USGS Science⌕ Search

SEARCH · USGS Science

Results for “Water Resources Research”

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
USGS2002

Multispecies reactive tracer test in an aquifer with spatially variable chemical conditions, Cape Cod, Massachusetts: Dispersive transport of bromide and nickel

Dispersive transport of groundwater solutes was investigated as part of a multispecies reactive tracer test conducted under spatially variable chemical conditions in an unconfined, sewage‐contaminated sand and gravel aquifer on Cape Cod, Massachusetts. Transport of the nonreactive tracer bromide (Br) reflected physical and hydrologic processes. Transport of the reactive tracer nickel (Ni) complexed with an organic ligand (NiEDTA) varied in response to pH and other chemical conditions within the aquifer. A loss of about 14% of the Ni mass was calculated from the distribution of tracers through time. This loss is consistent with reversible adsorption of NiEDTA onto the iron and aluminum oxyhydroxide coatings on the aquifer sediments. The Ni consistently lagged behind Br with a calculated retardation coefficient of 1.2. Longitudinal dispersivities reached constant values of 2.2 and 1.1 m for Br and Ni, respectively, by at least 69 m of travel. The smaller dispersivity for Ni possibly was due to nonlinear or spatially variant adsorption of NiEDTA. In the upper, uncontaminated zone of the aquifer, longitudinal dispersion of Ni was greater than that of Br early in the test as a result of reversible adsorption of NiEDTA. In general, transverse dispersivities were much smaller (horizontal: 1.4–1.5 × 10 −2 m; vertical: 0.5–3.8 × 10 −3 m) than the longitudinal dispersivities. The Br results are similar to those from a test conducted eight years earlier, suggesting that transport parameters are spatially stationary within the aquifer at the scale of 300 m covered by the spatially overlapping tests. A significant difference between the two tests was the travel distance (69 and 26 m) needed to reach a constant longitudinal dispersivity.

Massachusetts↗
USGS2002

Consistency of patterns in concentration‐discharge plots

Concentration‐discharge (c‐Q) plots have been used to infer how flow components such as event water, soil water, and groundwater mix to produce the observed episodic hydrochemical response of small catchments. Because c‐Q plots are based only on observed streamflow and solute concentration, their interpretation requires assumptions about the relative volume, hydrograph timing, and solute concentration of the streamflow end‐members. Evans and Davies [1998] present a taxonomy of c‐Q loops resulting from three‐component conservative mixing. Their analysis, based on a fixed template of end‐member hydrograph volume, timing, and concentration, suggests a unique relationship between c‐Q loop form and the rank order of end‐member concentrations. Many catchments exhibit variability in component contributions to storm flow in response to antecedent conditions or rainfall characteristics, but the effects of such variation on c‐Q relationships have not been studied systematically. Starting with a “baseline” condition similar to that assumed by Evans and Davies [1998] , we use a simple computer model to characterize the variability in c‐Q plot patterns resulting from variation in end‐member volume, timing, and solute concentration. Variability in these three factors can result in more than one c‐Q loop shape for a given rank order of end‐member solute concentrations. The number of resulting hysteresis patterns and their relative frequency depends on the rank order of solute concentrations and on their separation in absolute value. In ambiguous cases the c‐Q loop shape is determined by the relative “prominence” of the event water versus soil water components. This “prominence” is broadly defined as a capacity to influence the total streamflow concentration and may result from a combination of end‐member volume, timing, or concentration. The modeling results indicate that plausible hydrological variability in field situations can confound the interpretation of c‐Q plots, even when fundamental end‐member mixing assumptions are satisfied.

Water Resources Research↗
USGS2003

Physical and hydraulic properties of volcanic rocks from Yucca Mountain, Nevada

A database of physical and hydraulic properties was developed for rocks in the unsaturated zone at Yucca Mountain, Nevada, a site under consideration as a geologic repository for high‐level radioactive waste. The 5320 core samples were collected from 23 shallow (<100 m) and 10 deep (500–1000 m) vertical boreholes. Hydrogeologic units have been characterized in the unsaturated zone [ Flint , 1998] that represent rocks with ranges of welding, lithophysae, and high and low temperature alteration (as a result of the depositional, cooling, and alterational history of the lithostratigraphic layers). Lithostratigraphy, the hydrogeologic unit, and the corresponding properties are described. In addition, the physical properties of bulk density, porosity, and particle density; the hydraulic properties of saturated hydraulic conductivity and moisture retention characteristics; and the field water content were measured and compiled for each core sample.

Nevada↗
USGS2003

Regional flood probabilities

The T ‐year annual maximum flood at a site is defined to be that streamflow, that has probability 1/ T of being exceeded in any given year, and for a group of sites the corresponding regional flood probability (RFP) is the probability that at least one site will experience a T ‐year flood in any given year. The RFP depends on the number of sites of interest and on the spatial correlation of flows among the sites. We present a Monte Carlo method for obtaining the RFP and demonstrate that spatial correlation estimates used in this method may be obtained with rank transformed data and therefore that knowledge of the at‐site peak flow distribution is not necessary. We examine the extent to which the estimates depend on specification of a parametric form for the spatial correlation function, which is known to be nonstationary for peak flows. It is shown in a simulation study that use of a stationary correlation function to compute RFPs yields satisfactory estimates for certain nonstationary processes. Application of asymptotic extreme value theory is examined, and a methodology for separating channel network and rainfall effects on RFPs is suggested. A case study is presented using peak flow data from the state of Washington. For 193 sites in the Puget Sound region it is estimated that a 100‐year flood will occur on the average every 4.5 years.

Water Resources Research↗
USGS2002

Steady shape analysis of tomographic pumping tests for characterization of aquifer heterogeneities

Hydraulic tomography, a procedure involving the performance of a suite of pumping tests in a tomographic format, provides information about variations in hydraulic conductivity at a level of detail not obtainable with traditional well tests. However, analysis of transient data from such a suite of pumping tests represents a substantial computational burden. Although steady state responses can be analyzed to reduce this computational burden significantly, the time required to reach steady state will often be too long for practical applications of the tomography concept. In addition, uncertainty regarding the mechanisms driving the system to steady state can propagate to adversely impact the resulting hydraulic conductivity estimates. These disadvantages of a steady state analysis can be overcome by exploiting the simplifications possible under the steady shape flow regime. At steady shape conditions, drawdown varies with time but the hydraulic gradient does not. Thus transient data can be analyzed with the computational efficiency of a steady state model. In this study, we demonstrate the value of the steady shape concept for inversion of hydraulic tomography data and investigate its robustness with respect to improperly specified boundary conditions.

Water Resources Research↗
USGS2003

Interpreting tracer breakthrough tailing from different forced-gradient tracer experiment configurations in fractured bedrock

Conceptual and mathematical models are presented that explain tracer breakthrough tailing in the absence of significant matrix diffusion. Model predictions are compared to field results from radially convergent, weak-dipole, and push-pull tracer experiments conducted in a saturated crystalline bedrock. The models are based upon the assumption that flow is highly channelized, that the mass of tracer in a channel is proportional to the cube of the mean channel aperture, and the mean transport time in the channel is related to the square of the mean channel aperture. These models predict the consistent −2 straight line power law slope observed in breakthrough from radially convergent and weak-dipole tracer experiments and the variable straight line power law slope observed in push-pull tracer experiments with varying injection volumes. The power law breakthrough slope is predicted in the absence of matrix diffusion. A comparison of tracer experiments in which the flow field was reversed to those in which it was not indicates that the apparent dispersion in the breakthrough curve is partially reversible. We hypothesize that the observed breakthrough tailing is due to a combination of local hydrodynamic dispersion, which always increases in the direction of fluid velocity, and heterogeneous advection, which is partially reversed when the flow field is reversed. In spite of our attempt to account for heterogeneous advection using a multipath approach, a much smaller estimate of hydrodynamic dispersivity was obtained from push-pull experiments than from radially convergent or weak dipole experiments. These results suggest that although we can explain breakthrough tailing as an advective phenomenon, we cannot ignore the relationship between hydrodynamic dispersion and flow field geometry at this site. The design of the tracer experiment can severely impact the estimation of hydrodynamic dispersion and matrix diffusion in highly heterogeneous geologic media.

Water Resources Research↗
USGS2002

Analysis of temperature profiles for investigating stream losses beneath ephemeral channels

Continuous estimates of streamflow are challenging in ephemeral channels. The extremely transient nature of ephemeral streamflows results in shifting channel geometry and degradation in the calibration of streamflow stations. Earlier work suggests that analysis of streambed temperature profiles is a promising technique for estimating streamflow patterns in ephemeral channels. The present work provides a detailed examination of the basis for using heat as a tracer of stream/groundwater exchanges, followed by a description of an appropriate heat and water transport simulation code for ephemeral channels, as well as discussion of several types of temperature analysis techniques to determine streambed percolation rates. Temperature‐based percolation rates for three ephemeral stream sites are compared with available surface water estimates of channel loss for these sites. These results are combined with published results to develop conclusions regarding the accuracy of using vertical temperature profiles in estimating channel losses. Comparisons of temperature‐based streambed percolation rates with surface water‐based channel losses indicate that percolation rates represented 30% to 50% of the total channel loss. The difference is reasonable since channel losses include both vertical and nonvertical component of channel loss as well as potential evapotranspiration losses. The most significant advantage of the use of sediment‐temperature profiles is their robust and continuous nature, leading to a long‐term record of the timing and duration of channel losses and continuous estimates of streambed percolation. The primary disadvantage is that temperature profiles represent the continuous percolation rate at a single point in an ephemeral channel rather than an average seepage loss from the entire channel.

Water Resources Research↗
USGS2003

Geochemical response to variable streamflow conditions in contaminated and uncontaminated streams

Seasonal variations in stream inorganic geochemistry are not well documented or understood. We sampled two mining-impacted and two relatively pristine streams in western Montana over a 12-month period, collecting samples every 4 weeks, with supplemental sampling (at least weekly) during spring runoff. We analyzed all samples for dissolved (operationally defined as <0.2 ??m) and total recoverable concentrations. Generally, the trace elements (Al, As, Cu, Fe, Mn, and Zn) did not correlate linearly with streamflow, while the major elements (e.g., Ca, K, and Mg) did. Suspended sediment, total recoverable metals, and H+ followed clockwise hysteresis rotations, driven by short-term flushing events during the very early stages of spring runoff. Mining-impacted sites had higher concentrations of many trace elements than did relatively pristine sites. One of the mining impacted sites exhibited strong geochemical responses to spring rain events in the basin. The results underscore the need to sample streams frequently during changing hydrologic and climatic conditions in order to accurately monitor surface water quality and to determine solute and particulate loads (both contaminant and noncontaminant).

Water Resources Research↗
USGS2003

Inverse modeling of interbed storage parameters using land subsidence observations, Antelope Valley, California

We use land-subsidence observations from repeatedly surveyed benchmarks and interferometric synthetic aperture radar (InSAR) in Antelope Valley, California, to estimate spatially varying compaction time constants, ??, and inelastic specific skeletal storage coefficients, Skv*, in a previously calibrated regional groundwater flow and subsidence model. The observed subsidence patterns reflect both the spatial distribution of head declines and the spatially variable inelastic skeletal storage coefficient. Using the nonlinear parameter estimation program UCODE we estimate compaction time constants between 3.8 and 285 years. The Skv* values are estimated by linear estimation and range from 0 to almost 0.09. We find that subsidence observations over long time periods are necessary to constrain estimates of the large compaction time constants in Antelope Valley. The InSAR data used in this study cover only a three-year period, limiting their usefulness in constraining these time constants. This problem will be alleviated as more SAR data become available in the future or where time constants are small. By incorporating the resulting parameter estimates in the previously calibrated regional model of groundwater flow and land subsidence we can significantly improve the agreement between simulated and observed land subsidence both in terms of magnitude and spatial extent. The sum of weighted squared subsidence residuals, a common measure of model fit, was reduced by 73% with respect to the original model. However, the ability of the model to adequately reproduce the subsidence observed over only a few years is impaired by the fact that the simulated hydraulic heads over small time periods are often not representative of the actual aquifer hydraulic heads. Errors in the simulated hydraulic aquifer heads constitute the primary limitation of the approach presented here.

Water Resources Research↗
USGS2003

Methods for using groundwater model predictions to guide hydrogeologic data collection, with application to the Death Valley regional groundwater flow system

Calibrated models of groundwater systems can provide substantial information for guiding data collection. This work considers using such models to guide hydrogeologic data collection for improving model predictions by identifying model parameters that are most important to the predictions. Identification of these important parameters can help guide collection of field data about parameter values and associated flow system features and can lead to improved predictions. Methods for identifying parameters important to predictions include prediction scaled sensitivities (PSS), which account for uncertainty on individual parameters as well as prediction sensitivity to parameters, and a new “value of improved information” (VOII) method presented here, which includes the effects of parameter correlation in addition to individual parameter uncertainty and prediction sensitivity. In this work, the PSS and VOII methods are demonstrated and evaluated using a model of the Death Valley regional groundwater flow system. The predictions of interest are advective transport paths originating at sites of past underground nuclear testing. Results show that for two paths evaluated the most important parameters include a subset of five or six of the 23 defined model parameters. Some of the parameters identified as most important are associated with flow system attributes that do not lie in the immediate vicinity of the paths. Results also indicate that the PSS and VOII methods can identify different important parameters. Because the methods emphasize somewhat different criteria for parameter importance, it is suggested that parameters identified by both methods be carefully considered in subsequent data collection efforts aimed at improving model predictions.

California, Nevada↗
USGS2003

Direct assessment of groundwater vulnerability from single observations of multiple contaminants

Groundwater vulnerability is a central concept in pollution risk assessment, yet its estimation has been largely a matter of expert judgment. This work applies a method for the direct calculation of vulnerability from monitoring well observations of pesticide concentrations. The method has two major advantages: it is independent of the compounds being examined, and it has a direct probabilistic interpretation making it ideal for risk assessment. The methodology was applied to data from a groundwater monitoring program in the midwestern United States. The distribution of the vulnerabilities was skewed toward zero. Spatial distribution of the vulnerabilities shows them to be controlled by both regional and local factors. Methods are presented for estimating the necessary sample sizes for vulnerability studies. The further application of the approach developed in this study to understanding groundwater pollution is discussed.

Midwest region↗
USGS2003

Use of hydraulic head to estimate volumetric gas content and ebullition flux in northern peatlands

Hydraulic head was overpressured at middepth in a 4.2‐m thick raised bog in the Glacial Lake Agassiz peatlands of northern Minnesota, and fluctuated in response to atmospheric pressure. Barometric efficiency (BE), determined by calculating ratios of change in hydraulic head to change in atmospheric pressure, ranged from 0.05 to 0.15 during July through November of both 1997 and 1998. The overpressuring and a BE response were caused by free‐phase gas contained primarily in the center of the peat column between two or more semielastic, semiconfining layers of more competent peat. Two methods were used to determine the volume of gas bubbles contained in the peat, one using the degree of overpressuring in the middepth of the peat, and the other relating BE to specific yield of the shallow peat. The volume of gas calculated from the overpressuring method averaged 9%, assuming that the gas was distributed over a 2‐m thick overpressured interval. The volume of gas using the BE method averaged 13%. Temporal changes in overpressuring and in BE indicate that the volume of gaseous‐phase gas also changed with time, most likely because of rapid degassing (ebullition) that allowed sudden loss of gas to the atmosphere. Estimates of gas released during the largest ebullition events ranged from 0.3 to 0.7 mol m −2 d −1 . These ebullition events may contribute a significant source of methane and carbon dioxide to the atmosphere that has so far largely gone unmeasured by gas‐flux chambers or tower‐mounted sensors.

Water Resources Research↗
USGS2003

Pumping‐induced leakage in a bounded aquifer: An example of a scale‐invariant phenomenon

A new approach is presented for calculation of the volume of pumping‐induced leakage entering an aquifer as a function of time. This approach simplifies the total leakage calculation by extending analytical‐based methods developed for infinite systems to bounded aquifers of any size. The simplification is possible because of the relationship between drawdown and leakage in aquifers laterally bounded by impermeable formations. This relationship produces a scale‐invariant total leakage; i.e., the volume of leakage as a function of time does not change with the size of the aquifer or with the location of the pumping well. Two examples and image well theory are used to demonstrate and prove, respectively, the generality of this interesting phenomenon.

Water Resources Research↗
USGS2003

Influence of transitional volcanic strata on lateral diversion at Yucca Mountain, Nevada

Natural hydraulic barriers exist at Yucca Mountain, Nevada, a potential high‐level nuclear waste repository, that have been identified as possible lateral diversions for reducing deep percolation through the waste storage area. Historical development of the conceptual model of lateral diversion has been limited by available field data, but numerical investigations presented the possibility of significant lateral diversion due to the presence of a thin, porous rock layer, the Paintbrush nonwelded tuffs. Analytical analyses of the influence of transitional changes in properties suggest that minimal lateral diversion is likely at Yucca Mountain. Numerical models, to this point, have not accounted for the gradual transition of properties or the existence of multiple layers that could inadvertently influence the simulation of lateral diversion as an artifact of numerical model discretization. Analyses were made of subsurface matric potential measurements, and comparisons were made of surface infiltration estimates with deeper percolation flux calculations using chloride‐mass‐balance calculations and simulations of measured temperature profiles. These analyses suggest that insignificant lateral diversion has occurred above the repository horizon and that water generally moves vertically through the Paintbrush nonwelded tuffs.

Water Resources Research↗
USGS2003

Diel cycles in dissolved metal concentrations in streams: Occurrence and possible causes

Substantial diel (24‐hour) cycles in dissolved (0.1‐μm filtration) metal concentrations were observed during low flow for 18 sampling episodes at 14 sites on 12 neutral and alkaline streams draining historical mining areas in Montana and Idaho. At some sites, concentrations of Cd, Mn, Ni, and Zn increased as much as 119, 306, 167, and 500%, respectively, from afternoon minimum values to maximum values shortly after sunrise. Arsenic concentrations exhibited the inverse temporal pattern with increases of up to 54%. Variations in Cu concentrations were small and inconsistent. Diel metal cycles are widespread and persistent, occur over a wide range of metal concentrations, and likely are caused primarily by instream geochemical processes. Adsorption is the only process that can explain the inverse temporal patterns of As and the divalent metals. Diel metal cycles have important implications for many types of water‐quality studies and for understanding trace‐metal mobility.

Water Resources Research↗
USGS2003

Variations in flow and transport in thick desert vadose zones in response to paleoclimatic forcing (0-90 kyr): Field measurements, modeling, and uncertainties

An understanding of unsaturated flow and potential recharge in interdrainage semiarid and arid regions is critical for quantification of water resources and contaminant transport. We evaluated system response to paleoclimatic forcing using water potential and Cl profiles and modeling of nonisothermal liquid and vapor flow and Cl transport at semiarid (High Plains, Texas) and arid (Chihuahuan Desert, Texas; Amargosa Desert, Nevada) sites. Infiltration in response to current climatic forcing is restricted to the shallow (∼0.3–3 m) subsurface. Subsurface Cl accumulations correspond to time periods of 9–90 kyr. Bulge-shaped Cl profiles generally represent accumulation during the Holocene (9–16 kyr). Lower Cl concentrations at depth reflect higher water fluxes (0.04–8.4 mm/yr) during the Pleistocene and earlier times. Low water potentials and upward gradients indicate current drying conditions. Nonisothermal liquid and vapor flow simulations indicate that upward flow for at least 1–2 kyr in the High Plains and for 12–16 kyr at the Chihuahuan and Amargosa desert sites is required to reproduce measured upward water potential gradients and that recharge is negligible (<0.1 mm/yr) in these interdrainage areas.

Water Resources Research↗
USGS2003

Stochastic analysis of transverse dispersion in density‐coupled transport in aquifers

Spectral perturbation techniques have been used previously to derive integral expressions for dispersive mixing in concentration‐dependent transport in three‐dimensional, heterogeneous porous media, where fluid density and viscosity are functions of solute concentration. Whereas earlier work focused on evaluating longitudinal dispersivity in isotropic media and incorporating the result in a mean one‐dimensional transport model, the emphasis of this paper is on evaluation of the complete dispersion tensor, including the more general case of anisotropic media. Approximate analytic expressions for all components of the macroscopic dispersivity tensor are derived, and the tensor is shown to be asymmetric. The tensor is separated into its symmetric and antisymmetric parts, where the symmetric part is used to calculate the principal components and principal directions of dispersivity, and the antisymmetric part of the tensor is shown to modify the velocity of the solute body compared to that of the background fluid. An example set of numerical simulations incorporating the tensor illustrates the effect of density‐coupled dispersivity on a sinking plume in an aquifer. The simulations show that the effective transverse vertical spreading in a sinking plume to be significantly greater than would be predicted by a standard density‐coupled transport model that does not incorporate the coupling in the dispersivity tensor.

Water Resources Research↗
USGS2003

Effect of isolated fractures on accelerated flow in unsaturated porous rock

Fractures that begin and end in the unsaturated zone, or isolated fractures, have been ignored in previous studies because they were generally assumed to behave as capillary barriers and remain nonconductive. We conducted a series of experiments using Berea sandstone samples to examine the physical mechanisms controlling flow in a rock containing a single isolated fracture. The input fluxes and fracture orientation were varied in these experiments. Visualization experiments using dyed water in a thin vertical slab of rock were conducted to identify flow mechanisms occurring due to the presence of the isolated fracture. Two mechanisms occurred: (1) localized flow through the rock matrix in the vicinity of the isolated fracture and (2) pooling of water at the bottom of the fracture, indicating the occurrence of film flow along the isolated fracture wall. These mechanisms were observed at fracture angles of 20 and 60 degrees from the horizontal, but not at 90 degrees. Pooling along the bottom of the fracture was observed over a wider range of input fluxes for low‐angled isolated fractures compared to high‐angled ones. Measurements of matrix water pressures in the samples with the 20 and 60 degree fractures also demonstrated that preferential flow occurred through the matrix in the fracture vicinity, where higher pressures occurred in the regions where faster flow was observed in the visualization experiments. The pooling length at the terminus of a 20 degree isolated fracture was measured as a function of input flux. Calculations of the film flow rate along the fracture were made using these measurements and indicated that up to 22% of the flow occurred as film flow. These experiments, apparently the first to consider isolated fractures, demonstrate that such features can accelerate flow through the unsaturated zone and should be considered when developing conceptual models.

Water Resources Research↗