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

Advection and diffusion in a variable-salinity confining layer

A numerical model that simulates groundwater flow and solute transport for cases in which fluid properties are variable was applied in one dimension (vertical) to the shallow, low-permeability, clayey, confining layer in Doñana National Park in southwestern Spain. The salinity in the 80-m-thick confining layer decreases from a brine near the land surface to fresh water near its base. Results of model simulations indicate that the system could be in or close to a steady state condition. The model calibration was very sensitive to small variations in individual model parameters and was nonunique in the sense that equally good calibrations could be achieved by compensatory joint perturbations in the permeability, diffusion coefficient, and overall governing hydraulic gradient. At present, there is probably an upward flow of the order of 1 mm yr −1 to 1 cm yr −1 and a balance in the solute flux between upward advection and downward diffusion. The time scale of calculated responses to changes in boundary conditions in this low-permeability system ranges from thousands to hundreds of thousands of years when considering extremes within the range of uncertainty of values of the evaluated parameters.

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↗
USGS2010

The use of simulation and multiple environmental tracers to quantify groundwater flow in a shallow aquifer

Measurements of the concentrations of chlorofluorocarbons (CFCs), tritium, and other environmental tracers can be used to calculate recharge ages of shallow groundwater and estimate rates of groundwater movement. Numerical simulation also provides quantitative estimates of flow rates, flow paths, and mixing properties of the groundwater system. The environmental tracer techniques and the hydraulic analyses each contribute to the understanding and quantification of the flow of shallow groundwater. However, when combined, the two methods provide feedback that improves the quantification of the flow system and provides insight into the processes that are the most uncertain. A case study near Locust Grove, Maryland, is used to investigate the utility of combining groundwater age dating, based on CFCs and tritium, and hydraulic analyses using numerical simulation techniques. The results of the feedback between an advective transport model and the estimates of groundwater ages determined by the CFCs improve a quantitative description of the system by refining the system conceptualization and estimating system parameters. The plausible system developed with this feedback between the advective flow model and the CFC ages is further tested using a solute transport simulation to reproduce the observed tritium distribution in the groundwater. The solute transport simulation corroborates the plausible system developed and also indicates that, for the system under investigation with the data obtained from 0.9-m-long (3-foot-long) well screens, the hydrodynamic dispersion is negligible. Together the two methods enable a coherent explanation of the flow paths and rates of movement while indicating weaknesses in the understanding of the system that will require future data collection and conceptual refinement of the groundwater system.

Water Resources Research↗
USGS2010

How permeable are clays and shales?

The permeability of argillaceous formations, although rarely measured and poorly understood, is commonly a critical parameter in analyses of subsurface flow. Data now available suggest a regular relation between permeability and porosity in clays and shales and permeabilities that, even at large scales, are significantly lower than usually assumed. Permeabilities between 10 −23 and 10 −17 m 2 have been obtained at porosities between 0.1 and 0.4 in both laboratory and regional studies. Although it is clear that transmissive fractures or other heterogeneities control the large-scale hydraulic behavior of certain argillaceous units, the permeability of many others is apparently scale independent. These results have significant implications for understanding fluid transport rates and abnormal pressure generation in basins, and could prove important for waste isolation efforts.

Water Resources Research↗
USGS2010

Hydrologic processes controlling sulfate mobility in a small forested watershed

Hydrologic controls on sulfate mobility were investigated in a forested catchment in the Georgia Piedmont using a watershed mass balance approach. Variations in annual sulfate export were governed primarily by differences in runoff rather than by differences in sulfate deposition or in total annual precipitation. However, 2 years with similar total runoff had substantially different sulfate export. A residual analysis indicated that a shift in the concentration‐discharge relationship accounted for 72% of the difference in sulfate export and that a change in the pattern of discharge accounted for the remainder of the difference. Stream water sulfate concentrations reflected past hydrologic conditions. Concentrations at the same discharge were higher following an extended dry period than following average periods. The elevation in stream water sulfate concentrations following dry periods persisted for several months. The influence of rainfall patterns on sulfate export underscores the need for long‐term records to adequately characterize the acidification status of the watershed and to understand trends in water quality.

Georgia↗
USGS2010

Interpretation of surface flux measurements in heterogeneous terrain during the Monsoon '90 experiment

A network of 9-m-tall surface flux measurement stations were deployed at eight sparsely vegetated sites during the Monsoon '90 experiment to measure net radiation, Q , soil heat flux, G , sensible heat flux, H (using eddy correlation), and latent heat flux, λ E (using the energy balance equation). At four of these sites, 2-m-tall eddy correlation systems were used to measure all four fluxes directly. Also a 2-m-tall Bowen ratio system was deployed at one site. Magnitudes of the energy balance closure ( Q + G + H + λ E ) increased as the complexity of terrain increased. The daytime Bowen ratio decreased from about 10 before the monsoon season to about 0.3 during the monsoons. Source areas of the measurements are developed and compared to scales of heterogeneity arising from the sparse vegetation and the topography. There was very good agreement among simultaneous measurements of Q with the same model sensor at different heights (representing different source areas), but poor agreement among different brands of sensors. Comparisons of simultaneous measurements of G suggest that because of the extremely small source area, extreme care in sensor deployment is necessary for accurate measurement in sparse canopies. A recently published model to estimate fetch is used to interpret measurements of H at the 2 m and 9 m heights. Three sites were characterized by undulating topography, with ridgetops separated by about 200–600 m. At these sites, sensors were located on ridgetops, and the 9-m fetch included the adjacent valley, whereas the 2-m fetch was limited to the immediate ridgetop and hillside. Before the monsoons began, vegetation was mostly dormant, the watershed was uniformly hot and dry, and the two measurements of H were in close agreement. After the monsoons began and vegetation fully matured, the 2-m measurements of H were significantly greater than the 9-m measurements, presumably because the vegetation in the valleys was denser and cooler than on the ridgetops and hillsides. At one lowland site with little topographic relief, the vegetation was more uniform, and the two measurements of H were in close agreement during peak vegetation. Values of λ E could only be compared at two sites, but the 9-m values were greater than the 2-m values, suggesting λ E from the dense vegetation in the valleys was greater than elsewhere.

Water Environment Research↗
USGS2010

Surface energy balance estimates at local and regional scales using optical remote sensing from an aircraft platform and atmospheric data collected over semiarid rangelands

Remotely sensed data in the visible, near-infrared, and thermal-infrared wave bands were collected from a low-flying aircraft during the Monsoon '90 field experiment. Monsoon '90 was a multidisciplinary experiment conducted in a semiarid watershed. It had as one of its objectives the quantification of hydrometeorological fluxes during the “monsoon” or wet season. The remote sensing observations along with micrometeprological and atmospheric boundary layer (ABL) data were used to compute the surface energy balance over a range of spatial scales. The procedure involved averaging multiple pixels along transects flown over the meteorological and flux (METFLUX) stations. Average values of the spectral reflectance and thermal-infrared temperatures were computed for pixels of order 10 −1 to 10 1 km in length and were used with atmospheric data for evaluating net radiation ( R n ), soil heat flux ( G ), and sensible ( H ) and latent ( LE ) heat fluxes at these same length scales. The model employs a single-layer resistance approach for estimating H that requires wind speed and air temperature in the ABL and a remotely sensed surface temperature. The values of R n and G are estimated from remote sensing information together with near-surface observations of air temperature, relative humidity, and solar radiation. Finally, LE is solved as the residual term in the surface energy balance equation. Model calculations were compared to measurements from the METFLUX network for three days having different environmental conditions. Average percent differences for the three days between model and the METFLUX estimates of the local fluxes were about 5% for R n , 20% for G and H , and 15% for LE . Larger differences occurred during partly cloudy conditions because of errors in interpreting the remote sensing data and the higher spatial and temporal variation in the energy fluxes. Minor variations in modeled energy fluxes were observed when the pixel size representing the remote sensing inputs changed from 0.2 to 2 km. Regional scale estimates of the surface energy balance using bulk ABL properties for the model parameters and input variables and the 10-km pixel data differed from the METFLUX network averages by about 4% for R n , 10% for G and H , and 15% for LE . Model sensitivity in calculating the turbulent fluxes H and LE to possible variations in key model parameters (i.e., the roughness lengths for heat and momentum) was found to be fairly significant. Therefore the reliability of the methods for estimating key model parameters and potential errors needs further testing over different ecosystems and environmental conditions.

Water Resources Research↗
USGS2010

Use of ground-based remotely sensed data for surface energy balance evaluation of a semiarid rangeland

An interdisciplinary field experiment was conducted to study the water and energy balance of a semiarid rangeland watershed in southeast Arizona during the summer of 1990. Two subwatersheds, one grass dominated and the other shrub dominated, were selected for intensive study with ground-based remote sensing systems and hydrometeorological instrumentation. Surface energy balance was evaluated at both sites using direct and indirect measurements of the turbulent fluxes (eddy correlation, variance, and Bowen ratio methods) and using an aerodynamic approach based on remote measurements of surface reflectance and temperature and conventional meteorological information. Estimates of net radiant flux density ( R n ), derived from measurements of air temperature, incoming solar radiation, and surface temperature and radiance compared well with values measured using a net radiometer (mean absolute difference (MAD) ≃ 50 W/m 2 over a range from 115 to 670 W/m 2 ). Soil heat flux density ( G ) was estimated using a relation between G / R n and a spectral vegetation index computed from the red and near-infrared surface reflectance. These G estimates compared well with conventional measurements of G using buried soil heat flux plates (MAD ≃ 20 W/m 2 over a range from −13 to 213 W/m 2 ). In order to account for the effects of sparse vegetation, semiempirical adjustments to the single-layer bulk aerodynamic resistance approach were required for evaluation of sensible heat flux density ( H ). This yielded differences between measurements and remote estimates of H of approximately 33 W/m 2 over a range from 13 to 303 W/m 2 . The resulting estimates of latent heat flux density, LE , were of the same magnitude and trend as measured values; however, a significant scatter was still observed: MAD ≃ 40 W/m 2 over a range from 0 to 340 W/m 2 . Because LE was solved as a residual, there was a cumulative effect of errors associated with remote estimates of R n , G , and H .

Water Resources Research↗
USGS2010

Local energy flux estimates for unstable conditions using variance data in semiarid rangelands

A network of meteorological stations was installed during the Monsoon '90 field campaign in the Walnut Gulch experimental watershed. The study area has a fairly complex surface. The vegetation cover is heterogeneous and sparse, and the terrain is mildly hilly, but dissected by ephemeral channels. Besides measurement of some of the standard weather data such as wind speed, air temperature, and solar radiation, these sites also contained instruments for estimating the local surface energy balance. The approach utilized measurements of net radiation ( R n ), soil heat flux ( G ) and Monin-Obukhov similarity theory applied to first- and second-order turbulent statistics of wind speed and temperature for determining the sensible heat flux ( H ). The latent heat flux ( LE ) was solved as a residual in the surface energy balance equation, namely, LE = −( R n + G + H ). This procedure (VAR-RESID) for estimating the energy fluxes satisfied monetary constraints and the requirement for low maintenance and continued operation through the harsh environmental conditions experienced in semiarid regions. Comparison of energy fluxes using this approach with more traditional eddy correlation techniques showed differences were within 20% under unstable conditions. Similar variability in flux estimates over the study area was present in the eddy correlation data. Hence, estimates of H and LE using the VAR-RESID approach under unstable conditions were considered satisfactory. Also, with second-order statistics of vertical velocity collected at several sites, the local momentum roughness length was estimated. This is an important parameter used in modeling the turbulent transfer of momentum and sensible heat fluxes across the surface-atmosphere interface.

Water Resources Research↗
USGS2010

Movements of water, solutes, and stable isotopes in the unsaturated zones of two sand plains in the upper Midwest

Four month-long field experiments investigated movements of water and solutes through unsaturated sand plains near Princeton, Minnesota, and Oakes, North Dakota. Atrazine and bromide were applied to bare soils and soils planted with corn. The field plots were irrigated according to local farming practices. At the end of each experiment, unsaturated soils were analyzed for atrazine and bromide concentrations and oxygen and hydrogen isotope compositions of soil water. Most soil water was affected by evaporation but groundwater beneath the plots had no evaporative isotopic signature. Therefore most recharge consisted of water that was unaffected by evaporation. Sources of such water may have included snowmelt, prolonged or high-intensity rainfalls that were not interrupted by periods of drying, and water that moved through preferential flow paths. Preferential flow also was suggested by the detection of atrazine, deethylatrazine, and bromide in groundwater shortly after each application of irrigation water at Princeton and by isolated concentrations of atrazine and bromide in soil well below the main masses of chemicals at Oakes.

Minnesota, North Dakota↗
USGS2010

Modeling of soil water retention from saturation to oven dryness

Most analytical formulas used to model moisture retention in unsaturated porous media have been developed for the wet range and are unsuitable for applications in which low water contents are important. We have developed two models that fit the entire range from saturation to oven dryness in a practical and physically realistic way with smooth, continuous functions that have few parameters. Both models incorporate a power law and a logarithmic dependence of water content on suction, differing in how these two components are combined. In one model, functions are added together (model “sum”); in the other they are joined smoothly together at a discrete point (model “junction”). Both models also incorporate recent developments that assure a continuous derivative and force the function to reach zero water content at a finite value of suction that corresponds to oven dryness. The models have been tested with seven sets of water retention data that each cover nearly the entire range. The three-parameter sum model fits all data well and is useful for extrapolation into the dry range when data for it are unavailable. The two-parameter junction model fits most data sets almost as well as the sum model and has the advantage of being analytically integrable for convenient use with capillary-bundle models to obtain the unsaturated hydraulic conductivity.

Water Resources Research↗