Landslide inventory and susceptibility assessment using multiple statistical approaches along the Karakoram highway,northern Pakistan

China-Pakistan Economic Corridor(CPEC)is a framework of regional connectivity,which will not only benefit China and Pakistan but will have positive impact on Iran,Afghanistan,India,Central Asian Republic,and the region.The surrounding area in CPEC is prone to frequent disruption by geological hazards mainly landslides in northern Pakistan.Comprehensive landslide inventory and susceptibility assessment are rarely available to utilize for landslide mitigation strategies.This study aims to utilize the high-resolution satellite images to develop a comprehensive landslide inventory and subsequently develop landslide susceptibility maps using multiple techniques.The very high-resolution(VHR)satellite images are utilized to develop a landslide inventory using the visual image classification techniques,historic records and field observations.A total of 1632 landslides are mapped in the area.Four statistical models i.e.,frequency ratio,artificial neural network,weights of evidence and logistic regression were used for landslide susceptibility modeling by comparing the landslide inventory with the topographic parameters,geological features,drainage and road network.The developed landslides susceptibility maps were verified using the area under curve(AUC)method.The prediction power of the model was assessed by the prediction rate curve.The success rate curves show 93%,92.8%,92.7%and 87.4%accuracy of susceptibility maps for frequency ratio,artificial neural network,weights of evidence and logistic regression,respectively.The developed landslide inventory and susceptibility maps can be used for land use planning and landslide mitigation strategies.     

Landslide inventory and susceptibility modelling using geospatial tools,in Hunza-Nagar valley,northern Pakistan

A comprehensive landslide inventory and susceptibility maps are prerequisite for developing and implementing landslide mitigation strategies. Landslide susceptibility maps for the landslides prone regions in northern Pakistan are rarely available. The Hunza-Nagar valley in northern Pakistan is known for its frequent and devastating landslides. In this paper, we have developed a landslide inventory map for Hunza-Nagar valley by using the visual interpretation of the SPOT-5 satellite imagery and mapped a total of 172 landslides. The landslide inventory was subsequently divided into modelling and validation data sets. For the development of landslide susceptibility map seven discrete landslide causative factors were correlated with the landslide inventory map using weight of evidence and frequency ratio statistical models. Four different models of conditional independence were used for the selection of landslide causative factors. The produced landslides susceptibility maps were validated by the success rate and area under curves criteria. The prediction power of the models was also validated with the prediction rate curve. The validation results shows that the success rate curves of the weight of evidence and the frequency models are 82% and 79%, respectively. The prediction accuracy results obtained from this study are 84% for weight of evidence model and 80% for the frequency ratio model. Finally, the landslide susceptibility index maps were classified into five different varying susceptibility zones. The validation and prediction result indicates that the weight of evidence and frequency ratio model are reliable to produce an accurate landslide susceptibility map, which may be helpful for landslides management strategies.     

Regional scale rainfall- and earthquake-triggered landslide susceptibility assessment in Wudu County, China

Wudu County in northwestern China frequently experiences large-scale landslide events.High-magnitude earthquakes and heavy rainfall events are the major triggering factors in the region.The aim of this research is to compare and combine landslide susceptibility assessments of rainfalltriggered and earthquake-triggered landslide events in the study area using Geographical Information System(GIS) and a logistic regression model.Two separate susceptibility maps were produced using inventories reflecting single landslide-triggering events,i.e.,earthquakes and heavy rain storms.Two groups of landslides were utilized: one group containing all landslides triggered by extreme rainfall events between 1995 and 2003 and the other group containing slope failures caused by the 2008 Wenchuan earthquake.Subsequently,the individual maps were combined to illustrate the locations of maximum landslide probability.The use of the resulting three landslide susceptibility maps for landslide forecasting,spatial planning and for developing emergency response actions are discussed.The combined susceptibility map illustrates the total landslide susceptibility in the study area.     

GIS-based Frequency Ratio and Logistic Regression Modelling for Landslide Susceptibility Mapping of Debre Sina Area in Central Ethiopia

Ethiopia has a mountainous landscape which can be divided into the Northwestern and Southeastern plateaus by the Main Ethiopian Rift and Afar Depression. Debre Sina area is located in Central Ethiopia along the escarpment where landslide problem is frequent due to steep slope, complex geology, rift tectonics, heavy rainfall and seismicity. In order to tackle this problem, preparing a landslide susceptibility map is very important. For this, GISbased frequency ratio(FR) and logistic regression(LR) models have been applied using landslide inventory and the nine landslide factors(i.e. lithology, land use, distance from river fault, slope, aspect, elevation, curvature and annual rainfall). Database construction, weighting each factor classes or factors, preparing susceptibility map and validation were the major steps to be undertaken. Both models require a rasterized landslide inventory and landslide factor maps. The former was classified into training and validation landslides. Using FR model, weights for each factor classes were calculated and assigned so that all the weighted factor maps can be added to produce a landslide susceptibility map. In the case of LR model, the entire study area is firstly divided into landslide and non-landslide areas using the training landslides. Then, these areas are changed into landslide and non-landslide points so as to extract the FR maps of the nine landslide factors. Then a linear relationship is established between training landslides and landslide factors in SPSS. Based on this relationship, the final landslide susceptibility map is prepared using LR equation. The success-rate and prediction-rate of FR model were 74.8% and 73.5%, while in case of LR model these were 75.7% and 74.5% respectively. A close similarity in the prediction and validation rates showed that the model is acceptable. Accuracy of LR model is slightly better in predicting the landslide susceptibility of the area compared to FR model.     

Landslide susceptibility mapping using an integrated model of information value method and logistic regression in the Bailongjiang watershed,Gansu Province,China

Bailongjiang watershed in southern Gansu province, China, is one of the most landslide-prone regions in China, characterized by very high frequency of landslide occurrence. In order to predict the landslide occurrence, a comprehensive map of landslide susceptibility is required which may be significantly helpful in reducing loss of property and human life. In this study, an integrated model of information value method and logistic regression is proposed by using their merits at maximum and overcoming their weaknesses, which may enhance precision and accuracy of landslide susceptibility assessment. A detailed and reliable landslide inventory with 1587 landslides was prepared and randomly divided into two groups, (i) training dataset and (ii) testing dataset. Eight distinct landslide conditioning factors including lithology, slope gradient, aspect, elevation, distance to drainages, distance to faults, distance to roads and vegetation coverage were selected for landslide susceptibility mapping. The produced landslide susceptibility maps were validated by the success rate and prediction rate curves. The validation results show that the success rate and the prediction rate of the integrated model are 81.7 % and 84.6 %, respectively, which indicate that the proposed integrated method is reliable to produce an accurate landslide susceptibility map and the results may be used for landslides management and mitigation.     

Landslide Susceptibility Assessment of the Youfang Catchment using Logistic Regression

A detailed landslide susceptibility map was produced in the Youfang catchment using logistic regression method with datasets developed for a geographic information system(GIS).Known as one of the most landslide-prone areas in China, the Youfang catchment of Longnan mountain region,which lies in the transitional area among QinghaiTibet Plateau, loess Plateau and Sichuan Basin, was selected as a representative case to evaluate the frequency and distribution of landslides.Statistical relationships for landslide susceptibility assessment were developed using landslide and landslide causative factor databases.Logistic regression(LR)was used to create the landslide susceptibility maps based on a series of available data sources: landslide inventory; distance to drainage systems, faults and roads; slope angle and aspect; topographic elevation and topographical wetness index, and land use.The quality of the landslide susceptibility map produced in this paper was validated and the result can be used fordesigning protective and mitigation measures against landslide hazards.The landslide susceptibility map is expected to provide a fundamental tool for landslide hazards assessment and risk management in the Youfang catchment.     

Mapping landslide susceptibility at the Three Gorges Reservoir,China, using gradient boosting decision tree,random forest and information value models

This work was to generate landslide susceptibility maps for the Three Gorges Reservoir(TGR) area, China by using different machine learning models. Three advanced machine learning methods, namely, gradient boosting decision tree(GBDT), random forest(RF) and information value(InV) models, were used, and the performances were assessed and compared. In total, 202 landslides were mapped by using a series of field surveys, aerial photographs, and reviews of historical and bibliographical data. Nine causative factors were then considered in landslide susceptibility map generation by using the GBDT, RF and InV models. All of the maps of the causative factors were resampled to a resolution of 28.5 m. Of the 486289 pixels in the area,28526 pixels were landslide pixels, and 457763 pixels were non-landslide pixels. Finally, landslide susceptibility maps were generated by using the three machine learning models, and their performances were assessed through receiver operating characteristic(ROC) curves, the sensitivity, specificity,overall accuracy(OA), and kappa coefficient(KAPPA). The results showed that the GBDT, RF and In V models in overall produced reasonable accurate landslide susceptibility maps. Among these three methods, the GBDT method outperforms the other two machine learning methods, which can provide strong technical support for producing landslide susceptibility maps in TGR.     

基于信息量模型和数据标准化的滑坡易发性评价

本文以北川曲山-擂鼓片区为研究区,将坡度、坡向、高程、地层、距断层的距离、距水系的距离和距道路的距离作为该区域滑坡易发性评价因子。采用信息量模型计算了各项评价因子的信息量值,并运用4种标准化模型对信息量值进行标准化处理。各评价因子的权重由层次分析法（AHP）确定。在GIS中将权重值和各评价因子的标准化信息量值,进行叠加计算得到区域滑坡总信息量值,并基于自然断点法对其进行重分类,将研究区划分为极高易发区、高易发区、中易发区、低易发区和极低易发区5级易发区。将基于4种标准化模型和信息量模型得到的滑坡易发性评价结果进行了对比分析,结果表明：基于最值标准化信息量模型的滑坡易发性评价结果的ROC曲线下面积AUC值为0.807,高于其余模型的AUC值,说明最值标准化信息量模型的滑坡易发性评价效果最好。极高易发区面积占研究区面积的20.03%,离断层和水系较近,主要分布地层为寒武系、志留系和三迭系。研究结果可为区内滑坡风险评价和灾害防治提供参考。     

Landslide Susceptibility Mapping along Bhalubang–Shiwapur Area of Mid-Western Nepal Using Frequency Ratio and Conditional Probability Models

Roads constructed in fragile Siwaliks are prone to large number of instabilities. Bhalubang–Shiwapur section of Mahendra Highway lying in Western Nepal is one of them. To understand the landslide causative factor and to predict future occurrence of the landslides, landslide susceptibility mapping(LSM) of this region was carried out using frequency ratio(FR) and weights-of-evidence(W of E) models. These models are easy to apply and give good results. For this, landslide inventory map of the area was prepared based on the aerial photo interpretation, from previously published/unpublished reposts, and detailed field survey using GPS. About 332 landslides were identified and mapped, among which 226(70%) were randomly selected for model training and the remaining 106(30%) were used for validation purpose. A spatial database was constructed from topographic, geological, and land cover maps. The reclassified maps based on the weight values of frequency ratio and weights-of-evidence were applied to get final susceptibility maps. The resultant landslide susceptibility maps were verified andcompared with the training data, as well as with the validation data. From the analysis, it is seen that both the models were equally capable of predicting landslide susceptibility of the region(W of E model(success rate = 83.39%, prediction rate = 79.59%); FR model(success rate = 83.31%, prediction rate = 78.58%)). In addition, it was observed that the distance from highway and lithology, followed by distance from drainage, slope curvature, and slope gradient played major role in the formation of landsides. The landslide susceptibility maps thus produced can serve as basic tools for planners and engineers to carry out further development works in this landslide prone area.     

GIS based landslide susceptibility mapping of Tevankarai Ar sub-watershed, Kodaikkanal, India using binary logistic regression analysis

Landslide susceptibility mapping is the first step in regional hazard management as it helps to understand the spatial distribution of the probability of slope failure in an area.An attempt is made to map the landslide susceptibility in Tevankarai Ar subwatershed,Kodaikkanal,India using binary logistic regression analysis.Geographic Information System is used to prepare the database of the predictor variables and landslide inventory map,which is used to build the spatial model of landslide susceptibility.The model describes the relationship between the dependent variable(presence and absence of landslide) and the independent variables selected for study(predictor variables) by the best fitting function.A forward stepwise logistic regression model using maximum likelihood estimation is used in the regression analysis.An inventory of 84 landslides and cells within a buffer distance of 10m around the landslide is used as the dependent variable.Relief,slope,aspect,plan curvature,profile curvature,land use,soil,topographic wetness index,proximity to roads and proximity to lineaments are taken as independent variables.The constant and the coefficient of the predictor variable retained by the regression model are used to calculate the probability of slope failure and analyze the effect of each predictor variable on landslide occurrence in thestudy area.The model shows that the most significant parameter contributing to landslides is slope.The other significant parameters are profile curvature,soil,road,wetness index and relief.The predictive logistic regression model is validated using temporal validation data-set of known landslide locations and shows an accuracy of 85.29 %.     

GIS based spatial data analysis for landslide susceptibility mapping

Landslide susceptibility map delineates the potential zones for landslides occurrence. The paper presents a statistical approach through spatial data analysis in GIS for landslide susceptibility mapping in parts of Sikkim Himalaya. Six important causative factors for landslide occurrences were selected and corresponding thematic data layers were prepared in GIS. Topographic maps,satellite image,field data and published maps constitute the input data for thematic layer preparation. Numerical weights for different categories of these factors were determined based on a statistical approach and the weighted thematic layers were integrated in GIS environment to generate the landslide susceptibility map of the area. The landslide susceptibility map classifies the area into five different landslide susceptible zones i.e.,very high,high,moderate,low and very low. This map was validated using the existing landslide distribution in the area.     

Hazard assessment of landslide disaster using information value method and analytical hierarchy process in highly tectonic Chamba region in bosom of Himalaya

The present study is focused on a comparative evaluation of landslide disaster using analytical hierarchy process and information value method for hazard assessment in highly tectonic Chamba region in bosom of Himalaya. During study, the information about the causative factors was generated and the landslide hazard zonation maps were delineated using Information Value Method (IV) and Analytical Hierarchy Process (AHP) using ArcGIS (ESRI). For this purpose, the study area was selected in a part of Ravi river catchment along one of the landslide prone Chamba to Bharmour road corridor of National Highway (NH-154A) in Himachal Pradesh, India. A numeral landslide triggering geoenvironmental factors i.e. slope, aspect, relative relief, soil, curvature, land use and land cover (LULC), lithology, drainage density, and lineament density were selected for landslide hazard mapping based on landslide inventory. Landslide hazard zonation map was categorized namely “very high hazard, high hazard, medium hazard, low hazard, and very low hazard”. The results from these two methods were validated using Area Under Curve (AUC) plots. It is found that hazard zonation map prepared using information value method and analytical hierarchy process methods possess the prediction rate of 78.87% and 75.42%, respectively. Hence, landslide hazard zonation map obtained using information value method is proposed to be more useful for the study area. These final hazard zonation maps can be used by various stakeholders like engineers and administrators for proper maintenance and smooth traffic flow between Chamba and Bharmour cities, which is the only route connecting these tourist places.     

Landslide susceptibility assessment using the certainty factor and analytic hierarchy process

A new approach combining the certainty factor (CF) and analytic hierarchy process (AHP) methods was proposed to assess landslide susceptibility in the Ziyang district, which is situated in the Qin-Ba Mountain region, China. Landslide inventory data were collected based on field investigations and remote sensing interpretations. A total of 791 landslides were identified. A total of 633 landslides were randomly selected from this data set as the training set, and the remaining landslides were used for validation as the test set. Nine factors, including the slope angle, slope aspect, slope curvature, lithology, distance to faults, distance to streams, precipitation, road network intensity degree and land use were chosen as the landslide causal factors for further susceptibility assessment. The weight of each factor and its subclass were calculated by AHP and CF methods. Landslide susceptibility was compared between the bivariate statistical method and the proposed CF-AHP method. The results indicate that the distance to streams, distance to faults and lithology are the most dominant causal factors associated with landslides. The susceptibility zonation was categorized into five classes of landslide susceptibility, i.e., very high, high, moderate, low and very low level. Lastly, the relative operating characteristics (ROC) curve was used to validate the accuracy of the new approach, and the result showed a satisfactory prediction rate of 78.3%, compared to 69.2% obtained with the landslide susceptibility index method. The results indicate that the CF-AHP combined method is more appropriate for assessing the landslide susceptibility in this area.     

Landslide initiation and runout susceptibility modeling in the context of hill cutting and rapid urbanization: a combined approach of weights of evidence and spatial multi-criteria

Rainfall induced landslides are a common threat to the communities living on dangerous hill-slopes in Chittagong Metropolitan Area, Bangladesh. Extreme population pressure, indiscriminate hill cutting, increased precipitation events due to global warming and associated unplanned urbanization in the hills are exaggerating landslide events. The aim of this article is to prepare a scientifically accurate landslide susceptibility map by combining landslide initiation and runout maps. Land cover, slope, soil permeability, surface geology, precipitation, aspect, and distance to hill cut, road cut, drainage and stream network factor maps were selected by conditional independence test. The locations of 56 landslides were collected by field surveying. A weight of evidence (WoE) method was applied to calculate the positive (presence of landslides) and negative (absence of landslides) factor weights. A combination of analytical hierarchical process (AHP) and fuzzy membership standardization (weighs from 0 to 1) was applied for performing a spatial multi-criteria evaluation. Expert opinion guided the decision rule for AHP. The Flow-R tool that allows modeling landslide runout from the initiation sources was applied. The flow direction was calculated using the modified Holmgren’s algorithm. The AHP landslide initiation and runout susceptibility maps were used to prepare a combined landslide susceptibility map. The relative operating characteristic curve was used for model validation purpose. The accuracy of WoE, AHP, and combined susceptibility map was calculated 96%, 97%, and 98%, respectively.     

Geospatial susceptibility mapping of earthquake-induced landslides in Nuweiba area,Gulf of Aqaba,Egypt

Earthquake induced landslides are one of the most severe geo-environmental hazards that cause enormous damage to infrastructure, property, and loss of life in Nuweiba area. This study developed a model for mapping the earthquake-induced landslide susceptibility in Nuweiba area in Egypt with considerations of geological, geomorphological, topographical, and seismological factors. An integrated approach of remote sensing and GIS technologies were applied for that target. Several data sources including Terra SAR-X and SPOT 5 satellite imagery, topographic maps, field data, and other geospatial resources were used to model landslide susceptibility. These data were used specifically to produce important thematic layers contributing to landslide occurrences in the region. A rating scheme was developed to assign ranks for the thematic layers and weights for their classes based on their contribution in landslide susceptibility. The ranks and weights were defined based on the knowledge from field survey and authors experiences related to the study area. The landslide susceptibility map delineates the hazard zones to three relative classes of susceptibility: high, moderate, and low. Therefore, the current approach provides a way to assess landslide hazards and serves for geo-hazard planning and prediction in Nuweiba area.     

GIS-based evaluation on the fault motion-induced coseismic landslides

Earthquake-induced potential landslides are commonly estimated using landslide susceptibility maps. Nevertheless, the fault location is not identified and the ground motion caused by it is unavailable in the map. Thus, potential coseismic landslides for a specific fault motion-induced earthquake could not be predicted using the map. It is meaningful to incorporate the fault location and ground motion characteristics into the landslide predication model. A new method for a specific fault motion-induced coseismic landslide prediction model using GIS (Geographic Information System) is proposed herein. Location of mountain ridges, slope gradients over 45 o , PVGA (Peak Vertical Ground Accelerations) exceeded 0.15 g, and PHGA (Peak Horizontal Ground Accelerations) exceeded 0.25 g of slope units were representing locations that initiated landslides during the 1999 Chi-Chi earthquake in Taiwan. These coseismic landslide characteristics were used to identify areas where landslides occurred during Meishan fault motion-induced strong ground motions in Chiayi County in Taiwan. The strong ground motion (over 8 Gal in the database, 1 Gal = 0.01 m/s 2 , and 1 g = 981 Gal) characteristics were evaluated by the fault length, site distance to the fault, and topography, and their attenuation relations are presented in GIS. The results of the analysis show that coseismic landslide areas could be identified promptly using GIS. The earthquake intensity and focus depth have visible effects on ground motion. The shallower the focus depth, the larger the magnitude increase of the landslides. The GIS-based landslide predication method is valuable combining the geomorphic characteristics and ground motion attenuation relationships for a potential region landslide hazard assessment and in disaster mitigation planning.     

Susceptibility assessment of landslides caused by the wenchuan earthquake using a logistic regression model

The Wenchuan earthquake on May 12, 2008 caused numerous collapses, landslides, barrier lakes, and debris flows. Landslide susceptibility mapping is important for evaluation of environmental capacity and also as a guide for post-earthquake reconstruction. In this paper, a logistic regression model was developed within the framework of GIS to map landslide susceptibility. Qingchuan County, a heavily affected area, was selected for the study. Distribution of landslides was prepared by interpretation of multi-temporal and multi-resolution remote sensing images (ADS40 aerial imagery, SPOT5 imagery and TM imagery, etc.) and field surveys. The Certainly Factor method was used to find the influencial factors, indicating that lithologic groups, distance from major faults, slope angle, profile curvature, and altitude are the dominant factors influencing landslides. The weight of each factor was determined using a binomial logistic regression model. Landslide susceptibility mapping was based on spatial overlay analysis and divided into five classes. Major faults have the most significant impact, and landslides will occur most likely in areas near the faults. Onethird of the area has a high or very high susceptibility, located in the northeast, south and southwest, including 65.3% of all landslides coincident with the earthquake. The susceptibility map can reveal the likelihood of future failures, and it will be useful for planners during the rebuilding process and for future zoning issues.     

GIS-based landslide susceptibility mapping using numerical risk factor bivariate model and its ensemble with linear multivariate regression and boosted regression tree algorithms

In this study, a novel approach of the landslide numerical risk factor(LNRF) bivariate model was used in ensemble with linear multivariate regression(LMR) and boosted regression tree(BRT) models, coupled with radar remote sensing data and geographic information system(GIS), for landslide susceptibility mapping(LSM) in the Gorganroud watershed, Iran. Fifteen topographic, hydrological, geological and environmental conditioning factors and a landslide inventory(70%, or 298 landslides) were used in mapping. Phased array-type L-band synthetic aperture radar data were used to extract topographic parameters. Coefficients of tolerance and variance inflation factor were used to determine the coherence among conditioning factors. Data for the landslide inventory map were obtained from various resources, such as Iranian Landslide Working Party(ILWP), Forestry, Rangeland and Watershed Organisation(FRWO), extensive field surveys, interpretation of aerial photos and satellite images, and radar data. Of the total data, 30% were used to validate LSMs, using area under the curve(AUC), frequency ratio(FR) and seed cell area index(SCAI).Normalised difference vegetation index, land use/land cover and slope degree in BRT model elevation, rainfall and distance from stream were found to be important factors and were given the highest weightage in modelling. Validation results using AUC showed that the ensemble LNRF-BRT and LNRFLMR models(AUC = 0.912(91.2%) and 0.907(90.7%), respectively) had high predictive accuracy than the LNRF model alone(AUC = 0.855(85.5%)). The FR and SCAI analyses showed that all models divided the parameter classes with high precision. Overall, our novel approach of combining multivariate and machine learning methods with bivariate models, radar remote sensing data and GIS proved to be a powerful tool for landslide susceptibility mapping.     

Relationships between landslide types and topographic attributes in a loess catchment, China

Topographic attributes have been identified as the most important factor in controlling the initiation and distribution of shallow landslides triggered by rainfall.As a result,these landslides influence the evolution of local surface topography.In this research,an area of 2.6 km 2 loess catchment in the Huachi County was selected as the study area locating in the Chinese Loess Plateau.The landslides inventory and landslide types were mapped using global position system(GPS) and field mapping.The landslide inventory shows that these shallow landslides involve different movement types including slide,creep and fall.Meanwhile,main topographic attributes were generated based on a high resolution digital terrain model(5 m × 5 m),including aspect,slope shape,elevation,slope angle and contributing area.These maps were overlaid with the spatial distributions of total landslides and each type of landslides in a geographic information system(GIS),respectively,to assess their spatial frequency distributions and relative failure potentials related to these selected topographic attributes.The spatial analysis results revealed that there is a close relation between the topographic attributes of the postlandsliding local surface and the types of landslide movement.Meanwhile,the types of landslide movement have some obvious differences in local topographic attributes,which can influence the relative failure potential of different types of landslides.These results have practical significance to mitigate natural hazard and understandgeomorphologic process in thick loess area.     

霍西煤矿区地表滑坡灾害敏感性数值建模与等级划分

本文以山西省霍西煤矿区为研究区,利用遥感和GIS方法对滑坡灾害的敏感性进行了数值建模与定量评价。利用交叉检验方法构建了径向基核函数支持向量机滑坡敏感性评价模型,并基于拟合精度对模型进行了定量评价;对各评价因子在模型中的重要性进行对比分析;基于空间分辨率为30m的评价因子,通过径向基核函数支持向量机模型获得了霍西煤矿区滑坡敏感性指数值,并利用分位数法将霍西煤矿区的滑坡敏感性分为极高、高、中和低4个等级。结果表明：拟合精度建模阶段和验证阶段分别为87.22%和70.12%;与滑坡敏感性关系最密切的5个评价因子依次是岩性、距道路距离、坡向、高程和土地利用类型;极高和高敏感区域分布了93.49%的滑坡点,面积占总面积的50.99%,是比较合理的分级方案。本研究不仅可以为研究区人工边坡调查和煤矿资源合理开采提供借鉴,对相似矿区的相关工作也具有参考价值。