USGS Science⌕ Search

SEARCH · USGS Science

Results for “Water Resources 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 451 records · Page 25Linked to original sources
USGS2006

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Simple predictions of maximum transport rate in unsaturated soil and rock

In contrast with the extreme variability expected for water and contaminant fluxes in the unsaturated zone, evidence from 64 field tests of preferential flow indicates that the maximum transport speed V max , adjusted for episodicity of infiltration, deviates little from a geometric mean of 13 m/d. A model based on constant‐speed travel during infiltration pulses of actual or estimated duration can predict V max with approximate order‐of‐magnitude accuracy, irrespective of medium or travel distance, thereby facilitating such problems as the prediction of worst‐case contaminant traveltimes. The lesser variability suggests that preferential flow is subject to rate‐limiting mechanisms analogous to those that impose a terminal velocity on objects in free fall and to rate‐compensating mechanisms analogous to Le Chatlier's principle. A critical feature allowing such mechanisms to dominate may be the presence of interfacial boundaries confined by neither solid material nor capillary forces.

Water Resources Research↗
USGS2007

Analysis of pumping‐induced unsaturated regions beneath a perennial river

The presence of an unsaturated region beneath a streambed during groundwater pumping near streams can reduce the pumping capacity, change flow paths, and alter the types of biological transformations in the streambed sediments. A three‐dimensional, multiphase flow model of two horizontal collector wells along the Russian River near Forestville, California, was developed to investigate the impact of varying the ratio of the aquifer to streambed permeability on (1) the formation of an unsaturated region beneath the stream, (2) the pumping capacity, (3) stream water fluxes through the streambed, and (4) stream water traveltimes to the collector wells. The aquifer to streambed permeability ratio at which the unsaturated region was initially observed ranged from 10 to 100. The size of the unsaturated region beneath the streambed increased as the aquifer to streambed permeability ratio increased. The simulations also indicated that for a particular aquifer permeability, decreasing the streambed permeability by only a factor of 2–3 from the permeability where desaturation initially occurred resulted in reducing the pumping capacity. In some cases, the stream water fluxes increased as the streambed permeability decreased. However, the stream water residence times increased and the fraction of stream water that reached that the wells decreased as the streambed permeability decreased, indicating that a higher streambed flux does not necessarily correlate to greater recharge of stream water around the wells.

Water Resources Research↗
USGS2007

Ecohydrological controls on soil moisture and hydraulic conductivity within a pinyon‐juniper woodland

The impact of pinyon‐juniper woodland encroachment on rangeland ecosystems is often associated with a reduction of streamflow and recharge and an increase in soil erosion. The objective of this study is to investigate vegetational control on seasonal soil hydrologic properties along a 15‐m transect in pinyon‐juniper woodland with biocrust. We demonstrate that the juniper tree controls soil water content (SWC) patterns directly under the canopy via interception, and beyond the canopy via shading in a preferred orientation, opposite to the prevailing wind direction. The juniper also controls the SWC and unsaturated hydraulic conductivity measured close to water saturation ( K ( h )) under the canopy by the creation of soil water repellency due to needle drop. We use this information to refine the hydrologic functional unit (HFU) concept into three interacting hydrologic units: canopy patches, intercanopy patches, and a transitional unit formed by intercanopy patches in the rain shadow of the juniper tree. Spatial autoregressive state‐space models show the close relationship between K ( h ) close to soil water saturation and SWC at medium and low levels, integrating a number of influences on hydraulic conductivity.

Water Resources Research↗