Search

Global high-resolution observations of atmospheric composition from satellites can greatly improve our understanding of surface emissions through inverse analyses. Variational inverse methods can optimize surface emissions at any resolution but do not readily quantify the error and information content of the posterior solution. The information content of satellite data may be much lower than its coverage would suggest because of failed retrievals, instrument noise, and error correlations that propagate through the inversion. Analytical solution of the inverse problem provides closed-form characterization of posterior error statistics and information content but requires the construction of the Jacobian matrix that relates emissions to atmospheric concentrations. Building the Jacobian matrix is computationally expensive at high resolution because it involves perturbing each emission element, typically individual grid cells, in the atmospheric transport model used as the forward model for the inversion. We propose and analyze two methods, reduced dimension and reduced rank, to construct the Jacobian matrix at greatly decreased computational cost while retaining information content. Both methods are two-step iterative procedures that begin from an initial native-resolution estimate of the Jacobian matrix constructed at no computational cost by assuming that atmospheric concentrations are most sensitive to local emissions. The reduced-dimension method uses this estimate to construct a Jacobian matrix on a multiscale grid that maintains a high resolution in areas with high information content and aggregates grid cells elsewhere. The reduced-rank method constructs the Jacobian matrix at native resolution by perturbing the leading patterns of information content given by the initial estimate. We demonstrate both methods in an analytical Bayesian inversion of Greenhouse Gases Observing Satellite (GOSAT) methane data with augmented information content over North America in July 2009. We show that both methods reproduce the results of the native-resolution inversion while achieving a factor of 4 improvement in computational performance. The reduced-dimension method produces an exact solution at a lower spatial resolution, while the reduced-rank method solves the inversion at native resolution in areas of high information content and defaults to the prior estimate elsewhere.

We quantify 2019 annual mean methane emissions in the contiguous US (CONUS) at 0.25° × 0.3125° resolution by inverse analysis of atmospheric methane columns measured by the Tropospheric Monitoring Instrument (TROPOMI). A gridded version of the US Environmental Protection Agency (EPA) Greenhouse Gas Emissions Inventory (GHGI) serves as the basis for the prior estimate for the inversion. We optimize emissions and quantify observing system information content for an eight-member inversion ensemble through analytical minimization of a Bayesian cost function. We achieve high resolution with a reduced-rank characterization of the observing system that optimally preserves information content. Our optimal (posterior) estimate of anthropogenic emissions in CONUS is 30.9 (30.0–31.8) Tg a -1 , where the values in parentheses give the spread of the ensemble. This is a 13 % increase from the 2023 GHGI estimate for CONUS in 2019. We find emissions for livestock of 10.4 (10.0–10.7) Tg a -1 , for oil and gas of 10.4 (10.1–10.7) Tg a -1 , for coal of 1.5 (1.2–1.9) Tg a -1 , for landfills of 6.9 (6.4–7.5) Tg a -1 , for wastewater of 0.6 (0.5–0.7), and for other anthropogenic sources of 1.1 (1.0–1.2) Tg a -1. The largest increase relative to the GHGI occurs for landfills (51 %), with smaller increases for oil and gas (12 %) and livestock (11 %). These three sectors are responsible for 89 % of posterior anthropogenic emissions in CONUS. The largest decrease (28 %) is for coal. We exploit the high resolution of our inversion to quantify emissions from 70 individual landfills, where we find emissions are on median 77 % larger than the values reported to the EPA's Greenhouse Gas Reporting Program (GHGRP), a key data source for the GHGI. We attribute this underestimate to overestimated recovery efficiencies at landfill gas facilities and to under-accounting of site-specific operational changes and leaks. We also quantify emissions for the 48 individual states in CONUS, which we compare to the GHGI's new state-level inventories and to independent state-produced inventories. Our posterior emissions are on average 27 % larger than the GHGI in the largest 10 methane-producing states, with the biggest upward adjustments in states with large oil and gas emissions, including Texas, New Mexico, Louisiana, and Oklahoma. We also calculate emissions for 95 geographically diverse urban areas in CONUS. Emissions for these urban areas total 6.0 (5.4–6.7) Tg a -1 and are on average 39 (27–52) % larger than a gridded version of the 2023 GHGI, which we attribute to underestimated landfill and gas distribution emissions. [ABSTRACT FROM AUTHOR]

Satellite observations of dry-column methane mixing ratios (XCH 4) from shortwave infrared (SWIR) solar backscatter radiation provide a powerful resource to quantify methane emissions in service of climate action. The TROPOspheric Monitoring Instrument (TROPOMI), launched in October 2017, provides global daily coverage at a 5.5 × 7 km 2 (nadir) pixel resolution, but its methane retrievals can suffer from biases associated with SWIR surface albedo, scattering from aerosols and cirrus clouds, and across-track variability (striping). The Greenhouse gases Observing SATellite (GOSAT) instrument, launched in 2009, has better spectral characteristics and its methane retrieval is much less subject to biases, but its data density is 250 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 April 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.3 to 10.4 ppb at a 0.25 ∘ × 0.3125 ∘ resolution. Validation with Total Carbon Column Observing Network (TCCON) ground-based column measurements shows reductions in variable bias compared with the original TROPOMI data from 4.7 to 4.4 ppb and in single-retrieval precision from 14.5 to 11.9 ppb. TCCON data are all in locations with a 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. [ABSTRACT FROM AUTHOR]

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 nonlinear 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 Harmonized Emissions Component (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 post-processing 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 ∘ × 2.5 ∘ spatial resolution for 2019 using TROPOspheric Monitoring Instrument (TROPOMI) satellite observations. CHEEREIO achieves a 50-fold improvement in computational performance compared to the equivalent analytical inversion of TROPOMI observations. [ABSTRACT FROM AUTHOR]

We quantify weekly methane emissions at 0.25 ∘ × 0.3125 ∘ (≈25 × 25 km 2) 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. [ABSTRACT FROM AUTHOR]

We quantify 2019 methane emissions in the contiguous U.S. (CONUS) at 0.25° × 0.3125° resolution by inverse analysis of atmospheric methane columns measured by the Tropospheric Monitoring Instrument (TROPOMI). A gridded version of the U.S. Environmental Protection Agency (EPA) Greenhouse Gas Emissions Inventory (GHGI) serves as the basis for the prior estimate for the inversion. We optimize emissions and quantify observing system information content for an eight-member inversion ensemble through analytical minimization of a Bayesian cost function. We achieve high resolution with a reduced-rank characterization of the observing system that optimally preserves information content. Our optimal (posterior) estimate of anthropogenic emissions in CONUS is 30.9 (30.0–31.8) Tg a-1, where the values in parentheses give the spread of the ensemble. This is a 13 % increase from the 2023 GHGI estimate for CONUS in 2019. We find livestock emissions of 10.4 (10.0–10.7) Tg a-1, oil and gas of 10.4 (10.1–10.7) Tg a-1, coal of 1.5 (1.2–1.9) Tg a-1, landfills of 6.9 (6.4–7.5) Tg a-1, wastewater of 0.6 (0.5–0.7), and other anthropogenic sources of 1.1 (1.0–1.2) Tg a-1. The largest increase relative to the GHGI occurs for landfills (51 %), with smaller increases for oil and gas (12 %) and livestock (11 %). These three sectors are responsible for 89 % of posterior anthropogenic emissions in CONUS. The largest decrease (28 %) is for coal. We exploit the high resolution of our inversion to quantify emissions from 73 individual landfills, where we find emissions are on median 77 % larger than the values reported to the EPA's Greenhouse Gas Reporting Program (GHGRP), a key data source for the GHGI. We attribute this underestimate to overestimated recovery efficiencies at landfill gas facilities and to under-accounting of site-specific operational changes and leaks. We also quantify emissions for the 48 individual states in CONUS, which we compare to the GHGI's new state-level inventories and to independent state-produced inventories. Our posterior emissions are on average 34 % larger than the 2022 GHGI in the largest 10 methane-producing states, with the biggest upward adjustments in states with large oil and gas emissions, including Texas, New Mexico, Louisiana, and Oklahoma. We also calculate emissions for 95 geographically diverse urban areas in CONUS. Emissions for these urban areas total 6.0 (5.4–6.7) Tg a-1 and are on average 39 (27–52) % larger than a gridded version of the 2023 GHGI, which we attribute to underestimated landfill and gas distribution emissions. [ABSTRACT FROM AUTHOR]

We use 2019 TROPOMI satellite observations of atmospheric methane in an analytical inversion to quantify methane emissions from the Middle East and North Africa at up to ∼25 km × 25 km resolution, using spatially allocated national United Nations Framework Convention on Climate Change (UNFCCC) reports as prior estimates for the fuel sector. Our resulting best estimate of anthropogenic emissions for the region is 35 % higher than the prior bottom-up inventories (+ 103 % for gas, + 53 % for waste, + 49 % for livestock, -14 % for oil) with large variability across countries. Oil and gas account for 38 % of total anthropogenic emissions in the region. TROPOMI observations can effectively optimize and separate national emissions by sector for most of the 23 countries in the region, with 6 countries accounting for most of total anthropogenic emissions including Iran (5.3 (5.0–5.5) Tg a -1 ; best estimate and uncertainty range), Turkmenistan (4.4 (2.8–5.1) Tg a -1) , Saudi Arabia (4.3 (2.4–6.0) Tg a -1) , Algeria (3.5 (2.4–4.4) Tg a -1) , Egypt (3.4 (2.5–4.0) Tg a -1) , and Turkey (3.0 (2.0–4.1) Tg a -1). Most oil–gas emissions are from the production (upstream) subsector, but Iran, Turkmenistan, and Saudi Arabia have large gas emissions from transmission and distribution subsectors. We identify a high number of annual oil–gas emission hotspots in Turkmenistan, Algeria, and Oman and offshore in the Persian Gulf. We show that oil–gas methane emissions for individual countries are not related to production, invalidating a basic premise in the construction of activity-based bottom-up inventories. Instead, local infrastructure and management practices appear to be key drivers of oil–gas emissions, emphasizing the need for including top-down information from atmospheric observations in the construction of oil–gas emission inventories. We examined the methane intensity, defined as the upstream oil–gas emission per unit of methane gas produced, as a measure of the potential for decreasing emissions from the oil–gas sector and using as reference the 0.2 % target set by the industry. We find that the methane intensity in most countries is considerably higher than this target, reflecting leaky infrastructure combined with deliberate venting or incomplete flaring of gas. However, we also find that Kuwait, Saudi Arabia, and Qatar meet the industry target and thus show that the target is achievable through the capture of associated gas, modern infrastructure, and the concentration of operations. Decreasing methane intensities across the Middle East and North Africa to 0.2 % would achieve a 90 % decrease in oil–gas upstream emissions and a 26 % decrease in total anthropogenic methane emissions in the region, making a significant contribution toward the Global Methane Pledge. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

We use 2019 TROPOMI satellite observations of atmospheric methane in an analytical inversion to quantify methane emissions from the Middle East and North Africa at up to ~25 km × 25 km resolution, using spatially allocated national UNFCCC reports as prior estimates for the fuel sector. Our resulting best estimate of anthropogenic emissions for the region is 35 % higher than the prior bottom-up estimate (+103 % for gas, +53 % for waste, +49 % for livestock, −14 % for oil) with large variability across countries. Oil and gas account for 38 % of total anthropogenic emissions in the region. TROPOMI observations can effectively optimize and separate national emissions by sector for most of the 23 countries in the region, with 6 countries accounting for most of total anthropogenic emissions including Iran (5.3 (5.0–5.5) Tg a−1; best estimate and uncertainty range), Turkmenistan (4.4 (2.8–5.1) Tg a−1), Saudi Arabia (4.3 (2.4–6.0) Tg a−1), Algeria (3.5 (2.4–4.4) Tg a−1), Egypt (3.4 (2.5–4.0) Tg a−1) , and Turkey (3.0 (2.0–4.1) Tg a−1). Most oil/gas emissions are from the production (upstream) subsector, but Iran, Turkmenistan, and Saudi Arabia have large gas emissions from transmission and distribution subsectors. We identify a high number of annual oil/gas emission hotspots in Turkmenistan, Algeria, Oman, and offshore in the Persian Gulf. We show that oil/gas methane emissions for individual countries are not related to production, invalidating a basic premise in the construction of activity-based bottom-up inventories. Instead, local infrastructure and management practices appear to be key drivers of oil/gas emissions, emphasizing the need for including top-down information from atmospheric observations in the construction of oil/gas emission inventories. We examined the methane intensity, defined as the upstream oil/gas emission per unit of methane gas produced, as a measure of the potential for decreasing emissions from the oil/gas sector, and using as reference the 0.2 % target set by industry. We find that the methane intensity in most countries is considerably higher than this target, reflecting leaky infrastructure combined with deliberate venting or incomplete flaring of gas. However, we also find that Kuwait, Saudi Arabia, and Qatar meet the industry target and thus show that the target is achievable through capture of associated gas, modern infrastructure, and concentration of operations. Decreasing methane intensities across the Middle East and North Africa to 0.2 % would achieve a 90 % decrease in oil/gas upstream emissions and a 26 % decrease of total anthropogenic methane emissions in the region, making a significant contribution toward the Global Methane Pledge. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

We quantify weekly methane emissions at 0.25°×0.3125° (≈25×25 km²) resolution from the Permian Basin, the largest oil production basin in the United States, 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 too-low prior emissions or biased 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 local 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 over the period 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 recent industry targets 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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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 a prior estimate the latest 2014 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 log-normally distributed errors on prior emissions. Errors and information content on the optimized estimates are obtained directly from the analytical solution and also through a 36-member inversion ensemble. Our best estimate for total anthropogenic emissions in China is 65.0 (57.7–68.4) Tg a -1 , where parentheses indicate the uncertainty range determined by the inversion ensemble. 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-for 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. [ABSTRACT FROM AUTHOR]

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.25o×0.3125o 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. [ABSTRACT FROM AUTHOR]

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) Tga-1 for CONUS, 5.3 (3.6–5.7) Tga-1 for Canada, and 6.0 (4.7–6.1) Tga-1 for Mexico. These are higher than the most recent reported national inventories of 26.0 Tga-1 for the US (EPA), 4.0 Tga-1 for Canada (ECCC), and 5.0 Tga-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) Tga-1 for CONUS, 9.9 (7.8–12.0) Tga-1 for Canada, and 0.6 (0.4–0.6) Tga-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. [ABSTRACT FROM AUTHOR]

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 FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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° x 2.5° horizontal resolution. We compare the results to an inversion using the sparser but more mature GOSAT satellite retrievals, as well as 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 optimization 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° x 2.5° scale of the inversion and large spatial variations of the number of observations make it less useful than GOSAT for quantifying emissions at that resolution. 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 FROM AUTHOR]

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 FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Global high-resolution observations of atmospheric composition from satellites can greatly improve our understanding of surface emissions through inverse analyses. Variational inverse methods can optimize surface emissions at any resolution but do not readily quantify the error and information content of the posterior solution. In fact, the information content of satellite data may be orders of magnitude lower than its coverage suggests because of failed retrievals, instrument noise, and error correlations that propagate through the inversion. Analytical solution to the inverse problem provides closed-form characterization of posterior error statistics and information content but requires the construction of the Jacobian matrix that relates emissions to atmospheric concentrations. Building the Jacobian matrix is computationally expensive at high resolution because it involves perturbing each emission element, typically individual grid cells, in the atmospheric transport model used as forward model for the inversion. We propose and analyze two methods, reduced-dimension and reduced-rank, to construct the Jacobian matrix at greatly decreased computational cost while retaining information content. Both methods begin from an initial native-resolution estimate of the Jacobian matrix constructed at no computational cost by assuming that atmospheric concentrations are most sensitive to local emissions. The reduced-dimension method uses this estimate to construct a Jacobian matrix on a multiscale grid that maintains high resolution in areas with high information content and aggregates grid cells elsewhere. The reduced-rank method constructs the Jacobian matrix at native resolution by perturbing the leading patterns of information content given by the initial estimate. We demonstrate both methods in an analytical Bayesian inversion of GOSAT methane satellite data with augmented information content over North America in July 2009. We show that both methods reproduce the results of the native-resolution inversion while achieving a factor of 4 improvement in computational performance. The reduced-dimension method produces an exact solution at lower spatial resolution while the reduced-rank method solves the inversion at native resolution in areas of high information content and defaults to the prior estimate elsewhere. [ABSTRACT FROM AUTHOR]

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). [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

The article discusses that Methane is a potent greenhouse gas with a relatively short average atmospheric residence time of about a decade and is also a precursor of tropospheric ozone; and mentions the rapid increase in oil and natural gas (O/G) production in the U.S. since around 2005, driven primarily by hydraulic fracturing and horizontal drilling, has led to major concerns about increasing methane emissions and adverse climate impacts.

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., (Copyright © 2020 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works. Distributed under a Creative Commons Attribution NonCommercial License 4.0 (CC BY-NC).)