ALBERT

Description: High-quality satellite-based measurements are crucial to the assessment of global stratospheric composition change. The Stratospheric Aerosol and Gas Experiment II (SAGE II) provides the longest, continuous data set of vertically resolved ozone and aerosol extinction coefficients to date and therefore remains a cornerstone of understanding and detecting long-term ozone variability and trends in the stratosphere. Despite its stability, SAGE II measurements must be screened for outliers that are a result of excessive aerosol emitted into the atmosphere and that degrade inferences of change. Current methods for SAGE II ozone measurement quality assurance consist of multiple ad hoc and sometimes conflicting rules, leading to too much valuable data being removed or outliers being missed. In this work, the SAGE II ozone data set version 7.00 is used to develop and present a new set of screening recommendations and to compare the output to the screening recommendations currently used. Applying current recommendations to SAGE II ozone leads to unexpected features, such as removing ozone values around zero if the relative error is used as a screening criterion, leading to biases in monthly mean zonal mean ozone concentrations. Most of these current recommendations were developed based on “visual inspection”, leading to inconsistent rules that might not be applicable at every altitude and latitude. Here, a set of new screening recommendations is presented that take into account the knowledge of how the measurements were made. The number of screening recommendations is reduced to three, which mainly remove ozone values that are affected by high aerosol loading and are therefore not reliable measurements. More data remain when applying these new recommendations compared to the rules that are currently being used, leading to more data being available for scientific studies. The SAGE II ozone data set used here is publicly available at https://doi.org/10.5281/zenodo.3710518 (Kremser et al., 2020). The complete SAGE II version 7.00 data set, which includes other variables in addition to ozone, is available at https://eosweb.larc.nasa.gov/project/sage2/sage2_v7_table (last access: December 2019), https://doi.org/10.5067/ERBS/SAGEII/SOLAR_BINARY_L2-V7.0 (SAGE II Science Team, 2012; Damadeo et al., 2013).

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: When computing climatological averages of atmospheric trace-gas mixing ratios obtained from satellite-based measurements, sampling biases arise if data coverage is not uniform in space and time. Homogeneous spatiotemporal coverage is essentially impossible to achieve. Solar occultation measurements, by virtue of satellite orbit and the requirement of direct observation of the sun through the atmosphere, result in particularly sparse spatial coverage. In this proof-of-concept study, a method is presented to adjust for such sampling biases when calculating climatological means. The method is demonstrated using carbonyl sulfide (OCS) measurements at 16 km altitude from the ACE-FTS (Atmospheric Chemistry Experiment Fourier Transform Spectrometer). At this altitude, OCS mixing ratios show a steep gradient between the poles and Equator. ACE-FTS measurements, which are provided as vertically resolved profiles, and integrated stratospheric OCS columns are used in this study. The bias adjustment procedure requires no additional information other than the satellite data product itself. In particular, the method does not rely on atmospheric models with potentially unreliable transport or chemistry parameterizations, and the results can be used uncompromised to test and validate such models. It is expected to be generally applicable when constructing climatologies of long-lived tracers from sparsely and heterogeneously sampled satellite measurements. In the first step of the adjustment procedure, a regression model is used to fit a 2-D surface to all available ACE-FTS OCS measurements as a function of day-of-year and latitude. The regression model fit is used to calculate an adjustment factor that is then used to adjust each measurement individually. The mean of the adjusted measurement points of a chosen latitude range and season is then used as the bias-free climatological value. When applying the adjustment factor to seasonal averages in 30∘ zones, the maximum spatiotemporal sampling bias adjustment was 11 % for OCS mixing ratios at 16 km and 5 % for the stratospheric OCS column. The adjustments were validated against the much denser and more homogeneous OCS data product from the limb-sounding MIPAS (Michelson Interferometer for Passive Atmospheric Sounding) instrument, and both the direction and magnitude of the adjustments were in agreement with the adjustment of the ACE-FTS data.

Description: Upper-air measurements of essential climate variables (ECVs), such as temperature, are crucial for climate monitoring and climate change detection. Because of the internal variability of the climate system, many decades of measurements are typically required to robustly detect any trend in the climate data record. It is imperative for the records to be temporally homogeneous over many decades to confidently estimate any trend. Historically, records of upper-air measurements were primarily made for short-term weather forecasts and as such are seldom suitable for studying long-term climate change as they lack the required continuity and homogeneity. Recognizing this, the Global Climate Observing System (GCOS) Reference Upper-Air Network (GRUAN) has been established to provide reference-quality measurements of climate variables, such as temperature, pressure, and humidity, together with well-characterized and traceable estimates of the measurement uncertainty. To ensure that GRUAN data products are suitable to detect climate change, a scientifically robust instrument replacement strategy must always be adopted whenever there is a change in instrumentation. By fully characterizing any systematic differences between the old and new measurement system a temporally homogeneous data series can be created. One strategy is to operate both the old and new instruments in tandem for some overlap period to characterize any inter-instrument biases. However, this strategy can be prohibitively expensive at measurement sites operated by national weather services or research institutes. An alternative strategy that has been proposed is to alternate between the old and new instruments, so-called interlacing, and then statistically derive the systematic biases between the two instruments. Here we investigate the feasibility of such an approach specifically for radiosondes, i.e. flying the old and new instruments on alternating days. Synthetic data sets are used to explore the applicability of this statistical approach to radiosonde change management.

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

Description: A method, based on climate pattern scaling, has been developed to expand a small number of projections of fields of a selected climate variable (X) into an ensemble that encapsulates a wide range of indicative model structural uncertainties. The method described in this paper is referred to as the Ensemble Projections Incorporating Climate model uncertainty (EPIC) method. Each ensemble member is constructed by adding contributions from (1) a climatology derived from observations that represents the time-invariant part of the signal; (2) a contribution from forced changes in X, where those changes can be statistically related to changes in global mean surface temperature (Tglobal); and (3) a contribution from unforced variability that is generated by a stochastic weather generator. The patterns of unforced variability are also allowed to respond to changes in Tglobal. The statistical relationships between changes in X (and its patterns of variability) and Tglobal are obtained in a training phase. Then, in an implementation phase, 190 simulations of Tglobal are generated using a simple climate model tuned to emulate 19 different global climate models (GCMs) and 10 different carbon cycle models. Using the generated Tglobal time series and the correlation between the forced changes in X and Tglobal, obtained in the training phase, the forced change in the X field can be generated many times using Monte Carlo analysis. A stochastic weather generator is used to generate realistic representations of weather which include spatial coherence. Because GCMs and regional climate models (RCMs) are less likely to correctly represent unforced variability compared to observations, the stochastic weather generator takes as input measures of variability derived from observations, but also responds to forced changes in climate in a way that is consistent with the RCM projections. This approach to generating a large ensemble of projections is many orders of magnitude more computationally efficient than running multiple GCM or RCM simulations. Such a large ensemble of projections permits a description of a probability density function (PDF) of future climate states rather than a small number of individual story lines within that PDF, which may not be representative of the PDF as a whole; the EPIC method largely corrects for such potential sampling biases. The method is useful for providing projections of changes in climate to users wishing to investigate the impacts and implications of climate change in a probabilistic way. A web-based tool, using the EPIC method to provide probabilistic projections of changes in daily maximum and minimum temperatures for New Zealand, has been developed and is described in this paper.

Description: In this study backward trajectories from the tropical lower stratosphere were calculated for the Northern Hemisphere (NH) winters 1995–1996, 1997–1998 (El Niño) and 1998–1999 (La Niña) and summers 1996, 1997 and 1999 using both ERA-40 reanalysis data of the European Centre for Medium-Range Weather Forecast (ECMWF) and coupled Chemistry-Climate Model (CCM) data. The calculated trajectories were analysed to determine the distribution of points where individual air masses encounter the minimum temperature and thus minimum water vapour mixing ratio during their ascent through the tropical tropopause layer (TTL) into the stratosphere. The geographical distribution of these dehydration points and the local conditions there determine the overall water vapour entry into the stratosphere. Results of two CCMs are presented: the ECHAM4.L39(DLR)/CHEM (hereafter: E39/C) from the German Aerospace Center (DLR) and the Freie Universität Berlin Climate Middle Atmosphere Model with interactive chemistry (hereafter: FUB-CMAM-CHEM). In the FUB-CMAM-CHEM model the minimum temperatures are overestimated by about 9 K in NH winter and about 3 K in NH summer, resulting in too high water vapour entry values compared to ERA-40. However, the geographical distribution of dehydration points is fairly similar to ERA-40 for NH winter 1995–1996 and 1998–1999. The distribution of dehydration points in the boreal summer 1996 suggests an influence of the Indian monsoon upon the water vapour transport. The E39/C model displays a temperature bias of about +5 K. Hence, the minimum water vapour mixing ratios are higher relative to ERA-40. The geographical distribution of dehydration points is fairly well in NH winter 1995–1996 and 1997–1998 with respect to ERA-40. The distribution is not reproduced for the NH winter 1998–1999 (La Niña event) compared to ERA-40. There is an excessive water vapour flux through warm regions e.g. Africa in the NH winter and summer. The possible influence of the Indian monsoon on the transport is not seen in the boreal summer 1996. Further, the residence times of air parcels in the TTL were derived from the trajectory calculations. The analysis of the residence times reveals that in both CCMs residence times in the TTL are lower compared to ERA-40 and the seasonal variation is hardly present.

Description: A method, based on climate pattern-scaling, has been developed to expand a small number of projections of fields of a selected climate variable (X) into an ensemble that encapsulates a wide range of model structural uncertainties. The method described in this paper is referred to as the Ensemble Projections Incorporating Climate model uncertainty (EPIC) method. Each ensemble member is constructed by adding contributions from (1) a climatology derived from observations that represents the time invariant part of the signal, (2) a contribution from forced changes in X where those changes can be statistically related to changes in global mean surface temperature (Tglobal), and (3) a contribution from unforced variability that is generated by a stochastic weather generator. The patterns of unforced variability are also allowed to respond to changes in Tglobal. The statistical relationships between changes in X (and its patterns of variability) with Tglobal are obtained in a "training" phase. Then, in an "implementation" phase, 190 simulations of Tglobal are generated using a simple climate model tuned to emulate 19 different Global Climate Models (GCMs) and 10 different carbon cycle models. Using the generated Tglobal time series and the correlation between the forced changes in X and Tglobal, obtained in the "training" phase, the forced change in the X field can be generated many times using Monte Carlo analysis. A stochastic weather generator model is used to generate realistic representations of weather which include spatial coherence. Because GCMs and Regional Climate Models (RCMs) are less likely to correctly represent unforced variability compared to observations, the stochastic weather generator takes as input measures of variability derived from observations, but also responds to forced changes in climate in a way that is consistent with the RCM projections. This approach to generating a large ensemble of projections is many orders of magnitude more computationally efficient than running multiple GCM or RCM simulations. Such a large ensemble of projections permits a description of a Probability Density Function (PDF) of future climate states rather than a small number of individual story lines within that PDF which may not be representative of the PDF as a whole; the EPIC method corrects for such potential sampling biases. The method is useful for providing projections of changes in climate to users wishing to investigate the impacts and implications of climate change in a probabilistic way. A web-based tool, using the EPIC method to provide probabilistic projections of changes in daily maximum and minimum temperatures for New Zealand, has been developed and is described in this paper.