Search USGS⌕ Search

SEARCH · Search USGS

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

Errors and parameter estimation in precipitation-runoff modeling: 2. Case study

A case study is presented which illustrates some of the error analysis, sensitivity analysis, and parameter estimation procedures reviewed in the first part of this paper. It is shown that those procedures, most of which come from statistical nonlinear regression theory, are invaluable in interpreting errors in precipitation-runoff modeling and in identifying appropriate calibration strategies.

Water Resources Research↗
USGS2008

Type curve analysis of inertial effects in the response of a well to a slug test.

The water level response to a slug or bailer test in a well completed in a confined aquifer has been evaluated taking into account well-bore storage and inertial effects of the water column in the well. The response range, from overdamped with negligible inertial effects to damped oscillation, was covered employing numerical inversions of the Laplace-transform solution. By scaling the time with respect to the undamped natural period of the well-aquifer system and by using the damping parameter for a second-order damped, inertial-elastic system, a set of type curves was constructed that enables water level response data from a slug or bailer test to be analyzed under conditions where the inertial parameter is large. Values of transmissivity and effective static water column length can be determined when an estimate of storage coefficient is available. The numerical solution and resulting type curves cover the transition range between the limiting cases of negligible inertial effects and of damped oscillation that have been treated by others. Two examples of slug test analysis show that precise results depend on accurate measurements of water level displacement (±1% of initial value).

Water Resources Research↗
USGS2008

A comparison of several methods of solving nonlinear regression groundwater flow problems

Computational efficiency and computer memory requirements for four methods of minimizing functions were compared for four test nonlinear-regression steady state groundwater flow problems. The fastest methods were the Marquardt and quasi-linearization methods, which required almost identical computer times and numbers of iterations; the next fastest was the quasi-Newton method, and last was the Fletcher-Reeves method, which did not converge in 100 iterations for two of the problems. The fastest method per iteration was the Fletcher-Reeves method, and this was followed closely by the quasi-Newton method. The Marquardt and quasi-linearization methods were slower. For all four methods the speed per iteration was directly related to the number of parameters in the model. However, this effect was much more pronounced for the Marquardt and quasi-linearization methods than for the other two. Hence the quasi-Newton (and perhaps Fletcher-Reeves) method might be more efficient than either the Marquardt or quasi-linearization methods if the number of parameters in a particular model were large, although this remains to be proven. The Marquardt method required somewhat less central memory than the quasi-linearization metilod for three of the four problems. For all four problems the quasi-Newton method required roughly two thirds to three quarters of the memory required by the Marquardt method, and the Fletcher-Reeves method required slightly less memory than the quasi-Newton method. Memory requirements were not excessive for any of the four methods.

Water Resources Research↗
USGS2008

Geochemistry of groundwater in Cretaceous sediments of the southeastern coastal plain of eastern Mississippi and western Alabama

Geochemical samples of waters along two hydrologic flow paths in four Upper Cretaceous aquifers of northeastern Mississippi and western Alabama indicate similar geochemical evolution of their respective waters. The waters of the Coker, Gordo, and Eutaw-McShan aquifers, noncalcareous sands, increase downgradient in dissolved solids and p H, and are dominated by sodium and bicarbonate ions, which generally result from a calcite dissolution-cation exchange process. Increases in dissolved iron from oxidation reduction reactions followed by decreases in total inorganic carbon from siderite precipitation occur along the flow paths. As the total inorganic carbon increases, carbon 13 (δ 13 C) generally is enriched in the moving waters, indicating the addition of a predominantly heavy source of carbon, most likely dissolving calcite. In the Coker aquifer δ 13 C values in the waters become more negative downgradient, resulting from lignite oxidation, followed by δ 13 C values becoming more positive, resulting from dissolving calcite and perhaps some mixing with brines. In northeastern Mississippi the Ripley aquifer, a calcareous sand, initially contains calcium-bicarbonate dominated water that evolves to a sodium- bicarbonate dominated water downgradient, primarily from the calcite dissolution-cation exchange process. Feldspar hydrolysis to kaolinite dominates aluminosilicate reactions in the upgradient parts of the aquifers. Authigenesis of smectite clay may be occurring in the deeper, downgradient parts of the aquifers.

Alabama, Mississippi↗
USGS2008

Assessment of long-term salinity changes in an irrigated stream-aquifer system

Changes in salinity in groundwater and surface water in the Arkansas River valley of southeastern Colorado are primarily related to irrigation practices. A solute transport model was applied to an 11-mile reach of the valley to compute salinity changes in response to spatially and temporally varying stresses. The model was calibrated in 1973 using detailed field measurements made during 1971 and 1972. In 1973 the calibrated model was used to predict that a gradual long-term increase in groundwater salinity of about 2–3% per year would occur if the observed irrigation practices continued. The study area was resampled during the winter of 1982 to help evaluate if any long-term changes in salinity are actually occurring. Nonparametric and parametric statistical tests were used to help assess the significance of observed changes in groundwater salinity. These tests indicate that a statistically significant increase in salinity occurred between the winters of 1971 and 1972 (the model calibration period). However, a comparison of the winter 1972 and winter 1982 data indicates that no significant net change in salinity has occurred during this 10-year period. An analysis of the few available historical data (1895, 1923, 1959–1961, and 1964) supports the hypothesis that groundwater salinity in this irrigated area has reached a long-term dynamic equilibrium in response to irrigation practices. The model predictions of long-term salinity increases were invalid probably because the calibration period occurred during a short-term annual trend of increasing salinity in the river (and hence in leaky irrigation canals and in applied irrigation water), which was not representative of the long-term trend.

Colorado↗
USGS2008

Field determination of the three-dimensional hydraulic conductivity tensor of anisotropic media: 1. Theory

A field method is proposed for determining the three-dimensional hydraulic conductivity tensor and specific storage of an anisotropic porous or fractured medium. The method, known as cross-hole testing (to distinguish it from conventional single-hole packer tests), consists of injecting fluid into (or withdrawing fluid from) packed-off intervals in a number of boreholes and monitoring the transient head response in similar intervals in neighboring boreholes. The directions of the principal hydraulic conductivities need not be known prior to the test, and the boreholes may have arbitrary orientations (e.g., they can all be vertical). An important aspect of the proposed method is that it provides direct field information on whether it is proper to regard the medium as being uniform and anisotropic on the scale of the test. The first paper presents theoretical expressions describing transient and steady state head response in monitoring intervals of arbitrary lengths and orientations, to constant-rate injection into (or withdrawal from) intervals having similar or different lengths and orientations. The conditions under which these intervals can be treated mathematically as points are investigated by an asymptotic analysis. The effect of planar no-flow and constant-head boundaries on the response is analyzed by the theory of images. The second paper describes the field methodology and shows how the proposed approach works in the case of fractured granitic rocks.

Water Resources Research↗
USGS2008

Field determination of the three-dimensional hydraulic conductivity tensor of anisotropic media: 2. Methodology and application to fractured rocks

The analytical solutions developed in the first paper can be used to interpret the results of cross-hole tests conducted in anisotropic porous or fractured media. In the particular case where the injection and monitoring intervals are short relative to the distance between them, the test results can be analyzed graphically. From the transient variation of hydraulic head in a given monitoring interval, one can determine the directional hydraulic diffusivity, K d ( e )/ S s , and the quantity D / S s , by curve matching. (Here K d ( e ) is directional hydraulic conductivity parallel to the unit vector, e , pointing from the injection to the monitoring interval, S s is specific storage, and D is the determinant of the hydraulic conductivity tensor, K .) The principal values and directions of K , together with S s , can then be evaluated by fitting an ellipsoid to the square roots of the directional diffusivities. Ideally, six directional measurements are required. In practice, a larger number of measurements is often necessary to enable fitting an ellipsoid to the data by least squares. If the computed [ K d ( e )/ s s ] ½ values fluctuate so severely that a meaningful least squares fit is not possible, one has a direct indication that the subsurface does not behave as a uniform anisotropic medium on the scale of the test. Test results from a granitic rock near Oracle in southern Arizona are presented to illustrate how the method works for fractured rocks. At the site, the Oracle granite is shown to respond as a near-uniform, anisotropic medium, the hydraulic conductivity of which is strongly controlled by the orientations of major fracture sets. The cross-hole test results are shown to be consistent with the results of more than 100 single-hole packer tests conducted at the site.

Water Resources Research↗
USGS2008

Assessment of the instantaneous unit hydrograph derived from the theory of topologically random networks

An instantaneous unit hydrograph (iuh) based on the theory of topologically random networks (topological iuh) is evaluated in terms of sets of basin characteristics and hydraulic parameters. Hydrographs were computed using two linear routing methods for each of two drainage basins in the southeastern United States and are the basis of comparison for the topological iuh's. Elements in the sets of basin characteristics for the topological iuh's are the number of first-order streams only, ( N ), or the nuber of sources together with the number of channel links in the topological diameter ( N, D ); the hydraulic parameters are values of the celerity and diffusivity constant. Sensitivity analyses indicate that the mean celerity of the internal links in the network is the critical hydraulic parameter for determining the shape of the topological iuh, while the diffusivity constant has minimal effect on the topological iuh. Asymptotic results (source-only) indicate the number of sources need not be large to approximate the topological iuh with the Weibull probability density function.

Water Resources Research↗
USGS1985

A policy evaluation tool: Management of a multiaquifer system using controlled stream recharge

A model for the optimal allocation of water resources was developed for a multiaquifer groundwater and surface water system near Livermore, California. The complex groundwater system was analyzed using a transient, quasi-three-dimensional model that considers the nonlinear behavior of the unconfined aquifer. The surface water system consists of a reservoir that discharges water to three streams which in turn recharge the upper aquifer. Nonlinear streamflow-recharge relationships were developed based upon synoptic field measurements of streamflow. The management model uses constrained optimization to minimize the cost of allocating surface water subject to physical and economic restrictions. Results indicate that a combined hydrologic and economic management model can be used to evaluate management practices of a complex hydrogeologic system. Questions can be posed which either would be impossible or extremely difficult to solve without the management model. We demonstrate the utility of such a model in three areas. First, the efficiency of intra-basin water allocations is evaluated. Second, critical factors that control management decisions of the basin are identified. Third, the influence of economic incentives that can best satisfy the conflicting objectives of various water users is explored.

Water Resources Research↗
USGS2010

Estimation of distributional parameters for censored trace level water quality data: 1. Estimation techniques

A recurring difficulty encountered in investigations of many metals and organic contaminants in ambient waters is that a substantial portion of water sample concentrations are below limits of detection established by analytical laboratories. Several methods were evaluated for estimating distributional parameters for such censored data sets using only uncensored observations. Their reliabilities were evaluated by a Monte Carlo experiment in which small samples were generated from a wide range of parent distributions and censored at varying levels. Eight methods were used to estimate the mean, standard deviation, median, and interquartile range. Criteria were developed, based on the distribution of uncensored observations, for determining the best performing parameter estimation method for any particular data set. The most robust method for minimizing error in censored-sample estimates of the four distributional parameters over all simulation conditions was the log-probability regression method. With this method, censored observations are assumed to follow the zero-to-censoring level portion of a lognormal distribution obtained by a least squares regression between logarithms of uncensored concentration observations and their z scores. When method performance was separately evaluated for each distributional parameter over all simulation conditions, the log-probability regression method still had the smallest errors for the mean and standard deviation, but the lognormal maximum likelihood method had the smallest errors for the median and interquartile range. When data sets were classified prior to parameter estimation into groups reflecting their probable parent distributions, the ranking of estimation methods was similar, but the accuracy of error estimates was markedly improved over those without classification.

Water Resources Research↗
USGS2010

Estimation of distributional parameters for censored trace level water quality data: 2. Verification and applications

Estimates of distributional parameters (mean, standard deviation, median, interquartile range) are often desired for data sets containing censored observations. Eight methods for estimating these parameters have been evaluated by R. J. Gilliom and D. R. Helsel (this issue) using Monte Carlo simulations. To verify those findings, the same methods are now applied to actual water quality data. The best method (lowest root-mean-squared error (rmse)) over all parameters, sample sizes, and censoring levels is log probability regression (LR), the method found best in the Monte Carlo simulations. Best methods for estimating moment or percentile parameters separately are also identical to the simulations. Reliability of these estimates can be expressed as confidence intervals using rmse and bias values taken from the simulation results. Finally, a new simulation study shows that best methods for estimating uncensored sample statistics from censored data sets are identical to those for estimating population parameters. Thus this study and the companion study by Gilliom and Helsel form the basis for making the best possible estimates of either population parameters or sample statistics from censored water quality data, and for assessments of their reliability.

Water Resources Research↗
USGS2010

Modeling contamination of shallow unconfined aquifers through infiltration beds

We model the transport of a simply reactive contaminant through an infiltration bed and underlying shallow, one-dimensional, unconfined aquifer with a plane, steeply sloping bottom in the assumed absence of dispersion and downgradient dilution. The effluent discharge and ambient groundwater flow under the infiltration beds are presumed to form a vertically mixed plume marked by an appreciable radial velocity component in the near field flow region. The near field analysis routes effluent contamination as a single linear reservoir whose output forms a source plane for the one-dimensional, far field flow region downgradient of the facility; the location and width of the source plane reflect the relative strengths of ambient flow and effluent discharge. We model far field contaminant transport, using an existing method of characteristics solution with frame speeds modified by recharge, bottom slope, and linear adsorption, and concentrations reflecting first-order reaction kinetics. The near and far field models simulate transport of synthetic detergents, chloride, total nitrogen, and boron in a contaminant plume at the Otis Air Force Base sewage treatment plant in Barnstable County, Massachusetts, with reasonable accuracy.

Massachusetts↗
USGS2010

Groundwater transport of strontium 90 in a glacial outwash environment

As part of the investigation of groundwater contamination at a uranium-scrap recovery plant at Wood River Junction, Rhode Island, laboratory experiments led to the development of a model for predicting the transport of strontium 90 in glacial outwash sediments based on an approximate mechanism for ion exchange. The multicomponent system was simplified to two components by regarding all exchangeable cations other than strontium 90 as a single component. The binary ion-exchange parameter was a function of the variable, total ion concentration. A one-dimensional solute transport model was formulated to evaluate the time necessary for natural groundwater flow to remove the strontium 90 contamination plume from the groundwater system to the Pawcatuck River. The finite difference transport equations were solved sequentially for total ion concentrations, then strontium 90 concentrations. Clay-free quartz and feldspar sands at the study site have little potential for strontium 90 sorption, and high calcium, magnesium, and sodium concentrations compete for the few ion exchange sites. As the total ion concentration plume moves out of the system, ion exchange of strontium 90 increases, reducing the strontium 90 concentration in the groundwater. Cleanout times predicted using the binary ion exchange mechanism were about two thirds of those predicted using a constant distribution coefficient. It is suggested that this type of model can simulate solute transport more realistically in many groundwater systems where the total ion concentration is not constant.

Rhode Island↗
USGS2010

Stochastic analysis of three-dimensional flow in a bounded domain

A commonly accepted first-order approximation of the equation for steady state flow in a fully saturated spatially random medium has the form of Poisson's equation. This form allows for the advantageous use of Green's functions to solve for the random output (hydraulic heads) in terms of a convolution over the random input (the logarithm of hydraulic conductivity). A solution for steady state three- dimensional flow in an aquifer bounded above and below is presented; consideration of these boundaries is made possible by use of Green's functions to solve Poisson's equation. Within the bounded domain the medium hydraulic conductivity is assumed to be a second-order stationary random process as represented by a simple three-dimensional covariance function. Upper and lower boundaries are taken to be no-flow boundaries; the mean flow vector lies entirely in the horizontal dimensions. The resulting hydraulic head covariance function exhibits nonstationary effects resulting from the imposition of boundary conditions. Comparisons are made with existing infinite domain solutions.

Water Resources Research↗
USGS2010

The problem of complex eigensystems in the semianalytical solution for advancement of time in solute transport simulations: a new method using real arithmetic

In the numerical modeling of groundwater solute transport, explicit solutions may be obtained for the concentration field at any future time without computing concentrations at intermediate times. The spatial variables are discretized and time is left continuous in the governing differential equation. These semianalytical solutions have been presented in the literature and involve the eigensystem of a coefficient matrix. This eigensystem may be complex (i.e., have imaginary components) due to the asymmetry created by the advection term in the governing advection-dispersion equation. Previous investigators have either used complex arithmetic to represent a complex eigensystem or chosen large dispersivity values for which the imaginary components of the complex eigenvalues may be ignored without significant error. It is shown here that the error due to ignoring the imaginary components of complex eigenvalues is large for small dispersivity values. A new algorithm that represents the complex eigensystem by converting it to a real eigensystem is presented. The method requires only real arithmetic.

Water Resources Research↗
USGS2010

Ground-water flow in low permeability environments

Certain geologic media are known to have small permeability; subsurface environments composed of these media and lacking well developed secondary permeability have groundwater flow sytems with many distinctive characteristics. Moreover, groundwater flow in these environments appears to influence the evolution of certain hydrologic, geologic, and geochemical systems, may affect the accumulation of pertroleum and ores, and probably has a role in the structural evolution of parts of the crust. Such environments are also important in the context of waste disposal. This review attempts to synthesize the diverse contributions of various disciplines to the problem of flow in low-permeability environments. Problems hindering analysis are enumerated together with suggested approaches to overcoming them. A common thread running through the discussion is the significance of size- and time-scale limitations of the ability to directly observe flow behavior and make measurements of parameters. These limitations have resulted in rather distinct small- and large-scale approaches to the problem. The first part of the review considers experimental investigations of low-permeability flow, including in situ testing; these are generally conducted on temporal and spatial scales which are relatively small compared with those of interest. Results from this work have provided increasingly detailed information about many aspects of the flow but leave certain questions unanswered. Recent advances in laboratory and in situ testing techniques have permitted measurements of permeability and storage properties in progressively “tighter” media and investigation of transient flow under these conditions. However, very large hydraulic gradients are still required for the tests; an observational gap exists for typical in situ gradients. The applicability of Darcy's law in this range is therefore untested, although claims of observed non-Darcian behavior appear flawed. Two important nonhydraulic flow phenomena, osmosis and ultrafiltration, are experimentally well established in prepared clays but have been incompletely investigated, particularly in undisturbed geologic media. Small-scale experimental results form much of the basis for analyses of flow in low-permeability environments which occurs on scales of time and size too large to permit direct observation. Such large-scale flow behavior is the focus of the second part of the review. Extrapolation of small-scale experimental experience becomes an important and sometimes controversial problem in this context. In large flow systems under steady state conditions the regional permeability can sometimes be determined, but systems with transient flow are more difficult to analyze. The complexity of the problem is enhanced by the sensitivity of large-scale flow to the effects of slow geologic processes. One-dimensional studies have begun to elucidate how simple burial or exhumation can generate transient flow conditions by changing the state of stress and temperature and by burial metamorphism. Investigation of the more complex problem of the interaction of geologic processes and flow in two and three dimensions is just beginning. Because these transient flow analyses have largely been based on flow in experimental scale systems or in relatively permeable systems, deformation in response to effective stress changes is generally treated as linearly elastic; however, this treatment creates difficulties for the long periods of interest because viscoelastic deformation is probably significant. Also, large-scale flow simulations in argillaceous environments generally have neglected osmosis and ultrafiltration, in part because extrapolation of laboratory experience with coupled flow to large scales under in situ conditions is controversial. Nevertheless, the effects are potentially quite important because the coupled flow might cause ultra long lived transient conditions. The difficulties associated with analysis are matched by those of characterizing hydrologic conditions in tight environments; measurements of hydraulic head and sampling of pore fluids have been done only rarely because of the practical difficulties involved. These problems are also discussed in the second part of this paper.

Water Resources Research↗
USGS2010

Groundwater Seepage Vectors and the Potential for Hillslope Failure and Debris Flow Mobilization

Insight for understanding the effect of groundwater flow on the potential for hillslope failure and liquefaction is provided by a novel limit‐equilibrium analysis of infinite slopes with steady, uniform Darcian seepage of arbitrary magnitude and direction. Normalization of the limit‐equilibrium solution shows that three dimensionless parameters govern completely the Coulomb failure potential of saturated, cohesionless, infinite homogeneous hillslopes: (1) the ratio of seepage force magnitude to gravitational body force magnitude; (2) the angle θ − Φ, where θ is the surface slope angle and Φ is the angle of internal friction of the soil; and (3) the angle λ + Φ, where λ is the angle of the seepage vector measured with respect to an outward‐directed surface‐normal vector. An additional dimensionless parameter affects the solution if soil cohesion is included in the analysis. Representation of the normalized solution as a single family of curves shows that minimum slope stability universally occurs when the seepage direction is given by λ = 90° − Φ. It also shows that for some upward seepage conditions, slope stability is limited by static liquefaction rather than by Coulomb failure. Close association between these liquefaction conditions and certain Coulomb failure conditions indicates that slope failure in such instances could be responsible for nearly spontaneous mobilization of destructive flowing soil masses on hillslopes.

Water Resources Research↗