Decision-Support System for LID Footprint Planning and Urban Runo ﬀ Mitigation in the Lower Rio Grande Valley of South Texas

: To address regional ﬂooding in the United States, federal and state agencies are adopting strict drainage policies in any large-scale commercial development within the watershed boundary. The conventional approach of implementing a wet detention pond (WP) reduces the land cover and causes operation and maintenance challenges eventually. The present study developed a decision-support system (DSS) in the Lower Rio Grande Valley region of South Texas for optimal selection of Best Management Practices (BMPs) by substituting a portion of the WP footprint with three regionally promising low-impact development practices, namely, porous concrete pavement (PCP), bioretention (BR), and bioswale (BS). Source Load Assessment and Management Model for Windows (WinSLAMM) was used as the foundation for the DSS database and algorithm development. This tool suggested that the implementation of bioswale alone can considerably reduce the footprint and construction cost. Less than 0.95 ha of installation of BR and BS can mitigate 79–91% of runo ﬀ from a maximum of 5 ha of commercial development. A combination of BR, BS, and WP was found to reduce runo ﬀ signiﬁcantly (~100%), which suggests that the successful adoption of DSS might support better planning of the urban stormwater management in the Lower Rio Grande Valley (LRGV). evaporation rate was considered to be 0.25 mm / h. A 5% impervious cover was assumed for the predevelopment condition of the site. The DSS was simulated to determine the footprint area for each control practice required to meet the 10-year storm discharge goal as a function of the size of the development (0.1–6 ha). The shaded area represents the recommended footprint of control practices estimated by the tool. The DSS also estimated the outfall peak discharge values for controlled conditions to check the technical feasibility of the selected BMP option, which is shown by the connected rectangular dots. The dashed dots represent the targeted discharge from the 10-year 24-h storm event. As per the rainfall-runo ﬀ algorithm written for the rational method, the DSS-simulated 10-year peak discharge followed a linear trend concerning the size of the development. The DSS reveals that the commercial development of 2 ha is expected to release a runo ﬀ at a predeveloped rate of fewer than 0.06 m 3 / s to meet the 10-year storm discharge. Since no parameters were altered, these results were the same for all control practice scenarios (WP, PP, BR, and BS). The ANOVA results suggest that there is a signiﬁcant di ﬀ erence ( p > 0.05) in the footprint recommendation for di ﬀ erent BMP options for the same size of commercial development. The values from Pearson multiple correlation (R) analysis (Table S4) showed that the runo ﬀ reduction and runo ﬀ coe ﬃ cient for each BMP have a strong inverse correlation (R = − 0.99) with the peak discharge.


Introduction
Rapid urbanization and industrialization have significantly fueled the total percentage of impervious land cover (such as traditional parking lots driveways, sidewalks, roofs, etc.) within the major watersheds of the United States (US) [1][2][3][4].Urban land developments increased by four times between the years 1945 and 2010 [5], which eventually altered the watershed hydrology and triggered an uncontrolled runoff towards major rivers in the US [6].The first flush of urban runoff carries non-point source (NPS) pollutant load, which eventually adversely affects the water quality of the watershed [7][8][9].Since the 1980s, federal and state agencies have been working to update the watershed protection plans and adjust existing drainage policy measures to address the impairment of the US rivers, with the involvement of cities and the Stormwater Task Force (STF) partners [10,11].
In most US states, the conventional approach to meet the stormwater discharge goal through the detention of runoff is to design a wet or dry detention pond with a large footprint and maintenance burdens [12].The flat terrain of the watershed provides stormwater engineers with complicated flow, detention, and flood design problems [13].Although designing detention ponds is somewhat effective and the calculation of storage is simple using Excel spreadsheets, this approach may cause the reduction of land cover; problems with safety, operation, and maintenance; and aesthetic issues in the future [14].Communities have been stressing successful planning of innovative stormwater management techniques to reduce non-point source pollutants; to mitigate localized flooding in urban, colonial, and rural settings; and to improve water quality [15].
The implementation of Low-Impact Development (LID)-Best Management Practices (BMPs)-has been encouraged to promote onsite mitigation of polluted runoff before discharging into the local storm sewer [10].With the use of innovative LID strategies, the traditional design of a large detention pond footprint can be reduced or eliminated by allowing the LID footprint to be incorporated into design detention calculations [16].LID is an ecologically based engineered landscape that helps to return the urbanized watersheds to predevelopment hydrological conditions.LID differs from the conventional approach because it has the function of returning the runoff to the hydrological cycle, including runoff volume reduction and improvement in water quality [17].It also provides environmental services, territorial cohesion, preservation of biodiversity, and sustainable development [18,19].However, LID is not typically integrated at the beginning of the planning process because of the high capital and retrofitting cost for some BMPs [12,20], more specifically in rapidly grown urban areas.Even though it is recommended on some occasions, it is not mandated in the drainage ordinances [21].The recent challenge is to decide how to plan LIDs efficiently within a commercial development to meet a certain storm discharge goal from the drainage site.
Several decision-support systems were developed as tools to compare and evaluate different LID scenarios within the watershed boundary [22].One study developed a tool for the placement of LIDs based on the defined precipitation and soil moisture conditions.Results showed the areas that are the most suitable locations for LID installments and the optimal facility with higher retention abilities of the accumulated runoff [23].Another study developed an algorithm to allow the selection of LIDs with the lowest cost, the highest runoff reduction, and the highest likelihood of private-owner maintenance in various sub-watersheds [24].Another similar study developed a decision-making tool by evaluating various alternatives to determine the most cost-effective types and combinations of LIDs that maximize the water quality benefits [25].In dealing with the complex issues in stormwater management, the use of decision-support systems supported by modeling tools is becoming increasingly popular [26].Analyzing LID through coupling modeling software with a multi-objective optimization using genetic algorithms has been a common method [22,27].Ercolani et al. used the MOBIDIC-U hydrologic model equipped with SMART-GREEN, which is a Quantum Geographic Information System (QGIS) plugin interface, to study the effects of green roofs in an urban watershed [28].Their study revealed that green roofs might work better under more frequent storms of smaller magnitude, and their efficiency can be increased with spatially heterogeneous implementation.Martin-Mikle et al. presented a spatially explicit method using Geographic Information System (GIS) for the placement of LIDs within the urban watershed to increase cost-effectiveness and ecological benefits [5].Meeorw and Newell developed a Green Infrastructure Spatial Planning (GISP) model to aid in identifying trade-offs, synergies, and hotspots for future green infrastructures (GI) with maximum ecosystem services [29].Charoenkit and Piyathamrongchai developed a framework for the comparative assessment of multifunctionality among different urban ecosystem structures [30].A previous study evaluated a machine learning approach using artificial neural network (ANN), adaptive network-based fuzzy inference system (ANFIS) algorithms, and statistical modeling to predict the probability of green or grey transformation for vacant places, poor condition sites, and waterway corridors [31].Previously, the Storm Water Management Model (SWMM)-calibrated green infrastructure models were used to assess runoff reductions via different location strategies within the watershed, both spatially known and unknown [32].
Previous studies evaluated the performance of LID systems to be location-dependent.Several factors can affect the LID effectiveness; efficacy depends on the type, design, and local conditions such as the topography of the site, the soil type/conditions, the rainfall patterns, and other types of hydrological and meteorological properties [33].Due to the variation of each watershed's character, there are many disciplinary factors, such as hydrology (stormwater runoff reduction), ecology, health (air and water pollution), urban landscape (social activity), etc.The guidelines for BMP planning are not consistent, even between neighboring cities, due to the hydrogeological variation.For example, stormwater runoff generated from a 25-year frequency storm event in a new commercial development in San Benito, Texas was generally required to be detained on-site to be released into a receiving system at a 25-year predevelopment discharge rate.On the other hand, the city of Weslaco followed more stringent drainage policies, with a 10-year predevelopment discharge rate because of its lower elevation [34].Most of the previous studies focused on the development of a framework for identifying the optimal distribution of BMPs within a mixed-use watershed for maximizing ecosystem services [5,32].One previous study also prioritized economic evaluation as a decision-support tool for strategic spatial planning of green infrastructures [35].Very few studies have attempted to develop a region-specific decision-support framework to assess the LID transformation in a small-scale or commercial development boundary to reveal both maximum ecosystem services and minimum costs.Still, there is a lack of a well-established tool to enforce the local drainage policy in optimal selection and evaluation of combined BMP scenarios for maximizing hydrologic benefits in urban development [32,36,37].The proper selection of a BMP or BMPs in series through innovative planning, management, and engineering approaches will not only benefit communities but will also protect the water quality and ecosystems of the watershed [7].
Source Load Assessment and Management Model for Windows (WinSLAMM) is a Stormwater Quality Modeling Tool, which can evaluate stormwater controls based on the field environment and actual design values [38].It was designed to provide relatively complex outputs, such as runoff volume, pollutant mass loadings, and the effects of critical conditions on control practices.It is mostly used as a planning tool to better understand the sources of urban runoff pollutants and their control strategies.It has been widely used and extensively reviewed in many LID implementation studies in North America and has been proven reasonable in predicting stormwater flows and pollutant characteristics from a wide range of rainfall events, urban development characteristics, and control practices [38][39][40].WinSLAMM predicts runoff based on the results of field data from the National Stormwater Database and actual design values of BMPs [38,41].Therefore, the model is capable of predicting results even when there are no field-input data by using the complete package of default-calibrated files generated from previous LID data collection within the USA.It uses a simple rainfall-runoff equation (Rational Method) to predict runoff volumes, with errors ranging from 10% to 30% [40].By using 3-5 years of rainfall-runoff data, it is possible to calibrate the model for future runoff and even more reliable water quality prediction.WinSLAMM is highly advanced in its LID control features [42].The model can use BMPs as external control practices beyond the boundary of land use development.The model outputs can be summarized per month and incorporated into the SWAT model and Geographic Information System (GIS) Platforms [39,[43][44][45].
To address the water quality issues at the US-Mexico border associated with non-point source pollution in the Lower Rio Grande Valley (LRGV) region within the Arroyo Colorado watershed, this study was designed to develop an innovative Decision-Support System (DSS) for the planning of a BMP(s) footprint to replace a portion of conventional detention practices with LID echo-technologies.The preliminary objective of this study is to generate calibrated-BMP models of three LID BMPs (porous concrete pavement, bioretention, and bioswale) using WinSLAMM for the semi-arid climatic region of LRGV in South Texas.The primary objective is to apply BMP models to establish a novel methodology for the DSS development and to assess its feasibility in determining the footprint of an individual BMP or BMPs in series, as well as their hydrologic performance (peak flow, runoff volume reduction), for a range of sizes of commercial developments.In a broader picture, this novel algorithm is expected to evaluate existing BMP footprint requirements or seek LID options as an alternative to existing or conventional facilities (e.g., wet detention pond).If the tool proves adequate, a similar algorithm can be applied to other impaired watersheds in the US and the rest of the world.

Study Case
This study focused on the Arroyo Colorado River, which passes through the heart of the Lower Rio Grande Valley (LRGV) of South Texas, which extends up to 145 km from Mission, TX, to the Lower Laguna Madre.The total reach of the river has been divided into two segments by the Texas Commission on Environmental Quality (TCEQ): the tidal segment (Segment 2201) and the above-tidal segment (Segment 2202) [4,13,46].TCEQ identified both of these segments as impaired due to high nutrient and bacterial concentration and low Dissolved Oxygen (DO) level.Rapid urbanization and industrialization adversely replaced the existing vegetative cover [13,46,47], which eventually triggered the non-point source (NPS) pollution within the Arroyo Colorado watershed (Figure 1).The impairment inspired the development of the Arroyo Colorado Watershed Protection Plan in 2007 by local, state, and federal stakeholders to improve the surface water quality of the watershed.
Under the TCEQ-sponsored Clean Water Act 319 LRGV LID implementation program, flow data were monitored from six LID BMPs initially.Flow data were monitored from BMPs installed in different locations in the LRGV from August 2014 through March 2017 [48].From our field assessment, three LID BMPs (porous pavement, bioretention, and bioswale) were found "most promising" in terms of peak flow and runoff volume reduction in the semi-arid climatic region of LRGV in South Texas.Therefore, the following BMPs were prioritized in this study to be incorporated into the DSS algorithm:

Study Case
This study focused on the Arroyo Colorado River, which passes through the heart of the Lower Rio Grande Valley (LRGV) of South Texas, which extends up to 145 km from Mission, TX, to the Lower Laguna Madre.The total reach of the river has been divided into two segments by the Texas Commission on Environmental Quality (TCEQ): the tidal segment (Segment 2201) and the abovetidal segment (Segment 2202) [4,13,46].TCEQ identified both of these segments as impaired due to high nutrient and bacterial concentration and low Dissolved Oxygen (DO) level.Rapid urbanization and industrialization adversely replaced the existing vegetative cover [13,46,47], which eventually triggered the non-point source (NPS) pollution within the Arroyo Colorado watershed (Figure 1).The impairment inspired the development of the Arroyo Colorado Watershed Protection Plan in 2007 by local, state, and federal stakeholders to improve the surface water quality of the watershed.
Under the TCEQ-sponsored Clean Water Act 319 LRGV LID implementation program, flow data were monitored from six LID BMPs initially.Flow data were monitored from BMPs installed in different locations in the LRGV from August 2014 through March 2017 [48].From our field assessment, three LID BMPs (porous pavement, bioretention, and bioswale) were found "most promising" in terms of peak flow and runoff volume reduction in the semi-arid climatic region of LRGV in South Texas.Therefore, the following BMPs were prioritized in this study to be incorporated into the DSS algorithm:

DSS Framework for LID Footprint Planning and Evaluation
The DSS was developed as an innovative calculator for the LID footprint planning in a proposed commercial development to meet the 10-year discharge goal from a 50-year rainfall event.As a planning tool, the DSS algorithm has been written in the Visual Basic programming language and formulated as a macro-enabled Excel spreadsheet to support the discharge calculation as a basis for BMP planning.The user has control over the selection of a city, which is linked to the local soil characterization parameters (runoff coefficient, surface infiltration rate, etc.).Feedback from local engineers/planners/stakeholders was gathered to assess city-specific hydrologic and geologic information.The tool was successfully designed to accept input parameters, such as development size, predeveloped land cover, evaporation rate, type of BMP(s), and per unit construction cost.The

DSS Framework for LID Footprint Planning and Evaluation
The DSS was developed as an innovative calculator for the LID footprint planning in a proposed commercial development to meet the 10-year discharge goal from a 50-year rainfall event.As a planning tool, the DSS algorithm has been written in the Visual Basic programming language and formulated as a macro-enabled Excel spreadsheet to support the discharge calculation as a basis for BMP planning.The user has control over the selection of a city, which is linked to the local soil characterization parameters (runoff coefficient, surface infiltration rate, etc.).Feedback from local engineers/planners/stakeholders was gathered to assess city-specific hydrologic and geologic information.The tool was successfully designed to accept input parameters, such as development size, predeveloped land cover, evaporation rate, type of BMP(s), and per unit construction cost.The database and algorithm of the DSS were generated from the output resulting from the assessment, evaluation, and field-calibrated WinSLAMM model of three BMPs (permeable pavements, bioretention, and bioswale).The tool was designed to provide the outcomes, such as required BMP footprint, detention volume, size, and the number of BMP units needed to accomplish the objectives.The tool can predict the uncontrolled and controlled performance of the site in terms of peak discharge, total runoff volume, percentage of runoff reduction, runoff coefficient, etc.The tool is also able to estimate the total lifecycle cost of BMP(s).Figure 2 shows the overall framework of the DSS for the accomplishment of the objective of this study.database and algorithm of the DSS were generated from the output resulting from the assessment, evaluation, and field-calibrated WinSLAMM model of three BMPs (permeable pavements, bioretention, and bioswale).The tool was designed to provide the outcomes, such as required BMP footprint, detention volume, size, and the number of BMP units needed to accomplish the objectives.
The tool can predict the uncontrolled and controlled performance of the site in terms of peak discharge, total runoff volume, percentage of runoff reduction, runoff coefficient, etc.The tool is also able to estimate the total lifecycle cost of BMP(s).Figure 2 shows the overall framework of the DSS for the accomplishment of the objective of this study.

BMP Data Acquisition and Analysis
The raw flow data were retrieved from the ISCO Signature flow meter and analyzed for each rainfall event to calculate the flow volume at the BMP outfall.The total inflow volume onto the BMPs was calculated by summing up the direct rainfall volume (calculated from the ISCO 674 Rain Gauge data) onto the BMP surface and the runoff generated in the drainage areas.Later, volume reduction was calculated and converted to a normalized unit by dividing it by the area of BMP installation.Normalized runoff volume reduction (ΔV) was calculated to evaluate the performance of each BMP by using the following equation (1): where Vi = total inflow volume (m 3 ), V0 = total outflow volume (m 3 ), and Ap = area of BMP installation (m 2 ).

WinSLAMM Model Development
A site-specific rainfall file (.RAN) was created for each type of monitored BMP.The particulate solids concentration (.pscx) and pollutant probability distribution (.ppdx) files were selected from the default standardized files for the Central US zone.The street delivery file (.std) was selected for the other urban and commercial land use developments.The particle size distribution and peak-toaverage ratio files (.csv) were selected from the NURP database.The "as-built" design information of each BMP was assimilated and entered as model input parameters.All cross-sectional design values were collected from the drawing sheets provided by the contractors.The subgrade seepage rate was

BMP Data Acquisition and Analysis
The raw flow data were retrieved from the ISCO Signature flow meter and analyzed for each rainfall event to calculate the flow volume at the BMP outfall.The total inflow volume onto the BMPs was calculated by summing up the direct rainfall volume (calculated from the ISCO 674 Rain Gauge data) onto the BMP surface and the runoff generated in the drainage areas.Later, volume reduction was calculated and converted to a normalized unit by dividing it by the area of BMP installation.Normalized runoff volume reduction (∆V) was calculated to evaluate the performance of each BMP by using the following Equation (1): where V i = total inflow volume (m 3 ), V 0 = total outflow volume (m 3 ), and A p = area of BMP installation (m 2 ).

WinSLAMM Model Development
A site-specific rainfall file (.RAN) was created for each type of monitored BMP.The particulate solids concentration (.pscx) and pollutant probability distribution (.ppdx) files were selected from the default standardized files for the Central US zone.The street delivery file (.std) was selected for the other urban and commercial land use developments.The particle size distribution and peak-to-average ratio files (.csv) were selected from the NURP database.The "as-built" design information of each BMP was assimilated and entered as model input parameters.All cross-sectional design values were collected from the drawing sheets provided by the contractors.The subgrade seepage rate was obtained from the geotechnical reports provided by the lead engineers.Within a designated land use, the type and area of source contributions to runoff generation were defined.For the first simulation, default runoff coefficient files (.rsvx) were selected for the Central US zone.The flow routing process was defined in the model window, which includes the delineation of land use, conduits, junctions, BMP, and outfall.Local calibration of BMPs was emphasized based on the degree of discrepancy between observed and simulated runoff reduction results.For significant deviation, the default runoff coefficient (R v ) values were adjusted to bring the simulated results close to the observed results.The goodness of fit between both results was statistically justified through determining the Kling-Gupta efficiency (KGE) and Root Mean Square Error (RMSE) [49].Furthermore, the p-value for supporting the correlation was reported through a regression hypothesis test.Later, all calibrated BMP models were validated with a separate observed dataset.The following Equations ( 2) and (3) were used to calculate statistical parameters to assess the correlation between observed and model-simulated results: where r = linear correlation between observations and simulations, σ obs = standard deviation in observations, σ sim = standard deviation in simulations, µ obs = observation mean, µ sim = simulation mean, n = number of observations, x i = model-simulated runoff reduction for a particular rainfall event (i) (mm), and y i = observed runoff reduction for a particular rainfall event (i) (mm).A higher value of KGE (>0.8), lower value of RMSE (<30%), and a p-value less than 0.05 were preferred to support the adequacy of the best-fitted correlation between observed and model-predicted results [50-53].

DSS Algorithm Development
The DSS uses the Rational Method (uncontrolled conditions), and WinSLAMM transformed BMP equations (controlled conditions) to calculate the peak discharge from proposed commercial developments for the Arroyo Colorado watershed.For the uncontrolled or baseline condition, the peak discharge from a rainfall event was calculated by taking the difference between the discharge calculated by the rational method and the volumetric evaporation rate (rate of evaporation, mm/h × total area of commercial development, m 2 ) [54].In the discharge calculation, the rainfall intensity was assumed to be the seasonal average for each of the three LRGV counties in Texas: Hidalgo, Willacy, and Cameron.From the Intensity-Duration-Frequency (IDF) curve for these counties, the rainfall intensity corresponding to the 50 and 10-year 24-h storm events was estimated to be 12.7 and 7.62 mm/h, respectively.The intensity is multiplied by a seasonal variation factor of 3.8 to simulate a wet weather condition.Within the algorithm, the total volume of runoff was calculated by multiplying the rainfall duration by the peak discharge as per flow-volume relationship [55].The following Equations ( 4)- (7) were written for the DSS algorithm to calculate 50-and 10-year peak discharge and runoff volume: V p,50 = Q p,50 × t (5) V p,10 = Q p,10 × t where Q p,10 , Q p,50 = peak discharge from 10-and 50-year 24-h storm events (m 3 /s); V p,10 , V p,50 = peak volume of runoff from 10 and 50-year 24-h storm events (m 3 ); i 10 , i 50 = average rainfall intensity of 10and 50-year 24-h storm events; A = total area of proposed commercial development (ha); A i = total area of existing impervious cover at the predeveloped site (ha); C i = runoff coefficient for existing impervious cover = 0.98 (assumed); A p = total area of impervious cover at the predeveloped site (ha); C p = runoff coefficient for impervious cover = 0.98 (ha); SVF = seasonal variation factor, for wet seasons = 3.8 (assumed); E = evaporation rate (mm/h); t = time.
In the development of the DSS algorithm for the controlled condition, WinSLAMM was used to create a database for four BMPs.The database was generated in terms of peak discharge as an output resulting from the WinSLAMM-calibrated BMP models for varying development sizes and installation areas.The information gathered from the database was transformed into equations to build the algorithm for the DSS in the Visual Basic programming language.The calculated volumetric evaporation rate was subtracted from the peak discharge calculated from the WinSLAMM-translated BMP equation to estimate the discharge after the installation of a BMP that covers a specific footprint area.The following Equation ( 8) is the generic form of the algorithm written to calculate the peak discharge with a single control practice: where i = BMP footprint trial number; Q p,CP (i) = peak discharge with single control practice for an ith footprint trial (m 3 /s); b i , d i = regression coefficients, which depend on the size of the BMP footprint; A = total area of proposed commercial development (ha); E = evaporation rate (mm/h).
For BMPs in series, the summation of peak discharge reductions for all control practices was subtracted from the 50-year 24-h peak discharge to calculate the final peak discharge.The following Equation ( 10) is the generic form of the algorithm: Equation ( 10) can be simplified and written in the following form: where Q p,comb, (i) = peak discharge when a commercial development is incorporated with the combination of ith footprint trial of all BMPs in series (m 3 /s); Q p,CP1 = peak discharge when a commercial development is incorporated with LID control practice 1; Q p,CP2 = peak discharge when a commercial development is incorporated with LID control practice 2; Q p,CP3 = peak discharge when a commercial development is incorporated with LID control practice 3; Q p,50 = peak discharge from the 50-year 24-h storm event (m 3 /s).
Later, a convergence algorithm in Visual Basic was developed for BMP size determination.If the calculated Q p, CP (i) or Q p, comb, (i) is found to be less than Q p,10 for a specific footprint of the BMP(s), the algorithm ceases and considers the footprint size to be finalized.Otherwise, the algorithm goes for a further trial to look for the next available size and check for convergence again.

WinSLAMM-Calibrated BMP Model Development
Figure 3 shows the WinSLAMM calibration plots for three different types of BMPs: permeable pavements, bioretention, and bioswale.The plots show the runoff reduction behavior of these BMPs as a function of a range of rainfall magnitudes.The scattered markers represent the observed runoff reduction for a certain rainfall depth.However, only rainfall events greater than 1.0 mm were found to be sufficient to trigger a significant volume of runoff.For calibration purposes, the default runoff coefficient (R v ) values were adjusted to fit the model-simulated results to be closer to the observed results when needed.The model-simulated results were fitted with a second-order polynomial trendline since the observed reduction pattern was not found to be linear to rainfall depth.Previous studies observed that the outflow volume followed a curvilinear trend as rainfall magnitude increased.The goodness of fit and the statistical correlation between observed and simulated results were determined through the Kling-Gupta Efficiency (KGE), Root Mean Square Error (RMSE), and regression hypothesis test (p-value).Table 1 summarizes the monitoring period, the number of rainfall events, max rainfall depth, calibrated parameters, and performance indexes of the period for the three BMPs.rainfall events, max rainfall depth, calibrated parameters, and performance indexes of the period for the three BMPs.A total of 14 significant rainfall events were monitored from (September 2014 to November 2014) for the modeling of the porous concrete pavement (PCP), which is located at Monte Bella Park, Brownsville, TX (Figure 3A).The highest depth of rainfall was observed to be 67 mm within the dataset.The runoff reduction behavior seemed to follow an increasing trend against rainfall depths.Within this monitoring timeframe, the normalized volume reduction achieved was 0.054 ± 0.047 m 3 /m 2 .The model-simulated results showed a reasonable fit (KGE = 0.833, RMSE = 0.012 m 3 /m 2 , p < 0.05) with the observed data with default R v values.Therefore, the PCP model did not undergo further calibration trials.For the modeling of bioretention, 43 rainfall events (March 2016 to March 2017) were monitored from the South Texas college parking lot (Figure 3B).The highest depth of rainfall was observed to be 38 mm inches within the dataset.Observed reductions followed a steadily increasing trend, which was a good indication of significant volume reduction performance (0.25 ± 0.21 m 3 /m 2 ) from the bioretention.The model was fitted (KGE = 0.866, RMSE = 0.054 m 3 /m 2 , p < 0.05) to the observed results by decreasing the default R v values by 66% for pervious surfaces.Eleven rainfall events (August 2014 through February 2015) were monitored for the modeling of bioswale installed in the Cascade Park parking lot (Figure 3C).The highest depth of rainfall was 67 mm within the observed dataset.The model-simulated results with default R v values were a good fit (KGE = 0.893, RMSE = 0.033 m 3 /m 2 , p < 0.05) with the observed data.Thus, no further calibration was necessary for the bioswale model.Within the monitoring timeframe, the normalized reduction was observed at 0.13 ± 0.14 m 3 /m 2 .Overall, bioretention (BR) and bioswale (BS) showed somewhat better performance than porous concrete pavement (PCP), perhaps because of its improved infiltration rate through the vegetative cover.Our validation results (Tables S1-S3) showed that the WinSLAMM-calibrated PCP, BR, and BS models might over or under-predict within a 30% error range.Later, these calibrated models were used to calculate the controlled peak discharge for different sizes of commercial developments and BMP installations to create a BMP-specific database and equations for the integration into the DSS algorithm.

Application of DSS
The application of the DSS was evaluated to identify a suitable control practice(s) with a recommendation on the footprint area and number of units to meet the 10-year storm discharge at the outfall.According to the user-provided site and BMP information, our DSS predicted the size and units of the structure required to meet the 10-year discharge goal.During the DSS simulation, the algorithm utilized BMP-specific calibrated equations as the basis for the estimation of hydrologic outputs (peak discharge, runoff volume, percentage runoff reduction, and runoff coefficient).The area of the BMP footprint was determined after several iteration trials (convergence algorithm) until the calculated peak discharge met the targeted 10-year peak discharge.Apart from the footprint, the tool determined the total detention volume, geometric recommendations, and the number of BMP units required, none of which were the focus of our analysis.The tool also predicted the total lifecycle cost of LID implementation.
An illustration of the simulated results of the DSS was created to assess the effects of changing the development size, ranging between 0.1 and 6.0 ha, in Brownsville, TX, on the LID-controlled peak discharge from the site, as shown in Figure 4.In the preparation of the input data, the typical evaporation rate was considered to be 0.25 mm/h.A 5% impervious cover was assumed for the predevelopment condition of the site.The DSS was simulated to determine the footprint area for each control practice required to meet the 10-year storm discharge goal as a function of the size of the development (0.1-6 ha).The shaded area represents the recommended footprint of control practices estimated by the tool.The DSS also estimated the outfall peak discharge values for controlled conditions to check the technical feasibility of the selected BMP option, which is shown by the connected rectangular dots.The dashed dots represent the targeted discharge from the 10-year 24-h storm event.
As per the rainfall-runoff algorithm written for the rational method, the DSS-simulated 10-year peak discharge followed a linear trend concerning the size of the development.The DSS reveals that the commercial development of 2 ha is expected to release a runoff at a predeveloped rate of fewer than 0.06 m 3 /s to meet the 10-year storm discharge.Since no parameters were altered, these results were the same for all control practice scenarios (WP, PP, BR, and BS).The ANOVA results suggest that there is a significant difference (p > 0.05) in the footprint recommendation for different BMP options for the same size of commercial development.The values from Pearson multiple correlation (R) analysis (Table S4) showed that the runoff reduction and runoff coefficient for each BMP have a strong inverse correlation (R = −0.99)with the peak discharge.An illustration of the simulated results of the DSS was created to assess the effects of changing the development size, ranging between 0.1 and 6.0 ha, in Brownsville, TX, on the LID-controlled peak discharge from the site, as shown in Figure 4.In the preparation of the input data, the typical evaporation rate was considered to be 0.25 mm/h.A 5% impervious cover was assumed for the predevelopment condition of the site.The DSS was simulated to determine the footprint area for each control practice required to meet the 10-year storm discharge goal as a function of the size of the development (0.1-6 ha).The shaded area represents the recommended footprint of control practices estimated by the tool.The DSS also estimated the outfall peak discharge values for controlled conditions to check the technical feasibility of the selected BMP option, which is shown by the connected rectangular dots.The dashed dots represent the targeted discharge from the 10-year 24-h storm event.As per the rainfall-runoff algorithm written for the rational method, the DSS-simulated 10-year peak discharge followed a linear trend concerning the size of the development.The DSS reveals that the commercial development of 2 ha is expected to release a runoff at a predeveloped rate of fewer than 0.06 m 3 /s to meet the 10-year storm discharge.Since no parameters were altered, these results were the same for all control practice scenarios (WP, PP, BR, and BS).The ANOVA results suggest that there is a significant difference (p > 0.05) in the footprint recommendation for different BMP options for the same size of commercial development.The values from Pearson multiple correlation (R) analysis (Table S4) showed that the runoff reduction and runoff coefficient for each BMP have a strong inverse correlation (R = -0.99)with the peak discharge.Figure 4A reveals the size of the wet detention pond (WP), which might be adequate in runoff management for a range of development sizes (0.1-6.0 ha).Results suggested that the 2 ha development size might require the building of a 0.8 ha wet pond that would induce the least peak discharge of 0.034 m 3 /s and a maximum runoff reduction of 94% from the site.The installation of 0.8 ha of WP might cover the development size ranging between 0.1 and 6.0 ha.The DSS suggests that 1.22 ha of WP alone can be promising in terms of successful mitigation of runoff (94%) from 4 ha of development.The DSS also did the consultation to make a recommendation on LID BMPs. Figure 4B suggests that it might take 1.25 ha of porous concrete pavement (PCP) installation from 2.0 ha of development.The cross-sectional profile of the modeled PCP consists of a clean, washed, and angular aggregate, typically #57 (as a base) and #2 (as sub-base).The interim layers of PCP form 35% of the void space between the aggregates, which may only be counted as detention.Therefore, the DSS estimated a higher footprint for PCP as compared to WP.The DSS results suggest that 1.25 ha of PCP installation should be enough to handle runoff from a maximum of 2.4 ha of development.After that, the PCP size might differ along with the size of the development.Installation of 4 ha of PCP might serve a development size as large as 6.1 ha.These results reveal the idea that PCP might be capable to detain 34-50% of runoff generated outside its footprint boundary.The DSS also estimated the footprint of bioretention control structures to be 0.45 ha for the same 2 hectares of development, as shown in Figure 4C.The more preferable size would be 0.92 ha, which might cover a development size as large as 5 ha with 90-95% runoff reduction.For the same development size, the DSS recommends a comparatively smaller bioswale (BS) footprint (0.35 ha) than bioretention, with a similar peak discharge (0.1 m 3 /s) since they are more linear systems with a greater length than width.Only 0.7 ha of BS installation might be adequate to withstand the runoff received from the maximum 4.5 ha of commercial development (84% bigger than the BS footprint itself).The results summary for individual BMP cases is given in Table 2.The DSS module also provides the user with an option to choose BMPs in series.The simulation through the DSS recommends a footprint for all probable BMP configurations when they are installed in series.The comparison of different BMP combinations in terms of footprint requirements for 4 ha of development is shown in Figure 5. DSS inputs were fixed to a 4 ha commercial development in the city of Brownsville, TX.The figure distinguishes the distribution of footprint area between BMPs in the form of color-coded bar graphs.During the simulation, the 50-year uncontrolled peak discharge was used as an input to estimate the total inflow volume and runoff reduction (by subtracting the LID-controlled runoff volume).The figure also depicts the performance in terms of percentage runoff reduction after the implementation of the BMPs in the series.
When LIDs are installed in series with a wet detention pond, they may raise the total land use to be higher than that expected from a WP alone, as suggested by the DSS.The maximum 30% of the wet pond footprint can be replaced when combined with other BMPs in series.The DSS estimates a maximum BMP installation area of 2.5 ha when the exfiltration from the PCP is designed to be detained in the hydrologically connected downstream WP.It is important to note that the involvement of PCP in series can cover a major portion (>50%) of the total areas of BMP installation.The implementation of bioswale (0.72 ha, $1.25 M) or bioretention (0.92 ha, $3.5 M) alone can considerably reduce the construction cost (>75%) and promote the aesthetics of the site.A combined total of 0.47 ha of bioretention (incorporated as an island or median), 0.36 ha of bioswale (on the periphery of the site), and 0.85 ha of large wet detention pond (at the downstream) is the most promising configuration, with almost 100% runoff reduction from the site.Such a result indicates that the post-development condition that contributes to flooding can be reduced to zero under this configuration.When LIDs are installed in series with a wet detention pond, they may raise the total land use to be higher than that expected from a WP alone, as suggested by the DSS.The maximum 30% of the wet pond footprint can be replaced when combined with other BMPs in series.The DSS estimates a maximum BMP installation area of 2.5 ha when the exfiltration from the PCP is designed to be detained in the hydrologically connected downstream WP.It is important to note that the involvement of PCP in series can cover a major portion (>50%) of the total areas of BMP installation.The implementation of bioswale (0.72 ha, $1.25M) or bioretention (0.92 ha, $3.5M) alone can considerably reduce the construction cost (>75%) and promote the aesthetics of the site.A combined total of 0.47 ha of bioretention (incorporated as an island or median), 0.36 ha of bioswale (on the periphery of the site), and 0.85 ha of large wet detention pond (at the downstream) is the most promising configuration, with almost 100% runoff reduction from the site.Such a result indicates that the post-development condition that contributes to flooding can be reduced to zero under this configuration.

Discussion
Our LID decision-support algorithm predicted that a single wet detention pond (WP) of 1.22 ha might efficiently handle a 50-year frequency storm event from 4 ha of commercial development (peak inflow ~0.53 m 3 /s) in Brownsville, Texas, which can release the least peak outflow of 0.034 m 3 /s with a maximum flow reduction of 94%.Although WP can be promising in terms of runoff detention (94%), its implementation can sometimes produce both unpleasant and unwanted liabilities.The main concern for a developer is the high amount of land utilization and the maintenance and

Discussion
Our LID decision-support algorithm predicted that a single wet detention pond (WP) of 1.22 ha might efficiently handle a 50-year frequency storm event from 4 ha of commercial development (peak inflow ~0.53 m 3 /s) in Brownsville, Texas, which can release the least peak outflow of 0.034 m 3 /s with a maximum flow reduction of 94%.Although WP can be promising in terms of runoff detention (94%), its implementation can sometimes produce both unpleasant and unwanted liabilities.The main concern for a developer is the high amount of land utilization and the maintenance and construction costs associated with WP [12,14].Due to improper maintenance, the wetland can cause health and safety issues, declination of waterfowl population, clogging, mosquito breeding, etc. [56].The DSS estimated the total cost (including construction, operation & maintenance (O&M), and others) of 1.22 ha of WP in a 4 ha development to be almost 16 million, considering $236 per cubic meter of construction cost.
The DSS results suggest that the wet pond requirements can easily be avoided by using porous concrete pavement (PCP), bioretention (BR), or bioswale (BS) alone.Although PCP might require more space, the highly durable concrete or asphalt surfaces pose great serviceability by providing additional benefits when utilized in parking lots and areas with light traffic (driveway, shoulder, etc.) [57].Like other conventional pavements, it can also support high axle loads from fire trucks, garbage trucks, and emergency vehicles without any problem [58].Our present study in the semi-arid climatic LRGV observed a comparatively higher normalized runoff volume reduction from PCP than in the eastern sub-tropical, humid part of the US (e.g., North Carolina), where the reduction was calculated to be as low as 0.008 ± 0.006 m 3 /m 2 [59].This comparison indicates that our established DSS algorithm might predict higher space requirements for the same BMPs in sub-tropical climates.Our DSS estimated the total cost of 2.03 ha of PCP in a 4 ha development to be almost 3.5 million, considering $638 per cubic meter of construction cost.The DSS predicted the most promising results from the BS in terms of land use utilization to meet 10-year frequency targeted discharge (0.12 m 3 /s), with the peak outflow rate of 0.113 m 3 /s and a smaller footprint (0.72 ha) than other types.The BS structure is often designed for better conveyance, collection, and infiltration of stormwater runoff [60].The structure is highly vegetated with grass throughout its length and sloped sides to improve the surface runoff and infiltration, which was also encountered in our modeled bioswale.Before infiltration, the structure also causes the runoff water to pool on the surface temporarily [61,62].On the other hand, BR alone might induce a similar outlet peak discharge (0.046 m 3 /s) and runoff volume reduction (90-95%) to a wet pond, with much less land utilization (0.93 ha) than porous concrete pavement and wet pond in a 4-5 ha commercial development in Brownsville, TX.Bioretention can be described as a porous landscape detention basin since it captures a huge volume of runoff through its highly porous media [63].The top layer of our modeled BR cell was vegetated, deep-rooted native plants and local soil, which together build a strong foundation to enhance the on-site infiltration and evapotranspiration of the stormwater runoff, with less chance of potential clogging [64].Our modeled BR structure also provides an essential ponding depth (up to 0.3 m), which might temporarily hold the runoff from the parking lots and roadway [65].Hence, our field data also support the fact that the normalized volume reduction from the bioretention (0.25 ± 0.21 m 3 /m 2 ) was higher than that of PCP (0.054 ± 0.047 m 3 /m 2 ).The performance of our monitored bioretention in the semi-arid LRGV was even higher than that for the sub-tropical state of Virginia, US, where normalized volume reduction was calculated as 0.11 ± 0.23 m 3 /m 2 [66].A study in Italy showed that a single green roof with 100% green conversion could achieve a runoff volume reduction of up to 35% [67].Our DSS predicted that WP, PP, BR, or BS alone might achieve more than 75% runoff reduction, with a green conversion of 31%, 51%, 23%, and 18%, respectively.These results indicate a greater hydrologic potential of these four BMPs than the green roof.Our DSS estimated $3.6 and 1.6 million as the total life cycle cost for BR and BS implementation, respectively, in a 4 ha commercial development in Brownsville, TX, considering $546 per cubic meter of construction cost.Considering the DSS-predicted hydrological benefits and cost, BR alone might be the optimal option for any commercial development in the LRGV region.However, Zhang and Hu [68] observed the most promising runoff volume reduction (up to 100%) from a rainwater harvesting cistern in China in the case of a critical rainfall storm (50 mm).
In terms of volume reduction, a combination of PCP, BR, BS, and WP in series might require the least construction cost of $7.5 million, with a runoff reduction of 78%.By incorporating this combination, the serviceability and property value of the site can be increased.The DSS also predicted a combination of BR (0.47 ha), BS (0.36 ha), and WP (0.84 ha) in series, which might be the most promising option to handle 50-year frequency storm event runoff reduction of up to 100% from a 4 ha commercial development with a total cost of $9 million approximately.Besides, this combination can improve land use management and foster healthy aesthetics on the site.Considering all hydrological benefits and costs, this configuration might be optimal for combined BMP cases for any commercial development in the LRGV region.
With more flexible drainage policy options, the recommended size of BMPs can be reduced.However, the adopted stringent drainage policy in the DSS algorithm provided an opportunity for more sustainable runoff management.Also, it is possible to improve DSS recommendations through the integration of innovative BMP design models with higher storage capacity into the DSS algorithm.
Overall, the results of the presented analysis may support planners, stakeholders, and stormwater task force partners for optimized planning, design, and implementation of LID facilities in any commercial development in the region.Considering DSS-predicted hydrologic results (peak flow, runoff volume reduction) and cost, the most promising BMP or BMPs in series can be suggested at the beginning of the land development process in the semi-arid climatic region of LRGV.These results can also help to evaluate any existing/conventional facilities (e.g., wet detention pond) or to seek LID alternatives to maximize ecosystem services with a minimized cost.Furthermore, the outcomes of the tool can be incorporated into the enhancement of property development options and values in the region.Sustainable urban stormwater management will eventually lead to surface water quality improvement in the Arroyo Colorado watershed.

Conclusions
The present study established a novel methodology of a decision-support system for the planning of LID footprint using calibrated BMP models with the estimation of hydrologic performances (peak flow, runoff volume reduction) in any commercial development in the Lower Rio Grande Valley of South Texas.The database and algorithm of the DSS were developed based on WinSLAMM simulations of three BMPs (permeable pavements, bioretention, and bioswale) of different installation sizes, within a range of sizes of commercial development.Later DSS results were evaluated with specific site information to identify the suitable control practice(s) along with the footprint recommendation to meet the 10-year storm discharge from the site.The following are highlights of some important findings gathered from the DSS recommendations in the present study:

•
A wet detention pond (WP) alone can be promising to provide hydrologic benefits, but it is not optimal.WP can achieve the least peak discharge (0.034 m 3 /s) and maximum runoff reduction (94%) from a 4 ha development.However, WP can cause safety, aesthetic, and water quality issues in the long run.• Bioswale (BS) showed the smallest footprint requirement among all BMPs analyzed.Only 0.7 ha of BS installation can be adequate to hold the runoff generated from a maximum of 4.5 ha of commercial development (84% bigger than the BS footprint itself).Therefore, the implementation of bioswale alone can considerably reduce the footprint and construction cost.• Bioretention (BR) can also be promising with a preferable footprint of 0.93 ha, which might cover a commercial development as large as 5 ha with a 90-95% runoff reduction.

•
Apart from providing good serviceability and enhancing property value, porous concrete pavement (PCP) might be capable of detaining 34-50% of runoff generated outside its footprint boundary.For example, 4 ha of PCP installation might serve a commercial development as large as 6.1 ha.

•
Considering all hydrological benefits and costs for a single BMP case, BR alone might be the optimal option from any commercial development in the LRGV region.

•
Considering all hydrological benefits, land use management, aesthetics, and costs, a combination of BR, BS, and WP may be the optimal option, with almost 100% runoff reduction from the site.However, when LIDs are coupled with a wet detention pond, it may require a larger footprint area than that expected from WP alone.
The DSS will benefit the public and the environment by providing innovative and critical information to the region to meet the challenging requirements of balancing floods, water quality demands, and public safety needs.Reducing localized flooding and identifying sustainable, economical, and innovative stormwater management strategies are key benefits projected to result from the development of the DSS.The successful implementation and adoption of the DSS in urban stormwater management should improve the stormwater runoff water quality within the LRGV.
Further studies for the improvement of the DSS may include a detailed cost-benefit analysis of BMPs, based on the analyzed costs in terms of construction and O&M costs of the BMPs, and benefits in terms of runoff reduction, property value increases, leasehold revenue, flood damage mitigation, and groundwater recharge.Future studies can be recommended to incorporate more BMP options (e.g., green roof, dry pond, subsurface wetland, rainwater harvesting cistern, etc.) into the DSS algorithm.Enhancement of the DMT tool can be suggested in terms of providing user-defined drainage policy constraints.Moreover, the current DSS convergence algorithm determines BMP(s) footprint by dealing with a 50-year frequency storm event (12.7 mm/h) at a post-development rate as input and with a 10-year frequency storm (7.62 mm/h) discharge at a predevelopment rate as an output.Further studies can be done to increase its robustness so that the users can select their design storm frequency from the drop-down menu of the tool.This finding might help to assess long-term runoff reductions by dealing with varying rainfall characteristics from multiple events.Since WinSLAMM approximates a composite peak flow from the first flush runoff generated with complex drainage surfaces (impervious, pervious, semi-pervious, etc.), our DSS might not reflect a continuous rainfall-runoff characteristic (with double or more peaks) to determine the optimal BMP alternative.For the development of BMP models, future works can be recommended to emphasize the calibration of the hydrograph based on the user-defined hyetograph to capture the characteristics of continuous rainfall.The findings of our developed DSS will only be applicable to the LRGV region since the algorithm allows hydrogeological input from 14 cities in the region.However, a similar methodology and convergence algorithm can be applied to other impaired watersheds in the US and the rest of the world only after a regional calibration and validation of BMP models.

•Figure 1 .
Figure 1.Map of the Arroyo Colorado Watershed and site locations of monitored porous concrete pavement, bioretention, and bioswale within the Lower Rio Grande Valley in Texas [13].

Figure 1 .
Figure 1.Map of the Arroyo Colorado Watershed and site locations of monitored porous concrete pavement, bioretention, and bioswale within the Lower Rio Grande Valley in Texas [13].

Figure 2 .
Figure 2. The framework of the decision-support system (DSS) demonstrating input parameters, methodology, and outcomes of the tool (note: all parallelogram-shaped boxes represent input/output parameters, and all square-shaped boxes represent a process).

Figure 2 .
Figure 2. The framework of the decision-support system (DSS) demonstrating input parameters, methodology, and outcomes of the tool (note: all parallelogram-shaped boxes represent input/output parameters, and all square-shaped boxes represent a process).

Figure 4 .
Figure 4. Graphical representation of DSS-simulated results to assess the effect of the size of commercial development on the recommendation of LID footprints and controlled peak discharge: (a) wet detention pond (WP), (b) porous concrete pavement (PCP), (c) bioretention, and (d) bioswale.

Figure
Figure4Areveals the size of the wet detention pond (WP), which might be adequate in runoff management for a range of development sizes (0.1-6.0 ha).Results suggested that the 2 ha development size might require the building of a 0.8 ha wet pond that would induce the least peak discharge of 0.034 m 3 /s and a maximum runoff reduction of 94% from the site.The installation of 0.8

Figure 4 .
Figure 4. Graphical representation of DSS-simulated results to assess the effect of the size of commercial development on the recommendation of LID footprints and controlled peak discharge: (A) wet detention pond (WP), (B) porous concrete pavement (PCP), (C) bioretention, and (D) bioswale.

Sustainability 2020 ,
12,  x FOR PEER REVIEW 12 of 19 in series.The comparison of different BMP combinations in terms of footprint requirements for 4 ha of development is shown in Figure5.DSS inputs were fixed to a 4 ha commercial development in the city of Brownsville, TX.The figure distinguishes the distribution of footprint area between BMPs in the form of color-coded bar graphs.During the simulation, the 50-year uncontrolled peak discharge was used as an input to estimate the total inflow volume and runoff reduction (by subtracting the LID-controlled runoff volume).The figure also depicts the performance in terms of percentage runoff reduction after the implementation of the BMPs in the series.

Figure 5 .
Figure 5.The comparison of different BMP configurations in terms of footprint requirements in 4 ha of development.

Figure 5 .
Figure 5.The comparison of different BMP configurations in terms of footprint requirements in 4 ha of development.

Table 1 .
Summary of the monitoring period, the number of rainfall events, max rainfall depth, calibrated parameters, and performance indexes of the period for the three BMPs.

Table 2 .
DSS-simulated results summary for different types of control practices and sizes of commercial developments.