Bayesian estimation of non-linear centroid moment tensors using multiple seismic data sets

Abstract

SUMMARY Centroid moment tensor (CMT) parameters of earthquakes are routinely estimated to gain information on structures and regional tectonics. However, for small earthquakes (M < 4), it is still challenging to determine CMTs due to the lack of high-quality waveform data. In this study, we propose to improve solutions for small earthquakes by incorporating multiple seismic data types in Bayesian joint inversion: polarities picked on broad-band signals, amplitude spectra for intermediate frequency bands (0.2–2.0 Hz), and waveforms at low frequencies (0.05–0.2 Hz). Both measurement and theory errors are accounted for by iterative estimation of non-Toeplitz covariance matrices, providing objective weightings for the different data types in the joint parameter estimation. Validity and applicability of the method are demonstrated using simulated and field data. Results demonstrate that combinations of data, such as a single high-quality waveform, a few amplitude spectra and many waveform polarities, are able to resolve CMT parameters to comparable quality as if many high-quality waveforms were available. Results of 10 induced seismic events that occurred in northeastern British Columbia, Canada, between January 2020 and February 2022 indicate predominantly strike-slip focal mechanisms with low non-double-couple components. These events appear to be located at shallow depths with short time duration, as expected for induced seismicity. These results are consistent with previous studies, indicating that this method reduces the dependence of source inversion on high-quality waveforms, and can provide resolution of CMT parameters for earthquakes as small as ML 1.6.