Conditioning and preconditioning of the variational data assimilation problem.

*(English)*Zbl 1433.86007Summary: Numerical weather prediction (NWP) centres use numerical models of the atmospheric flow to forecast future weather states from an estimate of the current state. Variational data assimilation (VAR) is used commonly to determine an optimal state estimate that miminizes the errors between observations of the dynamical system and model predictions of the flow. The rate of convergence of the VAR scheme and the sensitivity of the solution to errors in the data are dependent on the condition number of the Hessian of the variational least-squares objective function. The traditional formulation of VAR is ill-conditioned and hence leads to slow convergence and an inaccurate solution. In practice, operational NWP centres precondition the system via a control variable transform to reduce the condition number of the Hessian. In this paper we investigate the conditioning of VAR for a single, periodic, spatially-distributed state variable. We present theoretical bounds on the condition number of the original and preconditioned Hessians and hence demonstrate the improvement produced by the preconditioning. We also investigate theoretically the effect of observation position and error variance on the preconditioned system and show that the problem becomes more ill-conditioned with increasingly dense and accurate observations. Finally, we confirm the theoretical results in an operational setting by giving experimental results from the Met Office variational system.

##### Keywords:

variational data assimilation; error covariances; least-squares optimization; condition number; Hessian
PDF
BibTeX
XML
Cite

\textit{S. A. Haben} et al., Comput. Fluids 46, No. 1, 252--256 (2011; Zbl 1433.86007)

Full Text:
DOI

##### References:

[1] | Andersson, E.; Fisher, M.; Munro, R.; McNally, A., Diagnosis of background errors for radiances and other observable quantities in a variational data assimilation scheme, and the explanation of a case of poor convergence, Quart J roy met soc, 126, 1455-1472, (2000) |

[2] | Gauthier, P.; Charette, C.; Fillion, L.; Koclas, P.; Laroche, S., Implementation of a 3D variational data assimilation system at the Canadian meteorological centre. part 1: the global analysis, Atmos Ocean, 37, 103-156, (1999) |

[3] | Golub, G.H.; Van Loan, C.F., Matrix computations, (1996), Johns Hopkins University Press · Zbl 0865.65009 |

[4] | Gray, R.M., Toeplitz and circulant matrices: a review, Found trendsÂ® commun inform theory, 2, 3, 155-239, (2006) |

[5] | Haben SA, Lawless AS, Nichols NK. Conditioning of the 3DVar data assimilation problem. University of Reading, Dept. of Mathematics, Math Report Series 3/2009; 2009. <http://www.reading.ac.uk/maths-and-stats/research/maths-report-series.aspx>. |

[6] | Lawless, A.S.; Gratton, S.; Nichols, N.K., An investigation of incremental 4D-var using non-tangent linear models, Quart J roy met soc, 131, 459-476, (2005) |

[7] | Lorenc, A., Optimal nonlinear objective analysis, Quart J roy met soc, 114, 205-240, (1988) |

[8] | Lorenc, A.C., Development of an operational variational assimilation scheme, J met soc jpn, 75, 339-346, (1997) |

[9] | Lorenc, A.C.; Ballard, S.P.; Bell, R.S.; Ingleby, N.B.; Andrews, P.L.F., The met office global three-dimensional variational data assimilation scheme, Quart J roy met soc, 126, 2991-3012, (2000) |

[10] | Rabier, F.; Jarvinen, H.; Klinker, E.; Mahfouf, J.F.; Simmons, A., The ECMWF operational implementation of four-dimensional variational assimilation. I: experimental results with simplified physics, Quart J roy met soc, 126, 1143-1170, (2000) |

[11] | Rawlins, F.; Ballard, S.P.; Bovis, K.J.; Clayton, A.M.; Li, D.; Inverarity, G.W., The met office global four-dimensional variational data assimilation scheme, Quart J roy met soc, 133, 347-362, (2007) |

[12] | Tremolet, Y., Incremental 4D-var convergence study, Tellus A, 59, 706-718, (2007) |

This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.