Next Article in Journal
Snow Cover Reconstruction in the Brunswick Peninsula, Patagonia, Derived from a Combination of the Spectral Fusion, Mixture Analysis, and Temporal Interpolation of MODIS Data
Previous Article in Journal
Spatial Prediction of Fluvial Flood in High-Frequency Tropical Cyclone Area Using TensorFlow 1D-Convolution Neural Networks and Geospatial Data
Previous Article in Special Issue
Pointwise Modelling and Prediction for Ground Surface Uplifts in Abandoned Coal Mines from InSAR Observations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Landslide Dynamic Susceptibility Mapping Base on Machine Learning and the PS-InSAR Coupling Model

1
Faculty of Engineering, China University of Geosciences, Wuhan 430074, China
2
College of Artificial Intelligence, Southwest University, Chongqing 400715, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2023, 15(22), 5427; https://doi.org/10.3390/rs15225427
Submission received: 24 September 2023 / Revised: 15 November 2023 / Accepted: 18 November 2023 / Published: 20 November 2023
(This article belongs to the Special Issue Ground Deformation Source Modeling Using Remote Sensing Techniques)

Abstract

:
Complex and fragile geological conditions combined with periodic fluctuations in reservoir water levels have led to frequent landslide disasters in the Three Gorges Reservoir area. With the development of remote sensing technology, many scholars have applied it to landslide susceptibility assessment to improve model accuracy; however, how to couple these two to obtain the optimal susceptibility assessment model remains to be studied. Based on Sentinel-1 data, relevant data, and existing research results, the information value method (IV), random forest (RF), support vector machine (SVM), and convolutional neural network (CNN) models were selected to analyze landslide susceptibility in the urban area of Wanzhou. Models with superior performance will be coupled with PS-InSAR deformation data using two methods: joint training and weighted overlay. The accuracy of different models was assessed and compared with the aim of determining the optimal coupling model and the role of InSAR in the model. The results indicate that the accuracy of different landslide susceptibility prediction models is ranked as RF > SVM > CNN > IV. Among the coupled dynamic models, the performance ranking was as follows: InSAR jointly trained RF (IJRF) > InSAR weighted overlay RF (IWRF) > InSAR jointly trained SVM (IJSVM) > InSAR weighted overlay SVM (IWSVM). Notably, the IJRF model, which combines InSAR deformation data through joint training, exhibited the highest accuracy, with an AUC value of 0.995. In the factor importance analysis within the IJRF model, InSAR deformation data ranked third after hydrological distance (0.210) and elevation (0.163), with a value of 0.154. A comparison between landslide dynamic susceptibility mapping (LDSM) and landslide susceptibility mapping (LSM) revealed that the inclusion of InSAR deformation data effectively reduced false positives around the landslide areas. The results suggest that joint training is the most suitable coupling method, allowing for the optimal expression of InSAR deformation data and enhancing the predictive accuracy of the model. This study serves as a reference for future research and provides a foundation for landslide risk management.

1. Introduction

Landslides are a prevalent geologic disaster worldwide [1,2]. Reservoir water storage can significantly change the hydrogeological conditions of bank slopes, leading to landslides and bank slope deformation [3,4,5]. As one of the largest reservoirs, numerous landslides have developed in the Three Gorges Reservoir area [6], posing a serious threat to the life safety of local residents [7,8,9]. As an important means for landslide risk control, landslide susceptibility assessment analyzes the nonlinear relationship between landslides and related environmental factors [10], providing an accurate spatial probability distribution of landslide susceptible areas [11,12]. This assessment holds great significance for landslide disaster prevention and mitigation efforts.
The key to landslide susceptibility assessment lies in the construction of a susceptibility assessment model, the calculation of the spatial probability distribution of the landslide susceptible area [10], and the conduction of landslide susceptibility mapping (LSM). Three types of landslide susceptibility assessment models are commonly used in existing studies, including heuristic models, mathematical statistical models, and machine learning (ML) models. Heuristic models such as the analytic hierarchy process (AHP) [13] require researchers to judge the weights and thresholds of the factors based on their expertise [14]. However, this subjective experience leads to the uncertainty of the model [15]. Mathematical statistical models [16] including the information value method (IV) [17] and deterministic coefficient method, among others, rely on the engineering analogy method and superimpose factors in different ways to express the nonlinear relationship between factors and landslides. Machine learning models such as logistic regression [18], SVM [19], and RF [16] can efficiently capture the relationship between factors and landslides, which are widely used in landslide susceptibility assessment on account of their excellent performance and efficient modeling process [20], although partial models mat be challenging to interpret due to the black-box analysis process. Deep learning models, represented by CNN, being developed and derived from machine learning models [14], are able to mine deep features and relationships in data through their neural network structure, enabling them to address complex problems [21,22]. Although different models exhibit advantages in different aspects, there is lack of consensus on the optimal model for predicting landslide susceptibility [23]. To verify the applicability and performance of different models, four models were employed for landslide susceptibility prediction research in this study: IV from the mathematical statistical models, SVM and RF from the machine learning models, and CNN from the deep learning models.
The evolution of a landslide often initiates with the creep stage [24]; monitoring ground deformation can be used to determine the extent and stability of landslides. [3,25,26]. However, the current landslide susceptibility assessment typically constructs models through static factors without considering dynamic features such as ground deformation, which leads models to misjudge certain landslide areas [23]. With the development of remote sensing technology, the accurate and long-time surface deformation monitoring data from InSAR (interferometric synthetic aperture radar) [27,28,29] enables the construction of a landslide dynamic susceptibility model [30,31]. Utilizing the persistent scatterer interferometry (PSI) technique, Mishra and Jain [32] retrieved the displacements of the Baglihar Dam Reservoir slope. Hussain et al. [33] updated the landslide inventory for susceptibility mapping. Zhu et al. [34] constructed the landslide dynamic susceptibility assessment model by means of an empirical matrix. Liu et al. [35] and Cao et al. [36] used the weighted overlay method to construct a landslide dynamic susceptibility assessment model which achieved LDSM (landslide dynamic susceptibility mapping) by weighting the InSAR deformation data and LSM. Various landslide susceptibility assessment models and different InSAR coupling methods were used in these studies to prove the significance of InSAR deformation data in landslide susceptibility assessment. However, it remains uncertain which coupling method can be used to obtain the best landslide dynamic susceptibility assessment model.
In summary, this study aimed to compare the performance of different coupled models, find the optimal method for updating the LSM with InSAR deformation data, and obtain the best landslide dynamic susceptibility assessment model. The study in this paper was conducted in three stages: (1) Based on the landslide-related factor data in the study area, we constructed IV, SVM, RF, and CNN models to obtain the LSM; validation assessment and performance comparison of the selected models indicated that RF and SVM exhibit the best performance. (2) Using weighted overlay and joint training methods, the SVM and RF models are coupled with the InSAR deformation data to construct dynamic susceptibility models and generate the LDSM. (3) We compared the accuracy of the applied dynamic susceptibility models based on various statistical indicators, appraised the applicable coupling method, and analyzed the role of the InSAR deformation data in model through LDSM.

2. Study Area and Dataset

2.1. Study Area

The study area is in the urban area of Wanzhou District, within the Three Gorges Reservoir area. Its geographic coordinates range from 108°19′50″E to 108°31′10″E and from 30°43′20″N to 30°51′30″N (Figure 1), covering an approximate area of 289 km2. Wanzhou is situated in the Yangtze River valley of the eastern Sichuan Basin, characterized by a subtropical monsoon climate with features of a warm and humid climate, abundant rainfall, four distinct seasons, and notable vertical stratification [37]. It experiences an average annual temperature of 18.2 °C and an annual precipitation of 1155.8 mm. The region belongs to the Yangtze River system and features an intricate network of rivers with significant elevation differences and a dendritic distribution.
Within the study area, there is a distribution of erosional and depositional landforms along the Yangtze River and its tributaries, characterized by floodplains and terraces [38], in which the overall elevation is relatively low, with the lowest altitude at 149.2 m above sea level. The surrounding area consists of low mountains and hills, with elevations mostly in the range of 440~620 m and reaching a maximum of 825.2 m.
The outcropping strata in the study area belong to the monoclinic Middle Jurassic period and Quaternary [39]. The J2s lithology comprises sandstone, siltstone, silty claystone, and mudstone, accounting for approximately 73% of the total; the J3s lithology consists of purple–red mudstone, muddy siltstone, and brownish-red siltstone with varying thicknesses, interbedded with purple–gray fine-grained feldspathic sandstone, making up around 23% of the total area. Additionally, small quantities of J3p, J2xs, J2x are distributed in the eastern and southeastern parts of the study area, as shown in Figure 2f.

2.2. SAR Dataset

The Sentinel-1A SAR C-band data were obtained from the European Space Agency (ESA). The SAR dataset comprises 24 ascending scenes and each image covers an area of approximately 40 km × 85 km along the Wanzhou–Yunyang section of the Yangtze River (as shown in Figure 3), with a total coverage area of 3380.8 km2. The image acquisition time spans from 17 July 2019 to 12 June 2021. Detailed data acquisition dates are included in Table 1.
The SAR images are single-look complex data generated using dual HH+HV polarization and the interferometric wide (IW) swath mode, with an incidence angle of 34° and a spatial resolution of 5 × 20 m. Additionally, the 12.5 m DEM covering the study area was prepared for the PS processing.

2.3. Landslide-Related Factors and Dataset

This paper considers various conditional factors related to landslide occurrence combined with existing studies on landslide susceptibility and available data in the study area [6,38,39]. A basic dataset for landslide susceptibility assessment in the study area was constructed from aspects of stratum lithology, topography, and geological structure [40]. This dataset comprises eleven factors: elevation, slope, aspect, distance to roads, distance to river, lithology, plan curvature, profile curvature, NDVI, TWI, and land use type.
The digital elevation model (DEM) serves as the primary tool for representing topographic features. Topographic and geomorphic factors including elevation, slope, aspect, plan curvature, profile curvature, and TWI are derived from the DEM data. These factors provide valuable insights into the topographic relief, convexity, and other characteristics which are significant causative factors for landslide occurrence [14].
Water flow can accelerate soil erosion and weaken soil strength [41], while vehicular traffic on roads and surrounding structures can exert additional loads [42]. Collectively, these factors contribute to the risk of slope instability. Therefore, the distance to the river and the distance to roads were calculated by Euclidean distance tools in ArcGIS based on vector data of rivers and roads in the study area and subsequently incorporated into the dataset.
Geological lithology factors are extracted from geological maps of the study area. Using Landsat 5/8 remote sensing imagery from the Google Earth Engine (GEE) platform, spatial data with a resolution of 30 m for NDVI and land use types can be obtained. Vegetation loss due to human activities is one of the major factors causing shallow landslides [14]. NDVI and land use types provide insights into vegetation density and human activity in the study area [43].
Standardizing all factors to a consistent spatial scale, cropping them to the study area’s extent, and dividing them into grid of 30 × 30 m with a total of 434,005 uniform raster cells, we obtained a raster dataset for the study area which can be used for further landslide susceptibility analysis. During the training process for machine learning and deep learning, the dataset was divided: 70% was used for the training set and 30% was used for the validation set. A total of 10,000 positive and 20,000 negative samples were randomly selected from the training set and used for training landslide susceptibility assessment models so that the susceptibility prediction could be conducted for the whole study area.

3. Methodology

3.1. Multicollinearity Analysis

In landslide susceptibility assessment, multicollinearity of influencing factors will lead to redundancy and the mutual interference of inputs, which has a negative influence on the predictive ability of the model [36]. To assess the independence of factors, researchers often employ methods such as principal component analysis, the variance inflation factor (VIF), regression analysis, and the Pearson correlation coefficient [14,44,45,46]. In this study, VIF and TOL (tolerance) are used to measure the linear correlation of factors and can be calculated using Equation (1).
V I F = 1 1 R 2 = 1 T O L
Here, R2 represents the variance between influencing factors. When VIF > 5 or TOL < 0.2, the multicollinearity in the dataset will be considered.

3.2. Landslide Susceptibility Assessment Models

3.2.1. Information Value (IV) Model

The information value (IV) method is a mathematical statistical analysis method based on information theory and the engineering analogy method [16]. It was proposed by Claude Shannon in the 1940s during his study of information theory and has since been widely applied in the assessment of geohazard susceptibility. The occurrence of geohazards is influenced a multitude of factors, with each factor contributing to varying degrees [47]. The IV model constructs the landslide susceptibility assessment model by analyzing the distribution relationship between the known hazard points and factors, quantitatively calculating the weight of the factor on the occurrence of the hazard, and standardizes the scales of each factor [13,48,49,50]. The formula is as follows:
I = i = 1 n I i = i = 1 n ln N i / N S i / S
where: I represents the information value of the assessment raster cells; Ii is the information value provided by factor i for geohazards; Ni is the number of raster cells in which geohazards have occurred for influence factor i; N is the total number of geohazard raster cells in study area; Si is the total number of raster cells for influence factor i; and S is the total number of raster cells in the study area.

3.2.2. Random Forest (RF)

The random forest (RF) model is an integrated learning model based on decision trees [50]. It constructs multiple decision trees to make individual predictions and combines the results by a majority vote to arrive at the final prediction [49,51]. Leveraging bootstrap aggregating and the random subspace method, the RF method randomly samples training data for each decision tree and selects a subset of features for node splitting [14]. This random selection process can control both the bias and variance, reduce inter-feature correlations, and enhance the model’s robustness against outliers and noise [52]. Additionally, RF offers advantages such as not requiring dimensionality reduction for high-dimensional problems and the ability to evaluate the importance of each feature in classification tasks [53].

3.2.3. Support Vector Machine (SVM)

The support vector machine (SVM) model, introduced by Vapnik et al. [54], is a machine learning model that maps data into high-dimensional feature space using kernel functions and identifies the optimal hyperplane for data classification with the aim of maximizing the distance between data samples [14,19,24]. Based on support vector theory and structural risk minimization, it can be applied to address issues such as classification, regression, and pattern recognition [55]. The most generalized kernel function is the Gaussian kernel (also known as the radial basis function, RBF) [51]; based on the RBF, the final decision function can be composed as follows:
g x = s i g n i = 1 n γ i α i K x i , x j + b
where K(xi, xj) is the RBF kernel function. Since its inception, SVM has rapidly emerged as an indispensable tool in the field of machine learning; its exceptional classification capabilities, adaptability, and interpretability have led to its widespread application across various domains.

3.2.4. Convolutional Neural Network (CNN)

The convolutional neural network (CNN) model is a deep learning model initially proposed by LeCun for addressing image recognition tasks [56]. It is a NN (neural network) architecture model consisting of convolution layers, pooling layers, and fully connected layers [57]. By performing convolution operations with appropriately shaped kernels, transmitting and mapping data through the neurons, CNN can extract features from data progressively. The convolution algorithm is as follows:
C j = i = 0 N w j x b + b j               j = 1,2 , , k  
Here, k represents the number of convolution kernels; Cj represents the output of the j-th convolution kernel; i signifies the spatial position; xi signifies the input within the convolution window; wj refers the weights; bj refers the bias. CNNs gradually improve their model performance to fit training data and accomplish classification tasks through multiple iterations by calculating the loss function, backpropagating gradients, and updating weights.
In this study, a CNN model was constructed for landslide susceptibility assessment, comprising six convolution layers, six pooling layers, and two fully connected layers. Max-pooling layers were added to reduce the feature dimensions of the convolution layer outputs, aiding the model in learning different features and deep relationships from the input data [58]. Subsequently, the learned features were combined, transformed, and mapped to various categories through the fully connected layers. Adding the ReLU and tanh functions as activation functions for the convolutional and fully connected layers enables nonlinear transformations of the output results, thereby enhancing the model’s ability to express intricate relationships [59]. The ReLU and tanh functions are as follows:
f x = x       i f   x > 0 0       i f   x 0       max 0 , x
f x = t a n h x = e x e x e x + e x
Finally, the model combines the extracted features and maps them to the (0, 1) interval through the sigmoid activation function, generating probability values as the predictions results for landslide classification.

3.3. Model Accuracy Verification Methods

We employed ROC (receiver operating characteristic) curves, a confusion matrix, and the FR (frequency ratio) to quantitatively evaluate the accuracy and performance of models. The ROC curve is one of the most used accuracy evaluation methods in probabilistic models. The curve is constructed using specificity and sensitivity based on the model’s prediction results and evaluates the model accuracy through the AUC (area under the ROC curve) [60]. AUC values fall within the interval of [0.5, 1], and a higher AUC value approaching 1 indicates a better model accuracy, while a value near 0.5 suggests the model lacks predictive capability [61].
The confusion matrix categorizes binary classification results into four classes [36]: true positive (TP), true negative (TN), false positive (FP), and false negative (FN). Based on the parameters derived from the confusion matrix, the following four model performance evaluation metrics can be calculated [14]:
Accuracy measures the proportion of correctly classified samples in the total sample, with the formula:
A c c u r a c y = T P + T N T P + F P + T N + F N
Recall assesses the model’s ability to identify positive samples, represented by the proportion of correctly predicted positive samples out of all true positive samples, and can be calculated using:
R e c a l l = T P T P + F N
Precision is the proportion of correctly predicted positive samples out of all predicted positive samples, measuring the accuracy of the model’s predictions for positive samples. Its formula is as follows:
P r e c i s i o n = T P T P + F P
F1 score is a composite metric that considers both recall and precision, computed as the harmonic mean of both. It is used to comprehensively evaluated the performance of the classification model and its formula is defined as:
F 1   S c o r e = 2 × P r e c i s i o n R e c a l l P r e c i s i o n + R e c a l l
The frequency ratio (FR) represents the ratio of the landslide area in a specific attribute interval of a factor to the total area of that attribute interval within the entire study area [10,62]. In different susceptibility partitions generated by various models, calculating their corresponding FR values can indict the partition’s ability to express positive and negative landslide samples [13,49]. The formula is as follows:
F R = N i / N S i / S
FR values greater than 1 indicate that the corresponding attribute interval is more conducive for landslide occurrence and vice versa.

3.4. PS-InSAR

PS-InSAR employs interferometry synthetic aperture radar technology in conjunction with persistent scatterers (PS) to monitor subtle surface deformations [63,64]. Based on time-series InSAR data and the single master interferogram method, we conduct operations of registration, clipping, and noise reduction on multiple PS points within the image area of study, resulting in the generation of time-series deformation data for these PS points [65]. Subsequently, an interpolation method was applied to generate a ground deformation map (GDM) by integrating the PS points within the study area. This map provides precise measurements of surface subsidence, uplift, and displacement, thus serving as a valuable basis for the timely identification of potential geological hazards [66].

3.5. Dynamic Evaluation of Landslide Susceptibility

The application of InSAR technology in landslide susceptibility assessment enhances the consideration of dynamic factors such as ground uplift and subsidence, enabling the identification of potential landslide areas during the creep phase, leading to more accurate assessment results [67]. Common methods for integrating InSAR deformation data with landslide susceptibility models include joint training [23], weighted overlays [36], and constructing evaluation matrices [35,64]. In this study, we employed the first two coupling methods, selecting the top-performing two models to couple with InSAR deformation data in IV, SVM, RF, and CNN, resulting in dynamic landslide susceptibility assessment models. Subsequently, we calculated the probability distribution of landslide occurrence based on these coupled models, generated the LDSM, and conducted a comparative analysis of model performance.
The joint training method integrates InSAR deformation data as influencing factors into the base dataset. After undergoing the same preprocessing as other factors, InSAR deformation data is jointly trained in the model. This process extracted features from InSAR deformation data by leveraging the nonlinear function of models, revealing its relationship with landslide disasters and constructing a susceptibility assessment model that considers dynamic features.
The weighted overlay method considers both geological conditions and landslide evolution. It assigns weights to the GDM derived from InSAR data and the LSM obtained from landslide susceptibility assessment models after normalizing the absolute values of InSAR deformation data [36]. These two weighted components are overlaid to create a landslide dynamic susceptibility assessment model. The weighted overlay formula is as follows:
L D S M = 0.61 L S M + 0.39 G D M
The flow chart of the research process is shown in Figure 4.

4. Results

4.1. The Multicollinearity Analysis of Related Factors

Table 2 presents the results of multicollinearity testing among the 11 factors related to landslides. Among all the factors, elevation has the highest VIF value (3.011) and the lowest TOL value (0.332), while the distance to roads displays the lowest VIF value (1.080) and the highest TOL value (0.926). The results of the multicollinearity test show that the VIF and TOL values of selected factors are better than the critical value of insignificant multicollinearity (VIF > 5 and TOL < 0.2), indicating the absence of significant multicollinearity. This suggests that the selected factors have good independence and can be used for training models to achieve highly accurate predictive results.

4.2. Landslide Susceptibility Mapping

Based on the prepared dataset, the IV, SVM, RF, and CNN methods were applied to assess landslide susceptibility in the study area. Segmented into high (H), moderate (M), low (L), and very low (VL) landslide susceptibility using the natural breakpoint method, the probability predictions of selected models were visualized as LSM in ArcGIS (Figure 5). As depicted in Figure 5, the H susceptibility areas are primarily concentrated along the Yangtze River and its tributaries, especially in the sloped areas along the banks of the Yangtze River in Tangjiao Village and Tiancheng District, the regions along the bank of the Zhuxi River, and the valley areas at the origins of the Wuqiao River and Longbao River in the southern part of the study area. These areas are distributed along riverbanks and transitional zones between low hills and riverbanks. The erosion caused by fluctuations in reservoir water levels has diminished the stability of slopes in these areas. Furthermore, as these regions are interconnected with densely populated areas, frequent engineering activities and the depletion of vegetation exacerbate the susceptibility to landslides.
Among these models, the LSM of the IV model demonstrates a strong correlation with the factor of distance to rivers. All areas along the Yangtze River and its tributaries are identified as H susceptibility landslide areas. Overall, the susceptibility tends to decrease as the distance to rivers increases. While the IV model effectively partitions susceptibility areas to some extent, it struggles to distinguish varying landslide susceptibility when regions have similar distances to rivers, primarily due to the high weighting of the distance to rivers factor, which results from the predominant distribution of landslides along watercourses. Therefore, the IV model is considered to have the issue of insufficient differentiation. Table 3 illustrates the landslide grid ratios and information value in each class of IV model.
The susceptibility results of the other three prediction models showed that the areas with H susceptibility were consistent with the landslide distribution area, but there were still some differences in the details. The RF model achieved its best performance with parameters set at max_depth = 20, min_samples_leaf = 1, min_samples_split = 6, and n_estimators = 46. Among the factors, the distance to rivers carried the highest weight of 0.225, while lithology had the lowest weight of 0.028. The LSM generated based on RF-predicted probabilities shows adequate continuity and a reasonable distribution at different susceptibility levels. The SVM model performed best with parameters set at C = 1, kernel = ‘rbf’, and gamma = 0.8. The susceptibility partitions in the LSM produced by the SVM have a similar distribution to the RF model overall. However, the results of the SVM exhibit fragmentation and lower continuity. The CNN model achieved optimal results with parameters set at learning_rate = 0.01 and batch_size = 64. The LSM demonstrates moderate continuity, and areas closer to watercourses were accurately identified as high susceptibility areas. However, some landslide areas far from watercourses were not correctly labeled as high susceptibility areas.

4.3. Model Accuracy Verification

Various statistical model evaluation metrics enable a comprehensive assessment of model accuracy from different perspectives, facilitating a quantitative performance comparison. The ROC curve is a valuable tool for assessing model performance, and the AUC value can provide an intuitive representation of predictive accuracy. According to Figure 6, the AUC values for the IV, SVM, RF, and CNN models are 0.865, 0.953, 0.980, and 0.888, respectively. All models have AUC values above 0.8, indicating their competence in landslide prediction. Notably, the RF and SVM models outperform the others in terms of the AUC value, with the RF model achieving the highest AUC value of 0.980.
The confusion matrix is used to compare the classification results with the actual circumstances and measure the classification performance of models by deriving the parameters of accuracy, recall, precision and F1 score. Models with higher accuracy, F1 score, and a balance in recall and precision are typically regarded as performing better. Figure 7 portrays the accuracy, recall, precision, F1 score, and training speed of the selected models, respectively. For the sake of convenient comparison, the training speed is represented by the natural logarithm of the training time (except for the IV model). As shown in Figure 7, most models perform well in accuracy but show subpar performance in precision. The RF model exhibited the best performance metrics, with an accuracy of 0.976, a recall of 0.647, a precision of 0.625, an F1 score of 0.636, and the shortest training time compared to the other ML models. These indicators suggest that the model performance is ranked as: RF > SVM > CNN > IV.
FR is a reliable technique for quantifying the ability to express landslides within different susceptibility zones. To mitigate the differences in magnitudes and facilitate a more convenient comparison, natural logarithm processing was applied to FR, and the results are presented in Figure 8. Models with a higher FR in the H susceptibility zone and a lower FR in the M, L, and VL susceptibility zones are regarded as possessing enhanced classification capabilities. As shown in Figure 8, the prediction results of the SVM and RF exhibit higher landslide raster ratios in the H susceptibility zone and lower values in other susceptibility zones. In contrast, the IV and CNN models show lower landslide raster ratios in the H susceptibility zone but higher values in other susceptibility zones. The results indicate the superior ability of SVM and RF to distinguish landslide-prone areas.
Based on the comparative analysis conducted, it is evident that the SVM model and RF model exhibit the best performance in analyzing the relationship between landslide disasters and landslide-related factors, providing the most accurate landslide susceptibility predictions in the study area. Consequently, we have selected the SVM and RF models as the foundation for coupling with InSAR deformation data and conducting further analysis.

4.4. Result of PS-InSAR

Based on the SAR dataset, persistent scatterer interferometry (PSI) was conducted in ENVI 5.6. The scenes were paired to create 23 connection graphs, and the baseline is illustrated in Figure 9. Subsequently, interferometric processing, two inversions, and geocoding were performed, resulting in a deformation dataset containing 260,039 PS points. The parameters used in the PS process are detailed in Table 4.
After removing noisy PS points, we converted the PS points into a raster image in ArcGIS using their deformation rates as a field, producing the deformation rates in PS points for the study area (as shown in Figure 10a). Then, the “Kriging” tool was employed to interpolate the PS points, generating the deformation rates for the entire study area, as depicted in Figure 10b.
As shown in Figure 10, the areas adjacent to rivers in the study area exhibit an uplift trend, while the low hills and ridges regions further from rivers display subsidence trends. The maximum annual average uplift rate within the study area, measured at 40.18 mm/year, was observed in the urban area on the south bank of the Zhuxi River. Conversely, the highest subsidence rate, reaching 40.25 mm/year, was identified around Wanzhou Wuqiao Airport and the hilly region north of the Yangtze River.

4.5. Landslide Dynamic Susceptibility Mapping

The SVM and RF models were coupled with InSAR in two different methods, resulting in four landslide dynamic susceptibility assessment models: InSAR weighted overlay SVM (IWSVM), InSAR jointly trained SVM (IJSVM), InSAR weighted overlay RF (IWRF), and InSAR jointly trained RF (IJRF). The landslide probability predictions for each model were segmented through the natural breakpoint method in ArcGIS and visualized as LDSM (as shown in Figure 11). Overall, the LDSMs of the coupled models show a similar susceptibility distribution to the uncoupled models. With the introduction of the deformation characteristic factor, the jointly trained models demonstrate a decrease in M and L susceptibility areas, while VL susceptibility areas increase; in the weighted overlay models, the M and L susceptibility areas increase and VL susceptibility areas decrease.
The LDSM of the IJSVM model has partially alleviated the issue of fragmentation in susceptibility partitions. However, there is still some intermingling in certain H and M susceptibility areas, resulting in insufficient continuity. The IWSVM model predicts more areas as H susceptibility areas, with more dispersion-prone areas. The IJRF model achieved a more reasonable distribution of susceptibility partitions, accurately identifying landslide areas as H susceptibility areas. The susceptibility of non-landslide areas decreased after considering the deformation information feature, leading to a reduction in H, M, and L susceptibility areas and an increase in VL susceptibility areas. The prediction results of the IWRF model exhibit more M and L susceptibility areas and a concentration of H susceptibility areas around the watercourses.

5. Discussion

Landslide susceptibility mapping is a crucial method in the risk assessment and control of landslide disasters. However, due to the lack of dynamic landslide characteristics, which results in the absence of landslide movement state features, traditional models fail to provide reliable susceptibility assessment results. In our previous research, four models were applied for landslide susceptibility mapping in the study area, and a comparative analysis of these models was conducted using various evaluation parameters. Then, we used PS-InSAR technology to obtain the GDM of the study area and coupled the GDM with the SVM and RF models using two methods, resulting in four dynamic models and the generation of LDSMs for the study area. The difference between the LDSM and LSM indicates the effect of InSAR deformation data, demonstrating the feasibility of susceptibility assessment models that consider dynamic features. However, different coupling methods expressed the InSAR deformation data in separate ways. Therefore, we employed the same evaluation methods as before, validating the performance of these dynamic susceptibility models. Through further comparative analysis, we aimed to evaluate the performance of different dynamic models and assess the effect of InSAR deformation data in these models.

5.1. Performance Comparison of Landslide Dynamic Susceptibility Models

The original models and the four dynamic models before and after considering dynamic features were evaluated based on ROC curves and AUC values. The comparative results are illustrated in Figure 12. The AUC values for IJRF, IWRF, IJSVM, and IWSVM were 0.995, 0.973, 0.972, and 0.947. IJRF and IJSVM exhibited improvements of 1.53% and 1.99% compared with RF and SVM. Conversely, the AUC values of IWRF and IWSVM showed reductions of 0.71% and 0.63%. Among these models, the IJRF model achieved the best performance, with an AUC value of 0.995. The results of the ROC indicated that the coupling method of joint training enhances model accuracy, whereas the weighted overlay method diminishes model accuracy.
Similarly, the computed confusion matrix derived further parameters of dynamic susceptibility assessment models, as shown in Figure 13. It is evident from the figure that all models exhibit an accuracy exceeding 0.95, indicating their high precision. However, except for IJRF, the precision values of the models are significantly lower than their recall values, indicating that these models fail to strike a balance between recall and precision. This suggests these models could predict landslides but struggle in ensuring the accuracy of samples identified as landslides. In contrast, the IJRF model displays an accuracy, recall, precision, and F1 score of 0.991, 0.847, 0.885, and 0.865, respectively, and has the highest performance of all examined models in these terms.
FR values represent the landslide proportions in the VL, L, M, and H susceptibility levels for each model. These FR values are naturally logarithmically scaled and are presented in Figure 14. From the figure, it can be observed that SVM and its coupled models have higher landslide proportions in the VL and L susceptibility areas but lower proportions in the H susceptibility areas. That indicates its poor discrimination ability in distinguishing landslide susceptibility areas. On the other hand, the RF and its coupled models have both lower landslide proportions in the VL and L susceptibility areas and higher proportions in M and H partitions. Notably, the IJRF model, while having a lower proportion of landslide areas in the VL susceptibility area, has the highest proportion of landslide areas in the H susceptibility area, indicating the superiority of the IJRF model in discriminating landslide susceptibility areas from non-landslide susceptibility areas.

5.2. Effect Analysis of InSAR Deformation Data

The importance of landslide-related factors in the RF and IJRF models is depicted in Figure 15. Factor importance indicates the extent to which each factor influences the model’s classification performance, signifying how each factor helps the model explain the relationships of landslides and related factors. According to Figure 15, these two models share a similar distribution of factor importance, with the distance to rivers being the most crucial factor in both cases, while lithology, distance to roads, land use type have lower factor importance. After the inclusion of InSAR deformation data as a factor, the importance of other factors in the IJRF model slightly decreased, with InSAR deformation data having an importance score of 0.154, ranking third among the twelve factors. In the IJRF model, InSAR deformation data as a factor provides discriminative information between landslide and non-landslide samples, playing a significant role in the RF model’s classification decisions, thereby enhancing the model’s performance.
Furthermore, we discuss the impact of introducing InSAR deformation data by analyzing the changes in LSM within typical landslide areas. The Sifangbei landslide and Cizhuyuan landslide areas were chosen as representative cases; their LSMs generated by RF and GDM and LDSMs obtained through IJRF are depicted in Figure 16.
For the Sifangbei landslide area, situated on the northeast bank of the Yangtze River (a), the RF model identified the riverbank portion as a high susceptibility area for landslide disasters, but the LSM (b1) failed to accurately delineate the landslide boundaries. The GDM (b2) revealed substantial subsidence deformation in the Sifangbei landslide area compared to its surroundings. The inclusion of InSAR deformation data in the IJRF model led to improved differentiation between landslide and non-landslide areas along the riverbank area, as evidenced by the LDSM (b3). It effectively outlined the regions of the Sifangbei landslide while reducing the susceptibility of other surrounding non-landslide areas, enhancing the model’s performance.
The Cizhuyuan landslide area is located at the source of the Longbao River on the west bank of the Yangtze River (a). In this area, the LSM generated by the RF model (c1) predicted a significant number of false-positive areas within the H and M susceptibility zones. The introduction of the InSAR deformation factor (c2) reduced the H and M susceptibility areas for non-landslide areas in the LSM, resulting in a more accurate judgment of non-landslide regions in the LDSM (c3).
Overall, these comparisons demonstrate the enhanced performance of the IJRF model with the consideration of deformation features, which aids in more accurately identifying landslide areas and reducing susceptibility ratings for non-landslide regions around the landslide susceptible area.

6. Conclusions

Remote sensing images and the time-series InSAR technique provide critical data for the early identification and stability assessment of landslides [31,68,69], but are often underutilized in existing landslide susceptibility evaluation studies. Focusing on the Wanzhou urban area within the Three Gorges Reservoir area (TGRA), this paper presented landslide susceptibility mapping (LSM) and time series deformation analysis based on geological and Sentinel-1 data and explored the landslide dynamic susceptibility mapping (LDSM) method coupled with machine learning and the PS-InSAR model in different approaches.
LSM shows that high landslide susceptibility areas in the Wanzhou urban area are primarily concentrated along the riverside. The SVM and RF show superior performance in various evaluation indicators than other models (IV, CNN). The subsidence of Wanzhou urban area obtained using InSAR techniques ranges from −40.25 to 40.18 mm/year, with some deformed regions correlated with potential landslides. After incorporating InSAR deformation data as a training factor, the results of the IJRF model exhibited fewer false positives and a more precise delineation of landslide areas. Additionally, the importance of factors indicated the critical role of InSAR deformation data in the IJRF model.
Summarily, the integration of time series InSAR analysis with the RF machine learning model effectively enhanced the LSM model’s precision, contributing to disaster prevention and mitigation in the TGRA. However, the C-band (wavelength 5.6 cm) data used in this study was prone to decorrelation with the presence of vegetation, leading to inaccurate estimations of ground deformation [32]. Therefore, longer-band, multi-spatial view SAR data and multi-temporal InSAR (MT-InSAR) methods will be considered to obtain comprehensive ground deformation data and assist in landslide susceptibility mapping.

Author Contributions

Writing—original draft, F.M.; Data curation, Q.R.; Writing—review and editing, Y.W. and Z.Q. (Zhao Qian); Methodology, Z.K.; Investigation, Z.Q. (Zhangkui Qin). All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the National Key R&D Program of China (2023YFC3007201), National Natural Science Foundation of China (42377161, 41977244), Natural Science Foundation of Hubei Province (2023AFB580), and Guizhou Provincial Science and Technology Project (QKHZC[2023]YB127). The authors thank the colleagues in our laboratory for their constructive comments and assistance.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

This study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jari, A.; Khaddari, A.; Hajaj, S.; Bachaoui, E.M.; Mohammedi, S.; Jellouli, A.; Mosaid, H.; El Harti, A.; Barakat, A. Landslide Susceptibility Mapping Using Multi-Criteria Decision-Making (MCDM), Statistical, and Machine Learning Models in the Aube Department, France. Earth 2023, 4, 698–713. [Google Scholar] [CrossRef]
  2. Chen, C.; Shen, Z.; Weng, Y.; You, S.; Lin, J.; Li, S.; Wang, K. Modeling Landslide Susceptibility in Forest-Covered Areas in Lin’an, China, Using Logistical Regression, a Decision Tree, and Random Forests. Remote Sens. 2023, 15, 4378. [Google Scholar] [CrossRef]
  3. Macciotta, R.; Hendry, M.T. Remote Sensing Applications for Landslide Monitoring and Investigation in Western Canada. Remote Sens. 2021, 13, 366. [Google Scholar] [CrossRef]
  4. Michoud, C.; Baumann, V.; Lauknes, T.R.; Penna, I.; Derron, M.-H.; Jaboyedoff, M. Large slope deformations detection and monitoring along shores of the Potrerillos dam reservoir, Argentina, based on a small-baseline InSAR approach. Landslides 2016, 13, 451–465. [Google Scholar] [CrossRef]
  5. Urgilez Vinueza, A.; Handwerger, A.L.; Bakker, M.; Bogaard, T. A new method to detect changes in displacement rates of slow-moving landslides using InSAR time series. Landslides 2022, 19, 2233–2247. [Google Scholar] [CrossRef]
  6. Miao, F.; Zhao, F.; Wu, Y.; Li, L.; Török, Á. Landslide susceptibility mapping in Three Gorges Reservoir area based on GIS and boosting decision tree model. Stoch. Environ. Res. Risk Assess. 2023, 37, 2283–2303. [Google Scholar] [CrossRef]
  7. Guo, L.; Miao, F.; Zhao, F.; Wu, Y. Data mining technology for the identification and threshold of governing factors of landslide in the Three Gorges Reservoir area. Stoch. Environ. Res. Risk Assess. 2022, 36, 3997–4012. [Google Scholar] [CrossRef]
  8. Ma, J.; Wang, Y.; Niu, X.; Jiang, S.; Liu, Z. A comparative study of mutual information-based input variable selection strategies for the displacement prediction of seepage-driven landslides using optimized support vector regression. Stoch. Environ. Res. Risk Assess. 2022, 36, 3109–3129. [Google Scholar] [CrossRef]
  9. Tang, H.; Wasowski, J.; Hsein, J.C. Geohazards in the three Gorges Reservoir Area, China—Lessons learned from decades of research. Eng. Geol. 2019, 261, 105267. [Google Scholar] [CrossRef]
  10. Huang, H.; Yan, J.; Fan, X.; Yao, C.; Huang, J.; Chen, W.; Hong, H. Uncertainty pattern in landslide susceptibility prediction modelling: Effects of different landslide boundaries and spatial shape expressions. Geosci. Front. 2022, 13, 101317. [Google Scholar] [CrossRef]
  11. Chen, W.; Zhang, S.; Li, R.; Shahabi, H. Performance evaluation of the GIS-based data mining techniques of best-first decision tree, random forest, and naïve Bayes tree for landslide susceptibility modeling. Sci. Total Environ. 2018, 644, 1006–1018. [Google Scholar] [CrossRef] [PubMed]
  12. Huang, F.; Cao, Z.; Jiang, S.H.; Zhou, C.; Huang, J.; Guo, Z. Landslide susceptibility prediction based on a semi-supervised multiple-layer perceptron model. Landslides 2022, 17, 2919–2930. [Google Scholar] [CrossRef]
  13. Wubalem, A.; Tesfaw, G.; Dawit, Z.; Getahun, B.; Mekuria, T.; Jothimani, M. Comparison of statistical and analytical hierarchy process methods on flood susceptibility mapping: In a case study of the Lake Tana sub-basin in northwestern Ethiopia. Open Geosci. 2021, 13, 1668–1688. [Google Scholar] [CrossRef]
  14. Saha, S.; Majumdar, P.; Bera, B. Deep learning and benchmark machine learning based landslide susceptibility investigation, Garhwal Himalaya (India). Quat. Sci. Adv. 2023, 10, 100075. [Google Scholar] [CrossRef]
  15. Kavzoglu, T.; Sahin, E.K.; Colkesen, I. Landslide susceptibility mapping using GIS-based multi-criteria decision analysis, support vector machines, and logistic regression. Landslides 2014, 11, 425–439. [Google Scholar] [CrossRef]
  16. Polat, A. An innovative, fast method for landslide susceptibility mapping using GIS-based LSAT toolbox. Environ. Earth Sci. 2021, 80, 217. [Google Scholar] [CrossRef]
  17. Wang, Z.; Ma, C.; Qiu, Y.; Xiong, H.; Li, M. Refined Zoning of Landslide Susceptibility: A Case Study in Enshi County, Hubei, China. Int. J. Environ. Res. Public Health 2022, 19, 9412. [Google Scholar] [CrossRef] [PubMed]
  18. Broeckx, J.; Vanmaercke, V.; Duchateau, R.; Poesen, J. A data-based landslide susceptibility map of Africa. Earth-Sci. Rev. 2018, 185, 102–121. [Google Scholar] [CrossRef]
  19. Yu, C.; Chen, J. Landslide Susceptibility Mapping Using the Slope Unit for Southeastern Helong City, Jilin Province, China: A Comparison of ANN and SVM. Symmetry 2020, 12, 1047. [Google Scholar] [CrossRef]
  20. Kavzoglu, T.; Colkesen, I.; Sahin, E.K. Machine Learning Techniques in Landslide Susceptibility Mapping: A Survey and a Case Study. In Landslides: Theory, Practice and Modelling; Advances in Natural and Technological Hazards Research; Springer: Cham, Swizerland, 2019; p. 50. [Google Scholar] [CrossRef]
  21. Saha, S.; Saha, A.; Hembram, T.K.; Mandal, K.; Sarkar, R.; Bhardwaj, D. Prediction of spatial landslide susceptibility applying the novel ensembles of CNN, GLM and random forest in the Indian Himalayan region. Stoch. Environ. Res. Risk Assess. 2022, 36, 3597–3616. [Google Scholar] [CrossRef]
  22. Nikoobakht, S.; Azarafza, M.; Akgün, H.; Derakhshani, R. Landslide Susceptibility Assessment by Using Convolutional Neural Network. Appl. Sci. 2022, 12, 5992. [Google Scholar] [CrossRef]
  23. Yao, J.; Yao, X.; Zhao, Z.; Liu, X. Performance comparison of landslide susceptibility mapping under multiple machine-learning based models considering InSAR deformation: A case study of the upper Jinsha River. Nat. Hazards Risk 2023, 14, 1. [Google Scholar] [CrossRef]
  24. Miao, F.; Wu, Y.; Xie, Y.; Li, Y. Prediction of landslide displacement with step-like behavior based on multialgorithm optimization and a support vector regression model. Landslides 2018, 15, 475–488. [Google Scholar] [CrossRef]
  25. Reyes-Carmona, C.; Galve, J.P.; Moreno-Sánchez, M.; Riquelme, A.; Ruano, P.; Millares, A.; Teixidó, T.; Sarro, R.; Pérez-Peña, J.V.; Barra, A.; et al. Rapid characterisation of the extremely large landslide threatening the Rules Reservoir (Southern Spain). Landslides 2021, 18, 3781–3798. [Google Scholar] [CrossRef]
  26. Li, Y.; Wang, X.; Mao, H. Influence of human activity on landslide susceptibility development in the Three Gorges area. Nat. Hazards 2020, 104, 2115–2151. [Google Scholar] [CrossRef]
  27. Shi, X.; Chen, C.; Dai, K.; Deng, J.; Wen, N.; Yin, Y.; Dong, X. Monitoring and Predicting the Subsidence of Dalian Jinzhou Bay International Airport, China by Integrating InSAR Observation and Terzaghi Consolidation Theory. Remote Sens. 2022, 14, 2332. [Google Scholar] [CrossRef]
  28. Xie, M.; Zhao, W.; Ju, N.; He, C.; Huang, H.; Cui, Q. Landslide evolution assessment based on InSAR and real-time monitoring of a large reactivated landslide, Wenchuan, China. Eng. Geol. 2020, 277, 105781. [Google Scholar] [CrossRef]
  29. Dai, K.; Li, Z.; Xu, Q. Entering the Era of Earth Observation-Based Landslide Warning Systems: A Novel and Exciting Framework. Geosci. Remote Sens. 2020, 8, 136–153. [Google Scholar] [CrossRef]
  30. Dai, K.; Chen, C.; Shi, X.; Wu, M.; Feng, W.; Xu, Q.; Liang, R.; Zhuo, G.; Li, Z. Dynamic landslides susceptibility evaluation in Baihetan Dam area during extensive impoundment by integrating geological model and InSAR observations. Int. J. Appl. Earth Obs. Geoinf. 2023, 116, 103157. [Google Scholar] [CrossRef]
  31. Zhou, C.; Cao, Y.; Hu, X.; Yin, K.; Wang, Y.; Catani, F. Enhanced dynamic landslide hazard mapping using MT-InSAR method in the Three Gorges Reservoir Area. Landslide 2022, 19, 1585–1597. [Google Scholar] [CrossRef]
  32. Mishra, V.; Jain, K. Satellite based assessment of artificial reservoir induced landslides in data scarce environment: A case study of Baglihar reservoir in India. J. Appl. Geophys. 2022, 205, 104754. [Google Scholar] [CrossRef]
  33. Hussain, S.; Pan, B.; Afzal, Z.; Ali, M.; Zhang, X.; Shi, X. Landslide detection and inventory updating using the time-series InSAR approach along the Karakoram Highway, Northern Pakistan. Sci. Rep. 2023, 13, 7485. [Google Scholar] [CrossRef] [PubMed]
  34. Zhu, Z.; Gan, S.; Yuan, X.; Zhang, J. Landslide Susceptibility Mapping with Integrated SBAS-InSAR Technique: A Case Study of Dongchuan District, Yunnan (China). Sensors 2022, 22, 5587. [Google Scholar] [CrossRef] [PubMed]
  35. Liu, W.; Zhang, Y.; Liang, Y.; Sun, P.; Li, Y.; Su, X.; Wang, A.; Meng, X. Landslide Risk Assessment Using a Combined Approach Based on InSAR and Random Forest. Remote Sens. 2022, 14, 2131. [Google Scholar] [CrossRef]
  36. Cao, C.; Zhu, K.; Xu, P.; Shan, B.; Yang, G.; Song, S. Refined landslide susceptibility analysis based on InSAR technology and UAV multi-source data. J. Clean. Prod. 2022, 368, 133146. [Google Scholar] [CrossRef]
  37. Yu, X.; Wang, Y.; Niu, R.; Hu, Y. A Combination of Geographically Weighted Regression, Particle Swarm Optimization and Support Vector Machine for Landslide Susceptibility Mapping: A Case Study at Wanzhou in the Three Gorges Area, China. Int. J. Environ. Res. Public Health 2016, 13, 487. [Google Scholar] [CrossRef]
  38. Zhang, H.; Song, Y.; Xu, S.; He, T.; Li, Z.; Yu, X.; Liang, Y.; Wu, W.; Wang, Y. Combining a class-weighted algorithm and machine learning models in landslide susceptibility mapping: A case study of Wanzhou section of the Three Gorges Reservoir, China. Comput. Geosci. 2022, 158, 104966. [Google Scholar] [CrossRef]
  39. Song, Y.; Niu, R.; Xu, S.; Ye, R.; Peng, L.; Guo, T.; Li, S.; Chen, T. Landslide Susceptibility Mapping Based on Weighted Gradient Boosting Decision Tree in Wanzhou Section of the Three Gorges Reservoir Area (China). ISPRS Int. J. Geo-Inf. 2019, 8, 4. [Google Scholar] [CrossRef]
  40. Cheng, J.; Dai, X.; Wang, Z.; Li, J.; Qu, G.; Li, W.; She, J.; Wang, Y. Landslide Susceptibility Assessment Model Construction Using Typical Machine Learning for the Three Gorges Reservoir Area in China. Remote Sens. 2022, 14, 2257. [Google Scholar] [CrossRef]
  41. Pradhan, A.M.S.; Kim, Y.-T. Rainfall-Induced Shallow Landslide Susceptibility Mapping at Two Adjacent Catchments Using Advanced Machine Learning Algorithms. ISPRS Int. J. Geo-Inf. 2020, 9, 569. [Google Scholar] [CrossRef]
  42. Hua, Y.; Wang, X.; Li, Y.; Xu, P.; Xia, W. Dynamic development of landslide susceptibility based on slope unit and deep neural networks. Landslides 2021, 18, 281–302. [Google Scholar] [CrossRef]
  43. Glade, T. Landslide occurrence as a response to land use change: A review of evidence from New Zealand. CATENA 2003, 51, 297–314. [Google Scholar] [CrossRef]
  44. Sun, X.; Chen, J.; Li, Y.; Rene, N.N. Landslide Susceptibility Mapping along a Rapidly Uplifting River Valley of the Upper Jinsha River, Southeastern Tibetan Plateau, China. Remote Sens. 2022, 14, 1730. [Google Scholar] [CrossRef]
  45. Arabameri, A.; Karimi-Sangchini, E.; Pal, S.C.; Saha, A.; Chowdhuri, I.; Lee, S.; Tien Bui, D. Novel Credal Decision Tree-Based Ensemble Approaches for Predicting the Landslide Susceptibility. Remote Sens. 2020, 12, 3389. [Google Scholar] [CrossRef]
  46. Roy, P.; Chakrabortty, R.; Chowdhuri, I.; Malik, S.; Das, B.; Pal, S.C. Development of Different Machine Learning Ensemble Classifier for Gully Erosion Susceptibility in Gandheswari Watershed of West Bengal, India. In Machine Learning for Intelligent Decision Science; Algorithms for Intelligent Systems; Springer: Singapore, 2020. [Google Scholar] [CrossRef]
  47. Miao, F.; Zhao, F.; Wu, Y.; Li, L.; Xue, Y.; Meng, J. A novel seepage device and ring-shear test on slip zone soils of landslide in the Three Gorges Reservoir Area. Eng. Geol. 2022, 307, 106779. [Google Scholar] [CrossRef]
  48. Rahman, G.; Bacha, A.S.; Ul Moazzam, M.F.; Rahman, A.U.; Mahmood, S.; Almohamad, H.; Al Dughairi, A.A.; Al-Mutiry, M.; Alrasheedi, M.; Abdo, H.G. Assessment of landslide susceptibility, exposure, vulnerability, and risk in shahpur valley, eastern hindu kush. Front. Earth Sci. 2022, 10, 953627. [Google Scholar] [CrossRef]
  49. Islam, F.; Riaz, S.; Ghaffar, B.; Tariq, A.; Shah, S.U.; Nawaz, M.; Hussain, M.L.; Amin, N.U.; Li, Q.; Lu, L.; et al. Landslide susceptibility mapping (LSM) of Swat District, Hindu Kush Himalayan region of Pakistan, using GIS-based bivariate modeling. Front. Environ. Sci. 2022, 10, 1027423. [Google Scholar] [CrossRef]
  50. Yu, C.; Chen, J. Application of a GIS-Based Slope Unit Method for Landslide Susceptibility Mapping in Helong City: Comparative Assessment of ICM, AHP, and RF Model. Symmetry 2020, 12, 1848. [Google Scholar] [CrossRef]
  51. Fu, Z.; Li, C.; Yao, W. Landslide susceptibility assessment through TrAdaBoost transfer learning models using two landslide inventories. CATENA 2023, 222, 106799. [Google Scholar] [CrossRef]
  52. Dou, J.; Yunus, A.P.; Bui, D.T.; Merghadi, A.; Sahana, M.; Zhu, Z.; Chen, C.-W.; Han, Z.; Pham, B.T. Improved landslide assessment using support vector machine with bagging, boosting, and stacking ensemble machine learning framework in a mountainous watershed, Japan. Landslides 2020, 17, 641–658. [Google Scholar] [CrossRef]
  53. Xing, Y.; Huang, S.; Yue, J.; Chen, Y.; Xie, W.; Wang, P.; Xiang, Y.; Peng, Y. Patterns of influence of different landslide boundaries and their spatial shapes on the uncertainty of landslide susceptibility prediction. Nat. Hazards 2023, 118, 709–727. [Google Scholar] [CrossRef]
  54. Vapnik, V. Support-vector networks. Mach. Learn. 1995, 20, 273–297. [Google Scholar]
  55. Yu, H.; Lu, Z. Review on landslide susceptibility mapping using support vector machines. CATENA 2018, 165, 520–529. [Google Scholar] [CrossRef]
  56. LeCun, Y.; Bengio, Y. Convolutional networks for images, speech, and time series. In The Handbook of Brain Theory and Neural Networks; MIT Press: Boston, MA, USA, 1995; Volume 3361, pp. 255–258. [Google Scholar]
  57. Thi Ngo, P.T.; Panahi, M.; Khosravi, K.; Ghorbanzadeh, O.; Kariminejad, N.; Cerda, A.; Lee, S. Evaluation of deep learning algorithms for national scale landslide susceptibility mapping of Iran. Geosci. Front. 2021, 12, 505–519. [Google Scholar] [CrossRef]
  58. Hakim, W.; Rezaie, F.; Nur, A.; Panahi, M.; Khosravi, K.; Lee, G.; Lee, S. Convolutional neural network (CNN) with metaheuristic optimization algorithms for landslide susceptibility mapping in Icheon, South Korea. J. Environ. Manag. 2022, 305, 114367. [Google Scholar] [CrossRef]
  59. Saha, S.; Bhattacharjee, S.; Shit, P.K.; Sengupta, N.; Bera, B. Deforestation probability assessment using integrated machine learning algorithms of Eastern Himalayan foothills (India). Resour. Conserv. Recycl. Adv. 2022, 14, 200077. [Google Scholar] [CrossRef]
  60. Bianchini, S.; Solari, L.; Del Soldato, M.; Raspini, F.; Montalti, R.; Ciampalini, A.; Casagli, N. Ground Subsidence Susceptibility (GSS) Mapping in Grosseto Plain (Tuscany, Italy) Based on Satellite InSAR Data Using Frequency Ratio and Fuzzy Logic. Remote Sens. 2019, 11, 2015. [Google Scholar] [CrossRef]
  61. Merghadi, A.; Yunus, A.P.; Dou, J.; Whiteley, J.; ThaiPham, B.; Bui, D.; Avtar, R.; Abderrahmane, B. Machine learning methods for landslide susceptibility studies: A comparative overview of algorithm performance. Earth-Sci. Rev. 2020, 207, 103225. [Google Scholar] [CrossRef]
  62. Guo, Z.; Shi, Y.; Huang, F.; Fan, X.; Huang, J. Landslide susceptibility zonation method based on C5.0 decision tree and K-means cluster algorithms to improve the efficiency of risk management. Geosci. Front. 2021, 12, 101249. [Google Scholar] [CrossRef]
  63. Shen, C.; Feng, Z.; Xie, C.; Fang, H.; Zhao, B.; Ou, W.; Zhu, Y.; Wang, K.; Li, H.; Bai, H.; et al. Refinement of Landslide Susceptibility Map Using Persistent Scatterer Interferometry in Areas of Intense Mining Activities in the Karst Region of Southwest China. Remote Sens. 2019, 11, 2821. [Google Scholar] [CrossRef]
  64. Ciampalini, A.; Raspini, F.; Lagomarsino, D.; Catani, F.; Casagli, N. Landslide susceptibility map refinement using PSInSAR data. Remote Sens. Environ. 2016, 184, 302–315. [Google Scholar] [CrossRef]
  65. Hakim, W.; Fadhillah, M.; Park, S.; Pradhan, B.; Won, J.; Lee, C. InSAR time-series analysis and susceptibility mapping for land subsidence in Semarang, Indonesia using convolutional neural network and support vector regression. Remote Sens. Environ. 2023, 287, 113453. [Google Scholar] [CrossRef]
  66. Zhang, G.; Wang, S.; Chen, Z.; Liu, Y.; Xu, Z.; Zhao, R. Landslide susceptibility evaluation integrating weight of evidence model and InSAR results, west of Hubei Province, China. Egypt. J. Remote Sens. Space Sci. 2023, 26, 95–106. [Google Scholar] [CrossRef]
  67. Fobert, M.-A.; Singhroy, V.; Spray, J.G. InSAR Monitoring of Landslide Activity in Dominica. Remote Sens. 2021, 13, 815. [Google Scholar] [CrossRef]
  68. Zhou, C.; Cao, Y.; Yin, K.; Intrieri, E.; Catani, F.; Wu, L. Characteristic comparison of seepage-driven and buoyancy-driven landslides in Three Gorges Reservoir area, China. Eng. Geol. 2022, 301, 106590. [Google Scholar] [CrossRef]
  69. Zhou, C.; Cao, Y.; Yin, K.; Wang, Y.; Shi, X.; Catani, F.; Ahmed, B. Landslide Characterization Applying Sentinel-1 Images and InSAR Technique: The Muyubao Landslide in the Three Gorges Reservoir Area, China. Remote Sens. 2020, 12, 3385. [Google Scholar] [CrossRef]
Figure 1. (a) Geographical location of the Three Gorges Reservoir area; (b) Geological map and geological hazards (red dots) in the Three Gorges Reservoir area; (c) DEM of study area.
Figure 1. (a) Geographical location of the Three Gorges Reservoir area; (b) Geological map and geological hazards (red dots) in the Three Gorges Reservoir area; (c) DEM of study area.
Remotesensing 15 05427 g001
Figure 2. Landslide-related factors in the study area: (a) elevation, (b) slope, (c) aspect, (d) distance to roads, (e) distance to waters, (f) lithology, (g) plan curvature, (h) profile curvature, (i) land use type, (j) NDVI, and (k) TWI.
Figure 2. Landslide-related factors in the study area: (a) elevation, (b) slope, (c) aspect, (d) distance to roads, (e) distance to waters, (f) lithology, (g) plan curvature, (h) profile curvature, (i) land use type, (j) NDVI, and (k) TWI.
Remotesensing 15 05427 g002aRemotesensing 15 05427 g002b
Figure 3. Sentinel-1A SAR data coverage and study area wrapped on the 12.5 m DEM used for the PS processing.
Figure 3. Sentinel-1A SAR data coverage and study area wrapped on the 12.5 m DEM used for the PS processing.
Remotesensing 15 05427 g003
Figure 4. Flow chart of the proposed method.
Figure 4. Flow chart of the proposed method.
Remotesensing 15 05427 g004
Figure 5. Landslide susceptibility maps: (a) LSM by IV, (b) LSM by RF, (c) LSM by SVM, and (d) LSM by CNN.
Figure 5. Landslide susceptibility maps: (a) LSM by IV, (b) LSM by RF, (c) LSM by SVM, and (d) LSM by CNN.
Remotesensing 15 05427 g005
Figure 6. ROC curve and AUC values of the four susceptibility models.
Figure 6. ROC curve and AUC values of the four susceptibility models.
Remotesensing 15 05427 g006
Figure 7. Comparison analysis of the applied models based on confusion matrix-derived parameters and training speed.
Figure 7. Comparison analysis of the applied models based on confusion matrix-derived parameters and training speed.
Remotesensing 15 05427 g007
Figure 8. Landslide raster ratio analysis and comparison based on natural logarithmically scaled FR.
Figure 8. Landslide raster ratio analysis and comparison based on natural logarithmically scaled FR.
Remotesensing 15 05427 g008
Figure 9. Baseline distribution of Sentiel-1A acquisitions.
Figure 9. Baseline distribution of Sentiel-1A acquisitions.
Remotesensing 15 05427 g009
Figure 10. (a) Mean vertical deformation rate PS-points map; (b) Ground deformation map of the study area after Kriging.
Figure 10. (a) Mean vertical deformation rate PS-points map; (b) Ground deformation map of the study area after Kriging.
Remotesensing 15 05427 g010
Figure 11. Landslide dynamic susceptibility maps generated by coupled models: (a) LDSM by IJSVM, (b) LDSM by IWSVM, (c) LDSM by IJRF, and (d) LDSM by IWRF.
Figure 11. Landslide dynamic susceptibility maps generated by coupled models: (a) LDSM by IJSVM, (b) LDSM by IWSVM, (c) LDSM by IJRF, and (d) LDSM by IWRF.
Remotesensing 15 05427 g011
Figure 12. ROC curves and AUC values of RF, SVM, and landslide dynamic susceptibility models.
Figure 12. ROC curves and AUC values of RF, SVM, and landslide dynamic susceptibility models.
Remotesensing 15 05427 g012
Figure 13. Comparison analysis of RF, SVM, and dynamic models based on confusion matrix-derived parameters.
Figure 13. Comparison analysis of RF, SVM, and dynamic models based on confusion matrix-derived parameters.
Remotesensing 15 05427 g013
Figure 14. Landslide raster ratio analysis and comparison of dynamic model based on natural logarithmically scaled FR.
Figure 14. Landslide raster ratio analysis and comparison of dynamic model based on natural logarithmically scaled FR.
Remotesensing 15 05427 g014
Figure 15. The importance of factors in the RF and IJRF.
Figure 15. The importance of factors in the RF and IJRF.
Remotesensing 15 05427 g015
Figure 16. Detailed comparative analysis of typical landslide areas (a) based on the LSM, GDM, and LDSM of the Sifangbei Landslide area (b1b3) and the Cizhuyuan Landslide area (c1c3).
Figure 16. Detailed comparative analysis of typical landslide areas (a) based on the LSM, GDM, and LDSM of the Sifangbei Landslide area (b1b3) and the Cizhuyuan Landslide area (c1c3).
Remotesensing 15 05427 g016
Table 1. Detailed data acquisition dates of Sentinel-1A.
Table 1. Detailed data acquisition dates of Sentinel-1A.
No.Date (yyyy-mm-dd)No.Date (yyyy-mm-dd)
12019-07-17132020-07-11
22019-08-10142020-08-16
32019-09-15152020-09-21
42019-10-09162020-10-15
52019-11-14172020-11-08
62019-12-08182020-12-02
72020-01-13192021-01-07
82020-02-18202021-02-12
92020-03-13212021-03-08
102020-04-18222021-04-13
112020-05-12232021-05-07
122020-06-17242021-06-12
Table 2. Results of multicollinearity assessment (VIF and tolerance value).
Table 2. Results of multicollinearity assessment (VIF and tolerance value).
Impact FactorsVIFTOL
Elevation3.0110.332
Slope1.7550.570
Slope aspect1.0840.923
Plan curvature1.2100.827
Profile curvature1.1130.898
Lithologic1.1860.843
Distance to rivers2.0640.484
Distance to roads1.0800.926
TWI2.1260.470
NDVI2.4430.409
Land use type1.8290.547
Table 3. Detailed parameters of each class in the IV model.
Table 3. Detailed parameters of each class in the IV model.
Assessment FactorsValueLandslide GridsTotal Study AreaInformation Value
CountPercentage (%)CountPercentage (%)
Elevation
(m)
149~250541354.5480,74736.870.708
250~450371037.38129,43643.06−0.142
450~6508028.0883,53527.79−1.235
650~82500.0068462.28−5.421
Slope (°)0~10407741.08129,85743.20−0.505
10~20350835.3592,70530.840.136
20~30198820.0358,59219.490.027
30~403243.2616,7425.57−0.534
>40280.2826680.89−1.146
Slope aspectFlat20.0220,1056.69−5.805
North204020.5552,72917.541.588
Northeast166416.7734,76311.750.371
East101710.2529,6379.860.039
Southeast99710.0528,4779.470.059
South130413.1432,85910.930.184
Southwest6706.7533,01610.98−0.487
West99510.0332,16310.70−0.065
Northwest123612.4536,81512.250.017
Plan curvature−1<458146.16129,10242.950.072
−1~15395.4335,89511.94−0.788
>1480548.41135,56745.100.071
Profile curvature−1<450645.40129,11242.960.055
−1~13553.5830,52510.16−1.044
>1506451.02140,92746.890.085
LithologicJ3s5145.1866,64022.17−1.45
J2s941194.82227,52275.70.225
J2xs00.0011740.39−3.658
J3p00.0050831.69−5.123
J2x00.001450.05−1.566
Distance to river (m)0100.10 30,76210.23−4.62
0~100452745.61 37,44112.461.298
100~200287028.9242,63614.190.712
200~300131613.2671,58223.82−0.586
300~4009159.2257,63819.18−0.732
>4002872.8960,50520.13−1.940
Distance to road (m)0~20103610.4419,9876.650.451
20~50151515.2627,3849.110.516
50~100174017.5335,84111.920.385
100~200214621.6246,37815.430.337
>200348835.14170,97456.88−0.482
TWI0~58918.9837,10612.35−0.319
5~7345834.84118,05139.28−0.112
7~10459546.396,03531.950.371
>109819.8849,37216.43−0.508
NDVI0~0.48108.1646,90415.61−0.648
0.4~0.6140814.1932,15510.700.282
0.6~0.75169017.0335,95311.960.353
0.75~0.9376937.97100,45133.420.128
>0.9224822.6585,10128.31−0.223
Land use typeAgricultural land637764.25171,66057.110.118
Forest1981.9916,9525.64−1.039
Shrubland770.7847831.59−0.718
River1711.72 25,5218.49−1.59
Artificial surface310231.25 81,64827.160.140
Table 4. Main set parameters for PS processing.
Table 4. Main set parameters for PS processing.
ParameterValue
Orbit configurationAscending
Size of scenes40 × 85 km
Number of scenes24
Look azimuth-angle80.45°
Max. temporal baseline396 days
Max. normal baseline119.23 m
Coherence thresholds0.35
Subarea for single reference point25 km2
Overlap for subarea30%
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Miao, F.; Ruan, Q.; Wu, Y.; Qian, Z.; Kong, Z.; Qin, Z. Landslide Dynamic Susceptibility Mapping Base on Machine Learning and the PS-InSAR Coupling Model. Remote Sens. 2023, 15, 5427. https://doi.org/10.3390/rs15225427

AMA Style

Miao F, Ruan Q, Wu Y, Qian Z, Kong Z, Qin Z. Landslide Dynamic Susceptibility Mapping Base on Machine Learning and the PS-InSAR Coupling Model. Remote Sensing. 2023; 15(22):5427. https://doi.org/10.3390/rs15225427

Chicago/Turabian Style

Miao, Fasheng, Qiuyu Ruan, Yiping Wu, Zhao Qian, Zimo Kong, and Zhangkui Qin. 2023. "Landslide Dynamic Susceptibility Mapping Base on Machine Learning and the PS-InSAR Coupling Model" Remote Sensing 15, no. 22: 5427. https://doi.org/10.3390/rs15225427

APA Style

Miao, F., Ruan, Q., Wu, Y., Qian, Z., Kong, Z., & Qin, Z. (2023). Landslide Dynamic Susceptibility Mapping Base on Machine Learning and the PS-InSAR Coupling Model. Remote Sensing, 15(22), 5427. https://doi.org/10.3390/rs15225427

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop