Abstract. This study presents an updated assessment of stratospheric and total column ozone trends over the 2000–2024 period using six merged limb-profile datasets and six merged total ozone datasets. Long-term changes were quantified using a multiple linear regression framework that accounts for dynamical and chemical variability. In addition to standard regressors (solar cycle, QBO, ENSO, stratospheric aerosol optical depth), we include Arctic and Antarctic Oscillation indices and the eddy heat flux in each hemisphere as proxies for dynamic variability. Volcanic (and wildfire) aerosol forcing is represented by separate proxies for three periods dominated by the major volcanic events of El Chichón, Pinatubo, and post-2000 volcanic eruptions, including Hunga-Tonga. These period-specific proxies are employed to better account for varying dynamical ozone responses that largely depend on the season and location of the eruptions. All profile datasets consistently show positive trends in the upper stratosphere, with the strongest ozone recovery in southern mid-latitudes, in agreement with other studies. In the lower stratosphere, trends remain weak, spatially heterogeneous, and predominantly negative. A comparison of stratospheric column trends derived from profile data with total ozone trends shows close agreement across latitude bands. Within the trend uncertainties, total column trends since 2000 are largely driven by stratospheric ozone changes, while tropospheric contributions to zonal-mean total ozone trends (the difference between total and stratospheric column trends) appear negligible. The extended regression framework improves the representation of recent dynamical variability and provides an updated perspective on stratospheric ozone recovery through 2024.

Abstract. A variety of chemical and dynamical processes in the troposphere and stratosphere affect tropical total column ozone (TCO), the net effect of which may cause changes in surface UV radiation and impact human and ecosystem health. We use dynamical linear modeling to estimate tropical trends in TCO and partial column ozone (PCO) in the troposphere and three stratospheric layers to assess agreement between satellite observational composites and chemistry–climate model (CCM) simulations from two multi-model experiments (CCMI-1 and CCMI-2022). While both model experiments show tropical TCO increases over 2000–2021, CCMI–2022 trends (+2.5 DU) agree slightly better with observations than CCMI-1 (+1.6 DU). However, this overall agreement obscures multiple systematic differences in PCO trends between the models and observations across atmospheric layers. For example, since 2000 tropical tropospheric PCO increased significantly in CCMI-2022 (+1.5 DU) but not in CCMI–1 (+0.3 DU), largely explaining the difference in TCO trends. Also, despite nearly identical stratospheric PCO trends, CCMI-2022 trends are slightly more negative in the lower stratospheric (by ~0.5 DU), compensated by more positive middle/upper stratosphere trends compared to CCMI-1. Crucially, substantial differences exist across observational PCO trends, particularly in the troposphere and middle/upper stratosphere, and these disagreements limit the ability to evaluate CCM fidelity. Furthermore, while the inter-model correlation between late and early 21st century trends is suggestive of a potential emergent constraint on future ozone trends, the spread in observational trends precludes its observational implementation.

In recent years, the need for high-quality long-term mesospheric ozone records has become increasingly evident, as they are essential for understanding chemical, dynamical, and radiative processes in the middle and upper atmosphere and their coupling with the lower layers. Here, we present a new merged dataset of ozone profiles in the middle atmosphere (METEOR-O3), created from several limb-viewing satellite instruments: HALOE, GOMOS, MIPAS, ACE-FTS, MLS, and SOFIE. The merged dataset covers the period from 1991 to 2023 and provides deseasonalized ozone anomalies in 10 degrees latitude bins between 80 degrees S and 80 degrees N, from approximately 22 to 100 km. The deseasonalized ozone anomalies are used for global and seasonal trend analysis. The results show positive upper stratospheric ozone trends in both hemispheres, with magnitudes of 1 %-2 % per decade between 35 and 45 km, indicating continued ozone recovery consistent with previous assessments. In contrast, mesospheric ozone (above similar to 60 km) exhibits negative trends of -1 % to -3 % per decade, with the strongest decreases of about -8 % to -12 % per decade between 80 and 90 km. Seasonal analyses confirm positive trends in the upper stratosphere across all seasons and persistent negative trends in the upper mesosphere, strongest at high latitudes above 80 km. The METEOR-O3 dataset provides the first global, long-term merged record suitable for detailed studies of mesospheric/lower thermospheric ozone variability and trend evaluation, providing valuable information for model validation and assessments of upper atmospheric changes.

Anthropogenic greenhouse gas emissions significantly impact the middle and upper atmosphere. They cause cooling and thermal shrinking and affect the atmospheric structure. Atmospheric contraction results in changes in key atmospheric features, such as the stratopause height or the peak ionospheric electron density, and also results in reduced thermosphere density. These changes can impact, among others, the lifespan of objects in low Earth orbit, refraction of radio communication and GPS signals, and the peak altitudes of meteoroids entering the Earth's atmosphere. Given this, there is a critical need for observational capabilities to monitor the middle and upper atmosphere. Equally important is the commitment to maintaining and improving long‐term, homogeneous data collection. However, capabilities to observe the middle and upper atmosphere are decreasing rather than improving.

This paper presents an intercomparison between existing tropospheric ozone column (TrOC) datasets obtained using combined limb and nadir observations, i.e., exploiting collocated stratospheric profile and total column information retrieved from limb and nadir satellite observations, respectively. In particular, seven datasets have been considered, covering the past 2 decades and consisting of monthly-averaged time series with nearly global coverage. We perform a comparison in terms of climatology and seasonality, investigate the tropopause height used for the construction of each dataset and the related biases, and finally discuss long-term TrOC drifts and trends. The overall goal of the study is to assess the consistency between the datasets and explore possible strategies to reconcile the differences between them. Despite uncertainties associated with the limb–nadir residual methodology and large biases between the mean values of the considered datasets, we identify an overall agreement of TrOC distribution patterns. The different tropopause height definitions used to construct the datasets did not show a relevant role in explaining the biases between them. We demonstrate that a thorough investigation of the drifts with respect to ground-based observations is needed to evaluate TrOC trends from satellite data and that long-term trends in specific regions can be consistently detected, e.g., a positive trend of up to 1.5 DU per decade over China for the 2005–2021 period.

This work uses multispectral measurements of vertically resolved aerosol extinction coefficients from the Stratospheric Aerosol and Gas Experiment (SAGE) III on the International Space Station (ISS) to demonstrate how the use of the Ångström exponent for interpolation of aerosol data between two different wavelengths creates a bias. An empirical relationship is derived between the magnitude of this bias and the Ångström exponent at several different SAGE wavelengths. This relationship can thus be used as a correction factor for other studies, such as multi-instrument intercomparisons or merging, that wish to convert aerosol data from one wavelength to another using the Ångström exponent and is applicable to all stratospheric non-cloud aerosol except highly aged particles that are evaporating at altitudes above the Junge layer.

In this paper, we present the new multi-wavelength dataset of aerosol extinction profiles, which are retrieved from the averaged transmittance spectra by the Global Ozone Monitoring by Occultation of Stars instrument aboard the Envisat satellite. Using monthly and zonally averaged transmittances as a starting point for the retrievals enables us to improve the signal-to-noise ratio and eliminate possible modulation of transmittance spectra by uncorrected scintillations. The two-step retrieval method is used: the spectral inversion is followed by the vertical inversion. The spectral inversion relies on the removal of contributions from ozone, NO2, NO3 and Rayleigh scattering from the optical depth spectra for each ray perigee altitude. In the vertical inversion, the profiles of aerosol extinction coefficients at several wavelengths are retrieved from the collection of slant aerosol optical depth profiles. The retrieved aerosol extinction profiles (FMI-GOMOSaero dataset v1) are provided in the altitude range 10–40 km at wavelengths of 400, 440, 452, 470, 500, 525, 550, 672 and 750 nm for the whole GOMOS operating period from August 2002 to March 2012. Extensive intercomparisons of the retrieved FMI-GOMOSaero aerosol profiles with aerosol profile data from other satellite instruments at several wavelengths have been performed. It is found that the average difference between FMI-GOMOSaero and other datasets is within 20 %–40 % in the lower and middle stratosphere, the standard deviation is ∼ 20 %–50 %, and the correlation coefficient of the time series is 0.65–0.85. The created FMI-GOMOSaero dataset can be used in merged datasets of stratospheric aerosols. It might be also used as a priori information for satellite retrievals during 2002–2012.

Climate-related studies need information about the distribution of stratospheric aerosols, which influence the energy balance of the Earth's atmosphere. In this work, we present a merged dataset of vertically resolved stratospheric aerosol extinction coefficients, which is derived using data from six limb and occultation satellite instruments: SAGE (Stratospheric Aerosol and Gas Experiment) II on ERBS (Earth Radiation Budget Satellite), GOMOS (Global Ozone Monitoring by Occultation of Stars) and SCIAMACHY (Scanning Imaging Spectrometer for Atmospheric Chartography) on Envisat, OSIRIS (Optical Spectrograph and InfraRed Imaging System) on Odin, OMPS (Ozone Monitor Profiling Suite Limb Profiler) on Suomi NPP, and SAGE III on the ISS (International Space Station). The merging of aerosol profiles is performed via the transformation of the aerosol datasets from individual satellite instruments to the same wavelength (750 nm) and their de-biasing and homogenization by adjusting the seasonal cycles. After such homogenization, the data from individual satellite instruments are in good agreement. The merged aerosol extinction coefficient is computed as the median of the adjusted data from the individual instruments. The merged time series of vertically resolved monthly mean aerosol extinction coefficients at 750 nm is provided in 10 degrees latitudinal bins from 90 degrees S to 90 degrees N, in the altitude range from 8.5 to 39.5 km. The time series of the stratospheric aerosol optical depth (SAOD) is created via the integration of aerosol extinction profiles from the tropopause to 39.5 km; it is also provided as monthly mean data in 10 degrees latitudinal bins. The created aerosol climate record covers the period from October 1984 until December 2023, and it is intended to be extended in the future. The merged CREST aerosol dataset (v2) is available at 10.57707/fmib2share.dfe14351fd8548bcaca3c2956b17f665 (Sofieva et al., 2024a). It can be used in various climate-related studies.

In this paper, we present the updated SAGE-CCI-OMPS+ climate data record of monthly zonal mean ozone profiles. This dataset covers the stratosphere and combines measurements by nine limb and occultation satellite instruments – SAGE II (Stratospheric Aerosol and Gases Experiment II), OSIRIS (Optical Spectrograph and InfraRed Imaging System), MIPAS (Michelson Interferometer for Passive Atmospheric Sounding), SCIAMACHY (SCanning Imaging Spectrometer for Atmospheric CHartographY), GOMOS (Global Ozone Monitoring by Occultation of Stars), ACE-FTS (Atmospheric Chemistry Experiment Fourier Transform Spectrometer), OMPS-LP (Ozone Monitor Profiling Suite Limb Profiler), POAM (Polar Ozone and Aerosol Measurement) III, and SAGE III/ISS (Stratospheric Aerosol and Gases Experiment III on the International Space Station). Compared to the original version of the SAGE-CCI-OMPS dataset (Sofieva et al., 2017b), the update includes new versions of MIPAS, ACE-FTS, and OSIRIS datasets and introduces data from additional sensors (POAM III and SAGE III/ISS) and retrieval processors (OMPS-LP). In this paper, we show detailed intercomparisons of ozone profiles from different instruments and data versions, with a focus on the detection of possible drifts in the datasets. The SAGE-CCI-OMPS+ dataset has a better coverage of polar regions and of the upper troposphere and the lower stratosphere (UTLS) than the previous dataset. We also studied the influence of including new datasets on ozone trends, which are estimated using multiple linear regression. The changes in the merged dataset do not change the overall morphology of post-1997 ozone trends; statistically significant trends are observed in the upper stratosphere. The largest changes in ozone trends are observed in polar regions, especially in the Southern Hemisphere. The updated SAGE-CCI-OMPS+ dataset contains profiles of deseasonalized anomalies and ozone concentrations from 1984 to 2021, in 10∘ latitude bins from 90∘ S to 90∘ N and in the altitude range from 10 to 50 km. The dataset is open access and available at https://climate.esa.int/en/projects/ozone/data/ (last access: 9 March 2023) and at ftp://cci_web@ftp-ae.oma.be/esacci (ESA Climate Office; last access: 9 March 2023).

<p>Total ozone is a measure of the protection of the biosphere from UV radiation. Extratropical total ozone recovery trends of about +0.5%/decade are consistent with the continuous decline in stratospheric halogen loading since the middle 1990s as a consequence of the Montreal Protocol and its Amendments on phasing out ODS. Nevertheless, the recovery (or chemistry-related) trends in the northern hemisphere (NH) are compensated by changes in atmospheric dynamics leading to stable ozone levels since about 2000 (apart from year-to-year variability). The near-global (60°S-60°N) average total ozone level is currently about 2-3% below the 1964-1980 mean.</p> <p>Lower stratospheric ozone, the dominating contributor to the total ozone column and derived from limb satellite observations, shows globally (below 60° latitude) a slightly negative (but mostly statistically insignificant) trend of about -1.5%/decade since 2000, which is not consistent with the stable total ozone levels in the NH.  Some studies suggested that increases in tropospheric ozone compensate for negative lower stratospheric trends leading to stable column levels. The trend regression models applied to ozone profiles and total ozone also differ. For instance, proxies representing changes in atmospheric dynamics and transport have been included for total ozone, but not for ozone profiles, which may also lead to inconsistencies between total and stratospheric column trends.</p> <p>This presentation will report updated total ozone trends from five merged total ozone datasets up to 2022. The same regression model, including proxies representing atmospheric dynamics and transport, will be then applied to ozone profiles from the merged limb dataset (SAGE II-SCIAMACHY-OMPS) to evaluate the consistency between column and profile trends.</p> <p><!-- Deuteranopia filter --> <!-- Tritanopia filter --></p>

3.5/3.6(panel (a)) and version 4.1/4.2(panel (b)).

Satellite measurements in nadir and limb viewing geometry provide a complementary view of the atmosphere. An effective combination of the limb and nadir measurements can give new information about atmospheric composition. In this work, we present tropospheric ozone column datasets that have been created using a combination of total ozone columns from OMI (Ozone Monitoring Instrument) and TROPOMI (TROPOspheric Monitoring Instrument) with stratospheric ozone column datasets from several available limb-viewing instruments: MLS (Microwave Limb Sounder), OSIRIS (Optical Spectrograph and InfraRed Imaging System), MIPAS (Michelson Interferometer for Passive Atmospheric Sounding), SCIAMACHY (SCanning Imaging Spectrometer for Atmospheric CHartographY), OMPS-LP (Ozone Mapping and Profiles Suite – Limb Profiler), and GOMOS (Global Ozone Monitoring by Occultation of Stars). We have developed further the methodological aspects of the assessment of tropospheric ozone using the residual method supported by simulations with the chemistry transport model SILAM (System for Integrated modeLling of Atmospheric coMposition). It has been shown that the accurate assessment of ozone in the upper troposphere and the lower stratosphere (UTLS) is of high importance for detecting the ground-level ozone patterns. The stratospheric ozone column is derived from a combination of ozone profiles from several satellite instruments in limb-viewing geometry. We developed a method for the data homogenization, which includes the removal of biases and a posteriori estimation of random uncertainties, thus making the data from different instruments compatible with each other. The high-horizontal- and vertical-resolution dataset of ozone profiles is created via interpolation of the limb profiles from each day to a 1∘×1∘ horizonal grid. A new kriging-type interpolation method, which takes into account data uncertainties and the information about natural ozone variations from the SILAM-adjusted ozone field, has been developed. To mitigate the limited accuracy and coverage of the limb profile data in the UTLS, a smooth transition to the model data is applied below the tropopause. This allows for the estimation of the stratospheric ozone column with full coverage of the UTLS. The derived ozone profiles are in very good agreement with collocated ozonesonde measurements. The residual method was successfully applied to OMI and TROPOMI clear-sky total ozone data in combination with the stratospheric ozone column from the developed high-resolution limb profile dataset. The resulting tropospheric ozone column is in very good agreement with other satellite data. The global distributions of tropospheric ozone exhibit enhancements associated with the regions of high tropospheric ozone production. The main datasets created are (i) a monthly 1∘×1∘ global tropospheric ozone column dataset (from ground to 3 km below the tropopause) using OMI and limb instruments, (ii) a monthly 1∘×1∘ global tropospheric ozone column dataset using TROPOMI and limb instruments, and (iii) a daily 1∘×1∘ interpolated stratospheric ozone column from limb instruments. Other datasets, which are created as an intermediate step of creating the tropospheric ozone column data, are (i) a daily 1∘×1∘ clear-sky and total ozone column from OMI and TROPOMI, (ii) a daily 1∘×1∘ homogenized and interpolated dataset of ozone profiles from limb instruments, and (iii) a daily 1∘×1∘ dataset of ozone profiles from SILAM simulations with adjustment to satellite data. These datasets can be used in various studies related to variability and trends in ozone distributions in both the troposphere and the stratosphere. The datasets are processed from the beginning of OMI and TROPOMI measurements until December 2020 and are planned to be regularly extended in the future.

Caption: Lightning discharges appear in various colours depending on the scatter of light inside the thundercloud and in the atmosphere.The intracloud lightning discharges in the centre of the thundercloud appear to be white with a bluish tint, and the cloud-to-ground discharge below appears to be orange.The right hand side of the thundercloud exhibits a green tint that is attributed to the unique composition of hydrometeors inside the thundercloud.The photo was taken in the late evening of 10 September 2013, near Tarragona in northeastern Spain.

Validation of the Integrated Water Vapour (IWV) from Sentinel-3 Ocean and Land Colour Instrument (OLCI) was performed as a part of the “ESA/Copernicus Space Component Validation for Land Surface Temperature, Aerosol Optical Depth and Water Vapour Sentinel-3 Products” (LAW) project. High-spatial-resolution IWV observations in the near-infrared spectral region from the OLCI instruments aboard the Sentinel-3A and Sentinel-3B satellites provide continuity with observations from MERIS (Medium Resolution Imaging Spectrometer). The IWV was compared with reference observations from two networks: GNSS (Global Navigation Satellite System) precipitable water vapour from the SuomiNet network and integrated lower tropospheric columns from radio-soundings from the IGRA (Integrated Radiosonde Archive) database. Results for cloud-free matchups over land show a wet bias of 7 %–10 % for OLCI, with a high correlation against the reference observations (0.98 against SuomiNet and 0.90 against IGRA). Both OLCI-A and OLCI-B instruments show almost identical results, apart from an anomaly observed in camera 3 of the OLCI-B instrument, where observed biases are lower than in other cameras in either instrument. The wavelength drift in sensors was investigated, and biases in different cameras were found to be independent of wavelength. Effect of cloud proximity was found to have almost no effect on observed biases, indicating that cloud flagging in the OLCI IWV product is sufficiently reliable. We performed validation of random uncertainty estimates and found them to be consistent with the statistical a posteriori estimates, but somewhat higher.

This study presents an updated evaluation of stratospheric ozone profile trends in the 60∘ S–60∘ N latitude range over the 2000–2020 period using an updated version of the Long-term Ozone Trends and Uncertainties in the Stratosphere (LOTUS) regression model that was used to evaluate such trends up to 2016 for the last WMO Ozone Assessment (2018). In addition to the derivation of detailed trends as a function of latitude and vertical coordinates, the regressions are performed with the datasets averaged over broad latitude bands, i.e. 60–35∘ S, 20∘ S–20∘ N and 35–60∘ N. The same methodology as in the last assessment is applied to combine trends in these broad latitude bands in order to compare the results with the previous studies. Longitudinally resolved merged satellite records are also considered in order to provide a better comparison with trends retrieved from ground-based records, e.g. lidar, ozonesondes, Umkehr, microwave and Fourier transform infrared (FTIR) spectrometers at selected stations where long-term time series are available. The study includes a comparison with trends derived from the REF-C2 simulations of the Chemistry Climate Model Initiative (CCMI-1). This work confirms past results showing an ozone increase in the upper stratosphere, which is now significant in the three broad latitude bands. The increase is largest in the Northern and Southern Hemisphere midlatitudes, with ∼2.2 ± 0.7 % per decade at ∼2.1 hPa and ∼2.1 ± 0.6 % per decade at ∼3.2 hPa respectively compared to ∼1.6 ± 0.6 % per decade at ∼2.6 hPa in the tropics. New trend signals have emerged from the records, such as a significant decrease in ozone in the tropics around 35 hPa and a non-significant increase in ozone in the southern midlatitudes at about 20 hPa. Non-significant negative ozone trends are derived in the lowermost stratosphere, with the most pronounced trends in the tropics. While a very good agreement is obtained between trends from merged satellite records and the CCMI-1 REF-C2 simulation in the upper stratosphere, observed negative trends in the lower stratosphere are not reproduced by models at southern and, in particular, at northern midlatitudes, where models report an ozone increase. However, the lower-stratospheric trend uncertainties are quite large, for both measured and modelled trends. Finally, 2000–2020 stratospheric ozone trends derived from the ground-based and longitudinally resolved satellite records are in reasonable agreement over the European Alpine and tropical regions, while at the Lauder station in the Southern Hemisphere midlatitudes they show some differences.

In this paper, we present the MErged GRIdded Dataset of Ozone Profiles (MEGRIDOP) in the stratosphere with a resolved longitudinal structure, which is derived from data from six limb and occultation satellite instruments: GOMOS, SCIAMACHY and MIPAS on Envisat, OSIRIS on Odin, OMPS on Suomi-NPP, and MLS on Aura. The merged dataset was generated as a contribution to the European Space Agency Climate Change Initiative Ozone project (Ozone_cci). The period of this merged time series of ozone profiles is from late 2001 until the end of 2018. The monthly mean gridded ozone profile dataset is provided in the altitude range from 10 to 50 km in bins of 10∘ latitude × 20∘ longitude. The merging is performed using deseasonalized anomalies. The created MEGRIDOP dataset can be used for analyses that probe our understanding of stratospheric chemistry and dynamics. To illustrate some possible applications, we created a climatology of ozone profiles with resolved longitudinal structure. We found zonal asymmetry in the climatological ozone profiles at middle and high latitudes associated with the polar vortex. At northern high latitudes, the amplitude of the seasonal cycle also has a longitudinal dependence. The MEGRIDOP dataset has also been used to evaluate regional vertically resolved ozone trends in the stratosphere, including the polar regions. It is found that stratospheric ozone trends exhibit longitudinal structures at Northern Hemisphere middle and high latitudes, with enhanced trends over Scandinavia and the Atlantic region. This agrees well with previous analyses and might be due to changes in dynamical processes related to the Brewer–Dobson circulation.

In this paper, we discuss the method for validation of random uncertainties in the remote sensing measurements based on evaluation of the structure function, i.e., root-mean-square differences as a function of increasing spatiotemporal separation of the measurements. The limit at the zero mismatch provides the experimental estimate of random noise in the data. At the same time, this method allows probing of the natural variability of the measured parameter. As an illustration, we applied this method to the clear-sky total ozone measurements by the TROPOspheric Monitoring Instrument (TROPOMI) on board the Sentinel-5P satellite. We found that the random uncertainties reported by the TROPOMI inversion algorithm, which are in the range 1–2 DU, agree well with the experimental uncertainty estimates by the structure function. Our analysis of the structure function has shown the expected results on total ozone variability: it is significantly smaller in the tropics compared to mid-latitudes. At mid-latitudes, ozone variability is much larger in winter than in summer. The ozone structure function is anisotropic (being larger in the latitudinal direction) at horizontal scales larger than 10–20 km. The structure function rapidly grows with the separation distance. At mid-latitudes in winter, the ozone values can differ by 5 % at separations 300–500 km. The method discussed is a powerful tool in experimental estimates of the random noise in data and studies of natural variability, and it can be used in various applications.