Advancing Discharge Ratings: A Novel Approach Based on Observed and Derivable GIS Factors in Alluvial Systems

: Depth–discharge rating is required at gauged and ungauged locations for hydrologic modeling of alluvial systems to evaluate streamﬂow and manage regional water resources. Spanning low to high-ﬂow conditions, manual ﬁeld measurements are used to develop discharge ratings at gauged locations, producing continuous ﬂow data from automated water depth measurements. The discharge rating is dependent on channel geometry, stream slope, vegetation, roughness coefﬁcient, sediment load, and bank stability. To construct discharge ratings for many locations within larger model domains (hundreds to thousands of km 2 ), intensive GIS and manual (spreadsheet) data manipulation are often required. In this analysis, available USGS gauging stations and readily available GIS coverages were used to learn and implement a novel method to characterize the depth–discharge relationships for hydrologic modeling of larger or complex areas using commonly available data and normalization techniques. The improved procedure simply uses drainage area, channel slope, and channel width, readily derivable GIS data, to develop discharge ratings for gauged and ungauged sections. The discharge rating curves for 70 USGS streamﬂow gauges were reproduced using the procedure. Then, the produced and observed discharge rating curves were compared to evaluate the accuracy of the method. In the analysis of streamﬂow depth predictions, the average Root Mean Squared Error was recorded at approximately 0.38 m ( ≈ 1.24 ft), with an interquartile range between 0.21 m and 0.49 m. The Mean Error remained centered around 0 m, with interquartile values ranging from − 0.24 m to 0.24 m.


Introduction
Managing water resources is challenging under the impacts of climate change and the ever-increasing demand from global population growth, which also leads to increased urbanization [1].Managing water resources depends on quantifying the available water resources and understanding the impacts of these and other stresses.In many cases, decisions concerning the best water management practice require using hydrologic simulation models [2].Hydrologic simulation models simplify complex real-world hydrology and estimate the hydrologic responses to alternative water resource management actions.These models can be used for extended periods of continuous simulation (years) or evaluating single storm event responses [3].The objective of developing hydrologic models can be to assess changes in streamflow due to changes to input stresses (e.g., climate variability, land use change, diversion options, mining effects, and stream-aquifer interaction from pumping stresses) to provide decision makers with the knowledge that enables better water resources management [4].Hydrologic models are conceptualized based on the model objectives (questions asked) and available budgets [5].
Water 2023, 15, 4152 2 of 20 Streamflow discharge is the quantity of water transported through a channel during a specific time, and the stage is the water surface height in a waterbody above a known datum [3].A depth-discharge relationship or discharge rating is a relationship between the flow rate and stream stage or water depth at a channel outlet [6].Within a hydrologic model, a discharge rating is required at the outlet of each simulation unit representing a portion of an alluvial wetland/stream system [7].For example, combined input of stage-discharge storage tabulated discharge rating is required for popular hydrologic models such as HSPF [8] and SWMM [9], as well as many others.
The hydrologic response within a model is very sensitive to discharge rating relationships, and these data are time-consuming for modelers to develop.The discharge rating depends on channel geometry, stream slope, vegetation, roughness coefficient, sediment load, and bank stability [10].Intensive GIS and manual (spreadsheet) data manipulation is required to construct discharge ratings for a model domain.This is especially problematic when the outlets of most stream reaches (i.e., model simulation units) do not have measured cross-sections, discharge ratings, or continuous discharge records.Moreover, the required work becomes daunting when hydrologic models are required for large areas with extensive alluvial stream and wetland systems.
In most applications, the available data to develop discharge ratings include GIS coverages and a few U.S. Geological Survey (USGS) streamflow gauging stations.The available stage-discharge data most likely originate from USGS stage gauging stations [11].Generally, USGS stage-discharge rating curves are available at established USGS streamflow measurement stations.However, they are relative to fixed local or National Geodetic Vertical Datum of 1929 (NGVD 29) or North American Vertical Datum of 1988 (NAVD 88) datums and can shift over time [12].These local or fixed datum stage-discharge rating relationships must be converted to depth-discharge ratings, considered more fixed over time.However, these data only represent one stream cross-section and do little to help populate discharge ratings when a large regional model may possess several hundred or thousands of stream sections that all require discharge ratings [4].
The standard approach to developing a stream's stage-discharge relationship or discharge rating is to visit the field several times to measure the flow for different stages, covering a series of low-to-high flows [13].These observations are then used to develop a curve that can convert the continuous stream-stage measurements to discharge estimates [13].Normally, to build a measured stage-discharge relationship in natural streams, observed stages must be plotted against the observed discharges [14].Alternatively, empirical methods to build discharge ratings can be used, such as assuming the streamflow is in uniform open channel flow with a prismatic channel section and then using Manning's or similar equations to predict or calculate flow estimates [3].These empirical methods require measured or estimated channel properties such as the cross-section area, slope, and subjectively interpreted roughness coefficients [4,13,15].Channel properties are seldom available online or from GIS databases.They are expensive and time-consuming to obtain, especially for larger regional assessments with hundreds to thousands of stream sections that must be parameterized with discharge rating tables.
Depth-discharge relationships are more consistent than stage-discharge in live beds and vegetated natural sections [6].They are thus more useful input for the discharge rating of alluvial stream/wetland systems [6].This study sought to develop a method to characterize the depth-discharge relationship for alluvial channels by using readily available GIS data, building on the work of Mueses et al. [6], Lewelling [16], Emmet [17], and Leopold [18].
As pointed out in the discussion above, discharge rating curves of natural streams tied to a fixed datum are not fixed; they shift over time.The changes in the live-bed channels and floodplains over time result in shifting discharge ratings in natural streams [19].Vegetation changes, fallen trees, live (variable) bed materials, and anthropogenic activities cause discharge rating curves to shift up or down over time.As a result, discharge rating data tied to the stage will not have a fixed singular one-to-one relationship, and a particular discharge threshold may be experienced for many different stages.Hydrologic models generally require a fixed time, one-to-one (depth-discharge pairs) discharge rating as an input.Mueses et al. [6] developed a simple model to produce an approximated fixed discharge rating curve using readily available GIS data.Mueses et al. [6] used 35 streamflow stations in southwest Florida to reproduce depth-discharge curves by normalizing the daily discharge and depth using the 10% discharge exceedance (Q 10 ), the corresponding depth (d 10 ), and a power law exponent (m Mueses ) to describe rating curve behavior.
In the expected application of the Mueses et al. [6] approach, Q 10 , d 10 , and m Mueses (exponent) values would be derived from GIS data, and the exponents would be found from nearby gauge settings or using average fitted exponents.Using the drainage area to predict Q 10 proved to be an efficient approach for predicting that variable.However, it was noted that there would be uncertain outcomes in other areas, and deriving Q 10 relationships could be regionally or locally dependent.However, searching for a simplistic method to predict d 10 using a similar analysis was not as productive.A method based on the square root of drainage area was hypothesized and tested but resulted in high errors and pore R 2 (R 2 ≤ 0.6) in depth prediction.There was another dilemma: how best to estimate the m Mueses exponent for ungauged sections?When analyzing gauged sites, each site was found to have a different exponent.It is suspected that each site (and thus each reach) possesses a different exponent based on the impact of the channel's shape, friction, and slope on the discharge rating behavior.
Lewelling [16] discussed one of the fundamentals of streamflow behavior that should be considered in this approach: discharge rating behavior is different for in-channel stages than out-of-bank flows.Generally, flows greater than 10% exceedance flows (less than 10% exceedance) begin to approach out-of-bank flows (flood stage).In comparison, the majority of the discharges (greater than 10% exceedance) are contained within the channel bank in natural and man-made sections [16].In addition, Leopold et al. [18] introduced the idea of normalizing the depth and discharge using the bankfull values to create a regional dimensionless depth-discharge relationship for ungauged locations.According to Mueses et al. [6], when their approach was compared to USGS observations, the results generally presented only small errors for intermediate and low discharges, with generally larger errors reserved for high (and presumably more out-of-bank) discharges.
It is well-recognized that natural channels exhibit different flow behavior (control) with characteristically different frictional conditions and geometries for in-bank and outof-bank flows [18].Thus, each control phase would exhibit different discharge rating behavior, especially the transition from channel control to floodplain control.Therefore, using a single exponent (m Mueses -value) to reproduce the depth-discharge curve for low and extreme discharges may be expected to yield relatively higher errors for extreme discharges, which was the hypothesized reason for the relatively high errors in larger flows in Mueses et al. [6].
Alluvial wetlands are those floodplain-type wetlands that encompass (surround) the main channel (thalweg), the storage of which greatly affects the out-of-bank flows through storage attenuation and differing floodplain frictional conditions [20].Alluvial wetlands, therefore, have incised channels representing the main thalweg and/or subchannels or branches in a dendritic stream system [21].These wetland channels are formed naturally and are in relation (proportion) to the dominant discharges strongly controlled by the cumulative drainage area, upland, and climatologic factors (e.g., soil types and mean annual rainfall).Thus, the channel size naturally increases with the high flow discharge index (e.g., Q 10 ) [17].Consequently, natural streams seem to achieve a balance between the channel geometry and the conveyance of the periodic but infrequent larger flows experienced in the channel [3].
Lawlor [22] investigated the relationship between the bankfull discharge, drainage area, and channel morphology to estimate the bankfull discharge at ungauged stations.The results showed that relating the drainage area to channel morphologic characteristics and the drainage area to the bankfull discharge using regression analysis can be useful to approximate channel morphology and/or bankfull discharge [22].Based on the findings of [17,23], the size of the drainage area within the same region can be reliably used to determine the bankfull channel discharge, bankfull mean depth, bankfull width, and bankfull cross-section area.
As pointed out previously, the only way to build a depth-discharge relationship in natural streams is by plotting computed depths vs. observed discharges [14].Plotting the depths against discharges on arithmetic paper over a single control phase or section yields a parabolic line, while plotting them via log-log paper or relations produces a straight line with a constant slope [14].Plotting a long period of depths and discharges in log-log relations yields straight lines with similar slopes but different offsets that enable hydrologists to dictate the change in control sections, and especially the transition from in-channel control to floodplain control, which is, by definition, beyond bankfull discharge.
The bankfull discharge of an incised channel is the dominant discharge controlling the evolution of the channel geomorphology.Bankfull flow is defined as the channel discharge rate just before the streamflow starts to overflow into floodplain wetlands [18].Periodic bankfull discharge dictates the channel's shape, size and invert and thus the depth model [14].In field measurements, the bankfull depth can often be recognized from stain lines or scour marks and thus measured from the channel bed, vegetation boundaries, transitions between the bed and bank materials, and/or sudden fluctuations on side slopes between channels and floodplains [14].As a result, the bankfull depth measurement is frequently subjectively interpreted or opinion-based and thus prone to error [14].
The studies on bankfull discharge estimation using reoccurrence interval, which is the exceedance probability of a bankfull discharge in a given year, are inconsistent.Lawlor [22] found that the bankfull discharge reoccurrence interval is between 1 year and 4.4 years, with a median of a 1.5-year recurrence threshold overall.Further supporting this threshold, previous studies showed that the recurrence interval was approximately 1.5-year exceedance [24][25][26][27][28].However, based on another study by Parrett and Johnson [29], the bankfull discharge at ungauged stations could be calculated through regression of the bankfull discharge to the two-year reoccurrence interval peak discharge.According to Lawlor [22] and previous studies, using a 1.5-year reoccurrence interval to calculate the bankfull discharge for ungauged stations might yield better estimates than the 2-year exceedance threshold.In addition, Mueses et al. [6] sought to develop a non-dimensional discharge rating model by normalizing with respect to bankfull discharge and depth; however, challenges in predicting the bankfull discharge and depth based on the drainage area rendered this approach unviable.
The main goal of this study was to develop a generic model and parameter guideline to produce discharge rating curves for vast ungauged alluvial wetland networks within a continuous hydrologic simulation model.The model seeks to present procedures and parameters for both in-channel and out-of-bank rating.This approach reduces the need for extensive field surveys.It requires only sufficiently resolved Digital Elevation Model(s) (DEMs); commonly available, delineated basins; and cumulative drainage areas.This method enables modelers to establish depth-discharge relationships for any area without direct gauge data, offering a cost-effective and adaptable solution for expansive terrains.

Data and Geographic Scope
USGS maintains streamflow gauging stations for numerous streams in southwest Florida.These stations are often monitored through cooperation between USGS and the local water management authorities, such as the Southwest Florida Water Management District (SWFWMD) (Figure 1).USGS measures daily stages for each gauge location relative to a particular datum, either locally established, such as NGVD 1929 or NAVD 88.Also, USGS staff visit gauge sites periodically to report field measurements such as the stage, cross-section area, vegetation condition, width, section control, and discharge.Discharge, stage, and field measurements were downloaded for 90 locations in southwest Florida, including 6 major rivers in the study area: Peace, Manatee, Myakka, Alafia, Withlacoochee, and Hillsborough.These data span from January 1980 to December 2022, covering a full range of flow conditions, including low, medium, and high flows.Some of these stations had no published stages or field measurements and/or were located in areas without DEM data.Data were used based on the availability of information from these stations.Therefore, some of these stations could not be used for all parts of this analysis or to reliably estimate depths or cross-section areas (please refer to Table 1 for more details).
Water 2023, 15, x FOR PEER REVIEW 5 of 21 Discharge, stage, and field measurements were downloaded for 90 locations in southwest Florida, including 6 major rivers in the study area: Peace, Manatee, Myakka, Alafia, Withlacoochee, and Hillsborough.These data span from January 1980 to December 2022, covering a full range of flow conditions, including low, medium, and high flows.Some of these stations had no published stages or field measurements and/or were located in areas without DEM data.Data were used based on the availability of information from these stations.Therefore, some of these stations could not be used for all parts of this analysis or to reliably estimate depths or cross-section areas (please refer to Table 1 for more details).Notes: Q10, Q50, and Q1 are 10th-, 50th-, and 1st-percentile exceedance discharges, respectively.A10 is the cross-section area at Q10, and d10, d50, and d1 are 10th-, 50th-, and 1st-percentile exceedance depths, respectively.AD is the drainage area, and SC is the channel slope.Notes: Q 10 , Q 50 , and Q 1 are 10th-, 50th-, and 1st-percentile exceedance discharges, respectively.A 10 is the crosssection area at Q 10 , and d 10 , d 50 , and d 1 are 10th-, 50th-, and 1st-percentile exceedance depths, respectively.A D is the drainage area, and S C is the channel slope.

Transforming Water Stage Measurements into Depth Values
For reference, USGS does not report the depth of a stream or flow station but monitors the stage and maintains a relationship between the measured stage and discharge.These stage-discharge records and field measurements were used to convert stage-discharge relationships to the depth-discharge relationships required for this analysis.Effective depth (d e ) of flow is relative to a derived depth that is representative of zero flow (even if not measured), which is termed the stream invert or zero flow depth.In the case of perennial stream flow, where no zero flow has been measured, effective depth must be derived from logarithmic extrapolation [6].The depth used in this analysis is the sum of the invert depth and the difference between the invert stage and the recorded stage.This methodology (referred to as effective depth) was previously presented in Mueses et al. [6].By calculating the difference between the estimated minimum flow stage (S 0 ) and the recorded stage (S) and adding the hydraulic mean depth for minimum discharge from the field measurement (d min ) (see also Figure 2), effective depth is where de is effective depth (L), d min is minimum hydraulic mean depth (L), S is stage (L), and S 0 is stage at minimum flow (L).It is important to note that the stage for d min and S 0 often changes with time.Recognizing this change, epics must be identified in the longer-term records to adjust the stage record with time [1,30].
Water 2023, 15, x FOR PEER REVIEW 6 of 21

Transforming Water Stage Measurements into Depth Values
For reference, USGS does not report the depth of a stream or flow station but monitors the stage and maintains a relationship between the measured stage and discharge.These stage-discharge records and field measurements were used to convert stage-discharge relationships to the depth-discharge relationships required for this analysis.Effective depth (de) of flow is relative to a derived depth that is representative of zero flow (even if not measured), which is termed the stream invert or zero flow depth.In the case of perennial stream flow, where no zero flow has been measured, effective depth must be derived from logarithmic extrapolation [6].The depth used in this analysis is the sum of the invert depth and the difference between the invert stage and the recorded stage.This methodology (referred to as effective depth) was previously presented in Mueses et al. [6].By calculating the difference between the estimated minimum flow stage (S0) and the recorded stage (S) and adding the hydraulic mean depth for minimum discharge from the field measurement (dmin) (see also Figure 2), effective depth is where de is effective depth (L), dmin is minimum hydraulic mean depth (L), S is stage (L), and S0 is stage at minimum flow (L).It is important to note that the stage for dmin and S0 often changes with time.Recognizing this change, epics must be identified in the longerterm records to adjust the stage record with time [1,30].

Improved Non-Dimensional Discharge Rating Curve
Mueses et al. [6] developed a nondimensional model to characterize the depth-discharge rating curve from the reference stage and discharge.Mueses et al. [6] normalized the discharge and depth using the 10th-percentile exceedance discharge, Q10, and the corresponding 10th-percentile depth, d10.Note that these exceedances are daily flows and stages and not annual exceedances.Thus, the 10th-percentile exceedance should be exceeded, on average, 36 days/year.The simple model input included Q10, d10, and mMueses which describe the relationship between discharges and depths related to the cross-sectional shape and frictional condition of the channel.Therefore, to develop a discharge rating for an ungauged simulation unit (waterbody), values of Q10, corresponding d10, and mMueses must be derived by empirical fitting techniques.The Mueses et al. [6] model is simply

Improved Non-Dimensional Discharge Rating Curve
Mueses et al. [6] developed a nondimensional model to characterize the depthdischarge rating curve from the reference stage and discharge.Mueses et al. [6] normalized the discharge and depth using the 10th-percentile exceedance discharge, Q 10 , and the corresponding 10th-percentile depth, d 10 .Note that these exceedances are daily flows and stages and not annual exceedances.Thus, the 10th-percentile exceedance should be exceeded, on average, 36 days/year.The simple model input included Q 10 , d 10 , and m Mueses which describe the relationship between discharges and depths related to the cross-sectional shape and frictional condition of the channel.Therefore, to develop a discharge rating for an ungauged simulation unit (waterbody), values of Q 10 , corresponding d 10 , and m Mueses must be derived by empirical fitting techniques.The Mueses et al. [6] model is simply where Q 10 is the 10th-percentile daily discharge exceedance (L 3 /T), d 10 is the 10th-percentile daily depth exceedance (L), and m Mueses is Mueses et al.'s [6] discharge behavior exponent.This method has been used with reasonable success to populate and calibrate large regional models in West-Central Florida [31].
M Mueses was calculated by finding the slope of a best-fit line through the logarithm of normalized depth and discharge values on an arithmetic plot.Thus, the method found the best fit log-linear equation, which is a power law expression of the form where b is the intercept.Using a construction method that involved the hydraulic depth of the lowest flow measurement to establish the depth at zero flow, even if it did not exist in the observations, the intercept (b) was eliminated.
A modified version of the Mueses et al. [6] non-dimensional model was formulated to better reproduce discharge rating curves in this study using GIS deliverable variables.The hypothesis behind the modified approach is that when the streamflow depth is within the channel bank (in-bank flows), the discharge rating curve follows one exponent.That exponent changes after the transition to out-of-bank discharges (floodplain stages).Smoothness in the discharge rating is ensured by using a common reference index.Breaks in discharge rating curves for flood flows (out-of-bank flows) have been well documented [14].The modified method applies two exponents.The first exponent reproduces the discharge rating curve for discharges less than Q 10 , and the second reproduces the discharge rating curve for discharges higher than Q 10 .To help adapt the method to be used for ungauged hydro-features, the model low-and medium-flow (in-bank) exponent (m 1 ) was calculated from the 50th-percentile exceedance flow and depth, Q 50 and d 50 , respectively, and the exponent for out-of-bank flows (m 2 ) was calculated from the 1st-percentile exceedance flow and depth, Q 1 and d 1 .
The relationships found for the exponents are For m 2 , the out-of-bank exponent is where Q 10 and d 10 were previously defined; Q 50 is the 50th-percentile daily exceedance discharge (L 3 /T); Q 1 is the 1st-percentile exceedance discharge (L 3 /T); d 50 and d 1 are the 50th-and 1st-percentiles exceedance effective depths (L), respectively; m 1 is the exponent for in-bank discharges; and m 2 is the exponent for out-of-bank discharges.
When the rating curves of USGS flow stations were studied, a trend was observed in these rating curves.Figure 3 shows four normalized rating curves of the studied gauge stations, each characterized by varying drainage areas ranging from 222 to 4325 km 2 .This figure is pivotal in highlighting the behaviors of m 1 and m 2 , demonstrating variations in the slope of the rating curve before and after the Q 10 (0,0) point.
gested that Q50 is the most reliably predictable low-flow percentile discharge.
Furthermore, Q1 was selected to characterize out-of-bank behavior.The methodology focused on using a flow metric that could both accurately predict and represent high-flow conditions.Attempts to model extremely high flows, specifically the 0.01st-percentile exceedance flow (Q0.01), resulted in an R 2 of 0.5.Consequently, Q1 was chosen as the representative metric for high flows.Therefore, Q10 was used to delineate the shift from in-bank behavior (m1) to out-of-bank behavior (m2), as further detailed in Figure 3. Normalized rating curves at 4 USGS gauge stations to demonstrate variations in m1 and m2 behavior centered around Q10 transition point (0,0).Note Q10 and d10 are 10th-percentile exceedance discharge (L 3 /T) and depth (L), respectively.m1 is the exponent for in-bank discharges, and m2 is the exponent for out-of-bank discharges.

Estimations of Q10, Q50, and Q1
The drainage area size, AD, and physical watershed characteristics are the key variables for regulating hydrologic fluxes within an environment with similar climatic and topographic conditions [17,23].While the AD is by far the most significant, other physical characteristics such as the watershed slope, wetland and upland vegetative cover, urban landform percentage (and imperviousness), soil type, and proximity of the water table to land surface play a lesser role in determining the high-flow (index) discharge rate within a single environment, especially for larger areas [1,30,32].This is because the variability of these influences for larger watersheds in many environments, including the study area in southwest Florida, is small.Mueses et al. [6] and Mueses [30] examined variables Observations revealed that Q 10 , typically not associated with out-of-bank discharge, remained within bank limits for these streams.A closer examination of the rating curve showed that Q 10 is a critical transition point for the onset of out-of-bank behavior.Mathematically, Q 10 serves as an ideal marker for identifying the beginning of out-of-bank flow, primarily due to its position at the lower end of this transition phase.
This study also utilized Q 50 to represent both mean and minimum flows.As Q 50 is the median flow, it was deemed a reliable indicator.However, the accuracy of low-flow measurements was often compromised by factors like deposition, erosion, or channel obstructions, which could alter the rating curve and affect measurement precision.It was suggested that Q 50 is the most reliably predictable low-flow percentile discharge.
Furthermore, Q 1 was selected to characterize out-of-bank behavior.The methodology focused on using a flow metric that could both accurately predict and represent high-flow conditions.Attempts to model extremely high flows, specifically the 0.01st-percentile exceedance flow (Q 0.01 ), resulted in an R 2 of 0.5.Consequently, Q 1 was chosen as the representative metric for high flows.Therefore, Q 10 was used to delineate the shift from in-bank behavior (m 1 ) to out-of-bank behavior (m 2 ), as further detailed in Figure 3.The drainage area size, A D , and physical watershed characteristics are the key variables for regulating hydrologic fluxes within an environment with similar climatic and topographic conditions [17,23].While the A D is by far the most significant, other physical characteristics such as the watershed slope, wetland and upland vegetative cover, urban landform percentage (and imperviousness), soil type, and proximity of the water table to land surface play a lesser role in determining the high-flow (index) discharge rate within a single environment, especially for larger areas [1,30,32].This is because the variability of these influences for larger watersheds in many environments, including the study area in southwest Florida, is small.Mueses et al. [6] and Mueses [30] examined variables derivable from GIS coverages influencing Q 10 .The conclusion was A D is the only parameter showing a strong correlation.
The purpose of this simplified analysis was to identify relationships that can estimate these reference discharges for ungauged hydro-features (where no discharge rating data Water 2023, 15, 4152 9 of 20 exist) and in vast hydrography settings, typical of larger regional models.It is important to point out that these exceedances follow a non-normal distribution and are skewed in generally predictable ways [33].For these data, the accuracy of parametric tests, which measure the fit of linear regression models, is invariably poor [33].The normality of these reference discharges was tested graphically using a Q-Q plot [33], and the results show the expected non-normality distribution (Figure 4).Therefore, the discharge datasets were logarithmically transformed (log Q values) to make the distribution symmetrical and normal (Figure 5).Then, the Q 1 , Q 50 , and Q 10 for 90 gauge stations were randomly divided into training and test datasets, with a heavier weight toward the training dataset (80% of the data), while the test dataset contained the remaining 20%.The training dataset was linearly regressed, and the resultant estimators were tested using parametric tests against the test dataset to evaluate the fit quality.The results were then evaluated using R 2 , which is the variation of the discharges that the drainage area can explain.
derivable from GIS coverages influencing Q10.The conclusion was AD is the only p ter showing a strong correlation.
The purpose of this simplified analysis was to identify relationships that can e these reference discharges for ungauged hydro-features (where no discharge rati exist) and in vast hydrography settings, typical of larger regional models.It is im to point out that these exceedances follow a non-normal distribution and are ske generally predictable ways [33].For these data, the accuracy of parametric tests measure the fit of linear regression models, is invariably poor [33].The normality reference discharges was tested graphically using a Q-Q plot [33], and the resul the expected non-normality distribution (Figure 4).Therefore, the discharge datase logarithmically transformed (log Q values) to make the distribution symmetrical a mal (Figure 5).Then, the Q1, Q50, and Q10 for 90 gauge stations were randomly d into training and test datasets, with a heavier weight toward the training dataset the data), while the test dataset contained the remaining 20%.The training data linearly regressed, and the resultant estimators were tested using parametric tests the test dataset to evaluate the fit quality.The results were then evaluated using R 2 is the variation of the discharges that the drainage area can explain.

Estimation of A 10
As pointed out previously, USGS visits streamflow gauges for low, medium, and high discharges multiple times annually and takes field measurements such as velocity, stage, and cross-section area.These data are generally published online as field notes for the gauging stations.The time series of cross-section areas from 70 streamflow gauge stations with suitable high-resolution DEM data in southwest Florida were used to derive a relationship for estimating the cross-section area at Q 10 , A 10 , for each station.
SWFWMD produced high-resolution 1.5 m horizontal, vertical accuracy < 0.1 m (5 ft horizontal, vertical accuracy < 0.3 ft) DEM for southwest Florida regions to estimate the local longitudinal channel slopes, S c , and channel width, W C .The S c of 70 gauge stations was estimated by determining the elevation at the gauge station and approximately a thousand meters upstream of the gauge station.W C was also calculated by measuring the width of the channel from the DEM using a GIS measure tool at the gauge locations.

Estimation of A10
As pointed out previously, USGS visits streamflow gauges for low, medium, and high discharges multiple times annually and takes field measurements such as velocity, stage, and cross-section area.These data are generally published online as field notes for the gauging stations.The time series of cross-section areas from 70 streamflow gauge stations with suitable high-resolution DEM data in southwest Florida were used to derive a relationship for estimating the cross-section area at Q10, A10, for each station.
SWFWMD produced high-resolution 1.5 m horizontal, vertical accuracy < 0.1 m (5 ft horizontal, vertical accuracy < 0.3 ft) DEM for southwest Florida regions to estimate the local longitudinal channel slopes, Sc, and channel width, WC.The Sc of 70 gauge stations was estimated by determining the elevation at the gauge station and approximately a thousand meters upstream of the gauge station.WC was also calculated by measuring the width of the channel from the DEM using a GIS measure tool at the gauge locations.
The AD and the Sc 0.5 for these streamflow gauging stations were logarithmically transformed because of non-normality (Figures 4 and 5) and divided into training and test datasets, with 80% allocated to the training set and 20% to the test set.Then, the training data were used in a multiple-variables linear regression model to estimate the A10.Then, the model was tested against the test dataset, again by regression.The multiple-variables linear regression model was evaluated using R 2 .

Estimations of d10, d50, and d1
Three linear regression models were developed: one to predict d10 using the hydraulic mean depth at Q10, dh10, and the other two models to predict d1 and d50 using d10.dh10 for 83 streamflow gauge stations was calculated from observed discharge rating curves by dividing A10 (L 2 ) by WC (L) and then regressed to estimate d10 (L).d10 was estimated using the Mueses et al. [6] method for effective depth.The stage at Q10 was subtracted from the The A D and the S c 0.5 for these streamflow gauging stations were logarithmically transformed because of non-normality (Figures 4 and 5) and divided into training and test datasets, with 80% allocated to the training set and 20% to the test set.Then, the training data were used in a multiple-variables linear regression model to estimate the A 10 .Then, the model was tested against the test dataset, again by regression.The multiple-variables linear regression model was evaluated using R 2 .

Estimations of d 10 , d 50 , and d 1
Three linear regression models were developed: one to predict d 10 using the hydraulic mean depth at Q 10 , d h10 , and the other two models to predict d 1 and d 50 using d 10 .d h10 for 83 streamflow gauge stations was calculated from observed discharge rating curves by dividing A 10 (L 2 ) by W C (L) and then regressed to estimate d 10 (L).d 10 was estimated using the Mueses et al. [6] method for effective depth.The stage at Q 10 was subtracted from the minimum flow stage and added to the hydraulic mean depth of the lowest field measurement.A similar method was employed to estimate d 1 and d 50 .Then, the data were divided into training and test datasets of 80% and 20%, respectively.Finally, d 10 for 66 of the 83 streamflow gauge stations was regressed to estimate d 1 and d 50 .Then, the results were evaluated against the test dataset.The regression model accuracy of predicting d1 and d 50 was evaluated using R 2 .
where d h10 is the hydraulic-mean depth of Q 10 (L), A 10 is the cross-section area of Q 10 (L 2 ), and W C is the width of Q 10 (L).

m Mueses vs. m 1 (in Channel Flow Exponent)
As discussed earlier, a single Mueses' exponent performed poorly for discharges higher than Q 10 .Therefore, at a minimum, there is a need to use two exponents to describe the full range of discharge ratings, including discharges below and above Q 10 (in-and out-of-bank flows).The required work to generate a m Mueses value is only possible with extensive flow measurement data.Hence, as the method is desired to populate the vast ungauged hydrography in a regional model domain, the m 1 , derived from Equation ( 4) using observed Q 10 , Q 50 , d 10 , and d 50 , was compared against m Mueses .Then, m Mueses and m 1 were evaluated using R 2 .On the other hand, the exponent for discharges higher than Q 10 , m 2 , was produced using Equation ( 6) but was not used in the comparison since m 2 was introduced within this study.

Testing and Verification: Addressing Uncertainty
Regression models were formulated to estimate discharge rating parameters using three readily derivable variables for any ungauged simulation unit: A D , S c , and W C .These determinants were applied within the regression models to generate predictions for Q 10 , Q 50 , Q 1 , d 10 , d 50 , and d 1 .Following this, m 1 and m 2 were derived from Equations ( 4) and ( 6).The precision of these forecasts was gauged via the Root Mean Squared Error (RMSE) and Mean Error (ME).Subsequently, these predictions were employed in Equations ( 5) and ( 7) to reproduce the USGS discharge rating curves.Then, the replicated discharge rating curves were compared to the observed USGS discharge rating curves.Also, the original model from Mueses et al. [6] for ungauged simulation units was compared to the observed USGS discharge rating curves.Then, the accuracy was evaluated through RMSE and ME.

Predicting Q 10 , Q 50 , and Q 1
Predicting the exceedance discharges for Q 50 , Q 10 , and Q 1 using only A D showcased strong performance.Specifically, R 2 values for the training dataset were 0.78, 0.9, and 0.84 for Q 50 , Q 10 , and Q 1 , respectively.For the test dataset, the R 2 values for Q 50 , Q 10 , and Q 1 were statistically comparable, with values of 0.76, 0.89, and 0.81, respectively (see Figure 6).The derived linear equations are as follows: Log(Q 1 ) = 0.67 Log(A D ) − 0.32 (11) where Q 50 , Q 10 , and Q 1 were previously defined (m 3 /s), and A D is the drainage area in (km 2 ).

Predicting A 10
The multi-variable linear regression model using A D and S c 0.5 demonstrated a robust relationship for the A 10 .For the training and test datasets, the R 2 values were 0.87 and 0.92, respectively.Notably, the p-values for A D and S c 0.5 were less than 0.0000.A 10 exhibited a positive relationship with A D and a negative relationship with S c 0.5 (see Figure 7).The results indicated that, statistically, A D and S c 0.5 are reliable estimators for A 10 , as Log(A 10 ) = 0.69 Log(A D ) − 0.72 Log S c 0.5 − 1.62 (12) where A 10 (m 2 ) and A D (km 2 ) were previously defined, and S c 0.5 is the channel slope (m/m).

Predicting A10
The multi-variable linear regression model using AD and Sc 0.5 demonstrated a robust relationship for the A10.For the training and test datasets, the R 2 values were 0.87 and 0.92, respectively.Notably, the p-values for AD and Sc 0.5 were less than 0.0000.A10 exhibited a positive relationship with AD and a negative relationship with Sc 0.5 (see Figure 7).The results indicated that, statistically, AD and Sc 0.5 are reliable estimators for A10, as ( 10 ) = 0.69 (  ) − 0.72 (  0.5 ) − 1.62 (12) where A10 (m 2 ) and AD (km 2 ) were previously defined, and Sc 0.5 is the channel slope (m/m).

Predicting d 10 , d 50 , and d 1
The model's ability to predict the different percentile exceedance depths was evaluated.The R 2 for the training datasets of d 10 , d 50 , and d 1 were 0.85, 0.77, and 0.88, respectively, and R 2 for the test datasets for these same variables was substantially equivalent at 0.87, Water 2023, 15, 4152 13 of 20 0.78, and 0.9, respectively (see Figure 8).The results conclusively showed that d h10 can be used to estimate d 10 , d 50 , and d 1 , and the resulting equations (for units of meter) are Log(d 10 ) = 0.84 Log(d h10 ) + 0.16 ( 13) Log(d 1 ) = 0.86 Log(d 10 ) + 0.21 (15) where d h10 , d 10 , d 50 , and d 1 were previously defined (m).

mMueses vs. m1 (in Channel Flow Exponent)
The exponent, mMueses, proposed by Mueses et al. [6], represents the best-fitted line through normalized depths and discharges of the discharge rating curves.When linearly regressed against the in-bank exponent of the refined model, m1, calculated from Q10, Q50, d10, and d50, it yielded an R 2 value and slope close to 1 (≈45°), as depicted in Figure 9.

m Mueses vs. m 1 (in Channel Flow Exponent)
The exponent, m Mueses , proposed by Mueses et al. [6], represents the best-fitted line through normalized depths and discharges of the discharge rating curves.When linearly regressed against the in-bank exponent of the refined model, m 1 , calculated from Q 10 , Q 50 , d 10 , and d 50 , it yielded an R 2 value and slope close to 1 (≈45 • ), as depicted in Figure 9.

Assessing Uncertainty in Parameter Estimation and Rating Curve Reproduction
When the AD of 70 gauge stations was used in the proposed linear regression models to predict Q10, Q50, and Q1, the RMSE values were 4.98 m 3 /s, 1.04 m 3 /s, and 13.37 m 3 /s, respectively.While the ME exhibited small bias in all predicted discharges, the ME values were −1.38 m 3 /s, −0.35 m 3 /s, and −4.21 m 3 /s, respectively.Also, when AD, Sc 0.5 , and WC were used to predict A10, d10, d50, and d1, the RMSE values were 18.81 m 2 , 0.64 m, 0.34 m, and 0.88 m, respectively.While the ME exhibited a small bias as well, the values were −2.75 m 2 , −0.02 m, −0.01 m, and −0.04 m, respectively.The observed bias in the results was expected, given the diverse streamflow gauge stations utilized to predict each parameter within the linear regression models.These predictions were subsequently tested across 70 streamflow gauge stations, some of which were part of the training dataset while others were not (please refer to Table 2 for more details).Compared with the Mueses et al.'s [6] model, the improved model performance was superior across different RMSE ranges.For an RMSE of ≤0.3 m (≈≤1 ft), the improved model showed that 40% of the stations (28 out of 70) achieved this accuracy, compared to 25.7% (18 out of 70 stations) in the Mueses et al. [6] model.In the RMSE range of 0.3 to 0.6 m (≈1 to 2 ft), 50% of stations fell under the improved model, whereas only 32.86% did so under the Mueses model.In the improved model, 90% of the stations achieved an RMSE ≤0.6 m (≤2 ft), compared to only 58% of stations achieving this range with the Mueses

Assessing Uncertainty in Parameter Estimation and Rating Curve Reproduction
When the A D of 70 gauge stations was used in the proposed linear regression models to predict Q 10 , Q 50 , and Q 1 , the RMSE values were 4.98 m 3 /s, 1.04 m 3 /s, and 13.37 m 3 /s, respectively.While the ME exhibited small bias in all predicted discharges, the ME values were −1.38 m 3 /s, −0.35 m 3 /s, and −4.21 m 3 /s, respectively.Also, when A D , S c 0.5 , and W C were used to predict A 10 , d 10 , d 50 , and d 1 , the RMSE values were 18.81 m 2 , 0.64 m, 0.34 m, and 0.88 m, respectively.While the ME exhibited a small bias as well, the values were −2.75 m 2 , −0.02 m, −0.01 m, and −0.04 m, respectively.The observed bias in the results was expected, given the diverse streamflow gauge stations utilized to predict each parameter within the linear regression models.These predictions were subsequently tested across 70 streamflow gauge stations, some of which were part of the training dataset while others were not (please refer to Table 2 for more details).Compared with the Mueses et al.'s [6] model, the improved model performance was superior across different RMSE ranges.For an RMSE of ≤0.3 m (≈≤1 ft), the improved model showed that 40% of the stations (28 out of 70) achieved this accuracy, compared to 25.7% (18 out of 70 stations) in the Mueses et al. [6] model.In the RMSE range of 0.3 to 0.6 m (≈1 to 2 ft), 50% of stations fell under the improved model, whereas only 32.86% did so under the Mueses model.In the improved model, 90% of the stations achieved an RMSE ≤ 0.6 m (≤2 ft), compared to only 58% of stations achieving this range with the Mueses model.A more detailed performance analysis for higher RMSE ranges, along with the corresponding percentages of stations, is outlined in Table 2 and illustrated in Figure 10.The average RMSE was notably lower in the improved model at 0.38 m (≈1.24 ft), compared to 0.67 m (≈2.2 ft) in the Mueses model, reflecting a significant enhancement in prediction accuracy.For a detailed comparison of these models across different RMSE ranges, refer to Table 3 and Figure 10.

RMSE (m)
USGS Percentage of Stations within RMSE Range (% Improved Model Mueses Mo 0 to 0.3 (≈0 to 1 ft) 40% 25.71% 0.3 to 0.6 (≈1 to 2 ft) 50% 32.86% 0.6 to 0.9 (≈2 to 3 ft) 7.14% 24.29% 0.9 to 1.2 (≈3 to 4 ft) 1.43% 5.71% 1.2 to 1.5 (≈4 to 5 ft) 1.43% 2.86% 1.5 to 1.8 (≈5 to 6 ft) 0% 1.43% 1.8 to 2.1 (≈6 to 7 ft) 0% 1.43% The discharge rating curves generated by the improved model (d′), Mueses e model (dMueses), and the observed data from USGS streamflow gauge stations (d) we ted for three different locations, each characterized by a unique AD.The improved demonstrated robust performance in simulating observed discharge rating curves and high-flow conditions at all three locations.This enhanced accuracy can be attr to applying dual exponents and integrating both Sc and Wc within the model fram On the other hand, the Mueses model, which predominantly relies on drainage the predictor, was found to underestimate the observed discharges at AD of 530 km ure 11c) generally.At a larger drainage area of 3540 km 2 , the model tended to overes the observed behavior, underscoring the limitations inherent in using only AD as dependent variable for developing discharge rating curves.For further details, r Figure 11.The discharge rating curves generated by the improved model (d ), Mueses et al. [6] model (d Mueses ), and the observed data from USGS streamflow gauge stations (d) were plotted for three different locations, each characterized by a unique A D .The improved model demonstrated robust performance in simulating observed discharge rating curves for low and high-flow conditions at all three locations.This enhanced accuracy can be attributed to applying dual exponents and integrating both Sc and Wc within the model framework.On the other hand, the Mueses model, which predominantly relies on drainage area as the predictor, was found to underestimate the observed discharges at A D of 530 km 2 (Figure 11c) generally.At a larger drainage area of 3540 km 2 , the model tended to overestimate the observed behavior, underscoring the limitations inherent in using only A D as the independent variable for developing discharge rating curves.For further details, refer to Figure 11.

Discussion, Recommendations, and Limitations
This study's objective was to develop a more reliable method to characterize the depth-discharge relationship for ungauged alluvial stream/wetland systems.Using commonly available GIS data and simple relationships, the method has been tested for southwest Florida hydrology and demonstrated improved performance compared to the earlier procedures presented by Mueses et al. [6].The method uses minimal calibration parameters and can be used or adapted for other environments.Since the drainage area size is the most important variable and previous studies indicated much lower sensitivity to other variables such as upland land use and watershed slope, it is believed that this variable, together with a GIS determination of longitudinal channel slope and transverse width, can be used to estimate the important controlling factors of A10, d10, Q10, and in-bank and

Discussion, Recommendations, and Limitations
This study's objective was to develop a more reliable method to characterize the depthdischarge relationship for ungauged alluvial stream/wetland systems.Using commonly available GIS data and simple relationships, the method has been tested for southwest Florida hydrology and demonstrated improved performance compared to the earlier procedures presented by Mueses et al. [6].The method uses minimal calibration parameters and can be used or adapted for other environments.Since the drainage area size is the most important variable and previous studies indicated much lower sensitivity to other variables such as upland land use and watershed slope, it is believed that this variable, together with a GIS determination of longitudinal channel slope and transverse width, can be used to estimate the important controlling factors of A 10 , d 10 , Q 10 , and in-bank and out-of-bank discharge rating exponents, m 1 and m 2 , for applications where measured cross-sections and discharge rating conditions are sparse.
Providing separate exponents for in-bank and out-of-bank discharge ratings greatly increases the accuracy and physical basis of the approach.However, for ungauged sections, doubling the number of uncertain parameters would lead to overparameterization if an objective means to define both were unavailable.This approach indicated that adding another flow index, Q 50 , could be used to help predict the in-bank flow exponent, m 1 .Since Q 50 was found to be strongly correlated to drainage area and, by definition, represents the median flow behavior, it was shown to suitably provide for the determination of the low-flow (in-bank) discharges and the appropriate exponent, m 1 .
For higher flows, Q 1 was also shown to be strongly correlated to the drainage area, and this flow represents a likely out-of-bank or near-bankfull flow condition.This flow index has a lower probability of occurrence and is observed for approximately only 3 days per year on average.Since it is appropriate for high flows and for convenience, it was selected to provide a method to estimate the high-flow behavior, m 2 , and showed acceptable predictive performance.However, where data exist for out-bank behavior, discharge exceedances with extremely low probability, such as the 0.01 percentile conditions, can be utilized.It should be noted that these percentile exceedances refer to daily, rather than yearly, occurrences.
Not surprisingly, the 90-percentile exceedance discharge was not suitably dependent on the drainage area (Figure 12), perhaps due to the strong sensitivity to baseflow and groundwater processes, and thus provided an index that was not useful for this method.Furthermore, natural streams experience live bed processes, shifting bed materials, and sensitivity to fallen debris, making the discharge rating curve behavior for low flows unreliable or even unstable.As a result, Q 50 is believed to be the lowest index that is suitable for describing the medium-and low-flow discharge rating behavior.
Water 2023, 15, x FOR PEER REVIEW 18 of 21 out-of-bank discharge rating exponents, m1 and m2, for applications where measured cross-sections and discharge rating conditions are sparse.Providing separate exponents for in-bank and out-of-bank discharge ratings greatly increases the accuracy and physical basis of the approach.However, for ungauged sections, doubling the number of uncertain parameters would lead to overparameterization if an objective means to define both were unavailable.This approach indicated that adding another flow index, Q50, could be used to help predict the in-bank flow exponent, m1.Since Q50 was found to be strongly correlated to drainage area and, by definition, represents the median flow behavior, it was shown to suitably provide for the determination of the lowflow (in-bank) discharges and the appropriate exponent, m1.
For higher flows, Q1 was also shown to be strongly correlated to the drainage area, and this flow represents a likely out-of-bank or near-bankfull flow condition.This flow index has a lower probability of occurrence and is observed for approximately only 3 days per year on average.Since it is appropriate for high flows and for convenience, it was selected to provide a method to estimate the high-flow behavior, m2, and showed acceptable predictive performance.However, where data exist for out-bank behavior, discharge exceedances with extremely low probability, such as the 0.01 percentile conditions, can be utilized.It should be noted that these percentile exceedances refer to daily, rather than yearly, occurrences.
Not surprisingly, the 90-percentile exceedance discharge was not suitably dependent on the drainage area (Figure 12), perhaps due to the strong sensitivity to baseflow and groundwater processes, and thus provided an index that was not useful for this method.Furthermore, natural streams experience live bed processes, shifting bed materials, and sensitivity to fallen debris, making the discharge rating curve behavior for low flows unreliable or even unstable.As a result, Q50 is believed to be the lowest index that is suitable for describing the medium-and low-flow discharge rating behavior.The southwest Florida region is mostly characterized as a low-slope environment (generally <0.3%) with protected alluvial floodplain wetlands, where natural streams tend toward uniform flow more often than higher-gradient systems.This may be one reason the concept of a fixed and predictable discharge rating behavior prevails and seems appropriate.
Relating the drainage area to depth assumes that when the drainage area is increasing, the depth is also, which is not always the case considering channel width variability.However, switching to a cross-sectional area dependence (A10), assuming that the crosssectional area is directly related to the magnitude of Q10, is conceptually sound and has The southwest Florida region is mostly characterized as a low-slope environment (generally <0.3%) with protected alluvial floodplain wetlands, where natural streams tend toward uniform flow more often than higher-gradient systems.This may be one reason the concept of a fixed and predictable discharge rating behavior prevails and seems appropriate.
Relating the drainage area to depth assumes that when the drainage area is increasing, the depth is also, which is not always the case considering channel width variability.However, switching to a cross-sectional area dependence (A 10 ), assuming that the crosssectional area is directly related to the magnitude of Q 10 , is conceptually sound and has been shown to be a better basis for defining d 10 .Hence, GIS procedures or other interpretations will be needed to define the hydraulic (or top-of-bank) width for Q 10 .
The primary limitation of this study is that the empirical relationships may differ in other environmental settings, particularly in highly variable alluvial systems (e.g., an-thropogenetic and hardened channel settings).This is because it can be challenging to ascertain whether channels in such settings are fully developed and contain any relation to flow return probability (index).The difficulty lies in determining if these channels have evolved between regular flow states and the erosion susceptibility of channel morphology.In sandy-terrain Florida, channels are highly variable in degree of incision, some deeply incised, others not, but are in relative dynamic equilibrium with probabilistic flows and cross-section relation, leading to a strong relationship between drainage area and crosssectional area.This would likely not be the case in artificially hardened or natural rock-bed channels in non-depositional environments.However, even with different soil types, land use, and hydrology in coastal plain environments, these empirical relationships will likely still be applicable, but some local adjustments to the coefficients may be necessary.To apply this methodology in a new environment, it will be essential to have a sufficient number of gauge stations with extensive records of in-bank and out-of-bank discharges to validate the approach and coefficient behavior.Additionally, observational verification of the method in different settings may prove to indicate yet another variable or process that becomes important.In the absence of a sufficiently resolved DEM, channel width and slope must be manually determined from topographic data or estimated by the user.

Conclusions
This study sought to develop a more reliable method for characterizing the depthdischarge relationship in dendritic alluvial stream/wetland systems for hydrologic modeling where limited or sparse gauge or survey data are available.The method uses simple GIS and normalized analyses based on derivable flow indexes defined for each channel section.The improved procedure simply requires cumulative drainage area, local channel slope, channel width, and readily derivable GIS data to develop discharge ratings for gauged and ungauged sections.The discharge rating curves for 70 USGS streamflow gauges were reproduced using the procedure, and then the produced and observed discharge rating curves were compared to evaluate the accuracy of the method.The methodology developed demonstrated improved performance compared to previous methods with minimal calibration parameters.
The comparison showed robust performance in predicting depth with minimal errors.In the analysis of streamflow depth predictions, the average was 0.38 m (≈1.24 ft), with an interquartile range between 0.21 m and 0.49 m.The ME showed no bias (centered around 0 m), with manageable interquartile values ranging +/−0.24 m.
The approach relied heavily on the cumulative drainage area to each reach (A D ), which was found to be the most significant variable.By combining GIS-determined channel slopes (S c ) and width (W C ), as well as key flow and depth indexes of d 10 , d 50 , and d 1 , and Q 10 , Q 50 , and Q 1 , discharge rating exponents m 1 and m 2 can be estimated, which were shown to be reliable for predicting in-bank and flood flow (out-of-bank) discharge ratings.This study found strong relationships between Q 50 and A D for low-flow discharges and between Q 1 and A D for high flows, aiding in estimating the respective exponents m 1 and m 2 .
The success of the approach in the studied Florida region, characterized by similar vegetation and slopes in fully developed (depositional equilibrium) floodplain wetlands, suggests potential applicability in other coastal plain environments.Adjustments may be necessary in different environmental settings, particularly where channels have high variability or are not in full depositional equilibrium.However, success and confidence in other settings will require sufficient gauging stations and adequate GIS data before validating the approach.
Funding: The authors gratefully acknowledge the funding support provided by Tampa Bay Water, under grant number 2014/1325/00, and project manager Jeffrey Geurink.Geurink provided invaluable expertise and technical and financial assistance in this research.
Data Availability Statement: Please contact the corresponding author for data.

Figure 1 .
Figure 1.Map of USGS streamflow gauge stations across southwest Florida.

Figure 1 .
Figure 1.Map of USGS streamflow gauge stations across southwest Florida.

Figure 2 .
Figure 2. Definition of effective depth, de, from measured stage.S is the measured stage (L), S0 is the minimum measures stage (L), and dmin is the hydraulic mean depth of S0 (L).

Figure 2 .
Figure 2. Definition of effective depth, d e , from measured stage.S is the measured stage (L), S 0 is the minimum measures stage (L), and d min is the hydraulic mean depth of S 0 (L).

Figure 3 .
Figure 3.Normalized rating curves at 4 USGS gauge stations to demonstrate variations in m1 and m2 behavior centered around Q10 transition point (0,0).Note Q10 and d10 are 10th-percentile exceedance discharge (L 3 /T) and depth (L), respectively.m1 is the exponent for in-bank discharges, and m2 is the exponent for out-of-bank discharges.

Figure 3 .
Figure 3. Normalized rating curves at 4 USGS gauge stations to demonstrate variations in m 1 and m 2 behavior centered around Q 10 transition point (0,0).Note Q 10 and d 10 are 10th-percentile exceedance discharge (L 3 /T) and depth (L), respectively.m 1 is the exponent for in-bank discharges, and m 2 is the exponent for out-of-bank discharges.

Figure 4 .
Figure 4. Testing the normality of pre-logarithmically normalized (a) Q1, (b) Q10, (c) Q50, and where Q10 was previously defined, Q1 and Q50 represent the first and fiftieth percentile exc discharge (L 3 /T), and A10 is the cross-sectional area of Q10 (L 2 ).The red line represents th quantiles from a normal distribution, and the blue dots represent the actual quantiles.

Figure 4 .
Figure 4. Testing the normality of pre-logarithmically normalized (a) Q 1 , (b) Q 10 , (c) Q 50 , and (d) A 10 , where Q 10 was previously defined, Q 1 and Q 50 represent the first and fiftieth percentile exceedance discharge (L 3 /T), and A 10 is the cross-sectional area of Q 10 (L 2 ).The red line represents theoretical quantiles from a normal distribution, and the blue dots represent the actual quantiles.

Figure 5 .
Figure 5. Testing the normality of post-logarithmically normalized (a) log(Q 1 ), (b) log(Q 10 ), (c) log(Q 50 ), and (d) log(A 10 ).The red line represents theoretical quantiles from a normal distribution, and the blue dots represent the actual quantiles.

Figure 7 .
Figure 7. Multiple linear regression model for A 10 (m 2 ) using A D (km 2 ) and channel slope (S c ) (m/m).

Figure 10 .
Figure 10.Histogram to compare the RMSE between the improved model and Mueses e Model.The gray area represents the overlap range between the improved and Mueses mod

Figure 10 .
Figure 10.Histogram to compare the RMSE between the improved model and Mueses et al. [6] Model.The gray area represents the overlap range between the improved and Mueses models.

Figure 11 .
Figure 11.Evaluation of discharge rating curves using improved method (d ) (m) with A D , S c , and channel width (W C ); Mueses et al. [6] model (d Mueses ) utilizing only A D ; and observed USGS discharge rating curve (d) at three gauge stations: (a) station No. 02296750, (b) station No. 02295637, (c) station No. 02294161.

Figure 12 .
Figure 12.The relationship between AD and Q90.

Figure 12 .
Figure 12.The relationship between A D and Q 90 .

Table 1 .
Number of USGS stations and corresponding statistics for each variable (with units beside each variable).

Table 2 .
The Root Mean Squared Error (RMSE) and Mean Error (ME) of the developed regression models (with units beside each variable).

Table 2 .
The Root Mean Squared Error (RMSE) and Mean Error (ME) of the developed regression models (with units beside each variable).Notes: Q 10 , Q 50 , Q 1 , A 10 , d 10 , d 50 , and d 1 were previously defined.