An EPIC Tikhonov Regularization: Application to Quasi-Static Fault Slip Inversion
Abstract
Imaging subsurface fault slip from surface observations is essential to improving our understanding of the physics of earthquakes and tsunamis. As the estimation of subsurface fault slip is inherently ill-posed, common inversion methods usually require a regularization term to counteract instabilities. Such regularization introduces biases in inferred slip estimates. Here, we discuss the effects that prior information, implied by a given regularization scheme, can have on fault slip estimates. We propose a novel Equal Posterior Information Condition (EPIC)—based Tikhonov regularization that generalizes the concept of prior information. The EPIC determines variances of prior information based on a chosen form of the structure of the posterior covariance matrix. In the context of subduction zone earthquakes, use of the EPIC counterbalances the spatial heterogeneity of observational constraints on fault slip, improving stability, quality and interpretability of fault slip estimates. We investigate the efficiency of the EPIC in the context of various synthetic fault slip distributions. We also demonstrate the methodology by inferring co-seismic slip from geodetic data for the 2011 (M_w 9.0) Tohoku-Oki, Japan, earthquake, obtaining robust slip estimates that are similar to those inferred using a unregularized fully Bayesian approach.
Additional Information
© 2021 American Geophysical Union. Issue Online: 28 June 2021; Version of Record online: 28 June 2021; Accepted manuscript online: 04 June 2021; Manuscript accepted: 02 June 2021; Manuscript revised: 20 May 2021; Manuscript received: 08 October 2020. F. Ortega-Culaciati acknowledges support from Proyecto Fondecyt 11140904 and 1181479, as well as from ANID PIA ACT 192169. J. Ruiz acknowledges support from Proyecto Fondecyt 1200679, and CONICYT PIA grant ACT172002. Powered@NLHPC: This research was partially supported by the supercomputing infrastructure of the NLHPC (ECM-02). We thank the editors and anonymous reviewers for providing insightful comments that helped to improve the clarity and quality of this manuscript. Data Availability Statement: Data sets for this research are included in this paper (and its supplementary information files): Minson et al. (2014). Python codes to solve the general linear least squares inversion problem using EPIC Tikhonov are available at https://doi.org/10.5281/zenodo.4922933, as well as in https://github.com/frortega/EPIC_LS.Attached Files
Published - 2020JB021141.pdf
Supplemental Material - 2020jb021141-sup-0002-supporting_information_si-s01.pdf
Files
Name | Size | Download all |
---|---|---|
md5:6f5ac123e65c32e3f22341ab82d88750
|
4.2 MB | Preview Download |
md5:7db7ba43a3606ae03cb3cfd5f8e82e4f
|
56.5 MB | Preview Download |
Additional details
- Eprint ID
- 109413
- Resolver ID
- CaltechAUTHORS:20210607-115053217
- Fondo Nacional de Desarrollo Científico y Tecnológico (FONDECYT)
- 11140904
- Fondo Nacional de Desarrollo Científico y Tecnológico (FONDECYT)
- 1181479
- Comisión Nacional de Investigación Científica y Tecnológica (CONICYT)
- ACT192169
- Fondo Nacional de Desarrollo Científico y Tecnológico (FONDECYT)
- 1200679
- Comisión Nacional de Investigación Científica y Tecnológica (CONICYT)
- ACT172002
- Created
-
2021-06-07Created from EPrint's datestamp field
- Updated
-
2023-10-03Created from EPrint's last_modified field
- Caltech groups
- Seismological Laboratory, Division of Geological and Planetary Sciences (GPS)