Landslide susceptibility evaluation and hazard zonation techniques – a review

Landslides are the most destructive geological hazard in the hilly regions. For systematic landslide mitigation and management, landslide evaluation and hazard zonation is required. Over the past few decades several techniques have been developed that can be used for landslide evaluation and zonation. These techniques can broadly be classified into qualitative and quantitative approaches. Qualitative approaches include geomorphological analysis and heuristic techniques whereas quantitative approaches include statistical, artificial intelligence and deterministic techniques. In quantitative techniques prediction for landslide susceptibility is based on the actual realistic data and interpretations. Further, the quantitative techniques also overcome the subjectivity of qualitative approaches. Each of these techniques may consider different causative factors and utilizes various means for factor evaluation and analysis. When compared, each of these techniques has its own advantage and disadvantage over other techniques. The selection of appropriate technique for landslide hazard evaluation and zonation is very crucial. The factors that need to be considered to adopt an appropriate approach are; investigation purpose, the extent of the area to be covered, the type of mapping units, the scale of map to be produced, type of data to be used, type of landslides, availability of resources, capability and skill set of an evaluator and the accessibility to the study area. The main aim of this article is to present a comprehensive review on various techniques and approaches available for landslide susceptibility and hazard zonation mapping. Further, attempt is also made to assess the effectiveness of these techniques in landslide hazard zonation studies.


Introduction
Landslides are considered to be the most damaging geological hazard in mountainous regions (Mengistu et al. 2019;Hamza and Raghuvanshi 2017;Girma et al. 2015;Raghuvanshi et al. 2014a;Pan et al. 2008;Kanungo et al. 2006;Crozier and Glade 2005;Parise and Jibson 2000;Varnes 1996). In the landslide study there are three basic components, these are landslide susceptibility, landslide hazard and landslide risk. Landslide susceptibility mapping or zonation is the subdivision of the terrain in to zones that have a different likelihood for landslide occurrence. It includes spatial distribution, size, location and displacement of the landslide deposit (Fell et al. 2008;Guzzetti et al. 1999;Varnes 1984). Landslide hazard mapping/zonation is a division of terrain into zones that are basically characterized by the spatial and temporal probability of landslide occurrence that includes description of location, volume, and prediction of future landslide occurrence in an area (Fell et al. 2008;Fell et al. 2007;AGS 2000).
In order to minimize the damage due to landslides it is necessary to evaluate the factors that are responsible for the landslides. These factors are related to geology, geomorphology, landuse and landcover, rainfall, seismicity, manmade activities etc. (Raghuvanshi et al. 2014a;Anbalagan 1992). For landslide studies it is generally assumed that the combination of these factors may possibly lead to landslides in a given area. Therefore, evaluation of these factors and their relation with the past landslides in an area may form the basis for the prediction of future landslides (Chimidi et al. 2017;Negassa and Kala 2015;Girma et al. 2015).
In order to evaluate and zone the area for potential landslides, landslide hazard evaluation and zonation techniques may be applied (Hamza and Raghuvanshi 2017;Girma et al. 2015;Raghuvanshi et al. 2014a;Bisson et al. 2014;Pan et al. 2008;Anbalagan 1992). In past several techniques have been developed by various researchers that can broadly be classified into inventory based approaches, expert evaluation (Raghuvanshi et al. 2014a;Guzzetti et al. 1999;Turrini and Visintainer 1998;Sarkar et al. 1995;Anbalagan 1992;Pachauri and Pant 1992 etc.), statistical (Negassa and Kala 2015; Girma et al. 2015;Kanungo et al. 2006;Carrara et al. 1992 etc.), deterministic (Fall et al. 2006), probabilistic and distribution free approaches (Kanungo et al. 2006). Several studies have attempted to list and review the techniques those are applied for landslide susceptibility and hazard zonation (Lee 2015;Thiebes et al. 2012;Michoud et al. 2012;Jaboyedoff et al. 2012;Jongmans and Garambois 2007;Angeli et al. 2000;Turner and Schuster 1996, etc.). Each of these techniques considers various factors and uses different means for factor evaluation and analysis. When compared, each of these techniques have its own characteristics and may have certain advantage and disadvantage over the other techniques (Raghuvanshi et al. 2014a;Fall et al. 2006;Kanungo et al. 2006;Casagli et al. 2004;Guzzetti et al. 1999;Leroi 1997). Further, some degree of uncertainty always exists in these techniques owing to factors that are considered or the methods that are followed to generate the factor data (Negassa and Kala 2015; Carrara et al. 1995).
Further, these landslide hazard techniques can be applied at different scales depending on the area to be covered, the geological and geo-morphological factors to be considered, the methods by which the concerned factor data will be generated and the capability and skill set of an evaluator (Ermias et al. 2017;Chimidi et al. 2017;Fall et al. 2006;Kanungo et al. 2006;Casagli et al. 2004;Guzzetti et al. 1999;Leroi 1997;Carrara et al. 1992). The main objective of the present research work is to present a comprehensive review on various techniques and approaches that are available for landslide susceptibility and hazard zonation mapping. Further, it is also attempted to assess the effectiveness of these techniques in landslide hazard zonation studies.

Data types and software for landslide studies
Data collected from the field visit and/or remote sensing image analyses are the initial steps for landslide studies. The technique, by which data pertaining to landslide study is generated, primarily depends on the scale of the study, purpose of the study and the accessibility to the study area. The data that is needed for landslide susceptibility and hazard zonation studies mainly includes data on; landslide inventory, environmental and triggering factors (Sreedevi and Yarrakula 2016;Metternicht et al. 2005;Soeters and Van Westen 1996). The environmental factors mainly consist of slope angle, slope aspect, flow accumulation, lithology, geological structures/ proximity to faults, soil type, soil depth/ thickness, slope hydrology, geomorphology, proximity to streams, road cuts, land use types and changes etc. The triggering factors include daily rainfall, rainfall intensity/ duration, earthquake and volcanic eruption. In general, it is not possible to define a uniform list of such causal factors that would be required for the landslide studies. The selection of such causal factors may differ, depending on the scale of analysis, the characteristics of the study area, the landslide type, and on the failure mechanisms. The main difference between the data collection for landslide susceptibility and hazard zonation is in the collection of landslide inventory data; for landslide hazard spatial and temporal landslide data from archive, image analysis etc. is required whereas for landslide susceptibility only spatial data is sufficient. Further, based on the type of the selected technique, the data listed above is changed either in to grid cell, terrain units, uniquecondition units, slope units or in to topographic units (Reichenbach et al. 2018;Guzzetti et al. 1999).
The commonly used software for landslide studies are ArcGIS package that is used for mapping landslide influenced area, factor maps preparation, overlay analysis and interpretations. Erdas Imagine is used for image processing that allows processing for geospatial factor layers used for landslide susceptibility and hazard studies. Further, Geostudio, Geoslope's and SLOPE/W are used for geotechnical slope stability studies. Also, SPSS/real statistics and MS Excel are used for data management and statistical analysis.

Methods of landslide susceptibility and hazard zonation
In order to evaluate landslide susceptibility and hazard zonation several approaches and techniques have been proposed. These includes landslide inventories, heuristic terrain and susceptibility zoning, statistical methods, deterministic methods, probabilistic methods etc. (Reichenbach et al. 2018;Corominas et al. 2014;Raghuvanshi et al. 2014a;Negassa and Kala 2015;Kanungo et al. 2006;Aleotti and Chowdhury 1999;Guzzetti et al. 1999). Broadly, all these techniques or approaches may be further classified into qualitative and quantitative approaches (Fig. 1). The qualitative approaches include distribution analysis or inventory, geomorphic analysis and the expert (heuristic) evaluation techniques which are based on the knowledge and experience of the evaluator (Corominas et al. 2014;Raghuvanshi et al. 2014a). The quantitative approaches mainly include statistical, deterministic, probabilistic and distribution free techniques (Raghuvanshi et al. 2014a;Kanungo et al. 2006). The qualitative approaches are considered to be subjective in nature (Girma et al. 2015;Raghuvanshi et al. 2014a;Kanungo et al. 2006;Fall et al. 2006;Casagli et al. 2004) whereas quantitative approaches are considered to be objective in nature (Girma et al. 2015;Fall et al. 2006). In recent times quantitative approaches have been widely used for landslide susceptibility and hazard evaluation.

Direct approaches Geomorphic and landslide inventory techniques
The geomorphic techniques are the direct methods of landslide susceptibility and hazard zonation mapping. Most of the data pertaining to various factors responsible for landslides are obtained through remote sensing and the data procurement from direct field work is rarely employed in these methods. The geomorphic mapping of landslides susceptibility is carried out through aerial photographs (Mandaglio et al. 2016;Kanungo et al. 2009;Guzzetti et al. 1999;Verstappen 1983) or from satellite images (Nossin 1989). These methods most commonly, but not necessarily are associated with the production of a landslide inventory map. Landslide inventory includes mapping of the past landslides in the area. This method also includes collection and recording of data on location, type and dimensions of landslides. Beside, data/ information on responsible causative factors, triggering mechanism and landslides occurrence frequency is also collected (Girma et al. 2015;Fall et al. 2006;. Landslide inventory mapping is considered to be straight forward and is required for most of the susceptibility and hazard zonation techniques  either to formulate the general rules for the hazard prediction or to validate the predicted model. Generally, these geomorphic techniques do not consider the factors that have initiated or triggered the landslides. The principal factors that are considered for landslide hazard zonation are lithology, geomorphology, landuse and landcover and presence or absence of landslides in the area.

Indirect approaches Expert evaluation
These techniques are indirect method of landslide evaluation and zonation. These techniques are based on various causative factors that are selected by the evaluator. In these techniques expert makes a decision on the degree and type of hazard for the given area, for which he may use either an indirect or a direct mapping approach. The landslide hazard is evaluated based on quasi-static variables (Fall et al. 2006; which are primarily based on the judgment and experience of an evaluator (Girma et al. 2015;Raghuvanshi et al. 2014a;Fell et al. 2008;Gorsevski et al. 2003;Guzzetti et al. 1999;Aleotti and Chowdhury 1999;Van Westen et al. 1997;Soeters and Van Westen 1996;Varnes and IAEG 1984;Carrara et al. 1995;Hutchinson 1995). Since these methods are based on the evaluation of an expert therefore they are considered to be subjective in nature (Girma et al. 2015;Ruff and Czurda 2008;Kanungo et al. 2006;Fall et al. 2006;Casagli et al. 2004).
The landslide susceptibility and hazard evaluated by heuristic technique for a given area may vary considerably if evaluated by different experts. However, these techniques are popular because of their simplicity in application. These techniques are based on data, primarily acquired from the field and are well supported by the judgment and experience of an evaluator (Raghuvanshi et al. 2014a(Raghuvanshi et al. , 2014b. Some examples of expert evaluation techniques are; techniques proposed by Raghuvanshi et al. (2014a), Ayenew and Barbieri (2005), Guzzetti et al. (1999), Turrini and Visintainer (1998), Sarkar et al. (1995), Anbalagan (1992), Pachauri and Pant (1992) etc.

Multi-criteria decision analysis methods
These methods are semi-quantitative approaches which are mostly used for landslide susceptibility evaluation (Abija et al. 2020;Erener et al. 2016;Ahmed 2015;Feizizadeh et al. 2014;Kavzoglu et al. 2013;Feizizadeh and Blaschke 2012;Gorsevski and Jankowski 2010) and hazard zonation studies (Bera et al. 2019). The methods which are categorized under multi-criteria decision analysis are analytical hierarchy process (AHP), fuzzy set based analysis, weighted linear combination and ordered weighted average (Bera et al. 2019;Ahmed 2015;Feizizadeh and Blaschke, 2012). Recently to minimize the subjectivity during factors weighting fuzzy linguistic based approaches are also being used. Some aforementioned subgroups of multi-criteria decision methods are discussed in the following paragraphs. These are AHP, the concept of fuzzy logic and weighted overlay method.
Analytic hierarchy process (AHP) approach The AHP is a multi-criteria decision making approach which is utilized to evaluate the landslide hazard susceptibility and hazard zonation mapping. It is a systematic approach that involves problem definition, goals and alternatives determination, formulation of pair wise comparison matrix, weight determination and finding an overall priority (Pardeshi et al. 2013;Saaty 2008). Landslide is a complex process that is resulted by combination of various causative factors (Raghuvanshi, 2019;Chimidi et al. 2017). The AHP approach can be applied in absolute or a relative measurement of relationship between causative factors and landslides. In absolute measurement each alternative is compared with one ideal alternative whereas, in relative measurement each alternative is compared with many other alternatives. The absolute measurement approach is normative that is conditioned by what is known to be the best. In contrast, the relative measurement approach is conditioned by the experience and ability to judge observations by an evaluator thus, it is descriptive in nature (Pardeshi et al. 2013;Saaty and Vargas 2006;Saaty 2008). Each of the landslide causative factors can be taken as alternatives. Further, these causative factors are assigned with absolute numbers (1-9) based on their relative significance in inducing instability to the slope (Yagi 2003). Thus, comparison matrices are developed to determine Consistency Ratio (CR) and Consistency Index (CI) (Pardeshi et al. 2013).
Weighted overlay method Weighted overlay method is a simple and direct method that can be used for the evaluation of landslide hazard in the given area (Ayele et al. 2014;Erener and Duzgun 2012;Bachri and Shresta 2010;Intarawichian and Dasananda 2010). This method is based on the assumption that the conditions that has resulted into the past landslides if reoccur in future in other areas, again landslides can occur. In weighted overlay method, for each considered causative factor layers, every cell is reclassified based on the preferred scale as per the significance with respect to the contribution that particular factor class can have on landslide occurrence. Such preferred scale is generally taken from 1 to 9 with 9 being the most significant with respect to the landslide occurrence. Thus, all causative factor layers are combined in GIS environment by using Weighted Overlay Model equation (Eq. 1). Besides, respective weights are also assigned to the causative factors on 0-9 ordinal scale. The numeric numbers assigned to factor class are termed as ratings and the numbers assigned to respective causative factors are known as weight (Kanungo et al. 2006;Lee et al. 2004). The higher weight or rating to a causative factor or to its respective class represents its greater significance for landslide occurrence (Ayele et al. 2014;Kanungo et al. 2006).
Where; 'W i ' is the weight of ith factor map, 'S ij ' is the ith spatial class weight of 'jth' factor map, 'S' is the spatial unit value in output map.
Fuzzy logic method The processes, events or function that change continuously may possibly be un-defined in terms of true or false events, thus these events or activities are expressed in fuzzy manner. The things which are vague or not clear are known as fuzzy. The Fuzzy logic is considered to be a superset of Boolean logic which is a extension to address partial truth values that are completely false or completely true (Lee 2007;Chacón et al. 2006;Gorsevski et al. 2003). Fuzzy logic method is applied in landslide studies to avoid subjectivity owing to selection of various multiple causative factors (Meten et al. 2015;Chacón et al. 2006;Zadeh 1978). The Fuzzy logic method is based on the fuzzy set theory in which membership degrees of elements can have varying degrees of confidence in an interval of [0, 1] (Ross 2010;Lee 2007;Kanungo et al. 2006). There are various fuzzy operators that can be used, these are; 'fuzzy or', 'fuzzy and', 'fuzzy algebraic product', 'fuzzy algebraic sum' and 'gamma operator'.
Where 'μ Combination ' is the fuzzy membership function, 'μ A ' represents the membership value of map 'A' at specified location, and 'μ B ' represents the value of map 'B' at specified location, likewise 'μ C ' and so on.
Similarly, the 'fuzzy or' is equivalent to 'Boolean OR' (logical union) that represents the values of output membership, controlled by the maximum values for any of the input maps. Thus, the 'fuzzy or' can be represented as (Lee 2007;Zimmermann 1996) Further, the 'fuzzy algebraic product' operator can be expressed as (Zimmermann 1996;Bonham-Carter 1994) Eq.4; Where; 'μ i ' represents the fuzzy membership function for the 'i th ' map and 'i' represents 1, 2... n maps that are to be combined.
Similarly, the 'fuzzy algebraic sum' compliments to 'fuzzy algebraic product' and is expressed as (Lee 2007;Zimmermann 1996; Bonham-Carter 1994) (Eq. 5); Where; 'μ i ' represents the fuzzy membership function for the 'i th ' map and 'i' represents 1, 2... n maps that are to be combined.
Finally, the 'gamma operation' can be expressed in terms of the 'fuzzy algebraic product' and the 'fuzzy algebraic sum' as (Anbalagan et al. 2015;Zimmermann 1996;Bonham-Carter 1994) Where; 'λ'express a parameter that is chosen in the range of (0,1), the fuzzy algebraic sum is calculated through Eq. 5 and fuzzy algebraic product is calculated by using Eq. 4.
When λ is 1 in 'fuzzy gamma operation' the combination is identical to 'fuzzy algebraic sum' and when the value of λ is 0 the combination is identical to 'fuzzy algebraic product'. Thus, sensible selection of 'λ' may provide values as output that will make sure flexible negotiations among decreasing trends in fuzzy algebraic product and increasing tendencies in 'fuzzy algebraic sum' (Lee 2007). For landslide susceptibility study Lee (2007) and Ahmed et al. (2014) have assigned 'λ' values randomly in between 0 to 1; as 0.025, 0.05, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 0.95, and 0.975.

Statistical approaches
Statistical approaches are the most commonly used methods for landslide susceptibility and hazard zonation (Mengistu et al. 2019;Hamza and Raghuvanshi 2017;Girma et al. 2015;Kanungo et al. 2006). Reviews on various methods developed in past and subsequent evaluation and development of methods has been presented by Reichenbach et al. (2018), Corominas et al. (2014) and Fell et al. (2008). Some of the techniques, which were developed earlier, have been discarded, some methods were modified and new techniques are always being developed (Ayalew and Yamagishi 2005). A thorough literature review made during the present work revealed that the statistical techniques that were used for landslide susceptibility evaluation and zonation by previous researchers can be broadly classified into; Bivariate and Multivariate statistical approaches (Kanungo et al. 2006;Gorsevski et al. 2006aGorsevski et al. , 2006bGorsevski et al. 2000).
Each of these statistical methods showed significant heterogeneity in adopted thematic factors, selected scale, statistical modeling tools used and the methods used for the model performance evaluation and validation. In terms of performance and validation some of these statistical methods may produce better results; however no individual method has proved to be perfect in its performance. In fact integration of these methods into optimal model by utilizing the same thematic factors data and landslide inventory may ideally produce better results than a single model. It needs to be realized that the capability, skill set and experience of an evaluator in applying a particular statistical model is more important than the method by itself (Reichenbach et al. 2018;Corominas et al. 2014). In general, for landslide susceptibility evaluation and zonation approach through statistical modeling as such no established standards and code of practices are available. Thus, it is always a challenge in adopting an appropriate method for landslide susceptibility evaluation. Further, uncertainty always exists on account of its credibility and the results thus produced by the adopted method (Reichenbach et al. 2018;Guzzetti et al. 2006). General limitation of these statistical methods is that they do not propose mechanisms that control slope failure, but rather it assumes that the prediction of future landslide areas can be assessed by measuring the combinations of variables that have led to landslide occurrence in the past (Safaei et al. 2011;Lee et al. 2004;Guzzetti et al. 1999).
Bivariate statistical approach The bivariate statistical approach is based on an inductive logic, which suggests that if a situation holds in all observed cases than the situation holds in all cases. Thus, these techniques are based on a general assumption that "past and the present are the key for future" . For landslide studies it is generally assumed that the combination of conditions pertaining to various causative factors may possibly lead to landslide in a given area. Therefore, evaluation of these factors and their relation with the past landslides in the area may form the basis for the prediction of potential areas where landslides may occur in future (Chimidi et al. 2017;Negassa and Kala 2015;Girma et al. 2015;Lan et al. 2004;. In bivariate statistical techniques contribution of individual causative factors in inducing landslide in an area is evaluated separately. Such contribution of various causative factors in inducing landslides is assessed quantitatively through density ratio. The density ratio is derived through overlay analysis where each of the causative factor layers is overlaid over landslide inventory distribution layer. The overlay analysis will thus provide a ratio between landslide that 'did' occur to the landslide that 'did not' occur within each causative factor class (Hamza and Raghuvanshi 2017;Lee and Min 2001). Further, combining these density ratios for various factor classes will form the basis for landslide hazard evaluation and zonation in the given area. The commonly used methods for landslide susceptibility and hazard zonation that fall under bivariate statistical approach are; Frequency Analysis (likelihood ratio) method (Lee and Min 2001) (Lee 2007;Chacón et al. 2006;Kanungo et al. 2006;Gorsevski et al. 2003) etc.

(i) Frequency Analysis (likelihood ratio) method
Among the bivariate statistical approaches Frequency analysis is the most commonly used method (Chimidi et al. 2017;Hamza and Raghuvanshi 2017;Girma et al. 2015;Lee and Min 2001). This method utilizes the correlation between each of the responsible causative factor classes and the spatial distribution of past landslides in the area (Chimidi et al. 2017;Girma et al. 2015;Moung-Jin et al. 2014;Akgun et al. 2012;Pradhan and Lee 2009;Lee 2005). The frequency analysis method is based on the ratio of landslides in a factor class as a percent of landslides to the area of the factor class as a percent to the entire area. This frequency ratio is normalized to '1'. Thus, a frequency ratio for a parameter class greater than one shows strong correlation to the occurrence of landslide whereas a frequency ratio less than one indicates lower correlation of the factor class with the landslides (Chimidi et al. 2017;Girma et al. 2015;Lee and Min 2001). The frequency ratio (FR d ) can be expressed as (Eq. 7); Where; 'FR d ' is the Frequency ratio for the causative factor class, '%Ls' is the percentage of landslides in a causative factor class and '%Am' is the area of the causative factor class as a percentage of the entire map.
Further, the landslide susceptibility index (LSI) for each pixel is the summation of total overlapped pixels and is given by Eq. 8; If the LSI value is high the landslide hazard is considered to be high. The LSI can also be represented as landslide hazard index (LHI) (Pradhan and Lee 2009).

(ii) Weights of Evidence (WOE) method
Weight of evidence (WOE) method is widely used for the prediction of potential susceptible areas for landslides (Süzen and Doyuran 2004;Van Westen et al. 2003;Van Westen 1993). This method is used to combine the dataset of landslides through a quantitative data-driven approach. With the WOE method prior probability (Prior P ), Conditional probabilities (Cond P ) and negative and positive weights of landslide susceptibility can be calculated. For 'Prior P ' the landslide density has to be computed as the total number of pixels within the area covered by the landslides to the total number of pixels in the area under study. 'Prior P ' can be expressed as Eq. 9; Where; 'Prior P ' is the prior probability, P{S} is the conditional probability of having a landslide 'S'; N Pix (slide) is the total number of pixels within the landslides in the area of study and NP ix (total) is the total number of pixels in the area of study.
The Conditional probabilities (Cond P ) is the ratio between the total number of pixels of a corresponding factor class within the landslides to the total number of pixels of that factor class within the area of study (Van Western 2002). This can be expressed by Eq. 10; Where; 'Cond p ' is the conditional probability of having landslide 'S' in a given parameter class 'B', 'N Pix {SnB}' is the total number of pixels of a parameter class 'B' within the landslides and 'N Pix {B}' is the total number of pixels of that factor class 'B' within the total area of study.
Weight of evidence (WOE) method determines the weight for all considered predictive causative factors (B) with respect to absence or presence of the landslides in the area (Mohammady et al. 2012;Bonham-Carter 1994).
Where, 'In' is the natural log, 'P' is the probability, 'Bi' is the presence of potential landslide predictive causative factor, ' B i ' is the absence of the potential landslide predictive causative factor, 'S' is the presence of landslide and ' S ' is the absence of landslide, ' W þ i ' represents a positive weight which implies that predictable factor is present in landslide, ' W − i ' represents a negative weight which suggests absence of predictable factor in the landslide.

(iii)Information Value (IV) method
In information value (IV) method the weighted class value is computed through density of landslides with respect to each causative factor. For spatial landslide prediction and zonation, information value method has been utilized by many researchers such as; Mengistu et al. (2019), Balasubramani and Kumaraswamy (2013), Yalcin (2008), Lin and Tung (2003), Zêzere et al. (2004), Jade and Sarkar (1993), Yin and Yan (1988) etc.
The information values of various causative factors can be used to determine the possible areas of landslide occurrence which may facilitate in landslide hazard zonation. Based on the presence or absence of the causative factor classes within the past landslides the Information values can be determined. In order to calculate the weights for respective classes the landslide map can be combined with the causative factor maps. Thus, landslide density for causative factor sub-classes can be determined by overlaying the causative factors map on the inventoried landslide map (Mengistu et al. 2019). According to Yin and Yan (1988), if the information value is positive the causative factor class represents strong relationship with the landslides in the area. According to Van Westen et al. (1997) the weighted value of a causative factor class can be represented as the natural logarithm of density of landslide in a factor class, divided by landslide density in the total map area. The Information value method is also named as landslide index (W i ) method. The information value can be computed by utilizing eqs. 13 to 16 (Yin and Yan 1988); Weight of factor class WBi ð Þ¼ CP PP ð15Þ Where; 'N Pix {SBi}' is the number of landslide pixels within the factor class, N Pix {Bi}is the number of pixels of a factor class, 'N Pix {TS}' is the total sum of pixels of landslide of the whole study area and N Pix {A}is the total pixels of the whole study area.
Once the Information values are obtained, the respective IVs can be assigned to each factor class and the weighted causative factor maps can be prepared. These causative factor maps can further be processed in GIS environment and landslide susceptibility index (LSI) can be computed (Eq. 17) for each pixel with the help of raster calculated. Thus, this LSI will form the basis to produce the landslide hazard zonation map (Mengistu et al. 2019).
Where; LSI is the landslide susceptibility index and ' ∑IVB′ is the sum total of information values for all causative factors.
Multivariate statistical approach The landslide susceptibility and hazard zonation through multivariate statistical approach is based on the relationship of relative contribution of each of the causative factor to the total landslide susceptibility in the area (Nandi and Shakoor 2009;Kanungo et al. 2009;Ayalew and Yamagishi 2005;Süzen and Doyuran 2004;. In multivariate methods for landslide susceptibility analysis, percentage of landslides for each pixel is determined and data layer on landslide presence or absence is developed through statistical analysis. The commonly used multivariate statistical methods for landslide susceptibility and hazard zonation are; Logistic regression model, multiple regression models and Discriminant analysis (Kanungo et al. 2006;Guzzetti et al. 1999;Chung and Fabbri 1995;Van Westen 1993;van Westen 1994;Yin and Yan 1988).

(i) Logistic regression method
Among various multivariate statistical methods Logistic regression is the most popular method used for the spatial prediction of landslide susceptibility and hazard zonation (Wu et al. 2015;Schicker and Moon 2012;Gorsevski et al. 2000). The logistic regression method is capable of predicting a binary response variable such as presence or absence of landslides through categorical and continuously scaled predictors (Steger et al. 2016;Goetz et al. 2015;Felicísimo et al. 2012;Pradhan and Lee 2009;Van Den Eeckhaut et al. 2006;Atkinson and Massari 1998). This method is applied to predict the probability of occurrence of landslides following the logistic regression statistical analysis. The relationship between the occurrence of landslides in the area and its dependency on several variables (causative factors) can be defined through Eq. 18 (Ayalew and Yamagishi 2005;Menard 1995).
Where; 'P' is the estimated probability of a landslide occurrence which varies from 0 to 1 on a'S shaped' curve. The term 'z' represents a linear combination. The logistic regression utilizes fitting of an equation (19) to the data set; Where; 'b 0 ' represents the intercept of the model, the b i (i = 0, 1, 2, …, n) are the slope coefficient of the logistic regression model and the x i (i = 0, 1, 2, 3………..,n) represents the independent variable.
Logistic regression is applicable for the situations where either the data is continuous or discrete or combination of both. This is the main advantage of logistic regression method. However, the results of logistic regression cannot identify the influence of different classes on landslide occurrence. To apply logistic regression the dependent variable should be binary, for example yes or no, zero or one, absence or presence etc. (Lee 2015;Chen and Wang 2007).

(ii) Discriminant method
Discriminant method is a commonly used method for landslide susceptibility and hazard zonation. It is a multivariate statistical model in which the dependent variable is considered to be categorical rather than continuous (Gorsevski et al. 2000). The discriminant analysis facilitates to know the maximum difference for each landslide causative factors (independent variable) among non-landslide group and the landslide group. Thus, it forms the basis to calculate the weights for these causative factors (Pardeshi et al. 2013;Lee et al. 2008). The slope units can be classified to account for landslide free and the landslide affected classes. Further, by determining the Standardized Discriminant Function Coefficient (SDFC) relative significance of each variable can be expressed in terms of discriminant function as a predictor of potential slope instability. Variable that shows high coefficient are strongly correlated to absence or presence of the landslide (Pardeshi et al. 2013;Guzzetti et al. 2005b).

Artificial intelligence (AI) methods
Artificial intelligence (AI) method uses some of the statistical concepts. However, these methods are based on assumptions, predetermined algorithms and output. AI methods are suitable when a direct mathematical relationship cannot be established between cause and effect (Chowdhury and Sadek 2012).There are number of AI methods or machine learning methods (Kavzoglu et al. 2018) that can be used for the landslide studies. These can be categorized as; artificial neural network (ANN), fuzzy based, hybrid, kernel based and tree based (decision tree) methods. Further, there are number of sub methods under each category, for instant, ANN-adaptive neuro-fuzzy inference system, back-propagation neural network etc.; fuzzy based-fuzzy clustering; hybrid-SVMparticle swarm optimization; kernel based-kernel logistic regression, SVM (support machine vector); tree basedbagging, decision trees, random forest etc. (Đurić et al. 2019;Kavzoglu et al. 2018;Goetz et al. 2015) etc. AI methods are effective regardless of data dimension (i.e., the number of conditioning factors) and data type (applicable to deal with discrete and continuous data) (Đurić et al. 2019;Kavzoglu et al. 2018). Furthermore, they can reveal good generalization performance on many real life issues, they have few parameters to adjust and provide the architecture of learning machines without experimentation (Pawley et al. 2017). Hence, they are more suitable for the analysis of high dimensional data and complex systems.
Artificial neural network (ANN) method ANN method facilitates to obtain, represent and perform mapping of landside susceptibility and hazard from one multivariate space of information into another by providing a set of data or information relating to representative mapping (Pradhan and Lee 2010;Nefeslioglu et al., 2008;Garrett 1994). The landslide is a complex process which is resulted due to combination of various causative and triggering factors. Also, the relationships between the landslides and the causative and triggering factors are considered to be non-linear in nature (Ercanoglu 2005). Therefore, in order to address such non-linear complex relationships between the factors and the landslides ANN method is used which includes back propagation learning algorithm and is capable of developing rules for weight assignment for the respective factors (Kanungo et al. 2006;Yesilnacar and Topal 2005;Gómez and Kavzoglu 2005;Lee et al. 2004;Zhou et al. 2002;Haykin 1999;Ripley 1996).
In ANN, back propagation is used to determine a gradient (Kanungo et al. 2006) which is required to calculate the weights that are used in the network. For landslide susceptibility and hazard zonation back propagation is used as shorthand to the backward propagation of the errors. For this, at the output an error is determined which is distributed backward all along the layers of the network (Pradhan and Lee 2010;Paola and Schowengerdt 1995).The ANN facilitates the database management with qualitative or quantitative information having rules (if/ then) that can infer results. The entry nodes can have quantitative, qualitative, reclassified or direct data. For landslide susceptibility study, landslide densities of different types or determinant causative factor data can be assigned to the entry nodes. The data is processed in the hidden layer nodes by following rules for weighting that depends on the considered importance (Chacón et al. 2006;Zell et al. 1993;Lipmann 1987). Further, the ANN method can manage imperfect or incomplete data and it also allows analyzing nonlinear complex variable relationships between the causative factors and the landslides (Chacón et al. 2006;Lee et al. 2003).The major limitations with ANN method are in its computing time, required for data format conversion to GIS environment and the application complexity owing to the process in the hidden layers (Chacón et al. 2006). The ANN method is only applicable for landslides which are shallow and it has limitation for deep seated landslides (Basma and Kallas 2004;Neaupane and Achet 2004).

Support vector machines (SVM)
This method is used to solve complex classification and regression problems (Tien Bui et al. 2015;Kavzoglu et al. 2013;Ballabio and Sterlacchini 2012;Kavzoglu and Colkesen 2009;Brenning 2005). SVMs are originally developed as a binary classifier aiming to find a linear hyperplane that separates two classes optimally (Vapnik 1999). For binary classification problems, SVMs attempt to find a separating hyperplane in the feature space such that the distance between the positive and negative samples is maximized for the linearly separable case (Kavzoglu et al. 2013). The hyperplane providing maximum margin between two classes is called optimum hyperplane and the points that constrain the width of the margin are called support vectors. In many classification and regression problems, it is difficult to separate data linearly. In such cases, the technique can be extended to allow for nonlinear decision surfaces (Cortes and Vapnik 1995).The main idea behind SVM is to find the boundary line that separates the two classes, but in such a way that the boundary line creates a maximum separation between the classes (Han et al. 2019;Kavzoglu et al. 2018;He et al., 2012).

Probabilistic approach
In order to evaluate landslide susceptibility and hazard zonation the degree of relationship between the past landslide distribution and the causative factors is converted to a value which is based on a probability distribution function. Probabilistic approach helps in prediction of spatial and temporal landslide distribution probability in the given area (Lari et al. 2014;Guzzetti et al. 2005b). The approach utilizes the comparison of spatial landslide distribution with the considered causative factors (explanatory variables) in a framework of probabilistic theory (Lari et al. 2014;Kanungo et al. 2009;Straub and Schubert 2008). The probabilistic approach though considered quantitative however it has certain degree of subjectivity in assignment of weights to various causative factors (Kanungo et al. 2006). Thus, probabilistic approach may be considered as semiquantitative.

Validation methods
All statistical methods, machine learning or artificial intelligence and some of the semi-quantitative approaches are tested by field observations and statistical tests. Sample field observation may be applicable for every type of landslide studies except for small area or single landslide studies. However, statistical validations are used for statistical, AI and semi-quantitative approaches such as AHP, fuzzy logic etc. Receiver operator characteristic (ROC) is used for statistical, AI and semiquantitative approaches. For multivariate statistical approaches there are different validation methods starting from data reliability to model accuracy in addition to ROC. In order to check the performance of the prediction ability of the above mentioned methods the most popular technique is receiver operator characteristic (ROC) curve. The ROC curve represents a plot of the probability with true positive identified landslides against the probability of false identified landslides, as the cutoff probability varies (Gorsevski et al. 2006a). With the selection of decision parameter that adjust the tradeoff between the proportion of correctly identified pixels and the incorrectly identified pixels, the ROC performance curve can be varied systematically (Kavzoglu et al. 2014;Swets 1988). Other validation methods that can be used are landslide density analysis, landslide percentage comparison, relative landslide density index; relative error etc.

Advantage and limitation of each method
The different methods used for landslide susceptibility analysis has their own advantage and disadvantage owing to its application, data procurement and scales at which these can be applied (Raghuvanshi et al. 2014a(Raghuvanshi et al. , 2014bFall et al. 2006;Kanungo et al. 2006;Casagli et al. 2004;Guzzetti et al. 1999;Leroi 1997).
Geomorphological methods can be applied in relatively short period of time at low cost. These methods can be applied in relatively large areas (small scale map) and can cover inaccessible area through various multitemporal satellite image interpretations (Chelli et al. 2015;Delmonaco et al. 2013;Guzzetti et al. 1999). Since these methods are applied at small scale therefore geomorphological details considered in these methods are limited (Guzzetti et al. 2005a(Guzzetti et al. , 2005bIbsen and Brunsden 1996). Thus, the reliability and application of results is limited for specific purpose. Similarly, the limitations of landslide inventories refer to their subjectivity and to the difficulty of assessing their reliability. The reliability of archive inventories depends largely on the quality and abundance of information sources (Raghuvanshi et al. 2014a;Guzzetti et al. 2005aGuzzetti et al. , 2005bCruden and Varnes 1996;Ibsen and Brunsden 1996;).
Further, the advantage in heuristic evaluation technique is in its simplicity. The weight assigned by the expert to various causative factors do not need historical data base and are simply based on the judgment of relative contribution of each parameter class. These techniques are popular because of their simplicity in application. These techniques are based on data, primarily acquired from the field and are well supported by the judgment and experience of an evaluator (Raghuvanshi et al. 2014a(Raghuvanshi et al. , 2014b. However, the major disadvantage of these techniques is in its subjectivity in assigning weight and ratings to the parameter classes (Raghuvanshi et al. 2014a;Delmonaco et al. 2013;Guzzetti et al. 1999). Thus, the landslide susceptibility evaluated by heuristic technique may vary considerably if evaluated by different experts.
The statistical approaches are primarily based on statistical relationship between past landslides and various causative factors. In these techniques the weights for various causative factors are determined statistically . These techniques are based on the analysis of functional relations between instability factors (thematic variables) and the past and present distribution of slope failures (landslide inventory) (Negassa and Kala 2015; Guzzetti et al. 2000). The advantage of statistical techniques is that these techniques can be applied over a large area and the ratings for various factor classes can be deduced statistically from the past landslide data. However, the limitation of these methods is in collection of past landslide inventory data over large areas. Since the general rules for landslide susceptibility are formulated from the past landslide data therefore well distributed past landslides in the study area will ensure good quality results. Thus, the type, extent and distribution of past landslide data in the study area is a necessary input for all statistical techniques (Negassa and Kala 2015).
Additionally, the statistical models require large efforts to collect and validate the necessary input data, which are often not readily available. Besides, for better results, interaction is required between geomorphologists and the statisticians to process the landslide and geoenvironmental data to avoid unrealistic results. Further, statistical models are influenced (negatively) by the extent of the study area, making it difficult to compare susceptibility classes from different locations. Also, often susceptibility maps prepared by statistical methods are difficult to understand by a non-specialists, including planners and policymakers (Reichenbach et al. 2018;Raghuvanshi et al. 2014aRaghuvanshi et al. , 2014bFall et al. 2006;Van Westen et al. 2003;Van Westen et al. 1997).
Artificial intelligence (AI) applications are utilized to simulate human intelligence for either solving a problem or making a decision. AI provides the advantages of permanency, reliability, and cost-effectiveness while also addressing uncertainty and speed in either solving a problem or reaching a decision (Chowdhury and Sadek 2012).
The main advantage of using ANN methods is that there is no need for statistical variables in the process. The ANN methods allow the target classes to be defined in relation to their distribution in each source data set, facilitating the integration of data emerging from remote sensing, GIS, etc. Further, the time involved is shorter than required for most statistical techniques and calculation pixel by pixel is possible. The method allows incomplete or imperfect data to be managed as well as the analysis of interactions between non-linear or complex variables (Chacón et al. 2006;Lee et al. 2003). According to Chacón et al. (2006) ANN has limitations that are complexity for application due to the internal process within the hidden layers and the amount of computing time needed for changes in data format to be used in GIS. Also it has limited application in deep landslides, rather it can mainly be used for shallow landslide (Basma and Kallas 2004;Neaupane and Achet 2004).
Geotechnical modeling or deterministic approach relies upon the understanding of the physical laws that control the slope instability (Guzzetti et al. 2000). Geotechnical modeling are useful for quantitative investigation for landslides, however these techniques are useful only in areas of limited extent due to difficulty of collecting geotechnical data with appropriate resolution over larger regions (Negassa and Kala 2015; Kanungo et al. 2009;Aleotti and Chowdhury 1999).

Selection criteria
As discussed, there are several techniques available for evaluation of landslide susceptibility and hazard zonation. The choice of selection of appropriate technique for landslide study will depend on the criteria such as; investigation purpose, the extent of the area to be covered, the types of mapping units, the scale of map to be produced, type of data to be used, type of landslides, availability of resources, capability and skill set of evaluator, accessibility of the study area etc. (Guzzetti et al. 2000(Guzzetti et al. , 2012Van Westen et al. 2005).
The landslide studies may be carried out for various purposes these may be primarily for construction sites, regional land-use planning, and evaluation of susceptibility of an area, for hazard zonation and mitigation measures. For each of these purposes specific type of technique has to be followed. Landslide studies for new construction sites requires techniques that may provide quantitative results of slope stability that may directly be utilized for design purpose. These studies are site specific and needs to be carried out at large scale (Raghuvanshi et al. 2014a(Raghuvanshi et al. , 2014bFall et al. 2006). The most appropriate techniques which may provide required results are deterministic techniques (Chimidi et al. 2017;Negassa and Kala 2015;Kanungo et al. 2006). Further, if the purpose is for regional planning, depending on the area coverage, heuristic or statistical techniques may be followed. The heuristic techniques generally require considerable field data (Ermias et al. 2017;Raghuvanshi et al. 2014aRaghuvanshi et al. , 2014b therefore these techniques may be conveniently applied at local scale. However, statistical techniques may be conveniently applied at medium scale. In some cases, depending on the purpose, integration of techniques may also be followed to get the desired results. The scale of mapping will depend on the extent of the area to be covered, the purpose for which landslide study has to be carried out and the availability of data and resources (Huabin et al. 2005;Aleotti and Chowdhury 1999). The different scales proposed for landslide studies has been classified as; detail scale (> 1:5000), large scale (< 1:5000-1:10,000), Medium scale (1:25,000-1:50,000) and Regional scale (< 1:250,000) (Huabin et al. 2005). Similarly, IAEG, Commission on Landslides (1990) and CEOS (2001) have classified landslide mapping scales as; Large scale (1:5000-1:15,000), Medium scale (1:25,000-1:100,000), Regional scale (1:100,000-1:500,000) and National scale (< 1:500,000). The detail scale or large scale are applicable for construction sites where deterministic techniques are applied for slope stability evaluation and where only single element at risk is considered (Corominas et al. 2014). The medium scale is applicable for land-use planning where evaluation of susceptibility of the area and hazard zonation is required. The landslide studies at medium scale are required at the preliminary stage of regional developmental plans or for the evaluation of the possible instability problems related to the development of large engineering projects such as; dam and reservoir projects.
The maps prepared at medium scale delineate the areas that are susceptible for landslides and may be a constraint for the development of engineering projects. However, they do not provide required data for design purpose (Soeters and Van Westen 1996). Further, landslide studies at national mapping scale are carried out covering very large areas. Studies at such scale are basically carried out to facilitate national policy makers. These maps will provide low level of details as the investigation is carried out by following general rules (Corominas et al. 2014;Soeters and Van Westen 1996). Table 1 shows the application of various landslide susceptibility and hazard zonation techniques at different scales, as sourced and modified from the works of Soeters and van Westen (1996) and Aleotti and Chowdhury 1999).
Other important criteria related to selection of appropriate technique for landslide studies is related to the type of data to be used. The method by which the required primary data on causative factors and landslide inventory will be collected and the secondary data that is available often forms the basis to adopt the appropriate technique for landslide studies. If systematic multi-stage historical landslide registered data in good quality and quantity is available statistical approaches may be applied effectively. However, such systematic data, in most of the cases, is not maintained by the concerned authorities (Girma et al. 2015;Negassa and Kala 2015). Thus, the landslide inventory data is produced through remote sensing and field based investigations. However, the manifestations for past landslides may not always be visible due to vegetation cover or due to displaced soils or due to erosional processes. Thus, the landslide inventory data generated through remote sensing or through field investigation may not be complete and the landslide susceptibility evaluated may not always produce actual results. Also, if systematic multi-temporal data on past landslides and other governing factors is amply available both statistical and probabilistic approaches can be applied effectively for landslide studies. The deterministic techniques are scale dependent and may require good quality of geological and geotechnical data (Raghuvanshi et al. 2014a). Therefore, deterministic techniques needs to be adopted only when the study area is relatively small and the possible resource to generate good quality geological and geotechnical data is available. Further, most of the heuristic approaches do not require detailed landslide inventory data however; they may require considerable field data (Girma et al. 2015;Raghuvanshi et al. 2014a) on various parameters, particularly on geological factors. Thus, the heuristic techniques may only be considered when the study area is relative small, area is well accessible for field investigation and sufficient time and resource is available.

Discussion
Landslide is a complex process that is resulted due to combination of several inherent causative and external triggering factors (Raghuvanshi et al. 2014a;Anbalagan 1992). For landslide susceptibility evaluation it is important to know these factors and to asses there possible influence in inducing instability to the slopes. A basic assumption which is followed in most of the techniques in landslide studies, particularly for quantitative techniques is that; the combination of causative factors that has lead to the past landslides in an area if reoccur in other areas, again landslides can be expected. Therefore, it is important to evaluate causative factors and assess their possible relationship with the past landslides in the area; this will form the basis for the prediction of the potential areas where landslides may occur in future (Chimidi et al. 2017;Girma et al. 2015;Negassa and Kala 2015).
The qualitative approaches can effectively be applied at small or large scale and are used to delineate the potential landslide areas. The landslide inventory approach, which is direct mapping method, are used to show the landslide distribution in the area and they also provide data for landslide density analysis (Reichenbach et al. 2018;Guzzetti et al. 1999;Soeters and Van Westen 1996). The landslide inventory mapping is the most straight-forward and are basic input for most of the quantitative approaches (Chimidi et al. 2017;Negassa and Kala 2015;Lan et al. 2004;. Landslide inventory approaches are focused to analyze distribution of landslides, to classify the landslides, to evaluate the temporal change in pattern of landslides and to determine density of landslides in various terrain units (Guzzetti et al. 2012). Further, the indirect method includes geomorphological mapping and techniques that are based on qualitative map combinations. These techniques include assessment of parameters that have significant influence on landslide process. These techniques involve combination of various causative factor maps. For this, based on the expert opinion weight and ratings are assigned to causative factors and the factor classes (Reichenbach et al. 2018;Raghuvanshi et al. 2014a;Ayalew and Yamagishi 2005;Anbalagan 1992). In general, the indirect approaches; heuristic approaches are subjective in nature and greatly depend on the expert opinion and evaluation (Guzzetti et al. 1999;Aleotti and Chowdhury 1999). However, they are simple in application and can utilize more actual field observed data well supported by the experience of an evaluator (Raghuvanshi et al. 2014a).
Most of the statistical techniques can effectively be applied at medium scale (Girma et al. 2015;Fell et al. 2008;Lee and Min 2001). Also, statistical techniques can be applied at small scale however the results thus produced are less reliable. Since, statistical techniques are based on historical landslide inventory data therefore collection of such data over very large areas is practically not feasible. Thus, because of this limitation at small scale the statistical techniques may not provide reliable results.
Artificial intelligence (AI) methods have an advantage of permanency, reliability, and cost-effectiveness and they also address uncertainty and rapid means of solving a problem for quick decision (Chowdhury and Sadek 2012).Further, the ANN methods allow the target classes to be defined in relation to their distribution in each source data set, facilitating the integration of data produced from various sources. The process requires less time as compared to other statistical techniques. The method allows analysis of interactions between nonlinear or complex variables (Chacón et al. 2006;Lee et al. 2003). ANN has limited application in deep landslides and can conveniently be used for shallow landslide evaluation (Basma and Kallas 2004;Neaupane and Achet 2004). The deterministic techniques are most suitable when applied at large scale, as these techniques require considerable field data on geological and geotechnical parameters. Ideally, they are most suitable for construction sites and these techniques may produce results that can directly be used for the design purpose (Chimidi et al. 2017;Negassa and Kala 2015;Kanungo et al. 2006). Moreover, the quantitative techniques are reliable as these are objective in nature and the predictions for landslide susceptibility or slope instability is based on the actual realistic data and interpretations. The quantitative techniques overcome the subjectivity owing to expert evaluations and are capable of producing repeatable results (Safaei et al. 2011;Carrara et al. 1991;Carrara et al. 1992;Carrara et al. 1995).
Often it is difficult to make a choice on the selection of an appropriate approach for landslide susceptibility evaluation and hazard zonation mapping. It is reasonable to consider factors that may guide to adopt an appropriate approach. These factors are related to investigation purpose, the extent of the area to be covered, the types of mapping units, the scale of map to be produced, type of data to be used, type of landslides, availability of resources, capability and skill set of an evaluator, accessibility of the study area etc. (Guzzetti et al. 2012;Guzzetti et al. 2000;Van Westen et al. 2005).
The purpose for which the landslide studies needs to be carried out is the most important factor in deciding about the technique. For regional planning purpose statistical or heuristic techniques are most appropriate approaches. The statistical techniques for this purpose may conveniently be applied at medium scale whereas the heuristic techniques may be applied in relatively small area at both medium and large scales. At regional scales heuristic techniques are not practically feasible as these techniques generally requires considerable field data on various causative factors (Ermias et al. 2017;Raghuvanshi et al. 2014a). Therefore, heuristic techniques may be conveniently applied at a local scale. The most appropriate techniques that may be applied at new construction site are the deterministic techniques. These techniques can be applied at large scales, most effectively applied for individual slopes and are capable of providing quantitative results that can directly be applied in engineering design (Chimidi et al. 2017;Raghuvanshi et al. 2014a;Fall et al. 2006;Kanungo et al. 2006). Landslide and related slope stability studies may be required for new project sites that may include road construction sites, canal alignments, building sites on slopes or adjoining to slopes and valley slopes along the dam reservoirs. For all such construction sites the landslide studies needs to be carried out by deterministic or probabilistic approaches, even if the study area is large. For such construction sites the required results must be of quantitative nature that can directly be utilized for the design purpose. The data generated through heuristic or statistical techniques may not be useful for the construction sites as these techniques are only capable of delineating potential hazard zone of landslide susceptibility and the results thus produced cannot be utilized for design purpose. Moreover, for the preliminary study for the road design or for the town planning, inventory or statistical approaches may be adopted to identify the potential landslide and related unstable areas. However, at detailed investigation stage deterministic or probabilistic approaches can be followed.
For landslide studies selection of appropriate mapping units is also important. Different type of mapping units such as; grid cell (pixels), slope unit, terrain units, unique condition units, geo-hydrological units, topographical units and political or administrative units may be adopted (Lee 2015;Guzzetti et al. 1999). For statistical and physical based approaches, grid cell or pixels mapping units in the raster data sets are most appropriate (Reichenbach et al. 2018;Lee 2015;Guzzetti et al. 2006). Similarly, unique condition units are commonly used in probabilistic approaches in both raster and vector data sets (Guzzetti et al. 1999;Carrara et al. 1995;Chung and Fabbri 1995). The primary advantage of adopting unique condition units is in their simplicity by which they can be obtained in GIS environment and in fact relatively they reduces some of the conceptual and operational problems of the grid cells (Reichenbach et al. 2018;Lee 2015) .
Further, the basis to adopt appropriate technique for landslide studies is also guided by the method by which the required data on causative factors will be collected. The collection of required data on causative factors and landslide inventory for various techniques may differ considerably. For statistical approaches systematic multistage historical landslide registered data in good quality and quantity is mandatory, as this data will form the basis for the prediction of potential landslides in the given area. Further, landslide inventory data is often not well maintained by the concerned authorities. Therefore, as a general practice the landslide inventory data is produced through remote sensing and field based investigations. This consumes considerable time and resources. Therefore, availability of data on various causative factors and systematic landslide inventory data will not only save the time but it will also ensure good quality results. The other aspect related to availability of data is related to the area coverage. If secondary data on causative factors and the landslide inventory is available in good quality and quantity, relatively large area can be covered for the landslide study. In contrast non availability of required secondary data for statistical approach will force to restrict the landslide study to limited coverage area. Similarly, the deterministic techniques are scale dependent and may require good quality of geological and geotechnical data (Raghuvanshi et al. 2014a). These techniques can only be adopted when the study area is relatively small and the required resource and time to generate good quality of geological and geotechnical data is available.
For relatively small areas, qualitative approaches particularly heuristic approaches are effective for the delineation of the potential landslide susceptible areas. Most of the data collections in the direct method or geomorphological analysis are through remote sensing sources, supplemented and verified from the field observations. In general, medium scale mapping is used to investigate landslide susceptibility and hazard zonation by using statistical, heuristic and probabilistic approaches. Large scale mapping is used to investigate specific area particularly for construction purpose or for landslide hazard mitigation.

Conclusion
For proper landslide mitigation and management, susceptibility evaluation and hazard zonation is very important. Landslides are influenced by several inherent causative and external triggering factors that vary significantly from place to place. Proper evaluation of these factors is essential for landslide studies. Each of these factors may influence landslide process and in combination they result into landslide activity. Further, various approaches that are used for landslide evaluation studies can broadly be divided into qualitative and quantitative approaches. Qualitative approaches include geomorphological analysis (inventory approach) and heuristic techniques (knowledge based approach) whereas quantitative approach includes statistical and deterministic techniques. The qualitative approaches can effectively be applied at small or large scale and are used to delineate the potential landslide areas. Statistical approaches are the most commonly used methods for landslide susceptibility and hazard zonation. The statistical techniques can be classified into; Bivariate and Multivariate statistical approaches. In bivariate statistical techniques contribution of individual causative factors in inducing landslide in an area is evaluated separately. Such contribution of various causative factors in inducing landslides is assessed quantitatively through density ratio. Combining these density ratios for various factor classes will form the basis for landslide hazard evaluation and zonation in the given area. The commonly used methods under bivariate statistical approach are; Frequency Analysis (likelihood ratio) method, Weighted overlay model, Weights of Evidence Model, Information Value Model (IV) and Fuzzy logic methods. The multivariate statistical approach is based on the relationship of relative contribution of each of the causative factor to the total landslide susceptibility in the area. In multivariate methods, percentage of landslides for each pixel is determined and data layer on landslide presence or absence is developed through statistical analysis. The commonly used multivariate statistical methods are; Logistic regression model, multiple regression models, Discriminant analysis, conditional analysis and Artificial Neural Networks (ANN).
The quantitative techniques are reliable as these are objective in nature and the predictions for landslide susceptibility or slope instability is based on the actual realistic data and interpretations. The quantitative techniques overcome the subjectivity, owing to the expert evaluations and are capable of producing repeatable results. Further, for the selection of appropriate technique for landslide hazard evaluation and zonation the factors that may be considered to adopt an appropriate approach are; investigation purpose, the extent of the area to be covered, the types of mapping units, the scale of map to be produced, type of data to be used, type of landslides, availability of resources, capability and skill set of an evaluator and accessibility of the study area. Finally, landslide is a complex natural geo-hazard and it requires good understanding on the governing factors. It is necessary to adopt appropriate technique for the landslide evaluation and zonation, for which above mentioned selection criteria needs to be followed.