The Error-Subspace Transform Kalman Filter
Ensemble square-root Kalman filters are currently the computationally most efficient ensemble-based Kalman filter methods. In particular, the Ensemble Transform Kalman Filter (ETKF) is known to provide a minimum ensemble transformation in a very efficient way. A similar filter algorithm is the Singular Evolutive Interpolated Kalman (SEIK) filter. In contrast to the ETKF, the SEIK filter solves the estimation problem of the Kalman filter directly in the in the error-subspace that is represented by the ensemble. Further, the matrix square roots computed in the analysis step of the filters are distinct. Comparing both filters shows small differences between their ensemble transformations, while their state estimates are identical if the same forecast ensemble is used for the analysis. To obtain the minimum transformation also with a filter using directly the error subspace, a new projection onto the error subspace is introduced. It leads to the Error-Subspace Transform Kalman Filter (ESTKF). The ESTKF shows a slightly lower computational cost compared to the ETKF. Using numerical experiments with the parallel data assimilation system PDAF [http://pdaf.awi.de], we show that also the SEIK filter can profit from the use of a symmetric matrix square root. In addition, the effect of the new projection in the ESTKF is shown.
AWI Organizations > Infrastructure > Scientific Computing