Dynamical seasonal predictions with the COLA atmospheric model

Seasonal hindcasts were made using the COLA (Center for Ocean-Land-Atmosphe re Studies) model for 16 winter seasons (mid-December through March for 198 1/82 to 1996/97). For each season, a nine-member ensemble was generated usi ng observed initial conditions in mid-December and observed global sea surf ace temperatures (SSTs). It is found that in the presence of large tropical SST anomalies, the model is quite successful in simulating seasonal-mean h eight anomalies over the Pacific-North America (PacNA) region. A local spatial pattern correlation field is computed for the ensemble seas onal mean of 500 hPa height from the General Circulation Model (GCM) and th e seasonal means from the re-analyses of the National Centers for Environme ntal Prediction, hereafter 'observations'. This field exceeds 0.6 over the eastern tropical Pacific and western North America; the maximum values are greatly enhanced during El Nino Southern Oscillation (ENSO) years. Similar results are obtained for 200 hPa winds. The observed enhancement of intrase asonal low-pass (10-90 day) variability of 500 hPa height during cold event s is simulated, as is the shift in the storm tracks (2-10 day variability o f 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 Oscilla tion is removed from both GCM and observations. The leading EOF mode explai ns about 50% of the variance for both GCM and observations; the two pattern s 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 th e GCM and observations. The great similarity between the height patterns in the first SVD mode and the corresponding leading mode EOF pattern indicate s that the leading height variations are forced by the 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) record of observed heights and SSTs. The pattern of percentage va riance explained in the GCM and in the observations are very similar to eac h other land to the EOF described above). The higher absolute values in the GCM case reflect the effectiveness of the ensemble in filtering out variab ility unrelated to SST forcing. SVD analysis was applied to 100 GCM 'samples' coupled with the observations ; a sample is defined as a single 16-year record obtained by picking one en semble member randomly for each of the 16 years. Probability distribution f unctions (p.d.f.s) of the pattern correlation for the leading SVD patterns, the percentage explained squared covariance, and the time series correlati ons all indicated sharp peaks at values of 0.87, 87%, and 0.82 respectively . The p.d.f. 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 p.d.f. 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 corr elation coefficient (ACC) for each of 36 possible intra-ensemble pairs. For eight years in which histograms of the ACC indicate predominantly positive values, the ACC of the ensemble mean with the observed seasonal mean is al so relatively high. Brier skill scores and reliability diagrams were computed for the 'event' o f the 500 hPa height being one standard deviation above (or below) the norm al, with all such events pooled over the entire: northern hemisphere or Nor th America only. All skill scores are positive and statistically significan t at the 99% level. The North American scores are higher than the whole hem isphere scores; the ENSO year scores are higher than those for all years.