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.

Southeast Australia experienced intensive and geographically extensive wildfires during the 2019-2020 summer season.sup.1,2. The fires released substantial amounts of carbon dioxide into the atmosphere.sup.3. However, existing emission estimates based on fire inventories are uncertain.sup.4, and vary by up to a factor of four for this event. Here we constrain emission estimates with the help of satellite observations of carbon monoxide.sup.5, an analytical Bayesian inversion.sup.6 and observed ratios between emitted carbon dioxide and carbon monoxide.sup.7. We estimate emissions of carbon dioxide to be 715 teragrams (range 517-867) from November 2019 to January 2020. This is more than twice the estimate derived by five different fire inventories.sup.8-12, and broadly consistent with estimates based on a bottom-up bootstrap analysis of this fire episode.sup.13. Although fires occur regularly in the savannas in northern Australia, the recent episodes were extremely large in scale and intensity, burning unusually large areas of eucalyptus forest in the southeast.sup.13. The fires were driven partly by climate change.sup.14,15, making better-constrained emission estimates particularly important. This is because the build-up of atmospheric carbon dioxide may become increasingly dependent on fire-driven climate-carbon feedbacks, as highlighted by this event.sup.16. The amount of carbon dioxide released by the Australian wildfires of 2019-2020 is uncertain, but is estimated here using satellite observations of carbon monoxide to be more than twice the amount suggested by fire inventories., Author(s): Ivar R. van der Velde [sup.1] [sup.2] , Guido R. van der Werf [sup.2] , Sander Houweling [sup.1] [sup.2] , Joannes D. Maasakkers [sup.1] , Tobias Borsdorff [sup.1] , [...]

Venezuela has long been identified as an area with large methane emissions and intensive oil exploitation, especially in the Lake Maracaibo region, but production has strongly decreased in recent years. The area is notoriously difficult to observe from space due to its complex topography and persistent cloud cover. We use the unprecedented coverage of the TROPOspheric Monitoring Instrument (TROPOMI) methane observations in analytical inversions with the Integrated Methane Inversion (IMI) framework at the national scale and at the local scale with the Weather Research and Forecasting model with chemistry (WRF-Chem). In the IMI analysis, we find Venezuelan emissions of 7.5 (5.7–9.3) Tga-1 in 2019, where about half of emissions can be informed by TROPOMI observations, and emissions from oil exploitation are a factor of ∼ 1.6 higher than in bottom-up inventories. Using WRF, we find emissions of 1.2 (1.0–1.5) Tga-1 from the Lake Maracaibo area in 2019, close to bottom-up estimates. Our WRF estimate is ∼ 40 % lower than the result over the same region from the IMI due to differences in the meteorology used by the two models. We find only a small, non-significant trend in emissions between 2018 and 2020 around the lake, implying the area's methane emission intensity expressed against oil and gas production has doubled over the time period, to ∼ 20 %. This value is much higher than what has previously been found for other oil and gas production regions and indicates that there could be large emissions from abandoned infrastructure. [ABSTRACT FROM AUTHOR]

Venezuela has long been identified as an area with large methane emissions, especially in the Lake Maracaibo region with intensive oil exploitation, but production has strongly decreased in recent years. The area is notoriously difficult to observe from space due to its complex topography and persistent cloud cover. We use the unprecedented coverage of the TROPOspheric Monitoring Instrument (TROPOMI) methane observations in analytical inversions at the national scale with the Integrated Methane Inversion (IMI) framework and at regional scale with the Weather Research and Forecasting (WRF) model. In the IMI analysis, we find 2019 Venezuelan emissions of 7.5 (5.7–9.3) Tg a-1, where about half of emissions can be informed by TROPOMI observations and emissions from oil exploitation are a factor ~1.6 higher than in bottom-up inventories. Using WRF, we find 2019 emissions of 1.2 (1.0–1.5) Tg a-1 from the Lake Maracaibo area, close to bottom-up estimates. Our WRF estimate is ~40 % lower than the regional result from the IMI due to differences in the meteorology used by the two models. We only find a small, non-significant, trend in emissions between 2018 and 2020 around the lake, implying the area's methane emission intensity expressed against oil/gas production has doubled over the time period to ~20 %. This value is much higher than what has previously been found for other oil/gas production regions and indicates there could be large emissions from abandoned infrastructure. [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]

Carbon monoxide (CO) is an air pollutant that plays an important role in atmospheric chemistry and is mostly emitted by forest fires and incomplete combustion in, for example, road transport, residential heating, and industry. As CO is co-emitted with fossil fuel CO2 combustion emissions, it can be used as a proxy for CO2. Following the Paris Agreement, there is a need for independent verification of reported activity-based bottom-up CO2 emissions through atmospheric measurements. CO can be observed daily at a global scale with the TROPOspheric Monitoring Instrument (TROPOMI) satellite instrument with daily global coverage at a resolution down to 5.5 × 7 km 2. To take advantage of this unique TROPOMI dataset, we develop a cross-sectional flux-based emission quantification method that can be applied to quantify emissions from a large number of cities, without relying on computationally expensive inversions. We focus on Africa as a region with quickly growing cities and large uncertainties in current emission estimates. We use a full year of high-resolution Weather Research and Forecasting (WRF) simulations over three cities to evaluate and optimize the performance of our cross-sectional flux emission quantification method and show its reliability down to emission rates of 0.1 Tg CO yr -1. Comparison of the TROPOMI-based emission estimates to the Dynamics–Aerosol–Chemistry–Cloud Interactions in West Africa (DACCIWA) and Emissions Database for Global Atmospheric Research (EDGAR) bottom-up inventories shows that CO emission rates in northern Africa are underestimated in EDGAR, suggesting overestimated combustion efficiencies. We see the opposite when comparing TROPOMI to the DACCIWA inventory in South Africa and Côte d'Ivoire, where CO emission factors appear to be overestimated. Over Lagos and Kano (Nigeria) we find that potential errors in the spatial disaggregation of national emissions cause errors in DACCIWA and EDGAR respectively. Finally, we show that our computationally efficient quantification method combined with the daily TROPOMI observations can identify a weekend effect in the road-transport-dominated CO emissions from Cairo and Algiers. [ABSTRACT FROM AUTHOR]

We apply atmospheric methane column retrievals from two different satellite instruments (Greenhouse gases Observing SATellite – GOSAT; TROPOspheric Monitoring Instrument – TROPOMI) to a regional inversion framework to quantify East Asian methane emissions for 2019 at 0.5 ∘ × 0.625 ∘ horizontal resolution. The goal is to assess if GOSAT (relatively mature but sparse) and TROPOMI (new and dense) observations inform consistent methane emissions from East Asia with identically configured inversions. Comparison of the results from the two inversions shows similar correction patterns to the prior inventory in central northern China, central southern China, northeastern China, and Bangladesh, with less than 2.6 Tg a -1 differences in regional posterior emissions. The two inversions, however, disagree over some important regions, particularly in northern India and eastern China. The methane emissions inferred from GOSAT observations are 7.7 Tg a -1 higher than those from TROPOMI observations over northern India but 6.4 Tg a -1 lower over eastern China. The discrepancies between the two inversions are robust against varied inversion configurations (i.e., assimilation window and error specifications). We find that the lower methane emissions from eastern China inferred by the GOSAT inversion are more consistent with independent ground-based in situ and total column (TCCON) observations, indicating that the TROPOMI retrievals may have high XCH 4 biases in this region. We also evaluate inversion results against tropospheric aircraft observations over India during 2012–2014 by using a consistent GOSAT inversion of earlier years as an intercomparison platform. This indirect evaluation favors lower methane emissions from northern India inferred by the TROPOMI inversion. We find that in this case the discrepancy in emission inference is contributed by differences in data coverage (almost no observations by GOSAT vs. good spatial coverage by TROPOMI) over the Indo-Gangetic Plain. The two inversions also differ substantially in their posterior estimates for northwestern China and neighboring Kazakhstan, which is mainly due to seasonally varying biases between GOSAT and TROPOMI XCH 4 data that correlate with changes in surface albedo. [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]

The United States is the world's largest oil/gas methane emitter according to current national reports. Reducing these emissions is a top priority in the US government's climate action plan. Here, we use a 2010 to 2019 high-resolution inversion of surface and satellite observations of atmospheric methane to quantify emission trends for individual oil/gas production regions in North America and relate them to production and infrastructure. We estimate a mean US oil/gas methane emission of 14.8 (12.4 to 16.5) Tg a-1 for 2010 to 2019, 70% higher than reported by the US Environmental Protection Agency. While emissions in Canada and Mexico decreased over the period, US emissions increased from 2010 to 2014, decreased until 2017, and rose again afterward. Increases were driven by the largest production regions (Permian, Anadarko, Marcellus), while emissions in the smaller production regions generally decreased. Much of the year-to-year emission variability can be explained by oil/gas production rates, active well counts, and new wells drilled, with the 2014 to 2017 decrease driven by reduction in new wells and the 2017 to 2019 surge driven by upswing of production. We find a steady decrease in the oil/gas methane intensity (emission per unit methane gas production) for almost all major US production regions. The mean US methane intensity decreased from 3.7% in 2010 to 2.5% in 2019. If the methane intensity for the oil/gas supply chain continues to decrease at this pace, we may expect a 32% decrease in US oil/gas emissions by 2030 despite projected increases in production. [ABSTRACT FROM AUTHOR]