The future ozone trends in changing climate simulated with SOCOLv4

Abstract

This study evaluates the future evolution of atmospheric ozone simulated with the Earth system model (ESM) SOCOLv4. Simulations have been performed based on two potential shared socioeconomic pathways (SSPs): the middle-of-the-road (SSP2-4.5) and fossil-fueled (SSP5-8.5) scenarios. The future changes in ozone, as well as in chemical drivers (NO${}_x$ and CO) and temperature, were estimated between 2015 and 2099 and for several intermediate subperiods (i.e., 2015–2039, 2040–2069, and 2070–2099) via dynamic linear modeling. In both scenarios, the model projects a decline in tropospheric ozone in the future that starts in the 2030s in SSP2-4.5 and after the 2060s in SSP5-8.5 due to a decrease in concentrations of NO${}_x$ and CO. The results also suggest a very likely ozone increase in the mesosphere and the upper and middle stratosphere, as well as in the lower stratosphere at high latitudes. Under SSP5-8.5, the ozone increase in the stratosphere is higher because of stronger cooling ($>\mathrm{1}$ K per decade) induced by greenhouse gases (GHGs), which slows the catalytic ozone destruction cycles. In contrast, in the tropical lower stratosphere, ozone concentrations decrease in both experiments and increase over the middle and high latitudes of both hemispheres due to the speeding up of the Brewer–Dobson circulation, which is stronger in SSP5-8.5. No evidence was found of a decline in ozone levels in the lower stratosphere at mid-latitudes. In both future scenarios, the total column ozone is expected to be distinctly higher than present in middle to high latitudes and might be lower in the tropics, which causes a decrease in the mid-latitudes and an increase in the tropics in the surface level of ultraviolet radiation. The results of SOCOLv4 suggest that the stratospheric-ozone evolution throughout the 21st century is strongly governed not only by a decline in halogen concentration but also by future GHG forcing. In addition, the tropospheric-ozone column changes, which are mainly due to the changes in anthropogenic emissions of ozone precursors, also have a strong impact on the total column. Therefore, even though the anthropogenic halogen-loading problem has been brought under control to date, the sign of future ozone column changes, globally and regionally, is still unclear and largely depends on diverse future human activities. The results of this work are, thus, relevant for developing future strategies for socioeconomic pathways.

1 Introduction

The stratospheric-ozone layer plays an essential role in the Earth’s atmosphere. It shields the ecosystem from dangerous ultraviolet radiation, shapes the vertical temperature profiles, and thus affects the general circulation of the atmosphere. The tropospheric ozone is one of the most potent greenhouse gases (GHGs; e.g., ), contributing to the rise in near-surface temperature, as well as a toxic air pollutant that is harmful to human health and vegetation. Thus, ozone contributes not only to climate change but also to human, agriculture, and ecosystem development (e.g., ).

A serious challenge for humanity is the consequences of stratospheric-ozone depletion caused by man-made halogenated ozone-depleting substances (hODSs). This prompted nations to ratify the Montreal Protocol in 1987, an international treaty to phase out hODSs. The Montreal Protocol and its Amendments and Adjustments (MPA) allows the ozone layer to recover from the hODS effect. Various studies show that the total ozone column decrease has been reversed at most latitudes, which is attributed to the decline in hODS concentrations, highlighting the success of the MPA in protecting the ozone layer . Different projections of the future ozone layer evolution during the 21st century suggest that the decrease in hODSs facilitates ozone recovery in the stratosphere . Ozone abundances are expected to return to the pre-1960 level in most atmospheric layers by the middle to late century, except in the lower stratosphere . Thus, it is believed that declining hODSs will gradually lose their leading role in determining the evolution of the ozone layer throughout the 21st century .

Studies claim that GHGs, such as carbon dioxide (CO${}_{\mathrm{2}}$), methane (CH${}_{\mathrm{4}}$), and nitrous oxide (N${}_{\mathrm{2}}$O), will largely control ozone changes in the 21st century . CO${}_{\mathrm{2}}$ facilitates the stratospheric-ozone enhancement due to direct radiative cooling of the stratosphere, slowing down the gas-phase ozone destruction rate . Therefore, for some parts of the stratosphere, even a super-recovery (i.e., ozone levels well above the pre-1980s level) is expected .

Whilst N${}_{\mathrm{2}}$O is mainly inert in the troposphere, the growth of its concentration will hamper the increase of the stratospheric ozone in the future due to the increased production of nitrogen oxides (NO${}_x$ $=$ NO $+$ NO${}_{\mathrm{2}}$), which catalytically destroy ozone . Yet, the GHG-related cooling of the stratosphere may reduce the efficiency of catalytic cycles involving NO${}_x$. This is due to the fact that more NO${}_x$ is converted to inactive N${}_{\mathrm{2}}$; i.e., the N${}_{\mathrm{2}}$O contribution to ozone destruction can be somewhat lowered .

CH${}_{\mathrm{4}}$ plays an ambivalent role in ozone change, as it may have both negative and positive effects on ozone. The negative effect of increased CH${}_{\mathrm{4}}$ on stratospheric ozone is that it increases the efficiency of the hydroxyl oxide (HO${}_x$) catalytic cycle of ozone destruction, since CH${}_{\mathrm{4}}$ is the main source of H${}_{\mathrm{2}}$O in the middle atmosphere . However, it should be noted that additional HO${}_x$ and NO${}_x$ radicals would also partly compensate for the negative effects of each other in the stratosphere through the production of the reservoir species HNO${}_{\mathrm{3}}$ (OH $+$ NO${}_{\mathrm{2}}$ $+$ $M$ $\to$ HNO${}_{\mathrm{3}}$ $+$ $M$). CH${}_{\mathrm{4}}$ also has a positive effect on ozone, as it causes an additional chlorine deactivation (CH${}_{\mathrm{4}}$ $+$ Cl $\to$ CH${}_{\mathrm{3}}$ $+$ HCl) throughout the stratosphere and promotes an increase in tropospheric ozone by being a source of CO , which is a precursor for ozone formation in the lower atmosphere.

The future evolution of tropospheric ozone will be strongly driven by the changes in CO and NO${}_x$, leading to large differences in projections of tropospheric ozone for distinct climate scenarios . In addition, the projections indicate that the future ozone changes in the troposphere are even more non-linear than in the stratosphere .

Most chemistry–climate models (CCMs) project that the ozone layer will continue to thin in the tropical lower stratosphere throughout the 21st century . The speed of this thinning depends on the climate scenario for GHGs . GHG-induced temperature changes in the lower atmosphere strengthen the meridional transport via the shallow branch of the Brewer–Dobson circulation (BDC) due to an increase in the temperature gradient between tropical and mid-latitudes. This raises the tropopause, alters the wave propagation and dissipation, and extends the subtropical transport barriers upward . The faster atmospheric upwelling decreases the ozone production in the ascending air parcel . The intensified transport also increases the stratosphere–troposphere exchange, with more ozone-poor tropospheric air being transported to the lower stratosphere . Models also exhibit significant differences in the magnitude of the simulated GHG-induced acceleration of the BDC .

Projections of the ozone layer and, hence, of the future surface UV levels strongly depend on the GHG scenarios applied, especially by the end of the 21st century . The current Intergovernmental Panel on Climate Change (IPCC) Coupled Model Intercomparison Project Phase 6 (CMIP6) activities have developed GHG emission scenarios based on shared socioeconomic pathways (SSP), which take economic, demographic, and technological perspectives into account . Therefore, an important task is to examine the sensitivity of the ozone evolution to these contemporary GHG scenarios applied. By analyzing simulations with CMIP6 models under various SSP scenarios,  showed that, under SSP5-8.5, the total ozone column is expected to be 10 DU higher than its 1960 level by the end of the 21st century. On the contrary, total tropical column ozone is not predicted to return to 1960 levels in most of the SSP scenarios due to either tropospheric or lower-stratospheric-ozone decreases . In , an intercomparison of three CMIP6 models under several SSP scenarios (SSP1-2.6, SSP2-4.5, SSP3-7.0, and SSP5-8.5) was presented. The general ozone increase in the global stratosphere has been demonstrated for all employed scenarios. Also, all GHG scenarios contribute positively to closing the Antarctic ozone hole. However, the projected changes in the tropical stratospheric-ozone column are shown to scale non-linearly with the growth of social development, i.e., with incrementing GHG emissions. In addition, showed that, due to the decline in lower-stratospheric ozone, the tropical ozone column is expected to be largely determined by tropospheric-ozone abundance, which might be higher if the SSP5-8.5 scenario plays out. showed the importance of simulating stratospheric ozone accurately for Southern Hemisphere climate change projections, particularly of wind, by comparing CMIP6 model simulations performed with and without interactive chemistry under moderate (SSP2-4.5) and high (SSP5-8.5) scenarios. Their results demonstrate inconsistency between simulations with and without interactive chemistry, showing differences in temperature and westerly wind patterns in the Southern Hemisphere that are driven by differences in Antarctic springtime ozone. This underscores the importance of accurately modeling ozone changes for future climate projections.

Despite the future evolution of atmospheric ozone and its trends on a global and regional scale from various CCMs based on SSP scenarios that have been recently evaluated , the assessment was made without performing a robust statistical or multivariate regression analysis, i.e., excluding the well-known natural forcings to derive future ozone trends. The quantitative analysis of ozone changes can be promoted by state-of-the-art regression models, utilizing a complex and robust statistical approach to diagnose ozone trends. Applying such tools may increase the accuracy of the trend estimation, especially if the tools can handle variables that have a non-linear time-varying change. In the past, it was found that one of the most suitable tools for analyzing ozone evolution and for estimating ozone trends is an advanced type of regression modeling, namely the method of dynamic linear modeling (DLM).

Using DLM to analyze space-borne ozone measurements, provided evidence for an ongoing ozone decrease in the mid-latitude lower stratosphere despite the ozone recovery from the decline in hODSs. CCMs are still incapable of fully reproducing these trends, yet they exhibit some marginally significant signs of ozone decline, which are not completely consistent with observations . The model projections also show no evidence of future lower-stratospheric-ozone decreases at mid-latitudes, whereas they do project the ozone decline in the tropics . This brings into question their ability to accurately simulate future ozone evolution at mid-latitudes, including the most densely populated regions. In essence, asserting the statistical significance and robustness of ozone trends in the lower stratosphere is not straightforward due to the large uncertainties induced by natural variability . Yet, DLM has proven itself to be a flexible regression tool for quantifying highly variable ozone changes and the natural variability contribution to these changes, giving us a reason to expect a higher level of accuracy for trend calculation and estimation of the statistical significance than via conventional multi-linear regression . This motivates an application of DLM to properly evaluate the anthropogenic impact on future atmospheric-ozone trends under modern SSP scenarios.

In this study, we assess future atmospheric-ozone evolution simulated with the SOCOLv4 Earth system model for the period 2015–2099 and for several subperiods (i.e., 2015–2039, 2040–2069, and 2070–2099). To provide the estimates for ozone trends, we carried out two sets of simulations, where the prescribed future GHG evolution and tropospheric-ozone precursors follow either the SSP2-4.5 or SSP5-8.5 scenario. Changes are derived and evaluated by employing the advanced dynamic linear-modeling algorithm . Section 2 outlines the computational methods and experiment design. The results of this study are provided in Sect. 3, followed by the discussion and conclusions summarized in Sect. 4.

2 Computational methods

2.1 The SOCOLv4 ESM description

In this study, simulations were performed with the Earth system model (ESM) SOCOLv4.0 (SOlar Climate Ozone Links, version 4; hereinafter SOCOLv4). SOCOLv4 consists of the Max Planck Institute for Meteorology (MPI-M) ESM version 1.2 (MPI-ESM1.2; ), the chemical module MEZON , and the size-resolving sulfate aerosol microphysical module AER . MPI‐ESM1.2 contains the general circulation model MA-ECHAM6 (the middle-atmosphere version of the European Centre-Hamburg Model version 6) to compute atmospheric transport, physics, and radiation transfer; the Hamburg Ocean Carbon Cycle (HAMOCC) model; the Max Planck Institute for Meteorology Ocean Model (MPIOM); and the Jena Scheme for Biosphere-‐Atmosphere Coupling in Hamburg (JSBACH). A chemical solver is based on the Newton–Raphson implicit iterative method that includes approximately 100 chemical compounds and 216 gas-phase, 72 photochemical, and 16 stratospheric heterogeneous reactions on polar stratospheric cloud particles and in aqueous sulfuric acid aerosols. It is worth saying that updates for MEZON in SOCOLv4, compared to its previous version SOCOLv3 used in CCMI-1, also include several newly discovered and unregulated hODSs, as well as additional chlorine- and bromine-containing very-short-lived substances uncontrolled by the MPA (see ). The advection scheme of operates the transport of chemical species. Photolysis rates are calculated using a lookup-table approach , including the effects of the solar irradiance variability. MA-ECHAM6, MEZON, and AER are interactively coupled, exchanging gas concentrations, sulfate aerosol properties, and meteorological fields.

SOCOLv4 is formulated on the T63 horizontal resolution, which corresponds to $\sim$ 1.9${}^{\circ }$ $\times$ 1.9${}^{\circ }$ and uses 47 vertical levels in hybrid pressure coordinates between Earth's surface and 0.01 hPa ($\sim \mathrm{80}$ km). The 15 min time step is used in SOCOLv4 to calculate dynamic processes, while chemistry and radiation calculations are performed every 2 h. SOCOLv4 reproduces well the distribution of atmospheric tracers, climatology, and the variability of the temperature and circulation fields. Details of the SOCOLv4 model description and validation can be found in .

2.2 Experiment design

Here, we analyze two types of transient simulations spanning the 2015–2099 period and based on projections of GHG emissions from the up-to-date climate scenarios under the shared socioeconomic pathways (SSPs; ). In our study, simulations are performed using two selected SSP scenarios representing pathways of middle-of-the-road (SSP2-4.5) and fossil-fueled (SSP5-8.5) development . Under these scenarios, the surface temperature is expected to rise by about 3 and 5 ${}^{\circ }$C by around 2100, respectively .

The SOCOLv4 simulations are conducted under standard conditions. This means that runs were initiated from MPI-ESM 1.2 restart files for 1970 and that chemistry was initiated from SOCOLv3 runs . This experiment was carried out starting from the year 1949. In 1980, the experiment was divided into ensemble members, which were initialized with slightly changing initial conditions, namely with a small (about 0.1 %) perturbation of the first-month CO${}_{\mathrm{2}}$ concentration. From 2015, all historical climate forcings (following the recommendations of CMIP6; ) are branched to either SSP2-4.5 or SSP5-8.5 scenarios using projected GHG concentrations. The future solar irradiance projection is provided by HEPPA-SOLARIS, as it is also recommended for CMIP6 . Each experiment consists of three ensemble members in order to properly address the impact of internal model variability on ozone evolution and to assess the level of statistical significance of the obtained results. In this study, we analyze trends in the ensemble mean ozone time series, as well as in chemical drivers and temperature.

2.3 Dynamic linear modeling (DLM)

We employ DLM to quantify long-term changes in the variables under study. DLM is a stochastic model to explain the natural or anthropogenic variability in times series using explanatory and proxy variables. Its application for the historical ozone trends and a detailed description can be found in previous studies .

In this study, the DLM setup includes time series of several statistically independent explanatory variables, which attribute to the known natural variability of ozone and which are commonly used for regression analysis of ozone time series . These include the projection of total solar irradiance (TSI – W m${}^{-\mathrm{2}}$; ); the El Niño–Southern Oscillation (ENSO) variability, represented by ENSO's 3.4 index (K) and calculated from the sea surface temperature field; equatorial zonal winds at 30 and 50 hPa, which are two principal components of the quasi-biennial oscillation variability (QBO30 and QBO50, m s${}^{-\mathrm{1}}$); a stratospheric aerosol optical depth (SAOD – dimensionless), which is determined by the aerosol extinction at 300–500 nm band; and the Arctic and Antarctic Oscillation indices (AO and AAO – hPa), calculated from the geopotential height fields at 1000 and 700 mb pressure levels. These proxies are prepared for each ensemble member of both experiments (except for TSI, which is the same for all simulations; see Fig.  in the Appendix). Although future volcanic eruptions were not considered in the simulations, SAOD is also included in the analysis. This was done because aerosol fields are calculated interactively in SOCOLv4 and are slightly different between SSP scenarios (see Fig. ), as they depend in particular on temperature and atmospheric dynamic changes, driven by GHGs. The DLM might be sensitive to these changes.

All used proxies are orthogonal, have admissible covariance, and can be used in the regression analysis (see Fig. ). The DLM also accounts for a first-order autoregressive (AR1) process . In addition, DLM estimates 6- and 12-month harmonics for the seasonal cycle.

The advantage of DLM against conventional multiple linear regressions is that DLM accounts for the level of trend nonlinearity as a free parameter, allowing the trend to evolve over time. This nonlinearity parameter is inferred from the data along with the trend term, seasonal cycle, proxy amplitudes, and the AR1 process . In principle, this makes the DLM method more accurate for capturing the ozone variability, especially for the after-turnaround period (post-1997) of the ozone evolution .

The long-term evolution of the dependent variable, excluding the effects of several independent proxies, is characterized in DLM by the “trend term” or background level. Consequently, we extracted the background level from the DLM output that, in this case, represents the evolution, excluding the effect of natural variability, induced by proxies. We inferred the posterior distributions on the background level by means of the Markov chain Monte Carlo sampling . The DLM was applied for each individual ensemble member of both experiments, using appropriate proxies for each calculation. We have drawn 200 samples from the DLM results, which describe the uncertainty in the posterior distribution. The resulting trends were estimated from the sample mean background levels at each grid point by means of the Mann–Kendall test for the entire 2015–2099 period, as well as for several subperiods, namely 2015–2039, 2040–2069, and 2070–2099, respectively. This was done to properly trace the evolution of trends during the considered period. This is essential, especially in the context of the clarity of ozone change prediction. Following this, the trend estimates from all individual ensemble members are averaged to get the mean trends in the ensemble of each experiment. The statistical significance of the calculated ensemble mean trend is estimated by applying the Student's $t$ test using the standard deviation of trends between individual ensemble members.

3 Results

3.1 Evolution of drivers of ozone change

The temporal evolution of several critical drivers of ozone changes is displayed in Fig.  and demonstrates a considerable difference between the SSP2-4.5 and SSP5-8.5 scenarios. CH${}_{\mathrm{4}}$ starts to decrease in the mid-2040s in SSP2-4.5, whereas it occurs only in the 2070s under SSP5-8.5. Since CO is partially a product of methane, its evolution in the lower atmosphere resembles the change in CH${}_{\mathrm{4}}$ but with a decrease during the first decades in SSP5-8.5. In contrast, near-global NO${}_x$ in SSP5-8.5 increases during this period, and after 2045 it starts to decrease, similarly to the RCP6.0 scenario . Yet, under SSP2-4.5, NO${}_x$ gradually decreased during the entire period. As such, the decline in tropospheric NO${}_x$ and CO columns relates to the air quality change and the decline in CH${}_{\mathrm{4}}$. In addition, NO${}_x$ in the troposphere is produced by lightning activity and airplanes; i.e., future changes in convective activity due to climate change and the growth of aircraft use may contribute to NO${}_x$ production. The resilient increase in CO${}_{\mathrm{2}}$ and N${}_{\mathrm{2}}$O is observed in both scenarios, with a higher and abrupt increase in SSP5-8.5 but with a sharp slowdown in the growth of CO${}_{\mathrm{2}}$ and N${}_{\mathrm{2}}$O concentrations in the last decades of the century, according to SSP2-4.5. Chlorine-containing ODSs (red line in panel e of Fig. ) decrease throughout the whole period. In its turn, a decline in bromine-containing hODSs (dark-red line in panel e of Fig. ) is decelerated by the end of the century. Biogenic isoprene (C${}_{\mathrm{5}}$H${}_{\mathrm{8}}$) evolves with a steady increase in SSP2-4.5, whilst in SSP5-8.5, C${}_{\mathrm{5}}$H${}_{\mathrm{8}}$ decreases till the 2060s and increases slightly by the end of the century. The global mean temperature changes relative to the present time show a stable increase in mean tropospheric and a decrease in mean stratospheric temperatures, with a more intense change under SSP5-8.5 that is in line with expectations discussed in the IPCC report . Under SSP2-4.5, the temperature changes become less pronounced in the late century due to a significant slowdown in CO${}_{\mathrm{2}}$ growth.

Figure 1

Annual mean evolution of drivers of ozone changes between 2015 and 2099 from both SSP2-4.5 (faded lines) and SSP5-8.5 (bold lines) – ODSs are presented as a single line, since their amounts are identical in both considered scenarios. This includes (a) global surface methane (CH${}_{\mathrm{4}}$) concentration [ppbv], (b) global surface nitrous oxide (N${}_{\mathrm{2}}$O) concentration [ppbv], (c) global surface carbon dioxide (CO${}_{\mathrm{2}}$) concentration [ppmv], (d) near-global (60${}^{\circ }$ N–60${}^{\circ }$ S) surface isoprene (C${}_{\mathrm{5}}$H${}_{\mathrm{8}}$) concentration [ppbv], (e) surface total organic chlorine (red line) and bromine ($\times \mathrm{100}$; dark-red line) ODS concentrations [pptv], (f) near-global (60${}^{\circ }$ N–60${}^{\circ }$ S) tropospheric carbon monoxide (CO) column [molec. cm${}^{-\mathrm{2}}$], (g) near-global (60${}^{\circ }$ N–60${}^{\circ }$ S) tropospheric nitrogen oxide (NO${}_x$) column [molec. cm${}^{-\mathrm{2}}$], and (h) global mean changes in averaged tropospheric (red line) and averaged stratospheric (blue line) temperature ($\mathrm{\Delta }T$) regarding its mean for the period 2015–2019.

[Figure omitted. See PDF]

3.2 Ozone anomalies for the period 2015–2099 relative to the present-day ozone concentration

Figure  shows the annual mean partial and total column ozone changes in the near-global region throughout the 21st century with respect to the period 2015–2019 in different atmospheric layers. We calculated changes in O${}_{\mathrm{3}}$ relative to the period 2015–2019 to estimate the future modeled ozone change regarding its current concentration. It was noted that the evolution of tropospheric ozone is largely determined by changes in CH${}_{\mathrm{4}}$, CO, and NO${}_x$ (see Fig. ). The contribution of CO seems to play a larger role in both scenarios, which might be due to less abundant NO${}_x$. Nevertheless, in SSP5-8.5, the sharp decrease in O${}_{\mathrm{3}}$, starting after 2065, resulted mainly from the decrease in NO${}_x$, as both CH${}_{\mathrm{4}}$ and CO start to decrease later. On the other hand, the projected sharp decline in NO${}_x$ makes CO a more important driver of tropospheric-ozone evolution in the last part of the century, especially under SSP5-8.5. Note that the decline in tropospheric-ozone concentration in SSP2-4.5 starts in the 2030s, much earlier than in SSP5-8.5. A steady increase in tropospheric ozone in SSP5-8.5 is observed by the 2060s; afterwards, it starts to sharply decrease during the last decades of the century, similarly to the RCP6.0 scenario . Albeit, in the late century, tropospheric ozone will be lower than it is now in both scenarios – the difference in the zero-crossing point time is about 50 years between scenarios. In SSP2-4.5, the tropospheric-ozone concentration already becomes lower than the present-day one after 2050, while in SSP5-8.5, it is lower only around the end of the century.

Figure 2

Near-global (60${}^{\circ }$ N–60${}^{\circ }$ S) annual mean anomaly ($\mathrm{\Delta }$O${}_{\mathrm{3}}$) of column ozone and DLM fits (both in Dobson units – DU) between 2015 and 2099, presented regarding the O${}_{\mathrm{3}}$ mean for the 2015–2019 period. The red line indicates $\mathrm{\Delta }$O${}_{\mathrm{3}}$ under the SSP2-4.5 scenario; the blue line indicates $\mathrm{\Delta }$O${}_{\mathrm{3}}$ under the SSP5-8.5 scenario. $\mathrm{\Delta }$O${}_{\mathrm{3}}$ is presented for (a) the mesosphere, (b) the upper stratosphere, (c) the middle stratosphere, (d) the lower stratosphere, (e) the entire model atmosphere, (f) the entire stratosphere, and (g) the troposphere. Shadings represent the $\mathrm{1}\mathit{\sigma }$ standard deviation between ensemble members of the experiment.

[Figure omitted. See PDF]

The lower-stratospheric ozone on a near-global scale shows signs of a slight increase until the mid-century in SSP2-4.5. However, over the last half of the century, it began to gradually decline, showing a moderate reduction of about $-\mathrm{1}$ DU (Dobson units) by 2099. In SSP5-8.5, the gradual decrease in ozone is visible during the whole considered period, showing a decrease of about $-\mathrm{4}$ DU by the end of the century. In fact, this ozone decrease is mainly induced by the intensification of transport from the tropics toward the mid-latitudes. In addition, the decline in averaged ozone over 60${}^{\circ }$ N–60${}^{\circ }$ S indicates that the tropical ozone decrease in the lower stratosphere starts to prevail over the ozone recovery from the effects of hODSs on a near-global scale, as seen in Fig. . In contrast, increased NO${}_x$ might still contribute to ozone production in the lower stratosphere via smog reactions (e.g., ), and when it starts to decline, the ozone abundance also decreases more strongly. It should be also mentioned that the expansion of the ozone hole in SOCOLv4 is larger than in observations (see ), and the near-global averaged future ozone decline in the lower stratosphere can be slightly underestimated. In the middle and upper stratosphere, ozone recovers throughout the period due to a decline in hODS level, with a growth of 1 DU (in the middle stratosphere) and 4 DU (in the upper stratosphere) according to SSP2-4.5 and about 4.5 DU (in the middle stratosphere) and 9 DU (in the upper stratosphere) according to SSP5-8.5 by the end of the century. In SSP5-8.5, a much more intense growth after the 2040s is observed when the discrepancy in CO${}_{\mathrm{2}}$ evolution between scenarios becomes larger (see Fig. ); i.e., the stratospheric temperature is lower in SSP5-8.5. In both scenarios, the evolution of the near-global averaged mesospheric ozone also increases. By the end of the century, the ozone content in the mesosphere will be higher than its modern level by $\sim$ 0.12 DU under SSP2-4.5 and by $\sim$ 0.22 DU under SSP5-8.5. This larger ozone enhancement in SSP5-8.5 might be due to lower temperatures and some influence of decreasing NO${}_x$ in the mesosphere. The total ozone column in SSP2-4.5 increases until the mid-century, after which it starts to slowly decline. In SSP5-8.5, the sharp increase in total ozone column transitions into a gradual decrease after the 2060s, which correlates well with the timing of both the sharp decline in tropospheric ozone and the intensification of the ozone decrease in the lower stratosphere. Even so, extra-polar mean total column ozone content by the end of the century will definitely be higher than present, wherein the magnitude of the increase is $\sim$ 2 to 3 times higher in SSP5-8.5 than in SSP2-4.5, which agrees well with previous studies (e.g., ).

3.3 Ozone and driver trends development during the period 2015–2099

The understanding of ozone evolution requires knowing the changes in the driving agents such as temperature and important gas species involved in the ozone production and destruction cycles. The evolution of the CO, NO${}_x$, temperature, and O${}_{\mathrm{3}}$ trends between 2015 and 2099 from both experiments is presented in Fig. .

Figure 3

Profiles of trends in O${}_{\mathrm{3}}$, NO${}_x$, CO, and temperature for the 2015–2099 period and different subperiods from both SSP2-4.5 and SSP5-8.5 simulations. The names of the corresponding scenario and the period are indicated in the upper left corner of each panel. The dashed line is the delimiter of the region, with significance at the 90 % level for positive or negative changes; the solid line is the same at the 95 % level.

[Figure omitted. See PDF]

Carbon monoxide is produced via CO${}_{\mathrm{2}}$ photolysis by solar irradiance in the upper atmosphere (e.g., ) and can be transported downward, mostly over the high latitudes during cold seasons. Therefore, its abundance in these areas strongly reflects CO${}_{\mathrm{2}}$ behavior, mimicking a steady increase in SSP5-8.5 and stabilization in SSP2-4.5. The increase in CO in the stratosphere should not strongly contribute to the ozone changes; however, some slight effect can be expected from the removal of OH caused by the CO $+$ OH $\to$ CO${}_{\mathrm{2}}$ $+$ H reaction . In the troposphere, the CO source is driven by methane and biogenic volatile organic compounds (VOCs). Therefore, we observe a steady CO decline after 2040 in SSP2-4.5 following the drop in methane emissions (see Fig. ). For the SSP5-8.5, the change of signs appears in 2070 after flattening and a small decline of the methane mixing ratio. An initial negative tendency for the 2015–2039 subperiod is related to a small decrease in VOCs. The CO tendencies in the stratosphere are defined by the upward transport and mixing of the tropospheric air. Carbon monoxide can be considered as a proxy for the level of organic species, which are a necessary part of the tropospheric-ozone production mechanism. The concentration of NO${}_x$ is the second part participating in this process.

In the mesosphere, NO${}_x$ (NO $+$ NO${}_{\mathrm{2}}$) is mostly produced by N${}_{\mathrm{2}}$O oxidation, energetic particles, and influx from the thermosphere. They can be destroyed by solar irradiance via NO${}_{\mathrm{2}}$ photolysis followed by a cannibalistic N $+$ NO $\to$ N${}_{\mathrm{2}}$ $+$ O reaction. Because the thermospheric source is the same for both cases and is partly accounted for by solar proxies, the NO${}_x$ trend in the mesosphere depends on the available N${}_{\mathrm{2}}$O and temperature, which regulate the efficiency of the cannibalistic reaction, making it faster for the cooler environment in the future. Despite a steady N${}_{\mathrm{2}}$O increase (see Fig. ), the N${}_{\mathrm{2}}$O in the mesosphere is less available due to its higher destruction by enhanced ozone and O(${}^{\mathrm{1}}$D) concentration in the stratosphere. Thus, less N${}_{\mathrm{2}}$O abundance and cooler temperature lead to a general decrease of the mesospheric NO${}_x$. For the 2070–2099 subperiod, however, the NO${}_x$ depletion for the SSP2-4.5 case is not so pronounced due to mesospheric cooling that is probably very small.

Stratospheric NO${}_x$ concentration is mostly regulated by the production via N${}_{\mathrm{2}}$O $+$ O(${}^{\mathrm{1}}$D) $\to$ NO $+$ NO and the conversion to reservoir species, which depend on the temperature and availability of hydrogen- and halogen-containing species, which deactivate NO${}_x$-building reservoir species like HNO${}_{\mathrm{3}}$ or ClONO${}_{\mathrm{2}}$. Therefore, the stratospheric NO${}_x$ increase is more substantial in the SSP2-4.5 case when the cooling and water vapor increases are not so pronounced as in the SSP5-8.5 case (see for trends in H${}_{\mathrm{2}}$O).

NO${}_x$ mostly declines in the lower troposphere due to improved air quality. Also, most periods in both scenarios show the permanent increase in free tropospheric NO${}_x$ over the Northern Hemisphere's upper troposphere, which is maintained by aircraft emissions.

The temperature trend patterns look as expected (e.g., ). A continuous increase of greenhouse gases leads to tropospheric warming and stratospheric cooling (e.g., ), and both are substantially more pronounced in the SSP5-8.5 scenario due to more intensive anthropogenic activity. The tropospheric warming in this case is more prominent over the Northern Hemisphere due to Arctic amplification (e.g., ) and in the lower stratosphere over the southern high latitude, where the ozone concentration is increasing due to the recovery from the halogen loading to pre-ozone hole conditions. The radiative cooling by greenhouse gases dominates in the stratosphere in relation to some warming caused by the stratospheric-ozone increase, and it agrees with the time evolution shown in Fig. . During the first 2015–2039 subperiod, the quadrupole structure of stratospheric temperature trends is observed in both scenarios, which is dynamically induced and is barely observed in later subperiods.

The ozone change patterns substantially differ between layers. In the troposphere, the ozone decrease is observed for the SSP2-4.5 scenario starting from 2040, as well as for the entire period. This behavior is explained by the continuous decrease of the ozone precursors related to the improvement of air quality. For the SSP5-8.5, a similar process occurs only after 2070 when the NO${}_x$ atmospheric abundance decline, together with the decline in CO, is the most prominent (see Fig. g). A similar decrease in the tropospheric ozone resembles the results obtained by  using the RCP6.0 scenario. Some increase in NO${}_x$ level before 2070 leads to positive tropospheric-ozone trends, which makes the ozone trend positive for the entire period. The pattern and magnitude of obtained statistically significant tropospheric-ozone trends over the entire period are consistent with those in the multi-model mean given in  and for WACCM and IPSL models from  for corresponding SSP scenarios.

In the upper stratosphere and southern lower stratosphere, the ozone increase is very persistent because it is driven by a steady decline of the halogen loading (see Fig. ). The ozone increase in the upper stratosphere is stronger for the SSP5-8.5 case because a more pronounced stratospheric cooling leads to less intensive catalytic ozone destruction cycles. Another area with a persistent trend appears in the tropical lower stratosphere, where intensification in the warmer-climate Brewer–Dobson circulation drives negative ozone trends (e.g., ). This feature is more pronounced for the SSP5-8.5 scenario after 2070 because of the stronger warming. Before 2070 and for the entire period, the magnitude of the ozone decline in this area is virtually the same for both cases due to compensation of the dynamical loss by means of increased tropospheric ozone obtained for SSP5-8.5. Overall, stratospheric-ozone trends are mostly statistically significant for the entire period and are consistent with previous findings . In MRI-ESM2, the pattern of future ozone trends (see ) differs from that modeled with SOCOLv4, but this was anticipated due to limitations identified in MRI-ESM2 .

In the upper mesosphere, ozone decreases until 2070 under SSP5-8.5 due to an increase in CH${}_{\mathrm{4}}$ causing an increase in mesospheric abundance of H${}_{\mathrm{2}}$O and, hence, an enhancement of HO${}_x$ radicals. Under SSP2-4.5, mesospheric ozone has generally increased over the entire period of 2015–2099, since CH${}_{\mathrm{4}}$ only increases slightly until the 2040s and then begins to decline.

3.4 Total column ozone trend development during the period 2015–2099

One way or another, changes in ozone in different layers of the atmosphere contribute to a change in total column ozone. It is essential for humanity to know the future evolution of total ozone because it affects changes in ground-level UV radiation. Figure  shows the evolution of trends in total column ozone as a function of month and latitude over the period 2015–2099 and intermediate subperiods.

Figure 4

Trends in total column ozone as a function of month and latitude for the 2015–2099 period and different subperiods from both SSP2-4.5 and SSP5-8.5 simulations. The names of the corresponding scenario and the period are indicated in the upper left corner of each panel. The dashed line is the delimiter of the region, with significance at the 90 % level for positive or negative changes; the solid line is the same at the 95 % level.

[Figure omitted. See PDF]

The total ozone recovery in austral spring over the Southern Hemisphere is generally similar between both scenarios, since it is driven by the phase-out of hODSs emissions, which are identical in both scenarios. However, the ozone increase is slightly higher in SSP5-8.5 owing to a lower temperature in the stratosphere. It is also seen that, during 2070–2099, in both scenarios, the ozone increase is slowed down. This might be because of the slower hODS decline and the GHG increase (see Fig. ) and due to the contribution of tropospheric-ozone decline (see Fig. ). In mid-latitudes, the ozone increase is also higher in SSP5-8.5 due to both temperature and more intense transport from the tropics. In contrast, in the tropics, trends in total ozone largely differ between scenarios. In SSP2-4.5, the tropical total ozone tends to reduce during the entire period by about $-\mathrm{2}$ DU per decade due to ozone decreases in the lower stratosphere and troposphere. A strong decline in total ozone of about $-\mathrm{4}$ DU per decade between 2070 and 2099 is observed in SSP5-8.5, but in other subperiods, the trend in tropical total column ozone is generally near zero due to an increase in tropospheric ozone that partly compensates for the ozone decline in the lower-stratospheric ozone. During boreal spring, the total ozone also increases in the Northern Hemisphere in both scenarios, with a higher increase in SSP5-8.5. The presented statistically significant total ozone column change distributions for the entire 2015–2099 period are highly compatible with the multi-model mean given in  for both considered SSP scenarios.

Thus, some increase in surface UV level over the tropics and a decrease over the middle and high latitudes can be expected in both scenarios, but in SSP2-4.5, it will be higher in the tropics throughout the entire period, while in SSP5-8.5, it will be higher in the late century only. On the contrary, the decrease in surface UV level at middle and high latitudes is expected to be greater in SSP5-8.5 than in SSP2-4.5 due to higher ozone.

4 Discussion and conclusions

In this paper, we have evaluated atmospheric-ozone trends based on two sets of ensemble simulations using SOCOLv4, covering the period from 2015 to 2099. One simulation is based on the SSP2-4.5 scenario, and the other is based on the SSP5-8.5 scenario; these scenarios differ in terms of greenhouse gas emissions and ozone precursors. The trends in ozone, as well as in non-hODS drivers of ozone evolution such as NO${}_x$, CO, and temperature, are derived using DLM. The ozone layer is expected to increase on a near-global scale throughout the entire century because of the ban on the production of hODSs by the Montreal Protocol. However, the evolution of atmospheric ozone in different atmospheric layers differs greatly between the two SSP scenarios. The tropospheric-ozone evolution, driven mainly by CO and NO${}_x$ changes, shows a difference in the time of inflection point when tropospheric ozone begins to decline. In SSP2-4.5, it began to be observed after 2040, while in SSP5-8.5, it began to be observed after 2065. In the mesosphere and the upper and middle stratosphere, a resilient increase in ozone is about 2 to 3 times higher in SSP5-8.5 than in SSP2-4.5, since the negative temperature trends in these regions in SSP5-8.5 are more than 1 K per decade stronger, which retards the catalytic ozone loss. In the lower stratosphere, the near-global ozone content tends to decline after 2040 in SSP2-4.5 and during the whole considered period in SSP5-8.5; by the end of the century, the decrease in SSP5-8.5 is more than 3 times higher due to faster meridional transport of ozone to the poles. The obtained ozone trends for different regions are consistent with those presented in previous studies .

In general, it is difficult to establish trends in ozone in the extratropical lower stratosphere due to the large uncertainty associated with natural variability . However, for the long-term periods, statistically robust projection in this part of the atmosphere is possible, as we show in this study. In addition, there are other factors which might contribute to the uncertainty of future ozone evolution. For instance, we have no future volcanic activity considered in our study because it is hard to predict volcanic eruptions. However, severe implications for the ozone layer in the future are expected if strong volcanic eruptions occur . In addition, no less important for the future ozone evolution might be the projected decline in solar activity throughout the 21st century , which also has not been considered in our study. Yet, it is well known that solar activity mainly drives photochemical and dynamical processes in the stratosphere and is responsible for the ozone formation and radiation budget . Therefore, a decline in solar activity might lead to a decrease in atmospheric-ozone production, causing some negative implications for its future evolution .

Nevertheless, total column ozone is expected to increase almost everywhere, except in the tropics. In both polar regions, the total ozone increases with a slightly higher intensity in SSP5-8.5. In the mid-latitudes, the total ozone also increases thanks to the upper-stratospheric-ozone increase and transport from the tropics. Conversely, in the tropics, it generally declines in SSP2-4.5 due to both tropospheric- and lower-stratospheric-ozone decreases, equating to about $-\mathrm{2}$ DU per decade; it changes in SSP5-8.5, with a sharp decrease of about $-\mathrm{4}$ DU per decade only during the last decades of the century due to a severe reduction in both tropospheric- and lower-stratospheric-ozone content. We showed that, besides changes in the stratospheric-ozone column, it is also essential to consider the tropospheric column ozone evolution, since it may seriously contribute to total column ozone evolution, especially in the tropics.

A much stronger ozone increase in the upper part of the middle atmosphere and the middle to high latitudes of the lower stratosphere might also be expected in the SSP5-8.5 scenario. In this regard, it may seem that the more-greenhouse-gases scenario is better because, despite higher near-surface temperatures, it will be more favorable for ozone increases over the most populated areas. However, the excessive increase in ozone over middle to high latitudes may also have negative consequences for human well-being. Exceeding the required level of total ozone content, especially over the most inhabited areas, means more UV absorption and, consequently, less surface-level UV radiation than what is required for human health. It causes less vitamin D synthesis and therefore increases the risk of diseases related to vitamin D deficiency, like rickets and osteomalacia . In addition, it is worth paying attention to the evolution of ozone in the tropics. There is a risk of a decrease in total ozone content, leading to an increase in surface UV level abnormalities, which also has negative effects on human health, like an increased risk of skin cancer and cataracts .

The important message in this regard is to find a way to bring the ozone content in the atmosphere to an equilibrium state, where it is neither lower nor higher than necessary. Thus, we emphasize that the findings presented in this study will be useful for further improvement of socioeconomic-pathway policies to determine the route to maintain the global total ozone content that is favorable for the sustainable development of human civilization.

Appendix A

Figure A1

Input quantities (proxy variables) for the forcing of the SOCOLv4 simulations. Faded colors are for SSP2-4.5; bright colors are for SSP 5-8.5.

[Figure omitted. See PDF]

Figure A2

Pearson’s linear correlation coefficients for different covariate variables (proxy variables) of the SOCOLv4 simulations for the 2015–2099 period: (a) for SSP2-4.5 and (b) for SSP5-8.5. Error bars represent the $\mathrm{1}\mathit{\sigma }$ standard deviation of the correlation coefficients between ensemble members.

[Figure omitted. See PDF]

Data availability

The SOCOLv4 simulations of future ozone evolution based on the SSP2-4.5 and SSP5-8.5 emission scenarios can be accessed from 10.5281/zenodo.7318315 .

Author contributions

AKD processed the data, visualized the results, and prepared the original draft of the paper. ER and TP supervised this research. ER originated the idea for this study. TS, JS, and AKD designed the experiments and performed simulations. TE refined the draft of the paper and contributed to the analysis of the results. All the authors participated in editing the paper and discussing the results.

Competing interests

The contact author has declared that none of the authors has any competing interests.

Disclaimer

Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Acknowledgements

Arseniy Karagodin-Doyennel, Eugene Rozanov, Timofei Sukhodolov, Tatiana Egorova, and Jan Sedlacek are grateful to the Swiss National Science Foundation for supporting this research through the no. 200020-182239 project POLE (Polar Ozone Layer Evolution). The work of Eugene Rozanov and Timofei Sukhodolov has been partly performed in the SPbSU “Ozone Layer and Upper Atmosphere Research” laboratory, supported by the Ministry of Science and Higher Education of the Russian Federation under agreement no. 075-15-2021-583. Calculations were supported by a grant from the Swiss National Supercomputing Centre (CSCS) under projects S-901 (ID no. 154), S-1029 (ID no. 249), and S-903.

Financial support

This research has been supported by the Schweizerischer Nationalfonds zur Förderung der Wissenschaftlichen Forschung (project POLE (Polar Ozone Layer Evolution; grant no. 200020-182239)) and the Ministry of Science and Higher Education of the Russian Federation (grant no. 075-15-2021-583). Publisher's note: Copernicus Publications has not received any payments from Russian or Belarusian institutions for this paper.

Review statement

This paper was edited by Jens-Uwe Grooß and reviewed by two anonymous referees.