Harmonized Dataset of Ozone Profiles from Satellite Limb and Occultation Measurements

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 similar to 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

Relative drifts and biases between six ozone limb satellite measurements from the last decade

As part of European Space Agency’s (ESA) climate change initiative, high vertical resolution ozone profiles from three instruments all aboard ESA’s Envisat (GOMOS, MIPAS, SCIAMACHY) and ESA’s third party missions (OSIRIS, SMR, ACE-FTS) are to be combined in order to create an essential climate variable data record for the last decade. A prerequisite before combining data is the examination of differences and drifts between the data sets. In this paper, we present a detailed analysis of ozone profile differences based on pairwise collocated measurements, including the evolution of the differences with time. Such a diagnosis is helpful to identify strengths and weaknesses of each data set that may vary in time and introduce uncertainties in long-term trend estimates. The analysis reveals that the relative drift between the sensors is not statistically significant for most pairs of instruments. The relative drift values can be used to estimate the added uncertainty in physical trends. The added drift uncertainty is estimated at about 3% decade$^{-1}$ (1σ). Larger differences and variability in the differences are found in the lowermost stratosphere (below 20 km) and in the mesosphere

Relative Drifts and Biases Between Six Ozone Limb Satellite Measurements From the Last Decade

As part of European Space Agency\u27s (ESA) climate change initiative, high vertical resolution ozone profiles from three instruments all aboard ESA\u27s Envisat (GOMOS, MIPAS, SCIAMACHY) and ESA\u27s third party missions (OSIRIS, SMR, ACE-FTS) are to be combined in order to create an essential climate variable data record for the last decade. A prerequisite before combining data is the examination of differences and drifts between the data sets. In this paper, we present a detailed analysis of ozone profile differences based on pairwise collocated measurements, including the evolution of the differences with time. Such a diagnosis is helpful to identify strengths and weaknesses of each data set that may vary in time and introduce uncertainties in long-term trend estimates. The analysis reveals that the relative drift between the sensors is not statistically significant for most pairs of instruments. The relative drift values can be used to estimate the added uncertainty in physical trends. The added drift uncertainty is estimated at about 3% decade-1 (1σ). Larger differences and variability in the differences are found in the lowermost stratosphere (below 20 km) and in the mesosphere

Validation of water vapour profiles (version 13) retrieved by the IMK/IAA scientific retrieval processor based on full resolution spectra measured by MIPAS on board Envisat

Vertical profiles of stratospheric water vapour measured by the Michelson Interferometer for Passive Atmospheric Sounding (MIPAS) with the full resolution mode between September 2002 and March 2004 and retrieved with the IMK/IAA scientific retrieval processor were compared to a number of independent measurements in order to estimate the bias and to validate the existing precision estimates of the MIPAS data. The estimated precision for MIPAS is 5 to 10% in the stratosphere, depending on altitude, latitude, and season. The independent instruments were: the Halogen Occultation Experiment (HALOE), the Atmospheric Chemistry Experiment Fourier Transform Spectrometer (ACE-FTS), the Improved Limb Atmospheric Spectrometer-II (ILAS-II), the Polar Ozone and Aerosol Measurement (POAM III) instrument, the Middle Atmospheric Water Vapour Radiometer (MIAWARA), the Michelson Interferometer for Passive Atmospheric Sounding, balloon-borne version (MIPAS-B), the Airborne Microwave Stratospheric Observing System (AMSOS), the Fluorescent Stratospheric Hygrometer for Balloon (FLASH-B), the NOAA frostpoint hygrometer, and the Fast In Situ Hygrometer (FISH). For the in-situ measurements and the ground based, air- and balloon borne remote sensing instruments, the measurements are restricted to central and northern Europe. The comparisons to satellite-borne instruments are predominantly at mid- to high latitudes on both hemispheres. In the stratosphere there is no clear indication of a bias in MIPAS data, because the independent measurements in some cases are drier and in some cases are moister than the MIPAS measurements. Compared to the infrared measurements of MIPAS, measurements in the ultraviolet and visible have a tendency to be high, whereas microwave measurements have a tendency to be low. The results of χ<sup>2</sup>-based precision validation are somewhat controversial among the comparison estimates. However, for comparison instruments whose error budget also includes errors due to uncertainties in spectrally interfering species and where good coincidences were found, the χ<sup>2</sup> values found are in the expected range or even below. This suggests that there is no evidence of systematically underestimated MIPAS random errors

Reconciliation of halogen-induced ozone loss with the total-column ozone record

The observed depletion of the ozone layer from the 1980s onwards is attributed to halogen source gases emitted by human activities. However, the precision of this attribution is complicated by year-to-year variations in meteorology, that is, dynamical variability, and by changes in tropospheric ozone concentrations. As such, key aspects of the total-column ozone record, which combines changes in both tropospheric and stratospheric ozone, remain unexplained, such as the apparent absence of a decline in total-column ozone levels before 1980, and of any long-term decline in total-column ozone levels in the tropics. Here we use a chemistry–climate model to estimate changes in halogen-induced ozone loss between 1960 and 2010; the model is constrained by observed meteorology to remove the effects of dynamical variability, and driven by emissions of tropospheric ozone precursors to separate out changes in tropospheric ozone. We show that halogen-induced ozone loss closely followed stratospheric halogen loading over the studied period. Pronounced enhancements in ozone loss were apparent in both hemispheres following the volcanic eruptions of El Chichon and, in particular, Mount Pinatubo, which significantly enhanced stratospheric aerosol loads. We further show that approximately 40% of the long-term non-volcanic ozone loss occurred before 1980, and that long-term ozone loss also occurred in the tropical stratosphere. Finally, we show that halogen-induced ozone loss has declined by over 10% since stratospheric halogen loading peaked in the late 1990s, indicating that the recovery of the ozone layer is well underway

Merged SAGE II, Ozone_cci and OMPS ozone profile dataset and evaluation of ozone trends in the stratosphere

In this paper, we present a merged dataset of ozone profiles from several satellite instruments: SAGE II on ERBS, GOMOS, SCIAMACHY and MIPAS on Envisat, OSIRIS on Odin, ACE-FTS on SCISAT, and OMPS on Suomi-NPP. The merged dataset is created in the framework of the European Space Agency Climate Change Initiative (Ozone_cci) with the aim of analyzing stratospheric ozone trends. For the merged dataset, we used the latest versions of the original ozone datasets. The datasets from the individual instruments have been extensively validated and intercompared; only those datasets which are in good agreement, and do not exhibit significant drifts with respect to collocated ground-based observations and with respect to each other, are used for merging. The long-term SAGE–CCI–OMPS dataset is created by computation and merging of deseasonalized anomalies from individual instruments. The merged SAGE–CCI–OMPS dataset consists of deseasonalized anomalies of ozone in 10° latitude bands from 90° S to 90° N and from 10 to 50 km in steps of 1 km covering the period from October 1984 to July 2016. This newly created dataset is used for evaluating ozone trends in the stratosphere through multiple linear regression. Negative ozone trends in the upper stratosphere are observed before 1997 and positive trends are found after 1997. The upper stratospheric trends are statistically significant at midlatitudes and indicate ozone recovery, as expected from the decrease of stratospheric halogens that started in the middle of the 1990s and stratospheric cooling

HEPPA-II model-measurement intercomparison project : EPP indirect effects during the dynamically perturbed NH winter 2008-2009

We compare simulations from three high-top (with upper lid above 120 km) and five medium-top (with upper lid around 80 km) atmospheric models with observations of odd nitrogen (NOx D NO+NO2), temperature, and carbon monoxide from seven satellite instruments (ACE-FTS on SciSat, GOMOS, MIPAS, and SCIAMACHY on Envisat, MLS on Aura, SABER on TIMED, and SMR on Odin) during the Northern Hemisphere (NH) polar winter 2008/2009. The models included in the comparison are the 3-D chemistry transport model 3dCTM, the ECHAM5/MESSy Atmospheric Chemistry (EMAC) model, FinROSE, the Hamburg Model of the Neutral and Ionized Atmosphere (HAMMO-NIA), the Karlsruhe Simulation Model of the Middle Atmosphere (KASIMA), the modelling tools for SOlar Climate Ozone Links studies (SOCOL and CAO-SOCOL), and the Whole Atmosphere Community Climate Model (WACCM4). The comparison focuses on the energetic particle precipitation (EPP) indirect effect, that is, the polar winter descent of NOx largely produced by EPP in the mesosphere and lower thermosphere. A particular emphasis is given to the impact of the sudden stratospheric warming (SSW) in January 2009 and the subsequent elevated stratopause (ES) event associated with enhanced descent of mesospheric air. The chemistry climate model simulations have been nudged toward reanalysis data in the troposphere and stratosphere while being unconstrained above. An odd nitrogen upper boundary condition obtained from MIPAS observations has further been applied to medium-top models. Most models provide a good representation of the mesospheric tracer descent in general, and the EPP indirect effect in particular, during the unperturbed (pre-SSW) period of the NH winter 2008/2009. The observed NOx descent into the lower mesosphere and stratosphere is generally reproduced within 20 %. Larger discrepancies of a few model simulations could be traced back either to the impact of the models' gravity wave drag scheme on the polar wintertime meridional circulation or to a combination of prescribed NOx mixing ratio at the uppermost model layer and low vertical resolution. In March-April, after the ES event, however, modelled mesospheric and stratospheric NOx distributions deviate significantly from the observations. The too-fast and early downward propagation of the NO x tongue, encountered in most simulations, coincides with a temperature high bias in the lower mesosphere (0.2-0.05 hPa), likely caused by an overestimation of descent velocities. In contrast, upper-mesospheric temperatures (at 0.05-0.001 hPa) are generally underestimated by the high-top models after the onset of the ES event, being indicative for too-slow descent and hence too-low NOx fluxes. As a consequence, the magnitude of the simulated NOx tongue is generally underestimated by these models. Descending NOx amounts simulated with mediumtop models are on average closer to the observations but show a large spread of up to several hundred percent. This is primarily attributed to the different vertical model domains in which the NOx upper boundary condition is applied. In general, the intercomparison demonstrates the ability of state-of- the-art atmospheric models to reproduce the EPP indirect effect in dynamically and geomagnetically quiescent NH winter conditions. The encountered differences between observed and simulated NOx, CO, and temperature distributions during the perturbed phase of the 2009 NH winter, however, emphasize the need for model improvements in the dynamical representation of elevated stratopause events in order to allow for a better description of the EPP indirect effect under these particular conditions.Peer reviewe

HEPPA-II model–measurement intercomparison project: EPP indirect effects during the dynamically perturbed NH winter 2008-2009

We compare simulations from three high-top (with upper lid above 120 km) and five medium-top (with upper lid around 80 km) atmospheric models with observations of odd nitrogen (NOx  =  NO + NO2), temperature, and carbon monoxide from seven satellite instruments (ACE-FTS on SciSat, GOMOS, MIPAS, and SCIAMACHY on Envisat, MLS on Aura, SABER on TIMED, and SMR on Odin) during the Northern Hemisphere (NH) polar winter 2008/2009. The models included in the comparison are the 3-D chemistry transport model 3dCTM, the ECHAM5/MESSy Atmospheric Chemistry (EMAC) model, FinROSE, the Hamburg Model of the Neutral and Ionized Atmosphere (HAMMONIA), the Karlsruhe Simulation Model of the Middle Atmosphere (KASIMA), the modelling tools for SOlar Climate Ozone Links studies (SOCOL and CAO-SOCOL), and the Whole Atmosphere Community Climate Model (WACCM4). The comparison focuses on the energetic particle precipitation (EPP) indirect effect, that is, the polar winter descent of NOx largely produced by EPP in the mesosphere and lower thermosphere. A particular emphasis is given to the impact of the sudden stratospheric warming (SSW) in January 2009 and the subsequent elevated stratopause (ES) event associated with enhanced descent of mesospheric air. The chemistry climate model simulations have been nudged toward reanalysis data in the troposphere and stratosphere while being unconstrained above. An odd nitrogen upper boundary condition obtained from MIPAS observations has further been applied to medium-top models. Most models provide a good representation of the mesospheric tracer descent in general, and the EPP indirect effect in particular, during the unperturbed (pre-SSW) period of the NH winter 2008/2009. The observed NOx descent into the lower mesosphere and stratosphere is generally reproduced within 20 %. Larger discrepancies of a few model simulations could be traced back either to the impact of the models\u27 gravity wave drag scheme on the polar wintertime meridional circulation or to a combination of prescribed NOx mixing ratio at the uppermost model layer and low vertical resolution. In March–April, after the ES event, however, modelled mesospheric and stratospheric NOx distributions deviate significantly from the observations. The too-fast and early downward propagation of the NOx tongue, encountered in most simulations, coincides with a temperature high bias in the lower mesosphere (0.2–0.05 hPa), likely caused by an overestimation of descent velocities. In contrast, upper-mesospheric temperatures (at 0.05–0.001 hPa) are generally underestimated by the high-top models after the onset of the ES event, being indicative for too-slow descent and hence too-low NOx fluxes. As a consequence, the magnitude of the simulated NOx tongue is generally underestimated by these models. Descending NOx amounts simulated with medium-top models are on average closer to the observations but show a large spread of up to several hundred percent. This is primarily attributed to the different vertical model domains in which the NOx upper boundary condition is applied. In general, the intercomparison demonstrates the ability of state-of-the-art atmospheric models to reproduce the EPP indirect effect in dynamically and geomagnetically quiescent NH winter conditions. The encountered differences between observed and simulated NOx, CO, and temperature distributions during the perturbed phase of the 2009 NH winter, however, emphasize the need for model improvements in the dynamical representation of elevated stratopause events in order to allow for a better description of the EPP indirect effect under these particular conditions

Chapter 4: The LOTUS regression model

One of the primary motivations of the LOTUS effort is to attempt to reconcile the discrepancies in ozone trend results from the wealth of literature on the subject. Doing so requires investigating the various methodologies employed to derive long-term trends in ozone as well as to examine the large array of possible variables that feed into those methodologies and analyse their impacts on potential trend results. Given the limited amount of time, the LOTUS group focused on the most common methodology of multiple linear regression and performed a number of sensitivity tests with the goal of trying to establish best practices and come to a consensus on a single regression model to use for this study. This chapter discusses the details and results of the sensitivity tests before describing the components of the final single model that was chosen and the reasons for that choice