An inexpensive technique for using past forecast errors to improve future forecast skill Part II, A quasi-inverse linear method
Advanced Search
Select up to three search categories and corresponding keywords using the fields to the right. Refer to the Help section for more detailed instructions.

Search our Collections & Repository

All these words:

For very narrow results

This exact word or phrase:

When looking for a specific result

Any of these words:

Best used for discovery & interchangable words

None of these words:

Recommended to be used in conjunction with other fields



Publication Date Range:


Document Data


Document Type:






Clear All

Query Builder

Query box

Clear All

For additional assistance using the Custom Query please check out our Help Page


An inexpensive technique for using past forecast errors to improve future forecast skill Part II, A quasi-inverse linear method

Filetype[PDF-2.67 MB]


  • Alternative Title:
    A quasi-inverse linear method
  • Description:
    A quasi-inverse linear method has been developed to study the sensitivity forecast errors to initial conditions for the National Centers for Environmental Prediction (NCEP)'s global spectral model. The inverse is approximated by running the tangent linear model (TLM) of the nonlinear forecast model with a negative time step, but reversing the sign of friction and diffusion terms, in order to avoid the computational instability that would be associated with these terms if they were run backwards. As usually done using the adjoint model integrations, the quasi-inverse TLM is started at the time of the verified forecast error and integrated backwards to the corresponding initial time. First, a numerical experiment shows that this quasi-inverse linear estimation is able to trace back the differences between two perturbed forecasts from the NCEP ensemble forecasting system, and recover with good accuracy the known difference between the two forecasts at the initial time. This result shows that both the linear estimation and the quasi-inverse linear estimation are quite close to the nonlinear evolution of the perturbation in the nonlinear forecast model, suggesting that we should be able to apply the method to study of the sensitivity of forecast errors to initial conditions. We then calculate the perturbation field at the initial time (linear sensitivity perturbation) by tracing back a one-day forecast error using the TLM quasi-inverse estimation. As could be expected from the previous experiment, when the estimated error is subtracted from the original analysis, the new initial conditions lead to an almost perfect one-day forecast. The forecasts beyond day one are also considerably improved, indicating that the initial conditions have indeed been improved. In the remainder of the paper, this quasi-inverse linear sensitivity method is compared with the adjoint sensitivity method (Pu et al. 1995; Rabier et al. 1996) for medium range weather forecasting. We find that both methods are able to trace back the forecast error to sensitivity perturbations which improve the initial conditions. However, the forecast improvement obtained by the quasi-inverse linear method is considerably better than that obtained with a single adjoint iteration, and similar to the one obtained using 5 iterations of the adjoint method, even though each adjoint iteration requires at least twice the computer resources of the quasi-inverse TLM estimation. Whereas the adjoint forecast sensitivities are closely related to singular vectors, the quasi-inverse linear sensitivities are associated with the bred (Lyapunov) vectors used for ensemble forecasting at NCEP (Toth and Kalnay 1993). The features of the two type sensitivity perturbations are also compared in this study. Finally, the possibility of the use of the sensitivity perturbation to improve future forecast skill is discussed, and preliminary experiments encourage us to further test this rather inexpensive method for possible operational use. The model used in this study is the NCEP operational global spectral model at a resolution of T62/L28. The corresponding TLM, and its adjoint, are based on an adiabatic version of the model but include both horizontal and vertical diffusion.
  • Content Notes:
    Zhao-Xia Pu, Eugenia Kalnay.

    "March 1996."

    "This is an unreviewed manuscript, primarily intended for informal exchange of information among NCEP staff members."

    System requirements: Adobe Acrobat Reader.

    Includes bibliographical references (pages 21-23).

  • Document Type:
  • Place as Subject:
  • Rights Information:
    Public Domain
  • Compliance:
  • Main Document Checksum:
  • File Type:

Supporting Files

  • No Additional Files

More +

You May Also Like

Checkout today's featured content at

Version 3.26