Search

We quantify weekly methane emissions at 0.25∘ × 0.3125∘ (≈25 × 25 km2) resolution from the Permian Basin, the largest oil production basin in the US, by inverse analysis of satellite observations from the TROPOspheric Monitoring Instrument (TROPOMI) from May 2018 to October 2020. The mean oil and gas emission from the region (± standard deviation of weekly estimates) was 3.7 ± 0.9 Tg a−1, higher than previous TROPOMI inversion estimates that may have used biased prior emissions or background assumptions. We find strong week-to-week variability in emissions superimposed on longer-term trends, and these are consistent with independent inferences of temporal emission variability from tower, aircraft, and multispectral satellite data. New well development and natural gas spot price were significant drivers of variability in emissions over our study period but the concurrent 50 % increase in oil and gas production was not. The methane intensity (methane emitted per unit of methane gas produced) averaged 4.6 % ± 1.3 % and steadily decreased from 5 %–6 % in 2018 to 3 %–4 % in 2020. While the decreasing trend suggests improvement in operator practices during the study period, methane emissions from the Permian Basin remained high, with methane intensity an order of magnitude above the industry target of <0.2 %. Our success in using TROPOMI satellite observations for weekly estimates of emissions from a major oil production basin shows promise for application to near-real-time monitoring in support of climate change mitigation efforts.

We quantify methane emissions in China and the contributions from different sectors by inverse analysis of 2019 TROPOMI satellite observations of atmospheric methane. The inversion uses as prior estimate the national sector-resolved anthropogenic emission inventory reported by the Chinese government to the United Nations Framework Convention on Climate Change (UNFCCC) and thus serves as a direct evaluation of that inventory. Emissions are optimized with a Gaussian mixture model (GMM) at up to 0.25° × 0.3125° resolution. The optimization is done analytically assuming lognormally distributed errors on prior emissions. Errors and information content on the optimal estimates are obtained directly from the analytical solution and also through a 36-member inversion ensemble. Our optimal estimate for total anthropogenic emissions in China is 65.0 (57.7–68.4) Tg a-1, where parentheses indicate uncertainty range. Contributions from individual sectors include 16.6 (15.6–17.6) Tg a-1 for coal, 2.3 (1.8–2.5) for oil, 0.29 (0.23–0.32) for gas, 17.8 (15.1–21.0) for livestock, 9.3 (8.2–9.9) for waste, 11.9 (10.7–12.7) for rice paddies, and 6.7 (5.8–7.1) for other sources. Our estimate is 21 % higher than the Chinese inventory reported to the UNFCCC (53.6 Tg a-1), reflecting upward corrections to emissions from oil (+147 %), gas (+61 %), livestock (+37 %), waste (+41 %), and rice paddies (+34 %), but downward correction for coal (-15 %). It is also higher than previous inverse studies (43–62 Tg a-1) that used the much sparser GOSAT satellite observations and were conducted at coarser resolution. We are in particular better able to separate coal and rice emissions. Our higher livestock emissions are attributed largely to northern China where GOSAT has little sensitivity. Our higher waste emissions reflect at least in part a rapid growth in wastewater treatment in China. Underestimate of oil emissions in the UNFCCC report appears to reflect unaccounted super-emitting facilities. Gas emissions in China are mostly from distribution, in part because of low emission factors from production and in part because 42 % of the gas is imported. Our estimate of emissions per unit of domestic gas production indicates a low life-cycle loss rate of 1.7 (1.3–1.9) %, which would imply net climate benefits from the current coal-to-gas energy transition in China. However, this small loss rate is somewhat misleading considering China’s high gas imports, including from Turkmenistan where emission per unit of gas production is very high.

We present a versatile, powerful, and user-friendly chemical data assimilation toolkit for simultaneously optimizing emissions and concentrations of chemical species based on atmospheric observations from satellites or suborbital platforms. The CHemistry and Emissions REanalysis Interface with Observations (CHEEREIO) exploits the GEOS-Chem chemical transport model and a localized ensemble transform Kalman filter algorithm (LETKF) to determine the Bayesian optimal (posterior) emissions and/or concentrations of a set of species based on observations and prior information, using an easy-to-modify configuration file with minimal changes to the GEOS-Chem or LETKF code base. The LETKF algorithm readily allows for non-linear chemistry and produces flow-dependent posterior error covariances from the ensemble simulation spread. The object-oriented Python-based design of CHEEREIO allows users to easily add new observation operators such as for satellites. CHEEREIO takes advantage of the HEMCO modular structure of input data management in GEOS-Chem to update emissions from the assimilation process independently from the GEOS-Chem code. It can seamlessly support GEOS-Chem version updates and is adaptable to other chemical transport models with similar modular input data structure. A postprocessing suite combines ensemble output into consolidated NetCDF files and supports a wide variety of diagnostic data and visualizations. We demonstrate CHEEREIO’s capabilities with an out-of-the-box application, assimilating global methane emissions and concentrations at weekly temporal resolution and 2°x2.5° spatial resolution for 2019 using TROPOMI satellite observations. CHEEREIO achieves a 50-fold improvement in computational performance compared to the equivalent analytical inversion of TROPOMI observations.

Satellite observations of dry column methane mixing ratios (XCH4) from shortwave infrared (SWIR) solar backscatter radiation provide a powerful resource to quantify methane emissions in service of climate action. The TROPOMI instrument launched in October 2017 provides global daily coverage at 5.5 × 7 km2 nadir pixel resolution but its retrievals can suffer from biases associated with SWIR surface albedo, scattering from aerosols and cirrus clouds, and across-track variability (striping). The GOSAT instrument launched in 2009 uses a retrieval method that is less subject to biases, but its data density is 200 times sparser than TROPOMI. Here we present a blended TROPOMI+GOSAT methane product obtained by training a machine learning (ML) model to predict the difference between TROPOMI and GOSAT co-located measurements, using only predictor variables included in the TROPOMI retrieval, and then applying the correction to the complete TROPOMI record from January 2018 to present. We find that the largest corrections are associated with coarse aerosol particles, high SWIR surface albedo, and across-track pixel index. Our blended product corrects a systematic difference between TROPOMI and GOSAT over water, and it features corrections exceeding 10 ppb over arid land, persistently cloudy regions, and high northern latitudes. It reduces the TROPOMI spatially variable bias over land (referenced to GOSAT data) from 14.0 to 10.7 ppb at 0.25° × 0.3125° resolution. Validation with TCCON ground-based column measurements shows reductions in variable bias compared to the original TROPOMI data from 6.0 to 5.2 ppb and in single-retrieval precision from 13.8 to 11.7 ppb. TCCON data are all in locations of SWIR surface albedo below 0.4 (where TROPOMI biases tend to be relatively low), but they confirm the dependence of TROPOMI biases on SWIR surface albedo and coarse aerosol particles, as well as the reduction of these biases in the blended product. Fine-scale inspection of the Arabian Peninsula shows that a number of hotspots in the original TROPOMI data are removed as artifacts in the blended product. The blended product also corrects striping and aerosol/cloud biases in single-orbit TROPOMI data, enabling better detection and quantification of ultra-emitters. Residual coastal biases can be removed by applying additional filters. The ML method presented here can be applied more generally to validate and correct data from any new satellite instrument by reference to a more established instrument.

We quantify methane emissions and their 2010–2017 trends by sector in the contiguous United States (CONUS), Canada, and Mexico by inverse analysis of in situ (GLOBALVIEWplus CH4 ObsPack) and satellite (GOSAT) atmospheric methane observations. The inversion uses as a prior estimate the national anthropogenic emission inventories for the three countries reported by the US Environmental Protection Agency (EPA), Environment and Climate Change Canada (ECCC), and the Instituto Nacional de Ecología y Cambio Climático (INECC) in Mexico to the United Nations Framework Convention on Climate Change (UNFCCC) and thus serves as an evaluation of these inventories in terms of their magnitudes and trends. Emissions are optimized with a Gaussian mixture model (GMM) at 0.5∘×0.625∘ resolution and for individual years. Optimization is done analytically using lognormal error forms. This yields closed-form statistics of error covariances and information content on the posterior (optimized) estimates, allows better representation of the high tail of the emission distribution, and enables construction of a large ensemble of inverse solutions using different observations and assumptions. We find that GOSAT and in situ observations are largely consistent and complementary in the optimization of methane emissions for North America. Mean 2010–2017 anthropogenic emissions from our base GOSAT + in situ inversion, with ranges from the inversion ensemble, are 36.9 (32.5–37.8) Tg a−1 for CONUS, 5.3 (3.6–5.7) Tg a−1 for Canada, and 6.0 (4.7–6.1) Tg a−1 for Mexico. These are higher than the most recent reported national inventories of 26.0 Tg a−1 for the US (EPA), 4.0 Tg a−1 for Canada (ECCC), and 5.0 Tg a−1 for Mexico (INECC). The correction in all three countries is largely driven by a factor of 2 underestimate in emissions from the oil sector with major contributions from the south-central US, western Canada, and southeastern Mexico. Total CONUS anthropogenic emissions in our inversion peak in 2014, in contrast to the EPA report of a steady decreasing trend over 2010–2017. This reflects offsetting effects of increasing emissions from the oil and landfill sectors, decreasing emissions from the gas sector, and flat emissions from the livestock and coal sectors. We find decreasing trends in Canadian and Mexican anthropogenic methane emissions over the 2010–2017 period, mainly driven by oil and gas emissions. Our best estimates of mean 2010–2017 wetland emissions are 8.4 (6.4–10.6) Tg a−1 for CONUS, 9.9 (7.8–12.0) Tg a−1 for Canada, and 0.6 (0.4–0.6) Tg a−1 for Mexico. Wetland emissions in CONUS show an increasing trend of +2.6 (+1.7 to +3.8)% a−1 over 2010–2017 correlated with precipitation.

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).

We use satellite methane observations from the Tropospheric Monitoring Instrument (TROPOMI), for May 2018 to February 2020, to quantify methane emissions from individual oil and natural gas (O/G) basins in the US and Canada using a high-resolution (∼25 km) atmospheric inverse analysis. Our satellite-derived emission estimates show good consistency with in situ field measurements (R=0.96) in 14 O/G basins distributed across the US and Canada. Aggregating our results to the national scale, we obtain O/G-related methane emission estimates of 12.6±2.1 Tg a−1 for the US and 2.2±0.6 Tg a−1 for Canada, 80 % and 40 %, respectively, higher than the national inventories reported to the United Nations. About 70 % of the discrepancy in the US Environmental Protection Agency (EPA) inventory can be attributed to five O/G basins, the Permian, Haynesville, Anadarko, Eagle Ford, and Barnett basins, which in total account for 40 % of US emissions. We show more generally that our TROPOMI inversion framework can quantify methane emissions exceeding 0.2–0.5 Tg a−1 from individual O/G basins, thus providing an effective tool for monitoring methane emissions from large O/G basins globally.

We present a user-friendly, cloud-based facility for quantifying methane emissions with 0.25∘ × 0.3125∘ (≈ 25 km × 25 km) resolution by inverse analysis of satellite observations from the TROPOspheric Monitoring Instrument (TROPOMI). The facility is built on an Integrated Methane Inversion optimal estimation workflow (IMI 1.0) and supported for use on the Amazon Web Services (AWS) cloud. It exploits the GEOS-Chem chemical transport model and TROPOMI data already resident on AWS, thus avoiding cumbersome big-data download. Users select a region and period of interest, and the IMI returns an analytical solution for the Bayesian optimal estimate of period-average emissions on the 0.25∘ × 0.3125∘ grid including error statistics, information content, and visualization code for inspection of results. The inversion uses an advanced research-grade algorithm fully documented in the literature. An out-of-the-box inversion with rectilinear grid and default prior emission estimates can be conducted with no significant learning curve. Users can also configure their inversions to infer emissions for irregular regions of interest, swap in their own prior emission inventories, and modify inversion parameters. Inversion ensembles can be generated at minimal additional cost once the Jacobian matrix for the analytical inversion has been constructed. A preview feature allows users to determine the TROPOMI information content for their region and time period of interest before actually performing the inversion. The IMI is heavily documented and is intended to be accessible by researchers and stakeholders with no expertise in inverse modelling or high-performance computing. We demonstrate the IMI's capabilities by applying it to estimate methane emissions from the US oil-producing Permian Basin in May 2018.

We quantify methane emissions and their 2010–2017 trends by sector in the contiguous United States (CONUS), Canada, and Mexico by inverse analysis of in situ (GLOBALVIEWplus CH4 ObsPack) and satellite (GOSAT) atmospheric methane observations. The inversion uses as prior estimate the national anthropogenic emission inventories for the three countries reported by the US Environmental Protection Agency (EPA), Environment and Climate Change Canada (ECCC), and the Instituto Nacional de Ecologia y Cambio Climatico (INECC) in Mexico to the United Nations Framework Convention on Climate Change (UNFCCC), and thus serves as an evaluation of these inventories in terms of their magnitudes and trends. Emissions are optimized with a Gaussian mixture model (GMM) at 0.5° × 0.625° resolution and for individual years. Optimization is done analytically using log-normal error forms. This yields closed-form statistics of error estimates and information content on the posterior (optimized) estimates, allows better representation of the high tail of the emission distribution, and enables construction of a large ensemble of inverse solutions using different observations and assumptions. We find that GOSAT and in situ observations are largely consistent and complementary in the optimization of methane emissions for North America. Mean 2010–2017 anthropogenic emissions from our base GOSAT + in situ inversion, with ranges from the inversion ensemble, are 36.9 (32.5–37.8) Tg a−1 for CONUS, 5.3 (3.6–5.7) Tg a−1 for Canada, and 6.0 (4.7–6.1) Tg a−1 for Mexico. These are higher than the most recent reported national inventories of 26.0 Tg a−1 for the US (EPA), 4.0 Tg a−1 for Canada (ECCC), and 5.0 Tg a−1 for Mexico (INECC). The correction in all three countries is largely driven by a factor of 2 underestimate in emissions from the oil sector with major contributions from the south-central US, western Canada, and southeast Mexico. Total CONUS anthropogenic emissions in our inversion peak in 2014, in contrast to the EPA report of a steady decreasing trend over 2010–2017. This reflects combined effects of increases in emissions from the oil and landfill sectors, decrease from the gas, and flat emissions from the livestock and coal sectors. We find decreasing trends in Canadian and Mexican anthropogenic methane emissions over the 2010–2017 period, mainly driven by oil and gas emissions. Our best estimates of mean 2010–2017 wetland emissions are 8.4 (6.4–10.6) Tg a−1 for CONUS, 9.9 (7.8–12.0) Tg a−1 for Canada, and 0.6 (0.4–0.6) Tg a−1 for Mexico. Wetland emissions in CONUS show an increasing trend of 2.6 (1.7–3.8) % a−1 over 2010–2017 correlated with precipitation.

We use 2010–2015 GOSAT satellite observations of atmospheric methane columns over North America in a high- resolution inversion of methane emissions, including contributions from different sectors and long-term trends. The inversion involves 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. 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. Prior estimates for the inversion include a gridded version of the EPA Inventory of U.S. 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 a factor 2 lower than our estimate. 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) and calling 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/gas production in the eastern US. We also find that oil/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.

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 datasets. 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 mid-latitudes and for relaxing global error correlations between methane emissions and the main methane sink (oxidation by OH radicals). The GOSAT observations achieve 212 independent pieces of information (DOFS) for quantifying mean 2010–2017 anthropogenic emissions on 1009 global model grid elements, and a DOFS of 122 for 2010–2017 emission trends. Adding the in situ data increases the DOFS by about 20–30 %, to 262 and 161 respectively for mean emissions and trends. Our joint inversion finds that oil/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, while coal emissions in China are overestimated. Wetland emissions in North America are much lower than in the mean WetCHARTs inventory used as prior estimate. Oil/gas emissions in the US increase over the 2010–2017 period but decrease in Canada and Europe. Our joint GOSAT+in situ 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).

Satellite observes large methane emissions from the Permian Basin in the United States., 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.

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 log-normal. Inversion results show large overestimates of Chinese coal emissions and Middle East oil/gas emissions in the EDGAR v4.3.2 inventory, but little error in the US where we use a new gridded version of the EPA national greenhouse gas inventory as prior estimate. Oil/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.