Debris Flow Formation Process and Critical Hydrodynamic Conditions in the Meizoseismal Area of the Wenchuan Earthquake

The Wenchuan earthquake generated strong surface disturbances and triggered a large number of loose deposits, resulting in the disasterprone environment with special watershed hydrological characteristics. This paper was to propose a debris flow formation process and explore the permeability characteristics and critical hydrodynamic conditions of the loose deposits triggered by the earthquake. The Guo Juanyan gully(31°05′27″ N to 31°05′46″ N, 103°36′58″ E to 103°37′09″ E) in Du Jiangyan City, located in the meizoseismal areas of the Wenchuan earthquake, was chosen as the study area and the disaster-prone environment was analyzed. The formation process of the debris flow was first proposed using a stability analysis, and then, the permeability characteristics of loose deposits were determined via in situ permeability experiments. Finally, the critical 1 h rainfall was simulated through a distributed hydrological model and verified by field observations. The formation process of debris flow could be divided into three stages based on the relationship between the hydrodynamic force and loose deposit resistance. The critical 1 h rainfall amounts under three antecedent moisture conditions(I-dry, II-normal and III-wet) were 52 mm/h, 43 mm/h and 34 mm/h, respectively. This study proposed a debris flow formation process in the meizoseismal areas of the Wenchuan earthquake based on the stability analysis and defined the rainfall threshold for debris flow early warning at the local level, which is significant for debris flow mitigation and risk management.     

Influences of the Wenchuan Earthquake on sediment supply of debris flows

The 5.12 Wenchuan Earthquake and the subsequent rainstorms induced a large number of landslides, which later were transformed into debris flows. To evaluate the effect of the earthquake on the sediment supply of debris flows, eight debris flow basins near Beichuan City, Sichuan Province, China were chosen as the study area. The area variations of the debris flow source after the Wenchuan Earthquake and the subsequent rainstorm are analyzed and discussed in this paper. Interpretations of aerial photographs (after the 5.12 Wenchuan Earthquake) and SPOT5 images (after the rainstorm event of September 24, 2008) as well as field investigations were compared to identify the transformation of landslide surface in the study area, indicating that the landslide area in the eight debris flow basins significantly increased. The loose sediment area on the channel bed increased after the rainstorm event. In order to estimate the relationship of the landslide area with the rainfall intensity in different return periods, a model proposed by Uchihugi was adopted. Results show that new landslide area induced by heavy rainfall with 50-year and 100-year return period will be 0.87 km2 and 1.67 km2, respectively. The study results show the Wenchuan earthquake had particular influences on subsequent rainfall-induced debris flow occurrence.     

Debris flow susceptibility assessment using a probabilistic approach: A case study in the Longchi area,Sichuan province,China

The Longchi area with the city of Dujiangyan, in the Sichuan province of China, is composed of Permian stone and diorites and Triassic sandstones and mudstones intercalated with slates. An abundance of loose co-seismic materials were present on the slopes after the May 12, 2008 Wenchuan earthquake, which in later years served as source material for rainfall-induced debris flows or shallow landslides. A total of 48 debris flows, all triggered by heavy rainfall on 13th August 20l0, are described in this paper. Field investigation, supported by remote sensing image interpretation, was conducted to interpret the co-seismic landslides in the debris flow gullies. Specific characteristics of the study area such as slope, aspect, elevation, channel gradient, lithology, and gully density were selected for the evaluation of debris flow susceptibility. A score was given to all the debris flow gullies based on the probability of debris flow occurrence for the selected factors. In order to get the contribution of the different factors, principal component analyses were applied. A comprehensive score was obtained for the 48 debris flow gullies which enabled us to make a susceptibility map for debris flows with three classes. Twenty-two gullies have a high susceptibility, twenty gullies show a moderate susceptibility and six gullies have a low susceptibility for debris flows.     

Damming of large river by debris flow: Dynamic process and particle composition

The frequency and extent of debris flows have increased tremendously due to the extreme weather and the Wenchuan earthquake on May 12, 2008. Previous studies focused on the debris flow from gullies damming the mountain streams. In this paper, an equation for the run-out distance of debris flow in the main river is proposed based on the dynamic equation of debris flow at different slopes given by Takahashi. By undertaking field investigations and flume experiments, a new calculation method of the volume of debris flow damming large river is obtained. Using the percolation theory and the renormalization group theory it was deduced that the large particles should comprise more than 50% for forming a stable debris flow dam. Hence, the criteria of damming large river by debris flow is presented in terms of run-out distance and grain composition which was then validated through the event of damming river by debris flow at Gaojia gully, the upper reaches of the Minjiang River, Sichuan, China, on July 3, 2011.     

Curvature derived from LiDAR digital elevation models as simple indicators of debris-flow susceptibility

To mitigate the damage caused by debris flows resulting from heavy precipitation and to aid in evacuation plan preparation, areas at risk should be mapped on a scale appropriate for affected individuals and communities. We tested the effectiveness of simply identifying debris-flow hazards through automated derivation of surface curvatures using LiDAR digital elevation models. We achieved useful correspondence between plan curvatures and areas of existing debris-flow damage in two localities in Japan using the analysis of digital elevation models(DEMs). We found that plan curvatures derived from 10 m DEMs may be useful to indicate areas that are susceptible to debris flow in mountainous areas. In residential areas located on gentle sloping debris flow fans, the greatest damage to houses was found to be located in the elongated depressions that are connected to mountain stream valleys. Plan curvaturederived from 5 m DEM was the most sensitive indicators for susceptibility to debris flows.     

Changes in runout distances of debris flows over time in the Wenchuan earthquake zone

A large number of debris flows occurred in the Wenchuan earthquake zone after the 12 May 2008 earthquake.The risks posed by these debris flows were rather high.An appropriate model is required to predict the possible runout distance and impacted area.This paper describes a study on the runout characteristics of the debris flows that occurred in the Wenchuan earthquake zone over the past four years.A total of 120 debris flows are analyzed.Separate multivariate regression models are established for the runout distances of hill-slope debris flows and channelized debris flows.The control variables include type of debris flow,debris flow volume,and elevation difference.Comparison of the debris flows occurring before and after the earthquake shows that the runout distance increased after the earthquake due to sufficient material supply and increased mobility of the source materials.In addition,the runout distances of annual debris flow events in 2008,2010 and 2011 are analyzed and compared.There is a tendency that the runout distance decreases over time due to the decreasing source material volume and possible changes of debris flow type.Comparison between the debris flows in the earthquake zone and the debris flows in Swiss Alps,Canada,Austria,and Japan shows that the former have a smaller mobility.     

Experimental analysis for the dynamic initiation mechanism of debris flows

Debris flow is one of the major secondary mountain hazards following the earthquake. This study explores the dynamic initiation mechanism of debris flows based on the strength reduction of soils through static and dynamic triaxial tests. A series of static and dynamic triaxial tests were conducted on samples in the lab. The samples were prepared according to different grain size distribution, degree of saturation and earthquake magnitudes. The relations of dynamic shear strength, degree of saturation, and number of cycles are summarized through analyzing experimental results. The findings show that the gravelly soil with a wide and continuous gradation has a critical degree of saturation of approximately 87%, above which debris flows will be triggered by rainfall, while the debris flow will be triggered at a critical degree of saturation of about 73% under the effect of rainfall and earthquake(M>6.5). Debris flow initiation is developed in the humidification process, and the earthquake provides energy for triggering debris flows. Debris flows are more likely to be triggered at the relatively low saturation under dynamic loading than under static loading. The resistance of debris flow triggering relies more on internal frication angle than soil cohesion under the effect of rainfall and earthquake. The conclusions provide an experimental analysis method for dynamic initiation mechanism of debris flows.     

Prediction of debris-flow danger area by combining hydrological and inundation simulation methods

Debris flows have caused serious human casualties and economic losses in the regions strongly affected by the Ms8.0 Wenchuan earthquake of 2008. Debris flow mitigation and risk assessment is a key issue for reconstruction.The existing methods of inundation simulation are based on historical disasters and have no power of prediction.The rain-flood method can not yield detailed flow hydrograph and does not meet the need of inundation simulation. In this paper,the process of water flow was studied by using the Arc-SCS model combined with hydraulic method,and then the debris flow runoff process was calculated using the empirical formula combining the result from Arc-SCS.The peak discharge and runoff duration served as input of inundation simulation. Then,the dangerous area is predicted using kinematic wave method and Manning equation.Taking the debris flow in Huashiban gully in Beichuan County,Sichuan Province,China on 24 Sep.2008 as example,the peak discharge of water flow and debris flow were calculated as 35.52 m3·s-1 and 215.66 m3·s-,with error of 4.15%compared to the measured values.The simulated area of debris-flow deposition was 161,500 m2,vs.the measured area of 144,097 m2,in error of 81.75%.The simulated maximum depth was 12.3 m,consistent with the real maximum depth between 10 and 15 m according to the field survey.The minor error is mainly due to the flow impact on buildings and variations in cross-section configuration.The present methodology can be applied to predict debris flow magnitude and evaluate its risk in other watersheds inthe earthquake area.     

Types and Causes of Debris Flow Damage to Drainage Channels in the Wenchuan Earthquake Area

Debris flows are among the most common geological disasters in China,and have been particularly frequent in Sichuan Province since the Wenchuan earthquake on 12 May 2008.The construction of debris flow drainage channels is a countermeasure used to distribute debris flow fans,and these channels play a critical role in the mitigation and prevention of damage resulting from debris flows.Under field conditions,the useful life of drainage channels can be greatly shortened as a result of strong abrasions to the drainage structure caused by the debris flow.Field investigations have shown that the types of damage to drainage channels include（a） erosion caused by hyper-concentrated silt flow,（b） impact fractures and foundation scour at the groundsills of the drainage channel,（c） destruction of the drainage channel outlet,and（d） destruction of the drainage channel caused by debris flow abrasion.In addition,based on the destruction of the drainage channel during the debris flow drainage process,a new type of drainage channel with energy dissipation components was proposed and applied in a steep,narrow gully for debris flow mitigation.Moreover,design and engineering repair recommendations for drainage channels are provided as a reference for repairing the damage to the channel.The results can provide an important reference for the effective repair and optimal design of drainage channels.     

Risk assessments of debris flow based on improved analytic hierarchy process and efficacy coefficient method

In order to evaluate the danger of debris flow properly, eight factors were selected as the risk assessment indexes of the debris flow, namely the vertical slope, valley relative difference, hillside slope, area of basin, loose solid material reserves, the path length of sediment supply probability, silting and scouring derricking and vegetation coverage. The improved Analytic Hierarchy Process (AHP) method was used to obtain the weights of the factors; and the efficacy coefficient method was adopted to evaluate the risks of six typical debris flow gullies. According to the research, the improved AHP method not only avoids the subjectivity in the individual factor valuation by comparing two factors of each layer, but also makes the subsequent consistency check unnecessary.     

Effectiveness analysis of the prediction of regional debris flow susceptibility in post-earthquake and drought site

Debris-flow disasters occurred frequently after the Mw 8.0 Wenchuan earthquake on 12 May 2008 in Sichuan Province, China. Based on historical accounts of debris-flow disaster events, it found that debris flow occurrence is closely related to the impact of earthquakes and droughts, because earthquakedrought activities can increase the loose solid materials, which can transform into debris flows under the effect of rainstorms. Based on the analysis of historical earthquake activity(frequency, magnitude and location), drought indexes and the trend of climate change(amount of rainfall), a prediction method was established, and the regional debris flow susceptibility was predicted. Furthermore, in a debris flow-susceptible site, effective warning and monitoring are essential not only from an economicpoint of view but are also considered as a frontline approach to alleviate hazards. The advantages of the prediction and early monitoring include(1) the acquired results being sent to the central government for policy making;(2) lives and property in mountainous areas can be protected, such as the 570 residents in the Aizi valley, who evacuated successfully before debris flows in 2012;(3) guiding the government to identify the areas of disasters and the preparation for disaster prevention and mitigation, such as predicting disasters in high-risk areas in the period 2012-2017, helping the government to recognize the development trend of disasters;(4) the quantitative prediction of regional debris-flow susceptibility, such as after the Wenchuan earthquake, can promote scientific and sustainable development and socioeconomic planning in earthquake-struck areas.     

Impacts of future climate change (2030-2059) on debris flow hazard: A case study in the Upper Minjiang River basin,China

An increase in extreme precipitation events due to future climate change will have a decisive influence on the formation of debris flows in earthquake-stricken areas. This paper aimed to describe the possible impacts of future climate change on debris flow hazards in the Upper Minjiang River basin in Northwest Sichuan of China, which was severely affected by the 2008 Wenchuan earthquake. The study area was divided into 1285 catchments, which were used as the basic assessment units for debris flow hazards. Based on the current understanding of the causes of debris flows, a binary logistic regression model was used to screen key factors based on local geologic, geomorphologic, soil, vegetation, and meteorological and climatic conditions. We used the weighted summation method to obtain a composite index for debris flow hazards, based on two weight allocation methods: Relative Degree Analysis and rough set theory. Our results showed that the assessment model using the rough set theory resulted in better accuracy. According to the bias corrected and downscaled daily climate model data, future annual precipitation (2030-2059) in the study area are expected to decrease, with an increasing number of heavy rainfall events. Under future climate change, areas with a high-level of debris flow hazard will be even more dangerous, and 5.9% more of the study area was categorized as having a high-level hazard. Future climate change will cause an increase in debris flow hazard levels for 128 catchments, accounting for 10.5% of the total area. In the coming few decades, attention should be paid not only to traditional areas with high-level of debris flow hazards, but also to those areas with an increased hazard level to improve their resilience to debris flow disasters.     

Dynamics of loose granular flow and its subsequent deposition in a narrow mountainous river

A large amount of loose debris materials were deposited on the slope of mountainous areas after the 2008 Ms 8.0 Wenchuan earthquake. During and after the earthquake, these loose debris deposits collapsed and slide into valleys or rivers, changing river sediment supply condition and channel morphology. To investigate the mechanisms of granular flow and deposition, the dynamics of slope failure and sediment transportation in typical mountainous rivers of different intersection angles were analyzed with a coupling model of Computational Fluid Dynamics and Discrete Element Method(CFD-DEM). The numerical results show that the change of intersection angle between the granular flow flume and the river channel can affect the deposit geometry and the fluid flow field significantly. As the intersection angle increases, the granular velocity perpendicular to the river channel increases, while the granular velocity parallel to the river channel decreases gradually. Compared to the test of dry granular flow, the CFD-DEM coupling tests show much higher granular velocity and larger volume of sediments entrained in the river. Due to the river flow, particles located at the edge of the deposition will move downstream gradually and the main section of sediments deposition moves from the center to the edge of the river channel. As a result, sediment supply in the downstream river will distribute unevenly. Under the erosion of fluid flow, the proportion of fine particles increases, while the proportion of coarse particles decreases gradually in the sediment deposition. The build-up of accumulated sediment mass will cause a significant increase in water level in the river channel, thus creating serious flooding hazard in mountainous rivers.     

The critical rainfall characteristics for torrents and debris flows in the Wenchuan earthquake stricken area

Critical rainfall assessment is a very important tool for hazard management of torrents and debris flows in mountainous areas. The Wenchuan Earthquake 2008 caused huge casualties and property damages in the earthquake-stricken area, which also generated large quantities of loose solid materials and increased occurrence probabilities of debris flows. There is an urgent need to quantify the critical rainfall distribution in the area so that better hazard management could be planned and if real time rainfall forecast is available, torrent and debris flow early-warning could be issued in advance. This study is based on 49-year observations （1954-2003） of up to 678 torrent and debris flow events. Detailed contour maps of 1 hour and 24 hour critical rainfalls have been generated （Due to the data limitation, there was insufficient 10 minute critical rainfall to make its contour map）. Generally, the contour maps from 1 hour and 24 hours have similar patterns. Three zones with low, medium and high critical rainfalls have been identified. The characteristics of the critical rainfall zones are linked with the local vegetation cover and land forms. Further studies and observations are needed to validate the finding and improve the contour maps.     

Activity evolution of landslides and debris flows after the Wenchuan earthquake in the Qipan catchment,Southwest China

The Wenchuan earthquake that occurred on 12 May 2008 induced numerous landslides. Loose landslide materials were deposited on hillslopes, and deep channels were easily remobilized and transformed into debris flows by extreme rainstorms. Twelve years after the Wenchuan earthquake, debris flows were still active in the Qipangou Ravine in the quake-hit area. In this paper, we continuously tracked the spatiotemporal evolution of the landslides and vegetation restoration and evaluated the evolution of debris flow activity in the Qipan catchment with the aid of a GIS platform and field investigations from 2008 to 2019. We observed that the area with active landslides increased sharply immediately following the earthquake, and then decreased with time; however, the total area of landslides continued to increase from 6.93 km2 in 2008 to 10.55 km2 in 2019. The active landslides shifted towards lower angles and higher elevations after 2013. Since 2009, the vegetation coverage has been gradually increasing and approaching the coverage present before the earthquake as of 2019. The landslide activity was high and the vegetation recovery rates were rapidly rising during the first five years after the earthquake; the recovery rates then slowed over time. Therefore, we divided the evolution that occurred during the post landslide period into an active period(2008-2013), a self-adjustment period(2013-2026) and a stable period(after 2026). We then proposed a quantitative model to determine the trends of landslide activity rates and NDVI values in the catchment, which indicated that the landslide activities and postseismic vegetation restoration rates in this catchment will return to preseismic levels within approximately two decades. We also analysed the runout volumes of the debris flows after the earthquakes(Diexi and Wenchuan) and the standard deviation of the vegetation coverage and predicted that the debris flow activities will last for an additional 50 years or more.     

Effects of loose deposits on debris flow processes in the Aizi Valley,southwest China

Loose deposits, rainfall and topography are three key factors that triggering debris flows.However, few studies have investigated the effects of loose deposits on the whole debris flow process.On June 28, 2012, a catastrophic debris flow occurred in the Aizi Valley, resulting in 40 deaths.The Aizi Valley is located in the Lower Jinsha River,southwestern Sichuan Province, China. The Aizi Valley debris flow has been selected as a case for addressing loose deposits effects on the whole debris flow process through remote sensing, field investigation and field experiments. Remote sensing interpretation and laboratory experiments were used to obtain the distribution and characteristics of the loose deposits, respectively. A field experiment was conducted to explore the mechanics of slope debris flows, and another field investigation was conducted to obtain the processes of debris flow formation, movement and amplification. The results showed that loose deposits preparation, slope debris flow initiation,gully debris flow confluence and valley debris flow amplification were dominated by the loose deposits.Antecedent droughts and earthquake activities may have increased the potential for loose soil sources in the Aizi Valley, which laid the foundation for debris flow formation. Slope debris flow initiated under rainfall, and the increase in the water content as well as the pore water pressure of the loose deposits were the key factors affecting slope failure. The nine gully debris flows converged in the valley, and the peak discharge was amplified 3.3 times due to a blockage and outburst caused by a large boulder. The results may help in predicting and assessing regional debris flows in dry-hot and seismic-prone areas based on loose deposits, especially considering large boulders.     

四川省小流域泥石流危险性评价

泥石流危险性评价是泥石流防灾减灾的重要内容。本文以四川省为研究区,以DEM为数据源,通过提取水流方向,计算汇流累积量,实现四川省小流域划分。基于收集的已查明泥石流流域资料,分析了泥石流孕灾环境与成灾特点,选择流域高差、流域面积为指标,建立基于能量条件的潜势泥石流流域判识模型,对划分的小流域进行判识,识别出7798个小流域具备泥石流发生所需能量条件,面积为31.1×10⁴ km²,占四川省总面积的64.18 %。进而建立了泥石流危险性评价指标体系和可拓物元模型,开展了小流域泥石流危险性评价,划分了危险度等级,得到中度、高度、极高危险区的小流域个数分别为1946、1725和1002个,面积分别为9.1×10⁴、7.7×10⁴和3.4×10⁴ km²,中度以上危险区面积共20.2×10⁴ km2,占四川省总面积的41.67%。最后对评价结果可靠性和各等级泥石流危险区在各地市级行政区、各大流域的分布进行了分析。其结果对促进泥石流判识与危险性评价理论,区域泥石流防灾减灾与山区可持续发展等具有重要的理论和现实意义。