1. Introduction
Today, rapid population growth, urbanization, industrialization, and technological developments have revealed the problem of pollution. Solid waste, which constitutes a large part of the pollution problem, is a situation that needs to be managed and improved in all countries of the world. Solid wastes have polluted and continue to pollute ground and underground waters. After the damage caused by ground and groundwater pollution became visible, the importance of this type of pollution was understood, and laws and regulations were prepared to eliminate the pollution that had occurred. If waste is not stored in a way that causes minimum damage to the environment, it is inevitable that negative effects will occur in the future. Since it is not possible to completely destroy waste, the main objective should be to reduce its mass, that is, to reduce its volume and to store it regularly in a way that will allow for reuse with technological developments that may occur in the following years.
Slope stability analyses in MSW landfills are very important both to know how high it can be filled without any movement on the slopes, and to determine the angle of the slopes that will be reconstructed if old or abandoned MSW sites are reused as construction sites.
In general, stability problems on slopes in landfills are similar to movements seen on normal slopes. Due to the non-homogeneous materials and non-uniform distribution in a landfill, movement can be observed on the slopes ranging from surface ruptures to large landslides. Similar to all natural or artificial slopes, the stability of slopes in MSW landfills against slipping and collapse under various loads is analyzed by the limit equilibrium method (LEM). In this method, the equilibrium between resisting and sliding forces is investigated on a known or accepted critical slip surface. In slope stability analyses, a certain factor of safety (FS) is determined, and necessary analyses are made. The FS calculated as a result of these analyses reveals whether the slope is stable against failure or the degree of stability.
With the finite element method (FEM), which is another method for slope stability analysis, slopes can be analyzed under different soil conditions and load combinations. By designing the slopes in two or three dimensions and using the appropriate soil model, realistic stress and displacement values can be obtained. In addition, necessary analyses can be made in cases where the slope is strengthened with different materials. Two approaches are generally used when performing slope stability analyses with FEM. In the first approach, the gravity increase method, the gravitational acceleration is increased until slope failure. In the second approach, the strength reduction method, cohesion and internal friction angles are reduced until slope failure.
There are numerous studies in the literature on the stability of slopes in MSW landfills [
1,
2,
3,
4,
5,
6,
7,
8]. In these studies, the slope stability analyses were carried out using both conventional LEM and FEM.
Seed et al. [
1] analyzed the causes of slope failure at a landfill in California. According to the results obtained from the analyses performed in the study, the FS was determined to be between 0.85 and 1.25 at the failure time.
Eid et al. [
2] investigated the stability of MSW landfill slopes. According to the results obtained from laboratory tests, field experiments, and back-calculations, it has been shown that the shear strength parameters for landfills can be considered as 35° and 0–50 kpa for the internal friction angle and cohesion, respectively.
In the study carried out by Koerner and Soong [
3], ten landfill failures were discussed. According to the results obtained, numerical analyses for stability calculations depend more on shear strength than other parameters.
Gharabaghi et al. [
4] carried out slope stability analyses at two existing solid waste landfills in Brazil. As a result of the study, it was stated that the determination of suitable shear strength parameters and the solid waste composition, grain diameter, degree of deterioration, and moisture content in the field are very important in the analysis.
Chang [
5] carried out three-dimensional numerical analyses to model the slope failure occurring in the Kettleman Hills Landfill. It was shown that the internal friction angle and cohesion values obtained as a result of the back-analysis carried out for the post-failure slope and the values obtained by the laboratory experiments were compatible.
Stark et al. [
6] investigated the shear strength of MSW landfills by back-analysis of failed waste slopes using laboratory tests and field experiments. Using the results obtained from the study, suggestions were made for modeling the shear strength of landfills in analyses.
Hossain and Haque [
7] analyzed the stability of landfill slopes as a function of time and decomposition. In the study, numerical analyses were carried out using FEM and LEM. At the end of the study, they compared the obtained results from the numerical analyses.
In the study carried out by Huvaj-Sarihan and Stark [
8], the shear strength parameters of waste landfills were investigated using back-analysis of failed waste slopes. They presented the back analyses of the failed waste slopes in four landfills. Each of the landfill slope failures was reviewed, and the results of the back-analyses were presented. At the end of the study, the obtained and recommended parameters were compared with other studies.
Because the tensile strength of the soils is low, the slopes can be made more stable with materials such as metal reinforcements, geotextiles, or geogrids placed inside the slope. High-strength elements are used in these structures, which are defined as reinforced slopes. While metal strips or reinforcements were used in the first applications, today such materials have been replaced by geosynthetic materials. In slopes that are reinforced to increase stability, a solution is realized by including the geogrid forces in the limit equilibrium or finite element analysis as known forces. Geogrids are high-strength, polymer-structured, geosynthetic materials with sufficient space and different grid structures to ensure interlocking with the soil. In its chemical structure, polyethylene-, polyester-, or polypropylene-type polymers are available. Geogrids have an important role in the reinforcement of soils with their high tensile strength and are often used to reinforce slopes.
Although there are many studies in the literature on the stability of slopes in landfills, investigations on the stability behavior of reinforced landfill slopes are limited. Today, increasing the stability of slopes using geogrid materials has become a routine and standard practice. However, their use for the reinforcement of landfill slopes is rare. The few studies that have been carried out on reinforcing solid waste slopes using geogrids usually include field applications.
Zornberg and Kavazanjian [
9] investigated the integrity of geogrid-reinforced landfill slopes subjected to differential settlements and seismic loadings experimentally and numerically. According to the results obtained from the study, it showed that the critical reinforced region corresponding to different loading mechanisms occurred at different elevations within the reinforced soil structure.
Jiang et al. [
10] performed finite element analysis to investigate the failure mechanism and FS values of a geotextile-reinforced extended berm in a MSW landfill slope of 45°.
Ke Han et al. [
11] performed a one-year field monitoring of a MSW landfill slope reinforced by geogrids. Using the results obtained, suggestions were made for future studies on the subject.
While the reinforcement of MSW landfill slopes using geogrids is available at existing landfills, the interactions between the solid waste material and geogrids have not yet been clearly demonstrated. In particular, studies on field observations of the results obtained by numerical and theoretical studies are quite limited. Subjects such as the long-term behavior of geogrids, the effect of mechanical creep, and the friction between geogrids and MSWs still need to be investigated. Friction characteristics between synthetic materials and geotechnical materials also have an important impact on slope stability. Interface shear behavior of synthetic and geotechnical materials has been studied recently due to slope failures in landfills.
Bergado et al. [
12] investigated the properties of different interfaces in landfills using laboratory experiments.
In the study by Kim and Frost [
13], the geotextile/geomembrane interfacial shear behavior was investigated experimentally and generalized interfacial shear mechanisms were proposed based on the results obtained.
Cen et al. [
14] performed a study on the cyclic interface shear between geomembranes and sandy gravel. At the end of the study, a new empirical correlation was suggested.
In the experimental study carried out by Shi et al. [
15], shear tests were performed at the geomembrane/geotextile and geomembrane/geocomposite/sand liner interfaces where the normal stress changes at different shear stages. According to the results obtained from the study, it was stated that the effects of normal stress changes on the lining interface strength should be considered in the slope stability analysis of MSW landfills.
A study was carried out experimentally investigating the interface creep behavior between geomembranes and geotextiles in the lining system by Lu et al. [
16]. Conventional direct shear tests were conducted using different combinations of geomembranes and geotextiles. It was recommended to conduct longer-term geomembrane–geotextile interface creep tests in order to understand the interface creep mechanism in depth.
This study presents the slope stability analysis of MSW landfills carried out using finite element and limit equilibrium methods. Finite element analyses were performed using PLAXIS [
17], and limit equilibrium analyses were performed using the SLOPE/W [
18] computer programs. The stability behavior of the landfill slopes was analyzed for both unreinforced and geogrid-reinforced conditions. A landfill slope model was selected for the analysis, and model parameters were obtained from the studies available in the literature. In this study, the optimum value of the vertical distance between the geogrid reinforcements and the number of reinforcements, which increase the FS of the MSW slope, are obtained, and the necessary recommendations for practical applications are presented.
2. Materials and Methods
For the slope stability analyses in landfills, a model with a width of 210 m and a height of 45 m was created (
Figure 1). The geometry was reconstructed for different slope angles in the analyses. The symmetry of the geometric model was used for convenience in the analyses. Slope stability analyses of the model were performed with PLAXIS using FEM, and SLOPE/W using LEM. Unreinforced and reinforced analysis results obtained by both methods were compared.
In the analysis, the waste material and soil were modeled with the Mohr–Coulomb model (MCM). The input parameters were obtained from previous experimental studies. The most important issue in the slope stability analysis of MSW landfills is the realistic estimation of the parameters of the waste material including unit weight, internal friction angle, and cohesion. Due to the natural properties of solid waste materials, such as moisture content, degree of degradation, waste composition, and particle size, it is very difficult to determine these parameters in the laboratory. When the existing studies in the literature are examined, it is seen that the mentioned shear strength parameters and unit weight values are in a wide range.
Many investigations have been performed on the geotechnical properties of landfills to examine the settlement and stability of landfills [
19,
20,
21,
22,
23,
24].
Average in-place unit weights, γ, used by owners and operators for landfill capacity estimates are typically in the range of 8.6 to 10.2 kN/m
3. Values in this range have also been used for seismic analyses by [
25,
26,
27]. Fassett et al. [
28], in their study, presented a summary of the unit weight values obtained in situ by various researchers for MSW. Accordingly, the unit weight values were between 2.9 kN/m
3 to 14.4 kN/m
3.
The shear strength values of the MSW landfills presented in the studies carried out on the subject are very variable, with internal friction angles, ϕ, as low as 10 degrees and as high as 53 degrees and cohesion values, c, varying from 0 to 67 kPa [
29].
In the study performed by Sharma et al. [
30] in an existing MSW landfill, Poisson’s ratio value was obtained as 0.49. Karimpour-Fard and Machado [
31] carried out an experimental study to estimate the deformation properties of MSW materials. The results obtained from this study show that Poisson’s ratio value varies between 0.35 and 0.49.
The modulus of elasticity, E, of waste materials is not constant and depends on the average stress [
28,
32]. Therefore, the modulus of elasticity value is in a wide range.
In the stability analysis, the material parameters presented in
Table 1, obtained as a result of the literature study, were used for the modeling of the solid waste material and soil.
The geometric model of the MSW landfill was created in the PLAXIS computer program as two-dimensional. Slope analyses were performed according to plane strain conditions. In the modeling, triangular elements with 15 nodes were chosen to obtain more accurate and sensitive results.
In the analyses, the standard boundary conditions available in the PLAXIS program were selected. Vertical and horizontal displacements at the base of the geometric model (ux = 0, uy = 0) were assumed to be zero, while on the vertical side of the slope, only vertical displacements were taken into account (ux = 0, uy = free).
In the analyses, geogrid elements were used to model the geogrid reinforcement layers. The EA value of 1100 kN/m was entered into the program as the material property of the geogrid reinforcement layers.
While creating the finite element mesh, the most suitable mesh structure that gives the most accurate result in the optimum time where the analysis results do not change much was investigated. For this purpose, the effect of the mesh structure on the FS was investigated by performing analyses for medium, fine, and very fine mesh. As a result of the analyses carried out, no significant change was observed in the results in the case of the fine mesh, and the model was created in a fine mesh structure.
In the unreinforced case, the analyses were carried out using two phases. In the first phase, the initial stresses due to the soil self-weight were created, and in the second stage, the phi-c reduction analysis was performed. In the reinforced case, the geogrid reinforcement was activated, and the analyses were carried out. For the finite element analysis, the slope geometries considered in the reinforced and unreinforced conditions are shown in
Figure 2.
In this study, limit equilibrium analyses were performed using the SLOPE/W computer program. With the LEM, the behavior of the soils was modeled based on certain assumptions, and geotechnical problems could be solved according to these assumptions. In this method, assuming that the slip circle was formed on a certain surface, FS values were obtained by comparing the forces sliding along this surface with the forces resisting the slip. With these analysis methods, the static equilibrium equation of the slope was attempted to be obtained using the Mohr–Coulomb stress criteria. There are three main types of limit equilibrium analyses used in slope stability calculations. These are: the slice method, the wedge method, and the infinite slope method. The most widely used method in slope stability analysis is the slice method. The slice methods used in practice are based on dividing the slip surface into sufficient vertical slices. Many approaches have been developed using the slice method, such as those by Bishop [
33], Janbu [
34], Spencer [
35], and Morgenstern and Price [
36]. Although there are some differences in practice between these methods, the common feature is the investigation of the slip mass equilibrium on a known or accepted critical slip surface. The slice method was used in the analyses performed in this study. The Bishop and half-sine function methods were selected in the calculations.
Figure 3 shows the geometry of the unreinforced and reinforced models considered in the limit equilibrium analysis.
The Mohr–Coulomb model (MCM) was used as the material model in the analysis. The values in
Table 1 were used for the material parameters required in this model. After the material properties were entered into the program, the slip surface model was defined. There are two types of failures in the SLOPE/W program: circular and block slip. The circular slip surface model was used in the slope analysis. Although there are many methods to determine the circular slip surface, the entry and exit method was used. The reason for using this method was that the boundaries of the slip surface could be determined clearly.