Articles | Volume 19, issue 13
Atmos. Chem. Phys., 19, 8669–8686, 2019
Atmos. Chem. Phys., 19, 8669–8686, 2019

Research article 09 Jul 2019

Research article | 09 Jul 2019

Impact of El Niño–Southern Oscillation on the interannual variability of methane and tropospheric ozone

Impact of El Niño–Southern Oscillation on the interannual variability of methane and tropospheric ozone
Matthew J. Rowlinson1, Alexandru Rap1, Stephen R. Arnold1, Richard J. Pope1,2, Martyn P. Chipperfield1,2, Joe McNorton3, Piers Forster4, Hamish Gordon1, Kirsty J. Pringle1, Wuhu Feng1,5, Brian J. Kerridge6,7, Barry L. Latter6,7, and Richard Siddans6,7 Matthew J. Rowlinson et al.
  • 1Institute for Climate and Atmospheric Science, School of Earth and Environment, University of Leeds, Leeds, LS2 9JT, UK
  • 2National Centre for Earth Observation, University of Leeds, Leeds, LS2 9JT, UK
  • 3European Centre for Medium-Range Weather Forecasts, Reading, RG2 9AX, UK
  • 4Priestley International Centre for Climate, University of Leeds, LS2 9JT, Leeds, UK
  • 5National Centre for Atmospheric Science, University of Leeds, LS2 9JT, Leeds, UK
  • 6Remote Sensing Group, STFC Rutherford Appleton Laboratory, Harwell, Oxfordshire, OX11 0QX, UK
  • 7National Centre for Earth Observation, Harwell, Oxfordshire, OX11 0QX, UK

Correspondence: Matthew J. Rowlinson (


The interannual variability of the greenhouse gases methane (CH4) and tropospheric ozone (O3) is largely driven by natural variations in global emissions and meteorology. The El Niño–Southern Oscillation (ENSO) is known to influence fire occurrence, wetland emission and atmospheric circulation, affecting sources and sinks of CH4 and tropospheric O3, but there are still important uncertainties associated with the exact mechanism and magnitude of this effect. Here we use a modelling approach to investigate how fires and meteorology control the interannual variability of global carbon monoxide (CO), CH4 and O3 concentrations, particularly during large El Niño events. Using a three-dimensional chemical transport model (TOMCAT) coupled to a sophisticated aerosol microphysics scheme (GLOMAP) we simulate changes to CO, hydroxyl radical (OH) and O3 for the period 1997–2014. We then use an offline radiative transfer model to quantify the climate impact of changes to atmospheric composition as a result of specific drivers.

During the El Niño event of 1997–1998, there were increased emissions from biomass burning globally, causing global CO concentrations to increase by more than 40 %. This resulted in decreased global mass-weighted tropospheric OH concentrations of up to 9 % and a consequent 4 % increase in the CH4 atmospheric lifetime. The change in CH4 lifetime led to a 7.5 ppb yr−1 increase in the global mean CH4 growth rate in 1998. Therefore, biomass burning emission of CO could account for 72 % of the total effect of fire emissions on CH4 growth rate in 1998.

Our simulations indicate that variations in fire emissions and meteorology associated with El Niño have opposing impacts on tropospheric O3 burden. El Niño-related changes in atmospheric transport and humidity decrease global tropospheric O3 concentrations leading to a −0.03 W m−2 change in the O3 radiative effect (RE). However, enhanced fire emission of precursors such as nitrogen oxides (NOx) and CO increase O3 and lead to an O3 RE of 0.03 W m−2. While globally the two mechanisms nearly cancel out, causing only a small change in global mean O3 RE, the regional changes are large – up to −0.33 W m−2 with potentially important consequences for atmospheric heating and dynamics.

1 Introduction

In terms of radiative forcing, methane (CH4) is the second most important anthropogenically emitted greenhouse gas after CO2 (Myhre et al., 2013). Concentrations of CH4 have risen from approximately 722 ppb in 1750 to over 1850 ppb in 2018, an increase of more than 150 % (Dlugokencky, 2019). During this time period, CH4 has contributed an estimated radiative forcing (RF) of 0.48±0.05 W m−2, around 20 % of the total direct anthropogenic RF from greenhouse gases (Myhre et al., 2013). Furthermore, CH4 is a precursor of tropospheric ozone (O3), which is also a greenhouse gas responsible for a RF of 0.4±0.2 W m−2 since the pre-industrial era (Myhre et al., 2013), as well as a harmful pollutant that damages human health (Anenberg et al., 2010) and ecosystems (Sitch et al., 2007). While anthropogenic emissions have driven the long-term increase in CH4 concentrations, CH4 is also emitted from a range of natural sources, leading to strong interannual variability (IAV) (Bousquet et al., 2006; Dlugokencky et al., 2011; Nisbet et al., 2016). Understanding the mechanisms driving IAV is important for accurate predictions of future CH4 concentrations, especially in the context of anthropogenic emission reductions.

Previous studies indicate that although anthropogenic sources may contribute to seasonal variations in atmospheric CH4, natural sources are the primary drivers of IAV (Bousquet et al., 2006; Meng et al., 2015). Emissions from natural wetlands have been shown to be the dominant process, with emissions from fires and changes to the atmospheric sink also playing important roles (Bousquet et al., 2006; Chen and Prinn, 2006; Dlugokencky et al., 2011; Kirschke et al., 2013; McNorton et al., 2016a, 2018; Corbett et al., 2017). These natural sources are climate sensitive, so interannual changes to temperature and precipitation affect the amount of CH4 emitted into the atmosphere, as well as the spatial distribution (Zhu et al., 2017). A number of studies have found that biomass burning emissions are largely responsible for the IAV of carbon monoxide (CO) and also affect O3 concentrations (Granier et al., 2000; Monks et al., 2012; Voulgarakis et al., 2015); however, Szopa et al. (2007) suggested that meteorology is a more important driver of IAV for CO, explaining 50 %–90 % of IAV.

A major driver of climatic IAV is the El Niño–Southern Oscillation (ENSO) – a mode of climate variability originating in the Pacific Ocean with alternating warm (El Niño) and cold (La Niña) modes (McPhaden et al., 2006). Positive-phase El Niño events lead to warmer and drier conditions in much of the tropics, disrupting global circulation patterns and leading to widespread changes in fire occurrence, wetland emissions and atmospheric transport (Feely et al., 1987; Jones et al., 2001; McPhaden et al., 2006). These influences occur most strongly in the tropics but have global consequences (Jones et al., 2001). Global CH4 concentrations have been observed to increase significantly during El Niño events, with an especially strong signal during the 1997–1998 event when the CH4 growth rate was 12 ppb yr−1, almost triple the 1750–2018 mean annual growth rate (Rigby et al., 2008; Hodson et al., 2011). Due to the wide-ranging effects of El Niño and varied sources of CH4, there are multiple factors which could trigger the increase in CH4 growth rate. Chen and Prinn (2006) attributed the increase to anomalies in global wetland emissions; however, Zhu et al. (2017) estimated that although 49 % of the interannual variation in wetland emissions can be explained by ENSO, wetland emissions were significantly lower during El Niño, including the 1997–1998 event. Conversely, Schaefer et al. (2018) estimated that ENSO is responsible for up to 35 % of global CH4 variability, but the effect of wetland and biomass burning emission changes are dwarfed by processes affecting the OH sink. Bousquet et al. (2006) suggested that the increased CH4 growth rate during the 1997–1998 El Niño was primarily caused by abnormally large peat fires in Indonesia emitting huge amounts of CH4 while wetlands emissions remained stable (van der Werf et al., 2004; Butler et al., 2005; Bousquet et al., 2006).

In addition to direct emissions of CH4 from fires, it has been proposed that anomalously large CO emissions during enhanced El Niño fire events could explain the changes to CH4 growth rate (Butler et al., 2005; Bousquet et al., 2006). CO is emitted from biomass burning in much larger quantities than CH4 (20× larger) and its reaction with the hydroxyl radical (OH) is its primary atmospheric sink (Voulgarakis and Field, 2015). Abnormal increases in CO concentrations may suppress the availability of OH, thereby extending CH4 lifetime and increasing its growth rate during and following large fire events (Butler et al., 2005; Manning et al., 2005). The reaction of CH4 with OH is the largest term in the global CH4 budget, accounting for ∼90 % of its sink (McNorton et al., 2016a); therefore, even minor changes to OH caused by the presence of other compounds or changes to atmospheric transport and photolysis rates could have a large impact on CH4 growth rate (Dlugokencky et al., 2011). Butler et al. (2005) found that CO emissions suppressed OH concentrations by 2.2 % in 1997–1998, which accounted for 75 % of the observed change in CH4 concentration. Bousquet et al. (2006) also reported a weakened OH sink during this El Niño event.

Here we use a modelling approach to investigate how El Niño events affect global CH4, CO and tropospheric O3 concentrations through changes to fire occurrence and atmospheric transport. Using long-term simulations spanning multiple El Niño and La Niña events, we quantify the relative influence of changes to fire emissions and dynamical transport. We also differentiate between the effect of direct CH4 emissions from fires and the indirect effect via CO emissions and atmospheric chemistry changes.

2 Models and simulations

2.1 Model description

For this study we use the three-dimensional chemical transport model (TOMCAT) (Chipperfield, 2006) coupled to the GLOMAP global aerosol microphysics scheme (Mann et al., 2010). The version of TOMCAT-GLOMAP used here is a further development of that described by Monks et al. (2017). Cloud fields are now provided from the European Centre for Medium-Range Weather Forecasts (ECMWF) reanalyses (Dee et al., 2011), replacing the climatological clouds fields used previously from the International Satellite Cloud Climatology Project (ISCCP) (Rossow and Schiffer, 1999), leading to improved representation of photolysis. Other developments include updated emission inventories, the inclusion of CERN Cosmic Leaving Outdoor Droplets (CLOUD)-based new particle formation and the introduction of Mårtensson sea spray emissions (Gordon et al., 2017; Monks et al., 2017). The model is run at 2.8×2.8 horizontal resolution with 31 vertical levels from the surface to 10 hPa, driven by 6-hourly ECMWF ERA-Interim reanalyses. The planetary boundary layer (PBL) scheme is based on Holtslag and Boville (1993) and sea surface temperatures are from ECMWF reanalyses. ECMWF ERA-Interim reanalyses have been shown to have good skill in capturing Madden–Julian Oscillation (MJO) events, which in turn impact the onset of ENSO events (Dee et al., 2011), giving confidence that the model competently simulates El Niño meteorological conditions.

The tropospheric chemistry scheme used is as described in Monks et al. (2017), with anthropogenic emissions from the Monitoring Atmospheric Composition and Climate (MACCity) emissions inventories (Lamarque et al., 2010). Annually varying emission inventories are included for all fire-emitted gas-species and aerosol emissions, such as black carbon (BC). The Global Fire Emissions Database (GFED) used by TOMCAT-GLOMAP has been updated to version 4, with CO, nitrogen oxides (NOx) and volatile organic compound (VOC) emissions from fires (Randerson et al., 2017; Reddington et al., 2018). Monthly varying biogenic VOC emissions are from the MEGAN-MACC emissions inventory for reference year 2000, calculated from the Model of Emissions of Gases and Aerosols from Nature (MEGAN version 2) (Sindelarova et al., 2014). The CH4 inventory was produced by McNorton et al. (2016b), with wetland emissions derived from the Joint UK Land Environment Simulator (JULES) and biomass burning emissions from GFEDv4 (Randerson et al., 2017). These are then combined with anthropogenic emissions from EDGAR version 3.2, paddy field emissions from Yan et al. (2009) and termite, wild animal, mud volcano, hydrate and ocean emissions from Matthews and Fung (1987) (McNorton et al., 2016b). The global mean surface CH4 mixing ratio is scaled in TOMCAT-GLOMAP to a best estimate based on observed global surface mean concentration (McNorton et al., 2016a; Dlugokencky, 2019).

2.2 Radiative transfer model

Radiative effects of O3 changes are calculated using the O3 radiative kernel approach derived by Rap et al. (2015) using an offline version of the Edwards and Slingo (1996) radiative transfer model. This considers six bands in the shortwave (SW), nine bands in the longwave (LW) and uses a delta-Eddington two-stream scattering solver at all wavelengths (Rap et al., 2015). This version has been used extensively in conjunction with TOMCAT-GLOMAP for calculating radiative forcing from simulated distributions of several short-lived climate pollutants (SLCPs) including BC, O3 and CH4 (Spracklen et al., 2011; Riese et al., 2012; Rap et al., 2013, 2015; Richards et al., 2013).

2.3 Simulations

All simulations are performed for 1997–2014 with a 4-year spin-up through 1993–1996. The control run (CTRL) allows all emissions and meteorology to vary throughout the modelled period. GFED biomass burning emission inventories began in 1997; therefore, the 1993–1996 spin-up simulation uses repeating 1999 emissions instead, as the closest year of “average” emissions, having excluded 1997 and 1998 due to the exceptionally high emissions in those years (Schultz et al., 2008).

To test the impact of El Niño events on atmospheric chemistry, we also performed three perturbed simulations (Table 1). Where model simulations used “fixed” parameters in Table 1, the year 2013 emissions or meteorology are specified as invariant throughout the simulation. This year is chosen as the ENSO-neutral case, due to it being the least active ENSO year during 1997–2014, with a maximum bimonthly multivariate ENSO index (MEI) magnitude of −0.4 and the only year without a single MEI value that could be considered an active El Niño or La Niña (Wolter and Timlin, 1993, 1998). Throughout this study, an El Niño event was taken to be ongoing if the MEI was greater than +1.0. We perform a factorial analysis based on perturbed simulations in which we fix global biomass burning emissions (FIREFIX) or global meteorology (METFIX) to the “ENSO-neutral” case. An additional perturbed simulation was performed in order to examine the secondary impact of CO on CH4 via oxidation changes, where only CO emissions from biomass burning were fixed (COFIX).

Table 1Details of TOMCAT model simulations. All simulations are run for 1997–2014.

Download Print Version | Download XLSX

3 Model evaluation

We have conducted a comprehensive evaluation of the coupled TOMCAT-GLOMAP model using aircraft observations and data from ozone-sondes and satellites. In general, the model is able to capture absolute concentrations, global distribution and seasonal variations in major species including O3, CO and CH4. MOPITT satellite retrievals have been used to evaluate CO at 800 and 500 hPa (Emmons et al., 2004) and are shown in Figs. S1 and S2 in the Supplement, respectively, along with a description of the satellite product and the averaging kernels applied to the model output. TOMCAT performs similarly here, as in Monks et al. (2017), underestimating CO concentrations in the Northern Hemisphere (NH) while overestimating peak concentrations in biomass burning regions, with a maximum difference of ∼75 ppb (Figs. S1 and S2). However, TOMCAT is able to reproduce seasonal variations in CO and locates peak CO accurately over East Asia and Central Africa.

Simulated O3 concentrations from TOMCAT were also compared with satellite observations of lower tropospheric (0–6 km) O3 from the Ozone Monitoring Instrument (OMI). These data were provided by the Rutherford Appleton Laboratory (RAL; data version fv0214) using an optimal estimation retrieval scheme, which resolves O3 in the 0–6 km layer by exploiting information in the Hartley and Huggins UV bands. The scheme derives from the one discussed by Miles et al. (2015) for another UV sounder GOME-2. TOMCAT representation of O3 concentrations between 0 and 6 km in NH winter are slightly improved on the Monks et al. (2017) version, particularly in tropical and Southern Hemisphere (SH) concentrations (Fig. S3). However, there remains a general low bias in global O3 of up to 10 Dobson Units (DU) in winter in regions such as the southern Atlantic Ocean.

TOMCAT O3 has also been evaluated using sonde observations (Figs. 1 and S4) (Tilmes et al., 2012), with the model generally representing the vertical profiles, seasonal variation and absolute concentrations of O3 very well, with a normalized mean bias (NMB) of 1.1 % across all sites at 700–1000 hPa and 2.1 % at 300–700 hPa. The model capably simulates the seasonality of tropospheric O3 (Fig. 1), with a maximum seasonal bias of 6.3 % at 300–700 hPa in March–May. There is no apparent regional or latitudinal bias, although simulated concentrations are overestimated in India (Fig. S4). In addition, the TOMCAT-simulated global tropospheric burden of O3 in 2000 is 342 Tg, which falls within the range of published values (Table 2).

Figure 1Comparison of seasonal mean simulated O3 concentrations (ppb) against mean ozone-sonde observations from Tilmes et al. (2012), for the period 1995–2011. Panels (a–d) show mean concentrations at 700–1000 hPa across all sites, while panels (e–h) show mean concentrations at 300–700 hPa. Values in each panel are seasonal means in (from left to right) December–February (DJF), March–May (MAM), June–August (JJA) and September–November (SON). The red line represents the linear regression. Normalized mean bias (NMB) values between the model and observations are also shown.


Table 2Present day (2000) TOMCAT model diagnostics compared to previous model version from Monks et al. (2017) and other published values.

Download Print Version | Download XLSX

We have also assessed the capability of TOMCAT-GLOMAP to simulate observed responses to El Niño events. Ziemke et al. (2010) derived an O3 ENSO index using satellite observations, finding that for a +1 K change in the Niño 3.4 index, there was a 2.4 DU increase in the Ozone ENSO Index (OEI). In TOMCAT-GLOMAP, we calculate a 2.8 DU increase per +1 K in the Niño 3.4, indicating a slightly larger but comparable response to El Niño events. The regional response of tropospheric O3 to El Niño was evaluated against an analysis using various observations and a chemistry–climate model in Zhang et al. (2015). That study observed increased total O3 column in the North Pacific, southern USA, northeastern Africa and East Asia, with decreases over central Europe and the North Atlantic. All of these observed responses were present in TOMCAT-GLOMAP simulations, except with a slight increase in tropospheric ozone column (TOC) in central Europe and a simulated decrease in Western Europe and the eastern Atlantic (Fig. S5).

3.1 Aircraft observations

We compare annual mean simulated gas-phase species for 1999 against a climatological dataset of aircraft observations from 16 campaigns conducted with a broad spatial and temporal range from 1992 to 2001 (Emmons et al., 2010). While the comparison of observational data from intermittent aircraft campaigns does not offer a perfect comparison with the model simulated long-term mean concentrations, it allows evaluation of broad characteristics of a number of species over vertical profiles in many global regions. Figure 2 shows the comparison of simulated annual mean global concentrations of CO, CH4 and Peroxyacetyl nitrate (PAN), with aircraft observations at 0–2, 2–6 and 6–10 km. We have also calculated the normalized mean bias between the model and observations (Fig. S6). Full details of the aircraft measurement campaigns used can be found in Table S1 in the Supplement.

Figure 2Global mean volume mixing ratios of CO (ppb), CH4 (ppb) and PAN (ppt) from TOMCAT for the period 1993–2001 at 0–2 km (left panels), 2–6 km (middle panels) and 6–10 km (right panels). The filled circles show mean values from aircraft observation campaigns that took place between 1992 and 2001 (Table S1) (Emmons et al., 2010).


The model captures broad characteristics of spatial distribution for all species, simulating higher concentrations in polluted urban or biomass burning regions, with lower concentrations over ocean and in the SH. CO concentrations decrease with altitude but the largest values still occur around urban areas and burning regions, which can be seen in both model and aircraft concentrations. Consistent with the comparison with MOPITT satellite retrievals (Figs. S1 and S2), the model underestimates CO concentrations particularly near the surface, with a NMB of −11.1 %, −9.93 % and −0.25 % at 0–2, 2–6 and 6–10 km, respectively. Absolute concentrations of CH4 in TOMCAT simulations match aircraft data very well, although given the global mean surface concentration scaling, we expect the magnitude of CH4 to be well simulated. The latitudinal and vertical distributions are also well captured, giving confidence in the model transport and OH simulation. Aircraft observations show CH4 also decreases with altitude and the hemispheric disparity becomes more pronounced, with higher concentrations in the NH. For PAN concentrations, the simulated spatial distribution is broadly well captured, as is the increased concentration with altitude. There is a general low bias in absolute concentrations near the surface (NMB =-12.3 %), with a better comparison at 2–6 km (NMB =1.68 %) and overestimation at 6–10 km (NMB =18.17 %).

3.2 OH evaluation

Due to its very short lifetime, it is challenging to evaluate model-simulated OH over representative spatial and temporal scales. Here we follow the evaluation methodology recommended by Lawrence et al. (2001) of dividing tropospheric OH into 12 subdomains, from the surface to a climatologically derived tropopause. This method was also used to evaluate a previous version of TOMCAT(version 1.76) by Monks et al. (2017), allowing direct comparison. The evaluation is performed for the year 2000. Figure 3 shows our simulated OH compared to Monks et al. (2017), the ACCMIP model mean (Naik et al., 2013) and the Spivakovsky et al. (2000) OH dataset estimated from methyl chloroform observations.

Figure 3Annual zonal mean hydroxyl radical (OH) concentrations (×106 molec. cm−3) divided into 12 subdomains, as recommended by Lawrence et al. (2001). The simulated OH from this study is compared to a dataset estimated from methyl chloroform observations (Spivakovsky et al., 2000) and the Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP) multi-model mean (Naik et al., 2013). Results from a previous version of TOMCAT from Monks et al. (2017) are also shown. A climatological tropopause, indicated by the smooth black line near the top of each panel, has been used to remove stratospheric OH.


The models and observationally constrained distribution broadly agree in terms of the latitudinal spread of OH concentrations with a minimum in the SH and a maximum in the tropics; however, there is disagreement over the exact altitude of the maximum OH concentrations. In both versions of TOMCAT the highest concentration is between the surface and 750 hPa, while ACCMIP and Spivakovsky et al. (2000) find peak OH in the upper and mid-level troposphere, respectively. The updated cloud fields used in the current TOMCAT-GLOMAP version have slightly increased OH concentrations in the mid-level and upper domains compared to Monks et al. (2017) but concentrations remain significantly higher in the NH and surface domains than in other studies. In addition, our simulated NH:SH ratio of 1.48 in the current TOMCAT version remains substantially higher than in the ACCMIP models (1.28±0.1), indicating that TOMCAT photolysis rates and OH production in the NH are larger.

The total global tropospheric average OH in this version of TOMCAT is 1.04×106 molec. cm−3, a decrease from Monks et al. (2017) and within the range of other published values (Table 2). This is primarily due an updated treatment of clouds, in which climatological cloud fields have been replaced with cloud fraction from ECMWF reanalyses data, affecting photolysis rates. The tropospheric O3 burden of 342 Tg has increased relative to Monks et al. (2017) (331 Tg) and is within the range found in Wild (2007) (335±10 Tg) and ACCMIP models (337±23 Tg) (Young et al., 2013). Due to the simplified treatment of CH4, the scaling applied and its relatively long atmospheric lifetime, the total atmospheric lifetime cannot be determined from TOMCAT simulations. Instead a chemical lifetime due to reaction with OH is calculated from CH4 and OH burdens, disregarding stratospheric sinks and soil sinks (Fuglestvedt et al., 1999; Berntsen et al., 2005; Voulgarakis et al., 2013). The lifetime diagnosed from TOMCAT is 8.0 years, compared to the multi-model mean and range of 9.3±0.9 years from Voulgarakis et al. (2013). The shorter lifetime in TOMCAT is due to the overestimation of OH at the surface, particularly in the NH where CH4 concentrations are highest due to anthropogenic emissions.

4 Results and discussion

4.1 Impact of meteorology and fire emissions on trace gas interannual variability

First we examine the mechanisms controlling interannual variability of simulated tropospheric CO, O3 and mean OH. We use the difference between the control (CTRL) and the perturbed simulations with fixed fires (FIREFIX) and fixed meteorology (METFIX) to determine the driving cause of IAV. Of particular interest is the effect of the 1997–1998 El Niño event (henceforth referred to as 1997 El Niño) and how the prevailing mechanisms controlling IAV change during such events. To define El Niño events, we use the bimonthly multivariate ENSO index, which is calculated from six observed variables and standardized to accurately monitor ENSO occurrence (Wolter and Timlin, 1998, 2011).

Previous studies examining the dominant factor controlling global CO IAV have found contrasting results. Szopa et al. (2007) suggested that meteorology was the main driver, accounting for 50 %–90 % of IAV in the tropics. Conversely, a study by Monks et al. (2012) considered CO IAV in the Arctic, finding that biomass burning was the dominant driver with a strong correlation to El Niño. Voulgarakis et al. (2015) also suggested that biomass burning was the more important driver of IAV with only a small effect from meteorology. Some of these differences in results can be explained by Szopa et al. (2007) considering only surface CO, rather than the whole troposphere as in Voulgarakis et al. (2015). Here we also consider the whole of tropospheric CO and our results are in line with those from Voulgarakis et al. (2015). We find the dominant source of IAV across the entire period is emissions from biomass burning – indicated by the large difference between simulations CTRL and FIREFIX (Fig. 4a), with a small effect from meteorological changes (CTRL – METFIX). This effect was largest during the 1997 El Niño, where an increase in fire events increased CO concentrations by more than 40 %. Smaller increases of 5.8 % and 7.6 % occur during less extreme El Niño events of 2002/03 and 2006, respectively, with only a 1.8 % increase during the 2009/10 El Niño, indicating that El Niño only significantly impacts CO concentrations when there is an associated increase in global fire events.

Figure 4Time series of simulated differences (%) between the control and the fixed meteorology (CTRL – METFIX, blue line) and fixed fire emissions (CTRL – FIREFIX, purple line) simulations for the global tropospheric burden of (a) CO, (b) OH and (c) O3. The ENSO bimonthly mean multivariate index is plotted with the dashed red line using the right-hand y axis in each panel.


Expanding on the work of Voulgarakis et al. (2015), we analysed IAV using the coefficient of variation (CV), calculated as the multi-year standard deviation normalized by the mean (Fig. 5). The global annual mean CO IAV over the whole period is 11.0 % for the whole troposphere and 14.3 % for surface concentrations. This is in very good agreement with Voulgarakis et al. (2015), who calculated 10 % IAV; in fact, the comparison is even better when we consider the same time period (2005–2009) with our corresponding IAV estimate at 9.7 %. The slightly lower estimate here may be a result of the fixed-year biogenic volatile organic compound (BVOC) emissions, removing the effect of the IAV of biogenic emissions on CO IAV. BVOC oxidation is estimated to contribute 15 % of the total source of CO (Duncan et al., 2007); however, the IAV of BVOC emissions has been found to be relatively small, ∼2–4 % (Naik et al., 2004; Lathière et al., 2005). Despite good global comparison with Voulgarakis et al. (2015), there are regional differences; CO IAV from TOMCAT is much larger in high-latitude boreal regions. This is likely due to the difference in the period studied, meaning this study includes additional extreme events including unusually large Russia boreal wildfires in 2010 and 2012 (Gorchakov et al., 2014; Kozlov et al., 2014). Infrequent and extreme events such as these significantly increase IAV.

Figure 5The calculated interannual variability (coefficient of variation) of CO over the period 1997–2014 for September–October (a, c, e) and March–April (b, d, f) from (a,  b) the control simulation (CTRL), (c, d) fixed meteorology (METFIX) and (e, f) fixed fire emissions (FIREFIX).


CO IAV is significantly greater in September–October, with peaks in known fire regions such as tropical South America, Africa, Southeast Asia and boreal forests. This indicates a strong contribution of fire emissions to IAV, especially from Indonesia (Fig. 4a), as also suggested by previous studies (Monks et al., 2012; Huang et al., 2014; Voulgarakis et al., 2015). In the FIREFIX simulation IAV is ∼55 % of the CTRL value, showing a large reduction in variability when interannual variability in fire emissions is removed. The IAV in March–April is significantly smaller than September–October as this period is outside the primary fire season for South America and Eurasia, although hotspots remain in Southeast Asia and Africa where fires commonly occur in March–April (van der Werf et al., 2017). Meteorology and atmospheric transport changes are most important in Africa in September–October and Indonesia in March–April (Fig. 5c, d). Fire emissions occur in these regions but the meteorological effects are important sources of IAV. This is in good agreement with Voulgarakis et al. (2015), who found that with fixed biomass burning emissions, high IAV remained over Africa during December–January, and Huang et al. (2014), who found CO over Central Africa correlated more closely with ice water content than CO emissions due to increased convective transport. However, the overall effect of meteorology on global IAV found here is much smaller than the 50 %–90 % suggested by Szopa et al. (2007): when we consider only surface CO over the same period, fixing meteorology decreases the mean CO IAV by just 5 %.

The IAV of OH and O3 have more complex contributions from fire emissions and meteorology (Fig. 4b, c). For both species, meteorology is the dominant cause of variability for the majority of the period, indicated by, on average, greater deviation from CTRL in METFIX simulation than FIREFIX, including during El Niño events outside of the 1997 El Niño, such as in 2006. Our results compare well to Inness et al. (2015), who also found that changes to tropospheric O3 during El Niño were driven by a combination of emissions and atmospheric dynamics. This is also in agreement with Doherty et al. (2006), where a strong correlation was found between ENSO meteorology and global O3 burden, albeit with a lag period of several months. Various meteorological variables are known to affect OH and O3 variability, including humidity, clouds and temperature (Stevenson et al., 2005; Holmes et al., 2013; Nicely et al., 2018). OH variability is particularly sensitive to changes in lightning NOx production which decreases during El Niño conditions (Turner et al., 2018). Murray et al. (2014) also examined factors affecting OH variability since the last glacial maximum, finding tropospheric water vapour, overhead stratospheric O3 and lightning NOx to be key controlling factors. Furthermore, circulation changes during El Niño events have been linked to lower stratospheric O3 variability (Zhang et al., 2015; Manatsa and Mukwada, 2017), which in turn influences tropospheric OH and O3 concentrations (Holmes et al., 2013; Murray et al., 2014). Despite the importance of meteorological drivers, we find that fire emissions are the dominant cause of variation in both OH and O3 during the 1997 El Niño, increasing global tropospheric O3 burden by up to ∼7 % and decreasing tropospheric OH by up to ∼6 %. This result is supported by several other studies, which have found that during large fire events such as that caused by the 1997 El Niño, fire emissions substantially decrease tropospheric OH and increase tropospheric O3 (Hauglustaine et al., 1999; Sudo and Takahashi, 2001; Holmes et al., 2013). Our results indicate that while meteorology is generally the most important driver of IAV in global tropospheric OH and O3, fire emissions can also play a key role and become the dominant driver when there are particularly large fire emissions related to El Niño.

Figure 6 shows the IAV of O3, supporting the analysis of Fig. 4 that also suggests meteorology is the dominant process in controlling IAV. METFIX-simulated IAV differs substantially from the CTRL, with much lower IAV in September–October (33 % decrease) and in March–April (42 % decrease) when meteorology is repeated. However, in the METFIX run there remain peaks in variability in close proximity to regions with large biomass burning emissions, demonstrating the significant contribution from fire emissions. In the FIREFIX simulation the distribution of IAV is broadly similar to the CTRL simulation and only shows a small change in global mean CV, indicating that fire emissions have less control on O3 IAV. These results are again comparable to Voulgarakis et al. (2015) as the distribution of O3 IAV in both CTRL and FIREFIX simulations is similar, despite slightly larger values of variation due to differing time period.

Figure 6The calculated interannual variability (coefficient of variation) of ozone over the period 1997–2014 for September–October (a, c, e) and March–April (a, c, e) from (a, b) the control simulation (CTRL), (c, d) fixed meteorology (METFIX) and (e, f) fixed fire emissions (FIREFIX).


4.2 Indirect effect of CO on oxidation and lifetime of CH4

The COFIX sensitivity experiment was conducted to determine the indirect influence of CO emissions on CH4 variability through changes in tropospheric OH concentrations. Figure 7a shows the difference in COFIX monthly mean OH concentrations from the control experiment, compared to that from the METFIX and FIREFIX simulations. When CO emissions from biomass burning are fixed, OH concentrations are consistently higher than in the CTRL simulation. This indicates that high CO emissions decrease global mean tropospheric OH. The greatest impact is during the 1997 El Niño, where CO emissions were abnormally large, suppressing mass weighted global monthly mean OH concentrations by up to 9 %. The mean effect on OH over the 1997 El Niño of −3.6 % is comparable to that simulated by Butler et al. (2005), who also found an increase in CO resulted in a change in OH of −2.2 %. Duncan et al. (2003) found a similar magnitude response in OH to the Indonesian wildfires in 1997 of between −2.1 % and −6.8 %. The suppression of OH concentrations due to CO emission is also simulated to a lesser degree in the 2003 and 2006 El Niño events but is absent in the 2010 El Niño as this event had little impact on global fire occurrence (Randerson et al., 2017). The effect of fixing only CO from fires is greater than the effect of fixing all fire emissions due to co-emitted species such as NOx, which act to increase OH concentrations.

Figure 7Time series of (a) the change (%) in mass-weighted tropospheric OH, (b) change (%) in CH4 lifetime and (c) resultant change (ppb) in annual CH4 growth rate calculated using an offline box model. The ENSO bimonthly mean multivariate index is plotted with the dashed red line using the right-hand y axis in (a).


As OH is also the primary sink of CH4 (∼90 %) (McNorton et al., 2016a), another effect of the decrease in OH due to CO emissions is to weaken the sink of CH4, increasing its atmospheric lifetime. The magnitude of this can be seen in Fig. 7b; the COFIX simulation indicates that CO emissions from fires extended CH4 atmospheric lifetime by more than 4 % during the 1997 El Niño. Fixing all fire emissions also enhances CH4 lifetime by around 2 %. Increasing the lifetime of a species increases its concentration in steady-state equilibrium. Due to the scaling applied to CH4 in TOMCAT we are unable to directly calculate the response in CH4 growth rate from TOMCAT, as simulated global mean surface CH4 concentrations are nudged to the observed value. Therefore, to determine the impact of the change to OH on CH4 concentrations, we used a simple global box model. This box model is similar to that described in McNorton et al. (2016a), which was found to compare well with other global and 12-box CH4 models (Rigby et al., 2013; McNorton et al., 2016a). In this case, the box model used monthly mean tropospheric OH concentrations and CH4 emissions for each simulation, while assuming constant temperature to calculate the effect of changing OH on global mean surface CH4. A fixed temperature was used as varying the temperature has been found to have a relatively small impact on derived CH4 concentrations (McNorton et al., 2016a). The impact of fire emissions on the CH4 growth rate is greatest in 1998, where all emissions from fires increased global CH4 by 10.5 ppb (Fig. 7c). Analysis of the COFIX simulation demonstrates that up to 7.5 ppb (72 %) of that change could have been caused by the release of CO alone and its role as a sink for OH. The effect on growth rate in the FIREFIX simulation is larger than in the COFIX despite a greater effect on CH4 lifetime from the COFIX, due to directly emitted CH4 varying with El Niño conditions in the COFIX simulation and not in FIREFIX. The influence of CO emissions on CH4 growth rate calculated here is smaller than in Butler et al. (2005), despite a much larger effect on tropospheric OH. The radiative effect of the change to CH4 from CO emitted from biomass burning alone in 1998 is 0.004 W m−2, calculated using updated expressions from Etminan et al. (2016).

4.3 Limiting factors of O3 production

In this section we examine trends and the impact of El Niño on the production of tropospheric O3. El Niño is known to have a large effect on tropospheric O3 precursors such as CO and NOx; therefore, examining O3 production regimes during El Niño can provide insights into the main mechanism responsible for the observed changes in tropospheric O3. The ratio between formaldehyde (HCHO) and nitrogen dioxide (NO2) concentrations can be used to indicate the limiting factor for tropospheric O3 production (Duncan et al., 2010). Ratios smaller than 1 indicate that removing VOCs will decrease tropospheric O3 formation (i.e. a VOC-limited regime), while ratios larger than 2 indicate that removing NOx will reduce O3 (i.e. a NOx-limited regime). Ratios of 1–2 indicate that both NOx and VOC reductions could decrease O3 (i.e. a “both-limited” regime). Here we apply this methodology to determine the changes to this ratio from 1997 to 2014 and dependence of O3 formation during the 1997 El Niño event. We compare the early period mean (1999–2003) to the end period mean (2010–2014) to determine whether significant changes have occurred over the 18-year period and compared mean El Niño conditions to both.

Figure 8Mean ratio of simulated tropospheric column HCHO to NO2 amounts for (a) the beginning of model period (1999–2003), (b) the end of model period (2010–2014) and (c) during all El Niño events. Panels (d, e) show the difference during El Niño from the 5-year mean values in panels (a, b), respectively.


In general, the SH and tropical regions have very high ratios, meaning they are strongly NOx-limited (Fig. 8). The NH is also predominantly NOx-limited, although less robustly, and polluted regions tend to be either VOC-limited or both-limited regimes. The ratio is largely constant across the modelled period; however, there are some significant shifts, such as in India, which was once solely NOx-limited, becoming increasing VOC-limited due to increased NOx pollution (Hilboll et al., 2017). This shift in the spatial distribution of O3 precursor emissions to lower latitudes leads to increased tropospheric O3 production proportional to total emissions (Zhang et al., 2016).

During El Niño there are large changes, increasing the ratio and therefore the NOx limitation by more than 40 % in the tropical Pacific. Significant changes to the ratio were also found in biomass burning regions of South America and Southeast Asia. This is due to the increase in NOx emissions in larger fire seasons associated with El Niño. However, these regions are already very heavily NOx-limited due to high VOC emissions in forest regions, meaning that although the shift in HCHO∕NO2 ratio during El Niño is large, it is not substantial enough to alter the limiting factor for formation of tropospheric O3 from one regime to another. Over India, El Niño conditions inhibit the trend towards a both-limited regime, as the NOx-limited regime continues to dominate throughout.

4.4 Impact on tropospheric ozone and radiative effects

The 1997 El Niño significantly altered the vertical distribution of O3 in the troposphere, increasing O3 concentrations in the NH while decreasing in the SH and tropics with an overall decrease in tropospheric O3 of −0.82 % compared to the 1997–2014 mean (Fig. 9a). In the CTRL simulation there is decreased O3 in the tropical upper troposphere, possibly related to increased convection over the eastern Pacific (Oman et al., 2013; Neu et al., 2014). We also simulate large increases in the midlatitude upper troposphere of both hemispheres in the CTRL and FIREFIX simulations but not in METFIX, implying that this is produced by El Niño-associated meteorological processes which promote intrusion of stratospheric air into the troposphere. These positive anomalies were also observed in Oman et al. (2013) and Zeng and Pyle (2005), attributed to El Niño influence on circulation patterns and enhanced stratosphere–troposphere exchange.

Figure 9Latitude–pressure cross sections of the percentage difference in O3 concentrations during the 1997 El Niño event compared to 1997–2014 period mean for the TOMCAT simulations: (a) CTRL, (b) METFIX and (c) FIREFIX simulations.


In general, the METFIX run simulates higher O3 concentrations in the NH than the period mean and lower concentrations in the SH (Fig. 9b). This hemispherical shift is also present in the CTRL and FIREFIX simulations but with greater negative O3 anomalies in the SH. The simulated NH increases in the CTRL simulation correspond to other studies of the 1997 El Niño (Koumoutsaris et al., 2008), while Oman et al. (2013) similarly reported negative O3 anomalies in the SH during El Niño. Large increases in tropospheric O3 in the western Pacific, Indian Ocean and Europe contribute to the increase in O3 in the NH, despite decreased O3 in the eastern Pacific (Chandra et al., 1998; Koumoutsaris et al., 2008; Oman et al., 2011).

Figure 10Tropospheric O3 radiative effects (W m−2) from the TOMCAT simulations: (a) control (CTRL), (b) fixed meteorology and fire emissions (BOTHFIX), (c) fixed meteorology only (METFIX) and (d) fixed fire emissions only (FIREFIX). Panels (e–g) show percentage differences between the control and the three perturbed simulations.


There is an overall increase in O3 (∼2 %) when meteorology was fixed to an ENSO-neutral year (i.e. 2013), meaning that meteorology during the 1997 El Niño caused a decrease in tropospheric O3 concentrations despite large increases in O3 in regions of the upper troposphere due to stratospheric intrusion. During the 1997 El Niño we find a 0.4 % increase in global tropospheric humidity compared to the period mean. This is likely partly responsible for the general decrease in O3 due to meteorology, as increased humidity enhances O3 loss (Stevenson et al., 2000; Isaksen et al., 2009; Kawase et al., 2011). Changes to transport and distribution of O3 will also impact how efficiently tropospheric O3 is produced and lost.

The similarities between the tropospheric O3 distribution in the CTRL and FIREFIX simulations show that fire emissions have a relatively small impact on the global distribution of O3 but do affect absolute values, as concentrations in the FIREFIX run are significantly lower in the tropics. This is likely because of the removal of large emissions of O3 precursors in that latitude band when fire emissions are fixed to a non-El Niño year, as several studies have found that enhanced fires in 1997 El Niño increased tropospheric O3 in the region (Chandra et al., 1998; Thompson et al., 2001; Doherty et al., 2006; Oman et al., 2013).

Figure 10 shows the tropospheric O3 radiative effect (RE) during the 1997 El Niño in each TOMCAT simulation, calculated using the Rap et al. (2015) tropospheric O3 radiative kernel. Consistent with the relative changes in O3 concentration, fire emissions and meteorology have contrasting effects on O3 RE. When isolated, these effects are opposite and almost equal: fire emissions increase O3 RE by 0.031 W m−2, while meteorology decreases by −0.030 W m−2. We performed an additional simulation to determine the effect of these factors occurring simultaneously (BOTHFIX) and found the increasing effect from fire emissions to be dominant over the decreasing effect from meteorology, leading to an overall increase in global mean O3 RE of 0.015 W m−2. The effect of fire emissions occurs almost entirely over Indonesia and the eastern Indian Ocean where the large influx of NOx, CO and CH4 from fire emissions during the 1997 El Niño causes large regional increases in tropospheric O3. This increase, also observed in Chandra et al. (1998), causes a regional RE of up to 0.17 W m−2. Meteorology has more varied impacts during El Niño, causing large decreases in O3 RE over the central Pacific Ocean (-0.36 W m−2) but also increases in the midlatitudes of the Pacific Ocean (∼0.33 W m−2). Globally the mean change to O3 RE is small, around 0.015 W m−2, but large regional changes have the potential to significantly alter atmospheric heating and dynamics.

5 Summary and conclusions

Global model simulations using annually invariant meteorology and fire emissions were performed for the period 1997–2014 in order to determine their relative impacts on the IAV of O3 and CH4, particularly during El Niño events. The TOMCAT-GLOMAP model used has been updated from that described by Monks et al. (2017), with improved cloud and photolysis representation and the introduction of Mårtensson sea spray emissions (Gordon et al., 2017). Model simulations were evaluated for a number of chemical species (O3, CH4, NOx, CO), with observations from aircraft, satellites and ozone-sondes. In general, the model shows a good agreement with observed values, although with some regional biases. Differences between the model and observations may be due to a number of factors, such as the relatively coarse model resolution, uncertainties in the model emission inventories and errors in observations. However, good overall agreement of model simulations with different observations, including the ability of the model to simulate the observed atmospheric responses to El Niño events (i.e. OEI change of 2.8 DU compared to 2.4 DU in Ziemke et al., 2010), provides confidence in model performance and results.

We find that the IAV of global CO concentrations is large and is primarily controlled by fire emissions over the modelled period. Exceptionally large CO emissions linked to El Niño in 1997 led to a decrease in OH concentrations of ∼9 %, which subsequently increased CH4 lifetime by ∼4 %. Using a box model we quantify the isolated impact of this change in atmospheric chemistry on global CH4 growth rate to be 7.75 ppb, ∼75 % of the total effect of fires. This effect, combined with concurrent direct CH4 emission from fires, explains the observed changes to CH4 growth rate during the 1997 El Niño.

Variability of oxidants O3 and OH is far more dependent on meteorology than fire emissions, except during very large El Niño events, such as in 1997 and 1998, when fires become dominant in terms of total tropospheric burden, although meteorology still controls distribution. The change to tropospheric O3 concentrations during El Niño has increased O3 RE by 0.17 W m−2 over Southeast Asia and decreased by 0.36 W m−2 over the central Pacific. The global mean O3 RE change due to 1997 El Niño meteorology and fires is an increase of 0.015 W m−2, as emissions of O3 precursors from fires causes increased O3. El Niño also causes significant shifts in the ratio of HCHO:NOx – an indicator of O3 production regime – but most significantly in the tropics, which are heavily NOx-limited, so this change does not cause a regime shift.

This work has shown that El Niño events significantly affect the variability of two important drivers of anthropogenic climate change. Further research into how El Niño events, with their associated effect on fire emissions, are likely to change in a warming climate is required to understand how these links between ENSO, CH4 and O3 may influence future climate change mitigation attempts.

Data availability

All data can be accessed upon request to the corresponding author.


The supplement related to this article is available online at:

Author contributions

MJR, AR and SRA conceptualized the study and planned the experiments. MPC, KJP, HG, WF and MJR developed and evaluated the version of TOMCAT-GLOMAP used here. BJK, BLL and RS provided the satellite retrievals for the O3 comparison, which was conducted by RJP. MPC and JM provided assistance and advice for the CH4 box model. MJR performed the TOMCAT model runs, SOCRATES and box model calculations. MJR analysed the results with help from AR and SRA. MJR compiled results and prepared the manuscript. All co-authors contributed to the final version with comments.

Competing interests

The authors declare that they have no conflict of interest.


This work was undertaken using the ARCHER UK National Supercomputing Service ( and ARC3, part of the High Performance Computing facilities at the University of Leeds, UK.

Financial support

This research has been supported by the Leeds York NERC DTP (grant no. NE/L002574/1).

Review statement

This paper was edited by Tim Butler and reviewed by two anonymous referees.


Anenberg, S. C., Horowitz, L. W., Tong, D. Q., and West, J. J.: An Estimate of the Global Burden of Anthropogenic Ozone and Fine Particulate Matter on Premature Human Mortality Using Atmospheric Modeling, Environ. Health Persp., 118, 1189–1195,, 2010. 

Berntsen, T. K., Fuglestvedt, J. S., Joshi, M. M., Shine, K. P., Stuber, N., Ponater, M., Sausen, R., Hauglustaine, D. A., and Li, L.: Response of climate to regional emissions of ozone precursors: sensitivities and warming potentials, Tellus B, 57, 283–304,, 2005. 

Bousquet, P., Hauglustaine, D. A., Peylin, P., Carouge, C., and Ciais, P.: Two decades of OH variability as inferred by an inversion of atmospheric transport and chemistry of methyl chloroform, Atmos. Chem. Phys., 5, 2635–2656,, 2005. 

Bousquet, P., Ciais, P., Miller, J. B., Dlugokencky, E. J., Hauglustaine, D. A., Prigent, C., Van der Werf, G. R., Peylin, P., Brunke, E. G., Carouge, C., Langenfelds, R. L., Lathière, J., Papa, F., Ramonet, M., Schmidt, M., Steele, L. P., Tyler, S. C., and White, J.: Contribution of anthropogenic and natural sources to atmospheric methane variability, Nature, 443, 439–443,, 2006. 

Butler, T. M., Rayner, P. J., Simmonds, I., and Lawrence, M. G.: Simultaneous mass balance inverse modeling of methane and carbon monoxide, J. Geophys. Res., 110, D21310,, 2005. 

Chandra, S., Ziemke, J. R., Min, W., and Read, W. G.: Effects of 1997–1998 El Niño on tropospheric ozone and water vapor, Geophys. Res. Lett., 25, 3867–3870,, 1998. 

Chen, Y.-H. and Prinn, R. G.: Estimation of atmospheric methane emissions between 1996 and 2001 using a three-dimensional global chemical transport model, J. Geophys. Res., 111, D10307,, 2006. 

Chipperfield, M. P.: New version of the TOMCAT/SLIMCAT off-line chemical transport model: Intercomparison of stratospheric tracer experiments, Q. J. Roy. Meteorol. Soc., 132, 1179–1203,, 2006. 

Corbett, A., Jiang, X., Xiong, X., Kao, A., and Li, L.: Modulation of midtropospheric methane by El Niño, Earth and Space Science, 4, 590–596,, 2017. 

Dee, D. P., Uppala, S. M., Simmons, A. J., Berrisford, P., Poli, P., Kobayashi, S., Andrae, U., Balmaseda, M. A., Balsamo, G., Bauer, P., Bechtold, P., Beljaars, A. C. M., van de Berg, L., Bidlot, J., Bormann, N., Delsol, C., Dragani, R., Fuentes, M., Geer, A. J., Haimberger, L., Healy, S. B., Hersbach, H., Hólm, E. V., Isaksen, L., Kållberg, P., Köhler, M., Matricardi, M., McNally, A. P., Monge-Sanz, B. M., Morcrette, J.-J., Park, B.-K., Peubey, C., de Rosnay, P., Tavolato, C., Thépaut, J.-N., and Vitart, F.: The ERA-Interim reanalysis: configuration and performance of the data assimilation system, Q. J. Roy. Meteorol. Soc., 137, 553–597,, 2011. 

Dlugokencky, E.: NOAA/ESRL, available at:, last access: 29 April 2019. 

Dlugokencky, E. J., Nisbet, E. G., Fisher, R., and Lowry, D.: Global atmospheric methane: budget, changes and dangers, Philos. T. Roy. Soc. A, 369, 2058–2072,, 2011. 

Doherty, R. M., Stevenson, D. S., Johnson, C. E., Collins, W. J., and Sanderson, M. G.: Tropospheric ozone and El Niño–Southern Oscillation: Influence of atmospheric dynamics, biomass burning emissions, and future climate change, 111,, 2006. 

Duncan, B., Yoshida, Y., R. Olson, J., Sillman, S., Martin, R., Lamsal, L., Hu, Y., E. Pickering, K., Retscher, C., and J. Allen, D.: Application of OMI observations to a space-based indicator of NOx and VOC controls on surface ozone formation, Atmos. Environ., 44, 2213–2223, 2010. 

Duncan, B. N., Bey, I., Chin, M., Mickley, L. J., Fairlie, T. D., Martin, R. V., and Matsueda, H.: Indonesian wildfires of 1997: Impact on tropospheric chemistry, J. Geophys. Res., 108, 4458,, 2003. 

Duncan, B. N., Logan, J. A., Bey, I., Megretskaia, I. A., Yantosca, R. M., Novelli, P. C., Jones, N. B., and Rinsland, C. P.: Global budget of CO, 1988–1997: Source estimates and validation with a global model, J. Geophys. Res., 112, D22301,, 2007. 

Edwards, J. M. and Slingo, A.: Studies with a flexible new radiation code. 1. Choosing a configuration for a large-scale model, Q. J. Roy. Meteor. Soc., 122, 689–719,, 1996. 

Emmons, L. K., Deeter, M. N., Gille, J. C., Edwards, D. P., Attié, J.-L., Warner, J., Ziskin, D., Francis, G., Khattatov, B., Yudin, V., Lamarque, J.-F., Ho, S.-P., Mao, D., Chen, J. S., Drummond, J., Novelli, P., Sachse, G., Coffey, M. T., Hannigan, J. W., Gerbig, C., Kawakami, S., Kondo, Y., Takegawa, N., Schlager, H., Baehr, J., and Ziereis, H.: Validation of Measurements of Pollution in the Troposphere (MOPITT) CO retrievals with aircraft in situ profiles, J. Geophys. Res., 109, D03309,, 2004. 

Emmons, L. K., Walters, S., Hess, P. G., Lamarque, J.-F., Pfister, G. G., Fillmore, D., Granier, C., Guenther, A., Kinnison, D., Laepple, T., Orlando, J., Tie, X., Tyndall, G., Wiedinmyer, C., Baughcum, S. L., and Kloster, S.: Description and evaluation of the Model for Ozone and Related chemical Tracers, version 4 (MOZART-4), Geosci. Model Dev., 3, 43–67,, 2010. 

Etminan, M., Myhre, G., Highwood, E. J., and Shine, K. P.: Radiative forcing of carbon dioxide, methane, and nitrous oxide: A significant revision of the methane radiative forcing, Geophys. Res. Lett., 43, 12614–612623,, 2016. 

Feely, R. A., Gammon, R. H., Taft, B. A., Pullen, P. E., Waterman, L. S., Conway, T. J., Gendron, J. F., and Wisegarver, D. P.: Distribution of chemical tracers in the eastern equatorial Pacific during and after the 1982–1983 El Niño/Southern Oscillation event, J. Geophys. Res., 92, 6545–6558,, 1987. 

Fuglestvedt, J. S., Berntsen, T. K., Isaksen, I. S. A., Mao, H., Liang, X.-Z., and Wang, W.-C.: Climatic forcing of nitrogen oxides through changes in tropospheric ozone and methane; global 3D model studies, Atmos. Environ., 33, 961–977,, 1999. 

Gorchakov, G. I., Sitnov, S. A., Sviridenkov, M. A., Semoutnikova, E. G., Emilenko, A. S., Isakov, A. A., Kopeikin, V. M., Karpov, A. V., Gorchakova, I. A., Verichev, K. S., Kurbatov, G. A., and Ponomareva, T. Y.: Satellite and ground-based monitoring of smoke in the atmosphere during the summer wildfires in European Russia in 2010 and Siberia in 2012, Int. J. Remote Sens., 35, 5698–5721,, 2014. 

Gordon, H., Kirkby, J., Baltensperger, U., Bianchi, F., Breitenlechner, M., Curtius, J., Dias, A., Dommen, J., Donahue, N. M., Dunne, E. M., Duplissy, J., Ehrhart, S., Flagan, R. C., Frege, C., Fuchs, C., Hansel, A., Hoyle, C. R., Kulmala, M., Kürten, A., Lehtipalo, K., Makhmutov, V., Molteni, U., Rissanen, M. P., Stozkhov, Y., Tröstl, J., Tsagkogeorgas, G., Wagner, R., Williamson, C., Wimmer, D., Winkler, P. M., Yan, C., and Carslaw, K. S.: Causes and importance of new particle formation in the present-day and preindustrial atmospheres, J. Geophys. Res.-Atmos., 122, 8739–8760,, 2017. 

Granier, C., Müller, J.-F., and Brasseur, G.: The Impact of Biomass Burning on the Global Budget of Ozone and Ozone Precursors, in: Biomass Burning and Its Inter-Relationships with the Climate System, edited by: Innes, J. L., Beniston, M., and Verstraete, M. M., Springer, Dordrecht, 2000. 

Hauglustaine, D. A., Brasseur, G. P., and Levine, J. S.: A sensitivity simulation of tropospheric ozone changes due to the 1997 Indonesian fire emissions, Geophys. Res. Lett., 26, 3305–3308,, 1999. 

Hilboll, A., Richter, A., and Burrows, J. P.: NO2 pollution over India observed from space – the impact of rapid economic growth, and a recent decline, Atmos. Chem. Phys. Discuss.,, in review, 2017. 

Hodson, E. L., Poulter, B., Zimmermann, N. E., Prigent, C., and Kaplan, J. O.: The El Niño–Southern Oscillation and wetland methane interannual variability, Geophys. Res. Lett., 38, L08810,, 2011. 

Holmes, C. D., Prather, M. J., Søvde, O. A., and Myhre, G.: Future methane, hydroxyl, and their uncertainties: key climate and emission parameters for future predictions, Atmos. Chem. Phys., 13, 285–302,, 2013. 

Holtslag, A. A. M. and Boville, B. A.: Local Versus Nonlocal Boundary-Layer Diffusion in a Global Climate Model, J. Climate, 6, 1825–1842,<1825:Lvnbld>2.0.Co;2, 1993. 

Huang, L., Fu, R., and Jiang, J. H.: Impacts of fire emissions and transport pathways on the interannual variation of CO in the tropical upper troposphere, Atmos. Chem. Phys., 14, 4087–4099,, 2014. 

Inness, A., Benedetti, A., Flemming, J., Huijnen, V., Kaiser, J. W., Parrington, M., and Remy, S.: The ENSO signal in atmospheric composition fields: emission-driven versus dynamically induced changes, Atmos. Chem. Phys., 15, 9083–9097,, 2015. 

Isaksen, I. S. A., Granier, C., Myhre, G., Berntsen, T. K., Dalsøren, S. B., Gauss, M., Klimont, Z., Benestad, R., Bousquet, P., Collins, W., Cox, T., Eyring, V., Fowler, D., Fuzzi, S., Jöckel, P., Laj, P., Lohmann, U., Maione, M., Monks, P., Prevot, A. S. H., Raes, F., Richter, A., Rognerud, B., Schulz, M., Shindell, D., Stevenson, D. S., Storelvmo, T., Wang, W. C., van Weele, M., Wild, M., and Wuebbles, D.: Atmospheric composition change: Climate–Chemistry interactions, Atmos. Environ., 43, 5138–5192,, 2009. 

Jones, C. D., Collins, M., Cox, P. M., and Spall, S. A.: The Carbon Cycle Response to ENSO: A Coupled Climate–Carbon Cycle Model Study, J. Climate, 14, 4113–4129,<4113:tccrte>2.0.CO;2, 2001. 

Kawase, H., Nagashima, T., Sudo, K., and Nozawa, T.: Future changes in tropospheric ozone under Representative Concentration Pathways (RCPs), Geophys. Res. Lett., 38, L05801,, 2011. 

Kirschke, S., Bousquet, P., Ciais, P., Saunois, M., Canadell, J. G., Dlugokencky, E. J., Bergamaschi, P., Bergmann, D., Blake, D. R., Bruhwiler, L., Cameron-Smith, P., Castaldi, S., Chevallier, F., Feng, L., Fraser, A., Heimann, M., Hodson, E. L., Houweling, S., Josse, B., Fraser, P. J., Krummel, P. B., Lamarque, J.-F., Langenfelds, R. L., Le Quéré, C., Naik, V., O'Doherty, S., Palmer, P. I., Pison, I., Plummer, D., Poulter, B., Prinn, R. G., Rigby, M., Ringeval, B., Santini, M., Schmidt, M., Shindell, D. T., Simpson, I. J., Spahni, R., Steele, L. P., Strode, S. A., Sudo, K., Szopa, S., van der Werf, G. R., Voulgarakis, A., van Weele, M., Weiss, R. F., Williams, J. E., and Zeng, G.: Three decades of global methane sources and sinks, Nat. Geosci., 6, 813–823,, 2013. 

Koumoutsaris, S., Bey, I., Generoso, S., and Thouret, V.: Influence of El Niño–Southern Oscillation on the interannual variability of tropospheric ozone in the northern midlatitudes, J. Geophys. Res., 113, D19301,, 2008. 

Kozlov, V. S., Yausheva, E. P., Terpugova, S. A., Panchenko, M. V., Chernov, D. G., and Shmargunov, V. P.: Optical-microphysical properties of smoke haze from Siberian forest fires in summer 2012, Int. J. Remote Sens., 35, 5722–5741, 2014. 

Krol, M. and Lelieveld, J.: Can the variability in tropospheric OH be deduced from measurements of 1,1,1-trichloroethane (methyl chloroform)?, J. Geophys. Res., 108, 4125,, 2003. 

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,, 2010. 

Lathière, J., Hauglustaine, D. A., De Noblet-Ducoudré, N., Krinner, G., and Folberth, G. A.: Past and future changes in biogenic volatile organic compound emissions simulated with a global dynamic vegetation model, Geophys. Res. Lett., 32, L20818,, 2005. 

Lawrence, M., Jöckel, P., and Kuhlmann, R.: What does the global mean OH concentration tell us?, Atmos. Chem. Phys., 1, 37–49,, 2001. 

Manatsa, D. and Mukwada, G.: A connection from stratospheric ozone to El Niño-Southern Oscillation, Scientific Reports, 7, 5558,, 2017. 

Mann, G. W., Carslaw, K. S., Spracklen, D. V., Ridley, D. A., Manktelow, P. T., Chipperfield, M. P., Pickering, S. J., and Johnson, C. E.: Description and evaluation of GLOMAP-mode: a modal global aerosol microphysics model for the UKCA composition-climate model, Geosci. Model Dev., 3, 519–551,, 2010. 

Manning, M. R., Lowe, D. C., Moss, R. C., Bodeker, G. E., and Allan, W.: Short-term variations in the oxidizing power of the atmosphere, Nature, 436, 1001–1004,, 2005. 

Matthews, E. and Fung, I.: Methane emission from natural wetlands: Global distribution, area, and environmental characteristics of sources, Global Biogeochem. Cy., 1, 61–86,, 1987. 

McNorton, J., Chipperfield, M. P., Gloor, M., Wilson, C., Feng, W., Hayman, G. D., Rigby, M., Krummel, P. B., O'Doherty, S., Prinn, R. G., Weiss, R. F., Young, D., Dlugokencky, E., and Montzka, S. A.: Role of OH variability in the stalling of the global atmospheric CH4 growth rate from 1999 to 2006, Atmos. Chem. Phys., 16, 7943–7956,, 2016a. 

McNorton, J., Gloor, E., Wilson, C., Hayman, G. D., Gedney, N., Comyn-Platt, E., Marthews, T., Parker, R. J., Boesch, H., and Chipperfield, M. P.: Role of regional wetland emissions in atmospheric methane variability, Geophys. Res. Lett., 43, 11433–411444,, 2016b. 

McNorton, J., Wilson, C., Gloor, M., Parker, R. J., Boesch, H., Feng, W., Hossaini, R., and Chipperfield, M. P.: Attribution of recent increases in atmospheric methane through 3-D inverse modelling, Atmos. Chem. Phys., 18, 18149–18168,, 2018. 

McPhaden, M. J., Zebiak, S. E., and Glantz, M. H.: ENSO as an Integrating Concept in Earth Science, Science, 314, 1740–1745,, 2006. 

Meng, L., Paudel, R., Hess, P. G. M., and Mahowald, N. M.: Seasonal and interannual variability in wetland methane emissions simulated by CLM4Me' and CAM-chem and comparisons to observations of concentrations, Biogeosciences, 12, 4029–4049,, 2015. 

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,, 2015. 

Monks, S. A., Arnold, S. R., and Chipperfield, M. P.: Evidence for El Niño–Southern Oscillation (ENSO) influence on Arctic CO interannual variability through biomass burning emissions, Geophys. Res. Lett., 39, L14804,, 2012. 

Monks, S. A., Arnold, S. R., Hollaway, M. J., Pope, R. J., Wilson, C., Feng, W., Emmerson, K. M., Kerridge, B. J., Latter, B. L., Miles, G. M., Siddans, R., and Chipperfield, M. P.: The TOMCAT global chemical transport model v1.6: description of chemical mechanism and model evaluation, Geosci. Model Dev., 10, 3025–3057,, 2017. 

Murray, L. T., Mickley, L. J., Kaplan, J. O., Sofen, E. D., Pfeiffer, M., and Alexander, B.: Factors controlling variability in the oxidative capacity of the troposphere since the Last Glacial Maximum, Atmos. Chem. Phys., 14, 3589–3622,, 2014. 

Myhre, G., Shindell, D., Breì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, 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., Cambridge University Press, Cambridge, United Kingdom and New York, NY, USA, 659–740, 2013. 

Naik, V., Delire, C., and Wuebbles, D. J.: Sensitivity of global biogenic isoprenoid emissions to climate variability and atmospheric CO2, J. Geophys. Res., 109, D06301,, 2004. 

Naik, V., Voulgarakis, A., Fiore, A. M., Horowitz, L. W., Lamarque, J.-F., Lin, M., Prather, M. J., Young, P. J., Bergmann, D., Cameron-Smith, P. J., Cionni, I., Collins, W. J., Dalsøren, S. B., Doherty, R., Eyring, V., Faluvegi, G., Folberth, G. A., Josse, B., Lee, Y. H., MacKenzie, I. A., Nagashima, T., van Noije, T. P. C., Plummer, D. A., Righi, M., Rumbold, S. T., Skeie, R., Shindell, D. T., Stevenson, D. S., Strode, S., Sudo, K., Szopa, S., and Zeng, G.: Preindustrial to present-day changes in tropospheric hydroxyl radical and methane lifetime from the Atmospheric Chemistry and Climate Model Intercomparison Project (ACCMIP), Atmos. Chem. Phys., 13, 5277–5298,, 2013. 

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,, 2014. 

Nicely, J. M., Canty, T. P., Manyin, M., Oman, L. D., Salawitch, R. J., Steenrod, S. D., Strahan, S. E., and Strode, S. A.: Changes in Global Tropospheric OH Expected as a Result of Climate Change Over the Last Several Decades, J. Geophys. Res.-Atmos., 123, 10774–710795,, 2018. 

Nisbet, E. G., Dlugokencky, E. J., Manning, M. R., Lowry, D., Fisher, R. E., France, J. L., Michel, S. E., Miller, J. B., White, J. W. C., Vaughn, B., Bousquet, P., Pyle, J. A., Warwick, N. J., Cain, M., Brownlow, R., Zazzeri, G., Lanoisellé, M., Manning, A. C., Gloor, E., Worthy, D. E. J., Brunke, E.-G., Labuschagne, C., Wolff, E. W., and Ganesan, A. L.: Rising atmospheric methane: 2007–2014 growth and isotopic shift, Global Biogeochem. Cy., 30, 1356–1370,, 2016. 

Oman, L. D., Ziemke, J. R., Douglass, A. R., Waugh, D. W., Lang, C., Rodriguez, J. M., and Nielsen, J. E.: The response of tropical tropospheric ozone to ENSO, Geophys. Res. Lett., 38, L13706,, 2011. 

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., 118, 965–976,, 2013. 

Prinn, R. G., Huang, J., Weiss, R. F., Cunnold, D. M., Fraser, P. J., Simmonds, P. G., McCulloch, A., Harth, C., Salameh, P., O'Doherty, S., Wang, R. H. J., Porter, L., and Miller, B. R.: Evidence for Substantial Variations of Atmospheric Hydroxyl Radicals in the Past Two Decades, Science, 292, 1882–1888,, 2001. 

Randerson, J. T., van der Werf, G. R., Giglio, L., Collatz, G. J., and Kasibhatla, P. S.: Global Fire Emissions Database, Version 4.1 (GFEDv4), ORNL Distributed Active Archive Center,, 2017. 

Rap, A., Scott, C. E., Spracklen, D. V., Bellouin, N., Forster, P. M., Carslaw, K. S., Schmidt, A., and Mann, G.: Natural aerosol direct and indirect radiative effects, Geophys. Res. Lett., 40, 3297–3301,, 2013. 

Rap, A., Richards, N. A. D., Forster, P. M., Monks, S. A., Arnold, S. R., and Chipperfield, M. P.: Satellite constraint on the tropospheric ozone radiative effect, Geophys. Res. Lett., 42, 5074–5081,, 2015. 

Reddington, C. L., Morgan, W. T., Darbyshire, E., Brito, J., Coe, H., Artaxo, P., Marsham, J., and Spracklen, D. V.: Biomass burning aerosol over the Amazon: analysis of aircraft, surface and satellite observations using a global aerosol model, Atmos. Chem. Phys. Discuss.,, in review, 2018. 

Richards, N. A. D., Arnold, S. R., Chipperfield, M. P., Miles, G., Rap, A., Siddans, R., Monks, S. A., and Hollaway, M. J.: The Mediterranean summertime ozone maximum: global emission sensitivities and radiative impacts, Atmos. Chem. Phys., 13, 2331–2345,, 2013. 

Riese, M., Ploeger, F., Rap, A., Vogel, B., Konopka, P., Dameris, M., and Forster, P.: Impact of uncertainties in atmospheric mixing on simulated UTLS composition and related radiative effects, J. Geophys. Res., 117, D16305,, 2012. 

Rigby, M., Prinn, R. G., Fraser, P. J., Simmonds, P. G., Langenfelds, R. L., Huang, J., Cunnold, D. M., Steele, L. P., Krummel, P. B., Weiss, R. F., O'Doherty, S., Salameh, P. K., Wang, H. J., Harth, C. M., Mühle, J., and Porter, L. W.: Renewed growth of atmospheric methane, Geophys. Res. Lett., 35, L22805,, 2008. 

Rigby, M., Prinn, R. G., O'Doherty, S., Montzka, S. A., McCulloch, A., Harth, C. M., Mühle, J., Salameh, P. K., Weiss, R. F., Young, D., Simmonds, P. G., Hall, B. D., Dutton, G. S., Nance, D., Mondeel, D. J., Elkins, J. W., Krummel, P. B., Steele, L. P., and Fraser, P. J.: Re-evaluation of the lifetimes of the major CFCs and CH3CCl3 using atmospheric trends, Atmos. Chem. Phys., 13, 2691–2702,, 2013. 

Rossow, W. B. and Schiffer, R. A.: Advances in Understanding Clouds from ISCCP, B. Am. Meteorol. Soc., 80, 2261–2288,<2261:Aiucfi>2.0.Co;2, 1999. 

Schaefer, H., Smale, D., Nichol, S. E., Bromley, T. M., Brailsford, G. W., Martin, R. J., Moss, R., Englund Michel, S., and White, J. W. C.: Limited impact of El Niño–Southern Oscillation on variability and growth rate of atmospheric methane, Biogeosciences, 15, 6371–6386,, 2018. 

Schultz, M. G., Heil, A., Hoelzemann, J. J., Spessa, A., Thonicke, K., Goldammer, J. G., Held, A. C., Pereira, J. M. C., and van het Bolscher, M.: Global wildland fire emissions from 1960 to 2000, Global Biogeochem. Cy., 22, GB2002,, 2008. 

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,, 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–794,, 2007. 

Spivakovsky, C. M., Logan, J. A., Montzka, S. A., Balkanski, Y. J., Foreman-Fowler, M., Jones, D. B. A., Horowitz, L. W., Fusco, A. C., Brenninkmeijer, C. A. M., Prather, M. J., Wofsy, S. C., and McElroy, M. B.: Three-dimensional climatological distribution of tropospheric OH: Update and evaluation, J. Geophys. Res., 105, 8931–8980,, 2000. 

Spracklen, D. V., Jimenez, J. L., Carslaw, K. S., Worsnop, D. R., Evans, M. J., Mann, G. W., Zhang, Q., Canagaratna, M. R., Allan, J., Coe, H., McFiggans, G., Rap, A., and Forster, P.: Aerosol mass spectrometer constraint on the global secondary organic aerosol budget, Atmos. Chem. Phys., 11, 12109–12136,, 2011. 

Stevenson, D., Doherty, R., Sanderson, M., Johnson, C., Collins, B., and Derwent, D.: Impacts of climate change and variability on tropospheric ozone and its precursors, Faraday Discuss., 130, 41–57,, 2005. 

Stevenson, D. S., Johnson, C. E., Collins, W. J., Derwent, R. G., and Edwards, J. M.: Future estimates of tropospheric ozone radiative forcing and methane turnover – The impact of climate change, Geophys. Res. Lett., 27, 2073–2076,, 2000. 

Sudo, K. and Takahashi, M.: Simulation of tropospheric ozone changes during 1997–1998 El Niño: Meteorological impact on tropospheric photochemistry, Geophys. Res. Lett., 28, 4091–4094,, 2001. 

Szopa, S., Hauglustaine, D. A., and Ciais, P.: Relative contributions of biomass burning emissions and atmospheric transport to carbon monoxide interannual variability, Geophys. Res. Lett., 34, L18810,, 2007. 

Thompson, A. M., Witte, J. C., Hudson, R. D., Guo, H., Herman, J. R., and Fujiwara, M.: Tropical Tropospheric Ozone and Biomass Burning, Science, 291, 2128–2132,, 2001. 

Tilmes, S., Lamarque, J.-F., Emmons, L. K., Conley, A., Schultz, M. G., Saunois, M., Thouret, V., Thompson, A. M., Oltmans, S. J., Johnson, B., and Tarasick, D.: Technical Note: Ozonesonde climatology between 1995 and 2011: description, evaluation and applications, Atmos. Chem. Phys., 12, 7475–7497,, 2012. 

Turner, A. J., Fung, I., Naik, V., Horowitz, L. W., and Cohen, R. C.: Modulation of hydroxyl variability by ENSO in the absence of external forcing, P. Natl. Acad. Sci. USA, 115, 8931–8936,, 2018. 

van der Werf, G. R., Randerson, J. T., Collatz, G. J., Giglio, L., Kasibhatla, P. S., Arellano, A. F., Olsen, S. C., and Kasischke, E. S.: Continental-Scale Partitioning of Fire Emissions During the 1997 to 2001 El Niño/La Niña Period, Science, 303, 73–76,, 2004. 

van der Werf, G. R., Randerson, J. T., Giglio, L., van Leeuwen, T. T., Chen, Y., Rogers, B. M., Mu, M., van Marle, M. J. E., Morton, D. C., Collatz, G. J., Yokelson, R. J., and Kasibhatla, P. S.: Global fire emissions estimates during 1997–-2016, Earth Syst. Sci. Data, 9, 697–720,, 2017. 

Voulgarakis, A. and Field, R. D.: Fire Influences on Atmospheric Composition, Air Quality and Climate, Current Pollution Reports, 1, 70–81,, 2015. 

Voulgarakis, A., Naik, V., Lamarque, J.-F., Shindell, D. T., Young, P. J., Prather, M. J., Wild, O., Field, R. D., Bergmann, D., Cameron-Smith, P., Cionni, I., Collins, W. J., Dalsøren, S. B., Doherty, R. M., Eyring, V., Faluvegi, G., Folberth, G. A., Horowitz, L. W., Josse, B., MacKenzie, I. A., Nagashima, T., Plummer, D. A., Righi, M., Rumbold, S. T., Stevenson, D. S., Strode, S. A., Sudo, K., Szopa, S., and Zeng, G.: Analysis of present day and future OH and methane lifetime in the ACCMIP simulations, Atmos. Chem. Phys., 13, 2563–2587,, 2013. 

Voulgarakis, A., Marlier, M. E., Faluvegi, G., Shindell, D. T., Tsigaridis, K., and Mangeon, S.: Interannual variability of tropospheric trace gases and aerosols: The role of biomass burning emissions, J. Geophys. Res.-Atmos., 120, 7157–7173,, 2015. 

Wang, J. S., McElroy, M. B., Logan, J. A., Palmer, P. I., Chameides, W. L., Wang, Y., and Megretskaia, I. A.: A quantitative assessment of uncertainties affecting estimates of global mean OH derived from methyl chloroform observations, J. Geophys. Res., 113, D12302,, 2008. 

Wild, O.: Modelling the global tropospheric ozone budget: exploring the variability in current models, Atmos. Chem. Phys., 7, 2643–2660,, 2007. 

Wolter, K. and Timlin, M.: Monitoring ENSO in COADS with a seasonally adjusted principal component index, 17th Climate Diagnostics Workshop, 18-23 October 1992, Norman, OK, 52–57, 1993. 

Wolter, K. and Timlin, M. S.: Measuring the strength of ENSO events: How does 1997/98 rank?, Weather, 53, 315–324,, 1998. 

Wolter, K. and Timlin, M. S.: El Niño/Southern Oscillation behaviour since 1871 as diagnosed in an extended multivariate ENSO index (MEI.ext), Int. J. Climatol., 31, 1074–1087,, 2011. 

Yan, X., Akiyama, H., Yagi, K., and Akimoto, H.: Global estimations of the inventory and mitigation potential of methane emissions from rice cultivation conducted using the 2006 Intergovernmental Panel on Climate Change Guidelines, Global Biogeochem. Cy., 23, GB2002,, 2009.  

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,, 2013. 

Zeng, G., and Pyle, J. A.: Influence of El Niño Southern Oscillation on stratosphere/troposphere exchange and the global tropospheric ozone budget, Geophys. Res. Lett., 32, L01814,, 2005. 

Zhang, J., Tian, W., Wang, Z., Xie, F., and Wang, F.: The Influence of ENSO on Northern Midlatitude Ozone during the Winter to Spring Transition, J. Climate, 28, 4774–4793,, 2015. 

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,, 2016. 

Zhu, Q., Peng, C., Ciais, P., Jiang, H., Liu, J., Bousquet, P., Li, S., Chang, J., Fang, X., Zhou, X., Chen, H., Liu, S., Lin, G., Gong, P., Wang, M., Wang, H., Xiang, W., and Chen, J.: Interannual variation in methane emissions from tropical wetlands triggered by repeated El Niño Southern Oscillation, Glob. Change Biol., 23, 4706–4716,, 2017. 

Ziemke, J. R., Chandra, S., Oman, L. D., and Bhartia, P. K.: A new ENSO index derived from satellite measurements of column ozone, Atmos. Chem. Phys., 10, 3711–3721,, 2010. 

Short summary
Wildfires and meteorology have a substantial effect on atmospheric concentrations of greenhouse gases such as methane and ozone. During the 1997 El Niño event, unusually large fire emissions indirectly increased global methane through carbon monoxide emission, which decreased the oxidation capacity of the atmosphere. There were also large regional changes to tropospheric ozone concentrations, but contrasting effects of fire and meteorology resulted in a small change to global radiative forcing.
Final-revised paper