Damage Model and the Inﬂuence Factors of Mitigation Engineering against Glacial Debris Flow in the Parlung River Basin, SE Tibetan Plateau

: Understanding the damage mechanism of glacial debris ﬂow mitigation systems is crucial for the risk prevention and assessment of the 200 km trafﬁc corridor in the Southeast of the Tibetan Plateau, where the Sichuan-Tibet railway and expressway have been planned. Based on the phenomena, position and residual efﬁciency of damaged engineering, our analysis of satellite imagery and ﬁeld investigation in multitype spatial reveals the damage types and inﬂuencing factors of glacial debris ﬂow mitigation engineering. An evaluation model which can be used to estimate the engineering damage grade is established by using the relationship between mono engineering works and mitigation systems. In the new model, the engineering damage is divided into ﬁve grades: undamaged, slightly damaged, relatively damaged, seriously damaged, and totally damaged. For glacial debris ﬂow in the Parlung river basin, the ﬁve grades of damage of mitigation works account for 8.70%, 34.78%, 21.74%, 13.04% and 21.74%, respectively. Furthermore, the soil source type and channel proﬁle gradient are the key factors in engineering damage. Design defect of proﬁle gradient is the controlling factor of damage in drainage channel engineering. Based on those results, an engineering damage model is established, which can provide an important reference for risk reduction and prevention of hazards due to the increasing development of trafﬁc engineering.


Introduction
The area of the SE Tibetan Plateau contains the largest number of marine glaciers in China [1].With the confirmation of global warming [2], these glaciers are experiencing significant shrinkage, which has increased year on year [3].This has resulted in abundant numbers of large-scale glacial debris flows in the surrounding area [4].These glacial debris flows can be classified according to the multi-disaster chain by their uncertain initiation mechanism, their various types, and their complex dynamic processes.Since the 1950s, there have been several famous debris flow hazards occurring in the Parlung River (PLR), which have done a great deal of damage to the river and ecosystem.These geohazards dammed the river, such as in the debris flows occurring in Guxiang gully [5] and Peilong gully [6].Geohazards also destroyed the highway and villages, as shown in the debris flows that occurred in Dongru gully [7], Midui gully [8] and Tianmo gully [9].Over the past 60 years, the Chinese government has carried out debris flow mitigation projects continuously to protect the unobstructed parts of the Sichuan-Tibet highway and the villages along the traffic corridor.Along the section of G318 Sichuan Tibet highway from Ranwu town to Tongmai town, there were 23 glacial debris flows that have been controlled by mitigation engineering, accounting for 35.9% of the total number of 64 glacial debris flows that directly harmed the highway.The mitigation system options comprise a drainage channel, check dam, aqueduct, and retention basin.Among them, drainage channel engineering is the most widely used.
Although the evaluation of hazard mitigation engineering has been a popular topic in recent years, the research has not been carried out systematically.The damage grade of engineering has two criteria: Firstly, the damage to the engineering structure, such as fracture, displacement and erosion of the concrete structure.Second is the damage to future mitigation ability, such as the drainage channel losing its diversion functionality, the check dam losing the function of its reservoir, and the reduction of discharge.A lot of studies have been carried out on debris flow mitigation engineering in the southwest mountainous areas of China, especially after the Wenchuan earthquake [10].These works include the evaluation method for determining the benefits of mitigation [11][12][13], the damage types of check dam and drainage channel [14][15][16][17], the laws of the characteristic dynamic changes due to mitigation engineering [18][19][20][21] and the influencing mechanics of earthquake damage engineering [22].In addition to field investigation and experimental research, some numerical methods have been used to real the interaction of debris flow and structure.These works include the Discrete Element Method [23], the Finite Element Method [24]), the SPH method [25], the coupled SPH-FEM [26][27][28], and the Arbitrary Lagrangian-Eulerian method [29].These results have great significance in guiding the disaster mitigation of rainfall-triggered debris flow (RDF).As these studies are mainly involved in RDF, there is a lack of work on the systematic investigation and analysis of glacial debris flow mitigation.
In general, the research on the mitigation of glacier debris flow is insufficient in both theory and practice.Specifically, there is no effective calculation model for the dynamic parameters, such as density, discharge, velocity, total volume, et cetera.Additionally, the interaction mechanisms of engineering structures in response to complex operation conditions are unclear.These ongoing scientific problems lead to glacial debris flow mitigation systems existing in a repeating cycle of building and destruction, followed by rebuilding.For example, the mitigation engineering in Bitong gully and Guxiang gully need to be repaired after every debris flow.This cycle often leads to damage to the highway and causes interruptions to traffic.According to the construction plans of The Sichuan-Tibet Railway, this project passes through the PLR, which should generate impetus for the study and implementation of the glacial debris flow mitigation [30].Therefore, there is a need to summarize the operational experience of glacier debris flow mitigation engineering and identify the key scientific problems in hazard prevention so as to cope with the changing situation of regional glacial debris flow disaster events.Based on the evaluation of the methods by which mitigation engineering resists RDF, this study, through the continuous investigation of glacial debris flow mitigation systems in the PLR, looks at the damage types, establishes a classification method for the degree of engineering damage, explores the key factors affecting engineering damage under initial debris flow conditions and reveals the damage mechanisms addressed by glacial debris flow mitigation engineering.Thus, the paper puts forward some practical suggestions for the improvement of glacial debris flow hazard mitigation.

Formation Conditions of Glacial Debris Flow
The Parlung River basin, crossed by the Sichuan-Tibet railway, is one of the most active regions experiencing glacial debris flow, where debris flow mitigation engineering plays a critical role in the reduction and prevention of damage resulting from glacial geohazards.
The PLR basin is located southeast of the Tibetan Plateau, near the great bend of the Yarlung Zangbo River (Figure 1).In this study area, the geohazard, which is represented by glacial debris flow, is caused by the internal and external dynamics arising from the uplift of the Qinghai-Tibet Plateau [31].
passages, like the Galongqu gully, are distributed in the lower elevation pass of Gangrigabu Mountain.
Influenced by deglaciation and snow avalanche, glacial debris flow can be characterized by fast velocity, large peak discharge, huge total volume, low activity frequency, and high boulder content.These dynamic characteristics, which are different from the characteristics of RDF, present a significant challenge to the function of traditional mitigation engineering.When traditional engineering encounters glacial debris flow, it will be subjected to larger impact forces, much higher discharge and greater scouring effects than specified in the codes to which they have been designed.

Glacial Debris Flow Mitigation
About half of the debris flow gullies with mitigation engineering are in the Ranwu-Yupu section of PLR (Figure 1b), which is only 1/3 of the total length of the river.The main reason is that the debris flow fan is lost due to erosion by the PLR, such that the highway can hardly avoid glacial debris flow disaster.Due to the large flat fans of the others gullies, the highway usually crosses the debris flow gullies with bridges in less hazardous areas.The main purpose of the mitigation design of glacial debris flow in the PLR is to ensure the safety of the Sichuan-Tibet highway, so it primarily considers drainage, along with check dams for disaster dynamic adjustment.Among these engineering projects, drainage channel work is performed in 95.6% of debris flow gullies, with only the Wulian gully not having a designed drainage channel but relying on the retention basin with a check dam.The climate of PLR is mainly affected by the passage of moisture along the Yarlung Zangbo River [32,33].As the summer monsoon from the Indian Ocean enters the valley of PLR, it produces the rainfall gradients in this basin, with a decreasing trend from the lower reaches to the upper reaches [34].At the same time, several branching moisture passages, like the Galongqu gully, are distributed in the lower elevation pass of Gangrigabu Mountain.
Influenced by deglaciation and snow avalanche, glacial debris flow can be characterized by fast velocity, large peak discharge, huge total volume, low activity frequency, and high boulder content.These dynamic characteristics, which are different from the characteristics of RDF, present a significant challenge to the function of traditional mitigation engineering.When traditional engineering encounters glacial debris flow, it will be subjected to larger impact forces, much higher discharge and greater scouring effects than specified in the codes to which they have been designed.

Glacial Debris Flow Mitigation
About half of the debris flow gullies with mitigation engineering are in the Ranwu-Yupu section of PLR (Figure 1b), which is only 1/3 of the total length of the river.The main reason is that the debris flow fan is lost due to erosion by the PLR, such that the highway can hardly avoid glacial debris flow disaster.Due to the large flat fans of the others gullies, the highway usually crosses the debris flow gullies with bridges in less hazardous areas.The main purpose of the mitigation design of glacial debris flow in the PLR is to ensure the safety of the Sichuan-Tibet highway, so it primarily considers drainage, along with check dams for disaster dynamic adjustment.Among these engineering projects, drainage channel work is performed in 95.6% of debris flow gullies, with only the Wulian gully not having a designed drainage channel but relying on the retention basin with a check dam.The second-largest number of engineering projects are check dams.There are nine gullies with 21 check dams, accounting for 39.1% of the total.
These mitigation engineering projects have played a significant role in mitigating glacial debris flow hazards along the Sichuan-Tibet highway.In addition, they have also had benefits in ecological protection and landscape improvement throughout the study area.Even with the trend of increasing frequency of debris flow hazards, those construction works still have an important benefit for the safety of the highway.

Evaluation Method of Engineering Damage
Although the evaluation of hazard mitigation engineering has been a popular topic in recent years, the research has not been carried out systematically.The damage grade of engineering has two criteria: Firstly, the damage to the engineering structure, such as fracture, displacement and erosion of the concrete structure.Second is the damage to future mitigation ability, such as the drainage channel losing its diversion functionality, the check dam losing the function of its reservoir, and the reduction of discharge.This paper refers to the earthquake damage evaluation method [22,35] and builds a new model based on the damage grade of structural reliability and functionality of mono engineering works, as shown in Table 1.In this new model, damage grade is divided into five levels from I to V. Encompassing the association and function combination relationship between mono engineering works and mitigation systems, the classification method is established, as shown in Figure 2. The second-largest number of engineering projects are check dams.There are nine gullies with 21 check dams, accounting for 39.1% of the total.These mitigation engineering projects have played a significant role in mitigating glacial debris flow hazards along the Sichuan-Tibet highway.In addition, they have also had benefits in ecological protection and landscape improvement throughout the study area.Even with the trend of increasing frequency of debris flow hazards, those construction works still have an important benefit for the safety of the highway.

Evaluation Method of Engineering Damage
Although the evaluation of hazard mitigation engineering has been a popular topic in recent years, the research has not been carried out systematically.The damage grade of engineering has two criteria: Firstly, the damage to the engineering structure, such as fracture, displacement and erosion of the concrete structure.Second is the damage to future mitigation ability, such as the drainage channel losing its diversion functionality, the check dam losing the function of its reservoir, and the reduction of discharge.This paper refers to the earthquake damage evaluation method [22,35] and builds a new model based on the damage grade of structural reliability and functionality of mono engineering works, as shown in Table 1.In this new model, damage grade is divided into five levels from I to V. Encompassing the association and function combination relationship between mono engineering works and mitigation systems, the classification method is established, as shown in Figure 2.

Data of Mitigation Engineering
Almost all drainage channels are damaged after debris flow events.Through the field investigation of mitigation engineering in the Ranwu-Tongmai section of PLR from 2016 to 2020, the data of design standard, dynamic characteristics, engineering location, implementation period and treatment method are obtained, and the debris flow hazard characteristics are collected, shown in Table 2.

Damage Type of Mitigation Engineering
The statistical results show that there are five damage types resulting from both dynamic characteristics and design defects, shown in Table 3.Among those five types, the main damage type is sedimentary damage, which occurs in 16 engineering works of 14 debris flow gullies, followed by scour damage, which occurs in 16 engineering works of 12 debris flow gullies.Impact damage is the least common, only occurring in one engineering work.At the same time, the damage caused by design defects is very significant, proving that there are still great deficiencies in the design method of glacial debris flow mitigation in China.We summarize the engineering damage to drainage channels with four classifications: abrasion damage, scour damage, sedimentary damage, and brim-over damage.Abrasion damage and scour damage are the most common failure forms in the study area, such as erosion of the engineering surface (Figure 3a,b), scour pit in the channel bed (Figure 3d), and side-wall collapse (Figure 3c).Sedimentary damage mainly occurs at the section where the channel gradient changes and most typically occurs in the channel under a bridge (Figure 3e,f).Brim-over damage is mainly caused by design deficiency when the discharge is greater than the drainage capacity, and giant rocks contained in the fluid are the key factor.It usually occurs at the entrance section and bend section of the drainage channel.
Based on the damage location, as with the drainage channel, the damage of check dams can also be divided into four types: dam abutment damage, dam body damage, dam foundation damage, and auxiliary structure damage, shown in Figure 4.
Dam abutment failure is caused by seepage and scour, which erode the soil and loosen the abutment foundation (Figure 4a,b).Dam body damage is not as common as abutment damage in the study area, which is mainly characterized by cracks caused by impact force (Figure 4c) and notches caused by abrasion (Figure 4d).Dam foundation damage is the most common failure mode.In general, all check dams would be damaged with varying degrees after experiencing debris flow.However, when the scouring depth is less than the foundation depth, its impact on dam stability is limited.Auxiliary structure damage usually includes the erosion and deposition of the overflow structure (Figure 4e), which affects the reliability and function of the check dam.In addition, the damage to the anti-erosion structure is widely distributed (Figure 4f).

Evaluation of Mono Engineering Work Damage
Applying the damage evaluation model established in Table 1 and Figure 2 to the mitigation engineering in PLR, the glacial debris flow drainage channel is statistically analyzed.The results are shown in Figure 5.There are two undamaged engineering works (I grade damage), seven slightly damaged engineering works (II grade damage), three relatively damaged engineering works (III grade damage), five seriously damaged engineering works (IV grade damage), and five totally damaged engineering works (V grade damage).The percentages of these damage grades were 9.09%, 31.82%,13.63%, 22.73%, and 22.73%, respectively.The evaluation results of the check dam damage are shown in Table 4.There are two undamaged engineering works (I grade damage), four slightly damaged engineering works (II grade damage), six relatively damaged engineering works (III grade damage), five seriously damaged engineering works (IV grade damage), and four totally damaged engineering works (V grade damage).The percentages of these damage grades were 9.53%, 19.05%, 28.57%, 23.81%, and 19.05%, respectively.These damaged engineering works indicate that check dams are not suitable for low-frequency glacial debris flow.

Evaluation of Mono Engineering Work Damage
Applying the damage evaluation model established in Table 1 and Figure 2 to the mitigation engineering in PLR, the glacial debris flow drainage channel is statistically analyzed.The results are shown in Figure 5.There are two undamaged engineering works (I grade damage), seven slightly damaged engineering works (II grade damage), three relatively damaged engineering works (III grade damage), five seriously damaged engineering works (IV grade damage), and five totally damaged engineering works (V grade damage).The percentages of these damage grades were 9.09%, 31.82%,13.63%, 22.73%, and 22.73%, respectively.The evaluation results of the check dam damage are shown in Table 4.There are two undamaged engineering works (Ⅰ grade damage), four slightly damaged engineering works (II grade damage), six relatively damaged engineering works (III grade damage), five seriously damaged engineering works (IV grade damage), and four totally damaged engineering works (V grade damage).The percentages of these damage grades were

Evaluation of Mitigation System Damage
We have shown the evaluation results of mitigation systems in Figure 2 in Table 5.Among the 23 mitigated debris flow gullies in the PLR basin, there are two undamaged mitigation systems (I grade damage), eight slightly damaged mitigation systems (II grade damage), five relatively damaged mitigation systems (III grade damage), three seriously damaged mitigation systems (IV grade damage), and five totally damaged mitigation systems (V grade damage).The percentages of these damage grades were 8.70%, 34.78%, 21.74%, 13.04%, and 21.74%, respectively.Table 5. Damage evaluation of the mitigation systems of single debris flow gullies.(The legend of damage grade is the same as Figure 5).The force causing the engineering damage shows a pulse distribution on a time scale, with the force position and vector direction also being random.This makes it difficult to quantitatively analyze the mechanism of structural damage.Besides, the engineering damage has the characteristics of multi-pathogenesis and the same symptoms, which also makes the analysis of the mechanical processes difficult.

Damage
(1) Influence of debris flow formation condition on drainage channel damage In this paper, we summarize the formation conditions of glacial debris flow into four categories: basin area, channel gradient, glacier area and soil source type.We analyze the relationship between engineering damage and formation conditions, as shown in Figures 6 and 7.
The results show that the channel gradient has a positive correlation with the damage grade and has the biggest impact on sedimentary damage.When the soil sources of debris flow are moraine (M) plus landslide (L) or snow avalanche (A) plus slope erosion (E), it usually results in a seriously damaged grade and a sedimentary damage type.Such types of glacial debris flow should be paid significant attention to when designing and constructing hazard mitigation projects.Broadly speaking, the watershed area has a twolevel differentiation effect on the damage grade and damage type.The large area and small area are both closely related to the damage grade, although the exact relationship is unclear.The impact of the glacier area is consistent with the watershed area due to the hydrodynamic conditions of glacial debris flow, mainly rainfall and glacier melt, which are closely related to the area parameter.For the damage type, when the watershed area is larger than 10 km 2 , there is mainly scour damage and abrasion damage, while when less than 10 km 2 , there is mainly brim-over damage and sedimentary damage.
(2) Influence of design defects on drainage channel damage According to previous damage analysis, engineering design defects mainly lead to scour damage of the overflow structure of check dams, caused by a lack of efficiency of flow; it could also cascade by triggering damage to the dam abutment and dam foundation.Through field investigation, the two flow areas of debris flow events, both in the reservoir entrance and the check dam, are compared with the designed area of the overflow structure, as shown in Figure 8.It can be seen from the figure that the most serious abutment damage is to the dam of 22-2 (No.2 dam in Guxiang gully), also shown in Figure 4b and Table 3.However, dam 15-3 (No.3 dam in Xinkongdong gully) experienced no damage because the reservoir of the check dam was empty, and the debris flow was blocked by the reservoir.The check dam greatly reduces the peak discharge so that the overflow structure can meet the hazard requirements.As the reservoir is now full of soil materials, if it experiences a similar debris flow in the future, the overflow capacity will be utterly insufficient.Such an event will trigger damage to the overflow structure when the overflow capacity is less than the size of the debris flow hazard.This kind of issue usually arises from the discharge calculating model of glacial debris flow being inaccurate.The results show that the channel gradient has a positive correlation with the damage grade and has the biggest impact on sedimentary damage.When the soil sources of debris by the reservoir.The check dam greatly reduces the peak discharge so that the overflow structure can meet the hazard requirements.As the reservoir is now full of soil materials, if it experiences a similar debris flow in the future, the overflow capacity will be utterly insufficient.Such an event will trigger damage to the overflow structure when the overflow capacity is less than the size of the debris flow hazard.This kind of issue usually arises from the discharge calculating model of glacial debris flow being inaccurate.The most common defect of drainage channel design is the channel gradient design, and the engineering work that is built is often designed according to the terrain of the debris flow fan.This leads to several types of damage to the channel engineering work.The drainage channel experiences sedimentary damage when the gradient is less than 200.00‰; otherwise, it experiences abrasion damage and scour damage (Figure 9).There is an outlier in the form of No.21 in Figure 9, which is the drainage channel in Guxiang gully, which experiences abrasion damage when the gradient is less than 200.00‰.The reason for this is that there are two check dams upstream of the drainage channel.After most of the soil materials of debris flow are deposited in the reservoir, the debris flow changes to a torrential flood, which will hardly deposit any material but will produce strong abrasion and scour forces under the higher velocity.The most common defect of drainage channel design is the channel gradient design, and the engineering work that is built is often designed according to the terrain of the debris flow fan.This leads to several types of damage to the channel engineering work.The drainage channel experiences sedimentary damage when the gradient is less than 200.00‰; otherwise, it experiences abrasion damage and scour damage (Figure 9).There is an outlier in the form of No.21 in Figure 9, which is the drainage channel in Guxiang gully, which experiences abrasion damage when the gradient is less than 200.00‰.The reason for this is that there are two check dams upstream of the drainage channel.After most of the soil materials of debris flow are deposited in the reservoir, the debris flow changes to a torrential flood, which will hardly deposit any material but will produce strong abrasion and scour forces under the higher velocity.The gradient change of a drainage channel is the most important factor for sedimentary damage.The ratio of damage length is related to the initial gradient of the channel, the amplitude of gradient change, and the length ratio between two sections of the channel with different gradients (Figure 10).On the whole, the length ratio and amplitude of gradient change are inversely proportional to the sedimentary damage grade, while the initial gradient has less effect on the sedimentary damage grade.Based on the four influence factors of the ratio of damage length, through mathematical function fitting, the formula of sedimentary damage length is established as in Equa- The gradient change of a drainage channel is the most important factor for sedimentary damage.The ratio of damage length is related to the initial gradient of the channel, the amplitude of gradient change, and the length ratio between two sections of the channel with different gradients (Figure 10).On the whole, the length ratio and amplitude of gradient change are inversely proportional to the sedimentary damage grade, while the initial gradient has less effect on the sedimentary damage grade.
Based on the four influence factors of the ratio of damage length, through mathematical function fitting, the formula of sedimentary damage length is established as in Equation ( 1) where the Y is sedimentary damage length; L 0 is the initial length divided by the gradient change length; η is the amplitude of gradient change, itself calculated by J 0 /J 1 , where J 0 is the initial gradient, and J 1 is the changed gradient; e is the natural number.
The test of Equation ( 1) is shown in Figure 11.The test result shows that there is a good correlation between the fitting value of Equation ( 1) and the measured value, such  The gradient change of a drainage channel is the most important factor for sedimentary damage.The ratio of damage length is related to the initial gradient of the channel, the amplitude of gradient change, and the length ratio between two sections of the channel with different gradients (Figure 10).On the whole, the length ratio and amplitude of gradient change are inversely proportional to the sedimentary damage grade, while the initial gradient has less effect on the sedimentary damage grade.
where the Y is sedimentary damage length; L0 is the initial length divided by the gradient change length; η is the amplitude of gradient change, itself calculated by J0/J1, where J0 is the initial gradient, and J1 is the changed gradient; e is the natural number.The test of Equation ( 1) is shown in Figure 11.The test result shows that there is a good correlation between the fitting value of Equation ( 1) and the measured value, such that the empirical formula is suitable for the study of drainage channels subject to glacial debris flow in the PLR basin.

Analysis of Engineering Damage Model
Based on the above analysis, the damage models of the drainage channel and check dam can be established, shown in Figures 12 and 13, respectively.
The four factors of the debris flow formation conditions influence the disaster characteristic of discharge, velocity, density, and boulder content, causing brim-over damage, scour damage, sedimentary damage, and abrasion damage, respectively.In engineering design, the design of horizontal bend and gradient should consider sedimentary damage on the concave bank.The design of energy dissipation engineering should consider the scour and abrasion damage caused by boulders.The damage to the check dam is more complex, such that we cannot obtain the force mechanism based on our investigations.For glacial debris flow hazard mitigation, check dams are commonly used for auxiliary engineering in the mitigation system because of the restraint of the protection object and the terrain of the debris fan.The key factors in check dam damage are the interaction between engineering and debris flow, designing of the overflow structure and designing of the anti-erosion structure (Figure 13).None of these problems are thoroughly studied in the mitigation of rainfall-triggered debris flow.

Analysis of Engineering Damage Model
Based on the above analysis, the damage models of the drainage channel and check dam can be established, shown in Figures 12 and 13, respectively.
The four factors of the debris flow formation conditions influence the disaster characteristic of discharge, velocity, density, and boulder content, causing brim-over damage, scour damage, sedimentary damage, and abrasion damage, respectively.In engineering design, the design of horizontal bend and gradient should consider sedimentary damage on the concave bank.The design of energy dissipation engineering should consider the scour and abrasion damage caused by boulders.The damage to the check dam is more complex, such that we cannot obtain the force mechanism based on our investigations.For glacial debris flow hazard mitigation, check dams are commonly used for auxiliary engineering in the mitigation system because of the restraint of the protection object and the terrain of the debris fan.The key factors in check dam damage are the interaction between engineering and debris flow, designing of the overflow structure and designing of the anti-erosion structure (Figure 13).None of these problems are thoroughly studied in the mitigation of rainfall-triggered debris flow.
complex, such that we cannot obtain the force mechanism based on our investigations.For glacial debris flow hazard mitigation, check dams are commonly used for auxiliary engineering in the mitigation system because of the restraint of the protection object and the terrain of the debris fan.The key factors in check dam damage are the interaction between engineering and debris flow, designing of the overflow structure and designing of the anti-erosion structure (Figure 13).None of these problems are thoroughly studied in the mitigation of rainfall-triggered debris flow.

Discussion
(1) Complexity of interaction between glacial debris flow and engineering The characteristics of glacial debris flow are quite different from that of RDF.The particle size distribution of the giant block is larger, resulting in greater heterogeneity of physical properties.The wear process and mechanical mechanism of the engineering surface under the alternating action of water flow and debris flow is also an important difference.Under the effect of various types of debris flows, there is no essential difference in the damage phenomenon of mitigation engineering, but some differences consist in the apparent details.Compared with the mitigation engineering of RDF, for the GRDF, the pit abrasion damage and the dam body damage are more serious [16].The high content rate and velocity of the giant rock in GRDF are the triggering factors, which lead to strong impact energy on the engineering.The damage to the dam foundation is relatively slight, which may be caused by the higher strength of the moraine soil [36].The sedimentary damage in the drainage of GRDF is caused by the change in the hydraulic gradient of the channel bed, which is different from that of RDF due to the fluid properties [15].The damage to the check dam mainly occurs from the abutment part but not the center part, which

Discussion
(1) Complexity of interaction between glacial debris flow and engineering The characteristics of glacial debris flow are quite different from that of RDF.The particle size distribution of the giant block is larger, resulting in greater heterogeneity of physical properties.The wear process and mechanical mechanism of the engineering surface under the alternating action of water flow and debris flow is also an important difference.Under the effect of various types of debris flows, there is no essential difference in the damage phenomenon of mitigation engineering, but some differences consist in the apparent details.Compared with the mitigation engineering of RDF, for the GRDF, the pit abrasion damage and the dam body damage are more serious [16].The high content rate and velocity of the giant rock in GRDF are the triggering factors, which lead to strong impact energy on the engineering.The damage to the dam foundation is relatively slight, which may be caused by the higher strength of the moraine soil [36].The sedimentary damage in the drainage of GRDF is caused by the change in the hydraulic gradient of the channel bed, which is different from that of RDF due to the fluid properties [15].The damage to the check dam mainly occurs from the abutment part but not the center part, which is completely different from the numerical simulation results [13].The reason is that the vulnerability of the joint between soil and engineering is overlooked in the numerical study.The damage characteristics obtained based on the operation of the structure can provide a great reference value for the design of mitigation engineering of GRDF.
(2) Some suggestions on mitigation systems for glacial debris flow When the possibility of blocking the main river by debris flow is slight, the drainage channel solution should be used as far as possible.The drainage channel has the advantages of established technology and clearly understood influence mechanisms.Unlike a check dam, even in the case of partial damage, a drainage channel would rarely trigger a disaster chain.Before there is a new calculation model of the dynamic parameters of glacial debris flow, it is not suitable to use a check dam as the main mitigation engineering solution, especially a gravity check dam.On the other hand, due to the high content of boulders in glacial debris flow, the grid-type dam can be considered to retain boulders in the wider and gentler sections, which would be of great value to the safety of the drainage channel.
(3) Suggestions on the mono-mitigation engineering The analysis of influencing factors reveals that the gradient, bend and boulders have a significant impact on the drainage channel.Therefore, some suggestions are put forward for engineering design.The gradient of the channel should be controlled at 200~300‰.When the drainage channel passes through a bridge, the gradient under the bridge should be greater than that upstream so as to reduce sedimentary damage.When the gradient change cannot be avoided, Equation (1) suggested in this paper can be used to calculate the sedimentary damage length.
When it is necessary to build a check dam, it should be designed as a grid-type dam, and there should be an increase in the size of the overflow structure.In the absence of the calculation method of debris flow overflow, the overflow area can be based on the flow area at the entrance of the reservoir to ensure that the overflow area of the dam is larger than that of the channel.The influence of boulders should be considered in the design of auxiliary engineerings, such as anti-scour engineering.
(4) The Sichuan-Tibet railway and the Sichuan-Tibet Expressway both cross the areas where glacial debris flows are the most common geohazard.Disaster mitigation in project construction creates new requirements for the theoretical research of glacial debris flow.There are two points to be made closely related to mitigation technology.First, it is necessary to clarify the hydrological process of water recharge by glacier and snow avalanche in debris flow processes, such as to understand the influence of ice/snow on dynamic parameters.Second, it is urgent to study the relationship between the initiation mechanism of glacial debris flow and hazard characteristics and establish an evaluation method to determine the extreme conditions of the geohazard.

Conclusions
To evaluate the operational state of the mitigation engineering, which was constructed to mitigate glacial debris flow along the Sichuan-Tibet highway in the PLR basin, a field investigation of the engineering was carried out.Based on the field investigation results, the damage type and influence factors are summarized, and some suggestions for mitigation system design are made for the construction of the Sichuan-Tibet railway.
(1) Five damage types to the mitigation engineering of glacial debris flow are summarized: scour damage, abrasion damage, impact damage, sedimentary damage and damage triggered by insufficient function.Among those five types, the driving forces of the first three types are dynamic characteristics, with design deficiency causing the last two types.
(2) A new evaluation model of mitigation engineering for glacial debris flow is established based on the damage to function and damage to engineering structure.The five damage grades are obtained by the field investigation, which is: undamaged (I grade), slightly damaged (II grade), relatively damaged (III grade), seriously damaged (IV grade), and totally damaged (V grade).
(3) By analyzing the key factors of damage type and damage grade, it is found that the geohazard formation conditions, which are characterized by channel gradient, watershed area and soil source type, are related to the damage to the drainage channel.On this basis, the damage model for engineering works is established.
(4) It is revealed that the LJ gradient for determining the damage type of the drainage channel is 20.00%.When the gradient is below 20.00%, the drainage channel will experience sedimentary damage.Otherwise, it will experience scour damage and abrasion damage.A formula to calculate the length of sedimentary damage is built using the three parameters of the ratio of channel length with different gradients, amplitudes of gradient change, and the initial gradient.

Figure 1 .
Figure 1.Sketch map of the investigated area.

Figure 1 .
Figure 1.Sketch map of the investigated area.

Figure 2 .
Figure 2. The classification method of damage grade for debris flow mitigation engineering.

Figure 2 .
Figure 2. The classification method of damage grade for debris flow mitigation engineering.

Water 2023 , 1 Figure 3 .
Figure 3. Damage types for drainage channels subject to glacial debris flow.

Figure 3 .
Figure 3. Damage types for drainage channels subject to glacial debris flow.

Figure 4 .
Figure 4. Damage types for check dams subject to glacial debris flow.

Figure 4 . 19 Figure 5 .
Figure 4. Damage types for check dams subject to glacial debris flow.Water 2023, 15, x FOR PEER REVIEW 10 of 19

Figure 5 .
Figure 5. Damage grade division of drainage channels in PLR.

Figure 6 .
Figure 6.The relationship between debris flow formation conditions and damage grade.

Figure 7 .
Figure 7.The relationship between debris flow formation conditions and damage type.B: brimover damage, S: scour damage, A: abrasion damage, D: sedimentary damageThe results show that the channel gradient has a positive correlation with the damage grade and has the biggest impact on sedimentary damage.When the soil sources of debris

Figure 6 . 19 Figure 6 .
Figure 6.The relationship between debris flow formation conditions and damage grade.

Figure 7 .
Figure 7.The relationship between debris flow formation conditions and damage type.B: brimover damage, S: scour damage, A: abrasion damage, D: sedimentary damage

Figure 7 .
Figure 7.The relationship between debris flow formation conditions and damage type.B: brim-over damage, S: scour damage, A: abrasion damage, D: sedimentary damage.

Figure 8 .
Figure 8. Comparative analysis of overflow area of typical check dam.

Figure 8 .
Figure 8. Comparative analysis of overflow area of typical check dam.

Water 2023 , 19 Figure 9 .
Figure 9.The relationship between the gradient of the drainage channel and damage type.

Figure 10 .
Figure 10.Relationship between gradient change and the ratio of sedimentary damage length.

Figure 9 .
Figure 9.The relationship between the gradient of the drainage channel and damage type.
formula is suitable for the study of drainage channels subject to glacial debris flow in the PLR basin.

Figure 9 .
Figure 9.The relationship between the gradient of the drainage channel and damage type.

Figure 10 .
Figure 10.Relationship between gradient change and the ratio of sedimentary damage length.Based on the four influence factors of the ratio of damage length, through mathematical function fitting, the formula of sedimentary damage length is established as in Equation(1)

Figure 10 . 19 Figure 11 .
Figure 10.Relationship between gradient change and the ratio of sedimentary damage length.

Figure 12 .
Figure 12.The damage models of drainage channels under glacial debris flow.Figure 12.The damage models of drainage channels under glacial debris flow.

Figure 12 . 19 Figure 13 .
Figure 12.The damage models of drainage channels under glacial debris flow.Figure 12.The damage models of drainage channels under glacial debris flow.Water 2023, 15, x FOR PEER REVIEW 16 of 19

Figure 13 .
Figure 13.The damage models of a check dam under glacial debris flow.

Table 1 .
The classification method of damage grade of engineering structure.

Table 2 .
The investigation data of mitigation engineering of glacial debris flow in the PLR.

Table 3 .
Statistics of damage type of mitigation engineering of glacial debris flow.

Table 4 .
Statistics of damage grade division of check dam.(The legend of damage grade is the same as Figure5).