Developing a Water Quality Assessment Framework for Southwest Florida Tidal Creeks

The three contiguous National Estuary Programs of southwest Florida, along with partners from six coastal counties, have developed a tidal creek water quality assessment framework to help prioritize natural resource investigations across a large population of tidal creeks between Tampa Bay and Estero Bay, Florida. The assessment framework is based on outcomes of a multidisciplinary study and includes a nutrient based report card that characterizes nutrient conditions relative to regional numeric nutrient criteria developed for contributing freshwater streams, identification of site-specific water quality indicators of tidal creek condition, and an online open science dashboard to display the assessment framework and provide access to all information relevant to its implementation. Application of the assessment framework has provided an actionable list of southwest Florida tidal creeks prioritized for further research and potential management action along with a host of site-specific indicator results that provide insights into drivers of tidal creek condition. The open science dashboard provides a platform for dissemination of this information in a readily accessible and reproducible format and a means to incorporate new data and indicators as they become available. Local resource managers are in need of tools to help prioritize natural resource investigations and management actions that achieve the greatest resource benefit with limited available resources. This assessment framework informs these efforts and builds capacity for future research to identify and refine management tools for these creeks where management resources, data, and sentinel biological response endpoints are limited.


Introduction
Tidal creeks are abundant in estuarine ecosystems, particularly in the low-energy coastal areas of the southeastern USA, yet their ecological relevance remains undervalued compared to larger, better-known estuarine systems (Mallin and Lewitus 2004). Tidal creeks provide critical habitat for many important estuarine dependent nekton species as a complementary ecotone to the larger estuary Holland et al. 2004;Greenwood et al. 2008;Krebs et al. 2010;Sherwood 2010) and their collective importance in material transfer and other ecological processes may equal or exceed that of larger estuaries (Dame et al. 2000;). The introduction of excess nutrients to coastal systems is a leading cause of observed water quality problems (Nixon 1995;National Research Council 2000;Crain 2007;Gilbert et al. 2010), and tidal creeks are critical pathways for nutrients between the landscape and the estuary. Tidal creeks receive and deliver nutrients from both anthropogenic sources (Mallin et al. 1999Mallin and Lewitus 2004;Holland et al. 2004;Sanger et al. 2015) and via decomposition of natural sources (Cloern 2001;Castro et al. 2003;Wolanski 2007;Deinenger and Frigstad 2019). The differences between these sources are important for management, yet tidal creeks remain understudied. In Florida, tidal creeks tend to be locally underrepresented in ambient water quality and environmental monitoring programs. This is due in part to difficulties in accessing what are often small tributaries surrounded by private lands or inaccessible floodplain areas. In addition, resource managers have tended to concentrate efforts on characterizing water quality delivered to the tidal portion of the system to estimate nutrient loads to the estuary rather than characterizing nutrient inputs and processes within the tidal portion of these creeks. Consequently, management and regulatory efforts to protect Florida's tidal creeks have lagged behind other waterbody types governed by the Clean Water Act. Southwest Florida includes several hundred tidal creeks that are not currently protected by defined numeric nutrient criteria but drain into prominent "estuaries of national significance" including Tampa Bay, Sarasota Bay, and Charlotte Harbor, many of which have overlapping designations as Outstanding Florida Waters or are within expansive state aquatic preserves. While numeric nutrient criteria have been developed for Florida springs, lakes, freshwater streams, and estuaries, tidal tributaries (including rivers and creeks) remain governed by a narrative criterion to protect waters from "an imbalance of flora and fauna" as defined in Florida Administrative Code (Rule 62-302.530(48)(b) F.A.C.) unless site-specific criteria have been established. Evaluation of the narrative criterion is based on chlorophyll a and dissolved oxygen standards developed principally for open bay estuaries. When violations of these numeric interpretations of the narrative criterion occur, nutrients are inferred as the causative pollutant. Local resource managers would benefit from water quality standards and numeric nutrient criteria developed using data collected within tidal creeks.
The need to develop a management strategy for southwest Florida tidal creeks was first promoted by the Tampa Bay Regional Planning Council in 1986 (Clark 1987;Clark 1991). Since that time, several important studies have been conducted on the value, function, and drivers of ecological integrity of southwest Florida tidal tributaries, including both larger tidal river systems and smaller tidal creeks. Sherwood (2008) and collaborators conducted a mutlidisciplinary tidal tributary habitat initiative in Tampa Bay, Florida, and concluded among other things that benthic microalgae were the primary source of food for fishes in tidal tributaries and that flashy stormwater runoff to tidal tributaries from more developed watersheds threatens benthic microalgae production by reducing water quality and light availability to benthic habitats. Burghart et al. (2013) used multivariate analyses of zooplankton community structure to develop breakpoints in the ratio of pelagic and demersal species composition as indicators of eutrophication based on sampling 4 spring fed and 4 surface water fed tidal rivers in southwest Florida. Peebles (2015) used zooplankton and hyperbenthos community metrics to develop indictors of trophic state for 6 tidal creeks in Tampa Bay, Florida, and found standard community diversity metrics (e.g. Shannon Weiner, Pielou's Evenness) may serve as valuable indicators of tidal creek trophic state. These studies, as well as several others Yeager et. al. 2007;Greenwood et. al. 2008;Sherwood 2010;Malkin et. al. 2010; Barbour and Adams 2012;Ley and Rolls 2018), were synoptic surveys designed to assess tidal creek conditions beneficial for estuarine dependent and resident fish, zooplankton, or hyperbenthos communities as biological response endpoints. These studies have contributed substantially to the identification of potential biological response endpoints for these systems as well as identifying drivers related to their success. However, these studies were not directly focused on identifying nutrient concentrations that could be used to list and prioritize creeks for management actions and gaps remain in defining the linkage between water column nutrient concentrations and biological response endpoints in a manner that can be used to develop discrete numeric water quality thresholds indicative of degradation of these systems. The three southwest Florida National Estuary Programs (NEPs) and their partners were interested in expanding research on southwest Florida tidal creeks to fulfill their own conservation and management objectives with respect to their proper stewardship and with the hope that this information may be used to inform future numerically based nutrient criteria for these systems.
This paper describes the development of a tidal creek water quality assessment framework to evaluate tidal creek condition across a large study area using a set of metrics derived from commonly collected water quality data. The goal of the effort was to coordinate and accelerate stewardship of more than 300 tidal creeks between Tampa Bay and Estero Bay in southwest Florida. Our assessment framework does not have any regulatory implications. Rather, the framework is designed to provide an openly available and continually updated assessment of tidal creek water quality condition that resource managers can use to screen and prioritize tidal creeks for nutrient management initiatives and/ or future research. We collected water quality and biological response endpoint data for a subset of creeks across a large study area, evaluated those data to identify water quality indicators and potential thresholds related to creek condition, and then applied the thresholds and indicators as an assessment framework to the identified population of creeks in southwest Florida using a large state and federal water quality repository. The objectives of using the assessment framework were to (1) report on tidal creek water quality conditions across a large population of tidal creeks between Tampa Bay and Estero Bay, Florida, (2) complement existing narrative-based regulatory standards for nutrients by providing additional context on site-specific creek conditions that may affect observed outcomes from current regulatory evaluations, and (3) prioritize investigations and potential management actions across the population of tidal creeks. The study involved identifying the population of creeks, conducting analyses of field data collections to develop a report card and indicator metrics, applying and testing metrics using a large long-term regulatory water quality database, and implementing the assessment tools with open source software in an online, interactive dashboard to directly support proactive local environmental management actions for southwest Florida tidal creeks.

Data Collection
The study area encompassed the watershed boundaries of the three southwest Florida NEPs between Tampa Bay and Estero Bay, Florida (Fig. 1). The National Hydrography Datasets (i.e., NHD-Plus;McKay et al. 2012), at 1:24 K and 1:100 K scales, were used in Geographic Information System (GIS) to identify the potential population of tidal creeks within the study area. Meetings were held with project partners to review the hydrography dataset, gather local knowledge on the creek population, and identify inclusion criteria which we defined as any natural or physically altered system (typically Strahler (1952) 3rd order or less) receiving freshwater inputs and exhibiting a range in salinities during a typical hydrologic year. The 1:100 K scale was initially chosen for identifying the population, which was then supplemented by additional creeks at the 1:24 K scale based on local knowledge and expertise of project partners. For instance, a small 2nd-order river (Estero River) was included in the population since its watershed area and creek length were smaller than other tributaries identified as "creeks" in the study area. We use the term "tidal creeks" to differentiate these systems from larger tributaries (i.e. higher order rivers) within the study area. The final identified population of 306 tidal creeks for this study is shown as black lines in Fig. 1. The tidal boundary (shading in Fig. 1) was defined using the Florida Department of Environmental Protection's (FDEP) Marine Water Body Identification Number (WBID) and FDEP saltwater interface boundaries along with a 1.66-km (1 mile) buffer which defined the landward most perimeter of the tidal portion of the creek.
Southwest Florida tidal creeks range from relatively natural systems to heavily modified ditches and urban stormwater conveyances, and a classification scheme was developed to balance the sampling effort across a range of potential conditions. First, a matrix of watershed attributes was developed that included (1) existing information on hydrologic soil groups and soil phosphorus concentrations (U.S. Department of Agriculture (USDA), Natural Resources Conservation Service (NRCS), National Soil Survey Handbook, n.d.)), (2) a landscape development index (LDI: Brown and Vivas 2005) for the creekshed and for a 50-m buffer from the creek centerline, (3) the percent of the creekshed that was classified as impervious surface based on the 2011 National Land Cover Database classification scheme (Wickham 2020), (4) elevation change within the tidal portion of the creek based on Digital Elevation Model data, (5) estimated nutrient loadings from previous NEP watershed loading models, and (6) several other metrics characterizing landscape attributes to further refine a classification scheme. Hierarchical agglomerative clustering was used to develop a cluster dendogram based on a resemblance matrix of watershed attributes and simple matching using normalized data values (Clarke and Gorley 2006). Three principal categories were separated at a cutoff value of 60% similarity (Class B, C, and D). To limit the size of the smallest class, the creeks with the lowest group average similarity relative to the three other classes were combined into a single additional class (A) and 16 Fig. 1 Geographic boundaries of the study area for defining the population of tidal creeks with the 16 sampled creeks bolded and labeled creeks were selected (4 from each class) for sampling. Class A creeks tended to be the longest creeks with the highest nutrient loadings, class B creeks tended to have the highest proportion of type A soils and the highest LDI scores, class C creeks tended to have the lowest LDI scores and the highest proportion of USDA NRCS B soils, and class D creeks tended to have the lowest soil P concentrations, but otherwise appeared to be similar to either B or C classes depending upon which watershed characteristic was compared. These classes were merely used to distribute our sampling across the range of watershed attribute and were not directly used to develop the assessment framework.
Two sampling surveys were conducted to collect data on water column nutrients, chlorophyll a, physical chemistry, and biological response endpoints in southwest Florida tidal creeks. Water quality analytes included constituents commonly collected in routine water quality monitoring programs while biological response endpoints included those identified in previous regional studies mentioned above. First, 16 creeks (labeled as bold lines in Fig. 1) were sampled every other month between November 2013 and September 2014. This survey (Janicki Environmental and Mote 2016) was designed to obtain unbiased estimates of the annual geometric mean nutrient concentrations relevant to typical Florida regulatory assessments, while also evaluating instream nutrient dynamics and biological response endpoints including benthic chlorophyll a content and nekton (fish) communities utilizing the creeks. The chlorophyll content of the top 1 cm of the creek bed (i.e. benthic chlorophyll a) was sampled using a 60-cm 3 syringe as a core based on outcomes of Sherwood (2008). The sample was collected either as a sub-sample of a petite ponar grab of the sediment or, if water levels were shallow enough, directly from the creek bed. Five replicate samples were collected at evenly spaced intervals along the left and right bank of the creek at the same randomly generated sampling point where water quality and fish samples were collected. Samples for each bank were then composited into a single sample for analysis. Benthic chlorophyll a content was determined using the modified fluorometric method of Whitney and Darley (1979) and Pinkney et al., (1994). Fish were collected according to Sherwood (2008) and others using a 9.1-m (3 mm mesh) raft seine (1.2 m high with a 1.2-m × 0.9-m × 0.9-m bag placed in the center) along the shoreline in depths 1.0 m or less to evaluate fish community structure. All species were identified to the lowest taxonomic level, measured in standard length, and enumerated.
For the second study (Janicki Environmental and Mote 2019), 6 of the previously selected 16 creeks (creeks with asterisk in Fig. 1) were chosen for sampling by our County partners (a single creek within each county) in consultation with the project team and sampled monthly during Spring and Fall months for 2 years beginning in April 2017. Benthic chlorophyll a and fish sampling was omitted from the second survey to focus on water quality dynamics over multiple years in a consensus-based decision with project partners given the fixed allocation of effort and results of the previous survey. The spring and fall seasons were chosen to bracket the traditional summer wet season when ca. 60% of the annual rainfall occurs in SW Florida and because nutrient transformations and phytoplankton uptake were observed to be most dynamic within the system during this time based on outcomes from the first study.
Both surveys collected data to evaluate water quality dynamics within the creek from the upstream extent of the creek's tidal segment to near the mouth of the creek by sampling four locations along the length of the creek. For each event, a fixed station located just upstream of the anticipated tidal head, and three sites located within three 200 m strata located along the length of the tidal segment were sampled (Fig. 2). The strata were not necessarily contiguous based on the length of the tidal portion of the creek. For the first survey, probabilistic site selection was used within the three tidal strata while for the second survey; a fixed site was established within each of the three strata at the center point. Sites were accessed either by small boat or by wading, dependent upon access to the location and depth of the creek. All creeks were sampled within a two-day window across the study area for each event and all samples were conducted during daylight hours, typically prior to noon and irrespective of tidal stage as the area exhibits a mixed semidiurnal tidal pattern (McPherson and Hammett 1991). All water quality samples were conducted in accordance with FDEP standard operating procedures (http:// www. dep. state. fl. us/ water/ sas/ sop/ sops. htm). Although the analytical methods defining the water quality constituents varied slightly between the two surveys (laboratory analysis for the second study included dissolved fractions and oceanographic methods with lower detection limits), analytical chemistry for both surveys was performed by the same NELAC certified laboratory, compared for reliability, and approved as part of a USEPA Quality Assurance Project Plan. Detailed information on analytical methods can be found in the supplementary materials for this paper (Tables S1 and S2).

Data Analysis
The field data were analyzed to assess potential stressorresponse relationships between nutrient concentrations and water quality and biological response endpoints that could be used to identify potential unifying numeric nutrient thresholds to protect tidal creek integrity. Data analysis methods included general linear models with tests for influence of random intercepts and residual correlation among samples collected within creek by date using Akaike Information Criteria for small sample size (AICC: Littell et al. 2006), analysis of variance with Tukey post hoc multiple comparison adjustment for fish catch statistics, robust regression (Rousseeuw and Leroy 2003), and Monte Carlo simulation (Davidson and Hinkley 2003) implemented using SAS software (SAS Institute, Inc. 2019). In addition, conditional probability analysis was conducted using classification and regression trees (Rokach and Maimon 2008) for exploratory data analysis and indicator development and implemented in R (R Core Team 2020). Specifically, conditional inference trees and ensemble models were implemented using the PARTY (Hothorn et al. 2006) and RandomForest (Breiman 2001) packages, respectively. These tools provide a flexible approach to evaluate nonlinear or threshold responses and can provide insight on conditional relationships among a host of potential explanatory drivers. Potential predictors included a host of landscape attributes including impervious surface, landscape development intensity index, natural buffer area, land use type, etc. (Table S3), as well as collected instream nutrient and physical chemistry water quality data. The response measures analyzed for this paper included water column and benthic chlorophyll a, dissolved oxygen, and fish/benthos community descriptors indicative of potential adverse impacts from nutrient pollution. The Standardized Precipitation Index was calculated using the SPEI package in R (Vincente-Serrano 2010) to evaluate the sampling effort relative to rainfall anomalies and assess potential relationships between rainfall anomalies and water quality.

Report Card and Indicator Development
Data from our field surveys were also used to inform development of a water quality "Report Card" and "Indicators" of tidal creek condition and these tools were then applied to the entire population of tidal creeks using data from the Florida statewide regulatory Watershed Information Network (WIN) database. The report card was derived as a screening tool to characterize creeks based on concentrations of the limiting nutrient (i.e. nitrogen) in the system. The report card utilized a "reference condition" approach with the reference being the established numeric nutrient criteria for the contributing freshwater stream. The current regulatory standards for total nitrogen in contributing freshwater streams in our study area are based on annual geometric means (AGMs) of 1.65 mg/l for the West Central region which includes Hillsborough, Manatee, Sarasota, and Charlotte Counties, and 1.54 mg/l for the Peninsula region which includes Pinellas and Lee counties. A 5-level classification scheme was developed based on these reference values and each year of data was categorized into one of 5 levels ("Prioritize," "Investigate," "Caution," "Monitor," "No Data"). A 10 year period analogous to Florida's water quality assessment "planning period" was evaluated and each year was assigned category. Creeks exceeding the freshwater numeric thresholds were assigned the "Prioritize" level. A proactive margin of safety was used to define the "Investigate" level category using Monte Carlo simulation and the standard deviation based on the pooled data from our first study. The simulation process used the expected value and standard deviation of log transformed observed data to derive a large simulated data pool (10,000 observations). Assuming six sampling events per year, a distribution of AGMs was generated based on 1000 replicates and the number of AGMs exceeding the regulatory threshold was calculated. The expected value was then incrementally shifted until a distribution with a 5% probability of exceeding the regulatory threshold for the contributing freshwater portion of the creek was identified. The expected value of that distribution was then chosen as a threshold AGM value for the Investigate level category.
The Prioritize and Investigate levels are regionally based and proactively identify creeks that may be approaching impaired status for further study and potential management action. However, an additional goal of this study was to provide guidance for proper stewardship of individual creeks that may require more site-specific nutrient management strategies. As such, a "Caution" category was developed to identify creeks where nutrient concentrations were considered elevated given site-specific conditions. Preliminary analyses demonstrated that several creek watershed attributes were correlated with the AGM concentrations for total nitrogen and total phosphorus, including creek length, watershed size, buffer acreage, area of natural wetland vegetation (predominantly mangroves) within tidal reaches, and the AGM concentration of colored dissolved organic matter. A simple linear regression relationship between nitrogen concentration and the estuarine total creek length (Table S4) was chosen as a scaling factor for defining the lower limit of the Caution category, due to its broad applicability for all creeks whether or not water quality data or other watershed attributes were readily available. The Caution level was constructed to protect smaller creeks from elevated nutrient concentrations as a stewardship goal. Finally, a "Monitor" category was designated for any creeks with existing data that fell below the Caution level criterion value, and a "No Data" category was assigned to those creeks where data were unavailable for evaluation.
A scoring algorithm was derived to define the final report card outcome for each creek using the entire ten year record of available data based on the following criteria. A single exceedance of the Prioritize and Investigate categories in any year of the ten year record would result in a classification of that creek into the respective category unless at least three other years of data were below the threshold level for that category. Creeks were assigned the next lower category if only one AGM for total nitrogen (TN) was above a given level while multiple other years (i.e., more than two) were below the given levels defining the cutoff points for each category. For example, a creek with at least 4 years of data and only a single exceedance of the Prioritize threshold would be assigned the Investigate category. The WIN database includes assigned waterbody identifiers (WBIDs) intended to represent discrete hydrologic units. Since it was possible for multiple WBIDs to exist within the tidal portion of the creek, the report card evaluation was carried out by WBID. In rare cases where multiple WBIDs exist within the tidal creek, the creek was assigned the highest observed category for the report card display. A schematic representation of the data synthesis process is provided in Fig. 3 where the raw data timeseries is provided (Fig. 3a); each year's data is assessed relative to the categories (Fig. 3b). The 10 years of data are summarized by the frequency counts in each category (Fig. 3c) and the final report card outcome score is derived (Fig. 3d).
Water quality Indicators were developed based on outcomes of our study as described in detail in the results section. The results for each indicator were calculated on an annual basis and then synthesized for the 10-year period by calculating the percentage of annual outcomes above the identified threshold indicator values out of the total number of years with available data. An integrative summary for all indicators was presented using a standardized polar coordinate system and radar charts to provide a single multimetric summary plot of the results across indicators.

Tidal Creeks Dashboard
An online, interactive dashboard was created for the assessment framework to help automate and distribute summary information on tidal creek condition to regional environmental managers. The dashboard is also freely available for the general public to better understand the importance of tidal creeks as environmental resources that support healthy estuarine ecosystems. All dashboard components were created using the open source R Statistical Programming Language (R Core Team 2020). Specifically, the R Shiny package (Chang et al. 2020) was used to create the interactive web interface and the "tbeptools" R package (Beck et al. 2021) created by the authors was used to read, analyze, and graph data for the tidal creek indicators.

Study Results
A total of 672 water quality samples were collected between the two studies (Table 1). TN concentrations ranged from 0.41 to 2.32 mg/l. Total phosphorus (TP) ranged from 0.02 to 3.19 mg/l. Corrected chlorophyll a (Chla) concentrations ranged from 0.37 to 229.0 µg/l and were highly variable. Episodic spikes over 100 µg/l were observed in the tidal portion of the creek, though average concentrations were typically less than 20 µg/l. Salinity ranged from 0 to 32 psu with the fresh water stations mostly near zero though higher salinities were recorded at the most upstream station of a few creeks for some dry season sampling events during the first study. Turbidity ranged from 0.6 to 340 ntu, and dissolved oxygen (DO) ranged from 0.17 to 12.50 mg/l. The tidal portion of the creeks were principally nitrogen-limited, with median inorganic N/P ratios well below 16 (Redfield limitation) for all creeks except Estero River at the southern end of the study area (Fig. S1).
Our sampling effort generally coincided with standardized precipitation index anomalies within 1 standard deviation of long term average monthly conditions (Fig. S2a); however, differences in rainfall deviations across creeks for a particular month were evident (Fig. S2b). While ca. 60% of the rain typically falls within 4 months of the year (June-September) and Spring months (i.e. April and May) tend to be the driest months of the year, convective activity and the potential for early season tropical storms modulate these general patterns. Seasonality in rainfall is a principal forcing function for water quality in these creeks, and seasonal dynamics in water quality generally conformed to the expected patterns. Dissolved oxygen concentrations and percent saturation values were lower in the summer, while TN and TP concentrations tended to be highest in summer months. Chlorophyll a concentrations tended to be highest and most variable in the spring dry season, but Chla results were particularly site-specific and were affected by antecedent rainfall conditions throughout the year. A representative example of the dynamic nature of instream tidal creek water quality dynamics from our study is provided in Fig. 4 for samples collected in the spring of 2017 in Estero River. A high antecedent rainfall condition (31.5 cm of rain in the 14 days prior to a June 2017 sampling event resulted in water quality with near zero psu salinities throughout the sampling area (4a), high color concentrations throughout the tidal portion of the creek (4b), relatively high TN throughout the tidal portion of the creek (4c), and low Chla, as residence times were presumably insufficient for phytoplankton production (4d). Alternatively, the April and May 2017 sampling events, with 0.5 cm and 4.3 cm of rain in the 14 days prior to sampling, respectively, exhibited more typical dry season conditions. During the dry season, a strong salinity gradient, lower color, lower TN, and higher overall Chla concentrations were observed in the tidal portion of the stream. Interestingly, the highest Chla concentrations were associated with some antecedent rainfall but not enough to have a marked effect on reducing salinity. These events all occurred at or near high slack tide with consistent depths of between 1.2 and 1.6 m.
Bivariate stressor-response relationships between nutrients and water quality response endpoints were observed when grouping all data among creeks, though standard linear regression model fit was poor (Fig. 5). Significant (α = 0.05) positive relationships were observed between Chla and both TN and TP (Fig. 5a, b). Dissolved oxygen was negatively related to TN concentrations (Fig. 5c). These results conform to expectations under the current management paradigm; however, the majority of variation was left unexplained by these models making inferences based on these simple relationships unreliable for identifying protective criteria for the population of creeks. In addition, there was a lack of relationship between Chla concentrations and DO (Fig. 5d). Using random intercepts to account for between creek effects generally improved linear regression models by absorbing inter-creek variation from the error. Incorporating a random effect for residual date correlation yielded little improvement in AICC, indicating little, if any, effect of event-based correlation introduced by the sampling scheme (Table 2). Evaluating creek specific bivariate relationships suggested that in many cases the assumption of a constant slope would be erroneous. For example, in some creeks, Chla was positively associated with DO, while in other creeks, Chla was negatively associated or not associated at all with DO. This contributed to our inability to identify a unifying relationship between these two constituents.
This study confirmed that southwest Florida tidal creeks serve as critical nursery areas for important estuarine sport fish including the common snook (aka snook; Centropomus undecimalis), a species of considerable economic value in southwest Florida. Impressively, snook was the ninth most abundant species collected. Six other recreationally important species were found in the top 25 most numerically dominant species collected including the red drum (Sciaenops ocellatus), blue crab (Callinectes sapidus), pink shrimp (Farfantepenaeus duorarum), striped mullet (Mugil cephalus), and sheepshead (Archosargus probatocephalus) as well as a number of species that serve as prey items for estuarine dependent sport fish (Janicki Environmental and Mote 2016). Catch per unit effort of snook aggregated across the sampled creeks was significantly higher in the tidal creeks we sampled compared to larger tidal river systems in Tampa Bay routinely sampled by partners over the same period of time (Fig. 6), confirming that these creeks are critical habitats that serve as nurseries supporting snook early life stages. The length frequency distribution of snook captured in this study included both young of the year and age 1 snook which are underrepresented in data currently being utilized for snook stock assessments. However, snook catch, including density, length frequency or biomass, was not significantly correlated with nutrient concentrations in our study and assessments of the fish community composition suggested that while there were substantial differences between creeks, fish diversity metrics (after removing exotic species) were either not significantly or positively correlated with observed nutrient concentrations (Table 3). The results of our field surveys suggested that underlying stressor-response relationships were likely affected by characteristics of the watershed, nutrient inputs from both natural and anthropogenic sources, and physical characteristics of the creek itself. Even our most natural creeks had some evidence of anthropogenic impacts and the interaction among many types of anthropogenic influences limited our ability to identify unifying water quality thresholds that would provide both sensitivity in identifying creeks with nutrient related impairments to the response endpoints we tested while also providing specificity in identifying creeks where similar nutrient concentrations were not resulting in imbalances of flora and fauna that would be deleterious to creek condition. The conditional nature of these relationships was confirmed by conditional inference models which suggested that a large combination of habitat-related factors were required to explain substantial proportions of the variation in water quality response endpoints. The most influential drivers for nutrient concentrations included creek length and watershed buffer area characteristics (including agricultural areas within the buffer) while Chla responses were more influenced by other instream water quality characteristics (Table 4). The degree of impervious surface in the watershed was a significant but not dominant predictor of instream nutrient of chlorophyll concentrations in our study; perhaps due to the heavily urbanized nature of much of our study area. Detailed results evaluating the many aspects of this multidisciplinary study including investigating relationships between fish catch and Chla and DO conditions, water quality comparisons between creeks, and water quality dynamics along the salinity gradient can be found in Janicki Environmental and Mote (2016 and 2019).
The remainder of this paper focuses on the development and application of our assessment framework. Instead of proposing a single criterion value that would likely be mis-specified for a large proportion of creeks within the systems, we developed an assessment framework that was completely outside of any regulatory a b c d  framework but that could be used by the NEPs and their partners as a screening and prioritization tool to target creeks at risk of nutrient related adverse effects for future investigation and research. The assessment framework includes a nutrient report card and evaluation of a number of water quality indicators developed to provide more site-specific insights on contributing factors to observed adverse effects related to alterations in the system that reduce assimilative capacity and resiliency to nutrient inputs.

Report Card
Report card outcomes for the period of record 2009 to 2018 are mapped in Fig. 7. Ninety-eight creeks (108 WBIDs) within the population had data collected within the 10 year assessment period. Of those WBIDs with data, 13 were assigned as Prioritize, 16 as Investigate, and 10 as Caution (Fig. 8). The remaining 69 WBIDs with data were assigned to the Monitor category indicating that nutrient conditions in these creeks were likely not adverse to the biological endpoints we measured. The majority of WBIDs had no data reported during the assessment period, highlighting the need for more monitoring to support stewardship of these resources. The flexibility of the scoring algorithm to classify a creek can be seen by the shading of cells in Fig. 8 which depicts the frequency of years classified into each category. While most creeks had multiple years of data, 27 creeks had data only for a single year; however, the utility of evaluating even a single year was highlighted by our findings that several of these creeks exhibited nitrogen concentrations above our Caution criteria. This information is critical in data limited systems. Several creeks in the Investigate or prioritize categories had only episodic spikes in nitrogen concentrations and several years of data where TN concentrations were below the margin of safety providing important information on nutrient inputs to the system while other systems were consistently above the threshold for the Investigate category.
The cumulative distributions of the TN AGMs for each report card category were visually distinct across the range of observations (Fig. S3). The distribution of AGMs for creeks in the Prioritize category was right-skewed with the majority of AGMs exceeding their respective contributing freshwater nitrogen criterion values. However, several of these tidal creeks are not listed as impaired under the traditional, narrative regulatory assessments that rely on indirect assessments of nutrients to identify an "imbalance of flora and fauna." The Caution and Investigate categories were similar to each other for the lower quartiles of their respective distribution, but differed in the upper quartiles where the Investigate distribution was slightly right-skewed. The Investigate distribution had, by definition, at least one AGM concentration for TN above the prescribed margin of safety, explaining the difference at the upper tails. The Caution category was site-specific to protect smaller creeks from elevated TN concentrations, and an observed overlap with the Investigate category for larger creek systems was not unexpected, as larger creeks grading category thresholds were higher.

Tidal Creek Indicators
The report card is one aspect of the larger assessment framework to evaluate tidal creek condition. The development of indicators to evaluate site-specific water quality attributes was based on results from our study and the generalization of those results to the larger population of creeks. For example, despite weak linear relationships, conditional probability analysis suggested a highly significant difference (p < 0.001) in Chla based on a TN concentration threshold of 1.1 mg/l. When TN concentrations were at or below 1.1 mg/l, the median Chla was below the current regulatory numeric interpretation of the narrative criterion (i.e. Chla AGM of 11 µg/l) applied to tidal creeks. This threshold TN  (Paulic et al. 1996). Therefore, an AGM TN of 1.1 mg/l was identified as one indicator of tidal creek condition. Carstensen et al. (2011) used the ratio of Chla per unit TN (Chla/TN) as an indicator of a shifting ecological baseline in the relationship between the two water quality constituents and we considered this ratio as a useful indicator of phytoplankton response to nutrient inputs in tidal creeks. In our study, tidal creeks with median Chla/TN ratios higher than 15 were observed to have a greater degree of physical alteration to either creek morphology or source water hydrology (West Spring, Forked, and Buck Creeks; Fig. S4) though the degree of physical alterations was difficult to quantify. Physical alterations to creek morphology are common in southwest Florida and can alter residence times and disrupt the natural hydrologic regime, affecting creek biogeochemistry and ecology. A high Chla/TN ratio may indicate more autotrophic conditions, while extremely low ratios may indicate dysfunction in expected nutrient assimilation pathways that manifest as Chla. Based on our observation that tidal creeks with median Chla/TN ratios higher than 15 tended to Fig. 8 Summary of prioritized creeks displaying categorized annual outcomes (first four columns) and the final report card score (final column). Cells in the first four columns are shaded (opacity) by the number of years that were assigned to each category within the 10-year evaluation period have a greater degree of physical alteration to the system, an annual median Chla/TN ratio of 15 was identified as another indicator of creek condition.
Inorganic nitrogen is readily assimilated by biological processes in nitrogen-limited systems and is also a common indicator of anthropogenic inputs (Gowen, 1992). Inorganic nitrogen was typically assimilated rapidly along the length of the creeks in our study, indicating that the dilution or uptake of inorganic nitrogen by biological processes could be a useful indicator of creek condition. However, there were cases where inorganic nitrogen inputs within the tidal portion of the creek were observed, indicating potential anthropogenic inputs. A plot of the source water and tidal inorganic nitrogen concentrations relative to a 1:1 line (Fig. S5) was used to identify creeks that may have anthropogenic inputs within the tidal portion of the creek. In our study, Estero River was a consistent example having an inorganic source of nitrogen of currently unknown origin within the tidal portion of the system. We included this ratio with a threshold value of 1 as an indicator of inorganic nitrogen inputs into the tidal portion of the creek though its assessment requires sampling at locations in both the tidal and freshwater portions of the creek.
Florida's Trophic State Index (TSI) (Paulic et al. 1996) was used as a water quality screening tool for Florida lakes and estuaries prior to the development of more site-specific criteria and was included as a tidal creek indicator to evaluate the utility of a composite index, combining both Chla and TN (and in some cases TP, depending on the N/P ratio). The Florida TSI is a modification of Carlson's (1977) Trophic State Index that excludes water clarity due to the blackwater characteristics of Florida waters but combines measurements of Chla, TN, and TP concentrations into a single value ranging from 0 to 100 to describe a waterbody's trophic state, with higher values suggesting more eutrophic, degraded conditions. For Florida lakes, values above 70 are rated a "poor" condition, values from 60 to 69 are "fair," and values below 60 are rated "good." Category breakpoints for Florida estuaries are ten points lower to reflect a lower productivity relative to lakes (Paulic et al. 1996). Application of the Florida TSI to the population of tidal creeks in our study resulted in a median value of 47 suggesting the majority of creek years were in the "good" estuary category. Scores ranged from 17 and 87 with only 5% of the scores above the threshold used to denote a poor estuary (TSI 60, Fig. S6). We selected a threshold value to a TSI threshold score of 55, the midpoint of a "fair" estuary, to indicate that the productivity of a creek may be shifting towards a more eutrophic condition.
An example synthesis of our multimetric indicator results, which include the 4 indicators described above as well as the regulatory Chla and DO standards currently used to evaluate tidal creeks under the narrative criterion, is provided for four example creeks in Fig. 9 using radar charts. The radar charts present the percentage of annual evaluations that exceed the threshold values for each indicator. For example, Fig. 9a displays results for a tidal creek with a high percentage of TN indicator exceedances, high percentage of Chla exceedences, a high number of dissolved oxygen exceedances (i.e. dissolved oxygen below the threshold criteria), and the majority of years above the selected TSI threshold value. This creek conforms to the traditional management paradigm where high nutrient concentrations are correlated with high Chla and low DO conditions. Figure 9b depicts an outcome with no TN indicator scores above the 1.1 mg/l threshold, but the majority of years exceeding the regulatory standards for both Chla and DO. The nitrate ratio threshold in this creek was consistently exceeded indicating that despite relatively low TN, nitrate inputs into the tidal portion of the creek were evident. This combined with the high exceedance of the Chla/TN ratio thresholds indicates that this creek may be particularly susceptible to inorganic nutrient inputs and/or be highly physically altered, requiring a more site-specific management approach. Figure 9c depicts results for a creek where the TN indicator score and TSI score were consistently exceeded but other metrics were not. The creek represented by Fig. 9d had high nitrate ratio score and TSI exceedences despite TN AGMs typically less than 1.1 mg/l. The Chla AGMs were routinely exceeded, and the Chla/TN ratio was exceeded in the same years indicating this creek may also be susceptible due to physical alterations or inorganic nutrient inputs. Together, these results highlight the complex interaction among individual water quality indicators and the need for additional information to supplement current narrative based tidal creek water quality standards. The radar charts provide an easy way to summarize multi-metric indicators and provide insights into sitespecific attributes that may govern water quality responses in our tidal creeks. This information can be used to inform site-specific investigations into tidal creek condition. The application of these tools using the open science dashboard provides an efficient platform for synthesizing and disseminating this information and allows the underlying science to advance over time.

The Management Framework as an Open Science Product
The tidal creeks dashboard (https:// shiny. tbep. org/ tidal creekdash/; Beck and Wessel 2021) was created to leverage open science concepts that enable transparency, reproducibility, and accessibility of methods. The dashboard provides a simple user interface to the tidal creeks assessment results and the underlying methods and source code are freely available in the supporting "tbeptools" R package (Beck et al. 2021;https:// tbep-tech. github. io/ tbept ools/ artic les/ tidal creeks. html). This approach facilitates understanding of the assessment framework among stakeholders and can also encourage others to apply similar techniques by providing access to underlying code behind the analysis methods. Such approaches have been shown to promote "better science in less time" (Lowndes 2017). There are several components of the dashboard that facilitate understanding of the assessment framework. The main components include an overview page that describes how to use the dashboard and a brief summary of the methods, a report card page that provides mapped and summary outcomes for the population of creeks, and an indicator page to show detailed results for each creek. A useful component of the dashboard for local resource managers is the ability to select broad regions or specific creeks of interest. For example, an interactive map for the report card page allows the user to draw a custom polygon over an area of interest to select creeks within a particular jurisdiction. Results for creeks within the polygon are then shown in tabular format for download and also summarized by percentage in each reporting category, either by the number of creeks or by total creek length. The indicator page allows a user to select a particular creek to view detailed results including a complete time series of the ten years of data for TN and supporting indicators described above (e.g., Chla, TSI, DO). This provides managers additional context to understand site-specific creek condition within the context of the population. The ability to interact with and extract results of interest for specific creeks provides unlimited customization to retrieve results of interest, helping facilitate an understanding of larger patterns while also informing more detailed assessments that can better guide discrete and/or regional management activities. The open science aspects also facilitate broader adoption and application of these methods beyond southwest Florida.

Discussion
The assessment framework was developed to provide actionable outcomes for southwest Florida tidal creeks based directly on commonly collected water quality data reported in state and federal repositories. Outcomes can be used to identify and prioritize natural resource investigations and management actions for these tidal creeks, while also promoting and accelerating additional research that may be needed to better understand these systems. The adaptive management paradigm is promoted and facilitated by the open science application of the dashboard. The nitrogen report card characterizes concentrations relative to criteria established for contributing freshwater streams while the indicators provide context to the report card outcomes by evaluating other water quality metrics as well as nutrient dynamics between the fresh and tidal portion of the creek where data are available. It is important to note that the assessment framework is not a regulatory tool but is designed to complement existing regulatory assessments related to nutrient impacts in southwest FL tidal creeks and, if pursued in a proactive manner, may even preclude the need for more stringent regulatory interventions in some cases.
The report card outcome categories of Prioritize, Investigate, Caution, and Monitor provide local natural resource managers needed information from which to allocate limited resources to prioritize and pursue potential nutrient management activities. Several creeks identified as Prioritize or Investigate by our report card have current or legacy anthropogenic nutrient sources that have/had the potential to deliver large nutrient loads to the downstream estuary. In several cases, these tidal creeks would pass an "impairment" assessment utilizing the current narrative regulatory standards. Although the Caution category is based on a somewhat weak relationship between creek length and nitrogen concentrations, it was chosen as a precautionary measure to protect smaller creeks from elevated nutrients concentrations. As such, follow-up assessments using the indicator outcomes are designed to support the ultimate decision making process developed from initial report card grade. The report card utilizes data from both long-term monitoring programs and synoptic studies and has proven useful as a screening tool to characterize tidal creek condition even when only very limited data are available. We observed that even using a single year's worth of data in the assessment could be very informative. In one case, a creek was categorized into the Investigate category based on only 4 samples and subsequent investigations identified additional data associated with a point source discharge in this creek that was not represented in the statewide data repository. Resource managers can use the dashboard application to easily facilitate these types of desktop evaluations.
Our study results suggested nitrogen limitation for the majority of southwest Florida tidal creeks, and therefore, our report card was based on TN (the form of nitrogen used to establish numeric nutrient criteria for regulatory assessments in Florida with the exception of springs). However, the report card is easily adaptable to cases where phosphorus may be identified as the limiting nutrient or when colimitation of both nutrients may be important for a particular system. A large portion of our study area is characterized by a phosphorus rich sediment layer in west central Florida known as the "Bone Valley," resulting in high surface water phosphorus concentrations relative to other surface waters of the State; though, we did observe some creeks south of the Caloosahatchee River (Estero River and Spring Creek) that exhibited N/P ratios indicative of potential for phosphorus limitation at times. Consistent with Cloern et al. (2020), our results also suggested that there was substantial variability in N/P ratios dependent upon location along the tidal gradient, seasonality, and which elemental form was evaluated. The range of nutrient and Chla concentrations we observed was extremely similar to those reported for North Carolina tidal creeks ) where nutrient addition bioassay experiments suggested nitrogen may be limiting even with N/P ratios are variable. However, we do not completely rule out the potential for phosphorus limitation under certain conditions in some locations within our study area.
Our field surveys supported indicator development useful in identifying site-specific attributes related to the water quality dynamics and the expression of nutrient inputs in these systems. In lakes and open bay estuarine systems, stressor-response relationships between nutrients and chlorophyll a can be reliable as residence times are generally sufficient for phytoplankton reproduction if a sufficient mass of limiting nutrients is available. In fluvial systems including tidal creeks, nutrient-chlorophyll relationships are more complex. In fact, Florida's freshwater stream numeric nutrient criteria (62-302.531 F.A.C.) were based on a reference system approach due to an inability to define a stressorresponse relationship that would serve as unifying protective criteria for freshwater streams. The regulatory requirements for freshwater streams now include a biological confirmation step using stream condition index scores and periphyton indices to validate or supersede outcomes based on application of regulatory water quality criteria. Evidence from our tidal creeks study suggested that nutrient-chlorophyll relationships are mediated by physical alterations to the system and the form of nitrogen available, as well as riparian buffer and landscape attributes. Together, these site-specific confounding factors limited the ability to derive a single nutrient threshold that would be protective of all creeks without being unnecessarily restrictive for many creeks with a robust fish community. Instead, a host of indicators was used to provide site-specific context to observed water quality and the report card outcomes.
The indicator outcomes we observed exemplify the broad range of responses even among creeks with similar nitrogen concentrations. We found the radar charts to be an intuitive way to facilitate evaluating these multi-indicator outcomes. Some creek responses appeared directly related to total nitrogen concentrations while others may be related to inputs of inorganic nitrogen (as measured by nitrate) in the tidal portion of the creek as well as physical alterations to the system that may make certain creeks more susceptible to nutrient pollution than others. Typically, nitrates are rapidly assimilated by biological processes and the finding of exceedances in our nitrate ratio indicator was used to indicate the presence of a source (e.g. septics) in the tidal portion of the system. This outcome often correlated with exceedances of the Chla/TN ratio, though fewer creeks had nitrate measurements in both the tidal and contributing freshwater portion of the stream. The Chla/TN ratio adds further context as an indication of how nutrients are assimilated by water column phytoplankton in these systems. The Chla/TN ratio is also likely influenced by water residence time and may be a good surrogate descriptor of hydrological or physical alteration of a creek. Resource managers can use these metrics to pursue evaluations of nutrient sources and identify potential limiting physical conditions based on these outcomes. The TSI indicator was included to evaluate a composite water quality index within the assessment framework, but its utility as an indicator relative to the individual constituents that make up the TSI is still being evaluated.
Physical alterations of the landscape are critical components of tidal creek ecological condition (Holland et al. 2004;Sanger et al. 2015), and the southwest Florida landscape is dominated by highly urbanized areas but no single landscape level indicator stood out as a differentiator for tidal creek condition in our analysis. Anthropogenic effects were evident even in our most natural systems which could have affected these results and limited our ability to define a "reference system". Historically, many southwest Florida tidal creeks have been physically altered through dredging and channelizing for stormwater conveyance, flood control, recreation, and or navigation but have the same regulatory "designated use" for evaluating compliance with water quality standards under the Clean Water Act. These physical alterations can either dramatically increase the water volume relative to watershed size or disconnect the natural hydrology. Both of these conditions can significantly affect residence times, and it is interesting to note that hydrologic alteration is now being considered as a source of pollution within the USEPA Clean Water Act (Novak et al. 2016). Ellis et al. (2017) suggested that creeks with relatively large discharge volumes could more efficiently dilute contaminant concentrations from terrestrial runoff compared to smaller, more urbanized creeks, suggesting the latter are more vulnerable to increased runoff from developed landscapes. Very few creeks in our study had both gaged flows and bathymetric data to estimate residence times, which is a current data gap critical to better understanding the role of freshwater discharge as a physical driver of water quality dynamics. The pursuit of additional indicators related to tidal creek condition, including the role of physical alterations on system function, the role of emergent wetland vegetation as a nutrient source and habitat, and the pursuit of additional biological response endpoints, is a recommended direction for future study.
Establishing biological response endpoints indicative of adverse effects to tidal creeks due to nutrient pollution was a challenging endeavor. The decision to use existing freshwater criteria values as the basis for the nutrient report card was based on the lack of evidence from our study that the response endpoints we evaluated were more sensitive to nutrient pollution than those identified for contributing freshwater criteria. Our study found that these creeks are critical habitats for several important fish species including the snook which supports similar findings by other studies in the region Yeager et. al. 2007;Greenwood et. al. 2008;Sherwood 2010;Barbour and Adams 2012;Ley and Rolls 2018). Fish dose-response curves form the basis of the marine dissolved oxygen standard in Florida (FDEP 2013) that is currently applied to tidal creeks providing precedent for the use of fish as an indicator of adverse effects. However, neither individual species catch statistics nor community diversity metrics appeared to be adversely affected by the nutrient conditions we observed, even when those concentrations approached regulatory criteria for their contributing freshwater streams (Janicki Environmental and Mote 2016). This indicated that the faunal endpoints we evaluated may not be most sensitive to the effects of nutrient pollution. This finding does not diminish the importance of monitoring these creeks as fish habitat as our initial collaboration meetings established that a consensus of local resource managers view a healthy fish community as a primary goal for these systems which may also increase their assimilative capacity. While our goals are to continue to focus on developing robust biological response endpoints as indicators of tidal creek condition, our framework is currently based on commonly collected water quality data to provide an actionable list of creeks from which future research can focus to improve our understanding of the many aspect that constitute tidal creek biological integrity.
The three southwest Florida NEPs currently have dedicated management objectives in their individual Comprehensive Conservation and Management Plans that seek to identify indicators of tidal creek integrity and increase stewardship and management efforts in these systems. Development of this framework and its open science tools helps achieve those objectives. Currently, the assessment framework and dashboard allow natural resource managers to prioritize investigations into their tidal creeks and researchers to easily evaluate the distribution of nutrient and indicator conditions for these creeks to aid in developing future studies. The open science tools included in this study are expected to accelerate this effort both locally and in other regions where a water quality assessment framework for tidal creeks is desired. Extrapolation of these methods to other areas would require an accessible database with tidal creek water quality data, spatially referenced flow lines for tidal creeks, and numeric nutrient standards or management level thresholds for their contributing freshwater reaches. The assessment framework and the resulting open science tools are intended to assist resource managers in prioritizing investigative and proactive management actions in southwest Florida tidal creeks, thereby improving and accelerating the proper stewardship of these critical coastal habitats facing increased pressure from coastal landscape development, anthropogenic nutrient inputs, and climate change.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.