  We study the Electrical Impedance Tomography Bayesian inverse problem for
recovering the conductivity given noisy measurements of the voltage on some
boundary surface electrodes. The uncertain conductivity depends linearly on a
countable number of uniformly distributed random parameters in a compact
interval, with the coefficient functions in the linear expansion decaying at an
algebraic rate. We analyze the surrogate Markov Chain Monte Carlo (MCMC)
approach for sampling the posterior probability measure, where the multivariate
sparse adaptive interpolation, with interpolating points chosen according to a
lower index set, is used for approximating the forward map. The forward
equation is approximated once before running the MCMC for all the realizations,
using interpolation on the finite element (FE) approximation at the parametric
interpolating points. When evaluation of the solution is needed for a
realization, we only need to compute a polynomial, thus cutting drastically the
computation time. We contribute a rigorous error estimate for the MCMC
convergence. In particular, we show that there is a nested sequence of
interpolating lower index sets for which we can derive an interpolation error
estimate in terms of the cardinality of these sets, uniformly for all the
parameter realizations. An explicit convergence rate for the MCMC sampling of
the posterior expectation of the conductivity is rigorously derived, in terms
of the interpolating point number, the accuracy of the FE approximation of the
forward equation, and the MCMC sample number. We perform numerical experiments
using an adaptive greedy approach to construct the sets of interpolation
points. We show the benefits of this approach over the simple MCMC where the
forward equation is repeatedly solved for all the samples and the non-adaptive
surrogate MCMC with an isotropic index set treating all the random parameters
equally.
