WATERSHED WEIGHTING OF EXPORT COEFFICIENTS TO MAP CRITICAL PHOSPHOROUS LOADING AREAS1

ABSTRACT: The Export Coefficient model (ECM) is capable of generating reasonable estimates of annual phosphorous loading simply from a watershed's land cover data and export coefficient values (ECVs). In its current form, the ECM assumes that ECVs are homogeneous within each land cover type, yet basic nutrient runoff and hydrological theory suggests that runoff rates have spatial patterns controlled by loading and filtering along the flow paths from the upslope contributing area and downslope dispersal area. Using a geographic information system (GIS) raster, or pixel, modeling format, these contributing area and dispersal area (CADA) controls were derived from the perspective of each individual watershed pixel to weight the otherwise homogeneous ECVs for phosphorous. Although the CADA‐ECM predicts export coefficient spatial variation for a single land use type, the lumped basin load is unaffected by weighting. After CADA weighting, a map of the new ECVs addressed the three fundamental criteria for targeting critical pollutant loading areas: (1) the presence of the pollutant, (2) the likelihood for runoff to carry the pollutant offsite, and (3) the likelihood that buffers will trap nutrients prior to their runoff into the receiving water body. These spatially distributed maps of the most important pollutant management areas were used within New York's West Branch Delaware River watershed to demonstrate how the CADA‐ECM could be applied in targeting phosphorous critical loading areas.     

PARTIAL AREA RESPONSE ON SMALL SEMIARID WATERSHEDS1

ABSTRACT: Significant errors in estimating surface runoff and erosion rates are possible if a watershed is assumed to contribute runoff uniformly over the entire area, when actually only a portion of the entire area may be contributing. Generation of overland flow on portions of small semiarid watersheds was analyzed by three methods: an average loss rate procedure, a lumped-linear model, and a distributed-nonlinear model. These methods suggested that, on the average, 45, 60, and 50% of the drainage area was contributing runoff at the watershed outlet. Infiltrometer data support the partial area concept and indicate that the low infiltration zones are the runoff source areas as simulated with the distributed-nonlinear model.     

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.     

A DISTRIBUTED DYNAMIC WATERSHED MODEL1

ABSTRACT: A distributed watershed model was developed to mathematically simulate overland and channel flow for a single-event storm. The modeled watersheds in the study were subdivided into rectangular grid elements. All hydrologically significant parameters, such as land slope, rainfall and precipitation excess, were assumed to be uniform within each element. The Green-Ampt method was adopted to generate precipitation excess for each element during the simulation period. A two-dimensional diffusion wave model was used for overland flow routing and an iterative Alternative Direction Implicit scheme was used to solve the simultaneous overland flow equations. Once the overland flow became inflow to the channel, a one-dimensional dynamic wave flood routing technique, based on a four-point, implicit, non-linear finite difference solution of the St. Venant equation of unsteady flow, was applied. A limited number of comparisons were made between simulated and observed hydrographs for areas of about one square mile. Given the appropriate parameters, the model was able to accurately simulate runoff for single-event storms. This paper describes a distributed watershed model developed to simulate overland and channel flow. Comparisons were made between simulated and observed hydrographs for three watersheds. The model was able to accurately simulate the runoff for single-event storms using 61-m by 61-m (200 ft by 200 ft) watershed grid elements.     

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.     

Phosphorus transport in overland flow in response to position of manure application

Phosphorus (P) loss in overland flow varies with spatial distribution of soil P, management, and hydrological pathways. The effect of flow time, flowpath length, and manure position on P loss in overland flow from two central Pennsylvania soils packed in boxes of varying length (0.5, 1.0, 1.5, 2.75, and 4.0 m long x 15 cm wide x 5 cm deep) were examined by collecting flow samples at 5-min intervals for 30 min (50 mm h(-1) rainfall) without and with 75 kg P ha(-1) applied as swine (Sus scrofa) manure over 0.5 m of the box slope length at distances of 0 to 3.5 m from the downslope collection point. Dissolved reactive P concentration was more closely related to the proportion of clay in sediment of overland flow before (r = 0.98) than after (r = 0.56) manure application. This was attributed to the transport of larger, low-density particles after applying manure. The concentration of dissolved and particulate P fractions decreased with increasing flowpath length, due to dilution rather than sorption of P by surface soil during overland flow. Total P loss (mainly as particulate P) from the Watson channery silt loam (fine-loamy, mixed, active, mesic Typic Fragiudult) was more than from Berks channery silt loam (loamy-skeletal, mixed, active, mesic Typic Dystrudept), even with manure applied. Thus, while P loss in overland flow is affected by where manure is applied relative to flowpath length, initial soil P concentration should not be discounted when looking at areas of potential P loss within a watershed.     

Wetland Flowpaths Mediate Nitrogen and Phosphorus Concentrations across the Upper Mississippi River Basin

Eutrophication, harmful algal blooms, and human health impacts are critical environmental challenges resulting from excess nitrogen and phosphorus in surface waters. Yet we have limited information regarding how wetland characteristics mediate water quality across watershed scales. We developed a large, novel set of spatial variables characterizing hydrological flowpaths from wetlands to streams, that is, “wetland hydrological transport variables,” to explore how wetlands statistically explain the variability in total nitrogen (TN) and total phosphorus (TP) concentrations across the Upper Mississippi River Basin (UMRB) in the United States. We found that wetland flowpath variables improved landscape-to-aquatic nutrient multilinear regression models (from R² = 0.89 to 0.91 for TN; R² = 0.53 to 0.84 for TP) and provided insights into potential processes governing how wetlands influence watershed-scale TN and TP concentrations. Specifically, flowpath variables describing flow-attenuating environments, for example, subsurface transport compared to overland flowpaths, were related to lower TN and TP concentrations. Frequent hydrological connections from wetlands to streams were also linked to low TP concentrations, which likely suggests a nutrient source limitation in some areas of the UMRB. Consideration of wetland flowpaths could inform management and conservation activities designed to reduce nutrient export to downstream waters.     

Soil, water and nutrient conservation in mountain farming systems: case-study from the Sikkim Himalaya

The Khanikhola watershed in Sikkim is agrarian with about 50% area under rain-fed agriculture representing the conditions of the middle mountains all over the Himalaya. The study was conducted to assess overland flow, soil loss and subsequent nutrient losses from different land uses in the watershed, and identify biotechnological inputs for management of mountain farming systems. Overland flow, soil and nutrient losses were very high from open agricultural (cropped) fields compared to other land uses, and more than 72% of nutrient losses were attributable to agriculture land use. Forests and large cardamom agroforestry conserved more soil compared to other land uses. Interventions, like cultivation of broom grass upon terrace risers, N2-fixing Albizia trees for maintenance of soil fertility and plantation of horticulture trees, have reduced the soil loss (by 22%). Soil and water conservation values (> 80%) of both large cardamom and broom grass were higher compared to other crops. Use of N2-fixing Albizia tree in large cardamom agroforestry and croplands contributed to soil fertility, and increased productivity and yield. Bio-composting of farm resources ensured increase in nutrient availability specially phosphorus in cropped areas. Agricultural practices in mountain areas should be strengthened with more agroforestry components, and cash crops like large cardamom and broom grass in agroforestry provide high economic return and are hydroecologically sustainable.     

Phosphorus transport pathways to streams in tile-drained agricultural watersheds

Agriculture is a major nonpoint source of phosphorus (P) in the Midwest, but how surface runoff and tile drainage interact to affect temporal concentrations and fluxes of both dissolved and particulate P remains unclear. Our objective was to determine the dominant form of P in streams (dissolved or particulate) and identify the mode of transport of this P from fields to streams in tile-drained agricultural watersheds. We measured dissolved reactive P (DRP) and total P (TP) concentrations and loads in stream and tile water in the upper reaches of three watersheds in east-central Illinois (Embarras River, Lake Fork of the Kaskaskia River, and Big Ditch of the Sangamon River). For all 16 water year by watershed combinations examined, annual flow-weighted mean TP concentrations were >0.1 mg L(-1), and seven water year by watershed combinations exceeded 0.2 mg L(-1). Concentrations of DRP and particulate P (PP) increased with stream discharge; however, particulate P was the dominant form during overland runoff events, which greatly affected annual TP loads. Concentrations of DRP and PP in tiles increased with discharge, indicating tiles were a source of P to streams. Across watersheds, the greatest DRP concentrations (as high as 1.25 mg L(-1)) were associated with a precipitation event that followed widespread application of P fertilizer on frozen soils. Although eliminating this practice would reduce the potential for overland runoff of P, soil erosion and tile drainage would continue to be important transport pathways of P to streams in east-central Illinois.     

RASTER-BASED HYDROLOGIC MODELING OF SPATIALLY-VARIED SURFACE RUNOFF1

ABSTRACT: The proliferation of watershed databases in raster Geographic Information System (GIS) format and the availability of radar-estimated rainfall data foster rapid developments in raster-based surface runoff simulations. The two-dimensional physically-based rainfall-runoff model CASC2D simulates spatially-varied surface runoff while fully utilizing raster GIS and radar-rainfall data. The model uses the Green and Ampt infiltration method, and the diffusive wave formulation for overland and channel flow routing enables overbank flow storage and routing. CASC2D offers unique color capabilities to display the spatio-temporal variability of rainfall, cumulative infiltrated depth, and surface water depth as thunderstorms unfold. The model has been calibrated and independently verified to provide accurate simulations of catchment response to moving rainstorms on watersheds with spatially-varied infiltration. The model can accurately simulate surface runoff from flashfloods caused by intense thunderstorms moving across partial areas of a watershed.     

Uptake and release of phosphorus from overland flow in a stream environment

Phosphorus runoff from agricultural fields has been linked to fresh-water eutrophication. However, edge-of-field P losses can be modified by benthic sediments during stream flow by physiochemical processes associated with Al, Fe, and Ca, and by biological assimilation. We investigated fluvial P when exposed to stream-bed sediments (top 3 cm) collected from seven sites representing forested and agricultural areas (pasture and cultivated), in a mixed-land-use watershed. Sediment was placed in a 10-m-long, 0.2-m-wide fluvarium to a 3-cm depth and water was recirculated over the sediment at 2 L s(-1) and 5% slope. When overland flow (4 mg dissolved reactive phosphorus [DRP] and 9 mg total phosphorus [TP] L(-1)) from manured soils was first recirculated, P uptake was associated with Al and Fe hydrous oxides for sediments from forested areas (pH 5.2-5.4) and by Ca for sediments from agricultural areas (pH 6.5-7.2). A large increase (up to 200%) in readily available P NH4Cl fraction was noted. After 24 h, DRP concentration in channel flow was related to sediment solution P concentration at which no net sorption or desorption of P occurs (EPC0) (r2 = 0.77), indicating quasi-equilibrium. When fresh water (approximately 0.005 mg P L(-1) mean base flow DRP at seven sites) was recirculated over the sediments for 24 h, P release kinetics followed an exponential function. Microbial biomass P accounted for 34 to 43% of sediment P uptake from manure-rich overland flow. Although abiotic sediment processes played a dominant role in determining P uptake, biotic process are clearly important and both should be considered along with the location and management of landscape inputs for remedial strategies to be effective.     

DELIVERY OF SUSPENDED SEDIMENT AND POLLUTANTS FROM NONPOINT SOURCES DURING OVERLAND FLOW1

ABSTRACT: Delivery of sediment and particulate pollutants from diffuse sources is shown to be related to the loss of sediment carrying energy of runoff during the overland flow phase. The loss is caused by the termination of rainfall and by reduction of flow energy during the recession phase of the overland flow hydrograph. It has been demonstrated both by theoretical analyses and experimental measurements that the saturated sediment concentration in overland flow is a function of rainfall erosivity and the runoff flow rate. The hypotheses were verified by field measurements from a small homogeneous watershed.     

INTEGRATING LANDSCAPE ASSESSMENT AND HYDROLOGIC MODELING FOR LAND COVER CHANGE ANALYSIS1

ABSTRACT: Significant land cover changes have occurred in the watersheds that contribute runoff to the upper San Pedro River in Sonora, Mexico, and southeast Arizona. These changes, observed using a series of remotely sensed images taken in the 1970s, 1980s, and 1990s, have been implicated in the alteration of the basin hydrologic response. The Cannonsville subwatershed, located in the Catskill/Delaware watershed complex that delivers water to New York City, provides a contrast in land cover change. In this region, the Cannonsville watershed condition has improved over a comparable time period. A landscape assessment tool using a geographic information system (GIS) has been developed that automates the parameterization of the Soil and Water Assessment Tool (SWAT) and KINEmatic Runoff and EROSion (KINEROS) hydrologic models. The Automated Geospatial Watershed Assessment (AGWA) tool was used to prepare parameter input files for the Upper San Pedro Basin, a subwatershed within the San Pedro undergoing significant changes, and the Cannonsville watershed using historical land cover data. Runoff and sediment yield were simulated using these models. In the Cannonsville watershed, land cover change had a beneficial impact on modeled watershed response due to the transition from agriculture to forest land cover. Simulation results for the San Pedro indicate that increasing urban and agricultural areas and the simultaneous invasion of woody plants and decline of grasslands resulted in increased annual and event runoff volumes, flashier flood response, and decreased water quality due to sediment loading. These results demonstrate the usefulness of integrating remote sensing and distributed hydrologic models through the use of GIS for assessing watershed condition and the relative impacts of land cover transitions on hydrologic response.     

INTEGRATION OF GEOGRAPHIC INFORMATION SYSTEMS AND A COMPUTER MODEL TO EVALUATE IMPACTS OF AGRICULTURAL RUNOFF ON WATER QUALITY1

ABSTRACT: This study integrates an Agricultural Non-Point Source Pollution Model (AGNPS), the Geographic Resource Analysis Support System (GRASS) (U.S. Army Corps of Engineers, 1987), and GRASS WATERWORKS (a hydrologic modeling tool box being developed at the Michigan State University Center for Remote Sensing) to evaluate the impact of agricultural runoff on water quality in the Cass River, a subwatershed of Saginaw Bay. AGNPS is used to estimate the amounts, origin, and distribution of sediment, nitrogen (N), and phosphorus (P) in the watershed. GRASS and GRASS WATERWORKS are used to generate parameters needed for AGNPS from digital maps, which include soil association, land use, watershed boundaries, water features, and digital elevation. Outputs of the model include spatially distributed estimates of volume and peak runoff, overland and channel erosion, sediment yields, and concentrations of nitrogen and phosphorus. Management scenarios are explored in the AGNPS model to minimize sedimentation and nutrient loading. Scenarios evaluated include variations in crop cover, tillage methods, and other agricultural management practices. In addition, areas vulnerable to erosion are identified for best management practices.     

SIMULATION OF SEDIMENT AND PLANT NUTRIENT LOSSES BY THE CREAMS WATER QUALITY MODEL1

ABSTRACT: CREAMS was applied to a field-sized watershed planted to cotton in the Limestone Valley region of northern Alabama. The field was cultivated for three years with conventional tillage (CvT) followed by three years of conservation tillage (CsT). CREAMS is composed of three components: hydrology, erosion, and chemistry. Surface runoff and losses of sediment, N and P were simulated and results were compared with the observed data from the watershed. Curve numbers recommended in the CREAMS user's guide were not adequate for the watershed conditions. The hydrology submodel improved runoff simulation from CvT and CsT when field-data based curve numbers were used. The erosion submodel demonstrated that CsT reduced sediment loss more than CvT, even though CsT had higher runoff than CvT. The nutrient submodel based on the simulated runoff and sediment underpredicted N loss for both CvT and CsT. This submodel, however, accurately predicted P loss for CvT, but underpredicted for CsT (50 percent lower than the observed). The results of CREAMS simulation generally matched the observed order of magnitude for higher runoff, lower sediment, and higher N and P losses from CsT than from CvT.     

Managing stormwater for urban sustainability: an evaluation of local comprehensive plans in the Chesapeake Bay watershed region

This study uses a developed plan coding protocol in evaluating the quality of 76 comprehensive plans to examine whether local comprehensive plans have adequately integrated the concepts of sustainable stormwater management. The Chesapeake Bay watershed was chosen for the investigation because degraded stormwater runoff from nearby urban and suburban jurisdictions have critically polluted the watershed. The findings indicate that the majority of local governments have not sufficiently incorporated the sustainable stormwater management principles into their comprehensive plans. Five plan components (factual basis, goals and objectives, inter-organizational cooperation, policies, tools and strategies, and implementation) appear weak in realizing the concepts. The current study concludes by providing policy implications and recommendations to increase awareness and understanding of sustainable stormwater management concepts and to produce better implementation plans that integrate stormwater, ecosystem, and environmental planning comprehensively.     

Integrating contributing areas and indexing phosphorus loss from agricultural watersheds

Most states in the USA have adopted P Indexing to guide P-based management of agricultural fields by identifying the relative risk of P loss at farm and watershed scales. To a large extent, this risk is based on hydrologic principles that frequently occurring storms can initiate surface runoff from fields. Once initiated, this hydrological pathway has a high potential to transport P to the stream. In regions where hydrologically active areas of watersheds vary in time and space, surface runoff generation by "saturation excess" has been linked to distance from stream, with larger events resulting in larger contributing distances. Thus, storm-return period and P loss from a 39.5-ha mixed-land-use watershed in Pennsylvania was evaluated to relate return-period thresholds and distances contributing P to streams. Of 248 storm flows between 1997 and 2006, 93% had a return period of 1 yr, contributing 47% of total P (TP) export, while the largest two storms (10-yr return period) accounted for 23% of TP export. Contributing distance thresholds for the watershed were determined (50-150 m) for a range of storm-return periods (1-10 yr) from hydrograph analysis. By modifying storm-return period thresholds in the P Index and thereby contributing distance, it is possible to account for greater risk of P loss during large storms. For instance, increasing return period threshold from 1 (current P indices) to 5 yr, which accounted for 67% of TP export, increased the P-management restricted area from 20 to 58% of the watershed. An increase in impacted area relative to a decreased risk of P loss creates a management-policy dilemma that cannot be ignored.     

GENERATING DESIGN HYDROGRAPHS BY DEM ASSISTED GEOMORPHIC RUNOFF SIMULATION: A CASE STUDY1

ABSTRACT: The geomorphic instantaneous unit hydrograph (GIUH) may be one of the most successful methodologies for predicting flow characteristics in ungauged watersheds. However, one difficulty in applying the GIUH model is determination of travel time, and the other difficulty is the large amount of geomorphologic information required in the study watershed. Recently, using the kinematic-wave theory Lee and Yen (1997) have analytically determined the travel times for overland and channel flows in watersheds. The limitation of using an empirical velocity equation to estimate the runoff travel time for a specified watershed is then relaxed. To simplify the time-consuming work involved in geomorphic parameter measurement on topographic maps, the GIUH model is linked with geographic information systems to obtain geomorphic parameters from digital elevation models. In this paper, a case study performed for peak flow analysis in an ungauged watershed is presented. The geomorphic characteristics of the study watershed were analyzed using a digital elevation model and were used to construct the runoff simulation model. The design storm was then applied to the geomorphic runoff simulation model to obtain the design hydrograph. The analytical procedures proposed in this study can provide a convenient way for hydrologists to estimate hydrograph characteristics based on limited hydrologic information.     

RUNOFF RESPONSE FROM TWO RECLAIMED WATERSHEDS1

ABSTRACT: Detailed studies of the surface hydrology of reclaimed surface-mined watersheds for both rainfall and snowmelt events are non-existent for central Alberta yet this information is crucial for design of runoff conveyance and storage structures. A study was initiated in 1992 with principal objectives of quantifying surface runoff for both summer rainfall and spring snowmelt events and identifying the dominant flow processes occurring in two reclaimed watersheds. Snowmelt accounted for 86 and 100% of annual watershed runoff in 1993 and 1994, respectively. The highest instantaneous peak flow was recorded during a summer rainfall event with a return period of greater than 50 years. Infiltration-excess overland flow was identified as the dominant flow process occurring within the Sandy Subsoil Watershed, whereas saturation overland flow was the principal runoff process occurring within the West Watershed.     

Phosphorus concentrations in overland flow from diverse locations on a New York dairy farm

The National Phosphorus Project rainfall simulator was used to quantify overland flow and P transport from nine sites distributed throughout the watershed of a New York City Watershed Agriculture Program collaborating dairy farm. Observed concentrations of total dissolved phosphorus (TDP) were low (0.007-0.12 mg L(-1)) in flow from deciduous forest, extensively managed pasture, and hillside seeps; moderate (0.18-0.64 mg L(-1)) in flow from intensively managed pastures, a hayfield, and a cow path; and extremely high (11.6 mg L(-1)) in flow from a manured barnyard. Concentrations of TDP from sites without fresh manure were strongly correlated with soil test P (TDP [mg L(-1)] = 0.0056 + 0.0180 x Morgan's soil test phosphorus [STP, mg kg(-1)]; R2 = 84%). Observed concentrations of suspended solids were low (16-137 mg L(-1)) in flow from vegetated sites, but were higher (375-615 mg L(-1)) in flow from sites with little ground cover (barnyard, cow path, plowed field). Under dry summer conditions the time to observed overland flow was shorter (<18 min) for nonfield areas (seeps, barnyard, cow path) than for field and forest areas (27-93 min), indicating that hydrologically active nonfield areas of minor spatial extent but with high soil P (e.g., cow paths and barnyards) can play a significant role in summertime P loading. When soils started from field capacity (second-day) time to overland flow was uniformly less than 23 min, indicating that under wet watershed conditions low-P source areas can dilute overland flow from concentrated sources.