An MCMC Bayesian full moment tensor inversion constrained by first-motion polarities and double couple percent

Monte Carlo Markov chain (MCMC) samplings can obtain a set of samples by directed random walk, mapping the posterior probability density of the model parameters in Bayesian framework. We perform earthquake waveform inversion to retrieve focal angles or the elements of moment tensor and source location using a Bayesian MCMC method with the constraints of first-motion polarities and double couple percentage using full Green functions and data covariance matrix. The algorithm tests the compatibility with polarities and also checks the double couple percentage of every site before the time-consuming synthetic seismogram computation for every sample of moment tensor of every trial source position. Other than large earthquakes, the method is especially suitable for weak events (M