A Comparative Study of Frequency Ratio, Shannon’s Entropy and Analytic Hierarchy Process (AHP) Models for Landslide Susceptibility Assessment

Landslide susceptibility maps are very important tools in the planning and management of landslide prone areas. Qualitative and quantitative methods each have their own advantages and dis-advantages in landslide susceptibility mapping. The aim of this study is to compare three models, i.e., frequency ratio (FR), Shannon’s entropy and analytic hierarchy process (AHP) by implementing them for the preparation of landslide susceptibility maps. Shimla, a district in Himachal Pradesh (H.P.), India was chosen for the study. A landslide inventory containing more than 1500 landslide events was prepared using previous literature, available historical data and a field survey. Out of the total number of landslide events, 30% data was used for training and 70% data was used for testing purpose. The frequency ratio, Shannon’s entropy and AHP models were implemented and three landslide susceptibility maps were prepared for the study area. The final landslide susceptibility maps were validated using a receiver operating characteristic (ROC) curve. The frequency ratio (FR) model yielded the highest accuracy, with 0.925 fitted ROC area, while the accuracy achieved by Shannon’s entropy model was 0.883. Analytic hierarchy process (AHP) yielded the lowest accuracy, with 0.732 fitted ROC area. The results of this study can be used by engineers and planners for better management and mitigation of landslides in the study area.


Introduction
Hilly regions are frequently affected by landslide disasters. Landslides are responsible for huge loss of life and damage to infrastructure [1,2]. A landslide is the result of movement of a slope under the effect of gravity, but the movement can be triggered by the factors such as rainfall, earthquakes, erosion, etc. [3,4]. The accurate prediction of landslides can help in the planning and management of landslide hazards, and can ultimately be used for the reduction of risk [5]. Landslide susceptibility maps help in the identification of landslide prone areas. Landslide susceptibility maps can be very efficient tools for planners and risk managers [5]. The occurrence of landslides is a complex phenomenon which depends upon various factors. Geological factors, drainage characteristics, land-use of the region, construction activities, etc., can all be responsible for the occurrence of landslides [6]. Anthropogenic activities and the development of infrastructure in landslide prone areas can disturb natural slopes, resulting in massive disaster [7,8].
The study area selected in this study is the Shimla district in Himachal Pradesh (H.P.), India. The study area is developing at a very rapid rate. Anthropogenic activities have been increasing in the study area, resulting in the instability of natural slopes. Landslide occurrence is very frequent in the study area due to rugged topographic and the typical climatic conditions in the region. Therefore, there is a need for the preparation of a landslide susceptibility map for the region which can be used by planners and engineers. 2 of 20 In different regions of the world, landslide susceptibility mapping has been undertaken by various researchers, using geographic information systems (GIS) and remote sensing [2,4,6,[9][10][11][12][13][14]. Qualitative and quantitative approaches are both widely used. Qualitative techniques consider expert opinion, while quantitative approaches establish a mathematical relationship between causative factors and the occurrence of landslide [2,5,13,[15][16][17]. Landslide inventory-based methods utilize a direct approach in which susceptibility maps are prepared by considering the geomorphological characteristics of the region. Ranking and weightage are introduced to improve direct susceptibility mapping, through analytic hierarchy processes (AHP), analytic network processes (ANP), fuzzy logic, etc. [15][16][17][18]. These methods can be categorized as semi-qualitative methods. Quantitative methods such as the bivariate method [19][20][21], multi-variate method [22][23][24], weight of evidence [25,26], etc., are widely used with the development of GIS techniques. GIS platforms help in the calculation and visualization of the cumulative effects of causative factors on landslides. Qualitative techniques are based on the subjective judgment of experts, while the output in quantitative techniques depends on the availability of historical landslide data [27].
Landslide susceptibility mapping has not been undertaken for the study area until now. There is a need for the preparation of an accurate landslide susceptibility map for the region. In this study, we used three techniques for landslide susceptibility mapping, i.e., frequency ratio (FR), Shannon's entropy model and analytic hierarchy process (AHP). Shannon's entropy model is an improvement over the frequency ratio model and is rarely used in the Indian continent for landslide susceptibility mapping. The results of this study can help in the delineation of landslide prone regions in the study area. Construction planners, environmental engineers and risk managers can use the results of this study for the planning and mitigation of landslide hazards in the study area. Figure 1a shows the methodology adopted in this study. In the first step, the study area is selected. Shimla is the capital of Himachal Pradesh (H.P.) province in India. The region lies between longitudes 77 • 0 E and 78 • 19 and latitudes 30 • 45 and 31 • 44 . The elevation of the region varies from 300 m to 6000 m and shows variation in the vegetation [28,29]. The study area lies in the Kashmir and Western Himalayas seismic region. Shimla district lies in a high damage risk zone in the case of seismic activity (MSK VIII) [28,29]. Major earthquakes that have occurred in Himachal Pradesh (H.P.) are the Kangra earthquake 1905 (M = 8), the Chamba earthquake 1945 (M = 6.5), the Kinnaur earthquake (M = 6.8), etc. The study area is predominantly covered by rock from the Jutogh, Shali, Simla and Rampur groups. The study area consists of lithological units of shale, siltstones, quartzite, dolomites, phyllite, schist, conglomerate, etc. The region is predominantly covered by the Proterozoic age group, Plaeo-proterozoic age group and Neo-proterozoic age group. The study area is mainly forest and cultivated land. Some of the study area in upper Shivaliks consists of hill tops with snow. The temperature varies from 4 • C to 31 • C [28,29]. Therefore, the region is cool in winter and moderately warm in summer. Figure 1b shows the study area and landslide inventory.

Landslide Inventory
A landslide inventory plays an important role in the mathematical modelling and validation of the final output [13]. Landslide inventories consist of information regarding the location of landslides, area of landslides, activity and date of landslides, etc. [30]. A landslide inventory can be simple to complex, depending upon the requirement of study. A landslide inventory is prepared from previous literature, Google Earth and field visits. The field visits cannot be undertaken frequently due to restricted movement in the current COVID-19 pandemic. The landslide inventory consists of more than 1500 landslides in the study area. The landslides are represented by polygons. Figure 1 shows the landslide inventory.

Landslide Inventory
A landslide inventory plays an important role in the mathematical modelling and validation of the final output [13]. Landslide inventories consist of information regarding the location of landslides, area of landslides, activity and date of landslides, etc. [30]. A landslide inventory can be simple to complex, depending upon the requirement of study. A landslide inventory is prepared from previous literature, Google Earth and field visits. The field visits cannot be undertaken frequently due to restricted movement in the current COVID-19 pandemic. The landslide inventory consists of more than 1500 landslides in the study area. The landslides are represented by polygons. Figure 1 shows the landslide inventory.

Landslide Conditioning Factors
Understanding the causative factors of landslides for the preparation of landslide susceptibility maps is very important. The data for the preparation of thematic layers of causative factors were obtained from different sources. An ASTER digital elevation model (DEM) was obtained from the USGS website. The resolution of the DEM used in this study is 1 arc second, which is approximately equal to 30 × 30 meters. The average route mean square error (RMSE) in the measurement of elevation is ±13.31 m. In terms of absolute positioning, the ASTER DEM models have a tendency for displacement by 5.3 m in our study. Landsat 8 data was used for the preparation of a land-use map for the study area. The resolution for Landsat bands 1-7 is 30 × 30 meters. Hydrogeological and topograph-

Landslide Conditioning Factors
Understanding the causative factors of landslides for the preparation of landslide susceptibility maps is very important. The data for the preparation of thematic layers of causative factors were obtained from different sources. An ASTER digital elevation model (DEM) was obtained from the USGS website. The resolution of the DEM used in this study is 1 arc second, which is approximately equal to 30 × 30 m. The average route mean square error (RMSE) in the measurement of elevation is ±13.31 m. In terms of absolute positioning, the ASTER DEM models have a tendency for displacement by 5.3 m in our study. Landsat 8 data was used for the preparation of a land-use map for the study area. The resolution for Landsat bands 1-7 is 30 × 30 m. Hydrogeological and topographical maps were used for the digitization of faults, lithology and the road network. The land cover of the study area was obtained from Landsat 8 data.
Selection of condition factors depends on the availability of data, ease of use and the results of a literature survey. In many studies, DEM based factors are considered to be major causative factors, but in recent studies, land use, drainage networks and faults or lineaments are also considered to be responsible for landslide occurrence. The factors considered in this study are explained as follows. The relative importance of the causative factors can be found using Shannon's entropy and the AHP model.

Slope Gradient
Slope gradient may be defined as the rate of change in elevation with respect to horizontal surface. When the shear stresses in a downward direction exceed the shear strength of soil or rock mass, the failure of slopes occur. Therefore, the slope angle directly affects the stability of slope [31]. The slope of the study area was extracted from a digital elevation model (DEM) of the study area. The study area is comprised predominantly of steep slopes greater than 60 • . As the slope angle increases, the possibility of slope failure increases. Figure 2a shows slope gradient of the study area.
Selection of condition factors depends on the availability of data, ease of use and the results of a literature survey. In many studies, DEM based factors are considered to be major causative factors, but in recent studies, land use, drainage networks and faults or lineaments are also considered to be responsible for landslide occurrence. The factors considered in this study are explained as follows. The relative importance of the causative factors can be found using Shannon's entropy and the AHP model.

Slope Gradient
Slope gradient may be defined as the rate of change in elevation with respect to horizontal surface. When the shear stresses in a downward direction exceed the shear strength of soil or rock mass, the failure of slopes occur. Therefore, the slope angle directly affects the stability of slope [31]. The slope of the study area was extracted from a digital elevation model (DEM) of the study area. The study area is comprised predominantly of steep slopes greater than 60°. As the slope angle increases, the possibility of slope failure increases. Figure 2a shows slope gradient of the study area.

Slope Aspect
Slope aspect gives the direction of maximum slope with respect to magnetic north. Slope aspect affects erosion, evapotranspiration, desertification, solar heating and soil moisture [5]. Slope aspect indirectly shows the impact of water on the mass. Figure 2b shows the slope aspect of the study area.

Slope Aspect
Slope aspect gives the direction of maximum slope with respect to magnetic north. Slope aspect affects erosion, evapotranspiration, desertification, solar heating and soil moisture [5]. Slope aspect indirectly shows the impact of water on the mass. Figure 2b shows the slope aspect of the study area.

Relative Relief
Relative relief shows the change in elevation. The change in elevation indicates variation in vegetation type and topographical characteristics of the region [13]. Vegetation roots hold the rock and soil mass together and make it more stable. The relative relief varies from 0 to 254 m in the present study. Figure 2c shows the relative relief of the region.

Topographic Wetness Index (TWI)
The topographic wetness index (TWI) is the tendency of terrain to accumulate water. It indicates variation in soil moisture [32]. Figure 2d shows the TWI of the region. If SCA is specific catchment area and ϕ is the angle of slope, TWI is given by the following formula:

Lithology
Lithology is an important causative factor of landslides. Variation in rock characteristics and age prominently affect landslide events [33]. Detailed information about the lithological and geological character of the region is given in Table 1. Figure 3a shows the lithology of the region. The drainage characteristics of the region significantly affect the occurrence of landslides. River water erodes the soil at the toe of the slope and weakens the base, which may result in slope failure [33,34]. Drainage density is divided into four categories for the study area. Higher drainage density increases the risk of landslides. Figure 3b shows the drainage density of study area.

Distance from Roads
Road construction disturbs the natural slopes along highway corridors. Road construction activities involve significant cutting and filling activities. The bottom section of a slope is weakened by the soil and rock cutting activities during highway construction [35]. Vehicular vibrations, excessive excavations and disturbed natural slopes all induce instability and increase the risk of landslide occurrence. Slopes existing near roads are more prone to experience landslides [13,35]. A detailed analysis of slope stability during the planning of construction of projects can reduce the risk of landslide occurrence. Distance from roads indirectly represents the effects of human interference on natural slopes. If proper slope protection is provided and the slopes along highways are planned properly, the risk of landslides along highways can be reduced. Figure 3c shows the distance from roads.

Drainage Density
The drainage characteristics of the region significantly affect the occurrence of landslides. River water erodes the soil at the toe of the slope and weakens the base, which may result in slope failure [33,34]. Drainage density is divided into four categories for the study area. Higher drainage density increases the risk of landslides. Figure 3b shows the drainage density of study area.

Distance from Roads
Road construction disturbs the natural slopes along highway corridors. Road construction activities involve significant cutting and filling activities. The bottom section of a slope is weakened by the soil and rock cutting activities during highway construction [35]. Vehicular vibrations, excessive excavations and disturbed natural slopes all induce instability and increase the risk of landslide occurrence. Slopes existing near roads are more prone to experience landslides [13,35]. A detailed analysis of slope stability during the planning of construction of projects can reduce the risk of landslide occurrence. Distance from roads indirectly represents the effects of human interference on natural slopes. If proper slope protection is provided and the slopes along highways are planned properly, the risk of landslides along highways can be reduced. Figure 3c shows the distance from roads.

Distance from Faults
Cracks, faults, fissures and lineament can increase pore pressure in the rock and soil mass, and can induce instability in slopes [36,37]. In this study, faults were digitized from

Distance from Faults
Cracks, faults, fissures and lineament can increase pore pressure in the rock and soil mass, and can induce instability in slopes [36,37]. In this study, faults were digitized from NRSA data. Proximity to faults is considered a causative factor in this study, and Euclidean distance from faults was calculated in ArcGIS. Figure 3d shows distance from faults.

Land Cover
Land cover is an important factor that affects the occurrence of landslides. Vegetation area, forest area, barren land, etc. are shown in a land cover map, and may affect the occurrence of landslides due to their diverse characteristics. A land cover map was extracted from Landsat 8 imagery. The land use was classified using a supervised learning technique of image classification in ArcGIS. The study area is predominantly covered with the cropland and forests. In the upper elevations, there are glaciers and snow covered hill tops as well. Figure 4 shows the land cover map of the region. area, forest area, barren land, etc. are shown in a land cover map, and may affect the occurrence of landslides due to their diverse characteristics. A land cover map was extracted from Landsat 8 imagery. The land use was classified using a supervised learning technique of image classification in ArcGIS. The study area is predominantly covered with the cropland and forests. In the upper elevations, there are glaciers and snow covered hill tops as well. Figure 4 shows the land cover map of the region.

Frequency Ratio
Frequency ratio models are based upon the relationship between the occurrence of landslide events and the causative factors [38,39]. A relationship between the landslide inventory and causative factors can be established in a GIS environment. A landslide will occur under the same conditions as past landslides [5]. The frequency ratio is the ratio between the percent areas where landslides have occurred in a class to the percent area of the influencing class relative to the whole study area as shown in the following formula [40].

Frequency Ratio
Frequency ratio models are based upon the relationship between the occurrence of landslide events and the causative factors [38,39]. A relationship between the landslide inventory and causative factors can be established in a GIS environment. A landslide will occur under the same conditions as past landslides [5]. The frequency ratio is the ratio between the percent areas where landslides have occurred in a class to the percent area of the influencing class relative to the whole study area as shown in the following formula [40].
where Fr is frequency ratio and n is number of causative factors. Table 2 shows the frequency ratio (FR) calculation.

Shannon's Entropy
Shannon's entropy model is an improvement on the frequency ratio model. The frequency ratio model only considers the weightage of sub-factors, and does not consider the weightage of causative factors. Shannon's entropy measures the uncertainty or instability of a system [41]. With landslide susceptibility mapping, it measures the influence of causative factors on the occurrence of landslides [42][43][44][45]. Weightage of causative factors is calculated using the following procedure: where Pij is the probability density and FR is the frequency ratio of sub-factors. Wj is the weightage of causative factors obtained from Shannon's entropy technique. These values can then be used for assigning weight to causative factors, while frequency ratio values are used for sub-factors. Table 2 shows Shannon's entropy weightage of causative factors.

Analytic Hierarchy Process (AHP)
Analytic hierarchy process (AHP) is a decision-making tool which helps to solve complex problems with simple criteria. AHP is based upon three principles i.e., decomposition of the problem, comparative judgment, and synthesis of relative importance or rankings [43]. In AHP, the problem is broken down into hierarchical criteria. These criteria are compared to each other. This process of relative comparison is called pair-wise com-parison. The eigen vector method is used to calculate the rankings and the consistency of the solution is then also checked by calculating the consistency ratio [43]. The consistency of the weights for relative importance assigned during the pairwise comparison can be checked for consistency using the equation given below.
Consistency Ratio (CR) = CI RI (8) where, CI is consistency index while RI is randomness index. CI is calculated as follows: where, λmax = Major eigen value and n = order of matrix.
Randomness index values are given by Saaty, which is dependent on the value of n. RI is the result of extensive experimentation on large dataset samples. If the CR values are less than 10%, the pairwise comparison is considered consistent. If the CR value is more than 10%, the solution is considered inconsistent and weights are reassigned in a pairwise comparison matrix. Tables 3 and 4 show the relative comparison matrix of AHP for causative factors and sub-factors.

Frequency Ratio
The frequency ratio (FR) of the causative sub-factors is shown in the Table 2. The relative importance of the sub-factors can be evaluated from the weightage calculated in the Table 2. The impact of slopes less than 45 • is very low, but as the slope increases above 45 • , ISPRS Int. J. Geo-Inf. 2021, 10, 603 13 of 20 the frequency ratio increases sharply. More than 39% of landslides occur in regions where the slope is more than 75 • , while such slopes only cover 18% of the study area. Higher slope values trigger the effect of gravity and also increase shear stress. The maximum frequency ratio (FR) value for slopes is 2.163. East, south-east, south, and North Slope aspects affect the occurrence of landslides more compared other directions. The frequency ratio (FR) is highest for the east direction (1.673), followed by the south east direction (1.486). These aspects receive more rainfall and are also subjected to erosion, which makes them more prone to landslides. The areas of low (0-50 m) and high (100-150 m) relative relief were more prone to landslides. The frequency ratio value for relative relief was highest (1.262) for 0-50 m relief, followed by 1.073 for 100-150 m relief. Lower values of topographic wetness index (TWI) have a significant effect on landslide occurrence, with a frequency ratio of 1.043, while very high values of TWI (FR = 0.037) had minimal effect on landslide occurrence. As TWI increases, FR decreases. Neoproterozoic age deposits had the highest FR values (1.252), followed by Mesoproterozoic deposits (FR = 1.653). The frequency ratios follow an increasing trend as the drainage density increases. For very high drainage density, the FR value was 1.685. Smaller distance from roads had a significant effect on the occurrence of landslides. Small distances from roads have the highest FR value i.e., 1.253, and similarly, smaller distances from faults also showed maximum FR (0.994) values in its category. Cutting activities performed during road construction reduced strength at the toe of slopes, making them more prone to landslides. Snow covered land (FR = 2.059) was more prone to landslides, while agricultural land (FR = 0.562) had the least effect on landslide occurrence. After snow cover, barren land due to less vegetative cover was found to be more prone to landslides, with an FR value of 1.420.
The landslide susceptibility index (LSI) for the frequency ratio model varied from 3.173 to 13.297. The landslide susceptibility map is shown in Figure 5. The map is divided using natural breaks in ArcGIS showing very low, low, moderate, high and very high susceptibility zones. It is observed that high and very high landslide susceptibility zones cover around 45% of the study area, but more than 70% landslides occur in these regions. The very high susceptibility zone covers only 16.62% of the total study area but it is subjected to 40.7% of the total landslide area. 11% of the total landslide area occur in very low and low susceptibility zones while these zones cover more than 25% area.

Shannon's Entropy Model
The weightage of sub-factors in Shannon's entropy model was based on frequency ratio (FR) values. The weightage of causative factors was evaluated from FR values of sub-factors. It is found that distance from faults and distance from roads are the major causative factors which have a very high impact on the occurrence of landslides. The weightage of distance from faults is 0.159 and weightage of vicinity to roads is found to be 0.148. The anthropogenic activities of road cutting and development have significantly increased landslide activities in the study area. The topographic wetness index (TWI) is also an important factor, with a weightage of 0.135 as per Shannon's entropy model. Slope has a weightage of 0.120. Lithology, drainage density, and relative relief have almost equal weightage. The very high susceptibility zone covers only 16.62% of the total study area but it is subjected to 40.7% of the total landslide area. 11% of the total landslide area occur in very low and low susceptibility zones while these zones cover more than 25% area.

Shannon's Entropy Model
The weightage of sub-factors in Shannon's entropy model was based on frequency ratio (FR) values. The weightage of causative factors was evaluated from FR values of subfactors. It is found that distance from faults and distance from roads are the major causative factors which have a very high impact on the occurrence of landslides. The weightage of distance from faults is 0.159 and weightage of vicinity to roads is found to be 0.148. The anthropogenic activities of road cutting and development have significantly increased landslide activities in the study area. The topographic wetness index (TWI) is also an important factor, with a weightage of 0.135 as per Shannon's entropy model. Slope has a weightage of 0.120. Lithology, drainage density, and relative relief have almost equal weightage.
The results of the FR model have been changed a little following the implementation of Shannon's entropy weightage to major causative factors. The simplest landslide susceptibility equation for this model is given as follows: The results of the FR model have been changed a little following the implementation of Shannon's entropy weightage to major causative factors. The simplest landslide susceptibility equation for this model is given as follows: LSI varies from 0.402 to 1.412 for Shannon's entropy model. The final landslide susceptibility map is divided into five categories, i.e., very low, low, moderate, high and very high, using natural breaks. The landslide susceptibility map by Shannon's entropy is shown in the Figure 6. Shannon's entropy model shows that 12.06% of the area has very high susceptibility, which contains 31.93% landslide area of the region. The model shows that 23.66% of the area has high susceptibility, which contains 34.45% landslide area. Very low and low landslide susceptibility zones cover 9.84% and 24.85% of study area, respectively.
shown in the Figure 6. Shannon's entropy model shows that 12.06% of the area has high susceptibility, which contains 31.93% landslide area of the region. The model s that 23.66% of the area has high susceptibility, which contains 34.45% landslide area. low and low landslide susceptibility zones cover 9.84% and 24.85% of study area, re tively.

Analytic Hierarchy Process (AHP) Approach
According to the AHP model, distance from roads and drainage density are the m factors which have high weightage, i.e., 0.280 and 0.205 respectively. Higher weighta distance from roads shows the significance of anthropogenic activities towards occur of landslides. Lithology has a weightage of 0.096 in the AHP model. Neoproterozoi mesoproterozoic deposits have been more responsible for occurrence of landslides pared to other lithological units, as per the AHP model. Distance from faults has a we age of 0.161, which states the high significance of this factor in the occurrence of slides, while aspect has the lowest weightage, i.e., 0.016. Slope (weightage = 0.105 (weightage = 0.081) are important DEM-based factors which affect the landslide phe enon. Lower TWI values are more responsible for the occurrence of landslides comp

Analytic Hierarchy Process (AHP) Approach
According to the AHP model, distance from roads and drainage density are the major factors which have high weightage, i.e., 0.280 and 0.205 respectively. Higher weightage of distance from roads shows the significance of anthropogenic activities towards occurrence of landslides. Lithology has a weightage of 0.096 in the AHP model. Neoproterozoic and mesoproterozoic deposits have been more responsible for occurrence of landslides compared to other lithological units, as per the AHP model. Distance from faults has a weightage of 0.161, which states the high significance of this factor in the occurrence of landslides, while aspect has the lowest weightage, i.e., 0.016. Slope (weightage = 0.105) and (weightage = 0.081) are important DEM-based factors which affect the landslide phenomenon. Lower TWI values are more responsible for the occurrence of landslides compared to higher values. Tables 3 and 4 show the AHP weightage for all the factors and sub-factors, respectively.
The landslide susceptibility index (LSI) for the AHP approach is given by following equation: The landslide susceptibility index (LSI) for the landslide susceptibility map obtained from the AHP approach varies from 0.061 to 0.648. The landslide susceptibility map is divided by natural breaks based on the landslide susceptibility index shown in Figure 7. ( ) = 0.105 × + 0.016 × + 0.036 × + 0.081 × + 0.096 × ℎ + 0.205 × + 0.280 × + 0.168 × + 0.020 × The landslide susceptibility index (LSI) for the landslide susceptibility map obtained from the AHP approach varies from 0.061 to 0.648. The landslide susceptibility map is divided by natural breaks based on the landslide susceptibility index shown in Figure 7. It is found that 23.05% area lies under very low landslide susceptibility according to the AHP approach, and contains around 10% landslide area. Around 24% area lies under high susceptibility and contains around 29% landslide area, while around 8% area lies under very high susceptibility containing 9.24% of the landslides. It seems evident that the AHP model remains less efficient in the classification of landslide susceptibility as compared to the previous two models. The relative importance of the causative factors has been selected based on the subjective approach of experts. There is no mathematical relationship for the comparison of the relative importance of causative factors and sub- It is found that 23.05% area lies under very low landslide susceptibility according to the AHP approach, and contains around 10% landslide area. Around 24% area lies under high susceptibility and contains around 29% landslide area, while around 8% area lies under very high susceptibility containing 9.24% of the landslides. It seems evident that the AHP model remains less efficient in the classification of landslide susceptibility as compared to the previous two models. The relative importance of the causative factors has been selected based on the subjective approach of experts. There is no mathematical relationship for the comparison of the relative importance of causative factors and subfactors. In qualitative techniques such as AHP, the input can be varied for better output. The lower accuracy in the AHP may be due to underestimation and overestimation of the impacts of some causative factors and sub-factors.

Validation of Results
The results were validated using a receiver operating curve (ROC). Figure 8 shows the ROCs for the landslide susceptibility maps prepared with different techniques. The ROC curves show a relationship between sensitivity and (1-specificity). The area under curve (AUC) value for the frequency ratio model is 0.925, which shows very high accuracy. The frequency ratio model shows very encouraging results in the study area. Shannon's entropy model shows an accuracy of 0.883. The analytic hierarchy process (AHP) based landslide susceptibility map shows an AUC value of 0.732, which is the lowest among all the models.

Validation of Results
The results were validated using a receiver operating curve (ROC). Figure 8 sho the ROCs for the landslide susceptibility maps prepared with different techniques. T ROC curves show a relationship between sensitivity and (1-specificity). The area und curve (AUC) value for the frequency ratio model is 0.925, which shows very high ac racy. The frequency ratio model shows very encouraging results in the study area. Sh non's entropy model shows an accuracy of 0.883. The analytic hierarchy process (AH based landslide susceptibility map shows an AUC value of 0.732, which is the low among all the models.

Conclusions and Future Scope
Landslides are one of the most disastrous phenomena in the Himalayan region. this study, statistical approaches and expert based approaches are compared for landsl susceptibility mapping in a GIS environment. The frequency ratio (FR) model shows b ter accuracy i.e., 0.925 compared to Shannon's entropy model and the AHP. Shanno entropy model had an accuracy of 0.883, and AHP model had an accuracy of 0.732 us the ROC technique. The FR model and Shannon's entropy model have inputs based landslide inventories which cannot be varied, while in the AHP model, the input can

Conclusions and Future Scope
Landslides are one of the most disastrous phenomena in the Himalayan region. In this study, statistical approaches and expert based approaches are compared for landslide susceptibility mapping in a GIS environment. The frequency ratio (FR) model shows better accuracy i.e., 0.925 compared to Shannon's entropy model and the AHP. Shannon's entropy model had an accuracy of 0.883, and AHP model had an accuracy of 0.732 using the ROC technique. The FR model and Shannon's entropy model have inputs based on landslide inventories which cannot be varied, while in the AHP model, the input can be varied to improve accuracy. The relative importance of each factor and sub-factor can be observed from these models. It can be concluded that frequency ratio is the most efficient method, and is the simplest to implement in quantitative methods. An accurate and extensive landslide inventory is required for the efficient implementation of frequency ratio models. The results from qualitative methods like AHP depend upon the selection of weights, so need additional research. AHP can be used in the regions for which historical landslide data is not available. There is evidence of rock fall, debris flows and landslides in high and very high susceptibility zones. Landslide susceptibility maps using these models have not previously been prepared in this study area. Therefore, the results of this study can provide a very useful input for the planner and risk managers.
The study could be extended to predict landslide hazard and risk assessment of the study area. The seasonal variation in causative factors such as rainfall, land-use variation and vegetation cover variation can be studied and their impact on the occurrence of landslides can be observed.