Deriving 3-D Surface Deformation Time Series with Strain Model and Kalman Filter from GNSS and InSAR Data

Abstract

This study proposes a new set of processing procedures based on the strain model and the Kalman filter (SM-Kalman) to obtain high-precision three-dimensional surface deformation time series from interferometric synthetic aperture radar (InSAR) and global navigation satellite system (GNSS) data. Implementing the Kalman filter requires the establishment of state and observation equations. In the time domain, the state equation is generated by fitting the pre-existing deformation time series based on a deformation model containing linear and seasonal terms. In the space domain, the observation equation is established with the assistance of the strain model to realize the spatial combination of InSAR and GNSS observation data at each moment. Benefiting from the application of the Kalman filter, InSAR and GNSS data at different moments can be synchronized. The time and measurement update steps are performed dynamically to generate a 3-D deformation time series with high precision and a high resolution in the temporal and spatial domains. Sentinel-1 SAR and GNSS datasets in the Los Angeles area are used to verify the effectiveness of the proposed method. The datasets include twenty-seven ascending track SAR images, thirty-four descending track SAR images and the daily time series of forty-eight GNSS stations from January 2016 to November 2018. The experimental result demonstrates that the proposed SM-Kalman method can produce high-precision deformation results at the millimeter level and provide two types of 3-D deformation time series with the same temporal resolution as InSAR or GNSS observations according to the needs of users. The new method achieves a high degree of temporal and spatial fusion of GNSS and InSAR data.