Search

There is now substantial literature on climate model studies of equatorial or tropical stratospheric SO2 injections that aim to counteract the surface warming produced by rising concentrations of greenhouse gases. Here we present the results from the first systematic intercomparison of climate responses in three Earth system models wherein the injection of SO2 occurs at different latitudes in the lower stratosphere: CESM2-WACCM6, UKESM1.0 and GISS-E2.1-G. The first two use a modal aerosol microphysics scheme, while two versions of GISS-E2.1-G use a bulk aerosol (One-Moment Aerosol, OMA) and a two-moment (Multiconfiguration Aerosol TRacker of mIXing state, MATRIX) microphysics approach, respectively. Our aim in this work is to determine commonalities and differences between the climate model responses in terms of the distribution of the optically reflective sulfate aerosols produced from the oxidation of SO2 and in terms of the surface response to the resulting reduction in solar radiation. A focus on understanding the contribution of characteristics of models transport alongside their microphysical and chemical schemes, and on evaluating the resulting stratospheric responses in different models, is given in the companion paper (Bednarz et al., 2023). The goal of this exercise is not to evaluate these single-point injection simulations as stand-alone proposed strategies to counteract global warming; instead we determine sources and areas of agreement and uncertainty in the simulated responses and, ultimately, the possibility of designing a comprehensive intervention strategy capable of managing multiple simultaneous climate goals through the combination of different injection locations. We find large disagreements between GISS-E2.1-G and the CESM2-WACCM6 and UKESM1.0 models regarding the magnitude of cooling per unit of aerosol optical depth (AOD) produced, which varies from 4.7 K per unit of AOD in CESM2-WACCM6 to 16.7 K in the GISS-E2.1-G version with two-moment aerosol microphysics. By normalizing the results with the global mean response in each of the models and thus assuming that the amount of SO2 injected is a free parameter that can be managed independently, we highlight some commonalities in the overall distributions of the aerosols, in the inter-hemispheric surface temperature response and in shifts to the Intertropical Convergence Zone, as well as some areas of disagreement, such as the extent of the aerosol confinement in the equatorial region and the efficiency of the transport to polar latitudes. In conclusion, we demonstrate that it is possible to use these simulations to produce more comprehensive injection strategies in multiple climate models. However, large differences in the injection magnitudes can be expected, potentially increasing inter-model spreads in some stratospheric quantities (such as aerosol distribution) while reducing the spread in the surface response in terms of temperature and precipitation; furthermore, the selection of the injection locations may be dependent on the models' specific stratospheric transport.

We present here results from the Geoengineering Model Intercomparison Project (GeoMIP) simulations for the experiments G6sulfur and G6solar for six Earth system models participating in the Climate Model Intercomparison Project (CMIP) Phase 6. The aim of the experiments is to reduce the warming that results from a high-tier emission scenario (Shared Socioeconomic Pathways SSP5-8.5) to that resulting from a medium-tier emission scenario (SSP2-4.5). These simulations aim to analyze the response of climate models to a reduction in incoming surface radiation as a means to reduce global surface temperatures, and they do so either by simulating a stratospheric sulfate aerosol layer or, in a more idealized way, through a uniform reduction in the solar constant in the model. We find that over the final two decades of this century there are considerable inter-model spreads in the needed injection amounts of sulfate (29 ± 9 Tg-SO2/yr between 2081 and 2100), in the latitudinal distribution of the aerosol cloud and in the stratospheric temperature changes resulting from the added aerosol layer. Even in the simpler G6solar experiment, there is a spread in the needed solar dimming to achieve the same global temperature target (1.91 ± 0.44 %). The analyzed models already show significant differences in the response to the increasing CO2 concentrations for global mean temperatures and global mean precipitation (2.05 K ± 0.42 K and 2.28 ± 0.80 %, respectively, for SSP5-8.5 minus SSP2-4.5 averaged over 2081–2100). With aerosol injection, the differences in how the aerosols spread further change some of the underlying uncertainties, such as the global mean precipitation response (−3.79 ± 0.76 % for G6sulfur compared to −2.07 ± 0.40 % for G6solar against SSP2-4.5 between 2081 and 2100). These differences in the behavior of the aerosols also result in a larger uncertainty in the regional surface temperature response among models in the case of the G6sulfur simulations, suggesting the need to devise various, more specific experiments to single out and resolve particular sources of uncertainty. The spread in the modeled response suggests that a degree of caution is necessary when using these results for assessing specific impacts of geoengineering in various aspects of the Earth system. However, all models agree that compared to a scenario with unmitigated warming, stratospheric aerosol geoengineering has the potential to both globally and locally reduce the increase in surface temperatures.

The paper constitutes part 2 of a study performing a first systematic inter-model comparison of the atmospheric responses to stratospheric sulfate aerosol injections (SAI) at various latitudes as simulated by three state-of-the-art Earth System Models – CESM2(WACCM6), UKESM1.0, and GISS-E2.1-G. We use a set of five sensitivity experiments with constant annual injections of SO2 in the lower stratosphere at either 30° S, 15° S, 0°, 15° N or 30° N. We identify the similarities and differences in the simulated responses amongst the models as well as demonstrate the role of biases in the climatological circulation and specific aspects of the model microphysics and chemistry in driving the inter-model differences. Building on part 1 (Visioni et al., 2022), we explain the simulated differences in the aerosol spatial distribution between the models: CESM2 shows significantly faster shallow branches of the Brewer Dobson circulation facilitating transport of the relatively larger-sized aerosol to higher latitudes; UKESM shows a relatively isolated tropical pipe and older tropical age-of-air confining the relatively smaller-sized aerosols to the tropics; and the two GISS versions with either bulk or modal aerosol microphysics show elevated sulfate levels at higher latitudes as the result of smaller aerosol sizes and relatively stronger horizontal mixing (thus very young stratospheric age-of-air). We then elucidate the role of these factors in driving the stratospheric responses to SAI. We find a large spread in the magnitudes of the tropical lower stratospheric warming amongst the models, which can be partially attributed to the differences in aerosol distribution and sizes. Regarding the stratospheric ozone responses, we find a good agreement in the tropics between the models with modal 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 to the relative agreement at low latitudes, 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. Finally, we also find large differences in stratospheric water vapour responses amongst the models, with CESM2 and GISS with modal microphysics both showing significant increases in stratospheric water vapour under SAI consistent with the increase in cold point temperatures that were largely not reproduced in UKESM. For the GISS runs with bulk microphysics, the SAI simulations show contrastingly different stratospheric responses to the models using the modal aerosol treatment, including the absence of lower stratospheric warming as well as significant reductions in stratospheric water vapour and ozone. The results point towards the importance of detailed treatment of aerosol processes, although some problems in halogen chemistry in GISS are identified that require further attention. Overall, our results contribute to an increased understanding of the underlying physical mechanisms as well as the sources of uncertainty in model projections of climate impacts from SAI.

Stratospheric aerosol geoengineering has been proposed as a potential solution to reduce climate change and its impacts. Here, we explore the responses of the Hadley circulation (HC) intensity and the intertropical convergence zone (ITCZ) using the strategic stratospheric aerosol geoengineering, in which sulfur dioxide was injected into the stratosphere at four different locations to maintain the global-mean surface temperature and the interhemispheric and equator-to-pole temperature gradients at present-day values (baseline). Simulations show that, relative to the baseline, strategic stratospheric aerosol geoengineering generally maintains northern winter December–January–February (DJF) HC intensity under RCP8.5, while it overcompensates for the greenhouse gas (GHG)-forced southern winter June–July–August (JJA) HC intensity increase, producing a 3.5 ± 0.4% weakening. The residual change of southern HC intensity in JJA is mainly associated with stratospheric heating and tropospheric temperature response due to enhanced stratospheric aerosol concentrations. Geoengineering overcompensates for the GHG-driven northward ITCZ shifts, producing 0.7° ± 0.1° and 0.2° ± 0.1° latitude southward migrations in JJA and DJF, respectively relative to the baseline. These migrations are affected by tropical interhemispheric temperature differences both at the surface and in the free troposphere. Further strategies for reducing the residual change of HC intensity and ITCZ shifts under stratospheric aerosol geoengineering could involve minimizing stratospheric heating and restoring and preserving the present-day tropical tropospheric interhemispheric temperature differences.

Solar geoengineering is an emerging topic in climate-change discussions. To support future decisions on the deployment of this technology, society requires better estimates of its environmental impacts and limitations. As solar geoengineering has never been deployed, conclusions about its climatic effects are primarily obtained through models and natural analogues. As such, our confidence in projections of solar geoengineering, the basis for that confidence and how our confidence can be improved is limited. In this Perspective, we review our current understanding of uncertainty and risk in solar geoengineering via stratospheric aerosols. Using a risk-register framework, we illustrate key uncertainties, such as sub-grid-scale mixing or effects of stratospheric heating, investigations of which should be prioritized to transition the field to a mission-driven research agenda. We conclude with recommendations for possible avenues of research, including targeted model intercomparisons and appropriately governed small-scale field experiments. Discussions surrounding stratospheric aerosol geoengineering have increased in climate-change discourse, but our understanding of its impacts is relatively limited. This Perspective discusses the uncertainty and risk associated with stratospheric aerosol geoengineering, offering suggestions on research topics that should be prioritized.

Consideration of solar geoengineering as a potential response to climate change will demand complex decisions. These include not only the choice of whether to deploy solar engineering, but decision...

Stratospheric aerosol injection (SAI) is a prospective climate intervention technology that would seek to abate climate change by deflecting back into space a small fraction of the incoming solar radiation. While most consideration given to SAI assumes a global intervention, this paper considers an alternative scenario whereby SAI might be deployed only in the subpolar regions. Subpolar deployment would quickly envelope the poles as well and could arrest or reverse ice and permafrost melt at high latitudes. This would yield global benefit by retarding sea level rise. Given that effective SAI deployment could be achieved at much lower altitudes in these regions than would be required in the tropics, it is commonly assumed that subpolar deployment would present fewer aeronautical challenges. An SAI deployment intended to reduce average surface temperatures in both the Arctic and Antarctic regions by 2 °C is deemed here to be feasible at relatively low cost with conventional technologies. However, we do not find that such a deployment could be undertaken with a small fleet of pre-existing aircraft, nor that relegating such a program to these sparsely populated regions would obviate the myriad governance challenges that would confront any such deployment. Nevertheless, given its feasibility and potential global benefit, the prospect of subpolar-focused SAI warrants greater attention.

Deliberately blocking out a small portion of the incoming solar radiation would cool the climate. One such approach would be injecting SO$_2$ into the stratosphere, which would produce sulfate aero...

Most socio-economic activities in Africa depend on the continent’s river basins, but effectively managing drought risks over the basins in response to climate change remains a big challenge. While studies have shown that the stratospheric aerosol injection (SAI) intervention could mitigate temperature-related climate change impacts over Africa, there is a dearth of information on how the SAI intervention could influence drought characteristics and drought risk managements over the river basins. The present study thus examines the potential impacts of climate change and the SAI intervention on droughts and drought management over the major river basins in Africa. Multi-ensemble climate simulation datasets from the Stratospheric Aerosol Geoengineering Large Ensemble (GLENS) Project were analysed for the study. The Standardized Precipitation Evapotranspiration Index (SPEI) and the Standardized Precipitation Index (SPI) were used to characterize the upper and lower limits of future drought severity, respectively, over the basins. The SPEI is a function of rainfall and potential evapotranspiration, whereas the SPI is only a function of rainfall, so the difference between the two indices is influenced by atmospheric evaporative demand. The results of the study show that, while the SAI intervention, as simulated in GLENS, may offset the impacts of climate change on temperature and atmospheric evaporative demand, the level of SAI that compensates for temperature change would overcompensate for the impacts on precipitation and therefore impose a climate water balance deficit in the tropics. SAI would narrow the gaps between SPEI and SPI projections over the basins by reducing SPEI drought frequency through reduced temperature and atmospheric evaporative demand while increasing SPI drought frequency through reduced rainfall. The narrowing of this gap lowers the level of uncertainty regarding future changes in drought frequency, but nonetheless has implications for future drought management in the basins, because while SAI lowers the upper limit of the future drought stress, it also raises the lower limit of the drought stress.

This paper describes the Stratospheric Aerosol Geoengineering Large Ensemble (GLENS) project, which promotes the use of a unique model dataset, performed with the Community Earth System Model, with the Whole Atmosphere Community Climate Model as its atmospheric component [CESM1(WACCM)], to investigate global and regional impacts of geoengineering. The performed simulations were designed to achieve multiple simultaneous climate goals, by strategically placing sulfur injections at four different locations in the stratosphere, unlike many earlier studies that targeted globally averaged surface temperature by placing injections in regions at or around the equator. This advanced approach reduces some of the previously found adverse effects of stratospheric aerosol geoengineering, including uneven cooling between the poles and the equator and shifts in tropical precipitation. The 20-member ensemble increases the ability to distinguish between forced changes and changes due to climate variability in global and regional climate variables in the coupled atmosphere, land, sea ice, and ocean system. We invite the broader community to perform in-depth analyses of climate-related impacts and to identify processes that lead to changes in the climate system as the result of a strategic application of stratospheric aerosol geoengineering.

We describe here changes in stratospheric dynamics and chemistry in a first century‐long sulfate aerosol geoengineering simulation in which the mean surface temperature and the interhemispheric and equator‐to‐pole surface temperature gradients were kept near their 2020 levels despite the RCP8.5 emission scenario. Simulations were carried out with the Community Earth System Model, version 1 with the Whole Atmosphere Community Climate Model as its atmospheric component [CESM1(WACCM)] coupled to a feedback algorithm controlling the magnitude of sulfur dioxide (SO_2) injections at four injection latitudes. We find that, throughout the entire geoengineering simulation, the lower stratospheric temperatures increase by ∼0.19 K per Tg SO_2 injection per year or ∼10 K with ∼40 Tg SO_2/year total SO_2 injection. These temperature changes are associated with a strengthening of the polar jets in the stratosphere and weakening of the mean zonal wind in the lower stratosphere subtropics and throughout the troposphere, associated with weaker storm track activity. In the geoengineering simulation the quasi‐biennial oscillation of the tropical lower stratospheric winds remains close to the presently observed quasi‐biennial oscillation, even for large amounts of SO2 injection. Water vapor in the stratosphere increases substantially: by 25% with ∼20 Tg SO_2/year annual injection and by up to 90% with a ∼40 Tg SO_2/year injection. Stratospheric column ozone in the geoengineering simulation is predicted to recover to or supersede preozone hole conditions by the end of the century.

Strategically applied geoengineering is proposed to reduce some of the known side effects of stratospheric aerosol modifications. Specific climate goals could be reached depending on design choices of stratospheric sulfur injections by latitude, altitude, and magnitude. Here we explore in detail the stratospheric chemical and dynamical responses to injections at different altitudes using a fully coupled Earth System Model. Two different scenarios are explored that produce approximately the same global cooling of 2°C over the period 2042–2049, a high‐altitude injection case using 24 Tg SO_2/year at 30 hPa (≈25‐km altitude) and a low‐altitude injection case using 32 Tg SO_2/year injections at 70 hPa (between 19‐ and 20‐km altitude), with annual injections divided equally between 15°N and 15°S. Both cases result in a warming of the lower tropical stratosphere up to 10 and 15°C for the high‐ and low‐altitude injection case and in substantial increases of stratospheric water vapor of up to 2 and 4 ppm, respectively, compared to no geoengineering conditions. Polar column ozone in the Northern Hemisphere is reduced by up to 18% in March for the high‐altitude injection case and up to 8% for the low‐altitude injection case. However, for winter middle and high northern latitudes, low‐altitude injections result in greater column ozone values than without geoengineering. These changes are mostly driven by dynamics and advection. Antarctic column ozone in 2042–2049 does not recover from present‐day (2002–2009) values for both cases.

Geoengineering methods could potentially offset aspects of greenhouse gas‐driven climate change. However, before embarking on any such strategy, a comprehensive understanding of its impacts must be obtained. Here, a 20‐member ensemble of simulations with the Community Earth System Model with the Whole Atmosphere Community Climate Model as its atmospheric component is used to investigate the projected hydroclimate changes that occur when greenhouse gas‐driven warming, under a high emissions scenario, is offset with stratospheric aerosol geoengineering. Notable features of the late 21st century hydroclimate response, relative to present day, include a reduction in precipitation in the Indian summer monsoon, over much of Africa, Amazonia and southern Chile and a wintertime precipitation reduction over the Mediterranean. Over most of these regions, the soil desiccation that occurs with global warming is, however, largely offset by the geoengineering. A notable exception is India, where soil desiccation and an approximate doubling of the likelihood of monsoon failures occurs. The role of stratospheric heating in the simulated hydroclimate change is determined through additional experiments where the aerosol‐induced stratospheric heating is imposed as a temperature tendency, within the same model, under present day conditions. Stratospheric heating is found to play a key role in many aspects of projected hydroclimate change, resulting in a general wet‐get‐drier, dry‐get‐wetter pattern in the tropics and extratropical precipitation changes through midlatitude circulation shifts. While a rather extreme geoengineering scenario has been considered, many, but not all, of the precipitation features scale linearly with the offset global warming.

Multidecadal variability in the Atlantic meridional overturning circulation (AMOC) is shown to differ significantly between the 4 × CO2 and preindustrial control simulations of the GFDL Earth System Model, version 2M (ESM2M) general circulation model (GCM). In the preindustrial simulation, this model has a peak in the power spectrum of both AMOC and northward heat transport at latitudes between 26° and 50°N. In the 4 × CO2 simulation, the only significant spectral peak is near 60°N. Understanding these differences is important for understanding the effect of future climate change on climate variability, as well as for providing insight into the physics underlying AMOC variability. Transfer function analysis demonstrates that the shift is predominantly due to a shift in the internal ocean dynamics rather than a change in stochastic atmospheric forcing. Specifically, the reduction in variance from 26° to 45°N is due to an increased stratification east of Newfoundland that results from the shallower and weaker mean overturning. The reduced AMOC variance that accompanies the reduced mean value of the AMOC at 4 × CO2 differs from predictions of simple box models that predict a weaker circulation to be closer to a stability bifurcation point and, therefore, be accompanied by amplified variability. The high-latitude variability in the 4 × CO2 simulation is related to the advection of anomalies by the subpolar gyre, distinct from the variability mechanism in the control simulation at lower latitudes. The 4 × CO2 variability has only a small effect on midlatitude meridional heat transport, but does significantly affect sea ice in the northern North Atlantic.

Solar Radiation Management (SRM) geoengineering is a proposed response to anthropogenic global warming (AGW). Stratospheric aerosol injection (SAI) is one proposed method, reliant on lofting material into the stratosphere. Engineering reviews related to this technology approach have been sparse, with most major primary analyses now at least five years old. We attempt to bridge this gap—with a short, qualitative review of recent developments in various fields of engineering that have potential applicability to SAI. Our analysis shows that a new conventional aircraft design is still likely to be the most dependable and affordable technology solution (cost estimates start around $1000–1500 per ton lofted), with hybrid or vacuum airships a potential challenger. Rockets, gas guns and MAGLEV/coilguns show some potential—although they lack the inherent level-flight capability that would be needed for direct aerosol distribution (versus distribution of gaseous precursors), without substantial additional engineering. Should very high-altitude access be required, rockets, jet-hybrid rockets, and various guns (especially light-gas guns) potentially offer the required capability. Costs and performance for tethered balloons remain highly uncertain. Towers are not found to be promising. The extreme accessibility of free balloons suggests that this method may be used primarily for reasons of political leverage, as opposed to being an optimal engineering solution.

The authors describe a new method of comparing different climate forcing agents (e.g., CO2 concentration, CH4 concentration, and total solar irradiance) in climate models that circumvents many of the difficulties associated with explicit calculations of efficacy. This is achieved by introducing an explicit feedback loop external to a climate model that adjusts one forcing agent to balance another while keeping global-mean surface temperature constant. The convergence time of this feedback loop can be adjusted, allowing for comparisons of forcing agents to be achieved with relatively short simulations. Comparisons between forcing agents are highly linear in concordance with predicted scaling relationships; for example, the global-mean climate response to a doubling of the CO2 concentration is equivalent to that of a 2.1% change in total solar irradiance. This result is independent of the magnitude of the forcing agent (within the range of radiative forcings considered here) and is consistent across two different climate models.

Injections of sulfur dioxide into the stratosphere are among several proposed methods of solar radiation management. Such injections could cool the Earth's climate. However, they would significantly alter the dynamics of the stratosphere. We explore here the stratospheric dynamical response to sulfur dioxide injections ∼ 5 km above the tropopause at multiple latitudes (equator, 15° S, 15° N, 30° S and 30° N) using a fully coupled Earth system model, Community Earth System Model, version 1, with the Whole Atmosphere Community Climate Model as its atmospheric component (CESM1(WACCM)). We find that in all simulations, the tropical lower stratosphere warms primarily between 30° S and 30° N, regardless of injection latitude. The quasi-biennial oscillation (QBO) of the tropical zonal wind is altered by the various sulfur dioxide injections. In a simulation with a 12 Tg yr−1 equatorial injection, and with fully interactive chemistry, the QBO period lengthens to ∼ 3.5 years, but never completely disappears. However, in a simulation with specified (or non-interactive) chemical fields, including O3 and prescribed aerosols taken from the interactive simulation, the oscillation is virtually lost. In addition we find that geoengineering does not always lengthen the QBO. We further demonstrate that the QBO period changes from 24 to 12 - 17 months in simulations with sulfur dioxide injections placed poleward of the equator. Our study points to the importance of understanding and verifying of the complex interactions between aerosols, atmospheric dynamics, and atmospheric chemistry as well as understanding the effects of sulfur dioxide injections placed away from the Equator on the QBO.

We describe the first simulations of stratospheric sulfate aerosol geoengineering using multiple injection locations to meet multiple simultaneous surface temperature objectives. Simulations were performed using CESM1(WACCM), a coupled atmosphere-ocean general circulation model with fully interactive stratospheric chemistry, dynamics (including an internally generated quasi-biennial oscillation), and a sophisticated treatment of sulfate aerosol formation, microphysical growth, and deposition. The objectives are defined as maintaining three temperature features at their 2020 levels against a background of the RCP8.5 scenario over the period 2020–2099. These objectives are met using a feedback mechanism in which the rate of sulfur dioxide injection at each of the four locations is adjusted independently every year of simulation. Even in the presence of uncertainties, nonlinearities, and variability, the objectives are met, predominantly by SO_2 injection at 30°N and 30°S. By the last year of simulation, the feedback algorithm calls for a total injection rate of 51 Tg SO_2 per year. The injections are not in the tropics, which results in a greater degree of linearity of the surface climate response with injection amount than has been found in many previous studies using injection at the equator. Because the objectives are defined in terms of annual mean temperature, the required geongineering results in “overcooling” during summer and “undercooling” during winter. The hydrological cycle is also suppressed as compared to the reference values corresponding to the year 2020. The demonstration we describe in this study is an important step toward understanding what geoengineering can do and what it cannot do.

We present new insights into the evolution and interactions of stratospheric aerosol using an updated version of the Whole Atmosphere Community Climate Model (WACCM). Improved horizontal resolution, dynamics, and chemistry now produce an internally generated quasi-biennial oscillation and significant improvements to stratospheric temperatures and ozone compared to observations. We present a validation of WACCM column ozone and climate calculations against observations. The prognostic treatment of stratospheric sulfate aerosols accurately represents the evolution of stratospheric aerosol optical depth and perturbations to solar and longwave radiation following the June 1991 eruption of Mount Pinatubo. We confirm the inclusion of interactive OH chemistry as an important factor in the formation and initial distribution of aerosol following large inputs of sulfur dioxide (SO_2) to the stratosphere. We calculate that depletion of OH levels within the dense SO_2 cloud in the first weeks following the Pinatubo eruption significantly prolonged the average initial e-folding decay time for SO_2 oxidation to 47 days. Previous observational and model studies showing a 30 day decay time have not accounted for the large (30–55%) losses of SO_2 on ash and ice within 7–9 days posteruption and have not correctly accounted for OH depletion. We examine the variability of aerosol evolution in free-running climate simulations due to meteorology, with comparison to simulations nudged with specified dynamics. We assess calculated impacts of volcanic aerosols on ozone loss with comparisons to observations. The completeness of the chemistry, dynamics, and aerosol microphysics in WACCM qualify it for studies of stratospheric sulfate aerosol geoengineering.

Injection of SO_2 into the stratosphere has been proposed as a method to, in part, counteract anthropogenic climate change. So far, most studies investigated injections at the equator or in a region in the tropics. Here we use Community Earth System Model version 1 Whole Atmosphere Community Climate Model (CESM1(WACCM)) to explore the impact of continuous single grid point SO_2 injections at seven different latitudes and two altitudes in the stratosphere on aerosol distribution and climate. For each of the 14 locations, 3 different constant SO_2 emission rates were tested to identify linearity in aerosol burden, aerosol optical depth, and climate effects. We found that injections at 15°N and 15°S and at 25 km altitude have equal or greater effect on radiation and surface temperature than injections at the equator. Nonequatorial injections transport SO_2 and sulfate aerosols more efficiently into middle and high latitudes and result in particles of smaller effective radius and larger aerosol burden in middle and high latitudes. Injections at 15°S produce the largest increase in global average aerosol optical depth and increase the change in radiative forcing per Tg SO_2/yr by about 15% compared to equatorial injections. High-altitude injections at 15°N produce the largest reduction in global average temperature of 0.2° per Tg S/yr for the last 7 years of a 10 year experiment. Injections at higher altitude are generally more efficient at reducing surface temperature, with the exception of large equatorial injections of at least 12 Tg SO_2/yr. These findings have important implications for designing a strategy to counteract global climate change.

By injecting different amounts of SO_2 at multiple different latitudes, the spatial pattern of aerosol optical depth (AOD) can be partially controlled. This leads to the ability to influence the climate response to geoengineering with stratospheric aerosols, providing the potential for design. We use simulations from the fully coupled whole-atmosphere chemistry climate model CESM1(WACCM) to demonstrate that by appropriately combining injection at just four different locations, 30°S, 15°S, 15°N, and 30°N, then three spatial degrees of freedom of AOD can be achieved: an approximately spatially uniform AOD distribution, the relative difference in AOD between Northern and Southern Hemispheres, and the relative AOD in high versus low latitudes. For forcing levels that yield 1–2°C cooling, the AOD and surface temperature response are sufficiently linear in this model so that the response to different combinations of injection at different latitudes can be estimated from single-latitude injection simulations; nonlinearities associated with both aerosol growth and changes to stratospheric circulation will be increasingly important at higher forcing levels. Optimized injection at multiple locations is predicted to improve compensation of CO_2-forced climate change relative to a case using only equatorial aerosol injection (which overcools the tropics relative to high latitudes). The additional degrees of freedom can be used, for example, to balance the interhemispheric temperature gradient and the equator to pole temperature gradient in addition to the global mean temperature. Further research is needed to better quantify the impacts of these strategies on changes to long-term temperature, precipitation, and other climate parameters.

If solar radiation management (SRM) were ever implemented, feedback of the observed climate state might be used to adjust the radiative forcing of SRM in order to compensate for uncertainty in either the forcing or the climate response. Feedback might also compensate for unexpected changes in the system, e.g. a nonlinear change in climate sensitivity. However, in addition to the intended response to greenhouse-gas induced changes, the use of feedback would also result in a geoengineering response to natural climate variability. We use a box-diffusion dynamic model of the climate system to understand how changing the properties of the feedback control affect the emergent dynamics of this coupled human–climate system, and evaluate these predictions using the HadCM3L general circulation model. In particular, some amplification of natural variability is unavoidable; any time delay (e.g., to average out natural variability, or due to decision-making) exacerbates this amplification, with oscillatory behavior possible if there is a desire for rapid correction (high feedback gain). This is a challenge for policy as a delayed response is needed for decision making. Conversely, the need for feedback to compensate for uncertainty, combined with a desire to avoid excessive amplification of natural variability, results in a limit on how rapidly SRM could respond to changes in the observed state of the climate system.

Solar radiation management (SRM) has been proposed as a possible option for offsetting some anthropogenic radiative forcing, with the goal of reducing some of the associated climatic changes. There are clearly significant uncertainties associated with SRM, and even small-scale experiments that might reduce uncertainty would carry some risk. However, there are also natural and anthropogenic analogs to SRM, such as volcanic eruptions in the case of stratospheric aerosol injection and ship tracks in the case of marine cloud albedo modification. It is essential to understand what we can learn from these analogs in order to validate models, particularly because of the problematic nature of outdoor experiments. It is also important to understand what we cannot learn, as this might better focus attention on what risks would need to be solely examined by numerical models. Stratospheric conditions following a major volcanic eruption, for example, are not the same as those to be expected from intentional geoengineering, both because of confounding effects of volcanic ash and the differences between continuous and impulsive injection of material into the stratosphere. Nonetheless, better data would help validate models; we thus recommend an appropriate plan be developed to better monitor the next large volcanic eruption. Similarly, more could be learned about cloud albedo modification from careful study not only of ship tracks, but of ship and other aerosol emission sources in cloud regimes beyond the narrow conditions under which ship tracks form; this would benefit from improved satellite observing capabilities.

Climate emulators trained on existing simulations can be used to project project the climate effects that result from different possible future pathways of anthropogenic forcing, without further relying on general circulation model (GCM) simulations. We extend this idea to include different amounts of solar geoengineering in addition to different pathways of greenhouse gas concentrations, by training emulators from a multi-model ensemble of simulations from the Geoengineering Model Intercomparison Project (GeoMIP). The emulator is trained on the abrupt 4 × CO2 and a compensating solar reduction simulation (G1), and evaluated by comparing predictions against a simulated 1 % per year CO2 increase and a similarly smaller solar reduction (G2). We find reasonable agreement in most models for predicting changes in temperature and precipitation (including regional effects), and annual-mean Northern Hemisphere sea ice extent, with the difference between simulation and prediction typically being smaller than natural variability. This verifies that the linearity assumption used in constructing the emulator is sufficient for these variables over the range of forcing considered. Annual-minimum Northern Hemisphere sea ice extent is less well predicted, indicating a limit to the linearity assumption.

We introduce system identification techniques to climate science wherein multiple dynamic input-output relationships can be simultaneously characterized in a single simulation. This method, involving multiple small perturbations (in space and time) of an input field while monitoring output fields to quantify responses, allows for identification of different timescales of climate response to forcing without substantially pushing the climate far away from a steady state. We use this technique to determine the steady state responses of low cloud fraction and latent heat flux to heating perturbations over 22 regions spanning Earth's oceans. We show that the response characteristics are similar to those of step-change simulations, but in this new method, the responses for 22 regions can be characterized simultaneously. Furthermore, we can estimate the timescale over which the steady state response emerges. The proposed methodology could be useful for a wide variety of purposes in climate science, including characterization of teleconnections and uncertainty quantification to identify the effects of climate model tuning parameters.

Systematic and compensating errors can lead to degraded predictive skill in climate models. Such errors may be identified by comparing different models in an analysis of individual physical processes. We examine model simulations of El Nino–Southern Oscillation (ENSO) in five Coupled Model Intercomparison Project (CMIP) models, using transfer functions to analyze nine processes critical to ENSO's dynamics. The input and output of these processes are identified and analyzed, some of which are motivated by the recharge oscillator theory. Several errors and compensating errors are identified. The east-west slope of the equatorial thermocline is found to respond to the central equatorial Pacific zonal wind stress as a damped driven harmonic oscillator in all models. This result is shown to be inconsistent with two different formulations of the recharge oscillator. East Pacific sea surface temperature (SST) responds consistently to changes in the thermocline depth in the eastern Pacific in the five CMIP models examined here. However, at time scales greater than 2 years, this consistent model response disagrees with observations, showing that the SST leads thermocline depth at long time scales. Compensating errors are present in the response of meridional transport of water away from the equator to SST: two different models show different response of the transport to off-equatorial wind curl and wind curl response to East Pacific SST. However, these two models show the same response of meridional transport to East Pacific SST. Identification of errors in specific physical processes can hopefully lead to model improvement by focusing model development efforts on these processes.

The dynamics of the Atlantic meridional overturning circulation (AMOC) vary considerably among different climate models; for example, some models show clear peaks in their power spectra while others do not. To elucidate these model differences, transfer functions are used to estimate the frequency domain relationship between surface forcing fields, including sea surface temperature, salinity, and wind stress, and the resulting AMOC response. These are estimated from the outputs of the Coupled Model Intercomparison Project phase 5 (CMIP5) and phase 3 (CMIP3) control runs for eight different models, with a specific focus on Geophysical Fluid Dynamics Laboratory Climate Model, version 2.1 (GFDL CM2.1), and the Community Climate System Model, version 4 (CCSM4), which exhibit rather different spectral behavior. The transfer functions show very little agreement among models for any of the pairs of variables considered, suggesting the existence of systematic model errors and that considerable uncertainty in the simulation of AMOC in current climate models remains. However, a robust feature of the frequency domain analysis is that models with spectral peaks in their AMOC correspond to those in which AMOC variability is more strongly excited by high-latitude surface perturbations that have periods corresponding to the frequency of the spectral peaks. This explains why different models exhibit such different AMOC variability. These differences would not be evident without using a method that explicitly computes the frequency dependence rather than a priori assuming a particular functional form. Finally, transfer functions are used to evaluate two proposed physical mechanisms for model differences in AMOC variability: differences in Labrador Sea stratification and excitation by westward-propagating subsurface Rossby waves.

Solar geoengineering has been receiving increased attention in recent years as a potential temporary solution to offset global warming. One method of approximating global-scale solar geoengineering in climate models is via solar reduction experiments. Two generations of models in the Geoengineering Model Intercomparison Project (GeoMIP) have now simulated offsetting a quadrupling of the CO