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.     

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.     

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.     

牛眠沟流域泥石流形成条件、发展趋势及其防治探讨

汶川地震使牛眠沟流域变成了一个多泥石流灾害的区域,到目前为止,已经发生7次大规模泥石流灾害。经研究,该流域诱发泥石流灾害的3个基本条件非常发育,而目前正处于发展阶段,在震后5~10年内的雨季,极易发生危害性较大、冲击力较强的大规模降雨型泥石流。如果考虑直接利用主沟内沉积的泥石流堆积物,在主沟堆积区的中上游左右两侧修建约为沟宽一半的交替式拦挡土石坝、在中下游修建与沟宽大致相等的拦挡土石坝及在相邻两坝间修建与沟向相对且具有一定坡度斜坡的土木防治工程,可实现消能与耗能双重目的。     

Characteristics of earthquake-triggered landslides and post-earthquake debris flows in Beichuan County

Field investigations and aerial photography after the earthquake of May 12,2008 show a large number of geo-hazards in the zone of extreme earthquake effects.In particular,landslides and debris flows,the geo-hazards that most threaten post-disaster reconstruction,are widely distributed.We describe the characteristics of these geo-hazards in Beichuan County using high-resolution remote sensing of landslide distribution,and the relationships between the area and volume of landslides and the peak-discharges of debris flows both pre-and post-earthquake.The results show:1) The concentration(defined as the number of landslide sources per unit area:Lc) of earthquaketriggered landslides is inversely correlated with distance from the earthquake(DF) fault.The relationship is described by the following equation:Lc = 3.2264exp(-0.0831DF)(R2 = 0.9246);2) 87 % of the earthquake-triggered landslides were less than 15×104 m2 in area,and these accounted only for 50% of the total area;84% of the landslide volumes were less than 60×104 m3,and these accounted only for 50% of the total volume.The probability densities of the area and volume distributions are correlated:landslide abundance increases with landslide area and volume up to maximum values of 5 × 104 m2 and 30 × 104 m3,respectively,and then decreases exponentially.3) The area(AL) and volume(VL) of earthquake-triggered landslides are correlated as described with the following equation:VL=6.5138AL1.0227(R2 = 0.9131);4) Characteristics of the debris flows changed after the earthquake because of the large amount of landslide material deposited in the gullies.Consequently,debris flow peak-discharge increased following the earthquake as described with the following equation:Vpost = 0.8421Vpre1.0972(R2 = 0.9821)(Vpre is the peak discharge of pre-earthquake flows and the Vpost is the peak discharge of post-earthquake flows).We obtained the distribution of the landslides based on the above analyses,as well as the magnitude of both the landslides and the post-earthquake debris flows.The results can be useful for guiding post-disaster reconstruction and recovery efforts,and for the future mitigation of these geo-hazards.However,the equations presented are not recommended for use in site-specific designs.Rather,we recommend their use for mapping regional seismic landslide hazards or for the preliminary,rapid screening of sites.     

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.     

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.     

Comparison of the entrainment rate of debris flows in distinctive triggering conditions

Debris flows can be extremely destructive because they can increase in magnitude via progressive entrainment. In this paper, a total of 18 landslide-type debris flows and 268 channelized debris flows in Wenchuan earthquake and Taiwan region, as well as other regions were collected to analyze the entrainment rate of debris flows in each triggering condition. Results show that there is a power relationship between volume of initial triggered mass and final deposited debris for landslide type debris flow. The debris flows during 2008 and 2013 in Wenchuan earthquake-region have smaller entrainment rate than that from 2001 t0 2009 in Taiwan. The entrainment rate of debris flow events from 2001 to 2009 in Taiwan shows a decaying tendency as elapsed time. Comparison of the entrainment rate in the two earthquake-hit regions with other regions proves that entrainment rate has a close relation with major sediment availability and secondary rainstorm conditions.     

Hybrid simulation of the initiation and runout characteristics of a catastrophic debris flow

On 13 August 2010, a catastrophic debris flow with a volume of 1.17 million m³ occurred in Xiaojiagou Ravine near Yingxiu town of Wenchuan county in Sichuan Province, China. The main source material was the landslide deposits retained in the ravine during the 2008 Wenchuan earthquake. This paper describes a two-dimensional hybrid numerical method that simulates the entire process of the debris flow from initiation to transportation and finally to deposition. The study area is discretized into a grid of square zones. A two dimensional finite difference method is then applied to simulate the rainfall-runoff and debris flow runout processes. The analysis is divided into three steps; namely, rainfall-runoff simulation, mixing water and solid materials, and debris flow runout simulation. The rainfall-runoff simulation is firstly conducted to obtain the cumulative runoff near the location of main source material and at the outlet of the first branch. The water and solid materials are then mixed to create an inflow hydrograph for the debris flow runout simulation. The occurrence time and volume of the debris flow can be estimated in this step. Finally the runout process of the debris flow is simulated. When the yield stress is high, it controls the deposition zone. When the yield stress is medium or low, both yield stress and viscosity influence the deposition zone. The flow velocity is largely influenced by the viscosity. The estimated yield stress by the equation, τ _y = ρghsin θ, and the estimated viscosity by the equation established by Bisantino et al. (2010) provide good estimates of the area of the debris flow fan and the distribution of deposition depth.     

Risk assessment of highways affected by debris flows in Wenchuan earthquake area

Since the Wenchuan earthquake in China on May 12^th, 2008, highways in earthquake-affected areas have been frequently interrupted by debris flows. We analyzed the hazard effect modes and damage processes along highways and developed three key indexes, scale of debris flows, deposits on highways and river blockage, to describe quantitatively the highway disasters. By combining the empirical methods and the actual terrain conditions, we proposed new methods to determine the value of hazard indexes. In addition, we used the economic value and resistance of highway as vulnerability assessment indexes, then determined the specific subindexes for the subgrade, bridges and culverts, and developed a way for the quantified vulnerability zoning. Moreover, we proposed the assessment and mapping methods for highway risk. The risk is described into 5 grades: extremely low risk, low risk, middle risk, high risk and extremely high risk. We applied these methods in a case study carried out on provincial highway S303 from Yingxiu Town to Wolong Town, in Wenchuan County. Analysis of debris flow risk for the whole highway, showed that the total length of highway in extremely low risk area was 28.26 km, 4.83 km in low risk area, 8.0 km in middle risk area, 3.65 km in high risk area, and 3.06 km in extremely high risk area. The assessment results are consistent with the field survey data which reflected the disaster situation. This risk method can be used objectively to evaluate the debris-flow risk along highways, and is useful for highway reconstruction in mountainous areas suffering from active debris flows.     

A catastrophic debris flow in the Wenchuan Earthquake area,July 2013: characteristics,formation, and risk reduction

In the Wenchuan Earthquake area,many co-seismic landslides formed blocking-dams in debris flow channels. This blocking and bursting of landslide dams amplifies the debris flow scale and results in severe catastrophes. The catastrophic debris flow that occurred in Qipan gully(Wenchuan,Southwest China) on July 11,2013 was caused by intense rainfall and upstream cascading bursting of landslide dams. To gain an understanding of the processes of dam bursting and subsequent debris flow scale amplification effect,we attempted to estimate the bursting debris flow peak discharges along the main gully and analyzed the scale amplification process. The results showed that the antecedent and triggering rainfalls for 11 July debris flow event were 88.0 mm and 21.6 mm,respectively. The event highlights the fact that lower rainfall intensity can trigger debris flows after the earthquake. Calculations of the debris flow peak discharge showed that the peak discharges after the dams-bursting were 1.17–1.69 times greater than the upstream peak discharge. The peak discharge at the gully outlet reached 2553 m~3/s which was amplified by 4.76 times in comparison with the initial peak discharge in the upstream. To mitigate debris flow disasters,a new drainage channel with a trapezoidal V-shaped cross section was proposed. The characteristic lengths(h1 and h2) under optimal hydraulic conditions were calculated as 4.50 m and 0.90 m,respectively.     

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.     

Debris flow susceptibility analysis based on the combined impacts of antecedent earthquakes and droughts: a case study for cascade hydropower stations in the upper Yangtze River,China

The upper Yangtze River region is one of the most frequent debris flow areas in China. The study area contains a cascade of six large hydropower stations located along the river with total capacity of more than 70 million kilowatts. The purpose of the study was to determine potential and dynamic differences in debris flow susceptibility and intensity with regard to seasonal monsoon events. We analyzed this region’s debris flow history by examining the effective peak acceleration of antecedent earthquakes, the impacts of antecedent droughts, the combined effects of earthquakes and droughts, with regard to topography, precipitation, and loose solid material conditions. Based on these factors, we developed a debris flow susceptibility map. Results indicate that the entire debris flow susceptibility area is 167,500 km², of which 26,800 km² falls within the high susceptibility area, with 60,900 km² in medium and 79,800 km² are in low susceptibility areas. Three of the six large hydropower stations are located within the areas with high risk of debris flows. The synthetic zonation map of debris flow susceptibility for the study area corresponds with both the investigation data and actual distribution of debris flows. The results of debris flow susceptibility provide base-line data for mitigating, assessing, controlling and monitoring of debris flows hazards.     

Assessment of prospective hazards resulting from the 2017 earthquake at the world heritage site Jiuzhaigou Valley,Sichuan, China

On August 8, 2017, a Ms = 7.0 magnitude earthquake occurred in the Jiuzhaigou Valley, in Sichuan Province, China(N: 33.20°, E: 103.82°). Jiuzhaigou Valley is an area recognized and listed as a world heritage site by UNESCO in 1992. Data analysis and field survey were conducted on the landslide, collapse, and debris flow gully, to assess the coseismic geological hazards generated by the earthquake using an unmanned aerial vehicle(UAV), remote-sensing imaging, laser range finders, geological radars, and cameras. The results highlighted the occurrence of 13 landslides, 70 collapses, and 25 potential debris flow gullies following the earthquake. The hazards were classified on the basis of their size and the potential property loss attributable to them. Consequently, 14 large-scale hazards, 30 medium-sized hazards, and 64 small hazards accounting for 13%, 28%, and 59% of the total hazards, respectively, were identified. Based on the variation tendency of the geological hazards that ensued in areas affected by the Kanto earthquake(Japan), Chi-chi earthquake(Taiwan China), and Wenchuan earthquake(Sichuan China), the study predicts that, depending on the rain intensity cycle, the duration of geological hazard activities in the Jiuzhaigou Valley may last over ten years and will gradually decrease for the following five to ten yearsbefore returning to pre-earthquake levels. Thus,necessary monitoring and early warning systems must be implemented to ensure the safety of residents,workers and tourists during the construction of engineering projects and reopening of scenic sites to the public.     

Case history of the disastrous debris flows of Tianmo Watershed in Bomi County,Tibet, China: Some mitigation suggestions

Debris flows and landslides, extensively developing and frequently occurring along Parlung Zangbo, seriously damage the Highway from Sichuan to Tiebt(G318) at Bomi County. The disastrous debris flows of the Tianmo Watershed on Sept. 4, 2007, July 25, 2010 and Sept. 4, 2010, blocked Parlung Zangbo River and produced dammed lakes, whose outburst flow made 50 m high terrace collapse at the opposite bank due to intense scouring on the foot of the terrace. As a result, the traffic was interrupted for 16 days in 2010 because that 900 m highway base was destructed and 430 m ruined. These debris flows were initiated by the glacial melting which was induced by continuous higher temperature and the following intensive rainfall, and expanded by moraines along channels and then blocked Parlung Zangbo. At the outlet of watershed,the density, velocity and peak discharge of debris flow was 2.06 t/m3, 12.7 m/s and 3334 m3/s, respectively. When the discharge at the outlet and the deposition volume into river exceeds 2125 m3/s and 126×103 m3, respectively, debris flow will completely blocked Parlung Zangbo. Moreover,if the shear stress of river flow on the foot of terrace and the inclination angel of terrace overruns 0. 377 N/m2 and 26°, respectively, the unconsolidated terrace will be eroded by outburst flow and collapse. It was strongly recommended for mitigation that identify and evade disastrous debris flows, reduce the junction angel of channels between river and watershed, build protecting wall for highway base and keep appropriate distance between highway and the edge of unconsolidated terrace.     

Artificial Neural Network-based prediction of glacial debris flows in the ParlungZangbo Basin,southeastern Tibetan Plateau,China

Accurate prediction on geological hazards can prevent disaster events in advance and greatly reduce property losses and life casualties.Glacial debris flows are the most serious hazards in southeastern Tibet in China due to their complexity in formation mechanism and the difficulty in prediction.Data collected from 102 glacier debris flow events from 31 gullies since 1970 and regional meteorological data from 1970 to 2019 in ParlungZangbo River Basin in southeastern Tibet were used for Artificial Neural Network(ANN)-based prediction of glacial debris flows.The formation mechanism of glacial debris flows in the ParlungZangbo Basin was systematically analyzed,and the calculations involving the meteorological data and disaster events were conducted by using the statistical methods and two layers fully connected neural networks.The occurrence probabilities and scales of glacial debris flows(small,medium,and large)were predicted,and promising results have been achieved.Through the proposed model calculations,a prediction accuracy of 78.33%was achieved for the scale of glacial debris flows in the study area.The prediction accuracy for both large-and medium-scale debris flows are higher than that for small-scale debris flows.The debris flow scale and the probability of occurrence increase with increasing rainfall and temperature.In addition,the K-fold cross-validation method was used to verify the reliability of the model.The average accuracy of the model calculated under this method is about 93.3%,which validates the proposed model.Practices have proved that the combination of ANN and disaster events can provide sound prediction on geological hazards under complex conditions.     

Characteristics of clustering debris flows in Wenchuan earthquake zone

Clustering debris-flow events, namely many debris flows simultaneously triggered by a regional rainstorm in a large-scale mountainous area, occurred in four regions of Wenchuan earthquake stricken areas in 2008 and 2010. The characteristics of the clustering debris flows are examined with regard to triggering rainfall, formation process, and relationship with the earthquake by field survey and remote sensing interpretation. It is found that the clustering events occurred nearly at the same time with the local peak rainstorms, and the rainfall intensity-duration bottom limit line for clustering debris flows is higher than the worldwide line. It means that more rainfall is needed for the occurrence of the clustering debris flows. Four kinds of major formation processes for these debris flows are summarized： tributary-dominated, mainstream- dominated, transformation from slope failures, and mobilization or liquefaction of landslide. The four regions has a spatial correlation with the strong- quake-influenced zone with the peak ground acceleration = 0.2 g and the seismic intensity 〉 X.     

Slope debris flows in the Wenchuan Earthquake area

Avalanches and landslides, induced by the Wenchuan Earthquake on May 12, 2008, resulted in a lot of disaggregated, solid material on slopes that could be readily mobilized as source material for debris flows. Rainstorms triggered numerous slope debris flows with great damage to highways and rivers over the subsequent two years. Slope debris flows (as opposed to channelized debris flows) are defined as phenomena in which high-concentration mixtures of debris and water flow down slopes for short distances to highways and river banks. Based on field investigations and measurements of 19 slope debris flows, their main characteristics and potential mitigation strategies were studied. High rainfall intensity is the main triggering factor. Critical rainfall intensities for simultaneous occurrence of single, several and numerous slope debris flow events were 20 mm/day, 30mm/day, and 90 mm/day, respectively. Field investigations also revealed that slope debris flows consist of high concentrations of cobbles, boulders and gravel. They are two-phase debris flows. The liquid phase plays the role of lubrication instead of transporting medium. Solid particles collide with each other and consume a lot of energy. The velocities of slope debris flows are very low, and their transport distances are only several tens of meters. Slope debris flows may be controlled by construction of drainage systems and by reforestation.     

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.     

A Modified Certainty Coefficient Method（M-CF） for Debris Flow Susceptibility Assessment：A Case Study for the Wenchuan Earthquake Meizoseismal Areas

In the meizoseismal areas hit by the China Wenchuan earthquake on May 12, 2008, the disasterprone environment has changed dramatically, making the susceptibility assessment of debris flow more complex and uncertain. After the earthquake, debris flow hazards occurred frequently and effective susceptibility assessment of debris flow has become extremely important. Shenxi gully in Du Jiangyan city, located in the meizoseismal areas, was selected as the study area. Based on the research of disaster-prone environment and the main factors controlling debris flow, the susceptibility zonations of debris flow were mapped using factor weight method（FW）, certainty coefficient method（CF） and geomorphic information entropy method（GI）. Through comparative analysis, the study showed that these three methods underestimated susceptible degree of debris flow when used in the meizoseismal areas of Wenchuan earthquake. In order to solve this problem, this paper developed a modified certainty coefficient method（M-CF） to reflect the impact of rich loose materials on the susceptible degree of debris flow. In the modified method, the distribution and area of loose materials were obtained by field investigations and postearthquake remote sensing image, and four data sets, namely, lithology, elevation, slop and aspect, wereused to calculate the CF values. The result of M-CF method is in agreement with field investigations and the accuracy of the method is satisfied. The method has a wide application to the susceptibility assessment of debris flow in the earthquake stricken areas.