Journal cover Journal topic
Atmospheric Chemistry and Physics An interactive open-access journal of the European Geosciences Union
Journal topic
Atmos. Chem. Phys., 18, 6867–6885, 2018
https://doi.org/10.5194/acp-18-6867-2018

Special issue: Quadrennial Ozone Symposium 2016 – Status and trends...

Atmos. Chem. Phys., 18, 6867–6885, 2018
https://doi.org/10.5194/acp-18-6867-2018

Research article 16 May 2018

Research article | 16 May 2018

# Decrease in tropospheric O3 levels in the Northern Hemisphere observed by IASI

Decrease in tropospheric O3 levels in the Northern Hemisphere observed by IASI
Catherine Wespes1, Daniel Hurtmans1, Cathy Clerbaux1,2, Anne Boynard2, and Pierre-François Coheur1 Catherine Wespes et al.
• 1Spectroscopie de l'Atmosphère, Service de Chimie Quantique et Photophysique, Faculté des Sciences, Université Libre de Bruxelles (ULB), Brussels, Belgium
• 2LATMOS/IPSL, UPMC Univ. Paris 06 Sorbonne Universités, UVSQ, CNRS, Paris, France

Correspondence: Catherine Wespes (cwespes@ulb.ac.be)

Abstract

In this study, we describe the recent changes in the tropospheric ozone (O3) columns measured by the Infrared Atmospheric Sounding Interferometer (IASI), onboard the Metop satellite, during the first 9 years of operation (January 2008 to May 2017). Using appropriate multivariate regression methods, we differentiate significant linear trends from other sources of O3 variations captured by IASI. The geographical patterns of the adjusted O3 trends are provided and discussed on the global scale. Given the large contribution of the natural variability in comparison with that of the trend (25–85 % vs. 15–50 %, respectively) to the total O3 variations, we estimate that additional years of IASI measurements are generally required to detect the estimated O3 trends with high precision. Globally, additional 6 months to 6 years of measurements, depending on the regions and the seasons, are needed to detect a trend of $|\mathrm{5}|$ DU decade−1. An exception is interestingly found during summer at mid- and high latitudes of the Northern Hemisphere (NH;  40 to  75 N), where the large absolute fitted trend values ($|\mathrm{0.5}|$ DU yr−1 on average) combined with the small model residuals ( 10 %) allow for detection of a band-like pattern of significant negative trends. Despite no consensus in terms of tropospheric O3 trends having been reached from the available independent datasets (UV or IR satellites, O3 sondes, aircrafts, ground-based measurements, etc.) for the reasons that are discussed in the text, this finding is consistent with the reported decrease in O3 precursor emissions in recent years, especially in Europe and USA. The influence of continental pollution on that latitudinal band is further investigated and supported by the analysis of the O3–CO relationship (in terms of correlation coefficient, regression slope and covariance) that we found to be the strongest at northern midlatitudes in summer.

1 Introduction

O3 plays a key role throughout the whole troposphere where it is produced by the photochemical oxidation of carbon monoxide (CO), non-methane volatile organic compounds (NMVOCs) and methane (CH4) in the presence of nitrogen oxides (NOx) (e.g., Logan et al., 1981). O3 sources in the troposphere consist of the in situ photochemical production from anthropogenic and natural precursors and the downwards transport of stratospheric O3. Being a strong pollutant, a major reactive species and an important greenhouse gas in the upper troposphere, O3 is of greatest interest regarding air quality, atmospheric chemistry and radiative forcing studies. Thanks to its long lifetime (several weeks) relative to transport timescales in the free troposphere (Fusco and Logan, 2003), O3 also contributes to large-scale transport of pollution far from source regions with further impacts on global air quality (e.g., Stohl et al., 2002; Parrish et al., 2012) and climate. Monitoring and understanding the time evolution of tropospheric O3 at a global scale is, therefore, crucial to apprehend future climate changes. Nevertheless, a series of limitations make O3 trends particularly challenging to retrieve and to interpret.

While the O3 precursors anthropogenic emissions have increased and shifted equatorward in the developing countries (Zhang et al., 2016) since the 1980s, extensive campaigns and routine in situ and remote measurements at specific urban and rural sites have provided long-term but sparse datasets of tropospheric O3 (e.g., Cooper et al., 2014 and references therein). Ultraviolet and visible (UV/VIS) atmospheric sounders onboard satellites provide tropospheric O3 measurements with much wider coverage, but they result either from indirect methods (e.g., Fishman et al., 2005) or from direct retrievals which are limited by coarse vertical resolution (Liu et al., 2010). All these datasets also suffer from a lack of homogeneity in terms of measurement methods (instrument and algorithm) and spatio-temporal samplings (e.g., Doughty et al., 2011; Heue et al., 2016; Leventidou et al., 2017). Those limitations, in addition to the large natural interannual variability (IAV) and decadal variations in tropospheric O3 levels (due to large-scale dynamical modes of O3 variations and to changes in stratospheric O3, in stratosphere–troposphere exchanges, in precursor emissions and in their geographical patterns), introduce strong biases in trends determined from independent studies and datasets (e.g., Zbinden et al., 2006; Thouret et al., 2006; Logan et al., 2012; Parrish et al., 2012 and references therein). As a consequence, determining accurate trends requires a long period of high density and homogeneous measurements (e.g., Payne et al., 2017).

Such long-term datasets are now becoming obtainable with the new generation of nadir-looking and polar-orbiting instruments measuring in the thermal infrared region. In particular, about one decade of O3 profile measurements, with good sensitivity in the troposphere independent from the layers above, is now available from the IASI (Infrared Atmospheric Sounding Interferometer) sounder aboard the European Metop platforms, allowing the monitoring of regional and global variations in tropospheric O3 levels (e.g., Dufour et al., 2012; Safieddine et al., 2014; Wespes et al., 2016).

In this study, we examine the tropospheric O3 changes behind the natural IAV as measured by IASI over January 2008–May 2017. To that end, we use the approach described in Wespes et al. (2017), which relies on a multi-linear regression (MLR) procedure, for accurately differentiating trends from other sources of O3 variations, the latter being robustly identified and quantified in the companion study. In Sect. 2, we briefly review the IASI mission and the tropospheric O3 product, and we shortly describe the multivariate models (annual or seasonal) that we use for fitting the daily O3 time series. In Sect. 3, after verifying the performance of the MLR models over the available IASI dataset, we evaluate the feasibility of capturing and retrieving significant trend parameters, apart from natural O3 dependencies, by performing trend sensitivity studies. In Sect. 4, we present and discuss the global distributions of the O3 trends estimated from IASI in the troposphere. The focus is given in summer over and downwind anthropogenic polluted areas of the NH where the possibility of inferring significant trends from the first  9 years of available IASI measurements is demonstrated. Finally, the O3–CO correlations, enhancement ratios and covariance are examined for characterizing the origin of the air masses in regions of positive and negative trends.

2 IASI O3 measurements and multivariate regression

The IASI instrument is a nadir-viewing Fourier transform spectrometer that records the thermal infrared emission of the Earth–atmosphere system between 645 and 2760 cm−1 from the polar Sun-synchronous orbiting meteorological Metop series of satellites. Metop-A and -B were successively launched in October 2006 and September 2012. The third and last launch, of Metop-C, is planned for 2018, which will ensure homogeneous long-term IASI measurements. The measurements are taken every 50 km along the track of the satellite at nadir and over a swath of 2200 km across track, with a field of view of four simultaneous footprints of 12 km at nadir, which provides global coverage of the Earth twice a day (at 09:30 and 21:30 mean local solar time). The instrument has good spectral resolution and low radiometric noise, which allows for retrieval of numerous gas-phase species in the troposphere (e.g., Clerbaux et al., 2009, and references therein; Hilton et al., 2012; Clarisse et al., 2011).

In this paper, we use the FORLI O3 profiles (Fast Optimal Retrievals on Layers for IASI processing chain set up at ULB; v20151001) retrieved from the IASI-A (aboard Metop-A) daytime measurements (defined with a solar zenith angle to the sun of < 80), which are characterized by good spectral fit (determined here by quality flags on biased or sloped residuals, suspect averaging kernels, maximum number of iteration exceeded, etc.) and which correspond to clear or almost-clear scenes (a filter based on a fractional cloud cover below 13 % has been applied; see Clerbaux et al., 2009; Hurtmans et al., 2012). These profiles are characterized by good vertical sensitivity in the troposphere and the stratosphere (e.g., Wespes et al., 2017). The FORLI algorithm relies on a fast radiative transfer and retrieval methodology based on the optimal estimation method (Rodgers, 2000) and is fully described in Hurtmans et al. (2012). The FORLI O3 profiles, which are retrieved at 40 constant vertical layers from surface up to 40 km and an additional 40–60 km one, have already undergone thorough characterization and validation exercises (e.g., Anton et al., 2011; Dufour et al., 2012; Gazeaux et al., 2013; Hurtmans et al., 2012; Parrington et al., 2012; Pommier et al., 2012; Scannell et al., 2012; Oetjen et al., 2014; Boynard et al., 2016, 2018; Wespes et al., 2016; Keppens et al., 2018). They demonstrated a good degree of accuracy, precision and vertical sensitivity, and no instrumental drift, for capturing large-scale dynamical modes of O3 variability in the troposphere independently from the layers above (Wespes et al., 2017), with the possibility of further differentiating long-term O3 changes in the troposphere (Wespes et al., 2016). Note, however, that a drift in the NH middle–low troposphere (MLT) O3 over the whole IASI dataset is reported in Keppens et al. (2018) and Boynard et al. (2018) from comparison with O3 sondes. This drift ( 2.8 DU decade−1 in the NH) is shown in Boynard et al. (2018) to result from a discontinuity (called “jump” by the author) in September 2010 in the IASI O3 time series, for reasons that are unclear at present. Furthermore, the drift strongly decreases (< $|\mathrm{1}|$ DU decade−1 on average) after the jump and it becomes even non-significant for most of the stations (significant positive drift is also found for some stations) over the periods before or after the jump, separately.

For the purpose of this work, we focus on a tropospheric column ranging from ground to 300 hPa (MLT) that includes the altitude of maximum sensitivity of IASI in the troposphere (usually between 4 and 8 km altitude), which limits the influences of the stratospheric O3 as much as possible and which was shown in Wespes et al. (2017) to exhibit independent deseasonalized anomalies/dynamical processes from those in the stratospheric layers. The stratospheric contribution into the tropospheric O3 columns have been previously estimated in Wespes et al. (2016) as ranging between 30 and 65 % depending on the region and the season with the smallest contribution and the largest sensitivity in the northern midlatitudes in spring–summer, where the O3 variations, thus, mainly originate from the troposphere. We use almost the same MLR model (in its annual or its seasonal formulation) as the one developed in the companion paper (see Eqs. 1 and 2; Sect. 2.2 in Wespes et al., 2017), which includes a series of geophysical variables in addition to a linear trend (LT) term. In order to take account of the observed jump properly, we modified the previously used MLR model so that the constant term is split into two components covering the periods before and after the September 2010 jump, separately. The MLR which is performed using an iterative stepwise backward elimination approach to retain the most relevant explanatory variables (called “proxies”) at the end of the iterations (e.g. Mäder et al., 2007) is applied to the daily IASI O3 time series. The main proxies selected to account for the natural variations in O3 are namely the QBO (quasi-biennial oscillation), the NAO (North Atlantic Oscillation) and the ENSO (El Niño–Southern Oscillation; see Table 1 in Wespes et al., 2017 for the exhaustive list of the proxies used). Their associated standard errors are estimated from the covariance matrix of the regression coefficients and are corrected to take into account the uncertainty due to the autocorrelation of the noise residual (see Eq. 3 in Wespes et al., 2016). The common rule that the regression coefficients are significant if they are greater in magnitude than 2 times their standard errors is applied (95 % confidence limits defined by 2σ level). The MLR model was found to give a good representation of the IASI O3 records in the troposphere over 2008–2016, allowing us to identify/quantify the main O3 drivers with marked regional differences in the regression coefficients. Time lags of 2 and 4 months for ENSO are also included hereafter in the MLR model to account for a large but delayed impact of ENSO on mid- and high latitudes O3 variations far from the equatorial Pacific, where the ENSO signal originates (Wespes et al., 2017).

Figure 1Seasonal distribution of O3 tropospheric columns (in DU, integrated from ground to 300 hPa) measured by IASI and averaged over January 2008–May 2017 (a), of the RMSE of the regression fits (in DU, b) and of the fraction of the variation in IASI data explained by the regression model, calculated as $\mathrm{100}×\left(\mathit{\sigma }\left({{\mathbit{O}}_{\mathbf{3}}}^{\text{Fitted_ model}}\left(t\right)\right)/\mathit{\sigma }\left({\mathbit{O}}_{\mathbf{3}}\left(t\right)\right)\right)$ (in percent, c). Data are averaged over a 2.5× 2.5 grid box.

3 Regression performance and sensitivity to trend

In this section, we first verify the performance of the MLR models (annual and seasonal; in terms of residual errors and variation explained by the model) to globally reproduce the time evolution of O3 records over the entire studied period (January 2008–May 2017). Based on this, we then investigate the statistical relevance of a trend study from IASI measurements in the troposphere by examining the sensitivity of the pair IASI-MLR to the retrieved LT term.

Figure 1 presents the seasonal distributions of tropospheric O3 measured by IASI averaged over January 2008–May 2017 (Fig. 1a), along with the RMSE of the seasonal regression fit (in DU; Fig. 1b) and the contribution of the fitted seasonal model into the IASI O3 time series (in percent; Fig. 1c), calculated as $\frac{\mathit{\sigma }\left({{\mathbit{O}}_{\mathbf{3}}}^{\mathrm{Fitted}\mathrm{_}\mathrm{model}}\left(t\right)\right)}{\mathit{\sigma }\left({\mathbit{O}}_{\mathbf{3}}\left(t\right)\right)}$, where σ is the standard deviation relative to the regression models and to the IASI O3 time series. These two statistical parameters help to evaluate how well the fitted model explains the variability in the IASI O3 observations. The seasonal patterns of O3 measurements are close to those reported in Wespes et al. (2017) for a shorter period (see Sects. 2.1 and 3.1 in Wespes et al., 2017 for a detailed description of the distributions) and they clearly show, for instance, high O3 values over the highly populated areas of Asia in summer. The distributions from Fig. 1 show that the model reproduces between 35 and 90 % of the daily O3 variation captured by IASI and that the residual errors vary between 0.01 and 5 DU (i.e., the RMSEs relative to the IASI O3 time series are of  15 % on global average and vary between 10 % in the NH in summer and 30 % in specific tropical regions). On an annual basis (data not shown), the model explains a large fraction of the variation in the IASI O3 dataset (from  45 to  85 %) and the RMSEs are lower than 4.5 DU everywhere ( 3 DU on the global average). The relative RMSE is less than 1 % in almost all situations indicating the absence of bias.

Figure 2Seasonal distributions of the contribution from the seasonal and explanatory variables into the IASI O3 variations estimated as $\left[\mathrm{100}×\mathit{\sigma }\left(\sum _{n=\mathrm{1};j=\mathrm{2}}^{\mathrm{4};m}\left[{a}_{n};{b}_{n};{x}_{j}\right]\left[\mathrm{cos}\left(n\mathit{\omega }t\right);\mathrm{sin}\left(n\mathit{\omega }t\right);{\mathbit{X}}_{\mathrm{norm},\mathbit{j}}\left(t\right)\right]\right)/\mathit{\sigma }\left({\mathbit{O}}_{\mathbf{3}}\left(t\right)\right)\right]$ (in percent, a) and of the contribution from the linear trend calculated as $\left[\mathrm{100}×\mathit{\sigma }\left({x}_{j=\mathrm{1}}\cdot \mathit{\text{trend}}\right)/\mathit{\sigma }\left({\mathbit{O}}_{\mathbf{3}}\left(t\right)\right)\right]$ (in percent, b). The grey areas and crosses refer to the non-significant grid cells in the 95 % confidence limits (2σ level). Note that the scales are different.

Figure 3(a) Example of daily time series of IASI O3 measurements (dark blue) and of the fitted seasonal regression models with (light blue) and without (orange) the linear term in the troposphere. Daily time series (b) of the deseasonalized O3 (observations and regression models) and (c) of the difference of the fitted models with and without the linear trend term as well as the adjusted annual trend (pink and grey lines, respectively; c; given in DU). The RMSE (annual and for the JJA period in DU) and the trend values (annual and for the JJA period in DU yr−1) are also indicated.

The seasonal distributions of the contribution to the total variation in the MLT from the adjusted harmonic terms and explanatory variables, which account for the “natural” variability, and from the LT term are shown in Fig. 2a and b, respectively. The grey areas in the LT panels refer to the LT terms rejected by the stepwise backward elimination process. The crosses indicate that the trend estimate in the grid cell is non-significant at the 95 % confidence limits (2σ level) when accounting for the autocorrelation in the noise residual at the end of the elimination procedure. While the large influences of the seasonal variations and of the main drivers – namely ENSO, NAO and QBO – on the IASI O3 records have been clearly attested in Wespes et al. (2017), we demonstrate in Fig. 2 that the LT also contributes considerably to the O3 variations detected by IASI in the troposphere. The LT contribution generally ranges from 15 to 50 %, with the largest values ( 30 to  50 %) being observed at mid- and high latitudes in the SH (30–70 S) and in the NH (45–70 N) in summer. In the SH, they are associated with the smallest tropospheric O3 columns (Fig. 1a) and the smallest natural contributions (< 25 %; Fig. 2a), while in the NH summer, they interestingly correspond to large MLT O3 columns, large natural contributions ( 50 to  60 %) and the smallest RMSE (< 12 % or < 3 DU). From the annual regression model, the natural variation and the trend respectively contribute 30–85 % and up to 40 % to the total variation in the MLT.

Figure 4Seasonal distribution of the differences between the RMSE of the regression fits with and without the linear trend term [(RMSE_w/o_LT RMSE_with_LT)  RMSE_with_LT×100] (in percent). The blue circle in the JJA panel refers to the case presented in Fig. 3.

In Fig. 3, we further investigate the robustness of the estimated trends by performing sensitivity tests in regions of significant trend contributions (e.g., in the NH midlatitudes in summer; see Fig. 2). The  9-year time series of IASI O3 daily averages (dark blue) along with the results from the seasonal regression model with and without the LT term included in the model (light blue and orange lines, respectively) are represented in Fig. 3a for one specific location (highlighted by a blue circle in the JJA panel in Fig. 4). Figure 3b provides the deseasonalized IASI (dark blue line) and fitted time series (calculated by subtracting the adjusted seasonal cycle from the time series) resulting from the adjustment with and without the LT term included in the MLR model (light blue and orange lines, respectively). The differences between the fitted models with and without LT are shown in Fig. 3c (pink lines). They match fairly well the adjusted trend over the IASI period (Fig. 3c, grey lines; the trend and the RMSE values are also indicated) and the adjustment without LT leads to larger residuals (e.g., RMSE_JJA_w/o_LT=3.37 DU vs. RMSE_JJA_with_LT=3.21 DU in summer). This result demonstrates the possibility of capturing trend information from  9 years of IASI-MLR with only some compensation effects from the other explanatory variables, contrary to what was observed when considering a shorter period of measurements or less-frequent sampling (i.e., monthly dataset; e.g., Wespes et al., 2016). It is also worth mentioning that the O3 change calculated over the whole IASI dataset in summer is larger than the RMSE of the model residuals (increase of 6.50 ± 2.35 DU vs. RMSE of 3.21 DU), underlying the statistical relevance of trend estimates.

Figure 5(a) Annual and (b) seasonal distributions of the adjusted trends in DU yr−1 from the multi-linear regression models. The grey areas and crosses refer to the non-significant grid cells in the 95 % confidence limits (2σ level).

The robustness of the adjusted trend is verified at the global scale in Fig. 4 which represents the seasonal distributions of the relative differences in the RMSE with and without the LT included in the MLR model, calculated as (RMSE_w/o_LT RMSE_with_LT)  RMSE_with_LT×100 (in percent). An increase in the RMSE when excluding LT from the MLR is observed almost everywhere in regions of significant trend contributions (Fig. 2), especially in mid- and high latitudes of the SH and of the NH in summer, where it reaches 10 %. This result indicates that adjusting LT improves the performance of the model and, hence, that a trend signal is well captured by IASI at a regional scale in the troposphere. From the annual model, the increase in the RMSE only reaches 5 % at mid- and high latitudes of the SH (data not shown). In regions of weak or non-significant trend contribution (see crosses in Fig. 2), no improvement is logically found.

4 O3 trend over 2008–2017

## 4.1 Annual and seasonal trends

The annual and the seasonal distributions of the fitted LT terms which are retained in the annual and the seasonal MLR models by the stepwise elimination procedure are respectively represented in Fig. 5a and b (in DU yr−1). Generally, the mid- and high latitudes of both hemispheres and, more particularly, the NH midlatitudes in summer reveal significant negative trends, while the tropics are mainly characterized by non-significant or weak significant trends. Even if trends in emissions have already been able to qualitatively explain measured tropospheric O3 trends over specific regions, the magnitude and the trend estimates considerably vary between independent measurement datasets (e.g., Cooper et al., 2014; the TOAR-Climate report coordinated by the International Global Atmospheric Chemistry Project and available soon at http://www.igacproject.org/activities/TOAR, last access: 8 May 2018, of Gaudel et al., 2018, and references therein) for the reasons discussed in Sect. 1, and they are not reproduced/explained by model simulations (e.g., Jonson et al., 2006; Cooper et al., 2010; Logan et al., 2012; Wilson et al., 2012; Hess and Zbinden, 2013; and references therein). As a result, comparing/reconciling the adjusted trends with independent measurements, even on a qualitative basis, remains difficult. Nevertheless, several of the statistically significant features observed in Fig. 5 show, interestingly, qualitative consistency with respect to recent published findings:

• The SH tropical region extending from the Amazon to tropical eastern Indian Ocean seems to indicate a general increase with, for example, a DJF trend of  0.23 ± 0.18 DU yr−1 (i.e., 2.09 ± 1.70 DU over the IASI measurement period), despite the large IAV in the MLT which characterizes the tropics and which likely explains the high frequency of non-significant trends. Enhanced O3 levels over that region have already been analyzed for previous periods (e.g., Logan, 1985; Logan and Kirchhoff, 1986; Fishman et al., 1991; Moxim and Levy, 2000; Thompson and Wallace, 2000; Thompson et al., 2007; Sauvage et al., 2006, 2007; Heue et al., 2016; Ebojie et al., 2016; Archibald et al., 2018; Leventidou et al., 2017). For instance, the large O3 enhancement of  0.36 ± 0.25 DU yr−1 (i.e., 3.3 ± 2.3 DU over the whole IASI period) stretching from southern Africa to Australia over the northeast of Madagascar during the austral winter–spring likely originates from large IAV in the subtropical jet-related stratosphere–troposphere exchanges, which have been found to primarily contribute to the tropospheric O3 trends over that region (Liu et al., 2016, 2017). Nevertheless, this finding should be mitigated by the fact that the trend value in the SH tropics is of the same magnitude as the RMSE of the regression residuals ( 2 to  4.5 DU; see Fig. 1).

Conversely, the tropical Pacific region exhibits significant negative trends that are similar to those reported from UV sounders in Ebojie et al. (2016) and in Leventidou et al. (2017) over previous periods, while Heue et al. (2016) mainly report a significant positive trend over that region.

• The trends over Southeast Asia are mostly non-significant and vary by season. In spring–summer, some grid cells in India, in mainland China and eastern downwind China exhibit significant positive trends reaching  0.45 DU yr−1 (i.e.,  4.2 DU over the IASI measurement period). This tends to indicate that the tropospheric O3 increases, which have been shown to mainly result from a strong positive trend in the Asian emissions over the past decades (e.g., Zhao et al., 2013; Cooper et al., 2014; Zhang et al., 2016; Cohen et al., 2017; Tarasick et al., 2018; and references therein) but also from a substantial change in the stratospheric contribution (Verstraeten et al., 2015), persist through 2008–2017 despite the recent decrease in O3 precursor emissions recorded in China after 2011 (e.g., Duncan et al., 2016; Krotkov et al., 2016; Miyazaki et al., 2017; van der A et al., 2017). This would indicate that this decrease is probably too recent/weak to recover the 2008 O3 levels over the entire region. Significant positive trends over Southeast Asia have also been reported from UV sounders over previous periods (e.g., Ebojie et al., 2016). Note, however, that this finding has to be taken carefully given the large model residuals (RMSE of  2 to  4 DU; see Sect. 3, Fig. 1) over that region. Finally, the large uncertainty in trend estimates over Southeast Asia might reflect the large IAV in the biomass-burning emissions and lightning NOx sources, in addition to the recent changes in emissions.

• The mid- and high latitudes of the SH show clear patterns of negative trends, throughout the year and in a more pronounced manner during winter–spring, with larger amplitudes than those of the RMSE values (0.33 ± 0.14 DU yr−1 on average in the 35–65 S band; i.e., a trend amplitude of $|\mathrm{3.1}|$± 1.3 DU over the studied period vs. an RMSE value of  2.5 DU). These significant negative trends in the SH are hard to explain, but considering the stratospheric contribution into the tropospheric columns (natural and artificial due to the limited IASI vertical sensitivity) in the mid- and high latitudes of the SH ( 40 to  60 %; see Supplement in Wespes et al., 2016) and the negative significant trends previously reported from IASI in the UTLS/low stratosphere in the 30–50 S band, they could be in line with those derived by Zeng et al. (2017) in the UTLS for a clean rural site of the SH (Lauder, New Zealand), which mainly originate from increasing tropopause height and O3-depleting substances. Significant negative change in tropospheric O3 over these regions were also reported in Ebojie et al. (2016).

• In the NH, a band-like pattern of negative trends is observed in the 40–75 N band covering Europe and North America, especially during summer. Averaged annual trend of 0.31 ± 0.17 DU yr−1 and summer trend of 0.47 ± 0.22 DU yr−1 (i.e., 2.87 ± 1.57 and 4.36 ± 2.02 DU, respectively, from January 2008 to May 2017) are estimated in that latitudinal band. These trend values are significantly larger than the RMSEs of the MLR model (< 3.5 DU in JJA; see Sect. 3, Fig. 1). Interestingly, both the annual and summer trends are amplified relative to the ones calculated in the midlatitudes of the NH over the 2008–2013 period of IASI measurements (0.19 ± 0.05 DU yr−1 and 0.30 ± 0.10 DU yr−1 for the annual and the summer trends, respectively, calculated in the 30–50 N band; see Wespes et al., 2016). This finding is line with previous studies which point out a possible leveling off of tropospheric O3 in summer due to the decline of anthropogenic O3 precursor emissions observed since 2010-2011 in North America, in western Europe and also in some regions of China (e.g., Cooper et al., 2010, 2012; Logan et al., 2012; Parrish et al., 2012; Oltmans et al., 2013; Simon et al., 2015; Ebojie et al., 2016; Archibald et al., 2018; Miyazaki et al., 2017). The present study even goes a step further by suggesting a possible decrease in the tropospheric O3 levels. Archibald et al. (2018) recently reported a net decrease of  5 % in the global anthropogenic NOx emissions in the 30–90 N latitude band, which is consistent with the annual significant negative trend of 0.27 ± 0.15 DU yr−1 for O3 estimated from IASI in that band. We should also note that, even if these latitudes are characterized by the lowest stratospheric contribution ( 30 to  45 %; see Supplement in Wespes et al., 2016), it might partly mask/attenuate the variability in the tropospheric O3 levels.

## 4.2 Expected year for trend detection

In this section, we further verify that it is indeed possible to infer, from the studied IASI period, the significant negative trend derived in the 40–75 N band in summer ($|\mathrm{0.5}|$ DU yr−1 on average; see Sect. 4.1) by determining the expected year from which such a trend amplitude would be detectable at a global scale. This is achieved by estimating the minimum duration (with probability 0.90) of the IASI O3 measurements that would be required to detect a trend of a specified magnitude, and its 95 % confidence level, following the formalism developed in Tiao et al. (1990) and in Weatherhead et al. (1998):

$\begin{array}{}\text{(1)}& & {N}^{\ast }\approx {\left[\frac{\mathrm{3.3}\cdot {\mathit{\sigma }}_{\mathit{\epsilon }}}{\left|{\mathit{\tau }}_{\mathrm{yr}}\right|}\cdot \sqrt{\frac{\mathrm{1}+\mathit{\varphi }}{\mathrm{1}-\mathit{\varphi }}}\right]}^{\frac{\mathrm{2}}{\mathrm{3}}},\text{(2)}& & {\mathrm{CL}}_{{N}^{\ast }}=\left[{N}^{\ast }\cdot {e}^{-B};{N}^{\ast }\cdot {e}^{+B},\right]\end{array}$

where Nis the number of the required years, σε is the standard deviation of the autoregressive noise residual εt, τyr is the magnitude of the trend per year and ϕ is the lag-1 autocorrelation of the noise. The magnitudes of the variation and the autocorrelation in the noise residuals are taken into account for reaching a better trend estimate precision. Given that large variance (${\mathit{\sigma }}_{\mathit{\epsilon }}^{\mathrm{2}}\right)$ and large positive autocorrelation ϕ of the noise induce small signal-to-noise ratio and long trend-like segments in the dataset, respectively, these two parameters increase the number of years that would be required for detecting a specified trend. CL${}_{{N}^{\ast }}$ is the 95 % confidence limit, which is not symmetric around Nand depends on B, an estimated uncertainty factor calculated as $\frac{\mathrm{4}}{\mathrm{3}\sqrt{D}}\sqrt{\frac{\mathrm{1}+\mathit{\varphi }}{\mathrm{1}-\mathit{\varphi }}}$, with D the number of days in the IASI datasets, which accounts for the uncertainty in ϕ (the uncertainty in σε being negligible given that only a few years of data are needed to estimate it; see Weatherhead et al., 1998). As a result, based on the available IASI-A and proxies datasets and assuming that the MLR model used in this study is accurate, we estimate, in Fig. 6a and b, the expected year when an O3 trend amplitude of $|\mathrm{5}|$ DU decade−1 (i.e., ${\mathit{\tau }}_{\mathrm{yr}}=|\mathrm{0.5}|$ DU yr−1 which corresponds to the averaged absolute value of the fitted negative trends in the NH summer; see Fig. 5b) is detectable, and its associated maximal confidence limit, respectively. The results in Fig. 6 clearly demonstrate the possibility of inferring, from the available IASI dataset, such significant trends in the mid- and high latitudes of the NH in summer and fall (trend detectable from  2016 to 2017 with an uncertainty of  6 to  9 months; see Fig. 6b). Conversely, the tropical regions and the NH in winter–spring would require additional  6 months to  6 years of measurements to detect an amplitude of $|\mathrm{0.5}|$ DU yr−1 (trend significant only from  2017 to 2023 or after depending on the location and the season). Note also that the strongest negative trends (up to 0.85 DU yr−1, i.e., ${\mathit{\tau }}_{\mathrm{yr}}=|\mathrm{0.85}|$ DU yr−1; see Fig. 5b) observed in specific regions of the NH midlatitudes would only require  6 years of IASI measurements for being detected. The mid- and high latitudes of the SH would require the shortest period of IASI measurement for detecting a specified trend, with only  7 years ± 6 months of IASI measurements to detect a $|\mathrm{0.5}|$ DU yr−1 trend (trend detectable from  2015). That band-like pattern in the SH corresponds to the region with the weakest IAV and contribution from large-scale dynamical modes of variability in the IASI MLT columns (see Sect. 3, Figs. 1 and 2), which translates into small ${\mathit{\sigma }}_{\mathit{\epsilon }}^{\mathrm{2}}$ and ϕ. Note however that an additional few months of IASI data are required to confirm the smaller negative trend of 0.35 DU yr−1 measured on average in the SH (see Fig. 5; a period  9 years ± 6 months being necessary to detect a trend amplitude of $|\mathrm{3.5}|$ DU decade−1). Given that large σε means large noise residual in the IASI data, the regions of short (long) required measurement period coincide, as expected, well with the small (high) RMSE values of the regression residuals (see Sect. 3, Fig. 1).

Figure 6(a) Estimated year of tropospheric IASI O3 trend detection (with a probability of 90 %) for a given trend of $|\mathrm{5}|$ DU decade−1 starting at the beginning of the studied period (20080101) and (b) associated maximal confidence limits from the annual (top panel) and the seasonal (bottom panels) regression models.

The regions of the longest measurement periods required in the tropics for the $|\mathrm{5}|$ DU decade−1 trend detection (up to  16 years of IASI data) correspond to known patterns of widespread high O3: (a) above intense biomass burning in Amazonia and eastwards across the tropical Atlantic (Logan and Kirchhoff, 1986; Fishman et al., 1991; Moxim and Levy, 2000; Thompson and Wallace, 2000; Thompson et al., 2007; Sauvage et al., 2007), (b) eastwards of Africa across the southern Indian Ocean which is subject to large variations in the stratospheric influences during the winter–spring austral period (JJA–SON; Liu et al., 2016, 2017), (c) eastwards of Africa across the northern Indian Ocean to India likely due to large lightning NOx emissions above central Africa during the wet season associated with the northeastward jet conducting a so-called “O3 river” (Tocquer et al., 2015) and (d) above regions of positive ENSO “chemical” effect in equatorial Asia/Australia and eastwards above northern and southern tropical regions (Wespes et al., 2016) explained by reduced rainfall and biomass fires during El Niño conditions (e.g., Worden et al., 2013). In fact, most of these patterns (a, b and d) are closely connected with strong El Niño events which extend the duration of the dry season and cause severe droughts, producing intense biomass-burning emissions, for instance, over South America (e.g., Chen et al., 2011; Lewis et al., 2011) and South Asia/Australia (e.g., Oman et al., 2013; Valks et al., 2014; Ziemke et al., 2015), and which alter the tropospheric circulation by increasing the transport of stratospheric O3 into the troposphere (e.g., Voulgarakis et al., 2011b; Neu et al., 2014) and the transport of biomass-burning air masses to the Indian Ocean (Zhang et al., 2012). In summary, these large-scale indirect ENSO-related variations in tropospheric O3 and the lightning NOx impact on O3, which are not accounted for in the MLR by specific representative proxies, are misrepresented in the regression models. They induce large noise residuals, i.e., large σε, and hence extend the time period needed to detect a trend of any given magnitude.

Figure 6, finally, suggests that a long duration is also required, especially in summer, above and east of China to quantify the anthropogenic impact on the local changes in the MLT: additional 3 ± 1.5 or 5 ± 1.5 years for a given $|\mathrm{5}|$ or $|\mathrm{3.5}|$ DU decade−1 trend are respectively calculated. This result could be explained by large perturbations in the MLT columns induced by recent decreases after decades of almost constant increases in surface emissions in China (e.g., Cohen et al., 2017; Miyazaki et al., 2017).

## 4.3 Multi-linear vs. single linear model

Figure 7Seasonal distributions of the fitted linear term trends (given in DU yr−1) derived from a single linear regression model. The crosses refer to the non-significant grid cells in the 95 % confidence limits (2σ level).

Even if MLR models have already been used for extracting trends in stratospheric and total O3 columns (e.g., Mäder et al., 2007; Frossard et al., 2013; Rieder et al., 2013; Knibbe et al., 2014), single linear regression (SLR) models, which do not differentiate the natural (chemical and dynamical) factors describing the O3 variability, are still commonly used (e.g., Cooper et al., 2014; the TOAR-Climate report of Gaudel et al., 2018; and references therein). They are, however, suspected to contribute to trend biases retrieved between independent measurements. In addition to the time-varying instrumental biases, trend biases can also be related to differences in the spatial and the temporal samplings (e.g., Doughty et al., 2011; Saunois et al., 2012; Lin et al., 2015), in the measurement period, in the upper boundary of the O3 columns, in the algorithm and in the vertical sensitivity of the measurements. The latter artificially alters the characteristics of the sounded layer by contaminations from the above and below layers leading to a mixing of the trend and also of the natural characteristics originating from these different layers (e.g., troposphere and stratosphere). The differences in the studied period, in the tropopause definition and in the spatio-temporal sampling might also imply differences in the natural influence on the measured O3 variations. While the impact of the natural contribution is taken into account in the MLR model, it might introduce an additional bias in the trend determined from SLR, making further challenging to compare trends estimated from a series of inhomogeneous independent measurements.

Substantial efforts in homogenizing independent tropospheric O3 column (TOC) datasets have been made in the TOAR-Climate assessment report (Gaudel et al., 2018), but large SLR trend biases remain between the TOAR datasets, in particular, between the satellite datasets where the lack of homogeneity in terms of tropopause calculation (same tropopause definition but different temperature profiles are used), of instrument vertical sensitivities and of spatial sampling have been specifically raised as possible causes for the trend divergence.

Figure 8Daily time series of (a) O3 measured by IASI (dark blue) and adjusted by the multivariate annual regression model (light blue), (b) the deseasonalized O3 and (c) the fitted signal of ENSO proxy (one of the main retained proxies in the multivariate regression model) calculated as xjXnorm,j (given in DU) along with the adjusted trends derived from the single and the multivariate linear regressions (SLR and MLR) over the equatorial central Pacific (region of negative ENSO “dynamical” effect). The RMSE of the MLR fit and the fitted SLR and MLR trend values are also indicated.

If reconciling the trend biases between the datasets by applying the vertical sensitivity of each measurement type to a common platform, as proposed in the TOAR-Climate assessment report is beyond the scope of this study and if, at this stage, there is no consensus in determining tropospheric O3 trends, the improvement in using a MLR instead of a SLR model for determining more accurate/realistic trends is explored here by comparing the seasonal distributions of the trends estimated from MLR (see Fig. 5b in Sect. 4.1.) and from SLR (presented in Fig. 7). Note that the constant term in the SLR model is split into two components (covering the periods before and after the September 2010 jump) to take account of the observed jump (see Sect. 2). The highest differences in the fitted trends derived from the two methods are found in the tropics and in some regions of the midlatitudes of the NH. They likely result from overlaps between the LT term and other covariates. For instance, the regions with high significant SLR trends ( 0.3 to  0.6 DU yr−1 over the tropical western and middle Pacific) during the period extending from September to May match the regions with strong El Niño/Southern Oscillation influence (see Figs. 8 and 12 in Wespes et al., 2016). Conversely, the MLR model shows generally weak significant negative or non-significant trends in the Pacific El Niño region during that period and it would even need additional  3 to  4 years of IASI measurements for detecting the fitted SLR trends (see Sect. 4.2 above). The effect of ENSO in overestimating the fitted SLR trend is further illustrated in Fig. 8 which represents the time series of O3 observed by IASI and adjusted by the annual MLR model (Fig. 8a) along with the deseasonalized times series (Fig. 8b) and the fitted SLR and MLR trends (Fig. 8c). The fitted signal of the ENSO proxy from the MLR model (calculated as xjXnorm,j following Wespes et al., 2017) is also represented (Fig. 8c). That example clearly shows that the ENSO influence is considerably compensated for by the adjustment of the linear trend in the SLR (annual trend of 0.48 ± 0.06 DU yr−1 from SLR vs. 0.23 ± 0.16 DU yr−1 from MLR for that example). Finally, differences between the SLR and the MLR models are also observed in the region with strong positive NAO influence over the Icelandic/Arctic region during MAM (see Wespes et al., 2016, for a description of the NAO-related O3 changes). Conversely, the sub-tropical SH exhibits similar seasonal patterns of negative trends from both the SLR and the MLR. This results from the weak natural IAV and contributions in tropospheric O3 above that region (see Sect. 3, Figs. 1 and 2), which, hence, limits the compensation effects.

In summary, despite the fact that considering a long period of measurements is usually recommended in SLR studies to overcome the dynamical cycles and, hence, to help in differentiating their influences from trends, we show that, considering that some dynamics have irregular or no particular periodicity (e.g., NAO, ENSO), it is not accurate enough. Furthermore, accurate satellite measurements of tropospheric O3 at a global scale are quite recent, limiting the period of available and comparable datasets (e.g., Payne et al., 2017). As a consequence, we support here that using a reliable multivariate regression model based on geophysical parameters and adapted for each specific sounded layer is a robust method for differentiating a “true” trend from any other sources of variability and, hence, that it should help in resolving trend differences between independent datasets.

Figure 9Global distributions of (a) the correlation coefficients (${R}_{{\mathrm{O}}_{\mathrm{3}}-\mathrm{CO}}\right)$, (b) the regression slope (dO3 dCO in mol cm−2/mol cm−2) and (c) the covariances (COV${}_{{\mathrm{O}}_{\mathrm{3}}-\mathrm{CO}}$ in 1033 mol cm−2× mol cm−2) of daily median IASI tropospheric O3 and CO over January 2008–May 2017. Data are averaged over a 2.5× 2.5 grid box. Crosses in ${R}_{{\mathrm{O}}_{\mathrm{3}}-\mathrm{CO}}$ (a) refer to the non-significant grid cells in the 95 % confidence intervals (2σ level).

## 4.4 Continental influence

In this section, we use the capabilities of IASI to simultaneously measure O3 and CO in order (1) to differentiate tropospheric and stratospheric air masses, (2) to identify the regions influenced by the continental export/intercontinental transport of O3 pollution and (3) to evaluate that continental influence on tropospheric O3 trends as observed by IASI. Similar tracer correlations between CO and O3 have already been used to give insights into the photochemical O3 enhancement in air pollution transport (e.g., Parrish et al., 1993; Bertschi and Jaffe, 2005). However, there are only a few studies using global satellite data for this purpose (Zhang et al., 2006; Voulgarakis et al., 2011a; Kim et al., 2013) and the analysis of the O3–CO relationship for better understanding the origin of O3 trends in the troposphere has, to the best of our knowledge, never been explored.

Figure 9a and b show the seasonal patterns of the O3–CO correlations (referred as ${R}_{{\mathrm{O}}_{\mathrm{3}}-\mathrm{CO}}\right)$ and of the dO3 dCO regression slopes calculated in the troposphere (from the surface to 300 hPa) over the studied IASI period (January 2008–May 2017). The dO3 dCO regression slopes, which represent the so-called O3–CO enhancement ratio, are used to evaluate the photochemical O3 production in continental outflow regions. The ${R}_{{\mathrm{O}}_{\mathrm{3}}-\mathrm{CO}}$ and the dO3 dCO distributions are similar and clearly show regional and seasonal differences in the strength of the O3–CO relationships. The patterns of positive and negative correlations allows differentiating the outflow regions characterized by photochemical O3 production from precursors (including CO) or CO destruction (both identified by positive ${R}_{{\mathrm{O}}_{\mathrm{3}}-\mathrm{CO}}\right)$ from the regions characterized by O3 loss (chemical destruction or surface deposition) or by strong stratospheric contaminations (both identified by negative ${R}_{{\mathrm{O}}_{\mathrm{3}}-\mathrm{CO}}\right)$. Negative ${R}_{{\mathrm{O}}_{\mathrm{3}}-\mathrm{CO}}$ and dO3 dCO are measured at high latitudes of both hemispheres all over the year, but more specifically at high latitudes of the SH in summer–fall (with ${R}_{{\mathrm{O}}_{\mathrm{3}}-\mathrm{CO}}$ < 0.25 on averages in DJF and MAM). Given that high latitudes experience more O3 destruction than the low latitudes due to a lack of sunlight, the negative correlations for the high latitude regions might also reflect air masses originating from/characterizing the stratosphere due to natural intrusion or to artificial mixing with the troposphere introduced by the limited vertical sensitivity of IASI in the highest latitudes (stratospheric contribution varying between  40 and  65 %; see Supplement in Wespes et al., 2016). These processes are likely at the origin of the band-like pattern of negative trends in the SH discussed in Sects. 3 and 4.1. Negative ${R}_{{\mathrm{O}}_{\mathrm{3}}-\mathrm{CO}}$ and dO3 dCO are also found above the Caribbean, the Arabic Peninsula and the northern Indian Ocean in JJA/SON and the South Atlantic in DJF. They are in line with Kim et al. (2013) and they likely reflect the influence of lightning NOx which produce O3 but also OH oxidizing CO (e.g., Sauvage et al., 2007; Labrador et al., 2004).

Strong positive correlations are identified in both ${R}_{{\mathrm{O}}_{\mathrm{3}}-\mathrm{CO}}$ and dO3 dCO patterns over the tropical regions and for midlatitudes of both hemispheres during the peak of photochemistry in summer. Maxima (${R}_{{\mathrm{O}}_{\mathrm{3}}-\mathrm{CO}}$ > 0.8 and dO3 dCO > 0.5) are detected in continental pollution outflow regions in the NH, especially downwind Southeast Asia and over the southern Africa/Amazonia/South Atlantic region. These O3–CO correlation patterns from IASI are fully consistent with those measured by TES (Zhang et al., 2006, 2008; Voulgarakis et al., 2011a) and by OMI/AIRS (Kim et al., 2013), which have been interpreted with global chemical transport models (CTMs) as originating from Asian pollution influences and combustion sources including biomass burning, respectively. The high positive ${R}_{{\mathrm{O}}_{\mathrm{3}}-\mathrm{CO}}$ found in JJA at midlatitudes of the NH are detected to a lesser extent in DJF reflecting the decreasing photochemistry. It is also worth pointing out the clear hemispheric differences in Fig. 9 in the ${R}_{{\mathrm{O}}_{\mathrm{3}}-\mathrm{CO}}$ and dO3 dCO values at mid- and high latitudes and, in particular, the shift of positive ${R}_{{\mathrm{O}}_{\mathrm{3}}-\mathrm{CO}}$ and dO3 dCO towards higher latitudes of the NH during summer (e.g., ${R}_{{\mathrm{O}}_{\mathrm{3}}-\mathrm{CO}}=\mathrm{0.24}$ in summer vs. 0.038 in spring in the 35–55 N band). As a consequence, these results suggest that the band-like pattern of negative trends measured by IASI in summer might substantially reflect the continental pollution influence and, hence, that it could result from the decline of anthropogenic O3 precursor emissions. Nevertheless, interpreting O3–CO correlations in the free troposphere, especially in photochemically aged pollution plumes far from the emission sources towards the highest latitudes, remains complicated by the mixing of the continental combustion outflow with stratospheric air masses, in addition to the background dynamic and photochemistry (e.g., Liang et al., 2007).

Finally, we also provide in Fig. 9c the seasonal patterns of O3–CO covariances (COV${}_{{\mathrm{O}}_{\mathrm{3}}-\mathrm{CO}}\right)$. They confirm the band-like pattern of the weak natural variation captured in the SH midlatitudes (see Sects. 3 and 4.1) and help identifying the region downwind eastern China, the northern midlatitudes outflow region and the southern tropical region as the ones with the highest pollution variability, in addition to the strongest O3–CO correlations. To conclude, the particularly strong positive O3–CO relationship in terms of ${R}_{{\mathrm{O}}_{\mathrm{3}}-\mathrm{CO}}$, dO3 dCO and COV${}_{{\mathrm{O}}_{\mathrm{3}}-\mathrm{CO}}$ measured over and downwind northeast India/eastern China in summer in comparison with the ones measured downwind eastern USA and over Europe indicate that Southeast Asia experiences the most of the intense pollution episodes of the NH with the largest O3–CO variability (COV${}_{{\mathrm{O}}_{\mathrm{3}}-\mathrm{CO}}$ > 40 × 1033 mol2 cm−4) and the largest O3 enhancement (dO3 dCO > 0.5) over the last decade. The strong O3–CO relationship in that region is associated with the significant increase that is detected in the IASI O3 levels downwind east of Asia (see Sect. 4.1) despite the net decrease in O3 precursor emissions recorded in China after 2011 (e.g., Cohen et al., 2017; Miyazaki et al., 2017).

5 Conclusions

In this study, we have explored, for the first time, the possibility of inferring significant trends in tropospheric O3 from the first  10 years (January 2008–May 2017) of IASI daily measurements at a global scale. To this end, MLR analyses have been performed by applying a multivariate regression model (annual and seasonal formulations), including a linear trend term in addition to chemical and dynamical proxies, on gridded mean tropospheric ozone time series. This work follows on the analysis of the main dynamical drivers of O3 variations measured by IASI, which was recently published in Wespes et al. (2017). We have first verified the performance of the MLR models in explaining the variations in daily time series over the whole studied period. In particular, we have shown that the model reproduces a large part of the O3 variations (> 70 %) with small residuals errors (RMSE of  10 %) at northern latitudes in summer. We have then performed O3 trend sensitivity tests to verify the possibility of capturing trend characteristics independently from natural variations. Despite the weak contribution of trends to the total variation in the MLT O3 columns at a global scale, the results demonstrate the possibility of differentiating significant trends from the explanatory variables, especially in summer at mid- and high latitudes of the NH ( 40–75 N) where the contribution and the sensitivity of trends are the largest (contribution of  30 to  50 % and a  10 % increase in the RMSE excluding the LT in the model). We then focused on the interpretation of the global trend estimates. We have found interesting significant positive trends in the SH tropical region extending from the Amazon to the tropical eastern Indian Ocean and over Southeast Asia in spring–summer, which should however be carefully considered given the high RMSEs of the regression residuals in these regions. The MLR analysis reveals a band-like pattern of high significant negative trends in the NH mid- and high latitudes in summer (0.47 ± 0.22 DU yr−1 on average in the 40–75 N band). The statistical significance of such trend estimates is further verified by estimating, based on the autocorrelation and on the variance of the noise residuals, the minimum number of years of IASI measurements that are required to detect a trend of a | 5| DU decade−1 magnitude. The results clearly demonstrate the possibility of determining such a trend amplitude from the available IASI dataset and the used MLR model at northern mid- and high latitudes in summer, while much larger measurement periods are necessary elsewhere. In particular, the regions with the longest required period, in the tropics, highlight a series of known processes that are closely related to the El Niño dynamic, which underlies the lack of associated parameterizations in the MLR model. The importance of using reliable MLR models in understanding large-scale O3 variations and in determining trends is further explored by comparing the trends inferred from MLR and from SLR, the latter still being commonly used by the international community. The comparison has clearly highlighted the advantage of MLR in attributing the trend-like segments in natural variations, such as ENSO, to the right processes and, hence, in avoiding misinterpretation of “apparent” trends in the measurement datasets. Nevertheless, it is worth noting that there could be a possible impact of the sampling (because of the cloud and quality filters applied) and of the jump in September 2010 that has been identified in the IASI dataset (see Sect. 2), in both MLR and SLR estimations. Finally, by exploiting the simultaneous and vertically resolved O3 and CO measurements from IASI, we have provided and used the O3–CO correlations in the troposphere to help determine the origins/characteristics of patterns of negative or positive trends. The distributions have allowed us to identify, in particular, strong positive O3–CO correlations, regression slopes and covariance in the NH midlatitudes and northward during summer, which suggest a continental pollution influence in the NH band-like pattern of high significant negative trends recorded by IASI and, hence, a direct effect of the policy measures taken to reduce emissions of O3 precursor species.

Overall, this study supports the importance of using (1) high density and long-term homogenized satellite records, such as those provided by IASI and (2) complex models with predictor functions that describe the O3 regressors dependencies for more accurate determination of trends in tropospheric O3 – as required by the scientific community, e.g., in the Intergovernmental Panel on Climate Change (IPCC, 2013) – and for further resolving trend biases between independent datasets (Payne et al., 2017; the TOAR-Climate assessment report of Gaudel et al., 2018). Currently, no consensus in terms of O3 trends in the troposphere is reached from the available measurements (UV or IR satellites, O3 sondes, aircrafts, ground-based measurements, etc.) for several reasons (time-varying instrumental biases, differences in the methodology used for calculating trends, in the measurement period, in the upper boundary of the O3 columns, in the retrieval algorithm, in the spatio-temporal sampling, in the vertical sensitivity of the instrument, etc.; see Sect. 4.3 and the TOAR-Climate report of Gaudel et al., 2018). However, determination, with IASI, of robust trends in tropospheric O3 at the global scale will be achievable in the near future by merging homogeneous O3 profiles from the three successive instruments onboard Metop-A (2006); -B (2012) and –C (2018) platforms and from the IASI next-generation instrument onboard the Metop Second Generation series of satellites (Clerbaux and Crevoisier, 2013; Crevoisier et al., 2014). A long record of tropospheric O3 measurements will be also provided by the Cross-track Infrared Sounder (CrIS) onboard the Joint Polar Satellite System series of satellites.

Data availability
Data availability.

The IASI O3 and CO data respectively processed with FORLI-O3 and FORLI-CO v0151001 can be downloaded from the Aeris portal (http://iasi.aeris-data.fr/O3/ and http://iasi.aeris-data.fr/CO/, respectively; Aeris, 2017a, b).

Competing interests
Competing interests.

The authors declare that they have no conflict of interest.

Special issue statement
Special issue statement.

Acknowledgements
Acknowledgements.

IASI was developed and built under the purview of the Centre National d'Etudes Spatiales (CNES, France). It is flown onboard the Metop satellites as part of the EUMETSAT Polar System. The IASI L1 data are received through the EUMETCast near-realtime data distribution service. We acknowledge the financial support from the ESA O3-CCI and the Copernicus O3-C3S projects. The AC SAF project is acknowledged for supporting the IASI FORLI-O3 and FORLI-CO implementation at EUMETSAT. The research in Belgium is also funded by the Belgian State Federal Office for Scientific, Technical and Cultural Affairs and the European Space Agency (ESA Prodex IASI Flow and AC SAF).

Edited by: Mark Weber
Reviewed by: two anonymous referees

References

Aeris: The IASI O3 products processed with FORLI-O3 v20151001, available at: http://iasi.aeris-data.fr/O3/ (last access: 14 May 2018), 2017a.

Aeris: The IASI CO products processed with FORLI-CO v20151001, available at: http://iasi.aeris-data.fr/CO/ (last access: 14 May 2018), 2017b.

Anton, M., Loyola, D., Clerbaux, C., Lopez, M., Vilaplana, J., Banon, M., Hadji-Lazaro, J., Valks, P., Hao, N., Zimmer, W., Coheur, P., Hurtmans, D., and Alados-Arboledas, L.: Validation of the Metop-A total ozone data from GOME-2 and IASI using reference ground-based measurements at the Iberian peninsula, Remote Sens. Environ., 115, 1380–1386, 2011.

Archibald, A., Elshorbany, Y., et al.: Tropospheric Ozone Assessment Report: Critical review of the present-day and near future tropospheric ozone budget, Elem. Sci. Anth., in review, 2018.

Bertschi, I. T. and Jaffe, D. A.: Long-range transport of ozone, carbon monoxide, and aerosols to the NE Pacific troposphere during the summer of 2003: Observations of smoke plumes from Asian boreal fires, J. Geophys. Res., 110, D05303, https://doi.org/10.1029/2004JD005135, 2005.

Boynard, A., Hurtmans, D., Koukouli, M. E., Goutail, F., Bureau, J., Safieddine, S., Lerot, C., Hadji-Lazaro, J., Wespes, C., Pommereau, J.-P., Pazmino, A., Zyrichidou, I., Balis, D., Barbe, A., Mikhailenko, S. N., Loyola, D., Valks, P., Van Roozendael, M., Coheur, P.-F., and Clerbaux, C.: Seven years of IASI ozone retrievals from FORLI: validation with independent total column and vertical profile measurements, Atmos. Meas. Tech., 9, 4327–4353, https://doi.org/10.5194/amt-9-4327-2016, 2016.

Boynard, A., Hurtmans, D., Garane, K., Goutail, F., Hadji-Lazaro, J., Koukouli, M. E., Wespes, C., Keppens, A., Pommereau, J.-P., Pazmino, A., Balis, D., Loyola, D., Valks, P., Coheur, P.-F., and Clerbaux, C.: Validation of the IASI FORLI/Eumetsat ozone products using satellite (GOME-2), ground-based (Brewer-Dobson, SAOZ) and ozonesonde measurements, Atmos. Meas. Tech. Discuss., https://doi.org/10.5194/amt-2017-461, in review, 2018.

Chen, Y., Randerson, J. T., Morton, D. C., DeFries, R. S., Collatz, G. J., Kasibhatla, P. S., Giglio, L., Jin, Y., and Marlier, M. E.: Forecasting Fire Season Severity in South America Using Sea Surface emperature Anomalies, Science, 334, 787–791, https://doi.org/10.1126/science.1209472, 2011.

Clarisse, L., R'Honi, Y., Coheur, P.-F., Hurtmans, D., and Clerbaux, C.: Thermal infrared nadir observations of 24 atmospheric gases, Geophys. Res. Lett., 38, L10802, https://doi.org/10.1029/2011GL047271, 2011.

Clerbaux, C. and Crevoisier, C.: New Directions: Infrared remote sensing of the troposphere from satellite: Less, but better, Atmos. Environ., 72, 24–26, 2013.

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.

Cohen, Y., Petetin, H., Thouret, V., Marécal, V., Josse, B., Clark, H., Sauvage, B., Fontaine, A., Athier, G., Blot, R., Boulanger, D., Cousin, J.-M., and Nédélec, P.: Climatology and long-term evolution of ozone and carbon monoxide in the UTLS at northern mid-latitudes, as seen by IAGOS from 1995 to 2013, Atmos. Chem. Phys. Discuss., https://doi.org/10.5194/acp-2017-778, in review, 2017.

Cooper, O., Parrish, D., Stohl, A., Trainer, M., Nédélec, P., Thouret, V., Cammas, J.-P., Oltmans, S., Johnson, B., and Tarasick, D.: Increasing springtime ozone mixing ratios in the free troposphere over western North America, Nature, 463, 344–348, https://doi.org/10.1038/nature08708, 2010.

Cooper, O. R., Gao, R.-S., Tarasick, D., Leblanc, T., and Sweeney, C.: Long-term ozone trends at rural ozone monitoring sites across the United States, 1990–2010, J. Geophys. Res., 117, D22307, https://doi.org/10.1029/2012JD018261, 2012.

Cooper, O. R., Parrish, D. D., Ziemke, J., Balashov, N. V., Cupeiro, M., Galbally, I. E., Gilge, S., Horowitz, L., Jensen, N. R., Lamarque, J.-F., Naik, V., Oltmans, S. J., Schwab, J., Shindell, D. T., Thompson, A. M., Thouret, V., Wang, Y., and Zbinden, R. M.: Global distribution and trends of tropospheric ozone: An observation-based review, Elem. Sci. Anth., 2, 000029, https://doi.org/10.12952/journal.elementa.000029, 2014.

Crevoisier, C., Clerbaux, C., Guidard, V., Phulpin, T., Armante, R., Barret, B., Camy-Peyret, C., Chaboureau, J.-P., Coheur, P.-F., Crépeau, L., Dufour, G., Labonnote, L., Lavanant, L., Hadji-Lazaro, J., Herbin, H., Jacquinet-Husson, N., Payan, S., Péquignot, E., Pierangelo, C., Sellitto, P., and Stubenrauch, C.: Towards IASI-New Generation (IASI-NG): impact of improved spectral resolution and radiometric noise on the retrieval of thermodynamic, chemistry and climate variables, Atmos. Meas. Tech., 7, 4367004385, https://doi.org/10.5194/amt-7-4367-2014, 2014.

Doughty, D. C., Thompson, A. M., Schoeberl, M. R., Stajner, I., Wargan, K., and Hui, W. C. J.: An intercomparison of tropospheric ozone retrievals derived from two Aura instruments and measurements in western North America in 2006, J. Geophys. Res., 116, D06303, https://doi.org/10.1029/2010JD014703, 2011.

Dufour, G., Eremenko, M., Griesfeller, A., Barret, B., LeFlochmoën, E., Clerbaux, C., Hadji-Lazaro, J., Coheur, P.-F., and Hurtmans, D.: Validation of three different scientific ozone products retrieved from IASI spectra using ozonesondes, Atmos. Meas. Tech., 5, 611–630, https://doi.org/10.5194/amt-5-611-2012, 2012.

Duncan, B. N., Lamsal, L. N., Thompson, A. M., Yoshida, Y., Lu, Z., Streets, D. G., Hurwitz, M. M., and Pickering, K. E.: A space-based, high-resolution view of notable changes in urban NOx pollution around the world (2004–2005), J. Geophys. Res., 121, 976–96, 2016.

Ebojie, F., Burrows, J. P., Gebhardt, C., Ladstätter-Weißenmayer, A., von Savigny, C., Rozanov, A., Weber, M., and Bovensmann, H.: Global tropospheric ozone variations from 2003 to 2011 as seen by SCIAMACHY, Atmos. Chem. Phys., 16, 417–436, https://doi.org/10.5194/acp-16-417-2016, 2016.

Fishman, J., Fakhruzzaman, K., Cros, B., and Nganda, D.: Identification of widespread pollution in the southern-hemisphere deduced from satellite analyses, Science, 252, 1693–1696, 1991.

Fishman, J., Creilson, J. K., Wozniak, A. E., and Crutzen, P. J.: Interannual variability of stratospheric and tropospheric ozone determined from satellite measurements, J. Geophys. Res., 110, D20306, https://doi.org/10.1029/2005JD005868, 2005.

Frossard, L., Rieder, H. E., Ribatet, M., Staehelin, J., Maeder, J. A., Di Rocco, S., Davison, A. C., and Peter, T.: On the relationship between total ozone and atmospheric dynamics and chemistry at mid-latitudes – Part 1: Statistical models and spatial fingerprints of atmospheric dynamics and chemistry, Atmos. Chem. Phys., 13, 147–164, https://doi.org/10.5194/acp-13-147-2013, 2013.

Fusco, A. C. and Logan, J. A.: Analysis of 1970–1995 trends in tropospheric ozone at northern hemisphere midlatitudes with the GEOSCHEM model, J. Geophys. Res., 108, 4449, https://doi.org/10.1029/2002JD002742, 2003.

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., Garcia, O., Granados-Muñoz, M. J., Hannigan, J., Hase, F., Hassler, B., Huang, G., Hurtmans, D., Jaffe, D., Jones, N., Kalabokas, P., Kerridge, B., Kulawik, S., Latter, B., Leblanc, T., Le Flochmoën, E., Lin, W., Liu, J., Liu, X., Mahieu, E., McClure-Begley, A., Neu, J., 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., Thouret, V., Thompson, A. M., Trickl, T., Weatherhead, E., Wespes, C., Worden, H., 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, in review, 2018.

Gazeaux, J., Clerbaux, C., George, M., Hadji-Lazaro, J., Kuttippurath, J., Coheur, P.-F., Hurtmans, D., Deshler, T., Kovilakam, M., Campbell, P., Guidard, V., Rabier, F., and Thépaut, J.-N.: Intercomparison of polar ozone profiles by IASI/MetOp sounder with 2010 Concordiasi ozonesonde observations, Atmos. Meas. Tech., 6, 613–620, https://doi.org/10.5194/amt-6-613-2013, 2013.

Hess, P. G. and Zbinden, R.: Stratospheric impact on tropospheric ozone variability and trends: 1990–2009, Atmos. Chem. Phys., 13, 649–674, https://doi.org/10.5194/acp-13-649-2013, 2013.

Heue, K.-P., Coldewey-Egbers, M., Delcloo, A., Lerot, C., Loyola, D., Valks, P., and van Roozendael, M.: Trends of tropical tropospheric ozone from 20 years of European satellite measurements and perspectives for the Sentinel-5 Precursor, Atmos. Meas. Tech., 9, 5037–5051, https://doi.org/10.5194/amt-9-5037-2016, 2016.

Hilton, F., Armante, R., August, T., Barnet, C., Bouchard, A., Camy-Peyret, C., Capelle, V., Clarisse, L., Clerbaux, C., Coheur, P.-F., Collard, A., Crevoisier, C., Dufour, G., Edwards, D., Faijan, F., Fourrié, N., Gambacorta, A., Goldberg, M., Guidard, V., Hurtmans, D., Illingworth, S., Jacquinet-Husson, N., Kerzenmacher, T., Klaes, D., Lavanant, L., Masiello, G., Matricardi, M., McNally, A., Newman, S., Pavelin, E., Payan, S., Péquignot, E., Peyridieu, S., Phulpin, T., Remedios, J., Schlüssel, P., Serio, C., Strow, L., Stubenrauch, C., Taylor, J., Tobin, D., Wolf, W., and Zhou, D: Hyperspectral Earth Observation from IASI: Five Years of Accomplishments, B. Am. Meteorol. Soc., 93, 347–370, 2012.

Hurtmans, D., Coheur, P., Wespes, C., Clarisse, L., Scharf, O., Clerbaux, C., Hadji-Lazaro, J., George, M., and Turquety, S.: FORLI radiative transfer and retrieval code for IASI, J. Quant. Spectrosc. Ra., 113, 1391–1408, 2012.

Intergovernmental Panel on Climate Change (IPCC): The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, edited by: Stocker, T. F., Qin, D., Plattner, G.-K., Tignor, M., Allen, S. K., Boschung, J., Nauels, A., Xia, Y., Bex, V., and Midgley, P. M., 164–270, Cambridge Univ. Press, Cambridge, UK and New York, USA, 2013.

Jonson, J. E., Simpson, D., Fagerli, H., and Solberg, S.: Can we explain the trends in European ozone levels?, Atmos. Chem. Phys., 6, 51–66, https://doi.org/10.5194/acp-6-51-2006, 2006.

Keppens, A., Lambert, J.-C., Granville, J., Hubert, D., Verhoelst, T., Compernolle, S., Latter, B., Kerridge, B., Siddans, R., Boynard, A., Hadji-Lazaro, J., Clerbaux, C., Wespes, C., Hurtmans, D. R., Coheur, P.-F., van Peet, J., van der A, R., Garane, K., Koukouli, M. E., Balis, D. S., Delcloo, A., Kivi, R., Stübi, R., Godin-Beekmann, S., Van Roozendael, M., and Zehner, C.: Quality assessment of the Ozone_cci Climate Research Data Package (release 2017): 2. Ground-based validation of nadir ozone profile data products, in preparation, 2018.

Kim, P. S., Jacob, D. J., Liu, X., Warner, J. X., Yang, K., Chance, K., Thouret, V., and Nedelec, P.: Global ozone–CO correlations from OMI and AIRS: constraints on tropospheric ozone sources, Atmos. Chem. Phys., 13, 9321–9335, https://doi.org/10.5194/acp-13-9321-2013, 2013.

Knibbe, J. S., van der A, R. J., and de Laat, A. T. J.: Spatial regression analysis on 32 years of total column ozone data, Atmos. Chem. Phys., 14, 8461–8482, https://doi.org/10.5194/acp-14-8461-2014, 2014.

Krotkov, N. A., McLinden, C. A., Li, C., Lamsal, L. N., Celarier, E. A., Marchenko, S. V., Swartz, W. H., Bucsela, E. J., Joiner, J., Duncan, B. N., Boersma, K. F., Veefkind, J. P., Levelt, P. F., Fioletov, V. E., Dickerson, R. R., He, H., Lu, Z., and Streets, D. G.: Aura OMI observations of regional SO2 and NO2 pollution changes from 2005 to 2015, Atmos. Chem. Phys., 16, 4605–4629, https://doi.org/10.5194/acp-16-4605-2016, 2016.

Labrador, L. J., von Kuhlmann, R., and Lawrence, M. G.: Strong sensitivity of the global mean OH concentration and the tropospheric oxidizing efficiency to the source of NOx from lightning, Geophys. Res. Lett., 31, L06102, https://doi.org/10.1029/2003GL019229, 2004.

Leventidou, E., Weber, M., Eichmann, K.-U., and Burrows, J. P.: Harmonisation and trends of 20-years tropical tropospheric ozone data, Atmos. Chem. Phys. Discuss., https://doi.org/10.5194/acp-2017-815, in review, 2017.

Lewis, S. L., Brando, P. M., Phillips, O. L., van der Heijden, G. M. F., and Nepstad, D.: The 2010 Amazon Drought, Science, 331, 554–554, https://doi.org/10.1126/science.1200807, 2011.

Liang, Q., Jaegle, L., Hudman, R. C., Turquety, S., Jacob, D. J., Avery, M. A., Browell, E. V., Sachse, G. W., Blake, D. R., Brune, W., Ren, X., Cohen, R. C., Dibb, J. E., Fried, A., Fuelberg, H., Porter, M., Heikes, B. G., Huey, G., Singh, H. B., and Wennberg, P. O.: Summertime influence of Asian pollution in the free troposphere over North America, J. Geophys. Res., 112, D12S11, https://doi.org/10.1029/2006JD007919, 2007.

Lin, M., Horowitz, L. W., Cooper, O. R., Tarasick, D., Conley, S., Iraci, L. T., Johnson, B., Leblanc, T., Petropavlovskikh, I., and Yates, E. L.: Revisiting the evidence of increasing springtime ozone mixing ratios in the free troposphere over western North America, Geophys. Res. Lett., 42, 8719–8728, https://doi.org/10.1002/2015GL065311, 2015.

Liu, J., Rodriguez, J. M., Thompson, A. M., Logan, J. A., Douglass, A. R., Olsen, M. A., Steenrod, S. D., and Posny, F.: Origins of tropospheric ozone interannual variation over Reunion: A model investigation, J. Geophys. Res.-Atmos., 121, 521–537, https://doi.org/10.1002/2015jd023981, 2016.

Liu, J., Rodriguez, J. M., Steenrod, S. D., Douglass, A. R., Logan, J. A., Olsen, M. A., Wargan, K., and Ziemke, J. R.: Causes of interannual variability over the southern hemispheric tropospheric ozone maximum, Atmos. Chem. Phys., 17, 3279–3299, https://doi.org/10.5194/acp-17-3279-2017, 2017.

Liu, X., Bhartia, P. K., Chance, K., Spurr, R. J. D., and Kurosu, T. P.: Ozone profile retrievals from the Ozone Monitoring Instrument, Atmos. Chem. Phys., 10, 2521–2537, https://doi.org/10.5194/acp-10-2521-2010, 2010.

Logan, J. A.: Tropospheric Ozone: Seasonal behaviour, Trends, and Anthropogenic Influence, J. Geophys. Res., 90, 10463–10482, 1985.

Logan, J. A. and Kirchhoff, V. W. J. H.: Seasonal variations of tropospheric ozone at Natal, Brazil, J. Geophys. Res., 91, 7875–7881, https://doi.org/10.1029/JD091iD07p07875, 1986.

Logan, J. A., Prather, M. J., Wofsy, S. C., and McElroy, M. B.: Tropospheric chemistry: A global perspective, J. Geophys. Res., 86, 7210–7254, https://doi.org/10.1029/JC086iC08p07210, 1981.

Logan, J. A., Staehelin, J., Megretskaia, I. A., Cammas, J.-P., Thouret, V., Claude, H., De Backer, H., Steinbacher, M., Scheel, H.-E., Stübi, R., Fröhlich, M., and Derwent, R.: Changes in ozone over Europe: Analysis of ozone measurements from sondes, regular aircraft (MOZAIC) and alpine surface sites, J. Geophys. Res., 117, D09301, https://doi.org/10.1029/2011JD016952, 2012.

Mäder, J. A., Staehelin, J., Brunner, D., Stahel, W. A., Wohltmann, I., and Peter, T.: Statistical modelling of total ozone: Selection of appropriate explanatory variables, J. Geophys. Res., 112, D11108, https://doi.org/10.1029/2006JD007694, 2007.

Miyazaki, K., Eskes, H., Sudo, K., Boersma, K. F., Bowman, K., and Kanaya, Y.: Decadal changes in global surface NOx emissions from multi-constituent satellite data assimilation, Atmos. Chem. Phys., 17, 807–837, https://doi.org/10.5194/acp-17-807-2017, 2017.

Moxim, W. J. and Levy II, H.: A model analysis of the tropical South Atlantic Ocean tropospheric ozone maximum: The interaction of transport and chemistry, J. Geophys. Res., 105, 17393–17415, https://doi.org/10.1029/2000JD900175, 2000.

Neu, J. L., Flury, T., Manney, G. L., Santee, M. L., Livesey, N. J., and Worden, J.: Tropospheric ozone variations governed by changes in stratospheric circulation, Nat. Geosci., 7, 340–344, https://doi.org/10.1038/ngeo2138, 2014.

Oetjen, H., Payne, V. H., Kulawik, S. S., Eldering, A., Worden, J., Edwards, D. P., Francis, G. L., Worden, H. M., Clerbaux, C., Hadji-Lazaro, J., and Hurtmans, D.: Extending the satellite data record of tropospheric ozone profiles from Aura-TES to MetOp-IASI: characterisation of optimal estimation retrievals, Atmos. Meas. Tech., 7, 4223–4236, https://doi.org/10.5194/amt-7-4223-2014, 2014.

Oltmans, S. J, Lefohn, A. S., Shadwick, D., Harris, J. M., Scheel, H.-E., Galbally, I., Tarasick, D. W., Johnson, B. J., Brunke, E. G., Claude, H., Zeng, G., Nichol, S., Schmidlin, F. J., Davies, J., Cuevas, E., Redondas, A., Naoe, H., Nakano, T., and Kawasato, T.: Recent Tropospheric Ozone Changes – A Pattern Dominated by Slow or No Growth, Atmos. Environ., 67, 331–351, https://doi.org/10.1016/j.atmosenv.2012.10.057, 2013.

Oman, L. D., Douglass, A. R., Ziemke, J. R., Rodriguez, J. M., Waugh, D. W., and Nielsen, J. E.: The ozone response to ENSO in Aura satellite measurements and a chemistry–climate simulation, J. Geophys. Res.-Atmos., 118, 965–976, 2013.

Parrington, M., Palmer, P. I., Henze, D. K., Tarasick, D. W., Hyer, E. J., Owen, R. C., Helmig, D., Clerbaux, C., Bowman, K. W., Deeter, M. N., Barratt, E. M., Coheur, P.-F., Hurtmans, D., Jiang, Z., George, M., and Worden, J. R.: The influence of boreal biomass burning emissions on the distribution of tropospheric ozone over North America and the North Atlantic during 2010, Atmos. Chem. Phys., 12, 2077–2098, https://doi.org/10.5194/acp-12-2077-2012, 2012.

Parrish, D. D., Holloway, J. S., Trainer, M., Murphy, P. C., Forbes, G. L., and Fehsenfeld, F. C.: Export of North American Ozone Pollution to the North Atlantic Ocean, Science, 259, 1436–1439, 1993.

Parrish, D. D., Law, K. S., Staehelin, J., Derwent, R., Cooper, O. R., Tanimoto, H., Volz-Thomas, A., Gilge, S., Scheel, H.-E., Steinbacher, M., and Chan, E.: Long-term changes in lower tropospheric baseline ozone concentrations at northern mid-latitudes, Atmos. Chem. Phys., 12, 11485–11504, https://doi.org/10.5194/acp-12-11485-2012, 2012.

Payne, V. H., Neu, J. L., and Worden, H. M.: Satellite observations for understanding the drivers of variability and trends in tropospheric ozone, J. Geophys. Res., 122, 6130–6134, https://doi.org/10.1002/2017JD026737, 2017.

Pommier, M., Clerbaux, C., Law, K. S., Ancellet, G., Bernath, P., Coheur, P.-F., Hadji-Lazaro, J., Hurtmans, D., Nédélec, P., Paris, J.-D., Ravetta, F., Ryerson, T. B., Schlager, H., and Weinheimer, A. J.: Analysis of IASI tropospheric O3 data over the Arctic during POLARCAT campaigns in 2008, Atmos. Chem. Phys., 12, 7371–7389, https://doi.org/10.5194/acp-12-7371-2012, 2012.

Rieder, H. E., Frossard, L., Ribatet, M., Staehelin, J., Maeder, J. A., Di Rocco, S., Davison, A. C., Peter, T., Weihs, P., and Holawe, F.: On the relationship between total ozone and atmospheric dynamics and chemistry at mid-latitudes – Part 2: The effects of the El Niño/Southern Oscillation, volcanic eruptions and contributions of atmospheric dynamics and chemistry to long-term total ozone changes, Atmos. Chem. Phys., 13, 165–179, https://doi.org/10.5194/acp-13-165-2013, 2013.

Rodgers, C. D.: Inverse Methods for Atmospheric Sounding: Theory and Practice, World Scientific, Series on Atmospheric, Oceanic and Planetary Physics, 2, edited by: Hackensack, N. J., World Scientific, Oxford, UK, 2000.

Safieddine, S., Boynard, A., Coheur, P.-F., Hurtmans, D., Pfister, G., Quennehen, B., Thomas, J. L., Raut, J.-C., Law, K. S., Klimont, Z., Hadji-Lazaro, J., George, M., and Clerbaux, C.: Summertime tropospheric ozone assessment over the Mediterranean region using the thermal infrared IASI/MetOp sounder and the WRF-Chem model, Atmos. Chem. Phys., 14, 10119–10131, https://doi.org/10.5194/acp-14-10119-2014, 2014.

Saunois, M., Emmons, L., Lamarque, J.-F., Tilmes, S., Wespes, C., Thouret, V., and Schultz, M.: Impact of sampling frequency in the analysis of tropospheric ozone observations, Atmos. Chem. Phys., 12, 6757–6773, https://doi.org/10.5194/acp-12-6757-2012, 2012.

Sauvage, B., Thouret, V., Thompson, A. M., Witte, J. C., Cammas, J. P., Nedelec, P., and Athier, G.: Enhanced view of the “tropical Atlantic ozone paradox” and “zonal wave one” from the in situ MOZAIC and SHADOZ data, J. Geophys. Res.-Atmos., 111, D01301, https://doi.org/10.1029/2005jd006241, 2006.

Sauvage, B., Martin, R. V., van Donkelaar, A., and Ziemke, J. R.: Quantification of the factors controlling tropical tropospheric ozone and the South Atlantic maximum, J. Geophys. Res.-Atmos., 112, D11309, https://doi.org/10.1029/2006jd008008, 2007.

Scannell, C., Hurtmans, D., Boynard, A., Hadji-Lazaro, J., George, M., Delcloo, A., Tuinder, O., Coheur, P.-F., and Clerbaux, C.: Antarctic ozone hole as observed by IASI/MetOp for 2008–2010, Atmos. Meas. Tech., 5, 123–139, https://doi.org/10.5194/amt-5-123-2012, 2012.

Simon, H., Reff, A., Wells, B., Xing, J., and Frank, N.: Ozone Trends Across the United States over a Period of Decreasing NOx and VOC Emissions, Environ. Sci. Technol., 49, 186–195, https://doi.org/10.1021/es504514z, 2015.

Stohl, A., Eckhardt, S., Forster, C., James, P., and Spichtinger, N.: On the pathways and timescales of intercontinental air pollution transport, J. Geophys. Res., 107, 4684, https://doi.org/10.1029/2001JD001396, 2002.

Tarasick, D., Galbally, I., Cooper, O. R., Ancellet, G., Leblanc, T., Wallington, T. J., Ziemke, J., Liu, X., Steinbacher, M., Stähelin, J., Vigouroux, C., Hannigan, J., García, O., Foret, G., Zanis, P., Weatherhead, E., Petropavlovskikh, I., Worden, H., Osman, M., Liu, J., Lin, M., Schultz, M., Granados-Muñoz, M., Thompson, A. M., Oltmans, S. J., Cuesta, J., Dufour, G., Thouret, V., Hassler, B., and Trickl, T.: Tropospheric Ozone Assessment Report: Tropospheric ozone observations – How well do we know tropospheric ozone changes?, Elem. Sci. Anth., in review, 2018.

Thompson, A. M., Witte, J. C., Smit, H. G. J., Oltmans, S. J., Johnson, B. J., Kirchhoff, V. W. J. H., and Schmidlin, F. J.: Southern Hemisphere Additional Ozonesondes (SHADOZ) 1998–2004 tropical ozone climatology: 3. Instrumentation, station-to-station variability, and evaluation with simulated flight profiles, J. Geophys. Res., 112, D03304, https://doi.org/10.1029/2005JD007042, 2007.

Thompson, D. W. J. and Wallace, J. M.: Annular modes in the extratropical circulation. Part I: month-to month variability, J. Climate, 13, 1000–1016, 2000.

Thouret, V., Cammas, J.-P., Sauvage, B., Athier, G., Zbinden, R., Nédélec, P., Simon, P., and Karcher, F.: Tropopause referenced ozone climatology and inter-annual variability (1994–2003) from the MOZAIC programme, Atmos. Chem. Phys., 6, 1033001051, https://doi.org/10.5194/acp-6-1033-2006, 2006.

Tiao, G. C., Reinsel, G. C., Xu, D., Pedrick, J. H., Zhu, X., Miller, A. J., DeLuisi, J. J., Mateer, C. L., and Wuebbles, D. J.: Effects of autocorrelation and temporal sampling schemes on estimates of trend and spatial correlation, J. Geophys. Res., 95, 20507–20517, 1990.

Tocquer, F., Barret, B., Mari, C., Le Flochmoen, E., Cammas, J. P., and Sauvage, B.: An upper tropospheric “ozone river” from Africa to India during the 2008 Asian post-monsoon season, Tellus B, 67, 25350, https://doi.org/10.3402/tellusb.v67.25350, 2015.

Valks, P., Hao, N., Gimeno Garcia, S., Loyola, D., Dameris, M., Jöckel, P., and Delcloo, A.: Tropical tropospheric ozone column retrieval for GOME-2, Atmos. Meas. Tech., 7, 2513–2530, https://doi.org/10.5194/amt-7-2513-2014, 2014.

van der A, R. J., Mijling, B., Ding, J., Koukouli, M. E., Liu, F., Li, Q., Mao, H., and Theys, N.: Cleaning up the air: effectiveness of air quality policy for SO2 and NOx emissions in China, Atmos. Chem. Phys., 17, 1775–1789, https://doi.org/10.5194/acp-17-1775-2017, 2017.

Verstraeten, W. W., Neu, J. L., Williams, J. E., Bowman, K. W., Worden, J. R., and Boersma, K. F.: Rapid increases in tropospheric ozone production and export from China, Nat. Geosci., 8, 690–695, https://doi.org/10.1038/NGEO2493, 2015.

Voulgarakis, A., Telford, P. J., Aghedo, A. M., Braesicke, P., Faluvegi, G., Abraham, N. L., Bowman, K. W., Pyle, J. A., and Shindell, D. T.: Global multi-year O3–CO correlation patterns from models and TES satellite observations, Atmos. Chem. Phys., 11, 5819–5838, https://doi.org/10.5194/acp-11-5819-2011, 2011a.

Voulgarakis, A., Hadjinicolaou, P., and Pyle, J. A.: Increases in global tropospheric ozone following an El Nino event: examining stratospheric ozone variability as a potential driver, Atmos. Sci. Lett., 12, 228–232, https://doi.org/10.1002/asl.318, 2011b.

Weatherhead, E. C., Reinsel, G. C., Tiao, G. C., Meng, X.-L., Choi, D., Cheang, W.-K., Keller, T., DeLuisi, 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.-Atmos., 103, 17149–17161, 1998.

Wespes, C., Hurtmans, D., Emmons, L. K., Safieddine, S., Clerbaux, C., Edwards, D. P., and Coheur, P.-F.: Ozone variability in the troposphere and the stratosphere from the first 6 years of IASI observations (2008–2013), Atmos. Chem. Phys., 16, 5721–5743, https://doi.org/10.5194/acp-16-5721-2016, 2016.

Wespes, C., Hurtmans, D., Clerbaux, C., and Coheur, P.-F.: O3 variability in the troposphere as observed by IASI over 2008–2016 – Contribution of atmospheric chemistry and dynamics, J. Geophys. Res.-Atmos., 122, 2429–2451, https://doi.org/10.1002/2016JD025875, 2017.

Wilson, R. C., Fleming, Z. L., Monks, P. S., Clain, G., Henne, S., Konovalov, I. B., Szopa, S., and Menut, L.: Have primary emission reduction measures reduced ozone across Europe? An analysis of European rural background ozone trends 1996–2005, Atmos. Chem. Phys., 12, 437–454, https://doi.org/10.5194/acp-12-437-2012, 2012.

Worden, J., Jiang, Z., Jones, D. B. A., Alvarado, M., Bowman, K., Frankenberg, C., Kort, E. A., Kulawik, S. S., Lee, M., Liu, J., Payne, V., Wecht, K., and Worden, H.: El Niño, the 2006 Indonesian peat fires, and the distribution of atmospheric methane, Geophys. Res. Lett., 40, 1–6, https://doi.org/10.1002/grl.50937, 2013.

Zbinden, R. M., Cammas, J.-P., Thouret, V., Nédélec, P., Karcher, F., and Simon, P.: Mid-latitude tropospheric ozone columns from the MOZAIC program: climatology and interannual variability, Atmos. Chem. Phys., 6, 1053–1073, https://doi.org/10.5194/acp-6-1053-2006, 2006.

Zeng, G., Morgenstern, O., Shiona, H., Thomas, A. J., Querel, R. R., and Nichol, S. E.: Attribution of recent ozone changes in the Southern Hemisphere mid-latitudes using statistical analysis and chemistry–climate model simulations, Atmos. Chem. Phys., 17, 10495–10513, https://doi.org/10.5194/acp-17-10495-2017, 2017.

Zhang, L., Jacob, D. J., Bowman, K. W., Logan, J. A., Turquety, S., Hudman, R. C., Li, Q., Beer, R., Worden, H. M., Worden, J. R., Rinsland, C. P., Kulawik, S. S., Lampel, M. C., Shephard, M. W., Fisher, B. M., Eldering, A., and Avery, M. A.: Ozone-CO correlations determined by the TES satellite instrument in continental outflow regions, Geophys. Res. Lett., 33, L18804, https://doi.org/10.1029/2006GL026399, 2006.

Zhang, L., Jacob, D. J., Boersma, K. F., Jaffe, D. A., Olson, J. R., Bowman, K. W., Worden, J. R., Thompson, A. M., Avery, M. A., Cohen, R. C., Dibb, J. E., Flock, F. M., Fuelberg, H. E., Huey, L. G., McMillan, W. W., Singh, H. B., and Weinheimer, A. J.: Transpacific transport of ozone pollution and the effect of recent Asian emission increases on air quality in North America: an integrated analysis using satellite, aircraft, ozonesonde, and surface observations, Atmos. Chem. Phys., 8, 6117–6136, https://doi.org/10.5194/acp-8-6117-2008, 2008.

Zhang, L., Li, Q. B., Murray, L. T., Luo, M., Liu, H., Jiang, J. H., Mao, Y., Chen, D., Gao, M., and Livesey, N.: A tropospheric ozone maximum over the equatorial Southern Indian Ocean, Atmos. Chem. Phys., 12, 4279–4296, https://doi.org/10.5194/acp-12-4279-2012, 2012.

Zhang, Y., Cooper, O. R., Gaudel, A., Thompson, A. M., Nédélec, P., Ogino, S.-Y., and West, J. J.: Tropospheric ozone change from 1980 to 2010 dominated by equatorward redistribution of emissions, Nat. Geosci., 9, 875–879, https://doi.org/10.1038/NGEO2827, 2016.

Zhao, B., Wang, S. X., Liu, H., Xu, J. Y., Fu, K., Klimont, Z., Hao, J. M., He, K. B., Cofala, J., and Amann, M.: NOx emissions in China: historical trends and future perspectives, Atmos. Chem. Phys., 13, 9869–9897, https://doi.org/10.5194/acp-13-9869-2013, 2013.

Ziemke, J. R., Douglass, A. R., Oman, L. D., Strahan, S. E., and Duncan, B. N.: Tropospheric ozone variability in the tropics from ENSO to MJO and shorter timescales, Atmos. Chem. Phys., 15, 8037–8049, https://doi.org/10.5194/acp-15-8037-2015, 2015.