Exciting news! We're transitioning to the Statewide California Earthquake Center. Our new website is under construction, but we'll continue using this website for SCEC business in the meantime. We're also archiving the Southern Center site to preserve its rich history. A new and improved platform is coming soon!

Group B, Poster #064, Tectonic Geodesy

Joint Inversion of GNSS and InSAR Data for Continuous 3-D Velocity and Strain Rate Fields in Southern California

Mradula Vashishtha, Jeonghyeop Kim, & William E. Holt
Poster Image: 

Poster Presentation

2023 SCEC Annual Meeting, Poster #064, SCEC Contribution #13235 VIEW PDF
GNSS and InSAR complement each other to provide more information on the spatial gradients of crustal motions. GNSS stations hold precision for the three components of crustal velocities but are spatially limited, while InSAR data provide a more continuous spatial coverage but with inherent ambiguities in the three components owing to surface motions measured along a single line of sight for each pixel.

We present a new joint inversion algorithm for GNSS and InSAR data to obtain a 3-D surface velocity field for Southern California as a product for SCEC CGM. Contrary to prior treatments we generate basis function responses to horizontal body-force equivalent inputs in a spheric...
al treatment. These responses satisfy force-balance equations. Vertical basis function responses are weakly coupled to the horizontal strain rates through the Vr/R term, where Vr is the vertical velocity and R is the radius of the Earth; these vertical responses are applicable to represent motions linked to both elastic and poroelastic processes. This suggests that these basis functions can also serve as fitting functions, and thus we can apply this algorithm to study various surface processes such as groundwater signals. We use UCERF3 velocities as boundary conditions around 1° by 1° areas. The boundary condition solution is added to the optimal set of basis function responses as a linear superposition. Together, these represent an estimate of steady-state deformation. Solutions for 1° by 1° areas are combined to represent a region-wide solution for Southern California.

We investigate both isotropic and anisotropic responses to the body-force equivalent inputs. We use Ridge (L2-norm) regularization and Akaike Bayesian Information Criterion to estimate an optimal damping level as well as a relative weighting of InSAR with respect to GNSS data. We show through tests with synthetic data that the joint inversion provides improved estimates for the surface velocity gradient tensor field, with anisotropic treatment yielding the highest quality match to horizontal velocity gradients. We jointly invert ascending and descending tracks of InSAR data from Xu et al. (2021) with GNSS in ITRF frame to obtain a solution for steady-state horizontal strain rates, rotation rates, vertical velocities, and horizontal gradients in vertical velocities in Southern California. We align adjoining regions together using simple interpolation techniques.