Nothing Special   »   [go: up one dir, main page]

You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Sign in to use this feature.

Years

Between: -

Subjects

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Journals

Article Types

Countries / Regions

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Search Results (284)

Search Parameters:
Keywords = rainfall induced landslide

Order results
Result details
Results per page
Select all
Export citation of selected articles as:
16 pages, 4171 KiB  
Article
Study on the Impact of Seepage Filtration Under Wet–Dry Cycles on the Stability of Mudstone Limestone Slopes
by Rui Li, Puyi Wang, Xiang Lu, Wei Zhou, Yihan Guo, Rongbo Lei, Zixiong Zhao, Ziyu Liu and Yu Tian
Water 2025, 17(4), 592; https://doi.org/10.3390/w17040592 - 18 Feb 2025
Abstract
Open-pit mining often exposes weak rock layers, the strength of which significantly affects the stability of slopes. If these rock layers are also prone to disintegration and expansion, cyclic rainfall can exacerbate instability. Rainfall-induced changes in the seepage field also indirectly threaten the [...] Read more.
Open-pit mining often exposes weak rock layers, the strength of which significantly affects the stability of slopes. If these rock layers are also prone to disintegration and expansion, cyclic rainfall can exacerbate instability. Rainfall-induced changes in the seepage field also indirectly threaten the stability of slopes. Therefore, investigating the characteristics of mudstone limestone and the impact of the seepage field on slope instability under different wet–dry cycles is of great significance for the safe mining of open-pit mines. This paper takes the mudstone limestone slope of a certain open-pit mine in the southwest as the starting point and conducts experiments on saturated density, water absorption rate, permeability coefficient, compressive strength, and variable angle shear strength. Combined with scanning electron microscopy and phase analysis of X-ray diffraction analysis, the macroscopic and microscopic characteristics of the samples are comprehensively analyzed. FLAC3D software is used to explore the changes in the seepage field and the mechanism of instability. Our research found that for the preparation of mudstone limestone samples, a particle size of less than 1 mm and a drying temperature of 50 °C are optimal, with specific values for initial natural and saturated density, and natural water content. As the number of wet–dry cycles increases, the saturated density of mudstone limestone increases; the water absorption rate first rises sharply and then rises slowly; the permeability coefficient first rises sharply and then stabilizes, finally dropping sharply; the compressive and shear strength decreases slowly, and the internal friction angle changes little; frequent cycles also lead to mudification and seepage filtration. At the microscopic level, pores become larger and more regular, and the distribution is more concentrated; changes in mineral content weaken the strength. Combined with numerical simulation, the changes in the seepage field at the bottom of the slope exceed those at the slope surface and top, the transient saturated area expands, and the overall and local slope stability coefficients gradually decrease. During the third cycle, the local stability is lower than the overall stability, and the landslide trend shifts. In conclusion, wet–dry cycles change the pores and mineral content, affecting the physical and mechanical properties, leading to the deterioration of the transient saturated area, a decrease in matrix suction, and an increase in surface gravity, eventually causing slope instability. Full article
Show Figures

Figure 1

Figure 1
<p>Diagram of the pilot study programme.</p>
Full article ">Figure 2
<p>Flow chart of the test.</p>
Full article ">Figure 3
<p>Flowchart of the wet and dry cycle scheme.</p>
Full article ">Figure 4
<p>Slope status and model establishment.</p>
Full article ">Figure 5
<p>Variation curve of specimen shear strength parameters.</p>
Full article ">Figure 6
<p>Effect of the number of wet and dry cycles on porosity and pore diameter.</p>
Full article ">Figure 7
<p>X-ray diffraction physical image analysis.</p>
Full article ">Figure 8
<p>Changes in seepage field in argillaceous limestone slopes after different rainfall cycles.</p>
Full article ">Figure 9
<p>Pore water pressure distribution of argillaceous limestone slope detection line under dry and wet cycle.</p>
Full article ">Figure 10
<p>Patterns of change between stability coefficients and depths of infiltration lines.</p>
Full article ">
24 pages, 15593 KiB  
Article
Study on Shallow Landslide Induced by Extreme Rainfall: A Case Study of Qichun County, Hubei, China
by Yousheng Li, Echuan Yan and Weibo Xiao
Water 2025, 17(4), 530; https://doi.org/10.3390/w17040530 - 12 Feb 2025
Viewed by 410
Abstract
In light of the increasing frequency of extreme rainfall events, there has been a concomitant rise in landslides triggered by such precipitation. Despite the extensive research conducted on rainfall-induced landslides, the practical implementation of these findings is constrained by geological and environmental factors. [...] Read more.
In light of the increasing frequency of extreme rainfall events, there has been a concomitant rise in landslides triggered by such precipitation. Despite the extensive research conducted on rainfall-induced landslides, the practical implementation of these findings is constrained by geological and environmental factors. Notably, there is a paucity of research on rainfall-induced shallow landslides in Hubei Province, China. Therefore, this study analyzes the fundamental characteristics and rainfall characteristics of landslides induced by multiple rounds of extreme rainfall in Qichun County in June and July 2016. The study explores the influence of five variables—namely, altitude, slope, slope aspect, stratum lithology, and rainfall—on landslides. The study uses numerical analysis to reveal the initiation mechanism of landslides. The research conclusions are as follows: The landslides within the study area are closely related to its natural topography, stratum lithology, and human activities. The majority of landslides are triggered by short-term extreme rainfall, while a smaller number are related to long-term continuous rainfall. The formation mechanism of landslides is primarily driven by dynamic water seepage, and the destruction process often lags behind the rainfall process. The conclusions can provide theoretical guidance for risk prevention and early warning of rainfall-induced landslides in the region. Full article
(This article belongs to the Special Issue Water-Related Landslide Hazard Process and Its Triggering Events)
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>) Hubei Province, China, (<b>b</b>) Qichun County, Huanggang City, (<b>c</b>) topographic map of Qichun County.</p>
Full article ">Figure 2
<p>Precipitation statistics in the study area from 2000 to 2019: (<b>a</b>) annual average rainfall, (<b>b</b>) monthly average rainfall.</p>
Full article ">Figure 3
<p>(<b>a</b>) The average monthly rainfall in 2016; (<b>b</b>) the average monthly rainfall in June and July from 2011 to 2019.</p>
Full article ">Figure 4
<p>The average daily and cumulative rainfall at 13 rainfall monitoring stations from June to July 2016.</p>
Full article ">Figure 5
<p>Comparison of daily rainfall at five different rainfall monitoring stations.</p>
Full article ">Figure 6
<p>Landslide volume statistics.</p>
Full article ">Figure 7
<p>Statistics of landslide area.</p>
Full article ">Figure 8
<p>Thickness statistics of sliding body.</p>
Full article ">Figure 9
<p>Disaster situation: (<b>a</b>) a house has been destroyed, (<b>b</b>) a house has been damaged, (<b>c</b>) a house has been destroyed, (<b>d</b>) an electric pole has been tilted, (<b>e</b>) a house has been damaged, (<b>f</b>) a house has been damaged.</p>
Full article ">Figure 10
<p>Number of landslides and altitude.</p>
Full article ">Figure 11
<p>Number of landslides and slope.</p>
Full article ">Figure 12
<p>Number of landslides and slope direction.</p>
Full article ">Figure 13
<p>Landslides and different stratum lithology.</p>
Full article ">Figure 14
<p>(<b>a</b>) Daily rainfall; (<b>b</b>) number of landslides per day.</p>
Full article ">Figure 15
<p>The relationship between the number of landslides and cumulative rainfall under different rainfall days.</p>
Full article ">Figure 16
<p>(<b>a</b>) Landslide deformation and damage, (<b>b</b>) 2D slope model, (<b>c</b>) rainfall conditions.</p>
Full article ">Figure 17
<p>The variations in pore water pressure: (<b>a</b>) initial, (<b>b</b>) 4th day, (<b>c</b>) 6th day, (<b>d</b>) 7th day, (<b>e</b>) 8th day, (<b>f</b>) 10th day.</p>
Full article ">Figure 18
<p>The variations in saturation: (<b>a</b>) initial, (<b>b</b>) 4th day, (<b>c</b>) 6th day, (<b>d</b>) 7th day, (<b>e</b>) 8th day, (<b>f</b>) 10th day.</p>
Full article ">Figure 19
<p>Changes in slope stability coefficient over time.</p>
Full article ">Figure 20
<p>Rainfall seepage damage to shallow landslide.</p>
Full article ">
28 pages, 5914 KiB  
Article
Predicting Landslide Deposit Zones: Insights from Advanced Sampling Strategies in the Ilopango Caldera, El Salvador
by Laura Paola Calderon-Cucunuba, Abel Alexei Argueta-Platero, Tomás Fernández, Claudio Mercurio, Chiara Martinello, Edoardo Rotigliano and Christian Conoscenti
Land 2025, 14(2), 269; https://doi.org/10.3390/land14020269 - 27 Jan 2025
Viewed by 483
Abstract
In landslide susceptibility modeling, research has predominantly focused on predicting landslides by identifying predisposing factors, often using inventories primarily based on the highest points of landslide crowns. However, a significant challenge arises when the transported mass impacts human activities directly, typically occurring in [...] Read more.
In landslide susceptibility modeling, research has predominantly focused on predicting landslides by identifying predisposing factors, often using inventories primarily based on the highest points of landslide crowns. However, a significant challenge arises when the transported mass impacts human activities directly, typically occurring in the deposition areas of these phenomena. Therefore, identifying the terrain characteristics that facilitate the transport and deposition of displaced material in affected areas is equally crucial. This study aimed to evaluate the predictive capability of identifying where displaced material might be deposited by using different inventories of specific parts of a landslide, including the source area, intermediate area, and deposition area. A sample segmentation was conducted that included inventories of these distinct parts of the landslide in the hydrographic basin of Lake Ilopango, which experienced debris flows and debris floods triggered by heavy rainfall from Hurricane Ida in November 2009. Given the extensive variables extracted for this evaluation (20 variables), the Induced Smoothed (IS) version of the Least Absolute Shrinkage and Selection Operator (LASSO) methodology was employed to determine the significance of each variable within the datasets. Additionally, the Multivariate Adaptive Regression Splines (MARS) algorithm was used for modeling. Our findings revealed that models developed using the deposition area dataset were more effective compared with those based on the source area dataset. Furthermore, the accuracy of models using deposition area data surpassed that of that using data from both the source and intermediate areas. Full article
Show Figures

Figure 1

Figure 1
<p>A northern basin of Ilopango Lake, El Salvador.</p>
Full article ">Figure 2
<p>Lithologic map of the study area. The map illustrates the geological composition of the study area, featuring pyroclastic rocks from Formation Tierra Blanca (s4) in the upslope, and Formation Comalapa (c1) in the middle slope. In the gentile and plain areas, alluvial deposits (Qf) and gravity deposits (Qd) are predominantly present.</p>
Full article ">Figure 3
<p>(<b>A1</b>,<b>A2</b>) Natural-color and false-color images from February 2009 (before the hurricane). (<b>B1</b>,<b>B2</b>) Natural-color and false-color images from November 2009 (after the hurricane). Source: Natural-color imagery obtained via Google Earth by Maxar Technologies; false-color imagery obtained via Google Earth Engine using data from the Aster sensor.</p>
Full article ">Figure 4
<p>Illustration of the four different approaches adopted for selecting cells within a landslide polygon: (<b>1</b>) a conceptual representation and (<b>2</b>) the actual subdivision in the dataset. The approaches were (<b>1a</b>) MAX: the cell with the highest elevation; (<b>1b</b>) SUP: randomly selected cells from the upper 10% of the landslide area; (<b>1c</b>) INF: randomly selected cells from the lower 10% of the landslide area; and (<b>1d</b>) BODY: randomly selected cells from across the entire landslide area.</p>
Full article ">Figure 5
<p>Landslides inventories. (<b>A</b>) “Debris Flow” dataset, which comprised landslides with similar geometries and sizes. (<b>B</b>) “Debris Flood” dataset, characterized by extensive deposit areas.</p>
Full article ">Figure 6
<p>Landslide mapping of debris floods (<b>A</b>) and debris flows (<b>B</b>).</p>
Full article ">Figure 7
<p>Significance of the variables observed on each dataset (MAX, SUP, INF, and BODY) for the debris flow database according to the IS-lasso method.</p>
Full article ">Figure 8
<p>Significance of the variables observed on each dataset (MAX, SUP, INF, and BODY) for the debris flood database according to the IS-lasso method.</p>
Full article ">Figure 9
<p>Percentage of each lithology category within landslide areas, calculated for the debris flow and debris flood databases. (1) Fm. Comalapa (acid pyroclastic). (2) Fm. Comalapa (acid effusive). (3) Fm. Tierras Blancas (pyroclastics). (4) Alluvial deposits. (5) Gravity deposits.</p>
Full article ">Figure 10
<p>Variable importance scores derived from the MARS models for each debris flow dataset (MAX, SUP, INF, and BODY).</p>
Full article ">Figure 11
<p>Variable importance scores derived from the MARS models for each debris flood dataset (MAX, SUP, INF, and BODY).</p>
Full article ">Figure 12
<p>Boxplots depicting the distributions of the AUC values for each model (BODY, INF, MAX, and SUP) applied to the debris flow testing datasets.</p>
Full article ">Figure 13
<p>Boxplots depicting the distributions of the AUC values for each model (BODY, INF, MAX, and SUP) applied to the debris flood testing datasets.</p>
Full article ">Figure 14
<p>Susceptibility map for debris flow deposits, obtained using the INF model from the debris flow database.</p>
Full article ">Figure 15
<p>Susceptibility map for debris flow source areas, obtained using the SUP model from the debris flow database.</p>
Full article ">Figure 16
<p>Susceptibility map for debris flow deposits, obtained using the INF model from the debris flood database.</p>
Full article ">Figure 17
<p>Boxplots illustrating the behavior of the variables SLO, ELE, TPI, and LSF in the MAX, SUP, INF, and BODY segments. Each boxplot displays the data distribution in the groups that represented areas prone to landslides compared with groups without landslides on debris flows.</p>
Full article ">Figure 18
<p>Boxplots illustrating the behavior of the variables CVI, RSP, TPI, and TWI in the MAX, SUP, INF, and BODY segments. Each boxplot displays data distribution in groups representing the areas prone to landslides compared with groups without landslides based on debris floods.</p>
Full article ">
22 pages, 18807 KiB  
Article
Development of a New Method for Debris Flow Runout Assessment in 0-Order Catchments: A Case Study of the Otoishi River Basin
by Ahmad Qasim Akbar, Yasuhiro Mitani, Ryunosuke Nakanishi, Hiroyuki Honda and Hisatoshi Taniguchi
Geosciences 2025, 15(2), 41; https://doi.org/10.3390/geosciences15020041 - 25 Jan 2025
Viewed by 656
Abstract
Debris flows are rapid, destructive landslides that pose significant risks in mountainous regions. This study presents a novel algorithm to simulate debris flow dynamics, focusing on sediment transport from 0-order basins to depositional zones. The algorithm integrates the D8 flow direction method with [...] Read more.
Debris flows are rapid, destructive landslides that pose significant risks in mountainous regions. This study presents a novel algorithm to simulate debris flow dynamics, focusing on sediment transport from 0-order basins to depositional zones. The algorithm integrates the D8 flow direction method with an adjustable friction coefficient to enhance the accuracy of debris flow trajectory and deposition modeling. Its performance was evaluated on three real-world cases in the Otoishi River basin, affected by rainfall-induced debris flows in July 2017, and the Aso Bridge landslide triggered by the 2016 Kumamoto Earthquake. By utilizing diverse friction coefficients, the study effectively captured variations in debris flow behavior, transitioning from fluid-like to more viscous states. Simulation results demonstrated a precision of 88.9% in predicting debris flow paths and deposition areas, emphasizing the pivotal role of the friction coefficient in regulating mass movement dynamics. Additionally, Monte Carlo (MC) simulations enhanced the identification of critical slip surfaces within 0-order basins, increasing the accuracy of debris flow source detection. This research offers valuable insights into debris flow hazards and risk mitigation strategies. The algorithm’s proven effectiveness in simulating real-world scenarios highlights its potential for integration into disaster risk assessment and prevention frameworks. By providing a reliable tool for hazard identification and prediction, this study supports proactive disaster management and aligns with the goals of sustainable development in regions prone to debris flow disasters. Full article
(This article belongs to the Special Issue Landslides Runout: Recent Perspectives and Advances)
Show Figures

Figure 1

Figure 1
<p>Schematic of debris flow and sediment transportation [<a href="#B5-geosciences-15-00041" class="html-bibr">5</a>].</p>
Full article ">Figure 2
<p>The location map of the study area is shown, overlaid with a Hillshade extracted from the 1 m resolution LiDAR DEM and the associated data used in this study. The red square area represents the study area located in Kyushu Island, Japan.</p>
Full article ">Figure 3
<p>Process map of pre-processor algorithm.</p>
Full article ">Figure 4
<p>Results of the preprocessor algorithm, illustrating the identified sources of debris flow as black polygons.</p>
Full article ">Figure 5
<p>D8 Flow Direction schematic.</p>
Full article ">Figure 6
<p>Schematic of slope and sediment particle distribution.</p>
Full article ">Figure 7
<p>Source of debris flow extracted from the results of MC calculations for the Otoishi River catchment [<a href="#B30-geosciences-15-00041" class="html-bibr">30</a>].</p>
Full article ">Figure 8
<p>Location map of case study areas where black polygons are the actual landslide after intense rain fall of July 2017.</p>
Full article ">Figure 9
<p>Illustration of every 20 steps of the debris flow runout simulation results for three case studies, as shown in <a href="#geosciences-15-00041-f008" class="html-fig">Figure 8</a>. <b>Case Study 1</b>: Result of debris flow movement from the source to its depositional zone with friction coefficients of 0.0001, 0.05, and 0.9. <b>Case Study 2:</b> Result of debris flow movement from the source to its depositional zone with friction coefficients of 0.0001, 0.05, and 0.9. <b>Case Study 3:</b> Result of debris flow movement from the source to its depositional zone with friction coefficients of 0.0001, 0.05, and 0.9.</p>
Full article ">Figure 9 Cont.
<p>Illustration of every 20 steps of the debris flow runout simulation results for three case studies, as shown in <a href="#geosciences-15-00041-f008" class="html-fig">Figure 8</a>. <b>Case Study 1</b>: Result of debris flow movement from the source to its depositional zone with friction coefficients of 0.0001, 0.05, and 0.9. <b>Case Study 2:</b> Result of debris flow movement from the source to its depositional zone with friction coefficients of 0.0001, 0.05, and 0.9. <b>Case Study 3:</b> Result of debris flow movement from the source to its depositional zone with friction coefficients of 0.0001, 0.05, and 0.9.</p>
Full article ">Figure 9 Cont.
<p>Illustration of every 20 steps of the debris flow runout simulation results for three case studies, as shown in <a href="#geosciences-15-00041-f008" class="html-fig">Figure 8</a>. <b>Case Study 1</b>: Result of debris flow movement from the source to its depositional zone with friction coefficients of 0.0001, 0.05, and 0.9. <b>Case Study 2:</b> Result of debris flow movement from the source to its depositional zone with friction coefficients of 0.0001, 0.05, and 0.9. <b>Case Study 3:</b> Result of debris flow movement from the source to its depositional zone with friction coefficients of 0.0001, 0.05, and 0.9.</p>
Full article ">Figure 9 Cont.
<p>Illustration of every 20 steps of the debris flow runout simulation results for three case studies, as shown in <a href="#geosciences-15-00041-f008" class="html-fig">Figure 8</a>. <b>Case Study 1</b>: Result of debris flow movement from the source to its depositional zone with friction coefficients of 0.0001, 0.05, and 0.9. <b>Case Study 2:</b> Result of debris flow movement from the source to its depositional zone with friction coefficients of 0.0001, 0.05, and 0.9. <b>Case Study 3:</b> Result of debris flow movement from the source to its depositional zone with friction coefficients of 0.0001, 0.05, and 0.9.</p>
Full article ">Figure 10
<p>Simulation results of the Aso Bridge landslide, Kumamoto. The yellow boundary, as referenced from [<a href="#B36-geosciences-15-00041" class="html-bibr">36</a>], indicates the actual landslide area.</p>
Full article ">Figure 11
<p>The debris movement from the source to the depositional zone: (<b>a</b>) a friction coefficient of 0.0001 applied with 2 m of mass moving from the source. (<b>b</b>) A friction coefficient of 0.5 applied with 2 m of mass moving from the source. (<b>c</b>) A friction coefficient of 0.9 applied with 2 m of mass moving from the source.</p>
Full article ">
20 pages, 22339 KiB  
Article
Evaluation of Rainfall-Induced Accumulation Landslide Susceptibility Based on Remote Sensing Interpretation
by Zhen Wu, Runqing Ye, Jue Huang, Xiaolin Fu and Yao Chen
Remote Sens. 2025, 17(2), 339; https://doi.org/10.3390/rs17020339 - 20 Jan 2025
Viewed by 558
Abstract
Landslide susceptibility evaluation is an indispensable part of disaster prevention and mitigation work. Selecting effective evaluation methods and models for landslide susceptibility assessment is of significant importance. This study focuses on selected areas in Yunyang County, Chongqing City. By interpreting high-resolution satellite remote [...] Read more.
Landslide susceptibility evaluation is an indispensable part of disaster prevention and mitigation work. Selecting effective evaluation methods and models for landslide susceptibility assessment is of significant importance. This study focuses on selected areas in Yunyang County, Chongqing City. By interpreting high-resolution satellite remote sensing images from before and after heavy rainfall on 31 August 2014, the distribution of rainfall-induced accumulation landslides was obtained. To evaluate the susceptibility of accumulation landslides, we have equated evaluation factors to accumulation distribution prediction factors. Eight evaluation factors were extracted using multi-source data, including lithology, elevation, slope, remote sensing image texture features, and the normalized difference vegetation index (NDVI). Various machine learning models, such as Random Forest (RF), Support Vector Machine (SVM), and BP Neural Network models, were employed to assess the susceptibility of rainfall-induced accumulation landslides in the study area. Subsequently, the accuracy of the evaluation models was compared and verified using the Receiver Operating Characteristic (ROC) curve, and the evaluation results were analyzed. Finally, the developed Random Forest model was applied to Gongping Town in Fengjie County to verify its applicability in other regions. The findings indicate that the complex geological conditions and the unique tectonic erosion landform patterns in the northeastern region of Chongqing not only make this area a center of heavy rainfall but also lead to frequent and recurrent rainfall-induced landslides. The Random Forest model effectively reflects the development characteristics of accumulation landslides in the study area. High and very high susceptibility zones are concentrated in the northern and central regions of the study area, while low and moderate susceptibility zones predominantly occupy the mountainous and riverside areas. Landslide susceptibility mapping in the study area shows that the Random Forest model yields reasonably graded results. Elevation, remote sensing image texture features, and lithology are highly significant factors in the evaluation system, indicating that the development factors of slope geological disasters in the study area are mainly related to topography, geomorphology, and lithology. The landslide susceptibility evaluation results in Gongping Town, Fengjie County, validate the applicability of the Random Forest model developed in this study to other regions. Full article
Show Figures

Figure 1

Figure 1
<p>Technical flowchart.</p>
Full article ">Figure 2
<p>Study area.</p>
Full article ">Figure 3
<p>Remote sensing interpretation distribution map of landslides induced by extreme rainfall on 31 August 2014.</p>
Full article ">Figure 4
<p>Remote sensing interpretation photos of accumulation landslides induced by rainfall around 31 August 2014.</p>
Full article ">Figure 5
<p>Correlation heatmap of evaluation factors.</p>
Full article ">Figure 6
<p>Evaluation factors for susceptibility to accumulation landslides.</p>
Full article ">Figure 6 Cont.
<p>Evaluation factors for susceptibility to accumulation landslides.</p>
Full article ">Figure 7
<p>Zoning maps of accumulation landslide susceptibility in the study area using different machine learning models.</p>
Full article ">Figure 8
<p>ROC curves and AUC values of three machine learning models.</p>
Full article ">Figure 9
<p>Comparison of landslide location distribution interpreted from remote sensing and susceptibility zoning results of accumulation landslides.</p>
Full article ">Figure 10
<p>Accumulation landslide susceptibility evaluation and remote sensing interpreted location map of Gongping Town, Fengjie County.</p>
Full article ">Figure 11
<p>Comparison of Guiba Landslide (<b>left</b>) interpreted image and Random Forest model accumulation landslide susceptibility evaluation results.</p>
Full article ">
28 pages, 24510 KiB  
Article
A Case Study of Using Numerical Analysis to Assess the Slope Stability of National Freeways in Northern Taiwan
by Hao-Wei Chiu, Yi-Hao Tsai, Chao-Wei Tang, Chih-Yu Chu and Shong-Loong Chen
Appl. Sci. 2025, 15(2), 635; https://doi.org/10.3390/app15020635 - 10 Jan 2025
Viewed by 487
Abstract
Taiwan is located at a junction of tectonic plates, which results in frequent earthquakes. Its terrain is mostly hilly, and its rainfall ranks among the highest in the world. Each of these elements affects the stability of slopes in various regions of Taiwan. [...] Read more.
Taiwan is located at a junction of tectonic plates, which results in frequent earthquakes. Its terrain is mostly hilly, and its rainfall ranks among the highest in the world. Each of these elements affects the stability of slopes in various regions of Taiwan. Several slopes along Taiwan’s Freeway 1 and 5 have experienced landslides and rockfalls. It is imperative that the slope stability of these national freeways be analyzed to avoid future slope collapses brought on by precipitation or other outside factors. Thus, three sites on Taiwan’s Freeway 1 and 5 were chosen for numerical slope stability analysis in this study. PLAXIS 2D CE (Version: 24.02.00.1144) finite element software was used in this study to simulate and analyze the safety of freeway slope protection projects. Displacements induced by normal and high groundwater levels were discussed. Moreover, a pseudo-static study of slope displacements under seismic conditions was performed. According to the results of the numerical study, the force operating on the slope was centered on the sliding surface when the groundwater level was normal, and it extended to the top when the groundwater level was high. By comparison, under seismic conditions, the force acting on the slope extended to the whole slope. Furthermore, the slope safety factor of Site 1 was greater than the design specification value in three different scenarios. This confirms that the slope protection project at Site 1 is effective. Full article
Show Figures

Figure 1

Figure 1
<p>Schematic diagram of Site 1.</p>
Full article ">Figure 2
<p>Schematic diagram of Site 2 and Site 3.</p>
Full article ">Figure 3
<p>Geological map of Site 1.</p>
Full article ">Figure 4
<p>Drilling location and core photos of Site 1.</p>
Full article ">Figure 5
<p>Slope cross-section and slope protection engineering facilities at Site 1.</p>
Full article ">Figure 6
<p>Geological map of Site 2.</p>
Full article ">Figure 7
<p>Slope cross-section and slope protection engineering facilities at Site 2.</p>
Full article ">Figure 8
<p>Geological map of Site 3.</p>
Full article ">Figure 9
<p>Slope cross-section and slope protection engineering facilities at Site 3.</p>
Full article ">Figure 10
<p>Setting of model boundary conditions.</p>
Full article ">Figure 11
<p>Schematic diagram of “create borehole”.</p>
Full article ">Figure 12
<p>Settings of soil parameters and groundwater level.</p>
Full article ">Figure 13
<p>PLAXIS 2D model of Site 1.</p>
Full article ">Figure 14
<p>Contour line diagram of total displacement with normal groundwater level at Site 1.</p>
Full article ">Figure 15
<p>Contour line diagram of total displacement with high groundwater level at Site 1.</p>
Full article ">Figure 16
<p>Contour line diagram of total displacement with normal groundwater level and pseudo-static analysis for Site 1.</p>
Full article ">Figure 17
<p>The potential sliding surface based on safety analysis in Scenario 1 at Site 1.</p>
Full article ">Figure 18
<p>The potential sliding surface based on safety analysis in Scenario 2 at Site 1.</p>
Full article ">Figure 19
<p>The potential sliding surface based on safety analysis in Scenario 3 at Site 1.</p>
Full article ">Figure 20
<p>The shading distribution diagram of the potential sliding surface based on safety analysis in Scenario 3 at Site 1.</p>
Full article ">Figure 21
<p>Site 1’s Σ<span class="html-italic">Msf</span> in various scenarios.</p>
Full article ">Figure 22
<p>PLAXIS 2D model for Site 2.</p>
Full article ">Figure 23
<p>Contour line diagram of total displacement with normal groundwater level for Site 2.</p>
Full article ">Figure 24
<p>Contour line diagram of total displacement with high groundwater level for Site 2.</p>
Full article ">Figure 25
<p>Contour line diagram of total displacement with normal groundwater level and pseudo-static analysis for Site 2.</p>
Full article ">Figure 26
<p>The potential sliding surface based on safety analysis in Scenario 1 at Site 2.</p>
Full article ">Figure 27
<p>The potential sliding surface based on safety analysis in Scenario 2 at Site 2.</p>
Full article ">Figure 28
<p>The potential sliding surface based on safety analysis in Scenario 3 at Site 2.</p>
Full article ">Figure 29
<p>The shading distribution diagram of the potential sliding surface based on safety analysis in Scenario 3 at Site 2.</p>
Full article ">Figure 30
<p>Site 2’s Σ<span class="html-italic">Msf</span> in various scenarios.</p>
Full article ">Figure 31
<p>PLAXIS 2D model for Site 3.</p>
Full article ">Figure 32
<p>Contour line diagram of total displacement with normal groundwater level for Site 3.</p>
Full article ">Figure 33
<p>Contour line diagram of total displacement with high groundwater level for Site 3.</p>
Full article ">Figure 34
<p>Contour line diagram of total displacement with normal groundwater level and pseudo-static analysis for Site 3.</p>
Full article ">Figure 35
<p>The potential sliding surface based on safety analysis in Scenario 1 at Site 3.</p>
Full article ">Figure 36
<p>The potential sliding surface based on safety analysis in Scenario 2 at Site 3.</p>
Full article ">Figure 37
<p>The potential sliding surface based on safety analysis in Scenario 3 at Site 3.</p>
Full article ">Figure 38
<p>The shading distribution diagram of the potential sliding surface based on safety analysis in Scenario 3 at Site 3.</p>
Full article ">Figure 39
<p>Site 3’s Σ<span class="html-italic">Msf</span> in various scenarios.</p>
Full article ">
18 pages, 3009 KiB  
Article
Ecological Sensitivity of the Mata Allo Sub-Watershed, South Sulawesi: A Spatial Analysis Using Principal Component Analysis
by Syamsu Rijal, Samsuri, Heni Masruroh, Munajat Nursaputra, Chairil A and Nur Zamzam Putri Ardi
Sustainability 2025, 17(2), 447; https://doi.org/10.3390/su17020447 - 9 Jan 2025
Viewed by 639
Abstract
Watersheds are critical ecosystems that provide essential services, but they face increasing threats from deforestation, land use changes, and climate variability. The Mata Allo Sub-Watershed, which is characterized by steep topography and high rainfall, is particularly vulnerable to erosion, landslides, and habitat loss, [...] Read more.
Watersheds are critical ecosystems that provide essential services, but they face increasing threats from deforestation, land use changes, and climate variability. The Mata Allo Sub-Watershed, which is characterized by steep topography and high rainfall, is particularly vulnerable to erosion, landslides, and habitat loss, necessitating robust conservation strategies. This study used principal component analysis (PCA) to assess ecological sensitivity, focusing on slope, rainfall, vegetation density, and land cover. The PCA results identified land cover as the most influential positive factor in F1 (loading value: 0.588), increasing sensitivity due to human-induced land use changes, while rainfall contributed most negatively (−0.638) by potentially mitigating extreme ecological risks. These contrasting roles underscore the complexity of interactions shaping watershed sensitivity. Slope strongly influenced F2 (−0.795), explaining 26.48% of the variance and highlighting the critical role of steep slopes in exacerbating erosion risks. Vegetation density in F3 (−0.679) and rainfall in F4 (−0.724) played significant roles in stabilizing soil and mitigating ecological risks, emphasizing their importance in reducing watershed sensitivity. The “Extremely Sensitive” class covers 48.79% of the watershed, primarily in areas with steep slopes and sparse vegetation, while “High Sensitivity” areas occupy 34.93%. Projections for 2032 suggest a reduction in “Extremely Sensitive” zones to 41.00%, reflecting improvements from targeted management interventions. These findings provide a foundation for promoting sustainable watershed management, enhancing climate resilience, and supporting biodiversity conservation efforts in vulnerable regions. Full article
(This article belongs to the Section Sustainable Forestry)
Show Figures

Figure 1

Figure 1
<p>Map of study area.</p>
Full article ">Figure 2
<p>Land cover map of Mata Allo Sub-Watershed in 2014 (<b>a</b>), 2018 (<b>b</b>), and 2023 (<b>c</b>).</p>
Full article ">Figure 3
<p>Projected land cover map of Mata Allo Sub-Watershed in 2032.</p>
Full article ">Figure 4
<p>Comparison of monthly average rainfall and percentage change in rainfall in the Mata Allo Sub-Watershed: actual (2014–2023) vs. projection (2032).</p>
Full article ">Figure 5
<p>Map of average annual rainfall comparison in the Mata Allo Sub-Watershed: actual (2014–2023) (<b>a</b>) vs. projection (2032) (<b>b</b>).</p>
Full article ">Figure 6
<p>Maps of actual (<b>a</b>) and predicted (<b>b</b>) sensitivity index of Mata Allo Watershed.</p>
Full article ">
22 pages, 4716 KiB  
Article
Global Sensitivity Analysis of Slope Stability Considering Effective Rainfall with Analytical Solutions
by Chuan-An Xia, Jing-Quan Zhang, Hao Wang and Wen-Bin Jian
Water 2025, 17(2), 141; https://doi.org/10.3390/w17020141 - 7 Jan 2025
Viewed by 559
Abstract
Rainfall-induced landslides are widely distributed in many countries. Rainfall impacts the hydraulic dynamics of groundwater and, therefore, slope stability. We derive an analytical solution of slope stability considering effective rainfall based on the Richards equation. We define effective rainfall as the total volume [...] Read more.
Rainfall-induced landslides are widely distributed in many countries. Rainfall impacts the hydraulic dynamics of groundwater and, therefore, slope stability. We derive an analytical solution of slope stability considering effective rainfall based on the Richards equation. We define effective rainfall as the total volume of rainfall stored within a given range of the unsaturated zone during rainfall events. The slope stability at the depth of interest is provided as a function of effective rainfall. The validity of analytical solutions of system states related to effective rainfall, for infinite slopes of a granite residual soil, is verified by comparing them with the corresponding numerical solutions. Additionally, three approaches to global sensitivity analysis are used to compute the sensitivity of the slope stability to a variety of factors of interest. These factors are the reciprocal of the air-entry value of the soil α, the thickness of the unsaturated zone L, the cohesion of soil c, the internal friction angle ϕ related to the effective normal stress, the slope angle β, the unit weights of soil particles γs, and the saturated hydraulic conductivity Ks. The results show the following: (1) The analytical solutions are accurate in terms of the relative differences between the analytical and the numerical solutions, which are within 5.00% when considering the latter as references. (2) The temporal evolutions of the shear strength of soil can be sequentially characterized as four periods: (i) strength improvement due to the increasing weight of soil caused by rainfall infiltration, (ii) strength reduction controlled by the increasing pore water pressure, (iii) strength reduction due to the effect of hydrostatic pressure in the transient saturation zone, and (iv) stable strength when all the soil is saturated. (3) The large α corresponds to high effective rainfall. (4) The factors ranked in descending order of sensitivity are as follows: α > L > c > β > γs > Ks > ϕ. Full article
Show Figures

Figure 1

Figure 1
<p>Diagram showing infinite slope with slope angle <span class="html-italic">β</span> under the effect of rainfall, where <span class="html-italic">z</span> = 0 and <span class="html-italic">L</span> indicate the locations of the water table and the ground surface, respectively.</p>
Full article ">Figure 2
<p>Temporal evolution of (<b>a</b>) <span class="html-italic">S<sub>e</sub></span>, (<b>b</b>) <span class="html-italic">S<sub>r</sub></span>, (<b>c</b>) <span class="html-italic">φ</span>, and (<b>d</b>) <span class="html-italic">θ</span> obtained through the analytical and the numerical methods.</p>
Full article ">Figure 3
<p>Temporal evolution of <span class="html-italic">q<sub>w</sub></span> obtained through the analytical and numerical methods.</p>
Full article ">Figure 4
<p>Temporal evolution of the values <math display="inline"><semantics> <mrow> <msubsup> <mi>R</mi> <mi>e</mi> <mo>*</mo> </msubsup> </mrow> </semantics></math> and <span class="html-italic">Q<sub>w</sub></span> (<b>a</b>) and the corresponding values of <math display="inline"><semantics> <mrow> <msubsup> <mi>I</mi> <mi>p</mi> <mo>*</mo> </msubsup> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mfenced> <mrow> <mn>1</mn> <mo>−</mo> <msubsup> <mi>I</mi> <mi>p</mi> <mo>*</mo> </msubsup> </mrow> </mfenced> </mrow> </semantics></math> (<b>b</b>) when <span class="html-italic">L</span> = 2, 3, and 8.</p>
Full article ">Figure 5
<p>Temporal evolution in the values of <math display="inline"><semantics> <mrow> <msubsup> <mi>R</mi> <mi>e</mi> <mo>*</mo> </msubsup> </mrow> </semantics></math> and <span class="html-italic">Q<sub>w</sub></span> (<b>a</b>) and the corresponding results of <math display="inline"><semantics> <mrow> <msubsup> <mi>I</mi> <mi>p</mi> <mo>*</mo> </msubsup> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mfenced> <mrow> <mn>1</mn> <mo>−</mo> <msubsup> <mi>I</mi> <mi>p</mi> <mo>*</mo> </msubsup> </mrow> </mfenced> </mrow> </semantics></math> (<b>b</b>) when <span class="html-italic">α</span> = 0.32, 1.28, and 5.12.</p>
Full article ">Figure 6
<p>Temporal evolution of the values of <math display="inline"><semantics> <mrow> <msubsup> <mi>R</mi> <mi>e</mi> <mo>*</mo> </msubsup> </mrow> </semantics></math> and <span class="html-italic">Q<sub>w</sub></span> (<b>a</b>) and the corresponding results of <math display="inline"><semantics> <mrow> <msubsup> <mi>I</mi> <mi>p</mi> <mo>*</mo> </msubsup> </mrow> </semantics></math> and <math display="inline"><semantics> <mrow> <mfenced> <mrow> <mn>1</mn> <mo>−</mo> <msubsup> <mi>I</mi> <mi>p</mi> <mo>*</mo> </msubsup> </mrow> </mfenced> </mrow> </semantics></math> (<b>b</b>) when with <math display="inline"><semantics> <mrow> <msub> <mi>q</mi> <mrow> <mi>B</mi> <mo>*</mo> </mrow> </msub> </mrow> </semantics></math> = 0.1, 0.5, and 1.0.</p>
Full article ">Figure 7
<p>Temporal evolution of the values of <math display="inline"><semantics> <mi>τ</mi> </semantics></math> (left column) and <math display="inline"><semantics> <mrow> <mo>Δ</mo> <mi>τ</mi> <mo>%</mo> </mrow> </semantics></math> (right column) when <span class="html-italic">α</span> = 0.2 (<b>a</b>,<b>b</b>), 1.28 (<b>c</b>,<b>d</b>), and 5.12 (<b>e</b>,<b>f</b>) with (solid lines) and without (cross symbols) effective rainfall.</p>
Full article ">Figure 8
<p>Temporal evolution of the values of <math display="inline"><semantics> <mi>τ</mi> </semantics></math> (left column) and <math display="inline"><semantics> <mrow> <mo>Δ</mo> <mi>τ</mi> <mo>%</mo> </mrow> </semantics></math> (right column) when <math display="inline"><semantics> <mrow> <msub> <mi>q</mi> <mrow> <mi>B</mi> <mo>*</mo> </mrow> </msub> </mrow> </semantics></math> = 0.1 (<b>a</b>,<b>b</b>), 0.5 (<b>c</b>,<b>d</b>), and 1.0 (<b>e</b>,<b>f</b>) with (solid lines) and without (cross symbols) effective rainfall.</p>
Full article ">Figure 9
<p><span class="html-italic">R<sub>Ti</sub></span>, <math display="inline"><semantics> <mrow> <mover accent="true"> <mrow> <mover accent="true"> <mi>δ</mi> <mo stretchy="true">^</mo> </mover> </mrow> <mo stretchy="true">^</mo> </mover> </mrow> </semantics></math><span class="html-italic"><sub>i</sub></span>, and <span class="html-italic">S<sub>i</sub></span> (with <span class="html-italic">i</span> = <span class="html-italic">c</span>, <span class="html-italic">φ</span>, and <span class="html-italic">L</span>) versus <span class="html-italic">N</span> when <span class="html-italic">t</span> = 0.</p>
Full article ">Figure 10
<p>Temporal evolution of <span class="html-italic">F<sub>s</sub></span> corresponding to the first 100 realizations (gray lines) of factors when considering effective rainfall (<b>a</b>). The corresponding temporal changes in the relative differences of <span class="html-italic">F<sub>s</sub></span> (gray lines) with and without effective rainfall (<b>b</b>). The red lines indicate the means of the values presented as gray lines.</p>
Full article ">Figure 11
<p><span class="html-italic">R<sub>Ti</sub></span>, <math display="inline"><semantics> <mrow> <mover accent="true"> <mrow> <mover accent="true"> <mi>δ</mi> <mo stretchy="true">^</mo> </mover> </mrow> <mo stretchy="true">^</mo> </mover> </mrow> </semantics></math><span class="html-italic"><sub>i</sub></span>, and <span class="html-italic">S<sub>i</sub></span> (with <span class="html-italic">i</span> = <span class="html-italic">c</span>, <span class="html-italic">φ</span>, <span class="html-italic">β</span>, <span class="html-italic">γ</span>, <span class="html-italic">K<sub>s</sub></span>, <span class="html-italic">α</span>, and <span class="html-italic">L</span>) when <span class="html-italic">t</span> = 24 h and <span class="html-italic">N</span> = 5000, with (solid lines) and without (dashed lines) effective rainfall.</p>
Full article ">Figure 12
<p>Temporal evolution of the arithmetic mean of the three indexes (MI) for each factor when <span class="html-italic">N</span> = 5000.</p>
Full article ">
20 pages, 3362 KiB  
Article
Stability Prediction Model of Transmission Tower Slope Based on ISCSO-SVM
by Zilong Zhang, Xiaoliang Liu, Yanhai Wang, Enyang Li and Yuhao Zhang
Electronics 2025, 14(1), 126; https://doi.org/10.3390/electronics14010126 - 31 Dec 2024
Viewed by 455
Abstract
Landslides induced by heavy rainfall are common in southern China and pose significant risks to the safe operation of transmission lines. To ensure the reliability of transmission line operations, this paper presents a stability prediction model for transmission tower slopes based on the [...] Read more.
Landslides induced by heavy rainfall are common in southern China and pose significant risks to the safe operation of transmission lines. To ensure the reliability of transmission line operations, this paper presents a stability prediction model for transmission tower slopes based on the Improved Sand Cat Swarm Optimization (ISCSO) algorithm and Support Vector Machine (SVM). The ISCSO algorithm is enhanced with dynamic reverse learning and triangular wandering strategies, which are then used to optimize the kernel and penalty parameters of the SVM, resulting in the ISCSO-SVM prediction model. In this study, a typical transmission tower slope in southern China is used as a case study, with the transmission tower slope database generated through orthogonal experimental design and Geo-studio simulations. In addition to traditional input features, an additional input—transmission tower catchment area—is incorporated, and the stable state of the transmission tower slope is set as the predicted output. The results demonstrate that the ISCSO-SVM model achieves the highest prediction accuracy, with the smallest errors across all metrics. Specifically, compared to the standard SVM, the MAPE, MAE, and RMSE values are reduced by 70.96%, 71.41%, and 57.37%, respectively. The ISCSO-SVM model effectively predicts the stability of transmission tower slopes, thereby ensuring the safe operation of transmission lines. Full article
Show Figures

Figure 1

Figure 1
<p>Test function for F1–F6.</p>
Full article ">Figure 2
<p>ISCSO-SVM predictive model flow chart.</p>
Full article ">Figure 3
<p>(<b>a</b>) Results of buckling condition. (<b>b</b>) Results of understable state. (<b>c</b>) Results of basically stable condition. (<b>d</b>) Results of steady condition.</p>
Full article ">Figure 4
<p>(<b>a</b>) Prediction results of SSA-SVM. (<b>b</b>) Prediction results of SVM. (<b>c</b>) Prediction results of SCSO-SVM. (<b>d</b>) Prediction results of ISCSO-SVM.</p>
Full article ">
20 pages, 9663 KiB  
Article
Research on the Failure Mechanism and Treatment Technology of Landslides in Typical Accumulation Bodies Along Highways in Qinghai Province
by Yunfei Yang, Zixuan Yang, Wanzhong Xu, Fayou A, Yinghang Guo and Jieru Zheng
Water 2025, 17(1), 34; https://doi.org/10.3390/w17010034 - 26 Dec 2024
Viewed by 507
Abstract
Landslides on the Jiaxi Highway in Qinghai Province threaten construction safety and quality. The on-site data analysis shows that excavation at the foot of the slope and heavy rainfall are the key factors causing the displacement of the Q1 monitoring point by 1825 [...] Read more.
Landslides on the Jiaxi Highway in Qinghai Province threaten construction safety and quality. The on-site data analysis shows that excavation at the foot of the slope and heavy rainfall are the key factors causing the displacement of the Q1 monitoring point by 1825 mm. This article uses numerical simulation methods combined with the strength reduction method to study the stability changes of slopes under different working conditions. Numerical simulations identified the landslide location and predicted a 1960 mm slip and a safety factor of 1.26 under natural conditions, indicating risks. The study adopted a strategy combining slope cutting, load reduction, and sheet pile wall reinforcement. After the first treatment, the safety factor rose to 1.83 with a 40 mm displacement; after the second, it reached 2.36 with a 37 mm displacement. Continuous monitoring showed a 50 mm displacement over six months, indicating stability. Rainfall simulations before and after treatment explained the stability evolution and local slope stability. Treatments increased the safety factor to 2.16 with a 17.6 mm displacement. This study significantly improved highway landslide stability and verified treatment effectiveness, providing a reference for similar geological conditions. Full article
Show Figures

Figure 1

Figure 1
<p>Geographic location of the study area.</p>
Full article ">Figure 2
<p>Borehole information for the study area. (<b>a</b>) Upper Triassic Nanying’er Formation sandstone; (<b>b</b>) Quaternary Holocene floodplain sand pebbles.</p>
Full article ">Figure 3
<p>Slope displacement monitoring points.</p>
Full article ">Figure 4
<p>Slope displacement monitoring data.</p>
Full article ">Figure 5
<p>Panoramic image of the landslide.</p>
Full article ">Figure 6
<p>H1 landslide boundary.</p>
Full article ">Figure 7
<p>Geological profile of section 1-1′.</p>
Full article ">Figure 8
<p>Numerical simulation model diagram.</p>
Full article ">Figure 9
<p>Displacement and plastic zone variations under natural working conditions: (<b>a</b>) x-displacement; (<b>b</b>) plastic zone.</p>
Full article ">Figure 10
<p>First cut slope + support.</p>
Full article ">Figure 11
<p>Effectiveness of governance: (<b>a</b>) x-displacement; (<b>b</b>) plastic zone.</p>
Full article ">Figure 12
<p>Second slope cutting + support.</p>
Full article ">Figure 13
<p>Effectiveness of governance: (<b>a</b>) x-displacement; (<b>b</b>) plastic zone.</p>
Full article ">Figure 14
<p>Post-disposal displacement monitoring data.</p>
Full article ">Figure 15
<p>Pore water pressure.</p>
Full article ">Figure 16
<p>Variation of safety factor with seepage time.</p>
Full article ">Figure 17
<p>Slope displacement monitoring: (<b>a</b>) monitoring point 1; (<b>b</b>) monitoring point 2.</p>
Full article ">Figure 18
<p>Variation of pore water pressure in slope: (<b>a</b>) monitoring point 1; (<b>b</b>) monitoring point 2.</p>
Full article ">Figure 19
<p>Pore water pressure (A1, A2 and A3 are slope displacement monitoring points).</p>
Full article ">Figure 20
<p>Slope displacement monitoring.</p>
Full article ">
15 pages, 10783 KiB  
Article
Evaluation of the Effects of Rainfall Infiltration Boundaries on the Stability of Unsaturated Soil Slopes Using the Particle Flow Code
by Jian Zhang, Fangrui Hu, Qi Zhang, Jun Wang, Wenting Deng, Li Zhang and Xiaoquan Shao
Water 2024, 16(24), 3704; https://doi.org/10.3390/w16243704 - 22 Dec 2024
Viewed by 623
Abstract
Rainfall infiltration is the primary triggering factor for the instability of unsaturated slopes. At present, rainfall-induced landslides are mainly considered to be influenced by the overall infiltration conditions, while few investigations have been conducted on the influence of infiltration boundaries on slope instability. [...] Read more.
Rainfall infiltration is the primary triggering factor for the instability of unsaturated slopes. At present, rainfall-induced landslides are mainly considered to be influenced by the overall infiltration conditions, while few investigations have been conducted on the influence of infiltration boundaries on slope instability. This study proposes a rainfall infiltration method using a discrete element model (DEM), which is based on saturated–unsaturated seepage theory. The influence of three infiltration boundaries on the instability of homogeneous unsaturated soil slopes was studied. The results showed that the infiltration rate of a rainfall-covered slope crest was faster than that of rainfall-covered slope surfaces. A transient saturated zone was formed on the slope surface after a certain duration of rainfall. Rain infiltration boundary conditions significantly impact the saturation distribution, seepage field, failure mode, and failure period. The safety and stability factors for the rainfall-covered slope crest and full area decreased monotonically with the increase in rainfall duration, while there was a brief increase at the initial stage of rainfall before a quick decline for rainfall-covered slope surfaces. This research provides a preliminary exploration of the impact of rainfall boundary conditions on the instability of slopes, offering a reference basis for DEM simulations that consider slope stability under the influence of rainfall infiltration. Full article
Show Figures

Figure 1

Figure 1
<p>Flowchart of the analysis of seepage behavior with the DEM.</p>
Full article ">Figure 2
<p>The geometry and boundary conditions of a simple slope.</p>
Full article ">Figure 3
<p>The balls in the DEM models and monitoring points.</p>
Full article ">Figure 4
<p>The stress–strain curves of soil under different confining pressures.</p>
Full article ">Figure 5
<p>The variation in the degree of saturation during the rainfall process at different monitoring points: (<b>a</b>) rainfall-covered slope crest; (<b>b</b>) rainfall-covered slope surface; (<b>c</b>) full rainfall-covered slope area.</p>
Full article ">Figure 6
<p>Vector field of the seepage velocity of the slope in different rainfall infiltration scenarios: (<b>a</b>) rainfall-covered slope crest; (<b>b</b>) rainfall-covered slope surface; (<b>c</b>) full rainfall-covered slope area.</p>
Full article ">Figure 7
<p>Crack distribution, landslide debris, and slope displacement during slope failure: (<b>a</b>) rainfall-covered slope crest; (<b>b</b>) rainfall-covered slope surface; (<b>c</b>) full rainfall-covered slope area.</p>
Full article ">Figure 8
<p>Failure pattern obtained from the model test of a homogeneous sandy slope in [<a href="#B29-water-16-03704" class="html-bibr">29</a>]: (<b>a</b>) rainfall-covered slope crest; (<b>b</b>) rainfall-covered slope surface; (<b>c</b>) full rainfall-covered slope area.</p>
Full article ">Figure 9
<p>Evolution law of crack numbers in different rainfall infiltration scenarios.</p>
Full article ">Figure 10
<p>Evolution law of the slope safety factor during rainfall.</p>
Full article ">
16 pages, 4207 KiB  
Article
Calibration and Performance Evaluation of Cost-Effective Capacitive Moisture Sensor in Slope Model Experiments
by Muhammad Nurjati Hidayat, Hemanta Hazarika and Haruichi Kanaya
Sensors 2024, 24(24), 8156; https://doi.org/10.3390/s24248156 - 20 Dec 2024
Viewed by 688
Abstract
Understanding the factors that contribute to slope failures, such as soil saturation, is essential for mitigating rainfall-induced landslides. Cost-effective capacitive soil moisture sensors have the potential to be widely implemented across multiple sites for landslide early warning systems. However, these sensors need to [...] Read more.
Understanding the factors that contribute to slope failures, such as soil saturation, is essential for mitigating rainfall-induced landslides. Cost-effective capacitive soil moisture sensors have the potential to be widely implemented across multiple sites for landslide early warning systems. However, these sensors need to be calibrated for specific applications to ensure high accuracy in readings. In this study, a soil-specific calibration was performed in a laboratory setting to integrate the soil moisture sensor with an automatic monitoring system using the Internet of Things (IoT). This research aims to evaluate a low-cost soil moisture sensor (SKU:SEN0193) and develop calibration equations for the purpose of slope model experiment under artificial rainfall condition using silica sand. The results indicate that a polynomial function is the best fit, with a coefficient of determination (R2) ranging from 0.918 to 0.983 and a root mean square error (RMSE) ranging from 1.171 to 2.488. The calibration equation was validated through slope model experiments, with soil samples taken from the models after the experiment finished. Overall, the moisture content readings from the sensors showed approximately a 12% deviation from the actual moisture content. The findings suggest that the cost-effective capacitive soil moisture sensor has the potential to be used for the development of landslide early warning system. Full article
(This article belongs to the Section Electronic Sensors)
Show Figures

Figure 1

Figure 1
<p>Equipment used: (<b>a</b>) moisture sensor, (<b>b</b>) soil moisture meter, (<b>c</b>) M5Stack Core2, and (<b>d</b>) ENV III sensor.</p>
Full article ">Figure 2
<p>Overview of calibration experiment setup.</p>
Full article ">Figure 3
<p>Schematic of the soil box for sensor calibration: (<b>a</b>) 3D view of soil box and six soil moisture sensors, and (<b>b</b>) soil samples were collected from the blue dashed line area near the moisture sensors to measure actual water content using a soil sampler.</p>
Full article ">Figure 4
<p>Particle size distribution of K7 sand.</p>
Full article ">Figure 5
<p>Overview of slope model and moisture sensor layout.</p>
Full article ">Figure 6
<p>Calibration of moisture sensor using a linear equation.</p>
Full article ">Figure 7
<p>Calibration of moisture sensor using logarithmic equation.</p>
Full article ">Figure 8
<p>Calibration of moisture sensor using polynomial equation.</p>
Full article ">Figure 9
<p>Room temperature, humidity and air pressure during the calibration process.</p>
Full article ">Figure 10
<p>Moisture content under soil temperature variation.</p>
Full article ">Figure 11
<p>Comparison of sensor moisture content with actual moisture content on different rainfall intensities.</p>
Full article ">Figure 12
<p>Temperature, humidity and air pressure during slope model experiment (100 mm/h).</p>
Full article ">Figure 13
<p>Comparison of derived moisture content with actual moisture content during the calibration process.</p>
Full article ">Figure 14
<p>Calibration equation of Moisture Sensor number 2.</p>
Full article ">
22 pages, 22575 KiB  
Article
Back Analysis of Rainfall-Induced Landslide in Cimanggung District of Sumedang Regency in West Java Using Deterministic and Probabilistic Analyses
by Dwi Sarah, Zulfahmi Zulfahmi, Moch Hilmi Zaenal Putra, Nendaryono Madiutomo, Gunawan Gunawan, Sumaryadi Sumaryadi and Deden Agus Ahmid
Geosciences 2024, 14(12), 347; https://doi.org/10.3390/geosciences14120347 - 17 Dec 2024
Viewed by 844
Abstract
Rainfall-induced landslides are widespread in Indonesia, particularly in West Java, where volcanic residual soils are typically stable but may become unstable during heavy rainfall. This study aims to back analyze the geotechnical factors contributing to the Cimanggung landslide in 2021. The methods applied [...] Read more.
Rainfall-induced landslides are widespread in Indonesia, particularly in West Java, where volcanic residual soils are typically stable but may become unstable during heavy rainfall. This study aims to back analyze the geotechnical factors contributing to the Cimanggung landslide in 2021. The methods applied in this study include site investigations, laboratory testing, and numerical modeling. We performed deterministic, coupled seepage-slope stability analysis and Monte Carlo probabilistic analysis to assess the slope performance prior to and after rainfall infiltration. The results reveal that the initial water level significantly affects slope stability, and heavy rainfall infiltration triggered the landslide’s initiation. The deep water table (over 20 m below ground level) maintains the slope stability, and increasing the water table to 16 m compromises its stability. Heavy rainfall infiltration reduces suction in the unsaturated zone, decreasing the shear strength and triggering landslides. The heavy rainfall infiltration did not penetrate deep enough to raise the water table; rather, poor urban drainage on the upper slope caused it. Rainfall infiltration caused wetting in the upper zone, weakening the slope and causing loss of support. It is recommended that effective drainage management and integrated slope monitoring be applied to mitigate landslide risks in this region. Full article
(This article belongs to the Section Natural Hazards)
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>) Topographic map of the Cimanggung study area in Sumedang Regency, West Java derived from DEMNAS data, showing major fault systems and volcanic features. (<b>b</b>) Geological map of the Bandung-Garut Quadrangle, highlighting lithological units, including young volcanic products (Qyu), basaltic lava (Qyl), and lake deposits (Ql). The study area is marked with a red rectangle.</p>
Full article ">Figure 2
<p>Maximum annual daily rainfall in Cimanggung District from 2003 to 2017.</p>
Full article ">Figure 3
<p>(<b>a</b>) The 3D satellite view of the landslide area at Mount Aseupan’s foothills, marked by a red polygon, with urbanized flat terrain and vegetated hills. (<b>b</b>) The 3D geological map showing the dominance of undifferentiated young volcanic products (Qyu). (<b>c</b>,<b>d</b>) Illustrations of the conditions before (December 2020) and after (May 2021) the landslide, when the Pondok Daud Housing Complex was entirely destroyed and five houses in the SBG Housing Complex were severely damaged, with some carried downslope. (<b>e</b>,<b>f</b>) Slope profiles from boreholes S-01 and S-02, showing subsurface layers of silty sand, silty clay, and clay.</p>
Full article ">Figure 4
<p>(<b>a</b>) Soil water characteristic curve and (<b>b</b>) hydraulic conductivity function.</p>
Full article ">Figure 5
<p>Slope model with boundary conditions.</p>
Full article ">Figure 6
<p>(<b>a</b>) Factor of safety and probability of failure vs. initial water table. (<b>b</b>) Factor of safety and reliability index vs. initial water table.</p>
Full article ">Figure 7
<p>Pore water pressure profiles observed at specific times (t3, t4, t7, t8, and t24) for different rainfall intensities and durations. In each subplot, t3 represents the pore water pressure profile at 3 h, t4 represents that at 4 h, t7 represents that at 7 h, t8 represents that at 8 h, and t24 represents that at 24 h. The corresponding rainfall conditions are (<b>a</b>) 100 mm/day for 3 h, (<b>b</b>) 142 mm/day for 3 h, (<b>c</b>) 192 mm/day for 3 h, (<b>d</b>) 100 mm/day for 4 h, (<b>e</b>) 142 mm/day for 4 h, (<b>f</b>) 192 mm/day for 4 h, (<b>g</b>) 100 mm/day for 5 h, (<b>h</b>) 142 mm/day for 5 h, (<b>i</b>) 192 mm/day for 5 h, (<b>j</b>) 100 mm/day for 24 h, (<b>k</b>) 142 mm/day for 24 h, and (<b>l</b>) 192 mm/day for 24 h.</p>
Full article ">Figure 7 Cont.
<p>Pore water pressure profiles observed at specific times (t3, t4, t7, t8, and t24) for different rainfall intensities and durations. In each subplot, t3 represents the pore water pressure profile at 3 h, t4 represents that at 4 h, t7 represents that at 7 h, t8 represents that at 8 h, and t24 represents that at 24 h. The corresponding rainfall conditions are (<b>a</b>) 100 mm/day for 3 h, (<b>b</b>) 142 mm/day for 3 h, (<b>c</b>) 192 mm/day for 3 h, (<b>d</b>) 100 mm/day for 4 h, (<b>e</b>) 142 mm/day for 4 h, (<b>f</b>) 192 mm/day for 4 h, (<b>g</b>) 100 mm/day for 5 h, (<b>h</b>) 142 mm/day for 5 h, (<b>i</b>) 192 mm/day for 5 h, (<b>j</b>) 100 mm/day for 24 h, (<b>k</b>) 142 mm/day for 24 h, and (<b>l</b>) 192 mm/day for 24 h.</p>
Full article ">Figure 8
<p>(<b>a</b>) Factor of safety over 3 h of rainfall. (<b>b</b>) Factor of safety over 4 h of rainfall. (<b>c</b>) Factor of safety over 5 h of rainfall. (<b>d</b>) Factor of safety over 24 h of rainfall.</p>
Full article ">
23 pages, 22589 KiB  
Article
Landslide Prediction Validation in Western North Carolina After Hurricane Helene
by Sophia Lin, Shenen Chen, Ryan A. Rasanen, Qifan Zhao, Vidya Chavan, Wenwu Tang, Navanit Shanmugam, Craig Allan, Nicole Braxtan and John Diemer
Geotechnics 2024, 4(4), 1259-1281; https://doi.org/10.3390/geotechnics4040064 - 14 Dec 2024
Viewed by 961
Abstract
Hurricane Helene triggered 1792 landslides across western North Carolina and has caused damage to 79 bridges to date. Helene hit western North Carolina days after a low-pressure system dropped up to 254 mm of rain in some locations of western North Carolina (e.g., [...] Read more.
Hurricane Helene triggered 1792 landslides across western North Carolina and has caused damage to 79 bridges to date. Helene hit western North Carolina days after a low-pressure system dropped up to 254 mm of rain in some locations of western North Carolina (e.g., Asheville Regional Airport). The already waterlogged region experienced devastation as significant additional rainfall occurred during Helene, where some areas, like Asheville, North Carolina received an additional 356 mm of rain (National Weather Service, 2024). In this study, machine learning (ML)-generated multi-hazard landslide susceptibility maps are compared to the documented landslides from Helene. The landslide models use the North Carolina landslide database, soil survey, rainfall, USGS digital elevation model (DEM), and distance to rivers to create the landslide variables. From the DEM, aspect factors and slope are computed. Because recent research in western North Carolina suggests fault movement is destabilizing slopes, distance to fault was also incorporated as a predictor variable. Finally, soil types were used as a wildfire predictor variable. In total, 4794 landslides were used for model training. Random Forest and logistic regression machine learning algorithms were used to develop the landslide susceptibility map. Furthermore, landslide susceptibility was also examined with and without consideration of wildfires. Ultimately, this study indicates heavy rainfall and debris-laden floodwaters were critical in triggering both landslides and scour, posing a dual threat to bridge stability. Field investigations from Hurricane Helene revealed that bridge damage was concentrated at bridge abutments, with scour and sediment deposition exacerbating structural vulnerability. We evaluated the assumed flooding potential (AFP) of damaged bridges in the study area, finding that bridges with lower AFP values were particularly vulnerable to scour and submersion during flood events. Differentiating between landslide-induced and scour-induced damage is essential for accurately assessing risks to infrastructure. The findings emphasize the importance of comprehensive hazard mapping to guide infrastructure resilience planning in mountainous regions. Full article
Show Figures

Figure 1

Figure 1
<p>Study area with location map illustrating North Carolina’s mountain area. (<b>a</b>) North Carolina’s distinct physiographic region distribution, (<b>b</b>) Blue Ridge Mountain area, and (<b>c</b>) hypothetical Appalachian Mountain formation during the Alleghenian orogeny.</p>
Full article ">Figure 2
<p>Path of Hurricane Helene moving through the Gulf of Mexico and landing near Perry, Florida as a Category 4 storm. Note the mountainous topography of western North Carolina.</p>
Full article ">Figure 3
<p>A composite representation of damaged bridges and landslide locations after Hurricane Helene.</p>
Full article ">Figure 4
<p>Some of the landslide locations after Hurricane Helene. (Photo credit: Shen-En Chen, Sophia Lin, and Qifan Zhao).</p>
Full article ">Figure 5
<p>A schematic of the calculation workflow for the probability of multi-hazard (wildfire, landslide, earthquake, and flooding) occurrence map, the probability of wildfire occurrence map, and of bridges of average flooding potential (AFP). Note that L+W+E represents landslides, wildfires, and earthquakes, and L+E represents landslides and earthquakes.</p>
Full article ">Figure 6
<p>Multi-hazard (without wildfire effect) risk map of North Carolina.</p>
Full article ">Figure 7
<p>Multi-hazard (with wildfire effect) susceptibility map of North Carolina.</p>
Full article ">Figure 8
<p>Multi-hazard susceptibility map in North Carolina with reported landslide locations: (<b>a</b>) landslide, wildfire, and earthquake; (<b>b</b>) landslide and earthquake.</p>
Full article ">Figure 9
<p>Analysis of reported landslides with the corresponding susceptibility probabilities: (<b>a</b>) multi-hazard scenario L+W+E; (<b>b</b>) multi-hazard scenario L+E; (<b>c</b>) difference between L+W+E and L+E; and (<b>d</b>) bar chart comparing the two scenarios by number of slides. Note that L+W+E represents landslides, wildfires, and earthquakes and L+E represents landslides and earthquakes.</p>
Full article ">Figure 10
<p>Hurricane Helene landslide damage to transportation structures and facilities: (<b>a</b>) by a roadside near Lake Lure; (<b>b</b>) by a parking space near Chimney Rock; (<b>c</b>) near a parking lot in Chimney Rock Village; and (<b>d</b>) below a county highway in Henderson County (Photo credit: Shen-En Chen, Sophia Lin, and Qifan Zhao).</p>
Full article ">Figure 11
<p>Hurricane Helene landslide damage to bridge structures: (<b>a</b>) Main Street bridge over a railroad, Saluda, NC; (<b>b</b>) bridge near Lake Lure; (<b>c</b>) the Big Hungry Road Bridge, Flat Rock; and (<b>d</b>) dam crossing, Lake Lure. (Photo credit: Shen-En Chen, Sophia Lin, and Qifan Zhao).</p>
Full article ">Figure 12
<p>Hurricane Helene flood-battered region in Chimney Rock Village, NC: (<b>a</b>) washed away bridge on the Chimney Rock Scenic Road over the Broad River, Chimney Rock Village, NC; (<b>b</b>) view from Main Street looking over Broad River; (<b>c</b>) scoured Broad River valley in front of Burnshirt Vineyards Bistro on Main Street, Chimney Rock Village, NC; and (<b>d</b>) the parking lot in front of Burnshirt Vineyards Bistro on Main Street, Chimney Rock Village, NC. (Photo credit: Shen-En Chen, Sophia Lin, and Qifan Zhao).</p>
Full article ">Figure 13
<p>Helene landslides and the associated susceptibility values as an accumulated function. Susceptibility values for the following multi-hazard scenarios: L+E (landslide and earthquake) and L+W+E (landslide, wildfire, and earthquake).</p>
Full article ">Figure 14
<p>Landslides with zero and 99~100% predictions for (<b>a</b>) without wildfire effects and (<b>b</b>) with wildfire effects.</p>
Full article ">Figure 15
<p>Conditioning factors used in this study, including reported landslides: (<b>a</b>) Elevation; (<b>b</b>) slope; (<b>c</b>) aspect; (<b>d</b>) soil type; (<b>e</b>) rainfall; (<b>f</b>) temperature; (<b>g</b>) forest cover; (<b>h</b>) distance to rivers; (<b>i</b>) distance to faults; (<b>j</b>) distance to roads; (<b>k</b>) distance to high population density; and (<b>l</b>) probability of wildfire occurrence.</p>
Full article ">Figure 16
<p>Typical bridge scour damage mechanism, including the formation of scour holes (local scour) around bridge piers, which can result in increased stress in the supporting geo-medium (riverbed material): (<b>a</b>) typical scour mechanism; (<b>b</b>) geo-medium stressing due to scour hole formation; (<b>c</b>) scour depths due to clear water scour vs. live-bed scour.</p>
Full article ">Figure 17
<p>Debris slide and scour combined mass waste mechanism of the Big Hungry River: (<b>a</b>) whole view of the Big Hungry Road (County route 1889) landslide, Flat Rock, NC; and (<b>b</b>) closeup of the slide and the river deposits, and (<b>c</b>) landslide assumption by [<a href="#B36-geotechnics-04-00064" class="html-bibr">36</a>].</p>
Full article ">Figure 18
<p>Reconstruction of the Big Hungry Road Bridge: (<b>a</b>) on the Flat Rock side; (<b>b</b>) on the Flat Rock side; (<b>c</b>) on the Flat Rock side, and (<b>d</b>) opposite to Flat Rock.</p>
Full article ">
19 pages, 8336 KiB  
Article
Analysis of the Differences Between Two Landslides on One Slope in Yongguang Village Based on Physical Models and Groundwater Identification
by Fucun Lu, Kun Liu, Shunhua Xu, Jianyu Zhang and Dingnan Guo
Water 2024, 16(24), 3591; https://doi.org/10.3390/w16243591 - 13 Dec 2024
Viewed by 635
Abstract
In 2013, a Ms 6.6 earthquake occurred at the boundary of Min County and Zhang County, triggering numerous landslides. Notably, two landslides with significantly different sliding characteristics emerged less than 100 m apart in Yongguang Village, Min County. The eastern landslide was characterized [...] Read more.
In 2013, a Ms 6.6 earthquake occurred at the boundary of Min County and Zhang County, triggering numerous landslides. Notably, two landslides with significantly different sliding characteristics emerged less than 100 m apart in Yongguang Village, Min County. The eastern landslide was characterized by instability induced by seismic inertial forces, whereas the western landslide exhibited flow slides triggered by liquefaction in loess. To further analyze the causes of these landslides, this study employed a 1 m depth ground temperature survey to probe the shallow groundwater in the area, aiming to understand the distribution of shallow groundwater. Based on the results from the 1 m depth ground temperature survey, a random forest model was applied to regressively predict the initial groundwater levels. The TRIGRS model was utilized to evaluate the influence of pre-earthquake rainfall conditions on landslide stability, and the pore water pressure outputs from TRIGRS were integrated with the Scoops3D model to analyze landslide stability under seismic effects. The results indicate that the combination of the 1 m depth ground temperature survey with high-density electrical methods and random forest approaches effectively captures the initial groundwater levels across the region. Notably, the heavy rainfall occurring one day prior to the earthquake did not significantly reduce the stability of the landslide in Yongguang Village. Instead, the abundant groundwater in the source area of the western landslide, combined with several months of pre-earthquake rainfall, resulted in elevated groundwater levels that created favorable conditions for its occurrence. While the primary triggering factor for both landslides in Yongguang Village was the earthquake, the distinct topographic and groundwater conditions led to significantly different sliding characteristics under seismic influence at the same slope. Full article
Show Figures

Figure 1

Figure 1
<p>Topographic and Google Earth imagery of Yongguang Village.</p>
Full article ">Figure 2
<p>Flowchart showing the methodology of this study.</p>
Full article ">Figure 3
<p>Schematic diagram of 1 m depth ground temperature survey.</p>
Full article ">Figure 4
<p>Schematic diagram of cell grid (cited from Li et al. [<a href="#B34-water-16-03591" class="html-bibr">34</a>]).</p>
Full article ">Figure 5
<p>Results of 1 m depth ground temperature survey: (<b>a</b>) distribution of temperature measurement points; (<b>b</b>) distribution of 1 m ground temperature and inferred direction of groundwater flow vein in Yongguang Village.</p>
Full article ">Figure 6
<p>Results of the high-density electrical method and distribution of measurement lines.</p>
Full article ">Figure 7
<p>Initial groundwater level in the study area.</p>
Full article ">Figure 8
<p>Prediction accuracy of the random forest (RF) model.</p>
Full article ">Figure 9
<p>Illustration of the auxiliary variable used in co-kriging interpolation of soil layer thickness and the resulting interpolation outcomes: (<b>a</b>) schematic diagram illustrating the relative position P<sub>1</sub>, P<sub>2</sub> at a point on the slope surface; (<b>b</b>) soil thickness in the study area.</p>
Full article ">Figure 10
<p>Results of the TRIGRS model: (<b>a</b>) dry conditions, (<b>b</b>) actual rainfall conditions, and (<b>c</b>) heavy rainfall conditions.</p>
Full article ">Figure 11
<p>Results of the Scoops3D model: (<b>a</b>) conditions without groundwater; (<b>b</b>) conditions with actual rainfall and groundwater.</p>
Full article ">Figure 12
<p>Groundwater storage capacity.</p>
Full article ">Figure 13
<p>Historical landslide distribution in the study area.</p>
Full article ">
Back to TopTop