Search NASA⌕ Search

SEARCH · Search NASA

Results for “Dynamic 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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

The Integrated Virtual Blast Furnace: Enabling Physics-Based Operational Guidance

As part of a DOE-supported research effort, Purdue University Northwest researchers are collaborating with Oak Ridge National Laboratory and United States Steel Corporation to develop a tool to provide blast furnace operators and engineers with process performance insight comparable to high-fidelity computational fluid dynamics modeling, accelerated to provide “what-if” scenarios at near-real-time speed. This is accomplished by pre-simulating a baseline case and a range of potential operating scenarios to establish how the furnace responds to changing inputs, then training a neural network-based Reduced Order Model to accelerate the speed at which predictions of key parameters can be generated.

Okosun, Tyamo↗
DOE OSTI2024

In-vivo neuronal dysfunction by Aβ and tau overlaps with brain-wide inflammatory mechanisms in Alzheimer’s disease

The molecular mechanisms underlying neuronal dysfunction in Alzheimer’s disease (AD) remain uncharacterized. Here, we identify genes, molecular pathways and cellular components associated with whole-brain dysregulation caused by amyloid-beta (Aβ) and tau deposits in the living human brain. We obtained in-vivo resting-state functional MRI (rs-fMRI), Aβ- and tau-PET for 47 cognitively unimpaired and 16 AD participants from the Translational Biomarkers in Aging and Dementia cohort. Adverse neuronal activity impacts by Aβ and tau were quantified with personalized dynamical models by fitting pathology-mediated computational signals to the participant’s real rs-fMRIs. Then, we detected robust brain-wide associations between the spatial profiles of Aβ-tau impacts and gene expression in the neurotypical transcriptome (Allen Human Brain Atlas). Within the obtained distinctive signature of in-vivo neuronal dysfunction, several genes have prominent roles in microglial activation and in interactions with Aβ and tau. Moreover, cellular vulnerability estimations revealed strong association of microglial expression patterns with Aβ and tau’s synergistic impact on neuronal activity (q < 0.001). These results further support the central role of the immune system and neuroinflammatory pathways in AD pathogenesis. Neuronal dysregulation by AD pathologies also associated with neurotypical synaptic and developmental processes. In addition, we identified drug candidates from the vast LINCS library to halt or reduce the observed Aβ-tau effects on neuronal activity. Top-ranked pharmacological interventions target inflammatory, cancer and cardiovascular pathways, including specific medications undergoing clinical evaluation in AD. Our findings, based on the examination of molecular-pathological-functional interactions in humans, may accelerate the process of bringing effective therapies into clinical practice.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Experimental and Computational Characterization of a Modified Sioutas Cascade Impactor for Respirable Radioactive Aerosols

Oak Ridge National Laboratory is collecting and characterizing aerosols released when spent nuclear fuel (SNF) rods are fractured in bending. An aerosol collection system was designed and tested to collect respirable sized (<10 μm aerodynamic diameter [AED]) particulates inside a hot cell facility. The setup is a modified version of the commercially available Sioutas cascade impactor, to which additional stages were added to expand the aerosol collection range from 2.5 to ~15 μm AED. To accommodate the additional stages and specific test conditions, the operating flow rate for aerosol collection was reduced, and testing was conducted by using pressure drop measurements, surrogate dust collection, and particle size characterization. The fluid flow distribution within the cascade and its stages was simulated in STAR-CCM+, and the stage-wise pressure drops obtained using the computational fluid dynamics model were then compared to experimental data. Lagrangian particle simulations were also performed, and stage-wise collection statistics were obtained from the simulation for comparison with the experimental data obtained using SNF-surrogate dust particles. The results provide valuable insights into the stage-wise particle collection characteristics of the modified cascade impactor and can also be used to improve the prediction accuracy of the manufacturer-determined analytical correlations.

aerosol modeling↗
DOE OSTI2025

Powering the Woods Hole X-Spar Buoy with Ocean Wave Energy—A Control Co-Design Feasibility Study

Despite its success in measuring air–sea exchange, the Woods Hole Oceanographic Institution’s (WHOI) X-Spar Buoy faces operational limitations due to energy constraints, motivating the integration of an energy harvesting apparatus to improve its deployment duration and capabilities. This work explores the feasibility of an augmented, self-powered system in two parts. Part 1 presents the collaborative design between X-Spar developers and wave energy researchers translating user needs into specific functional requirements. Based on requirements like desired power levels, deployability, survivability, and minimal interference with environmental data collection, unsuitable concepts are pre-eliminated from further feasibility study consideration. In part 2, we focus on one of the promising concepts: an internal rigid body wave energy converter. We apply control co-design methods to consider commercial of the shelf hardware components in the dynamic models and investigate the concept’s power conversion capabilities using linear 2-port wave-to-wire models with concurrently optimized control algorithms that are distinct for every considered hardware configuration. During this feasibility study we utilize two different control algorithms, the numerically optimal (but acausal) benchmark and the optimized damping feedback. We assess the sensitivity of average power to variations in drive-train friction, a parameter with high uncertainty, and analyze stroke limitations to ensure operational constraints are met. Our results indicate that a well-designed power take-off (PTO) system could significantly extend the WEC-Spar’s mission by providing additional electrical power without compromising data quality.

autonomous systems↗
DOE OSTI2025

DELVE-ing into the Milky Way’s Globular Clusters: Assessing Extratidal Features in NGC 5897, NGC 7492, and Testing Detectability with Deeper Photometry

Extratidal features around globular clusters (GCs) are tracers of their disruption, stellar stream formation, and their host’s gravitational potential. However, these features remain challenging to detect due to their low surface brightness. We conduct a systematic search for such features around 19 GCs in the DECam Local Volume Exploration (DELVE) survey Data Release 2, discovering a new extra-tidal envelope around NGC 5897 and find tentative evidence for an extended envelope surrounding NGC 7492. Through a combination of dynamical modeling and analyzing synthetic stellar populations, we demonstrate these envelopes may have formed through tidal disruption. We use these models to explore the detectability of these features in the upcoming Legacy Survey of Space and Time (LSST), finding that while LSST’s deeper photometry will enhance detection significance, additional methods for foreground removal like proper motions or metallicities may be important for robust stream detection. Our results both add to the sample of globular clusters with extratidal features and provide insights on interpreting similar features in current and upcoming data.

Chiti, A. [Univ. of Chicago, IL (United States); S↗
DOE OSTI2026

Applications of Multibody for Everybody (M4E) in Marine Energy

Multibody for Everybody (M4E) is an open-source symbolic dynamics modeling framework designed to automate the derivation of equations of motion and simulation of constrained multibody systems using the joint coordinate formulation.

16 TIDAL AND WAVE POWER↗
DOE OSTI2026

Summary of GFL Inverter Capability and Performance: Malawi [Slides]

This presentation provides a comprehensive comparison of inverter-based resource (IBR) technical requirements across multiple international grid codes - including IEEE 2800, ESO, EirGrid, and Malawi's evolving regulatory framework - to support capacity building for electricity-sector institutions in Malawi. It highlights key performance expectations for grid-following (GFL) inverters across critical domains such as primary and fast frequency response, ramping behavior, frequency and voltage ride-through capabilities, reactive power and voltage control, ROCOF withstand, phase jump tolerance, system strength considerations, protection schemes, and dynamic modeling requirements. Emphasizing a "Do no harm, do something, do more" philosophy, the document synthesizes best practices and operational benchmarks to guide the integration and stability of inverter-based resources within Malawi's power system. The presentation forms part of the SABESS Center of Excellence initiative supported by the Global Energy Alliance for People and the Planet.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Chemical clocks: using otolith geochemistry to enhance estimation of age and growth of white hake (Urophycis tenuis)

The white hake (Urophycis tenuis) is a groundfish distributed throughout the Gulf of Maine. Catch advice is based on stock assessments done with age-based population dynamics models; however, otolith aging is challenging because of unclear growth increments. To address this concern, we compared the consistency of aging with counts of visual annuli to that of aging with cycles of elemental concentrations measured by using laser ablation inductively coupled plasma mass spectrometry. We tested the hypothesis that oscillations in both environmental conditions and internal physiology through time influence uptake of elements during otolith mineralization. Concentrations of manganese, in comparison with those of the other investigated trace elements (magnesium, strontium, and barium), had the most promising correlation with visual growth increments (~100% age agreement, ±1 year), offering an additional tool to enhance increment identification. In our examination of 550 otoliths collected during 2007–2021, we found that white hake lived a maximum of 10.3 years and exhibited sexual dimorphism in maximum length and age. By using generated von Bertalanffy growth functions, L(t)=110(1–e–0.113(t+0.45)) for males and L(t)=140(1–e–0.113(t–0.30)) for females (where L(t) is length at time t), size and age at maturity were calculated for males (37.4 cm in total length (TL), 3.3 years) and females (47.4 cm TL, 4.2 years). These results demonstrate that otolith geochemistry can be used to improve the accuracy and precision of the estimation of fish age and maturity, even for challenging species.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Circular Economy for Photovoltaics in Service of Energy Transition

The challenge of energy transition is immediate and immense; current projections target 75 TW of photovoltaics (PV) capacity by 2050. While any transition to renewable energy technology is preferable to the current fossil-based system, it is ideal to improve the sustainability of PV to minimize negative environmental and social impacts. Circular economy (CE) has been proposed as a method to improve the sustainability of PV, especially for emerging materials like perovskites. CE is a set of actions, principles, and systems which aim to design out waste and keep products and materials in use, to reduce environmental impacts and enable sustainable development. At the most basic level, CE is "reduce, reuse, recycle", the R-actions, in ranked order. CE of a PV technology can be metricized in a variety of ways, such as the Material Circularity Indicator (Smith and Jones, Ellen MacArthur Foundation, 2019) or recycling rates. Unfortunately, standard CE metrics have several shortcomings for measuring renewable energy technologies in the context of deployment for energy transition (Figge 2018, Saidani 2019): 1) Only measure mass flows; 2) De-prioritization of the use phase in favor of mass circularity when scoring; and 3) Tight focus on a single product scale The use phase and energy flows of PV are key to energy transition, and therefore need to be quantified. Additionally, correlating product-scale to system-scale is necessary for quantifying the environmental impacts of energy transition. Life Cycle Assessment (LCA) can address some of these concerns, but also focuses on a single product scale and has trouble capturing the dynamics of system-scale energy transition, such as the interaction of module lifetime with manufacturing demands for energy transition deployment schedules. Therefore, we developed an open-source Python-based system dynamics model to quantify the mass, energy and carbon impacts of CE R-actions for PV technologies in the energy transition; PV in the CE (PV ICE) (Ovaitt & Mirletz 2021). The tool captures supply chains from material extraction through end of life, incorporating 5 circular end of life pathways. PV ICE takes in any evolving bill of materials, module properties and deployment schedule to support researchers and decision makers with data-backed insights. In this work, we quantify and compare proposed CE sustainable PV module designs and lifecycle management strategies, spanning currently commercialized technologies, government and industry technology targets, and several low Technology Readiness Level (TRL) emerging PV technologies, including perovskites. Our analyses capture the projected evolutions of lifetime, efficiency and material circularity of these PV technologies, as well as their material supply chains. Our analyses emphasize the importance of examining a suite of metrics to identify priorities and tradeoffs, and inform design or lifecycle management decisions holistically. Previous analyses have demonstrated the central importance of PV module lifetime to support energy transition while minimizing impacts. High levels of material circularity (>90%) enable minimizing lifecycle wastes, can reduce virgin material demands if paired with improving efficiency, but demonstrate tradeoffs in energy return on investment. In the fervor of new material and technology development, it is important to remember that CE is not the end goal; decarbonization and energy transition are the end goal. CE should be used in service to improve the sustainability of PV, and R-actions evaluated for their usefulness and efficacy to this end.

carbon↗
DOE OSTI2025

Novel Results Visualization for Dynamic PSA and New Modeling Features in EMRALD

The Event Modeling Risk Assessment Linked Diagram (EMRALD) tool, developed at the Idaho National Laboratory (INL), was designed to simplify the creation of dynamic models and support various research projects. One of the primary goals of EMRALD was to provide visual methods for modeling. EMRALD consists of two main components: a web-based user interface for model development and a solve engine for running model simulations. Over time, it has evolved to meet the diverse needs of its users. Initially, EMRALD's results were simple text outputs with final key state percentages and uncertainty bounds. However, because EMRALD utilizes a three-phase discrete event simulation and tracks the paths of each simulation run leading to a key state, there is significant potential to analyze large sets of path results data, including state paths, events, and timing. Visualizing this data meaningfully posed a challenge. To address this, a novel time-based Sankey diagram was developed. EMRALD exports results data in a format that can be opened in this Sankey viewer, allowing users to visualize paths, occurrences, events, and probability data for the entire simulation run in a single diagram. Moreover, when EMRALD was first created, there were limited tools capable of meeting its graphical requirements, many of which are no longer supported. In 2024, a new web-based interface was developed using modern graphing tools, enabling additional modeling features. This paper discusses the new dynamic PSA results visualization capability and the enhanced modeling tools available in EMRALD.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2025

How initial conditions-, structural-, and parameter-based model uncertainty interact and influence predictions in permafrost ecosystems: Modeling Archive

This dataset contains model output and input data, as well as source code examples for the Terrestrial Ecosystem Model with the Dynamic Vegetation Model and Dynamic Organic Soil (DVM-DOS-TEM) for the field sites Imnavait creek and the Bonanza creek Long Term Ecological Research Network (LTER). The data covers simulations from the last glacial maximum (LGM) until 2100 for a selection of paleo scenarios, setting the mean temperature of the LGM up to 10°C lower than pre-industrial conditions. The model structure was modulated to represent various model versions, and this dataset contains the relevant changes in the source code. The raw output data, the processed statistical data, the setup and processing scripts as well as parameter value distribution files from a parameter sensitivity analysis are included as well. Model outputs include active layer depth, organic soil carbon, soil layer depths, gross primary productivity (GPP) with and without nitrogen limitation, net primary productivity (NPP), soil liquid water content, heterotrophic, maintenance, and growth respiration, soil temperature, and vegetation carbon (*.nc files). The Next-Generation Ecosystem Experiments in the Arctic (NGEE Arctic) project is a research effort to reduce uncertainty in the Department of Energy’s Energy Exascale Earth System Model (E3SM) by developing a predictive understanding of Arctic tundra ecosystems underlain by permafrost and to quantify feedbacks from the Arctic tundra to the Earth system. NGEE Arctic is supported by the Department of Energy's Office of Biological and Environmental Research.Over Phases 1–3, observations made by the NGEE Arctic team across a gradient of permafrost landscapes in Arctic Alaska improved the representation of tundra processes in the land surface component of E3SM (the E3SM Land Model, ELM). Model improvements emphasized unique aspects of permafrost environments and explored reductions in model complexity while retaining predictive power. The Arctic-informed ELM developed by NGEE Arctic has been used to make novel predictions on processes ranging from permafrost thaw to soil biogeochemical cycling to Earth system feedbacks associated with the unique characteristics of tundra plants. In Phase 4, the NGEE Arctic team is evaluating our new predictive understanding under novel conditions across the Arctic domain. In collaboration with partners at long-term pan-Arctic research sites we are examining whether an Arctic-informed ELM can faithfully simulate interactions among surface and subsurface processes at site, regional, and pan-Arctic scales. In turn, we are using variety of tools to dynamically extend and evaluate ELM inference, with an emphasis on data synthesis and pan-Arctic model evaluation, reintegration of code with an evolving E3SM, scaling across heterogeneous Arctic landscapes, and the appropriate representation of the impacts of increasingly frequent Arctic disturbances.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Modeling of Macroscopic Dynamics in Three-Dimensional Plasma Configurations (Final Technical Report)

A numerical simulation-based study to investigate the stability properties of inherently three-dimensional plasma magnetic configurations is presented. These configurations, called stellarators and heliotrons, use external helical coils to provide the twisting magnetic field that is required to confine ionized plasma particles. In contrast to the tokamak configuration, which nominally has the property of axisymmetry, electrical current does not need to be induced in stellarator and heliotron plasma. Experiments have shown that these three-dimensional configurations are remarkably robust to macroscopic instability, even when driven to relatively high levels of plasma pressure. Our understanding of how the plasma is able to withstand macroscopic instability will be enhanced by the work proposed here. The effort will generalize existing numerical modeling capability, and it will apply the new capability to model macroscopic dynamics in these three-dimensional configurations. The development work has been applied to the NIMROD simulation code (https://nimrodteam.org) that was originally based on axisymmetric problem domains. It has been generalized to accommodate the three-dimensional shaping of the plasma volume in stellarators and heliotrons. Its use of equilibrium data has also been generalized, and the implementation allows different numerical representation of the magnetic field. Analysis has been performed to inform the selection of numerical methods, and a method to accelerate the iterative solution of the necessary algebraic systems has been developed and implemented. Verification of the new NIMSTELL code on benchmark problems has also been an important part of the effort.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Comparison of Global Aboveground Biomass Estimates From Satellite Observations and Dynamic Global Vegetation Models

The global forest carbon stocks represent the amount of carbon stored in woody vegetation and are important for quantifying the ability of the global forests to sequester atmospheric CO 2 and to provide ecosystem services (e.g., timber) under climate change. The forest ecosystem carbon pool estimates are highly variable and poorly quantified in areas lacking forest inventory estimates. Here, we compare and analyze aboveground biomass (AGB) estimates from five satellite-based global data sets and nine dynamic global vegetation models (DVGMs). We find that across the data sets, mean AGB exhibits the largest variability around the tropical area. In addition, AGB shows a similar latitudinal trend but large variability among the data sets. Satellite-based AGB estimates are lower than those simulated by DVGMs. The divergence among the satellite-based AGB estimates can be driven by the methodology, input satellite products, and the forested areas used to estimate AGB. The modeled NPP, autotrophic respiration, and carbon allocation mostly drive the variability of AGB simulated by DGVMs. The future availability of a high-quality global forest area map is anticipated to improve AGB estimate accuracy and to reduce the discrepancies among different satellite- and model-based AGB estimates. Furthermore, we suggest the carbon-modeling community reexamine the methodology used to estimate AGB and forested areas for a more robust global forest carbon stock estimation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Dynamic Line Rating Models and Their Potential for a Cost‐Effective Transition to Carbon‐Neutral Power Systems

Most transmission system operators (TSOs) currently use seasonally steady-state models considering limiting weather conditions that serve as reference to compute the transmission capacity of overhead power lines. The use of dynamic line rating (DLR) models can avoid the construction of new lines, market splitting, false congestions, and the degradation of lines in a cost-effective way. DLR can also be used in the long run in grid extension and new power capacity planning. In the short run, it should be used to help operate power systems with congested lines. The operation of the power systems is planned to have the market trading into account; thus, it computes transactions hours ahead of real-time operation, using power flow forecasts affected by large errors. In the near future, within a “smart grid” environment, in real-time operation conditions, TSOs should be able to rapidly compute the capacity rating of overhead lines using DLR models and the most reliable weather information, forecasts, and line measurements, avoiding the current steady-state approach that, in many circumstances, assumes ampacities above the thermal limits of the lines. Here, this work presents a review of the line rating methodologies in several European countries and the United States. Furthermore, it presents the results of pilot projects and studies considering the application of DLR in overhead power lines, obtaining significant reductions in the congestion of internal networks and cross-border transmission lines.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2025

A theoretical model for sheath dynamics in the two-plasma mode of a cylindrical filament discharge

We present a refined understanding of cylindrical plasma-facing emitter cathodes by examining the limitations of the conventional space charge limited model, which overlooks ion trapping dynamics within the virtual cathode. While conventional theory distinguishes between temperature-limited and space charge limited regimes based on electron current constraints, our findings reveal that trapped ions can form a quasi-neutral “second plasma,” expanding upstream and interacting with the primary plasma. In cylindrical and other non-planar geometries, an expanding trapped-ions layer produces current enhancement by converting more of the emitted electron flux to passing flux. The current enhancement leads to a complex coupling between the trapped-ions plasma and upstream plasma through ionizing collisions in the upstream region and charge exchange collisions in the virtual cathode. The coupling results in mutual density enhancement and at the same time, antagonistic expansion dynamics between the two plasma layers. The process is understood using a new “aid-compete” model. Numerical solutions, validated against simulations, reveal the model's capacity to predict system evolution and complex discharge behaviors, albeit with qualitative approximations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Molecular Modeling and Molecular Dynamics Simulation of a Packed and Intact Bacterial Microcompartment

Bacterial microcompartments (BMCs) are protein-bound organelles found in some bacteria which encapsulate enzymes for enhanced catalytic activity. These compartments spatially sequester enzymes within semipermeable shell proteins and are packed full of enzyme cargoes and metabolites as they fulfill their function. Coupling together recent SAXS and proteomics work, it is possible to develop molecular models for these microcompartments and interrogate enzyme and metabolite dynamics within. Our primary goal of this study is to quantify the permeability of metabolite glyceraldehyde-3-phosphate (G3P) and dihydroxyacetone phosphate (DHAP) across the BMC shell through classical molecular dynamics simulation. The Haliangium ochraceum model of BMC shell (PDB: 6MZX) was used to model an intact BMC of approximately 10 million atoms. Working at this scale presented its own challenges in managing large data sets, with multiple challenges and hardware advances discussed that facilitated this work. Over approximately 750 ns of aggregate simulation, we see multiple permeation events for these metabolites that were added at high concentration through the pores present within BMC shell tiles. When compared to independent permeability estimates for the same metabolites determined through replica exchange umbrella sampling simulations, the permeabilities varied by approximately 3 orders of magnitude. Regardless, the permeability coefficients for both G3P and DHAP are highly similar and very high, such that only very small concentration gradients can be maintained across the BMC shell between the cytosol and BMC interior. The large simulation systems also facilitated comparisons for molecular diffusivity in the crowded environment within the BMC shell. By our estimates, the viscosity within a packed BMC shell is at least 10-fold higher than it would be in neat solution and is the real driver for varying permeability estimates we obtained through simulation. These findings will be used as design inputs for future bioengineering efforts to make products from BMCs, highlighting how permeable BMC shells can be.

Diffusion↗
DOE OSTI2025

Global warming potential estimates of mass timber constructions beyond the first life: A dynamic radiative forcing modeling approach

In the face of a warming planet, steps must be taken to reduce the greenhouse gas emissions (GHG) associated with our building industry, which is a significant contributor to global emissions. Large, prefabricated wood elements such as mass timber panels (MTP) have great potential to achieve these reductions as they help displace high-embodied‑carbon materials like concrete and steel. Furthermore, storing the wood's biogenic carbon in buildings benefits the climate because it delays the eventual release of the carbon into the atmosphere. While these climate impacts have been assessed for the construction phase of mass timber buildings, relatively few life cycle assessment (LCA) studies have evaluated the climate impacts for the buildings' end-of-life (EOL) phase. This research estimates the climate impacts of four EOL scenarios for MTP: reusing as MTP, recycling into particleboard, incinerating, and landfilling. Using dynamic radiative forcing modeling and factoring in temporal GHG emissions and biogenic carbon storage, the global warming potential impacts are calculated for construction, deconstruction, and EOL processing of hybrid mass timber buildings in the U.S. Pacific Northwest for 160 years (GWP 160 ). The 160-year temporal scale used in this paper is an arbitrary scale, with the first 80 years being the assumed life of the building, followed by a series of reuse, recycle, or disposal scenarios over the second 80 years of that temporal scale. Of the four EOL scenarios considered in this paper, the ‘reuse’ scenario has the lowest net GWP 160 impact (calculated by summing the GWP 160 and carbon storage benefits, i.e., GWP$^{bioCS}_{160}$ ), emerging as a climate-preferred scenario, followed by ‘landfill’, ‘incinerate’, and ‘recycle’ scenarios. The lower net GWP 160 impact associated with the reuse scenario is due to the low fossil carbon emissions during EOL processing, as well as the biogenic carbon storage benefits. The results of this study also highlight the importance of efficient reuse and recycling strategies for wood in MTP.

42 ENGINEERING↗
DOE OSTI2024

Observational benchmarks inform representation of soil organic carbon dynamics in land surface models

Abstract. Representing soil organic carbon (SOC) dynamics in Earth system models (ESMs) is a key source of uncertainty in predicting carbon–climate feedbacks. Machine learning models can help identify dominant environmental controllers and establish their functional relationships with SOC stocks. The resulting knowledge can be integrated into ESMs to reduce uncertainty and improve predictions of SOC dynamics over space and time. In this study, we used a large number of SOC field observations (n=54 000), geospatial datasets of environmental factors (n=46), and two machine learning approaches (namely random forest, RF, and generalized additive modeling, GAM) to (1) identify dominant environmental controllers of global and biome-specific SOC stocks, (2) derive functional relationships between environmental controllers and SOC stocks, and (3) compare the identified environmental controllers and predictive relationships with those in models used in Phase 6 of the Coupled Model Intercomparison Project (CMIP6). Our results showed that the diurnal temperature, drought index, cation exchange capacity, and precipitation were important observed environmental predictors of global SOC stocks. While the RF model identified 14 environmental factors that describe climatic, vegetation, and edaphic conditions as important predictors of global SOC stocks (R2=0.61, RMSE = 0.46 kg m−2), current ESMs oversimplify the relationships between environmental factors and SOC, with precipitation, temperature, and net primary productivity explaining > 96 % of the variability in ESM-modeled SOC stocks. Further, our study revealed notable disparities among the functional relationships between environmental factors and SOC stocks simulated by ESMs compared with observed relationships. To improve SOC representations in ESMs, it is imperative to incorporate additional environmental controls, such as the cation exchange capacity, and refine the functional relationships to align more closely with observations.

54 ENVIRONMENTAL SCIENCES↗