Landslide Detection in the Linzhi–Ya’an Section along the Sichuan–Tibet Railway Based on InSAR and Hot Spot Analysis Methods

: Construction of the 998.64-km Linzhi–Ya’an section of the Sichuan–Tibet Railway has been inﬂuenced by landslide disasters, threatening the safety of Sichuan–Tibet railway projects. Landslide identiﬁcation and deformation analysis in this area are urgently needed. In this context, it was the ﬁrst time that 164 advanced land-observing satellite-2 (ALOS-2) phased array type L-band synthetic aperture radar-2 (PALSAR-2) images were collected to detect landslide disasters along the entire Linzhi–Ya’an section. Interferogram stacking and small baseline interferometry methods were used to derive the deformation rate and time-series deformation from 2014–2020. After that, the hot spot analysis method was introduced to conduct spatial clustering analysis of the annual deformation rate, and the effective deformation area was quickly extracted. Finally, 517 landslide disasters along the Linzhi–Ya’an route were detected by integrating observed deformation, Google Earth optical images, and external geological data. The main factors controlling the spatial landslide distribution were analyzed. In the vertical direction, the spatial landslide distribution was mainly concentrated in the elevation range of 3000–5000 m, the slope range of 10–40 ◦ , and the aspect of northeast and east. In the horizontal direction, landslides were concentrated near rivers, and were also closely related to earthquake-prone areas, fault zones, and high-precipitation areas. In short, rainfall, freeze–thaw weathering, seismic activity, and fault zones are the main factors inducing landslides along this route. This research provides scientiﬁc support for the construction and operation of the Linzhi–Ya’an section of the Sichuan–Tibet Railway. and analysis.


Introduction
The Sichuan-Tibet Railway, China, travels from eastern Chengdu, Sichuan to Lhasa, Xizang, passing through the Dadu River, Yalong River, Jinsha River, Lancang River, and Nujiang River, and, finally, travelling along the Yarlung Zangbo River to Lhasa, with a total length of 1533.64 km. It is a bridge connecting Tibet and the central region of China, which demonstrates immense strategic significance for the development of the Tibetan economy and maintenance of border stability [1]. Railway construction of the Chengdu-Ya'an section and Linzhi-Lhasa section have been completed, while the Linzhi-Ya'an section is under the project implementation stage. The planned Linzhi-Ya'an section is approximately 998.64 km long and passes through 21 tunnels, accounting for approximately 80.4% of the total length. There are many active fault zones along the route with extremely strong geological tectonic activities [2]. Moreover, the stratum lithology is complex, and there exists considerable variation in the elevation of the terrain, which mostly belongs the InSAR method. First, we explored and used targeted methods to correct errors in InSAR data processing, such as atmospheric error, geometric distortion, unwrapping error, and digital elevation model (DEM) error, and the accuracy of the InSAR method for surface deformation detection was effectively improved. Then, the hot spot analysis method was used to realize automatic, rapid, and accurate identification of deformation regions with spatial clustering properties of the InSAR deformation rate. Landslide disasters were screened by superimposing and analyzing optical remote sensing images and terrain data. Then, a disaster inventory map was generated. Finally, the spatial distribution and main controlling factors of landslides were analyzed based on the catalogue results, which provided a reference for the study of landslide disasters along the Sichuan-Tibet Railway, and the information could be applied to engineering construction.

Study Area
The Linzhi-Ya'an section is currently under construction as part of the Sichuan-Tibet Railway, and spans from the west edge of the Sichuan Basin to the eastern part of the Qinghai-Tibet Plateau, crossing the Himalayas, Nianqing Tanggula Mountains, and Hengduan Mountains, among other regions, as shown in Figure 1. It demonstrates the following features: (1) the terrain is steep and fluctuates considerably, with an elevation range of 627-5013 m, and is high in the west and low in the east, and alpine valleys are widespread [1]. (2) It passes through Dadu River, Yalong River, Jinsha River, Lancang River, and Nujiang River. Therefore, a dense river distribution leads to a remarkable cutting effect [15]. (3) The climate along the route transits from a temperate monsoon climate to a plateau alpine climate, which can be divided into six climate types: subtropical humid climate zone, western Sichuan plateau temperate humid climate zone, eastern Tibet plateau temperate semi-humid climate zone, southeast Tibet subtropical mountainous humid climate zone, southern Qiangtang plateau subtropical semi-arid climate zone, and southern Tibet plateau temperate semi-arid climate zone [16]. The period of rainfall and high temperature is mostly concentrated in summer, which is also a period of high landslide occurrence [17]. (4) The temperature difference between day and night in this region is substantial, and freeze-thaw weathering is extremely strong, which further reduces the stability of the rock mass [18,19]. (5) There are many active faults along the Sichuan-Tibet Railway and the geological tectonic activity is intense. It passes through 24 major faults, such as the Xianshui River, Yushu-Ganzi, Litang-Dewu, Lancang River, Nu River, and Jinsha River faults. Additionally, since the Quaternary, the Indian plate has been subducted and squeezed to the Eurasian plate and the Qinghai-Tibet Plateau has risen rapidly at a rate of 9.5 mm/year; therefore, seismic activity is frequent [20], which constitutes a background for surface deformation in this region. (6) There are substantial stratigraphic differences and weak lithologies along this route; the stratigraphic direction is consistent with the trend of the fault zone. Based on the interaction of the above-mentioned natural characteristics, fractures in the rock mass are widespread, which is conducive to rainwater infiltration [21]. These factors not only lead to frequent landslide disasters, but also induce large-scale, high-level, and long-distance landslides. For example, the Yigong high-level and huge landslide occurred in the year 2000. The sliding time of this landslide lasted for 10 min, and the sliding distance reached 8 km. The Yigong Zangbo River was markedly blocked and the highway was interrupted for several months. It also resulted in large-scale casualties and property losses upstream and downstream [22,23]. In 2018, a large-scale high-level rock landslide occurred in Baige Village, which caused a break in the Jinsha River for 2 days [24]. The economic loss caused by this landslide amounted to CNY 15 billion. Unmanned aerial vehicle (UAV) images of these two landslides are shown in Figure 2. These landslides demonstrated the characteristics of high concealment, rapid sliding speed, long sliding distance, and huge destructive force, which are deemed major challenges for the planning and construction of the Sichuan-Tibet Railway.

Data
To better overcome the decoherence caused by vegetation and ice-snow cover along the railway during InSAR processing, we collected 164 L-band ALOS/PALSAR-2 ascending images, which demonstrate good penetration ability for surface covering, to detect the hidden landslides along the Sichuan-Tibet Railway. Considering the low temporal resolution of ALOS/PALSAR-2 satellite (46 days), the time period of data from 2014 (satellite launching time) to 2020 was selected in this study. The cover range and corresponding parameters of the 16 orbits in the study area are shown in Figure 1 and Table 1, respectively. The visibility, layover, and shadow area of the images were calculated and masked according to the flight parameters of the satellite and external SRTM DEM with a 30-m resolution. Although the research area topography was mostly alpine canyon, according to the overall visibility map of the research area shown in Figure 3, the visibility area accounted for 91.5%, which proved that the selected data could be considered for disaster detection in most study areas. Additionally, SRTM DEM was also used for InSAR data processing. sliding speed, long sliding distance, and huge destructive force, which are deemed major challenges for the planning and construction of the Sichuan-Tibet Railway.

Data
To better overcome the decoherence caused by vegetation and ice-snow cover along the railway during InSAR processing, we collected 164 L-band ALOS/PALSAR-2 ascending images, which demonstrate good penetration ability for surface covering, to detect the hidden landslides along the Sichuan-Tibet Railway. Considering the low temporal resolution of ALOS/PALSAR-2 satellite (46 days), the time period of data from 2014 (satellite launching time) to 2020 was selected in this study. The cover range and corresponding parameters of the 16 orbits in the study area are shown in Figure 1 and Table 1, respectively. The visibility, layover, and shadow area of the images were calculated and masked according to the flight parameters of the satellite and external SRTM DEM with a 30-m resolution. Although the research area topography was mostly alpine canyon, according to the overall visibility map of the research area shown in Figure 3, the visibility area accounted for 91.5%, which proved that the selected data could be considered for disaster detection in most study areas. Additionally, SRTM DEM was also used for InSAR data processing.

Data
To better overcome the decoherence caused by vegetation and ice-snow cover along the railway during InSAR processing, we collected 164 L-band ALOS/PALSAR-2 ascending images, which demonstrate good penetration ability for surface covering, to detect the hidden landslides along the Sichuan-Tibet Railway. Considering the low temporal resolution of ALOS/PALSAR-2 satellite (46 days), the time period of data from 2014 (satellite launching time) to 2020 was selected in this study. The cover range and corresponding parameters of the 16 orbits in the study area are shown in Figure 1 and Table 1, respectively. The visibility, layover, and shadow area of the images were calculated and masked according to the flight parameters of the satellite and external SRTM DEM with a 30-m resolution. Although the research area topography was mostly alpine canyon, according to the overall visibility map of the research area shown in Figure 3, the visibility area accounted for 91.5%, which proved that the selected data could be considered for disaster detection in most study areas. Additionally, SRTM DEM was also used for InSAR data processing.

Methods
The collected 16-orbit data were preprocessed to obtain single-look complex (SLC) images, and then the combinations of SLC pairs were interfered and unwrapped. Because the terrain along the Sichuan-Tibet Railway fluctuates greatly, the vegetation and icesnow cover in some areas is serious, and many errors are inevitable during data processing (such as decoherence, atmospheric error and geometric distortion). In this study, errors in the interference and unwrapping process were explored and corrected according to the characteristics of the study area. After selecting high-quality unwrapping maps, the annual deformation rate was calculated based on the stacking-InSAR method. Based on the results, the effective deformation area was directly extracted by the hot spot analysis method. Then, the identification of potential landslides was completed based on effective deformation combined with optical images, geological factors, and time-series deformation methods. Finally, the main controlling factors of landslides were analyzed using a spatial analysis method. A flowchart of the specific method is provided in Figure 4. Detailed methods of stacking, SBAS-InSAR, and hot spot analysis are introduced below.

Methods
The collected 16-orbit data were preprocessed to obtain single-look complex (SLC) images, and then the combinations of SLC pairs were interfered and unwrapped. Because the terrain along the Sichuan-Tibet Railway fluctuates greatly, the vegetation and ice-snow cover in some areas is serious, and many errors are inevitable during data processing (such as decoherence, atmospheric error and geometric distortion). In this study, errors in the interference and unwrapping process were explored and corrected according to the characteristics of the study area. After selecting high-quality unwrapping maps, the annual deformation rate was calculated based on the stacking-InSAR method. Based on the results, the effective deformation area was directly extracted by the hot spot analysis method. Then, the identification of potential landslides was completed based on effective deformation combined with optical images, geological factors, and time-series deformation methods. Finally, the main controlling factors of landslides were analyzed using a spatial analysis method. A flowchart of the specific method is provided in Figure 4. Detailed methods of stacking, SBAS-InSAR, and hot spot analysis are introduced below. Remote Sens. 2021, 13, x 6 of 21

Stacking-InSAR
The stacking-InSAR method can be used in conjunction with the differential interferometric phase after unwrapping to obtain the annual deformation rate. The atmospheric phase error in the superimposed phase can be weakened, which can effectively be used to suppress atmospheric delay and to obtain high-precision deformation monitoring results.
In the present study, interferograms with residual error and low coherence were eliminated, and a certain number of high-quality unwrapping maps were used to obtain the phase deformation rate, as follows [25]: where represents the weight of the unwrapped phase; the time intervals are considered as the weight, i.e., = ∆ ; ℎ represents the th interference unwrapping phase; and represents the number of interferograms.

SBAS-InSAR
SBAS-InSAR is a short baseline differential interferometry method proposed by Berardino et al. in 2002, which can be used to obtain a continuous deformation time series [26]. This method exhibits a low requirement for SAR data. It also can weaken the influence of decoherence and topographic changes by controlling the time and spatial baselines [27].
In the present study, all SAR images were divided into different short baseline subsets according to a certain time and spatial threshold, and then each combinatorial pair was subjected to interference. After interference processing, the minimum cost flow

Stacking-InSAR
The stacking-InSAR method can be used in conjunction with the differential interferometric phase after unwrapping to obtain the annual deformation rate. The atmospheric phase error in the superimposed phase can be weakened, which can effectively be used to suppress atmospheric delay and to obtain high-precision deformation monitoring results.
In the present study, interferograms with residual error and low coherence were eliminated, and a certain number of high-quality unwrapping maps were used to obtain the phase deformation rate, as follows [25]: where w i represents the weight of the unwrapped phase; the time intervals are considered as the weight, i.e., w i = ∆T i ; ph i represents the ith interference unwrapping phase; and N represents the number of interferograms.

SBAS-InSAR
SBAS-InSAR is a short baseline differential interferometry method proposed by Berardino et al. in 2002, which can be used to obtain a continuous deformation time series [26]. This method exhibits a low requirement for SAR data. It also can weaken the influence of decoherence and topographic changes by controlling the time and spatial baselines [27].
In the present study, all SAR images were divided into different short baseline subsets according to a certain time and spatial threshold, and then each combinatorial pair was subjected to interference. After interference processing, the minimum cost flow (MCF) unwrapping method was used for phase unwrapping [28,29], and the observation equation was established, combined with high coherence points. Finally, the least square (LS) or singular value decomposition (SVD) method was used to obtain the minimum norm solution of the velocity vector, and the cumulative deformation of each time period could be obtained by integration in the time domain [26]. To overcome the spatial-temporal decorrelation and consider the low temporal resolution of ALOS-2 data, the temporal and spatial baseline thresholds were set as 500 d and ±700 m, respectively, to obtain more interferograms. Owing to the influence of satellite geometric parameters and terrain factors, the shortened, layover, and shadow errors result in the distortion of SAR images during the detection process [30], a phenomenon which affects the registration of SAR images and leads to the generation of certain noise areas. In the present study, the geometric parameters of satellites and DEM data were used to calculate and effectively remove the shadows and layover areas in the GAMMA software. The error areas could be eliminated, and the accuracy of the interferograms was improved. We also introduced the SRTM DEM data to remove the terrain phase. The adaptive filtering method was used to remove the noise existing in the interferograms [31]. Notably, during filtering, a small filtering window (32 or 16) was selected to avoid excessive information loss. The atmospheric error was corrected by adopting the empirical linear model constructed between the InSAR phase and elevation [32]. In short, the errors were corrected in different situations to obtain the optimal deformation results during data processing.

Automatic Extraction of Effective Deformation Data via Hot Spot Analysis
As the spatial range of the deformation rate is markedly considerable, and as few discrete error deformation points remain, direct landslide identification is deemed timeconsuming and laborious. Therefore, this study introduced a hot spot analysis method for large-scale landslide detection. Through the spatial statistical analysis of the obtained large-scale deformation values, the high and low points with evident clustering properties in space were obtained.
The hot spot statistical method used in the study was the Getis−Ord G * i method. G i (d) statistics is a method used to monitor whether the feature points possess significant spatial clustering properties by checking each element in the adjacent element environment within a certain distance [33,34]. The calculation formula is as follows: where w ij (d) represents the spatial weight matrix; n represents the total number of elements; and x represents the surface displacement recorded by InSAR. Assuming that G i (d) obeys laws of normal distribution, the Z score of each element point can be calculated using Equation (3); the higher the Z score, the more evident the spatial clustering properties. The following equation is considered: Among them, represents the expectation. In the present study, p < 0.01 (confidence level is 99%) and z > ±|2.58| of the deformation point was a reliable deformation with evident spatial clustering properties, and p was considered the significance level in the hypothesis test.

Deformation Rate Map Acquired via InSAR
Annual line of sight (LOS) deformation rate maps of the Ya'an-Linzhi railway section were generated using the stacking-InSAR method, as shown in Figure 5. The region showing an annual deformation rate value of <±10 mm was considered a stable region (represented by using green in Figure 5). Notably, the negative values (represented by using red) indicated that the surface moved away from the satellite, whereas the positive values (represented by using blue) indicated that the surface moved toward the satellite sensor ( Figure 5) [35]. It could be inferred that the coherence of the Ya'an-Kangding and Luolong-Linzhi sections was relatively low compared to other regions, which might have been caused by the complex environment [36,37]. Meanwhile, according to the annual deformation rate map, we found that large-scale errors (such as atmospheric and unwrapping errors) could be effectively removed. The spatial landslide distribution was uneven; there were more deformation regions in the east and west, whereas the middle part was relatively stable. Additionally, most deformation regions were distributed along expansive rivers. (represented by using green in Figure 5). Notably, the negative values (represented by using red) indicated that the surface moved away from the satellite, whereas the positive values (represented by using blue) indicated that the surface moved toward the satellite sensor ( Figure 5) [35]. It could be inferred that the coherence of the Ya'an-Kangding and Luolong-Linzhi sections was relatively low compared to other regions, which might have been caused by the complex environment [36,37]. Meanwhile, according to the annual deformation rate map, we found that large-scale errors (such as atmospheric and unwrapping errors) could be effectively removed. The spatial landslide distribution was uneven; there were more deformation regions in the east and west, whereas the middle part was relatively stable. Additionally, most deformation regions were distributed along expansive rivers.

Landslide Deformation Rate Map Obtained via Hot Spot Analysis
Based on the annual deformation rate map obtained by InSAR, a hot spot analysis was performed to automatically extract the effective deformation area. In the present study, five areas with a dense distribution of large-scale landslides were selected for illustration in Figure 6 (note: the locations are shown as yellow rectangles in Figure 5). It could be inferred from Figure 6 that the stable area and the discrete noise points removed after conduction of the hot spot analysis were related to the original deformation rate map, and that only data on the effective deformation area remained. These effective deformation points were more prominent [33], which could markedly be considered to reduce the difficulties of landslide interpretation and rapid identification of landslide hazards.
To further illustrate the landslide mapping results, Figure 7 has been depicted, which shows the original results, hot spot analysis results, and optical images of 17 typical largescale landslides in the key area. The specific locations of the landslide points are shown in Figure 6. By comparing the results, we found that the hot spot analysis could be performed to effectively eliminate the error area, and that the extracted landslide boundary was in good agreement with the optical image. It could also be inferred that some detected landslides were not easily distinguished on optical images because of their high concealment characteristics. However, the landslide could be identified by InSAR, which showed that the InSAR method could not only be used to identify landslides with evident deformation in optical images, but also be used to identify deformed landslides with no remarkable characteristics on optical images. Particularly, the hot spot analysis method can be used to determine the real hidden dangers posed by potential landslides and to provide scientific support for railway planning and construction.

Landslide Deformation Rate Map Obtained via Hot Spot Analysis
Based on the annual deformation rate map obtained by InSAR, a hot spot analysis was performed to automatically extract the effective deformation area. In the present study, five areas with a dense distribution of large-scale landslides were selected for illustration in Figure 6 (note: the locations are shown as yellow rectangles in Figure 5). It could be inferred from Figure 6 that the stable area and the discrete noise points removed after conduction of the hot spot analysis were related to the original deformation rate map, and that only data on the effective deformation area remained. These effective deformation points were more prominent [33], which could markedly be considered to reduce the difficulties of landslide interpretation and rapid identification of landslide hazards.
To further illustrate the landslide mapping results, Figure 7 has been depicted, which shows the original results, hot spot analysis results, and optical images of 17 typical large-scale landslides in the key area. The specific locations of the landslide points are shown in Figure 6. By comparing the results, we found that the hot spot analysis could be performed to effectively eliminate the error area, and that the extracted landslide boundary was in good agreement with the optical image. It could also be inferred that some detected landslides were not easily distinguished on optical images because of their high concealment characteristics. However, the landslide could be identified by InSAR, which showed that the InSAR method could not only be used to identify landslides with evident deformation in optical images, but also be used to identify deformed landslides with no remarkable characteristics on optical images. Particularly, the hot spot analysis method can be used to determine the real hidden dangers posed by potential landslides and to provide scientific support for railway planning and construction.

Active Landslide Inventory
In addition to landslides, large-scale disasters, such as ice collapse, glacier lake outburst, and debris flow, also occurred due to the complex geological and geomorphological conditions along the Sichuan-Tibet Railway. Therefore, landslide disasters must be selected by considering reliable surface deformation data that strictly include information based on geological factors, such as DEM and slope aspect. Furthermore, we mainly identified landslides that posed a threat to the construction of tunnels and bridges within 100 km, and did not consider small landslides and disasters that posed no threat to the railway, such as small collapses and landslides along the river banks. Based on the annual deformation rate after conduction of hot spot analysis, 517 landslides with evident deformation and certain threats to railways were interpreted (Figure 8). Many detected landslides have been validated by field investigation and published articles, such as those pertaining to the Yigong landslide, Lancang River, Jinsha River, and Yalong River landslide groups [11,12,[38][39][40][41][42]. By comparing the identification results of hot spot analysis and traditional manual visual interpretation, it was found that traditional visual interpretation resulted in the exclusion of data on 178 missed landslides and 58 misjudged landslides. This finding proves that hot spot analysis can be performed to solve the problem of unintentional exclusion of missed and misjudged landslides during identification to a considerable extent.

Time-Series Monitoring for Typical Landslides
Based on the landslide identification results, certain high-risk landslides can be monitored using time-series deformation. Therefore, we selected two landslide bodies, G and H, as shown in Figure 9, both of which demonstrate an evident continuous deformation trend. As shown in Figure 9b, the maximum cumulative deformation of landslide G reaches −400 mm, and a continuous deformation trend exists at all three points. This landslide is located on the right bank of the Jinsha River, and there are multiple tensile cracks in the back edge of the slope, and findings indicate that the range of cracks gradually expands with the continuous increase in deformation. Therefore, human activities on the slope, river cuts, and rainfall may cause failure of this slope. If sliding occurs, it may block the Jinsha River and endanger the upstream bridge, a phenomenon which may also induce failure of other unstable slopes and may expand the scale of damage. Therefore, continuous monitoring is recommended in the future. Landslide H is an old landslide, and there is less vegetation cover on this slope. Optical imaging results reveal that it is a freeze-thaw landslide. According to the time-series deformation of this landslide, the maximum cumulative deformation of P3 on this landslide can reach −296 mm (Figure 9d), and the three points demonstrate a state of continuous deformation, a phenomenon which will also pose a grave threat to the construction and operation of nearby buildings and structures. It could be inferred from the time-series deformation of these two disaster points that although ALOS-2 data exhibited certain disadvantages in terms of time resolution, the overall cumulative deformation trend could also be studied appropriately; therefore, such landslides can be monitored effectively by considering InSAR time-series deformation. Landslide analysis showed that the overall spatial landslide distribution was uneven, and was relatively dense in Luolong-Basu, Changdu-Tuoba, Jinsha River coast, and Kangding-Ya'an, wherein a total of 358 landslides were identified, accounting for 69.2% of the total number of landslides. It could be inferred from Figures 5 and 8 that the large-scale landslides were mainly concentrated in Lulang-Guxiang, Yuka-Xiangdui, Yongba-Juebo, Jinsha River coast, and Guza-Tianquan County. The large-scale landslides in the Lulang-Guxiang section (Figure 6a) mainly developed along the Palalung Zangbo River. The river flow rate in this area is high and annual runoff can reach 700 m 3 /s [15]. River erosion, significant tectonic activities, and freeze-thaw weathering are the main triggering factors for the frequent occurrence of landslides in this area [4]. Moreover, the annual deformation rate of most landslides is >100 mm/year and the volume is >1 km 2 . For example, the famous Yigong landslide occurred along the Palalung Zangbo River [22,23]. Regarding the coastal area of the Jinsha River (Figure 6b), several large-scale landslides were distributed in this area and the deformation rates of these slopes were all >100 mm/year, with a maximum of 180 mm/year. This area presents with substantial terrain fluctuations, evident differences in rainfall, extremely complex climate, and low vegetation coverage [4]. It is also located through four north-south fault zones, and tectonic activity is strong, factors which collectively result in the breakage of the slope and frequent occurrence of landslides. The Yuka-Xiangdui area (Figure 6c) is mainly affected by glaciers, weathering, and crustal uplift, which collectively lead to the wide distribution of the loose deposits and bare bedrock, resulting in the occurrence of large-scale landslides [5]. The Guza-Tianquan section (Figure 6d) belongs to the eastern edge of the Sichuan-Tibet line and is exposed to strong fault and earthquake activities. Therefore, numerous landslides have occurred in this region. The Yongba-Juebo section (Figure 6e) is strongly subjected to erosion and is cut by the Nujiang River. The ice-snow cover area is large and is affected by the subtropical warm and humid climate, and the freeze-thaw effect is remarkable. Landslides occur relatively easily in this environment. Generally, these five sections warrant attention, with a focus on continuous deformation monitoring. In addition, the inventory of 17 typical landslides marked in Figures 6 and 7, which may pose a great threat to railway construction, is shown in Appendix A, Table A1.

Time-Series Monitoring for Typical Landslides
Based on the landslide identification results, certain high-risk landslides can be monitored using time-series deformation. Therefore, we selected two landslide bodies, G and H, as shown in Figure 9, both of which demonstrate an evident continuous deformation trend. As shown in Figure 9b, the maximum cumulative deformation of landslide G reaches −400 mm, and a continuous deformation trend exists at all three points. This landslide is located on the right bank of the Jinsha River, and there are multiple tensile cracks in the back edge of the slope, and findings indicate that the range of cracks gradually expands with the continuous increase in deformation. Therefore, human activities on the slope, river cuts, and rainfall may cause failure of this slope. If sliding occurs, it may block the Jinsha River and endanger the upstream bridge, a phenomenon which may also induce failure of other unstable slopes and may expand the scale of damage. Therefore, continuous monitoring is recommended in the future. Landslide H is an old landslide, and there is less vegetation cover on this slope. Optical imaging results reveal that it is a freeze-thaw landslide. According to the time-series deformation of this landslide, the maximum cumulative deformation of P3 on this landslide can reach −296 mm (Figure 9d), and the three points demonstrate a state of continuous deformation, a phenomenon which will also pose a grave threat to the construction and operation of nearby buildings and structures. It could be inferred from the time-series deformation of these two disaster points that although ALOS-2 data exhibited certain disadvantages in terms of time resolution, the overall cumulative deformation trend could also be studied appropriately; therefore, such landslides can be monitored effectively by considering InSAR time-series deformation.

Effectiveness of the Hot Spot Analysis Method
To verify the effectiveness of hot spot analysis, we considered the annual deformation rate containing many error regions as an example to explore the performance of

Effectiveness of the Hot Spot Analysis Method
To verify the effectiveness of hot spot analysis, we considered the annual deformation rate containing many error regions as an example to explore the performance of hot spot analysis in extracting effective deformation. This area is located in Tianquan County and many errors exist in the annual deformation rate, owing to the high vegetation coverage. Figure 10a shows the original annual deformation rate. Figure 10b shows the deformation area obtained only by removing data on the stable region from Figure 10a using a deformation threshold, which contained many discontinuous deformation points. Interpretation of landslides directly from the information presented in Figure 10a,b is difficult. Figure 10c shows the results of extracting data on the effective deformation area by hot spot analysis. It could be inferred that most error points were eliminated and the effective deformation area was remarkably clear; thus, landslide interpretation could be performed rapidly and accurately, whereby data on 15 effective landslide areas were effectively extracted in this area (Figure 10c). Figure 10d shows the difference between Figure 10b,c, which indicates the error deformation area removed via hot spot analysis. It could be inferred that the hot spot analysis could be performed to remove many discrete and invalid deformation points that caused significant interference to interpretation.

Controlling Factors of Spatial Landslide Distribution
To determine the main factors controlling the spatial landslide distribution along the Sichuan-Tibet Railway, the elevation, slope, aspect, river, fault zone, seismic peak acceleration, and rainfall were selected and their relationships with the spatial landslide distribution were analyzed. In the vertical direction, the relationships between landslide and elevation, slope, and aspect were analyzed using spatial analysis statistical tools. As shown in Figure 11a, landslide disasters gradually increased with an increase in elevation,

Controlling Factors of Spatial Landslide Distribution
To determine the main factors controlling the spatial landslide distribution along the Sichuan-Tibet Railway, the elevation, slope, aspect, river, fault zone, seismic peak acceleration, and rainfall were selected and their relationships with the spatial landslide distribution were analyzed. In the vertical direction, the relationships between landslide and elevation, slope, and aspect were analyzed using spatial analysis statistical tools. As shown in Figure 11a, landslide disasters gradually increased with an increase in elevation, and decreased markedly when the elevation exceeded 5000 m. Therefore, elevation is an important factor affecting spatial landslide distribution [43]. Additionally, 77% of the landslide hazard points were located at altitudes ranging from 3000-5000 m. Figure 11b depicts the relationship between the landslide points and the slope. It could be inferred that a high slope almost developed along the rivers. Statistical analysis showed that 81% of the landslides were distributed on slopes of 10-40 • . The reason for this phenomenon is that low slopes do not present with sufficiently dynamic conditions, and high slopes are not conducive to the accumulation of landslide materials [44]. Therefore, the slope in the middle was beneficial for the occurrence of landslides. Furthermore, the terrain slope affects the rainfall and heat distribution, thereby affecting vegetation coverage. Figure 11c illustrates the relationship between landslides and the slope aspect. According to the statistical diagram, it could be inferred that the landslides were mostly developed in the northeastern and eastern areas, accounting for 53% of the total landslides. This may be attributable to the presence of less sunshine in these two directions, underdeveloped vegetation, and poor slope stability, thereby rendering these areas more prone to landslides [45].
In the horizontal direction, the relationships between landslides and fault zone, seismic peak acceleration, and rainfall were statistically analyzed. There are many rivers present along the line, and river cutting and lateral erosion seriously affect the slope stability. After seven buffer zones were created on the river, the relationship between the spatial landslide distribution and the distance from the river was analyzed. The number of landslides decreased with an increased distance from the river (Figure 12a), and the highest landslide density was observed to be in the area 300-600 m from the river (i.e., 0.91 landslides/km). Meanwhile, these areas are considered more suitable for human habitation; therefore, they have also demonstrated vulnerability due to human activities [43]. In general, there are close relationships between rivers and landslides. The fault zone changes the rock mass structure and causes crushing, which is conducive to rainwater infiltration and induces natural disasters. Figure 12b shows the statistical analysis of the distance between the landslide points and the fault zone. The number of landslides showed a slowly decreasing trend, with an increase in the distance between the landslide and the fault zone, a finding which was consistent with the theory and proved that the fault zone exerted a certain influence on landslide development. As one of the main triggering factors of landslides, earthquakes can not only directly induce landslide failure, but can also reduce the threshold of landslide failure by reducing slope stability [20,46]. Figure 12c shows that the landslide density reached a substantial value when the seismic peak acceleration was 0.3 g and 0.4 g, a finding which proved that it affected the spatial landslide distribution. Precipitation is another important factor inducing landslides [47]. In this study, information on annual average precipitation (GPM) [48] from 2016-2021 was used to explore the relationship between landslides and precipitation. It could be inferred from Figure 12d that rainfall was more distributed in the east and less in the west. When the rainfall was >700 mm/year, the landslide density gradually increased with increasing rainfall. When it reached a value >1000 mm/year, the landslide density reached a maximum of 0.98 landslides/km. When the rainfall was <700 mm/year, it exerted minimal effect on landslide development, and the landslide in the region was dominated by other factors; therefore, there were no evident statistical relationships with rainfall. In general, in the horizontal direction, the seismic peak acceleration, rainfall, river, and fault zone correlated with the landslide density, which largely determined the spatial landslide distribution along the Sichuan-Tibet Railway.
fects the rainfall and heat distribution, thereby affecting vegetation coverage. Figure 11c illustrates the relationship between landslides and the slope aspect. According to the statistical diagram, it could be inferred that the landslides were mostly developed in the northeastern and eastern areas, accounting for 53% of the total landslides. This may be attributable to the presence of less sunshine in these two directions, underdeveloped vegetation, and poor slope stability, thereby rendering these areas more prone to landslides [45]. In the horizontal direction, the relationships between landslides and fault zone, seismic peak acceleration, and rainfall were statistically analyzed. There are many rivers present along the line, and river cutting and lateral erosion seriously affect the slope stability. The above-mentioned analysis shows that landslides along the Sichuan-Tibet Railway are mainly affected by the following factors: (1) the river erosion cutting and river cutting will render the slope steep and unstable; therefore, most landslides in the region are affected by this factor and mainly occur on both sides of the river. (2) Fault zones and earthquakes: there are many fault zones along the Sichuan-Tibet Railway, and seismic activity is frequent. These activities destroy the integrity of the rock mass, enable easy infiltration of the rainwater, and lead to the production of a weak sliding surface, resulting in the occurrence of landslides. Especially in the east, landslides are concentrated near the Baqing-Leiniaoqi, Lancang River, Batang, Xianshui River, and Wenchuan-Maowen fault zones. (3) Change in freeze-thaw: the analysis shows that landslides may occur frequently at elevations of 3000-5000 m. The temperature difference between day and night is considerable, and repeated freezing and thawing may weaken the slope integrity. (4) Precipitation: owing to the simultaneous rainfall and heat along the line, extreme precipitation is common [49], and ice-snow melting water and rainfall will increase the slope weight and sliding force, phenomena which will also lead to the occurrence of landslides.

Conclusions
In the present study, using multi-temporal InSAR methods and 164 scenes of ALOS-2 ascending data, preliminary detection of landslide hazards in the Linzhi-Ya'an section of the Sichuan-Tibet Railway was conducted. After obtaining the annual deformation rate, the hot spot analysis method was introduced to rapidly extract data on the effective deformation. Then, the optical images and other geological factors were combined for landslide interpretation. A total of 517 landslide points were detected, which posed a certain threat to railway construction and safe operation. Large-scale landslides were mainly concentrated in Lulang-Guxiang, Yuka-Xiangdui, Yongba-Juebo, Jinsha River coast, and Guza-Tianquan County. By comparing the results of the original InSAR and the hot spot analysis, it was found that the hot spot analysis method could be applied to extract data on the effective deformation area, which considerably improved the speed of visual interpretation. Most importantly, it reduced the number of misjudgments and missing judgments to a certain extent. Finally, according to the detection results of landslides, the correlations between spatial landslide distribution and geological factors (such as slope, aspect, elevation, fault zone, river, and precipitation) were explored. The spatial landslide distribution was strongly and positively correlated with the distance between rivers and fault zones, and the main influencing factors were rainfall, river, freeze-thaw weathering, and earthquakes. This study provides important information that may be considered as vital technical support for the construction and safe operation of the Linzhi-Ya'an section of the Sichuan-Tibet Railway. Although the reliability of the hot spot analysis method has been proved in this study, it is limited by the determination of distance and stable deformation thresholds. Generally, these two thresholds are defined by human experience, which may result in the missing judgments of efficient deformation. Thus, accurate determination of distance and stable deformation thresholds will be studied in our subsequent research.