Next Article in Journal
Application of Different Weighting Schemes and Stochastic Simulations to Parameterization Processes Considering Observation Error: Implications for Climate Change Impact Analysis of Integrated Watershed Models
Previous Article in Journal
Natural Swimming Ponds as an Application of Treatment Wetlands—A Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

A Review of Research Methods and Evolution Mechanisms of Landslide-Induced Tsunamis

1
State Key Laboratory of Geohazard Prevention and Geoenvironment Protection, Chengdu University of Technology, Chengdu 610059, China
2
College of Environment and Civil Engineering, Chengdu University of Technology, Chengdu 610059, China
*
Author to whom correspondence should be addressed.
Water 2023, 15(10), 1879; https://doi.org/10.3390/w15101879
Submission received: 22 March 2023 / Revised: 28 April 2023 / Accepted: 11 May 2023 / Published: 16 May 2023
(This article belongs to the Section Water-Energy Nexus)

Abstract

:
The research on landslide surge has always been a hot issue in the research of engineering geology, geological disasters and hydraulics in both domestic and foreign circumstances. Based on the research status and achievements of landslide surge, this paper mainly reviews the research methods of landslide surge including numerical simulation methods, theoretical analysis methods and physical methods. The first part summarized the examples of geological disaster surge and the classification, characteristics and influencing factors of geological disaster surge. In the second part, the research methods of landslide surge are summarized. Firstly, the theoretical analysis progress of landslide surge in recent years is summarized. Secondly, the present situation of physical model test is summarized from four aspects: model, surge generation, surge propagation and surge run-up. Finally, the numerical simulation methods are summarized from the generation and propagation of surge. The third part discussed the influence factors of landslide surge. In the next part, the difficult problems and research status of landslide surge are evaluated. Finally, the study of landslide surge is prospected from the perspective of various factors of landslide surge.

1. Introduction

The vast territory of China benefits us with abundant material resources. Meanwhile, there are also various geological environments that bring a lot of geological problems. Among them, landslides are one of the common geological problems, which has been studied by domestic and foreign scholars over the years. With a lot of research carried out to deepen the research on landslide-induced issues, attention shifting towards landslide surges has been increased, and it is found that the damage caused by landslide surge is overwhelming [1,2].
The research on surges has always been a hot issue in the field of engineering geology [3,4]. In the beginning, it was associated with earthquake-induced tsunamis. For example, on 28 December 1908 in Messina, Italy, an earthquake with a magnitude of 7.5 triggered a tsunami, resulting in the death of more than 80,000 people. On 11 March 2011, a tsunami, induced by a magnitude 9.0 earthquake that occurred in the eastern sea area of Japan, inflicted huge casualties and economic losses. These cases demonstrate that the occurrence of tsunamis is uncertain and random, which often means devastating disaster.
A landslide surge is the surge developed from the vicinity water resources that are impacted by a landslide or collapse of masses, which threaten the security of buildings and people along its propagation path. Landslide surges mostly occur on the bank slope of a reservoir area. Compared to tsunamis, these landslide-prone areas have more buildings, people and infrastructure. As a result, landslide surges have a greater chance of causing property damage and human casualties. In addition, the development of the economy increased the construction of large reservoirs, especially for the continuous improvement of large water conservancy facilities in recent years, leaving the possibility of landslide surges. Therefore, extensive research has been carried out on the topic of landslide surges. Actually, there are many factors affecting landslide surges [5,6,7,8,9,10,11]. It is a complicated geological problem, which is not only influenced by complicated environment and large-scale solid landslide material but also complicated fluid material. To this end, it is challenging to study landslide surges.
The common research methods of landslide surges mainly include theoretical analysis, a physical model test and numerical simulation. Each method has its advantages in different aspects. In terms of theoretical analysis method, it can give an exact solution to a mathematical problem. Meanwhile in the study of the action mode of landslide surge damage, it can be demonstrated more clearly through physical model test. Compared with the model test, the numerical simulation method is low in cost and capable of adapting to the landslide surge damage of complex conditions, which satisfies multi-condition research and repeated simulation. On the basis of consulting a lot of literature, this paper uses real cases of landslide surges to illustrate the destruction of landslide surges; the research methods and progress of slope surge are summarized. Some difficulties and prospects in the study of landslide surge are prospected in the hope of providing some help for the follow-up study of landslide surges.

2. Landslide Surge Cases

The surges have caused awe in humans since ancient times. The records of surges can be sourced from more than a thousand years ago, in which there are many cases of landslide surges recorded hundreds of years ago. With the improvement and development of monitoring systems, more and more landslide surge cases have been recorded. Table 1 lists the major cases of landslide surge worldwide. As can be seen from Table 1, landslide surges are common all over the world. They usually have large landslide volumes and generate surges ranging from tens of meters to hundreds of meters. The landslide surges caused massive casualties.

2.1. Classification of Landslide Surges

Previous studies have concluded that the important reason for the formation of landslide surge waves can be attributed to the surface gravity wave, which is caused by the external force-induced disequilibrium between the fluid inertia and restoring force. The research about the surface gravity wave can be sourced from two kinds of literature. Herein, the first type of literature, represented by ‘Basic wave mechanics: for coastal and ocean engineers’ [27], focuses on the further research of water wave dynamics and wave theory, whilst the second kind of literature is represented by ‘fluid mechanics’ [28]. Fluid mechanics, a branch of mechanics, is the study of the related mechanical behavior of fluids (including gases, liquids and plasmas). Fluid mechanics can be divided into hydrostatics and fluid dynamics according to the motion mode of the object of study. The former studies the fluid at rest, and the latter studies the influence of force on fluid motion. Fluid mechanics according to the scope of application is divided into hydraulics and aerodynamics and so on. Fluid mechanics is a branch of continuum mechanics, which considers the characteristics of a system from a macroscopic perspective rather than a microscopic one [29]. Fluid mechanics mainly use nonlinear theoretical analysis and a modern number sequence method to describe the characteristics of water flow. The theory of water wave dynamics holds that the wave of landslide surges belongs to the category of gravity waves, which satisfies or roughly satisfies the constraint of wave theory formula. However, the wave property of landslide surges is ignored in fluid mechanics, which directly calculates the water particles.
According to the relative position of the landslide and water surface, the landslide surge waves can be classified into three types [30]: (1) landslide surges induced by the collapse masses above water surface, which involves the movement of solid, water and air, (2) landslide surges induced by collapse masses that are partially under water surface and (3) landslide surges induced by underwater collapse masses. Three positional relationships are shown in Figure 1.

2.2. Study on the Mechanism of Landslide Surge

According to the classification of the description parameter of wave theory, the characteristics of landslide surges are more complex. Most landslide surge waves are non-periodic waves and have strong nonlinearity, which are between moderate water waves and shallow water waves [31,32].
Through 211 groups of physical similarity tests, Valentin Heller [33] found that the height of surge wave is generally 5/4 times of the wave amplitude, and the wave velocity is generally close to the theoretical velocity of solitary waves.
Along the same line, Noda, Huber, Panizzo, Zweife et al. [34,35] found that different forms of surge waves are formed when the Froude coefficient ( F ) , the relative thickness of the block ( S ) , the relative mass ( M ) and the effective impact angle of the landslide ( β ) are different.
Panizzo, Zweifel, Yin et al. [36] carried out a large number of similarity tests of the landslide surge and found that the maximum height of surge is controlled by seven factors, including the water-entry velocity of the collapse masses of landslide, the thickness and width of the collapse masses, the volume and density of the collapse masses, the impact angle and the water depth. The test results also demonstrate that the wave height approximately satisfies a specific empirical formula. Besides, the wave height of the propagation wave is a more complicated problem. It is influenced by more factors, which are not only affected by the maximum height of surge but also the river surface and landform along the propagation path.
While Kamphuis et al. [9] deemed that the effect of sliding impact angles is negligible, Heller et al. [28] believed that wave height is related to cos ( 6 / 7 α ) and negatively correlated with α . These research results show that the influence multiple of different impact angles is varied between 0.5 and 1.95. However, in many physical similarity tests, velocity and sliding impact angle are interrelated, which leads to the difference of the above research conclusions. In addition, even different shapes of masses entering the water will induce different surges.
The study of Xu [37] used the coupled Eulerian–Lagrangian (CEL) algorithm and a variety of orthogonal numerical experiments and found that the influence of landslide shape on landslide surge, which is more significant within the near field. Besides, the volume of the landslide has a great influence on the height of the surge, and the influence decreases with the increase of the propagation distance.
Watts [38] pointed out that the deformation of the moving deformable masses is controlled by the motion of its centroid, and the difference between the surge generated by the deformable masses and the rigid body (with the same initial shape as the deformable masses) is inconspicuous.
The above research and investigation indicate that the occurrence of landslide surge is caused by the interaction of many factors. Due to the frequent occurrence of landslide surge in recent years and the complexity and diversified characteristics of landslide surge, many scholars have used a variety of research methods to study landslide surge from multiple angles, and the research on landslide surge has gradually become mature and perfect. The following is a summary of current landslide surge research methods from the perspective of different research methods.

3. Landslide Surge Research Methods

3.1. Theoretical Analysis of Landslide Surge

3.1.1. Surge Generation Stage

Compared with numerical simulation and physical model experiments, theoretical analysis can attain an exact solution of a closed idealized mathematical problem. Based on the energy conservation law, Noda [35] assumed that the reservoir area is a semi-infinite water body and the landslide body is a rigid body. In this case, the landslide body movement can be regarded as the motion of its centroid. Figure 2 is a schematic diagram of the landslide
Therefore, the water-entry velocity of landslide body can be expressed by Equation (1):
v = 1 cot α tan φ c l m g sin α 2 g H
where v is the water-entry velocity of landslide body; α is the sliding angle; φ , c is the shear strength parameter of sliding surface; m is the weight of the landslide body; H is the distance from the centroid of landslide body to the water surface; l is the length of the contact surface between the landslide body and sliding surface; g is the local gravity.
The study of Noda [35] simplified the landslide body into a rectangular rigid box, and the height H was set to be greater than the depth of dead water. This is to eliminate the impact of the height of landslide body on surge. The process of surge caused by the impact of rigid box is expressed as:
η ( x , t ) λ = 2 π × 0 d u u 0 T sin h u 1 s ( τ ) cos ( u x ) V ( τ ) d τ [ 1 s ( τ ) ] cos ( h u )
where η ( x , t ) denotes the elevation of the wavefront at the coordinate x and time t , λ is the width of the rigid box, T is the time duration for velocity to decrease from the water-entry velocity to 0, V is the velocity of the rigid box in water, s is the moving distance of the rigid box in water, σ = tan ( h u ) u , u = k h , k is the wave number, h is the depth of dead water and τ denotes an extremely short time interval. The limitation of this analytical solution is that the effects during the box impacting the free surface of water and the stage of underwater motion on the surge are not considered.
Based on Noda’s research, Di et al. [39] proposed a new analytical solution of initial surge under the framework of linear theory, which considered the limitation of Noda’s method.
By using the one-dimensional linear shallow water equation, Liu et al. [40] derived the analytical solution of the surge generated by landslide in two-dimensional model. The expression of wavefront elevation generated under the uniform sliding velocity of landslide body is:
η ( ξ , t ) = 1 3 h 0 ξ h 0 ξ + 0 ω a ( ω ) J 0 ( ω ξ ) cos ( ω t ) d ω + 0 ω b ( ω ) J 0 ( ω ξ ) sin ( ω t ) d ω
where J 0 is the zero-order Bessel function of the first kind; a ( ω ) and b ( ω ) are spectral functions determined by initial conditions; ω is wave frequency; ξ is a wavefront function.
Pan [35] combined the study of Noda [35] with the linear wave theory. This study divided the landslide into several blocks, and the interaction force between the blocks is assumed to be ignored; then, the landslide velocity of the first block is expressed as:
V x i = 2 a x ( i 1 ) Δ L + V x ( i 1 ) 2
where V x i is the landslide velocity of the i th block; a x ( i 1 ) is the landslide velocity of i blocks; Δ L is the width of the i th block.
Based on the momentum theorem, Wang et al. [30] analyzed the landslide surge by the slice method and obtained the calculated height of the initial surge by the momentum theorem:
η = η i + H i 2 h 0
η = m i v i H i B i t i ρ g
where η is the calculated height of the initial surge, m i is the weight of the landslide body, v i is the velocity of the landslide body, t i is the time duration of the landslide body moving underwater, H i is the height of the landslide body, B i is the width of the landslide and η i is the distance from the centroid of the landslide body to the highest point of the surge.
Based on the Noda method, Wang et al. [36] believed that the underwater movement of the landslide body has no impact on the water surface and would not cause a surge. Therefore, the H in the expression for calculating the water-entry velocity of the landslide body should be modified as H , which should be the distance between the centroid of the underwater part of the landslide body and the water surface. The expression is modified as:
v = 1 cot α tan φ c l m g sin α 2 g H
Dai et al. [19] deemed that the Noda method does not consider the resistance of water on the landslide body after entering the water. They consider the water resistance is R = 1 2 c w ρ f v 2 S , and the calculation expression of the initial sliding velocity after considering the water resistance is derived:
v = 2 H ( m g sin α m g cos α tan φ c L ) ( m sin α + H c W ρ f S )
Additionally, Dai et al. [19] also believed that the slice method proposed by Pan could not accurately describe the shape of the landslide surface, so they proposed using the motion equation method to calculate the initial sliding velocity of the landslide, and the acceleration calculation is improved by the motion equation method:
a i = { i k ( W i sin θ i μ i cos θ i ) + i k P w i [ cos ( α i θ i ) μ i sin ( α i θ i ) ] i n c i L i } / i = 1 n M i
where W i is the weight of the i th block; c i is the cohesive force of the slip band of the ith block; L i is the ground length of the i th block; θ i is the inclination angle of the i th block; μ i is the dynamic friction coefficient of the i th block; k is the number of blocks.
Some scholars have explored the generation mechanism of other waves except landslide surge wave through theoretical analysis [41,42]. Haugen et al. [43] studied the basic mechanism of tsunami generated by underwater mass flow under ideal geometric conditions through a two-dimensional linear analysis model. Ursell et al. [44] proposed the classical small amplitude wave theory of the piston generates wave, but this theory is only applicable to the small amplitude motion of plate. Kennard [45] developed a linear theory of waves, and Das et al. [46] determined that the linear theory proposed by Kennard has to satisfy the following three conditions: (1) the displacement of the vertical wall must be smaller than the dead water depth; (2) the acceleration of the wall must be less than the gravity; (3) the sliding constant must be less than unit 1.

3.1.2. Propagation Stage

By assuming that the coastline is a parallel steep face, Pan [47] ignored the energy attenuation in the process of surge propagation and the nonlinear influence of boundary conditions. The surge propagation processes are regarded as the linear superposition of a series of wavelets generated from the source point, where each wavelet is a solitary wave and has a total reflection on the opposite bank of the reservoir. Assuming that the reflection coefficient is known, the calculation expression of landslide surge height is obtained:
ζ = 2 ζ 0 π ( 1 + k ) × n = 1 , 3 , 5 n { k 2 ( n 1 ) × ln l ( 2 n 1 ) B + 1 + l 2 ( 2 n 1 ) 2 B 2 }
where ζ is the surge height on the opposite bank of landslide; ζ 0 is the initial surge height on the opposite bank of the landslide; k is the reflection coefficient of surge; B is the width of the river; n are the times of reflections. The calculation expression of landslide surge height at downstream dam site is expressed as:
ξ = ξ 0 π × n = 1 , 3 , 5 n ( 1 + k cos θ n ) k n 1 × ln 1 + ( n B x 0 L ) 2 1 x 0 x 0 L 1 + ( n B x 0 ) 2 1
where ξ is the surge height on the opposite bank of landslide; ξ 0 is the initial surge height at the downstream dam site; x 0 is the distance from the centroid of the landslide to the downstream dam site; θ n is the intersection angle between the nth incident line and normal line of bank slope.
Wang et al. [30] studied the attenuation of surge in the propagation process and deemed that the propagation stage of surge in the river is divided into a sharp attenuation stage and a slow attenuation stage. Based on the continuity equation and motion equation of unsteady flow of open channel, the calculation expression of landslide surge in the sharp attenuation stage is obtained:
h s ( x , t ) = h 1 e k 1 x g h j t
where h s is the height of surge with a distance of x from the landslide occurrence point, h 1 is the wave height and h j is the initial surge height; k 1 = T 0 a 2 g k T 0 a A 0 g , T 0 is the initial width of the water surface; x is the ratio of actual propagation distance to water depth; t is propagation time in sharp attenuation stage.
The calculation expression of landslide surge in slow attenuation stage is:
h s ( x ) = h 1 2 g ( h 0 + h 1 ) 2 n 2 ( 2 h 0 + h 1 ) R 4 3 x
where h 1 is the height of surge with a distance of x from the landslide occurrence point; h 0 is water depth; n is the roughness coefficient; x is the propagation distance in the slow attenuation stage; R is the hydraulic radius.
Li et al. [48] regarded the propagation process of landslide surge as a process of the interaction between a shock wave and a rarefaction wave. The interaction between the shock wave and rarefaction wave is qualitatively obtained by adopting the shallow water wave equation, mass coordinate and velocity wave height curve method:
When
u r < u l g ( h r + h l ) ( h r h l ) 2 2 h r h l
an incident shock wave will be generated.
When
u r > u l g ( h r + h l ) ( h r h l ) 2 2 h r h l
an incident sparse wave is generated.
u r and h r are, respectively, the right side flow velocity and water depth before the interaction between rarefaction wave and shock wave, while u l and h l denote the flow velocity and water depth on the left side before the interaction. According to the mass conservation equation, Li calculated the wave velocity of the shock wave produced by the interaction between the rarefaction wave and the shock wave:
D = u 0 + g h 0 1 + h 1 / h 0 2
where D is the velocity of shock wave; h 0 is the water depth of shock wave front; u 0 is the flow velocity of the wave front of shock wave; h 1 is the water depth of the wave back of shock wave; u 1 is the flow velocity of the wave back of shock wave.
Based on the study of Pan, Ha et al. [49] considered the phase difference of wave superposition during propagation, and the calculation expression of surge height at the opposite bank and dam site of landslide reservoir area was modified:
ζ = ( 1 + K ) ζ ζ = ( 1 + K cos θ ) ζ
where ζ is the wave height before reflection, 2 is the reflection wave height, 3 is the secondary reflection wave height, K is the reflection coefficient, and θ is the incident angle.
From the perspective of energy, Law [50] made a qualitative study on the propagation process of the landslide surge. The study found that the energy of the initial wave accounts for 98% of the total energy of the landslide surge, and the attenuation of the initial wave is inversely proportional to the square of the propagation distance along the propagation path.
Meanwhile, Hunt [51] made a judgment on the selection of the calculation theory of landslide surge: when the landslide occurs rapidly and the transverse size of the underwater part of the landslide body does not exceed twice of the water depth, the generated surge is short-wave. If the transverse size of the landslide body exceeds 10 times of the water depth, the generated surge wave is long-wave.

3.2. Physical Model Test Method

According to the development of the landslide surge, the evolution process can be divided into three stages, namely, the generation stage, the propagation and attenuation stage and the climbing stage. Figure 3 shows the stages of landslide surge.

3.2.1. Physical Model

Landslide surge is affected by many factors such as landslide body, water body and terrain. Therefore, there are certain advantages to generalizing the landslide body and water body into a simple physical model. However, this may impact the reliability of test results. Domestic and foreign scholars have optimized the physical model method constantly to attain the robust reflection of the real damage.
In 1955, Wiegel [52] carried out physical model tests by adopting steel plates to simulate landslide blocks.
Noda [35] and Pan [47] focused on the surge induced by horizontal and vertical landslides. The model was generalized into a cuboid thereby ignoring the influence of the shape of the landslide.
Ashtiani [53] considered the influence of landslide shape by using sand-gravel particle material to make granular model.
Considering that the particle model can only simulate the loose particle landslide, Yin et al. [5] applied the river physical model to simulate the water surface, while the granular materials and cement were adopted to fabricate different sizes of landslide body. Meanwhile, in order to reconstruct the real process of landslide surge and fully consider the complexity of landslide surge, the landslide inclination and water-entry speed were controlled by means of the test control system, which significantly improved the flexibility of the test.
Other scholars have used different materials to simulate slides. For example, Han et al. [6] used multiple rigid blocks with sizes controlled by structural planes in the three directions (length, width and thickness) to simulate the landslide body. Viroulet et al. [2] used glass beads of different sizes and shapes to simulate the landslide body.
With the deepening of the research on the simulation of physical model test slope, some scholars consider the similarity between the test model and the prototype. Huang et al. [8] investigated the Gongjiafang landslide. The marble rock which is high in similarity in terms of density with the main compositions of the prototype landslide was adopted to simulate the landslide body, in order to ensure the similarity principle. Wang et al. [11] simulated the water body by using the river generalized model and used concrete, steel bars, plastic plates, plastic pipes and other materials to simulate the wharf in Wanzhou Port. The parameters including the hull shape, full load displacement and draft depth were designed to meet the similarity relationship.
Some scholars not only simulated the disaster factors but also restored the disaster process of the carrier when landslide surge occurred. Chen et al. [18] studied the safety of mooring ships under landslide surge by replicating the river section of Tuokou wharf in Jiangnan, Wanzhou, and a deck barge of 3000 t was selected as a test ship. In addition, some scholars designed a pneumatic landslide-generating device to simulate landslide surge, which can greatly improve the test repeatability and accurately, and all the important parameters of surge can be controlled.

3.2.2. Surge Generation Stage

As an important part of the landslide surge disaster, the generation stage of surge has been studied by many scholars through physical model tests.
Noda [35] carried out theoretical derivation through physical model test and first proposed the initial wave equation under the impact of horizontal landslide:
H max h = 1.17 v g h
Based on Noda’s research, Pan [42] continued to carry out physical model tests and proposed the initial wave calculation equation of vertical landslide surge.
H max = v g h , ( 0 < v g h < 0.5 ) f v g h , ( 0 < v g h < 2 ) 1 , ( v g h > 2 )
Wang et al. [11] studied the landslide surge of mountain channel reservoir by a 1:70 scale model test and proposed the equation of total wave energy of initial wave:
E w = P T S = 1 8 ρ g H 2 L S n
where
n = 1 2 1 + 2 k h sin ( 2 k h )
H is the initial wave height, S is the length of the peak line of the initial wave, c is the wave velocity, and L is the wavelength. When the water level is high, n = 1/2; when the water level is low, n 1 .
Ke et al. [54] carried out physical model tests based on a river of Wushan County and Fengjie County. The correction coefficient K shown in Equation (22) was proposed. Then, the correction coefficient was validated by taking the Gongjiafang landslide as an example.
K = 2.30 ( v s g h ) 0.93 ( s b h 2 ) 0.43
Through the physical model test, Peng et al. [50] designed an orthogonal test scheme comprised of six factors which affect the landslide surge. The results demonstrated that the water level of the landslide body entering the water and the thickness of the landslide body are important factors to control the initial wave height of the landslide surge. The accordingly fitting empirical equation of the initial wave is obtained:
s h 0.462 × tan α 0.379
Yin et al. [5] transformed the experimental data into dimensionless form and optimized the maximum initial wave equation of landslide obtained by a previous study which is expressed as Equation (24):
H max h = 1.17 v g h ( sin 2 α + m cos 2 α ) g 1 ( l t b h ) g 2 ( w b )
The equation for calculating wave propagation was obtained by the multivariate nonlinear regression method:
H p h = 1.47 H max h ( x h ) 0.5 , ( x h > 2.13 )
While the equation for calculating wave at climbing stage is
R h = 2.3 H c h ( 90 β ) 0.2
Kamphuis et al. [7] conducted a landslide surge test in a two-dimensional flume, in which the angle of the landslide was controlled between 20° and 90°. Besides, the dimensionless relation of initial wave height is determined by dimensional analysis:
H h = f ( F r , L , S , α , φ , n , D , X , T )
where F r is the sliding Froude number; L is the relative landslide length; S is the relative landslide thickness; α is the sliding dip angle; φ is the angle of the front of landslide; n is the porosity of landslide; D is relative landslide density; x is relative propagation distance; T is relative time. Herein, the influence of F r , L and S is the most conspicuous.
Walder et al. [55] studied the near-field surge generation area through a two-dimensional flume model, in which the motion characteristics of the landslide model were recorded, and the dimensionless expression of the maximum amplitude of the near-field surge was accordingly derived:
a m h = 1.32 ( T s V ) 0.68
where T s is the time duration of the relative underwater motion of the landslide; V is the relative landslide volume, and a m is the maximum amplitude.
Based on the momentum balance theory, Han et al. [6] derived the theoretical expression of the maximum near-field amplitude of surge under hydrostatic and fluctuating conditions:
a m 1 = h 2 + 2 ρ s s b v s 2 cos α L ρ g K h
a m 2 = 2 ρ s s b v s 2 cos α ρ g h T ( s + b / 2 )
where a m 1 and a m 2 are the maximum amplitude under hydrostatic conditions and fluctuating conditions, respectively; b and s are the width and thickness of the landslide; v s is the velocity of the landslide along the slope direction; h is the dead water depth; α is the angle of landslide body; ρ s is the density of the landslide body; ρ is the density of the water body, and T is elliptic coefficient.
Huang et al. [56] established a 1:150 experimental model based on the Lechangxia reservoir. The characteristic of landslide surge with respect to different water-entry velocities and water level conditions was investigated, and the test results were compared with the theoretical results of Pan Jiazheng method, the method recommended by the American Civil Engineering Society and the empirical equation method of the Academy of Water Sciences. The results demonstrated that the Pan Jiazheng method is more accurate in the case of applying the empirical equation method to predict the surge.
Huang et al. studied the influence of water entry velocity and water level conditions on the generation of surge wave but ignored the influence of landslide volume on the generation of surge wave. Therefore, some scholars have explored the influence of landslide thickness and width on the generation and propagation of surge waves. Yuan et al. [3] studied the generation and propagation law of surges subjected to landslides with different widths and thicknesses within the curved channel. The results illustrated that the surge height is proportional to the width and thickness of landslide.
In order to explore the influence of more factors on landslide surge wave, the landslide surge wave height formula was obtained. Based on the channel model of a hydropower station on the Lancang River, Ding et al. [57] used the orthogonal test method to carry out the landslide surge test. The effects of influencing factors such as landslide size, water-entry velocity, water depth and river width on the maximum initial wave height, as well as the regularity of wave propagation, were studied. The equation for calculating the maximum initial wave and the propagating wave were fitted by the method of multiple regression. In addition, Li et al. [58] took the reservoir area near the dam as the research object to analyze the variation law of surge height under different landslide factors such as shape, water depth, volume and falling height in the narrow channel. He concluded that the height of the surge is proportional to the sliding height and the volume of the initial drainage and inversely proportional to the water depth. Xiao et al. [59] established a 1:200 river channel model to study the influence of different factors on the characteristics of the water tongue and the height of the initial wave in the near-field region. The empirical equation for calculating the height of the water tongue and the length of the water tongue were, respectively, obtained by the fitting variance method and the regression method. The reliability of the equations was verified by the cases of the Dayantang landslide and the Qianjiangping landslide.
In order to explore the damage pattern of surge in different scenes and bank slope shapes, McFall [60] studied the characteristics of surge in more scenarios based on model experiments, including open seas, narrow fiords, curved straits and conical islands. The results show that the initial wave crest amplitude generated by the landslide along the plane is greater than that along the conical island slope. However, the amplitude of the initial wave trough and the second wave crest is smaller than that of the conical island. In addition, the energy conversion rate of the landslide kinetic energy was within the range of 1% and 24%, and the energy conversion rate of pebble landslide was 43% larger than that of gravel landslide on average.
In addition, some scholars have adopted some special methods to make the physical model test of landslide surge wave better reflect the real damage situation. Mohammed et al. [61] conducted a series of granular landslide surge experiments in a three-dimensional wave basin. The model landslide body adopted natural circular river gravel with a particle size of 12.7–19 mm. The landslide movement was controlled by a new pneumatic landslide surge generator fixed on a slope. Based on the experimental results, the empirical equation of near-field wave was obtained by using multivariate regression analysis method, which is mainly related to sliding Froude number, relative landslide thickness, relative landslide width and relative landslide length.
Some scholars design experiments based on some field investigation data to propose risk analysis methods. Liu et al. [62] established a 1:200 physical model based on the Baishui River landslide channel, and the geological data of more than 100 potential landslides in the Three Gorges Reservoir area were statistically analyzed. Generally, the three-dimensional physical model test was performed within the scheme of the orthogonal test. The velocity calculation model of wading landslide after instability was accordingly proposed, which was fed as the calculation model for the landslide surge in the Three Gorges Reservoir area after the sensitivity analysis method. Heller et al. [31] discussed the scale effect in the landslide surge model experiment and deemed that when the water depth of the model is less than 0.2 m, the scale effect caused by fluid viscosity and surface tension cannot be ignored. Besides, the landslide surge tests with respect to four different granular materials were conducted by using the Fritz experimental device. In the experiment, the deformation of the granular landslide and the velocity of its centroid before entering water was captured by the two laser distance sensors arranged above the slide. By analyzing the maximum surge in the near field, it was found that the surge wave height is about 5/4 times of the wave crest amplitude, indicating that the nonlinearity of the near field surge is conspicuous, and the wave velocity can be approximately described by the solitary wave velocity. Moreover, the test also demonstrated that the empirical equation of surge characteristic parameters caused by granular landslide is related to a dimensionless parameter P, which is called surge generation parameter.
P = F r S 1 / 2 M 1 / 4 ( cos β ) 1 / 2
In this formula, F r is the sliding Froude number; S is the relative landslide thickness; M is the relative landslide mass, β = 6 / 7 α ; α is the angle of landslide.

3.2.3. Propagation Process Stage

After the landslide body enters the water for a period of time, due to the propagation of surge energy in space, frequency diffusion and nonlinear effects, the surge wave height will gradually decrease as the propagation distance increases [63]. The purpose of studying the propagation stage of landslide surge wave by physical model test is to analyze the attenuation rule and put forward the empirical equation through the experimental data. Li et al. [64] established a flume experiment model based on the Three Gorges Reservoir, and the variation law of initial surge height subjected to different angles of landslide slope, solid-fluid effective contact surfaces and water depths was investigated. The regression analysis was applied to analyze the link between the wave height and the three factors. The empirical equation of relative initial surge height and wave height of different landslide water-entry points was, respectively, fitted, and the data error verification analysis was accordingly carried out. Finally, two empirical equations for calculating surge height along the propagation path were proposed, which provides a basis for calculating surge height at different positions. Yang et al. [65] studied the attenuation law of the initial surge height along the river channel by taking the steep rock landslide in the Three Gorges Reservoir area as the research object through model test and summarized the wave height attenuation relationship. The wave height attenuation coefficient was divided into two stages by the propagation distance of 3 m. The empirical equations for calculating the wave height attenuation coefficient at different stages were proposed. Finally, the applicability of the empirical equations was verified by the experimental results.

3.2.4. Surge Climbing Stage

The study of surge climbing process through physical model test is mainly aimed at the factors which affecting the climbing process, such as bank slope angle, bank slope permeability and so on. In response to these factors, scholars have adopted many different wave-making methods to generate surge waves.
Hall and Watts [66] simulated the generation of solitary waves through a piston wavemaker at the Waterway Experiment Station (WES) of the U.S. Army Corps of Engineers and studied the climbing process of solitary waves on impervious bank slopes. Through experiments, it was found that when β   <   12 ° ( β is the bank slope angle) the maximum wave height is mainly affected by friction, while when β   <   12 ° , the maximum wave height is mainly affected by gravity and inertial force. The empirical formula of the maximum solitary wave height was obtained by fitting the experimental data:
r h = 11 β 0.67 ( H h ) 1.9 β 0.35 ( 5 ° β   <   12 ° ) r h = 3.05 β 0.13 ( H h ) 1.15 β 0.02 ( 12 ° β 45 ° )
Synolakis [67] proposed an approximate theory of unbroken solitary wave propagation under the condition of constant water depth. The model test was carried out at a slope of 1:19.85, and the theoretical derivation results were compared with the experimental data. It was found that the theoretical derivation results were basically consistent with the experimental data when the solitary wave was not broken. The unbroken solitary wave climbing equation is called the “climbing law”.
r h = 2.831 ( cot β ) 1 / 2 ( H h ) 5 / 4
Müller [68] used Russell’s wave-making method to generate surges through falling heavy objects. The process of surge climbing and dam overflow with respect to slopes of 18.4°, 45° and 90° was, respectively, investigated. The prediction equation of surge climbing is expressed as:
r h = 1.25 ( 90 ° β ) 0.2 ( H h ) 1.25 ( H L ) 0.15

3.3. Numerical Method

3.3.1. Surge Generation Stage

The study of surge generation involves the interaction between landslides, reservoir water and air. Therefore, it is an important problem to simulate shock generation in numerical simulation of landslide surge [69,70].
The tracking of free surface is an important research topic in the simulation of surge generation. Accurately capturing the free surface of the fluid can truly simulate the generation process of surge. At present, the VOF method, SPH method, MAC method and LEVEL SET method have been applied to the study of this problem and have achieved good results. The LEVELSET method is more convenient to solve and easier to deal with complex interfaces and dramatic changes in free surfaces. Some scholars use this method to study the influence of surge and verify the applicability of CEL method in the study of landslide surge by comparing with the experimental results. However, the research is mainly aimed at a single block landslide and does not consider the granular landslide. The smooth particle hydrodynamics method (SPH) was originally applied to the calculation of irregular motion in three-dimensional space. Compared with the finite difference analysis method, SPH uses particle to represent the medium fluid to estimate the partial differential equation based on multiple variables [71].
< f ( x ) > = E f ( x ) W ( x x , h ) d x
In the formula, < f ( x ) > denotes the kernel estimate at coordinates; f ( x ) is the field value; W ( x , y ) is the kernel function. In recent years, more and more scholars have applied the SPH method to the field of fluid mechanics, which is suitable for simulating large surges caused by high-speed landslides.
The VOF method has strong advantages in dealing with free boundary problems and has good adaptability to large deformation. At the same time, it requires a small amount of storage, which is widely used in research.
Li et al. [72] used the finite volume method to establish the calculation model and discussed the mechanism of surge generation. The fluid calculation software FLUENT combined with dynamic mesh technology and VOF method was adopted to compare and analyze the surge height and pressure field generated by different landslide models.
Based on the Navier–Stokes equations, Song et al. [73] used the VOF method to simulate the generation process of surges. The simulation results were in good agreement with the experimental results and the observation results. Yuk et al. [74] used the VOF method to simulate the process of landslide surge generation, and the free surface was effectively tracked. Abadie et al. [75] described the air, water and landslide with Newtonian fluid and Navier–Stokes equation and simulated the surge by combining them with the VOF method. This research pointed out that the air also had a significant influence on the motion of landslides, but the research works of Abadie et al. were limited to rigid landslides.
Liu et al. [76] used the VOF method and LES model to study the phenomenon of water climbing and falling caused by the sliding of three-dimensional regular blocks. Based on the shallow water wave equation, the generation and propagation of landslide surges under a variable grid were simulated by using dynamic grid technology under the premise of known landslide movement. The results were compared with the experimental data, which showed that the dynamic grid technology can simulate the surge accurately.

3.3.2. Propagation Process Stage

There have been a lot of research results on surge propagation simulation, and most of these numerical models can directly be applied in research. In terms of the research method, the empirical equations have a certain impact on the accuracy of the calculation results since they are obtained based on idealized assumptions. Therefore, the numerical calculation can simulate the surge propagation more accurately than the empirical equations.
The numerical models of surge propagation mainly include the Navier–Stokes equation model, Boussinesq equation model, the model based on shallow water wave equation, etc.
Some scholars applied the N-S equation of fluid mechanics and established finite element and finite difference methods to analyze specific landslides; the height of the surge was obtained, and the law of energy exchange between landslides and water bodies was studied [77]. Ren et al. [78] derived the DIF equation for simulating landslide surge and used the irregular grid finite product method and the explicit MacCormack prediction-correction numerical method to solve the equation. The numerical model of landslide surge was established and verified by the relevant data of Xintan landslide. Xu et al. [79] adopted the Navier–Stokes equation model in the FLOW3D software and the VOF method to simulate the generation and propagation of surges in three-dimensional wide waters. Compared with the experimental results, it is indicated that the numerical model can well simulate the generation and propagation of surges, but it ignored the influence of slider form on surges, and the simulation accuracy needs to be improved.
By adopting the shallow water equation, Bosa et al. [80] established a 2DH numerical model to simulate the surge caused by the Vajont reservoir landslide. The results show that the shock wave caused by the large landslide can be simulated by the two-dimensional model.
Based on the two-dimensional fourth-order Boussinesq equation model, Ashtiani et al. [81] simulated the surge height, wave climbing and maximum wave height of Maku and Shafa-Roud reservoir areas, which provided a reference for further evaluating the surge disaster caused by potential landslides in the reservoir area.
In addition, through the secondary development of GEO-WAVE software, Wang et al. [82] established FAST system to study the propagation and attenuation law of surge and established a surge warning system, which has a good application prospect. Huang et al. [15] used GEOFLOW to simulate the propagation process of surges and discussed the propagation law of surges. Heidarzadeh et al. [83] used the TUNAMI-N2 numerical model to simulate the propagation of tsunami caused by submarine landslides.
In order to verify the accuracy of the numerical model, some scholars use the landslide prototype to verify. Jiang et al. [84] established a numerical model of landslide surge by solving the DIF equation and verified it by the Xintan landslide. However, it assumed that the landslide was a single block with a fixed shape, which ignored the impact of the shape of rock landslide on the surge. Shen et al. [85] simulated the propagation process of elliptic cosine waves, in which the VOF method was used to deal with the problem of wave free surface, and the simulation results were consistent with the experimental results.

4. Study on Influence Factors and Mechanism of Landslide Surge

The formation of landslide surge is a multidisciplinary problem involving landslide dynamics, rock and soil mechanics and fluid mechanics. Many scholars use numerical simulation, a physical model test and other methods to study the factors affecting landslide surge [86], and some patterns were found, including the following research results.

4.1. Sliding Body Factor

The study of sliding body factors mainly involves the landslide shape, volume of the sliding body and friction angle of the sliding surface. Through the study of the relationship between landslide shape and surge, it is found that the surge formed by wedge-shaped sliding mass is the largest, subsequently the rectangle, while the surge provoked by the oval block is the smallest. The surges induced by the rectangular and wedge-shaped sliding masses are generally the same when the water-entry point is far, but they are larger than those induced by elliptical blocks. Additionally, through the study of the relationship between the volume of the sliding mass and the surge, it is found that the volume of the sliding mass has a great influence on the surge height, especially in the near field condition. Scholars believe that the volume of the landslide is the most important factor controlling surge size. Through the study of the relationship between the friction angle of the sliding surface and the surge, it is found that the internal friction angle of the sliding surface has a great influence on the surge height under the same height. The friction angle of sliding surface affects the surge size by changing the velocity of sliding body entering water [38].

4.2. Water Body Factor

Among the water factors, the main factors that have a great influence are the width of the water surface and the distance between the sliding body and the water. Through the study of the relationship between water surface width and surge, it is found that when the surge reaches the opposite bank, the effect of the opposite bank on the surge gradually increases due to the blocking effect. With the decrease of water surface width, the surge height shows a slight increase trend, especially the climbing height. While studying the height of a landslide in water, it was found that with the increase of the distance from the water-entry point, the influence of the front edge shape and the volume of the sliding on the landslide surge gradually decreases [47].

5. Difficulties and Problems of Landslide Surge Research

Although there are many research methods on landslide surge, there are still many problems and difficulties. This section summarizes the difficulties in the current research on landslide surge from three perspectives: theoretical analysis, a physical model test and numerical simulation.

5.1. Theoretical Analysis Method

The theoretical analysis method is a simple and fast surge estimation method based on mathematics and mechanics. The selection of equations is restrained by the landslide failure mode and the water body. Due to the complexity of the surge, each calculation method needs to meet certain assumptions. For example, the foreign empirical calculation method assumes that surge is a one-way flow. In practical engineering, except for a small part of rock and soil mass falling vertically, most of the unstable sliding bodies were oblique sliding along the coastline, and the fluctuation of waves is mostly dynamic, which cannot be linearly superimposed. In addition, in order to simplify the calculation equation, some researchers assume that the problem of landslide surge is a two-dimensional plane problem, ignoring the influence of landslide thickness and scale on landslide surge. Furthermore, some researchers introduce the concept of centroid or rigid body and regard the water-entry rock soil as a complete block. However, in the actual process, these landslide bodies are often deformed or even broken, ignoring the influence of internal parameters will cause large errors in the prediction of some landslide surges. Moreover, the interaction between sliding body and water body is complex, which involves the fluid mechanics, rock mechanics and kinematics, while the influence factors of resistance such as fluid viscous force are often insufficiently considered. Therefore, it is considered that the continuity of surge calculated by empirical equation method is poor, which cannot well predict the influence of landslide failure process and wave evolution on surge disaster. Although the research on landslide surges has a long history, there are still some difficulties and problems.

5.2. Physical Model Test Method

In order to simulate the complexity of geological problems and solve the problem of theoretical analysis, scholars use a physical model test to study problems, but the physical model test method also has some difficulties and problems. The physical simulation method is based on the similarity theory to replicate the landslide body and water body in scaled dimension. This method can directly observe the wave body shape and the whole landslide surge formation, propagation and even climbing process. It is a relatively intuitive and reliable research method. However, the preparation period is long and requires a lot of labor and material resources. Mostly, the physical model test is generally aimed at a certain case, so it is only applicable to a single case or phenomenon, resulting in poor versatility. Meanwhile, there are many influencing factors of landslide surge; to summarize the variation law from the model test is the focus and difficulty of the research, and the influence on the size effect is not clear at present. Therefore, determining the appropriate scale is also a key point of the model test.

5.3. Numerical Method

The numerical simulation method can solve the problems of poor generality and the high cost of the physical model test. The numerical analysis method is a kind of flexible and intelligent solution method based on the continuous development of mathematics and computer science, which is low in cost and easy to control and can better replicate the whole process of landslide surge. In the early time, the numerical calculation is mostly for simplifying the complex multi-phase fluid–solid coupling problem into the unidirectional channel flow problem after determining the scale of the landslide body and the water-entry speed and then using the Saint-Venant equation to determine the boundary conditions, so as to calculate the simulation results. However, assumptions are too idealized to attain precise results. With the development of computational fluid dynamics, the numerical analysis method begins to be applied in solving the mathematical problems of initial conditions and surrounding boundary conditions. However, due to the complexity of the actual situation, researchers often question the accuracy of the obtained initial conditions and surrounding boundary conditions. Therefore, the numerical simulation method is mostly used as an auxiliary fitting method, which needs to be verified by other methods.

6. Conclusions

At present, the study of landslide surge cannot clearly include all factors. In the case of complex influencing factors, different research methods have certain shortcomings. This section summarizes and forecasts the landslide surge from three perspectives: theoretical analysis, a physical model test and numerical simulation.

6.1. Theoretical Analysis Method

The generation stage of surge is essentially the process of transforming the potential energy of the landslide into wave energy [87]. A landslide surge is highly nonlinear, high-order and complex. The energy conversion rate is not only related to the initial position of the landslide but also to the geometric shape, motion state, water boundary conditions, water depth and other factors of the landslide. Energy conversion rate affects the magnitude and scope of the disaster. At present, there is little research on energy conversion, and more research on landslide surges is needed in the future. In addition, the propagation and attenuation of landslide surge are hotspots of current research. However, it mainly focuses on the spatial attenuation law of wave amplitude or wave height along the radial direction at present and rarely involves the change of surge with time. In fact, the attenuation stage of surge has both space attenuation and time attenuation. The study of the spatial attenuation law is helpful to determine the spatial range of the disaster. The attenuation law in time dimension can provide theoretical guidance for early warning and prevention of the disaster.

6.2. Physical Model Test Method

At present, the physical model test of landslide surge mainly simulates block landslide and granular landslide, but there are a lot of rock landslides in the reservoir slope. The broken rock mass is neither a block model nor a granular model, so it is necessary to study this type of landslide surge. In addition, Through the analysis of the physical model of the landslide, it can be seen that the current research mainly considers the density, saturation and water content of sandstone to design the landslide model and pays less attention to mudstone. However, the particle size of mudstone is relatively small, and the water absorption is good; it is easy to soften in water and become discrete or break during the sliding process. Therefore, more attention should be paid to the physical model test of the mudstone landslide in the future.

6.3. Numerical Method

Although the development of computing technology and the computer has made it possible to use the numerical simulation technology of the N-S equation to study the landslide wave, the complex model and the decomposition, collision and deformation of complex rock and soil mass cannot be fully considered. At the same time, using the numerical simulation technology of the N-S equation to simulate requires too much computer performance and time. Therefore, a breakthrough should be made in the calculation equation and computer performance of the numerical simulation of landslide surge in the future.

Author Contributions

C.X. wrote the original draft preparation; Z.H. collected the literature related to swells; W.F. edited and corrected the article. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the National Natural Science Foundation of China (Grant No. 52108361 and 41977252); Sichuan Provincial Department of Science and Technology key R & D project (Grant No. 2023YFS0436); Sichuan Provincial Department of Science and Technology project (Grant No. 2023NSFSC0376).

Data Availability Statement

The original data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Acknowledgments

We acknowledge the anonymous reviewers for their comments.

Conflicts of Interest

The authors declare that they have no known competing financial interest or personal relationships that could have appeared to influence the work reported in this paper.

References

  1. Wu, S. Theory and Technology of Landslide Risk Assessment; Science Press: Beijing, China, 2012; pp. 8–75. [Google Scholar]
  2. Li, H.; Tan, Y.; Li, E. Research on quantitative of landslide assessment disaster risk in reservoir area of Kala hydropower station. Chin. J. Undergr. Space Eng. 2013, 9 (Suppl. S2), 2040–2046. [Google Scholar]
  3. Yuan, P.Y.; Wang, P.Y.; Zhao, Y.; Wang, M.L. Experimental study on the nonlinear behavior of a sailing container ship under landslide induced surges. Adv. Civ. Eng. 2019, 2019, 9081586. [Google Scholar] [CrossRef]
  4. Wang, P.; Han, L.; Yu, T.; Meng, C. Effects of landslide generated impulse waves on ship impact force for pile wharf. J. Harbin Eng. Univ. 2016, 37, 878–884. [Google Scholar]
  5. Yin, K.; Liu, Y.; Wang, Y.; Jiang, Z. Physical model experiments of landslide-induced surge in Three Gorge Reservoir. Earth Sci. (J. China Univ. Geosci.) 2012, 37, 1067–1074. [Google Scholar]
  6. Han, L.; Wang, P. Prediction of the maximum near-field wave amplitude of impulse waves generated by three dimensional landslides based on momentum balance. Chin. J. Rock Mech. Eng. 2018, 37, 2584–2592. [Google Scholar]
  7. Viroulet, S.; Sauret, A.; Kimmoun, O. Tsunami generated by a granular collapse down a rough inclined plane. EPL (Europhys. Lett.) 2014, 105, 34004. [Google Scholar] [CrossRef]
  8. Huang, B.; Wang, S.; Chen, X.; Yin, Y.; Jiang, Z. Prototype physical similarity experimental study on surgegenerated by instability of cataclastic rock mass. Chin. J. Rock Mech. Eng. 2013, 32, 1417–1425. [Google Scholar]
  9. Zweifel, A.; Hager, W.H.; Minor, H.E. Plane impulse waves in reservoirs. J. Waterw. Port Coast. Ocean Eng. 2006, 132, 358–368. [Google Scholar] [CrossRef]
  10. Yue, S.; Diao, M.; Wang, L. Research on initial formation and attenuation of landslide-generated waves. J. Hydraul. Eng. 2016, 47, 816–825. [Google Scholar]
  11. Wang, M.; Zu, F.; Wang, P.; Han, L. Head Wave Energy Analysis of Landslide Surge in Mountainous Channel Reservoir. Waterw. Eng. 2020, 143, 79–83. [Google Scholar]
  12. Wang, J.; Wang, L.; Chen, X. Calculation of lenggu hydroelectric station landslide surge height based on pan jiazheng calculation method of speed and surge. Water Resour. Power 2010, 9, 95–98. [Google Scholar]
  13. Wang, Y.; Liu, J.; Zhang, Y.; Yin, K.; Huo, Z. Review of wave amplitude prediction generatedbylandslide based on physical experiments. Geol. Miner. Resour. South China 2018, 34, 279–288. [Google Scholar]
  14. Yin, K.; Du, J.; Wang, Y. Analysis of surge triggered by Dayantang landslide in Shuibuya reservoir of Qingjiang river. Rock Soil Mech. 2008, 29, 3266–3270. [Google Scholar]
  15. Dai, Y.; Yin, K.; Wang, Y. Discussion on method of landslide velocity calculation and surge prediction. Rock Soil Mech. 2008, 29, 407–411. [Google Scholar]
  16. Du, B. Tangyanguang landslide of tuoxi reservior-first large landslide induced by reservoir storage in China. In Proceedings of the Symposium of the Second Conference of Chinese Geotechnical and Engineering, Wuhan, China, 28–31 October 2006; China Science and Technology Press: Beijing, China, 2006; Volume 1, pp. 918–922. [Google Scholar]
  17. Zhong, L. Enlightenment from vajont landslide in Italy. Chin. J. Geol. Hazard Control 1993, 5, 77–84. [Google Scholar]
  18. Hermann, M.F. Initial Phase of Landslide Generated Impulse Waves; VAW Publikationen: Zurich, Switzerland, 2002. [Google Scholar]
  19. Chen, L. Destroying tsunamis generated by landslides in history. Yangtze River 1984, 2, 88. (In Chinese) [Google Scholar]
  20. Huang, R. Typical Catastrophic Landslide in China; Science Press: Beijing, China, 2008. [Google Scholar]
  21. Huang, B.; Yin, Y.; Liu, G.; Wang, S.; Chen, X.; Huo, Z. Analysis of waves generated by Gongjiafang landslide in Wu Gorge, three Gorges reservoir, on 23 November 2008. Landslides 2012, 9, 395–405. [Google Scholar] [CrossRef]
  22. Liao, Q.; Li, X.; Li, S.; Dong, Y. Occurrence, geology and geomorphology of Qianjiangping landslide in Three Gorges Reservoir area Study on characteristics, cause and landslide criterion. J. Rock Mech. Eng. 2005, 17, 3146–3153. [Google Scholar]
  23. Huang, B.; Yin, Y.; Du, C. Risk management study on impulse waves generated by Hongyanzi landslide in Three Gorges Reservoir of China on 24 June 2015. Landslides 2016, 13, 603–616. [Google Scholar] [CrossRef]
  24. Zhang, H.; Lin, F. Preliminary study on the characteristics and formation mechanism of small dam landslide in Fuquan, Guizhou. Sci. Technol. Eng. 2015, 15, 1617–1815. [Google Scholar]
  25. Du, J.; Wang, Y.; Peng, G.; Yin, K. Surge prediction of dashiban landslide in three gorges reservoir. Saf. Environ. Eng. 2007, 14, 92–95. [Google Scholar]
  26. Jin, D.; Wang, G. Tangyanguang landslide in Tuoxi Reservoir. In Typtical Landslides in China; Science Press: Beijing, China, 1986. (In Chinese) [Google Scholar]
  27. Sorensen, R.M. Basic Wave Mechanics: For Coastal and Ocean Engineers; John Wiley: New York, NY, USA, 1993. [Google Scholar]
  28. Kong, L. Fluid Mechanics; Higher Education Press: Beijing, China, 2004. [Google Scholar]
  29. Wang, S. Fluid Mechanics; China Electric Power Press: Beijing, China, 2007. [Google Scholar]
  30. Wang, Y.; Yin, K. Perturbation method of superposing initial surge height of landslide along reservoir shoreline. Chin. J. Rock Mech. Eng. 2004, 23, 717–720. [Google Scholar]
  31. Heller, V. Landslide Generated Impulse Waves: Prediction of Near Field Characteristics. Ph.D. Thesis, University of Nottingham, Nottingham, UK, 2007. [Google Scholar]
  32. Heller, V.; Hager, H.W.; Minor, H.-E. Landslide Generated Impulse Waves in Reservoirs: Basics and Computation; VAW Publikationen: Zürich, Switzerland, 2009. [Google Scholar]
  33. Heller, V.; Hager, W.H. Impulse Product Parameter in Landslide Generated Impulse Waves. J. Waterw. Port Coast. Ocean. Eng. 2010, 136, 145–155. [Google Scholar] [CrossRef]
  34. Panizzo, A.; De Girolamo, P.; Petaccia, A. Forecasting impulse waves generated by subaerial landslides. J. Geophys. Res. Ocean. 2005, 110, C12025. [Google Scholar] [CrossRef]
  35. Noda, E. Water waves generated by landslides. J. Waterw. Harb. Coast. Eng. 1970, 96, 835–853. [Google Scholar] [CrossRef]
  36. Mou, P.; Wang, P.; Han, L.; Meng, C. China’s safety production science and technology, review of physical model test research on landslide and surge disaster. J. Saf. Sci. Technol. 2020, 16, 43–49. [Google Scholar]
  37. Xu, W. Research on factors affecting landslide surge. J. Eng. Geol. 2012, 20, 491–507. [Google Scholar]
  38. Watts, P.; Grilli, S.T.; Kirby, J.T.; Fryer, G.J.; Tappin, D. Landslide tsunami case studies using a Boussinesq model and a fully nonlinear tsunami generation model. Nat. Hazards Earth Syst. Sci. 2003, 3, 391–402. [Google Scholar] [CrossRef]
  39. Di Risio, M.; Sammarco, P. Analytical modeling of landslide generated waves. J. Waterw. Port Coast. Ocean Eng. 2008, 134, 53–60. [Google Scholar] [CrossRef]
  40. Liu, P.L.-F.; Lynett, P.; Synolakis, C.E. Analytical solutions for forced long waves on a sloping beach. J. Fluid Mech. 2003, 478, 101–109. [Google Scholar] [CrossRef]
  41. Qin, Y.; Jiang, Q.; Guo, H. Discussion on two methods for predicting velocity of landslides. Rock Soil Mech. 2008, 29 (Suppl. S1), 373–378. [Google Scholar]
  42. Hu, J.; Wang, D.; Hu, B. Analysis on reservoir landslides and its surge. J. East China Jiaotong Univ. 2003, 20, 26–30. [Google Scholar]
  43. Haugen, K.B.; Løvholt, F.; Harbitz, C.B. Fundamental mechanisms for tsunami generation by submarine mass flows in idealized geometries. Mar. Pet. Geol. 2005, 22, 209–217. [Google Scholar] [CrossRef]
  44. Ursell, F.; Dean, R.G.; Yu, Y.S. Forced small-amplitude water waves: A comparison of theory and experiment. J. Fluid Mech. 1960, 7, 33–52. [Google Scholar] [CrossRef]
  45. Kennard, E.H. Generation of surface waves by a moving partition. Q. Appl. Math. 1949, 7, 303–312. [Google Scholar] [CrossRef]
  46. Das, M.M.; Wiegel, R.L. Waves generated by horizontal motion of a wall. J. Waterw. Port Coast. Ocean Eng. 1972, 98, 49–65. [Google Scholar] [CrossRef]
  47. Pan, J. Stability against Sliding of Buildings and Analysis of Landslide; China Water Power Press: Beijing, China, 1980. [Google Scholar]
  48. Li, M.; Wang, R.; Liu, X. On landslide surge and waveform spreading. J. Zhejiang Water Conserv. Hydropower Coll. 2003, 15, 5–7. [Google Scholar]
  49. Ha, Q.; Hu, W. The calculation of landslide surge in reservoir. Yellow River 1980, 2, 30–36. [Google Scholar]
  50. Law, L.; Brebner, A. On Water Waves Generated by Landslides. In Proceedings of the 3rd Australasian Conference on Hydraulics and Fluid Mechanics, Sydney, Australia, 25–29 November 1968; The Institution of Engineers: Sydney, Australia, 1968; pp. 155–159. [Google Scholar]
  51. Bruce, H. Water waves generated by distant landslides. J. Hydraul. Res. 1988, 26, 307–322. [Google Scholar]
  52. Wiegel, R.L. Laboratory studies of gravity waves generated by the movement of a submerged body. Trans. Am. Geo Phys. Union 1955, 36, 759–774. [Google Scholar] [CrossRef]
  53. Ataie-Ashtiani, B.; Nik-Khah, A. Impulsive waves caused by sub-aerial landslides. J. Environ. Fluid Mech 2008, 8, 263–280. [Google Scholar] [CrossRef]
  54. Ke, C.; Wang, Y.; Huo, Z.; Zhang, Y.; Liu, J. Study on the prediction model of initial surge amplitude of high and steep reservoir bank landslide. Saf. Environ. Eng. 2021, 28, 164–169. [Google Scholar]
  55. Kamphuis, J.W.; Bowering, R.J. Impulse waves generated by landslides. In Proceedings of the 12th International Conference on Coastal Engineering, Washington, DC, USA, 13–18 September 1970; American Society of Civil Engineers: Reston, VA, USA, 2016. [Google Scholar]
  56. Huang, J.; Zhang, T.; Li, J. Comparative analysis of empirical estimation methods for landslide surge on reservoir bank. Geotech. Mech. 2014, 35, 133–140. [Google Scholar]
  57. Ding, J.; Deng, H.; Wu, J.; Yong, M. Physical model test of landslide surge in a reservoir area of Lancang River. J. Yangtze River Acad. Sci. 2017, 34, 39–44. [Google Scholar]
  58. Li, R.; Jiang, C.; Deng, B.; Liu, Y. Surge height and propagation law of landslide in narrow channel near dam reservoir area. Traffic Sci. Eng. 2016, 32, 79–84. [Google Scholar] [CrossRef]
  59. Xiao, L.; Yin, K.; Wang, J.; Liu, Y. Reservoir bank landslide impact surge based on physical simulation test. J. Cent. South Univ. (Nat. Sci. Ed.) 2014, 45, 1618–1626. [Google Scholar]
  60. McFall, B.C.; Fritz, H.M. Physical modelling of tsunamis generated by three-dimensional deformable granular landslides on planar and conical island slopes. Proc. R. Soc. A Math. 2016, 472, 20160052. [Google Scholar] [CrossRef] [PubMed]
  61. Mohammed, F. Physical Modeling of Tsunamis Generated by Three-Dimensional Deformable Granular Landslides; Georgia Institute of Technology: Atlanta, GA, USA, 2010. [Google Scholar]
  62. Liu, Y. Research on Landslide-Induced Surge in Three Gorges Reservoir Area; China University of Geosciences: Wuhan, China, 2013. [Google Scholar]
  63. Xue, H.; Ma, Q.; Diao, M.; Jiang, L. Propagation characteristics of sub aerial landslide-generated impulse waves. Environ. Fluid Mech. 2018, 19, 203–230. [Google Scholar] [CrossRef]
  64. Li, Y.; Wang, P.; Hu, X. Calculation and research of landslide surge in mountain river-type reservoirs. J. Chongqing Jiaotong Univ. (Nat. Sci. Ed.) 2011, 30, 295–299. [Google Scholar]
  65. Yang, Q.; Wang, P.; Yu, T.; Chen, L. Analysis of surge climbing test of steep rock landslide in the Three Gorges reservoir area. China J. Geol. Hazards Prev. 2014, 25, 43–48+55. [Google Scholar]
  66. Hall, J.V.; Watts, G.M. Laboratory Investigation of the Vertical Rise of Solitary Wave on Impermeable Slopes; U. S. Army Corps of Engineers Beach Erosion Board: Washington, DC, USA, 1953. [Google Scholar]
  67. Synolakis, C.E. The Runup of Solitary Waves. J. Fluid Mech. 1987, 185, 523–545. [Google Scholar] [CrossRef]
  68. Müller, D.R. Auflaufen Undüberschwappen von Impulswellen an Talsperren (Run-Up and Run-Over of Impulse Waves at Dams); VAW-Mitteilung 137; ETH Zurich: Zurich, Switzerland, 1995. [Google Scholar]
  69. Biscarini, C. Computational fluid dynamics modelling of landslide generated water waves. Landslides 2010, 7, 117–124. [Google Scholar] [CrossRef]
  70. Serrano-Pacheco, A.; Murillo, J.; García-Navarro, P. A finite volume method for the simulation of the waves generated by landslides. J. Hydrol. 2009, 373, 273–289. [Google Scholar] [CrossRef]
  71. Xu, W. CEL algorithm study of reservoir surge Induced by landslide. J. Eng. Geol. 2012, 20, 350–354. [Google Scholar]
  72. Li, W.; Yue, G.; Wang, H. Dynamic mechanism of high speed landslide surge. J. Nat. Disasters 2013, 22, 127–133. [Google Scholar]
  73. Song, X.; Xing, A.; Chen, L. Numerical simulation of two-dimensional water waves due to landslide based on FLUENT. Hydrogeol. Eng. Geol. 2009, 36, 90–94. (In Chinese) [Google Scholar]
  74. Yuk, D.; Yim, S.C.; Liu, P.F. Numerical modeling of submarine mass-movement generated waves using RANS model. Comput. Geosci. 2006, 32, 927–935. [Google Scholar] [CrossRef]
  75. Abadie, S.; Morichon, D.; Grilli, S.; Glockner, S. Numerical simulation of waves generated by landslides using a multiple-fluid Navier-Stokes model. Coast. Eng. 2010, 57, 779–794. [Google Scholar] [CrossRef]
  76. Liu, P.F.; Wu, T.R.; Raichlen, F.; Synolakis, C.E.; Borrero, J.C. Run-up and rundown generated by three-dimensional sliding masses. J. Fluid Mech. 2005, 536, 107–144. [Google Scholar] [CrossRef]
  77. Liu, G.; Liu, M. Smoothed Particle Hydrodynamics—A Meshfree Particle Method; World Scientific Publishing Company: Singapore, 2003; pp. 26–32. [Google Scholar]
  78. Ren, K.; Han, J. Experimental research on primary wave height generated by loose earth landslide. Yangtze River 2011, 42, 69–72. [Google Scholar]
  79. Xu, B.; Jiang, C.; Deng, B. Three dimensional numerical simulations of water waves generated by landslides and its propagation process. J. Transp. Sci. Eng. 2011, 27, 39–45. [Google Scholar]
  80. Bosa, S.; Petti, M. A numerical model of the wave that overtopped the Vajont Dam in 1963. Water Resour. Manag. 2013, 27, 1763–1779. [Google Scholar] [CrossRef]
  81. Ataie-Ashtiani, B.; Yavari-Ramshe, S. Numerical simulation of wave generated by landslide incidents in dam reservoirs. Landslides 2011, 8, 417–432. [Google Scholar] [CrossRef]
  82. Wang, S.; Huang, B.; Liu, G.; Chen, X. Numerical simulation of tsunami due to slope failure at Gongjiafang on Three Gorges Reservoir. Rock Soil Mech. 2015, 36, 212–218. (In Chinese) [Google Scholar]
  83. Heidarzadeh, M.; Pirooz, M.D.; Zaker, N.H.; Yalciner, A.C.; Mokhtari, M.; Esmaeily, A. Historical tsunami in the Makran Subduction Zone off the southern coasts of Iran and Pakistan and results of numerical modeling. Ocean Eng. 2008, 35, 774–786. [Google Scholar] [CrossRef]
  84. Jiang, Z.B.; Jin, F.; Sheng, J. Numerical simulations of water waves due to landslide. J. Yangtze River Sci. Res. Inst. 2005, 22, 1–3. [Google Scholar]
  85. Shen, Y.M.; Ng, C.O.; Zheng, Y.H. Simulation of wave propagation over a submerged bar using the VOF method with a two-equation k-εturbulence modeling. Ocean Eng. 2004, 31, 87–95. [Google Scholar] [CrossRef]
  86. Teng, G.; Zhang, Y.; Chen, J. Comparison of the Application of the Empirical Estimation Method for Reservoir Bank Landslide Surge. Henan Sci. Technol. 2021, 40, 69–72. [Google Scholar]
  87. Deng, C.; Dang, F.; Chen, X. Study on the surge wave propagation in the reservoir area and its interaction mechanism with the dam. J. Hydraul. Eng. 2019, 50, 815–823. [Google Scholar]
Figure 1. Relative position of landslide and water (a) landslide surges induced by the collapse masses above water surface, which involves the movement of solid, water and air. (b) landslide surges induced by collapse masses that are partially under water surface. (c) landslide surges induced by underwater collapse masses.
Figure 1. Relative position of landslide and water (a) landslide surges induced by the collapse masses above water surface, which involves the movement of solid, water and air. (b) landslide surges induced by collapse masses that are partially under water surface. (c) landslide surges induced by underwater collapse masses.
Water 15 01879 g001
Figure 2. Schematic diagram of landslide.
Figure 2. Schematic diagram of landslide.
Water 15 01879 g002
Figure 3. Schematic diagram for three stages of landslide-generated impulse waves.
Figure 3. Schematic diagram for three stages of landslide-generated impulse waves.
Water 15 01879 g003
Table 1. Major cases of landslide surges worldwide.
Table 1. Major cases of landslide surges worldwide.
NumberTimeLocationLandslide Volume
(104 m3)
Failure ModeNumber of Deaths
11792Senyun, Japan [12]535The formation of 10 m high surge.15,000
21933Diexi, Minjiang River, China [13]--8800
31956Lanfjord, Norway [14]12The formation of 140 m high surge.32
41958Lituya landslide in Alaska, USA [15]30The formation of 30 m high surge.2
51958Italy Pontesei arch dam reservoir area [15]-A 20 m high surge over the dam.1
61961Zishui Zhexi Reservoir, Hunan Province, China [16]1.65Damage caused by 21 m high surge over dam.40
71963Italy vajont reservoir [17]240A 175 m surge overtop3000
81971Peru chungar lake shore [15]0.1The formation of 30 m high surge.300–500
91980Mount St. Helens Spirit Lake [18] 2500The landslide whipped up a more than 200 m high surge-
101982Jibazi landslide in Yunyang County, Three Gorges Reservoir Area, China [19]-The river shoreline advanced inward more than 50 m, the riverbed silt height of 30 m.-
111985Xintan Landslide, Zigui County, Three Gorges Reservoir Area, China [20,21]30Initial waves of up to 49 m destroyed nearby ships.10
122002Stromboli, Italy [15]-The formation of 10.9 m high surge.2
132003Qianjiangping Village in Three Gorges Reservoir Area of China [20]24Up to 30 m surges destroyed 120 acres of farmland nearby.14
142007Shuibuya Dam, China [22]3The formation of 50 m high surge.1
152008GongJiaFang, Wushan, Chongqing, China [23]0.38Waves about 15 m high formed on the opposite bank.-
162009The landslide in Xiaowan1The landslide whipped up a 30 m high surge24
172014Xiaoba landslide in Guizhou province [24]0.33The landslide whipped up a 20 m high surge12
182015Chongqing Wushan Daning River Jiangdong Temple North Shore [25,26]-A 6 m high surge toppled nearby boats.1
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Xin, C.; He, Z.; Feng, W. A Review of Research Methods and Evolution Mechanisms of Landslide-Induced Tsunamis. Water 2023, 15, 1879. https://doi.org/10.3390/w15101879

AMA Style

Xin C, He Z, Feng W. A Review of Research Methods and Evolution Mechanisms of Landslide-Induced Tsunamis. Water. 2023; 15(10):1879. https://doi.org/10.3390/w15101879

Chicago/Turabian Style

Xin, Chunlei, Zhiqian He, and Wenkai Feng. 2023. "A Review of Research Methods and Evolution Mechanisms of Landslide-Induced Tsunamis" Water 15, no. 10: 1879. https://doi.org/10.3390/w15101879

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop