Three-dimensional and long-term landslide displacement estimation by fusing C- and L-band SAR observations: A case study in Gongjue County, Tibet, China

Recently, a large number of synthetic aperture radar (SAR) images has been introduced into landslide investigations with the growing launch of new SAR satellites, such as ALOS/PALSAR-2 and Sentinel-1. Therefore, it is appropriate to develop new approaches to retrieve three-dimensional (3D) displacements and long-term (> 10 years) displacement time series to investigate the spatio-temporal evolution and creep behavior of landslides. In this study, a new approach for the estimation of 3D and long-term displacement time series of landslides, based on the fusion of C- and L-band SAR observations, is presented. This method is applied to map 3D and long-term displacements (nearly 12 years) of the landslides in Gongjue County, Tibet in China; four sets of SAR images from different platforms (i.e., L-band ascending ALOS/PALSAR-1, C-band descending ENVISAT, and C-band ascending and descending Sentinel-1 SAR datasets) covering the period of January 2007 to November 2018 were collected and exploited. First, the assumption that the landslide moves parallel to its ground surface is used to produce 3D displacement rates and time series by fusing ascending and descending Sentinel-1 SAR images, from which the optimal sliding direction for each pixel of the slope is well estimated. Then, the long-term displacement time-series of the landslide between January 2007 and October 2018 in the estimated sliding direction is recovered by fusing L-band ALOS/PALSAR-1 and C-band Sentinel-1 SAR images. In order to fill the time gap of nearly four years between ALOS/PALSAR-1 and Sentinel-1 SAR images, the Tikhonov regularization (TR) method is developed to establish the observational equation. Moreover, to solve the problem arising from ALOS/PALSAR-1 and Sentinel-1 images with different wavelengths, incidence angles and flight directions, the measurements from ALOS/PALSAR-1 and Sentinel-1 images are both projected to the estimated optimal sliding direction to achieve a unified displacement datum. Our results from ascending and descending Sentinel-1 images suggest that the maximum displacement rates of the study area in the vertical and east-west directions from December 2016 to October 2018 were greater than 70 and 80 mm/year, respectively, and 2D displacement results reveal that the displacement patterns and movement characteristics of all the detected landslides are not identical in the study area. Specifically, the 3D displacement results successfully revealed the spatiotemporal dis