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]

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]

In the Arctic stratosphere, the combination of chemical ozone depletion by halogenated ozone-depleting substances (hODSs) and dynamic fluctuations can lead to severe ozone minima. These Arctic ozone minima are of great societal concern due to their health and climate impacts. Owing to the success of the Montreal Protocol, hODSs in the stratosphere are gradually declining, resulting in a recovery of the ozone layer. On the other hand, continued greenhouse gas (GHG) emissions cool the stratosphere, possibly enhancing the formation of polar stratospheric clouds (PSCs) and, thus, enabling more efficient chemical ozone destruction. Other processes, such as the acceleration of the Brewer–Dobson circulation, also affect stratospheric temperatures, further complicating the picture. Therefore, it is currently unclear whether major Arctic ozone minima will still occur at the end of the 21st century despite decreasing hODSs. We have examined this question for different emission pathways using simulations conducted within the Chemistry-Climate Model Initiative (CCMI-1 and CCMI-2022) and found large differences in the models' ability to simulate the magnitude of ozone minima in the present-day climate. Models with a generally too-cold polar stratosphere (cold bias) produce pronounced ozone minima under present-day climate conditions because they simulate more PSCs and, thus, high concentrations of active chlorine species (ClOx). These models predict the largest decrease in ozone minima in the future. Conversely, models with a warm polar stratosphere (warm bias) have the smallest sensitivity of ozone minima to future changes in hODS and GHG concentrations. As a result, the scatter among models in terms of the magnitude of Arctic spring ozone minima will decrease in the future. Overall, these results suggest that Arctic ozone minima will become weaker over the next decades, largely due to the decline in hODS abundances. We note that none of the models analysed here project a notable increase of ozone minima in the future. Stratospheric cooling caused by increasing GHG concentrations is expected to play a secondary role as its effect in the Arctic stratosphere is weakened by opposing radiative and dynamical mechanisms. [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]

Ozone in the Arctic stratosphere is subject to large interannual variability, driven by both chemical ozone depletion and dynamical variability. Anomalies in Arctic stratospheric ozone become particularly important in spring, when returning sunlight allows them to alter stratospheric temperatures via shortwave heating, thus modifying atmospheric dynamics. At the same time, the stratospheric circulation undergoes a transition in spring with the final stratospheric warming (FSW), which marks the end of winter. A causal link between stratospheric ozone anomalies and FSWs is plausible and might increase the predictability of stratospheric and tropospheric responses on sub-seasonal to seasonal timescales. However, it remains to be fully understood how ozone influences the timing and evolution of the springtime vortex breakdown. Here, we contrast results from chemistry climate models with and without interactive ozone chemistry to quantify the impact of ozone anomalies on the timing of the FSW and its effects on surface climate. We find that ozone feedbacks increase the variability in the timing of the FSW, especially in the lower stratosphere. In ozone-deficient springs, a persistent strong polar vortex and a delayed FSW in the lower stratosphere are partly due to the lack of heating by ozone in that region. High-ozone anomalies, on the other hand, result in additional shortwave heating in the lower stratosphere, where the FSW therefore occurs earlier. We further show that FSWs in high-ozone springs are predominantly followed by a negative phase of the Arctic Oscillation (AO) with positive sea level pressure anomalies over the Arctic and cold anomalies over Eurasia and Europe. These conditions are to a significant extent (at least 50 %) driven by ozone. In contrast, FSWs in low-ozone springs are not associated with a discernible surface climate response. These results highlight the importance of ozone–circulation coupling in the climate system and the potential value of interactive ozone chemistry for sub-seasonal to seasonal predictability. [ABSTRACT FROM AUTHOR]

Ozone in the Arctic stratosphere is subject to large interannual variability, driven by both chemical ozone depletion and dynamical variability. Anomalies in Arctic stratospheric ozone become particularly important in spring, when returning sunlight allows them to alter stratospheric temperatures via shortwave heating, thus modifying atmospheric dynamics. At the same time, the stratospheric circulation undergoes a transition in spring with the Stratospheric Final Warming (FSW), which marks the end of winter. A causal link between stratospheric ozone anomalies and FSWs is plausible and might increase the predictability of stratospheric and tropospheric responses on sub-seasonal to seasonal timescales. However, it remains to be fully understood how ozone influences the timing and evolution of the springtime vortex breakdown. Here, we contrast results from chemistry climate models with and without interactive ozone chemistry to quantify the impact of ozone anomalies on the timing of the FSW and its effects on surface climate. We find that ozone feedbacks increase the variability in the timing of the FSW, especially in the lower stratosphere. In ozone-deficient springs, a persistent strong polar vortex and a delayed FSWin the lower stratosphere are partly due to lacking heating by ozone in that region. High ozone anomalies, on the other hand, result in additional shortwave heating in the lower stratosphere, where the FSW therefore occurs earlier. We further show that FSWs in high ozone springs are predominantly followed by a negative phase of the Arctic Oscillation (AO) with positive sea level pressure anomalies over the Arctic and cold anomalies over Eurasia and Europe. These conditions are to a significant extent (at least 50%) driven by ozone. In contrast, FSWs in low ozone springs are not associated with a discernible surface climate response. These results highlight the importance of ozone-circulation coupling in the climate system and the potential value of interactive ozone chemistry for sub-seasonal to seasonal predictability. [ABSTRACT FROM AUTHOR]

This paper features the new atmosphere–ocean–aerosol–chemistry–climate model, SOlar Climate Ozone Links (SOCOL) v4.0, and its validation. The new model was built by interactively coupling the Max Planck Institute Earth System Model version 1.2 (MPI-ESM1.2) (T63, L47) with the chemistry (99 species) and size-resolving (40 bins) sulfate aerosol microphysics modules from the aerosol–chemistry–climate model, SOCOL-AERv2. We evaluate its performance against reanalysis products and observations of atmospheric circulation, temperature, and trace gas distribution, with a focus on stratospheric processes. We show that SOCOLv4.0 captures the low- and midlatitude stratospheric ozone well in terms of the climatological state, variability and evolution. The model provides an accurate representation of climate change, showing a global surface warming trend consistent with observations as well as realistic cooling in the stratosphere caused by greenhouse gas emissions, although, as in previous model versions, a too-fast residual circulation and exaggerated mixing in the surf zone are still present. The stratospheric sulfur budget for moderate volcanic activity is well represented by the model, albeit with slightly underestimated aerosol lifetime after major eruptions. The presence of the interactive ocean and a successful representation of recent climate and ozone layer trends make SOCOLv4.0 ideal for studies devoted to future ozone evolution and effects of greenhouse gases and ozone-destroying substances, as well as the evaluation of potential solar geoengineering measures through sulfur injections. Potential further model improvements could be to increase the vertical resolution, which is expected to allow better meridional transport in the stratosphere, as well as to update the photolysis calculation module and budget of mesospheric odd nitrogen. In summary, this paper demonstrates that SOCOLv4.0 is well suited for applications related to the stratospheric ozone and sulfate aerosol evolution, including its participation in ongoing and future model intercomparison projects. [ABSTRACT FROM AUTHOR]

As part of the Model Intercomparison Project on the climatic response to Volcanic forcing (VolMIP), several climate modeling centers performed a coordinated pre-study experiment with interactive stratospheric aerosol models simulating the volcanic aerosol cloud from an eruption resembling the 1815 Mt. Tambora eruption (VolMIP-Tambora ISA ensemble). The pre-study provided the ancillary ability to assess intermodel diversity in the radiative forcing for a large stratospheric-injecting equatorial eruption when the volcanic aerosol cloud is simulated interactively. An initial analysis of the VolMIP-Tambora ISA ensemble showed large disparities between models in the stratospheric global mean aerosol optical depth (AOD). In this study, we now show that stratospheric global mean AOD differences among the participating models are primarily due to differences in aerosol size, which we track here by effective radius. We identify specific physical and chemical processes that are missing in some models and/or parameterized differently between models, which are together causing the differences in effective radius. In particular, our analysis indicates that interactively tracking hydroxyl radical (OH) chemistry following a large volcanic injection of sulfur dioxide (SO2) is an important factor in allowing for the timescale for sulfate formation to be properly simulated. In addition, depending on the timescale of sulfate formation, there can be a large difference in effective radius and subsequently AOD that results from whether the SO2 is injected in a single model grid cell near the location of the volcanic eruption, or whether it is injected as a longitudinally averaged band around the Earth. [ABSTRACT FROM AUTHOR]

Polar stratospheric clouds (PSCs) contribute to catalytic ozone destruction by providing surfaces for the conversion of inert chlorine species into active forms and by denitrification. The latter describes the removal of HNO3 from the stratosphere by sedimenting PSC particles, which hinders chlorine deactivation by the formation of reservoir species. Therefore, an accurate representation of PSCs in chemistry–climate models (CCMs) is of great importance to correctly simulate polar ozone concentrations. Here, we evaluate PSCs as simulated by the CCM SOCOLv3.1 for the Antarctic winters 2006, 2007 and 2010 by comparison with backscatter measurements by CALIOP on board the CALIPSO satellite. The year 2007 represents a typical Antarctic winter, while 2006 and 2010 are characterized by above- and below-average PSC occurrence. The model considers supercooled ternary solution (STS) droplets, nitric acid trihydrate (NAT) particles, water ice particles and mixtures thereof. PSCs are parameterized in terms of temperature and partial pressures of HNO3 and H2O , assuming equilibrium between the gas and particulate phase. The PSC scheme involves a set of prescribed microphysical parameters, namely ice number density, NAT particle radius and maximum NAT number density. In this study, we test and optimize the parameter settings through several sensitivity simulations. The choice of the value for the ice number density affects simulated optical properties and dehydration, while modifying the NAT parameters impacts stratospheric composition via HNO3 uptake and denitrification. Depending on the NAT parameters, reasonable denitrification can be modeled. However, its impact on ozone loss is minor. The best agreement with the CALIOP optical properties and observed denitrification was for this case study found with the ice number density increased from the hitherto used value of 0.01 to 0.05 cm-3 and the maximum NAT number density from 5×10-4 to 1×10-3 cm-3. The NAT radius was kept at the original value of 5 µm. The new parameterization reflects the higher importance attributed to heterogeneous nucleation of ice and NAT particles following recent new data evaluations of the state-of-the-art CALIOP measurements. A cold temperature bias in the polar lower stratosphere results in an overestimated PSC areal coverage in SOCOLv3.1 by up to 40 %. Offsetting this cold bias by + 3 K delays the onset of ozone depletion by about 2 weeks, which improves the agreement with observations. Furthermore, the occurrence of mountain-wave-induced ice, as observed mainly over the Antarctic Peninsula, is continuously underestimated in the model due to the coarse model resolution (T42L39) and the fixed ice number density. Nevertheless, we find overall good temporal and spatial agreement between modeled and observed PSC occurrence and composition. This work confirms previous studies indicating that simplified PSC schemes, which avoid nucleation and growth calculations in sophisticated but time-consuming microphysical process models, may also achieve good approximations of the fundamental properties of PSCs needed in CCMs. [ABSTRACT FROM AUTHOR]

Water vapor (H 2 O) is the source of reactive hydrogen radicals in the middle atmosphere, whereas carbon monoxide (CO), being formed by CO 2 photolysis, is suitable as a dynamical tracer. In the mesosphere, both H 2 O and CO are sensitive to solar irradiance (SI) variability because of their destruction/production by solar radiation. This enables us to analyze the solar signal in both models and observed data. Here, we evaluate the mesospheric H 2 O and CO response to solar irradiance variability using the Chemistry-Climate Model Initiative (CCMI-1) simulations and satellite observations. We analyzed the results of four CCMI models (CMAM, EMAC-L90MA, SOCOLv3, and CESM1-WACCM 3.5) operated in CCMI reference simulation REF-C1SD in specified dynamics mode, covering the period from 1984–2017. Multiple linear regression analyses show a pronounced and statistically robust response of H 2 O and CO to solar irradiance variability and to the annual and semiannual cycles. For periods with available satellite data, we compared the simulated solar signal against satellite observations, namely the GOZCARDS composite for 1992–2017 for H 2 O and Aura/MLS measurements for 2005–2017 for CO. The model results generally agree with observations and reproduce an expected negative and positive correlation for H 2 O and CO, respectively, with solar irradiance. However, the magnitude of the response and patterns of the solar signal varies among the considered models, indicating differences in the applied chemical reaction and dynamical schemes, including the representation of photolyzes. We suggest that there is no dominating thermospheric influence of solar irradiance in CO, as reported in previous studies, because the response to solar variability is comparable with observations in both low-top and high-top models. We stress the importance of this work for improving our understanding of the current ability and limitations of state-of-the-art models to simulate a solar signal in the chemistry and dynamics of the middle atmosphere. [ABSTRACT FROM AUTHOR]

As part of the Model Intercomparison Project on the climatic response to Volcanic forcing (VolMIP), several climate modeling centers performed a coordinated pre-study experiment with interactive stratospheric aerosol models simulating the volcanic aerosol cloud from an eruption resembling the 1815 Mt Tambora eruption (VolMIP-Tambora ISA ensemble). The pre-study provided the ancillary ability to assess intermodel diversity in the radiative forcing for a large stratospheric-injecting equatorial eruption when the volcanic aerosol cloud is simulated interactively. An initial analysis of the VolMIP-Tambora ISA ensemble showed large disparities between models in the stratospheric global mean aerosol optical depth (AOD). In this study, we now show that stratospheric global mean AOD differences among the participating models are primarily due to differences in aerosol size, which we track here by effective radius. We identify specific physical and chemical processes that are missing in some models and/or parameterized differently between models, which are together causing the differences in effective radius. In particular, our analysis indicates that interactively tracking hydroxyl radical (OH) chemistry following a large volcanic injection of sulfur dioxide (SO2) is an important factor in allowing for the timescale for sulfate formation to be properly simulated. In addition, depending on the timescale of sulfate formation, there can be a large difference in effective radius and subsequently AOD that results from whether the SO2 is injected in a single model gridcell near the location of the volcanic eruption, or whether it is injected as a longitudinally averaged band around the Earth. [ABSTRACT FROM AUTHOR]

Water vapor (H2O) is the source of reactive hydrogen radicals in the middle atmosphere, whereas carbon monoxide (CO), being formed by CO2 photolysis, is suitable as a dynamical tracer. In the mesosphere, both H2O and CO are sensitive to solar irradiance variability because of their destruction/production by solar radiation. This enables to analyze the solar signal in both, models and observed data. Here, we evaluate the mesospheric H2O and CO response to solar irradiance variability using the Chemistry-Climate Model Initiative (CCMI-1) simulations and satellite data. We analyzed the results of four CCMI models (CMAM, EMAC-L90MA, SOCOLv3, CESM1-WACCM 3.5) operated in CCMI reference simulation REF-C1SD in specified dynamics mode, covering the period from 1984 to 2017. Multiple linear regression analysis shows a pronounced and statistically robust response of H2O and CO to solar irradiance variability, and to the annual and semiannual cycles. For periods with available satellite data, we compared the simulated solar signal against satellite observations, namely during 1992-2017 for H2O and 2005-2017 for CO. The model results generally agree with observations and reproduce an expected negative and positive correlation for H2O and CO, respectively, with solar irradiance. However, the magnitude of the response and patterns of the solar signal varies among the considered models, indicating differences in the applied chemical reaction and dynamical schemes including the representation of photolyses. We suggest that there is no dominating thermospheric influence of solar irradiance in CO, as reported in previous studies because the response to solar variability is comparable with observations in both, low-top and high-top models. We stress the importance of this work for improving our understanding of the current ability and limitations of state-of-the-art models to simulate a solar signal in the chemistry and dynamic of the middle atmosphere. [ABSTRACT FROM AUTHOR]

The stratospheric ozone layer shields surface life from harmful ultraviolet radiation. Following the Montreal Protocol ban on long-lived ozone-depleting substances (ODSs), rapid depletion of total column ozone (TCO) ceased in the late 1990s, and ozone above 32 km is now clearly recovering. However, there is still no confirmation of TCO recovery, and evidence has emerged that ongoing quasi-global (60 ∘ S–60 ∘ N) lower stratospheric ozone decreases may be responsible, dominated by low latitudes (30 ∘ S–30 ∘ N). Chemistry–climate models (CCMs) used to project future changes predict that lower stratospheric ozone will decrease in the tropics by 2100 but not at mid-latitudes (30–60 ∘). Here, we show that CCMs display an ozone decline similar to that observed in the tropics over 1998–2016, likely driven by an increase in tropical upwelling. On the other hand, mid-latitude lower stratospheric ozone is observed to decrease, while CCMs that specify real-world historical meteorological fields instead show an increase up to present day. However, these cannot be used to simulate future changes; we demonstrate here that free-running CCMs used for projections also show increases. Despite opposing lower stratospheric ozone changes, which should induce opposite temperature trends, CCMs and observed temperature trends agree; we demonstrate that opposing model–observation stratospheric water vapour (SWV) trends, and their associated radiative effects, explain why temperature changes agree in spite of opposing ozone trends. We provide new evidence that the observed mid-latitude trends can be explained by enhanced mixing between the tropics and extratropics. We further show that the temperature trends are consistent with the observed mid-latitude ozone decrease. Together, our results suggest that large-scale circulation changes expected in the future from increased greenhouse gases (GHGs) may now already be underway but that most CCMs do not simulate mid-latitude ozone layer changes well. However, it is important to emphasise that the periods considered here are short, and internal variability that is both intrinsic to each CCM and different to observed historical variability is not well-characterised and can influence trend estimates. Nevertheless, the reason CCMs do not exhibit the observed changes needs to be identified to allow models to be improved in order to build confidence in future projections of the ozone layer. [ABSTRACT FROM AUTHOR]

Polar Stratospheric Clouds (PSCs) contribute to catalytic ozone destruction by providing surfaces for the conversion of inert chlorine species into active forms and by denitrification of the stratosphere. Therefore, an accurate representation of PSCs in chemistry-climate models (CCMs) is of great importance to correctly simulate polar ozone concentrations. Here, we evaluate PSCs as simulated by the CCM SOCOLv3.1 for the Antarctic winter 2007 by comparison with backscatter measurements by CALIOP onboard the CALIPSO satellite. The model considers supercooled ternary solution (STS) droplets, nitric acid trihydrate (NAT) particles, water ice particles, and mixtures thereof. PSCs are parametrized in terms of temperature and partial pressures of HNO3 and H2O, assuming equilibrium between gas and particulate phase. We use the CALIOP measurements to optimize three prescribed microphysical parameters of the PSC scheme, namely ice number density, NAT particle radius and maximum NAT number density. The choice of the prescribed value of the ice number density affects simulated optical properties and dehydration, while modifying the maximum NAT number density or the NAT particle radius impacts stratospheric composition by enhancing the HNO3-uptake and denitrification. Best agreement with the CALIOP optical properties and observed denitrification was for this case study found with the ice number density increased from the hitherto used value of 0.01 to 0.05 cm−3 and the maximum NAT number density from 5×10−4 to 1×10−3 cm−3. The NAT radius was kept at the original value of 5 µm. The new parametrization reflects the higher importance attributed to heterogeneous nucleation of ice and NAT particles, e.g. on meteoric dust, following recent new data evaluations of the state-of-the-art CALIOP measurements. A cold temperature bias in the polar lower stratosphere results in an overestimated PSC areal coverage in SOCOLv3.1 by up to 100 %. Furthermore, the occurrence of mountain-wave induced ice, as observed mainly over the Antarctic Peninsula, is continuously underestimated in the model due to the coarse model resolution and the fixed ice number density. However, overall we find a good temporal and spatial agreement between modeled and observed PSC occurrence and composition, as well as reasonable modeled denitrification and ozone loss. Based on constraining three important parameters by means of the CALIOP measurements, this work demonstrates that also a simplified PSC scheme, which describes STS, NAT, ice and mixtures thereof with equilibrium assumptions and avoids nucleation and growth calculations in sophisticated, but time-consuming microphysical process models, may achieve good approximations of fundamental properties of PSCs needed in CCMs. [ABSTRACT FROM AUTHOR]

An estimated 0.5–1 billion people globally have inadequate intakes of selenium (Se), due to a lack of bioavailable Se in agricultural soils. Deposition from the atmosphere, especially through precipitation, is an important source of Se to soils. However, very little is known about the atmospheric cycling of Se. It has therefore been difficult to predict how far Se travels in the atmosphere and where it deposits. To answer these questions, we have built the first global atmospheric Se model by implementing Se chemistry in an aerosol–chemistry–climate model, SOCOL-AER (modeling tools for studies of SOlar Climate Ozone Links – aerosol). In the model, we include information from the literature about the emissions, speciation, and chemical transformation of atmospheric Se. Natural processes and anthropogenic activities emit volatile Se compounds, which oxidize quickly and partition to the particulate phase. Our model tracks the transport and deposition of Se in seven gas-phase species and 41 aerosol tracers. However, there are large uncertainties associated with many of the model's input parameters. In order to identify which model uncertainties are the most important for understanding the atmospheric Se cycle, we conducted a global sensitivity analysis with 34 input parameters related to Se chemistry, Se emissions, and the interaction of Se with aerosols. In the first bottom-up estimate of its kind, we have calculated a median global atmospheric lifetime of 4.4 d (days), ranging from 2.9 to 6.4 d (2nd–98th percentile range) given the uncertainties of the input parameters. The uncertainty in the Se lifetime is mainly driven by the uncertainty in the carbonyl selenide (OCSe) oxidation rate and the lack of tropospheric aerosol species other than sulfate aerosols in SOCOL-AER. In contrast to uncertainties in Se lifetime, the uncertainty in deposition flux maps are governed by Se emission factors, with all four Se sources (volcanic, marine biosphere, terrestrial biosphere, and anthropogenic emissions) contributing equally to the uncertainty in deposition over agricultural areas. We evaluated the simulated Se wet deposition fluxes from SOCOL-AER with a compiled database of rainwater Se measurements, since wet deposition contributes around 80 % of total Se deposition. Despite difficulties in comparing a global, coarse-resolution model with local measurements from a range of time periods, past Se wet deposition measurements are within the range of the model's 2nd–98th percentiles at 79 % of background sites. This agreement validates the application of the SOCOL-AER model to identifying regions which are at risk of low atmospheric Se inputs. In order to constrain the uncertainty in Se deposition fluxes over agricultural soils, we should prioritize field campaigns measuring Se emissions, rather than laboratory measurements of Se rate constants. [ABSTRACT FROM AUTHOR]

The top-of-atmosphere (TOA) outgoing longwave flux over the 9.6 µm ozone band is a fundamental quantity for understanding chemistry–climate coupling. However, observed TOA fluxes are hard to estimate as they exhibit considerable variability in space and time that depend on the distributions of clouds, ozone (O3), water vapor (H2O), air temperature (Ta), and surface temperature (Ts). Benchmarking present-day fluxes and quantifying the relative influence of their drivers is the first step for estimating climate feedbacks from ozone radiative forcing and predicting radiative forcing evolution. To that end, we constructed observational instantaneous radiative kernels (IRKs) under clear-sky conditions, representing the sensitivities of the TOA flux in the 9.6 µm ozone band to the vertical distribution of geophysical variables, including O3 , H2O , Ta , and Ts based upon the Aura Tropospheric Emission Spectrometer (TES) measurements. Applying these kernels to present-day simulations from the Chemistry-Climate Model Initiative (CCMI) project as compared to a 2006 reanalysis assimilating satellite observations, we show that the models have large differences in TOA flux, attributable to different geophysical variables. In particular, model simulations continue to diverge from observations in the tropics, as reported in previous studies of the Atmospheric Chemistry Climate Model Intercomparison Project (ACCMIP) simulations. The principal culprits are tropical middle and upper tropospheric ozone followed by tropical lower tropospheric H2O. Five models out of the eight studied here have TOA flux biases exceeding 100 mW m -2 attributable to tropospheric ozone bias. Another set of five models have flux biases over 50 mW m -2 due to H2O. On the other hand, Ta radiative bias is negligible in all models (no more than 30 mW m -2). We found that the atmospheric component (AM3) of the Geophysical Fluid Dynamics Laboratory (GFDL) general circulation model and Canadian Middle Atmosphere Model (CMAM) have the lowest TOA flux biases globally but are a result of cancellation of opposite biases due to different processes. Overall, the multi-model ensemble mean bias is -133±98 mW m -2 , indicating that they are too atmospherically opaque due to trapping too much radiation in the atmosphere by overestimated tropical tropospheric O3 and H2O. Having too much O3 and H2O in the troposphere would have different impacts on the sensitivity of TOA flux to O3 and these competing effects add more uncertainties on the ozone radiative forcing. We find that the inter-model TOA outgoing longwave radiation (OLR) difference is well anti-correlated with their ozone band flux bias. This suggests that there is significant radiative compensation in the calculation of model outgoing longwave radiation. [ABSTRACT FROM AUTHOR]

SOCOL-AERv1 was developed as an aerosol–chemistry–climate model to study the stratospheric sulfur cycle and its influence on climate and the ozone layer. It includes a sectional aerosol model that tracks the sulfate particle size distribution in 40 size bins, between 0.39 nm and 3.2 µm. showed that SOCOL-AERv1 successfully matched observable quantities related to stratospheric aerosol. In the meantime, SOCOL-AER has undergone significant improvements and more observational datasets have become available. In producing SOCOL-AERv2 we have implemented several updates to the model: adding interactive deposition schemes, improving the sulfate mass and particle number conservation, and expanding the tropospheric chemistry scheme. We compare the two versions of the model with background stratospheric sulfate aerosol observations, stratospheric aerosol evolution after Pinatubo, and ground-based sulfur deposition networks. SOCOL-AERv2 shows similar levels of agreement as SOCOL-AERv1 with satellite-measured extinctions and in situ optical particle counter (OPC) balloon flights. The volcanically quiescent total stratospheric aerosol burden simulated in SOCOL-AERv2 has increased from 109 Gg of sulfur (S) to 160 Gg S, matching the newly available satellite estimate of 165 Gg S. However, SOCOL-AERv2 simulates too high cross-tropopause transport of tropospheric SO2 and/or sulfate aerosol, leading to an overestimation of lower stratospheric aerosol. Due to the current lack of upper tropospheric SO2 measurements and the neglect of organic aerosol in the model, the lower stratospheric bias of SOCOL-AERv2 was not further improved. Model performance under volcanically perturbed conditions has also undergone some changes, resulting in a slightly shorter volcanic aerosol lifetime after the Pinatubo eruption. With the improved deposition schemes of SOCOL-AERv2, simulated sulfur wet deposition fluxes are within a factor of 2 of measured deposition fluxes at 78 % of the measurement stations globally, an agreement which is on par with previous model intercomparison studies. Because of these improvements, SOCOL-AERv2 will be better suited to studying changes in atmospheric sulfur deposition due to variations in climate and emissions. [ABSTRACT FROM AUTHOR]

We perform the first multi-model intercomparison of the impact of nudged meteorology on the stratospheric residual circulation using hindcast simulations from the Chemistry–Climate Model Initiative (CCMI). We examine simulations over the period 1980–2009 from seven models in which the meteorological fields are nudged towards a reanalysis dataset and compare these with their equivalent free-running simulations and the reanalyses themselves. We show that for the current implementations, nudging meteorology does not constrain the mean strength of the stratospheric residual circulation and that the inter-model spread is similar, or even larger, than in the free-running simulations. The nudged models generally show slightly stronger upwelling in the tropical lower stratosphere compared to the free-running versions and exhibit marked differences compared to the directly estimated residual circulation from the reanalysis dataset they are nudged towards. Downward control calculations applied to the nudged simulations reveal substantial differences between the climatological lower-stratospheric tropical upward mass flux (TUMF) computed from the modelled wave forcing and that calculated directly from the residual circulation. This explicitly shows that nudging decouples the wave forcing and the residual circulation so that the divergence of the angular momentum flux due to the mean motion is not balanced by eddy motions, as would typically be expected in the time mean. Overall, nudging meteorological fields leads to increased inter-model spread for most of the measures of the mean climatological stratospheric residual circulation assessed in this study. In contrast, the nudged simulations show a high degree of consistency in the inter-annual variability in the TUMF in the lower stratosphere, which is primarily related to the contribution to variability from the resolved wave forcing. The more consistent inter-annual variability in TUMF in the nudged models also compares more closely with the variability found in the reanalyses, particularly in boreal winter. We apply a multiple linear regression (MLR) model to separate the drivers of inter-annual and long-term variations in the simulated TUMF; this explains up to ∼75 % of the variance in TUMF in the nudged simulations. The MLR model reveals a statistically significant positive trend in TUMF for most models over the period 1980–2009. The TUMF trend magnitude is generally larger in the nudged models compared to their free-running counterparts, but the intermodel range of trends doubles from around a factor of 2 to a factor of 4 due to nudging. Furthermore, the nudged models generally do not match the TUMF trends in the reanalysis they are nudged towards for trends over different periods in the interval 1980–2009. Hence, we conclude that nudging does not strongly constrain long-term trends simulated by the chemistry–climate model (CCM) in the residual circulation. Our findings show that while nudged simulations may, by construction, produce accurate temperatures and realistic representations of fast horizontal transport, this is not typically the case for the slower zonal mean vertical transport in the stratosphere. Consequently, caution is required when using nudged simulations to interpret the behaviour of stratospheric tracers that are affected by the residual circulation. [ABSTRACT FROM AUTHOR]

An estimated 0.5–1 billion people globally have inadequate intakes of selenium (Se), due to a lack of bioavailable Se in agricultural soils. Deposition from the atmosphere, especially through precipitation, is an important source of Se to soils. However, very little is known about the atmospheric cycling of Se. It has therefore been difficult to predict how far Se travels in the atmosphere and where it deposits. To answer these questions, we have built the first global atmospheric Se model by implementing Se chemistry into an aerosol–chemistry–climate model, SOCOL-AER. In the model, we include information from the literature about the emissions, speciation, and chemical transformation of atmospheric Se. Natural processes and anthropogenic activities emit volatile Se compounds, which oxidize quickly and partition to the particulate phase. Our model tracks the transport and deposition of Se in 7 gas-phase species and 41 aerosol tracers. However, there are large uncertainties associated with many of the model's input parameters. In order to identify which model uncertainties are the most important for understanding the atmospheric Se cycle, we conducted a global sensitivity analysis with 34 input parameters related to Se chemistry, Se emissions, and the interaction of Se with aerosols. In the first bottom-up estimate of its kind, we have calculated a median global atmospheric lifetime of 4.4 d (days), ranging from 2.9–6.4 d (2nd–98th percentile) given the uncertainties of the input parameters. The uncertainty in the Se lifetime is mainly driven by the uncertainty in the carbonyl selenide (OCSe) oxidation rate and the lack of tropospheric aerosol species other than sulfate aerosols in SOCOL-AER. In contrast to uncertainties in Se lifetime, the uncertainty in deposition flux maps are governed by Se emission factors, with all four Se sources (volcanic, marine biosphere, terrestrial biosphere, and anthropogenic emissions) contributing equally to the uncertainty in deposition over agricultural areas. We evaluated the simulated Se wet deposition fluxes from SOCOL-AER with a compiled database of rainwater Se measurements, since wet deposition contributes around 80 % of total Se deposition. Despite difficulties in comparing a global, coarse resolution model with local measurements from a range of time periods, past Se wet deposition measurements are within the range of the model's 2nd–98th percentile at 79 % of background sites. This agreement validates the application of the SOCOL-AER model to identifying regions which are at risk of low atmospheric Se inputs. In order to constrain the uncertainty in Se deposition fluxes over agricultural soils we should prioritize field campaigns measuring Se emissions, rather than laboratory measurements of Se rate constants. [ABSTRACT FROM AUTHOR]

We have derived values of the ultraviolet index (UVI) at solar noon using the Tropospheric Ultraviolet Model (TUV) driven by ozone, temperature and aerosol fields from climate simulations of the first phase of the Chemistry-Climate Model Initiative (CCMI-1). Since clouds remain one of the largest uncertainties in climate projections, we simulated only the clear-sky UVI. We compared the modelled UVI climatologies against present-day climatological values of UVI derived from both satellite data (the OMI-Aura OMUVBd product) and ground-based measurements (from the NDACC network). Depending on the region, relative differences between the UVI obtained from CCMI/TUV calculations and the ground-based measurements ranged between -5.9 % and 10.6 %. We then calculated the UVI evolution throughout the 21st century for the four Representative Concentration Pathways (RCPs 2.6, 4.5, 6.0 and 8.5). Compared to 1960s values, we found an average increase in the UVI in 2100 (of 2 %–4 %) in the tropical belt (30 ∘ N–30 ∘ S). For the mid-latitudes, we observed a 1.8 % to 3.4 % increase in the Southern Hemisphere for RCPs 2.6, 4.5 and 6.0 and found a 2.3 % decrease in RCP 8.5. Higher increases in UVI are projected in the Northern Hemisphere except for RCP 8.5. At high latitudes, ozone recovery is well identified and induces a complete return of mean UVI levels to 1960 values for RCP 8.5 in the Southern Hemisphere. In the Northern Hemisphere, UVI levels in 2100 are higher by 0.5 % to 5.5 % for RCPs 2.6, 4.5 and 6.0 and they are lower by 7.9 % for RCP 8.5. We analysed the impacts of greenhouse gases (GHGs) and ozone-depleting substances (ODSs) on UVI from 1960 by comparing CCMI sensitivity simulations (1960–2100) with fixed GHGs or ODSs at their respective 1960 levels. As expected with ODS fixed at their 1960 levels, there is no large decrease in ozone levels and consequently no sudden increase in UVI levels. With fixed GHG, we observed a delayed return of ozone to 1960 values, with a corresponding pattern of change observed on UVI, and looking at the UVI difference between 2090s values and 1960s values, we found an 8 % increase in the tropical belt during the summer of each hemisphere. Finally we show that, while in the Southern Hemisphere the UVI is mainly driven by total ozone column, in the Northern Hemisphere both total ozone column and aerosol optical depth drive UVI levels, with aerosol optical depth having twice as much influence on the UVI as total ozone column does. [ABSTRACT FROM AUTHOR]

Energetic particle precipitation (EPP) affects the chemistry of the polar middle atmosphere by producing reactive nitrogen (NOy) and hydrogen (HOx) species, which then catalytically destroy ozone. Recently, there have been major advances in constraining these particle impacts through a parametrization of NOy based on high-quality observations. Here we investigate the effects of low (auroral) and middle (radiation belt) energy range electrons, separately and in combination, on reactive nitrogen and hydrogen species as well as on ozone during Southern Hemisphere winters from 2002 to 2010 using the SOCOL3-MPIOM chemistry-climate model. Our results show that, in the absence of solar proton events, low-energy electrons produce the majority of NOy in the polar mesosphere and stratosphere. In the polar vortex, NOy subsides and affects ozone at lower altitudes, down to 10 hPa. Comparing a year with high electron precipitation with a quiescent period, we found large ozone depletion in the mesosphere; as the anomaly propagates downward, 15 % less ozone is found in the stratosphere during winter, which is confirmed by satellite observations. Only with both low- and middle-energy electrons does our model reproduce the observed stratospheric ozone anomaly. [ABSTRACT FROM AUTHOR]

The top-of-atmosphere (TOA) outgoing longwave flux over the 9.6-μm ozone band is a fundamental quantity for understanding chemistry-climate coupling. However, observed TOA fluxes are hard to estimate as they exhibit considerable variability in space and time that depend on the distributions of clouds, ozone (O3), water vapor (H2O), air temperature (Ta), and surface temperature (Ts). Benchmarking present day fluxes and quantifying the relative influence of their drivers is the first step for estimating climate feedbacks from ozone radiative forcing and predicting its evolution. To that end, we construct observational instantaneous radiative kernels (IRKs) representing the sensitivities of the TOA flux in the 9.6-μm ozone band to the vertical distribution of geophysical variables, including O3, H2O, Ta, and Ts based upon the Aura Tropospheric Emission Spectrometer (TES) measurements. Applying these kernels to present-day simulations from the Chemistry-Climate Model Initiative (CCMI) project as compared to a 2006 reanalysis assimilating satellite observations, we show that the models have large differences in TOA flux, attributable to different geophysical variables. In particular, model simulations continue to diverge from observations in the tropics, as reported in previous studies of the Atmospheric Chemistry Climate Model Inter-comparison Project (ACCMIP) simulations. The principal culprits are tropical mid and upper tropospheric ozone followed by tropical lower tropospheric H2O. Five models out of the eight studied here have TOA flux biases exceeding 100 mWm−2 attributable to tropospheric ozone bias. Another set of five models flux biases over 50 mWm−2 due to H2O. On the other hand, Ta radiative bias is negligible in all models (no more than 30 mWm−2). We found that AM3 and CMAM have the lowest TOA flux biases globally but are a result of cancellation of difference processes. Overall, the multi-model ensemble mean bias is −132.9 ± 98 mWm−2, indicating that they are too atmospherically opaque thereby reducing sensitivity of TOA flux to ozone and potentially an underestimate of ozone radiative forcing. We find that the inter-model TOA OLR difference is well anti-correlated with their ozone band flux bias. This suggests that there is significant radiative compensation in the calculation of model outgoing longwave radiation. [ABSTRACT FROM AUTHOR]

Stratospheric sulfate geoengineering (SSG) could contribute to avoiding some of the adverse impacts of climate change. We used the SOCOL-AER global aerosol–chemistry–climate model to investigate 21 different SSG scenarios, each with 1.83 Mt S yr -1 injected either in the form of accumulation-mode H2SO4 droplets (AM H2SO4), gas-phase SO2 or as combinations of both. For most scenarios, the sulfur was continuously emitted at an altitude of 50 hPa (≈20 km) in the tropics and subtropics. We assumed emissions to be zonally and latitudinally symmetric around the Equator. The spread of emissions ranged from 3.75 ∘ S–3.75 ∘ N to 30 ∘ S–30 ∘ N. In the SO2 emission scenarios, continuous production of tiny nucleation-mode particles results in increased coagulation, which together with gaseous H2SO4 condensation, produces coarse-mode particles. These large particles are less effective for backscattering solar radiation and have a shorter stratospheric residence time than AM H2SO4 particles. On average, the stratospheric aerosol burden and corresponding all-sky shortwave radiative forcing for the AM H2SO4 scenarios are about 37 % larger than for the SO2 scenarios. The simulated stratospheric aerosol burdens show a weak dependence on the latitudinal spread of emissions. Emitting at 30 ∘ N–30 ∘ S instead of 10 ∘ N–10 ∘ S only decreases stratospheric burdens by about 10 %. This is because a decrease in coagulation and the resulting smaller particle size is roughly balanced by faster removal through stratosphere-to-troposphere transport via tropopause folds. Increasing the injection altitude is also ineffective, although it generates a larger stratospheric burden, because enhanced condensation and/or coagulation leads to larger particles, which are less effective scatterers. In the case of gaseous SO2 emissions, limiting the sulfur injections spatially and temporally in the form of point and pulsed emissions reduces the total global annual nucleation, leading to less coagulation and thus smaller particles with increased stratospheric residence times. Pulse or point emissions of AM H2SO4 have the opposite effect: they decrease the stratospheric aerosol burden by increasing coagulation and only slightly decrease clear-sky radiative forcing. This study shows that direct emission of AM H2SO4 results in higher radiative forcing for the same sulfur equivalent mass injection strength than SO2 emissions, and that the sensitivity to different injection strategies varies for different forms of injected sulfur. [ABSTRACT FROM AUTHOR]

Previous multi-model intercomparisons have shown that chemistry-climate models exhibit significant biases in tropospheric ozone compared with observations. We investigate annual-mean tropospheric column ozone in 15 models participating in the SPARC-IGAC (Stratosphere-troposphere Processes And their Role in Climate-International Global Atmospheric Chemistry) Chemistry-Climate Model Initiative (CCMI). These models exhibit a positive bias, on average, of up to 40%-50% in the Northern Hemisphere compared with observations derived from the Ozone Monitoring Instrument and Microwave Limb Sounder (OMI/MLS), and a negative bias of up to ∼ 30% in the Southern Hemisphere. SOCOLv3.0 (version 3 of the Solar-Climate Ozone Links CCM), which participated in CCMI, simulates global-mean tropospheric ozone columns of 40.2DU - approximately 33% larger than the CCMI multi-model mean. Here we introduce an updated version of SOCOLv3.0, SOCOLv3.1, which includes an improved treatment of ozone sink processes, and results in a reduction in the tropospheric column ozone bias of up to 8DU, mostly due to the inclusion of N2O5 hydrolysis on tropospheric aerosols. As a result of these developments, tropospheric column ozone amounts simulated by SOCOLv3.1 are comparable with several other CCMI models. We apply Gaussian process emulation and sensitivity analysis to understand the remaining ozone bias in SOCOLv3.1. This shows that ozone precursors (nitrogen oxides (NOx), carbon monoxide, methane and other volatile organic compounds, VOCs) are responsible for more than 90% of the variance in tropospheric ozone. However, it may not be the emissions inventories themselves that result in the bias, but how the emissions are handled in SOCOLv3.1, and we discuss this in the wider context of the other CCMI models. Given that the emissions data set to be used for phase 6 of the Coupled Model Intercomparison Project includes approximately 20% more NOx than the data set used for CCMI, further work is urgently needed to address the challenges of simulating sub-grid processes of importance to tropospheric ozone in the current generation of chemistry-climate models. [ABSTRACT FROM AUTHOR]

Energetic particle precipitation (EPP) affects the chemistry of the polar middle atmosphere by producing reactive nitrogen (NOy) and hydrogen (HOx) species, which then catalytically destroy ozone. Recently, there have been major advances in constraining these particle impacts through a parametrization based on high quality observations. Here we investigate the effects of low (auroral) and middle (radiation belt) energy range electrons, separately and in combination, on reactive nitrogen and hydrogen species as well as on ozone during Southern Hemisphere winters from 2002 to 2010 using the chemistry-climate model SOCOL3-MPIOM. Our results show that, in absence of solar proton events, low energy electrons produce the majority of NOy in the polar mesosphere and stratosphere. In the polar vortex, NOy subsides and affects ozone at lower altitudes, down to 10hPa. Comparing a year with high electron precipitation with a quiescent period, we found large ozone depletion in the mesosphere; as the anomaly propagates downward, 15% less ozone is found in the stratosphere during winter, which is confirmed by satellite observations. Only with both low and middle energy electrons, our model reproduces the observed stratospheric ozone anomaly. [ABSTRACT FROM AUTHOR]

Stratospheric sulfate geoengineering (SSG) could contribute to avoiding some of the adverse impacts of climate change. We used the global 3D-aerosol-chemistry-climate model, SOCOL-AER, to investigate 21 different SSG scenarios, each with 1.83MtSyr−1 injected either in the form of accumulation-mode-H2SO4 droplets (AM-H2SO4), gas-phase SO2, or as combinations of both. For most scenarios, the sulfur was continuously emitted at 50hPa (≈20km) altitude in the tropics and subtropics, zonally and latitudinally symmetric about the equator (ranging from ±3.75° to ±30°). In the SO2 emission scenarios, continuous production of tiny nucleation mode particles results in increased coagulation, which together with condensation produces larger coarse mode particles. These larger particles are less effective for backscattering solar radiation and sedimentation out of the stratosphere is faster. On average, AM-H2SO4 injection increases stratospheric aerosol residence times by 32% and stratospheric aerosol burdens 37–41% when comparing to SO2 injection. The modelled all-sky (clear-sky) short-wave radiative forcing for AM-H2SO4 injection scenarios is up to 17–70% (44%–57%) larger than is the case for SO2. Aerosol burdens have a surprisingly week dependence on the latitudinal spread of emissions with emission in the stratospheric surf zone (>15°N–15°S) decreasing burdens by only about 10%. This is because the faster removal through stratosphere-to-troposphere transport via tropopause folds found when injection is spread farther from the equator is roughly balanced by a decrease in coagulation. Increasing injection altitude is also surprisingly ineffective because the increase in burden is compensated by an increase in large aerosols due to increased condensation. Increasing the local SO2 flux in the injection region by pulse or point emissions reduces the total global annual nucleation. Coagulation is also reduced due to the interruption of the continuous flow of freshly formed particles. The net effect of pulse or point emission of SO2 is to increase stratospheric aerosol residence time and radiative forcing. Pulse or point emissions of AM-H2SO4 has the opposite effect—decreasing stratospheric aerosol burden and radiative forcing by increasing coagulation. In summary, this study corroborates previous studies with uncoupled aerosol and radiation modules, suggesting that, compared to SO2 injection, the direct emission of AM-H2SO4 results in more radiative forcing for the same sulfur equivalent mass injection strength and that sensitivities to different injection strategies may vary for different forms of injected sulfur. [ABSTRACT FROM AUTHOR]

Climate models consistently predict an acceleration of the Brewer-Dobson circulation (BDC) due to climate change in the 21st century. However, the strength of this acceleration varies considerably among individual models, which constitutes a notable source of uncertainty for future climate projections. To shed more light upon the magnitude of this uncertainty and on its causes, we analyze the stratospheric mean age of air (AoA) of 10 climate projection simulations from the Chemistry Climate Model Initiative phase 1 (CCMI-I), covering the period between 1960 and 2100. In agreement with previous multi-model studies, we find a large model spread in the magnitude of the AoA trend over the simulation period. Differences between future and past AoA are found to be predominantly due to differences in mixing (reduced aging by mixing and recirculation) rather than differences in residual mean transport. We furthermore analyze the mixing efficiency, a measure of the relative strength of mixing for given residual mean transport, which was previously hypothesized to be a model constant. Here, the mixing efficiency is found to vary not only across models, but also over time in all models. Changes in mixing efficiency are shown to be closely related to changes in AoA and quantified to roughly contribute 10 % to the long-term AoA decrease over the 21st century. Additionally, mixing efficiency variations are shown to considerably enhance model spread in AoA changes. To understand these mixing efficiency variations, we also present a consistent dynamical framework based on diffusive closure, which highlights the role of basic state potential vorticity gradients in controlling mixing efficiency and therefore aging by mixing. [ABSTRACT FROM AUTHOR]

The impact of changes in incoming solar irradiance on stratospheric ozone abundances should be included in climate simulations to aid in capturing the atmospheric response to solar cycle variability. This study presents the first systematic comparison of the representation of the 11-year solar cycle ozone response (SOR) in chemistry--climate models (CCMs) and in pre-calculated ozone databases specified in climate models that do not include chemistry, with a special focus on comparing the recommended protocols for the Coupled Model Intercomparison Project Phase 5 and Phase 6 (CMIP5 and CMIP6). We analyse the SOR in eight CCMs from the Chemistry--Climate Model Initiative (CCMI-1) and compare these with results from three ozone databases for climate models: the Bodeker Scientific ozone database, the SPARC/Atmospheric Chemistry and Climate (AC&C) ozone database for CMIP5 and the SPARC/CCMI ozone database for CMIP6. The peak amplitude of the annual mean SOR in the tropical upper stratosphere (1-5 hPa) decreases by more than a factor of 2, from around 5 to 2 %, between the CMIP5 and CMIP6 ozone databases. This substantial decrease can be traced to the CMIP5 ozone database being constructed from a regression model fit to satellite and ozonesonde measurements, while the CMIP6 database is constructed from CCM simulations. The SOR in the CMIP6 ozone database therefore implicitly resembles the SOR in the CCMI-1 models. The structure in latitude of the SOR in the CMIP6 ozone database and CCMI-1 models is considerably smoother than in the CMIP5 database, which shows unrealistic sharp gradients in the SOR across the middle latitudes owing to the paucity of long-term ozone measurements in polar regions. The SORs in the CMIP6 ozone database and the CCMI-1 models show a seasonal dependence with enhanced meridional gradients at mid- to high latitudes in the winter hemisphere. The CMIP5 ozone database does not account for seasonal variations in the SOR, which is unrealistic. Sensitivity experiments with a global atmospheric model without chemistry (ECHAM6.3) are performed to assess the atmospheric impacts of changes in the representation of the SOR and solar spectral irradiance (SSI) forcing between CMIP5 and CMIP6. The larger amplitude of the SOR in the CMIP5 ozone database compared to CMIP6 causes a likely overestimation of the modelled tropical stratospheric temperature response between 11-year solar cycle minimum and maximum by up to 0.55 K, or around 80% of the total amplitude. This effect is substantially larger than the change in temperature response due to differences in SSI forcing between CMIP5 and CMIP6. The results emphasize the importance of adequately representing the SOR in global models to capture the impact of the 11-year solar cycle on the atmosphere. Since a number of limitations in the representation of the SOR in the CMIP5 ozone database have been identified, we recommend that CMIP6 models without chemistry use the CMIP6 ozone database and the CMIP6 SSI dataset to better capture the climate impacts of solar variability. The SOR coefficients from the CMIP6 ozone database are published with this paper. [ABSTRACT FROM AUTHOR]

We evaluate how the coupled aerosol-chemistry-climate model SOCOL-AERv1.0 represents the influence of the 1991 eruption of Mt. Pinatubo on stratospheric aerosol properties and atmospheric state. The aerosol module is coupled to the radiative and chemical modules and includes comprehensive sulfur chemistry and microphysics, in which the particle size distribution is represented by 40 size bins with radii spanning from 0.39 nm to 3.2 µm. SOCOL-AER simulations are compared with satellite and in situ measurements of aerosol parameters, temperature reanalyses, and ozone observations. In addition to the reference model configuration, we performed series of sensitivity experiments looking at different processes affecting the aerosol layer. An accurate sedimentation scheme is found to be essential to prevent particles from diffusing too rapidly to high and low altitudes. The aerosol radiative feedback and the use of a nudged quasibiennial oscillation help to keep aerosol in the tropics and significantly affect the evolution of the stratospheric aerosol burden, which improves the agreement with observed aerosol mass distributions. The inclusion of van der Waals forces in the particle coagulation scheme suggests improvements in particle effective radius, although other parameters (such as aerosol longevity) deteriorate. Modification of the Pinatubo sulfur emission rate also improves some aerosol parameters, while it worsens others compared to observations. Observations themselves are highly uncertain and render it difficult to conclusively judge the necessity of further model reconfiguration. The model revealed problems in reproducing aerosol sizes above 25 km and also in capturing certain features of the ozone response. Besides this, our results show that SOCOL-AER is capable of predicting the most important global-scale atmospheric effects following volcanic eruptions, which is also a prerequisite for an improved understanding of solar geoengineering effects from sulfur injections to the stratosphere. [ABSTRACT FROM AUTHOR]

Previous multi-model intercomparisons have shown that chemistry-climate models exhibit significant biases in tropospheric ozone compared with observations. We investigate annual-mean tropospheric column ozone in 15 models participating in the SPARC/IGAC (Stratosphere-troposphere Processes and their Role in Climate/International Global Atmospheric Chemistry) Chemistry-Climate Model Initiative (CCMI). These models exhibit a positive bias, on average, of up to 40-50 % in the Northern Hemisphere compared with observations derived from the Ozone Monitoring Instrument and Microwave Limb Sounder (OMI/MLS), and a negative bias of up to ~ 30 % in the Southern Hemisphere. SOCOLv3.0 (version 3 of the Solar-Climate Ozone Links CCM), which participated in CCMI, simulates global-mean tropospheric ozone columns of 40.2 DU - approximately 33 % larger than the CCMI multi-model mean. Here we introduce an updated version of SOCOLv3.0, "SOCOLv3.1, which includes an improved treatment of ozone sink processes, and results in a reduction in the tropospheric column ozone bias of up to 8 DU, mostly due to the inclusion of N2O5 hydrolysis on tropospheric aerosols. As a result of these developments, tropospheric column ozone amounts simulated by SOCOLv3.1 are comparable with several other CCMI models. We apply Gaussian process emulation and sensitivity analysis to understand the remaining ozone bias in SOCOLv3.1. This shows that ozone precursors (nitrogen oxides (NOx, carbon monoxide, methane and other volatile organic compounds) are responsible for more than 90 % of the variance in tropospheric ozone. However, it may not be the emissions inventories themselves that result in the bias, but how the emissions are handled in SOCOLv3.1, and we discuss this in the wider context of the other CCMI models. Given that the emissions data set to be used for phase 6 of the Coupled Model Intercomparison Project includes approximately 20 % more NOx than the data set used for CCMI, further work is urgently needed to address the challenges of simulating sub-grid processes of importance to tropospheric ozone in the current generation of chemistry-climate models. [ABSTRACT FROM AUTHOR]

We have derived values of the Ultraviolet Index (UVI) at solar noon from the Tropospheric Ultraviolet Model (TUV) driven by ozone, temperature and aerosol fields from the first phase of the Chemistry-Climate Model Initiative (CCMI-1). Since clouds remain one of the largest uncertainties in climate projections, we simulated only clear-sky UVI. We compared the UVI climatologies obtained from CCMI and TUV against present-day climatological values of UVI derived from satellite data (the OMI-Aura OMUVBd product) and ground-based measurements (from the NDACC network). Depending on the region, relative differences between the UVI obtained from CCMI and TUV and ground-based measurements ranged between -4 % and 11 %. We calculated the UVI evolution throughout the 21st century for the four Representative Concentration Pathways (RCPs 2.6, 4.5, 6.0 and 8.5). Compared to 1960s values, we found an average increase in UVI in 2100 (of 2-4 %) in the tropical belt (30° N-30° S). For the mid-latitudes, we observed a 1.8 to 3.4 % increase in the Southern Hemisphere for RCP 2.6, 4.5 and 6.0, and found a 2.3 % decrease in RCP 8.5. Higher UV indices are projected in the Northern Hemisphere except for RCP 8.5. At high latitudes, ozone recovery is well identified and induces a complete return of mean UVI levels to 1960 values for RCP 8.5 in the Southern Hemisphere. In the Northern Hemisphere, UVI levels in 2100 are higher by 0.5 to 5.5 % for RCP 2.6, 4.5 and 6.0 and they are lower by 7.9 % for RCP 8.5. We analysed the impacts of greenhouse gases (GHGs) and ozone-depleting substances (ODSs) on UVI from 1960 by comparing CCMI sensitivity simulations (1960-2100) with fixed GHGs or ODSs at their respective 1960 levels. As expected with ODS fixed at their 1960 levels, there is no large decrease in ozone levels and consequently no sudden increase in UVI levels. With fixed GHG, we observed a delayed return of ozone to 1960 values, the same signal is observed on UVI, and looking at the UVI difference between 2090s values and 1960s values, we found an 8 % increase in the tropical belt during the summer of each hemisphere. Finally, we show that, while in the Southern Hemisphere UVI is mainly driven by total ozone column, in the Northern Hemisphere both total ozone column and aerosol optical depth drive UVI levels, with aerosol optical depth having twice as much influence on UVI as total column does. [ABSTRACT FROM AUTHOR]

The stratospheric age of air (AoA) is a useful measure of the overall capabilities of a general circulation model (GCM) to simulate stratospheric transport. Previous studies have reported a large spread in the simulation of AoA by GCMs and coupled chemistry–climate models (CCMs). Compared to observational estimates, simulated AoA is mostly too low. Here we attempt to untangle the processes that lead to the AoA differences between the models and between models and observations. AoA is influenced by both mean transport by the residual circulation and two-way mixing; we quantify the effects of these processes using data from the CCM inter-comparison projects CCMVal-2 (Chemistry–Climate Model Validation Activity 2) and CCMI-1 (Chemistry–Climate Model Initiative, phase 1). Transport along the residual circulation is measured by the residual circulation transit time (RCTT). We interpret the difference between AoA and RCTT as additional aging by mixing. Aging by mixing thus includes mixing on both the resolved and subgrid scale. We find that the spread in AoA between the models is primarily caused by differences in the effects of mixing and only to some extent by differences in residual circulation strength. These effects are quantified by the mixing efficiency, a measure of the relative increase in AoA by mixing. The mixing efficiency varies strongly between the models from 0.24 to 1.02. We show that the mixing efficiency is not only controlled by horizontal mixing, but by vertical mixing and vertical diffusion as well. Possible causes for the differences in the models' mixing efficiencies are discussed. Differences in subgrid-scale mixing (including differences in advection schemes and model resolutions) likely contribute to the differences in mixing efficiency. However, differences in the relative contribution of resolved versus parameterized wave forcing do not appear to be related to differences in mixing efficiency or AoA. [ABSTRACT FROM AUTHOR]

Continued anthropogenic greenhouse gas (GHG) emissions are expected to cause further global warming throughout the 21st century. Understanding the role of natural forcings and their influence on global warming is thus of great interest. Here we investigate the impact of a recently proposed 21st century grand solar minimum on atmospheric chemistry and climate using the SOCOL3-MPIOM chemistry--climate model with an interactive ocean element. We examine five model simulations for the period 2000- 2199, following the greenhouse gas concentration scenario RCP4.5 and a range of different solar forcings. The reference simulation is forced by perpetual repetition of solar cycle 23 until the year 2199. This reference is compared with grand solar minimum simulations, assuming a strong decline in solar activity of 3.5 and 6.5 W m-2, respectively, that last either until 2199 or recover in the 22nd century. Decreased solar activity by 6.5 W m-2 is found to yield up to a doubling of the GHG-induced stratospheric and mesospheric cooling. Under the grand solar minimum scenario, tropospheric temperatures are also projected to decrease compared to the reference. On the global scale a reduced solar forcing compensates for at most 15% of the expected greenhouse warming at the end of the 21st and around 25% at the end of the 22nd century. The regional effects are predicted to be significant, in particular in northern high-latitude winter. In the stratosphere, the reduction of around 15% of incoming ultraviolet radiation leads to a decrease in ozone production by up to 8 %, which overcompensates for the anticipated ozone increase due to reduced stratospheric temperatures and an acceleration of the Brewer--Dobson circulation. This, in turn, leads to a delay in total ozone column recovery from anthropogenic halogen-induced depletion, with a global ozone recovery to the pre-ozone hole values happening only upon completion of the grand solar minimum. [ABSTRACT FROM AUTHOR]

Ozone fields simulated for the first phase of the Chemistry-Climate Model Initiative (CCMI-1) will be used as forcing data in the 6th Coupled Model Intercomparison Project. Here we assess, using reference and sensitivity simulations produced for CCMI-1, the suitability of CCMI-1 model results for this process, investigating the degree of consistency amongst models regarding their responses to variations in individual forcings. We consider the influences of methane, nitrous oxide, a combination of chlorinated or brominated ozone-depleting substances, and a combination of carbon dioxide and other greenhouse gases. We find varying degrees of consistency in the models' responses in ozone to these individual forcings, including some considerable disagreement. In particular, the response of total-column ozone to these forcings is less consistent across the multi-model ensemble than profile comparisons. We analyse how stratospheric age of air, a commonly used diagnostic of stratospheric transport, responds to the forcings. For this diagnostic we find some salient differences in model behaviour, which may explain some of the findings for ozone. The findings imply that the ozone fields derived from CCMI-1 are subject to considerable uncertainties regarding the impacts of these anthropogenic forcings. We offer some thoughts on how to best approach the problem of generating a consensus ozone database from a multi-model ensemble such as CCMI- 1. [ABSTRACT FROM AUTHOR]

We evaluate how the coupled aerosol-chemistry-climate model SOCOL-AER represents the influence of the 1991 eruption of Mt. Pinatubo on stratospheric aerosol loading, aerosol microphysical processes, radiative effects, and atmospheric chemistry. The aerosol module includes comprehensive sulfur chemistry and microphysics, in which the particle size distribution is represented by 40 size bins spanning radii from 0.39nm to 3.2μm. Radiative forcing is computed online using aerosol optical properties calculated according to Mie theory. SOCOL-AER simulations are compared with satellite and in situ measurements of aerosol parameters, temperature reanalyses, and ozone observations. In addition to the reference model configuration, we performed a series of sensitivity experiments looking at different processes affecting the aerosol layer. An accurate sedimentation scheme is found to be essential to prevent particles diffusing too rapidly to high and low altitudes. The aerosol radiative feedback and the use of a nudged quasi-biennial oscillation help to keep aerosol in the tropics and significantly affect the evolution of the stratospheric aerosol burden, which improves the agreement with observed aerosol mass distributions. Changes in the aerosol distribution affected by an inclusion of Van der Waals forces to the particle coagulation scheme suggest improvements in particle effective radius, although other parameters (such as aerosol longevity) deteriorate. Modification of the Pinatubo emission rate also improves some aerosol parameters, while worsens others compared to observations. Observations themselves are highly uncertain and render it difficult to conclusively judge the necessity of further model reconfiguration. In conclusion, our results show that SOCOL-AER is capable of predicting the most important global-scale atmospheric and climate effects following volcanic eruptions, which is also a prerequisite for improved understanding of anthropogenic effects from sulfur emissions. [ABSTRACT FROM AUTHOR]

Stratospheric age of air (AoA) is a useful measure of the overall capabilities of a general circulation model (GCM) to simulate stratospheric transport. Previous studies have reported a large spread in the simulation of AoA by GCMs and coupled chemistry-climate models (CCMs). Compared to observational estimates simulated AoA is mostly too low. Here we attempt to untangle the processes that lead to the AoA differences between the models and between models and observations. AoA is influenced by both, mean transport along the residual circulation and two-way mixing; we quantify the effects of these processes using data from the CCM inter-comparison projects CCMVal-2 and CCMI-1. Transport along the residual circulation is measured by the residual circulation transit time (RCTT). We interpret the difference between AoA and RCTT as additional aging by mixing. Aging by mixing thus includes mixing on both the resolved and subgrid scale. We find that the spread in AoA between the models is primarily caused by differences in the effects of mixing, and only to some extent by differences in residual circulation strength. These effects are quantified by the mixing efficiency, a measure of the relative increase of AoA by mixing. The mixing efficiency varies strongly between the models from 0.21 to 0.99. We show that the mixing efficiency is not only controlled by horizontal mixing, but by vertical mixing and vertical diffusion as well. Possible causes for the differences in the models' mixing efficiencies are discussed. Differences in subgrid scale mixing (including differences in advection schemes and model resolutions) likely contribute to the differences in mixing efficiency. However, differences in the relative contribution of resolved versus parametrized wave forcing do not appear to be related to differences in mixing efficiency or AoA. [ABSTRACT FROM AUTHOR]

Ozone forms in the Earth's atmosphere from the photodissociation of molecular oxygen, primarily in the tropical stratosphere. It is then transported to the extratropics by the Brewer-Dobson circulation (BDC), forming a protective ozone layer around the globe. Human emissions of halogen-containing ozone-depleting substances (hODSs) led to a decline in stratospheric ozone until they were banned by the Montreal Protocol (MP), and since 1998 ozone in the upper stratosphere shows a likely recovery. Total column ozone (TCO) measurements of ozone between the Earth's surface and the top of the atmosphere, indicate that the ozone layer has stopped declining across the globe, but no clear increase has been observed at latitudes outside the polar regions (60-90). Here we report evidence from multiple satellite measurements that ozone in the lower stratosphere between 60° S and 60° N has declined continuously since 1985. We find that, even though upper stratospheric ozone is recovering in response to the MP, the lower stratospheric changes more than compensate for this, resulting in the conclusion that, globally (60° S 60° N), stratospheric column ozone (StCO) continues to deplete. We find that globally, TCO appears to not have decreased because tropospheric column ozone (TrCO) increases, likely the result of human activity and harmful to respiratory health, are compensating for the stratospheric decreases. The reason for the continued reduction of lower stratospheric ozone is not clear, models do not reproduce these trends, and so the causes now urgently need to be established. Reductions in lower stratospheric ozone trends may partly lead to a small reduction in the warming of the climate, but a reduced ozone layer may also permit an increase in harmful ultra-violet (UV) radiation at the surface and would impact human and ecosystem health. [ABSTRACT FROM AUTHOR]

Continued anthropogenic greenhouse gas (GHG) emissions are expected to cause further global warming throughout the 21st century. Understanding potential interferences with natural forcings is thus of great interest. Here we investigate the impact of a recently proposed 21st century grand solar minimum on atmospheric chemistry and climate using the SOCOL3-MPIOM chemistry-climate model with interactive ocean. We examine several model simulations for the period 2000-2199, following the greenhouse gas scenario RCP4.5, but with different solar forcings: the reference simulation is forced by perpetual repetition of solar cycle 23 until the year 2199, whereas the grand solar minimum simulations assume strong declines in solar activity of 3.5 and 6.5 W m-2 with different durations. Decreased solar activity is found to yield up to a doubling of the GHG induced stratospheric and mesospheric cooling. Under the grand solar minimum scenario tropospheric temperatures are also projected to decrease. On the global scale the reduced solar forcing compensates at most 15% of the expected greenhouse warming at the end of 21st and around 25% at the end of 22nd century. The regional effects are predicted to be stronger, in particular in northern high latitude winter. In the stratosphere, the reduced incoming ultraviolet radiation leads to less ozone production by up to 8%, which overcompensates the anticipated ozone increase due to reduced stratospheric temperatures and an acceleration of the Brewer-Dobson circulation. This, in turn, leads to a delay in total ozone column recovery from anthropogenic chlorine-induced depletion, with a global ozone recovery to the pre-ozone hole values happening only upon completion of the grand solar minimum in the 22nd century or later. [ABSTRACT FROM AUTHOR]

To simulate the impacts of volcanic eruptions on the stratosphere, chemistry-climate models that do not include an online aerosol module require temporally and spatially resolved aerosol size parameters for heterogeneous chemistry and aerosol radiative properties as a function of wavelength. For phase 1 of the Chemistry-Climate Model Initiative (CCMI-1) and, later, for phase 6 of the Coupled Model Intercomparison Project (CMIP6) two such stratospheric aerosol data sets were compiled, whose functional capability and representativeness are compared here. For CCMI-1, the SAGE-4λ data set was compiled, which hinges on the measurements at four wavelengths of the SAGE (Stratospheric Aerosol and Gas Experiment) II satellite instrument and uses ground-based Lidar measurements for gap-filling immediately after the Mt. Pinatubo eruption, when the stratosphere was optically opaque for SAGE II. For CMIP6, the new SAGE-3λ data set was compiled, which excludes the least reliable SAGE II wavelength and uses CLAES (Cryogenic Limb Array Etalon Spectrometer) measurements on UARS, the Upper Atmosphere Research Satellite, for gap-filling following the Mt. Pinatubo eruption instead of ground-based Lidars. Here, we performed SOCOLv3 (Solar Climate Ozone Links version 3) chemistry-climate model simulations of the recent past (1986-2005) to investigate the impact of the Mt. Pinatubo eruption in 1991 on stratospheric temperature and ozone and how this response differs depending on which aerosol data set is applied. The use of SAGE-4λ results in heating and ozone loss being overestimated in the lower stratosphere compared to observations in the post-eruption period by approximately 3 K and 0.2 ppmv, respectively. However, less heating occurs in the model simulations based on SAGE-3λ, because the improved gap-filling procedures after the eruption lead to less aerosol loading in the tropical lower stratosphere. As a result, simulated temperature anomalies in the model simulations based on SAGE-3λ for CMIP6 are in excellent agreement with MERRA and ERA-Interim reanalyses in the post-eruption period. Less heating in the simulations with SAGE-3λ means that the rate of tropical upwelling does not strengthen as much as it does in the simulations with SAGE-4λ, which limits dynamical uplift of ozone and therefore provides more time for ozone to accumulate in tropical mid-stratospheric air. Ozone loss following the Mt. Pinatubo eruption is overestimated by 0.1 ppmv in the model simulations based on SAGE-3λ, which is a better agreement with observations than in the simulations based on SAGE-4λ. Overall, the CMIP6 stratospheric aerosol data set, SAGE-3λ, allows SOCOLv3 to more accurately simulate the post-Pinatubo eruption period. [ABSTRACT FROM AUTHOR]

Ozone fields simulated for the Chemistry-Climate Model Initiative (CCMI) will be used as forcing data in the 6th Coupled Model Intercomparison Project (CMIP6). Here we assess, using reference and sensitivity simulations produced for phase 1 of CCMI, the suitability of CCMI-1 model results for this process, investigating the degree of consistency amongst models regarding their responses to variations in individual forcings. We consider the influences of methane, nitrous oxide, a combination of chlorinated or brominated ozone-depleting substances (ODSs), and a combination of carbon dioxide and other greenhouse gases (GHGs). We find varying degrees of consistency in the models' responses in ozone to these individual forcings, including some considerable disagreement. In particular, the response of total-column ozone to these forcings is less consistent across the multi-model ensemble than profile comparisons. The likely cause of this is lower-stratospheric transport and dynamical responses exhibiting substantial inter-model differences. The findings imply that the ozone fields derived from CCMI-1 are subject to considerable uncertainties regarding the impacts of these anthropogenic forcings. [ABSTRACT FROM AUTHOR]

We find that wintertime temperature anomalies near 4 hPa and 50° N/S are related, through dynamics, to anomalies in ozone and temperature, particularly in the tropical stratosphere but also throughout the upper stratosphere and mesosphere. These mid-latitude anomalies occur on timescales of up to a month, and are related to changes in wave forcing. A change in the meridional Brewer-Dobson circulation extends from the middle stratosphere into the mesosphere and forms a temperature-change quadrupole from Equator to pole. We develop a dynamical index based on detrended, deseasonalised mid-latitude temperature. When employed in multiple linear regression, this index can account for up to 60% of the total variability of temperature, peaking at ~5 hPa and dropping to 0 at ~50 and ~0.5 hPa, respectively, and increasing again into the mesosphere. Ozone similarly sees up to an additional 50% of variability accounted for, with a slightly higher maximum and strong altitude dependence, with zero improvement found at 10 hPa. Further, the uncertainty on all equatorial multiple-linear regression coefficients can be reduced by up to 35 and 20% in temperature and ozone, respectively, and so this index is an important tool for quantifying current and future ozone recovery. [ABSTRACT FROM AUTHOR]

Stratospheric water vapour (SWV) is an important component of the Earth's atmosphere as it affects both radiative balance and the chemistry of the atmosphere. Key processes driving changes in SWV include dehydration of air masses transiting the cold-point tropopause (CPT) and methane oxidation. We use a chemistry-climate model to simulate changes in SWV through the 21st century following the four canonical representative concentration pathways (RCPs). Furthermore, we quantify the contribution that methane oxidation makes to SWV following each of the RCPs. Although the methane contribution to SWV maximizes in the upper stratosphere, modelled SWV trends are found to be driven predominantly by warming of the CPT rather than by increasing methane oxidation. SWV changes by -5 to 60% (depending on the location in the atmosphere and emissions scenario) and increases in the lower stratosphere in all RCPs through the 21st century. Because the lower stratosphere is where water vapour radiative forcing maximizes, SWV's influence on surface climate is also expected to increase through the 21st century. [ABSTRACT FROM AUTHOR]

We find that wintertime temperature anomalies near 4 hPa and 50° N/S are related, through dynamics, to anomalies in ozone and temperature, particularly in the tropical stratosphere, but also throughout the upper stratosphere and mesosphere. These mid-latitude anomalies occur on timescales of up to a month, and are related to changes in wave-forcing. A change in the meridional circulation extends from the middle stratosphere into the mesosphere and forms a temperature-change quadrupole from equator to pole. We develop a dynamical index based on detrended, deseasonalised mid-latitude temperature. When employed in multiple linear regression, this index can account for up to 40 % of the total variability of temperature and ozone and a doubling of the total coefficient of determination in the equatorial stratosphere above 20 hPa. Further, the uncertainty on all multiple-linear regression coefficients can be reduced by up to 45 % and 25 % in temperature and ozone, respectively, and so this index is an important tool for quantifying current and future ozone recovery. [ABSTRACT FROM AUTHOR]