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 955 records · Page 53Linked to original sources
USGS2010

Factors influencing wood mobilization in Minnesota streams

Natural pieces of wood provide a variety of ecosystem functions in streams including habitat, organic matter retention, increased hyporheic exchange and transient storage, and enhanced hydraulic and geomorphic heterogeneity. Wood mobilization is a critical process in determining the residence time of wood. We documented the characteristics and locations of 865 natural wood pieces (>0.05 m in diameter for a portion >1 m in length) in nine streams along the north shore of Lake Superior in Minnesota. We determined the locations of the pieces again after an overbank stormflow event to determine the factors that influenced mobilization of stationary wood pieces in natural streams. Seven of 11 potential predictor variables were identified with multiple logistic regression as significant to mobilization: burial, effective depth, ratio of piece length to effective stream width (length ratio), bracing, rootwad presence, downstream force ratio, and draft ratio. The final model ( P < 0.001, r 2 = 0.39) indicated that wood mobilization under natural conditions is a complex function of both mechanical factors (burial, length ratio, bracing, rootwad presence, draft ratio) and hydraulic factors (effective depth, downstream force ratio). If stable pieces are a goal for stream management then features such as partial burial, low effective depth, high length relative to channel width, bracing against other objects (e.g., stream banks, trees, rocks, or larger wood pieces), and rootwads are desirable. Using the model equation from this study, stewards of natural resources can better manage in-stream wood for the benefit of stream ecosystems.

Water Resources Research↗
USGS2010

Hillslope hydrologic connectivity controls riparian groundwater turnover: Implications of catchment structure for riparian buffering and stream water sources

Hydrologic connectivity between catchment upland and near stream areas is essential for the transmission of water, solutes, and nutrients to streams. However, our current understanding of the role of riparian zones in mediating landscape hydrologic connectivity and the catchment scale export of water and solutes is limited. We tested the relationship between the duration of hillslope‐riparian‐stream (HRS) hydrologic connectivity and the rate and degree of riparian shallow groundwater turnover along four HRS well transects within a set of nested mountain catchments (Tenderfoot Creek Experimental Forest, MT). Transect HRS water table connectivity ranged from 9 to 123 days during the annual snowmelt hydrograph. Hillslope water was always characterized by low specific conductance (∼27 μ S cm −1 ). In transects with transient hillslope water tables, riparian groundwater specific conductance was elevated during base flow conditions (∼127 μ S cm −1 ) but shifted toward hillslope signatures once a HRS groundwater connection was established. The degree of riparian groundwater turnover was proportional to the duration of HRS connectivity and inversely related to the riparian: hillslope area ratios (buffer ratio; r 2 = 0.95). We applied this relationship to the stream network in seven subcatchments within the Tenderfoot Creek Experimental Forest and compared their turnover distributions to source water contributions measured at each catchment outlet. The amount of riparian groundwater exiting each of the seven catchments was linearly related ( r 2 = 0.92) to their median riparian turnover time. Our observations suggest that the size and spatial arrangement of hillslope and riparian zones along a stream network and the timing and duration of groundwater connectivity between them is a first‐order control on the magnitude and timing of water and solutes observed at the catchment outlet.

Water Resources Research↗
USGS2010

Mixing effects on apparent reaction rates and isotope fractionation during denitrification in a heterogeneous aquifer

Gradients in contaminant concentrations and isotopic compositions commonly are used to derive reaction parameters for natural attenuation in aquifers. Differences between field‐scale (apparent) estimated reaction rates and isotopic fractionations and local‐scale (intrinsic) effects are poorly understood for complex natural systems. For a heterogeneous alluvial fan aquifer, numerical models and field observations were used to study the effects of physical heterogeneity on reaction parameter estimates. Field measurements included major ions, age tracers, stable isotopes, and dissolved gases. Parameters were estimated for the O 2 reduction rate, denitrification rate, O 2 threshold for denitrification, and stable N isotope fractionation during denitrification. For multiple geostatistical realizations of the aquifer, inverse modeling was used to establish reactive transport simulations that were consistent with field observations and served as a basis for numerical experiments to compare sample‐based estimates of “apparent” parameters with “true“ (intrinsic) values. For this aquifer, non‐Gaussian dispersion reduced the magnitudes of apparent reaction rates and isotope fractionations to a greater extent than Gaussian mixing alone. Apparent and true rate constants and fractionation parameters can differ by an order of magnitude or more, especially for samples subject to slow transport, long travel times, or rapid reactions. The effect of mixing on apparent N isotope fractionation potentially explains differences between previous laboratory and field estimates. Similarly, predicted effects on apparent O 2 threshold values for denitrification are consistent with previous reports of higher values in aquifers than in the laboratory. These results show that hydrogeological complexity substantially influences the interpretation and prediction of reactive transport.

Water Resources Research↗
USGS2011

Spatial connectivity in a highly heterogeneous aquifer: From cores to preferential flow paths

This study investigates connectivity in a small portion of the extremely heterogeneous aquifer at the Macrodispersion Experiment (MADE) site in Columbus, Mississippi. A total of 19 fully penetrating soil cores were collected from a rectangular grid of 4 m by 4 m. Detailed grain size analysis was performed on 5 cm segments of each core, yielding 1740 hydraulic conductivity ( K ) estimates. Three different geostatistical simulation methods were used to generate 3‐D conditional realizations of the K field for the sampled block. Particle tracking calculations showed that the fastest particles, as represented by the first 5% to arrive, converge along preferential flow paths and exit the model domain within preferred areas. These 5% fastest flow paths accounted for about 40% of the flow. The distribution of preferential flow paths and particle exit locations is clearly influenced by the occurrence of clusters formed by interconnected cells with K equal to or greater than the 0.9 decile of the data distribution (10% of the volume). The fraction of particle paths within the high‐ K clusters ranges from 43% to 69%. In variogram‐based K fields, some of the fastest paths are through media with lower K values, suggesting that transport connectivity may not require fully connected zones of relatively homogenous K . The high degree of flow and transport connectivity was confirmed by the values of two groups of connectivity indicators. In particular, the ratio between effective and geometric mean K (on average, about 2) and the ratio between the average arrival time and the arrival time of the fastest particles (on average, about 9) are consistent with flow and advective transport behavior characterized by channeling along preferential flow paths.

Water Resources Research↗
USGS2011

Historical summer base flow and stormflow trends for New England rivers

River base flow is important to aquatic ecosystems, particularly because of its influence on summer water temperatures. Summer (June through September) daily mean streamflows were separated into base flow and stormflow components by use of an automated method at 25 stations in the New England region of the United States that drain predominantly natural basins. Summer monthly mean base flows increased from 1950–2006 at most stations in western New England with many large increases (>20%) and some very large increases (>50%) in and near New Hampshire and Vermont. The same was true for increases in summer 7 day low base flows in and near New Hampshire and Vermont during this same period; in contrast, there were small and large decreases in 7 day low base flows in northern and coastal areas of Maine. Summer stormflows increased from 1950–2006 by more than 50% at many stations in New England, particularly in and near New Hampshire and Vermont. The increases in base flows and stormflows at many stations in and near New Hampshire and Vermont were likely driven by the large increases in summer precipitation recorded at weather stations in this area. Summer rainfall increased at most weather stations in New England from 1950–2006 with many increases of more than 20% in western New England. Summer air temperature increased on average by 1.1°C from 1950–2006 in New England and may have played a role in the decreased base flows in northern and coastal Maine through increased evapotranspiration. Many variables increased less from 1930–2006 than from 1950–2006.

Connecticut, Maine, Massachusetts, New Hampshire, ↗
USGS2010

Use of electrical imaging and distributed temperature sensing methods to characterize surface water–groundwater exchange regulating uranium transport at the Hanford 300 Area, Washington

We explored the use of continuous waterborne electrical imaging (CWEI), in conjunction with fiber‐optic distributed temperature sensor (FO‐DTS) monitoring, to improve the conceptual model for uranium transport within the Columbia River corridor at the Hanford 300 Area, Washington. We first inverted resistivity and induced polarization CWEI data sets for distributions of electrical resistivity and polarizability, from which the spatial complexity of the primary hydrogeologic units was reconstructed. Variations in the depth to the interface between the overlying coarse‐grained, high‐permeability Hanford Formation and the underlying finer‐grained, less permeable Ringold Formation, an important contact that limits vertical migration of contaminants, were resolved along ∼3 km of the river corridor centered on the 300 Area. Polarizability images were translated into lithologic images using established relationships between polarizability and surface area normalized to pore volume ( S por ). The FO‐DTS data recorded along 1.5 km of cable with a 1 m spatial resolution and 5 min sampling interval revealed subreaches showing (1) temperature anomalies (relatively warm in winter and cool in summer) and (2) a strong correlation between temperature and river stage (negative in winter and positive in summer), both indicative of reaches of enhanced surface water–groundwater exchange. The FO‐DTS data sets confirm the hydrologic significance of the variability identified in the CWEI and reveal a pattern of highly focused exchange, concentrated at springs where the Hanford Formation is thickest. Our findings illustrate how the combination of CWEI and FO‐DTS technologies can characterize surface water–groundwater exchange in a complex, coupled river‐aquifer system.

Washington↗
USGS2011

Predictive uncertainty analysis of a saltwater intrusion model using null-space Monte Carlo

Because of the extensive computational burden and perhaps a lack of awareness of existing methods, rigorous uncertainty analyses are rarely conducted for variable-density flow and transport models. For this reason, a recently developed null-space Monte Carlo (NSMC) method for quantifying prediction uncertainty was tested for a synthetic saltwater intrusion model patterned after the Henry problem. Saltwater intrusion caused by a reduction in fresh groundwater discharge was simulated for 1000 randomly generated hydraulic conductivity distributions, representing a mildly heterogeneous aquifer. From these 1000 simulations, the hydraulic conductivity distribution giving rise to the most extreme case of saltwater intrusion was selected and was assumed to represent the "true" system. Head and salinity values from this true model were then extracted and used as observations for subsequent model calibration. Random noise was added to the observations to approximate realistic field conditions. The NSMC method was used to calculate 1000 calibration-constrained parameter fields. If the dimensionality of the solution space was set appropriately, the estimated uncertainty range from the NSMC analysis encompassed the truth. Several variants of the method were implemented to investigate their effect on the efficiency of the NSMC method. Reducing the dimensionality of the null-space for the processing of the random parameter sets did not result in any significant gains in efficiency and compromised the ability of the NSMC method to encompass the true prediction value. The addition of intrapilot point heterogeneity to the NSMC process was also tested. According to a variogram comparison, this provided the same scale of heterogeneity that was used to generate the truth. However, incorporation of intrapilot point variability did not make a noticeable difference to the uncertainty of the prediction. With this higher level of heterogeneity, however, the computational burden of generating calibration-constrained parameter fields approximately doubled. Predictive uncertainty variance computed through the NSMC method was compared with that computed through linear analysis. The results were in good agreement, with the NSMC method estimate showing a slightly smaller range of prediction uncertainty than was calculated by the linear method. Copyright 2011 by the American Geophysical Union.

Water Resources Research↗
USGS2011

Transient groundwater chemistry near a river: Effects on U(VI) transport in laboratory column experiments

In the 300 Area of a U(VI)-contaminated aquifer at Hanford, Washington, USA, inorganic carbon and major cations, which have large impacts on U(VI) transport, change on an hourly and seasonal basis near the Columbia River. Batch and column experiments were conducted to investigate the factors controlling U(VI) adsorption/desorption by changing chemical conditions over time. Low alkalinity and low Ca concentrations (Columbia River water) enhanced adsorption and reduced aqueous concentrations. Conversely, high alkalinity and high Ca concentrations (Hanford groundwater) reduced adsorption and increased aqueous concentrations of U(VI). An equilibrium surface complexation model calibrated using laboratory batch experiments accounted for the decrease in U(VI) adsorption observed with increasing (bi)carbonate concentrations and other aqueous chemical conditions. In the column experiment, alternating pulses of river and groundwater caused swings in aqueous U(VI) concentration. A multispecies multirate surface complexation reactive transport model simulated most of the major U(VI) changes in two column experiments. The modeling results also indicated that U(VI) transport in the studied sediment could be simulated by using a single kinetic rate without loss of accuracy in the simulations. Moreover, the capability of the model to predict U(VI) transport in Hanford groundwater under transient chemical conditions depends significantly on the knowledge of real-time change of local groundwater chemistry.

Washington↗
USGS2010

Influence of hummocks and emergent vegetation on hydraulic performance in a surface flow wastewater treatment wetland

A series of tracer experiments were conducted biannually at the start and end of the vegetation growing season in a surface flow wastewater treatment wetland located near Phoenix, AZ. Tracer experiments were conducted prior to and following reconfiguration and replanting of a 1.2 ha treatment wetland from its original design of alternating shallow and deep zones to incorporate hummocks (shallow planting beds situated perpendicular to flow). Tracer test data were analyzed using analysis of moments and the one‐dimensional transport with inflow and storage numerical model to evaluate the effects of the seasonal vegetation growth cycle and hummocks on solute transport. Following reconfiguration, vegetation coverage was relatively small, and minor changes in spatial distribution influenced wetland hydraulics. During start‐up conditions, the wetland underwent an acclimation period characterized by small vegetation coverage and large transport cross‐sectional areas. At the start of the growing season, new growth of emergent vegetation enhanced hydraulic performance. At the end of the growing season, senescing vegetation created short‐circuiting. Wetland hydrodynamics were associated with high volumetric efficiencies and velocity heterogeneities. The hummock design resulted in breakthrough curves characterized by multiple secondary tracer peaks indicative of varied flow paths created by bottom topography.

Arizona↗
USGS2011

Assessing the detail needed to capture rainfall-runoff dynamics with physics-based hydrologic response simulation

Concept development simulation with distributed, physics-based models provides a quantitative approach for investigating runoff generation processes across environmental conditions. Disparities within data sets employed to design and parameterize boundary value problems used in heuristic simulation inevitably introduce various levels of bias. The objective was to evaluate the impact of boundary value problem complexity on process representation for different runoff generation mechanisms. The comprehensive physics-based hydrologic response model InHM has been employed to generate base case simulations for four well-characterized catchments. The C3 and CB catchments are located within steep, forested environments dominated by subsurface stormflow; the TW and R5 catchments are located in gently sloping rangeland environments dominated by Dunne and Horton overland flows. Observational details are well captured within all four of the base case simulations, but the characterization of soil depth, permeability, rainfall intensity, and evapotranspiration differs for each. These differences are investigated through the conversion of each base case into a reduced case scenario, all sharing the same level of complexity. Evaluation of how individual boundary value problem characteristics impact simulated runoff generation processes is facilitated by quantitative analysis of integrated and distributed responses at high spatial and temporal resolution. Generally, the base case reduction causes moderate changes in discharge and runoff patterns, with the dominant process remaining unchanged. Moderate differences between the base and reduced cases highlight the importance of detailed field observations for parameterizing and evaluating physics-based models. Overall, similarities between the base and reduced cases indicate that the simpler boundary value problems may be useful for concept development simulation to investigate fundamental controls on the spectrum of runoff generation mechanisms. Copyright 2011 by the American Geophysical Union.

Water Resources Research↗
USGS2011

Rethinking hyporheic flow and transient storage to advance understanding of stream-catchment connections

Although surface water and groundwater are increasingly referred to as one resource, there remain environmental and ecosystem needs to study the 10 m to 1 km reach scale as one hydrologic system. Streams gain and lose water over a range of spatial and temporal scales. Large spatial scales (kilometers) have traditionally been recognized and studied as river-aquifer connections. Over the last 25 years hyporheic exchange flows (1-10 m) have been studied extensively. Often a transient storage model has been used to quantify the physical solute transport setting in which biogeochemical processes occur. At the longer 10 m to 1 km scale of stream reaches it is now clear that streams which gain water overall can coincidentally lose water to the subsurface. At this scale, the amounts of water transferred are not necessarily significant but the exchanges can, however, influence solute transport. The interpretation of seemingly straightforward questions about water, contaminant, and nutrient fluxes into and along a stream can be confounded by flow losses which are too small to be apparent in stream gauging and along flow paths too long to be detected in tracer experiments. We suggest basic hydrologic approaches, e.g., measurement of flow along the channel, surface and subsurface solute sampling, and routine measurements of the water table that, in our opinion, can be used to extend simple exchange concepts from the hyporheic exchange scale to a scale of stream-catchment connection.

Water Resources Research↗
USGS2011

Hydrologic response of catchments to precipitation: Quantification of mechanical carriers and origins of water

Precipitation-induced overland and groundwater flow and mixing processes are quantified to analyze the temporal (event and pre-event water) and spatial (groundwater discharge and overland runoff) origins of water entering a stream. Using a distributed-parameter control volume finite-element simulator that can simultaneously solve the fully coupled partial differential equations describing 2-D Manning and 3-D Darcian flow and advective-dispersive transport, mechanical flow (driven by hydraulic potential) and tracer-based hydrograph separation (driven by dispersive mixing as well as mechanical flow) are simulated in response to precipitation events in two cross sections oriented parallel and perpendicular to a stream. The results indicate that as precipitation becomes more intense, the subsurface mechanical flow contributions tend to become less significant relative to the total pre-event stream discharge. Hydrodynamic mixing can play an important role in enhancing pre-event tracer signals in the stream. This implies that temporally tagged chemical signals introduced into surface-subsurface flow systems from precipitation may not be strong enough to detect the changes in the subsurface flow system. It is concluded that diffusive/dispersive mixing, capillary fringe groundwater ridging, and macropore flow can influence the temporal sources of water in the stream, but any sole mechanism may not fully explain the strong pre-event water discharge. Further investigations of the influence of heterogeneity, residence time, geomorphology, and root zone processes are required to confirm the conclusions of this study. Copyright 2011 by the American Geophysical Union.

Water Resources Research↗
USGS2011

Simulating adsorption of U(VI) under transient groundwater flow and hydrochemistry: Physical versus chemical nonequilibrium model

Coupled intragrain diffusional mass transfer and nonlinear surface complexation processes play an important role in the transport behavior of U(VI) in contaminated aquifers. Two alternative model approaches for simulating these coupled processes were analyzed and compared: (1) the physical nonequilibrium approach that explicitly accounts for aqueous speciation and instantaneous surface complexation reactions in the intragrain regions and approximates the diffusive mass exchange between the immobile intragrain pore water and the advective pore water as multirate first-order mass transfer and (2) the chemical nonequilibrium approach that approximates the diffusion-limited intragrain surface complexation reactions by a set of multiple first-order surface complexation reaction kinetics, thereby eliminating the explicit treatment of aqueous speciation in the intragrain pore water. A model comparison has been carried out for column and field scale scenarios, representing the highly transient hydrological and geochemical conditions in the U(VI)-contaminated aquifer at the Hanford 300A site, Washington, USA. It was found that the response of U(VI) mass transfer behavior to hydrogeochemically induced changes in U(VI) adsorption strength was more pronounced in the physical than in the chemical nonequilibrium model. The magnitude of the differences in model behavior depended particularly on the degree of disequilibrium between the advective and immobile phase U(VI) concentrations. While a clear difference in U(VI) transport behavior between the two models was noticeable for the column-scale scenarios, only minor differences were found for the Hanford 300A field scale scenarios, where the model-generated disequilibrium conditions were less pronounced as a result of frequent groundwater flow reversals.

Water Resources Research↗
USGS2011

The significance of turbulent flow representation in single-continuum models

Karst aquifers exhibit highly conductive features caused from rock dissolution processes. Flow within these structures can become turbulent and therefore can be expressed by nonlinear gradient functions. One way to account for these effects is by coupling a continuum model with a conduit network. Alternatively, turbulent flow can be considered by adapting the hydraulic conductivity within the continuum model. Consequently, the significance of turbulent flow on the dynamic behavior of karst springs is investigated by an enhanced single-continuum model that results in conduit-type flow in continuum cells (CTFC). The single-continuum approach CTFC represents laminar and turbulent flow as well as more complex hybrid models that require additional programming and numerical efforts. A parameter study is conducted to investigate the effects of turbulent flow on the response of karst springs to recharge events using the new CTFC approach, existing hybrid models, and MODFLOW-2005. Results reflect the importance of representing (1) turbulent flow in karst conduits and (2) the exchange between conduits and continuum cells. More specifically, laminar models overestimate maximum spring discharge and underestimate hydraulic gradients within the conduit. It follows that aquifer properties inferred from spring hydrographs are potentially impaired by ignoring flow effects due to turbulence. The exchange factor used for hybrid models is necessary to account for the scale dependency between hydraulic properties of the matrix continuum and conduits. This functionality, which is not included in CTFC, can be mimicked by appropriate use of the Horizontal Flow Barrier package for MODFLOW. Copyright 2011 by the American Geophysical Union.

Water Resources Research↗
USGS2011

Characterization of the intragranular water regime within subsurface sediments: pore volume, surface area, and mass transfer limitations

Although "intragranular" pore space within grain aggregates, grain fractures, and mineral surface coatings may contain a relatively small fraction of the total porosity within a porous medium, it often contains a significant fraction of the reactive surface area, and can thus strongly affect the transport of sorbing solutes. In this work, we demonstrate a batch experiment procedure using tritiated water as a high-resolution diffusive tracer to characterize the intragranular pore space. The method was tested using uranium-contaminated sediments from the vadose and capillary fringe zones beneath the former 300A process ponds at the Hanford site (Washington). Sediments were contacted with tracers in artificial groundwater, followed by a replacement of bulk solution with tracer-free groundwater and the monitoring of tracer release. From these data, intragranular pore volumes were calculated and mass transfer rates were quantified using a multirate first-order mass transfer model. Tritium-hydrogen exchange on surface hydroxyls was accounted for by conducting additional tracer experiments on sediment that was vacuum dried after reaction. The complementary ("wet" and "dry") techniques allowed for the simultaneous determination of intragranular porosity and surface area using tritium. The Hanford 300A samples exhibited intragranular pore volumes of ~1% of the solid volume and intragranular surface areas of ~20%&ndash;35% of the total surface area. Analogous experiments using bromide ion as a tracer yielded very different results, suggesting very little penetration of bromide into the intragranular porosity.

Water Resources Research↗
USGS2011

How restructuring river connectivity changes freshwater fish biodiversity and biogeography

Interbasin water transfer projects, in which river connectivity is restructured via man-made canals, are an increasingly popular solution to address the spatial mismatch between supply and demand of fresh water. However, the ecological consequences of such restructuring remain largely unexplored, and there are no general theoretical guidelines from which to derive these expectations. River systems provide excellent opportunities to explore how network connectivity shapes habitat occupancy, community dynamics, and biogeographic patterns. We apply a neutral model (which assumes competitive equivalence among species within a stochastic framework) to an empirically derived river network to explore how proposed changes in network connectivity may impact patterns of freshwater fish biodiversity. Without predicting the responses of individual extant species, we find the addition of canals connecting hydrologically isolated river basins facilitates the spread of common species and increases average local species richness without changing the total species richness of the system. These impacts are sensitive to the parameters controlling the spatial scale of fish dispersal, with increased dispersal affording more opportunities for biotic restructuring at the community and landscape scales. Connections between isolated basins have a much larger effect on local species richness than those connecting reaches within a river basin, even when those within-basin reaches are far apart. As a result, interbasin canal projects have the potential for long-term impacts to continental-scale riverine communities.

Water Resources Research↗
USGS2011

Effects of dynamically variable saturation and matrix-conduit coupling of flow in karst aquifers

Well-developed karst aquifers consist of highly conductive conduits and a relatively low permeability fractured and/or porous rock matrix and therefore behave as a dual-hydraulic system. Groundwater flow within highly permeable strata is rapid and transient and depends on local flow conditions, i.e., pressurized or nonpressurized flow. The characterization of karst aquifers is a necessary and challenging task because information about hydraulic and spatial conduit properties is poorly defined or unknown. To investigate karst aquifers, hydraulic stresses such as large recharge events can be simulated with hybrid (coupled discrete continuum) models. Since existing hybrid models are simplifications of the system dynamics, a new karst model (ModBraC) is presented that accounts for unsteady and nonuniform discrete flow in variably saturated conduits employing the Saint-Venant equations. Model performance tests indicate that ModBraC is able to simulate (1) unsteady and nonuniform flow in variably filled conduits, (2) draining and refilling of conduits with stable transition between free-surface and pressurized flow and correct storage representation, (3) water exchange between matrix and variably filled conduits, and (4) discharge routing through branched and intermeshed conduit networks. Subsequently, ModBraC is applied to an idealized catchment to investigate the significance of free-surface flow representation. A parameter study is conducted with two different initial conditions: (1) pressurized flow and (2) free-surface flow. If free-surface flow prevails, the systems is characterized by (1) a time lag for signal transmission, (2) a typical spring discharge pattern representing the transition from pressurized to free-surface flow, and (3) a reduced conduit-matrix interaction during free-surface flow. Copyright 2011 by the American Geophysical Union.

Water Resources Research↗
USGS2011

O2 reduction and denitrification rates in shallow aquifers

O 2 reduction and denitrification rates were determined in shallow aquifers of 12 study areas representing a wide range in sedimentary environments and climatic conditions. Zero‐ and first‐order rates were determined by relating reactant or product concentrations to apparent groundwater age. O 2 reduction rates varied widely within and between sites, with zero‐order rates ranging from <3 μmol L −1 yr −1 to more than 140 μmol L −1 yr −1 and first‐order rates ranging from 0.02 to 0.27 yr −1 . Moderate denitrification rates (10–100 μmol N L −1 yr −1 ; 0.06–0.30 yr −1 ) were observed in most areas with O 2 concentrations below 60 μmol L −1 , while higher rates (>100 μmol N L −1 yr −1 ; >0.36 yr −1 ) occur when changes in lithology result in a sharp increase in the supply of electron donors. Denitrification lag times (i.e., groundwater travel times prior to the onset of denitrification) ranged from <20 yr to >80 yr. The availability of electron donors is indicated as the primary factor affecting O 2 reduction rates. Concentrations of dissolved organic carbon (DOC) and/or sulfate (an indicator of sulfide oxidation) were positively correlated with groundwater age at sites with high O 2 reduction rates and negatively correlated at sites with lower rates. Furthermore, electron donors from recharging DOC are not sufficient to account for appreciable O 2 and nitrate reduction. These relations suggest that lithologic sources of DOC and sulfides are important sources of electrons at these sites but surface‐derived sources of DOC are not. A review of published rates suggests that denitrification tends to occur more quickly when linked with sulfide oxidation than with carbon oxidation.

Water Resources Research↗