Application of the El Nino-Southern Oscillation CLImatology and PERsistence (CLIPER) Forecasting Scheme

contributed by Christopher W. Landsea¹ and John A. Knaff²

¹NOAA/AOML/Hurricane Research Division Miami, Florida

²Cooperative Institute for Research in the Atmosphere Colorado State University Fort Collins, Colorado

To provide a baseline of skill in seasonal ENSO forecasting, a multiple regression has been used to take best advantage of CLImatology, PERsistence and trend of initial conditions the ENSO-CLIPER (Knaff and Landsea 1997). This replaces simple persistence as a skill threshold. "Skill" is then redefined as the ability to outforecast the ENSO-CLIPER - a more difficult task.

This statistical prediction method is based entirely on the optimal combination of persistence, month-to-month trend of initial conditions and climatology. Multiple least squares regression is employed to test a total of fourteen possible predictors for the selection of the bestpredictors, based upon 1950-1994 developmental data. A range of zero to four predictors were chosen in developing twelve separate regression models, developed separately for each initial calendar month. The predictands to be forecast include the Southern Oscillation (pressure) Index (SOI) and the Nino1+2, Nino 3, Nino 4 and Nino 3.4 SST indices for the equatorial eastern and central Pacific at lead times ranging from zero seasons (0 - 2 months) through seven seasons (21 - 23 months). Though hindcast ability is strongly seasonally dependent, substantial improvement is achieved over simple persistence wherein largest gains occur for two to seven season (6 to 23 months) lead times. The ENSO-CLIPER model thus not only offers a baseline "no-skill" forecast of ENSO variability, but a practical forecast based upon the CLIPER premise.

ENSO-CLIPER captures the climatological aspects (Huschke 1959, WMO 1992) of the whole ENSO complex capturing both mean conditions and propagation of those features in time. In essence, this model given initial conditions of ENSO (SSTs in Nino regions 1&2, 3, 3.4, and 4, and the SOI) and the recent past valid at a particular time will fit, using regression techniques, the best evolution from those initial conditions. The method has been frozen following its development (42 years), and yields the mean climatological evolutions for that period (1951-1992 to1953-1994 depending on lead-time). This procedure is analogous to statistical tropical cyclone track forecasting using a CLIPER approach (see Neumann 1972).

The regression design called leaps and bounds (IMSL 1987) is used to develop optimal models (the best subsets of a prescribed number of predictors). Predictors include 1, 3 or 5-month averages of initial predictor anomalies as well as their recent trends. Predictors are the predictands themselves at earlier times. Some limits on predictor selection were imposed to reduce overfitting (Aczel 1989). Skills are degraded from dependent sample results to reflect estimated independent forecast skill following Davis (1979) and Shapiro (1984). Final skill estimates reflect levels comparable to those of more sophisticated statistical and dynamical models. More details about the ENSO-CLIPER model, including its skill and its predictor selection rules, are given in the June 1997 issue of this Bulletin (p. 55). The program to run ENSO-CLIPER is also available upon request.

Employing the chosen predictors in the ENSO-CLIPER model on a 1 December 2000 initialization date yields forecasts for Dec.-Jan.-Feb. 2000 (0 season lead) out through Sep.-Oct.-Nov. 2002 (7 season lead). Results for the SSTs in Nino 3.4 and Nino 3 as well as the SOI are shown in Fig.1. These forecasts indicate that mild La Nina to neutral may continue through early spring of 2001, but there are large regional differences in the evolution of SST anomalies. This coming spring, ENSO-CLIPER develops warming in the western equatorial Pacific (Nino 4 and Nino 3.4) and cooling in the eastern Pacific (Nino 3 and Nino 1&2). All regions show warming from summer 2001 to winter 2001-2002. A return to an El Nino condition - perhaps a moderate event - is being forecast in western regions to begin in Sep.-Oct.-Nov 2001 and peak in Dec.-Jan.-Feb 2001-02. The forecast of conditions in Nino 3 produces a similar warming, but SST anomalies remain negative or neutral through next winter. The short leads are based primarily upon persistence of initial conditions and recent trends in initial conditions, while the change to El Nino-like conditions in over one year's time depend upon the climatological preference to change phase on the biennial timescale.

ENSO-CLIPER predictions made over the last several seasons have verified reasonably well, for lead times less than a year (Table 1, and past ELLFB issues). The longer leads are strongly biennial an thus have poorly forecast the persistent cool event. The 0 and 1 season lead forecasts in particular have performed very well, being within 0.6C of the actual anomaly all the way back to those issued in late 1998 with the exception of the DJF 1999-2000 verification. Over the last couple of years the performance of this model has been very competitive with both statistical and numerical ENSO forecast models. This fact suggests that there is little skill associated with the present ENSO forecast schemes (Landsea and Knaff 2000).

Acknowledgments: The authors wish to thank William Gray, Tony Barnston, John Sheaffer, Dave Enfield Dennis Mayer, Barb Brumit, Amie Hedstrom, Bill Thorson and Rick Taft for all their help and comments concerning this work. The second author is being supported by NOAA under contract NA67RJ0152 and is employed with the Regional Mesoscale Meteorology Team at CIRA.

References:

Aczel, A. D., 1989: Complete Business Statistics. Richard D. Irwin, Inc., 1056 pp.

Davis, R. E., 1979: A search for short-range climate productivity. Dyn. Atmos. Oceans, 3, 485-497.

Huschke, R. E., 1959: _Glossary of Meteorology. Second Printing, American Meteorological Society, 45 Beacon St., Boston, MA, 02108, 638 pp.

IMSL, 1987: FORTRAN subroutines for statistical analysis. International Mathematical & Statistical FORTRAN Library, 1232 pp.

Knaff, J. A. and C. W. Landsea, 1997: An El Nino-Southern Oscillation CLImatology and PERsistence (CLIPER) forecasting scheme. Wea. Forecasting, 12, 633-652.

Landsea, C. W., and J. A. Knaff, 2000: How much skill was there in forecasting the very strong 1997-98 El Nino? Bull. Amer. Meteor. Soc. 81, 2107-2119.

Neumann, C. J., 1972: An alternative to the HURRAN tropical cyclone model system. NOAA Tech Memo. NWS SR-62. 22 pp.

Shapiro, L. J., 1984: Sampling errors in statistical models of tropical cyclone motion: A comparison of predictor screening and EOF techniques. Mon. Wea. Rev., 112, 1378-1388.

WMO, 1992: International Meteorological Vocabulary. Second Edition, Secretariat of the World Meteorological Organization, Geneva, Switzerland, 784pp.

Figure 1: Forecast of Nino 3.4, Nino 3 and SOI using data available through 1 June 2000. Forecasts are valid for Jun.-Aug. (JJA) 2000, Sep.-Nov. (SON) 2000, Dec.-Feb. (DJF) 2000-01, Mar.-May (MAM) 2001, JJA 2001, SON 2001, DJF 2001-02 and MAM 2002. Actual numerical forecast values for these times are shown on each plot along with estimated RMSE bars. These anomalies are based upon a 1950-1979 mean.

Table 1: Recent history of ENSO-CLIPER forecasts and corresponding observations for the Nino 3.4 SST region (in degree C).

Target Period	Forecast Made 1 Jun 1999	Forecast Made 1 Sep 1999	Forecast Made 1 Dec 1999	Forecast Made 1 Mar 2000	Forecast Made 1 Jun 2000	Forecast Made 1 Sep 2000	Forecast Made 1 Dec 2000	Observed Anomaly
JJA99	-0.4							-0.8
SON99	-0.4	-0.7						-0.9
DJF99-00	-0.5	-0.9	-0.9					-1.6
MAM00	-0.1	-0.1	-0.3	-0.7				-0.7
JJA00	0.1	0.1	0.5	-0.1	-0.4			-0.3
SON00	0.5	0.2	0.2	-0.3	-0.3	-0.1		-0.5
DJF00-01	0.7	0.8	0.7	-0.3	-0.1	-0.7	-0.5
MAM01	0.4	0.3	0.4	0.3	0.1	0.1	-0.6