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:
 
 
VeryDry/Dry Dry/Average Average/Wet Wet /Very Wet
74% 86% 102% 124%

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
AGCM simulation skill correlation=0.65

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.
 
Year 1994 1995 1996 1997 1998 1999 2000
Forecast
Category

Very Wet

Dry

Average

Wet
Dry or
Average

Dry
Dry or Average
Observed
Category

Wet

Dry

Very Dry 

Very Wet

Very Dry

Average
Average

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:
 
 
Very Dry Dry Average Wet  Very Wet
0.33 0.34 0.30 0.03 0.00

 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:
 
Very Dry Dry Average Wet  Very Wet
0.55 0.22 0.13 0.10 0.00

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:

  • Correlation between independent hindcasts and observations over 1949-1998 are significant at the 5% level (shown in figure 1)
  • Correlation between independent hindcasts of this years forecast tercile and observations during 1949-1998 are significant at the 5% level

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.
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.
 

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 A2:
 
 


 
 

Figure A3: