Limited impact of El Niño – Southern Oscillation on variability and growth rate of atmospheric methane

The El Niño–Southern Oscillation (ENSO) has been suggested as a strong forcing in the methane cycle and as a driver of recent trends in global atmospheric methane mole fractions [CH4]. Such a sensitivity of the global CH4 budget to climate events would have important repercussions for climate change mitigation strategies and the accuracy of projections for future greenhouse forcing. Here, we test the impact of ENSO on atmospheric CH4 in a correlation analysis. We use local and global records of [CH4], as well as stable carbon isotopic records of atmospheric CH4 (δCH4), which are particularly sensitive to the combined ENSO effects on CH4 production from wetlands and biomass burning. We use a variety of nominal, smoothed, and detrended time series including growth rate records. We find that at most 36 % of the variability in [CH4] and δCH4 is attributable to ENSO, but only for detrended records in the southern tropics. Trend-bearing records from the southern tropics, as well as all studied hemispheric and global records, show a minor impact of ENSO, i.e. < 24 % of variability explained. Additional analyses using hydrogen cyanide (HCN) records show a detectable ENSO influence on biomass burning (up to 51 %– 55 %), suggesting that it is wetland CH4 production that responds less to ENSO than previously suggested. Dynamics of the removal by hydroxyl likely counteract the variation in emissions, but the expected isotope signal is not evident. It is possible that other processes obscure the ENSO signal, which itself indicates a minor influence of the latter on global CH4 emissions. Trends like the recent rise in atmospheric [CH4] can therefore not be attributed to ENSO. This leaves anthropogenic methane sources as the likely driver, which must be mitigated to reduce anthropogenic climate change.


Introduction
Attributing recent changes in the methane budget, and the associated impact on its growth rate, to specific natural or anthropogenic causes is essential for climate change mitigation.The impact of climatic variability on methane emissions is particularly important to assess the potential for CH 4 release under future climate scenarios (e.g., from permafrost and wetland environments, as well as gas hydrates) in a reinforcing feedback.Atmospheric methane mole fractions [CH 4 ] have increased by 140 % over pre-industrial levels (MacFarling Meure et al., 2006).The associated increase in radiative forcing makes CH 4 the second most important anthropogenic greenhouse gas (Shindell et al., 2009).The long-term [CH 4 ] increase until the late 1990s can be attributed to increasing emissions from fossil fuel production (Ferretti et al., 2005;Schaefer et al., 2016), as well as sources from agriculture (enteric fermentation in livestock, rice production), waste management, and anthropogenic burning (van Aardenne et al., 2001;Saunois et al., 2016).After a plateau in the early 2000s, [CH 4 ] has been rising again since 2007.Considering recent reconstructions of methane's dominant atmospheric sink, i.e. the hydroxyl radical OH, we consider it likely that increasing emissions contribute to (Rigby et al., 2017), if not dominate (Naus et al., 2018), the [CH 4 ] rise.If so, the methane source type that varied can be investigated with measurements of stable carbon isotope ratios in atmospheric methane (δ 13 CH 4 ).The latter are influenced by the relative source contributions from 13 C-depleted biogenic, 13 C-rich pyrogenic, and thermogenic methane with intermediate δ 13 C. Isotope studies suggest that biogenic methane sources make either a dominant (Schaefer et al., 2016;Nisbet et al., 2016) Figure 1.Regions of ENSO impacts and monitoring stations used in this study.The map indicates the locations of the atmospheric monitoring stations on Ascension Island (ASC), Samoa (SMO), Baring Head (BHD), and Lauder (LAU).General precipitation anomalies during Northern Hemisphere El Niño conditions for December-February are taken from https://www.climate.gov/news-features/featured-images/global-impacts-el-nino-and-la-nina (last access: 5 December 2017).El Niño dry regions in June-August are similar for southern Asia and South America; during La Niña events opposite patterns for wet-and dryness develop in roughly the same regions.
or strong (Worden et al., 2017) contribution to the recent [CH 4 ] rise.Biogenic methane comes predominantly from wetlands and agriculture.Schaefer et al. (2016) suggested agriculture as the more likely cause, primarily because satellite data place the increased emissions in Southeast Asia, India, and China (Houweling et al., 2014).However, this geographic footprint from an inversion of satellite data is also consistent with fluxes from one particular wetland emissions model (Houweling et al., 2014).Other studies also assume a stronger role of wetlands due to drier conditions during the plateau years (Bousquet et al., 2006) and higher wetland emissions afterwards, which are attributed to a switch to predominant La Niña conditions around 2007 (Bousquet et al., 2011;Nisbet et al., 2016).La Niña is the cold phase of El Niño-Southern Oscillation (ENSO) cycles, which have a strong impact on precipitation anomalies in tropical regions (Ropelewski and Halpert, 1987;Lyon and Barnston, 2005) (Fig. 1) that are key source areas for methane production from wetlands and biomass burning (Kirschke et al., 2013).ENSO impacts are strongest in the tropics, generally from December to February.During El Niño (La Niña) events in the December to February period, it tends to be drier (wetter) in the Indonesian region, north-east Brazil, and south-eastern Africa, whereas it tends to be wetter (drier) in the southern USA and Mexico, eastern China and Taiwan, and east-central Africa (Fig. 1).During El Niño (La Niña) events in the June to August period, it tends to be drier (wetter) in the Indonesian region, central America, and India.
The generally drier conditions during El Niños suppress global wetland emissions in models by up to 19 Tg yr −1 in the 1990s (Hodson et al., 2011).Several anthropogenic sources are subject to the same ENSO forcing and are expected to vary in concert with wetlands (e.g., rice agriculture, possibly livestock).At the same time, dry El Niño phases enhance CH 4 emissions from both natural and anthropogenic biomass burning (van der Werf et al., 2006).Wet La Niña conditions have the opposite effect; summed across the globe they increase wetland emissions and lower biomass burning CH 4 .As tropical wetland fluxes are considerably larger than biomass burning emissions (Saunois et al., 2017), the expected net effect is a lower [CH 4 ] growth rate caused by El Niño conditions and a higher one due to La Niñas.The ENSO impact on δ 13 CH 4 should be more pronounced than the one on [CH 4 ] because changes in wetland and biomass burning emissions combine to enrich atmospheric CH 4 in 13 C during El Niños and deplete it during La Niñas.Biogenic methanogenesis in wetlands discriminates strongly against 13 C and creates methane that is 13 C-depleted (δ 13 C = −58 ‰ for tropical wetlands) relative to the plant precursor material (δ 13 C of −12 ‰ to −28 ‰ ) and to the combined total of global emissions (δ 13 C ∼ −53.5 ‰).In contrast, during burning the isotope ratios of the precursor plant material are essentially conserved and lead to δ 13 C ∼ −22 ‰ for CH 4 emissions from fires (Schwietzke et al., 2016).The simultaneous suppression of 13 C-depleted wetland CH 4 and enhancement of very 13 C-rich pyrogenic emissions (and vice versa) act in the same direction on the δ 13 CH 4 of the combined source.The latter should be detectable in atmospheric δ 13 CH 4 records if the impact of ENSO on the CH 4 cycle is sufficiently large, as is predicted by the emission anomalies in wetland emission models (Hodson et al., 2011), reconstructed from satellite observations of burned area (van der Werf et al., 2010), and observed through variability in hydrogen cyanide (HCN) (Pumphrey et al., 2018), which is an indicator of biomass burning.
Varying contributions from wetlands dominated by C 3 and C 4 plants, which differ in the δ 13 CH 4 of their emissions, may be part of the ENSO-CH 4 signal or work to obscure it if controlled by other drivers.In general, we assume that δ 13 CH 4 of the various emission sources has not changed over the ∼ 35year period of our study.Although such changes, correlated to atmospheric CO 2 mole fractions, have been reported to occur over centuries to millennia in ice core studies (Möller et al., 2013), they are likely negligible over the short duration and > 20 % CO 2 change of our study period.
Changes in OH have also been suggested as partial or dominant drivers in recent CH 4 trends, both for the onset of the 1999-2006plateau (McNorton et al., 2016;;Schaefer et al., 2016) and for the post-2007 [CH 4 ] increase (Rigby et al., 2017;Turner et al., 2017).A chemistry-climate model suggests that ENSO modulates tropical OH (where hydroxyl levels are highest) via changes in NO x production through lightning, ozone availability, and specific humidity, as well as emissions of reactive carbon (Turner et al., 2018).Resulting changes in methane removal could create their own signal in atmospheric records of [CH 4 ] and δ 13 CH 4 .They could also either reinforce or dampen the emission impacts discussed above.
We conduct correlation analyses between ENSO variability and [CH 4 ], as well as δ 13 CH 4 records, to quantify how much ENSO anomalies in emissions and sinks affect atmospheric CH 4 .Specifically, we explore how much of the yearto-year variability in atmospheric methane can be attributed to ENSO and how large the ENSO-CH 4 signal is in dependence of latitude.We test if recent trends in methane growth rate can be attributed to wetland emissions controlled by ENSO dynamics or if agricultural sources are more likely drivers.ENSO is quantified by four different indices, which are based on ocean temperature, sea level pressure gradients, and a multivariate combination.[CH 4 ] and δ 13 CH 4 time series from four different locations were used, two from stations in the southern tropics (Samoa, SMO, and Ascension Island, ASC), the southern mid-latitudes (Baring Head, NZ; BHD), taken as representative of the Southern Hemisphere, and global average time series of [CH 4 ] and δ 13 CH 4 calculated from a network of global stations (Dlugokencky et al., 2011;Schaefer et al., 2016).We also investigate ENSO's impact on HCN data measured in Lauder, NZ (LAU), to quantify the biomass burning contribution separately.The aim is to detect the impact of ENSO on atmospheric CH 4 on various spatial scales.

Data
For access to all data sets used in this study, see Sect.10.

ENSO indices
We used four different indices in our analysis to cover various climatic effects of the ENSO cycle (Figs.2a and 3a).The Southern Oscillation Index (SOI) is calculated from the gradient in mean sea-level pressure observations at Tahiti and Darwin, Australia (Troup, 1965).Further information on the SOI is given by Horel and Wallace (1981) and Trenberth (1976).The Ocean Niño Index (ONI) uses sea surface temperature (SST) anomalies in the eastern Pacific Niño 3.4.region (5 • N-5 • S, 120-170 • W), which show smaller intraseasonal variability than pressure and are further smoothed by using 3-month running means (Barnston et al., 1997;Kousky and Higgins, 2007).
The El Niño Modoki Index (EMI) is based on SST anomalies in the central Pacific (Ashok et al., 2007) rather than the eastern Pacific (the canonical El Niño).Events with the largest SST anomalies in the Modoki region show differences in the climate teleconnections to canonical El Niño events.The tropical precipitation differences are modest, but large differences in tropospheric circulation and wind anomalies (Yeh et al., 2009) can produce large extratropical differences in precipitation and temperature.The EMI has also been shown to be a significant predictor of tropical atmospheric ozone variations (Xie et al., 2014).
Variability in both atmospheric pressure and SST anomalies informs the Multivariate ENSO Index (MEI) (Wolter andTimlin, 1993, 1998).The various indices correlate highly with each other (r 2 = 0.85 and higher), except the EMI (r 2 between 0.33 and 0.52 for correlations with SOI, ONI, and MEI), which deviates from the others during the strong 1997-1998 El Niño event.Excluding the latter brings the correlation to r 2 between 0.74 and 0.79.
An ENSO index based on precipitation data, the ESPI, (Curtis and Adler, 2000) correlates very highly with the MEI, the ONI, and the SOI (r 2 of 0.902, 0.909, and 0.839, respectively).Therefore, we did not conduct separate calculations for the ESPI.(Lowe et al., 1991).Both data sets are on the same international scale (Dlugokencky et al., 2005), although for the presented analysis, internal consistency of the time series is the relevant criterion; interlaboratory offsets do not affect the findings.The individual time series (Fig. 2b-e) show seasonal cycles, inter-annual variability (IAV), and long-term trends.To investigate ENSO effects on these different timescales we derived the following seven records from the individual measurements at each sta-tion (Table 1): first, the nominal monthly mean values to capture the full variability in the data ("nom"), second, 12-month running means to represent IAV and trends ("run"), third, monthly resolved growth rate defined as the difference between the following 12 months and the preceding 12 months ("gro"), and fourth, a residual ("res") as calculated by seasonal trend analysis by loess (STL; Cleveland et al., 1990).The seasonal window was set at 120 months, which forces a uniform seasonal cycle over the duration of the record.
The residual therefore represents IAV in the expression of the seasonal cycle as well as other short-term anomalies.The fifth, sixth, and seventh records are the detrended time series for which the STL trend component is subtracted from the monthly means with subsequent determination of detrended monthly means, 12-month running means, and growth rate ("det-nom", "det-run", "det-gro").

δ 13 CH 4 time series
The δ 13 CH 4 time series used in this study were measured at three different laboratories, i.e., the Institute of Arctic and Alpine Research (INSTAAR), USA; the Institute for Environmental Physics (IUP) at Heidelberg University, Germany; In this analysis, we use the measurements covering 1992-2016 (Fig. 3c).For the global δ 13 CH 4 data set we conducted the analysis for the nomi-nal annual means (nom) and growth rate, i.e. the difference between two subsequent yearly values (gro).We also detrended the time series by subtracting linear trends for the sub-periods 1992-1999 and 2007-2016 ("det") and then calculating a detrended growth rate (det-gro).For the singlestation δ 13 CH 4 records of BHD, ASC, and SMO, we derived the same seven records as described for the [CH 4 ] data (Fig. 3d-f).

HCN time series
HCN retrievals were computed from mid-infrared solar spectra measured at LAU (45.04 • S, 169.68 • E) as part of the Network for the Detection of Atmospheric Composition Change (NDACC).The time series has been described by Zeng et al. (2012), but the data used here are from updated retrievals using the improved SFIT4 algorithm (NDACC, 2014).The HCN data show strong seasonality that is even more pronounced in the updated retrievals.Zeng et al. (2012) found a significant negative trend for 1997-2009 and attributed it to variations in biomass burning.A similar deseasonalized trend is apparent in the updated record.HCN values are reported here as petamolecules cm −2 .Measurements cover the period 1998-2017 when combined for two different instruments with a change-over point in 2000.We conducted our analyses for total column values (0-100 km).The latter signal is dominated by the tropospheric burden as measured in the 0-12 km height partial column; the correlation between total and tropospheric HCN is r 2 = 0.997.In addition, we investigated whether the stratospheric HCN burden is differently impacted by ENSO.To that end, we used the 12-100 km partial column, which holds ∼ 22 % of the total HCN burden.This layer shows lower correlation with the total column record (r 2 = 0.45).Analogous to the monthly resolved methane records, we constructed monthly means, 12-month running means, growth rates, and STL residuals for the total column and stratospheric HCN data (Table 1, Fig. 3b).No detrended records other than STL residual were considered.

Analysis
We conducted correlation analyses between the time series of a chosen ENSO index and either a [CH 4 ], δ 13 CH 4 , or HCN record as the dependent variable.The degree of correlation is quantified by the square (r 2 value) of the Pearson correlation coefficient or, alternatively, of the Spearman ranking coefficient.The Pearson coefficient is more commonly used, but it assumes linear relationships between the variables and may underestimate non-linear correlations.We therefore also used the Spearman rank, which does not require linearity.Note that not all correlation combinations were tested using both coefficients.
A lag time between ENSO forcing and detection of resulting δ 13 CH 4 or HCN variability at the measurement site, (or in the global average) is likely, due to a variety of factors that may lead to lags of unknown length and some of which may be cumulative, e.g., hydrology, plant growth and decay, microbial response, and seasonal triggers for methanogenesis or burning, as well as atmospheric chemistry, mixing, and transport between source regions and sampling sites.Therefore, it is difficult to define a cut-off for lags.Literature estimates of specific lags range from days (Chamberlain et al., 2016) to 7 months (Zhang et al., 2018;Zhu et al., 2017), not counting atmospheric transport.Given ENSO variability with a periodicity of 2-7 years (McPhaden et al., 2006), our analysis therefore allows for lag times of up to 5 years in monthly increments in the calculations and reports the maximum r 2 and lag time (in months) for a given ENSO-[CH 4 ] / δ 13 CH 4 / HCN combination.We conducted the analysis for all permutations of the four ENSO indices as monthly means and their 12-month running means as well as the [CH 4 ], δ 13 CH 4 , and HCN data products listed in Sect.2.1 and Table 1.For all [CH 4 ], δ 13 CH 4 , and HCN parameters we used the period 1998-2016, except for ASC, where data are only available from late 2000.Using the same period for all time series avoids differing correlation results due to varying data coverage.The period includes the strong El Niños of 1998 and 2015, as well as the strong La Niñas of 1999, 2007, and 2010.We also calculated correlations for the period 1983-2016 ([CH 4 ([CH 4 ] of SMO, ASC, and global) and 1992-2016 (δ 13 CH 4 at BHD and global).

Results
Most combinations have r 2 values < 0.1 when comparing one dependent data set to the different ENSO time series (Tables 2-4).In the following, we only summarize results for the highest r 2 for each dependent time series (across all the nominal, smoothed, and detrended records for a station).Given that Pearson coefficient and Spearman rank give comparable results (Tables 3 and 4), we quote the Spearman results, unless otherwise mentioned.The p values for the Spearman ranks indicate that all results for r 2 > 0.1 are significant (p< 0.001), except for global δ 13 CH 4 correlations, for which no p values below 0.05 occur.Although the analysis provides r 2 values for lags up to 60 months (Tables 2-4), we consider it likely that lags of > 3 years indicate spurious correlations, given that individual ENSO events last 1-2 years and global atmospheric mixing times are on the order of 1 year.Therefore, we also report the highest r 2 for lags < 3 years in the following sections.For other cases with lags > 3 years in Tables 2-4, the highest relevant r 2 value is lower than the reported value, which places an upper limit on the influence of ENSO.
Methane mole fractions show correlations with ENSO of r 2 values up to 0.36 at SMO, but only for detrended time series (Table 2).The highest values are from (detrended) growth rates, which can be more indicative of dynamics within an ENSO event, rather than its overall emissions impact (Zhang et al., 2018).For SMO detrended [CH 4 ] series, lag times are fairly consistent across the various ENSO indices and generally shorter than 1 year.For other [CH 4 ] records at SMO and ASC, the highest correlations are r 2 < 0.23 and have lags of over 3 years (r 2 < 0.19 for lags < 3 years).The global running mean [CH 4 ] time series shows r 2 = 0.24 (lag: 4.5 years; r 2 = 0.04 for lag < 3 years) with the SOI running mean for the period 1998-2016.However, for the full length of available data, as well as all BHD records, all correlations are below r 2 = 0.20, with lag times that are variable, extremely short (zero or 1 month), or over 3 years.
The highest correlations are between HCN total column running means, as well as stratospheric growth rates, and 12-month running mean ENSO records (Table 3).Here, ENSO accounts for 30 %-51 % of the observed variability, depending on the ENSO index.For both total and stratospheric HCN, lag times for maximum correlation are generally shorter than 1 year and are consistent (≤ 6 months difference) between the various ENSO indices, with exception of the EMI.
The δ 13 CH 4 records from the stations SMO, ASC, and BHD all have r 2 values below 0.24 (Table 3).Variability in lag times between different ENSO indices for the same dependent record is generally high.
None of the global δ 13 CH 4 series produced statistically robust correlations with ENSO; all p values were higher than 0.05.The following findings are therefore not relevant for further interpretation.The highest correlation is between global detrended δ 13 CH 4 and SOI monthly means with r 2 = 0.37.Global δ 13 CH 4 , is the only parameter for which ENSO monthly means produce higher correlations than the smoothed (12-month running mean) record.Because the correlation calculation between annual δ 13 CH 4 and ENSO monthly means is specific to the month of year, this indicates that global δ 13 CH 4 is more sensitive to the seasonality of ENSO than its IAV.The actual ENSO influence on global δ 13 CH 4 is shown through correlation with running ENSO indices, which is highest between nominal δ 13 CH 4 values and SOI, with Pearson r 2 = 0.25 for 1998-2016.For the period 1992-2016 this value drops to Pearson r 2 = 0.20.The lack of statistical robustness for global δ 13 CH 4 -ENSO correlations may stem from the different resolution of the two sets of time series.In this case, the southern hemispheric record from BHD may represent the extratropical impact of ENSO on δ 13 CH 4 .
The full BHD record for 1992-2016 gives very similar results as the 1998-2016 subset used for comparison with the other stations (as discussed above).However, the shorter subset for 1998-2014 produces larger Pearson r 2 values (0.26    (Nisbet et al., 2016) this should be manifested in sizeable correlations.

Contrasting correlation patterns for [CH 4 ] and δ 13 CH 4 versus HCN
In all [CH 4 ] and δ 13 CH 4 records, ENSO cycles explain at most about one-third of the variability in detrended records and less than one-quarter in others.This is true even for the southern tropics, where ENSO has strong climatic impacts and where the majority of low-latitude wetland emissions and of biomass burning emissions originate (Kirschke et al., 2013).Correlations found for ASC and SMO, which represent this latitude band in our study, exceed those for the southern mid-latitudes or the global record only by a limited margin and only for detrended records.Further, inconsistent lag times, lags of more than 3 years, and higher correlation coefficients for the exclusion of major ENSO events point to spurious correlations.
In contrast, we find a prominent influence of ENSO on the biomass burning proxy HCN.ENSO impacts on HCN have been reported before, e.g., by Pumphrey et al. (2018), who observe suppression of HCN levels during La Niña events and enhancement during El Niños, particularly in equatorial Asia.That study found a rather confined geographical impact of El Niño events with strongly enhanced HCN emissions around Malaysia, Indonesia, and Papua New Guinea, as well as generally rapid transport eastward and to the stratosphere.We speculate that the fast, upward transport (although not observed for all El Niño events) explains why stratospheric growth rates are the most sensitive data set to ENSO.For the total column, the HCN burden is concentrated in lower tropospheric levels and may be subjected to more mixing of different air parcels.According to the results of Pumphrey et al. (2018), data from LAU in the southern mid-latitudes are outside the region of the strongest HCN signal.This is also evident in the zonal mean HCN climatologies of Sheese et al. (2017).However, ENSO accounts for up to 51 % of the variability in our biomass burning proxy record.One explanation for the lower combined wetland-pyrogenic δ 13 CH 4 signal is low sensitivity of wetland CH 4 production to ENSO events.This is consistent with r 2 values of 0.12-0.26 between modelled wetland methane emissions (using different climate data sets as drivers) and MEI as reported by Zhang et al. (2018).Alternatively, other processes in the CH 4 cycle obscure the ENSO impacts.

Impact of ENSO on methane emission rates
In a correlation analysis by Zhu et al. (2017), ENSO explained 49 % of IAV in modelled tropical wetland CH 4 emissions.This is far higher than the combined effect with biomass burning on δ 13 CH 4 in this study, and therefore seems to be an overestimate.Even so, the magnitude of the modelled emission changes is 6 Tg yr −1 at most.The modelling study of Hodson et al. (2011) finds slightly larger anomalies in global wetland emissions due to ENSO with mean reductions of −9 ± 3 Tg yr −1 and mean gains of +8 ± 4 Tg yr −1 for El Niño and La Niña events, respectively.Pandey et al. (2017) found in a comprehensive inversion study that the net effect of the strong 2011 La Niña on tropical and northern extratropical CH 4 emissions was a global increase of +6.6 Tg yr −1 .The wetland emission anomalies are expected to be partly compensated by changes in biomass burning that are of opposite sign.We are not aware of studies that quantify biomass burning anomalies for specific ENSO events.Assuming that ENSO is the main control of biomass burning emissions of CH 4 , the IAV in the Global Fire Emissions Database (GFED) data (van der Werf et al., 2010) may serve as an indication for possible ENSO impacts.In that case, the standard deviation of 2.4 Tg yr −1 for 1998-2014 would approximate the average impact, with maximum anomalies of up to 4 Tg yr −1 .We use these numbers together with results from Hodson et al. (2011) in the following proof-of-concept discussions.The combined wetland-biomass burning anomalies are ∼ 6 Tg yr −1 for average ENSO events and ∼ 8 Tg yr −1 for extreme ones, restricted to 1-2-year-long individual events.This is well short of the sustained increase after 2007 when yearly emissions were ∼ 20 Tg higher than during the 1999-2006 plateau period and the 9 Tg yr −1 reduction during the 1990s (Schaefer et al., 2016).Previous findings that modelled tropical (Zhu et al., 2015) and global (Zhang et al., 2018) wetland CH 4 emissions can explain at most 25 % and 14 %, respectively, of the variation in atmospheric methane growth rates therefore agree with our results that ENSO exerts only a minor control on atmospheric CH 4 .

Process-based understanding of ENSO impact on wetlands
A major contribution of ENSO to the recent [CH 4 ] increase is inconsistent with independent assessments of wetland response, as shown above, but our findings do not detect any clear minor contribution of ENSO to [CH 4 ] and δ 13 CH 4 time series, either.Several reasons may explain the lack of correlation, whereby we assume that wetlands respond less than proposed.The main ENSO forcing on tropical wetland CH 4 production is thought to be via wetland extent, which is driven by precipitation (Hodson et al., 2011, and Holmes et al., 2015, in contrast to Zhu et al., 2017, who find temperature to be dominant).However, a case study in the eastern Amazon finds that precipitation changes explain only 21 % of wetland CH 4 emission variance during the wet season and 7 % over the whole year (Basso et al., 2016).The lack of a direct link between precipitation and wetland CH 4 production is also evident in the large range in output from various wetland models even when forced with the same meteorological conditions (Melton et al., 2013), although the disagree-ment between models could also be due to an incomplete understanding of influences on the wetland cycle other than precipitation (Turetsky et al., 2014;Bridgham et al., 2013;Parker et al., 2018).Zhang et al. (2018) report an evolving response of wetland emissions to El Niños, where an initial reduction due to decreased wetland extend is counteracted by increased microbial activity under higher temperatures during the later stages of the event.A complex response of wetland CH 4 production is not only seen in models, however.The inversion study of Pandey et al. (2017) found a global increase of +6.6 Tg yr −1 for the strong 2011 La Niña, but a reduction by −6.1 Tg yr −1 during the 2012 weak La Niña.Similarly, Liu et al. (2017) found that El Niño conditions produced opposing weather forcing and carbon cycle responses between various tropical regions, as well as differing ones between the 1998 and 2015 events.Another example of this is flooding in the Amazon region during La Niña events, while flooding in the wetlands of the Paraná region occurs during El Niños (Parker et al., 2018).Depending on the strength and geographical expression of the climate anomaly, ENSO may thus cause regional or global emission anomalies that are opposite to the expected pattern.

Evaluating the consistency of ENSO impacts throughout the record
The atmospheric [CH 4 ] history shows global emission reductions in the 1990s and increases after 2007 (Schaefer et al., 2016).This would be consistent with ENSO forcing of the methane cycle, whereby the 1990s were dominated by drier El Niño periods, whereas the recent years of predominant La Niña conditions were wetter.Given that the magnitude of the low-latitude wetland CH 4 source exceeds pyrogenic emissions rates, the expected emissions history would qualitatively match atmospheric trends.Also, for a short period between 2008 and 2011, Schaefer et al. (2016) observed the activation of CH 4 emissions with an extremely 13 C-depleted cumulative δ 13 CH 4 (∼ −75 ‰).Such a value on the global scale is hard to match by a single source type.The cumulative effect of wetland enhancement and fire suppression forced by the 2008 La Niña event would provide an excellent explanation.However, the isotopic signal of the emissions reductions in the 1990s should be similar if ENSO forcing was the cause.In contrast, Schaefer et al. (2016) found that the "lost emissions" during that period are quite 13 Crich and rather indicate a reduction in fossil fuel methane.
An alternative interpretation of these isotope trends by Rice et al. (2016) requires simultaneous reductions of pyrogenic and biogenic emissions, which is also inconsistent with the expected ENSO forcing.A consistent match between ENSO phases and global δ 13 CH 4 is therefore neither evident in the dominant δ 13 CH 4 trends nor in the correlation analysis presented in this study.

Using isotopes to attribute emission changes
The impact of an ENSO emissions "perturbation" (i.e. the combined emissions anomaly of an event) on atmospheric δ 13 CH 4 can be assessed in isotope mass balance calculations according to where, for a given source, S and δ represent emission rate and δ 13 CH 4 , respectively (note that for scenarios discussed here, S may be negative, i.e., a reduction in emissions).Using generic isotope source signatures for biogenic, fossil fuel, and pyrogenic methane emissions from Schwietzke et al. (2016), we find that the average La Niña perturbations proposed in Sect.4.3.have an effective δ 13 CH 4 of −79 ‰, with −83 ‰ for extreme ones.As expected, the combined isotope leverage of wetland enhancement and fire reductions on the global source is strong, equalling the leverage of much larger source anomalies (20 Tg yr −1 ), with lower δ 13 CH 4 of ∼ −60 ‰ after 2007 as calculated by Schaefer et al. (2016).
In addition to the assumed 6 Tg yr −1 ENSO perturbation, another ∼ 14 Tg yr −1 of emissions with δ 13 CH 4 = −52 ‰ would be necessary to produce the observed [CH 4 ] and δ 13 CH 4 trends.The isotope mass balance then shows that the non-ENSO additional emissions are a roughly equal mix of fossil fuel and biogenic methane.Noting that the assumption that all years after 2007 experienced average La Niña conditions is unrealistic; these findings therefore show the following three points: (i) ENSO effects alone cannot explain the post-2007 [CH 4 ] rise.(ii) There was an increase in biogenic sources in addition to ENSO-driven wetland anomalies.Other wetland variability may have contributed to the rise (Zhang et al., 2018); given the range in wetland model output (Melton et al., 2013) this stands to be confirmed by ensemble runs.In the absence of boreal emission increases (Sweeney et al., 2016), the only biogenic source large enough to accommodate the required changes is agriculture (Saunois et al., 2016).(iii) Any ENSO-driven reduction in biomass burning after 2007 allows for, or requires, growing fossil fuel emissions.The latter has recently been proposed by Worden et al. (2017), who reconstructed larger biomass burning reductions after 2007 than recorded by GFED, although without assigning the reductions to ENSO or other causes.

Role of other methane cycle processes
There is an alternative explanation for the lack of correlation between ENSO and the methane records.ENSO could affect CH 4 emissions from tropical wetlands and biomass burning as predicted by Hodson et al. (2011)  ] and δ 13 CH 4 compared to southern mid-latitudes and global average would be expected for several of these scenarios.This is because both biomass burning and wetland emissions show strong maxima in the southern tropics and should be the dominant sources in this latitudinal band (Kirschke et al., 2013).The detrended [CH 4 ] records from SMO show such a signal, but one that explains only one-third of the IAV and does not seem to have significant impact on the trends.Further, we do not find higher ENSO forcing of the δ 13 CH 4 variability, even in the core region of its climatic impact.Corbett et al. (2017)  Higher emissions of methane and carbon monoxide from biomass burning will draw down OH and weaken the sink.Emission factors from fires for CO are between 10-and 30fold higher than for CH 4 (van der Werf et al., 2017), so that the biomass burning dynamics dominate the source of reactive carbon, leaving less OH during El Niños and more during La Niñas to draw down CH 4 .This would provide a negative feedback for the emissions' [CH 4 ] signal from ENSO forcing.In contrast, the feedback on the ENSO emissions' δ 13 CH 4 signal would be positive due to varying enrichment of 13 C-methane through sink fractionation (less removal leads to less 13 C enrichment of relatively 13 C-depleted wetland emissions during La Niñas; more removal increases the 13 C enrichment from biomass burning emissions during El Niños further).In addition to the reactive carbon effect, Turner et al. (2018) found a further OH increase during La Niñas due to higher lightning rates with NO x production.Turner et al. (2018) could attribute 17 % of OH variability that is forced by climate cycles (rather than emissions of other atmospheric compounds) to ENSO.This is a minor part of the variability, but in consequence the dampening effect on [CH 4 ] and the reinforcing feedback on δ 13 CH 4 would add further to the reactive carbon feedbacks.In our correlation results these sink impacts are not apparent, as the [CH 4 ] correlations for the tropical stations are higher than δ 13 CH 4 correlations (Tables 2 and 3).Nevertheless, the OH dynamics provide a possible explanation for the limited ENSO impact on [CH 4 ] variability and trends.They also make δ 13 CH 4 a conservative proxy for the influence that ENSO exerts on tropical methane.
Whether ENSO has less influence on CH 4 emissions than assumed or whether such an impact is overwhelmed by atmospheric removal or other CH 4 cycle processes, our results suggest that global atmospheric trends in [CH 4 ] and δ 13 CH 4 are dominated by other components in the methane budget.

Conclusions
To study the impact of natural climate variability on recent trends in atmospheric methane concentration, we investigated the correlation between ENSO cycles and records of the mole fractions and stable carbon isotopes of methane, as well as HCN as a biomass burning indicator.As δ 13 CH 4 is subject to a mutually reinforcing signal from ENSO suppression of wetland emissions and enhancement of biomass burning CH 4 (or vice versa), as well as positive feedbacks from OH dynamics, it is particularly suited to study the role of ENSO in the CH 4 cycle.
We find a sizeable effect of ENSO on biomass burning, as indicated by HCN variability in southern mid-latitudes.In contrast, ENSO explains a smaller fraction (≤ 37 %) of [CH 4 ] IAV, even in the southern tropics, where the expected effect should be greatest.Trends in [CH 4 ] and δ 13 CH 4 in these latitudes are far less influenced by ENSO (≤ 23 %).On hemispheric and global scales the ENSO signal in the methane records is similarly weak.Our results do not rule out that ENSO influences CH 4 emissions from wetlands and biomass burning through temperature, enhanced precipitation, or droughts in key regions, but any such impacts are overwhelmed by OH dynamics or other source and sink processes.We review literature estimates of ENSO-driven emissions and find them too small and sporadic to account for the post-2007 rise.Counteracting OH dynamics are expected to further dampen any influence ENSO may have on methane growth rates.Our findings suggest that ENSO is not an important driver for recent global trends in methane, including the [CH 4 ] plateau and the increase in [CH 4 ] since 2007.The latter must therefore have different causes.Our results do not rule out that wetland production is a contributor to the post-2007 [CH 4 ] rise if driven by environmental controls other than ENSO.This is suggested by an increase in wetland CH 4 production between the periods 2000-2006 and 2006-2014, although with the limited confidence of a single wetland emissions model (Zhang et al., 2018).The longer the atmospheric [CH 4 ] and δ 13 CH 4 trends persist, the less probable processes that impact IAV and short-lived cyclical events like ENSO as the driver are.Therefore, we consider increasing anthropogenic sources as the more likely cause of the [CH 4 ] rise.Changes in removal rates via OH have been suggested as an additional (Rigby et al., 2017) or alternative (Turner et al., 2017) driver of the increase, but recent work suggests that sink impacts are not dominant (Naus et al., 2018).There is evidence for additional methane emissions from agriculture (Wolf et al., 2017) and from fossil fuel sources (Hausmann et al., 2016); both may contribute to the current rise in [CH 4 ] (Worden et al., 2017).Further identification of these processes is necessary to inform climate change mitigation policies and climate projections.

Figure 2 .
Figure 2. Selected time series of ENSO indices and [CH 4 ].Panels from top to bottom: (a) Multivariate ENSO Index (MEI), Southern Oscillation Index (SOI), Ocean Niño Index (ONI), and El Modoki Index (EMI) shown for nominal literature data and their 12-month running means.(b) Global [CH 4 ] records; monthly means, 12month running mean, detrended 12-month running mean, as well as nominal and detrended growth rates.(c) [CH 4 ] records from BHD (d) [CH 4 ] records from ASC.(e) [CH 4 ] records from SMO. BHD, ASC, and SMO display the same records as for global time series.

Figure 3 .
Figure 3. Selected time series of ENSO indices, HCN, and δ 13 CH 4 .Panels from top to bottom: (a) multivariate ENSO Index (MEI), Southern Oscillation Index (SOI), Ocean Niño Index (ONI), and El Modoki Index (EMI) shown for nominal literature data and their 12-month running means.(b) HCN records as 12-month running means from LAU for total atmospheric column and stratosphere (12-100 km) and respective growth rates.(c) Global annually averaged δ 13 CH 4 according to Schaefer et al. (2016) updated to end of 2016; nominal and detrended values and their respective growth rates.(d) δ 13 CH 4 from BHD; monthly means, 12-month running mean, and detrended 12-month running mean, as well as nominal and detrended growth rates.(e) δ 13 CH 4 from ASC. (f) δ 13 CH 4 from SMO. ASC and SMO display same records as for BHD.Note scale differences between all δ 13 C axes to accentuate variability for comparison with ENSO.

Table 1 .
Description of time series products used in the correlation analyses.

Table 2 .
Spearman correlation of methane mole fraction with ENSO variability.Correlations (r 2 values) for the Spearman ranking coefficient between [CH 4 ] time series from various sites and ENSO indices with lag times (in months) for optimum results.Colour backgrounds indicate r 2 values in 10 % classes.Grey background indicates correlations with p values > 0.001.

Table 2 : Spearman correlation of methane mole fraction with ENSO variability.
Correlations (r 2 -values) for the Spearman ranking coefficient between [CH4] time series from various sites and ENSO indices with lag times (in months) for optimum results.Colour backgrounds indicate r 2 -values in 10% classes.Grey background indicates correlations with p-values > 0.001.

Table 3 .
Spearman correlation of δ 13 CH 4 and HCN with ENSO variability.Correlations (r 2 values) for the Spearman ranking coefficient between dependent variables, i.e. δ 13 CH 4 and HCN time series from various sites, and ENSO indices with lag times (in months) for optimum results.Colour backgrounds indicate r 2 values in 10 % classes.Grey background indicates correlations with p values > 0.001.
Prather et al., 2012)nriched (or vice versa), the two influences partly cancel out the combined emission rates, i.e. their impact on[CH 4].However, they reinforce each other's impact on total source δ 13 CH 4 .It is possible that biomass burning and wetland CH 4 production have different response times to ENSO forcing, which would weaken their cumulative impact on δ 13 CH 4 .Similarly, a longer atmospheric residence time of CH 4 (∼ 9 years,Prather et al., 2012)overHCN (∼ 3 months, Li et al., 2000)and a smaller relative portion of ENSO-sensitive emissions in the global methane source may lead to dampening effects

Table 4 .
Pearson correlation of δ 13 CH 4 and HCN with ENSO variability.Correlations (r 2 values) for the Pearson correlation coefficient between dependent variables, i.e. δ 13 CH 4 and HCN time series from various sites, and ENSO indices with lag times (in months) for optimum results.Colour backgrounds indicate r 2 values in 10 % classes.Results have not been screened for p values.ENSO-triggered plume.However, if ENSO is invoked as a main cause of recent trends in [CH 4 ] and δ 13 CH 4 correlation with ENSO indices compared to HCN.The available records for HCN and δ 13 CH 4 from ASC and SMO only cover a small number of ENSO events, which could affect the results.However, when analysing sub-periods of global and BHD[CH 4] and δ 13 CH 4 records, we find larger correlations for shorter periods, particularly when strong ENSO events are excluded.This shows that the results are not biased against the detection of ENSO influences because records are too short.We also note that all stations measure background air; they are set up to detect broad spatial and temporal trends and not specific emission events such as an www.biogeosciences.net/15/6371/2018/Biogeosciences, 15, 6371-6386, 2018 show that during La Niña events high surface temperatures over the western Pacific lead to upward transport over the Indonesian region (a CH 4 source area from wetlands and rice paddies) and negative CH 4 anomalies in the mid-troposphere (tropical surface air with relatively low[CH 4] replaces air from the Northern Hemisphere with higher [CH 4 ]).This mechanism would dampen the signal of higher La Niña emissions in surface records like SMO and ASC.However, the corresponding El Niño anomalies in midtropospheric CH 4 over the central Pacific are smaller.This indicates that central Pacific surface air, where there are no CH 4 sources, is closer in[CH 4] to mid-tropospheric levels than surface air from the western Pacific.Unless there were strong longitudinal differences in mid-tropospheric[CH 4], this is inconsistent with a scenario in which high concentrations of CH 4 are generated over the western Pacific in La Niñas but transported upwards and away from the surface stations used in this study.On hemispheric or global scales, transport processes are unlikely to play a strong role, given the short mixing time of methane relative to its atmospheric turnover.The low correlations of [CH 4 ] and δ 13 CH 4 with ENSO rule out a dominant role for ENSO-triggered sink changes in atmospheric methane records.Removal processes could lead to either amplification or dampening of source signals.