1. Introduction
A landslide is a geological phenomenon that occurs in hilly terrains, resulting in significant ramifications on human lives and financial resources. It possesses the potential to inflict substantial damage on both the ecosystem and society within a specific geographical area [
1,
2,
3]. As the world’s population grows, due to the increase in human engineering activities, and the complex and changeable world climate, the impact of landslides on human beings is progressively growing, and its frequency of attention has also increased year by year. At present, the landslide has become the second-largest geological disaster in the world after the earthquake disaster, which seriously endangers human engineering activities and the safety of human life and property [
4,
5,
6,
7,
8,
9]. Landslide hazard assessment is very important for disaster prevention and mitigation. The Tuoding section of the study area is located in the upper reaches of the Jinsha River, northwest of Yunnan and the southeast edge of the Qinghai–Tibet Plateau. Under the special tectonic background in the region, there is a complex high ground stress field, strong tectonic activity, frequent occurrence of geological disasters such as landslides, and complex and changeable climatic conditions. These geological hazards frequently lead to river blockages, triggering landslides and earthquakes, as well as the formation of dammed lakes that can result in devastating floods. These events pose a significant threat to the safety of residents living downstream of the affected river [
10,
11]. The Baige landslide and river blocking occurred twice in the upper reaches of Jinsha River, resulting in the destruction of a large number of houses and a direct economic loss of about 963.5 million US dollars [
11]. In addition, there exist numerous historical landslide formations on both banks of the Jinsha River [
12,
13], threatening the lives of downstream residents and buildings. Therefore, it is necessary to carry out the study on the susceptibility assessment of potential landslides along the Jinsha River.
At present, domestic and foreign scholars mainly use artificial neural network [
14], logistic regression [
15], frequency ratio [
16], weight of evidence [
17], and other methods to evaluate landslide susceptibility. The United States and some European countries were the first to study landslide susceptibility. The United States Geological Survey (USGS) commenced landslide susceptibility assessment during the 1970s, and developed a variety of assessment models, such as the weighted stability index model (WIS), the physical stability index model (PSI), etc. These models take into account multiple factors such as topography, geology, and climate, and have high accuracy. In addition, European countries have also conducted in-depth research on landslide hazard zoning, and countries such as Germany and France have also proposed corresponding assessment methods, such as the landform evolution influence factor method (FOSM). The earliest systematic exploration results can be traced back to 1964. Dobrovony analyzed and interpreted the geological map, topographic map, and geomorphic features of the Anchorage area, and classified the sensitivity of the local landslide [
18]. Regmi et al. evaluated the susceptibility to landslides through GIS techniques using Bayes’ theorem based on weights of evidence [
19]. During the 1976 International Engineering Geological Congress, the attribution of natural disasters to geological phenomena was established, and the term “geological disaster” was formally defined for the first time. The analysis and research work of landslide disasters has also become more clear. In the 1980s, the geographic information system officially entered the stage of landslide hazard classification. The powerful data analysis function of the system made the study of landslide hazards gradually transform from qualitative research to quantitative research, making the evaluation methods more diverse and scientific. Brand E W used quantitative methods to try to calculate the susceptibility of disasters, expounded the development of landslide geological disasters in Hong Kong, and put forward suggestions for susceptibility management and prevention [
20]. Wei chen combined GIS technology to use a new hybrid computational intelligence model on different mapping units to evaluate landslide susceptibility [
21]. Bourenane validated and compared the landslide susceptibility maps (LSMs) produced by applying four geographic information system (GIS)-based statistical approaches including frequency ratio (FR), statistical index (SI), weights of evidence (WoE), and logistic regression (LR) for the urban area of Azazga [
22]. Based on the distribution characteristics and scale of landslides in the field survey, the study area was divided into susceptibility zones by using the superposition analysis function of GIS. This study made the landslide susceptibility assessment more systematic for the first time. In the 1990s, geographic information technology developed rapidly. GIS technology, with its powerful spatial analysis ability, provides a powerful means for the extraction of landslide disaster characteristics and the calculation of landslide susceptibility. It has been combined with various types of machine learning and artificial intelligence to make landslide susceptibility zoning more and more scientific [
23,
24,
25]. Basofi used the analytic hierarchy process combined with the natural discontinuity method in the process of landslide susceptibility assessment in Indonesia, and achieved good results [
26].
In recent years, the focus of research on landslide susceptibility assessment has shifted from exploring new models to utilizing a combination of GIS and various mathematical statistical models, along with advanced technologies such as InSAR, satellite and aerial photography, and machine learning. Skillful integration of these technologies enables more scientific and rational discovery and interpretation of landslides (including potential ones), enhances the reliability of data, and improves the accuracy of susceptibility assessment results [
27]. For example, Liu Jing successfully employed GIS in combination with UAV technology to assess landslide susceptibility in Zhouqu County, Gansu Province, yielding positive outcomes [
28]. Novellino A utilized InSAR technology and machine-learning methods in the process of landslide susceptibility assessment within the study area, showcasing their effectiveness [
29]. In the regional landslide susceptibility assessment conducted by Fang Ranke, machine-learning methods are comprehensively summarized, and novel insights for susceptibility assessment are proposed [
30].
The advancements in InSAR technology have opened up possibilities for landslide susceptibility assessment [
31]. InSAR technology has been utilized for surface deformation monitoring since 1969 and has proven effective in identifying various geological disasters such as earthquakes, surface subsidence, and landslides. This technology plays a crucial role in our efforts toward disaster prevention and mitigation [
32,
33,
34,
35,
36]. Currently, domestic research primarily focuses on identifying indicators for landslide susceptibility assessment and improving assessment methods. For instance, Wu Shuren et al. successfully evaluated landslide susceptibility in Fengdu County, Chongqing, using a geological disaster information system and information quantity method [
37]. Fan Linfeng and other scholars employed a weighted information model to assess landslide susceptibility in Enshi City, Hubei Province [
38]. Li Wenyan and colleagues compared the frequency ratio and information quantity models in their study of landslide susceptibility in certain areas of Gaolan County, Gansu Province [
39]. Foreign research generally possesses a mature theoretical framework and diverse evaluation models, but the applicability of these models is limited by geological and climatic conditions. Domestic research, on the other hand, places more emphasis on practical application and has achieved significant results through the continuous improvement and optimization of evaluation methods. When combined with optical imagery and on-site field investigations, this approach can effectively verify and validate the accuracy of the findings [
40,
41,
42,
43].
The objective of this study is to utilize ALOS-2 and Sentinel-1A data to detect potential landslides in the Tuoding section of the Jinsha River and assess the suitability of the data. The study is conducted in three stages: (a) SBAS-InSAR technology is employed to analyze the long-term phase changes, while multi-period Google Earth images are utilized to identify potential landslides. This approach helps reduce result uncertainties. (b) Field investigations are carried out to verify the identification results and analyze the genetic mechanisms of typical landslides. (c) Based on the findings from ALOS-2 and Sentinel-1A, the deformation characteristics of specific landslides such as Shentingla (STL), Kongzhigong (KZG), Dingzhui (DZ), and Duila (DL) are analyzed. Finally, the reliability and applicability of the ALOS-2 and Sentinel-1A data are validated.
2. Study Area
The study region is situated within Tuoding Lisu Township, Deqin County, Diqing Tibetan Autonomous Prefecture, northern Yunnan Province (
Figure 1a). The area is 28 km east of Shangri-La City, 13 km west of Xiaruo Lisu Township, 19 km south of Tacheng Town, and 30 km north of Nixi Township. Jinsha River passes through the study area from north to south. The left bank is Wujing Township, and the right bank is Tuoding Lisu Township. The study section flows through Shengli, Baishugong, Chigutang, Guilong, Badong, and other places. The specific range is (99°18′41″–99°36′14″ E), (27°18′23″–27°58 ′30″ N) rectangular area, with an area of 1230.36 square kilometers. The residents in the area are mainly Lisu people, and there are Tibetan, Han, and other population distributions. Tuoding Township is connected with Shangri-La City and surrounding major towns through national highways G214 and G215.
The study site is situated within the tectonically active region where the Indian plate subducts beneath the Eurasian continental plate. The tectonic compression is strong, and the vegetation in the area is developed, which is a heavy forest area. The region’s robust topographical uplift has resulted in the formation of a V-shaped, deeply incised valley, characterized by prominent denudation processes and glacial erosion, thus rendering it an alpine canyon landform [
42]. The river width of the study area is usually between 60 m and 100 m. The study area is an incised river section, with high mountains and deep slopes. The villages in the area are distributed in strips on both sides of the Jinsha River. The elevation of the study area is 1862~4502 m, and the height difference is 2640 m.
The geological formations exposed within the investigated region comprise Cenozoic Quaternary (Q) Holocene and Pleistocene strata, Tertiary (E) Eocene Lumeiyi Formation strata, Mesozoic Triassic (T) Upper Zhongwo Formation, Middle Beiya Formation and Hewanjie Formation, Lower Lamei Formation, Permian (T) Upper Basalt Formation, Paleozoic Carboniferous (C) Lower Formation, Devonian (D) Lower Xiaoyangpo Formation, Middle Lingpaishan Formation and Guangtoupo Formation, Lower Ranjiawan Formation and Hailuo Formation, and Cambrian (∈) Lower Tacheng Formation and Yangpo Formation. The stratigraphic age spans from the Cenozoic to the Paleozoic, and there are many types of strata. The missing strata are Cretaceous (K), Jurassic (J), and Ordovician (O). The distribution patterns of the geological strata within the study area exhibit discernible regularities, and the strata gradually become older from W to E. Different series are usually in fault contact. The fault is a sub-fault of the eastern branch fault zone of Jinsha River in F1 and the Zhongdian–Longpan–Qiaohou fault zone in F2, while the strata in the same series are generally in integrated contact (
Figure 1b).
The seasonal performance of the study area is that the rainy season of the basin is from May to October every year. Influenced by both the southwest and southeast monsoons, the study area experiences an abundance of atmospheric moisture, resulting in concentrated precipitation. The mean annual precipitation measures 954.0 mm, while the average annual evaporation stands at 2179 mm. Furthermore, the region exhibits an annual mean runoff of 1360 m
3/s, accompanied by an average annual temperature of 12.6 °C. According to the “China ground motion parameter zoning map”, the peak acceleration of the ground motion of the class II site in Tuoding Township, Deqin County, Diqing City, Yunnan Province is 0.20 g. The ground motion response spectrum within the study area exhibits a distinctive peak at a period of 0.40 s, corresponding to a seismic intensity of VIII on the scale. The predisposition to slope instability on both sides of the river, as well as the recurring incidents of landslides, can be attributed to a confluence of factors, including robust tectonic activity, ongoing fluvial incision erosion, prolonged weathering processes, a complex high ground stress field, and the effects of free surface unloading [
43].
4. Discussion
InSAR technology is a reliable method for identifying potential landslides. Previous studies have utilized Sentinel-1 and ALOS-2 satellite data for surface deformation monitoring, which we also used in this study to identify potential landslides in the study area. However, while interpretation methods such as visible light remote sensing technology and InSAR technology are generally scientific, their measurement accuracy can be impacted by factors such as cloud cover and surface reflectance. This can result in individual landslide data being accidental and unreliable. Therefore, we should not rely solely on remote sensing technology for landslide data survey and instead aim to collect data through multiple channels. Furthermore, we need to enhance the processing capacity of InSAR technology to reduce interference and improve accuracy, ensuring that our survey work is more precise. By doing so, we can obtain more reliable data to inform our analysis and decision-making related to landslides.
The examination of factors contributing to landslide susceptibility zoning reveals that regions characterized by very low susceptibility and low susceptibility are primarily clustered in the western sector, and in the southeast part of the region, these areas are characterized by high altitudes and dense vegetation, which in turn contribute to enhanced slope protection through the presence of vegetation on the slopes, and because of its high altitude, rivers, and roads. It is difficult to affect these areas, making these places safer; the high-susceptibility and higher-susceptibility areas are mainly concentrated on both sides of the Jinsha River and near the highways along the river. This phenomenon can likely be attributed to the persistent erosive forces exerted by the river along the base of the slopes on both flanks, coupled with the adverse impacts stemming from construction activities along the adjacent highways, which have resulted in damage and disruption to the slopes. The slope is more fragile, and because the elevation of the region is lower, the lower the vegetation in the low-altitude area closer to the foot of the slope on both sides is, the rock mass is bare and vulnerable to erosion, so the landslide susceptibility is greater. As a result, it is crucial to prioritize the protection of landslides on both sides of the river and implement appropriate preventive measures to reduce the dangerous consequences of landslides. By doing so, we can minimize the impact of landslides and ensure the safety of the local population.
We organized and calculated the correlation between nine evaluation factors: elevation, slope, aspect, distance from river, vegetation coverage, distance from road, distance from fault, lithology, and precipitation. This allowed us to have a more intuitive understanding of the relationship between each factor’s range and landslide occurrence, facilitating a deeper understanding of local landslides. Based on the results of the landslide hazard analysis and single-factor data in the study area, it is known that landslides in the study area are mainly concentrated in areas with elevations below 2600 m, slopes of 18–30 degrees, aspects of 0–90 degrees, distances from rivers of 500–800 m, distances from roads of 200–500 m, vegetation coverage below 25%, distances from faults less than 500 m, and rainfall exceeding 980 mm. Areas with a higher overall landslide risk in the study area are concentrated along the rivers and roads running north–south, as well as near geological faults. There is still room for development in the selection of factors. Some factor data only reflect the moment of collection and do not consider the impact of time. In future research, in addition to ensuring the screening capability of factors, the timeliness of factors should also be considered.
Based on the results of the landslide sensitivity analysis, the study area should pay special attention to the banks of the Jinsha River, Jirenshui, Shenta, Lulu Ge, Guiba, and the easternmost part of the study area near the Shangri-La region. Monitoring measures should be strengthened and early preventive measures taken to reduce or avoid significant losses caused by landslides to the local area. Local governments should develop prevention plans, enhance mountain monitoring during heavy rainfall periods, and establish emergency overtime systems during flood seasons. At the same time, it is necessary to regulate construction work and strengthen supervision of the construction process to minimize disturbance to local slopes caused by construction activities. This can enhance the understanding of the hazards posed by rainfall and human engineering activities as landslide triggering factors, effectively protect the interests of the local people, and minimize harm.
5. Conclusions
This paper employs InSAR technology to acquire surface deformation monitoring data for the Topping River section in the study area. These data are utilized to identify potential landslide boundaries, which are further validated through field investigations. Vulnerability zoning of the nine influencing factors is conducted using the frequency ratio model and information quantity model. The accuracy of the results is verified using the ROC curve and density method, with the information quantity model outperforming the frequency ratio model. The information quantity model’s calculation results, combined with the analytic hierarchy process and superposition point density map, are utilized to perform landslide susceptibility zoning in the study area. The results reveal five susceptibility zones, ranging from low to high, with respective areas of 197.48, 455.72, 408.21, 152.66, and 16.22 km2. The corresponding proportion of landslides in each zone is 0.17%, 1.60%, 3.88%, 8.41%, and 16.65%. As expected, the proportion of landslides significantly increases with higher susceptibility levels, confirming the accuracy of the susceptibility results. Finally, this study analyzes the single-factor classification map of the study area based on the susceptibility zoning results. This analysis provides a reasonable explanation for the varying distribution of local landslide susceptibility. The insights gained from this study contribute to the understanding of landslide susceptibility zoning and can inform the development of effective preventive measures.