USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Hydrology”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 361 records · Page 20Linked to original sources

Fast ground-water mixing and basal recharge in an unconfined, alluvial aquifer, Konza LTER Site, Northeastern Kansas

Ground-water chemistry and water levels at three levels in a well nest were monitored biweekly for two and a half years in a shallow unconfined floodplain aquifer in order to study the dynamics of such shallow aquifers. The aquifer, in northeastern Kansas, consists of high porosity, low hydraulic conductivity fine-grained sediments dominated by silt and bounded by fractured limestone and shale bedrock. Results show that the aquifer underwent chemical stratification followed by homogenization three times during the study period. The length of time between maximum stratification and complete homogenization was 3-5 months. The chemical parameters most useful for demonstrating the mixing trends were dissolved nitrate and sulfate. Higher nitrate concentrations were typical of unsaturated zone water and were sourced from fertilizer applied to the cultivated fields on the floodplain. Variations in sulfate concentrations are attributed to dissolution of rare gypsum in limestone bedrock and variable evapoconcentration in the unsaturated zone. The mixing of three chemically different waters (entrained, unsaturated-zone water; water entering the base of the floodplain aquifer; and water in residence before each mixing event) was simulated. The resident water component for each mixing event was a fixed composition based on measured water chemistry in the intermediate part of the aquifer. The entrained water composition was calculated using a measured composition of the shallow part of the aquifer and measurements of soil-water content in the unsaturated zone. The incoming basal water composition and the fractions of each mixing component were fitted to match the measured chemistry at the three levels in the aquifer. A conceptual model for this site explains: (1) rapid water-level rises, (2) water-chemistry changes at all levels in the aquifer coincident with the water-level rises, (3) low measured hydraulic conductivity of the valley fill and apparent lack of preferential flow pathways, (4) minuscule amounts of unsaturated-zone recharge, and (5) dissolved oxygen peaks in the saturated zone lagging water-level peaks. We postulate that rainfall enters fractures in bedrock adjacent to the floodplain. This recharge water moves rapidly through the fractured bedrock into the base of the floodplain aquifer. The recharge event through the bedrock causes a rapid rise in water level in the floodplain aquifer, and the chemistry of the deepest water in the floodplain aquifer changes at that time. The rising water also entrains slow-moving, nitrate-rich, unsaturated-zone water, altering the chemistry of water in the shallow part of the aquifer. Vertical chemical stratification in the aquifer is thus created by the change in water chemistry in the upper and lower parts of the saturated zone. As the water level begins to decline, the aquifer undergoes mixing that eventually results in homogeneous water chemistry. The rise in water level from the recharge event also displaces gas from the unsaturated zone that is then replaced as the water level declines following the recharge event. This new, oxygen-rich vadose-zone air equilibrates rapidly with saturated-zone water, resulting in a dissolved oxygen pulse in the ground water that peaks one-half to 2 months after the water-level peak. This oxygen pulse subsequently declines over a period of 2-6 months. ?? 2003 Elsevier B.V. All rights reserved.

Journal of Hydrology↗

Water table fluctuations near an incised stream, Walnut Creek, Iowa

Incised channels are common features in many agricultural watersheds, but the effects of channel incision on riparian water table conditions have been poorly documented. In this study, we evaluate the water table fluctuations in the floodplain near an incised stream (Walnut Creek, Iowa) and investigate the roles that channel incision and variable recharge play in modifying the water table configuration in the floodplain. Groundwater flows from higher landscape positions towards Walnut Creek under hydraulic gradients that were steepest near the upland/floodplain contact and in the near-stream riparian zone. Annually, water table fluctuations on the floodplain were greatest in wells located 30 m from the creek, midway between the creek and upland. Water levels monitored continuously during a runoff event indicated that bank storage was confined to a narrow zone adjacent to the channel. A steady-state, one-dimensional analytical model was developed to describe the shape of the water table surface near an incised stream and evaluate how variable groundwater recharge and channel bed lowering has affected the shape of the water table surface. Results from this study have implications for managing the riparian buffers of incised streams with successful establishment dependent upon matching buffer vegetation to riparian water table conditions. ?? 2003 Elsevier B.V. All rights reserved.

Journal of Hydrology↗

The rotating movement of three immiscible fluids - A benchmark problem

A benchmark problem involving the rotating movement of three immiscible fluids is proposed for verifying the density-dependent flow component of groundwater flow codes. The problem consists of a two-dimensional strip in the vertical plane filled with three fluids of different densities separated by interfaces. Initially, the interfaces between the fluids make a 45??angle with the horizontal. Over time, the fluids rotate to the stable position whereby the interfaces are horizontal; all flow is caused by density differences. Two cases of the problem are presented, one resulting in a symmetric flow field and one resulting in an asymmetric flow field. An exact analytical solution for the initial flow field is presented by application of the vortex theory and complex variables. Numerical results are obtained using three variable-density groundwater flow codes (SWI, MOCDENS3D, and SEAWAT). Initial horizontal velocities of the interfaces, as simulated by the three codes, compare well with the exact solution. The three codes are used to simulate the positions of the interfaces at two times; the three codes produce nearly identical results. The agreement between the results is evidence that the specific rotational behavior predicted by the models is correct. It also shows that the proposed problem may be used to benchmark variable-density codes. It is concluded that the three models can be used to model accurately the movement of interfaces between immiscible fluids, and have little or no numerical dispersion. ?? 2003 Elsevier B.V. All rights reserved.

Journal of Hydrology↗

Use of chemical and isotopic tracers to assess nitrate contamination and ground-water age, Woodville Karst Plain, USA

Concerns regarding ground-water contamination in the Woodville Karst Plain have arisen due to a steady increase in nitrate-N concentrations (0.25-0.90 mg/l) during the past 30 years in Wakulla Springs, a large regional discharge point for water (9.6 m3/s) from the Upper Floridan aquifer (UFA). Multiple isotopic and chemical tracers were used with geochemical and lumped-parameter models (exponential mixing (EM), dispersion, and combined exponential piston flow) to assess: (1) the sources and extent of nitrate contamination of ground water and springs, and (2) mean transit times (ages) of ground water. Delta 15N-NO3 values (1.7-13.8???) indicated that nitrate in ground water originated from localized sources of inorganic fertilizer and human/animal wastes. Nitrate in spring waters (??15N-NO3=5.3-8.9???) originated from both inorganic and organic N sources. Nitrate-N concentrations (<0.02-16 mg/l) were highly variable both spatially and vertically in the oxic UFA, based on water samples from 46 wells and four springs collected from 1997 to 2000. During high-flow conditions, spring waters had decreased nitrate and increased DOC concentrations that resulted from mixtures of 20-95% surface water. Although higher nitrate-N concentrations (>1.0 mg/l) were associated with shallow wells (open intervals less than 15 m below land surface), elevated nitrate concentrations in deeper wells are consistent with mixtures of water from shallow and deep zones in the UFA as indicated from geochemical mixing models and the distribution of mean transit times (5-90 years) estimated using lumped-parameter flow models. Ground water with mean transit times of 10 years or less tended to have higher dissolved organic carbon concentrations, lower dissolved solids, and lower calcite saturation indices than older waters, indicating mixing with nearby surface water that directly recharges the aquifer through sinkholes. Significantly higher values of pH, magnesium, dolomite saturation index, and phosphate in springs and deep water (>45 m) relative to a shallow zone (<45 m) were associated with longer ground-water transit times (50-90 years). Chemical differences with depth in the aquifer result from deep regional flow of water recharged through low permeability sediments (clays and clayey sands of the Hawthorn Formation) that overlie the UFA upgradient from the karst plain.

Journal of Hydrology↗

Nitrate and herbicide loading in two groundwater basins of Illinois' sinkhole plain

This investigation was designed to estimate the mass loading of nitrate (NO3-) and herbicides in spring water discharging from groundwater basins in an agriculturally dominated, mantled karst terrain. The loading was normalized to land use and NO3- and herbicide losses were compared to estimated losses in other agricultural areas of the Midwestern USA. Our study area consisted of two large karst springs that drain two adjoining groundwater basins (total area of 37.7 km2) in southwestern Illinois' sinkhole plain, USA. The springs and stream that they form were monitored for almost 2 years. Nitrate-nitrogen (NO3-N) concentrations at three monitoring sites were almost always above the background concentration (1.9 mg/l). NO3-N concentrations at the two springs ranged from 1.08 to 6.08 with a median concentration of 3.61 mg/l. Atrazine and alachlor concentrations ranged from <0.01 to 34 ??g/l and <0.01 to 0.98 ??g/l, respectively, with median concentrations of 0.48 and 0.12 ??g/l, respectively. Approximately 100,000 kg/yr of NO3-N, 39 kg/yr of atrazine, and 2.8 kg/yr of alachlor were discharged from the two springs. Slightly more than half of the discharged NO3- came from background sources and most of the remainder probably came from fertilizer. This represents a 21-31% loss of fertilizer N from the groundwater basins. The pesticide losses were 3.8-5.8% of the applied atrazine, and 0.05-0.08% of the applied alachlor. The loss of atrazine adsorbed to the suspended solid fraction was about 2 kg/yr, only about 5% of the total mass of atrazine discharged from the springs. ?? 2004 Elsevier B.V. All rights reserved.

Journal of Hydrology↗

Identifying areas of basin-floor recharge in the Trans-Pecos region and the link to vegetation

Comparative water potential and chloride profiles (∼10 m deep) collected from four vegetation communities in the Trans-Pecos region of the Chihuahuan Desert were assessed to evaluate the potential for using vegetation patterns as a means of efficiently improving large-scale estimates of basin-floor recharge in semiarid and arid regions. Analytical solutions and multiphase flow and transport modeling constrained flux histories and current fluxes across the water table at each site. Chloride bulge profiles containing ∼12–15 kyr of atmospheric deposition and long-term drying water potential profiles typified most desertscrub and grassland sites. In contrast, evidence of episodic sub-root zone percolation and chloride profiles containing <250 yr of atmospheric deposition characterized the woodland site. The results suggested that the desertscrub and grassland areas support small upward fluxes across the water table (nonrecharge), whereas the woodland site supports significant downward fluxes across the water table (recharge). A nonrecharge–recharge transition was identified to be collocated with a grassland–woodland ecotone. The establishment of vegetation–recharge relationships such as this will improve estimates of basin-scale recharge by identifying regions where no recharge is expected and regions where recharge is expected and point measurements should be concentrated. An approach integrating remotely sensed spatial distributions of vegetation and indicator relationships to recharge is both timely and warranted, although several caveats, as revealed in this study, should be noted. For example, the relative importance and distribution of vertical conduits that permit percolation to the water table merits future investigation.

Journal of Hydrology↗

Comparison of ground-water flow model particle-tracking results and isotopic data in the Mojave River ground-water basin, southern California, USA

Flow-path and time-of-travel results for the Mojave River ground-water basin, southern California, calculated using the ground-water flow model MODFLOW and particle-tracking model MODPATH were similar to flow path and time-of-travel interpretations derived from delta-deuterium and carbon-14 data. Model and isotopic data both show short flow paths and young ground-water ages throughout the floodplain aquifer along most the Mojave River. Longer flow paths and older ground-water ages as great as 10,000 years before present were measured and simulated in the floodplain aquifer near the Mojave Valley. Model and isotopic data also show movement of water between the floodplain and regional aquifer and subsequent discharge of water from the river to dry lakes in some areas. It was not possible to simulate the isotopic composition of ground-water in the regional aquifer away from the front of the San Gabriel and San Bernardino Mountains - because recharge in these areas does not occur under the present-day climatic conditions used for calibration of the model.

Journal of Hydrology↗

Aquifer vulnerability to pesticide pollution: Combining soil, land-use and aquifer properties with molecular descriptors

This study uses an extensive survey of herbicides in groundwater across the midwest United States to predict occurrences of a range of compounds across the region from a combination of their molecular properties and the properties of the catchment of a borehole. The study covers 100 boreholes and eight pesticides. For each of the boreholes its catchment the soil, land-use and aquifer properties were characterized. Discriminating boreholes where pollution occurred from those where no pollution occurred gave a model that was 74% correct with organic carbon content, percentage sand content and depth to the water table being significant properties of the borehole catchment. Molecular topological descriptors as well as K oc , solubility and half-life were used to characterize each compound included in the study. Inclusion of molecular properties makes it possible to discriminate between occurrence and non-occurrence of each compound in each well. The best-fit model combines: organic carbon content, percentage sand content and depth to the water table with molecular descriptors representing molecular size, molecular branching and functional group composition of the herbicides.

Midwest region↗

Baseflow contribution to nitrate-nitrogen export from a large, agricultural watershed, USA

Nitrate-nitrogen export from the Raccoon River watershed in west-central Iowa is among the highest in the United State and contributes to impairment of downstream water quality. We examined a rare long-term record of streamflow and nitrate concentration data (1972-2000) to evaluate annual and seasonal patterns of nitrate losses in streamflow and baseflow from the Raccoon River. Combining hydrograph separation with a load estimation program, we estimated that baseflow contributes approximately two-thirds (17.3 kg/ha) of the mean annual nitrate export (26.1 kg/ha). Baseflow transport was greatest in spring and late fall when baseflow contributed more than 80% of the total export. Herein we propose a 'baseflow enrichment ratio' (BER) to describe the relation of baseflow water with baseflow nitrate loads. The long-term ratio of 1.23 for the Raccoon River suggests preferential leaching of nitrate to baseflow. Seasonal patterns of the BER identified the strong link between the baseflow nitrate loads and seasonal crop nitrogen requirements. Study results demonstrate the utility of assessing the baseflow contribution to nitrate loads to identify appropriate control strategies for reducing baseflow delivery of nitrate. ?? 2004 Elsevier B.V. All rights reserved.

Journal of Hydrology↗

Linear model describing three components of flow in karst aquifers using 18O data

The stable isotope of oxygen, 18O, is used as a naturally occurring ground-water tracer. Time-series data for ??18O are analyzed to model the distinct responses and relative proportions of the conduit, intermediate, and diffuse flow components in karst aquifers. This analysis also describes mathematically the dynamics of the transient fluid interchange between conduits and diffusive networks. Conduit and intermediate flow are described by linear-systems methods, whereas diffuse flow is described by mass-balance methods. An automated optimization process estimates parameters of lognormal, Pearson type III, and gamma distributions, which are used as transfer functions in linear-systems analysis. Diffuse flow and mixing parameters also are estimated by these optimization methods. Results indicate the relative proximity of a well to a main conduit flowpath and can help to predict the movement and residence times of potential contaminants. The three-component linear model is applied to five wells, which respond to changes in the isotopic composition of point recharge water from a sinking stream in the Madison aquifer in the Black Hills of South Dakota. Flow velocities as much as 540 m/d and system memories of as much as 71 years are estimated by this method. Also, the mean, median, and standard deviation of traveltimes; time to peak response; and the relative fraction of flow for each of the three components are determined for these wells. This analysis infers that flow may branch apart and rejoin as a result of an anastomotic (or channeled) karst network.

Journal of Hydrology↗

Investigating surface water-well interaction using stable isotope ratios of water

Because surface water can be a source of undesirable water quality in a drinking water well , an understanding of the amount of surface water and its travel time to the well is needed to assess a well's vulnerability. Stable isotope ratios of oxygen in river water at the City of La Crosse, Wisconsin, show peak-to-peak seasonal variation greater than 4‰ in 2001 and 2002. This seasonal signal was identified in 7 of 13 city municipal wells, indicating that these 7 wells have appreciable surface water contributions and are potentially vulnerable to contaminants in the surface water. When looking at wells with more than 6 sampling events, a larger variation in δ 18 O compositions correlated with a larger fraction of surface water, suggesting that samples collected for oxygen isotopic composition over time may be useful for identifying the vulnerability to surface water influence even if a local meteoric water line is not available. A time series of δ 18 O from one of the municipal wells and from a piezometerlocated between the river and the municipal well showed that the travel time of flood water to the municipal well was approximately 2 months; non-flood arrival times were on the order of 9 months. Four independent methods were also used to assess time of travel. Three methods (groundwater temperature arrival times at the intermediate piezometer, virus-culture results, and particle tracking using a numerical groundwater-flow model) yielded flood and non-flood travel times of less than 1 year for this site. Age dating of one groundwater sample using 3 H– 3 He methods estimated an age longer than 1 year, but was likely confounded by deviations from piston flow as noted by others. Chlorofluorocarbons and SF 6 analyses were not useful at this site due to degradation and contamination, respectively. This work illustrates the utility of stable hydrogen and oxygen isotope ratios of water to determine the contribution and travel time of surface water in groundwater, and demonstrates the importance of using multiple methods to improve estimates for time of travel of 1 year or less.

Journal of Hydrology↗

Utility of Penman-Monteith, Priestley-Taylor, reference evapotranspiration, and pan evaporation methods to estimate pasture evapotranspiration

Actual evapotranspiration (ETa) was measured at 30-min resolution over a 19-month period (September 28, 2000-April 23, 2002) from a nonirrigated pasture site in Florida, USA, using eddy correlation methods. The relative magnitude of measured ETa (about 66% of long-term annual precipitation at the study site) indicates the importance of accurate ET a estimates for water resources planning. The time and cost associated with direct measurements of ETa and the rarity of historical measurements of ETa make the use of methods relying on more easily obtainable data desirable. Several such methods (Penman-Monteith (PM), modified Priestley-Taylor (PT), reference evapotranspiration (ET 0), and pan evaporation (Ep)) were related to measured ETa using regression methods to estimate PM bulk surface conductance, PT ??, ET0 vegetation coefficient, and Ep pan coefficient. The PT method, where the PT ?? is a function of green-leaf area index (LAI) and solar radiation, provided the best relation with ET a (standard error (SE) for daily ETa of 0.11 mm). The PM method, in which the bulk surface conductance was a function of net radiation and vapor-pressure deficit, was slightly less effective (SE=0.15 mm) than the PT method. Vegetation coefficients for the ET0 method (SE=0.29 mm) were found to be a simple function of LAI. Pan coefficients for the Ep method (SE=0.40 mm) were found to be a function of LAI and Ep. Historical or future meteorological, LAI, and pan evaporation data from the study site could be used, along with the relations developed within this study, to provide estimates of ETa in the absence of direct measurements of ETa. Additionally, relations among PM, PT, and ET0 methods and ETa can provide estimates of ETa in other, environmentally similar, pasture settings for which meteorological and LAI data can be obtained or estimated. ?? 2004 Elsevier B.V. All rights reserved.

Journal of Hydrology↗

Effects of suburban development on runoff generation in the Croton River basin, New York, USA

The effects of impervious area, septic leach-field effluent, and a riparian wetland on runoff generation were studied in three small (0.38–0.56 km 2 ) headwater catchments that represent a range of suburban development (high density residential, medium density residential, and undeveloped) within the Croton River basin, 70 km north of New York City. Precipitation, stream discharge, and groundwater levels were monitored at 10–30 min intervals for 1 year, and stream water and groundwater samples were collected biweekly for δ 18 O, NO 3 − , and SO 4 2− analysis for more than 2 years during an overlapping period in 2000–2002. Data from 27 storms confirmed that peak magnitudes increased and recession time decreased with increasing development, but lags in peak arrival and peak discharge/mean discharge were greatest in the medium density residential catchment, which contains a wetland in which storm runoff is retained before entering the stream. Baseflow during a dry period from Aug. 2001–Feb. 2002 was greatest in the high-density residential catchment, presumably from the discharge of septic effluent through the shallow groundwater system and into the stream. In contrast, moderate flows during a wet period from Mar.–Aug. 2002 were greatest in the undeveloped catchment, possibly as a result of greater subsurface storage or greater hydraulic conductivity at this site. The mean residence time of baseflow was about 30 weeks at all three catchments, indicating that human influence was insufficient to greatly affect the groundwater recharge and discharge properties that determine catchment residence time. These results suggest that while suburban development and its associated impervious surfaces and storm drains accelerate the transport of storm runoff into streams, the combined effects of remnant natural landscape features such as wetlands and human alterations such as deep groundwater supply and septic systems can change the expected effects of human development on storm runoff and groundwater recharge.

New York↗

Coupled inverse modeling of vadose zone water, heat, and solute transport: Calibration constraints, parameter nonuniqueness, and predictive uncertainty

In this study, an inverse methodology is presented and used to evaluate the effect that calibration of a synthetic artificial recharge model, constrained by different combinations of measurements (pressure head, temperature, and concentration), has on estimated vadose zone model parameter-value nonuniqueness and predictive water, heat, and solute transport uncertainty. Several findings are arrived at following model calibration and predictive analysis. First, composite scaled sensitivities revealed that all calibration measurement combinations contributed to the estimation of 30 water, heat, and solute transport parameters by inverting a set of vadose zone transport equations that were coupled explicitly through dependent variables and implicitly through parameters and fluid properties. Second, despite excellent model quality and perfect match of simulated-to-measured dependent field variables, the limitations in information content of field measurements used to constrain the calibration process promoted correlation among parameters; correlation among parameters promoted parameter nonuniqueness; and parameter nonuniqueness promoted predictive uncertainty. Consequently, simulations by transport models calibrated against field information represent a single realization associated with some quantifiable range of predictive uncertainty. Third, a primary reduction in uncertainty was achieved by increasing the number of calibration-constraint measurements, but reductions in uncertainty appeared restricted implying a practical limit to parameterization detail. Fourth, for a fixed number of measurements, a less prominent reduction in the range of predictive uncertainty could be realized through selective use of measurement types to constrain the calibration process. Therefore, field measurement types used to constrain the calibration process should be matched to target predictions. Fifth, because correlation among parameters contributes to predictive uncertainty, it may be possible to further reduce predictive uncertainty by estimating parameters that also minimize the largest eigenvalue in the normalized eigenvector matrix.

Journal of Hydrology↗

Comparison of the lognormal and beta distribution functions to describe the uncertainty in permeability

The permeability of a single hydrostratigraphic unit is associated with considerable uncertainty due to measurement errors and significant spatial variability. Historically this uncertainty is characterized by a lognormal distribution. This distribution is generally heavy tailed, so using this distribution to describe the permeability has the limitation that all positive values of permeability for a given hydrostratigraphic unit have positive (not zero) probability of occurrence. The beta distribution, with its bounded domain, is explored as an alternative to the lognormal distribution in describing the uncertainty of permeability. The lognormal distribution and the beta distribution are both fit to the historic data provided by Jan Law in 1944 that was used to generalize the statement that the uncertainty in the permeability follows a lognormal distribution. The lognormal distribution and the beta distribution are also fit to an extensive permeability data set taken from regions within the Dakota Sandstone that have been shown to represent a single hydrostratigraphic unit. The results of this curve fitting exercise indicate that the beta distribution is a viable alternative to the lognormal distribution to characterize the uncertainty in permeability. ?? 2005 Elsevier B.V. All rights reserved.

Journal of Hydrology↗

Repeated surveys by acoustic Doppler current profiler for flow and sediment dynamics in a tidal river

A strategy of repeated surveys by acoustic Doppler current profiler (ADCP) was applied in a tidal river to map velocity vectors and suspended-sediment indicators. The Sacramento River at the junction with the Delta Cross Channel at Walnut Grove, California, was surveyed over several tidal cycles in the Fall of 2000 and 2001 with a vessel-mounted ADCP. Velocity profiles were recorded along flow-defining survey paths, with surveys repeated every 27 min through a diurnal tidal cycle. Velocity vectors along each survey path were interpolated to a three-dimensional Cartesian grid that conformed to local bathymetry. A separate array of vectors was interpolated onto a grid from each survey. By displaying interpolated vector grids sequentially with computer animation, flow dynamics of the reach could be studied in three-dimensions as flow responded to the tidal cycle. Velocity streamtraces in the grid showed the upwelling of flow from the bottom of the Sacramento River channel into the Delta Cross Channel. The sequential display of vector grids showed that water in the canal briefly returned into the Sacramento River after peak flood tides, which had not been known previously. In addition to velocity vectors, ADCP data were processed to derive channel bathymetry and a spatial indicator for suspended-sediment concentration. Individual beam distances to bed, recorded by the ADCP, were transformed to yield bathymetry accurate enough to resolve small bedforms within the study reach. While recording velocity, ADCPs also record the intensity of acoustic backscatter from particles suspended in the flow. Sequential surveys of backscatter intensity were interpolated to grids and animated to indicate the spatial movement of suspended sediment through the study reach. Calculation of backscatter flux through cross-sectional grids provided a first step for computation of suspended-sediment discharge, the second step being a calibrated relation between backscatter intensity and sediment concentration. Spatial analyses of ADCP data showed that a strategy of repeated surveys and flow-field interpolation has the potential to simplify computation of flow and sediment discharge through complex waterways. The use of trade, product, industry, or firm names in this report is for descriptive purposes only and does not constitute endorsement of products by the US Government. ?? 2005 Elsevier B.V. All rights reserved.

Journal of Hydrology↗

Estimating recharge using relations between precipitation and yield in a mountainous area with large variability in precipitation

Estimates of recharge to bedrock aquifers from infiltration of precipitation can be difficult to obtain, especially in areas with large spatial and temporal variability in precipitation. In the Black Hills area of western South Dakota and eastern Wyoming, streamflow yield is highly influenced by annual precipitation, with yield efficiency (annual yield divided by annual precipitation) increasing with increasing annual precipitation. Spatial variability in annual yield characteristics for Black Hills streams is predictably influenced by precipitation patterns. Relations between precipitation and yield efficiency were used to estimate annual recharge from long-term records of annual precipitation. A series of geographic information system algorithms was used to derive annual estimates for 1000- by 1000-m grid cells. These algorithms were composited to derive estimates of annual recharge rates to the Madison and Minnelusa aquifers in the Black Hills area of western South Dakota and eastern Wyoming during water years 1931-1998 and an estimate of average recharge for water years 1950-1998. This approach provides a systematic method of obtaining consistent and reproducible estimates of recharge from infiltration of precipitation. Resulting estimates of average annual recharge (water years 1950-1998) ranged from 1 cm in the southern Black Hills to 22 cm in the northwestern Black Hills. Recharge rates to these aquifers from infiltration of precipitation on outcrops was estimated to range from 0.9 m 3/s in 1936 to 18.8 m3/s in 1995.

Journal of Hydrology↗