Constraining Earthquake Source Inversions with GPS Data 1: Resolution Based Removal of Artifacts

Morgan T. Page, Susana Custodio, Ralph J. Archuleta, & Jean M. Carlson

Published 2009, SCEC Contribution #1127

We present a resolution analysis of an inversion of GPS data from the 2004 M<sub>w</sub> 6.0 Parkfield Earthquake. This earthquake was recorded at 13 1-Hz GPS receivers, which provides for a truly co-seismic dataset that can be used to infer the static-slip field. We find that the resolution of our inverted slip model is poor at depth and near the edges of the modeled fault plane that are far from GPS receivers. The spatial heterogeneity of the model resolution in the static field inversion leads to artifacts in poorly resolved areas of the fault plane. These artifacts look qualitatively similar to asperities commonly seen in the final slip models of earthquake source inversions, but in this inversion they are caused by a surplus of free parameters. The location of the artifacts depends on the station geometry and the assumed velocity structure. We demonstrate that a nonuniform gridding of model parameters on the fault can remove these artifacts from the inversion. We generate a nonuniform grid with a grid spacing that matches the local resolution length on the fault, and show that it outperforms uniform grids, which either generate spurious structure in poorly resolved regions or lose recoverable information in well-resolved areas of the fault. In a synthetic test, the nonuniform grid correctly averages slip in poorly resolved areas of the fault while recovering small-scale structure near the surface. Finally, we present an inversion of the Parkfield GPS dataset on the nonuniform grid and analyze the errors in the final model.

Citation
Page, M. T., Custodio, S., Archuleta, R. J., & Carlson, J. M. (2009). Constraining Earthquake Source Inversions with GPS Data 1: Resolution Based Removal of Artifacts. Journal of Geophysical Research: Solid Earth, 114. doi: 10.1029/2007JB005449.