Geology Reports⌕ Search

SEARCH · Geology Reports

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 253 records · Page 14Linked to original sources

Investigating runoff efficiency in upper Colorado River streamflow over past centuries

With increasing concerns about the impact of warming temperatures on water resources, more attention is being paid to the relationship between runoff and precipitation, or runoff efficiency. Temperature is a key influence on Colorado River runoff efficiency, and warming temperatures are projected to reduce runoff efficiency. Here, we investigate the nature of runoff efficiency in the upper Colorado River (UCRB) basin over the past 400 years, with a specific focus on major droughts and pluvials, and to contextualize the instrumental period. We first verify the feasibility of reconstructing runoff efficiency from tree-ring data. The reconstruction is then used to evaluate variability in runoff efficiency over periods of high and low flow, and its correspondence to a reconstruction of late runoff season UCRB temperature variability. Results indicate that runoff efficiency has played a consistent role in modulating the relationship between precipitation and streamflow over past centuries, and that temperature has likely been the key control. While negative runoff efficiency is most common during dry periods, and positive runoff efficiency during wet years, there are some instances of positive runoff efficiency moderating the impact of precipitation deficits on streamflow. Compared to past centuries, the 20th century has experienced twice as many high flow years with negative runoff efficiency, likely due to warm temperatures. These results suggest warming temperatures will continue to reduce runoff efficiency in wet or dry years, and that future flows will be less than anticipated from precipitation due to warming temperatures.

Colorado River↗

Spatial patterns of development drive water use

Water availability is becoming more uncertain as human populations grow, cities expand into rural regions and the climate changes. In this study, we examine the functional relationship between water use and the spatial patterns of developed land across the rapidly growing region of the southeastern United States. We quantified the spatial pattern of developed land within census tract boundaries, including multiple metrics of density and configuration. Through non‐spatial and spatial regression approaches we examined relationships and spatial dependencies between the spatial pattern metrics, socio‐economic and environmental variables and two water use variables: a) domestic water use, and b) total development‐related water use (a combination of public supply, domestic self‐supply and industrial self‐supply). Metrics describing the spatial patterns of development had the highest measure of relative importance (accounting for 53% of model's explanatory power), explaining significantly more variance in water use compared to socio‐economic or environmental variables commonly used to estimate water use. Integrating metrics characterizing the spatial pattern of development into water use models is likely to increase their utility and could facilitate water‐efficient land use planning.

North Carolina, South Carolina↗

On the sensitivity of annual streamflow to air temperature

Although interannual streamflow variability is primarily a result of precipitation variability, temperature also plays a role. The relative weakness of the temperature effect at the annual time scale hinders understanding, but may belie substantial importance on climatic time scales. Here we develop and evaluate a simple theory relating variations of streamflow and evapotranspiration ( E ) to those of precipitation ( P ) and temperature. The theory is based on extensions of the Budyko water‐balance hypothesis, the Priestley‐Taylor theory for potential evapotranspiration ( ), and a linear model of interannual basin storage. The theory implies that the temperature affects streamflow by modifying evapotranspiration through a Clausius‐Clapeyron‐like relation and through the sensitivity of net radiation to temperature. We apply and test (1) a previously introduced “strong” extension of the Budyko hypothesis, which requires that the function linking temporal variations of the evapotranspiration ratio ( E / P ) and the index of dryness ( / P ) at an annual time scale is identical to that linking interbasin variations of the corresponding long‐term means, and (2) a “weak” extension, which requires only that the annual evapotranspiration ratio depends uniquely on the annual index of dryness, and that the form of that dependence need not be known a priori nor be identical across basins. In application of the weak extension, the readily observed sensitivity of streamflow to precipitation contains crucial information about the sensitivity to potential evapotranspiration and, thence, to temperature. Implementation of the strong extension is problematic, whereas the weak extension appears to capture essential controls of the temperature effect efficiently.

Water Resources Research↗

Regional variability of nitrate fluxes in the unsaturated zone and groundwater, Wisconsin, USA

Process-based modeling of regional NO3− fluxes to groundwater is critical for understanding and managing water quality, but the complexity of NO3− reactive transport processes make implementation a challenge. This study introduces a regional vertical flux method (VFM) for efficient estimation of reactive transport of NO3− in the vadose zone and groundwater. The regional VFM was applied to 443 well samples in central-eastern Wisconsin. Chemical measurements included O2, NO3−, N2 from denitrification, and atmospheric tracers of groundwater age including carbon-14, chlorofluorocarbons, tritium, and tritiogenic helium. VFM results were consistent with observed chemistry, and calibrated parameters were in-line with estimates from previous studies. Results indicated that (1) unsaturated zone travel times were a substantial portion of the transit time to wells and streams (2) since 1945 fractions of applied N leached to groundwater have increased for manure-N, possibly due to increased injection of liquid manure, and decreased for fertilizer-N, and (3) under current practices and conditions, approximately 60% of the shallow aquifer will eventually be affected by downward migration of NO3−, with denitrification protecting the remaining 40%. Recharge variability strongly affected the unsaturated zone lag times and the eventual depth of the NO3− front. Principal components regression demonstrated that VFM parameters and predictions were significantly correlated with hydrogeochemical landscape features. The diverse and sometimes conflicting aspects of N management (e.g. limiting N volatilization versus limiting N losses to groundwater) warrant continued development of large-scale holistic strategies to manage water quality and quantity.

Wisconsin↗

The role of frozen soil in groundwater discharge predictions for warming alpine watersheds

Climate warming may alter the quantity and timing of groundwater discharge to streams in high alpine watersheds due to changes in the timing of the duration of seasonal freezing in the subsurface and snowmelt recharge. It is imperative to understand the effects of seasonal freezing and recharge on groundwater discharge to streams in warming alpine watersheds as streamflow originating from these watersheds is a critical water resource for downstream users. This study evaluates how climate warming may alter groundwater discharge due to changes in seasonally frozen ground and snowmelt using a 2‐D coupled flow and heat transport model with freeze and thaw capabilities for variably saturated media. The model is applied to a representative snowmelt‐dominated watershed in the Rocky Mountains of central Colorado, USA, with snowmelt time series reconstructed from a 12 year data set of hydrometeorological records and satellite‐derived snow covered area. Model analyses indicate that the duration of seasonal freezing in the subsurface controls groundwater discharge to streams, while snowmelt timing controls groundwater discharge to hillslope faces. Climate warming causes changes to subsurface ice content and duration, rerouting groundwater flow paths but not altering the total magnitude of future groundwater discharge outside of the bounds of hydrologic parameter uncertainties. These findings suggest that frozen soil routines play an important role for predicting the future location of groundwater discharge in watersheds underlain by seasonally frozen ground.

Colorado↗

The timing and magnitude of changes to Hortonian overland flow at the watershed scale during the post-fire recovery process

Extreme hydrologic responses following wildfires can lead to floods and debris flows with costly economic and societal impacts. Process-based hydrologic and geomorphic models used to predict the downstream impacts of wildfire must account for temporal changes in hydrologic parameters related to the generation and subsequent routing of infiltration-excess overland flow across the landscape. However, we lack quantitative relationships showing how parameters change with time-since-burning, particularly at the watershed scale. To assess variations in best-fit hydrologic parameters with time, we used the KINEROS2 hydrological model to explore temporal changes in hillslope saturated hydraulic conductivity ( K sh ) and channel hydraulic roughness ( n c ) following a wildfire in the upper Arroyo Seco watershed (41.5 km 2 ), which burned during the 2009 Station fire in the San Gabriel Mountains, California, USA. This study explored runoff-producing storms between 2008 and 2014 to infer watershed hydraulic properties by calibrating the model to observations at the watershed outlet. Modelling indicates K sh is lowest in the first year following the fire and then increases at an average rate of approximately 4.2 mm/h/year during the first 5 years of recovery. The estimated values for K sh in the first year following the fire are similar to those obtained in previous studies on smaller watersheds (<1.5 km 2 ) following the Station fire, suggesting hydrologic changes detected here can be applied to lower-order watersheds. Hydraulic roughness, n c , was lowest in the first year following the fire, but increased by a factor of 2 after 1 year of recovery. Post-fire observations suggest changes in n c are due to changes in grain roughness and vegetation in channels. These results provide quantitative constraints on the magnitude of fire-induced hydrologic changes following severe wildfires in chaparral-dominated ecosystems as well as the timing of hydrologic recovery.

California↗

Characterizing particle-scale equilibrium adsorption and kinetics of uranium(VI) desorption from U-contaminated sediments

Rates of U(VI) release from individual dry-sieved size fractions of a field-aggregated, field-contaminated composite sediment from the seasonally saturated lower vadose zone of the Hanford 300-Area were examined in flow-through reactors to maintain quasi-constant chemical conditions. The principal source of variability in equilibrium U(VI) adsorption properties of the various size fractions was the impact of variable chemistry on adsorption. This source of variability was represented using surface complexation models (SCMs) with different stoichiometric coefficients with respect to hydrogen ion and carbonate concentrations for the different size fractions. A reactive transport model incorporating equilibrium expressions for cation exchange and calcite dissolution, along with rate expressions for aerobic respiration and silica dissolution, described the temporal evolution of solute concentrations observed during the flow-through reactor experiments. Kinetic U(VI) desorption was well described using a multirate SCM with an assumed lognormal distribution for the mass-transfer rate coefficients. The estimated mean and standard deviation of the rate coefficients were the same for all <2 mm size fractions but differed for the 2–8 mm size fraction. Micropore volumes, assessed using t-plots to analyze N2 desorption data, were also the same for all dry-sieved <2 mm size fractions, indicating a link between micropore volumes and mass-transfer rate properties. Pore volumes for dry-sieved size fractions exceeded values for the corresponding wet-sieved fractions. We hypothesize that repeated field wetting and drying cycles lead to the formation of aggregates and/or coatings containing (micro)pore networks which provided an additional mass-transfer resistance over that associated with individual particles. The 2–8 mm fraction exhibited a larger average and standard deviation in the distribution of mass-transfer rate coefficients, possibly caused by the abundance of microporous basaltic rock fragments.

Water Resources Research↗

Seawater circulation in sediments driven by interactions between seabed topography and fluid density

Measurements of submarine groundwater discharge (SGD) in coastal areas often show that the saltwater discharge component is substantially greater than the freshwater discharge. Several mechanisms have been proposed to explain these high saltwater discharge values, including saltwater circulation driven by wave and tidal pumping, wave and tidal setup in intertidal areas, currents over bedforms, and density gradients resulting from mixing along the freshwater-saltwater interface. In this study, a new mechanism for saltwater circulation and discharge is proposed and evaluated. The process results from interaction between bedform topography and buoyancy forces, even without flow or current over the bedform. In this mechanism, an inverted salinity (and density) profile in the presence of both a bedform on the seafloor and an upward flow of fresher groundwater from depth induces a downward flow of saline pore water under the troughs and upward flow under the adjacent crest of the bedform. The magnitude and occurrence of the mechanism were tested using numerical methods. The results indicate that this mechanism could drive seawater circulation under a limited range of conditions and contribute 20%–30% of local SGD when and where the process is operative. Bedform shape, hydraulic conductivity, hydraulic head, and salinity at depth in the porous media, aquifer thickness, effective porosity, and hydrodynamic dispersion are among the factors that control the occurrence and magnitude of the circulation of seawater by this mechanism.

Water Resources Research↗

Balancing practicality and hydrologic realism: a parsimonious approach for simulating rapid groundwater recharge via unsaturated-zone preferential flow

The impact of preferential flow on recharge and contaminant transport poses a considerable challenge to water-resources management. Typical hydrologic models require extensive site characterization, but can underestimate fluxes when preferential flow is significant. A recently developed source-responsive model incorporates film-flow theory with conservation of mass to estimate unsaturated-zone preferential fluxes with readily available data. The term source-responsive describes the sensitivity of preferential flow in response to water availability at the source of input. We present the first rigorous tests of a parsimonious formulation for simulating water table fluctuations using two case studies, both in arid regions with thick unsaturated zones of fractured volcanic rock. Diffuse flow theory cannot adequately capture the observed water table responses at both sites; the source-responsive model is a viable alternative. We treat the active area fraction of preferential flow paths as a scaled function of water inputs at the land surface then calibrate the macropore density to fit observed water table rises. Unlike previous applications, we allow the characteristic film-flow velocity to vary, reflecting the lag time between source and deep water table responses. Analysis of model performance and parameter sensitivity for the two case studies underscores the importance of identifying thresholds for initiation of film flow in unsaturated rocks, and suggests that this parsimonious approach is potentially of great practical value.

Water Resources Research↗

Variations in surface water-ground water interactions along a headwater mountain stream : comparisons between transient storage and water balance analyses

The accumulation of discharge along a stream valley is frequently assumed to be the primary control on solute transport processes. Relationships of both increasing and decreasing transient storage, and decreased gross losses of stream water have been reported with increasing discharge; however, we have yet to validate these relationships with extensive field study. We conducted transient storage and mass recovery analyses of artificial tracer studies completed for 28 contiguous 100 m reaches along a stream valley, repeated under four base-flow conditions. We calculated net and gross gains and losses, temporal moments of tracer breakthrough curves, and best fit transient storage model parameters (with uncertainty estimates) for 106 individual tracer injections. Results supported predictions that gross loss of channel water would decrease with increased discharge. However, results showed no clear relationship between discharge and transient storage, and further analysis of solute tracer methods demonstrated that the lack of this relation may be explained by uncertainty and equifinality in the transient storage model framework. Furthermore, comparison of water balance and transient storage approaches reveals complications in clear interpretation of either method due to changes in advective transport time, which sets a the temporal boundary separating transient storage and channel water balance. We have little ability to parse this limitation of solute tracer methods from the physical processes we seek to study. We suggest the combined analysis of both transient storage and channel water balance more completely characterizes transport of solutes in stream networks than can be inferred from either method alone.

Montana↗

Colloid transport in saturated porous media: Elimination of attachment efficiency in a new colloid transport model

A colloid transport model is introduced that is conceptually simple yet captures the essential features of colloid transport and retention in saturated porous media when colloid retention is dominated by the secondary minimum because an electrostatic barrier inhibits substantial deposition in the primary minimum. This model is based on conventional colloid filtration theory (CFT) but eliminates the empirical concept of attachment efficiency. The colloid deposition rate is computed directly from CFT by assuming all predicted interceptions of colloids by collectors result in at least temporary deposition in the secondary minimum. Also, a new paradigm for colloid re-entrainment based on colloid population heterogeneity is introduced. To accomplish this, the initial colloid population is divided into two fractions. One fraction, by virtue of physiochemical characteristics (e.g., size and charge), will always be re-entrained after capture in a secondary minimum. The remaining fraction of colloids, again as a result of physiochemical characteristics, will be retained “irreversibly” when captured by a secondary minimum. Assuming the dispersion coefficient can be estimated from tracer behavior, this model has only two fitting parameters: (1) the fraction of the initial colloid population that will be retained “irreversibly” upon interception by a secondary minimum, and (2) the rate at which reversibly retained colloids leave the secondary minimum. These two parameters were correlated to the depth of the Derjaguin-Landau-Verwey-Overbeek (DLVO) secondary energy minimum and pore-water velocity, two physical forces that influence colloid transport. Given this correlation, the model serves as a heuristic tool for exploring the influence of physical parameters such as surface potential and fluid velocity on colloid transport.

Water Resources Research↗

Temporal variability of exchange between groundwater and surface water based on high-frequency direct measurements of seepage at the sediment-water interface

Seepage at the sediment-water interface in several lakes, a large river, and an estuary exhibits substantial temporal variability when measured with temporal resolution of 1 min or less. Already substantial seepage rates changed by 7% and 16% in response to relatively small rain events at two lakes in the northeastern USA, but did not change in response to two larger rain events at a lake in Minnesota. However, seepage at that same Minnesota lake changed by 10% each day in response to withdrawals from evapotranspiration. Seepage increased by more than an order of magnitude when a seiche occurred in the Great Salt Lake, Utah. Near the head of a fjord in Puget Sound, Washington, seepage in the intertidal zone varied greatly from −115 to +217 cm d−1 in response to advancing and retreating tides when the time-averaged seepage was upward at +43 cm d−1. At all locations, seepage variability increased by one to several orders of magnitude in response to wind and associated waves. Net seepage remained unchanged by wind unless wind also induced a lake seiche. These examples from sites distributed across a broad geographic region indicate that temporal variability in seepage in response to common hydrological events is much larger than previously realized. At most locations, seepage responded within minutes to changes in surface-water stage and within minutes to hours to groundwater recharge associated with rainfall. Likely implications of this dynamism include effects on water residence time, geochemical transformations, and ecological conditions at and near the sediment-water interface.

Water Resources Research↗

How runoff begins (and ends): characterizing hydrologic response at the catchment scale

Improved understanding of the complex dynamics associated with spatially and temporally variable runoff response is needed to better understand the hydrology component of interdisciplinary problems. The objective of this study was to quantitatively characterize the environmental controls on runoff generation for the range of different streamflow-generation mechanisms illustrated in the classic Dunne diagram. The comprehensive physics-based model of coupled surface-subsurface flow, InHM, is employed in a heuristic mode. InHM has been employed previously to successfully simulate the observed hydrologic response at four diverse, well-characterized catchments, which provides the foundation for this study. The C3 and CB catchments are located within steep, forested terrain; the TW and R5 catchments are located in gently sloping rangeland. The InHM boundary-value problems for these four catchments provide the corner-stones for alternative simulation scenarios designed to address the question of how runoff begins (and ends). Simulated rainfall-runoff events are used to systematically explore the impact of soil-hydraulic properties and rainfall characteristics. This approach facilitates quantitative analysis of both integrated and distributed hydrologic responses at high-spatial and temporal resolution over the wide range of environmental conditions represented by the four catchments. The results from 140 unique simulation scenarios illustrate how rainfall intensity/depth, subsurface permeability contrasts, characteristic curve shapes, and topography provide important controls on the hydrologic-response dynamics. The processes by which runoff begins (and ends) are shown, in large part, to be defined by the relative rates of rainfall, infiltration, lateral flow convergence, and storage dynamics within the variably saturated soil layers.

Water Resources Research↗

Coupled hydrogeomorphic and woody-seedling responses to controlled flood releases in a dryland river

Interactions among flow, geomorphic processes, and riparian vegetation can strongly influence both channel form and vegetation communities. To investigate such interactions, we took advantage of a series of dam-managed flood releases that were designed in part to maintain a native riparian woodland system on a sand-bed, dryland river, the Bill Williams River, Arizona, USA. Our resulting multiyear flow experiment examined differential mortality among native and nonnative riparian seedlings, associated flood hydraulics and geomorphic changes, and the temporal evolution of feedbacks among vegetation, channel form, and hydraulics. We found that floods produced geomorphic and vegetation responses that varied with distance downstream of a dam, with scour and associated seedling mortality closer to the dam and aggradation and burial-induced mortality in a downstream reach. We also observed significantly greater mortality among nonnative tamarisk (Tamarix) seedlings than among native willow (Salix gooddingii) seedlings, reflecting the greater first-year growth of willow relative to tamarisk. When vegetation was small early in our study period, the effects of vegetation on flood hydraulics and on mediating flood-induced channel change were minimal. Vegetation growth in subsequent years resulted in stronger feedbacks, such that vegetation's stabilizing effect on bars and its drag effect on flow progressively increased, muting the geomorphic effects of a larger flood release. These observations suggest that the effectiveness of floods in producing geomorphic and ecological changes varies not only as a function of flood magnitude and duration, but also of antecedent vegetation density and size.

Arizona↗

Using heat as a tracer to estimate spatially distributed mean residence times in the hyporheic zone of a riffle-pool sequence

Biochemical reactions that occur in the hyporheic zone are highly dependent on the time solutes that are in contact with sediments of the riverbed. In this investigation, we developed a 2-D longitudinal flow and solute-transport model to estimate the spatial distribution of mean residence time in the hyporheic zone. The flow model was calibrated using observations of temperature and pressure, and the mean residence times were simulated using the age-mass approach for steady-state flow conditions. The approach used in this investigation includes the mixing of different ages and flow paths of water through advection and dispersion. Uncertainty of flow and transport parameters was evaluated using standard Monte Carlo and the generalized likelihood uncertainty estimation method. Results of parameter estimation support the presence of a low-permeable zone in the riffle area that induced horizontal flow at a shallow depth within the riffle area. This establishes shallow and localized flow paths and limits deep vertical exchange. For the optimal model, mean residence times were found to be relatively long (9–40.0 days). The uncertainty of hydraulic conductivity resulted in a mean interquartile range (IQR) of 13 days across all piezometers and was reduced by 24% with the inclusion of temperature and pressure observations. To a lesser extent, uncertainty in streambed porosity and dispersivity resulted in a mean IQR of 2.2 and 4.7 days, respectively. Alternative conceptual models demonstrate the importance of accounting for the spatial distribution of hydraulic conductivity in simulating mean residence times in a riffle-pool sequence.

Water Resources Research↗

Impacts of an ethanol-blended fuel release on groundwater and fate of produced methane: simulation of field observations

In a field experiment at Vandenberg Air Force Base (VAFB) designed to mimic the impact of a small-volume release of E10 (10% ethanol and 90% conventional gasoline), two plumes were created by injecting extracted groundwater spiked with benzene, toluene, and o-xylene, abbreviated BToX (No-Ethanol Lane) and BToX plus ethanol (With-Ethanol Lane) for 283 days. We developed a reactive transport model to understand processes controlling the fate of ethanol and BToX. The model was calibrated to the extensive field dataset and accounted for concentrations of sulfate, iron, acetate, and methane along with iron-reducing bacteria, sulfate-reducing bacteria, fermentative bacteria, and methanogenic archaea. The benzene plume was about 4.5 times longer in the With-Ethanol Lane than in the No-Ethanol Lane. Matching this different behavior in the two lanes required inhibiting benzene degradation in the presence of ethanol. Inclusion of iron reduction with negligible growth of iron-reducers was required to reproduce the observed constant degradation rate of benzene. Modeling suggested that vertical dispersion and diffusion of sulfate from an adjacent aquitard were important sources of sulfate in the aquifer. Matching of methane data required incorporating initial fermentation of ethanol to acetate, methane loss by outgassing, and methane oxidation coupled to sulfate and iron reduction. Simulation of microbial growth using dual Monod kinetics, and including inhibition by more favorable electron acceptors, generally resulted in reasonable yields for microbial growth of 0.01-0.05.

California↗

A generalized Grubbs-Beck test statistic for detecting multiple potentially influential low outliers in flood series

he Grubbs-Beck test is recommended by the federal guidelines for detection of low outliers in flood flow frequency computation in the United States. This paper presents a generalization of the Grubbs-Beck test for normal data (similar to the Rosner (1983) test; see also Spencer and McCuen (1996)) that can provide a consistent standard for identifying multiple potentially influential low flows. In cases where low outliers have been identified, they can be represented as “less-than” values, and a frequency distribution can be developed using censored-data statistical techniques, such as the Expected Moments Algorithm. This approach can improve the fit of the right-hand tail of a frequency distribution and provide protection from lack-of-fit due to unimportant but potentially influential low flows (PILFs) in a flood series, thus making the flood frequency analysis procedure more robust.

Water Resources Research↗

Hysteresis in suspended sediment to turbidity relations due to changing particle size distributions

Turbidity (T) is the most ubiquitous of surrogate technologies used to estimate suspended-sediment concentration (SSC). The effects of sediment size on turbidity are well documented; however, effects from changes in particle size distributions (PSD) are rarely evaluated. Hysteresis in relations of SSC-to-turbidity (SSC~T) for single stormflow events was observed and quantified for a data set of 195 concurrent measurements of SSC, turbidity, discharge, velocity, and volumetric PSD collected during five stormflows in 2009–2010 on Yellow River at Gees Mill Road in metropolitan Atlanta, Georgia. Regressions of SSC-normalized turbidity (T/SSC) on concurrently measured PSD percentiles show an inverse, exponential influence of particle size on turbidity that is not constant across the size range of the PSD. The majority of the influence of PSD on T/SSC is from particles of fine-silt and smaller sizes (finer than 16 microns). This study shows that small changes in the often assumed stability of the PSD are significant to SSC~T relations. Changes of only 5 microns in the fine silt and smaller size fractions of suspended sediment PSD can produce hysteresis in the SSC~T rating that can increase error and produce bias. Observed SSC~T hysteresis may be an indicator of changes in sediment properties during stormflows and of potential changes in sediment sources. Trends in the PSD time series indicate that sediment transport is capacity-limited for sand-sized sediment in the channel and supply-limited for fine silt and smaller sediment from the hillslope.

Georgia↗