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,297 records · Page 72Linked to original sources
USGS1996

The use of streambed temperature profiles to estimate the depth, duration, and rate of percolation beneath arroyos

Temporal variations in a streambed temperature profile between 30 and 300 cm beneath Tijeras Arroyo, New Mexico, were analyzed at 30-min intervals for 1990 to estimate the depth, duration, and rate of percolation during streamflows. The depth of percolation was clearly documented by the rapid response of the streambed temperature profile to streamflows. Results indicate that the streambed possessed small thermal gradients with significant diurnal variations from late November to late May, indicating that ephemeral streamflows created continuous, advection-dominated heat transport to depths below 300 cm during this period. Timing and duration of percolation suggested by temporal variations in the temperature profile were verified by comparison with measured streamflow records for the study reach over 1990. Percolation rates were estimated using a technique based on the travel time of the daily maximum temperature into the streambed. Percolation rates were compared with streambed seepage rates determined from measurements of streamflow loss, stream surface area, and stream evaporative loss for the entire study reach. Travel time estimates of streambed percolation rates ranged from 9 to 40 cm/hr, while streamflow estimates of streambed seepage rates ranged from 6 to 26 cm/hr during the study period. Discrepancies between streambed percolation and seepage rates may be caused by differences in the areal extent of measurements for percolation versus seepages rates. In summary, the depth, timing, and duration of streamflow-induced percolation were well documented by temporal variations in a single streambed temperature profile, while rates of percolation based on the temperature profile were about double the seepage rates based on streamflow records for the entire study reach.

Water Resources Research↗
USGS1997

Effects of solution mining of salt on wetland hydrology as inferred from tree rings

Radial growth and concentrations of selected elements within rings were studied in white pine ( Pinus strobus ) trees from a wetland in central New York approximately 5 km north of a salt-solution mining field that operated from 1889 to 1988. Trees seemingly document three sequential episodes of mine-induced alterations of groundwater discharge irrigating the wetland during the 100-year period. The radial growth of trees established before the onset of mining declined abruptly in the early 1890s and remained suppressed until about 1960, as did growth of numerous other trees that became established after the onset of mining. Suppressed pre-1960 radial growth coincided with the interval that surface water was injected into the saltbeds, suggesting that losses of injected water to the bedrock and/or unconsolidated deposits increased groundwater flow into the wetland. An abrupt and sustained enhancement of radial growth beginning about 1960 indicates that the wetland became drier, and thus more conducive to tree growth, when injection practices were discontinued in the late 1950s despite the continued pumping of brine. Following the cessation of mining in the late 1980s, head pressures again increased in the upper valley, driving chloride-enriched flow northward along regional bedding-plane fractures and into the wetland. Large concentrations of chloride were detected within the most recently formed rings of some trees. As the result of chloride-enriched irrigation, the radial growth of some trees declined, and some trees died. Thus trees have preserved evidence of a century of hydrologic alterations, unobtainable by other means, where the effects of brine mining have not been documented previously.

Water Resources Research↗
USGS1997

Deep well injection of brine from Paradox Valley, Colorado: Potential major precipitation problems remediated by nanofiltration

Groundwater brine seepage into the Dolores River in Paradox Valley, Colorado, increases the dissolved solids load of the Colorado River annually by ∼2.0 × 10 8 kg. To abate this natural contamination, the Bureau of Reclamation plans to pump ∼3540 m 3 /d of brine from 12 shallow wells located along the Dolores River. The brine, with a salinity of 250,000 mg/L, will be piped to the deepest (4.9 km) disposal well in the world and injected mainly into the Mississippian Leadville Limestone. Geochemical modeling indicates, and water-rock experiments confirm, that a huge mass of anhydrite (∼1.0 × 10 4 kg/d) likely will precipitate from the injected brine at downhole conditions of 120°C and 500 bars. Anhydrite precipitation could increase by up to 3 times if the injected brine is allowed to mix with the highly incompatible formation water of the Leadville Limestone and if the Mg in this brine dolomitizes the calcite of the aquifer. Laboratory experiments demonstrate that nanofiltration membranes, which are selective to divalent anions, provide a new technology that remediates the precipitation problem by removing ∼98% of dissolved SO 4 from the hypersaline brine. The fluid pressure used (50 bars) is much lower than would be required for traditional reverse osmosis membranes because nanofiltration membranes have a low rejection efficiency (5–10%) for monovalent anions. Our results indicate that the proportion of treatable brine increases from ∼60% to >85% with the addition of trace concentrations of a precipitation inhibitor and by blending the raw brine with the effluent stream.

Water Resources Research↗
USGS1997

Regional streamflow regimes and hydroclimatology of the United States

The dominant regions of interannual streamflow variability in the United States are defined, and their seasonality and persistence characteristics identified, using an orthogonally rotated principal components analysis (RPCA) of a climatically sensitive network of 559 stream gages for the period 1941–1988. This classification of streamflow regimes is comprehensive and unique in that separate analyses of the streamflow record, for each month of the year, are carried out to detail the month-to-month changes in the dominant streamflow patterns. Streamflow variations, or anomalies, in the Upper Mississippi, South Atlantic/Gulf, Far West, Ohio Valley, Northeast, and Eastern/Mid- Atlantic regions, as well as a pattern of opposing streamflow anomalies in the West, are observed in all seasons of the year. Anomalies in the Southern Plains and New England regions are observed in autumn, winter, and spring; those in the Rocky Mountains and Middle Mississippi regions occur in late spring and summer.

Water Resources Research↗
USGS1997

Flow to a well of finite diameter in a homogeneous, anisotropic water table aquifer

A Laplace transform solution is presented for the problem of flow to a partially penetrating well of finite diameter in a slightly compressible water table aquifer. The solution, which allows for evaluation of both pumped well and observation piezometer data, accounts for effects of well bore storage and skin and allows for the noninstantaneous release of water from the unsaturated zone. For instantaneous release of water from the unsaturated zone the solution approaches the line source solution derived by Neuman as the diameter of the pumped well approaches zero. Delayed piezometer response, which is significant during times of rapidly changing hydraulic head, is included in the theoretical treatment and shown to be an important factor in accurate evaluation of specific storage. By means of a hypothetical field example it is demonstrated that evaluations of specific storage ( S s ) using classical line source solutions may yield values of S s that are overestimated by a factor of 100 or more, depending upon the location of the observation piezometers and whether effects of delayed piezometer response are included in the analysis. Theoretical responses obtained with the proposed model are used to suggest methods for evaluating specific storage.

Water Resources Research↗
USGS1997

Strontium 87/strontium 86 as a tracer of mineral weathering reactions and calcium sources in an alpine/subalpine watershed, Loch Vale, Colorado

Sr isotopic ratios of atmospheric deposition, surface and subsurface water, and geologic materials were measured in an alpine/subalpine watershed to characterize weathering reactions and identify sources of dissolved Ca in stream water. Previous studies have noted an excess of Ca in stream water above that expected from stoichiometric weathering of the dominant bedrock minerals. Mixing calculations based on 87 Sr/ 86 Sr indicate that on an annual basis, 26 ± 7% of Ca export in streams is atmospherically derived, 23 ± 1% is from weathering of plagioclase, and the remainder is from weathering of calcite present in trace amounts in the bedrock. A potential source of error when applying Sr isotopes in catchment studies is determination of the 87 Sr/ 86 Sr of Sr released by mineral weathering, which is complicated by the wide range of mineral isotopic compositions, particularly in older rocks, and the variable rates at which the minerals weather. In this study, base-flow stream chemistry was used to represent the 87 Sr/ 86 Sr of Sr derived from mineral weathering because it effectively integrates the potentially variable isotopic composition of Sr released by weathering in the alpine environment.

Colorado↗
USGS1997

Multiphase flow modeling of a crude-oil spill site with a bimodal permeability distribution

Fluid saturation, particle-size distribution, and porosity measurements were obtained from 269 core samples collected from six boreholes along a 90-m transect at a subregion of a crude-oil spill site, the north pool, near Bemidji, Minnesota. The oil saturation data, collected 11 years after the spill, showed an irregularly shaped oil body that appeared to be affected by sediment spatial variability. The particle-size distribution data were used to estimate the permeability ( k ) and retention curves for each sample. An additional 344 k estimates were obtained from samples previously collected at the north pool. The 613 k estimates were distributed bimodal lognormally with the two population distributions corresponding to the two predominant lithologies: a coarse glacial outwash deposit and fine-grained interbedded lenses. A two-step geostatistical approach was used to generate a conditioned realization of k representing the bimodal heterogeneity. A cross-sectional multiphase flow model was used to simulate the flow of oil and water in the presence of air along the north pool transect for an 11-year period. The inclusion of a representation of the bimodal aquifer heterogeneity was crucial for reproduction of general features of the observed oil body. If the bimodal heterogeneity was characterized, hysteresis did not have to be incorporated into the model because a hysteretic effect was produced by the sediment spatial variability. By revising the relative permeability functional relation, an improved reproduction of the observed oil saturation distribution was achieved. The inclusion of water table fluctuations in the model did not significantly affect the simulated oil saturation distribution.

Water Resources Research↗
USGS1997

Experimental design for estimating parameters of rate-limited mass transfer: Analysis of stream tracer studies

Tracer experiments are valuable tools for analyzing the transport characteristics of streams and their interactions with shallow groundwater. The focus of this work is the design of tracer studies in high-gradient stream systems subject to advection, dispersion, groundwater inflow, and exchange between the active channel and zones in surface or subsurface water where flow is stagnant or slow moving. We present a methodology for (1) evaluating and comparing alternative stream tracer experiment designs and (2) identifying those combinations of stream transport properties that pose limitations to parameter estimation and therefore a challenge to tracer test design. The methodology uses the concept of global parameter uncertainty analysis, which couples solute transport simulation with parameter uncertainty analysis in a Monte Carlo framework. Two general conclusions resulted from this work. First, the solute injection and sampling strategy has an important effect on the reliability of transport parameter estimates. We found that constant injection with sampling through concentration rise, plateau, and fall provided considerably more reliable parameter estimates than a pulse injection across the spectrum of transport scenarios likely encountered in high-gradient streams. Second, for a given tracer test design, the uncertainties in mass transfer and storage-zone parameter estimates are strongly dependent on the experimental Damkohler number, DaI , which is a dimensionless combination of the rates of exchange between the stream and storage zones, the stream-water velocity, and the stream reach length of the experiment. Parameter uncertainties are lowest at DaI values on the order of 1.0. When DaI values are much less than 1.0 (owing to high velocity, long exchange timescale, and/or short reach length), parameter uncertainties are high because only a small amount of tracer interacts with storage zones in the reach. For the opposite conditions ( DaI ≫ 1.0), solute exchange rates are fast relative to stream-water velocity and all solute is exchanged with the storage zone over the experimental reach. As DaI increases, tracer dispersion caused by hyporheic exchange eventually reaches an equilibrium condition and storage-zone exchange parameters become essentially nonidentifiable.

Water Resources Research↗
USGS1997

Geochemical mole-balance modeling with uncertain data

Geochemical mole-balance models are sets of chemical reactions that quantitatively account for changes in the chemical and isotopic composition of water along a flow path. A revised mole-balance formulation that includes an uncertainty term for each chemical and isotopic datum is derived. The revised formulation is comprised of mole-balance equations for each element or element redox state, alkalinity, electrons, solvent water, and each isotope; a charge-balance equation and an equation that relates the uncertainty terms for p H, alkalinity, and total dissolved inorganic carbon for each aqueous solution; inequality constraints on the size of the uncertainty terms; and inequality constraints on the sign of the mole transfer of reactants. The equations and inequality constraints are solved by a modification of the simplex algorithm combined with an exhaustive search for unique combinations of aqueous solutions and reactants for which the equations and inequality constraints can be solved and the uncertainty terms minimized. Additional algorithms find only the simplest mole-balance models and determine the ranges of mixing fractions for each solution and mole transfers for each reactant that are consistent with specified limits on the uncertainty terms. The revised formulation produces simpler and more robust mole-balance models and allows the significance of mixing fractions and mole transfers to be evaluated. In an example from the central Oklahoma aquifer, inclusion of up to 5% uncertainty in the chemical data can reduce the number of reactants in mole-balance models from seven or more to as few as three, these being cation exchange, dolomite dissolution, and silica precipitation. In another example from the Madison aquifer, inclusion of the charge-balance constraint requires significant increases in the mole transfers of calcite, dolomite, and organic matter, which reduce the estimated maximum carbon 14 age of the sample by about 10,000 years, from 22,700 years to 12,600 years.

Water Resources Research↗
USGS1997

Soil-water movement under natural-site and waste-site conditions: A multiple-year field study in the Mojave Desert, Nevada

Soil-water movement under natural-site and simulated waste-site conditions were compared by monitoring four experimental sites in the Mojave Desert, Nevada, during a 5-year period: one vegetated soil profile, one soil profile where vegetation was removed, and two nonvegetated test trenches. Precipitation ranged from 14 to 162 mm/yr. Temporal changes in water content measured by neutron probe were limited to the upper 0.5–1 m; values ranged from 0.01 to 0.19 m 3 /m 3 . Water potential and temperature were measured by thermocouple psychrometers; 77% remained operable for ≥4.5 years. For vegetated soil, precipitation that accumulated in the upper 0.75 m of soil was removed by evapotranspiration: water potentials decreased seasonally by 4 to >8 MPa. During 2 years with below-average precipitation, water potentials below the app arent root zone decreased by 2.3 (1.2-m depth) to 0.4 MPa (5-m depth), and the gradients became predominantly upward. Water potentials then rebounded during 2 years with near- and above-average precipitation, and seasonally variant water potential gradients were reestablished above the 4.2-m depth. Under nonvegetated waste-site conditions, data indicated the long-term accumulation and shallow, but continued, penetration of precipitation: water potentials showed moisture penetration to depths of 0.75−1.85 m. The method of simulated-waste drum placement (stacked versus random) and the associated differences in subsidence showed no measurable influence on the water balance of the trenches: subsidence totaled ≤13 mm during the study. Water potentials below the trenches and below the 2-m depth for the nonvegetated soil remained low (≈−5.5 to −7.5 MPa) and indicated the persistence of typically upward driving forces for isothermal water flow. Water fluxes estimated from water potential and temperature data suggested that isothermal liquid, isothermal vapor, and nonisothermal vapor flow need to be considered in the conceptualization of unsaturated flow at the field sites. Below the depth of temporal water content change, the estimated liquid fluxes ranged from 10 −10 to 10 −15 cm/s, isothermal vapor fluxes ranged from 10 −10 to 10 −13 cm/s, and the nonisothermal vapor fluxes ranged from 10 −8 to 10 −10 cm/s.

Nevada↗
USGS1997

Methods for predicting peak discharge of floods caused by failure of natural and constructed earthen dams

Floods from failures of natural and constructed dams constitute a widespread hazard to people and property. Expeditious means of assessing flood hazards are necessary, particularly in the case of natural dams, which may form suddenly and unexpectedly. We revise statistical relations (derived from data for past constructed and natural dam failures) between peak discharge ( Q p ) and water volume released ( V 0 ) or drop in lake level ( d ) but assert that such relations, even when cast into a dimensionless form, are of limited utility because they fail to portray the effect of breach-formation rate. We then analyze a simple, physically based model of dam-breach formation to show that the hydrograph at the breach depends primarily on a dimensionless parameter η= kV 0 / g l/2 d 7/2 , where k is the mean erosion rate of the breach and g is acceleration due to gravity. The functional relationship between Q p and η takes asymptotically distinct forms depending on whether η ≪ 1 (relatively slow breach formation or small lake volume) or η ≫ 1 (relatively fast breach formation or large lake volume). Theoretical predictions agree well with data from dam failures for which k , and thus η, can be estimated. The theory thus provides a rapid means of predicting the plausible range of values of peak discharge at the breach in an earthen dam as long as the impounded water volume and the water depth at the dam face can be estimated.

Water Resources Research↗
USGS1997

An algorithm for computing moments-based flood quantile estimates when historical flood information is available

This paper presents the expected moments algorithm (EMA), a simple and efficient method for incorporating historical and paleoflood information into flood frequency studies. EMA can utilize three types of at-site flood information: systematic stream gage record; information about the magnitude of historical floods; and knowledge of the number of years in the historical period when no large flood occurred. EMA employs an iterative procedure to compute method-of-moments parameter estimates. Initial parameter estimates are calculated from systematic stream gage data. These moments are then updated by including the measured historical peaks and the expected moments, given the previously estimated parameters, of the below-threshold floods from the historical period. The updated moments result in new parameter estimates, and the last two steps are repeated until the algorithm converges. Monte Carlo simulations compare EMA, Bulletin 17B 's [ United States Water Resources Council , 1982] historically weighted moments adjustment, and maximum likelihood estimators when fitting the three parameters of the log-Pearson type III distribution. These simulations demonstrate that EMA is more efficient than the Bulletin 17B method, and that it is nearly as efficient as maximum likelihood estimation (MLE). The experiments also suggest that EMA has two advantages over MLE when dealing with the log-Pearson type III distribution: It appears that EMA estimates always exist and that they are unique, although neither result has been proven. EMA can be used with binomial or interval-censored data and with any distributional family amenable to method-of-moments estimation.

Water Resources Research↗
USGS1997

Evaluation of process errors in bed load sampling using a Dune Model

Reliable estimates of the streamwide bed load discharge obtained using sampling devices are dependent upon good at-a-point knowledge across the full width of the channel. Using field data and information derived from a model that describes the geometric features of a dune train in terms of a spatial process observed at a fixed point in time, we show that sampling errors decrease as the number of samples collected increases, and the number of traverses of the channel over which the samples are collected increases. It also is preferable that bed load sampling be conducted at a pace which allows a number of bed forms to pass through the sampling cross section. The situations we analyze and simulate pertain to moderate transport conditions in small rivers. In such circumstances, bed load sampling schemes typically should involve four or five traverses of a river, and the collection of 20–40 samples at a rate of five or six samples per hour. By ensuring that spatial and temporal variability in the transport process is accounted for, such a sampling design reduces both random and systematic errors and hence minimizes the total error involved in the sampling process.

Water Resources Research↗
USGS1997

Constant-concentration boundary condition: Lessons from the HYDROCOIN variable-density groundwater benchmark problem

In a solute-transport model, if a constant-concentration boundary condition is applied at a node in an active flow field, a solute flux can occur by both advective and dispersive processes. The potential for advective release is demonstrated by reexamining the Hydrologic Code Intercomparison (HYDROCOIN) project case 5 problem, which represents a salt dome overlain by a shallow groundwater system. The resulting flow field includes significant salinity and fluid density variations. Several independent teams simulated this problem using finite difference or finite element numerical models. We applied a method-of-characteristics model (MOCDENSE). The previous numerical implementations by HYDROCOIN teams of a constant-concentration boundary to represent salt release by lateral dispersion only (as stipulated in the original problem definition) was flawed because this boundary condition allows the release of salt into the flow field by both dispersion and advection. When the constant-concentration boundary is modified to allow salt release by dispersion only, significantly less salt is released into the flow field. The calculated brine distribution for case 5 depends very little on which numerical model is used, as long as the selected model is solving the proper equations. Instead, the accuracy of the solution depends strongly on the proper conceptualization of the problem, including the detailed design of the constant-concentration boundary condition. The importance and sensitivity to the manner of specification of this boundary does not appear to have been recognized previously in the analysis of this problem.

Water Resources Research↗
USGS1997

Regional interpretation of water-quality monitoring data

We describe a method for using spatially referenced regressions of contaminant transport on watershed attributes (SPARROW) in regional water-quality assessment. The method is designed to reduce the problems of data interpretation caused by sparse sampling, network bias, and basin heterogeneity. The regression equation relates measured transport rates in streams to spatially referenced descriptors of pollution sources and land-surface and stream-channel characteristics. Regression models of total phosphorus (TP) and total nitrogen (TN) transport are constructed for a region defined as the nontidal conterminous United States. Observed TN and TP transport rates are derived from water-quality records for 414 stations in the National Stream Quality Accounting Network. Nutrient sources identified in the equations include point sources, applied fertilizer, livestock waste, nonagricultural land, and atmospheric deposition (TN only). Surface characteristics found to be significant predictors of land-water delivery include soil permeability, stream density, and temperature (TN only). Estimated instream decay coefficients for the two contaminants decrease monotonically with increasing stream size. TP transport is found to be significantly reduced by reservoir retention. Spatial referencing of basin attributes in relation to the stream channel network greatly increases their statistical significance and model accuracy. The method is used to estimate the proportion of watersheds in the conterminous United States (i.e., hydrologic cataloging units) with outflow TP concentrations less than the criterion of 0.1 mg/L, and to classify cataloging units according to local TN yield (kg/km 2 /yr).

Water Resources Research↗
USGS1997

Hydrological and chemical estimates of the water balance of a closed-basin lake in north central Minnesota

Chemical mass balances for sodium, magnesium, chloride, dissolved organic carbon, and oxygen 18 were used to estimate groundwater seepage to and from Williams Lake, Minnesota, over a 15-month period, from April 1991 through June 1992. Groundwater seepage to the lake and seepage from the lake to groundwater were determined independently using a flow net approach using data from water table wells installed as part of the study. Hydrogeological analysis indicated groundwater seepage to the lake accounted for 74% of annual water input to the lake; the remainder came from atmospheric precipitation, as determined from a gage in the watershed and from nearby National Weather Service gages. Seepage from the lake accounted for 69% of annual water losses from the lake; the remainder was removed by evaporation, as determined by the energy budget method. Calculated annual water loss exceeded calculated annual water gain, and this imbalance was double the value of the independently measured decrease in lake volume. Seepage to the lake determined from oxygen 18 was larger (79% of annual water input) than that determined from the flow net approach and made the difference between calculated annual water gain and loss consistent with the independently measured decrease in lake volume. Although the net difference between volume of seepage to the lake and volume of seepage from the lake was 1% of average lake volume, movement of water into and out of the lake by seepage represented an annual exchange of groundwater with the lake equal to 26–27% of lake volume. Estimates of seepage to the lake from sodium, magnesium, chloride, and dissolved organic carbon did not agree with the values determined from flow net approach or oxygen 18. These results indicated the importance of using a combination of hydrogeological and chemical approaches to define volume of seepage to and from Williams Lake and identify uncertainties in chemical fluxes.

Minnesota↗
USGS1998

A nonlinear model for analysis of slug-test data

While doing slug tests in high-permeability aquifers, we have consistently seen deviations from the expected response of linear theoretical models. Normalized curves do not coincide for various initial heads, as would be predicted by linear theories, and are shifted to larger times for higher initial heads. We have developed a general nonlinear model based on the Navier-Stokes equation, nonlinear frictional loss, non-Darcian flow, acceleration effects, radius changes in the well bore, and a Hvorslev model for the aquifer, which explains these data features. The model produces a very good fit for both oscillatory and nonoscillatory field data, using a single set of physical parameters to predict the field data for various initial displacements at a given well. This is in contrast to linear models which have a systematic lack of fit and indicate that hydraulic conductivity varies with the initial displacement. We recommend multiple slug tests with a considerable variation in initial head displacement to evaluate the possible presence of nonlinear effects. Our conclusion is that the nonlinear model presented here is an excellent tool to analyze slug tests, covering the range from the underdamped region to the overdamped region.

Water Resources Research↗
USGS1998

High-resolution Monte Carlo simulation of flow and conservative transport in heterogeneous porous media: 2. Transport results

In this, the second of two papers concerned with the use of numerical simulation to examine flow and transport parameters in heterogeneous porous media via Monte Carlo methods, results from the transport aspect of these simulations are reported on. Transport simulations contained herein assume a finite pulse input of conservative tracer, and the numerical technique endeavors to realistically simulate tracer spreading as the cloud moves through a heterogeneous medium. Medium heterogeneity is limited to the hydraulic conductivity field, and generation of this field assumes that the hydraulic-conductivity process is second-order stationary. Methods of estimating cloud moments, and the interpretation of these moments, are discussed. Techniques for estimation of large-time macrodispersivities from cloud second-moment data, and for the approximation of the standard errors associated with these macrodispersivities, are also presented. These moment and macrodispersivity estimation techniques were applied to tracer clouds resulting from transport scenarios generated by specific Monte Carlo simulations. Where feasible, moments and macrodispersivities resulting from the Monte Carlo simulations are compared with first- and second-order perturbation analyses. Some limited results concerning the possible ergodic nature of these simulations, and the presence of non-Gaussian behavior of the mean cloud, are reported on as well.

Water Resources Research↗