Bayeslands Workshop
Channel: Rohitash Chandra
Date: 2018-10-20

https://github.com/badlands-model/Bayeslands_Workshop

Bayesian inference provides a principled approach towards uncertainty quantification of free parameters in geophysical forward models. This provides advantages over optimization methods that provide single point estimates as solutions, which lack uncertainty quantification. Badlands (basin and landscape dynamics model) is geophysical forward model that simulates topography development at various space and time scales. Badlands consists of a number of geophysical parameters that need to be estimated with appropriate uncertainty quantification, given the observed ground truth such as surface topography, sediment thickness and stratigraphy through time. This is challenging due to the scarcity of data, sensitivity of the parameters and complexity of the Badlands model. In this paper, we take a Bayesian approach to provide inference using Markov chain Monte Carlo sampling (MCMC). Hence, we present \textit{BayesLands}, a Bayesian framework for Badlands that fuses information obtained from complex forward models with observational data and prior knowledge. As a proof-of-concept, we consider a synthetic and real-world topography with two free parameters, namely precipitation and erodibility, that we need to estimate through BayesLands. The results of the experiments shows that BayesLands yields a promising distribution of the parameters. Moreover, the challenge in sampling due to multi-modality is presented through visualizing a likelihood surface that has a range of suboptimal modes. https://arxiv.org/abs/1805.03696

In recent years, Bayesian inference has become a popular methodology for the estimation and uncertainty quantification of parameters in geological and geophysical forward models. Badlands is a basin and landscape evolution forward model for simulating topography evolution at a large range of spatial and time scales. Previously, Bayesian inference has been used for parameter estimation and uncertainty quantification in Badlands, an extension known as Bayeslands. It was demonstrated that the posterior surface of these parameters could exhibit highly irregular features such as multi-modality and discontinuities making standard Markov Chain Monte Carlo (MCMC) sampling difficult. Parallel tempering (PT) is an advanced MCMC method suited for irregular and multi-modal distributions. Moreover, PT is more suitable for multi-core implementations that can speed up computationally expensive geophysical models. In this paper, we present a multi-core PT algorithm implemented in a high performance computing architecture for enhancing Bayeslands. The results show that PT in Bayeslands not only reduces the computation time over a multi-core architecture, but also provides a means to improve the sampling process in a multi-modal landscape. This motivates its usage in larger-scale problems in basin and landscape evolution models.

https://arxiv.org/abs/1806.10939