Open Access Article
Sergio
Ibarra-Espinosa†
*ab,
Lei
Hu
b,
Colin
Harkins
ac,
Brian C.
McDonald
c,
Scot M.
Miller
d,
Youmi
Oh
ab,
Lori
Bruhwiler
b,
Colm
Sweeney
b and
Arlyn
Andrews
b
aCooperative Institute for Research in Environmental Sciences, University of Colorado-Boulder, Boulder, CO, USA. E-mail: sibarrae@umd.edu
bNOAA Global Monitoring Laboratory, Boulder, CO, USA
cNOAA Chemical Sciences Laboratory, Boulder, CO, USA
dDepartment of Environmental Health and Engineering, Whiting School of Engineering, Johns Hopkins University, Baltimore, MD 21218, USA
First published on 11th May 2026
The coronavirus disease 2019 (COVID-19) pandemic disrupted normal human activities worldwide, and mobility reductions resulted in reduced levels of air pollutants and greenhouse gases emissions. Here, we examine the impact of these disruptions on a potent greenhouse gas, methane (CH4), over the U.S. In this study, we quantified CH4 emissions from the contiguous U.S. between 2019 and 2021 by analyzing inverse modeling results derived from atmospheric measurements made at 35 sites across the country. Our estimates indicate emission reductions of −2.5 (±0.43) Tg year−1 CH4 in 2020 and −2.9 (±1.63) Tg year−1 in 2021, relative to 2019. The respective percentage change was a −4.3 (−5.1 to −3.5) % reduction in 2020 and −4.8 (−8.3 to −0.7) % in 2021, relative to 2019. Combining with process-based inventory emission datasets, we found that this reduction was primarily due to decreased fossil fuel and agricultural emissions; however, record-breaking forest fires resulted in an increase of 0.4 (0.1 to 0.8) Tg year−1 in 2020–2019, equal to a 20 (3 to 46) % increase in CH4 emissions from the western U.S.
Environmental significanceThis study leverages the COVID-19 pandemic as a unique “natural experiment” to assess the sensitivity of atmospheric monitoring systems to abrupt changes in anthropogenic activity. By applying high-resolution inverse modeling to atmospheric methane (CH4) observations across the U.S., we detected and quantified subtle emission reductions in 2020 and 2021 that process-based inventories partly underestimated. Specifically, our results reveal that while fossil fuel reductions were consistent with drilling activity declines, agricultural emission reductions—likely driven by disruptions in the livestock sector—were more pronounced than reported in the EPA Greenhouse Gas Inventory. Furthermore, the analysis captures the opposing signal of record-breaking wildfire emissions in the Western U.S. These findings demonstrate the critical role of top–down atmospheric inversions in verifying bottom–up inventories and refining our understanding of how specific sectors drive regional methane budgets, a key step for designing effective climate mitigation policies. |
In this study, we investigated the impact of COVID-19 on U.S. emissions of CH4, a potent greenhouse gas with a ∼11.8 year atmospheric lifetime and a Global Warming Potential of 27 to 30 over 100 years.8 Its current atmospheric abundance contributes a total radiative forcing of 0.54 (0.43 to 0.65) W m−2 relative to 1750.8 On average, between 2008 and 2017, global CH4 emissions have been 577 (550–594) Tg year−1, with the primary sources being agriculture, fossil fuel production and use, and wetlands.9 The main sink of atmospheric CH4 is the reaction with hydroxyl radicals (or OH).10 Atmospheric mole fractions of CH4 increased in the decades before 2000, were relatively stable between 1999 and 2006, and have since increased rapidly.11–13 Some studies have proposed that the recent increases in global atmospheric CH4 abundance are due to intensified emissions from microbial sources such as wetlands and agriculture, in combination with reduced atmospheric mole fractions of OH.14–18
CH4 emission estimates for the U.S. show discrepancies across studies and inventories. For example, the U.S. Environmental Protection Agency Greenhouse Gas Inventory (USEPA-GHGI), a bottom–up inventory based on activity data and emission factors from various source categories, shows a small negative trend in total anthropogenic CH4 emissions, with 33 Tg year−1 in 1990 and 29 Tg year−1 in 2019.19 In 2021, the energy and agriculture sectors were reported to emit 11 Tg year−1 and 10 Tg year−1, or 38% and 36% of total U.S. anthropogenic CH4 emissions, respectively. Other studies using different approaches provide varying estimates. NOAA's CarbonTracker-CH4, an atmosphere-based inversion model that uses CH4 and CH4 isotope observations to constrain sectors, estimated 29 Tg year−1 fossil emissions in 2021 for the contiguous United States (CONUS).20 A fuel-based inventory focusing on direct measurements and modeling specific to the oil and gas sector reported 15 Tg year−1 for oil and gas emissions in the CONUS in 2015.21 For natural sources in the U.S., estimates also differ, ranging from 9.5 Tg year−1, average 2010 to 2017,22 14.3 Tg year−1, average 2008 to 2017,9 to 12.5 Tg year−1, average 2010 to 2019.23 These variations underscore inconsistencies in the current understanding of U.S. CH4 emissions, arising from uncertainties inherent in the different estimation methods.
The marked reduction in anthropogenic activity during the COVID-19 pandemic provided a unique opportunity to investigate its impact on emissions of key pollutants and greenhouse gases, such as CH4. However, this area of research remains a complex area of study, evidenced by discrepancies in quantitative estimates across different analyses and scales. The EPA's GHGI reported that U.S. anthropogenic CH4 emissions decreased by 0.6 Tg year−1 (2.1%) from 2019 to 2020 and by 1.3 Tg year−1 (4.3%) from 2019 to 2021.19 In CarbonTracker-CH4, a near-zero change was found for the U.S. between 2019 and 2020 and a reduction of 2.2 Tg year−1 (−3.5%) was found by comparing 2019 and 2021.20 Nevertheless, both estimates agree with the fact that there was a more substantial reduction in 2021 than in 2020,24 using the satellite-based TROPOspheric Monitoring Instrument (TROPOMI) observations and a 4D-Var Integrated Forecasting System (4D-Var IFS), compared the first six months of 2020 with 2019 and found a global increase of 1.6%, and a 2.2% increase in the U.S. In California, aircraft measurements in 2020 revealed a 35% reduction relative to 2017 in CH4 emissions from the oil and gas sector.25 Here, we analyze U.S. CH4 emissions from 2019 to 2021 derived from a high-resolution regional inverse modeling study.26 We used the notation 2020–2019 and 2021–2019 to represent the anomaly period of the pandemic years, 2020 and 2021, versus the pre-pandemic year 2019. Using these atmospheric observation-based emission estimates, we investigate changes in U.S. national and regional CH4 emissions during and after the COVID-19 pandemic relative to the pre-pandemic year of 2019. Then, we compare these changes with process-based inventories to improve the characterization of changes in CH4 emissions during and after the COVID-19 pandemic.
The presented uncertainties for the annual emissions below include two measures of uncertainty as described in (ref. 29). One is the analytical (internal) uncertainty that quantifies the random errors within a single inversion resulted from prior error covariance and model-data mismatch errors. This quantity reflects how well the atmospheric observations constrain emissions given a specific set of prior assumptions and background estimate. Hence, observational constraints should result in a decrease in uncertainty relative to the prior uncertainty. To characterize the features of our inversion, we define the external uncertainty as one standard deviation across the results of the six ensemble members. This metric captures the sensitivity of the methane flux estimates to variations in prior assumptions, background and transport model configurations. A narrow ensemble spread, therefore, increases confidence that our findings are robust and not an artifact of a single set of modeling choices.
The pandemic anomalies were calculated from the difference in annual CH4 emissions between the pre-pandemic year 2019 and the pandemic years 2020 and 2021 for each of the six ensemble members. Our best estimate of the pandemic CH4 anomalies discussed below were the mean anomaly among the six members and their one standard deviation as its uncertainty. To determine if the emission changes between the pre-pandemic (2019) and pandemic (2020, 2021) periods were statistically significant, we used z-scores for each ensemble member and its associated uncertainty, as shown in Fig. S5. The z-score quantifies the difference between the groups in units of standard deviation, where a larger absolute z-score and a corresponding p-value < 0.05 indicate a statistically significant difference.30 In addition, we used the non-parametric Mann–Whitney U Test, appropriate for comparing two independent groups, in this case the ensemble group for each year and region, as shown in Table S1.
| Product | Sources | Timespan | Mean 2019–2021 (Tgy−1) | References |
|---|---|---|---|---|
| Carbon tracker – CH4 priors | Wetlands | 1984–2021 | 16.3 (±0.4) | 33 |
| Soil sink | −6.3 (±< 0.1) | |||
| Termites | 0.4 (±< 0.1) | 34 | ||
| Wild animals | 0.2 (±< 0.1) | 32 | ||
| Geological seep | 6.6 (±< 0.1) | 35 | ||
| Kaplan | Wetlands | 2006–2021 | 6.6 (±< 0.1) | 26 |
| Zhang | Wetlands CRU | 2001–2021 | 9.9 (±0.1) | 18 |
| Wetlands MERRA | 2001–2021 | 10.8 (±0.2) |
![]() | ||
| Fig. 1 Differences in 2020–2019 and 2021–2019 CH4 emissions over the U.S. from the inversion for (a) net emissions, (b) the natural component, (c) the anthropogenic component, (d) from the U.S. EPA inventory,19 (e) from the FOG oil and gas inventory,21 and (f) global fire emissions database.43 Error bars in panel (a) represent 1-sigma of the anomalies. Panel (g) shows the total number of cattle (heads) in the U.S.41 EPA and FOG columns show the net results in dashed lines. The points represent the average anomaly for the priors. | ||
To investigate the processes that may have caused the differences in CH4 emissions in 2020–2021 relative to 2019, we analyzed process-based information from anthropogenic emission inventories. The U.S. EPA inventory-based emissions (Fig. 1d) are also lower in 2020 and 2021 compared to 2019, consistent with our finding, although smaller in magnitude. The 2020 reduction in the EPA inventory was −0.4 Tgy−1, and, for 2021, it was −1.1 Tgy−1, or −1.5% and −3.8%, respectively, relative to 2019. The EPA data suggest that this reduction was mainly caused by the energy sector, whose anomalies were −0.8 Tgy−1 and −1.1 Tgy−1, or −6.7% and −9.7% for 2020 and 2021. The EPA-reported reductions of CH4 emissions in the energy sector in 2020–2021 are also consistent with the change reported by the FOG inventory, but with lower magnitudes. The net decrease of CH4 emissions in the oil and gas sector, estimated by the FOG inventory, is −1.1 Tgy−1 in 2020 and −1.6 Tgy−1 in 2021 relative to 2019 (Fig. 1e), representing −6.1% and −9.0% reductions. The reductions in FOG were dominated by reduced emissions from drilling, while emissions from other production activities were increased in 2020–2021. The CH4 emissions reduction associated with drilling was −2.8 Tgy−1 in 2020 and −1.9 Tgy−1 in 2021, similar to our optimized posterior anomalies. Besides reduction of CH4 emissions in the energy sector, there were also likely reductions of CH4 emissions from agriculture (Fig. 1g). Comparing U.S. Department of Agriculture (USDA) cattle inventory anomalies to EPA agriculture emission anomalies for 2020 and 2021, both relative to a 2019 baseline, reveals a consistent decreasing trend. USDA data show cattle inventory anomalies of −1.0% in 2020 and -1.2% in 2021 compared to 2019. The EPA's 2023 report indicates U.S. agriculture emission anomalies of 0.02 Tgy−1 (0.23%) in 2020 and −0.07 Tgy−1 (−0.76%) in 2021 relative to 2019. While the EPA's reported emission changes are smaller, particularly the slight increase in 2020, the decrease in both cattle inventory and overall agriculture emissions in 2021 relative to 2019 is consistent with the expected relationship between livestock numbers and a component of agricultural emissions, acknowledging that total agricultural emissions encompass various sources beyond just cattle. Wildfire CH4 emissions shown in Fig. 1f) increased during this period, similarly to Land Use and Land Use Forestry (LULUF) in Fig. 1d), discussed in more detail in next section.
The inversion system and spatial and temporal density of the measurements allow emissions to be derived on regional scales, which provide a spatial understanding of the changes during the COVID period. Hence, Fig. 2 shows the six regions considered in this study and Fig. 3 shows annual changes in total posterior emission changes by region. Consistent negative CH4 emission anomalies observed across multiple datasets in the oil and gas-dominated Central South (CS) region point strongly to reduced drilling activity as the primary driver during the COVID period. We found negative anomalies in the optimized posterior fluxes, and given that the natural emissions in this region are insignificant, the anthropogenic component is very similar to the posterior values.
![]() | ||
| Fig. 2 Regions in the U.S.: northeast (NE), southeast (SE), central north (CN), central south (CS), mountain (M), west (W), as defined by (ref. 44). The figure was generated using R package ggmap45 and ggplot2.46 | ||
![]() | ||
| Fig. 3 Difference 2020–2019 and 2021–2019 in CH4 emissions over the US Tgy−1, by regions in the US north east (NE), south east (SE), central north (CN), central south (CS), mountain (M), west (W), as defined by ref. 44 (a) year, (b) natural component, (c) anthropogenic component, (d) from U.S. EPA inventory,19 (e) by the FOG oil and gas inventory,21 and (f) Global Fire Emissions Database,43 and (g) by the total number of cattle (heads) in the US.41 Error bars in panel (a) represent 1-sigma of the anomalies. The points represent the average anomaly for the priors. | ||
The negative anomalies in the U.S. EPA19 are similar to our estimates and primarily arise from changes in the energy sector. The net anomaly in the FOG inventory also agrees with our results and those of the US. EPA.19 Specifically, in the region CS, for the periods 2020–2019 and 2021–2019, we found reductions of −0.7 (±0.31) Tgy−1 and −1.6 (±0.51) Tgy−1, in our posterior results. These changes are statistically significant after applying z-score test between the estimates of ensemble members, and Mann–Whitney U test, as shown in Fig. S5 and Table S1. In case of US EPA,19 we found −0.6 Tgy−1 and −0.9 Tgy−1 for 2020–2019 and 2021–2019; and for FOG21 −1.1 Tgy−1. There are some reductions in cattle in the CS during these years. In this region, we also see a significant reduction in drilling and increment of methane emissions associated with oil and gas production. Indeed, a study found that during the pandemic, lower demand for oil and gas resulted in turned-off wells, reducing CH4 emissions.42 Other studies also found substantial CH4 reductions in oil-dominated regions.23,47 There are similarities between negative FOG emissions in the drilling sector and negative optimized estimates. Fig. S6 shows correlations with interval confidence between posterior emissions and FOG emissions by production and drilling. Here, we see that the region CS has significant correlations between posterior and drilling CH4 emissions in 2020, corroborating that the main driver in reduced CH4 emissions in 2020 was likely due to reduced drilling activities. The percentage change in the FOG-drilling inventory in 2020 and 2021, about 2019, was −46.6% and −31.4%, respectively, while the same change for the FOG-production sector was 15.1% and 2.7%. This is because the COVID lockdown resulted in fewer drilling activities and associated emissions, but the cause of production emissions increase is unclear. The total change was −6.1% from 2019 to 2020 and −9.0% from 2019 to 2021.
The Central North (CN) region exhibited substantial negative anomalies in the U.S. CH4 emissions during the COVID period, representing one of the largest regional changes (Fig. 3). Specifically, for the periods 2020–2019 and 2021–2019, we found reductions of −3.4 (±0.62) Tgy−1 and −3.3 (±0.69) Tgy−1, in our posterior results (see Fig. S5 and Table S1). The FOG shows negative anomalies for both periods, and production showed a slight increase in 2020–2019 and a reduction in 2021–2019. Even though these anomalies agree with the emission decrease in our results, they do not explain the higher magnitude of the changes found in our study. The EPA shows slight negative anomalies, which also are not consistent with the values we found. According to the U.S. EPA Greenhouse Gas Inventory report,19 70% of the CH4 agricultural emissions come from enteric fermentation. Furthermore, among the agricultural categories (field burning of agricultural residues, enteric fermentation, manure management, and rice cultivation) only enteric fermentation exhibited strong and negative anomalies (Fig. S7), with a stronger anomaly in 2021 than 2020, similar to our results. Also, from the enteric fermentation sub-sector in the inventory,19 beef cattle are responsible for 71% of the totals for 2021 and show a more significant negative anomaly than other animals, as shown in Fig. S8 SI, which is consistent with the anomalies presented by the U.S. EPA, but our findings suggest a potential under-representation. The methodology for the agriculture sector is Tier 2 of IPCC.8 For the enteric fermentation, they use the Cattle Enteric Fermentation Model (CEFM), as shown in the Common Reporting Format (CRF) Source Category 3A.48 This methodology uses estimates of cattle populations by animal type and state, as shown in table A-126, Annex 3, greenhouse gases report in the US.48 Research has shown that 3-nitrooxypropanol (3-NOP), included in ruminants' diets, reduces methane emissions by about 30%.49,50 However, 3-NOP (BovaerTM) was approved by the Food and Drug Administration just in 2024,51 well after the period discussed here. Therefore, changes in the diet is not the cause to the reduction we observed at CN. During COVID-19 and the associated lockdowns, the livestock sector was particularly impacted because of a shortage of animal feed and health resources, workers, and a decrease in purchasing power.52 Also, sharp decreases in meat, milk, and egg production have been documented during this period, which is explained by mobility restrictions of workers, including veterinary services, among others.53 A map of the anomalies in livestock can be seen in Fig. S9. Furthermore, posterior emissions derived by multiple independent atmospheric inversions, i.e., from ref. 26 CarbonTracker-CH4, Copernicus Atmosphere Monitoring Service, and54 (Fig. S9 in ref. 26) in the CN region are higher than EPA inventory estimates, suggesting likely underestimation of agricultural emissions, and therefore reductions in agricultural emissions, by the EPA inventory.
In the other regions, the anomalies had smaller magnitudes. In the North East (NE) region, there is only a significant increment of emissions for 2020–2019 (see also Fig. S5 and Table S1) and not 2021–2019, which agrees with production in FOG and US EPA showing small negative anomalies. The posterior anomalies in the South East (SE) region were only significant for 2021–2019 and not 2020–2019. Coincidentally, the cattle also showed a small increase of around 20 thousand heads. Since the SE region is a primarily sub-tropical region with high precipitation and humidity levels, natural emissions may play a more pivotal role here. The natural anomaly for 2020–2019 is significant for this region. Wetlands are an important source of natural CH4 emissions,9,16 and over the US, the southeast has a high density of this source. Our model captures the natural variability since the optimized posterior values agree with the natural anomalies derived independently, as shown in Fig. 2, panel (b) for the difference 2021–2019 in the region NE.
In the Mountain (M) region, our posterior total emission anomalies were +0.3 (±0.33) Tgy−1 for 2020 relative to 2019 and +0.1 (±0.28) Tgy−1 for 2021 relative to 2019. Nevertheless, the anomaly was not significant after checking scores and Mann–Whitney U test (p.value > 0.05) in Fig. S5 and Table S1, respectively. For comparison, US EPA inventory anomalies for the same periods were +0.11 Tgy−1 and +0.02 Tgy−1, primarily attributed to the LULUF sector. The GFED fire emissions dataset also indicated positive anomalies of +0.14 Tgy−1 and +0.09 Tgy−1 for 2020–2019 and 2021–2019, respectively. A slight negative anomaly in cattle was also observed in Fig. 3g.
In the West (W) region, our analysis reveals positive posterior anomalies of +0.4 (±0.29) Tgy−1 (2020–2019) and +0.3 (±0.37) Tgy−1 (2021–2019). These results align well with19 inventory anomalies derived from the LULUF sector, reported as +0.4 Tgy−1 for both periods. Furthermore, the change was significant only for the period 2020–2019, also corroborated by the z-scores and Mann–Whitney U test Fig. S5 and Table S1. While the baseline LULUF sector emissions (per EPA definitions) are dominated by sources like flooded land remaining flooded (FLRF, 58.4%), which exhibit relatively stable emissions, our results indicate the observed positive anomaly in this region is driven primarily by Forest Fires (Fig. S10). This aligns with Forest Fires constituting 31.4% of baseline LULUF. California displayed the most significant state-level positive anomaly (Fig. S11), consistent with the record number of large fires in 2020 propagated by severe drought.55 Therefore, our optimized fluxes effectively capture the influence of wildfire activity on CH4 emissions, corroborating inventory trends.19 Further, GFED biomass burning anomalies for the W region were substantial at 0.49 Tgy−1 (2020–2019) and 0.37 Tgy−1 (2021–2019). Indeed, the annual time series from GFED confirms that 2020 saw record-breaking wildfire emissions, particularly impacting the M and W regions in Fig. S12 and S13.
To investigate the influence of natural emissions further, we calculated the anthropogenic anomalies at 1-degree resolution by subtracting the natural emissions from our posterior emissions. As shown in Fig. S14, we then compared these anthropogenic anomalies with the posterior anomalies. The strong correlation observed suggests that the posterior emissions are primarily driven by the anthropogenic component, allowing us to focus on the posterior results for further analysis. No statistical differences and high correlations were inferred after applying Mann–Whitney U and correlation tests between total optimized and anthropogenic components, as shown in panels (c) and (f).
Even though19 considered CH4 losses related to its reactions with OH, the OH field considered were monthly climatological values that do not contain interannual variability. Studies found that atmospheric OH concentration was reduced globally and over regions like the US during the 2020 pandemic period, largely due to decreased NOX emissions.17 Globally, tropospheric OH concentrations decreased by 1.6 ± 0.2% relative to 2019, a reduction that could explain 53 ± 10% of the coincident increase in global atmospheric methane abundance in 2020.17 Regional methane emission estimates are sensitive to interannually varying OH fields.56 Over the U.S., the tropospheric OH reduction was approximately 5% at maximum due to pandemic-related NOX emission decrease.57 This reduction in OH has implications for methane budget calculations. For example, if less OH is present in the atmosphere, the atmospheric lifetime of CH4 becomes longer. This increased lifetime enhances the contribution of background CH4 arriving from remote regions to the concentrations observed at a measurement site and smaller concentration enhancements would be attributable to emissions. For instance, assuming a linear sensitivity to the pandemic-related 5% reduction in tropospheric OH to our 2020 posterior estimate of 56.7 Tgy−1 could increase the reported reduction of −2.5 Tgy−1 to −2.8 (2.6 to 3.0) Tgy−1. This further supports our results as a conservative estimate of the pandemic's impact. Therefore, if reduction in OH was considered in ref. 26, a slightly larger reduction of CH4 emission would be expected to be derived.
While our regional inversion26 used a monthly climatological OH field, we now explicitly acknowledge that a global drop in OH during 2020–2021, driven by reduced NOX and CO emissions, likely explained a dominant portion of the global methane surge.11,12 Our finding of reduced U.S. fossil fuel emissions (particularly in the Central South) aligns with latest isotopic evidence showing that the record global surge in 2020–2022 was driven almost entirely by microbial sources, while fossil fuel emissions remained stable or declined in many regions.14,15
Supplementary information (SI): assimilated methane monitoring network (aircraft and surface sites); regional optimization summaries for total, natural, and anthropogenic emissions alongside 2020 and 2021 anomalies; statistical evaluations (z-scores and Mann–Whitney U tests) of the regional inversions; and sector-specific data covering fossil oil and gas activity, agriculture, and biomass burning. See DOI: https://doi.org/10.1039/d6ea00004e.
Footnote |
| † Now at Cooperative Institute for Satellite Earth System Studies (CISESS)/Earth System Science Interdisciplinary Center (ESSIC), University of Maryland. |
| This journal is © The Royal Society of Chemistry 2026 |