Transdimensional Markov Chain Monte Carlo joint inversion of direct current resistivity and transient electromagnetic data

Abstract

SUMMARYJoint inversion of multiple geophysical data sets with complementary information content can significantly reduce the non-uniqueness inherent to each individual data set and, therefore, can improve subsurface characterization. Gradient-based joint inversion methods depend on the choice of model regularization and usually produce one single optimal model, and rely on linearization to estimate model parameter uncertainty. However, a quantitative evaluation of the parameter uncertainty of the derived model parameters is crucial for reliable data interpretation. In this study, we present a transdimensional Markov Chain Monte Carlo (MCMC) method for the joint inversion of direct current resistivity and transient electromagnetic data, which provides a rigorous assessment of the uncertainty associated with the derived model. The transdimensional property of the algorithm allows the number of unknown model parameters to be determined adaptively by the data. This usually favours models with fewer parameters through the parsimony criterion of the Bayesian method by choosing suitable prior distributions. In this paper, we demonstrate that the transdimensional MCMC method combines complementary information contained in each data set and reduces the overall uncertainty using synthetic examples. Furthermore, we successfully applied the new joint inversion scheme to field data from Azraq, Jordan. The transdimensional MCMC inversion results are in good agreement with the results obtained by deterministic inversion techniques. From the MCMC inversion results we identified the thickness of a basalt formation and a conductive zone, which were uncertain and not interpreted in prior studies, adding to the geological interpretation.