Optimal solution error covariance in highly nonlinear problems of variational data assimilation

Shutyaev, V. and Gejadze, I and Copeland, G.J.M and Le Dimet, F.X. (2012) Optimal solution error covariance in highly nonlinear problems of variational data assimilation. Nonlinear Processes in Geophysics, 19 (2). pp. 177-184. ISSN 1023-5809 (https://doi.org/10.5194/npg-19-177-2012)

[thumbnail of npg_19_177_2012.pdf]
Preview
PDF. Filename: npg_19_177_2012.pdf
Final Published Version

Download (877kB)| Preview

Abstract

The problem of variational data assimilation for a nonlinear evolution model is formulated as an optimal control problem (see, e.g.[1]) to find the initial condition, boundary conditions or model parameters. The input data contain observation and background errors, hence there is an error in the optimal solution. For mildly nonlinear dynamics, the covariance matrix of the optimal solution error can be approximated by the inverse Hessian of the cost functional of an auxiliary data assimilation problem ([2], [3]). The relationship between the optimal solution error covariance matrix and the Hessian of the auxiliary control problem is discussed for different degrees of validity of the tangent linear hypothesis. For problems with strongly nonlinear dynamics a new statistical method based on computation of a sample of inverse Hessians is suggested. This method relies on the efficient computation of the inverse Hessian by means of iterative methods (Lanczos and quasi-Newton BFGS) with preconditioning. The method allows us to get a sensible approximation of the posterior covariance matrix with a small sample size. Numerical examples are presented for the model governed by Burgers equation with a nonlinear viscous term. The first author acknowledges the funding through the project 09-01-00284 of the Russian Foundation for Basic Research, and the FCP program "Kadry".