On the Orthogonality of Bred Vectors

J. D. Annan Frontier Research System for Global Change, Yokohama, Kanagawa, Japan, and Proudman Oceanographic Laboratory, Liverpool, United Kingdom

Search for other papers by J. D. Annan in
Current site
Google Scholar
PubMed
Close
Restricted access

Abstract

The “breeding method” is a well-established and computationally inexpensive method for generating perturbations for ensemble integrations. One feature of bred vectors is that, although arbitrary initial perturbations generally converge toward the fastest growing directions, distinct bred vectors for global numerical weather prediction (NWP) models are generally quasi-orthogonal. However, an examination of the local structure of the vectors indicates that this conclusion may be somewhat misleading since the perturbations frequently have very similar shapes over local areas (consistent with the low local dimension of atmospheric physics). Therefore, there may be substantial redundancy when multiple independent breeding cycles are performed in parallel, and the vectors can be inefficient in spanning the range of locally growing perturbations. Experiments with the 3-variable Lorenz63 model indicate that orthogonalizing the bred vectors can result in significantly improved performance. A simple local orthogonalization technique has been applied to the spatially resolved 40-variable Lorenz96 model with encouraging results. It is suggested that a similar treatment to bred vectors might be helpful for global NWP or other applications.

Corresponding author address: J. D. Annan, Frontier Research System for Global Change, 3173-25 Showa-machi Kanazawa-ku, Yokohama, Kanagawa 236-0001, Japan. Email: jdannan@jamstec.go.jp

Abstract

The “breeding method” is a well-established and computationally inexpensive method for generating perturbations for ensemble integrations. One feature of bred vectors is that, although arbitrary initial perturbations generally converge toward the fastest growing directions, distinct bred vectors for global numerical weather prediction (NWP) models are generally quasi-orthogonal. However, an examination of the local structure of the vectors indicates that this conclusion may be somewhat misleading since the perturbations frequently have very similar shapes over local areas (consistent with the low local dimension of atmospheric physics). Therefore, there may be substantial redundancy when multiple independent breeding cycles are performed in parallel, and the vectors can be inefficient in spanning the range of locally growing perturbations. Experiments with the 3-variable Lorenz63 model indicate that orthogonalizing the bred vectors can result in significantly improved performance. A simple local orthogonalization technique has been applied to the spatially resolved 40-variable Lorenz96 model with encouraging results. It is suggested that a similar treatment to bred vectors might be helpful for global NWP or other applications.

Corresponding author address: J. D. Annan, Frontier Research System for Global Change, 3173-25 Showa-machi Kanazawa-ku, Yokohama, Kanagawa 236-0001, Japan. Email: jdannan@jamstec.go.jp

Save
  • Anderson, J. L., 2001: An ensemble adjustment Kalman filter for data assimilation. Mon. Wea. Rev., 129 , 28842902.

  • Bishop, C. H., B. J. Etherton, and S. J. Majumdar, 2001: Adaptive sampling with the ensemble transform Kalman filter. Part I: Theoretical aspects. Mon. Wea. Rev., 129 , 420436.

    • Search Google Scholar
    • Export Citation
  • Cohn, S. E., and R. Todling, 1996: Approximate data assimilation for stable and unstable dynamics. J. Meteor. Soc. Japan, 74 , 6375.

  • Corazza, M., E. Kalnay, D. J. Patil, E. Ott, J. Yorke, I. Szunyogh, and M. Cai, 2002: Use of the breeding technique in the estimation of the background error covariance matrix for a quasigeostrophic model. Preprints, Symp. on Observations, Data Assimilation, and Probabilistic Prediction, Orlando, FL, Amer. Meteor. Soc., 154–157.

    • Search Google Scholar
    • Export Citation
  • Evensen, G., 1994: Sequential data assimilation with a nonlinear quasi-geostrophic model using Monte Carlo methods to forecast error statistics. J. Geophys. Res., 99 (C5) 1014310162.

    • Search Google Scholar
    • Export Citation
  • Heemink, A. W., K. Bolding, and M. Verlaan, 1997: Storm surge forecasting using Kalman filtering. J. Meteor. Soc. Japan, 75 , 305318.

  • Houtekamer, P. L., and H. L. Mitchell, 1998: Data assimilation using an ensemble Kalman filter technique. Mon. Wea. Rev., 126 , 796811.

    • Search Google Scholar
    • Export Citation
  • Kalnay, E., and Z. Toth, 1994: Removing growing errors in the analysis. Preprints, 10th Conf. on Numerical Weather Prediction, Portland, OR, Amer. Meteor. Soc., 212–215.

    • Search Google Scholar
    • Export Citation
  • Kalnay, E., M. Corazza, and M. Cai, 2002: Are bred vectors the same as Lyapunov vectors? Preprints, Symp. on Observations, Data Assimilation, and Probabilistic Prediction, Orlando, FL, Amer. Meteor. Soc., 173–177.

    • Search Google Scholar
    • Export Citation
  • Kincaid, D. R., and E. W. Cheney, 1990: Numerical Analysis: The Mathematics of Scientific Computing. Brooks/Cole, 690 pp.

  • Lorenz, E. N., 1963: Deterministic non-periodic flow. J. Atmos. Sci., 20 , 130141.

  • Lorenz, E. N., 1996: Predictability: A problem partly solved. Proc. ECMWF Seminar on Predictability, Reading, United Kingdom, ECMWF, 1–18.

    • Search Google Scholar
    • Export Citation
  • Lorenz, E. N., and K. A. Emmanuel, 1998: Optimal sites for supplementary weather observations: Simulation with a small model. J. Atmos. Sci., 55 , 399414.

    • Search Google Scholar
    • Export Citation
  • Patil, D. J., B. R. Hunt, E. Kalnay, J. A. Yorke, and E. Ott, 2001: Local low dimensionality of atmospheric dynamics. Phys. Rev. Lett., 86 , 58785881.

    • Search Google Scholar
    • Export Citation
  • Peña, M., and E. Kalnay, 2003: Separating fast and slow modes in coupled chaotic systems. Geophys. Res. Abstr., 5 , 07540.

  • Pham, D. T., J. Verron, and M. C. Roubaud, 1998: A singular evolutive extended Kalman filter for data assimilation in oceanography. J. Mar. Syst., 16 , 323340.

    • Search Google Scholar
    • Export Citation
  • Toth, Z., and E. Kalnay, 1993: Ensemble forecasting at NMC: The generation of perturbations. Bull. Amer. Meteor. Soc., 74 , 23172330.

  • Toth, Z., and E. Kalnay, 1997: Ensemble forecasting at NCEP and the breeding method. Mon. Wea. Rev., 125 , 32973319.

  • van Leeuwen, P. J., 1999: Comments on “Data assimilation using an ensemble Kalman filter.”. Mon. Wea. Rev., 127 , 13741377.

  • Wang, X., and C. H. Bishop, 2003: A comparison of breeding and ensemble transform Kalman filter ensemble forecast schemes. J. Atmos. Sci., 60 , 11401158.

    • Search Google Scholar
    • Export Citation
  • Yamane, S., and S. Yoden, 2001: Finite-time evolution of small perturbations superposed on a chaotic solution: Experiments with an idealized barotropic model. J. Atmos. Sci., 58 , 10661078.

    • Search Google Scholar
    • Export Citation
All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 1271 926 29
PDF Downloads 185 66 12