Impacts of Logistically Constrained High Latitude, Low Altitude Stratospheric Aerosol Injection Scenarios
Poster presentation for ICARP IV Summit (25 - 28 March 2025).
SEARCH · Search NASA
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.
Poster presentation for ICARP IV Summit (25 - 28 March 2025).
Explore the source record for details and available documents.
The 2022 National Defense Strategy of the United States listed climate change as a serious threat to national security. Climate intervention methods, such as stratospheric aerosol injection, have been proposed as mitigation strategies, but the downstream effects of such actions on a complex climate system are not well understood. The development of algorithmic techniques for quantifying relationships between source and impact variables related to a climate event (i.e., a climate pathway) would help inform policy decisions. Data-driven deep learning models have become powerful tools for modeling highly nonlinear relationships and may provide a route to characterize climate variable relationships. In this paper, we explore the use of an echo state network (ESN) for characterizing climate pathways. ESNs are a computationally efficient neural network variation designed for temporal data, and recent work proposes ESNs as a useful tool for forecasting spatiotemporal climate data. However, ESNs are noninterpretable black-box models along with other neural networks. The lack of model transparency poses a hurdle for understanding variable relationships. We address this issue by developing feature importance methods for ESNs in the context of spatiotemporal data to quantify variable relationships captured by the model. We conduct a simulation study to assess and compare the feature importance techniques, and we demonstrate the approach on reanalysis climate data. In the climate application, we consider a time period that includes the 1991 volcanic eruption of Mount Pinatubo. This event was a significant stratospheric aerosol injection, which acts as a proxy for an anthropogenic stratospheric aerosol injection. Furthermore, we are able to use the proposed approach to characterize relationships between pathway variables associated with this event that agree with relationships previously identified by climate scientists.
Three Coupled Model Intercomparison Project 5 (CMIP5) models that simulated the G4 experiment of the Geoengineering Model Intercomparison Project (GeoMIP) were used to investigate the impact of stratospheric aerosol injection (SAI) on combined temperature and precipitation extremes in Africa that can have greater negative impacts on human and the environment than individual rainfall or temperature extremes. The examined compound extremes included the dry (R warm|dry and R cold|dry ) and wet (R warm|wet and R cold|wet ) modes assessed during the injection (SAI, 2050-2069) and post-injection (postSAI, 2070-2089) periods compared with the historical period (1986-2005). We found a significant projected change in the occurrence of both wet and dry modes during SAI and postSAI related to the historical period. The magnitude and sign of this change depend on the season and the geographical location. During the SAI and postSAI, the wet (R warm|wet and R cold|wet ) modes are projected to be significantly lower while the dry modes are noted to increase in a large part of African continent depending on the season and the geographical location and may consequently leads to an increase of the droughts prone areas. Furthermore, the termination effect is noted to re duce the occurrence of dry modes, which may reduce the potential negative effects of the injection after halting. As the effect may vary from one region to another and according to the season, it suggested assessing the key sector impacts of SAI. Thus, this change in dry modes due to SAI could affect all activities which depend on water resources such as water supply, agriculture and food production, energy demand, and production with adverse effects on health, security, and sustainable development, but this needs to be assessed and quantified at regional scales.
Solar radiation modification (SRM) is a possible deliberate approach to decrease or reflect incoming solar radiation with the goal of reducing global temperatures, which have increased over the last decades due to high atmospheric greenhouse gas concentrations. Stratospheric aerosol injection, specifically, has shown potential for successfully reducing global temperatures in climate model simulations. Despite the growing literature in the areas of climate change and SRM, their combined effects on renewable energy generation, a climate change mitigation strategy, have not been addressed. In this review paper, we synthesize previous literature on the possible effects of climate change and SRM on renewable energy resources (i.e., wind energy, solar energy, biomass energy, and hydropower), review the status of climate change and SRM research, and explore potential effects of SRM on renewable energy primarily in the Continental United States (CONUS), but with global perspectives as well. We discuss the research challenges and impacts of SRM on renewable energy and conclude by discussing the potential implications of SRM for renewables for SRM governance and policy. This work is not advocating for or against SRM. It is highlighting an important potential impact for future decision makers.
The National Academies of Sciences, Engineering, and Medicine (NASEM) report "Reflecting Sunlight: Recommendations for Solar Geoengineering Research and Research Governance" focused on three atmospheric-based Solar Geoengineering strategies: Stratospheric Aerosol Injection, Marine Cloud Brightening, Cirrus Cloud Thinning. The study recommends that a national (and internationally-coordinated) Solar Geoengineering research program be established and coordinated by the U.S. Global Change Research Program—with a focus on developing information needed by decision makers, not advancing a path to deployment. The study recommends a detailed framework for governance of this research, at both national and international levels. The study also recommends a detailed research agenda for this program, encompassing over a dozen transdisciplinary research topics that can be loosely grouped into the categories of (i) technical dimensions and impacts, (ii) social dimensions, and (iii) context and goals for solar geoengineering research itself.
The paper constitutes Part 2 of a study performing a first systematic inter-model comparison of the atmospheric responses to stratospheric aerosol injection (SAI) at various single latitudes in the tropics, as simulated by three state-of-the-art Earth system models – CESM2-WACCM6, UKESM1.0, and GISS-E2.1-G. Building on Part 1 (Visioni et al., 2023) we demonstrate the role of biases in the climatological circulation and specific aspects of the model microphysics in driving the inter-model differences in the simulated sulfate distributions. We then characterize the simulated changes in stratospheric and free-tropospheric temperatures, ozone, water vapor, and large-scale circulation, elucidating the role of the above aspects in the surface SAI responses discussed in Part 1. We show that the differences in the aerosol spatial distribution can be explained by the significantly faster shallow branches of the Brewer–Dobson circulation in CESM2, a relatively isolated tropical pipe and older tropical age of air in UKESM, and smaller aerosol sizes and relatively stronger horizontal mixing (thus very young stratospheric age of air) in the two GISS versions used. We also find a large spread in the magnitudes of the tropical lower-stratospheric warming amongst the models, driven by microphysical, chemical, and dynamical differences. These lead to large differences in stratospheric water vapor responses, with significant increases in stratospheric water vapor under SAI in CESM2 and GISS that were largely not reproduced in UKESM. For ozone, good agreement was found in the tropical stratosphere amongst the models with more complex microphysics, with lower stratospheric ozone changes consistent with the SAI-induced modulation of the large-scale circulation and the resulting changes in transport. In contrast, we find a large inter-model spread in the Antarctic ozone responses that can largely be explained by the differences in the simulated latitudinal distributions of aerosols as well as the degree of implementation of heterogeneous halogen chemistry on sulfate in the models. The use of GISS runs with bulk microphysics demonstrates the importance of more detailed treatment of aerosol processes, with contrastingly different stratospheric SAI responses to the models using the two-moment aerosol treatment; however, some problems in halogen chemistry in GISS are also identified that require further attention. Overall, our results contribute to an increased understanding of the underlying physical mechanisms as well as identifying and narrowing the uncertainty in model projections of climate impacts from SAI.
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.
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.
The June 1991 Mt. Pinatubo eruption resulted in a massive increase of sulfate aerosols in the atmosphere, absorbing radiation and leading to global changes in surface and stratospheric temperatures. A volcanic eruption of this magnitude serves as a natural analog for stratospheric aerosol injection, a proposed solar radiation modification method to combat a warming climate. The impacts of such an event are multifaceted and region-specific. Our goal is to characterize the multivariate and dynamic nature of the atmospheric impacts following the Mt. Pinatubo eruption. We developed a multivariate space-time dynamic linear model to understand the full extent of the spatially- and temporally-varying impacts. Specifically, spatial variation is modeled using a flexible set of basis functions for which the basis coefficients are allowed to vary in time through a vector autoregressive (VAR) structure. This novel model is cast in a Dynamic Linear Model (DLM) framework and estimated via a customized MCMC approach. We demonstrate how the model quantifies the relationships between key atmospheric parameters prior to and following the Mt. Pinatubo eruption with reanalysis data from MERRA-2 and highlight when such a model is advantageous over univariate models.
As global temperatures continue to rise, climate mitigation strategies such as stratospheric aerosol injections (SAI) are increasingly discussed, but the downstream effects of these strategies are not well understood. As such, there is interest in developing statistical methods to quantify the evolution of climate variable relationships during the time period surrounding an SAI. Feature importance applied to echo state network (ESN) models has been proposed as a way to understand the effects of SAI using a data-driven model. This approach depends on the ESN fitting the data well. If not, the feature importance may place importance on features that are not representative of the underlying relationships. Typically, time series prediction models such as ESNs are assessed using out-of-sample performance metrics that divide the times series into separate training and testing sets. However, this model assessment approach is geared towards forecasting applications and not scenarios such as the motivating SAI example where the objective is using a data driven model to capture variable relationships. Here, in this paper, we demonstrate a novel use of climate model replicates to investigate the applicability of the commonly used repeated hold-out model assessment approach for the SAI application. Simulations of an SAI are generated using a simplified climate model, and different initialization conditions are used to provide independent training and testing sets containing the same SAI event. The climate model replicates enable out-of-sample measures of model performance, which are compared to the single time series hold-out validation approach. For our case study, it is found that the repeated hold-out sample performance is comparable, but conservative, to the replicate out-of-sample performance when the training set contains enough time after the aerosol injection.
SAND2025-11608O Conditional multi-step attribution for climate forcings and its associated dataset present a framework for attributing a climate impacts, such as stratospheric warming or surface cooling, to a specific magnitude of climate forcing, such as stratospheric aerosol injection. The objective is to improve on existing methods that struggle to attribute highly noisy climate signals. This has potential broader use in climate security applications to identify and understand geoengineering and climate tipping point effects. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.
Abstract Dimension reduction techniques are an essential part of the climate analyst’s toolkit. Due to the enormous scale of climate data, dimension reduction methods are used to identify major patterns of variability within climate dynamics, to create compelling and informative visualizations, and to quantify major named modes such as El Niño–Southern Oscillation. Principal components analysis (PCA), also known as the method of empirical orthogonal functions (EOFs), is the most commonly used form of dimension reduction, characterized by a remarkable confluence of attractive mathematical, statistical, and computational properties. Despite its ubiquity, PCA suffers from several difficulties relevant to climate science: high computational burden with large datasets, decreased statistical accuracy in high dimensions, and difficulties comparing across multiple datasets. In this paper, we introduce several variants of PCA that are likely to be of use in climate sciences and address these problems. Specifically, we introduce non-negative , sparse , and tensor PCA and demonstrate how each approach provides superior pattern recognition in climate data. We also discuss approaches to comparing PCA-family results within and across datasets in a domain-relevant manner. We demonstrate these approaches through an analysis of several runs of the E3SM climate model from 1991 to 1995, focusing on the simulated response to the Mt. Pinatubo eruption; our findings are consistent with a recently identified stratospheric warming fingerprint associated with this type of stratospheric aerosol injection.
Extreme weather events have been demonstrated to be increasing in frequency and intensity across the globe and are anticipated to increase further with projected changes in climate. Solar climate intervention strategies, specifically stratospheric aerosol injection (SAI), have the potential to minimize some of the impacts of a changing climate while more robust reductions in greenhouse gas emissions take effect. However, to date little attention has been paid to the possible responses of extreme weather and climate events under climate intervention scenarios. We present an analysis of 16 extreme surface temperature and precipitation indices, as well as associated vegetation responses, applied to the Geoengineering Large Ensemble (GLENS). GLENS is an ensemble of simulations performed with the Community Earth System Model (CESM1) wherein SAI is simulated to offset the warming produced by a high-emission scenario throughout the 21st century, maintaining surface temperatures at 2020 levels. GLENS is generally successful at maintaining global mean temperature near 2020 levels; however, it does not completely offset some of the projected warming in northern latitudes. Some regions are also projected to cool substantially in comparison to the present day, with the greatest decreases in daytime temperatures. The differential warming–cooling also translates to fewer very hot days but more very hot nights during the summer and fewer very cold days or nights compared to the current day. Extreme precipitation patterns, for the most part, are projected to reduce in intensity in areas that are wet in the current climate and increase in intensity in dry areas. We also find that the distribution of daily precipitation becomes more consistent with more days with light rain and fewer very intense events than currently occur. In many regions there is a reduction in the persistence of long dry and wet spells compared to present day. However, asymmetry in the night and day temperatures, together with changes in cloud cover and vegetative responses, could exacerbate drying in regions that are already sensitive to drought. Overall, our results suggest that while SAI may ameliorate some of the extreme weather hazards produced by global warming, it would also present some significant differences in the distribution of climate extremes compared to the present day.
An ozone lidar system was installed at the National Institute for Environmental Studies (NIES) in Tsukuba, Japan in March 1988 and has been measuring vertical profiles of ozone (15 - 45 km) since September 1988. The lidar system consists of a XeCl (308 nm) excimer laser, its deuterium Raman shifter (339 nm), a XeF excimer laser (351 nm), a 2 m telescope, a receiving system and a data processing system. The precision of the derived ozone concentration is about 10 percent of an altitude of 40 km for a 4 hr observation. Temperature profiles (30 - 80 km) are also obtained from the Rayleigh scattering signals at 351 nm. Approximate 50 ozone measurements are carried out in a year and variations of vertical profiles of ozone such as seasonal variations and shorter-term variations are observed. Systematic errors due to aerosols had been negligible until the arrival of the stratospheric aerosols injected by the eruption of Mt. Pinatubo. Effects of the volcanic aerosols on ozone measurements depend on the differences between wavelengths used as the on- and off-resonance.
Stratospheric aerosol geoengineering focused on the Arctic could substantially reduce local and worldwide impacts of anthropogenic global warming. Because the Arctic receives little sunlight during the winter, stratospheric aerosols present in the winter at high latitudes have little impact on the climate, whereas stratospheric aerosols present during the summer achieve larger changes in radiative forcing. Injecting SO 2 in the spring leads to peak aerosol optical depth (AOD) in the summer. We demonstrate that spring injection produces approximately twice as much summer AOD as year-round injection and restores approximately twice as much September sea ice, resulting in less increase in stratospheric sulfur burden, stratospheric heating, and stratospheric ozone depletion per unit of sea ice restored. Furthermore, we also find that differences in AOD between different seasonal injection strategies are small compared to the difference between annual and spring injection.
Simulations of stratospheric geoengineering with black carbon (BC) aerosols using a general circulation model with fixed sea surface temperatures show that the climate effects strongly depend on aerosol size and altitude of injection. 1 Tg BC/a injected into the lower stratosphere would cause little surface cooling for large radii but a large amount of surface cooling for small radii and stratospheric warming of over 60 C. With the exception of small particles, increasing the altitude of injection increases surface cooling and stratospheric warming. Stratospheric warming causes global ozone loss by up to 50% in the small radius case. The Antarctic shows less ozone loss due to reduction of polar stratospheric clouds, but strong circumpolar winds would enhance the Arctic ozone hole. Using diesel fuel to produce the aerosols is likely prohibitively expensive and infeasible. Although studying an absorbing aerosol is a useful counterpart to previous studies involving sulfate aerosols, black carbon geoengineering likely carries too many risks to make it a viable option for deployment.
The Pacific Northwest pyrocumulonimbus Event (PNE) took place in British Columbia during the evening and night-time hours between the 12th and 13th of August 2017. Several pyroconvective clouds erupted in this occasion, and released in the upper troposphere/lower stratosphere unprecedented amounts of carbonaceous aerosols (300 ktn). Only a few years later, an even larger pyroCb injection took place over Australia. This event, named the Australian New Year (ANY) event, injected up to 1100 ktn of aerosol between December 29th 2019 and January 4th 2020. Such large injections of carbonaceous aerosol modify the stratospheric radiative budgets, locally perturbing stratospheric temperatures and winds. In this study, we use the Goddard Earth Observing System Chemistry Climate Model (GEOS CCM) to study the pertubations on the stratospheric meteorology induced by an aerosol injection of the magnitude of the PNE. Our simulations include the radiative interactions of aerosols, so that their impact on temperatures and winds are explicitly simulated. We show how the presence of the carbonaceous aerosols from the pyroCb causes the formation and maintenance of a synoptic scale stratospheric anticyclone. We follow this disturbance considering the potential vorticity anomaly and the brown carbon aerosol loading and we describe its dynamical and thermodynamical structure and its evolution in time. The analysis presented here shows that the simulated anticyclone undergoes daily expansion-compression cycles governed by the radiative heating, which are directly related to the vertical motion of the plume, and that the aerosol radiative heating is essential in maintaining the anticyclone itself.