UPLAND SOIL EROSION SIMULATION FOR AGRICULTURAL WATERSHEDS1

ABSTRACT: A soil erosion simulation model that considered the physical conditions of agricultural watersheds and that interfaced with the modified USDAHL-74 watershed hydrology model was developed. The erosion model simulates the detachment and transport of soil particles caused by raindrop impact and overland flow from rill and interrill areas. The model considers temporal and spatial variation of plant residue, crop canopy cover, snow cover, and the moisture content of surface soil as modifying factors of the erosive forces of raindrop impact and overland flow. The hydrology model simulates overland flow and some of the physical parameters that are used in the erosion model. The simulation is executed in the time interval determined by the rainfall rate or snowmelt rate. The erosion model compares the transport capacity of the overland flow and the sediment loaded in the overland flow to determine the fate account for the free soil particles that have already been detached and are readily available to be transported by the overland flow. The model was tested with data from two small agricultural watersheds in the Palouse region of the Pacific Northwest dryland. The model was calibrated by trial-and-error to determine the coefficients of the model.     

SIMULATING THREE-DIMENSIONAL GROUND WATER RESPONSE IN A SMALL MOUNTAINOUS WATERSHED1

ABSTRACT: Snowmelt from deep mountainous snowpacks is seldom rapid enough to exceed infiltration rates; thus, the source of streamflow in many mountainous watersheds is snowmelt recharge through shallow ground water systems. The hydrologic response and interaction between surface and sub-surface flow processes in these watersheds, which is controlled by basin structure, the spatial distribution of snowmelt, and the hydrogeology of the subsurface, are not well understood. The purpose of this study was to test a three-dimensional ground water model using simulated snowmelt input to simulate ground water response to spatially distributed snowmelt on the Upper Sheep Creek Watershed located within the Reynolds Creek Experimental Watershed in Southwestern Idaho. The model was used to characterize the mountainous aquifer and to delineate the subsurface flow mechanisms. Difficulty in finding a reasonable combination of grid spacing and time stepping within the model was encountered due to convergence problems with the Picard solution to the non-linear variably saturated ground water flow equations. Simulation results indicated that flow may be either unconfined or confined depending on inflow rate and hydrogeologic conditions in the watershed. The flow mechanism had a much faster response time when confined flow occurred. Response to snowmelt from a snow drift approximately 90 m away took only a few hours when flow was confined. Simulated results showed good agreement with piezometer measurements both in magnitude and timing; however, convergence problems with the Picard solution limited applicability of the model.     

Variable Source Area Hydrology Modeling with the Water Erosion Prediction Project Model

In nondegraded watersheds of humid climates, subsurface flow patterns determine where the soil saturates and where surface runoff is occurring. Most models necessarily use infiltration‐excess (i.e., Hortonian) runoff for predicting runoff and associated constituents because subsurface flow algorithms are not included in the model. In this article, we modify the Water Erosion Prediction Project (WEPP) model to simulate subsurface flow correctly and to predict the spatial and temporal location of saturation, the associated lateral flow and surface runoff, and the location where the water can re‐infiltrate. The modified model, called WEPP‐UI, correctly simulated the hillslope drainage data from the Coweeta Hydrologic Laboratory hillslope plot. We applied WEPP‐UI to convex, concave, and S‐shaped hillslope profiles, and found that multiple overland flow elements are needed to simulate distributed lateral flow and runoff well. Concave slopes had the greatest runoff, while convex slopes had the least. Our findings concur with observations in watersheds with saturation‐excess overland flow that most surface runoff is generated on lower concave slopes, whereas on convex slopes runoff infiltrates before reaching the stream. Since the WEPP model is capable of simulating both saturation‐excess and infiltration‐excess runoff, we expect that this model will be a powerful tool in the future for managing water quality.     

SAGEBRUSH RANGELAND HYDROLOGY AND EVALUATION OF THE SPUR HYDROLOGY MODEL1

ABSTRACT: An excellent hydrologic record on sagebrush range-land has been developed at the Reynolds Creek Experimental Watershed in southwestern Idaho. The objectives of this paper were two-fold: (1) to analyze and describe the hydrologic record (8–18 years) from four sagebrush watersheds (1–83 ha); and (2) to evaluate the hydrology component of SPUR, a comprehensive rangeland model. The watersheds represent a gradient in elevation (1180–1658 m) and precipitation (240–350 mm/yr). Runoff was a small fraction (> 2 percent) of the total water budget for all of the watersheds. It occurred very infrequently at the three lower elevation watersheds: Summit, Flats, and Nancy Gulch. At Lower Sheep, the highest elevation watershed, runoff occurred most years for a period of 1 to 17 weeks in the winter. Frozen soil combined with rainfall or snowmelt was associated with most of the runoff from Flats and Nancy Gulch. At Summit summertime thunderstorms produced all of the runoff. The average annual sediment yield from all of the watersheds was low (17–950 kg/ha). It was highest from Summit, which had well developed alluvial channels and very steep slopes. SPUR was able to simulate runoff with reasonable accuracy only at Summit, where frozen soils were not a factor. There was poor correlation between predicted and actual annual 8ediment loss. The model tended to overpredict evapotranspiration early in the growing season and underpredict it in the late summer.     

Curve Number Derivation for Watersheds Draining Two Headwater Streams in Lower Coastal Plain South Carolina,USA

The objective of this study was to assess curve number (CN) values derived for two forested headwater catchments in the Lower Coastal Plain (LCP) of South Carolina using a three‐year period of storm event rainfall and runoff data in comparison with results obtained from CN method calculations. Derived CNs from rainfall/runoff pairs ranged from 46 to 90 for the Upper Debidue Creek (UDC) watershed and from 42 to 89 for the Watershed 80 (WS80). However, runoff generation from storm events was strongly related to water table elevation, where seasonally variable evapotranspirative wet and dry moisture conditions persist. Seasonal water table fluctuation is independent of, but can be compounded by, wet conditions that occur as a result of prior storm events, further complicating flow prediction. Runoff predictions for LCP first‐order watersheds do not compare closely to measured flow under the average moisture condition normally associated with the CN method. In this study, however, results show improvement in flow predictions using CNs adjusted for antecedent runoff conditions and based on water table position. These results indicate that adaptations of CN model parameters are required for reliable flow predictions for these LCP catchments with shallow water tables. Low gradient topography and shallow water table characteristics of LCP watersheds allow for unique hydrologic conditions that must be assessed and managed differently than higher gradient watersheds.     

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.     

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.     

Hydrologic Calibration and Validation of SWAT in a Snow‐Dominated Rocky Mountain Watershed,Montana, U.S.A.1

Abstract: The Soil and Water Assessment Tool (SWAT) has been applied successfully in temperate environments but little is known about its performance in the snow‐dominated, forested, mountainous watersheds that provide much of the water supply in western North America. To address this knowledge gap, we configured SWAT to simulate the streamflow of Tenderfoot Creek (TCSWAT). Located in central Montana, TCSWAT represents a high‐elevation watershed with ～85% coniferous forest cover where more than 70% of the annual precipitation falls as snow, and runoff comes primarily from spring snowmelt. Model calibration using four years of measured daily streamflow, temperature, and precipitation data resulted in a relative error (RE) of 2% for annual water yield estimates, and mean paired deviations (Dv) of 36 and 31% and Nash‐Sutcliffe (NS) efficiencies of 0.90 and 0.86 for monthly and daily streamflow, respectively. Model validation was conducted using an additional four years of data and the performance was similar to the calibration period, with RE of 4% for annual water yields, Dv of 43% and 32%, and NS efficiencies of 0.90 and 0.76 for monthly and daily streamflow, respectively. An objective, regression‐based model invalidation procedure also indicated that the model was validated for the overall simulation period. Seasonally, SWAT performed well during the spring and early summer snowmelt runoff period, but was a poor predictor of late summer and winter base flow. The calibrated model was most sensitive to snowmelt parameters, followed in decreasing order of influence by the surface runoff lag, ground water, soil, and SCS Curve Number parameter sets. Model sensitivity to the surface runoff lag parameter reflected the influence of frozen soils on runoff processes. Results indicated that SWAT can provide reasonable predictions of annual, monthly, and daily streamflow from forested montane watersheds, but further model refinements could improve representation of snowmelt runoff processes and performance during the base flow period in this environment.     

Effect of varying the phosphorus content of dairy cow diets on losses of phosphorus in overland flow following surface applications of manure

The increasing use of concentrate feedstuffs within Northern Ireland dairy systems has resulted in significant farm gate phosphorus (P) surpluses, and these have contributed to increased soil P levels and risk of P loss to overland flow. However, the P content of feed concentrates can be lowered without compromising animal performance. This study focuses on P losses from grassland and evaluates how adjusting the P content of manure impacts on the P composition and concentration in overland flow. Dairy cows were offered diets containing 5.3 to 3.0 g P kg(-1) dry matter (DM) and produced manures with a range of P contents. Manure was applied at a rate of 50 m3 ha(-1) to 0.5-m2 grassland plots, and simulated rainfall (40 mm h(-1)) was applied repeatedly 2, 9, 28, and 49 d after during the summer, winter, and spring. Decreasing the P content in the diet, from the highest to the lowest P treatment (43%), produced a proportionately greater reduction in manure TP content (63%), but reductions were not exclusively in the water-soluble fraction. Following surface applications of manure, P concentrations in overland flow increased in all seasons (P < or = 0.001), while the greatest impact of varying the manure P content was most evident during the first simulated overland flow event. When diet P content was reduced from 5.4 to 3.0 g P kg(-1) DM, a statistically significant reduction in runoff P concentration was observed in all seasons. Elevated P concentrations in overland flow were observed for 28 d in spring and 9 d in summer and winter. The large drop in P concentrations between simulated rainfall events on Day 2 and Day 9 suggests that increasing the time interval between manure application and the generation of overland flow has a greater impact on P losses than does varying the dietary P content.     

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.     

THE EFFECTS OF CLIMATE CHANGE ON STREAM FLOW AND NUTRIENT LOADING'

ABSTRACT: This study assesses the potential impact of climate change on stream flow and nutrient loading in six watersheds of the Susquehanna River Basin using the Generalized Watershed Loading Function (GWLF). The model was used to simulate changes in stream flow and nutrient loads under a transient climate change scenario for each watershed. Under an assumption of no change in land cover and land management, the model was used to predict monthly changes in stream flow and nutrient loads for future climate conditions. Mean annual stream flow and nutrient loads increased for most watersheds, but decreased in one watershed that was intensively cultivated. Nutrient loading slightly decreased in April and late summer for several watersheds as a result of early snowmelt and increasing evapotranspiration. Spatial and temporal variability of stream flow and nutrient loads under the transient climate scenario indicates that different approaches for future water resource management may be useful.     

Investigation of the Curve Number Method For Surface Runoff Estimation In Tropical Regions

This study tests the applicability of the curve number (CN) method within the Soil and Water Assessment Tool (SWAT) to estimate surface runoff at the watershed scale in tropical regions. To do this, surface runoff simulated using the CN method was compared with observed runoff in numerous rainfall‐runoff events in three small tropical watersheds located in the Upper Blue Nile basin, Ethiopia. The CN method generally performed well in simulating surface runoff in the studied watersheds (Nash‐Sutcliff efficiency [NSE] > 0.7; percent bias [PBIAS] < 32%). Moreover, there was no difference in the performance of the CN method in simulating surface runoff under low and high antecedent rainfall (PBIAS for both antecedent conditions: ~30%; modified NSE: ~0.4). It was also found that the method accurately estimated surface runoff at high rainfall intensity (e.g., PBIAS < 15%); however, at low rainfall intensity, the CN method repeatedly underestimated surface runoff (e.g., PBIAS > 60%). This was possibly due to low infiltrability and valley bottom saturated areas typical of many tropical soils, indicating that there is scope for further improvements in the parameterization/representation of tropical soils in the CN method for runoff estimation, to capture low rainfall‐intensity events. In this study the retention parameter was linked to the soil moisture content, which seems to be an appropriate approach to account for antecedent wetness conditions in the tropics.     

RUNOFF AND SEDIMENT YIELD UNDER GRAZING IN FOOTHILLS FESCUE GRASSLANDS OF ALBERTA1

ABSTRACT: Grazing can have a profound impact on infiltration and thus runoff and erosion. The objectives of this study were to quantify the effects of select grazing systems on rainfall and snowmelt induced runoff and sediment yield from sloped areas of the foothills fescue grasslands of Alberta, Canada. The effects of two grazing intensities (heavy and very heavy) for two durations (short duration and continuous throughout the growing season) were compared to an ungrazed control between June 1988 and April 1991. Runoff was measured using 1-rn² runoff frames and collection bucket systems. Sediment yields were then determined on samples from the collected runoff. Snowmelt was the dominant source of runoff. Snowmelt runoff was higher from the heavily grazed areas than from the very heavily grazed areas, due to the higher standing vegetation which accumulated snow in the former areas. Sediment yields as a result of snowmelt were generally low in all areas. Only a few summer storms caused runoff. Runoff volumes and sediment yields from summer rainstorms were low, due to low rainfall and to generally dry antecedent soil moisture conditions. The greatest risk of summer runoff, and thus sediment yield, appears to occur in August.     

A Hydrologic Data Screening Procedure for Exploring Monotonic Trends and Shifts in Rainfall and Runoff Patterns

A thorough understanding of past and present hydrologic responses to changes in precipitation patterns is crucial for predicting future conditions. The main objectives of this study were to determine temporal changes in rainfall‐runoff relationship and to identify significant trends and abrupt shifts in rainfall and runoff time series. Ninety‐year rainfall and runoff time series datasets from the Gasconade and Meramec watersheds in east‐central Missouri were used to develop data screening procedure to assess changes in the rainfall and runoff temporal patterns. A statistically significant change in mean and variance was detected in 1980 in the rainfall and runoff time series within both watersheds. In addition, both the rainfall and runoff time series indicated the presence of nonstationary attributes such as statistically significant monotonic trends and/or change in mean and variance, which should be taken into consideration when using the time series to predict future scenarios. The annual peak runoff and the annual low flow in the Meramec watershed showed significant temporal changes compared to that in the Gasconade watershed. Water loss in both watersheds was found to be significantly increasing which is potentially due to the increase in groundwater pumping for water supply purposes.     

Rate of fall-applied liquid swine manure: effects on runoff transport of sediment and phosphorus

Reducing the delivery of phosphorus (P) from land-applied manure to surface water is a priority in many watersheds. Manure application rate can be controlled to manage the risk of water quality degradation. The objective of this study was to evaluate how application rate of liquid swine manure affects the transport of sediment and P in runoff. Liquid swine manure was land-applied and incorporated annually in the fall to runoff plots near Morris, Minnesota. Manure application rates were 0, 0.5, 1, and 2 times the rate recommended to supply P for a corn (Zea mays L.)-soybean [Glycine max (L.) Merr.] rotation. Runoff volume, sediment, and P transport from snowmelt and rainfall were monitored for 3 yr. When manure was applied at the highest rate, runoff volume and sediment loss were less than the control plots without manure. Reductions in runoff volume and soil loss were not observed for spring runoff when frozen soil conditions controlled infiltration rates. The reduced runoff and sediment loss from manure amended soils compensated for addition of P, resulting in similar runoff losses of total P among manure application rates. However, losses of dissolved P increased with increasing manure application rate for runoff during the spring thaw period. Evaluation of water quality risks from fall-applied manure should contrast the potential P losses in snowmelt runoff with the potential that incorporated manure may reduce runoff and soil loss during the summer.     

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.     

Projections of 21st Century Sierra Nevada Local Hydrologic Flow Components Using an Ensemble of General Circulation Models1

Abstract: Sierra Nevada snowmelt and runoff is a key source of water for many of California’s 38 million residents and nearly the entire population of western Nevada. The purpose of this study was to assess the impacts of expected 21st Century climatic changes in the Sierra Nevada at the subwatershed scale, for all hydrologic flow components, and for a suite of 16 General Circulation Models (GCMs) with two emission scenarios. The Soil and Water Assessment Tool (SWAT) was calibrated and validated at 35 unimpaired streamflow sites. Results show that temperatures are projected to increase throughout the Sierra Nevada, whereas precipitation projections vary between GCMs. These climatic changes drive a decrease in average annual streamflow and an advance of snowmelt and runoff by several weeks. The largest streamflow reductions were found in the mid‐range elevations due to less snow accumulation, whereas the higher elevation watersheds were more resilient due to colder temperatures. Simulation results showed that decreases in snowmelt affects not only streamflow, but evapotranspiration, surface, and subsurface flows, such that less water is available in spring and summer, thus potentially affecting aquatic and terrestrial ecosystems. Declining spring and summer flows did not equally affect all subwatersheds in the region, and the subwatershed perspective allowed for identification for the most sensitive basins throughout the Sierra Nevada.     

Overland and near-surface transport of Cryptosporidium parvum from vegetated and nonvegetated surfaces

Understanding microbial pathogen transport patterns in overland flow is important for developing best management practices for limiting microbial transport to water resources. Knowledge about the effectiveness of vegetative filter strips (VFS) to reduce pathogen transport from livestock confinement areas is limited. In this study, overland and near-surface transport of Cryptosporidium parvum has been investigated. Effects of land slopes, vegetation, and rainfall intensities on oocyst transport were examined using a tilting soil chamber with two compartments, one with bare ground and the other with brome (Bromus inermis Leyss.) vegetation. Three slope conditions (1.5, 3.0, and 4.5%) were used in conjunction with two rainfall intensities (25.4 and 63.5 mm/h) for 44 min using a rainfall simulator. The vegetative surface was very effective in reducing C. parvum in surface runoff. For the 25.4 mm/h rainfall, the total percent recovery of oocysts in overland flow from the VFS varied from 0.6 to 1.7%, while those from the bare ground condition varied from 4.4 to 14.5%. For the 63.5 mm/h rainfall, the recovery percentages of oocysts varied from 0.8 to 27.2% from the VFS, and 5.3 to 59% from bare-ground conditions. For all slopes and rainfall intensities, the total (combining both surface and near-surface) recovery of C. parvum oocysts was considerably less from the vegetated surface than those from the bare-ground conditions. These results indicate that the VFS can be a best management practice for controlling C. parvum in runoff from animal production facilities.     

Characterization of Storm Flow Dynamics of Headwater Streams in the South Carolina Lower Coastal Plain1

Epps, Thomas H., Daniel R. Hitchcock, Anand D. Jayakaran, Drake R. Loflin, Thomas M. Williams, and Devendra M. Amatya, 2012. Characterization of Storm Flow Dynamics of Headwater Streams in the South Carolina Lower Coastal Plain. Journal of the American Water Resources Association (JAWRA) 1‐14. DOI: 10.1111/jawr.12000 Abstract: Hydrologic monitoring was conducted in two first‐order lower coastal plain watersheds in South Carolina, United States, a region with increasing growth and land use change. Storm events over a three‐year period were analyzed for direct runoff coefficients (ROC) and the total storm response (TSR) as percent rainfall. ROC calculations utilized an empirical hydrograph separation method that partitioned total streamflow into sustained base flow and direct runoff components. ROC ratios ranged from 0 to 0.32 on the Upper Debidue Creek (UDC) watershed and 0 to 0.57 on Watershed 80 (WS80); TSR results ranged from 0 to 0.93 at UDC and 0.01 to 0.74 at WS80. Variability in event runoff generation was attributed to seasonal trends in water table elevation fluctuation as regulated by evapotranspiration. Groundwater elevation breakpoints for each watershed were identified based on antecedent water table elevation, streamflow, ROCs, and TSRs. These thresholds represent the groundwater elevation above which event runoff generation increased sharply in response to rainfall. For effective coastal land use decision making, baseline watershed hydrology must be understood to serve as a benchmark for management goals, based on both seasonal and event‐based surface and groundwater interactions.