Nonlinear stability and ergodicity of ensemble based Kalman filters

被引:67
作者
Tong, Xin T. [1 ,2 ]
Majda, Andrew J. [1 ,2 ]
Kelly, David [1 ]
机构
[1] NYU, Courant Inst Math Sci, 251 Mercer St, New York, NY 10012 USA
[2] NYU, Ctr Atmosphere Ocean Sci, 550 1St Ave, New York, NY 10012 USA
关键词
ensemble Kalman filter; filter stability; geometric ergodicity; Lyapunov function; catastrophic filter divergence; DISCRETE;
D O I
10.1088/0951-7715/29/2/657
中图分类号
O29 [应用数学];
学科分类号
070104 ;
摘要
The ensemble Kalman filter (EnKF) and ensemble square root filter (ESRF) are data assimilation methods used to combine high dimensional, nonlinear dynamical models with observed data. Despite their widespread usage in climate science and oil reservoir simulation, very little is known about the long-time behavior of these methods and why they are effective when applied with modest ensemble sizes in large dimensional turbulent dynamical systems. By following the basic principles of energy dissipation and controllability of filters, this paper establishes a simple, systematic and rigorous framework for the nonlinear analysis of EnKF and ESRF with arbitrary ensemble size, focusing on the dynamical properties of boundedness and geometric ergodicity. The time uniform boundedness guarantees that the filter estimate will not diverge to machine infinity in finite time, which is a potential threat for EnKF and ESQF known as the catastrophic filter divergence. Geometric ergodicity ensures in addition that the filter has a unique invariant measure and that initialization errors will dissipate exponentially in time. We establish these results by introducing a natural notion of observable energy dissipation. The time uniform bound is achieved through a simple Lyapunov function argument, this result applies to systems with complete observations and strong kinetic energy dissipation, but also to concrete examples with incomplete observations. With the Lyapunov function argument established, the geometric ergodicity is obtained by verifying the controllability of the filter processes; in particular, such analysis for ESQF relies on a careful multivariate perturbation analysis of the covariance eigen-structure.
引用
收藏
页码:657 / 691
页数:35
相关论文
共 43 条
[1]  
Anderson JL, 2001, MON WEATHER REV, V129, P2884, DOI 10.1175/1520-0493(2001)129<2884:AEAKFF>2.0.CO
[2]  
2
[3]  
Anderson JL, 1999, MON WEATHER REV, V127, P2741, DOI 10.1175/1520-0493(1999)127<2741:AMCIOT>2.0.CO
[4]  
2
[5]  
[Anonymous], 1972, STOCHASTIC PROCESSES
[6]  
[Anonymous], 2005, INFORM THEORY STOCHA
[7]  
[Anonymous], 2003, Ocean Dynamics, DOI [10.1007/s10236-003-0036-9, DOI 10.1007/S10236-003-0036-9]
[8]  
[Anonymous], 2003, ATMOSPHERIC MODELING
[9]  
[Anonymous], DISCRETE CONTINUOUS
[10]  
[Anonymous], 2000, ITO CALCULUS