Assimilation of Standard and Targeted Observations within the Unstable Subspace of the Observation–Analysis–Forecast Cycle System

Anna Trevisan CNR-ISAC, Bologna, Italy

Search for other papers by Anna Trevisan in
Current site
Google Scholar
PubMed
Close
and
Francesco Uboldi LEGOS/SHOM, Toulouse, France

Search for other papers by Francesco Uboldi in
Current site
Google Scholar
PubMed
Close
Restricted access

Abstract

In this paper it is shown that the flow-dependent instabilities that develop within an observation–analysis–forecast (OAF) cycle and that are responsible for the background error can be exploited in a very simple way to assimilate observations. The basic idea is that, in order to minimize the analysis and forecast errors, the analysis increment must be confined to the unstable subspace of the OAF cycle solution. The analysis solution here formally coincides with that of the classical three-dimensional variational solution with the background error covariance matrix estimated in the unstable subspace.

The unstable directions of the OAF system solution are obtained by breeding initially random perturbations of the analysis but letting the perturbed trajectories undergo the same process as the control solution, including assimilation of all the available observations. The unstable vectors are then used both to target observations and for the assimilation design.

The approach is demonstrated in an idealized environment using a simple model, simulated standard observations over land with a single adaptive observation over the ocean. In the application a simplified form is adopted of the analysis solution and a single unstable vector at each analysis time whose amplitude is determined by means of the adaptive observation. The remarkable reduction of the analysis and forecast error obtained by this simple method suggests that only a few accurately placed observations are sufficient to control the local instabilities that take place along the cycle.

The stability of the system, with or without forcing by observations, is studied and the growth rate of the leading instability of the different control solutions is estimated. Whereas the model has more than one positive Lyapunov exponent, the solution of the OAF scheme that includes the adaptive observation is stable. It is suggested that a negative exponent can be considered a necessary condition for the convergence of a particular OAF solution to the truth, and that the estimate of the degree of stability of the control trajectory can be used as a simple criterion to evaluate the efficiency of data assimilation and observation strategies.

The present findings are in line with previous quantative observability results with more realistic models and with recent studies that indicate a local low dimensionality of the unstable subspace.

Corresponding author address: Dr. Anna Trevisan, CNR-ISAC, Via Gobetti 101, 40129 Bologna, Italy. Email: a.trevisan@isac.cnr.it

Abstract

In this paper it is shown that the flow-dependent instabilities that develop within an observation–analysis–forecast (OAF) cycle and that are responsible for the background error can be exploited in a very simple way to assimilate observations. The basic idea is that, in order to minimize the analysis and forecast errors, the analysis increment must be confined to the unstable subspace of the OAF cycle solution. The analysis solution here formally coincides with that of the classical three-dimensional variational solution with the background error covariance matrix estimated in the unstable subspace.

The unstable directions of the OAF system solution are obtained by breeding initially random perturbations of the analysis but letting the perturbed trajectories undergo the same process as the control solution, including assimilation of all the available observations. The unstable vectors are then used both to target observations and for the assimilation design.

The approach is demonstrated in an idealized environment using a simple model, simulated standard observations over land with a single adaptive observation over the ocean. In the application a simplified form is adopted of the analysis solution and a single unstable vector at each analysis time whose amplitude is determined by means of the adaptive observation. The remarkable reduction of the analysis and forecast error obtained by this simple method suggests that only a few accurately placed observations are sufficient to control the local instabilities that take place along the cycle.

The stability of the system, with or without forcing by observations, is studied and the growth rate of the leading instability of the different control solutions is estimated. Whereas the model has more than one positive Lyapunov exponent, the solution of the OAF scheme that includes the adaptive observation is stable. It is suggested that a negative exponent can be considered a necessary condition for the convergence of a particular OAF solution to the truth, and that the estimate of the degree of stability of the control trajectory can be used as a simple criterion to evaluate the efficiency of data assimilation and observation strategies.

The present findings are in line with previous quantative observability results with more realistic models and with recent studies that indicate a local low dimensionality of the unstable subspace.

Corresponding author address: Dr. Anna Trevisan, CNR-ISAC, Via Gobetti 101, 40129 Bologna, Italy. Email: a.trevisan@isac.cnr.it

Save
  • Bennett, A. F., 2002: Inverse Modeling of the Ocean and Atmosphere. Cambridge University Press, 256 pp.

  • Bergot, T., G. Hello, A. Joly, and S. Marlardel, 1999: Adaptive observations: A feasibility study. Mon. Wea. Rev., 127 , 743765.

  • Berliner, L. M., Z-Q. Lu, and C. Snyder, 1999: Statistical design for adaptive weather observations: A feasibility study. J. Atmos. Sci., 56 , 25362552.

    • Search Google Scholar
    • Export Citation
  • Bishop, C. H., and Z. Toth, 1999: Ensemble transformation and adaptive observations. J. Atmos. Sci., 56 , 17481765.

  • Daley, R., 1991: Atmospheric Data Analysis. Cambridge University Press, 457 pp.

  • Dee, D., S. E. Cohn, A. Dalcher, and M. Ghil, 1985: An efficient algorithm for estimating covariances in distributed systems. IEEE Trans. Autom. Control, 30AC , 10571065.

    • Search Google Scholar
    • Export Citation
  • Evensen, G., and P. J. van Leewen, 1996: Assimilation of Geosat altimeter data for the Agulhas Current using the ensemble Kalman filter with a quasigeostrophic model. Mon. Wea. Rev., 124 , 8596.

    • Search Google Scholar
    • Export Citation
  • Gelaro, R., R. Langland, G. D. Rohaly, and T. E. Rosmond, 1999: An assessment of the singular vector approach to targeted observing using the FASTEX dataset. Quart. J. Roy. Meteor. Soc., 125 , 32993328.

    • Search Google Scholar
    • Export Citation
  • Ghil, M., 1989: Meteorological data assimilation for oceanographers. Part I: Description and theoretical framework. Dyn. Atmos. Oceans, 13 , 171218.

    • Search Google Scholar
    • Export Citation
  • Ghil, M., 1997: Advances in sequential estimation for atmospheric and oceanic flows. J. Meteor. Soc. Japan, 75 , 289304.

  • Ghil, M., and P. Malanotte-Rizzoli, 1991: Data assimilation in meteorology and oceanography. Adv. Geophys., 33 , 141266.

  • Ghil, M., and R. Todling, 1996: Tracking atmospheric instabilities with the Kalman filter. Part II: Two-layer results. Mon. Wea. Rev., 124 , 23402352.

    • Search Google Scholar
    • Export Citation
  • Ghil, M., S. Cohn, J. Tavantzis, K. Bube, and E. Isaacson, 1981: Applications of estimation theory to numerical weather prediction. Dynamic Meteorology: Data Assimilation Methods, L. Bengtsson, M. Ghil, and E. Kallen, Eds., Springer Verlag, 139–224.

    • Search Google Scholar
    • Export Citation
  • Hamill, T. M., and C. Snyder, 2000: A hybrid ensemble Kalman filter-3D-variational analysis scheme. Mon. Wea. Rev., 128 , 29052919.

  • Hansen, J. A., and A. L. Smith, 2000: The role of operational constraints in selecting supplementary observations. J. Atmos. Sci., 57 , 28592871.

    • Search Google Scholar
    • Export Citation
  • Houtekamer, P. L., and H. L. Mitchell, 2001: A sequential ensemble Kalman filter for atmospheric data assimilation. Mon. Wea. Rev., 129 , 796811.

    • Search Google Scholar
    • Export Citation
  • Ide, K., P. Courtier, M. Ghil, and A. C. Lorenc, 1997: Unified notation for data assimilation: Operational, sequential and variational. J. Meteor. Soc. Japan, 75 , 181189.

    • Search Google Scholar
    • Export Citation
  • Kalnay, E., 2002: Atmospheric Modeling, Data Assimilation, and Predictability. Cambridge University Press, 314 pp.

  • 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
  • Lacarra, J., and O. Talagrand, 1988: Short range evolution of small perturbations in a barotropic model. Tellus, 40A , 8195.

  • Langland, R. H., R. Gelaro, G. D. Rohaly, and M. A. Shapiro, 1999: Targeted observations in FASTEX: Adjoint-based targeting procedures and data impact experiments in IOP 17 and IOP 18. Quart. J. Roy. Meteor. Soc., 125 , 32413270.

    • Search Google Scholar
    • Export Citation
  • Lorenz, E. N., 1984: The local structure of a chaotic attractor in four dimensions. Physica D, 13 , 90104.

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

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

    • Search Google Scholar
    • Export Citation
  • Montani, A., A. J. Thorpe, R. Buizza, and P. Unden, 1999: Forecast skill of the ECMWF model using targeted observations during FASTEX. Quart. J. Roy. Meteor. Soc., 125 , 32193240.

    • Search Google Scholar
    • Export Citation
  • Morss, R. E., K. A. Emanuel, and C. Snyder, 2001: Idealized adaptive observation strategies for improving numerical weather prediction. J. Atmos. Sci., 58 , 210232.

    • Search Google Scholar
    • Export Citation
  • Navon, I. M., 1997: Practical and theoretical aspects of adjoint parameter estimation and identifiability in meteorology and oceanography. Dyn. Atmos. Oceans, 27 , 5579.

    • Search Google Scholar
    • Export Citation
  • Palmer, T. N., R. Gelaro, J. Barkmeijer, and R. Buizza, 1998: Singular vectors, metrics, and adaptive observations. J. Atmos. Sci., 55 , 633653.

    • 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
  • Pu, Z. X., and E. Kalnay, 1999: Targeting observations with the quasi-linear inverse and adjoint NCEP global models: Performance during FASTEX. Quart. J. Roy. Meteor. Soc., 125 , 33293337.

    • Search Google Scholar
    • Export Citation
  • Pu, Z. X., E. Kalnay, J. Sela, and I. Szuniog, 1997: Sensitivity of forecast errors to initial conditions with a quasi-inverse linear method. Mon. Wea. Rev., 125 , 24792503.

    • Search Google Scholar
    • Export Citation
  • Snyder, C., 1996: Summary of an informal workshop on adaptive observations and FASTEX. Bull. Amer. Meteor. Soc., 77 , 953961.

  • Szunyog, I., Z. Toth, K. A. Emanuel, C. Bishop, J. Woolen, T. Marchok, R. Morss, and C. Snyder, 1999: Ensemble based targeting experiments during FASTEX: The impact of dropsonde data from the Lear jet. Quart. J. Roy. Meteor. Soc., 125 , 31893218.

    • Search Google Scholar
    • Export Citation
  • Todling, R., and M. Ghil, 1994: Tracking atmospheric instabilities with the Kalman filter. Part I: Methodology and one-layer results. Mon. Wea. Rev., 122 , 183204.

    • Search Google Scholar
    • Export Citation
  • Trevisan, A., and F. Pancotti, 1998: Periodic orbits, Lyapunov vectors, and singular vectors in the Lorenz system. J. Atmos. Sci., 55 , 390398.

    • Search Google Scholar
    • Export Citation
  • Uboldi, F., and M. Kamachi, 2000: Time–space weak-constraint data assimilation for nonlinear models. Tellus, 52A , 412421.

All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 881 350 18
PDF Downloads 185 50 9