WATERSHED SCALE MODELING OF CRITICAL SOURCE AREAS OF RUNOFF GENERATION AND PHOSPHORUS TRANSPORT1

ABSTRACT: A curve number based model, Soil and Water Assessment Tool (SWAT), and a physically based model, Soil Moisture Distribution and Routing (SMDR), were applied in a headwater watershed in Pennsylvania to identify runoff generation areas, as runoff areas have been shown to be critical for phosphorus management. SWAT performed better than SMDR in simulating daily streamflows over the four‐year simulation period (Nash‐Sutcliffe coefficient: SWAT, 0.62; SMDR, 0.33). Both models varied streamflow simulations seasonally as precipitation and watershed conditions varied. However, levels of agreement between simulated and observed flows were not consistent over seasons. SMDR, a variable source area based model, needs further improvement in model formulations to simulate large peak flows as observed. SWAT simulations matched the majority of observed peak flow events. SMDR overpredicted annual flow volumes, while SWAT underpredicted the same. Neither model routes runoff over the landscape to water bodies, which is critical to surface transport of phosphorus. SMDR representation of the watershed as grids may allow targeted management of phosphorus sources. SWAT representation of fields as hydrologic response units (HRUs) does not allow such targeted management.     

EFFECTS OF WATERSHED REPRESENTATION ON RUNOFF AND SEDIMENT YIELD MODELING1

ABSTRACT: This paper evaluates the effects of watershed geometric representation (i.e., plane and channel representation) on runoff and sediment yield simulations in a semiarid rangeland watershed. A process based, spatially distributed runoff erosion model (KINEROS2) was used to explore four spatial representations of a 4.4 ha experimental watershed. The most complex representation included all 96 channel elements identifiable in the field. The least complex representation contained only five channel elements. It was concluded that oversimplified watershed representations greatly influence runoff and sediment yield simulations by inducing excessive infiltration on hillslopes and distorting runoff patterns and sediment fluxes. Runoff and sediment yield decrease systematically with decreasing complexity in watershed representation. However, less complex representations had less impact on runoff and sediment‐yield simulations for small rainfall events. This study concludes that the selection of the appropriate level of watershed representation can have important theoretical and practical implications on runoff and sediment yield modeling in semiarid environments.     

AN EVALUATION OF FLOOD FREQUENCY ESTIMATES BASED ON RAINFALL/RUNOFF MODELING1

ABSTRACT: An evaluation of flood frequency estimates simulated from a rainfall/runoff model is based on (1) computation of the equivalent years of record for regional estimating equations based on 50 small stream sites in Oklahoma and (2) computation of the bias for synthetic flood estimates as compared to observed estimates at 97 small stream sites with at least 20 years of record in eight eastern states. Because of the high intercorrelation of synthetic flood estimates between watersheds, little or no regional (spatial) information may be added to the network as a result of the modeling activity. The equivalent years of record for the regional estimating equations based totally on synthetic flood discharges is shown to be considerably less than the length of rainfall record used to simulate the runoff. Furthermore, the flood estimates from the rainfall/runoff model consistently underestimate the flood discharges based on observed record, particularly for the larger floods. Depending on the way bias is computed, the synthetic estimate of the 100-year flood discharge varies from 11 to 29 percent less than the value based on observed record. In addition, the correlation between observed and synthetic flood frequency estimates at the same site is also investigated. The degree of correlation between these estimates appears to vary with recurrence interval. Unless the correlation between these two estimates is known, it is not possible to compute a weighted estimate with minimum variance.     

PARTITIONING SMALL SCALE SPATIAL VARIABILITY OF RUNOFF AND EROSION ON SAGEBRUSH RANGELAND1

ABSTRACT: Most hydrologic models require input parameters which represent the variability found across an entire landscape. The estimation of such parameters is very difficult, particularly on rangeland. Improved model parameter estimation procedures are needed which incorporate the small-scale and temporal variability found on rangeland. This study investigates the use of a surface soil classification scheme to partition the spatial variability in hydrologic and interrill erosion processes in a sagebrush plant community. Four distinct microsites were found to exist within the sagebrush coppice-dune dune-interspace complex. The microsites explained the majority of variation in hydrologic and interrill erosion response found on the site and were discernable based on readily available soil and vegetation information. The variability within each microsite was quite low and was not well correlated with soil and vegetation properties. The surface soil classification scheme defined in this study can be quite useful for defining sampling procedures, for understanding hydrologic and erosion processes, and for parameterizing hydrologic models for use on sagebrush range-land.     

EVALUATION OF EROSION CONTROL PRODUCTS WITH AND WITHOUT ADDED POLYACRYLAMIDE1

ABSTRACT: Polyacrylamide (PAM) has been demonstrated to greatly reduce erosion in furrow irrigation, but much less is known about its effectiveness on the much steeper slopes typical of construction sites. The purpose of this study was to determine if anionic PAM would enhance erosion control either alone on bare soil or in combination with four types of ground covers commonly used for grass establishment: straw, straw erosion control blanket (ECB), wood fiber, and mechanically bonded fiber matrix (MBFM). Tests were conducted under natural rainfall and vegetation on a 4 percent slope (bare soil, straw, ECB, and MBFM) or using a rainfall simulator (bare soil, straw, wood fiber, MBFM) on either 10 percent or 20 percent slope on three different soil substrates. All ground cover treatments were evaluated with and without PAM applied in solution at 19 kg/ha. The straw, ECB, and MBFM significantly reduced runoff volume, average turbidity, and total sediment lost over five rainfall events on the vegetated plots. The addition of PAM to ground covers only occasionally had significant effects on runoff parameters but did significantly increase vegetative coverage overall. The rainfall simulator tests produced similar results after four events, with the straw, wood fiber, and MBFM all having significantly lower turbidity than the bare soil. The PAM significantly reduced turbidity for both the first and second events but did not consistently improve runoff quality after multiple rainfall events for any ground cover‐soil combinations tested. Separate tests of PAM applied before or after straw did not indicate a clear advantage of either approach, but runoff turbidity was often significantly reduced with PAM, especially at the 20 percent slope. Turbidity reductions were attributed to flocculation of eroded sediment.     

EFFECTS OF SOIL TYPE,PLOT CONDITION,AND SLOPE ON RUNOFF AND INTERRILL EROSION OF TWO SOILS IN THE LAKE TAHOE BASIN1

ABSTRACT: Few studies have addressed sediment discharge due to interrill erosion from natural and minimally disturbed alpine and subalpine forested watersheds. Infiltration, runoff, and surface erosion of two Tahoe Basin soils under several conditions were investigated using rainfall simulation. A significant three-way interaction among soil type, plot condition, and slope was identified. Although high erodibiity was commonly associated with disturbance and/or high slope, this was not always the case. Soil type, plot condition, slope, and duration of the event were all found to be important factors in determining the amount of erosion. Decreased water clarity in Lake Tahoe has been partly attributed to increased algal growth associated with surface runoff and erosion from adjacent watersheds. Interpretive evaluation for resource management planning should be event based and carefully delineated on a sitespecific basis.     

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.     

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.     

RUNOFF MODELING OF A MOUNTAINOUS CATCHMENT USING TOPMODEL: A CASE STUDY1

ABSTRACT: The rainfall‐runoff response of the Tygarts Creek Catchment in eastern Kentucky is studied using TOPMODEL, a hydrologic model that simulates runoff at the catchment outlet based on the concepts of saturation excess overland flow and subsurface flow. Unlike the traditional application of this model to continuous rainfall‐runoff data, the use of TOPMOEL in single event runoff modeling, specifically floods, is explored here. TOPMODEL utilizes a topographic index as an indicator of the likely spatial distribution of rainfall excess generation in the catchment. The topographic index values within the catchment are determined using the digital terrain analysis procedures in conjunction with digital elevation model (DEM) data. Select parameters in TOPMODEL are calibrated using an iterative procedure to obtain the best‐fit runoff hydrograph. The calibrated parameters are the surface transmissivity, T_O, the transmissivity decay parameter, m, and the initial moisture deficit in the root zone, S_r0. These parameters are calibrated using three storm events and verified using three additional storm events. Overall, the calibration results obtained in this study are in general agreement with the results documented from previous studies using TOPMODEL. However, the parameter values did not perform well during the verification phase of this study.     

MODELING UPLAND EROSION1

ABSTRACT: We have developed a computer model of soil loss on an upland watershed from the continuity considerations for sediment transport and from equations describing rill and interrill erosion. The model is based on dividing the upland area into a grid containing rill and interill zones, on the Universal Soil Loss Equation (USLE), and on equations describing detachment and transport capacity of rill flow. The USLE estimates the sediment load from the contributing areas. The location and amount of total erosion and deposition are determined by comparing the transport and detachment capacity of rill flow for specific storms. The model considers the mechanics of erosion process and can serve as basis for reservoir and channel design and land use planning.     

COMPARATIVE EVALUATION OF LAND COVER DATA SOURCES FOR EROSION PREDICTION1

ABSTRACT: A fundamental problem in protecting surface drinking water supplies is the identification of sites highly susceptible to soil erosion and other forms of nonpoint source (NPS) pollution. The New York City Department of Environmental Protection is trying to identify erodible sites as part of a program aimed at avoiding costly filtration. New York City's 2,000 square mile watershed system is well suited for analysis with geographic information systems (GIS); an increasingly important tool to determine the spatial distribution of sensitive NPS pollution areas. This study used a GIS to compare three land cover sources for input into the Modified Universal Soil Loss Equation (MUSLE), a model estimating soil loss from rangeland and forests, for a tributary watershed within New York City's water supply system. Sources included both conventional data (aerial photography) and Landsat data (MSS and TM images). Although land cover classifications varied significantly across these sources, location-specific and aggregate watershed predictions of the MUSLE were very similar. We conclude that using Landsat TM imagery with a hybrid classification algorithm provides a rapid, objective means of developing large area land cover databases for use in the MUSLE, thus presenting an attractive alternative to photo interpretation.     

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.     

RUNOFF AND EROSION RESPONSE OF SIMULATED WASTE BURIAL COVERS IN A SEMI-ARID ENVIRONMENT1

ABSTRACT: Control of runoff (reducing infiltration) and erosion at shallow land burials is necessary in order to assure environmentally safe disposal of low-level radioactive-waste and other waste products. This study evaluated the runoff and erosion response of two perennial grass species on simulated waste burial covers at Idaho National Engineering and Environmental Laboratory (INEEL). Rainfall simulations were applied to three plots covered by crested wheatgrass [Agropyron desertorum(Fischer ex Link) Shultes], three plots covered by streambank wheatgrass [Elymus lanceolatus(Scribner and Smith) Gould spp. lanceolaus], and one bare plot. Average total runoff for rainfall simulations in 1987, 1989, and 1990 was 42 percent greater on streambank wheatgrass plots than on crested wheatgrass plots. Average total soil loss for rainfall simulations in 1987 and 1990 was 105 percent greater on streambank wheatgrass plots than on crested wheatgrass plots. Total runoff and soil loss from natural rainfall and snowmelt events during 1987 were 25 and 105 percent greater, respectively, on streambank wheatgrass plots than on crested wheatgrass plots. Thus, crested wheatgrass appears to be better suited in revegetation of waste burial covers at INEEL than streambank wheatgrass due to its much lower erosion rate and only slightly higher infiltration rate (lower runoff rate).     

COMPARISON OF 1-D AND 2-D MODELING OF OVERLAND RUNOFF AND SEDIMENT TRANSPORT1

ABSTRACT: One-dimensional and two-dimensional modeling approaches were compared for their abilities in predicting overland runoff and sediment transport. Both 1-D and 2-D models were developed to test the hypothesis that the 2-D modeling approach could improve the model predictions over the 1-P approach, based on the same mathematical representations of physical processes for runoff and sediment transport. The models developed in this study were applied to overland areas with cross slopes. A hypothetical case and an experimental study reported by Storm (1991) were used. Based on the simulation results from the selected hypothetical case and experimental study, the 2-D model provided better representation of spatial distribution of flow depths and sediment concentrations than the 1-D model. However, no significant differences in predictions of total runoff volume and sediment yield at the outlet area were found between the 1-D and 2-D models.     

MODELING NORTH CASCADES ANNUAL RUNOFF DISTRIBUTION USING DRAINAGE BASIN ELEVATION DISTRIBUTION1

ABSTRACT: The annual distribution of flow in a drainage basin within a given region is a function of many factors. These may include annual distribution of rainfall, basin orientation, ground cover, or presence of glaciers. Since the North Cascades region of northern Washington State has little variation in precipitation distribution by month, and the region has significant snowpack, one would predict that in an unregulated basin, basin elevation would be one of the most important factors impacting an annual hydrograph distribution. Such a prediction can be made since the higher a drainage basin is, the larger the portion of runoff that would occur as late spring snowmelt. Given that there is a relationship between elevation distribution and annual hydrograph, the problem becomes one of how to use this relationship to model an ungaged basin's hydrograph. This study concludes that, within the North Cascades region and perhaps within other regions, an effective method of determining annual flow distribution is to model ungaged flows in the same manner as flows from a gaged basin with an elevation distribution similar to that of the subject basin.     

COMPARATIVE EVALUATION OF THREE URBAN RUNOFF MODELS1

ABSTRACT: Three urban runoff models, namely, the Road Research Laboratory Model (RRLM), the Storm Water Management Model (SWMM) and the University of Cincinnati Urban Runoff Model (UCURM), were examined by comparing the model simulated hydrographs with the hydrographs measured on several instrumented urban watersheds. This comparison was done for the hydrograph peak points as well as for the entire hydrographs using such statistical measures as the correlation coefficient, the special correlation coefficient and the integral square error. The results of the study indicated that, when applying the three selected non-calibrated models on small urban catchments, the SWM model performed marginally better than the RRL model and both these models were more accurate than the UCUR model. On larger watersheds, the comparisons between the SWM model and the other two models would be likely even more favourable for the SWM model, because it has the most advanced flow routing scheme among the studied models.     

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.     

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.     

RELATIONSHIP BETWEEN ANNUAL RUNOFF AND WATERSHED AREA FOR THE EASTERN UNITED STATES1

ABSTRACT: As part of the U.S. Environmental Protection Agency's effort to determine the long-term effects of acidic deposition on surface water chemistry, annual runoff was estimated for about 1000 ungaged sites in the eastern U.S. using runoff contour maps. One concern in using contour maps was that a bias may be introduced in the runoff estimates due to the size of the 1000 ungaged sites relative to the size of the watersheds used in developing the maps. To determine if a bias was present the relationship between the annual runoff (expressed as depth) and the watershed area for the Northeast (NE) and Southern Blue Ridge Province (SBRP) was tested using five regional data bases. One short-term data base (1984 Water Year, n = 531) and two long-term data bases (1940–57, n = 134 and 1951–80, n = 342) were used in the NE. In the SBRP one short-term database (1984 Water Year, n = 531) and one long-term data base (1951–80, n = 60) were used. For the NE and the SBRP, runoff was not directly correlated with watershed area using the five regional databases. Also, runoff normalized by precipitation was not related to watershed area.     

PROJECTION OF URBANIZATION EFFECTS ON RUNOFF USING CLARK INSTANTANEOUS UNIT HYDROGRAPH PARAMETERS1

ABSTRACT: To alleviate serious flooding problems brought upon by rapid urbanization in the Beargrass Creek watershed, located in Louisville, Kentucky, the U.S. Army Corps of Engineers undertook a major flood study in 1973. In order to predict flood conditions in 1990, the year when the watershed was expected to undergo complete urbanization, trends in the Clark Instantaneous Unit Hydrograph (Clark IUH) parameters were utilized to determine the 1990 unit hydrograph and flood conditions. Based on the results from this flood study, this paper demonstrates the applicability of using projected Clark IUH parameters for modeling future runoff conditions in an urbanizing watershed. Values of these parameters, as estimated from maximum annual historical flood data, are used to develop regression models for predicting future Clark IUH parameters. Using the projected parameters, selected annual flood events since 1973 are simulated in order to verify the accuracy of these projections. Results show a close correspondence between the simulated and observed flood characteristics. Hence, the use of projected Clark IUH parameters is an appropriate procedure for modeling future runoff conditions in an urbanizing watershed.