Simulation of interactions among multiple debris flows

Adjacent debris flows may interact in many ways: two or more concurrent debris flows may merge; one debris flow can run out over an existing debris flow fan. Such interactions may cause debris flow properties to change in the mixing process as well as more severe adverse effects than those caused by a single debris flow. This paper aims to investigate the interactions among channelized debris flows originated from adjacent catchments. Both concurrent and successive debris flows are considered. If several debris flows originate from different locations concurrently and merge, the volumetric sediment concentration (i.e., the ratio of the volume of solid material to the total volume of debris flow), C _v, is a good index to capture the mixing process of these debris flows. The change in C _v reflects where mixing occurs and the mixing degree. The debris flow properties (e.g., yield stress and dynamic viscosity) evolve in the mixing process and can be captured by the change in C _v. The debris flow with a larger volume dominates the mixing process, and the properties of the mixed debris flow are more similar to those of the larger debris flow. The inundated areas and runout distances of successive debris flows are smaller than those of concurrent debris flows of the same total volume due to the smaller scales of the individual events and blockage by the earlier debris flows. However, the deposit depth in the interacting part of the debris flow fans of successive debris flows can be much larger than that of concurrent debris flows, leading to more destructive cascading hazards (e.g., the formation of debris barrier lakes). The sequence of successive debris flows not only significantly influences the runout characteristics of the debris flows but also substantially affects the cascading hazards.     

Mechanisms and runout characteristics of the rainfall-triggered debris flow in Xiaojiagou in Sichuan Province, China

The 2008 Wenchuan earthquake induced a large number of landslides, and a vast amount of loose landslide materials deposited on steep hill slopes or in channels. Such loose materials can become sources of deadly debris flows once triggered by storms. On 13 August 2010, a storm swept Yingxiu and its vicinity, triggering a catastrophic debris flow with a volume of 1.17?million?m³ in Xiaojiagou Ravine. The debris flow buried 1,100?m of road, blocked a river and formed a debris flow barrier lake. A detailed field study was conducted to understand the initiation mechanisms and runout characteristics of this debris flow. Two types of debris flows are identified, namely hill-slope debris flow and channelized debris flow. The hill-slope debris flow was triggered in the forms of firehose effect, rilling and landsliding, whereas the channelized debris flow was triggered in the form of channel-bed failure. This debris flow was a water?Crock flow since most particles were gravel, cobble or larger rocks and the fraction of silt and clay was less than 2%. Grain contact friction, pore-pressure effects and inertial grain collision were the three most important physical interactions within the debris flow. Such interactions yielded a smaller runout distance (593?m) compared with those of mud?Crock flows of similar size.     

Linking rainfall-induced landslides with predictions of debris flow runout distances

Rapid debris flows are among the most destructive natural hazards in steep mountainous terrains. Prediction of their path and impact hinges on knowledge of initiation location and the size and constitution of the released mass. To better link mass release initiation with debris flow paths and runout lengths, we propose to capitalize on a newly developed model for rainfall-induced landslide initiation (“Catchment-scale Hydro-mechanical Landslide Triggering” CHLT model, von Ruette et al. 2013) and couple it with simple estimates of debris flow runout distances and pathways. Landslide locations and volumes provided by the CHLT model are used as inputs to simulate debris flow runout distances with two empirical- and two physically-based models. The debris flow runout models were calibrated using two landslide inventories in the Swiss Alps obtained following a large rainfall event in 2005. We first fitted and tested the models for the “Prättigau” inventory, where detailed information on runout path was available, and then applied the models to landslides inventoried from a different catchment (“Napf”). The predicted debris flow runout distances (emanating from CHLT simulated landslide positions) were well in the range of observed values for the physically-based approaches. The empirical approaches tend to overestimate runout distances relative to observations. These preliminary results demonstrate the added value of linking shallow landslide triggering models with predictions of debris flow runout pathways for a range of soil states and triggering events, thus providing a more complete hazard assessment picture for debris flow exposure at the catchment scale.     

Initiation process of debris flows on different slopes due to surface flow and trigger-specific strategies for mitigating post-earthquake in old Beichuan County, China

The 12 May 2008 Wenchuan earthquake (Ms 8.0) in China, produced an estimated volume of 28 × 10⁸ m³ loosened material, which led to debris flows after the earthquake. Debris flows are the dominant mountain hazards, and serious threat to lives, properties, buildings, traffic, and post-earthquake reconstruction in the earthquake-hit areas. It is very important to understand the debris flow initiation processes and characteristics, for designing debris flow mitigation. The main objective of this article is to examine the different debris flow initiation processes in order to identify suitable mitigation strategies. Three types of debris flow initiation processes were identified (designated as Types A, B, and C) by field survey and experiments. In “A” type initiation, the debris flow forms as a result of dam failure in the process of rill erosion, slope failure, landslide dam, or dam failure. This type of debris flow occurs at the slope of 10 ± 2°, with a high bulk density, and several surges following dam failure. “B” type initiation is the result of a gradual increase in headward down cutting, bank and lateral erosion, and then large amount of loose material interfusion into water flow, which increases the bulk density, and forms the debris flow. This type of debris flow occurs mainly on slopes of 15 ± 3° without surges. “C” type debris flow results from slope failures by surface flow, infiltration, loose material crack, slope failure, and fluidization. This type of debris flow occurs mainly on slopes of 21 ± 4°, and has several surges of debris flow following slope failure, and a high bulk density. To minimize the hazards from debris flows in areas affected by the Wenchuan earthquake, the erosion control measures, such as the construction of grid dams, slope failure control measures, the construction of storage sediment dams, and the drainage measures, such as construction of drainage ditches are proposed. Based on our results, it is recommend that the control measures should be chosen based on the debris flow initiation type, which affects the peak discharge, bulk density and the discharge process. The mitigation strategies discussed in this paper are based on experimental simulations of the debris flows in the Weijia, Huashiban, and Xijia gullies of old Beichuan city. The results are useful for post-disaster reconstruction and recovery, as well as for preventing similar geohazards in the future.     

Design and performance of a novel multi-function debris flow mitigation system in Wenjia Gully,Sichuan

The post-earthquake debris flows in the Wenjia Gully led to the exposure of the shortcomings in the design of the original conventional debris flow mitigation system. A predicament for the Wenjia mitigation system is a large amount of loose material (est. 50 × 10⁶ m³) that has been deposited in the gully by the co-seismic landslide, providing abundant source material for debris flows under saturation. A novel design solution for the replacement mitigation system was proposed and constructed, and has exhibited excellent performance and resilience in subsequent debris flows. The design was governed by the three-phase philosophy of controlling water, sediment, and erosion. An Early Warning System (EWS) for debris flow that uses real-time field data was developed; it issues alerts based on the probabilistic and empirical correlations between rainfall and debris flows. This two-fold solution reduces energy of the debris flow by combining different mitigation measures while minimizing the impact through event forecasting and rapid public information sharing. Declines in the number and size of debris flows in the gully, with increased corresponding rainfall thresholds and mean rainfall intensity-duration (I-D) thresholds, indicate the high efficacy of the new mitigation system and a lowered debris flow susceptibility. This paper reports the design of the mitigation system and analyzes the characteristics of rainfall and debris flow events that occurred before and after implementation of the system; it evaluates the effectiveness of one of the most advanced debris flow mitigation systems in China.     

The formation and development of debris flows in large watersheds after the 2008 Wenchuan Earthquake

The Wenchuan earthquake has caused abundance of loose materials supplies for debris flows. Many debris flows have occurred in watersheds in area beyond 20 km², presenting characteristics differing from those in small watersheds. The debris flows yearly frequency decreases exponentially, and the average debris flow magnitude increases linearly with watershed size. The rainfall thresholds for debris flows in large watersheds were expressed as I?=?14.7 D ^?0.79 (2 h?<?D?<?56 h), which is considerably higher than those in small watersheds as I?=?4.4 D ^?0.70 (2 h?<?D?<?37 h). A case study is conducted in Ergou, 39.4 km² in area, to illustrate the formation and development processes of debris flows in large watersheds. A debris flow develops in a large watershed only when the rainfall was high enough to trigger the wide-spread failures and erosions on slope and realize the confluence in the watershed. The debris flow was supplied by the widely distributed failures dominated by rill erosions (14 in 22 sources in this case). The intermittent supplying increased the size and duration of debris flow. While the landslide dam failures provided most amounts for debris flows (57 % of the total amount), and amplified the discharge suddenly. During these processes, the debris flow velocity and density increased as well. The similar processes were observed in other large watersheds, indicating this case is representative.     

Numerical runout simulation of debris avalanches in the Faroe Islands, North Atlantic Ocean

The Faroe Islands in the North Atlantic Ocean are susceptible to flow-type landslides in coarse-grained highly organic colluvium. Following several hazardous debris avalanche events, research work has been initiated to quantify landslide risk. A central task in this work is to predict landslide runout behavior. From numerical simulation of four debris avalanches, this study provides a first screening of which rheology and appertaining input parameters best predict runout behavior of debris avalanches in the Faroe Islands. Three rheologies (frictional, Voellmy, and Bingham) are selected and used for individual back analysis of the events in the numerical models BING and DAN3D. A best fit rheology is selected from comparing predicted and observed landslide runout behavior. General back analysis to identify the optimal input parameters for the chosen rheology is performed by cross validation, where each debris avalanche is modeled with input parameters from the three other events. Optimal input parameters are found from the model run producing the most accurate runout length and velocity. The Bingham is selected as the best fit rheology, a result differing from similar studies of coarse-grained landslides. A reason for why particularly the frictional rheology proves unsuitable is its tendency to produce too long runout lengths of the low-weight runout material, a result showing important limitations for using the frictional rheology in DAN3D. Optimal Bingham input parameters are τ _y ?=?980 Pa and μ _b ?=?117 Pa/s. However, future studies performed in 2D models are needed for precise parameterization before results can be used for landslide risk assessment.     

Flow behaviour of the submarine glacigenic debris flows on the Bear Island Trough Mouth Fan, western Barents Sea

Using 3·5 kHz high-resolution seismic data, gravity cores and side-scan sonar imagery, the flow behaviour of submarine, glacigenic debris flows on the Bear Island Trough Mouth Fan, western Barents Sea was studied. During their downslope movement, the sediments within the uppermost part of the debris flows (<3 m) are inferred to have been deformed as a result of the shear stress at the debris–water interface. Thus, the uppermost part of the flow did not move downslope as a rigid plug. If present, a rigid part of the flow was located at least some metres below the surface. At c . 1000 to at least 1600 m water depth, the debris flows eroded and probably incorporated substrate debris. Further downslope, the debris flows moved passively over substrate sediments. The hypothesis of hydroplaning of the debris flow front may explain why the debris flows moved across the lower fan without affecting the underlying sediments. Detailed morphological information from the surface of one of the debris flow deposits reveals arcuate ridges. These features were probably formed by flow surge. Hydroplaning of the debris flow front may also explain the formation of flow surge. The long runout distance of some of the large debris flows could be due to accretion of material to the base of the debris flow, thereby increasing in volume during flow, and/or to hydroplaning suppressing deceleration of the flow.     

A method to develop the input parameter database for site-specific debris flow hazard prediction under extreme rainfall

The predictive hazard analysis at a detailed scale for debris flow runout analysis can be improved significantly through reliable estimation of the input parameters. In this study, a method for database establishment of input parameters at a site-specific scale was laid out for the predictive-based debris flow hazard assessment under extreme rainfall. The adoption of the DAN-3D code necessitated the estimation of three main input parameters: initial volume, bulk basal frictional angle, and growth rate. The initial volume was assessed using a 3D coupled finite element seepage and limit equilibrium-based slope stability analysis. An artificial neural network-based model was developed using 27 debris flow events for predicting the basal bulk frictional angle and consisted of eight factors: plan curvature, profile curvature, percentage of fine content, D₅₀, initial unit weight, initial volume, relative relief ratio, and channel length. Finally, the growth rate was estimated using the previously assessed initial volume, soil depth, and the approximate runout length. The proposed method was validated by application to the Raemian slope in the Woomyeon mountain region, Seoul, for the extreme rainfall event of 27 July 2011. The analysis yielded a final volume of 53,067.9 m³, a velocity upon arrival on the road of 26.81 m/s, and an approximately 0.5-m debris thickness concentrated near the Raemian apartments. The comparison of the predicted debris flow path and debris flow velocity with the actual event demonstrates good similarity and provides a conservative estimate of the volume. This study therefore illustrates the importance of an input parameter database in providing a reliable debris flow runout hazard assessment.     

泥石流的二维数学模型

泥石流是在重力作用下,由砂粒石块和水等组成的固液混合物,是一种发生于山区的复杂的地质灾害现象。泥石流主要是由暴雨诱发引起的,它沿着复杂的三维地形高速流动,具有流体流动的特性。为了模拟泥石流的运动规律,预测降雨诱发的泥石流的到达距离和泛滥范围,减少和避免泥石流引起的灾害,把泥石和雨水组成的固液混合物假定为遵循均匀、连续、不可压缩的、非定常的牛顿流体运动规律。基于质量守恒方程和Naiver-stokes方程,采用深度积分方法,推导出了一个模拟泥石流运动的二维数学模型。所有方程式可用有限差分法来求解。结合GIS,该模型可用于预测泥石流的流动距离和泛滥范围,以及泛滥范围内的危险房屋和路段,也可以用于泥石流灾害的风险性分析。     

The Saharan debris flow: an insight into the mechanics of long runout submarine debris flows

New 3·5 kHz profiles and a series of piston cores from the north-west African margin provide evidence that the Saharan debris flow travelled for more than 400 km on a highly fluid, low-friction layer of poorly sorted sediment. Data suggest that the Saharan debris flow is a two-phase event, consisting of a basal, volcaniclastic debris flow phase overlain by a pelagic debris flow phase. Both phases were emplaced on the lower continental rise by a single large debris flow at around 60 ka. The volcaniclastic flow left a thin deposit less than 5 m thick. This contrasts with the much thicker (over 25 m) deposit left by the pelagic debris flow phase. We suggest that pelagic sediment, sourced and mobilized as debris flow from the African continental margin, loaded and destabilized volcaniclastic material in the vicinity of the western Canaries. When subjected to this loading, the volcaniclastic material appears to have formed a highly fluid sandy debris flow, capable of transporting with it the huge volumes of pelagic debris, and contributing to a runout distance extending over 400 km downslope of the Canary Islands on slopes that decrease to as little as 0·05°. It is likely that the pelagic debris formed a thick impermeable slab above the volcanic debris, thus maintaining high pore pressures generated by loading and giving rise to low apparent friction conditions. The distribution of the two debris phases indicates that the volcaniclastic debris flow stopped within a few tens of kilometres after escaping from beneath the pelagic debris flow, probably because of dissipation of excess pore pressure when the seal of pelagic material was removed.     

Catastrophic debris flows triggered by a 4 July 2013 rainfall in Shimian,SW China: formation mechanism,disaster characteristics and the lessons learned

On 4 July 2013, three catastrophic debris flows occurred in the Hougou, Majingzi, and Xiongjia gullies in Shimian county and produced debris dams and river blockages, resulting in serious casualties and huge economic loss. Though debris flows have been identified prior to the catastrophic events, their magnitudes and destructive power were far beyond early recognition and hazard assessment. Our primary objective for this study was to explore the formation mechanism and typical characteristics and to summarize the lessons learned from these disastrous events in order to avoid the repeat of such disasters in the future. Based on field investigation and imagery interpretation of remote sensing carried out following the catastrophic events, four conclusions were drawn: (1) The catastrophic debris flows were initiated from surface-water runoff, and the triggering factor was attributed to the local intensive rainfall with an hourly intensity of more than 46.7 mm. (2) Entrainment was the most important sediment-supplying method for the debris flow occurrence, and the source materials transported by debris flows from the three gullies were estimated to be about 97?×?10⁴ m³ in volume altogether. (3) As surface-water runoff eroded and entrained hillslope and channel materials persistently, debris flows were characterized by intensive incision at upper or middle reaches and significant magnification effect in flow discharge and volume downstream. Corresponding peak discharge surveyed at the outlets of the Hougou, Majingzi, and Xiongjia gullies was estimated up to 751.0 m³/s, 870.1 m³/s, and 758.7 m³/s, respectively. (4) Debris flows that occurred from the three gullies all belonged to viscous ones and the bulk densities were calculated more than 1.80 g/cm³, indicating a huge carrying capacity and destructive impacting power. In addition, the lessons learned from the catastrophic events were summarized, including recognition and assessment on debris flow hazard and utilization of deposition fan. In this paper, prevention suggestions on debris flow prone valleys with high-vegetation coverage and low occurrence frequency were also put forward. The results of this study contribute to a better understanding on the initiation mechanism, dynamic characteristics, and disaster mitigation of debris flows initiated from intense rainfall and surface-water runoff in mountainous areas.     

四川泸定MS6.8级地震区湾东河流域泥石流活动性预测

2022年9月5日四川泸定县发生M_S 6.8级地震,地震诱发大量同震崩滑体,并导致湾东河断流。基于现场调查、影像解译和区域地质资料分析,采用空间统计和水文计算的方法,对湾东河流域同震崩滑体分布特征和潜在泥石流危险性进行了研究。结果表明：湾东河流域内同震崩滑体主要分布在地震烈度Ⅸ度区,规模以中小型为主,主要沿沟道两侧展布,尤其是单薄山脊两侧临空面发育密度较大,距断层距离和坡度对其分布具有明显的控灾效应;未来湾东河流域暴发溃决型泥石流的冲出量可能为同等触发条件下震前泥石流的约两倍。依此提出了加强流域内溃决型泥石流风险防范,尽快通过综合监测预警获取泥石流发生的临界雨量值,在泥石流防治工程设计中应充分考虑泥石流规模放大系数等防灾减灾建议,为泸定地震后泥石流灾害防灾减灾提供科学参考。     

Driftwood deposition from debris flows at slit-check dams and fans

Experience shows that debris flows containing large woody debris (driftwood) can be more damaging than debris flows without driftwood. In this study, the deposition process of debris flows carrying driftwood was investigated using numerical simulations and flume experiments. Debris-flow trapping due to driftwood jamming in a slit-check dam was also investigated. A numerical model was developed with an interacting combination of Eulerian expression of the debris flow and Lagrangian expression of the driftwood, in which the fluctuating coordinates and rotation of the driftwood were treated stochastically. The calculated shapes and thicknesses of a debris-flow fan and the positions and orientations of the deposited driftwood on a debris-flow fan were consistent with experimental flume results. The jamming of driftwood in a slit-check dam was evaluated based on geometry and probability. The simulated results of outflow discharge and the proportion of driftwood passed through the slit-check dam also agreed with the experimental results.     

A hypoplastic constitutive model for debris materials

Debris flow is a very common and destructive natural hazard in mountainous regions. Pore water pressure is the major triggering factor in the initiation of debris flow. Excessive pore water pressure is also observed during the runout and deposition of debris flow. Debris materials are normally treated as solid particle–viscous fluid mixture in the constitutive modeling. A suitable constitutive model which can capture the solid-like and fluid-like behavior of solid–fluid mixture should have the capability to describe the developing of pore water pressure (or effective stresses) in the initiation stage and determine the residual effective stresses exactly. In this paper, a constitutive model of debris materials is developed based on a framework where a static portion for the frictional behavior and a dynamic portion for the viscous behavior are combined. The frictional behavior is described by a hypoplastic model with critical state for granular materials. The model performance is demonstrated by simulating undrained simple shear tests of saturated sand, which are particularly relevant for the initiation of debris flows. The partial and full liquefaction of saturated granular material under undrained condition is reproduced by the hypoplastic model. The viscous behavior is described by the tensor form of a modified Bagnold’s theory for solid–fluid suspension, in which the drag force of the interstitial fluid and the particle collisions are considered. The complete model by combining the static and dynamic parts is used to simulate two annular shear tests. The predicted residual strength in the quasi-static stage combined with the stresses in the flowing stage agrees well with the experimental data. The non-quadratic dependence between the stresses and the shear rate in the slow shear stage for the relatively dense specimens is captured.     

Hazard assessment of debris flows in the Wenchuan earthquake-stricken area,South West China

The Longxi river basin with the city of Dujiangyan, in the Sichuan province of South West China, belongs to the seismic area of the May 12, 2008 Wenchuan earthquake. Lots of loose co-seismic materials were present on the slopes, which in later years served as source material for rainfall-induced debris flows. A total of 12 debris flows, were triggered by heavy rainfall on August 13, 2010 in the study area. The FLO-2D numerical analysis software was adopted to simulate debris flows intensity, including movement velocities and maximum flow depths. A comparison of the measured fan spreading with the simulation results, the evaluation parameter Ω was used to verify accuracy of simulation, the results show Ω values ranging between 1.37 and 1.65 indicating relative good simulation results. This study also estimated the flood hydrograph for various recurrence intervals (20, 100, and 200 years, respectively) to perform scenario simulations of debris flows, and followed Swiss and Austrian standards to establish a debris flow hazard classification model on the basis of a combination of the debris flow intensity and the recurrence period. This study distinguishes three hazard classes: low, medium, and high. This proposed approach generated a debris flow hazard distribution map that could be used for disaster prevention in the Wenchuan earthquake-stricken area, South West China.     

Effects of the morphology of sediment-transporting channels on the erosion and deposition of debris flows

A series experiments are conducted to investigate the effects of streambed profile on the erosion and deposition of debris flows. It is found that straight channel can increase the run out of debris flows by 10–25%, compared to that of surfaces without channels, and that travel distance was positively correlated with the hydraulic radius of the channel. In addition, the presence of straight channels caused the volume of debris flow deposition to become normally distributed with respect to travel distance. In the case of curved channels, increases in the sinuosity index resulted in significant blockage and obstruction. In the deposition zone, the maximum deposition volume for a channel with a comparatively low sinuosity index (1.05) was <?50% of the minimum deposition volume for a straight channel. Furthermore, the channel curvature affected not only the positions of deposition peaks along the travel distance but also the debris flow magnitudes in each unit interval (0.5 m). This study demonstrates the effects of differences in channel morphology on the erosional and depositional processes of gully debris flows. These findings are of significant importance for guiding debris flow risk assessment and for the restoration and reconstruction of downstream regions.     

四川省九龙县娃娃沟泥石流运动特征与堵河成灾机制分析

2006年7月16日娃娃沟流域暴发的大规模泥石流,给下游3个电站造成巨大经济损失,是大渡河流域一次典型的灾害性泥石流。分析得出,娃娃沟泥石流重度高、搬运能力强,泥石流固体物质砂、石混杂,粗大砾石含量高;暴发频率低、规模大,流速及峰值流量分别高达10．78m／s及798．5m^3／s;在汇口处,泥石流堆积物堵塞河道是引起下游电站受灾的重要原因,高重度、粗颗粒、大流量的组合是此次泥石流堵江的重要原因。堵河判别计算结果显示在发生百年一遇泥石流时,该断面均有发生堵河的可能。娃娃沟泥石流表明：①在大渡河支流的泥石流沟周边的中小电站极有可能在泥石流暴发时受到破坏。因此,电站建设过程中应加强对周边泥石流沟的防灾减灾工作;②虽然娃娃沟流域植被良好,但仍然发生了大规模泥石流。表明植被不能完全避免泥石流的发生,对于此类泥石流沟不能疏忽大意。     

Analysis of the triggering conditions and erosion of a runoff-triggered debris flow in Miyun County,Beijing, China

The occurrence of debris flow from channel-bed failure is occasionally noted in small and steeply sloping watersheds where channelized water flow dominates debris flow initiation. On August 12, 2016, a debris flow from channel-bed failures occurred in the Caozhuangzi Watershed of the Longtan Basin, Miyun, Beijing. Rainfall records over 10-min intervals and field investigations including channel morphology measurements were used to study the triggering conditions and erosion process. The results indicated that the occurrence of this event lagged the peak 10-min rainfall interval and that the cumulative rainfall prior to the occurrence time played an important role in its formation. A mean 10-min rainfall intensity–duration expression in the form of I₁₀?=?5.0?×?D^?0.21, where I₁₀ denotes the mean 10-min rainfall intensity and D is the rainfall duration ranging from 10 to 60 h, was proposed. The debris flows have low proportions of grain size fractions <?0.1 mm and higher fractions of grains 0.1–2 mm in size, indicating that the flow had low viscosity and was coarse-grain dominated. Channel morphology analysis revealed that abrupt changes in topography in the study area, including a steep section, a concave stream bank area, and a partial concave stream section were eroded more extensively than other sites. The maximum sediment erosion volume and erosion depth were not proportional to the variation in stream gradient. Consideration of the degree of erosion in the channel at sites with abrupt morphology changes, the maximum sediment erosion volume, and the erosion depth and volume at the initial channel site and downstream region of forest area together showed that the prime factor controlling erosion was entrained sediment volume. This work, thus, provides a case study regarding the triggering conditions of runoff-triggered debris flows and the topographical changes by debris flow erosion.     

Characteristics and numerical runout modeling of the heavy rainfall-induced catastrophic landslide–debris flow at Sanxicun,Dujiangyan, China,following the Wenchuan Ms 8.0 earthquake

The 2008 Ms 8.0 Wenchuan earthquake triggered a large number of extensive landslides. It also affected geologic properties of the mountains such that large-scale landslides followed the earthquake, resulting in the formation of a disaster chain. On 10 July 2013, a catastrophic landslide–debris flow suddenly occurred in the Dujiangyan area of Sichuan Province in southeast China. This caused the deaths of 166 people and the burying or damage of 11 buildings along the runout path. The landslide involved the failure of ≈1.47 million m³, and the displaced material from the source area was ≈0.3 million m³. This landslide displayed shear failure at a high level under the effects of a rainstorm, which impacted and scraped an accumulated layer underneath and a heavily weathered rock layer during the release of potential and kinetic energies. The landslide body entrained a large volume of surface residual diluvial soil, and then moved downstream along a gully to produce a debris flow disaster. This was determined to be a typical landslide–debris flow disaster type. The runout of displaced material had a horizontal extent of 1200 m and a vertical extent of 400 m. This was equivalent to the angle of reach (fahrböschung angle) of 19° and covered an area of 0.2 km². The background and motion of the landslide are described in this study. On the basis of the above analysis, dynamic simulation software (DAN3D) and rheological models were used to simulate the runout behavior of the displaced landslide materials in order to provide information for the hazard zonation of similar types of potential landslide–debris flows in southeast China following the Wenchuan earthquake. The simulation results of the Sanxicun landslide revealed that the frictional model had the best performance for the source area, while the Voellmy model was most suitable for the scraping and accumulation areas. The simulations estimated that the motion could last for ≈70 s, with a maximum speed of 47.7 m/s.