On optimal solution error covariances in variational data assimilation problems

I.Y. Gejadze, F.X. Le Dimet, V. Shutyaev, Scottish Founding Council via GRPE (Funder)

Research output: Contribution to journalArticlepeer-review

38 Citations (Scopus)
56 Downloads (Pure)


The problem of variational data assimilation for a nonlinear evolution model is formulated as an optimal control problem to find unknown parameters such as distributed model coefficients or boundary conditions. The equation for the optimal solution error is derived through the errors of the input data (background and observation errors), and the optimal solution error covariance operator through the input data error covariance operators, respectively. The quasi-Newton BFGS algorithm is adapted to construct the covariance matrix of the optimal solution error using the inverse Hessian of an auxiliary data assimilation problem based on the tangent linear model constraints. Preconditioning is applied to reduce the number of iterations required by the BFGS algorithm to build a quasi-Newton approximation of the inverse Hessian. Numerical examples are presented for the one-dimensional convection-diffusion model.
Original languageEnglish
Pages (from-to)2159-2178
Number of pages20
JournalJournal of Computational Physics
Issue number6
Publication statusPublished - 20 Mar 2010


  • variational data assimilation
  • parameter estimation
  • optimal solution error covariances
  • hessian
  • preconditioning
  • mathematics


Dive into the research topics of 'On optimal solution error covariances in variational data assimilation problems'. Together they form a unique fingerprint.

Cite this