Comparison and combination of EAKF and SIR-PF in the Bayesian filter framework  

作　　者：SHEN Zheqi ZHANG Xiangming TANG Youmin 

机构地区：[1]State Key Laboratory of Satellite Ocean Environment Dynamics, Second Institute of Oceanography,State Oceanic Administration [2]Environmental Science and Engineering, University of Northern British Columbia

出　　处：《Acta Oceanologica Sinica》2016年第3期69-78,共10页海洋学报（英文版）

基　　金：The National Natural Science Foundation of China under contract Nos 41276029 and 41321004;the Project of State Key Laboratory of Satellite Ocean Environment Dynamics,Second Institute of Oceanography under contract No.SOEDZZ1404;the National Basic Research Program(973 Program)of China under contract No.2013CB430302

摘　　要：Bayesian estimation theory provides a general approach for the state estimate of linear or nonlinear and Gaussian or non-Gaussian systems. In this study, we first explore two Bayesian-based methods： ensemble adjustment Kalman filter（EAKF） and sequential importance resampling particle filter（SIR-PF）, using a well-known nonlinear and non-Gaussian model（Lorenz ＇63 model）. The EAKF, which is a deterministic scheme of the ensemble Kalman filter（En KF）, performs better than the classical（stochastic） En KF in a general framework. Comparison between the SIR-PF and the EAKF reveals that the former outperforms the latter if ensemble size is so large that can avoid the filter degeneracy, and vice versa. The impact of the probability density functions and effective ensemble sizes on assimilation performances are also explored. On the basis of comparisons between the SIR-PF and the EAKF, a mixture filter, called ensemble adjustment Kalman particle filter（EAKPF）, is proposed to combine their both merits. Similar to the ensemble Kalman particle filter, which combines the stochastic En KF and SIR-PF analysis schemes with a tuning parameter, the new mixture filter essentially provides a continuous interpolation between the EAKF and SIR-PF. The same Lorenz ＇63 model is used as a testbed, showing that the EAKPF is able to overcome filter degeneracy while maintaining the non-Gaussian nature, and performs better than the EAKF given limited ensemble size.Bayesian estimation theory provides a general approach for the state estimate of linear or nonlinear and Gaussian or non-Gaussian systems. In this study, we first explore two Bayesian-based methods： ensemble adjustment Kalman filter（EAKF） and sequential importance resampling particle filter（SIR-PF）, using a well-known nonlinear and non-Gaussian model（Lorenz ＇63 model）. The EAKF, which is a deterministic scheme of the ensemble Kalman filter（En KF）, performs better than the classical（stochastic） En KF in a general framework. Comparison between the SIR-PF and the EAKF reveals that the former outperforms the latter if ensemble size is so large that can avoid the filter degeneracy, and vice versa. The impact of the probability density functions and effective ensemble sizes on assimilation performances are also explored. On the basis of comparisons between the SIR-PF and the EAKF, a mixture filter, called ensemble adjustment Kalman particle filter（EAKPF）, is proposed to combine their both merits. Similar to the ensemble Kalman particle filter, which combines the stochastic En KF and SIR-PF analysis schemes with a tuning parameter, the new mixture filter essentially provides a continuous interpolation between the EAKF and SIR-PF. The same Lorenz ＇63 model is used as a testbed, showing that the EAKPF is able to overcome filter degeneracy while maintaining the non-Gaussian nature, and performs better than the EAKF given limited ensemble size.

关 键 词：data assimilation ensemble adjustment Kalman filter particle filter Bayesian estimation ensemble adjustment Kalman particle filter 

分 类 号：P714[天文地球—海洋科学]