2010–2016 methane trends over Canada, the United States, and Mexico observed by the GOSAT satellite: contributions from different source sectors

Abstract

We use 7 years (2010–2016) of methane column observations from the Greenhouse Gases Observing Satellite (GOSAT) to examine trends in atmospheric methane concentrations over North America and infer trends in emissions. Local methane enhancements above background are diagnosed in the GOSAT data on a $\mathrm{0.5}{}^{\circ }\times \mathrm{0.5}{}^{\circ }$ grid by estimating the local background as the low (10th–25th) percentiles of the deseasonalized frequency distributions of the data for individual years. Trends in methane enhancements on the $\mathrm{0.5}{}^{\circ }\times \mathrm{0.5}{}^{\circ }$ grid are then aggregated nationally and for individual source sectors, using information from state-of-science bottom-up inventories. We find that US methane emissions increased by $\mathrm{2.5}\pm \mathrm{1.4}\mathit{\%}$ a${}^{-\mathrm{1}}$ (mean $\pm$ 1 standard deviation) over the 7-year period, with contributions from both oil–gas systems (possibly unconventional oil–gas production) and from livestock in the Midwest (possibly swine manure management). Mexican emissions show a decrease that can be attributed to a decreasing cattle population. Canadian emissions show year-to-year variability driven by wetland emissions and correlated with wetland areal extent. The US emission trends inferred from the GOSAT data account for about 20 % of the observed increase in global methane over the 2010–2016 period.

Introduction

Methane is an important greenhouse gas with a calculated climate impact as important as carbon dioxide over a 10-year time horizon . Livestock, oil–gas, and waste are the leading anthropogenic sources. Wetlands are the dominant natural source. Contributions from different source sectors and regions remain poorly quantified . Atmospheric methane concentrations leveled off in the 1990s but have been increasing again since 2007 . Interpretations of atmospheric observations from surface networks have reached conflicting conclusions as to the cause of the renewed increase, with attributions to (1) natural gas production based on correlation with ethane , (2) agriculture and/or wetlands based on isotopic information , (3) reduced biomass burning to reconcile the ethane and isotopic constraints , and (4) declining concentrations of the OH radical (the main methane sink) based on the methylchloroform proxy .

Satellite-based observations of atmospheric methane columns have been available from the TANSO-FTS instrument aboard the Greenhouse Gases Observing Satellite (GOSAT) continuously since May 2009 . found that the GOSAT data had limited success in detecting regional year-to-year trends for 2009–2011. used a linear regression of GOSAT trends from January 2010 to January 2014 to infer a 2.8 % a${}^{-\mathrm{1}}$ increase in methane emissions from the contiguous United States (CONUS). Their analysis was based on the trend in the CONUS enhancement of methane relative to the Pacific Ocean taken as background. argued that such a trend inference could have been biased by the brevity of the GOSAT record, atmospheric transport variability, seasonal bias in GOSAT sampling frequency, and the use of Pacific data as background. They also pointed out that global inversions of the surface network data for 2000–2012 from the North American Carbon Program (NACP) reveal no significant CONUS emission trend.

Here we reexamine the trend in CONUS emissions implied by statistical analysis of the GOSAT data, addressing the concerns expressed by . We use a longer record (January 2010–December 2016), an improved definition of the background that accounts for atmospheric transport variability, and we also evaluate the trends for consistency with an inversion using the NACP surface network data. We further extend the trend analysis to Canada and Mexico, attribute the trends to specific methane source sectors using information from new emission inventories, and relate the trends to independent data on sectoral activities.

Methods

GOSAT was launched in January 2009 in a Sun-synchronous low Earth orbit. It retrieves the atmospheric methane column by nadir measurements of solar backscatter (1.65 $\mathrm{\mu }$m absorption band). There has been no degradation of retrieval accuracy since the beginning of the record . Observations in the standard mode are made at three circular pixels of 10 km diameter across the orbit track 260 km apart, separated by 260 km along the track. The same locations are sampled every 3 days, making for a temporally dense data set at those locations. The observations often switch from the standard mode to focus on targets and this affects the regularity of the sampling.

Frequency distributions and 2010–2016 trends of methane column average dry mole fractions $X_{{\mathrm{CH}}_{\mathrm{4}}}$ at Lamont, Oklahoma (36.6${}^{\circ }$ N, 97.4${}^{\circ }$ W), as measured by TCCON and GOSAT. The upper panels show the deseasonalized 2015 frequency distributions from TCCON and GOSAT. The percentiles (blue plus signs) are plotted on a normal probability scale such that a normal distribution would plot as a straight line (green). The local background is defined by the 10th–25th percentile range and the mean annual local enhancement relative to this background is defined by the difference with the mean of the distribution. Lower panels compare TCCON and GOSAT backgrounds and enhancements for 2010–2016, with error standard deviations on the enhancements as described in the text.

[Figure omitted. See PDF]

Here we use the version 7.0 proxy nadir retrievals of GOSAT methane data from . The proxy method uses prior knowledge of carbon dioxide columns, based on the MACC-II inversion product v13r2; accounting for seasonal and interannual variations, to infer methane column average dry mole fractions $X_{{\mathrm{CH}}_{\mathrm{4}}}$ (in ppb) from the ratio of retrieved methane and carbon dioxide columns. The proxy method takes advantage of the much larger variability in methane than in carbon dioxide mixing ratios . The resulting GOSAT $X_{{\mathrm{CH}}_{\mathrm{4}}}$ data have been validated against the ground-based Total Carbon Column Observing Network (TCCON) and found to be of high quality with a single-scene precision of 0.7 % (random error) and a systematic error of 4–6 ppb . GOSAT observes in all seasons with near-uniform frequency south of 45${}^{\circ }$ N (CONUS and Mexico), but observations further north (Canada) are biased toward summer. The number of successful retrievals over Canada is 2–3 times less in winter than in summer (see Supplement).

From a simple mass balance perspective, enhancements of column methane above the surrounding background in a strong source region can be linearly related to the emissions in that region . estimated the CONUS background by using glint mode retrievals from GOSAT over the Pacific Ocean for the corresponding latitudes. pointed out that changes in large-scale meridional transport could alias trends in this background estimate onto trends in the emissions.

Here we define local background methane for a given CONUS location ($\mathrm{0.5}{}^{\circ }\times \mathrm{0.5}{}^{\circ }$ grid cell, typically including a single repeated GOSAT measurement location) and for a given year as the low (10th–25th) percentiles of the deseasonalized GOSAT methane observations within the given $\mathrm{0.5}{}^{\circ }\times \mathrm{0.5}{}^{\circ }$ grid cell, with seasonality removed using the seasonal-trend-loess (STL) decomposition method . This approach assumes that the low percentiles of concentrations reflect meteorological conditions where local sources have relatively little effect on methane concentrations due to rapid ventilation. Low percentiles are a standard approach for estimating the regional background at a measurement location . By choosing the 10th–25th percentile rather than a lower extreme we guard against the effect of measurement noise (random error). A permutation resampling test shows that GOSAT observations across North America are sufficiently precise that observations $\ge \mathrm{10}$th percentile are not affected by measurement noise (see Supplement). We use the range defined by the 10th–25th percentile range as a measure of uncertainty in the background for the purpose of determining the enhancement.

The 2010–2016 trends in GOSAT methane enhancements over North America. (a) Ordinary least-squares linear regression trends for $\mathrm{0.5}{}^{\circ }\times \mathrm{0.5}{}^{\circ }$ grid cells with sufficient GOSAT observations, where the deseasonalized annual mean methane enhancements are defined relative to a local low-percentile background as described in the text. The trends are not statistically significant at that resolution (see text). (b, c) Spatial frequency distributions for the $\mathrm{0.5}{}^{\circ }\times \mathrm{0.5}{}^{\circ }$ grid cells over the contiguous United States (CONUS) of mean methane and local background (b), and local methane enhancements computed by difference (c). The dashed black line in the lower-right panel indicates the mean trend in CONUS enhancements. Also shown in the lower-left panel are the 2010–2013 trend distributions from .

[Figure omitted. See PDF]

Systematic errors of 4–6 ppb in GOSAT observations do not affect the enhancement because the bias can be expected to similarly affect all percentiles of the methane observations. Local enhancements are inversely proportional to wind speed , but we find no significant trends in wind speeds over the 2010–2016 period that would contribute to our aggregated trends in methane enhancements (see Supplement). Any trends in OH concentrations would also not affect the enhancement because the lifetime of methane against oxidation is 9–10 years , which is very long compared to the timescale for ventilation from the source region.

We examined the validity of our approach by comparing frequency distributions of GOSAT methane columns and related trends to continuous ground-based column observations available from the TCCON network site at Lamont, Oklahoma (36.6${}^{\circ }$ N, 97.4${}^{\circ }$ W). Figure  shows the frequency distributions of the deseasonalized GOSAT and TCCON observations at Lamont. The GOSAT background defined by the 10th–25th percentiles is consistent with TCCON; we see that the repeated observation strategy of GOSAT at its discrete sampling locations makes for a sufficiently dense data set for defining the 10th–25th percentiles with little effect from instrument noise. The local annual mean background increases between 2010 and 2016 in a consistent way in the GOSAT and TCCON data sets, reflecting the global increase in the methane background. The enhancements above background also show comparable 2010–2016 trends between the two data sets, although the error standard deviations defined by the ranges of the 10th–25th percentiles are large and the trends at this single site are marginally significant ($p=\mathrm{0.07}$). Below we will use enhancement statistics aggregated over a large number of sites in order to reduce that uncertainty and quantify trends.

To aggregate trends in methane enhancements for individual source sectors, we use bottom-up annual mean sectoral information with $\mathrm{0.1}{}^{\circ }\times \mathrm{0.1}{}^{\circ }$ spatial resolution from the gridded 2012 US EPA inventory of , the 2013 Canadian and 2010 Mexican oil–gas emission inventories of , and the EDGAR v4.2 global inventory for 2008 for other Canadian and Mexican sources. Compared to EDGAR v4.2, the more recent EDGAR v4.3.2 has similar national totals and spatial patterns for non-oil–gas anthropogenic methane emissions in North America. For wetlands, we use multiyear annual mean values from two climatological inventories with $\mathrm{0.5}{}^{\circ }\times \mathrm{0.5}{}^{\circ }$ spatial resolution: (1) the mean of inventories contributing to the Wetland ${\mathrm{CH}}_{\mathrm{4}}$ Intercomparison of Models Project (WETCHIMP) and (2) the 2010–2015 mean of the WetCHARTs extended ensemble wetland methane emissions inventory by . From these inventories we select high-emitting grid cells at $\mathrm{0.5}{}^{\circ }\times \mathrm{0.5}{}^{\circ }$ resolution dominated by a particular source sector. The high-emitting grid cells are defined as having emissions larger than 0.5 t h${}^{-\mathrm{1}}$, encompassing 80 %–90 % of anthropogenic and wetland emissions in all three countries. A high-emitting grid cell is identified as dominated by a given source sector if that source sector accounts for more than 70 % of the total emissions in the cell. This allows us to define grid cells dominated specifically by oil–gas, livestock, waste, and wetland emissions. Contributions from other sectors (up to 30 %) may lead to some smoothing of results. Wetland-dominated areas determined by the WETCHIMP mean and WetCHARTs inventories differ significantly (see Supplement). Using either of the two inventories alone may bias our results, and thus we conservatively require wetland-dominated areas to be determined as such in both inventories.

Methane emissions in North America and contributions from different source sectors. (a) shows $\mathrm{0.5}{}^{\circ }\times \mathrm{0.5}{}^{\circ }$ grid cells with high emissions dominated by a particular sector as identified by the bottom-up inventories (see text for details). High-emitting wetland areas are those identified by both the WETCHIMP mean inventory and the mean inventory. Livestock includes enteric fermentation and manure management. Oil/gas includes the complete systems from production to distribution. Waste includes landfills and wastewater plants. (b) shows national emissions for 2008–2013 from the bottom-up inventories. “Other” includes smaller sources from coal, rice, combustion, petrochemical production, ferroalloy production, and biomass burning. Total emissions in the high-emitting wetland areas are indicated by gridded areas.

[Figure omitted. See PDF]

We define a total methane enhancement $\mathrm{\Delta }$ for a given year, source sector, and country as $$\mathrm{\Delta }=\sum _i\left({\overline{X}}_{{\mathrm{CH}}_{\mathrm{4}},i}-X_{{\mathrm{CH}}_{\mathrm{4}},\mathrm{b},i}\right),$$ where ${\overline{X}}_{{\mathrm{CH}}_{\mathrm{4}},i}$ is the annual mean value of the deseasonalized column average dry mole fractions in the $\mathrm{0.5}{}^{\circ }\times \mathrm{0.5}{}^{\circ }$ grid cell $i$ for the given year, $X_{{\mathrm{CH}}_{\mathrm{4}},\mathrm{b},i}$ is the corresponding local background value, and the summation is over all high-emitting grid cells for that sector and country. We require grid cells to have at least eight valid retrievals for a given year (to avoid bias in the annual mean value and percentiles), and about 70 % of grid cells meet this requirement. To account for local background variation due to atmospheric transport, the summation in Eq. (1) is conducted for 1000 Monte Carlo realizations where the background $X_{{\mathrm{CH}}_{\mathrm{4}},\mathrm{b},i}$ for each grid cell and for individual years is obtained by random sampling of percentiles in the 10th–25th range. Results are only weakly sensitive to the choice of that range (see Supplement). The resulting summation statistics define the probability density function of the total enhancement $\mathrm{\Delta }$, and this is used in what follows to test the statistical significance of year-to-year trends in $\mathrm{\Delta }$.

National trends in methane emissions since 2010 inferred from GOSAT, and contributions from specific source sectors where sufficient data are available. The trends are defined by relative year-to-year changes in the summed methane enhancements $\mathrm{\Delta }$ relative to the local backgrounds as computed from Eq. (1), and vertical bars are standard deviations derived from uncertainty in the local background (see text).

[Figure omitted. See PDF]

Results and discussion

Figure (upper panel) shows the spatial distribution of GOSAT methane trends in local enhancements over North America at $\mathrm{0.5}{}^{\circ }\times \mathrm{0.5}{}^{\circ }$ spatial resolution from January 2010 to December 2016 (7 years of data). The $\mathrm{0.5}{}^{\circ }\times \mathrm{0.5}{}^{\circ }$ trends are inferred from ordinary least-squares linear regression of the enhancements for individual years. The trends are not statistically significant at that resolution. We will aggregate grid cells in what follows to increase statistical significance. Some areas are sparsely sampled, such as California, while the central US is more densely observed due to a more regular schedule of standard measurements. This sampling bias in California is unlikely to affect the trend estimate on the national scale because California accounts for less than 5 % of the US national total according to the bottom-up inventories. Spatial averaging to $\mathrm{4}{}^{\circ }\times \mathrm{4}{}^{\circ }$ as in does not improve significance (see Supplement) because methane emissions are not correlated on that scale. A major reason for the weaker statistical significance of our results relative to is the choice of background. Enhancements defined relative to the Pacific background, as in , are larger than in our approach where the background is defined locally.

The 2010–2016 changes in methane-emitting activities. (a) Oil and natural gas production in CONUS . (b) Cattle population in Iowa, Kansas, Missouri, and Nebraska , and animal units of swine in Iowa . One animal unit accounts for 3–5 heads of swine depending on body weight. (c) Total cattle population in Mexico .

[Figure omitted. See PDF]

We improve the statistical significance of the CONUS enhancement trends by taking national statistics over all $\mathrm{0.5}{}^{\circ }\times \mathrm{0.5}{}^{\circ }$ grid cells. This is shown in the lower panels of Fig.  with the CONUS frequency distribution of trends in mean methane, local background, and the enhancements computed by difference. The mean 2010–2016 trend in methane enhancements over CONUS is $\mathrm{0.25}\pm \mathrm{0.48}$ ppb a${}^{-\mathrm{1}}$ (mean $\pm$ 1 standard deviation), which is statistically significant (sample size $n=\mathrm{241}$ and $p$ value $<$ 0.01). The mean 2010 methane enhancement for high-emitting grid cells in CONUS relative to local background is 10.8 ppb, which is comparable to that found by . If this mean enhancement is taken as a measure of CONUS emissions, then a $\mathrm{0.25}\pm \mathrm{0.48}$ ppb a${}^{-\mathrm{1}}$ trend implies a $\mathrm{2.3}\pm \mathrm{4}\mathit{\%}$ a${}^{-\mathrm{1}}$ increase in emissions for 2010–2016. The frequency distributions, shown in the lower-left panel, are much broader than ours because they did not use annual averaging of the data. Their Pacific background distribution is similarly broader and is also lower than our local background, which is appropriately elevated by continental influences. Below we will use the aggregated enhancements by source sectors (Eq. 1) to infer the trends and reduce the uncertainty.

Figure shows the locations of high-emitting $\mathrm{0.5}{}^{\circ }\times \mathrm{0.5}{}^{\circ }$ grid cells dominated by different sectors as identified by the bottom-up inventories of Sect. 2. Also shown are national emission totals from these inventories. Wetland-dominated areas in Fig.  are those identified by both the WETCHIMP mean and inventories in order to avoid false positives. There is clear separation of grid cells dominated by wetlands, oil–gas, and livestock source sectors. Waste emissions dominate in urban areas but are more localized. Offshore oil–gas emissions over the Gulf of Mexico account for more than 50 % of Mexican oil–gas total , but are not directly detectable by GOSAT nadir measurements over land. Glint observations are available over the ocean but are much sparser.

Figure shows GOSAT methane enhancement trends for 2010–2016 (expressed as percent change since 2010) over Canada, CONUS, and Mexico, along with contributions from the sector-resolved high-emitting grid cells. Here the trends are calculated for the summed enhancement $\mathrm{\Delta }$ in Eq. (1) calculated for individual years and for high-emitting grid cells of individual countries or high-emitting sectors. Inferring significant trends for a given source sector generally requires $\sim \mathrm{50}$ contributing $\mathrm{0.5}{}^{\circ }\times \mathrm{0.5}{}^{\circ }$ grid cells. The largest source of uncertainty is the selection of the local background within the 10th–25th percentile range, and this is reflected by the error bars in the figure.

The Canadian methane emissions show no significant 7-year trend but large year-to-year variability driven by wetlands. The 2014 maximum can be explained by a maximum of wetland areal extent (see Fig. S6 in Supplement). Observations in the oil–gas-dominated region of Canada (mainly natural gas in Alberta) are too sparse for inferring a significant oil–gas emission trend and are not shown here.

Mexican national emissions (excluding oil–gas offshore emissions) show a 5 %–10 % decrease over the 2010 to 2016 period that appears to be largely driven by livestock. The decrease in livestock emissions ($\mathrm{4.0}\pm \mathrm{1.6}$ % a${}^{-\mathrm{1}}$) is consistent with the 17 % decrease in the Mexican cattle population over that period as reported by the Foreign Agriculture Service of the US Department of Agriculture () and shown in Fig. . The slight increase in Mexican emissions from 2012 on suggests an increasing source to compensate for the declining livestock emissions but GOSAT observations are too sparse to identify that source.

The CONUS data imply a significant increase in methane emissions from 2010 to 2016, with a trend of $\mathrm{2.5}\pm \mathrm{1.4}$ % a${}^{-\mathrm{1}}$ derived from linear regression that is consistent with our previously calculated mean trend of 2.3 % a${}^{-\mathrm{1}}$ averaged over the $\mathrm{0.5}{}^{\circ }\times \mathrm{0.5}{}^{\circ }$ gridded trends in Fig. . Breakdown by sector suggests that US oil–gas emissions increased at a marginally significant level (2.9 % a${}^{-\mathrm{1}}$, $p=\mathrm{0.03}$) from 2010 to 2016. Oil and unconventional (hydraulic fracturing) gas production grew by 15 % a${}^{-\mathrm{1}}$ and 19 % a${}^{-\mathrm{1}}$, respectively, during that period (Fig. ), though production rate is not necessarily a predictor of emissions .

The US livestock emissions show a $\mathrm{3.5}\pm \mathrm{1.8}\mathit{\%}$ a${}^{-\mathrm{1}}$ increase in our analysis, largely reflecting the agricultural Midwest where high-emitting grid cells are concentrated (Fig. ). These grid cells emit 0.95 Tg ${\mathrm{CH}}_{\mathrm{4}}$ a${}^{-\mathrm{1}}$ from enteric fermentation (cattle) and 0.55 Tg ${\mathrm{CH}}_{\mathrm{4}}$ a${}^{-\mathrm{1}}$ from manure management (swine) according to the gridded EPA inventory . The cattle population in that region does not show a significant trend (Fig. ), but the swine population in Iowa (accounting for most of the swine population in the Midwest) increased by two million heads from 2010 to 2016 (Fig. ). This would increase swine manure management emissions by 0.02–0.1 Tg ${\mathrm{CH}}_{\mathrm{4}}$ a${}^{-\mathrm{1}}$ over the 2010–2016 period, assuming the emission factor of 10–45 kg ${\mathrm{CH}}_{\mathrm{4}}$ head${}^{-\mathrm{1}}$ a${}^{-\mathrm{1}}$. Here a larger value of the emission factor is more likely. The emission factor may have increased during that time due to an increase in swine body weight and a 30 % rise in concentrated animal feeding operations (CAFOs) with more than 1000 animal units . Those CAFOs tend to use liquid manure storage and have extended manure storage time , which lead to greater methane emissions. A recent bottom-up study from found a steady increasing trend since the 1990s in US methane emissions from manure management.

US wetland emissions do not show a significant trend over 2010–2016 but large year-to-year variability, which contributes in part to the total national trend after 2012. Correlation with driving variables in the WetCHARTs yearly ensemble of suggests that this year-to-year variability is related to wetland areal extent, same as for Canada (see Fig. S6 in Supplement), though the definition of wetland areal extent may vary significantly . Here the WetCHARTs extended ensemble used (i) GlobCover land cover data and the Global Lakes and Wetlands Database GLWD; to represent spatial wetland extent and (ii) ERA-Interim precipitation to account for temporal wetland extent .

Time series of the residuals (observed minus simulated methane concentrations) of the CarbonTracker-Lagrange (CT-L) ${\mathrm{CH}}_{\mathrm{4}}$ transport model simulations driven by posterior emissions optimized for NACP data (green) and scaled to GOSAT-inferred emission trends (purple) for three surface sites particularly sensitive to emissions from different sectors: LEF (45.9${}^{\circ }$ N, 90.3${}^{\circ }$ W), WBI (41.7${}^{\circ }$ N, 91.4${}^{\circ }$ W), and WKT (31.3${}^{\circ }$ N, 97.3${}^{\circ }$ W). Solid lines show the medians of NACP and GOSAT trends, and shaded areas show the 25th–75th percentile envelope.

[Figure omitted. See PDF]

Inverse analyses of methane concentrations in surface air measured as part of the North American Carbon Program for 2010–2014 reveal no significant trends in US emissions over that period . We examined whether the trends inferred from this work (significant trends after 2012) are consistent with the information provided by NACP surface data. For this purpose, we examined the residuals (observed minus simulated methane concentrations) of the CarbonTracker-Lagrange (CT-L) methane transport model driven with two sets of emissions: (1) the CT-L posterior emissions for 2010–2014 that are optimized to match all NACP data and show no significant trend and (2) a scaled version of the CT-L posterior emissions that matches the sector-resolved trends derived in this work. Figure 6 shows annual statistics and trends of the residuals for both simulations at three NACP sites included in the CT-L inversion: LEF (Park Falls, Wisconsin; 45.9${}^{\circ }$ N, 90.3${}^{\circ }$ W), WBI (West Branch, Iowa; 41.7${}^{\circ }$ N, 91.4${}^{\circ }$ W), and WKT (Moody, Texas; 31.3${}^{\circ }$ N, 97.3${}^{\circ }$ W). These sites are strongly influenced by large livestock and/or wetlands, livestock, and oil–gas sources, respectively . There is no significant trend in the residuals of the CT-L simulation driven by either our GOSAT-inferred emission trends or CT-L posterior emissions, and the two sets of residuals are statistically indistinguishable. We find similar results for other NACP sites that are less sensitive to source regions. This implies that the trends found in this work are compatible with the constraints provided by NACP data. This also suggests that the surface data may be spatially too sparse to adequately infer trends of the magnitude detected by GOSAT.

Conclusions

In conclusion, analysis of 7 years (2010–2016) of GOSAT methane trends over Canada, the contiguous US (CONUS), and Mexico suggests a significant increase in total US methane emissions after 2012 and decrease in Mexican emissions. The Mexican decreasing trend appears to be due to a declining cattle population. Canada shows no significant long-term trend but large year-to-year variability associated with wetlands and correlated with variations in wetland areal extent, though this trend is weighted toward summer because of the seasonal bias in observation frequency (less observations in winter). The US trend is $+\mathrm{2.5}\pm \mathrm{1.4}\mathit{\%}$ a${}^{-\mathrm{1}}$ for the period and appears to reflect contributions from both oil–gas and livestock. Assuming 38–55 Tg ${\mathrm{CH}}_{\mathrm{4}}$ a${}^{-\mathrm{1}}$ for the CONUS emissions, including 29–40 Tg ${\mathrm{CH}}_{\mathrm{4}}$ a${}^{-\mathrm{1}}$ from anthropogenic sources and 9–15 Tg ${\mathrm{CH}}_{\mathrm{4}}$ a${}^{-\mathrm{1}}$ from wetlands , we deduce an increasing emission trend of 0.9–1.3 Tg ${\mathrm{CH}}_{\mathrm{4}}$ a${}^{-\mathrm{1}}$ over the 2010–2016 period, which would account for about 20 % of the global increase in atmospheric methane . Our analysis is mainly limited by the length of the GOSAT record, and a longer record can provide more reliable results. The definition of local background may also not fully account for the variation in atmospheric transport. Our trend analysis should be compared to trends inferred from inverse modeling , which better account for the role of atmospheric transport but have their own errors notably in the prior assumptions of emission patterns . Future inversions combining GOSAT and surface network data with improved bottom-up estimates are needed to provide more robust trend analyses.

The GOSAT Proxy $X_{{\mathrm{CH}}_{\mathrm{4}}}$ data (Parker et al., 2015) used in this publication are available from http://www.esa-ghg-cci.org/sites/default/files/documents/public/documents/GHG-CCI_DATA.html (last access: 1 December 2017). The TCCON $X_{{\mathrm{CH}}_{\mathrm{4}}}$ data (Wunch et al., 2011) used in this publication are from https://tccondata.org/ (lass access: 1 December 2017). The North American Carbon Program (NACP) data (Wofsy and Harris, 2002) used in this publication are from https://www.nacarbon.org/ (last access: 1 December 2017).

The supplement related to this article is available online at: https://doi.org/10.5194/acp-18-12257-2018-supplement.

JXS and DJJ designed the research; JXS and DJJ performed the research; JXS analyzed the GOSAT data; JDM and AAB provided the bottom-up ${\mathrm{CH}}_{\mathrm{4}}$ emission inventories; JB performed CarbonTracker-Lagrange model simulations; HB and RJP provided the GOSAT data; JXS and DJJ wrote the paper; all the authors discussed the results and contributed to the manuscript.

The authors declare that they have no conflict of interest.

Acknowledgements

This work was supported by the NASA Earth Science Division and by the Environmental Defense Fund. Part of the funding for this study was provided through NASA grant no. NNH14ZDA001N-CMS. Jian-Xiong Sheng and Claudia Arndt were partially funded by the Kravis Scientific Research Fund at the Environmental Defense Fund. Funding for EDF's work on livestock methane was provided by Sue and Steve Mandel. Alexander J. Turner was supported by a Department of Energy (DOE) Computational Science Graduate Fellowship (CSGF). Part of this research was carried out at the Jet Propulsion Laboratory, California Institute of Technology, under a contract with the National Aeronautics and Space Administration. Robert J. Parker was funded via an ESA Living Planet Fellowship with additional funding from the UK National Centre for Earth Observation (NCEO) and the ESA Greenhouse Gas Climate Change Initiative (GHG-CCI). We thank the Japanese Aerospace Exploration Agency, National Institute for Environmental Studies, and the Ministry of Environment for the GOSAT data and their continuous support as part of the Joint Research Agreement. This research used the ALICE High Performance Computing Facility at the University of Leicester. TCCON data were obtained from the TCCON Data Archive, hosted by CaltechData (http://tccondata.org). Edited by: Martin Heimann Reviewed by: three anonymous referees