1. Introduction
The natural hazard chain refers to a series of secondary hazards induced by a single natural hazard [
1]. Compared with a single hazard, a hazard chain has the characteristics of a long time scale, broad scope, and immense destructiveness and results in enormous losses of life and economic losses, and harmful social and environmental effects [
2]. Heavy rains, debris flows, earthquakes, landslides, landslide dams, and outburst floods are all individual hazards, but sometimes, they can be in a hazard chain and cause more severe hazards. A landslide is usually caused by heavy rains or geological effects such as earthquakes and occurs suddenly, with huge impacts on human production and life [
3]. When there is a river in the direction of the landslide, the landslide can easily block the river and form a landslide dam. The landslide dam causes the upstream water level to keep rising and also causes the submersion of the upstream residential area. Once the landslide dam breaches, the large amount of water accumulated in the barrier lake will be released quickly and will cause enormous hazards downstream.
In recent years, frequent earthquakes and heavy rains have led to a global outbreak of landslides [
4]. The Tangjiashan landslide dam formed in the 2008 Wenchuan earthquake and collapsed 29 days later. The dam forced nearly 200,000 people downstream to evacuate [
5]. In 2009, heavy rainfall caused by a typhoon led to landslides and debris flows in Hsiaolin Village. The landslides buried the village and blocked the gorge of the Cishan River, and the peak outflow rate reached 70,649 m
3/s after the landslide dam breached [
6]. In 2018, landslides blocked the Jinsha River twice at the same location. The outburst flood caused road damage and extensive damage to large areas of farmland and houses [
7].
Therefore, a risk assessment method needs to be proposed to reduce the economic loss, loss of life, and social and environmental effects of the hazard chain. In the existing research, scholars from China and overseas have mostly used mathematical methods, such as the analytic hierarchy process, the Poisson probability model, and information models, to evaluate the vulnerability and risk of landslide disasters. However, the existing risk assessments have focused primarily on single hazards, including forecasts of landslide movements [
8,
9], dam formation [
10,
11,
12,
13], dam stability [
14,
15,
16,
17], and outburst floods [
3,
10,
18,
19]. With the development of science and technology, more and more advanced methods are used for forecasting hazards. Ji et al. [
20] proposed the inverse FORM algorithm and, based on this, a geographic information system extension tool was developed for probabilistic physical modeling and forecasting landslides [
21]. Panagoulia [
22,
23] analyzed the responses of the medium-sized mountainous Mesochora catchment to climate change and proposed a multi-stage method for selecting input variables for ANN forecasts of river flows.
For assessing the risk of a hazard chain, a conceptual risk assessment model of a regional hazard chain was proposed based on Newmark’s permanent deformation model and applied to predicting a disaster chain induced by the Wenchuan earthquake [
24]. Dong [
25] summarized the hazard chain mode of river blockages induced by avalanches on the basis of a hazard chain structure. A risk assessment of the Layue landslide dam was carried out from the aspects of risk and vulnerability. These studies have mainly focused on the concept of a hazard chain, but have yet to emphasize a systematic risk assessment model. Although the remedy for a hazard chain lies in blocking it before expansion/transformation, sometimes, it is inevitable [
26]. Therefore, a risk assessment model is needed to reduce the losses caused by the hazard chain.
This study aimed to explore a risk assessment model of hazard chains and to apply it to the case of landslides. A database was established for statistical analysis. Based on existing research, formulas were established for assessing the risk of river blockages caused by landslides, the stability of landslide dams, and the peak flood discharge after dam breakage. A risk assessment model of the hazard chain of landslide–river blockage–outburst flood was proposed, with explicit formulas. The case of the Baige landslide was investigated to verify the applicability of the model.
2. Statistical Analyses of Individual Hazards
The evolution of the investigated hazard chain caused by landslides can be divided into three stages: river blockage, dam breakage, and the outburst flood. The risks of these three hazards are accessed individually in this section. A database including 57 landslide events, 50 historically documented landslide dams, and 34 landslide dam breaches was established with relatively complete and accurate data. Some improved risk assessment models were proposed and calibrated on the basis of the database.
2.1. Risk Assessment of River Blockage Caused by Landslides
The river blockage caused by landslides can be classified into complete blockage and partial blockage, as shown in
Figure 1. In the case of a partial blockage, the landslide deposit blocks part of the river channel and the river still flows through the unblocked channel. In the case of a complete blockage, the landslide deposit completely blocks the river channel and forms a landslide dam, which promotes the formation of a barrier lake. Therefore, assessing the degree of river blockage is the primary task of risk assessment.
Tacconi Stefanelli et al. [
13] investigated the relationship between the formation of landslide dams and two critical conditions, namely landslide volume and valley width, and proposed an assessment index, named the morphological obstruction index (MOI), as follows:
where
VL denotes the landslide’s volume (m
3) and
WV denotes the valley’s width (m). When
, a landslide dam forms; when
, a landslide dam does not form; and landslide dams are uncertain when
. The advantage of this formula is that the valley’s width and the landslide’s volume can easily be measured, so it is suitable for rapid assessments of the risk of river blockage.
However, the accuracy of the MOI is poor. Only 57.9% of the cases in the database could be assessed correctly (
Table A1) by this method. This is because the risk of a landslide blocking the river is also related to many other factors, such as the valley’s shape, the riverbank’s slope, and the landslide’s materials. The accuracy needs to be improved by taking more factors into consideration.
A landslide formed on a steeper slope is more likely to induce a landslide dam, as the landslide generates more kinetic energy [
10]. According to the database, the relationship between the steepness of the slope and the degree of river blockage can be statistically analyzed, as shown in
Figure 2. This illustrates that a landslide with a steeper slope is more likely to completely block the river and form a landslide dam. Moreover, landslide materials with larger particles are also more likely to accumulate directly at the bottom of the river and form a landslide dam. It is also easier for a landslide dam to form in a V-shaped valley than in a U-shaped valley [
27].
Hu et al. [
28] and Tacconi Stefanelli et al. [
29] collected many cases of river blockages caused by landslides in China and Italy. Some cases with complete data were selected to establish a database of 57 events. According to the database (
Table A1), an index based on the MOI was proposed, named the Landslide Blocking River Index (LBRI). The valley’s shape, the landslide’s materials, and the steepness of the slope were considered in the LBRI to make it more accurate and enhance its applicability. Because some cases lacked a description of the valley’s shape, the shapes of these cases were considered to be V-shape valleys for a conservative evaluation:
where
φ denotes the steepness of the landslide slope (°) and
AC denotes the “accumulation coefficient”, which is related to the valley’s shape and the particle size of landslide materials. According to the analysis results of database, different values of
AC are suggested for various landslide materials and different valley shapes, as shown in
Table 1.
In
Figure 3, the results calculated for the LBRI were divided into two different domains, i.e., complete blockage and partial blockage. According to the statistical analysis of these data, when
, the river is partly blocked; when
, the river is completely blocked. When the number of correct cases was divided by the total number of cases, the formula’s accuracy reached 96.5%.
2.2. Assessing the Stability of Landslide Dams
A landslide dam’s stability increases with the dam’s volume or the ratio of the dam’s width to its height. Obviously, the longer the dam is, with the higher hydraulic thrust, the lower the stability of the dam will be. Meanwhile, as the ratio of the reservoir capacity of the barrier lake to the dam’s volume increases, the landslide dam’s stability decreases [
30]. Canuti et al. [
14] proposed the blockage index (BI) to assess a dam’s stability. The formulation of BI is expressed as follows:
where
VD is the dam’s volume (m
3) and
SC is the catchment area (10
6 m
2). Canuti et al. [
14] suggested that when
, the dam is stable; when
, the stability is unsure; and when
, the dam is unstable.
The dam volume is easily measured, and the catchment area is easily calculated. The BI is instructive for the rapid assessment of the stability of landslide dams, so it is widely used [
13,
16,
31]. However, the accuracy of the BI is poor. Only 50.0% of the cases in the database could be assessed correctly (
Table A2).
Xu [
17] developed a rapid assessment model of the stability of landslide dams by collecting the geometric information and breach times of 110 landslide dams worldwide. The model consists of two parts [
17], i.e., Fisher’s discriminant model
YA and the logistics regression model
ZB.
where
Vl is the reservoir capacity of the barrier lake (m
3) and
Hd is the dam’s height (m).
Fisher’s discriminant model considered the influences of the dam’s height Hd, the landslide dam’s volume VD, the reservoir capacity of the barrier lake Vl, and the catchment area SC on the stability of landslide dams. The logistic regression model considered the influences of the dam’s height Hd, the landslide dam’s volume VD, and the catchment area SC on the stability of landslide dams. In Fisher’s discriminant model, if YA < 0, the landslide dam is stable; if YA > 0, the landslide dam is unstable. In the logistic regression model, if ZB < 0.5, the landslide dam is stable; if ZB > 0.5, the landslide dam is unstable.
Dam materials originate from the landslide’s materials, so they have almost the same material composition. Dams made up of earth or soft rock have a low shear strength and are prone to breaching, and dams made up of rock or debris are more stable [
32]. Therefore, it is also necessary to consider the influence of the dam materials.
2.2.1. The Landslide Dam Stability Index
Shan et al. [
16] established a database for rapid predictions of dam stability. In the database, 50 documented historical landslide dams with complete data were chosen for assessing the stability of landslide dams, as listed in
Table A2. Based on Xu’s model and the database, a new index named the landslide dam stability index (LDSI) was proposed, which is expressed as follows:
where
Wd is the dam’s width (m),
Ld is the dam’s length (m), and
M is the material coefficient, which describes the influence of the material on the dam’s stability. If the geometries of two landslide dams are exactly the same, a landslide dam formed by rock is more stable than one formed by earth, so the deposited material significantly affects the stability. Different values of
M have been suggested for various landslide materials according to the results of analyzing the database. The suggested values of
M are listed in
Table 2.
Figure 4 shows that when
, the landslide dam is stable, and when
, the landslide dam is unstable. When the number of correct cases was divided by the total number of cases, the formula’s accuracy reached 86%.
Before a landslide occurs, it is difficult to precisely evaluate the dam’s length, height, width, and volume; the reservoir capacity of barrier lake; catchment area; etc. Thus, to assess the risk of the hazard chain triggered by a landslide, it is necessary to estimate these parameters in advance. Some empirical methods are introduced in the following section.
2.2.2. Some Methods of Calculating the Parameters
Equation (6) involves the dam’s height, width, length, and volume; the reservoir capacity of the barrier lake; the catchment area; and the material coefficient. Only the material coefficient is known. Chen et al. [
33] studied the landslide events in Taiwan and proposed a formula for the volume of a dam. However, the formulas were proposed on the basis of rainfall and earthquake conditions. On the basis of the form of this formula and the database, a new formula for calculating the volume of a dam in all cases is as follows:
where
VD has a unit of 10
6 m and
VL has a unit of 10
6 m. The fitting results show that Equation (7) has a correlation coefficient of 0.675.
According to the results when analyzing the database, the dam’s width is mainly affected by the dam’s volume and has little correlation with the dam’s height. A formula for the dam’s width can be fitted as follows:
The fitting results of the database show that Equation (8) has a correlation coefficient of 0.850.
Hu et al. [
28] proposed an empirical formula based on 86 cases to calculate the length of a dam:
where
h is the depth of water in the river (m).
After the volume, width, and length of the dam have been determined, the height of the dam can be easily calculated as follows:
The fitting results of the database show that Equation (10) has a correlation coefficient of 0.799.
The geometric relationship for estimating the reservoir capacity of the barrier lake and the catchment area is shown in
Figure 5. The reservoir capacity of the barrier lake and the catchment area can be calculated as follows:
where
iu is the upstream riverbed’s inclination along the river (rad).
2.3. Assessment of the Peak Flood Discharge after a Dam Breakage
After the dam has been breached, the peak outflow rate occurs at the dam site, and the peak flood discharge decreases with an increase in the distance of the flow. The degree of risk is assessed by comparing the relationship between the peak flood discharge, the check flood discharge, and the average annual discharge in a downstream area.
Based on 12 cases of landslide dams, Costa and Schuster [
10] established the relationship between the peak outflow rate and potential energy:
where
PE is the potential energy (N·m) and
γw is the specific gravity of water (N/m
3).
2.3.1. The Method of Calculating the Peak Outflow Rate
The higher the landslide dam or the larger the reservoir capacity of the dam lake, the more potential energy will be stored in the dam. Because the shape of the breach is difficult to predict, the development process of the breach can be ignored, and the relationship between the influences of some related factors on the peak outflow rate can be estimated on the basis of a statistical analysis of the cases of dam breaches.
The estimation model proposed by Costa and Schuster [
10] should be improved by considering the influencing factors of the dam’s volume and its erodibility. According to the 34 dam breaching cases listed in
Table A3, a new formula is proposed:
where
β is the coefficient of the erodibility related to the dam materials. According to the results of analyzing the database, the suggested values of
β are listed in
Table 3.
Figure 6 shows a comparison between the measured peak outflow rate and the calculated peak outflow rate, and presents a good degree of fitting with a correlation coefficient of 0.934. However, the correlation coefficient of Costa and Schuster’s formula only reached 0.787 in this database.
2.3.2. The Method of Calculating the Downstream Peak Flood Discharge
Li [
34] proposed an empirical formula called the attenuation formula of the peak outflow rate to predict the peak flood discharge somewhere downstream:
where
L0 is the distance from the dam site to somewhere downstream (m),
is the peak flood discharge at
L0 from the dam site (m
3/s), and
Vmax is the maximum average flow velocity during the flood period (m/s). The historical maximum velocity can be used in the areas with detailed data. If there are no data, Li [
34] suggests that 3.0–5.0 m/s can be used in general mountainous areas, 2.0–3.0 m/s in semi-mountainous areas, and 1.0–2.0 m/s in plains.
K is an empirical coefficient. Li [
34] also suggests that
K equals 1.1–1.5 in mountainous areas, 1.0 in semi-mountainous areas, and 0.8–0.9 in plains.
If we compare
to the check flood discharge
and the average annual discharge
, the safety index (
Si) can be obtained. The formulation of
Si is expressed as follows:
5. Discussion
For the case of the Baige landslide, the proposed assessment model was compared with other models, such as Equations (1), (3), and (15), for calculating the risks of river blockage, the dam’s stability, and the peak outflow rate. The results of this comparison are shown in
Table 10. The analyzed results show that the proposed model had a more reasonable dam stability and peak outflow rate than the other two models, so it could have a higher accuracy in assessments of the risk of a landslide hazard chain.
The high accuracy was because more factors were considered in the assessment formula. In this study, many factors were added, but some factors were not included, such as the difference in height between the landslide and the riverbed and the distance between the landslide and the river. These factors are very important in assessments of the risk of a landslide hazard chain, but have seldom been recorded in the existing database. Although the proposed model could fit the database in
Appendix A very well, its applicability still needs to be verified with more detailed documented cases. More factors will be collected and added to the formulas to increase the accuracy of the assessments of the risk of a landslide hazard chain.
6. Conclusions
The hazard chain caused by a landslide is usually more harmful than the landslide itself. In this study, the most common hazard chain, namely landslide–river blockage–outburst flood, was studied.
According to 57 cases of landslide events, the formula for assessing the degree of river blockage was improved, and the accuracy rate in the studied cases reached 96.5%. According to 50 documented historical landslide dams, the discriminant formula of landslide dam stability was put forward, and the accuracy rate in the cases reached 86%. According to 34 cases of dam breaches, a formula for flood peak flow at the dam site was put forward, and the correlation coefficient R2 reached 0.934. A risk assessment model was proposed that combined the peak outflow rate attenuation formula and the improved vulnerability assessment index. The applicability of the proposed model was verified by the Baige landslide, and the calculated results were in good agreement with the measured results.
In practical applications, the advantage of this model is that it can be used to improve the efficiency of the emergency evacuation once a landslide has happened, as well as for increasing the level of engineering security in advance. However, the model has some limitations. Because of the lack of the corresponding measured data, factors such as the water content of landslide material, the river’s velocity, rainfall, and geological action after the landslide’s occurrence cannot be taken into account. The model should be improved regarding these aspects in subsequent research.