January 2023: | 1921.74 ppb |
January 2022: | 1907.44 ppb |
The graphs show globally-averaged, monthly mean atmospheric methane abundance determined from marine surface sites. The first graph shows monthly means for the last four years plus the current year, and the second graph shows the full NOAA time-series starting in 1983. Values for the last year are preliminary, pending recalibrations of standard gases and other quality control steps. Other impacts on the latest few months of data are described below.
The Global Monitoring Division of NOAA’s Earth System Research Laboratory has measured methane since 1983 at a globally distributed network of air sampling sites (Dlugokencky et al., 1994). A global average is constructed by first smoothing the data for each site as a function of time, and then smoothed values for each site are plotted as a function of latitude for 48 equal time steps per year. Global means are calculated from the latitude plot at each time step (Masarie and Tans, 1995). Go here for more details on how global means are calculated.
Methane is reported as a “dry air mole fraction”, defined as the number of molecules of methane divided by the total number of molecules in the sample, after water vapor has been removed. The mole fraction is expressed as nmol mol^{-1}, abbreviated “ppb” (for parts per billion; 1 ppb indicates that one out of every billion molecules in an air sample is CH_{4}).
In the figures, the red lines and circles are globally averaged monthly mean values centered on the middle of each month. The black line and squares show the long-term trend (in principle, similar to a 12-month running mean) where the average seasonal cycle has been removed.
Year | Growth (ppb) | Uncertainty (ppb) |
---|---|---|
1984 | 12.82 | 0.88 |
1985 | 12.30 | 0.98 |
1986 | 12.90 | 0.72 |
1987 | 11.40 | 0.75 |
1988 | 10.79 | 0.47 |
1989 | 11.08 | 0.59 |
1990 | 8.75 | 0.51 |
1991 | 14.06 | 0.59 |
1992 | 2.36 | 0.50 |
1993 | 3.83 | 0.52 |
1994 | 7.22 | 0.49 |
1995 | 3.96 | 0.50 |
1996 | 2.50 | 0.37 |
1997 | 6.27 | 0.74 |
1998 | 12.14 | 0.74 |
1999 | 2.40 | 0.63 |
2000 | -1.48 | 0.62 |
2001 | -0.72 | 0.48 |
2002 | 3.24 | 0.55 |
2003 | 4.79 | 0.56 |
2004 | -4.86 | 0.36 |
2005 | 0.37 | 0.48 |
2006 | 1.91 | 0.57 |
2007 | 7.81 | 0.61 |
2008 | 6.56 | 0.43 |
2009 | 4.66 | 0.49 |
2010 | 5.20 | 0.67 |
2011 | 4.87 | 0.64 |
2012 | 5.01 | 0.56 |
2013 | 5.70 | 0.60 |
2014 | 12.77 | 0.50 |
2015 | 10.02 | 0.66 |
2016 | 7.10 | 0.60 |
2017 | 6.84 | 0.72 |
2018 | 8.70 | 0.52 |
2019 | 9.70 | 0.60 |
2020 | 15.21 | 0.41 |
2021 | 17.70 | 0.47 |
2022 | 14.07 | 0.58 |
The table and plot summarize annual increases in atmospheric CH_{4} based on globally averaged marine surface data.
The annual increase in atmospheric CH_{4} in a given year is the increase in its abundance (mole fraction) from January 1 in that year to January 1 of the next year, after the seasonal cycle has been removed (as shown by the black lines in the figure above). It represents the sum of all CH_{4} added to, and removed from, the atmosphere during the year by human activities and natural processes. Our first preliminary estimate for the annual increase of a particular year is produced during April of the following year, using available data from the previous year. It is important to recognize that the initial, April estimate of the annual increase is likely to change significantly as more data are added to the analysis. That estimate will be updated in subsequent months as more samples are measured for CH_{4} and included in the analysis. By autumn of the following year the annual increase will typically converge toward a “final” value.
Estimates of the globally-averaged CH_{4} abundance (monthly- and annually-averaged means), and the annual increase, are updated every month as new samples are returned to Boulder, measured for CH_{4}, and added to the analysis. Adding new, more recent data improves the accuracy of the initial estimate by increasing the spatial density of data and eliminating “end effects” of the curve fitting procedures used. We’ve investigated the impacts of adding new data to the parameters reported here, and a summary of the results follows:
Initial estimates of the CH_{4} annual increase made in April for the previous year are biased compared to those that follow using additional data. The average bias in the initial estimate is +1±0.8 ppb yr^{-1} (1 standard deviation shown). Over the next few months, the average bias slowly decreases until it is negligible by July or August. In any given year though, bias in the initial estimate of the annual increase can be much larger than the average, with bias up to ±3 ppb yr^{-1}; that is, it can be positive or negative. In other words, until late in a year, bias in the annual increase can be much larger than the uncertainty reported based on the bootstrap method described below.
Behavior of initial annually-averaged means and monthly-averaged means are similar (see links to files below). For monthly mean CH_{4}, the initial value is typically too high, by up to 7.6 ppb.
The estimated uncertainty in the global annual CH_{4} increase varies by year. It is estimated using two terms: The first is a "bootstrap" (resampling) method that varies the sites in our network. Each bootstrap realization of the network is constructed by randomly picking sites, with restitution, from the existing marine boundary layer sites in the NOAA/GML cooperative global air sampling network (Dlugokencky et al., 1994). Each member of the ensemble of networks has the same number of sites as the real network, but some sites are missing, while others are represented more than once. An additional condition is that at least one site is present from high southern latitudes, one from the tropics, and one from high northern latitudes, because we have always maintained broad latitude coverage in the real network. Temporal data gaps at individual sites are present in the bootstrap networks. The second term is a Monte Carlo method that randomly modifies the data to account for measurement uncertainty. The modifications are based on assessment of random uncertainty in the measurements, and it varies over time. In both cases, 100 globally-averaged time series are created. We calculate the mean and standard deviation for each year's annual increase from the ensemble members, and one standard deviation from the two terms (network and analytical) is taken in quadrature to give the reported uncertainty at each time step. As mentioned earlier, bias in our first estimates of annual increase, monthly mean, and annual mean can be significantly greater than the stated uncertainty.