Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Control variates”

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 865 records · Page 48

Dynamic phase transformations in additively manufactured Ti-6Al-4V during thermo-mechanical gyrations

A complex interaction of process parameters, geometry and scan strategies in Additive Manufacturing (AM), can bring about spatial and temporal transients, i.e., Σ T ( x, y, z, time ), within a part. Published literature focusses on fluctuating thermal cycles on the microstructure evolution. However, the microstructural variations have not been correlated to dynamic flow behavior due to the macro- and micro-scale phenomena, i.e., accumulated plastic strains brought about by large thermal gradients, transformational strains and crystallographic misfit strains. Therefore, here we studied the mechanical response of Ti6Al4V alloys produced by AM under externally imposed controlled thermo-mechanical reversals in a Gleeble® thermo-mechanical simulator. The stress-strain behaviors were correlated to phase fractions, lattice strains, and also limited information on crystallographic texture using neutron diffraction techniques at the VULCAN Beamline at SNS, ORNL and also metallographic studies. The results are discussed and rationalized based on theories of static and dynamic phase transformations.

36 MATERIALS SCIENCE↗

Nodal superconductivity in miassite Rh 17 S 15

Solid state chemistry has produced a plethora of materials with properties not found in nature. For example, high-temperature superconductivity in cuprates is drastically different from the superconductivity of naturally occurring metals and alloys and is frequently referred to as unconventional. Unconventional superconductivity is also found in other synthetic compounds, such as iron-based and heavy-fermion superconductors. Here, we report compelling evidence of unconventional nodal superconductivity in synthetic samples of Rh 17 S 15 (T c = 5.4 K), which is also found in nature as the mineral miassite. We investigated the temperature-dependent variation of the London penetration depth Δλ(T) and the disorder evolution of the critical superconducting temperature T c and the upper critical field H c2 (T) in single crystalline Rh 17 S 15 . We found a T - linear temperature variation of Δλ(T) below 0.3T c , which is consistent with the presence of nodal lines in the superconducting gap of Rh 17 S 15 . The nodal character of the superconducting state is supported by the observed suppression of T c and H c2 (T) in samples with a controlled level of non-magnetic disorder introduced by 2.5 MeV electron irradiation. We propose a nodal sign-changing superconducting gap in the A 1g irreducible representation, which preserves the cubic symmetry of the crystal and is in excellent agreement with the derived superfluid density. To the best of our knowledge, this establishes miassite as the only mineral known so far that reveals unconventional superconductivity in its clean synthetic form, though it is unlikely that it is present in natural crystals because of unavoidable impurities that quickly destroy nodal superconductivity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Gyrokinetic simulations of the effects of magnetic islands on microturbulence in KSTAR

Gyrokinetic simulations are utilized to study effects of magnetic islands (MIs) on the ion temperature gradient (ITG) turbulence in the KSTAR tokamak with resonant magnetic perturbations. Simulations show that the transport is controlled by the nonlinear interactions between the ITG turbulence and self-generated vortex flows and zonal flows, leading to an anisotropic structure of fluctuation and transport on the poloidal plane and in the toroidal direction. MIs greatly enhance turbulent transport of both particle and heat. The turbulent transport exhibits variations in the toroidal direction, with transport through the resonant layer near the island X-point being enhanced when the X-point is located at the outer mid-plane. A quantitative agreement is shown between simulations and KSTAR experiments in terms of time frequency and perpendicular wavevector spectrum.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Latent Twins

Over the past decade, scientific machine learning has transformed the development of mathematical and computational frameworks for analyzing, modeling, and predicting complex systems. From inverse problems to numerical partial differential equations (PDEs), dynamical systems, and model reduction, these advances have pushed the boundaries of what can be simulated. Yet they have often progressed in parallel, with representation learning and algorithmic solution methods evolving largely as separate pipelines. With Latent Twins, we propose a unifying mathematical framework that creates a hidden surrogate in latent space for the underlying equations. Whereas digital twins mirror physical systems in the digital world, Latent Twins mirror mathematical systems in a learned latent space governed by operators. Through this lens, classical modeling, inversion, model reduction, and operator approximation all emerge as special cases of a single principle. We establish the fundamental approximation properties of Latent Twins for both ordinary differential equations (ODEs) and PDEs and demonstrate the framework across three representative settings: (i) canonical ODEs, capturing diverse dynamical regimes; (ii) a PDE benchmark using the shallow-water equations, contrasting Latent Twin simulations with deep operator network and forecasts with a four-dimensional variational method baseline; and (iii) a challenging real-data geopotential reanalysis dataset, reconstructing and forecasting from sparse, noisy observations. Latent Twins provide a compact, interpretable surrogate for solution operators that evaluate across arbitrary time gaps in a single-shot, while remaining compatible with scientific pipelines such as assimilation, control, and uncertainty quantification. Looking forward, this framework offers scalable, theory-grounded surrogates that bridge data-driven representation learning and classical scientific modeling across disciplines.

Latent Twins↗

Heterosynaptic plasticity in memristive and memcapacitive lipid bilayers: A snapshot review

Synaptic plasticity refers to activity-dependent synaptic strengthening or weakening between neurons. It is usually associated with homosynaptic plasticity, which refers to a synaptic junction controlled by interactions between specific neurons. Heterosynaptic plasticity, on the other hand, lacks this specificity. It involves much larger populations of synapses and neurons and can be associated with changes in synaptic strength due to nonlocal alterations in the ambient electrochemical environment. Here, this paper presents specific examples demonstrating how variations in the ambient electrochemical environment of lipid membranes can impact the nonlinear dynamical behaviors of memristive and memcapacitive systems in droplet interface bilayers (DIBs). Examples include the use of pH as a modulatory factor that alters the voltage-dependent memristive behavior of alamethicin ion channels in DIB lipid bilayers, and the discovery of long-term potentiation (LTP) in a lipid bilayer-only system after application of electrical stimulation protocols.

36 MATERIALS SCIENCE↗

Genetics and Genomics of Pathogen Resistance in Switchgrass (Final Report)

This project was funded by DOE under Grant no. DE-SC0016108. Originally approved for the 2016-2019 period, two no-cost extensions were solicited and approved, which prolonged the lifespan through July 2021. This final report informs on the results obtained so far from the research implemented. The research hinged on integrating genomics (genomic selection, RNAseq, virus-plant interactions) with classical genetics (conventional breeding) to incorporate durable resistance to fungal (rust) and viral (mosaic) diseases in switchgrass (Panicum virgatum) populations being bred for bioenergy. Higher biomass yield, higher quality (low lignin content), and durable disease resistance are key features to make lignocellulosic switchgrass feedstocks economically competitive and sustainable. Genomic selection is being applied on three generations of a switchgrass population derived from crossing two ecotypes (Kanlow as lowland female and Summer as upland male) with differential performance in terms of biomass yield and quality, disease resistance, and winter survivability. Target populations were screened for rust and mosaic in field and/or lab and phenotyped for biomass yield and quality traits. Genetic analyses were applied across generations to capture the joint inheritance of the targeted traits and predict breeding values for parents and progeny with greater accuracy. Parental and a panel of different switchgrass populations were genotyped with the DArTseq technology to develop SNP (0, 1, 2) and in-silico (presence/absence) DArT markers. Rust inoculations techniques were developed and applied successfully on switchgrass. The original populations (Kanlow and Summer) were sequenced with RNAseq to capture the gene expression profiles across sequential time-points and appraise the basis of greater resistance in the Kanlow vs the Summer ecotype. Constructs of PMV and sPMV mosaic virus were assembled and tested first on proso millet to find the best protocol to use later on switchgrass. Results from the preliminary analyses indicate that 1) ample additive genetic variation is available for selection and improving this inter-ecotypic population for yield, quality, and disease traits, 2) significant gains are to be expected with the genetic correlations being favorable between yield and lignin content and between yield and disease ratings, 3) substantial differences exist in the genetic regions controlling rust resistance in the two ecotypes, 4) co-infection with PMV isolates from Nebraska and its satellite from Kansas elicit severe mosaic symptoms, and 5) two different genetic systems are responsible for imparting resistance to rust and virus in switchgrass.

59 BASIC BIOLOGICAL SCIENCES↗

Analysis of Warped April Tag Impacts on Detection and Pose Estimation

This report evaluates the impact of geometric deformation on an April Tag, particularly when warped due to attachment on a curved surface, on its detectability and pose estimation performance. A comparative analysis was conducted using a flat April Tag as a control under identical experimental conditions, which involved recording video sequences with varying viewing angles. For detectability, the warped tag exhibited consistent detection failures at viewing angles beyond 40° and complete failures beyond 60°, whereas the flat tag maintained reliable detection across all angles. For pose estimation, measured by pose jitter (variation in rotation and translation), differences between the warped and flat tags were minimal and statistically insignificant, indicating robust performance even for the deformed tag. These findings suggest that while geometric warping reduces an April Tag’s detectability, its pose estimation accuracy remains relatively unaffected under the tested conditions.

42 ENGINEERING↗

SDSS-IV MaNGA: The Effect of Stellar Mass and Halo Mass on the Assembly Histories of Satellite Galaxies

We combine an unprecedented MaNGA sample of over 3000 passive galaxies in the stellar mass range 10 9 –10 12 M ⊙ with the Sloan Digital Sky Survey group catalog by Tinker to quantify how central and satellite formation, quantified by radial profiles in stellar age, [Fe/H], and [Mg/Fe], depends on the stellar mass of the galaxy (M * ) and the mass of the host halo (M h ). After controlling for M * and M h , the stacked spectra of centrals and satellites beyond the effective radius (r e ) show small, yet significant differences in multiple spectral features at the 1% level. According to spectral fitting with the code alf , a primary driver of these differences appears to be [Mg/Fe] variations, suggesting that stellar populations in the outskirts of satellites formed more rapidly than the outer populations of centrals. To probe the physical mechanisms that may be responsible for this signal, we examined how satellite stellar populations depend on M h . We find that satellites in high-M h halos show older stellar ages, lower [Fe/H], and higher [Mg/Fe] compared to satellites in low-M h halos, especially for M * = 10 9.5 –10 10.5 M ⊙ . These signals lend support to environmentally driven processes that quench satellite galaxies, although variations in the merger histories of central and satellite galaxies also emerge as a viable explanation.

79 ASTRONOMY AND ASTROPHYSICS↗

Optimizing power system restoration with damaged communications

Utility procedures for power system blackstart and restoration typically assume that energization decisions can be reliably communicated across the grid. In reality, the communications and control network would likely also be affected in power outages, such as those caused by extreme weather events or cyber-attacks. This paper studies the effect of damage to the power system communications and control infrastructure on restoration operations following a blackout. We model the communications infrastructure as a graph, overlaying the power grid, and imposing the requirement that every energized element in the power grid be observable from a control center. We expand on a specialized branch-and-bound algorithm from the literature to optimize the restoration process and devise an initialization heuristic and a rounding heuristic to improve solution speed. We perform numerical experiments on synthetic systems for Illinois and Texas with outages based on a solar flare or hurricane. We compare the results of our specialized branch-and-bound algorithm to the results from (i) the initialization heuristic alone, (ii) a variation of this heuristic that we use as a baseline, and (iii) the restoration optimization for the power system without communications constraints. Here, we find that damage to the communications infrastructure significantly increases the time required to re-energize the grid. Moreover, by simultaneously optimizing communications repairs and grid energization decisions, we are able to re-energize the grid significantly faster than if communications repairs and energization decisions were made independently or with partial coordination, motivating improvements to current industry practice.

97 MATHEMATICS AND COMPUTING↗

Correlating Surface Processing of Nb Superconducting RF Cavities with the Evolution of Surface Electronic States

Superconducting-radio frequency (RF) cavities provide an efficient way to accelerate particle beams with extremely high acceleration gradients while generating very small power dissipation. The few nanometers of the surface play a critical role in defining the RF performance of superconducting Nb based cavities. Over the past two decades, several pioneering surface treatment and processing methods have emerged, enabling remarkable improvements in superconducting cavity performance by simultaneously achieving high quality factors with increasing maximum acceleration gradients. These processing approaches include chemical polishing, distinct multi-step thermal treatments under ultra-high vacuum (UHV) conditions over low to high temperature regimes, as well as high-temperature treatments under controlled nitrogen atmospheres. Beyond their macroscopic impact on RF performance, these methods produce distinct surface oxide configurations characterized by different valence states, oxide thicknesses, chemical uniformity, and oxygen concentration profiles extending into the near-surface bulk of niobium. In this work, we are trying to understand how the surface-processing methods and the resulting oxide/oxygen profiles affect the electronic structure of surface and the mechanism of superconductivity. Using a combination of X-ray photoemission and X-ray absorption spectroscopies, we investigate how the valence-band structure and the electronic density of states (DoS) near the Fermi level evolve under different surface treatments. By employing tunable photon energies across multiple elemental absorption edges, we use resonant photoemission to disentangle and identify the elemental contributions to specific valence-band features. Our observations show that different surface-processing methods lead to distinct temperature evolutions of the DoS and valence-band states near the Fermi level. Our results suggest variations in different Nb-O orbital hybridizations in distinct processes and point towards the possibility of different underlying mechanisms of superconductivity governed by surface chemistry and oxide configuration. We also correlate these distinct superconducting mechanisms with RF cavity performance, specifically focusing on measured surface resistance, the nature of the Q-slope, and quench fields observed in superconducting RF measurements. These results will enable us to identify the potential limiting factors and relevant controllable parameters that can be further optimized to improve the performance of superconducting RF cavities.

Tripathi, Malvika [Fermilab] (ORCID:00000001989251↗

Synthetic Calcium Carbonate Production by Carbon Dioxide (CO 2 ) Mineralization of Industrial Waste Brines

The global scale of CO 2 emission has crossed 36 B tons, and the United State represents 14% of the total emissions. In light of the high cost associated with current CO 2 capture processes ($60-to-150 per ton of CO 2 ), carbon capture and utilization (CCU), wherein CO 2 is converted to beneficial products, provides a pragmatic path to overcome the economic barrier for CO 2 emissions control. In particular, CO 2 mineralization offers an attractive route as it produces high-value mineral carbonates that sequester CO 2 in a stable form. For instance, fine precipitated calcium carbonate, a valuable product with unit price in the range of $230-280/t and a global market projected to reach 99 M tons in 2020, can be produced by capturing CO 2 within aqueous Ca 2+ solution. However, two critical challenges need to be overcome—(i) the need for costly processes such as electrolysis or addition of alkali hydroxides to maintain alkalinity during mineralization, and (ii) the geographic availability of Ca-rich solutions (or brines) which often renders efficient integration with power plants impractical. To simultaneously address these challenges, this project developed two novel routes to enable carbonate production based on efficient CO 2 mineralization. In one variation, the project developed a new alkaline carbonation process to capture CO 2 and produce calcite precipitates from coal ashes. Herein, coal ash is carbonated first in a sodium carbonate solution. The carbonation reaction produces sodium hydroxide and raises the solution pH. This high-pH hydroxide solution is then used for CO 2 capture, which converts the sodium hydroxide solution back to sodium carbonate to repeat the carbonation cycle. The carbonated ash residue is refined with a CO 2 pressure swing step to produce high-purity precipitated calcium carbonate. In another variation, Ca-rich produced water brine serves as the Ca-source to mineralize CO 2 . An H + /Na + ion-exchange cycle was designed to provision alkalinity during mineralization and regenerate the ion-exchange reagent in brines with high salinity. Ca-depleted brines can then be treated within a centralized water treatment facility. Both variations beneficially utilize reject streams—such as coal ashes and brines from oil and gas extraction or CO 2 storage operations—that are available at substantial quantities in the vicinity of coal power plants within the U.S. In addition to technology development, this project included techno-economic and life-cycle analysis to identify technically and economically appropriate solutions for power plants in different geographic locations. Taken together, the project provides a unique route to integrate CO 2 emissions control and waste handling/treatment for coal power plants, while producing a high-value product to offset the economic burden associated with waste management. The developed processes utilize post-desulfurization flue gas from coal-fired power plants as is. The CO 2 conversion reactions are performed in alkaline brines at ambient pressure from the flue gas, thereby minimizing the energy burden. By beneficiation of industrial waste streams with low energy input, the processes offers significant technical advantages in energy and CO 2 footprint over the current paradigm of precipitated calcium carbonate production, which involves calcination of limestone at processing temperatures in excess of 800 °C. Overall, this project developed viable CO 2 mineralization processes while strategically maximizing the economic benefit.

01 COAL, LIGNITE, AND PEAT↗

Tailoring Cu-Zr gradient nanoglass structures: Influence of nanoparticle size and cooling rates on glass-glass interfaces

The study of gradient nanoglasses (GNGs) has gained attention due to their unique mechanical properties and potential applications in advanced materials. This study employs molecular dynamics simulations to synthesize a GNG using Cu-Zr metallic glass nanoparticles (NPs) sized from 3 to 15 nm. The NPs were produced by melting and quenching metallic clusters at a relatively slow quench rate of 10 9 K/s. The synthesis of GNG is elucidated along with the characterization of its heterogeneous metallic glass nanostructure. A seamless GNG structure is formed through cold compression of Cu 64 Zr 36 amorphous NPs of varying sizes. The influence of NP size on the GNG structure is investigated, utilizing deeply relaxed NPs, which exhibit a characteristic Cu segregation pattern on their surfaces. The results highlight an increase in structural heterogeneity due to heterogeneous mass transport and the development of local composition and density variations caused by Cu segregation at glass-glass interfaces (GGIs). A reduction in NP size is correlated with decreased Cu atomic displacements and local density at GGIs, suggesting that larger NPs may produce stronger GGIs. This research presents a novel methodology for synthesizing heterogeneous metallic glasses, demonstrating the capacity to control and customize nanostructure heterogeneity through the manipulation of NP sizes and cooling rates. Furthermore, these findings enhance our understanding of structural evolution during nanoglass synthesis and lay the foundation for further exploration in nanomaterial synthesis and characterization.

36 MATERIALS SCIENCE↗

Impact of Irradiation-Induced Filter Heating on Calibration of NIR-Longpass-Filtered Reference Solar Cells

Reference solar cells play a crucial role in determining the performance of photovoltaic (PV) devices. In the performance calibration chain of mainstream single-junction PV technologies, common reference cell types include Si, KG-filtered Si, and GaAs cells. For emerging multijunction (MJ) PV technologies, such as CdTe/Si and perovskite/Si, Si reference cells with colored glass near-infrared-longpass (LP) filters like Schott RG-715 and 850 glass have been proposed. They offer a better spectral response that matches the bottom junctions of the emerging MJs, which could lead to lower uncertainties in performance measurements. However, this article reveals a prominent decrease in short-circuit current (ISC) during National Renewable Energy Laboratory's primary calibration over the course of minutes when using this type of LP-filtered reference cell, which could result in unacceptable measurement errors. Unlike quartz or KG filter glasses, LP colored glass filters demonstrate temperature-sensitive cut-on wavelength. When incident irradiance reaches these LP-filtered reference cells, the increased temperature due to light absorption causes a shift of the cut-on to longer wavelength. As a result, the device I SC exhibits a continuous decrease (approximately a 2.5% drop with RG850 LP filter in a 5-min duration) even when the device temperature is controlled at 25 degrees C. To address the temperature gradient issue between colored glass LP filter and the reference cell, we propose the direct integration of thin-film LP semiconductor layer on reference cells. This type of LP-filter-integrated cell has shown minimal temperature-related I SC variation and can serve as a more reliable reference cell source for accurate performance measurements.

14 SOLAR ENERGY↗

Library of Advanced Models of large-scale PV (LAMP) (Final Technical Report)

A suite of dynamic models of photovoltaic (PV) systems is developed in this project to accurately capture PV system dynamics during different contingencies. The following models are included: (i) quasi-dynamic, (ii) high-fidelity electromagnetic transient (EMT), (iii) advanced, and (iv) baseline. The models have the flexibility to consider different control functionalities at different controller levels. The study’s requirements determine the type of model needed. For example, high-fidelity EMT models are recommended for contingency studies. Advanced models are recommended for online or real-time implementations with resource constraints, whereas quasi-dynamic models are useful for studying the impact of weather variations.

14 SOLAR ENERGY↗

Model form and sensitivity analysis of CALPHAD-based nucleation models in b-stabilized Ti alloys

Accurate prediction of α-phase nucleation and growth in β-stabilized titanium alloys is crucial for designing heat treatments to optimize mechanical properties in additively manufactured lightweight components. Ideally, predictions of nucleation and growth would incorporate both top-down observations of past experimental heat treatments and bottom-up modeling of phase transformations; however, the appropriate method of combining these information sources is not self-evident. Combining top-down and bottom-up information requires a unified form of model that can connect between spatiotemporal scales, as well as sets of fitting parameters that can be identified by each data source. The selection of which parameters to fit to which data source can be made based on expert opinion, or by performing a sensitivity analysis. In solid-solid nucleation, direct observation of the nucleation and growth process is challenging. Most data on the heat treatment-controlled phase transformations are not in-situ. To predict the process and outcome of the nucleation, growth and coarsening of precipitates, theoretical models of the nucleation pathway are used to bridge the gap. Many sources of uncertainty affect the modeling of this nucleation process. It can be influenced by small variations in the thermomechanical processing history, chemical composition, and initial microstructure. If molecular dynamics (MD) simulations are used to determine thermodynamic quantities and inform CALPHAD modeling, additional uncertainty can be introduced and accounted for using Bayesian methods. Top-down uncertainties require additional steps to quantify. The influence of nucleation model form on the sensitivity of predictions to input parameters and physical conditions is the focus of this study. Classical nucleation theory (CNT) allows modeling to formulate the nucleation as homogeneous or, more commonly, heterogeneous. Non-classical nucleation models are also increasingly explored as a means of reconciling top-down and bottom-up data. In this study, the sensitivity of the intragranular nucleation of α in a β-annealed, slow-cooled aging (BASCA) heat treatment of β-stabilized Ti5553 alloy is explored using CNT and both heterogeneous and homogeneous assumptions. The Kampmann-Wagner Numerical model of precipitate nucleation and growth is employed. Using open-source tools (pyCalphad and thermodynamic modeling of TiMo as a surrogate system, a sensitivity analysis is performed to measure variations in key parameters, including chemical driving force, interfacial energy, and diffusivity, as they relate to predictions of precipitate number density. The inclusion of top-down and bottom-up data in selection of nucleation model form is discussed.

Rodriguez Negron, A. M.↗

Measurement of the QE Map

113 MHz SRF gun utilized CsK 2 Sb and Na 2 KSb cathodes, which were irradiated with a pulsed green laser with a wavelength of 0.53 microns. The drive laser generated 78 kHz pulses with an acousto-optic modulator selecting the desired number of pulses reaching the cathode. The operator can select the aperture size, defining the laser spot diameter on the cathode, and the power level, which is controlled by a rotating polarizer. The laser power reaching the cathode can be measured with a meter based on a photodiode. The quantum efficiency (QE) of a photocathode serves as a primary indicator of its state. Small QE renders the cathode unusable when the desired bunch charge cannot be achieved. Strong variations in the quantum efficiency on the cathode surface might lead to the growth of the beam emittance due to the space charge effects. We have created a routine to measure the distribution of the QE on the cathode. The operator selects the aperture size, the central spot location, and the scan ranges for both planes. The scan starts from the central spots and proceeds along a spiral-like trajectory, making N rounds. The number of rounds can be chosen, and the total number of measured points is (2N+1) 2 . The laser spot movement is shown in Fig. 1.

43 PARTICLE ACCELERATORS↗

Twinac: A Universal Framework for Virtual Accelerator Controls

We propose a universal framework for a system of virtual accelerator controls (Twinac), a standard toolkit for research institutions to design, maintain, and use a real-time, end-to-end "digital twin" of their particle accelerator facility. This virtual counterpart will mirror any physical accelerator to provide (1) predictive maintenance; (2) surveillance of hidden environmental factors, such as seasonal temperature variations, which could impact performance of power supplies, magnets, and other instruments; and (3) a capability to model novel ways to operate the accelerator without risking equipment damage. Twinac is envisioned as accelerator facility agnostic, allowing institutions to share and reuse myriad simulation approaches (analytics-based, A.I.-driven, or combinations of the like) across facilities. The Twinac system lays the groundwork for a collaborative network of institutions to maintain and update this shared virtual accelerator technology.

Miceli, Tia M. [Fermilab] (ORCID:0000000265577789)↗

Divergent controls of soil organic carbon between observations and process-based models

Abstract The storage and cycling of soil organic carbon (SOC) are governed by multiple co-varying factors, including climate, plant productivity, edaphic properties, and disturbance history. Yet, it remains unclear which of these factors are the dominant predictors of observed SOC stocks, globally and within biomes, and how the role of these predictors varies between observations and process-based models. Here we use global observations and an ensemble of soil biogeochemical models to quantify the emergent importance of key state factors – namely, mean annual temperature, net primary productivity, and soil mineralogy – in explaining biome- to global-scale variation in SOC stocks. We use a machine-learning approach to disentangle the role of covariates and elucidate individual relationships with SOC, without imposing expected relationships a priori . While we observe qualitatively similar relationships between SOC and covariates in observations and models, the magnitude and degree of non-linearity vary substantially among the models and observations. Models appear to overemphasize the importance of temperature and primary productivity (especially in forests and herbaceous biomes, respectively), while observations suggest a greater relative importance of soil minerals. This mismatch is also evident globally. However, we observe agreement between observations and model outputs in select individual biomes – namely, temperate deciduous forests and grasslands, which both show stronger relationships of SOC stocks with temperature and productivity, respectively. This approach highlights biomes with the largest uncertainty and mismatch with observations for targeted model improvements. Understanding the role of dominant SOC controls, and the discrepancies between models and observations, globally and across biomes, is essential for improving and validating process representations in soil and ecosystem models for projections under novel future conditions.

58 GEOSCIENCES↗