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.     

Mechanisms involved in triggering debris flows within a cohesive gravel soil mass on a slope: a case in SW China

The triggering mechanisms of debris flows were explored in the field using artificial rainfall experiments in two gullies, Dawazi Gully and Aizi Gully, in Yunnan and Sichuan Provinces, China, respectively. The soils at both sites are bare, loose and cohesive gravel-dominated. The results of a direct shear test, rheological test and back-analysis using soil mass stability calculations indicate that the mechanisms responsible for triggering debris flows involved the decreases in static and dynamic resistance of the soil. The triggering processes can be divided into 7 stages: rainfall infiltration, generation of excess runoff, high pore water pressure, surface erosion, soil creep, soil slipping, debris flow triggering and debris flow increment. In addition, two critical steps are evident: (i) During the process of the soil mass changing from a static to a mobile state, its cohesion decreased sharply (e.g., the cohesion of the soil mass in Dawazi Gully decreased from 0.520 to 0.090 kPa, a decrease of 83%). This would have reduced the soil strength and the kinetic energy during slipping, eventually triggered the debris flow. (ii) When the soil mass began to slip, the velocity and the volume increment of the debris flow fluctuated as a result of the interaction of soil resistance and the sliding force. The displaced soil mass from the source area of the slope resulted in the deposition of a volume of soil more than 7 - 8 times greater than that in the source area.     

Characteristic rainfall for warning of debris flows

A characteristic rainfall is introduced to overcome the difficulties encountered in determining a critical rainfall value for triggering debris flow.The characteristic value is defined as the rainfall at which debris-flow occurrence probability shows a rapid increase,and can be used as a warning rainfall threshold for debris flows.Investigation of recorded debris flows and 24-hour rainfall data at Jiangjia basin,Yunnan Province,in southwestern China,demonstrates the existence of such a characteristic rainfa...     

Magnitude-frequency relationship of debris flows in the Jiangjia Gully,China

The magnitude-frequency(MF) relationship of debris flows is the basis for engineering designs and risk quantification. However, because of the lack of debris flow monitoring data, research progress in this area has been relatively slow. The MF relationship of debris flows in Jiangjia Gully, Yunnan Province was evaluated based on a regression analysis of 178 debris flow events that occurred from 1987-2004. The magnitude-cumulative frequency(MCF) relationship of the debris flows in the Jiangjia Gully is consistent with the linear logarithmic transformation function. Moreover, observed data for debris flows in Hunshui Gully of Yunnan Province and Huoshao Gully, Liuwan Gully, and Niwan Gully of Gansu Province were used to verify the function. The results showed that the MCF relationship of highfrequency debris flows is consistent with the power law equation, although the regression coefficients in the equation are considerably different. Further analysis showed a strong correlation between the differences in the constants and the drainage area and daily maximum precipitation.     

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.     

Influence of fine sediment on the fluidity of debris flows

Debris flows include a great diversity of grain sizes with inherent features such as inverse grading, particle size segregation, and liquefaction of fine sediment. The liquefaction of fine sediment affects the fluidity of debris flows, although the behavior and influence of fine sediment in debris flows have not been examined sufficiently. This study used flume tests to detect the effect of fine sediment on the fluidity of laboratory debris flows consisting of particles with various diameters. From the experiments, the greatest sediment concentration and flow depth were observed in the debris flows mixed with fine sediment indicating increased flow resistance. The experimental friction coefficient was then compared with the theoretical friction coefficient derived by substituting the experimental values into the constitutive equations for debris flow. The theoretical friction coefficient was obtained from two models with different fine-sediment treatments: assuming that all of the fine sediments were solid particles or that the particles consisted of a fluid phase involving pore water liquefaction. From the comparison of the friction coefficients, a fully liquefaction state was detected for the fine particle mixture. When the mixing ratio and particle size of the fine sediment were different, some other cases were considered to be in a partially liquefied transition state. These results imply that the liquefaction of fine sediment in debris flows was induced not only by the geometric conditions such as particle sizes, but also by the flow conditions.     

Experimental study of entrainment behavior of debris flow over channel inflexion points

On-spot observation and field reconnaissance of debris flows have revealed that inflexion points in the longitudinal profile of a movable channel may easily become unstable points that significantly affect their entrainment behavior. In this study, small-scale flume experiments were performed to investigate the entrainment characteristics of debris flows over two types of inflexion points, namely, a convex point, which has an upslope gradient that is less than the downslope gradient; and a concave point, which has an upslope gradient that is greater than the downslope gradient. It was observed that when debris flowed over a convex point, the entrainment developed gradually and progressively from the convex point in the downstream direction, and the primary control factors were the slope gradient and friction angle. Conversely, when debris flowed over a concave point, the entrainment was characterized by impacting and impinging erosion rather than traditional hydraulic erosion, and the impingement angle of the flow significantly determined the maximum erosion depth and outflow exit angle. An empirical relationship between the topography change and the control factors was obtained from the experimental data.     

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.     

Grain composition and erosive equilibrium of debris flows

Debris flows consist of grains of various sizes ranging from 10~(-6) m ～1 m. Field observations in the Jiangjia Gully (JJG) and other sites throughout China indicate that the grain size distribution of sediment in debris flows can be characterized by an exponential function fit to the cumulative distribution. The exponent value for the function varies by location and may be useful in distinguishing between debris flows from different valleys. For example, minimum values and ranges of the exponent are associated with the high frequency of debris flows in the JJG. Furthermore, the distribution presents piecewise fractality (i.e. scaling laws hold in various ranges of the grain size) and we propose that the fractal structure determines the matrix and that the fractal dimension plays a crucial role in material exchange between a debris flow and the substrate it flows over. Finally, the empirical data support an exponential relation between grain composition and non-dimensional shear stress for the critical state of the channel. Overall we propose a material-determinism approach to studying debris flows which contrasts with the enviro-determinism that has dominated much recent work in this field.     

Gravel accumulation in deposits of viscous debris flows with hyper-concentration

According to the observational data of viscous debris flows with hyper-concentration, debris flows can be classified into three types： high-viscous, viscous, and sub-viscous debris flows. Distinct formation mechanism of different graded bedding structures in deposits of viscous debris flows was analyzed in this paper by using their yield-stress ratio and flow plug ratio. This paper specially analyzed the effect of Weissenberg which the gravels in squirm condition of hyper-concentration viscous flows would tend to move vertically, and the formation mechanism of the gravels accumulated at surface was also studied. The analysis in this paper can establish a foundation for the studies on differentiation of bedding structures of debris flow deposits and studies on dynamic parameters of debris flows.     

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.     

Potential sediment sources identification of debris flows in the Jiangjia Gully,China

It is of great significance for gully prevention and management to identify the potential sediment source of debris flow. Debris flow in a gully always originates from tributaries that have different gravity potential energies and sediment condition. In this study, tributaries of the Jiangjia Gully(JJG) in Yunnan province, China, are taken as the study area to determine the possible sediment sources of debris flow. It was found that tributaries with a high evolution index(EI, the integral of the hypsometric curve) always had high gravity potential energy, which favors the occurrence of landslide activity. Furthermore, the relationship between sediment distribution, gravity potential energy, and EI is compared, respectively. The results showed that the EI had a greater influence on the occurrence of landslides, and sediments were concentrated in tributaries with EI between 0.5 and 0.6. Accordingly, tributaries with EI 0.5 were identified as the sediment sources of debris flow. In addition, the shape of a tributary was related to EI and can reflect the condition of water and sediment storage.     

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.     

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.     

Two-dimensional numerical model for debris flows in the Jiangjia Gully, Yunnan Province

Debris flows are recurrent natural hazards in many mountainous regions.This paper presents a numerical study on the propagation of debris flows in natural erodible open channels,in which the bed erosion and sedimentation processes are important.Based on the Bingham fluid theory,a mathematical model of the two-dimensional non-constant debris flow is developed.The governing equations include the continuity and momentum conservation equations of debris flow,the sediment convection-diffusion equation,the bed erosion-deposition equation and the bed-sediment size gradation adjustment equation.The yield stress and shear stress components are included to describe the dynamic rheological properties.The upwind control-volume Finite Volume Method (FVM) is applied to discretize the convection terms.The improved SIMPLE algorithm with velocity-free-surface coupled correction is developed to solve the equations on non-orthogonal,quadrilateral grids.The model is applied to simulate a debris flow event in Jiangjia Gully,Yunnan Province and to predict the flow pattern and bed erosion-deposition processes.The results show the effectiveness of the proposed numercial model in debris flow simulation and potential hazard analysis.     

Measuring the internal velocity of debris flows using impact pressure detecting in the flume experiment

Measuring the internal velocity of debris flows is very important for debris flow dynamics research and designing debris flow control works. However, there is no appropriate method for measuring the internal velocity because of the destructive power of debris flow process. In this paper, we address this problem by using the relationship between velocity and kinetic pressure, as described by surface velocity and surface kinetic pressure data. Kinetic pressure is the difference of impact pressure and static pressure. The former is detected by force sensors installed in the flow direction at the sampling section. Observations show that static pressure can be computed using the formula for static water pressure by simply substituting water density for debris flow density. We describe the relationship between surface velocity and surface kinetic pressure using data from seven laboratory flume experiments. It is consistent with the relationship for single phase flow, which is the measurement principle of the Pitot tube.     

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.     

Geomorphologic analysis and physico-dynamic characteristics of Zhatai-Gully debris flows in SW China

Zhatai gully is a typical debris flow channel in Butuo county of Sichuan province, southwestern China. The geomorphologic features are analyzed and the physical-dynamic characteristics are discussed on the basis of field investigation and laboratory tests. Geomorphologic analysis indicates that Zhatai-gully drainage in relation to debris flow can be divided into source area, transport area, and deposition area. The source area has a steep slope and has very limited vegetation cover, which favors runoff, allowing loose solid materials to be mobilized easily and rapidly. In the transport area, there are many small landslides, lateral lobes, and loose materials distributed on both banks. These landslides are active and constantly providing abundant source of soils for the debris flows. In the deposition area, three old debris-flow deposits of different ages can be observed. The dynamic calculation shows that within the recurrence intervals of 50 and 100 years, debris flow discharges are 155.77m³/s and 1y8.19m³/s and deposition volumes are 16.39 x 10⁴ m³ and 18.14 x 10⁴ m³, respectively. The depositional fan of an old debris flow in the outlet of the gully can be subdivided into six layers. There are three debris flow deposits on left and two on the right side of the gully. Grain-size tests of sediments from the soil, gulley bed deposits, and the fresh and old debris flow deposits showed that high amounts of clay and fine gravel were derived from the soil in the source area whereas much of the gravel fraction were sourced from the gully bed deposits. Comprehensive analysis indicates that Zhatai gully is viscous debris-flow gully with moderate to high frequency and moderate to large magnitude debris flows. The risk of a debris flow disaster in Zhatai-gully is moderate and poses a potential threat to the planned hydroelectric dam. Appropriate engineering measures are suggested in the construction and protection of the planned hydroelectric station.     

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.