ALBERT

Description: The discovery of atmospheric micro(nano)plastic transport and ocean–atmosphere exchange points to a highly complex marine plastic cycle, with negative implications for human and ecosystem health. Yet, observations are currently limited. In this Perspective, we quantify the processes and fluxes of the marine-atmospheric micro(nano)plastic cycle, with the aim of highlighting the remaining unknowns in atmospheric micro(nano)plastic transport. Between 0.013 and 25 million metric tons per year of micro(nano)plastics are potentially being transported within the marine atmosphere and deposited in the oceans. However, the high uncertainty in these marine-atmospheric fluxes is related to data limitations and a lack of study intercomparability. To address the uncertainties and remaining knowledge gaps in the marine-atmospheric micro(nano)plastic cycle, we propose a future global marine-atmospheric micro(nano)plastic observation strategy, incorporating novel sampling methods and the creation of a comparable, harmonized and global data set. Together with long-term observations and intensive investigations, this strategy will help to define the trends in marine-atmospheric pollution and any responses to future policy and management actions.

Description: We quantify the stratospheric injection of brominated very short-lived substances (VSLS) based on aircraft observations acquired in winter 2014 above the Tropical Western Pacific during the CONvective TRansport of Active Species in the Tropics (CONTRAST) and the Airborne Tropical TRopopause EXperiment (ATTREX) campaigns. The overall contribution of VSLS to stratospheric bromine was determined to be 5.0 ± 2.1 ppt, in agreement with the 5 ± 3 ppt estimate provided in the 2014 World Meteorological Organization (WMO) Ozone Assessment report (WMO 2014), but with lower uncertainty. Measurements of organic bromine compounds, including VSLS, were analyzed using CFC-11 as a reference stratospheric tracer. From this analysis, 2.9 ± 0.6 ppt of bromine enters the stratosphere via organic source gas injection of VSLS. This value is two times the mean bromine content of VSLS measured at the tropical tropopause, for regions outside of the Tropical Western Pacific, summarized in WMO 2014. A photochemical box model, constrained to CONTRAST observations, was used to estimate inorganic bromine from measurements of BrO collected by two instruments. The analysis indicates that 2.1 ± 2.1 ppt of bromine enters the stratosphere via inorganic product gas injection. We also examine the representation of brominated VSLS within 14 global models that participated in the Chemistry-Climate Model Initiative. The representation of stratospheric bromine in these models generally lies within the range of our empirical estimate. Models that include explicit representations of VSLS compare better with bromine observations in the lower stratosphere than models that utilize longer-lived chemicals as a surrogate for VSLS. ©2018. American Geophysical Union. All Rights Reserved.

Description: Simulated stratospheric temperatures over the period 1979–2016 in models from the Chemistry-Climate Model Initiative are compared with recently updated and extended satellite data sets. The multimodel mean global temperature trends over 1979–2005 are −0.88 ± 0.23, −0.70 ± 0.16, and −0.50 ± 0.12 K/decade for the Stratospheric Sounding Unit (SSU) channels 3 (~40–50 km), 2 (~35–45 km), and 1 (~25–35 km), respectively (with 95% confidence intervals). These are within the uncertainty bounds of the observed temperature trends from two reprocessed SSU data sets. In the lower stratosphere, the multimodel mean trend in global temperature for the Microwave Sounding Unit channel 4 (~13–22 km) is −0.25 ± 0.12 K/decade over 1979–2005, consistent with observed estimates from three versions of this satellite record. The models and an extended satellite data set comprised of SSU with the Advanced Microwave Sounding Unit-A show weaker global stratospheric cooling over 1998–2016 compared to the period of intensive ozone depletion (1979–1997). This is due to the reduction in ozone-induced cooling from the slowdown of ozone trends and the onset of ozone recovery since the late 1990s. In summary, the results show much better consistency between simulated and satellite-observed stratospheric temperature trends than was reported by Thompson et al. (2012, https://doi.org/10.1038/nature11579) for the previous versions of the SSU record and chemistry-climate models. The improved agreement mainly comes from updates to the satellite records; the range of stratospheric temperature trends over 1979–2005 simulated in Chemistry-Climate Model Initiative models is comparable to the previous generation of chemistry-climate models. ©2018. The Authors.

Description: With low concentrations of tropospheric aerosol, the Southern Ocean offers a “natural laboratory” for studies of aerosol–cloud interactions. Aerosols over the Southern Ocean are produced from biogenic activity in the ocean, which generates sulfate aerosol via dimethylsulfide (DMS) oxidation, and from strong winds and waves that lead to bubble bursting and sea spray emission. Here, we evaluate the representation of Southern Ocean aerosols in the Hadley Centre Global Environmental Model version 3, Global Atmosphere 7.1 (HadGEM3-GA7.1) chemistry–climate model. Compared with aerosol optical depth (AOD) observations from two satellite instruments (the Moderate Resolution Imaging Spectroradiometer, MODIS-Aqua c6.1, and the Multi-angle Imaging Spectroradiometer, MISR), the model simulates too-high AOD during winter and too-low AOD during summer. By switching off DMS emission in the model, we show that sea spray aerosol is the dominant contributor to AOD during winter. In turn, the simulated sea spray aerosol flux depends on near-surface wind speed. By examining MODIS AOD as a function of wind speed from the ERA-Interim reanalysis and comparing it with the model, we show that the sea spray aerosol source function in HadGEM3-GA7.1 overestimates the wind speed dependency. We test a recently developed sea spray aerosol source function derived from measurements made on a Southern Ocean research voyage in 2018. In this source function, the wind speed dependency of the sea spray aerosol flux is less than in the formulation currently implemented in HadGEM3-GA7.1. The new source function leads to good agreement between simulated and observed wintertime AODs over the Southern Ocean; however, it reveals partially compensating errors in DMS-derived AOD. While previous work has tested assumptions regarding the seawater climatology or sea–air flux of DMS, we test the sensitivity of simulated AOD, cloud condensation nuclei and cloud droplet number concentration to three atmospheric sulfate chemistry schemes. The first scheme adds DMS oxidation by halogens and the other two test a recently developed sulfate chemistry scheme for the marine troposphere; one tests gas-phase chemistry only, while the second adds extra aqueous-phase sulfate reactions. We show how simulated sulfur dioxide and sulfuric acid profiles over the Southern Ocean change as a result and how the number concentration and particle size of the soluble Aitken, accumulation and coarse aerosol modes are affected. The new DMS chemistry scheme leads to a 20 % increase in the number concentration of cloud condensation nuclei and cloud droplets, which improves agreement with observations. Our results highlight the importance of atmospheric chemistry for simulating aerosols and clouds accurately over the Southern Ocean.

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: 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 1991 Mt Pinatubo eruption, when the stratosphere was too 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 measurements from CLAES (Cryogenic Limb Array Etalon Spectrometer) 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 tropical 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 tropical 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 up to 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.

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

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

Description: 〉We analyse simulations performed for the Chemistry-Climate Model Initiative (CCMI) to estimate the return dates of the stratospheric ozone layer from depletion caused by anthropogenic stratospheric chlorine and bromine. We consider a total of 155 simulations from 20 models, including a range of sensitivity studies which examine the impact of climate change on ozone recovery. For the control simulations (unconstrained by nudging towards analysed meteorology) there is a large spread (±20 DU in the global average) in the predictions of the absolute ozone column. Therefore, the model results need to be adjusted for biases against historical data. Also, the interannual variability in the model results need to be smoothed in order to provide a reasonably narrow estimate of the range of ozone return dates. Consistent with previous studies, but here for a Representative Concentration Pathway (RCP) of 6.0, these new CCMI simulations project that global total column ozone will return to 1980 values in 2049 (with a 1σ uncertainty of 2043–2055). At Southern Hemisphere mid-latitudes column ozone is projected to return to 1980 values in 2045 (2039–2050), and at Northern Hemisphere mid-latitudes in 2032 (2020–2044). In the polar regions, the return dates are 2060 (2055–2066) in the Antarctic in October and 2034 (2025–2043) in the Arctic in March. The earlier return dates in the Northern Hemisphere reflect the larger sensitivity to dynamical changes. Our estimates of return dates are later than those presented in the 2014 Ozone Assessment by approximately 5–17 years, depending on the region, with the previous best estimates often falling outside of our uncertainty range. In the tropics only around half the models predict a return of ozone to 1980 values, around 2040, while the other half do not reach the 1980 value. All models show a negative trend in tropical total column ozone towards the end of the 21st century. The CCMI models generally agree in their simulation of the time evolution of stratospheric chlorine and bromine, which are the main drivers of ozone loss and recovery. However, there are a few outliers which show that the multi-model mean results for ozone recovery are not as tightly constrained as possible. Throughout the stratosphere the spread of ozone return dates to 1980 values between models tends to correlate with the spread of the return of inorganic chlorine to 1980 values. In the upper stratosphere, greenhouse gas-induced cooling speeds up the return by about 10–20 years. In the lower stratosphere, and for the column, there is a more direct link in the timing of the return dates of ozone and chlorine, especially for the large Antarctic depletion. Comparisons of total column ozone between the models is affected by different predictions of the evolution of tropospheric ozone within the same scenario, presumably due to differing treatment of tropospheric chemistry. Therefore, for many scenarios, clear conclusions can only be drawn for stratospheric ozone columns rather than the total column. As noted by previous studies, the timing of ozone recovery is affected by the evolution of N2O and CH4. However, quantifying the effect in the simulations analysed here is limited by the few realisations available for these experiments compared to internal model variability. The large increase in N2O given in RCP 6.0 extends the ozone return globally by ∼ 15 years relative to N2O fixed at 1960 abundances, mainly because it allows tropical column ozone to be depleted. The effect in extratropical latitudes is much smaller. The large increase in CH4 given in the RCP 8.5 scenario compared to RCP 6.0 also lengthens ozone return by ∼ 15 years, again mainly through its impact in the tropics. Overall, our estimates of ozone return dates are uncertain due to both uncertainties in future scenarios, in particular those of greenhouse gases, and uncertainties in models. The scenario uncertainty is small in the short term but increases with time, and becomes large by the end of the century. There are still some model–model differences related to well-known processes which affect ozone recovery. Efforts need to continue to ensure that models used for assessment purposes accurately represent stratospheric chemistry and the prescribed scenarios of ozone-depleting substances, and only those models are used to calculate return dates. For future assessments of single forcing or combined effects of CO2, CH4, and N2O on the stratospheric column ozone return dates, this work suggests that it is more important to have multi-member (at least three) ensembles for each scenario from every established participating model, rather than a large number of individual models.