Joint gravity and magnetic inversion in 3D using Monte Carlo methods

We jointly invert gravity and magnetic data following a Monte Carlo method that provides estimation for a 3D model of the structure and physical properties of the medium. In particular, the model layer geometry and the density and magnetic susceptibility fields within layers are estimated, and their uncertainties are described with posterior probabilities. This method combines the gravity and magnetic data with prior information of the mass density and magnetic susceptibility statistics, and statistical constraints on the model interface positions. The resulting model realizations jointly comply with the observations and the prior statistical information.