Abstract. We use 2010–2015 observations of atmospheric methane columns from the GOSAT satellite instrument in a global inverse analysis to improve estimates of methane emissions and their trends over the period, as well as the global concentration of tropospheric OH (the hydroxyl radical, methane's main sink) and its trend. Our inversion solves the Bayesian optimization problem analytically including closed-form characterization of errors. This allows us to (1) quantify the information content from the inversion towards optimizing methane emissions and its trends, (2) diagnose error correlations between constraints on emissions and OH concentrations, and (3) generate a large ensemble of solutions testing different assumptions in the inversion. We show how the analytical approach can be used, even when prior error standard deviation distributions are lognormal. Inversion results show large overestimates of Chinese coal emissions and Middle East oil and gas emissions in the EDGAR v4.3.2 inventory but little error in the United States where we use a new gridded version of the EPA national greenhouse gas inventory as prior estimate. Oil and gas emissions in the EDGAR v4.3.2 inventory show large differences with national totals reported to the United Nations Framework Convention on Climate Change (UNFCCC), and our inversion is generally more consistent with the UNFCCC data. The observed 2010–2015 growth in atmospheric methane is attributed mostly to an increase in emissions from India, China, and areas with large tropical wetlands. The contribution from OH trends is small in comparison. We find that the inversion provides strong independent constraints on global methane emissions (546 Tg a−1) and global mean OH concentrations (atmospheric methane lifetime against oxidation by tropospheric OH of 10.8±0.4 years), indicating that satellite observations of atmospheric methane could provide a proxy for OH concentrations in the future.

Using new satellite observations and atmospheric inverse modeling, we report methane emissions from the Permian Basin, which is among the world’s most prolific oil-producing regions and accounts for >30% of total U.S. oil production. Based on satellite measurements from May 2018 to March 2019, Permian methane emissions from oil and natural gas production are estimated to be 2.7 ± 0.5 Tg a−1, representing the largest methane flux ever reported from a U.S. oil/gas-producing region and are more than two times higher than bottom-up inventory-based estimates. This magnitude of emissions is 3.7% of the gross gas extracted in the Permian, i.e., ~60% higher than the national average leakage rate. The high methane leakage rate is likely contributed by extensive venting and flaring, resulting from insufficient infrastructure to process and transport natural gas. This work demonstrates a high-resolution satellite data–based atmospheric inversion framework, providing a robust top-down analytical tool for quantifying and evaluating subregional methane emissions.

Abstract. We evaluate the global atmospheric methane column retrievals from the new TROPOMI satellite instrument and apply them to a global inversion of methane sources for 2019 at 2∘ × 2.5∘ horizontal resolution. We compare the results to an inversion using the sparser but more mature GOSAT satellite retrievals and to a joint inversion using both TROPOMI and GOSAT. Validation of TROPOMI and GOSAT with TCCON ground-based measurements of methane columns, after correcting for retrieval differences in prior vertical profiles and averaging kernels using the GEOS-Chem chemical transport model, shows global biases of −2.7 ppbv for TROPOMI and −1.0 ppbv for GOSAT and regional biases of 6.7 ppbv for TROPOMI and 2.9 ppbv for GOSAT. Intercomparison of TROPOMI and GOSAT shows larger regional discrepancies exceeding 20 ppbv, mostly over regions with low surface albedo in the shortwave infrared where the TROPOMI retrieval may be biased. Our inversion uses an analytical solution to the Bayesian inference of methane sources, thus providing an explicit characterization of error statistics and information content together with the solution. TROPOMI has ∼ 100 times more observations than GOSAT, but error correlation on the 2∘ × 2.5∘ scale of the inversion and large spatial inhomogeneity in the number of observations make it less useful than GOSAT for quantifying emissions at that scale. Finer-scale regional inversions would take better advantage of the TROPOMI data density. The TROPOMI and GOSAT inversions show consistent downward adjustments of global oil–gas emissions relative to a prior estimate based on national inventory reports to the United Nations Framework Convention on Climate Change but consistent increases in the south-central US and in Venezuela. Global emissions from livestock (the largest anthropogenic source) are adjusted upward by TROPOMI and GOSAT relative to the EDGAR v4.3.2 prior estimate. We find large artifacts in the TROPOMI inversion over southeast China, where seasonal rice emissions are particularly high but in phase with extensive cloudiness and where coal emissions may be misallocated. Future advances in the TROPOMI retrieval together with finer-scale inversions and improved accounting of error correlations should enable improved exploitation of TROPOMI observations to quantify and attribute methane emissions on the global scale.

Abstract. We use satellite (GOSAT) and in situ (GLOBALVIEWplus CH4 ObsPack) observations of atmospheric methane in a joint global inversion of methane sources, sinks, and trends for the 2010–2017 period. The inversion is done by analytical solution to the Bayesian optimization problem, yielding closed-form estimates of information content to assess the consistency and complementarity (or redundancy) of the satellite and in situ data sets. We find that GOSAT and in situ observations are to a large extent complementary, with GOSAT providing a stronger overall constraint on the global methane distributions, but in situ observations being more important for northern midlatitudes and for relaxing global error correlations between methane emissions and the main methane sink (oxidation by OH radicals). The in-situ-only and the GOSAT-only inversions alone achieve 113 and 212 respective independent pieces of information (DOFS) for quantifying mean 2010–2017 anthropogenic emissions on 1009 global model grid elements, and respective DOFS of 67 and 122 for 2010–2017 emission trends. The joint GOSAT+ in situ inversion achieves DOFS of 262 and 161 for mean emissions and trends, respectively. Thus, the in situ data increase the global information content from the GOSAT-only inversion by 20 %–30 %. The in-situ-only and GOSAT-only inversions show consistent corrections to regional methane emissions but are less consistent in optimizing the global methane budget. The joint inversion finds that oil and gas emissions in the US and Canada are underestimated relative to the values reported by these countries to the United Nations Framework Convention on Climate Change (UNFCCC) and used here as prior estimates, whereas coal emissions in China are overestimated. Wetland emissions in North America are much lower than in the mean WetCHARTs inventory used as a prior estimate. Oil and gas emissions in the US increase over the 2010–2017 period but decrease in Canada and Europe. The joint inversion yields a global methane emission of 551 Tg a−1 averaged over 2010–2017 and a methane lifetime of 11.2 years against oxidation by tropospheric OH (86 % of the methane sink).

Abstract. We use 2010–2015 Greenhouse Gases Observing Satellite (GOSAT) observations of atmospheric methane columns over North America in a high-resolution inversion of methane emissions, including contributions from different sectors and their trends over the period. The inversion involves an analytical solution to the Bayesian optimization problem for a Gaussian mixture model (GMM) of the emission field with up to 0.5∘×0.625∘ resolution in concentrated source regions. The analytical solution provides a closed-form characterization of the information content from the inversion and facilitates the construction of a large ensemble of solutions exploring the effect of different uncertainties and assumptions in the inverse analysis. Prior estimates for the inversion include a gridded version of the Environmental Protection Agency (EPA) Inventory of US Greenhouse Gas Emissions and Sinks (GHGI) and the WetCHARTs model ensemble for wetlands. Our best estimate for mean 2010–2015 US anthropogenic emissions is 30.6 (range: 29.4–31.3) Tg a−1, slightly higher than the gridded EPA inventory (28.7 (26.4–36.2) Tg a−1). The main discrepancy is for the oil and gas production sectors, where we find higher emissions than the GHGI by 35 % and 22 %, respectively. The most recent version of the EPA GHGI revises downward its estimate of emissions from oil production, and we find that these are lower than our estimate by a factor of 2. Our best estimate of US wetland emissions is 10.2 (5.6–11.1) Tg a−1, on the low end of the prior WetCHARTs inventory uncertainty range (14.2 (3.3–32.4) Tg a−1), which calls for better understanding of these emissions. We find an increasing trend in US anthropogenic emissions over 2010–2015 of 0.4 % a−1, lower than previous GOSAT-based estimates but opposite to the decrease reported by the EPA GHGI. Most of this increase appears driven by unconventional oil and gas production in the eastern US. We also find that oil and gas production emissions in Mexico are higher than in the nationally reported inventory, though there is evidence for a 2010–2015 decrease in emissions from offshore oil production.