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

Journal metrics

  • IF value: 1.329 IF 1.329
  • IF 5-year<br/> value: 1.394 IF 5-year
  • CiteScore<br/> value: 1.27 CiteScore
  • SNIP value: 0.903 SNIP 0.903
  • SJR value: 0.709 SJR 0.709
  • IPP value: 1.455 IPP 1.455
  • h5-index value: 20 h5-index 20
© Author(s) 2017. This work is distributed under
the Creative Commons Attribution 4.0 License.
Research article
16 Nov 2017
Review status
This discussion paper is a preprint. A revision of this manuscript was accepted for the journal Nonlinear Processes in Geophysics (NPG) and is expected to appear here in due course.
Quasi static ensemble variational data assimilation
Anthony Fillion1,2, Marc Bocquet1, and Serge Gratton3 1CEREA, Joint Laboratory École des Ponts ParisTech and EDF R&D, Université Paris-Est, Champs-sur-Marne, France
2CERFACS, Toulouse, France
3INPT-IRIT, Toulouse, France
Abstract. The analysis in nonlinear variational data assimilation is the solution of a non-quadratic minimization. Thus, the analysis efficiency relies on its ability to locate a global minimum of the cost function. If this minimization uses a Gauss-Newton (GN) method, it is critical for the starting point to be in the attraction basin of a global minimum. Otherwise the method may converge to a local extremum, which degrades the analysis. With chaotic models, the number of local extrema often increases with the temporal extent of the data assimilation window, making the former condition harder to satisfy. This is unfortunate because the assimilation performance also increases with this temporal extent. However, a quasi-static (QS) minimization may overcome these local extrema. It consists in gradually injecting the observations in the cost function. This method was introduced by Pires et al. (1996) in a 4D-Var context.

We generalize this approach to four-dimensional nonlinear EnVar methods, which are based on both a nonlinear variational analysis and the propagation of dynamical error statistics via an ensemble. This forces to consider the cost function minimizations in the broader context of cycled data assimilation algorithms. We adapt this QS approach to the iterative ensemble Kalman smoother (IEnKS), an exemplar of nonlinear deterministic 4D EnVar methods. Using low-order models, we quantify the positive impact of the QS approach on the IEnKS, especially for long data assimilation windows. We also examine the computational cost of QS implementations and suggest cheaper algorithms.

Citation: Fillion, A., Bocquet, M., and Gratton, S.: Quasi static ensemble variational data assimilation, Nonlin. Processes Geophys. Discuss.,, in review, 2017.
Anthony Fillion et al.
Anthony Fillion et al.


Total article views: 372 (including HTML, PDF, and XML)

HTML PDF XML Total BibTeX EndNote
263 87 22 372 13 24

Views and downloads (calculated since 16 Nov 2017)

Cumulative views and downloads (calculated since 16 Nov 2017)

Viewed (geographical distribution)

Total article views: 371 (including HTML, PDF, and XML)

Thereof 371 with geography defined and 0 with unknown origin.

Country # Views %
  • 1



Latest update: 23 Mar 2018
Publications Copernicus
Short summary
This study generalizes a paper by Pires et al. (1996) to state-of-the-art data assimilation techniques, such as the iterative ensemble Kalman smoother (IEnKS). We show that the longer the time window over which observations are assimilated, the better the accuracy of the IEnKS. Beyond a critical time length that we estimate, we show that this accuracy finally degrades. We show that the use of the quasi-static minimizations, but generalized to the IEnKS yields a significantly improved accuracy.
This study generalizes a paper by Pires et al. (1996) to state-of-the-art data assimilation...