Assessment of a turfgrass sod best management practice on water quality in a suburban watershed

The disposal of manure on agricultural land has caused water quality concerns in many rural watersheds, sometimes requiring state environmental agencies to conduct total maximum daily load (TMDL) assessments of stream nutrients, such as nitrogen (N) and phosphorus (P). A best management practice (BMP) has been developed in response to a TMDL that mandates a 50% reduction of annual P load to the North Bosque River (NBR) in central Texas. This BMP exports composted dairy manure P through turfgrass sod from the NBR watershed to urban watersheds. The manure-grown sod releases P slowly and would not require additional P fertilizer for up to 20 years in the receiving watershed. This would eliminate P application to the sod and improve the water quality of urban streams. The soil and water assessment tool (SWAT) was used to model a typical suburban watershed that would receive the sod grown with composted dairy manure to assess water quality changes due to this BMP. The SWAT model was calibrated to simulate historical flow and estimated sediment and nutrient loading to Mary's Creek near Fort Worth, Texas. The total P stream loading to Mary's Creek was lower when manure-grown sod was transplanted instead of sod grown with inorganic fertilizers. Flow, sediment and total N yield were the same for both cases at the watershed outlet. The SWAT simulations indicated that the turfgrass BMP can be used effectively to import manure P into an urban watershed and reduce in-stream P levels when compared to sod grown with inorganic fertilizers.     

GENERALIZED WATERSHED LOADING FUNCTIONS FOR STREAM FLOW NUTRIENTS1

ABSTRACT: Loading functions are proposed as a general model for estimating monthly nitrogen and phosphorus fluxes in stream flow. The functions have a simple mathematical structure, describe a wide range of rural and urban nonpoint sources, and couple surface runoff and ground water discharge. Rural runoff loads are computed from daily runoff and erosion and monthly sediment yield calculations. Urban runoff loads are based on daily nutrient accumulation rates and exponential wash off functions. Ground water discharge is determined by lumped parameter unsaturated and saturated zone soil moisture balances. Default values for model chemical parameters were estimated from literature values. Validation studies over a three-year period for an 850 km² watershed showed that the loading functions explained at least 90 percent of the observed monthly variation in dissolved and total nitrogen and phosphorus fluxes in stream flow. Errors in model predictions of mean monthly fluxes were: dissolved phosphorus - 4 percent; total phosphorus - 2 percent; dissolved nitrogen - 18 percent; and total nitrogen - 28 percent. These results were obtained without model calibration.     

STREAM DISCHARGE PREDICTION VIA A GRID BASED SOIL WATER ROUTING WITH PADDY FIELDS1

ABSTRACT: A grid based daily hydrologic model for a watershed with paddy fields was developed to predict the stream discharge. ASCII formatted elevation, soil, and land use data supported by the GRASS Geographic Information System are used to generate distributed results such as surface runoff and subsurface flow, soil water content, and evapotranspiration. The model uses a single flow path algorithm and simulates a water balance at each grid element. A linear reservoir assumption was used to predict subsurface runoff components. The model was applied to a 75.6 km² watershed located in the middle of South Korea, and observed stream flow hydrographs from 1995 and 1996 were compared to model predictions. The stream flow predictions of 1995 and 1996 generally agreed with the observed flow, resulting in a Nash‐Sutcliffe efficiency R² of 0.60 and 0.62, respectively. The hydraulic conductivity for percolating water through the saturated layer affected baseflow generation. The levee height of the paddy influenced the time and magnitude of the surface runoff, depending on irrigation management. The model will be used for making low flow management decisions by evaluating the role of each land use to stream flow, especially in case of paddy decrease by gradual urbanization of a watershed.     

MODELING FRAMEWORK FOR MANAGING COPPER RUNOFF IN URBAN WATERSHEDS1

ABSTRACT: A modeling framework was developed for managing copper runoff in urban watersheds that incorporates water quality characterization, watershed land use areas, hydrologic data, a statistical simulator, a biotic ligand binding model to characterize acute toxicity, and a statistical method for setting a watershed specific copper loading. The modeling framework is driven by export coefficients derived from water quality parameters and hydrologic inputs measured in an urban watershed's storm water system. This framework was applied to a watershed containing a copper roof built in 1992. A series of simulations was run to predict the change in receiving stream water chemistry caused by roof aging and to determine the maximum copper loading (at the 99 percent confidence level) a watershed could accept without causing acute toxicity in the receiving stream. Forecasting the amount of copper flux responsible for exceeding the assimilation capacity of a watershed can be directly related to maximum copper loadings responsible for causing toxicity in the receiving streams. The framework developed in this study can be used to evaluate copper utilization in urban watersheds.     

Modeling the relationship between development and storm water and nutrient runoff

The EPA Storm Water Management Model was used to model the effects of urban and agricultural development on storm water runoff from uplands bordering a Louisiana swamp forest. Using this model, we examined the effects of changing land use patterns. By 1995 it is projected that urban land on the uplands bordering the swamp will increase by 321 percent, primarily at the expense of land currently in agriculture. Simulation results indicate that urbanization will cause storm water runoff rates to be up to 4.2 times greater in 1995 than in 1975. Nutrient runoff will increase 28 percent for nitrogen (N) and 16 percent for phosphorus (P) during the same period. The environmental effects of these changes in the receiving swamp forest are examined.     

LOADING FUNCTIONS FOR PREDICTING NUTRIENT LOSSES FROM COMPLEX WATERSHEDS1

ABSTRACT: A loading function methodology is presented for predicting runoff, sediment, and nutrient losses from complex watersheds. Separate models are defined for cropland, forest, urban and barnyard sources, and procedures for estimating baseflow nutrients are provided. The loading functions are designed for use as a preliminary screening tool to isolate the major contributors in a watershed. Input data sources are readily available and the functions do not require costly calibrations. Data requirements include watershed land use and soil information, daily precipitation and temperature records and rainfall erosivities. Comparison of predicted and measured water, sediment, and nutrient runoff fluxes for the West Branch Deleware River in New York, indicated that runoff was underpredicted by about 14 percent while dissolved nutrients were within 30 percent of observed values. Sediment and solid-phase nutrients were overpredicted by about 50 percent. An annual nutrient budget for the West Branch Delaware River showed that cornland was the major source of sediment, solid phase nutrients, and total phosphorus. Waste water treatment plants and ground water discharge contributed the most dissolved phosphorus and dissolved nitrogen, respectively.     

GIS BASED LONG TERM HYDROLOGIC IMPACT EVALUATION FOR WATERSHED URBANIZATION1

ABSTRACT: To adequately manage impacts of ongoing or future land use changes in a watershed, the magnitude of their hydrologic impacts needs to be assessed. A grid based daily streamflow model was calibrated with two years of observed streamflow data, using time periods when land use data are available and verified by comparison of model predictions with observed streamflow data. Streamflow data were separated into direct runoff and baseflow to estimate the impacts of urbanization on each hydrologic component. Analysis of the ratio between direct runoff and total runoff from 30 years of simulation results and the change in these ratios with urbanization shows that estimated annual direct runoff increased from 49.2 percent (1973) to 63.1 percent (1984) and 65.0 percent (1991), indicating the effects of urbanization are greater on direct runoff than on total runoff. The direct runoff ratio also varies with annual rainfall, with dry year ratios larger than those for wet years. This suggests that the impact of urbanization on areas that are sensitive to runoff ratios, such as stream ecosystems, might be more serious during drier years than in wetter years in terms of water quality and water yield. This indicates that sustainable base‐flow is important to maintaining sound stream ecosystems.     

MOBILITY AND AQUATIC TOXICITY OF COPPER IN AN URBAN WATERSHED1

ABSTRACT: Abundant use of copper based products has resulted in increased violation of copper water quality criteria in runoff from urban storm water systems. The objectives of this work were to understand the mobility and toxicity of copper in an urban watershed and to apportion the amount of copper entering the freshwater receiving stream from different urban land covers using a mass balance approach. Sixteen rainfall events collected from the University of Connecticut study watershed between August 1998 and September 2000 were analyzed to assess copper flux in an urban storm water system. Mean flow weighted dissolved copper concentrations observed in the study for copper based architectural material runoff, pervious area runoff, impervious area runoff, and in the receiving stream were 1210 ± 840, 9 ± 3, 8 ± 2, and 14 ± 7 μg/L, respectively. Mean dissolved copper concentrations in the receiving stream exceeded Connecticut's water quality criteria. Despite exceeding the dissolved concentration based criteria, cupric ion concentrations at the system outlet remained below 0.05 μg/L for all storms analyzed, and no acute toxicity (using Daphnia pulex as the test organism) was measured in samples collected from the stream.     

THE IMPACT OF RUNOFF GENERATION MECHANISMS ON THE LOCATION OF CRITICAL SOURCE AREAS1

ABSTRACT: Identifying phosphorus (P) source areas and transport pathways is a key step in decreasing P loading to natural water systems. This study compared the effects of two modeled runoff generation processes ‐ saturation excess and infiltration excess ‐ on total phosphorus (TP) and soluble reactive phosphorus (SRP) concentrations in 10 catchment streams of a Catskill mountain watershed in southeastern New York. The spatial distribution of runoff from forested land and agricultural land was generated for both runoff processes; results of both distributions were consistent with Soil Conservation Service‐Curve Number (SCS‐CN) theory. These spatial runoff distributions were then used to simulate stream concentrations of TP and SRP through a simple equation derived from an observed relation between P concentration and land use; empirical results indicate that TP and SRP concentrations increased with increasing percentage of agricultural land. Simulated TP and SRP stream concentrations predicted for the 10 catchments were strongly affected by the assumed runoff mechanism. The modeled TP and SRP concentrations produced by saturation excess distribution averaged 31 percent higher and 42 percent higher, respectively, than those produced by the infiltration excess distribution. Misrepresenting the primary runoff mechanism could not only produce erroneous concentrations, it could fail to correctly locate critical source areas for implementation of best management practices. Thus, identification of the primary runoff mechanism is critical in selection of appropriate models in the mitigation of nonpoint source pollution. Correct representation of runoff processes is also critical in the future development of biogeochemical transport models, especially those that address nutrient fluxes.     

IMPACTS OF CALIFORNIA'S CLIMATIC REGIMES AND COASTAL LAND USE CHANGE ON STREAMFLOW CHARACTERISTICS1

ABSTRACT: To investigate the impacts of urbanization and climatic fluctuations on stream flow magnitude and variability in a Mediterranean climate, the HEC‐HMS rainfall/runoff model is used to simulate stream flow for a 14‐year period (October 1, 1988, to September 30, 2002) in the Atascadero Creek watershed located along the southern coast of California for 1929, 1998, and 2050 (estimated) land use conditions (8, 38 and 52 percent urban, respectively). The 14‐year period experienced a range of climatic conditions caused mainly by El Nino‐Southern Oscillation variations. A geographic information system is used to delineate the watershed and parameterize the model, which is calibrated using data from two stream flow and eight rainfall gauges. Urbanization is shown to increase peak discharges and runoff volume while decreasing stream flow variability. In all cases, the annual and 14‐year distributions of stream flow are shown to be highly skewed, with the annual maximum 24 hours of discharge accounting for 22 to 52 percent of the annual runoff and the maximum ten days of discharge from an average El Nino year producing 10 to 15 percent of the total 14‐year discharge. For the entire period of urbanization (1929 to 2050), the average increase in annual maximum discharges and runoff was 45 m³/s (300 percent) and 15 cm (350 percent), respectively. Additionally, the projected increase in urbanization from 1998 to 2050 is half the increase from 1929 to 1998; however, increases in runoff (22 m³/s and 7 cm) are similar for both scenarios because of the region's spatial development pattern.     

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.     

LAKE CHAMPLAIN BASIN NONPOINT SOURCE PHOSPHORUS ASSESSMENT1

ABSTRACT: Existing land use data were used to estimate nonpoint source phosphorus loads to Lake Champlain (Vermont/New York/Quebec) in a loading function model that combined P concentration coefficients with regional hydrologic data. The estimates were verified against monitored loading data, then used to assess the relative magnitudes of contributions from major land uses and regions of the Lake Champlain Basin. The Basin is comprised of 62 percent forest, 28 percent agricultural land, 3 percent urban land, and 7 percent water. The best-fit model estimated an annual total P load of 457 mt/year, which did not differ significantly from the 458 metric tons/year measured for an average hydrologic year, and accurately predicted loads from major tributaries. Agriculture contributes 66 percent of the annual nonpoint source P load to Lake Champlain; urban and forest land contribute 18 percent and 16 percent, respectively. Because agricultural land contributes most nonpoint source P to Lake Champlain, load reduction effort must deal with agricultural sources. However, because the urban 3 percent of the basin contributes 18 percent of the estimated load, high load reduction efficiencies might be achieved by addressing urban sources. This assessment clearly demonstrated the relationship between land use and P loads in the Lake Champlain Basin, a prerequisite for policy-makers to endorse a P management strategy requiring changes in land use and management.     

ANALYSIS AND PREDICTION OF SURFACE RUNOFF IN AN URBANIZING WATERSHED USING SATELLITE IMAGERY1

ABSTRACT: This paper demonstrates how satellite image data [e.g., from Landsat 5 Thematic Mapper (TM)], in conjunction with an urban growth model and simple runoff calculations, can be used to estimate future surface runoff and, by implication, water quality within a watershed. To illustrate the method, predictions of land use change and surface runoff are shown for Spring Creek Watershed, a medium sized urbanizing watershed in Central Pennsylvania. Land cover classifications for this watershed were created from images for summertime 1986 and 1996 and subsequently used as input to the Clarke urban growth model, called SLEUTH, to predict land use changes to the year 2025. Simulations with this model show a progressive growth in the percentage of urban pixels and in impervious surface area in the watershed but also an increase in woodland, primarily in previously clear‐cut areas. Given that woodland area will continue to increase in area, surface runoff into Spring Creek is predicted to remain only slightly above present level. However, should the woodland amount fail to increase, surface runoff is then predicted to increase more significantly during the next 25 years. Finally, the concept of urban sprawl is addressed within the context of predicted increases in urbanization by relating the implied increase in impervious surface area to population density within the watershed.     

MINIMIZING THE IMPACT OF URBANIZATION ON LONG TERM RUNOFF1

Increasing concern about the problems caused by urban sprawl has encouraged development and implementation of smart growth approaches to land use management. One of the goals of smart growth is water resources protection, in particular minimizing the runoff impact of urbanization. To investigate the magnitude of the potential benefits of land use planning for water resources protection, possible runoff impacts of historical and projected urbanization were estimated for two watersheds in Indiana and Michigan using a long term hydrological impact analysis model. An optimization component allowed selection of land use change placements that minimize runoff increase. Optimizing land use change placement would have reduced runoff increase by as much as 4.9 percent from 1973 to 1997 in the Indiana study watershed. For nonsprawl and sprawl scenarios in the Michigan watershed for 1978 to 2040, optimizing land use change placement would have reduced runoff increase by 12.3 percent and 20.5 percent, respectively. The work presented here illustrates both an approach to assessing the magnitude of the impact of smart growth and the significant potential scale of smart growth in moderating runoff changes that result from urbanization. The results of this study have significant implications for urban planning.     

POTENTIAL PHOSPHORUS LOAD REDUCTIONS UNDER THE LAKE OKEECHOBEE REGULATORY PROGRAM1

ABSTRACT: Nutrient loading from beef pastures located within the northern Lake Okeechobee watershed in Florida, has been identified as a source of phosphorus contributing to the accelerated eutrophication of the lake. Since 1989 within the watershed, 557 agricultural drainage sites, mainly beef pasture, have been monitored for compliance under a regulatory program. Of those sites, 154 were actively monitored for phosphorus concentrations from October 1, 1998, to September 30, 1999. Of these 154 sites, 77 were considered to be out of compliance (OOC). An OOC site is defined as having runoff with a 12‐month average phosphorus concentration exceeding the permitted discharge limit. The average annual phosphorous load from the 77 OOC sites for an eight‐year study period from October 1, 1991, to September 30, 1999, was estimated using measured concentration values and simulated runoff obtained from an agricultural nonpoint source pollution model, CREAMS‐WT. The 77 OOC sites produced an estimated average annual 46 metric tonnes of phosphorus load, of which an estimated 22 tonnes of phosphorus reached Lake Okeechobee on an average annual basis. The remaining estimated average annual 24 tonnes of phosphorus load was retained by streams and wetlands in the discharge transport system between the sites and the lake. The estimated average annual load reaching Lake Okeechobee from the OOC sites represented 11 percent of the phosphorus load above a five‐year average annual target load for the lake. However, the OOC site drainage areas represented only 3 percent of the northern watershed that drains into the lake. Of the 77 OOC sites, 12 sites had an average annual phosphorus loading rate equal to or greater than 3.0 kg/ha and were placed on the priority list for the Critical Restoration Project in the Lake Okeechobee watershed. To estimate the possible phosphorus load reductions from the 77 sites, two scenarios were modeled. The first scenario reduced phosphorus concentrations in runoff to the permitted discharge limits under the Lake Okeechobee regulatory program. The second scenario changed current land uses to native rangeland with an estimated annual offsite total phosphorus areal loading rate of 0.114 kg/ha. These two scenarios are hypothetical with assumed concentration values and loading rate. Model results showed that the first management scenario reduced the average annual phosphorus load to the lake by an estimated 15 tonnes. The second scenario reduced the average annual phosphorus load to the lake by an estimated 21 tonnes.     

Weighting Nitrogen and Phosphorus Pixel Pollutant Loads to Represent Runoff and Buffering Likelihoods

Watershed models often estimate annual nitrogen (N) or phosphorus (P) pollutant loads in rural areas with export coefficient (EC) (kg/ha/yr) values based on land cover, and in urban areas as the product of spatially uniform event mean concentration (EMC) (mg/L) values and runoff volume. Actual N and P nonpoint source (NPS) pollutant loading has more spatial complexity due to watershed variation in runoff likelihood and buffering likelihood along surface and subsurface pathways, which can be represented in a contributing area dispersal area (CADA) NPS model. This research develops a CADA NPS model to simulate how watershed properties of elevation, land cover, and soils upslope and downslope of each watershed pixel influence nutrient loading. The model uses both surface and subsurface runoff indices (RI), and surface and subsurface buffer indices (BI), to quantify the runoff and buffering likelihood for each watershed pixel, and generate maps of weighted EC and EMC values that identify NPS pollutant loading hotspots. The research illustrates how CADA NPS model maps and pixel loading values are sensitive to the spatial resolution and accuracy of elevation and land cover data, and model predictions can represent the lower and upper bounds of NPS loading. The model provides managers with a tool to rapidly visualize, rank, and investigate likely areas of high nutrient export.     

ALAWAT: A spatially allocated watershed model for approximating stream,sediment, and pollutant flows in Hawaii,USA

The Ala Wai Canal Watershed Model (ALAWAT) is a planning-level watershed model for approximating direct runoff, streamflow, sediment loads, and loads for up to five pollutants. ALAWAT uses raster GIS data layers including land use, SCS soil hydrologic groups, annual rainfall, and subwatershed delineations as direct model parameter inputs and can use daily total rainfall from up to ten rain gauges and streamflow from up to ten stream gauges. ALAWAT uses a daily time step and can simulate flows for up to ten-year periods and for up to 50 subwatersheds. Pollutant loads are approximated using a user-defined combination of rating curve relationships, mean event concentrations, and loading/washoff parameters for specific subwatersheds, land uses, and times of year. Using ALAWAT, annual average streamflow and baseflow relationships and urban suspended sediment loads were approximated for the Ala Wai Canal watershed (about 10,400 acres) on the island of Oahu, Hawaii. Annual average urban suspended sediments were approximated using two methods: mean event concentrations and pollutant loading and washoff. Parameters for the pollutant loading and washoff method were then modified to simulate the effect of various street sweeping intervals on sediment loads.     

EFFECTS OF PRECIPITATION AND LAND USE ON STORM RUNOFF1

ABSTRACT: Storm-runoff quantity and quality were studied in three watersheds located near St. Paul in Ramsey County, Minnesota, from April 15 through September 15 of 1984, 1985, and 1986 to qualitatively determine the effects of precipitation and selected land uses on storm runoff. In respect to precipitation effects, differences in stormrunoff quantity between years in an urban watershed that lacks wetlands appear to be related to the average storm size (amount of precipitation) during the study period of each year. In contrast, the differences in storm-runoff quantity from watersheds that contain wetlands appear to be related to total precipitation during study period of each year. In respect to land use, the differences in storm-runoff quantity appear to be related to the amounts of impervious and wetland area. The watershed that contains the largest amount of impervious area and smallest amount of wetland area has the largest amount of storm runoff. Differences in storm-runoff quality appear to be related to the amounts of wetland and lake area. The watershed that contains the largest amounts of wetland and lake area has the smallest storm-runoff loading of suspended solids, phosphorus, and nitrogen. The wetland and lake areas likely retain the loading and, subsequently, lower the amount of storm-runoff loading exported from a watershed.     

INTEGRATING PHOSPHORUS AND NITROGEN DECISION MANAGEMENT AT WATERSHED SCALES1

ABSTRACT: The persistence of water quality problems has directed attention towards the reduction of agricultural nonpoint sources of phosphorus (P) and nitrogen (N). We assessed the practical impact of three management scenarios to reduce P and N losses from a mixed land use watershed in central Pennsylvania, USA. Using Scenario 1 (an agronomic soil P threshold of 100 mg Mehlich‐3 P kg^‐1, above which no crop response is expected), 81 percent of our watershed would receive no P as fertilizer or manure. Under Scenario 2 (an environmental soil P threshold of 195 mg Mehlich‐3 P kg^‐1, above which the loss of P in surface runoff and subsurface drainage increases greatly), restricts future P inputs in only 51 percent of the watershed. Finally, using scenario 3 (P and N indices that account for likely source and transport risks), 25 percent of the watershed was at high risk or greater of P loss, while 60 percent of the watershed was classified as of high risk of nitrate (NO₃) leaching. Areas at risk of P loss were near the stream channel, while areas at risk of NO₃ leaching were near the boundaries of the watershed, where freely draining soils and high manure and fertilizer N applications coincide. Remedial measures to minimize P export should focus on critical source areas, while remedial measures to reduce N losses should be source based, concentrating on more efficient use of N by crops.