Next Article in Journal
A Comparative Approach of Fuzzy Object Based Image Analysis and Machine Learning Techniques Which Are Applied to Crop Residue Cover Mapping by Using Sentinel-2 Satellite and UAV Imagery
Next Article in Special Issue
Land Subsidence Susceptibility Mapping Using Persistent Scatterer SAR Interferometry Technique and Optimized Hybrid Machine Learning Algorithms
Previous Article in Journal
Sea Ice Thickness Retrieval Based on GOCI Remote Sensing Data: A Case Study
Previous Article in Special Issue
Rapid Assessment of Hillslope Erosion Risk after the 2019–2020 Wildfires and Storm Events in Sydney Drinking Water Catchment
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

A New Approach for Identification of Potential Rockfall Source Areas Controlled by Rock Mass Strength at a Regional Scale

1
Key Laboratory of Shale Gas and Geoengineering, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 100029, China
2
Innovation Academy for Earth Sciences, Chinese Academy of Sciences, Beijing 100029, China
3
University of Chinese Academy of Sciences, Beijing 100049, China
*
Author to whom correspondence should be addressed.
Submission received: 7 February 2021 / Revised: 24 February 2021 / Accepted: 26 February 2021 / Published: 3 March 2021
(This article belongs to the Special Issue Remote Sensing for Natural Hazards Assessment and Control)

Abstract

:
The identification of rockfall source areas is a fundamental work for rockfall disaster prevention and mitigation. Based on the Culmann model, a pair of important indicators to estimate the state of slope stability is the relief and slope angles. Considering the limit of field survey and the increasing requirements for identification over a large area, a new approach using the relief–slope angle relationship to identify rockfall source areas controlled by rock mass strength at a regional scale is proposed in this paper. Using data from helicopter-based remote sensing imagery, a digital elevation model of 10 m resolution, and field work, historical rockfalls in the Wolong study area of Tibet where frequent rockfalls occur are identified. A clear inverse relationship between the relief and slope angles of historical rockfalls enables us to calculate the rock mass strength of the landscape scale by the Culmann model and the relief–slope angle relationship curve. Other parameters used in our proposed approach are calculated by ArcGIS and statistic tools. By applying our approach, the potential rockfall source areas in the study are identified and further zoned into three susceptibility classes that could be used as a reference for a regional rockfall susceptibility study. Using the space partition of historical rockfall inventory, our prediction result is validated. Most of the rockfall source areas (i.e., 71.92%) identified in the validation area are occupied by historical rockfalls, which proves the good prediction of our approach. The dominant uncertainty in this paper is derived from the process of calculating rock mass strength, defining the specific area for searching potential rockfall source areas, and the resolution of the digital elevation model.

Graphical Abstract

1. Introduction

The identification of potential rockfall source areas is the first step in assessing rockfall susceptibility, hazard, risk, and determining rockfall disaster prevention and mitigation [1,2,3,4,5,6,7]. However, it is very difficult work to carry out in mountain areas, especially the steep and high-relief slopes, which are not accessible on site. Field investigation provides the most effective method to survey the distribution of potential rockfall source areas at a specific site [8]. Through field investigation, the engineering geology conditions controlling rockfall distribution, including the rock mass strength, orientation of structures, joint density, slope angle, relief, and the activity of tectonic faults, could be carefully studied on site [7,9,10]. Recent technologies, including unmanned aerial vehicles, terrestrial laser scanning, monitoring systems, photogrammetry, and point cloud analysis software tools (e.g., AgiSoft, Photoscan, and Coltop) [11,12,13,14], help researchers to acquire detailed information of the above conditions.
For the identification of potential rockfall source areas at a regional scale, the traditional field investigation is not as effective as that at a specific site because of its limited investigation scope and because it consumes much time and human resources. Hence, some researchers have developed regional rockfall susceptibility models based on ArcGIS to estimate the spatial distribution of rockfall using the causal factors of rockfall including lithology, terrain, elevation, faults, rainfall, and earthquakes [15,16,17]. Based on the results of rockfall susceptibility assessments, the whole area could be zoned into different areas with different susceptible degrees, which provides useful guidance for regional land use and rockfall disaster prevention plans. Alternatively, some researchers have identified rockfall sources at a regional scale by remote sensing interpretation technologies using multi-temporal aerial photos, helicopter-based remote sensing imagery, and high-resolution digital elevation model (DEM) [5,6,18]. Subsequently, regional locations of rockfall source areas could be identified for further rockfall kinematic modeling and predicting regional rockfall hazards [19].
Rock mass strength is thought to be the basic controlling factor of slope stability [20,21]. A rock slope with a low rock mass strength fractured by different types of fractures is prone to rockfall [7,22]. Based on Culmann’s two-dimensional slope stability model [23], a hillslope is susceptible to rockfall if its relief, slope angle, or both are larger than the threshold values [24,25,26,27]. This means that the relief and slope angles are a pair of important indicators that could be used to identify rockfall source areas on the slopes whose stability is dominantly controlled by the rock mass strength. Until now, previous studies rarely focused on the approaches combining the relief and slope angles to identify rockfall source areas controlled by the rock mass strength at a regional scale [5,28]. Focusing on this issue, comprehensive technologies, including helicopter-based remote sensing imagery, a DEM with 10 m resolution, images from Google Earth, and field work, were adopted in this study. Lastly, a new approach, the procedures, and the application criteria for identifying rockfall source areas at a regional scale are proposed and were applied in the study area.

2. Materials and Methods

The Wolong (WL) area of Tibet where ideal geological conditions exist to investigate the characteristics of rock mass outcrops and the distribution of rockfall source areas was selected for this study (Figure 1). To identify potential rockfall source areas at a regional scale, a new approach combining the relief and slope angles based on the rock mass strength is proposed.

2.1. Study Area

The lithology of the area is mainly diorite and granite with a small component of gneiss [29,30]. Because of the steep terrain and the widely distributed tectonic structures, the main type of slope failure in the WL area is rockfall. Based on the rockfall scars left on the steep slopes and the rockfall deposits distributed widely, rockfalls in the study area occur frequently, and the dominant kinematic failure modes are toppling and planer sliding (Figure 2).
In this study, based on a power grid project, the helicopter-based remote sensing imagery obtained in 2017 and the 10 m resolution DEM over the study area were acquired with the help of the research group State Grid Corporation of China. Using the data of the complete study area, images of Google Earth, and field work, a rockfall inventory including 235 rockfall scars on bedrocks and 109 rockfall deposits was prepared (Figure 1). The rockfall scars were identified based on the fresh bedrock color left on the scars (Figure 2). Because many of the rockfall deposits were removed by the Yarlung Tsangpo River, and some of the rockfalls are adjacent, it was not possible to track each rockfall’s deposits. Rockfall deposits at the foot of the slopes were identified based on the shape of the deposit (e.g., pyramid) and the identifiable rockfall blocks (e.g., meters) left on the deposits.

2.2. Framework of Identifying Potential Rockfall Source Areas

By using ArcGIS and statistic tools, our new approach proposed in this paper includes the framework as follows:
(1) Calculate the rock mass strength parameters (c1, φ1) of the bedrock slope at the landscape scale and build the relief–slope (R–S) angle relationship curve in the study area (Figure 3a).
(2) Measure the present relief (H) and slope angle (β) of each specific area (A) (Figure 3b) with the potential rockfall over the slope areas.
(3) Calculate the limit relief ( H c ) of the specific area (A) by the Culmann model [23] (Figure 3a). The Culmann model indicates that the relief of the slope is controlled by the rock mass strength, and the slope angle (β) has the following relationship with the limit relief (maximum height) (Hc):
H c = 4 C γ s i n β c o s φ [ 1 cos ( β φ ) ]
where c is cohesion, γ is the bulk density, and φ is the internal friction angle.
(4) Compare the present relief (H) and the limit relief ( H c ) of the specific area to estimate its state of stability. When the present relief of the bedrock is larger than the limit relief (Figure 3b), the bedrock is prone to generate rockfalls. Hence, the area whose relief exceeds its limit relief is identified as the potential rockfall source area.
(5) Using the procedure from steps (1) to (4), all areas of the slopes over the study area are searched, and their states of stability are estimated. Eventually, all the potential rockfall source are identified.
Different from the availability of rock mass strength tested in the laboratory, the rock mass strength of bedrocks at the landscape scale is very difficult to test on site because of the lack of suitable approaches [20,21]. Schmidt and Montgomery [24] proposed an approach to estimate the rock mass strength parameters (c, φ) using the relief and slope angles of historical rockfalls with the Culmann model. Previous studies have applied the approach to calculate the parameters of rock mass strength using the relief and slope angles of historical landsides or rockfall scars in some cases [25,26,27,31]. In reference to previous studies, data of the relief and slope angles of 235 historical rockfall scars were first extracted by ArcGIS. Then, the parameters of the rock mass strength (c, φ) at the landscape scale in our study area were calculated under the precondition that the bedrock relief is controlled by the rock mass strength.

3. Results

The geometrical characteristics of historical rockfalls were analyzed to define the specific area (A) used in the new approach for searching rockfall source areas. Applying our proposed approach, the rockfall source areas were identified in the study area and were zoned into three susceptibility classes.

3.1. Geometrical Characteristics of Historical Rockfalls

Before calculating the relief and slope angles of each specific area, a suitable value of the specific area (A) should be first determined. To define the value of A in this study, the geometrical characteristics of 235 historical rockfalls were first analyzed by ArcGIS. The values of the relief and slope angles and the areas of historical rockfalls were measured separately. According to the statistical results, the relief of historical rockfalls is mainly distributed between 40 and 130 m, the slope angle is generally larger than 45°, and the area of each historical rockfall scar is generally less than 9000 m2. All three groups of data show Gaussian distribution characteristics (Figure 4). The mean area of historical rockfalls is 5217 m2, which was adopted as the specific area (A) for searching rockfall source areas of the slopes over the study area.

3.2. Locations of Identified Rockfall Sources

A clear inverse relationship between the relief and slope angles of 235 historical rockfalls (Figure 5) enabled us to calculate the rock mass strength at the landscape scale in this study. According to previous studies [25,26,27,31], the minimum and maximum parameters of rock mass strength can be estimated using the lower envelope and the upper envelope of the R–S curves obtained by data of the relief and slope angles of failed slopes. The upper envelope of the R–S curve (Figure 5) represents the maximum strength of the rock mass, and the lower envelope represents the minimum.
Applying the results of 235 historical rockfalls to Equation (1) by the optimization algorithm (Figure 5), the minimum and maximum rock mass cohesions (c) in the study area are 28 Pa and 270 kPa, respectively, and both internal friction angles (φ) are 23°. Using the upper envelope and the mean value of the upper envelope and the lower envelope, with the lower envelope corresponding to each slope angle, three R–S relationship curves were built as the threshold to determine if each specific area is stable or unstable (Figure 5). The R–S curves could be regarded as three different estimates: the aggressive estimate, the moderate estimate, and the conservative estimate, corresponding to the upper envelope, the mean value, and the lower envelope curves, which could be selected by different aims or rockfall disaster prevention and mitigation strategies.
Based on our approach, to identify the rockfall source area is to identify the area whose present relief exceeds the limit relief corresponding to its slope angle (Figure 3). Hence, we calculated the limit relief corresponding to each slope angle (Table 1) using the three R–S curves presented in Figure 5. Considering the actual slope distribution of the historical rockfalls (Figure 4b), we mainly focused on the slopes with angles larger than 45°. Comparing the present relief and the limit relief in each specific area of the slopes in the study area, we obtained the rockfall source areas in ArcGIS (Figure 6). More rockfall source areas are distributed at the lower parts of the slopes whose slope angles are relatively bigger, which is probably affected by the intense incision of the Yarlung Tsangpo River [32].

3.3. Zoning Map of Rockfall Susceptibility

In rockfall risk analysis, the rockfall susceptibility assessment is the first step to carry out [16,28,33]. The rockfall susceptibility map helps to highlight the spatial distribution of potentially unstable slopes [34], which is usually zoned into different susceptibility classes, e.g., high, medium, and low susceptibility, to represent different rockfall susceptible degrees of slopes. Based on the sketch in Figure 3, if H > Hc, the area A on the slope is unstable, which means it is prone to rockfall in the future. The larger H is than Hc, the higher the possibility of rockfall in area A, and hence the higher rockfall susceptibility of area A. To provide a reference for the susceptibility study following the identification of the rockfall source areas by our approach, the rockfall source areas were zoned into three susceptibility classes, high (I), medium (II), and low susceptibility (III) areas, using the upper envelope, the mean envelope, and the lower envelope in Figure 5. In this way, the regional rockfall susceptibility maps were produced in the study area (Figure 6).
The slope angle and the elevation of all the susceptibility classes were analyzed, and their frequency distributions were obtained. According to the results (Figure 7), the rockfall source areas with different susceptibility classes have different ranges of slope angle. The rockfall source areas within the high susceptibility class are mainly distributed on the slopes with the angles of 60–66°, those of medium susceptibility are distributed on the slopes with the angles of 54–61°, and those of low susceptibility are distributed on the slopes with the angles of 46–55°. However, based on the distribution statistics of elevation in Figure 7, no obvious relationship between elevation and susceptibility classes was observed in this study.

4. Discussion

It is important to validate the prediction results in prediction modeling [35]. Hence, the validation of our proposed new approach is discussed here. The uncertainty in our new approach is mainly derived from the process of calculating the rock mass strength parameters, defining the specific area (A) used for searching the rockfall source areas, and the resolution of the DEM.

4.1. Validation

Proper validation should be based on the comparison between the prediction results and the actual characteristics affected by future rockfalls [16,35]. The selection of approaches depends on the characteristics of dataset. In this study, the space partition to separate our rockfall inventory into two groups was chosen since information was lacking related to the time of the 235 historical rockfalls (Figure 8). To validate the prediction result, the rockfall inventory was partitioned into two groups. One group was used for prediction (Figure 1a) and the other was used for validation (Figure 1b).
The first group with 118 historical rockfalls (Figure 1a) was adopted to calculate the rock mass strength and build the R–S relationship curve. According to the fitting result of the selected 118 rockfall sources, the minimum and maximum rock mass cohesions (c) are 28 and 230 kPa, respectively, and both internal friction angles (φ) are 22°. To present a common case of application, we used the mean value of the upper envelope and the lower envelope to build the R–S relationship curve in this paper. Based on the procedure of our approach (Figure 3), the rockfall source areas were obtained. The second group with 117 historical rockfalls (Figure 1b) was used to validate the prediction result (Figure 8). The validation result demonstrated that 117 historical rockfalls occupy 71.92% of the rockfall source areas in the validation area predicted by our new approach, which proves to be a good prediction.
Besides, the prediction result was validated with field work. Most of the rockfall source areas identified by our approach are distributed on slopes with high relief and steep terrain. This is consistent with the observation of rockfall distribution on the maps of DEM and hillshade (Figure 8). Most of the historical rockfall deposits are distributed at the foot of the slopes with identified rockfall source areas, proving that the identified rockfall source areas are distributed rockfalls that occurred in the past, and that unstable slopes are prone to rockfall in the future [19].

4.2. Uncertainty Analysis

An accurate calculation of rock mass strength parameters is the fundamental work in our new proposed approach. For example, for a slope with a specific relief (H) and slope angle (β), the bigger the rock mass strength estimated, the higher the value of the limit relief (Hc) calculated by Equation (1), and hence the bigger the difference between Hc and H. For the stability of the slope, more stable (i.e., the lower rockfall susceptibility) slopes were predicted. In other words, if the rock mass strength parameters are overestimated, the slope would be prone to rockfall with an incorrect prediction of low susceptibility. In our new approach, the more accurate estimated rock mass strength parameters will greatly improve the accuracy of the R–S relationship curve, and hence the prediction results of rockfall source areas. However, it is very difficult to quantitatively estimate the rock mass strength at the landscape. An in situ test is recognized as one of the most reliable methods to obtain rock mass strength [10,20], which proves to be very difficult to carry out in high mountain areas. In this paper, the back analysis for rock mass strength using historical rockfalls based on the Culmann model is adopted, whose accuracy depends on the reliability of the identification, the boundary of historical rockfalls, and the uncertainty in fitting the data.
According to the procedure of our approach and the sketch in Figure 3, the specific area (A) is also important in affecting the accuracy of the prediction result. For example, if a smaller area than the real one was defined to search the potential rockfall source areas, a lower relief (H) than the real value of a potential source area would be obtained. For the same slope angle and corresponding estimated Hc, the area would mistakenly be regarded as more stable (Figure 3) than its real stability state. In this paper, the mean area of historical rockfalls is defined as the specific area (A) for searching the potential rockfall source areas. The definition of the parameter could be further studied in the future.
The resolution of DEM plays an important role in controlling the accuracy of the result [5] because it is the basic data in almost every data process of our new approach. The higher the resolution of the DEM, the higher the accuracy of the values of the relief and slope angles of historical rockfalls and the potential rockfall source areas. However, it is not easy to acquire a high-quality DEM in a large study area currently.

5. Conclusions

The main type of rock mass failure in the Wolong area of Tibet is rockfall. Using data from helicopter-based remote sensing imagery, a DEM with 10 m resolution of the study area, images from Google Earth, and field work, a rockfall inventory including 235 rockfalls scars on bedrocks and 109 rockfall deposits was prepared. According to the statistical results, the relief of historical rockfalls is mainly distributed between 40 and 130 m, the slope angle is generally larger than 45°, and the area of each historical rockfall scar is generally less than 9000 m2. A clear inverse relationship between the relief and slope angles of historical rockfalls enabled us to calculate the rock mass strength at the landscape scale base on the Culmann model, obtaining the minimum and maximum rock mass cohesions (c) in the study area from 28 to 270 kPa, respectively, and the internal friction angle of 23°.
Required by the actual needs of identification on high and steep slopes, this paper proposes a new approach using the relief–slope angle relationship to identify the rockfall source areas controlled by the rock mass strength on a regional scale. Based on historical rockfalls and a high-resolution DEM, we obtained the parameters used in our proposed approach. By applying our approach, the potential rockfall source areas in the study area were identified and further zoned into three susceptibility classes that could be used as a reference for the study of regional rockfall susceptibility assessment. According to the results, rockfall source areas within the high susceptibility class are mainly distributed on the slopes with the angles of 60–66°, those of medium susceptibility are distributed on the slopes with the angles of 54–61°, and those of low susceptibility are distributed on the slopes with the angles of 46–55°.
By the space partition and the field work, our prediction result was validated. Most of the rockfall source areas (i.e., 71.92%) identified in the validation area are occupied by actual historical rockfalls, which proves the accurate prediction ability of our approach. The locations of the rockfall source areas obtained in this paper could provide reference for actual rockfall disaster prevention and mitigation in the study area. Our proposed approach could be used to identify the rockfall source areas in the regional areas that are not accessible. In the paper, the dominant uncertainty is derived from the process of calculating the rock mass strength parameters, the process of defining the specific area (A) that is used for searching the rockfall source areas, and the resolution of the DEM. Many more studies estimating the rock mass strength at the landscape scale and defining the specific area (A) are necessary in the future.

Author Contributions

Conceptualization, X.W.; methodology X.W., H.L.; formal analysis, H.L., J.S.; writing—original draft preparation, X.W., H.L.; writing—review and editing, X.W.; visualization, H.L.; supervision, X.W.; project administration, J.S.; funding acquisition, X.W. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Second Tibetan Plateau Scientific Expedition and Research Program (STEP) (Grant No. 2019QZKK0904), the Strategic Priority Research Program of the Chinese Academy of Sciences (Grant No. XDA23090402), the Application of Synthetic Aperture Radar-Based Geological Hazard Analysis Technology on the Strategic Electricity Transmission Passage of Sichuan-Tibet Plateau (Grant No. 52199918000C).

Acknowledgments

We thank the program research group of the State Grid Corporation of China for providing the 10 m resolution topographic DEM data.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Mazzoccola, D.F.; Hudson, J.A. A comprehensive method of rock mass characterization for indicating natural slope instability. Q. J. Eng. Geol. Hydrogeol. 1996, 29, 37–56. [Google Scholar] [CrossRef]
  2. Bunce, C.M.; Cruden, D.M.; Morgenstern, N.R. Assessment of the hazard from rock fall on a highway. Can. Geotech. J. 1997, 34, 344–356. [Google Scholar] [CrossRef]
  3. Guzzetti, F.; Carrara, A.; Cardinali, M.; Reichenbach, P. Landslide hazard evaluation: A review of current techniques and their application in a multi-scale study, Central Italy. Geomorphology 1999, 31, 181–216. [Google Scholar] [CrossRef]
  4. Gokceoglu, C.; Sonmez, H.; Ercanoglu, M. Discontinuity controlled probabilistic slope failure risk maps of the Altindag (settlement) region in Turkey. Eng. Geol. 2000, 55, 277–296. [Google Scholar] [CrossRef]
  5. Loye, A.; Jaboyedoff, M.; Pedrazzini, A. Identification of potential rockfall source areas at a regional scale using a DEM-based geomorphometric analysis. Nat. Hazards Earth Syst. Sci. 2009, 9, 1643–1653. [Google Scholar] [CrossRef]
  6. Liu, H.; Wang, X.; Liao, X.; Sun, J.; Zhang, S. Rockfall Investigation and Hazard Assessment from Nang County to Jiacha County in Tibet. Appl. Sci. 2019, 10, 247. [Google Scholar] [CrossRef] [Green Version]
  7. Wang, X.; Frattini, P.; Stead, D.; Sun, J.; Liu, H.; Valagussa, A.; Li, L. Dynamic rockfall risk analysis. Eng. Geol. 2020, 272, 105622. [Google Scholar] [CrossRef]
  8. Wang, X.; Zhang, L.; Wang, S.; Agliardi, F.; Frattini, P.; Crosta, G.B.; Yang, Z. Field investigation and rockfall hazard zonation at the Shijing Mountains Sutra caves cultural heritage (China). Environ. Earth Sci. 2011, 66, 1897–1908. [Google Scholar] [CrossRef]
  9. Stead, D.; Wolter, A. A critical review of rock slope failure mechanisms: The importance of structural geology. J. Struct. Geol. 2015, 74, 1–23. [Google Scholar] [CrossRef]
  10. Saroglou, C.; Qi, S.; Guo, S.; Wu, F. ARMR, a new classification system for the rating of anisotropic rock masses. Bull. Int. Assoc. Eng. Geol. 2019, 78, 3611–3626. [Google Scholar] [CrossRef]
  11. Jaboyedoff, M.; Metzger, R.; Oppikofer, T.; Couture, R.; Derron, M.H.; Locat, J.; Turmel, D. New insight techniques to analyze rock-slope relief using DEM and 3D-imaging cloud points: COLTOP 3D software. In Proceedings of the 1st Canada-US Rock mechanics Symposium, Vancouver, BC, Canada, 27–31 May 2007; pp. 27–31. [Google Scholar] [CrossRef]
  12. DiMatteo, L.; Romeo, S.; Kieffer, D.S. Rockfall analysis in Alpine area by using a reliable integrated monitoring system: Results from the Ingelsberg slope (Salzburg Land, Austria). Bull. Eng. Geol. Environ. 2017, 76, 413–420. [Google Scholar] [CrossRef]
  13. Mazzanti, P.; Schilirò, L.; Martino, S.; Antonielli, B.; Brizi, E.; Brunetti, A.; Margottini, C.; Mugnozza, G.S. The Contribution of Terrestrial Laser Scanning to the Analysis of Cliff Slope Stability in Sugano (Central Italy). Remote Sens. 2018, 10, 1475. [Google Scholar] [CrossRef] [Green Version]
  14. Kovanič, Ľ.; Blistan, P.; Urban, R.; Štroner, M.; Blišťanová, M.; Bartoš, K.; Pukanská, K. Analysis of the Suitability of High-Resolution DEM Obtained Using ALS and UAS (SfM) for the Identification of Changes and Monitoring the Development of Selected Geohazards in the Alpine Environment—A Case Study in High Tatras, Slovakia. Remote Sens. 2020, 12, 3901. [Google Scholar] [CrossRef]
  15. Cascini, L. Applicability of landslide susceptibility and hazard zoning at different scales. Eng. Geol. 2008, 102, 164–177. [Google Scholar] [CrossRef]
  16. Wang, X.L.; Zhang, L.Q.; Wang, S.J.; Lari, S. Regional landslide susceptibility zoning with considering the aggregation of landslide points and the weights of factors. Landslides 2014, 11, 399–409. [Google Scholar] [CrossRef]
  17. Budetta, P.; De Luca, C.; Nappi, M. Quantitative rockfall risk assessment for an important road by means of the rockfall risk management (RO.MA.) method. Bull. Int. Assoc. Eng. Geol. 2016, 75, 1377–1397. [Google Scholar] [CrossRef]
  18. Duarte, R.M.; Marquínez, J. The influence of environmental and lithologic factors on rockfall at a regional scale: An evaluation using GIS. Geomorphology 2002, 43, 117–136. [Google Scholar] [CrossRef]
  19. Frattini, P.; Crosta, G.; Carrara, A.; Agliardi, F. Assessment of rockfall susceptibility by integrating statistical and physically-based approaches. Geomorphology 2008, 94, 419–437. [Google Scholar] [CrossRef]
  20. Hoek, E. Strength of jointed rock masses. Géotechnique 1983, 33, 187–223. [Google Scholar] [CrossRef] [Green Version]
  21. Gudmundsson, A. Rock Fractures in Geological Processes; Cambridge University Press (CUP): Cambridge, UK, 2011. [Google Scholar]
  22. Francioni, M.; Antonaci, F.; Sciarra, N.; Robiati, C.; Coggan, J.; Stead, D.; Calamita, F. Application of Unmanned Aerial Vehicle Data and Discrete Fracture Network Models for Improved Rockfall Simulations. Remote Sens. 2020, 12, 2053. [Google Scholar] [CrossRef]
  23. Culmann, C. Die Graphische Statik; Meyer and Zeller: Zurich, Switzerland, 1875; 644p. [Google Scholar]
  24. Schmidt, K.M.; Montgomery, D.R. Limits to Relief. Science 1995, 270, 617–620. [Google Scholar] [CrossRef]
  25. Montgomery, D.R.; Brandon, M.T. Topographic controls on erosion rates in tectonically active mountain ranges. Earth Planet. Sci. Lett. 2002, 201, 481–489. [Google Scholar] [CrossRef]
  26. Crosta, G.B.; Utili, S.; De Blasio, F.V.; Castellanza, R. Reassessing rock mass properties and slope instability triggering conditions in Valles Marineris, Mars. Earth Planet. Sci. Lett. 2014, 388, 329–342. [Google Scholar] [CrossRef]
  27. DiBiase, R.A.; Rossi, M.W.; Neely, A.B. Fracture density and grain size controls on the relief structure of bedrock landscapes. Geology 2018, 46, 399–402. [Google Scholar] [CrossRef] [Green Version]
  28. Crosta, G.B.; Agliardi, F. A methodology for physically based rockfall hazard assessment. Nat. Hazards Earth Syst. Sci. 2003, 3, 407–422. [Google Scholar] [CrossRef]
  29. Quidelleur, X.; Grove, M.; Lovera, O.M.; Harrison, T.M.; Yin, A.; Ryerson, F.J. Thermal evolution and slip history of the Renbu Zedong Thrust, southeastern Tibet. J. Geophys. Res. Space Phys. 1997, 102, 2659–2679. [Google Scholar] [CrossRef]
  30. Xu, Z.-Q.; Dilek, Y.; Yang, J.-S.; Liang, F.-H.; Liu, F.; Ba, D.-Z.; Cai, Z.-H.; Li, G.-W.; Dong, H.-W.; Ji, S.-C. Crustal structure of the Indus–Tsangpo suture zone and its ophiolites in southern Tibet. Gondwana Res. 2015, 27, 507–524. [Google Scholar] [CrossRef]
  31. Burbank, D.W.; Leland, J.; Fielding, E.; Anderson, R.S.; Brozovic, N.; Reid, M.R.; Duncan, C. Bedrock incision, rock uplift and threshold hillslopes in the northwestern Himalayas. Nat. Cell Biol. 1996, 379, 505–510. [Google Scholar] [CrossRef]
  32. Finnegan, N.J.; Hallet, B.; Montgomery, D.R.; Zeitler, P.K.; Stone, J.O.; Anders, A.M.; Liu, Y.P. Coupling of rock uplift and river incision in the Namche Barwa–Gyala Peri massif, Tibet. Geol. Soc. Am. Bull. 2008, 120, 142–155. [Google Scholar] [CrossRef]
  33. Fell, R.; Cororninas, J.; Bonnard, C.; Cascini, L.; Leroi, E.; Savage, W.Z. Guidelines for landslide susceptibility, hazard and risk-zoning for land use planning. Eng. Geol. 2008, 102, 85–98. [Google Scholar] [CrossRef] [Green Version]
  34. Ayalew, L.; Yamagishi, H.; Maruib, H.; Kannoc, T. Landslides in Sado Island of Japan: Part II. GIS-based susceptibility mapping with comparisons of results from two methods and verifications. Eng. Geol. 2005, 81, 432–445. [Google Scholar] [CrossRef]
  35. Chung, C.F.; Fabbri, A.G. Validation of spatial prediction models for landslide hazard mapping. Nat. Hazard. 2003, 30, 451–472. [Google Scholar] [CrossRef]
Figure 1. Simplified geology map of the study area. Two zones are separated by the space partition for validation of the prediction result, with (a) the back calculation zone and (b) the verification zone. (c) Sample area presented in Figure 2.
Figure 1. Simplified geology map of the study area. Two zones are separated by the space partition for validation of the prediction result, with (a) the back calculation zone and (b) the verification zone. (c) Sample area presented in Figure 2.
Remotesensing 13 00938 g001
Figure 2. Samples of the rockfall scars and deposits on the helicopter-based remote sensing imagery with the location c in Figure 1b and the historical rockfalls with kinematic failure modes of toppling and planer sliding and their deposits at the foot of the slopes.
Figure 2. Samples of the rockfall scars and deposits on the helicopter-based remote sensing imagery with the location c in Figure 1b and the historical rockfalls with kinematic failure modes of toppling and planer sliding and their deposits at the foot of the slopes.
Remotesensing 13 00938 g002
Figure 3. Sketch map describing the approach. (a) Slope angle versus relief for hillslope. (b) Sketch map of estimating the state of stability for the specific area (A).
Figure 3. Sketch map describing the approach. (a) Slope angle versus relief for hillslope. (b) Sketch map of estimating the state of stability for the specific area (A).
Remotesensing 13 00938 g003
Figure 4. The geometrical statistics of historical rockfalls.
Figure 4. The geometrical statistics of historical rockfalls.
Remotesensing 13 00938 g004
Figure 5. Rock mass strength fitting result and the relief–slope (R–S) angle relationship in the study area.
Figure 5. Rock mass strength fitting result and the relief–slope (R–S) angle relationship in the study area.
Remotesensing 13 00938 g005
Figure 6. The identification results of the rockfall source areas in the study area. The results of (ac) correspond to high (I), medium (II) and low (III) rockfall susceptibility classes, respectively.
Figure 6. The identification results of the rockfall source areas in the study area. The results of (ac) correspond to high (I), medium (II) and low (III) rockfall susceptibility classes, respectively.
Remotesensing 13 00938 g006aRemotesensing 13 00938 g006b
Figure 7. The distribution of slope and elevation of the identified rockfall source areas within different susceptibility classes.
Figure 7. The distribution of slope and elevation of the identified rockfall source areas within different susceptibility classes.
Remotesensing 13 00938 g007aRemotesensing 13 00938 g007b
Figure 8. The validation of the prediction result obtained by using one group of historical rockfalls separated by the space partition (Figure 1a).
Figure 8. The validation of the prediction result obtained by using one group of historical rockfalls separated by the space partition (Figure 1a).
Remotesensing 13 00938 g008
Table 1. Limit relief of the specific areas corresponding to the slope angle.
Table 1. Limit relief of the specific areas corresponding to the slope angle.
Slope
Angle (°)
Limit Relief (m)Slope
Angle (°)
Limit Relief (m)Slope
Angle (°)
Limit Relief (m)
C-UC-MC-LC-UC-MC-LC-UC-MC-L
46.035.5189.1342.761.016.286.2156.276.09.650.992.3
47.033.2176.7320.362.015.682.8150.077.09.349.489.5
48.031.2165.8300.463.015.079.6144.278.09.047.986.9
49.029.3155.8282.464.014.476.6138.879.08.746.584.3
50.027.6146.9266.165.013.973.7133.680.08.545.281.9
51.026.1138.7251.466.013.471.1128.881.08.343.979.6
52.024.7131.3238.067.012.968.5124.282.08.042.777.3
53.023.4124.6225.768.012.466.2119.983.07.841.575.2
54.022.2118.4214.569.012.063.9115.884.07.640.373.1
55.021.2112.7204.170.011.661.8111.985.07.439.271.1
56.020.2107.4194.671.011.259.7108.286.07.238.269.2
57.019.3102.5185.872.010.957.8104.787.07.037.267.3
58.018.498.0177.673.010.555.9101.488.06.836.265.5
59.017.693.8170.074.010.254.298.289.06.635.263.8
60.016.989.9162.975.09.952.595.290.06.534.362.2
Note: for the slope angle, the front part is open interval, and latter part is closed interval. For example, 46 represents the range of (45, 46]. C-U: R–S curve by upper envelope; C-M: R–S curve by mean value; C-L: R–S curve by lower envelope.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, X.; Liu, H.; Sun, J. A New Approach for Identification of Potential Rockfall Source Areas Controlled by Rock Mass Strength at a Regional Scale. Remote Sens. 2021, 13, 938. https://0-doi-org.brum.beds.ac.uk/10.3390/rs13050938

AMA Style

Wang X, Liu H, Sun J. A New Approach for Identification of Potential Rockfall Source Areas Controlled by Rock Mass Strength at a Regional Scale. Remote Sensing. 2021; 13(5):938. https://0-doi-org.brum.beds.ac.uk/10.3390/rs13050938

Chicago/Turabian Style

Wang, Xueliang, Haiyang Liu, and Juanjuan Sun. 2021. "A New Approach for Identification of Potential Rockfall Source Areas Controlled by Rock Mass Strength at a Regional Scale" Remote Sensing 13, no. 5: 938. https://0-doi-org.brum.beds.ac.uk/10.3390/rs13050938

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