Multicore Parallel Tempering Bayeslands for Basin and Landscape Evolution. (16th November 2019)

• Record Type: Journal Article
• Authors: Chandra, Rohitash; Müller, R. Dietmar; Azam, Danial; Deo, Ratneel; Butterworth, Nathaniel; Salles, Tristan; Cripps, Sally
• Is Part Of: Geochemistry, geophysics, geosystems. Volume 20:Number 11(2019)
• Page Start: 5082
• Page End: 5104
• Publication Date: 2019-11-16
• Subjects: landscape evolution models -- Bayesian inference -- MCMC -- parallel computing -- parallel tempering; Geochemistry -- Periodicals; Geophysics -- Periodicals; Earth sciences -- Periodicals; 550.5
• DOI: 10.1029/2019GC008465
• Languages: English
• ISSNs: 1525-2027

Abstract

Abstract: The Bayesian paradigm is becoming an increasingly popular framework for estimation and uncertainty quantification of unknown parameters in geophysical inversion problems. Badlands is a landscape evolution model for simulating topography evolution at a broad range of spatial and temporal scales. Our previous work presented Bayeslands that used the Bayesian inference for estimating unknown parameters in the Badlands model using Markov chain Monte Carlo sampling. Bayeslands faced challenges in terms of computational issues and convergence due to multimodal posterior distributions. Parallel tempering is an advanced Markov chain Monte Carlo method suited for irregular and multimodal posterior distributions. In this paper, we extend Bayeslands using parallel tempering with high‐performance computing to address previous limitations in Bayeslands. Our results show that parallel tempering Bayeslands not only reduces the computation time‚ but also provides an improvement in sampling multimodal posterior distributions, which motivates future application to continental scale landscape evolution models. Key Points: We enhance the Bayeslands framework for estimation of free parameters in landscape evolution models We use high‐performance computing with parallel tempering MCMC for enhancing Bayeslands The enhanced Bayeslands framework not only reduces the computation time, but also provides better exploration of multimodal distributions