Abstract

The analysis of landslide evolution using archived optical remote-sensing images is common, but it is often limited by the acquisition frequency, cloud cover, and resolution. With the development of space-borne radar observation technology, small baseline subset interferometric synthetic aperture radar (SAR) technology provides a new technical approach for detecting landslide deformation before disasters. The Sentinel-1A SAR datasets (20170219–20210330) were used to study the time-series deformation characteristics of the Wangjiashan landslide in the Baihetan Reservoir area before its impoundment. The time-series results show that the Wangjiashan landslide was in an initial deformation state in the prior four years, and the deformation first occurred in the middle part and then expanded to the landslide toe and crown retrogressive movement characteristics. Combined with an analysis of field deformation signs, these findings suggest that the upper landslide mass formed a local sliding surface, which caused serious deformation of the road. An analysis of historical rainfall data revealed that the Wangjiashan landslide is sensitive to rainfall, and the deformation is not only significantly correlated with cumulative rainfall but also influenced by concentrated heavy rainfall. The research in this study provides a basis for the monitoring, early warning, and risk management of the Wangjiashan landslide during the impoundment period. This work also provides a useful reference for investigations using space-borne SAR Earth observations in geological disaster prevention and control.

1. Introduction

The experience from the construction and impoundment of many large hydropower plants has shown that reservoir impoundment induces a large number of nascent landslide hazards while leading to increasing deformation of ancient landslides, resulting in many casualties and property losses [14]. Therefore, it is particularly important to understand the historical deformation of ancient landslides in reservoir inundation areas and then evaluate their stability. Due to the lack of monitoring equipment, the analyses of the deformation of landslides have mainly been performed using archived optical remote-sensing images [5, 6]. Optical remote sensing is characterized by wide coverage and a large spectral range but a long acquisition time interval; moreover, the images are affected by clouds and fog, and the spatial resolution is inadequate.

In recent years, with the development of Earth observation technology, interferometric synthetic aperture radar (InSAR) technology has provided a new technical way to study the historical deformation of landslides. The InSAR technique uses two or more synthetic aperture radar (SAR) images from the same area and the same orbit to perform interferometric differential processing to remove residual terrain errors, atmospheric effects, and other related errors to obtain high-precision three-dimensional information on the terrain and surface deformations. Compared with traditional ground measurement technology and space-borne optical remote sensing, InSAR technology has the characteristics of less weather influence, large-scare coverage, and high measurement accuracy. It has played an important role in landslide deformation monitoring [714] and large-scale landslide detection [1518]. With the continuous improvement of hardware and software of space-borne radar Earth observation technology, InSAR technology will play an even more important role in the future [1922]. However, due to the variability of signal propagation and scattering characteristics, the measurement results of conventional differential InSAR (DInSAR) methods are mainly affected by systematic and random errors, i.e., satellite orbit errors, atmospheric delays, elevation errors, and surface scatterer variations. These may lead to temporal and spatial discorrelation and reduce the deformation monitoring accuracy [2325]. In addition, these phenomena may also lead to the inability to extract surface deformation information in low-coherence regions. In contrast, the time-series InSAR technique can well compensate for the shortcomings of the traditional DInSAR technique, i.e., it uses a certain error reduction method to process multiple SAR images of long time series and then extracts the surface deformation information with high accuracy. Among the many time-series InSAR methods, the small baseline subset-InSAR (SBAS-InSAR) method has the advantage of obtaining small deformation information and long time series of slow surface deformation [18, 26, 27] and is therefore the most widely used in landslide monitoring.

Based on the above background, the SBAS-InSAR method is used to obtain the deformation rates of the Wangjiashan ancient landslide in the radar line-of-sight (LOS) and slope aspect directions based on Sentinel-1A radar images and to extract historical time-series deformation information of this landslide. Furthermore, the landslide deformation characteristics and influencing factors are analyzed in conjunction with a field investigation. On this basis, the evolution trend of the Wangjiashan landslide in the impoundment period is discussed. The research in this study provides a basis for the monitoring, early warning, and risk management of the Wangjiashan landslide during the impoundment period. This work also provides a useful reference for investigations using space-borne SAR Earth observations in geological disaster prevention and control.

2. Wangjiashan Landslide

Due to increasing demands for clean energy and the unique geographical conditions of Southwest China, several hydropower stations, including the Baihetan Hydropower Station, have been built in the lower reaches of the Jinsha River in this region. Within the Baihetan Reservoir area is the Wangjiashan landslide, which is located on the right bank of the Xiaojiang River and is only 1.3 km from the Xiangbiling resettlement area (Figure 1). This valley is characterized by a subtropical dry and warm climate with an average annual temperature ranging from 12 C to 20 C. The annual precipitation in the study area is approximately 600–800 mm, with precipitation falling more frequently from May to October, and this period accounts for 90% of the annual rainfall; in contrast, from November to April of the following year, the weather is sunny and dry with sparse rainfall.

The Wangjiashan landslide is 800 m long and 90–500 m wide and covers approximately 23.5 × 104 m2 (Figure 2). The elevation of the landslide toe on the right bank of the Xiaojiang River is approximately 765 m, while the crown is at an elevation of 1125 m, resulting in an elevation difference of 360 m between the landslide crown and toe. Gullies bound the left and right sides of the landslide, and these two gullies extended upward, ultimately intersecting at the landslide crown. The 870–950 m elevation section of the landslide features flat-topped hills and landslide depressions, and the terrain slope is 15–20°, whereas from an elevation of 870 m to the toe of the landslide, the terrain slope is 35–40°, and the local slope ranges from 45–50°.

The longitudinal profile of the landslide taken from the borehole data (see Figure 2, for borehole location) of the Power China Huadong Engineering Corporation (EHEC) is shown in Figure 3. The deposits on the upper landslide are mainly gravel soil that is gray-yellow, loose, and 14.0–49.9 m thick, whereas the lower landslide deposits are silty gravel that is also gray-yellow to dark gray in color, dense, and 3.04–13.30 m thick. The soil in the sliding zone comprises gray-black to gray-brown gravel clay distributed in interlayers with silty gravel in the lower part of the landslide. The landslide bed is bedrock composed of Ordovician, Devonian, and Carboniferous limestone, dolomite, and quartz sandstone intercalated with siltstone, mudstone, and shale. The strata are oriented N0-10 W and dip NE∠40–46°, and the strata are affected by the tectonic effect; moreover, the joint network is dense, and the rock mass is broken and highly weathered. According to drilling data, the landslide deposit is 14.0–87.6 m thick, with a total volume of approximately 611 × 104 m3. The landslide has long exhibited creep deformation, and thus, the stability and disaster risk of the Wangjiashan landslide after the impoundment of the Baihetan Reservoir has aroused concern among hydropower station authorities. Now, it is urgent to understand the deformation characteristics of the Wangjiashan landslide before impoundment to reduce or eliminate the risk of geological disasters during impoundment.

3. Data and Methods

3.1. Data

Because the ancient Wangjiashan landslide tends to move toward the SW, the terrain slope is smaller than the complementary angle of the central incident angle (36.8°) of the Sentinel-1A descending track [28]. Hence, Sentinel-1A descending track data are conducive to analyzing the deformation the landslide. The coverage of the Sentinel-1A descending track data is shown in Figure 1. In this study, 77 scenes from the descending track of Sentinel-1A are used to investigate the historical deformation of the Wangjiashan landslide before impoundment from February 19, 2017, to March 30, 2021, with a time interval of 12 days. To eliminate the influence of terrain relief on the InSAR results, the Advanced Land Observing Satellite (ALOS) World 3D (AW3D30) digital elevation model (DEM) with a spatial resolution of 12.5 m is used as the external DEM data to simulate the terrain phase.

Moreover, to reduce the phase error caused by the orbit error, before processing the Sentinel-1A data, we download precise orbit data at the corresponding time and import the relevant data to correct the orbit information. In addition, monthly precipitation data and seismic data for the study area from January 2017 to March 2021 are acquired from the China Meteorological Data Service Center and the China Earthquake Networks Center, respectively. Finally, the engineering geological information of the Wangjiashan landslide is obtained from the geotechnical investigation report of the EHEC.

3.2. SBAS-InSAR Method

The SBAS-InSAR approach was first proposed by Italian scholars Berardino et al. [29]. SBAS-InSAR technology overcomes the decoherence phenomenon caused by an excessively long spatiotemporal baseline and alleviates atmospheric effects and terrain influences on the interferometric difference [29, 30]. The SBAS-InSAR technique obtains time series of surface deformation by the least-squares method while using singular value decomposition (SVD) to fit the deformation in both time and space; the InSAR surface deformation monitoring accuracy reaches the millimeter level, and thus, this technique exhibits a high measurement accuracy in landslide deformation monitoring research [31].

N + 1 scenes of SAR images are acquired in the same area and arranged in chronological order; then, one scene is chosen as the main image for registration, while the other SAR images are registered to the main image by setting appropriate temporal and spatial baseline thresholds. Ultimately, M scene interferograms are obtained. Assuming that t is the image acquisition time, after removing the flat terrain effect, terrain phase, and atmospheric phase error, the interference phase at pixel point (x, r) (x is the azimuth coordinate and r is the range coordinate) can be obtained from the phases and in the images generated at times tA and tB, respectively [29, 32, 33].

During SBAS-InSAR data processing, M interferograms are generated. is adopted to represent the matrix composed of the image deformation phases at point (x, r) at N times; then, the matrix composed of M interferometric pair phase values is , where A is an M × N matrix. When M ≥ N, the deformation time series can be obtained by the least-squares method. When M < N, the result obtained by the least-squares method is not unique, SVD is implemented to jointly solve multiple small baselines, and the phase velocities obtained for each period are integrated into the time domain to obtain the deformation time series of the entire observation period.

In this study, ENVI SARscape software is used to analyze the historical deformation of the Wangjiashan landslide. The temporal baseline threshold is set to 72 days, and the critical baseline percentage is 2%. The 20201025 SAR image is used as the main image, while images from other periods are used as auxiliary images for interferometric pairing, and 213 sets of interferometric pairs are generated. The relationship between the specific acquisition dates of the descending track data and the spatiotemporal baseline is shown in Figure 4. Synthetic interferogram inversion, geocoding, and raster vector conversion are performed to generate the mean LOS deformation rate, and finally, the time-series deformation results are obtained.

3.3. Slope Aspect Velocity Fields

The sliding direction of a landslide is generally downward along the slope. In the process of radar interferometry, the obtained deformation is usually the deformation value along the radar LOS direction, which often cannot accurately reflect the true deformation of the slope surface [34, 35]. According to the geometric relationship among the radar LOS direction, the slope aspect, and the vertical settlement direction, the deformation in the radar LOS direction is transformed into slope aspect deformation [36, 37]. The specific conversion formula is as follows:

In the above equations, θ represents the incident angle of the radar, α is the satellite flight direction angle, β is the slope direction, and φ is the slope. Among them, β and φ can be calculated from the DEM.

(1) demonstrates that if the value of is close to 0, will tend toward infinity, but the ground unit at this time is completely decoherent. To correct this anomaly, according to the study of Herrera et al. [38], can be set as a threshold so that the acceptable projected velocity value can be distinguished from the data with unreliable projection values. That is, when , , and when , . In addition, when , the landslide mass moves upward along the slope. Even if expansion causes a bulge on the front edge of the landslide or the landslide mass accumulates at the foot of the slope, the displacement vector in the horizontal direction remains oriented in the downward direction of the slope surface. These deformation points with need to be eliminated consistently. In this way, the average deformation rate of the surface displacement along the slope direction in the study area can be obtained (Figure 5).

4. Results

4.1. Mean Deformation Velocities

In this study, the SBAS-InSAR method is utilized to analyze the mean deformation rate of the Wangjiashan landslide in the radar LOS and slope aspect directions. As shown in Figure 6, the mean deformation rate of the Wangjiashan landslide in the LOS direction ranges from −69 to −14 mm/year, while the mean deformation rate in the slope aspect direction reaches −220 mm/year. The landslide deformation is concentrated in the white dashed circle in the middle with a longitudinal length of approximately 600-700 m. The mean VLOS ranges from −41 to −69 mm/year, and the mean VSlope ranges from −120 to −220 mm/year. The field investigation revealed deformation phenomena such as fracturing (Figure 7(a)) and toppling of the retaining wall (Figure 7(b)) in the area of the white dashed circle. Furthermore, as many as 12 cracks are found in the slope (Figure 2); the longest crack is approximately 17 m long, and the tensile crack with the maximum width is 30 cm wide (Figure 7(c)). These cracks are mostly located around the road and are dominated by tension cracks, which may be the result of landslide deformation caused by road excavation. Some deformation phenomena, such as local collapse in front of the landslide, are also observed (Figure 7(d)).

4.2. Time-Series Cumulative Deformation

Based on SBAS-InSAR technology, nine time-series cumulative deformation graphs of the Wangjiashan landslide are obtained (Figure 8). It should be noted that all accumulated deformation is referenced to the first image on February 19, 2017. In 2018, the deformation area of the landslide was in the middle of the landslide area (red part). In 2019, the landslide deformation expanded to the front and rear edges. Over time, the deformation of the Wangjiashan landslide gradually increased, and the maximum cumulative deformation in the LOS direction reached 294 mm. As shown in Figure 9, the deformation of the Wangjiashan landslide exhibits obvious spatial heterogeneity, with the largest deformation in the middle of the landslide, followed by the landslide toe, and the smallest deformation at the landslide crown. According to the field investigation, the cracks are mainly developed in the middle of the landslide (Figure 2 and Table 1), consistent with the area with the largest cumulative deformation.

5. Interpretation of Landslide Motion

5.1. Kinematic Behavior

To further understand and analyze the temporal and spatial characteristics of landslide deformation, the mean LOS deformation rate curve along section B–B′ in Figure 6(a) and the LOS cumulative deformation curve along section B–B′ in Figure 9 are superimposed onto the geological profile, and the result is shown in Figure 5. A comparison of the data revealed the following:(i)The largest deformation of the Wangjiashan landslide is in zones B and C. The average LOS cumulative deformation is 220 mm, and the maximum LOS cumulative deformation is 289 mm. The deformation of area A is relatively low; the average LOS cumulative deformation in zone A is 90 mm, and the maximum LOS cumulative deformation is 120 mm. Likewise, the deformation in area E is very weak; the average LOS deformation rate is only eight mm/year. The cumulative deformation and LOS deformation rate in zone D increase along the profile. The deformation rate at the crown of the Wangjiashan landslide is the smallest, while the deformation rates of the central zone (zone C) and the toe (zone B) are larger. As shown in Figure 5, the middle of the landslide deformed first and then expanded to the front and rear edges. The landslide exhibits the deformation characteristics of a retrogressive landslide [39].(ii)The data of zone F show a decreasing deformation rate but increasing cumulative deformation. A field investigation revealed that the retaining wall on the side of the road in this area is broken (Figure 7(b)), and the road surface was bulging (Figure 10) and cracked (Figure 7(c)). These observations indicate that the sliding force of the landslide mass exceeded the slip resistance of the retaining wall, resulting in damage to both the retaining wall and the ground. Considering the deformation rate in conjunction with the topographic characteristics [40, 41], due to excavation of the road, the upper landslide seems to have formed a local sliding surface, leading to serious road deformation. Of course, the final determination of this local sliding surface requires further detailed research.

5.2. Influencing Factors and Evolution Trend

Gravity-driven retrogressive sliding is the main deformation mechanism of the Wangjiashan landslide. However, other factors, such as changes in the groundwater level, heavy rainfall, and earthquakes, can accelerate deformation [42]. Nevertheless, no groundwater was discovered in the boreholes in the Wangjiashan landslide, which indicates that the groundwater level within the landslide is much lower than the sliding surface (Figure 3). Therefore, groundwater is not the main factor affecting landslide deformation. To understand the response relationship between deformation and rainfall, four monitoring points (P1, P2, P3, and P4 in Figure 6(b)) with different elevations are selected from the slope deformation rate map to further analyze the correlation between the slope time-series deformation and the precipitation observed by adjacent meteorological stations.

Figure 11 shows the slope aspect time-series cumulative deformation and deformation rate of the Sentinel-1A data from February 2017 to March 2021 at the four observation points P1, P2, P3, and P4. The deformations at P1, P2, P3, and P4 increase linearly, and the deformations increase obviously during the rainy season. In the rainy season of 2017, the deformations at the four observation points increased. However, in the 2018 rainy season, only the deformations at P2, P3, and P4 changed significantly with seasonal precipitation. In July 2019, the monthly precipitation reached 334.5 mm, and the cumulative displacements at the four observation points increased significantly. Through the deformation rate curve of the landslide time series, we can see that the displacement rate of landslides in the rainy season is increased compared with that in the dry season. This shows that heavy rainfall will accelerate the deformation of a landslide. Because the Wangjiashan landslide is a reactivated ancient landslide, there is also deformation in the dry season, but the deformation rate is less than that in the rainy season. These findings indicate that rainfall is one of the main reasons for the reactivation and deformation of the Wangjiashan landslide.

In contrast, during the rainy season of 2020, the deformations at the four observation points did not increase significantly because the rainfall was low. In addition, the cumulative rainfall from June to August in the study area during 2019 and 2018 was 504.5 mm and 584.5 mm, respectively, and the cumulative deformation at point P2 was 61.86 mm and 50.38 mm, respectively. The above phenomenon indicates that concentrated heavy rainfall is more likely to cause large deformation of the landslide [43, 44]. The above data analysis shows that the deformation of the Wangjiashan landslide in the rainy season is related to the accumulation of rainfall and concentrated heavy rainfall.

In addition, seismic data within a 100 km radius of the Wangjiashan landslide from 2017 to 2020 were downloaded from the China Earthquake Networks Center (http://www.ceic.ac.cn/) (Figure 2). Based on its epicentral location, focal depth, magnitude, and intensity, the Ms5.0 earthquake at 21 : 47 on May 18, 2020, caused many houses to collapse and damaged several roads in Qiaojia County. The epicenter of the earthquake was 74 km from the Wangjiashan landslide, and the focal depth was 8 km. However, comparative data show that the earthquake did not accelerate the deformation of the Wangjiashan landslide (Figure 11).

6. Discussion

According to the theory of creep deformation [45], before a landslide occurs, it usually undergoes an initial deformation stage, a constant deformation stage, and an accelerated deformation stage (which in turn is divided into three substages: an initial acceleration stage, a moderate acceleration stage, and an imminent sliding stage). As shown in Figure 12, Xu et al. [47, 48] proposed using the tangential angle (α) as the basis for dividing these three stages. When the tangent angle is much less than 45°, the landslide is considered to be in its initial deformation stage; when the tangent angle is approximately equal to 45°, the landslide is assumed to undergo constant velocity deformation; when the tangent angle greatly exceeds 45°, the landslide is believed to be in its accelerated deformation stage (with 80° and 85° representing the boundaries between the first and second substages and the second and third substages, respectively). The gravity landslide warning model shown in Figure 12 has been widely used to monitor and warn landslides in Northwest and Southwest China [46, 4953].

The long-term displacement curve in Figure 11 shows that the Wangjiashan landslide has exhibited a stable deformation rate in the past four years, except for a briefly accelerated rate during the rainfall period and average tangential angle α between 23° and 32°, which is less than 45° (Figure 13). Therefore, from the perspective of its evolution, the Wangjiashan landslide was in the initial deformation stage before impoundment and maintained slow deformation without strong interference from external factors.

The signs of deformation during the impoundment period of the Wangjiashan landslide and studies on areas such as the Three Gorges Reservoir have shown that fluctuations in the water level of large reservoirs may lead to the reactivation of ancient landslides, resulting in considerable losses [3, 4, 54, 55]. Given the rainfall concentration characteristics of the study area and the sensitivity of the deformation of the Wangjiashan landslide to rainfall, its future development deserves further attention. Therefore, it is urgent to systematically identify and analyze potential landslides in the Baihetan Reservoir area and establish a monitoring and early warning system to reduce the losses of life and property.

7. Conclusions

Sentinel-1A SAR data were used to study the time-series deformation characteristics of the Wangjiashan landslide in the Baihetan Reservoir area before its impoundment. We found that the Wangjiashan landslide has been in an initial deformation state in recent years. Deformation first occurred in the middle of the landslide and then expanded to the landslide toe and crown, thereby demonstrating the characteristics of retrogressive movement. Combined with an analysis of field deformation evidence, the upper landslide mass seems to have formed a local sliding surface, which has caused serious deformation of the road. An analysis of historical rainfall data revealed that the Wangjiashan landslide is sensitive to rainfall, and the deformation is not only significantly correlated with cumulative rainfall but also influenced by concentrated heavy rainfall. Considering the sensitivity of the Wangjiashan landslide to rainfall, coupled with the environmental changes caused by impoundment, establishing a landslide monitoring and early warning system in this area is particularly important.

Data Availability

The precipitation data and historical earthquake data were provided by China Meteorological Data Service Center (http://data.cma.cn/) and China Seismic Network Center (http://www.ceic.ac.cn/), respectively. Sentinel-1 A/B images were from the European Space Agency (ESA) project, and these data were downloaded from the Sentinel-1 Scientific Data Hub. DEM data were provided by the Japan Aerospace Exploration Agency (JAXA).

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This research was funded by the National Key R&D Plan (Grant no. 2018YFC1505503), the State Key Laboratory of Geohazard Prevention and Geoenvironment Protection Independent Research Project (Grant no. SKLGP2020Z001), and the Zhejiang Huadong Construction Engineering Co., Ltd Research Project (KY2020-HDJS-19).