ALBERT

Description: Journal of Climate, Ahead of Print. 〈br/〉

Description: 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.

Description: 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.

Description: 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.

Description: 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 analyse 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 analyse the mixing efficiency, a measure of the relative strength of mixing for given residual mean transport, which was previously hypothesised 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.

Description: 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, and since 1998 ozone in the upper stratosphere is rising again, likely the recovery from halogen-induced losses. Total column 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 between 60° S and 60° N 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 indeed continued to decline since 1998. We find that, even though upper stratospheric ozone is recovering, the continuing downward trend in the lower stratosphere prevails, resulting in a downward trend in stratospheric column ozone between 60° S and 60° N. We find that total column ozone between 60° S and 60° N appears not to have decreased only because of increases in tropospheric column ozone that compensate for the stratospheric decreases. The reasons for the continued reduction of lower stratospheric ozone are not clear; models do not reproduce these trends, and thus the causes now urgently need to be established.

Description: 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.

Description: Major mid-winter stratospheric sudden warmings (SSWs) are the largest instance of wintertime variability in the Arctic stratosphere. Because SSWs are able to cause significant surface weather anomalies on intra-seasonal timescales, several previous studies have focused on their potential future change, as might be induced by anthropogenic forcings. However, a wide range of results have been reported, from a future increase in the frequency of SSWs to an actual decrease. Several factors might explain these contradictory results, notably the use of different metrics for the identification of SSWs and the impact of large climatological biases in single-model studies. To bring some clarity, we here revisit the question of future SSW changes, using an identical set of metrics applied consistently across 12 different models participating in the Chemistry–Climate Model Initiative. Our analysis reveals that no statistically significant change in the frequency of SSWs will occur over the 21st century, irrespective of the metric used for the identification of the event. Changes in other SSW characteristics – such as their duration, deceleration of the polar night jet, and the tropospheric forcing – are also assessed: again, we find no evidence of future changes over the 21st century.

Description: 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.

Description: The eruption of Mt. Tambora in 1815 was the largest volcanic eruption of the past 500 years. The eruption had significant climatic impacts, leading to the 1816 year without a summer, and remains a valuable event from which to understand the climatic effects of large stratospheric volcanic sulfur dioxide injections. The eruption also resulted in one of the strongest and most easily identifiable volcanic sulfate signals in polar ice cores, which are widely used to reconstruct the timing and atmospheric sulfate loading of past eruptions. As part of the Model Intercomparison Project on the climatic response to Volcanic forcing (VolMIP), five state-of-the-art global aerosol models simulated this eruption. We analyse both simulated background (no Tambora) and volcanic (with Tambora) sulfate deposition to polar regions and compare to ice core records. The models simulate overall similar patterns of background sulfate deposition, although there are differences in regional details and magnitude. However, the volcanic sulfate deposition varies considerably between the models with differences in timing, spatial pattern and magnitude. Mean simulated deposited sulfate on Antarctica ranges from 19 to 264 kg km−2 and on Greenland from 31 to 194 kg km−2, as compared to the mean ice-core-derived estimates of roughly 50 kg km−2 for both Greenland and Antarctica. The ratio of the hemispheric atmospheric sulfate aerosol burden after the eruption to the average ice sheet deposited sulfate varies between models by up to a factor of 15. Sources of this inter-model variability include differences in both the formation and the transport of sulfate aerosol. Our results suggest that deriving relationships between sulfate deposited on ice sheets and atmospheric sulfate burdens from model simulations may be associated with greater uncertainties than previously thought.