KRIGED ESTIMATES OF TRANSMISSIVITY IN THE MESILLA BOLSON,NEW MEXICO1

ABSTRACT: Kriging utilizes a statistically based procedure of spatial interpolation that incorporates the spatial correlation structure of the phenomenon, and provides an error estimate. Kriging was applied to a total of 141 transmissivity values in an attempt to quantify the transmissivity distribution of the Santa Fe aquifer in Mesilla Bolson. New Mexico. The analysis produced contour maps of estimated transmissivity values and associated estimation variances. Through variogram analysis and fitting of an exponential variogrsm to 141 natural log of transmissivity (InT) values, the range was determined to be 3 miles, the average variance 2.74 (σ_InT= 1.65) with a mean of 8.65. Kriged estimates were generally lower when compared to estimates based on available transmissivity maps.     

Estimating Hydraulic Properties of Volcanic Aquifers Using Constant‐Rate and Variable‐Rate Aquifer Tests1

Abstract: In recent years the ground‐water demand of the population of the island of Maui, Hawaii, has significantly increased. To ensure prudent management of the ground‐water resources, an improved understanding of ground‐water flow systems is needed. At present, large‐scale estimations of aquifer properties are lacking for Maui. Seven analytical methods using constant‐rate and variable‐rate withdrawals for single wells provide an estimate of hydraulic conductivity and transmissivity for 103 wells in central Maui. Methods based on constant‐rate tests, although not widely used on Maui, offer reasonable estimates. Step‐drawdown tests, which are more abundantly used than other tests, provide similar estimates as constant‐rate tests. A numerical model validates the suitability of analytical solutions for step‐drawdown tests and additionally provides an estimate of storage parameters. The results show that hydraulic conductivity is log‐normally distributed and that for dike‐free volcanic rocks it ranges over several orders of magnitude from 1 to 2,500 m/d. The arithmetic mean, geometric mean, and median values of hydraulic conductivity are respectively 520, 280, and 370 m/d for basalt and 80, 50, and 30 m/d for sediment. A geostatistical approach using ordinary kriging yields a prediction of hydraulic conductivity on a larger scale. Overall, the results are in agreement with values published for other Hawaiian islands.     

HYDROGEOLOGY AND PATHS OF FLOW IN THE CARBONATE BEDROCK AQUIFER,NORTHWESTERN INDIANA1

ABSTRACT: The U.S. Geological Survey (USGS) is assessing the ground-water resources of the carbonate bedrock aquifers in Indiana and Ohio as part of their Regional Aquifer Systems Analysis program. Part of this assessment includes the determination of unknown aspects of the hydraulic characteristics, boundaries, and flow paths of the carbonate aquifer. To accomplish this, the USGS drilled three wells through the carbonate aquifer near the Kankakee River in northwestern Indiana. Geophysical logs were used to help determine depths and thicknesses for testing and to help describe geology at the three wells. Packer tests were used to determine direction of ground-water flow and to provide data for an analysis of the distribution of transmissivity in the carbonate aquifer. Transmissivity of the carbonates is associated with two physical characteristics of the rocks: fractures and interconnected porosity. Almost all of the transmissivity is derived from horizontal fracturing; however, only a few of the fractures present in the carbonate are transmissive. Some transmissivity is associated with a zone of fossiliferous, vuggy dolomite, which yields water from the rock matrix. Most of the transmissivity is associated with large fractures and solution crevices in the upper 30 feet of the bedrock; less transmissivity is associated with the deeper vuggy reef material, even where extensively fractured. Transmissivity of individual fractures and fossiliferous zones ranges from 300 to 27,000 feet squared per day. The aquifer bottom is defined by a lack of transmissive fractures and an increased shale content near the contact of the Silurian and Ordovician sections. Water-level data from the three wells indicate that flow is horizontal at well site 1 north of the Kankakee River, upward at well site 2 near the river, and downward at well site 3 south of the river. Most of the flow occurs in the upper part of the carbonate bedrock where fracturing and solution-enlarged crevices are most developed. Water levels indicate the the Kankakee River is a hydrologic boundary for the regional carbonate bedrock aquifer.     

Enhancing Hyrdological Simulation Program – FORTRAN Model Channel Hydraulic Representation1

Abstract: The hydrological simulation program – FORTRAN (HSPF) is a comprehensive watershed model that employs depth‐area‐volume‐flow relationships known as the hydraulic function table (FTABLE) to represent the hydraulic characteristics of stream channel cross‐sections and reservoirs. An accurate FTABLE determination for a stream cross‐section site requires an accurate determination of mean flow depth, mean flow width, roughness coefficient, longitudinal bed slope, and length of stream reach. A method that uses regional regression equations to estimate mean flow depth, mean flow width, and roughness coefficient is presented herein. FTABLES generated by the proposed method (Alternative Method) and FTABLES generated by Better Assessment Science Integrating Point and Nonpoint Sources (BASINS) were compared. As a result, the Alternative Method was judged to be an enhancement over the BASINS method. First, the Alternative Method employs a spatially variable roughness coefficient, whereas BASINS employs an arbitrarily selected spatially uniform roughness coefficient. Second, the Alternative Method uses mean flow width and mean flow depth estimated from regional regression equations whereas BASINS uses mean flow width and depth extracted from the National Hydrography Dataset (NHD). Third, the Alternative Method offers an option to use separate roughness coefficients for the in‐channel and floodplain sections of compound channels. Fourth, the Alternative Method has higher resolution in the sense that area, volume, and flow data are calculated at smaller depth intervals than the BASINS method. To test whether the Alternative Method enhances channel hydraulic representation over the BASINS method, comparisons of observed and simulated streamflow, flow velocity, and suspended sediment were made for four test watersheds. These comparisons revealed that the method used to estimate the FTABLE has little influence on hydrologic calibration, but greatly influences hydraulic and suspended sediment calibration. The hydrologic calibration results showed that observed versus simulated daily streamflow comparisons had Nash‐Sutcliffe efficiencies ranging from 0.50 to 0.61 and monthly comparisons had efficiencies ranging from 0.61 to 0.84. Comparisons of observed and simulated suspended sediments concentrations had model efficiencies ranging from 0.48 to 0.56 for the daily, and 0.28 to 0.70 for the monthly comparisons. The overall results of the hydrological, hydraulic, and suspended sediment concentration comparisons show that the Alternative Method yielded a relatively more accurate FTABLE than the BASINS method. This study concludes that hydraulic calibration enhances suspended sediment simulation performance, but even greater improvement in suspended sediment calibration can be achieved when hydrological simulation performance is improved. Any improvements in hydrological simulation performance are subject to improvements in the temporal and spatial representation of the precipitation data.     

USE OF HISTORIC FLOOD INFORMATION IN ESTIMATING FLOOD PEAKS ON UNGAGED WATERSHEDS1

Regional procedures to estimate flood magnitudes for ungaged watersheds typically ignore available site-specific historic flood information such as high water marks and the corresponding flow estimates, otherwise referred to as limited site-specific historic (LSSH) flood data. A procedure to construct flood frequency curves on the basis of LSSH flood observations is presented. Simple inverse variance weighting is employed to systematically combine flood estimates obtained from the LSSH data base with those from a regional procedure to obtain improved estimtes of flood peaks on the ungaged watershed. For the region studied, the variance weighted estimates of flow had a lower logarithmic standard error than either the regional or the LSSH flow estimates, when compared to the estimates determined by three standard distributions for gaged watersheds investigated in the development of the methodology. Use of the simple inverse variance weighting procedure is recommended when “reliable” estimates of LSSH floods for the ungaged site are available.     

MODELING GROUND WATER FLOW USING FLUX BOUNDARY CONDITIONS1

ABSTRACT: Determination of the boundary conditions for modeling ground water flow is a critical point especially in regional models. Normally the regional models require model areas that are greater than the given area of interest. This work focuses on the prediction of hydraulic heads in regional models using flux boundary conditions. The model uses flux boundary conditions that were estimated using a radial flow analog and Darcy's law. The regional model that is presented uses no parameter identification (inverse estimation) procedures. In the present work, the Houston area was used. The simulation of the hydrological conditions of the Chicot and Evangeline Aquifers that underlie the Houston area were made using the available information about the geological profile in the Houston region and the current information about the existing production wells. The regional model works as a forward problem. The system parameters such as hydraulic conductivity, specific storage, and hydrological stresses were specified, and the model predicts the hydraulic head. Actual data from piezometers operated by the U.S. Geological Survey (USGS) in many places throughout Houston were used as initial conditions. Some piezometric head data were generated using the regional variable theory called kriging to supply head estimates in areas where data were unavailable. The Modular Three Dimensional Finite Difference Groundwater Flow Model developed by the USGS was used to predict the hydraulic heads. The predicted ground water heads are compared to the actual data. The results show that the model performs well for locations where data were available.     

AN APPLICATION OF THE ANALYTIC ELEMENT METHOD IN MODELING FLORIDA EVERGLADES HYDROLOGY1

ABSTRACT: The large volumes of ground water that are discharged from the Everglades toward the Miami metropolitan area have historically posed a significant environmental water supply problem. In order to analyze the effects of seepage barriers on these subsurface outflows, the analytic element modeling code GFLOW was used to construct a ground water flow model of a region that includes a portion of the Everglades along with adjacent developed areas. The hydrology of this region can be characterized by a highly transmissive surficial aquifer in hydraulic contact with wetlands and canals. Calibration of the model to both wet and dry season conditions yielded satisfactory results, and it was concluded that the analytic element method is a suitable technique for modeling ground water flow in the Everglades environment. Finally, the model was used to evaluate the potential effectiveness of a subsurface barrier approximately two miles long for increasing water levels within the adjacent fringes of the Everglades National Park. It was found that the barrier had a negligible effect on water levels due to both its relatively short length and the high transmissivity of the surficial aquifer.     

BANKFULL DISCHARGE RECURRENCE INTERVALS AND REGIONAL HYDRAULIC GEOMETRY RELATIONSHIPS: PATTERNS IN THE PACIFIC NORTHWEST,USA1

ABSTRACT: The model bankfull discharge recurrence interval (annual series) (T_a) in streams has been approximated at a 1.5‐year flow event. This study tests the linkage between regional factors (climate, physiography, and ecoregion) and the frequency of bank‐full discharge events in the Pacific Northwest (PNW). Patterns of Ta were found to be significant when stratified by EPA Ecoregion. The mean value for T_a in the PNW is 1.4 years; however, when the data is stratified by ecoregion, the humid areas of western Oregon and Washington have a mean value of 1.2 years, while the dryer areas of Idaho and eastern Oregon and Washington have a mean value of 1.4 to 1.5 years. Among the four factors evaluated, vegetation association and average annual precipitation are the primary factors related to channel form and T_a. Based on the results of the T_a analyses, regional hydraulic geometry relationships of streams were developed for the PNW, which relate variables, such as bank‐full cross‐sectional area, width, depth, and velocity, to bankfull discharge and drainage area. The verification of T_a values, combined with the development of regional hydraulic geometry relationships, provides geographically relevant information that will result in more accurate estimates of hydraulic geometry variables in the PNW.     

ESTIMATING TRANSMISSIVITY AND HYDRAULIC CONDUCTIVITY OF CHICOT AQUIFER FROM SPECIFIC CAPACITY DATA1

ABSTRACT: Specific capacity data obtained from Well Construction reports which are available from USGS offices, can provide useful estimates of tranamissivity (T), and hydraulic conductivity (K), of an aquifer. The Chicot Aquifer in Louisiana is one of the largest sources of fresh ground water in North America. Hydrologic data collected for the Chicot Aquifer indicate that specific capacity tests can be used in estimating local and regional values for T and K, if the Cooper-Jacob equation for transient flow is used with proper corrections for well loss and partial penetration. Where full scale pumping test data are scarce, specific capacity test data that are adequately distributed spatially can be used to map changes in T and K values and can be summarized statistically to indicate applicable regional values. A computer program called “TGUESS” which is available from International Ground Water Modeling Center, Holcomb Research Institute, was used in this study. The contour maps for T and K values are prepared for different well depth intervals to avoid wide variation of values.     

MASS LOAD ESTIMATION ERRORS UTILIZING GRAB SAMPLING STRATEGIES IN A KARST WATERSHED1

ABSTRACT: Developing a mass load estimation method appropriate for a given stream and constituent is difficult due to inconsistencies in hydrologic and constituent characteristics. The difficulty may be increased in flashy flow conditions such as karst. Many projects undertaken are constrained by budget and manpower and do not have the luxury of sophisticated sampling strategies. The objectives of this study were to: (1) examine two grab sampling strategies with varying sampling intervals and determine the error in mass load estimates, and (2) determine the error that can be expected when a grab sample is collected at a time of day when the diurnal variation is most divergent from the daily mean. Results show grab sampling with continuous flow to be a viable data collection method for estimating mass load in the study watershed. Comparing weekly, biweekly, and monthly grab sampling, monthly sampling produces the best results with this method. However, the time of day the sample is collected is important. Failure to account for diurnal variability when collecting a grab sample may produce unacceptable error in mass load estimates. The best time to collect a sample is when the diurnal cycle is nearest the daily mean.     

USING SENSITIVITY ANALYSIS TO ASSIST PARAMETER ZONATION IN GROUND WATER FLOW MODEL1

ABSTRACT: For numerical modeling of ground water movement in a real aquifer system, the aquifer is usually divided into hydrogeologically defined zones, each with its own parameter values. The responses of the system, such as head or drawdown, are often available only in some of the zones. The estimated parameters of all the zones are based on the measured response in these limited zones. However, the estimates for some of the zones may be very uncertain, and these zones are therefore not justified by the data. In this paper, an approach is presented to understand which zone may produce uncertain parameter values and should be lumped with its neighbor. This approach is demonstrated using a regional numerical model for pumping test analysis in the Nottinghamshire aquifer, UK. A step-by-step process is used in identifying the aquifer zones and estimating their parameters based on the principle of using the smallest possible numbers of zones and parameters for adequate representation of the drawdown response. After the parameters of each zone are estimated, the sensitivity features of these parameters are examined. The results show that the parameters in one zone can be estimated properly by the drawdown in another zone only when there is significant sensitivity. For transmissivity, sensitivity between zones occurs when there is significant flow between them. For storativity, sufficient sensitivity can occur without large flows between the zones, provided that one zone causes significant drawdown in the other. This idea can be extended to the flow model for a large aquifer system. If the aquifer is divided in such a way that aquifer responses are not sensitive to the parameters in some of the zones, the parameters in those zones cannot be estimated properly and should be lumped into their neighboring zones. In this way, a simple but more reasonable model can be built.     

UNCERTAINTIES IN ESTIMATING THE WATER BALANCE OF LAKES1

ABSTRACT: Evaluation of hydrologic methodology used in a number of water balance studies of lakes in the United States shows that most of these studies calculate one or more terms of the budget as the residual. A literature review was made of studies in which the primary purpose was error analysis of hydrologic measurement and interpretation. Estimates of precipitation can have a wide range of error, depending on the gage placement, gage spacing, and areal averaging technique. Errors in measurement of individual storms can be as high as 75 percent. Errors in short term averages are commonly in the 15-30 percent range, but decrease to about 5 percent or less for annual estimates. Errors in estimates of evaporation can also vary widely depending on instrumentation and methodology. The energy budget is the most accurate method of calculating evaporation; errors are in the 10–15 percent range. If pans are used that are located a distance from the lake of interest, errors can be considerable. Annual pan-to-lake coefficients should not be used for monthly estimates of evaporation because they differ from the commonly used coefficient of 0.7 by more than 100 percent. Errors in estimates of stream discharge are often considered to be within 5 percent. If the measuring section, type of flow profile, and other considerations, such as stage discharge relationship, are less than ideal errors in estimates of stream discharge can be considerably greater than 5 percent. Errors in estimating overland (nonchannelized) flow have not been evaluated, and in most lake studies this component is not mentioned. Comparison of several lake water balances in which the risdual consists solely of errors in measurement, shows that such a residual, if interpreted as ground water, can differ from an independent estimate of ground water by more than 100 percent.     

EVALUATION OF THE MAXEY-EAKIN METHOD FOR ESTIMATING RECHARGE TO GROUND-WATER BASINS IN NEVADA1

AESTRACT An evaluation of the Maxey-Eakin method for calculating recharge to ground-water basins in Nevada was performed. The evaluation consisted of comparing Maxey-Eakin estimates with independent estimates of recharge, and analyzing the nature of the differences between the groups of estimates. In the comparison with the Maxey-Eakin estimates, two different groups of independent estimates were used: (1) 40 recharge estimates that were identified from water budgets contained in reports by the Nevada Department of Conservation and Natural Resources and (2) 27 recharge estimates that were identified from previous studies that used models. The results of the comparisons indicate generally good agreement between the Maxey-Eakin estimates and both groups of independent estimates. To quantify this agreement, an analysis was conducted to estimate the uncertainty in the Maxey-Eakin method. The analysis produced an upper bound on the standard deviation of the Maxey-Eakin estimate for a given basin. For the group of 40 water-budget estimates, the upper bound on the standard deviation for an individual basin is 4,800 acre-ft/yr, and the corresponding coefficient of variation of the Maxey-Eakin estimate is no greater than 44 percent. For the group of 27 model estimates, the upper bound on the standard deviation is 4,100 acre-ft/yr, and the corresponding coefficient of variation is no greater than 24 percent.     

Assessment of Coarse Sediment Mobility in the Black Canyon of the Gunnison River,Colorado

The Gunnison River in the Black Canyon of the Gunnison National Park (BCNP) near Montrose, Colorado is a mixed gravel and bedrock river with ephemeral side tributaries. Flow rates are controlled immediately upstream by a diversion tunnel and three reservoirs. The management of the hydraulic control structures has decreased low-frequency, high-stage flows, which are the dominant geomorphic force in bedrock channel systems. We developed a simple model to estimate the extent of sediment mobilization at a given flow in the BCNP and to evaluate changes in the extent and frequency of sediment mobilization for flow regimes before and after flow regulation in 1966. Our methodology provides a screening process for identifying and prioritizing areas in terms of sediment mobility criteria when more precise systematic field data are unavailable. The model uses the ratio between reach-averaged bed shear stress and critical shear stress to estimate when a particular grain size is mobilized for a given reach. We used aerial photography from 1992, digital elevation models, and field surveys to identify individual reaches and estimate reach-averaged hydraulic geometry. Pebble counts of talus and debris fan deposits were used to estimate regional colluvial grain-size distributions. Our results show that the frequency of flows mobilizing river bank sediment along a majority of the Gunnison River in the BCNP has significantly declined since 1966. The model results correspond well to those obtained from more detailed, site-specific field studies carried out by other investigators. Decreases in the frequency of significant sediment-mobilizing flows were more pronounced for regions within the BCNP where the channel gradient is lower. Implications of these results for management include increased risk of encroachment of vegetation on the active channel and long-term channel narrowing by colluvial deposits. It must be recognized that our methodology represents a screening of regional differences in sediment mobility. More precise estimates of hydraulic and sediment parameters would likely be required for dictating quantitative management objectives within the context of sediment mobility and sensitivity to changes in the flow regime.     

GEOSTATISTICAL ESTIMATION OF HORIZONTAL HYDRAULIC CONDUCTIVITY FOR THE KIRKWOOD‐COHANSEY AQUIFER1

ABSTRACT: The Kirkwood‐Cohansey aquifer has been identified as a critical source for meeting existing and expected water supply needs for southern New Jersey. Several contaminated sites exist in the region; their impact on the aquifer has to be evaluated using ground water flow and transport models. Ground water modeling depends on availability of measured hydrogeologic data (e.g., hydraulic conductivity, for parameterization of the modeling runs). However, field measurements of such critical data have inadequate spatial density, and their locations are often clustered. The goal of this study was to research, compile, and geocode existing data, then use geostatistics and advanced mapping methods to develop a map of horizontal hydraulic conductivity for the Kirkwood‐Cohansey aquifer. Spatial interpolation of horizontal hydraulic conductivity measurements was performed using the Bayesian Maximum Entropy (BME) Method implemented in the BMELib code library. This involved the integration of actual measurements with soft information on likely ranges of hydraulic conductivity at a given location to obtain estimate maps. The estimation error variance maps provide an insight into the uncertainty associated with the estimates, and indicate areas where more information on hydraulic conductivity is required.     

Calibrating a Basin‐Scale Groundwater Model to Remotely Sensed Estimates of Groundwater Evapotranspiration

Remotely sensed vegetation indices correspond to canopy vigor and cover and have been successfully used to estimate groundwater evapotranspiration (ET_g) over large spatial and temporal scales. However, these data do not provide information on depth to groundwater (dtgw) necessary for groundwater models (GWM) to calculate ET_g. An iterative approach is provided that calibrates GWM to ET_g derived from Landsat estimates of the Enhanced Vegetation Index (EVI). The approach is applied to different vegetation groups in Mason Valley, Nevada over an 11‐year time span. An uncertainty analysis is done to estimate the resulting mean and 90% confidence intervals in ET_g to dtgw relationships to quantify errors associated with plant physiologic complexity, species variability, and parameter smoothing to the 100 m GWM‐grid, temporal variability in soil moisture and nonuniqueness in the solution. Additionally, a first‐order second moment analysis shows ET_g to dtgw relationships are almost exclusively sensitive to estimated land surface, or maximum, ET_g despite relatively large uncertainty in extinction depths and hydraulic conductivity. The EVI method of estimating ET_g appears to bias ET_g during years with exceptionally wet spring/summer conditions. Excluding these years improves model performance significantly but highlights the need to develop a methodology that accounts not only on quantity but timing of annual precipitation on phreatophyte greenness.     

ECONOMIC RESOURCE MULTIPLIERS FOR REGIONAL IMPACT ANALYSIS1

ABSTRACT: Tradeoffs between regional economic development and resource use is a question often confronting local decisionmakers. A resource-interindustry model can be used to depict the interrelationships between regional economic sectors as to household income and/or employment and resource use. A resource-interindustry model was developed for Humboldt and Lander Counties in Nevada which shows the tradeoffs between regional household income (wages, salaries, profits, and rents) and/or employment and water usage. Water income and water employment multipliers can be ranked, enabling decisionmakers to realize sectors which require greatest regional water usage to regional household income and/or employment.     

ASSESSMENT OF HYDRAULIC PROPERTIES IN AN UNCONFINED ALLUVIAL AQUIFER NEAR GRAND ISLAND,NEBRASKA1

ABSTRACT: A method is presented to analyze time-drawdown data from one or more observation wells for the calculation of four hydraulic parameters for unconfined aquifers: vertical hydraulic conductivity, horizontal hydraulic conductivity, storage coefficient, and specific yield. The hydraulic parameter results are further analyzed for reliability and the possible ranges of the actual parameter values. After verification using a theoretical example, the method was used to analyze pumping test data from 22 observation wells in an unconfined alluvial aquifer near Grand Island, Nebraska. Results indicate that this method can be used to efficiently calculate the four hydraulic parameters in this type of aquifer. The method can also identify the impact of measurement errors on the parameter estimates, and provide ranges of the actual parameter values. The parameter values calculated using this method were compared to those determined using other theories. It is found that this method is very useful for calculating the hydraulic properties from pumping test data and for analyzing the parameter reliability.     

Management of the Spring Snowmelt Recession in Regulated Systems

In an effort to restore predictable ecologically relevant spring snowmelt recession flow patterns in rivers regulated by dams, this study defined a methodology by which spring flow regimes can be modeled in regulated systems from the quantifiable characteristics of spring snowmelt recessions in unregulated rivers. An analysis of eight unregulated rivers across the Sierra Nevada mountain range in California found that unregulated systems behaved similarly with respect to seasonal spring patterns and recession limb curvature, and thus prescribed flows could be designed in a manner that mimics those predictable characteristics. Using the methodology to quantify spring recession flows in terms of a daily percent decrease in flow, a series of flow recession scenarios were created for application in an existing hydrodynamic model for the regulated Rubicon River. The modeling results showed that flow recessions with slow ramping rates similar to those observed in unregulated rivers (less than 10% per day) were likely to be protective of native aquatic species, such as the Foothill yellow‐legged frog, while flows that receded at greater rates would likely result in desiccation of egg masses and potential stranding of tadpoles and fry. Furthermore, recession rates of less than 10% per day provided the most spatially diverse hydraulic habitat in the modeled domain for an appropriate duration in spring to support all native species guilds and maximize aquatic biodiversity.     

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.