1. Introduction
Landfills are the most straightforward, cost-effective, and commonly used waste disposal systems worldwide. Despite recent advancements in reuse, recycling, and recovery methods, these efforts have not been able to halt the global expansion of landfills [
1]. According to UN projections [
2], a significant increase is anticipated in the total mass of waste generated, as well as in the amount of waste deposited in Municipal Solid Waste (MSW) landfills, as illustrated in
Figure 1.
From 2.1 billion tonnes in 2023, the generation of MSW is expected to grow to 3.8 billion tonnes by 2050. Of this, the mass of waste in MSW landfills is projected to increase from 0.64 to 1.09 billion tonnes, representing an approximate growth of 70%. While the direct cost of global waste management in 2020 was estimated at USD 252 billion, the inclusion of indirect costs associated with pollution, health impacts, and climate change from poor disposal practices raises this figure to USD 361 billion. The report warns that inaction on waste management could nearly double these costs by 2050, reaching a staggering USD 640.3 billion. These findings highlight the urgent need for improved waste management strategies to mitigate the economic and environmental burdens of MSW [
2].
MSW landfills are geotechnical structures; however, unlike traditional geotechnical structures, they are composed of waste rather than soil. According to existing research [
3,
4], the geotechnical properties of MSW relevant to slope stability are significantly more heterogeneous than those of soil. However, stability analyses often treat them as if they were soil. This approach can result in slopes that are insufficiently stable and robust, where even minor operational deficiencies during the lifespan of MSW landfills can lead to failure, as evidenced by a series of historical catastrophic MSW landfill slope failures. Such events pose significant risks to the environment, public health, and the economy.
An example is the Payatas landfill failure, which occurred on July 10, 2000, in the Philippines, with severe consequences for society and the environment [
5]. A total of 13,000–16,000 cubic meters of municipal solid waste (MSW), weighing 15,800–19,400 tons, was involved in the landslide. The causes of the failure were linked to the slope geometry (a slope that was too steep), large amounts of rainfall, and operational deficiencies. A series of MSW landfill failures have occurred worldwide, such as in Sarajevo, Yugoslavia, in 1977 [
6]; Bogota, Colombia, in 1997 [
7]; and Bandung, Indonesia, in 2005 [
8].
Ensuring adequate slope stability in MSW landfills is one of the most critical engineering challenges in landfill operations [
9]. Various authors have conducted back-analyses of MSW landfill failures. Zhang et al. [
10] analyzed 62 slope failures of MSW landfills, categorizing them into three groups based on the location of the slip surface. In the first group, the entire slip surface lies within the waste pile. In the second group, it passes through the foundation soil, and in the third group, it forms at the interface between the bottom liner and the MSW body. Blight [
11] analyzed six large-scale MSW landfill failures that occurred between 1997 and 2005. The study highlights that the causes of failure in the cases reviewed include the absence of engineered designs for landfills, lack of engineering oversight throughout the construction and operation phases, and the most important factor overlooked in all these cases: the influence of both leachate levels and slope angles on the overall stability of the landfills. Various authors have conducted back-analyses of the stability and reliability of MSW landfill slopes using different methods. The most common methods employed are the Limit Equilibrium Method (LEM), the Finite Element Method (FEM), and the Finite Difference Method (FDM).
The Limit Equilibrium Method (LEM) utilizes the concept of static equilibrium, calculating the factor of safety based on the principles of statics by considering force and/or moment equilibrium. This method is widely used in geotechnical engineering practice. However, its major shortcoming is the disregard for strain and displacement compatibility [
12]. In contrast, discrete element formulations offer significant advantages when modeling failure and post-failure behavior. They allow individual particles, blocks, or elements to separate, move, and rotate independently, potentially forming new contacts with other blocks. This approach is particularly valuable for analyzing failure mechanisms and large post-failure deformations. However, their drawback is that they do not provide a direct indication of the margin of safety [
13].
A study by Vinai et al. [
14] demonstrates the application of the Distinct Element Method (DEM) to describe the behavior of geogrid-reinforced embankments under rock impact. The researchers simulate various types of rock impacts on an embankment using Itasca’s PFC2D software, version 3.0. The modeling outputs include diagrams that illustrate displacement and damage, such as failures and tension cracks. These results are validated against full-scale impact test data. The study concludes that PFC2D modeling satisfactorily describes the real-world behavior of embankments under rock impact; however, such an approach does not directly quantify the margin of safety. Medizadeh et al. [
1] employ the Finite Difference Method (FDM) in conjunction with the Monte Carlo method to determine the reliability of MSW landfill slopes. In reliability analyses, they select soil unit weight, cohesion, and friction angle as random variables and assume a normal statistical distribution for the safety factor, determined through 2000 simulations. They conclude that this number of simulations is appropriate for assessing the reliability of MSW slopes and that higher coefficients of variation (COV) of random variables result in lower reliability. Chropenova et al. [
9] analyze the stability of two MSW landfills near the town of Handlová in Slovakia. They conduct analyses under both static and seismic conditions, with a horizontal acceleration factor of Kh = 0.07 (pseudo-static analyses). These analyses utilize the Bishop and Spencer methods (Limit Equilibrium Methods), assuming both circular and polygonal slip surfaces that pass through the waste body and foundation soil. The Spencer method and the polygonal slip surface result in lower safety factors, ranging from 1.70 to 2.01 under static conditions and from 1.41 to 1.59 under seismic conditions.
Research Purpose
The aim of this study is to quantify the safety margins and reliability of five MSW landfill slopes that have experienced failures by applying Eurocode standards. The following cases were analyzed: Cincinnati, OH, USA [
12]; Xiaping, Shenzhen, China [
10,
15]; Maoershan, Sichuan, China [
10]; Meethotamulla, Sri Lanka [
10,
16,
17]; and Xerolakka, Greece [
18,
19]. Relevant data for conducting stability and reliability analyses were sourced from the literature [
10,
15,
16,
17,
18,
19,
20], where back analyses determined the pre-failure geometry, failure mode, material properties, leachate levels, and other geotechnical parameters. From the literature review, it was determined that in these MSW landfill failures, authors primarily quantified slope stability using the traditional factor of safety, which was around a unit value in all analyzed cases. However, no data were found regarding the stability or reliability of these MSW landfills according to Eurocode 7 or any other building code. The landfills in this study were selected based on the availability of detailed pre-failure geometry and geotechnical data, enabling comprehensive analyses of stability and reliability. Additionally, these cases represent the most common causes of failure, including poor compaction, inadequate foundation soil, high leachate levels, steep slopes, and operational deficiencies. In all the cases considered, the slip surface passes through the waste body, which accounts for 70% of all landfill slope failures [Qian, Translational Failure Analysis of Landfills].
The research hypothesis posits that in all analyzed cases, the reliability of the MSW landfill slopes is significantly lower than typical values, indicating a much higher probability of failure. By conducting reliability analyses, the risk of failure could have been identified in advance, potentially allowing for the implementation of preventive measures to avert the failures.
In this study, the margins of stability were assessed using the factor of safety (FS) and the over-design factor (ODF), while reliability was quantified using the reliability index calculated by the Monte Carlo method. The study also analyses the influence of seismic loading on slope stability. Additionally, sensitivity analyses were conducted to investigate the effects of individual parameters on the stability of MSW slopes.
The findings of this study indicate that in all analyzed cases, under static loading, the ODF was significantly lower than the prescribed value, failing to meet the criteria set by Eurocode 7. Seismic loading further significantly reduced the ODF value. Reliability analyses indicated insufficient levels of the reliability index in all cases, with values for seismic loading being lower compared to static loading.
3. Results and Discussion
A total of 40 stability analyses were conducted for MSW landfills under both static and seismic loading conditions, along with an additional 191 sensitivity analyses. Additionally, ten reliability analyses were performed using the Monte Carlo method, involving a total of 20,000 simulations to determine the reliability index. Both the reliability and stability analyses were conducted for the same predefined slip surfaces where failure had occurred.
3.1. Stability of MSW Landfill Slopes under Static and Seismic Loading
The results of the stability analyses for static loading are presented in
Figure 7,
Figure 8,
Figure 9,
Figure 10,
Figure 11 and
Figure 12. The Xerolakka MSW landfill exhibits the lowest factor of safety (FS = 0.98), while the highest factor of safety is observed at the Maoershan MSW landfill (FS = 1.08). The higher factor of safety at the Maoershan MSW landfill can be attributed to the retaining structure located at the edge of the slope, which significantly enhances stability due to its shear strength and the gradual slope angle (14°, 1:4.05). Conversely, the lower factor of safety at the Xerolakka MSW landfill is primarily due to the very steep slope angle (31°, 1:1.67) and the very high leachate level, which increases pore pressures, thereby reducing shear strength on the slip surface due to decreased effective stress. The factors of safety for the other analyzed slopes range from 1.01 to 1.03, which are consistent with the observed slope failures in all the analyzed sites. The results of the slope stability analyses in all cases are in good agreement with the results from previous studies [
10,
15,
17,
18,
19,
20].
Table 12 presents a summary of the stability analyses conducted under both static and seismic loading conditions. The seismic analyses were performed using horizontal and vertical seismic coefficients of
kh = 0.1 and
kv = 0.05, respectively, which correspond approximately to an earthquake intensity of VI on the Modified Mercalli Scale. This intensity is characterized by the widespread perception of the earthquake, with many objects falling from shelves, plaster detaching from buildings, broken glass in windows, damaged chimneys, and localized incidents such as landslides, rockfalls, and liquefaction. The results of the analyses conducted under seismic loading indicate that even a relatively low-intensity earthquake is likely to trigger slope failure in the analyzed cases. Due to the significant reduction in the factor of safety observed across all analyses, no further analyses were conducted with higher seismic loading.
In all analyzed cases, there is a significant reduction in the factor of safety, ranging from 17% to 33%. The greatest reductions, 33%, and 31%, were observed in the cases of the Xiaping and Maoershan MSW landfills in China. These cases involve a combination of relatively large sliding mass areas and high leachate levels, which contribute to increased pore pressures. In contrast, the MSW landfill in Cincinnati has the largest sliding mass area among all the cases considered, but its relatively low leachate level results in a smaller reduction in the factor of safety due to seismic loading. In the case of Meethotamula, the low leachate level leads to a smaller reduction in the factor of safety, which is 17%. In the Xerolakka case, the reduction in the factor of safety due to seismic activity is smaller because of the relatively small sliding mass area. The results of the analyses indicate that a combination of a high leachate level and a large sliding mass area leads to the greatest reduction in the factor of safety due to seismic loading.
The significant decrease in the factor of safety under seismic conditions highlights the vulnerability of MSW landfills to earthquake-induced slope failures. These findings underscore the importance of comprehensive seismic hazard assessments and the implementation of risk mitigation measures for MSW landfills, even in regions where only low-intensity earthquakes are expected.
3.2. Stability of MSW Landfill Slopes under Static and Seismic Loading According to Eurocode 7
The summary of slope stability analyses conducted according to Eurocode 7 Design Approach 1, Combination 1 (DA1, C1), Design Approach 2 (DA2), and Design Approach 3 (DA3) for both static and seismic loading is presented in
Table 13,
Table 14 and
Table 15. Slope stability is quantified by the Over-Design Factor (ODF), which must be greater than or equal to one for the slope to satisfy the ultimate limit state criteria according to Eurocode 7. As shown in
Table 13, the stability analysis results for DA1 and C1 fully coincide with those of the classical stability analysis presented in
Table 12. This alignment is due to the specific combination of partial factors (
Table 7,
Table 8 and
Table 9), wherein DA1 and C1, the partial factors for materials and resistances, are set to one. Although the partial factor for external actions is 1.35, since there are no external actions on the slope, the results completely match those of the classical stability analysis.
Table 14 presents the summary of stability analysis results conducted according to Eurocode 7, Design Approach 2 (DA2). The results indicate that, under both static and seismic loading, none of the MSW landfills meet the ultimate limit state criteria specified by Eurocode 7, DA2. For all analyzed MSW landfills, ODF values are approximately 10% lower under both static and seismic loading compared to those in DA1 and C1. Additionally, the differences between the ODF values for static and seismic loading in DA2 are nearly identical to those observed in the previously analyzed design approach.
Table 15 presents the summary of stability analysis results conducted according to Eurocode 7, Design Approach 3 (DA3). The ODF values are further reduced by approximately 10% compared to DA2, and none of the analyzed MSW landfills meet the ultimate limit state criteria according to DA3. The differences between the ODF values for static and seismic loading are nearly identical to those observed in the previous design approaches.
According to Design Approach 1, Combination 1, all considered MSW landfills, except for the one in Xerolakka, meet the ultimate limit state criteria specified by Eurocode 7. However, none of the MSW landfills meet the ultimate limit state criteria under Design Approach 2 or Design Approach 3. Seismic loading with an intensity of VI on the Modified Mercalli intensity scale reduces the ODF by 16–33% compared to static loading.
The results of the stability analyses conducted in accordance with Eurocode 7 indicate that Design Approach 3 (DA3) yields the most conservative designs. This outcome arises from the application of partial factors from group M2 (
Table 8) to the shear strength parameters, whereas the other analyzed design approaches employ values from group M1. Consequently, this results in lower values for resisting forces (R
H,d, R
V,d) and moments (M
R,d), thereby reducing the shear strength of the MSW along the slip surface. In the absence of external actions, identical partial factors are applied to the design values of forces (V
H,d, V
V,d) and moments (M
Sd), contributing to sliding across all design approaches. This leads to the most conservative ODF value in the case of DA3. Design Approach 2 produces ODF values that are approximately 10% less conservative than those of DA3. This reduction is attributable to the smaller values of partial factors applied to the MSW shear strength parameters from group M2. Although DA2 uses the same partial factors for MSW shear strength parameters as DA1, C1, the ODF values in DA2 are about 10% lower. This difference is because DA2 applies partial factors for resistances from group R2, which are set at 1.1, while DA1, C1 uses values from group R1, set at 1.0. From these analyses, it becomes evident that the application of partial factors for MSW shear strength parameters from group M2 results in a 20% reduction in the ODF value compared to group M1. Similarly, the application of partial factors for resistances from group R2 reduces the ODF value by 10% compared to groups R1 and R3.
Given the potential consequences of MSW landfill slope failure and the significant heterogeneity of MSW, the findings of this study suggest that DA3 provides the safest and most robust designs.
However, a notable limitation of applying Eurocode 7 to the design of MSW landfills is the uniform value of partial factors for unit weight, set at 1.0 for both groups M1 and M2. The results of this study indicate that unit weight significantly influences the stability of MSW landfills. Therefore, it is recommended that further research be conducted to determine the appropriate values of partial factors for MSW unit weight, which, when used in conjunction with other partial factors, will ensure the desired reliability of MSW landfill structures.
3.3. Reliability of MSW Landfill Slopes under Static and Seismic Loading
Table 16 presents the summary of reliability analysis results.
The Monte Carlo method converged after 2000 simulations, which is consistent with the findings of other studies [
1,
27]. Reliability is quantified by the reliability index (β) and failure probability (Pf). For static loading, the Xiaping and Xerolakka landfills have the lowest reliability index values, β = 1, while Maoershan has the highest, β = 2.8. The reliability index values for the Cincinnati and Meethotamula landfills fall between these two extremes. The minimum target reliability index values according to Eurocode 7 for a reference period of 50 years are 3.8 for consequence class CC2 and 4.3 for CC3 (
Table 10), corresponding to failure probabilities of 7.23 × 10
−5 and 8.54 × 10
−6, respectively.
Table 16 shows that the reliability indexes and failure probabilities of all analyzed MSW landfills do not meet the minimum target values for consequence classes CC2 and CC3, which are typically applied to MSW landfills. According to the US Army Corps of Engineers (
Table 11), the expected performance level for the Maoershan landfill under static loading is below average; for Meethotamula it is poor, and for Cincinnati, Xiaping, and Xerolakka it is unsatisfactory. For seismic loading, the expected performance level is hazardous for all analyzed MSW landfills.
Based on the conducted analyses, it can be concluded that the calculated reliability of the slopes for all analyzed MSW landfills is significantly below the levels prescribed by Eurocode and the U.S. Army Corps of Engineers. As a result, slope failure could have been considered a plausible outcome for these landfills.
From
Table 12,
Table 13,
Table 14,
Table 15 and
Table 16, it can be observed that in the considered cases, there is no correlation between the safety factors of the ODF and the reliability index. This outcome is expected and is due to the fact that the results of reliability analyses predominantly depend on the statistical characteristics (coefficient of variation, statistical distribution) of the strength parameters, while these characteristics have less impact on the results of stability analyses [
33]. The lack of correlation between the ODF and the reliability index suggests that traditional stability analyses may not fully capture the risk associated with landfill slope stability. While ODF provides a measure of safety based on deterministic factors, the reliability index offers a probabilistic assessment that considers the inherent variability and uncertainty in material properties. This discrepancy underscores the importance of incorporating reliability-based approaches in the design and evaluation of MSW landfills to better account for the potential for failure under varying conditions.
3.4. Sensitivity Analyses Results
Figure 12 shows the results of the sensitivity analyses. The vertical axis represents the factor of safety, while the horizontal axis represents the sensitivity range, which reflects the minimum and maximum values of the input parameters considered in the analysis. These minimum and maximum values are determined according to Equation (5). In all cases, there is a positive correlation between increases in unit weight, cohesion, friction angle, and the factor of safety. For the Xiaping and Xerolakka landfills, the friction angle has the greatest influence on the factor of safety, while cohesion has the greatest impact in the case of the Maoershan landfill. In the Cincinnati and Maoershan landfills, unit weight is the most significant factor influencing the safety factor.
The results of the sensitivity analyses are consistent with expectations for shear strength parameters (cohesion and friction angle). In contrast to slopes constructed in soil, where soil unit weight has a negligible effect on the factor of safety, sensitivity analyses conducted in this study show that unit weight has a more significant impact on the stability of MSW landfill slopes. This is attributed to the relatively high coefficient of variation of MSW unit weight compared to that of soil, as well as the significantly lower value of the unit weight of MSW relative to soil unit weight. Therefore, in designing safe and reliable MSW landfill slopes, it is crucial to pay particular attention to accurately determining the characteristic values of both the strength parameters and the unit weight of MSW.
The findings from these sensitivity analyses underscore the importance of accounting for the variability in MSW properties during the design phase. The significant influence of unit weight on slope stability suggests that even small errors in estimating MSW unit weight could lead to considerable differences in the predicted safety factor.
These observations also point to the potential benefits of incorporating probabilistic methods into the design and analysis of MSW landfills. By considering the variability and uncertainty of MSW properties in a more structured way, engineers can better anticipate the range of possible outcomes and design more resilient structures.
3.5. Limitations and Implications for Future Research
The stability and reliability analyses conducted in this study were performed on simplified two-dimensional geotechnical models using the Spencer method, which is part of the limit equilibrium methods group. This method was chosen because it considers both moment equilibrium and force equilibrium. However, a limitation of the limit equilibrium method is that it is based solely on the principles of statics, neglecting strains and displacements. Despite this drawback, its main advantage is that it quantifies the safety margin of the slope and is routinely used in geotechnical engineering practice [
13], making it a suitable choice for this study. Another possible limitation of this study is the omission of the potential impact of geosynthetics on the stability and reliability of MSW landfill slopes. The analyses were based on literature data, which did not provide information on the influence of geosynthetics on MSW slope failures. Additionally, the literature lacks details on whether the design of these slopes was conducted according to any specific building codes.
The study revealed that the safety and reliability margins of all analyzed landfills are significantly lower than those prescribed by Eurocode 7. Although Eurocode 7 is primarily intended for the design of geotechnical structures involving soil, it is also applied to structures involving MSW based on the authors’ experience. Further research is needed to adapt Eurocode 7 for the design of structures involving MSW, particularly in terms of calibrating partial factors for loads, material properties, and resistances. According to the results of the sensitivity analyses, special attention should be given to the unit weight of MSW, as its impact on the stability of MSW landfill slopes is significant.
4. Conclusions
In this study, stability and reliability analyses were conducted on the slopes of five MSW landfills that experienced failures. Actual slip surfaces, failure modes, material properties, and leachate levels were obtained from the literature. As there were no available data on geosynthetics used in these MSW landfills, their influence on stability was not considered. The main factors leading to failure included high leachate levels due to inadequate drainage systems, poor MSW compaction, and overly steep slopes. The factors of safety for the analyzed slopes are around unity, which is expected given that slope failures occurred.
The primary goal of the study was to assess the safety and reliability margins of these slopes and compare them with the values prescribed in Eurocode. Additionally, the slopes of the considered MSW landfills were classified according to the US Army Corps of Engineers based on the reliability index. The study found that reliability indexes in all considered cases for static and seismic loading did not meet the target values for consequence classes CC2 and CC3 as prescribed in Eurocode 0. According to the U.S. Army Corps of Engineers classification, the expected performance level is below average for the Maoershan landfill, poor for Meethotamulla, and unsatisfactory for Cincinnati, Xiaping, and Xerolakka. For seismic loading, the expected performance level is hazardous for all analyzed MSW landfills. The results of the reliability analysis indicate an increased probability of failure in all the cases considered.
The study showed that the safety margins for the ultimate limit state of all considered MSW landfill slopes are below those prescribed in Eurocode 7, Design Approaches 2 and 3. In the case of Design Approach 1, Combination 1, 4 out of 5 analyzed landfills meet these criteria. Therefore, it can be concluded that the application of Design Approach 1, Combination 1, for the design of MSW landfill slopes can result in potentially unstable slopes. Design Approach 3 produces the most conservative designs, while Design Approach 2 results in designs that are 10–15% less conservative.
Seismic loading, equivalent to an earthquake of intensity VI according to the Modified Mercalli intensity scale, reduced the FS and ODF values by 10% in all considered cases. In contrast, the reduction in the reliability index was more pronounced, ranging from 20% to 74%. Based on these findings, it can be concluded that seismic loading is a critical factor to consider in the design of MSW landfill slopes, as it can significantly influence their stability.
The sensitivity analysis results showed that, in addition to the shear strength parameters of MSW, unit weight has a significant influence on the stability of MSW landfill slopes, which is not typically the case for slopes constructed in soil. Therefore, it is essential to accurately determine the characteristic values of MSW unit weight in the design process, which can be achieved, for example, by applying statistical methods.
To enhance the design of MSW landfill slopes, it is recommended to develop a comprehensive design framework that incorporates both traditional and advanced analytical methods. This framework could include stability analyses for both static and seismic loading conditions, reliability-based design, and sensitivity analyses. Implementing such a framework will facilitate the construction of more robust MSW landfill designs, capable of withstanding potential operational issues, such as drainage system blockages, inadequate material compaction, and other deficiencies.