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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Hydrological Controls on Riverbed Methane Emissions: A Numerical Investigation of Hydrodynamic and Ebullitive Mechanisms from Site to Basin Scales

Methane production and emission from riverbed sediments constitute a significant yet underexplored component of the river methane budget. Despite their importance, a mechanistic and quantitative understanding of these processes under dynamic flow conditions remains elusive, particularly at the basin scale. In this study, we investigate the multiscale (site-to-basin) mechanisms governing methane emissions from riverbeds, focusing on the interplay between biogeochemical processes and hydrological dynamics. We developed a reactive transport model that integrates methane production with hydrodynamic and ebullitive pathways and coupled it with a basin-scale groundwater flow model to simulate methane emissions across scales. By incorporating key factors such as vertical hydrological exchange flow (VHEF), particulate organic carbon availability, sediment properties, and temperature sensitivity, our model captures the spatiotemporal variability of methane ebullition at both the site and basin scales. Our results reveal that methane ebullition is strongly modulated by VHEF, with lower flow rates promoting methane accumulation and subsequent ebullition. Sensitivity analyses underscore nonlinear responses of methane fluxes to environmental drivers. Furthermore, these findings emphasize the critical roles of sediment composition, hydrological exchange, and environmental conditions in shaping methane dynamics. Overall, this work advances our understanding of riverine methane emissions and offers insights for guiding climate models and mitigation strategies.

Fluxes↗
DOE OSTI2023

Sea–Level Rise Impacts on Tidal Marshes and Estuarine Biogeochemical Processes

We used a numerical model to investigate the effects of sea-level rise (SLR) on the biogeochemical processes in the York River Estuary with extensive tidal marshes. The fully-coupled hydrodynamic-water quality-marsh model accounts for the spatial and temporal variations of physical-biogeochemical interactions between the tidal marshes and surrounding waters. This study focuses on an SLR scenario where the vertical accretion of tidal marshes keeps pace with the rising sea levels. Results show that SLR amplifies the tidal range and prolongs flooding duration, which results in enhanced porewater exchanges of materials between the tidal marshes and the surrounding waters. The increased availability of shallow-water habitats and enhanced light utilization in the shallow areas under SLR promote phytoplankton production (PP) in the shallow-water regions of the York River. Consequently, the organic carbon in the open water is fueled by the contributions from shallow waters and the enhanced export of organic carbon from the marshes. The change in the dissolved oxygen (DO) budget in the York River Estuary is attributed to changes in water column respiration, net metabolism of the benthic layer, reaeration, PP, and increased stratification under SLR. The net DO flux out of the York River increases at the York River mouth. Diel DO variation, especially in the marshes in the upper estuary, promotes phosphorus release from the sediment. The changes in dissolved nitrogen under SLR are relatively minimal.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Model Data Archive Associated with Manuscript "Fire-altered Carbon Pools Create Disturbance Memory in Stream Dissolved Organic Carbon"

This data package supports the publication “Fire-altered Carbon Pools Create Disturbance Memory in Stream Dissolved Organic Carbon” by Li et al. (2026). The package contains processed model inputs, configuration files, restart files, simulation outputs, scripts, and visualization products used to evaluate post-fire dissolved organic carbon (DOC) dynamics in the Naches River Watershed, Washington, USA, following the 2021 Schneider Springs Fire. The modeling workflow couples ELM-BGC, the biogeochemistry-enabled Energy Exascale Earth System Model Land Model; ATS, the Advanced Terrestrial Simulator for integrated surface-subsurface hydrology; and PFLOTRAN, a reactive transport model for multicomponent aqueous geochemistry. Together, these models simulate how wildfire-induced changes in vegetation, litter, coarse woody debris, and soil organic matter influence DOC production, transport, and reaction from burned hillslopes to stream networks. The archive includes preprocessed meteorological, geospatial, hydrologic, and biogeochemical forcing data; ELM-BGC-derived DOC source terms; ATS mesh files; PFLOTRAN reactive-transport inputs; model configuration files; spin-up and transient restart files; watershed-scale diagnostic outputs; stream concentration time series; and figures or visualization files used to inspect and reproduce key results. File types include Hierarchical Data Format 5 (HDF5) files for gridded forcing and model-coupling data, model input and configuration files for ELM-BGC, ATS, and PFLOTRAN, restart and simulation-output files generated by the modeling workflow, tabular or time-series diagnostic outputs, scripts for post-processing and figure generation, and image or visualization products associated with the manuscript. Use of the package depends on the intended task. Re-running the simulations requires the relevant modeling software, including ELM-BGC, ATS, and PFLOTRAN as ATS's geochemical engine. Inspecting outputs and reproducing figures requires Python with scientific plotting libraries such as Matplotlib, and three-dimensional model outputs may be viewed with ParaView. Geographic information system files or maps may be inspected with ArcGIS Pro or comparable GIS software. The data package is intended to enable traceability, reuse, and partial reproduction of the coupled land-to-watershed hydro-biogeochemical modeling workflow used to test how wildfire disturbance affects terrestrial carbon pools and downstream DOC dynamics.

ATS↗
DOE OSTI2023

Seasonal drivers of dissolved oxygen across a tidal creek–marsh interface revealed by machine learning

Abstract Dissolved oxygen (DO) is a key biogeochemical control in coastal systems, and its concentration and drivers vary markedly through time and space. This makes it difficult to accurately represent coastal DO and associated biogeochemical processes in models, limiting our ability to predict how these systems will respond to global change. We obtained high‐frequency (5‐min) in situ measurements of DO collected at three locations across the interface of a tidal creek and coastal marsh in the Pacific Northwest, USA. Random Forest machine learning models quantified the importance of three categories of environmental drivers (Aquatic, Climatic, and Terrestrial) of DO variability across the creek–marsh interface. We selected two 4‐month datasets representing Summer and Winter seasonal periods to test two hypotheses on the dominant drivers of DO at the coastal interface. We found that the Terrestrial driver—characterized by long periods of anaerobic conditions and episodic pulses in DO after floods—was most important during the Winter, whereas the Aquatic driver—characterized by variability over tidal, diel, and lunar cycles—was most important during the Summer. We explored how future climate change scenarios could alter the drivers of DO variability using a cumulative sums driver–response framework. Our results suggest that under climate change, Aquatic and Climatic drivers may increase in importance during the Summer, potentially linked to changing metabolic regimes and sea level, with Terrestrial driver importance potentially increasing during the Winter. Our approach highlights useful methods for understanding the spatiotemporal complexity of oxygen across coastal interfaces and quantifying the relative importance of distinct environmental drivers.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Agrivoltaic Grazing Systems for a Sustainable Future: A Multi‐Disciplinary Review & Gap Analysis

Solar photovoltaics (PV) is the fastest growing source of electricity in the world, however, its expansion has raised concerns about the displacement of agriculture and the degradation of rangeland ecosystems. This has prompted research and development of multi-land-use systems that co-prioritize energy production alongside agriculture, commonly referred to as agrivoltaics. Currently, agrivoltaic grazing (also called solar grazing) is one of the most common forms of these systems. However, peer-reviewed research on these systems is limited. This review synthesizes research on agrivoltaic grazing systems, drawing upon literature addressing agrivoltaics broadly and considering relevant literature in adjacent fields of study. Based on this review, we identify six key gaps and priority directions for future research. These include foundational fieldwork to characterize the layered ecological impacts of solar PV and grazing and understand PV-livestock integration from both an animal welfare and a PV infrastructure perspective. This will facilitate the development of biogeochemical and economic models that improve our predictive capacity and ability to compare system designs. We also identify gaps in understanding the human aspects of these systems and emphasize the importance of utilizing collaborative research methods and increasing research on the social dimensions of agrivoltaic grazing systems.

14 SOLAR ENERGY↗
DOE OSTI2020

Ecosystem–scale biogeochemical fluxes from three bioenergy crop candidates: How energy sorghum compares to maize and miscanthus

Perennial crops have been the focus of bioenergy research and development for their sustainability benefits associated with high soil carbon (C) and reduced nitrogen (N) requirements. However, perennial crops mature over several years and their sustainability benefits can be negated through land reversion. A photoperiod–sensitive energy sorghum (Sorghum bicolor) may provide an annual crop alternative more ecologically sustainable than maize (Zea mays) that can more easily integrate into crop rotations than perennials, such as miscanthus (Miscanthus × giganteus). This study presents an ecosystem–scale comparison of C, N, water and energy fluxes from energy sorghum, maize and miscanthus during a typical growing season in the Midwest United States. Gross primary productivity (GPP) was highest for maize during the peak growing season at 21.83 g C m –2 day –1 , followed by energy sorghum (17.04 g C m –2 day –1 ) and miscanthus (15.57 g C m –2 day –1 ). Maize also had the highest peak growing season evapotranspiration at 5.39 mm day –1 , with energy sorghum and miscanthus at 3.81 and 3.61 mm day –1 , respectively. Energy sorghum was the most efficient water user (WUE), while maize and miscanthus were comparatively similar (3.04, 1.75 and 1.89 g C mm –1 H 2 O, respectively). Maize albedo was lower than energy sorghum and miscanthus (0.19, 0.26 and 0.24, respectively), but energy sorghum had a Bowen ratio closer to maize than miscanthus (0.12, 0.13 and 0.21, respectively). Nitrous oxide (N 2 O) flux was higher from maize and energy sorghum (8.86 and 12.04 kg N ha –1 , respectively) compared with miscanthus (0.51 kg N ha –1 ), indicative of their different agronomic management. Furthermore, these results are an important first look at how energy sorghum compares to maize and miscanthus grown in the Midwest United States. This quantitative assessment is a critical component for calibrating biogeochemical and ecological models used to forecast bioenergy crop growth, productivity and sustainability.

09 BIOMASS FUELS↗
DOE OSTI2020

Biogeochemical Fluxes for Three Bioenergy Crops

Perennial crops have been the focus of bioenergy research and development for their sustainability benefits associated with high soil carbon (C) and reduced nitrogen (N) requirements. However, perennial crops mature over several years and their sustainability benefits can be negated through land reversion. A photoperiod‐sensitive energy sorghum (Sorghum bicolor) may provide an annual crop alternative more ecologically sustainable than maize (Zea mays) that can more easily integrate into crop rotations than perennials, such as miscanthus (Miscanthus × giganteus). This study presents an ecosystem‐scale comparison of C, N, water and energy fluxes from energy sorghum, maize and miscanthus during a typical growing season in the Midwest United States. Gross primary productivity (GPP) was highest for maize during the peak growing season at 21.83 g C m−2 day−1, followed by energy sorghum (17.04 g C m−2 day−1) and miscanthus (15.57 g C m−2 day−1). Maize also had the highest peak growing season evapotranspiration at 5.39 mm day−1, with energy sorghum and miscanthus at 3.81 and 3.61 mm day−1, respectively. Energy sorghum was the most efficient water user (WUE), while maize and miscanthus were comparatively similar (3.04, 1.75 and 1.89 g C mm−1 H2O, respectively). Maize albedo was lower than energy sorghum and miscanthus (0.19, 0.26 and 0.24, respectively), but energy sorghum had a Bowen ratio closer to maize than miscanthus (0.12, 0.13 and 0.21, respectively). Nitrous oxide (N2O) flux was higher from maize and energy sorghum (8.86 and 12.04 kg N ha−1, respectively) compared with miscanthus (0.51 kg N ha−1), indicative of their different agronomic management. These results are an important first look at how energy sorghum compares to maize and miscanthus grown in the Midwest United States. This quantitative assessment is a critical component for calibrating biogeochemical and ecological models used to forecast bioenergy crop growth, productivity and sustainability.

Field Data↗
DOE OSTI2021

Reduction of Vanadium(V) by Iron(II)-Bearing Minerals

Fe(II)-bearing minerals (magnetite, siderite, green rust, etc.) are common products of microbial Fe(III) reduction, and they provide a reservoir of reducing capacity in many subsurface environments that may contribute to the reduction of redox active elements such as vanadium; which can exist as V(V), V(IV), and V(III) under conditions typical of near-surface aquatic and terrestrial environments. To better understand the redox behavior of V under ferrugenic/sulfidogenic conditions, we examined the interactions of V(V) (1 mM) in aqueous suspensions containing 50 mM Fe(II) as magnetite, siderite, vivianite, green rust, or mackinawite, using X-ray absorption spectroscopy at the V K-edge to determine the valence state of V. Two additional systems of increased complexity were also examined, containing either 60 mM Fe(II) as biogenic green rust (BioGR) or 40 mM Fe(II) as a mixture of biogenic siderite, mackinawite, and magnetite (BioSMM). Within 48 h, total solution-phase V concentrations decreased to <20 µM in all but the vivianite and the biogenic BiSMM systems; however, >99.5% of V was removed from solution in the BioSMM and vivianite systems within 7 and 20 months, respectively. The most rapid reduction was observed in the mackinawite system, where V(V) was reduced to V(III) within 48 h. Complete reduction of V(V) to V(III) occurred within 4 months in the green rust system, 7 months in the siderite system, and 20 months in the BioGR system. Vanadium(V) was only partially reduced in the magnetite, vivianite, and BioSMM systems, where within 7 months the average V valence state stabilized at 3.7, 3.7, and 3.4, respectively. The reduction of V(V) in soils and sediments has been largely attributed to microbial activity, presumably involving direct enzymatic reduction of V(V); however the reduction of V(V) by Fe(II)-bearing minerals suggests that abiotic or coupled biotic–abiotic processes may also play a critical role in V redox chemistry, and thus need to be considered in modeling the global biogeochemical cycling of V.

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

Comparing Model Representations of Physiological Limits on Transpiration at a Semi-arid Ponderosa Pine Site

Mechanistic representations of biogeochemical processes in ecosystem models are rapidly advancing, requiring advancements in model evaluation approaches. Here we quantify multiple aspects of model functional performance to evaluate improved process representations in ecosystem models. We compare semi-empirical stomatal models with hydraulic constraints against more mechanistic representations of stomatal and hydraulic functioning at a semi-arid pine site using a suite of metrics and analytical tools. We find that models generally perform similarly under unstressed conditions, but performance diverges under atmospheric and soil drought. The more empirical models better capture synergistic information flows between soil water potential and vapor pressure deficit to transpiration, while the more mechanistic models are overly deterministic. Although models can be parameterized to yield similar functional performance, alternate parameterizations could not overcome structural model constraints that underestimate the unique information contained in soil water potential about transpiration. Additionally, both multilayer canopy and big-leaf models were unable to capture the magnitude of canopy temperature divergence from air temperature, and we demonstrate that errors in leaf temperature can propagate to considerable error in simulated transpiration. This study demonstrates the value of merging underutilized observational data streams with emerging analytical tools to characterize ecosystem function and discriminate among model process representations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Dataset and scripts for manuscript "Using Neural Network Ensembles to Separate Ocean Biogeochemical and Physical Drivers of Phytoplankton Biogeography in Earth System Models"

Please note: The title of this version contains an updated title for the manuscript compared to the previous version of this dataset. This is only due to title updates during the peer review process for the manuscript. The zip file contains the scripts, functions, and source files for the manuscript titled "Using Neural Network Ensembles to Separate Ocean Biogeochemical and Physical Drivers of Phytoplankton Biogeography in Earth System Models." The manuscript has been submitted for peer review. Please consult the README file for information on the specifications of the files. These files may occasionally be updated to add annotations to the scripts to make them more user friendly and to correct any errors.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Machine Learning in Environmental Chemistry: Application to Surface Complexation Modeling

Environmental chemistry – or biogeochemistry – is the scientific discipline typically invoked when examining and quantifying groundwater or surface water contamination, and nutrient cycling in the environment. Over the last three decades, there have been significant advances in mechanistic model development to describe and predict these complex biogeochemical processes. In particular, surface complexation models (SCMs) have been developed to describe the rock/soil surface reactions of metals and radionuclides, and their partitioning between various mobile species in the aqueous phase or immobile species sorbed on solid surfaces. Often represented by a simplified linear isotherm constant – Kd – in reactive transport models, these reactions play a critical role in many environmental science applications; particularly in contamination risk assessments and nuclear waste disposal performance assessments. In the past several decades, efforts by various institutions across the world have focused on developing SCMs based on datasets from laboratory measurements, including the identification of key parameters such as equilibrium constants.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Carbon–concentration and carbon–climate feedbacks in CMIP6 models and their comparison to CMIP5 models

Results from the fully and biogeochemically coupled simulations in which CO 2 increases at a rate of 1 % yr –1 (1pctCO2) from its preindustrial value are analyzed to quantify the magnitude of carbon–concentration and carbon–climate feedback parameters which measure the response of ocean and terrestrial carbon pools to changes in atmospheric CO 2 concentration and the resulting change in global climate, respectively. The results are based on 11 comprehensive Earth system models from the most recent (sixth) Coupled Model Intercomparison Project (CMIP6) and compared with eight models from the fifth CMIP (CMIP5). The strength of the carbon–concentration feedback is of comparable magnitudes over land (mean ± standard deviation = 0.97 ± 0.40 PgC ppm –1 ) and ocean (0.79 ± 0.07 PgC ppm –1 ), while the carbon–climate feedback over land (–45.1 ± 50.6 PgC °C –1 ) is about 3 times larger than over ocean (–17.2 ± 5.0 PgC °C –1 ). The strength of both feedbacks is an order of magnitude more uncertain over land than over ocean as has been seen in existing studies. These values and their spread from 11 CMIP6 models have not changed significantly compared to CMIP5 models. The absolute values of feedback parameters are lower for land with models that include a representation of nitrogen cycle. The transient climate response to cumulative emissions (TCRE) from the 11 CMIP6 models considered here is 1.77 ± 0.37 °C EgC –1 and is similar to that found in CMIP5 models (1.63 ± 0.48 °C EgC –1 ) but with somewhat reduced model spread. The expressions for feedback parameters based on the fully and biogeochemically coupled configurations of the 1pctCO2 simulation are simplified when the small temperature change in the biogeochemically coupled simulation is ignored. Decomposition of the terms of these simplified expressions for the feedback parameters is used to gain insight into the reasons for differing responses among ocean and land carbon cycle models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Plankton energy flows using a global size-structured and trait-based model

Plankton community models are critical tools for understanding the processes that shape marine plankton communities, how plankton communities impact biogeochemical cycles, and the feedbacks between community structure and function. Here, using the flexible Marine Biogeochemistry Library (MARBL), we present the Size-based Plankton ECological TRAits (MARBL-SPECTRA) model, which is designed to represent a diverse plankton community while remaining computationally tractable. MARBL-SPECTRA is composed of nine phytoplankton and six zooplankton size classes represented using allometric scaling relationships for physiological traits and interactions within multiple functional types. MARBL-SPECTRA is embedded within the global ocean component of the Community Earth System Model (CESM) and simulates large-scale, emergent patterns in phytoplankton growth limitation, plankton phenology, plankton generation time, and trophic transfer efficiency. The model qualitatively reproduces observed global patterns of surface nutrients, chlorophyll biomass, net primary production, and the biogeographies of a range of plankton size classes. In addition, the model simulates how predator:prey dynamics and trophic efficiency vary across gradients in total ecosystem productivity. Shorter food chains that export proportionally more carbon from the surface to the ocean interior occur in productive, eutrophic regions, whereas in oligotrophic regions, the food chains are relatively long and export less organic matter from the surface. The union of functional type modeling with size-resolved, trait-based modeling approaches allows MARBL-SPECTRA to capture both large-scale elemental cycles and the structure of planktonic food webs affecting trophic transfer efficiency.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Effective kinetics driven by dynamic concentration gradients under coupled transport and reaction

Biogeochemical reaction kinetics are generally established from batch reactors where concentrations are uniform. In natural systems, many biogeochemical processes are characterized by spatially and temporally variable concentration gradients that often occur at scales which are not resolved by field measurements or biogeochemical and reactive transport models. Yet, it is not clear how these sub-scale chemical gradients affect reaction kinetics compared to batch kinetics. Here we investigate this question by studying the paradigmatic case of localized pulses of solute reacting with a solid or a dissolved species in excess. Additionally, we consider non-linear biogeochemical reactions, representative of mineral dissolution, adsorption and redox reactions, which we quantify using simplified power-law kinetics. The combined effect of diffusion and reaction leads to effective kinetics that differ quantitatively and qualitatively from the batch kinetics. Depending on the nonlinearity (reaction order) of the local kinetics, these effects lead to either enhancement or decrease of the overall reaction rate, and result in a rich variety of reaction dynamics. We derive analytical results for the effective kinetics, which are validated by comparison to direct numerical simulations for a broad range of Damköhler numbers and reaction order. Our findings provide new insights into the interpretation of imperfectly mixed lab experiments, the effective kinetics of field systems characterized by intermittent reactant release and the integration of sub-scale concentration gradients in reactive transport models.

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

A new bioenergy model that simulates the impacts of plant‐microbial interactions, soil carbon protection, and mechanistic tillage on soil carbon cycling

Abstract Advancing our predictive understanding of bioenergy systems is critical to design decision tools that can inform which feedstock to plant, where to plant it, and how to manage its production to provide both energy and ecosystem carbon (C) benefits. Here, we lay the foundation for that advancement by integrating recent developments in the science of belowground processes in shaping the C cycle into a new bioenergy model, FUN‐BioCROP (Fixation and Uptake of Nitrogen‐Bioenergy Carbon, Rhizosphere, Organisms, and Protection). We show that FUN‐BioCROP can approximate the historical trajectory of soil C dynamics as natural ecosystems were successively converted into intensive agriculture and bioenergy systems. This ability relies in part on a novel tillage representation that mechanistically models tillage as a process that increases microbial access to C. Importantly, the impacts of tillage and feedstock choice also influence FUN‐BioCROP simulations of warming responses with no‐till perennial feedstocks, miscanthus, and switchgrass, having more C that is unprotected and susceptible to warming than tilled annual feedstocks like corn–corn–soybean. However, this susceptibility to warming is balanced by a greater potential for increases in belowground C allocation to enhance soil C stocks in perennial systems. Collectively, our model results highlight the importance of belowground processes in evaluating the ecosystem C benefits of bioenergy production.

09 BIOMASS FUELS↗
DOE OSTI2021

Toward more mechanistic representations of biogeochemical processes in stream and river networks: Implementation and demonstration of a multiscale model : Modeling Archive

This Modeling Archive supports the ORNL-SFA and IDEAS-Waterhsed publication submitted to Environmental Modeling and Software. A recently introduced multiscale model for representing the combined effect of hyporheic exchange flows and small-scale hyporheic-zone biogeochemical processes is extended from the reach scale to river network scales. The ADELS (Advection Dispersion Equation with Lagrangian Subgrid) model uses a one-dimensional advection-dispersion-reaction equation for the channel network and couples that equation at each channel location to a one-dimensional advection-reaction subgrid model representing an ensemble of streamlines that are diverted into the hyporheic zone before returning to the channel. The subgrid model is written in a computationally advantageous Lagrangian form with hyporheic age replacing the hyporheic travel distance. In the paper, we summarized implementation in the integrated surface/subsurface hydrology modeling system Advanced Terrestrial Simulator (ATS).

East Fork Poplar Creek↗
DOE OSTI2023

Simulations of ocean deoxygenation in the historical era: insights from forced and coupled models

Ocean deoxygenation due to anthropogenic warming represents a major threat to marine ecosystems and fisheries. Challenges remain in simulating the modern observed changes in the dissolved oxygen (O 2 ). Here, we present an analysis of upper ocean (0-700m) deoxygenation in recent decades from a suite of the Coupled Model Intercomparison Project phase 6 (CMIP6) ocean biogeochemical simulations. The physics and biogeochemical simulations include both ocean-only (the Ocean Model Intercomparison Project Phase 1 and 2, OMIP1 and OMIP2) and coupled Earth system (CMIP6 Historical) configurations. We examine simulated changes in the O 2 inventory and ocean heat content (OHC) over the past 5 decades across models. The models simulate spatially divergent evolution of O 2 trends over the past 5 decades. The trend (multi-model mean and spread) for upper ocean global O 2 inventory for each of the MIP simulations over the past 5 decades is 0.03 ± 0.39×10 14 [mol/decade] for OMIP1, -0.37 ± 0.15×10 14 [mol/decade] for OMIP2, and -1.06 ± 0.68×10 14 [mol/decade] for CMIP6 Historical, respectively. The trend in the upper ocean global O2 inventory for the latest observations based on the World Ocean Database 2018 is -0.98×10 14 [mol/decade], in line with the CMIP6 Historical multi-model mean, though this recent observations-based trend estimate is weaker than previously reported trends. A comparison across ocean-only simulations from OMIP1 and OMIP2 suggests that differences in atmospheric forcing such as surface wind explain the simulated divergence across configurations in O 2 inventory changes. Additionally, a comparison of coupled model simulations from the CMIP6 Historical configuration indicates that differences in background mean states due to differences in spin-up duration and equilibrium states result in substantial differences in the climate change response of O 2 . Finally, we discuss gaps and uncertainties in both ocean biogeochemical simulations and observations and explore possible future coordinated ocean biogeochemistry simulations to fill in gaps and unravel the mechanisms controlling the O 2 changes.

54 ENVIRONMENTAL SCIENCES↗