An extensive observational data set, consisting of more than 10⁶ SF₆ vertical profiles from MIPAS measurements distributed over the whole globe has been condensed into monthly zonal means of mean age of air for the period September 2002 to January 2010, binned at 10° latitude and 1–2 km altitude. The data were analysed with respect to their temporal variation by fitting a regression model consisting of a constant and a linear increase term, 2 proxies for the QBO variation, sinusoidal terms for the seasonal and semi-annual variation and overtones for the correction of the shapes to the observed data set. The impact of subsidence of mesospheric SF₆-depleted air and in-mixing into non-polar latitudes on mid-latitudinal absolute age of air and its linear increase was assessed and found to be small. The linear increase of mean age of stratospheric air was found to be positive and partly larger than the trend derived by Engel et al. (2009) for most of the Northern mid-latitudes, the middle stratosphere in the tropics, and parts of the Southern mid-latitudes, as well as for the Southern polar upper stratosphere. Multi-year decrease of age of air was found for the lowermost and the upper stratospheric tropics, for parts of Southern mid-latitudes, and for the Northern polar regions. Analysis of the amplitudes and phases of the seasonal variation shed light on the coupling of stratospheric regions to each other. In particular, the Northern mid-latitude stratosphere is well coupled to the tropics, while the Northern lowermost mid-latitudinal stratosphere is decoupled, confirming the separation of the shallow branch of the Brewer-Dobson circulation from the deep branch. We suggest an overall increased tropical upwelling, together with weakening of mixing barriers, especially in the Northern Hemisphere, as a hypothetical model to explain the observed pattern of linear multi-year increase/decrease, and amplitudes and phase shifts of the seasonal variation.

Drifts, trends and periodic variations were calculated from monthly zonally averaged ozone profiles. The ozone profiles were derived from level-1b data of the Michelson Interferometer for Passive Atmospheric Sounding (MIPAS) by means of the scientific level-2 processor run by the Karlsruhe Institute of Technology (KIT), Institute for Meteorology and Climate Research (IMK). All trend and drift analyses were performed using a multilinear parametric trend model which includes a linear term, several harmonics with period lengths from 3 to 24 months and the quasi-biennial oscillation (QBO). Drifts at 2-sigma significance level were mainly negative for ozone relative to Aura MLS and Odin OSIRIS and negative or near zero for most of the comparisons to lidar measurements. Lidar stations used here include those at Hohenpeissenberg (47.8° N, 11.0° E), Lauder (45.0° S, 169.7° E), Mauna Loa (19.5° N, 155.6° W), Observatoire Haute Provence (43.9° N, 5.7° E) and Table Mountain (34.4° N, 117.7° W). Drifts against the Atmospheric Chemistry Experiment Fourier Transform Spectrometer (ACE-FTS) were found to be mostly insignificant. The assessed MIPAS ozone trends cover the time period of July 2002 to April 2012 and range from −0.56 ppmv decade⁻¹ to +0.48 ppmv decade⁻¹ (−0.52 ppmv decade⁻¹ to +0.47 ppmv decade⁻¹ when displayed on pressure coordinates) depending on altitude/pressure and latitude. From the empirical drift analyses we conclude that the real ozone trends might be slightly more positive/less negative than those calculated from the MIPAS data, by conceding the possibility of MIPAS having a very small (approximately within −0.3 ppmv decade⁻¹) negative drift for ozone. This leads to drift-corrected trends of −0.41 ppmv decade⁻¹ to +0.55 ppmv decade⁻¹ (−0.38 ppmv decade⁻¹ to +0.53 ppmv decade⁻¹ when displayed on pressure coordinates) for the time period covered by MIPAS Envisat measurements, with very few negative and large areas of positive trends at mid-latitudes for both hemispheres around and above 30 km (~10 hPa). Negative trends are found in the tropics around 25 and 35 km (~25 and 5 hPa), while an area of positive trends is located right above the tropical tropopause. These findings are in good agreement with the recent literature. Differences of the trends compared with the recent literature could be explained by a possible shift of the subtropical mixing barriers. Results for the altitude–latitude distribution of amplitudes of the quasi-biennial, annual and the semi-annual oscillation are overall in very good agreement with recent findings.

Vertically resolved distributions of sulfur dioxide (SO₂) with global coverage in the height region from the upper troposphere to ~20 km altitude have been derived from observations by the Michelson Interferometer for Passive Atmospheric Sounding (MIPAS) on Envisat for the period July 2002 to April 2012. Retrieved volume mixing ratio profiles representing single measurements are characterized by typical errors in the range of 70–100 pptv and by a vertical resolution ranging from 3 to 5 km. Comparison with observations by the Atmospheric Chemistry Experiment Fourier transform spectrometer (ACE-FTS) revealed a slightly varying bias with altitude of −20 to 50 pptv for the MIPAS data set in case of volcanically enhanced concentrations. For background concentrations the comparison showed a systematic difference between the two major MIPAS observation periods. After debiasing, the difference could be reduced to biases within −10 to 20 pptv in the altitude range of 10–20 km with respect to ACE-FTS. Further comparisons of the debiased MIPAS data set with in situ measurements from various aircraft campaigns showed no obvious inconsistencies within a range of around ±50 pptv. The SO₂ emissions of more than 30 volcanic eruptions could be identified in the upper troposphere and lower stratosphere (UTLS). Emitted SO₂ masses and lifetimes within different altitude ranges in the UTLS have been derived for a large part of these eruptions. Masses are in most cases within estimations derived from other instruments. From three of the major eruptions within the MIPAS measurement period – Kasatochi in August 2008, Sarychev in June 2009 and Nabro in June 2011 – derived lifetimes of SO₂ for the altitude ranges 10–14, 14–18 and 18–22 km are 13.3 ± 2.1, 23.6 ± 1.2 and 32.3 ± 5.5 days respectively. By omitting periods with obvious volcanic influence we have derived background mixing ratio distributions of SO₂. At 10 km altitude these indicate an annual cycle at northern mid- and high latitudes with maximum values in summer and an amplitude of about 30 pptv. At higher altitudes of about 16–18 km, enhanced mixing ratios of SO₂ can be found in the regions of the Asian and the North American monsoons in summer – a possible connection to an aerosol layer discovered by Vernier et al. (2011b) in that region.

In this paper, we present a HARMonized dataset of OZone profiles (HARMOZ) based on limb and occultation measurements from Envisat (GOMOS, MIPAS and SCIAMACHY), Odin (OSIRIS, SMR) and SCISAT (ACE-FTS) satellite instruments. These measurements provide high-vertical-resolution ozone profiles covering the altitude range from the upper troposphere up to the mesosphere in years 2001–2012. HARMOZ has been created in the framework of the European Space Agency Climate Change Initiative project. The harmonized dataset consists of original retrieved ozone profiles from each instrument, which are screened for invalid data by the instrument teams. While the original ozone profiles are presented in different units and on different vertical grids, the harmonized dataset is given on a common pressure grid in netCDF (network common data form)-4 format. The pressure grid corresponds to vertical sampling of ~ 1 km below 20 km and 2–3 km above 20 km. The vertical range of the ozone profiles is specific for each instrument, thus all information contained in the original data is preserved. Provided altitude and temperature profiles allow the representation of ozone profiles in number density or mixing ratio on a pressure or altitude vertical grid. Geolocation, uncertainty estimates and vertical resolution are provided for each profile. For each instrument, optional parameters, which are related to the data quality, are also included. For convenience of users, tables of biases between each pair of instruments for each month, as well as bias uncertainties, are provided. These tables characterize the data consistency and can be used in various bias and drift analyses, which are needed, for instance, for combining several datasets to obtain a long-term climate dataset. This user-friendly dataset can be interesting and useful for various analyses and applications, such as data merging, data validation, assimilation and scientific research. The dataset is available at http://www.esa-ozone-cci.org/?q=node/161 or at doi:10.5270/esa-ozone_cci-limb_occultation_profiles-2001_2012-v_1-201308.

This study investigates the abrupt and severe water vapour decline in the stratosphere beginning in the year 2000 (the "millennium water vapour drop") and other similarly strong stratospheric water vapour reductions by means of various simulations with the state-of-the-art Chemistry-Climate Model (CCM) EMAC (ECHAM/MESSy Atmospheric Chemistry Model). The model simulations differ with respect to the prescribed sea surface temperatures (SSTs) and whether nudging is applied or not. The CCM EMAC is able to most closely reproduce the signature and pattern of the water vapour drop in agreement with those derived from satellite observations if the model is nudged. Model results confirm that this extraordinary water vapour decline is particularly obvious in the tropical lower stratosphere and is related to a large decrease in cold point temperature. The drop signal propagates under dilution to the higher stratosphere and to the poles via the Brewer–Dobson circulation (BDC). We found that the driving forces for this significant decline in water vapour mixing ratios are tropical sea surface temperature (SST) changes due to a coincidence with a preceding strong El Niño–Southern Oscillation event (1997/1998) followed by a strong La Niña event (1999/2000) and supported by the change of the westerly to the easterly phase of the equatorial stratospheric quasi-biennial oscillation (QBO) in 2000. Correct (observed) SSTs are important for triggering the strong decline in water vapour. There are indications that, at least partly, SSTs contribute to the long period of low water vapour values from 2001 to 2006. For this period, the specific dynamical state of the atmosphere (overall atmospheric large-scale wind and temperature distribution) is important as well, as it causes the observed persistent low cold point temperatures. These are induced by a period of increased upwelling, which, however, has no corresponding pronounced signature in SSTs anomalies in the tropics. Our free-running simulations do not capture the drop as observed, because a) the cold point temperature has a low bias and thus the water vapour variability is reduced and b) because they do not simulate the appropriate dynamical state. Large negative water vapour declines are also found in other years and seem to be a feature which can be found after strong combined El Niño/La Niña events if the QBO west phase during La Niña changes to the east phase.