4.1 Influence of the erodible layer on the sliding mass motion
From the velocity and displacement of the sliding masses, it is clear that the presence of the erodible layer on the horizontal basin had a significant restraining effect on the mobility of the sliding mass. Under the restraining effect, not only was the movement distance of the sliding mass reduced, but also were the length and width of the deposits. The main reason for this was that the density and depth of the erodible layer increased due to interactions with the sliding mass. It was also for this reason that the topography in the front edge of the deposits elevated and the shear strength of the erodible layer was fully exerted under the force of the sliding mass, which increased its resistance during motion. After the collision between the sliding mass and erodible layer, the kinetic energy of the sliding mass was consumed primarily by shearing and colliding among particles. On the other hand, the friction coefficient in the interior regions of the sand was 0.613, which is larger than the friction coefficient for the sand–plate interface (0.403). Therefore, the mobility of the laboratory landslides without an erodible layer was larger than that with an erodible layer (Crosta et al., 2017; Liu et al., 2020; Yuan et al., 2014; Zhou et al., 2019).
To reduce the influence of material differences on the experimental results, the materials of the sliding mass and erodible layer were both medium-fine sand, with the same physical and mechanical attributes. This type of sand possessed small rheological strength. In some model experiments using dry particles as the erodible layer, experimental results have shown that the erodible layer had an inhibitory effect on the mobility of the sliding mass (Crosta, 2013; Crosta et al., 2017; Dufresne, 2012). It should be noted that all of these studies set the erodible layer on the nearly horizontal basin. Nevertheless, when the erodible layer was set on a downward sliding path, such as an inclined plate, opposing results were found(Mangeney et al., 2010). Further, the mobility of the sliding mass would increase due to entrainment during collisions with the erodible layer, irrespective of whether the underlying layer was dry or wet, especially when the slope angle of the erodible layer was close to the internal friction angle of the material (Crosta et al., 2015a; Mangeney et al., 2010).
In this study, when the De was less than 20 mm, the larger the De was, the more significant the hindering effect was on the mobility of the sliding mass, especially in the planform of the deposits. However, when the De was 24 mm, the length and width of the deposits increased, from 287.59 mm and 630.00 mm to 291.34 mm and 650.60 mm, respectively. At the same time, the erosion volume decreased from 905.14×102 mm3 to 190.40×102 mm3 due to the decreasing kinetic energy of the sliding mass and its weakening erosion ability.
The position of the erodible layer may play an important role on the mobility of the sliding mass. Both the mobility and duration of the sliding mass will increase as the increasing thickness of the erodible layer when it is at an inclined plane (Mangeney et al. (2010), Fig. 11). However, the increasing mobility will be maximum at De = 4d to 24d (particle diameter) when the sliding mass is compacted (Farin et al. (2014), Fig. 10). With the propagating of the sliding mass, its velocity decreases, and consequently decreasing in erosion depth. During its slowly propagating, there is almost no erosion (Farin et al. (2014), Fig. 17). The larger the velocity of the sliding mass is, the stronger its erosion ability is, and the deeper the erosion depth is, which is also true for the present study. The difference of this study with those previous studies is that the mobility of the sliding mass is inhibited when the sliding mass is on the horizontal plane (Fig. 7). It can be seen that the position of the erodible layer on the landslides’ motion path has significant influences on its mobility and erosion ability.
4.2 Interaction form and structural characteristics
In physical model experiments, it is common to analyse the interaction between the sliding mass and erodible layer through the internal structural characteristics observed from a profile of the deposit (Deng et al., 2015; Paola et al., 2009). The structures of buckles, recumbent folds, thrust shear, horizontal shear, and convoluted bedding revealed in the profile of the deposits in this study indicated that there was an intense interaction between the sliding mass and erodible layer during sliding mass motion (Fig. 9, 10).
With increases in the De, the effect of the sliding mass on erodible layers changed from pushing to covering (Fig. 10). According to the profile features, the sliding mass in zone II mainly acted by pushing, scraping, and carrying the erodible layer to zone Ib. The structures formed in zone II are relatively simple, where some ridges formed on the deposit surface (Fig. 10e, h). However, the ridges disappeared at thicker erodible layers (Fig. 10k). The range of zone II along the y axis became shorter with increases in De. Zone Ib was mainly characterized by compression deformation, demonstrated by a relatively high compression of the strata and complex internal structures, such as buckles, recumbent folds, thrust shear, horizontal shear, and convoluted bedding. There were small bulges (Fig. 10j) and no fully covered area by the sliding mass (Fig. 10f) in zone Ib. These phenomena were due to that the resistance between the erodible layer and horizontal plane is greater than its internal one when the sliding mass interacted with it. The erodible layer in zone Ia is mainly subjected to thrusting. There are obvious folds formed at the front of the deposits (Fig. 10d). With continuous increases in De, the extent of compression of the erodible layer in zone Ia decreased, and the extension of the y axis gradually decreased until the area disappeared (Fig. 10c). When the erodible layer is thick, the pushing effect of landslide is weakened, and the front sliding mass will over the erodible layer without strong disturbance. The upheaval formed in the deposits separated zone II from zone Ib, which is usually caused by compression and accompanied by rapid decline in velocity of sliding mass. When the inner shear stress of the debris in motion is not sufficient to cause discrete shear and inverse thrust, longitudinal compression would take place, forming the transverse ridges in the debris mass (Strom, 2006). This phenomenon is also often found in the field, such as in scenarios where suddenly decreases in the topographic gradient; encounters with soft and/or deformable substrates; or slide into water (Clavero et al., 2004; Johnson, 1978; Shea et al., 2007; Siebert, 2002).
Folds are common geological phenomena, which form not only under tectonic stress (Deng et al., 2022), but also under sudden loading caused by turbulence or mountain floods (Rana et al., 2016). The friction resistance or shear stress applied to weak strata is the basis for the formation of folds (Owen, 1996). The folds in this study were distributed in the erodible layer close to the upheaval and front edge of the deposit (Fig. 10, Ib). Furthermore, a series of compressive deformation characteristics are also found in zone Ia (Fig. 10). Through model experiments with coloured sand, Dufresne et al. (2009)reproduced the folds sustaining its initial colour sequences and reported that there are bulldozing effects in the substrata when being covered rapidly by the propagating avalanche. A resistance was faced by the lower part of the sliding mass, while the upper part continued to move forward and into the core of the fold when the sliding mass scraped the erodible layer. Under the influence of inertial force, the sliding mass convoluted into the fold moved forward along with the folded erodible layer under the force of the trailing sliding mass, leading to the convolution phenomenon in zone Ib. The formation process of this phenomenon can explain the thrust nappe structure in loess landslides (Peng et al., 2018).
Thrust shear tends to occur in positions where the stress is most concentrated. These shear stresses may combine with the force of gravity to further deform folds in complex ways (Lowe, 1975). In this study, the shear position in the deposits mainly developed in the core and limbs of the recumbent folds.
The secondary shear surfaces in the deposits mentioned in Section 3.3 were caused by the stress concentration and velocity differences in the sliding mass. The secondary sliding surfaces formed by thrust shear primarily existed in the interior of the deformed erodible layer and the contact area between the erodible layer and sliding mass. At the thin erodible layer, there were few secondary sliding surfaces with a shallow burial depth, but a number of complex secondary sliding surfaces in the medium thickness erodible layers (Fig. 10a, b and c). At De = 24 mm, the increase in resistance led to an upward movement of the secondary sliding surfaces, which reduced their number. Further, thrust shear is common in landslides, and is observed in the Ancient Longwu Xishan Landslide No.2 (Tian et al., 2022)
Although many structural phenomena were reproduced in this study, there are still some others not observed, such as flash heating and consequently friction weakening (Habib, 1975; Kent, 1966; Shreve, 1966), and fluidization of substrata due to the impact of the sliding mass that results in rearranging of the soil skeletal particles (Sassa, 1998). In the study of (Zhou et al., 2016) to the entrainment of substrata containing water, the sliding mass entrains in the manner of impact and scour, and the water plays an important role during the process. However, our experiments is conducted without the engagement of water, which makes the changing pore water pressure and consequently fluidization of substrata impossible (Duan et al., 2021; Iverson et al., 2011). In addition, diapiric intrusions, which are considered due to compression (Wang et al., 2018) or friction shearing (Phillips et al., 2013) in the substrata, are also not observed. Therefore, further experiments need to be performed to reproduce these phenomena.
In fact, during the interaction between the sliding mass and erodible layer, the sliding mass impacted the erodible layer with a certain velocity. The erodible layer was shovelled, scraped, and pushed by the sliding mass. The erodible layer dissipated the kinetic energy of the sliding mass through the collision and friction among particles, thereby restraining the velocity and displacement of the sliding mass. The base of the moving mass accommodated transport by large amounts of simple shear (Roche et al., 2013). The erosion ability of the sliding mass with a certain volume is finite. When De was small, the interaction between the sliding mass and the erodible layer was gentle. Therefore, the kinetic energy of the sliding mass dissipated due to fewer collisions and less friction among particles. The sliding mass could push most of the erodible layer on its path a longer distance along the bare board. As a consequence, the position of the fold core was far from the slope break. Owing to the hindrance of the erodible layer on the sliding mass being small, the sliding mass could easily climb over the erodible layer and move horizontally on top of the erodible layer. Consequently, the shear surface was relatively gentle. In contrast, when De was large, the interaction between the sliding mass and the erodible layer was significant; therefore, there was a large fraction of kinetic energy loss in the sliding mass. The sliding mass had to move upward along the erodible layer under the thrust of the sliding mass and the resistance of the erodible layer. As a consequence, a steeper shear surface was formed in the deposit. The geological structures formed in this study are highly consistent with the structural evidence found in the field. Field evidences suggest a bulldozing action occurred at the front of the Nixu rock avalanche (Zeng et al., 2021). The highspeed avalanche debris in motion ploughed the erodible alluvial substrates along a basal decollement like a bulldozer. Complex structural phenomena such as entrainment and diapiric structure were formed in the deposits. This work reproduced inner structures seen in natural cases and reflected the kinematics in the final deposits.
4.3 Generalization of the sliding mass movement process
The motion process of the sliding mass was divided into three stages, according to differences in their velocity and motion characteristics.
In the first stage, the velocity of the sliding mass increased sharply to a peak velocity before colliding with the erodible layer (Fig. 11a).
In the second stage, there was a strong collision between the sliding mass and the erodible layer. The length of the sliding mass became significantly shorter, and the velocity of the sliding mass was significantly reduced. While squeezing, scraping, and pushing the erosion layer forward, the sliding mass was also subjected to increasing resistance due to the erodible layer (Fig. 11b). The disturbed erodible layer moved forward along with the sliding mass, which resulted in an increase in the erosion range and a further compressed erodible layer (Fig. 11c). The horizontal thrust of the sliding mass reached a maximum when all of the sliding mass slid on the horizontal basin. Therefore, secondary shear surfaces were formed at the position of the stress concentration (Fig. 11d). The strong interaction between sliding mass and erodible layer led to the sliding surfaces formed not only on the interface between the sliding mass and horizontal basin, but also at the contact surface between the sliding mass and erodible layers and the inner regions of the erodible layers in the second stage. Therefore, the movement mode of landslide was relatively complex and rich structural features were formed.
In the third stage, the interaction between the sliding mass and the erodible layer became gentle, and the disturbed erodible layer became dense under the pressure of the sliding mass. The erodible layer entrained by the sliding mass eventually stopped moving forward at a lower velocity. During this motion process, there was a small relative displacement between the front sliding mass and the erodible layer. Furthermore, a series of buckles formed at the front of the erodible layer under the squeezing of the sliding mass. In this stage, a secondary acceleration was observed when the De values were 4 mm, 12 mm, 16 mm, and 20 mm because the trailing sliding mass, under the action of an inertial force, would climb over and surpass the front sliding mass obstructed by the disturbed erodible layer.
4.4 Comparisons with previous studies
In sandbox experiments carried out by Dufresne (2012), there was an arc-shaped joint between the inclined and horizontal basins. The angles of the plate in the sliding mass zone and the accumulation zone were set 60° and 0°. Experimental results showed that the effects due to the sliding mass interacting with the erodible layer were mainly horizontal shearing. When the sliding mass pushed and collided with the erodible layer, the deep layer of the erodible layer was disturbed. However, the erodible layer was not completely eroded.
In a study by Crosta et al. (2017), the experimental device used was similar to that in this study. Structures such as convolutions, buckles, and thrust shearing were also found in the deposits (Fig. 12a-c), which further demonstrated that these structural phenomena tend to form under the condition that there was an erodible layer on the horizontal basin. However, with an erodible layer of 21mm it is hard to answer how the internal structures vary at different De. In addition, the influence of the De and the characteristics of the amount of eroded layer have not been quantified. Found in this study, With the increase of De, the action form of sliding mass on erodible layer changes gradually from pushing to covering. At the same time, there is a certain relation between internal structure and surface structure.
The material of the erodible layer used in this study was dry medium-fine quartz sand. This kind of sand is very loose under conditions without compaction; therefore, it can easily be disturbed. In actual landslides, although the erodible layers differ due to different geo-environments, they all have certain structural properties after long-term consolidation. It is difficult for a landslide to cause a large-scale or deep scraping effect when the landslide impacts the erodible layer. Therefore, recumbent folds are rarely observed in actual landslides. Conversely, the folds, thrusts, and pushover phenomenon on the erodible layer of landslides are more common and significant in actual landslides.
In addition, the volume of landslide also profoundly affects its erosional ability. For instance, the Yigong landslide primarily comprises weathered granite, and the erodible layer is sandy alluvial deposits (Wu et al., 2020). The landslide, with a volume of about 91×106 m3 along the gully moving downwards, squeezes the erodible layer, with a volume of about 24×106 m3, to form pushover structures. During this process, the erodible layer is scraped and entrained by the sliding mass, which makes the volume of the landslide continuously increase. The large impact force generated by the landslide takes away the erodible layers in the gully, which results in an outcrop of the bedrock in the eroded area (Delaney and Evans, 2015; Zhou et al., 2016). The underlying bedrock prevents deeper erosion and the flow obviously do not lose momentum in these sections (Dietrich and Krautblatter, 2019). This field observation is reflected in the experiment conducted herein, when the De is 4 mm. The impact force of the sliding mass makes the erodible layer at the slope break completely erode, resulting in direct contact between the sliding mass and the plexiglass plate. This experiment reproduces the geological phenomena observed in the field.
Another case, the Nixu landslide is also mainly composed of granite, and its erodible layer is mainly horizontal alluvial deposits, composed of coarse sand and gravel (Zeng et al., 2020). The landslide exhibits an intense erosion ability when moving downwards with a volume of about 47×106 m3, and the erodible layer undergoes strong deformation and distortion. Structures such as crushed substrate clasts, convoluted lamination, diapiric intrusion, decollement, and sand boils are manifestations of deformation and distortion, as shown in Fig. 12d–h (Wang et al., 2018; Zeng et al., 2021). Under this condition, the interaction between the landslide and the erodible layer, as well as the structures in the deposit, are similar to those at De is 24 mm in this study. These structures reflect the large resistance of the erodible layer to the sliding mass. However, due to the differences in the material properties of the landslide and the erodible layer in the field, the structural phenomenon is more widespread than that in laboratory experiments.
Finally, the Xingyuan landslide consists of silty clay, and the erodible layer is alluvial deposits comprising silty clay, sandy silt, and gravel (Duan et al., 2021; Peng et al., 2017). The landslide, with a volume of 0.17×106 m3, impacts the erodible layer consisting of silty clay, resulting its liquefaction. The position of the sliding surface in the landslide moves down 5–6 m of the ground surface when the landslide falls onto the erodible layer, then the landslide covers the undisturbed erodible layer (Fig. 12i and j). This is similar to the experimental conditions of De is 24mm in this study.