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,027 records · Page 57Linked to original sources
USGS2010

Use of output from high‐resolution atmospheric models in landscape‐scale hydrologic models: An assessment

In this paper we investigate the feasibility of coupling regional climate models (RCMs) with landscape‐scale hydrologic models (LSHMs) for studies of the effects of climate on hydrologic systems. The RCM used is the National Center for Atmospheric Research/Pennsylvania State University mesoscale model (MM4). Output from two year‐round simulations (1983 and 1988) over the western United States is used to drive a lake model for Pyramid Lake in Nevada and a streamfiow model for Steamboat Creek in Oregon. Comparisons with observed data indicate that MM4 is able to produce meteorologic data sets that can be used to drive hydrologic models. Results from the lake model simulations indicate that the use of MM4 output produces reasonably good predictions of surface temperature and evaporation. Results from the streamflow simulations indicate that the use of MM4 output results in good simulations of the seasonal cycle of streamflow, but deficiencies in simulated wintertime precipitation resulted in underestimates of streamflow and soil moisture. Further work with climate (multiyear) simulations is necessary to achieve a complete analysis, but the results from this study indicate that coupling of LSHMs and RCMs may be a useful approach for evaluating the effects of climate change on hydrologic systems.

Water Resources Research↗
USGS2010

Groundwater as a nonpoint source of atrazine and deethylatrazine in a river during base flow conditions

Alluvial groundwater adjacent to the main stem river is the principal nonpoint source of atrazine and deethylatrazine in the Cedar River of Iowa after the river has been in base flow conditions for 5 days. Between two sites along a 116-km reach of the Cedar River, tributaries contributed about 25% of the increase in the atrazine and deethylatrazine load, whereas groundwater from the alluvial aquifer contributed at least 75% of the increase in load. Within the study area, tributaries aggregate almost all of the discharge from tile drains, and yet the tributaries still only contribute 25% of the increase in loads in the main stem river. At an unfamned study site adjacent to the Cedar River, the sources of atrazine and deethylatrazine in the alluvial groundwater are bank storage of river water and groundwater recharge from areas distant from the river. Atrazine and deethylatrazine associated with bank storage water will provide larger concentrations to the river during early base flow conditions. After the depletion of bank storage, stable and smaller concentrations of atrazine and deethylatrazine, originating from groundwater recharge, continue to be discharged from the alluvial aquifer to the river; thus these results indicate that alluvial aquifers are an important nonpoint source of atrazine and deethylatrazine in rivers during base flow.

Iowa↗
USGS2010

Comparison of Penman-Monteith, Shuttleworth-Wallace, and modified Priestley-Taylor evapotranspiration models for wildland vegetation in semiarid rangeland

Eddy correlation measurements of sensible and latent heat flux are used with measurements of net radiation, soil heat flux, and other micrometeorological variables to develop the Penman-Monteith, Shuttleworth-Wallace, and modified Priestley-Taylor evapotranspiration models for use in a sparsely vegetated, semiarid rangeland. The Penman-Monteith model, a one-component model designed for use with dense crops, is not sufficiently accurate ( r 2 = 0.56 for hourly data and r 2 = 0.60 for daily data). The Shuttleworth-Wallace model, a two-component logical extension of the Penman-Monteith model for use with sparse crops, performs significantly better ( r 2 = 0.78 for hourly data and r 2 = 0.85 for daily data). The modified Priestley-Taylor model, a one-component simplified form of the Penman potential evapotranspiration model, surprisingly performs as well as the Shuttle worth-Wallace model. The rigorous Shuttleworth-Wallace model predicts that about one quarter of the vapor flux to the atmosphere is from bare-soil evaporation. Further, during daylight hours, the small leaves are sinks for sensible heat produced at the hot soil surface.

Water Resources Research↗
USGS2010

Simulation of fluid distributions observed at a crude oil spill site incorporating hysteresis, oil entrapment, and spatial variability of hydraulic properties

Subsurface oil, water, and air saturation distributions were determined using 146 samples collected from seven boreholes along a 120-m transect at a crude oil spill site near Bemidji, Minnesota. The field data, collected 10 years after the spill, show a clearly defined oil body that has an oil saturation distribution that appears to be influenced by sediment heterogeneities and water table fluctuations. The center of the oil body has depressed the water-saturated zone boundary and the oil appears to have migrated laterally within the capillary fringe. A multiphase cross-sectional flow model was developed and used to simulate the movement of oil and water at the spill site. Comparisons between observed and simulated oil saturation distributions serve as an indicator of the appropriateness of using such models to predict the actual spread of organic immiscible liquids at spill sites. Sediment hydraulic properties used in the model were estimated from particle size data. The general large-scale features of the observed oil body were reproduced only when hysteresis with oil entrapment and representations of observed spatial variability of hydraulic properties were incorporated into the model. The small-scale details of the observed subsurface oil distribution were not reproduced in the simulations. The discrepancy between observed and simulated oil distributions reflects the considerable uncertainty in model parameter estimates and boundary conditions, three-phase capillary pressure-saturation-relative permeability functions, representations of spatial variability of hydraulic properties, and hydrodynamics of the groundwater flow system at the study site.

Water Resources Research↗
USGS2010

Transport and accumulation of radionuclides and stable elements in a Missouri River Reservoir

Several long sediment cores from the Cheyenne River Embayment of Lake Oahe, a 250-km-long Missouri River reservoir in South Dakota, have been analyzed for radionuclides and stable elements. The combination of fine-scale sampling and rapid sedimentation produces radionuclide distributions that can be used to estimate the detailed chronology of particle transport processes in the Oahe reservoir system. A self-consistent and quantitative treatment of the 137 Cs data suggests processes to which characteristic times may be associated. Times that characterize system-wide processes include (1) an integration time of several years reflecting retention of the sediment-bound tracer in regions within or external to the reservoir, (2) a relaxation time of approximately 15 years reflecting a decreasing rate of sediment accumulation ascribed to shoreline stabilization, (3) a time of a few months characterizing the breadth of riverine signatures in cores due to integration effects in the Cheyenne River system and deltaic deposits, and (4) times of a few years associated with propagation of riverine load signatures along the embayment. The distribution of total sedimentary arsenic confirms the validity of the variable sedimentation model. In 1977, a tailings retention facility was built at the Homestake Mine site, and the unrestricted input of As ceased. As a result of this remedial action, the concentration of sedimentary As decreased dramatically. In the upper section of the core, above the depth represented by the year 1976, the concentration of As decreases tenfold. In this same core the distribution of lithologically discriminating chemical elements, calcium and vanadium, relate to major flow events in the Cheyenne River basin. Because there is minimal diagenesis of chemical constituents in these rapidly accumulating sediments, stable element signatures, in addition to radiotracers, may be used to reconstruct hydrologic events in drainage basins that contribute sediment to lakes and reservoirs.

South Dakota↗
USGS2010

A finite-volume Eulerian-Lagrangian Localized Adjoint Method for solution of the advection-dispersion equation

A new mass-conservative method for solution of the one-dimensional advection-dispersion equation is derived and discussed. Test results demonstrate that the finite-volume Eulerian-Lagrangian localized adjoint method (FVELLAM) outperforms standard finite-difference methods, in terms of accuracy and efficiency, for solute transport problems that are dominated by advection. For dispersion-dominated problems, the performance of the method is similar to that of standard methods. Like previous ELLAM formulations, FVELLAM systematically conserves mass globally with all types of boundary conditions. FVELLAM differs from other ELLAM approaches in that integrated finite differences, instead of finite elements, are used to approximate the governing equation. This approach, in conjunction with a forward tracking scheme, greatly facilitates mass conservation. The mass storage integral is numerically evaluated at the current time level, and quadrature points are then tracked forward in time to the next level. Forward tracking permits straightforward treatment of inflow boundaries, thus avoiding the inherent problem in backtracking, as used by most characteristic methods, of characteristic lines intersecting inflow boundaries. FVELLAM extends previous ELLAM results by obtaining mass conservation locally on Lagrangian space-time elements. Details of the integration, tracking, and boundary algorithms are presented. Test results are given for problems in Cartesian and radial coordinates.

Water Resources Research↗
USGS2010

Measurement of variation in soil solute tracer concentration across a range of effective pore sizes

Solute transport concepts in soil are based on speculation that solutes are distributed nonuniformly within large and small pores. Solute concentrations have not previously been measured across a range of pore sizes and examined in relation to soil hydrological properties. For this study, modified pressure cells were used to measure variation in concentration of a solute tracer across a range of pore sizes. Intact cores were removed from the site of a field tracer experiment, and soil water was eluted from 10 or more discrete classes of pore size. Simultaneous changes in water content and unsaturated hydraulic conductivity were determined on cores using standard pressure cell techniques. Bromide tracer concentration varied by as much as 100% across the range of pore sizes sampled. Immediately following application of the bromide tracer on field plots, bromide was most concentrated in the largest pores; concentrations were lower in pores of progressively smaller sizes. After 27 days, bromide was most dilute in the largest pores and concentrations were higher in the smaller pores. A sharp, threefold decrease in specific water capacity during elution indicated separation of two major pore size classes at a pressure of 47 cm H 2 O and a corresponding effective pore diameter of 70 μm. Variation in tracer concentration, on the other hand, was spread across the entire range of pore sizes investigated in this study. A two-porosity characterization of the transport domain, based on water retention criteria, only broadly characterized the pattern of variation in tracer concentration across pore size classes during transport through a macroporous soil.

Water Resources Research↗
USGS2010

Evaluation of the energy budget method of determining evaporation at Williams Lake, Minnesota, using alternative instrumentation and study approaches

Best estimates of evaporation at Williams Lake, north central Minnesota, were determined by the energy budget method using optimum sensors and optimum placement of sensors. These best estimates are compared with estimates derived from using substitute data to determine the effect of using less accurate sensors, simpler methods, or remotely measured data. Calculations were made for approximately biweekly periods during five open water seasons. For most of the data substitutions that affected the Bowen ratio, new values of evaporation differed little from best estimates. The three data substitution methods that caused the largest deviations from the best evaporation estimates were (1) using changes in the daily average surface water temperature as an indicator of the lake heat storage term, (2) using shortwave radiation, air temperature, and atmospheric vapor pressure data from a site 110 km away, and (3) using an analog surface water temperature probe. Recalculations based on these data substitutions resulted in differences from the best estimates as much as 89%, 21%, and 10%, respectively. The data substitution method that provided evaporation values that most closely matched the best estimates was measurement of the lake heat storage term at one location in the lake, rather than at 16 locations. Evaporation values resulting from this substitution method usually were within 2% of the best estimates.

Minnesota↗
USGS2010

Estimating discharge of shallow groundwater by transpiration from greasewood in the Northern Great Basin

Evapotranspiration from bare soil and phreatophytes is a principal mechanism of groundwater discharge in arid and semiarid regions of the midwestern and western United States including the Great Basin. The imbalance between independent estimates of groundwater recharge from precipitation and of groundwater discharge based on estimates of groundwater evapotranspiration leads to large uncertainties in groundwater budgets. Few studies have addressed this problem. Energy budget micrometeorological field studies were conducted in a stand of sparse-canopy greasewood growing in an area of shallow groundwater in the western Great Basin during the summer of 1989. The data were used to calculate above-canopy fluxes of sensible and latent heat using the energy budget-Bowen ratio method. The calculated energy budget fluxes were used, with soil surface and plant canopy temperature measurements, to calibrate and apply a two-component, energy-combination model that partitions the energy and heat fluxes between bare soil and the canopy. This permitted the separation of evaporation from the soil and transpiration from greasewood. The calibrated model was used to estimate daily transpiration of groundwater by greasewood growing in an area with a depth to water of about 2 m. The daily rate of groundwater discharge by transpiration during July and August was estimated to be 2.4 mm. A period of 100 days for groundwater discharge at this rate was assumed to estimate an annual discharge of groundwater of 24 cm at the study site.

Great Basin↗
USGS2010

Role of physical heterogeneity in the interpretation of small-scale laboratory and field observations of bacteria, microbial-sized microsphere, and bromide transport through aquifer sediments

The effect of physical variability upon the relative transport behavior of microbial-sized microspheres, indigenous bacteria, and bromide was examined in field and flow-through column studies for a layered, but relatively well sorted, sandy glaciofluvial aquifer. These investigations involved repacked, sieved, and undisturbed aquifer sediments. In the field, peak abundance of labeled bacteria traveling laterally with groundwater flow 6 m downgradient from point of injection was coincident with the retarded peak of carboxylated microspheres (retardation factor, RF = 1.7) at the 8.8 m depth, but preceded the bromide peak and the retarded microsphere peak (RF = 1.5) at the 9.0 m depth. At the 9.5 m depth, the bacterial peak was coincident with both the bromide and the microsphere peaks. Although sorption appeared to be a predominant mechanism responsible for immobilization of microbial-sized microspheres in the aquifer, straining appeared to be primarily responsible for their removal in 0.6-m-long columns of repacked, unsieved aquifer sediments. The manner in which the columns were packed also affected optimal size for microsphere transport, which in one experiment was near the size of the small (∼2 μm) groundwater protozoa (flagellates). These data suggest that variability in aquifer sediment structure can be important in interpretation of both small-scale field and laboratory experiments examining microbial transport behavior.

Water Resources Research↗
USGS2010

A Fast Fourier transform stochastic analysis of the contaminant transport problem

A three-dimensional stochastic analysis of the contaminant transport problem is developed in the spirit of Naff (1990). The new derivation is more general and simpler than previous analysis. The fast Fourier transformation is used extensively to obtain numerical estimates of the mean concentration and various spatial moments. Data from both the Borden and Cape Cod experiments are used to test the methodology. Results are comparable to results obtained by other methods, and to the experiments themselves.

Water Resources Research↗
USGS2010

Simulating the volatilization of solvents in unsaturated soils during laboratory and field infiltration experiments

This paper describes laboratory and field experiments which were conducted to study the dynamics of trichloroethylene (TCE) as it volatilized from contaminated groundwater and diffused in the presence of infiltrating water through the unsaturated soil zone to the land surface. The field experiments were conducted at the Picatinny Arsenal, which is part of the United States Geological Survey Toxic Substances Hydrology Program. In both laboratory and field settings the gas and water phase concentrations of TCE were not in equilibrium during infiltration. Gas-water mass transfer rate constants were calibrated to the experimental data using a model in which the water phase was treated as two phases: a mobile water phase and an immobile water phase. The mass transfer limitations of a volatile organic compound between the gas and liquid phases were described explicitly in the model. In the laboratory experiment the porous medium was nonsorbing, and water infiltration rates ranged from 0.076 to 0.28 cm h −1 . In the field experiment the water infiltration rate was 0.34 cm h −1 , and sorption onto the soil matrix was significant. The laboratory-calibrated gas-water mass transfer rate constant is 3.3×10 −4 h −1 for an infiltration rate of 0.076 cm h −1 and 1.4×10 −3 h −1 for an infiltration rate of 0.28 cm h −1 . The overall mass transfer rate coefficients, incorporating the contribution of mass transfer between mobile and immobile water phases and the variation of interfacial area with moisture content, range from 3×10 −4 h −1 to 1×10 −2 h −1 . A power law model relates the gas-water mass transfer rate constant to the infiltration rate and the fraction of the water phase which is mobile. It was found that the results from the laboratory experiments could not be extrapolated to the field. In order to simulate the field experiment the very slow desorption of TCE from the soil matrix was incorporated into the mathematical model. When desorption from the soil matrix was added to the model, the calibrated gas-water mass transfer rate constant is 2 orders of magnitude lower than that predicted using the power law model developed for the nonsorbing laboratory soil.

Water Resources Research↗
USGS2010

Correction of stream quality trends for the effects of laboratory measurement bias

We present a statistical model relating measurements of water quality to associated errors in laboratory methods. Estimation of the model allows us to correct trends in water quality for long-term and short-term variations in laboratory measurement errors. An illustration of the bias correction method for a large national set of stream water quality and quality assurance data shows that reductions in the bias of estimates of water quality trend slopes are achieved at the expense of increases in the variance of these estimates. Slight improvements occur in the precision of estimates of trend in bias by using correlative information on bias and water quality to estimate random variations in measurement bias. The results of this investigation stress the need for reliable, long-term quality assurance data and efficient statistical methods to assess the effects of measurement errors on the detection of water quality trends.

Water Resources Research↗
USGS2010

Mean flow and turbulence fields over two-dimensional bed forms

Detailed laser-Doppler velocity and Reynolds stress measurements over fixed two-dimensional bed forms are used to investigate the coupling between the mean flow and turbulence and to examine effects that play a role in producing the bed form instability and finite amplitude stability. The coupling between the mean flow and the turbulence is explored in both a spatially averaged sense, by determining the structure of spatially averaged velocity and Reynolds stress profiles, and a local sense, through computation of eddy viscosities and length scales. The measurements show that there is significant interaction between the internal boundary layer and the overlying wake turbulence produced by separation at the bed form crest. The interaction produces relatively low correlation coefficients in the internal boundary layer, which suggests that using local bottom stress to predict bed load flux may not only be erroneous, it may also disregard the essence of the bed form instability mechanism. The measurements also indicate that topographically induced acceleration over the bed form stoss slope has a more significant effect in damping the turbulence over bed forms than was previously supposed, which is hypothesized to play a role in the stabilization of fully developed bed forms.

Water Resources Research↗
USGS2010

An analytic solution of the stochastic storage problem applicable to soil water

The accumulation of soil water during rainfall events and the subsequent depletion of soil water by evaporation between storms can be described, to first order, by simple accounting models. When the alternating supplies (precipitation) and demands (potential evaporation) are viewed as random variables, it follows that soil-water storage, evaporation, and runoff are also random variables. If the forcing (supply and demand) processes are stationary for a sufficiently long period of time, an asymptotic regime should eventually be reached where the probability distribution functions of storage, evaporation, and runoff are stationary and uniquely determined by the distribution functions of the forcing. Under the assumptions that the potential evaporation rate is constant, storm arrivals are Poisson-distributed, rainfall is instantaneous, and storm depth follows an exponential distribution, it is possible to derive the asymptotic distributions of storage, evaporation, and runoff analytically for a simple balance model. A particular result is that the fraction of rainfall converted to runoff is given by (1 - R −1 )/( e α(1− R −1) − R −1 ), in which R is the ratio of mean potential evaporation to mean rainfall and a is the ratio of soil water-holding capacity to mean storm depth. The problem considered here is analogous to the well-known problem of storage in a reservoir behind a dam, for which the present work offers a new solution for reservoirs of finite capacity. A simple application of the results of this analysis suggests that random, intraseasonal fluctuations of precipitation cannot by themselves explain the observed dependence of the annual water balance on annual totals of precipitation and potential evaporation.

Water Resources Research↗
USGS2010

Chlorofluorocarbons (CCl3F and CCl2F2) as dating tools and hydrologic tracers in shallow groundwater of the Delmarva Peninsula, Atlantic Coastal Plain, United States

Concentrations of the Chlorofluorocarbons (CFCs) CFC-11 and CFC-12 were determined in groundwater from coastal plain sediments of the Delmarva Peninsula. CFC-modeled ages were calculated independently for CFC-11 and CFC-12, and agreed to within 2–3 years in the majority of the waters. Recharge temperatures, determined from dissolved nitrogen and argon concentrations, varied from 9±2°C over most of the peninsula to 14±2°C at the southernmost tip of the peninsula in Virginia. The CFC-modeled ages were examined in relation to the known hydrogeologic environment, both on regional scales and in more intensively sampled local scale networks. The CFC-modeled recharge years and measured tritium concentrations were used to reconstruct a tritium input function that was compared to the modeled tritium plus 3 He distribution. Most of the present distribution of tritium in Delmarva groundwater is consistent with low dispersivities. The results of the study strongly support the use of CFCs for dating shallow, aerobic groundwater.

Delaware, Maryland↗
USGS2010

Numerical investigation of steady liquid water flow in a variably saturated fracture network

Numerical simulation was used to study steady liquid water movement in a 5-m by 5-m vertical section containing a hypothetical fracture network under conditions of variable saturation. The fracture network was assumed to be embedded within an impermeable rock matrix. Three variations of a network were considered. The “mixed” network consisted of two fracture sets, a subvertical set containing five 125 μm average aperture fractures and a subhorizontal set containing four 25 μm average aperture fractures. The other two networks had identical fracture orientation and contained either all 125 μm or all 25 μm average aperture fractures. The TOUGH simulator was used to calculate the total steady liquid water flux through the network, the flux through individual fracture segments, and the pressure head at each fracture segment. A unit hydraulic gradient was imposed on the network by applying fixed pressure head boundaries (ranging from −0.25 to 0.0 m of water) of equal value to the top and bottom. Saturation and permeability versus pressure head relations for the two sets of fractures were determined with the VSFRAC model, which assumed that aperture was variable within an individual fracture. Results showed that the spatial distributions of pressure head and flux within the network, as well as the location of the dominant pathways, depended strongly on the prescribed boundary pressure head. For the mixed network, both pressure head and flux tended to become more spatially uniform when the boundary pressure head approached the pressure head at which the permeability thickness products of the large- and small-aperture fractures are equal (the crossover pressure head). These results imply that for systems similar to the one considered here, interpretation of actual measurements of pressure head and flux may be quite complex, and that representation of variably saturated fracture networks as an equivalent continuum may be more valid for some ranges in pressure head than for others. Equivalent permeability as a function of pressure head was calculated for the fracture network, illustrating how information collected on individual fractures may be used to estimate the flow properties of rock at larger scales.

Water Resources Research↗
USGS2010

Adaptive estimation of the log fluctuating conductivity from tracer data at the Cape Cod Site

An adaptive estimation scheme is used to obtain the integral scale and variance of the log-fluctuating conductivity at the Cape Cod site based on the fast Fourier transform/stochastic model of Deng et al. (1993) and a Kalmanlike filter. The filter incorporates prior estimates of the unknown parameters with tracer moment data to adaptively obtain improved estimates as the tracer evolves. The results show that significant improvement in the prior estimates of the conductivity can lead to substantial improvement in the ability to predict plume movement. The structure of the covariance function of the log-fluctuating conductivity can be identified from the robustness of the estimation. Both the longitudinal and transverse spatial moment data are important to the estimation.

Massachusetts↗