On-line version ISSN 1996-7489
S. Afr. j. sci. vol.105 n.5-6 Pretoria May./Jun. 2009
Instituto Nacional de Pesquisas Espaciais, Caixa Postal 515, 12245-970 São José dos Campos, SP, Brazil. E-mail: firstname.lastname@example.org
The precipitation data for some regions in South Africa were studied for the period 19001998. From the 11 regions examined, 8 in South Africa had maximum precipitations in the austral summer months (December, January, February, March), while 3 had maxima in autumn and winter. Annual values showed considerable year-to-year fluctuations (50% to 200% of the mean), while five-year running means showed long-term fluctuations (75% to 150% of the mean). A spectrum analysis indicated periodicities in the ranges 23 (quasi-biennial oscillation, QBO), 34 (quasi-triennial oscillation, QTO), 611, 1721, 2326, 3235 and 5566 years, some common to, and some different in different regions. The QBO and QTO accounted for a substantial fraction (3050%) of the total variance. In five-year running means, the effects of QBO and QTO were suppressed considerably. The plots showed distinct peaks, but the spacings varied in a wide range, indicating that predictions based on extrapolation of single peaks are not likely to come true even for decadal averages. El Niño effects for the giant event of 1982/83 were as expected but those for 1997/98 were obscure, almost absent. Running means over 21 years did not indicate linear trends, upwards or downwards. Instead, considerable oscillations were seen, with magnitudes different in different regions (525%). On average, high values during 1915/16 decreased considerably (58%) up to 1935, oscillated upwards thereafter and recouped by 1980, but decreased considerably thereafter.
Key words: South Africa, rainfall, ENSO, periodicities, predictions, trends
Rainfall regimes in different parts of Africa are different, even in smaller regions, though some regional coherence is reported broadly in Sahel, East Africa and South Africa.1 African rainfall studies fall into several categories such as (a) teleconnections between regional rainfall anomalies and transport, (b) rainfall cycles (periodicities), (c) regional sea-surface temperature (SST) forcing, (d) relationships between atmospheric pressure and circulation patterns and atmospheric modelling, (e) relationship between 23-year or quasi-biennial oscillations (QBO) and 34-year or quasi-triennial oscillations (QTO), and (f) trends and predictions.
(a) Teleconnections between regional rainfall anomalies and transport
Some major rainfall anomalies are reported to occur simultaneously throughout Africa2,3 and global teleconnections have been invoked to explain the similarity in rainfall anomalies.4 Nicholson2 observed a certain degree of interhemispheric connectivity, but the strongest spatial associations were reported within the northern hemisphere (Sahel) and within equatorial Africa and South Africa. There seemed to be two preferred configurations, namely, precipitation of the same sign over all the African continent, and opposite signs between tropical and subtropical latitudes. Tyson5 showed that the Sahel and Kalahari had synchronous wet and dry spells while equatorial latitudes and South Africa had opposite rainfalls. Washington and Todd6 explored the tropical-temperate links in southern Africa and the southwest Indian Ocean daily rainfall, while Todd et al.7 reported on water vapour transport associated with tropical-temperate trough systems over southern Africa and the southwest Indian Ocean.
(b) Rainfall cycles (periodicities)
For South Africa, Dyer8 noted regional differences and assigned rainfall stations into homogeneous groups. Dyer and Tyson9 delimited regions of South Africa where specific oscillations were predominant and showed that a quasi 20-year oscillation was predominant in the summer-rainfall region covering the northeastern half of South Africa. This oscillation was used by them for estimating possible future rainfall conditions. Earlier, Tyson et al.10 showed that a 10-year to 12-year oscillation affected the southern coastal region and the adjacent inland area predominantly, while Dyer11 and Dyer and Gosnell12 attempted to link the 10-year to 12-year oscillation to the 11-year sunspot cycle. Vines13 examined data in the southern and northeastern part of South Africa and reported periodicities of 67, 1012 and 1620 years. Nicholson and Entekhabi14 found peaks in four bands at 2.22.4, 2.62.8, 3.33.8 and 5.06.3 years, common throughout equatorial and southern Africa but only weakly evident in northern Africa. (These authors did observe peaks in the 1625-year range, but dismissed these on the grounds that the resolution at low frequencies, i.e. high periodicities, was poor.)
(c) Regional SST forcing
A relationship with the El Niño-Southern Oscillation (ENSO) was reported for eastern and southern Africa.1517 Janowiak17 reported that positive rainfall departures in equatorial East Africa and negative departures in South Africa followed ENSO events. Nicholson and Kim18 made a comprehensive study of the rainfall response to ENSO episodes over Africa. The strongest signals were reported to have appeared in southern, eastern and far northern Africa and the weakest in the Sahel. Nicholson and Entekhabi19 studied the relationship between rainfall in equatorial and southern Africa and the SST along the southwestern coast of Africa. When comparing the rainfall results with those of SST obtained in a previous study,20 Nicholson and Kim18 concluded that the ENSO episodes that influenced rainfall over Africa were those that manifested as SST fluctuations in the low-latitude Atlantic and western Indian Oceans. Mason2123 described the temporal variability of SST around South Africa. In recent years, Reason24 examined the relationship between subtropical Indian Ocean SST dipole events and southern African rainfall, while Reason and Rouault25 reported an ENSO-like decadal variability in South African rainfall. Rouault et al.26 examined the relationship between southeast tropical Atlantic warm events and southern African rainfall. Reason et al.27 reported results of interannual variability in rainfall over the Limpopo Province of South Africa (see also Kruger28).
(d) Relationships between atmospheric pressure and circulation patterns and atmospheric modelling
For South Africa, statistical models for seasonal forecasts of southern Africa summer rainfall have been earlier invoked.2931 Recently, Engelbrecht32 reported simulations of climate and climate change with the conformal-cubic atmospheric model. Hewitson and Crane33 studied the general circulation model (GCM) climate change projections with empirical downscaling.
(e) Relationship between QBO and QTO
A factor complicating the analysis of ENSO and solar activity effects on South African rainfall is their apparent, unexplained modulation by a QBO. Earlier, Labitzke and van Loon34,35 and van Loon and Labitzke36 reported that the atmospheric response to solar activity may be governed by the phase of the QBO as measured by the stratospheric zonal wind component over the equator. Mason and Tyson37 explored this idea for southern African rainfall and showed that solar activity effects on southern African rainfall were modulated by the initial QBO phase (easterly or westerly), and that, although solar forcing is not always a dominant effect, its contribution cannot be discounted and depends on the phase of the QBO (for further details and review see Mason and Jury38 and Mason and Tyson39).
(f) Trends and predictions
Kruger and Shongwe40 studied the temperature trends in South Africa during 19602003, while Kruger41 reported on trends in daily precipitation indices in South Africa during 19102004. Tadross et al.42 made regional climate model (RCM)-based projections of changes in southern African summer climate.
This article reports on the analysis of rainfall cycles and trends and ENSO effects over 11 regions in South Africa, and also, in part, revisits some of the earlier predictions.
Precipitation data reported are the monthly precipitation values from the high-resolution gridded (5° × 5°) data sets,4345 available on the website www.cru.uea.ac.uk/cru/data/hrg. The climate surfaces were constructed from a new data set of station climatological normals, numbering 19 800 for precipitation. The station data were interpolated as a function of latitude, longitude, and elevation using thin-plate splines. The accuracy of the interpolations was assessed using cross-validation and by comparison with other climatologies (more details are given in New et al.43). South African data were available for the locations A (19001997), B, C (19001998), D (19001997), E, F, G (19001998), H, I (19001996), J (19561998) and K (19491998), with geographical coordinates as shown in Table 1, and on a map of South Africa in Fig. 1.
Data for ENSO parameters, namely the Southern Oscillation Index (SOI) represented by the Tahiti (T) minus Darwin (D) atmospheric pressure difference (TD), and SST of the Niño 1+2 region in the east Pacific (010°S, 9080°W) were obtained from the Climate Prediction Center, National Oceanic and Atmospheric Administration website (www.cpc.ncep.noaa.gov/ data/indices/). Data for global temperatures46 were obtained from the website http://lwf.ncdc.noaa.gov/oa/climate/research/ anomalies/
Climatology of monthly rainfall totals
Figure 2 shows the climatology (average monthly precipitation values, 19001996) for these locations. The rainfall data available from the website (www.cru.uea.ac.uk) are given in units of (mm × 10), probably to avoid a decimal point. Hence, in Fig. 2, the ordinate is in units of (mm × 10) and the numbers indicating monthly average rainfalls (horizontal lines, not annual rainfall) are also in units of (mm × 10). Thus, in region A, the January rainfall was 1 483 (mm × 10), i.e. 148.3 mm, the monthly average was 580 (mm × 10), i.e. 58 mm, while the annual rainfall (not shown here) for region A was 6 961 (mm × 10), i.e. 696 mm. As noted by earlier researchers,5,37 the mean rainfall differs considerably from region to region; regions AG have peak rainfall in the southern hemisphere summer months (mostly December to March), region H seems to have very little seasonal variation (all season rainfall), and regions IK have peak rainfall in winter (MayJune). Rainfall in regions J and K is very high (J: monthly average ~270 mm, annual ~3 240 mm; K: monthly average ~196 mm, annual ~2 352 mm) compared to the other regions (monthly averages 9.593.1 mm, annual 1141 117 mm).
Rainfall cycles and trends
Annual values of every series were expressed as percentages of the overall (19001996) mean of that series. Figure 3a shows the plots (percentage of average monthly mean of each region) for series AI. (Again, the values are average monthly means in units of (mm × 10), thus 580 is actually 58 mmthe average monthly rainfall for region A. The annual value for region A (not shown) would be 696 mm.) Peaks are marked with dots and troughs with crosses. As can be seen, there are considerable year-to-year fluctuations (50200% of the mean), more in some series than in others. To suppress short-term fluctuations, five-year running means were calculated. These are shown in Fig. 3b. The plots are now smoother and the fluctuations are smaller (70150% of the mean), but the peak spacings (indicated by numbers) vary in a wide range (928 years). Similar peaks are joined by vertical lines and the peaks during 19741976 seem to be almost simultaneous. Dyer and Tyson9 reported that a quasi 20-year oscillation predominated in the summer rainfall region covering the northeastern half of South Africa and estimated above-normal rainfall for the intervals 19721981 and 19912000, and drier conditions for the in-between interval 19811990. In Fig. 3b, 19741976 had above-normal rainfall and 19801984 had below-normal rainfall, but the predicted above-normal rainfall during 19912000 does not seem to have occurred, at least until 1996 when the data series ends. However, it seems that whereas 19971999 also had low or normal rainfall, 2000 was a highly anomalously wet year with devastating floods, so much so that Mason47 commented that a return to a wet phase appeared to have occurred. So, in a restricted sense (for only one year during 19912000), the prediction of Dyer and Tyson seems to have come true.
Recently, Reason and Rouault25 reported the connection of ENSO-like decadal variability with South African rainfall. The results reported above are not directly comparable to theirs, because they are annual values whereas they studied summer- and winter-rainfall series separately and their main emphasis was on the comparison of the smoothed summer and winter series with SST-MSLP (sea-surface temperaturemean sea level pressure). However, their plot of seven-point running means of summer-rainfall values showed maxima in years 1910, 1920, 1937, 1941, 1956, 1975 and 2000, yielding spacings in successive maxima of 10, 17, 4, 15, 19 and 25 years (range 425 years), by no means near-decadal.25 Their winter-rainfall maxima were in 1905, 1914, 1942, 1954, 1974 and 1990, yielding spacings in successive maxima of 9, 28, 12, 20 and 16 years (range 928 years), again by no means near-decadal.
Thus the spacings in their plots were just as erratic as those shown in Fig. 3b (928 years). In any case, periodicities (decadal or otherwise) can be ascertained accurately only by a spectral analysis, as is shown in the following section.
The time series of annual values in Fig. 3a were subjected to spectral analysis. The method used was the Maximum Entropy Method48,49 (MEM), which locates peaks more accurately than the conventional Blackman and Tukey50 (BK) method. However, the amplitude (power) estimates in MEM are not very reliable.5153 Hence, MEM was used only for detecting all the possible peaks Tk (k = 1 to n), using length of the prediction error filter (LPEF) as 50% of the data length. These Tk were then used in the expression:
where f(t) is the observed series and E the error factor. The amplitude rk is obtained as a square root of the sum of the squares of ak and bk, while the phase øk is obtained as tan(inverse) of the ratio (bk/ak). A multiple regression analysis (MRA)54 was then carried out to estimate Ao (ak, bk), and their standard errors (by a least-squares fit). From these, amplitudes rk and their standard error sk (common for all rk in this method, which assumes white noise) were calculated. Any rk exceeding 2s is significant at a 95% (a priori) confidence level. This simple method will be termed as MEM-MRA.
Figure 4 shows the spectra (amplitudes versus periodicities detected by MEM-MRA of each series). As can be seen, strong periodicities are in the QBO (23 year periodicity) region. These contribute about 3050% of the total variance. Other periodicities are near 3.5, 4.2 (QTO, 34-year periodicity), 6.5, 11 (period of the sunspot cycle?), 20, 35 and 60 years, but none of these contributes more than 10% variance. That should be the reason why predictions based on the extrapolation of any one of these may not come true.
Reason and Rouault25 supposedly indicated that the 20-year periodicity is of great importance. This is not true. As reported above, their summer and winter rainfall data are shown as seven-point running means, hence the QBO, QTO and the 6.5-year periodicities are negated. In the smoothed plots too, the 20-year wave is not at all evident. In their plot of summer rainfall, the spacings in successive peaks are 10, 17, 4, 15, 19 and 25 years (a wide range of 425 years, average 15 years), and in the winter plot, the spacings in successive peaks are 9, 28, 12, 20, 16 years (a wide range of 928 years, average 17 years). Thus, a striking 20-year periodicity does not exist, and if it did, it would be closer to 1516 years. The values for years 19002000 of these plots were subjected to a spectral analysis using MEM-MRA. There was no 20-year periodicity (or for that matter, any other periodicity) significant at a 2s level, either in the summer or in the winter series.
Quasi-biennial and quasi-triennial oscillations
Since a large fraction of the total variance of any of these series is in the QBO and QTO regions, a detailed examination of the data at a fine resolution (say three-monthly) would be interesting. Figure 5 shows the plots for the restricted 20-year interval 19781998. Figure 5 (top) shows the three-monthly rainfalls (JFM, AMJ, JAS, OND, four values per year) for region A. Some QBO and QTO are seen but with much scatter. Hence, running means were calculated for four successive three-monthly values. The smoothed three-monthly values for region A are shown in the second plot. Now, smooth QBO and QTO waves are seen, but the spacings between successive peaks (indicated by numbers, in months) vary in a wide range (1551 months), indicating that the QBO and QTO are irregular (almost random) and beyond the reach of predictions. The other plots show the smoothed three-monthly running means for the regions BK. All show irregular QBO and QTO. Some peaks are similar for some regions (indicated by vertical lines), but the spacings are smaller in regions EK (many near 24 months) compared to spacings in regions AD (30 months or more), indicating considerable regional differences in QBO and QTO characteristics. This might have important implications for the origins of the QBO and QTO, which are most probably not random occurrences. The irregularity of the spacings is probably due to the presence of more than one QBO or QTO frequency.
El Niño effects
Trenberth55 mentions 'different flavours of El Niño'. During 19001998, there were several El Niño events, but two of them (19821983 and 19971998) were most prominent and occurred in the last few decades. Figure 6 shows the plots for some years before, during and after the El Niño of 19821983. The two top plots are for the Southern Oscillation Index represented by the Tahiti (T) minus Darwin (D) atmospheric pressure difference (TD), and SST anomalies of the El Niño 1+2 regions in the east Pacific (010°S, 9080°W). The El Niño interval from May 1982 to November 1983 is marked by thick vertical lines (the other vertical lines mark the month January, J) and the negative values of (TD) and the positive values of SST anomalies for this interval are shaded black. The other plots are de-seasoned (climatology subtracted) monthly precipitation deviations for the eleven regions AK as follows: de-seasoned means the average values for January, February etc. are calculated separately for the whole interval 19001998, thus yielding 12 average values for the 12 months JanuaryDecemberthis is climatology. Then, from the January value of every year, the average January value (climatology) is subtracted. The difference is considered the rainfall anomaly for January, and similarly for February to December. Thus, a continuous series of rainfall anomalies is obtained for each region, plotted as plots 311 for regions AK. (Note that the ordinate scale is mm × 10. Thus, 1 000 means 100 mm). These can be compared with any other parameters, such as SST anomalies. In Fig. 6, negative deviations are shaded black. The following may be noted:
(i) In region A (third plot), there are prominent rainfall decreases (droughts, shaded black) during October 1982 to April 1983, coinciding with the El Niño. However, there are similar decreases before the El Niño during November 1981 to May 1982 and after the El Niño during October 1983 to May 1984. Thus severe droughts occurred before and after the El Niño, indicating that there may be other physical processes causing droughts in the region.
(ii) Regions BH also show droughts of more or less extent during the El Niño, indicating that the 19821983 El Niño had a forcing effect on the rainfall of many regions of South Africa.
(iii) Regions AD (most of the summer rainfall region of South Africa) did not experience periods of significant above-normal rainfall during the 19821983 El Niño. In the remainder of the summer rainfall region of South Africa (regions E and G) and also region F over the Indian Ocean, anomalously high rainfall during the El Niño only occurred towards the end of the El Niño period. These results strengthen the general idea that El Niño is related with dry conditions over the summer rainfall region of South Africa.
(iv) Quiroz56 mentions the El Niño of 19821983 as a season of extraordinary climate anomalies. Most of the summer rainfall region of South Africa clearly experienced generally dry conditions during the El Niño.
These results using annual values are in general agreement with those of Reason et al.27 using summer values for the Limpopo region of northern South Africa, which state that El Niño 3.4 SST anomalies show a robust relationship with dry spell frequency during 19792002 (which includes the 19821983 event). On the other hand, the decadal relationship mentioned by Reason and Rouault25 is not very meaningful because they have used seven-point running means, so the major QBO and QTO effects in rainfall are lost. What remains is a mixture of insignificant higher periodicities (decadal range, if decadal means 10 and more, not just 10 and near). The SST-MSLP has very prominent variations in the QBO and QTO ranges (more in QTO, see details in Kane57), but Reason and Rouault25 have carried out an empirical orthogonal function (EOF) analysis and selected only the EOF representing the decadal to multidecadal band of SST-MSLP. Thus, whereas the main relationship between rainfall and SST-MSLP would be in the QBO and QTO regions, these have been negated in the decadal averaging and the results give the relationship between two comparatively insignificant series. It is no surprise then that their correlation was only about 0.50, which they claimed as significant at a 2σ level (probably true), but in my view is too low for considering as a meaningful relationship, and certainly is not useful for any meaningful predictions (a correlation of 0.50 implies a variance of only 0.25 or 25% of real relationship and 75% of random component).
Figure 7 shows similar plots for 19951998 where the El Niño interval (March 1997 to July 1998) is enclosed in thick vertical lines and negative (TD) values and positive SST anomalies are shaded black. (The ordinate scale is in units of mm × 10, such that 1 000 means 100 mm.) For some regions, data for 1998 or 19971998 were missing and only regions AC and HK are considered. Except for region A, where droughts seem to be associated with the El Niño, the El Niño effects are dubious, with positive and negative rainfall deviations spread throughout 19951998. Reason et al.27 reported somewhat similar results and some strange characteristics of the 19971998 El Niño are documented in Kane.58 An interesting aspect has been recently pointed out by Landman and Mason.59 They mentioned:
The El NiñoSouthern Oscillation (ENSO) signal in Indian Ocean sea-surface temperatures has weakened since the late-1970s, resulting in important changes in the association between tropical western Indian Ocean sea-surface temperatures and DecemberFebruary rainfall over South Africa and Namibia. Prior to the late-1970s, warm events in the tropical Indian Ocean occurred frequently in association with ENSO warm events, and both were typically associated with anomalously dry conditions over much of South Africa and Namibia. Over the most recent two decades, however, sea-surface temperature variability in the tropical western Indian Ocean has become significantly less dependent upon ENSO. At the same time, warm (cold) events in the tropical western Indian Ocean have become associated with wet (dry) conditions over the northeastern half of South Africa and northern Namibia since the late-1970s. These changes in rainfallsea-surface temperature associations can be successfully simulated using a general circulation model. Such changes in the climate system have important implications for the predictability of southern African rainfall, and highlight the need for an operational forecast system to predict Indian Ocean sea-surface temperatures.
In a recent review, Mason47 had the following to say:
Since the late-1970s, El Niño episodes have been unusually recurrent, while the frequency of strong La Niña events has been low. Prolonged/recurrent warm event conditions of the first half of the 1990s were the result of the persistence of an anomalously warm pool near the date line, which, in turn, may be part of an abrupt warming trend in tropical sea-surface temperatures that occurred in the late-1970s. The abrupt warming of tropical sea-surface temperatures has been attributed to the enhanced-greenhouse effect, but may be indicative of inter-decadal variability: earlier changes in the frequency of ENSO events and earlier persistent El Niño and La Niña sequences have occurred. Most forecasts of ENSO variability in a doubled-CO2 climate suggest that the recent changes in the tropical Pacific are anomalous. Of potential concern, however, is a possible reduction in the predictability of ENSO events given a warmer background climate.
The El Niño influence on rainfall over southern Africa occurs largely because of a weakening of tropical convection over the subcontinent. A warming of the Indian Ocean during El Niño events appears to be important in providing a teleconnection from the equatorial Pacific Ocean. The abrupt warming of the tropical Pacific and Indian oceans in the late-1970s is probably partly responsible for increasing air temperatures over southern Africa, and may have contributed to a prolongation of predominantly dry conditions. A return to a wet phase appears to have occurred, despite the persistence of anomalously high sea-surface temperatures associated with the late-1970s warming, and a record-breaking El Niño in 1997-1998.
It is obvious, therefore, that drastic imbalances have occurred in recent decades, and that ENSO and other relationships established before the 1970s are no longer valid. A new approach with changed inputs is necessary for estimating the future climate behaviour of South Africa, as is being done by Landman and Mason31 and others.
Rainfall trends over South Africa
In view of the multiple periodicities seen in the precipitation series, long-term trends will be difficult to identify. Some small periodicities can be suppressed by calculating running means over a fairly long interval, for example, 21 years. Figure 8 shows the plots of 21-year running means for (a) rainfall percentages for regions AI, (b) average rainfall for AI and for the whole of South Africa, and (c) global temperatures for land (full line) and sea (crosses).46 The following may be noted:
(i) None of these trends is monotonically upwards or downwards. Hence, fitting linear trends is meaningless.
(ii) The numbers are percentages and indicate long-term oscillations different in different regions, with ranges from 926%, with some spacings of 30 years or more.
(iii) After 1980, there seems to be a significant decrease in regions AG. (The plots are 21-year running means and the error is much less than 0.5%, while the decreases from the 1980 level to the 19891990 level are more than 5%).
(iv) The bottom plots (b) show the average percentages for regions AI and for South Africa as a whole (an independent estimate, given in the website of precipitation www.cru.uea. ac.uk). The two plots are very similar, show oscillations between 19141936, 19361948, 19481961, 19611967, 19671974 and 19741980, and a large decrease from 1980 to 1989, with magnitudes of 510%. Recently, Kruger41 reported that 'some significant changes in (precipitation) indices, averaged over specific areas in South Africa, could be identified. These include areas with significant increases and decreases in annual precipitation, increases in the longest annual dry spell indicating more extreme dry seasons, increases in the longest annual wet spells indicating more extreme wet seasons, and increases in high daily precipitation amounts. The conclusion is that, while in the largest part of South Africa there has been no real evidence of changes in precipitation over the past century, there are, however, some identifiable areas where significant changes in certain characteristics of precipitation have occurred over the period 1910 to 2004'. Kruger28 mentioned previously that for the period 19551991 (the latter half of Fig. 8), the 1960s and 1980s had below-normal rainfall while the 1950s and 1970s had above-normal rainfall, the turning points being 19591960, 19701971, 19801981 and 19911992. This is partly evident in Fig. 8, but the spacing of about 10 years (implying a 20-year full cycle) is not evident (the spacings of 19 and 14 years during 19481990 and 19101950 are equivalent to a long-term swing of about 33 years). However, Kruger28 made an interesting observation, namely, during the epochs of above-normal rainfall, the El Niño effects are less severe than those during epochs of below-normal rainfall (the opposite is true for La Niña). Kripalani and Kulkarni60,61 have also made a similar comment for the all-India summer monsoon rainfall.
(v) The bottom plot (c) of land temperature (full line) shows a steady increase during 19101944 (0.230°C to +0.074°C = +0.304°C), a slight decrease during 19441966 (+0.074°C to 0.037°C = 0.111°C), a slow increase during 19661974 (0.037°C to +0.029°C = + 0.066°C), and a rapid increase during 19741992 (+0.029°C to +0.469°C = +0.440°C). The bottom plot (c) of sea temperature (line of crosses) is almost the same as for land temperature except that from 1974 onwards, sea temperatures increased less (+0.073°C to +0.263°C = +0.190°C). In any case, the temperature patterns in (c) are not similar to those of the average rainfall patterns in (b). However, these are global temperatures, while a more relevant comparison would be with temperature trends over South Africa itself. Kruger and Shongwe40 have shown such temperature trends for South Africa but only for 19602003, so only 19601990 in Fig. 8 can be compared. Kruger and Shongwe40 report a positive trend in the annual average temperatures of 24 stations, such that a 21-year running mean would also show a smooth positive (upward) trend, just like the global trend during 19601990. Thus, neither the global temperature trend nor the localised South African temperature trend (both smoothly upward from 1960 onwards) match the fluctuations of rainfalls depicted in the plots (b).
Using the data available until the mid-1970s, Tyson et al.,10 apart from postulating regions of rainfall variability based on spectral analysis, concluded that there were no significant trends. Here, it is shown that there are long-term trends but that they are not linear. Rather, they fluctuate downwards in some intervals, upwards in others and decrease greatly after 1980.
Regarding alternating extended wet and dry spells during the 20th century in South Africa, the series of anomaly maps showing the temporal and spatial variability of annual rainfall in Tyson,5 updated in Tyson and Preston-Whyte,62 give a convincing illustration of the extent of the 20th-century variability patterns until the end of the 1980s. Obviously, the environment has changed considerably in recent decades.
General circulation model experiments have been suggesting that the unprecedented global warming in the 20th century may result in a more intense hydrological cycle, with an associated increase in the frequency and/or magnitude of heavy precipitation. Testing for South Africa, Fauchereau et al.63 mentioned that there was no significant trend in the Southern Africa Rainfall Index across the century, but parts of South Africa experienced a significant shift toward increasing probabilities of extreme rainfall. However, they also mentioned that the increase in rainfall variability over the latter decades of the 20th century was merely a return to similar conditions that existed in the early decades of the century, when it was significantly cooler than it is currently. They concluded that over the course of what climate alarmists call a century of unprecedented global warming, and contrary to the implications of nearly all 2 × CO2 GCM experiments, there has been no net change in either South African rainfall variability or the mean value of the Southern African Rainfall Index. The results shown in the bottom of Fig. 8 are slightly different. From a high level in 19151916, values oscillated downwards and decreased by about 510% but recovered in 1980 before decreasing considerably thereafter (~8%). This pattern is not similar to the pattern of land-sea global temperature variations.46
Using extrapolations of long-term periodicities, predictions have been attempted. Dyer and Tyson9 used a quasi 20-year oscillation and estimated that the intervals 19721981 and 19912000 would possibly have above-normal rainfall, while 19811990 would be drier than normal. Whereas estimates for 19721981 and 19811990 came true, the above-normal rainfall during 19911996 did not occur. (In 19971998, the strong El Niño should have caused below-normal rainfall, which did not occur; however, this cannot be interpreted as above-normal rainfall.) However, 2000 was a year of abnormally-high rainfall, so Dyer and Tyson's prediction is partially vindicated.
Predictions are made using two general methodsstatistical methods and dynamically-based (GCM) methods. Workers have often reported relationships between rainfall and parameters like Atlantic SST.24,26,64 Rouault et al.26 stated that, whereas warm events may be associated with unwelcome floods along the Angolan and Nambian coast, increased rainfall elsewhere can sometimes alleviate droughts in other regions of southern Africa. Recently, several workers reported seasonal forecasts for western, eastern and central Sahel, and for East Africa and South Africa. These involve integration of three European models (UKMO unified model, the ECMWF T63 model and the ARPEGE model) in a nine-member ensemble, neural networks, or multivariate regression analyses, where ENSO is only one of the parameters (not necessarily an important one) and global SST and/or SST in nearby oceans is generally an important parameter, the anomalies of which may or may not be related to Pacific SST anomalies. Several other parameters are also used. For South Africa, Jury29 used as many as 49 predictors (e.g. SSTs, upper winds, surface pressure and wind, QBO, satellite cloud cover, cloud albedo), one-third of which were related to ENSO, another one-third to transition of Indian monsoon and one-fourth to conditions in the south Atlantic. Nevertheless, the prediction skills are only modest to moderate. For example, for South Africa during the giant El Niño event of 19971998, Jury,29 Greischar et al.,65 Landman,66 and Thiaw and Barnston67 gave a forecast of below-average rainfall in the austral summer of 19971998, which turned out to be only partly true; rainfall was deficit in the west (Namibia), almost normal in the south, and excessive in the southeast of South Africa. For the 19981999 summers, Greischar et al.,65,69 Thiaw and Barnston68 and Landman and Mason70 predicted excess rainfall for South Africa, but droughts seem to have occurred there. Thus, these attempts cannot be considered as successful. We suspect that presently unpredictable, almost random, components are very prominent and vitiate the predictions.
The precipitation data for some regions in South Africa were studied for 19001998. The following was noted:
(1) From the 11 regions of South Africa examined, 8 had maximum precipitations in the summer months (January, February, March), whereas 3 had maxima in autumn and winter.
(2) Annual values showed considerable year-to-year fluctuations (50% to 200% of the mean).
(3) When smoothed by calculating five-year running means, long-term fluctuations were still large (70% to 150% of the mean).
(4) A spectrum analysis indicated periodicities in the ranges 23 (QBO), 34 (QTO), 611, 1721, 2326, 3235 and 5566 years, some common to, and some different in different regions.
(5) The QBO and QTO accounted for a substantial fraction (3050%) of the total variance. Since there was generally more than one peak in QBO and QTO, the net effect was almost a random-like (unpredictable) component. With such a large random component, statistical predictions of annual values to any reasonable degree of accuracy cannot be expected. (This does not apply to seasonal forecasts based on GCMs). In five-year running means, the effect of QBO and QTO got suppressed considerably. The plots did show distinct peaks, but the spacings varied in a wide range, indicating that predictions based on extrapolation of single periodicities (obtained in the spectral analysis) are not likely to come true even for decadal averages. Dyer and Tyson9 reported that a quasi 20-year oscillation predominated in the northeastern half of South Africa (probably regions A and B) and estimated that the intervals 19721981 and 19912000 may possibly have above-normal rainfall, while 19811990 may be drier than normal. Whereas the estimates for 19721981 and 19811990 came true, the predicted above-normal rainfall during 19911999 did not occur. However, 2000 was a year of abnormally-high rainfall, partially vindicating Dyer and Tyson's9 prediction.
(6) Regarding El Niño effects for the two giant events of 19821983 and 19971998most of the summer rainfall region of South Africa had droughts during 19821983. For 19971998, effects were obscure, almost absent.
(7) Running means over 21 years did not indicate linear trends, upwards or downwards. Instead, considerable oscillations were seen, different in different regions (525% of the mean). On average, high values during 19151916 decreased considerably (58%) up to 1935, oscillated upwards thereafter and recouped by 1980, but decreased considerably thereafter (210%, statistical uncertainty less than 0.5%).
In short-term fluctuations, the (quantitatively) strong influence of QBO and QTO is somewhat disconcerting. There is a plethora of periodicities even in narrow ranges of 23 years and 34 years and the peak separations are irregular. Random series do yield a plethora of periodicities, most of which are statistically insignificant or of borderline significance. Hence, a substantial part of the QBO and QTO appears to have an almost random nature, beyond the reach of predictions made a few years before a specific wet or dry spell. However, many peaks seemed to be similar in the different regions, indicating the possibility of at least a partial common physical cause. This needs further exploration. The long-term fluctuations involve many periodicities and except for the 11-year periodicity (possibly corresponding to the sunspot cycle), these do not have ready explanations. Also, since each explains only about 10% variance and stationarity is not guaranteed, predictions based on its extrapolation are likely to prove erroneous. El Niño effects, though prominent in earlier years, were mostly obscure for the giant event of 19971998, indicating that some drastic environmental changes have occurred in recent decades. Since 1960, the trends in temperatures were consistently upward, but rainfalls had fluctuations. Running means over 21-yearly values could not have data uncertainties, and yet, the trends even in nearby locations do not seem to be similar. Can local effects be so overpowering? It is likely that global changes may not keep up to the climatological expectations, particularly from model studies;71,72 but the changes may induce changes in the synoptic circulation over South Africa, which in turn may cause different trends over different regions.32,33,42
I thank Tim D. Mitchell (Tyndall Centre for Climate Change Research, U.K.) for useful suggestions. This work was partially supported by FNDCT, Brazil, under contract no. FINEP-537/CT.
1. Hulme M. (1992). Rainfall changes in Africa: 19311960 to 19611990. Int. J. Climatol. 12, 685690. [ Links ]
2. Nicholson S.E. (1986). The spatial coherence of African rainfall anomalies: interhemispheric teleconnections. J. Climate Appl. Meteor. 25, 13651381. [ Links ]
3. Nicholson S.E. (1986). The nature of rainfall variabilty in Africa south of equator. J. Climatol. 6, 515530. [ Links ]
4. Stoeckinus T. (1981). Interannual variations of tropical precipitation patterns. Mon. Wea. Rev. 109, 12331247. [ Links ]
5. Tyson P.D. (1986). Climate Change Variability in South Africa. Oxford University Press, Cape Town. [ Links ]
6. Washington R. and Todd M.C. (1999). Tropical temperate links in southern Africa and southwest Indian Ocean daily rainfall. Int. J. Climatol. 19, 16011616. [ Links ]
7. Todd M.C., Washington R. and Palmer P.I. (2004). Water vapour transport associated with tropicaltemperate trough systems over southern Africa and the southwest Indian Ocean. Int. J. Climatol. 24, 555568. [ Links ]
8. Dyer T.G.J. (1975). The assignment of rainfall stations into homogeneous groups: an application of principal component analysis. Quart. J. Roy. Meteor. Soc. 101, 10051013. [ Links ]
9. Dyer T.G.J. and Tyson P.D. (1977). Estimating above and below normal rainfall periods over South Africa, 19722000. J. Appl. Meteor. 16, 145147. [ Links ]
10. Tyson P.D., Dyer T.G.J. and Mametse M.N. (1975). Secular changes in South African rainfall: 18801972. Quart. J. Roy. Meteor. Soc. 101, 817833. [ Links ]
11. Dyer T.G.J. (1977). On the application of some stochastic models to precipitation forecasting. Quart. J. Roy. Meteor. Soc. 103, 177189. [ Links ]
12. Dyer T.G.J. and Gosnell J.M. (1978). Long-term rainfall trends in the South African sugar industry. S. Afr. Sugar Technol. Assoc. I, 18. [ Links ]
13. Vines R.G. (1980). Analyses of South African rainfall. S. Afr. J. Sci. 76, 404409. [ Links ]
14. Nicholson S.E. and Entekhabi D. (1986). The quasi-periodic behaviour of rainfall variability in Africa and its relationship to the Southern Oscillation. J. Climate Appl. Meteor. 34, 331348. [ Links ]
15. Dyer T.G.J. (1979). Rainfall along the east coast of Southern Africa, the Southern Oscillation and the latitude of the subtropical high pressure belt. Quart. J. Roy. Meteor. Soc. 105, 445452. [ Links ]
16. Lindesay J.A., Harrison M.S.J. and Haffner M.P. (1986). The Southern Oscillation and South African rainfall. S. Afr. J. Sci. 82, 196198. [ Links ]
17. Janowiak J. (1988). An investigation of interannual variability in Africa, J. Climate 1, 240255. [ Links ]
18. Nicholson S.E. and Kim J. (1997). The relationship of the El Niño/Southern Oscillation to African rainfall. Int. J. Climatol. 17, 117135. [ Links ]
19. Nicholson S.E. and Entekhabi D. (1987). Rainfall variability in equatorial and southern Africa. Part 1: Relationship with sea surface temperatures along the southwestern coast of Africa. J. Climate Appl. Meteor. 26, 561568. [ Links ]
20. Nicholson S.E. (1997). An analysis of the ENSO signal in the tropical Atlantic and western Indian Oceans. Int. J. Climatol. 17, 345375. [ Links ]
21. Mason S.J. (1990). Temporal variability of sea surface temperature around South Africa: a possible forcing mechanism for the 18-year rainfall oscillation? S. Afr. J. Sci. 86, 243252. [ Links ]
22. Mason S.J. (1995). Sea surface temperature-South African association 1910 1989. Int. J. Climatol. 15, 119135. [ Links ]
23. Mason S.J. (1998). Seasonal forecasting of south African rainfall using a non-linear discriminant analysis model. Int. J. Climatol. 18, 147164. [ Links ]
24. Reason C.J.C. (2001). Subtropical Indian Ocean SST dipole events and southern African rainfall. Geophys. Res. Lett. 28, 22252227. [ Links ]
25. Reason C.J.C. and Rouault M. (2002). ENSO-like decadal variability and South African rainfall. Geophys. Res. Lett. 29, doi:10.1029/2002GL014663. [ Links ]
26. Rouault M., Florenchie P., Fauchereau N. and Reason C.J.C. (2003). South east tropical Atlantic warm events and southern African rainfall. Geophys. Res. Lett. 30, doi:10.1029/2002GL014840. [ Links ]
27. Reason C.J.C., Hachigonta S. and Phaladi R.F. (2005). Interannual variability in rainy season characteristics over the Limpopo Province of South Africa. Int. J. Climatol. 25, 18351853. [ Links ]
28. Kruger A.C. (1999). The influence of decadal-scale variability of summer rainfall on the impact of El Niño and La Niña events in South Africa. Int. J. Climatol. 19, 5968. [ Links ]
29. Jury M. (1997). Statistical models for seasonal forecasts of southern Africa summer rainfall. CPC Experimental Long-lead Forecast Bulletin 6, no. 4, Decem- ber 1997. Online at: www.cpc.noaa.gov/products/predictions/experimental/ bulletin/Dec97/a30.html [ Links ]
30. Jury M.R., Mulenga H.M. and Mason S.J. (1999). Development of statistical long range models to predict summer climate variability over southern Africa. J. Climate 12, 18921899. [ Links ]
31. Landman W.A. and Mason S.J. (1999). Operational long-lead prediction of South African rainfall using canonical correlation analysis. Int. J. Climatol. 19, 10731090. [ Links ]
32. Engelbrecht F. (2005). Simulations of climate and climate change over southern and tropical Africa with the conformal-cubic atmospheric model. In Climate Change and Water Resources in Southern Africa: Studies on Scenarios, Impacts, Vulnerabilities and Adaptation, ed. R.E. Schulze, chap. 4, pp. 5774. Water Research Commission, Pretoria. [ Links ]
33. Hewitson B.C. and Crane R.G. (2006). Consensus between GCM climate change projections with empirical downscaling: precipitation downscaling over South Africa. Int. J. Climatol. 26, 13151337. [ Links ]
34. Labitzke K. and van Loon H. (1988). Associations between the 11-year solar cycle, the QBO, and the atmosphere. Part I: The troposphere and stratosphere in the northern hemisphere in winter. J. Atmos. Terres. Phys. 50, 197206. [ Links ]
35. Labitzke K. and van Loon H. (1989). Association between the 11-year solar cycle, the QBO, and the atmosphere. Part III: Aspects of the association. J. Climate 2, 554565. [ Links ]
36. Van Loon H. and Labitzke K. (1988). Association between the 11-year solar cycle, the QBO, and the atmosphere. Part II: Surface and 700 mb in the northern hemisphere in winter. J. Climate 1, 905920. [ Links ]
37. Mason S.J. and Tyson P.D. (1992). The modulation of sea surface temperature and rainfall associations over southern Africa with solar activity and the quasi-biennial oscillations. J. Geophys. Res. 97, 58475856. [ Links ]
38. Mason S.J. and Jury M.R. (1997). Climatic change and inter-annual variability over southern Africa: a reflection on underlying processes. Prog. Phys. Geogr. 21, 2350. [ Links ]
39. Mason S.J. and Tyson P.D. (2000). The occurrence and predictability of droughts over southern Africa. In Droughts, a Global Assessment, vol. 1, ed. D.A. Wilhite, pp. 113134. Routledge, New York. [ Links ]
40. Kruger A.C. and Shongwe S. (2004). Temperature trends in South Africa: 19602003. Int. J. Climatol. 24, 19291945. [ Links ]
41. Kruger A.C. (2006). Observed trends in daily precipitation indices in South Africa: 19102004. Int. J. Climatol. doi:10.1002/joc.1368. [ Links ]
42. Tadross M.A., Jack C. and Hewitson B.C. (2006). On RCM-based projections of change in southern African summer climate. Geophys. Res. Lett. 32, L23713, doi: 10.1029/2005GL024460. [ Links ]
43. New M., Hulme M. and Jones P.D. (1999). Representing twentieth century spacetime climate variability. Part 1: development of a 196190 mean monthly terrestrial climatology. J. Climate 12, 829856. [ Links ]
44. New M., Hulme M. and Jones P.D. (2000). Representing twentieth century space-time climate variability. Part 2: development of 190196 monthly grids of terrestrial surface climate. J. Climate 13, 22172238. [ Links ]
45. Mitchell T.D., Hulme M. and New M. (2002). Climate data for political areas. Area 34, 109112. [ Links ]
46. Jones P.D., New M., Parker D.E., Martin S. and Rigor I.G. (1999). Surface air temperature and its changes over the past 150 years. Rev. Geophys. 37, 173199. [ Links ]
47. Mason S.J. (2001). El Niño, climate change, and southern African climate. Environmetrics 12, 327345. [ Links ]
48. Burg J.P. (1967). Maximum Entropy Spectral Analysis. In Proceedings, 37th Meeting, Society of Exploration Geophysics, Oklahoma City. [ Links ]
49. Ulrych T.J. and Bishop T.N. (1975). Maximum Entropy Spectral Analysis and autoregressive decomposition. Rev. Geophys. 13, 183200. [ Links ]
50. Blackman R.B. and Tukey J.W. (1958). The Measurement of Power Spectra. Dover, New York. [ Links ]
51. Kane R.P. (1977). Power spectrum analysis of solar and geophysical parameters. J. Geomag. Geoelect. 29, 471495. [ Links ]
52. Kane R.P. (1979). Maximum Entropy Spectral Analysis of some artificial samples. J. Geophys. Res. 84, 965966. [ Links ]
53. Kane R.P. and Trivedi N.B. (1982). Comparison of maximum entropy spectral analysis (MESA) and least-square linear prediction (LSLP) methods for some artificial samples. Geophysics 47, 17311736. [ Links ]
54. Bevington P.R. (1969) Data Reduction and Error Analysis for the Physical Sciences, pp. 164176. McGraw-Hill, New York. [ Links ]
55. Trenberth K.E. (1993). The different flavors of El Niño. In Proceedings of the 18th Annual Climate Diagnostics Workshop, Boulder, CO, pp. 5053. [ Links ]
56. Quiroz R.S. (1983). Seasonal climate summary: the climate of the 'El Niño' winter of 198283 a season of extraordinary climate anomalies. Mon. Wea. Rev. 111, 16851706. [ Links ]
57. Kane R.P. (2004). Comparison of stratospheric zonal winds and El Niño/Southern Oscillation in recent decades. Int. J. Climatol. 24, 525532. [ Links ]
58. Kane R.P. (1999). Some characteristics and precipitation effects of the El Niño of 19971998. J. Atmos. Solar-Terr. Phys. 61, 13251346. [ Links ]
59. Landman W.A. and Mason S.J. (1999). Changes in the association between Indian Ocean sea-surface temperatures and South African summer rainfall. Int. J. Climatol. 19, 14771492. [ Links ]
60. Kripalani R.H. and Kulkarni A. (1997). Climate impact of El Niño/La Niña on Indian monsoon: a new perspective. Weather 52, 3946. [ Links ]
61. Kripalani R.H. and Kulkarni A. (1997). Rainfall variability over south-east Asia connections with Indian monsoon and ENSO extremes: new perspectives. Int. J. Climatol. 17, 11551168. [ Links ]
62. Tyson P.D. and Preston-Whyte R.A. (2000). The Weather and Climate of Southern Africa, 2nd edn. Oxford University Press, Cape Town. [ Links ]
63. Fauchereau N., Trzaska S., Roualt M. and Richard Y. (2003). Rainfall variability and changes in southern Africa during the 20th century in the global warming context. Nat. Hazards 29, 139154. [ Links ]
64. Servain J. and Seva M. (1987). On relationship between tropical Atlantic sea surface temperature, wind stress and regional precipitation indices: 19641984. Ocean Air Interaction 1, 183190. [ Links ]
65. Greischar L., Hastenrath S. and van Heerden J. (1997). Prediction of summer rainfall over eastern South Africa using empirical methods. CPC Experimental Long-lead Forecast Bulletin 6, no. 4, December 1997. Online at: http://www. cpc.noaa.gov/products/predictions/experimental/bulletin/Dec97/a39.html [ Links ]
66. Landman W. (1997). A Canonical Correlation Analysis Model to Predict South African 199798 Summer Rainfall and Temperature. CPC Experimental Long-lead Forecast Bulletin 6, no. 4, December 1997. Online at: http://www. cpc.noaa.gov/products/predictions/experimental/bulletin/Dec97/a36.html [ Links ]
67. Thiaw W. and Barnston A. (1997). CCA Forecast for Southern African Rainfall in JanFebMar 1998. CPC Experimental Long-lead Forecast Bulletin 6, no. 3, September 1997, pp. 3031. Online at: www.cpc.noaa.gov/products/predictions/ experimental/bulletin/Sep97/art11.html [ Links ]
68. Thiaw W. and Barnston A. (1998). CCA Forecast for Southern African Rainfall in JanFebMar 1999. COLA Experimental Long-lead Forecast Bulletin 7, no. 4, December 1998. Online at: http://grads.iges.org/ellfb/Dec98/thiaw.htm [ Links ]
69. Greischar L., Hastenrath S. and van Heerden J. (1998). Prediction of Summer Rainfall over Eastern South Africa Using Empirical Methods. COLA Experimental Long-lead Forecast Bulletin 7, no. 4, December 1998. Online at: http://grads. iges.org/ellfb/Dec98/greishar.htm [ Links ]
70. Landman W. and Mason S.J. (1998). Prediction of the January to March Rainfall of South Africa, Botswana and Namibia Using CCA. COLA Experimental Long-lead Forecast Bulletin 7, no. 4, December 1998. Online at: http://grads.iges.org/ellfb/ Dec98/landman.htm [ Links ]
71. Barnston A.G., Glantz M.H. and He Y.X. (1999). Predictive skill of statistical and dynamical climate models in SST forecasts during the 199798 El Niño episode and the 1998 La Niña onset. Bull. Amer. Meteorol. Soc. 80, 217243. [ Links ]
72. Bell G.D. and Halpert M.S. (1998). Climate assessment for 1997. Bull. Amer. Meteorol. Soc. 79, S1S50. [ Links ]
Received 15 February 2007. Accepted 5 March 2009.