Moment Tensor Inversion
The functionality of this module is based on the Bayes-ISOLA code (Vackář et al., 2017), a tool for automated CMT inversion in a Bayesian framework. This module a function to calculate an optional data covariance matrix from pre-event noise that yields an automatic weighting function for the receiver components according to their noise levels. This also serves as an automated frequency filter that suppresses noisy frequency ranges. During the inversion, a grid search is performed over time and space, combined with an analytical (least squares) moment tensor inversion at each grid point, which accelerates the inversion process. The more time-consuming tasks like the Green’s function computation and the grid search are parallelized for efficiency. The results provided by the module contain the best-fit solution, as well as the full posterior probability density function (PPDF), which allows the user to plot the marginal probability density functions (PDFs) for any of the CMT parameters.
Create Earth velocity model
To start you need to load a Earth model. Go to Build a 1D Earth Model and open the Earth Model Form, the fill it and save it.
Alternatively, Yo can find examples of regional models in ./isp/mti/input/iberia.dat and ./isp/mti/input/AK135.txt.
Bayesian Inversion
The user needs to fill the following fields to go on with the inversion:
Metadata >> Path to the file where you have the metadata. This is necessary to extract the stations coordinates and the instrument information to do the deconvolution (this is not totally necessary, the user can manage the seismograms already in velocity units).
-
Files >> path to the seismogram files
-
Earth Model >>> Path to the Earth model file (see previous section for see how to make this file with ISP, the user can also do it manually).
-
Time Window >>> Set the window in which the seismograms will be used (StartTime, EndTime and CHECK trim).
-
Macro >> You can set a Macro to make basic processing over the seismograms
Now Follow this steps,
Process and plot Seismograms Open Stations info. This action will allow you select which station/component will be used in the inversion
Plot stations Map (First, fill Latitude Longitude and depth) This is optional - Fill the info: - Hipocenter location and Uncertainity (m) - Origin Time and Uncertainity - Magnitude - Frequency range for the inversion (see Fig 2 for example) - Distance filter. Avoid statations out of the distance filter (Distance from epicenter to station) - Option Deviatoric, Circle, Shape, Covariance (please because it uses previous noise window for inversion, select long enough prevent window (i.e., 10 min)) and use precalculated Green functions
Pre-calculated Green Functions is useful because, once you have run an inversion then you may filter the inversion selecting different components or distances
