Bayesian approach to magnetotelluric tensor decomposition
Main Article Content
Abstract
Magnetotelluric directional analysis and impedance tensor decomposition are basic tools to validate a local/regional composite electrical model of the underlying structure. Bayesian stochastic methods approach the problem of the parameter estimation and their uncertainty characterization in a fully probabilistic fashion, through the use of posterior model probabilities.We use the standard GroomBailey 3D local/2D regional composite model in our bayesian approach. We assume that the experimental impedance estimates are contamined with the Gaussian noise and define the likelihood of a particular composite model with respect to the observed data. We use noninformative, flat priors over physically reasonable intervals for the standard GroomBailey decomposition parameters. We apply two numerical methods, the Markov chain Monte Carlo procedure based on the Gibbs sampler and a singlecomponent adaptive Metropolis algorithm. From the posterior samples, we characterize the estimates and uncertainties of the individual decomposition parameters by using the respective marginal posterior probabilities. We conclude that the stochastic scheme performs reliably for a variety of models, including the multisite and multifrequency case with up to several hundreds of parameters. Though the Monte Carlo samplers are computationally very intensive, the adaptive Metropolis algorithm increase the speed of the simulations for largescale problems.
Article Details
Open-Access License
No Permission Required
Istituto Nazionale di Geofisica e Vulcanologia applies the Creative Commons Attribution License (CCAL) to all works we publish.
Under the CCAL, authors retain ownership of the copyright for their article, but authors allow anyone to download, reuse, reprint, modify, distribute, so long as the original authors and source are cited. No permission is required from the authors or the publishers.