Methane (CH4) is a potent greenhouse gas (GHG) with significantly greater heat-trapping capacity than carbon dioxide (CO2). Landfill gas (LFG), generated by anaerobic degradation of organic waste in municipal solid waste (MSW), consists mainly of equimolar CH4 and CO2. Landfills, which encompass both uncontrolled dumpsites and sanitary landfills, represent a substantial source of anthropogenic CH4 emissions (Cai et al., 2018; Zhu et al., 2023). The global annual CH4 emissions from landfills are estimated to be on the order of several tens of Tg (Fig. 1a), accounting for ~ 20% of all annual anthropogenic CH4 emissions (Saunois et al., 2020). In certain megacities, the proportions of CH4 emissions from landfills even exceed 50% (Maasakkers et al., 2022).
There are approximately 300,000-500,000 existing landfills worldwide, receiving > 1.5 billion tonnes of MSW annually (Kaza et al., 2018). Some huge landfills emit CH4 at rates comparable to those of “ultra-emitters” in oil and gas industry (Lauvaux et al., 2022), not to mention the vast number of small landfills with continuous CH4 leakage. Nonetheless, landfill CH4 emissions receive disproportionally low attention in research and mitigation compared to other sources with similar emission levels. For example, the energy sector emits a comparable amount of CH4 annually, but requires 10 ~ 100 USD/t CO2-eq more mitigation cost than the waste sector (EPA, 2019). Furthermore, landfill emissions have an estimated uncertainty range of 10 ~ 20 Tg/year, the value of which already exceeds the total emission of wastewater treatment that received considerable attention recently (Du et al., 2023; Song et al., 2023) (Fig. 1a). A more accurate assessment of landfill CH4 emissions is essential to formulate effective waste management strategies and emission mitigation strategies (Fei et al., 2021) to realize the targeted 30% reduction in global annual CH4 emissions by 2030, as pledged in the 2021 United Nations Climate Change Conference (COP26).
The current estimation of global landfill CH4 emissions relies on the reporting of national greenhouse gas inventories (NGHGI) to the United Nations Framework Convention on Climate Change (UNFCCC), which follows the guideline set by the Intergovernmental Panel on Climate Change (IPCC). Landfill managers and government officers most often use a first-order decay (FOD) model (Eq. 1) to calculate site- or region-specific annual landfill CH4 emissions based on waste disposal amount and landfill operating time (IPCC, 2006).
$$\begin{array}{c}{r}_{CH4}={L}_{0}\times W\times k\times {e}^{-k\times (t-{t}_{0})}\#\left(1\right)\end{array}$$
where rCH4 is the gravimetric CH4 generation rate (kg CH4/year), W is the waste disposal amount (dry or wet mass, kg), t is the time since waste disposal (year), and t0 is the lag time between waste disposal and the initiation of CH4 generation (year).
The FOD model (Eq. 1) necessitates another two critical input parameters, namely first-order decay constant (also known as waste decay rate, k, year− 1) and methane generation potential (L0, L CH4/kg dry waste) (IPCC, 2006). The model was initially adopted by IPCC in 2006 and remained largely unchanged thereafter, except for a minor refinement in 2019 (IPCC, 2019). The values of k and L0 are subject to significant variability due to waste heterogeneity and landfill conditions. The value of L0, an inherent property of waste components, has been comprehensively reviewed by previous research (Krause et al., 2016), and can be accurately determined based on waste composition (IPCC, 2006; Krause et al., 2016). However, the IPCC’s guideline only provides four default k values based on different temperature and precipitation conditions. This fails to consider the diversity of landfills with different geographical locations, waste compositions, and potential climate change impacts, resulting in substantially biased k estimations. However, 198 countries and regions under UNFCCC framework report their landfill CH4 emission inventories employing IPCC’s default k values, thereby underlining the significant magnitude of potential CH4 emission misestimations and importance of k value refinement.
Besides IPCC’s inventory method (bottom-up) using a generalized model to calculate landfill CH4 emissions, remote sensing of CH4 concentration and atmospheric inversion modeling (top-down) provides direct and independent quantifications of CH4 fluxes from landfills (Deng et al., 2022; Erland et al., 2022). The biases of the conventional “bottom-up” inventory estimations have been acknowledged by previous studies (Spokas et al., 2015; Wang et al., 2015) and further corroborated by the recent “top-down” inversions leveraging state-of-the-art satellites (Maasakkers et al., 2022) and aircrafts (Duren et al., 2019), which reveal up to 2.6 times more CH4 emissions from certain landfills compared to the inventory estimations.
Nevertheless, the extensive data processing requirement, high data acquisition cost, and high personnel expertise threshold of atmospheric inversion limit its applicability to a large number of landfills. In contrast, as the IPCC’s inventory method is widely accessible and user-friendly, we seek to improve it by revising the k estimation method for the FOD model. Here, contrasting the previous studies relying on the k values measured from a single landfill (De la Cruz et al., 2016; Delgado et al., 2023) or a small number of laboratory tests (Karanjekar et al., 2015), we survey and analyze exhaustively reported k (kr) values obtained from landfill measurements (N = 196 in 18 countries) and laboratory tests (N = 80 in 12 countries) in the literature (Fig. 1b). The operating and environmental conditions of each landfill are also collected to establish a comprehensive landfill database.
Thereafter, a new k estimation method based on the three most influential factors, site-specific waste composition, moisture, and temperature corrections (CMT method) is developed. The superior performance of the CMT method is showcased by comparing the corrected k (kCMT) with the IPCC’s default k (kIPCC) following a two-stage approach. In the first stage, we compare the k estimation error between CMT method and IPCC’s method by referencing the reported kr values. In the second stage, we contrast the CH4 prediction results based on both kCMT and kIPCC, utilizing the atmospheric inversion results as independent benchmarks. Subsequently, we apply the kCMT to predict the cumulative CH4 emissions between 2010 and 2030 from 12 major landfills in different climate zones and income-level countries worldwide. The biases between the predictions made using kCMT and kIPCC are gauged. Lastly, we explore the novel insights on landfill CH4 emission trends uncovered by the kCMT predictions and potential mitigation strategies (Fig. 1c and 1d).
Survey of reported waste decay rates worldwide
This study conducts a comprehensive survey and analysis of measured and reported waste decay rates (kr) on a global scale. We exhaustively compile available kr values and distinguish them by the sources into field measurements and laboratory tests (Fig. 2a and 2b). Compared to kIPCC, the kr values are considered to be more accurate for individual landfills, which distribute broadly across all major climate zones (except the polar zone) with an annual average temperature of 3.1 ~ 28.1 ℃ and annual average precipitation of 79 ~ 2,344 mm. Most field measurements and laboratory tests have been conducted in the United States (U.S.), Europe, and China. Brazil represents the sole South American country with trackable CH4 emission data from landfills, while no reliable data exists for African countries. In developing countries, where field measurements or laboratory experiments are not practicable, kIPCC are extensively utilized, with over 100 studies employing kIPCC to estimate site- or city-level CH4 emissions from landfills.
Field and laboratory kr values complement each other to improve CH4 emission estimation and understanding of influential factors. Laboratory tests typically create controlled and ideal environments for MSW degradation, yielding high kr and nearly complete time series data of CH4 generation (He and Fei, 2020). Conversely, field kr measurements are the most valuable first-hand data to quantify landfill CH4 emissions. Nonetheless, field measurements encompass considerably higher uncertainties than laboratory tests and time series data of CH4 generation data is often incomplete. 92% of the field kr are between 0.01 ~ 0.3, with a few outliers resulting from unique operating conditions, e.g., kr = 0.87 for a pilot-scale bioreactor test cell (Yazdani et al., 2012). The laboratory kr values are more dispersed between 0.8 ~ 25.6. This is attributed to the widely varied operating conditions of laboratory tests, e.g., kr = 25.6 for a small test column with intensive leachate recirculation (Pezzolla et al., 2017) and kr = 0.8 for a larger column with a lower recirculation frequency (Huang et al., 2012). The average kr of the laboratory tests (5.0) substantially surpasses that of field measurements (0.12) due mainly to scale effect and degradation-promoting conditions (Fig. 2b) (Fei et al., 2016). Factors contributing to the scale effect between field and laboratory kr include operating conditions (moisture, temperature, inoculation, vertical stress, oxygen availability, measurement accuracy, etc.) and waste properties (size of interest, heterogeneity, etc.). Both the scale effect and degradation-promoting conditions on k estimation are considered in this study.
The k value is influenced by factors such as moisture, temperature, waste and leachate properties, and landfill operations, among which the IPCC’s guideline identifies moisture, temperature, and waste composition as the three primary factors. Our survey reveals that the kr values exhibit substantial variation between annual temperature of < 20 and > 20 ℃ (p = 0.03, Fig. 2c), but statistically insignificant variation between annual precipitation of < 1,000 and > 1,000 mm (p = 0.13, Fig. 2d). This does not imply that precipitation has a negligible impact on kr values; rather, it suggests that the dichotomous classification adopted by the IPCC’s guideline (precipitation: wet or dry, temperature: tropical or temperate) cannot reflect quantitatively the impacts of the influencing factors. Consequently, it is crucial to conduct a quantitative analysis to distinctly characterize the impacts of temperature and moisture on k.
Moreover, significant disparities in kr values are observed among the countries (Fig. 2e), primarily due to pronounced differences in waste composition among them. The difference arises from distinct income levels, lifestyles, product availabilities, and waste management strategies. The Annex I countries of UNFCCC, predominantly Organisation for Economic Co-operation and Development (OECD) countries, typically generate lower proportions of fast degradable waste (Bf in %, including food and yard waste) than the Non-Annex I countries (IPCC, 2019). In many European countries that advocate waste incineration, such as Germany and Denmark, Bf even approaches 0% in landfiled waste (Kaza et al., 2018). Conversely, Bf constitutes substantial portions of disposed waste in many developing countries, e.g., Bf = 59% in China and 53% in India (IPCC, 2019). A gap remains in quantifying the influence of waste composition on estimated k, which we aim to address in this study. Additionally, LFG recovery efficiency and type and placement time of landfill cover vary across countries, which also impact kr and estimated k. Overall, the significant discrepancy between the kr and kIPCC highlights the urgent need for improving k estimation to be used in the IPCC’s method.
Improved kCMT estimations are frequently closer to kr than kIPCC
To improve the estimation accuracy of k and CH4 emission and reduce the potential bias induced by kIPCC, we develop a CMT correction method based on the analysis of available kr values worldwide. The CMT method quantifies the influences of moisture, temperature, and waste composition on k (Eq. 2), rather than categorizing k value dichotomously as has been done in IPCC’s method. The details of the development and application procedures of the CMT method are provided in the Methods.
$$\begin{array}{c}{k}_{CMT}={k}_{0}{\times f}_{C}\times {f}_{T}{\times f}_{M}\#\left(2\right)\end{array}$$
where kCMT refers to k after CMT corrections; k0 refers to the k under baseline moisture, temperature, and waste composition; fC, fM, and fT refer to the correction factors for waste composition (C), moisture (M), and temperature (T), respectively. k0 represents the maximum k value under optimum conditions where fC = fM = fT = 1.
We first determine the kCMT and kIPCC for the 196 landfills with reported kr as described in the previous section. The relative errors of kIPCC (δ(kIPCC)) and kCMT (δ(kCMT)) as compared to the corresponding kr are calculated. The kr values are considered to be the most accurate available data for the landfills, while the relative error serves as an indicator of the disparity between the estimated k and the actual k value. A higher relative error denotes a diminished accuracy of the k estimation method. The average δ(kCMT) is 56.2%, which is significantly lower than the average δ(kIPCC) of 78.7% (p = 0.05). The δ(kCMT) values are < 50% (taken as the acceptable error) in 134 out of the 196 landfills (Fig. 4a), while 104 landfills show δ(kIPCC) < 50% (Fig. 4b). The error analysis reveals that the kCMT values are superior in 105 out of the 196 landfills and equivalent in 45 out of the 196 landfills as compared to the kIPCC values, giving kCMT a 77% likelihood of being equally or more accurate than kIPCC (Fig. 4c). While retaining the essence of kIPCC, that is simplicity and a priori estimation, kCMT improves significantly the accuracy of k estimation.
We recognize that high deviations between kCMT and kr still exist, e.g., kCMT = 0.09 and kr = 0.01 for a landfill in Maryland, U.S. (Jain et al., 2021). One explanation for the high δ(kCMT) could be the uncertainties in determining fT, fM, and fC for kCMT calculation. In addition, due to limited data availability, some landfills reported the waste composition during collection instead of disposal, which becomes a contributor to δ(kCMT). The accuracy of kCMT can be improved if more landfills with complete site information and time-dependent CH4 generation data are included in the analysis, whereas the accuracy of kIPCC is unaffected per its standard calculation procedure. The generally lower δ(kCMT) than δ(kIPCC) highlights that more accurate CH4 emission estimations can be made by using kCMT than kIPCC, as demonstrated in the next section.
Inversed atmospheric CH 4 fluxes are better estimated by QCMT than QIPCC
Recent advancements in high-resolution remote sensing technology using satellites and unmanned aerial vehicles have enabled atmospheric inversion as a powerful tool for measuring CH4 emissions from single landfills (Maasakkers et al., 2022; Tu et al., 2022). We gather 32 landfills (Fig. 2a) with available inversed atmospheric CH4 fluxes (Qr, Gg/y) and sufficient operating information as benchmarks against the corresponding CH4 emissions estimated using kCMT and kIPCC. For each landfill, we estimate the annual CH4 emissions using both kCMT (QCMT, Gg/y) and kIPCC (QIPCC, Gg/y), and compared these results with the corresponding Qr to yield the respective errors δ(QCMT) and δ(QIPCC). These landfills are not purposely chosen but represent the limited number of sites meeting the selection criteria (Fig. S2c). The comparisons underscore the significant amount of CH4 underestimated by the current IPCC’s method. Apart from one landfill with a Qr lower than the QIPCC (Victorville, U.S., Fig. S1), all the other Qr exceed QIPCC by 4-737%. The average δ(QCMT) and root mean square error (RMSE) of QCMT are 36% and 7.3 for the 32 landfills, while the δ(QIPCC) and RMSE of QIPCC are 45% and 14.5. The QCMT values are closer to the Qr than QIPCC in 26 out of the 32 landfills (Fig. 4i), further demonstrating that adopting kCMT generates more accurate estimations than using kIPCC. It is noteworthy that k only controls the annual emission rate, but not the cumulative emission amount, which is instead governed by L0 (Krause et al., 2016).
Although the QCMT is more accurate than the QIPCC for most of the landfills, the improvements by QCMT as compared to the QIPCC ((QCMT-QIPCC)/QIPCC×100%) vary between 9 ~ 177%. The difference is mainly attributed to landfill age and climate conditions. The CH4 generation and emission rates of a new landfill increase exponentially according to the FOD model, thus the difference between kCMT and kIPCC leads to an increasingly widened gap between the QIPCC and QCMT (e.g., Fig. 4b). As the annual CH4 emission reaches a peak while assuming a relatively stable annual waste disposal rate, the QCMT and QIPCC values tend to intersect (e.g., Fig. 4a, 4c, and 4d). Considering the climate conditions, the landfills in tropical wet zones have comparable kIPCC and kCMT values. Thus, the QIPCC and QCMT are also comparable. For instance, Kanjurmarg landfill in Mumbai, India, has kIPCC = 0.17 year− 1 and a waste decay half-life (t1/2,IPCC) of 4.1 years. In contrast, the kCMT = 0.3 year− 1 and t1/2,CMT = 2.3 years for this landfill. With an increase of 0.13 year− 1 from kIPCC to kCMT, the t1/2 is shortened by only 1.8 years, and the QCMT and QIPCC in 2021 differ by just 25% (Fig. 4c). Conversely, the West Miramar landfill in San Diego, U.S., in temperate dry zone has kIPCC = 0.02 year− 1 and kCMT = 0.06 year− 1, which is a three-time difference. This leads to the t1/2,CMT to be 23.1 years lower than the t1/2,IPCC and 70% higher QCMT than QIPCC in 2017 (Fig. 4f). The other landfills (Fig. 4b, 4e, 4g, and 4h,) have a similar trend as in this case.
Overall, the differences between kCMT and kIPCC and QCMT and QIPCC are particularly evident in new landfills in arid, temperate, and continental climate zones. This is a major finding that has been overlooked so far, as it highlights the sensitivity of time-dependent Q values to the accuracy of k values. In fact, the biased QIPCC are not confined to these 32 landfills, but most landfills worldwide. We apply the k correction method to additional landfills worldwide and discuss new perspectives and mitigation strategies for landfill CH4 emissions in the next section.
ΣQ CMT predictions are equal to or much higher than ΣQIPCC by 2030
Table 1
Underestimation (Qdiff%)a of CH4 emissions by current inventories in 12 typical landfills situated in different climate zones and the annual MSW generation amount and share in each climate zone.
Climate zone | A: tropical | B: arid | C: temperate | D: continental | E: polar |
Annual MSW generation in 2018 (tons/year) | 546,099,396 | 677,757,999 | 629,788,663 | 218,461,487 | 275,270 |
Share of global annual MSW generation in 2018 | 26% | 33% | 30% | 11% | 0.0% |
Qdiff% for 12 major landfills | Brasilia, Brazil (1%); Mumbai, India (2%); Guangzhou, China (11%); Jakarta, Indonesia (12%); | Pietermartzburg, South Africa (40%); Arequipa, Peru (49%); Bennett, U. S. (209%); | Shanghai, China (22%); Randleman, U.S. (30%); London, U.K. (47%); | Belgrade, Serbia (25%); Morrisville, U.S. (182%); | N.A. |
a Qdiff% = (ΣQCMT - ΣQIPCC) / ΣQIPCC × 100% (by year 2030) |
b Estimated based on the annual MSW generation rate (Kaza et al., 2018) and primary climate zone of each country, except for differentiating different provinces in China and different states in the U.S. (Kottek et al., 2006).
The higher accuracies of kCMT and QCMT compared to kIPCC and QIPCC establish a foundation for extending the CMT correction method and QCMT to different landfills and countries globally. Given the unique site conditions of each landfill, it is overwhelming and unfeasible to reassess all landfills worldwide using the CMT method in a single study. Alternatively, we select 12 representative major landfills in various climate zones and countries to predict their cumulative CH4 emissions since the operation to 2030 using kIPCC (ΣQIPCC, Gg CH4) and kCMT (ΣQCMT, Gg CH4) (Fig. 5). These chosen landfills are among the largest sanitary landfills and uncontrolled dumpsites worldwide (Fei et al., 2022), situating across all major geographical regions (East Asia, South Asia, Southeast Asia, East Europe, West Europe, Africa, South America, and North America) and climate zones (A: tropical, B: arid, C: temperate, D: continental, and E: polar).
The ΣQCMT values of the 12 landfills are on average 52% higher than the ΣQIPCC values. The landfills in arid, temperate, and continental climate zones have higher differences between ΣQIPCC and ΣQCMT (Qdiff%) than those in the tropical zone. The four landfills in the tropical zone show only 1–12% of Qdiff%, the reason for which has been discussed in the previous section. The eight landfills in arid, temperate, and continental zones reveal considerable underestimations by ΣQIPCC as compared to ΣQCMT. In the arid zone, where 33% of global MSW is generated annually, the Qdiff% reaches 209% in a 6-year-old landfill in Bennete, U.S. In temperate zone, where 30% of global MSW is generated annually, the Qdiff% reaches 47% in a landfill in London, the United Kingdom (U.K.). In the continental zone with 11% of global annual MSW generation, the Qdiff% reaches 182% in a 6-year-old landfill in Morrisville, U.S. (Table 1). If we extrapolate such level of underestimation to all landfills in arid, temperate, and continental climate zones, the potentially overlooked annual CH4 emissions in 2030 amounts to 10 ~ 20 Tg/yr, which is already commensurate with the total annual CH4 emission from the wastewater sector.
The overlooked CH4 emissions from landfills are informative to the operations and monitoring throughout their lifespans. A considerable portion of underestimated ΣQ occurs during early landfilling stage, when LFG collection is not implemented. To achieve the best economic and environmental benefits, landfill managers should determine flexibly the installation and operation time of LFG collection systems based on the kCMT values instead of the kIPCC values or a fixed delay time. Moreover, the faster waste degradation rates as represented by the higher kCMT values than the kIPCC values indicate shorter waste stabilization time after closure, which suggest that the recommended post-closure care periods of at least 30 years may be shortened on site-specific cases (EPA, 2016). Once the post-closure care period is over, the landfill can be redeveloped to alleviate land scarcity. The higher QCMT than QIPCC of most landfills indicates that the unit costs of implementing CH4 emission mitigation measures in landfills may be even lower than the current expectations of 0–10 USD/t CO2-eq (IPCC, 2023). On the other hand, the benefit of shifting waste management from disposal to recycling and incineration is also underestimated, as the current cost-benefit analysis of waste management is based on the kIPCC and QIPCC of a landfill, which underestimate its CH4 emission compared to reality. Overall, the newly introduced kCMT uncovers new mitigation opportunities obscured by the current kIPCC.