Journal cover Journal topic
Solid Earth An interactive open-access journal of the European Geosciences Union
Journal topic

Journal metrics

Journal metrics

  • IF value: 2.380 IF 2.380
  • IF 5-year value: 3.147 IF 5-year
    3.147
  • CiteScore value: 3.06 CiteScore
    3.06
  • SNIP value: 1.335 SNIP 1.335
  • IPP value: 2.81 IPP 2.81
  • SJR value: 0.779 SJR 0.779
  • Scimago H <br class='hide-on-tablet hide-on-mobile'>index value: 32 Scimago H
    index 32
  • h5-index value: 31 h5-index 31
Volume 5, issue 2
Solid Earth, 5, 1055–1069, 2014
https://doi.org/10.5194/se-5-1055-2014
© Author(s) 2014. This work is distributed under
the Creative Commons Attribution 3.0 License.
Solid Earth, 5, 1055–1069, 2014
https://doi.org/10.5194/se-5-1055-2014
© Author(s) 2014. This work is distributed under
the Creative Commons Attribution 3.0 License.

Method article 17 Nov 2014

Method article | 17 Nov 2014

Fully probabilistic seismic source inversion – Part 1: Efficient parameterisation

S. C. Stähler1 and K. Sigloch2,1 S. C. Stähler and K. Sigloch
  • 1Dept. of Earth and Environmental Sciences, Ludwig-Maximilians-Universität (LMU), Munich, Germany
  • 2Dept. of Earth Sciences, University of Oxford, South Parks Road, Oxford OX1 3AN, UK

Abstract. Seismic source inversion is a non-linear problem in seismology where not just the earthquake parameters themselves but also estimates of their uncertainties are of great practical importance. Probabilistic source inversion (Bayesian inference) is very adapted to this challenge, provided that the parameter space can be chosen small enough to make Bayesian sampling computationally feasible. We propose a framework for PRobabilistic Inference of Seismic source Mechanisms (PRISM) that parameterises and samples earthquake depth, moment tensor, and source time function efficiently by using information from previous non-Bayesian inversions. The source time function is expressed as a weighted sum of a small number of empirical orthogonal functions, which were derived from a catalogue of >1000 source time functions (STFs) by a principal component analysis. We use a likelihood model based on the cross-correlation misfit between observed and predicted waveforms. The resulting ensemble of solutions provides full uncertainty and covariance information for the source parameters, and permits propagating these source uncertainties into travel time estimates used for seismic tomography. The computational effort is such that routine, global estimation of earthquake mechanisms and source time functions from teleseismic broadband waveforms is feasible.

Publications Copernicus
Download
Citation