Show simple item record

dc.contributor.authorVasyura-Bathke, Hannes
dc.contributor.authorDettmer, J
dc.contributor.authorDutta, Rishabh
dc.contributor.authorMai, Paul Martin
dc.contributor.authorJonsson, Sigurjon
dc.date.accessioned2021-02-07T08:27:29Z
dc.date.available2021-02-07T08:27:29Z
dc.date.issued2021-01-25
dc.identifier.citationVasyura-Bathke, H., Dettmer, J., Dutta, R., Mai, P. M., & Jónsson, S. (2021). Accounting for theory errors with empirical Bayesian noise models in nonlinear centroid moment tensor estimation. Geophysical Journal International. doi:10.1093/gji/ggab034
dc.identifier.issn0956-540X
dc.identifier.issn1365-246X
dc.identifier.doi10.1093/gji/ggab034
dc.identifier.urihttp://hdl.handle.net/10754/667237
dc.description.abstractSummary Centroid moment-tensor (CMT) parameters can be estimated from seismic waveforms. Since these data indirectly observe the deformation process, CMTs are inferred as solutions to inverse problems which are generally under-determined and require significant assumptions, including assumptions about data noise. Broadly speaking, we consider noise to include both theory and measurement errors, where theory errors are due to assumptions in the inverse problem and measurement errors are caused by the measurement process. While data errors are routinely included in parameter estimation for full CMTs, less attention has been paid to theory errors related to velocity-model uncertainties and how these affect the resulting moment-tensor (MT) uncertainties. Therefore, rigorous uncertainty quantification for CMTs may require theory-error estimation which becomes a problem of specifying noise models. Various noise models have been proposed, and these rely on several assumptions. All approaches quantify theory errors by estimating the covariance matrix of data residuals. However, this estimation can be based on explicit modelling, empirical estimation, and/or ignore or include covariances. We quantitatively compare several approaches by presenting parameter and uncertainty estimates in non-linear full CMT estimation for several simulated data sets and regional field data of the Ml 4.4, 13 June 2015 Fox Creek, Canada, event. While our main focus is at regional distances, the tested approaches are general and implemented for arbitrary source model choice. These include known or unknown centroid locations, full MTs, deviatoric MTs, and double-couple MTs. We demonstrate that velocity-model uncertainties can profoundly affect parameter estimation and that their inclusion leads to more realistic parameter uncertainty quantification. However, not all approaches perform equally well. Including theory errors by estimating non-stationary (non-Toeplitz) error covariance matrices via iterative schemes during Monte Carlo sampling performs best and is computationally most efficient. In general, including velocity-model uncertainties is most important in cases where velocity structure is poorly known.
dc.description.sponsorshipSeismic waveforms were originally downloaded from Incorporated Research Institutions for Seismology (IRIS) https://ds.iris.edu. We thank Babak Hejrani, Simon St¨ahler and the editor Andrew Valentine for their comments that helped to improve the quality of this manuscript. Special thanks go to Sebastian Heimann, Frank Kr¨uger, Mehdi Nikkhoo and Olaf Zielke for valuable discussions. Plots were produced with Matplotlib and the Generic Mapping Tools (Hunter 2007; Wessel et al. 2013). This work employed the open source library py rocko (www.pyrocko.org) (Heimann et al. 2017) and the Bayesian Earthquake Analysis Tool (https://github.com/hvasbath/beat). The research was supported by King Abdullah University of Science and Technology (KAUST), under award numbers BAS/1/1353-01-01 and BAS/1/1339-01-1. H.V-B was partially supported by Geo.X, the Research Network for Geosciences in Berlin and Potsdam under the project number SO 087 GeoX.
dc.publisherOxford University Press (OUP)
dc.relation.urlhttps://academic.oup.com/gji/advance-article/doi/10.1093/gji/ggab034/6119912
dc.rightsThis is an Open Access article distributed under the terms of the Creative Commons Attribution License , which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited.
dc.rights.urihttp://creativecommons.org/licenses/by/4.0/
dc.titleAccounting for theory errors with empirical Bayesian noise models in nonlinear centroid moment tensor estimation
dc.typeArticle
dc.contributor.departmentComputational Earthquake Seismology (CES) Research Group
dc.contributor.departmentCrustal Deformation and InSAR Group
dc.contributor.departmentEarth Science and Engineering Program
dc.contributor.departmentPhysical Science and Engineering (PSE) Division
dc.identifier.journalGeophysical Journal International
dc.eprint.versionPublisher's Version/PDF
dc.contributor.institutionNow at: University of Potsdam, D-14476 Potsdam, Germany
dc.contributor.institutionUniversity of Calgary, Calgary, Canada
kaust.personVasyura-Bathke, Hannes
kaust.personDutta, Rishabh
kaust.personMai, Paul Martin
kaust.personJonsson, Sigurjon
kaust.grant.numberBAS/1/1339-01-1
kaust.grant.numberBAS/1/1353-01-01
refterms.dateFOA2021-02-07T08:29:20Z
dc.date.published-online2021-01-25
dc.date.published-print2021-03-12


Files in this item

Thumbnail
Name:
ggab034.pdf
Size:
3.370Mb
Format:
PDF
Description:
Published version

This item appears in the following Collection(s)

Show simple item record

This is an Open Access article distributed under the terms of the Creative Commons Attribution License , which permits unrestricted reuse,
distribution, and reproduction in any medium, provided the original work is properly cited.
Except where otherwise noted, this item's license is described as This is an Open Access article distributed under the terms of the Creative Commons Attribution License , which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited.