唐山市土地利用防灾适宜度评价方法分析

基于GIS的防灾适宜度多准则评价(MCE)是土地利用防灾规划的关键。根据唐山市地质灾害资料,建立基于距离的防灾适宜度评价准则并量化评价指标。依据决策风险指标计算次序权重,应用层次分析程序(AHP)构建比较矩阵并计算准则权重,分析基于GIS的OWA方法、布尔决策和权重线性叠加(WLC)等多准则评价方法的决策风险,确定唐山市土地利用防灾评价策略。基于决策风险和指标补偿原则计算次序权重、准则权重和一致性比率,得到唐山市土地利用防灾适宜度评价结果,据此提出唐山市土地资源合理利用建议。     

Spatially and temporally distributed modeling of landslide susceptibility

Mapping of landslide susceptibility in forested watersheds is important for management decisions. In forested watersheds, especially in mountainous areas, the spatial distribution of relevant parameters for landslide prediction is often unavailable. This paper presents a GIS-based modeling approach that includes representation of the uncertainty and variability inherent in parameters. In this approach, grid-based tools are used to integrate the Soil Moisture Routing (SMR) model and infinite slope model with probabilistic analysis. The SMR model is a daily water balance model that simulates the hydrology of forested watersheds by combining climate data, a digital elevation model, soil, and land use data. The infinite slope model is used for slope stability analysis and determining the factor of safety for a slope. Monte Carlo simulation is used to incorporate the variability of input parameters and account for uncertainties associated with the evaluation of landslide susceptibility. This integrated approach of dynamic slope stability analysis was applied to the 72-km² Pete King watershed located in the Clearwater National Forest in north-central Idaho, USA, where landslides have occurred. A 30-year simulation was performed beginning with the existing vegetation covers that represented the watershed during the landslide year. Comparison of the GIS-based approach with existing models (FSmet and SHALSTAB) showed better precision of landslides based on the ratio of correctly identified landslides to susceptible areas. Analysis of landslide susceptibility showed that (1) the proportion of susceptible and non-susceptible cells changes spatially and temporally, (2) changed cells were a function of effective precipitation and soil storage amount, and (3) cell stability increased over time especially for clear-cut areas as root strength increased and vegetation transitioned to regenerated forest. Our modeling results showed that landslide susceptibility is strongly influenced by natural processes and human activities in space and time; while results from simulated outputs show the potential for decision-making in effective forest planning by using various management scenarios and controlling factors that influence landslide susceptibility. Such a process-based tool could be used to deal with real-dynamic systems to help decision-makers to answer complex landslide susceptibility questions.     

Multicriteria Geologic Data Analysis for Mineral Favorability Mapping: Application to a Metal Sulphide Mineralized Area,Ribeira Valley Metallogenic Province,Brazil

GIS techniques have been used in the evaluation of favorability for base-metal mineralization in an area comprising the Cerro Azul and Apiaí quadrangles (SG.22-X-B-IV and V, scale 1:100.000), Ribeira Valley, São Paulo and Paraná States, Brazil. Methods have been employed for selection and weighting of prospective variables when applying GIS techniques to a digital database consisting of geological, geochemical and airborne geophysics, and mineral occurrence information. The exploration variable selection and analysis were based on two mineralization models: (1) Panelas type, vein-type carbonate hosted, and (2) Perau type, sedimentary-exhalative. The overlay was performed by weighted linear combination (WLC) and order weighted average (OWA) methods. Both methods proved suitable for the study area, yielding similar results. The ordered weighted averaging analysis provided the best results, with favorability maps showing a large number of classes occupying relatively minor areas. In comparison, the weighted linear combination analysis produced more coherent results but without details for minor areas. The prospective parameters obtained are considered suitable for both Perau and Panelas types. Both methods are inexpensive, and are suitable for selection of prospective areas during geological surveys in areas similar to the studied one.     

Exploring multi-criteria decision strategies in GIS with linguistic quantifiers: an extension of the analytical network process using ordered weighted averaging operators

Many real-world spatial planning and management problems give rise to a geographical information system (GIS)-based multi-criteria decision-making. Analytical network process (ANP) provides a comprehensive methodology for representing complex multi-criteria decision-making problems as a network of criteria and alternatives, where feedback and interdependence relationships may exist within and between all the criteria and alternatives. Experts’ experiences are used to estimate relative magnitudes of tangible and intangible factors through paired comparisons in order to make rational and consistent decisions. However, the GIS-based ANP, an adoption of weighted linear aggregation rule, typically employed a high trade-off decision strategy and neglected other decision strategies. This paper develops a novel GIS-based multi-criteria evaluation (MCE) procedure by extending the ANP using fuzzy quantifiers-guided ordered weighted averaging (OWA) operators. This extension, which generalizes the aggregation process used in the ANP, would provide a generic powerful decision-making tool that allows decision-makers to define a decision strategy on a continuum between pessimistic (risk-averse) and optimistic (risk-taking) strategies. By changing the linguistic quantifiers, the GIS-based ANP–OWA can generate a wide range of decision strategies taking into accounts the level of risk the decision-makers wish to assume in their MCE. A land-use suitability analysis in a region of Saudi Arabia is presented to demonstrate the application of the proposed procedure.     

Uncertainty due to DEM error in landslide susceptibility mapping

Terrain attributes such as slope gradient and slope shape, computed from a gridded digital elevation model (DEM), are important input data for landslide susceptibility mapping. Errors in DEM can cause uncertainty in terrain attributes and thus influence landslide susceptibility mapping. Monte Carlo simulations have been used in this article to compare uncertainties due to DEM error in two representative landslide susceptibility mapping approaches: a recently developed expert knowledge and fuzzy logic-based approach to landslide susceptibility mapping (efLandslides), and a logistic regression approach that is representative of multivariate statistical approaches to landslide susceptibility mapping. The study area is located in the middle and upper reaches of the Yangtze River, China, and includes two adjacent areas with similar environmental conditions – one for efLandslides model development (approximately 250 km²) and the other for model extrapolation (approximately 4600 km²). Sequential Gaussian simulation was used to simulate DEM error fields at 25-m resolution with different magnitudes and spatial autocorrelation levels. Nine sets of simulations were generated. Each set included 100 realizations derived from a DEM error field specified by possible combinations of three standard deviation values (1, 7.5, and 15 m) for error magnitude and three range values (0, 60, and 120 m) for spatial autocorrelation. The overall uncertainties of both efLandslides and the logistic regression approach attributable to each model-simulated DEM error were evaluated based on a map of standard deviations of landslide susceptibility realizations. The uncertainty assessment showed that the overall uncertainty in efLandslides was less sensitive to DEM error than that in the logistic regression approach and that the overall uncertainties in both efLandslides and the logistic regression approach for the model-extrapolation area were generally lower than in the model-development area used in this study. Boxplots were produced by associating an independent validation set of 205 observed landslides in the model-extrapolation area with the resulting landslide susceptibility realizations. These boxplots showed that for all simulations, efLandslides produced more reasonable results than logistic regression.     

城市居住用地防灾适宜度评价的OWA方法

城市居住用地的安全受到多种灾害威胁,防灾决策的主要问题之一是开发基于GIS的多准则评价 (MCE)方法。文章介绍了在基于GIS的次序权重平均法 (OWA)中,依据重要性等级计算次序权重和应用层次分析程序 (AHP)构建比较矩阵计算准则权重的方法;通过OWA评价方法与布尔决策和权重线性叠加 (WLC)等多准则评价方法的比较,分析了不同方法的决策风险和指标等级重要次序的影响,并以唐山市为例,建立了各灾害因子的适宜度评价准则,计算了唐山市居住用地的防灾适宜度,对唐山市居住用地的合理开发提出了几点建议。     

Comparison of satellite and air photo based landslide susceptibility maps

Landslide susceptibility maps can be prepared in a variety of ways. Many geoscientists favour the use of an overlay model approach in which several map layers are combined by some arithmetic rules to determine the potential for sliding in an area or region. The resulting susceptibility maps, although based on a subjective weighting of relevant factors, can often be of high accuracy and utility. In order to obtain the relevant input data for this type of analysis, remotely sensed data are often used. To date, susceptibility mapping, just as the mapping of historic and individual landslides, has tended to require higher-resolution imagery. This has somewhat limited the application of landslide susceptibility mapping. While high-resolution air photo or satellite imagery is superior to lower resolution imagery for the purpose of mapping of historic and individual landslides, such higher levels of resolution may not be required for the development of landslide susceptibility maps. In order to determine if medium-resolution satellite imagery, such as SPOT or ASTER, could provide the needed data for landslide susceptibility mapping, a comparison was undertaken of landslide susceptibility model output resulting from the use of stereo NAPP aerial photography versus the use of data obtained from stereo SPOT imagery. The test area selected for this study consisted of two watersheds, Pena Canyon and Big Rock Canyon, situated west of Santa Monica, California, USA, along the Pacific Coast Highway. Both watersheds have a long and well-documented history of landslide activity and sufficient geologic variability and complexity to provide a good test site. The specific overlay model used in this evaluation required input data consistent with the needs of many other models of this type. The model output derived from the two different data sources and presented here in the form of susceptibility maps were virtually identical. Statistical and difference analysis confirmed that both methods of obtaining input data provide similar results and successfully identified landslide prone areas. These results suggest that satellite imagery, in this instance, SPOT images, could potentially be used in lieu of conventional air photos, to evaluate landslide susceptibility. In many situations, especially in the case of remote locations and/or developing countries, this capability should result in substantial savings in terms of time, financial resources, and overall viability.     

The impact of spatial scale: exploring urban butterfly abundance and richness patterns using multi-criteria decision analysis and principal component analysis

ABSTRACT

Modelling changes in biodiversity have become a necessary component of smart urban planning practices. However, concepts such as biodiversity are often evaluated using area-based composite indices, the results of which are heavily reliant on specific parameters chosen. This paper explores the design and implementation of a butterfly biodiversity index by comparing two widely accepted modelling techniques: principal component analysis and spatial multi-criteria decision analysis (MCDA). A high degree of scale dependency has been demonstrated in previous studies exploring the use of area-based composite measures. To evaluate the impact of scale, each model was assessed at two different spatial resolutions. The outcomes were analyzed, mapped and compared using ordinary least squares, geographically weighted regression and global Moran’s I to evaluate relative biodiversity patterns across the City of Toronto, Canada. Findings indicate that the impact of spatial scale was significant, whereby the coarser resolution models were found to be more highly correlated with biodiversity, compared to the finer resolution models. The results of this study contribute to a growing body of literature that explores key conceptual questions regarding the robustness of GIS-based MCDA, the impact of scale in urban ecology studies, and the use of composite indices to manage spatial ecological data.     

Risk-Based Analysis in Mineral Potential Mapping: Application of Quantifier-Guided Ordered Weighted Averaging Method

Natural Resources Research - In this work, a quantifier-guided ordered weighted averaging (OWA) method was employed for mineral potential mapping (MPM) in Nowchun Cu–Mo prospect, SE Iran. The...     

DEM误差对滑坡危险性评价模型的影响

基于数字高程模型(DEM)计算得到的坡度、坡向等地形属性是滑坡危险性评价模型的重要输入数据, DEM误差会导致地形属性计算结果不确定性, 进而影响滑坡危险性评价模型的结果。本文选择基于专家知识的滑坡危险性评价模型和逻辑斯第回归模型, 采用蒙特卡洛模拟方法, 研究DEM误差所导致的滑坡危险性评价模型结果不确定性。研究区位于长江中上游的重庆开县, 采用5 m分辨率的DEM, 以序贯高斯模拟方法模拟了不同大小(误差标准差为1 m、7.5 m、15 m)和空间自相关性(变程为0 m、30 m、60 m、120 m)的12 类DEM误差场参与滑坡危险性评价。每次模拟包括100 个实现, 通过对每次模拟分别计算滑坡危险性评价结果的标准差图层和分类一致性百分比图层, 用以评价结果不确定性。评价结果表明, 在不同的DEM精度下, 两个滑坡危险性评价模型所得结果的总体不确定性随空间自相关程度的变化趋势并不相同。当DEM空间自相关性程度不同时, 基于专家知识的滑坡危险性评价模型的评价结果总体不确定随着DEM误差增加而呈现不同的变化趋势, 而逻辑斯第回归模型的评价结果总体不确定性随着DEM误差大小增加而单调增加。从评价结果总体不确定性角度而言, 总体上逻辑斯第回归模型比基于专家知识的滑坡危险性评价模型更加依赖于DEM数据质量。     

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.     

一种基于地理环境相似度的滑坡负样本可信度度量方法

滑坡负样本在统计型滑坡危险度制图中具有重要作用,能抑制统计模型对滑坡危险度的高估。当前滑坡负样本采样方法采集的负样本可信度未知,在负样本采样过程中,极有可能将那些潜在滑坡点错选为负样本,这些假的负样本会降低负样本集的质量和训练样本集的质量,进而影响统计模型的精度。本文基于“地理环境越相似、地理特征越相似”的地理学常识,认为与正样本有着相似地理环境的点极有可能是未来发生滑坡的点;与正样本的地理环境越不相似的点,则越有可能是负样本。基于此假设提出一种基于地理环境相似度的负样本可信度度量方法,将该方法应用于滑坡灾害频发的陇南山区油房沟流域,对油房沟进行滑坡负样本可信度评价制图;使用油房沟流域的滑坡发生初始面来验证该方法的有效性。结果发现：滑坡发生初始面上所有栅格点的负样本可信度平均值为0.26,超过95%的栅格点的负样本可信度都小于0.5,说明本文提出的负样本可信度度量方法合理。     

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.     

Web-based geospatial multiple criteria decision analysis using open software and standards

The emerging ubiquity of geospatial information is providing an unprecedented opportunity to apply Geographical Information Systems (GIS)-based multiple criteria decision analysis (MCDA) to a broad spectrum of use cases. Volunteered geographic information, open GIS software, geoservice-based tools, cloud-based virtualized platforms, and worldwide collaboration of both domain experts and general users have greatly increased the quantity and accessibility of geospatially referenced data resources. Currently, there is a lack of GIS-based MCDA tools that integrate this decision-driven process within a widely accessible, robust geoframework environment, designed for user-friendly interaction. In this contribution, we present a conceptual workflow and proof-of-concept software application, Geocentric Environment for Analysis and Reasoning (GEAR), which provides a viable transition path to enhance geospatial MCDA in the age of open GIS. We propose a Web-based platform that leverages open-source geotechnologies to incorporate a wide variety of geospatial data formats in a common solution space to allow for spatially enhanced and time-relevant decision analysis. Through the proposed workflow, a user can ingest and modify heterogeneous data formats, exploit temporally tagged data sources, create multicriteria decision analysis models, and visualize the results in an iterative and collaborative workspace. A sample case study applied to disaster relief is used to demonstrate the prototype and workflow. This proof-of-concept Web-based application provides a notional pathway of how to connect open-source data to open-source analysis through a geospatially enabled MCDA workflow that could be virtually accessible to many levels of decision makers from individuals to entire organizations.     

Landslide susceptibility mapping based on Support Vector Machine: A case study on natural slopes of Hong Kong, China

The Support Vector Machine (SVM) is an increasingly popular learning procedure based on statistical learning theory, and involves a training phase in which the model is trained by a training dataset of associated input and target output values. The trained model is then used to evaluate a separate set of testing data. There are two main ideas underlying the SVM for discriminant-type problems. The first is an optimum linear separating hyperplane that separates the data patterns. The second is the use of kernel functions to convert the original non-linear data patterns into the format that is linearly separable in a high-dimensional feature space. In this paper, an overview of the SVM, both one-class and two-class SVM methods, is first presented followed by its use in landslide susceptibility mapping. A study area was selected from the natural terrain of Hong Kong, and slope angle, slope aspect, elevation, profile curvature of slope, lithology, vegetation cover and topographic wetness index (TWI) were used as environmental parameters which influence the occurrence of landslides. One-class and two-class SVM models were trained and then used to map landslide susceptibility respectively. The resulting susceptibility maps obtained by the methods were compared to that obtained by the logistic regression (LR) method. It is concluded that two-class SVM possesses better prediction efficiency than logistic regression and one-class SVM. However, one-class SVM, which only requires failed cases, has an advantage over the other two methods as only “failed” case information is usually available in landslide susceptibility mapping.     

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

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

Automated landslide mapping using spectral analysis and high-resolution topographic data: Puget Sound lowlands, Washington, and Portland Hills, Oregon

Landslide inventory maps are necessary for assessing landslide hazards and addressing the role slope stability plays in landscape evolution over geologic timescales. However, landslide inventory maps produced with traditional methods — aerial photograph interpretation, topographic map analysis, and field inspection — are often subjective and incomplete. The increasing availability of high-resolution topographic data acquired via airborne Light Detection and Ranging (LiDAR) over broad swaths of terrain invites new, automated landslide mapping procedures. We present two methods of spectral analysis that utilize LiDAR-derived digital elevation models of the Puget Sound lowlands, Washington, and the Tualatin Mountains, Oregon, to quantify and automatically map the topographic signatures of deep-seated landslides. Power spectra produced using the two-dimensional discrete Fourier transform and the two-dimensional continuous wavelet transform identify the characteristic spatial frequencies of deep-seated landslide morphologic features such as hummocky topography, scarps, and displaced blocks of material. Spatial patterns in the amount of spectral power concentrated in these characteristic frequency bands highlight past slope instabilities and allow the delineation of landslide terrain. When calibrated by comparison with detailed, independently compiled landslide inventory maps, our algorithms correctly classify an average of 82% of the terrain in our five study areas. Spectral analysis also allows the creation of dominant wavelength maps, which prove useful in analyzing meter-scale topographic expressions of landslide mechanics, past landslide activity, and landslide-modifying geomorphic processes. These results suggest that our automated landslide mapping methods can create accurate landslide maps and serve as effective, objective, and efficient tools for digital terrain analysis.     

Evaluation and comparison of landslide susceptibility mapping methods: a case study for the Ulus district,Bartın,northern Turkey

The purpose of this study was to investigate the capabilities of different landslide susceptibility methods by comparing their results statistically and spatially to select the best method that portrays the susceptibility zones for the Ulus district of the Bart?n province (northern Turkey). Susceptibility maps based on spatial regression (SR), linear discriminant analysis (LDA), quadratic discriminant analysis (QDA), logistic regression (LR) method, and artificial neural network method (ANN) were generated, and the effect of each geomorphological parameter was determined. The landslide inventory map digitized from previous studies was used as a base map for landslide occurrence. All of the analyses were implemented with respect to landslides classified as rotational, active, and deeper than 5 m. Three different sets of data were used to produce nine explanatory variables (layers). The study area was divided into grids of 90 m × 90 m, and the ‘seed cell’ technique was applied to obtain statistically balanced population distribution over landslide inventory area. The constructed dataset was divided into two datasets as training and test. The initial assessment consisted of multicollinearity of explanatory variables. Empirical information entropy analysis was implemented to quantify the spatial distribution of the outcomes of these methods. Results of the analyses were validated by using success rate curve (SRC) and prediction rate curve (PRC) methods. Additionally, statistical and spatial comparisons of the results were performed to determine the most suitable susceptibility zonation method in this large-scale study area. In accordance with all these comparisons, it is concluded that ANN was the best method to represent landslide susceptibility throughout the study area with an acceptable processing time.     

Automatic selection of weights for GIS-based multicriteria decision analysis: site selection of transmission towers as a case study

Transmission line (TL) siting consists of finding suitable land to build transmission towers. This is just one of the numerous complex geographical problems often solved using GIS-based multicriteria decision analysis (MCDA), which is a set of techniques that weight several geographical features to identify suitable locations. This technique is mostly employed using expert knowledge to identify the correct set of weights; thus adding a certain amount of subjectivity to the analysis, meaning that for the same problem if we change the experts involved, we may reach different results.This research is a first attempt to try and solve this issue. We employed a statistical analysis to quantitatively calculate these weights and we tested our method on a case study about transmission line siting in Switzerland. We compared the distances between each sample in our dataset, in this case study these are location of transmission towers, with each geographical feature, e.g. distance from water features. Then we calculate the same distances but for random points, sampled throughout the study area. The reasoning behind this method is that if samples present a distance from a geographic feature statistically different from the random, it means that the feature played an important role in dictating the location of the sample. In this case for instance, high-voltage transmission towers are purposely built as far away as possible from urban areas. Random points are on the contrary by definition sampled without any constraint. Therefore, when comparing the two datasets, we should find that transmission towers have a larger average distance from urban areas than random points. This allows us to determine that this criterion (i.e. distance from urban centers) is important for planning new TL.The results indicate that this method can successfully weight and rank the most important criteria to be considered for an MCDA analysis, in line with weights proposed in the literature. The advantage of the proposed technique is that it completely excludes human factors, thus potentially increasing the social acceptance of the MCDA results.     

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.