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.     

EQUIVALENCE BETWEEN TOPMODEL AND THE NRCS CURVE NUMBER METHOD IN PREDICTING VARIABLE RUNOFF SOURCE AREAS1

An equivalence is proposed between two rainfall‐runoff methods with a long history of use in the United States and Europe. In watersheds where variable source areas dominate runoff, the two methods can have comparable probability distribution functions of moisture deficit, and therefore predict similar saturated runoff source areas. A novel approach is introduced to determine the S parameter in the Natural Resources Conservation Service (NRCS) method. This approach constrains S by the physical soil and topography characteristics of the catchment and depth to water table. The NRCS curve number method is at the core of many rainfall‐runoff models in hydrology. As a simple lumped parameter method, it is often scrutinized because it is not obvious how to derive S from catchment hydromorphological characteristics. The novel approach provides a clear physical meaning for S, allowing better estimation of this parameter in humid shallow water table environments where the variable source area can be the dominant runoff mechanism.     

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.     

Effects of Watershed Impervious Cover on Dissolved Silica Loading in Storm Flow1

Abstract: Dissolved silica (DSi) availability is a factor that affects the composition of algal populations in aquatic ecosystems. DSi cycling is tightly linked to the hydrological cycle, which is affected by human alterations of the landscape. Development activities that increase impervious cover change watershed hydrology and may increase the discharge of DSi‐poor rainwater and decrease the discharge of DSi‐rich ground water into aquatic ecosystems, possibly shifting algal community composition toward less desirable assemblages. In this study, DSi loadings from two adjacent coastal watersheds with different percent impervious cover were compared during four rain and five nonrain events. Loadings in the more impervious watershed contained a significantly larger proportion of surface runoff than base flow (ground‐water discharge) and had lower [DSi] water during rain events than the less impervious watershed. Application of the Soil Conservation Service Curve Number (CN) method showed that the minimum rainfall height necessary to yield runoff was significantly lower for the more impervious watershed, implying that runoff volumes increase with impervious cover as well as the frequency of runoff‐yielding events. Empirical data collected during this study and estimates derived from the CN method suggest that impervious cover may be responsible for both short‐term DSi limitation during rain events as well as long‐term reduction of DSi inputs into aquatic ecosystems.     

Runoff Curve Numbers at the Agricultural Field‐Scale and Implications for Continuous Simulation Modeling

This study used monitoring in the waterways of agricultural fields to understand the use of the runoff curve number (CN) in continuous simulation models. The CN has a long history as a design tool for estimating runoff volumes for large, single storms on small watersheds, but its use in continuous simulation models to describe runoff from smaller storms and relatively small areas is more recent and controversial. We examined 788 nonwinter rainfall events on four agricultural fields over five years (2004‐2008) during which runoff was generated in 87 events. The largest 20 runoff events on each field generated approximately 90% of the total runoff volume. The runoff event CNs showed an inverse correlation with storm depth that could not consistently be explained by previous precipitation. We review how small areas of higher runoff generation within larger areas will systematically increase the apparent CN of the larger area as the storm size decreases. If this variation is not incorporated into a model explicitly, continuous simulation modelers must understand that when source areas are aggregated or when runoff generation is spatially variable, the overall CN is not unique when smaller storms are included in the calibration set.     

GIS‐Based Predictive Models of Hillslope Runoff Generation Processes1

Abstract: Successful nonpoint source pollution control using best management practice placement is a complex process that requires in‐depth knowledge of the locations of runoff source areas in a watershed. Currently, very few simulation tools are capable of identifying critical runoff source areas on hillslopes and those available are not directly applicable under all runoff conditions. In this paper, a comparison of two geographic information system (GIS)‐based approaches: a topographic index model and a likelihood indicator model is presented, in predicting likely locations of saturation excess and infiltration excess runoff source areas in a hillslope of the Savoy Experimental Watershed located in northwest Arkansas. Based on intensive data collected from a two‐year field study, the spatial distributions of hydrologic variables were processed using GIS software to develop the models. The likelihood indicator model was used to produce probability surfaces that indicated the likelihood of location of both saturation and infiltration excess runoff mechanisms on the hillslope. Overall accuracies of the likelihood indicator model predictions varied between 81 and 87% for the infiltration excess and saturation excess runoff locations respectively. On the basis of accuracy of prediction, the likelihood indicator models were found to be superior (accuracy 81‐87%) to the predications made by the topographic index model (accuracy 69.5%). By combining statistics with GIS, runoff source areas on a hillslope can be identified by incorporating easily determined hydrologic measurements (such as bulk density, porosity, slope, depth to bed rock, depth to water table) and could serve as a watershed management tool for identifying critical runoff source areas in locations where the topographic index or other similar methods do not provide reliable results.     

STORM RUNOFF PREDICTION BASED ON A SPATIALLY DISTRIBUTED TRAVEL TIME METHOD UTILIZING REMOTE SENSING AND GIS1

ABSTRACT: In this study, remotely sensed data and geographic information system (GIS) tools were used to estimate storm runoff response for Simms Creek watershed in the Etonia basin in northeast Florida. Land cover information from digital orthophoto quarter quadrangles (DOQQ), and enhanced thematic mapper plus (ETM+) were analyzed for the years 1990, 1995, and 2000. The corresponding infiltration excess runoff response of the study area was estimated using the U.S. Department of Agriculture (USDA), Natural Resources Conservation Service Curve Number (NRCS‐CN) method. A digital elevation model (DEM)/GIS technique was developed to predict stream response to runoff events based on the travel time from each grid cell to the watershed outlet. A comparison of predicted to observed stream response shows that the model predicts the total runoff volume with an efficiency of 0.98, the peak flow rate at an efficiency of 0.85, and the full direct runoff hydrograph with an average efficiency of 0.65. The DEM/GIS travel time model can be used to predict the runoff response of ungaged watersheds and is useful for predicting runoff hydrographs resulting from proposed large scale changes in the land use.     

Hydrologic Landscape Classification to Estimate Bristol Bay,Alaska Watershed Hydrology

Hydrologic landscapes (HLs) have proven to be a useful tool for broad scale assessment and classification of landscapes across the United States as they help organize larger geographical areas into areas of similar hydrologic characteristics. We developed a HL classification for the Bristol Bay watershed of southwest Alaska that incorporates indices of annual climate and seasonality, terrain, geology, and the influences of large lakes and glaciers. A HL classification is particularly useful in this large watershed because of its hydrologic and landscape variability, important salmon fishery, variety of environmental and potential anthropogenic stressors, and lack of widespread hydrologic data. Following creation of Bristol Bay basin‐wide HL classes, we compared the HL distributions within watersheds grouped by two calculated runoff parameters derived from available long‐term streamflow records and found HL distributions within these groups provided predictive insight on hydrologic behavior. Using these developed runoff groups, we estimated expected hydrologic behavior in watersheds across the larger Bristol Bay watershed that lacked gauged streamflow records. The HL approach provides a scientific basis for estimating the first‐order hydrologic behavior of watersheds and landscapes that lack detailed hydrologic information.     

COMPUTATION OF RUNOFF REDUCTION CAUSED BY FARM PONDS1

ABSTRACT: A two-parameter farm pond storage index, FPSI, was Used to adjust computed surface. runoff using the partial area runoff contribution resulting from runoff captured by farm ponds. The validity of the index method was tested by fitting a continuous accounting version of the Soil Conservation Service curve number procedure to surface runoff data from each of three watersheds, first with and then without the FPSI routine. Evapotranspiration computed with the Jensen-Haise method and rainfall were input to the model. A linear relationship was assumed between the storage index and the portion of the controlled drainage area that was contributing to runoff. Adjusting the computed runoff with the FPSI reduced the coefficient of variation of monthly measured versus computed surface runoff for each of the three watersheds. The correlation coefficients for the same comparisons were increased. The annual predicted surface runoff Was improved for 12 of the 17 station years of data tested. The farm pond storage index could be used with any surface runoff model to improve the prediction of runoff from watersheds with drainage areas greater than 1 square mile and with about 20 percent or more of the drainage area controlled by farm ponds.     

VALIDATION OF AGNPS FOR SMALL WATERSHEDS USING AN INTEGRATED AGNPS/GIS SYSTEM1

ABSTRACT: The AGNPS (AGricultural NonPoint Source) model was evaluated for predicting runoff and sediment delivery from small watersheds of mild topography. Fifty sediment yield events were monitored from two watersheds and five nested subwater-sheds in East Central Illinois throughout the growing season of four years. Half of these events were used to calibrate parameters in the AGNPS model. Average calibrated parameters were used as input for the remaining events to obtain runoff and sediment yield data. These data were used to evaluate the suitability of the AGNPS model for predicting runoff and sediment yield from small, mild-sloped watersheds. An integrated AGNPS/GIS system was used to efficiently create the large number of data input changes necessary to this study. This system is one where the AGNPS model was integrated with the GRASS (Geographic Resources Analysis Support System) GIS (Geographical Information System) to develop a decision support tool to assist with management of runoff and erosion from agricultural watersheds. The integrated system assists with the development of input GIS layers to AGNPS, running the model, and interpretation of the results.     

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.     

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.     

HYDROLOGIC PARAMETERIZATION OF WATERSHEDS FOR RUNOFF PREDICTION USING SWAT1

ABSTRACT: The use of continuous time, distributed parameter hydrologic models like SWAT (Soil and Water Assessment Tool) has opened several opportunities to improve watershed modeling accuracy. However, it has also placed a heavy burden on users with respect to the amount of work involved in parameterizing the watershed in general and in adequately representing the spatial variability of the watershed in particular. Recent developments in Geographical Information Systems (GIS) have alleviated some of the difficulties associated with managing spatial data. However, the user must still choose among various parameterization approaches that are available within the model. This paper describes the important parameterization issues involved when modeling watershed hydrology for runoff prediction using SWAT with emphasis on how to improve model performance without resorting to tedious and arbitrary parameter by parameter calibration. Synthetic and actual watersheds in Indiana and Mississippi were used to illustrate the sensitivity of runoff prediction to spatial variability, watershed decomposition, and spatial and temporal adjustment of curve numbers and return flow contribution. SWAT was also used to predict stream runoff from actual watersheds in Indiana that have extensive subsurface drainage. The results of this study provide useful information for improving SWAT performance in terms of stream runoff prediction in a manner that is particularly useful for modeling ungaged watersheds wherein observed data for calibration is not available.     

EFFECT OF RAINFALL-EXCESS DETERMINATION ON RUNOFF COMPUTATION1

ABSTRACT: Four methods of determining rainfall-excess are considered. They include Ø-index and equations of Horton, Kostyakov and Philip. These methods are utilized in a nonlinear kinematic wave model to predict surface runoff from two natural agricultural watersheds. The effect of determining rainfall-excess on surface runoff response is then examined. It is observed that errors in rainfall-excess constitute a major source of error in runoff prediction. The choice of a method of determining rainfall-excess is thus crucial to runoff computation.     

Annual Estimates of Recharge,Quick‐Flow Runoff,and Evapotranspiration for the Contiguous U.S. Using Empirical Regression Equations

This study presents new data‐driven, annual estimates of the division of precipitation into the recharge, quick‐flow runoff, and evapotranspiration (ET) water budget components for 2000‐2013 for the contiguous United States (CONUS). The algorithms used to produce these maps ensure water budget consistency over this broad spatial scale, with contributions from precipitation influx attributed to each component at 800 m resolution. The quick‐flow runoff estimates for the contribution to the rapidly varying portion of the hydrograph are produced using data from 1,434 gaged watersheds, and depend on precipitation, soil saturated hydraulic conductivity, and surficial geology type. Evapotranspiration estimates are produced from a regression using water balance data from 679 gaged watersheds and depend on land cover, temperature, and precipitation. The quick‐flow and ET estimates are combined to calculate recharge as the remainder of precipitation. The ET and recharge estimates are checked against independent field data, and the results show good agreement. Comparisons of recharge estimates with groundwater extraction data show that in 15% of the country, groundwater is being extracted at rates higher than the local recharge. These maps of the internally consistent water budget components of recharge, quick‐flow runoff, and ET, being derived from and tested against data, are expected to provide reliable first‐order estimates of these quantities across the CONUS, even where field measurements are sparse.     

Relationship of Water Quality and Pollutant Loads to Land Uses in adjoining watersheds1

ABSTRACT: The Grand and Saugeen Rivers in southern Ontario were chosen for study as pilot watersheds under the Pollution From Land Use Activities Reference Group (PLUARG) study. The pilot watersheds have adjacent headwater areas and are physically similar in geology, physiography, and climate. Significant differences in water quality between the watersheds at their outlets are attributed to land use and population differences. The major pollutant sources in the two pilot watersheds were identified as trace elements from urban runoff and point source discharges; phosphorus from agricultural and urban runoff and private waste disposal; chloride from transportation corridors; and sediment and nitrogen from agricultural runoff. Yields at the watershed outlets were similar for suspended sediment and two to three times as high in the Grand River for phosphours, nitrogen, chloride, and lead. The higher phosphorus and nitrogen levels were attributed to larger point source inputs and the higher proportion of agricultural activity, comprising 75 percent of the Grand River basin compared to 64 percent in the Saugeen River basin. Similarly, the higher chloride and lead levels were attributed to an order of magnitude larger population and three times as much urban land in the Grand River basin compared to the Saugeen River basin.     

SUSCEPTIBILITY OF INDIANA WATERSHEDS TO HERBICIDE CONTAMINATION1

ABSTRACT: An index of watershed susceptibility to surface water contamination by herbicides could be used to improve source water assessments for public drinking water supplies, prioritize watershed restoration projects, and direct funding and educational efforts to areas where the greatest environmental benefit can be realized. The goal of this study is to use streamflow and herbicide concentration data to develop and evaluate a method for estimating comparative watershed susceptibility to herbicide loss. United States Geological Survey (USGS) concentration data for five relatively water soluble herbicides (alachlor, atrazine, cyanazine, metolachlor, and simazine) were analyzed for 16 Indiana watersheds. Correlation was assessed between observed herbicide losses and: (1) a herbicide runoff index using GIS‐based land use, soil type, SCS runoff curve number, tillage practice, herbicide use estimates, and combinations of these factors; and (2) predicted herbicide losses from a non‐point source pollution model (NAPRA‐Web, an Internet‐based interface for GLEAMS). The highest adjusted R²value was found between herbicide concentration and the runoff curve number alone, ranging from 0.25 to 0.56. Predictions from the simulation model showed a poorer correlation with observed herbicide loss. This indicates potential for using the runoff curve number as a simple herbicide contamination susceptibility index.     

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.     

HYDROLOGIC RESPONSE OF MASSACHUSETTS WATERSHEDS1

Hydrologic response, defined as the annual direct runoff divided by the annual precipitation, was computed for twenty-one watersheds in or near western Massachusetts, using a total of 232 years of hydrologic records. Variability of the results over the period of analysis was greater than is desirable to inspire confidence in the usefulness of the hydrologic response function; however, the results do suggest that the hydrologic response concept, with appropriate refinements, could be applied successfully to the problem of delineating hydrologic provinces and determination of drainage and storage in unregulated watersheds.     

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.