Search

Recent studies have suggested that injection of solid particles such as alumina and calcite particles for stratospheric aerosol injection (SAI) instead of sulfur-based injections could reduce some of the adverse side effects of SAI such as ozone depletion and stratospheric heating. Here, we present a version of the global aerosol–chemistry–climate model SOCOL-AERv2 and the Earth system model (ESM) SOCOLv4 which incorporate a solid-particle microphysics scheme for assessment of SAI of solid particles. Microphysical interactions of the solid particle with the stratospheric sulfur cycle were interactively coupled to the heterogeneous chemistry scheme and the radiative transfer code (RTC) for the first time within an ESM. Therefore, the model allows simulation of heterogeneous chemistry at the particle surface as well as feedbacks between microphysics, chemistry, radiation and climate. We show that sulfur-based SAI results in a doubling of the stratospheric aerosol burden compared to the same mass injection rate of calcite and alumina particles with a radius of 240 nm. Most of the sulfuric acid aerosol mass resulting from SO2 injection does not need to be lifted to the stratosphere but is formed after in situ oxidation and subsequent water uptake in the stratosphere. Therefore, to achieve the same radiative forcing, larger injection rates are needed for calcite and alumina particle injection than for sulfur-based SAI. The stratospheric sulfur cycle would be significantly perturbed, with a reduction in stratospheric sulfuric acid burden by 53 %, when injecting 5 Mtyr-1 (megatons per year) of alumina or calcite particles of 240 nm radius. We show that alumina particles will acquire a sulfuric acid coating equivalent to about 10 nm thickness if the sulfuric acid is equally distributed over the whole available particle surface area in the lower stratosphere. However, due to the steep contact angle of sulfuric acid on alumina particles, the sulfuric acid coating would likely not cover the entire alumina surface, which would result in available surface for heterogeneous reactions other than the ones on sulfuric acid. When applying realistic uptake coefficients of 1.0, 10-5 and 10-4 for H2SO4 , HCl and HNO3 , respectively, the same scenario with injections of calcite particles results in 94 % of the particle mass remaining in the form of CaCO3. This likely keeps the optical properties of the calcite particles intact but could significantly alter the heterogeneous reactions occurring on the particle surfaces. The major process uncertainties of solid-particle SAI are (1) the solid-particle microphysics in the injection plume and degree of agglomeration of solid particles on the sub-ESM grid scale, (2) the scattering properties of the resulting agglomerates, (3) heterogeneous chemistry on the particle surface, and (4) aerosol–cloud interactions. These uncertainties can only be addressed with extensive, coordinated experimental and modelling research efforts. The model presented in this work offers a useful tool for sensitivity studies and incorporating new experimental results on SAI of solid particles. [ABSTRACT FROM AUTHOR]

An accurate estimate of global hydroxyl radical (OH) abundance is important for projections of air quality, climate, and stratospheric ozone recovery. As the atmospheric mixing ratios of methyl chloroform (CH3CCl3) (MCF), the commonly used OH reference gas, approaches zero, it is important to find alternative approaches to infer atmospheric OH abundance and variability. The lack of global bottom-up emission inventories is the primary obstacle in choosing a MCF alternative. We illustrate that global emissions of long-lived trace gases can be inferred from their observed mixing ratio differences between the Northern Hemisphere (NH) and Southern Hemisphere (SH), given realistic estimates of their NH-SH exchange time, the emission partitioning between the two hemispheres, and the NH versus SH OH abundance ratio. Using the observed long-term trend and emissions derived from the measured hemispheric gradient, the combination of HFC-32 (CH2F2), HFC-134a (CH2FCF3, HFC-152a (CH3CHF2), and HCFC-22 (CHClF2), instead of a single gas, will be useful as a MCF alternative to infer global and hemispheric OH abundance and trace gas lifetimes. The primary assumption on which this multispecies approach relies is that the OH lifetimes can be estimated by scaling the thermal reaction rates of a reference gas at 272 K on global and hemispheric scales. Thus, the derived hemispheric and global OH estimates are forced to reconcile the observed trends and gradient for all four compounds simultaneously. However, currently, observations of these gases from the surface networks do not provide more accurate OH abundance estimate than that from MCF.

Reductions in sulfur (S) atmospheric deposition in recent decades have been attributed to S deficiencies in crops. Similarly, global soil selenium (Se) concentrations were predicted to drop, particularly in Europe, due to increases in leaching attributed to increases in aridity. Given its international importance in agriculture, reductions of essential elements, including S and Se, in European soils could have important impacts on nutrition and human health. Our objectives were to model current soil S and Se levels in Europe and predict concentration changes for the 21st century. We interrogated four machine-learning (ML) techniques, but after critical evaluation, only outputs for linear support vector regression (Lin-SVR) models for S and Se and the multilayer perceptron model (MLP) for Se were consistent with known mechanisms reported in literature. Other models exhibited overfitting even when differences in training and testing performance were low or non-existent. Furthermore, our results highlight that similarly performing models based on RMSE or R2 can lead to drastically different predictions and conclusions, thus highlighting the need to interrogate machine learning models and to ensure they are consistent with known mechanisms reported in the literature. Both elements exhibited similar spatial patterns with predicted gains in Scandinavia versus losses in the central and Mediterranean regions of Europe, respectively, by the end of the 21st century for an extreme climate scenario. The median change was −5.5% for S (Lin-SVR) and −3.5% (MLP) and −4.0% (Lin-SVR) for Se. For both elements, modeled losses were driven by decreases in soil organic carbon, S and Se atmospheric deposition, and gains were driven by increases in evapotranspiration. [ABSTRACT FROM AUTHOR]

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 SO 2 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 Wm-2 , depending on the microphysical configuration. Nucleation first shifts the size distribution towards radii better suited for solar scattering (0.3 µm

Stratospheric aerosols play important roles in Earth's radiative budget and in heterogeneous chemistry. Volcanic eruptions modulate the stratospheric aerosol layer by injecting particles and particle precursors like sulfur dioxide (SO 2) into the stratosphere. Beginning on 9 April 2021, La Soufrière erupted, injecting SO 2 into the tropical upper troposphere and lower stratosphere, yielding a peak SO 2 loading of 0.3–0.4 Tg. The resulting volcanic aerosol plumes dispersed predominately over the Northern Hemisphere (NH), as indicated by the CALIOP/CALIPSO satellite observations and model simulations. From June to August 2021 and May to July 2022, the NASA ER-2 high-altitude aircraft extensively sampled the stratospheric aerosol layer over the continental United States during the Dynamics and Chemistry of the Summer Stratosphere (DCOTSS) mission. These in situ aerosol measurements provide detailed insights into the number concentration, size distribution, and spatiotemporal variations of particles within volcanic plumes. Notably, aerosol surface area density and number density in 2021 were enhanced by a factor of 2–4 between 380–500 K potential temperature compared to the 2022 DCOTSS observations, which were minimally influenced by volcanic activity. Within the volcanic plume, the observed aerosol number density exhibited significant meridional and zonal variations, while the mode and shape of aerosol size distributions did not vary. The La Soufrière eruption led to an increase in the number concentration of small particles (<400 nm), resulting in a smaller aerosol effective diameter during the summer of 2021 compared to the baseline conditions in the summer of 2022, as observed in regular ER-2 profiles over Salina, Kansas. A similar reduction in aerosol effective diameter was not observed in ER-2 profiles over Palmdale, California, possibly due to the values that were already smaller in that region during the limited sampling period in 2022. Additionally, we modeled the eruption with the SOCOL-AERv2 aerosol–chemistry–climate model. The modeled aerosol enhancement aligned well with DCOTSS observations, although the direct comparison was complicated by issues related to the model's background aerosol burden. This study indicates that the La Soufrière eruption contributed approximately 0.6 % to Arctic and Antarctic ozone column depletion in both 2021 and 2022, which is well within the range of natural variability. The modeled top-of-atmosphere 1-year global average radiative forcing was -0.08 W m -2 clear-sky and -0.04 W m -2 all-sky. The radiative effects were concentrated in the tropics and NH midlatitudes and diminished to near-baseline levels after 1 year. [ABSTRACT FROM AUTHOR]

Stratospheric ozone is expected to recover by the mid-century due to the success of the Montreal Protocol in regulating the emission of ozone-depleting substances (ODSs). In the Arctic, ozone abundances are projected to surpass historical levels due to the combined effect of decreasing ODSs and elevated greenhouse gases (GHGs). While long-term changes in stratospheric ozone have been shown to be a major driver of future surface climate in the Southern Hemisphere during summertime, the dynamical and climatic impacts of elevated ozone levels in the Arctic have not been investigated. In this study, we use two chemistry climate models (the SOlar Climate Ozone Links – Max Planck Ocean Model (SOCOL-MPIOM) and the Community Earth System Model – Whole Atmosphere Community Climate Model (CESM-WACCM)) to assess the climatic impacts of future changes in Arctic ozone on stratospheric dynamics and surface climate in the Northern Hemisphere (NH) during the 21st century. Under the high-emission scenario (RCP8.5) examined in this work, Arctic ozone returns to pre-industrial levels by the middle of the century. Thereby, the increase in Arctic ozone in this scenario warms the lower Arctic stratosphere; reduces the strength of the polar vortex, advancing its breakdown; and weakens the Brewer–Dobson circulation. The ozone-induced changes in springtime generally oppose the effects of GHGs on the polar vortex. In the troposphere, future changes in Arctic ozone induce a negative phase of the Arctic Oscillation, pushing the jet equatorward over the North Atlantic. These impacts of future ozone changes on NH surface climate are smaller than the effects of GHGs, but they are remarkably robust among the two models employed in this study, canceling out a portion of the GHG effects (up to 20 % over the Arctic). In the stratosphere, Arctic ozone changes cancel out a much larger fraction of the GHG-induced signal (up to 50 %–100 %), resulting in no overall change in the projected springtime stratospheric northern annular mode and a reduction in the GHG-induced delay of vortex breakdown of around 15 d. Taken together, our results indicate that future changes in Arctic ozone actively shape the projected changes in the stratospheric circulation and their coupling to the troposphere, thereby playing an important and previously unrecognized role as a driver of the large-scale atmospheric circulation response to climate change. [ABSTRACT FROM AUTHOR]

Solar radiation management as 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 side effects of hypothesised strat-SRM scenarios. A recent strat-SRM model intercompar-ison study for composition-climate models with interactive stratospheric aerosol suggests that the modelled climate response 5 to a particular assumed injection strategy, depends on the type of aerosol microphysical scheme used (e.g., modal or sectional representation), alongside also host model resolution and transport. Compared to short-duration volcanic SO2 emission, the continuous SO2 injections in strat-SRM scenarios may pose a greater challenge to the numerical implementation of of microphysical processes such as nucleation, condensation, and coagulation. This study explores how changing the timesteps and sequencing of microphysical processes in the sectional aerosol-chemistry-climate model SOCOL-AERv2 (40 size bins) 10 affect model predicted climate and ozone layer impacts considering strat-SRM SO2 injections of of 5 and 25 Tg(S) yr-1 at 20 km altitude between 30°S and 30°N. The model experiments consider year 2040 boundary conditions for ozone depleting substances and green house gases. We focus on the length of the microphysical timestep 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 15 scenarios injecting 25 Mt/yr of sulfur with a default microphysical timesetp 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 nucle-ation 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 timestep 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 20 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 W m-2 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 ultra-fine aerosol particles, increasing the surface area density, resulting in 10 DU less ozone (about 3% of total column) in 25 the northern midlatitudes 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 minutes or less must be applied to accurately capture the magnitude of the H2SO2 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. [ABSTRACT FROM AUTHOR]

Stratospheric aerosols play important roles in Earth's radiative budget and in heterogeneous chemistry. Volcanic eruptions modulate the stratospheric aerosol layer by injecting particles and particle precursors like sulfur dioxide (SO2) into the stratosphere. Beginning on April 9th, 2021, La Soufrière erupted injecting SO2 into the tropical upper troposphere and lower stratosphere, yielding a peak SO2 loading of 0.3-0.4 Tg. The resulting volcanic aerosol plumes dispersed predominately over the northern hemisphere (NH), as indicated by the CALIOP/CALIPSO satellite observations and model simulations. From June to August 2021 and May to July 2022, the NASA ER-2 high-altitude aircraft extensively sampled the stratospheric aerosol layer over the continental United States during the Dynamics and Chemistry of the Summer Stratosphere (DCOTSS) mission. These in situ aerosol measurements provide detailed insights into the number concentration, size distribution, and spatiotemporal variations of particles within volcanic plumes. Notably, aerosol surface area density and number density in 2021 were enhanced by a factor of 2-4 between 380-500 K potential temperature compared to the 2022 DCOTSS observations, which were minimally influenced by volcanic activity. Within the volcanic plume, the observed aerosol number density exhibited significant meridional and zonal variations while the mode and shape of aerosol size distributions did not vary. The La Soufrière eruption led to an increase in the number concentration of small particles (<400 nm), resulting in a smaller aerosol effective diameter during the summer of 2021 compared to the baseline conditions in the summer of 2022, as observed in regular ER-2 profiles over Salina, Kansas. A similar reduction in aerosol effective diameter was not observed in ER-2 profiles over Palmdale, California, possibly due to the already smaller values in that region during the limited sampling period in 2022. The La Soufrière eruption was modeled with the SOCOL-AERv2 aerosol-chemistry-climate model. The modeled aerosol enhancement aligned well with DCOTSS observations, although the direct comparison was complicated by issues related to the model's background aerosol burden. This study indicates that the La Soufrière eruption contributed at most 0.6% to Arctic and Antarctic ozone column depletion in both 2021 and 2022, which is well within the range of natural variability. The modeled top-of-atmosphere one-year global average radiative forcing was -0.08 W/m2 clear-sky and -0.04 W/m2 all-sky. The radiative effects were concentrated in the tropics and NH midlatitudes and diminished to near-baseline levels after one year. [ABSTRACT FROM AUTHOR]

Stratospheric ozone is expected to recover by mid-century due to the success of the Montreal Protocol in regulating the emission of ozone-depleting substances (ODSs). In the Arctic, ozone abundances are projected to surpass historical levels due to the combined effect of decreasing ODSs and elevated greenhouse gases (GHGs). While ozone recovery has been shown to be a major driver of future surface climate in the Southern Hemisphere during summertime, the dynamical and climatic impacts of elevated ozone levels in the Arctic have not been investigated. In this study, we use two chemistry climate models (SOCOL-MPIOM and CESM-WACCM) to assess the climatic impacts of Arctic ozone recovery on stratospheric dynamics and surface climate in the Northern Hemisphere (NH) during the 21st century. Under the high-emission scenario (RCP8.5) examined in this work, Arctic ozone returns to pre-industrial levels by the middle of the century. Thereby, it warms the lower Arctic stratosphere, reduces the strength of the polar vortex, advancing its breakdown, and weakening the Brewer-Dobson circulation. In the troposphere, Arctic ozone recovery induces a negative phase of the Arctic Oscillation, pushing the jet equatorward over the Atlantic. These impacts of ozone recovery in the NH are smaller than the effects of GHGs, but they are remarkably robust among the two models employed in this study, cancelling out some of the GHG effects. Taken together, our results indicate that Arctic ozone recovery actively shapes the projected changes in the stratospheric circulation and their coupling to the troposphere, thereby playing an important and previously unrecognized role as driver of the large-scale atmospheric circulation response to climate change [ABSTRACT FROM AUTHOR]