Efficient hydrological model parameter optimization with Sequential Monte Carlo sampling

  • Authors:

  • Erwin Jeremiah;Scott A. Sisson;Ashish Sharma;Lucy Marshall

  • Affiliations:

  • School of Civil and Environmental Engineering, University of New South Wales, Sydney, NSW 2052, Australia;School of Mathematics and Statistics, University of New South Wales, Sydney, NSW, Australia;School of Civil and Environmental Engineering, University of New South Wales, Sydney, NSW 2052, Australia;Land Resources and Environmental Science, Montana State University, Bozeman, Montana, United States

  • Venue:
  • Environmental Modelling & Software
  • Year:
  • 2012

Quantified Score

Hi-index 0.00

Visualization

Abstract

Bayesian inference provides an ideal platform for assessing parameter uncertainty for complex physical models, such as conceptual hydrological models. Sequential Monte Carlo (SMC) samplers are well suited for its implementation as they are effective in sampling from posterior distributions with the non-linear dependency structures and multiple modes often present in hydrological models. A challenge in implementing SMC samplers is in the construction of a suitable sequence of intermediary distributions leading to the posterior distribution, in such a way that the sampler performs both efficiently in that minimal computation is needed, and robustly in that the samplers particle population is effectively maintained. In this article we demonstrate that naive implementation of an SMC sampler in a hydrological model can cause the sampler to collapse. To address this, we propose a new method of dynamically constructing the transition path between the intermediary distributions that effectively increases the amount of computation at the point where it is needed to avoid this collapse. We analyse the performance of our approach through the analysis of real and simulated hydrological data, and discuss related SMC sampler implementation issues. In addition we consider the question of the appropriate number of particles in the sampler as the dimensionality of the model increases.