MTfit: Bayesian Source InversionΒΆ
Bayesian Moment Tensor Inversion Code by David J Pugh
MTfit is a Bayesian forward model inversion code for moment tensor and double-couple source inversion using different data types, based on the Bayesian approach presented in Pugh et al, 2016a and Pugh, 2015. The code has been developed as part of a PhD project (Pugh, 2015). The solutions are estimated using polarity and amplitude ratio data, although the code is extensible (see Extending MTfit) so it is possible to include other data-types in this framework. MTfit can incorporate uncertainty estimates both in the data (noise etc.) and the model (and location) in the resultant posterior PDF. There are three sampling approaches that have been developed, with different advantages (Pugh et al, 2015c, and it is also possible to use the approach for relative amplitude inversion as well (Pugh et al, 2015e).
MTfit also works with the automated Bayesian polarity approach described in Pugh et al, 2016b as an alternative method of estimating polarity probabilities. This is written as the autopol Python module, and may be available on request.
The source code is available from GitHub.
Issues are tracked at the GitHub repository. Please raise a new issue or feature request here.
This documentation is available as a PDF and an epub file
Please note that this code is provided as-is, and no guarantee is given that this code will perform in the desired way. Additional development and support is carried out in the developer's free time.
Restricted: For Non-Commercial Use Only This code is protected intellectual property and is available solely for teaching and non-commercially funded academic research purposes. Applications for commercial use should be made to Schlumberger or the University of Cambridge.