Investigating the effects of random data errors on the waveform-based moment tensor inversion
Open access
Date
2022-04Type
- Journal Article
ETH Bibliography
yes
Altmetrics
Abstract
The linear Gauss-Markov model for waveform-based moment tensor inversion often relies on the overdetermined least-squares method. It needs a proper stochastic model of the observables for accurate and precise estimates of the unknown parameters. Furthermore, estimating the level and distribution of random errors in the observed waveforms is challenging due to assessing the minimum-variance unbiased estimator (MVUE). Hence, according to the considerable effects of random data errors in assessing the uncertainty of the moment tensor components, this paper aims to describe an MVUE of the data covariance matrix and its application on uncertainty quantification of the moment tensor. The used mathematical prescription allows us to use the covariance matrix for the three-component noise records at every station and all possible cross-correlations among the recorded noise wavefield. To illustrate the proposed method's performance, we conducted tests with synthetic data using configuration of the 2018 M-w 6.8 Zakynthos (Ionian Sea, Greece) earthquake. Both uncorrelated and correlated random noise traces were added to the synthetic waveform data in amounts between 5 and 20 per cent of the maximum amplitude. In order to test the efficiency of the method, we considered three different structures of covariance matrix: (i) diagonal matrix (contains a variance of individual measurements at seismic stations), (ii) block-diagonal matrix (considering cross-covariance among three components at each station), and (iii) full covariance matrix. Test results are presented by comparison of the moment tensor inversion outcomes with known noise levels of generated synthetic data and with synthetic focal mechanisms, the ability of the estimated full covariance matrix in illustrating the minimum variance of parameters (namely, minimum posterior uncertainties), unbiased of the parameters, and values of the cross-correlations between the components of each station and also among stations. Finally, we applied the method to the real waveforms of the Zakynthos earthquake having inferred focal mechanism of strike/dip/rake angles 13/40/171 (deg) with 33 per cent double couple (DC) and -61 per cent compensated linear vector dipole component (CLVD). The focal mechanism solution has strike/dip/rake angles 19/34/177 (deg) with 69 per cent DC and -23 per cent CLVD when using our estimated full covariance matrix. Show more
Permanent link
https://doi.org/10.3929/ethz-b-000530593Publication status
publishedExternal links
Journal / series
Geophysical Journal InternationalVolume
Pages / Article No.
Publisher
Oxford University PressSubject
Statistical methods; Time-series analysis; Waveform inversion; Statistical seismologyOrganisational unit
02818 - Schweiz. Erdbebendienst (SED) / Swiss Seismological Service (SED)
More
Show all metadata
ETH Bibliography
yes
Altmetrics