Estimating potential landslide sites of an upland sub-watershed in Western Ghat’s of Kerala (India) through frequency ratio and GIS

The purpose of this study is to assess the susceptibility of landslides in parts of Western Ghats, Kerala, India, using a geographical information system (GIS). Landslide inventory of the area was made by detailed field surveys and the analysis of the topographical maps. The landslide triggering factors are considered to be slope angle, slope aspect, slope curvature, slope length, distance from drainage, distance from lineaments, lithology, land use and geomorphology. ArcGIS version 8.3 was used to manipulate and analyse all the collected data. Probabilistic-likelihood ratio was used to create a landslide susceptibility map for the study area. The result was validated using the Area under Curve (AUC) method and temporal data of landslide occurrences. The validation results showed satisfactory agreement between the susceptibility map and the existing data on landslide locations. As the result, the success rate of the model was (84.46%) and the prediction rate of the model was (82.38%) shows high prediction accuracy. In the reclassified final landslide susceptibility zone map, 5.68% of the total area is classified as critical in nature. The landslide susceptibility map thus produced can be used to reduce hazards associated with landslides and to land cover planning.     

The application of GIS-based bivariate statistical methods for landslide hazards assessment in the upper Tons river valley,Western Himalaya,India

Landslides are one of the most frequent and common natural hazards in many parts of Himalaya. To reduce the potential risk, the landslide susceptibility maps are one of the first and most important steps in the landslide hazard mitigation. Earth observation satellite and geographical information system-based techniques have been used to derive and analyse various geo-environmental parameters significant to landslide hazards. In this study, a bivariate statistics method was used for spatial modelling of landslide susceptibility zones. For this purpose, thematic layers including landslide inventory, geology, slope angle, slope aspect, geomorphology, slope morphology, drainage density, lineament and land use/land cover were used. A large number of landslide occurrences have been observed in the upper Tons river valley area of Western Himalaya. The result has been used to spatially classify the study area into zones of very high, high, moderate, low and very low landslide susceptibility zones. About 72% of active landslides have been observed to occur in very high and high hazard zones. The result of the analysis was verified using the landslide location data. The validation result shows significant agreement between the susceptibility map and landslide location. The result can be used to reduce landslide hazards by proper planning.     

Landslide susceptibility assessment in the İzmir (West Anatolia,Turkey) city center and its near vicinity by the logistic regression method

A landslide susceptibility assessment for İzmir city (Western Turkey), which is the third biggest city of Turkey, was performed by a logistic regression method. A database of landslide characteristics was prepared using detailed field surveys. The major landslides in the study area are generally observed in the field, dominated by weathered volcanics, and 39.63% of the total landslide area is in this unit. The parameters of lithology, slope gradient, slope aspect, distance to drainage, distance to roads and distance to fault lines were used as variables in the logistic regression analysis. The effect of each parameter on landslide occurrence was assessed from the corresponding coefficients that appear in the logistic regression function. On the basis of the obtained coefficients, lithology plays the most important role in determining landslide occurrence and distribution. Slope gradient has a more significant effect than the other geomorphological parameters, such as slope aspect and distance to drainage. Using a predicted map of probability, the study area was classified into five categories of landslide susceptibility: very low, low, moderate, high and very high. Whereas 49.65% of the total study area has very low susceptibility, very high susceptibility zones make up 11.69% of the area.     

Landslide susceptibility mapping of Tevankarai Ar sub-watershed,Kodaikkanal taluk,India, using weighted similar choice fuzzy model

This paper deals with the landslide susceptibility zonation of Tevankarai Ar sub-watershed using weighted similar choice fuzzy method in a GIS environment. There has been a rapid increase in landslide occurrences in the Kodaikkanal town and area surrounding the town specially in the settlements around the town and road links leading to and from the town. This necessitates a detailed study of slope instability problems in this area. It is observed that these incidences occur frequently during the monsoon and summer showers. Rainfall is identified as the prime triggering factor. Eleven physical factors that cause instability are identified as causative factors from the field investigations and landslide occurrences. Land use pattern, slope gradient, curvature and aspect, weathering index which are evaluated from the weathering ratios of different chemical constituents of the three major lithological variations, soil type, hydraulic conductivity of soil and soil thickness, geomorphology, drainage, and lineament have been utilized to prepare the spatial variation. A weighted similar choice fuzzy model which ranks a set of alternatives by identifying the similarity between the outcome of alternatives and outcome of ideal alternatives is used to rank the causative factors. Each causative factor is classified into sub-categories and rated based on their effect on stimulating the landslide event using qualitative judgment derived from field studies and landslide history. The prepared thematic maps of causative factors are integrated, utilizing the GIS software Arcmap. The outcome has projected the low, moderate, high, and very high landslide susceptibility zones. The high-hazard and very high-hazard areas fall in the northwestern part characterized by croplands and agricultural plantations, while the moderate hazard zones are seen in prominent settlements and low-hazard zones are observed in the sparse settlements and zones of less agricultural activity. The model is verified using the relative landslide density (R) index, and the susceptibility map is found to be consistent with the mapped landslide incidences. The results from this study illustrate that the use of weighted similar choice fuzzy method is suitable for landslide susceptibility mapping on regional scale in growing hill towns as Kodaikkanal town.     

Landslide temporal analysis and susceptibility assessment as bases for landslide mitigation, Machu Picchu, Peru

Multi-temporal landslide occurrence information acquired through aerial photo interpretation and field mapping was used to assess occurrence frequencies on the slopes around the UNESCO cultural world heritage site of Machu Picchu, Peru. This showed that the coarse time resolution of the historical landslide information may lead to inaccurate interpretations regarding landslide occurrence frequencies in some parts of the study area. In addition, the assumption that the past landslide frequency can be used to describe the future landslide occurrence was not proved in the study area. Thereafter, unique conditional analyses were undertaken to assess landslide susceptibility using a limited number of preparatory factor maps. It showed that large majority of the Inca City is located on least susceptible areas within the region. The results of the susceptibility assessment combined with landslide occurrence frequencies may serve as a basis for the landslide hazard mitigation in the studied area. For these purposes, pixel-based susceptibility maps were generalized into expert-defined landslide management units. These units provide site managers with easily understandable and applicable hence reliable information about future landslide occurrences. An approach describing usage of the resulting susceptibility maps for onsite mitigation purposes was described with respect to the needs of Machu Picchu site managers.     

GIS-based landslide susceptibility mapping using bivariate statistical analysis in Devrek (Zonguldak-Turkey)

Devrek town with increasing population is located in a hillslope area where some landslides exist. Therefore, landslide susceptibility map of the area is required. The purpose of this study was to generate a landslide susceptibility map using a bivariate statistical index and evaluate and compare the results of the statistical analysis conducted with three different approaches in seed cell concept resulting in different data sets in Geographical Information Systems (GIS) based landslide susceptibility mapping applied to the Devrek region. The data sets are created from the seed cells of (a) crowns and flanks, (b) only crowns, and (c) only flanks of the landslides by using ten different causative parameters of the study area. To increase the data dependency of the analysis, all parameter maps are classified into equal frequency classes based directly on the percentile divisions of each corresponding seed cell data set. The resultant maps of the landslide susceptibility analysis indicate that all data sets produce fairly acceptable results. In each data set analysis, elevation, lithology, slope, aspect, and drainage density parameters are found to be the most contributing factors in landslide occurrences. The results of the three data sets are compared using Seed Cell Area Indexes (SCAI). This comparison shows that the crown data set produces the most accurate and successful landslide susceptibility map of the study area.     

Landslides susceptibility mapping based on geographical information system,GuiZhou, south-west China

The purpose of this study is to assess the susceptibility of landslides around the area of Guizhou province, in south-west of China, using a geographical information system (GIS). The base map is prepared by visiting the field area and mapping individual landslide at a scale of 1:500,000 topographic maps. In the study, slope, lithology, landslide inventory, tectonic activity, drainage distribution and annual precipitation were taken as independent causal factors. Therefore, six causal factors maps are prepared by collecting information from various authorized sources and converting them in to GIS maps. The susceptibility assessment is based on the qualitative map combination model and trapezoidal fuzzy number weighting (TFNW) approach. Using a predicted map of probability, the study area was classified into four categories of landslide susceptibility: low, moderate, high and very high. In addition, the weighting procedure showed that the TFNW is an efficient method for landslide causal factors weighting.     

Support Vector Machines for Landslide Susceptibility Mapping: The Staffora River Basin Case Study, Italy

The aim of this study is the application of support vector machines (SVM) to landslide susceptibility mapping. SVM are a set of machine learning methods in which model capacity matches data complexity. The research is based on a conceptual framework targeted to apply and test all the procedural steps for landslide susceptibility modeling from model selection, to investigation of predictive variables, from empirical cross-validation of results, to analysis of predicted patterns. SVM were successfully applied and the final susceptibility map was interpreted via success and prediction rate curves and receiver operating characteristic (ROC) curves, to support the modeling results and assess the robustness of the model. SVM appeared to be very specific learners, able to discriminate between the informative input and random noise. About 78% of occurrences was identified within the 20% of the most susceptible study area for the cross-validation set. Then the final susceptibility map was compared with other maps, addressed by different statistical approaches, commonly used in susceptibility mapping, such as logistic regression, linear discriminant analysis, and naive Bayes classifier. The SVM procedure was found feasible and able to outperform other techniques in terms of accuracy and generalization capacity. The over-performance of SVM against the other techniques was around 18% for the cross-validation set, considering the 20% of the most susceptible area. Moreover, by analyzing receiver operating characteristic (ROC) curves, SVM appeared to be less prone to false positives than the other models. The study was applied in the Staffora river basin (Lombardy, Northern Italy), an area of about 275 km² characterized by a very high density of landslides, mainly superficial slope failures triggered by intense rainfall events.     

Macro-zonation of landslide susceptibility in Garamaura-Swarghat-Gambhar section of national highway 21, Bilaspur District, Himachal Pradesh (India)

Landslides the most common geo-hazard in hilly terrain are short lived phenomena but cause extraordinary landscape changes and destruction of life and property. The frequency and intensity of landslides occurrences along NH-21 during the rainy season not only disrupts traffic movement but also misbalance the agro-economic and developmental activities of the region frittering away thousand crores of rupees from the exchequer. An assessment of landslide susceptibility is, therefore, a prerequisite for sustainable development of the region. The present study deals with the preparation of macro-zonation maps of landslide susceptibility in an area of about 100 sq km on 1:50,000 scale across Garamaura-Swarghat section of National Highway-21. The map has been prepared by superimposing the terrain evaluation maps in a particular zone such as lithological map, structural map, slope morphometry map, relative relief map, land use and land cover map and hydrological condition map using landslide susceptibility evaluation factor rating scheme and calculating the total estimated susceptibility as per the guidelines of IS: 14496 (Part-2) 1998). Numerical weightages are assigned to the prime causative factors of slope instability such as lithology, structure, slope morphometery, relative relief, land use and groundwater conditions as per the scheme approved by Bureau of Indian Standard for the purpose of landslide susceptibility zonation. The area depicts zones of different instability. The identified susceptibility zones compared with landslide intensity in the area show some congruence with the weightages of the inputs. The incongruence in intensity and frequency of landslide occurrences and the inferred susceptibility zones of BIS scheme allow other geotechnical considerations and causative factors to be incorporated for the landslide susceptibility zonation.     

Improvement of statistical landslide susceptibility mapping by using spatial and global regression methods in the case of More and Romsdal (Norway)

Statistical models are one of the most preferred methods among many landslide susceptibility assessment methods. As landslide occurrences and influencing factors have spatial variations, global models like neural network or logistic regression (LR) ignore spatial dependence or autocorrelation characteristics of data between the observations in susceptibility assessment. However, to assess the probability of landslide within a specified period of time and within a given area, it is important to understand the spatial correlation between landslide occurrences and influencing factors. By including these relations, the predictive ability of the developed model increases. In this respect, spatial regression (SR) and geographically weighted regression (GWR) techniques, which consider spatial variability in the parameters, are proposed in this study for landslide hazard assessment to provide better realistic representations of landslide susceptibility. The proposed model was implemented to a case study area from More and Romsdal region of Norway. Topographic (morphometric) parameters (slope angle, slope aspect, curvature, plan, and profile curvatures), geological parameters (geological formations, tectonic uplift, and lineaments), land cover parameter (vegetation coverage), and triggering factor (precipitation) were considered as landslide influencing factors. These influencing factors together with past rock avalanche inventory in the study region were considered to obtain landslide susceptibility maps by using SR and LR models. The comparisons of susceptibility maps obtained from SR and LR show that SR models have higher predictive performance. In addition, the performances of SR and LR models at the local scale were investigated by finding the differences between GWR and SR and GWR and LR maps. These maps which can be named as comparison maps help to understand how the models estimate the coefficients at local scale. In this way, the regions where SR and LR models over or under estimate the landslide hazard potential were identified.     

基于GIS的白龙江流域甘肃段滑坡易发性评价

在甘肃省白龙江流域地质灾害资料收集及现场调查的基础上, 统计分析了该区滑坡发育与地层岩性、坡度、坡向、高程、断裂、植被等因素之间的关系, 建立了白龙江流域滑坡易发性评价指标体系。采用基于GIS的层次分析法评价模型, 完成了滑坡易发性分区评价, 将研究区滑坡按易发程度划分为高易发区、中易发区、低易发区和极低易发区, 其中, 高易发区占研究区总面积的13.59%, 主要分布在断裂带、白龙江两侧以及软弱岩土体分布的区域; 中易发区占27.85%;主要分布在白龙江支流以及主要道路两侧的一定范围内; 低易发区占33.09%, 主要分布在海拔相对较高、植被覆盖度较高、基本上无断裂带通过的区域; 其余区域为极低易发区, 占25.46%。对比分析显示评价结果与实际滑坡发育情况吻合, 可以较好地反映区内滑坡灾害发育的总体特征。      

Landslide susceptibility mapping based on GIS and information value model for the Chencang District of Baoji,China

The main objective of this study was to apply a statistical (information value) model using geographic information system (GIS) to the Chencang District of Baoji, China. Landslide locations within the study area were identified using reports and aerial photographs, and a field survey. A total of 120 landslides were mapped, of which 84 (70 %) were randomly selected for building the landslide susceptibility model. The remaining 36 (30 %) were used for model validation. We considered a total of 10 potential factors that predispose an area to a landslide for the landslide susceptibility mapping. These included slope degree, altitude, slope aspect, plan curvature, geomorphology, distance from faults, lithology, land use, mean annual rainfall, and peak ground acceleration. Following an analysis of these factors, a landslide susceptibility map was produced using the information value model with GIS. The resulting landslide susceptibility index was divided into five classes (very high, high, moderate, low, and very low) using the natural breaks method. The corresponding distribution area percentages were 29.22, 25.14, 15.66, 15.60, and 14.38 %, respectively. Finally, landslide locations were used to validate the results of the landslide susceptibility map using areas under the curve (AUC). The AUC plot showed that the susceptibility map had a success rate of 81.79 % and a prediction accuracy of 82.95 %. Based on the results of the AUC evaluation, the landslide susceptibility map produced using the information value model exhibited good performance.     

Landslide susceptibility mapping and factor effect analysis using frequency ratio in a catchment scale: a case study from Garuwa sub-basin,East Nepal

For assessing landslide susceptibility, the spatial distribution of landslides in the field is essential. The landslide inventory map is prepared on the basis of historical information of individual landslide events from different sources such as previously published reports, satellite imageries, aerial photographs and interview with local inhabitants. Then, the distribution of landslides in the study area is verified with field surveys. However, the selection of contributing factors for modelling landslide susceptibility is an inhibit task. The previous studies show that the factors are chosen as per availability of data. This paper documents the landslide susceptibility mapping in the Garuwa sub-basin, East Nepal using frequency ratio method. Nine different contributing factors are considered: slope aspect, slope angle, slope shape, relative relief, geology, distance from faults, land use, distance from drainage and annual rainfall. To analyse the effect of contributing factors, the landslide susceptibility index maps are generated four times using (a) topographical factors and geological factors, (b) topographical factors, geological factors and land use, (c) topographical factors, geological factors, land use and drainage and (d) all nine causative factors. By comparing with the pre-existing landslides, the fourth case (considering all nine causative factors) yields the best success rate accuracy, i.e. 81.19 %, which is then used to produce the final landslide susceptibility zonation map. Then, the final landslide susceptibility map is validated through chi-square test. The standard chi-square value with 3 degrees of freedom at the 0.001 significance level is 16.3, whereas the calculated chi-square value is 7,125.79. Since the calculated chi-square value is greater than the standard chi-square value, it can be concluded that the landslide susceptibility map is considered as statistically significant. Moreover, the results show that the predicted susceptibility levels are found to be in good agreement with the past landslide occurrences.     

Landslide Susceptibility Mapping Using Fuzzy Logic System and Its Influences on Mainlines in Lashgarak Region,Tehran, Iran

Landslide susceptibility mapping is among the useful tools applied in disaster management and planning development activities in mountainous areas. The susceptibility maps prepared in this research provide valuable information for landslide hazard management in Lashgarak region of Tehran. This study was conducted to, first, prepare landslide susceptibility maps for Lashgarak region and evaluate landslide effect on mainlines and, second, to analyze the main factors affecting landslide hazard increase in the study area in order to propose efficient strategies for landslide hazard mitigation. A GIS-based multi-criteria decision analysis model (fuzzy logic) is used in the present work for scientific evaluation of landslide susceptible areas in Lashgarak region. To this end, ArcGIS, PCIGeomatica, and IDIRISI software packages were used. Eight information layers were selected for information analysis: ground strength class, slope angle, terrain roughness, normalized difference moisture index, normalized difference vegetation index, distance from fault, distance from the river, and distance from the road. Next, eight different scenarios were created to determine landslide susceptibility of the study area using different operators (intersection (AND), union (OR), algebraic sum (SUM), multiplication (PRODUCT), and different fuzzy gamma values) of fuzzy overlay approach. After that, the performance of various fuzzy operators in landslide susceptibility mapping was empirically compared. The results revealed the excellent consistency of landslide susceptibility map prepared using the fuzzy union (OR) operator with landslide distribution map in the study area. Eventually, the accuracy of landslide susceptibility map prepared using the fuzzy union (OR) operator was evaluated using the frequency ratio diagram. The results showed that frequency values of the landslides gradually increase from “low susceptibility” to high “susceptibility” as 88.34% of the landslides are categorized into two “high” and “very high” susceptibility classes, implying the satisfactory consistency between the landslide susceptibility map prepared using fuzzy union (OR) operator and landslide distribution map.     

A comparison of the GIS based landslide susceptibility assessment methods: multivariate versus bivariate

The purpose of this study is to evaluate and to compare the results of multivariate (logical regression) and bivariate (landslide susceptibility) methods in Geographical Information System (GIS) based landslide susceptibility assessment procedures. In order to achieve this goal the Asarsuyu catchment in NW Turkey was selected as a test zone because of its well-known landslide occurrences interfering with the E-5 highway mountain pass.Two methods were applied to the test zone and two separate susceptibility maps were produced. Following this a two-fold comparison scheme was implemented. Both methods were compared by the Seed Cell Area Indexes (SCAI) and by the spatial locations of the resultant susceptibility pixels.It was found that both of the methods converge in 80% of the area; however, the weighting algorithm in the bivariate technique (landslide susceptibility method) had some severe deficiencies, as the resultant hazard classes in overweighed areas did not converge with the factual landslide inventory map. The result of the multivariate technique (logical regression) was more sensitive to the different local features of the test zone and it resulted in more accurate and homogeneous susceptibility maps.     

Landslide detection and susceptibility mapping using LiDAR and an artificial neural network approach: a case study in the Cuyahoga Valley National Park,Ohio

The purpose of this study was to detect shallow landslides using hillshade maps derived from light detection and ranging (LiDAR)-based digital elevation model (DEM) derivatives. The landslide susceptibility mapping used an artificial neural network (ANN) approach and backpropagation method that was tested in the northern portion of the Cuyahoga Valley National Park (CVNP) located in northeast Ohio. The relationship between landslides and predictor attributes, which describe landform classes using slope, profile and plan curvatures, upslope drainage area, annual solar radiation, and wetness index, was extracted from LiDAR-based DEM using geographic information system (GIS). The approach presented in this paper required a training study area for the development of the susceptibility model and a validation study area to test the model. The results from the validation showed that within the very high susceptibility class, a total of 42.6 % of known landslides that were associated with 1.56 % of total area were correctly predicted. In contrast, the very low susceptibility class that represented 82.68 % of the total area was associated with 1.20 % of known landslides. The results suggest that the majority of the known landslides occur within a small portion of the study area, consistent with field investigation and other studies. Sample probabilistic maps of landslide susceptibility potential and other products from this approach are summarized and presented for visualization to help park officials in effective management and planning.     

基于GIS技术和频率比模型的三峡地区秭归向斜盆地滑坡敏感性评价

以三峡地区秭归向斜盆地为研究区,在野外滑坡编图的基础上,选择地形坡度、地层岩性、河流、道路、坡向、坡型结构6个因子参与滑坡敏感性评价,运用GIS空间分析技术,引入频率比模型,分别计算评价因子的贡献率,并进行叠加分析,最终划分为4个敏感性分区。评价结果表明:极高敏感性区面积占21.39%,滑坡发生面积占61.44%;高敏感性分区面积占24.99%,滑坡发生面积占21.67%;中敏感性分区面积占30.66%,滑坡发生面积占13.19%;低敏感性分区面积占22.93%,滑坡发生面积占3.69%。滑坡敏感性面积累计百分比曲线图表明评价结果具有较高的准确度和可靠性。研究成果可为政府部门减灾防灾工程提供科学支持,为滑坡灾害的预测和管理提供科学依据。     

Landslide susceptibility mapping using GIS and digital photogrammetric techniques: a case study from Ardesen (NE-Turkey)

Ardesen is a settlement area which has been significantly damaged by frequent landslides which are caused by severe rainfalls and result in many casualties. In this study a landslide susceptibility map of Ardesen was prepared using the Analytical Hierarchy Process (AHP) with the help of Geographical Information Systems (GIS) and Digital Photogrametry Techniques (DPT). A landslide inventory, lithology–weathering, slope, aspect, land cover, shear strength, distance to the river, stream density and distance to the road thematics data layers were used to create the map. These layer maps are produced using field, laboratory and office studies, and by the use of GIS and DPT. The landslide inventory map is also required to determine the relationship between these maps and landslides using DPT. In the study field in the Hemsindere Formation there are units that have different weathering classes, and this significantly affects the shear strength of the soil. In this study, shear strength values are calculated in great detail with field and laboratory studies and an additional layer is evaluated with the help of the stability studies used to produce the landslide susceptibility map. Finally, an overlay analysis is carried out by evaluating the layers obtained according to their weight, and the landslide susceptibility map is produced. The study area was classified into five classes of relative landslide susceptibility, namely, very low, low, moderate, high, and very high. Based on this analysis, the area and percentage distribution of landslide susceptibility degrees were calculated and it was found that 28% of the region is under the threat of landslides. Furthermore, the landslide susceptibility map and the landslide inventory map were compared to determine whether the models produced are compatible with the real situation resulting in compatibility rate of 84%. The total numbers of dwellings in the study area were determined one by one using aerial photos and it was found that 30% of the houses, with a total occupancy of approximately 2,300 people, have a high or very high risk of being affected by landslides.     

Use of fuzzy relations to produce landslide susceptibility map of a landslide prone area (West Black Sea Region, Turkey)

Preparation of landslide susceptibility maps is important for engineering geologists and geomorphologists. However, due to complex nature of landslides, producing a reliable susceptibility map is not easy. For this reason, many procedures have been used to produce such maps. In this study, a new attempt is tried to produce landslide susceptibility map of a part of West Black Sea Region of Turkey. To obtain the fuzzy relations for producing the susceptibility map, a landslide inventory database is compiled by both field surveys and airphoto studies. A total of 266 landslides are identified in the study area, and dominant mode of failure is rotational slide while the other mode of failures are soil flow and shallow translational slide. The landslide inventory and the parameter maps are analyzed together using a computer program (FULLSA) developed in this study. The computer program utilizes the fuzzy relations and produces the landslide susceptibility map automatically. According to this map, 9.6% of the study area is classified as very high susceptibility, 10.3% as high susceptibility, 8.9% as moderate susceptibility, 27.5% as low susceptibility and 43.8% as very low susceptibility or nonsusceptible areas. The prediction performance of the susceptibility map is checked by considering actual landslides in the study area. For this purpose, strength of the relation (r_ij) and the root mean square error (RMSE) values are calculated as 0.867 and 0.284, respectively. These values show that the produced landslide susceptibility map in the present study has a sufficient reliability. It is believed that the approach employed in this study mainly prevents the subjectivity sourced from the parameter selection and provides a support to improve the landslide susceptibility mapping studies.     

A subjective and objective integrated weighting method for landslides susceptibility mapping based on GIS

The purpose of this study is to present a weighting method, integrating subjective weight with objective weight, for landslides susceptibility mapping based on geographical information system (GIS). First, the landslide inventory, aspect, slope, proximity to streams of drainage network, proximity to railway, proximity to road, topography, elevation, lithology, tectonic activity and annual precipitation, including their subclasses, were taken as independent landslide causal factors. Second, objective weights of the causal factors were calculated according to the landslide area density based on entropy weighting method, and key factors were selected according to the rank of the objective weights. Third, trapezoidal fuzzy number weighting approach was used to assess the sub-classes of each key factor. Finally, a case study was carried out in Guizhou province, China. A landslide susceptibility map was created using weighted linear combination model based on GIS. Using a predicted map of probability, the study area was classified into four categories of landslide susceptibility: low, moderate, moderate-high, and high.