Modeling watershed-scale effectiveness of agricultural best management practices to reduce phosphorus loading

Planners advocate best management practices (BMPs) to reduce loss of sediment and nutrients in agricultural areas. However, the scientific community lacks tools that use readily available data to investigate the relationships between BMPs and their spatial locations and water quality. In rural, humid regions where runoff is associated with saturation-excess processes from variable source areas (VSAs), BMPs are potentially most effective when they are located in areas that produce the majority of the runoff. Thus, two critical elements necessary to predict the water quality impact of BMPs include correct identification of VSAs and accurate predictions of nutrient reduction due to particular BMPs. The objective of this research was to determine the effectiveness of BMPs using the Variable Source Loading Function (VSLF) model, which captures the spatial and temporal evolutions of VSAs in the landscape. Data from a long-term monitoring campaign on a 164-ha farm in the New York City source watersheds in the Catskills Mountains of New York state were used to evaluate the effectiveness of a range of BMPs. The data spanned an 11-year period over which a suite of BMPs, including a nutrient management plan, riparian buffers, filter strips and fencing, was installed to reduce phosphorus (P) loading. Despite its simplicity, VSLF predicted the spatial distribution of runoff producing areas well. Dissolved P reductions were simulated well by using calibrated reduction factors for various BMPs in the VSLF model. Total P losses decreased only after cattle crossings were installed in the creek. The results demonstrated that BMPs, when sited with respect to VSAs, reduce P loss from agricultural watersheds, providing useful information for targeted water quality management.     

Flood Risk Reduction from Agricultural Best Management Practices

Best management practices (BMPs) play an important role in improving impaired water quality from conventional row crop agriculture. In addition to reducing nutrient and sediment loads, BMPs such as fertilizer management, reduced tillage, and cover crops could alter the hydrology of agricultural systems and reduce surface water runoff. While attention is devoted to the water quality benefits of BMPs, the potential co‐benefits of flood loss reduction are often overlooked. This study quantifies the effects of selected commonly applied BMPs on expected flood loss to agricultural and urban areas in four Iowa watersheds. The analysis combines a watershed hydrologic model, hydraulic model outputs, and a loss estimation model to determine relationships between hydrologic changes from BMP implementations and annual economic flood loss. The results indicate a modest reduction in peak discharge and economic loss, although loss reduction is substantial when urban centers or other high‐value assets are located downstream in the watershed. Among the BMPs, wetlands, and cover crops reduce losses the most. The research demonstrates that watershed‐scale implementation of agricultural BMPs could provide benefits of flood loss reduction in addition to water quality improvements.     

MODELING RAINFALL-RUNOFF PROCESS INCLUDING LOSSES1

ABSTRACT: Existing discrete, linear rainfall-runoff models generally require the effective rainfall of a given storm as the input for computing the runoff hydrograph. This paper proposes a method for estimating, simultaneously, the optimal values of model parameters and the rainfall losses frem the measured rainfall hyetograph and the runoff hydrograph. The method involves an ARMA model for the rainfall-runoff process and a nonlinear iterative technique. The number of model parameters to be estimated for the ARMA model is much less than the unit hydrograph model. Applications of the model to three different watersheds show that the computed runoff hydrographs agree well with the measurements.     

Economic analysis of best management practices to reduce watershed phosphorus losses

In phosphorus-limited freshwater systems, small increases in phosphorus (P) concentrations can lead to eutrophication. To reduce P inputs to these systems, various environmental and agricultural agencies provide producers with incentives to implement best management practices (BMPs). In this study, we examine both the water quality and economic consequences of systematically protecting saturated, runoff-generating areas from active agriculture with selected BMPs. We also examine the joint water quality/economic impacts of these BMPs-specifically BMPs focusing on barnyards and buffer areas. Using the Variable Source Loading Function model (a modified Generalized Watershed Loading Function model) and net present value analysis (NPV), the results indicate that converting runoff-prone agricultural land to buffers and installing barnyard BMPs are both highly effective in decreasing dissolved P loss from a single-farm watershed, but are also costly for the producer. On average, including barnyard BMPs decreases the nutrient loading by about 5.5% compared with only implementing buffers. The annualized NPV for installing both buffers on only the wettest areas of the landscape and implementing barnyard BMPs becomes positive only if the BMPs lifetime exceeds 15 yr. The spatial location of the BMPs in relation to runoff producing areas, the time frame over which the BMPs are implemented, and the marginal costs of increasing buffer size were found to be the most critical considerations for water quality and profitability. The framework presented here incorporates estimations of nutrient loading reductions in the economic analysis, and is applicable to farms facing BMP adoption decisions.     

PARAMETER UNCERTAINTY IN RAINFALL-RUNOFF MODELING FOR POLDER SYSTEM DESIGN1

ABSTRACT: An approach is developed for incorporating the uncertainty of parameters for estimating runoff in the design of polder systems in ungaged watersheds. Monte Carlo Simulation is used to derive a set of realizations of streamflow hydrographs for a given design rainstorm using the U. S. Soil Conservation Service (SCS) unit hydrograph model. The inverse of the SCS curve number, which is a function of the antecedent runoff condition in the SCS model, is the random input in the Monte Carlo Simulation. Monte Carlo realizations of streamfiow hydrographs are used to simulate the performance of a polder flood protection system. From this simulation the probability of occurrence of flood levels for a particular hydraulic design may be used to evaluate its effectiveness. This approach is demonstrated for the Pluit Polder flood protection system for the City of Jakarta, Indonesia. While the results of the application indicate that uncertainty in the antecedent runoff condition is important, the effects of uncertainty in rainfall data, in additional runoff parameters, such as time to peak, in the hydraulic design, and in the rainfall-runoff model selected should also be considered. Although, the SCS model is limited to agricultural conditions, the approach presented herein may be applied to other flood control systems if appropriate storm runoff models are selected.     

Optimal nonpoint source pollution control strategies for a reservoir watershed in Taiwan

The purpose of this study is to develop a model for optimal nonpoint source pollution control for the Fei-Tsui Reservoir watershed in Northern Taiwan. Several structural best management practices (BMPs) are selected to treat stormwater runoff. The complete model consists of two interacting components: an optimization model based on discrete differential dynamic programming (DDDP) and a zero-dimensional reservoir water quality model. A predefined procedure is used to locate suitable sites for construction of various selected BMPs in the watershed. In the optimization model, the objective function is to find the best combination of BMP type and placement, which minimizes the total construction and operation, maintenance, and repair (OMR) costs of the BMPs. The constraints are the water quality standards for total phosphorus (TP) and total suspended solids (TSS) concentrations in the reservoir. A zero-dimensional reservoir water quality model of the Vollenweider type is embedded in the optimization framework to simulate pollutant concentrations in Fei-Tsui Reservoir. The resulting optimal cost and benefit of water quality improvement are depicted by the model-derived trade-off curves. The modeling framework developed in the present study could be used as an efficient tool for planning a watershed-wide implementation of BMPs for mitigating stormwater pollution impact on the receiving water bodies.     

UNCERTAINTY ANALYSIS OF BMP EFFECTWENESS FOR CONTROLLING NITROGEN FROM URBAN NONPOINT SOURCES1

ABSTRACT: The effectiveness of urban Best Management Practices (BMPs) in achieving the No-Net-Increase Policy (NNTP), a policy designed to limit nonpoint nitrogen loading to Long Island Sound (US), is analyzed. A unit loading model is used to simulate annual nitrogen exported from the Norwalk River watershed (Connecticut) under current and future conditions. A probabilistic uncertainty analysis is used to incorporate uncertainty in nitrogen export coefficients and BMP nitrogen removal effectiveness. The inclusion of uncertainty in BMP effectiveness and nitrogen export coefficients implies that additional BMPs, or BMPs with a greater effectiveness in nitrogen removal, will be required to achieve the NNIP. Even though including uncertainty leads to an increase in BMP implementation rates or BMP effectiveness, this type of analysis provides the decision maker with a more realistic assessment of the likelihood that implementing BMPs as a management strategy will be successful. Monte Carlo simulation results indicate that applying BMPs to new urban developments alone will not be sufficient to achieve the NNIP since BMPs are not 100 percent effective in removing the increase in nitrogen caused by urbanization. BMPs must also be applied to selected existing urban areas. BMPs with a nitrogen removal effectiveness of 40–60 percent, probably the highest level of removal that can be expected over an entire watershed, must be applied to at least 75 percent of the existing urban area to achieve the NNIP This high rate of application is not likely to be achieved in urbanized watersheds in the LIS watershed; therefore, additional point source control will be necessary to achieve the NNIP     

Farm factors associated with reducing Cryptosporidium loading in storm runoff from dairies

A systems approach was used to evaluate environmental loading of Cryptosporidium oocysts on five coastal dairies in California. One aspect of the study was to determine Cryptosporidium oocyst concentrations and loads for 350 storm runoff samples from dairy high use areas collected over two storm seasons. Selected farm factors and beneficial management practices (BMPs) associated with reducing the Cryptosporidium load in storm runoff were assessed. Using immunomagnetic separation (IMS) with direct fluorescent antibody (DFA) analysis, Cryptosporidium oocysts were detected on four of the five farms and in 21% of storm runoff samples overall. Oocysts were detected in 59% of runoff samples collected near cattle less than 2 mo old, while 10% of runoff samples collected near cattle over 6 mo old were positive. Factors associated with environmental loading of Cryptosporidium oocysts included cattle age class, 24 h precipitation, and cumulative seasonal precipitation, but not percent slope, lot acreage, cattle stocking number, or cattle density. Vegetated buffer strips and straw mulch application significantly reduced the protozoal concentrations and loads in storm runoff, while cattle exclusion and removal of manure did not. The study findings suggest that BMPs such as vegetated buffer strips and straw mulch application, especially when placed near calf areas, will reduce environmental loading of fecal protozoa and improve stormwater quality. These findings are assisting working dairies in their efforts to improve farm and ecosystem health along the California coast.     

Modeling Watershed‐Scale Impacts of Stormwater Management with Traditional versus Low Impact Development Design

Stormwater runoff and associated pollutants from urban areas in the greater Chesapeake Bay Watershed (CBW) impair local streams and downstream ecosystems, despite urbanized land comprising only 7% of the CBW area. More recently, stormwater best management practices (BMPs) have been implemented in a low impact development (LID) manner to treat stormwater runoff closer to its source. This approach included the development of a novel BMP model to compare traditional and LID design, pioneering the use of comprehensively digitized storm sewer infrastructure and BMP design connectivity with spatial patterns in a geographic information system at the watershed scale. The goal was to compare total watershed pollutant removal efficiency in two study watersheds with differing spatial patterns of BMP design (traditional and LID), by quantifying the improved water quality benefit of LID BMP design. An estimate of uncertainty was included in the modeling framework by using ranges for BMP pollutant removal efficiencies that were based on the literature. Our model, using Monte Carlo analysis, predicted that the LID watershed removed approximately 78 kg more nitrogen, 3 kg more phosphorus, and 1,592 kg more sediment per square kilometer as compared with the traditional watershed on an annual basis. Our research provides planners a valuable model to prioritize watersheds for BMP design based on model results or in optimizing BMP selection.     

Combined monitoring and modeling indicate the most effective agricultural best management practices

Although water quality problems associated with agricultural nonpoint source (NPS) pollution have prompted the rapid and widespread adoption of a variety of so called "best management practices" (BMPs), there have been few realistic efforts to assess their combined effectiveness in reducing NPS pollution. This study used the Variable Source Loading Function (VSLF) model, a distributed watershed model, to simulate phosphorus (P) loading from an upstate New York dairy farm before and after the implementation of a suite of BMPs. With minimal calibration, the model calculates the dissolved P (DP) losses from impervious surfaces (e.g., barnyards), the plant/soil complex, field-applied manure, and loads associated with baseflow conditions. The simulated DP loads agreed well with measured loads for both the pre-BMP and post-BMP periods. More importantly, results showed that BMPs reduced DP loads by 35%, which is over half of the expected reduction if all manure was removed from the watershed, i.e., approximately 50% reduction. The model results indicate that had no BMPs been installed DP loads would be approximately 37% greater than observed at the watershed outlet. The most effective BMPs were those that disassociated pollutant loading areas from areas prone to generating runoff, i.e., hydrologically sensitive areas. By contrast, attempts to reduce P content in manure were somewhat less effective. This study demonstrates that a combination of distributed, mechanistic modeling and long-term monitoring provides better insights into the effectiveness of water quality protection efforts than either individually.     

Effects of Impervious Area and BMP Implementation and Design on Storm Runoff and Water Quality in Eight Small Watersheds

The effects of increases in effective impervious area (EIA) and the implementation of water quality protection designed detention pond best management practices (BMPs) on storm runoff and stormwater quality were assessed in Gwinnett County, Georgia, for the period 2001‐2008. Trends among eight small watersheds were compared, using a time trend study design. Significant trends were detected in three storm hydrologic metrics and in five water quality constituents that were adjusted for variability in storm characteristics and climate. Trends in EIA ranged from 0.10 to 1.35, and changes in EIA treated by BMPs ranged from 0.19 to 1.32; both expressed in units of percentage of drainage area per year. Trend relations indicated that for every 1% increase in watershed EIA, about 2.6, 1.1, and 1.5% increases in EIA treated by BMPs would be required to counteract the effects of EIA added to the watersheds on peak streamflow, stormwater yield, and storm streamflow runoff, respectively. Relations between trends in EIA, BMP implementation, and water quality were counterintuitive. This may be the result of (1) changes in constituent inputs in the watersheds, especially downstream of areas treated by BMPs; (2) BMPs may have increased the duration of stormflow that results in downstream channel erosion; and/or (3) spurious relationships between increases in EIA, BMP implementation, and constituent inputs with development rates.     

Application of a Multi-Objective Optimization Method to Provide Least Cost Alternatives for NPS Pollution Control

Nonpoint source (NPS) pollutants such as phosphorus, nitrogen, sediment, and pesticides are the foremost sources of water contamination in many of the water bodies in the Midwestern agricultural watersheds. This problem is expected to increase in the future with the increasing demand to provide corn as grain or stover for biofuel production. Best management practices (BMPs) have been proven to effectively reduce the NPS pollutant loads from agricultural areas. However, in a watershed with multiple farms and multiple BMPs feasible for implementation, it becomes a daunting task to choose a right combination of BMPs that provide maximum pollution reduction for least implementation costs. Multi-objective algorithms capable of searching from a large number of solutions are required to meet the given watershed management objectives. Genetic algorithms have been the most popular optimization algorithms for the BMP selection and placement. However, previous BMP optimization models did not study pesticide which is very commonly used in corn areas. Also, with corn stover being projected as a viable alternative for biofuel production there might be unintended consequences of the reduced residue in the corn fields on water quality. Therefore, there is a need to study the impact of different levels of residue management in combination with other BMPs at a watershed scale. In this research the following BMPs were selected for placement in the watershed: (a) residue management, (b) filter strips, (c) parallel terraces, (d) contour farming, and (e) tillage. We present a novel method of combing different NPS pollutants into a single objective function, which, along with the net costs, were used as the two objective functions during optimization. In this study we used BMP tool, a database that contains the pollution reduction and cost information of different BMPs under consideration which provides pollutant loads during optimization. The BMP optimization was performed using a NSGA-II based search method. The model was tested for the selection and placement of BMPs in Wildcat Creek Watershed, a corn dominated watershed located in northcentral Indiana, to reduce nitrogen, phosphorus, sediment, and pesticide losses from the watershed. The Pareto optimal fronts (plotted as spider plots) generated between the optimized objective functions can be used to make management decisions to achieve desired water quality goals with minimum BMP implementation and maintenance cost for the watershed. Also these solutions were geographically mapped to show the locations where various BMPs should be implemented. The solutions with larger pollution reduction consisted of buffer filter strips that lead to larger pollution reduction with greater costs compared to other alternatives.     

A systematic approach for the comparative assessment of stormwater pollutant removal potentials

This paper describes the development of a methodology to theoretically assess the stormwater pollutant removal performances of structural best management practices (BMPs). The method combines the categorisation of the relative importance of the primary removal processes within 15 different BMPs with an evaluation of the ability of each process to remove a pollutant in order to generate a value representing the pollutant removal potential for each BMP. The methodology is demonstrated by applying it separately to a set of general water quality indicators (total suspended solids, biochemical and chemical oxygen demand, nitrates, phosphates and faecal coliforms) to produce a ranked list of BMP pollutant removal efficiencies. Given the limited amount of available monitoring data relating to the differential pollutant removal capabilities of BMPs, the resulting prioritization will support stakeholders in making urban drainage decisions from the perspective of pollutant removal. It can also provide inputs to existing urban hydrology models, which aim to predict the treatment performances of BMPs. The level of resilience of the proposed approach is tested using a sensitivity analysis and the limitations in terms of BMP design and application are discussed.     

EXPORT COEFFICIENT MODELING TO ASSESS PHOSPHORUS LOADING IN AN URBAN WATERSHED1

ABSTRACT: An export coefficient modeling approach was used to assess the influence of land use on phosphorus loading to a Southern Ontario stream. A model was constructed for the 1995–1996 water year and calibrated within ± 3 percent of the observed mean concentration of total phosphorus. It was found that runoff from urban areas contributed most to the loading of phosphorus to the stream. When the model was assessed by running it for the 1977–1978 water year, using water quality and land use data collected independently, agreement within ± 7 percent was obtained. The model was then used to forecast the impact of future urban development proposed for the watershed, in terms of phosphorus loading, and to evaluate the reduction in loading resulting from several urban best management practices (BMP). It was determined that phosphorus removal will have to be applied to all the urban runoff from the watershed to appreciably reduce stream phosphorus concentration. Of the BMP designs assessed, an infiltration pond system resulted in the greatest phosphorus load reduction, 50 percent from the 1995–1996 baseline.     

Management of Urban Road Runoff Containing PAHs: Probabilistic Modeling and Its Application in Beijing,China1

Abstract: Pollutant loading from storm runoff is considered to be an important component of nonpoint source pollution in urban areas. In developing countries, because of the accelerated urbanization and motorization, storm runoff pollution has become a challenge for improving aquatic environmental quality. An effective storm runoff management plan needs to be developed, and questions concerning how much and which proportion of a storm should be treated need to be answered. In this study, a model is developed to determine the fraction of storm runoff that needs to be treated to meet the discharge standard within a given probability. The model considers that the pollutants can be mobilized during the early stage of a storm. The model is applied to a field study of polycyclic aromatic hydrocarbons (PAHs) in road runoff in Beijing, China. In this case, the probability that the PAH load will be mobilized with suspended sediments by the earlier portion of the flush is 73%. Given the high PAH loading in the study area and the referenced discharge standard, the probability that the entire runoff should be captured and treated is 94%. Thus, urban planners need to consider treatment systems for the majority of the storms in this area, whether the PAH load is in the first flush or not. This methodology can be applied to other regions where PAH loads may result in different management outcomes.     

WATER LOSS EQUATIONS AND COEFFICIENTS TO ESTIMATE RUNOFF FROM RAINFALL1

Infiltration models are based on physical characteristics of the soil and initial soil moisture. For a given soil it is based on the initial soil moisture distribution. A computer simulation model for flood runoff systems (FH-Model) was used to analyze 39 sets of rainfall-runoff data on four small watersheds ranging in size from 17 to 342 square kilometers located in the Yamaska River basin in Quebec. From these analyses, parameters and coefficients have been determined for a water loss (infiltration) equation. A method for determining the loss parameters, using a nonlinear least square curve fitting technique, is presented. Expressions were made to relate the loss parameters to antecedent precipitation. The equations were tested on 11 storm rainfall and runoff events on a watershed located in the same region and close agreements were found.     

Simulation of Combined Best Management Practices and Low Impact Development for Sustainable Stormwater Management1

Damodaram, Chandana, Marcio H. Giacomoni, C. Prakash Khedun, Hillary Holmes, Andrea Ryan, William Saour, and Emily M. Zechman, 2010. Simulation of Combined Best Management Practices and Low Impact Development for Sustainable Stormwater Management. Journal of the American Water Resources Association (JAWRA) 1-12. DOI: 10.1111/j.1752-1688.2010.00462.x Abstract: Urbanization causes increased stormwater runoff volumes, leading to erosion, flooding, and the degradation of instream ecosystem health. Although Best Management Practices (BMPs) are used widely as a means for controlling flood runoff events, Low Impact Development (LID) options have been proposed as an alternative approach to better mimic the natural flow regime by using decentralized designs to control stormwater runoff at the source, rather than at a centralized location in the watershed. For highly urbanized areas, LID practices such as rainwater harvesting, green roofs, and permeable pavements can be used to retrofit existing infrastructure and reduce runoff volumes and peak flows. This paper describes a modeling approach to incorporate these LID practices in an existing hydrologic model to estimate the effects of LID choices on streamflow. The modeling approach has been applied to a watershed located on the campus of Texas A&M University in College Station, Texas, to predict the stormwater reductions resulting from retrofitting existing infrastructure with LID technologies. Results demonstrate that use of these LID practices yield significant stormwater control for small events and less control for flood events. A combined BMP-LID approach is tested for runoff control for both flood and frequent rainfall events.     

The Ability of Urban Residential Lawns to Disconnect Impervious Area from Municipal Sewer Systems1

Abstract: Runoff from urban catchments depends largely on the amount of impervious surface and the connectivity of these surfaces to the storm sewer drainage system. In residential areas, pervious lawns can be used to help manage stormwater runoff by intercepting and infiltrating runoff from impervious surfaces. The goal of this research was to develop and evaluate a simple method for estimating the reduction in stormwater runoff that results when runoff from an impervious surface (e.g., rooftop) is directed onto a pervious surface (e.g., lawn). Fifty‐two stormwater runoff reduction tests were conducted on six residential lawns in Madison, Wisconsin during the summer of 2004. An infiltration‐loss model that requires inputs of steady‐state infiltration rate, abstraction (defined here as surface storage, vegetation interception and cumulative total infiltration minus steady‐state infiltration during the period prior to steady‐state), and inundated area was evaluated using experimental data. The most accurate results were obtained using the observed steady‐state infiltration rates and inundated areas for each test, combined with a constant abstraction for all tests [root mean squared (RMS) difference = 1.0 cm]. A second case utilized lawn‐averaged steady‐state infiltration rates, a regression estimate of inundated area based on flow‐path length, and lawn‐specific abstractions based on infiltration rate (RMS difference = 2.2 cm). In practice, infiltration rates will likely be determined using double‐ring infiltration measurements (RMS difference = 3.1 cm) or soil texture (RMS difference = 5.7 cm). A generalized form of the model is presented and used to estimate annual stormwater runoff volume reductions for Madison. Results indicate the usefulness of urban lawns as a stormwater management practice and could be used to improve urban runoff models that incorporate indirectly connected impervious areas.     

A STRATEGY FOR IMPLEMENTING BMPs FOR CONTROLLING NONPOINT SOURCE POLLUTION: THE CASE OF THE FEI‐TSUI RESERVOIR WATERSHED IN TAIWAN1

ABSTRACT: This paper describes a concerted effort by the Taiwan Water Resources Bureau, the City of Taipei, and the Bureau of Fei‐tsui Reservoir Management to protect the water quality in the Fei‐tsui Reservoir.The reservoir is the major source of water supply for over two million people in the metropolitan area of Taipei. Over the years the reservoir has suffered from siltation and more recently eutrophication. The sources of the pollution are traced to the hundreds of tea gardens, rice fields and other agricultural areas in the watershed and to urban sources such as construction sites. Large amounts of nutrients enter the reservoir by way of storm water runoff during storm or typhoon events. Since 1999, various agencies have worked to initiate an effort to reduce nonpoint pollution in the Fei‐Tsui Reservoir watershed. Practices being considered include nonstructural measures such as nutrient management, and structural measures such as swales, detention basins, and wetlands, in addition to erosion and sediment control methods. A number of field tests have been completed on the performance of selected best management practices (BMPs). A strategy for implementing the BMPs at the watershed scale has been developed based on a total maximum daily load (TMDL) analysis that is reported in this paper.