Parallel tempering for strongly nonlinear geoacoustic inversion

Stan E. Dosso, Charles W. Holland, Malcolm Sambridge

Research output: Contribution to journalArticle

51 Scopus citations

Abstract

This paper applies parallel tempering within a Bayesian formulation for strongly nonlinear geoacoustic inverse problems. Bayesian geoacoustic inversion consists of sampling the posterior probability density (PPD) of seabed parameters to estimate integral properties, such as marginal probability distributions, based on ocean acoustic data and prior information. This sampling is usually carried out using the Markov-chain Monte Carlo method of Metropolis-Hastings sampling. However, standard sampling methods can be very inefficient for strongly nonlinear problems involving multi-modal PPDs with the potential to miss important regions of the parameter space and to significantly underestimate parameter uncertainties. Parallel tempering achieves efficient/effective sampling of challenging parameter spaces with the ability to transition freely between multiple PPD modes by running parallel Markov chains at a series of increasing sampling temperatures with probabilistic interchanges between chains. The approach is illustrated for inversion of (simulated) acoustic reverberation data for which the PPD is highly multi-modal. While Metropolis-Hastings sampling gives poor results even with very large sample sizes, parallel tempering provides efficient, convergent sampling of the PPD. Methods to enhance the efficiency of parallel tempering are also considered.

Original languageEnglish (US)
Pages (from-to)3030-3040
Number of pages11
JournalJournal of the Acoustical Society of America
Volume132
Issue number5
DOIs
StatePublished - Nov 1 2012

    Fingerprint

All Science Journal Classification (ASJC) codes

  • Arts and Humanities (miscellaneous)
  • Acoustics and Ultrasonics

Cite this