Accounting for theory errors with empirical Bayesian noise models in nonlinear centroid moment tensor estimation

Abstract

Summary 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.

A self-similar dynamic rupture model based on the simplified wave-rupture analogy

“…Although several Bayesian inversion approaches were recently developed to handle the uncertainties in finite fault inversion (e.g. Minson et al 2013;Shimizu et al 2020, Heimann et al 2018, including uncertainties related to Green's function (Yagi & Fukahata 2011, Vasyura-Bathke et al 2021, the effort of such methods is usually large which limits the approaches to the later detailed analysis long after the events.…”

Analytical solutions for gravity changes caused by triaxial volumetric sources

“…Some complexities such as layering or viscoelasticity can be accounted for (Amoruso et al, 2008) by using appropriate Green's functions for point dislocations (Okubo, 1993;Sun & Okubo, 1993;Wang et al, 2006). Besides, theory errors, arising from ignoring real Earth complexities, can be estimated in terms of noise covariance matrices within a Bayesian framework (see Duputel et al, 2014;Minson et al, 2013;Vasyura-Bathke et al, 2021).…”

Estimation of Seismic Moment Tensors Using Variational Inference Machine Learning

“…Finally, errors in the Earth structure model can be taken into account for robust inference on the estimated source mechanism (Vasyura‐Bathke et al., 2021). The theory error from the choice of the 1‐D Earth structure model can be included in our framework by training and evaluating BNNs on each grid point for each Earth structure.…”