Articles | Volume 18, issue 16
Atmos. Chem. Phys., 18, 12257–12267, 2018
Atmos. Chem. Phys., 18, 12257–12267, 2018

Research article 24 Aug 2018

Research article | 24 Aug 2018

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

2010–2016 methane trends over Canada, the United States, and Mexico observed by the GOSAT satellite: contributions from different source sectors
Jian-Xiong Sheng1,2, Daniel J. Jacob1, Alexander J. Turner1, Joannes D. Maasakkers1, Joshua Benmergui1, A. Anthony Bloom3, Claudia Arndt2, Ritesh Gautam2, Daniel Zavala-Araiza2, Hartmut Boesch4,5, and Robert J. Parker4,5 Jian-Xiong Sheng et al.
  • 1School of Engineering and Applied Sciences, Harvard University, Cambridge, MA, USA
  • 2Environmental Defense Fund, Austin, TX, USA
  • 3Jet Propulsion Laboratory, California Institute of Technology, Pasadena, CA, USA
  • 4Leicester Institute for Space and Earth Observation, University of Leicester, Leicester, UK
  • 5NERC National Center for Earth Observation, UK

Correspondence: Jian-Xiong Sheng (


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 0.5×0.5 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 0.5×0.5 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 2.5±1.4  % a−1 (mean ± 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.

1 Introduction

Methane is an important greenhouse gas with a calculated climate impact as important as carbon dioxide over a 10-year time horizon (Myhre et al.2013; Etminan et al.2016). 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 (Kirschke et al.2013; Saunois et al.2016). Atmospheric methane concentrations leveled off in the 1990s but have been increasing again since 2007 (Dlugokencky et al.2009). 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 (Franco et al.2016; Hausmann et al.2016; Helmig et al.2016), (2) agriculture and/or wetlands based on isotopic information (Nisbet et al.2016; Schaefer et al.2016), (3) reduced biomass burning to reconcile the ethane and isotopic constraints (Worden et al.2017), and (4) declining concentrations of the OH radical (the main methane sink) based on the methylchloroform proxy (Rigby et al.2017; Turner et al.2017).

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 (Kuze et al.2016). Cressot et al. (2016) found that the GOSAT data had limited success in detecting regional year-to-year trends for 2009–2011. Turner et al. (2016) used a linear regression of GOSAT trends from January 2010 to January 2014 to infer a 2.8 % a−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. Bruhwiler et al. (2017) 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 Bruhwiler et al. (2017). 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.

2 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 µm absorption band). There has been no degradation of retrieval accuracy since the beginning of the record (Kuze et al.2016). 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.

Figure 1Frequency distributions and 2010–2016 trends of methane column average dry mole fractions XCH4 at Lamont, Oklahoma (36.6 N, 97.4 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.


Here we use the version 7.0 proxy nadir retrievals of GOSAT methane data from Parker et al. (2011, 2015). The proxy method uses prior knowledge of carbon dioxide columns, based on the MACC-II inversion product (v13r2; Chevallier et al.2010) accounting for seasonal and interannual variations, to infer methane column average dry mole fractions XCH4 (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 (Frankenberg et al.2006; Parker et al.2015). The resulting GOSAT XCH4 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 (Parker et al.2015; Buchwitz et al.2015, 2016). GOSAT observes in all seasons with near-uniform frequency south of 45 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 (Jacob et al.2016; Buchwitz et al.2017). Turner et al. (2016) estimated the CONUS background by using glint mode retrievals from GOSAT over the Pacific Ocean for the corresponding latitudes. Bruhwiler et al. (2017) 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 (0.5×0.5 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 0.5×0.5 grid cell, with seasonality removed using the seasonal-trend-loess (STL) decomposition method (Cleveland et al.1990). 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 (Goldstein et al.1995). 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 ≥10th 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.

Figure 2The 2010–2016 trends in GOSAT methane enhancements over North America. (a) Ordinary least-squares linear regression trends for 0.5×0.5 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 0.5×0.5 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 Turner et al. (2016).


Systematic errors of 4–6 ppb in GOSAT observations (Buchwitz et al.2016) 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 (Jacob et al.2016), 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 (Prather et al.2012; Kirschke et al.2013), 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 (Wunch et al.2011) network site at Lamont, Oklahoma (36.6 N, 97.4 W). Figure 1 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=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 0.1×0.1 spatial resolution from the gridded 2012 US EPA inventory of Maasakkers et al. (2016), the 2013 Canadian and 2010 Mexican oil–gas emission inventories of Sheng et al. (2017), and the EDGAR v4.2 global inventory for 2008 (European Commission2011) for other Canadian and Mexican sources. Compared to EDGAR v4.2, the more recent EDGAR v4.3.2 (Janssens-Maenhout et al.2017) 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 0.5×0.5 spatial resolution: (1) the mean of inventories contributing to the Wetland CH4 Intercomparison of Models Project (WETCHIMP) (Melton et al.2013) and (2) the 2010–2015 mean of the WetCHARTs extended ensemble wetland methane emissions inventory by Bloom et al. (2017). From these inventories we select high-emitting grid cells at 0.5×0.5 resolution dominated by a particular source sector. The high-emitting grid cells are defined as having emissions larger than 0.5 t h−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.

Figure 3Methane emissions in North America and contributions from different source sectors. (a) shows 0.5×0.5 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 Bloom et al. (2017) 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.


We define a total methane enhancement Δ for a given year, source sector, and country as

(1) Δ = i X ¯ CH 4 , i - X CH 4 , b , i ,

where X¯CH4,i is the annual mean value of the deseasonalized column average dry mole fractions in the 0.5×0.5 grid cell i for the given year, XCH4,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 XCH4,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 Δ, and this is used in what follows to test the statistical significance of year-to-year trends in Δ.

Figure 4National 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 Δ 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).


3 Results and discussion

Figure 2 (upper panel) shows the spatial distribution of GOSAT methane trends in local enhancements over North America at 0.5×0.5 spatial resolution from January 2010 to December 2016 (7 years of data). The 0.5×0.5 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 4×4 as in Turner et al. (2016) 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 Turner et al. (2016) is the choice of background. Enhancements defined relative to the Pacific background, as in Turner et al. (2016), are larger than in our approach where the background is defined locally.

Figure 5The 2010–2016 changes in methane-emitting activities. (a) Oil and natural gas production in CONUS (Drillinginfo2016). (b) Cattle population in Iowa, Kansas, Missouri, and Nebraska (USDA National Agricultural Statistics Service2015a), and animal units of swine in Iowa (Iowa Department of Natural Resources2017). One animal unit accounts for 3–5 heads of swine depending on body weight(USDA National Agricultural Statistics Service1995). (c) Total cattle population in Mexico (USDA Foreign Agricultural Service2015).


We improve the statistical significance of the CONUS enhancement trends by taking national statistics over all 0.5×0.5 grid cells. This is shown in the lower panels of Fig. 2 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 0.25±0.48 ppb a−1 (mean ± 1 standard deviation), which is statistically significant (sample size n=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 Janardanan et al. (2017). If this mean enhancement is taken as a measure of CONUS emissions, then a 0.25±0.48 ppb a−1 trend implies a 2.3±4 % a−1 increase in emissions for 2010–2016. The Turner et al. (2016) 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 3 shows the locations of high-emitting 0.5×0.5 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. 3 are those identified by both the WETCHIMP mean and Bloom et al. (2017) 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 (Sheng et al.2017), but are not directly detectable by GOSAT nadir measurements over land. Glint observations are available over the ocean but are much sparser.

Figure 4 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 Δ 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 ∼50 contributing 0.5×0.5 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 (Bloom et al.2017) (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 (4.0±1.6 % a−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 (2015) and shown in Fig. 5. 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 2.5±1.4 % a−1 derived from linear regression that is consistent with our previously calculated mean trend of 2.3 % a−1 averaged over the 0.5×0.5 gridded trends in Fig. 2. Breakdown by sector suggests that US oil–gas emissions increased at a marginally significant level (2.9 % a−1, p=0.03) from 2010 to 2016. Oil and unconventional (hydraulic fracturing) gas production grew by 15 % a−1 and 19 % a−1, respectively, during that period (Fig. 5), though production rate is not necessarily a predictor of emissions (Peischl et al.2015).

The US livestock emissions show a 3.5±1.8 % a−1 increase in our analysis, largely reflecting the agricultural Midwest where high-emitting grid cells are concentrated (Fig. 3). These grid cells emit 0.95 Tg CH4 a−1 from enteric fermentation (cattle) and 0.55 Tg CH4 a−1 from manure management (swine) according to the gridded EPA inventory (Maasakkers et al.2016). The cattle population in that region does not show a significant trend (Fig. 5), 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 (USDA National Agricultural Statistics Service2015b; Iowa Department of Natural Resources2017) (Fig. 5). This would increase swine manure management emissions by 0.02–0.1 Tg CH4 a−1 over the 2010–2016 period, assuming the IPCC (2006) emission factor of 10–45 kg CH4 head−1 a−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 (Iowa Department of Natural Resources2017). Those CAFOs tend to use liquid manure storage (US EPA2016) and have extended manure storage time (Iowa Department of Natural Resources2011), which lead to greater methane emissions. A recent bottom-up study from Wolf et al. (2017) 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 Bloom et al. (2017) 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 (Poulter et al.2017). Here the WetCHARTs extended ensemble used (i) GlobCover land cover data (Bontemps et al.2011) and the Global Lakes and Wetlands Database (GLWD; Lehner and Dölla2004) to represent spatial wetland extent and (ii) ERA-Interim precipitation to account for temporal wetland extent (Bloom et al.2017).

Figure 6Time series of the residuals (observed minus simulated methane concentrations) of the CarbonTracker-Lagrange (CT-L) CH4 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 N, 90.3 W), WBI (41.7 N, 91.4 W), and WKT (31.3 N, 97.3 W). Solid lines show the medians of NACP and GOSAT trends, and shaded areas show the 25th–75th percentile envelope.


Inverse analyses of methane concentrations in surface air measured as part of the North American Carbon Program (Wofsy and Harris2002) for 2010–2014 reveal no significant trends in US emissions over that period (Benmergui et al.2015). 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 (Benmergui et al.2015) 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 N, 90.3 W), WBI (West Branch, Iowa; 41.7 N, 91.4 W), and WKT (Moody, Texas; 31.3 N, 97.3 W). These sites are strongly influenced by large livestock and/or wetlands, livestock, and oil–gas sources, respectively (Benmergui et al.2015). 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.

4 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 +2.5±1.4% a−1 for the period and appears to reflect contributions from both oil–gas and livestock. Assuming 38–55 Tg CH4 a−1 for the CONUS emissions, including 29–40 Tg CH4 a−1 from anthropogenic sources (Miller et al.2013; Wecht et al.2014; Turner et al.2015; Maasakkers et al.2016) and 9–15 Tg CH4 a−1 from wetlands (Melton et al.2013; Bloom et al.2017), we deduce an increasing emission trend of 0.9–1.3 Tg CH4 a−1 over the 2010–2016 period, which would account for about 20 % of the global increase in atmospheric methane (Rigby et al.2017). 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 (Bruhwiler et al.2017), which better account for the role of atmospheric transport but have their own errors notably in the prior assumptions of emission patterns (Maasakkers et al.2016). Future inversions combining GOSAT and surface network data with improved bottom-up estimates are needed to provide more robust trend analyses.

Data availability

The GOSAT Proxy XCH4 data (Parker et al., 2015) used in this publication are available from (last access: 1 December 2017). The TCCON XCH4 data (Wunch et al., 2011) used in this publication are from (lass access: 1 December 2017). The North American Carbon Program (NACP) data (Wofsy and Harris, 2002) used in this publication are from (last access: 1 December 2017).


The supplement related to this article is available online at:

Author contributions

JXS and DJJ designed the research; JXS and DJJ performed the research; JXS analyzed the GOSAT data; JDM and AAB provided the bottom-up CH4 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.

Competing interests

The authors declare that they have no conflict of interest.


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 (

Edited by: Martin Heimann
Reviewed by: three anonymous referees


Benmergui, J., Andrews, A., Thoning, K., Trudeau, M., Miller, S., Dlugokencky, E., Bruhwiler, L., Masarie, K., Worthy, D., Sweeney, C., and others: Integrating diverse observations of North American CH4 into flux inversions in CarbonTrackerLagrange-CH4, AGU Fall Meeting Abstracts, 2015. a, b, c

Bloom, A. A., Bowman, K. W., Lee, M., Turner, A. J., Schroeder, R., Worden, J. R., Weidner, R., McDonald, K. C., and Jacob, D. J.: A global wetland methane emissions and uncertainty dataset for atmospheric chemical transport models (WetCHARTs version 1.0), Geosci. Model Dev., 10, 2141–2156,, 2017. a, b, c, d, e, f, g

Bontemps, S., Defourny, P., Bogaert, E., Arino, O., Kalogirou, V., and Perez, J.: GLOBCOVER 2009: Products Description and Validation Report, Tech. rep., ESA, 2, 2011. a

Bruhwiler, L. M., Basu, S., Bergamaschi, P., Bousquet, P., Dlugokencky, E., Houweling, S., Ishizawa, M., Kim, H.-S., Locatelli, R., Maksyutov, S., Montzka, S., Pandey, S., Patra, P. K., Petron, G., Saunois, M., Sweeney, C., Schwietzke, S., Tans, P., and Weatherhead, E. C.: US CH4 Emissions from Oil and Gas Production: Have Recent Large Increases Been Detected?, J. Geophys. Res.-Atmos., 122, 4070–4083,, 2017. a, b, c, d

Buchwitz, M., Reuter, M., Schneising, O., Boesch, H., Guerlet, S., Dils, B., Aben, I., Armante, R., Bergamaschi, P., Blumenstock, T., Bovensmann, H., Brunner, D., Buchmann, B., Burrows, J. P., Butz, A., Chédin, A., Chevallier, F., Crevoisier, C. D., Deutscher, N. M., Frankenberg, C., Hase, F., Hasekamp, O. P., Heymann, J., Kaminski, T., Laeng, A., Lichtenberg, G., De Mazière, M., Noël, S., Notholt, J., Orphal, J., Popp, C., Parker, R., Scholze, M., Sussmann, R., Stiller, G. P., Warneke, T., Zehner, C., Bril, A., Crisp, D., Griffith, D. W. T., Kuze, A., O'Dell, C., Oshchepkov, S., Sherlock, V., Suto, H., Wennberg, P., Wunch, D., Yokota, T., and Yoshida, Y.: The Greenhouse Gas Climate Change Initiative (GHG-CCI): Comparison and quality assessment of near-surface-sensitive satellite-derived CO2 and CH4 global data sets, Remote Sens. Environ., 162, 344–362,, 2015. a

Buchwitz, M., Dils, B., Boesch, H., Crevoisier, C., Detmers, D., Frankenberg, C., Hasekamp, O., Hewson, W., Laeng, A., Noël, S., and Notholt, J.: ESA Climate Change Initiative (CCI) Product Validation and Intercomparison Report (PVIR) for the Essential Climate Variable (ECV) Greenhouse Gases (GHG) for data set Climate Research Data Package No. 4 (CRDP# 4), Version, 5, 253, available at: (last access: 20 April 2018), 2016. a, b

Buchwitz, M., Schneising, O., Reuter, M., Heymann, J., Krautwurst, S., Bovensmann, H., Burrows, J. P., Boesch, H., Parker, R. J., Somkuti, P., Detmers, R. G., Hasekamp, O. P., Aben, I., Butz, A., Frankenberg, C., and Turner, A. J.: Satellite-derived methane hotspot emission estimates using a fast data-driven method, Atmos. Chem. Phys., 17, 5751–5774,, 2017. a

Chevallier, F., Ciais, P., Conway, T. J., Aalto, T., Anderson, B. E., Bousquet, P., Brunke, E. G., Ciattaglia, L., Esaki, Y., Fröhlich, M., Gomez A., Gomez-Pelaez A. J., Haszpra L., Krummel P. B., Langenfelds, R. L., Leuenberger, M., Machida, T., Maignan, F., Matsueda, H., Morguí, J. A., Mukai, H., Nakazawa, T., Peylin, P., Ramonet, M., Rivier, L., Sawa, Y., Schmidt, M., Steele, L. P., Vay, S. A., Vermeulen, A. T., Wofsy, S., and Worthy, D.: CO2 surface fluxes at grid point scale estimated from a global 21 year reanalysis of atmospheric measurements, J. Geophys. Res.-Atmos., 115, D21307,, 2010. a

Cleveland, R. B., Cleveland, W. S., and Terpenning, I.: STL: A seasonal-trend decomposition procedure based on loess, Journal of Official Statistics, 6, 3–73, 1990. a

Cressot, C., Pison, I., Rayner, P. J., Bousquet, P., Fortems-Cheiney, A., and Chevallier, F.: Can we detect regional methane anomalies? A comparison between three observing systems, Atmos. Chem. Phys., 16, 9089–9108,, 2016. a

Dlugokencky, E. J., Bruhwiler, L., White, J. W. C., Emmons, L. K., Novelli, P. C., Montzka, S. A., Masarie, K. A., Lang, P. M., Crotwell, A. M., Miller, J. B., and Gatti, L. V.: Observational constraints on recent increases in the atmospheric CH4 burden, Geophys. Res. Lett., 36, L18803,, 2009. a

Drillinginfo: Drillinginfo’s Production Data Platform – HPDI, available at: (lass access: 1 December 2017), 2016. a

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

European Commission: Emission Database for Global Atmospheric Research (EDGAR), release version 4.2, 2011. a

Franco, B., Mahieu, E., Emmons, L. K., Tzompa-Sosa, Z. A., Fischer, E. V., Sudo, K., Bovy, B., Conway, S., Griffin, D., Hannigan, J. W., Strong, K., and Walker, K. A.: Evaluating ethane and methane emissions associated with the development of oil and natural gas extraction in North America, Environ. Res. Lett., 11, 044010,, 2016. a

Frankenberg, C., Meirink, J. F., Bergamaschi, P., Goede, A. P. H., Heimann, M., Körner, S., Platt, U., van Weele, M., and Wagner, T.: Satellite chartography of atmospheric methane from SCIAMACHY on board ENVISAT: Analysis of the years 2003 and 2004, J. Geophys. Res., 111, D07303,, 2006. a

Goldstein, A. H., Wofsy, S. C., and Spivakovsky, C. M.: Seasonal variations of nonmethane hydrocarbons in rural New England: Constraints on OH concentrations in northern midlatitudes, J. Geophys. Res., 100, 21023–21033,, 1995. a

Hausmann, P., Sussmann, R., and Smale, D.: 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. a

Helmig, D., Rossabi, S., Hueber, J., Tans, P., Montzka, S. A., Masarie, K., Thoning, K., Plass-Duelmer, C., Claude, A., Carpenter, L. J., Lewis, A. C., Punjabi, S., Reimann, S., Vollmer, M. K., Steinbrecher, R., Hannigan, J. W., Emmons, L. K., Mahieu, E., Franco, B., Smale, D., and Pozzer, A.: Reversal of global atmospheric ethane and propane trends largely due to US oil and natural gas production, Nature Geosci.,, in press, 2016. a

Iowa Department of Natural Resources: Manure on Frozen and Snow-covered Ground: Report to the Governor and General Assembly, available at: (last access: 1 December 2017), 2011. a

Iowa Department of Natural Resources: The online Animal Feeding Operations database, available at: (last access: 1 December 2017), 2017. a, b, c

IPCC: 2006 IPCC Guidelines for National Greenhouse Gas Inventories (Volume 4): Agriculture, Forestry and Other Land Use, available at: (last access: 1 December 2017), 2006. a

Jacob, D. J., Turner, A. J., Maasakkers, J. D., Sheng, J., Sun, K., Liu, X., Chance, K., Aben, I., McKeever, J., and Frankenberg, C.: Satellite observations of atmospheric methane and their value for quantifying methane emissions, Atmos. Chem. Phys., 16, 14371–14396,, 2016. a, b

Janardanan, R., Maksyutov, S., Ito, A., Yukio, Y., and Matsunaga, T.: Assessment of Anthropogenic Methane Emissions over Large Regions Based on GOSAT Observations and High Resolution Transport Modeling, Remote Sensing, 9, 941,, 2017. a

Janssens-Maenhout, G., Crippa, M., Guizzardi, D., Muntean, M., Schaaf, E., Dentener, F., Bergamaschi, P., Pagliari, V., Olivier, J. G. J., Peters, J. A. H. W., van Aardenne, J. A., Monni, S., Doering, U., and Petrescu, A. M. R.: EDGAR v4.3.2 Global Atlas of the three major Greenhouse Gas Emissions for the period 1970–2012, Earth Syst. Sci. Data Discuss.,, 2017. a

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

Kuze, A., Suto, H., Shiomi, K., Kawakami, S., Tanaka, M., Ueda, Y., Deguchi, A., Yoshida, J., Yamamoto, Y., Kataoka, F., Taylor, T. E., and Buijs, H. L.: Update on GOSAT TANSO-FTS performance, operations, and data products after more than 6 years in space, Atmos. Meas. Tech., 9, 2445–2461,, 2016. a, b

Lehner, B. and Dölla, P.: Development and validation of a global database of lakes, reservoirs and wetlands, J. Hydrol., 296, 1–22,, 2004. a

Maasakkers, J. D., Jacob, D. J., Sulprizio, M. P., Turner, A. J., Weitz, M., Wirth, T., Hight, C., DeFigueiredo, M., Desai, M., Schmeltz, R., Hockstad, L., Bloom, A. A., Bowman, K. W., Jeong, S., and Fischer, M. L.: Gridded National Inventory of U.S. Methane Emissions, Environ. Sci. Technol., 50, 13123–13133,, 2016. a, b, c, d

Melton, J. R., Wania, R., Hodson, E. L., Poulter, B., Ringeval, B., Spahni, R., Bohn, T., Avis, C. A., Beerling, D. J., Chen, G., Eliseev, A. V., Denisov, S. N., Hopcroft, P. O., Lettenmaier, D. P., Riley, W. J., Singarayer, J. S., Subin, Z. M., Tian, H., Zürcher, S., Brovkin, V., van Bodegom, P. M., Kleinen, T., Yu, Z. C., and Kaplan, J. O.: Present state of global wetland extent and wetland methane modelling: conclusions from a model inter-comparison project (WETCHIMP), Biogeosciences, 10, 753–788,, 2013. a, b

Miller, S. M., Wofsy, S. C., Michalak, A. M., Kort, E. A., Andrews, A. E., Biraud, S. C., Dlugokencky, E. J., Eluszkiewicz, J., Fischer, M. L., Janssens-Maenhout, G., Miller, B. R., Miller, J. B., Montzka, S. A., Nehrkorn, T., and Sweeney, C.: Anthropogenic emissions of methane in the United States, P. Natl. Acad. Sci., 110, 20018–20022,, 2013. a

Myhre, G., Shindell, D., Bréon, F.-M., Collins, W., Fuglestvedt, J., Huang, J., Koch, D., Lamarque, J.-F., Lee, D., Mendoza, B., and Nakajima, T.: Anthropogenic and natural radiative forcing, Climate Change, 423, 658–40, 2013. a

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

Parker, R., Boesch, H., Cogan, A., Fraser, A., Feng, L., Palmer, P. I., Messerschmidt, J., Deutscher, N., Griffith, D. W. T., Notholt, J., Wennberg, P. O., and Wunch, D.: Methane observations from the Greenhouse Gases Observing SATellite: Comparison to ground-based TCCON data and model calculations, Geophys. Res. Lett., 38, L15807,, 2011. a

Parker, R. J., Boesch, H., Byckling, K., Webb, A. J., Palmer, P. I., Feng, L., Bergamaschi, P., Chevallier, F., Notholt, J., Deutscher, N., Warneke, T., Hase, F., Sussmann, R., Kawakami, S., Kivi, R., Griffith, D. W. T., and Velazco, V.: Assessing 5 years of GOSAT Proxy XCH4 data and associated uncertainties, Atmos. Meas. Tech., 8, 4785–4801,, 2015. a, b, c

Peischl, J., Ryerson, T. B., Aikin, K. C., de Gouw, J. A., Gilman, J. B., Holloway, J. S., Lerner, B. M., Nadkarni, R., Neuman, J. A., Nowak, J. B., Trainer, M., Warneke, C., and Parrish, D. D.: Quantifying atmospheric methane emissions from the Haynesville, Fayetteville, and northeastern Marcellus shale gas production regions, J. Geophys. Res.-Atmos., 120, 2119–2139,, 2015. a

Poulter, B., Bousquet, P., Canadell, J. G., Ciais, P., Peregon, A., Marielle Saunois, Arora, V. K., Beerling, D. J., Brovkin, V., Jones, C. D., Joos, F., Nicola Gedney, Ito, A., Kleinen, T., Koven, C. D., McDonald, K., Melton, J. R., Peng, C., Shushi Peng, Prigent, C., Schroeder, R., Riley, W. J., Saito, M., Spahni, R., Tian, H., Lyla Taylor, Viovy, N., Wilton, D., Wiltshire, A., Xu, X., Zhang, B., Zhang, Z., and Zhu, Q.: Global wetland contribution to 2000–2012 atmospheric methane growth rate dynamics, Environ. Res. Lett., 12, 094013,, 2017. a

Prather, M. J., Holmes, C. D., and Hsu, J.: Reactive greenhouse gas scenarios: Systematic exploration of uncertainties and the role of atmospheric chemistry, Geophys. Res. Lett., 39, L09803,, 2012. a

Rigby, M., Montzka, S. A., Prinn, R. G., White, J. W. C., Young, D., O'Doherty, S., Lunt, M. F., Ganesan, A. L., Manning, A. J., Simmonds, P. G., Salameh, P. K., Harth, C. M., Mühle, J., Weiss, R. F., Fraser, P. J., Steele, L. P., Krummel, P. B., McCulloch, A., and Park, S.: Role of atmospheric oxidation in recent methane growth, P. Natl. Acad. Sci., 114, 5373–5377,, 2017. a, b

Saunois, M., Bousquet, P., Poulter, B., Peregon, A., Ciais, P., Canadell, J. G., Dlugokencky, E. J., Etiope, G., Bastviken, D., Houweling, S., Janssens-Maenhout, G., Tubiello, F. N., Castaldi, S., Jackson, R. B., Alexe, M., Arora, V. K., Beerling, D. J., Bergamaschi, P., Blake, D. R., Brailsford, G., Brovkin, V., Bruhwiler, L., Crevoisier, C., Crill, P., Covey, K., Curry, C., Frankenberg, C., Gedney, N., Höglund-Isaksson, L., Ishizawa, M., Ito, A., Joos, F., Kim, H.-S., Kleinen, T., Krummel, P., Lamarque, J.-F., Langenfelds, R., Locatelli, R., Machida, T., Maksyutov, S., McDonald, K. C., Marshall, J., Melton, J. R., Morino, I., Naik, V., O'Doherty, S., Parmentier, F.-J. W., Patra, P. K., Peng, C., Peng, S., Peters, G. P., Pison, I., Prigent, C., Prinn, R., Ramonet, M., Riley, W. J., Saito, M., Santini, M., Schroeder, R., Simpson, I. J., Spahni, R., Steele, P., Takizawa, A., Thornton, B. F., Tian, H., Tohjima, Y., Viovy, N., Voulgarakis, A., van Weele, M., van der Werf, G. R., Weiss, R., Wiedinmyer, C., Wilton, D. J., Wiltshire, A., Worthy, D., Wunch, D., Xu, X., Yoshida, Y., Zhang, B., Zhang, Z., and Zhu, Q.: The global methane budget 2000–2012, Earth Syst. Sci. Data, 8, 697–751,, 2016. a

Schaefer, H., Fletcher, S. E. M., Veidt, C., Lassey, K. R., Brailsford, G. W., Bromley, T. M., Dlugokencky, E. J., Michel, S. E., Miller, J. B., Levin, I., Lowe, D. C., Martin, R. J., Vaughn, B. H., and White, J. W. C.: A 21st-century shift from fossil-fuel to biogenic methane emissions indicated by 13CH4, Science, 352, 80–84,, 2016. a

Sheng, J.-X., Jacob, D. J., Maasakkers, J. D., Sulprizio, M. P., Zavala-Araiza, D., and Hamburg, S. P.: A high-resolution (0.1×0.1) inventory of methane emissions from Canadian and Mexican oil and gas systems, Atmos. Environ., 158, 211–215,, 2017. a, b

Turner, A. J., Jacob, D. J., Wecht, K. J., Maasakkers, J. D., Lundgren, E., Andrews, A. E., Biraud, S. C., Boesch, H., Bowman, K. W., Deutscher, N. M., Dubey, M. K., Griffith, D. W. T., Hase, F., Kuze, A., Notholt, J., Ohyama, H., Parker, R., Payne, V. H., Sussmann, R., Sweeney, C., Velazco, V. A., Warneke, T., Wennberg, P. O., and Wunch, D.: Estimating global and North American methane emissions with high spatial resolution using GOSAT satellite data, Atmos. Chem. Phys., 15, 7049–7069,, 2015. a

Turner, A. J., Jacob, D. J., Benmergui, J., Wofsy, S. C., Maasakkers, J. D., Butz, A., Hasekamp, O., and Biraud, S. C.: 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. a, b, c, d, e, f, g

Turner, A. J., Frankenberg, C., Wennberg, P. O., and Jacob, D. J.: Ambiguity in the causes for decadal trends in atmospheric methane and hydroxyl, P. Natl. Acad. Sci., 114, 5367–5372,, 2017. a

USDA Foreign Agricultural Service: Mexico: Livestock and Products Annual, available at: (last access: 1 December 2017), 2015.  a, b

USDA National Agricultural Statistics Service: Animal manure management, RCA (Soil and Water Resources Conservation Act) issue brief #7, available at: (last access: 1 December 2017), 1995. a

USDA National Agricultural Statistics Service: Cattle Inventory, available at: (last access: 1 December 2017), 2015a. a

USDA National Agricultural Statistics Service: Hog Inventory, available at: (last access: 1 December 2017), 2015b. a

US EPA: Inventory of U.S. Greenhouse Gas Emissions and Sinks 1990–2014, 2016. a

Wecht, K. J., Jacob, D. J., Frankenberg, C., Jiang, Z., and Blake, D. R.: Mapping of North American methane emissions with high spatial resolution by inversion of SCIAMACHY satellite data, J. Geophys. Res.-Atmos., 119, 7741–7756,, 2014. a

Wofsy, S. C. and Harris, R. C.: The North American Carbon Program (NACP): Report of the NACP Committee of the U.S. Interagency Carbon Cycle Science Program, U.S. Global Change Res. Program, Washington, D. C., 2002. a

Wolf, J., Asrar, G. R., and West, T. O.: Revised methane emissions factors and spatially distributed annual carbon fluxes for global livestock, Carbon Balance and Management, 12, 16,, 2017. a

Worden, J. R., Bloom, A. A., Pandey, S., Jiang, Z., Worden, H. M., Walker, T. W., Houweling, S., and Röckmann, T.: Reduced biomass burning emissions reconcile conflicting estimates of the post-2006 atmospheric methane budget, Nature Commun., 8, 2227,, 2017. a

Wunch, D., Toon, G. C., Blavier, J.-F. L., Washenfelder, R. A., Notholt, J., Connor, B. J., Griffith, D. W. T., Sherlock, V., and Wennberg, P. O.: The Total Carbon Column Observing Network, Philosophical Transactions of the Royal Society of London A: Mathematical, Phys. Eng. Sci., 369, 2087–2112,, 2011. a

Short summary
Analysis of 7 years (2010–2016) of GOSAT methane trends over Canada, the contiguous US, and Mexico suggests that US methane emissions increased by 2.5 ± 1.4 % a−1 over the 7-year period, with contributions from both oil–gas systems and livestock in the Midwest. 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.
Final-revised paper