3-D Satellite Interferometry for Interseismic Velocity Fields

The global interseismic strain rate map is being accomplished rapidly with measurements of the space-based geodetic technique of InSAR. High-resolution measurements of crustal deformation from InSAR can provide crucial constraints on a region's active tectonics, geodynamics, and seismic hazard. However, space-based InSAR usually only provides good constraints on horizontal displacement in the east-west direction, with the north-south component typically provided by low-resolution GNSS measurements. Sentinel-1, on the other hand, has the potential to provide measurements that are sensitive to north-south motion, through exploitation of the burst overlap areas produced by the TOPS acquisition mode. However, the significant noise contributions from decorrelation and propagation through the ionosphere make it challenging to detect surface displacements associated with interseismic deformation needing millimeters per year accuracy. The ionospheric phase advance is a significant nuisance term that can bias InSAR measurements. Although methods have been developed to mitigate the effect, they are not always routinely applied when processing C-band SAR images, for which the effect is generally expected to be small. Nevertheless, the effect can be significant, especially when analyzing low deformation gradients over large areas using time-series analysis. Here, the work in Chapter 3 presents a time-series approach to ionospheric noise mitigation, which improves on existing methods. Firstly, I estimate the ionospheric contribution for each individual acquisition from multiple interferograms, which reduces noise. Secondly, this work improved the identification of unwrapping errors, which can bias the estimation. Thirdly, I introduce a new filtering approach, which gives better results, particularly at image edges and areas with variable density of coherent measurements. Furthermore, the approach is applicable when estimating along-track motion in burst overlap areas. The results show that applying the correction improves velocity accuracy significantly for both conventional line-of-sight and burst overlap interferometry techniques. The application of measuring long-term tectonic signals that concentrate in the north-south component with millimeters per year accuracy is essential to constrain interseismic strain globally. In Chapter 4, I also demonstrate a time-series approach with the burst overlap interferometry appropriate for extracting subtle long-term displacements. The approach includes mitigation of ionospheric noise, and I investigate different filtering approaches to optimize the reduction of decorrelation noise. I present the mean ground velocity in the azimuth direction from data acquired between 2014 and 2019 along the West-Lut Fault, a north-south striking fault in eastern Iran. The chi-square statistic defines a good agreement between the results and independent GNSS measurements. Moreover, the denser coverage of the technique allows to detect the variation in strain accumulation between northern and southern segments of the fault, with our modeling indicating a variation of slip rate from 9.2±0.5 mm/yr in the south to 4.3±0.5 mm/yr in the north. With current efforts to use InSAR to constrain strain rates globally, along-track measurements can fill a crucial gap in north-south sensitivity. With the achievement of that the burst overlap InSAR technique can measure azimuth motions across a slowly deforming area where the surface displacements are concentrated in the north-south component, this results in that, in the TOPS burst overlap region, the number of observations for a ground displacement can reach 3-4 times with different observational components. Measurement redundancy allows for the decomposition of observed velocities into three-dimensional components. In Chapter 5, I apply InSAR observations to estimate a deformation across the Chaman fault in both line-of-sight and along-track components using images from ascending and descending passes. I demonstrate an inversion to estimate the decomposed velocities. The algorithm employs a sparse GNSS network across the region to transform InSAR velocities to the GNSS reference frame. The results show that constraining the long-wavelength signal across the InSAR observations using GNSS data can mitigate the long-wavelength ionospheric disturbance that remains in the observations. The variation in slip rates across the Chaman fault is depicted by two transect profiles. The mean velocity profile at latitude 31˚N, where the Chaman fault is the only tectonic structure to accommodate strain, is consistent with 10.4±0.4 mm/yr of slip rate derived from the interseismic modeling. The optimal fault slip rate to fit with the mean velocity of the southern profile at latitude 29˚N is 5.5±0.8 mm/yr across the Chaman fault and 15.5±0.9 mm/yr across the parallel fault (the Ghazaband fault). I also demonstrate the benefits of high temporal sampling of InSAR observations with TOPS acquisition mode to study time-dependent surface deformation. I present the evolution of fault creeps, including seismic and aseismic fault slip along the Chaman fault during 2014-2018.