Comparative Use of Hydrologic Indicators to Determine the Effects of Flow Regimes on Water Quality in Three Channels across Southern Florida, USA

: This study determines the relationships between water ﬂow and water quality in three types of channels in southern Florida, USA: Shark River Slough, Peace River, and Hillsboro Canal. Peace River most resembles a natural channel with ﬂoodplain connectivity, sinuosity, and uninhibited ﬂow. Shark River Slough has a natural, shallow channel with sheet ﬂow, while the Hillsboro Canal is the most modiﬁed channel due to dredging, straightening, and regulated ﬂow. Hydrologic indices for each channel were estimated to characterize ﬂow regimes and ﬂow variability, while concentration– discharge ( C–Q ) relationships were determined to quantify the impact of ﬂow regime on water quality. The greatest variability in ﬂow occurred at the Hillsboro Canal, followed by Peace River and Shark River Slough. Connectivity to ﬂoodplains and long durations of low and high ﬂow pulses at Peace River and Shark River Slough contributed to the dilution of water quality constituent concentrations at higher ﬂows. Conversely, the channelized characteristics of the Hillsboro Canal resulted in an enrichment of constituents, especially during high ﬂows. This study suggests that C–Q relationships can be used in canal discharge management to prevent water quality degradation of sensitive downstream wetland and aquatic ecosystems. heavily channelized Hillsboro Canal, whose ﬂow regime is regulated by water control structures, had the greatest ﬂow variability, and was dominated by enrichment behavior of its water quality constituents with higher ﬂows. This study found that C–Q relationships can be applied to canals and should be considered by water managers in discharge operations to reduce the transport of nutrient and pollutant loads downstream to ecologically delicate ecosystems.


Introduction
Human modifications of natural channels (rivers and streams) have resulted in channelization and deviation from natural flow regimes [1][2][3][4]. Channelization, or canalization, involves the straightening and shortening of existing natural channels to increase hydrologic transport capacity, and the dredging of land surfaces to create new channels for navigation, water conveyance, and drainage purposes [2,5]. Channelization includes the construction of levees and embankments along channel banks to restrict channel morphology and width, which, in some instances, results in the deepening of the channel [3]. Such changes in morphology disconnect channels from their floodplains, reduce sediment exchange capacities, prevent overbank flows, and change flow regimes [3,6,7]. Flow management via water control structures can further change the flow regimes of channels [5]. For example, pump stations, culverts, spillways, and weirs are used to regulate the flow of water within the canal networks of South Florida, USA [5,8]. Overall, flow regime changes affect the natural variability and cyclical patterns of low and high flow conditions of channels, often changing the magnitude and frequency of low and high flows [3,9,10].
A variety of hydrologic metrics have been developed to quantify deviations from natural flow regimes that are caused by human interventions. The Indicators of Hydrologic Alteration (IHA) are widely used hydrologic indices that characterize the magnitudes, timings, frequencies, and rates of changes in water flow, all of which are components of a flow regime [10]. These components are important because they regulate bedload transport, which influences channel morphology and facilitates exchanges of water quality constituents (e.g., nutrients and organic matter) between a channel and its floodplain [3]. Furthermore, flow regime changes in modified channels can cause water quality degradation by decreasing in-stream nutrient retention and increasing pollutant exports downstream to receiving ecosystems [11]. Low flow conditions can increase the residence time of water quality constituents [12], allowing biogeochemical processes such as biotic uptake, microbial reduction (e.g., denitrification), and sorption to regulate concentration variability [12,13]. Conversely, high flows can dilute the concentrations of point-sourced water quality constituents or mobilize concentrations of multi-sourced constituents, including particulates, which can then be transported downstream [12,[14][15][16].
Regression models have been developed to predict the behavior of the concentrations of constituents that affect water quality for specified flow regimes [17]. Specifically, several studies have used concentration-discharge (C-Q) relationship models to quantify the effect of flow on the behavior of water quality constituents within natural channel networks [12,15,16,[18][19][20]. A C-Q relationship is commonly expressed as a power-law function: where a and b are model coefficients that represent the intercept (same unit as concentration) and slope (unitless), respectively, on a logarithmic scale [12]. The slope (b) represents the pattern or behavior of the constituent [18]. A chemostatic or constant pattern, where b = 0 [18], occurs when the concentration of a constituent is not influenced by discharge, but by other factors such as biogeochemical cycling. For instance, chemostasis can result from the abundance of a constituent in the watershed [21], stemming from nutrient legacies from agriculture [22] or active geogenic weathering [18], which act as biogeochemical controls and regulate the concentration and loading of such constituent. Conversely, a chemodynamic pattern occurs when there is an effect of discharge on the concentration of a constituent as denoted by a non-zero slope. A negative slope (b < 0) represents a dilution behavior in which concentration varies inversely with discharge, while a positive slope (b > 0) represents an enrichment or mobilization behavior where a constituent's concentration increases with increasing discharge [12,16]. Oftentimes, there are changes in C-Q patterns, observable as inflections in slopes that further delineate the processes dominating the behavior of water quality constituents. These inflections in slopes, also known as discharge thresholds, reflect the interactions of hydrologic and biogeochemical forcing on the behavior of water quality constituents, whereby hydrology is more likely to dominate at high flow and biogeochemical processes at low flow [12]. Furthermore, the presence of slope inflections in C-Q relationships can identify the ranges of discharges in a flow regime where hydrology increases the concentration and export of constituents [12,15,20]. Theoretically, there are nine classifications (archetypes) of C-Q relationship patterns based on the presence of slope inflections ( Figure 1) [12]. Previous research has used the flow regime characteristics of channels to assess the impacts of flow changes caused by dams on river ecology [9,10,23,24], while few studies have addressed the effects of channelization [25]. This study expands on previous research by determining the impacts of flow regime changes on the behavior of water quality constituents in three channels in southern Florida, USA, whose channel morphologies and/or flow characteristics deviate from the flow regimes of natural channels. The objectives of this study were to (1) characterize and compare the flow regimes of a wetland slough, a river, and a flow-regulated, man-made canal; and (2) test the applicability of C-Q relationship models that have been traditionally used in channels with natural flow regimes to other channels that are different morphologically (channel) and hydrologically (flow regime). Predictions include the following: (1) the flow regime in the man-made canal will have the greatest variability compared to the slough and river; and (2) C-Q relationship models will be applicable in non-riverine channels to quantify the effect of flow on the behavior of water quality constituents.  [12]. Blue, red, and green represent dilution, chemostatic, and enrichment, respectively. Previous research has used the flow regime characteristics of channels to assess the impacts of flow changes caused by dams on river ecology [9,10,23,24], while few studies have addressed the effects of channelization [25]. This study expands on previous research by determining the impacts of flow regime changes on the behavior of water quality constituents in three channels in southern Florida, USA, whose channel morphologies and/or flow characteristics deviate from the flow regimes of natural channels. The objectives of this study were to (1) characterize and compare the flow regimes of a wetland slough, a river, and a flow-regulated, man-made canal; and (2) test the applicability of C-Q relationship models that have been traditionally used in channels with natural flow regimes to other channels that are different morphologically (channel) and hydrologically (flow regime). Predictions include the following: (1) the flow regime in the man-made canal will have the greatest variability compared to the slough and river; and (2) C-Q relationship models will be applicable in non-riverine channels to quantify the effect of flow on the behavior of water quality constituents.

Study Locations
Three channels were selected for this study: Shark River Slough, Peace River, and Hillsboro Canal, which are all located in southern Florida ( Figure 2). This region broadly has a sub-tropical to tropical climate, relatively flat topography and is underlain by carbonate (karst) formations that hold groundwater [26,27].
Shark River Slough is within Everglades National Park (ENP), a protected ecosystem comprising fresh-water wetlands, marshes, prairies, tidal flats, and mangroves [4]. The slough is a shallow, slow-moving, broad river, with channel depths in the order of centimeters [28] and is the primary pathway for surface-water flow through ENP [4]. Everglades National Park is underlain by porous marine limestone formations, including the Miami Limestone [29,30]. As part of a larger effort to drain the larger Everglades region for agriculture and urban development, discharges into Shark River Slough have been reduced by changes in the volume, timing, and distribution of water through regional canals, which diverted freshwater flows into the Atlantic Ocean [4].  [12]. Blue, red, and green represent dilution, chemostatic, and enrichment, respectively.

Study Locations
Three channels were selected for this study: Shark River Slough, Peace River, and Hillsboro Canal, which are all located in southern Florida ( Figure 2). This region broadly has a sub-tropical to tropical climate, relatively flat topography and is underlain by carbonate (karst) formations that hold groundwater [26,27].
Shark River Slough is within Everglades National Park (ENP), a protected ecosystem comprising fresh-water wetlands, marshes, prairies, tidal flats, and mangroves [4]. The slough is a shallow, slow-moving, broad river, with channel depths in the order of centimeters [28] and is the primary pathway for surface-water flow through ENP [4]. Everglades National Park is underlain by porous marine limestone formations, including the Miami Limestone [29,30]. As part of a larger effort to drain the larger Everglades region for agriculture and urban development, discharges into Shark River Slough have been reduced by changes in the volume, timing, and distribution of water through regional canals, which diverted freshwater flows into the Atlantic Ocean [4].
Peace River starts in south-central Florida and flows southwest until discharging into the Charlotte Harbor Estuary. The river drainage basin is 6086.47 sq km, with a length of about 168.98 km, and the basin geology consists of clastic sediments that overly carbonate rocks (Table 1) [31]. The river, although a natural watercourse, has had its flow regime modified by land use activities such as phosphate mining, urbanization, and agricultural development [32]. This has resulted in the loss and channelization of some of the headwater tributaries, and reduced surface inflows into the river [32]. However, much of the river corridor remains undeveloped; it has floodplains that support riparian vegetation, and the river flows unrestricted over its entire reach [32]. For example, the river reach at the town of Bartow, Florida, is within the upper river watershed (2139.33 sq km) and has abundant riparian vegetation and wetlands, including cypress and hardwood swamps [31]. Additionally, the channel reach is characterized by low bank heights and shallow bank slope angles [31] and is sinuous. Given these channel and hydrologic characteristics, Peace River most resembles a channel with a natural flow regime. corridor remains undeveloped; it has floodplains that support riparian vegetatio the river flows unrestricted over its entire reach [32]. For example, the river reach town of Bartow, Florida, is within the upper river watershed (2139.33 sq km) a abundant riparian vegetation and wetlands, including cypress and hardwood sw [31]. Additionally, the channel reach is characterized by low bank heights and s bank slope angles [31] and is sinuous. Given these channel and hydrologic characte Peace River most resembles a channel with a natural flow regime. The Hillsboro Canal is a major canal that originates from Lake Okeechobee and southeast, draining agricultural, wetland, and urban basins before discharging i Atlantic Ocean [8]. The geology of the canal basin is made up of the Fort Thomps mation, which consists of alternating layers of limestone, shell, sand, and marl pun by solution holes, and the southern portion grades into the more porous Miami Lim formation [30]. The Hillsboro Canal is one of the regional canals, constructed b 1910 and 1920, that was dredged to drain the Everglades and reroute freshwater nally headed south, eastward into the ocean [4]. Like all South Florida canals, the boro Canal lacks the natural channel features necessary to dissipate high-ener charges, such as floodplains, but rather has levees, and steep slopes to confine the the channel and prevent overbank discharge [5,8]. Discharge in the canal is manag water control structures to prevent flooding, drain excess water from farmlands, watershed runoff, recharge the surficial aquifer, and prevent salt-water intrusio The Hillsboro Canal is a major canal that originates from Lake Okeechobee and flows southeast, draining agricultural, wetland, and urban basins before discharging into the Atlantic Ocean [8]. The geology of the canal basin is made up of the Fort Thompson formation, which consists of alternating layers of limestone, shell, sand, and marl punctured by solution holes, and the southern portion grades into the more porous Miami Limestone formation [30]. The Hillsboro Canal is one of the regional canals, constructed between 1910 and 1920, that was dredged to drain the Everglades and reroute freshwater, originally headed south, eastward into the ocean [4]. Like all South Florida canals, the Hillsboro Canal lacks the natural channel features necessary to dissipate high-energy discharges, such as floodplains, but rather has levees, and steep slopes to confine the flow to the channel and prevent overbank discharge [5,8]. Discharge in the canal is managed via water control structures to prevent flooding, drain excess water from farmlands, convey watershed runoff, recharge the surficial aquifer, and prevent salt-water intrusion [5,8]. The channelization and flow regulation in the Hillsboro Canal makes it the most modified of the three channels in this study.

Data Collection
Discharge and water quality concentration data for fifteen years (2004-2018) were used in this study. The sites selected on the channels were Shark River Slough at Bottle Creek (SRS), Peace River at Bartow (PRB), and Hillsboro Canal at the S6 structure (HC6). Mean daily discharge data at SRS (25 •  , the water quality data were adjacent to the discharge station and were also retrieved from DBHYDRO.

Flow Regime Characterization
The discharge data in each channel were analyzed using the IHA program developed by the Nature Conservancy [10,34]. This program is free and user-friendly. The selected groups of hydrologic indices were the magnitude of annual flow conditions, frequency and duration of high and low pulses, and rate of change in water flow (see Table 2 for definitions). The Richards-Baker Flashiness Index (R-B Index), developed to determine the fluctuations in discharge relative to the total discharge [35], was calculated as an additional index of flow variability [36]. Because hydrologic data is usually not normally distributed, non-parametric statistics were used in the IHA program to calculate the values of the hydrologic indices and the median values were reported, except for the magnitudes, which were computed as means [34].
Principal component analysis (PCA) was used to identify the dominant hydrological indices that best described the flow regimes of the three channels. Principal component analysis is a multivariate analysis tool that is primarily used to reduce the dimensionality of large data by transforming the original variables in a dataset into new uncorrelated variables known as principal components, and the new variables or axes are linear combinations of the starting variables [37]. However, PCA can also be used to determine the relationships between variables and highlight the similarities and differences between categories [38].

Concentration-Discharge (C-Q) Relationships
The linear expression of the concentration-discharge (C-Q) relationship is written as follows: log 10 C = log 10 a + b log 10 Q For each channel, discharge and constituent concentration data were paired by date. Since flow is regulated at HC6, the discharge data used were restricted to pumping events when the canal was flowing. The data pairs were then log-transformed and regressed to yield slopes (b) [12], which were then used to understand the influence of flow on the concentrations of the water quality constituents. Student's t-test was conducted to ascertain the statistical significance of b being different from zero, with a level of significance at p < 0.05 [39,40]. p values lower than 0.05 indicated a significant log (C)-log (Q) slope; otherwise, b was not different from zero, and this meant that concentration was not discharge-dependent [40]. Next, piecewise (segmented) linear regressions were conducted on C-Q relationships to detect the presence of inflections in the slopes, which would signify changes in the behavior of the water quality constituents with discharge [20]. To do this, the Davies' test from the 'segmented' package in the R programming language [41] was used to iteratively search across 10 quantiles of the explanatory variable (Q) for breakpoints in the slope (b) of the log (C)-log (Q) regression for each constituent, with p < 0.05 selected as the level of significance [19,39]. When the significance level was met, breakpoint analysis from the package was conducted to estimate two new regression models, separated by an inflection in b. This indicated a change in the linear relationship between concentration and discharge [42,43]. An ANOVA test of independence was then conducted on the slopes above and below the inflection points to determine if they were statistically different from zero [39]. Flow duration curves, which are one of the outputs of the IHA program, were used to calculate the probabilities of exceeding the discharge thresholds for the constituents that had slope inflections. All statistical analyses were conducted in the R programming language [41].

Flow Regime Characterization
Flow regime characteristics differed between the three sites ( Table 2). The 1-day minimum discharge was highest at PRB (river), while the 90-day minimum, and 1-and 90-day maximum discharges, were highest at HC6 (canal). The channels SRS (slough) and PRB had low frequencies of high pulse counts while HC6 had the highest frequency of high pulse counts. SRS had the longest low-pulse duration (31 days) while PRB had the longest high-pulse duration (12 days), which was about twice (6.25 days) and four times (3 days) that of SRS and HC6, respectively. Conversely, HC6 had the highest hydrograph rise rate, fall rate, and highest flashiness (R-B Index) of the three channels. With consistently lower differences between low-flow and high-flow indices, the flow regime of SRS had the least variability, while HC6 had the greatest flow variability between the low-and high-flow indices and had the highest values across most hydrologic indices.
In the PCA, SRS and PRB were noticeably distinguished from HC6 ( Figure 3). For instance, SRS and PRB were differentiated from HC6 by low pulse duration and low pulse count. Additionally, PRB was differentiated by the fall rate, high pulse duration, 1-day minimum discharge, and 90-day minimum and maximum discharges, while HC6 was differentiated by R-B Index, high pulse count, 1-day maximum discharge and rise rate (rise rate is highly correlated with R-B Index and is thus omitted for clearer illustration) (Figure 3).

Concentration-Discharge (C-Q) Relationships across the Channels
The flows in the channels had quantifiable impacts on the concentrations of water quality constituents (Table 3). At SRS and PRB, dilution patterns tended to dominate with higher flows (negative b), particularly for TP, NN, SC, and Cl (PRB only). Conversely, turbidity and TSS were enriched with higher flows in PRB. At HC6, enrichment patterns were observed with higher flows (positive b) for TP, NN, and turbidity, and chemostatic patterns (b~0) were observed for TSS, SC, and Cl.
Piecewise regression analysis further revealed the presence of statistically significant slope inflections in the linear log (C)-log (Q) regressions for some water quality constituents across the channels ( Table 4). Five of the nine possible (C-Q) archetypes were found across SRS, PRB, and HC6 combined. At SRS, NN exhibited the CA archetype; that is, chemostasis at low flow and dilution at high flow, while SC exhibited the AA archetype signifying continual dilution with discharge ( Figure 4).  (a) (b)    At PRB, four of the five water quality constituents with statistically significant slope inflections (NN, turbidity, TSS, and SC) exhibited the BA archetype, which is enrichment at low flow and dilution at high flow ( Figure 5). The fourth constituent, Cl, exhibited the AA archetype.
At HC6, TP exhibited BB archetype-enrichment at both high and low flow, while NN and turbidity exhibited the CB archetype, which is chemostasis at low flow and enrichment at high flow (Table 4, Figure 6). The levels of significance for a change in slope according to the Davies' test were marginal for turbidity (0.05) ( Table 4). For SC, the archetype was CA and although the p value for the Davies' test was marginally significant at 0.04, ANOVA revealed that the slopes were not statistically different from zero. Furthermore, visual inspection of the C-Q plot shows that there is not really a definite pattern (Figure 6).

Impacts of Flow Regime Characteristics on Water Quality in the Channels
The longer low pulse duration at SRS translates into a longer residence time of low flow conditions where biogeochemical transformation can dominate water quality constituent behavior [12]. Furthermore, the low magnitude and flow variability at SRS, combined with the low pulses, resulted in chemostatic and dilution behavior, and the absence of enrichment behavior.
At PRB, a recurring behavior of the water quality constituents was an initial enrichment in concentrations at low flow followed by dilution at high flow (BA in Table 4, Figure  5). The flow regime characteristic most responsible for this is the high pulse duration. Previous studies suggest that low flow concentration enrichment can result from the abundance and proximity of widespread constituent pools that are easily exported and mobilized to (or from) the channels [16,19,39], while the switch to dilution at high flows is indicative of constituent exhaustion (source-limitation) [12,16,19]. Key reasons for the ex-

Impacts of Flow Regime Characteristics on Water Quality in the Channels
The longer low pulse duration at SRS translates into a longer residence time of low flow conditions where biogeochemical transformation can dominate water quality constituent behavior [12]. Furthermore, the low magnitude and flow variability at SRS, combined with the low pulses, resulted in chemostatic and dilution behavior, and the absence of enrichment behavior.
At PRB, a recurring behavior of the water quality constituents was an initial enrichment in concentrations at low flow followed by dilution at high flow (BA in Table 4, Figure 5). The flow regime characteristic most responsible for this is the high pulse duration. Previous studies suggest that low flow concentration enrichment can result from the abundance and proximity of widespread constituent pools that are easily exported and mobilized to (or from) the channels [16,19,39], while the switch to dilution at high flows is indicative of constituent exhaustion (source-limitation) [12,16,19]. Key reasons for the exhaustion of constituents at PRB are the extensive high flow duration, the presence of riparian vegetation, and gentle channel slopes. As a result, there is ample opportunity for extended deposition of water quality constituents from the channels into the floodplains [9,39]. Additionally, the high frequency of low pulse at PRB indicates more in-channel sediment deposition [23,24].
The channel HC6 had the greatest flow variability that resulted in the strong enrichment behavior at high flow. The unnatural and regulated discharges at HC6 resulted in higher flashiness and steeper hydrographs that contributed to the increases in the concentrations of water quality constituents. Flashier hydrologic responses have been linked to higher concentrations of nutrients, turbidity, and TSS, from channel erosion and sediment mobilization [25]. Results from a comparative study between a rural incised (channelized) stream and an urban unchanneled stream revealed a flashier hydrologic response and higher concentrations of water quality constituents in the rural stream, for all measured components except nitrate [25]. Just as storm events are periods of high sediment mobilization [20], so are canal discharge events, which explain the predominant enrichment behavior at HC6. The higher frequency of high pulse flows at HC6 translates into a greater export of sediment loads and the constituents they contain, such as TP. Furthermore, the absence of floodplains at HC6 means that a greater volume of water quality constituents and pollutant loads can be transported downstream into ecologically sensitive ecosystems.
The clear differentiation of HC6 from PRB and SRS, as shown by the PCA, alludes to the fundamental differences between natural and man-made channels. The high discharge maximums, pulse count, rise rate, and flashiness reflect the regulated flow management in the canal as well as the channel homogeneity. This contrasts with the slough and the river that have more natural heterogenous features and generally lower values for their flow regime characteristics.

Behavior of Water Quality Constituents at Low and High Flow Conditions across the Channels
All three channels had discharge thresholds for NN. The low flow chemostasis for NN at SRS and HC6 can be attributed to greater biogeochemical control on NN availability. Nitrogen species are biogeochemically active and can have high turnover rates [44]. This rapid cycling of NN between biotic and abiotic pools can regulate availability and can lead to a chemostatic response at low flow. As for most of the water quality constituents at PRB, NN exhibited a shift from low flow enrichment to high flow dilution. In Florida, NN-rich groundwater discharge that dominates at low flow conditions may be diluted with NN-poor temporary pools (rainfall and soil water) that are activated during high flow conditions, resulting in a shift from enrichment to dilution [19,20]. Although the Upper Peace River watershed mainly recharges the groundwater, high flow periods can raise the groundwater head higher than the river stage, which then discharges into the channel, and is facilitated by the porous kart formations (sinkholes and sand-filled depressions) [27]. At HC6, the low flow chemostasis of NN reflects the biogeochemical influence from the buildup of nutrient legacy stores that is typical of managed agricultural systems [20,22]. Conversely, the high flow enrichment indicates that higher discharges are responsible for the internal release and subsequent transport of NN in the canal.
For Cl at PRB, the greater dilution at low flow and the lesser dilution at high flow can be attributed to a point source contribution from the watershed. Dilution is common for Cl because it is not significantly altered by biogeochemical processes, making it a good conservative tracer [45]. Specific conductance followed the dominant PRB pattern of dilution, but the discharge threshold occurred at a lower value compared to the other constituents, suggesting a relatively scarce pool of the constituent that is quickly diluted. Although Shark River Slough is tidal, SRS is within the upstream freshwater zone, which accounts for the dilution of SC with increasing discharge. At HC6, the high variability of SC with discharge could result from the multiple natural sources of ions in the Everglades region, including trapped seawater in groundwater intersected by canals [46]. In addition, the complexities associated with canal discharge operations may have contributed to its high concentration variability.
Of the three channels, HC6 had a discharge threshold for TP, with moderate enrichment occurring at low flow and high enrichment at high flow. At low flow conditions, more phosphorus (P) can be retained in the channel and cycled physiochemically through sorption-desorption reactions and mineral precipitation-dissolution processes, and biologically through biotic uptake and decomposition of organic matter [13]. However, as flow increases, P can be transported through advection, which can cause entrainment of sediment-bound particulate P and desorption of P from suspended particles [13,19]. Decades of continuous nutrient loading from agriculture to South Florida canals, including the Hillsboro Canal, has led to the sequestration of P, often referred to as "legacy P" in canal sediments [47]. These sediments are highly organic and have low bulk density that makes them easily entrained and susceptible to downstream transport [48,49].
At PRB, turbidity and TSS had identical behavior and discharge thresholds, indicating the similarity between the two sediment-associated constituents. The low flow enrichment of TSS results from outer bank erosion and entrainment of particles from the channel bed [20], while the change to dilution at high flow indicates an exhaustion of the eroded sediments and their deposition in the floodplains [16]. At HC6, the strong enrichment behavior at high flow for turbidity results from transport limitation, whereby higher flows allow more connectivity between the watershed and the channel, and lead to greater mobilization and transport of sediments [12].

Limitations on the Use of C-Q Relationship Models of Man-Made and Regulated Channels
Unlike natural channels that may be free flowing, the man-made canals of South Florida do not flow freely; rather, their discharges are managed to meet the water needs of the region. As a result, canals exhibit extended periods of zero flow with stagnant conditions punctuated by short durations of high flows. Therefore, the C-Q relationship models could only explain the behavior of water quality constituents for periods where there was flow. This resulted in the flow exceedance probabilities of the water quality constituents with discharge thresholds at HC6 being consistently lower than the median discharge (0.5 flow exceedance probability). The median discharge is commonly used in natural channels, especially when analyzing riverine networks, as opposed to individual channels [12,19,40]. Consequently, median discharge may not always be the appropriate metric to use as the discharge threshold in C-Q slopes for singular channels or non-riverine channels. Furthermore, at HC6, a direct interpretation for the low discharge threshold probabilities (especially for TP and NN) is that the potential of canal flows mobilizing higher concentrations of these constituents is low. This means that the flows at HC6 that can export nutrients occur infrequently; however, the discharge thresholds and their associated probabilities will vary at other locations along the Hillsboro Canal and other South Florida canals depending on discharge operations. Therefore, C-Q relationships in such canals are useful in the determination of discharge that can limit the transport of nutrients to sensitive ecosystems, which can be a potential benefit to water resource managers.

Conclusions
This study described the relationships between water quality and flow in three morphologically and hydrologically different channels in southern Florida: Shark River Slough, Peace River, and Hillsboro Canal. The unchanneled and shallow Shark River Slough exhibited the least flow variability, which when combined with its naturally nutrient-poor conditions, resulted in predominantly dilution behavior of its water quality constituents at higher flows. The meandering Peace River, whose characteristics most resemble a natural flow regime, exhibited intermediate flow variability. Furthermore, Peace River had the highest pulse duration that afforded an extended temporal connectivity to its floodplains and resulted in dilution behavior of all the water quality constituents at higher flows. The heavily channelized Hillsboro Canal, whose flow regime is regulated by water control structures, had the greatest flow variability, and was dominated by enrichment behavior of its water quality constituents with higher flows. This study found that C-Q relationships can be applied to canals and should be considered by water managers in discharge operations to reduce the transport of nutrient and pollutant loads downstream to ecologically delicate ecosystems. Data Availability Statement: Publicly available datasets were analyzed in this study. These data can be found here: (https://www.sfwmd.gov/science-data/dbhydro, https://waterdata.usgs.gov/ nwis, https://fce-lter.fiu.edu/data/core/, https://www18.swfwmd.state.fl.us/ResData/Search/ ExtDefault.aspx) (accessed on 1 June 2021).