Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Low-Complexity One-Bit DOA Estimation for Massive ULA with a Single Snapshot
Previous Article in Journal
Accuracy Evaluation and Analysis of GNSS Tropospheric Delay Inversion from Meteorological Reanalysis Data
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Wide Area Detection and Distribution Characteristics of Landslides along Sichuan Expressways

1
College of Geological Engineering and Geomatics, Chang’an University, Xi’an 710054, China
2
Big Data Center for Geosciences and Satellites, Chang’an University, Xi’an 710054, China
3
Key Laboratory of Western China’s Mineral Resource and Geological Engineering, Ministry of Education, Xi’an 710054, China
4
School of Highway, Chang’an University, Xi’an 710064, China
5
Engineering Research Center of Highway Infrastructure Digitalization, Ministry of Education, Xi’an 710064, China
6
Sichuan Expressway Construction and Development Group, Chengdu 610041, China
7
College of Transportation Engineering, Chang’an University, Xi’an 710064, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2022, 14(14), 3431; https://doi.org/10.3390/rs14143431
Submission received: 17 May 2022 / Revised: 2 July 2022 / Accepted: 13 July 2022 / Published: 17 July 2022
(This article belongs to the Section Remote Sensing in Geology, Geomorphology and Hydrology)

Abstract

:
Wide area landslide detection is a major international research hotspot in the field of geological hazards, and the integration of multi-temporal optical satellite images and spaceborne interferometric synthetic aperture radar (InSAR) appears to be an effective way to realize this. In this paper, a technical framework is presented for wide area landslide detection: (i) multi-temporal satellite optical images are used to detect landslides with distinguishable geomorphological features; (ii) Generic Atmospheric Correction Online Service (GACOS) assisted InSAR stacking is employed to generate annual surface displacement rate maps in radar line of sight using satellite SAR images from both ascending and descending tracks, which are in turn utilized to automatically detect active landslides from ground motion using hotspot analysis, and (iii) the distribution characteristics of the detected landslides are investigated by examining their relationships with topographic and hydrological factors. Three expressways in Sichuan Province, China—namely the Yakang (Ya’an-Kangding), Yaxi (Ya’an-Xichang), and Lushi (Luding-Shimian) expressways—and their surrounding regions (a total area of approximately 20,000 square kilometers) were chosen as the study area. A total of 413 landslides were detected, among which 320 were detected using multi-temporal satellite optical images, and 109 were detected using GACOS-assisted InSAR stacking. It should be noted that only 16 landslides were detected by both approaches; these landslides all exhibited not only obvious geomorphological features but also ground motion. A statistical analysis of the topographic and hydrological factors shows that of the detected landslides: 81% are distributed at elevations of 1000–2500 m, over 60% lie within the elevation range of 100~400 m, and 90% present with medium and steep slopes (20°~45°), and 80% are located within areas seeing an annual rainfall of 950~1050 mm. Nine landslides were found to pose potential safety hazards to the expressways. The research findings in this paper have directly benefitted the Sichuan expressways; equally important, it is believed that the technical framework presented in this paper will provide guidance for hazard mitigation and the prevention of transportation hazards in the future.

1. Introduction

In 2020, according to the Ministry of Natural Resources of the People’s Republic of China’s statistical report, 4810 landslide hazards occurred in mainland China, accounting for 61.3% of all geological disasters (http://www.mnr.gov.cn/dt/dzdc (accessed on 23 October 2021)). Landslides have become a prominent contributor to China’s geological hazards. The term “landslide” describes the downward movement of a mass of rock, debris or earth along a slope, which can shape the landform, destroy mountains and provide sediment for river networks [1,2,3]. They pose a serious threat to people, properties and the environment in many areas. Many researchers are committed to wide area landslide detection, providing essential data support for landslide hazard mitigation and prevention [3,4,5].
Landslides are mainly caused by geophysical and meteorological factors, including heavy rain, long-term rainfall, wet weather, slope erosion, earthquake vibration, volcanic activity, rapid melting of snow and ice, human activity and so on [3,6]. In some cases, before a landslide becomes destabilized, the whole slope or part of the slope moves slowly, and the moving speed can reach several centimeters per year. Such ground motion is an early indicator of slope instability and is also key to mapping landslide-prone areas and establishing a landslide early warning system [7]. Currently, the commonly used methods for landslide detection include multi-temporal satellite optical image interpretation technology, spaceborne repeat-pass Synthetic Aperture Radar Interferometry (InSAR), light detection, and ranging (LiDAR), and geophysical prospecting.
In the 1970s, satellite optical image interpretation technology was adopted for landslide detection and mapping [8]. An increasing number of high spatial resolution satellite optical images are available alongside improved computer performance and this technique has been widely used in wide area landslide detection [9,10]. However, it has inherited limitations. On the one hand, satellite optical images are often limited by the presence of clouds and fog, which leads to the failure of landslide detection in some regions. On the other hand, whilst this approach can better detect landslides with obvious geomorphological features or apparent signs of deformation, its detection of landslides with no obvious signs of surface deformation is far more challenging [10]. Interpreted signs of landslides are usually less than 1 m in the spatial dimension. Therefore, the spatial resolution needs to be greater than 1 m when selecting satellite optical images.
As a significant earth observation technology, spaceborne repeat-pass InSAR technology has been extensively used in the field of surface deformation monitoring, and especially landslide detection and deformation monitoring [11,12,13,14], because of its high precision, high-resolution, all-day, all-weather useability, alongside its enormous catalogue of historical data (since the early 1990s). As such, it can effectively detect landslides in inaccessible regions. Although InSAR technology plays a significant role in landslide detection and monitoring, it nevertheless presents certain challenges such as interference decoherence, atmospheric interference, geometric distortion of squint imaging, terrain effects, and insensitivity to north–south ground motion [12]. Previous work has shown that the globally open Sentinel-1 imagery, with its wide coverage and short revisiting time, is suitable for wide area landslide detection. Combined with the ascending and descending tracks of Sentinel-1 images, the potential for misjudgment of a landslide caused by the geometric distortion of SAR images can be significantly reduced [12,13,15,16].
Airborne LiDAR is a remote sensing technology that can be used to acquire a digital terrain model (DTM) with sub-meter accuracy [17]. Through the global navigation satellite system and flight navigation information, the geographic position of the airborne instruments can be reconstructed accurately. Compared with satellite optical images, airborne LiDAR has incomparable advantages: (1) it can penetrate a forest canopy and obtain terrain elevation information [17], and (2) historical damage such as ancient landslides, large deposits, and cracks can be interpreted using LiDAR data [18]. Therefore, this technology is suitable for landslide detection in vegetation-covered areas [19,20]. It should be pointed out that the acquisition cost of airborne LiDAR data is high, and as such, it is largely used only for the detailed investigation of key areas and significant landslide hazards at present [18].
Geophysical prospecting technology is one non-destructive geophysical method that can be used to detect landslides, which utilizes both natural and artificially created physical fields, (e.g., gravitational, magnetic, electrical, seismic, and radiometric fields) acquired on the surface of the Earth, in the air and/or underground (in wells and shafts) to determine the location of geological structures, ore bodies, and so forth and their fundamental characteristics [21,22]. However, the method has certain errors and is costly to implement and is highly affected by geological conditions.
The transition zone from the Qinghai–Tibet Plateau to the Sichuan Basin is characterized by clouds and fogs, dense vegetation, high elevation differences, and strong water vapor variations. Therefore, wide area landslide detection in this region is challenging. In this paper, a technical framework is presented for wide-area landslide detection. First, multi-temporal satellite optical images are used to detect landslides with distinguishable geomorphological features and qualitatively assess the threat to highways. Secondly, Sentinel-1 images from both ascending and descending tracks are interferometrically processed and then analyzed with GACOS-assisted InSAR stacking to generate annual surface displacement rate maps in the radar line of sight (LOS). Deforming landslides are automatically extracted by the hotspot analysis method. Finally, the main conditions that control and influence the development and distribution of landslides are analyzed through the statistics of features of all the detected landslides, such as elevation, elevation difference, area, slope angle, slope aspect, rainfall, and other characteristics. The results will be directly applied to landslide hazard mitigation and prevention of hazards along Sichuan expressways.

2. Materials and Methods

2.1. Study Area

The research area of this study encompasses approximately 20,000 km2 along and around the Yakang, Yaxi, and Lushi expressways in Sichuan Province, China (Figure 1), located in the transitional zone between the Sichuan Basin and the Qinghai–Tibet Plateau. The area comprises a large topographic terrain, complex geomorphological types, and complicated geological conditions. Land cover is dominated by forests, grasslands, and shrubs (Figure 2a). Because of the strong influence of rainfall, earthquakes, and human engineering activity, landslides often occur along and near the expressways in the study area. For instance, at 03:50 on 21 August 2020, the Zhonghaicun landslide occurred due to heavy rainfall, and road traffic was interrupted for a long period [23], which seriously threatened the construction and safe operation of the expressways. Therefore, it is necessary to carry out wide area landslide detection along and around the expressways in the region as a matter of urgency.

2.2. Data

Multi-source satellite optical images and Sentinel-1 radar images were used in this study, as shown in Table 1. The selected satellite optical images included Gaofen-2 satellite images and Google Earth archived images. The Gaofen-2 satellite was successfully launched on 19 August 2014. The spatial resolution of the panchromatic band is 0.8 m, the multi-spectral resolution is 3.2 m, the width is 45 km, and the revisiting time is five days, meeting the requirements of landslide detection. The archived images of Google Earth are mainly satellite images with a spatial resolution below 1 m, such as WorldView, Quickbird, and GeoEye, which provide abundant data resources for landslide detection. The European Space Agency’s Sentinel-1A/B satellites operate day and night, performing C-band microwave SAR imaging and providing radar imagery with wide coverage, (e.g., 250 km × 250 km) and a short repeat cycle (12 days for a single satellite and six days for both). In this study, a total of 521 Sentinel-1A/B images in four tracks (Table 1) were selected for the single-look complex (SLC) images of VV polarization in the interferometric wide swath (IW).

2.3. Methods

2.3.1. Classification of Landslides

Landslides can be classified into different types in terms of different criteria [24]. In the case of landslide identification, Li et al. (2022) proposed to divide landslides into four types in terms of their histories and movement: actively deforming slopes (I), reactivated historically deformed slopes (II), stabilized historically deformed slopes (III), and undeformed but potentially unstable slopes (IV) [25] (Figure 2). The actively deforming slopes (I) refer to the slopes without historical failures which are actively deforming but there has been no historical slope failure; the reactivated historically deformed slopes (II) refer to the slopes with historical failures that have been reactivated, (i.e., in motion); the stabilized historically deformed slopes (III) refer to the slopes with historical failures that are in a stable state, and the undeformed but potentially unstable slopes (IV) refer to the slopes without historical failures that are in a stable state but with a slide-prone structure [25]. Types I and II can often be identified through mapping their deformation signals since they are in motion, and geomorphologic features are usually utilized to identify Types II and III since they exhibit historical failures. In contrast, it is difficult (if not impossible) to identify Type IV since there is neither a deformation signal nor a geomorphologic feature with such slopes.

2.3.2. Multi-Temporal Satellite Optical Image Interpretation

The multi-temporal satellite optical image interpretation method has the characteristics of wide range, non-contact, low cost, and high precision, and is suitable for wide area landslide detection [18]. Utilizing spectral and textural differences in satellite optical images acquired at different times, combined with topographic and geomorphologic features, it is straightforward to identify actively deforming slopes with obvious deformation signals (Figure 3a) and historically deformed slopes (Figure 3b,c) [10]. The actively deforming slopes with obvious deformation signals can exhibit cracks at tops (heads) of the slope and small-scale collapses on the lower parts of the slope [2,26]. Because of historical failures, the historically deformed slopes often exhibit certain geomorphological features, such as scarps, flanks, cracks, ridges, etc. [25,26]. Hence, such geomorphological features can be used to identify the abovementioned landslides. Note that satellite optical images can be seriously affected by the presence of clouds, fogs, and their corresponding shadows on the ground, which in turn might result in false identification of landslides. To address this, multi-temporal satellite optical images are often employed to ensure there are sufficient cloud-free images to examine the evolution of the slopes.

2.3.3. InSAR Automatic Landslide Detection

InSAR technology can be employed to generate information on surface elevation and surface deformation from the interference signals of two or more SAR images covering the same area [27,28]. When the microwave signal from the SAR satellite sensor passes through the atmosphere, the propagation path will be affected to varying degrees with time. Any differences in the atmosphere at the time of SAR image acquisition during different time periods will directly affect the InSAR observations. This influence is usually known as the atmospheric effect, which is also one of the most significant error sources in InSAR deformation monitoring [29,30]. Most landslides in the study area are located in mountainous areas, where topography-related tropospheric delays are particularly significant, and atmospheric disturbances are not negligible where there are rivers, lakes, and other bodies of water nearby, or strong convective weather is present during satellite transits. Therefore, prior to landslide detection, it is necessary to carry out a reasonable and effective atmospheric correction on interference images. The GACOS (http://www.gacos.net/ (accessed on 25 October 2021)) system developed by Prof. Zhenhong Li’s team is based on HERS ECMWF data with a spatial resolution of 0.125° and a temporal resolution of six hours. It uses the Iterative Tropospheric Decomposition (ITD) model to separate the stratification and turbulence signals from the total tropospheric delay, and can provide near-real-time, global coverage of InSAR atmospheric delay error correction mapping for atmospheric delay error correction in InSAR measurements [31,32,33].
Considering the large topographic relief in the study area, the presence of strong convective weather, and the obvious atmospheric influence on InSAR interferometric images, the GACOS-assisted InSAR stacking method is an efficient means to generate annual surface displacement rate maps [34,35]. In this study, a multi-looking factor of 4 in range, (i.e., a resolution of approximately 17 m in the ground direction) and 1 in azimuth, (i.e., a resolution of approximately 14 m in the ground direction) were used to process the Sentinel-1 SAR images. The SAR images were first pre-processed to obtain single-looking complex SLC images, and the topographic influence was attenuated using the 30 m resolution ALOS Global Digital Surface Model (AW3D30) released by the Japan Aerospace Exploration Agency (JAXA). A time baseline of 60 days and a spatial baseline of 150 m were set to generate interferograms, and a series of unwrapping maps were obtained by adaptive filtering and phase unwrapping processing, after which atmospheric correction products of GACOS were introduced to attenuate atmospheric effects (Figure 4). It should be pointed out that the average coherence greater than 0.4 can be observed in 80% of the study area, suggesting the feasibility of InSAR for this area (Figure 2b).
The InSAR stacking technique is based on the assumption that surface deformation varies approximately linearly and the average rate per pixel is estimated from N interferometric images [36]. This has the advantages of low computational resource requirements, straightforward automated processing, simple, effective, and fast application, and the effective suppression of random errors. The least square method was used to solve the average unwrapping rate following GACOS correction, the interference phases were weighted by time intervals and, finally, the annual surface deformation rate map of LOS direction in the study area was generated [35], as in Equation (1).
V m e a n = i = 1 N φ i Δ T i i = 1 N Δ T i 2
where V mean is the mean velocity, φ i demonstrates the unwrapped phase and Δ T i is the temporal baseline of the ith interferogram.
The Getis–Ord G i * statistic ( G i * )—also known as hot-spot analysis [37,38]—is a method commonly used to analyze location-related tendencies (clustering) among the attributes of spatial data (points or areas). The method is an adaptation of the General G-statistic [37], a global method for quantifying the degree of spatial autocorrelation over an area. It is first necessary to determine whether a feature point has significant spatial clustering characteristics by looking at each pixel in the environment of adjacent elements within a certain distance. Pixels with significant spatial clustering characteristics have high value attributes and are surrounded by other elements that also have high values [33,34]. It is time-consuming and laborious work to visually interpret the annual average surface deformation rates of a large area in order to ascertain locations of landslide risk. Therefore, the Getis–Ord hotspot analysis method was introduced to automatically extract effective deforming areas from the annual surface deformation rate map. The calculation formula is as follows:
G i * ( d ) = j = 1 n ω i j ( d ) x j X ¯ j = 1 n ω i j ( d ) S n j = 1 n ω i j 2 ( j = 1 n ω i j ( d ) ) 2 n 1
with X ¯ = j = 1 n x j n and S = j = 1 n x j 2 n ( X ¯ ) 2 where G i * ( d ) is computed for feature (i) at a distance (d) standardized as a z-score. x j is the annual surface deformation value of each neighbor pixel. ω i j is the spatial weight for the target neighbor i and j pair and n is the total number of pixels. The Euclidean method was used to calculate the geographical distances from each feature to its neighboring features.
The degree of clustering of deformation points is assessed based on the z-score and p-Value. The higher the z value, the greater the likelihood of clustering at a particular location. In this study, deformation points satisfying p < 0.01 (99% confidence level) and z > ± | 2.58 | were selected as having significant spatial clustering properties [39].

2.3.4. Landslide Identification with Integrated Remote Sensing Techniques

After obtaining geomorphological and deformational information as demonstrated in the previous two sections, landslides can be identified (Table 2). Type I landslides (actively deforming slopes) are detected mainly on the basis of deformation signals; since some of them may exhibit geomorphological features such as cracks, scarps and collapses, multi-temporal satellite optical image interpretation may provide assistance in some cases. Type II landslides (reactivated historically deformed slopes) should be detected on the basis of deformation signals and geomorphological features, whilst Type III landslides (stabilized historically deformed slopes) should be detected by interpreting geomorphological features with optical remote sensing images. Type IV landslides (undeformed but potentially unstable slopes) have no deformation signal or geomorphological feature, and neither InSAR nor optical remote sensing imagery can be utilized to detect Type IV landslides. Such landslides can be identified according to their geological structures, and geophysical prospecting technology should be used for this purpose [21,22]. Geophysical prospecting technology is costly to implement, and its performance is highly affected by geological conditions. Hence, it is still difficult to detect Type IV landslides over wide regions at present.

3. Results and Analysis

3.1. Annual Surface Deformation Rate

The GACOS-assisted InSAR stacking method was used to process 521 SAR images acquired during the period from 2017 to 2021 to generate LOS annual deformation rate maps from both ascending and descending tracks encompassing appropriately 40,000 km2 of three expressways in Sichuan, as shown in Figure 5a,b, respectively. The annual surface deformation rate maps show that most of the landslides threatening the expressway are concentrated in the Luding–Shimian, Shimian–Hanyuan and Hanyuan–Yingjing sections, as shown in the three red areas in Figure 5, with the landslides in the area 1 (Figure 6a) being significantly obvious. A comparison of the location of the landslide concentrations with hydropower stations shows that the landslide concentrations are all in the vicinity of the upstream of the hydropower stations. The hydropower station in the middle of Zones 1 and 2 is the Pubugou Hydropower Station, with a dam height of 186 m and a reservoir capacity of 5.39 billion m3. Close to Zone 3 is the Dagangshan Hydropower Station, with a dam height of 210 m and a reservoir capacity of 0.072 million m3.
The large capacity of the Pubugou Hydropower Station is likely to be a particularly significant factor in the landslides in Zone 1. There are 16 landslides in Zones 1, 2 and 3 with an annual surface deformation rate greater than 20 mm/a, representing great safety hazards to the expressway (Figure 6). From Figure 6a,d, it can be seen that more landslides are detectable in the ascending track images than in the descending track images in Zone 1, whereas landslides A and H are detectable in both the ascending and descending track images. However, landslides B, C, D, E, F, G and I are only detectable in the ascending track images and landslide J is only detectable in the descending track images. The situation in Zones 2 and 3 is similar to that in Zone 1, where there are landslides that could be detected by both ascending and descending track images and by descending or ascending track images only.

3.2. Landslide Detection and Mapping

The combined multi-temporal satellite optical image visual interpretation and GACOS-assisted InSAR stacking methods detected a total of 413 landslides. As shown in Figure 7a, 320 landslides were detected through satellite optical image visual interpretation. A total of 109 landslides were automatically detected by the hotspot analysis method based on InSAR-derived LOS annual surface displacement rate maps, and the hotspot analysis results for Zones 1 and 2 are shown in Figure 8. It is worth noting that only 16 landslides with obvious deformation signals and geomorphological features were detected by both approaches, mainly because the landslide types that can be detected by the two methods are different [40]. Then, according to the classification of landslides in this study, the 413 detected landslides were divided into 371 actively deforming slopes, 23 reactivated historically deformed slopes and 19 stabilized historically deformed slopes (Figure 7b).
Where the landslides that can be detected using high-resolution satellite optical images include landslides that have occurred and landslides that are incubating with obvious deformation signals, and since the spatial resolution of the satellite optical images used in this paper was 0.8 m, landslides without geomorphological features cannot be detected by this method. In contrast, the GACOS-assisted InSAR stacking method can only be used to detect landslides in motion, (i.e., Types I and II). Although the two approaches can detect different types of landslides, they are nevertheless both landslides and, therefore, when it comes to landslide detection, it is necessary to combine the two technologies in order that they may complement and verify each other.
Landslide cataloging must determine not only the location of any landslide but also its extent. Both the geomorphological features of the satellite optical images and the annual surface deformation rate generated by the GACOS-assisted InSAR stacking method can be used as a basis for mapping the extent of landslides, as shown in Figure 9. Four cases are included in this research, combining the two methods. In case 1, landslides can be detected by both approaches, as shown in Figure 9a, and mapping the extent of landslides is based on the range of annual surface deformation rate map, and geomorphological features taken from satellite optical images. In case 2, landslides can only be identified by multi-temporal satellite optical image interpretation, as shown in Figure 9b, and the extent of the mapped landslide is determined based on high-resolution satellite images over multiple periods. In case 3, landslides can only be detected by InSAR, as shown in Figure 9c, and mapping the extent of landslides requires a combination of deformation extent, multi-temporal satellite optical images, DEM, slope and other multi-source data. In case 4, as shown in Figure 9d, neither of these two methods can detect landslides; this is also the main area in which landslides are missed due to the extent of dense vegetation. In these areas, the geographical features are affected by vegetation, meaning that satellite optical images cannot be interpreted, and InSAR results are incoherent or geometrically distorted. In such areas, airborne LiDAR is believed to be an effective tool to detect old landslides covered by dense vegetation [40].

3.3. Inventory of Landslides That Pose a High Risk to Expressways

By scrutinizing the multi-temporal satellite optical images, annual surface displacement rate map and distance from the expressway, it was concluded that nine of the 413 landslides detected were potential threats to the expressway, of which three, four and two were safety hazards to Lushi, Yaxi and Yakang expressways, respectively. The length, width, area, distance from the expressway, satellite optical image texture characteristics, and annual surface displacement rate of the landslides were combined to qualitatively evaluate these nine landslides, as detailed in Table 3.

4. Discussion

4.1. Impact Factors on Results of Landslide Mapping

4.1.1. Effects of Spaceborne SAR Imaging Geometry

Due to the inherent oblique observation geometry of all imaging radar systems, surface slopes and similar terrain features lead to geometric distortions of the data acquired by SAR systems [12], the most relevant of which are layovers and shadows. These phenomena are particularly evident in high mountain valley areas and will result in parts of the detection area being left without a valid signal, ultimately leading to missed landslide detection [41]. In order to better compensate for the geometric distortions caused by SAR imaging geometry, this study combined the ascending and descending tracks of SAR images to obtain the deformation characteristics of detection targets under different imaging geometries, thus effectively reducing the areas that would not be effectively observed due to geometric distortions.
Quantitatively evaluating the combined ascending and descending tracks of SAR images can effectively reduce geometric distortions. In this paper, the geometric distortion and effective observable area were counted, as shown in Figure 10. The effective observation area made up 89.4% of the total area of the ascending track, the layover area represented 10.3% and the shadow area was less than 0.3%. The observable area in the descending track comprised 83.7% of the total area, the layover area made up 15.8% and the shadow area was just 0.5%. In the combined ascending and descending imagery, up to 97.3% of the region could be observed effectively, and the geometric distortion region was reduced to 2.7%. This result was consistent with the analysis that combined observation of ascending and descending tracks along the Sichuan–Tibet Railway, effectively reducing the single-track geometric distortion area by 30% to 1.5% [16]. These results show that the area of geometric distortion can be effectively reduced by using the joint observation of the ascending and descending tracks, whereupon the rate of missed detections of landslides can be reduced.

4.1.2. Identifying Deformation without Landslides

The area of deformation in the LOS annual surface deformation rate map includes not only landslides that are deforming but also alluvial fans formed when tributaries enter the main river channel, subsidence caused by agricultural irrigation, urban ground subsidence, and deformation caused by human engineering activities, as shown in Figure 11a–d, respectively. Certain cases, such as Zone I in Figure 8c and Zones I, S, and T in Figure 8d, will also be considered valid deformation areas when hotspot analysis is performed on the LOS annual surface deformation rate map. To effectively detect landslides, it is necessary to combine DEM and satellite optical images for comparison and validation, to exclude any misclassification of automatic landslide detection by hotspot analysis.

4.2. Development Characteristics of Landslides along Sichuan Expressways

According to the landslide mapping method proposed in this paper, the surface area of the 413 landslides detected amounted to 47.32 km2 (the total area of the study area is about 20,000 km2), with minimum and maximum areas of 1176 m2 and 6.980 km2, respectively.

4.2.1. Distribution Pattern of Detected Landslides with Respect to Topographic Factors

Rugged mountain terrain is more prone to landslides than flat terrain [42]. In previous studies, topographic parameters, such as elevation, relative altitude, slope angle, and aspect, were considered essential factors for the formation, development, and occurrence of landslides [42], and were widely used to evaluate the distribution of landslides [42,43]. The AW3D30 data at 30 m spatial resolution was employed to generate topographic factors, such as elevation, slope angle, slope aspect, and relative altitude, as shown in Figure 12.
Elevation plays an important role in landslide distribution because it is related to land cover, roads and river valleys, slope material, and the seismic wave amplification effect [42]. From Figure 7 and Figure 13a, it can be seen that approximately 80.5% of landslides occur between 1000 and 2500 m above sea level, and the maximum landslide density is distributed between 1000 and 2000 m. Below 1500 m, the number of landslides increases with elevation, while above 1500 m, the number decreases with elevation. Because slopes above 1500 m are less affected by rainfall and river erosion, the number and development of landslides are relatively low. As shown in Figure 13f, more than 77% of the landslides are located in the area with a relative altitude below 400 m, and the densest area of landslides is located in the area with relative altitudes of 100–300 m. When the relative altitude is below 200 m, the number of landslides increases with increased relative altitude. However, when the relative altitude is higher than 200 m, the number of landslides decreases with increased relative altitude.
Generally, a larger slope angle will lead to a larger slope sliding force, making it more susceptible to landslides [44]. Slope angles for the whole study area were divided into 12 categories: <10°, 10°~15°, 15°~20°, 20°~25°, 25°~30°, 30°~35°, 35°~40°, 40°~45°, 45°~50°, 50°~55°, 55°~60° and >60°. Statistical calculation and correlation results are shown in Figure 13b. The study area was found to be dominated by medium to steep slope angles (20–45°), accounting for about 89.5% of the whole study area. In the range of slope angles from 30 to 40°, the landslides are the most intensive (54.3%). At slope angles below 40°, the number of landslides increases as the slope angle increases, and at slope angles greater than 40°, the number of landslides decreases as the slope angle increases.
The slope aspect may be the only topographic parameter with orientation in common topographic distribution analysis, and so it is usually combined with other directional parameters for further examination [42]. Previous studies have shown that the slope aspect influences weathering conditions, temperature, precipitation, deposition conditions, and land cover, and so plays a role in the susceptibility of an area to landslides [45,46]. As shown in Figure 12b, the slope aspect was divided into the following 16 classes: north (N, ~348.76°–360° and ~0°–11.25°), north-northeast (NNE, ~11.26°–33.75°), northeast (NE, ~33.76°–56.25°), east-northeast (ENE, ~56.26°–78.75°), east (E, ~78.76°–101.25°), southeast (SE, ~101.26°–123.75°), east-southeast (ESE, ~123.76°–146.25°), south-southeast (SSE, ~146.26°–168.75°), south (S, ~168.76°–191.25°), south-southwest (SSW, ~191.26°–213.75°), southwest (SW, 213.76°~236.25°), west-southwest (WSW, ~236.26°–258.75°), west (W, ~258.76°–281.25°), west-northwest (WNW, ~281.26°–303.75°), northwest (NW, ~303.76°–326.25°) and north-northwest (NNW, 326.26°~348.75°). As shown in Figure 13c, the slope aspects ENE, E, ESE, SE, SSE, SSW, SW, WSW, and W together account for 82.3% of the total number of landslides—consistent with the main directions of the high mountain valleys in the study area: ESE, SSE, and SW.

4.2.2. Distribution Pattern of Detected Landslides with Respect to Hydrological Factors

On an annual scale, precipitation is widely recognized as the driver of the seasonal acceleration and deceleration of slow-moving landslides [47,48,49]. However, precipitation may not be the only factor to initiate a landslide that has been constantly active for hundreds of years [50]. A comparison was made of 20-year average precipitation (2001–2020) relative to the landslide distribution (Figure 12c). Overall, 80% of the identified landslides are located in areas with an annual rainfall of 950 to 1050 mm/a. When rainfall is less than 960 mm/a, the number of landslides increases with rainfall, and when rainfall is greater than 960 mm/a, the number of landslides decreases with rainfall. In addition, the identified landslides were compared with the average excess precipitation between 2017 and 2021 (Figure 12d). Here, excess precipitation is defined as the difference between annual precipitation and the 20-year average, where a positive value indicates excess. As shown in Figure 13e, it can be found that the landslides being deformed (as identified by InSAR from 2017 to 2021) increased with excess rainfall, indicating that rainfall is one of the driving factors of landslide occurrence. However, regions with less rainfall may experience more landslide activity, partly depending on other factors such as bedrock formation, land uplift, and tectonic activity [51].

4.2.3. Four General Rules of Landslide Distribution

By summarizing the spatial distribution patterns of landslides in this study and other areas (Table 4), four general rules have been determined. The first is the general rule of the distribution of elevation and landslides. There is an elevation boundary that makes the number of landslides reach the peak at that elevation boundary, whereupon the phenomenon of decreasing above that elevation boundary will occur. For example, the elevation boundary of the China–Pakistan Karakoram Highway is 2100 m [52]; along the Jinsha River Basin, the elevation boundary is 2000 m [13]. With increased altitude, due to human activity, river erosion, and rock weathering, the number of landslides is small and the degree of development is slow. The second rule determines that landslides mainly occur in areas with a relative altitude below 500 m, but they are also related to the research area. For instance, in the Bailong River Basin, landslides are most densely distributed in areas with a relative elevation difference of 50 to 150 m [53]. In the Jinsha River Basin, landslides are most densely distributed in areas with a relative elevation difference of 500 to 700 m [13]. The third rule determines that landslides are most frequently developed in the slope range 10°–45° because, as the slope angle increases, the rock weathering layer subsequently becomes thinner or non-existent [54] and shear stress increases [55], resulting in a lower number of landslides on steep slopes. The fourth rule determines that landslides tend to occur in mountainous areas with abundant rainfall. For example, on the west coast of the US, 75% of landslides are located in mountainous areas with relatively abundant rainfall (≥2000 mm/a) [56]. However, there are many exceptions. For example, along the China–Pakistan Karakoram Highway area, the rainfall range is 162–941 mm/a, with 84.3% of landslides occurring in areas in the semi-arid region (<400 mm/a) [52]; along the Sichuan–Tibet Railway Ya’an-Linzhi). No significant statistical relationship emerged between rainfall and landslide distribution [39].

5. Conclusions

In this paper, a technical framework for wide area landslide detection is proposed and successfully applied to an area of approximately 20,000 km2 along and near three expressways in Sichuan, leading to the following conclusions:
(1)
A total of 320 landslides were detected by interpreting GF-2 images and Google Earth historical images. GACOS-assisted InSAR stacking was employed to obtain annual surface displacement rate maps in the radar LOS direction using satellite SAR images from both ascending and descending tracks, which were in turn utilized to automatically detect 109 potential landslides with ground motion using hotspot analysis. The combination of the two techniques resulted in the detection of 413 landslides covering 47.32 km2, of which 371 landslides were actively deforming slopes (I), 23 landslides were reactivated historically deformed slopes (II), and 19 landslides were stabilized historically deformed slopes (III).
(2)
Combined observations from Sentinel-1 ascending and descending tracks can effectively increase the visible area from 89.4% and 83.7% to 97.3%, to potentially significantly reduce the number of missed landslides.
(3)
According to the statistics of the topographical and hydrological factors of the 413 detected landslides, the main distribution area of landslides have the following characteristics: elevation of 1000–2500 m, relative altitude of 100–400 m, slope angle of 20°–45°, slope aspect of 101°–168° and 191°–281°, and annual rainfall of 950–1050 mm.
Although the technical framework can significantly reduce the rate of missed landslides, some landslides continue to be missed, (e.g., those with large gradient deformations, those in densely vegetated areas, mini-landslides, and undeformed but potentially unstable slopes). To solve this problem, combining various techniques and data is necessary. For example, high-resolution SAR or satellite optical image pixel offset tracking measurement technology can be used to detect landslides with large gradient deformation, LiDAR can be used to detect Type III landslides in densely vegetation areas, high-resolution SAR images can obtain annual surface deformation rates for mini-landslides detection, and (semi-)airborne geophysical prospecting technology can detect undeformed but potentially unstable slopes. Finally, it should be noted that, although the hotspot analysis method can automatically extract effective deformation areas, it presents a variety of situations requiring manual identification. Nevertheless, in the future, it will be improved to realize the automatic detection of landslides.

Author Contributions

Conceptualization, Z.L. (Zhenhong Li) and B.C.; methodology, B.C. and Z.L. (Zhenhong Li); validation, B.C., Z.L. (Zhenhong Li), C.Z. (Chenglong Zhang), B.H., J.D. and Y.C.; formal analysis, B.C., Z.L. (Zhenhong Li) and C.Z. (Chenglong Zhang); writing original draft preparation, B.C.; writing review and editing, Z.L. (Zhenhong Li), C.Z. (Chenglong Zhang), M.D., W.Z., S.Z. (Shuangcheng Zhang), and C.Z. (Chi Zhang); resources, Z.L. (Zhenhong Li), Z.L. (Zhiyong Liao), S.Z. (Shuanke Zhou), J.W. and J.P.; visualization, B.C.; supervision, Z.L. (Zhenhong Li) and J.P.; funding acquisition, Z.L. (Zhenhong Li). All authors have read and agreed to the published version of the manuscript.

Funding

The research was funded by the National Key Research and Development Program of China (2020YFC1512000), the National Natural Science Foundation of China Major Program (41941019), the Shaanxi Province Science and Technology Innovation team (2021TD–51), the Shaanxi Province Geoscience Big Data and Geohazard Prevention Innovation Team (2022), and the Fundamental Research Funds for the Central Universities, CHD (Refs. 300102260301 and 300102261108).

Data Availability Statement

The Sentinel-1 datasets were freely provided by Copernicus and ESA, https://search.asf.alaska.edu/#/ (accessed on 21 March 2022). The ALOS digital elevation model (DEM) can be found at https://www.eorc.jaxa.jp/ALOS/en/aw3d30/data/index.htm (accessed on 21 March 2022), and precipitation data were freely downloaded from https://pmm.nasa.gov/data-access/downloads/gpm (accessed on 21 March 2022).

Acknowledgments

We gratefully thank anonymous reviewers and editors for their constructive comments and suggestions. We also thank the China Resources Satellite Application Centre for providing the GF-2 images and the Google Earth platform for providing the optical remote sensing images.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Cruden, D.M.; Varnes, D.J. Landslide types and processes. In Landslides: Investigation and Mitigation; Turner, A.K., Schuster, R.L., Eds.; Special Report; National Academy Press: Washington, DC, USA, 1996; pp. 36–75. [Google Scholar]
  2. Hungr, O.; Serge, L.; Luciano, P. The Varnes classification of landslide types, an update. Landslides 2013, 11, 167–194. [Google Scholar] [CrossRef]
  3. Mondini, A.C.; Guzzetti, F.; Chang, K.T.; Monserrat, O.; Martha, T.R.; Manconi, A. Landslide failures detection and mapping using Synthetic Aperture Radar: Past, present and future. Earth Sci. Rev. 2021, 217, 103611. [Google Scholar] [CrossRef]
  4. Maes, J.; Parra, C.; Mertens, K.; Bwambale, B.; Jacobs, L.; Poesen, J.; Dewitte, O.; Vranken, L.; Hontheim, A.D.; Kabaseke, C.; et al. Questioning network governance for disaster risk management: Lessons learnt from Landslide risk management in Uganda. Environ. Sci. Policy 2018, 85, 163–171. [Google Scholar] [CrossRef]
  5. Raska, P. Contextualizing community-based landslide risk reduction: An evolutionary perspective. Landslides 2019, 16, 1747–1762. [Google Scholar] [CrossRef]
  6. Schumm, S.A. Geomorphic Thresholds: The Concept and Its Applications. Trans. Inst. Br. Geogr. 1979, 4, 485–515. [Google Scholar] [CrossRef]
  7. Van, A.; Th, W.J. Creep processes in landslides. Earth Surf. Process. Landf. 1984, 9, 573–583. [Google Scholar] [CrossRef]
  8. Sauchyn, D.J.; Trench, N.R. Landsat applied to landslide mapping. Photogramm. Eng. Remote Sens. 1978, 44, 735–741. [Google Scholar]
  9. Guzzetti, F.; Mondini, A.C.; Cardinali, M.; Fiorucci, F.; Santangelo, M.; Chang, K.T. Landslide inventory maps: New tools for an old problem. Earth Sci. Rev. 2012, 112, 42–66. [Google Scholar] [CrossRef] [Green Version]
  10. Lu, H.Y.; Li, W.L.; Xu, Q.; Dong, X.J.; Dai, C.; Wang, D. Early detection of landslides in the upstream and downstream areas of the Baige Landslide, the Jinsha River based on optical remote sensing and InSAR technologies. Geomat. Inform. Sci. Wuhan Univ. 2019, 44, 1342–1354. (In Chinese) [Google Scholar] [CrossRef]
  11. Zhang, L.; Liao, M.S.; Dong, J.; Xu, Q.; Gong, J.Y. Early Detection of Landslide Hazards in Mountainous Areas of West China Using Time Series SAR Interferometry–A Case Study of Danba, Sichuan. Geomat. Inform. Sci. Wuhan Univ. 2018, 43, 2039–2049. (In Chinese) [Google Scholar] [CrossRef]
  12. Li, Z.H.; Song, C.; Yu, C.; Xiao, R.Y.; Chen, L.F.; Luo, H.; Dai, K.R.; Ge, D.Q.; Ding, Y.; Zhang, Y.X.; et al. Application of satellite radar remote sensing to landslide detection and monitoring: Challenges and solutions. Geomat. Inform. Sci. Wuhan Univ. 2019, 44, 967–979. (In Chinese) [Google Scholar] [CrossRef]
  13. Liu, X.J.; Zhao, C.Y.; Zhang, Q.; Lu, Z.; Li, Z.H.; Yang, C.S.; Zhu, W.; Zeng, J.L.; Chen, L.Q.; Liu, C.J. Integration of Sentinel-1 and ALOS/PALSAR-2 SAR datasets for mapping active landslides along the Jinsha River corridor, China. Eng. Geol. 2021, 284, 106033. [Google Scholar] [CrossRef]
  14. Zhang, C.L.; Li, Z.H.; Yu, C.; Song, C.; Xiao, R.Y.; Peng, J.B. Landslide Detection of the Jinsha River Region Using GACOS Assisted InSAR Stacking. Geomat. Inform. Sci. Wuhan Univ. 2021, 46, 1649–1657. (In Chinese) [Google Scholar] [CrossRef]
  15. Raspini, F.; Bianchini, S.; Ciampalini, A.; Soldato, M.D.; Solari, L.; Novali, F.; Conte, S.D.; Rucci, A.; Ferretti, A.; Casagli, N. Continuous, semi-automatic monitoring of ground deformation using Sentinel-1 satellites. Sci. Rep. 2018, 8, 7253. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Dai, K.; Zhang, L.L.; Song, C.; Li, Z.L.; Zhuo, G.X.; Xu, Q. Quantitative Analysis of Sentinel-1 Imagery Geometric Distortion and Their Suitability along Sichuan-Tibet Railway. Geomat. Inform. Sci. Wuhan Univ. 2021, 46, 1450–1460. (In Chinese) [Google Scholar] [CrossRef]
  17. Jaboyedoff, M.; Oppikofer, T.; Abellan, A.; Derron, M.H.; Loye, A.; Metzger, R.; Pedrazzini, A. Use of LIDAR in landslide investigations: A review. Nat. Hazards 2012, 61, 5–28. [Google Scholar] [CrossRef] [Green Version]
  18. Xu, Q. Understanding and Consideration of Related Issues in Early Identification of Potential Geohazards. Geomat. Inform. Sci. Wuhan Univ. 2020, 45, 4–12. (In Chinese) [Google Scholar] [CrossRef]
  19. Booth, A.M.; Roering, J.J.; Perron, J.T. Automated landslide mapping using spectral analysis and high-resolution topographic data: Puget Sound lowlands, Washington, and Portland Hills, Oregon. Geomorphology 2009, 109, 132–147. [Google Scholar] [CrossRef]
  20. Guo, C.; Xu, Q.; Dong, X.J.; Liu, X.S.; She, J.X. Geohazard Recognition by Airborne LiDAR Technology in Complex Mountain Areas. Geomat. Inform. Sci. Wuhan Univ. 2021, 46, 1538–1547. (In Chinese) [Google Scholar] [CrossRef]
  21. Xiong, S. The Strategic Consideration of the Development of China’s Airborne Geophysical Technology. Geol. China 2009, 36, 1366–1374. (In Chinese) [Google Scholar]
  22. Langhammer, L.; Rabenstein, L.; Schmid, L.; Bauder, A.; Grab, M.; Schare, P.; Maurer, H. Glacier Bed Surveying with Helicopter-Borne Dual-polarization Ground-penetrating Radar. J. Glaciol. 2019, 65, 123–135. [Google Scholar] [CrossRef] [Green Version]
  23. Ye, F.; Fu, W.X.; Zhou, H.F.; Liu, Y.; Ba, R.J.; Zheng, S. The “8·21” rainfall-induced Zhonghaicun landslide in Hanyuan County of China: Surface features and genetic mechanisms. Landslides 2021, 18, 3421–3434. [Google Scholar] [CrossRef]
  24. Highland, L.M.; Bobrowsky, P. The Landslide Handbook—A Guide to Understanding Landslides; U.S. Geological Survey: Reston, VA, USA, 2008; pp. 1325–1454.
  25. Li, Z.H.; Zhang, C.L.; Chen, B.; Zhan, J.W.; Ding, M.T.; Lü, Y.; Li, X.L.; Peng, J.B. A Technical Framework of Landslide Prevention Based on Multi-Source Remote sensing and Its Engineering Application. Earth Sci. 2022, 45, 1901–1916. (In Chinese) [Google Scholar]
  26. Li, W.L.; Xu, Q.; Lu, H.Y.; Dong, X.J.; Zhu, Y.Q. Tracking the Deformation History of Large-Scale Rocky Landslides and Its Enlightenment. Geomat. Inform. Sci. Wuhan Univ. 2019, 44, 1043–1053. (In Chinese) [Google Scholar] [CrossRef]
  27. Zebker, H.A.; Goldstein, R.M. Topographic mapping from interferometric synthetic aperture radar observations. J. Geophysic. Res. 1986, 91, 4993–4999. [Google Scholar] [CrossRef]
  28. Massonnet, D.; Rossi, M.; Carmona, C.; Adragna, F.; Peltzer, G.; Feigl, K.; Rabaute, T. The displacement of the Landers earthquake mapped by radar interferometry. Nature 1993, 364, 138–142. [Google Scholar] [CrossRef]
  29. Li, Z.H.; Muller, J.P.; Cross, P. Interferometric synthetic aperture radar (InSAR) atmospheric correction: GPS, Moderate Resolution Imaging Spectroradiometer (MODIS), and InSAR integration. J. Geophysic. Res. 2005, 110, B03410. [Google Scholar] [CrossRef]
  30. Li, Z.H.; Fielding, E.J.; Cross, P.; Preusker, R. Advanced InSAR atmospheric correction: MERIS/MODIS combination and stacked water vapour models. Int. J. Remote Sens. 2009, 30, 3343–3363. [Google Scholar] [CrossRef]
  31. Yu, C.; Penna, N.T.; Li, Z.H. Generation of real-time mode high-resolution water vapor fields from GPS observations. J. Geophys. Res. Atmos. 2017, 122, 2008–2025. [Google Scholar] [CrossRef]
  32. Yu, C.; Li, Z.H.; Penna, N.T. Interferometric synthetic aperture radar atmospheric correction using a GPS-based iterative tropospheric decomposition model. Remote Sens. Environ. 2018, 204, 109–121. [Google Scholar] [CrossRef]
  33. Yu, C.; Li, Z.H.; Penna, N.T.; Crippa, P. Generic Atmospheric Correction Model for Interferometric Synthetic Aperture Radar Observations. J. Geophys. Res. Solid Earth 2018, 123, 9202–9222. [Google Scholar] [CrossRef]
  34. Yu, C.; Li, Z.H.; Penna, N.T. Triggered afterslip on the southern Hikurangi subduction interface following the 2016 Kaikōura earthquake from the InSAR time series with atmospheric corrections. Remote Sens. Environ. 2020, 251, 112097. [Google Scholar] [CrossRef]
  35. Xiao, R.Y.; Yu, C.; Li, Z.H.; Song, C.; He, X.F. General Survey of Large-scale Land Subsidence by GACOS-Corrected InSAR Stacking: Case Study in North China Plain. Proc. Int. Assoc. Hydrol. Sci. 2020, 382, 213–218. [Google Scholar] [CrossRef] [Green Version]
  36. Sandwell, D.T.; Price, E.J. Phase gradient approach to stacking interferograms. J. Geophys. Res. 1998, 103, 30183–30204. [Google Scholar] [CrossRef] [Green Version]
  37. Getis, A.; Ord, J.K. The Analysis of Spatial Association by Use of Distance Statistics. Geogr Anal. 1992, 24, 189–206. [Google Scholar] [CrossRef]
  38. Ord, J.K.; Getis, A. Local Spatial Autocorrelation Statistics: Distributional Issues and an Application. Geogr Anal. 1995, 27, 286–306. [Google Scholar] [CrossRef]
  39. Zhang, J.M.; Zhu, W.; Cheng, Y.Q.; Li, Z.H. Landslide Detection in the Linzhi-Ya’an Section along the Sichuan-Tibet Railway Based on InSAR and Hot Spot Analysis Methods. Remote Sens. 2021, 13, 3566. [Google Scholar] [CrossRef]
  40. Xu, Q.; Lu, H.Y.; Li, W.L.; Dong, X.J.; Guo, C. Types of Potential Landslide and Corresponding Identification Technologies. Geomat. Inform. Sci. Wuhan Univ. 2022, 47, 377–387. (In Chinese) [Google Scholar] [CrossRef]
  41. Dai, K.R.; Tie, Y.B.; Xu, Q.; Feng, Y.; Zhuo, G.C.; Shi, X.L. Early identification of potential landslide geohazards in alpine-canyon terrain based on SAR interferometry–A case study of the middle section of Yalong river. J. Radars 2020, 9, 554–568. (In Chinese) [Google Scholar]
  42. Wu, C.H.; Peng, C.; Yu, S.L.; Ayala, I.A.; Chao, H.; Jian, Y.S. Seismogenic fault and topography control the spatial patterns of landslides triggered by the 2017 Jiuzhaigou earthquake. J. Mt. Sci. 2018, 15, 793–807. [Google Scholar] [CrossRef]
  43. Zhang, S.; Li, R.; Wang, F.W.; Lio, A. Characteristics of landslides triggered by the 2018 Hokkaido Eastern Iburi earthquake, Northern Japan. Landslides 2019, 16, 1691–1708. [Google Scholar] [CrossRef]
  44. Cui, P.; Ling, Z.X.; Qiang, Z. Risk assessment of highways affected by debris flows in Wenchuan earthquake area. J. Mt. Sci. 2013, 10, 173–189. [Google Scholar] [CrossRef]
  45. Papathanassiou, G.; Valkaniotis, S.; Ganas, A.; Pavlides, S. GIS-based statistical analysis of the spatial distribution of earthquake-induced landslides in the island of Lefkada, Ionian Islands, Greece. Landslides 2012, 10, 771–783. [Google Scholar] [CrossRef]
  46. Kamp, U.; Growley, B.J.; Khattak, G.A.; Owen, L.A. GIS-based landslide susceptibility mapping for the 2005 Kashmir earthquake region. Geomorphology 2008, 101, 631–642. [Google Scholar] [CrossRef]
  47. Bennett, G.L.; Roering, J.J.; Mackey, B.H.; Handwerger, A.L.; Schmidt, D.A.; Guillod, B.P. Historic drought puts the brakes on earth flows in Northern California. Geophys. Res. Lett. 2016, 43, 5725–5731. [Google Scholar] [CrossRef] [Green Version]
  48. Handwerger, A.L.; Fielding, E.J.; Huang, M.H.; Bennett, G.L.; Liang, C.; Schulz, W.H. Widespread Initiation, Reactivation, and Acceleration of Landslides in the Northern California Coast Ranges due to Extreme Rainfall. J. Geophys. Res.-Earth. 2019, 124, 1782–1797. [Google Scholar] [CrossRef] [Green Version]
  49. Xu, Y.K.; Kim, J.; George, D.L.; Lu, Z. Characterizing Seasonally Rainfall-Driven Movement of a Translational Landslide using SAR Imagery and SMAP Soil Moisture. Remote Sens. 2019, 11, 2347. [Google Scholar] [CrossRef] [Green Version]
  50. Roering, J.J.; Mackey, B.H.; Handwerger, A.L.; Booth, A.M.; Schmidt, D.A.; Bennett, G.L.; Darriau, C.C. Beyond the angle of repose: A review and synthesis of landslide processes in response to rapid uplift, Eel River, Northern California. Geomorphology 2015, 236, 109–131. [Google Scholar] [CrossRef] [Green Version]
  51. Xu, Y.K.; Schulz, W.H.; Lu, Z.; Kim, J.; Baxtrom, K. Geologic controls of slow-moving landslides near the US West Coast. Landslides 2021, 18, 3353–3365. [Google Scholar] [CrossRef]
  52. Su, X.J.; Zhang, Y.; Meng, X.M.; Yue, D.X.; Ma, J.H.; Guo, F.Y.; Zhao, Z.Q.; Rehman, M.U.; Khalid, Z.; Chen, G.; et al. Landslide mapping and analysis along the China-Pakistan Karakoram Highway based on SBAS-InSAR detection in 2017. J. Mt. Sci. 2021, 18, 2540–2564. [Google Scholar] [CrossRef]
  53. Li, Y.X.; Zhang, Y.; Su, X.J.; Zhao, F.M.; Liang, Y.W.; Meng, X.M.; Jia, J. Early identification and characteristics of potential landslides in the Bailong River Basin using InSAR technique. Natl. Remote Sens. Bulletin. 2021, 25, 677–690. (In Chinese) [Google Scholar] [CrossRef]
  54. Ietto, F.; Perri, F.; Cella, F. Weathering Characterization for Landslides Modeling in Granitoid Rock Masses of the Capo Vaticano Promontory (Calabria, Italy). Landslides 2018, 15, 43–62. [Google Scholar] [CrossRef]
  55. Saha, A.K.; Gupta, R.P.; Arora, M.K. GIS-Based Landslide Hazard Zonation in the Bhagirathi (Ganga) Valley, Himalayas. Int J. Remote Sens. 2002, 23, 357–369. [Google Scholar] [CrossRef]
  56. Conforti, M.; Ietto, F. Influence of Tectonics and Morphometric Features on the Landslide Distribution: A Case Study from the Mesima Basin (Calabria, South Italy). J. Earth Sci. 2019, 32, 393–409. [Google Scholar] [CrossRef]
Figure 1. Location of the three expressways and coverage of image datasets. Pink rectangles indicate the coverage of Sentinel-1 ascending tracks, green rectangles show the coverage of Sentinel-1 descending tracks, and red rectangle indicates the coverage of satellite optical images. Upper right corner red rectangle indicates the location of the study area in China.
Figure 1. Location of the three expressways and coverage of image datasets. Pink rectangles indicate the coverage of Sentinel-1 ascending tracks, green rectangles show the coverage of Sentinel-1 descending tracks, and red rectangle indicates the coverage of satellite optical images. Upper right corner red rectangle indicates the location of the study area in China.
Remotesensing 14 03431 g001
Figure 2. Suitability analysis of InSAR: (a) land cover, and (b) map of the average coherence during the period from October 2017 to July 2021.
Figure 2. Suitability analysis of InSAR: (a) land cover, and (b) map of the average coherence during the period from October 2017 to July 2021.
Remotesensing 14 03431 g002
Figure 3. Examples of landslides with different types: (a) actively deforming slope (I), (b) reactivated historically deformed slope (II), (c) stabilized historically deformed slope (III), and (d) undeformed but potentially unstable slope (IV).
Figure 3. Examples of landslides with different types: (a) actively deforming slope (I), (b) reactivated historically deformed slope (II), (c) stabilized historically deformed slope (III), and (d) undeformed but potentially unstable slope (IV).
Remotesensing 14 03431 g003
Figure 4. Flowchart for landslide detection and analysis.
Figure 4. Flowchart for landslide detection and analysis.
Remotesensing 14 03431 g004
Figure 5. LOS annual surface deformation rate map along the Yakang, Yaxi, and Lushi expressways: (a) ascending track, and (b) descending track. Negative value (red color) and positive value (blue color) indicate that the Earth’s surface is moving away from and towards radar sensor, respectively.
Figure 5. LOS annual surface deformation rate map along the Yakang, Yaxi, and Lushi expressways: (a) ascending track, and (b) descending track. Negative value (red color) and positive value (blue color) indicate that the Earth’s surface is moving away from and towards radar sensor, respectively.
Remotesensing 14 03431 g005
Figure 6. Key areas of landslides: (a,d) represent LOS annual deformation rates in Zone 1 from ascending and descending tracks, respectively; (b,e) represent LOS annual deformation rates in Zone 2 from ascending and descending tracks, respectively, and (c,f) represent LOS annual deformation rates in Zone 3 from ascending and descending tracks, respectively. Red rectangles indicate a deforming landslide threatening the expressway.
Figure 6. Key areas of landslides: (a,d) represent LOS annual deformation rates in Zone 1 from ascending and descending tracks, respectively; (b,e) represent LOS annual deformation rates in Zone 2 from ascending and descending tracks, respectively, and (c,f) represent LOS annual deformation rates in Zone 3 from ascending and descending tracks, respectively. Red rectangles indicate a deforming landslide threatening the expressway.
Remotesensing 14 03431 g006
Figure 7. Distribution of landslides along the three expressways in Sichuan Province, China detected by multi-temporal satellite optical images and GACOS-assisted InSAR stacking. (a) Distribution map of landslides detected by different methods, and (b) distribution map of actively deforming slopes (I), reactivated historically deformed slopes (II), and stabilized historically deformed slopes (III). Pink circles indicate historical earthquakes that occurred during the period from 1900 to 2021, collected from the United States Geological Survey.
Figure 7. Distribution of landslides along the three expressways in Sichuan Province, China detected by multi-temporal satellite optical images and GACOS-assisted InSAR stacking. (a) Distribution map of landslides detected by different methods, and (b) distribution map of actively deforming slopes (I), reactivated historically deformed slopes (II), and stabilized historically deformed slopes (III). Pink circles indicate historical earthquakes that occurred during the period from 1900 to 2021, collected from the United States Geological Survey.
Remotesensing 14 03431 g007
Figure 8. Hotspot analysis results: (a,b) represent hotspot analysis results of Zone 1 from ascending and descending tracks, respectively, and (c,d) represent hotspot analysis results of Zone 2 from ascending and descending tracks, respectively.
Figure 8. Hotspot analysis results: (a,b) represent hotspot analysis results of Zone 1 from ascending and descending tracks, respectively, and (c,d) represent hotspot analysis results of Zone 2 from ascending and descending tracks, respectively.
Remotesensing 14 03431 g008
Figure 9. Landslide identification from four cases: (a) landslides can be detected by multi-temporal satellite optical image interpretation method and GACOS assisted InSAR stacking method; (b) landslides can only be identified by satellite optical image interpretation; (c) landslides can only be detected by InSAR, and (d) landslides are not detectable by either method. (i) Google optical images; (ii) LOS annual surface deformation rate maps from ascending tracks, and (iii) LOS annual surface deformation rate maps from descending tracks. Solid red line means landslide is detectable and dashed red line means landslide is not detectable.
Figure 9. Landslide identification from four cases: (a) landslides can be detected by multi-temporal satellite optical image interpretation method and GACOS assisted InSAR stacking method; (b) landslides can only be identified by satellite optical image interpretation; (c) landslides can only be detected by InSAR, and (d) landslides are not detectable by either method. (i) Google optical images; (ii) LOS annual surface deformation rate maps from ascending tracks, and (iii) LOS annual surface deformation rate maps from descending tracks. Solid red line means landslide is detectable and dashed red line means landslide is not detectable.
Remotesensing 14 03431 g009
Figure 10. Geometric distortion (layovers and shadows) analysis.
Figure 10. Geometric distortion (layovers and shadows) analysis.
Remotesensing 14 03431 g010
Figure 11. Exemplar areas in motion with misclassified landslides: (a) alluvial fan, (b) agricultural irrigation, (c) ground subsidence, and (d) human engineering activities.
Figure 11. Exemplar areas in motion with misclassified landslides: (a) alluvial fan, (b) agricultural irrigation, (c) ground subsidence, and (d) human engineering activities.
Remotesensing 14 03431 g011
Figure 12. Spatial relationships between landslide distribution and slope angle, slope aspect, precipitation and excess precipitation: (a) spatial relationship between slope angle and detected landslides; (b) spatial relationship between slope aspect and detected landslides; (c) spatial relationship between precipitation and detected landslides, and (d) spatial relationship between excess precipitation and detected landslides.
Figure 12. Spatial relationships between landslide distribution and slope angle, slope aspect, precipitation and excess precipitation: (a) spatial relationship between slope angle and detected landslides; (b) spatial relationship between slope aspect and detected landslides; (c) spatial relationship between precipitation and detected landslides, and (d) spatial relationship between excess precipitation and detected landslides.
Remotesensing 14 03431 g012
Figure 13. Relationships between distribution of detected landslides and elevation, slope angle, slope aspect, precipitation, excess precipitation, and relative altitude, respectively: (a) relationship between distribution of detected landslides and elevation; (b) relationship between distribution of detected landslides and slope angle; (c) relationship between distribution of detected landslides and slope aspect; (d) relationship between distribution of detected landslides and precipitation; (e) relationship between distribution of detected landslides and excess precipitation, and (f) relationship between distribution of detected landslides and relative altitude.
Figure 13. Relationships between distribution of detected landslides and elevation, slope angle, slope aspect, precipitation, excess precipitation, and relative altitude, respectively: (a) relationship between distribution of detected landslides and elevation; (b) relationship between distribution of detected landslides and slope angle; (c) relationship between distribution of detected landslides and slope aspect; (d) relationship between distribution of detected landslides and precipitation; (e) relationship between distribution of detected landslides and excess precipitation, and (f) relationship between distribution of detected landslides and relative altitude.
Remotesensing 14 03431 g013
Table 1. SAR and optical images used in this study.
Table 1. SAR and optical images used in this study.
SatelliteOrbital DirectionPathNumber of ImagesRange/Azimuth Pixel SpacingTemporal Coverage
Gaofen-2--280.8 m20180725, 20181125
20190422, 20200827
20201109, 20210201
20210112, 20190207
Google Earth---<1 m2010–2020
Sentinel-1Ascending261142.3 m/13.9 m20171015–20210726
2611420171015–20210726
12811520171010–20210721
Descending62932.3 m/13.9 m20180602–20210728
1358520180619–20210721
Table 2. Landslide classification and identification methods.
Table 2. Landslide classification and identification methods.
Landslide TypesExamplesFeaturesIdentification Methods
Actively deforming slopes (I) Remotesensing 14 03431 i001Deformation, possibly with geomorphological features, e.g., cracks, scarps and collapsesDeformation monitoring techniques including InSAR and SAR/optical pixel offset tracking (POT), assisted by multi-temporal satellite optical image interpretation
Reactivated historically deformed slopes (II) Remotesensing 14 03431 i002Deformation together with clear geomorphological features, e.g., cracks, scarps and collapsesDeformation monitoring techniques including InSAR and POT, plus multi-temporal satellite optical image interpretation
Stabilized historically deformed slopes (III) Remotesensing 14 03431 i003Clear geomorphological features, e.g., cracks, scarps and collapsesMulti-temporal satellite optical image interpretation
Undeformed but potentially unstable slopes (IV) Remotesensing 14 03431 i004No deformation signal or geomorphological featureGeophysical prospecting technology
Table 3. Statistics of nine detected landslides with potential risks to expressways.
Table 3. Statistics of nine detected landslides with potential risks to expressways.
ExpresswayNoLandslidesLongitudeLatitudeLength/mWidth/mArea/km2Maximum LOS Rate (mm/a)Distance from Expressway/kmComprehensive Judgement
Lushi1Pangou102.329.8350023008.42−935This is an old landslide and deformation is visible locally.
2Shanshugang102.129.510686940.68−840This is an unstable slope, which is in deformation stage at present; the expressway was built in front of this unstable slope.
3Huangcaoping102.229.59668310.88−830.3This old landslide is currently in deformation stage.
Yaxi4Waduanshan102.729.48904270.31−730.8This unstable slope is currently in deformation stage.
5Jianshibao102.429.29596450.49−561.5This potential landslide was caused by mining; there are continuous cracks in the back wall of the slope and signs of human activity on the front edge of the landslide, which could lead to overall instability.
6Haizi102.629.4408515877.521020.14This is an old landslide. Both the leading edge and trailing edges are deformed, which may happen again.
7Wazhaping102.629.424656101.35−911.0This is an old landslide. Both the leading edge and middle edges are deformed, which may happen again.
Yakang8Lianghekou102.430.03816860.19−570.07This is an unstable slope body in deformation stage, caused by human activity.
9Tuanjie102.2301782260.037−560This is a deforming landslide, which was locally deformed and eroded by the river; it has the possibility of overall instability.
Table 4. Statistics on the spatial distribution of detected landslides in other areas.
Table 4. Statistics on the spatial distribution of detected landslides in other areas.
NumberElevation/mRelative Altitude/mSlope Angle/°Slope Aspect/°Precipitation/mmRegionReference Articles
1550–2100
(81.7%)
-20–45
(87.7%)
E, SE, S
(22%, 20%, 17%)
<400(162–941)
(84.3%)
China-Pakistan Karakoram Highway[52]
21000–3000
(69%)
<900
(57%)
10–45
(83%)
--Jinsha River corridor, China[13]
3-<300
(76.3%)
20–40
(83.3%)
SE, S, SW-Bailong River Basin, China[53]
43000–5000
(77%)
10–40
(83%)
N, NE, E, SE
(82%)
No statistically significant relationshipSichuan–Tibet Railway, China[39]
52200–3700
(85%)
200–50010–40
(81.2%)
ENE~SSE
(68.8%)
-Jiuzhaigou, northwest Sichuan Province, China[42]
6-160–320
(75%)
20–45
(70%)
--Calabria, South Italy[56]
7--5–30
(97.1%)
->2000(0–4000)
(75%)
US West Coast[51]
81000~2500
(80.5%)
<400
(77%)
20–45
(89.5%)
ESE, SSE, SW950~1050
(80%)
Southwest Sichuan Province, ChinaThis study
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chen, B.; Li, Z.; Zhang, C.; Ding, M.; Zhu, W.; Zhang, S.; Han, B.; Du, J.; Cao, Y.; Zhang, C.; et al. Wide Area Detection and Distribution Characteristics of Landslides along Sichuan Expressways. Remote Sens. 2022, 14, 3431. https://doi.org/10.3390/rs14143431

AMA Style

Chen B, Li Z, Zhang C, Ding M, Zhu W, Zhang S, Han B, Du J, Cao Y, Zhang C, et al. Wide Area Detection and Distribution Characteristics of Landslides along Sichuan Expressways. Remote Sensing. 2022; 14(14):3431. https://doi.org/10.3390/rs14143431

Chicago/Turabian Style

Chen, Bo, Zhenhong Li, Chenglong Zhang, Mingtao Ding, Wu Zhu, Shuangcheng Zhang, Bingquan Han, Jiantao Du, Yanbo Cao, Chi Zhang, and et al. 2022. "Wide Area Detection and Distribution Characteristics of Landslides along Sichuan Expressways" Remote Sensing 14, no. 14: 3431. https://doi.org/10.3390/rs14143431

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