Slope Stability Research Articles

An investigation on anisotropic soil slope stability by LS-SVM and LEM approaches

Purpose This study aims to use regression Least-Square Support Vector Machine (LS-SVM) as a probabilistic model to determine the factor of safety (FS) and probability of failure (PF) of anisotropic soil slopes. Design/methodology/approach This research uses machine learning (ML) techniques to predict soil slope failure. Due to the lack of analytical solutions for measuring FS and PF, it is more convenient to use surrogate models like probabilistic modeling, which is suitable for performing repetitive calculations to compute the effect of uncertainty on the anisotropic soil slope stability. The study first uses the Limit Equilibrium Method (LEM) based on a probabilistic evaluation over the Latin Hypercube Sampling (LHS) technique for two anisotropic soil slope profiles to assess FS and PF. Then, using one of the supervised methods of ML named LS-SVM, the outcomes (FS and PF) were compared to evaluate the efficiency of the LS-SVM method in predicting the stability of such complex soil slope profiles. Findings This method increases the computational performance of low-probability analysis significantly. The compared results by FS-PF plots show that the proposed method is valuable for analyzing complex slopes under different probabilistic distributions. Accordingly, to obtain a precise estimate of slope stability, all layers must be included in the probabilistic modeling in the LS-SVM method. Originality/value Combining LS-SVM and LEM offers a unique and innovative approach to address the anisotropic behavior of soil slope stability analysis. The initiative part of this paper is to evaluate the stability of an anisotropic soil slope based on one ML method, the Least-Square Support Vector Machine (LS-SVM). The soil slope is defined as complex because there are uncertainties in the slope profile characteristics transformed to LS-SVM. Consequently, several input parameters are effective in finding FS and PF as output parameters.

Numerical and experimental pseudo-static study on block arrangement effects in traditional dry-stone walls for out-of-plane mechanical behavior evaluation

In Lima’s slopes, the absence of urban planning has led to housing built on backfills contained by traditional dry-stone retaining walls, called “pircas,” which pose high seismic risk. A previous study [1], which combined numerical and experimental analyses, evaluated pircas through pseudo-static tests that reflected common construction methods in Lima’s eastern region. One recommendation from that study was to explore effective stone arrangements to enhance the seismic performance of pircas. Consequently, this study aims to examine how the construction process, particularly the arrangement of stones, impacts the out-of-plane behavior of pircas. The research focuses on identifying construction-feasible configurations that can effectively improve their seismic performance. For this purpose, natural-scale pseudo-static experimental studies and numerical modeling using the Discrete Element Method (DEM) were carried out. The experimental campaign consisted of nine pseudo-static tests on different wall arrangements. The pseudo-static numerical study included seventy wall arrangements. The geometry of the blocks was simplified to eliminate this component from the analysis, which focused on the arrangement of the wall. The numerical and experimental results showed the same trend. The study has demonstrated that the current construction technique can be significantly improved by regularly incorporating through-stones—tie stones that pass through the entire cross-section of the wall—and considering a cross-section with overlap. This approach leads to an increase of up to 50% in lateral strength compared to current practices, verifying the effectiveness of through-stones. The study does not intend to endorse the use of pirca for urban habitation in high seismic sectors. However, these traditional walls can be recommended for slope stability, environmental mitigation, urban agriculture, and protection of natural areas.

Effect of geometry on the natural frequency and seismic response characteristics of slopes subjected to pulse-like ground motions

Near-fault regions are particularly vulnerable to seismic-induced landslides due to the intense energy pulses in near-fault ground motions (NFGMs). These pulses, shaped by terrain geometry and material properties, significantly influence seismic response and slope stability. This study investigates the impact of slope geometry on natural frequency and seismic response characteristics under both pulse-like ground motions (PLGMs) and non-pulse ground motions (non-PGMs). The results show that increasing slope height lowers natural frequency, making the slope more susceptible to resonance with seismic waves, thus amplifying ground motion and increasing instability. Similarly, steeper slopes also reduces the natural frequency, heightening instability by up to 0.17%. PLGMs generate seismic responses approximately 7% stronger than those induced by non-PLGMs. Furthermore, as the frequency of PLGMs rises, so does their destructive potential. Material analysis reveals that Rock Class A has a natural frequency 68% higher than Rock Class D, making it significantly resistant to seismic deformation. These insights are essential for designing more resilient slopes in seismic-prone regions.

INNOVATIVE EXPERIMENTAL TESTING PROGRAM OF DIRECT SHEAR TEST IN SOIL MECHANICS

The research work aims at analyzing for the first time the data set obtained on cohesive soil samples following the publication of the Romanian Invention Patent RO 134239. The standard test method for the direct shear test provides the shear strength parameter – internal friction angle in consolidated drained condition - of either undisturbed or remolded soil samples forcing the shear plane at the midsection of the sample in the horizontal direction. The samples are provided in parallelepipedal shape (6 cm x 6 cm x 2 cm) and the displacement rate in horizontal direction is 0.1 mm/min. The new equipment patented in Romania changes the direction of shearing, from horizontal to vertical, and the soil samples are of cubic shape (6 cm x 6 cm x 6 cm). The experimental program involves testing both the parallelepipedal and cubic samples using the same motorized mechanism, with simultaneous readings from their respective micro-comparators. The UU test is performed without allowing consolidation and shearing at 1.0 mm/min. For the CD test, samples are consolidated under vertical loads for 24 hours before shearing at 0.1 mm/min. The shear stresses for cubic samples were higher than those for parallelepipedal samples, with residual stresses reflecting this trend. For cubic samples, both the peak and residual shear stresses trend lines indicated higher cohesion (c) and lower internal friction angle (𝜙) for UU tests and CD tests in contrast to parallelepipedal samples in both testing conditions. The innovative testing program allows for variability in shear strength parameters along the soil failure surface in both natural and compacted soil structures. This differentiation divides the soil condition into drained and undrained states at the initiation, emergence points, and the point of maximum depth along the failure surface. This approach is significant for accurately assessing soil shear resistance and potential failure mechanisms. The study's findings suggest a nuanced approach to parameter selection for slope stability analysis, ensuring accurate representation of both cohesion and internal friction in stability models.

Slope Stability Analysis of Open-Pit Mine Considering Weathering Effects

Weathering processes gradually alter the physical and mechanical attributes of slope materials, weakening the structural integrity and stability of slopes. This paper presents an in-depth analysis of slope stability in an open-pit mine, emphasizing the pivotal role of weathering effects in determining slope stability. To accurately capture the impact of weathering on slope stability, a comprehensive analysis model was developed, incorporating field observations, laboratory testing, and numerical simulations. The effects of weathering on the mechanical properties of black shale were studied through extensive laboratory tests. The uniaxial compressive strength, shear strength, and modulus of elasticity significantly decreased with increasing weathering, indicating a heightened vulnerability to slope failure. The correlation function between mechanical parameters and weathering time was obtained, providing the basis for evaluating the stability of mine slopes. It was found that more severe weathering conditions were strongly correlated with elevated risks of slope failure, including landslides and collapses. Based on these findings, practical recommendations are provided for slope reinforcement and management strategies, aimed at mitigating slope failure risks and ensuring the safe and efficient operation of the mine. By incorporating weathering effects into slope stability analysis, mine operators can make informed decisions that account for the dynamic nature of slope materials and their susceptibility to weathering, thereby improving overall mine performance and safety.

A scientometrics review of conventional and soft computing methods in the slope stability analysis

Predicting slope stability is important for preventing and mitigating landslide disasters. This paper examines the existing approaches for analyzing slope stability. There are several established conventional approaches for slope stability analysis that can be applied in this context. However, in recent decades, soft computing methods has been extensively developed and employed in stochastic slope stability analysis, notably as surrogate models to improve computing efficiency in contrast to traditional approaches. Soft computing methods can deal with uncertainty and imprecision, which may be quantified using performance indices like coefficient of determination, in regression and accuracy in classification. This review study focuses on conventional methods such as the Bishop’s method and Janbu’s method, as well as soft computing models such as support vector machine, artificial neural network, Gaussian process regression, decision tree, etc. The advantages and limitations of soft computing techniques in relation to conventional methods have also been thoroughly covered in this paper. The achievements of soft computing methods are summarized from two aspects—predicting factor of safety and classification of slope stability. Key potential research challenges and future prospects are also given.

Development of Soft Computing-based Predictive Tools for Estimating the Young Modulus of Weak Rocks

The deformation characteristics of rocks are of vital importance in addressing most geomechanical issues as they are one of the most critical input parameters in rock engineering analyses. For this reason, robust forecasting models are required when analysing the stability of tunnels, slopes, mine galleries, and other underground excavations. In this research, novel predictive models are proposed to estimate the tangential Young modulus (Eti) of weak rocks. To achieve this, an extensive literature review is performed to obtain a comprehensive database including critical physico-mechanical properties of various weak rocks. Thanks to the advantages of soft computing methods such as genetic algorithm (GA), adaptive neuro-fuzzy inference system (ANFIS), artificial neural networks (ANN) and multivariate adaptive regression splines (MARS), novel predictive models are established. The effectiveness of the developed predictive models is investigated using various statistical measures and it is concluded that empirical models utilizing ANN and ANFIS methodologies are the most effective tools for estimating the Eti of weak rocks. In addition, a practical design chart is also developed for assessing the Eti of weak rocks.

Bayesian ensemble learning and Shapley additive explanations for fast estimation of slope stability with a physics-informed database

Bayesian ensemble learning and Shapley additive explanations for fast estimation of slope stability with a physics-informed database

Enhancing deep learning-based slope stability classification using a novel metaheuristic optimization algorithm for feature selection

The evaluation of slope stability is of crucial importance in geotechnical engineering and has significant implications for infrastructure safety, natural hazard mitigation, and environmental protection. This study aimed to identify the most influential factors affecting slope stability and evaluate the performance of various machine learning models for classifying slope stability. Through correlation analysis and feature importance evaluation using a random forest regressor, cohesion, unit weight, slope height, and friction angle were identified as the most critical parameters influencing slope stability. This research assessed the effectiveness of machine learning techniques combined with modern feature selection algorithms and conventional feature analysis methods. The performance of deep learning models, including recurrent neural networks (RNNs), long short-term memory (LSTM) networks, and generative adversarial networks (GANs), in slope stability classification was evaluated. The GAN model demonstrated superior performance, achieving the highest overall accuracy of 0.913 and the highest area under the ROC curve (AUC) of 0.9285. Integration of the binary bGGO technique for feature selection with the GAN model led to significant improvements in classification performance, with the bGGO-GAN model showing enhanced sensitivity, positive predictive value, negative predictive value, and F1 score compared to the classical GAN model. The bGGO-GAN model achieved 95% accuracy on a substantial dataset of 627 samples, demonstrating competitive performance against other models in the literature while offering strong generalizability. This study highlights the potential of advanced machine learning techniques and feature selection methods for improving slope stability classification and provides valuable insights for geotechnical engineering applications.

Engineering geological investigation of Gololcha dam for evaluation of leakage and abutment slope stability, Eastern Ethiopia

The Ethiopian economy is primarily dependent on agriculture, making the construction of water harvesting facilities, such as dams, crucial for improving the productivity of this sector. The ongoing construction of the Gololcha dam on the Kurkura River located in Eastern Ethiopia aims to enhance irrigation schemes in the region. However, the dam site's complex geological and structural conditions pose challenges related to leakage and slope instability. Hence, this study focuses on addressing the abutment slope stability and leakage condition of this dam. This study employed kinematic analysis and the Limit Equilibrium Method (LEM) to assess slope stability. Additionally, engineering geological mapping, discontinuity surveys, seismic refraction tomography (SRT), and in-situ permeability testing were used to evaluate the leakage condition of the dam site. Notably, the permeability and SRT survey results identified potential leakage zones to the depth of 35, 30, and 35 m at the left, right, and central foundations of the dam, respectively. The kinematic method revealed one planar and two wedge modes of failure in the slope section covered by slightly weathered and fractured basalt rock at the right abutment. Further stability analysis of these two modes of failures via LEM analysis indicated slope instability under saturated conditions, emphasizing the role of pore water pressure. Furthermore, LEM modeling was directly utilized using the Slide 6.0 software to analyze the slope stability condition of the left abutment of the dam. This modeling also uncovered instability under saturated conditions. Based on the study findings, this study recommended curtain grouting to address potential leakage, as well as slope flattening and removing unstable rock wedges and loose material to stabilize unstable slope sections.

The Drawdown of a Reservoir: Its Effect on Seepage Conditions and Stability of Earth Dams

This article addresses the reliability and safety of an earth dam in the case of a change in the reservoir water level. The water level must often be reduced to remove water or as a response to an emergency situation in the process of operation of a hydraulic structure. Lower water levels change seepage conditions, such as the surface of depression, values and directions of seepage gradients, seepage rates, and volumetric hydrodynamic loading. Practical hydraulic engineering shows that these changes can have a number of negative consequences. Higher seepage gradients can lead to seepage-triggered deformations in the vicinity of the upstream slope of a structure. Hydrodynamic loads, arising during drawdown, reduce the stability of an upstream slope of a dam and cause its failure. Potential consequences of a drawdown can be evaluated by solving the problem of drawdown seepage for the dam body and base. A numerical solution to this problem is based on the finite element method applied using the PLAXIS 2D software package. Results thus obtained are compared with those obtained using the finite element method in the locally variational formulation. A numerical experiment was conducted to analyze factors affecting the value of the maximum seepage gradient and stability of the earth dam slope. Recommendations were formulated to limit the drawdown parameters and to ensure the safe operation of a structure.

Landslide Analysis with Incomplete Data: A Framework for Critical Parameter Estimation

Landslides are one of the most common geohazards, posing significant risks to infrastructure, recreation, and human life. Slope stability analyses rely on detailed data, accurate materials testing, and careful model parameter selection. These factors are not always readily available, and estimations must be made, introducing uncertainty and error to the final slope stability analysis results. The most critical slope stability parameters that are often missing or incompletely constrained include slope topography, depth to water table, depth to failure plane, and material property parameters. Though estimation of these values is common practice, there is limited guidance or best practice instruction for this important step in the analysis. Guidance is provided for the estimation of: original and/or post-failure slope topography via traditional methods as well as the use of open-source digital elevation models, water table depth across variable hydrologic settings, and the iterative estimation of depth to failure plane and slope material properties. Workflows are proposed for the systematic estimation of critical parameters based primarily on slide type and scale. The efficacy of the proposed estimation techniques, uncertainty quantification, and final parameter estimation protocol for data-sparse landslide analysis is demonstrated via application at a landslide in Colorado, USA.

Investigation on the Instability Mechanism of Expansive Soil Slope With Weak Interlayer Based on Strain Softening

ABSTRACTExpansive soils are widespread in the world and coincide with areas of high human activity. The main cause of deep instability of expansive soil slopes is due to their softening caused by excavation and seepage. By developing a comprehensive numerical model based on the theory of unsaturated soil, this study examines the characteristics of stress and displacement distribution of expansive soil slopes through hydraulic‐mechanical coupled numerical simulation. This study analyzes the evolution patterns of slopes with excavation unloading and seepage of water storage to reveal the mechanisms of deep‐seated instability of expansive soil slopes. The findings demonstrate that: The instability of expansive soil slopes begins at the foot of the slope and propagates along the interlayer, affecting the entire slope. Excavation leads to the softening of the expansive soil interlayer and the transfer of shear stress. During water storage, the weakening of the soil strength results in slope instability along the weak interlayer slip. Softening of the expansive soil interlayer facilitates the redistribution of shear forces in the slope and alters the distribution law of the plastic zone in the deep layer. Overly slowing down the slope leads to significant excavation unloading, which is detrimental to the slope's stability.

Discrete-element method simulations of shallow landslides triggered by rainfall

Rainfall infiltration is a primary cause of slope failure. Analyzing the entire process of a landslide triggered by rainfall is crucial for the design of infrastructure. However, the complexity associated with simulating large deformations has often hindered previous studies from modeling the entire landslide process subsequent to rainfall infiltration. This study aims to provide a method to simulate the entire process of landslides considering rainfall and more insight into the mechanism of landslides. The discrete element method (DEM) is first applied to model the entire process of shallow landslides triggered by rainfall. The numerical results indicate that the intensity and duration of rainfall infiltration have significant impacts on the slope stability. The particle displacement field, resulting from landslides due to rainfall, corresponds with three distinct stages within the landslide process. The increased rainfall intensity causes the area of the particles within failure mass grows rapidly during the initial stage of the landslide. The closer to the slope surface the faster the particle velocity. The simulations demonstrate that the slope surface transitions from a linear to a concave downward geometry during the landslides. Larger particles display greater displacement and velocity compared to smaller particles.

Experimental Investigation on Shear Strength at the Permeable Concrete–Fine-Grained Soil Interface for Slope Stabilization Using Deep Socket Counterfort Drains

In slopes where high pore water pressure exists, deep counterfort drains (also called drainage trenches or trench drains) represent one of the most effective methods for improving stability or mitigating landslide risks. In the cases of deep or very deep slip surfaces, this method represents the only possible intervention. Trench drains can be realized by using panels or secant piles filled with coarse granular material or permeable concrete. If the trenches are adequately “socket” into the stable ground (for example sufficiently below the sliding surface of a landslide or below the critical slip surface of marginally stable slopes) and the filling material has sufficient shear strength and stiffness, like porous concrete, there is a further increase in shear strength due to the “shear keys” effect. The increase in shear strength is due both to the intrinsic resistance of the concrete on the sliding surface and the resistance at the concrete–soil interface (on the lateral surface of the trench). The latter can be very significant in relation to the thickness of the sliding mass, the “socket depth”, and the spacing between the trenches. The increase in shear strength linked to the “shear keys effect” depends on the state of the porous concrete–soil interface. For silty–clayey base soils, it is very significant and is of the same order of magnitude as the increase in shear resistance linked to the permanent reduction on the slip surface in pore water pressure (draining effect). This paper presents the results of an experimental investigation on the shear strength at the porous interface of concrete and fine-grained soils and demonstrates the high significance and effectiveness of the “shear keys” effect.

An Analytical Insight Into Stability Analysis of Unsaturated Multi‐Layered Slopes Subjected to Rainfall Infiltration

ABSTRACTSlopes in nature usually present layered characteristics, and its stability is susceptible to rainfall events. Considering that current analytical solutions are only suited to simulate the rainfall infiltration of double‐layered infinite unsaturated slopes, an analytical procedure is hence proposed in this study to tackle the consideration of multiple layers. The variable separation method and transfer matrix method are combined to derive the analytical solution of pore water pressure (PWP) for simulating rainfall infiltration in layered infinite unsaturated slopes. After having validated the proposed model and analytical solutions by comparing with existing literature and numerical simulation, the closed‐form solution of PWP is incorporated into the limit equilibrium for assessing slope stability. A three‐layer slope is selected as an example for further discussion. PWP distribution and factor of safety are calculated, considering the effects of saturated hydraulic conductivity and thickness of the upper layer, intensity of antecedent and subsequent rainfall, and varied soil unit weight along the depth. The slope stability subjected to rainfall effects is consistent with the variation of PWP. The proposed analytical solutions provide a simple and practical avenue for computing PWP distribution and evaluating the stability of multi‐layered slopes under rainfall conditions, which can also serve as a benchmark for numerical solutions.

Effect of slope angle on fractured rock masses under combined influence of variable rainfall infiltration and excavation unloading

Intense precipitation infiltration and intricate excavation processes are crucial factors that impact the stability and security of towering and steep rock slopes within mining sites. The primary aim of this research was to investigate the progression of cumulative failure within a cracked rock formation, considering the combined effects of precipitation and excavation activities. The study was conducted in the Huangniuqian eastern mining area of the Dexing Copper Mine in Jiangxi Province, China. An engineering geological investigation was conducted, a physical model experiment was performed, numerical calculations and theoretical analysis were conducted using the matrix discrete element method (MatDEM), and the deformation characteristics and the effect of the slope angle of a fractured rock mass under different scenarios were examined. The failure and instability mechanisms of the fractured rock mass under three slope angle models were analyzed. The experimental results indicate that as the slope angle increases, the combined effect of rainfall infiltration and excavation unloading is reduced. A novel approach to simulating unsaturated seepage in a rock mass, based on the van Genuchten model (VGM), has been developed. Compared to the vertical displacement observed in a similar physical experiment, the average relative errors associated with the slope angles of 45°, 50°, and 55° were 2.094%, 1.916%, and 2.328%, respectively. Accordingly, the combined effect of rainfall and excavation was determined using the proposed method. Moreover, the accuracy of the numerical simulation was validated. The findings contribute to the seepage field in a meaningful way, offering insight that can inform and enhance existing methods and theories for research on the underlying mechanism of ultra-high and steep rock slope instability, which can inform the development of more effective risk management strategies.

How water, temperature, and seismicity control the preconditioning of massive rock slope failure (Hochvogel)

Abstract. The anticipation of massive rock slope failures is a key mitigation strategy in a changing climate and environment requiring a precise understanding of pre-failure process dynamics. Here we exploit >4 years of multi-method high-resolution monitoring data from a large rock slope instability close to failure. To quantify and understand the effect of possible drivers (water from rain and snowmelt, internal rock fracturing, and earthquakes), we correlate slope displacements with environmental data, local seismic recordings, and earthquake catalogues. During the snowmelt phase, displacements are controlled by meltwater infiltration with high correlation and a time lag of 4–9 d. During the snow-free summer, rainfall induces accelerations with a time lag of 1–16 h for up to several days without a minimum activation rain sum threshold. Rock fracturing, linked to temperature and freeze–thaw cycles, is predominantly near the surface and unrelated to displacement rates. A classic Newmark analysis of recent and historic earthquakes indicates a low potential for immediate triggering of a major failure at the case site, unless it is already very close to failure. Seismic topographic amplification of the peak ground velocity (PGV) at the summit ranges from a factor of 2–11 and is spatially heterogeneous, indicating a high criticality of the slope. The presented in-depth monitoring data analysis enables a comprehensive rockfall driver evaluation and indicates where future climatic changes, e.g. in precipitation intensity and frequency, may alter the preconditioning of major rock slope failures.

The Uses of Lightweight Material in Civil Engineering : A Review

A lightweight aggregate (LWA) is a group of aggregates with a relative density lower than normal density aggregates (natural sand, gravel, and crushed stone), sometimes referred to as low density aggregate. It is utilized for its reduced weight and exceptional qualities in sound reduction, fire resistance, insulation, and geotechnical applications LWA, or lightweight aggregate, is a cutting-edge construction material that is employed to decrease the overall weight of tall structures. Despite being lighter, it possesses a strength that is comparable to that of regular concrete. In recent times, there has been a significant surge in the use of lightweight aggregate in the field of building, mostly because of its exceptional performance in seismic situations. In addition, substituting natural aggregate with other industrial by-products and trash enables us to minimize the adverse environmental consequences. Laser wavelength ablation (LWA) has also been employed to address a wide range of geotechnical engineering challenges, including the transformation of soft and unstable soil into a suitable and usable property. Lightweight fillings have a weight that is around 50% less than fills made using typical materials. The load reduction, along with the high internal friction angle of the lightweight aggregate, can decrease vertical and lateral pressures by almost 50%, potentially resulting in less settling. This literature study examines the use of lightweight materials, including LECA, Bonza, and Thermostone, in geotechnical engineering applications. The results demonstrate a significant emphasis on improving slope stability and minimizing stresses on susceptible soils. Nevertheless, several knowledge gaps exist regarding the efficacy of these materials under challenging conditions and when subjected to dynamic loads. A growing inclination towards using these lightweight aggregates is observed because of their advantageous environmental and mechanical properties. Further investigation is needed to examine the extended-term performance and interactions of these materials with different soil types in order to enhance their performance in diverse geotechnical situations.

Cascading Landslide: Kinematic and Finite Element Method Analysis through Remote Sensing Techniques

Cascading landslides are specific multi-hazard events in which a primary movement triggers successive landslide processes. Areas with dynamic and quickly changing environments are more prone to this type of phenomena. Both the kind and the evolution velocity of a landslide depends on the materials involved. Indeed, rockfalls are generated when rocks fall from a very steep slope, while debris flow and/or mudslides are generated by fine materials like silt and clay after strong water imbibition. These events can amplify the damage caused by the initial trigger and propagate instability along a slope, often resulting in significant environmental and societal impacts. The Morino-Rendinara cascading landslide, situated in the Ernici Mountains along the border of the Abruzzo and Lazio regions (Italy), serves as a notable example of the complexities and devastating consequences associated with such events. In March 2021, a substantial debris flow event obstructed the Liri River, marking the latest step in a series of landslide events. Conventional techniques such as geomorphological observations and geological surveys may not provide exhaustive information to explain the landslide phenomena in progress. For this reason, UAV image acquisition, InSAR interferometry, and pixel offset analysis can be used to improve the knowledge of the mechanism and kinematics of landslide events. In this work, the interferometric data ranged from 3 January 2020 to 24 March 2023, while the pixel offset data covered the period from 2016 to 2022. The choice of such an extensive data window provided comprehensive insight into the investigated events, including the possibility of identifying other unrecorded events and aiding in the development of more effective mitigation strategies. Furthermore, to supplement the analysis, a specific finite element method for slope stability analysis was used to reconstruct the deep geometry of the system, emphasizing the effect of groundwater-level flow on slope stability. All of the findings indicate that major landslide activities were concentrated during the heavy rainfall season, with movements ranging from several centimeters per year. These results were consistent with numerical analyses, which showed that the potential slip surface became significantly more unstable when the water table was elevated.