RS与GIS支持下的汶川县城周边地质灾害危险性评价

地质灾害危险性评价是防灾减灾工作的重要依据。本文以汶川县城周边64 km2为例,应用遥感信息提取技术与GIS空间分析方法,根据IKONOS遥感图像和地形图及野外调查资料,提取了崩塌和滑坡易发性评价因子,采用信息量法确定了因子分值,计算了崩塌和滑坡易发性,并分别提出崩塌和滑坡的危险性计算方法,形成了汶川地区崩塌和滑坡危险性分区图。研究结果表明:新的崩塌和滑坡危险性评价方法能够反映区内地质灾害危险程度,该方法可行,结果合理,这为中、大比例尺区域范围内地质灾害危险性研究提供了有益的思路。     

Landslide hazard zonation using analytical hierarchy process along National Highway-3 in mid Himalayas of Himachal Pradesh,India

Landslide zonation studies emphasize on preparation of landslide hazard zonation maps considering major instability factors contributing to occurrence of landslides. This paper deals with geographic information system-based landslide hazard zonation in mid Himalayas of Himachal Pradesh from Mandi to Kullu by considering nine relevant instability factors to develop the hazard zonation map. Analytical hierarchy process was applied to assign relative weightages over all ranges of instability factors of the slopes in study area. To generate landslide hazard zonation map, layers in geographic information system were created corresponding to each instability factor. An inventory of existing major landslides in the study area was prepared and combined with the landslide hazard zonation map for validation purpose. The validation of the model was made using area under curve technique and reveals good agreement between the produced hazard map and previous landslide inventory with prediction accuracy of 79.08%. The landslide hazard zonation map was classified by natural break classifier into very low hazard, low hazard, moderate hazard, high hazard and very high landslide hazard classes in geographic information system depending upon the frequency of occurrence of landslides in each class. The resultant hazard zonation map shows that 14.30% of the area lies in very high hazard zone followed by 15.97% in high hazard zone. The proposed model provides the best-fit classification using hierarchical approach for the causative factors of landslides having complex structure. The developed hazard zonation map is useful for landslide preparedness, land-use planning, and social-economic and sustainable development of the region.     

The performance evaluation of information value,density area,LNRF, and frequency ratio methods for landslide zonation at Miandarband area,Kermanshah Province,Iran

Landslide is one of the natural disasters which causes a lot of annual damage directly or indirectly in the world. Many planned areas, especially in hilly regions, are prone to different types of landslides; therefore, landslide susceptibility maps become an urgent issue, so that landslide damages and impact can be minimized. The best method for studying landslides, which has long been of interest to researchers, is hazard zonation. In this method, due to the affecting factors in landslide occurrence, study areas are classified into areas with low to very high risk. Different methods have been developed for this purpose. In this paper, the four bivariate statistical methods namely information value, density area, LNRF, and frequency ratio are used to investigate the hazard zonation of landslide in Miandarband located north of the Kermanshah Province. The density ratio (D _r) and Qs values for information value, density area, frequency ratio, and LNRF methods used in this study were calculated to be 2.245312, 0.98146; 2.857816, 1.071185; 2.858085, 0.783945; and 2.418375, 1.070928, respectively. The results indicate that although there are minor differences, the frequency ratio method compared to the density area method that was used for the study of landslide zonation presents better results.     

Rainfall event-based landslide susceptibility zonation mapping

Landslide susceptibility assessment is a major research topic in geo-disaster management. In recent days, various landslide susceptibility and landslide hazard assessment methodologies have been introduced with diverse thoughts of assessment and validation method. Fundamentally, in landslide susceptibility zonation mapping, the susceptibility predictions are generally made in terms of likelihoods and probabilities. An overview of landslide susceptibility zoning practices in the last few years reveals that susceptibility maps have been prepared to have different accuracies and reliabilities. To address this issue, the work in this paper focuses on extreme event-based landslide susceptibility zonation mapping and its evaluation. An ideal terrain of northern Shikoku, Japan, was selected in this study for modeling and event-based landslide susceptibility mapping. Both bivariate and multivariate approaches were considered for the zonation mapping. Two event-based landslide databases were used for the susceptibility analysis, while a relatively new third event landslide database was used in validation. Different event-based susceptibility zonation maps were merged and rectified to prepare a final susceptibility zonation map, which was found to have an accuracy of more than 77 %. The multivariate approach was ascertained to yield a better prediction rate. From this study, it is understood that rectification of susceptibility zonation map is appropriate and reliable when multiple event-based landslide database is available for the same area. The analytical results lead to a significant understanding of improvement in bivariate and multivariate approaches as well as the success rate and prediction rate of the susceptibility maps.     

Landslide vulnerability and risk assessment for multi-hazard scenarios using airborne laser scanning data (LiDAR)

Landslide hazard, vulnerability, and risk-zoning maps are considered in the decision-making process that involves land use/land cover (LULC) planning in disaster-prone areas. The accuracy of these analyses is directly related to the quality of spatial data needed and methods employed to obtain such data. In this study, we produced a landslide inventory map that depicts 164 landslide locations using high-resolution airborne laser scanning data. The landslide inventory data were randomly divided into a training dataset: 70 % for training the models and 30 % for validation. In the initial step, a susceptibility map was developed using logistic regression approach in which weights were assigned to every conditioning factor. A high-resolution airborne laser scanning data (LiDAR) was used to derive the landslide conditioning factors for the spatial prediction of landslide hazard areas. The resultant susceptibility was validated using the area under the curve method. The validation result showed 86.22 and 84.87 % success and prediction rates, respectively. In the second stage, a landslide hazard map was produced using precipitation data for 15 years. The precipitation maps were subsequently prepared and show two main categories (two temporal probabilities) for the study area (the average for any day in a year and abnormal intensity recorded in any day for 15 years) and three return periods (15-, 10-, and 5-year periods). Hazard assessment was performed for the entire study area. In the third step, an element at risk map was prepared using LULC, which was considered in the vulnerability assessment. A vulnerability map was derived according to the following criteria: cost, time required for reconstruction, relative risk of landslide, risk to population, and general effect to certain damage. These criteria were applied only on the LULC of the study area because of lack of data on the population and building footprint and types. Finally, risk maps were produced using the derived vulnerability and hazard information. Thereafter, a risk analysis was conducted. The LULC map was cross-matched with the results of the hazard maps for the return period, and the losses were aggregated for the LULC. Then, the losses were calculated for the three return periods. The map of the risk areas may assist planners in overall landslide hazard management.     

Assessing landslide exposure in areas with limited landslide information

Landslide risk assessment is often a difficult task due to the lack of temporal data on landslides and triggering events (frequency), run-out distance, landslide magnitude and vulnerability. The probability of occurrence of landslides is often very difficult to predict, as well as the expected magnitude of events, due to the limited data availability on past landslide activity. In this paper, a qualitative procedure for assessing the exposure of elements at risk is presented for an area of the Apulia region (Italy) where no temporal information on landslide occurrence is available. Given these limitations in data availability, it was not possible to produce a reliable landslide hazard map and, consequently, a risk map. The qualitative analysis was carried out using the spatial multi-criteria evaluation method in a global information system. A landslide susceptibility composite index map and four asset index maps (physical, social, economic and environmental) were generated separately through a hierarchical procedure of standardising and weighting. The four asset index maps were combined in order to obtain a qualitative weighted assets map, which, combined with the landslide susceptibility composite index map, has provided the final qualitative landslide exposure map. The resulting map represents the spatial distribution of the exposure level in the study area; this information could be used in a preliminary stage of regional planning. In order to demonstrate how such an exposure map could be used in a basic risk assessment, a quantification of the economic losses at municipal level was carried out, and the temporal probability of landslides was estimated, on the basis of the expert knowledge. Although the proposed methodology for the exposure assessment did not consider the landslide run-out and vulnerability quantification, the results obtained allow to rank the municipalities in terms of increasing exposure and risk level and, consequently, to identify the priorities for designing appropriate landslide risk mitigation plans.     

Landslide hazard mapping along national highway-154A in Himachal Pradesh,India using information value and frequency ratio

For the socio-economic development of a country, the highway network plays a pivotal role. It has therefore become an imperative to have landslide hazard assessment along these roads to provide safety. The current study presents landslide hazard zonation maps, based on the information value method and frequency ratio method using GIS on 1:50,000 scale by generating the information about the landslide influencing factors. The study was carried out in the year 2017 on a part of Ravi river catchment along one of the landslide-prone Chamba to Bharmour road corridor of NH-154A in Himachal Pradesh, India. A number of landslide triggering geo-environmental factors like “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. The landslide inventory has been developed using satellite imagery, Google earth and by doing exhaustive field surveys. A digital elevation model was used to generate slope gradient, slope aspect, curvature, and relative relief map of the study area. The other information, i.e., soil maps, geological maps, and toposheets, have been collected from various departments. The 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 have been validated using area under curve (AUC) method. It has been found that hazard zonation map prepared using frequency ratio model had a prediction rate of 75.37% while map prepared using information value method had prediction rate of 78.87%. Hence, on the basis of prediction rate, the landslide hazard zonation map, obtained using information value method, was experienced to be more suitable for the study area.     

Regional-scale landslide activity and landslide susceptibility zonation in the Nepal Himalaya

Landslide susceptibility zonation mapping is a fundamental procedure for geo-disaster management in tropical and sub-tropical regions. Recently, various landslide susceptibility zonation models have been introduced in Nepal with diverse approaches of assessment. However, validation is still a problem. Additionally, the role of various predisposing causative parameters for landslide activity is still not well understood in the Nepal Himalaya. To address these issues of susceptibility zonation and landslide activity, about 4,000 km² area of central Nepal was selected for regional-scale assessment of landslide activity and susceptibility zonation mapping. In total, 655 new landslides and 9,229 old landslides were identified with the study area with the help of satellite images, aerial photographs, field data and available reports. The old landslide inventory was “blind landslide database” and could not explain the particular rainfall event responsible for the particular landslide. But considering size of the landslide, blind landslide inventory was reclassified into two databases: short-duration high-intensity rainfall-induced landslide inventory and long-duration low-intensity rainfall-induced landslide inventory. These landslide inventory maps were considered as proxy maps of multiple rainfall event-based landslide inventories. Similarly, all 9,884 landslides were considered for the activity assessment of predisposing causative parameters. For the Nepal Himalaya, slope, slope aspect, geology and road construction activity (anthropogenic cause) were identified as most affective predisposing causative parameters for landslide activity. For susceptibility zonation, multivariate approach was considered and two proxy rainfall event-based landslide databases were used for the logistic regression modelling, while a relatively recent landslide database was used in validation. Two event-based susceptibility zonation maps were merged and rectified to prepare the final susceptibility zonation map and its prediction rate was found to be more than 82 %. From this work, it is concluded that rectification of susceptibility zonation map is very appropriate and reliable. The results of this research contribute to a significant improvement in landslide inventory preparation procedure, susceptibility zonation mapping approaches as well as role of various predisposing causative parameters for the landslide activity.     

A comparative evaluation of heuristic and bivariate statistical modelling for landslide susceptibility mappings in Ghurmi–Dhad Khola, east Nepal

Landslides are one of the major natural disasters that occur in the Himalayan range with recurring frequency, causing enormous loss of life and property every year. Preparation of landslide inventory maps and landslide susceptibility zonation maps are the important tasks to be taken into account initially for safe mitigation measures. The present paper focuses on landslide susceptibility maps of the Ghurmi–Dhad Khola area, east Nepal, using Geographic Information System. For this purpose, the landslide susceptibility maps are prepared by using the heuristic and bivariate statistical methods. The parameters considered for the study are slope angle, slope aspect, elevation, distance from drainage, geology, land cover, rock and soil type, and distance from faults and folds. The landslide susceptibility zonation map produced from the heuristic method shows that 42.59 % of the observed landslide falls under the very high susceptible zone and 33.00 % under the high susceptible zone. Likewise, the landslide susceptibility zonation map produced from the bivariate method depicts that 44.19 % of the observed landslide falls under the very high susceptible zone and 31.59 % under the high susceptible zone. Both the landslide susceptibility zonation maps are identical, and success rates of both the maps are above 80 %. While comparing the landslide susceptibility maps obtained from two different methods, about 78 % of the study area falls in the identical susceptible zones. Special attention should be taken into consideration for the construction works in the areas which have been spatially agreed as very high and high susceptible zones from both techniques. Moreover, these maps can be used for slope management, land use planning, disaster management planning, etc., by the concerned authorities.     

Inclusion of earthquake strong ground motion in a geographic information system-based landslide susceptibility zonation in Garhwal Himalayas

Garhwal Himalayas are seismically very active and simultaneously suffering from landslide hazards. Landslides are one of the most frequent natural hazards in Himalayas causing damages worth more than one billion US$ and around 200 deaths every year. Thus, it is of paramount importance to identify the landslide causative factors to study them carefully and rank them as per their influence on the occurrence of landslides. The difference image of GIS-derived landslide susceptibility zonation maps prepared for pre- and post-Chamoli earthquake shows the effect of seismic shaking on the occurrence of landslides in the Garhwal Himalaya. An attempt has been made to incorporate seismic shaking parameters in terms of peak ground acceleration with other static landslide causative factors to produce landslide susceptibility zonation map in geographic information system environment. In this paper, probabilistic seismic hazard analysis has been carried out to calculate peak ground acceleration values at different time periods for estimating seismic shaking conditions in the study area. Further, these values are used as one of the causative factors of landslides in the study area and it is observed that it refines the preparation of landslide susceptibility zonation map in seismically active areas like Garhwal Himalayas.     

基于GIS与ANN模型的地震滑坡易发性区划

基于遥感数据、地理信息系统(GIS)技术和人工神经网络(ANN)模型,开展地震滑坡易发性区划研究.2010年4月14日玉树地震后,基于航片与卫星影像目视解译,并辅以野外调查的方法,在地震区圈定了2036处地震诱发滑坡.选择高程、坡度、坡向、斜坡曲率、坡位、与水系距离、地层岩性、与断裂距离、与公路距离、归一化植被指数(NDVI)、与同震地表破裂距离、地震动峰值加速度(PGA)共12个因子作为地震滑坡易发性评价因子.这些因子均是应用GIS技术与遥感影像处理技术,基于地形数据、地质数据、遥感数据得到.训练样本中的滑动样本有两组,一组是滑坡区整个单滑坡体的质心位置,另一组是滑坡滑源区滑前的坡体高程最高的位置.应用这12个影响因子,分别采用这两组评价样本,基于ANN模型建立地震滑坡易发性索引图,基于GIS工具建立地震滑坡易发性分级图.分别应用训练样本中滑坡分布的点数据去检验各自的结果正确率,正确率分别为81.53％与81.29％,表明ANN模型是一种高效科学的地震滑坡易发性区划模型.     

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.     

Assessing the influence of morpho-structural setting on landslide abundance

Knowing the factors that influence landslide abundance and distribution is important to evaluate landslide susceptibility and hazard. Visual interpretation of aerial photographs (API) can be used to collect spatially distributed information on bedding attitude (BA), in an area. Where a map of the location of bedding traces (BTs), i.e. lines showing the intersection of bedding planes with the local topography, is available, the map can be used to obtain BA point data and to prepare maps showing morpho-structural domains. The possibility of using BA maps to investigate the influence of morpho-structural settings on landslide abundance is hampered by the lack of understanding of the influence of the length of the BTs, and of the parameters used to interpolate the BA data on the structural zonation. To investigate the problem, we used information on 207 BTs obtained through API in the Collazzone area, Central Italy, and we prepared 150 maps showing BA information. This was accomplished using 15 different values for the segmentation length of the BTs (S), and 10 different values for the tension parameter (T) used for the interpolation. We compare the results against previous results obtained for the same area adopting a heuristic approach to the segmentation of the same set of BTs. Next, we compare the geographical distribution of old deep-seated, deep-seated and shallow landslides in five morpho-structural domains in the study area, and we analyse the influence of the structural settings on the abundance of the different types of landslides.     

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.     

Landslide Hazard and Risk Assessment on the Northern Slope of Mt. Changbai, China

Landslide hazard and risk assessment on the northern slope of Mt. Changbai, a well-known tourist attraction near the North Korean-Chinese border, are assessed. This study is divided into two parts, namely, landslide hazard zonation and risk assessment. The 1992 Anbalagan and Singh method of landslide hazard zonation （LHZ） was modified and used in this area. In this way, an Associative Analysis Method was used in representative areas to get a measure for controlling factors （slope gradient, relative relief, vegetation, geology, discontinuity development, weak layer thickness and ground water）. For the membership degree of factor to slope failure, the middle range of limited values was used to calculate LHZ. Based on an estimation of the potential damage from slope failure, a reasonable risk assessment map was obtained using the relationship of potential damage and probable hazard to aid future planning and prediction and to avert loss of life.     

Use of quantitative landslide hazard and risk information for local disaster risk reduction along a transportation corridor: a case study from Nilgiri district, India

The objective of analyzing hazard and risk in an area is to utilize the result in selecting appropriate landslide risk reduction strategies. However, this does not happen always, and most often results of the hazard and risk analysis remain at an academic level. The under or non-utilization of results in pre-disaster planning could be due to several reasons, including difficulties in understanding the scientific content/meaning of the models, and lack of information on the practical significance and utility of the models. In this study, an attempt is made to highlight the uses of hazard and risk information in different landslide risk reduction strategies along a transportation corridor in Nilgiri, India. At first, a quantitative analysis of landslide hazard and risk was made. The obtained information was then incorporated in risk reduction options such as land use zoning, engineering solutions, and emergency preparedness. For emergency preparedness, the perception of the local Nilgiri communities toward landslide risk was evaluated and simplified maps were generated for the benefit and understanding of end users. A rainfall threshold-based early warning system was presented, which could be used in risk awareness programs involving public participation. The use of quantitative risk information in the cost-benefit analysis for the planning of structural measures to protect the road and railway alignments was also highlighted, and examples were shown how the transport organizations could implement these measures. Finally, the study provided examples of the utility of hazard and risk information for spatial planning and zoning, indicating areas where the landslide hazard is too high for planning future developments.     

Pathways for adaptive and integrated disaster resilience

The GIS-multicriteria decision analysis (GIS-MCDA) technique is increasingly used for landslide hazard mapping and zonation. It enables the integration of different data layers with different levels of uncertainty. In this study, three different GIS-MCDA methods were applied to landslide susceptibility mapping for the Urmia lake basin in northwest Iran. Nine landslide causal factors were used, whereby parameters were extracted from an associated spatial database. These factors were evaluated, and then, the respective factor weight and class weight were assigned to each of the associated factors. The landslide susceptibility maps were produced based on weighted overly techniques including analytic hierarchy process (AHP), weighted linear combination (WLC) and ordered weighted average (OWA). An existing inventory of known landslides within the case study area was compared with the resulting susceptibility maps. Respectively, Dempster-Shafer Theory was used to carry out uncertainty analysis of GIS-MCDA results. Result of research indicated the AHP performed best in the landslide susceptibility mapping closely followed by the OWA method while the WLC method delivered significantly poorer results. The resulting figures are generally very high for this area, but it could be proved that the choice of method significantly influences the results.     

Landslides in the Mailuu-Suu Valley, Kyrgyzstan—Hazards and Impacts

Mailuu-Suu is a former uranium mining area in Kyrgyzstan (Central Asia) at the northern border of the Fergana Basin. This region is particularly prone to landslide hazards and, during the last 50 years, has experienced severe landslide disasters in the vicinity of numerous nuclear waste tailing dams. Due to its critical situation, the Mailuu-Suu region was and still is the target area for several risk assessment projects. This paper provides a brief review of previous studies, past landslide events and a discussion on possible future risk scenarios. Various aspects of landslide hazard and related impacts in the Mailuu-Suu Valley are analyzed in detail: landslide susceptibility, historical evolution of landslide activity, size-frequency relationship, river damming and flooding as well as impacts on inhabited areas and nuclear waste storage zones. The study was carried out with standard remote sensing tools for the processing of satellite imagery and the construction of digital elevation models (DEMs). The processed inputs were combined on a GIS platform with digital landslide distribution maps of 1962, 1977, and 2003, digitized geological and geographic maps, and information from landslide monitoring and geophysical investigation.As a result, various types of landslide susceptibility maps based on conditional analysis (CA) are presented as well as predictions of future landslide activity and related damming potential and their possible impact on the population. For some risk scenarios, remediation and prevention measures are suggested.     

Application of an analytical hierarchical process approach for landslide susceptibility mapping in A Luoi district,Thua Thien Hue Province,Vietnam

A Luoi is a Vietnamese–Laotian border district situated in the western part of Thua Thien Hue province, central Vietnam, where landslides occur frequently and seriously affect local living conditions. This study focuses on the spatial analysis of landslide susceptibility in this 263-km² area. To analyze landslide manifestation in the study area, causative factor maps are derived of slope angle, weathering, land use, geomorphology, fault density, geology, drainage distance, elevation, and precipitation. The analytical hierarchical process approach is used to combine these maps for landslide susceptibility mapping. A landslide susceptibility zonation map with four landslide susceptibility classes, i.e. low, moderate, high, and very high susceptibility for landsliding, is derived based on the correspondence with an inventory of observed landslides. The final map indicates that about 37% of the area is very highly susceptible for landsliding and about 22% is highly susceptible, which means that more than half of the area should be considered prone to landsliding.     

Application of a RS- and GIS-based semi-quantitative approach (analytical hierarchy process – AHP) in landslide hazard risk assessment of the Shivkhola Watershed,Darjiling Himalaya

The present study deals with the application of analytical hierarchy process to prepare landslide hazard risk map of the Shivkhola Watershed applying remote sensing and geographic information system (GIS). Firstly, to integrate all the required thematic data layers and to prepare landslide susceptibility map, prioritised class rating value and prioritised factor rating value were obtained by developing couple-comparing matrix with a reasonable consistency and with the help of MATLAB software after Saaty. Three important risk factor/element maps, that is, weighted land use/land cover map, road contributing area map and settlement density map, were developed and their weighted linear combination was performed to prepare landslide risk exposure map. Then by integrating landslide susceptibility map and landslide risk exposure map, a classification was incorporated on ARC GIS Platform to prepare landslide hazard risk map. To evaluate the validity of the landslide hazard risk map, probability/chance of landslide hazard risk event has been estimated by means of frequency ratio between landslide hazard risk area (%) and number of risk events (%) for each landslide hazard risk class. Finally, an accuracy assessment was also made on ERDAS Imagine (8.5) which depicts that the classification accuracy of the landslide hazard risk map was 92.89 with overall Kappa statistics of 0.8929.