COLA Report 67

Dynamical Seasonal Predictions with the COLA Atmospheric model

J. Shukla, D.A. Paolino, D.M. Straus, D. DeWitt, M. Fennessy, J.L. Kinter, L. Marx and R. Mo

May 1999


Abstract

Seasonal hindcasts were made using the COLA model for 16 winter seasons (mid-December through March for 1981/82 - 1996/97). For each season, a nine member ensemble was generated using observed initial conditions in mid-December and observed global sea surface temperature (SST). It is found that in the presence of large tropical SST anomalies, the model is quite successful in simulating seasonal mean height anomalies over the Pacific North America (PacNA) region.

A local spatial pattern correlation field is computed for the ensemble, seasonal mean of 500 hPa height from the GCM and the seasonal means from the reanalyses of the National Centers for Environmental Prediction, hereafter observations. This field exceeds 0.6 over the eastern tropical Pacific and western North America; the maximum values are greatly enhanced during ENSO years. Similar results are obtained for 200 hPa u-winds. The observed enhancement of intraseasonal low pass (10-90 day) variability of 500 hPa height during cold events is simulated, as is the shift in the stormtracks (2-10 day variability of 850 hPa meridional heat flux) during warm events.

Empirical orthogonal function (EOF) analysis is applied in the PacNA region to the ensemble seasonal means of GCM 500 hPa height and the corresponding observed seasonal means; the global combined Arctic / North Atlantic Oscillation is removed from both GCM and observations. The leading EOF mode explains about 50% of the variance for both GCM and observations; the two patterns are nearly identical. Singular Value Decomposition (SVD) analysis between the tropical Pacific SST and 500 hPa height in the PacNA region shows that the nature of the coupling between SST and height is nearly the same in the GCM and observations. The great similarity between the height patterns in the first SVD mode and the corresponding leading mode EOF patterns indicates that the leading height variations are forced by SST.

The SST-forced variance of height was also estimated by regression analysis of (ensemble) seasonal means for the GCM and observations for the 16 years onto an index of tropical Pacific SSTs derived from SVD analysis of a long (30 year) observed record of heights and SSTs. The pattern of percentage variance explained in the GCM and observations are very similar to each other (and to the EOF described above). The higher absolute values in the GCM case reflect the effectiveness of the ensemble in filtering out variability unrelated to SST forcing.

SVD analysis was applied to100 GCM "samples" coupled with the observations; a sample is defined as a single 16-year record obtained by picking one ensemble member randomly for each of the 16 years. Probability distribution functions (pdf) of the pattern correlation for the leading SVD patterns, the percentage explained squared covariance, and the time series correlations all indicated sharp peaks at values of (0.87, 87%, and 0.82) respectively.

The pdf of the projection of individual 5-day means onto the leading EOF described above is quite dramatically shifted during strong warm and cold tropical SST events; the warm [cold] event pdf has almost all its weight in the negative [positive] EOF region. GCM and observations agree well.

The intra-ensemble spread was estimated by computing the PacNA anomaly correlation coefficient (ACC) for each of 36 possible intra-ensemble pairs. For nine years in which histograms of the ACC indicate predominantly positive values, the ACC of the ensemble mean with the observed seasonal mean is also relatively high.

Brier skill scores and reliability diagrams were computed for the "event" of the 500 hPa height being one standard deviation above (or below) the normal, with all such events pooled over the entire Northern Hemisphere or North America only. All skill scores are positive and statistically significant at the 99% level. The North America scores are higher than the whole hemisphere scores; the ENSO year scores are higher than those considering all years.

Complete copies of this report are available from:

Center for Ocean-Land-Atmosphere Studies
4041 Powder Mill Road, Suite 302
Calverton, MD 20705-3106
(301) 595-7000
(301) 595-9793 Fax
jperez@cola.iges.org

Full Report


last update: 9 June 1999
comments to: www@grads.iges.org