EXPERIMENTAL FORECAST OF EAST AFRICAN RAINFALL FOR OCTOBER-DECEMBER 2001
by Andrew Colman
Seasonal Forecasting Group, Ocean Applications, Met Office, Bracknell, UK
INTRODUCTION
The Met Office is conducting research into the effects of sea surface temperatures and other climatic variables on tropical rainfall. As part of this research, experimental seasonal rainfall forecasts have been made for the Sahel and adjacent regions in tropical NW Africa since 1986, and for the Nordeste region of Brazil since 1987. Using similar statistical methods, forecasts for tropical East Africa October-December rainfall (the 'short rains') have been issued since 1994 and appear in previous September issues of this bulletin .
A long-lead forecast for East African rainfall using observed data up to mid August has already been produced and was contributed to the Greater Horn of Africa Climate Outlook Forum ( GHACOF8 ). This forecast uses observed data up to the end of August.
The region covered by the East Africa prediction is between 5N and 15S and between 30E and the Indian Ocean coast. These forecasts for E Africa were produced using statistical methods and by using Met Office's Atmospheric General Circulation Model (AGCM). This year we have introduced forecasts for 2.5o latititude x 3.75 o longitude rectanglar grid box regions. Skill is not so high at this higher resolution but these forecasts give an indication of rainfall distribution.
The statistical forecast is made by using linear regression and discriminant analysis techniques, with three indices of global sea surface temperature (SST) anomaly patterns (Appendix, figures A1-A3 respectively). The forecast model is derived from historical rainfall and SST information.
The AGCM forecast was extracted from a nine member ensemble of AGCM predictions using sea temperatures and atmospheric conditions observed just prior to when the forecast was run (September 14th).
The historical rainfall record is divided into 5 equi-probable categories.
Based on 1961-1990 rainfall, the category boundaries (as percentages of mean rainfall) are: FORECAST SKILL PERFORMANCE OF TRIAL FORECASTS FOR 50 PAST YEARS The statistical and dynamical forecasts were tested using trial forecasts over the period 1948 to 1997. The assessment measure used is correlation. Statistical linear regression forecasts were assessed using a method where a trial forecast is made for each year using a regression equation
calculated using data for the remaining years. This assessment provides a good measure of forecast skill from minimal data. To provide an indication of AGCM skill, the performance of a long term AGCM run forced with observed SST in simulating rainfall is measured. Statistical forecast skill correlation=0.50 These correlations are statistically significant at the 5% level. PERFORMANCE OF REAL TIME EMPIRICAL FORECASTS Forecasts have been made for this region since 1994. The forecasts for 1994 and 1995 were strongly influenced by above and below average SST
in the NW Pacific respectively and the forecasts for 1997 and 1998 where influenced by the 1997 El Nino and the 1998 La Nina events. Note: The categories used for the 1994-1998 forecasts are based on a 1951-1980 climatology. For the 1999 and later forecasts, categories based
on the 1961-1990 climatology are used as 1961-1990 is the accepted WMO standard climatology period and is used by most forecasters. The
1961-1990 rainfall average is 104% of the 1951-1980 average. FORECASTS FOR THE 2001 SEASON STATISTICAL FORECAST Below average SST anomalies off the coast of Peru and in the tropical SW Atlantic are favouring below average rainfall in E Africa this year.
The regression forecast is 76% of the 1961-1990 average and is in the DRY category The discriminant analysis technique gives the following probabilities for the 5 (1961-1990 based) categories: AGCM DYNAMICAL FORECAST The ensemble mean prediction provides a best estimate rainfall forecast which is: 73% of the 1961-1990 model climatology. Based on the performance of AGCM ensemble simulations of rainfall from 1961 to 1990, the AGCM ensemble forecast is presented as
probabilities of 5 (1961-1990 based) observed rainfall categories which are: GRID BOX FORECASTS The grid box forecasts are expressed as probabilities of terciles which are climatologically equiprobable over 1961-1990. This is in order to make
the forecasts compatable with GHACOF forecasts which are expressed in the same way. Figure 1 shows the skill of the empirical forecasts. The
empirical and dynamical forecasts are shown in figures 2 and 3 respectively. Two forecast maps are shown for each category, one includes all grid boxes for which there is data, the second (skill mask) version includes only
gridboxes where independent test correlation skill is significant. To be included on the skill mask map, hindcasts for the box must pass at least 1
of these 2 tests: FIGURE 1: CORRELATION SKILL OF EMPIRICAL REGRESSION FORECASTS FIGURE 2: PROBABILITY FORECASTS BY EMPIRICAL METHOD FIGURE 3: PROBABILITY FORECASTS FROM AGCM DYNAMICAL FORECAST OVERALL BEST ESTIMATE: This year, the AGCM and the empirical forecasts are both indicating below average rainfall to be likely,
with the dynamical forecast being slightly drier than the empirical. REFERENCE: Mutai, C.C., Ward, M.N and Colman, A.W. Prediction of East Africa seasonal "short rainfall" rooted in evidence for widespread SST-forced
variability during October-December. I.J.Climatol.18 975-997 (1998). ACKNOWLEDGEMENTS: Thanks to David Rowell for providing output from the HADAM3 model. Thanks to Mark Naylor, Pete Mclean and Richard Graham for
supplying dynamical forecast output. APPENDIX : Predictor patterns used for empirical Forecast The pattern shown in figure 3 is the most important predictor contributing to over 50% of the forecast variance. Figure A1: Figure A3:
VeryDry/Dry
Dry/Average
Average/Wet
Wet /Very
Wet 74%
86%
102%
124%
AGCM simulation skill correlation=0.65
Year
1994
1995
1996
1997
1998
1999
2000 Forecast
Category
Very Wet
Dry
Average
WetDry or
Average
DryDry
or
Average Observed
Category
Wet
Dry
Very
Dry
Very
Wet
Very Dry
AverageAverage
Very Dry
Dry
Average
Wet
Very Wet 0.33
0.34
0.30
0.03
0.00
Very Dry
Dry
Average
Wet
Very Wet 0.55
0.22
0.13
0.10
0.00
Hence our best estimate is that rainfall will be in the DRY or VERY DRY category. The highest probabilities for below aberage rainfall
are for the east of the Region.
Figure A2: