Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Biogeochemical modeling”

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 523 records · Page 29

Do we have globally representative data to understand soil processes?

Understanding and modeling soils and soil organic matter (SOM) are central to a variety of human needs, from food production to ecosystem management. Soil data have been collected for over a century, but the global spatial and process representativeness of soil data remains unclear. We assessed the representativeness of currently available soil data that could be used to understand a variety of SOM processes. We used 16 open-source soil databases and data from over 281,000 unique locations globally, categorizing the databases into three main data types necessary to understand SOM processes: soil carbon stocks and fluxes, mechanistic drivers of these stocks and fluxes, and soil carbon gain or loss potential. We found that stock and driver data have extensive global coverage. However, data on soil carbon gain or loss potential, particularly data describing change in soils over time such as time series data, are severely limited in their global coverage. We conclude that while significant strides have been made in measuring soil carbon stocks and fluxes, and their drivers, we are limited in global data related to changes in soils over time. Our recommendations for soil data generators are to ensure precise metadata reporting and prioritizing sampling in underrepresented areas like tropical, arctic, mountainous, wetland and arid regions. We also encourage designing revisit schemes that explicitly support change detection and reporting multi-modal datasets that can aid in model development. Targeted measurement of low coverage soil data types and regions is necessary for a range of applications including current and future biogeochemical predictions, and their management and policy implications.

carbon fluxes↗

X-ray chemical imaging for assessing redox microsites within soils and sediments

Redox reactions underlie several biogeochemical processes and are typically spatiotemporally heterogeneous in soils and sediments. However, redox heterogeneity has yet to be incorporated into mainstream conceptualizations and modeling of soil biogeochemistry. Anoxic microsites, a defining feature of soil redox heterogeneity, are non-majority oxygen depleted zones in otherwise oxic environments. Neglecting to account for anoxic microsites can generate major uncertainties in quantitative assessments of greenhouse gas emissions, C sequestration, as well as nutrient and contaminant cycling at the ecosystem to global scales. However, only a few studies have observed/characterized anoxic microsites in undisturbed soils, primarily, because soil is opaque and microsites require µm-cm scale resolution over cm-m scales. Consequently, our current understanding of microsite characteristics does not support model parameterization. To resolve this knowledge gap, we demonstrate through this proof-of-concept study that X-ray fluorescence (XRF) 2D mapping can reliably detect, quantify, and provide basic redox characterization of anoxic microsites using solid phase “forensic” evidence. First, we tested and developed a systematic data processing approach to eliminate false positive redox microsites, i.e., artefacts, detected from synchrotron-based multiple-energy XRF 2D mapping of Fe (as a proxy of redox-sensitive elements) in Fe-“rich” sediment cores with artificially injected microsites. Then, spatial distribution of Fe II and Fe III species from full, natural soil core slices (over cm-m lengths/widths) were mapped at 1–100 µm resolution. These investigations revealed direct evidence of anoxic microsites in predominantly oxic soils such as from an oak savanna and toeslope soil of a mountainous watershed, where anaerobicity would typically not be expected. We also revealed preferential spatial distribution of redox microsites inside aggregates from oak savanna soils. We anticipate that this approach will advance our understanding of soil biogeochemistry and help resolve “anomalous” occurrences of reduced products in nominally oxic soils.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Development of observation-based global multilayer soil moisture products for 1970 to 2016

Abstract. Soil moisture (SM) datasets are critical to understanding the global water, energy, and biogeochemical cycles and benefit extensive societal applications. However, individual sources of SM data (e.g., in situ and satellite observations, reanalysis, offline land surface model simulations, Earth system model – ESM – simulations) have source-specific limitations and biases related to the spatiotemporal continuity, resolutions, and modeling and retrieval assumptions. Here, we developed seven global, gap-free, long-term (1970–2016), multilayer (0–10, 10–30, 30–50, and 50–100 cm) SM products at monthly 0.5∘ resolution (available at https://doi.org/10.6084/m9.figshare.13661312.v1; Wang and Mao, 2021) by synthesizing a wide range of SM datasets using three statistical methods (unweighted averaging, optimal linear combination, and emergent constraint). The merged products outperformed their source datasets when evaluated with in situ observations (mean bias from −0.044 to 0.033 m3 m−3, root mean square errors from 0.076 to 0.104 m3 m−3, Pearson correlations from 0.35 to 0.67) and multiple gridded datasets that did not enter merging because of insufficient spatial, temporal, or soil layer coverage. Three of the new SM products, which were produced by applying any of the three merging methods to the source datasets excluding the ESMs, had lower bias and root mean square errors and higher correlations than the ESM-dependent merged products. The ESM-independent products also showed a better ability to capture historical large-scale drought events than the ESM-dependent products. The merged products generally showed reasonable temporal homogeneity and physically plausible global sensitivities to observed meteorological factors, except that the ESM-dependent products underestimated the low-frequency temporal variability in SM and overestimated the high-frequency variability for the 50–100 cm depth. Based on these evaluation results, the three ESM-independent products were finally recommended for future applications because of their better performances than the ESM-dependent ones. Despite uncertainties in the raw SM datasets and fusion methods, these hybrid products create added value over existing SM datasets because of the performance improvement and harmonized spatial, temporal, and vertical coverages, and they provide a new foundation for scientific investigation and resource management.

54 ENVIRONMENTAL SCIENCES↗

Hydrology Outweighs Temperature in Driving Production and Export of Dissolved Carbon in a Snowy Mountain Catchment

Abstract Terrestrial production and export of dissolved organic and inorganic carbon (DOC and DIC) to streams depends on water flow and biogeochemical processes in and beneath soils. Yet, understanding of these processes in a rapidly changing climate is limited. Using the watershed‐scale reactive‐transport model BioRT‐HBV and stream data from a snow‐dominated catchment in the Rockies, we show deeper groundwater flow averaged about 20% of annual discharge, rising to ∼35% in drier years. DOC and DIC production and export peaked during snowmelt and wet years, driven more by hydrology than temperature. DOC was primarily produced in shallow soils (1.94 ± 1.45 gC/m 2 /year), stored via sorption, and flushed out during snowmelt. Some DOC was recharged to and further consumed in the deeper subsurface via respiration (−0.27 ± 0.02 gC/m 2 /year), therefore reducing concentrations in deeper groundwater and stream DOC concentrations at low discharge. Consequently, DOC was primarily exported from the shallow zone (1.62 ± 0.96 gC/m 2 /year, compared to 0.12 ± 0.02 gC/m 2 /year from the deeper zone). DIC was produced in both zones but at higher rates in shallow soils (1.34 ± 1.00 gC/m 2 /year) than in the deep subsurface (0.36 ± 0.02 gC/m 2 /year). Deep respiration elevated DIC concentrations in the deep zone and stream DIC concentrations at low discharge. In other words, deep respiration is responsible for the commonly‐observed increasing DOC concentrations (flushing) and decreasing DIC concentrations (dilution) with increasing discharge. DIC export from the shallow zone was ~66% of annual export but can drop to ∼53% in drier years. Numerical experiments suggest lower carbon production and export in a warmer, drier future, and a higher proportion from deeper flow and respiration processes. These results underscore the often‐overlooked but growing importance of deeper processes in a warming climate.

54 ENVIRONMENTAL SCIENCES↗

The carbon isotope biogeochemistry of (epsilon)CO2 production in a methanogenic marine sediment

To investigate the relationship between sigma(CO2) delta(C-13) values and rates of the dominant remineralization processes at the organic-rich field site of Cape Lookout Bight, NC, the isotopic composition of porewater sigma(CO2) was measured on a seasonal basis. The sigma(CO2) delta(C-13) values varies seasonally in response to changes in rates of sulfate reduction and methanogenesis, the dominant remineralization processes at this site. A tube incubation experiment was also performed to determine the isotopic signature of the sigma(CO2) produced by sulfate reduction and methanogenesis. The delta(C-13) of the sigma(CO2) produced in the sulfate reduction zone determined from the tube incubation was -14.3 plus or minus 1.9, a value enriched in C-13 relative to the labile organic fraction. The C-13-enrichment may be caused by low rates of methanogenesis occurring in the sulfate reduction zone. The delta(C-13) of the sigma(CO2) produced in the methanogenic zone was estimated to be +44 per mil, whereas the co-produced methane was -65 per mil. The fractionation factor for CO2 reduction was calculated to be 1.055, a value in agreement with previous estimates at this site. The measured concentration and delta(C-13) of the sigma(CO2) at Cape Lookout was closely reproduced by a diagenetic model using the measured rates of sulfate reduction and sigma(CO2) production, and the isotopic signature of the sigma(CO2) production in the two biogeochemical zones.

Boehme, Susan E.↗

The Emergence and Promise of Functional Chemogeography of Organic Matter

Organisms in ecosystems continuously release a myriad of organic matter molecules that undergo microbial and abiotic transformation, processes that critically influence carbon storage and climate feedbacks. Yet, a systematic understanding of what determines the transformation and persistence of organic matter across spatiotemporal scales remains elusive. We propose an emerging framework, termed “functional chemogeography,” to understand transformation and persistence of organic matter based on the chemical traits of molecules. This framework extends beyond a sole focus on intrinsic traits, which remain relatively constant across spatiotemporal scales, to emphasize extrinsic traits such as biochemical transformations and environmental responses, which vary spatiotemporally and are shaped by both intrinsic traits and the environment. When upscaled to the assemblage level using functional diversity indices, these extrinsic traits reveal a significant, and in some cases superior, capacity than intrinsic traits to explain biogeochemical processes, as demonstrated through a case study of dissolved organic matter in China's lakes. In conclusion, by integrating trait-based perspectives into predictive models, this framework helps bridge chemical complexity with ecosystem biogeochemistry, thereby advancing our ability to predict the fate of global organic carbon under environmental change.

chemical traits↗

Microbial ecology of extreme environments: Antarctic dry valley yeasts and growth in substrate limited habitats

The multiple stresses temperature, moisture, and for chemoheterotrophs, sources of carbon and energy of the Dry Valley Antarctica soils allow at best depauperate communities, low in species diversity and population density. The nature of community structure, the operation of biogeochemical cycles, the evolution and mechanisms of adaptation to this habitat are of interest in informing speculations upon life on other planets as well as in modeling the limits of gene life. Yeasts of the Cryptococcus vishniacil complex (Basidiobiastomycetes) are investigated, as the only known indigenes of the most hostile, lichen free, parts of the Dry Valleys. Methods were developed for isolating these yeasts (methods which do not exclude the recovery of other microbiota). The definition of the complex was refined and the importance of nitrogen sources was established as well as substrate competition in fitness to the Dry Valley habitats.

Vishniac, H. S.↗

Hybridizing Machine Learning and Physically-based Earth System Models to Improve Prediction of Multivariate Extreme Events (AI Exploration of Wildland Fire Prediction)

Focal Areas: This project responds to two focal areas identified in the DOE Call for AI4ESP White Papers: 1) Predictive modeling through the use of artificial intelligence (AI) techniques, and 2) insights gleaned from complex data using explainable AI and big data analytics. Science Challenge: Large wildland fires (hereafter wildfires) appearing as high-impact compound climate extreme events are closely related to hydroclimate and water cycle extremes that modulate surface fuel supply and combustibility. These compound events have multivariate climatic features (e.g., temperature, precipitation, relative humidity, wind, lightning) and societal drivers (e.g., forest management, land use change, human caused ignitions). Meanwhile, they induce strong feedbacks to the coupled atmosphere, biosphere, and hydrosphere by perturbing regional and global radiation budget as well as ecological, biogeochemical, and water cycles across multiple spatiotemporal scales. The nonlinear interactions between these natural and anthropogenic components of the Earth system are too complex to be completely and adequately represented in today’s Earth system models (ESMs). The inherent stochastic nature of fire activity at all scales further increases the difficulty of its prediction using ESMs that are usually developed from deterministic equations and parameterizations. Besides, concurrence of long-term (decadal to interdecadal) global climate change and fire regime shifts overlapping with short-term (intraseasonal to interannual) variations of regional fire weather and burning activity confound predictability of these compound extreme events. We propose to address the above scientific challenges by using machine learning (ML)-based data-driven modeling techniques to integrate observations and physically-based ESMs’ simulations in a computationally efficient hybrid prediction system. This prediction system is supposed to characterize the wildfire’s sensitivity to climate and exogenous drivers at high resolution (~ 0.25°) on subseasonal to seasonal (S2S) timescales providing improved predictability and explainability. We will use the system to help identify: (1) What are the computational elements of a hybrid system needed to predict compound climate extreme events such as global wildfires? (2) What are the key drivers (either natural or anthropogenic) that modulate short-term variations of multivariate fire weather and burning activity over different regions? How can one take advantage of those driver-response relationships to improve the predictability of large wildfires on S2S time scales? (3) What are the underlying physical mechanisms and sources of improved predictability? Which ML techniques are optimal in revealing and adapting these mechanisms?

54 ENVIRONMENTAL SCIENCES↗

Contribution of the World's Main Dust Source Regions to the Global Cycle of Desert Dust

Even though desert dust is the most abundant aerosol by mass in Earth's atmosphere, the relative contributions of the world's major source regions to the global dust cycle remain poorly constrained. This problem hinders accounting for the potentially large impact of regional differences in dust properties on clouds, the Earth's energy balance, and terrestrial and marine biogeochemical cycles. Here, we constrain the contribution of each of the world's main dust source regions to the global dust cycle. We use an analytical framework that integrates an ensemble of global aerosol model simulations with observationally informed constraints on the dust size distribution, extinction efficiency, and regional dust aerosol optical depth (DAOD). We obtain a dataset that constrains the relative contribution of nine major source regions to size-resolved dust emission, atmospheric loading, DAOD, concentration, and deposition flux. We find that the 22–29 Tg (1 standard error range) global loading of dust with a geometric diameter up to 20 µm is partitioned as follows: North African source regions contribute ∼ 50 % (11–15 Tg), Asian source regions contribute ∼ 40 % (8–13 Tg), and North American and Southern Hemisphere regions contribute ∼ 10 % (1.8–3.2 Tg). These results suggest that current models on average overestimate the contribution of North African sources to atmospheric dust loading at ∼ 65 %, while underestimating the contribution of Asian dust at ∼ 30 %. Our results further show that each source region's dust loading peaks in local spring and summer, which is partially driven by increased dust lifetime in those seasons. We also quantify the dust deposition flux to the Amazon rainforest to be ∼ 10 Tg yr−1, which is a factor of 2–3 less than inferred from satellite data by previous work that likely overestimated dust deposition by underestimating the dust mass extinction efficiency. The data obtained in this paper can be used to obtain improved constraints on dust impacts on clouds, climate, biogeochemical cycles, and other parts of the Earth system.

Dust↗

Laser Induced Fluorescence (LIF) as a Remote Sensing Tool: A Review

Vegetational changes are primary indicators of the present and future ecological status of the globe. These are changes which not only impact upon the primary productivity, but the total of the biogeochemical processes occurring on the planet. The impacts of global climatic and other environmental changes on vegetation must be monitored by some means in order to develop models which will allow us to predict long term effects. Large scale monitoring is now possible only with remote sensing systems, primarily passive reflectance, obtained by the use of satellite and aircraft platforms. However, passive reflectance techniques at this time are limited in their ability to detect subtle changes in the concentration and oxidation states of the many compounds involved in the light reactions of photosynthesis. Knowledge of these changes we consider to be fundamental in the remote assessment of both the rate and efficiency of photosynthesis and also the early detection of stress damage. The above factors pointed to the desirability of a sensing technique with the sensitivity and specificity necessary for detecting and quantifying those biological entities involved in photosynthesis. Another optical technique for vegetation monitoring is fluorescence. Previously, the lack of adequate excitation light sources and detector technologies have limited the use of fluorescence on intact plant leaves in the field. It is only recently with the advent of lasers with short pulse duration and advanced detector technologies that fluorescence measurements in the remote mode have become possible in the presence of ambient light.

Chappelle, E. W.↗

The Ocean Carbon States Database: A Proof-of-Concept Application of Cluster Analysis in the Ocean Carbon Cycle

In this paper, we present a database of the basic regimes of the carbon cycle in the ocean, the 'ocean carbon states', as obtained using a data mining/pattern recognition technique in observation-based as well as model data. The goal of this study is to establish a new data analysis methodology, test it and assess its utility in providing more insights into the regional and temporal variability of the marine carbon cycle. This is important as advanced data mining techniques are becoming widely used in climate and Earth sciences and in particular in studies of the global carbon cycle, where the interaction of physical and biogeochemical drivers confounds our ability to accurately describe, understand, and predict CO2 concentrations and their changes in the major planetary carbon reservoirs. In this proof-of-concept study, we focus on using well-understood data that are based on observations, as well as model results from the NASA Goddard Institute for Space Studies (GISS) climate model. Our analysis shows that ocean carbon states are associated with the subtropical-subpolar gyre during the colder months of the year and the tropics during the warmer season in the North Atlantic basin. Conversely, in the Southern Ocean, the ocean carbon states can be associated with the subtropical and Antarctic convergence zones in the warmer season and the coastal Antarctic divergence zone in the colder season. With respect to model evaluation, we find that the GISS model reproduces the cold and warm season regimes more skillfully in the North Atlantic than in the Southern Ocean and matches the observed seasonality better than the spatial distribution of the regimes. Finally, the ocean carbon states provide useful information in the model error attribution. Model air-sea CO2 flux biases in the North Atlantic stem from wind speed and salinity biases in the subpolar region and nutrient and wind speed biases in the subtropics and tropics. Nutrient biases are shown to be most important in the Southern Ocean flux bias.

carbon cycle↗

Substrate-explicit Modeling Tutorial

This tutorial Narrative demonstrates how to apply thermodynamic theory, also known as lambda theory, to convert molecular formulas of compounds (derived from FTICR-MS peaks using Formulatiry and R codes) into stoichiometric and kinetic forms of biogeochemical reactions and how to use the resulting kinetic equations to simulate dynamic conversion of compounds in batch and continuous stirred tank reactors. It will address under what conditions respiration rates are driven by thermodynamics, and how respiration rates respond to the variations in parameters and input variables and how to interpret the results.

54 ENVIRONMENTAL SCIENCES↗

Phage resistance mutations in a marine bacterium impact biogeochemically relevant cellular processes

Phage–bacteria interactions shape ecology and biogeochemistry across biomes. Resistance, arising from their evolutionary arms race, is well documented for receptor mutations, but other resistance mechanisms and their ecological implications remain unexplored. Here we isolated, sequenced and characterized 13 phage-resistant mutants of marine Cellulophaga baltica (Flavobacteriia). Mechanistically, mutations in surface proteins provided broad and complete extracellular resistance against multiple phages through decreased adsorption. Intracellular mutations affecting serine, glycine and threonine metabolism produced narrower resistance against a single phage, permitting viral DNA replication, and, in one mutant, were shown to be lipid mediated. Putative ecosystem impacts inferred from in vitro experiments include: (1) altered carbon utilization for all mutants, but especially by surface ones, (2) increased metabolite secretion for one modelled intracellular mutant (including experimentally verified acetate) and (3) increased ‘stickiness’ for all mutants, with surface mutants also sedimenting faster. Our findings highlight new resistance mechanisms and suggest that the phage–host arms race could result in ecosystem-level biogeochemical impacts in marine microorganisms.

Urvoy, Marion [Ohio State University]↗

Macroecological distributions of gene variants highlight the functional organization of soil microbial systems

Abstract The recent application of macroecological tools and concepts has made it possible to identify consistent patterns in the distribution of microbial biodiversity, which greatly improved our understanding of the microbial world at large scales. However, the distribution of microbial functions remains largely uncharted from the macroecological point of view. Here, we used macroecological models to examine how the genes encoding the functional capabilities of microorganisms are distributed within and across soil systems. Models built using functional gene array data from 818 soil microbial communities showed that the occupancy-frequency distributions of genes were bimodal in every studied site, and that their rank-abundance distributions were best described by a lognormal model. In addition, the relationships between gene occupancy and abundance were positive in all sites. This allowed us to identify genes with high abundance and ubiquitous distribution (core) and genes with low abundance and limited spatial distribution (satellites), and to show that they encode different sets of microbial traits. Common genes encode microbial traits related to the main biogeochemical cycles (C, N, P and S) while rare genes encode traits related to adaptation to environmental stresses, such as nutrient limitation, resistance to heavy metals and degradation of xenobiotics. Overall, this study characterized for the first time the distribution of microbial functional genes within soil systems, and highlight the interest of macroecological models for understanding the functional organization of microbial systems across spatial scales.

59 BASIC BIOLOGICAL SCIENCES↗

Matrix Diffusion Controls Mountain Hillslope Groundwater Ages and Inferred Storage Dynamics

Groundwater age distributions provide fundamental insights on coupled water and biogeochemical processes in mountain watersheds. Field-based studies have found mixtures of young and old-aged groundwater in mountain catchments underlain by bedrock; yet, the processes that dictate these groundwater age distributions are poorly understood. In this work, we use the coupled ParFlow-CLM integrated hydrologic and EcoSLIM particle tracking models to simulate groundwater age distributions on a lower montane hillslope in the East River Watershed, Colorado (USA). We develop a convolution-based approach to propagate fracture-matrix diffusion processes to the EcoSLIM advection-dominated age distributions. We compare observed 3 H and 4 He concentrations from two groundwater wells against model predictions that have varying advective transport times and matrix diffusion magnitudes. Based on a Monte Carlo analysis that considers uncertain matrix and fracture parameters, we find that matrix diffusion is needed to jointly predict 3 H and 4 He observations at both wells. The advection-dominated age distributions lack adequate mixing of young and old-aged water to capture the observed co-occurrence of 3 H and 4 He. The model scenario that best matches the 3 H, 4 He, and water level observations when considering both advective flowpath and matrix diffusion mixing processes has a dynamic bedrock groundwater reservoir that is susceptible to considerable storage losses during low-snow periods. This dynamic groundwater system amplifies the need to assimilate deeper bedrock groundwater into watershed hydro-biogeochemical predictions. This work further highlights the importance of considering matrix diffusion when interpreting environmental tracers in bedrock groundwater systems.

54 ENVIRONMENTAL SCIENCES↗

Biogeochemical cycling in the ocean. Part 1: Introduction to the effects of upwelling along the west coast of North America

Coastal upwelling is examined as it relates to the cycling of chemical species in coastal waters along the west coast of North America. The temporal and spatial features of upwelling phenomena in the Eastern boundary regions of the North Pacific Ocean are presented and discussed in terms of upwelling episodes. Climate conditions affecting upwelling include: thermal effects, wind-induced shear stress which moves surface layers, and the curl of the wind stress vector which is thought to affect the extent and nature of upwelling and the formation of offshore convergent downwelling fronts. These effects and the interaction of sunlight and upwelled nutrients which result in a biological bloom in surface waters is modeled analytically. The roles of biological and chemical species, including the effects of predation, are discussed in that context, and relevant remote sensing and in situ observations are presented. Climatological, oceanographic, biological, physical, chemical events, and processes that pertain to biogeochemical cycling are presented and described by a set of partial differential equations. Simple preliminary results are obtained and are compared with data. Thus a fairly general framework has been laid where the many facets of biogeochemical cycling in coastal upwelled waters can be examined in their relationship to one another, and to the whole, to whatever level of detail or approximation is warranted or desired.

Howe, John T.↗

Data and scripts associated with “When do Riverine Systems 'Feel the Burn'? Simulating How Burn Extent and Severity Modulate Hydrologic Controls on Biogeochemical Export” (v2)

This data package is associated with the publication “When do Riverine Systems 'Feel the Burn'? Simulating How Burn Extent and Severity Modulate Hydrologic Controls on Biogeochemical Export” published in Water Resources Research (Wampler et al. 2025; preprint: https://doi.org/10.22541/essoar.174438106.63564767/v1). This study used the Soil and Water Assessment Tool (SWAT), a processed based model to explore the impacts of area burned and burn severity on streamflow, nitrate, and dissolved organic carbon (DOC) in two test basins: a semi-arid, mixed land use basin and a humid, primarily forested basin. We developed 1800 wildfire scenarios that we ran in each basin: 20 different burn extents (5 to 100% by 5%), 3 different burn severities (low, moderate, and high), and 30 different post-fire precipitation scenarios. We also ran an additional 30 scenarios associated with no wildfire for the 30 post-fire precipitation scenarios. For each scenario we were interested in the change in runoff ratio (streamflow) and average concentration and annual loads (nitrate and DOC) across the wildfire scenarios. This data package contains the data and scripts required to build SWAT models for the two test basins, create and run the wildfire scenarios, and generate the data summaries and figures used in the associated manuscript. This data package was originally published in March 2025. It was updated in January 2026 (v2; new and modified files) to include the final files after the manuscript went through reviews. See the change history section below for more details. For details on how to navigate data packages generated by this project, see https://data.ess-dive.lbl.gov/portals/PNNLRiverCorridorSFA/About.

54 ENVIRONMENTAL SCIENCES↗

Biogeochemical Processes Across Aquatic Interfaces

The aquatic interfaces exposing terrestrial soils to oxic-anoxic regime shifts represent biogeochemical “hotspots” that are extremely sensitive to climate and environmental change. However, processes and interaction across theses aquatic interfaces are poorly understood and underrepresented in current Earth system models. In this project, we aim to develop predictive understanding of the feedbacks between microbial systems and geochemical environments that determine emergent ecosystem behaviors and resilience in response to disturbances. We use experimental, mechanistic modeling and meta-analysis tools to elucidate interactions among soil, water, geomorphology and microbiology that regulate the molecular transformations and fluxes of carbon, nutrients, and redox-sensitive compounds across aquatic interfaces.

58 GEOSCIENCES↗