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 susceptibility mapping using Bayesian approach in the Sultan Mountains (Akşehir,Turkey)

Landslides cause heavy damage to property and infrastructure, in addition to being responsible for the loss of human lives in many parts of the Turkey. The paper presents GIS-based spatial data analysis for landslide susceptibility mapping in the regions of the Sultan Mountains, West of Akşehir, and central part of Turkey. Landslides occur frequently in the area and seriously affect local living conditions. Therefore, spatial analysis of landslide susceptibility in the Sultan Mountains is important. The relationships between landslide distributions with the 19 landslide affecting parameters were analysed using a Bayesian model. In the study area, 90 landslides were observed. The landslides were randomly subdivided into 80 training landslides and 10 test landslides. A landslide susceptibility map was produced by using the training landslides. The test landslides were used in the accuracy control of the produced landslide susceptibility map. Approximately 9% of the study area was classified as high susceptibility zone. Medium, low and very low susceptibility zones covered 8, 23 and 60% of the study area, respectively. Most of the locations of the observed landslides actually fall into moderate (17.78%) and high (77.78. %) susceptibility zones of the produced landslide susceptibility map. This validates the applicability of proposed methods, approaches and the classification scheme. The high susceptibility zone is along both sides of the Akşehir Fault and at the north-eastern slope of the Sultan Mountains. It was determined that the surface area of the Harlak and Deresenek formations, which have attained lithological characteristics of clayey limestone with a broken and separated base, and where area landslides occur, possesses an elevation of 1,100–1,600 m, a slope gradient of 25°–35° and a slope aspect of 22.5°–157.5° facing slopes.     

A distributed hydrological–geotechnical model using satellite-derived rainfall estimates for shallow landslide prediction system at a catchment scale

This paper describes the potential applicability of a hydrological–geotechnical modeling system using satellite-based rainfall estimates for a shallow landslide prediction system. The physically based distributed model has been developed by integrating a grid-based distributed kinematic wave rainfall-runoff model with an infinite slope stability approach. The model was forced by the satellite-based near real-time half-hourly CMORPH global rainfall product prepared by NOAA-CPC. The method combines the following two model outputs necessary for identifying where and when shallow landslides may potentially occur in the catchment: (1) the time-invariant spatial distribution of areas susceptible to slope instability map, for which the river catchment is divided into stability classes according to the critical relative soil saturation; this output is designed to portray the effect of quasi-static land surface variables and soil strength properties on slope instability and (2) a produced map linked with spatiotemporally varying hydrologic properties to provide a time-varying estimate of susceptibility to slope movement in response to rainfall. The proposed hydrological model predicts the dynamic of soil saturation in each grid element. The stored water in each grid element is then used for updating the relative soil saturation and analyzing the slope stability. A grid of slope is defined to be unstable when the relative soil saturation becomes higher than the critical level and is the basis for issuing a shallow landslide warning. The method was applied to past landslides in the upper Citarum River catchment (2,310 km²), Indonesia; the resulting time-invariant landslide susceptibility map shows good agreement with the spatial patterns of documented historical landslides (1985–2008). Application of the model to two recent shallow landslides shows that the model can successfully predict the effect of rainfall movement and intensity on the spatiotemporal dynamic of hydrological variables that trigger shallow landslides. Several hours before the landslides, the model predicted unstable conditions in some grids over and near the grids at which the actual shallow landslides occurred. Overall, the results demonstrate the potential applicability of the modeling system for shallow landslide disaster predictions and warnings.     

Physically based modelling of shallow landslide sediment yield at a catchment scale

 A shallow landslide erosion and sediment yield component, applicable at the basin scale, has been incorporated into the physically based, spatially distributed, hydrological and sediment transport modelling system, SHETRAN. The component determines when and where landslides occur in a basin in response to time-varying rainfall and snowmelt, the volume of material eroded and released for onward transport, and the impact on basin sediment yield. Derived relationships are used to link the SHETRAN grid resolution (up to 1 km), at which the basin hydrology and final sediment yield is modelled, to a subgrid resolution (typically around 10–100 m) at which landslide occurrence and erosion is modelled. The subgrid discretization, landslide susceptibility and potential landslide impact are determined in advance using a geographic information system (GIS), with SHETRAN then providing information on temporal variation in the factors controlling landsliding. The ability to simulate landslide sediment yield is demonstrated by a hypothetical application based on a catchment in Scotland. Received: 30 October 1996 · Accepted: 25 June 1997     

基于有效降雨强度的滑坡灾害危险性预警

选取湖北省恩施地区1 000 km2区域作为典型研究区, 在全面分析该区域历史滑坡资料的基础上, 根据该区滑坡生成与地层岩性之间的关系, 将研究区地层划分为高、中、低3类易发性岩组.分岩组统计降雨监测数据与历史滑坡信息, 得出有效降雨强度与关键降雨持续时间的散点图, 由此确定不同滑坡发生概率的有效降雨强度阈值, 提出该区的滑坡灾害危险性预警判别模型.基于样本区统计数据建立滑坡预测指标体系, 运用GIS得出研究区域的滑坡空间易发性区划结果, 并根据不同易发岩组-有效降雨强度模型, 叠加滑坡灾害易发性分区结果与降雨危险性预警等级分级结果, 对研究区的滑坡灾害危险性进行了预测预警.结果表明: 不同易发岩组-有效降雨强度模型所得预警结果与实际情况吻合, 预警模型具有考虑全面和预警精度高的特点, 在实际预警中切实可用.      

浙江省永嘉县区域滑坡灾害人口易损性评价和伤亡风险预测

区域滑坡灾害人口易损性及人口伤亡风险预测研究是区域滑坡灾害预警预报工作的一个重要环节,该研究对提高预警预报工作的针对性和有效性具有关键作用.在对浙江省永嘉县有关资料进行分析的基础上,从研究区人口年龄结构、居民对滑坡灾害风险的防范意识、政府对滑坡灾害的重视程度及滑坡灾害预警预报体系的完善程度4个方面评价了研究区人口易损性,并给出了计算人口易损性的公式,据此得到了永嘉县人口易损性分布图.根据永嘉县的实际情况,提出了耕地人口密度的概念.综合人口易损性分布图、人口密度分布图和滑坡灾害易发性预测图得到了研究区受威胁人口伤亡风险预测图,为当地政府职能部门实施滑坡灾害风险的控制和管理提供决策依据.     

Canelles landslide: modelling rapid drawdown and fast potential sliding

A large landslide (40 × 10⁶ m³) was reactivated on the left bank of Canelles reservoir, Spain. The instability was made evident after a considerable reduction of the reservoir level. The drawdown took place during the summer of 2006 after several years of high water levels. The drawdown velocity reached values between 0.5 and 1.2 m/day (registered at low elevations). The paper reports the geological and geotechnical investigations performed to define the movement. The geometry of the slip surface was established from the detailed analysis of the continuous cores recovered in deep borings and from limited information provided by inclinometers. Deep piezometric records provided also valuable information on the pore water pressure in the vicinity of the failure surface. These data allowed validating a flow–deformation coupled calculation model, which takes into account the changes in water level that occurred 4  years previous to the failure as well as the average rainfall. The analysis indicates that the most likely reason for the instability is the rapid drawdown that took place during the summer of 2006. The potential sudden acceleration of the slide is also analysed in the paper introducing coupled thermal hydraulic and mechanical effects that may develop at the basal shearing surface of the sliding mass. The results indicate that the slide velocity may reach values around 16 m/s when displacement reaches 250 m.     

基于投影寻踪模型的滑坡危险性等级评价

危险性评价是滑坡灾害预防与减灾工作首要解决的重要内容.在地理信息系统技术支持下, 以山地灾害频发区——小江流域作为研究对象, 选取坡度、土体粘聚力和内摩擦角这3个评价指标构建滑坡危险性分级评价指标体系, 将投影寻踪技术运用到滑坡危险性等级评价中, 对评价样本的各指标因素进行线性投影, 以最优投影方向所对应的投影特征值作为评价依据, 建立了滑坡危险性等级综合评价模型, 绘制了滑坡危险性等级分布图.结果表明: 研究区极高危险区、高危险区、中等危险区、低危险区和极低危险区的面积比例为14.28∶9.41∶69.12∶7.00∶0.19;根据所建立的5级评价指标体系对研究区60个土质滑坡点资料进行了验证, 在占研究区总面积23.69%的高、极高危险区的小范围内, 实际发生土质滑坡数量45个, 占总土质滑坡数量的75.00%;中等危险性级别以上区域拥有的土质滑坡数量占全部土质滑坡的96.67%;不同危险性级别的滑坡体积方量统计结果表明, 滑坡体积方量密度随危险性级别的提高而迅速增加.对比评价结果及实测结果可知, 投影寻踪分级结果符合实际情况, 证实了该方法的正确性, 为滑坡危险性评价提供了一条新思路.      

Hazard assessment of the Tidal Inlet landslide and potential subsequent tsunami, Glacier Bay National Park, Alaska

An unstable rock slump, estimated at 5 to 10 × 10⁶ m³, lies perched above the northern shore of Tidal Inlet in Glacier Bay National Park, Alaska. This landslide mass has the potential to rapidly move into Tidal Inlet and generate large, long-period-impulse tsunami waves. Field and photographic examination revealed that the landslide moved between 1892 and 1919 after the retreat of the Little Ice Age glaciers from Tidal Inlet in 1890. Global positioning system measurements over a 2-year period show that the perched mass is presently moving at 3–4 cm annually indicating the landslide remains unstable. Numerical simulations of landslide-generated waves suggest that in the western arm of Glacier Bay, wave amplitudes would be greatest near the mouth of Tidal Inlet and slightly decrease with water depth according to Green’s law. As a function of time, wave amplitude would be greatest within approximately 40 min of the landslide entering water, with significant wave activity continuing for potentially several hours.     

降雨诱发区域性滑坡预报预警方法研究

论文以四川省雅安市雨城区为研究区，将逻辑回归模型引入区域降雨型滑坡预警预报，建立了同时考虑降雨强度和降雨过程的降雨临界值表达式。在滑坡危险性区划的基础上，研究提出了区域降雨型滑坡预警预报指标，包括临界值降雨指数R和滑坡发生指数L，并利用20台自动遥测雨量计和地质灾害群测群防网络，采用历史记录雨量和预报雨量，建立了区域降雨型滑坡预报预警体系。     

Geomorphology and landslide susceptibility assessment using GIS and bivariate statistics: a case study in southern Italy

In this article, the results of a study aimed to assess the landslide susceptibility in the Calaggio Torrent basin (Campanian Apennines, southern Italy) are presented. The landslide susceptibility has been assessed using two bivariate-statistics-based methods in a GIS environment. In the first method, widely used in the existing literature, weighting values (Wi) have been calculated for each class of the selected causal factors (lithology, land-use, slope angle and aspect) taking into account the landslide density (detachment zones + landslide body) within each class. In the second method, which is a modification of the first method, only the landslide detachment zone (LDZ) density has been taken into account to calculate the weighting values. This latter method is probably characterized by a major geomorphological coherence. In fact, differently from the landslide bodies, LDZ must necessarily occur in geoenvironmental classes prone to failure. Thus, the calculated Wi seem to be more reliable in estimating the propensity of a given class to generate failure. The thematic maps have been reclassified on the basis of the calculated Wi and then overlaid, with the purpose to produce landslide susceptibility maps. The used methods converge both in indicating that most part of the study area is characterized by a high–very high landslide susceptibility and in the location and extent of the low-susceptible areas. However, an increase of both the high–very high and moderate–high susceptible areas occurs in using the second method. Both the produced susceptibility maps have been compared with the geomorphological map, highlighting an excellent coherence which is higher using method-2. In both methods, the percentage of each susceptibility class affected by landslides increases with the degree of susceptibility, as expected. However, the percentage at issue in the lowest susceptibility class obtained using method-2, even if low, is higher than that obtained using method-1. This suggests that method-2, notwithstanding its major geomorphological coherence, probably still needs further refinements.     

Use of satellite remote sensing data in the mapping of global landslide susceptibility

Satellite remote sensing data has significant potential use in analysis of natural hazards such as landslides. Relying on the recent advances in satellite remote sensing and geographic information system (GIS) techniques, this paper aims to map landslide susceptibility over most of the globe using a GIS-based weighted linear combination method. First, six relevant landslide-controlling factors are derived from geospatial remote sensing data and coded into a GIS system. Next, continuous susceptibility values from low to high are assigned to each of the six factors. Second, a continuous scale of a global landslide susceptibility index is derived using GIS weighted linear combination based on each factor’s relative significance to the process of landslide occurrence (e.g., slope is the most important factor, soil types and soil texture are also primary-level parameters, while elevation, land cover types, and drainage density are secondary in importance). Finally, the continuous index map is further classified into six susceptibility categories. Results show the hot spots of landslide-prone regions include the Pacific Rim, the Himalayas and South Asia, Rocky Mountains, Appalachian Mountains, Alps, and parts of the Middle East and Africa. India, China, Nepal, Japan, the USA, and Peru are shown to have landslide-prone areas. This first-cut global landslide susceptibility map forms a starting point to provide a global view of landslide risks and may be used in conjunction with satellite-based precipitation information to potentially detect areas with significant landslide potential due to heavy rainfall.     

Probabilistic landslide hazard assessment using homogeneous susceptible units (HSU) along a national highway corridor in the northern Himalayas, India

The increased socio-economic significance of landslides has resulted in the application of statistical methods to assess their hazard, particularly at medium scales. These models evaluate where, when and what size landslides are expected. The method presented in this study evaluates the landslide hazard on the basis of homogenous susceptible units (HSU). HSU are derived from a landslide susceptibility map that is a combination of landslide occurrences and geo-environmental factors, using an automated segmentation procedure. To divide the landslide susceptibility map into HSU, we apply a region-growing segmentation algorithm that results in segments with statistically independent spatial probability values. Independence is tested using Moran’s I and a weighted variance method. For each HSU, we obtain the landslide frequency from the multi-temporal data. Temporal and size probabilities are calculated using a Poisson model and an inverse-gamma model, respectively. The methodology is tested in a landslide-prone national highway corridor in the northern Himalayas, India. Our study demonstrates that HSU can replace the commonly used terrain mapping units for combining three probabilities for landslide hazard assessment. A quantitative estimate of landslide hazard is obtained as a joint probability of landslide size, of landslide temporal occurrence for each HSU for different time periods and for different sizes.     

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.     

Advances in landslide nowcasting: evaluation of a global and regional modeling approach

The increasing availability of remotely sensed data offers a new opportunity to address landslide hazard assessment at larger spatial scales. A prototype global satellite-based landslide hazard algorithm has been developed to identify areas that may experience landslide activity. This system combines a calculation of static landslide susceptibility with satellite-derived rainfall estimates and uses a threshold approach to generate a set of ‘nowcasts’ that classify potentially hazardous areas. A recent evaluation of this algorithm framework found that while this tool represents an important first step in larger-scale near real-time landslide hazard assessment efforts, it requires several modifications before it can be fully realized as an operational tool. This study draws upon a prior work’s recommendations to develop a new approach for considering landslide susceptibility and hazard at the regional scale. This case study calculates a regional susceptibility map using remotely sensed and in situ information and a database of landslides triggered by Hurricane Mitch in 1998 over four countries in Central America. The susceptibility map is evaluated with a regional rainfall intensity–duration triggering threshold and results are compared with the global algorithm framework for the same event. Evaluation of this regional system suggests that this empirically based approach provides one plausible way to approach some of the data and resolution issues identified in the global assessment. The presented methodology is straightforward to implement, improves upon the global approach, and allows for results to be transferable between regions. The results also highlight several remaining challenges, including the empirical nature of the algorithm framework and adequate information for algorithm validation. Conclusions suggest that integrating additional triggering factors such as soil moisture may help to improve algorithm performance accuracy. The regional algorithm scenario represents an important step forward in advancing regional and global-scale landslide hazard assessment.     

基于机器学习的滑坡崩塌地质灾害气象风险预警研究

在划分气象风险等级时,传统地质灾害气象风险预警方法忽略了承灾体脆弱性因素,且气象风险预报等级整体偏高,导致高等级风险区空报率较高。基于此,提出基于机器学习的滑坡、崩塌灾害气象风险预警方法。利用信息量法,分析气象因素影响程度。选取坐标点、降雨量、易发生等级,将其作为机器学习人工神经网络的输入节点,判断是否发生崩塌、滑坡灾害;针对地质灾害区域,根据影响程度计算气象引发因子指数,结合滑坡、崩塌灾害潜势度G和承灾体脆弱性M,确定气象风险预警指数R,划分预警级别,完成滑坡、崩塌灾害气象风险预警。实验结果表明,设计方法有效降低了三级预报和四级预警空报率,提升了预警精细化程度。     

Anatexis at 700 to 1000 MPa in the aureole of the Marcy anorthosite, Adirondack Highlands, New York

High temperature (>900 °C) metamorphism affected the New Russia gneiss complex in the aureole of the Marcy anorthosite, Adirondack Highlands, New York. Dehydration melting of pargasitic hornblende and plagioclase in metagabbro during contact metamorphism produced garnet among other phases, an indicator that pressure exceeded 700 MPa during anatexis. Partial melting also occurred in mangerite and charnockite. Minerals that equilibrated during melting yield barometric estimates of 970 ± 100 MPa (garnet–orthopyroxene–plagioclase–quartz in metagabbro and mangerite) and 735 ± 100 and 985 ± 100 MPa (garnet–hornblende–plagioclase–quartz, metagabbro and mangerite, respectively). From these results we infer that the Marcy anorthosite was emplaced at a depth of at least 23 km and probably near 32 km. Received: 9 February 2000 / Accepted: 4 April 2000     

Experimental study of the phase and melting relations of homogeneous basalt + peridotite mixtures and implications for the petrogenesis of flood basalts

Flood basalt provinces may constitute some of the most catastrophic volcanic events in the Earth's history. A popular model to explain them involves adiabatic ascent of plumes of anomalously hot peridotite from a thermal boundary layer deep in the mantle, across the peridotite solidus. However, peridotitic plumes probably require unreasonably high potential temperatures to generate sufficient volumes of magma and high enough melting rates to produce flood volcanism. This lead to the suggestion that low melting eclogitic or pyroxenitic heterogeneities may be present in the source regions of the flood basalts. In order to constrain petrogenetic models for flood basalts generated in this way, an experimental investigation of the melting relations of homogeneous peridotite + oceanic basalt mixtures has been performed. Experiments were conducted at 3.5 GPa on a fertile peridotite (MPY90)–oceanic basalt (GA1) compositional join. The hybrid basalt + peridotite compositions crystallised garnet lherzolite at subsolidus temperatures plus quenched ne-normative picritic liquids at temperatures just above the solidus, over the compositional range MPY90 to GA1₅₀MPY90₅₀. The solidus temperature decreased slightly from ∼1500 °C for MPY90 to ∼1450 °C for GA1₅₀MPY90₅₀. Compositions similar to GA1₃₀MPY90₇₀ have 100% melting compressed into a melting interval which is approximately 50–60% smaller than that for pure MPY90, due to a liquidus minimum. During adiabatic ascent of hybrid source material containing a few tens of percent basalt in peridotite, the lower solidus and compressed solidus–liquidus temperature interval may conspire to substantially enhance melt productivity. Mixtures of recycled oceanic crust and peridotite in mantle plumes may therefore provide a viable source for some flood volcanics. Evidence for this would include higher than normal Fe/Mg values in natural primary liquids, consistent with equilibration with more Fe-rich olivine than normal pyrolitic olivine (i.e. <Fo_89–92). Modelling of fractionation trends in West Greenland picrites is presented to demonstrate that melts parental to the Greenland picrites were in equilibrium at mantle P–T conditions with olivine with Fo_84–86, suggesting an Fe-enriched source compared with normal peridotite, and consistent with the presence of a basaltic component in the source. Received: 29 October 1999 / Accepted: 3 February 2000     

A statistical approach to landslide risk modelling at basin scale: from landslide susceptibility to quantitative risk assessment

A procedure for landslide risk assessment is presented. The underlying hypothesis is that statistical relationships between past landslide occurrences and conditioning variables can be used to develop landslide susceptibility, hazard and risk models. The latter require also data on past damages. Landslides occurred during the last 50 years and subsequent damages were analysed. Landslide susceptibility models were obtained by means of Spatial Data Analysis techniques and independently validated. Scenarios defined on the basis of past landslide frequency and magnitude were used to transform susceptibility into quantitative hazard models. To assess vulnerability, a detailed inventory of exposed elements (infrastructures, buildings, land resources) was carried out. Vulnerability values (0–1) were obtained by comparing damages experienced in the past by each type of element with its actual value. Quantitative risk models, with a monetary meaning, were obtained for each element by integrating landslide hazard and vulnerability models. Landslide risk models showing the expected losses for the next 50 years were thus obtained for the different scenarios. Risk values obtained are not precise predictions of future losses but rather a means to identify areas where damages are likely to be greater and require priority for mitigation actions.     

Study of regional monsoonal effects on landslide hazard zonation in Cameron Highlands, Malaysia

In general, landslides in Malaysia mostly occurred during northeast and southwest periods, two monsoonal systems that bring heavy rain. As the consequence, most landslide occurrences were induced by rainfall. This paper reports the effect of monsoonal-related geospatial data in landslide hazard modeling in Cameron Highlands, Malaysia, using Geographic Information System (GIS). Land surface temperature (LST) data was selected as the monsoonal rainfall footprints on the land surface. Four LST maps were derived from Landsat 7 thermal band acquired at peaks of dry and rainy seasons in 2001. The landslide factors chosen from topography map were slope, slope aspect, curvature, elevation, land use, proximity to road, and river/lake; while from geology map were lithology and proximity to lineament. Landslide characteristics were extracted by crossing between the landslide sites of Cameron Highlands and landslide factors. Using which, the weighting system was derived. Each landslide factors were divided into five subcategories. The highest weight values were assigned to those having the highest number of landslide occurrences. Weighted overlay was used as GIS operator to generate landslide hazard maps. GIS analysis was performed in two modes: (1) static mode, using all factors except LST data; (2) dynamic mode, using all factors including multi-temporal LST data. The effect of addition of LST maps was evaluated. The final landslide hazard maps were divided into five categories: very high risk, high risk, moderate, low risk, and very low risk. From verification process using landslide map, the landslide model can predict back about 13–16% very high risk sites and 70–93% of very high risk and high risk combined together. It was observed however that inclusion of LST maps does not necessarily increase the accuracy of the landslide model to predict landslide sites.