RUS  ENG
Full version
JOURNALS // Numerical methods and programming // Archive

Num. Meth. Prog., 2018 Volume 19, Issue 4, Pages 507–515 (Mi vmp938)

Comparison of data assimilation methods based on the classical, ensemble and local Kalman filter by the example of the advection equation and Lorenz system

D. A. Rostilova, M. N. Kaurkinb, R. A. Ibrayevc

a Lomonosov Moscow State University, Faculty of Computational Mathematics and Cybernetics
b P. P. Shirshov Institute of Oceanology, Russian Academy of Sciences
c Institute of Numerical Mathematics of the Russian Academy of Sciences, Moscow

Abstract: The paper is devoted to the comparison of three data assimilation methods: the Kalman Filter (Kalman Filter, KF), the ensemble Kalman Filter (EnKF), and the local Kalman Filter (LKF). A number of numerical experiments on data assimilation by these methods are performed on two different models described by systems of differential equations. The first one is a simple one-dimensional linear equation of advection and the second one is the Lorenz system. The mean errors and the execution time of these assimilation methods are compared for different model sizes. The numerical results are consistent with the theoretical estimates. It is shown that the computational complexity of local and ensemble Kalman filters grows linearly with the size of the model, whereas in the classical Kalman Filter this complexity increases according to the cubic law. The efficiency of parallel implementation of the local Kalman filter is considered.

Keywords: Kalman Filter, Ensemble Kalman Filter, Local Kalman Filter, Lorenz equations, advection equation, data assimilation.

UDC: 004.942

Received: 23.06.2018



© Steklov Math. Inst. of RAS, 2024