Singular vectors and estimates of the analysis‐error covariance metric

An important ingredient of ensemble forecasting is the computation of initial perturbatinos. Various techniques exist to generate initial perturbations. All these aim to produce an ensemble that, at initial time, reflects the uncertainty in the initial condition. In this paper a method for computing singular vectors consistent with current estimates of the analysis-error statistics is proposed and studied. The singular-vector computation is constrained at initial time by the Hessian of the three-dimensional variational assimilation (3D-Var) cost function in a way which is consistent with the operational analysis procedure. The Hessian is affected by the approximations made in the implementation of 3D-Var; however, it provides a more objective representation of the analysis-error covariances than other metrics previously used to constrain singular vectors. Experiments are performed with a T21L5 Primitive-Equation model. To compute the singular vectors we solve a generalized eigenvalue problem using a recently developed algorithm. It is shown that use of the Hessian of the cost function can significantly influence such properties of singular vectors as horizontal location, vertical structure and growth rate. The impact of using statistics of observational errors is clearly visible in that the amplitude of the singular vectors reduces in data-rich areas. Finally, the use of an approximation to the Hessian is discussed.

[1]  E. Lorenz.  A study of the predictability of a 28-variable atmospheric model , 1965 .

[2]  Philippe Courtier,et al.  Sensitivity of forecast errors to initial conditions , 1996 .

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

[4]  Martin Ehrendorfer,et al.  Optimal Prediction of Forecast Error Covariances through Singular Vectors , 1997 .

[5]  F. Bouttier.  A Dynamical Estimation of Forecast Error Covariances in an Assimilation System , 1994 .

[6]  P. Courtier,et al.  The ECMWF implementation of three‐dimensional variational assimilation (3D‐Var). II: Structure functions , 1998 .

[7]  P. Courtier,et al.  The ECMWF implementation of three‐dimensional variational assimilation (3D‐Var). I: Formulation , 1998 .

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

[9]  A. K. Cline.  Several Observations on the Use of Conjugate Gradient Methods , 1978 .

[10]  P. L. Houtekamer,et al.  Prediction Experiments with Two-Member Ensembles , 1994 .

[11]  Roberto Buizza,et al.  Sensitivity of optimal unstable structures , 1994 .

[12]  Andrew C. Lorenc,et al.  Analysis methods for numerical weather prediction , 1986 .

[13]  H. V. D. Vorst,et al.  Jacobi-davidson type methods for generalized eigenproblems and polynomial eigenproblems , 1995 .

[14]  C. Leith.  Theoretical Skill of Monte Carlo Forecasts , 1974 .

[15]  P. L. Houtekamer,et al.  A System Simulation Approach to Ensemble Prediction , 1996 .

[16]  Gilbert Strang,et al.  Introduction to applied mathematics , 1988 .

[17]  F. Molteni,et al.  The ECMWF Ensemble Prediction System: Methodology and validation , 1996 .

[18]  E. Davidson.  The iterative calculation of a few of the lowest eigenvalues and corresponding eigenvectors of large real-symmetric matrices , 1975 .

[19]  E. Kalnay,et al.  Ensemble Forecasting at NCEP and the Breeding Method , 1997 .

[20]  A. Simmons,et al.  The ECMWF medium-range prediction models development of the numerical formulations and the impact of increased resolution , 1989 .

[21]  T. Palmer,et al.  Singular Vectors, Metrics, and Adaptive Observations. , 1998 .

[22]  Roberto Buizza,et al.  Optimal perturbation time evolution and sensitivity of ensemble prediction to perturbation amplitude , 1995 .

[23]  C. Lanczos.  An iteration method for the solution of the eigenvalue problem of linear differential and integral operators , 1950 .

[24]  Roberto Buizza,et al.  The Singular-Vector Structure of the Atmospheric Global Circulation , 1995 .

[25]  John Derber,et al.  The National Meteorological Center's spectral-statistical interpolation analysis system , 1992 .

[26]  R. A. Walt,et al.  The numerical solution of algebraic equations , 1979 .

[27]  H. M. van den Dool,et al.  On the Weights for an Ensemble-Averaged 6–10-Day Forecast , 1994 .

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

[29]  Gerard L. G. Sleijpen,et al.  A Jacobi-Davidson Iteration Method for Linear Eigenvalue Problems , 1996, SIAM J. Matrix Anal. Appl..

[30]  E. Epstein,et al.  Stochastic dynamic prediction1 , 1969 .