Search USGSSearch

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,045 records · Page 58Linked to original sources

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

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

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

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

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

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

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

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

Transport of chromium and selenium in the suboxic zone of a shallow aquifer: Influence of redox and adsorption reactions

Breakthrough of Cr(VI) (chromate), Se(VI) (selenate), and O 2 (dissolved oxygen) was observed in tracer tests conducted in a shallow, sand and gravel aquifer with mildly reducing conditions. Loss of Cr, probably due to reduction of Cr(VI) to Cr(III) and irreversible sorption of Cr(III), occurred along with slight retardation of Cr(VI), owing to reversible sorption. Reduction of Se(VI) and O 2 was thermodynamically feasible but did not occur, indicating conditions, were unfavorable to microbial reduction. Cr(VI) reduction by constituents of aquifer sediments did not achieve local equilibrium during transport. The reduction rate was probably limited by incomplete contact between Cr(VI) transported along predominant flow paths and reductants located in regions within aquifer sediments of comparatively low permeability. Scatter in the amount of Cr reduction calculated from individual breakthrough curves at identical distances downgradient probably resulted from heterogeneities in the distribution of reductants in the sediments. Predictive modeling of the transport and fate of redox-sensitive solutes cannot be based strictly on thermodynamic considerations; knowledge of reaction rates is critical. Potentially important mass transfer rate limitations between solutes and reactants in sediments as well as heterogeneities in the distribution of redox properties in aquifers complicate determination of limiting rates for use in predictive simulations of the transport of redox-sensitive contaminants in groundwater.

Water Resources Research

Dating of shallow groundwater: Comparison of the transient tracers 3H/3He, chlorofluorocarbons, and 85Kr

This paper describes a direct comparison of apparent ages derived from 3 H/ 3 He, chlorofluorocarbons (CCl 3 F and CCl 2 F 2 ), and 85 Kr measurements in shallow groundwater. Wells chosen for this study are completed in the unconfined surficial aquifers in late Cenozoic Atlantic Coastal Plain sediments of the Delmarva Peninsula, on the east coast of the United States. Most of the apparent tracer ages agree within 2 years of each other for recharge dates between 1965 and 1990. Discrepancies in apparent tracer ages usually can be explained by hydrological processes such as mixing in a discharge area. Recharge rate calculations based on apparent tracer age gradients at multilevel well locations agree with previous recharge estimates. High recharge rates on the Delmarva Peninsula result in nearly complete dissolved-gas confinement in the groundwater. The remarkable agreement between the different tracer ages indicates negligible mixing of waters of different ages, insignificant dispersion, minimal gas loss to the atmosphere, and insignificant sorption-desorption processes at this location.

Water Resources Research

Marginal bed load transport in a gravel bed stream, Sagehen Creek, California

Marginal bed load transport describes the condition when relatively few bed particles are moving at any time. Bed particles resting in the shallowest bed pockets will move when the dimensionless shear stress т * exceeds a value of about 0.020. As т * increases, the number of bed particles moving increases. Significant motion of bed particles, i.e., when a substantial fraction of the bed particles are moving, occurs when т * exceeds a value of about 0.060. Thus marginal bed load transport occurs over the domain 0.020 < т * < 0.060. Marginal bed load transport rates and associated hydraulic characteristics of Sagehen Creek, a small mountain gravel bed stream, were measured on 55 days at discharges ranging from slightly less than one half of the bank-full discharge to more than 4 times the bank-full discharge. Dimensionless shear stress varied from 0.032 to 0.042, and bed particles as large as the 80th percentile of the bed surface were transported. The relation between reference dimensionless shear stress and relative particle protrusion for Sagehen Creek was determined by varying т * ri to obtain the best fit of the Parker bed load function to the measured transport rates. During the period of record (water years 1954–1991), the mean annual quantity of bed load transported past the Sagehen Creek gage was 24.7 tons. Forty-seven percent of all bed load transported during the 38 years of record occurred in just 6 years. During 10 of the 38 years of record, essentially no bed load was transported. The median diameter of bed load was 26 mm, compared to 58 mm in the surface bed material.

California

Climate, soil water storage, and the average annual water balance

This paper describes the development and testing of the hypothesis that the long-term water balance is determined only by the local interaction of fluctuating water supply (precipitation) and demand (potential evapotranspiration), mediated by water storage in the soil. Adoption of this hypothesis, together with idealized representations of relevant input variabilities in time and space, yields a simple model of the water balance of a finite area having a uniform climate. The partitioning of average annual precipitation into evapotranspiration and runoff depends on seven dimensionless numbers: the ratio of average annual potential evapotranspiration to average annual precipitation (index of dryness); the ratio of the spatial average plant-available water-holding capacity of the soil to the annual average precipitation amount; the mean number of precipitation events per year; the shape parameter of the gamma distribution describing spatial variability of storage capacity; and simple measures of the seasonality of mean precipitation intensity, storm arrival rate, and potential evapotranspiration. The hypothesis is tested in an application of the model to the United States east of the Rocky Mountains, with no calibration. Study area averages of runoff and evapotranspiration, based on observations, are 263 mm and 728 mm, respectively; the model yields corresponding estimates of 250 mm and 741 mm, respectively, and explains 88% of the geographical variance of observed runoff within the study region. The differences between modeled and observed runoff can be explained by uncertainties in the model inputs and in the observed runoff. In the humid (index of dryness <1) parts of the study area, the dominant factor producing runoff is the excess of annual precipitation over annual potential evapotranspiration, but runoff caused by variability of supply and demand over time is also significant; in the arid (index of dryness >1) parts, all of the runoff is caused by variability of forcing over time. Contributions to model runoff attributable to small-scale spatial variability of storage capacity are insignificant throughout the study area. The consistency of the model with observational data is supportive of the supply-demand-storage hypothesis, which neglects infiltration excess runoff and other finite-permeability effects on the soil water balance.

Water Resources Research

Reply to “Comment on ‘An efficient numerical solution of the transient storage equations for solute transport in small streams" by R. L. Runkel and S. C. Chapra

We thank Dawes and Short [this issue] for presenting an alternate technique for the efficient solution of the transient storage solute transport equations. After reading their comment, it is clear that several points in our original manuscript are in need of clarification. In this reply we provide additional information on our solution technique and comment briefly on the alternate scheme of Dawes and Short.

Water Resources Research

Flow path studies in forested watersheds of headwater tributaries of Brush Brook, Vermont

An investigation was undertaken into how headwater tributaries of Brush Brook, Vermont, could have average p H differences of almost two units (4.75 and 6.7). Sampling along four tributaries revealed that most of one tributary, below an area of seeps, had consistently higher p H, Ca 2+ , Mg 2+ , and K + , and lower Al than other sites. Bedrock mapping showed numerous fractures in vicinity of the seeps. A portion of this tributary's watershed and a portion of an acid tributary's watershed were intensively mapped for soil depth. Sampling showed the widespread existence of dense basal till in the watershed of the acid tributary but none in that of the near-neutral stream. Lateral flow, found above the dense till, was chemically similar to that of the acid tributary and to solutions sampled from soil B horizons. There were no differences in the average p H of nonseep soils sampled from either watershed. Flow paths are hypothesized to be through the B horizons in the acid tributaries and from below the soil profile in the near-neutral tributary. The acid catchment should be more sensitive to environmental change.

Vermont

Influence of diurnal variations in stream temperature on streamflow loss and groundwater recharge

We demonstrate that for losing reaches with significant diurnal variations in stream temperature, the effect of stream temperature on streambed seepage is a major factor contributing to reduced afternoon streamflows. An explanation is based on the effect of stream temperature on the hydraulic conductivity of the streambed, which can be expected to double in the 0° to 25°C temperature range. Results are presented for field experiments in which stream discharge and temperature were continuously measured for several days over losing reaches at St. Kevin Gulch, Colorado, and Tijeras Arroyo, New Mexico. At St. Kevin Gulch in July 1991, the diurnal stream temperature in the 160-m study reach ranged from about 4° to 18°C, discharges ranged from 10 to 18 L/s, and streamflow loss in the study reach ranged from 2.7 to 3.7 L/s. On the basis of measured stream temperature variations, the predicted change in conductivity was about 38%; the measured change in stream loss was about 26%, suggesting that streambed temperature varied less than the stream temperature. At Tijeras Arroyo in May 1992, diurnal stream temperature in the 655-m study reach ranged from about 10° to 25°C and discharge ranged from 25 to 55 L/s. Streamflow loss was converted to infiltration rates by factoring in the changing stream reach surface area and streamflow losses due to evaporation rates as measured in a hemispherical evaporation chamber. Infiltration rates ranged from about 0.7 to 2.0 m/d, depending on time and location. Based on measured stream temperature variations, the predicted change in conductivity was 29%; the measured change in infiltration was also about 27%. This suggests that high infiltration rates cause rapid convection of heat to the streambed. Evapotranspiration losses were estimated for the reach and adjacent flood plain within the arroyo. On the basis of these estimates, only about 5% of flow loss was consumed via stream evaporation and stream-side evapotranspiration, indicating that 95% of the loss within the study reach represented groundwater recharge.

Water Resources Research

Long-term hydrologic and biogeochemical responses of a soft water seepage lake in north central Wisconsin

Hydrologic and chemical budgets were determined for the period October 1980 through September 1988 for Vandercook Lake, a 43-ha soft water seepage lake located in predominantly sandy outwash in northern Wisconsin. The 1951–1980 mean annual precipitation for the lake, based on nearby National Weather Service (NWS) stations, was 825 mm; volume-weighted p H of bulk atmospheric deposition during the 8-year study averaged 4.7. From October 1980 through September 1983, annual precipitation was 108% of the 30-year NWS average, annual recharge was 144% of the 8-year study period average, and annual groundwater inflow to the lake ranged from 5 to 9% of the total inflow. From October 1986 through September 1988, annual precipitation was 90% of the NWS average, annual recharge was 30% of the study period average, and the lake received no groundwater inflow. During the study, the lake changed from a system whose buffering mechanism was significantly influenced by mineral weathering in inflowing groundwater to one dominated by in-lake sulfate reduction. The functional differences exhibited by this lake during the 8 years of study demonstrate the tenuousness with which conclusions based on shorter-term studies of similar systems must be considered.

Wisconsin

An updated model of induced airflow in the unsaturated zone

Simulation of induced movement of air in the unsaturated zone provides a method to determine permeability and to design vapor extraction remediation systems. A previously published solution to the airflow equation for the case in which the unsaturated zone is separated from the atmosphere by a layer of lower permeability (such as a clay layer) has been superseded. The new solution simulates airflow through the layer of lower permeability more rigorously by defining the leakage in terms of the upper boundary condition rather than by adding a leakage term to the governing airflow equation. This note presents the derivation of the new solution. Formulas for steady state pressure, specific discharge, and mass flow in the domain are obtained for the new model and for the case in which the unsaturated zone is in direct contact with the atmosphere.

Water Resources Research

Deducing the distribution of terminal electron-accepting processes in hydrologically diverse groundwater systems

The distribution of microbially mediated terminal electron-accepting processes (TEAPs( was investigated in four hydrologically diverse groundwater systems by considering patterns of electron acceptor (nitrate, sulfate) consumption, intermediate product (hydrogen (H 2 )) concentrations, and final product (ferrous iron, sulfide, and methane) production. In each hydrologic system a determination of predominant TEAPs could be arrived at, but the level of confidence appropriate for each determination differed. In a portion of the lacustrine aquifer of the San Joaquin Valley, for example, all three indicators (sulfate concentrations decreasing, H 2 concentrations in the 1–2 nmol range, and sulfide concentrations increasing along flow paths identified sulfate reduction as the predominant TEAP, leading to a high degree of confidence in the determination. In portions of the Floridan aquifer and a petroleum hydrocarbon-contaminated aquifer, sulfate reduction and methanogenesis are indicated by production of sulfide and methane, and hydrogen oncentrations in the 1–4 nmol and 5–14 nmol range, respectively. However, because electron acceptor consumption could not be documented in these systems, less confidence is warranted in the TEAP determination. In the Black Creek aquifer, no pattern of sulfate consumption and sulfide production were observed, but H2 concentrations indicated sulfate reduction as the predominant TEAP. In this case, where just a single line of evidence is available, the least confidence in the TEAP diagnosis is justified. Because this methodology is based on measurable water chemistry parameters and upon the physiology of microbial electron transfer processes, it provides a better description of predominant redox processes in groundwater systems than more traditional Eh-based methods.

Water Resources Research