Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “biogeochemistry”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Perspective on identifying and characterizing the processes controlling iron speciation and residence time at the atmosphere-ocean interface

It is well recognized that the atmospheric deposition of iron (Fe) affects ocean productivity, atmospheric CO 2 uptake, ecosystem diversity, and overall climate. Despite significant advances in measurement techniques and modeling efforts, discrepancies persist between observations and models that hinder accurate predictions of processes and their global effects. Here, we provide an assessment report on where the current state of knowledge is and where future research emphasis would have the highest impact in furthering the field of Fe atmosphere-ocean biogeochemical cycle. These results were determined through consensus reached by diverse researchers from the oceanographic and atmospheric science communities with backgrounds in laboratory and in situ measurements, modeling, and remote sensing. We discuss i) novel measurement methodologies and instrumentation that allow detection and speciation of different forms and oxidation states of Fe in deliquesced mineral aerosol, cloud/rainwater, and seawater; ii) oceanic models that treat Fe cycling with several external sources and sinks, dissolved, colloidal, particulate, inorganic, and organic ligand-complexed forms of Fe, as well as Fe in detritus and phytoplankton; and iii) atmospheric models that consider natural and anthropogenic sources of Fe, mobilization of Fe in mineral aerosols due to the dissolution of Fe-oxides and Fe-substituted aluminosilicates through proton-promoted, organic ligand-promoted, and photo-reductive mechanisms. In addition, the study identifies existing challenges and disconnects (both fundamental and methodological) such as i) inconsistencies in Fe nomenclature and the definition of bioavailable Fe between oceanic and atmospheric disciplines, and ii) the lack of characterization of the processes controlling Fe speciation and residence time at the atmosphere-ocean interface. Such challenges are undoubtedly caused by extremely low concentrations, short lifetime, and the myriad of physical, (photo)chemical, and biological processes affecting global biogeochemical cycling of Fe. However, we also argue that the historical division (separate treatment of Fe biogeochemistry in oceanic and atmospheric disciplines) and the classical funding structures (that often create obstacles for transdisciplinary collaboration) are also hampering the advancement of knowledge in the field. Finally, the study provides some specific ideas and guidelines for laboratory studies, field measurements, and modeling research required for improved characterization of global biogeochemical cycling of Fe in relationship with other trace elements and essential nutrients. In conclusion, the report is intended to aid scientists in their work related to Fe biogeochemistry as well as program managers at the relevant funding agencies

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

Diurnal rainfall response to the physiological and radiative effects of CO 2 in tropical forests in the Energy exascale earth system model v1

Understanding how the connection between rainfall and tropical forests will respond to increasing CO$_2$ concentrations is a key element in understanding how the tropical water cycle will respond to increasing CO$_2$. The plant physiological and radiative impacts of CO$_2$ on rainfall patterns over tropical forest regions are examined in the Energy Exascale Earth System Model version 1.1 (E3SMv1.1-BGC) biogeochemistry experiments. Composite analysis reveals a dampening of the diurnal cycle of rainfall over the Amazon, Congo, and Maritime Continent in response to rising CO$_2$ levels, regardless of the sign of total rainfall change. A full factorial model experiment confirms that the CO$_2$ radiative and CO$_2$ plant physiological effects can individually or jointly reduce the magnitude of the rainfall diurnal cycle, though the physical pathway giving rise to the reduction differs between the two effects. For the physiological response, stomatal closure reduces evapotranspiration, which dries the boundary layer and raises the lifting condensation level. These effects combine to reduce deep convective rainfall during its peak occurrence in the late daytime to early nighttime period. For the radiative response, a relative reduction in daytime CAPE (consistent with a reduction in the diurnal temperature range) leads to less frequent triggering of deep convection and a reduction of rainfall diurnal amplitude. These diurnal rainfall changes are structurally similar across seasons, and show little sensitivity to representation of nutrient coupling for the land biogeochemistry. In agreement with previous findings, the physiological response has only minor impact on extreme rainfall relative to the radiative response.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Conceptualizing Biogeochemical Reactions With an Ohm's Law Analogy

Abstract In studying problems like plant‐soil‐microbe interactions in environmental biogeochemistry and ecology, one usually has to quantify and model how substrates control the growth of, and interaction among, biological organisms (and abiotic factors, e.g., adsorptive mineral soil surfaces). To address these substrate‐consumer relationships, many substrate kinetics and growth rules have been developed, including the famous Monod kinetics for single‐substrate‐based growth and Liebig's law of the minimum for multiple‐nutrient‐colimited growth. However, the mechanistic basis that leads to these various concepts and mathematical formulations and the implications of their parameters are often quite uncertain. Here, we show that an analogy based on Ohm's law in electric circuit theory is able to unify many of these different concepts and mathematical formulations. In this Ohm's law analogy, a resistor is defined by a combination of consumers’ and substrates’ kinetic traits. In particular, the resistance is equal to the mean first passage time that has been used to derive the Michaelis‐Menten kinetics under substrate replete conditions for a single substrate as well as the predation rate of individual organisms. We further show that this analogy leads to important insights on various biogeochemical problems, such as (a) multiple‐nutrient‐colimited biological growth, (b) denitrification, (c) fermentation under aerobic conditions, (d) metabolic temperature sensitivity, and (e) the legitimacy of Monod kinetics for describing bacterial growth. We expect that our approach will help both modelers and nonmodelers to better understand and formulate hypotheses when studying certain aspects of environmental biogeochemistry and ecology.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Intraspecific variability in plant and soil chemical properties in a common garden plantation of the energy crop Populus

Optimizing crops for synergistic soil carbon (C) sequestration can enhance CO 2 removal in food and bioenergy production systems. Yet, in bioenergy systems, we lack an understanding of how intraspecies variation in plant traits correlates with variation in soil biogeochemistry. This knowledge gap is exacerbated by both the heterogeneity and difficulty of measuring belowground traits. Here, we provide initial observations of C and nutrients in soil and root and stem tissues from a common garden field site of diverse, natural variant, Populus trichocarpa genotypes—established for aboveground biomass-to-biofuels research. Our goal was to explore the value of such field sites for evaluating genotype-specific effects on soil C, which ultimately informs the potential for optimizing bioenergy systems for both aboveground productivity and belowground C storage. To do this, we investigated variation in chemical traits at the scale of individual trees and genotypes and we explored correlations among stem, root, and soil samples. We observed substantial variation in soil chemical properties at the scale of individual trees and specific genotypes. While correlations among elements were observed both within and among sample types (soil, stem, root), above-belowground correlations were generally poor. We did not observe genotype-specific patterns in soil C in the top 10 cm, but we did observe genotype associations with soil acid-base chemistry (soil pH and base cations) and bulk density. Finally, a specific phenotype of interest (high vs low lignin) was unrelated to soil biogeochemistry. Our pilot study supports the usefulness of decade-old, genetically-variable, Populus bioenergy field test plots for understanding plant genotype effects on soil properties. Finally, this study contributes to the advancement of sampling methods and baseline data for Populus systems in the Pacific Northwest, USA. Further species- and region-specific efforts will enhance C predictability across scales in bioenergy systems and, ultimately, accelerate the identification of genotypes that optimize yield and carbon storage.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Dissolved oxygen sensor in an automated hyporheic sampling system reveals biogeochemical dynamics

Many river corridor systems frequently experience rapid variations in river stage height, hydraulic head gradients, and residence times. The integrated hydrology and biogeochemistry of such systems is challenging to study, particularly in their associated hyporheic zones. Here we present an automated system to facilitate 4-dimensional study of dynamic hyporheic zones. It is based on combining real-time in-situ and ex-situ measurements from sensor/sampling locations distributed in 3-dimensions. A novel dissolved oxygen (DO) sensor was integrated into the system during a small scale study. We measured several biogeochemical and hydrologic parameters at three subsurface depths in the riverbed of the Columbia River in Washington State, USA, a dynamic hydropeaked river corridor system. During the study, episodes of significant DO variations (~+/- 4 mg/l) were observed, with minor variation in other parameters (e.g., <~+/-0.15 mg/l NO 3 ). DO concentrations were related to hydraulic head gradients, showing both hysteretic and non-hysteretic relationships with abrupt (hours) transitions between the two types of relationships. The observed relationships provide a number of hypotheses related to the integrated hydrology and biogeochemistry of dynamic hyporheic zones. We suggest that preliminary high-frequency monitoring is advantageous in guiding the design of long term monitoring campaigns. The study also demonstrated the importance of measuring multiple parameters in parallel, where the DO sensor provided the key signal for identifying/detecting transient phenomena.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Seasonal changes in the drivers of water physico-chemistry variability of a small freshwater tidal river

Where rivers meet the sea, tides can exert a physical and chemical influence on the lower reaches of a river. How tidal dynamics in these tidal river reaches interact with upstream hydrological drivers such as storm rainfall, which ultimately determines the quantity and composition of material transferred from watersheds to estuaries, is currently unknown. We monitored a small freshwater tidal river in the Pacific Northwest, USA in high resolution over one year to evaluate the relative importance of tides versus upstream hydrological flows (i.e., base flow and precipitation events) on basic physico-chemical parameters (pH, dissolved oxygen, turbidity, specific conductivity, and temperature), and how these interactions relate to the downstream estuary. Tidal variability and diurnal cycles (i.e. solar radiation) dominated water physico-chemical variability in the summer, but the influence of these drivers was overshadowed by storm-driven sharp pulses in river physico-chemistry during the remainder of the year. Within such events, we found incidences of counterclockwise hysteresis of pH, counterclockwise hysteresis of dissolved oxygen, and clockwise hysteresis of turbidity, although systematic trends were not observed across events. The dominance of storm rainfall in the river’s physico-chemistry dynamics, and similar pulses of decreased pH observed in adjacent estuarine waters, suggest that the linkage between tidal streams and the broader system is variable throughout the year. High-frequency monitoring of tidal river biogeochemistry is therefore crucial to enable the assessment of how the relative strength of these drivers may change with future sea level rise and altered precipitation patterns to modulate biogeochemical dynamics across the land-ocean-atmosphere continuum.

aquatic↗
DOE OSTI2019

Spatial gradients in the characteristics of soil-carbon fractions are associated with abiotic features but not microbial communities

Coastal terrestrial-aquatic interfaces (TAIs) are dynamic zones of carbon biogeochemical cycling influenced by landscape-scale salinity gradients. However, there is significant cross-system heterogeneity in salinity influences over TAI soil biogeochemical function, such as gas fluxes and dissolved organic carbon concentrations. One potential explanation for the lack of consistency is unrecognized variation in soil organic carbon chemistry and microbial communities. To provide a foundation for understanding the role of C chemistry and microbial communities on biogeochemical cycling in coastal soils, we evaluated associations among a landscape scale soil salinity gradient, molecular-level soil carbon chemistry, and microbial community assembly processes. Based on observations in related systems we hypothesized systematic shifts along the salinity gradient in organic carbon thermodynamics, biochemical transformation profiles, and heteroatom content. While the thermodynamic favorability of water soluble organic compounds in shallow soils did decrease with increasing salinity (R2=0.78, p<0.00001), this was opposite the expected direction. The number of inferred biochemical transformations and total heteroatom content also decreased with increasing salinity. These patterns indicate lower microbial activity at higher salinity that is potentially constrained by less favorable organic carbon. Furthermore, organic compound classes appeared to be primarily marine/algal derived in the floodplains with greater terrestrial signatures in upland soil. Based on a recent simulation-based study we further hypothesized a relationship between microbial community assembly processes and carbon chemistry. Null modelling revealed particularly strong influences of dispersal limitation over microbial composition, which may be due limited hydrologic connectivity within the clay-rich soils. Signatures of dispersal limitation indicate stochastically assembled microbial communities, which was further reflected in the lack of an association between community assembly processes and carbon chemistry. This suggests a disconnect between microbial composition and carbon biogeochemistry, thereby indicating that the salinity-associated gradient in carbon chemistry was driven by a combination of spatially-structured inputs and salinity-associated metabolic responses of microbial communities that are independent of community composition. Our results collectively indicate that soil carbon chemistry varies systematically with salinity at the landscape scale. This variation has not been accounted for in studies examining impacts of salinity to biogeochemical rates and pools. We propose that impacts of salinity in coastal soil biogeochemistry need to be understood in context of carbon chemistry and microbial physiology, while microbial composition may have less influence.

soil carbon, Coastal Biogeochemistry, microbial co↗
DOE OSTI2025

Response of Subsurface Nitrogen-Cycling Microbial Communities to Environmental Fluctuations (Final Technical Report)

Riparian floodplains are dynamic ecosystems linking terrestrial and riverine systems. These floodplains experience hydrological shifts such as changes in water table height, flooding, and drought and can be ‘hotspots’ of biogeochemical cycling due to shifting sediment moisture (and saturation) and subsurface exchanges of water, nutrients, and other compounds across different sediment layers. Subsurface microbial communities are the primary drivers of biogeochemical processes in floodplains, and thus their structure and function can directly influence both surface and groundwater quality. The microbial nitrogen (N) cycle is particularly important in floodplains as it affects nutrient availability and removal. Two functional guilds of chemoautotrophic (i.e. CO2-fixing) microorganisms are responsible for the first oxidative step of the N cycle, nitrification: ammonia-oxidizing archaea (AOA) and bacteria (AOB) catalyze the oxidation of ammonia to nitrite, while nitrite-oxidizing bacteria (NOB) oxidize nitrite to nitrate. Despite the critical role nitrification plays in N-cycling in both terrestrial and aquatic ecosystems, our understanding of the diversity, ecophysiology, and activity of nitrifying organisms in subsurface floodplain soils/sediments is extremely limited. To help address this critical knowledge gap, the overarching goal of this project was to determine how shifts in key environmental parameters and gradients impact microbial N-cycling communities/processes, with particular emphasis on nitrification, within hydrologically-variable floodplain sediments in the Wind River Basin near Riverton, Wyoming. The three specific objectives of this project were to: (1) to associate in situ environmental drivers of N cycling with distinct functional guilds; (2) determine the guild response to variation in key ecosystem drivers; and (3) develop a dynamic ecosystem model of the microbial N cycle with the Riverton subsurface using community genomic and biogeochemical data collected in the first two objectives. Over the course of this project, we employed both 16S rRNA gene amplicon sequencing and genome-resolved metagenomics to examine the phylogenetic diversity and metabolic potential of subsurface nitrifier communities within 68 samples collected across multiple sites, depths, and time points within the Riverton floodplain, allowing for both spatial and temporal investigations at different scales. This project benefitted tremendously from recent advances in high-throughput sequencing technologies coupled with dramatic improvements in the computational tools and algorithms available for analyzing such large, complex genomic datasets. By pairing these cutting-edge genomic approaches with depth-resolved sampling and detailed geochemical analyses of the Riverton floodplain, we have gained novel insights into the structure and function of subsurface nitrifier communities in relation to both hydrology and biogeochemistry. This project resulted in the most detailed and comprehensive characterization of N-cycling floodplain microbial communities to date and will hopefully inspire and pave the way for future studies using similar approaches in other floodplains. Indeed, such information is critical for understanding subsurface biogeochemical cycling and how elemental stores are altered from perturbations initiated by the water cycle within floodplains. Finally, because of the terrestrial-aquatic nature of the Riverton floodplain, results from this project are also of relevance to disciplines such as soil science, estuarine science, limnology & oceanography, biogeochemistry, geobiology, environmental engineering, as well as genomics and data science.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Long-term hydrological, biogeochemical, and climatological data from Walker Branch Watershed, East Tennessee, USA

In 1967, the original Walker Branch Watershed (WBW) project was established to study elemental cycling and mass balances in a relatively unimpacted watershed. Over the next 50+ years, findings from additional experimental studies and long-term observations on WBW advanced understanding of catchment hydrology, biogeochemistry, and ecology and established WBW as a seminal site for catchment science. The 97.5-ha WBW is located in East Tennessee, USA, on the U.S. Department of Energy's Oak Ridge Reservation. Vegetation on the watershed is characteristic of an eastern deciduous, second-growth forest. The watershed is divided into two subcatchments: the West Fork (38.4 ha) and the East Fork (59.1 ha). Headwater streams draining these subcatchments are fed by multiple springs, and thus flow is perennial. Stream water is high in base cations due to weathering of dolomite bedrock and nutrient concentrations are low. Long-term observations of climate, hydrology, and biogeochemistry include daily (1969–2014) and 15-min (1994–2014) stream discharge and annual runoff (1969–2014); hourly, daily, and annual rainfall (1969–2012); daily climate and soil temperature (1993–2010); and weekly stream water chemistry (1989–2013). These long-term datasets are publicly available on the WBW website (https://walkerbranch.ornl.gov/long-term-data/). While collection of these data has ceased, related long-term measurements continue through the National Ecological Observatory Network (NEON), where WBW is the core terrestrial and aquatic site in the Appalachian and Cumberland Plateau region (NEON's Domain 7) of the United States. These long-term datasets have been and will continue to be important in evaluating the influence of climatic and environmental drivers on catchment processes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Validation of the Community Land Model Version 5 over the Contiguous United States (CONUS) using in situ and remote sensing data sets

The Community Land Model (CLM) is an effective tool to simulate the biophysical and biogeochemical processes and their interactions with the atmosphere. Although CLM Version 5 (CLM5) constitutes various updates in these processes, its performance in simulating energy, water and carbon cycles over the Contiguous United States (CONUS) at scales which land surface changes and hydrometeorological and hydroclimatological applications are more locally relevant is yet to be assessed. In this study, we conducted three simulations at 0.125? during 1979-2018 over the CONUS using different configurations of CLM, namely CLM5-biogeochemistry (CLM5BGC), CLM4.5BGC, and CLM5-satellite phenology (CLM5SP). We validated and compared their simulations against multiple remote-sensed and in-situ datasets. Overall, the parametric and structural updates (e.g., carbon cost for nitrogen uptake, variable soil thickness, dry surface layer) in CLM5 improve its ability in capturing terrestrial biogeochemical dynamics. The low evapotranspiration in CLM5BGC is associated with biases in simulating vegetation phenological characteristics rather than soil water limitations. The mismatch between CLM5BGC-simulated peak leaf area index and reference data can be attributed to CLM5BGC's inability in simulating phenology of trees and grasses. The differences between CLM-simulated irrigation and reference estimates can be attributed to differences between processes represented in models and in reality, and uncertainties in input and validation datasets. Evaluation against observations at small catchments suggest that hydrologic parameters needed to be calibrated to improve simulations of runoff, especially subsurface runoff. Additional efforts are needed to incorporate spatially-distributed plant phenology and physiology parameters and regional-specific agricultural management practices (e.g., planting, harvest).

Cheng, Yanyan↗
DOE OSTI2024

Machine Learning Driven Sensitivity Analysis of E3SM Land Model Parameters for Wetland Methane Emissions

Methane (CH 4 ) is globally the second most critical greenhouse gas after carbon dioxide, contributing to 16%–25% of the observed atmospheric warming. Wetlands are the primary natural source of methane emissions globally. However, wetland methane emission estimates from biogeochemistry models contain considerable uncertainty. One of the main sources of this uncertainty arises from the numerous uncertain model parameters within various physical, biological, and chemical processes that influence methane production, oxidation, and transport. Sensitivity Analysis (SA) can help identify critical parameters for methane emission and achieve reduced biases and uncertainties in future projections. This study performs SA for 19 selected parameters responsible for critical biogeochemical processes in the methane module of the Energy Exascale Earth System Model (E3SM) land model (ELM). The impact of these parameters on various CH 4 fluxes is examined at 14 FLUXNET- CH 4 sites with diverse vegetation types. Given the extensive number of model simulations needed for global variance-based SA, we employ a machine learning (ML) algorithm to emulate the complex behavior of ELM methane biogeochemistry. We found that parameters linked to CH 4 production and diffusion generally present the highest sensitivities despite apparent seasonal variation. Comparing simulated emissions from perturbed parameter sets against FLUXNET-CH 4 observations revealed that better performances can be achieved at each site compared to the default parameter values. This presents a scope for further improving simulated emissions using parameter calibration with advanced optimization techniques.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Transition Zones at the Changing Coastal Terrestrial‐Aquatic Interface

Coastal soils are a significant but highly uncertain component of global biogeochemical cycles. These systems experience spatial and temporal variability in biogeochemical processes, driven by marsh‐to‐upland gradients and hydrological fluctuations. These fluctuations make it difficult to understand and predict biogeochemical processes in these highly dynamic systems. We studied coastal soil biogeochemistry and its variability (a) at regional scales and (b) across transects from upland forest to marsh, in two contrasting regions—Lake Erie, a freshwater lacustrine system, and Chesapeake Bay, a saltwater estuarine system. Salinity‐related analytes were a key source of variability in soil biogeochemistry, not just in the saltwater system, but surprisingly, also in the freshwater system. We had hypothesized linear trends in biogeochemical parameters along the TAI—however, contrary to expectations, transition soils were not consistently intermediate between upland and marsh endmembers; the non‐monotonic trends of C, P, Fe along our transects suggest that these do not behave as expected and may be difficult to model and predict—thus these are key analytes to study in our regions. Rapidly changing soil factors across coastal gradients (e.g., Ca, K, CEC, and TS) may act as precursors to ecosystem shifts. Our comprehensive soil characterization represents a snapshot of a single timepoint of surface soils and provides essential data for mechanistic modeling of ecosystem dynamics across coastal transects.

Patel, Kaizad F. [Pacific Northwest National Labor↗
DOE OSTI2025

Implementation and Evaluation of Emission‐Driven Land‐Atmosphere Coupled Simulation in E3SMv2.1

Emissions-driven (prognostic CO 2 ) simulations are essential for representing two-way carbon-climate feedback in Earth System Models. We present an emissions-driven land–atmosphere coupled biogeochemistry (BGC) configuration (BGCLNDATM_progCO2) in version 2.1 of the Energy Exascale Earth System Model (E3SMv2.1). This is the first E3SM configuration that performs land-atmosphere emission-hindcasts. Here, we document its implementation, evaluate the model's performance against observations and other models, and propose a structured evaluation protocol for such emissions-driven simulations. We conducted transient historical simulations (1850–2014) with BGCLNDATM_progCO2 and compare them to reference simulations—a land-atmosphere coupled simulation without BGC and a standalone land simulation with BGC, both using prescribed CO 2 concentrations—and to observations. BGCLNDATM_progCO2 overestimates atmospheric CO 2 concentrations by 11–23 ppm yet stays within the 40-ppm spread CMIP6 emission-driven models and retains physical climate properties comparable to the reference runs. The CO 2 biases are partly attributed to underrepresented oceanic CO 2 uptake and inadequate representations of some terrestrial processes. In general, introducing prognostic CO 2 did not change physical climate metrics at the global scale but had larger regional effects, particularly over land where spatially heterogeneous CO 2 and prognostic leaf area index influenced surface energy balance. Finally, we propose a general evaluation protocol including spin-up assessment, atmospheric CO 2 benchmarking, physical climate evaluation, and land biogeochemical analysis to support scientific rigor and facilitate inter-model comparisons. The new configuration lays the groundwork for future enhancements, including improved terrestrial biogeochemical processes, integrated marine biogeochemistry, and additional human–Earth system interactions. These developments advance E3SM toward fully coupled emissions-driven simulations, enabling more accurate carbon–climate feedback projections and informing mitigation policy by providing physically consistent carbon-budget metrics for mitigation scenarios.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Hydrological and Thermal Dynamics of a Supra‐Permafrost Subterranean Estuary

Subterranean estuaries (STEs), where groundwater interacts with seawater, influence surface and subsurface coastal ecology and biogeochemistry. In Arctic-STEs overlying permafrost, groundwater flow and heat transport determine the fate of organic matter. Yet, direct observations of groundwater flow and heat and solute transport processes in Arctic STEs remain limited. This study characterized groundwater flow paths and fluxes and heat transport within an Arctic-STE along Alaska's Beaufort Sea coast during thawing, summer, and freeze-up. Intertidal seabed temperature-depth profiles collected along a 10-m transect captured the active groundwater flow period, from thaw and flow onset in mid-June to freeze-up in late-September. During this period, aquifer geometry evolved non-uniformly due to spatially varying thaw rates across the STE (mean (m) thaw depths–beach: 0.25 to 0.55–0.6 m on 20 June, 25 July, 1 October; seabed: 0.6–0.9 m from 25 July to 1 October). Groundwater and surface water levels, salinity, and subsurface temperature profiles measured over tidal time scales were interpreted alongside groundwater flow-heat transport numerical simulations. Fresh groundwater discharge was sporadic during thawing (m: 0.32 m 3 /day/m), abundant in summer (m: 0.45 m 3 /day/m), and was largely absent during freeze-up. During freeze-up, groundwater flow was driven exclusively by seawater recirculation via tidal pumping (from thawing to summer to freeze-up: 0.00025–0.15–0.5 m3/day/m) and convection. Heat advection dominated near aquatic interfaces (shaping intertidal ice), and conduction controlled vertical temperature gradients in low-flow and unsaturated sediments. These findings will help predict how prolonged summers will alter Arctic-STE cryo-hydrology and biogeochemistry.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Data for Soil Oxygen Dynamics: A Key Mediator of Tile Drainage Impacts on Coupled Hydrological, Biogeochemical, and Crop Systems

Tile drainage removes excess water and is an essential, widely adopted management practice to enhance crop productivity in the US Midwest and throughout the world. Tile drainage has been shown to significantly change hydrological and biogeochemical cycles by lowering the water table and reducing the residence time of soil water, although examining the complex interactions and feedbacks in an integrated hydrology–biogeochemistry–crop system remains elusive. Oxygen dynamics are critical to unraveling these interactions and have been ignored or oversimplified in existing models. Understanding these impacts is essential, particularly so because tile drainage has been highlighted as an adaptation under projected wetter springs and drier summers in the changing climate in the US Midwest. We used the ecosys model that uniquely incorporates first-principle soil oxygen dynamics and crop oxygen uptake mechanisms to quantify the impacts of tile drainage on hydrological and biogeochemical cycles and crop growth in corn–soybean rotation fields. The model was validated with data from a multi-treatment, multi-year experiment in Washington, IA. The relative root mean square error (rRMSE) for the corn and soybean yield in validation is 5.66 % and 12.57 %, respectively. The Pearson coefficient (r) of the monthly tile flow during the growing season is 0.78. Plant oxygen stress turns out as an emergent property of the equilibrium between the soil oxygen supply and biological demand. The impact of tile drainage on the system is achieved through a series of coupled feedback mechanisms. The model results show that tile drainage reduces the soil water content and enhances soil oxygenation. It additionally increases the subsurface discharge and elevates inorganic nitrogen leaching, with seasonal variations influenced by climate and crop phenology. The improved aerobic condition alleviates crop oxygen stress during wet springs, thereby promoting crop root growth during the early growth stage. The development of greater root density, in turn, mitigates water stress during dry summers, leading to an overall increase in the crop yield by ∼6 %. These functions indicate the potential of tile drainage in bolstering crop resilience to climate change and the use of this modeling tool for large-scale assessments of tile drainage. The model reveals the underlying causal mechanisms that drive the agroecosystem response to drainage on the coupled hydrology, biogeochemistry, and crop system dynamics.

Modeling↗
DOE OSTI2024

Untargeted metabolite data from a root surface in a rhizobox

Raw data is provided from samples analyzed using separate reverse phase chromatographic methods on a high performance liquid chromatograph with mass spectrometry. These porewater samples were collected from a microdialysis which generated samples along the surface of a growing A. sative root (all_hc). This data was used to answer questions connecting rhizosphere metabolite (putatively identified metabolites, hc_putative_norm) changes over time (root growth) with changes in the surrounding rhizosphere biogeochemistry (DOC, redox, pH). Rhizosphere biogeochemistry values are provided in the hc_putative_norm file as averages over their respective range of time that they were collected at. The hc_putative_norm file also contains all normalized values over only the intensity values collected for putatively identified metabolites.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Simulating physical and biogeochemical dynamics of the Arctic Ocean [Slides]

Ocean physics (currents, mixing, air-sea transfer of heat and gases) and biogeochemistry (carbon and oxygen cycles, ecosystem dynamics, nutrient distributions) are tightly coupled. We used the DOE-developed Energy Exascale Earth System Model (E3SM) to simulate physical and biogeochemical dynamics of the global ocean, with an emphasis on the western North Pacific. In particular, we used the grid refinement capability of E3SM’s ocean component (LANL-developed MPASOcean) to quantify the effect of model resolution on the carbon and oxygen cycles using a pair of century-long simulations that differ only in the western North Pacific grid. Regional grid refinement has a noticeable impact on biogeochemistry. Distribution of important nutrients (such as nitrate) changes the amount of photosynthesis by phytoplankton, which in turn affects air-sea transfer of CO2 and the local carbon cycle.

58 GEOSCIENCES↗
DOE OSTI2024

Geochemical and Microbial Dynamics of Hydrogen in a Methane Storage Reservoir

Hydrogen has been identified as a flexible energy carrier with zero or negative emission across multiple energy systems, and existing natural gas infrastructure could be leveraged if hydrogen gas (H2) was blended with methane (CH4). For example, subsurface methane storage reservoirs could be slightly modified to also store hydrogen if a methane/hydrogen blend were injected. However, the compatibility of methane storage reservoirs to include H2 injection has not been fully demonstrated, and this could lead to geochemical and microbiological reactions that alter the reservoir and stored gas content. It is essential that we understand the impact of H2 gas on the biogeochemistry of subsurface storage reservoirs before deploying large-scale H2-CH4 storage, We collected produced fluid from two separate methane storage reservoirs in the Southwestern US. First, we completed a baseline analysis of the biogeochemistry through qPCR, 16S rRNA sequencing, metagenomic sequencing, and geochemical analysis. Each reservoir was found to have unique geochemical conditions and a unique microbial community structure, with Site 1 having a higher TDS and an abundance of Shewanella and Site 2 having a lower TDS and high abundance of Eubacterium and Acetobacterium. Next, we ran a series of high pressure, high temperature reactors under hydrogen storage conditions with the biological sample from one of the storage reservoirs and a 20% H2-80% CH4 gas blend for up to 7 days. Our results show a decrease of hydrogen by 5% in reactors as early as 1-3 days. Previous hydrogen storage work has linked subsurface microorganisms with methanogenesis hydrogen sulfide production, acid production, and microbial corrosion. Our results show minimal change in the fluid chemistry, with the exception of a decrease in dissolved sulfate concentrations. Taxonomic sequencing demonstrated the presence of microorganisms capable of iron redox, acid generation, and hydrogen sulfide production throughout the reactors, suggesting microbial hydrogen consumption may occur through various metabolic pathways. This work demonstrates that site-specific geochemistry and microbiology may impact the efficiency of hydrogen storage in methane storage reservoirs.

environmental microbiology↗