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 865 records · Page 48Linked to original sources

Dynamics of CFCs in northern temperate lakes and adjacent groundwater

Three dimictic lakes and one meromictic lake in and near the Trout Lake, Wisconsin, watershed were sampled to determine the variation of chlorofluorocarbon (CFC) concentrations within the lakes. The lakes were sampled during stratified conditions, during fall turnover, and during ice cover. The results demonstrate a considerable variation in CFC concentrations and corresponding atmospheric mixing ratios in the lakes sampled, both with depth and season within a given lake, and across different lakes. CFC profiles and observed degradation were not related to the groundwater inflow rate and hence are likely the result of in‐lake processes influenced by CFC degradation in the (lake) water column, CFC degradation in the lake‐bed sediments, and gas exchange rates and the duration of turnover (turnover efficiency).

Water Resources Research↗

Time series and recurrence interval models to predict the vulnerability of streams to episodic acidification in Shenandoah National Park, Virginia

Acid rain affects headwater streams by temporarily reducing the acid‐neutralizing capacity (ANC) of the water, a process termed episodic acidification. The increase in acidic components in stream water can have deleterious effects on the aquatic biota. Although acidic deposition is uniform across Shenandoah National Park (SNP) in north central Virginia, the stream water quality response during rain events varies substantially. This response is a function of the catchment's underlying geology and topography. Geologic and topographic data for SNP's 231 catchments are readily available; however, long‐term measurements (tens of years) of ANC and accompanying discharge are not and would be prohibitively expensive to collect. Transfer function time series models were developed to predict hourly ANC from discharge for five SNP catchments with long‐term water‐quality and discharge records. Hourly ANC predictions over short time periods (≤1 week) were averaged, and distributions of the recurrence intervals of annual water‐year minimum ANC values were model‐simulated for periods of 6, 24, 72, and 168 hours. The distributions were extrapolated to the rest of the SNP catchments on the basis of catchment geology and topography. On the basis of the models, large numbers of SNP streams have 6‐ to 168‐hour periods of low‐ANC values, which may stress resident fish populations. Smaller catchments are more vulnerable to episodic acidification than larger catchments underlain by the same bedrock. Catchments with similar topography and size are more vulnerable if underlain by less basaltic/carbonate bedrock. Many catchments are predicted to have successive years of low‐ANC values potentially sufficient to extirpate some species.

Virginia↗

Stochastic uncertainty analysis for unconfined flow systems

A new stochastic approach proposed by Zhang and Lu (2004), called the Karhunen‐Loeve decomposition‐based moment equation (KLME), has been extended to solving nonlinear, unconfined flow problems in randomly heterogeneous aquifers. This approach is on the basis of an innovative combination of Karhunen‐Loeve decomposition, polynomial expansion, and perturbation methods. The random log‐transformed hydraulic conductivity field ( lnK S ) is first expanded into a series in terms of orthogonal Gaussian standard random variables with their coefficients obtained as the eigenvalues and eigenfunctions of the covariance function of lnK S . Next, head h is decomposed as a perturbation expansion series Σ h ( m ) , where h ( m ) represents the m th‐order head term with respect to the standard deviation of lnK S . Then h ( m ) is further expanded into a polynomial series of m products of orthogonal Gaussian standard random variables whose coefficients h i 1 , i 2 ,..., i m ( m ) are deterministic and solved sequentially from low to high expansion orders using MODFLOW‐2000. Finally, the statistics of head and flux are computed using simple algebraic operations on h i 1 , i 2 ,..., i m ( m ) . A series of numerical test results in 2‐D and 3‐D unconfined flow systems indicated that the KLME approach is effective in estimating the mean and (co)variance of both heads and fluxes and requires much less computational effort as compared to the traditional Monte Carlo simulation technique.

Water Resources Research↗

Quantifying surface water–groundwater interactions using time series analysis of streambed thermal records: Method development

We present a method for determining streambed seepage rates using time series thermal data. The new method is based on quantifying changes in phase and amplitude of temperature variations between pairs of subsurface sensors. For a reasonable range of streambed thermal properties and sensor spacings the time series method should allow reliable estimation of seepage rates for a range of at least ±10 m d −1 (±1.2 × 10 −2 m s −1 ), with amplitude variations being most sensitive at low flow rates and phase variations retaining sensitivity out to much higher rates. Compared to forward modeling, the new method requires less observational data and less setup and data handling and is faster, particularly when interpreting many long data sets. The time series method is insensitive to streambed scour and sedimentation, which allows for application under a wide range of flow conditions and allows time series estimation of variable streambed hydraulic conductivity. This new approach should facilitate wider use of thermal methods and improve understanding of the complex spatial and temporal dynamics of surface water–groundwater interactions.

Water Resources Research↗

Reach‐averaged sediment routing model of a canyon river

Spatial complexity in channel geometry indicates that accurate prediction of sediment transport requires modeling in at least two dimensions. However, a one‐dimensional model may be the only practical or possible alternative, especially for longer river reaches of practical concern in river management or landscape modeling. We have developed a one‐dimensional model of the Colorado River through upper Grand Canyon that addresses this problem by reach averaging the channel properties and predicting changes in sand storage using separate source and sink functions coupled to the sand routing model. The model incorporates results from the application of a two‐dimensional model of flow, sand transport, and bed evolution, and a new algorithm for setting the near‐bed sand boundary condition for sand transported over an exposed bouldery bed. Model predictions were compared to measurements of sand discharge during intermittent tributary inputs and varying discharges controlled by dam releases. The model predictions generally agree well with the timing and magnitude of measured sand discharges but tend to overpredict sand discharge during the early stages of a high release designed to redistribute sand to higher‐elevation deposits.

Water Resources Research↗

Modeling the movement of a pH perturbation and its impact on adsorbed zinc and phosphate in a wastewater‐contaminated aquifer

[1] Chemical conditions were perturbed in an aquifer with an ambient pH of 5.9 and wastewater-derived adsorbed zinc (Zn) and phosphate (P) contamination by injecting a pulse of amended groundwater. The injected groundwater had low concentrations of dissolved Zn and P, a pH value of 4.5 resulting from equilibration with carbon dioxide gas, and added potassium bromide (KBr). Downgradient of the injection, breakthrough of nonreactive Br and total dissolved carbonate concentrations in excess of ambient values (excess TCO 2 ) were accompanied by a decrease in pH values and over twentyfold increases in dissolved Zn concentrations above preinjection values. Peak concentrations of Br and excess TCO 2 were followed by slow increases in pH values accompanied by significant increases in dissolved P above preinjection concentrations. The injected tracers mobilized a significant mass of wastewater-derived Zn. Reactive transport simulations incorporating surface complexation models for adsorption of Zn, P, hydrogen ions, and major cations onto the aquifer sediments, calibrated using laboratory experimental data, captured most of the important trends observed during the experiment. These include increases in Zn concentrations in response to the pH perturbation, perturbations in major cation concentrations, attenuation of the pH perturbation with transport distance, and increases in alkalinity with transport distance. Observed desorption of P in response to chemical perturbations was not predicted, possibly because of a disparity between the range of chemical conditions in the calibration data set and those encountered during the field experiment. Zinc and P desorbed rapidly in response to changing chemical conditions despite decades of contact with the sediments. Surface complexation models with relatively few parameters in the form of logK values and site concentrations show considerable promise for describing the influence of variable chemistry on the transport of adsorbing contaminants.

Water Resources Research↗

Debris flow deposition and reworking by the Colorado River in Grand Canyon, Arizona

Flow regulation by large dams affects downstream flow competence and channel maintenance. Debris flows from 740 tributaries in Grand Canyon, Arizona, transport coarse‐grained sediment onto debris fans adjacent to the Colorado River. These debris fans constrict the river to form rapids and are reworked during river flows that entrain particles and transport them downstream. Beginning in 1963, flood control operations of Glen Canyon Dam limited the potential for reworking of aggraded debris fans. We analyzed change in debris fans at the mouths of 75‐Mile and Monument Creeks using photogrammetry of aerial photography taken from 1965 to 2000 and supplemented with ground surveys performed from 1987 to 2005. Our results quantify the debris fan aggradation that resulted from debris flows from 1984 to 2003. Volume, area, and river constriction increased at both debris fans. Profiles of the two debris fans show that net aggradation occurred in the middle of debris fans at stages above maximum dam releases, and surface shape shifted from concave to convex. Dam releases above power plant capacity partially reworked both debris fans, although reworking removed much less sediment than what was added by debris flow deposition. Large dam releases would be required to create additional reworking to limit the rate of debris fan aggradation in Grand Canyon.

Arizona↗

Potential effects of regional pumpage on groundwater age distribution

Groundwater ages estimated from environmental tracers can help calibrate groundwater flow models. Groundwater age represents a mixture of traveltimes, with the distribution of ages determined by the detailed structure of the flow field, which can be prone to significant transient variability. Effects of pumping on age distribution were assessed using direct age simulation in a hypothetical layered aquifer system. A steady state predevelopment age distribution was computed first. A well field was then introduced, and pumpage caused leakage into the confined aquifer of older water from an overlying confining unit. Large changes in simulated groundwater ages occurred in both the aquifer and the confining unit at high pumping rates, and the effects propagated a substantial distance downgradient from the wells. The range and variance of ages contributing to the well increased substantially during pumping. The results suggest that the groundwater age distribution in developed aquifers may be affected by transient leakage from low‐permeability material, such as confining units, under certain hydrogeologic conditions.

Water Resources Research↗

Impacts of the 2004 tsunami on groundwater resources in Sri Lanka

The 26 December 2004 tsunami caused widespread destruction and contamination of coastal aquifers across southern Asia. Seawater filled domestic open dug wells and also entered the aquifers via direct infiltration during the first flooding waves and later as ponded seawater infiltrated through the permeable sands that are typical of coastal aquifers. In Sri Lanka alone, it is estimated that over 40,000 drinking water wells were either destroyed or contaminated. From February through September 2005, a team of United States, Sri Lankan, and Danish water resource scientists and engineers surveyed the coastal groundwater resources of Sri Lanka to develop an understanding of the impacts of the tsunami and to provide recommendations for the future of coastal water resources in south Asia. In the tsunami‐affected areas, seawater was found to have infiltrated and mixed with fresh groundwater lenses as indicated by the elevated groundwater salinity levels. Seawater infiltrated through the shallow vadose zone as well as entered aquifers directly through flooded open wells. Our preliminary transport analysis demonstrates that the intruded seawater has vertically mixed in the aquifers because of both forced and free convection. Widespread pumping of wells to remove seawater was effective in some areas, but overpumping has led to upconing of the saltwater interface and rising salinity. We estimate that groundwater recharge from several monsoon seasons will reduce salinity of many sandy Sri Lankan coastal aquifers. However, the continued sustainability of these small and fragile aquifers for potable water will be difficult because of the rapid growth of human activities that results in more intensive groundwater pumping and increased pollution. Long‐term sustainability of coastal aquifers is also impacted by the decrease in sand replenishment of the beaches due to sand mining and erosion.

Water Resources Research↗

Modeling the influence of river rehabilitation scenarios on bed material sediment flux in a large river over decadal timescales

A stochastic flood generator and calibrated sediment transport formulae were used to assess the decadal impact of major river rehabilitation strategies on two fraction bed material sediment flux and net storage, first‐order indicators of aquatic riverine habitat, in a large river system. Model boundary conditions were modified to reflect the implementation of three major river rehabilitation strategies being considered in the Sacramento River Valley: gravel augmentation, setting back of levees, and flow alteration. Fifty 30‐year model simulations were used to compute probabilities of the response in sediment flux and net storage to these strategies. Total annual average bed material sediment flux estimates were made at six gauged river cross sections, and ∼60 km reach‐scale sediment budgets were evaluated between them. Gravel augmentation to improve spawning habitat induced gravel accumulation locally and/or downstream, depending on the added mixture. Levee setbacks to recreate the river corridor reduced flow stages for most flows and hence lowered sediment flux. Flow alteration to mimic natural flow regimes systematically decreased total annual average flux, suggesting that high‐magnitude low‐frequency transport events do not affect long‐term trends in bed material flux. The results indicate that each rehabilitation strategy reduces sediment transport in its target reaches and modulates imbalances in total annual bed material sediment budgets at the reach scale. Additional risk analysis is necessary to identify extreme conditions associated with variable hydrology that could affect rehabilitation over decades. Sensitivity analysis suggests that sorting of bed material sediment is the most important determinant of modeled transport and storage patterns.

Water Resources Research↗

A field assessment of the value of steady shape hydraulic tomography for characterization of aquifer heterogeneities

Hydraulic tomography is a promising approach for obtaining information on variations in hydraulic conductivity on the scale of relevance for contaminant transport investigations. This approach involves performing a series of pumping tests in a format similar to tomography. We present a field‐scale assessment of hydraulic tomography in a porous aquifer, with an emphasis on the steady shape analysis methodology. The hydraulic conductivity (K) estimates from steady shape and transient analyses of the tomographic data compare well with those from a tracer test and direct‐push permeameter tests, providing a field validation of the method. Zonations based on equal‐thickness layers and cross‐hole radar surveys are used to regularize the inverse problem. The results indicate that the radar surveys provide some useful information regarding the geometry of the K field. The steady shape analysis provides results similar to the transient analysis at a fraction of the computational burden. This study clearly demonstrates the advantages of hydraulic tomography over conventional pumping tests, which provide only large‐scale averages, and small‐scale hydraulic tests (e.g., slug tests), which cannot assess strata connectivity and may fail to sample the most important pathways or barriers to flow.

Water Resources Research↗

Water use regimes: Characterizing direct human interaction with hydrologic systems

The sustainability of human water use practices is a rapidly growing concern in the United States and around the world. To better characterize direct human interaction with hydrologic systems (stream basins and aquifers), we introduce the concept of the water use regime. Unlike scalar indicators of anthropogenic hydrologic stress in the literature, the water use regime is a two‐dimensional, vector indicator that can be depicted on simple x‐y plots of normalized human withdrawals ( h out ) versus normalized human return flows ( h in ). Four end‐member regimes, natural‐flow‐dominated (undeveloped), human‐flow‐dominated (churned), withdrawal‐dominated (depleted), and return‐flow‐dominated (surcharged), are defined in relation to limiting values of h out and h in . For illustration, the water use regimes of 19 diverse hydrologic systems are plotted and interpreted. Several of these systems, including the Yellow River Basin, China, and the California Central Valley Aquifer, are shown to approach particular end‐member regimes. Spatial and temporal regime variations, both seasonal and long‐term, are depicted. Practical issues of data availability and regime uncertainty are addressed in relation to the statistical properties of the ratio estimators h out and h in . The water use regime is shown to be a useful tool for comparative water resources assessment and for describing both historic and alternative future pathways of water resource development at a range of scales.

Water Resources Research↗

Controls on the variability of net infiltration to desert sandstone

As populations grow in arid climates and desert bedrock aquifers are increasingly targeted for future development, understanding and quantifying the spatial variability of net infiltration becomes critically important for accurately inventorying water resources and mapping contamination vulnerability. This paper presents a conceptual model of net infiltration to desert sandstone and then develops an empirical equation for its spatial quantification at the watershed scale using linear least squares inversion methods for evaluating controlling parameters (independent variables) based on estimated net infiltration rates (dependent variables). Net infiltration rates used for this regression analysis were calculated from environmental tracers in boreholes and more than 3000 linear meters of vadose zone excavations in an upland basin in southwestern Utah underlain by Navajo sandstone. Soil coarseness, distance to upgradient outcrop, and topographic slope were shown to be the primary physical parameters controlling the spatial variability of net infiltration. Although the method should be transferable to other desert sandstone settings for determining the relative spatial distribution of net infiltration, further study is needed to evaluate the effects of other potential parameters such as slope aspect, outcrop parameters, and climate on absolute net infiltration rates.

Water Resources Research↗

Stochastic uncertainty analysis for solute transport in randomly heterogeneous media using a Karhunen‐Loève‐based moment equation approach

A new approach has been developed for solving solute transport problems in randomly heterogeneous media using the Karhunen‐Loève‐based moment equation (KLME) technique proposed by Zhang and Lu (2004). The KLME approach combines the Karhunen‐Loève decomposition of the underlying random conductivity field and the perturbative and polynomial expansions of dependent variables including the hydraulic head, flow velocity, dispersion coefficient, and solute concentration. The equations obtained in this approach are sequential, and their structure is formulated in the same form as the original governing equations such that any existing simulator, such as Modular Three‐Dimensional Multispecies Transport Model for Simulation of Advection, Dispersion, and Chemical Reactions of Contaminants in Groundwater Systems (MT3DMS), can be directly applied as the solver. Through a series of two‐dimensional examples, the validity of the KLME approach is evaluated against the classical Monte Carlo simulations. Results indicate that under the flow and transport conditions examined in this work, the KLME approach provides an accurate representation of the mean concentration. For the concentration variance, the accuracy of the KLME approach is good when the conductivity variance is 0.5. As the conductivity variance increases up to 1.0, the mismatch on the concentration variance becomes large, although the mean concentration can still be accurately reproduced by the KLME approach. Our results also indicate that when the conductivity variance is relatively large, neglecting the effects of the cross terms between velocity fluctuations and local dispersivities, as done in some previous studies, can produce noticeable errors, and a rigorous treatment of the dispersion terms becomes more appropriate.

Water Resources Research↗

Carbon export and cycling by the Yukon, Tanana, and Porcupine rivers, Alaska, 2001-2005

Loads and yields of dissolved and particulate organic and inorganic carbon (DOC, POC, DIC, PIC) were measured and modeled at three locations on the Yukon River (YR) and on the Tanana and Porcupine rivers (TR, PR) in Alaska during 2001–2005. Total YR carbon export averaged 7.8 Tg C yr −1 , 30% as OC and 70% as IC. Total C yields (0.39–1.03 mol C m −2 yr −1 ) were proportional to water yields (139–356 mm yr −1 ; r 2 = 0.84) at all locations. Summer DOC had an aged component (fraction modern (FM) = 0.94–0.97), except in the permafrost wetland‐dominated PR, where DOC was modern. POC had FM = 0.63–0.70. DOC had high concentration, high aromaticity, and high hydrophobic content in spring and low concentration, low aromaticity, and high hydrophilic content in winter. About half of annual DOC export occurred during spring. DIC concentration and isotopic composition were strongly affected by dissolution of suspended carbonates in glacial meltwater during summer.

Water Resources Research↗

Spatial and temporal variability in sedimentation rates associated with cutoff channel infill deposits: Ain River, France

Floodplain development is associated with lateral accretion along stable channel geometry. Along shifting rivers, the floodplain sedimentation is more complex because of changes in channel position but also cutoff channel presence, which exhibit specific overflow patterns. In this contribution, the spatial and temporal variability of sedimentation rates in cutoff channel infill deposits is related to channel changes of a shifting gravel bed river (Ain River, France). The sedimentation rates estimated from dendrogeomorphic analysis are compared between and within 14 cutoff channel infills. Detailed analyses along a single channel infill are performed to assess changes in the sedimentation rates through time by analyzing activity profiles of the fallout radionuclides 137 Cs and unsupported 210 Pb. Sedimentation rates are also compared within the channel infills with rates in other plots located in the adjacent floodplain. Sedimentation rates range between 0.65 and 2.4 cm a −1 over a period of 10 to 40 years. The data provide additional information on the role of distance from the bank, overbank flow frequency, and channel geometry in controlling the sedimentation rate. Channel infills, lower than adjacent floodplains, exhibit higher sedimentation rates and convey overbank sediment farther away within the floodplain. Additionally, channel degradation, aggradation, and bank erosion, which reduce or increase the distance between the main channel and the cutoff channel aquatic zone, affect local overbank flow magnitude and frequency and therefore sedimentation rates, thereby creating a complex mosaic of sedimentation zones within the floodplain and along the cutoff channel infills. Last, the dendrogeomorphic and 137 Cs approaches are cross validated for estimating the sedimentation rate within a channel infill.

Ain River↗

Relating low‐flow characteristics to the base flow recession time constant at partial record stream gauges

Base flow recession information is helpful for regional estimation of low‐flow characteristics. However, analyses that exploit such information generally require a continuous record of streamflow at the estimation site to characterize base flow recession. Here we propose a simple method for characterizing base flow recession at low‐flow partial record stream gauges (i.e., sites with very few streamflow measurements under low‐streamflow conditions), and we use that characterization as the basis for a practical new approach to low‐flow regression. In a case study the introduction of a base flow recession time constant, estimated from a single pair of strategically timed streamflow measurements, approximately halves the root‐mean‐square estimation error relative to that of a conventional drainage area regression. Additional streamflow measurements can be used to reduce the error further.

Water Resources Research↗

Groundwater noble gas, age, and temperature signatures in an Alpine watershed: Valuable tools in conceptual model development

Bedrock groundwater in alpine watersheds is poorly understood, mainly because of a scarcity of wells in alpine settings. Groundwater noble gas, age, and temperature data were collected from springs and wells with depths of 3–342 m in Handcart Gulch, an alpine watershed in Colorado. Temperature profiles indicate active groundwater circulation to a maximum depth (aquifer thickness) of about 200 m, or about 150 m below the water table. Dissolved noble gas data show unusually high excess air concentrations (>0.02 cm 3 STP/g, ΔNe > 170%) in the bedrock, consistent with unusually large seasonal water table fluctuations (up to 50 m) observed in the upper part of the watershed. Apparent 3 H/ 3 He ages are positively correlated with sample depth and excess air concentrations. Integrated samples were collected from artesian bedrock wells near the trunk stream and are assumed to approximate flow‐weighted samples reflecting bedrock aquifer mean residence times. Exponential mean ages for these integrated samples are remarkably consistent along the stream, four of five being from 8 to 11 years. The tracer data in combination with other hydrologic and geologic data support a relatively simple conceptual model of groundwater flow in the watershed in which (1) permeability is primarily a function of depth; (2) water table fluctuations increase with distance from the stream; and (3) recharge, aquifer thickness, and porosity are relatively uniform throughout the watershed in spite of the geological complexity of the Proterozoic crystalline rocks that underlie it.

Water Resources Research↗