USGS Science⌕ Search

SEARCH · USGS Science

Results for “Hydrology”

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 757 records · Page 42Linked to original sources
USGS1999

Evaluating the use of “goodness‐of‐fit” measures in hydrologic and hydroclimatic model validation

Correlation and correlation‐based measures (e.g., the coefficient of determination) have been widely used to evaluate the “goodness‐of‐fit” of hydrologic and hydroclimatic models. These measures are oversensitive to extreme values (outliers) and are insensitive to additive and proportional differences between model predictions and observations. Because of these limitations, correlation‐based measures can indicate that a model is a good predictor, even when it is not. In this paper, useful alternative goodness‐of‐fit or relative error measures (including the coefficient of efficiency and the index of agreement) that overcome many of the limitations of correlation‐based measures are discussed. Modifications to these statistics to aid in interpretation are presented. It is concluded that correlation and correlation‐based measures should not be used to assess the goodness‐of‐fit of a hydrologic or hydroclimatic model and that additional evaluation measures (such as summary statistics and absolute error measures) should supplement model evaluation tools.

Water Resources Research↗
USGS2002

Using chemical, hydrologic, and age dating analysis to delineate redox processes and flow paths in the riparian zone of a glacial outwash aquifer‐stream system

A combination of chemical and dissolved gas analyses, chlorofluorocarbon age dating, and hydrologic measurements were used to determine the degree to which biogeochemical processes in a riparian wetland were responsible for removing NO 3 − from groundwaters discharging to the Otter Tail River in west central Minnesota. An analysis of river chemistry and flow data revealed that NO 3 − concentrations in the river increased in the lower half of the 8.3 km study reach as the result of groundwater discharge to the river. Groundwater head measurements along a study transect through the riparian wetland revealed a zone of groundwater discharge extending out under the river. On the basis of combined chemical, dissolved gas, age date, and hydrologic results, it was determined that water chemistry under the riparian wetland was controlled largely by upgradient groundwaters that followed flow paths up to 16 m deep and discharged under the wetland, creating a pattern of progressively older, more chemically reduced, low NO 3 − water the farther one progressed from the edge of the wetland toward the river. These findings pose challenges for researchers investigating biogeochemical processes in riparian buffer zones because the progressively older groundwaters entered the aquifer in earlier years when less NO 3 − fertilizer was being used. NO 3 − concentrations originally present in the groundwater had also decreased in the upgradient aquifer as a result of denitrification and progressively stronger reducing conditions there. The resulting pattern of decreasing NO 3 − concentrations across the riparian zone may be incorrectly interpreted as evidence of denitrification losses there instead of in the upgradient aquifer. Consequently, it is important to understand the hydrogeologic setting and age structure of the groundwaters being sampled in order to avoid misinterpreting biogeochemical processes in riparian zones.

Water Resources Research↗
USGS2009

Interactions between soil thermal and hydrological dynamics in the response of Alaska ecosystems to fire disturbance

Soil temperature and moisture are important factors that control many ecosystem processes. However, interactions between soil thermal and hydrological processes are not adequately understood in cold regions, where the frozen soil, fire disturbance, and soil drainage play important roles in controlling interactions among these processes. These interactions were investigated with a new ecosystem model framework, the dynamic organic soil version of the Terrestrial Ecosystem Model, that incorporates an efficient and stable numerical scheme for simulating soil thermal and hydrological dynamics within soil profiles that contain a live moss horizon, fibrous and amorphous organic horizons, and mineral soil horizons. The performance of the model was evaluated for a tundra burn site that had both preburn and postburn measurements, two black spruce fire chronosequences (representing space-for-time substitutions in well and intermediately drained conditions), and a poorly drained black spruce site. Although space-for-time substitutions present challenges in model-data comparison, the model demonstrates substantial ability in simulating the dynamics of evapotranspiration, soil temperature, active layer depth, soil moisture, and water table depth in response to both climate variability and fire disturbance. Several differences between model simulations and field measurements identified key challenges for evaluating/improving model performance that include (1) proper representation of discrepancies between air temperature and ground surface temperature; (2) minimization of precipitation biases in the driving data sets; (3) improvement of the measurement accuracy of soil moisture in surface organic horizons; and (4) proper specification of organic horizon depth/properties, and soil thermal conductivity.

Alaska↗
USGS2008

Sources, transformations, and hydrological processes that control stream nitrate and dissolved organic matter concentrations during snowmelt in an upland forest

We explored catchment processes that control stream nutrient concentrations at an upland forest in northeastern Vermont, USA, where inputs of nitrogen via atmospheric deposition are among the highest in the nation and affect ecosystem functioning. We traced sources of water, nitrate, and dissolved organic matter (DOM) using stream water samples collected at high frequency during spring snowmelt. Hydrochemistry, isotopic tracers, and end‐member mixing analyses suggested the timing, sources, and source areas from which water and nutrients entered the stream. Although stream‐dissolved organic carbon (DOC) and dissolved organic nitrogen (DON) both originated from leaching of soluble organic matter, flushing responses between these two DOM components varied because of dynamic shifts of hydrological flow paths and sources that supply the highest concentrations of DOC and DON. High concentrations of stream water nitrate originated from atmospheric sources as well as nitrified sources from catchment soils. We detected nitrification in surficial soils during late snowmelt which affected the nitrate supply that was available to be transported to streams. However, isotopic tracers showed that the majority of nitrate in upslope surficial soil waters after the onset of snowmelt originated from atmospheric sources. A fraction of the atmospheric nitrogen was directly delivered to the stream, and this finding highlights the importance of quick flow pathways during snowmelt events. These findings indicate that interactions among sources, transformations, and hydrologic transport processes must be deciphered to understand why concentrations vary over time and over space as well as to elucidate the direct effects of human activities on nutrient dynamics in upland forest streams.

Water Resources Research↗
USGS2010

A rapid method for the measurement of sulfur hexafluoride (SF6), trifluoromethyl sulfur pentafluoride (SF5CF3), and Halon 1211 (CF2ClBr) in hydrologic tracer studies

A rapid headspace method for the simultaneous laboratory determination of intentionally introduced hydrologic tracers, sulfur hexafluoride (SF6), trifluoromethyl sulfur pentafluoride (SF5CF3), Halon 1211 (CF2ClBr), and other halocarbons in water and gases is described. The high sensitivity of the procedure allows for introduction of minimal tracer mass (a few grams) into hydrologic systems with a large dynamic range of analytical detection (dilutions to 1:108). Analysis times by gas chromatography with electron capture detector are less than 1 min for SF6; about 2 min for SF6 and SF5CF3; and 4 min for SF6, SF5CF3, and Halon 1211. Many samples can be rapidly collected, preserved in stoppered septum bottles, and analyzed at a later time in the laboratory. Examples are provided showing the effectiveness of the gas tracer test studies in varied hydrogeological settings.

Geochemistry, Geophysics, Geosystems↗
USGS2018

Enhancement of a parsimonious water balance model to simulate surface hydrology in a glacierized watershed

The U.S. Geological Survey monthly water balance model (MWBM) was enhanced with the capability to simulate glaciers in order to make it more suitable for simulating cold region hydrology. The new model, MWBMglacier, is demonstrated in the heavily glacierized and ecologically important Copper River watershed in Southcentral Alaska. Simulated water budget components compared well to satellite‐based observations and ground measurements of streamflow, evapotranspiration, snow extent, and total water storage, with differences ranging from 0.2% to 7% of the precipitation flux. Nash Sutcliffe efficiency for simulated and observed streamflow was greater than 0.8 for six of eight stream gages. Snow extent matched satellite‐based observations with Nash Sutcliffe efficiency values of greater than 0.89 in the four Copper River ecoregions represented. During the simulation period 1949 to 2009, glacier ice melt contributed 25% of total runoff, ranging from 12% to 45% in different tributaries, and glacierized area was reduced by 6%. Statistically significant ( p < 0.05) decreasing and increasing trends in annual glacier mass balance occurred during the multidecade cool and warm phases of the Pacific Decadal Oscillation, respectively, reinforcing the link between climate perturbations and glacier mass balance change. The simulations of glaciers and total runoff for a large, remote region of Alaska provide useful data to evaluate hydrologic, cryospheric, ecologic, and climatic trends. MWBM glacier is a valuable tool to understand when, and to what extent, streamflow may increase or decrease as glaciers respond to a changing climate.

Journal of Geophysical Research F: Earth Surface↗
USGS2019

Hillslope hydrology in global change research and earth system modeling

Earth System Models (ESMs) are essential tools for understanding and predicting global change, but they cannot explicitly resolve hillslope‐scale terrain structures that fundamentally organize water, energy, and biogeochemical stores and fluxes at subgrid scales. Here we bring together hydrologists, Critical Zone scientists, and ESM developers, to explore how hillslope structures may modulate ESM grid‐level water, energy, and biogeochemical fluxes. In contrast to the one‐dimensional (1‐D), 2‐ to 3‐m deep, and free‐draining soil hydrology in most ESM land models, we hypothesize that 3‐D, lateral ridge‐to‐valley flow through shallow and deep paths and insolation contrasts between sunny and shady slopes are the top two globally quantifiable organizers of water and energy (and vegetation) within an ESM grid cell. We hypothesize that these two processes are likely to impact ESM predictions where (and when) water and/or energy are limiting. We further hypothesize that, if implemented in ESM land models, these processes will increase simulated continental water storage and residence time, buffering terrestrial ecosystems against seasonal and interannual droughts. We explore efficient ways to capture these mechanisms in ESMs and identify critical knowledge gaps preventing us from scaling up hillslope to global processes. One such gap is our extremely limited knowledge of the subsurface, where water is stored (supporting vegetation) and released to stream baseflow (supporting aquatic ecosystems). We conclude with a set of organizing hypotheses and a call for global syntheses activities and model experiments to assess the impact of hillslope hydrology on global change predictions.

Water Resources Research↗
USGS2019

Assessing the feasibility of satellite-based thresholds for hydrologically driven landsliding

Elevated soil moisture and heavy precipitation contribute to landslides worldwide. These environmental variables are now being resolved with satellites at spatiotemporal scales that could offer new perspectives on the development of landslide warning systems. However, the application of these data to hydro-meteorological thresholds (which account for antecedent soil moisture and rainfall) first need to be evaluated with respect to proven, direct measurement-based thresholds that use rain gauges and in situ soil moisture sensors. Here, we compare ground-based hydrologic data to overlapping satellite-based data before, during, and after a recent season of widespread shallow landsliding in the San Francisco Bay Area (California, USA). We then explore how the remotely sensed information could be used to empirically define hypothetical thresholds for shallow landsliding. We find that the ground-based thresholds developed with a single monitoring station show superior performance because the in situ soil saturation data better reflect the gravity-dominated subsurface flow conditions that are characteristic of hillslopes during the rainy season. Although the satellite-based thresholds can identify most of the landslide days, they include a greater number of false alarms due to overestimates of soil moisture between major storm events. To avoid the type of false alarms that are characteristic of our satellite-based thresholds, further post-processing of the near-surface hydrologic response data to better reflect gravity-dominated drainage should be integrated into satellite-based model outputs. Our results encourage further deployment of ground stations in landslide-prone terrain and cautious exploration of satellite-based hydro-meteorological thresholds where in situ networks are nonexistent.

California↗
USGS2020

Hydrologically induced deformation in Long Valley Caldera and adjacent Sierra Nevada

Vertical and horizontal components of GNSS displacements in the Long Valley Caldera and adjacent Sierra Nevada range show a clear correlation with hydrological trends at both multiyear and seasonal time scales. We observe a clear vertical and horizontal seasonal deformation pattern primarily attributable to the solid earth response to hydrological surface loading at large-to-regional (Sierra Nevada range) scales. Several GNSS sites, located at the eastern edge of the Sierra Nevada along the southwestern rim of Long Valley Caldera, also show significant horizontal deformation that cannot be explained by elastic deformation from surface loading. Due to the location of these sites and the strong correlation between their horizontal displacements and spring discharge, we hypothesize that this deformation reflects poroelastic processes related to snowmelt runoff water infiltrating into the Sierra Nevada slopes around Long Valley Caldera. Interestingly, this is also an area where water infiltrates to feed the local hydrothermal system, and where snowmelt-induced earthquake swarms have been recently detected.

California, Nevada↗
USGS2020

Twelve-year dynamics and rainfall thresholds for alternating creep and rapid movement of the Hooskanaden landslide from integrating InSAR, pixel offset tracking, and borehole and hydrological measurements

The Hooskanaden landslide is a large (~600 m wide × 1,300 m long), deep (~30 – 45 m) slide located in southwestern Oregon. Since 1958, it has had five moderate/major movements that catastrophically damaged the intersecting U.S. Highway 101, along with persistent slow wet‐season movements and a long‐term accelerating trend due to coastal erosion. Multiple remote sensing approaches, borehole measurements, and hydrological observations have been integrated to interpret the motion behaviors of the slide. Pixel offset tracking of both Sentinel‐1 and Sentinel‐2 images was carried out to reconstruct the 3‐D displacement field of the 2019 major event, and the results agree well with field measurements. A 12‐year displacement history of the landslide from 2007 to 2019 has been retrieved by incorporating offsets from Light Detection and Ranging (LiDAR) digital elevation model (DEM) gradients and Interferometric Synthetic Aperture Radar (InSAR) processing of ALOS and Sentinel‐1 images. Comparisons with daily/hourly ground precipitation reveal that the motion dynamics are predominantly controlled by intensity and temporal pattern of rainfall. A new empirical threefold rainfall threshold was therefore proposed to forecast the dates for the moderate/major movements. This threshold relies upon antecedent water‐year and previous 3‐day and daily precipitation and was able to represent observed movement periods well. Adaptation of our threshold methodology could prove useful for other large, deep landslides for which temporal forecasting has long been generally intractable. The averaged characteristic hydraulic conductivity and diffusivity were estimated as 6.6 × 10 −6 m/s and 6.6 × 10 −4 m 2 /s, respectively, based on the time lags between rainfall pulses and slide accelerations. Hydrologic modeling using these parameters helps to explain the ability of the new rainfall threshold.

Oregon↗
USGS2022

Heterogeneous patterns of aged organic carbon export driven by hydrologic flow paths, soil texture, fire, and thaw in discontinuous permafrost headwaters

Climate change is thawing and potentially mobilizing vast quantities of organic carbon (OC) previously stored for millennia in permafrost soils of northern circumpolar landscapes. Climate-driven increases in fire and thermokarst may play a key role in OC mobilization by thawing permafrost and promoting transport of OC. Yet, the extent of OC mobilization and mechanisms controlling terrestrial-aquatic transfer are unclear. We demonstrate that hydrologic transport of soil dissolved OC (DOC) from the active layer and thawing permafrost to headwater streams is extremely heterogeneous and regulated by the interactions of soils, seasonal thaw, fire, and thermokarst. Repeated sampling of streams in eight headwater catchments of interior Alaska showed that the mean age of DOC for each stream ranges widely from modern to ∼2,000 years B.P. Together, an endmember mixing model and nonlinear, generalized additive models demonstrated that Δ 14 C-DOC signature (and mean age) increased from spring to fall, and was proportional to hydrologic contributions from a solute-rich water source, related to presumed deeper flow paths found predominantly in silty catchments. This relationship was correlated with and mediated by catchment properties. Mean DOC ages were older in catchments with >50% burned area, indicating that fire is also an important explanatory variable. These observations underscore the high heterogeneity in aged C export and difficulty of extrapolating estimates of permafrost-derived DOC export from watersheds to larger scales. Our results provide the foundation for developing a conceptual model of permafrost DOC export necessary for advancing understanding and prediction of land-water C exchange in changing boreal landscapes.

Alaska↗
USGS2022

Hydrologic connectivity and residence time affect the sediment trapping efficiency and dissolved oxygen concentrations of the Atchafalaya River Basin

Little is known about water movement, volume, or residence time (RT), and how those characteristics affect sediment trapping efficiency (TE) and dissolved oxygen concentrations (DO) in the United States' largest remaining bottomland hardwood swamp, the Atchafalaya River Basin. To better understand these dynamics, this study used bathymetry, lidar, and stage records to determine volumes in the Basin's hydrologically distinct water management units (WMUs). Discharge measurements determined flow distribution and RT. Residence time was compared with DO to identify conditions that coincided with DO increases or decreases. Suspended sediment concentrations (SSC) were used to determine TE relative to calculated and measured discharge and RT. Discharge through units (85–2,200 m 3 /s) and RT (0.37–231 d) depended on connectivity and river stage. At high stages, with water temperatures >20°C, DO in the largest WMU declined by −0.21 mg/l/day. DO trends indicated less well-connected areas of the WMU contributed hypoxic waters as the flood wave lengthened and stages fell. In the two WMUs examined for TE, TE (−266% to 99% and up to 38 Gg/day) correlated with hydrologic connectivity, SSC, RT, water volume, and, in one WMU, discharge losses. Long RT and high TE indicated a high potential to process nutrients. These relationships varied among WMUs. Large volumes of sediment-laden water moving over the floodplain combined with long RT, high TE, and hypoxia indicate that this ecosystem has continental-scale importance in reducing nutrient loads to the northern Gulf of Mexico. Reports from other systems suggest similar processes may be operating on other large river floodplains globally.

Louisiana↗
USGS2022

Hydrologic and landscape controls on dissolved organic matter composition across western North American Arctic lakes

Northern high-latitude lakes are hotspots for cycling dissolved organic carbon (DOC) inputs from allochthonous sources to the atmosphere. However, the spatial distribution of lake dissolved organic matter (DOM) is largely unknown across Arctic-boreal regions with respect to the surrounding landscape. We expand on regional studies of northern high-latitude DOM composition by integrating DOC concentrations, optical properties, and molecular-level characterization from lakes spanning the Canadian Taiga to the Alaskan Tundra. Lakes were sampled during the summer from July to early September to capture the growing season. DOM became more optically processed and molecular-level aromaticity increased northward across the Canadian Shield to the southern Arctic and from interior Alaska to the Tundra, suggesting relatively greater DOM incorporation from allochthonous sources. Using water isotopes (δ 18 O-H 2 O), we report a weak overall trend of increasing DOC and decreasing aromaticity in lakes that were hydrologically isolated from the landscape and enriched in δ 18 O-H 2 O, while within-region trends were stronger and varied depending on the landscape. Finally, DOC correlated weakly with chromophoric dissolved organic matter (CDOM) across the study sites, suggesting that autochthonous and photobleached DOM were a major component of the DOC in these regions; however, some of the northernmost and wetland-dominated lakes followed pan-Arctic riverine DOC-CDOM relationships, indicating strong contributions from allochthonous inputs. As many lakes across the North American Arctic are experiencing changes in temperature and precipitation, we expect the proportions of allochthonous and autochthonous DOM to respond with aquatic optical browning with greater landscape connectivity and more internally produced DOM in hydrologically isolated lakes.

Global Biogeochemical Cycles↗
USGS2024

Quantifying compound and nonlinear effects of hurricane-induced flooding using a dynamically coupled hydrological-ocean model

We recently developed a dynamically coupled hydrological-ocean modeling system that provides seamless coverage across the land-ocean continuum during hurricane-induced compound flooding. This study introduced a local inertial equation and a diagonal flow algorithm to the overland routing of the coupled system’s hydrology model (WRF-Hydro). Using Hurricane Florence (2018) as a test case, the performance of the coupled model was significantly improved, evidenced by its enhanced capability of capturing backwater and increased water level simulation accuracy and stability. With four model experiments, we present a framework to detangle, define, and quantify compound and nonlinear effects. The results revealed that the flood peaks in the lower Cape Fear River Basin and the coastal waters were contributed by inland flooding and storm surge, respectively. These two processes had comparable contributions to the flooding in the Cape Fear River Estuary. The compound effect was identified when the flood levels resulting from the combination of land and ocean processes surpassed those caused by an individual process alone. The compound effect during Hurricane Florence exhibited limited impact on flood peaks, primarily due to the time lag between the peaks of the storm surge and the inland flooding. In the period between the two peaks, the compound effect was salient and significantly impacted the magnitude and variation of the flood level. The nonlinear effect, defined as the difference between the compound flood level and the superposition of storm surge and inland flooding water levels, reduced flood levels in the river channels while increasing flood levels on the floodplain.

Water Resources Research↗
USGS2024

Climate, hydrology, and nutrients control the seasonality of Si concentrations in rivers

The seasonal behavior of fluvial dissolved silica (DSi) concentrations, termed DSi regime , mediates the timing of DSi delivery to downstream waters and thus governs river biogeochemical function and aquatic community condition. Previous work identified five distinct DSi regimes across rivers spanning the Northern Hemisphere, with many rivers exhibiting multiple DSi regimes over time. Several potential drivers of DSi regime behavior have been identified at small scales, including climate, land cover, and lithology, and yet the large-scale spatiotemporal controls on DSi regimes have not been identified. We evaluate the role of environmental variables on the behavior of DSi regimes in nearly 200 rivers across the Northern Hemisphere using random forest models. Our models aim to elucidate the controls that give rise to (a) average DSi regime behavior, (b) interannual variability in DSi regime behavior (i.e., Annual DSi regime), and (c) controls on DSi regime shape (i.e., minimum and maximum DSi concentrations). Average DSi regime behavior across the period of record was classified accurately 59% of the time, whereas Annual DSi regime behavior was classified accurately 80% of the time. Climate and primary productivity variables were important in predicting Average DSi regime behavior, whereas climate and hydrologic variables were important in predicting Annual DSi regime behavior. Median nitrogen and phosphorus concentrations were important drivers of minimum and maximum DSi concentrations, indicating that these macronutrients may be important for seasonal DSi drawdown and rebound. Our findings demonstrate that fluctuations in climate, hydrology, and nutrient availability of rivers shape the temporal availability of fluvial DSi.

Journal of Geophysical Research Biogeosciences↗
USGS2025

Uncertainty quantification of geophysical and hydrologic parameters estimated from borehole nuclear magnetic resonance data

Borehole nuclear magnetic resonance (bNMR) data are typically used to infer in situ hydrologic properties. Partial water content as a function of pore size is estimated by fitting the measured NMR response to a multi-exponential T 2 distribution, and the sum of estimated T 2 amplitudes equals the total volumetric water content. From these estimated parameters, several empirical relationships are commonly used to infer hydraulic conductivity from the NMR-estimated water content and T 2 distribution. Often, parameters are estimated through deterministic inversion methods that produce a single best-fit estimate, but do not reflect uncertainties in model parameters. Here, a Bayesian Markov chain Monte Carlo (McMC) approach for analyzing bNMR data is developed that allows for comprehensive uncertainty quantification of NMR parameters and derived hydrologic properties. The underlying model that describes the T 2 distribution is defined by a set of spline interpolation points. The number of interpolation points is allowed to vary in a trans-dimensional algorithm that naturally favors simple models with fewer interpolation points, allowing the data to inform the necessary level of model complexity. Additionally, data error is estimated as an unknown parameter. Analysis of the ensemble of models output from the McMC algorithm provides useful details on the range of plausible T 2 distributions that can fit a measured bNMR decay curve, as well as uncertainty estimates of total water content. The ensemble of NMR parameters can also be propagated through commonly used relationships to produce uncertainty estimates on derived parameters such as bound/capillary/mobile water content or hydraulic conductivity.

JGR Machine Learning and Computation↗
USGS2010

Hydrologic processes controlling sulfate mobility in a small forested watershed

Hydrologic controls on sulfate mobility were investigated in a forested catchment in the Georgia Piedmont using a watershed mass balance approach. Variations in annual sulfate export were governed primarily by differences in runoff rather than by differences in sulfate deposition or in total annual precipitation. However, 2 years with similar total runoff had substantially different sulfate export. A residual analysis indicated that a shift in the concentration‐discharge relationship accounted for 72% of the difference in sulfate export and that a change in the pattern of discharge accounted for the remainder of the difference. Stream water sulfate concentrations reflected past hydrologic conditions. Concentrations at the same discharge were higher following an extended dry period than following average periods. The elevation in stream water sulfate concentrations following dry periods persisted for several months. The influence of rainfall patterns on sulfate export underscores the need for long‐term records to adequately characterize the acidification status of the watershed and to understand trends in water quality.

Georgia↗
USGS2010

Deducing the distribution of terminal electron-accepting processes in hydrologically diverse groundwater systems

The distribution of microbially mediated terminal electron-accepting processes (TEAPs( was investigated in four hydrologically diverse groundwater systems by considering patterns of electron acceptor (nitrate, sulfate) consumption, intermediate product (hydrogen (H 2 )) concentrations, and final product (ferrous iron, sulfide, and methane) production. In each hydrologic system a determination of predominant TEAPs could be arrived at, but the level of confidence appropriate for each determination differed. In a portion of the lacustrine aquifer of the San Joaquin Valley, for example, all three indicators (sulfate concentrations decreasing, H 2 concentrations in the 1–2 nmol range, and sulfide concentrations increasing along flow paths identified sulfate reduction as the predominant TEAP, leading to a high degree of confidence in the determination. In portions of the Floridan aquifer and a petroleum hydrocarbon-contaminated aquifer, sulfate reduction and methanogenesis are indicated by production of sulfide and methane, and hydrogen oncentrations in the 1–4 nmol and 5–14 nmol range, respectively. However, because electron acceptor consumption could not be documented in these systems, less confidence is warranted in the TEAP determination. In the Black Creek aquifer, no pattern of sulfate consumption and sulfide production were observed, but H2 concentrations indicated sulfate reduction as the predominant TEAP. In this case, where just a single line of evidence is available, the least confidence in the TEAP diagnosis is justified. Because this methodology is based on measurable water chemistry parameters and upon the physiology of microbial electron transfer processes, it provides a better description of predominant redox processes in groundwater systems than more traditional Eh-based methods.

Water Resources Research↗