Significance of geomorphological and subsurface drainage controls on failures of peat‐covered hillslopes triggered by extreme rainfall

On 19 September 2003, 40 landslides of 140–18 000 m³ volume occurred within 2·5 km² on the slopes of Dooncarton Mountain (Republic of Ireland) during a storm that may have exceeded 90 mm within 90 minutes. The landslides were investigated to determine the reasons for such a high density of slope failures. All of the landslides were surveyed within four months, and nine of them were investigated in detail. The six largest landslides, all peat failures, accounted for 57% of the more than 100 000 m³ of material displaced during the event. A consistent sequence of superficial materials was found on the failed hillslopes, including an extensive iron pan at the base of a buried soil horizon 0·3 m below the base of the peat. Morphologically, almost all of the landslides occurred on steep planar slopes or around sharp convexities, with the latter failures developing retrogressively upslope. The only significant relationship found from analysis of 371 subsurface pipes and 142 seepage cracks (defined here as contiguous fissures conducting concentrated subsurface flow) across all the failures was that the thinner the peat cover, the deeper the pipes and seepage cracks occurred below the base of peat. It is concluded that most of the landslides were probably caused by a combination of excess water pressures in the buried soil horizon and the thinner overburden of peat or peaty soil associated with the steeper slope segments. Pipes and seepage cracks formed on the iron pan probably existed prior to the failure event and may have contributed to the high water pressures as rainwater inputs exceeded their discharge capacities. One large peat slide was probably triggered by excess water pressures developed within and between artificial tine cuts. The properties of the blanket peat were generally of little consequence in the occurrence of the landslides, but relict desiccation cracks and other structural weaknesses through the peat mass were probably highly significant. Although several aspects of the peat failures correspond to previously published examples, the context of these failures in terms of the topography and upland catena is distinctive. Copyright © 2007 John Wiley & Sons, Ltd.     

Seismic loading impacts on excess pore‐water pressure maintain landslide triggered flowslides

During the 2003 Sanriku‐Minami earthquake, Japan, a flowslide was triggered on a slope of about 13.5º. The displaced landslide mass developed into a flowslide and deposited on a horizontal rice paddy after traveling approximately 130 m. To study the trigger and movement mechanisms of this landslide, field investigation and laboratory ring‐shear tests were performed. Field investigation revealed that the landslide originated from a fill slope, where a gully was buried for cultivation some decades ago, and shallow ground water was present. Undrained monotonic and cyclic ring‐shear tests on a sample (pyroclastic deposits) taken from the source area revealed that the soil is highly liquefiable, and its steady‐state shear strength can be little affected by overconsolidation. Using the seismic records of the earthquake, probable seismic loadings on the sliding surface were synthesized and applied to the samples in ring‐shear tests, which were performed under undrained or partially drained conditions. The undrained and partially drained tests revealed that shear failure can be triggered by the introduction of seismic loading and formation of excess pore‐water pressure. The generation of excess pore‐water pressure along with increase of shear displacement and the inhibited dissipation of excess pore‐water pressure due to the thickness of the saturated soil layer above the sliding surface probably enabled the continued post‐failure landsliding. Copyright © 2008 John Wiley & Sons, Ltd.     

Seepage erosion and its implication to the formation of amphitheatre valley heads: A case study at Obara,Japan

Seepage erosion was investigated in an amphitheatre with a semicircular valley head, steep slopes, and a flat bottom developed in granodiorite hills at Obara, Aichi prefecture, Japan. A high sediment yield occurred where the measuring sites were located at the base of the landslide debris in the base of the convex slopes, whereas sediment outflows were small where the measuring sites were located at the base of the strong convex slopes. This implies that the seepage erosion was an effective agent for removal of debris deposited at the base of the slope. Small landslides can be found at the lower slopes within the area of the observed amphitheatre. The slope stability analysis and subsurface water observation of the lower slope suggest that the small landslides in this amphitheatre are due to over-steepened slopes, and relatively insensitive to subsurface water status. Colluvium in the flat valley bottom thinly covers the bedrock surface. Therefore the topography of the amphitheatre was found to be formed by parallel retreat of slopes by the repetition of basal seepage erosion and subsequent small landslides.     

Snowmelt runoff and soil moisture dynamics on steep subalpine hillslopes

Snowmelt water supplies streamflow and growing season soil moisture in mountain regions, yet pathways of snowmelt water and their effects on moisture patterns are still largely unknown. This study examined how flow processes during snowmelt runoff affected spatial patterns of soil moisture on two steep sub‐alpine hillslope transects in Rocky Mountain National Park, CO, USA. The transects have northeast‐facing and east‐facing aspects, and both extend from high‐elevation bedrock outcrops down to streams in valley bottoms. Spatial patterns of both snow depth and near‐surface soil moisture were surveyed along these transects in the snowmelt and summer seasons of 2008–2010. To link these patterns to flow processes, soil moisture was measured continuously on both transects and compared with the timing of discharge in nearby streams. Results indicate that both slopes generated shallow lateral subsurface flow during snowmelt through near‐surface soil, colluvium and bedrock fractures. On the northeast‐facing transect, this shallow subsurface flow emerged through mid‐slope seepage zones, in some cases producing saturation overland flow, whereas the east‐facing slope had no seepage zones or overland flow. At the hillslope scale, earlier snowmelt timing on the east‐facing slope led to drier average soil moisture conditions than on the northeast‐facing slope, but within hillslopes, snow patterns had little relation to soil moisture patterns except in areas with persistent snow drifts. Results suggest that lateral flow and exfiltration processes are key controls on soil moisture spatial patterns in this steep sub‐alpine location. Copyright © 2014 John Wiley & Sons, Ltd.     

Seepage erosion in layered stream bank material

Current stream restoration practices often require anthropogenic manipulation of natural field soils to reconstruct stream banks in the absence of stabilizing vegetation. For this study, researchers conducted laboratory experiments on reconstructed, non‐vegetated stream banks with layered soils experiencing seepage. The objective of the study was to determine the effect of seepage, pore water pressure, and bank geometry on erosion and bank stability of layered streambanks. The experimental design consisted of an intermediate‐size soil lysimeter packed with a sandy clay loam top soil and an underlying fine sand layer at three bank slopes (90°, 45° and 26°). Shallow groundwater flow and seepage resulted in bank failure of geometrically stable banks. Pop out failures, liquid deformation, and piping were all observed failure mechanisms in the underlying sand material, dependent on the bank angle. Groundwater seepage processes created small‐scale failures of the underlying sand leading to larger‐scale failures of the overlying sandy clay loam. The underlying sand layer eroded according to the initial bank angle and change in overburden loading. The overlying loam layer failed along linear failure planes. The gradually sloped bank (i.e. 26° slope) failed faster, hypothesized to be due to less confining pressure and greater vertical seepage forces. Researchers analyzed the laboratory experiments using the Bank Stability and Toe Erosion Model, version 4·1. The model calculated an accurate shear surface angle similar to the failure angle observed in the lysimeter tests. The model predicted failure only for the undercut 90° bank slope, and indicated stable conditions for the other geometries. Steeper initial bank slopes and undercut banks decreased the bank factor of safety. The observed failure mechanisms and measured saturation data indicated an interaction between overburden pressure, seepage forces, and bank slope on bank stability. Future bank stability modeling would benefit by incorporating lateral seepage erosion and soil liquefaction prediction calculations. Copyright © 2009 John Wiley & Sons, Ltd.     

Temporal and spatial variation of infilling processes in a landslide scar in a steep mountainous region,Japan

The duration of the soil‐depth recovery needed for reoccurrence of shallow colluvial landslides at a given site in humid regions is much longer than the return period of rainfall needed to generate sufficient pore water pressure to initiate a landslide. Knowledge of the rate of change in soil depth in landslide scars is therefore necessary to evaluate return intervals of landslides. Spatial variation in sediment transport at the Kumanodaira landslide scar in central Japan was investigated by field observations. Spatial distribution of the rate of change in soil depth was estimated using sediment transport data and geographic information system (GIS) analysis. Observations revealed that the timing of sediment transport differed for shallow and deep soil layers. Near‐surface sediment transport (mostly dry ravel and some shallow soil creep at depths ≤0·05 m) measured in sediment traps was active in winter and early spring and was affected by freezing–thawing; soil creep of subsoil (i.e. >0·05 m), monitored by strain probes, was active in summer and autumn when precipitation was abundant. Near‐surface sediment flux was estimated by a power law function of slope gradient. Deeper soil creep was more affected by relative location to the landslide scar, which influences soil depth, than by slope gradient. Our study indicated that the rate of soil‐depth recovery is high just below the head scarp of the landslide. Abrupt changes in the longitudinal slope topography immediately above, within and just below the head scarp became smoother with time due to degradation proximate to the landslide head scarp and flanks, as well as aggradation just below the head scarp. Copyright © 2014 John Wiley & Sons, Ltd.     

Testing a model for predicting the timing and location of shallow landslide initiation in soil‐mantled landscapes

The growing availability of digital topographic data and the increased reliability of precipitation forecasts invite modelling efforts to predict the timing and location of shallow landslides in hilly and mountainous areas in order to reduce risk to an ever‐expanding human population. Here, we exploit a rare data set to develop and test such a model. In a 1·7 km² catchment a near‐annual aerial photographic coverage records just three single storm events over a 45 year period that produced multiple landslides. Such data enable us to test model performance by running the entire rainfall time series and determine whether just those three storms are correctly detected. To do this, we link a dynamic and spatially distributed shallow subsurface runoff model (similar to TOPMODEL) to an in?nite slope model to predict the spatial distribution of shallow landsliding. The spatial distribution of soil depth, a strong control on local landsliding, is predicted from a process‐based model. Because of its common availability, daily rainfall data were used to drive the model. Topographic data were derived from digitized 1 : 24 000 US Geological Survey contour maps. Analysis of the landslides shows that 97 occurred in 1955, 37 in 1982 and ?ve in 1998, although the heaviest rainfall was in 1982. Furthermore, intensity–duration analysis of available daily and hourly rainfall from the closest raingauges does not discriminate those three storms from others that did not generate failures. We explore the question of whether a mechanistic modelling approach is better able to identify landslide‐producing storms. Landslide and soil production parameters were ?xed from studies elsewhere. Four hydrologic parameters characterizing the saturated hydraulic conductivity of the soil and underlying bedrock and its decline with depth were ?rst calibrated on the 1955 landslide record. Success was characterized as the most number of actual landslides predicted with the least amount of total area predicted to be unstable. Because landslide area was consistently overpredicted, a threshold catchment area of predicted slope instability was used to de?ne whether a rainstorm was a signi?cant landslide producer. Many combinations of the four hydrological parameters performed equally well for the 1955 event, but only one combination successfully identi?ed the 1982 storm as the only landslide‐producing storm during the period 1980–86. Application of this parameter combination to the entire 45 year record successfully identi?ed the three events, but also predicted that two other landslide‐producing events should have occurred. This performance is signi?cantly better than the empirical intensity–duration threshold approach, but requires considerable calibration effort. Overprediction of instability, both for storms that produced landslides and for non‐producing storms, appears to arise from at least four causes: (1) coarse rainfall data time scale and inability to document short rainfall bursts and predict pressure wave response; (2) absence of local rainfall data; (3) legacy effect of previous landslides; and (4) inaccurate topographic and soil property data. Greater resolution of spatial and rainfall data, as well as topographic data, coupled with systematic documentation of landslides to create time series to test models, should lead to signi?cant improvements in shallow landslides forecasting. Copyright © 2003 John Wiley & Sons, Ltd.     

Slope control on the frequency distribution of shallow landslides and associated soil properties,North Island,New Zealand

The extrapolation of results from field trials to larger areas of land for purposes of regional impact assessment is an important issue in geomorphology, particularly for landform properties that show high stochastic variability in space and time, such as shallow landslide erosion. It is shown in this study, that by identifying the main driver for spatial variability in shallow landslide erosion at field scales, namely slope angle, it is possible to develop a set of generic functions for assessing the impact of landslides on selected soil properties at larger spatial scales and over longer time periods. Research was conducted within an area of pastoral soft‐rock Tertiary hill country in the North Island of New Zealand that is subject to infrequent high intensity rainfall events, producing numerous landslides, most of which are smaller than several hundred square metres in size and remove soil to shallow depths. All landslides were mapped within a 0·6 km² area and registered to a high resolution (2 m) slope map to show that few landslides occur on slopes < 20° and 95% were on slopes > 24°. The areal density of landslides from all historical events showed an approximately linear increase with slope above 24°. Integrating landslide densities with soil recovery data demonstrates that the average value of a soil property fluctuates in a ‘saw‐tooth’ fashion through time with the overall shape of the curve controlled by the frequency of landslide inducing storm events and recovery rate of the soil property between events. Despite such fluctuations, there are gradual declines of 7·5% in average total carbon content of topsoil and 9·5% in average soil depth to bedrock, since the time of forest clearance. Results have application to large‐scale sediment budget and water quality models and to the New Zealand Soil Carbon Monitoring System (CMS). Copyright © 2012 John Wiley & Sons, Ltd.     

基于人工边坡稳定性的降雨型滑坡预警预报模型设计与应用

由削坡建房遗留的人工边坡存在大量滑坡隐患问题,在降雨引发土质边坡自身动力变化分析条件下,以稳定性评价建模为基础,提出降雨型滑坡动力学预警预报模型。文中以广东省梅州市花岗岩地区为例,使用GIS技术构建了1 727个预警分析单元,并进行关键地质环境因子赋值及与气象站点数据关联;按坡高、坡度等参数,分别构建16个边坡失稳动力学预警模块,并根据降雨量变化,计算边坡稳定性系数,最终按其阈值确定风险等级并予以预警。本研究对于推动人工边坡诱发的滑坡地质灾害预警预报与预防均具有重要意义。     

The influence of water retention curve hysteresis on the stability of unsaturated soil slopes

Hysteresis is a common feature exhibited in hydraulic properties of an unsaturated soil. The movement of wetting front and the hysteresis effect are important factors which impact the shear strength of the unsaturated soil and the mechanics of shallow landslides. These failures are mainly triggered by the deepening of the wetting front accompanied by a decrease in matric suction induced by infiltration. This research establishes a method for determining a stability analysis of unsaturated infinite soil slopes, integrating the influence of infiltration and the water retention curve hysteresis. Furthermore, the present stability analysis method including the infiltration model and the advanced Mohr–Coulomb failure criterion calculates the variations of the safety factor (FS) in accordance with different slope angle, depth and hydrological processes. The experimentally measured data on the effect of hysteresis are also carried out for comparison. Numerical analyses, employing both wetting and drying hydraulic behaviour of unsaturated soil, are performed to study the difference in soil‐water content as observed in the experiments. The simulating approximations also fully responded to the experimental data of sand box. The results suggest that the hysteresis behaviour affect the distribution of soil‐water content within the slope indeed. The hysteresis made the FS values a remarkable recovery during the period of non‐rainfall in a rainfall event. The appropriate hydraulic properties of soil (i.e. wetting or drying) should be used in accordance with the processes that unsaturated soil actually experience. This method will enable us to acquire more accurate matric suction head and the unsaturated soil‐shear strength as it changes with the hysteretic flow, in order to calculate into the stability analysis of shallow landslides. An advanced understanding of the process mechanism afforded by this method is critical to realizing a reliable and appropriate design for slope stabilization. It also offers some immediate reference information to the disaster reduction department of the government. Copyright © 2011 John Wiley & Sons, Ltd.     

Landslide activation behaviour illuminated by electrical resistance monitoring

A common factor in landslide activation (or reactivation) is subsurface moisture and associated pore pressure variations linked to rainfall. Monitoring of these subsurface hydrogeological processes is necessary to improve our understanding of water‐induced landslide activation. Geophysical approaches, electrical methods in particular, are increasingly being applied to landslide monitoring because they provide non‐invasive spatial information in heterogeneous subsurface environments that can be difficult to characterise using surface observations or intrusive sampling alone. Electrical techniques are sensitive to changing subsurface moisture conditions, and have proven to be a useful tool for investigating the hydrogeology of natural and engineered slopes. The objectives of this investigation were to further develop electrical resistance monitoring for slope stability assessment, and to validate the approach at an intermittently‐active UK landslide system to advance the understanding of complex landslide activation mechanisms. A long‐term transfer resistance dataset was collected from a grid of electrodes to allow spatial monitoring of the landslide. These data were interpreted using a synthesis of rainfall, temperature, GPS and piezometric records. The resistance data were corrected for seasonal temperature variations and electrode movements were monitored, as these processes were shown to mask moisture related changes. Results reveal that resistance monitoring is sensitive to soil moisture accumulation, including changes in piezometric levels, and can be used to study the principal activation mechanism of slow‐moving shallow earthflows. Spatial monitoring using resistance maps was shown to be particularly valuable as it revealed the evolution of subsurface moisture distribution, in the lead up to landslide activation. Key benefits of this approach are that it provides a simple, rapid and non‐invasive means of spatially monitoring subsurface moisture dynamics linked to landslide activation at high‐temporal resolution. Crucially, it provides a means of monitoring subsurface hydraulic changes in the build‐up to slope failure, thereby contributing to early warning of landslide events. Copyright © 2018 John Wiley & Sons, Ltd.     

How does forest structure affect root reinforcement and susceptibility to shallow landslides?

Forests can decrease the risk of shallow landslides by mechanically reinforcing the soil and positively influencing its water balance. However, little is known about the effect of different forest structures on slope stability. In the study area in St Antönien, Switzerland, we applied statistical prediction models and a physically‐based model for spatial distribution of root reinforcement in order to quantify the influence of forest structure on slope stability. We designed a generalized linear regression model and a random forest model including variables describing forest structure along with terrain parameters for a set of landslide and control points facing similar slope angle and tree coverage. The root distribution measured at regular distances from seven trees in the same study area was used to calibrate a root distribution model. The root reinforcement was calculated as a function of tree dimension and distance from tree with the root bundle model (RBMw). Based on the modelled values of root reinforcement, we introduced a proxy‐variable for root reinforcement of the nearest tree using a gamma distribution. The results of the statistical analysis show that variables related to forest structure significantly influence landslide susceptibility along with terrain parameters. Significant effects were found for gap length, the distance to the nearest trees and the proxy‐variable for root reinforcement of the nearest tree. Gaps longer than 20 m critically increased the susceptibility to landslides. Root reinforcement decreased with increasing distance from trees and is smaller in landslide plots compared to control plots. Furthermore, the influence of forest structure strongly depends on geomorphological and hydrological conditions. Our results enhance the quantitative knowledge about the influence of forest structure on root reinforcement and landslide susceptibility and support existing management recommendations for protection against gravitational natural hazards. Copyright © 2015 John Wiley & Sons, Ltd.     

A low‐dimensional physically based model of hydrologic control of shallow landsliding on complex hillslopes

Hillslopes have complex three‐dimensional shapes that are characterized by their plan shape, profile curvature of surface and bedrock, and soil depth. To investigate the stability of complex hillslopes (with different slope curvatures and plan shapes), we combine the hillslope‐storage Boussinesq (HSB) model with the infinite slope stability method. The HSB model is based on the continuity and Darcy equations expressed in terms of storage along the hillslope. Solutions of the HSB equation account explicitly for plan shape by introducing the hillslope width function and for profile curvature through the bedrock slope angle and the hillslope soil depth function. The presented model is composed of three parts: a topography model conceptualizing three‐dimensional soil mantled landscapes, a dynamic hydrology model for shallow subsurface flow and water table depth (HSB model) and an infinite slope stability method based on the Mohr–Coulomb failure law. The resulting hillslope‐storage Boussinesq stability model (HSB‐SM) is able to simulate rain‐induced shallow landsliding on hillslopes with non‐constant bedrock slope and non‐parallel plan shape. We apply the model to nine characteristic hillslope types with three different profile curvatures (concave, straight, convex) and three different plan shapes (convergent, parallel, divergent). In the presented model, the unsaturated storage has been calculated based on the unit head gradient assumption. To relax this assumption and to investigate the effect of neglecting the variations of unsaturated storage on the assessment of slope stability in the transient case, we also combine a coupled model of saturated and unsaturated storage and the infinite slope stability method. The results show that the variations of the unsaturated zone storage do not play a critical role in hillslope stability. Therefore, it can be concluded that the presented dynamic slope stability model (HSB‐SM) can be used safely for slope stability analysis on complex hillslopes. Our results show that after a certain period of rainfall the convergent hillslopes with concave and straight profiles become unstable more quickly than others, whilst divergent convex hillslopes remain stable (even after intense rainfall). In addition, the relation between subsurface flow and hillslope stability has been investigated. Our analyses show that the minimum safety factor (FS) occurs when the rate of subsurface flow is a maximum. In fact, by increasing the subsurface flow, stability decreases for all hillslope shapes. Copyright © 2008 John Wiley & Sons, Ltd.     

Bank undercutting and tension failure by groundwater seepage: predicting failure mechanisms

Groundwater seepage can lead to the erosion and failure of streambanks and hillslopes. Two groundwater instability mechanisms include (i) tension failure due to the seepage force exceeding the soil shear strength or (ii) undercutting by seepage erosion and eventual mass failure. Previous research on these mechanisms has been limited to non‐cohesive and low cohesion soils. This study utilized a constant‐head, seepage soil box packed with more cohesive (6% and 15% clay) sandy loam soils at prescribed bulk densities (1.30 to 1.70 Mg m^?3) and with a bank angle of 90° to investigate the controls on failure mechanisms due to seepage forces. A dimensionless seepage mechanism (SM) number was derived and evaluated based on the ratio of resistive cohesion forces to the driving forces leading to instability including seepage gradients with an assumed steady‐state seepage angle. Tension failures and undercutting were both observed dependent primarily on the saturated hydraulic conductivity, effective cohesion, and seepage gradient. Also, shapes of seepage undercuts for these more cohesive soils were wider and less deep compared to undercuts in sand and loamy sand soils. Direct shear tests were used to quantify the geotechnical properties of the soils packed at the various bulk densities. The SM number reasonably predicted the seepage failure mechanism (tension failure versus undercutting) based on the geotechnical properties and assumed steady‐state seepage gradients of the physical‐scale laboratory experiments, with some uncertainty due to measurement of geotechnical parameters, assumed seepage gradient direction, and the expected width of the failure block. It is hypothesized that the SM number can be used to evaluate seepage failure mechanisms when a streambank or hillslope experiences steady‐state seepage forces. When prevalent, seepage gradient forces should be considered when analyzing bank stability, and therefore should be incorporated into commonly used stability models. Copyright © 2013 John Wiley & Sons, Ltd.     

Modeling the spatial distribution of landslide‐prone colluvium and shallow groundwater on hillslopes of Seattle,WA

Landslides in partially saturated colluvium on Seattle, WA, hillslopes have resulted in property damage and human casualties. We developed statistical models of colluvium and shallow‐groundwater distributions to aid landslide hazard assessments. The models were developed using a geographic information system, digital geologic maps, digital topography, subsurface exploration results, the groundwater flow modeling software VS2DI and regression analyses. Input to the colluvium model includes slope, distance to a hillslope–crest escarpment, and escarpment slope and height. We developed different statistical relations for thickness of colluvium on four landforms. Groundwater model input includes colluvium basal slope and distance from the Fraser aquifer. This distance was used to estimate hydraulic conductivity based on the assumption that addition of finer‐grained material from down‐section would result in lower conductivity. Colluvial groundwater is perched so we estimated its saturated thickness. We used VS2DI to establish relations between saturated thickness and the hydraulic conductivity and basal slope of the colluvium. We developed different statistical relations for three groundwater flow regimes. All model results were validated using observational data that were excluded from calibration. Eighty percent of colluvium thickness predictions were within 25% of observed values and 88% of saturated thickness predictions were within 20% of observed values. The models are based on conditions common to many areas, so our method can provide accurate results for similar regions; relations in our statistical models require calibration for new regions. Our results suggest that Seattle landslides occur in native deposits and colluvium, ultimately in response to surface‐water erosion of hillslope toes. Regional groundwater conditions do not appear to strongly affect the general distribution of Seattle landslides; historical landslides were equally dispersed within and outside of the area potentially affected by regional groundwater conditions. Published in 2007 by John Wiley & Sons, Ltd.     

Landslide Amplification by Liquefaction of Runout‐Path Material after the 2008 Wenchuan (M 8·0) Earthquake,China

Here, we propose that an earthquake can trigger the failure of a landslide mass while simultaneously triggering liquefaction of runout‐path materials before the arrival of the landslide mass, thus greatly increasing the size and mobility of an overriding landslide. During the 2008 Wenchuan earthquake, about 60 000 landslides were triggered, directly resulting in about 20 000 casualties. While these landslides mainly originated from steep slopes, some landslides with high mobility formed in colluvial valley deposits. Among these, the most catastrophic was the Xiejiadian landslide in Pengzhou city, which traveled hundreds of meters before coming to rest. Through field investigation and laboratory testing, we conclude that this landslide primarily formed from colluvial deposits in the valley and secondarily from failure of slopes in granitic rock located uphill. Much of the granitic slope failure was deposited in the upper part of the travel path (near the slide head); the remainder was dispersed throughout the main landslide deposit. Superposition of deposits at the landslide toe indicates that landslide debris derived from colluvial soil was deposited first. The deposits at the landslide toe displayed flow characteristics, such as fine materials comprising basal layers and large boulders covering the deposit surface. We hypothesize that the main part of the landslide resulted from seismogenic liquefaction of valley colluvium, rather than from liquefaction potentially caused by undrained loading from the granitic slope failures impacting the colluvium. To examine the likelihood that seismogenic liquefaction occurred, we took samples from different areas of the landslide deposit and performed undrained cyclic shear tests on them in the laboratory. The results showed that the sandy soils that comprise most of the deposit are highly liquefiable under seismic loading. Therefore, we conclude that liquefaction of the colluvium in the valley during the earthquake was the main reason for this rapid (~46 m/s) long‐runout (1·7 km) landslide. Copyright © 2012 John Wiley & Sons, Ltd.     

A finite element model for simulating surface run‐off and unsaturated seepage flow in the shallow subsurface

This work introduces water–air two‐phase flow into integrated surface–subsurface flow by simulating rainfall infiltration and run‐off production on a soil slope with the finite element method. The numerical model is formulated by partial differential equations for hydrostatic shallow flow and water–air two‐phase flow in the shallow subsurface. Finite element computing formats and solution strategies are presented to obtain a numerical solution for the coupled model. An unsaturated seepage flow process is first simulated by water–air two‐phase flow under the atmospheric pressure boundary condition to obtain the rainfall infiltration rate. Then, the rainfall infiltration rate is used as an input parameter to solve the surface run‐off equations and determine the value of the surface run‐off depth. In the next iteration, the pressure boundary condition of unsaturated seepage flow is adjusted by the surface run‐off depth. The coupling process is achieved by updating the rainfall infiltration rate and surface run‐off depth sequentially until the convergence criteria are reached in a time step. A well‐conducted surface run‐off experiment and traditional surface–subsurface model are used to validate the new model. Comparisons with the traditional surface–subsurface model show that the initiation time of surface run‐off calculated by the proposed model is earlier and that the water depth is larger, thus providing values that are closer to the experimental results.     

Triggering Mechanisms of Slope Instability and their Relationship to Earthquakes and Tsunamis

— Submarine and shoreline slope failures that accompany large earthquakes and large tsunamis are triggered by several mechanisms. Triggering mechanisms range from direct effects, such as inertial forces from earthquake shaking, to indirect effects, such as rapid drawdown that occurs when an earthquake-generated tsunami first approaches a shoreline. Soil shear strength also plays an important role in earthquake-related slope failures. Earthquakes change the shear strength of the soil by inducing excess pore water pressures. These excess pore water pressures change with time after the earthquake, resulting in changes in shear strength and slope stability with time. This paper reviews earthquake-related triggering mechanisms for submarine and shoreline slope failures. The variation in shear strength with time following an earthquake is examined and it is shown that delayed slope failures after an earthquake can occur as a result of changes in earthquake-induced excess pore water pressures and shear strength with time.     

Prediction of rainfall‐induced shallow landslides in the Loess Plateau,Yan'an,China, using the TRIGRS model

In this work, a transient rainfall infiltration and grid‐based regional slope‐stability model (TRIGRS) was implemented in a case study of Yan'an City, Northwest China. In this area, widespread shallow landslides were triggered by the 12 July 2013 exceptional rainstorm event. A high‐resolution DEM, soil parameters from in‐situ and laboratory measurements, water table depths, the maximum depth of precipitation infiltration and rain‐gauge‐corrected precipitation of the event, were used as inputs in the TRIGRS model. Shallow landslides triggered on the same day were used to evaluate the modeling results. The summarized results are as follows: (i) The characteristics and distribution of thirty‐five shallow landslides triggered by the 12 July 2013 rainfall event were identified in the study area and all were classified as shallow landslides with the maximum depth, area and volume less than 3 m, 200 m² and 1000 m³, respectively, (ii) Four intermediate factor of safety (F_S) maps were generated using the TRIGRS model to represent the scenarios 6, 12, 18 and 24 hours after the storm event. The area with F_S < 1 increased with the rainfall duration. The percentage of the area with F_S < 1 was 0.2%, 3.3%, 3.8% and 5.1% for the four stages, respectively. Twenty‐four hours after the rainstorm, TRIGRS predicted that 1255 grid cells failed, which is consistent with the field data. (iii) TRIGRS generated more satisfactory results at a given precipitation threshold than SINMAP, which is ideal for landslide hazard zoning for land‐use planning at the regional scale. Comparison results showed that TRIGRS is more useful for landslide prediction for a certain precipitation threshold, also in the regional scale. (iv) Analysis of the responses of loess slope prone to slope failure after different precipitation scenarios revealed that loess slopes are particularly sensitive to extended periods of heavy precipitation. Copyright © 2016 John Wiley & Sons, Ltd.