Journal cover Journal topic
Nonlinear Processes in Geophysics An interactive open-access journal of the European Geosciences Union
Journal topic

Journal metrics

Journal metrics

  • IF value: 1.699 IF 1.699
  • IF 5-year value: 1.559 IF 5-year
    1.559
  • CiteScore value: 1.61 CiteScore
    1.61
  • SNIP value: 0.884 SNIP 0.884
  • IPP value: 1.49 IPP 1.49
  • SJR value: 0.648 SJR 0.648
  • Scimago H <br class='hide-on-tablet hide-on-mobile'>index value: 52 Scimago H
    index 52
  • h5-index value: 21 h5-index 21
Volume 21, issue 5
Nonlin. Processes Geophys., 21, 955–970, 2014
https://doi.org/10.5194/npg-21-955-2014
© Author(s) 2014. This work is distributed under
the Creative Commons Attribution 3.0 License.
Nonlin. Processes Geophys., 21, 955–970, 2014
https://doi.org/10.5194/npg-21-955-2014
© Author(s) 2014. This work is distributed under
the Creative Commons Attribution 3.0 License.

Research article 23 Sep 2014

Research article | 23 Sep 2014

Improving the ensemble transform Kalman filter using a second-order Taylor approximation of the nonlinear observation operator

G. Wu1, X. Yi1, L. Wang2, X. Liang3, S. Zhang1, X. Zhang1, and X. Zheng1 G. Wu et al.
  • 1State Key Laboratory of Remote Sensing Science, College of Global Change and Earth System Science, Beijing Normal University, Beijing, China
  • 2Department of Statistics, University of Manitoba, Winnipeg, Canada
  • 3National Meteorological Information Center, China Meteorological Administration, Beijing, China

Abstract. The ensemble transform Kalman filter (ETKF) assimilation scheme has recently seen rapid development and wide application. As a specific implementation of the ensemble Kalman filter (EnKF), the ETKF is computationally more efficient than the conventional EnKF. However, the current implementation of the ETKF still has some limitations when the observation operator is strongly nonlinear. One problem in the minimization of a nonlinear objective function similar to 4D-Var is that the nonlinear operator and its tangent-linear operator have to be calculated iteratively if the Hessian is not preconditioned or if the Hessian has to be calculated several times. This may be computationally expensive. Another problem is that it uses the tangent-linear approximation of the observation operator to estimate the multiplicative inflation factor of the forecast errors, which may not be sufficiently accurate.

This study attempts to solve these problems. First, we apply the second-order Taylor approximation to the nonlinear observation operator in which the operator, its tangent-linear operator and Hessian are calculated only once. The related computational cost is also discussed. Second, we propose a scheme to estimate the inflation factor when the observation operator is strongly nonlinear. Experimentation with the Lorenz 96 model shows that using the second-order Taylor approximation of the nonlinear observation operator leads to a reduction in the analysis error compared with the traditional linear approximation method. Furthermore, the proposed inflation scheme leads to a reduction in the analysis error compared with the procedure using the traditional inflation scheme.

Publications Copernicus
Download
Citation