Variable-storage Quasi-Newton Operators as Inverse Forecast/Analysis Error Covariance Matrices in Variational Data Assimilation

Two approximations of the Hessian matrix as limited-memory operators are built from the limited-memory BFGS inverse Hessian approximation provided by the minimization code, in view of the specification of the inverse analysis/forecast error covariance matrix in variational data assimilation. Some numerical experiments and theoretical considerations lead to reject the limited-memory DFP Hessian approximation and to retain the BFGS one for the applications foreseen. Conditioning issues are explored and a preconditioning strategy via a change of control variable is proposed, based on a suitable Cholesky factorization of the limited-memory inverse Hessian matrix. This factorization is implemented as the composition of linear operators. The memory requirements and the number of floating-point operations required by the method are given and confirmed by numerical experiments. The method is found to have a strong potential for variational data assimilation systems using high resolution ocean or atmosphere general circulation models.

[1]  Monique Tanguay,et al.  Four‐dimensional data assimilation with a wide range of scales , 1995 .

[2]  Jorge Nocedal,et al.  Representations of quasi-Newton matrices and their use in limited memory methods , 1994, Math. Program..

[3]  Jacques Verron,et al.  Sensitivity Analysis in Variational Data Assimilation , 1997 .

[4]  Philippe Courtier,et al.  Data Assimilation in the Ocean and in the Atmosphere : What Should be Next? (gtSpecial IssueltData Assimilation in Meteology and Oceanography: Theory and Practice) , 1997 .

[5]  J. F. Bonnans Optimisation numérique : aspects théoriques et pratiques , 1997 .

[6]  Jorge Nocedal,et al.  Theory of algorithms for unconstrained optimization , 1992, Acta Numerica.

[7]  S. Cohn,et al.  An Introduction to Estimation Theory , 1997 .

[8]  Ionel M. Navon,et al.  Conjugate-Gradient Methods for Large-Scale Minimization in Meteorology , 1987 .

[9]  J. J. Moré,et al.  Quasi-Newton Methods, Motivation and Theory , 1974 .

[10]  Claude Lemaréchal,et al.  Some numerical experiments with variable-storage quasi-Newton algorithms , 1989, Math. Program..

[11]  P. Gauthier,et al.  A validation of the incremental formulation of 4D variational data assimilation in a nonlinear barotropic flow , 1998 .

[12]  F. Veerse,et al.  Multiple‐truncation incremental approach for four‐dimensional variational data assimilation , 1998 .

[13]  C.,et al.  Analysis methods for numerical weather prediction , 2022 .

[14]  M. J. D. Powell,et al.  Updating conjugate directions by the BFGS formula , 1987, Math. Program..

[15]  M. Zhu,et al.  The Quasi-Cauchy Relation and Diagonal Updating , 1999, SIAM J. Optim..

[16]  Paola Malanotte-Rizzoli,et al.  Assimilation studies of open‐ocean flows: 2. Error measures with strongly nonlinear dynamics , 1996 .

[17]  Philippe Courtier,et al.  Dynamical structure functions in a four‐dimensional variational assimilation: A case study , 1996 .

[18]  C. Wunsch The Ocean Circulation Inverse Problem , 1996 .

[19]  Pierre De Mey,et al.  Data Assimilation at the Oceanic Mesoscale: A Review , 1997 .

[20]  Chao Yang,et al.  ARPACK users' guide - solution of large-scale eigenvalue problems with implicitly restarted Arnoldi methods , 1998, Software, environments, tools.

[21]  Philippe Courtier,et al.  Four‐Dimensional Assimilation In the Presence of Baroclinic Instability , 1992 .

[22]  Jorge Nocedal,et al.  On the limited memory BFGS method for large scale optimization , 1989, Math. Program..

[23]  William Carlisle Thacker,et al.  The role of the Hessian matrix in fitting models to measurements , 1989 .

[24]  J. Nocedal Updating Quasi-Newton Matrices With Limited Storage , 1980 .

[25]  P. Courtier,et al.  A strategy for operational implementation of 4D‐Var, using an incremental approach , 1994 .