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.     

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.     

EVALUATION OF RUNOFF,EROSION, AND PHOSPHORUS MODELING SYSTEM—SIMPLE1

ABSTRACT: The purpose of this study was to evaluate the performance of Spatially Integrated Models for Phosphorus Loading and Erosion (SIMPLE) in predicting runoff volume, sediment loss, and phosphorus loading from two watersheds. The modeling system was applied to the 334 ha QOD subwatershed, part of the Owl Run watershed, located in Fauquier County, Virginia, and to the 2240 ha watershed, Battle Branch, located in Delaware County, Oklahoma. Simulation runs were conducted at cell and field scales, and simulation results were compared with observed data. Runoff volume and dissolved phosphorus loading were measured at the Battle Branch watershed. Runoff volume, sediment yield, and total phosphorus loading were measured at the QOD site. SIMPLE tended to underestimate runoff volumes during the dormant period, from November to March. The comparison between observed and predicted dissolved phosphorus showed better correlation than for observed and predicted total phosphorus loading. Cell level simulations provided similar estimates of runoff volume and phosphorus loading when compared to field level simulations for both watersheds. However, observed sediment yields better compared with the values predicted from the cell level simulation when compared to field level simulation. Finally, results of model evaluation indicated that SIMPLE's predictive ability is acceptable for screening applications but not for site-specific quantitative predictions.     

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.     

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.     

A FRAMEWORK FOR PHOSPHORUS TRANSPORT MODELING IN THE LAKE OKEECHOBEE WATERSHED1

ABSTRACT: A modeling framework was developed to determine phosphorus loadings to Lake Okeechobee from watersheds located north of the lake. This framework consists of the land-based model CREAMS-WT, the in-stream transport model QUAL2E, and an interface procedure to format the land-based model output for use by the in-stream model. QUAL2E hydraulics and water quality routines were modified to account for flow routing and phosphorus retention in both wetlands and stream channels. Phosphorus loadings obtained from previous applications of CREAMS-WT were used by QUAL2E, and calibration and verification showed that QUAL2E accurately simulated seasonal and annual phosphorus loadings from a watershed. Sensitivity and uncertainty analyses indicated that the accuracy of monthly loadings can be improved by using better estimates of in-stream phosphorus decay rates, ground water phosphorus concentrations, and runoff phosphorus concentrations as input to QUAL2E.     

MODELING OF NONPOINT SOURCE POLLUTION OF NITROGEN AT THE WATERSHED SCALE1

ABSTRACT: The Watershed Nutrient Transport and Transformation (NTT-Watershed) model is a physically based, energy-driven, multiple land use, distributed model that is capable of simulating water and nutrient transport in a watershed. The topographic features and subsurface properties of the watershed are refined into uniform, homogeneous square grids. The vertical discretization includes vegetation, overland flow, soil water redistribution and groundwater zones. The chemical submodel simulates the nitrogen dynamics in terrestrial and aquatic systems. Three chemical state variables are considered (NO₃^-_-, NH₄⁺, and Org-N). The NTT-Watershed model was used to simulate the fate and transport of nitrogen in the Muddy Brook watershed in Connecticut. The model was shown to be capable of capturing the hydrologic and portions of the nitrogen dynamics in the watershed. Watershed planners could use this model in developing strategies of best management practices that could result in maximizing the reductions of nitrogen export from a watershed.     

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.     

SWMHMS - SMALL WATERSHED MONTHLY HYDROLOGIC MODELING SYSTEM1

ABSTRACT: SWMHMS is a conceptual computer modeling program developed to simulate monthly runoff from a small nonurban watershed. The input needed to run model simulations include daily precipitation, monthly data for evapotranspiration determination (average temperature, crop consumptive coefficients, and percent daylight hours), and six watershed parameter values. Evapotranspiration was calculated with the Blaney-Criddle equation while surface runoff was determined using the Soil Conservation Service curve number procedure. For watershed parameter evaluation, SWMHMS provides options for both optimization and sensitivity analysis. Observed runoff data are required along with the model input previously mentioned in order to conduct parameter optimization. SWMEIMS was tested with data from six watersheds located in different regions of the United States. Model accuracy was generally found to be very good except on watersheds having substantial snowfall accumulation. In having only six watershed parameters, SWMHMS is less complex to use than many other computer programs that calculate monthly runoff. Consequently, SWMHMS may find its greatest application as an educational tool for students learning principles of hydrologic modeling, such as parameter evaluation procedures and the impacts of input data uncertainty on model results.     

IDENTIFICATION OF CRITICAL NONPOINT POLLUTION SOURCE AREAS USING GEOGRAPHIC INFORMATION SYSTEMS AND WATER QUALITY MODELING1

n integrated approach coupling water quality computer simulation modeling with a geographic information system (GIS) was used to delineate critical areas of nonpoint source (NPS) pollution at the watershed level. Two simplified pollutant export models were integrated with the Virginia Geographic Information System (VirGIS) to estimate soil erosion, sediment yield, and phosphorus (P) loading from the Nomini Creek watershed located in Westmoreland County, Virginia. On the basis of selected criteria for soil erosion rate, sediment yield, and P loading, model outputs were used to identily watershed areas which exhibit three categories (low, medium, high) of non-point source pollution potentials. The percentage of the watershed area in each category, and the land area with critical pollution problems were also identified. For the 1505-ha Nomini Creek watershed, about 15, 16, and 21 percent of the watershed area were delineated as sources of critical soil erosion, sediment, and phosphorus pollution problems, respectively. In general, the study demonstrated the usefulness of integrating GIS with simulation modeling for nonpoint source pollution control and planning. Such techniques can facilitate making priorities and targeting nonpoint source pollution control programs.     

MODELING RUNOFF FROM VARIABLE SOURCE AREAS IN HUMID,SHALLOW WATER TABLE ENVIRONMENTS1

ABSTRACT: Variable Source Areas (VSAs) are zones with water saturated soils in forested wetlands fringing streams and creeks. Runoff from these areas is generated by saturation excess after a shallow water table rises and inundates the ground surface. In humid regions, like Florida and the Southeast, VSAs are believed to produce most of the runoff in shallow water table environments. Modeling the spatial extent and temporal fluctuation of a VSA is difficult because the formation of a VSA depends on a number of hydrological and morphological factors like rainfall intensity, soil texture, water table depth, and topographic attributes of the terrain. In this paper, we couple a digital elevation model with a two‐dimensional variable saturation model to illustrate the formation of a VSA at the hillside scale. The topography derived from the digital elevation model forms the upper domain geometry for the two‐dimensional finite element simulations of variable saturated flow. The objectives are: (1) to model the spatial and dynamic fluctuation of a VSA, and (2) to understand the roles of rainfall variability and terrain attributes on the formation of a VSA. Results show that hillsides with shallow water table depths, low saturated hydraulic conductivity, mild slopes, and concave slope curvature were more susceptible to runoff from a variable source. Runoff from a variable source showed little sensitivity to rainfall intensity. In general, landscapes with steep slopes generated a small VSA and a seepage face that vanished rapidly with time. In contrast, flat terrains are more amenable to VSA and retain ground surface inundation for longer periods of time.     

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.     

PATTERNS OF WATERSHED MONTHLY RUNOFF1

ABSTRACT: Two dimensional sliding polynomials were adapted to pattern analysis of watershed monthly rainfall and runoff. Contours of runoff in the two-dimensional space of time and rainfall are constructed on a grid of 16 nodes whose values are determined by least squares. This method is form free, hence derived patterns are not biased to selected functional forms, but can directly represent the smoothed data. Values of the nodes are localized averages of the data constrained by required mathematical continuity across the grid of values. An advantage of the method is that the standard deviation can be calculated for each node, thus producing patterns of uncertainty of the deterministic component revealed by the data.     

APPLICATION OF THE GREEN-AMPT INFILTRATION EQUATION TO WATERSHED MODELING1

ABSTRACT: A computer model was developed, based on the Green-Ampt infiltration equation, to computed rainfall excess for a single precipitation event. The model requires an estimate of parameters related to hydraulic conductivity, wetting front section, and fillable porosity of the soil layers. Values of parameters were estimated from soil textural averages or regression equations based on percent sand, percent clay, and porosity. Average values of effective porosity and wetting front suction were largely acceptable due to the relatively low variability and low model sensitivity to the parameters. Hydraulic conductivity was the most erratic constituent of the loss rate computation due to the high variability and the high sensitivity of the computed infiltration to the parameter. The performance of the Green-Ampt infiltration model was tested through a comparison with the SCS curve number procedure. Seven watersheds and 23 storms with precipitation of one inch or greater were used in the comparison. For storms with less than one inch of rainfall excess, the SCS curve number procedure generally gave the best results; however, for six of the seven storms with precipitation excess greater than one inch, the Green-Ampt procedure delivered better results. In this comparison, both procedures used the same initial abstractions. The separation of rainfall losses into infiltration, interception, and surface retention is, in theory, an accurate method of estimating precipitation excess. In the second phase of the study using nine watersheds and 39 storms, interception and surface retention losses were computed by the Horton equations. Green-Ampt and interception parameters were estimated from value sin the literature, while the surface retention parameter was calibrated so that the computed runoff volumes matched observed volumes. A relationship was found between the surface retention storage capacity and the 15-day antecedent precipitation index, month of year, and precipitation amount.     

HYDROLOGICAL MODELING OF THE IROQUOIS RIVER WATERSHED USING HSPF AND SWAT1

ABSTRACT: The performance of two popular watershed scale simulation models — HSPF and SWAT — were evaluated for simulating the hydrology of the 5,568 km² Iroquois River watershed in Illinois and Indiana. This large, tile drained agricultural watershed provides distinctly different conditions for model comparison in contrast to previous studies. Both models were calibrated for a nine‐year period (1987 through 1995) and verified using an independent 15‐year period (1972 through 1986) by comparing simulated and observed daily, monthly, and annual streamflow. The characteristics of simulated flows from both models are mostly similar to each other and to observed flows, particularly for the calibration results. SWAT predicts flows slightly better than HSPF for the verification period, with the primary advantage being better simulation of low flows. A noticeable difference in the models' hydrologic simulation relates to the estimation of potential evapotranspiration (PET). Comparatively low PET values provided as input to HSPF from the BASINS 3.0 database may be a factor in HSPF's overestimation of low flows. Another factor affecting baseflow simulation is the presence of tile drains in the watershed. HSPF parameters can be adjusted to indirectly account for the faster subsurface flow associated with tile drains, but there is no specific tile drainage component in HSPF as there is in SWAT. Continued comparative studies such as this, under a variety of hydrologic conditions and watershed scales, provide needed guidance to potential users in model selection and application.     

ROLE OF WATERSHED SUBDIVISION ON MODELING THE EFFECTIVENESS OF BEST MANAGEMENT PRACTICES WITH SWAT1

Distributed parameter watershed models are often used for evaluating the effectiveness of various best management practices (BMPs). Streamflow, sediment, and nutrient yield predictions of a watershed model can be affected by spatial resolution as dictated by watershed subdivision. The objectives of this paper are to show that evaluation of BMPs using a model is strongly linked to the level of watershed subdivision; to suggest a methodology for identifying an appropriate subdivision level; and to examine the efficacy of different BMPs at field and watershed scales. In this study, the Soil and Water Assessment Tool (SWAT) model was calibrated and validated for streamflow, sediment, and nutrient yields at the outlet of the Dreisbach (623 ha) and Smith Fry (730 ha) watersheds in Maumee River Basin, Indiana. Grassed waterways, grade stabilization structures, field borders, and parallel terraces are the BMPs that were installed in the study area in the 1970s. Sediment and nutrient outputs from the calibrated model were compared at various watershed subdivision levels, both with and without implementation of these BMPs. Results for the study watersheds indicated that evaluation of the impacts of these BMPs on sediment and nutrient yields was very sensitive to the level of subdivision that was implemented in SWAT. An optimal watershed subdivision level for representation of the BMPs was identified through numerical simulations. For the study watersheds, it would appear that the average subwatershed area corresponding to approximately 4 percent of total watershed area is needed to represent the influence of these BMPs when using the SWAT model.     

APPLICATION OF AN ENVIRONMENTAL AND ECONOMIC MODELING SYSTEM FOR WATERSHED ASSESSMENTS1

ABSTRACT: A National Pilot Project (NPP) on Livestock and the Environment was initiated in 1992 to help provide solutions to environmental problems associated with livestock production. A major development of the NPP was the Comprehensive Economic and Environmental Optimization Tool‐Livestock and Poultry (CEEOT‐LP), an integrated modeling system designed to produce economic and environmental indicators for alternative policy scenarios applied to intensive livestock production watersheds. The system consists of a farm‐level economic model (FEM) and two environmental models: the field‐scale APEX model and the watershed‐level SWAT model. To date, CEEOT‐LP has been applied to two watersheds in Texas and one in Iowa. Predicted reductions in P losses for two P‐based manure application rate scenarios, relative to baseline conditions, ranged from ?4 to ?54 percent across the three watersheds; however, N loss impacts ranged from a decrease of 34 percent to an increase of 79 percent. For five other alternative scenarios that were simulated for only one watershed, N and P loss impacts ranged between a reduction of 78 percent to an increase of 20 percent. Aggregate watershed‐level economic impacts of the seven scenarios spanned a spectrum of a 27 percent decrease to a 25 percent increase in profit, relative to the baseline.     

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.     

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.     

UNCERTAINTY ANALYSIS OF RUNOFF ESTIMATES FROM A RUNOFF CONTOUR MAP1

ABSTRACT: The U.S. Environmental Protection Agency (EPA) in cooperation with the U.S. Geological Survey (USGS) conducted an analysis to quantify the uncertainty associated with interpolating runoff to specific sites using a runoff contour map. We interpolated runoff to 93 gaged watersheds from a runoff contour map using (1) hand interpolation to the watershed outlet, (2) a computer interpolation to the watershed outlet, and (3) hand interpolation to the watershed centroid. We compared the interpolated values to the actual gaged values and found that there was a bias in the average interpolated value for runoff estimated at basin outlets, with interpolated values being less than the actual. We found no significant difference between the hand interpolation method and the computer interpolation method except that the computer method tended to have higher variability due to factors inherent to the software used. There were no strong spatial correlations or regional patterns in the runoff interpolations, which indicates that there are no regional biases introduced in the development of the contour map. We determined that we could estimate runoff, on the average, within approximately 8.9 cm (3.5 in; 15 percent) of the measured value using the three methods. The results of this work indicate that runoff contour maps can he used in regional studies to estimate runoff to ungaged systems with quantifiable uncertainty.