Next Article in Journal
Grassland Conservation Effectiveness of National Nature Reserves in Northern China
Next Article in Special Issue
Fourfold Bounce Scattering-Based Reconstruction of Building Backs Using Airborne Array TomoSAR Point Clouds
Previous Article in Journal
Combinations of Feature Selection and Machine Learning Algorithms for Object-Oriented Betel Palms and Mango Plantations Classification Based on Gaofen-2 Imagery
Previous Article in Special Issue
A Comparative Study on Classification Features between High-Resolution and Polarimetric SAR Images through Unsupervised Classification Methods
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

InSAR Study of Landslides: Early Detection, Three-Dimensional, and Long-Term Surface Displacement Estimation—A Case of Xiaojiang River Basin, China

1
Space Microwave Remote Sensing System Department, Aerospace Information Research Institute, Chinese Academy of Sciences, Beijing 100190, China
2
School of Electronics, Electrical and Communication Engineering, University of Chinese Academy of Sciences, Beijing 100049, China
3
China Aero Geophysical Survey & Remote Sensing Center for Natural and Resources (AGRS), Beijing 100083, China
4
National Key Laboratory of Microwave Imaging Technology, Aerospace Information Research Institute, Chinese Academy of Sciences, Beijing 100190, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2022, 14(7), 1759; https://doi.org/10.3390/rs14071759
Submission received: 2 March 2022 / Revised: 30 March 2022 / Accepted: 1 April 2022 / Published: 6 April 2022
(This article belongs to the Special Issue Recent Progress and Applications on Multi-Dimensional SAR)

Abstract

:
Landslides, a major natural geohazard, obstruct municipal constructions and may destroy villages and towns, at worst causing significant casualties and economic losses. Interferometric Synthetic Aperture Radar (InSAR) technique offers distinct advantages on landslide detection and monitoring. In this paper, a more systematic workflow is designed for InSAR study of landslides, in terms of three levels: (i) early detection on regional scale, (ii) three-dimensional (3D) surface displacement rates estimation on detailed scale, and (iii) time series analysis on long-term temporal scale. The proposed workflow is applied for landslide research on the Xiaojiang River Basin, China, using ascending and descending Sentinel-1 images acquired from March 2017 to May 2019. First, the landslide inventory has been mapped and updated using InSAR stacking method, supporting geohazard prevention on a regional scale. A total of 22 active landslides are identified, ranging from medium to super large scale. Compared with the existing inventory, three unrecorded landslides are newly detected by our approach, and five recorded landslides are detected significant expansion of their boundaries. Then, specific to a detected landslide, Baobao landslide, a Total Least Squares–Kalman Filter-based approach is presented. Two outcomes are provided for further spatial-temporal pattern analysis: 3D displacement rates, providing an intuitive insight on the spatial characteristics and sliding direction of landslide, which are analyzed to deep the understanding of its kinematic mechanism, and long-term time series, which contribute to deduce the dynamic evolution of landslide, presenting benefits in landslide forecasting.

1. Introduction

As a major natural geological hazard, landslides interfere with the transport net-works and infrastructure development or may hold back the expansion of urbanization [1,2,3], and at worst, destroy villages and towns, causing significant casualties and eco-nomic losses [4,5,6,7]. There is a need for increased research on landslides, especially in terms of three levels, i.e., early detection on a regional scale, estimation of three-dimensional displacement rates on a detailed scale, and retrieval of displacement time series on a long-term temporal scale.
The detection of unstable slopes is an essential prerequisite for early warning and geohazard prevention [8], because of their extensive distribution and concealment. Interferometric Synthetic Aperture Radar (InSAR) performs quantitative remote sensing for small displacements over wide areas, with the advantages of centi- or even millimeter precision, high spatial-temporal resolutions, and all-day, all-weather working capabilities [9]. Different InSAR techniques [10,11,12] have demonstrated their distinctive capabilities in regional landslide detection, achieving valuable examples in Italy nationwide [13,14], the Alps [15,16], California [17], and the Jinsha River Catchment, China [18,19]. For instance, Novellino et al. [20] utilized the intermittent Small Baseline Subset (SBAS) to analyze the slope stability in northwestern Sicily, Italy; Bouali et al. [21] combined the Permanent Scatterers Interferometry (PSI) results with the landslide inventory to detect extremely slow (<16 mm/year) landslides across the Palos Verdes Peninsula, California; Liu et al. [19] employed InSAR stacking on Scan Synthetic Aperture Radar (SAR) images, achieving large-scale geohazard identification along the Jinsha River Corridor. Furthermore, the performance of different InSAR techniques for landslide detection has been compared in detail [22], revealing that the stacking method performs better on poor-coherence areas with the advantages of technical requirements and low computation labor. However, conventional InSAR techniques measure only one component of the surface displacement, i.e., in the line-of-sight (LOS) direction of satellite. This poses a challenge to interpretation and detection of landslides, and more than that, narrows the understanding of their displacement patterns and kinematic mechanisms.
Three-dimensional (3D) displacement rates provide intuitive insights on the spatial characteristics and sliding direction of landslides, which can be contributed to kinematic mechanism analysis [23]. To date, various strategies have been developed to estimate 3D displacement rates based on InSAR observations. The most straightforward strategy is to integrate multi-track InSAR observations acquired with different imaging geometries [24,25]. However, since most SAR satellite orbits are near-polar, this strategy may be applied solely for some specific cases in high latitude regions with higher azimuth angle diversity. Another strategy is to combine InSAR measurements with outputs from either offset tracking [26,27], Multiple Aperture Interferometry (MAI) [28,29], or other ancillary sources, e.g., Global Navigation Satellite System (GNSS) [30,31,32]. Moreover, a strategy is to impose a priori information as constraints on deformation process studied, which is typically applied for mining subsidence [33,34,35] and glacier or landslide motions [36,37,38,39,40,41]. Specific to landslide studies, a classical constraint is the surface-parallel flow model, i.e., movement along the steepest downslope, which has been accepted in many representative landslide cases, e.g., the Slumgullion Landslide, USA [38] and the Jiaju Landslide, China [39]. However, this constraint seems too strict for most landslides except translational landslides [40], since it does not allow consideration of the difference between ground surface and slip surface, forcible ignoring potential errors in model design. Overall, it is worth noting that the above strategies are typically developed to measure single deformation episodes with large gradient, e.g., earthquakes, fasting landslides, and volcanic activities. There are few studies on 3D displacement estimation of slow-moving landslides [40,41], particularly on a long-term temporal scale.
In contrast, as the entire landslide process can last from months to years, long-term time series of landslide movement is of great significance to deducing the dynamic evolution, presenting particular benefits in landslide forecasting and risk management [42,43,44]. Some efforts have been devoted for extending aforementioned strategies on temporal scale, e.g., the Pixel-Offset SBAS technique [45], the Multidimensional SBAS (MSBAS) technique [40,46], the Minimum Acceleration (MinA) approach [47], and the approach based on the Kalman Filter (KF) [48]. Among them, MSBAS and MinA are developed for simultaneous post processing of hundreds to thousands of individual interferograms. They regularize the rank-deficient and ill-posed inversion system by imposing further condition (e.g., Tikhonov regularization [40] or the assumption that displacement time series have minimum acceleration [47]), and then solve the 3D displacement time series utilizing the Singular Value Decomposition (SVD) method. In contrast, the KF-based approach is a sequential processing, which combines multi-source InSAR observations at the present epoch in a Kalman filter model with displacement information determined from previous interferograms. This approach allows the displacement time series dynamically update with an improved temporal resolution and limited computational cost. Its performance has been demonstrated in researches on volcanics and faults motions [48,49], but up to now, this approach is rarely applied to the retrieval of the time series of landslide displacements.
In this study, a more systematic workflow is designed for InSAR-based landslide study, including three levels, i.e., (i) early detection on regional scale, (ii) 3D displacement rates estimation on detailed scale, and (iii) time series analysis on long-term temporal scale. First, a standard InSAR stacking procedure is employed to investigate the distribution and activities of landslides over study area. A relatively complete landslide inventory is mapped and updated for geohazard prevention on regional scale. Then, specific to a detected landslide, 3D displacement rates and long-term time series are further estimated for the spatial-temporal pattern analysis in more detail. To this end, a Total Least Squares–Kalman Filter (TLS–KF)-based approach is presented, combining advantages of Total Least Squares (TLS) and advantages of the Kalman Filter (KF). The outcomes can provide intuitive description of sliding direction and displacement characteristics of landslides, deepening the understanding of its kinematic mechanism and dynamic evolution.
Here, the proposed workflow is applied for landslide research on the Xiaojiang River Basin, Eastern Yunnan, China, which is the most serious area of soil erosion and geological disaster in the upper reaches of Yangtze River [50]. This paper is organized as follows. A brief description of the study area, as well as the datasets, is given in Section 2. The proposed workflow is presented in Section 3, with its theory and processing elaborated in detail. In Section 4, a new landslide inventory has been mapped and updated by our approach, with additional information of LOS displacement rates. Moreover, one of the detected super large landslides, the Baobao landslide, is further studied in Section 5. The three-dimensional and long-term displacements are estimated by the TLS–KF-based approach and its spatial-temporal displacement pattern are discussed.

2. Study area and Datasets

2.1. Study Area

The study area is situated in the Xiaojiang River Basin, Eastern Yunnan, China (Figure 1), north from Menggu and south of and nearby Dongchuan, covering an area of approximately 2314 km2. It is the most serious area of soil erosion and geological disaster in the upper reaches of the Yangtze River [50], as well known by the Natural Museum of Landslides [51], under combined actions of multiple factors, e.g., complex natural geography, lithology, active tectonic movements, special climate, and human activities. The elevation in most parts of study area is higher than 2500 ma.s.l., reaching more than 4000 ma.s.l. in some regions. Due to violent river downward eroding, valleys feature strong V-shaped topography with height differences up to 3000 m, resulting in slope angles of greater than 25° in most slopes, thus providing favorable geological conditions for the occurrences of landslides. The tectonic setting is conditioned by a series of nearly north-south trending faults, including the Xiaojiang Faults, the Pudu River Faults, and the Yuanmou Faults (marked as F1, F2, F3 in Figure 1), resulting in frequent seismic activities. There have been approximately 41 earthquakes of Mw ≥ 3.0 in the study area and its surroundings in the past ten years (2009–2019), including three stronger earthquakes greater than Mw = 5.0. Generally, rock masses in active tectonic zones are more prone to develop joints and fissures and to fracture [52], presenting poor stabilities. Therefore, landslide susceptibility tightly correlates with the distribution and activities of faults. Specific to the study area, tectonic movements are mainly affected by the north section of the Xiaojiang Faults, which is the most active sinistral strike-slip fault zone in southeast margin of the Qinghai-Tibet Plateau [53].
Figure 2a shows the geological setting map of the study area. The outcrops are composed of almost all strata from Jurassic to Mesoprotero-zoic, in detail, including middle-lower Jurassic, Permian, Triassic, Carboniferous-Devonian, Silurian, Ordovician, Cambrian, Paleozoic, Sinian, and Kunyang Group strata. They mainly include basic volcanic rocks, dolomites, limestones, phyllites, mudstones, and sandstones, and diabase in local regions. The lithology of strata features as soft and hard rocks alternatively distributed, leading to poor compactness and low shear strength [54], and it is an inherent factor in the full development of landslides.
The climate belongs to subtropical plateau monsoon climate [55]. The annual temperature and precipitation are approximately 15 °C and 800 mm, respectively. As a typical hot-dry valley, the climate features relatively high temperature and low humidity and precipitation concentration (88% of precipitation is concentrated during June to October). The former leads to strong physical weathering on the slope surfaces; the latter leads to rain infiltration decreases shear strength of materials, thus driving landslide movements. For instance, Jiangjia Gully debris flow (in Figure 2a, mapped with white curve) is induced by strongly weathered limestones under concentrated heavy rainfalls. In addition, the study area is rich in mineral resources, among which Dongchuan is one of six famous copper productions [56]. Mining activities lasting more than 2000 years have broken the balance of groundwater and destroyed the initial stress states of some slopes, promoting the occurrence of landslides.
As noted, multiple factors work together to result in the extensive distribution and strong activity of large-scale landslides in the Xiaojiang River Basin. Potential landslides threaten the safety of cities and towns as well as traffic lines (e.g., Provincial Highway S23). Moreover, the landslides could dump large amounts of rock/soil mass into the Xiaojiang River, largely endangering the operation of Baihetan giant hydropower station on downstream of its mainstream. China Aero Geophysical Survey & Remote Sensing Center for Land and Resources (AGRS) has mapped the distribution of landslides in this area through optical remote sensing (see Figure 2) and verified them by field investigation. However, the detection only based on optical images may omit some active landslides, especially those without significant morphological features (e.g., scarps, sliding masses, and bulging toes). More importantly, complete investigations of landslides in terms of 3D displacements, spatial pattern, and kinematic evolution are absent [41].

2.2. Datasets

The detection of active landslides is an essential task for disaster prevention. More than that, the estimation of 3D landslide displacement rates and long-term time series is of great significant to deepen the understanding of its kinematic mechanism and dynamic evolution. To these ends, a total of 117 SAR images are collected from Sentinel-1 ascending Track 26 and descending Track 62, covering the period from March 2017 to May 2019. The spatial coverage of the SAR images provided by each track is illustrated in Figure 1. Other basic parameters of the SAR images are summarized in Table 1, in detail. Moreover, a 1-arc-second SRTM digital elevation model (DEM) is employed in the standard differential InSAR processing, and in the calculation of the LOS visibility maps of ascending and descending tracks [57]. The holes (non-value pixels) in the DEM are filled by bicubic spline interpolation.

3. Methodology

The InSAR stacking method is employed to detect active landslides in the study area. A new approach combining Total Least Squares (TLS) and Kalman Filter (KF) is presented to estimate 3D displacement rates and long-term time series of detected landslides. The workflow is shown in Figure 3, which consists of three steps, as follows.
Step 1. Generation of LOS displacement rate maps and detection of active landslides. First, each SAR dataset is handled with a standard differential InSAR procedure, mainly including co-registration, interferograms generation, phase unwrapping, errors correction, and quality screening. Then, high-quality unwrapped interferograms are outputted for the following processing, stacking interferograms [12]. Relative InSAR processing is elaborated in Section 3.1. As outputs, the LOS displacement rate maps are derived from ascending and descending datasets, with a spatial resolution of 15 m and in a Universal Transverse Mercator (UTM) geodetic system. These rate maps are superimposed onto optical images, combined with morphological features (e.g., scarps, sliding masses, and bulging toes), to interpret and identify active landslides. Finally, suspected landslides detected from each dataset are mosaiced [21] to map a relatively complete landslide inventory.
Step 2. Estimation of 2D/3D displacement rates and analysis of spatial pattern. Specific to a detected landslide, if the north-south displacement component can be neglected, the 2D displacement rates are derived by fusing ascending and descending unwrapping interferograms [40]. Otherwise, the 3D displacement rates are estimated in total least squares sense by imposing a relaxed constraint, i.e., movement along slip surface, which is more general than the ground-surface parallel flow constraint [40,58]. The methodology is presented in Section 3.2. Based on the outputs, the optimal sliding direction of landslide is determined for the analysis of the spatial displacement pattern.
Step 3. Retrieval of long-term displacement times series and temporal evolution. Ascending and descending unwrapped interferograms are fused to retrieve the long-term displacement time series in the optimal sliding direction based on the Kalman Filter [48]. The methodology is described in Section 3.3. Then, the temporal evolution of landslide is investigated based on the OSD displacement times series and the possible driving factor, precipitation.

3.1. InSAR Processing

A standard differential InSAR procedure is employed to process all SAR images, as follows. First, single look complex (SLC) images from ascending and descending Sentinel-1 are co-registered, separately, with the aid of auxiliary data (e.g., precise orbit files and SRTM DEM). Second, the small baseline subset (SBAS) combination strategy [10] is adopted to generate high-quality interferograms with less temporal and spatial decoherence. Specific to our datasets, the temporal baseline and spatial baseline are restricted to 40 days and 200 m, respectively. In processing, the multi-look factor is set to 4 (4 pixels in range and 1 pixel in azimuth), resulting in interferograms characterized by a ground resolution of approximately 15 m × 15 m. Similar resolution data have proven their performance in regional landslide mapping in many previous works [59,60]. For a given interferogram, the minimum cost flow algorithm [61] is applied to unwrap phases. A semi-adaptive procedure is performed to correct various items of phase errors. The unwrapping error is corrected by either a phase compensation algorithm or re-unwrapping at a high coherence point [19]. The error related to residual baseline is estimated and reduced based on a quadratic polynomial model [62]. The error caused by DEM is corrected using the method presented in [63]. The atmospheric phase is corrected by the spatially variable power law tropospheric correction technique [64], preliminarily. It is worth noting that the climate in the study area features complex and ever-changing local micro-climate, so the atmospheric phase needs further correction in the follow-up. Finally, a total number of 186 and 131 high-quality unwrapped interferograms are screened from ascending and descending Sentinel-1 datasets, respectively. Their spatiotemporal distributions are illustrated in Figure 4.
Then, the stacking interferograms method [12] is employed to retrieve the displacement rate of each dataset in LOS direction. Equation (1) shows its kernel, a linear combination of interferograms ( φ i ) with the temporal baseline ( Δ T i ) as weights:
V L O S = λ 4 π i = 1 N Δ T i φ i i = 1 N Δ T i 2  
While the atmospheric phase is random for independent interferograms, the related error was effectively corrected by stacking, and the signal-to-noise ratio was improved by a factor of N [65]. Specific to our datasets, compared with the descending dataset, the ascending dataset generates a larger number of interferograms (approximately 1.5 times) for stacking, thus inferring that atmospheric errors will be corrected more effectively and a higher accuracy will be achieved. In addition, the LOS displacement rate maps are geocoded into UTM grid system (e.g., World Geodetic System 1984), preparing for the integration of multi-source InSAR results to detect active landslides. Moreover, invalid measurement scatterers in flat areas and in areas presenting low degree of LOS visibility (e.g., layover and shadow) are masked to avoid possible misjudgments.

3.2. Total Least Squares Estimation of Three-Dimensional Landslide Displacement Rates

If only two independent SAR datasets are available and the north-south displacement component is non-negligible, a prior model should be imposed as the constraint on landslide movement to reduce the freedom degree for inversion of 3D displacement rates. A classical constraint is the surface-parallel flow model [58], which has proved its applicability on translational landslides. However, this constraint is too restrictive for most landslides, as their mechanisms and types are more complex in reality, and their buried slip surfaces are usually undeterminable to be parallel to topography. To generalize further, the constraint is relaxed for non-translational landslides, as follows:
V U = H + Δ h x V E + H + Δ h y V N  
where V U , V E , and V N are the displacement rate parameters in vertical, east-west, and north-south direction, respectively, H + Δ h x and H + Δ h y represent the first derivatives of slip surface in the east and north directions, respectively, H is the elevation of the ground surface, which was derived from a smoothed external DEM [66], and Δ h is the elevation difference between the ground surface and slip surface.
Under the relaxed constraint, i.e., movement along slip surface, the observation model is constructed as follows:
S E G S N G S U G H + Δ h x H + Δ h y 1 V E V N V U = L 0
where L is the observation vector in LOS direction derived from ascending and descending tracks, S E , S N , and S U are the east-west, north-south, and vertical components of the LOS unit vector S [40], and G represents the time intervals between consecutive SAR acquisitions. It is worth noting that the slip surface is roughly deduced using the strain model [67] and/or field exploration [68], and is sometimes even inaccessible; thus, Δ h is unknown with considerable error. Furthermore, whereas Δ h is generally small enough relative to H , the terms composed of Δ h are regarded as Errors-In-Variables (EIV) and the observation model is converted into the EIV frame, as follows:
L = G ^ + E G X + V  
where X = V E , V N T is the unknown parameter vector, G ^ = G S E + H x S U ,   S N + H y S U is the updated design matrix, and E G = S U G Δ h x , Δ h y and V represent errors of model design and observation vector, respectively, with stochastic properties characterized by:
V v e c E G ~ 0 0 , Σ V 0 0 Σ G  
where v e c indicates the operator that stacks one column of a matrix under the previous column and Σ V and Σ G are the variance-covariance matrixes of observation error and design error, respectively.
It is well known that errors of model design are ignored in the classical Least Squares (LS) estimation; thus, the solutions in LS sense are no longer optimal. To adjust to the EIV observation model, the Total Least Squares (TLS) estimation takes into account both errors in observations and model design, resolving based on the criterion of minimizing the objective function with Equation (6). Their performances have been proved in some similar problems, e.g., 3D coordinated transformation [69], in GNSS:
m i n V T Σ V 1 V + v e c E G T Σ G 1 v e c E G  
Furthermore, Σ G can be approximately expressed in form of the Kronecker-Zehfuss product ( ), as Σ G : = I 2 Σ 0 , where I 2 is the identity matrix and Σ 0 is the variance matrix of each column of E G , composed of variance of elevation difference ( σ Δ h 2 ), as Σ 0 σ Δ h 2   S U 2 d i a g G 1 , 1 2 G N , 1 2 . Thus, a fast convergent iterative algorithm presented in [70] is adopted to obtain the optimal solution in the TLS sense (see Table 2).

3.3. Kalman Filter for Long-Term Landslide Displacement Time Series

The retrieval of long-term time series of landslide displacement is particularly meaningful to the deduction of landslide dynamic evolution. Toward this end, the Kalman filter, an effective approach to fuse multi-source SAR observations obtained at different times [48], is employed to estimate the landslide displacement time series for each pixel in the optimal sliding direction (OSD). The OSD is determined in Step 2 using the TLS estimations of 3D displacement rates, V T L S , U :   V T L S , E : V T L S , N = a : b : c . The observation model and the state model of the Kalman filter are constructed for a sliding pixel at acquisition date t k , as follows:
L k = H X k + V k  
X k = A X k 1 + W k 1
where L k is the LOS observation, X k = d k v k T is the unknown state vector, d k and v k are the OSD displacement and velocity, respectively, H = a S U + b S E + c S N a 2 + b 2 + C 2 0 is the projection matrix, A = 1 Δ T 0 1 is the state transition matrix, Δ T is the time interval between t k 1 and t k , and V k and W k 1 are the vectors of measurement noise and process noise, respectively, with stochastic properties characterized by:
V k ~ N 0 ,   R k               W k ~ N 0 , Q k
where R k is the observation variance of the kth interferogram and Q k is the variance-covariance matrix of process noise, which can be recursed using Sage-Husa adaptive filter [71], as Equation (10). The OSD displacement time series can be dynamically estimated by the Kalman filter in near-real time. All relevant equations are listed in Table 3.
Sage-Husa adaptive filter of Q k :
Q k = 1 f k Q k 1 + f k K k L k H X k L k H X k T K k T + P k A P k 1 A T
where x ˜ k and x ˜ k are the prior and posteriori state estimations, with P k and P k being the corresponding variance-covariance matrixes, K k is the Kalman gain, f k = 1 b 1 b k is the forgetting factor, and b is a preset parameter, recommended to be set at 0.95–0.99.

4. Results: Line-of-Sight Displacement Rates between March 2007 to May 2019

Following Step 1, the LOS displacement rate maps are derived from ascending and descending Sentinel-1 SAR datasets, respectively, as shown in Figure 5. It is worth noting that the negative values (red color) indicate the movements away from the satellite, and the positive values (blue color) indicate the movements toward the satellite.
Measurement scatterers (MSs) with a total number of 9,584,080 and 9,179,644 are identified from the ascending dataset (Figure 5a) and descending dataset (Figure 5b), respectively, thus producing an overall spatial density of approximately 4000 MSs/km2. These scatterers are selected by a fixed standard (a coherence threshold of 0.2 and an intensity threshold of −1 dB) and generally distributed on the roads, bridges, artificialities, and rocks/soils with sparse vegetations. Compared with the ascending track, the heading angle of descending track is more likely to cause severe geometric distortions on eastward slopes, resulting in a relatively sparse MSs for landslide detection in the study area. Moreover, the distribution histograms of MSs in stable regions are calculated for accuracy assessment, as shown in Figure 6. The mean values (red line) are lower than 0.1 cm/year, and the standard deviations (blue line) are approximately 1 cm/year.
Refocusing on Figure 5a,b, the locations of active displacement regions are roughly consistent between ascending and descending results, while the displacement magnitude and their detailed extents are locally different. This is mainly attributed to their different LOS sensitivities to mass movements in each slope orientation. To complement each other, the ascending and descending LOS displacement rate maps are combined to map a relatively complete landslide inventory. Finally, a total of 22 active landslides, ranging from medium scale to super large scale, are identified, mapped, and labeled in Figure 5, and their detailed information are summarized in Table 4. These landslides are concentrated near the Xiaojiang River and some of them are re-active parts of historical debris flows. Compared with the recorded landslides (black curves) in AGRS inventory, three unrecorded landslides have been newly detected by our approach (No. 20–22, mapped within red curves): enlarged views of them are shown in Figure 5e–g. Furthermore, five recorded landslides have been detected with a significant expansion of their boundaries (No. 14, 16–19, updated within red dotted curves): enlarged views of them are shown in Figure 5c,d.
At this point, the landslide inventory of the study area has been mapped and updated by our approach with addition information of the LOS displacement rate, supporting the regional disaster prevention and risk control. However, for the detected landslides, their 3D displacement rates and long-term time series remain to be further studied, in order to deepen the understanding of their kinematic mechanism and dynamic evolution. In Section 5, the Baobao landslide (labeled No.14 in Table 4) is taken as an example case to retrieve its 3D displacement rates and long-term time series. Furthermore, the spatial-temporal displacement pattern of Baobao landslide will be discussed in the following sections.

5. Discussion: Three-Dimensional and Long-Term Displacements of Baobao Landslide

5.1. Three-Dimensional Displacement Rates for Spatial Displacement Pattern Analysis

The Baobao landslide is one of detected landslides on super large scale, the extent of which is ~0.76 km in width and ~1.57 km in length, located on the east bank of Pucha River and in the Xiaojiang Fault zone. Its optical remote sensing image acquired in April 2019 is shown in Figure 7a, where polygons with different colors represent four Blocks B1–B4 of the landslide, which are divided according to the geomorphological analysis. Following Step 2, three-dimensional displacement rates are estimated in the vertical, north-south, and east-west directions during March 2017 to May 2019, as shown in Figure 7b–d, respectively.
It is worth noting that negative values (red color) indicate downward movements in Figure 7b, southward movements in Figure 7c, and westward movements in Figure 7d, and positive values (blue color) indicate northward movements in Figure 7c. The maximum displacement rates in the vertical, south-north, and east-west directions are more than −6.6 cm/year, ±11.0 cm/year, and −16.2 cm/year, respectively. The two most active movements are observed at the middle of block B1, i.e., Region R1, and lower-left edge of block B2, i.e., Region R2, with averaging displacement rates of 13.5 cm/year and 14.9 cm/year in the optimal sliding direction (OSD), respectively. Overall, the entire landslide presents a trend of westward and downward movement. As evidence in Figure 7c, northward displacements mainly occurred in block B1, and are sliding toward the northwest direction, while southward displacements mainly occurred in Region R2 and the lower half of block B3, and are sliding toward the southwest direction. A distinct sliding boundary can be clearly observed between block B1 and blocks B2–B3.
For spatial displacement pattern analysis, displacement rates and elevation are extracted along two representative Profiles AA’ and BB’, as shown in Figure 8. The Profile AA’ is along the terrain gradient descent, approximately parallel to the main sliding direction of Blocks B1-B2, and the Profile BB’ transversely passes through Blocks B2–B3 at similar elevation (see Figure 7d). In Figure 8a, the sliding movement in OSD of Blocks B1–B2 (along AA’) presents a negative correlation with the elevation, that is, displacement in the front part (Block B2) is larger than the posterior part (Block B1). This evidence indicates that Blocks B1–B2 are developed under a pull-type mechanism [73], which can be preliminarily verified by geomorphological features (e.g., tensile cracks and scratches at back edge of B1) presented in Figure 7a. In Figure 8b, evidence illustrates that Block B3 is more stable than Block B2 (along BB’). Moreover, as evidenced in Figure 7b–d, 3D displacements in Blocks B3–B4 are significantly smaller than its above Block B1 and present positive correlations with elevation, which suggest that Blocks B3–B4 is very likely to be deformed with a thrust-type mechanism [73], under the continuous extrusion of Block B1.
Furthermore, Figure 9 shows the optimal sliding direction (OSD) for each pixel of landslide, determined using the estimated 3D displacement rates. It worth noting that OSD vectors indicate that Block B3 is sliding towards the southwest, imposing a direct threaten on the Baobao Village. Under the pressure of Block B1, the sliding movement of Block B3 will be more active. It is strongly recommended to enhance the monitoring, prevention, and control of the Baobao landslide, in order to avoid potential losses of life and property in the surroundings.
Overall, through Figure 7, Figure 8 and Figure 9, an intuitive description of the sliding direction and displacement characteristics of each block is provided for the study of Baobao landslide from a three-dimensional perspective. The spatial pattern and deformation mechanism is analyzed on a detailed scale. In particular in Figure 9, the OSD vectors superimposed on the optical image most intuitively exhibit the potential threat of landslide to surrounding villages, traffics, rivers, etc., playing a positive role in disaster risk management.

5.2. Long-Term Displacement Time Series for Temporal Displacement Pattern Analysis

In order to reconstruct the dynamic evolution of Baobao landslide, four points, P1–P4 (marked in Figure 7b), located in different blocks of the landslide are selected to retrieve the long-term time series of displacement in optimal sliding direction from ascending and descending Sentinel-1 images, according to Step 3.
Figure 10 shows the OSD displacement time series of Points P1-P4 during March 2017 to May 2019, where the triangle and error bar indicate the cumulative displacement at each epoch and its standard deviation calculated with the Kalman filter, respectively. Among them, the largest movement is measured at Point P1, with a magnitude of 0.39 m and a maximum standard deviation of 0.0024 m; the smallest movement is measured at Point P4, with a magnitude of 0.23 m and a maximum standard deviation of 0.0027 m. In addition, the cumulative displacement observed at Points P2 and P3 are 0.37 m and 0.29 m, respectively, with relatively large standard deviations of 0.0096 m and 0.0071 m.
Except for gravity, the primary driving factor, the acceleration of landslide movement may also be caused by some external environment factors, including earthquakes, heavy precipitation, and river and groundwater level fluctuations. Specific to the study area, two recorded earthquakes marked in Figure 1 occurred in 2009 and 2013, respectively, outside the temporal coverage of SAR acquisition. With this in mind, focusing on the correlation between displacement time series and monthly precipitations, the evolution process of Baobao landslide is analyzed in combination with the standard three-stage creep rupture curve model (refer to Figure 14a in [41]). From March 2017 to May 2018, the displacements of Points P1–P4 all present linear evolution trends. Landslide displacement during this period exhibits creep behavior of the secondary creep stage, i.e., constant state with a steady displacement rate. After June 2018, a significant signal of displacement acceleration is observed at all points P1–P4, which may be correlated with heavy rainfalls in flood season. Large amounts of rain infiltration can increase the pore-water pressures of slopes and, thus, reduce the shear strength of materials, resulting in an increase in sliding movement. The Baobao landslide has entered the tertiary creep stage. Then, with the arrival of dry season (from October 2018 to May 2019), precipitation gradually decreased and, thus, the Baobao landslide reached a new limit equilibrium.

6. Conclusions

In this paper, a full-scale workflow is designed for the current increasingly systematic and multi-level InSAR study of landslides, in terms of three levels: (i) early detection on a regional scale, in favor of landslide inventory and risk assessment, (ii) 3D displacement rates estimation on a detailed scale, enhancing understanding on kinematic mechanism and failure mode, and (iii) time series analysis on long-term temporal scale, presenting benefits in landslide forecasting. This workflow is applied for landslide research on the Xiaojiang River Basin, Eastern Yunnan, China, using ascending and descending Sentinel-1 images acquired from March 2017 to May 2019. A relatively complete landslide inventory is mapped and updated for geohazard prevention. New landslides are discovered, and known ones are re-evaluated on their activities and boundaries. Furthermore, three-dimensional and long-term displacement information of Baobao Landslide are estimated through the proposed Total Least Squares–Kalman Filter-based approach. Based on these outcomes, the spatial-temporal pattern is further analyzed to deduce the failure mechanism and current evolution stage.
Through this study, the potential of Sentinel-1 images in InSAR study of landslides has been fully explored and demonstrated. As a rare open resource, Sentinel-1 satellites provide a large amount of free valuable observations since their launch, greatly promoting the development of InSAR technology and its applications in landslide studies. However, Sentinel-1 data with C band are relatively sensitive to seasonal variation in vegetation coverage, often resulting in unsatisfactory performance on landslide monitoring in natural environments with lush flora. Moreover, medium ground resolution of Sentinel-1 data (~15 m) can hardly detect small to medium landslides, let alone to characterize motion mechanism in detail. These two points limit the capability of Sentinel-1 in the InSAR study of landslides to some extent. In contrast, ALOS-2 data with L band and high ground resolution (~3 m under ultra-fine mode) could have played a greater role in this filed. However, its unstable data acquisition results in that the actual observation interval is several times longer than the revisit period (14 days) in most regions. Previous experience indicates that stable data acquisition with high observation frequency is of great significance to closely monitor the landslide development and capture the omen of failure. A hot prospect of new and forthcoming SAR satellites (e.g., LT-1, Tandem-L), with longer wavelengths, higher resolution, and shorter revisit period, is predictable in InSAR study of landslides, focusing on two aspects: (i) regional early detection based on multi-source and multi-mode data integration and (ii) 3D displacement time series estimation and prediction combing high-resolution high-frequency observations and landslide models.

Author Contributions

Conceptualization, H.J., Y.W. and D.G.; methodology, H.J.; software, H.J.; validation, H.J., Y.W. and D.G.; formal analysis, H.J.; investigation, H.J. and D.G.; resources, D.G.; data curation, H.J.; writing—original draft preparation, H.J.; writing—review and editing, H.J. and Y.W.; visualization, H.J.; supervision, Y.W.; project administration, R.W. and Y.W.; funding acquisition, Y.D. and Y.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Funds, China, grant number 61901444.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We acknowledge the European Space Agency (ESA) for providing Sentinel-1 SAR data through open sources. We also acknowledge the China Aero Geophysical Survey & Remote Sensing Center for Land and Resources (AGRS) for providing the available landslide inventory of the study area. We are grateful for the data support from the Geographic Data Sharing Infrastructure, College of Urban and Environmental Science, Peking University. Additionally, we are very grateful to all reviewers, institutions, and studies for their help and advice on our work.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Necula, N.; Niculiță, M.; Fiaschi, S.; Genevois, R.; Riccardi, P.; Floris, M. Assessing Urban Landslide Dynamics through Multi-Temporal InSAR Techniques and Slope Numerical Modeling. Remote Sens. 2021, 13, 3862. [Google Scholar] [CrossRef]
  2. Zhu, Y.; Qiu, H.; Liu, Z.; Wang, J.; Yang, D.; Pei, Y.; Ma, S.; Du, C.; Sun, H.; Wang, L. Detecting Long-Term Deformation of a Loess Landslide from the Phase and Amplitude of Satellite SAR Images: A Retrospective Analysis for the Closure of a Tunnel Event. Remote Sens. 2021, 13, 4841. [Google Scholar] [CrossRef]
  3. Yang, D.; Qiu, H.; Zhu, Y.; Liu, Z.; Pei, Y.; Ma, S.; Du, C.; Sun, H.; Liu, Y.; Cao, M. Landslide Characteristics and Evolution: What We Can Learn from Three Adjacent Landslides. Remote Sens. 2021, 13, 4579. [Google Scholar] [CrossRef]
  4. Petley, D. Global patterns of loss of life from landslides. Geology 2012, 40, 927–930. [Google Scholar] [CrossRef]
  5. Froude, M.; Petley, D. Global fatal landslide occurrence from 2004 to 2016. Nat. Hazards Earth Syst. Sci. 2018, 18, 2161–2181. [Google Scholar] [CrossRef] [Green Version]
  6. Jia, H.; Wang, Y.; Ge, D.; Deng, Y.; Wang, R. Improved Offset Tracking for Predisaster Deformation Monitoring of the 2018 Jinsha River Landslide (Tibet, China). Remote Sens. Environ. 2020, 247, 111899. [Google Scholar] [CrossRef]
  7. Intrieri, E.; Raspini, F.; Fumagalli, A.; Lu, P.; del Conte, S.; Farina, P.; Allievi, J.; Ferretti, A.; Casagli, N. The Maoxian Landslide as Seen from Space: Detecting Precursors of Failure with Sentinel-1 Data. Landslides 2018, 15, 123–133. [Google Scholar] [CrossRef] [Green Version]
  8. Dai, K.; Li, Z.; Xu, Q.; Burgmann, R.; Milledge, D.G.; Tomas, R.; Fan, X.; Zhao, C.; Liu, X.; Peng, J.; et al. Entering the Era of Earth Observation-Based Landslide Warning Systems: A Novel and Exciting Framework. IEEE Geosci. Remote Sens. Mag. 2020, 8, 136–153. [Google Scholar] [CrossRef] [Green Version]
  9. Massonnet, D.; Feigl, K.L. Radar Interferometry and Its Application to Changes in the Earth’s Surface. Rev. Geophys. 1998, 36, 441–500. [Google Scholar] [CrossRef] [Green Version]
  10. Berardino, P.; Fornaro, G.; Lanari, R.; Sansosti, E. A New Algorithm for Surface Deformation Monitoring Based on Small Baseline Differential SAR Interferograms. IEEE Trans. Geosci. Remote Sens. 2002, 40, 2375–2383. [Google Scholar] [CrossRef] [Green Version]
  11. Ferretti, A.; Prati, C.; Rocca, F. Permanent Scatterers in SAR Interferometry. IEEE Trans. Geosci. Remote Sens. 2001, 39, 8–20. [Google Scholar] [CrossRef]
  12. Lyons, S.; Sandwell, D. Fault Creep along the Southern San Andreas from Interferometric Synthetic Aperture Radar, Permanent Scatterers, and Stacking: Southern San Andreas from Insar. J. Geophys. Res. 2003, 108, 108. [Google Scholar] [CrossRef]
  13. Rosi, A.; Tofani, V.; Tanteri, L.; Tacconi Stefanelli, C.; Agostini, A.; Catani, F.; Casagli, N. The New Landslide Inventory of Tuscany (Italy) Updated with PS-InSAR: Geomorphological Features and Landslide Distribution. Landslides 2018, 15, 5–19. [Google Scholar] [CrossRef] [Green Version]
  14. Solari, L.; del Soldato, M.; Raspini, F.; Barra, A.; Bianchini, S.; Confuorto, P.; Casagli, N.; Crosetto, M. Review of Satellite Interferometry for Landslide Detection in Italy. Remote Sens. 2020, 12, 1351. [Google Scholar] [CrossRef]
  15. Schlögel, R.; Doubre, C.; Malet, J.-P.; Masson, F. Landslide Deformation Monitoring with ALOS/PALSAR Imagery: A D-InSAR Geomorphological Interpretation Method. Geomorphology 2015, 231, 314–330. [Google Scholar] [CrossRef]
  16. Aslan, G.; Foumelis, M.; Raucoules, D.; de Michele, M.; Bernardie, S.; Cakir, Z. Landslide Mapping and Monitoring Using Persistent Scatterer Interferometry (PSI) Technique in the French Alps. Remote Sens. 2020, 12, 1305. [Google Scholar] [CrossRef] [Green Version]
  17. Zhao, C.; Lu, Z.; Zhang, Q.; de la Fuente, J. Large-Area Landslide Detection and Monitoring with ALOS/PALSAR Imagery Data over Northern California and Southern Oregon, USA. Remote Sens. Environ. 2012, 124, 348–359. [Google Scholar] [CrossRef]
  18. Zhao, C.; Kang, Y.; Zhang, Q.; Lu, Z.; Li, B. Landslide Identification and Monitoring along the Jinsha River Catchment (Wudongde Reservoir Area), China, Using the InSAR Method. Remote Sens. 2018, 10, 993. [Google Scholar] [CrossRef] [Green Version]
  19. Liu, X.; Zhao, C.; Zhang, Q.; Lu, Z.; Li, Z.; Yang, C.; Zhu, W.; Liu-Zeng, J.; Chen, L.; Liu, C. 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]
  20. Novellino, A.; Cigna, F.; Sowter, A.; Ramondini, M.; Calcaterra, D. Exploitation of the Intermittent SBAS (ISBAS) Algorithm with COSMO-SkyMed Data for Landslide Inventory Mapping in North-Western Sicily, Italy. Geomorphology 2017, 280, 153–166. [Google Scholar] [CrossRef]
  21. Bouali, E.H.; Oommen, T.; Escobar-Wolf, R. Mapping of Slow Landslides on the Palos Verdes Peninsula Using the California Landslide Inventory and Persistent Scatterer Interferometry. Landslides 2018, 15, 439–452. [Google Scholar] [CrossRef]
  22. Zhang, L.; Dai, K.; Deng, J.; Ge, D.; Liang, R.; Li, W.; Xu, Q. Identifying Potential Landslides by Stacking-InSAR in Southwestern China and Its Performance Comparison with SBAS-InSAR. Remote Sens. 2021, 13, 3662. [Google Scholar] [CrossRef]
  23. Colesanti, C.; Wasowski, J. Investigating Landslides with Space-Borne Synthetic Aperture Radar (SAR) Interferometry. Eng. Geol. 2006, 88, 173–199. [Google Scholar] [CrossRef]
  24. Gray, L. Using Multiple RADARSAT InSAR Pairs to Estimate a Full Three-Dimensional Solution for Glacial Ice Movement: Multiple Interferograms for 3-D Motion. Geophys. Res. Lett. 2011, 38, 38. [Google Scholar] [CrossRef]
  25. Wright, T.J. Toward Mapping Surface Deformation in Three Dimensions Using InSAR. Geophys. Res. Lett. 2004, 31, L01607. [Google Scholar] [CrossRef] [Green Version]
  26. Fialko, Y.; Sandwell, D.; Simons, M.; Rosen, P. Three-Dimensional Deformation Caused by the Bam, Iran, Earthquake and the Origin of Shallow Slip Deficit. Nature 2005, 435, 295–299. [Google Scholar] [CrossRef]
  27. Hu, J.; Li, Z.; Zhu, J.; Ren, X.; Ding, X. Inferring Three-Dimensional Surface Displacement Field by Combining SAR Interferometric Phase and Amplitude Information of Ascending and Descending Orbits. Sci. China Earth Sci. 2010, 53, 550–560. [Google Scholar] [CrossRef]
  28. Jung, H.S.; Lu, Z.; Won, J.S.; Poland, M.P.; Miklius, A. Mapping Three-Dimensional Surface Deformation by Combining Multiple-Aperture Interferometry and Conventional Interferometry: Application to the June 2007 Eruption of Kilauea Volcano, Hawaii. IEEE Geosci. Remote Sensing Lett. 2011, 8, 34–38. [Google Scholar] [CrossRef]
  29. Hu, J.; Li, Z.W.; Ding, X.L.; Zhu, J.J.; Zhang, L.; Sun, Q. 3D Coseismic Displacement of 2010 Darfield, New Zealand Earthquake Estimated from Multi-Aperture InSAR and D-InSAR Measurements. J. Geod. 2012, 86, 1029–1041. [Google Scholar] [CrossRef]
  30. Gudmundsson, S.; Sigmundsson, F.; Carstensen, J.M. Three-Dimensional Surface Motion Maps Estimated from Combined Interferometric Synthetic Aperture Radar and GPS Data: THREE-DIMENSIONAL SURFACE MOTION MAPS. J. Geophys. Res. 2002, 107, ETG 13-1–ETG 13-14. [Google Scholar] [CrossRef]
  31. Ji, P.; Lv, X.; Chen, Q.; Sun, G.; Yao, J. Applying InSAR and GNSS Data to Obtain 3-D Surface Deformations Based on Iterated Almost Unbiased Estimation and Laplacian Smoothness Constraint. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2021, 14, 337–349. [Google Scholar] [CrossRef]
  32. Ji, P.; Lv, X.; Yao, J.; Sun, G. A New Method to Obtain 3-D Surface Deformations from InSAR and GNSS Data With Genetic Algorithm and Support Vector Machine. IEEE Geosci. Remote Sens. Lett. 2022, 19, 1–5. [Google Scholar] [CrossRef]
  33. Dong, L.; Wang, C.; Tang, Y.; Tang, F.; Zhang, H.; Wang, J.; Duan, W. Time Series InSAR Three-Dimensional Displacement Inversion Model of Coal Mining Areas Based on Symmetrical Features of Mining Subsidence. Remote Sens. 2021, 13, 2143. [Google Scholar] [CrossRef]
  34. Yang, Z.; Li, Z.; Zhu, J.; Feng, G.; Wang, Q.; Hu, J.; Wang, C. Deriving Time-Series Three-Dimensional Displacements of Mining Areas from a Single-Geometry InSAR Dataset. J. Geod. 2018, 92, 529–544. [Google Scholar] [CrossRef]
  35. Li, Z.W.; Yang, Z.F.; Zhu, J.J.; Hu, J.; Wang, Y.J.; Li, P.X.; Chen, G.L. Retrieving Three-Dimensional Displacement Fields of Mining Areas from a Single InSAR Pair. J. Geod. 2015, 89, 17–32. [Google Scholar] [CrossRef]
  36. Kumar, V.; Venkataramana, G.; Høgda, K.A. Glacier Surface Velocity Estimation Using SAR Interferometry Technique Applying Ascending and Descending Passes in Himalayas. Int. J. Appl. Earth Obs. Geoinf. 2011, 13, 545–551. [Google Scholar] [CrossRef]
  37. Sun, Q.; Hu, J.; Zhang, L.; Ding, X. Towards Slow-Moving Landslide Monitoring by Integrating Multi-Sensor InSAR Time Series Datasets: The Zhouqu Case Study, China. Remote Sens. 2016, 8, 908. [Google Scholar] [CrossRef] [Green Version]
  38. Hu, X.; Bürgmann, R.; Fielding, E.J.; Lee, H. Internal Kinematics of the Slumgullion Landslide (USA) from High-Resolution UAVSAR InSAR Data. Remote Sens. Environ. 2020, 251, 112057. [Google Scholar] [CrossRef]
  39. Ao, M.; Zhang, L.; Shi, X.; Liao, M.; Dong, J. Measurement of the Three-Dimensional Surface Deformation of the Jiaju Landslide Using a Surface-Parallel Flow Model. Remote Sens. Lett. 2019, 10, 776–785. [Google Scholar] [CrossRef]
  40. Samsonov, S.; Dille, A.; Dewitte, O.; Kervyn, F.; d’Oreye, N. Satellite Interferometry for Mapping Surface Deformation Time Series in One, Two and Three Dimensions: A New Method Illustrated on a Slow-Moving Landslide. Eng. Geol. 2020, 266, 105471. [Google Scholar] [CrossRef]
  41. Liu, X.; Zhao, C.; Zhang, Q.; Yin, Y.; Lu, Z.; Samsonov, S.; Yang, C.; Wang, M.; Tomás, R. Three-Dimensional and Long-Term Landslide Displacement Estimation by Fusing C- and L-Band SAR Observations: A Case Study in Gongjue County, Tibet, China. Remote Sens. Environ. 2021, 267, 112745. [Google Scholar] [CrossRef]
  42. Aydan, Ö.; Ito, T.; Özbay, U.; Kwasniewski, M.; Shariar, K.; Okuno, T.; Özgenoğlu, A.; Malan, D.F.; Okada, T. ISRM Suggested Methods for Determining the Creep Characteristics of Rock. Rock Mech. Rock Eng. 2014, 47, 275–290. [Google Scholar] [CrossRef] [Green Version]
  43. Wasowski, J.; Bovenga, F. Investigating Landslides and Unstable Slopes with Satellite Multi Temporal Interferometry: Current Issues and Future Perspectives. Eng. Geol. 2014, 174, 103–138. [Google Scholar] [CrossRef]
  44. Intrieri, E.; Carlà, T.; Gigli, G. Forecasting the Time of Failure of Landslides at Slope-Scale: A Literature Review. Earth-Sci. Rev. 2019, 193, 333–349. [Google Scholar] [CrossRef]
  45. Casu, F.; Manconi, A.; Pepe, A.; Lanari, R. Deformation Time-Series Generation in Areas Characterized by Large Displacement Dynamics: The SAR Amplitude Pixel-Offset SBAS Technique. IEEE Trans. Geosci. Remote Sens. 2011, 49, 2752–2763. [Google Scholar] [CrossRef]
  46. Samsonov, S.; d’Oreye, N. Multidimensional time-series analysis of ground deformation from multiple insar data sets applied to virunga volcanic province. Geophys. J. Int. 2012, 191, 1095–1108. [Google Scholar]
  47. Pepe, A.; Solaro, G.; Calo, F.; Dema, C. A Minimum Acceleration Approach for the Retrieval of Multiplatform InSAR Deformation Time Series. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2016, 9, 3883–3898. [Google Scholar] [CrossRef]
  48. Hu, J.; Ding, X.-L.; Li, Z.-W.; Zhu, J.-J.; Sun, Q.; Zhang, L. Kalman-Filter-Based Approach for Multisensor, Multitrack, and Multitemporal InSAR. IEEE Trans. Geosci. Remote Sens. 2013, 51, 4226–4239. [Google Scholar] [CrossRef]
  49. Dalaison, M.; Jolivet, R. A Kalman Filter Time Series Analysis Method for InSAR. J. Geophys. Res. Solid Earth 2020, 125, 125. [Google Scholar] [CrossRef]
  50. Zhang, X.; Gan, S.; Yuan, X.; Wang, L.; Yang, M.; Chen, X. Analysis of land use patterns of debris flow fans in Xiaojiang River Basin of Yunnan Province. J. Yunnan Univ. Nat. Sci. Ed. 2020, 42, 499–506. [Google Scholar] [CrossRef]
  51. Wu, J.; Zeng, Y. Xiaojiang River Basin: Natural Museum of Landslides. Hum. Nat. 2002, 4, 21–25. [Google Scholar]
  52. Guo, C.; Montgomery, D.R.; Zhang, Y.; Wang, K.; Yang, Z. Quantitative Assessment of Landslide Susceptibility along the Xianshuihe Fault Zone, Tibetan Plateau, China. Geomorphology 2015, 248, 93–110. [Google Scholar] [CrossRef]
  53. Yang, J.; Li, Y.; Yun, N.; Zhou, S.; Yang, H.; Yang, Z.; Yao, Y. Dynamic earthquake trig-gering in the north of Xiaojiang fault zone, Yunnan. Chin. J. Geophys. 2021, 64, 3207–3219. (In Chinese) [Google Scholar] [CrossRef]
  54. Zhao, X. Study on the Risk Assessment of Landslide Hazards in Northeastern Yunnan Province. Master’s Thesis, Yunnan University, Kunming, China, 2018. Available online: https://kns.cnki.net/kcms/detail/detail.aspx?FileName=1018248736.nh&DbName=CMFD2020 (accessed on 1 May 2018). (In Chinese).
  55. Zhao, N.; Wang, Z.; Pan, B.; Li, Z.; Daun, X. Ecological Functions of Riverbed Structures with Different Strengths in the Xiaojiang River Basin. J. Tsinghua Univ. Sci. Technol. 2015, 54, 584–589. [Google Scholar]
  56. Zhou, D.; Zuo, X. Surface subsidence monitoring and prediction in mining area based on SBAS-InSAR and PSO-BP neural network algorithm. J. Yunnan Univ. Nat. Sci. Ed. 2021, 43, 895–905. [Google Scholar] [CrossRef]
  57. Dun, J.; Feng, W.; Yi, X.; Zhang, G.; Wu, M. Detection and Mapping of Active Landslides before Impoundment in the Baihetan Reservoir Area (China) Based on the Time-Series InSAR Method. Remote Sens. 2021, 13, 3213. [Google Scholar] [CrossRef]
  58. Gourmelen, N.; Kim, S.W.; Shepherd, A.; Park, J.W.; Sundal, A.V.; Björnsson, H.; Pálsson, F. Ice Velocity Determined Using Conventional and Multiple-Aperture InSAR. Earth Planet. Sci. Lett. 2011, 307, 156–160. [Google Scholar] [CrossRef] [Green Version]
  59. Cascini, L.; Fornaro, G.; Peduto, D. Advanced Low- and Full-Resolution DInSAR Map Generation for Slow-Moving Landslide Analysis at Different Scales. Eng. Geol. 2010, 112, 29–42. [Google Scholar] [CrossRef]
  60. Wang, Y.; Liu, D.; Dong, J.; Zhang, L.; Guo, J.; Liao, M.; Gong, J. On the Applicability of Satellite SAR Interferometry to Landslide Hazards Detection in Hilly Areas: A Case Study of Shuicheng, Guizhou in Southwest China. Landslides 2021, 18, 2609–2619. [Google Scholar] [CrossRef]
  61. Costantini, M. A Novel Phase Unwrapping Method Based on Network Programming. IEEE Trans. Geosci. Remote Sens. 1998, 36, 813–821. [Google Scholar] [CrossRef]
  62. Liao, M.; Zhang, L.; Shi, X.; Jiang, Y.; Dong, J.; Liu, Y. Radar Remote Sensing Monitoring Methods and Practices of Landslide Deformation; Science Press: Beijing, China, 2017. [Google Scholar]
  63. Zhao, C.; Zhang, Q.; He, Y.; Peng, J.; Yang, C.; Kang, Y. Small-Scale Loess Landslide Monitoring with Small Baseline Subsets Interferometric Synthetic Aperture Radar Technique—Case Study of Xingyuan Landslide, Shaanxi, China. J. Appl. Remote Sens. 2016, 10, 026030. [Google Scholar] [CrossRef] [Green Version]
  64. Bekaert, D.P.S.; Hooper, A.; Wright, T.J. A Spatially Variable Power Law Tropospheric Correction Technique for InSAR Data. J. Geophys. Res. Solid Earth 2015, 120, 1345–1356. [Google Scholar] [CrossRef]
  65. Biggs, J.; Wright, T.; Lu, Z.; Parsons, B. Multi-Interferogram Method for Measuring Interseismic Deformation: Denali Fault, Alaska. Geophys. J. Int. 2007, 170, 1165–1179. [Google Scholar] [CrossRef] [Green Version]
  66. Frattini, P.; Crosta, G.B.; Rossini, M.; Allievi, J. Activity and Kinematic Behaviour of Deep-Seated Landslides from PS-InSAR Displacement Rate Measurements. Landslides 2018, 15, 1053–1070. [Google Scholar] [CrossRef]
  67. Aryal, A.; Brooks, B.A.; Reid, M.E. Landslide Subsurface Slip Geometry Inferred from 3-D Surface Displacement Fields. Geophys. Res. Lett. 2015, 42, 1411–1417. [Google Scholar] [CrossRef]
  68. Rezaei, S.; Shooshpasha, I.; Rezaei, H. Reconstruction of Landslide Model from ERT, Geotechnical, and Field Data, Nargeschal Landslide, Iran. Bull. Eng. Geol. Environ. 2019, 78, 3223–3237. [Google Scholar] [CrossRef]
  69. Wang, B.; Li, J.; Liu, C.; Yu, J. Generalized Total Least Squares Prediction Algorithm for Universal 3D Similarity Transformation. Adv. Space Res. 2017, 59, 815–823. [Google Scholar] [CrossRef]
  70. Schaffrin, B.; Wieser, A. On Weighted Total Least-Squares Adjustment for Linear Regression. J. Geod. 2008, 82, 415–421. [Google Scholar] [CrossRef]
  71. Sage, A.P.; Husa, G.W. Adaptive filtering with unknown prior statistics. In Proceedings of the Joint Automatic Control Conference, Boulder, CO, USA, 5–7 August 1969; Volume 7, pp. 760–769. [Google Scholar]
  72. Welch, G.; Gary, B. An Introduction to the Kalman Filter; University of North Carolina at Chapel Hill: Chapel Hill, NC, USA, 1995; pp. 127–132. Available online: https://www.cs.unc.edu/~welch/media/pdf/kalman_intro.pdf (accessed on 24 July 2006).
  73. Lu, Y. Deformation and Failure Mechanism of Slope in Three Dimensions. J. Rock Mech. Geotech. Eng. 2015, 7, 109–119. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Location of the study area and coverage of the SAR images, with SRTM DEM as the base map. White and black rectangles illustrate the study area and the coverage of SAR images, respectively. Blue lines represent the Jinsha River and its tributary, Xiaojiang River, and black flag icons are hydropower stations. Red lines represent the active faults, where F1: Xiaojiang Fault, F2: Pudu River Fault, F3: Yuanmou Fault. Red dots are the earthquakes that occurred in the study area and vicinity during the period of 2009 to 2019 (information acquired from Unified earthquake catalogue of China Seismic Network).
Figure 1. Location of the study area and coverage of the SAR images, with SRTM DEM as the base map. White and black rectangles illustrate the study area and the coverage of SAR images, respectively. Blue lines represent the Jinsha River and its tributary, Xiaojiang River, and black flag icons are hydropower stations. Red lines represent the active faults, where F1: Xiaojiang Fault, F2: Pudu River Fault, F3: Yuanmou Fault. Red dots are the earthquakes that occurred in the study area and vicinity during the period of 2009 to 2019 (information acquired from Unified earthquake catalogue of China Seismic Network).
Remotesensing 14 01759 g001
Figure 2. (a) Geological setting of the study area (modified from data sharing by Geographic Data Sharing Infrastructure, College of Urban and Environmental Science, Peking University. White curves map the major historical debris flows. Black curves map the recorded landslides (data support from China Aero Geophysical Survey & Remote Sensing Center for Land and Resources (AGRS)). (b) and (c) are enlarged views of two local regions, with SRTM DEM as the base map.
Figure 2. (a) Geological setting of the study area (modified from data sharing by Geographic Data Sharing Infrastructure, College of Urban and Environmental Science, Peking University. White curves map the major historical debris flows. Black curves map the recorded landslides (data support from China Aero Geophysical Survey & Remote Sensing Center for Land and Resources (AGRS)). (b) and (c) are enlarged views of two local regions, with SRTM DEM as the base map.
Remotesensing 14 01759 g002
Figure 3. Workflow of InSAR-driven landslide detection and Total Least Squares–Kalman Filter-based approach for three-dimensional displacement rates and long-term time series estimation of landslides.
Figure 3. Workflow of InSAR-driven landslide detection and Total Least Squares–Kalman Filter-based approach for three-dimensional displacement rates and long-term time series estimation of landslides.
Remotesensing 14 01759 g003
Figure 4. Spatiotemporal distributions of accepted unwrapped interferograms. (a) Ascending Sentinel-1 dataset for Track 26 and (b) descending Sentinel-1 dataset for Track 62.
Figure 4. Spatiotemporal distributions of accepted unwrapped interferograms. (a) Ascending Sentinel-1 dataset for Track 26 and (b) descending Sentinel-1 dataset for Track 62.
Remotesensing 14 01759 g004
Figure 5. Line-of-sight (LOS) displacement rate maps for the study area derived from (a) ascending Sentinel-1 images during March 2017 to May 2019; (b) descending Sentinel-1 images during March 2017 to April 2019. White curves map the major historical debris flows, and black curves map the recorded landslides from AGRS. By our approach, five recorded landslides (No. 14, 16–19) have been detected significant expansion of their boundaries, redrawn within red dotted curves; three unrecorded landslides (No. 20–22) have been newly detected and mapped within red curves. (cg) show enlarged views of these landslides, within their ascending observations, descending observations, and optical images for interpretation.
Figure 5. Line-of-sight (LOS) displacement rate maps for the study area derived from (a) ascending Sentinel-1 images during March 2017 to May 2019; (b) descending Sentinel-1 images during March 2017 to April 2019. White curves map the major historical debris flows, and black curves map the recorded landslides from AGRS. By our approach, five recorded landslides (No. 14, 16–19) have been detected significant expansion of their boundaries, redrawn within red dotted curves; three unrecorded landslides (No. 20–22) have been newly detected and mapped within red curves. (cg) show enlarged views of these landslides, within their ascending observations, descending observations, and optical images for interpretation.
Remotesensing 14 01759 g005aRemotesensing 14 01759 g005b
Figure 6. Accuracy assessment of LOS displacement maps. Distribution histograms of measurement scatterers in stable area are calculated based on (a) ascending Sentinel-1 results and (b) descending Sentinel-1 results.
Figure 6. Accuracy assessment of LOS displacement maps. Distribution histograms of measurement scatterers in stable area are calculated based on (a) ascending Sentinel-1 results and (b) descending Sentinel-1 results.
Remotesensing 14 01759 g006
Figure 7. Optimal remote sensing image and three-dimensional displacement rate maps during March 2017 to May 2019 of the Baobao landslide. (a) Optical image acquired in April 2019. White curves map the boundary of the landslide movement and different colors represent four Blocks B1–B4. (b) Vertical displacement rate map, where Points P1–P4 are selected to analysis their time series of displacement in optimal sliding direction. (c) North-south displacement rate map, where black dotted curves mapped two most active sliding regions (R1 and R2). (d) East-west displacement rate map, where two lines indicate the locations of Profiles AA’ and BB’. Moreover, yellow dotted lines are topographic contour lines.
Figure 7. Optimal remote sensing image and three-dimensional displacement rate maps during March 2017 to May 2019 of the Baobao landslide. (a) Optical image acquired in April 2019. White curves map the boundary of the landslide movement and different colors represent four Blocks B1–B4. (b) Vertical displacement rate map, where Points P1–P4 are selected to analysis their time series of displacement in optimal sliding direction. (c) North-south displacement rate map, where black dotted curves mapped two most active sliding regions (R1 and R2). (d) East-west displacement rate map, where two lines indicate the locations of Profiles AA’ and BB’. Moreover, yellow dotted lines are topographic contour lines.
Remotesensing 14 01759 g007
Figure 8. Elevation, three-dimensional displacement rates, and movement in the optimal sliding direction along the Profiles AA’ and BB’ label in Figure 7d. (a) Profile AA’ and (b) Profile BB’, where B1, B2, and B3 indicate blocks of Baobao landslide labeled in Figure 7a, and where R1 and R2 indicate active regions labeled in Figure 7c.
Figure 8. Elevation, three-dimensional displacement rates, and movement in the optimal sliding direction along the Profiles AA’ and BB’ label in Figure 7d. (a) Profile AA’ and (b) Profile BB’, where B1, B2, and B3 indicate blocks of Baobao landslide labeled in Figure 7a, and where R1 and R2 indicate active regions labeled in Figure 7c.
Remotesensing 14 01759 g008
Figure 9. (a) Horizontal movement vectors of the Baobao landslide and (b) the enlarged horizontal movement vector over the Baobao Village. The base map is the optimal image acquired on April 2019.
Figure 9. (a) Horizontal movement vectors of the Baobao landslide and (b) the enlarged horizontal movement vector over the Baobao Village. The base map is the optimal image acquired on April 2019.
Remotesensing 14 01759 g009
Figure 10. The displacement time series with optimal sliding direction for points P1–P4 (marked in Figure 7b) of the Baobao landslide during March 2017 to May 2019 against monthly precipitation. The error bars indicate the standard deviations of the time series derived using the Kalman filter.
Figure 10. The displacement time series with optimal sliding direction for points P1–P4 (marked in Figure 7b) of the Baobao landslide during March 2017 to May 2019 against monthly precipitation. The error bars indicate the standard deviations of the time series derived using the Kalman filter.
Remotesensing 14 01759 g010
Table 1. Basic parameters of SAR images acquired for the study.
Table 1. Basic parameters of SAR images acquired for the study.
SensorsTrackOrbitHeading (°)Incidence
Angle (°)
Start DateEnd DateNo. of Images
Sentinel-126Ascending−12.520939.650313 March 201707 May 201965
Sentinel-162Descending192.525939.693315 March 201710 April 201952
Table 2. Equations of iterative total least squares estimation presented in [70].
Table 2. Equations of iterative total least squares estimation presented in [70].
Step   1 :   ξ 0 = 0 ;   X 0 = G ^ T Σ V 1 G ^ 1 G ^ T Σ V 1 L ;
Step 2:
λ i = ( Σ V + ( X i T I 2 X i ) Σ 0 ) 1 · L G ^ X i
ξ i = λ i T Σ 0 λ i
X i + 1 = G ^ T ( Σ V + ( X i T I 2 X i ) Σ 0 ) 1 G ^ ξ i I 2 1 · G ^ T ( Σ V + ( X i T I 2 X i ) Σ 0 ) 1 L
Step 3: Repeat Step 2 until
X i + 1 X i 2 < t h r e s 1
1thres is a set according to the required accuracy of 3D displacement rates.
Table 3. Equations of the Kalman filter presented in [72].
Table 3. Equations of the Kalman filter presented in [72].
Time update (“predict”):
x ˜ k = A x ˜ k 1
P k = A P k 1 A T + Q k 1
Measurement update (“correct”):
K k = P k H T H P k H T + R k 1
x ˜ k = x ˜ k + K k L k H x ˜ k
P k = I K k H P k
Table 4. Information of the detected active landslides 1.
Table 4. Information of the detected active landslides 1.
No.Location NameAspectSlope (°)Area (ha)ScaleDetected from SAR ImagesRisk LevelThreat Object
1West of XiaoyanjiaoN, NE14–4970Super largeS1A, S1DMediumRoads, Jinsha River
2East of XiaoyanjiaoN, NW17–5271Super largeS1A, S1DMediumRoads, Jinsha River
3lijialiangziNW12–4882Super largeS1AMediumVillages on back edge, farmlands, roads | debris flow source
4LaoliukouSE32–4614LargeS1AMediumRoads on back edge, buildings, farmlands and roads in gully
5Qingmen No.1SW23–4715LargeS1DHighVillages and roads on back edge
6Qingmen No.2SW22–396LargeS1DMediumFarmlands on back edge
7BainijingSW24–393MediumS1A, S1DLowDebris flow source
8Dingjiabaob-aoSW, W6–3531LargeS1A, S1DMediumFarmlands, and roads
9SongpingziNW12–59159Super largeS1A, S1DHighVillages, farmlands, and roads
10Laoqing GullyN28–4120LargeS1DMediumFarmlands on back edge
11Taiping VillageW18–3814Large MediumFarmlands, and roads
12Pingzi VillageN8–4145LargeS1AMediumFarmlands on back edge, roads in gully
13WujiawanE16–4320LargeS1AMediumRoads
14Baobao VillageNW17–47140Super largeS1A, S1DMediumVillages, and roads
15HeinaokeSE13–5167Super largeS1AMediumFarmlands on back edge, buildings, roads, and bridges in gully | debris flow source
16JilongaoNW, N9–49167Super largeS1A, S1DMediumFarmlands
17Sanjiacun-XiaogouqingE, SE8–41275Super largeS1A, S1DHighBuildings, roads, bridges and farmlands in gully | debris flow source
18
19LanshanaoNW9–36175Super largeS1A, S1DMediumRoads, farmland | debris flow source
20GelepingziW6–5638LargeS1A, S1DMediumRoads
21Dongfangho-ng GullySE6–5262Super largeS1AHighRoads, farmlands, and buildings in gully | debris flow source
22Laoao GullyW, SW6–45134Super largeS1A. S1DLowDebris flow source
1 Notes: S1A and S1D represent ascending and descending Sentinel-1 SAR images, respectively. Newly detected landslides (No. 20–22) and recorded landslides with an obvious expansion (No. 14, 16–19) are highlighted in red.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Jia, H.; Wang, Y.; Ge, D.; Deng, Y.; Wang, R. InSAR Study of Landslides: Early Detection, Three-Dimensional, and Long-Term Surface Displacement Estimation—A Case of Xiaojiang River Basin, China. Remote Sens. 2022, 14, 1759. https://doi.org/10.3390/rs14071759

AMA Style

Jia H, Wang Y, Ge D, Deng Y, Wang R. InSAR Study of Landslides: Early Detection, Three-Dimensional, and Long-Term Surface Displacement Estimation—A Case of Xiaojiang River Basin, China. Remote Sensing. 2022; 14(7):1759. https://doi.org/10.3390/rs14071759

Chicago/Turabian Style

Jia, Hongying, Yingjie Wang, Daqing Ge, Yunkai Deng, and Robert Wang. 2022. "InSAR Study of Landslides: Early Detection, Three-Dimensional, and Long-Term Surface Displacement Estimation—A Case of Xiaojiang River Basin, China" Remote Sensing 14, no. 7: 1759. https://doi.org/10.3390/rs14071759

APA Style

Jia, H., Wang, Y., Ge, D., Deng, Y., & Wang, R. (2022). InSAR Study of Landslides: Early Detection, Three-Dimensional, and Long-Term Surface Displacement Estimation—A Case of Xiaojiang River Basin, China. Remote Sensing, 14(7), 1759. https://doi.org/10.3390/rs14071759

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