TY - JOUR
T1 - Well-posedness and accuracy of the ensemble Kalman filter in discrete and continuous time
AU - Kelly, D. T B
AU - Law, Kody
AU - Stuart, Andrew M.
N1 - KAUST Repository Item: Exported on 2020-10-01
Acknowledgements: The authors are grateful to A J Majda for helpful discussions concerning this work. DTBK is supported by ONR grant N00014-12-1-0257. The work of AMS is supported by ERC, EPSRC, ESA and ONR. KJHL was supported by the King Abdullah University of Science and Technology (KAUST) and is a member of the KAUST Strategic Research Initiative Center for Uncertainty Quantification.
PY - 2014/9/22
Y1 - 2014/9/22
N2 - The ensemble Kalman filter (EnKF) is a method for combining a dynamical model with data in a sequential fashion. Despite its widespread use, there has been little analysis of its theoretical properties. Many of the algorithmic innovations associated with the filter, which are required to make a useable algorithm in practice, are derived in an ad hoc fashion. The aim of this paper is to initiate the development of a systematic analysis of the EnKF, in particular to do so for small ensemble size. The perspective is to view the method as a state estimator, and not as an algorithm which approximates the true filtering distribution. The perturbed observation version of the algorithm is studied, without and with variance inflation. Without variance inflation well-posedness of the filter is established; with variance inflation accuracy of the filter, with respect to the true signal underlying the data, is established. The algorithm is considered in discrete time, and also for a continuous time limit arising when observations are frequent and subject to large noise. The underlying dynamical model, and assumptions about it, is sufficiently general to include the Lorenz '63 and '96 models, together with the incompressible Navier-Stokes equation on a two-dimensional torus. The analysis is limited to the case of complete observation of the signal with additive white noise. Numerical results are presented for the Navier-Stokes equation on a two-dimensional torus for both complete and partial observations of the signal with additive white noise.
AB - The ensemble Kalman filter (EnKF) is a method for combining a dynamical model with data in a sequential fashion. Despite its widespread use, there has been little analysis of its theoretical properties. Many of the algorithmic innovations associated with the filter, which are required to make a useable algorithm in practice, are derived in an ad hoc fashion. The aim of this paper is to initiate the development of a systematic analysis of the EnKF, in particular to do so for small ensemble size. The perspective is to view the method as a state estimator, and not as an algorithm which approximates the true filtering distribution. The perturbed observation version of the algorithm is studied, without and with variance inflation. Without variance inflation well-posedness of the filter is established; with variance inflation accuracy of the filter, with respect to the true signal underlying the data, is established. The algorithm is considered in discrete time, and also for a continuous time limit arising when observations are frequent and subject to large noise. The underlying dynamical model, and assumptions about it, is sufficiently general to include the Lorenz '63 and '96 models, together with the incompressible Navier-Stokes equation on a two-dimensional torus. The analysis is limited to the case of complete observation of the signal with additive white noise. Numerical results are presented for the Navier-Stokes equation on a two-dimensional torus for both complete and partial observations of the signal with additive white noise.
UR - http://hdl.handle.net/10754/563763
UR - https://iopscience.iop.org/article/10.1088/0951-7715/27/10/2579
UR - http://www.scopus.com/inward/record.url?scp=84907246706&partnerID=8YFLogxK
U2 - 10.1088/0951-7715/27/10/2579
DO - 10.1088/0951-7715/27/10/2579
M3 - Article
SN - 0951-7715
VL - 27
SP - 2579
EP - 2603
JO - Nonlinearity
JF - Nonlinearity
IS - 10
ER -