Three-dimensional deformation and its uncertainty derived by integrating multiple SAR data analysis methods

Three-dimensional (3D) surface deformation data with high accuracy and resolution can help reveal the complex mechanisms and sources of subsurface deformation, both tectonic and anthropogenic. Detailed 3D deformation data are also beneficial for maintaining the position coordinates of existing ground features, which is critical for developing and advancing global positioning technologies and their applications. In seismically active regions, large earthquakes have repeatedly caused significant ground deformation and widespread damage to human society. However, the delay in updating position coordinates following deformation can hamper disaster recovery. Synthetic aperture radar (SAR) data allow high-accuracy and high-resolution 3D deformation measurements. Three analysis methods are currently available to measure 1D or 2D deformation: SAR interferometry (InSAR), split-bandwidth interferometry (SBI), and the pixel offset method. In this paper, we propose an approach to derive 3D deformation by integrating deformation data from the three methods. The theoretical uncertainty of the derived 3D deformations was also estimated using observed deformation data for each of these methods and the weighted least square (WLS) approach. Furthermore, we describe two case studies (the 2016 Kumamoto earthquake sequence and the 2016 Central Tottori earthquake in Japan) using L-band Advanced Land Observing Satellite 2 (ALOS-2) data. The case studies demonstrate that the proposed approach successfully retrieved 3D coseismic deformation with the standard error of similar to 1, similar to 4, and similar to 1 cm in the east-west, north-south, and vertical components, respectively, with sufficient InSAR data. SBI and the pixel offset method filled the gaps of the InSAR data in large deformation areas in the order of 10 cm accuracy. The derived standard errors for each pixel are also useful for subsequent applications, such as updating position coordinates and deformation source modeling. The proposed approach is also applicable to other SAR datasets. In particular, next-generation L-band SAR satellites, such as ALOS-4 and NASA-ISRO SAR (NISAR), which have a wider swath width, more frequent observation capabilities than the former L-band satellites, and exclusive main look directions (i.e., right and left) will greatly enhance the applicability of 3D deformation derivation and support the quick recovery from disasters with significant deformation.