Modeling of rainfall-triggered shallow landslide

By integrating hydrological modeling with the infinite slope stability analysis, a rainfall-triggered shallow landslide model was developed by Iverson (Water Resour Res 36:1897-1910, 2000). In Iverson’s model, the infiltration capacity is assumed to be equivalent to the saturated hydraulic conductivity for finding pressure heads analytically. However, for general infiltration process, the infiltration capacity should vary with time during the period of rain, and the infiltration rate is significantly related to the variable infiltration capacity. To avoid the unrealistically high pressure heads, Iverson employed the beta-line correction by specifying that the simulated pressure heads cannot exceed those given by the beta line. In this study, the suitability of constant infiltration capacity together with the beta-line correction for hydrological modeling and landslide modeling of hillslope subjected to a rainfall is examined. By amending the boundary condition at ground surface of hillslope in Iverson’s model, the modified Iverson’s model with considering general infiltration process is developed to conduct this examination. The results show that the unrealistically high pressure heads from Iverson’s model occur due to the overestimation of infiltration rate induced from the assumption that the infiltration capacity is identical to the saturated hydraulic conductivity. Considering with the general infiltration process, the modified Iverson’s model gives acceptable results. In addition, even though the beta-line correction is applied, the Iverson’s model still produces greater simulated pressure heads and overestimates soil failure potential as compared with the modified Iverson’s model. Therefore, for assessing rainfall-triggered shallow landslide, the use of constant infiltration capacity together with the beta-line correction needs to be replaced by the consideration of general infiltration process.     

Spatio-temporal estimation of shallow landslide hazard triggered by rainfall using a three-dimensional model

A shallow landslide triggered by rainfall can be forecast in real-time by modeling the relationship between rainfall infiltration and decrease of slope stability. This paper describes a promising approach that combines an improved three-dimensional slope stability model with an approximate method based on the Green and Ampt model, to estimate the time–space distribution of shallow landslide hazards. Once a forecast of rainfall intensity and slope stability-related data, e.g., terrain and geology data, are acquired, this approach is shown to have the ability to estimate the variation of slope stability of a wide natural area during rainfall and to identify the location of potential failure surfaces. The effectiveness of the estimation procedures described has been tested by comparison with a one-dimensional method and by application to a landslide-prone area in Japan.     

Evaluation of potential landslide damming: Case study of Urni landslide,Kinnaur, Satluj valley,India

This work aims to understand the process of potential landslide damming using slope failure mechanism,dam dimension and dam stability evaluation. The Urni landslide, situated on the right bank of the Satluj River, Himachal Pradesh(India) is taken as the case study. The Urni landslide has evolved into a complex landslide in the last two decade(2000-2016) and has dammed the Satluj River partially since year 2013,damaging ～200 m stretch of the National Highway(NH-05). The crown of the landslide exists at an altitude of ～2180-2190 m above msl, close to the Urni village that has a human population of about 500.The high resolution imagery shows ～50 m long landslide scarp and ～100 m long transverse cracks in the detached mass that implies potential for further slope failure movement. Further analysis shows that the landslide has attained an areal increase of 103,900 ± 1142 m^2 during year 2004-2016. About 86% of this areal increase occurred since year 2013. Abrupt increase in the annual mean rainfall is also observed since the year 2013. The extreme rainfall in the June, 2013; 11 June(～100 mm) and 16 June(～115 mm),are considered to be responsible for the slope failure in the Urni landslide that has partially dammed the river. The finite element modelling(FEM) based slope stability analysis revealed the shear strain in the order of 0.0-0.16 with 0.0-0.6 m total displacement in the detachment zone. Further, kinematic analysis indicated planar and wedge failure condition in the jointed rockmass. The debris flow runout simulation of the detached mass in the landslide showed a velocity of ～25 m/s with a flow height of ～15 m while it(debris flow) reaches the valley floor. Finally, it is also estimated that further slope failure may detach as much as 0.80 ±0.32 million m^3 mass that will completely dam the river to a height of 76±30 m above the river bed.     

Centrifuge modelling of landslides and landslide hazard mitigation: A review

Landslides are serious geohazards that occur under a variety of climatic conditions and can cause many casualties and significant economic losses. Centrifuge modelling, as a representative type of physical modelling, provides a realistic simulation of the stress level in a small-scale model and has been applied over the last 50 years to develop a better understanding of landslides. With recent developments in this technology, the application of centrifuge modelling in landslide science has significantly increased. Here, we present an overview of physical models that can capture landslide processes during centrifuge modelling. This review focuses on (i) the experimental principles and considerations, (ii) landslide models subjected to various triggering factors, including centrifugal acceleration, rainfall, earthquakes, water level changes, thawing permafrost, excavation, external loading and miscellaneous conditions, and (iii) different methods for mitigating landslides modelled in centrifuge, such as the application of nails, piles, geotextiles, vegetation, etc. The behaviors of all the centrifuge models are discussed, with emphasis on the deformation and failure mechanisms and experimental techniques. Based on this review, we provide a best-practice methodology for preparing a centrifuge landslide test and propose further efforts in terms of the seven aspects of model materials, testing design and equipment, measurement methods, scaling laws, full-scale test applications, landslide early warning, and 3D modelling to better understand the complex behaviour of landslides.     

Analysis of time-varying rainfall infiltration induced landslide

A case study of rainstorm-induced landslide is modeled following the hourly rainfall time series from the stage of infiltration caused by induced slope movement and soil saturated to excess pore pressures—Transient Rainfall Infiltration and Grid-based Regional Slope-Stability Analysis (TRIGRS). The grid-based landslide stability analysis was conducted to model the increased pore pressures and runoff in the study area under the specified rainfall conditions. The generated time variances of pore pressures help determine landslide characteristics and mechanisms under rainfall conditions. Inputs of soil properties and permeability parameters for landslide stability analysis in the study area were prepared by TRIGRS adopted for transient infiltration analysis. Results of the analyses show that under heavy rainfall conditions, the infiltrated slope is unstable and the time of debris masses movement initiated is correlated to the recorded time. In the initiated landslide, characteristics and effects are considered and reflected in the numerical modeling under combinations of topography, land use, climatic and geological conditions. Results reveal that there is a plane failure surface and a potential circular failure surface at the study site besides the rock topple failures in the crest. A grid-based slope-stability analysis incorporated with the GIS spatial functions is more advantageous than the traditional two-dimensional analysis for specified slope profiles to determine the whole behavior of a slope.     

Analytic network process model for landslide hazard zonation

Various controlling factors such as lithology, slope angle, slope aspect, landuse, channel proximity etc. are generally considered for the landslide hazard assessment. Although outer dependence of these parameters to a landslide is inevitably taken into account, inter-dependence among the factors is seldom addressed. Analytic Network Process (ANP) is the multi-criteria decision making (MCDM) tool which takes into account such a complex relationship among parameters. In this research, an ANP model for landslide susceptibility is proposed, priority weights for each parameter controlling the landslide were determined, and a hazard map was prepared of an area in a fragile mountainous terrain in the eastern part of Nepal. The data used in the example were derived from published sources, aerial photographs and a topographic map. However, the procedures developed can readily incorporate additional information from more detailed investigations.     

Permanent seismic deformation analysis of a landslide

A failed slope may not necessarily require a remedial treatment if it can be shown with confidence that the maximum movement of the slide mass will be within tolerable limits, i.e., not cause loss of life or property. A permanent displacement analysis of a landslide for static and seismic conditions is presented using a continuum mechanics approach. Computed values of displacement for static conditions compare favorably with field measurements and computed values of seismic displacements for a postulated earthquake motion appear reasonable. Also, the seismic displacements using the continuum mechanics approach compare favorably with those obtained using the Newmark sliding block procedure for assessing seismically-induced slope deformations.     

Implementing low-cost landslide risk reduction: a pilot study in unplanned housing areas of the Caribbean

Landslides pose a serious physical and environmental threat to vulnerable communities living in areas of unplanned housing on steep slopes in the Caribbean. Some of these communities have, in the past, had to be relocated, at costs of millions of dollars, because of major slides triggered by tropical storm rainfall. Even so, evidence shows that: (1) risk reduction is a marginal activity; (2) there has been minimal uptake of hazard maps and vulnerability assessments and (3) there is little on-the-ground delivery of construction for risk reduction. This article directly addresses these issues by developing a low-cost approach to the identification of the potential pore pressure changes that trigger such slides we seek to address these three commentaries directly. A complex 45–60° slope site in St Lucia, West Indies was selected as a pilot for a modelling approach that uses numerical models (FLAC and CHASM) to verify the need for surface water management to effectively reduce landslide risk. Following the model confirmation, a series of drains were designed and constructed at the site. Post-construction evidence indicates the methodology to be sound, in that the site was stable in subsequent 1-in-1 to 1-in-4 year rainfall events. A critical feature of the approach is that it is community-based from data acquisition through to community members participating in construction.     

Uncertainty pattern in landslide susceptibility prediction modelling: Effects of different landslide boundaries and spatial shape expressions

In some studies on landslide susceptibility mapping (LSM), landslide boundary and spatial shape characteristics have been expressed in the form of points or circles in the landslide inventory instead of the accurate polygon form. Different expressions of landslide boundaries and spatial shapes may lead to substantial differences in the distribution of predicted landslide susceptibility indexes (LSIs); moreover, the presence of irregular landslide boundaries and spatial shapes introduces uncertainties into the LSM. To address this issue by accurately drawing polygonal boundaries based on LSM, the uncertainty patterns of LSM modelling under two different landslide boundaries and spatial shapes, such as landslide points and circles, are compared. Within the research area of Ruijin City in China, a total of 370 landslides with accurate boundary information are obtained, and 10 environmental factors, such as slope and lithology, are selected. Then, correlation analyses between the landslide boundary shapes and selected environmental factors are performed via the frequency ratio (FR) method. Next, a support vector machine (SVM) and random forest (RF) based on landslide points, circles and accurate landslide polygons are constructed as point-, circle- and polygon-based SVM and RF models, respectively, to address LSM. Finally, the prediction capabilities of the above models are compared by computing their statistical accuracy using receiver operating characteristic analysis, and the uncertainties of the predicted LSIs under the above models are discussed. The results show that using polygonal surfaces with a higher reliability and accuracy to express the landslide boundary and spatial shape can provide a markedly improved LSM accuracy, compared to those based on the points and circles. Moreover, a higher degree of uncertainty of LSM modelling is present in the expression of points because there are too few grid units acting as model input variables. Additionally, the expression of the landslide boundary as circles introduces errors in measurement and is not as accurate as the polygonal boundary in most LSM modelling cases. In addition, the results under different conditions show that the polygon-based models have a higher LSM accuracy, with lower mean values and larger standard deviations compared with the point- and circle-based models. Finally, the overall LSM accuracy of the RF is superior to that of the SVM, and similar patterns of landslide boundary and spatial shape affecting the LSM modelling are reflected in the SVM and RF models.     

Micro-sheeting of granite and its relationship with landsliding specifically after the heavy rainstorm in June 1999, Hiroshima Prefecture, Japan

Certain types of granite in mountainous areas are microscopically sheeted to a depth of 50 m due to unloading under the stress field that reflects slope morphology. Micro-sheets generally strike parallel to major slope surfaces and gently dip downslope, forming cataclinal overdip slopes. The cataclinal overdip slope accelerates creep movement of micro-sheeted granite, which in turn loosens and disintegrates granite via the widening or neoformation of cracks, probably in combination with stress release, temperature change, and changes in water content near the ground surface. The surface portion of micro-sheeted granite is thus loosened with a well-defined basal front, which finally slides in response to heavy rain. Innumerable landslides of this type occurred in Hiroshima Prefecture, western Japan, following the heavy rainstorm of 29 June 1999. Following such landslides, the weathering of micro-sheeted granite exposed on the landslide scar recommences, setting the stage for future landslide.     

暴雨型滑坡降水入渗机理分析

暴雨能诱发大量的滑坡发生,在不同地区存在不同的临界降雨值,每当雨强大于该值时,滑坡就会大面积发生,这已经被多起滑坡事件证实。一般地,入渗理论认为,当雨强大于一定值时,入渗率是常量,过强的降水会转化为地表径流,这与雨强越大滑坡越易发生、群发性越强的观察结果不一致。在斜坡体上存在大量的裂隙,它们对降水入渗有着很大的贡献,应用有限单元法分别模拟了有裂隙和无裂隙时斜坡体内瞬态渗流场的变化,并用Bishop法对不同埋深滑动面的稳定性做了计算分析。模拟结果表明：考虑裂隙时的计算结果与宏观观察结果更一致,在评价降雨入渗对斜坡稳定性的影响时,一定要充分考虑裂隙的存在,只有雨强大于一个临界值时才会有充足的降水通过裂隙渗入到斜坡深部,裂隙的影响才显著地表现出来。     

Early warning system for shallow landslides using rainfall threshold and slope stability analysis

A combined cluster and regression analysis were performed for the first time to identify rainfall threshold that triggers landslide events in Amboori, Kerala, India. Amboori is a tropical area that is highly vulnerable to landslides. The 2, 3, and 5-day antecedent rainfall data versus daily rainfall was clustered to identify a cluster of critical events that could potentially trigger landslides. Further, the cluster of critical events was utilized for regression analysis to develop the threshold equations. The 5-day antecedent (x-variable) vs. daily rainfall (y-variable) provided the best fit to the data with a threshold equation of y = 80.7–0.1981x. The intercept of the equation indicates that if the 5-day antecedent rainfall is zero, the minimum daily rainfall needed to trigger the landslide in the Amboori region would be 80.7 mm. The negative coefficient of the antecedent rainfall indicates that when the cumulative antecedent rainfall increases, the amount of daily rainfall required to trigger monsoon landslide decreases. The coefficient value indicates that the contribution of the 5-day antecedent rainfall is ∼20% to the landslide trigger threshold. The slope stability analysis carried out for the area, using Probabilistic Infinite Slope Analysis Model (PISA-m), was utilized to identify the areas vulnerable to landslide in the region. The locations in the area where past landslides have occurred demonstrate lower Factors of Safety (FS) in the slope stability analysis. Thus, rainfall threshold analysis together with the FS values from slope stability can be suitable for developing a simple, cost-effective, and comprehensive early-warning system for shallow landslides in Amboori and similar regions.     

Downslope volume enlargement of a debris slide–debris flow in the 1999 Hiroshima, Japan, rainstorm

Following a heavy rainstorm on 29 June 1999, hundreds of slope failures occurred in granitic mountains in Hiroshima Prefecture, Japan. Among these events, a highly mobile landslide (termed the Kameyama landslide in this paper), which occurred in Kameyama area of Hiroshima city, was the most catastrophic, and was investigated in detail. The displaced soil mass from the source area of this landslide traveled about 300 m and deposited a volume more than 10 times as great as that in the source area. The landslide originated in and traversed a valley-shaped concave slope covered by pre-existing colluvial debris deposits. In addition, a spring was visible in the source area and very shallow ground water was observed in an observation pit dug in the source area. Thus, it is inferred that the ground-water table rose quickly during the rainfall, and that this rise triggered the slope failure in the source area. Based on a field survey along the landslide cross section, a possible explanation for the mechanism of the landslide was obtained: the displaced soil mass from the source area impacted the debris deposit in the path of the landslide, thus triggering liquefaction failure of the saturated part of debris. The original landslide and the liquefied debris then moved downslope as a single mass. To examine this assumption, ring-shear tests were performed on samples taken from the source area. Undrained ring-shear tests on the saturated samples showed that the sample is highly liquefiable, and liquefaction failure could have been triggered in the debris deposits by a very small impact from the displaced soil mass of the initial failure. In addition, laboratory tests simulating the impacts on the debris deposits at natural water content, i.e., unsaturated (at the survey time, 2 days after the failure) showed that although shear failure could be caused by the assumed impact force, the displaced soils stopped after a few centimeters displacement, indicating that existence of a saturated zone in debris deposits is prerequisite for this kind of failure.     

A simple definition of a landslide

This study makes a significant progress in addressing the challenges of short-term slope displacement prediction in the Universal Landslide Monitoring Program, an unprecedented disaster mitigation program in China, where lots of newly established monitoring slopes lack sufficient historical deformation data, making it difficult to extract deformation patterns and provide effective predictions which plays a crucial role in the early warning and forecasting of landslide hazards. A slope displacement prediction method based on transfer learning is therefore proposed. Initially, the method transfers the deformation patterns learned from slopes with relatively rich deformation data by a pre-trained model based on a multi-slope integrated dataset to newly established monitoring slopes with limited or even no useful data, thus enabling rapid and efficient predictions for these slopes. Subsequently, as time goes on and monitoring data accumulates, fine-tuning of the pre-trained model for individual slopes can further improve prediction accuracy, enabling continuous optimization of prediction results. A case study indicates that, after being trained on a multi-slope integrated dataset, the TCN-Transformer model can efficiently serve as a pre-trained model for displacement prediction at newly established monitoring slopes. The three-day average RMSE is significantly reduced by 34.6% compared to models trained only on individual slope data, and it also successfully predicts the majority of deformation peaks. The fine-tuned model based on accumulated data on the target newly established monitoring slope further reduced the three-day RMSE by 37.2%, demonstrating a considerable predictive accuracy. In conclusion, taking advantage of transfer learning, the proposed slope displacement prediction method effectively utilizes the available data, which enables the rapid deployment and continual refinement of displacement predictions on newly established monitoring slopes.     

台风暴雨条件下滑坡稳定性影响因素分析

以实验数据为依据,通过数值模拟对台风暴雨条件下滑坡体稳定性的主要影响因素进行模拟研究。台风暴雨条件下滑坡体的主要影响因素为特殊降雨雨型：(1) 在相同雨强程度下,台风暴雨的单峰型降雨水文响应速度、稳定性破坏速度及程度均大于普通均匀型降雨,单峰型降雨导致边坡短时间内达到失稳状态,呈“即雨即滑”情况;(2) 强台风多峰型降雨下,经过第一个峰值强降雨后浅层土质边坡已处于不稳定状态,当第二个强峰值降雨来临时,坡体的稳定性较单峰型降低更为迅速。这些为福建省台风暴雨这一特殊降雨类型条件下滑坡的发生机理探索及预警预报提供了新的思路和科学依据。     

Using an airborne laser scanner for the identification of shallow landslides and susceptibility assessment in an area of ignimbrite overlain by permeable pyroclastics

An airborne laser scanner can identify shallow landslides even when they are only several meters in diameter and are hidden by vegetation, if the vegetation is coniferous or deciduous trees in a season with fewer leaves. We used an airborne laser scanner to survey an area of the 1998 Fukushima disaster, during which more than 1,000 shallow landslides occurred on slopes of vapor-phase crystallized ignimbrite overlain by permeable pyroclastics. We identified landslides that have occurred at the 1998 event and also previous landslides that were hidden by vegetation. The landslide density of slopes steeper than 20° was 117 landslides/km² before the 1998 disaster. This event increased the density by 233 landslides/km² indicating that this area is highly susceptible to shallow landsliding.     

A time-space based approach for mapping rainfall-induced shallow landslide hazard

A rainfall-induced shallow landslide is a major hazard in mountainous terrain, but a time-space based approach is still an unsettled issue for mapping rainfall-induced shallow landslide hazards. Rain induces a rise of the groundwater level and an increase in pore water pressure that results in slope failures. In this study, an integrated infinite slope analysis model has been developed to evaluate the influence of infiltration on surficial stability of slopes by the limit equilibrium method. Based on this new integrated infinite slope analysis model, a time-space based approach has been implemented to map the distributed landslide hazard in a GIS (Geographic Information Systems) and to evaluate the shallow slope failure induced by a particular rainfall event that accounts for the rainfall intensity and duration. The case study results in a comprehensive time-space landslide hazard map that illustrates the change of the safety factor and the depth of the wetting front over time.     

Influence of rainfall-induced wetting on the stability of slopes in weathered soils

Shallow failures of slopes in weathered soil are caused by infiltration due to prolonged rainfall. These failures are mainly triggered by the deepening of the wetting band accompanied by a decrease in matric suction induced by the water infiltration. This paper reports trends of rainfall-induced wetting band depth in two types of weathered soils that are commonly found in Korea. Both theoretical and numerical analyses for wetting band depth are presented based on the soil–water characteristic curve obtained using filter paper as well as tensiometer tests. It is found that the magnitude of wetting front suction plays a key role in the stability of slopes in weathered soils. Theoretical analysis based on modified Green and Ampt model tends to underestimate the wetting band depth for typical Korean weathered soils. It was also deduced that for Korean weathered soils, the factor of safety drops rapidly once the wetting band depth of 1.2 m reached.     

SHALSTAB模型在浙南林溪流域滑坡预测中的应用

受台风暴雨影响,浙南林溪流域滑坡频发。针对该区域滑坡规模小、长度与厚度比值大的特点,采用浅层滑坡稳定性模型(SHALSTAB)对潜在滑坡进行了预测,以log(降雨量q/土壤的导水系数T)作为划分标准,结果显示随着log(q/T)值的提高,预测的滑坡区域逐渐扩大,预测捕获率升高的同时,误判率也随之上升。以log(q/T)≤-3.1作为预测滑坡的判别标准,模型效果较好,预测捕获率为62.50%,误判率(17.79%)较低。预测结果显示,滑坡潜在区域主要位于斜坡下部、土体厚度大和坡度陡峭的地区,山体顶部、土体厚度薄和地形平坦的区域斜坡稳定。