Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “aerosol_optical_depth”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Biomass burning emission analysis based on MODIS aerosol optical depth and AeroCom multi-model simulations: implications for model constraints and emission inventories

We assessed the biomass burning (BB) smoke aerosol optical depth (AOD) simulations of 11 global models that participated in the AeroCom phase III BB emission experiment. By comparing multi-model simulations and satellite observations in the vicinity of fires over 13 regions globally, we (1) assess model-simulated BB AOD performance as an indication of smoke source–strength, (2) identify regions where the common emission dataset used by the models might underestimate or overestimate smoke sources, and (3) assess model diversity and identify underlying causes as much as possible. Using satellite-derived AOD snapshots to constrain source strength works best where BB smoke from active sources dominates background non-BB aerosol, such as in boreal forest regions and over South America and southern hemispheric Africa. The comparison is inconclusive where the total AOD is low, as in many agricultural burning areas, and where the background is high, such as parts of India and China. Many inter-model BB AOD differences can be traced to differences in values for the mass ratio of organic aerosol to organic carbon, the BB aerosol mass extinction efficiency, and the aerosol loss rate from each model. The results point to a need for increased numbers of available BB cases for study in some regions and especially to a need for more extensive regional-to-global-scale measurements of aerosol loss rates and of detailed particle microphysical and optical properties; this would both better constrain models and help distinguish BB from other aerosol types in satellite retrievals. More generally, there is the need for additional efforts at constraining aerosol source strength and other model attributes with multi-platform observations.

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

Analysis of contrasting aerosol indirect effects in liquid water clouds over the northern part of Arabian Sea

The extensive daily statistics of aerosol properties, cloud properties, and their mutual correlations provide crucial information for better assessing future climate change. Here, in this paper, 14 years (2010–2023) of data from the Moderate Resolution Imaging Spectroradiometer (MODIS) are analyzed over the northern part of Arabian Sea (Latitude: 21°–25° N and Longitude: 62°–68° E) to assess the characteristics of aerosols and clouds and their relationships under different meteorological conditions. When aerosol optical depth (AOD) is less than ~0.7, the observations exhibit a positive correlation between AOD and cloud droplet effective radius (CDR) but negative correlations between AOD and cloud droplet number concentration (CDNC), between AOD and cloud optical depth (COD), between AOD and cloud liquid water path (CLWP), and between AOD and cloud geometrical thickness (H). The corresponding aerosol-cloud correlations change signs when the AOD values are larger than 0.7. However, the single folded positive AOD-cloud fraction (CF) relationship is observed in both AOD regimes. Similar correlations are also observed between precipitable water vapor (PWV) and CDR, CDNC, COD, H, CF and CLWP, together with a positive correlation between PWV and AOD. Further isolation of the environmental effects from aerosol effects by stratifying AOD and cloud data into different LTS and PWV bins shows that the signature of the well-known Twomey effect is observed under high LTS-high PWV conditions, while an opposite effect (anti-Twomey) is observed under low PWV conditions, regardless of LTS values. Additionally, negative correlations between AOD and COD, AOD and CLWP, and AOD and H are observed under low LTS, regardless of PWV conditions, with a slight positive correlation when AOD >0.4 under high LTS and PWV conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

MODIS Aerosol and Low‐Cloud Retrievals: Orographic Effects in the Wake of Macquarie Island

Abstract Visible/infrared imagery from passive satellites is commonly relied upon to study low cloud microphysics over oceanic regions, including for the Southern Ocean (SO), but relatively little validation has been undertaken for the SO. In this article, we compare low‐cloud effective radius ( r e ), cloud droplet number concentration ( N d ) and cloud liquid water path (LWP) retrievals from the NASA Moderate Imaging Spectroradiometer (MODIS) with surface measurements collected during the Macquarie Island Cloud and Radiation Experiment (MICRE). MODIS 3.7‐μm band retrievals show little bias and moderately good correlation relative to MICRE retrievals for liquid‐phase low clouds when restricted to Solar Zenith Angles <65° on spatial scales of 50–100 km. However, the low overall bias in 3.7‐μm band effective radius ( r e 3.7 ) retrievals partly results from cancellation of errors: r e 3.7 is overestimated for non‐to‐lightly precipitating clouds, and underestimated for heavier drizzling clouds by ∼1–1.5 μm. 1.6‐μm and 2.1‐μm band r e retrievals are biased high. N d may likewise be slightly under‐ or overestimated depending on the concentration, but there is insufficient data to provide confidence in this result. Interestingly, a composite of MODIS retrievals from 2002 to 2020 show a distinct region of enhanced cloud cover and N d (and lower r e ) in the wake of Macquarie Island associated with orographic cloud formation. MODIS retrievals of aerosol optical depth (AOD) and Angstrom Exponent (AE) upwind and downwind of the island do not differ significantly. Comparison with MICRE measurements suggests that MODIS Collection 6 AOD retrievals are reasonable, while AE is problematically large.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Tracing the impacts of Mount Pinatubo eruption on regional climate using spatially-varying changepoint detection

Significant events, such as volcanic eruptions, can have global and long-lasting impacts on climate. These global impacts, however, are not uniform across space and time. Understanding how the Mt. Pinatubo eruption affects global and regional climate is of great interest for predicting the impact on climate due to similar events as well as understanding the possible effect of the stratospheric aerosol injections proposed to combat climate change. While many studies illustrated the impact of the Pinatubo eruption on a global scale, studies at a fine regional scale are scarce. Here, we propose a Bayesian spatially-varying changepoint detection and estimation method to trace the impact of Mt. Pinatubo eruption on regional climate. Our approach takes into account the diffusing nature and spatial correlation of the climate changes attributed to the volcanic eruption. We illustrate our method and demonstrate its advantages over an existing changepoint detection method through simulations. Finally, we apply our method to monthly stratospheric aerosol optical depth and surface temperature data from 1985 to 1995 to detect and estimate changepoints following the 1991 Mt. Pinatubo eruption. Our results quantitatively characterize the spatial pattern of the eruption’s impact on regional climate, complementing the previous studies on the global impact of the Pinatubo eruption.

Aerosol optical depth↗
DOE OSTI2025

Sensitivity of Regional WRF‐Chem Air Quality and Weather Simulations to Biomass‐Burning Emission Data Sets: A Case Study of the Impact of Canadian Wildfire on the US°

This study focuses on the period from June 26 to 29, 2023, when record‐breaking Canadian wildfires severely impacted air quality in the Midwest United States. Using the Weather Research and Forecasting Model with Chemistry (WRF‐Chem) and four biomass‐burning data sets (Fire Inventory from NCAR version 1, Fire Inventory from NCAR version 2.5, Quick Fire Emissions Data set [QFED], and Regional ABI‐VIIRS Emission), we analyzed aerosol transport from Canada to the US and assessed the model's accuracy in predicting PM 2.5 , O 3 , CO and aerosol weather feedback. Model simulations were compared with ground‐based and remote sensing observations as well as field measurements from the Community Research on Climate and Urban Science (CROCUS) project. Our findings show that the movement of a low‐pressure system from the Great Lakes to the Atlantic, combined with the high‐pressure system over the Atlantic, caused the transport of aerosols from Canadian wildfires to the US. Results show WRF‐Chem significantly underestimated key atmospheric components: aerosol optical depth (AOD) by over 50%, PM 2.5 by 65%–90% and peak O 3 concentrations by 50%–55% across four biomass burning data sets. Additionally, CO and NO 2 concentrations were underpredicted. The substantial underestimation of PM 2.5 led to an overestimation of temperature by up to 3.6 °C primarily due to excessive downward shortwave radiation, which resulted from the underestimation of direct aerosol effects and an increase in sensible heat flux. Among the biomass‐burning data sets, QFED produced the most accurate AOD and PM 2.5 predictions due to improved wildfire emission estimates, leading to a 1.0 to 1.5 °C reduction in temperature overestimation during the daytime. These findings underscore the need for improving wildfire emission estimates for trace gases and aerosols to enhance air quality and weather feedback predictions.

WRF-chem model↗
DOE OSTI2024

Enhanced light absorption for solid-state brown carbon from wildfires due to organic and water coatings

Abstract Wildfires emit solid-state strongly absorptive brown carbon (solid S-BrC, commonly known as tar ball), critical to Earth’s radiation budget and climate, but their highly variable light absorption properties are typically not accounted for in climate models. Here, we show that from a Pacific Northwest wildfire, over 90% of particles are solid S-BrC with a mean refractive index of 1.49 + 0.056 i at 550 nm. Model sensitivity studies show refractive index variation can cause a ~200% difference in regional absorption aerosol optical depth. We show that ~50% of solid S-BrC particles from this sample uptake water above 97% relative humidity. We hypothesize these results from a hygroscopic organic coating, potentially facilitating solid S-BrC as nuclei for cloud droplets. This water uptake doubles absorption at 550 nm and the organic coating on solid S-BrC can lead to even higher absorption enhancements than water. Incorporating solid S-BrC and water interactions should improve Earth’s radiation budget predictions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Radiative impact of record-breaking wildfires from integrated ground-based data

The radiative effects of wildfires have been traditionally estimated by models using radiative transfer calculations. Assessment of model-predicted radiative effects commonly involves information on observation-based aerosol optical properties. However, lack or incompleteness of this information for dense plumes generated by intense wildfires reduces substantially the applicability of this assessment. Here we introduce a novel method that provides additional observational constraints for such assessments using widely available ground-based measurements of shortwave and spectrally resolved irradiances and aerosol optical depth (AOD) in the visible and near-infrared spectral ranges. We apply our method to quantify the radiative impact of the record-breaking wildfires that occurred in the Western US in September 2020. For our quantification we use integrated ground-based data collected at the Atmospheric Measurements Laboratory in Richland, Washington, USA with a location frequently downwind of wildfires in the Western US. We demonstrate that remarkably dense plumes generated by these wildfires strongly reduced the solar surface irradiance (up to 70% or 450 Wm -2 for total shortwave flux) and almost completely masked the sun from view due to extremely large AOD (above 10 at 500 nm wavelength). We also demonstrate that the plume-induced radiative impact is comparable in magnitude with those produced by a violent volcano eruption occurred in the Western US in 1980 and continental cumuli.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Two decades of aerosol trends over India: seasonal characteristics and urban-rural dynamics

Abstract India faces significant air quality challenges, with one of the highest air pollution levels of any country in the world. Here, we examine two decades (2001–2019) of both particulate matter (PM 2.5 ) concentration and aerosol optical depth (AOD) over the country. Increases are seen between the two decadal averages, for 2001–2010 and 2011–2019, in western India, particularly in the Indo-Gangetic Plain (IGP). IGP region, including Bihar, West Bengal, Jharkhand, and Uttar Pradesh, shows the highest increases in AOD (+0.03, 13%) and PM 2.5, s (+8 µ g m −3 ). Seasonal AOD patterns fluctuate, with the IGP experiencing the highest wintertime increase, especially in Bihar (+0.07). In summer, there are increases in AOD along the southern and eastern coastal areas. Monsoons cause a slight rise in AOD, except in Rajasthan. In the post-monsoon season, the IGP experiences a notable increase in AOD (+0.057, 25%), potentially driven by biomass burning in Bihar (+0.11) and Uttar Pradesh (+0.075). Dividing our study area into urban and peri-urban clusters ( n = 2791), AOD is found to be similar, possibly due to advective mixing. However, the differences between urban and rural areas become more noticeable, especially in the second decade. Correlations between AOD and PM 2.5, g vary across locations, with the highest found in Kanpur ( R 2 = 0.61) and weaker in Delhi ( R 2 = 0.42), highlighting the need for more ground monitoring. However, it suggests that satellite-derived AOD can generally be used to examine trends in PM 2.5 over longer time frames.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Intensifying heat extremes in China attributed to rising greenhouse gases and declining aerosols since the 2010s

Rising temperatures have increased the frequency, duration, and intensity of extreme temperature events over China in recent decades. The upward trends in heat extremes in China in the warm season (May–September) and their relationships with changes in aerosols and greenhouse gases are investigated using observations, reanalysis data and model results. Significant increasing trends in China are observed in daily maximum temperatures (TXx), heatwave frequency, and heatwave mean duration during 2011–2023, with increasing rates of 0.70 °C/decade, 3.77 d/decade and 0.31 d/event/decade, respectively. This study shows that 43 ± 3% of the TXx increases in China are attributed to the rising CO2 concentrations. Aerosol optical depth in China decreased at a rate of 0.054 per decade from 2011 to 2023 due to significant air quality improvements. The weakened aerosol cooling effect due to declining aerosols contributes 27 ± 3% to the observed TXx trend. In eastern China, where aerosol reductions were most significant, aerosol reduction even accounts for 79 ± 10% of the TXx increasing trend. The intensifying heat waves attributed to GHGs and aerosols are in accordance with the increasing extreme high temperatures. The results highlight the significant impacts of rising GHGs and decreasing aerosols on heat extreme events over China in recent years, emphasizing the need of considering both GHGs and aerosols to address the issue of intensification of heat extremes in a warming future.

Zou, Aili↗
DOE OSTI2025

Multiscale ACI Satellite Database

The SATELLITE_EAGLES_PNNL NetCDF dataset contains a suite of satellite- and reanalysis-derived atmospheric and surface parameters on a regular latitude–longitude grid. The dataset includes core geophysical fields such as land fraction, aerosol optical depth at multiple wavelengths (465, 550, 667, and 865 nm), sea surface temperature, estimated inversion strength, and various thermodynamic and dynamic quantities (e.g., relative humidity, vertical velocity, boundary-layer height, and surface fluxes) from both MERRA and ERA reanalysis products, provided as daily-mean and instantaneous values. A major component of the dataset consists of MODIS-retrieved cloud microphysical properties, including cloud droplet number concentration, cloud effective radius, optical thickness, and liquid water path, provided for three compositing regimes (“All,” “Q06,” and “G18”). Corresponding cloud-top parameters—temperature, height, and pressure—along with total and domain-mean cloud fraction fields are also included. The file further integrates additional satellite data from AMSR-E (for cloud water, rain water, and surface precipitation retrievals) and CERES (for top-of-atmosphere radiative fluxes, cloud fractions, and albedo). This dataset is designed to evaluate aerosol–cloud interactions in warm clouds, emphasizing the use of MODIS for deriving cloud droplet number concentration and liquid water path statistics. The complementary satellite and reanalysis fields are co-located and time-matched to the same instantaneous MODIS observations, enabling consistent comparisons between cloud properties, aerosol loading, and large-scale meteorological conditions. The dataset is recently featured in Christensen et al. (2025), Machine Learning Reveals Strong Grid-Scale Dependence in the Satellite Nd–LWP Relationship, EGUsphere [preprint], https://doi.org/10.5194/egusphere-2025-3850, 2025.

Christensen, Matthew [Pacific Northwest National L↗
DOE OSTI2025

Multiscale ACI Satellite Database

The SATELLITE_EAGLES_PNNL NetCDF dataset contains a suite of satellite- and reanalysis-derived atmospheric and surface parameters on a regular latitude–longitude grid. The dataset includes core geophysical fields such as land fraction, aerosol optical depth at multiple wavelengths (465, 550, 667, and 865 nm), sea surface temperature, estimated inversion strength, and various thermodynamic and dynamic quantities (e.g., relative humidity, vertical velocity, boundary-layer height, and surface fluxes) from both MERRA and ERA reanalysis products, provided as daily-mean and instantaneous values. A major component of the dataset consists of MODIS-retrieved cloud microphysical properties, including cloud droplet number concentration, cloud effective radius, optical thickness, and liquid water path, provided for three compositing regimes (“All,” “Q06,” and “G18”). Corresponding cloud-top parameters—temperature, height, and pressure—along with total and domain-mean cloud fraction fields are also included. The file further integrates additional satellite data from AMSR-E (for cloud water, rain water, and surface precipitation retrievals) and CERES (for top-of-atmosphere radiative fluxes, cloud fractions, and albedo). This dataset is designed to evaluate aerosol–cloud interactions in warm clouds, emphasizing the use of MODIS for deriving cloud droplet number concentration and liquid water path statistics. The complementary satellite and reanalysis fields are co-located and time-matched to the same instantaneous MODIS observations, enabling consistent comparisons between cloud properties, aerosol loading, and large-scale meteorological conditions. The dataset is recently featured in Christensen et al. (2025), Machine Learning Reveals Strong Grid-Scale Dependence in the Satellite Nd–LWP Relationship, EGUsphere [preprint], https://doi.org/10.5194/egusphere-2025-3850, 2025.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Quantifying the impacts of marine aerosols over the southeast Atlantic Ocean using a chemical transport model: implications for aerosol–cloud interactions

The southeast Atlantic region, characterized by persistent stratocumulus clouds, has one of the highest uncertainties in aerosol radiative forcing and significant variability across climate models. In this study, we analyze the seasonally varying role of marine aerosol sources and identify key uncertainties in aerosol composition at cloud-relevant altitudes over the southeast Atlantic using the GEOS-Chem chemical transport model. We evaluate simulated aerosol optical depth (AOD) and speciated aerosol concentrations against those collected from ground observations and aircraft campaigns such as LASIC, ORACLES, and CLARIFY, conducted during 2017. The model consistently underestimates AOD relative to AERONET, particularly at remote locations like Ascension Island. However, when compared with aerosol mass concentrations from aircraft campaigns during the biomass burning period, it performs adequately at cloud-relevant altitudes, with a normalized mean bias (NMB) between -3.5 % (CLARIFY) and -7.5 % (ORACLES). At these altitudes, in the model, organic aerosols (63 %) dominate during the biomass burning period, while sulfate (41 %) prevails during austral summer, when dimethylsulfide (DMS) emissions peak in the model. Our findings indicate that marine sulfate can account for up to 69 % of total sulfate during the high-DMS period. Sensitivity analyses indicate that refining DMS emissions and oxidation chemistry may increase sulfate aerosol produced from marine sources, highlighting that there remains large uncertainty as to the role of DMS emissions in the marine boundary layer. Additionally, we find marine primary organic aerosol emissions may substantially increase total organic aerosol concentrations, particularly during austral summer. This study underscores the imperative need to refine marine emissions and their chemical transformations, as aerosols from marine sources are a major component of total aerosols at cloud-relevant altitudes and may impact uncertainties in aerosol radiative forcing over the southeast Atlantic.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Size-resolved process understanding of stratospheric sulfate aerosol following the Pinatubo eruption

Stratospheric sulfate aerosol produced by volcanic eruptions plays important roles in atmospheric chemistry and the global radiative balance of the atmosphere. The simulation of stratospheric sulfate concentrations and optical properties is highly dependent on the chemistry scheme and microphysical treatment. In this work, we implemented a sophisticated gas-phase chemistry scheme (full chemistry, FC) and a 5-mode version of the Modal Aerosol Module with Prognostic Stratospheric Aerosol (MAM5-PSA) for the interactive treatment of stratospheric sulfate aerosol in the Department of Energy's Energy Exascale Earth System Model version 2 (E3SMv2) model to better simulate the chemistry-aerosol feedback following the Pinatubo eruption, and to compare it against a simulation using simplified chemistry (SC) and the default 4-mode version of the Modal Aerosol Module (MAM4). MAM5-PSA experiments were found to better capture the stratospheric sulfate burden from the eruption of the volcano to the end of 1992 as compared to the High-resolution Infrared Sounder (HIRS) observations, and the formation of sulfate in MAM5-PSA with FC (with an additional OH replenishment reaction) was significantly faster than in MAM4 with FC. Analyses of microphysical processes indicate that more sulfate aerosol mass was generated in total in FC experiments than in SC experiments. MAM5-PSA performs better than MAM4 in simulation of aerosol optical depth (AOD); AOD anomalies from the MAM5-PSA experiment have better agreement with observations. The simulated largest changes in global mean net radiative flux at the top of the atmosphere following the eruption were about −3 W m −2 in MAM5-PSA experiments and roughly −1.5 W m −2 in MAM4 experiments.

AEROSOL↗
DOE OSTI2026

Advancing the quantification of aerosol-cloud interactions with the CALIPSO-CloudSat-Aqua/MODIS record

Aerosol-cloud-precipitation interactions are assessed over the non-polar ocean using more than 11 years of combined Aqua-MODIS, CALIPSO-CALIOP, and CloudSat products. The analysis first shows the benefit of incorporating vertically resolved aerosol extinction coefficient (σext) in aerosol-cloud interactions (ACI) assessments, demonstrating that: σext vertically collocated with the cloud layer () correlates best with cloud droplet number concentration (Nd), column-integrated aerosol optical depth (AOD) cannot explain the Nd variability in the extratropics, and the S-shape of the AOD-Nd relationship reported in previous studies is not replicated when using instead of AOD, with a Nd- linearity more consistent with in-situ studies over the ocean. ACI metric, estimated as the log-scale regression between CALIOP and MODIS Nd reveals that the eastern Pacific is the region with the strongest ACI, followed by the Southern Ocean. The susceptibility of clouds to changes in their liquid water path (LWP) and frequency of precipitation followed a 2-step calculation by combining the Nd- regression (ACI) with the regression between these macrophysical variables and Nd. LWP susceptibility is negative (LWP decreases with aerosol loading) and statistically significant over the eastern Pacific, eastern Atlantic, and extratropics. In contrast, vast areas of the tropical and subtropical ocean feature negligible changes in LWP with aerosol. Precipitation frequency susceptibility is negative, but the values are only significant over the coastal eastern Pacific and Atlantic. The findings suggest that previous modeling assessments relying on AOD may need to be revisited by taking advantage of the synergy between passive and active sensors.

Li, Zhujun↗
DOE OSTI2026

Subsets of geostationary satellite data over international observing network sites for studying the diurnal dynamics of energy, carbon, and water cycles

The latest generation of geostationary satellites provide Earth observations similar to widely used polar-orbiting sensors but at intervals as frequently as every 5–10 min, making them ideal for studying the diurnal dynamics of land–atmosphere interactions. The NASA Earth Exchange (NEX) group created the GeoNEX datasets by collating data from several geostationary platforms, including GOES-16/17/18, Himawari-8/9, and GK-2A, and placing them on a common grid to facilitate use by the Earth science community. Here, we document the GeoNEX Coincident Ground Observations (GeCGO) dataset for terrestrial ecosystem studies and provide examples for its use. Currently, GeCGO provides GOES-16 Advanced Baseline Imager (ABI) data over a 10 km × 10 km area surrounding 1586 network sites across the Americas. GeCGO makes it easy to compare the time series of geostationary data with the diurnal ground observations, including carbon/water fluxes and aerosol optical depth, and is extensible to other regions. We also develop GeoNEXTools to facilitate analyses that require both GeoNEX data and other NASA satellite data. The objectives of this paper are to introduce GeCGO and GeoNEXTools and demonstrate their applications. First, we describe the details of GeCGO and GeoNEXTools. Second, we explain how GeCGO can be integrated with other satellite data. Finally, we showcase comparisons between GeCGO and observations from three ground-based networks. GeCGO is available at https://doi.org/10.25966/y5pe-xp41 (Hashimoto et al., 2025).

Hashimoto, Hirofumi [NASA Ames Research Center (AR↗
DOE OSTI2025

Using feature importance as an exploratory data analysis tool on Earth system models

Abstract. Machine learning (ML) models are commonly used to generate predictions, but these models can also support the discovery of new science. Generating accurate predictions necessitates that a model captures the structure of the underlying data. If the structure is properly extracted, ML could be a useful exploratory and evidential tool. In this paper, we present a case study that demonstrates the use of ML for exploratory data analysis (EDA) in the climate space. We apply the ML explainability method of spatiotemporal zeroed feature importance (stZFI) to understand how climate-variable associations evolve over space and time. Our analyses focus on data from ensembles of Earth system models (ESMs) which provide data on different climate states and conditions. We elect to work with ESM ensembles since they allow us to compare feature importance across alternative scenarios not available with observed data. The ensembles also account for natural variability so that we can distinguish between signal and noise due to natural climate variability when computing feature importance. The use of perturbed initial condition ensembles introduces variability mimicking the natural variability in the atmosphere; thus the signals emerging using feature importance (FI) can be evaluated against the natural variability in the climate system. For our analyses, we consider the 1991 volcanic eruption of Mount Pinatubo, which was a large stratospheric aerosol injection. We explore the climate pathway associated with the eruption from aerosols to radiation to temperature at both the near-surface and stratospheric levels. In addition to applying the method to data generated from two different ESMs, we apply stZFI to reanalysis data to compare the associations identified by stZFI. We show how stZFI tracks the importance of aerosol optical depth over time on forecasting temperatures. This case study illustrates usefulness of an ML tool (stZFI) for EDA on a well-studied climate exemplar.

Ries, Daniel (ORCID:0000000250294647)↗
DOE OSTI2026

Stratospheric aerosol forcing for CMIP7 – Part 1: optical properties for pre-industrial, historical, and scenario simulations

Stratospheric aerosols, most of which originate from explosive volcanic sulfur emissions into the stratosphere, are a key natural driver of climate variability. They are thus a forcing provided by the Coupled Model Intercomparison Project (CMIP) Climate Forcings Task Team to climate modelling groups participating in phase 7 of CMIP. For the historical period, we provide two datasets covering 1750–2023: (i) a volcanic upper tropospheric-stratospheric sulfur emission dataset, documented in a companion paper; and (ii) a stratospheric sulfate aerosol optical property dataset, documented in this paper. For the satellite era (from 1979 onwards), stratospheric aerosol optical properties are derived from the Global Space-based Stratospheric Aerosol Climatology (GloSSAC) dataset. For the pre-satellite era (1750–1978), optical properties are derived from our volcanic SO 2 emission dataset using a new version of the reduced-complexity volcanic aerosol model Easy Volcanic Aerosol (Height) (EVA_H v2). A background, non-volcanic stratospheric aerosol climatology is derived from the 1998–2001 period with a trend over 1850–1978 accounting for increasing anthropogenic aerosols. A monthly stratospheric aerosol climatology is derived from the 1850–2021 average for both pre-industrial and scenario (future) simulations, with a 9-year ramp over 2022–2030 for scenario simulations to ensure a smooth transition from the historical period. CMIP7 simulations should use version 2.2.1 of our pre-industrial and historical datasets, and version 2.2.2 of our scenario dataset. Our methodology to produce historical aerosol optical properties significantly differs from CMIP6 for the pre-satellite era, and the resulting forcings in turn largely differ. In particular, the CMIP6 dataset was mostly based on the sparse and uncertain pyrheliometer record, which resulted in strongly underrepresented emissions from small-to-moderate magnitude eruptions. The resulting bias is addressed in CMIP7, which is entirely emission-derived in the pre-satellite era and uses more recent ice-core-based volcanic sulfur emission inventories than CMIP6. Our approach results in an overall larger volcanic aerosol forcing for CMIP7, with the 1850–2014 mean mid-visible global mean stratospheric aerosol optical depth (SAOD) in CMIP7 (0.0138) being 29 % higher than in CMIP6 (0.0107). The pre-industrial mean of the same variable is 26 % higher in CMIP7 (0.0135, derived from the historical 1850–2021) than CMIP6 (0.0107, derived from the historical 1850–2014 mean). Using a reduced-complexity climate model, we simulate a global mean surface temperature that is 0.07 °C colder for 1850–1900 when using the CMIP7 dataset instead of CMIP6, whereas 2000–2014 is 0.03 °C warmer in CMIP7. Our dataset also exhibits lower forcing for 1960–1980, resulting in temperatures 0.06 °C warmer when averaged over 1960–1990, a period for which CMIP6 climate models exhibit a cold bias. Given the large uncertainties characterizing the dataset, in particular for the pre-satellite era, we advise against treating the CMIP7 or CMIP6 dataset as uniquely superior for any specific year and highlight the need for further evaluation. We conclude the study by discussing sources of uncertainty for the dataset, future research avenues to improve it, as well as requirements to operationalize the production of the dataset, i.e. extend it and update it on an annual basis instead of every 5–7 years following CMIP cycles.

Aubry, Thomas J. [Univ. of Exeter, Devon (United K↗