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.     

Decreased Runoff Response to Precipitation,Little Missouri River Basin,Northern Great Plains,USA

High variability in precipitation and streamflow in the semiarid northern Great Plains causes large uncertainty in water availability. This uncertainty is compounded by potential effects of future climate change. We examined historical variability in annual and growing season precipitation, temperature, and streamflow within the Little Missouri River Basin and identified differences in the runoff response to precipitation for the period 1976‐2012 compared to 1939‐1975 (n = 37 years in both cases). Computed mean values for the second half of the record showed little change (<5%) in annual or growing season precipitation, but average annual runoff at the basin outlet decreased by 22%, with 66% of the reduction in flow occurring during the growing season. Our results show a statistically significant (p < 0.10) 27% decrease in the annual runoff response to precipitation (runoff ratio). Surface‐water withdrawals for various uses appear to account for <12% of the reduction in average annual flow volume, and we found no published or reported evidence of substantial flow reduction caused by groundwater pumping in this basin. Results of our analysis suggest that increases in monthly average maximum and minimum temperatures, including >1°C increases in January through March, are the dominant driver of the observed decrease in runoff response to precipitation in the Little Missouri River Basin.     

Spatial Distribution of Water Supply in the Coterminous United States1

Abstract: Available water supply across the contiguous 48 states was estimated as precipitation minus evapotranspiration using data for the period 1953‐1994. Precipitation estimates were taken from the Parameter‐Elevation Regressions on Independent Slopes Model (PRISM). Evapotranspiration was estimated using two models, the Advection‐Aridity model and the Zhang model. The evapotranspiration models were calibrated using precipitation and runoff data for 655 hydrologically undisturbed basins, and then tested using estimates of natural runoff for the 18 water resource regions (WRR) of the 48 contiguous states. The final water supply coverage reflects a mixture of outputs from the two evapotranspiration models. Political, administrative, and land cover boundaries were mapped over the coverage of mean annual water supply. Across the entire study area, we find that 53% of the water supply originates on forested land, which covers only 29% of the surface area, and that 24% originates on federal lands, including 18% on national forests and grasslands alone. Forests and federal lands are even more important in the West (the 11 western contiguous states), where 65% of the water supply originates on forested land and 66% on federal lands, with national forests and grasslands contributing 51%.     

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.     

NETWORK DESIGN FOR WATER SUPPLY FORECASTING IN THE WEST1

ABSTRACT: The objective is to develop techniques to evaluate how changes in basic data networks can improve accuracy of water supply forecasts for mountainous areas. The approach used was to first quantify how additional data would improve our knowledge of winter precipitation, and second to estimate how this knowledge translates, quantitatively, into improvement in forecast accuracy. A software system called DATANET was developed to analyze each specific gage network alternative. This system sets up a fine mesh of grid points over the basin. The long-term winter mean precipitation at each grid point is estimated using a simple atmospheric model of the orographic precipitation process. The mean runoff at each grid point is computed from the long-term mean precipitation estimate. The basic runoff model is calibrated to produce the observed long-term runoff. The error analysis is accomplished by comparing the error in forecasts based on the best possible estimate of precipitation using all available data with the error in the forecasts based on the best possible estimate of winter precipitation using only the gaged data. Different data network configurations of gage sites can be compared in terms of forecast errors.     

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.     

Comparison and Evaluation of Gridded Precipitation Datasets in a Kansas Agricultural Watershed Using SWAT

Gridded precipitation datasets are becoming a convenient substitute for gauge measurements in hydrological modeling; however, these data have not been fully evaluated across a range of conditions. We compared four gridded datasets (Daily Surface Weather and Climatological Summaries [DAYMET], North American Land Data Assimilation System [NLDAS], Global Land Data Assimilation System [GLDAS], and Parameter‐elevation Regressions on Independent Slopes Model [PRISM]) as precipitation data sources and evaluated how they affected hydrologic model performance when compared with a gauged dataset, Global Historical Climatology Network‐Daily (GHCN‐D). Analyses were performed for the Delaware Watershed at Perry Lake in eastern Kansas. Precipitation indices for DAYMET and PRISM precipitation closely matched GHCN‐D, whereas NLDAS and GLDAS showed weaker correlations. We also used these precipitation data as input to the Soil and Water Assessment Tool (SWAT) model that confirmed similar trends in streamflow simulation. For stations with complete data, GHCN‐D based SWAT‐simulated streamflow variability better than gridded precipitation data. During low flow periods we found PRISM performed better, whereas both DAYMET and NLDAS performed better in high flow years. Our results demonstrate that combining gridded precipitation sources with gauge‐based measurements can improve hydrologic model performance, especially for extreme events.     

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.     

A COMPARISON OF DELTA CHANGE AND DOWNSCALED GCM SCENARIOS FOR THREE MOUNTAINOUS BASINS IN THE UNITED STATES1

ABSTRACT: Simulated daily precipitation, temperature, and runoff time series were compared in three mountainous basins in the United States: (1) the Animas River basin in Colorado, (2) the East Fork of the Carson River basin in Nevada and California, and (3) the Cle Elum River basin in Washington State. Two methods of climate scenario generation were compared: delta change and statistical downscaling. The delta change method uses differences between simulated current and future climate conditions from the Hadley Centre for Climate Prediction and Research (HadCM2) General Circulation Model (GCM) added to observed time series of climate variables. A statistical downscaling (SDS) model was developed for each basin using station data and output from the National Center for Environmental Prediction/National Center for Atmospheric Research (NCEPINCAR) reanalysis regridded to the scale of HadCM2. The SDS model was then used to simulate local climate variables using HadCM2 output for current and future conditions. Surface climate variables from each scenario were used in a precipitation‐runoff model. Results from this study show that, in the basins tested, a precipitation‐runoff model can simulate realistic runoff series for current conditions using statistically down‐scaled NCEP output. But, use of downscaled HadCM2 output for current or future climate assessments are questionable because the GCM does not produce accurate estimates of the surface variables needed for runoff in these regions. Given the uncertainties in the GCMs ability to simulate current conditions based on either the delta change or downscaling approaches, future climate assessments based on either of these approaches must be treated with caution.     

Effect of the Spatial Variability of Land Use,Soil Type,and Precipitation on Streamflows in Small Watersheds1

Abstract: The spatial variability of the data used in models includes the spatial discretization of the system into subsystems, the data resolution, and the spatial distribution of hydrologic features and parameters. In this study, we investigate the effect of the spatial distribution of land use, soil type, and precipitation on the simulated flows at the outlet of “small watersheds” (i.e., watersheds with times of concentration shorter than the model computational time step). The Soil and Water Assessment Tool model was used to estimate runoff and hydrographs. Different representations of the spatial data resulted in comparable model performances and even the use of uniform land use and soil type maps, instead of spatially distributed, was not noticeable. It was found that, although spatially distributed data help understand the characteristics of the watershed and provide valuable information to distributed hydrologic models, when the watershed is small, realistic representations of the spatial data do not necessarily improve the model performance. The results obtained from this study provide insights on the relevance of taking into account the spatial distribution of land use, soil type, and precipitation when modeling small watersheds.     

ROLE OF BASEFLOW IN RAINFALL-RUNOFF RELATIONS1

ABSTRACT. The role of initial baseflow, or the baseflow at the beginning of storm precipitation, in modifying mathematical rainfall-runoff relations is analyzed by using data from 95 storms over a drainage basin in Illinois. A regression model is set up with total runoff, surface runoff, baseflow runoff, and peak flow as dependent variables, and storm precipitation, initial baseflow, effective and total storm durations, and highest and lowest temperatures during the storm as independent variables. Stepwise regression analyses show that storm precipitation and initial baseflow are the most important variables for making dependent variable estimates. The standard error estimates using only storm precipitation and initial baseflow as predictors show a seasonal trend with a peak in July, August, or September. An understanding of the role of baseflow as an indicator of average soil moisture condition over the basin can be of great help in short-term reservoir regulation and flood warning.     

VERIFICATION OF THE RHEA-OROGRAPHIC-PRECIPITATION MODEL1

ABSTRACT: Observed April 1 snowpack accumulations within and near the Gunnison River basin in southwestern Colorado are compared with simulations from the Rhea-orographic-precipitation model to determine if the model simulates reliable magnitudes and temporal and spatial variability in winter precipitation for the basin. Twenty simulations of the Rhea model were performed using‘optimal’parameter sets determined for 10-kilometer (km) grids (10-km by 10-km grid cells) through stochastic calibration. Comparisons of Rhea-model simulations of winter precipitation with April 1 snowpack accumulations at 32 snowcourse stations were performed for the years 1972–1990. For most stations and most years the Rhea model reliably simulates the temporal and spatial variability in April 1 snowpack accumulations. However, in general, the Rhea-model underestimates April 1 snowpack accumulations in the Gunnison River basin area, and the underestimation is greatest for locations that receive the largest amount of snow. A significant portion of the error in Rhea-model simulations is due to the calibration of the Rhea model using gauge-catch precipitation measurements which can be as much as 50 percent below actual snowfall accumulations. Additional error in the Rhea-model simulations is a result of the comparison of gridded precipitation values to observed values measured at points.     

A COMPARISON OF HAND- AND SPLINE-DRAWN PRECIPITATION MAPS FOR MOUNTAINOUS MONTANA1

ABSTRACT: Average annual precipitation for the period 1961–1990 was estimated for a mountainous region in Montana with a Laplacian thin-plate spline (ANUSPLIN) and compared to a hand-drawn map. Input data included latitude, longitude, and elevation from a three-arc-second U.S. Geological Survey Digital Elevation Model of the Bozeman and Billings 1 × 2 topographic quadrangles and precipitation data at 96 stations. The two maps are similar in appearance. Digital comparison of the two maps with ARC/INFO's Grid tools shows that mean annual precipitation for the hand-contoured map is 22.9 inches and for the ANUSPLIN map is 23.7 inches. Of the 5,760,000 cells, 53 percent showed no difference between ANUSPLIN and hand-drawn maps; 19 percent showed a two-inch difference, and 28 percent showed more than 2 inches difference. Input data and model output at the same location are not different (standard deviation 1.77, p-value 0.76). Hand-drawn maps show two inches more precipitation during the 1961–90 period than during the 1941–1970 period. Similarly, measured data at 73 sites for the period 1961–1990 are on average 2.4 inches higher than the same stations during the 1941–1970 period. The difference is significant (p-value > 0.0001).     

ESTIMATED IMPACTS OF CLIMATE WARMING ON CALIFORNIA WATER AVAILABILITY UNDER TWELVE FUTURE CLIMATE SCENARIOS1

Spatially disaggregated estimates of over 131 stream‐flow, ground water, and reservoir evaporation monthly time series in California have been created for 12 different climate warming scenarios for a 72‐year period. Such disaggregated hydrologic estimates of multiple hydrologic cycle components are important for impact and adaptation studies of California's water system. A statewide trend of increased winter and spring runoff and decreased summer runoff is identified. Without operations modeling, approximate changes in water availability are estimated for each scenario. Even most scenarios with increased precipitation result in less available water because of the current storage systems' inability to catch increased winter streamflow in compensation for reduced summer runoff. The water availability changes are then compared with estimated changes in urban and agricultural water uses in California between now and 2100. The methods used in this study are relatively simple, but the results are qualitatively consistent with other studies focusing on the hydrologies of single basins or surface water alone.     

Illinois State Water Survey Bulletin 75: New Precipitation Frequency Study for Illinois

Storm frequency estimates and their temporal distributions are important in determining estimates of runoff or peak flow rates in many engineering and hydrological problems. Illinois State Water Survey Bulletin 70 has been serving as the design rainfall standard in Illinois since its publication in 1989. Although Bulletin 70 represented the best available data at the time of its publication, the standards needed to be reevaluated and updated after more than three decades and with the growing evidence of the nonstationary nature of heavy precipitation. The trends in heavy precipitation in Illinois prompted the creation of a new frequency study named Bulletin 75, providing precipitation frequencies for event durations ranging from 5 min to 10 days and for recurrence intervals ranging from 2 months to 500 years. The results are presented for the same 10 geographic sections in Illinois as in Bulletin 70 to maintain the continuity of hydrologic studies and compatibility with regulations. The primary goal of this paper is to outline some of the key methodological issues and challenges, to compare the results with the previous sources, and to highlight the effects of the changing precipitation standards on the development of amendments to existing ordinances. Lake County in Illinois, as one of the most affected urban areas with the highest change in heavy precipitation, was selected to illustrate issues related to the application of modified precipitation standards.     

Using a GIS Model to Identify Internally Drained Areas and Runoff Contribution in a Glaciated Watershed1

Macholl, Jacob A., Katherine A. Clancy, and Paul M. McGinley, 2011. Using a GIS Model to Identify Internally Drained Areas and Runoff Contribution in a Glaciated Watershed. Journal of the American Water Resources Association (JAWRA) 47(1):114‐125. DOI: 10.1111/j.1752‐1688.2010.00495.x Abstract: Glaciated watersheds are not easily delineated using geographic information systems’ elevation‐based algorithms, especially where stream networks are disconnected and there are large regions of internally drained areas. This paper presents the results of an analysis using the Potential Contributing Source Area (PCSA) model to identify potential contributing areas, defined as areas from which runoff is physically capable of reaching a drainage network. The investigation was conducted to define the potential contributing areas in a glaciated region of northwest Wisconsin. The curve number (CN) method was used to predict runoff volumes in the watershed. The streamflows of four tributaries were measured and the runoff portion of the hydrograph quantified to be compared with runoff estimates calculated using the potential contributing areas and the traditional catchment area. Runoff producing events occurred, but the use of area‐weighted CN values was unsuccessful in modeling runoff due to all precipitation depths during the study period falling below the initial abstraction. A distributed CN approach provided runoff estimates that were generally better using the potential contributing areas compared with using the traditional catchment area. The extent of the minimum contributing area, estimated for a range of precipitation events, was found to be substantially less than the potential contributing areas, suggesting that the PCSA model delimits the maximum boundary of potential contributing areas.     

WATER HARVESTING TREATMENT EVALUATION AT GRANITE REEF1

ABSTRACT: Yearly runoff efficiencies (total runoff/total precipitation), threshold retentions (precipitation needed to initiate runoff), and runoff-efficiencies-after-thresholds were determined for several water-harvesting catchment treatments at the Granite Reef test site. This information was found to be useful for showing (1) overall performance of catchments with time; (2) the distribution of the precipitation among runoff, surface retention, and infiltration; (3) why, how, and when certain treatments weathered and failed; (4) when to repair treatments; and (5) how to design catchments (size, site preparation, material selection, etc.). New impermeable membranes with smooth surfaces yielded nearly 100% of the precipitation. An asphalt-fiberglass treatment continues this high efficiency after 10 years of weathering – polyethylene after 8; efficiency of butyl sheeting was high initially but decreased rapidly after 9 years weathering. A standard roofing treatment retained up to 30% of the precipitation in the gravel covering. A concrete catchment lost as much as 50% of the total precipitation through micropores and surface cracks. Silicone treatments rapidly lost repellancy and efficiency, while paraffin treatments have weathered 5 years with little loss of efficiency.     

SPATIAL VARIABILITY IN WATER‐BALANCE MODEL PERFORMANCE IN THE CONTERMINOUS UNITED STATES1

ABSTRACT: A monthly water‐balance (WB) model was tested in 44 river basins from diverse physiographic and climatic regions across the conterminous United States (U.S.). The WB model includes the concepts of climatic water supply and climatic water demand, seasonality in climatic water supply and demand, and soil‐moisture storage. Exhaustive search techniques were employed to determine the optimal set of precipitation and temperature stations, and the optimal set of WB model parameters to use for each basin. It was found that the WB model worked best for basins with: (1) a mean elevation less than 450 meters or greater than 2000 meters, and/or (2) monthly runoff that is greater than 5 millimeters (mm) more than 80 percent of the time. In a separate analysis, a multiple linear regression (MLR) was computed using the adjusted R‐square values obtained by comparing measured and estimated monthly runoff of the original 44 river basins as the dependent variable, and combinations of various independent variables [streamflow gauge latitude, longitude, and elevation; basin area, the long‐term mean and standard deviation of annual precipitation; temperature and runoff; and low‐flow statistics (i.e., the percentage of months with monthly runoff that is less than 5 mm)]. Results from the MLR study showed that the reliability of a WB model for application in a specific region can be estimated from mean basin elevation and the percentage of months with gauged runoff less than 5 mm. The MLR equations were subsequently used to estimate adjusted R‐square values for 1,646 gauging stations across the conterminous U.S. Results of this study indicate that WB models can be used reliably to estimate monthly runoff in the eastern U.S., mountainous areas of the western U.S., and the Pacific Northwest. Applications of monthly WB models in the central U.S. can lead to uncertain estimates of runoff.     

Does the Temporal Resolution of Precipitation Input Influence the Simulated Hydrological Components Employing the SWAT Model?

This study aimed to evaluate the influence of sub‐daily precipitation time steps on model performance and hydrological components by applying the Green and Ampt infiltration method using the Soil and Water Assessment Tool (SWAT). Precipitation was measured at a resolution of 0.1 mm and aggregated to 5‐, 15‐, 30‐, and 60‐min time steps. Daily discharge data over a 10‐year period were used to calibrate and validate the model. Following a global sensitivity analysis, relevant parameters were optimized through an automatic calibration procedure using SWAT‐CUP for each time step. Daily performance statistics were almost equal among all four time steps (NSE ≈ 0.47). Discharge mainly consisted of groundwater flow (55%) and tile flow (42%), in reasonable proportions for the investigated catchment. In conclusion, model outputs were almost identical, showing simulations responded nearly independently of the chosen precipitation time step. This held true for (1) the selection of sensitive parameters, (2) performance statistics, (3) the shape of the hydrographs, and (4) flow components. However, a scenario analysis revealed that the precipitation time step becomes important when saturated hydraulic conductivities are low and curve numbers are high. The study suggests that there is no need in using precipitation time steps <1 h for lowland catchments dominated by soils with a low surface runoff potential if daily flow values are being considered. Editor's note : This paper is part of the featured series on SWAT Applications for Emerging Hydrologic and Water Quality Challenges. See the February 2017 issue for the introduction and background to the series.