Spatial distribution of block falls using volumetric GIS–decision-tree models

https://doi.org/10.1016/j.jag.2010.05.008Get rights and content

Abstract

Block falls are considered a significant aspect of surficial instability contributing to losses in land and socio-economic aspects through their damaging effects to natural and human environments. This paper predicts and maps the geographic distribution and volumes of block falls in central Lebanon using remote sensing, geographic information systems (GIS) and decision-tree modeling (un-pruned and pruned trees). Eleven terrain parameters (lithology, proximity to fault line, karst type, soil type, distance to drainage line, elevation, slope gradient, slope aspect, slope curvature, land cover/use, and proximity to roads) were generated to statistically explain the occurrence of block falls. The latter were discriminated using SPOT4 satellite imageries, and their dimensions were determined during field surveys. The un-pruned tree model based on all considered parameters explained 86% of the variability in field block fall measurements. Once pruned, it classifies 50% in block falls’ volumes by selecting just four parameters (lithology, slope gradient, soil type, and land cover/use). Both tree models (un-pruned and pruned) were converted to quantitative 1:50,000 block falls’ maps with different classes; starting from Nil (no block falls) to more than 4000 m3. These maps are fairly matching with coincidence value equal to 45%; however, both can be used to prioritize the choice of specific zones for further measurement and modeling, as well as for land-use management. The proposed tree models are relatively simple, and may also be applied to other areas (i.e. the choice of un-pruned or pruned model is related to the availability of terrain parameters in a given area).

Introduction

Block falls represent the most intensely geomorphic process of the cliff zone in mountainous areas (Luckman and Fiske, 1995). They are widely distributed in the Mediterranean region (including Lebanon), affecting life and properties severely. In Lebanon, estimations over the last 40 years indicate that mass movements (including block falls) cost 10–15 millions of US dollars per year causing numerous fatalities and injuries (Khawlie, 2000). Concrete examples can give an idea about the importance of mass movement occurrence damages in this small country (Abdallah, 2007). A mass movement in Mazrait Bani Assaf (central Lebanon) caused tremendous damages, i.e. complete destruction of several farmlands, partial destruction of the irrigation system, and occurrence of some fractures in the earth next to some houses with huge resettlement of the ground at the main road access to the village. On 18 March 1992, a huge movement occurred in Hazerta near Zahlé (in the Bekaa valley) which destroyed several houses and 20 people died. On 25th March 2000, a mass movement destroyed 10 houses and several ski resorts in Aayoun Orghouch (Mount Lebanon). Moreover, dense road networks often pass through rugged and unstable topographic features and hence are affected by block falls with high impact on the population and environment. Protection plans to prevent block falls are limited and are often not adequate to meet general safety requirements.

A number of methodologies for mass movement hazard/susceptibility zonation over large areas exist in the literature, which are belonging to qualitative, semi-quantitative and quantitative GIS-assisted approaches. Each of these approaches has its value for certain applications and disadvantages for other objectives, given the specific conditions of the Mediterranean environments. In this context, the GIS technology has been applied to examine landslide hazard in several Mediterranean countries, such as Jordan (Al-Homoud and Masanat, 1998, Husein et al., 2000); northern Spain (Duarte and Marquinez, 2002); Tuscany, Italy (Luzi et al., 2000), central Italy (Carrara et al., 1991), Southern Italy (Budetta et al., 1994, Refice and Capolongo, 2002), and Fabriano, Italy (Luzi and Pergalani, 1996).

Studies of rockfalls are often based on field surveys (Selby, 1980, Kirkaldie, 1985, Romana, 1985, Bieniawski, 1993, Habibagahi and Katebi, 2004, Lindsay et al., 2001), and hazard is estimated either by an empirical assessment of the susceptibility to failure, or by the calculation of a safety factor derived from models of rock mechanics (Carere et al., 2001, Crosta et al., 2001, Irigaray et al., 2003). Very few statistical/deterministic GIS models have been developed for predictive susceptibility modeling of rock falls at regional scale (Abbot et al., 1998, Ayala-Carcedo et al., 2003, Marquinez et al., 2003, Baillifard et al., 2004, Calcaterra et al., 2004, Guzzetti et al., 2004). They are essentially classifying a given area into units according to hazard degree based on information obtained from the intercorrelation between landscape factors and mass movement distribution. Our approach provides a systematic decision-tree model for quantifying and predicting volumes of block falls in a representative area of Lebanon, based on analyzing terrain parameters, in the frame of geographic information systems (GIS). The proposed model comprises a set of rules to classify (predict) a dependent target variable (volume of block falls in m3) using the values of independent predictor variables (terrain parameters). These predictors are both continuous (e.g. elevation) and categorical (e.g. lithology). The obtained block fall maps from converting various decision-tree models, at 1:50,000 cartographic scale, can serve as geological inventories useful for land-use management and environmental decision-making.

Section snippets

Study area

The study area (762 km2), located in the central part of the Lebanese Mediterranean slopes, is chosen because it comprises a mixture of lithological, morphological, pedological, hydrological and climatic diversity that witnesses major block falls. The area is bounded to the west by the Mediterranean Sea, to the north-west by Biaqout village, and to the south by Damour river (Fig. 1).

The outcropping stratigraphic sequence exposes rock formations spanning from the middle Jurassic to Recent.

Materials and methods

The volumetric mapping of block falls distribution was realized in several steps (Fig. 2). SPOT4 (HRVIR) satellite imageries were used together with field surveys for delineating existing block falls, and determining their dimensions. The obtained block fall layer information (point location) was then intersected with maps of predictor terrain parameters. Several decision-tree models were explored on the result of this intersection combining block falls locations and the corresponding terrain

Character of block falls’ occurrence

The sixty-two detected block falls were ranged from 234 to 5436 m3, with a mean volume being equal to 2139 m3 and a standard deviation of 1959. Most of them occur under the following conditions: on elevations ranging between 500 and 1000 m (37% of block falls), with moderately steep slopes oscillating between 12 and 18° (48%), facing southward (16%), in concave (50%) forested areas of coniferous type (21%), growing over highly fissured and jointed dolomites and dolomitic limestone of the Jurassic

Conclusion

The established decision-tree models enabled, for the first time, mapping of predicted block falls’ volumes in a region of Lebanon at a scale of 1:50,000, based on geo-environmental characteristics (e.g. topography, geology, soil and land cover/use). The modeling approach was easily implemented with available GIS software, and is suitable for data exploration and predictive block falls’ volumetric mapping. It is explicit and can be critically evaluated and revised when necessary. It can be also

Acknowledgments

This research is a part of a project entitled “Establishing Mass Movements Susceptibility and Hazard maps for the Lebanese territory using Remote Sensing and GIS techniques” (2008–2010) funded by the Lebanese National Council for Scientific Research (CNRS). We acknowledge all responsible at the CNRS. We also thank the editor and reviewers for helpful reviews of the manuscript.

References (56)

  • A.S. Al-Homoud et al.

    A classification system for the assessment of slope stability of terrain along highway routes in Jordan

    Environ. Geol.

    (1998)
  • F.J. Ayala-Carcedo et al.

    Large scales rock fall reach susceptibility maps in La Cabrera Sierra (Madrid) performed with GIS and dynamic analysis at 1:5,000

    Nat. Hazards

    (2003)
  • F. Baillifard et al.

    Towards a GIS-based rockfall hazard assessment along the Quebec City Promontory, Quebec, Canada

  • A.A. Basma et al.

    Modeling soil collapse by artificial neural networks

    Geotech. Geol. Eng.

    (2004)
  • Berk, R.A., 2003. An introduction to ensemble methods for data analysis. UCLA Department of Statistics Technical...
  • Z.T. Bieniawski

    Classification of rock masses for engineering: the R.M.R. system and future trends

  • E. Binaghi et al.

    Prediction of displacements in unstable areas using a neural model

    Nat. Hazards

    (2004)
  • R. Bou Kheir et al.

    Caractérisation morpho-pédologique des zones karstiques du Liban. Sensibilité des sols à l’érosion hydrique

    Sécheresse

    (2003)
  • R. Bou Kheir et al.

    Use of terrain variables for mapping gully erosion susceptibility in Lebanon

    Earth Surf. Proc. Land.

    (2007)
  • L. Breiman

    Decision-tree forests

    Mach. Learn.

    (2001)
  • L. Breiman et al.

    Classification and Regression Trees

    (1984)
  • P. Budetta et al.

    Engineering-geological zoning of potentially unstable rock slopes in Sorrentine Peninsula (South Italy)

  • T.H. Bui et al.

    Toward affective dialogue modeling using partially observable markov decision processes

  • CAL

    Climatic Atlas of Lebanon

    (1982)
  • D. Calcaterra et al.

    Rockfall hazard assessment at Mt. San Costanzo (Sorrento Peninsula, Italy)

  • A. Carrara et al.

    GIS techniques and statistical models in evaluating landslide hazard

    Earth Surf. Process. Land.

    (1991)
  • Carere, K., Ratto, S., Zanolini, F., 2001. Programme Intereg 2c – Falaises. Prévention des mouvements de versants et...
  • G. Crosta et al.

    Valutazione e gestione del rischio da frana

    (2001)
  • Cited by (6)

    • Risk at the margins: A natural hazards perspective on the Syrian refugee crisis in Lebanon

      2019, International Journal of Disaster Risk Reduction
      Citation Excerpt :

      Since the Hyogo Framework for Action [60], the GoL and international actors have made significant investments in landslide risk reduction initiatives throughout the country [25]. Local landslide studies include the susceptibility to rockfall in the central Mount Lebanon district [2], an investigation of rockfall hazard outside of Kfar Nahbrakh [16], and landslide hazard mapping in the Nahr Ibrahim watershed [4]. The Lebanese National Council for Scientific Research (CNRS) released 1/50,000 (25 m effective resolution) scale maps of landslide hazard zones for the entire country based on statistical methods [25,3].

    • Detection and analysis of morphotectonic features utilizing satellite remote sensing and GIS: An example in SW Jordan

      2016, Geomorphology
      Citation Excerpt :

      We achieved this by combining the higher-resolution panchromatic band 8 (spatial resolution of 15 m) with selected lower-resolution multispectral bands (spatial resolution of 30 m) to improve the spatial resolution of the Landsat 8 satellite imagery (Fig. S1 in the supplementary material available online). The latter procedure is known as pan-sharpening technique and is most widely used for improved geological and environmental information extraction (e.g., Schetselaar, 2001; Al Rawashdeh et al., 2006; Pavlic et al., 2008; Zhang, 2008; Abdallah, 2010; Argyriou, 2012). Using the DEM and its derived slope, we conducted longitudinal profile analyses for five streams or channels, with lengths ranging from 25 to 34 km, which flow on the eastern side of the Wadi Araba rift valley and cross a strong topographic transition between the low-gradient uplands and the deeply incised mountain front on the west side of Jordan (Fig. S1).

    • Multimodal method for coseismic landslide hazard assessment

      2016, Engineering Geology
      Citation Excerpt :

      We studied the coseismic landslide hazard across all of Lebanon. Previous work in Lebanon on landslide susceptibility and detection (e.g. Abdallah, 2010; Abdallah et al., 2007; Khawlie and Hassanain, 1984) has focused on the rugged western slopes of the Mount Lebanon Range, where near-vertical limestone cliffs, long talus slopes, and extensive hillslopes of relatively steep (> 20°) terrain cover the landscape. East of Ain Mouaaffaq, near Hammana, Khawlie and Hassanain (1984) identified active and former slope failures involving translational and rotational slides, block failures, toppling, and complex failures, while Abdallah (2010) identified slides and falls of unknown origin across the northwest of Lebanon.

    View full text