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 1,189 records · Page 66Linked to original sources
USGS2008

Evaluation of unsaturated zone air permeability through pneumatic tests

Predicting the steady state distribution of air pressure in the unsaturated zone resulting from a pneumatic test provides a method for determining air-phase permeability. This technique is analogous to the inverse problem of well hydraulics; however, air flow is more complicated than ground water flow because of air compressibility, the Klinkenberg effect, variations in air density and viscosity that result from temperature fluctuations in the unsaturated zone and the possibility of inducing water movement during the pneumatic test. An analysis of these complicating factors reveals that, when induced water movement can be neglected, a linear version of the airflow equation can provide an appropriate approximation for the purpose of determining air-phase permeability. Two analytical solutions for steady state, two-dimensional, axisymmetric airflow to a single well partially screened in the unsaturated zone are developed. One solution applies where there is a stratum of relatively low air permeability, separating the stratum in which the well is completed, from the atmosphere. The other solution applies where there is no separating stratum between the domain and atmosphere. In both situations the water table forms the lower horizontal boundary. Applications of both solutions to determine air permeability from data collected during pneumatic tests are presented.

Water Resources Research↗
USGS2010

Arrival times and temporal moments of breakthrough curves for an imperfectly stratified aquifer

Using concepts developed in an earlier study, a solution in Laplace transform space is obtained for transport of resident concentration in an imperfectly but yet highly stratified porous medium. The flow field, into which an instantaneous pulse of tracer is injected, is taken to be steady and mean uniform parallel to the direction of stratification. From this transform-space solution either temporal moments can be derived by taking derivatives with respect to the Laplace parameter, or the transform-space solution can be inverted numerically to obtain breakthrough curves for the mean concentration. When compared to an equivalent solution with a Fickian dispersive flux, these temporal moments indicate the extent to which transport in heterogeneous porous media deviates from classical Fickian behavior. The numerical inversion of the Laplace transform solution gives partial breakthrough curves for the mean concentration which have the appearance of conflicting with the derived moment information. A hypothesis is put forth which resolves this apparent conflict; this hypothesis is verified by adding a component of local dispersion to the governing transport equation. On the basis of the flux-averaged concentration a form for the expected probability density function for the arrival time of a tracer particle is derived; arrival time moments and an arrival time cumulative distribution function are available as a consequence. Arrival time moments, as derived from the flux-averaged concentration, do not differ significantly from the resident moments, as derived from the resident concentration.

Water Resources Research↗
USGS2008

Convergent radial dispersion: A note on evaluation of the Laplace transform solution

A numerical inversion algorithm for Laplace transforms that is capable of handling rapid changes in the computed function is applied to the Laplace transform solution to the problem of convergent radial dispersion in a homogeneous aquifer. Prior attempts by the author to invert this solution were unsuccessful for highly advective systems where the Peclet number was relatively large. The algorithm used in this note allows for rapid and accurate inversion of the solution for all Peclet numbers of practical interest, and beyond. Dimensionless breakthrough curves are illustrated for tracer input in the form of a step function, a Dirac impulse, or a rectangular input.

Water Resources Research↗
USGS2010

Gravity-driven groundwater flow and slope failure potential: 1. Elastic effective-stress model

Hilly or mountainous topography influences gravity-driven groundwater flow and the consequent distribution of effective stress in shallow subsurface environments. Effective stress, in turn, influences the potential for slope failure. To evaluate these influences, we formulate a two-dimensional, steady state, poroelastic model. The governing equations incorporate groundwater effects as body forces, and they demonstrate that spatially uniform pore pressure changes do not influence effective stresses. We implement the model using two finite element codes. As an illustrative case, we calculate the groundwater flow field, total body force field, and effective stress field in a straight, homogeneous hillslope. The total body force and effective stress fields show that groundwater flow can influence shear stresses as well as effective normal stresses. In most parts of the hillslope, groundwater flow significantly increases the Coulomb failure potential Φ, which we define as the ratio of maximum shear stress to mean effective normal stress. Groundwater flow also shifts the locus of greatest failure potential toward the slope toe. However, the effects of groundwater flow on failure potential are less pronounced than might be anticipated on the basis of a simpler, one-dimensional, limit equilibrium analysis. This is a consequence of continuity, compatibility, and boundary constraints on the two-dimensional flow and stress fields, and it points to important differences between our elastic continuum model and limit equilibrium models commonly used to assess slope stability.

Water Resources Research↗
USGS2010

Gravity-driven groundwater flow and slope failure potential: 2. Effects of slope morphology, material properties, and hydraulic heterogeneity

Hillslope morphology, material properties, and hydraulic heterogeneities influence the role of groundwater flow in provoking slope instability. We evaluate these influences quantitatively by employing the elastic effective stress model and Coulomb failure potential concept described in our companion paper (Iverson and Reid, this issue). Sensitivity analyses show that of four dimensionless quantities that control model results (i.e., Poisson's ratio, porosity, topographic profile, and hydraulic conductivity contrast), slope profiles and hydraulic conductivity contrasts have the most pronounced and diverse effects on groundwater seepage forces, effective stresses, and slope failure potentials. Gravity-driven groundwater flow strongly influences the shape of equilibrium hillslopes, which we define as those with uniform near-surface failure potentials. For homogeneous slopes with no groundwater flow, equilibrium hillslope profiles are straight; but with gravity-driven flow, equilibrium profiles are concave or convex-concave, and the largest failure potentials exist near the bases of convex slopes. In heterogeneous slopes, relatively slight hydraulic conductivity contrasts of less than 1 order of magnitude markedly affect the seepage force field and slope failure potential. Maximum effects occur if conductivity contrasts are of four orders of magnitude or more, and large hydraulic gradients commonly result in particularly large failure potentials just upslope from where low-conductivity layers intersect the ground surface.

Water Resources Research↗
USGS2010

Debris flow rheology: Experimental analysis of fine-grained slurries

The rheology of slurries consisting of ≤2-mm sediment from a natural debris flow deposit was measured using a wide-gap concentric-cylinder viscometer. The influence of sediment concentration and size and distribution of grains on the bulk rheological behavior of the slurries was evaluated at concentrations ranging from 0.44 to 0.66. The slurries exhibit diverse rheological behavior. At shear rates above 5 s −1 the behavior approaches that of a Bingham material; below 5 s −1 , sand exerts more influence and slurry behavior deviates from the Bingham idealization. Sand grain interactions dominate the mechanical behavior when sand concentration exceeds 0.2; transient fluctuations in measured torque, time-dependent decay of torque, and hysteresis effects are observed. Grain rubbing, interlocking, and collision cause changes in packing density, particle distribution, grain orientation, and formation and destruction of grain clusters, which may explain the observed behavior. Yield strength and plastic viscosity exhibit order-of-magnitude variation when sediment concentration changes as little as 2–4%. Owing to these complexities, it is unlikely that debris flows can be characterized by a single rheological model.

Water Resources Research↗
USGS2010

Flow resistance under conditions of intense gravel transport

A study of flow resistance was undertaken in a channelized reach of the North Fork Toutle River, downstream of Mount St. Helens, Washington. Hydraulic and sediment transport data were collected in flows with velocities up to 3 m/s and shear stresses up to 7 times the critical value needed for bed load transport. Details of the flow structure as revealed in vertical velocity profiles indicate that weak bed load transport over a plane gravel bed has little effect on flow resistance. The plane gravel bed persists up to stresses ∼3 times critical, at which point, irregular bed forms appear. Bed forms greatly increase flow resistance and cause velocity profiles to become distorted. The latter arises as an effect of flows becoming depth-limited as bed form amplitude increases. At very high rates of bed load transport, an upper stage plane bed appeared. Velocity profiles measured in these flows match the law of the wall closely, with the equivalent roughness being well represented by k s = 3 D 84 of the bed load. The effects noted here will be important in very large floods or in rivers that are not free to widen, such as those cut into bedrock.

Washington↗
USGS2010

A solute flux approach to transport in heterogeneous formations: 2. Uncertainty analysis

Uncertainty in the mass flux for advection dominated solute movement in heterogeneous porous media is investigated using the Lagrangian framework developed in paper 1 by Dagan et al. (this issue). Expressions for the covariance of the mass flux and cumulative mass flux are derived as functions of the injection volume and sampling area size relative to the scale of heterogeneity. The result is illustrated for solute advection in three types of heterogeneous porous media: stratified formations, two‐ and three‐dimensional porous media; small perturbation approximation is used for the two‐ and three‐dimensional cases. Variances of the mass flux and cumulative mass flux are evaluated as functions of the injection volume (area) scale versus log‐hydraulic conductivity integral scale. The greatest decrease in coefficient of variation (CV) of the mass flux is for the source scale 1–5 times the hydraulic conductivity integral scale; further increase in the source size decreases CV comparatively less. The variance of the cumulative mass flux (or total discharge) indicates that for the source size of 20 hydraulic conductivity integral scales, the transport conditions are almost ergodic. The present results also indicate that the cumulative mass flux is a relatively robust quantity for describing field‐scale solute transport.

Water Resources Research↗
USGS2010

A solute flux approach to transport in heterogeneous formations: 1. The general framework

It is common to represent solute tranport in heterogeneous formations in terms of the resident concentration C (x, t), regarded as a random space function. The present study investigates the alternative representation by q , the solute mass flux at a point of a control plane normal to the mean flow. This representation is appropriate for many field applications in which the variable of interest is the mass of solute discharged through a control surface. A general framework to compute the statistical moments of q and of the associated total solute discharge Q and mass M is established. With x the direction of the mean flow, a solute particle is crossing the control plane at y = η, z = ζ and at the travel (arrival) time τ. The associated expected solute flux value is proportional to the joint probability density function (pdf) g 1 of η, ζ and τ, whereas the variance of q is shown to depend on the joint pdf g 2 of the same variables for two particles. In turn, the statistical moments of η, ζ and τ depend on those of the velocity components through a system of stochastic ordinary differential equations. For a steady velocity field and neglecting the effect of pore‐scale dispersion, a major simplification of the problem results in the independence of the random variables η, ζ and τ. As a consequence, the pdf of η and ζ can be derived independently of τ. A few approximate approaches to derive the statistical moments of η, ζ and τ are outlined. These methods will be explored in paper 2 in order to effectively derive the variances of the total solute discharge and mass, while paper 3 will deal with the nonlinear effect of the velocity variance upon the moments of η, ζ and τ

Water Resources Research↗
USGS2010

Legal mechanisms for protecting riparian resource values

Riparian resources include the borders of rivers, lakes, ponds, and potholes. These border areas are very important for a number of reasons, including stream channel maintenance, flood control, aesthetics, erosion control, fish and wildlife habitat, recreation, and water quality maintenance. These diverse functions are not well protected by law or policy. We reviewed law and policies regarding endangered species habitat designation, land use planning, grazing management, water allocation, takings, and federal permits and licenses, along with the roles of federal, state, and local governments. We discuss the politics of implementing these policies, focusing on the difficulties in changing entrenched water and land use practices. Our review indicates a lack of direct attention to riparian ecosystem issues in almost all environmental and land use programs at every level of government. Protection of riparian resource values requires a means to integrate existing programs to focus on riparian zones.

Water Resources Research↗
USGS2010

Solute transport with multisegment, equilibrium-controlled, classical reactions: Problem solvability and feed forward method's applicability for complex segments of at most binary participants

The feed forward (FF) method derives efficient operational equations for simulating transport of reacting solutes. It has been shown to be applicable in the presence of networks with any number of homogeneous and/or heterogeneous, classical reaction segments that consist of three, at most binary participants. Using a sequential (network type after network type) exploration approach and, independently, theoretical explanations, it is demonstrated for networks with classical reaction segments containing more than three, at most binary participants that if any one of such networks leads to a solvable transport problem then the FF method is applicable. Ways of helping to avoid networks that produce problem insolvability are developed and demonstrated. A previously suggested algebraic, matrix rank procedure has been adapted and augmented to serve as the main, easy-to-apply solvability test for already postulated networks. Four network conditions that often generate insolvability have been identified and studied. Their early detection during network formulation may help to avoid postulation of insolvable networks.

Water Resources Research↗
USGS2010

Large-scale natural gradient tracer test in sand and gravel, Cape Cod, Massachusetts: 3. Hydraulic conductivity variability and calculated macrodispersivities

Hydraulic conductivity ( K ) variability in a sand and gravel aquifer on Cape Cod, Massachusetts, was measured and subsequently used in stochastic transport theories to estimate macrodispersivities. Nearly 1500 K measurements were obtained by borehole flowmeter tests and permeameter analyses of cores. The geometric mean for the flowmeter tests (0.11 cm/s) is similar to that estimated from other field tests. The mean for the permeameter tests (0.035 cm/s) is significantly lower, possibly because of compaction of the cores. The variance for the flowmeter (0.24) is also greater than that for the permeameter (0.14). Geostatistical analyses applying negative exponential models with and without nuggets reveal similar spatial correlation structures for the two data sets. Estimated correlation scales range from 2.9 to 8 m in the horizontal and from 0.18 to 0.38 m in the vertical. Estimates of asymptotic longitudinal dispersivity (b.35–0.78 m) are similar in magnitude to that observed in the natural gradient tracer test (0.96 m) previously conducted at this site.

Massachusetts↗
USGS2010

Solute transport with multiple equilibrium-controlled or kinetically controlled chemical reactions

A new approach is applied to the problem of modeling solute transport accompanied by many chemical reactions. The approach, based on concepts of the concentration space and its stoichiometric subspaces, uses elements of the subspaces as primary dependent variables. It is shown that the resulting model equations are compact in form, isolate the chemical reaction expressions from flow expressions, and can be used for either equilibrium or kinetically controlled reactions. The implications of the results on numerical algorithms for solving the equations are discussed. The application of the theory is illustrated throughout with examples involving a simple but broadly representative set of reactions previously considered in the literature. Numerical results are presented for four interconnected reactions: a homogeneous complexation reaction, two sorption reactions, and a dissolution/precipitation reaction. Three cases are considered: (1) four kinetically controlled reactions, (2) four equilibrium-controlled reactions, and (3) a system with two kinetically controlled reactions and two equilibrium-controlled reactions.

Water Resources Research↗
USGS2010

Relation of nickel concentrations in tree rings to groundwater contamination

Increment cores were collected from trees growing at two sites where groundwater is contaminated by nickel. Proton-induced X ray emission spectroscopy was used to determine the nickel concentrations in selected individual rings and in parts of individual rings. Ring nickel concentrations were interpreted on the basis of recent concentrations of nickel in aquifers, historical information about site use activities, and model simulations of groundwater flow. Nickel concentrations in rings increased during years of site use but not in trees outside the contaminated aquifers. Consequently, it was concluded that trees may preserve in their rings an annual record of nickel contamination in groundwater. Tulip trees and oaks contained higher concentrations of nickel than did sassafras, sweet gum, or black cherry. No evidence was found that nickel accumulates consistently within parts of individual rings or that nickel is translocated across ring boundaries.

Water Resources Research↗
USGS2010

Fat fractal scaling of drainage networks from a random spatial network model

An alternative quantification of the scaling properties of river channel networks is explored using a spatial network model. Whereas scaling descriptions of drainage networks previously have been presented using a fractal analysis primarily of the channel lengths, we illustrate the scaling of the surface area of the channels defining the network pattern with an exponent which is independent of the fractal dimension but not of the fractal nature of the network. The methodology presented is a fat fractal analysis in which the drainage basin minus the channel area is considered the fat fractal. Random channel networks within a fixed basin area are generated on grids of different scales. The sample channel networks generated by the model have a common outlet of fixed width and a rule of upstream channel narrowing specified by a diameter branching exponent using hydraulic and geomorphologic principles. Scaling exponents are computed for each sample network on a given grid size and are regressed against network magnitude. Results indicate that the size of the exponents are related to magnitude of the networks and generally decrease as network magnitude increases. Cases showing differences in scaling exponents with like magnitudes suggest a direction of future work regarding other topologic basin characteristics as potential explanatory variables.

Water Resources Research↗
USGS2010

Use of chlorofluorocarbons (CCl3F and CCl2F2) as hydrologic tracers and age‐dating tools: The alluvium and terrace system of central Oklahoma

The use of Chlorofluorocarbons (CFCs) as an age‐dating tool and tracer in shallow groundwaters has been investigated. New methodology for field sampling and preserving groundwaters containing parts per trillion concentrations of the CFCs, F‐1l and F‐12, is presented. Samples are analyzed by purge‐and‐trap gas chromatography with an electron capture detector. Physical and chemical processes that can alter natural concentrations (air‐water equilibrium) of CFCs were investigated to assess dating uncertainties. CFC model recharge ages appear to be defined within 2 years under optimum conditions. The method was applied to central Oklahoma to demonstrate the usefulness of CFCs as (1) an age‐dating tool of shallow groundwaters, (2) a tracer of sewage effluent in surface and shallow groundwaters, and (3) a tracer of shallow groundwater. Results of dating indicate two primary recharge periods in central Oklahoma over the past 45 years that correspond to the wet periods 1945–1960 and 1967–1975.

Water Resources Research↗
USGS2010

A geochemical transport model for redox-controlled movement of mineral fronts in groundwater flow systems: A case of nitrate removal by oxidation of pyrite

A one-dimensional prototype geochemical transport model was developed in order to handle simultaneous precipitation-dissolution and oxidation-reduction reactions governed by chemical equilibria. Total aqueous component concentrations are the primary dependent variables, and a sequential iterative approach is used for the calculation. The model was verified by analytical and numerical comparisons and is able to simulate sharp mineral fronts. At a site in Denmark, denitrification has been observed by oxidation of pyrite. Simulation of nitrate movement at this site showed a redox front movement rate of 0.58 m yr −1 , which agreed with calculations of others. It appears that the sequential iterative approach is the most practical for extension to multidimensional simulation and for handling large numbers of components and reactions. However, slow convergence may limit the size of redox systems that can be handled.

Water Resources Research↗
USGS2010

Evolution of coarse gravel bed forms: Field measurements at flood stage

Field measurements to investigate the origin and growth of mesoscale gravel bed forms in deep flows were made in the North Fork Toutle River, Washington. Sonar observations of the gravel streambed at a stationary point were recorded during two storm flows in December 1989 and January 1990 with concurrent bed load sampling and continuous velocity measurements. Mean diameter of bed load was about 3 cm, flow depths were 1.4–2.4 m, and bed shear stresses were 2–5 times the critical stress of mean bed load diameter, as computed from the depth-slope product. These records document the hydrodynamic conditions under which dunelike coarse gravel bed forms were observed. Coarse gravel dunes (height, 20 cm; length, 6–15 m) evolved more than 24 hours after peak stage, primarily by accretion, as inferred from bed form changes revealed in dual sonar records. Dune heights increased to 40 cm as mean trough elevation rose about 50 cm over several hours. Smaller dunes (wavelength, 1–3 m), transitional from bed load sheets, migrated on the backs of the large dunes. The superposed dunes finally became indistinguishable from the large dunes, which diminished in height by increasing the mean level of troughs. Gravel deposition occurred at the observation point in conjunction with migration of gravel dunes. The direct comparison of known bed form regimes and gravel bar facies provides alternative interpretations of gravelly deposits.

Water Resources Research↗