Particle filtering with path sampling and an application to a bimodal ocean current model

This paper introduces a recursive particle filtering algorithm designed to filter high dimensional systems with complicated non-linear and non-Gaussian effects. The method incorporates a parallel marginalization (PMMC) step in conjunction with the hybrid Monte Carlo (HMC) scheme to improve samples generated by standard particle filters. Parallel marginalization is an efficient Markov chain Monte Carlo (MCMC) strategy that uses lower dimensional approximate marginal distributions of the target distribution to accelerate equilibration. As a validation the algorithm is tested on a 2516 dimensional, bimodal, stochastic model motivated by the Kuroshio current that runs along the Japanese coast. The results of this test indicate that the method is an attractive alternative for problems that require the generality of a particle filter but have been inaccessible due to the limitations of standard particle filtering strategies.

[1]  Paul Krause,et al.  Dimensional reduction for a Bayesian filter. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[2]  R. E. Kalman,et al.  A New Approach to Linear Filtering and Prediction Problems , 2002 .

[3]  Shenn-Yu Chao,et al.  Bimodality of the Kuroshio , 1984 .

[4]  G. Kitagawa.  Monte Carlo Filter and Smoother for Non-Gaussian Nonlinear State Space Models , 1996 .

[5]  Goodman,et al.  Multigrid Monte Carlo method for lattice field theories. , 1986, Physical review letters.

[6]  Alexandre J. Chorin.  Conditional Expectations and Renormalization , 2003, Multiscale Model. Simul..

[7]  G. Eyink,et al.  Accelerated Monte Carlo for Optimal Estimation of Time Series , 2005 .

[8]  Jun S. Liu,et al.  Sequential Monte Carlo methods for dynamic systems , 1997 .

[9]  Andrew M. Stuart,et al.  A Bayesian approach to Lagrangian data assimilation , 2008 .

[10]  Jonathan Weare,et al.  Parallel marginalization Monte Carlo with applications to conditional path sampling , 2007, 0709.1721.

[11]  Robert N. Miller,et al.  Data assimilation into nonlinear stochastic models , 1999 .

[12]  J. Banavar,et al.  Computer Simulation of Liquids , 1988 .

[13]  Jun S. Liu,et al.  Monte Carlo strategies in scientific computing , 2001 .

[14]  W. Gilks,et al.  Following a moving target—Monte Carlo inference for dynamic Bayesian models , 2001 .

[15]  Berend Smit,et al.  Understanding Molecular Simulation , 2001 .

[16]  Geir Evensen,et al.  The Ensemble Kalman Filter: theoretical formulation and practical implementation , 2003 .

[17]  Nando de Freitas,et al.  Sequential Monte Carlo Methods in Practice , 2001, Statistics for Engineering and Information Science.

[18]  S. Duane,et al.  Hybrid Monte Carlo , 1987 .

[19]  Jonathan Weare,et al.  Efficient Monte Carlo sampling by parallel marginalization , 2007, Proceedings of the National Academy of Sciences.

[20]  Bo Qiu,et al.  Kuroshio Path Variations South of Japan: Bimodality as a Self-Sustained Internal Oscillation , 2000 .

[21]  A. Stuart,et al.  Conditional Path Sampling of SDEs and the Langevin MCMC Method , 2004 .

[22]  Berend Smit,et al.  Understanding molecular simulation: from algorithms to applications , 1996 .