TY - JOUR
T1 - A Particle Filter-based Adaptive Inflation Scheme for the Ensemble Kalman Filter
AU - Ait-El-Fquih, Boujemaa
AU - Hoteit, Ibrahim
N1 - KAUST Repository Item: Exported on 2020-04-23
Acknowledgements: Research reported in this publication was supported by KingAbdullah University of Science and Technology (KAUST). The authors would like to thank the anonymous reviewers and Franc ̧oisDesbouvries for the fruitful discussions about the resampling in the particle filter.
PY - 2020/1/24
Y1 - 2020/1/24
N2 - An adaptive covariance inflation scheme is proposed for the ensemble Kalman filter (EnKF) to mitigate for the loss of ensemble variance. Adaptive inflation methods are mostly based on a Bayesian approach, which considers the inflation factor as a random variable with a given prior probability distribution, and then combines it with the inflation likelihood through Bayes’ rule to obtain its posterior distribution. In this work, we introduce a numerical implementation of this generic Bayesian approach that uses a particle filter (PF) to compute a Monte Carlo approximation of the inflation posterior distribution. To alleviate the sample attrition issue, the proposed PF employs an artificial dynamical model for the inflation factor based on the well-known smoothing-kernel West and Liu model. The positivity constraint on the inflation factor is further imposed through an inverse-Gamma transition density, whose parameters suggest analytical expressions. The resulting PF-EnKF scheme is straightforward to implement, and can use different number of particles in its EnKF and PF components. Numerical experiments are conducted with the Lorenz-96 model to demonstrate the effectiveness of the proposed method under various experimental scenarios.
AB - An adaptive covariance inflation scheme is proposed for the ensemble Kalman filter (EnKF) to mitigate for the loss of ensemble variance. Adaptive inflation methods are mostly based on a Bayesian approach, which considers the inflation factor as a random variable with a given prior probability distribution, and then combines it with the inflation likelihood through Bayes’ rule to obtain its posterior distribution. In this work, we introduce a numerical implementation of this generic Bayesian approach that uses a particle filter (PF) to compute a Monte Carlo approximation of the inflation posterior distribution. To alleviate the sample attrition issue, the proposed PF employs an artificial dynamical model for the inflation factor based on the well-known smoothing-kernel West and Liu model. The positivity constraint on the inflation factor is further imposed through an inverse-Gamma transition density, whose parameters suggest analytical expressions. The resulting PF-EnKF scheme is straightforward to implement, and can use different number of particles in its EnKF and PF components. Numerical experiments are conducted with the Lorenz-96 model to demonstrate the effectiveness of the proposed method under various experimental scenarios.
UR - http://hdl.handle.net/10754/660550
UR - https://onlinelibrary.wiley.com/doi/abs/10.1002/qj.3716
UR - http://www.scopus.com/inward/record.url?scp=85078664126&partnerID=8YFLogxK
U2 - 10.1002/qj.3716
DO - 10.1002/qj.3716
M3 - Article
SN - 0035-9009
JO - Quarterly Journal of the Royal Meteorological Society
JF - Quarterly Journal of the Royal Meteorological Society
ER -