GIS

Aug 19, 2023

Scientific Reports volume 13, Article number: 7684 (2023) Cite this article

374 Accesses

Metrics details

Combined with the spatial data processing capability of the Geographic Information Systems (GIS), the Pan Jiazheng method is extended from two-dimensional (2D) to three-dimensional (3D), and a 3D landslide surge height calculation method is proposed based on grid column units. First, the data related to the landslide are rasterized to form grid columns, and a force analysis model of 3D landslides is established. Combining the vertical strip method with Newton's laws of motion, dynamic equilibrium equations are established to solve the surge height. Moreover, a 3D landslide surge height calculation expansion module is developed in the GIS environment, and the results are compared with those of the 2D Pan Jiazheng method. Comparisons showed that the maximum surge height obtained by the proposed method is 24.6% larger than that based on the Pan Jiazheng method. Compared with the traditional 2D method, the 3D method proposed in this paper better represent the actual spatial state of the landslide and is more suitable for risk assessment.

A landslide disaster refers to a disaster caused by the overall downhill slide of rock mass or soil mass under the action of gravity and is one of the main geological disasters in the world. Quick sliding speeds and long sliding distances of medium and large landslides result in massive losses of life and property every year. Among these, the reservoir bank landslide produces a bigger surge during its sliding into the water, and this causes great harm to the passing ships and the surrounding buildings, receiving much concern around the globe. There are many surging events caused by landslides in the world, such as the 100 m high surges triggered by the landslide in the Vajont Reservoir, Italy, in 1963, killing at least 2,000 people1; the red rock surges in Wushan, Chongqing, China in 20152; and the Dayantang surges in Hubei, China in 20143; all of which brought about large casualties and economic losses. Therefore, it is of vital significance to evaluate the surge wave hazard caused by the reservoir bank landslide.

Calculating the surge height is one of the key indexes to evaluate the surge hazard. The methods of calculating the surge height can mainly be divided into theoretical analysis method4,5,6,7,8, numerical simulation method10,11,12,13,14, and physical modelling method15,16,17. Among them, the Pan Jiazheng method in the theoretical analysis method is widely used in engineering applications because of its simple modelling processes, which has few requirements for engineers and high precision.

The Pan Jiazheng method originated from Noda. Noda4 proposed an approximate method to find the amplitude of the largest surge in the nonlinear region by utilizing the solutions obtained from linear theory. Since then, many scholars have conducted more in-depth research. On this basis, Academician Pan Jiazheng5 divided the landslide body into many two-dimensional (2D) vertical strips and calculated the surge height by considering the horizontal and vertical movement of the landslide. This method is called the Pan Jiazheng method. The method has been applied and improved over the years. For example, Dai et al.6 used Pan Jiazheng method to calculate the sliding speed of Xiaduling landslide in the Three Gorges Reservoir area. Huang et al.7 improved the Pan Jiazheng method by considering the resistance of water and the change in the friction coefficient. Miao et al.8 proposed a sliding block model based on the 2D vertical strip method to predict the maximum surge height.

Although the Pan Jiazheng method has made some improvements, it is still in the 2D stage, and the calculation results will be different if the 2D section is selected. However, the actual state of the landslide is three-dimensional (3D), and the 2D analysis method cannot simulate the real landslide state in a reasonable way. Hu18 proposed that the value obtained by 2D state analysis is about 70% of the 3D state value. The Pan Jiazheng method based on 2D model has been unable to obtain high-precision result values, so the study of the Pan Jiazheng method based on 3D state analysis is of great significance to improve its calculation accuracy.

Geographic Information Systems (GIS) is widely used in geotechnical engineering19,20,21. GIS features strong spatial analysis capabilities. It can carry out unified management and storage of spatial data, with functions such as spatial positioning, spatial database management, digital elevation model establishment, etc. And it can provide users with multi-factor spatial analysis, prediction and forecasting, simulation optimization, and other analysis methods of spatial data. Its biggest feature is that it can transform vector data into raster dataset based on grid column unit model. Because of its good 3D spatial data processing capability, many scholars have added geotechnical professional models to their geographic information systems22,23. Xie's team24,25,26 built a 3D limit equilibrium method in the GIS, and developed the slope stability analysis module. This paper, based on the previous study, will establish a set of 3D landslide surge height calculation methods grounded in the GIS.

Coinciding with the spatial data processing capability of the GIS, this paper extends Pan Jiazheng's method from 2 to 3D and proposes a 3D landslide surge height calculation method based on the GIS. Firstly, considering the 3D spatial relation of grid column units, the calculating expression of the required parameters is given. Through the force analysis of the grid column unit, combined with Newton's law of motion and considering the action of water, the dynamic equation based on the grid column unit is established, and the sliding speed of the sliding body during the sliding process is solved, and then the surge height is calculated. At the same time, an expansion module for calculating surge height is developed based on Component Object Model (COM) technology in the ArcGIS environment, and the calculation is carried out through the case of the Kaiding landslide of Houziyan Hydropower Station, and compared with the results of 2D Pan Jiazheng method, which verified the applicability of the module.

After this introduction, the second section introduces the GIS-based 3D landslide surge height calculation method, gives the specific surge height calculation formula, and introduces the secondary development of the extension module. The third section verifies the correctness and applicability of the extension module through case calculation. Finally, the fourth section comes to discussion and summary.

For a slope, the representation of data is mainly in the form of vectors. These data include but are not limited to slip surface, strata, groundwater, fault, slip, and other types of data. These vector data layers can be converted to raster data layers using the spatial analysis capabilities of GIS to form a grid data set. The grid data structure consists of rectangular units. Each rectangular unit has a corresponding row and column number and is assigned an attribute value that represents the grid unit27,28,29,30,31. Therefore, the slope can be divided into square columns based on the grid units to form a grid column unit model, as shown in Fig. 1.

3D space model of landslide ((a) 3D schematic diagram of landslide, (b) 3D schematic diagram of one column).

Figure 2 shows the spatial relationship between parameters. θ is the dip of the grid column at the slip surface; α is the dip direction of the grid column at the slip surface; β is the sliding direction of the landslide; θr is the apparent dip of the main inclination direction of the landslide; αx is the apparent dip of the X-axis; and αy is the apparent dip of the Y-axis. Of these six parameters, α, θ, and β are known and is calculated in a paper published by myself25. Other three parameters αx, αy and θr can be calculated according to the spatial relationship in Fig. 2b, which are calculated as follows.

The spatial relationship between parameters. (a) Is the spatial relationship diagram of each parameter in the grid column; (b) is the calculation relationship diagram between parameters.

The bottom area of one grid column is calculated by

where cellsize represents the length of the side of each grid column.

The weight W of the grid column is expressed as

where m is the number of strata, hm is the height of each stratum, and rm is the unit weight of each stratum. For the grid column units above the water, rm is calculated from the natural unit weight. For grid column units under water, rm is calculated from the buoyant unit weight.

The pore water pressure is obtained as follows32.

where R is the distance from the centre bottom of the grid column to the water surface.

When the sliding body enters the water, the resistance of the water is calculated as follows32.

where G is the resultant force of the resistance of the water to the sliding body (mN); cw is the viscous resistance coefficient, which is 0.15 to 0.18; ρf is the buoyant density (g/m3), taking the average of all stratum; v is the velocity of the landslide (m/s); and S is the surface area of the grid column in the water (m2).

As shown in Fig. 3, it is the force analysis of one grid column. In the 3D sliding body, one random grid column ABCDA1B1C1D1 is selected, and the force analysis is explained as follows26.

The weight of one grid column is W; the direction is the Z-axis; and the weight acts at the centroid of the grid column.

The resultant horizontal seismic force is kW, where k is the “seismic coefficient”; the direction of kW is the sliding direction of the landslide; and the resultant horizontal force acts at the centroid of the grid column.

The external loads on the ground surface are represented by P; the direction of P is the Z-axis, and these external loads act at the centre of the top of the grid column. The external loads represent loads caused by objects on the surface of the landslide, such as buildings, trees, and so on.

The normal and shear stresses on the slip surface are represented by σ and τ, respectively. The normal stress is perpendicular to the slip surface, and the shear stress is in the sliding direction of the landslide. The normal and shear stresses act at the centre of the bottom of the grid column.

The pore water pressure on the slip surface is u. The direction of u is directed as σ.

The horizontal tangential forces on the vertical face at y = 0 and vertical face at y = △y (△y represents the size of the grid column along Y-axes) are T and T + △T, respectively; the vertical tangential forces on the vertical face at y = 0 and vertical face at y = △y are R and R + △R, respectively; the normal forces on the vertical face at y = 0 and vertical face at y = △y are F and F + △F, respectively; the horizontal tangential forces on the vertical face at x = 0 and vertical face at x = △x are E and E + △E, respectively; the vertical tangential forces on the vertical face at x = 0 and vertical face at x = △x are V and V + △V, respectively; and the normal forces on the vertical face at x = 0 and vertical face at x = △x are H and H + △H, respectively. For convenience, the resultant force between columns in the sliding direction of the landslide is defined as ΔD.

Force analysis of one grid column.

We assume that all the grid column units move continuously, do not separate in the macroscopic dimension and remain vertical after sliding, as also assumed by Pan Jiazheng5. The force analysis of one grid column and the spatial relationships between parameters at the slip surface are shown in Figs. 2 and 3, respectively.

We arbitrarily selected a grid column unit (the grid column unit in row i and column j). According to Newton’s laws of motion, dynamic equilibrium equations are established in the sliding direction of the landslide and the vertical direction. The force analyses in the sliding direction of the landslide and vertical direction are shown in Fig. 4.

where ax and a are the horizontal acceleration and vertical acceleration of the grid column, respectively; φ is the effective friction angle of the grid column at the slip surface; g is gravitational acceleration; c is the effective cohesion of the grid column at the slip surface.

Force analysis in the vertical direction and sliding direction of the landslide.

According to this assumption, the horizontal acceleration ax of each grid column unit is the same, and the vertical acceleration a of each grid column unit varies. Pan Jiazheng5 suggested that there is a certain proportional relationship between ax and a, that is, aax = tanδ. δ is the horizontal inclination angle of the line connecting the centre bottom of the grid column to the centre bottom of the next grid column in the sliding direction of the landslide. The effect of vertical tangential forces is ignored, namely, ΔV-ΔR = 0; therefore, Eq. (8) can be transformed as follows.

The simultaneous Eqs. (9) and (10) can be obtained as follows.

For the entire sliding body, the forces between the grid columns are internal forces, that is, the resultant force is 0, yielding Eq. (12).

By summing all the grid column units, the horizontal acceleration ax can be determined by Eq. (7).

Substituting Eqs. (9) and (11) into Eq. (13) yields the following equation.

where

The steps in calculating the landslide sliding velocity are as follows.

Using the spatial analysis capability of GIS, the landslide body is rasterized, and the size of the grid column unit (i, j) can be set to an arbitrary square. A partitioning line is drawn from the bottom to the top of the landslide every ΔL in the sliding direction of the landslide, and the resulting regions are numbered zone 1, zone 2, zone 3, …, zone (n-1), zone n. Each partition includes numbers of grid column units, and the length of zone n is less than or equal to ΔL, as shown in Fig. 5. For a grid column unit that is not completely contained within a partition, if the area within the partition is greater than half of the total area, the unit is divided into that partition; otherwise, the unit is divided into the next partition.

For each grid column unit, the parameters required in Eq. (14) are calculated.

t0 is the starting point of when the landslide body begins to slide, and t0 = 0. When the landslide body moves distance ΔL sequentially in the sliding direction of the landslide, the corresponding time is recorded as t1, t2, t3…tn, and the corresponding velocity is expressed as vx1, vx2, vx3…vxn.

The horizontal acceleration at t0 can be calculated by Eq. (14) and is denoted as ax0, and the velocity at time t0 is zero. After sliding distance ΔL is reached, the following equations can be obtained.

At t = t1, the landslide body has horizontally moved by a distance ΔL in the sliding direction of the landslide, zone 1 has slipped from the sliding surface. The horizontal acceleration ax1 at t1 is still calculated by Eq. (14). Unlike t0, the weight for zone (n-1) changes to the weight for zone n, and the weight for zone (n-2) becomes the weight for zone (n-1), and so on (At this time, there is no grid column for zone n). After ax1 is calculated, the following can be established.

The calculation is continued in turn. When the obtained horizontal acceleration is negative, the maximum velocity can be obtained. Finally, ax and vx in the calculation process can be plotted as respective curves versus the sliding time.

Rasterization and partitioning of landslides.

The China Institute of Water Resources and Hydropower Research proposed an empirical formula for surge height calculation33. In the formula, the main factors that affect the surge height are the sliding velocity and volume of the landslide. The formula for calculating the maximum surge height is as follows.

where ξmax is the maximum surge height (m); d is the comprehensive influence coefficient, with an average value of 0.12; vm is the maximum sliding velocity (m/s); V is the volume of the landslide body in the water (m3); and g is gravitational acceleration, which equals 9.8 m/s2.

The formula for calculating the surge height at different distances from the landslide body is as follows.

where ξ is the surge height at a point from the landslide body L metres (m); n is the calculation coefficient, which is 1.4; and d1 is the influence coefficient related to distance L, which is determined by the following formula.

Combined with the surge height calculation method, an expansion module was developed based on component object model (COM) technology in the ArcGIS environment25,26.

First, rasterize the data related to slope in ArcGIS software to form a raster data set (including the information of elevation, strata, groundwater, fault, slip surface, etc.), and calculate the required parameters by using the extension module. Then, divide the sliding body into units, and calculate the sliding speed at each time point step by step with the extension module. And the calculation stops when the acceleration is less than 0. Finally, obtain the maximum sliding speed and calculate the maximum surge height. Fig. 6 illustrates the algorithmic steps.

The algorithmic steps.

The Kaiding landslide is approximately 14.5 km away from the dam of the Houziyan hydropower station in Sichuan, China. The length of the landslide along the river is approximately 490 m, the top elevation is 2080 m, the bottom elevation is 1754 m, and the volume is approximately 4.5 million m3. Plan and section views are shown in Figs. 7 and 8, respectively. The scale in Fig. 7 is 1:50,000, and the river direction is north to south.

The plan view of the Kaiding landslide.

The section view of the Kaiding landslide.

The unit size of a grid column is 5 m × 5 m, and ΔL = 10 m. The internal friction angle φ at the slip surface is 22.8°, the natural unit weight is 18.84 kN/m3, the buoyant unit weight is 19.43 kN/m3, the buoyant density is 2.11 × 106 g/m3, the viscous resistance coefficient is 0.18, and the elevation of the reservoir water level is 1810.3 m. When the landslide body slides, the effective cohesion c at the slip surface will decrease to 0, that is, c = 05. Using this method and Pan Jiazheng's 2D method, the acceleration and velocity curves with the sliding time can be obtained, as shown in Fig. 9a,b, respectively. The calculation results are shown in Table 1.

(a) Horizontal acceleration curve with the sliding time. (b) Sliding velocity curve with the sliding time.

The calculation results indicate that the maximum velocity obtained by the proposed method is 11.21 m/s, the starting acceleration is 1.25 m/s2, and the sliding time required to reach the maximum velocity is 8.67 s. In comparison, the maximum velocity obtained by the Pan Jiazheng method is 9.51 m/s, the starting acceleration is 0.84 m/s2, and the sliding time required to reach the maximum velocity is 9.73 s.

Comparing the results of the proposed method with those of the Pan Jiazheng method, the maximum velocity of the proposed method is 15.2% higher than that calculated by the Pan Jiazheng method, the starting acceleration is 32.8% higher, and the sliding time required to reach the maximum velocity is 1.06 s short.

According to the most dangerous working conditions, it is assumed that the landslide body all slips into the water. The volume V of the landslide body under water is 340 × 104 m3. According to Eqs. (25) and (26), the maximum surge height obtained by the proposed method is 9.66 m, and the surge height at the dam site is 0.56 m. The maximum surge height obtained by the Pan Jiazheng method is 7.28 m and the surge height at the dam site is 0.44 m. The calculation results are shown in Table 2.

The landslide is approximately 14.5 km from the dam, the crest elevation is 1847.02 m, and the elevation of the reservoir water level is maintained at 1810.3 m. When the surge height at the dam site is 0.56 m, water will not flow over the dam crest and the safe operation of the dam will not be affected.

The maximum surge height obtained by the proposed method is 24.6% larger than that based on the Pan Jiazheng method, and the surge height at the dam site obtained by the proposed method is 21.4% larger than that based on the Pan Jiazheng method.

The calculation results show that the difference between the results of the 2D method and the 3D method is 24.6%. Hu13 proposed that the value obtained by 2D state analysis is about 70% of the 3D state value. The result is consistent with Hu's prediction. Compared that of the 2D method, the computational model of the 3D method better represents the actual spatial state of the landslide. Therefore, the method in this paper is more applicable than Pan Jiazheng method in actual risk evaluation.

This paper proposes a 3D landslide surge height calculation method, and it divides the landslide into several grid column units. The method assumes that all the grid column units move continuously and do not separate in the macroscopic dimension, remaining vertical after sliding. That is, assume that the column units are rigid materials. But in the actual sliding process, the column units are unable to keep vertical regularly, especially in the soil landslide, while for the rock landslide, its column units can keep better integrity during the sliding process, so the method in this paper will have relatively good applicability for the rock landslide.

Combined with the spatial data processing capability of the GIS, the Pan Jiazheng method is extended from 2 to 3D, and a 3D landslide surge height calculation method is proposed for the first time. Combined with Newton's law of motion, the dynamic balance equation for calculating the sliding speed of a 3D sliding body is derived, and then the surge height is calculated.

This is the first time to combine the surge height calculation model with GIS. At the same time, an extension module is developed based on ArcGIS software, and the feasibility of the module is verified by a case study. The module has the advantages of unified data format and simple preparation process.

Because the Pan Jiazheng method is a calculation method focused on 2D sections, the calculation results will be different if different sections are selected. After the 3D landslide is carried out in rasterization, the 3D calculation model based on the grid column unit is established in this paper, which overcomes the above defects and makes the calculation model close to the actual situation, therefore the method in this paper is more applicable than Pan Jiazheng method in actual risk evaluation.

As the application of GIS in geotechnical engineering becomes increasingly extensive, the calculation method of surge height established in GIS in this paper will provide a theoretical basis for scholars to add surge height calculation modules in their respective geographic information Systems.

The datasets used and/or analysed during the current study available from the corresponding author on reasonable request.

Zhao, T., Utili, S. & Crosta, G. B. Rockslide and impulse wave modelling in the Vajont reservoir by DEM-CFD analyses. Rock Mech. Rock Eng. 49, 2437–2456 (2016).

Article ADS Google Scholar

Han, B. et al. The monitoring-based analysis on deformation-controlling factors and slope stability of reservoir landslide: Hongyanzi Landslide in the Southwest of China. Geofluids 35, 1–14 (2018).

Google Scholar

Wang, J. X., Xiao, L. L., Ward, S. N. & Du, J. Tsunami Squares modeling of the 2007 Dayantang landslide generated waves considering the effects in slide/water interactions. Eng. Geol. 284, 342–354 (2021).

Article Google Scholar

Risio, M. D., Bellotti, G., Panizzo, A. & Girolamo, P. D. Three-dimensional experiments on landslide generated waves at a sloping coast. Coast. Eng. 56, 659–671 (2009).

Article Google Scholar

Xiao, L., Yin, K., Wang, J. & Liu, Y. Impulse wave generated by reservoir landslides based on physical model experiments. J. Central South Univ. (Sci. Technol.) 45, 1618–1626 (2014).

Google Scholar

Dai, Y. X., Yin, K. L. & Wang, Y. Discussion on method of landslide velocity calculation and surge prediction. Rock Soil Mech. 29, 407–411 (2008).

Google Scholar

Huang, J. Z., Zhong, Z. H. & Zhang, M. F. Sliding velocity analysis of reservoir bank landslides based on the improved vertical slice method. J. Mt. Sci. 30, 555–560 (2012).

Google Scholar

Miao, T. D., Liu, Z. Y., Niu, Y. H. & Ma, C. W. A sliding block model for the runout prediction of high-speed landslides. Can. Geotech. J. 38, 217–226 (2001).

Article Google Scholar

Wang, M. L., Tian, Y. & Li, X. L. Experimental study on pressure distribution of spur dike under the combined action of landslide surge and water flow. AIP Adv. 12, 1221–1242 (2022).

Google Scholar

Sassa, K., Dang, K. H., Yanagisawa, H. & He, B. A new landslide-induced tsunami simulation model and its application to the 1792 Unzen-Mayuyama landslide-and-tsunami disaster. Landslides 13, 1–15 (2016).

Article Google Scholar

Yavari-Ramshe, S. & Ataie-Ashtiani, B. Numerical modeling of subaerial and submarine landslide-generated tsunami waves-recent advances and future challenges. Landslides 13, 1325–1368 (2016).

Article Google Scholar

Majd, M. S. & Sanders, B. F. The LHLLC scheme for two-layer and two-phase transcritical flows over a mobile bed with avalanching, wetting and drying. Adv. Water Resour. 67, 16–31 (2014).

Article ADS Google Scholar

Yin, Y. P., Huang, B. L., Liu, G. N. & Wang, S. C. Potential risk analysis on a Jianchuandong dangerous rockmass-generated impulse wave in the Three Gorges Reservoir, China. Environ. Earth Sci. 74, 2595–2607 (2015).

Article ADS Google Scholar

Brune, S., Babeyko, A. Y., Gaedicke, C. & Ladage, S. Hazard assessment of underwater landslide-generated tsunamis: A case study in the Padang region, Indonesia. Nat. Hazards 53, 205–218 (2010).

Article Google Scholar

Cui, P. & Zhu, X. Surge generation in reservoirs by landslides triggered by the wenchuan earthquake. J. Earthq. Tsunami 05, 461–474 (2011).

Article Google Scholar

Huang, B. L. et al. A physical similarity model of an impulsive wave generated by Gongjiafang landslide in Three Gorges Reservoir, China. Landslides 11, 513–525 (2014).

Article Google Scholar

Huang, B. L., Yin, Y. P. & Du, C. L. Risk management study on impulse waves generated by Hongyanzi landslide in Three Gorges Reservoir of China on June 24, 2015. Landslides 13, 603–616 (2016).

Article Google Scholar

Hu, G. T. Landslide dynamics (Geological Press) 105–108 (Beijing, 1995).

Xie, M., Esaki, T. & Cai, M. GIS-based implementation of three-dimensional limit equilibrium approach of slope stability. J. Geotech. Geoenviron. Eng. 132, 656–666 (2006).

Article Google Scholar

Senouci, R., Taibi, N. E., Teodoro, A. C., Duarte, L. & Meddah, R. Y. GIS-based expert knowledge for landslide susceptibility mapping (LSM): Case of Mostaganem Coast District, West of Algeria. Sustainability 13, 477–486 (2021).

Article Google Scholar

Jia, N., Mitani, Y., Xie, M. W., Tong, J. & Yang, Z. H. GIS deterministic model-based 3D large-scale artificial slope stability analysis along a highway using a new slope unit division method. Nat. Hazards 76, 873–890 (2015).

Article Google Scholar

Jia, N., Yang, Z., Xie, M. W., Mitani, Y. & Tong, J. X. GIS-based three-dimensional slope stability analysis considering rainfall infiltration. Bull Eng Geol Environ 74, 919–931 (2015).

Article Google Scholar

Mergili, M., Marchesini, I., Rossi, M., Guzzetti, F. & Fellin, W. Spatially distributed three-dimensional slope stability modelling in a raster GIS. Geomorphology 206, 178–195 (2014).

Article ADS Google Scholar

Jia, N., Yang, Z. H., Xie, M. W., Mitani, Y. & Tong, J. X. GIS-based three-dimensional slope stability analysis considering rainfall infiltration. Bull. Eng. Geol. Env. 74, 919–931 (2015).

Article Google Scholar

Yu, G., Bu, L., Wang, C. F. & Farooq, A. Composition analysis and distributed assumption GIS model of normal stress on the slope sliding surface. Front. Earth Sci. 10, 123–135 (2022).

Article ADS CAS Google Scholar

Yu, G., Xie, M. W., Liang, J., Asim, F. & Williams, E. J. A GIS-based 3D slope stability analysis method based on the assumed normal stress on the slip surface. Sci. Rep. 10, 291–298 (2020).

Google Scholar

Xie, M., Esaki, T. & Zhou, G. GIS-based probabilistic mapping of landslide hazard using a three-dimensional deterministic model. Nat. Hazards 33, 265–282 (2004).

Article Google Scholar

Bhagyaraj, U. & Kolathayar, S. Seismic hazard assessment and landslide vulnerability mapping for Ladakh, and Jammu & Kashmir using GIS technique. J. Geol. Soc. India 99, 377–382 (2023).

Article Google Scholar

Agrawal, N. & Dixit, J. Assessment of landslide susceptibility for Meghalaya (India) using bivariate (frequency ratio and Shannon entropy) and multi-criteria decision analysis (AHP and fuzzy-AHP) models. All Earth 34, 179–201 (2022).

Article Google Scholar

Dai, Y., Bai, X., Nie, G. & Huang, F. G. A rapid assessment method for earthquake-induced landslide casualties based on GIS and logistic regression model. Geomat. Nat. Haz. Risk 13, 222–248 (2022).

Article Google Scholar

Coco, L., Macrini, D., Piacentini, T. & Buccolini, M. Landslide susceptibility mapping by comparing GIS-Based bivariate methods: A focus on the geomorphological implication of the statistical results. Remote Sensing 13, 260–274 (2021).

Article Google Scholar

Habib, P. Production of gaseous pore pressure during rock slides. Rock Mech. 7, 193–197 (1975).

Article Google Scholar

Peng, M. J. et al. Stability analysis of landslide dams under surge action based on large-scale flume experiments. Eng. Geol. 259, 276–288 (2019).

Article Google Scholar

Download references

I would like to express my sincere gratitude to Prof. Xiaowen Zhou, Lei Bu, Chengfeng Wang, and Asim Farooq for their motivation and for providing me access to their immense knowledge during this research work. This work was supported by the National Natural Science Foundation of China [Grant Numbers 51508203].

School of Civil Engineering and Transportation, South China University of Technology, Guangzhou, 510640, China

Guo Yu & Chengfeng Wang

State Key Laboratory of Subtropical Building Science, South China University of Technology, Guangzhou, 510640, China

Xiaowen Zhou

China Coal Technology & Engineering Group, Nanjing Design & Research Institute Co., Ltd, Nanjing, 210031, China

Lei Bu

Center of Excellence in Transportation Engineering, Pak Austria Facshhoule Institute of Applied Sciences and Technology, Khanpur Road, Haripur, Pakistan

Asim Farooq

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

G.Y. was responsible for the implementation of the algorithm and programming work. X.Z. was responsible for the ideas of the article. L.B. and C.W. was responsible for helping with the programming work. A.F. was responsible for helping with the programming work.

Correspondence to Xiaowen Zhou.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Yu, G., Zhou, X., Bu, L. et al. GIS-based calculation method of surge height generated by three-dimensional landslide. Sci Rep 13, 7684 (2023). https://doi.org/10.1038/s41598-023-34798-1

Download citation

Received: 02 December 2022

Accepted: 08 May 2023

Published: 11 May 2023

DOI: https://doi.org/10.1038/s41598-023-34798-1

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.