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 523 records · Page 29Linked to original sources
USGS2009

Water velocity and the nature of critical flow in large rapids on the Colorado River, Utah

Rapids are an integral part of bedrock‐controlled rivers, influencing aquatic ecology, geomorphology, and recreational value. Flow measurements in rapids and high‐gradient rivers are uncommon because of technical difficulties associated with positioning and operating sufficiently robust instruments. In the current study, detailed velocity, water surface, and bathymetric data were collected within rapids on the Colorado River in eastern Utah. With the water surface survey, it was found that shoreline‐based water surface surveys may misrepresent the water surface slope along the centerline of a rapid. Flow velocities were measured with an ADCP and an electronic pitot‐static tube. Integrating multiple measurements, the ADCP returned velocity data from the entire water column, even in sections of high water velocity. The maximum mean velocity measured with the ADCP was 3.7 m/s. The pitot‐static tube, while capable of only point measurements, quantified velocity 0.39 m below the surface. The maximum mean velocity measured with the pitot tube was 5.2 m/s, with instantaneous velocities up to 6.5 m/s. Analysis of the data showed that flow was subcritical throughout all measured rapids with a maximum measured Froude number of 0.7 in the largest measured rapids. Froude numbers were highest at the entrance of a given rapid, then decreased below the first breaking waves. In the absence of detailed bathymetric and velocity data, the Froude number in the fastest‐flowing section of a rapid was estimated from near‐surface velocity and depth soundings alone.

Water Resources Research↗
USGS2010

Comparison of aquifer characterization approaches through steady state groundwater model validation: A controlled laboratory sandbox study

Groundwater modeling has become a vital component to water supply and contaminant transport investigations. An important component of groundwater modeling under steady state conditions is selecting a representative hydraulic conductivity (K) estimate or set of estimates which defines the K field of the studied region. Currently, there are a number of characterization approaches to obtain K at various scales and in varying degrees of detail, but there is a paucity of information in terms of which characterization approach best predicts flow through aquifers or drawdowns caused by some drawdown inducing events. The main objective of this paper is to assess K estimates obtained by various approaches by predicting drawdowns from independent cross-hole pumping tests and total flow rates through a synthetic heterogeneous aquifer from flow-through tests. Specifically, we (1) characterize a synthetic heterogeneous aquifer built in the sandbox through various techniques (permeameter analyses of core samples, single-hole, cross-hole, and flow-through testing), (2) obtain mean K fields through traditional analysis of test data by treating the medium to be homogeneous, (3) obtain heterogeneous K fields through kriging and steady state hydraulic tomography, and (4) conduct forward simulations of 16 independent pumping tests and six flowthrough tests using these homogeneous and heterogeneous K fields and comparing them to actual data. Results show that the mean K and heterogeneous K fields estimated through kriging of small-scale K data (core and single-hole tests) yield biased predictions of drawdowns and flow rates in this synthetic heterogeneous aquifer. In contrast, the heterogeneous K distribution or ?K tomogram? estimated via steady state hydraulic tomography yields excellent predictions of drawdowns of pumping tests not used in the construction of the tomogram and very good estimates of total flow rates from the flowthrough tests. These results suggest that steady state groundwater model validation is possible in this laboratory sandbox aquifer if the heterogeneous K distribution and forcing functions (boundary conditions and source/sink terms) are characterized sufficiently. ?? 2010 by the American Geophysical Union.

Water Resources Research↗
USGS2010

Quantification of surface water and groundwater flows to open‐ and closed‐basin lakes in a headwaters watershed using a descriptive oxygen stable isotope model

Accurate quantification of hydrologic fluxes in lakes is important to resource management and for placing hydrologic solute flux in an appropriate biogeochemical context. Water stable isotopes can be used to describe water movements, but they are typically only effective in lakes with long water residence times. We developed a descriptive time series model of lake surface water oxygen‐18 stable isotope signature ( δL ) that was equally useful in open‐ and closed‐basin lakes with very different hydrologic residence times. The model was applied to six lakes, including two closed‐basin lakes and four lakes arranged in a chain connected by a river, located in a headwaters watershed. Groundwater discharge was calculated by manual optimization, and other hydrologic flows were constrained by measured values including precipitation, evaporation, and streamflow at several stream gages. Modeled and observed δL were highly correlated in all lakes ( r = 0.84–0.98), suggesting that the model adequately described δL in these lakes. Average modeled stream discharge at two points along the river, 16,000 and 11,800 m 3 d −1 , compares favorably with synoptic measurement of stream discharge at these sites, 17,600 and 13,700 m 3 d −1 , respectively. Water yields in this watershed were much higher, 0.23–0.45 m, than water yields calculated from gaged streamflow in regional rivers, approximately 0.10 m, suggesting that regional groundwater discharge supports water flux through these headwaters lakes. Sensitivity and robustness analyses also emphasized the importance of considering hydrologic residence time when designing a sampling protocol for stable isotope use in lake hydrology studies.

Water Resources Research↗
USGS2010

Evaluation of aquifer heterogeneity effects on river flow loss using a transition probability framework

River-aquifer exchange is considered within a transition probability framework along the Rio Grande in Albuquerque, New Mexico, to provide a stochastic estimate of aquifer heterogeneity and river loss. Six plausible hydrofacies configurations were determined using categorized drill core and wetland survey data processed through the TPROGS geostatistical package. A base case homogeneous model was also constructed for comparison. River loss was simulated for low, moderate, and high Rio Grande stages and several different riverside drain stage configurations. Heterogeneity effects were quantified by determining the mean and variance of the K field for each realization compared to the root-mean-square (RMS) error of the observed groundwater head data. Simulation results showed that the heterogeneous models produced smaller estimates of loss than the homogeneous approximation. Differences between heterogeneous and homogeneous model results indicate that the use of a homogeneous K in a regional-scale model may result in an overestimation of loss but comparable RMS error. We find that the simulated river loss is dependent on the aquifer structure and is most sensitive to the volumetric proportion of fines within the river channel. Copyright 2010 by the American Geophysical Union.

Water Resources Research↗
USGS2009

Introduction to special section on impacts of land use change on water resources

Changes in land use have potentially large impacts on water resources, yet quantifying these impacts remains among the more challenging problems in hydrology. Water, food, energy, and climate are linked through complex webs of direct and indirect effects and feedbacks. Land use is undergoing major changes due not only to pressures for more efficient food, feed, and fiber production to support growing populations but also due to policy shifts that are creating markets for biofuel and agricultural carbon sequestration. Hydrologic systems embody flows of water, solutes, sediments, and energy that vary even in the absence of human activity. Understanding land use impacts thus necessitates integrated scientific approaches. Field measurements, remote sensing, and modeling studies are shedding new light on the modes and mechanisms by which land use changes impact water resources. Such studies can help deconflate the interconnected influences of human actions and natural variations on the quantity and quality of soil water, surface water, and groundwater, past, present, and future.

Water Resources Research↗
USGS2010

Impacts of precipitation seasonality and ecosystem types on evapotranspiration in the Yukon River Basin, Alaska

Evapotranspiration (ET) is the largest component of water loss from terrestrial ecosystems; however, large uncertainties exist when estimating the temporal and spatial variations of ET because of concurrent shifts in the magnitude and seasonal distribution of precipitation as well as differences in the response of ecosystem ET to environmental variabilities. In this study, we examined the impacts of precipitation seasonality and ecosystem types on ET quantified by eddy covariance towers from 2002 to 2004 in three ecosystems (grassland, deciduous broadleaf forest, and evergreen needleleaf forest) in the Yukon River Basin, Alaska. The annual precipitation changed greatly in both magnitude and seasonal distribution through the three investigated years. Observations and model results showed that ET was more sensitive to precipitation scarcity in the early growing season than in the late growing season, which was the direct result of different responses of ET components to precipitation in different seasons. The results demonstrated the importance of seasonal variations of precipitation in regulating annual ET and overshadowing the function of annual precipitation. Comparison of ET among ecosystems over the growing season indicated that ET was largest in deciduous broadleaf, intermediate in evergreen needleleaf, and lowest in the grassland ecosystem. These ecosystem differences in ET were related to differences in successional stages and physiological responses.

Water Resources Research↗
USGS2010

Emergency use of groundwater as a backup supply: Quantifying hydraulic impacts and economic benefits

Groundwater can play an important role in water‐supply emergency planning. A framework is presented for assessing the hydraulic impacts and associated costs of using groundwater as a backup supply when imported‐water deliveries are disrupted, and for quantifying the emergency benefits of groundwater management strategies that enable better response to such disruptions. Response functions are derived, which relate additional groundwater pumpage during water‐supply emergencies to impacts such as increased pumping costs, subsidence, and seawater intrusion. Monte Carlo analysis is employed to estimate the incremental costs of using groundwater as a backup supply. The emergency benefits of alternative groundwater management strategies are computed for different expected durations of imported water disruption, percentages of imported water replaced by groundwater, and threshold drawdowns for subsidence impacts. The methodology is applied to the coastal Los Angeles Basin. For this case study, emergency benefits of artificial recharge strategies are dominated by reduction of potential subsidence costs. The variance of the results also is primarily due to subsidence effects. Incorporation of probability distributions reflecting a larger expected use of groundwater during the imported‐water disruption results in higher estimated emergency benefits of artificial recharge strategies. The framework presented for quantifying incremental costs and economic benefits of using groundwater as a backup supply could be applied to a broad range of water emergency planning decisions.

California↗
USGS2010

Use of regression‐based models to map sensitivity of aquatic resources to atmospheric deposition in Yosemite National Park, USA

An abundance of exposed bedrock, sparse soil and vegetation, and fast hydrologic flushing rates make aquatic ecosystems in Yosemite National Park susceptible to nutrient enrichment and episodic acidification due to atmospheric deposition of nitrogen (N) and sulfur (S). In this study, multiple linear regression (MLR) models were created to estimate fall‐season nitrate and acid neutralizing capacity (ANC) in surface water in Yosemite wilderness. Input data included estimated winter N deposition, fall‐season surface‐water chemistry measurements at 52 sites, and basin characteristics derived from geographic information system layers of topography, geology, and vegetation. The MLR models accounted for 84% and 70% of the variance in surface‐water nitrate and ANC, respectively. Explanatory variables (and the sign of their coefficients) for nitrate included elevation (positive) and the abundance of neoglacial and talus deposits (positive), unvegetated terrain (positive), alluvium (negative), and riparian (negative) areas in the basins. Explanatory variables for ANC included basin area (positive) and the abundance of metamorphic rocks (positive), unvegetated terrain (negative), water (negative), and winter N deposition (negative) in the basins. The MLR equations were applied to 1407 stream reaches delineated in the National Hydrography Data Set for Yosemite, and maps of predicted surface‐water nitrate and ANC concentrations were created. Predicted surface‐water nitrate concentrations were highest in small, high‐elevation cirques, and concentrations declined downstream. Predicted ANC concentrations showed the opposite pattern, except in high‐elevation areas underlain by metamorphic rocks along the Sierran Crest, which had relatively high predicted ANC (>200 μ eq L −1 ). Maps were created to show where basin characteristics predispose aquatic resources to nutrient enrichment and acidification effects from N and S deposition. The maps can be used to help guide development of water‐quality programs designed to monitor and protect natural resources in national parks.

Yosemite Nataional Park↗
USGS2009

A new method for high-resolution characterization of hydraulic conductivity

A new probe has been developed for high-resolution characterization of hydraulic conductivity (K) in shallow unconsolidated formations. The probe was recently applied at the Macrodispersion Experiment (MADE) site in Mississippi where K was rapidly characterized at a resolution as fine as 0.015 m, which has not previously been possible. Eleven profiles were obtained with K varying up to 7 orders of magnitude in individual profiles. Currently, high-resolution (0.015-m) profiling has an upper K limit of 10 m/d; lower-resolution (???0.4-m) mode is used in more permeable zones pending modifications. The probe presents a new means to help address unresolved issues of solute transport in heterogeneous systems. Copyright 2009 by the American Geophysical Union.

Water Resources Research↗
USGS2010

Testing mixing models of old and young groundwater in a tropical lowland rain forest with environmental tracers

We tested three models of mixing between old interbasin groundwater flow (IGF) and young, locally derived groundwater in a lowland rain forest in Costa Rica using a large suite of environmental tracers. We focus on the young fraction of water using the transient tracers CFC‐11, CFC‐12, CFC‐113, SF 6 , 3 H, and bomb 14 C. We measured 3 He, but 3 H/ 3 He dating is generally problematic due to the presence of mantle 3 He. Because of their unique concentration histories in the atmosphere, combinations of transient tracers are sensitive not only to subsurface travel times but also to mixing between waters having different travel times. Samples fall into three distinct categories: (1) young waters that plot along a piston flow line, (2) old samples that have near‐zero concentrations of the transient tracers, and (3) mixtures of 1 and 2. We have modeled the concentrations of the transient tracers using (1) a binary mixing model (BMM) of old and young water with the young fraction transported via piston flow, (2) an exponential mixing model (EMM) with a distribution of groundwater travel times characterized by a mean value, and (3) an exponential mixing model for the young fraction followed by binary mixing with an old fraction (EMM/BMM). In spite of the mathematical differences in the mixing models, they all lead to a similar conceptual model of young (0 to 10 year) groundwater that is locally derived mixing with old (>1000 years) groundwater that is recharged beyond the surface water boundary of the system.

Water Resources Research↗
USGS2010

Mercury dynamics in relation to dissolved organic carbon concentration and quality during high flow events in three northeastern U.S. streams

Mercury (Hg) contamination is widespread in remote areas of the northeastern United States. Forested uplands have accumulated a large reservoir of Hg in soil from decades of elevated anthropogenic deposition that can be released episodically to stream water during high flows. The objective of this study was to evaluate spatial and temporal variations in stream water Hg species and organic matter fractions over a range of hydrologic conditions in three forested upland watersheds (United States). Mercury and organic matter concentrations increased with discharge at all three sites; however, the partitioning of Hg fractions (dissolved versus particulate) differed among sites and seasons. Associated with increased discharge, flow paths shifted from mineral soil under base flow to upper soil horizons. As flow paths shifted, greater concentrations of dissolved organic carbon (DOC) richer in aromatic substances were flushed from upper soil horizons to stream water. The hydrophobic organic matter associated with humic material from upper soils appears to have had a greater capacity to bind Hg. Because of the strong correlation between Hg and DOC, we hypothesize that there was a concurrent shift in the source of Hg with DOC from lower mineral soil to upper soil horizons. Our study suggests that stream discharge is an effective predictor of dissolved total Hg flux.

Water Resources Research↗
USGS2010

Map correlation method: Selection of a reference streamgage to estimate daily streamflow at ungaged catchments

Daily streamflow time series are critical to a very broad range of hydrologic problems. Whereas daily streamflow time series are readily obtained from gaged catchments, streamflow information is commonly needed at catchments for which no measured streamflow information exists. At ungaged catchments, methods to estimate daily streamflow time series typically require the use of a reference streamgage, which transfers properties of the streamflow time series at a reference streamgage to the ungaged catchment. Therefore, the selection of a reference streamgage is one of the central challenges associated with estimation of daily streamflow at ungaged basins. The reference streamgage is typically selected by choosing the nearest streamgage; however, this paper shows that selection of the nearest streamgage does not provide a consistent selection criterion. We introduce a new method, termed the map‐correlation method, which selects the reference streamgage whose daily streamflows are most correlated with an ungaged catchment. When applied to the estimation of daily streamflow at 28 streamgages across southern New England, daily streamflows estimated by a reference streamgage selected using the map‐correlation method generally provides improved estimates of daily streamflow time series over streamflows estimated by the selection and use of the nearest streamgage. The map correlation method could have potential for many other applications including identifying redundancy and uniqueness in a streamgage network, calibration of rainfall runoff models at ungaged sites, as well as for use in catchment classification.

New England↗
USGS2010

Effects of upstream dams versus groundwater pumping on stream temperature under varying climate conditions

The relative impact of a large upstream dam versus in‐reach groundwater pumping on stream temperatures was analyzed for humid, semiarid, and arid conditions with long dry seasons to represent typical climate regions where large dams are present, such as the western United States or eastern Australia. Stream temperatures were simulated using the CE‐QUAL‐W2 water quality model over a 110 km model grid, with the presence or absence of a dam at the top of the reach and pumping in the lower 60 km of the reach. Measured meteorological data from three representative locations were used as model input to simulate the impact of varying climate conditions on streamflow and stream temperature. For each climate condition four hypothetical streamflow scenarios were modeled: (1) natural (no dam or pumping), (2) large upstream dam present, (3) dam with in‐reach pumping, and (4) no dam with pumping, resulting in 12 cases. Dam removal, in the presence or absence of pumping, resulted in significant changes in stream temperature throughout the year for all three climate conditions. From March to August, the presence of a dam caused monthly mean stream temperatures to decrease on average by approximately 3.0°C, 2.5°C, and 2.0°C for the humid, semiarid, and arid conditions, respectively; however, stream temperatures generally increased from September to February. Pumping caused stream temperatures to warm in summer and cool in winter by generally less than 0.5°C because of a smaller pumping‐induced alteration in streamflow relative to the dam. Though the presence or absence of a large dam led to greater changes in stream temperature than the presence or absence of pumping, ephemeral conditions were increased both temporally and spatially because of pumping.

Water Resources Research↗
USGS2010

An approach for modeling sediment budgets in supply-limited rivers

Reliable predictions of sediment transport and river morphology in response to variations in natural and human-induced drivers are necessary for river engineering and management. Because engineering and management applications may span a wide range of space and time scales, a broad spectrum of modeling approaches has been developed, ranging from suspended-sediment "rating curves" to complex three-dimensional morphodynamic models. Suspended sediment rating curves are an attractive approach for evaluating changes in multi-year sediment budgets resulting from changes in flow regimes because they are simple to implement, computationally efficient, and the empirical parameters can be estimated from quantities that are commonly measured in the field (i.e., suspended sediment concentration and water discharge). However, the standard rating curve approach assumes a unique suspended sediment concentration for a given water discharge. This assumption is not valid in rivers where sediment supply varies enough to cause changes in particle size or changes in areal coverage of sediment on the bed; both of these changes cause variations in suspended sediment concentration for a given water discharge. More complex numerical models of hydraulics and morphodynamics have been developed to address such physical changes of the bed. This additional complexity comes at a cost in terms of computations as well as the type and amount of data required for model setup, calibration, and testing. Moreover, application of the resulting sediment-transport models may require observations of bed-sediment boundary conditions that require extensive (and expensive) observations or, alternatively, require the use of an additional model (subject to its own errors) merely to predict the bed-sediment boundary conditions for use by the transport model. In this paper we present a hybrid approach that combines aspects of the rating curve method and the more complex morphodynamic models. Our primary objective was to develop an approach complex enough to capture the processes related to sediment supply limitation but simple enough to allow for rapid calculations of multi-year sediment budgets. The approach relies on empirical relations between suspended sediment concentration and discharge but on a particle size specific basis and also tracks and incorporates the particle size distribution of the bed sediment. We have applied this approach to the Colorado River below Glen Canyon Dam (GCD), a reach that is particularly suited to such an approach because it is substantially sediment supply limited such that transport rates are strongly dependent on both water discharge and sediment supply. The results confirm the ability of the approach to simulate the effects of supply limitation, including periods of accumulation and bed fining as well as erosion and bed coarsening, using a very simple formulation. Although more empirical in nature than standard one-dimensional morphodynamic models, this alternative approach is attractive because its simplicity allows for rapid evaluation of multi-year sediment budgets under a range of flow regimes and sediment supply conditions, and also because it requires substantially less data for model setup and use.

Water Resources Research↗
USGS2010

Effect of surficial disturbance on exchange between groundwater and surface water in nearshore margins

Low‐permeability sediments situated at or near the sediment‐water interface can influence seepage in nearshore margins, particularly where wave energy or currents are minimal. Seepage meters were used to quantify flow across the sediment‐water interface at two lakes where flow was from surface water to groundwater. Disturbance of the sediment bed substantially increased seepage through the sandy sediments of both lakes. Seepage increased by factors of 2.6 to 7.7 following bed disturbance at seven of eight measurement locations at Mirror Lake, New Hampshire, where the sediment representing the greatest restriction to flow was situated at the sediment‐water interface. Although the veneer of low‐permeability sediment was very thin and easily disturbed, accumulation on the bed surface was aided by a physical setting that minimized wind‐generated waves and current. At Lake Belle Taine, Minnesota, where pre‐disturbance downward seepage was smaller than at Mirror Lake, but hydraulic gradients were very large, disturbance of a 20 to 30 cm thick medium sand layer resulted in increases in seepage of 2 to 3 orders of magnitude. Exceptionally large seepage rates, some exceeding 25,000 cm/d, were recorded following bed disturbance. Since it is common practice to walk on the bed while installing or making seepage measurements, disruption of natural seepage rates may be a common occurrence in nearshore seepage studies. Disturbance of the bed should be avoided or minimized when utilizing seepage meters in shallow, nearshore settings, particularly where waves or currents are infrequent or minimal.

Minnesota, New Hampshire↗
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↗