Version 1 (modified by lnerger, 9 years ago) (diff)


A new filter: The Error Subspace Transform Kalman Filter - ESTKF

Recently, we have developed a new filter formulation, the Error Subspace Transform Kalman Filter. Details on this filter will be available in the paper "A unification of ensemble square-root filters" by L. Nerger, T. Janjic, J. Schroeter, and W. Hiller to appear in the Monthly Weather Review in 2012.

While studying the detailed relationship of the SEIK filter and the Ensemble Transform Kalman Filter (ETKF), we found a possibility to obtain the same minimum transformation of the ensemble (i.e. the distance between the forecast and analysis ensembles is minimal in the Frobenius norm) at a slightly lower computational cost. The new filter results from a modification of the SEIK filter to use consistent projections between the state space and the error subspace represented by the ensemble of model states. (The SEIK filter itself already projects onto the error subspace. However, this projection is not fully consistent and will lead to small differences in the analysis ensemble that depend on the order of the states on the ensemble matrix.) As the new formulation is very similar to the ETKF, but

In the paper, we also tested the use of the symmetric square root in the SEIK filter instead of the Cholesky decomposition that is commonly applied. The symmetric square root is commonly used in the ETKF and will also be the default for the ESTKF. For the SEIK filter, the symmetric square-root improved state estimates in experiments with the Lorenz-96 model, for the case that the ensemble transformation was deterministic.

The ESTKF, it's localized variant LESTKF, as well as a formulation the SEIK filter with symmetric square-root and explicit ensemble transformation will soon become available with the release of Version 1.8 of PDAF.