A Dynamical One-Month Lead Seasonal Rainfall Prediction for
March to May 1997 for the North-eastern Area of South America
contributed by Mike Harrison1, Tony Evans1,
Ruth Evans1, Mike Davey2, and Andrew Colman2
1NWP Division 2Ocean Applications Branch
UK Meteorological Office, Bracknell, United Kingdom
One objective of the European PROVOST experiment (PRediction Of climate
Variations On Seasonal and interannual Timescales) is to estimate potential
dynamical seasonal predictability given ideal surface boundary conditions
on a global scale. To this end, four European models (the UKMO Unified
Model at climate resolution, the ECMWF T63 model and the ARPEGE model run
at T42 by Mätäo-France and T63 by EDF) have been integrated in
9-member ensembles initialized at 24-hour intervals for four months for
each season over 15 years from 1979 to 1993 (climatologies are calculated
over the same period). Common initial conditions and verifying analyses
obtained from the ECMWF reanalysis as well as common SST anomalies from
the UKMO GISST and Reynold's OI data sets were used in all experiments.
All initializations were at 0000Z finishing on the day prior to the start
of the season. The north-eastern region of South America has been identified
through these PROVOST experiments using the UKMO model as an area of relatively
high predictability and in this paper a real-time seasonal forecast at
one month lead is provided for the region based on these assessments.
The area with higher predictability broadly stretches through much of
eastern Brazil, the Amazon Basin, French Guiana, Surinam and Guyana, even
covering the southeastern parts of the Caribbean. Within this area correlations
over the 15 years between ensemble-mean rainfall anomalies and anomaly
values obtained from the gridded observed land-surface rainfall data set
of Hulme (1994) exceed 0.5, with correlations in excess of 0.8 in the vicinity
of the Nordeste region of Brazil (Fig. 1). The observed data set is gridded
to the same resolution as the model (2.5° x 3.75°), but in order
to reduce noise 4° x 4° blocks have been joined together to produce
Figure 1; only blocks with adequate data were retained. Time series of
en-semble mean rainfall and the Hulme data for selected regions illustrate
the fact that the dynamical model, although able to capture the interannual
variability reasonably well, has insufficient variability (Fig. 2). Hence
a variance inflation has been calculated using en-semble means for each
gridded area and applied to the forecasts from both the members and the
ensemble means given below. Also shown in Figure 2 are correlations between
ensemble mean predictions and some of the common rainfall indices for the
Nordeste (see Colman et al. 1997 for index details); again the model has
a high level of predictability (correlations approach 0.9) but requires
variance inflation (here done against the FQ index, but results for other
indices are almost identical). Note that observed anomalies frequently
lie within the ranges of the inflated ensembles or are close outliers to
those ranges.
Forecasts as produced for the 1997 March to May season are derived from
nine-member ensemble runs, but with two major design differences from the
predictability experiments outlined above. First, initialization for the
predictions is from late January rather than late February as in the PROVOST
runs. Secondly, the real-time experiments use persisted SST anomalies (from
January) throughout. Neither is thought likely to have significant negative
impacts on the model's ability to provide real-time predictions. Skill
over the region as deduced from the PROVOST runs remains high throughout
the year, whether for months 1 to 3 or 2 to 4 of the simulations. Indeed,
equivalent levels of skill tend to be present on a monthly time scale,
although with some drop-off into the fourth month. Hence the shift in the
start date is considered highly unlikely to affect potential predictability
for the region from that deduced for the "standard" seasons.
Concerning use of persisted SST anomalies, experiments for six winter seasons
have been carried out to date with persisted SST anomalies. While there
is some inevitable loss of predictability associated with the use of persisted
anomalies, this appears to be minimal in areas of relatively high predictability
such as considered here, and certainly does not eliminate predictability
in terms of the levels normally associated with seasonal forecasts (Fig.
3). Use of persisted anomalies fails, of course, during seasons in which
there is a substantial readjustment of SST anomalies over ocean areas related
to a given region's rainfall; experience has been gained of such failed
forecasts for the Nordeste in preliminary work with the model. Currently
there is no solution to this problem of rapid intraseasonal SST anomaly
distribution changes: the forecasts given below are conditional on the
assumption of continuity of the January anomalies.
The ensemble provides a consistent prediction of above-average rainfall
during March to May 1997 across the region except over the continental
interior (area E - below average in all members) and over the Nordeste
(area D), where five members are below-average (Fig. 4 and Table 1). A
southward displacement of the ITCZ from its model climatological location
is suggested. Inflated ensemble-mean, together with maximum and minimum,
rainfall anomalies for each of the gridded areas depicted in Fig. 1, plus
for various Nordeste rainfall indices, are listed in Table 1. Note that
the gridded prediction for the Nordeste is somewhat less than for the indices,
mainly because the block includes signal from the oceanic dry area representing
the shift in ITCZ location (Fig. 4).
For the first time since this research began, the dynamical model forecasts
for the Nordeste rainfall indices are entirely inconsistent with those
from the empirical techniques developed at UKMO (Colman et al. 1997; this
issue), techniques with an extended history of high skill. Examination
of the empirical methods indicates that most of the information this season
is being extracted from the Atlantic dipole, with a sign associated with
dry Nordeste conditions. An ENSO predictor is also included in the empirical
model, but this is currently weak. Scrutiny of Atlantic SST anomalies in
January 1997 reveals a region of warm anomalies extending eastward from
the Nordeste/Bahia coastal regions. Brief analysis suggests that anomalies
of this sign in this region are typically associated with above-average
Nordeste rainfall. It is hypothesized, therefore, that the difference in
the empirical and dynamical predictions may result from the dynamical model's
ability to respond to SST anomalies off Bahia whereas the UKMO empirical
methods are tuned only to basin-scale anomalies and hence do not give much
weight to the relatively localized effects from the Bahia region.
Colman, et al., 1997: Multiple regression and discriminant analysis
predictions of Mar-Apr-May 1997 rainfall in northeast Brazil. Experimental
Long-Lead Forecast Bulletin, Vol. 6, No. 1 (this issue).
Hulme, M., 1994: Validation of large-scale precipitation fields in general
circulation models. Global Precipitation and Climate Change, M.
Desbois and F. Desalmand, Eds., NATO ASI Series, Vol. 23, Springer-Verlag,
387-406.
Table 1. March to May 1997 seasonal forecast rainfall percentages of
normal for the unmodified Ensemble Mean (E Mean - with respect to the model
1979-1993 climate) and inflated Mean (I Mean - with respect to the Hulme
data for 1979-1983) and for the highest and lowest (inflated) ensemble
members for each of the 6 areas depicted in Figure 1. Also shown are equivalent
predictions for three Nordeste rainfall indices (FQ - Fortaleza-Quixeramobim;
N -Nobre; H - Hastenrath; see Colman et al., 1997), inflated using the
FQ index, estimated from the closest two and four model grid points; the
equivalent Hulme area is D. Note that the rainfall indices refer to different
periods although all have been predicted using March to May model rainfall
(FQ - March-May; Nobre - February-May; Hastenrath - March-April); thus
the Nobre prediction is not at long lead.
Forecast | A | B | C | D | E | F | 2 Grid Points | 4 Grid Points | ||||||||
(%) | FQ | N | H | FQ | N | H | ||||||||||
E Mean | 172 | 224 | 121 | 101 | 87 | 142 | 121 | 121 | 121 | 138 | 138 | 138 | ||||
I Mean | 299 | 377 | 146 | 103 | 57 | 198 | 146 | 145 | 145 | 166 | 149 | 153 | ||||
I Highest | 344 | 429 | 215 | 118 | 99 | 273 | 164 | 160 | 164 | 188 | 178 | 183 | ||||
I Lowest | 188 | 324 | 138 | 85 | 23 | 126 | 130 | 128 | 130 | 134 | 127 | 131 |
Fig. 1. Correlations for hindcasts over March
to May 1979-93 between ensemble mean rainfall and the Hulme gridded rainfall
set over 10 x 15 degree blocks.
Fig. 2. Time series of Ensemble Mean (EM) rainfall
(as % of normal) pre- and post-inflation of ensembles created using observed
SST and of the Hulme dataset for selected representative areas (see Fig.
1) and for the FQ Nordeste rainfall index. Results for the other indices
are similar. Bars indicate inflated range of ensembles and the median.
Capture Rate (CR) indicates the number of years out of 15 when the observations
lies within the range of the inflated ensemble.
Fig. 3. Ensemble mean January to March one
month lead hindcasts as % of normal for Area D (Fig. 1) for ensembles using
observed and persisted SSTs from 1989 to 1994. Correlation between observed
and persisted SST hindcasts is 0.78. Results for the other regions and
the Nordeste indices are similar. Symbols: circle-observed SSTs; triangle-persisted
SSTs; square-Hulme observed anomaly.