Search

Nationally reported greenhouse gas inventories are a core component of the Paris Agreements transparency framework. Comparisons with emission estimates derived from atmospheric observations help identify improvements to reduce uncertainties and increase the confidence in reported values. To facilitate comparisons over the contiguous United States, we present a 0.1° × 0.1° gridded inventory of annual 2012-2018 anthropogenic methane emissions, allocated to 26 individual source categories, with scale-dependent error estimates. Our inventory is consistent with the U.S. Environmental Protection Agency (EPA) Inventory of U.S. Greenhouse Gas Emissions and Sinks (GHGI), submitted to the United Nations in 2020. Total emissions and patterns (spatial/temporal) reflect the activity and emission factor data underlying the GHGI, including many updates relative to a previous gridded version of the GHGI that has been extensively compared with observations. These underlying data are not generally available in global gridded inventories, and comparison to EDGAR version 6 shows large spatial differences, particularly for the oil and gas sectors. We also find strong regional variability across all sources in annual 2012-2018 spatial trends, highlighting the importance of understanding regional- and facility-level activities. Our inventory represents the first time series of gridded GHGI methane emissions and enables robust comparisons of emissions and their trends with atmospheric observations.

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.

Global modeling of aerosol‐particle number and size is important for understanding aerosol effects on Earth's climate and air quality. Fine‐resolution global models are desirable for representing nonlinear aerosol‐microphysical processes, their nonlinear interactions with dynamics and chemistry, and spatial heterogeneity. However, aerosol‐microphysical simulations are computationally demanding, which can limit the achievable global horizontal resolution. Here, we present the first coupling of the TwO‐Moment Aerosol Sectional (TOMAS) microphysics scheme with the High‐Performance configuration of the GEOS‐Chem model of atmospheric composition (GCHP), a coupling termed GCHP‐TOMAS. GCHP's architecture allows massively parallel GCHP‐TOMAS simulations including on the cloud, using hundreds of computing cores, faster runtimes, more memory, and finer global horizontal resolution (e.g., 25 km × 25 km, 7.8 × 105 model columns) versus the previous single‐node capability of GEOS‐Chem‐TOMAS (tens of cores, 200 km × 250 km, 1.3 × 104 model columns). GCHP‐TOMAS runtimes have near‐ideal scalability with computing‐core number. Simulated global‐mean number concentrations increase (dominated by free‐tropospheric over‐ocean sub‐10‐nm‐diameter particles) toward finer GCHP‐TOMAS horizontal resolution. Increasing the horizontal resolution from 200 km × 200–50 km × 50 km increases the global monthly mean free‐tropospheric total particle number by 18.5%, and over‐ocean sub‐10‐nm‐diameter particles by 39.8% at 4‐km altitude. With a cascade of contributing factors, free‐tropospheric particle‐precursor concentrations increase (32.6% at 4‐km altitude) with resolution, promoting new‐particle formation and growth that outweigh coagulation changes. These nonlinear effects have the potential to revise current understanding of processes controlling global aerosol number and aerosol impacts on Earth's climate and air quality. Plain Language Summary: Small particles in the air have important effects on Earth's climate and air quality. Representing the number and size of these particles in global models is challenging because their processes are complex. This factor has often limited global‐model horizontal resolution because fine global resolution models (e.g., 25 km × 25 km or smaller) generally ran too slowly but would be useful for representing details missed at traditional coarse resolution (e.g., 200 km × 250 km). We start with a detailed particle scheme that previously only ran at coarse global resolution because fine resolution would take too long. We present the initial use of this scheme in an updated model version, with a structure allowing a fast‐running, high‐memory model with fine resolution, by using hundreds to thousands of computer cores. In the updated structure, model speed increases with the number of cores used. We find that the total number of particles in the model is more with fine compared to coarse model resolution. These increases are most in Earth's remote regions and for particles which come from gas. Using fine model resolution globally when representing particles could change our understanding of how they impact Earth's climate and air quality. Key Points: We couple aerosol microphysics with GEOS‐Chem's High‐Performance configuration for fine (25 km × 25 km) global‐resolution capabilityGlobal‐mean aerosol number increases with model resolution, dominated by particles smaller than 10 nm in the over‐ocean free troposphereToward finer horizontal resolution, enhanced particle precursor loading in the free troposphere promotes particle formation and growth [ABSTRACT FROM AUTHOR]

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]

Ozone pollution in the Southeast US involves complex chemistry driven by emissions of anthropogenic nitrogen oxide radicals (NOx ≡ NO + NO2) and biogenic isoprene. Model estimates of surface ozone concentrations tend to be biased high in the region and this is of concern for designing effective emission control strategies to meet air quality standards. We use detailed chemical observations from the SEAC4RS aircraft campaign in August and September 2013, interpreted with the GEOS-Chem chemical transport model at 0.25°×0.3125° horizontal resolution, to better understand the factors controlling surface ozone in the Southeast US. We find that the National Emission Inventory (NEI) for NOx from the US Environmental Protection Agency (EPA) is too high. This finding is based on SEAC4RS observations of NOx and its oxidation products, surface network observations of nitrate wet deposition fluxes, and OMI satellite observations of tropospheric NO2 columns. Our results indicate that NEI NOx emissions from mobile and industrial sources must be reduced by 30-60%, dependent on the assumption of the contribution by soil NOx emissions. Upper tropospheric NO2 from lightning makes a large contribution to satellite observations of tropospheric NO2 that must be accounted for when using these data to estimate surface NOx emissions. We find that only half of isoprene oxidation proceeds by the high-NOx pathway to produce ozone; this fraction is only moderately sensitive to changes in NOx emissions because isoprene and NOx emissions are spatially segregated. GEOS-Chem with reduced NOx emissions provides an unbiased simulation of ozone observations from the aircraft, and reproduces the observed ozone production efficiency in the boundary layer as derived from a regression of ozone and NOx oxidation products. However, the model is still biased high by 8±13 ppb relative to observed surface ozone in the Southeast US. Ozonesondes launched during midday hours show a 7 ppb ozone decrease from 1.5 km to the surface that GEOS-Chem does not capture. This bias may reflect a combination of excessive vertical mixing and net ozone production in the model boundary layer.

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 theWeather 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]

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 stateproduced 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]

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]

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]