Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hydrologic exchange flows”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Hyporheic Oxygen Dynamics in the East River, Colorado: Insights From an In–Situ, High Frequency Time Series During Two Distinct Flow Seasons

Dissolved oxygen (DO) is critical for aquatic ecosystems, however, few studies have focused on the long-term DO dynamics in hyporheic zones, which are a function of both transport (hydrologic exchange between river and hyporheic zone) and uptake by biogeochemical reactions or respiration. We explore the dynamics of temperature and DO at 10, 20, and 35 cm depth beneath the East River, Colorado, from July–October 2017 (relatively normal water year) and April to October 2018 (comparatively low flow year), enabled by distinctive, in-situ, high frequency (Δt = 5 min) sensors that provided continuous time-series from the undisturbed study site over 14 months. We expect that hyporheic DO, which has a regular daily fluctuation pattern, is supplied by the surface water (at all times we estimate downwelling) and that diurnal hyporheic DO temporal patterns should be aligned with diurnal hyporheic temperature patterns. However, this was not found to be the case. Hyporheic DO becomes depleted briefly at 20 and 35 cm depths in 2017, and at all three hyporheic depths for extended periods in 2018. Whereas diurnal temperature fluctuations have consistent timings of maxima and minima, hyporheic DO rarely has as regular a pattern, and daily ranges are inconsistent. Rainfall events caused some of these changes to diurnal hyporheic DO patterns without repeatable effects. Antecedent snowpack conditions influence streamflow dynamics and therefore hyporheic DO dynamics in this alpine river. These results also point to the strong and variable influence of hyporheic microbial communities regulating hyporheic DO.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Investigation of Scale-Dependent Groundwater/Surface-water Exchange in Rivers by Gradient Self-Potential Logging: Numerical Modeling and Field Experiments

Exchanges of groundwater and surface-water are fundamental to a wide range of water-supply and water-quality management issues but challenging to map beyond the reach scale. Waterborne gradient self-potential (SP) measurements are directly sensitive to water flow through riverbed sediments and can be used to infer exchange locations, direction (gain versus loss), scale, and relative changes, but to date applications to river corridor hydrology are limited. Numerical modeling and field experiments were therefore performed herein, each emphasizing waterborne gradient SP logging for identifying and locating focused vertical groundwater discharge (surface-water gain) and recharge (surface-water loss) in a river. Two and three-dimensional numerical models were constructed to simulate the polarities, appearances, and peak amplitudes of streaming-potential and electric-field anomalies on a riverbed and in the surface-water that were attributable to steady-state vertical fluxes of groundwater through high-permeability conduits in the riverbed. Effects of varied hydraulic length-scale of exchange and surface-water depth were tested through numerical modeling. Modeling results aided in data acquisition and interpretation for three repeated field experiments performed along a 1.5–2.0 km reach of the Quashnet River in Cape Cod, Massachusetts, where focused, meter-scale groundwater discharges occur at discrete locations within otherwise ubiquitous and more diffuse groundwater upwelling conditions. Strong gradient SP anomalies were repeatedly measured in the Quashnet River at previously confirmed locations of focused groundwater discharge, showing the efficacy of waterborne gradient SP logging in identifying and characterizing groundwater/surface water exchange dynamics at multiple river network scales.

Engineering↗
DOE OSTI2021

Can Simple Machine Learning Tools Extend and Improve Temperature-Based Methods to Infer Streambed Flux?

Temperature-based methods have been developed to infer 1D vertical exchange flux between a stream and the subsurface. Current analyses rely on fitting physically based analytical and numerical models to temperature time series measured at multiple depths to infer daily average flux. These methods have seen wide use in hydrologic science despite strong simplifying assumptions including a lack of consideration of model structural error or the impacts of multidimensional flow or the impacts of transient streambed hydraulic properties. We performed a “perfect-model experiment” investigation to examine whether regression trees, with and without gradient boosting, can extract sufficient information from model-generated subsurface temperature time series, with and without added measurement error, to infer the corresponding exchange flux time series at the streambed surface. Using model-generated, synthetic data allowed us to assess the basic limitations to the use of machine learning; further examination of real data is only warranted if the method can be shown to perform well under these ideal conditions. We also examined whether the inherent feature importance analyses of tree-based machine learning methods can be used to optimize monitoring networks for exchange flux inference.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1998

Development of a Coupled Hydrological/Sediment Yield Model for a Watershed at Regional Level

Development of a hydrologic model for the study of environmental conservation requires a comprehensive understanding of individual-storm affecting hydrologic and sedimentologic processes. The hydrologic models that we are currently coupling are the Simulator for Hydrology and Energy Exchange at the Land Surface (SHEELS) and the Distributed Runoff Model (DRUM). SHEELS runs continuously to estimate surface energy fluxes and sub-surface soil water fluxes, while DRUM operates during and following precipitation events to predict surface runoff and peak flow through channel routing. The lateral re-distribution of surface water determined by DRUM is passed to SHEELS, which then adjusts soil water contents throughout the profile. The model SHEELS is well documented in Smith et al. (1993) and Laymen and Crosson (1995). The model DRUM is well documented in Vieux et al. (1990) and Vieux and Gauer (1994). The coupled hydrologic model, SHEELS/DRUM, does not simulate sedimentologic processes. The simulation of the sedimentologic process is important for environmental conservation planning and management. Therefore, we attempted to develop a conceptual frame work for coupling a sediment yield model with SHEELS/DRUM to estimate individual-storm sediment yield from a watershed at a regional level. The sediment yield model that will be used for this study is the Universal Soil Loss Equation (USLE) with some modifications to enable the model to predict individual-storm sediment yield. The predicted sediment yield does not include wind erosion and erosion caused by irrigation and snow melt. Units used for this study are those given by Foster et al. (1981) for SI units.

Rajbhandaril, Narayan↗
DOE OSTI2020

A multirate mass transfer model to represent the interaction of multicomponent biogeochemical processes between surface water and hyporheic zones (SWAT-MRMT-R 1.0)

Abstract. Surface water quality along river corridors can be modulated by hyporheic zones (HZs) that are ubiquitous and biogeochemically active. Watershed management practices often ignore the potentially important role of HZs as a natural reactor. To investigate the effect of hydrological exchange and biogeochemical processes on the fate of nutrients in surface water and HZs, a novel model, SWAT-MRMT-R, was developed coupling the Soil and Water Assessment Tool (SWAT) watershed model and the reaction module from a flow and reactive transport code (PFLOTRAN). SWAT-MRMT-R simulates concurrent nonlinear multicomponent biogeochemical reactions in both the channel water and its surrounding HZs, connecting the channel water and HZs through hyporheic exchanges using multirate mass transfer (MRMT) representation. Within the model, HZs are conceptualized as transient storage zones with distinguished exchange rates and residence times. The biogeochemical processes within HZs are different from those in the channel water. Hyporheic exchanges are modeled as multiple first-order mass transfers between the channel water and HZs. As a numerical example, SWAT-MRMT-R is applied to the Hanford Reach of the Columbia River, a large river in the United States, focusing on nitrate dynamics in the channel water. Major nitrate contaminants entering the Hanford Reach include those from the legacy waste, irrigation return flows (irrigation water that is not consumed by crops and runs off as point sources to the stream), and groundwater seepage resulting from irrigated agriculture. A two-step reaction sequence for denitrification and an aerobic respiration reaction is assumed to represent the biogeochemical transformations taking place within the HZs. The spatially variable hyporheic exchange rates and residence times in this example are estimated with the basin-scale Networks with EXchange and Subsurface Storage (NEXSS) model. Our simulation results show that (1), given a residence time distribution, how the exchange fluxes to HZs are approximated when using MRMT can significantly change the amount of nitrate consumption in HZs through denitrification and (2) source locations of nitrate have a different impact on surface water quality due to the spatially variable hyporheic exchanges.

58 GEOSCIENCES↗
DOE OSTI2024

Sinuosity‐Driven Hyporheic Exchange: Hydrodynamics and Biogeochemical Potentials

Abstract Hydrologic exchange processes are critical for ecosystem services along river corridors. Meandering contributes to this exchange by driving channel water, solutes, and energy through the surrounding alluvium, a process called sinuosity‐driven hyporheic exchange. This exchange is embedded within and modulated by the regional groundwater flow (RGF), which compresses the hyporheic zone and potentially diminishes its overall impact. Quantifying the role of sinuosity‐driven hyporheic exchange at the reach‐to‐watershed scale requires a mechanistic understanding of the interplay between drivers (meander planform) and modulators (RGF) and its implications for biogeochemical transformations. Here, we use a 2D, vertically integrated numerical model for flow, transport, and reaction to analyze sinuosity‐driven hyporheic exchange systematically. Using this model, we propose a dimensionless framework to explore the role of meander planform and RGF in hydrodynamics and how they constrain nitrogen cycling. Our results highlight the importance of meander topology for water flow and age. We demonstrate how the meander neck induces a shielding effect that protects the hyporheic zone against RGF, imposing a physical constraint on biogeochemical transformations. Furthermore, we explore the conditions when a meander acts as a net nitrogen source or sink. This transition in the net biogeochemical potential is described by a handful of dimensionless physical and biogeochemical parameters that can be measured or constrained from literature and remote sensing. This work provides a new physically based model that quantifies sinuosity‐driven hyporheic exchange and biogeochemical reactions, a critical step toward their representation in water quality models and the design and assessment of river restoration strategies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Simplified Method for the In Situ Collection and Laboratory Analysis of Cosmogenic Tracers (Sulfur-35 and Sodium-22) to Determine Residence Time Distributions and Water Ages

The use of cosmogenically produced sulfur-35 (T 1/2 = 87 days) and sodium-22 (T 1/2 = 2.6 years) as intrinsic tracers can provide valuable information on catchment hydrology, flow paths, and subsurface storage. A new and straightforward method was created to determine the activities of both 35 S and 22 Na in various water sources by pumping large volumes (up to 1000 L) of water through cation- and anion-exchange resin columns in the field to collect sodium and sulfate ions and simple chemistry in the lab. Samples are counted for 35 S using liquid scintillation counting (LSC) and for 22 Na via γ spectroscopy. Our novel in situ method provides faster sample throughput as well as better counting statistics and lower detection limits. Both methods were successfully applied at the Southern Sierra Critical Zone Observatory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Hyporheic hydraulic geometry: Conceptualizing relationships among hyporheic exchange, storage, and water age

Hyporheic exchange is now widely acknowledged as a key driver of ecosystem processes in many streams. Yet stream ecologists have been slow to adopt nuanced hydrologic frameworks developed and applied by engineers and hydrologists to describe the relationship between water storage, water age, and water balance in finite hydrosystems such as hyporheic zones. Here, in the context of hyporheic hydrology, we summarize a well-established mathematical framework useful for describing hyporheic hydrology, while also applying the framework heuristically to visualize the relationships between water age, rates of hyporheic exchange, and water volume within hyporheic zones. Building on this heuristic application, we discuss how improved accuracy in the conceptualization of hyporheic exchange can yield a deeper understanding of the role of the hyporheic zone in stream ecosystems. Although the equations presented here have been well-described for decades, our aim is to make the mathematical basis as accessible as possible and to encourage broader understanding among aquatic ecologists of the implications of tailed age distributions commonly observed in water discharged from and stored within hyporheic zones. Our quantitative description of “hyporheic hydraulic geometry,” associated visualizations, and discussion offer a nuanced and realistic understanding of hyporheic hydrology to aid in considering hyporheic exchange in the context of river and stream ecosystem science and management.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Modeling framework for evaluating the impacts of hydrodynamic pressure on hydrologic exchange fluxes and residence time for a large-scale river section over a long-term period

Quantifying hydrologic exchange fluxes (HEF) at the river and subsurface interface and their residence times (RT) in subsurface are important for managing the water quality and ecosystem health in dynamic river corridor systems. In this study, a modeling framework is developed for coupling the three-dimensional (3D) multi-phase surface, subsurface flow transport, and numerical tracer model for RT in a large-scale river section over a long period. The framework is utilized to evaluate the impacts of hydrodynamic pressure on HEFs and RT for a 30 km section of the Columbia River in Washington State over a three-year period. Based on comparisons between model simulations with and without considering hydrodynamic pressure, we found that hydrodynamic pressure increases the net HEFs by 7% with river gaining water from the subsurface domain, and leads to slight reduction of RT.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Analysis of Hydrology and Interim Measure Performance in the Chromium Plume at Los Alamos National Laboratory - 20498

Hexavalent chromium, Cr(VI), is present in the regional aquifer beneath Los Alamos National Laboratory (LANL) at concentrations greater than the New Mexico groundwater standard (50 ppb). The source of the Cr is blowdown from a power plant where potassium dichromate was used as a corrosion inhibitor from 1956 to 1972. Blowdown comingled with effluent from another source and migrated approximately 3 miles down a canyon, primarily as surface flow, until reaching an alluvial infiltration zone. The infiltrating water moved vertically through a thick (approximately 900 ft), complex stratigraphic sequence of unsaturated tuffs, basalt, and unconsolidated alluvial deposits. Horizontal transport also took place within spatially limited perched groundwater horizons in the vadose zone (VZ). Breakthrough at the regional aquifer water table is thought to have occurred in multiple hydraulic windows or 'drip points' located approximately 3 miles from the initial release site at the power plant. The regional aquifer is highly heterogeneous and anisotropic, and displays unconfined behavior in an upper zone and confined behavior at depth. The Cr(VI) contamination has so far remained confined to the upper approximately 25 m of the aquifer. The ambient hydraulic gradient in the regional aquifer is toward the east/southeast, but the water table is extremely flat in the Cr plume area compared to surrounding areas. There is some local variability in flow direction, potentially related to effects of seasonal pumping from nearby Los Alamos county water-supply wells, which draw from the regional aquifer several hundred feet below the water table. An interim measure (IM) is currently operating at the site to achieve hydraulic control of plume migration using a series of extraction wells, ion-exchange treatment, and injection along the downgradient portion of the plume. The migration of Cr contamination in the aquifer has been previously evaluated using a calibrated probabilistic numerical model, which has also been used to design the spatial configuration of extraction and injection wells for the IM, and to identify locations for monitoring wells. The model is built using the Finite Element Heat and Mass (FEHM) transfer code, and the Model Analysis and Decision Support (MADS) software is used for analyses. To date, the portion of the IM system that has operated in a continuous manner has been very effective at hydraulic control. Modifications to the flow field associated with the influence of IM pumping and injection are evaluated here, with modeling constraints to be provided by recent observations of tracers injected into both monitoring wells and IM injection wells, as well as opportunistic observations of geochemical signatures of treated ion-exchange effluent injected into IM injection wells. The tracer and geochemical data provide insights into modified flow directions within the Cr plume area that cannot be otherwise obtained. Results of the modeling work are being used to inform adaptive management of IM operations and will be a valuable tool for development and refinement of long-term remediation strategies. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1992

Emissions of biogenic sulfur gases from northern bogs and fens

Sulfur gases are important components of the global cycle of S. They contribute to the acidity of precipitation and they influence global radiation balance and climate. The role of terrestrial sources of biogenic S and their effect on atmospheric chemistry remain as major unanswered questions in our understanding of the natural S cycle. The role of northern wetlands as sources and sinks of gaseous S by measuring rates of S gas exchange as a function of season, hydrologic conditions, and gradients in tropic status was investigated. Experiments were conducted in wetlands in New Hampshire, particularly a poor fen, and in Mire 239, a poor fen at the Experimental Lakes Area (ELA) in Ontario. Emissions were determined using Teflon enclosures, gas cryotrapping methods and gas chromatography (GC) with flame photometric detection. Dynamic (sweep flow) and static enclosures were employed which yielded similar results. Dissolved S gases and methane were determined by gas stripping followed by GC.

Demello, William Zamboni↗
DOE OSTI2020

Modeling Functional Organic Chemistry in Arctic Rivers: An Idealized Siberian System

Rivers of the Arctic will become ever more important for the global climate, since they carry a majority of continental dissolved organic carbon flux into the rapidly changing polar ocean. Aqueous organics comprise a wide array of functional groups, several of which are likely to impact coastal and open water biophysical properties. Light attenuation, interfacial films, aerosol formation, gas release and momentum exchange can all be cited. We performed Lagrangian kinetic modeling for the evolution of riverine organic chemistry as the molecules in question make their way from the highlands to Arctic outlets. Classes as diverse as the proteins, sugars, lipids, re-condensates, humics, bio-tracers and small volatiles are all included. Our reduced framework constitutes an idealized northward flow driving a major hydrological discharge rate and primarily representing the Russian Lena. Mountainous, high solute and tundra sources are all simulated, and they meet up at several points between soil and delta process reactors. Turnover rates are parameterized beginning with extrapolated coastal values imposed along a limited tributary network, with connections between different terrestrial sub-ecologies. Temporal variation of our total dissolved matter most closely resembles the observations when we focus on the restricted removal and low initial carbon loads, suggesting relatively slow transformation along the water course. Thus, channel combinations and mixing must play a dominant role. Nevertheless, microbial and photochemical losses help determine the final concentrations for most species. Chemical evolution is distinct for the various functionalities, with special contributions from pre- and post-reactivity in soil and delta waters. Several functions are combined linearly to represent the collective chromophoric dissolved matter, characterized here by its absorption. Tributaries carry the signature of lignin phenols to segregate tundra versus taiga sources, and special attention is paid to the early then marine behaviors of low molecular weight volatiles. Heteropolycondensates comprise the largest percentage of reactive carbon in our simulations due to recombination/accumulation, and they tend to be preeminent at the mouth. Outlet concentrations of individual structures such as amino acids and absorbers lie above threshold values for biophysical influence, on the monolayer and light attenuation. The extent of coastal spreading is examined through targeted regional box modeling, relying on salinity and color for calibration. In some cases, plumes reach the scale of peripheral arctic seas, and amplification is expected during upcoming decades. Conclusions are mapped from the Lena to other boreal discharges, and future research questions are outlined regarding the bonding type versus mass release as permafrost degrades. Dynamic aqueous organic coupling is recommended for polar system models, from headwaters to coastal diluent.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

2D reactive transport model of shale chemical weathering and biogeochemical fluxes along a mountainous hillslope, East River Watershed, Colorado: Input files and simulation results

This data package contains input files and simulation results for a two-dimensional (2D) reactive transport model used to quantitatively analyze the coupled hydrological and biogeochemical processes governing shale weathering and associated biogeochemical fluxes under realistic environmental conditions in the high-elevation East River Watershed. These data support the conclusions presented in Stolze et al. (Water Resources Research, under review), "Model-based interpretation of solute exports and carbon partitioning during shale weathering in a mountainous hillslope". The model simulates atmospheric-subsurface gas exchange, subsurface water flow, and shale weathering processes under dynamic, year-scale conditions along a shale-underlain hillslope located in the East River watershed. The simulations were performed using the PFLOTRAN flow and reactive transport code and executed on the Perlmutter supercomputer to leverage its large-scale parallel computing capabilities. The data package contains two zipped folders, "model_input_files" and "simulation_results", and one readme.txt file. "model_input_files" contains the necessary input files to run the calibrated base-base model presented in Stolze et al. (Water Resources Research, under review). "simulation_results" contains a single hdf5 file ("Output_2D_hillslope_model.h5") which includes the results of simulation performed using the base-case model. This file can be opened with HDFView 3.1.4, Python, or MATLAB. "readme.txt" contains relevant information about the base-case model and provides guidelines on how to run the associated input files provided in the folder "model_input_files". Furthermore, readme.txt provides information regarding the model results provided in "Output_2D_hillslope_model.h5" such as matrix dimensionality and output units. Field datasets used to evaluate model performance were collected at three monitoring wells located along a hillslope transect (PLM1, PLM2, and PLM3). Dissolved ion concentration data were collected from November 2016 to October 2021 for Ca, Mg, DIC, Na, K, SO4 (Dong et al., 2025 - dic_npoc_data_2014_2024.zip - DOI:10.15485/1660459; Williams et al., 2025 - anion_data_2014_2024.zip - DOI:10.15485/1668054; Dong et al., 2025 - cation_data_2014_2024.zip - DOI:10.15485/1668055). Note that we used the files named er_PLM1_xx_yy, er_PLM2_xx_yy, and er_PLM3_xx_yy where xx stands for the name of the aqueous species and yy stands for the depth where the measurements were performed. Soil water content ([0 - 1] m) and water table depth were collected from November 2016 to October 2021 (Wan et al., 2024 - Dynamic_water_table__depthsFig2b.csv and Soil_water_content_Fig4e.csv - DOI:10.15485/2322567). Gaseous CO2 concentration were collected from October 2020 to December 2021(Wan et al., 2024 - Soil_CO2_concentrations_Fig4h.csv - DOI:10.15485/2322567) Gaseous CO2 flux from the subsurface to the atmosphere were collected in the vicinity of PLM2 from October 2019 to May 2022 (Wu et al., 2025). Soil microbial biomass concentration was measured from August 2016 to June 2017 (Sorensen et al., 2019 - 2017_East_River_Pumphouse_Microbial_Biomass__1_.csv - DOI:10.15485/1577267) All field data are published as CSV files compatible with Microsoft Excel, MATLAB, and Python, or as text files. The coordinates of the monitoring wells and the CO2(g) flux sensor in the coordinate system WGS84 are: -PLM1: [38.9197710 ; -106.9492750] -PLM2: [38.9201580 ; -106.9487170] -PLM3: [38.9207843 ; -106.9483668] -PLM4: 38.9210060 ; -106.9479528] -CO2(g) flux sensor: [38.9199180 ; -106.9489906] ------------------------------------------------------------------------------------------- This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231. This research used resources of the National Energy Research Scientific Computing Center (NERSC), a Department of Energy User Facility using NERSC award BER-ERCAP 23980, BER-ERCAP 28550, and BER-ERCAP 33789.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Technical note: Using long short-term memory models to fill data gaps in hydrological monitoring networks

Abstract. Quantifying the spatiotemporal dynamics in subsurface hydrological flows over a long time window usually employs a network of monitoring wells. However, such observations are often spatially sparse with potential temporal gaps due to poor quality or instrument failure. In this study, we explore the ability of recurrent neural networks to fill gaps in a spatially distributed time-series dataset. We use a well network that monitors the dynamic and heterogeneous hydrologic exchanges between the Columbia River and its adjacent groundwater aquifer at the U.S. Department of Energy's Hanford site. This 10-year-long dataset contains hourly temperature, specific conductance, and groundwater table elevation measurements from 42 wells with gaps of various lengths. We employ a long short-term memory (LSTM) model to capture the temporal variations in the observed system behaviors needed for gap filling. The performance of the LSTM-based gap-filling method was evaluated against a traditional autoregressive integrated moving average (ARIMA) method in terms of error statistics and accuracy in capturing the temporal patterns of river corridor wells with various dynamics signatures. Our study demonstrates that the ARIMA models yield better average error statistics, although they tend to have larger errors during time windows with abrupt changes or high-frequency (daily and subdaily) variations. The LSTM-based models excel in capturing both high-frequency and low-frequency (monthly and seasonal) dynamics. However, the inclusion of high-frequency fluctuations may also lead to overly dynamic predictions in time windows that lack such fluctuations. The LSTM can take advantage of the spatial information from neighboring wells to improve the gap-filling accuracy, especially for long gaps in system states that vary at subdaily scales. While LSTM models require substantial training data and have limited extrapolation power beyond the conditions represented in the training data, they afford great flexibility to account for the spatial correlations, temporal correlations, and nonlinearity in data without a priori assumptions. Thus, LSTMs provide effective alternatives to fill in data gaps in spatially distributed time-series observations characterized by multiple dominant frequencies of variability, which are essential for advancing our understanding of dynamic complex systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Interactions between molecular-scale processes and hyporheic exchange for understanding Fe-S-C cycling in riparian wetlands (Final Report)

Wetlands represent some of the most productive ecosystems on the planet and critically influence global environmental health. Specifically, wetlands promote water quality by transforming nutrients and organic compounds and sequestering metals and contaminants. Riparian wetland hyporheic zones, where toxic surface water and anoxic groundwater mix, exhibit dynamic conditions that drive steep redox gradients and promote hotspots of diverse and fluctuating microbial activity. Changes in climate, water quality, and water quantity can disturb hydrologic flow and biogeochemical processing in these environments. Understanding how sulfate loading, impacted by hydrologic flux and anthropogenic inputs, influences iron and carbon cycling in wetlands will be crucial for predicting water quality issues driven by iron mineral precipitation and sorption, such as the release of heavy metals and other toxic elements. We used a fully integrated multi-scale and multi-method approach to develop a mechanistic understanding of how hydrologic flow influences coupled iron and sulfur cycles in riparian wetlands. This entailed hydrological, geochemical, and microbial observations at two locations: an anthropogenic sulfate-impacted riparian wetland in northern Minnesota, and a low-sulfate Fe-rich riparian wetland in Tims Branch at the Savannah River Site (SRS). These sites were characterized by hydrologically dynamic conditions where the stream and wetland systems oscillated between gaining (upward flow) and losing (downward flow) conditions that recharged the system occasionally with oxidants that fueled a variety of biogeochemical reactions. Aqueous geochemical measurements of surface water, groundwater, and porewater samples were made alongside solid-phase geochemical analyses of sediment gravity cores. Bulk X-ray absorption spectroscopy at the Advanced Photon Source (APS; Argonne) interrogated the speciation and distribution of Fe and S mineral phases of the sediments. Interestingly, it was discovered that, despite strongly reducing conditions, Fe(III) compounds and a variety of intermediate valence S compounds were stable in the subsurface. Indeed, compounds like thiosulfate, S(0), and intermediate valence organosulfur compounds were more prevalent than FeS and pyrite. The composition of the sediments did change with changes in hydrologic flow, showing their reactivity in changing redox conditions. The abundance of these intermediate S compounds was likely formed as a result of anaerobic oxidation by aqueous and solid-phase Fe(III) compounds, fueling a cryptic S cycle that is driving the breakdown of organic matter in the hyporheic zone. Microbiome surveys showed a core community that seemed stable across the landscape, but changed with increasing depth into the sediment. The community composition did not seem to change dramatically with changes in season or hydrologic flow, except for some organisms that have the potential to contribute to S cycling. More work is needed to confirm their functional activity. These fine process-scale analyses were placed within a dynamic field context using physical flow parameters from surface water and groundwater level measurements. These data helped shed light on sulfur-driven biogeochemical processes in hydrologically dynamic riparian wetlands, addressing a gap in our understanding about the impacts of pollution and other anthropogenic changes on ecologically sensitive environments.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1992

Sulfur gas exchange in Sphagnum-dominated wetlands

Sulfur gases are important components of the global cycle of S. They contribute to the acidity of precipitation and they influence global radiation balance and climate. The role of terrestrial sources of biogenic S and their effect on atmospheric chemistry remain as major unanswered questions in our understanding of the natural S cycle. The role of northern wetlands as sources and sinks of gaseous S was investigated by measuring rates of S gas exchange as a function of season, hydrologic conditions, and gradients in trophic status. The effects of inorganic S input on the production and emission of gaseous S were also investigated. Experiments were conducted in wetlands in New Hampshire, particularly a poor fen, fens within the Experimental Lakes Area (ELA) in Ontario, Canada and in freshwater and marine tundra. Emissions were determined using Teflon enclosures, gas cryotrapping methods, and gas chromatography (GC) with flame photometric detection. Dynamic (sweep flow) and static enclosures were employed. Dissolved gases were determined by gas stripping followed by GC.

Hines, Mark E.↗
DOE OSTI2025

Evaluating the Influence of Plants on Hydrologic Cycling: Quantifying and Validating the Role of Plant Processes and Stomatal Conductance (Final Report)

Plants can exert strong controls on water cycling, both locally and across the globe. The impacts are driven by how plants control the exchange of water and energy between the land and the atmosphere, which in turn are influenced by how plants photosynthesize and grow (e.g. biogeochemical cycling). Our team has shown previously that plant responses to increasing CO 2 in the atmosphere can influence the water cycle through changes in rainfall and the amount of water on land with implications for drought. Plant responses can also alter both average and extreme flow of water in rivers, with implications for freshwater availability and flood frequency. Although the impact of plant responses on water cycling has been demonstrated to exist, significant uncertainty remains in our understanding of the magnitude and form of the plant responses themselves, as well as their ultimate impact on water availability for people and ecosystems. In this project we worked to quantify the role of plant processes in regulating water cycling (e.g., precipitation, evapotranspiration, runoff, droughts, etc.). We addressed the following objectives: quantify the control of physiological vs. radiative effects on water cycling across many coupled models in idealized and realistic scenarios, as well as to quantify how our assumptions about leaf-level processes (e.g. coupling between stomatal conductance and photosynthesis), and organism-level processes (e.g. leaf area response to high CO 2 ) contributed to uncertainties and biases in water cycling. We focused on the impact of these processes on water cycling to assess the sensitivities of precipitation, evapotranspiration, and runoff to model assumptions about plant processes at the leaf, organism, and community level. As part of this project we used measurements of stable carbon isotopes to generate observationally based estimates of plant functioning during the historical period. This work is directly relevant to DOE’s primary scientific research questions in the RGMA topic area (b) on Biogeochemical processes, feedbacks, and interactions in the Earth system.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Coupling Surface Flow with High-performance Subsurface Reactive Flow and Transport Code PFLOTRAN

Water exchange between the surface and subsurface is important for both water resource management and environmental protection. In this paper, we develop coupled surface and subsurface flow simulation capability in a parallel subsurface flow and reactive transport code PFLOTRAN. We sequentially couple the diffusion wave-based surface flow with the subsurface flow governedby the Richards equation in PFLOTRAN. These two flow domains are linked with a boundary condition switching method that ensures continuity of pressure and flux at the surface-subsurface interface. We verify the coupled code against other existing hydrologic models and observation data using a number of numerical experiments. The coupled hydrological model exhibits good performance in strong parallel scaling tests. The new coupled surface and subsurface simulator significantly advance community simulation capability towards improving integrated hydrologic and biogeochemical understanding of complex systems such as watersheds and river corridors. Keywords: Surface flow, Integrated hydrological modeling, Boundary condition switching, Parallel computing

Wu, Runjian↗