the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
![](https://acp.copernicus.org/articles/24/9177/2024/acp-24-9177-2024-avatar-thumb150.png)
Investigation of the impact of satellite vertical sensitivity on long-term retrieved lower-tropospheric ozone trends
Richard J. Pope
Fiona M. O'Connor
Mohit Dalvi
Brian J. Kerridge
Richard Siddans
Barry G. Latter
Brice Barret
Eric Le Flochmoen
Anne Boynard
Martyn P. Chipperfield
Wuhu Feng
Matilda A. Pimlott
Sandip S. Dhomse
Christian Retscher
Catherine Wespes
Richard Rigby
Ozone is a potent air pollutant in the lower troposphere and an important short-lived climate forcer (SLCF) in the upper troposphere. Studies investigating long-term trends in the tropospheric column ozone (TCO3) have shown large-scale spatio-temporal inconsistencies. Here, we investigate the long-term trends in lower-tropospheric column ozone (LTCO3, surface–450 hPa sub-column) by exploiting a synergy of satellite and ozonesonde data sets and an Earth system model (UK's Earth System Model, UKESM) over North America, Europe, and East Asia for the decade 2008–2017. Overall, we typically find small LTCO3 linear trends with large uncertainty ranges using the Ozone Monitoring Instrument (OMI) and the Infrared Atmospheric Sounding Interferometer (IASI), while model simulations indicate a stable LTCO3 tendency. The satellite a priori data sets show negligible trends, indicating that any year-to-year changes in the spatio-temporal sampling of these satellite data sets over the period concerned have not artificially influenced their LTCO3 temporal evolution. The application of the satellite averaging kernels (AKs) to the UKESM simulated ozone profiles, accounting for the satellite vertical sensitivity and allowing for like-for-like comparisons, has a limited impact on the modelled LTCO3 tendency in most cases. While, in relative terms, this is more substantial (e.g. on the order of 100 %), the absolute magnitudes of the model trends show negligible change. However, as the model has a near-zero tendency, artificial trends were imposed on the model time series (i.e. LTCO3 values rearranged from smallest to largest) to test the influence of the AKs, but simulated LTCO3 trends remained small. Therefore, the LTCO3 tendencies between 2008 and 2017 in northern-hemispheric regions are likely to be small, with large uncertainties, and it is difficult to detect any small underlying linear trends due to interannual variability or other factors which require further investigation (e.g. the radiative transfer scheme (RTS) used and/or the inputs (e.g. meteorological fields) used in the RTS).
- Article
(3644 KB) - Full-text XML
-
Supplement
(3140 KB) - BibTeX
- EndNote
-
Satellite lower-tropospheric column ozone (LTCO3) records in the Northern Hemisphere show small trends with large uncertainty ranges between 2008 and 2017.
-
Modelled LTCO3 over that period is temporally stable, and application of the satellite averaging kernels (AKs), accounting for the satellite vertical sensitivity, to the model yields little impact on the simulated trends.
Tropospheric ozone (TO3) is a short-lived climate forcer (SLCF) and an important greenhouse gas (GHG; Myhre et al., 2013; Forster et al., 2021). TO3 is also a hazardous air pollutant, with adverse impacts on human health (Doherty et al., 2017; WHO, 2022) and on agricultural and natural vegetation (Sitch et al., 2007; Hollaway et al., 2012). Since the pre-industrial (PI) period, anthropogenic activities have increased the atmospheric loading of ozone (O3) precursor gases, most notably methane (CH4) and nitrogen oxides (NOx), resulting in an increase in TO3 of 25 %–50 % since 1900 (Gauss et al., 2006; Lamarque et al., 2010; Young et al., 2013). The PI to present-day (PD) radiative forcing (RF) from TO3 is estimated by the Intergovernmental Panel on Climate Change (IPCC) to be 0.47 Wm−2 (Forster et al., 2021), with an uncertainty range of 0.24–0.70 Wm−2.
During the satellite era (i.e. since the mid-1990s), extensive records of TO3 have been produced, e.g. by the European Space Agency Climate Change Initiative (ESA-CCI; ESA, 2019). However, the large presence of stratospheric O3, coupled with the different vertical sensitivities and sources of error associated with observations in different wavelength regions (e.g. Eskes and Boersma 2003; Ziemke et al., 2011; Miles et al., 2015), means large-scale inconsistencies in time and space exist between the records of satellite tropospheric column ozone (TCO3) (as shown by Gaudel et al., 2018).
The work by Gaudel et al. (2018) was part of the Tropospheric Ozone Assessment Report (TOAR), which represented a large global effort to understand spatio-temporal patterns and variability in TO3. Their investigation of ozonesondes (2003–2012) and products from nadir-viewing satellites in polar orbits (three from the Ozone Monitoring Instrument (OMI) (2005–2015 and 2016) and two from the Infrared Atmospheric Sounding Interferometer (IASI) (2008–2016)) displayed discrepancies in the spatial distribution, magnitude, direction, and significance of the TCO3 trends. They noted that the records cover slightly different time periods but were unable to provide any definitive reasons for these discrepancies beyond briefly suggesting that differences in measurement techniques and retrieval methods were likely to be causing the observed spatial inconsistencies. The range of potential definitions of the tropopause height used to derive TCO3 from these nadir-viewing profile products could also lead to differences between the satellite product absolute values and their temporal evolutions. While the five products discussed above use the same definition (i.e. World Meteorological Organization (WMO) 2 K km−1 lapse rate; WMO, 1957), several of the other products analysed by Gaudel et al. (2018) did use other definitions.
The vertical sensitivity of each retrieved product (function of measurement technique and retrieval methodology) used by Gaudel et al. (2018) will have had an impact on which part of the troposphere the O3 signal is weighted towards. This is potentially one of the drivers behind the different OMI and IASI TCO3 trends, where OMI showed predominantly positive trends between 60° S and 60° N, while the opposite was the case for IASI. The vertical sensitivity is represented by the averaging kernel (AK), which provides the relationship between perturbations at different levels in the retrieved and true profiles (Eskes and Boersma, 2003). Typically, for the products used by Gaudel et al. (2018), the peak AK sensitivities for TO3 are in the 0–6 km range for OMI (Miles et al., 2015) and around 11–12 km for IASI (Keim et al., 2009), while there is a secondary peak at approximately 5 km (Boynard et al., 2009). In the case of the Rutherford Appleton Laboratory (RAL) Space OMI data, used in Gaudel et al. (2018), TCO3 values were derived from retrieved surface–450 hPa layer average mixing ratios; this was also applied to the overlying 450 hPa–tropopause layer using ERA-Interim profiles. As the TO3 values were derived from different (UV and IR) sensors and methodologies whose vertical sensitivities differ, they were likely to represent O3 controlled by different contributions of atmospheric processes (e.g. precursor emissions from the surface and stratosphere–troposphere exchanges). Therefore, TCO3 trends from the different satellite products are not necessarily expected to be similar. The determination of the linear trend in a satellite TCO3 record(s) can also be difficult as many factors (e.g. chemistry, emissions, deposition, and transport) which control ozone interannual variability, especially for time periods of a decade or less (Barnes et al., 2016; Chang et al., 2020; Fiore et al., 2022).
In this study, we undertake the first assessment of spatio-temporal variability in satellite-derived lower-tropospheric column ozone (LTCO3, surface–450 hPa) from three satellite products over a consistent decade (2008–2017). In combination with an Earth system model (ESM), we aim to quantify the impact of year-to-year spatio-temporal sampling, the satellite instrument uncertainties, and the instrument vertical sensitivity on long-term LTCO3 trends. We focus our analysis on North America, Europe, and East Asia given their large emissions of ozone precursor gases and their temporal variability. In our article, Sect. 2 discusses the satellite and ozonesonde data sets and the model used, Sect. 3 presents our results, and our discussion and conclusions are summarized in Sects. 4 and 5.
2.1 Satellite data sets
The satellite products (see Table 1) used here are from nadir-viewing polar-orbiting platforms providing ozone sub-column profiles. This includes ozone profile data from the OMI product developed by the RAL Space and the IASI products from the Laboratoire d'Aérologie (IASI-SOFRID) and the Université Libre de Bruxelles, in collaboration with the Laboratoire Atmosphères, Observations Spatiales (ULB-LATMOS) (IASI-FORLI). OMI and IASI are on NASA's Aura and EUMETSAT's Metop-A satellites in sun-synchronous low Earth orbits at 13:30 and 09:30 in local solar overpass time, respectively. OMI and IASI are ultraviolet–visible (UV–Vis) and infrared (IR) sounders with spectral ranges of 270–500 nm (Boersma et al., 2008, 2011) and 645–2760 cm−1 (Illingworth et al., 2011), respectively. OMI has a spatial footprint at nadir of 24 km × 13 km, while IASI measures simultaneously in four fields of view (FOVs, each circular at nadir with a diameter of 12 km) in a 50 km × 50 km square; these are scanned across-track to sample a 2200 km wide swath (Clerbaux et al., 2009).
The OMI retrieval scheme is based on an optimal-estimation (OE) approach, produced by RAL Space, which is described in detail by Miles et al. (2015). The retrieval schemes for IASI-FORLI and IASI-SOFRID O3 are discussed in detail by Boynard et al. (2018) and Barret et al. (2020). The lowest sub-column in the OMI sub-column profile represents the surface–450 hPa layer (i.e. LTCO3). For the IASI products, there were several sub-columns spanning the surface–450 hPa range. Therefore, the IASI sub-columns were totalled up between the surface and the layer beneath or equal to the 450 hPa level. Where the 450 hPa level was located within a sub-column (i.e. was located between its bounding upper and lower pressure levels), the sub-column proportion between the lower pressure barrier and the 450 hPa level was determined and added to the sub-columns below (i.e. towards the surface). For the ozone a priori profile, the RAL Space and FORLI schemes use the ozone latitude vs. month of the year climatology of McPeters et al. (2007), while IASI-SOFRID uses the dynamical ozone climatology described in Sofieva et al. (2014). However, the FORLI scheme uses a single ozone profile (Boynard et al., 2018) derived from the McPeters et al. (2007) data set; thus, it has no seasonality or latitude dependence, unlike the other retrieval schemes.
In this work, the OMI data were filtered for good-quality retrievals where the geometric cloud fraction was < 0.2, the sub-column O3 values were > 0.0, the solar zenith angle was < 80.0°, the retrieval convergence flag was 1.0, and the normalized cost function was < 2.0. The IASI-FORLI data were filtered for a geometric cloud fraction < 0.13 (pre-filtered), degrees of freedom > 2.0, O3 values > 0.0, a solar zenith angle < 80.0°, and a ratio of LTCO3 to total column O3 < 0.085. The IASI-SOFRID data were provided on a 1.0° × 1.0° horizontal grid (i.e. level-3 product but at a daily temporal resolution – we use the daytime data in this study) with filtering already applied, as in Barret et al. (2020). Here, only O3 values > 0.0 were used. To remove systematic biases between the satellite records while maintaining the long-term interannual variability of each record, ozonesondes were used to generate bias correction offsets (BCOs) (2008–2017) to help harmonize the data sets (i.e. subtraction term in units of Dobson units, DU – as done in Russo et al., 2023; Pope et al., 2024); this is discussed in the Supplement (i.e. Sect. S1 in the Supplement). By applying the BCOs, this improves the robustness of the satellite data sets (in absolute terms). This is important when intercomparing the products but also when using them to evaluate the UK's Earth System Model (UKESM) and when determining the model's skill in simulating LTCO3 as used in this study (see Sect. S4 in the Supplement).
Here, each ozonesonde profile was co-located with the nearest satellite retrieval within 500 km and 6 h to reduce spatio-temporal sampling biases (e.g. Keppens et al., 2019). The ozonesonde profile was then interpolated in the vertical onto the satellite pressure grid where the sub-columns between pressure levels were determined. The ozonesonde sub-column profiles were then convolved by the satellite averaging kernels (AKs), which represent the sensitivity of the satellite's ability to retrieve ozone as a function of altitude, thus allowing for a robust like-for-like comparison between the ozonesondes and the retrieved LTCO3. The application of AKs to ozonesonde profiles to evaluate satellite ozone products is discussed in detail by Pope et al. (2023). The application of the AKs to the ozonesondes (and the model) is outlined in Eq. (1):
where sondeAK is the modified ozonesonde sub-column profile (Dobson units, DU), AK is the averaging-kernel matrix, sondeint is the sonde sub-column profile (DU) on the satellite pressure grid, and apr is the a priori (DU). The application of the AKs to the ozonesondes is discussed in more detail in Sect. S1 in the Supplement.
To investigate long-term trends over North America, Europe, and East Asia, the Hemispheric Transport of Air Pollution (HTAP) regional sea–land mask (European Commission, 2016; see Sect. S2 in the Supplement, Fig. S5) is used to sub-sample the gridded satellite data for the respective regions and then generate average monthly time series for each product over each region of interest. For the ozonesonde time series for each HTAP region investigated, only ozonesonde sites which are located within each HTAP region are selected. This results in 15, 13, and 6 ozonesonde sites for North America, Europe, and East Asia, respectively. As ozonesonde data for East Asia are all from Japan, Taiwan, and Hong Kong, trends in ozone LTCO3 will likely be different to satellite and/or model trends over all of East Asia.
In Sect. 3.2, where we discuss the impact of satellite retrieval errors on derived LTCO3 linear trends, the OMI and IASI-FORLI retrieval errors are provided in their product files, but these are not available for IASI-SOFRID. Therefore, while it is not a perfect metric to represent the error in the IASI-SOFRID data, we use the standard deviation in the monthly spatial average of the regional time series.
2.2 United Kingdom Earth System Model (UKESM)
The UK's Earth System Model, UKESM1.0, is a state-of-the-art ESM with fully interactive coupled-component models (e.g. atmosphere, ocean, land surface, atmospheric chemistry) which has been developed by the UK Met Office and the Natural Environment Research Council (NERC). The detailed coupling of all the Earth system components is described by Sellar et al. (2019). However, in this study, we run UKESM1.0 in an atmosphere-only configuration (e.g. similarly to Archibald et al., 2020). The aim is to use UKESM1.0 to investigate long-term trends in TO3 and to help explore inconsistencies between satellite records; thus, it is computationally more time efficient as only the atmospheric dynamics and chemistry components are simulated. Over the 2008–2017 time period (with a 1-year spin-up), the UKESM1.0 model tracers and diagnostics (e.g. ozone, pressure) are output as 3D fields at sub-daily (6-hourly) time steps to allow robust comparisons between the model and satellite data sets (i.e. model–satellite spatio-temporal co-location to reduce representation biases and application of the satellite AKs to map the instrument vertical sensitivity onto the model, yielding like-for-like comparisons). The satellite AKs from OMI and IASI-FORLI are provided in the level-2 files (i.e. an AK matrix per retrieval). However, the IASI-SOFRID AKs are provided from the gridded level-3 data product (i.e. an AK matrix for each 1° × 1° grid box).
Here, the UKESM1.0 land and atmosphere share a regular latitude–longitude grid with a resolution of 1.25° × 1.875°, with 85 vertical levels on a terrain-following hybrid height coordinate and with a model lid at 85 km above sea level (50 levels are below 18 km). All the key inputs to the model from other Earth system components (e.g. sea surface temperature (SST) and land surface vegetation) were prescribed from ancillary files. The ocean and ice forcings are represented by the monthly Reynolds sea ice and SST data from the National Oceanic and Atmospheric Administration (NOAA, https://climatedataguide.ucar.edu/climate-data/, last access: 1 July 2022). Solar forcings are provided by Phase 6 of the Coupled Model Intercomparison Project (CMIP6; Matthes et al., 2017; Eyring et al., 2016), as is the stratospheric aerosol climatology in order to represent contributions from volcanic eruptions (Sellar et al., 2019). The land cover is provided by the output from the land surface component of the ESM (JULES; Wiltshire et al., 2021) from a fully coupled historical simulation. Anthropogenic and biomass burning emissions from Hoesly et al. (2018) and van Marle et al. (2017) are prescribed for the period 2008–2014. After 2014, anthropogenic and biomass burning emissions are from the Shared Socioeconomic Pathway (SSP, Rao et al., 2017) 2–4.5 (i.e. a middle-of-the-road climate and emission scenario).
Biological emissions are a climatology between 2001 and 2010 from the MEGAN-MACC database (Sindelarova et al., 2014), while natural emissions are from the Precursors of Ozone and their Effects in the Troposphere (POET, http://accent.aero.jussieu.fr/database_table_inventories.php, last access: 1 July 2022) based on the year 1990. Dry deposition of O3 to the land surface is represented by the Wesley scheme, which is applied as in O'Connor et al. (2014). The model is also in a nudged or “specified-dynamics” configuration (i.e. meteorological analyses are used to “nudge” the model's meteorological variables, i.e. u- and v-wind components and potential temperature, towards reality; Telford et al., 2008) using 6-hourly reanalysis data from the European Centre for Medium-Range Weather Forecasts (ECMWF) ERA-Interim product. A similar configuration of UKESM1.0 was used by Archibald et al. (2020), in which a thorough evaluation against multiple observations (e.g. surface, aircraft, and satellite) was carried out.
2.3 Trend approach
LTCO3 trends are calculated using the linear least-squares fit approach formulated by van der A et al. (2006, 2008) and utilized by Pope et al. (2018), who investigated LTCO3 trends. Here, the monthly LTCO3 time series are represented by the following function:
where Yt is the observed monthly LTCO3 for month t, Xt is the number of months since the start of the record, C is the first monthly mean LTCO3 value of the record, B is the monthly linear trend, and Asin(ωXt+ϕ) is the seasonal model component (Weatherhead et al., 1998). A is the amplitude, ω is the frequency (set to 1 year; ), and ϕ is the phase shift. C, B, A, and ϕ are the fit parameters from the linear least-squares fit. Nt represents the model errors and/or residuals. The linear trend uncertainty, σB, represents the trend precision and is calculated as follows:
where n is the number of years, α is the autocorrelation in the residuals (Nt), and σN is the standard deviation in the residuals. As in van der A et al. (2006) and Pope et al. (2018), we calculate the autocorrelation for each time series using a lag of one time step (i.e. 1 month). The autocorrelation in Eq. (2) is not accounted for directly; thus, it is factored into the trend uncertainty (Eq. 3), as used and discussed by van der A et al. (2006) and Weatherhead et al. (1998), respectively.
A detailed evaluation of UKESM1.0 LTCO3 through comparisons with the three satellite products and ozonesondes is presented in Sect. S4 in the Supplement. Overall, UKESM1.0 robustly simulates LTCO3 spatially and seasonally in comparison to the ozonesondes and satellite instruments (i.e. typically within the ozonesonde variability and satellite uncertainty range).
3.1 UKESM1.0 and satellite LTCO3 trends
3.1.1 North America
LTCO3 trends from OMI, IASI-FORLI, IASI-SOFRID, and ozonesondes are derived between 2008 and 2017 (i.e. consistent time record for all instruments) using the linear–seasonal trend model (Eq. 2). For each satellite product, the corresponding UKESM1.0 time series (with and without AKs) are analysed, along with the satellite a priori. For the North America OMI metrics (Fig. 1a, Table 2), there is clear seasonality in the a priori, ranging between approximately 17.0 and 22.0 DU. As this is based on the climatology of McPeters et al. (2007), there is no trend, and there is a very good model fit (i.e. R2=1.0). The key point is that, as a climatology, the a priori will have no trend, but if there are substantial temporal sampling differences between years then an artificial trend could be introduced. OMI LTCO3 ranges between 20.0 and 27.0 DU, with substantial variability. There is a drop in LTCO3 to 19.0 DU in 2009 before peaking at 25.0–27.0 DU between 2010 and 2015. Peak LTCO3 then drops to 22.0–24.0 DU in 2016 and 2017. As a result, the linear–seasonal trend model, which does not account for interannual variations such as this, only has a fit skill of R2=0.59. The corresponding OMI LTCO3 trend is −0.79[−7.07, 5.48] DU per decade, showing a negligible trend with a large uncertainty range. Here, −0.79 DU per decade is the trend, while [−7.07, 5.48] DU per decade constitute the 95 % confidence interval. The UKESM1.0 LTCO3 time series ranges between 17.0 and 22.0 DU, with clear seasonality though somewhat less interannual variation than OMI, and the linear–seasonal trend model therefore has a considerably better fit with R2=0.95. The model trend has the opposite sign at 0.21[−0.37, 0.78] DU per decade. Here, the model trend is near-zero, with a relatively large uncertainty range (though not as sizable as OMI). When the AKs are applied to the model, the trend switches sign to −0.57[−1.58, 0.45] DU per decade, and the linear–seasonal trend model fit decreases in skill to R2=0.90. The trend switch in sign, though small, is potentially linked to the application of the AKs, which also increases LTCO3 by 2.0–3.0 DU in general.
Table 2LTCO3 trends (DU per decade) for the satellite trend (trend), the satellite − uncertainty trend (trend error 1), the satellite + uncertainty trend (trend error 2), the satellite a priori trend (a priori trend), the UKESM trend (UKESM trend), the trend of UKESM with AKs applied (UKESM + AKs trend), the forced UKESM trend (forced UKESM trend) and the trend of the forced UKESM with AKs applied (forced UKESM + AKs trend). The “trend lower” and “trend upper” represent the trend 95 % confidence interval based on the trend precision calculated from Eq. (3). R2 is the trend fit skill (i.e. correlation squared), and the p value is also shown.
![](https://acp.copernicus.org/articles/24/9177/2024/acp-24-9177-2024-t02-thumb.png)
We also investigated the satellite degrees of freedom of signal (DOFS) over the lower troposphere (i.e. surface to 450 hPa), which provides an estimate of the number of independent pieces of information in the LTCO3. The DOFS are calculated by taking the trace of the AK matrix over the lower-tropospheric levels in the satellite vertical grid. Overall, we found that the products for the three regions had negligible trends in their time series (i.e. within ± 1.0 % yr−1), meaning that the information content of satellite LTCO3 remained stable with time (see Sect. S3 in the Supplement).
![https://acp.copernicus.org/articles/24/9177/2024/acp-24-9177-2024-f01](https://acp.copernicus.org/articles/24/9177/2024/acp-24-9177-2024-f01-thumb.png)
Figure 1Lower-tropospheric column ozone (LTCO3, surface to 450 hPa, DU) regional time series for North America, based on the HTAP land mask, from OMI (a), IASI-FORLI (b), IASI-SOFRID (c), and ozonesondes (d), are shown by the black lines in the respective panels. UKESM1.0 simulations without and with satellite averaging kernels (AKs) applied are shown in red and blue lines. Green lines show the satellite a priori. Dashed lines show the LTCO3 linear trends, which are labelled at the top of each panel. The R2-squared values show the linear–seasonal trend model fit to the corresponding LTCO3 time series (i.e. correlation squared).
The IASI-FORLI LTCO3 time series (Fig. 1b) tends to be lower than that of OMI and ranges between 17.0 and 22.0 DU. There is a substantial negative IASI-FORLI trend (−1.42[−2.35, −0.50] DU per decade; Table 2), although, as stated by Boynard et al. (2018) and Wespes et al. (2018), the IASI level-1 data sets input into the FORLI retrieval are not consistent with time; they suffer from a specific discontinuity in September 2010, which degrades the robustness of this trend. While we are aware of the artificial trend in the IASI-FORLI data set, it is still a valuable long-term product, allowing us to quantify multiple factors (e.g. impact of AKs on model tendencies and absolute values and year-to-year spatio-temporal sampling stability – i.e. near-zero trend in the a priori). The a priori has a negligible trend, but there is no clear seasonality in the a priori time series. As a result, the linear–seasonal trend model has a more limited fit skill (i.e. R2=0.67). The impact of the satellite AKs appears to be less for IASI-FORLI as both UKESM1.0 and UKESM1.0 + AKs have time series ranging between approximately 17.0 and 21.0 (though the UKESM1.0 + AKs range is slightly smaller) and linear–seasonal trend model fits of R2=0.93 and R2=0.92, respectively. The corresponding trends are small at −0.13[−0.75, 0.49] DU per decade and −0.32 [−0.82, 0.20] DU per decade, but the introduction of the AKs does move the UKESM1.0 trend slightly towards that of the satellite. Interestingly, while the application of the IASI-FORLI AKs to UKESM marginally pushes the convolved model trend in LTCO3 towards that of the satellite (which has a substantial negative trend), the IASI-FORLI DOFS have small positive trends (0.37 % yr−1–0.57 % yr−1 – see Sect. S3 in the Supplement). Therefore, there is a minor-scale, yet contrasting, discrepancy in how the vertical sensitivity influences the long-term LTCO3 trends.
For IASI-SOFRID (Fig. 1c), there is little difference between any of the time series as they all range between 16.0 and 21.0 DU, with corresponding linear–seasonal trend model fits of R2=0.94 to 0.98 and negligible trends. The IASI-SOFRID and a priori trends are 0.12[−0.59, 0.82] DU per decade and 0.11[−0.17, 0.39] DU per decade (Table 2), respectively, with the model showing near-zero trends in both cases. Given the close agreement between the satellite and a priori time series and fit metrics, it is suggestive that IASI-SOFRID TO3 is more closely confined to the a priori profile than are the other products.
The ozonesondes show a substantial trend of −1.15[−2.0, −0.10] DU per decade, while the model trend sampled as the sondes is −0.16[−1.67, 1.35] DU per decade. The co-located model and ozonesonde linear–seasonal trend model fits are R2=0.62 and 0.64, respectively. The noise and lack of seasonality in the ozonesonde time series are slightly unexpected given the reasonable density of stations across North America, though the spatial coverage and temporal sampling are much less than for the satellite products.
3.1.2 Europe
In Europe, the OMI LTCO3 values are larger than in North America, ranging between 19.0 and 30.0 DU (Fig. 2a). The same interannual variability exists, peaking between 2010 and 2015, with the minimum in 2009. Hence, the linear–seasonal trend model, which does not represent interannual variation, has moderate skill and R2=0.72. The corresponding trend is −0.80[−7.29, 5.69] DU per decade and so has a similar direction and magnitude to that for North America, though this is not substantial. The a priori ranges between 17.0 and 22.5 DU, with a trend of −0.12[−0.26, 0.03] DU per decade (Table 2). Given the relatively small trend and uncertainty range, unlike for the OMI equivalent, this suggests that there is unlikely to be an artificial trend arising through year-to-year spatio-temporal sampling changes in geographical sampling across the European region. UKESM1.0 LTCO3 ranges between approximately 19.0 and 22.0 DU, with a good linear–seasonal trend model fit of R2=0.99 and a trend of −0.11[−0.50, 0.29] DU per decade. As for North America, when the OMI AKs are applied, the UKESM LTCO3 values systematically increase by 2.0–3.0 DU, move further away from the satellite a priori, and more closely follow the variability of OMI (R2 decreases slightly to 0.95). The trend tends towards that of OMI at −0.72[−1.77, 0.32] DU per decade.
![https://acp.copernicus.org/articles/24/9177/2024/acp-24-9177-2024-f02](https://acp.copernicus.org/articles/24/9177/2024/acp-24-9177-2024-f02-thumb.png)
Figure 2LTCO3 (DU) regional time series for Europe, based on the HTAP land mask, from OMI (a), IASI-FORLI (b), IASI-SOFRID (c), and ozonesondes (d), are shown by the black lines in the respective panels. UKESM1.0 simulations without and with satellite AKs applied are shown in red and blue lines. Green lines show the satellite a priori. Dashed lines show the LTCO3 linear trends, which are labelled at the top of each. The R2-squared values show the linear–seasonal trend model fit to the corresponding LTCO3 time series (i.e. correlation squared).
As in the case of North America, the European IASI-FORLI a priori has no seasonal cycle (and moderate R2 of 0.48 in the linear–seasonal trend model fit) with a near-zero trend (0.09[−0.09, 0.27] DU per decade) (Fig. 2b, Table 2). The IASI-FORLI data exhibit a substantial negative trend of −1.83[−2.78, −0.89] DU per decade, again due to step changes in the IASI level-1 processor, with a good linear–seasonal trend model fit of R2=0.92. UKESM1.0 LTCO3 trends, without and with AKs applied, are −0.28[−0.77, 0.20] DU per decade and −0.43[−1.21, 0.35] DU per decade. Again, though a small change, the application of the AKs introduces a slight perturbation of the model trend compared to IASI-FORLI.
The IASI-SOFRID a priori, ranging between 17.0 and 21.0 DU, has a trend of 0.17[−0.12, 0.45] DU per decade, with a good fit skill of R2=0.98 (Fig. 2c). The IASI-SOFRID and UKESM1.0 metrics, with and without averaging kernels applied, are similar, with LTCO3 trends of 0.05[−0.91, 1.01] DU per decade, −0.27 [−0.72, 0.19] DU per decade, and 0.08[−0.33, 0.49] DU per decade, respectively, and with R2 values between 0.93 and 0.98.
The ozonesonde monthly regional means (Fig. 2d) have a more pronounced time series than for North America, yielding a less noisy time series of LTCO3. Here, there is clear seasonality ranging between 17.0 and 24.0 DU, with a large R2 value of 0.95. The ozonesonde trend is relatively small at −0.6[−1.39, 0.17] DU per decade, while the UKESM1.0 equivalent is more substantial at −0.96[−1.56, 0.35] DU per decade.
3.1.3 East Asia
For East Asia, OMI LTCO3 again has both a pronounced seasonal cycle and interannual variability (19.0–27.0 DU), which is consistent with the other two regions discussed above (Fig. 3a, Table 2). This yields a moderate skill fit to the linear–seasonal trend model of R2=0.52 and a near-zero trend (−0.09[−7.88, 7.70] DU per decade). The a priori has a trend of −0.25[−0.71, 0.22] DU per decade; thus, year-to-year spatio-temporal sampling changes could be influencing the robustness of OMI retrieved time series in this region. However, both the instrument and a priori trend uncertainties intersect with 0.0. UKESM1.0 LTCO3 ranges between approximately 16.0 and 22.0 DU, with a good fit of R2=0.98. Like the other regions, the application of the OMI AKs increases the model values systematically by several Dobson units. The UKESM1.0 LTCO3 trend is −0.16[−0.94, 0.62] DU per decade, which is small, but the AKs increase the trend magnitude to −0.62[−2.24, 1.00] DU per decade, which moves it away from the OMI trend.
![https://acp.copernicus.org/articles/24/9177/2024/acp-24-9177-2024-f03](https://acp.copernicus.org/articles/24/9177/2024/acp-24-9177-2024-f03-thumb.png)
Figure 3LTCO3 (DU) regional time series for East Asia, based on the HTAP land mask, from OMI (a), IASI-FORLI (b), IASI-SOFRID (c), and ozonesondes (d), are shown by the black lines in the respective panels. UKESM1.0 simulations without and with satellite AKs applied are shown in red and blue lines. Green lines show the satellite a priori. Dashed lines show the LTCO3 linear trends, which are labelled at the top of each panel. The R2-squared values show the linear–seasonal trend model fit to the corresponding LTCO3 time series (i.e. correlation squared).
IASI-FORLI (Fig. 3b, Table 2), like the other two regions, has a substantial negative trend of −1.52[−2.16, −0.88] DU per decade. The a priori again exhibits virtually no seasonal cycle (low fit skill of R2=0.21) and negligible year-to-year spatio-temporal sampling differences yielding a near-zero trend of −0.03[−0.22, 0.16] DU per decade. For UKESM1.0, the East Asian seasonal range is much larger than that of other regions, ranging between 17.0 and 27.0 DU (i.e. seasonal amplitude of approximately ± 5.0 DU). When the AKs are applied, this range shrinks to approximately 19.0 to 23.0 DU, which is more in line with the IASI-FORLI LTCO3 values. The corresponding model trends are −0.03[−0.62, 0.56] DU per decade and −0.29[−0.80, 0.22] DU per decade; thus, the AKs push the model tendency towards that of the instrument, though the impact is small in absolute terms (large in relative terms).
IASI-SOFRID and its a priori LTCO3 seasonality are, again, very similar, ranging between 16.0 and 21.0 DU, with very little interannual variability and with linear–seasonal trend model fit skills of R2=0.96 and 0.98 (Fig. 3c, Table 2). The IASI-SOFRID and a priori linear trends are therefore also consistent at −0.19[−1.01, 0.63] DU per decade and −0.15[−0.73, 0.58] DU per decade. The UKESM1.0 seasonal variability is, again, large, between 17.0 and 26.0 DU, and, as in the case of IASI-FORLI, when the instrument AKs are applied to the model, the seasonal range shrinks (i.e. 16.0–22.0 DU) to be much closer to the ranges of the retrieval and its prior. The model trends are −0.42[−0.97, 0.13] DU per decade and −0.24[−0.67, 0.20] DU per decade (with AKs), where there is a minor shift in the model tendency towards that of IASI-SOFRID and its prior.
For the ozonesondes (Fig. 3d), there is a substantial LTCO3 trend of 3.17[0.16, 6.17] DU per decade, with a fit skill of R2=0.79, which is larger than those for North America and Europe. LTCO3 increases from 18.0–25.0 in 2008 to 21.0–28.0 in 2011. This remains similar in 2012 and 2013 before dropping by several Dobson units between 2014 and 2017. The UKESM1.0 sampled as the ozonesondes has considerably less interannual variability, with a smaller trend of 0.37[−0.90, 1.64] DU per decade. Therefore, UKESM1.0 and the satellite product trends are generally smaller (in magnitude) than the ozonesonde tendencies. However, it is worth considering that there are only a few sites (e.g. Hong Kong and Taiwan) where ozonesonde data are available in East Asia.
3.2 Influence of satellite averaging kernels on UKESM1.0 LTCO3
To investigate the impact of applying the satellite averaging kernels to UKESM1.0 – and, thus, to learn something about the influence of vertical sensitivity on retrieved LTCO3 – three different metrics are considered for the 2008–2017 time-period. These are the absolute LTCO3 value, the amplitude of the LTCO3 seasonal cycle, and the linear trend. These metrics are compared for the satellite, the satellite ± the error term, the a priori, UKESM1.0, and UKESM1.0 + AKs for the three regions discussed above.
From Fig. 4, the average OMI LTCO3 is approximately 22.0, 24.0, and 23.0 DU for North America, Europe, and East Asia, respectively. This represents a substantial deviation from the a priori values of 17.5, 20.0, and 16.0 DU, respectively. However, the average error term for OMI LTCO3 is sizeable at approximately ± 8.0 to ± 9.0 DU for all regions. The average UKESM1.0 values for each region are approximately 19.5, 21.5, and 19.0 DU, but the application of the AKs increases this by several Dobson units to 22.0, 24.0, and 21.0 DU. In comparison, mean values for both IASI products vary less between the three geographical areas: IASI-FORI (IASI-SOFRID) LTCO3 values are 20.0 (18.5), 19.0 (18.5), and 22.0 (18.0) DU. The corresponding error ranges, in comparison with OMI, are smaller between 17.0 and 23.0 (16.0 and 21.5), 16.0 and 21.5 (16.0 and 21.0), and 18.0 and 23.5 (14.5 and 21.5) DU for North America, Europe, and East Asia, respectively. With the IASI-FORLI AKs applied to UKESM1.0, LTCO3 decreases from 19.5–19.25, 21.25–19.5, and 22.75–21.25 DU for the three regions. For IASI-SOFRID, there is a decrease from 21.0–19.5 DU in Europe and a decrease from 22.0–19.5 DU in East Asia, while no change occurs in North America. Overall, OMI has the largest error range, and the application of the AKs to UKESM1.0 systematically increases the model LTCO3 time series by several Dobson units. The opposite occurs for the IASI products where there is a smaller decrease to UKESM1.0 LTCO3 of 1.0–2.0 DU. The error ranges are also smaller than that of OMI.
![https://acp.copernicus.org/articles/24/9177/2024/acp-24-9177-2024-f04](https://acp.copernicus.org/articles/24/9177/2024/acp-24-9177-2024-f04-thumb.png)
Figure 4Average LTCO3 (DU) values across the 2008–2017 time period for the satellite (black), lower satellite (cyan), upper satellite (grey), a priori (blue), UKESM1.0 (red), and UKESM1.0 + AKs (orange). The lower-satellite and upper-satellite values are the average of the satellite ± its error term time series (note: these values do not always fit in the y-axis range). O, F, and S represent OMI, IASI-FORLI, and IASI-SOFRID for North America (a), Europe (b), and East Asia (c).
In terms of the LTCO3 seasonal amplitude (Fig. 5), OMI (including the error terms) is approximately 2.6 DU (for all), 3.3–3.8, and 2.3–2.6 DU for North America, Europe, and East Asia. The a priori seasonal amplitude ranges from 2.7 to 2.9 DU across the regions. The IASI-FORLI averages (including the error terms) tend to be lower than OMI but have similar seasonal ranges. North America, Europe, and East Asia have amplitudes of 2.3–2.5, 2.3–2.5, and 1.6–1.8 DU, respectively. It is noteworthy that this seasonal cycle occurs despite the IASI-FORLI prior exhibiting virtually no seasonal cycle at all. IASI-SOFRID has a European range of 2.4–2.6 DU and comparable ranges for North America and East Asia at 1.8–2.5 and 2.3–3.0 DU. Therefore, the seasonal amplitude in IASI-SOFRID is more sensitive to the error metric, but as the “error” term is based on the LTCO3 standard deviation, given the lack of an error term in the product, it is unsurprising that there is more variability in the seasonal amplitude. For the OMI comparisons, the application of the AKs to UKESM1.0 shifts the simulated amplitude slightly upwards from 2.0–2.1, 3.1–3.3, and 4.0–4.4 DU for the respective regions. The IASI-FORLI AK impacts are a decrease from 1.9–1.4, 3.0–2.1, and 4.2–1.9 DU. For IASI-SOFRID, the corresponding impact on UKESM1.0 is 2.2–2.4, 3.3–2.9, and 4.5–3.2 DU. Therefore, the OMI AKs have a minimal impact, increasing the model seasonal amplitude by 0.1–0.3 DU, but the IASI products suppress the simulated amplitude by 1.0–2.0 DU at the most extreme.
![https://acp.copernicus.org/articles/24/9177/2024/acp-24-9177-2024-f05](https://acp.copernicus.org/articles/24/9177/2024/acp-24-9177-2024-f05-thumb.png)
Figure 5Average LTCO3 seasonal-cycle amplitude (DU) values across the 2008–2017 time period for the satellite (black), lower satellite (cyan), upper satellite (grey), a priori (blue), UKESM1.0 (red), and UKESM1.0 + AKs (orange). The lower-satellite and upper-satellite values are the average of the satellite ± its error term time series (note: these values do not always fit in the y-axis range). O, F, and S represent OMI, IASI-FORLI, and IASI-SOFRID for North America (a), Europe (b), and East Asia (c).
The impacts of the satellite LCTO3 error terms on the derived linear trends are shown in Fig. 6. For OMI, the range in trends calculated (i.e. satellite ± error term) is approximately −1.50[−7.04, 4.04] DU per decade to −0.09 [−6.98, 6.81] DU per decade, −1.65[−6.92, 3.62] DU per decade to 0.05[−7.44, 7.53] DU per decade, and −1.05[−6.61, 4.52] DU per decade to 0.87[−8.24, 9.98] DU per decade for North America, Europe, and East Asia, respectively. The IASI-FORLI trends (i.e. satellite ± error term) are substantial, ranging from −1.50[−2.51, −0.50] DU per decade to −1.34[−2.21, −0.47] DU per decade, −1.87[−2.87, −0.87] DU per decade to −1.80[−2.72, −0.88] DU per decade, and −1.62[−2.27, −0.98] DU per decade to −1.42[−2.06, −0.78] DU per decade for the three regions. The corresponding IASI-SOFRID trends were 0.09[−0.48, 0.66] DU per decade to 0.14[−0.59, 0.88] DU per decade, −0.07[−0.91, 0.78] DU per decade to 0.16[−0.74, 1.07] DU per decade, and −0.30[−1.02, 0.42] DU per decade to −0.08[−0.73, 0.58] DU per decade. Therefore, only the IASI-FORLI trends (i.e. satellite ± error term) are substantially different from zero (i.e. p < 0.05). However, that is due in part to discontinuities in the input meteorological data used to generate this version of the product (Boynard et al., 2018).
![https://acp.copernicus.org/articles/24/9177/2024/acp-24-9177-2024-f06](https://acp.copernicus.org/articles/24/9177/2024/acp-24-9177-2024-f06-thumb.png)
Figure 6Average LTCO3 linear trends (DU per decade) across the 2008–2017 time period for the satellite (black), lower satellite (cyan), upper satellite (grey), a priori (blue), UKESM1.0 (red), UKESM1.0 + AKs (orange), forced UKESM1.0 (dark green), and forced UKESM1.0 + AKs (light green). The lower-satellite and upper-satellite values are the average of the satellite ± its error term time series (note: these values do not always fit in the y-axis range). O, F, and S represent OMI, IASI-FORLI, and IASI-SOFRID for North America (a), Europe (b), and East Asia (c). Diamond and circular symbols represent linear trends with p values > 0.05 or p ≤ 0.05, respectively.
The application of the OMI AKs to UKESM1.0 had the largest impacts on the simulated trends, with changes in a negative direction from 0.21[−0.37, 0.78] DU per decade to −0.57[−1.58, 0.45] DU per decade, −0.11[−0.50, 0.29] DU per decade to −0.72[−1.77, 0.32] DU per decade, and −0.16[−0.94, 0.62] DU per decade to −0.62[−2.24, 1.00] DU per decade for the respective regions. IASI-FORLI AKs introduced small decreases from −0.13[−0.75, 0.49] DU per decade to −0.32[−0.82, 0.20] DU per decade, −0.28[−0.77, 0.20] DU per decade to −0.43[−1.21, 0.35] DU per decade, and −0.03[−0.62, 0.56] DU per decade to −0.29[−0.80, 0.22] DU per decade. IASI-SOFRID AKs introduced small increases in the LTCO3 trend from −0.24[−0.85, 0.37] DU per decade to −0.04[−0.53, 0.45] DU per decade, −0.27[−0.72, 0.19] DU per decade to 0.08[−0.33, 0.49] DU per decade, and −0.42[−0.97, 0.13] DU per decade to −0.24[−0.67, 0.20] DU per decade.
As the absolute model trends are small, it is difficult to determine the impact of the AKs on the simulated trends. In relative terms, it can have impacts of several hundred percent, but the model and model + AK trend ranges (95 % confidence interval) always intersect. Therefore, in an attempt to derive more substantial UKESM1.0 LTCO3 trends (without and with AKs applied) in order to assess the maximum impact the AKs can have on UKESM LTCO3 trends, the modelled data were sorted from lowest to highest, and the trend was re-calculated. In North America, this approach forced positive model trends, sub-sampled to OMI, IASI-FORLI, and IASI-SOFRID, of 0.73[0.22, 1.25] DU per decade, 0.64[−3.50, 4.77] DU per decade, and 0.80 (0.41, 1.19) DU per decade. When the AKs were applied, it yielded trends of −0.74[−1.89, 0.40] DU per decade, 0.55[0.08, 1.03] DU per decade, and 0.58[0.24, 0.92] DU per decade. In Europe, this forced positive model trends of 0.62[0.14, 1.10] DU per decade, 0.37[−0.05, 0.79] DU per decade, and 0.46[0.09, 0.84] DU per decade. With the AKs applied, the trends become 0.47[−0.51, 1.44] DU per decade, 0.28[−0.38, 0.94] DU per decade, and 0.10[−0.32, 0.51] DU per decade. Finally, in East Asia, the forced model trends are 0.90[0.34, 1.47] DU per decade, 0.66[0.15, 1.17] DU per decade, and 0.63[0.26, 1.00] DU per decade. The application of the AKs introduced model trends of 1.02[−0.04, 2.09] DU per decade, 0.08[−0.44, 0.61] DU per decade, and 0.20[−0.20, 0.61] DU per decade.
Even with forced trends in the UKESM1.0 regional time series, the trends are relatively small (i.e. typically less than 1.0 DU per decade in magnitude). Therefore, the application of the AKs to the forced UKESM LTCO3 time series still yields small-scale changes in tendencies, and there is an overlap in the two model trend uncertainty ranges (i.e. 95 % confidence levels). However, in relative terms, the trend changes are larger (e.g. > 100 % in multiple cases), and there is often a shift in the modelled LTCO3 trend uncertainty range so that it is either intersecting or no longer intersecting with zero (i.e. a shift in the p-value regime from < 0.05 to > 0.05). Therefore, in modelled and satellite data sets with more substantial trends, the impacts of the AKs, and thus of the satellite vertical sensitivity, on LTCO3 trends would be much greater and could potentially help pinpoint sources of differences between satellite products in terms of their TO3 temporal evolution.
3.3 Diurnal variability in regional LTCO3 and temporal evolution
As TO3 varies diurnally due to meteorological and photochemical processes (e.g. Gaudel et., 2018), the different satellite overpass times (i.e. Aura and Metop-A daytime overpasses are around 13:30 and 09:30 local time, respectively) will likely influence the spatial distributions of TO3 which OMI and IASI will retrieve. In principle, this could therefore explain some differences between the two sensors and their long-term LTCO3 trends. Here, the model is a useful tool to help investigate this, and we used the 6-hourly output to derive the UKESM simulated LTCO3 spatial distributions at the Aura (13:30 LT) and Metop-A (09:30 LT) daytime overpasses. These model fields were then used to calculate regional time series for North America, Europe, and East Asia. For each region and month, between 2008 and 2017, we calculated the regional average absolute difference (i.e. from the selection of model grid cells which fell within the HTAP-2 mask for a specific month) and the standard deviation of the absolute differences between the overpass times. Here, across all months and regions, we found the peak average absolute difference (13:30–09:30 LT) and the standard deviation to be small at 2.03 % and 2.56 %, respectively. For the long-term trends, across all regions and overpass times, all of the UKESM trends were smaller than ± 0.5 DU per decade. Therefore, the model LTCO3 regional trends are negligibly different between overpass times. This might not be surprising given the negligible model trends in the satellite spatio-temporal trend comparisons (see Sect. 3.1), but the actual absolute differences (average and range) in simulated LTCO3 are also small, supporting the argument that, on the regional scale, the daytime diurnal-cycle differences between satellite overpass times have a limited influence on the reported satellite trend discrepancies (e.g. in Gaudel et al., 2018).
This investigation of satellite LTCO3 focussed on 2008–2017, representing a decade of overlap of the OMI and IASI records. The analysis focussed on North America, Europe, and East Asia as these regions are subject to large emissions of and temporal changes in O3 precursor gases. LTCO3 is typically spatially homogeneous, with shallow gradients between background and source-induced O3 concentrations. Secondly, individual retrievals of LTCO3 are often associated with large uncertainties (e.g. random and systematic uncertainties). There are multiple contributory factors concerning both instrumental attributes (notably, spectroradiometric noise and calibration accuracy) and variability in geophysical variables which influence radiative transfer and vertical sensitivity (e.g. stratospheric ozone, cloud and aerosol, water vapour, surface spectral reflectivity and/or emissivity, and pressure and temperature profiles) and which can result in LTCO3 time series with substantial variability and/or noise when derived at a high spatial resolution (e.g. when deriving time series from data gridded at 0.5° or 1.0°). Therefore, we undertake our analysis at the regional (e.g. continental) scale, where more satellite retrievals are included in time series monthly means, yielding a reduction in the random-error component of the sample.
Ideally, this analysis would have utilized several more records (e.g. several UV–Vis and IR products) to quantify long-term trends in LTCO3 and to investigate the potential reasons for any discrepancies, as shown by Gaudel et al. (2018) for TCO3. While RAL Space and other providers have generated UV–Vis profile O3 products for more instruments, e.g. the Global Ozone Monitoring Experiment 1 and 2 (GOME-1 and GOME-2) and the SCanning Imaging Absorption spectroMeter for Atmospheric CHartographY (SCIAMACHY), the GOME-1 and SCIAMACHY records do not overlap with the IASI and OMI records either at all or for a sufficient time period. Secondly, step changes in the GOME-2A level-1 processing scheme used to produce the available LTCO3 level-2 version mean it is not sufficiently homogeneous (see Pope et al., 2023). For the IR instruments, other potential sensors include the Tropospheric Emissions Spectrometer (TES; Richards et al., 2008) and the RAL Space IASI Extended Infrared Microwave Sounding (IMS; Pimlott et al., 2022) scheme applied to IASI. Unfortunately, the TES record only covers 2005–2013, with decreasing spatial coverage with time, and, at the time of this work, the IASI-IMS product had only been processed on a sub-sampled basis of 1 in 10 d.
In this work, we some find discrepancies in the observed long-term tendencies from the utilized LTCO3 products in these northern-hemispheric regions. The OMI product is subject to large-scale interannual variability over the 2008–2017 decade, in comparison with which the underlying linear trends are small in absolute terms, with large confidence ranges (i.e. 95 % confidence intervals) intersecting with zero. However, the OMI LTCO3 product has been shown to be stable over this period relative to ozonesondes by Pope at el. (2023). IASI-FORLI has substantial negative LTCO3 tendencies, but this is driven by a specific discontinuity in 2010 due to an inhomogeneity in EUMETSAT (water vapour, temperature) data used in IASI-FORLI level-2 processing (Boynard et al., 2018; Wespes et al., 2018). This inhomogeneity induces an artificial drift that explains the substantial negative LTCO3 trends reported here and in Gaudel et al. (2018). The IASI-SOFRID LTCO3 and the a priori are very similar, with little interannual variability, which suggests that the IASI-SOFRID O3 retrieval in this height range is more constrained by the a priori (i.e. less TO3 sensitivity than the other products – see Sect. S3 in the Supplement). Importantly, analysis of the three products' a priori LTCO3 records show negligible trends, meaning that year-to-year spatio-temporal sampling differences (i.e. the number of retrievals used in the spatio-monthly regional averages) are not skewing long-term satellite trends. In summary: any underlying linear trend in LTCO3 occurring during the decade 2008–2017 was masked by interannual variability in the OMI retrieval and by constraint to the a priori in the IASI-SOFRID retrieval; although substantial for IASI-FORLI retrieval, this is due to changing meteorological inputs to the data processing (Boynard et al., 2018; Wespes et al., 2018).
For UKESM1.0, the model exhibits negligible long-term temporal variability in LTCO3 for all regions and for all the instruments' samplings. Modelled LTCO3 trends never exceeded 1.0 DU per decade in magnitude and were all deemed to be insignificant due to large associated p values using the linear–seasonal trend model detailed in Sect. 2.3 and Eqs. (2) and (3). The ozonesondes for each region were included to ground-truth the model and satellite trends. The North American sites' LTCO3 time series were relatively noisy and exhibited considerable interannual variability in their seasonal cycles. The comparatively low level of interannual variability in the European UKESM1.0 record of LTCO3 was in good agreement with the ozonesondes and so was its low trend, providing confidence in the model over that region. For East Asia, the interannual variability differed substantially between UKESM1.0 and ozonesondes, and the reported ozonesonde trend was significantly much larger than for UKESM1.0. Therefore, when considering UKESM1.0 and the ozonesondes, no consistent LTCO3 trends can be determined for any of the regions. Overall, taking all data sets into account, LTO3 appears to have neither increased nor decreased markedly over these three regions between the beginning and end of the study decade (i.e. 2008–2017).
One key aspect of this work was to exploit UKESM1.0 to determine the importance of vertical sensitivity for retrieved LTO3 and how this influences the reported long-term tendency. In terms of the absolute model trends (with and without the satellite AKs), the impact on LTCO3 was small, typically with near-zero tendencies and large uncertainty ranges (i.e. the 95 % confidence interval). In relative terms, the changes in model trend values were more substantial on the order of 100 %. To explore this further, the UKESM1.0 LTCO3 time series (with and without the satellite AKs) were sorted from lowest to highest (based on annual averages) to impose the most substantial trend in the model data. When the trends were re-calculated, the largest model LTCO3 trends ranged between 0.37 and 0.90 DU per decade. When the AKs were applied, the LTCO3 trends ranged from −0.74 to 1.02 DU per decade. Again, in relative terms, this represents a large impact of the AKs on simulated LTO3 tendencies, but in absolute terms, these are small changes. However, it should be noted that many of the 95% confidence intervals shifted between statistical regimes (i.e. a p value ≤ 0.05 vs. a p value > 0.05) once the AKs were applied to the model (i.e. the confidence interval initially did not intersect with zero but then did, or the interval did intersect with zero before shifting to a non-zero range).
Gaudel et al. (2018) suggested two potential reasons for the TCO3 trend discrepancies in their study:
-
time-varying instrument biases or drift
-
the impact of satellite vertical sensitivity.
A further two important reasons are as follows:
-
Changes over time in latitude–longitude domains sampled by satellite sensors (e.g. GOME-1 has substantial issues after 2003)
-
the time period used for the trend analysis.
As stated by Boynard et al. (2018) and Wespes et al. (2018), the IASI-FORLI-v20151001 product has an artificial negative drift with time, explained by a discontinuity found in the level-2 meteorological inputs taken from EUMETSAT. However, in the near-future, a new consistent IASI-FORLI ozone climate data record will be available using homogeneous level-1 and level-2 EUMETSAT meteorological data. The analysis of OMI LTCO3 by Pope et al. (2023) showed OMI LTCO3 to be temporally stable against ozonesondes. A similar analysis (not shown here) indicates that IASI-SOFRID LTCO3 is also temporally stable, with near-zero drift in bias. For the satellite vertical sensitivity, some of our results were unexpected. While the application of the AKs to UKESM1.0 can substantially shift the simulated absolute LTO3 values and squash or stretch the seasonal amplitude, the impacts on the simulation LTCO3 tendencies are small in absolute terms. In relative terms, the impacts can be large (e.g. 100 % change in trend rate). However, as the UKESM1.0 simulated LTCO3 trends are generally near zero, it is difficult to confidently say either way if the vertical sensitivity, when retrieving LTCO3, is important for influencing long-term tendencies, even when a more substantial trend is forced upon UKESM1.0. Future work on this would probably need to look at artificial model data which already contain substantial TO3 trends (e.g. 5.0 or 10.0 DU per decade). This will obviously not match reality but would provide some further quantification on how important vertical sensitivity is in relation to different instruments or sounders in LTO3 trend determination.
As for year-to-year spatio-temporal sampling, our results suggest negligible trends for the product LTCO3 a priori time series; thus, monthly sampling biases are unlikely to introduce artificial trends as the a priori data sets are trendless. Finally, the time period over which the trend analysis is undertaken is critically important. Gaudel et al. (2018), using the available data at the time, focussed on 2005–2015 and 2016 and 2008–2015 and 2016 for the OMI and IASI products they used. For the IASI products, using a slightly extended time period, the trends show similar tendencies. However, for OMI, 2016 and 2017 represent years of lower TO3. As a result, this dampens the strong significant positive trends in TCO3 reported by Gaudel et al. (2018). It is notable that the substantial positive increase in tropical LTO3 between 1995 and 2017 reported by Pope et al. (2023) from a series of UV–Vis sounders, including the same OMI global data set that is used here, further suggests that the selection of the time period and geographical region is crucial with regard to the role of interannual variability in linear trend detection.
Gaudel et al. (2018) undertook a multi-satellite analysis of long-term trends in tropospheric column ozone (TCO3). They found large-scale differences between these products with no clear consensus on the signs or drivers of these TCO3 trends. To avoid complications with regard to tropopause definition and to reduce the influence of stratospheric ozone on retrieved values, this study has undertaken a detailed follow-up assessment of decadal trends in LTCO3 (surface–450 hPa layer) rather than TCO3, exploiting ozonesonde records and model simulations and accounting carefully for satellite O3 metrics (e.g. averaging kernels (AKs), a priori information, and satellite uncertainties). We have focussed on LTCO3 data sets from the Ozone Monitoring Instrument (OMI) produced by the RAL Space scheme and from the Infrared Atmospheric Sounding Interferometer (IASI) produced by the IASI-FORLI and IASI-SOFRID schemes, for which there were consistent records from the period 2008–2017.
The evaluation of satellite LTO3 from these three products over the North American, European, and East Asian regions resulted in linear trends which varied over a small range close to zero and with confidence intervals intersecting with zero. This was consistent with simulations from the UK Earth System Model (UKESM1.0). There were no large-scale trends in the a priori information; thus, changes in satellite year-to-year spatio-temporal sampling have not been driving inconsistencies between products. When convolving UKESM1.0 with the satellite AKs (i.e. to assess the impact of the satellite vertical sensitivity), it did change the size of the model trend and, in some instances, the direction of the trend, but as the simulated LTO3 trends were small and insignificant, they had limited influence. Overall, our results show that changes in LTO3 during the decade 2008–2017 in North America, Europe, and East Asia were dominated by variability in processes which control LTO3 on shorter timescales.
In the near-future, the new European polar-orbiting mission Metop Second Generation will include IASI Next Generation and Sentinel-5 UV–VIS sounders to provide height-resolved ozone products to extend current missions through to the mid-2040s. This will be supplemented by the new USA Near Earth Orbit Network (NEON) series as a replacement for the Joint Polar Satellite System (JPSS). The Geostationary Environment Monitoring Spectrometer (GEMS) and Tropospheric Emissions: Monitoring of Pollution (TEMPO) have also recently been launched, and there will be new geostationary platforms: the Infrared Sounder (IRS) and Sentinel-4 UV–VIS sounder on Europe's Meteosat Third Generation (MTG-S), again through to the mid-2040s, and the USA Geostationary Extended Observations (GeoXO) series. Overall, these platforms will provide large volumes of data (e.g. diurnal observations) on tropospheric ozone over a long timescale for future regional trend analyses.
The IASI-FORLI and IASI-SOFRID data can be obtained from https://iasi.aeris-data.fr/O3 (IASI-SOFRID, 2022) and https://iasi-sofrid.sedoo.fr/ (IASI-FORLI, 2020). The RAL OMI data are available via the NERC Centre for Environmental Data Analysis (CEDA) through JASMIN platform, subject to data requests. However, the RAL Space OMI and UKESM1.0 LTCO3 data used for the analysis in this work can be found on the Zenodo open-access platform (https://doi.org/10.5281/zenodo.13342181; Pope, 2024). The ozonesonde data for WOUDC, SHADOZ, and NOAA are available from https://woudc.org/ (WOUDC, 2023), https://tropo.gsfc.nasa.gov/shadoz/ (SHADOZ, 2023), and https://gml.noaa.gov/ozwv/ozsondes/ (NOAA, 2023).
The supplement related to this article is available online at: https://doi.org/10.5194/acp-24-9177-2024-supplement.
RJP conceptualized, planned, and undertook the research study. BB, ELF, BJK, RS, BGL, AB, and CW provided the OMI and IASI ozone data and advice on using the products and their analysis. FMO'C and MD provided advice and expertise on using and running UKESM. CR provided advice and help during RP's ESA CCI fellowship. Scientific and technical contributions came from MPC, WF, MAP, SSD, and RR. RJP prepared the article with input from all the co-authors.
The contact author has declared that none of the authors has any competing interests.
Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. While Copernicus Publications makes every effort to include appropriate place names, the final responsibility lies with the authors.
This article is part of the special issue “Tropospheric Ozone Assessment Report Phase II (TOAR-II) Community Special Issue (ACP/AMT/BG/GMD inter-journal SI)”. It is a result of the Tropospheric Ozone Assessment Report, Phase II (TOAR-II, 2020–2024).
This work was funded by the UK Natural Environment Research Council (NERC) through the provision of funding for the National Centre for Earth Observation (NCEO, award reference no. NE/R016518/1), the NERC-funded UKESM Earth system modelling project (award no. reference NE/N017978/1), and the European Space Agency (ESA) Climate Change Initiative (CCI) post-doctoral fellowship scheme (award reference no. 4000137140). For UKESM1.0 model runs, we acknowledge the use of the Monsoon2 system, a collaborative facility supplied under the Joint Weather and Climate Research Programme, a strategic partnership between the Met Office and NERC. IASI is a joint mission of EUMETSAT and the Centre National d'Etudes Spatiales (CNES, France). The IASI-SOFRID research was conducted at LAERO, with some financial support from the CNES French spatial agency (TOSCA–IASI project). The authors acknowledge the AERIS data infrastructure for providing access to the IASI-FORLI data, ULB-LATMOS for the development of the FORLI retrieval algorithm, and the AC SAF project of EUMETSAT for providing the IASI-FORLI data used in this paper. Anna Maria Trofaier (ESA Climate Office) provided support and advice throughout the fellowship.
This research has been supported by the National Centre for Earth Observation (grant no. NE/R016518/1), the NERC-funded UKESM Earth system modelling project (award reference no. NE/N017978/1), and the European Space Agency (grant no. 4000137140).
This paper was edited by Jianzhong Ma and reviewed by Helen Worden and one anonymous referee.
Archibald, A. T., O'Connor, F. M., Abraham, N. L., Archer-Nicholls, S., Chipperfield, M. P., Dalvi, M., Folberth, G. A., Dennison, F., Dhomse, S. S., Griffiths, P. T., Hardacre, C., Hewitt, A. J., Hill, R. S., Johnson, C. E., Keeble, J., Köhler, M. O., Morgenstern, O., Mulcahy, J. P., Ordóñez, C., Pope, R. J., Rumbold, S. T., Russo, M. R., Savage, N. H., Sellar, A., Stringer, M., Turnock, S. T., Wild, O., and Zeng, G.: Description and evaluation of the UKCA stratosphere–troposphere chemistry scheme (StratTrop vn 1.0) implemented in UKESM1, Geosci. Model Dev., 13, 1223–1266, https://doi.org/10.5194/gmd-13-1223-2020, 2020.
Barnes, E. A., Fiore, A. M., and Horowitz, L. W.: Detection of trends in surface ozone in the presence of climate variability, J. Geophys. Res.-Atmos., 121, 6112–6129, 2016.
Barret, B., Emili, E., and Le Flochmoen, E.: A tropopause-related climatological a priori profile for IASI-SOFRID ozone retrievals: improvements and validation, Atmos. Meas. Tech., 13, 5237–5257, https://doi.org/10.5194/amt-13-5237-2020, 2020.
Boersma, K. F., Jacob, D. J., Eskes, H. J., Pinder, R. W., Wang, J., and van der A, R. J.: Intercomparison of SCIAMACHY and OMI tropospheric NO2 columns: Observing the diurnal evolution of chemistry and emissions from space, J. Geophys. Res.-Atmos., 113, D16S26, https://doi.org/10.1029/2007JD008816, 2008.
Boersma, K. F., Eskes, H. J., Dirksen, R. J., van der A, R. J., Veefkind, J. P., Stammes, P., Huijnen, V., Kleipool, Q. L., Sneep, M., Claas, J., Leitão, J., Richter, A., Zhou, Y., and Brunner, D.: An improved tropospheric NO2 column retrieval algorithm for the Ozone Monitoring Instrument, Atmos. Meas. Tech., 4, 1905–1928, https://doi.org/10.5194/amt-4-1905-2011, 2011.
Boynard, A., Clerbaux, C., Coheur, P.-F., Hurtmans, D., Turquety, S., George, M., Hadji-Lazaro, J., Keim, C., and Meyer-Arnek, J.: Measurements of total and tropospheric ozone from IASI: comparison with correlative satellite, ground-based and ozonesonde observations, Atmos. Chem. Phys., 9, 6255–6271, https://doi.org/10.5194/acp-9-6255-2009, 2009.
Boynard, A., Hurtmans, D., Garane, K., Goutail, F., Hadji-Lazaro, J., Koukouli, M. E., Wespes, C., Vigouroux, C., Keppens, A., Pommereau, J.-P., Pazmino, A., Balis, D., Loyola, D., Valks, P., Sussmann, R., Smale, D., Coheur, P.-F., and Clerbaux, C.: Validation of the IASI FORLI/EUMETSAT ozone products using satellite (GOME-2), ground-based (Brewer–Dobson, SAOZ, FTIR) and ozonesonde measurements, Atmos. Meas. Tech., 11, 5125–5152, https://doi.org/10.5194/amt-11-5125-2018, 2018.
Chang, K.-L., Cooper, O. R., Gaudel, A., Petropavlovskikh, I., and Thouret, V.: Statistical regularization for trend detection: an integrated approach for detecting long-term trends from sparse tropospheric ozone profiles, Atmos. Chem. Phys., 20, 9915–9938, https://doi.org/10.5194/acp-20-9915-2020, 2020.
Clerbaux, C., Boynard, A., Clarisse, L., George, M., Hadji-Lazaro, J., Herbin, H., Hurtmans, D., Pommier, M., Razavi, A., Turquety, S., Wespes, C., and Coheur, P.-F.: Monitoring of atmospheric composition using the thermal infrared IASI/MetOp sounder, Atmos. Chem. Phys., 9, 6041–6054, https://doi.org/10.5194/acp-9-6041-2009, 2009.
Doherty, R. M., Heal, M. R., and O'Connor, F. M.: Climate change impacts on human health over Europe through its effect on air quality, Environ. Health, 16, 33–44, https://doi.org/10.1186/s12940-017-0325-2, 2017.
ESA: Climate Change Initiative, http://cci.esa.int/ozone (last access: 1 September 2022), 2019.
Eskes, H. J. and Boersma, K. F.: Averaging kernels for DOAS total-column satellite retrievals, Atmos. Chem. Phys., 3, 1285–1291, https://doi.org/10.5194/acp-3-1285-2003, 2003.
European Commission: Joint Research Centre, Dentener F, et al. 2016, Hemispheric Transport of Air Pollution (HTAP): specification of the HTAP2 experiments: ensuring harmonized modelling, Publications Office, https://data.europa.eu/doi/10.2788/725244 (last access: 15 July 2022), 2016.
Eyring, V., Bony, S., Meehl, G. A., Senior, C. A., Stevens, B., Stouffer, R. J., and Taylor, K. E.: Overview of the Coupled Model Intercomparison Project Phase 6 (CMIP6) experimental design and organization, Geosci. Model Dev., 9, 1937–1958, https://doi.org/10.5194/gmd-9-1937-2016, 2016.
Gaudel, A., Cooper, O. R., Ancellet, G., Barret, B., Boynard, A., Burrows, J. P., Clerbaux, C., Coheur, P. F., Cuesta, J., Cuevas, E., Doniki, S., Dufour, G., Ebojie, F., Foret, G., Garia, O., GranadosMunoz, M. J., Hannigan, J. W., Hase, F., Hassler, B., Huang, G., Hurtmans, D., Jaffe, D., Jones, N., Kalabokas, P., Kerridge, B., Kulwaik, S., Latter, B., Leblanc, T., Le Flochmoen, E., Lin, W., Liu, J., Liu, X., Mahieu, E., McClure-Begley, A., Neu, J. L., Osman, M., Palm, M., Petetin, H., Petropavlovskikh, I., Querel, R., Rahpoe, N., Rozanov, A., Schultz, M. G., Schwab, J., Siddans, R., Smale, D., Steinbacher, M., Tanimoto, H., Tarasick, D. W., Thouret, V., Thompson, A. M., Trickl, T., Weatherhead, E., Wespes, C., Worden, H. M., Vigouroux, C., Xu, X., Zeng, G., and Ziemke, J.: Tropospheric Ozone Assessment Report: Present day distribution and trends of tropospheric ozone relevant to climate and global atmospheric chemistry model evaluation, Elementa, 6, 1–58, https://doi.org/10.1525/elementa.291, 2018.
Fiore, F. M, Hancock, S. E., Lamarque, J.-F., Correa, G. P., Chang, K.-L., Ru, M., Cooper, O., Gaudel, A., Polvani, L. M., Sauvage, B., and Ziemke, J. R.: Understanding recent tropospheric ozone trends in the context of large internal variability: A new perspective from chemistry-climate model ensembles, Environ. Res. Clim., 1, 025008, https://doi.org/10.1088/2752-5295/ac9cc2, 2022.
Forster, P., Storelvmo, T., Armour, K., Collins, W., Dufresne, J.- L., Frame, D., Lunt, D. J., Mauritsen, T., Palmer, M. D., Watanabe, M., Wild, M., and Zhang, H.: The Earth's Energy Budget, Climate Feedbacks, and Climate Sensitivity, in: Climate Change 2021: The Physical Science Basis, Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change, edited by: Masson-Delmotte, V., Zhai, P., Pirani, A., Connors, S. L., Péan, C., Berger, S., Caud, N., Chen, Y., Goldfarb, L., Gomis, M. I., Huang, M., Leitzell, K., Lonnoy, E., Matthews, J. B. R.,Maycock, T. K., Waterfield, T., Yelekçi, O., Yu, R., and Zhou, B., Cambridge University Press, Cambridge, United Kingdom and New York, NY, USA, 923–1054, https://doi.org/10.1017/9781009157896.009, 2021.
Gaudel, A., Cooper, O. R., Ancellet, G., Barret, B., Boynard, A., Burrows, J. P., Clerbaux, C., Coheur, P. F., Cuesta, J., Cuevas, E., Doniki, S., Dufour, G., Ebojie, F., Foret, G., Garia, O., Granados-Munoz, M. J., Hannigan, J. W., Hase, F., Hassler, B., Huang, G., Hurtmans, D., Jaffe, D., Jones, N., Kalabokas, P., Kerridge, B., Kulwaik, S., Latter, B., Leblanc, T., Le Flochmoen, E., Lin, W., Liu, J., Liu, X., Mahieu, E., McClure-Begley, A., Neu, J. L., Osman, M., Palm, M., Petetin, H., Petropavlovskikh, I., Querel, R., Rahpoe, N., Rozanov, A., Schultz, M.G., Schwab, J., Siddans, R., Smale, D., Steinbacher, M., Tanimoto, H., Tarasick, D.W., Thouret, V., Thompson, A. M., Trickl, T., Weatherhead, E., Wespes, C., Worden, H. M., Vigouroux, C., Xu, X., Zeng, G., and Ziemke, J.: Tropospheric Ozone Assessment Report: Present day distribution and trends of tropospheric ozone relevant to climate and global atmospheric chemistry model evaluation, Elementa, 6, 1–58, https://doi.org/10.1525/elementa.291, 2018.
Gauss, M., Myhre, G., Isaksen, I. S. A., Grewe, V., Pitari, G., Wild, O., Collins, W. J., Dentener, F. J., Ellingsen, K., Gohar, L. K., Hauglustaine, D. A., Iachetti, D., Lamarque, F., Mancini, E., Mickley, L. J., Prather, M. J., Pyle, J. A., Sanderson, M. G., Shine, K. P., Stevenson, D. S., Sudo, K., Szopa, S., and Zeng, G.: Radiative forcing since preindustrial times due to ozone change in the troposphere and the lower stratosphere, Atmos. Chem. Phys., 6, 575–599, https://doi.org/10.5194/acp-6-575-2006, 2006.
Hoesly, R. M., Smith, S. J., Feng, L., Klimont, Z., Janssens-Maenhout, G., Pitkanen, T., Seibert, J. J., Vu, L., Andres, R. J., Bolt, R. M., Bond, T. C., Dawidowski, L., Kholod, N., Kurokawa, J.-I., Li, M., Liu, L., Lu, Z., Moura, M. C. P., O'Rourke, P. R., and Zhang, Q.: Historical (1750–2014) anthropogenic emissions of reactive gases and aerosols from the Community Emissions Data System (CEDS), Geosci. Model Dev., 11, 369–408, https://doi.org/10.5194/gmd-11-369-2018, 2018.
Hollaway, M. J., Arnold, S. R., Challinor, A. J., and Emberson, L. D.: Intercontinental trans-boundary contributions to ozone-induced crop yield losses in the Northern Hemisphere, Biogeosciences, 9, 271–292, https://doi.org/10.5194/bg-9-271-2012, 2012.
IASI-SOFRID: Welcome to the IASI-SOFRID database (vn3.5) [data set], http://thredds.sedoo.fr/iasi-sofrid-o3-co/ (last access: 1 December 2022), 2022.
IASI-FORLI: Daily IASI/Metop-A ULB-LATMOS ozone (O3) L2 product (total column and vertical profile) (v20151001) [data set], https://iasi.aeris-data.fr/catalog/ (last access: 15 December 2022), 2020.
Illingworth, S. M., Remedios, J. J., Boesch, H., Moore, D. P., Sembhi, H., Dudhia, A., and Walker, J. C.: ULIRS, an optimal estimation retrieval scheme for carbon monoxide using IASI spectral radiances: sensitivity analysis, error budget and simulations, Atmos. Meas. Tech., 4, 269–288, https://doi.org/10.5194/amt-4-269-2011, 2011.
Keim, C., Eremenko, M., Orphal, J., Dufour, G., Flaud, J.-M., Höpfner, M., Boynard, A., Clerbaux, C., Payan, S., Coheur, P.-F., Hurtmans, D., Claude, H., Dier, H., Johnson, B., Kelder, H., Kivi, R., Koide, T., López Bartolomé, M., Lambkin, K., Moore, D., Schmidlin, F. J., and Stübi, R.: Tropospheric ozone from IASI: comparison of different inversion algorithms and validation with ozone sondes in the northern middle latitudes, Atmos. Chem. Phys., 9, 9329–9347, https://doi.org/10.5194/acp-9-9329-2009, 2009.
Keppens, A., Compernolle, S., Verhoelst, T., Hubert, D., and Lambert, J.-C.: Harmonization and comparison of vertically resolved atmospheric state observations: methods, effects, and uncertainty budget, Atmos. Meas. Tech., 12, 4379–439, https://doi.org/10.5194/amt-12-4379-2019, 2019.
Lamarque, J.-F., Bond, T. C., Eyring, V., Granier, C., Heil, A., Klimont, Z., Lee, D., Liousse, C., Mieville, A., Owen, B., Schultz, M. G., Shindell, D., Smith, S. J., Stehfest, E., Van Aardenne, J., Cooper, O. R., Kainuma, M., Mahowald, N., McConnell, J. R., Naik, V., Riahi, K., and van Vuuren, D. P.: Historical (1850–2000) gridded anthropogenic and biomass burning emissions of reactive gases and aerosols: methodology and application, Atmos. Chem. Phys., 10, 7017–7039, https://doi.org/10.5194/acp-10-7017-2010, 2010.
Matthes, K., Funke, B., Andersson, M. E., Barnard, L., Beer, J., Charbonneau, P., Clilverd, M. A., Dudok de Wit, T., Haberreiter, M., Hendry, A., Jackman, C. H., Kretzschmar, M., Kruschke, T., Kunze, M., Langematz, U., Marsh, D. R., Maycock, A. C., Misios, S., Rodger, C. J., Scaife, A. A., Seppälä, A., Shangguan, M., Sinnhuber, M., Tourpali, K., Usoskin, I., van de Kamp, M., Verronen, P. T., and Versick, S.: Solar forcing for CMIP6 (v3.2), Geosci. Model Dev., 10, 2247–2302, https://doi.org/10.5194/gmd-10-2247-2017, 2017.
McPeters, R. D., Labow, G. J., and Logan, J. A.: Ozone climatological profiles for satellite retrieval algorithms, J. Geophys. Res., 112, D05308, https://doi.org/10.1029/2005JD006823, 2007.
Miles, G. M., Siddans, R., Kerridge, B. J., Latter, B. G., and Richards, N. A. D.: Tropospheric ozone and ozone profiles retrieved from GOME-2 and their validation, Atmos. Meas. Tech., 8, 385–398, https://doi.org/10.5194/amt-8-385-2015, 2015.
Myhre, G., Shindell, D., Bréon, F.-M, Collins, W., Fuglestvedt, J., Huang, J., Koch, D., Lamarque, J.-F., Lee, D., Mendoza, B., Nakajima, T., Robock, A., Stephens, G., Takemura, T., and Zhang, H.: Anthropogenic and Natural Radiative Forcing, in: Climate Change 2013: The Physical Science Basis, Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, Cambridge University Press, Cambridge, United Kingdom and New York, NY, USA, 659–740, https://doi.org/10.1017/CBO9781107415324.018., 2013.
NOAA: ESRL/GML Ozonesondes [data set], https://gml.noaa.gov/ozwv/ozsondes/ (last access: 1 June 2022), 2023.
O'Connor, F. M., Johnson, C. E., Morgenstern, O., Abraham, N. L., Braesicke, P., Dalvi, M., Folberth, G. A., Sanderson, M. G., Telford, P. J., Voulgarakis, A., Young, P. J., Zeng, G., Collins, W. J., and Pyle, J. A.: Evaluation of the new UKCA climate-composition model – Part 2: The Troposphere, Geosci. Model Dev., 7, 41–91, https://doi.org/10.5194/gmd-7-41-2014, 2014.
Pimlott, M. A., Pope, R. J., Kerridge, B. J., Latter, B. G., Knappett, D. S., Heard, D. E., Ventress, L. J., Siddans, R., Feng, W., and Chipperfield, M. P.: Investigating the global OH radical distribution using steady-state approximations and satellite data, Atmos. Chem. Phys., 22, 10467–10488, https://doi.org/10.5194/acp-22-10467-2022, 2022.
Pope, R.: Investigation of satellite vertical sensitivity on long-term retrieved lower tropospheric ozone trends, in: Atmospheric Chemistry and Physics, Zenodo [data set], https://doi.org/10.5281/zenodo.13342181, 2024.
Pope, R. J., Arnold, S. R., Chipperfield, M. P., Latter, B. G., Siddans, R., and Kerridge, B. J.: Widespread changes in UK air quality observed from space, Atmos. Sci. Lett., 19, e817, https://doi.org/10.1002/asl.817, 2018.
Pope, R. J., Kerridge, B. J., Siddans, R., Latter, B. G., Chipperfield, M. P., Feng, W., Pimlott, M. A., Dhomse, S. S., Retscher, C., and Rigby, R.: Investigation of spatial and temporal variability in lower tropospheric ozone from RAL Space UV–Vis satellite products, Atmos. Chem. Phys., 23, 14933–14947, https://doi.org/10.5194/acp-23-14933-2023, 2023.
Pope, R. J., Rap, A., Pimlott, M. A., Barret, B., Le Flochmoen, E., Kerridge, B. J., Siddans, R., Latter, B. G., Ventress, L. J., Boynard, A., Retscher, C., Feng, W., Rigby, R., Dhomse, S. S, Wespes, C. and Chipperfield, M. P.: Quantifying the tropospheric ozone radiative effect and its temporal evolution in the satellite era, Atmos. Chem. Phys., 24, 3613–3626, https://doi.org/10.5194/acp-24-3613-2024, 2024.
Rao, S., Klimont, Z., Smith, S. J., Van Dingenen, R., Dentener, F., Bouwman, L., Riahi, K., Amann, M., Bodirsky, B. L., van Vuuren, D.P., Reus, L. R., Calvin, K., Drouet, L., Fricko, O., Fujimori, S., Gernaat, D., Havlik, P., Harmsen, M., Hasegawa, T., Heyes, C., and Tavoni, M.: Future air pollution in the shared socio-economic pathways, Glob. Environ. Change, 42, 346–358, https://doi.org/10.1016/j.gloenvcha.2016.05.012, 2017.
Richards, N. A. D, Osterman, G. B., Browell, E. V., Hair, J. W., Avery, M., and Li, Q.: Validation of tropospheric emission spectrometer ozone profiles with aircraft observations during the intercontinental chemical transport experiment – B, J. Geophys. Res., 113, D16S29, https://doi.org/10.1029/2007JD008815, 2008.
Russo, M. R., Kerridge, B. J., Abraham, N. L., Keeble, J., Latter, B. G., Siddans, R., Weber, J., Griffiths, P. T., Pyle, J. A., and Archibald, A. T.: Seasonal, interannual and decadal variability of tropospheric ozone in the North Atlantic: comparison of UM-UKCA and remote sensing observations for 2005–2018, Atmos. Chem. Phys., 23, 6169–6196, https://doi.org/10.5194/acp-23-6169-2023, 2023.
Sellar, A. A., Jones, C. G., Mulcahy, J. P., Tang, Y., Yool, A., Wiltshire, A., O'Connor, F. M., Stringer, M., Hill, R., Palmieri, J., Woodward, S., de Mora, L., Kuhlbrodt, T., Rumbold, S. T., Kelley, D. I., Ellis, R., John, C. E., Walton, J., Abraham, N. L., Andrews, M. B., Andrews, T., Archibald, A. T., Berthou, S., Burke, E., Blockley, E., Carslaw, K., Dalvi, M., Edwards, J., Folbert, G. A., Gedney, N., Griffiths, P. T., Harper, A. B., Hendry, M. A., Hewitt, A. J., Johnson, B., Jones, A., Jones, C. D., Keebie, J., Liddicoat, S., Morgenstern, O., Parker, R. J., Predoi, V., Robertson, E., Siahaan, A., Smith, R. S., Swaminathan, R., Woodhouse, M. T., Zeng, G., and Zerroukat, M.: Description and Evaluation of the UK Earth System Model, J. Adv. Model. Earth Sy., 11, 4513–4558, https://doi.org/10.1029/2019MS001739, 2019.
SHADOZ: SHADOZ Data Archive [data set], https://tropo.gsfc.nasa.gov/shadoz/Archive.html (last access: 1 June 2022), 2023.
Sindelarova, K., Granier, C., Bouarar, I., Guenther, A., Tilmes, S., Stavrakou, T., Müller, J.-F., Kuhn, U., Stefani, P., and Knorr, W.: Global data set of biogenic VOC emissions calculated by the MEGAN model over the last 30 years, Atmos. Chem. Phys., 14, 9317–9341, https://doi.org/10.5194/acp-14-9317-2014, 2014.
Sitch, S., Cox, P. M., Collins, W. J., and Huntingford, C.: Indirect radiative forcing of climate change through ozone effects on the land carbon sink, Nature, 448, 791–795, https://doi.org/10.1038/nature06059, 2007.
Sofieva, V. F., Tamminen, J., Kyrölä, E., Mielonen, T., Veefkind, P., Hassler, B., and Bodeker, G. E.: A novel tropopause-related climatology of ozone profiles, Atmos. Chem. Phys., 14, 283–299, https://doi.org/10.5194/acp-14-283-2014, 2014.
Telford, P. J., Braesicke, P., Morgenstern, O., and Pyle, J. A.: Technical Note: Description and assessment of a nudged version of the new dynamics Unified Model, Atmos. Chem. Phys., 8, 1701–1712, https://doi.org/10.5194/acp-8-1701-2008, 2008.
van der A, R. J., Peters, D. H. M. U., Eskes, H., Boersma, K. F., Van Roozendael, M., De Smedt, I., and Kelder, H. M.: Detection of the trend and seasonal variation in tropospheric NO2 over China, J. Geophys. Res., 11, D12317, https://doi.org/10.1029/2005JD006594, 2006.
van der A, R. J., Eskes, H. J., Boersma, K. F., van Noije, T. P. C., Van Roozendael, M., De Smedt, I., Peters, D. H. M. U., and Meijer, E. W.: Trends, seasonal variability and dominant NOx sources derived from a ten year record of NO2 measured from space, J. Geophys. Res., 113, 1–12, https://doi.org/10.1029/2007JD009021, 2008.
van Marle, M. J. E., Kloster, S., Magi, B. I., Marlon, J. R., Daniau, A.-L., Field, R. D., Arneth, A., Forrest, M., Hantson, S., Kehrwald, N. M., Knorr, W., Lasslop, G., Li, F., Mangeon, S., Yue, C., Kaiser, J. W., and van der Werf, G. R.: Historic global biomass burning emissions for CMIP6 (BB4CMIP) based on merging satellite observations with proxies and fire models (1750–2015), Geosci. Model Dev., 10, 3329–3357, https://doi.org/10.5194/gmd-10-3329-2017, 2017.
Weatherhead, E. C., Reinsel, G. C., Tiao, G. C., Meng, X., Choi, D., Cheang, W., Keller, T., De Luisi, J., Wuebbles, D. J., Kerr, J. B., Miller, A. J., Oltmans, S. J., and Frederick, J. E.: Factors affecting the detection of trends: statistical considerations and applications to environmental data, J. Geophys. Res., 103, 17149–17161, https://doi.org/10.1029/98JD00995, 1998.
Wespes, C., Hurtmans, D., Clerbaux, C., Boynard, A., and Coheur, P.-F.: Decrease in tropospheric O3 levels in the Northern Hemisphere observed by IASI, Atmos. Chem. Phys., 18, 6867–6885, https://doi.org/10.5194/acp-18-6867-2018, 2018.
WHO (World Health Organisation): Ambient (outdoor) air pollution, https://www.who.int/news-room/fact-sheets/detail/ambient-(outdoor)-air-quality-and-health (last access: 2 January 2024), 2022.
Wiltshire, A. J., Burke, E. J., Chadburn, S. E., Jones, C. D., Cox, P. M., Davies-Barnard, T., Friedlingstein, P., Harper, A. B., Liddicoat, S., Sitch, S., and Zaehle, S.: JULES-CN: a coupled terrestrial carbon–nitrogen scheme (JULES vn5.1), Geosci. Model Dev., 14, 2161–2186, https://doi.org/10.5194/gmd-14-2161-2021, 2021.
WMO: Meteorology – A three-dimensional science, World Meteorological Organisation, Bulletin, 6, 134–138, 1957.
WOUDC: Data Search/Download [data set], https://woudc.org/data/explore.php (last access: 1 June 2022), 2023.
Young, P. J., Archibald, A. T., Bowman, K. W., Lamarque, J.-F., Naik, V., Stevenson, D. S., Tilmes, S., Voulgarakis, A., Wild, O., Bergmann, D., Cameron-Smith, P., Cionni, I., Collins, W. J., Dalsøren, S. B., Doherty, R. M., Eyring, V., Faluvegi, G., Horowitz, L. W., Josse, B., Lee, Y. H., MacKenzie, I. A., Nagashima, T., Plummer, D. A., Righi, M., Rumbold, S. T., Skeie, R. B., Shindell, D. T., Strode, S. A., Sudo, K., Szopa, S., and Zeng, G.: Pre-industrial to end 21st century projections of tropospheric ozone from the Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP), Atmos. Chem. Phys., 13, 2063–2090, https://doi.org/10.5194/acp-13-2063-2013, 2013.
Ziemke, J. R., Chandra, S., Labow, G. J., Bhartia, P. K., Froidevaux, L., and Witte, J. C.: A global climatology of tropospheric and stratospheric ozone derived from Aura OMI and MLS measurements, Atmos. Chem. Phys., 11, 9237–9251, https://doi.org/10.5194/acp-11-9237-2011, 2011.