Adaptive Observations And Multilevel Optimization In Data Assimilation

We propose to use a decomposition of large-scale incremental four dimensional (4D-Var) data assimilation problems in order to make their numerical solution more efficient. This decomposition is based on exploiting an adaptive hierarchy of the observations. Starting with a low-cardinality set and the solution of its corresponding optimization problem, observations are adaptively added based on a posteriori error estimates. The particular structure of the sequence of associated linear systems allows the use of a variant of the conjugate gradient algorithm which effectively exploits the fact that the number of observations is smaller than the size of the vector state in the 4D-Var model. The method proposed is justified by deriving the relevant error estimates at different levels of the hierarchy and a practical computational technique is then derived. The new algorithm is tested on a 1D-wave equation and on the Lorenz-96 system, the latter one being of special interest because of its similarity with Numerical Weather Prediction (NWP) systems.

[1]  Philippe Courtier,et al.  Dual formulation of four‐dimensional variational assimilation , 1997 .

[2]  E. Lorenz.  Predictability of Weather and Climate: Predictability – a problem partly solved , 2006 .

[3]  Anders Logg,et al.  Automated Solution of Differential Equations by the Finite Element Method: The FEniCS Book , 2012 .

[4]  A. Brandt.  Multi-level adaptive technique (MLAT) for fast numerical solution to boundary value problems , 1973 .

[5]  Klaus Böhmer,et al.  Defect Correction Methods , 1984, Computing Supplementum.

[6]  A. Karimi,et al.  Extensive chaos in the Lorenz-96 model. , 2009, Chaos.

[7]  Erik Andersson,et al.  Influence‐matrix diagnostic of a data assimilation system , 2004 .

[8]  K. Emanuel,et al.  Optimal Sites for Supplementary Weather Observations: Simulation with a Small Model , 1998 .

[9]  W. Dörfler.  A convergent adaptive algorithm for Poisson's equation , 1996 .

[10]  Leonard A. Smith,et al.  The role of operational constraints in selecting supplementary observations , 2000 .

[11]  Chris Snyder,et al.  Statistical Design for Adaptive Weather Observations , 1999 .

[12]  Samuel Buis,et al.  Intercomparison of the primal and dual formulations of variational data assimilation , 2008 .

[13]  S. McCormick,et al.  Fast Adaptive Composite Grid (FAC) Methods: Theory for the Variational Case , 1984 .

[14]  Takemasa Miyoshi,et al.  Analysis sensitivity calculation in an ensemble Kalman filter , 2009 .

[15]  Serge Gratton,et al.  Approximate Gauss-Newton Methods for Nonlinear Least Squares Problems , 2007, SIAM J. Optim..

[16]  Dacian N. Daescu,et al.  Adaptive observations in the context of 4D-Var data assimilation , 2004 .

[17]  Serge Gratton,et al.  Preconditioning and globalizing conjugate gradients in dual space for quadratically penalized nonlinear-least squares problems , 2013, Comput. Optim. Appl..

[18]  Serge Gratton,et al.  An observation‐space formulation of variational assimilation using a restricted preconditioned conjugate gradient algorithm , 2009 .

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

[20]  Ricardo H. Nochetto,et al.  Data Oscillation and Convergence of Adaptive FEM , 2000, SIAM J. Numer. Anal..

[21]  M. Hestenes,et al.  Methods of conjugate gradients for solving linear systems , 1952 .

[22]  M. Leutbecher.  Adaptive observations, the Hessian metric and singular vectors , 2003 .

[23]  Shinfield Park.  Adaptive observations , the Hessian metric and singular vectors , 2003 .

[24]  Gene H. Golub,et al.  Matrix computations , 1983 .

[25]  M. Leutbecher.  A Reduced Rank Estimate of Forecast Error Variance Changes due to Intermittent Modifications of the Observing Network. , 2003 .