SO2 Injection Research Articles

Observing the SO2 and Sulfate Aerosol Plumes From the 2022 Hunga Eruption With the Infrared Atmospheric Sounding Interferometer (IASI)

AbstractThe Hunga volcano violently erupted on 15 January 2022, producing the largest perturbation of the stratospheric aerosol layer since Pinatubo 1991, despite the initially estimated modest injection of SO2. This study presents novel SO2 and sulfate aerosol (SA) co‐retrievals from the Infrared Atmospheric Sounding Interferometer, and uses them to quantify the initial progression of the Hunga plume. These observations are consistent with rapid conversion of SO2 (e‐folding time: 17.1 ± 4.3 days) to SA, with an injected burden of >1.0 Tg SO2. This points at larger SO2 injections than previously thought. A long‐lasting SA plume was observed, with two separate build‐up phases, and with a meridional dispersion of marked anomalies from the tropics to the higher southern hemispheric latitudes. A limited (∼20%) SA removal was observed after 1‐year dispersion. The total injected SA mass burden was estimated at 1.6 ± 0.5 Tg in the total atmospheric column, with a build‐up e‐folding time of about 2 months.

Microphysical Simulation of the 2022 Hunga Volcano Eruption Using a Sectional Aerosol Model

AbstractApproximately 150 Tg of water vapor and 0.42 Tg of sulfur dioxide were injected directly into the stratosphere by the January 2022 Hunga volcanic eruption, which represents the largest water vapor injection in the satellite era. A comparison of numerical simulations to balloon‐borne and satellite observations of the water‐rich plume suggests that particle coagulation contributed to the Hunga aerosol's effective dry radius increase from 0.2 μm in February to around 0.4 μm in March. Our model suggests that the stratospheric aerosol effective radius is persistently perturbed for years by moderate and large‐magnitude volcanic events, whereas extreme wildfire events show limited impact on the stratospheric background particle size. Our analysis further suggests that both the particle optical efficiency and the aerosols' stratospheric lifetime explain Hunga's unusually large aerosol optical depth per unit of the SO2 injection, as compared with the Pinatubo eruption.

Carbon Cycle Response to Stratospheric Aerosol Injection With Multiple Temperature Stabilization Targets and Strategies

AbstractWe analyze the global carbon cycle response to a set of stratospheric aerosol injection (SAI) simulations performed by the CESM2(WACCM6‐MA) model. The simulations are performed under the specified SSP2‐4.5 CO2 concentration pathway. It is found that both the temperature stabilization target and the SO2 injection strategy have important effects on the global carbon sink. Relative to the SSP2‐4.5 scenario, averaged over the last 20 years of our simulations (year 2050–2069), simultaneous multi‐location SO2 injection causes an increase in cumulative land carbon uptake of 45 and 23 PgC, and an increase in cumulative ocean carbon uptake of 6 and 2 PgC for temperature stabilization targets of 0.5°C and 1.5°C respectively. For a temperature stabilization target of 1.0°C, SO2 injections increase land and ocean carbon sinks by 22–42 PgC and 4–7 PgC, respectively, depending on the strategies of SO2 injections (low latitude, mid‐to‐high latitude, and multi‐objective injection). Relative to SSP2‐4.5, by year 2069, SAI increases diagnosed cumulative CO2 emissions by 25–53 PgC (3%–6%), implying a decrease in atmospheric CO2 if SO2 injections were performed under a prescribed CO2 emission pathway. Stratospheric SO2 injections slow permafrost thaw, but do not restore permafrost to the previous extent at the same warming level for all injection strategies. An abrupt termination of SO2 injection weakens both the ocean and land carbon sink, and causes a rapid decline of permafrost extent. A gradual phaseout of SO2 injection slows sharp decline of permafrost and delays the rebound of carbon sink.

Importance of microphysical settings for climate forcing by stratospheric SO2 injections as modeled by SOCOL-AERv2

Abstract. Solar radiation modification by a sustained deliberate source of SO2 into the stratosphere (strat-SRM) has been proposed as an option for climate intervention. Global interactive aerosol–chemistry–climate models are often used to investigate the potential cooling efficiencies and associated side effects of hypothesized strat-SRM scenarios. A recent model intercomparison study for composition–climate models with interactive stratospheric aerosol suggests that the modeled climate response to a particular assumed injection strategy depends on the type of aerosol microphysical scheme used (e.g., modal or sectional representation) alongside host model resolution and transport. Compared to short-duration volcanic SO2 emissions, the continuous SO2 injections in strat-SRM scenarios may pose a greater challenge to the numerical implementation of microphysical processes such as nucleation, condensation, and coagulation. This study explores how changing the time steps and sequencing of microphysical processes in the sectional aerosol–chemistry–climate model SOCOL-AERv2 (40 mass bins) affects model-predicted climate and ozone layer impacts considering strat-SRM by SO2 injections of 5 and 25 Tg(S) yr−1 at 20 km altitude between 30° S and 30° N. The model experiments consider the year 2040 to be the boundary conditions for ozone-depleting substances and greenhouse gases (GHGs). We focus on the length of the microphysical time step and the call sequence of nucleation and condensation, the two competing sink processes for gaseous H2SO4. Under stratospheric background conditions, we find no effect of the microphysical setup on the simulated aerosol properties. However, at the high sulfur loadings reached in the scenarios injecting 25 Tg(S) yr−1 of SO2 with a default microphysical time step of 6 min, changing the call sequence from the default “condensation first” to “nucleation first” leads to a massive increase in the number densities of particles in the nucleation mode (R<0.01 µm) and a small decrease in coarse-mode particles (R>1 µm). As expected, the influence of the call sequence becomes negligible when the microphysical time step is reduced to a few seconds, with the model solutions converging to a size distribution with a pronounced nucleation mode. While the main features and spatial patterns of climate forcing by SO2 injections are not strongly affected by the microphysical configuration, the absolute numbers vary considerably. For the extreme injection with 25 Tg(S) yr−1, the simulated net global radiative forcing ranges from −2.3 to −5.3 W m−2, depending on the microphysical configuration. Nucleation first shifts the size distribution towards radii better suited for solar scattering (0.3 µm <R< 0.4 µm), enhancing the intervention efficiency. The size distribution shift, however, generates more ultrafine aerosol particles, increasing the surface area density and resulting in 10 DU (Dobson units) less ozone (about 3 % of the total column) in the northern mid-latitudes and 20 DU less ozone (6 %) over the polar caps compared to the condensation first approach. Our results suggest that a reasonably short microphysical time step of 2 min or less must be applied to accurately capture the magnitude of the H2SO4 supersaturation resulting from SO2 injection scenarios or volcanic eruptions. Taken together, these results underscore how structural aspects of model representation of aerosol microphysical processes become important under conditions of elevated stratospheric sulfur in determining atmospheric chemistry and climate impacts.

The 2019 Raikoke eruption as a testbed used by the Volcano Response group for rapid assessment of volcanic atmospheric impacts

Abstract. ​​​​​​​The 21 June 2019 Raikoke eruption (48° N, 153° E) generated one of the largest amounts of sulfur emission to the stratosphere since the 1991 Mt. Pinatubo eruption. Satellite measurements indicate a consensus best estimate of 1.5 Tg for the sulfur dioxide (SO2) injected at an altitude of around 14–15 km. The peak Northern Hemisphere (NH) mean 525 nm stratospheric aerosol optical depth (SAOD) increased to 0.025, a factor of 3 higher than background levels. The Volcano Response (VolRes) initiative provided a platform for the community to share information about this eruption which significantly enhanced coordination efforts in the days after the eruption. A multi-platform satellite observation subgroup formed to prepare an initial report to present eruption parameters including SO2 emissions and their vertical distribution for the modeling community. It allowed us to make the first estimate of what would be the peak in SAOD 1 week after the eruption using a simple volcanic aerosol model. In this retrospective analysis, we show that revised volcanic SO2 injection profiles yield a higher peak injection of the SO2 mass. This highlights difficulties in accurately representing the vertical distribution for moderate SO2 explosive eruptions in the lowermost stratosphere due to limited vertical sensitivity of the current satellite sensors (±2 km accuracy) and low horizontal resolution of lidar observations. We also show that the SO2 lifetime initially assumed in the simple aerosol model was overestimated by 66 %, pointing to challenges for simple models to capture how the life cycle of volcanic gases and aerosols depends on the SO2 injection magnitude, latitude, and height. Using a revised injection profile, modeling results indicate a peak NH monthly mean SAOD at 525 nm of 0.024, in excellent agreement with observations, associated with a global monthly mean radiative forcing of −0.17 W m−2 resulting in an annual global mean surface temperature anomaly of −0.028 K. Given the relatively small magnitude of the forcing, it is unlikely that the surface response can be dissociated from surface temperature variability.

Formation integrity evaluation for geosequestration of CO2 in depleted petroleum reservoirs under cyclic stress conditions

The geological storage of carbon dioxide (CO2), also referred to as CO2 geosequestration, represents one of the most promising options for reducing greenhouse gases in the atmosphere. However, most of the time, CO2 is captured with small amounts of other industrial gases such as sulphur dioxide (SO2) and hydrogen sulphide (H2S), which might be compressed together and stored in depleted petroleum reservoirs or aquifers. Moreover, during CO2 geosequestration in reservoirs, pressure variations during injection could force some amount of CO2 (with or without other acid gas impurities) into the caprock; thereby, altering the petrophysical, geochemical, and geomechanical properties of the caprock. Thus, the brittleness index of the reservoir and caprock might be impacted during CO2 geosequestration due to the chemical reactions between the rock minerals and the formation fluid. Furthermore, to meet the world net-zero carbon target, the promotion of CO2 utilization is paramount. This could be possible by developing an effective technology for cyclic CO2 geosequestration (with or without gas impurities). Therefore, studies on the co-injection of CO2 with other acid gases from industrial emissions, their withdrawal from the porous medium, and their impact on reservoir and caprock integrity are paramount. In this study, a dual-tubing string well completion technology was designed for cyclic injection and withdrawal of CO2 (with or without another acid gas), and numerical simulations were performed using TOUGHREACT codes, to model the cyclic process and investigate the co-injection of SO2 and H2S (separately) with CO2 in sandstone formations overlain by shale caprock. A novel technique of converting the volume fraction of minerals to their weight fraction was developed in this study, to evaluate the brittleness index of the sandstone reservoir and shale caprock during CO2 geosequestration. The findings of the study indicate that the porosity and permeability increase for the CO2 only and CO2–H2S injection cases, in the shale caprock; while for the CO2–SO2 injection case, porosity and permeability only decreased in the layers of the shale caprock contacted by SO2 and due to anhydrite precipitation. In all the injection cases, the porosity and permeability of the sandstone reservoir decreased in a few layers directly below the perforation interval of the production zone. However, in other regions in the sandstone reservoir, the porosity and permeability increased for the CO2 only and CO2–H2S injection cases. In contrast, for the CO2–SO2 co-injection case, porosity and permeability decreased in the layers of the sandstone rock contacted by SO2. In all the CO2 geosequestration cases, the brittleness of the shale and sandstone rocks investigated decreased slightly, except in the CO2–SO2 co-injection case where the brittleness of the sandstone rock decreased significantly. Based on the mineralogical composition of the formations in this study, co-injection of SO2 gas with CO2 gas, only decreased the brittleness index of the shale caprock slightly, but significantly decreased the brittleness of the sandstone reservoir.

Mechanisms of CO2-rich industrial waste gas enhanced shale oil recovery in kerogen slit based on adsorption behavior, gas flooding and surfactant synergy study

Enhanced oil recovery (EOR) by injection of CO2-rich industrial waste gas has been proposed as an energy-efficient and environmentally friendly method for shale oil reservoirs, in which the waste gas can be utilized and sequestered to reduce greenhouse gas emissions. In this work, the adsorption behaviors of shale oil and waste gas components in kerogen slits and the effects of gas flooding on EOR were investigated using molecular dynamics (MD) method, which can provide technical guide for shale oil development. In terms of the adsorption behaviors, an intensity decrease of oil adsorption is observed in kerogen slits after the injection of the waste gas components. Desorption of shale oil during injection of SO2 and CO2 is more pronounced in shallow reservoirs, whereas it is more evident for NO and N2 in deep reservoirs. The desorption efficiency can be 10.91 % for pure CO2 injection and 9.09 % for NO injection at a depth of 1000 m, and the relative stability ratios of CO2 and NO sequestration are 77.09 % and 76.06 %, respectively. For gas flooding, the existence of pressure difference provides favorable conditions for gas molecules to pass through the oil zone, and waste gas can be more effective than pure CO2 as driving gas. For the waste gas injection, the percentage of CO2 retained in the kerogen nanopores is much higher than that of pure CO2 injection, suggesting that more CO2 molecules can enter into kerogen matrix to take part in the displacement process. Furthermore, the effect of surfactants on oil desorption and displacement was also investigated and the interaction energy of the surfactant/oil/gas system was calculated and analyzed, finding that addition of the oil-soluble surfactants can be beneficial to shale oil displacement. This work can provide a comprehensive evaluation of the feasibility of CO2-rich industrial waste gas application for EOR, which is of great significance in the effective production of shale oil and the sequestration of waste gas in shale reservoirs.

How well are aerosol–cloud interactions represented in climate models? – Part 1: Understanding the sulfate aerosol production from the 2014–15 Holuhraun eruption

Abstract. For over 6 months, the 2014–2015 effusive eruption at Holuhraun, Iceland, injected considerable amounts of sulfur dioxide (SO2) into the lower troposphere with a daily rate of up to one-third of the global emission rate, causing extensive air pollution across Europe. The large injection of SO2, which oxidises to form sulfate aerosol (SO42-), provides a natural experiment offering an ideal opportunity to scrutinise state-of-the-art general circulation models' (GCMs) representation of aerosol–cloud interactions (ACIs). Here we present Part 1 of a two-part model inter-comparison using the Holuhraun eruption as a framework to analyse ACIs. We use SO2 retrievals from the Infrared Atmospheric Sounding Interferometer (IASI) instrument and ground-based measurements of SO2 and SO42- mass concentrations across Europe, in conjunction with a trajectory analysis using the Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) model, to assess the spatial and chemical evolution of the volcanic plume as simulated by five GCMs and a chemical transport model (CTM). IASI retrievals of plume altitude and SO2 column load reveal that the volcanic perturbation is largely contained within the lower troposphere. Compared to the satellite observations, the models capture the spatial evolution and vertical variability of the plume reasonably well, although the models often overestimate the plume altitude. HYSPLIT trajectories are used to attribute to Holuhraun emissions 111 instances of elevated sulfurous surface mass concentrations recorded at European Monitoring and Evaluation Programme (EMEP) stations during September and October 2014. Comparisons with the simulated concentrations show that the modelled ratio of SO2 to SO42- during these pollution episodes is often underestimated and overestimated for the young and mature plume, respectively. Models with finer vertical resolutions near the surface are found to better capture these elevated sulfurous ground-level concentrations. Using an exponential function to describe the decay of observed surface mass concentration ratios of SO2 to SO42- with plume age, the in-plume oxidation rate constant is estimated as 0.032 ± 0.002 h−1 (1.30 ± 0.08 d e-folding time), with a near-vent ratio of 25 ± 5 (µg m−3 of SO2 / µg m−3 of SO42-). The majority of the corresponding derived modelled oxidation rate constants are lower than the observed estimate. This suggests that the representation of the oxidation pathway/s in the simulated plumes is too slow. Overall, despite their coarse spatial resolutions, the six models show reasonable skill in capturing the spatial and chemical evolution of the Holuhraun plume. This capable representation of the underlying aerosol perturbation is essential to enable the investigation of the eruption's impact on ACIs in the second part of this study.

The Emissions Model Intercomparison Project (Emissions-MIP): quantifying model sensitivity to emission characteristics

Abstract. Anthropogenic emissions of aerosols and precursor compounds are known to significantly affect the energy balance of the Earth–atmosphere system, alter the formation of clouds and precipitation, and have a substantial impact on human health and the environment. Global models are an essential tool for examining the impacts of these emissions. In this study, we examine the sensitivity of model results to the assumed height of SO2 injection, seasonality of SO2 and black carbon (BC) particulate emissions, and the assumed fraction of SO2 emissions that is injected into the atmosphere as particulate phase sulfate (SO4) in 11 climate and chemistry models, including both chemical transport models and the atmospheric component of Earth system models. We find large variation in atmospheric lifetime across models for SO2, SO4, and BC, with a particularly large relative variation for SO2, which indicates that fundamental aspects of atmospheric sulfur chemistry remain uncertain. Of the perturbations examined in this study, the assumed height of SO2 injection had the largest overall impacts, particularly on global mean net radiative flux (maximum difference of −0.35 W m−2), SO2 lifetime over Northern Hemisphere land (maximum difference of 0.8 d), surface SO2 concentration (up to 59 % decrease), and surface sulfate concentration (up to 23 % increase). Emitting SO2 at height consistently increased SO2 and SO4 column burdens and shortwave cooling, with varying magnitudes, but had inconsistent effects across models on the sign of the change in implied cloud forcing. The assumed SO4 emission fraction also had a significant impact on net radiative flux and surface sulfate concentration. Because these properties are not standardized across models this is a source of inter-model diversity typically neglected in model intercomparisons. These results imply a need to ensure that anthropogenic emission injection height and SO4 emission fraction are accurately and consistently represented in global models.

Tonga volcanic eruption triggered anomalous Arctic warming in early 2022

On 15 January 2022, the submarine volcano on Hunga Tonga-Hunga Ha'apai Island erupted and injected a small amount of SO2 but an unprecedent amount of water vapor into the stratosphere. The impact of Tonga eruption on climate was dynamically simulated using a fully coupled climate prediction system FIO-CPS v2.0 by considering both SO2 and water vapor injections into the stratosphere, which is the first attempt to incorporate the volcanic water vapor injection in a climate prediction system to evaluate its effect on surface climate response. Model results indicated that the surface air temperature (SAT) over the Arctic region was significantly warmer in early 2022 as a response to the co-injection of SO2 and water vapor from Tonga eruption, although its effect on global mean SAT was found to be insignificant. In February 2022, the average warming anomaly in the Arctic area reached up to 1.1 °C, with zonal mean warming of 2.0 °C near the North Pole. The Tonga eruptions induced anomalous heat transport which primarily contributed to the Arctic warming in February. However, it had a cooling effect on the Arctic area in January and March. In addition to heat transport, the eruption also triggered anomalous northward transport of water vapor. This led to an increase in the water vapor content, which further enhanced the downward longwave radiation and cloud-radiative feedback over the Arctic region, ultimately contributing to the Arctic SAT warming in early 2022. In the simulation only SO2 injection was considered without water vapor injection, the pattern of SAT response showed warming over the Eurasian continent and cooling in the north over Arctic Ocean, which is opposite to that in the simulation with co-injection of SO2 and water vapor. This indicated that volcanic water vapor injection of Tonga eruption played a dominant role in the Arctic warming. Consequently, submarine volcano should receive more attention due to its potential climate effects, especially when it involves significant water vapor releases.

Potential impact of stratospheric aerosol geoengineering on projected temperature and precipitation extremes in South Africa

Stratospheric aerosol injection (SAI) is the theoretical deployment of sulphate particles into the stratosphere to reflect incoming solar radiation and trigger a cooling impact at the Earth’s surface. This study assessed the potential impact of SAI geoengineering on temperature and precipitation extremes over South Africa (SAF) and its climatic zones in the future (2075–2095) using simulations from the Stratospheric Aerosol Geoengineering Large Ensemble (GLENS) project. We analyse three different experiments from the GLENS project, each of which simulate stratospheric SO2 injection under the representative concentration pathway 8.5 (RCP8.5) emissions scenario: (i) tropical injection around 22.8–25 km altitude (GLENS), (ii) tropical injection around 1 km above the tropopause (GLENS_low), and (iii) injection near the equator around 20–25 km (GLENS_eq). The study used a set of the Expert Team on Climate Change Detection and Indices describing temperature and rainfall extremes to assess the impact of the three SAI experiments on extreme weather in the future over SAF. The results of this study indicate that, relative to the baseline period (2010–2030), all three SAI experiments are mostly over-effective in offsetting the projected RCP8.5 increase in the frequency of hot (up to −60%) and decrease (up to +10%) in cold temperature extremes over SAF and its climatic zones. These findings suggest that SAI could cause over-cooling in SAF. However, SAI impact on precipitation extremes is less linear and varies across the country’s climatic zones. For example, SAI could reinforce the projected decrease in precipitation extremes across most of SAF, although it could exacerbate heavy precipitation over the KwaZulu-Natal Coast. These findings are consistent across SAI experiments except in magnitude, as GLENS_eq and GLENS_low could cause larger decreases in precipitation extremes than GLENS. These findings imply that SAI could alleviate heat stress on human health, agriculture, and vulnerable communities while simultaneously decreasing infrastructure and crops’ vulnerability to flooding. It is, however, essential to interpret these findings cautiously as they are specific to the SAI experiments and modelling settings considered in the GLENS project.

Mechanistic and performance insights into low-temperature NH3-SCR based on Ce-modified Mn-Ti catalysts

In this study, Mn-Ti supported Ce catalysts (derived from manganese tailings) were synthesized by coprecipitation method. The denitration performance and SO2 resistance of the catalysts were investigated. The results showed that introducing Ce oxides into Mn-Ti catalyst improved the low temperature catalytic performance, with almost 100% NO conversion at 225 °C. Additionally, the N2 selectivity was improved due to the decrease in lattice oxygen. The modified catalyst also presented excellent SO2 resistance, with around 80% NO conversion after 10 h of SO2 injection. The modified catalyst had higher ratios of Mn4+ and Mn3+ species, stronger surface acidity, and more redox properties than the Mn-Ti catalyst, which improved NH3-SCR performance. After SO2 poisoning, the surface Lewis acidity decreased and the surface Brønsted acidity became stronger, indicating that the loss of Lewis acidity might cause the decrease of NO conversion. Moreover, the nitrate species on Mn-Ti catalyst were less than those on Mn-Ce-Ti catalyst in the presence of SO2, suggesting that doping Ce into Mn-Ti catalyst retained more oxidation performance after being poisoned by SO2.

Volcanic Drivers of Stratospheric Sulfur in GFDL ESM4

AbstractStratospheric injections of sulfur dioxide from major volcanic eruptions perturb the Earth's global radiative balance and dominate variability in stratospheric sulfur loading. The atmospheric component of the GFDL Earth System Model (ESM4.1) uses a bulk aerosol scheme and previously prescribed the distribution of aerosol optical properties in the stratosphere. To quantify volcanic contributions to the stratospheric sulfur cycle and the resulting climate impact, we modified ESM4.1 to simulate stratospheric sulfate aerosols prognostically. Driven by explicit volcanic emissions of aerosol precursors and non‐volcanic sources, we conduct ESM4.1 simulations from 1989 to 2014, with a focus on the Mt. Pinatubo eruption. We evaluate our interactive representation of the stratospheric sulfur cycle against data from Moderate Resolution Imaging Spectroradiometer, Multi‐angle Imaging SpectroRadiometer, Advanced Very High Resolution Radiometer, High Resolution Infrared Radiation Sounder, and Stratospheric Aerosol and Gas Experiment II. To assess the key processes associated with volcanic aerosols, we performed a sensitivity analysis of sulfate burden from the Mt. Pinatubo eruption by varying injection heights, emission amount, and stratospheric sulfate's dry effective radius. We find that the simulated stratospheric sulfate mass burden and aerosol optical depth in the model are sensitive to these parameters, especially volcanic SO2 injection height, and the optimal combination of parameters depends on the metric we evaluate.

Efficacy and sorption of sulfur dioxide as a fumigant for control of navel orangeworm (Amyelois transitella) on stored pistachios

Sulfur dioxide (SO2) fumigation was evaluated for postharvest control of navel orangeworm (NOW), Amyelois transitella, on stored pistachios in a laboratory study. Sorption of SO2 on pistachios were measured at 5, 15, and 25 °C. Because of rapid sorption of SO2 on pistachios after SO2 injection, 3 h short fumigations were used to determine effective treatments against NOW eggs, larvae, and pupae. SO2 fumigation was effective against all of the life stages tested. However, large larvae in infested pistachios were most tolerant and eggs were most susceptible to SO2 fumigation. Smaller larvae were more susceptible to SO2 fumigation than larger larvae. Complete controls of eggs, larvae, and pupae were achieved in 3 h fumigations with 0.2, 2.0, and 1.0% SO2, respectively. Large SO2 fumigation trials each with 4.5 kg pistachios were conducted in a 19.8L chamber and complete control of 6th instar larvae were achieved in 3 h fumigation with about 1.6–1.8% SO2 at ambient temperatures of 18–21 °C. The results showed that SO2 fumigation has potential to control NOW on stored pistachios as well as to control other postharvest pests on stored products.

Interactive stratospheric aerosol models' response to different amounts and altitudes of SO2 injection during the 1991 Pinatubo eruption

Abstract. A previous model intercomparison of the Tambora aerosol cloud has highlighted substantial differences among simulated volcanic aerosol properties in the pre-industrial stratosphere and has led to questions about the applicability of global aerosol models for large-magnitude explosive eruptions prior to the observational period. Here, we compare the evolution of the stratospheric aerosol cloud following the well-observed June 1991 Mt. Pinatubo eruption simulated with six interactive stratospheric aerosol microphysics models to a range of observational data sets. Our primary focus is on the uncertainties regarding initial SO2 emission following the Pinatubo eruption, as prescribed in the Historical Eruptions SO2 Emission Assessment experiments (HErSEA), in the framework of the Interactive Stratospheric Aerosol Model Intercomparison Project (ISA-MIP). Six global models with interactive aerosol microphysics took part in this study: ECHAM6-SALSA, EMAC, ECHAM5-HAM, SOCOL-AERv2, ULAQ-CCM, and UM-UKCA. Model simulations are performed by varying the SO2 injection amount (ranging between 5 and 10 Tg S) and the altitude of injection (between 18–25 km). The comparisons show that all models consistently demonstrate faster reduction from the peak in sulfate mass burden in the tropical stratosphere. Most models also show a stronger transport towards the extratropics in the Northern Hemisphere, at the expense of the observed tropical confinement, suggesting a much weaker subtropical barrier in all the models, which results in a shorter e-folding time compared to the observations. Furthermore, simulations in which more than 5 Tg S in the form of SO2 is injected show an initial overestimation of the sulfate burden in the tropics and, in some models, in the Northern Hemisphere and a large surface area density a few months after the eruption compared to the values measured in the tropics and the in situ measurements over Laramie. This draws attention to the importance of including processes such as the ash injection for the removal of the initial SO2 and aerosol lofting through local heating.

Climate response to off-equatorial stratospheric sulfur injections in three Earth system models – Part 1: Experimental protocols and surface changes

Abstract. 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.

Investigations on the Global Spread of the Hunga Tonga-Hunga Ha’apai Volcanic Eruption Using Space-Based Observations and Lagrangian Transport Simulations

On 15 January 2022, the Hunga Tonga-Hunga Ha’apai (HTHH) (175.38° W, 20.54° S) volcano erupted explosively. It is considered the most explosive volcanic eruption during the past 140 years. The HTHH volcanic eruption caused intense ripples, Lamb waves, and gravity waves in the atmosphere which encircled the globe several times, as reported by different studies. In this study, using OMI, SAGE-III/ISS, and CALIPSO satellite observations, we investigated the spread of the volcanic SO2 cloud due to the HTHH eruption and subsequent formation of sulfuric acid clouds in the stratosphere. It took about 19–21 days for the stratospheric SO2 injections of the HTHH to encircle the globe longitudinally due to a dominant westward jet with wind speeds of ~2500 km/day, and it slowly dispersed over the whole globe within several months due to poleward spread. The formation of sulfuric acid clouds intensified after about a month, causing the more frequent occurrence of high aerosol optical depth elevated layers in the stratosphere at an altitude of about 20–26 km. This study deals with the dynamics of volcanic plume spread in the stratosphere, knowledge of which is essential in estimating the accurate radiative effects caused by perturbations in the earth–atmosphere system due to a volcanic eruption. In addition, this knowledge provides important input for studies related to the geo-engineering of the earth’s atmosphere by injecting particulates and gases into the stratosphere.

Tracking the 2022 Hunga Tonga-Hunga Ha'apai Aerosol Cloud in the Upper and Middle Stratosphere Using Space-Based Observations.

On 15 January 2022, the submarine Hunga Tonga volcanic eruption lofted materials high into the upper stratosphere, reaching a record-breaking altitude of ∼58km, unprecedented in the satellite observations era. Within two weeks, the bulk of the injected material circulated the globe between 20-30km altitude, as observed by satellite instruments. We estimate that the stratospheric aerosol optical depth (sAOD) is the largest since the Pinatubo eruption and is at least twice as great as the sAOD after the 2015 Calbubo eruption despite the similar SO2 injection from that eruption. We use space-based observations to monitor the Hunga-Tonga volcanic plume evolution and transport at different altitudes as it circulates the globe. While the main aerosol layer remains trapped in the tropical pipe, small parts have already made it to both the northern and southern hemisphere poles by April, which is almost certain to influence this year's ozone hole.

Impact of the Latitude of Stratospheric Aerosol Injection on the Southern Annular Mode

AbstractThe impacts of Stratospheric Aerosol Injection (SAI) strategies on the Southern Annular Mode (SAM) are analyzed with the Community Earth System Model. Using a set of simulations with fixed single‐point SO2 injections we demonstrate the first‐order dependence of the SAM response on the latitude of injection, with the northern hemispheric and equatorial injections driving a response corresponding to a positive phase of SAM and the southern hemispheric injections driving a negative phase of SAM. We further demonstrate that the results can to first order explain the differences in the SAM responses diagnosed from the two recent large ensembles of geoengineering simulations utilizing more complex injection strategies – Geoengineering Large Ensemble and Assessing Responses and Impacts of Solar climate intervention on the Earth system with Stratospheric Aerosol Injection (GLENS and ARISE‐SAI) – as driven by the differences in the simulated sulfate aerosol distributions. Our results point to the meridional extent of aerosol‐induced lower stratospheric heating as an important driver of the sensitivity of the SAM response to the injection location.

The Hunga Tonga-Hunga Ha'apai Hydration of the Stratosphere.

Following the 15 January 2022 Hunga Tonga‐Hunga Ha'apai eruption, several trace gases measured by the Aura Microwave Limb Sounder (MLS) displayed anomalous stratospheric values. Trajectories and radiance simulations confirm that the H2O, SO2, and HCl enhancements were injected by the eruption. In comparison with those from previous eruptions, the SO2 and HCl mass injections were unexceptional, although they reached higher altitudes. In contrast, the H2O injection was unprecedented in both magnitude (far exceeding any previous values in the 17‐year MLS record) and altitude (penetrating into the mesosphere). We estimate the mass of H2O injected into the stratosphere to be 146 ± 5 Tg, or ∼10% of the stratospheric burden. It may take several years for the H2O plume to dissipate. This eruption could impact climate not through surface cooling due to sulfate aerosols, but rather through surface warming due to the radiative forcing from the excess stratospheric H2O.