This is software for inversion of seismic Receiver Functions for 1-D Shear wave velocity structure using Reversible jump Markov chain Monte Carlo algorithm. Full details appear in Bodin et al. (2012). In this algorithm the number of layers in the model is unknown as are the interface depths between layers and the Vs wave velocity in each layer. The algorithm results in an ensemble of solutions and properties of the ensemble may be generated for inspection. The point wise spatial mean of the ensemble of Earth models is used as a representative image of the solution. The variance of the ensemble of models at each point in the medium provides information on uncertainty.

Outputs include an ensemble averaged Vs velocity profile and probability distributions on the position of interfaces between layers.


The Fortran 90 source code can be downloaded here. You will need to register with iEarth prior to download. It has been run successfully on unix workstations. Instructions for use and some examples are included. For examples of applications see the paper below. Enquires should be directed to the author.

The contents of this package are provided under the terms of the GPL licence


The software requires the rj-MCMC library and if this prerequisite is in a standard place then you can compile this software with the following commands:

    > ./configure
    > make
and optionally:
    > sudo make install

The Rf software uses pkg-config to configure itself while building. If the rjmcmc library is installed into a non-standard directory (eg a users home directory) or to a path that is not by default in pkg-config's search path then you will need to set the PKG_CONFIG_PATH environment variable, eg for bash users:

    > export PKG_CONFIG_PATH=/opt/rjmcmc/lib/pkg-config
    > ./configure
    > make
It is recommended that new users run through this tutorial as a starting point.


    Transdimensional inversion of receiver functions and surface wave dispersion Bodin, T., Sambridge, M., Tkalčić, H., Arroucau, P., Gallagher, K., and Rawlinson, N., J. Geophys. Res., 117, B02301. 2012. doi:10.1029/2011JB008560.

    "Crustal and Uppermost Mantle structure variation beneath La Réunion hotspot track", Fontaine, F.R., Barruol, G., H. Tkalčić, Wölbern, I., Rumpker, G., Bodin, T. and Haugmard, M., Geophys. J. Int., 203, 107-126, 2015.