HEC Montréal, Canada, May 6  8, 2013
2013 Optimization Days
HEC Montréal, Canada, 6 — 8 May 2013
TBP Séance plénière IV/ Plenary Session IV
May 7, 2013 02:00 PM – 03:00 PM
Location: Amphithéâtre IBM
Chaired by Dominique Orban
1 Presentation

02:00 PM  03:00 PM
Convergence Properties of the Primal and Dual Forms of Variational Data Assimilation
Our understanding of the fundamental processes governing the atmosphere and their interactions is achieved through observations and modeling based on first principles. Numerical modeling is the laboratory in which our knowledge is encompassed within a mathematical and numerical framework, the results of which are then confronted to observations through data assimilation, in which observations and model predictions are combined based on their relative accuracy to produce analyses representing our best estimate of the current state of the system. This paper discusses the variational data assimilation problem which can be solved in either its primal (3D/4DVar) or dual (3D/4DPSAS) form. The methods are equivalent at convergence but the dual method exhibits a spurious behaviour at the beginning of the minimization which leads to less probable states than the background state. This is a serious concern when using the dual method in operational implementations when only a finite number of iterations can be afforded. Two classes of minimization algorithms are examined in this article: the conjugate gradient (CG) and the minimum residual (MINRES) methods. While the CG algorithms ensure a monotonic reduction of the cost function, those based on the MINRES enforce instead a monotonic decrease of the norm of the gradient. In this article, it is shown that when applied to the minimization of the dual problem, the MINRES algorithms also lead to iterates for which their ‘image’ in physical space leads to a monotonic decrease of the primal cost function. A relationship is established showing that the primal objective function is related to the value of the dual cost function and the norm of its gradient. This holds for the incremental forms of both the three and fourdimensional cases. A new convergence criterion is introduced based on the error norm in model space to make sure that, for the dual problem, the same accuracy is obtained in the analysis when only a finite number of iterations is completed.