USGS Science⌕ Search

SEARCH · USGS Science

Results for “Water 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 775 records · Page 43Linked to original sources
USGS2002

Partial entrainment of gravel bars during floods

Spatial patterns of bed material entrainment by floods were documented at seven gravel bars using arrays of metal washers (bed tags) placed in the streambed. The observed patterns were used to test a general stochastic model that bed material entrainment is a spatially independent, random process where the probability of entrainment is uniform over a gravel bar and a function of the peak dimensionless shear stress τ 0 * of the flood. The fraction of tags missing from a gravel bar during a flood, or partial entrainment, had an approximately normal distribution with respect to τ 0 * with a mean value (50% of the tags entrained) of 0.085 and standard deviation of 0.022 (root‐mean‐square error of 0.09). Variation in partial entrainment for a given τ 0 * demonstrated the effects of flow conditioning on bed strength, with lower values of partial entrainment after intermediate magnitude floods (0.065 < τ 0 *< 0.08) than after higher magnitude floods. Although the probability of bed material entrainment was approximately uniform over a gravel bar during individual floods and independent from flood to flood, regions of preferential stability and instability emerged at some bars over the course of a wet season. Deviations from spatially uniform and independent bed material entrainment were most pronounced for reaches with varied flow and in consecutive floods with small to intermediate magnitudes.

Water Resources Research↗
USGS2002

Weathering reactions and hyporheic exchange controls on stream water chemistry in a glacial meltwater stream in the McMurdo Dry Valleys

In the McMurdo Dry Valleys, Antarctica, dilute glacial meltwater flows down well‐established streambeds to closed basin lakes during the austral summer. During the 6–12 week flow season, a hyporheic zone develops in the saturated sediment adjacent to the streams. Longer Dry Valley streams have higher concentrations of major ions than shorter streams. The longitudinal increases in Si and K suggest that primary weathering contributes to the downstream solute increase. The hypothesis that weathering reactions in the hyporheic zone control stream chemistry was tested by modeling the downstream increase in solute concentration in von Guerard Stream in Taylor Valley. The average rates of solute supplied from these sources over the 5.2 km length of the stream were 6.1 × 10 −9 mol Si L −1 m −1 and 3.7 × 10 −9 mol K L −1 m −1 , yielding annual dissolved Si loads of 0.02–1.30 mol Si m −2 of watershed land surface. Silicate minerals in streambed sediment were analyzed to determine the representative surface area of minerals in the hyporheic zone subject to primary weathering. Two strategies were evaluated to compute sediment surface area normalized weathering rates. The first applies a best linear fit to synoptic data in order to calculate a constant downstream solute concentration gradient, dC / dx (constant weathering rate contribution, CRC method); the second uses a transient storage model to simulate dC / dx , representing both hyporheic exchange and chemical weathering (hydrologic exchange, HE method). Geometric surface area normalized dissolution rates of the silicate minerals in the stream ranged from 0.6 × 10 −12 mol Si m −2 s −1 to 4.5 × 10 −12 mol Si m −2 s −1 and 0.4 × 10 −12 mol K m −2 s −1 to 1.9 × 10 −12 mol K m −2 s −1 . These values are an order of magnitude lower than geometric surface area normalized weathering rates determined in laboratory studies and are an order of magnitude greater than geometric surface area normalized weathering rates determined in a warmer, wetter setting in temperate basins, despite the cold temperatures, lack of precipitation and lack of organic material. These results suggest that the continuous saturation and rapid flushing of the sediment due to hyporheic exchange facilitates weathering in Dry Valley streams.

Water Resources Research↗
USGS2002

Estimating the sources and transport of nutrients in the Waikato River Basin, New Zealand

We calibrated SPARROW (Spatially Referenced Regression on Watershed Attributes) surface water‐quality models using measurements of total nitrogen and total phosphorus from 37 sites in the 13,900‐km 2 Waikato River Basin, the largest watershed on the North Island of New Zealand. This first application of SPARROW outside of the United States included watersheds representative of a wide range of natural and cultural conditions and water‐resources data that were well suited for calibrating and validating the models. We applied the spatially distributed model to a drainage network of nearly 5000 stream reaches and 75 lakes and reservoirs to empirically estimate the rates of nutrient delivery (and their levels of uncertainty) from point and diffuse sources to streams, lakes, and watershed outlets. The resulting models displayed relatively small errors; predictions of stream yield (kg ha −1 yr −1 ) were typically within 30% or less of the observed values at the monitoring sites. There was strong evidence of the accuracy of the model estimates of nutrient sources and the natural rates of nutrient attenuation in surface waters. Estimated loss rates for streams, lakes, and reservoirs agreed closely with experimental measurements and empirical models from New Zealand, North America, and Europe as well as with previous U.S. SPARROW models. The results indicate that the SPARROW modeling technique provides a reliable method for relating experimental data and observations from small catchments to the transport of nutrients in the surface waters of large river basins.

Waikato River↗
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↗