Variational assimilation of meteorological observations with the adjoint vorticity equation. I: Theory

Abstract

The following variational approach is taken to the problem of assimilation of meteorological observations: find the solution of the assimilating model which minimizes a given scalar function measuring the ‘distance’ between a model solution and the available observations. It is shown how the ‘adjoint equations’ of the model can be used to compute explicitly the ‘gradient’ of the distance function with respect to the model's initial conditions. the computation of one gradient requires one forward integration of the full model equations over the time interval on which the observations are available, followed by one backward integration of the adjoint equations. Successive gradients thus computed are introduced into a descent algorithm in order to determine the initial conditions which define the minimizing model solution. The theory is applied to the vorticity equation. Successful numerical experiments performed on a Haurwitz wave are described.