Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “climate observations”

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 253 records · Page 14

Toward data assimilation of ship-induced aerosol–cloud interactions

Satellite imagery can detect temporary cloud trails or ship tracks formed from aerosols emitted from large ships traversing our oceans, a phenomenon that global climate models cannot directly reproduce. Ship tracks are observable examples of marine cloud brightening, a potential solar climate intervention that shows promise in helping combat climate change. In this paper, we demonstrate a simulation-based approach in learning the behavior of ship tracks based upon a novel stochastic emulation mechanism. Our method uses wind fields to determine the movement of aerosol–cloud tracks and uses a stochastic partial differential equation (SPDE) to model their persistence behavior. This SPDE incorporates both a drift and diffusion term which describes the movement of aerosol particles via wind and their diffusivity through the atmosphere, respectively. We first present our proposed approach with examples using simulated wind fields and ship paths. We then successfully demonstrate our tool by applying the approximate Bayesian computation method-sequential Monte Carlo for data assimilation.

54 ENVIRONMENTAL SCIENCES↗

Mesoscale Cellular Convection Detection and Classification Using Convolutional Neural Networks: Insights From Long-Term Observations at ARM Eastern North Atlantic Site

Marine boundary layer clouds are crucial in Earth's climate system. They frequently manifest as closed or open cell mesoscale cellular convection (MCC). MCC clouds are challenging to represent accurately in current climate models, highlighting the need for detailed observational data sets and in-depth analyses. This study utilizes over 8 years of observations from the U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) User Facility Eastern North Atlantic (ENA) site at Graciosa Island, Azores, to investigate these clouds. We first apply a convolutional neural network with a U-Net architecture to classify open and closed cells, marking the first application of such an approach for automatically detecting MCC patterns from ground-based radar measurements. This method addresses some observational gaps in satellite data related to low temporal resolution, nighttime challenges, and limited vertical structure capture. The analysis of the MCC cases shows clear differences between closed and open MCCs: Closed MCC clouds are characterized by lower cloud tops and bases, shallower cloud geometrical depth, weaker horizontal wind speeds, stronger atmospheric stability, and a more homogeneous liquid water path than open MCCs. Finally, we demonstrate two potential applications of our radar-based MCC classifications: (a) facilitating the investigation of aerosol-cloud interactions and (b) exploring meteorological factors along with MCC's evolution by integrating satellite imagery and back-trajectory analysis. The identified MCC cases offer a valuable resource for the scientific community to study MCC processes further and improve climate model accuracy.

54 ENVIRONMENTAL SCIENCES↗

Warming homogenizes apparent temperature sensitivity of ecosystem respiration

Warming-induced carbon loss through terrestrial ecosystem respiration (Re) is likely getting stronger in high latitudes and cold regions because of the more rapid warming and higher temperature sensitivity of Re (Q 10 ). However, it is not known whether the spatial relationship between Q 10 and temperature also holds temporally under a future warmer climate. Here, we analyzed apparent Q 10 values derived from multiyear observations at 74 FLUXNET sites spanning diverse climates and biomes. We found warming-induced decline in Q 10 is stronger at colder regions than other locations, which is consistent with a meta-analysis of 54 field warming experiments across the globe. We predict future warming will shrink the global variability of Q 10 values to an average of 1.44 across the globe under a high emission trajectory (RCP 8.5) by the end of the century. Therefore, warming-induced carbon loss may be less than previously assumed because of Q 10 homogenization in a warming world.

54 ENVIRONMENTAL SCIENCES↗

Snow thermal conductivity and conductive flux in the Central Arctic: Estimates from observations and implications for models

During the Arctic winter, the conductive heat flux through the sea ice and snow balances the radiative and turbulent heat fluxes at the surface. Snow on sea ice is a thermal insulator that reduces the magnitude of the conductive flux. The thermal conductivity of snow, that is, how readily energy is conducted, is known to vary significantly in time and space from observations, but most forecast and climate models use a constant value. This work begins with a demonstration of the importance of snow thermal conductivity in a regional coupled forecast model. Varying snow thermal conductivity impacts the magnitudes of all surface fluxes, not just conduction, and their responses to atmospheric forcing. Given the importance of snow thermal conductivity in models, we use observations from sea ice mass balance buoys installed during the Multidisciplinary drifting Observatory for the Study of Arctic Climate expedition to derive the profiles of thermal conductivity, density, and conductive flux. From 13 sites, median snow thermal conductivity ranges from 0.33 W m −1 K −1 to 0.47 W m −1 K −1 with a median from all data of 0.39 W m −1 K −1 from October to February. In terms of surface energy budget closure, estimated conductive fluxes are generally smaller than the net atmospheric flux by as much as 20 W m −2 , but the average residual during winter is −6 W m −2 , which is within the uncertainties. The spatial variability of conductive heat flux is highest during clear and cold time periods. Higher surface temperature, which often occurs during cloudy conditions, and thicker snowpacks reduce temporal and spatial variability. These relationships are compared between observations and the coupled forecast model, emphasizing both the importance and challenge of describing thermodynamic parameters of snow cover for modeling the Arctic as a coupled system.

54 ENVIRONMENTAL SCIENCES↗

Estimates of Southern Hemispheric Gravity Wave Momentum Fluxes across Observations, Reanalyses, and Kilometer-Scale Numerical Weather Prediction Model

Abstract Gravity waves (GWs) are among the key drivers of the meridional overturning circulation in the mesosphere and upper stratosphere. Their representation in climate models suffers from insufficient resolution and limited observational constraints on their parameterizations. This obscures assessments of middle atmospheric circulation changes in a changing climate. This study presents a comprehensive analysis of stratospheric GW activity above and downstream of the Andes from 1 to 15 August 2019, with special focus on GW representation ranging from an unprecedented kilometer-scale global forecast model (1.4 km ECMWF IFS), ground-based Rayleigh lidar (CORAL) observations, modern reanalysis (ERA5), to a coarse-resolution climate model (EMAC). Resolved vertical flux of zonal GW momentum (GWMF) is found to be stronger by a factor of at least 2–2.5 in IFS compared to ERA5. Compared to resolved GWMF in IFS, parameterizations in ERA5 and EMAC continue to inaccurately generate excessive GWMF poleward of 60°S, yielding prominent differences between resolved and parameterized GWMFs. A like-to-like validation of GW profiles in IFS and ERA5 reveals similar wave structures. Still, even at ∼1 km resolution, the resolved waves in IFS are weaker than those observed by lidar. Further, GWMF estimates across datasets reveal that temperature-based proxies, based on midfrequency approximations for linear GWs, overestimate GWMF due to simplifications and uncertainties in GW wavelength estimation from data. Overall, the analysis provides GWMF benchmarks for parameterization validation and calls for three-dimensional GW parameterizations, better upper-boundary treatment, and vertical resolution increases commensurate with increases in horizontal resolution in models, for a more realistic GW analysis. Significance Statement Gravity wave–induced momentum forcing forms a key component of the middle atmospheric circulation. However, complete knowledge of gravity waves, their atmospheric effects, and their long-term trends are obscured due to limited global observations, and the inability of current climate models to fully resolve them. This study combines a kilometer-scale forecast model, modern reanalysis, and a coarse-resolution climate model to first compare the resolved and parameterized momentum fluxes by gravity waves generated over the Andes, and then evaluate the fluxes using a state-of-the-art ground-based Rayleigh lidar. Our analysis reveals shortcomings in current model parameterizations of gravity waves in the middle atmosphere and highlights the sensitivity of the estimated flux to the formulation used.

Meteorology & Atmospheric Sciences↗

Rapid attribution analysis of the extraordinary heat wave on the Pacific coast of the US and Canada in June 2021

Towards the end of June 2021, temperature records were broken by several degrees Celsius in several cities in the Pacific Northwest areas of the US and Canada, leading to spikes in sudden deaths and sharp increases in emergency calls and hospital visits for heat-related illnesses. Here we present a multi-model, multi-method attribution analysis to investigate the extent to which human-induced climate change has influenced the probability and intensity of extreme heat waves in this region. Based on observations, modelling and a classical statistical approach, the occurrence of a heat wave defined as the maximum daily temperature (TXx) observed in the area 45–52 °N, 119–123 °W, was found to be virtually impossible without human-caused climate change. The observed temperatures were so extreme that they lay far outside the range of historical temperature observations. This makes it hard to state with confidence how rare the event was. Using a statistical analysis that assumes that the heat wave is part of the same distribution as previous heat waves in this region led to a first-order estimation of the event frequency of the order of once in 1000 years under current climate conditions. Using this assumption and combining the results from the analysis of climate models and weather observations, we found that such a heat wave event would be at least 150 times less common without human-induced climate change. Also, this heat wave was about 2 °C hotter than a 1-in-1000-year heat wave would have been in 1850–1900, when global mean temperatures were 1.2 °C cooler than today. Looking into the future, in a world with 2 °C of global warming (0.8 °C warmer than today), a 1000-year event would be another degree hotter. Our results provide a strong warning: our rapidly warming climate is bringing us into uncharted territory with significant consequences for health, well-being and livelihoods. Adaptation and mitigation are urgently needed to prepare societies for a very different future.

54 ENVIRONMENTAL SCIENCES↗

Multiscale Wildfire Simulation Framework and Remote Sensing

Wildfire as one type of climate extreme events causes huge socioeconomic losses and damages. Large wildfires (i.e., generated pyrocumulonimbus (PyroCb)) can inject tremendous amounts of smoke into the stratosphere, where black carbon and organic carbon aerosols can persist months to years and influence climate by imposing a significant reduction in the radiative forcing like that associated with large volcanic eruptions or proposed via climate interventions such as geoengineering. Both observations and numerical modeling results clearly indicate an increasing trend in wildfire frequency and intensity in many regions during the recent decades with climate change. However, current understanding of wildfire remains largely uncertain owing to limitations of modeling capabilities in representing the multiscale wildfire physics and dynamics and a scarcity of observations constraining important wildfire and environmental variables. This study primarily aims to improve the wildfire simulation capabilities in the state-of-the-art climate model by filling in two major gaps: (1) model resolution is typically too coarse to resolve fine scale processes associated with fires, and (2) chemistry and aerosol processes in fire smoke are poorly represented.

54 ENVIRONMENTAL SCIENCES↗

Pan-Pacific low-frequency modes of sea level and climate variability

Tide gauges provide a long observational record that can inform the nature of satellite-era basin-scale sea level trends. However, common signals must be extracted from geographically sparse records. Here, by applying low-frequency component analysis (LFCA) to tide gauge records and surface climate reconstructions, we isolate three coherent modes of Pacific Ocean variability that we ascribe to: a secular, greenhouse gas–driven climate change (LFC1); a nonlinear mode of variability with a reversal around 1980, potentially linked to aerosols (LFC2); and the Interdecadal Pacific Oscillation (LFC3). Although sea level trend patterns reflect the superimposed contribution of all modes, satellite-era trends are dominated by an increasing phase of LFC2: They are thus potentially unrepresentative of both longer-term historical patterns and those expected in the future.

Science & Technology - Other Topics↗

Development of the ARM Lagrangian Large-Scale Forcing Data (ARMLAGTRAJ) Value-Added Product Based on the lagtraj Framework

The Atmospheric Radiation Measurement (ARM) large-scale forcing data developed based on the constrained variational analysis (VARANAL) value-added product (VAP) (Zhang and Lin 1997, Zhang et al. 2001, Xie et al. 2004, Tang et al. 2019) has been widely used for single-column models (SCMs), cloud-resolving models (CRMs), and large-eddy simulation models (LESs) to understand and improve physical processes in models. Recently, the U.S. Department of Energy (DOE) ARM user facility conducted several major field campaigns using ship-based moving observational platforms. For example, the Marine ARM GPCI Investigation of Clouds (MAGIC) field campaign focused on the role of subtropical marine-boundary layer (MBL) clouds, and the Multidisciplinary Drifting Observatory for the Study of Arctic Climate (MOSAiC) field campaign aimed to improve understanding of the coupled climate systems in the Arctic. Observations from moving platforms are critical to provide a comprehensive characterization of coupled-system processes associated with all stages of the cloud and/or sea-ice life cycle. Traditional ARM large-scale forcing data have been developed at fixed locations. They need to be extended to include these moving platforms to address data needs for ship-based field campaigns or to support LES modeling in a Lagrangian framework. With these considerations in mind, we develop ARM-type Lagrangian large-scale forcing data sets based on the lagtraj framework (Boeing et al. 2020) with notable enhancements in generating forcings that are more suitable for ARM field campaigns. The lagtraj is a novel tool that generates forcings for LES and SCM simulation in both Lagrangian and Eulerian perspective. This technical report focuses on the major changes we performed on the lagtraj algorithm and provides an overview of the ARM Lagrangian Large-Scale Forcing Data (ARMLAGTRAJ) value-added products.

54 ENVIRONMENTAL SCIENCES↗

Network of networks: Time series clustering of AmeriFlux sites

Environmental observation networks, such as AmeriFlux, are foundational for monitoring ecosystem response to climate change, management practices, and natural disturbances; however, their effectiveness depends on their representativeness for the regions or continents. We proposed an empirical, time series approach to quantify the similarity of ecosystem fluxes across AmeriFlux sites. We extracted the diel and seasonal characteristics (i.e., amplitudes, phases) from carbon dioxide, water vapor, energy, and momentum fluxes, which reflect the effects of climate, plant phenology, and ecophysiology on the observations, and explored the potential aggregations of AmeriFlux sites through hierarchical clustering. While net radiation and temperature showed latitudinal clustering as expected, flux variables revealed a more uneven clustering with many small (number of sites < 5), unique groups and a few large (> 100) to intermediate (15–70) groups, highlighting the significant ecological regulations of ecosystem fluxes. Many identified unique groups were from under-sampled ecoregions and biome types of the International Geosphere-Biosphere Programme (IGBP), with distinct flux dynamics compared to the rest of the network. At the finer spatial scale, local topography, disturbance, management, edaphic, and hydrological regimes further enlarge the difference in flux dynamics within the groups. Nonetheless, our clustering approach is a data-driven method to interpret the AmeriFlux network, informing future cross-site syntheses, upscaling, and model-data benchmarking research. Finally, we highlighted the unique and underrepresented sites in the AmeriFlux network, which were found mainly in Hawaii and Latin America, mountains, and at under-sampled IGBP types (e.g., urban, open water), motivating the incorporation of new/unregistered sites from these groups.

54 ENVIRONMENTAL SCIENCES↗

A Methodology for Estimating the Energy and Moisture Budget of the Convective Boundary Layer Using Continuous Ground-Based Infrared Spectrometer Observations

Abstract Land–atmosphere interactions play a critical role in both the atmospheric water and energy cycles. Changes in soil moisture and vegetation alter the partitioning of surface water and energy fluxes, influencing diurnal evolution of the planetary boundary layer (PBL). The mixing-diagram framework has proven useful in understanding the evolution of the heat and moisture budget within the convective boundary layer (CBL). We demonstrate that observations from the Department of Energy Atmospheric Radiation Measurement (ARM) Southern Great Plains site provide all of the needed inputs needed for the mixing-diagram framework, allowing us to quantify the impact from the surface fluxes, advection, radiative heating, encroachment, and entrainment on the evolution of the CBL. Profiles of temperature and humidity retrieved from the ground-based infrared spectrometer [Atmospheric Emitted Radiance Interferometer (AERI)] are a critical component in this analysis. Large-eddy simulation results demonstrate that mean mixed-layer values derived are shown to be critical to close the energy and moisture budgets. A novel approach demonstrated here is the use of network of AERIs and Doppler lidars to quantify the advective fluxes of heat and moisture. The framework enables the estimation of the entrainment fluxes as a residual, providing a way to observe the entrainment fluxes without using multiple lidar systems. The high temporal resolution of the AERI observations enables the morning, midday, and afternoon evolution of the CBL to be quantified. This work provides a new way to use observations in this framework to evaluate weather and climate models. Significance Statement The energy and moisture budget of the planetary boundary layer (PBL) is influenced by multiple sources, and accurately representing this evolution in numerical models is critical for weather forecasts and climate predictions. The mixing-diagram approach, driven by profiling observations as illustrated here, provides a powerful way to quantify the contributions from each of these sources. In particular, the energy and moisture mixed into the PBL from above the PBL can be determined accurately from ground-based remote sensors using this approach.

54 ENVIRONMENTAL SCIENCES↗

Biome-scale temperature sensitivity of ecosystem respiration revealed by atmospheric CO2 observations

Abstract The temperature sensitivity of ecosystem respiration regulates how the terrestrial carbon sink responds to a warming climate but has been difficult to constrain observationally beyond the plot scale. Here we use observations of atmospheric CO 2 concentrations from a network of towers together with carbon flux estimates from state-of-the-art terrestrial biosphere models to characterize the temperature sensitivity of ecosystem respiration, as represented by the Arrhenius activation energy, over various North American biomes. We infer activation energies of 0.43 eV for North America and 0.38 eV to 0.53 eV for major biomes therein, which are substantially below those reported for plot-scale studies (approximately 0.65 eV). This discrepancy suggests that sparse plot-scale observations do not capture the spatial-scale dependence and biome specificity of the temperature sensitivity. We further show that adjusting the apparent temperature sensitivity in model estimates markedly improves their ability to represent observed atmospheric CO 2 variability. This study provides observationally constrained estimates of the temperature sensitivity of ecosystem respiration directly at the biome scale and reveals that temperature sensitivities at this scale are lower than those based on earlier plot-scale studies. These findings call for additional work to assess the resilience of large-scale carbon sinks to warming.

54 ENVIRONMENTAL SCIENCES↗

CMIP5 Climate Models Overestimate Cooling by Volcanic Aerosols

We compare projections of the observed hemispherical mean surface temperature (HadCRUT4.6.0.0) and the ensemble mean of CMIP5 climate models' simulations on a set of standard regression model forcing variables. We find that the volcanic aerosol regression coefficients of the CMIP5 simulations are consistently significantly larger (by 40–49%) than the volcanic aerosol coefficients of the observed temperature. The probability that the observed differences are caused just by chance is much less than 0.01. The overestimate is due to the climate models' response to volcanic aerosol radiative forcing. The largest overestimate occurs in the winter season of each hemisphere. We hypothesize that the models' parameterization of aerosol-cloud interactions within ice and mixed phase clouds is a likely source of this discrepancy. Furthermore, the models significantly underestimate the effect of solar variability on temperature for both hemispheres.

54 ENVIRONMENTAL SCIENCES↗

Impacts of long-range transport of aerosols on marine-boundary-layer clouds in the eastern North Atlantic

Vertical profiles of aerosols are inadequately observed and poorly represented in climate models, contributing to the current large uncertainty associated with aerosol–cloud interactions. The US Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) Aerosol and Cloud Experiments in the Eastern North Atlantic (ACE-ENA) aircraft field campaign near the Azores islands provided ample observations of vertical distributions of aerosol and cloud properties. Here we utilize the in situ aircraft measurements from the ACE-ENA and ground-based remote-sensing data along with an aerosol-aware Weather Research and Forecast (WRF) model to characterize the aerosols due to long-range transport over a remote region and to assess their possible influence on marine-boundary-layer (MBL) clouds. The vertical profiles of aerosol and cloud properties measured via aircraft during the ACE-ENA campaign provide detailed information revealing the physical contact between transported aerosols and MBL clouds. The European Centre for Medium-Range Weather Forecasts Copernicus Atmosphere Monitoring Service (ECMWF-CAMS) aerosol reanalysis data can reproduce the key features of aerosol vertical profiles in the remote region. The cloud-resolving WRF sensitivity experiments with distinctive aerosol profiles suggest that the transported aerosols and MBL cloud interactions (ACIs) require not only aerosol plumes to get close to the marine-boundary-layer top but also large cloud top height variations. Based on those criteria, the observations show that the occurrence of ACIs involving the transport of aerosol over the eastern North Atlantic (ENA) is about 62 % in summer. For the case with noticeable long-range-transport aerosol effects on MBL clouds, the susceptibilities of droplet effective radius and liquid water content are -0.11 and +0.14, respectively. When varying by a similar magnitude, aerosols originating from the boundary layer exert larger microphysical influence on MBL clouds than those entrained from the free troposphere.

54 ENVIRONMENTAL SCIENCES↗

Role of the Terrestrial Biosphere in Atmospheric Chemistry and Climate

The terrestrial biosphere–atmosphere interface provides a key chemical, biological, and physical lower boundary for the atmosphere. The presence of vegetation itself modifies the physical boundary, or the biogeophysical aspects of the system, by controlling important climate drivers such as soil moisture, light environment, and temperature. The leaf surface area of the terrestrial biosphere provides additional surface area for emissions, and it can be up to 55% of the total Earth’s surface area during the boreal summer. Vegetation also influences the biogeochemical aspects of the system by emitting a broad suite of reactive trace gases such as biogenic volatile organic compound (BVOC) emissions and climate-relevant primary biological aerosol particles (PBAP). Many of these emissions are a function of meteorological and climatological conditions at the surface, including temperature, light environment, soil moisture, and winds. Once emitted, they can be processed in the troposphere through a suite of chemical reactions. BVOC can contribute to the formation of ozone and secondary organic aerosols (SOA), and PBAP can rupture to form smaller particles with climatic relevance. These emissions and subsequent aerosol products can influence atmospheric processes that affect the surface climate, such as the attenuation of radiation, the formation of greenhouse gases such as ozone that can feedback to surface air temperature, and the alteration of clouds and subsequent precipitation. These atmospheric changes can then feedback to the land surface and emissions themselves, creating positive or negative feedback loops that can dampen or amplify the emission response. For the dominant BVOC isoprene, the feedback response to temperature can be positive or negative depending on ambient temperatures that drive isoprene emissions. The feedback response to soil moisture and precipitation can be positive, negative, or uncoupled depending on the moisture content of the soil and the total atmospheric aerosol loading. For light, the isoprene response can be positive or negative depending on the role of diffuse light. Overall, these feedbacks highlight the dynamical response of the biosphere to changing atmospheric conditions across a range of time scales, from minutes for trace gases and aerosols, to months for phenological changes, to years for land cover and land use change. Here, the dynamic aspect of this system requires us to understand, simulate, and predict the complex feedbacks between the biosphere and atmosphere and understand their role in the simulation and understanding of climate and global change. From the observational perspective, these feedbacks are challenging to identify in observations, and predictive modeling tools provide a crucial link for understanding how these feedbacks will change under warming climate scenarios.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Recent Weakening of the Global Radiative Feedback

Earth's climate stability, characterized by the global radiative feedback parameter (𝜆), varies decadally due to changing surface temperature patterns. Recent variations in 𝜆 are poorly understood as coordinated model simulations typically end in 2014. We apply a convolutional neural network trained on climate model simulations to observation-based surface temperature reconstructions to estimate variations in 𝜆 up to 2025. We find that 𝜆 reached a minimum (maximum stability) around the mid 1990s (𝜆 ≃ −3 Wm −2 /K), but has since weakened significantly (𝜆 ≃ −2 Wm −2 /K). We confirm these results with climate model simulations extended to 2022. The recent 𝜆 weakening is not significantly affected by El Niño Southern Oscillation or Pacific Decadal Oscillation. Attribution reveals that warming in the subtropical Northeast Pacific is an important driver of the recently weakened feedback, confirmed by targeted experiments in E3SMv2. Our approach enables near real-time monitoring of Earth's climate stability.

54 ENVIRONMENTAL SCIENCES↗

Hemispheric Asymmetry of Phase Partition in Mixed‐Phase Clouds Based on Near Global‐Scale Airborne Observations

Mixed-phase clouds contribute to substantial uncertainties in global climate models due to their complex microphysical properties. Former model evaluations almost exclusively rely on satellite observations to assess cloud phase distributions globally. This study investigated mixed-phase cloud properties using near global-scale in situ observation data sets from 14 flight campaigns in combination with collocated output from a global climate model. The Southern Hemisphere (SH) shows significantly higher occurrence frequencies and higher mass fractions of supercooled liquid water than Northern Hemisphere (NH) based on observations at 0.2 and 100 km horizontal scales. Such hemispheric asymmetry is not captured by the model. The model also consistently overestimates liquid water content (LWC) in all cloud phases but shows ice water content (IWC) biases that vary with phase. Key processes contributing to model biases in phase partition can be identified through the combination of evaluation of phase frequency, liquid mass fraction, LWC and IWC.

Yang, Ching An [San Jose State University, CA (Uni↗

Emergent temperature sensitivity of soil organic carbon driven by mineral associations

Abstract Soil organic matter decomposition and its interactions with climate depend on whether the organic matter is associated with soil minerals. However, data limitations have hindered global-scale analyses of mineral-associated and particulate soil organic carbon pools and their benchmarking in Earth system models used to estimate carbon cycle–climate feedbacks. Here we analyse observationally derived global estimates of soil carbon pools to quantify their relative proportions and compute their climatological temperature sensitivities as the decline in carbon with increasing temperature. We find that the climatological temperature sensitivity of particulate carbon is on average 28% higher than that of mineral-associated carbon, and up to 53% higher in cool climates. Moreover, the distribution of carbon between these underlying soil carbon pools drives the emergent climatological temperature sensitivity of bulk soil carbon stocks. However, global models vary widely in their predictions of soil carbon pool distributions. We show that the global proportion of model pools that are conceptually similar to mineral-protected carbon ranges from 16 to 85% across Earth system models from the Coupled Model Intercomparison Project Phase 6 and offline land models, with implications for bulk soil carbon ages and ecosystem responsiveness. To improve projections of carbon cycle–climate feedbacks, it is imperative to assess underlying soil carbon pools to accurately predict the distribution and vulnerability of soil carbon.

54 ENVIRONMENTAL SCIENCES↗