GIS and ANN model for landslide susceptibility mapping

Landslide hazard is as the probability of occurrence of a potentially damaging landslide phenomenon within specified period of time and within a given area. The susceptibility map provides the relative spatial probability of landslides occurrence. A study is presented of the application of GIS and artificial neural network model to landslide susceptibility mapping, with particular reference to landslides on natural terrain in this paper. The method has been applied to Lantau Island, the largest outlying island within the territory of Hong Kong. A three-level neural network model was constructed and trained by the back-propagate algorithm in the geographical database of the study area. The data in the database includes digital elevation modal and its derivatives, landslides distribution and their attributes, superficial geological maps, vegetation cover, the raingauges distribution and their 14 years 5-minute observation. Based on field inspection and analysis of correlation between terrain variables and landslides frequency, lithology, vegetation cover, slope gradient, slope aspect, slope curvature, elevation, the characteristic value, the rainstorms corresponding to the landslide, and distance to drainage line are considered to be related to landslide susceptibility in this study. The artificial neural network is then coupled with the ArcView3.2 GIS software to produce the landslide susceptibility map, which classifies the susceptibility into three levels: low, moderate, and high. The results from this study indicate that GIS coupled with artificial neural network model is a flexible and powerful approach to identify the spatial probability of hazards.     

GIS and ANN model for landslide susceptibility mapping

1 IntroductionThe population growth and the expansion of settlements and life-lines over hazardous areas exert increasingly great impact of natural disasters both in the developed and developing countries. In many countries, the economic losses and casualties due to landslides are greater than commonly recognized and generate a yearly loss of property larger than that from any other natural disasters, including earthquakes, floods and windstorms. Landslides in mountainous terrain often occur a…     

Landslide characteristics and slope instability modeling using GIS, Lantau Island, Hong Kong

Steep terrain and high a frequency of tropical rainstorms make landslide occurrence on natural terrain a common phenomenon in Hong Kong. This paper reports on the use of a Geographical Information Systems (GIS) database, compiled primarily from existing digital maps and aerial photographs, to describe the physical characteristics of landslides and the statistical relations of landslide frequency with the physical parameters contributing to the initiation of landslides on Lantau Island in Hong Kong. The horizontal travel length and the angle of reach, defined as the angle of the line connecting the head of the landslide source to the distal margin of the displaced mass, are used to describe runout behavior of landslide mass. For all landslides studied, the horizontal travel length of landslide mass ranges from 5 to 785 m, with a mean value of 43 m, and the average angle of reach is 27.7°. This GIS database is then used to obtain a logistic multiple regression model for predicting slope instability. It is indicated that slope gradient, lithology, elevation, slope aspect, and land-use are statistically significant in predicting slope instability, while slope morphology and proximity to drainage lines are not important and thus excluded from the model. This model is then imported back into the GIS to produce a map of predicted slope instability. The results of this study demonstrate that slope instability can be effectively modeled by using GIS technology and logistic multiple regression analysis.     

Objective landslide detection and surface morphology mapping using high-resolution airborne laser altimetry

A map of extant slope failures is the most basic element of any landslide assessment. Without an accurate inventory of slope instability, it is not possible to analyze the controls on the spatial and temporal patterns of mass movement or the environmental, human, or geomorphic consequences of slides. Landslide inventory maps are tedious to compile, difficult to make in vegetated terrain using conventional techniques, and tend to be subjective. In addition, most landslide inventories simply outline landslide boundaries and do not offer information about landslide mechanics as manifested by internal deformation features. In an alternative approach, we constructed accurate, high-resolution DEMs from airborne laser altimetry (LIDAR) data to characterize a large landslide complex and surrounding terrain near Christchurch, New Zealand. One-dimensional, circular (2-D) and spherical (3-D) statistics are used to map the local topographic roughness in the DEMs over a spatial scale of 1.5 to 10 m. The bedrock landslide is rougher than adjacent unfailed terrain and any of the statistics can be employed to automatically detect and map the overall slide complex. Furthermore, statistics that include a measure of the local variability of aspect successfully delineate four kinematic units within the gently sloping lower half of the slide. Features with a minimum size of surface folds that have a wavelength of about 11 to 12 m and amplitude of about 1 m are readily mapped. Two adjacent earthflows within the landslide complex are distinguished by a contrast in median roughness, and texture and continuity of roughness elements. The less active of the earthflows has a surface morphology that presumably has been smoothed by surface processes. The Laplacian operator also accurately maps the kinematic units and the folds and longitudinal levees within and at the margins of the units. Finally, two-dimensional power spectra analyses are used to quantify how roughness varies with length scale. These results indicate that no dominant length scale of roughness exists for smooth, unfailed terrain. In contrast, zones with different styles of landslide deformation exhibit distinctive spectral peaks that correspond to the scale of deformation features, such as the compression folds. The topographic-based analyses described here may be used to objectively delineate landslide features, generate mechanical inferences about landslide behavior, and evaluate relatively the recent activity of slides.     

A landslide susceptibility model using the Analytical Hierarchy Process method and multivariate statistics in perialpine Slovenia

Landslides cause damage to property and unfortunately pose a threat even to human lives. Good landslide susceptibility, hazard, and risk models could help mitigate or even avoid the unwanted consequences resulted from such hillslope mass movements. For the purpose of landslide susceptibility assessment the study area in the central Slovenia was divided to 78 365 slope units, for which 24 statistical variables were calculated. For the land-use and vegetation data, multi-spectral high-resolution images were merged using Principal Component Analysis method and classified with an unsupervised classification. Using multivariate statistical analysis (factor analysis), the interactions between factors and landslide distribution were tested, and the importance of individual factors for landslide occurrence was defined. The results show that the slope, the lithology, the terrain roughness, and the cover type play important roles in landslide susceptibility. The importance of other spatial factors varies depending on the landslide type. Based on the statistical results several landslide susceptibility models were developed using the Analytical Hierarchy Process method. These models gave very different results, with a prediction error ranging from 4.3% to 73%. As a final result of the research, the weights of important spatial factors from the best models were derived with the AHP method. Using probability measures, potentially hazardous areas were located in relation to population and road distribution, and hazard classes were assessed.     

Comparison of landslide susceptibility based on a decision-tree model and actual landslide occurrence: The Akaishi Mountains, Japan

This paper proposes a statistical decision-tree model to analyze landslide susceptibility in a wide area of the Akaishi Mountains, Japan. The objectives of this study were to validate the decision-tree model by comparing landslide susceptibility and actual landslide occurrence, and to reveal the relationships among landslide occurrence, topography, and geology. Landslide susceptibility was examined through ensemble learning with a decision tree. Decision trees are advantageous in that estimation processes and order of important explanatory variables are explicitly represented by the tree structures. Topographic characteristics (elevation, slope angle, profile curvature, plan curvature, and dissection and undissection height) and geological data were used as the explanatory variables. These topographic characteristics were calculated from digital elevation models (DEMs). The objective variables were landslide occurrence and reactivation data between 1992 and 2002 that were depicted by satellite image analysis. Landslide susceptibility was validated by comparing actual data on landslides that occurred and reactivated after the model was constructed (between 2002 and 2004).This study revealed that, from 2002 to 2004, landslides tended to occur and reactivate in catchments with high landslide susceptibility. The landslide susceptibility map thus depicts the actual landslide occurrence and reactivation in the Akaishi Mountains. This result indicates that the decision-tree model has appropriate accuracy for estimating the probabilities of future landslides. The tree structure indicates that landslides occurred and reactivated frequently in the catchments that had an average slope angle exceeding ca. 29° and a mode of slope angle exceeding 33°, which agree well with previous studies. A decision tree also quantitatively expresses important explanatory variables at the higher order of the tree structure.     

Landslide susceptibility mapping using geological data, a DEM from ASTER images and an Artificial Neural Network (ANN)

An efficient and accurate method of generating landslide susceptibility maps is very important to mitigate the loss of properties and lives caused by this type of geological hazard. This study focuses on the development of an accurate and efficient method of data integration, processing and generation of a landslide susceptibility map using an ANN and data from ASTER images. The method contains two major phases. The first phase is the data integration and analysis, and the second is the Artificial Neural Network training and mapping. The data integration and analysis phase involve GIS based statistical analysis relating landslide occurrence to geological and DEM (digital elevation model) derived geomorphological parameters. The parameters include slope, aspect, elevation, geology, density of geological boundaries and distance to the boundaries. This phase determines the geological and geomorphological factors that are significantly correlated with landslide occurrence. The second phase further relates the landslide susceptibility index to the important geological and geomorphological parameters identified in the first phase through ANN training. The trained ANN is then used to generate a landslide susceptibility map. Landslide data from the 2004 Niigata earthquake and a DEM derived from ASTER images were used. The area provided enough landslide data to check the efficiency and accuracy of the developed method. Based on the initial results of the experiment, the developed method is more than 90% accurate in determining the probability of landslide occurrence in a particular area.     

Qualitative landslide susceptibility assessment by multicriteria analysis: A case study from San Antonio del Sur, Guantánamo, Cuba

Geomorphological information can be combined with decision-support tools to assess landslide hazard and risk. A heuristic model was applied to a rural municipality in eastern Cuba. The study is based on a terrain mapping units (TMU) map, generated at 1:50,000 scale by interpretation of aerial photos, satellite images and field data. Information describing 603 terrain units was collected in a database. Landslide areas were mapped in detail to classify the different failure types and parts. Three major landslide regions are recognized in the study area: coastal hills with rockfalls, shallow debris flows and old rotational rockslides denudational slopes in limestone, with very large deep-seated rockslides related to tectonic activity and the Sierra de Caujerí scarp, with large rockslides. The Caujerí scarp presents the highest hazard, with recent landslides and various signs of active processes. The different landforms and the causative factors for landslides were analyzed and used to develop the heuristic model. The model is based on weights assigned by expert judgment and organized in a number of components such as slope angle, internal relief, slope shape, geological formation, active faults, distance to drainage, distance to springs, geomorphological subunits and existing landslide zones. From these variables a hierarchical heuristic model was applied in which three levels of weights were designed for classes, variables, and criteria. The model combines all weights into a single hazard value for each pixel of the landslide hazard map. The hazard map was then divided by two scales, one with three classes for disaster managers and one with 10 detailed hazard classes for technical staff. The range of weight values and the number of existing landslides is registered for each class. The resulting increasing landslide density with higher hazard classes indicates that the output map is reliable. The landslide hazard map was used in combination with existing information on buildings and infrastructure to prepare a qualitative risk map. The complete lack of historical landslide information and geotechnical data precludes the development of quantitative deterministic or probabilistic models.     

GIS techniques for regional-scale landslide susceptibility assessment: the Sicily (Italy) case study

This study describes the assessment of landslide susceptibility in Sicily (Italy) at a 1:100,000 scale using a multivariate logistic regression model. The model was implemented in a GIS environment by using the ArcSDM (Arc Spatial Data Modeller) module, modified to develop spatial prediction through regional data sets. A newly developed algorithm was used to automatically extract the detachment area from mapped landslide polygons. The following factors were selected as independent variables of the logistic regression model: slope gradient, lithology, land cover, a curve number derived index and a pluviometric anomaly index. The above-described configuration has been verified to be the best one among others employing from three to eight factors. All the regression coefficients and parameters were calculated using selected landslide training data sets. The results of the analysis were validated using an independent landslide data set. On an average, 82% of the area affected by instability and 79% of the not affected area were correctly classified by the model, which proved to be a useful tool for planners and decision-makers.     

Analysis of LiDAR-derived topographic information for characterizing and differentiating landslide morphology and activity

This study used airborne laser altimetry (LiDAR) to examine the surface morphology of two canyon-rim landslides in southern Idaho. The high resolution topographic data were used to calculate surface roughness, slope, semivariance, and fractal dimension. These data were combined with historical movement data (Global Positioning Systems (GPS) and laser theodolite) and field observations for the currently active landslide, and the results suggest that topographic elements are related to the material types and the type of local motion of the landslide. Weak, unconsolidated materials comprising the toe of the slide, which were heavily fractured and locally thrust upward, had relatively high surface roughness, high fractal dimension, and high vertical and lateral movement. The body of the slide, which predominantly moved laterally and consists mainly of undisturbed, older canyon floor materials, had relatively lower surface roughness than the toe. The upper block, consisting of a down-dropped section of the canyon rim that has remained largely intact, had a low surface roughness on its upper surface and high surface roughness along fractures and on its west face (unrelated to landslide motion). The upper block also had a higher semivariance than the toe and body. The topographic data for a neighboring, older and larger landslide complex, which failed in 1937, are similarly used to understand surface morphology, as well as to compare to the morphology of the active landslide and to understand scale-dependent processes. The morphometric analyses demonstrate that the active landslide has a similar failure mechanism and is topographically more variable than the 1937 landslide, especially at scales > 20 m. Weathering and the larger scale processes of the 1937 slide are hypothesized to cause the lower semivariance values of the 1937 slide. At smaller scales (< 10 m) the topographic components of the two landslides have similar roughness and semivariance. Results demonstrate that high resolution topographic data have the potential to differentiate morphological components within a landslide and provide insight into the material type and activity of the slide. The analyses and results in this study would not have been possible with coarser scale digital elevation models (10-m DEM). This methodology is directly applicable to analyzing other geomorphic surfaces at appropriate scales, including glacial deposits and stream beds.     

GIS支持下三峡库区秭归县滑坡灾害空间预测

基于GIS空间分析和统计模型相结合进行区域评价与空间预测是滑坡灾害研究的重要方向之一。以三峡库区秭归县为研究区,选择坡度、坡向、边坡结构、工程岩组、排水系统、土地利用和公路开挖作为评价因子。为提高模型的预测精度、可信度和推广能力,利用窗口采样规则降低训练样本之间的空间相关性。建立Logistic回归模型,对滑坡灾害与评价因子进行定量相关性分析。计算研究区滑坡灾害易发性指数,对其进行聚类分析,绘制滑坡易发性分区图,其中高、中易发区占整个研究区面积的38.9%,主要分布在人类工程活动频繁和靠近排水系统的区域。经过验证,该模型的预测精度达到77.57%。     

中国地质灾害伤亡事件的空间格局及影响因素

对中国2000-2012年造成人员伤亡的地质灾害事件进行分析,其空间分布格局受地形等自然环境要素的影响,南多北少,主要位于川西山区和云贵高原地区,东南丘陵地区,北方黄土丘陵,以及祁连山脉和天山山脉等地区,但局部地区的分布格局表明其还受到人为因素影响。构建基于二元Logistic回归的中国地质灾害伤亡事件发生概率模型（CELC）,定量分析自然、人为因素的影响程度,结果表明GDP增长率是仅次于地形起伏度的第二大影响因素,GDP增长率每增加2.72%,地质灾害伤亡事件发生的概率变为原来的2.706倍。此外还有多年平均降水、植被覆盖度、岩性、土壤类型、断裂带、产业类型和人口密度等因素。将CELC模型应用于中国县域,计算各个县的地质灾害伤亡事件概率,发现尚未发生但概率较高的县有27个,或为贫困县、或为矿产工业县域,或为房产过度开发县,它们是未来中国需要重点防范地质灾害的县域。     

Probabilistic landslide hazard assessment at the basin scale

We propose a probabilistic model to determine landslide hazard at the basin scale. The model predicts where landslides will occur, how frequently they will occur, and how large they will be. We test the model in the Staffora River basin, in the northern Apennines, Italy. For the study area, we prepare a multi-temporal inventory map through the interpretation of multiple sets of aerial photographs taken between 1955 and 1999. We partition the basin into 2243 geo-morpho-hydrological units, and obtain the probability of spatial occurrence of landslides by discriminant analysis of thematic variables, including morphological, lithological, structural and land use. For each mapping unit, we obtain the landslide recurrence by dividing the total number of landslide events inventoried in the unit by the time span of the investigated period. Assuming that landslide recurrence will remain the same in the future, and adopting a Poisson probability model, we determine the exceedance probability of having one or more landslides in each mapping unit, for different periods. We obtain the probability of landslide size by analysing the frequency–area statistics of landslides, obtained from the multi-temporal inventory map. Assuming independence, we obtain a quantitative estimate of landslide hazard for each mapping unit as the joint probability of landslide size, of landslide temporal occurrence and of landslide spatial occurrence.     

Spatial pattern and influencing factors of landslide casualty events

Analysis of casualties due to landslides from 2000 to 2012 revealed that their spatial pattern was affected by terrain and other natural environmental factors, which resulted in a higher distribution of landslide casualty events in southern China than in northern China. Hotspots of landslide-generated casualties were in the western Sichuan mountainous area and Yunnan-Guizhou Plateau region, southeast hilly area, northern part of the loess hilly area, and Tianshan and Qilian Mountains. However, local distribution patterns indicated that landslide casualty events were also influenced by economic activity factors. To quantitatively analyse the influence of natural environment and human-economic activity factors, the Probability Model for Landslide Casualty Events in China (LCEC) was built based on logistic regression analysis. The results showed that relative relief, GDP growth rate, mean annual precipitation, fault zones, and population density were positively correlated with casualties caused by landslides. Notably, GDP growth rate ranked only second to relative relief as the primary factors in the probability of casualties due to landslides. The occurrence probability of a landslide casualty event increased 2.706 times with a GDP growth rate increase of 2.72%. In contrast, vegetation coverage was negatively correlated with casualties caused by landslides. The LCEC model was then applied to calculate the occurrence probability of landslide casualty events for each county in China. The results showed that there are 27 counties with high occurrence probability but zero casualty events. The 27 counties were divided into three categories: poverty-stricken counties, mineral-rich counties, and real-estate overexploited counties; these are key areas that should be emphasized in reducing landslide risk.     

中国地质灾害伤亡事件的空间格局及影响因素（英文）

Analysis of casualties due to landslides from 2000 to 2012 revealed that their spatial pattern was affected by terrain and other natural environmental factors, which resulted in a higher distribution of landslide casualty events in southern China than in northern China. Hotspots of landslide-generated casualties were in the western Sichuan mountainous area and Yunnan-Guizhou Plateau region, southeast hilly area, northern part of the loess hilly area, and Tianshan and Qilian Mountains. However, local distribution patterns indicated that landslide casualty events were also influenced by economic activity factors. To quantitatively analyse the influence of natural environment and human-economic activity factors, the Probability Model for Landslide Casualty Events in China(LCEC) was built based on logistic regression analysis. The results showed that relative relief, GDP growth rate, mean annual precipitation, fault zones, and population density were positively correlated with casualties caused by landslides. Notably, GDP growth rate ranked only second to relative relief as the primary factors in the probability of casualties due to landslides. The occurrence probability of a landslide casualty event increased 2.706 times with a GDP growth rate increase of 2.72%. In contrast, vegetation coverage was negatively correlated with casualties caused by landslides. The LCEC model was then applied to calculate the occurrence probability of landslide casualty events for each county in China. The results showed that there are 27 counties with high occurrence probability but zero casualty events. The 27 counties were divided into three categories: poverty-stricken counties, mineral-rich counties, and real-estate overexploited counties; these are key areas that should be emphasized in reducing landslide risk.     

Terrain Pattern Recognition and Spatial Decision Making for Regional Slope Stability Studies

A terrain partition scheme is presented that allows the identification of regions with high landslide risk in natural terrain zones on the basis of geomorphometric criteria from moderate resolution DEMs. The key factor being the terrain segmentation to aspect regions (regions formed by points preserving the same aspect direction) instead of using an artificial regular-grid terrain partition scheme. The study area is in western Greece (NW Peloponnesus) whereas a moderate resolution digital elevation model with spacing 75 m is used. Landslide inventory analysis and knowledge conceptualization identified that the landslide susceptibility of a particular aspect region is high, if the mean elevation is low and the mean gradient is high. Each aspect region was parametrically represented on the basis of its mean gradient and elevation. The domain of each parameter was divided to seven slices (classes) on the basis of the observed density. Subsequent knowledge based mapping identified aspect regions with high landslide susceptibility for the following spatial rule: (a) “mean slope in class 6 or 7” and (b) “mean elevation in class 1 to 5”. Alternatively the rule is expressed as mean slope to be equal or greater than 15^∘ whereas mean elevation to be in the range 0 to 750 m. These identified zones correspond to regions where historical landslides occurred (populated coastal areas in the North) as well as to south regions (natural terrain zone) where no landslide record is available, because of the limitations posed by the natural terrain landslide mapping program in Greece. The presented terrain segmentation technique combined to the spatial decision-making process, provided both an object framework for integrating geomorphometric parameters and a method for landslide risk analysis in natural terrain zones.     

A comparison of landslide rates following helicopter and conventional cable-based clear-cut logging operations in the Southwest Coast Mountains of British Columbia

A comparison of landslide rates following helicopter and conventional, cable-based, clear-cut logging was carried out using results from two independent terrain attribute studies in the Eldred and Lois River watersheds in the Southwest Coast Mountains of British Columbia. Landslides initiating from directly within a road prism were excluded from the study in order to focus the comparison on landslides related primarily to conventional versus helicopter yarding methods. A landslide rate of 0.02 landslides/ha was observed in 162 terrain polygons logged by helicopter 8 years prior to this study. Landslide rates in 38 gullied polygons were 0.06 landslides/ha. No landslides were observed in 124 open-slope polygons. Over a similar 8-year average period, 0.03 landslides/ha were observed in 142 cable-yarded terrain polygons; 0.06 and 0.02 landslides/ha occurred in gullied and open-slope polygons, respectively. t-Tests indicate that total landslide rates are not significantly different following helicopter and conventional logging; however, a dichotomy exists between gullied and open-slope terrain polygons. Landslide rates are not significantly different in gullied terrain but are significantly higher on open-slopes following conventional cable logging. Consequently, landslides appear to have a greater potential to occur in open-slope terrain following conventional logging, but differences in gullied polygons are less likely. Increased post-logging landslide rates in conventionally logged, open slopes are more likely the result of undetected road-related drainage changes than differences between helicopter and conventional yarding-related ground disturbance.     

An integrated TIN and Grid method for constructing multi‐resolution digital terrain models

Multi‐resolution terrain models are an efficient approach to improve the speed of three‐dimensional (3D) visualizations, especially for terrain visualization in Geographical Information Systems (GIS). As a further development to existing algorithms and models, a new model is proposed for the construction of multi‐resolution terrain models in a 3D GIS. The new model represents multi‐resolution terrains using two major methods for terrain representation: Triangulated Irregular Network (TIN) and regular grid (Grid). In this paper, first, the concepts and formal definitions of the new model are presented. Second, the methodology for constructing multi‐resolution terrain models based on the new model is proposed. Third, the error of multi‐resolution terrain models is analysed, and a set of rules is proposed to retain the important features (e.g. boundaries of man‐made objects) within the multi‐resolution terrain models. Finally, several experiments are undertaken to test the performance of the new model. The experimental results demonstrate that the new model can be applied to construct multi‐resolution terrain models with good performance in terms of time cost and maintenance of the important features. Furthermore, a comparison with previous algorithms/models shows that the speed of rendering for 3D walking/flying through has been greatly improved by applying the new model.     

Landslide susceptibility mapping by correlation between topography and geological structure: the Janghung area, Korea

The purpose of this study is to develop and apply the technique for landslide susceptibility analysis using geological structure in a Geographic Information System (GIS). In the study area, the Janghung area of Korea, landslide locations were detected from Indian Remote Sensing (IRS) satellite images by change detection, where the geological structure of foliation was surveyed and analysed. The landslide occurrence factors (location of landslide, geological structure and topography) were constructed into a spatial database. Then, strike and dip of the foliation and the aspect and slope of the topography were compared and the results, which were verified using landslide location data, show that foliation of gneiss has a geometrical relation to the joint or fault that leads to a landslide. Using the geometrical relations, the landslide susceptibility was assessed and verified. The verification results showed satisfactory agreement between the susceptibility map and the landslide location data.