Evaluation of different machine learning models and novel deep learning-based algorithm for landslide susceptibility mapping

The losses and damage caused by landslide are countless in the world every year. However, the existing approaches of landslide susceptibility mapping cannot fully meet the requirement of landslide prevention, and further excavation and innovation are also needed. Therefore, the main aim of this study is to develop a novel deep learning model namely landslide net (LSNet) to assess the landslide susceptibility in Hanyin County, China, meanwhile, support vector machine model (SVM) and kernel logistic regression model (KLR) were employed as reference model. The inventory map was generated based on 259 landslides, the training dataset and validation dataset were, respectively, prepared using 70% landslides and the remaining 30% landslides. The variance inflation factor (VIF) was applied to optimize each landslide predisposing factor. Three benchmark indices were used to evaluate the result of susceptibility mapping and area under receiver operating characteristics curve (AUROC) was used to compare the models. Result demonstrated that although the processing speed of LSNet model is the slowest, it still significantly outperformed its corresponding benchmark models with validation dataset, and has the highest accuracy (0.950), precision (0.951), F1 (0.951) and AUROC (0.941), which reflected excellent predictive ability in some degree. The achievements obtained in this study can improve the rapid response capability of landslide prevention for Hanyin County.


Instruction
Landslide is defined as the special geological phenomenon that is threatening to mankind triggered by human activities or natural factors.Under the dual background of human activities and natural transmutations, the occurrence rate of landslides in the world increased rapidly (Sun et al. 2020).Depending on the latest statistical report of Ministry of Natural Resources of China (http:// www.mnr.gov.cn/), in total 6,181 geological hazards occurred in 2019, during that year, 211 people were killed and direct economic loss valued at 397 million dollars.Due to complex terrain, tectonic development and human activities, there are more than 230,000 potential geological hazards in China and landslides account for 53.50% of it.Therefore, the prevention of landslide development in China is crucial.In the face of increasingly serious landslide threats, the development of disaster prevention and mitigation work can effectively reduce the threat posed by landslides.In order to plan and construct the city safely and effectively, and to carry out the work of disaster prevention and mitigation successfully, it is necessary to quantitatively assess the landslide susceptibility on the regional scale.
Over the past years, various techniques and methods for forecasting the landslides have been applied to landslide susceptibility assessment (LSA).At first, with the help of historical landslide data and the geological Open Access *Correspondence: whysxdj2021@163.comenvironment background, the scope of landslide prediction was directly delineated, but this method relies heavily on experience, resulting in low reliability of the results (Guzzetti et al. 2012;Youssef and Pourghasemi 2021).With the development of geographic information system (GIS) and satellite remote sensing technology within each subject area, the more commonly used solutions can be summarized as a few steps.The first step of LSA is to collect the development characteristics and spatial distribution features of historical and hidden danger landslides (Pradhan and Lee 2010).Then the predisposing factors of landslide occurrence are selected from the geological and environment background.Subsequently, the linear or non-linear mapping relationship between predisposing factors and the degree of landslide susceptibility is analyzed by using evaluation model (qualitative or quantitative), and the contribution rate of each landslide predisposing factor is determined.In the end, some techniques of analysis and comparison are used to choose the suitable model for the study area (Carrara et al. 1995).
In principle, the evaluation models used in LSA could be crude categorized as two classes: statistical model and machine learning model.Usually, the statistical models were conjugated to GIS for spatial variety prediction of landslide disaster.The frequently used models are index of entropy (IOE) (Constantin et al. 2011;Youssef et al. 2015), analytical hierarchy process (AHP) (Kayastha et al. 2013), frequency ratio (FR) (Umar et al. 2014;Razavizadeh et al. 2017), certainty factor (CF) (Fan et al. 2017;Li and Zhang 2017), logistic regression (LR) (Pourghasemi et al. 2013;Aditian et al. 2018) and evidential belief function (EBF) (Carranza 2015;Li and Wang 2019).As the predictive ability of statistical model is still deficient and with the development of information technology, machine learning models could be a better alternative to solve the susceptibility assessment problem, such as artificial neural networks (ANN) (Aditian et al. 2018;Polykretis and Chalkias 2018;SOMA et al. 2019), support vector machines (SVM) (Bui et al. 2016;Pandey and Pourghasemi 2020), adaptive neuro-fuzzy inference systems (ANFIS) (Aghdam et al. 2017), decision trees (DT) (Pham et al. 2016), fuzzy logic (FL) (Saadoud et al. 2018) and multivariate adaptive regression spines (MARS) (Conoscenti et al. 2014).Due to data quality, factor selection, model parameter adjustment and other factors, some low accuracy, over fitting, and under fitting problems often appear (Bui et al. 2018).In order to solve these problems, hybrid model was developed in recent years, such as reduced error pruning trees (REPT) (Pham et al. 2019), kernel logistic regression model integrated with fractal dimension (KLR box-counting ) (Zhang et al. 2019), support vector regression model integrated with gray wolf optimization algorithm (SVR-GWO) (Balogun et al. 2021), adaptive neuro-fuzzy inference system model integrated with satin bowerbird optimizer algorithms (ANFIS-SBO) (Chen et al. 2021).However, these machine learning models still cannot avoid some disadvantages, for example, (1) It still requires a lot of prior knowledge and erection; (2) the existing network cannot fully extract potential landslide features; (3) the models are sensitive to missing data, and prone to fall into local optima.
Compared with machine learning, deep learning (DL) does not need to manually construct and select feature layers when dealing with object features, and at the same time, deep learning accepts a larger sample size, which is gradually applied in various fields.For example, Panahi (2020) used convolutional neural networks and recurrent neural networks to predict the probability of flash flood (Panahi et al. 2020); Kumar (2020) used deep learning model to complete the prediction of ground water depth (Kumar et al. 2020); Benzekri (2020) employed the deep learning model to construct an early forest fire detection system (Benzekri et al. 2020).In general, DL model performed a satisfactory ability of classification and regression.The main reason is that DL is completely a data-driven feature learning method, and has multi-level non-linear operations, which can abstractly represent classification features from a large amount of data, and combines gradient transfer method to optimize its end-to-end network structure (Zhu et al. 2020).Recently, more and more DL models have been successively applied in the field of LSA (Xiao et al. 2018;Huang et al. 2020;Li et al. 2021), among the different DL models, convolutional neural network (CNN) plays a significant role in landslide recognition and prediction, for example, Wang (2019) applied three novel CNN architectures in LSA, which achieved higher prediction accuracy than conventional methods (Wang et al. 2019); Sameen (2020) developed a DL-based technique for LSA through a 1-dimensional CNN, which performed better than ANN and SVM (Sameen et al. 2020); Fang et al. (2021) constructed four heterogeneous ensemble-learning techniques combining with CNN, RNN, SVM, and LR, the final results also showed that the DL model performed best (Fang et al. 2021).
However, these previous studies have some shortcomings, for example, in the process of building DL model, each convolutional layer will bring different prediction effects, the construction of multi-channel networks is ignored for landslide predisposing factor.Moreover, due to China's vast territory and diverse geological environment, there is still a lack of databases and studies related to landslides, and more reasonable, reliable and more accurate DL models still need to be explored.
Therefore, this study proposed a novel deep learning network named LSNet composed of multiple convolutional layers to predict the landslide susceptibility in Hanyin County, Shaanxi Province, China.The patches of landslide predisposing factor maps were used as the input data to train the LSNet, meanwhile the LSI was regarded as the output to predict the landslide susceptibility.In addition, the support vector machine model (SVM) and kernel logistic regression model (KLR) were employed to compare with LSNet.The primary difference here between this study and the literature mentioned is that approaches existed in this paper are seldom used and compared in landslide susceptibility assessment, especially LSNet and KLR.Another point is that the three models were first applied in Hanyin County and LSNet with multiple channels was established combined with the data structure representing the landslide, with the aim to improve the accuracy of LSA in the study area.Finally, all the results may help the government to make efficient decisions about landslide prevention and provide prevention references for landslide risk.

Sample description of study area
Hanyin County belongs to the hilly area in southern Shaanxi Province, the geographical coordinates are 32°68'-33°09'north latitude and 108°11'-108°44'east longitude (Fig. 1).The study area is about 51 km wide from east to west, 58 km long from north to south, and covers an area of about 1347 Km 2 .The climate type of study area is continental tropical monsoon climate and the temperature varies greatly.According to the local meteorological statistics, the mean annual precipitation in the past 50 years is about 920 mm, and the rainfall in the northern region is significantly less than that in the southern region.The water resources in the study area are very abundant, and there are 4 rivers in total, all of which belong to Yangtze River system.There are three types of groundwater in the study area, including loose rock pore water, carbonate fissure water, and bedrock fissure water.
The geomorphology of study area is dominated by low and middle mountains, with valleys, hills and basins, and the area of mountains accounts for 87%.The exposed strata and main lithology in the study area are shown in Table 1.Since the geotectonic location of the study area is located in the core zone of the Qinling microplate, there are many faults and folds in this area.In fact, there are a total of 5 faults that have been proven.Besides, according to the historical records, there have been 16 earthquakes in the study area, with an average magnitude of 4, but these earthquakes did not cause major damage.

Landslide inventory
Before carrying out the LSA, it is critical to verify about the information of landslides in the study area.Landslide inventory is to integration of landslide boundaries, locations, types and so on, which is the subsequent basis of data analysis and model construction.Based on the historical landslide data (SBGMR 1989;PRC 2020), remote sensing image (Cloud 2020), literatures (Liu and Huang 2006) and field survey, a total of 267 landslide locations were identified from 1989-2020 as the reference.These landslide locations were then imported into the GF-2 remote sensing images to delineated landslide boundary in ArcGIS software.In order to generate the landslide inventory map of study area (Fig. 1), all landslide boundaries were converted into polygons and resampled with the resolution of 30 m × 30 m. Figure 2 shows the preparation process of landslide inventory map.

Data preparation
In order to prepare the input dataset for model construction, 267 landslide samples were separated into two parts according to the ratio of 7/3 (Zhao and Chen 2020).Among them, 187 landslide samples were used as the training dataset to train the model, and the remaining 80 landslide samples were applied as the validation dataset to finish the validation purpose.

Analysis and quantification of landslide predisposing factors
There are no fixed guidelines for selecting the predisposing factors for LSA.After the comprehensive analysis of regional geo-environment characteristics and previous researches (Zhou and Fang 2015;Wang et al. 2020;Wu et al. 2020), we proposed altitude, slope angle, slope aspect, normalized difference vegetation index (NDVI), distance to rivers, distance to roads, distance to faults, mean annual precipitation (MAP) and lithology as the landslide predisposing factors.
Among them, altitude, slope angle, slope aspect are commonly used terrain factors, which can describe the influence of terrain on landslide with multi-dimension.For instance, altitude can affect the topographic properties, vegetation distribution and temperature differences at different heights.Slope angle is regarded as one major impact factor affecting instability and deformation of landslides, for example, the larger the slope angle, the less stable the slope is.
NDVI is a commonly used to indicate the vegetation coverage condition and it ranges from -1 to 1. Negative value expresses that there is water, snow and cloud  coverage.Positive value expresses that the earth coverage is vegetation, and the larger the coverage extent of the vegetation lead to higher NDVI positive value.At present, it is generally believed that the mechanism of landslide occurrence is the increase of pore water pressure in slope compositions by using the rainwater infiltration and accumulation during the rainy period.The high pore water pressure will lead to the increase of the effective load stress and the decrease of the shear strength of slope components, which are main triggering factors caused landslide.Therefore, MAP also occupies vital status of the landslide occurrence.
Rivers can affect the hydrogeology characteristics of slopes, and rivers usually corrode the toe of slope, which may decrease the anti-slide force.In mountainous area, it is common that numerous of landslide hazards are triggered by road constructions.Besides, in areas with frequent tectonic movements, the dislocation of faults can also cause landslide occurrence.Hence, distance to rivers, distance to roads, and distance to faults were regarded as predisposing factors in this study.
The lithology of the slope is the material foundation of the landslide.Some slopes are made up of hard rocks, some of which are made up of soft rocks, and some of them are made up of soil.Because of the difference between the lithologies, their shear strength varies.Often, slopes made of hard rocks do not easily fall, while slopes material of soft rocks or soil are easier to be destabilization deformation.
Since the original attribute data of each predisposing factor are very different, the frequency ratio (FR) is introduced to unify the dimension of each predisposing factor.The calculation process of FR value is shown in Eq. ( 1): where Sam ij stands for the percentage of landslides in each landslide predisposing factor class, and Are ij is the area percentage of each landslide predisposing factor class (Siahkamari et al. 2017).
Additionally, in order to calculate the FR value, it is necessary to classify the predisposing factors, and the data sources, resolution and classification method of each predisposing factor are listed in Table 2.

Methodologies
The main research contents include 4 parts: (1) using the data that are already available to complete the landslide inventory; (2) using FR value to quantify the landslide (1)

Factor optimization method
Since the assumption of machine learning modeling is that the variables are independent of each other, it needs to detect whether there is strong correlation between the factors.This strong correlation relationship is called multicollinearity which may cause the over-fitting or underfitting problems (Hong et al. 2018).In this study, the variance inflation factor (VIF) and tolerances (TOL) were applied to reflect the multicollinearity problem, which can be calculated by constructing a linear regression model based on the training dataset.When VIF > 10 and TOL < 0.1, it indicates that the predisposing factor has a multicollinearity problem and needs to be eliminated, vice versa (Pham et al. 2019).

Support vector machine model (SVM)
The basic principle of SVM is to search the optimal separating hyperplane that can maximize the interval between positive and negative samples in training dataset (Wang and Brenning 2021).Initially, SVM model was used as the supervised learning algorithm to solve binary classification problem, while the non-linear classification problem can be solved after introducing the kernel function.Therefore, the SVM model was applied in many researches about landside susceptibility assessment.In addition, there are three parameters namely penalty factor (C 0 ), non-sensitive loss function ( ε ), and kernel func- tion parameter ( γ ) that need to be adjusted appropriately in the process of constructing the SVM model (Xie et al. 2021).The main steps of SVM model construction can be described as below.
At first, the landslide predisposing factors are defined as the dataset of instance label pairs (s i , t i , i = 1, 2, …, n), where s i stands for the input data, t i is the output classes (landslide and non-landslide), and n is the number of training samples (Kumar et al. 2017).The training samples are mapped in to a n-dimensional hyperplane by using the RBF kernel function which can be defined as: Then mathematical expression of the n-dimensional hyperplane L needs to satisfy the following condition: where w denotes for the norm of normal hyperplane, and b is the constant.The maximum interval between vector and hyperplane can be derived by applying the (2) K s i , s j = −γ s i − s j , γ > 0.
(3) 3 The flowchart of the study Lagrangian multiplier (Abedini et al. 2019), and cost function can be expressed as:

Kernel logistic regression model (KLR)
In statistical learning, when there are phenomena such as non-linear estimation, non-normal estimation, and uneven variance, it may cause invalid estimation by using the ordinary regression method (Chen et al. 2018).These problems were overcome after the introduction of logistic regression, and logistic regression is widely used to solve binary classification problem.However, the structure of original logistic regression model is relatively simple, the flexibility is relatively low, and it still has defects in dealing with non-linear classification problems (Chen et al. 2019).While the kernel function can help to solve these problems effectively in constructing logistic regression model.Therefore, the hybrid model, namely kernel logistic regression is created.In order to be consistent with the SVM model above, the RBF kernel function is determined to build KLR model.The expression of KLR model is as follows: where p i is the probability of landslide occurrence, k i stands for the i th row of K(s i , s j ), and α is a constant for the intercept (Thai and Indra 2018).

Landslide net model (LSNet)
The deep learning has been widely used in the field of remote sensing image processing, including change detection, land use classification, image registration and so on.The deep belief networks, convolutional neural network (CNN), and auto coder are the three most commonly used network models in deep learning.The operating principle of these networks is to stack multiple layers within the model, and use the output of the previous item as the input of the next item, so that the features of each layer in the network can be converted into higher-dimensional features (Bui et al. 2020).Among them, the CNN has robust feature extraction capabilities and has been successfully applied in the field of image processing.
LSNet is a multi-layer feedforward neural network, the advantage of which is that it can process large-scale data in the form of multiple arrays from the local and global input data.The structure of LSNet is consist of multiple layers, which are related to each other through a set of learnable weights and biases.The local and global-scale (5) features can be captured by these convolutional blocks using scanning of the entire image.Meanwhile, the pooling layer and rectified linear unit (ReLU) layer are used for generalization to improve the non-linear fitting ability of the network (Li et al. 2021).Additionally, each convolutional layer contains feature maps obtained by multiple convolution kernels, and these feature maps share the node weights of the convolution kernels, so features can be extracted from different parts (Fig. 4).Specifically, the main operation performing in CNN can be generalized as follows: where O l−1 denotes the input feature map in lth layer, W l and b l , respectively, represent the weight and deviation of input feature layer convoluting by linear convolution, and σ is a non-linear function outside the convolution layer.

Result assessment method
In order to assess the accuracy of classification result and compare the performance of each model, statistical indexes are purposed to finish this work.A matrix (Table 3) is constructed by true positive (TP), false positive (FP), true negative (TN), and false negative (FN) calculating from training dataset (Pham et al. 2021).The accuracy and precision are calculated according to Eqs. ( 7) and (8) for accuracy assessment, meanwhile, the consistency of the results is verified with F1.The calculation process is as follows:

Model comparison method
In this study, the work of model comparison is purposed to carry out from three indicators including the running speed of the model, the classification ability for landslide and non-landslide, and the generalized performance of the model.Among them, based on the validation dataset, the running speed of the model is quantitative expressed by time, and the sensitivity and specificity are, respectively, used to reflect the classification ability for landslide and non-landslide (Eqs.( 9) and ( 10)) (Yanar et al. 2020).Additionally, the receiver operating characteristics curve ( 6) (ROC) is used for assessing the generalized performance, and in general, the larger the area under ROC curve (AUROC), the stronger the generalization ability of the model (Dang et al. 2020):

The quantification results of FR for landslide predisposing factors
In this study, the FR value was employed to quantify each landslide predisposing factor according to the classification result.It can be observed from Table 4 that the interval of Tertiary from the lithology factor has the highest FR value (FR = 2.32), followed by the range of < 100 from the distance to roads factor (FR = 1.85), and the range (10) of 278-548 from the altitude factor (FR = 1.81).On the contrary, the lowest FR value appears in both the 1432-2107 interval of the altitude factor (FR = 0.00) and the flat interval of the slope factor (FR = 0.00).

The optimization result of landslide predisposing factors
The VIF and TOL values of each landslide predisposing factor were calculated based on the quantified landslide predisposing factors, and the calculation results are shown in Table 5.As can be seen from the results, the largest VIF value and the smallest TOL value appear in NDVI (VIF = 1.433,TOL = 0.698), followed by the altitude (VIF = 1.293,TOL = 0.773) and the aspect (VIF = 1.268,TOL = 0.789).By contrast, the distance to roads has the smallest VIF value and the largest TOL value (VIF = 1.019,TOL = 0.981).Since the VIF and TOL values of all landslide predisposing factors are not inside the critical range (VIF > 10 and TOL < 0.1), all factors are retained and used to prepare the dataset.
Based on the optimized landslide predisposing factors, the training and validation datasets were prepared according to aforementioned partition principle.Subsequently, the training dataset was used as the input data to implement the following three models.

Implementation of SVM model
In this study, the training dataset was used to construct the SVM model.Since the parameters of RBF kernel function are significant for model construction, the   tenfold cross-validation method was used to search the most suitable parameter set (C 0 ,γ ).The optimized parameter set is (241, 0.02).Then run the trained SVM model in the python platform, and adjust the output range of the model to 0.000-1.000which also represents the LSI.In the end, the natural break (Jenks) method was used to divide the LSI into five ranges which, respectively, represent the very low susceptibility area (0.0899-0.2084), low susceptibility area (0.2085-0.4646), moderate susceptibility area (0.4647-0.6228), high susceptibility area (0.6229-0.7893) and very high susceptibility area (0.7894-0.9224), furthermore the LSM was generated by converting these areas to image in ArcGIS software (Fig. 5).

Implementation of KLR model
The construction progress of KLR model is similar to the SVM model.For the purpose of comparison, the parameter set (C 0 ,γ ) was consistent with that of the SVM model.Subsequently, the training dataset was used as the input data for KLR model construction in the python platform, and the output range of the LSI was adjusted to 0.000-1.000.Finally, the LSI was divided into five ranges by using the natural break (Jenks) method.These five ranges, respectively, represent the very low susceptibility area (0.0145-0.2459), low susceptibility area (0.2460-0.3695), moderate susceptibility area (0.3696-0.5161), high susceptibility area (0.5162-0.6974) and very high susceptibility area (0.6975-0.9983), moreover the LSM corresponding to KLR model was generated in Arc-GIS software (Fig. 6).

Implementation of LSNet model
The LSNet was coded using tensorflow 2.0 under the python environment, and running on a personal computer with Intel(R) Core(TM) i7-7700 k CPU, RTX 3080Ti GPU, 32 GB RAM, and the Windows 10 operating system.The LSNet had multi-layer structure, the size of input window was designed as 224 × 224.
In fact, many researchers in image processing have   AlexNet can implement more than 1000 categories of classification.In contrast, the landslide susceptibility mapping is a binary classification problem, which does not require deep network design.For this reason, in this study, the size of input layer in convolution kernel of LSNet was set to 5 × 5, the size of the convolution kernel for the other layers was set to 3 × 3, the number of feature maps for each layer was set to 64, 128, 128, 256, 256, respectively.At the same time, a pooling layer, non-linear activation function ReLU and batch normalized BN were set after each convolutional layer.Based on computer graphics vision, all other parameters of the LSNet were empirically optimized, for instance, the learning rate and epoch are set as 0.0001 and 600 to learn the depth features through back propagation.Subsequently, the number of neurons in the fully connected layers was set to 1024, 256, 128, 2, respectively (Table 6).Since the classical network, VGG and UNet have proved the validity of 3 × 3 convolution kernel and 64,128,256,1024 feature image parameters in various regression and classification tasks.Therefore, on the basis of these parameters, we design a network structure with multiple channel inputs and fuse them at the end of the network.In particular, we use the hyperparameter search method to obtain the fully connected layers and training parameters of the multi-channel fusion structure, such as learning rate, batch size and epoch.Finally, softmax was used to estimate the probability of landslide occurrence to output confidence, namely LSI.Similarly, the output range of LSI for LSNet was adjusted to 0.000-1.000,and, respectively, represents the very low susceptibility area (0.0045-0.2021), low susceptibility area (0.2022-0.3458), moderate susceptibility area (0.3459-0.4814), high susceptibility area (0.4815-0.8033) and very high susceptibility area (0.8034-0.9972) (Fig. 7).

The result of accuracy assessment
After mapping the LSMs of these three models, it is necessary to assess the quality of results.In this study, the matrix has been organized based on the validation dataset, then the accuracy, precision, and F1 values for each LSM were calculated (Table 7).As shown in Table 7, the LSNet model gets the highest accuracy value and precision value (accuracy = 0.950, precision = 0.951), by contrast, the SVM model gets the lowest accuracy value and precision value (accuracy = 0.825, precision = 0.850), while the performance of the KLR model is moderate.From the value of F1, the LSNet also gets the highest value (F1 = 0.951), followed by the KLR model and SVM model, which is also consistent with the ordering of accuracy and precision values.

The result of model comparison
In order to compare the running speed, classification and generalization performance, the run time, sensitivity, specificity and AUROC values were introduced to finish this work.As the results shown in Table 7, the largest sensitivity and specificity values belong to the LSNet model, indicating that the LSNet model has the best landslide and non-landslide classification abilities among these three models.On the contrary, the smallest sensitivity and specificity values belong to the SVM model, indicating that the landslide and non-landslide classification abilities of SVM model are the weakest among these three models.
For AUROC values (Fig. 8), the LSNet model also obtains the largest AUROC value (AUROC = 0.941), followed by the KLR model (AUROC = 0.899) and SVM model (AUROC = 0.835), and the results show that the LSNet model has the best generalization ability.

Discussion
In this paper, we show the progress and results of landslide susceptibility mapping based on SVM model, KLR model, and LSNet model in Hanyin County, Shaanxi Province, China.In terms of the model performance, although the classification accuracy of the three models is higher, the accuracy of LSNet and other statistical indexes are higher than that of SVM and KLR, which fully shows that the LSNet performs best in the study area.
Since both SVM and KLR are developed based on statistical theory, the quality of input data and the adjustment of model parameters in the process of model construction may affect the final result.Before preparation of input datasets, three classification methods, i.e., natural break (Jenks), equal interval, and custom interval were all used to grade FR-quantified landslide predisposing factors.However, the classification methods and results of landslide predisposing factors are inevitably affected by human factors, which may lead to over-fitting or under-fitting (Yacine and Pourghasemi 2019).For this reason, it is necessary to deeply analyze the impact of classification methods on data quality.Besides, this study only used two machine learning models for comparison, therefore, more models should be added for reference in subsequent research, so that the advantages and disadvantages of deep learning and machine learning in landslide susceptibility mapping can be more comprehensively compared.
In contrast, as a deep learning model, the input data of LSNet is a complete remote sensing image containing all the information.In order to distinguish landslide and non-landslide from image data, not only the objects in the image patch need to be characterized as landslides, but also need to accurately and reliably represent the contextual information of the landslide space background.The advantage of LSNet is to derive the category of the object at the image block level, and learn the spatial distribution through the CNN network with hierarchical representation, and finally obtain the probability of each object's category through multiple fully connected layers and softmax.It is different from machine learning in principle, and its specific advantages include: (1) LSNet can classify based on object blocks in a deep learning network of convolutional structure, and output the category probability; (2) LSNet uses the CNN model to learn the internal and overall spatial information of the object block to represent the contextual spatial semantic information of the category.LSNet represents the probability of the category at the object block level, which can avoid pixel-level misfits and improve the accuracy of classification (Dimililer et al. 2021).(3) LSNet can directly read remote sensing images without destroying the data structure, and can obtain richer and multi-source data sets, making the trained model robust.Interestingly, the running time of LSNet is significantly longer than that of SVM and KLR, which may be limited by the hardware performance of the computer, resulting in slower calculations.Nevertheless, this does not mean that the LSNet is not a state-of-the-art model and other studies have reached similar conclusions in their researches.However, since the occurrence of landslides is affected by multi-source factors, the characteristics of the landslides themselves are also very complex and cannot be described solely by human-represented features.Therefore, we try our best to use deep learning methods that are completely data-driven.In addition, in this study, the final network architecture of LSNet was determined by gradually increasing the ablation experiments, and at the same time, an LSNet with multiple channels was established combined with the data structure representing the landslide, which could well fulfill the requirements of LSA.Furthermore, because other network structures are generated based on their respective research objects and purposes, although we cannot conduct experiments on these networks one by one, we have designed a network structure that matches the confidence of the current landslide research.We have also demonstrated the effectiveness and accuracy of LSNet in our experiments.
On the other hand, as a black box model, DL cannot intuitively reflect the spatial distribution features of landslides in the study area during data preparation.On the contrary, in machine learning modeling, because FR is used to quantify the graded landslide predisposing factors, the spatial distribution of the landslide under the conditions of each predisposing factor can be intuitively reflected from the quantified results (Zhang et al. 2020).For instance, from the view of distance to rivers and roads, as the distance from roads and rivers increase, the FR value decreases, indicating that the closer to the river and the road, the more landslides are distributed.This is because the exposed rock and soil in study area have low mechanical strength, the surface is easily weathered and eroded, and the joints and fissures are very developed.Moreover, due to the scouring action from the river and excavation of the slope toe during road construction, the original stress structure of the slope was destroyed, which resulted in the instability of the slope and generated a large number of potential landslides.This consistent with the phenomenon we observed in the field, and is similar to the results of geological hazard studies in similar areas of the study area (Wang et al. 2016;Liu et al. 2020).

Conclusion
Landslide susceptibility mapping is a key step for landslide prevention work.This study used Hanyin County, Shaanxi Province, China as the study area to complete the work of landslide susceptibility mapping by building the LSNet model, SVM model, and KLR model, and generated the LSM.Then various statistical indexes were applied for the accuracy assessment, and the ROC curves was employed to compare the performance and classification ability of the models.In summary, the main conclusions are as follows: (1) in the process of dataset preparation and parameter adjustment, the machine learning model will inevitably be affected by human factors, resulting in unstable classification results.However, LSNet can overcome human interference and generate objective classification results.(2) LSNet can avoid the problems of over-fitting and under-fitting.The classification accuracy in the study area is high, moreover the generalization is stronger than the SVM model and the KLR model.The LSNet can be promoted and used in the study area.
This study introduced the construction method of LSNet model in detail, and compared the performance of LSNet model (deep learning), SVM model (machine learning), and KLR model (hybrid model), which enriched the landslide database of the study area and can provide reference for the application of deep learning model in landslide prevention in the future.Furthermore, the results of this study can improve the efficiency of landslide prevention for government decision-making in similar study areas, which is conducive to rapid response of landslide warning.Comprehensive risk assessors and land use planner can benefit from our study findings.Additionally, the proposed approach is an innovative method that may also help other scientists to develop landslide susceptibility maps in other areas, but also as an approach that could be used in geo-environmental problems besides natural hazard assessments.

Fig. 1
Fig. 1 The study area a location and landslide inventory map; b forest destruction caused by landslide; c road damage caused by landslides

Fig. 2
Fig.2The preparation process of landslide inventory map

Fig. 5
Fig. 5 Landslide susceptibility map of study area derived by SVM model

Fig. 6
Fig. 6 Landslide susceptibility map of study area derived by KLR model

Fig. 7
Fig. 7 Landslide susceptibility map of study area derived by LSNet model

Table 1
The main lithology information of the study area

Table 2
The information and data source of landslide predisposing factors

Table 3
Discriminant matrix of statistical indexes

Table 4
The FR calculation result for each class of landslide predisposing factors

Table 5
The VIF and TOL values of each landslide predisposing factor

Table 6
The architecture of LSNet

Table 7
Calculation results of statistical indexes for landslide susceptibility mapping The ROC curves of each landslide susceptibility model based on validation datasetLastly, we measured the running speed of each model, the results show that the running speed of the SVM model (32 s) and the KLR model (27 s) are relatively close, while the running speed of the LSNet model (107 s) is significantly slower than the first models.