<iframe src="https://www.googletagmanager.com/ns.html?id=GTM-KCV32QR" height="0" width="0" style="display:none;visibility:hidden">
Quick search in Citations
*Required field
Quick Search in Journals
Quick Search in Journals
AnywhereCitationPNAS NexusProceedings of the National Academy of Sciences
Advanced Search
SUGGESTED SEARCHES:
Covid-19 Artificial Intelligence Mpox Climate Change Gun Violence
Research Article
Environmental Sciences
Open access

Ambiguity in the causes for decadal trends in atmospheric methane and hydroxyl

Alexander J. Turner [email protected], Christian Frankenberg https://orcid.org/0000-0002-0546-5857 [email protected], Paul O. Wennberg, and Daniel J. JacobAuthors Info & Affiliations
Edited by Mark H. Thiemens, University of California, San Diego, La Jolla, CA, and approved December 28, 2016 (received for review September 26, 2016)
April 17, 2017
114 (21) 5367-5372
https://doi.org/10.1073/pnas.1616020114
Commentary
May 15, 2017
Overexplaining or underexplaining methane’s role in climate change
Michael J. Prather, Christopher D. Holmes
Research Article
April 17, 2017
Role of atmospheric oxidation in recent methane growth
PDF/EPUB

Significance

Recent trends in atmospheric methane are not well understood as evidenced by multiple hypotheses proposed to explain the stabilization of methane concentrations in the early 2000s and the renewed growth since 2007. Here we use a multispecies inversion to determine the cause of these decadal trends. The most likely explanation for the renewed growth in atmospheric methane involves a decrease in hydroxyl (OH), the main sink for atmospheric methane, that is partially offset by a decrease in methane emissions. However, we also demonstrate that the problem of attributing methane trends from the current surface observation network, including isotopes, is underdetermined and does not allow unambiguous attribution of decadal trends.

Abstract

Methane is the second strongest anthropogenic greenhouse gas and its atmospheric burden has more than doubled since 1850. Methane concentrations stabilized in the early 2000s and began increasing again in 2007. Neither the stabilization nor the recent growth are well understood, as evidenced by multiple competing hypotheses in recent literature. Here we use a multispecies two-box model inversion to jointly constrain 36 y of methane sources and sinks, using ground-based measurements of methane, methyl chloroform, and the C13/C12 ratio in atmospheric methane (δ13CH4) from 1983 through 2015. We find that the problem, as currently formulated, is underdetermined and solutions obtained in previous work are strongly dependent on prior assumptions. Based on our analysis, the mathematically most likely explanation for the renewed growth in atmospheric methane, counterintuitively, involves a 25-Tg/y decrease in methane emissions from 2003 to 2016 that is offset by a 7% decrease in global mean hydroxyl (OH) concentrations, the primary sink for atmospheric methane, over the same period. However, we are still able to fit the observations if we assume that OH concentrations are time invariant (as much of the previous work has assumed) and we then find solutions that are largely consistent with other proposed hypotheses for the renewed growth of atmospheric methane since 2007. We conclude that the current surface observing system does not allow unambiguous attribution of the decadal trends in methane without robust constraints on OH variability, which currently rely purely on methyl chloroform data and its uncertain emissions estimates.
Atmospheric methane ( CH 4 ) is the second strongest anthropogenic greenhouse gas (1) and concentrations have been increasing for much of the past century (2) due, primarily, to increasing anthropogenic emissions. Atmospheric concentrations stabilized in the early 2000s (3) (hereafter referred to as the “methane stabilization”) and began increasing again in 2007 (4, 5) (hereafter referred to as the “renewed growth”). There has been much speculation about the cause of these trends (3–24). Attribution has proved to be a difficult task in part because this period of renewed growth is characterized by a methane growth rate of ∼ 6 ppb/y, which represents a source–sink imbalance of only 3% [or an increase of 20 Tg/y compared with an estimated annual source of 550 Tg/y (13)].
Previous work investigating the trends in atmospheric methane has generally used observations of either atmospheric ethane or bulk carbon isotope ratios in atmospheric methane ( δ 13 CH4), in conjunction with methane observations, to provide additional constraints on the sources of methane. This is because ethane is coemitted with methane from fossil-fuel sources, which represent ∼ 62% of the ethane budget (25), and has been used to infer changes in methane emissions from fossil-fuel sources. Similarly, δ 13 CH4 has been used to determine the sources governing atmospheric methane concentrations because different methane sources and sinks have distinct isotopic signatures.
However, previous works using ethane and δ 13 CH4 have come to differing conclusions about the causes of the stabilization in the early 2000s and the renewed growth since 2007. For example, Kai et al. (9) used isotope measurements and attributed the methane stabilization to a reduction in microbial sources whereas Aydin et al. (8) and Simpson et al. (12) used ethane observations and attributed it to a reduction in fossil-fuel sources. Levin et al. (11) found isotope measurements to be inconclusive. Similarly, Hausmann et al. (17), Franco et al. (18), and Helmig et al. (19) used ethane observations and attributed part of the renewed growth to fossil-fuel sources whereas Schaefer et al. (16), Rice et al. (22), Nisbet et al. (23), and Schwietzke et al. (24) reached differing conclusions using isotope measurements. Schaefer et al. (16) concluded that fossil fuels did not contribute to the renewed growth, Nisbet et al. (23) concluded that fossil fuels were not a dominant factor, Rice et al. (22) concluded that fugitive fossil-fuel sources have increased since 2000, and Schwietzke et al. (24) concluded that fossil-fuel sources have not increased. This prompts the question: Why do many, apparently, plausible explanations disagree with each other?
Quantitative attribution of methane emissions to fossil-fuel sources at global scales using ethane is hampered by the large variability in methane-to-ethane emission ratios and recent increases in ethane sources that emit little methane (26, 27). Similarly, there is a large overlap in the signatures from fossil-fuel and nonfossil methane sources. Part of this overlap is because fossil-fuel sources are not strictly thermogenic in origin with more than 20% of the world’s natural gas reserves coming from microbial activity (28–30). This overlap makes it difficult to draw quantitative conclusions about the methane sources, using atmospheric measurements of δ 13 CH4 on global scales (SI Appendix, section 1).
In addition, changes in the hydroxyl radical (OH), the main sink for atmospheric methane, complicate the issue. Previous work has used observations of methyl chloroform ( CH 3 CCl 3 ; abbreviated as “ MCF ”) to provide constraints on global-to-hemispheric OH concentrations because the loss of methyl chloroform is controlled by OH (31–34). Recent work has shown how small increases in the OH sink can explain most of the methane stabilization (21).
Here we present a simple two-box (Northern and Southern Hemisphere) model to investigate the cause of the methane stabilization and renewed growth. Fig. 1 shows a schematic of the two-box model. The model simulates annual hemispheric concentrations of CH 4 12 , CH 4 13 , methyl chloroform, and OH. We use atmospheric observations of methane, δ13CH4, and methyl chloroform to constrain annual hemispheric emissions of methane, the isotopic composition of the emissions, methyl chloroform emissions, and the OH abundance in a nonlinear, stochastic, Bayesian inversion.
Fig. 1.
Schematic of the two-box model. Inputs are annual hemispheric OH anomalies, methyl chloroform emissions, methane emissions, and δ13CH4 for the methane emissions. Outputs are annual hemispheric concentrations of methyl chloroform, methane, and the δ13CH4 of atmospheric methane. Interhemispheric exchange time is 1 y.

Results

Most Likely Solution.

The most likely solution, defined as the largest posterior probability sampled, found here is that the renewed growth is due to a decline in the OH sink, partially offset by a decrease in methane emissions. Similarly, the stabilization is explained by an increase in the OH sink offsetting an increase in methane emissions. Fig. 2 shows the most likely solution from our nonlinear inversion for the drivers of decadal trends in atmospheric methane and OH as well as the modeled methane, δ13CH4, and methyl chloroform concentrations. The posterior model accurately represents the observed concentrations of all three species. This shows how changes of a few percent in the methane sources and sinks can explain all of the observed hemispheric-scale variability in methane and methyl chloroform. The large uncertainty in isotopic signatures (SI Appendix, section 1) makes it difficult to draw quantitative conclusions from the isotopic composition of Northern and Southern Hemispheric methane sources.
Fig. 2.
Most likely solution. Left column shows observed (black triangles) and modeled (solid lines) concentrations of atmospheric CH4 (Top), δ 13 CH4 (Middle), and methyl chloroform (Bottom). The Northern Hemisphere is yellow and the Southern Hemisphere is blue. Right column shows the methane emissions (plotted as a deviation from the constant prior emissions; Top), the isotopic composition of the methane emissions (Middle), and the OH anomaly relative to a global mean concentration of 1 × 10 6 molecules (molec)/cm3 (Bottom).
In the most likely solution the renewed growth in methane is, counterintuitively, explained by a reduction in methane emissions. We find an ∼ 25-Tg/y decrease in global methane emissions from 2003 to 2016, during the period of renewed growth that started in 2007. This decrease in methane emissions is offset by an ∼ 7% decrease in global mean OH over the same period, which is constrained by methyl chloroform observations. The combination of this decrease in methane emissions and OH results in an increase in atmospheric methane concentrations, observed in both hemispheres.
In this solution, we find an ∼ 35-Tg/y increase in methane emissions from 1993 to 2003 and an ∼ 7% increase in global mean OH from 1991 to 2000, with the bulk of the methane emission increase occurring in the Northern Hemisphere. The rate of increase in the methane emissions slows after about 1998, resulting in relatively stable atmospheric methane concentrations. Northern Hemispheric isotope concentrations do not exhibit any systematic changes during this period.
The isotopic compositions of the Northern and Southern Hemispheric methane sources (Fig. 2, Middle Right) are, largely, decoupled from the changes seen in the methane emissions and OH. The Northern Hemispheric methane emissions undergo modest changes whereas the isotopic composition of those emissions fluctuate about −52.5‰. Conversely, the Southern Hemispheric methane emissions remain largely unchanged from 1980 to 2016 whereas the isotopic composition decreases by about 2‰  from 1990 (when publicly available isotope measurements began) to 2015.
The OH anomalies derived here are consistent with previous work examining global mean OH (21, 31–34). In particular, Montzka et al. (31), Rigby et al. (33), and McNorton et al. (21) used methyl chloroform observations to derive anomalies in global mean OH. Fig. 3, Top shows that their OH anomalies exhibit a similarity to the OH anomalies found here. Patra et al. (34) found that the interhemispheric ratio of OH has been roughly constant from 2004 to 2011 ( 0.97 ± 0.12 ). Prather et al. (32) used observations of methyl chloroform to derive a methane lifetime of 9.1 ± 0.9 y whereas our global mean methane lifetime from 1980 to 2016 is 9.2 ± 0.2 y. The consistency between the results presented here and past findings is not particularly surprising because our work uses many of the same methyl chloroform datasets and prior emissions estimates as previous work.
Fig. 3.
Analysis of OH anomalies and the methane lifetime from the most likely solution. Top is the same as Fig. 2, Bottom Right but includes the OH anomalies from Montzka et al. (31), Rigby et al. (33), and McNorton et al. (21) (black lines). OH anomalies from Montzka et al. (31), Rigby et al. (33), and McNorton et al. (21) are offset such that their mean matches the mean 1997–2007 anomaly found here. Middle is the ratio of Northern to Southern Hemispheric OH and the black line is from Patra et al. (34). Bottom is the methane lifetime in our two-box model and the black line is the lifetime from Prather et al. (32). OH is the only sink included in our two-box model so the methane lifetime shown here is more representative of the actual methane lifetime, not a lifetime due to OH loss.
This most likely solution is found to be robust to small perturbations in the prior error variance parameters for methane emissions and OH anomalies (SI Appendix, section 4.1), interhemispheric exchange times (SI Appendix, section 4.4), and alternate observation operators (SI Appendix, section 4.2). However, we find the amplitude of the changes in methane emissions and OH anomalies is strongly sensitive to the methyl chloroform reaction constant with OH (SI Appendix, section 4.5). Whereas the exact magnitude of our most likely solution changes in the different sensitivity tests, the general spatiotemporal pattern of increasing methane emissions and OH anomalies in the mid-1990s and decreasing methane emissions and OH anomalies from 2000 to present is robust to small perturbations but not large perturbations (as we present in the following sections).

Assuming Fixed OH Concentrations.

We performed a sensitivity test (Fig. 4) where the inversion assumed time-invariant OH concentrations (global mean OH concentration of 1 × 10 6 molec/cm3), thus linearizing the problem. The agreement with the observations is largely unchanged whereas methane emissions exhibit a fundamentally different temporal pattern. Here, we find that the renewed growth since 2007 is due to an ∼ 20-Tg/y increase in global methane emissions with each hemisphere contributing ∼ 10 Tg/y. The isotopic composition of the methane emissions is almost unchanged from the full nonlinear inversion, even though the methane emissions are radically different. MCF observations are reconciled through changes to the uncertain MCF emissions.
Fig. 4.
Sensitivity test with fixed OH concentrations. Details are the same as in Fig. 2. Dashed lines are from the most likely solution (Fig. 2).
The only discernible difference between the simulated concentrations in the full nonlinear case and this sensitivity test with fixed OH concentrations is in the first 5 y of the δ13CH4 concentrations when there are no observations (Fig. 4, Left). The consistency of the isotopic compositions indicates that simulating the δ13CH4 observations is largely unaffected by changes in the methane emissions or OH concentrations and that the δ13CH4 observations are providing constraints only on the isotopic compositions of the sources; it does not indicate that this spatiotemporal pattern in the isotopic compositions is a robust feature. This difference in the first 5 y of δ13CH4 concentrations is due to a slightly different treatment of the prior distribution for the initial conditions (Materials and Methods).

Assuming Fixed Methane Emissions.

As an extreme test, we performed an additional sensitivity study (Fig. 5) where the inversion assumed that methane emissions were time invariant (global methane emissions of 550 Tg/y), also linearizing the problem. Only modest changes to the OH concentrations are needed to explain the observed methane concentrations with fixed methane emissions and a small divergence in Northern and Southern Hemispheric OH can explain changes in the interhemispheric methane difference. The renewed growth is explained by an ∼ 3–5% decrease in global mean OH from 2005 to 2016. However, this sensitivity test does require larger changes to the methyl chloroform emissions (SI Appendix, Fig. S4), relative to the full nonlinear inversion. Again, the isotopic composition of the methane emissions is almost unchanged from the full nonlinear inversion and is identical to the isotopic composition of the methane emissions from the sensitivity test with fixed OH concentrations.
Fig. 5.
Sensitivity test with fixed methane emissions. Details are the same as in Fig. 2. Dashed lines are from the most likely solution (Fig. 2).
Assuming the methane emissions are fixed implicitly places a constraint on the magnitude of the interannual variability of the global mean OH concentration. This is because an increase (decrease) in the OH anomaly could be offset by an increase (decrease) in the methane emissions to satisfy the observational record. As such, assuming the methane emissions are a fixed parameter (as opposed to a parameter to be estimated) in an inversion solving for global mean OH limits the potential interannual variability of global mean OH.
There are no discernible differences between the simulated concentrations in the two sensitivity tests and, as with the fixed OH sensitivity test, the only difference between the simulated concentrations in the full nonlinear case and this sensitivity test is in the first 5 y of the δ13CH4 concentrations (Figs. 4, Left and 5, Left).

Discussion and Conclusions

We performed a nonlinear Bayesian inversion to infer the most likely set of drivers of decadal trends in atmospheric methane and OH. Based on our assumptions (Table 1), we find that decreasing OH concentrations is the most likely explanation for the renewed growth since 2007, with methane emissions actually decreasing during that period. This result is robust to small perturbations in our prior assumptions but not to large perturbations. The isotopic composition of the Southern Hemispheric methane sources in our most likely solution decreased by ∼ 1‰  during this period whereas the emissions decreased by ∼ 10 Tg/y. This would indicate that the isotopically heavy sources in the Southern Hemisphere (such as biomass burning) may have decreased whereas the isotopically light sources remained constant. There is evidence from the satellite record of CO pointing to a decrease in Southern Hemispheric biomass burning since 2001 (35) but some of the decrease in CO could be due to anthropogenic sources (36). The isotopic composition of the Northern Hemispheric sources remains relatively stable during this period, indicating negligible changes in the proportion of emissions coming from isotopically light and heavy sources. Finally, the most likely solution found here suggests that the recent increase in atmospheric methane may be a transient feature driven by fluctuations in OH concentrations.
Table 1.
Prior distributions for components of the state vector
Model input Distribution 𝐚 𝐛 𝝁 𝝈 τ , y
Annual components (emissions, isotopic compositions, and anomalies)            
NH CH 4 , Tg/y N B 300 500 412.5 20 5
SH CH 4 , Tg/y N B 100 250 137.5 20 5
NH δ 13 CH 4 , ‰ U −60.0 −45.0 — — —
SH δ 13 CH 4 , ‰ U −60.0 −45.0 — — —
NH MCF, Gg/y N B −1.0 β α γ 3
SH MCF, Gg/y U −1.0 1.0 — — —
NH OH, % N B −20 20 0 10 3
SH OH, % N B −20 20 0 10 3
Initial conditions (IC)
NH CH 4 , ppb U 1,540 1,620 — — —
SH CH 4 , ppb U 1,480 1,580 — — —
NH δ 13 CH 4 , ‰ U −48.2 −46.6 — — —
SH δ 13 CH 4 , ‰ U −48.2 −46.6 — — —
NH MCF, ppt U 15 135 — — —
SH MCF, ppt U 15 135 — — —
α is the updated annual MCF emissions from Prinn et al. (39), β =max ([1:0 Gg=y, 2α]), and γ=max ([0:2α, 1:5 Gg=y]). ppb, parts per billion; ppt, parts per trillion.
It is important to be cautious with source attribution based on isotope measurements. The isotopic composition of the Northern and Southern Hemispheric methane emissions remained largely unchanged in our three inversions (compare Figs. 2, Middle Right, 4, Middle Right, and 5, Middle Right) even though the spatiotemporal patterns in the methane emissions were radically different. As such, the interpretation of the sources driving the changes in methane emissions would differ. There is also a large overlap in the isotopic composition of different sources, further complicating the interpretation of the isotope measurements.
As for the methane stabilization, we find an ∼ 30-Tg/y increase in Northern Hemispheric methane emissions from 1992 to 2003 combined with an increasing OH anomaly is the most likely explanation for the methane trends from the early 1990s through the early 2000s; this generally agrees with the recent findings from McNorton et al. (21) who find that changes in OH can explain most of the methane stabilization. Based on our analysis, this seems more likely than previous work that has argued for the collapse of the former Soviet Union as a cause of the methane stabilization in the early 2000s (3, 8, 12, 16). However, we stress that our study relates only to aggregated (hemispheric-scale) drivers of the methane growth rate and does not preclude large subhemispheric-scale changes.
We performed two sensitivity tests where different potential drivers of decadal trends were held constant in the inversion. These sensitivity tests yielded two important conclusions: (i) Multiple (fundamentally different) scenarios can explain the observations and (ii) previous work that did not jointly estimate methane and OH aliased errors from one species to another.
For the former conclusion, both sensitivity tests are able to simulate the observations to within the observational uncertainties and the main difference between their likelihoods is due to the specification of the prior and assumptions in the analysis. As for the latter, previous work has rarely jointly estimated all parameters (e.g., methane emissions and OH concentrations) and is predisposed to a subset of solutions. For example, Schaefer et al. (16) derive a step increase of 19.7 Tg/y in methane emissions starting in 2007 and then attempt to explain the cause with isotope measurements. Their result is similar to our sensitivity test with fixed OH concentrations (Fig. 4) where we find an ∼ 20-Tg/y increase in methane emissions starting in 2007. However, this solution is not the most likely one if we allow OH concentrations to vary. A similar argument can be made against using fixed methane emissions when estimating global mean OH concentrations.
The methane stabilization and renewed growth can be reconciled through small changes to the sources and sinks (on the order of a few percent, relative to their global budgets). As such, small changes in the sources and sinks can have important implications for the observed atmospheric concentrations and make quantitative determination of the causes a difficult task. We find here that global methane emissions and OH likely changed by ± 7% over 36 y but multiple scenarios can explain the observed changes in atmospheric methane, δ13CH4, and methyl chloroform. As such, the apparent disagreement between past works on the causes of decadal trends in atmospheric methane (3–22) is, almost certainly, due to the problem being underdetermined.
Moving forward, stronger conclusions on the causes of decadal trends in atmospheric methane and OH could be drawn if we had other independent proxies for OH. Alternatively, a mechanistic explanation with supporting evidence for the potential changes in OH concentrations could allow us to draw stronger conclusions. For example, changes in J ( O 1 D ) due to stratospheric ozone could provide a large enough change in OH (37) but the changes in stratospheric ozone are only weakly correlated with the OH anomalies derived here (SI Appendix, section 5). Future work could also focus on obtaining higher spatially resolved constraints where the methane-to-ethane ratios and isotopic signatures of the sources are better constrained. This would also allow the use of more gradient information as could be discerned from satellite observations (15). More work is needed to discern the causes of decadal trends in atmospheric methane and OH.

Materials and Methods

The model and data used are available at https://github.com/alexjturner/BoxModel_PNAS_20161223.

Observational Records Used.

All datasets used are publicly available. Methane observations are from National Oceanic and Atmospheric Administration’s Earth System Research Laboratory (NOAA/ESRL). δ13CH4 observations are from NOAA/ESRL; University of Washington, Seattle; University of Heidelberg, Heidelberg; and University of California, Irvine, CA. Methyl chloroform observations are from NOAA/ESRL and the Global Atmospheric Gases Experiment (GAGE)/Advanced GAGE (AGAGE) network. See SI Appendix, section 2 for more information on the observations.

Bootstrapping Hemispheric Averages and Uncertainties.

We construct a hemispheric average atmospheric methane, δ13CH4, and methyl chloroform via bootstrapping from the different observational records. The observational records are deseasonalized with a site-specific stable seasonal filter and we require that each observational record has at least a 5 y of data. We then generate a hemispherically averaged observational record and uncertainty by randomly drawing n observational records from the population of possible records, where n is the total number of observational records for that hemisphere. These n records are then combined using a block-averaging scheme with a 1-y window. This process is repeated 50 times. We then compute the mean and uncertainty from the 50 different time series. We also impose minimum uncertainty of 2 ppb and 0.03‰  for methane and δ13CH4, respectively. Finally, we require that the uncertainty for older observations is greater than that for newer observations. See SI Appendix, section 2 for more information on the methodology.

Two-Box Model.

We use a two-box model with three species: CH 4 12 , CH 4 13 , and MCF . The model simulates annual hemispheric concentrations for each species and considers loss, through reaction with OH , and interhemispheric transport with a timescale of 1 y: τ NS = 1 y (shown schematically in Fig. 1).
This gives us a set of six coupled ordinary differential equations:
∂ [ X ] N ( t ) ∂ t = E X , N ( t ) − k [ X ] [ OH ] N ( t ) [ X ] N ( t ) + [ X ] S ( t ) − [ X ] N ( t ) τ NS ∂ [ X ] S ( t ) ∂ t = E X , S ( t ) − k [ X ] [ OH ] S ( t ) [ X ] S ( t ) + [ X ] N ( t ) − [ X ] S ( t ) τ NS ,
where [ X ] ( t ) and E X ( t ) are the hemispheric concentrations and emissions, respectively, for a given species. N and S subscripts denote Northern and Southern Hemispheres. We do not consider other minor loss terms for methane [e.g., methanotrophic bacteria in aerated soils, chlorine and atomic oxygen radicals in the stratosphere, or reactions with chlorine radicals from sea salt in the marine boundary layer (13)] or methyl chloroform [e.g., uptake by the ocean (38)]. Lifetimes and reaction rate constants ( k [ X ] ) used here should be viewed as total atmospheric lifetimes, not the lifetime with respect to OH loss. OH is plotted as an anomaly relative to a global mean concentration of 1 ×  10 6 molec/cm3.

Nonlinear, Stochastic, Bayesian Inversion.

The two-box model ( 𝐅 ) can be used to relate a state vector ( 𝐱 ) to the observations ( 𝐲 ),
𝐲 = 𝐅 ( 𝐱 ) + ϵ ,
but there will always be some error ( ϵ ) associated with the both the observations and the model. The state vector in our work is the annual hemispheric MCF emissions, OH anomalies, methane emissions, and their isotopic composition, as well as the initial conditions for methane, δ13CH4, and MCF. We can estimate the terms in this state vector using Bayesian inference,
P ( 𝐱 | 𝐲 ) ∝ P ( 𝐲 | 𝐱 ) ⋅ P ( 𝐱 ) ,
where P ( 𝐱 | 𝐲 ) is the posterior distribution, P ( 𝐲 | 𝐱 ) is the likelihood distribution, and P ( 𝐱 ) is the prior distribution. The two-box model is nonlinear due to the interaction between OH and methane: k [ CH 4 ] [ OH ] ( t ) [ CH 4 ] ( t ) . As such, we have adopted a stochastic method to infer the most likely solution.
We use the covariance matrix adaptation evolution strategy (CMA-ES) [Hansen (40) and references therein] to find the most likely solution. Typical sampling methods [such as Markov chain Monte Carlo (MCMC)] become prohibitively slow as the dimension of the state vector becomes large because they have trouble defining the proposal distribution. CMA-ES is an evolutionary algorithm that modifies the covariance matrix of the proposal distribution based on the fitness of multiple candidate solutions in a given generation. This allows CMA-ES to efficiently sample the posterior distribution. We restart CMA-ES with 10 different initializations and covariance matrices in an attempt to find a global minimum. In total, we draw 500,000,000 samples from the posterior distribution.
We assume the likelihood distribution is Gaussian with a diagonal covariance matrix populated by the uncertainties from the bootstrapping process. Because we are using a stochastic method, we can use non-Gaussian distributions that may be less restrictive and allow more flexible specification of the prior distribution. Our prior distribution is a convolution of uniform distributions: U ( 𝐚 , 𝐛 ) and bounded normal distributions (which can be written as the product of a normal distribution and a uniform distribution): N B ( 𝐚 , 𝐛 , 𝝁 , 𝝈 , τ ) = U ( 𝐚 , 𝐛 ) ⋅ N ( 𝝁 , 𝚺 ( 𝝈 , τ ) ) , where 𝐚 is the lower bound, 𝐛 is the upper bound, 𝝁 is the mean, 𝚺 is the covariance matrix, 𝝈 is the square root of the diagonal of 𝚺 , and τ is the temporal correlation length scale for 𝚺 . The prior distributions for each component of the state vector are shown in Table 1.

Linear, Gaussian, Bayesian Inversion for Sensitivity Tests.

The two-box model is nonlinear because of the interaction between OH and methane, as mentioned above. However, the model becomes linear if we assume that either methane or OH is fixed. As such, our sensitivity tests (presented in Figs. 4 and 5) have a linear response. For computational efficiency, we assume Gaussian errors in the sensitivity tests to obtain a closed-form solution for the posterior distribution [for example, Rodgers (41)],
𝐱 ^ = ( 𝐒 𝐚 − 1 + 𝐊 T 𝐒 𝐨 − 1 𝐊 ) − 1 𝐊 T 𝐒 𝐨 − 1 ( 𝐲 − 𝐊𝐱 ) ,
where 𝐱 ^ is the maximum a posteriori solution, 𝐒 𝐚 is the prior error covariance matrix, 𝐊 is the Jacobian matrix of the two-box model, and 𝐒 𝐨 is the observational error covariance matrix.

Data Availability

Data deposition: The model and data reported in this paper have been deposited in GitHub, https://github.com/alexjturner/BoxModel_PNAS_20161223.

Acknowledgments

We thank E. Dlugokencky for providing methane data; S. Montzka, R. Prinn, S. O’Doherty, and R. Weiss for providing MCF data; and I. Levin, C. Veidt, B. Vaughn, J. White, and S. Englund for providing δ13CH4 data. This work was supported by a Department of Energy Computational Science Graduate Fellowship (to A.J.T.) and by a NASA Carbon Monitoring System grant (to D.J.J.).

Supporting Information

Appendix (PDF)

References

1
; IPCC, Climate change 2013: The physical science basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, eds TF Stocker, D Qin, G-K Plattner, M Tignor, SK Allen, J Boschung, A Nauels, Y Xia, V Bex, PM Midgley (Cambridge Univ Press, Cambridge, UK, 2013).
Google Scholar
2
DM Etheridge, GI Pearman, PJ Fraser, Changes in tropospheric methane between 1841 and 1978 from a high accumulation-rate antarctic ice core. Tellus Ser B Chem Phys Meteorol 44, 282–294 (1992).
Crossref
Google Scholar
3
EJ Dlugokencky, Atmospheric methane levels off: Temporary pause or a new steady-state? Geophys Res Lett 30, 1992 (2003).
Crossref
Google Scholar
4
M Rigby, et al., Renewed growth of atmospheric methane. Geophys Res Lett 35, L22805 (2008).
Crossref
Google Scholar
5
EJ Dlugokencky, et al., Observational constraints on recent increases in the atmospheric CH4 burden. Geophys Res Lett 36, L18803 (2009).
Crossref
Google Scholar
6
JS Wang, et al., A 3-D model analysis of the slowdown and interannual variability in the methane growth rate from 1988 to 1997. Global Biogeochem Cycles 18, GB3011 (2004).
Crossref
Google Scholar
7
AM Fiore, LW Horowitz, EJ Dlugokencky, JJ West, Impact of meteorology and emissions on methane trends, 1990–2004. Geophys Res Lett 33, L12809 (2006).
Crossref
Google Scholar
8
M Aydin, et al., Recent decreases in fossil-fuel emissions of ethane and methane derived from firn air. Nature 476, 198–201 (2011).
Crossref
PubMed
Google Scholar
9
FM Kai, SC Tyler, JT Randerson, DR Blake, Reduced methane growth rate explained by decreased northern hemisphere microbial sources. Nature 476, 194–197 (2011).
Crossref
PubMed
Google Scholar
10
P Bousquet, et al., Source attribution of the changes in atmospheric methane for 2006–2008. Atmos Chem Phys 11, 3689–3700 (2011).
Crossref
Google Scholar
11
I Levin, et al., No inter-hemispheric δ 13 CH4 trend observed. Nature 486, E3–E4 (2012).
Crossref
PubMed
Google Scholar
12
IJ Simpson, et al., Long-term decline of global atmospheric ethane concentrations and implications for methane. Nature 488, 490–494 (2012).
Crossref
PubMed
Google Scholar
13
S Kirschke, et al., Three decades of global methane sources and sinks. Nat Geosci 6, 813–823 (2013).
Crossref
Google Scholar
14
I Pison, B Ringeval, P Bousquet, C Prigent, F Papa, Stable atmospheric methane in the 2000s: Key-role of emissions from natural wetlands. Atmos Chem Phys 13, 11609–11623 (2013).
Crossref
Google Scholar
15
AJ Turner, et al., A large increase in U.S. methane emissions over the past decade inferred from satellite data and surface observations. Geophys Res Lett 43, 2218–2224 (2016).
Crossref
Google Scholar
16
H Schaefer, et al., A 21st-century shift from fossil-fuel to biogenic methane emissions indicated by 13CH4. Science 352, 80–84 (2016).
Crossref
PubMed
Google Scholar
17
P Hausmann, R Sussmann, D Smale, Contribution of oil and natural gas production to renewed increase in atmospheric methane (2007–2014): Top-down estimate from ethane and methane column observations. Atmos Chem Phys 16, 3227–3244 (2016).
Crossref
Google Scholar
18
B Franco, et al., Evaluating ethane and methane emissions associated with the development of oil and natural gas extraction in North America. Environ Res Lett 11, 044010 (2016).
Crossref
Google Scholar
19
D Helmig, et al., Reversal of global atmospheric ethane and propane trends largely due to us oil and natural gas production. Nat Geosci 9, 490–495 (2016).
Crossref
Google Scholar
20
SB Dalsøren, et al., Atmospheric methane evolution the last 40 years. Atmos Chem Phys 16, 3099–3126 (2016).
Crossref
Google Scholar
21
J McNorton, et al., Role of OH variability in the stalling of the global atmospheric CH4 growth rate from 1999 to 2006. Atmos Chem Phys 16, 7943–7956 (2016).
Crossref
Google Scholar
22
AL Rice, et al., Atmospheric methane isotopic record favors fossil sources flat in 1980s and 1990s with recent increase. Proc Natl Acad Sci USA 113, 10791–10796 (2016).
Crossref
PubMed
Google Scholar
23
EG Nisbet, et al., Rising atmospheric methane: 2007-2014 growth and isotopic shift. Global Biogeochem Cycles 30, 1356–1370 (2016).
Crossref
Google Scholar
24
S Schwietzke, et al., Upward revision of global fossil fuel methane emissions based on isotope database. Nature 538, 88–91 (2016).
Crossref
PubMed
Google Scholar
25
Y Xiao, et al., Global budget of ethane and regional constraints on U.S. sources. J Geophys Res 113, D21306 (2008).
Crossref
Google Scholar
26
EA Kort, et al., Fugitive emissions from the Bakken shale illustrate role of shale production in global ethane shift. Geophys Res Lett 43, 4617–4623 (2016).
Crossref
Google Scholar
27
J Peischl, et al., Quantifying atmospheric methane emissions from oil and natural gas production in the Bakken shale region of North Dakota. J Geophys Res 121, 6101–6111 (2016).
Crossref
Google Scholar
28
DD Rice, GE Claypool, Generation, accumulation, and resource potential of biogenic gas. AAPG Bull 65, 5–25 (1981).
Google Scholar
29
AM Martini, JM Budai, LM Walter, M Schoell, Microbial generation of economic accumulations of methane within a shallow organic-rich shale. Nature 383, 155–158 (1996).
Crossref
Google Scholar
30
JB Curtis, Fractured shale-gas systems. AAPG Bull 86, 1921–1938 (2002).
Google Scholar
31
SA Montzka, et al., Small interannual variability of global atmospheric hydroxyl. Science 331, 67–69 (2011).
Crossref
PubMed
Google Scholar
32
MJ Prather, CD Holmes, J Hsu, Reactive greenhouse gas scenarios: Systematic exploration of uncertainties and the role of atmospheric chemistry. Geophys Res Lett 39, L09803 (2012).
Crossref
Google Scholar
33
M Rigby, et al., Re-evaluation of the lifetimes of the major CFCs and CH3CCl3 using atmospheric trends. Atmos Chem Phys 13, 2691–2702 (2013).
Crossref
Google Scholar
34
PK Patra, et al., Observational evidence for interhemispheric hydroxyl-radical parity. Nature 513, 219–23 (2014).
Crossref
PubMed
Google Scholar
35
HM Worden, et al., Decadal record of satellite carbon monoxide observations. Atmos Chem Phys 13, 837–850 (2013).
Crossref
Google Scholar
36
G Zeng, et al., Trends and variations in CO, C2H6, and HCN in the Southern hemisphere point to the declining anthropogenic emissions of CO and C2H6. Atmos Chem Phys 12, 7543–7555 (2012).
Crossref
Google Scholar
37
CD Holmes, MJ Prather, OA Søvde, G Myhre, Future methane, hydroxyl, and their uncertainties: Key climate and emission parameters for future predictions. Atmos Chem Phys 13, 285–302 (2013).
Crossref
Google Scholar
38
PO Wennberg, et al., Recent changes in the air-sea gas exchange of methyl chloroform. Geophys Res Lett 31, L16112 (2004).
Crossref
Google Scholar
39
RG Prinn, et al., A history of chemically and radiatively important gases in air deduced from ALE/GAGE/AGAGE. J Geophys Res 105, 17751–17792 (2000).
Crossref
Google Scholar
40
N Hansen, The CMA Evolution Strategy: A comparing review. Towards a New Evolutionary Computation. Advances on Estimation of Distribution Algorithms, eds JA Lozano, P Larraaga, I Inza, E Bengoetxea (Springer, Berlin), pp. 75–102 (2006).
Google Scholar
41
CD Rodgers Inverse Methods for Atmospheric Sounding (World Scientific, Singapore, 2000).
Crossref
Google Scholar

Information & Authors

Information

Published in

Go to Proceedings of the National Academy of Sciences
Go to Proceedings of the National Academy of Sciences
Proceedings of the National Academy of Sciences
Vol. 114 | No. 21
May 23, 2017
PubMed: 28416668

Classifications

  1. Physical Sciences
  2. Environmental Sciences

Copyright

Freely available online through the PNAS open access option.

Data Availability

Data deposition: The model and data reported in this paper have been deposited in GitHub, https://github.com/alexjturner/BoxModel_PNAS_20161223.

Submission history

Published online: April 17, 2017
Published in issue: May 23, 2017

Keywords

  1. methane
  2. renewed growth
  3. hydroxyl
  4. oxidative capacity
  5. troposphere

Acknowledgments

We thank E. Dlugokencky for providing methane data; S. Montzka, R. Prinn, S. O’Doherty, and R. Weiss for providing MCF data; and I. Levin, C. Veidt, B. Vaughn, J. White, and S. Englund for providing δ13CH4 data. This work was supported by a Department of Energy Computational Science Graduate Fellowship (to A.J.T.) and by a NASA Carbon Monitoring System grant (to D.J.J.).

Notes

This article is a PNAS Direct Submission.
See Commentary on page 5324.

Authors

Affiliations

Alexander J. Turner1 [email protected]
School of Engineering and Applied Sciences, Harvard University, Cambridge, MA 02138;
View all articles by this author
Christian Frankenberg1 https://orcid.org/0000-0002-0546-5857 [email protected]
Division of Geological and Planetary Sciences, California Institute of Technology, Pasadena, CA 91125;
Jet Propulsion Laboratory, California Institute of Technology, Pasadena, CA 21109
View all articles by this author
Paul O. Wennberg
Division of Geological and Planetary Sciences, California Institute of Technology, Pasadena, CA 91125;
View all articles by this author
Daniel J. Jacob
School of Engineering and Applied Sciences, Harvard University, Cambridge, MA 02138;
View all articles by this author

Notes

1
To whom correspondence may be addressed. Email: [email protected] or [email protected].
Author contributions: A.J.T. and C.F. designed research; A.J.T. and C.F. performed research; A.J.T. contributed new reagents/analytic tools; A.J.T., C.F., P.O.W., and D.J.J. analyzed data; and A.J.T., C.F., P.O.W., and D.J.J. wrote the paper.

Competing Interests

The authors declare no conflict of interest.

Metrics & Citations

Metrics

Note: The article usage is presented with a three- to four-day delay and will update daily once available. Due to ths delay, usage data will not appear immediately following publication. Citation information is sourced from Crossref Cited-by service.


Citation statements




Altmetrics

Citations

If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. Simply select your manager software from the list below and click Download.

Cited by

    View Options

    View options

    PDF format

    Download this article as a PDF file

    DOWNLOAD PDF

    Get Access

    Login options

    Check if you have access through your login credentials or your institution to get full access on this article.

    Personal login Institutional Login

    Recommend to a librarian

    Recommend PNAS to a Librarian
    Save for later Item saved, go to cart

    Purchase options

    Purchase this article to get full access to it.

    Single Article Purchase

    Ambiguity in the causes for decadal trends in atmospheric methane and hydroxyl
    Proceedings of the National Academy of Sciences
    • Vol. 114
    • No. 21
    • pp. 5319-E4317

    Restore content access

    Restore content access for purchases made as a guest

    Media

    Figures

    Tables

    Other

    Share

    Share

    Share article link

    Share on social media

    Further reading in this issue

    Trending