  Bayesian modelling and computational inference by Markov chain Monte Carlo
(MCMC) is a principled framework for large-scale uncertainty quantification,
though is limited in practice by computational cost when implemented in the
simplest form that requires simulating an accurate computer model at each
iteration of the MCMC. The delayed acceptance Metropolis--Hastings MCMC
leverages a reduced model for the forward map to lower the compute cost per
iteration, though necessarily reduces statistical efficiency that can, without
care, lead to no reduction in the computational cost of computing estimates to
a desired accuracy. Randomizing the reduced model for the forward map can
dramatically improve computational efficiency, by maintaining the low cost per
iteration but also avoiding appreciable loss of statistical efficiency.
Randomized maps are constructed by a posteriori adaptive tuning of a randomized
and locally-corrected deterministic reduced model. Equivalently, the
approximated posterior distribution may be viewed as induced by a modified
likelihood function for use with the reduced map, with parameters tuned to
optimize the quality of the approximation to the correct posterior
distribution. Conditions for adaptive MCMC algorithms allow practical
approximations and algorithms that have guaranteed ergodicity for the target
distribution. Good statistical and computational efficiencies are demonstrated
in examples of calibration of large-scale numerical models of geothermal
reservoirs and electrical capacitance tomography.
