Geological causes and geomorphological precursors of the Tsaoling landslide triggered by the 1999 Chi-Chi earthquake, Taiwan

The Tsaoling landslide, one of the largest landslide areas in Taiwan, has been affected by catastrophic events triggered by rain or earthquakes six times since 1862. These landslides, including that caused by the 1999 earthquake, have essentially not been reactivated old slides, but were sequential new ones that developed upslope, retrogressively. The landslide area is underlain by Pliocene sandstone and shale to form a dip slope with a bedding plane, dipping uniformly at 14°. The slip surface of the 1999 landslide was smooth and planar, parallel to the bedding plane with a slightly stepped profile; it formed within thinly alternated beds of fine sandstone and shale with ripple lamination or in a shale bed. The shale is weathered by slaking and probably by sulfuric acid, which is inferred to be one of the major causes of the intermittent retrogressive development of the landslides. The weathering was likely accelerated by the removal of overlying beds during earlier landslides in 1941 and 1942. The top margin of the 1999 landslide, in plan view, coincided with a V-shaped scarplet, which can be clearly recognized on aerial photographs taken before the landslide. This geomorphological feature indicates that this landslide had already moved slightly before its 1999 occurrence, providing precursory evidences.     

Geophysical characterization of the lithological control on the kinematic pattern in a large clayey landslide (Avignonet,French Alps)

Lithology variation is known to have a major control on landslide kinematics, but this effect may remain unnoticed due to low spatial coverage during investigation. The large clayey Avignonet landslide (French Alps) has been widely studied for more than 35 years. Displacement measurements at 38 geodetic stations over the landslide showed that the slide surface velocity dramatically increases below an elevation of about 700 m and that the more active zones are located at the bottom and the south of the landslide. Most of the geotechnical investigation was carried out in the southern part of the landslide where housing development occurred on lacustrine clay layers. In this study, new electrical prospecting all across the unstable area revealed the unexpected presence of a thick resistive layer covering the more elevated area and overlying the laminated clays, which is interpreted as the lower part of moraine deposits. The downslope lithological boundary of this layer was found at around 700 m asl. This boundary coincides with the observed changes in slide velocity and in surface roughness values computed from a LiDAR DTM acquired in 2006. This thick permeable upper layer constitutes a water reservoir, which is likely to influence the hydromechanical mechanism of the landslide. The study suggests a major control of vertical lithological variations on the landslide kinematics, which is highlighted by the relation between slide velocity and electrical resistivity.     

Movement of deep-seated rainfall-induced landslide at Hsiaolin Village during Typhoon Morakot

This study analyzes the mechanism of the landslide event at Hsiaolin Village during Typhoon Morakot in 2009. This landslide event resulted in 400 deaths. The extremely high intensity and accumulative rainfall events may cause large-scale and complex landslide disasters. To study and understand a landslide event, a combination of field investigations and numerical models is used. The landslide area is determined by comparing topographic information from before and after the event. Physiographic parameters are determined from field investigations. These parameters are applied to a numerical model to simulate the landslide process. Due to the high intensity of the rainfall event, 1,675 mm during the 80 h before the landslide event, the water content of soil was rapidly increased causing a landslide to occur. According to the survivors, the total duration of the landslide run out was less than 3 min. Simulation results indicated that the total duration was about 150 s. After the landslide occurrence, the landslide mass separated into two parts by a spur at EL 590 in about 30 to 50 s. One part passed the spur in about 30 to 60 s. One part inundated the Hsiaolin Village and the other deposited at a local river channel and formed a landslide dam. The landslide dam had height between 50 and 60 m and length between 800 and 900 m. The simulation result shows that the proposed model can be used to evaluate the potential areas of landslides induced by extremely high intensity rainfall events.     

Rainfall-induced landslide in the active frontal fold–thrust belt of Northwestern Himalaya,Jammu: dynamics inferred by geological evidences and Ground Penetrating Radar

Landslides are the most established geological hazards in the frontal fold–thrust belt of Northwestern Himalaya comprising of Siwaliks and Murree strata. The continuous rainfall from 2 to 6 September, 2014 caused a massive landslide at village Sadal in Udhampur district of Jammu and Kashmir state. The landslide occurred in the early morning of September 6, 2014, destroying entire Sadal habitation comprising 45 houses, and killing 41 people and more than 500 domestic animals. Google earth images of pre and post-landslide events along with the field measurements show the kinematics of upper and lower parts of the slide. Horizontal and vertical components of displacement and mode of failure suggest the landslide as of complex nature. The shallow subsurface geophysical imaging through Ground Penetrating Radar (GPR) survey shows the failure plane composed of friable mudstone bed underlain by massive mudstone and overlain by cross-bedded sandstone. The depth of debris material above the failure plane ranges from 6 m at Site S1a-b to 10 m at Site-S2b and 20 m at Site S3a. The velocity analysis of Site-3 shows four thick layers represented from bottom to surface by L1—sandstone (V?=?0.16 m/ns, travel time?=?356.36 ns), L2—mudstone (V?=?0.17 m/ns, travel time?=?288.48 ns), L3—massive mudstone (V?=?0.19 m/ns, travel time 220.68 ns), and L4—cross-laminated sandstone (V?=?0.20 m/ns, travel time?=?77.58 ns) overlaying the failure plane. The study shows the landslide occur along the western limb of a fold identified during the present work. We mapped an old landslide on the same limb which shows 5–6 m-thick subsurface debris material with thick rock fragments involved in the landslide process. The detailed geological and geophysical investigations suggest that both the landslides were triggered by extreme rain fall events.     

Field monitoring of groundwater responses to heavy rainfalls and the early warning of the Kualiangzi landslide in Sichuan Basin,southwestern China

The Kualiangzi landslide was triggered by heavy rainfalls in the “red beds” area of Sichuan Basin in southwestern China. Differing from other bedrock landslides, the movement of the Kualiangzi landslide was controlled by the subvertical cracks and a subhorizontal bedding plane (dip angle < 10°). The ingress of rainwater in the cracks formed a unique groundwater environment in the slope. Field measurement for rainfall, groundwater movement, and slope displacement has been made for the Kualiangzi landslide since 2013. The field monitoring system consists of two rainfall gauges, seven piezometers, five water-level gauges, and two GPS data loggers. The equipments are embedded near a longitudinal section of the landslide, where severe deformation has been observed in the past 3 years. The groundwater responses to four heavy rainfall events were analyzed between June 16 and July 24 in 2013 coincided with the flood season in Sichuan. Results showed that both of the water level and the pore-water pressure increased after each rainfall event with delay in the response time with respect to the precipitation. The maximum time lag reached 35 h occurred in a heavy rainfall event with cumulative precipitation of 127 mm; such lag effect was significantly weakened in the subsequent heavy rainfall events. In each presented rainfall event, longer infiltration period in the bedrock in the upper slope increased the response time of groundwater, compared to that of in the gravels in the lower slope. A translational landslide conceptual model was built for the Kualiangzi landslide, and the time lag was attributed to the gradual formation of the uplift pressure on the slip surface and the softening of soils at the slip surface. Another important observation is the effect on the slope movement which was caused by the water level (H _w) in the transverse tension trough developed at the rear edge of the landslide. Significant negative correlation was found for H _w and the slope stability factor (F _s), in particular for the last two heavy rainfall events, of which the drastic increase of water level caused significant deterioration in the slope stability. The rapid drop (Δ?=?22.5 kPa) of pore-water pressure in the deep bedrock within 1 h and the large increase (Δ?=?87.3 mm) of surficial displacement were both monitored in the same period. In the end, a four-level early warning system is established through utilizing H _w and the displacement rate D _r as the warning indicators. When the large deformation occurred in flood season, the habitants at the leading edge of the landslide can be evacuated in time.     

The July 1, 2017 Wangjiawan landslide in Ningxiang County,China

Many landslides were triggered by heavy rainfall from 29 June to 1 July 2017 in Ningxiang County, Hunan Province, China. A field investigation into one of the landslides, known as the Wangjiawan landslide, was undertaken on July 3, 2017, to understand the landslide mechanisms and the factors that triggered the event. The landslide is a translational and rotational slide that degraded downslope to an earth flow. It occurred on a steeply dipping mud-rich slate intercalation. Field investigation shows that the landslide had a movement rate of more than 25 m/s and had been triggered by up to 338 mm of continuous rainfall over 2 days. The landslide was fast-moving because it occurred on a steep slope and showed a rapid reduction in shear strength. The landslide resulted in nine fatalities and 19 injuries, most during initial rescue efforts. To mitigate a possible secondary landslide disaster, later emergency measures including evacuation of the survivors, setting up warning signs, and covering landslide cracks with plastic sheeting were adopted.     

Numerical modeling of the June 24, 2015, Hongyanzi landslide generated impulse waves in Three Gorges Reservoir,China

Subaerial landslides falling into confined water bodies often generate impulsive waves. Damaging landslide tsunamis in Three Gorges Reservoir, China, have struck several times in the last 15 years. On June 24, 2015, a 23?×?10⁴ m³ slope failure occurred on the east bank of the Daning River opposite Wushan Town. The sliding mass intruded into the Three Gorges Reservoir and initiated a reservoir tsunami that resulted in two deaths and significant damage to shipping facilities. A post-event survey revealed the landslide geometry and wave run-up distribution, while an eyewitness video captured most of the landslide motion. Employing these firm constraints, we applied the Tsunami Squares method to simulate the 2015 Hongyanzi landslide and tsunami. The simulation revealed that the landslide experienced a progressive failure in the first few seconds and impacted the water with a maximum velocity of ~?16 m/s. The initial wave propagated to the opposite shore in an arch shape, and the water surface reached a maximum amplitude of ~?11 m near the landslide. Wave amplitude-time curves at four points on the river cross section show that the initial wave reached Wushan town in about 50 s with an average wave velocity of ~?30 m/s. The maximum wave run-ups on the shoreline opposite the landslide are around 6 m and attenuate to less than 1 m beyond 2-km distance. The landslide simulation matches the observed geological profile and the eyewitness video, and the numerical results coincide with the observed wave run-up heights. Nearly 80% of landslide energy is lost due to frictional resistances, but the remaining fraction imparted to the tsunami carried catastrophic consequences to a large region. The numerical results emphasize the efficiency and accuracy of Tsunami Squares method for a “Quick Look” simulation of a potential landslide.     

Dynamic process of the Wenjiagou rock landslide in Sichuan Province,China

A rock avalanche is a geological event that is always sudden, rapid and with a long run-out, and can result in large loss of lives and property. The Wenjiagou rock avalanche was a high-speed rock landslide caused by a strong earthquake, in Mianzhu, Sichuan Province, southwest China. In this study, we reproduce the movement and deposition processes of the sliding mass by numerical simulation. We analyze the effects of the friction coefficient of each slip surface and the strength of the parallel bonds and contact stiffness between particles on the dynamic process and deposit features using three-dimensional particle flow code (PFC3D). The simulation results agree with the field measurements when the friction coefficient is 0.2, parallel bond strength is 2 MPa, and contact stiffness is 2?×?10⁸ kN/m. The landslide lasted about 115 s from the initial movement to the final deposition at the exit of the valley. The maximum velocity of the sliding mass was 114 m/s.     

Catastrophic rockslide-debris avalanche at St. Bernard, Southern Leyte, Philippines

On 17 February 2006, a rockslide-debris avalanche cascaded down the steep slope of Mt. Can-abag, burying the entire village of Guinsaugon in St. Bernard, Southern Leyte, Philippines. Casualties include 139 dead with 980 still missing and presumed dead, making it perhaps the most catastrophic landslide in Philippine history. The landslide started at the ridge top along a fault plane associated with the active Philippine Fault Zone. It started as a block slide that transformed into an avalanche. The entire event lasted for only a few minutes. Estimated maximum landslide velocity is 120–130 m/s. The landslide left behind a deep, wedge-shaped scarp. The central part of the deposit exhibits a hummocky topography typical of avalanches, in contrast to the flatter surface of the debris-flow-type marginal deposit. High amounts of soil in the matrix were derived from the scouring of ancient landslide deposits and rice fields in the valley. The landslide has a total area of 3.2 km² and a runout distance of 4.1 km. Estimated volume of debris is approximately 20 Mm³. At least four streams were dammed by the landslide debris. Intense precipitation and earthquakes preceding the landslide are the potential triggers. Preliminary back analyses assuming a planar and wedge slip surface yielded very low factors of safety even under dry conditions. A more rigorous analysis of the failure mechanism of the landslide is needed.     

Analysis of Baishuihe landslide influenced by the effects of reservoir water and rainfall

Since the impoundment of the Three Gorges Reservoir in June 2003, a number of new landslides have occurred and existing landslides have been made worse. The 1,260 × 10⁴ m³ Baishuihe landslide, located at 56 km west of the Three Gorges Dam, began to deform more noticeably after the first impoundment in early July 2003. The sliding of the two blocks comprising the landslide, one an active block and the other a relatively stable block, became apparent after approximately 5 years of monitoring. Field recordings show that the landslide displacement is affected by the combined effects of the rainfall and water level in the reservoir. These effects have been investigated in the present paper, including the deformation characteristics (movement pattern, direction, displacement and velocity) earmarking the temporal evolution of the active block. Based on a practical creep model of a large rock slide, alert velocity thresholds for pre-alert, alert and emergency phases have been computed corresponding to the imminence of failure. The alert velocity thresholds are being proposed to be included as a part of an early-warning system of an emergency plan drawn up to minimize the adverse impact in the event of landslide failure. The emergency plan is intended to be implemented as a risk management tool by the relevant authorities of the Three Gorges Reservoir in the near future.     

Landslide at Su-Hua Highway 115.9k triggered by Typhoon Megi in Taiwan

This study focused on the landslide case at Su-Hua Highway 115.9k, Taiwan. A preliminary investigation was conducted on geomorphologic features change and landslide mechanisms using digital elevation models, geographical maps, and remote sensing images at different times in conjunction with geological surveys and analysis results. Using the results of geological surveys and physical model experiments, we constructed a discrete element method to simulate the process of landslide movement. The results revealed deformation in the metamorphic rock slopes upstream of 115.9k. The slopes around the erosion gully upstream presented visible slope toes cutting and tension cracks at the crest as well as unstable rock masses. According to the results of numerical simulation for typhoon Megi event, intense rains could induce slippage in the rock debris/masses in the source area, initially at a speed of 5–20 m/s. Subsequently, steeper terrain could cause the rock debris/masses to accelerate to form a high-speed (>30 m/s) debris slide quickly moving downstream to form an alluvial fan downstream by the sea.     

Estimation of dynamic friction and movement history of large landslides

We performed seismic waveform inversions and numerical landslide simulations of deep-seated landslides in Japan to understand the dynamic evolution of friction of the landslides. By comparing the forces obtained from a numerical simulation to those resolved from seismic waveform inversion, the coefficient of friction during sliding was well-constrained between 0.3 and 0.4 for landslides with volumes of 2–8 ×10⁶ m³. We obtained similar coefficients of friction for landslides with similar scale and geology, and they are consistent with the empirical relationship between the volume and dynamic coefficient of friction obtained from the past studies. This hybrid method of the numerical simulation and seismic waveform inversion shows the possibility of reproducing or predicting the movement of a large-scale landslide. Our numerical simulation allows us to estimate the velocity distribution for each time step. The maximum velocity at the center of mass is 12–36 m/s and is proportional to the square root of the elevation change at the center of mass of the landslide body, which suggests that they can be estimated from the initial DEMs. About 20% of the total potential energy is transferred to the kinetic energy in our volume range. The combination of the seismic waveform inversion and the numerical simulation helps to obtain the well-constrained dynamic coefficients of friction and velocity distribution during sliding, which will be used in numerical models to estimate the hazard of potential landslides.     

Deformation and mechanism of landslide influenced by the effects of reservoir water and rainfall, Three Gorges, China

Since the impoundment of the Three Gorges Reservoir in June 2003, numerous preexisting landslides have been reactivated. This paper seeks to find the factors influencing landslide deformation and the relationship between displacement and fluctuation of the reservoir water level, while the displacement and the intensity of rainfall based on monitoring data; 6 years of monitoring were carried out on the Shiliushubao landslide, a old landslide, consisting of a deep-seated main block and two shallow blocks, with a volume of 1,180 × 10⁴ m³ and located on the left bank of the Yangtze River, 66 km upstream of the Three Gorges dam. This landslide was reactivated by the impoundment and since then the landslide body has been experiencing persistent deformation with an observed maximum cumulative displacement of 8,598.5 mm up to December 2009. Based on the monitoring data, we analyzed the relationship between the fluctuation of the reservoir water level and displacement, rainfall and displacement, and found that the rainfall is the major factor influencing deformation for two shallow blocks and the displacement has a positive correlation with the variation of rainfall intensity. The fluctuation of the reservoir water level is the primary factor for main block, and the deformation rate has a negative correlation with the variation of reservoir water level, declined with the rise of the water level and increased with the drawdown of the water level.     

四川茂县周场坪滑坡发育特征与变形监测分析

位于四川茂县南新镇的周场坪滑坡是一大型古滑坡,曾于1982年发生大规模快速复活,目前滑体半堰塞岷江。野外地质调查表明,周场坪滑坡在平面上呈不规则长舌形,长约850 m,滑体前后缘高差约350 m;在剖面上发育3级滑动,钻探揭露滑带埋深以50~70 m为主,推测潜在失稳滑坡体积为1 500×10 ⁴~2 000×10 ⁴ m ³。周场坪滑坡在平面上分为4个变形区,在滑体中后部和前缘坡脚发育大量拉裂缝与下错陡坎,拉裂缝宽度以0.2~3.0 m为主,陡坎下错高度2~10 m。在野外调查和钻探分析的基础上,对该滑坡开展了地表位移(GNSS)、深部蠕滑变形(钻孔测斜仪)和雨量等监测。监测分析表明,目前周场坪滑坡后缘变形速率达0.80 m/a,中部和前缘分别为0.69 m/a和0.51 m/a,呈推移式滑动变形,整体向NW310°方向滑动;地表位移速率在监测期内基本在1~3 mm/d之间波动,波动主要受降雨量影响,且略滞后于降雨量;滑移加速度基本在0~6 mm/d ²范围波动。ZK2钻孔测斜仪监测数据表明,滑坡在80 m深度内主要沿2层滑带蠕滑,其中浅层滑带埋深在22 m左右,深层滑带埋深在46 m左右,滑移速率在0~5 mm/d范围波动。综合研究认为,周场坪滑坡目前处于缓慢变形的深层蠕滑中,其变形速率受降雨和河流侵蚀等因素的影响,在极端内、外动力条件下,可能会加速滑动,并再次整体复活,造成堵塞岷江等危害。     

Deep-seated gravitational slope deformations triggered by extreme rainfall and agricultural practices (eastern Michoacan,Mexico)

The study of deep-seated gravitational slope deformations (DSGSD) in Mexico is scarce; therefore, their localization and causes are highly overlooked. The present paper examines the characterization of the DSGSD of Jungapeo and Las Pilas in eastern Michoacan state, currently active and endangering their inhabitants. An integrated study, including detailed lithology, morpho-structural inventories, analysis of land use, and pluviometric regime, was performed and complemented with differential global positioning system monitoring networks. Both landslides developed over highly weathered volcano-sedimentary rocks. On the one hand, the Jungapeo landslide has an estimated volume of 990,455 m³ with steady decreasing velocity rates from 41 to 15 cm/month in the first monitoring period to 13–3 cm/month in the last one. On the other hand, the Las Pilas landslide estimated volume is about 1,082,467 m³ with a stable velocity rate of 1.3 to 0.1 cm/month. Despite the multi-storeyed style of activity, two behaviors of instability were distinguished: slow deformation and secondary landslide stages. The conditioning factors for slow deformation in both DSGSD are the combination of weathered lithology with clay- and sand-rich content, and the shift toward intensive monoculture. The triggering factor is related to excess water produced by an inefficient flood-irrigation system that also generates an atypical acceleration behavior in both landslides during the dry season. The DSGSD activity thus predisposes the generation of tension cracks and secondary scarps from which the collateral landslides are triggered by atypical rainfall, such as that of 2010.     

Slope dynamics and streambed uplift during the Pergalla landslide reactivation in March 2016 and discussion of concurrent causes (Northern Apennines,Italy)

On March 28, 2016, the toe zone of the apparently dormant Pergalla earthslide-earthflow (Northern Apennines, Italy) had a paroxysmal reactivation. In the course of 2 days, displacements up to almost 8 m severely damaged several houses and roads. At the bottom of the slope, the emersion of rotational sliding surfaces determined the uplift of almost 3 m of the Nure river streambed that was consequently partially dammed. The paper describes the landslide event on the basis of field surveys and analysis of post-event aerial photos, as well as data from geophysical surveys and pre- to post-failure displacement monitoring. It also discusses the possible concurrent causes of the event, including antecedent rainfall, the migration of active streambed channels of Nure river toward the landslide toe in the previous year, and the existence of long-term pre-failure slow movements. It is concluded that these factors, together with the presence of sliding surfaces extending beneath the valley floor, should be primarily considered if a preventive assessment of river damming potential due to streambed uplift should be made for other similar landslides in the Apennines.     

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.     

Dynamic process simulation of rainfall-induced Sanxicun landslide based on a thermo-poro-elastic approach

The Sanxicun landslide occurred on July 10, 2013, in Sanxicun Village, which is located in Dujiangyan City, Sichuan Province, China. It travelled up to 1200 m, destroyed 11 houses, and killed 166 people in the village. To explain how this landslide could travel such a large distance and cause such serious damage, this study used a thermo-poro-elastic approach coupled with the Savage–Hutter model to simulate the dynamic process of the Sanxicun landslide. The simulated results were compared with the actual results, as well as those of other researchers. It showed that the simulated results for the landslide profile and mass accumulation scope were basically consistent with the actual results. The simulated landslide runout was 1242 m, which was quite close to the actual value. The simulated maximum mass accumulation thickness was 16.4 m. The maximum velocity was 32.6 m/s, which was between those calculated by Yin et al. (J Eng Geol 22(2):309–318, 2014), Yin et al. (Landslide 13:9–23, 2016), and the various trends were found to be consistent. The temperature change and the pore water pressure evolution in the shear zone during sliding are also obtained by simulation. This study had recreated the Sanxicun landslide motion process from the view of thermo-poro-elastic coupling within the shear zone.