USGS Science⌕ Search

SEARCH · USGS Science

Results for “Hydrology Journal”

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 559 records · Page 31Linked to original sources

Estimation of hectare-scale soil-moisture characteristics from aquifer-test data

Analysis of a 72-h, constant-rate aquifer test conducted in a coarse-grained and highly permeable, glacial outwash deposit on Cape Cod, Massachusetts revealed that drawdowns measured in 20 piezometers located at various depths below the water table and distances from the pumped well were significantly influenced by effects of drainage from the vadose zone . The influence was greatest in piezometers located close to the water table and diminished with increasing depth. The influence of the vadose zone was evident from a gap, in the intermediate-time zone, between measured drawdowns and drawdowns computed under the assumption that drainage from the vadose zone occurred instantaneously in response to a decline in the elevation of the water table. By means of an analytical model that was designed to account for time-varying drainage, simulated drawdowns could be closely fitted to measured drawdowns regardless of the piezometer locations. Because of the exceptional quality and quantity of the data and the relatively small aquifer heterogeneity, it was possible by inverse modeling to estimate all relevant aquifer parameters and a set of three empirical constants used in the upper-boundary condition to account for the dynamic drainage process. The empirical constants were used to define a one-dimensional (1D) drainage versus time curve that is assumed to be representative of the bulk material overlying the water table. The curve was inverted with a parameter estimation algorithm and a 1D numerical model for variably saturated flow to obtain soil-moisture retention curves and unsaturated hydraulic conductivity relationships defined by the Brooks and Corey equations. Direct analysis of the aquifer-test data using a parameter estimation algorithm and a two-dimensional, axisymmetric numerical model for variably saturated flow yielded similar soil-moisture characteristics. Results suggest that hectare-scale soil-moisture characteristics are different from core-scale predictions and even relatively small amounts of fine-grained material and heterogeneity can dominate the large-scale soil-moisture characteristics and aquifer response.

Journal of Hydrology↗

Chloride mass-balance method for estimating ground water recharge in arid areas: Examples from western Saudi Arabia

The chloride mass-balance method, which integrates time and aerial distribution of ground water recharge, was applied to small alluvial aquifers in the wadi systems of the Asir and Hijaz mountains in western Saudi Arabia. This application is an extension of the method shown to be suitable for estimating recharge in regional aquifers in semi-arid areas. Because the method integrates recharge in time and space it appears to be, with certain assumptions, particularly well suited for and areas with large temporal and spatial variation in recharge. In general, recharge was found to be between 3 to 4% of precipitation - a range consistent with recharge rates found in other arid and semi-arid areas of the earth.

Asir mountains, Hijaz mountains↗

Shallow ground-water quality beneath a major urban center: Denver, Colorado, USA

A survey of the chemical quality of ground water in the unconsolidated alluvial aquifer beneath a major urban center (Denver, Colorado, USA) was performed in 1993 with the objective of characterizing the quality of shallow ground-water in the urban area and relating water quality to land use. Thirty randomly selected alluvial wells were each sampled once for a broad range of dissolved constituents. The urban land use at each well site was sub- classified into one of three land-use settings: residential, commercial, and industrial. Shallow ground-water quality was highly variable in the urban area and the variability could be related to these land-use setting classifications. Sulfate (SO4) was the predominant anion in most samples from the residential and commercial land-use settings, whereas bicarbonate (HCO3) was the predominant anion in samples from the industrial land-use setting, indicating a possible shift in redox conditions associated with land use. Only three of 30 samples had nitrate concentrations that exceeded the US national drinking-water standard of 10 mg l-1 as nitrogen, indicating that nitrate contamination of shallow ground water may not be a serious problem in this urban area. However, the highest median nitrate concentration (4.2 mg l-1) was in samples from the residential setting, where fertilizer application is assumed to be most intense. Twenty-seven of 30 samples had detectable pesticides and nine of 82 analyzed pesticide compounds were detected at low concentrations, indicating that pesticides are widely distributed in shallow ground water in this urban area. Although the highest median total pesticide concentration (0.17 ??g l-1) was in the commercial setting, the herbicides prometon and atrazine were found in each land-use setting. Similarly, 25 of 29 samples analyzed had detectable volatile organic compounds (VOCs) indicating these compounds are also widely distributed in this urban area. The total VOC concentrations in sampled wells ranged from nondetectable to 23 442 ??g l-1. Widespread detections and occasionally high concentrations point to VOCs as the major anthropogenic ground-water impact in this urban environment. Generally, the highest VOC concentrations occurred in samples from the industrial setting. The most frequently detected VOC was the gasoline additive methyl tertbutyl ether (MTBE, in 23 of 29 wells). Results from this study indicate that the quality of shallow ground water in major urban areas can be related to land-use settings. Moreover, some VOCs and pesticides may be widely distributed at low concentrations in shallow ground water throughout major urban areas. As a result, the differentiation between point and non-point sources for these compounds in urban areas may be difficult.

Journal of Hydrology↗

Denitrification and mixing in a stream-aquifer system: Effects on nitrate loading to surface water

Ground water in terrace deposits of the South Platte River alluvial aquifer near Greeley, Colorado, USA, had a median nitrate concentration of 1857 μmol l −1 . Median nitrate concentrations in ground water from adjacent floodplain deposits (468 μmol l −1 ) and riverbed sediments (461 μmol l −1 ), both of which are downgradient from the terrace deposits, were lower than the median concentration in the terrace deposits. The concentrations and δ 15 N values of nitrate and N 2 in ground water indicated that denitrifying activity in the floodplain deposits and riverbed sediments accounted for 15–30% of the difference in nitrate concentrations. Concentrations of Cl − and SiO 2 indicated that mixing between river water and ground water in the floodplain deposits and riverbed sediments accounted for the remainder of the difference in nitrate concentrations. River flux measurements indicated that ground-water discharge in a 7.5 km segment of river had a nitrate load of 1718 kg N day − and accounted for about 18% of the total nitrate load in the river at the downstream end of that segment. This nitrate load was 70% less than the load predicted on the basis of the median nitrate concentration in the terrace deposits and assuming no denitrification or mixing in the aquifer. Water exchange between the river and aquifer caused ground water that originally discharged to the river to reenter denitrifying sediments in the riverbed and floodplain, thereby further decreasing the nitrate load in this stream—aquifer system. Results from this study indicated that denitrification and mixing within akluvial aquifer sediments may substantially decrease the nitrate load added to rivers by discharging ground water.

Journal of Hydrology↗

Dynamics of water-table fluctuations in an upland between two prairie-pothole wetlands in North Dakota

Data from a string of instrumented wells located on an upland of 55 m width between two wetlands in central North Dakota, USA, indicated frequent changes in water-table configuration following wet and dry periods during 5 years of investigation. A seasonal wetland is situated about 1.5 m higher than a nearby semipermanent wetland, suggesting an average ground water-table gradient of 0.02. However, water had the potential to flow as ground water from the upper to the lower wetland during only a few instances. A water-table trough adjacent to the lower semipermanent wetland was the most common water-table configuration during the first 4 years of the study, but it is likely that severe drought during those years contributed to the longevity and extent of the water-table trough. Water-table mounds that formed in response to rainfall events caused reversals of direction of flow that frequently modified the more dominant water-table trough during the severe drought. Rapid and large water-table rise to near land surface in response to intense rainfall was aided by the thick capillary fringe. One of the wettest summers on record ended the severe drought during the last year of the study, and caused a larger-scale water-table mound to form between the two wetlands. The mound was short in duration because it was overwhelmed by rising stage of the higher seasonal wetland which spilled into the lower wetland. Evapotranspiration was responsible for generating the water-table trough that formed between the two wetlands. Estimation of evapotranspiration based on diurnal fluctuations in wells yielded rates that averaged 3–5 mm day −1 . On many occasions water levels in wells closer to the semipermanent wetland indicated a direction of flow that was different from the direction indicated by water levels in wells farther from the wetland. Misinterpretation of direction and magnitude of gradients between ground water and wetlands could result from poorly placed or too few observation wells, and also from infrequent measurement of water levels in wells.

Journal of Hydrology↗

Effect of nitrate, organic carbon, and temperature on potential denitrification rates in nitrate-rich riverbed sediments

A study conducted in 1994 as part of the US Geological Survey's National Water-Quality Assessment Program, South Platte River Basin investigation, examined the effect of certain environmental factors on potential denitrification rates in nitrate-rich riverbed sediments. The acetylene block technique was used to measure nitrous oxide (N2O) production rates in laboratory incubations of riverbed sediments to evaluate the effect of varying nitrate concentrations, organic carbon concentrations and type, and water temperature on potential denitrification rates. Sediment incubations amended with nitrate, at concentrations ranging from 357 to 2142 ??mol l-1 (as measured in the field), produced no significant increase (P > 0.05) in N2O production rates, indicating that the denitrification potential in these sediments was not nitrate limited. In contrast, incubations amended with acetate as a source of organic carbon, at concentrations ranging from 0 to 624 ??mol l-1, produced significant increases (P < 0.05) in N2O production rates with increased organic carbon concentration, indicating that the denitrification potential in these sediments was organic carbon limited. Furthermore, N2O production rates also were affected by the type of organic carbon available as an electron donor. Acetate and surface-water-derived fulvic acid supported higher N2O production rates than groundwater-derived fulvic acid or sedimentary organic carbon. Lowering incubation temperatures from 22 to 4??C resulted in about a 77% decrease in the N2O production rates. These results help to explain findings from previous studies indicating that only 15-30% of nitrate in groundwater was denitrified before discharging to the South Platte River and that nitrate concentrations in the river generally were higher in winter than in summer.

Journal of Hydrology↗

Evaluation of unconfined-aquifer parameters from pumping test data by nonlinear least squares

Nonlinear least squares (NLS) with automatic differentiation was used to estimate aquifer parameters from drawdown data obtained from published pumping tests conducted in homogeneous, water-table aquifers. The method is based on a technique that seeks to minimize the squares of residuals between observed and calculated drawdown subject to bounds that are placed on the parameter of interest. The analytical model developed by Neuman for flow to a partially penetrating well of infinitesimal diameter situated in an infinite, homogeneous and anisotropic aquifer was used to obtain calculated drawdown. NLS was first applied to synthetic drawdown data from a hypothetical but realistic aquifer to demonstrate that the relevant hydraulic parameters (storativity, specific yield, and horizontal and vertical hydraulic conductivity) can be evaluated accurately. Next the method was used to estimate the parameters at three field sites with widely varying hydraulic properties. NLS produced unbiased estimates of the aquifer parameters that are close to the estimates obtained with the same data using a visual curve-matching approach. Small differences in the estimates are a consequence of subjective interpretation introduced in the visual approach.

Journal of Hydrology↗

Results of a prototype surface water network design for pesticides developed for the San Joaquin River Basin, California

A nested surface water monitoring network was designed and tested to measure variability in pesticide concentrations in the San Joaquin River and selected tributaries during the irrigation season. The network design an d sampling frequency necessary for determining the variability and distribution in pesticide concentrations were tested in a prototype study. The San Joaquin River Basin, California, was sampled from April to August 1992, a period during the irrigation season where there was no rainfall. Orestimba Creek, which drains a part of the western San Joaquin Valley, was sampled three times per week for 6 weeks, followed by a once per week sampling for 6 weeks, and the three times per week sampling for 6 weeks. A site on the San Joaquin River near the mouth of the basin, and an irrigation drain of the eastern San Joaquin Valley, were sampled weekly during the entire sampling period. Pesticides were most often detected in samples collected from Orestimba Creek. This suggests that the western valley was the principal source of pesticides to the San Joaquin River during the irrigation season. Irrigation drainage water was the source of pesticides to Orestimba Creek. Pesticide concentrations of Orestimba Creek showed greater temporal variability when sampled three times per week than when sampled once a week, due to variations in field management and irrigation. The implication for the San Joaquin River basin (an irrigation-dominated agricultural setting) is that frequent sampling of tributary sites is necessary to describe the variability in pesticides transported to the San Joaquin River.

Journal of Hydrology↗

222Rn transport in a fractured crystalline rock aquifer: Results from numerical simulations

Dissolved 222Rn concentrations in ground water from a small wellfield underlain by fractured Middle Proterozoic Pikes Peak Granite southwest of Denver, Colorado range from 124 to 840 kBq m-3 (3360-22700 pCi L-1). Numerical simulations of flow and transport between two wells show that differences in equivalent hydraulic aperture of transmissive fractures, assuming a simplified two-fracture system and the parallel-plate model, can account for the different 222Rn concentrations in each well under steady-state conditions. Transient flow and transport simulations show that 222Rn concentrations along the fracture profile are influenced by 222Rn concentrations in the adjoining fracture and depend on boundary conditions, proximity of the pumping well to the fracture intersection, transmissivity of the conductive fractures, and pumping rate. Non-homogeneous distribution (point sources) of 222Rn parent radionuclides, uranium and 226Ra, can strongly perturb the dissolved 222Rn concentrations in a fracture system. Without detailed information on the geometry and hydraulic properties of the connected fracture system, it may be impossible to distinguish the influence of factors controlling 222Rn distribution or to determine location of 222Rn point sources in the field in areas where ground water exhibits moderate 222Rn concentrations. Flow and transport simulations of a hypothetical multifracture system consisting of ten connected fractures, each 10 m in length with fracture apertures ranging from 0.1 to 1.0 mm, show that 222Rn concentrations at the pumping well can vary significantly over time. Assuming parallel-plate flow, transmissivities of the hypothetical system vary over four orders of magnitude because transmissivity varies with the cube of fracture aperture. The extreme hydraulic heterogeneity of the simple hypothetical system leads to widely ranging 222Rn values, even assuming homogeneous distribution of uranium and 226Ra along fracture walls. Consequently, it is concluded that 222Rn concentrations vary, not only with the geometric and stress factors noted above, but also according to local fracture aperture distribution, local groundwater residence time, and flux of 222Rn from parent radionuclides along fracture walls.

Journal of Hydrology↗

Statistical modeling of agricultural chemical occurrence in midwestern rivers

Agricultural chemicals in surface water may constitute a human health risk or have adverse effects on aquatic life. Recent research on unregulated rivers in the midwestern USA documents that elevated concentrations of herbicides occur for 1-4 months following application in late spring and early summer. In contrast, nitrate concentrations in unregulated rivers are elevated during fall, winter, and spring months. Natural and anthropogenic variables of fiver drainage basins, such as soil permeability, amount of agricultural chemicals applied, or percentage of land planted in corn, affect agricultural chemical concentration and mass transport in rivers. Presented is an analysis of selected data on agricultural chemicals collected for three regional studies conducted by the US Geological Survey. Statistical techniques such as multiple linear and logistic regression were used to identify natural and anthropogenic variables of drainage basins that have strong relations to agricultural chemical concentrations and mass transport measured in rivers. A geographic information system (GIS) was used to manage and analyze spatial data. Statistical models were developed that estimated the concentration, annual transport, and annual mean concentration of selected agricultural chemicals in midwestern rivers. Multiple linear regression models were not very successful (R2 from 0.162 to 0.517) in explaining the variance in observed agricultural chemical concentrations during post-planting runoff. Logistic regression models were somewhat more successful, correctly matching the observed concentration category in 61-80% of observations. Linear and multiple linear regression models were moderately successful (R2 from 0.522 to 0.995) in explaining the variance in observed annual transport and annual mean concentration of agricultural chemicals. Explanatory variables that were commonly significant in the regression models include estimates of agricultural chemical use, crop acreage, soil characteristics, and basin topography.

Journal of Hydrology↗

Modeling saltwater upconing in a freshwater aquifer in south-central Kansas

The Great Bend Prairie freshwater alluvial aquifer in south-central Kansas overlies a bedrock brine aquifer of Permian age. The continuous extraction of freshwater mainly for irrigation in this area has accelerated the upward movement of the saltwater, resulting in the deterioration of water quality. Predicting saltwater upconing is critical for maintaining a long-term supply of water of good quality to the Groundwater Management District No. 5. This paper uses a numerical model to predict the effect of saltwater upconing on the salinity of pumped water. A sensitivity analysis was conducted for estimating the effects of uncertainties in model parameters on the numerical results. The most important factors affecting the salinity of discharged water are found to be the location and nature of clay layers in the aquifer, the pumping rate, the location of the well screen, and the hydraulic conductivities of the medium. The effect of uncertainties in aquifer porosity and dispersivity on salinity is appreciable but not substantial, whereas that of recharge from precipitation in the study area is relatively insignificant. In addition, a numerical model based on the field conditions observed at the Siefkes site was constructed and calibrated to reproduce and project the variation of measured water levels and discharged groundwater concentration. Finally, a number of practical management recommendations based on this study are presented.

Journal of Hydrology↗

Multi-level slug tests in highly permeable formations: 2. Hydraulic conductivity identification, method verification, and field applications

Using the developed theory and modified Springer-Gelhar (SG) model, an identification method is proposed for estimating hydraulic conductivity from multi-level slug tests. The computerized algorithm calculates hydraulic conductivity from both monotonic and oscillatory well responses obtained using a double-packer system. Field verification of the method was performed at a specially designed fully penetrating well of 0.1-m diameter with a 10-m screen in a sand and gravel alluvial aquifer (MSEA site, Shelton, Nebraska). During well installation, disturbed core samples were collected every 0.6 m using a split-spoon sampler. Vertical profiles of hydraulic conductivity were produced on the basis of grain-size analysis of the disturbed core samples. These results closely correlate with the vertical profile of horizontal hydraulic conductivity obtained by interpreting multi-level slug test responses using the modified SG model. The identification method was applied to interpret the response from 474 slug tests in 156 locations at the MSEA site. More than 60% of responses were oscillatory. The method produced a good match to experimental data for both oscillatory and monotonic responses using an automated curve matching procedure. The proposed method allowed us to drastically increase the efficiency of each well used for aquifer characterization and to process massive arrays of field data. Recommendations generalizing this experience to massive application of the proposed method are developed.Using the developed theory and modified Springer-Gelhar (SG) model, an identification method is proposed for estimating hydraulic conductivity from multi-level slug tests. The computerized algorithm calculates hydraulic conductivity from both monotonic and oscillatory well responses obtained using a double-packer system. Field verification of the method was performed at a specially designed fully penetrating well of 0.1-m diameter with a 10-m screen in a sand and gravel alluvial aquifer (MSEA site, Shelton, Nebraska). During well installation, disturbed core samples were collected every 0.6 m using a split-spoon sampler. Vertical profiles of hydraulic conductivity were produced on the basis of grain-size analysis of the disturbed core samples. These results closely correlate with the vertical profile of horizontal hydraulic conductivity obtained by interpreting multi-level slug test responses using the modified SG model. The identification method was applied to interpret the response from 474 slug tests in 156 locations at the MSEA site. More than 60% of responses were oscillatory. The method produced a good match to experimental data for both oscillatory and monotonic responses using an automated curve matching procedure. The proposed method allowed us to drastically increase the efficiency of each well used for aquifer characterization and to process massive arrays of field data. Recommendations generalizing this experience to massive application of the proposed method are developed.

Journal of Hydrology↗

Multi-level slug tests in highly permeable formations: 1. Modification of the Springer-Gelhar (SG) model

A multi-level slug test model and a method for the evaluation of vertical profiles of hydraulic conductivity in highly permeable formations are developed. A double-packer system is employed to estimate local hydraulic conductivity. Depending on the formation and double-packer system parameters, the water level recovery in the tested well can exhibit a monotonic or oscillatory response. To discern information on aquifer properties from artifacts introduced by the measurement system, the theory of flow in an aquifer and a double-packer system is developed. The mathematical model incorporates features of the Springer and Gelhar (1991) model and reduces to the Bouwer and Rice (1976) model in a special case. The model involves equations of momentum and mass conservation for the double-packer system with quasi-steady well-aquifer interaction equations. The method is uniformly applicable for both monotonic and oscillatory well responses and can produce profiles of hydraulic conductivity for the tested well. The criterion is given to determine the type of well response for given slug test conditions.A multi-level slug test model and a method for the evaluation of vertical profiles of hydraulic conductivity in highly permeable formations are developed. A double-packer system is employed to estimate local hydraulic conductivity. Depending on the formation and double-packer system parameters, the water level recovery in the tested well can exhibit a monotonic or oscillatory response. To discern information on aquifer properties from artifacts introduced by the measurement system, the theory of flow in an aquifer and a double-packer system is developed. The mathematical model incorporates features of the Springer and Gelhar (1991) model and reduces to the Bouwer and Rice (1976) model in a special case. The model involves equations of momentum and mass conservation for the double-packer system with quasi-steady well-aquifer interaction equations. The method is uniformly applicable for both monotonic and oscillatory well responses and can produce profiles of hydraulic conductivity for the tested well. The criterion is given to determine the type of well response for given slug test conditions.

Journal of Hydrology↗

Effects of a beaver pond on runoff processes: comparison of two headwater catchments

Natural variations in concentrations of 18 O, D, and H 4 SiO 4 in two tributary catchments of Woods Lake in the west-central Adirondack Mountains of New York were measured during 1989&ndash;1991 to examine runoff processes and their implications for the neutralization of acidic precipitation by calcium carbonate treatment. The two catchments are similar except that one contained a 1.3 ha beaver pond. Evaporation from the beaver pond caused a seasonal decrease in the slope of the meteoric water line in stream water from the catchment with a beaver pond (WO2). No corresponding change in slope of the meteoric water line was evident in stream water from the other catchment (WO4), nor in ground water nor soil water from either catchment, indicating that evaporative fractionation was not significant. Application of a best-fit sine curve to &delta; 18 O data indicated that base flow in both catchments had a residence time of about 100 days. Ground water from a well finished in thick till had the longest residence time (160 days); soil water from the O-horizon and B-horizon had residence times of 63 and 80 days, respectively. Water previously stored within each catchment (pre-event water) was the predominant component of streamflow during spring snowmelt and during spring and autumn rainfall events, but the proportion of streamflow that consisted of pre-event water differed significantly in the two catchments. The proportion of event water (rain and snowmelt) in WO2 was smaller than at WO4 early in the spring snowmelt of March 13&ndash;17, 1990, but the proportions of source water components for the two catchments were almost indistinguishable by the peak flow on the third day of the melt. The event water was further separated into surface-water and subsurface-water components by utilizing measured changes in H 4 SiO 4 concentrations in stream water during the snowmelt. Results indicated that subsurface flow was the dominant pathway by which event water reached the stream except during the peak flow of a rain-on-snow event on the last day of the melt. Streamflow from a spring rain storm with dry antecendent conditions two months later (May 16&ndash;18, 1990), was less than 5% event water at peak flow in WO2 and 26% in WO4. This change from the runoff pattern in March is attributed to retention of event water in the beaver pond favored by relatively low pre-event storage and isothermal (nonstratified) conditions in the pond that allowed mixing. Streamflow during several autumn storms was about 15&ndash;25% event water at peak flow in WO4; the highest values for event water were associated with wet antecedent moisture conditions. These results indicate that a beaver pond can significantly affect the downstream delivery of event water through evaporation and mixing, but provides minimal retention during large runoff events such as snowmelt. Beaver ponds are expected to provide greater opportunity for neutralization of acidic waters during most of the year in catchments treated with calcium carbonate, but little neutralization effect during snowmelt.

Journal of Hydrology↗

Infiltration and hydraulic connections from the Niagara River to a fractured-dolomite aquifer in Niagara Falls, New York

The spatial distribution of hydrogen and oxygen stable-isotope values in groundwater can be used to distinguish different sources of recharge and to trace groundwater flow directions from recharge boundaries. This method can be particularly useful in fractured-rock settings where multiple lines of evidence are required to delineate preferential flow paths that result from heterogeneity within fracture zones. Flow paths delineated with stable isotopes can be combined with hydraulic data to form a more complete picture of the groundwater flow system. In this study values of ??D and ??18O were used to delineate paths of river-water infiltration into the Lockport Group, a fractured dolomite aquifer, and to compute the percentage of fiver water in groundwater samples from shallow bedrock wells. Flow paths were correlated with areas of high hydraulic diffusivity in the shallow bedrock that were delineated from water-level fluctuations induced by diurnal stage fluctuations in man-made hydraulic structures. Flow paths delineated with the stable-isotope and hydraulic data suggest that fiver infiltration reaches an unlined storm sewer in the bedrock through a drainage system that surrounds aqueducts carrying river water to hydroelectric power plants. This finding is significant because the storm sewer is the discharge point for contaminated groundwater from several chemical waste-disposal sites and the cost of treating the storm sewer's discharge could be reduced if the volume of infiltration from the river were decreased.The spatial distribution of hydrogen and oxygen stable-isotope values in groundwater can be used to distinguish different sources of recharge and to trace groundwater flow directions from recharge boundaries. This method can be particularly useful in fractured-rock settings where multiple lines of evidence are required to delineate preferential flow paths that result from heterogeneity within fracture zones. Flow paths delineated with stable isotopes can be combined with hydraulic data to form a more complete picture of the groundwater flow system. In this study values of ??D and ??18O were used to delineate paths of river-water infiltration into the Lockport Group, a fractured dolomite aquifer, and to compute the percentage of river water in groundwater samples from shallow bedrock wells. Flow paths were correlated with areas of high hydraulic diffusivity in the shallow bedrock that were delineated from water-level fluctuations induced by diurnal stage fluctuations in man-made hydraulic structures. Flow paths delineated with the stable-isotope and hydraulic data suggest that river infiltration reaches an unlined storm sewer in the bedrock through a drainage system that surrounds aqueducts carrying river water to hydroelectric power plants. This finding is significant because the storm sewer is the discharge point for contaminated groundwater from several chemical waste-disposal sites and the cost of treating the storm sewer's discharge could be reduced if the volume of infiltration from the river were decreased.

Journal of Hydrology↗

Analysis and simulation of reactive transport of metal contaminants in ground water in Pinal Creek Basin, Arizona

Large-scale mining activities have generated a plume of acidic ground water more than 15 km long in the regional aquifer of the Pinal Creek Basin. A one-dimensional reactive-transport model was developed using PHREEQC to aid in the analysis of transport and chemical processes in the plume and to determine the uses and limitations of this type of modeling approach. In 1984, the acidic part of the plume had a pH as low as 3.4 and contained milligram-per-liter concentrations of iron, copper, aluminum and other metals. From 1984 to 1994, concentrations of contaminants in the alluvial aquifer in Pinal Creek Basin, Arizona, decreased as a result of mixing, recharge, remedial pumping and chemical reactions. For reactions involving gypsum and rhodochrosite, the equilibrium modeling assumption of a local geochemical equilibrium was generally valid. From 1984 to 1990, water along the simulated flow path was at equilibrium or slightly supersaturated with gypsum, and gypsum equilibria controlled dissolved concentrations of calcium and sulfate. Beginning in 1991, water in the acidic part of the plume became increasingly undersaturated with respect to gypsum, indicating that the gypsum available for dissolution in the aquifer may have been completely consumed by about 1991. Rhodochrosite precipitation was thought responsible for the measured attenuation in dissolved manganese in the neutralized zone. For reactions involving calcite, the assumption of a local geochemical equilibrium was generally not valid. Dissolution of calcite in the transition zone was not sufficient to establish equilibrium although, following neutralization, the calcite saturation index decreased to −1.2 in 1986. Calcite undersaturation decreased along the flow path in the neutralized zone, and equilibrium was attained about 7 km downgradient of the transition zone. The assumption of a local geochemical equilibrium was not valid for oxidation–reduction reactions that involved iron oxides and manganese oxides. Kinetically controlled oxidation–reduction reactions continued in the acidic part of the flow path for years following the passage of the transition zone. Although the equilibrium approach helped to provide an increased understanding of contaminant transport at Pinal Creek, future work will require a kinetic modeling approach to more accurately simulate selected reactions between the plume and aquifer materials.

Arizona↗

Estimating formation properties from early-time recovery in wells subject to turbulent head losses

A mathematical model is developed to interpret the early-time recovering water level following the termination of pumping in wells subject to turbulent head losses. The model assumes that turbulent head losses dissipate immediately when pumping ends. In wells subject to both borehole storage and turbulent head losses, the early-time recovery exhibits a slope equal to 1/2 on log-log plots of the recovery versus time. This half-slope response should not be confused with the half-slope response associated with a linear flow regime during aquifer tests. The presence of a borehole skin due to formation damage or stimulation around the pumped well alters the early-time recovery in wells subject to turbulent head losses and gives the appearance of borehole storage, where the recovery exhibits a unit slope on log-log plots of recovery versus time. Type curves can be used to estimate the formation storafivity from the early-time recovery data. In wells that are suspected of having formation damage or stimulation, the type curves can be used to estimate the 'effective' radius of the pumped well, if an estimate of the formation storativity is available from observation wells or other information. Type curves for a homogeneous and isotropic dual-porosity aquifer are developed and applied to estimate formation properties and the effect of formation stimulation from a single-well test conducted in the Madison limestone near Rapid City, South Dakota.A mathematical model is developed to interpret the early-time recovering water level following the termination of pumping in wells subject to turbulent head losses. The model assumes that turbulent head losses dissipate immediately when pumping ends. In wells subject to both borehole storage and turbulent head losses, the early-time recovery exhibits a slope equal to 1/2 on log-log plots of the recovery versus time. This half-slope response should not be confused with the half-slope response associated with a linear flow regime during aquifer tests. The presence of a borehole skin due to formation damage or stimulation around the pumped well alters the early-time recovery in wells subject to turbulent head losses and gives the appearance of borehole storage, where the recovery exhibits a unit slope on log-log plots of recovery versus time. Type curves can be used to estimate the formation storativity from the early-time recovery data. In wells that are suspected of having formation damage or stimulation, the type curves can be used to estimate the `effective' radius of the pumped well, if an estimate of the formation storativity is available from observation wells or other information. Type curves for a homogeneous and isotropic dual-porosity aquifer are developed and applied to estimate formation properties and the effect of formation stimulation from a single-well test conducted in the Madison limestone near Rapid City, South Dakota.

Journal of Hydrology↗

Noble gases, stable isotopes, and radiocarbon as tracers of flow in the Dakota aquifer, Colorado and Kansas

A suite of chemical and isotope tracers (dissolved noble gases, stable isotopes of water, radiocarbon, and CI) have been analyzed along a flow path in the Dakota aquifer system to determine likely recharge sources, ground water residence times, and the extent of mixing between local and intermediate flow systems, presumably caused by large well screens. Three water types were distinguished with the tracers, each having a very different history. Two of the water types were found in south-eastern Colorado where the Dakota is poorly confined. The tracer data suggest that the first group recharged locally during the last few thousand years and the second group was composed of ground water that recharged earlier during a cooler climate, presumably during the last glacial period (LGP) and mixed aged water. The paleotemperature record archived in this groundwater system indicates that south-eastern Colorado was about 5??C cooler during the LGP than during the late Holocene. Similar temperature changes derived from dissolved noble gases in other aquifer systems have been reported earlier for the south-western United States. The third water type was located down gradient of the first two in the confined Dakota in western and central Kansas. Groundwater residence time of this water mass is on the order of 104-105 yrs and its recharge location is near the Colorado and Kansas border down gradient of the other water types. The study shows the importance of using multiple tracers when investigating ground water systems.A suite of chemical and isotope tracers (dissolved noble gases, stable isotopes of water, radiocarbon, and CL) were analyzed along a flow path in the Dakota aquifer system to determine likely recharge sources, ground water residence times, and the extent of mixing between local and intermediate flow systems. Three water types were distinguished with the tracers, each having a very different history. Two of the water types were located in south-eastern Colorado where the Dakota is poorly confined. The third water type was located down gradient of the first two in the confined Dakota in western and central Kansas.

Journal of Hydrology↗