Recent Advances in Stability and Failure Mechanisms of Landslide Dams

Numerous landslide dams have been induced in recent years as a result of frequent earthquakes and extreme climate hazards. Landslide dams present serious threats to lives and properties downstream due to potentially breaching floods from the impounded lakes. To investigate the factors influencing the stability of landslide dams, a large database has been established based on an in-depth investigation of 1,737 landslide dam cases. The effects of triggers, dam materials, and geomorphic characteristics of landslide dams on dam stability are comprehensively analyzed. Various evaluation indexes of landslide dam stability are assessed based on this database, and stability evaluation can be further improved by considering the dam materials. Stability analyses of aftershocks, surges, and artificial engineering measures on landslide dams are summarized. Overtopping and seepage failures are the most common failure modes of landslide dams. The failure processes and mechanisms of landslide dams caused by overtopping and seepage are reviewed from the perspective of model experiments and numerical analyses. Finally, the research gaps are highlighted, and pathways to achieve a more complete understanding of landslide dam stability are suggested. This comprehensive review of the recent advances in stability and failure mechanisms of landslide dams can serve as a key reference for stability prediction and emergency risk mitigation.


INTRODUCTION
Landslide dams are natural dams formed by river blockages with massive amounts of materials from avalanches, landslides, or debris flows (Swanson et al., 1986;Capra, 2006;Korup and Tweed, 2007;Hermanns, 2013;Shi et al., 2014). A large number of landslide dams have been induced by frequent earthquakes (Huang and Fan, 2013;Fan et al., 2014), extreme climate hazards Gariano and Guzzetti, 2016), or snowmelt (Costa and Schuster, 1991;Strom, 2010;Frey et al., 2018). For instance, multiple landslide dams were clustered along short river reaches induced by the 1999 Chi-Chi earthquake (Liao and Lee, 2000), the 2004 Chuetsu earthquake (Wang et al., 2007), and the 2008 Wenchuan earthquake (Fan et al., 2012). Eighteen landslide dams were formed by the 2009 Morakot typhoon in Taiwan (Chen and Chang, 2016). Landslide dams have frequently occurred in the rivers on the Tibetan Plateau, caused by slides and debris flows induced by snowmelt water or active glaciers (Fan et al., 2019;Wang et al., 2020).
Landslide dams present serious threats to lives and properties downstream from a potentially rapid release of the impounded water if the dams fail (Huang and Fan, 2013;Fan et al., 2014). This threat is greater for some landslide dams with large volumes: if loss of stability occurs, the breached floods could be hundreds of times the maximal recorded floods in history, and could destroy almost everything in their paths (Schuster, 2000;Hegan et al., 2001;Yin et al., 2009). The evaluation of landslide dam stability is thus crucial before taking any effective mitigation measures.
The stability of a landslide dam is usually difficult to determine. The formation and instability of landslide dams are regulated by the complex interactions between the geomorphological parameters of the valley (geometries of landslide dam and lake), the slide characteristics (dam material), and the hydrological parameters of the river (inflow rate and lake volume). Moreover, landslide dam stability is affected by aftershocks, surge actions, and human control measures. Korup (2002) reviews the geomorphic assessment indexes of landslide dam stability and points out that there is still a considerable lack of understanding regarding the stability and failure mechanisms of landslide dams. Since then, landslide dam cases on a regional and global scale have been widely collected to obtain their spatial morphology and distribution regularity (Fan et al., 2020). Based on these field cases, the effects of the characteristic parameters on landslide dam stability have been explored Evans et al., 2011). In addition, evaluation models of landslide dam stability have been established by means of statistical analysis (Dong et al., 2009;Shen et al., 2020). Subsequently, the failure mechanisms of landslide dams caused by overtopping and seepage have been explored using model experiments and numerical analyses (Shi et al., 2015a,b,c;Zhou G.G.D. et al., 2019;. However, there are few comprehensive reviews of the stability and failure mechanisms of landslide dams.
In this study, a large database with 1,737 landslide dam cases worldwide is compiled to give a comprehensive perspective of landslide dam stability. Based on this database, the effects of the characteristic parameters on landslide dam stability are first reviewed in Section "Influences of Characteristic Parameters on Landslide Dam Stability." Stability analyses of aftershocks, surges, and artificial engineering measures on the landslide dams are summarized in Section "Stability Analyses of Special Factors on Landslide Dams." Then, detailed investigations of failure mechanisms of landslide dams caused by overtopping and seepage are reviewed in Sections "Overtopping Failure Analyses of Landslide Dams" and "Seepage Failure Analyses of Landslide Dams, " respectively. Finally, research gaps and pathways are highlighted in Section "Future Research Directions." An exhaustive review of the recent advances in stability and failure mechanisms of landslide dams can provide a reliable basis for stability prediction, early warning, and emergency risk mitigation in mountainous areas.

INFLUENCES OF CHARACTERISTIC PARAMETERS ON LANDSLIDE DAM STABILITY
In this section, the effects of single and multiple characteristic parameters on landslide dam stability are reviewed. A database containing 1,737 landslide dams located in 45 countries and regions has been assembled over a decade through literature review (Shen et al., 2020), to provide detailed insight into the stability of landslide dams.
The distribution of landslide dams is closely related to seismic zones worldwide as shown in Figure 1. Landslide dams are widely distributed in China, Italy, Japan, and the United States, due to a wide range of mountainous regions and numerous seismic belts. By contrast, few landslide dams documented in the literature occur in Africa.
Due to uncertainties in the life periods and failure modes of landslide dams, their stability can be difficult to determine. A landslide dam might have existed for several decades or even hundreds of years, but its ultimate failure may be caused by a heavy rainstorm or other factors. In this paper, a landslide dam is called stable when a backwater lake still exists or is filled in with sediments at the time of its analysis (Fan et al., 2020). Those classified as "unstable" have encountered a catastrophic breach with the abrupt release of the impounded lake waters or piping having developed from upstream to downstream of the dam.
Landslide dam failure could be induced by overtopping, piping or downstream slope slide Fan et al., 2020). Overtopping failures appear where landslide dams undergo erosion or collapse leading to a catastrophic breach, with the subsequent release of the impounded lake waters (Figures 2A-C). Piping gradually develops when small particles are transported to free exits or into coarse openings driven by hydraulic gradient and lake waters are released through the leakage channel ( Figure 2D). Downstream slope failures occur where the buoyant weight of the downstream dam is lower than the uplift force with a rapid increase in the upstream water level.
Overtopping failure is the most common instability mode, accounting for approximately 90% of the inventoried landslide dams which fail (Figure 3), whereas seepage failure is less common. The failure situation is significantly different from that of embankment dams (Foster et al., 2000;Richards and Reddy, 2007). The width of landslide dams can be hundreds of meters or even thousands of meters, and the ratio of the dam width to the height is small, contributing to a low hydraulic gradient from the dam upstream to downstream sides (Shen et al., 2020). Therefore, seepage failures do not often develop in landslide dams. Nevertheless, seepage has a significant impact on the process of overtopping failure due to the hydraulic gradient from dam upstream to downstream (Gregoretti et al., 2010;.

Single-Factor Analyses of Landslide Dam Stability
Landslide dam stability is influenced by various parameters. According to case analyses in the literature (Costa and Schuster, 1988;Casagli et al., 2003), the triggers of slope failures and accumulated materials affect landslide dam stability. The geomorphic parameters characterizing both the dam and the blocked river channel also have an influence on the dam stability (Pirocchi, 1992;Casagli and Ermini, 1999). Based on the compiled landslide dam database, the effects of a single parameter on landslide dam stability are discussed below.

Landslide Dam Triggers
Landslide dams are triggered by earthquakes, rainfall, snowmelt, volcanic eruptions, and other factors . Most landslide dams are caused by earthquakes and rainfall; 50.4% of landslide dams (700 cases) are induced by earthquakes and 39.3% of landslide dams (546 cases) are caused by rainfall, as shown in Figure 4A.
In general, 65% of landslide dams cannot remain stable as shown in Figure 4B. Compared with landslide dams induced by earthquakes and other factors, a higher proportion of dams induced by rainfall and snowmelt are unstable: 67.4 and 66.7%, respectively. This trend is consistent with the landslide dams inventoried by Stefanelli et al. (2015) that indicate dam failures caused by intense rainfall and snowmelt make up 53 and 60% of the failures, respectively. It is speculated that the inflow rate of the backwater lake is larger when accompanied by rainfall and snowmelt, resulting in a higher flow velocity and erosion rate during the breach process. Moreover, the water content of landslide dams is high and the material strength of the wet dams is lower than that of relatively dry dams (Nian et al., 2018). Both of these factors could reduce the stability of a landslide dam.

Landslide Dam Material
The material composition of landslide dams could be divided into soil, debris, and rock . As shown in Figure 5A, the proportion of soil-type landslide dams among the inventoried landslide dams is the smallest (21%). The proportion of soiltype landslide dams that are unstable is higher than for the rock and debris types, reaching nearly 80%. Rock-type landslide dams maintain stability better than soil-type and debris-type dams.
Rock-type landslide dams are induced by concordant, discordant, or oblique failures at bedding, jointing, and other contact or bedrock structures (Cui et al., 2009). The Xiaojiaqiao landslide dam was a typical case of rock-type, which had a sliding surface along a bedding-parallel fault with a slickenside and a gouge with a thickness of a few centimeters in carbonate rock (Chigira et al., 2010). As a result of high material strength and stability, blasting was explored to reduce the potential risk of a backwater lake.
Debris-type landslide dams, which are the most prevalent among the inventoried dams (Figure 5A), are a mixture of soil, rock, coarse grains, and cuttings. These dams either consist of a relatively homogenous mixture of rock and soil or a distribution of these two materials in different zones (Shen et al., 2020). Taking the Tangjiashan landslide dam as an example, the soil was concentrated on the top of the dam and rock was on the  (Quenta et al., 2007). Overtopping failure occurred for Jure, Tangjiashan, and Baige landslide dams, whereas seepage failure occurred for Allpacoma landslide dam. bottom. The breaching process ceased when erosion induced by an outburst flooded into the rock layer .
Soil-type landslide dams are mostly caused by the failure of shallow slopes consisting of colluvium deposit, weathered sediment, loess, man-made fill, or materials displaced from previous man-made construction (Cui et al., 2009). Generally, this type of landslide dam has a low material strength and high erodibility. Taking the Hsiaolin landslide dam as an example, this landslide dam lasted less than 1 h, and it only took approximately 8 min to completely breach due to its low soil strength . This situation also arose for the New-Street soil-type landslide dam (Cui et al., 2009).

Geomorphic Characteristics
The geomorphic characteristics of landslide dams could be categorized as geometric and hydrological parameters. Peng and Zhang (2012) define these parameters as shown in Table 1. The geometric parameters are summarized as the dam height, dam length, dam width, and dam volume. The hydrological parameters are listed as the inflow rate and lake volume.
It is commonplace that the geometric parameters of landslide dams cannot be obtained immediately with field investigations due to their remote locations. Fan et al. (2012) employed postlandslide remote sensing images and the pre-landslide digital terrain model in the geographic information system (GIS) to estimate the geomorphic characteristics. This method has also been adopted by Kuo et al. (2011) and Dong et al. (2014). While the dam height, dam length, and lake volume could be obtained by GIS, the dam width in the upstream direction is difficult to determine from a satellite image because of the water stored in the backwater lake .

Dam width and dam height
Landslide dam stability increases with increasing dam width and the ratio of the dam width to the dam height as shown in Figures 5B,C. The proportion of failed landslide dams with dam widths smaller than 250 m is 88%, compared to 48% for landslide dams with dam widths larger than 500 m. The proportion of failed landslide dams for which the ratio of the dam width to the dam height is smaller than 15 is 87%, which reduces to 65% for landslide dams with a ratio of the dam width to the dam height larger than 30.
From a physical point of view, the dam height and dam width are important parameters to assess landslide dam stability against both overtopping and seepage failures. For the former, they determine the steepness of the dam slope downstream and the corresponding flow velocity and erosion ability of the overtopping waters. For the latter, they control the elevation of the water table through the landslide dam and especially its hydraulic gradient .

Dam volume
The global dam stability increases with increasing landslide dam volume (Figure 6). The proportions of failed landslide dams  (Peng and Zhang, 2012  with dam volumes smaller than 4 × 10 6 m 3 and larger than 20 × 10 6 m 3 are 82 and 48%, respectively. The volume could well identify the dam stability by affecting the dam geometry: width, height, and length. Dam volume is always considered as an influencing factor in the assessment indexes of landslide dam stability (Canuti et al., 1998;Ermini and Casagli, 2003). Evans et al. (2011) present a special review on natural and artificial rockslide dams and conclude that most recorded dams exceed 1 million m 3 in volume. This also applies to our inventoried landslide dams (Figure 6). Whether the landslide volume or landslide dam volume should be selected as the factor influencing landslide dam stability remains controversial. According to Stefanelli et al. (2016), the morphometric data processing of post-landslide dams leads to some errors, proportional to the amount of eroded material. The loss of relative volume and the percentage error are smaller if compared to the total landslide volume rather than to the dam volume.

Lake volume
There is no discernable relationship between landslide dam stability and lake volume ( Figure 6B). However, landslide dam stability decreases gradually as the ratio of the lake volume to the dam volume increases ( Figure 6C). According to the statistics, the lake volume does not have a significant correlation with the dam volume (Shen et al., 2020). The lake volume is determined by the dam height and valley geometry. Compared with the lake volume, the catchment area is usually more accessible and often used as the influencing factor in assessment indexes of landslide dam stability Stefanelli et al., 2016).
A stable landslide dam in New Zealand would require a 10-fold dam volume for a given lake volume (Korup, 2004). Nevertheless, some landslide dams among our inventoried landslide dams remain stable even though the lake volume is larger than the dam volume. For example, the Val Pola landslide dam in the Adda River is still stable due to a large dam width of 2.5 km (Crosta et al., 2004).

Inflow rate
There is a negative correlation between landslide dam stability and the inflow rate as shown in Figure 6D; 42% of landslide dams cannot remain stable when the inflow rate is smaller than 5 m 3 /s. Due to overtopping, landslide dams do not usually remain stable when the annual average inflow rate is larger than 10 m 3 /s. Some ancient landslide dams could exist for centuries with a small inflow rate from weak precipitation, and they are currently used for recreation and tourism, such as Valasht lake (Ehteshami-Moinabadi and Nasiri, 2019) and Sarez lake (Schuster and Alford, 2004).
The inflow rate of a landslide dam caused by rainfall and snowmelt is normally higher than the annual average river discharge, leading to a lower landslide dam stability. For example, the inflow rate at the blockage of the Hsiaolin landslide was 2,974 m 3 /s : significantly larger than the annual average inflow rate (30 m 3 /s). However, the actual inflow rate during the period of landslide dam formation is not usually recorded in historical cases, and thus, the average flow rate in the studied river is chosen to take its place (Shen et al., 2020).

Multifactor Analyses of Landslide Dam Stability
Based on the single-factor analyses discussed above, the following three factors are closely related to landslide dam stability: (1) the backwater lake; (2) the dam geometry; and (3) the dam material. Currently, the geomorphic approach is widely used to correlate the characteristic parameters of the dam and river with the stability of a landslide dam (Costa and Schuster, 1988;Casagli and Ermini, 1999;Ermini and Casagli, 2003;Korup, 2004;Dong et al., 2009). The parameters used in the indexes are all interlinked, and their relative significances are compared to fully evaluate landslide dam stability (Fan et al., 2020). Canuti et al. (1998) proposed the blockage index (BI): where V d is the dam volume (m 3 ) and A b is the upstream basin area at the point of blockage (km 2 ). BI reflects the contributions of the landslide volume and the drainage basin area to the stability of a landslide dam (Swanson et al., 1986).
Considering the dam height in the equation, Ermini and Casagli (2003) proposed a different formulation, the Dimensionless BI (DBI): where H d is the dam height (m). According to their study, the dam height is an important parameter to assess landslide dam stability against both overtopping and seepage failures. It influences the downstream slope in an overtopping failure and the hydraulic gradient in a piping failure . The DBI was obtained by analyzing 84 landslide dam cases. The two limits identifying the different domains were DBI = 2.75 for the stability domain and DBI = 3.08 for the instability domain . Chen et al. (2017) applied the DBI to assess the Attabad landslide dam stability: its DBI was 4.62-4.85 which confirmed the instability of the dam. Based on the dam height, Korup (2004) proposed three dimensionless indexes: the Backstow index I s , Basin index I a , and Relief index I r .
where V l is the lake volume. The values of I s < −3 and I s > 0 delimited the stable and unstable landslide dam domains, respectively. Similarly, I a = log H 2 d /A b . A landslide dam would be stable with I a > 3.
The relief index I r is expressed as follows: where H b is the relief upstream of the point of blockage. I r = −1 is used to distinguish stable from unstable landslide dams. Considering dam destabilization by the river using a simplified expression of stream power (A b S) per unit channel length (where S is the local longitudinal slope of the channel bed), Stefanelli et al. (2016) presented the hydromorphological dam stability index (HDSI): where V L is the landslide volume. Kumar et al. (2019) forecasted the damming process of the Urni landslide and employed the HDSI and DBI indexes to evaluate the stability of the potential landslide dam. The above geomorphometric parameters do not consider the effect of the dam material on the stability. Based on the DBI, the landslide dam stability is further discerned by use of the median particle size d 50 . The DBI was calculated by examining 229 landslide dam cases as shown in Figure 7A. A landslide dam was unstable with DBI > 3.6, whereas the data above this envelope remained inconclusive. Fifty landslide dams (Supplementary Appendix 1) with their grain compositions were extracted from these 229 cases. Three different domains of existence could be recognized (Figure 7B), as follows: Unstable domain, log (d 50 ) < 1.0. Below this threshold value, the dam materials mainly consist of debris grains smaller than 10 mm and have a marked tendency to erode. Therefore, the landslide dams are unstable. Uncertain domain, 1.0 < log (d 50 ) < 2.1. In this domain, the grain size distribution of the dam material is uncertain. It may be a gap-graded mixture or may be dominated by the median particle size, which makes it difficult to determine the landslide dam stability. Stable domain, log (d 50 ) > 2.1. Above the threshold value, the dam materials mainly consist of gravel and cobble with a high anti-erosion ability.
The database of landslide dam cases is used to test various evaluation indexes as shown in Table 2. It is difficult to discern the stability of a landslide dam using BI and I s . I a has an advantage to predict the instability of landslide dams. On the whole, DBI can be used to predict the stability and instability of landslide dams with an accuracy of approximately 60%. The prediction of landslide dam instability can be further improved by considering the dam materials.

STABILITY ANALYSES OF SPECIAL FACTORS ON LANDSLIDE DAMS
The stability of a landslide dam is determined not only by geomorphometric parameters, but also by special factors including aftershocks, surges, and engineering mitigation measures. These are discussed individually below.

Stability Analysis of Aftershocks on Landslide Dams
Landslide dams triggered by earthquakes are always subjected to many aftershocks. After the Wenchuan earthquake, more than 300 aftershocks with magnitudes higher than Ms 4.0 occurred after the formation of these landslide dams (Shen et al., 2013), which might significantly influence the soil properties and structures of these dams.
The effects of seismic action on embankment dams have been widely studied (Chen and Harichandran, 2001;Proulx et al., 2001;Wu, 2001;Swaisgood, 2003;Calayir and Karaton, 2005;Arabshahi and Lotfi, 2008;Chen et al., 2008;Bilici et al., 2009;Sevim et al., 2011;Zou et al., 2013). The stress concentration, soil liquefaction, ground deformation, and crushing damage in the zone of the slabs are the common failure patterns for embankment dams. These failure patterns are rarely observed for landslide dams due to the different dam materials and morphological characteristics.
A series of large-scale shaking table tests were conducted to investigate the dynamic behavior of landslide dams under various aftershocks (Shi et al., 2015b,c). Videogrammetry was employed to measure the dynamic deformation of the model dams. The seismic actions gave rise to settlements and horizontal deformations in the landslide dams due to loosened dam material. The settlements increased from the dam bottom to the dam crest. The displacements close to the slope surface of the dam were larger than those in the interior. Higher peak ground acceleration led to larger dam deformations due to higher seismic forces. Aftershocks may not directly lead to dam failure, but the cracks and settlements caused by the aftershocks may accelerate dam failure, accompanied by overtopping as shown in Figure 8.

Stability Analysis of Surges on Landslide Dams
The rising water level in backwater lakes may induce a series of landslides or avalanches (Xu et al., 2017). In addition, debris flows possibly migrate out of a gully and slide into a lake area (Hu et al., 2009). When a landslide or debris flow quickly rushes BI obtained from 255 landslide dams presented in the paper is smaller than 1.0 and it lies outside of the criteria proposed by Casagli and Ermini (1999). into a reservoir, it may result in huge surge waves (Hager et al., 2004;Koo and Kim, 2008). For example, a landslide of 300 million m 3 rapidly slid into the Vaiont Reservoir in Italy in 1963, and an enormous surge with a wave height as high as 300 m was generated (Tang and Lee, 1992). A glacial avalanche rushed into a moraine lake with a volume of 6.5 million m 3 in Nastetuku River, Canada, in 1983, and the moraine dam was completely breached in less than 5 h under the impact of large surge waves (Risley et al., 2006). The surges can significantly erode a landslide dam and reduce its width and height, causing more rapid overtopping failures than under normal conditions. The surge scale increases with an increase of the contact area between a landslide and the water surface and landslide height, directly reducing the stability of a landslide dam. Nevertheless, the effect of the distance from the entry point of a landslide to the dam site on the surge wave is limited (Wiegel, 1970;Xu et al., 2015). In addition, dynamic water pressure caused by the surge wave is applied on the upstream slope of a landslide dam (Chen H.Y. et al., 2015). The maximum pressure load increases with increasing sliding distance and significantly reduces the landslide dam stability by diminishing the effective stress.
As shown in Figure 9, the failure process of a landslide dam caused by a surge could be generalized as a three-stage process . In Stage I, a scour base plane is formed on the upstream slope of a landslide dam. Then, the scour surface intersects with the dam crest and gradually moves downstream under the action of waves. In Stage II, erosion occurs on both the upstream slope by surge waves and the downstream slope by overflow. In Stage III, an inclined erosion slope appears and much more intense erosion occurs during the overtopping process. Then, the erosion gradually slows down, leaving a residual dam behind. The wave height and water level play different roles during the process of dam failure . The wave height determines the erosion boundary above the water level, while the water level determines the erosion position on the upstream slope. The landslide dam stability is determined by the difference ( H) between the effective water level h w (the sum of water level h l and wave height h h ) and the effective dam height h d (the dam height after reduction due to local erosion and collapse). When H < 0, a stable eroded upstream surface eventually forms; when H > 0, the dam is overtopped and fails under the action of the subsequent wave loadings (Figure 9).

Stability Analysis of Engineering Mitigation Measures on Landslide Dams
If a landslide dam has a high tendency to fail, engineering mitigation measures might be employed to stabilize the blockage or accelerate its failure to reduce the potential risk . Yang et al. (2010) and Cui et al. (2012) introduced the experiences of handling some landslide dams induced by the Wenchuan earthquake. Schuster and Evans (2011) summarized the mitigation measures employed for 20 historical cases. Sattar and Konagai (2012) presented hazard-mitigation strategies for several large landslide dams. These measures could be categorized into short-term and long-term measures .
Short-term measures consist of diversion, drainage, and control of the erosion rate during overtopping failure. These measures are a temporary approach to control a landslide dam's stability. Diversion or drainage could be employed to stabilize the dam by controlling the rising reservoir level when surrounding hydraulic facilities, such as a reservoir, siphon, and pump, are available. This allows time to build the long-term measures. For example, diversion was applied to the Randa landslide dam formed in 1991 in the Vispa River, Switzerland, and a spillway across the dam was then constructed (Bonnard, 2004). Pumps were installed at the Higashi Takezawa landslide dam, and a 280-m-long concrete-faced spillway was then constructed (Sattar and Konagai, 2012). However, some measures may be used to accelerate a dam's instability to prevent continuous rising of the reservoir level. For example, excavation and blasting were employed for the Yanziyan landslide dam to eliminate large rocks to hinder further overtopping failure .
Long-term measures consist of drainage tunnels, spillways, and drainage conduits. These long-term measures are conducive to controlling the reservoir water level and are fundamental approaches to maintain the stability of the landslide dam. The Hongshiyan landslide dam triggered by the Ludian earthquake is a typical example (Shi et al., 2017). An existing drainage tunnel connecting the hydropower plant and reservoir became a diversion channel for the landslide dam. Spillways and drainage conduits are easier and cheaper to construct than drainage tunnels. Limited by available time and transportation facilities, constructed spillways are commonly designed to reduce the water level of the backwater lake rather than to serve as permanent structures . However, not all spillways are successful in preventing rapid overtopping failures. For example, the peak flow rate of the Yigong dam after constructing a spillway 1,000 m long, 24 m deep, and 150 m wide on the top and 20 m wide on the bottom was still as high as 124,000 m 3 /s (Shang et al., 2003).

Experimental Analyses of Overtopping Failures of Landslide Dams
Model experiments on the overtopping failures of landslide dams are valuable to gain more insight into the development process of dam failures and to calibrate and validate of the corresponding numerical models described below (Zhu et al., 2004).
Homogenous sand or gravel are currently considered as the model dam materials (Javadi and Mahdi, 2014;Liu et al., 2019). The silt and sand components below 1 mm are always excluded due to their cohesive properties, resulting in a complex overtopping failure mechanism (Schmocker and Hager, 2009). Cao et al. (2011b) observed that cohesive clays may act to mitigate seepage through the dam, and the gravel in the dam could appreciably depress the rate of the dam failure process. This is consistent with the observation of largescale model experiments by Morris et al. (2007) and Zhong et al. (2019). In contrast, coarser grains lead to a faster breach process compared with finer materials as indicated by Pickert et al. (2011). The reason proposed for the difference is that negative pore-water pressure is generated in the finer material with the apparent cohesion (Pickert et al., 2011). The grain size distributions (grain diameter greater than 0.75 mm) have a small effect on the overall overtopping failure process and no increased erosion resistance is observed for the soil mixtures (Schmocker et al., 2014). The results show that the mean grain diameter adequately describes the non-cohesive material characteristics and general overtopping failure features at laboratory scale can be investigated using uniform material.
According to the specifics of the breach process, overtopping failure can be divided into various stages (Chen R. et al., 2015). The overtopping failure process could be described as initial downstream surface erosion progressing to stair-stepped multiple overfalls, ultimately merging into a single upstream-migrating headcut (Hanson et al., 2002). Kakinuma and Shimizu (2014) categorized the overtopping failure process into the following four stages: downstream slope erosion, breach widens gradually, overflow rate increases significantly, and breach rate decreases. By comparing the field data with the overtopping erosion process, the final phase of overtopping erosion is overestimated. An armor layer appears on the downstream slope during the final stage of breach development in centrifugal model tests (Zhao et al., 2018(Zhao et al., , 2019. Commonly, the cross section of the breach during the overtopping process is considered as triangular, rectangular, trapezoidal, or parabolic (Fread, 1988;Pickert et al., 2011;Peng et al., 2014). The predominant presumption of a trapezoidal breach cross-section is based on observations of final breach shapes, like Tangjiashan and Xiaogangjian landslide dams (Chang and Zhang, 2010). However, the final breach shape is influenced by falling reservoir water levels and therefore does not represent the breach shape during the breach event (Coleman et al., 2002). Breach side walls are typically vertical and even overhanging during breach development (Morris et al., 2007;Pickert et al., 2011).
The overtopping failure patterns can be summarized as four types through the existing model experiments as shown in Figure 10. For Case (a), the downstream slope angle rapidly increases until a constant critical soil friction angle is achieved; thereafter, the angle is maintained to breach end (Guan, 2018). This model is consistent with the theoretical model reported by Powledge et al. (1989). For Case (b), the breach channel initially erodes the downstream face of the dam with an inverted slope parallel to the face; then, the breach inverted slope progressively flattens to a terminal value by rotating about a fixed pivot point along the dam base (Coleman et al., 2002;Schmocker and Hager, 2009). For Case (c), the channel first develops a stepped profile and the upstream migration of the steps coalesces into a headcut. Thereafter, the retreating headcut maintains a slope near the internal friction angle of the dam material (Walder et al., 2015). For Case (d), the erosion point moves from the downstream dam crest toward the upstream dam crest and the dam toe. Then, a spindle-like failure process is observed along the flow direction (Zhou G.G.D. et al., 2019;. The overtopping failure patterns are affected by various factors including the lake volume, downstream slope, and dam material. Compared with Cases (a) and (b), the lake volume is relatively small in Case (b), and the entrainment capacity is drastically reduced with the release of the lake volume during the failure process, resulting in a progressive flattening of the slope. The landslide dam is located in the flume with a large longitudinal slope in Case (d) (Zhou G.G.D. et al., 2019), whereas the flumes are horizontal in Cases (a), (b), and (c). As a result of a larger bed slope, the erosion stress and entrainment power of a beaching flood are larger in Case (d), leading to rapid erosion close to the crest. The dam materials in Cases (a) and (c) are fine-grained soils and the lake volumes are comparatively large; thus, a slope near the internal friction angle of the dam material could be maintained. The stepped profile in Case (d) is caused by the hydraulic jump when a breaching flood flows through the dam crest and downstream slope of the dam (Walder et al., 2015).
Multiple landslide dams triggered by strong shock or rainstorm are usually closely distributed along river reaches or gullies (Shi et al., 2015a). The outburst flood released by a landslide dam upstream can induce failures of downstream landslide dams one after another. Compared with a single landslide dam, the impact of cascading landslide dam failures is more complicated . For instance, Tangjiashan landslide dam and two other dams downstream (Shi et al., 2015a) as well as Xiaogangjian and Yibadao landslide dams (Niu et al., 2012) breached in succession. At least 19 landslide dams in the Sanyanyu gully were destroyed by upland flash floods and developed into a catastrophic debris flow in Zhouqu . The breach of cascading landslide dams may result in a sharp increase in the peak outflow rate and a more rapid breach process of any downstream dam (Cao et al., 2011a). In addition, multi-peak floods are very likely to develop at a downstream landslide dam, due to the overlapping effect of breach discharges upstream and downstream (Shi et al., 2015a). However, similarities and differences in the breach processes and failure patterns of single and cascading landslide dams are still unclear.

Numerical Analyses of the Overtopping Failure of Landslide Dams
To enable early warning and reduce risks and losses, it is absolutely necessary to develop a numerical model to predict the overtopping failure process of landslide dams (Peng and Zhang, 2013;Chen et al., 2017). Generally, the numerical models of overtopping failures of landslide dams can be classified as physically based models or coupled shallow water hydrodynamic models.

Physically Based Models
Physically based models adopt the principles of hydraulics, soil mechanics, and sediment erosion to calculate time-stepping solutions of the growth process of dam failure and the outflow hydrograph. Typical models of this type include BREACH (Fread, 1988;Fujisawa et al., 2009) and its modified models (Mohamed et al., 2002;Morris et al., 2009), BEED (Singh and Scarlatos, 1988), and BRES (Loukola and Huokuna, 1998). These models are widely employed to investigate overtopping failures of embankment dams. The geometrical parameters of the dam and initial breach as well as the material strength such as cohesion and internal friction angle are required before iterative calculation.
Compared with the abovementioned models, some modified models, such as DABA (Chang and Zhang, 2010), DB-IWHR (Chen Z.Y. et al., 2015;Zhong et al., 2018), and FIGURE 10 | Overtopping failure patterns of landslide dams: (A) downstream slope angle rapidly increases until a constant critical soil friction angle is achieved and then the angle is maintained to breach end; (B) downstream dam face is initially eroded with an slope parallel to the face, then the breach inverts slope, progressively flattening to a terminal value; (C) a stepped profile is first developed and then the retreating headcut maintains a slope near the internal friction angle of the dam material; and (D) a spindle-like failure process is observed along the flow direction.
DLBreach (Wu and Li, 2017), have been developed to analyze the overtopping failure process of landslide dams. The dam material is commonly assumed to be uniformly distributed along the depth for homogeneous embankment dams (Figures 11A,B). However, for a landslide dam, the erodibility changes significantly from the dam crest to the native foundation due to changes in the soil type, density, and grain composition (Chang et al., 2011). Variations in the soil erodibility along the depth are thus considered in the DABA model. The erosion processes of Tangjiashan and Xiaogangjian landslide dams have been successfully simulated by using the DABA model. Commonly, the lateral banks are assumed to uniformly widen (Figures 11C,D). If the slope stability of the lateral banks is considered, bank erosion occurs only within the water area at the initial stage. The lateral banks undergo undercutting by the breaching flood and the debris soil on the lateral banks above the water surface slides down to the breach bottom at a critical slope (Chen Z.Y. et al., 2015;Zhong et al., 2018). For the cascading breach of landslide dams or surge wave in the backwater lake, the sudden increase in the inflow rate would cause erosion both within and outside of the breach (Shi et al., 2015a;Wu and Li, 2017). The erosion on the dam crest and breach should be considered, and the overflow discharge outside and within the breach should be combined as the total outflow discharge (Figures 11E,F).

Coupled Shallow Water Hydrodynamic Models
Coupled shallow water hydrodynamic models are based on the mass conservation equations for the sediment and dam material, and the mass and momentum conservation equations for the water-sediment mixture flow (Cao et al., 2011c). This method has been widely employed to investigate the overtopping failure process (Huang W. et al., 2014;Chen S.C. et al., 2015;Do et al., 2016;Zheng et al., 2016;Liu and He, 2018;Wang, 2018;Wu and Lin, 2019). The calculation cost of coupled shallow water hydrodynamic models is significantly larger than that of physically based models due to the solution of conservation equations. Nevertheless, the former class of models offers a higher spatial resolution. To consider sediment transport and bed morphological evolution, a double layer-averaged model has been proposed by Li et al. (2013) and Ouyang et al. (2014). The influences of seepage flow through the dam material on the apparent cohesion and overtopping failure initiation have been investigated (Volz et al., 2017). In recent years, some novel solution algorithms and graphics processing units (GPUs) have been employed to achieve speedups of up to two orders of magnitude (Dazzi et al., 2019). This is very beneficial for predicting the overtopping failure processes of landslide dams with large volumes.
In addition, other computational fluid dynamics methods are also used to investigate overtopping failure of landslide dams like the smoothed particle hydrodynamics (SPH) approach (Memarzadeh et al., 2018) or discontinuous deformation analysis (DDA) coupled with SPH method . Two-phase formulations are used in the SPH numerical algorithms to examine the free surface and bed evolution profiles, in which the entrained sediments are treated as a different fluid component (Ran et al., 2015). Compared with these numerical methods, coupled shallow water hydrodynamic models offer the advantage of high computational efficiency. This  (Chang and Zhang, 2010); (C,D) uniformly and non-uniformly widening the breach, respectively ; (E) outflow discharge within the breach; and (F) outside and within the breach (Wu and Li, 2017). review will not introduce other computational fluid dynamics methods in detail.
Unfortunately, the overtopping failure predictions of landslide dams contain some uncertainty (Coleman et al., 2002;Zhu et al., 2004). The overtopping failure process of a landslide dam is believed to depend on the various factors introduced in Section "Influences of Characteristic Parameters on Landslide Dam Stability, " including the landslide triggers, geomorphic characteristics, and dam materials. These factors vary greatly from case to case. At present, the numerical models of overtopping failures in landslide dams have some disadvantages. After a landslide dam has formed, it is difficult to rapidly obtain the soil parameters of the dam. Almost all of the simulations available are limited to overtopping failures in homogeneous dams, which is not consistent with the heterogeneity of landslide dams. In addition, the coupled erosion rate equations in the shallow water hydrodynamic models are empirical, and the entrainment parameters need to be calibrated for every landslide dam case. Furthermore, the capacity of the breaching flood to transport sediment has not been included in most modeling.

SEEPAGE FAILURE ANALYSES OF LANDSLIDE DAMS
Piping and flowing soil are typical phenomena of seepage failure in landslide dams. Laboratory tests and numerical simulations are the main research approaches to investigate the seepage failure of landslide dams.

Experimental Analyses of Seepage Failure
The experimental research investigating the seepage failure of landslide dams could be categorized into element and model tests. The element tests of seepage failure are also known as permeability tests, and the specimen sizes used in these tests are small. Model tests carried out with appropriate proportions and similar materials resemble prototype testing.

Element Tests of Seepage Failure
The failure process of piping and flowing soil could be observed in various experiments. Based on the pore pressures in sandy soils, four stages of piping development were identified: initial movement, progressive heave, boil formation, and total heave Rice, 2013, 2014). For widely graded material, the finer fraction in the initial stage moves as uniform loss across the entire base and a distinctly preferential path (piping) develops in the last stage (Moffat et al., 2011;Huang et al., 2015). Influenced by the amount and types of fine particles present in the noncohesive and cohesive soils, three modes of piping behaviors were recognized: concentrated leak erosion, backward erosion, and suffusion (Richards and Reddy, 2012). The flowing soil failure of uniform sand is associated with an effective stress equal to zero, whereas a piping failure is caused by the internal erosion of fine particles in gap-graded sand (Ke and Takahashi, 2014;Van Beek et al., 2014;Yang and Wang, 2017).
The permeability coefficient and critical hydraulic gradient could be obtained from an element test to evaluate seepage stability (Richards and Reddy, 2007;Fleshman and Rice, 2013;Zhou et al., 2016;Zhu et al., 2019). The critical hydraulic gradient for seepage failure caused by flowing soil may be consistent with Terzaghi's theoretical value (Fontana, 2008;Huang et al., 2015;Yang and Wang, 2017). However, the critical hydraulic gradient for seepage failure caused by piping erosion is smaller than Terzaghi's theoretical value.
Based on grain size measures such as d 15 , d 85 , and grain sizes ranging from d to 4d, the Kezdi (1979), Sherard (1979), and Kenney and Lau (1986) criteria have been proposed to ascertain the seepage stability. These geometric criteria were verified and modified by subsequent researchers (Burenkova, 1993;Skempton and Brogan, 1994;Richards and Reddy, 2007;Li and Fannin, 2008;Wan and Fell, 2008;Shire et al., 2014;Zhou et al., 2016). Considering the high coarse fraction in the landslide dam materials, Chang and Zhang (2013) extended the internal stability criteria for gap-graded and well-graded soils based on a physical understanding of the microstructures of the soils.
Grain compositions of landslide dams differ significantly between various landslide dams or even for a specific dam, resulting in different seepage failures. Shi et al. (2018) conducted seepage tests on four typical dam materials: fine-grained, coarsegrained, well-graded, and gap-graded materials at various dry densities, based on Tangjiashan landslide dam (Chang et al., 2011;Zhao et al., 2013). The seepage failures of the fine-grained and gap-graded soils were flowing soil . By contrast, they were piping for the coarse-grained and well-graded soils.

Model Tests of Seepage Failure
The prediction index of seepage failure of landslide dams is discerned in various experiments. The critical hydraulic gradients corresponding to the onset of seepage erosion and collapse of the dam crest were found to increase with an increase in the uniformity coefficient of dam material (Okeke and Wang, 2016a). The potential to form a piping path through the dams is reduced with an increase in the soil density and the homogeneity of the dam materials (Okeke and Wang, 2016b). The premonitory factors during the seepage failure were identified by Wang et al. (2018) as self-potential change, pore-pressure change, and seepage-water turbidity.
The seepage failure of landslide dams can be classified into four patterns based on the dam materials as shown in Figure 12. For a fine-grained landslide dam, the crest width is shortened by the slide of the downstream slope, and then the sliding area on the downstream side increasingly expands to the entire downstream slope: a typical sliding failure Zhu et al., 2019). This failure pattern is a combination of a flowing soil failure and an overtopping failure. For a coarsegrained landslide dam, the inflow rate and seepage flow rate basically maintain a balance and dam failure is not observed (Guan, 2018). For a well-graded landslide dam, the flow erodes the crest and downstream surface after the water level rises above the crest: a typical erosion failure (Xiong et al., 2018). This failure pattern is a combination of internal erosion and overtopping failures.
Step-pool systems develop during the failure process, which is consistent with observations of the Tiger-leaping Gorge and Yujunmen landslide dams (Wang et al., , 2012. For a gap-graded landslide dam, some fine grains are removed from the cracks driven by hydraulic gradient and these cracks easily broaden and interconnect, which eventually contributes to the formation of the piping channel: a typical piping failure (Xiong et al., 2018). This failure pattern was applicable to the landslide dam in the La Paz river catchment (Quenta et al., 2007).

Numerical Analyses of Seepage Failure
Compared with element and model tests, numerical analyses of seepage failure can provide a great deal of insight into the seepage process and hence the failure mechanisms of landslide dams. Currently, continuous medium methods and discontinuous medium methods are employed to analyze the seepage stability of landslide dams.

Continuous Medium Method
The finite-element method (FEM) and finite-difference method (FDM) are typical continuous medium methods. The distributions of the seepage field, stress, and displacement can be obtained using these methods, and the tendency of seepage failure is then determined as shown in Figure 13.
Considering the detachment of the soil grains from the soil fabric, backward erosion and the development of the piping path are presented as they vary with the grain compositions, porosity, and pore pressure (Fujisawa et al., 2010). The critical hydraulic gradient and the progression of internal erosion could be obtained by capturing the main hydraulic characteristics of the turbulent flow that occur in an erodible pipe and the seepage flow in the remaining area of a dam foundation (Wang et al., 2014;Rotunno et al., 2019). The heterogeneity could accelerate the development of preferential flow paths and increase the likelihood of seepage failures by verifying the stochastic parameters of the dam material, such as the hydraulic conductivity, void ratio, and grain contents (Liang et al., 2017). A local critical gradient rather than an average critical gradient FIGURE 12 | Failure patterns of a landslide dam for different materials: (A) sliding failure for a fine-grained landslide dam Zhu et al., 2019); (B) no failure for a coarse-grained landslide dam; (C) erosion failure for a well-graded landslide dam (Xiong et al., 2018); (D) piping failure for a gap-graded landslide dam (Quenta et al., 2007;Xiong et al., 2018).
FIGURE 13 | Evolution of (A) saturation degree, (B) displacement field, and (C) shear strain during a seepage failure process (Xiong et al., 2018). should be used as the criterion for pipe progression to account for the scale effects observed between large-and small-scale experiments (Robbins, 2016). Based on a rational constitutive model for a saturated/unsaturated soil (Zhang and Ikariya, 2011), Xiong et al. (2018) correlated fine-grained, well-graded, and gapgraded dam materials with the sliding failure, erosion failure, and piping failure by analyzing the development of the seepage lines, displacement field, saturation degree, and shear strain.

Discontinuous Medium Method
Compared with the continuous medium method, the discontinuous medium method has the outstanding advantage of simulating the interactions among particles. The discreteelement method (DEM) is a typical method based on the discontinuous media theory.
The process of internal erosion in landslide dams could be simulated by progressively removing the finer particles, and the mechanical consequences of erosion can thus be analyzed. The removal of particles produces an increase in the specific volume because of the solid volume decrease and void volume increase. The removal of fine grains leads to a decrease of the sliding resistance of each interparticle contact and the occurrence of local sliding, which contributes to the instability of the granular assemblies at a shear stress level much lower than the critical state failure line (Muir et al., 2010;Hicher, 2013). The number of contacts between the fine particles is significantly reduced by the fine particle loss, while the contact forces gradually transfer to the coarse soil particles .
Alternatively, the migration process of the soil particles moved by the pore fluid could be directly simulated by coupling DEM with computational fluid dynamics Zhang et al., 2020). The representative size ratio of the soil skeleton has a great influence on the effectiveness of preventing the dam soil from being eroded (Huang Q.F. et al., 2014). The eroded percentage of soil particles gradually grows with the increase of the representative size ratio of the soil skeleton, and the effective vertical stresses reach zero when the hydraulic gradient reaches the estimated critical hydraulic gradient (Abdelhamid and El Shamy, 2015;. The internal erosion rate is proportional to the flow velocity for both spherical and non-spherical particles, and a critical velocity exists for angular particles owing to grain interlocking which is not observed for the spherical particles (Guo et al., 2018). As shown in Figure 14, the seepage field is distributed uniformly in the fine-grained dam material which results in a global seepage failure (flowing soil), whereas it is gradually deflected with the increase in the hydraulic gradient and local seepage failure (piping) occurs in the coarse-grained dam material (Cheng et al., 2018;Shi et al., 2018).

FUTURE RESEARCH DIRECTIONS
The research on landslide dam stability is of enormous significance for early warnings and emergency evacuations and, consequently, also for disaster mitigation after dam failure. Considerable progress has been made over more than 30 years in characteristic analysis and developing some understanding of the failure mechanisms. Nevertheless, a complete understanding of the stability evaluation and failure mechanisms of a landslide dam remain elusive, and the state-of-the-art in failure modeling is far from advanced. The suggestions are listed as follows: (1) The erosion rate of dam materials should be taken into account in the multifactor evaluation of the stability of a landslide dam. The erosion rate of dam materials can be obtained by in situ tests (i.e., Chang and Zhang, 2010;FIGURE 14 | Seepage failure modes in (A) fine-grained and (B) coarse-grained soils . The red dashed lines denote the seepage failure tendency. Flowing soil and piping failures occur for the fine-grained and coarse-grained soils, respectively. Chang et al., 2011) and laboratory model tests. However, most landslide dams last a short time, mainly because of the high erosion rate (Chang et al., 2011). It is necessary to establish the functional relationship between the grain composition and erosion rate by laboratory model tests to achieve the goal of rapid assessment. A similar approach is that the internal erosion of dam materials caused by seepage can be evaluated from the grain composition, like the Sherard (1979) and Kenney and Lau (1986) criteria. Moreover, the erosion function of the dam materials could be embedded in physically-based models and coupled shallow water hydrodynamic models. The overtopping failure process of a landslide dam could then be rapidly predicted with high accuracy. (2) The morphological characteristics of a landslide dam can be rapidly obtained by remote sensing or unmanned aerial vehicle (UAV). However, the grain composition and material distribution are difficult to determine. Traditionally, the grain size distribution of a landslide dam is obtained by sieve analysis. However, sieve analyses are impractical when dealing with materials often ranging in size from blocks tens of cubic meters in size to microscopic particles. The sampling method of grid by number analysis has been proposed to obtain the coarser part of the debris material . Considering landslide dams located in remote mountainous areas, hyperspectral remote sensing technology is suggested to obtain the grain size on the surface of landslide dams. The hyperspectral system is now able to cover the wavelength region from 0.4 to 2.5 µm using hundreds of spectral channels . The grain size can be determined by the hyperspectral technology in tandem with machine learning methods. (3) Many research studies have been conducted for single landslide dams, investigating the breach discharge and failure pattern (Coleman et al., 2002;Schmocker and Hager, 2009;Walder et al., 2015;Jiang et al., 2016;Zhou G.G.D. et al., 2019). The outburst flood released by a landslide dam upstream can induce failures of downstream landslide dams one after another (Shi et al., 2015a). The breach of cascading landslide dams may result in a sharp increase in the peak outflow rate and thus a more serious disaster downstream. These research results of single landslide dams may not be suitable for cascading landslide dams due to the complex breach process and unsteady inflow rate (Shi et al., 2015a). In addition, the research attention paid to the cascading failure of landslide dams has been relatively limited. The effects of grain composition and geomorphic characteristics (dam geometry and initial water level) of landslide dams upstream and downstream on the cascading failure process and the amplification effect of breach discharge have not been explored. (4) Significant differences in grain compositions occur in various dam zones because of grain segregation in the landslide dam accumulation process (Zhou Y.Y. et al., 2019). The grain composition near the bank may be different from that away from the bank. It is also applicable to the grain compositions in the dam upstream and downstream. The failure process and pattern caused by overtopping and seepage could be affected by this heterogeneity. However, the model dam is normally prepared manually by uniformly mixing the debris material. A sectionalized dam model should be considered in the corresponding model experiments and numerical analyses. (5) The evolution of the flow properties during an overtopping failure should be considered in the corresponding model experiments and numerical analyses. The released flow changes from a pure water flow into a sediment flow or even a debris flow by substantial entrainment of dam materials Zhou G.G.D. et al., 2019). The flow density and viscosity increase during the breach process which affects the erosion potential of dam material. Despite the importance of this, the enlargement process of debris flows during an overtopping failure is relatively unexplored (Chen et al., 2014;Chen H.Y. et al., 2015;Jiang et al., 2016;Zheng et al., 2016). In addition, the capacity of the overflow to transport sediment should be considered in the modeling. (6) The effects of seepage on the overtopping failure of landslide dams should be investigated. At present, the overtopping and seepage failures are separately analyzed. Actually, seepage has a significant effect on the breach process and breach mode during overtopping. On the one hand, significant positive pore pressure occurs in the dam below seepage line and reduces the effective stress in the downstream slope. The pore pressures display an obvious difference among various dam materials and thus affect the overtopping failure . On the other hand, negative pore pressure is generated in the finer material with the apparent cohesion above the seepage line (Pickert et al., 2011). However, it is difficult to observe in the coarsegrained material. Coarser grains with a higher strength lead to a faster breach process compared with finer materials with a smaller strength (Morris et al., 2007;Zhong et al., 2019). This seemingly contradictory conclusion is closely related with pore pressure.

AUTHOR CONTRIBUTIONS
HZ and ZS finished the manuscript. DS and CM collected the data and landslide dam cases. MP, KH, and LZ modified the manuscript. All authors contributed to the article and approved the submitted version.

FUNDING
The research reported in this manuscript was substantially supported by the Natural Science Foundation of China (Nos. 42007252 and 41731283).