Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “state dependent variable”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Optimal energy storage portfolio for high and ultrahigh carbon-free and renewable power systems

Achieving 100% carbon-free or renewable power systems can be facilitated by the deployment of energy storage technologies at all timescales, including short-duration, long-duration, and seasonal scales; however, most current literature focuses on cost assessments of energy storage for a given timescale or type of technology. In this work, we use an optimization framework with high spatial and temporal resolution to simultaneously assess the variable renewable power deployment and the optimal storage portfolio for seven independent system operators in the United States. Results indicate that achieving high (75–90%) and ultrahigh (>90%) energy mixes requires combining several flexibility options, including renewable curtailment, short-duration, long-duration, and seasonal storage. For instance, carbon-free and renewable energy mix targets of up to 80% are achieved with economic curtailment and a combination of short- and long-duration energy storage for the performance and cost assumptions used. After that, there is a point between 80% and 95% where seasonal storage becomes cost-competitive, depending on the specific power system. Moreover, our results indicate that storage-to-storage operation—one storage device used to charge another storage device—and the decoupling of charging and discharging storage power capacity are cost-effective options for the integration of high and ultrahigh shares of carbon-free or renewable power sources. Additionally, the results from this study show that an 85% carbon-free or renewable energy mix can be achieved at a cost of avoided CO 2 emissions of US$66.0 per tonne or less, regardless of the power system.

25 ENERGY STORAGE↗
DOE OSTI2024

Modeling Oxygen Partial Pressure in Solid Oxide Electrolysis Cells: The Microstructure Effect

Oxygen partial pressure is an important thermodynamic state variable that affects both the performance and degradation of solid oxide electrolysis cells. In this work, a 3D model developed from Virkar’s 1D model has been applied to reconstructed and synthetic microstructures of Ni-YSZ-GDC-LSCF cell. The effect of the microstructures, including the thickness of the YSZ and GDC layer, and the compositions of the hydrogen and oxygen electrode, on the distribution of oxygen partial pressure has been investigated. The results show that the maximum oxygen partial pressure may occur on the interface between oxygen electrode and GDC layer or between GDC layer and YSZ layer depending on the rate of oxygen ion exchange between GDC and YSZ. A thicker GDC layer lowers the maximum oxygen partial pressure in the cell, while a thicker YSZ layer lowers the maximum oxygen partial pressure in the hydrogen electrode. In addition, the Ni:YSZ ratio and porosity also affects the maximum partial pressure. These findings provide insights on mitigating degradation in solid oxide electrolysis cell by tuning the microstructures.

Lei, Yinkai↗
DOE OSTI2020

Dominant effect of anharmonicity on the equation of state and thermal conductivity of MgO under extreme conditions

Understanding heat transfer in the earth mantle region is of great scientific interest but highly challenging due to dissimilar conditions from our daily lives. Thermal conductivity of highly pressured and thermally activated materials can be completely different from that at ambient condition. In this study we calculate the equation of state and thermal conductivity of cubic crystalline MgO, a major component of the mantle, in such an environment. It is shown that the material properties are not accurately captured unless rigorous treatment of anharmonicity of phonon vibrational motion is explicitly applied. To predict renormalized phonon dispersion relations, interatomic force constants of high-order expansion are required, especially at elevated temperature. The anharmonicity plays an increasingly substantial role for regulating volume expansion as temperature rises, because the Gibbs free energy of MgO considerably deviates from what was predicted by the quasiharmonic approximation. We find that the lattice thermal conductivity of MgO and the heat transfer of the earth's mantle are reliably estimated only when both the temperature-dependent renormalization of phonon frequencies and four-phonon scattering effects are taken into consideration. We map the characteristics at Geotherm into contour diagrams in which equilibrium volumes and the thermal conductivities are described as functions of temperature and pressure. Furthermore, the thermal conductivity values of MgO in this study are smaller than those in previous theoretical studies due to the additional incorporation of four-phonon scatterings in the calculations. Specifically, deviation between the values in previous studies and in this study becomes larger as thermodynamic variables reach those near the core-mantle boundary. For instance, thermal conductivities of 52.15 and 23.67 Wm -1 K -1 are estimated at ambient and extreme conditions (T = 4000K and P = 135 GPa), respectively.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Microstructure-sensitive mechanical behavior of an additively manufactured psuedoelastic shape memory alloy

The additive manufacturing of shape memory alloys into complex geometries enables fabrication of advanced functional systems across a variety of fields and domains. This work presents results focused on the mechanical behavior of additively manufactured shape memory pseudoelastic NiTi. The deformation induced solid state phase transformation from austenite to martensite allows this system to accommodate large recoverable strains. This deformation behavior is fundamentally driven by crystal-scale transformation physics. Laser powder bed fusion processing reveals that the resulting microstructure, both grain morphology and crystallographic texture, is strongly dependent on the manufacturing processing history. Exhaustive mechanical testing demonstrates that these microstructural factors strongly impact both tensile and cyclic stress–strain behavior. Cyclic dissipative behavior, however, is similar across all tested microstructures following an initial transient period. Remarkably, analysis of spatial strain fields during tensile loading reveals two distinctly different localization “modes”. The first is initiation of localized deformation bands which continuously propagate through the tensile bar during loading. In the second mode localization is observed but lacks propagation; instead additional localization cites nucleate during subsequent loading. The latter phenomena is suspected to be driven by grain-scale deformation physics as the localized band morphologies coincide with grain morphologies. These phenomena strongly impact the resulting aggregate stress–strain behavior. Hence, manufacturers and designers of psuedoelastic functional components must at the very least consider the potential variability in properties when considering additive manufacturing processing. More ideally the process–structure–property relations can be used to further tailor and optimize final functional performance.

Additive manufacturing↗
DOE OSTI2022

Search for resonant production of strongly coupled dark matter in proton-proton collisions at 13 TeV

The first collider search for dark matter arising from a strongly coupled hidden sector is presented and uses a data sample corresponding to 138 fb$^{-1}$, collected with the CMS detector at the CERN LHC, at $ \sqrt{s} $ = 13 TeV. The hidden sector is hypothesized to couple to the standard model (SM) via a heavy leptophobic Z' mediator produced as a resonance in proton-proton collisions. The mediator decay results in two “semivisible” jets, containing both visible matter and invisible dark matter. The final state therefore includes moderate missing energy aligned with one of the jets, a signature ignored by most dark matter searches. No structure in the dijet transverse mass spectra compatible with the signal is observed. Assuming the Z' boson has a universal coupling of 0.25 to the SM quarks, an inclusive search, relevant to any model that exhibits this kinematic behavior, excludes mediator masses of 1.5–4.0 TeV at 95% confidence level, depending on the other signal model parameters. To enhance the sensitivity of the search for this particular class of hidden sector models, a boosted decision tree (BDT) is trained using jet substructure variables to distinguish between semivisible jets and SM jets from background processes. When the BDT is employed to identify each jet in the dijet system as semivisible, the mediator mass exclusion increases to 5.1 TeV, for wider ranges of the other signal model parameters. These limits exclude a wide range of strongly coupled hidden sector models for the first time.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Machine-learning-assisted deciphering of microstructural effects on ionic transport in composite materials: A case study of Li 7 La 3 Zr 2 O 12 -LiCoO 2

The effective diffusivity of ionic species in multiphase materials is critical for the design and function of composite materials for electrochemical energy storage. In practice, effective diffusivity depends sensitively not only on the intrinsic diffusivities of constituting materials but also on their topological arrangement; nevertheless, these coupled contributions are oversimplified in most analytical models. Here, we combine atomistically informed mesoscale modeling and machine learning (ML) analysis to unravel how such features affect effective diffusivity in two-phase composites. Using the Li 7 La 3 Zr 2 O 12 -LiCoO 2 composite solid-state battery cathode as a model system, we compute effective diffusivity for 600 distinct dense polycrystalline microstructures with different topological configurations of grains, grain boundaries, and heterointerfaces. We verify that in addition to atomic-scale variabilities, microstructural feature diversity can significantly impact effective transport properties. Across the ensemble of test microstructures, this often results in bimodal distributions of effective diffusivity that encompass two qualitatively distinct operating mechanisms, which we identify via flux analysis. An ML approach reveals that the most critical determining factors for effective diffusivity are the connectivity of bulk phases and their heterointerfaces. The role of ionic mobility at the heterointerfaces is also discussed. These insights highlight the combined importance of microstructure and interface engineering in tuning the transport properties of ionic species in composite materials. In conclusion, our framework can also be extended for understanding generic microstructure-property relationships in other complex multiphase materials.

25 ENERGY STORAGE↗
DOE OSTI2022

1 H chemical shift anisotropy: a high sensitivity solid-state NMR dynamics probe for surface studies?

We report dynamics play significant roles in chemistry and biochemistry-molecular motions impact both large- and small-scale chemical reactions in addition to biochemical processes. In many systems, including heterogeneous catalysts, the characterization of dynamics remains a challenge. The most common approaches involve the solid-state NMR measurement of anisotropic interactions, in particular 2 H quadrupolar coupling and 1 H–X dipolar coupling, which generally require isotope enrichment. Due to the high sensitivity of 1 H NMR, 1 H chemical shift anisotropy (CSA) is a particularly enticing, and underexplored, dynamics probe. We carried out 1 H CSA and 1 H– 13 C dipolar coupling measurements in a series of model supported complexes to understand how 1 H CSA can be leveraged to gain dynamic information for heterogeneous catalysts. Mathematical descriptions are given for the dynamic averaging of the CSA tensor, and its dependence on orientation and asymmetry. The variability of the orientation of the tensor in the molecular frame, in addition to its magnitude and asymmetry, negatively impacts attempts to extract quantitative dynamic information. Nevertheless, 1 H CSA measurements can reveal useful qualitative insights into the motions of a particularly dilute site, such as from a surface species.

1H chemical shift anisotropy↗
DOE OSTI2025

Turbulence statistical analysis of the L-H transition and RMPs in KSTAR

Here, we investigate the turbulence statistics associated with low-to-high confinement (L-H) transitions and externally applied resonant magnetic perturbations (RMPs) in KSTAR. Time-series fluctuations of electron density n e , electron temperature T e , and the time derivative of the poloidal magnetic field dB θ /dt (Mirnov coils) are analysed using information-geometric measures (information rate Γ and information length $\mathcal{L}$ = ∫ Γ dt), together with kurtosis κ and variance σ 2 . In low-density upper single-null plasmas (n e ~ 1.2 x 10 19 m -3 ), a ~80 kHz magnetic mode coupling n e , T e , dB θ /dt and emerges prior to the L-H transition and persists into the edge-localised modes H-mode. Edge-localised RMPs (ERMPs) suppress this coherent mode but enhance intermittency, producing frequent bursts that abruptly reshape the time-dependent probability density functions (PDFs) and generate large spikes in Γ (with smaller changes in κ), signalling ERMP-driven departures from quasi-stationarity. The impact of ERMPs on background fluctuation levels depends on density, radial location, and the fluctuating variable itself ($\tilde{n}$, $\tilde{T}$, $\dot{B}$ θ ), whereas $\mathcal{L}$ provides a robust, regime-agnostic measure of cumulative statistical reorganisation and spatial decorrelation. In particular, at low density we observe weaker coupling between $\tilde{n}$ and $\tilde{T}$, along with a tendency toward decreased radial correlation-most clearly for $\tilde{T}$-under ERMPs. Overall, information geometry cleanly captures intermittent events, quantifies non-equilibrium PDF evolution, and offers a compact, cross-diagnostic metric for assessing resonant magnetic perturbation effects on edge transport and correlation across densities, radial locations, and confinement states.

Kim, Eun-jin [Coventry Univ. (United Kingdom); Seo↗
DOE OSTI2026

Comparative Evaluation of Spectral Methods for Robust Reactor Noise Estimation

Reactor noise analysis provides a noninvasive means to determine neutron kinetic parameters from stochastic fluctuations in detector signals. However, standard cross-power spectral density (CPSD) analyses can be sensitive to numerical processing choices, which may introduce processing-dependent systematic shifts in estimates of the prompt neutron decay constant (α) and limit reproducibility. This study uses a hybrid multitaper–Welch spectral estimator to analyze subcritical noise measurements from a fast-spectrum critical assembly. The decay constant α was extracted using three frequency-domain methods: the CPSD, the magnitude-squared coherence (MSC), and the generalized magnitude-squared coherence (GMSC). These coherence-based estimators normalize detector auto-spectral structure and are expected to reduce the sensitivity of fitted α values to processing parameters. A Sobol global sensitivity analysis identified which numerical inputs most strongly influence the fitted values of α. All estimators produced a linear dependence of α on inverse count rate, with delayed-critical extrapolations near 1.7 × 10 4 s −1 , in agreement within 8% of MCNP6.3 KOPTS benchmark calculations. Sensitivity results show that while the CPSD depends on both time-bin width and taper selection, the MSC and GMSC are dominated by time-bin width alone, indicating reduced parameter coupling and greater robustness to processing variability. These findings demonstrate the feasibility and practical value of coherence-based spectral estimators for extracting α from reactor noise and support their broader application to multi-detector and irregular datasets in subcritical system characterization.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Advanced HVAC Humidity Control for Hot-Humid Climates

During this project we develop and validate a cost-effective, integrated control solution to improve humidity control and comfort for energy-efficient homes in hot-humid climates. This study focuses on developing a strategy that is effective, field tested, and practical for builders to install with minimal disruption to standard practices. A successful solution would simplify the transition to high-performance humidity control and be the basis for design and installation guidance. By relying on the central system as a starting point, the strategy employed minimizes system complexity and cost for builders, while improving comfort and operating cost for homeowners. The solution strategy was to coordinate the cooling, dehumidification, and ventilation functions of central, ducted HVAC systems to better control indoor humidity, improve occupant thermal comfort, and capture energy savings. The primary strategic goals were to: (1) optimize dehumidification by the central air-conditioning system, particularly during part-load conditions, using conventional equipment with modified control settings and lower system airflows; (2) maximize ventilation during heating/cooling on-cycles, to “bank” and condition outdoor ventilation air, and minimize ventilation during off-cycles; (3) quantify the effectiveness and energy impact of the dehumidification and ventilation strategies, while identifying a metric that would be useful to evaluate latent effectiveness. For the test houses in our study, located in Richmond Hill, Georgia; Houston, Texas; and Monroe, Louisiana we observed: (1) the indoor humidity did not exceed 60% RH during the monitored cooling season for 99% of the time in Richmond Hill, 96% of the time in Houston, and 90% of the time in Monroe; (2) the dehumidification strategy improved the steady-state latent capacity of the HVAC system at design conditions by 16% to 49% at the Houston test house and by 28% to 71% at the Monroe test house, depending on which mode the system was operating in; and (3) the good results at the test houses were primarily due to the amount of time the air-conditioning system operated in ramping or dehumidification modes, or both, particularly during the early cooling season. This study demonstrates that air conditioners or heat pumps with a single-stage compressor can provide good humidity control without the need for a two-stage or variable-stage compressor system. The airflow and control settings for ramping and dehumidification modes are critical to control indoor humidity in hot-humid climates, particularly during part-load and shoulder season conditions. The dehumidification strategy used in this study did not jeopardize the mechanical reliability of the cooling equipment. The strategies used in this study are applicable across various equipment brands, models, and efficiency levels, and also applicable to a broad range of homes in hot-humid climates. Results will vary by specific equipment, location, and house configuration and construction.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2023

Exploring the impacts of unprecedented climate extremes on forest ecosystems: hypotheses to guide modeling and experimental studies

Abstract. Climatic extreme events are expected to occur more frequently in the future, increasing the likelihood of unprecedented climate extremes (UCEs) or record-breaking events. UCEs, such as extreme heatwaves and droughts, substantially affect ecosystem stability and carbon cycling by increasing plant mortality and delaying ecosystem recovery. Quantitative knowledge of such effects is limited due to the paucity of experiments focusing on extreme climatic events beyond the range of historical experience. Here, we present a road map of how dynamic vegetation demographic models (VDMs) can be used to investigate hypotheses surrounding ecosystem responses to one type of UCE: unprecedented droughts. As a result of nonlinear ecosystem responses to UCEs that are qualitatively different from responses to milder extremes, we consider both biomass loss and recovery rates over time by reporting a time-integrated carbon loss as a result of UCE, relative to the absence of drought. Additionally, we explore how unprecedented droughts in combination with increasing atmospheric CO2 and/or temperature may affect ecosystem stability and carbon cycling. We explored these questions using simulations of pre-drought and post-drought conditions at well-studied forest sites using well-tested models (ED2 and LPJ-GUESS). The severity and patterns of biomass losses differed substantially between models. For example, biomass loss could be sensitive to either drought duration or drought intensity depending on the model approach. This is due to the models having different, but also plausible, representations of processes and interactions, highlighting the complicated variability of UCE impacts that still need to be narrowed down in models. Elevated atmospheric CO2 concentrations (eCO2) alone did not completely buffer the ecosystems from carbon losses during UCEs in the majority of our simulations. Our findings highlight the consequences of differences in process formulations and uncertainties in models, most notably related to availability in plant carbohydrate storage and the diversity of plant hydraulic schemes, in projecting potential ecosystem responses to UCEs. We provide a summary of the current state and role of many model processes that give way to different underlying hypotheses of plant responses to UCEs, reflecting knowledge gaps which in future studies could be tested with targeted field experiments and an iterative modeling–experimental conceptual framework.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

A Framework for Identifying Building Energy Models of Localized Utility Service Areas Using Smart Meter Data

Bottom-up load modeling of buildings offers a versatile approach to simulating baseline demand and scenarios of future technology evolution and adoption at the individual building level. This capability is essential to understanding how future load shapes may change with the adoption of electric equipment and vehicles, particularly as it relates to grid planning and infrastructure investments. Traditionally, grid planning techniques have used historical load data to predict future load and infrastructure needs. However, with the anticipated rise in adoption of electrification technologies such as heat pumps and electric vehicles, historical data become less reliable predictors of the future. By employing ResStock, a high-fidelity building stock modeling tool, we can fine-tune electrification scenarios and aggregate models to represent varying geographic resolutions of the grid system, while considering the underlying features of homes. This may enable a more accurate and responsive approach to anticipate and plan for the evolving landscape of energy demands. We present a new framework that leverages building stock energy modeling to identify building models that align with the load shapes and housing attributes of buildings with AMI data. This approach applies two model layers: (1) a classification step that identifies the presence of air conditioning, electric heating, and electric water heating, and (2) an optimization routine that identifies building energy models aligning with load profile data from advanced metering infrastructure meters. This report demonstrates one approach to deploying this framework, and presents results for three test cases that use both modeled and AMI data to assess performance. For a test case using AMI data in Fort Collins, Colorado, we observed a median monthly electricity load CV-RMSE of 16.6%, and a top ten daily heating and cooling median absolute percent error of 7.7% and 8.3%, respectively. For each AMI meter, we identify a set of potential energy models so that downstream use-cases can account for uncertainty driven by variability of baseline technologies and occupant behavior, which impact the response to electrification and energy efficiency scenarios. Our results indicate that ResStock has potential as a scalable solution for modeling residential energy demand at local grid resolutions. Its performance depends on location-specific factors, underlying building characteristics, and the level of aggregation, offering a path towards more precise and adaptive distribution grid planning for the evolving energy landscape.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Managing Solar Photovoltaic Integration in the Western United States: Resource Adequacy Considerations

This study examines the impact of reserve margin-based reliability assessment, as commonly used in capacity expansion models, on planning resource-adequate power systems under high penetrations of solar photovoltaics (PV). As a generation resource, PV is operationally different from the conventional dispatchable resources for which most capacity expansion models were designed. The question this study attempts to answer is whether large amounts of PV on a system (in this case, the Western Interconnection of North America) would bias the results of conventional reserve margin-based capacity expansion modeling towards an over- or under-provisioning of resource adequacy. This analysis used NREL’s Resource Planning Model (RPM) for capacity expansion modeling and NREL’s Probabilistic Resource Adequacy Suite (PRAS) for resource adequacy assessment. RPM uses a reserve margin requirement to enforce resource adequacy. PRAS, a collection of tools for studying the resource adequacy of power systems and the adequacy contributions of individual resources on a probabilistic basis, was used to compute multiple resource adequacy metrics across a number of simulated scenarios and system representations with differing regional detail. In all cases, including high PV penetrations (up to 33% annual generation from PV, interconnection-wide), RPM was able to produce resource-adequate systems as measured by normalized expected unserved energy and loss-of-load expectation results from PRAS. The accuracy of reserve margin approaches depends heavily on the underlying assumptions informing the capacity credit assigned to variable and energy-limited resources, particularly when such resources are abundant in the modeled system. RPM’s standard methodology for estimating variable and flexible resources’ capacity contributions, which is based on the top 100 hours of net load, did not appear to systematically undervalue or overvalue variable generation relative to a more rigorous equivalent firm capacity assessment using PRAS, although both over- and under-valuations were observed in specific scenarios. In the worst cases, the top 100 hour method underestimated the equivalent firm capacity of PV by two percentage points, and overestimated the equivalent firm capacity of PV by five percentage points. Calculating capacity contributions based on the top 10 hours of net load systematically underestimated equivalent firm capacities at more modest PV penetrations, but was often a better approximation of equivalent firm capacity than the existing 100-hour approach in scenarios with higher PV penetrations.

14 SOLAR ENERGY↗
DOE OSTI2021

Learning to Correct Climate Projection Biases

The fidelity of climate projections is often undermined by biases in climate models due to their simplification or misrepresentation of unresolved climate processes. While various bias correction methods have been developed to post-process model outputs to match observations, existing approaches usually focus on limited, low-order statistics, or break either the spatiotemporal consistency of the target variable, or its dependency upon model resolved dynamics. We develop a Regularized Adversarial Domain Adaptation (RADA) methodology to overcome these deficiencies, and enhance efficient identification and correction of climate model biases. Instead of pre-assuming the spatiotemporal characteristics of model biases, we apply discriminative neural networks to distinguish historical climate simulation samples and observation samples. The evidences based on which the discriminative neural networks make distinctions are applied to train the domain adaptation neural networks to bias correct climate simulations. We regularize the domain adaptation neural networks using cycle-consistent statistical and dynamical constraints. An application to daily precipitation projection over the contiguous United States shows that our methodology can correct all the considered moments of daily precipitation at approximately $1^\circ$ resolution, ensures spatiotemporal consistency and inter-field correlations, and can discriminate between different dynamical conditions. Our methodology offers a powerful tool for disentangling model parameterization biases from their interactions with the chaotic evolution of climate dynamics, opening a novel avenue toward big-data enhanced climate predictions.

58 GEOSCIENCES↗
DOE OSTI2020

Electricity Pricing aware Deep Reinforcement Learning based Intelligent HVAC Control

Recently, deep reinforcement learning (DRL) based intelligent control of Heating, Ventilation, and Air Conditioning (HVAC) has gained a lot of attention due to DRL's ability to optimally control HVAC for minimizing operational cost while maintaining resident's comfort. The success of such DRL-based techniques largely depends on the articulation of the problem in terms of states, actions, and reward function. Inclusion of the electricity pricing information in the problem formulation can play an important role in saving the cost of HVAC operation. However, less attention has been given in the literature on formulating well-crafted state features based on electricity pricing. In this work, we propose an approach for training the DRL model with a specific focus on feature engineering based on electricity pricing. During training, we generate random but sufficiently realistic electricity price signals so that the pre-trained DRL model is robust and adaptive to the dynamic and variable electricity prices. The validation results are encouraging and show the potential of ≈12%-15% savings in the one day cost of HVAC operation, proving the usefulness of including electricity pricing related features as state features.

Kurte, Kuldeep↗
DOE OSTI2022

Global soil profiles indicate depth-dependent soil carbon losses under a warmer climate

Soil organic carbon (SOC) changes under future climate warming are difficult to quantify in situ. Here we apply an innovative approach combining space-for-time substitution with meta-analysis to SOC measurements in 113,013 soil profiles across the globe to estimate the effect of future climate warming on steady-state SOC stocks. We find that SOC stock will reduce by 6.0 ± 1.6% (mean±95% confidence interval), 4.8 ± 2.3% and 1.3 ± 4.0% at 0–0.3, 0.3–1 and 1–2 m soil depths, respectively, under 1 °C air warming, with additional 4.2%, 2.2% and 1.4% losses per every additional 1 °C warming, respectively. The largest proportional SOC losses occur in boreal forests. Existing SOC level is the predominant determinant of the spatial variability of SOC changes with higher percentage losses in SOC-rich soils. Our work demonstrates that warming induces more proportional SOC losses in topsoil than in subsoil, particularly from high-latitudinal SOC-rich systems.

54 ENVIRONMENTAL SCIENCES↗