Fully probabilistic seismic source inversion – Part 2: Modelling errors and station covariances

Abstract. Seismic source inversion, a central task in seismology, is concerned with the estimation of earthquake source parameters and their uncertainties. Estimating uncertainties is particularly challenging because source inversion is a non-linear problem. In a companion paper, Stahler and Sigloch (2014) developed a method of fully Bayesian inference for source parameters, based on measurements of waveform cross-correlation between broadband, teleseismic body-wave observations and their modelled counterparts. This approach yields not only depth and moment tensor estimates but also source time functions. A prerequisite for Bayesian inference is the proper characterisation of the noise afflicting the measurements, a problem we address here. We show that, for realistic broadband body-wave seismograms, the systematic error due to an incomplete physical model affects waveform misfits more strongly than random, ambient background noise. In this situation, the waveform cross-correlation coefficient CC, or rather its decorrelation D = 1 − CC, performs more robustly as a misfit criterion than lp norms, more commonly used as sample-by-sample measures of misfit based on distances between individual time samples. From a set of over 900 user-supervised, deterministic earthquake source solutions treated as a quality-controlled reference, we derive the noise distribution on signal decorrelation D = 1 − CC of the broadband seismogram fits between observed and modelled waveforms. The noise on D is found to approximately follow a log-normal distribution, a fortunate fact that readily accommodates the formulation of an empirical likelihood function for D for our multivariate problem. The first and second moments of this multivariate distribution are shown to depend mostly on the signal-to-noise ratio (SNR) of the CC measurements and on the back-azimuthal distances of seismic stations. By identifying and quantifying this likelihood function, we make D and thus waveform cross-correlation measurements usable for fully probabilistic sampling strategies, in source inversion and related applications such as seismic tomography.

[2]  Steven M. Day,et al.  Misfit Criteria for Quantitative Comparison of Seismograms , 2006 .

[3]  M. Sambridge.  Geophysical inversion with a neighbourhood algorithm—II. Appraising the ensemble , 1999 .

[4]  P. Mahalanobis.  On the generalized distance in statistics , 1936 .

[5]  Samuel Kotz,et al.  The Laplace Distribution and Generalizations: A Revisit with Applications to Communications, Economics, Engineering, and Finance , 2001 .

[6]  Martin Schimmel,et al.  Phase Cross-Correlations: Design, Comparisons, and Applications , 1999 .

[7]  É. Stutzmann,et al.  Global climate imprint on seismic noise , 2009 .

[8]  J. Beck,et al.  Accounting for prediction uncertainty when inferring subsurface fault slip , 2014 .

[9]  Simon Stähler,et al.  Fully probabilistic seismic source inversion - Part 1: Efficient parameterisation , 2013 .

[10]  Karin Sigloch,et al.  ObsPyLoad: A Tool for Fully Automated Retrieval of Seismological Waveform Data , 2013 .

[11]  A. Owen.  Empirical likelihood ratio confidence intervals for a single functional , 1988 .

[12]  C. H. Chapman,et al.  A new method for computing synthetic seismograms , 1978 .

[13]  B. Kennett,et al.  Traveltimes for global earthquake location and phase identification , 1991 .

[14]  Masayuki Kikuchi,et al.  Inversion of complex body waves—III , 1991, Bulletin of the Seismological Society of America.

[15]  J. Peterson,et al.  Observations and modeling of seismic background noise , 1993 .

[16]  Carl Tape,et al.  Adjoint Tomography of the Southern California Crust , 2009, Science.

[17]  Vincent Douet,et al.  A new database of source time functions (STFs) extracted from the SCARDEC method , 2016 .

[18]  W. Menke,et al.  Polarization and coherence of 5 to 30 Hz seismic wave fields at a hard-rock site and their relevance to velocity heterogeneities in the crust , 1990, Bulletin of the Seismological Society of America.

[19]  A uniform parametrization of moment tensors , 2015 .

[20]  A.,et al.  Inverse Problems = Quest for Information , 2022 .

[21]  J. Kummerow.  Using the value of the crosscorrelation coefficient to locate microseismic events , 2010 .

[22]  M. Vallée,et al.  SCARDEC: a new technique for the rapid determination of seismic moment magnitude, focal mechanism and source time functions for large earthquakes using body‐wave deconvolution , 2011 .

[23]  D. Storchak,et al.  Improved location procedures at the International Seismological Centre , 2011 .

[24]  G. Nolet,et al.  Measuring finite‐frequency body‐wave amplitudes and traveltimes , 2006 .

[25]  B. Kennett,et al.  Seismic Event Location: Nonlinear Inversion Using a Neighbourhood Algorithm , 2001 .

[26]  P. Shearer,et al.  Shear and compressional velocity models of the mantle from cluster analysis of long‐period waveforms , 2008 .

[27]  A. Fichtner,et al.  Probabilistic full waveform inversion based on tectonic regionalization - development and application to the Australian upper mantle , 2013 .

[28]  T. Nissen‐Meyer,et al.  Triplicated P-wave measurements for waveform tomography of the mantle transition zone , 2012 .

[29]  Masayuki Kikuchi,et al.  Inversion of complex body waves , 1982 .

[30]  T. Bodin.  Transdimensional approaches to geophysical inverse problems , 2010 .

[31]  Karin Sigloch,et al.  Mantle provinces under North America from multifrequency P wave tomography , 2011 .

[32]  S. Kotz,et al.  The Laplace Distribution and Generalizations , 2012 .

[33]  Karin Sigloch,et al.  Intra-oceanic subduction shaped the assembly of Cordilleran North America , 2013, Nature.

[34]  Long-Path Interferometry through an Uncontrolled Atmosphere* , 1962 .

[35]  B. P. Bogert.  Correction of seismograms for the transfer function of the seismometer , 1962 .

[36]  D. L. Anderson,et al.  Preliminary reference earth model , 1981 .

[37]  B. Kennett,et al.  Source Depth and Mechanism Inversion at Teleseismic Distances Using a Neighborhood Algorithm , 2000 .

[38]  Alberto Malinverno,et al.  Expanded uncertainty quantification in inverse problems: Hierarchical Bayes and empirical Bayes , 2004 .

[39]  B. Romanowicz,et al.  Imaging anisotropic layering with Bayesian inversion of multiple data types , 2015 .

[40]  M. Sambridge,et al.  Transdimensional tomography with unknown data noise , 2012 .

[41]  Jan Dettmer,et al.  Trans-dimensional finite-fault inversion , 2014 .

[42]  Guust Nolet,et al.  Two-stage subduction history under North America inferred from multiple-frequency tomography , 2008 .

[43]  Hiroo Kanamori,et al.  Use of long-period surface waves for rapid determination of earthquake-source parameters , 1981 .

[44]  Simon C. Stähler,et al.  Instaseis: instant global seismograms based on a broadband waveform database , 2015 .

[45]  Eric Larose,et al.  Locating a small change in a multiple scattering environment , 2010 .

[46]  Hiroo Kanamori,et al.  Uncertainty estimations for seismic source inversions , 2011 .

[47]  K. Sigloch,et al.  Multifrequency measurements of core-diffracted P waves (Pdiff) for global waveform tomography , 2015 .

[48]  Krzysztof Podgórski,et al.  Multivariate generalized Laplace distribution and related random fields , 2013, J. Multivar. Anal..

[49]  Hrvoje Tkalcic,et al.  Point source moment tensor inversion through a Bayesian hierarchical model , 2016 .

[50]  William Menke,et al.  Using waveform similarity to constrain earthquake locations , 1999 .