On-line version ISSN 1816-7950
Water SA (Online) vol.36 no.3 Pretoria Apr. 2010
David LoveI, II,*; Stefan UhlenbrookIII,IV; Steve TwomlowII; Pieter van der ZaagIII,IV
IWaterNet, PO Box MP600, Mount Pleasant, Harare, Zimbabwe
IIICRISAT Bulawayo, Matopos Research Station, PO Box 776 Bulawayo, Zimbabwe
IIIUNESCO-IHE, Westvest 7, PO Box 3015, 2601 DA Delft, The Netherlands
IVDepartment of Water Resources, Delft University of Technology, PO Box 5048, 2600 GA Delft, The Netherlands
Changing regional and global trends in climate and discharge, such as global warming-related declines in annual rainfall in south-eastern Africa, are likely to have a strong influence on water resource availability, and to increase livelihood risk. It is therefore important to characterise such trends. Information can be obtained by examining and comparing the rainfall and runoff records at different locations within a basin. In this study, trends in various parameters of temperature (4 stations), rainfall (10 stations) and discharge (16 stations) from the northern part of the Limpopo Basin, Zimbabwe, were statistically analysed, using the Spearman rank test, the Mann-Kendall test and the Pettitt test. It was determined that rainfall and discharge in the study area have undergone a notable decline since 1980, both in terms of total annual water resources (declines in annual rainfall, annual unit runoff) and in terms of the temporal availability of water (declines in number of rainy days, increases in dry spells, increases in days without flow). Annual rainfall is negatively correlated to an index of the El Niño Southern Oscillation phenomenon. The main areas of rising risk are an increasing number of dry spells, which is likely to decrease crop yields, and an increasing probability of annual discharge below the long-term average, which could limit blue-water availability. As rainfall continues to decline, it is likely that a multiplier effect will be felt on discharge. Increasing food shortages are a likely consequence of the impact of this declining water resource availability on rain-fed and irrigated agriculture. Declining water resource availability will also further stress urban water supplies, notably those of Zimbabwe's second-largest city of Bulawayo, which depends to a large extent from these water resources and already experiences chronic water shortages.
Keywords: climate variability, climate change, discharge analysis, Pettitt test, rainfall analysis, water resources, Limpopo Basin, Zimbabwe, Southern Africa
Rainfall in south-eastern Africa is temporally and spatially intermittent (Unganai and Mason, 2002). Annual rainfall for some sites can vary by up to 1 000 mm·a-1 from year to year (Twomlow and Bruneau, 2000).A drought year may record less than 250 mm·a-1, such as the 2004-2005 season in the Limpopo Basin (Love et al., 2006). Rainfall variability in south-eastern Africa is strongly influenced by the coupled ocean-atmosphere El Niño Southern Oscillation phenomenon (ENSO) (Trenberth et al., 2007). Positive ENSO anomalies generally result in reduced rainfall in the region and are becoming more common (Makarau and Jury, 1998; Alemaw and Chaoka, 2006). Furthermore, there has been a general decline in rainfall in Southern Africa since 1961 (New et al., 2006), with the period 1986-1995 being the driest decade of the 20th century (Trenberth et al., 2007). General circulation models developed with the Intergovernmental Panel on Climate Change Special Report on Emission Scenarios (IPCC SRES) scenarios suggest that annual rainfall in south-eastern Africa will decline further under the impact of global warming (Christensen et al., 2007). This is expected to be between 10% and 20% below the 19001970 averages by 2050 (Milly et al., 2008) or by up to 10% below the 1980-1999 averages by 2099 (Christensen et al., 2007).
Declines in rainfall, as discussed above, may translate into more than proportional declines in discharge due to non-linear processes, including, for example, interception thresholds. By the years 2041-2060 in south-eastern Africa, runoff is expected to decline by between 10% and 40% compared to 1900-1970 averages (Milly et al., 2005: Fig. 4a).
These changing regional and global trends in climate and discharge are likely to have a strong influence on water resource availability, and increase livelihood risk. It has been shown that household food security in southern Africa is highly vulnerable to climate stress (Archer et al., 2007). Already in much of southern Africa there is a precarious balance between available water resources and water demand as a result of generally low conversion of rainfall to runoff and potential evaporation exceeding rainfall (e.g. Mazvimavi, 2004). Frequently, the water yield from the developed surface water resource falls short of the demand, deficits being more evident during the frequent droughts (e.g. Nyabeze, 2004). Furthermore, some catchments, especially within the Limpopo Basin, are already over-committed (Kabel, 1984; Basson and Rossouw, 2003), leading to water stress: a high ratio of water withdrawal or water use to discharge (Vörösmarty et al., 2000). Changes of this nature constitute a major challenge to water resources management (Milly et al., 2008).
Information on water resource availability, past and present trends and future predictions, can be obtained by examining and comparing the rainfall and runoff records at different locations within a basin. Multi-decadal temporal trends of discharge, rainfall and temperature can be analysed to assess risk to strategic water resource systems (e.g. Hundecha and Bárdossy, 2005; Chen et al., 2007) or to attribute causes of change (e.g. Anderson et al., 2001).
In this study, temperature, rainfall and discharge trends from a series of locations in the northern part of the Limpopo Basin, Zimbabwe, are analysed, in order to determine the:
- Stationarity of the time series and any observed changes on water resource availability
- Causes of any changes observed
- Relationships between rainfall and discharge in the study area and causes thereof
- Implications of any changes
The study area is the drainage basins of the Shashe (18 991 km2) and Mzingwane Rivers (15 695 km2), in the northern Limpopo Basin, and selected stations are within Zimbabwe (Fig. 1). The 2 sub-basins contribute 12.3% and 9.3% of the mean annual runoff of the transboundary Limpopo basin, respectively. This makes them the third- and fourth-largest contributors to the Limpopo basin, after the Olifants (32.7%, 68 450 km2 area) and Luvuhu (13.1%, 4 826 km2 area) sub-basins (Görgens and Boroto, 1997).
Occurrence of rainfall in Zimbabwe is mainly connected to the Inter Tropical Convergence Zone (ITCZ). The movement of the ITCZ from the equator marks the start of the rainy season in the southern hemisphere. In a normal year, it fluctuates half-way between Tanzania and Zimbabwe but never moves beyond the Limpopo River in the south. The ITCZ moves with the sun, southwards at the start of summer (October/November) and northwards in late summer (March/April) (Twomlow et al., 2006). Because of this, the wet season in the southern parts of Zimbabwe (including the northern part of the Limpopo Basin, also known as the Mzingwane Catchment) starts later and ends sooner than in the northern areas. Furthermore, the northern winds in the convergence are moister than the southern winds, leading to less frequent rainfall in the southern areas than the north, for a given air moisture level. Rainfall in southern Zimbabwe thus occurs over a limited period of time, and often a large portion of the annual rainfall can fall in a small number of events (Twomlow and Bruneau, 2000; Unganai and Mason, 2002; De Groen and Savenije, 2006). These short, intense rain spells, if they occur at critical stages during crop growth, can exercise a strong control on crop yields (Twomlow et al., 2006). Within the study area, annual rainfall generally decreases from the North to the South: from around 630 mm·a-1 at Esigodini, to 560 mm·a-1 at Filabusi, to 360 mm·a-1 at Beitbridge (see Table 2).
A statistical analysis was carried out of daily temperature, rainfall and runoff data from several stations (see Fig. 1). Stations with a minimum of 30 years continuous data were selected. Data were organised into the Southern African hydrological year of October to September.
All selected discharge stations were upstream of major dams and with limited upstream water users (see table 3 for water use upstream of the discharge stations and Fig. 1 for the locations of the discharge stations in relation to major dams) and did not require naturalisation.
The time series were visually inspected, along with supporting materials such as the station files. The following exclusions were made for each station, in order to remove unreliable data:
- Where data were missing for 2 months or more, the year was excluded
- Where data were missing for 2 weeks or more during the months of November to April (rainy season), the year was excluded
- Where a note had been made in the station file that readings were unreliable (e.g. due to siltation, security), the year was excluded
- Where the runoff coefficient for a given year was more than 1, the rainfall data were checked for correspondence with the next nearest climate station. If this was unsatisfactory, then the discharge data for that year were excluded.
- Where a daily average flow of more than zero but less than 0.02 m3/s was recorded for 1 month or more, the year's data for 'number of days of no flow' was excluded as being below the detection level of the gauge. The data selected for analysis, following these exclusions, are shown in Tables 2 and 3.
The Multivariate ENSO Index (MEI) used was developed by Wolter and Timlin (1998) from sea-level pressure, zonal and meridional components of the surface wind, sea surface temperature, surface air temperature and total cloudiness fraction of the sky. An average MEI for each hydrological year in the period 1951-2005 was obtained from published bimonthly MEI values (Wolter, 2007).
A series of non-parametric tests were carried out on several parameters that were derived from the temperature, rainfall, discharge and MEI dataset. For datasets which are not drawn from a population with specific statistical conditions, such as normal distribution, non-parametric tests are appropriate (McCuen, 2003). This is the case for the study data, which is not unusual for hydroclimatic data (Tilahun, 2006). To determine whether there was an overall trend in each time series studied, 2 trend analysis tests were used. The 1st trend analysis test, the Spearman rank test, is a simple measure of correlation between 2 data series and is a special case of the Pearson product-moment coefficient (Myers and Well, 2003). It has often been used for temporal analysis of climatic variables (e.g. Yu and Neil, 1993; González-Hidalgo et al., 2001). The Mann-Kendall test was also applied. This is a more elaborate test for identifying shift trends (e.g. Shao et al., 2009), and compares the relative magnitudes of the data, as opposed to their actual values (Gilbert, 1987). It has been widely used for assessing the significance of temporal trends in hydrological and climatic data (Hirsch et al., 1982), including in semi-arid regions (Zhang et al., 2008).
The time series from each parameter was also subjected to the Pettitt Test (Pettitt, 1979), which is often used to detect abrupt changes in hydrological series (e.g. Shao et al., 2009). The test determines the timing of a change in the distribution of a time series, known as a 'change point' (Zhang et al., 2008). The change point divides the series into 2 sub-series. The significance of the change point is then assessed by F- and t-tests on the change in the mean and the variance. The Pettitt test was used to identify change points in the time series, at a probability threshold of p = 0.8, followed by F- and t-tests at 2.5% significance level. This procedure has been used for identifying change points in hydroclimatic data in both humid and semi-arid environments (Tu et al., 2005; Ashagrie et al., 2006; Zhang et al., 2008). Tests were carried out using SPELL-stat v.188.8.131.52B (Guzman and Chu, 2004). To determine whether or not change points that were shown to be significant were also substantively important, the effect size was determined. The Cohen's d statistic was computed, with values of 0.8 and above considered to have a 'large' effect size (Cohen, 1988).
Using the Southern African hydrological year of October to September, the following time series were prepared from temperature records from Matopos Research Station (for which the longest time series in the region is available):
- Annual average of daily maximum temperature in °>C
- Annual average of daily minimum temperature in °>C
- Annual maximum (extreme) of daily maximum temperature in °>C
- Annual minimum (extreme) of daily minimum temperature in °>C
The following time series were prepared for rainfall for all climate stations:
- Annual rainfall in mm/a
- Number of wet days in d/a with rainfall >0.5 mm/d
- Number of heavy rain days per year with rainfall >10 mm/d
- Number of heavy rain days per year with rainfall >20 mm/d
- Number of heavy rain days per year with rainfall >30 mm/d
- Length of longest dry spell for the months November-March (days)
- Number of 5 to 7 d dry spells for the months November to March (-)
- Number of 8 to 14 d dry spells for the months November to March (-)
- Number of more than 14 d dry spells for the months November to March (-)
- Number of more than 20 d dry spells for the months November to March (-)
- Total number of dry spells of 5 d or more for the months November to March (-)
Annual rainfall anomalies were prepared after the method of Hulme et al. (2001) for each station's annual rainfall series. The Thiessen polygon method was applied to estimate catchment rainfall; given the small number of rainfall stations per catchment, a more complex interpolation method was not justified. Catchment rainfall and discharge were compared for each discharge station at an annual time step and for selected stations at monthly and daily time steps.
The following time series were prepared for all discharge stations:
- Annual unit runoff (mm/a)
- Annual maximum of monthly flows (m3/month)
- Days per annum no flow
- Maximum average daily flow recorded per year (m3/s)
- Maximum flood recorded (instantaneous flows) per year (m3/s)
- Runoff coefficient (-)
A simple risk analysis was carried out for each change point where a medium or large effect size had been determined:
P is the probability that load exceeds risk
l is the load: the behaviour of the system under an external stress
r is resistance: the capacity of the system to overcome the load (Ganoulis, 2004)
For each parameter, a value was taken for resistance based on expert knowledge and related studies; each parameter represented a livelihood or an ecological risk. Thus for the selected parameters, the risk is the probability of occurrence of the below conditions:
- Temperature more than 3.1°C above the long-term mean for Bulawayo, shown to result in reduction of maize yield by 16% (Dimes et al., 2009)
- Rainfall under 450 mm/a, the minimum rainfall recommended for rain-fed maize production (Shumba and Maposa, 1996; Hoffmann et al., 2002)
- Less than 60 and less than 30 wet days per year: representing 1 wet day per 3 days and per 6 days, respectively, in the 6-month rainy season
- Less than 20 days per year with rainfall of over 10 mm/d, often considered the baseline minimum of heavy rainfall events (Lebel et al., 2000; Douville et al., 2001). Heavy rainfall events are very important for runoff generation in an environment with high interception: a threshold of 5 mm/d has been suggested for the study area (Love at al., 2010).
- More than 10 days per year with rainfall of over 20 mm/d, which could result in waterlogging or flooding
- At least 1 dry spell per year lasting over 14 d: dry spells of this length are considered critical in reducing the yield of maize (Lal, 1997; Magombeyi and Taigbenu, 2008).
- At least 1 dry spell per year lasting over 20 d: more extreme than the preceding risk measure
- Annual unit runoff less than 22 mm/a: this is the mean annual unit runoff of the Mzingwane Catchment (Görgens and Boroto, 1997; Chibi et al., 2006)
- Days of no flow more than 219 d/a: this shows the river to be flowing less than 40% of the time.
- Maximum average daily flow recorded per year under 5·m3/s and maximum instantaneous flood recorded per year under 20 m3/s: the occurrence of sufficiently large floods annually is important for hydro-ecological processes such as habitat maintenance (King and Louw, 1998). The values selected are indicative as the ecological significance of a specific flow volume or return period is specific to each river, and a detailed analysis is beyond the scope of this paper.
Figure 2 summarises the minimum and maximum daily temperature trends observed for 4 locations in the study area, with peaks corresponding with the occurrence of El Niño events (high positive Multivariate ENSO Index - MEI). Similar trends were observed at the 4 locations, but our focus will be on the Matopos Research Station location, which had the longest data set, some 20 years longer. The annual average of maximum daily temperatures at Matopos Research Station (Table 4) showed a significantly increasing trend (Spearman rank test). The Multivariate ENSO Index also shows a significantly increasing trend (Spearman rank test).
Temperature parameters recorded at Matopos Research Station show change points in 1975 and the mid to late 1980s, with the average daily maximum temperature rising from 1975 and the annual average of maximum daily temperatures rising from 1986 (Table 4). A drop in minimum temperatures is recorded from around 1990 onwards.
The high inter-annual variability of rainfall in the northern Limpopo Basin is similar to that reported for much of Zimbabwe (Twomlow et al., 2006). The influence of the ENSO phenomenon is shown by the correspondence of Multivariate El Niño/Southern Oscillation Index (MEI) peaks with unusually low rainfall for seasons such as 1982/83, 1986/87, 1991/92 and 1997/98 (Fig. 3). For almost all stations there is a strong negative correlation between total annual rainfall and MEI: Mzingwane Dam (r=0.43), Bulawayo (r=0.46), Matopos Research Station (r=0.50), Filabusi (r=0.56), Mbalabala (r=0.56), Kezi (r=0.63), Beitbridge (r=0.70) and West Nicholson (r=0.71). Cyclicity of 17 to 20 years is apparent in the 5-year moving averages and to a lesser extent in the annual rainfall anomalies (Fig. 4 ).
Spearman rank and Mann-Kendall tests determined that only a limited number of parameters showed a significant drying trend for the whole time series (Table 5). These include the total annual rainfall and days of heavier rainfall at Matopos National Park, as well as the number of longer dry spells at Filabusi and Matopos Research Station. The Pettitt test identified numerous change points, with a significant difference in means confirmed by t-test. The majority showed a medium size effect (0.3 < d < 0.8) and several showed a large size effect (d ≥ 8) (Table 5). Data for Bulawayo and Esigodini did not show significant change points and are not displayed. The majority of these showed a change point to a drier regime from around 1980, with a smaller number of change points that showed a change to a drier regime in the 1960s, as well as around 1980 (see Fig. 6 for examples of the differences identified in Table 5). F-tests showed that most change points also showed a significant decrease in the variance of the rainfall parameter in question. All stations except Bulawayo and Matopos Research Station showed change to fewer wet days. Three stations showed change to an increased number of dry spells, a decreased annual rainfall or both. Only 1 station showed a change to fewer heavy storms. None of the stations showed any change to any wetter regime.
The Pettitt test identified numerous change points, with a significant difference in means confirmed by t-test for discharge data obtained from 16 locations in the basin (Table 6). The majority of change points had a large size effect (d > 0.8)Most stations on tributaries of the Mzingwane, Thuli and Shashe Rivers show a change point to a drier regime from around 1980 (see Fig. 6 for examples of the differences identified in Table 6). Most change points also showed a significant decrease in the variability of the parameter in question, confirmed by F-test although only about half of the stations show a decline in coefficient of variation. Almost all stations showed a change to a smaller number of days of flow and lower annual maximums of average monthly flows.
Many stations showed a consistently drying trend across the full time series for some parameters, including absolute declines in annual unit runoff at 6 stations, confirmed by Spearman rank and Mann-Kendall tests (Table 6). Data for B60 and B87 did not show significant change points and are not displayed. However, data from B87 showed continuous drying trends across the full time series, based on Spearman tests on the annual unit runoff, annual maximum of monthly flows, maximum daily flow per year and maximum flood recorded. Change points detected in B77 occur close to gaps in the time series and these results are also not displayed.
Only a few discharge stations showed a clear relationship between annual unit runoff and annual catchment rainfall (Fig. 7 and Table 7); these are predominantly tributaries of the Shashe and Thuli Rivers in the west: B77, B78, B80 and B83 (Fig. 1). All of these stations showed change to a drier regime in or around 1980 with at least 1 parameter showing an absolute declining trend across the time series (Table 6), although which parameter showed such a decline is not constant between the stations.
Change over time in the relationship between rainfall and runoff is shown by temporal change in the annual runoff coefficients (Table 7). Three stations show a significant decline in runoff coefficient, across the whole time series, 4 stations show a change to lower runoff coefficients in the early 1980s, and 1 station shows such a change in the mid 1980s. Examination of rainfall-runoff relationships at sub-annual time steps did not give acceptable levels of correlation, probably due to one of the following factors: differences in rainfall and discharge measurement days, delays between rainfall and discharge peaks at a short time step, or the high spatial variability of rainfall in Zimbabwe.
The most important climatic risk observed is the increased probability of dry spells, with the probability of dry spells over 14 d and over 20 d rising to above 0.7 in many stations (Table 8). The probability of the number of rainy days per year falling below 60 and below 30 also increased substantially after each change point (mainly in 1979/80). There is no substantive risk posed by temperature and little increase in risk posed by annual rainfall levels. The probability of heavy storms declined at the 2 stations for which statistically significant change was observed.
There is a high risk of low annual runoff occurring. The probability of the overall Mzingwane Catchment average unit runoff exceeding a study catchment's annual unit runoff doubled at most stations after the change points (Table 9). At all stations except B80, the probability of the river flowing less than 40% of time has increased substantially since the change points. In most cases, the probability of lower maximum flows and floods has increased.
Stationarity of the climatic and hydrological time series
The general trends in the time series, as identified by Spearman rank and Mann-Kendall tests, observed in the northern Limpopo Basin and the significant change points indicate:
- The relationship between rainfall and temperature and ENSO (Figs. 2 and 3); see also Alemaw and Chaoka (2006)
- An increase in maximum daily temperatures (Table 3); see also New et al. (2006)
- A decline in rainfall (Table 5 and Fig. 5) with high inter-decadal variability (Fig. 4); see also Trenberth et al. (2007)
- A 3- to 5-year cyclicity in rainfall (Fig. 4), similar to and Jury and Mpeta (2005) and possibly related to El Niño. This cycle could be contained within a 17- to 20-year cyclicity in rainfall (Fig. 4) which could be similar to Tyson (1986) and Mason et al. (1997). Thus, there is possibly a shorter cycle within a longer cycle.
The vast majority of climate stations showed a change towards a drier climate from the early 1980s (Table 5) the notable exception being Bulawayo, which is at a higher elevation and actually outside the Limpopo Basin. Decreased land rainfall from this period onwards was reported at global level by Trenberth et al. (2007). Decline in rainfall reflecting a new rainfall regime, as opposed to a gradual decrease, has also been reported from western Australia (Hope et al., 2006).
A corresponding change to a drier regime was also observed for most discharge stations (Table 6), notably from the early 1980s in the stations upstream of major dams (Mzingwane, Shashe and Thuli tributaries; see Fig. 1 for locations). Such change in runoff regime has not been previously reported in south-eastern Africa. In 6 of the 16 stations (there is no consistent geographical pattern), there was a significant decline in annual unit runoff for the full time series. The change to a regime of relatively lower or less frequent rainfall and discharge demonstrates that water resource availability has shown a decline in the northern Limpopo Basin in the last half-century, even if trends in parameters such as annual rainfall are obscured in many places by the observed cyclicity.
Causes of change in rainfall and discharge
The decline in rainfall and other drying trends in rainfall parameters (Table 5) could be related to:
- The increasing frequency and severity of El Niño events (Alemaw and Chaoka, 2006).
- Changes due to global warming: the more pronounced increase in number of dry spells with more limited changes in total rainfall is similar to those predicted from global circulation models by Christensen et al. (2007) for Southern Africa. The median precipitation response predicted is zero for December to May rainfall with a 13% decline in September to November rainfall (in southern Zimbabwe there is often very little rainfall during September and October). The number of dry spells are predicted to increase in most subtropical areas.
- Possible effects of land use change on moisture recycling. Large areas of Zimbabwe, including parts of the study area, have experienced significant population growth, agricultural intensification and loss of forest cover in the last 50 years (e.g. Mapedza et al., 2003; Zinyama and Whitlow, 2004). However, changing sea surface temperatures (and thus ENSO) are thought to be more important than changing land use patterns in controlling warm season rainfall variability and trends in Southern Africa (Christensen et al., 2007). This needs further research.
- Observed cyclicity in rainfall: Fig. 5 shows that the early 1980s are at the lower point of the 5-year moving average of annual rainfall. However, the cyclicity exhibited would suggest another change point (to a wetter regime) would be expected by the late 1990s. Such a change could not (yet) be observed, but may become apparent once more recent data become available.
The change to a drier regime in many discharge stations (Table 6) could be related to the following factors:
- Declining rainfall, as observed in Fig. 5, resulting in declining runoff generation. This explanation is favoured by the observed co-incidence of change points to decreasing rainfall and decreasing discharge in the early 1980s.
- Increased upstream withdrawals due to small dam construction. The Zimbabwe Government carried out an extensive dam-building programme; an estimated 15 000 small dams were built in Zimbabwe in the early to mid 1980s (e.g. Jewsbury and Imevbore, 1988; Chimomwa and Nugent, 1993).
- Post-1980 increases in agricultural water use, which will have occurred due to increased production since independence in 1980, at both smallholder and commercial scale (Thirtle et al., 1993), as well as due to changing population distribution, especially in areas which had experienced population displacement during the civil war (Zinyama and Whitlow, 2004).
- Possible effects of land-use change on hydrology: increasing population density and conversion of forest or rangeland to agricultural land, which has occurred in Zimbabwe over the past 50 years, see above, can cause a decrease in runoff (e.g. Calder et al., 1995) although it can also cause an increase in runoff (e.g. Bari and Smettem, 2004).
It should be noted that correlation of runoff with catchment rainfall will have been weakened by poor rainfall station distribution in the study area. Only a few discharge stations showed a good correlation (B77, B78, B80 and B83; see Table 7). These stations are all from the west of the study area (see Fig. 1), where small dams are fewer, farm sizes are larger and crop production is minimal (Surveyor General of Zimbabwe and Forestry Commission, 1996); thus abstractions by farmers are likely to be lower in these areas. Seven discharge stations spread across the study area (B15, B40, B61, B74, B77, B83 and B87, Table 7) have declining runoff coefficients, showing a possible change in the runoff generation regime. Field studies, using experimental catchments with rainfall measured from several stations within the catchment, rather than interpolated from the national rainfall stations used in this study, might assist in further understanding of this problem.
The principal livelihood risks posed by the changing hydroclimatic patterns observed are: risk of rainfed crop failure, due to an increasing probability of dry spells over 14 days and 20 days in length (Table 8), and risk to water supply (for agriculture and domestic use) from an increasing probability of annual unit runoff dropping below the Mzingwane Catchment average (Table 9). The changes in temperature and total annual rainfall do not appear to impose risk at this time. Changing flow regimes, with a higher probability of smaller-sized floods, may have ecological implications, but this cannot be assessed without a full flood analysis.
There is a complex variety of factors that may influence the temporal distribution of rainfall and discharge in southern Africa, some of which result in cyclicity. This complexity may be part of the reason that hydroclimatic studies using different statistical techniques make different findings. For example, Mazvimavi (2008) did not detect significant change in extreme high or low rainfall in 40 stations across Zimbabwe, using the Mann-Kendal test and quantile regression analysis. However, this study has shown that, in terms of the statistical tests selected in this study, total seasonal rainfall and annual discharges in the northern Limpopo Basin, Zimbabwe, have shown notable declines since 1980 at many measuring stations. The evidence from the time series analyses shows a declining rainfall trend that appears to be related to the increased incidence and severity of El Niño events and to changes associated with global warming. Discharge has been reduced by the declining rainfall input, but also by declining runoff generation in headwater catchments, likely due to impoundments in small dams and increased abstractions for irrigated crop production, especially in the north and east of the study area. If rainfall continues to decline, it is likely that a multiplier effect will be felt on discharge, due to non-linear processes resulting in more than proportional declines in discharge (De Wit and Stankiewicz, 2006), as well as increased agricultural water abstractions by farmers in the headwater catchments to compensate with irrigation for the reduced rainfall.
It has been seen that since 1980 water resource availability in the study area has declined, both in terms of total annual water available for storage (i.e. declines in annual rainfall, annual unit runoff) and in terms of the frequency of water availability (i.e. declines in number of rainy days, increases in dry spells, increases in days without flow). Although construction of dams can mitigate inter-annual variability and the declining frequency of water availability, the decline in total annual water available for storage will decrease the yield of dams. These trends make it unlikely that the construction on the Mzingwane River of the proposed Glassblock Dam (Chibi et al., 2006), located between Mbalabala and West Nicholson (see Fig. 1), would be beneficial, and will further stress urban water supplies in the study area, notably those of Zimbabwe's second-largest city of Bulawayo, which already experiences chronic water shortages (Gumbo, 2004). Further downstream, the proposed expansion of irrigation from Zhovhe Dam to provide for smallholder cereal production as well as cash crops may also be compromised, along with the proposed transboundary water supply to Musina in South Africa. Taken together, these trends suggest that water supply development for south-western Zimbabwe (including any water exports from Zimbabwe) should focus on alternatives to the Mzingwane Catchment and the Limpopo Basin. It would be strategic to link such alternatives to the long-awaited Matabeleland-Zambezi water carrier.
Drought and subsequent crop failure is becoming an all too common feature of south-eastern Africa (Richardson, 2007). This study has shown a rising risk of crop failure, due to an increasing probability of dry spells longer than 14 days occurring. Declining water resource availability in the region has had, and will undoubtedly continue to have, dramatic impacts on staple food production and food security (Du Toit and Prinsloo, 2001; Stige et al., 2006). The declining frequency of water availability will make the situation worse, and it has been shown that, at least for some parts of the study area, there is an increasing risk of annual runoff below the long-term catchment average. Unfortunately, many strategies to address hunger in southern Africa are water intensive, and some have been shown to result in water use conflicts (Love et al., 2006). Such initiatives should be complemented by interventions that maximise the use of the existing scarce water resources. These include soil water conservation techniques (e.g. Mupangwa et al., 2006), rainwater harvesting (e.g. Mwenge Kahinda et al., 2007), improving the management of existing irrigation schemes (e.g. Samakande et al., 2004) and the use of alternative crops such as switching from maize to sorghum when a drier season is to be expected based on the condition of the ENSO phenomenon (Stige et al., 2006).
This paper is an output of the CGIAR Challenge Programme on Water and Food Project 17 'Integrated Water Resource Management for Improved Rural Livelihoods: Managing risk, mitigating drought and improving water productivity in the water scarce Limpopo Basin' led by WaterNet, with additional funding provided by the International Foundation for Science (Grant W4029-1). The opinions and results presented in this paper are those of the authors and do not necessarily represent the opinions of the donors or participating institutions. Discharge data were provided by the Zimbabwe National Water Authority and climate data by the Department of Meteorological Services, Zimbabwe Ministry of Environment and Tourism. Dr Marieke de Groen and Dr Washington Nyabeze are thanked for useful discussions.
ALEMAW BF and CHAOKA TR (2006) The 1950-1998 warm ENSO events and regional implication to river flow variability in Southern Africa. Water SA 32 459-463. [ Links ]
ANDERSON ML, KAVVAS ML and MIERZWA MD (2001) Probabilistic/ensemble forecasting: a case study using hydrologic response distributions associated with El Niño Southern Oscillation (ENSO). J. Hydrol. 249 134-147. [ Links ]
ARCHER E, MUKHALA E, WALKER S, DILLEY M and MASAMVU K (2007) Sustaining agricultural production and food security in Southern Africa: an improved role for climate prediction? Clim. Change 83 287-300. [ Links ]
ASHAGRIE AG, DE LAAT PJM, DE WIT MJM, TU M and UHLENBROOK S (2006) Detecting the influence of land use changes on discharges and floods in the Meuse River Basin the predictive power of a ninety-year rainfall-runoff relation? Hydrol. Earth Syst. Sci.10 691-701. [ Links ]
BARI M and SMETTEM KRJ (2004) Modelling monthly runoff generation processes following land use changes: groundwater-surface runoff interactions. Hydrol. Earth Syst. Sci. 8 903-922. [ Links ]
BASSON MS and ROSSOUW JD (2003) Olifants Water Management Area: Overview of Water Resources Availability and Utilization. Report P WMA 04/000/00/0203, Department of Water Affairs and Forestry, Pretoria, South Africa. 55 pp. [ Links ]
CALDER IR, HALL RL, BASTABLE HG, GUNSTON HM, SHELA O, CHIRWA A and KAFUNDU R (1995) The impact of land use change on water resources in sub-Saharan Africa: a modelling study of Lake Malawi. J. Hydrol. 170 123-135. [ Links ]
CHEN H, GUO S, XU C and SINGH VP (2007) Historical temporal trends of hydro-climatic variables and runoff response to climate variability and their relevance in water resource management in the Hanjiang basin. J. Hydrol. 344 171-184. [ Links ]
CHIBI T, SAKAHUNI RC and MASVOPO TH (2006) Mzingwane Catchment Outline Plan. Zimbabwe National Water Authority, Bulawayo, Zimbabwe. 60 pp. [ Links ]
CHIMOMWA M and NUGENT C (1993) A fisheries GIS for Zimbabwe: an initial analysis of the numbers, distribution and size of Zimbabwe's small dams. Rep. FAO/UNDP ZIM/88/021. Fisheries and Aquaculture Department, Food and Agriculture Organisation of the United Nations. URL: http://www.fao.org/docrep/field/003/AB969E/AB969E00.htm (Accessed 22 December 2008). [ Links ]
CHRISTENSEN JH, HEWITSON B, BUSUIOC A, CHEN A, GAO X, HELD I, JONES R, KOLLI RK, KWON W-T, LAPRISE R, MAGAÑA RUEDA V, MEARNS L, MENÉNDEZ CG, RÄISÄNEN J, RINKE A, SARR A and WHETTON P (2007) Regional Climate Projections. In: Solomon SD, Qin M, Manning Z, Chen M, Marquis KB, Averyt M, Tigno R and Miller HL (eds.) Climate Change 2007: The Physical Science Basis. Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge University Press, Cambridge. 847-940. [ Links ]
COHEN J (1988) Statistical Power Analysis for the Behavioural Sciences Analysis (2nd edn.). Lawrence Erlbaum, Philadelphia. 567 pp. [ Links ]
DE GROEN MM and SAVENIJE HHG (2006) A monthly interception equation based on the statistical characteristics of daily rainfall. Water Resour. Res. 42 W12417, DOI:10.1029/2006WR005013. [ Links ]
DE WIT M and STANKIEWICZ J (2006) Changes in surface water supply across Africa with predicted climate change. Sci. 311 1917-1921. [ Links ]
DOUVILLE H, CHAUVIN F and BROQUA H (2001) Influence of Soil Moisture on the Asian and African Monsoons. Part I: Mean Monsoon and Daily Precipitation. J. Clim. 14 2381-2403. [ Links ]
DIMES J, COOPER P and RAO KPC (2009) Climate change impact on crop productivity in the semi-arid tropics of Zimbabwe in the 21st century. In: Humphreys L (ed.) Proceedings of the Workshop on Increasing the Productivity and Sustainability of Rainfed Cropping Systems of Poor, Smallholder Farmers, Tamale, Ghana, 22-25 September 2008. The CGIAR Challenge Program on Water and Food, Colombo. [ Links ]
DU TOIT AS and PRINSLOO MA (2001) El Niño-southern oscillation effects on maize production in South Africa: A preliminary study. In: Impact of El Nino and Climate Variability on Agriculture. ASA Spec. Pub. 63 77-86. [ Links ]
GANOULIS J (2004) Integrated Risk Analysis for Sustainable Water Resources Management. In: Linkov I and Ramadan ABD (eds.) Comparative Risk Assessment and Environmental Decision Making. NATO Science Series, Kluwer Academic Publishers, Dordrecht. 275-286. [ Links ]
GILBERT RO (1987) Statistical Methods for Environmental Pollution Monitoring. Wiley, New York, 320 pp. [ Links ]
GONZÁLEZ-HIDALGO JC, DE LUIS M, RAVENTÓS J and SÁNCHEZ JR (2001) Spatial distribution of seasonal rainfall trends in a western Mediterranean area. Int. J. Clim. 21 843-860. [ Links ]
GÖRGENS AHM and BOROTO RA (1997) Limpopo River: flow balance anomalies surprises and implications for integrated water resources management. Proc. 8th S. Afr. Natl. Hydrol. Symp. Pretoria, South Africa. URL: http://www.wrc.org.za/wrcpublications/wrcsanciahs/boroto.htm (Accessed 15 August 2005). [ Links ]
GUMBO B (2004) The status of water demand management in selected cities of southern Africa. Phys. Chem. Earth 29 1225-1231. [ Links ]
GUZMAN JA and CHU ML (2004) SPELL-stat v 15110 B. Grupo en Prediccion y Modelamiento Hidroclimatico Universidad Industrial de Santander, Colombia. Hydrological Modelling and Prediction Group, Industrial University of Santander, Colombia. [ Links ]
HEWITSON BC and CRANE RG (2006) Consensus between GCM climate change projections with empirical downscaling: precipitation downscaling over South Africa. Int. J. Clim. 26 1315-1337. [I cannot find this citation in the text?] [ Links ]
HIRSCH RM, SLACK JR and SMITH RA (1982) Techniques of trend analysis for monthly water quality data, Water Resour. Res. 18 107-121. [ Links ]
HOFFMANN S, CSITÁRI G and HEGEDÜS L (2002) The humus content and soil biological properties as a function of organic and mineral fertilization. Arch. Agron. Soil Sci. 48 141-146. [ Links ]
HOPE PK, DROSDOWSKY W and NICHOLLS N (2006) Shifts in the synoptic systems influencing southwest Western Australia. Clim. Dyn. 26 751-764. [ Links ]
HULME M, DOHERTY RM, NGARA T, NEW MG and LISTER D (2001) African climate change: 1900-2100. Clim. Res. 17 145-168. [ Links ]
HUNDECHA Y and BÁRDOSSY A (2005) Trends in daily precipitation and temperature extremes across western Germany in the second half of the 20th century. Int. J. Clim. 25 1189-1202. [ Links ]
JEWSBURY JM and IMEVBORE AMA (1988) Small dam health studies. Parasitol. Today 4 57-59. [ Links ]
JURY MR and MPETA EJ (2005) The annual cycle of African climate and its variability. Water SA 31 (1) 1-8. [ Links ]
KABEL TC (1984) An Assessment of Surface Water Resources of Zimbabwe and Guidelines for Development Planning. Government of Zimbabwe, Harare, Zimbabwe. 16 pp. [ Links ]
KING J and LOUW D (1998) Instream flow assessments for regulated rivers in South Africa using the Building Block Methodology. Aq. Eco. Health Man. 1 109-124. [ Links ]
LAL R (1997) Long-term tillage and maize monoculture effects on a tropical Alfisol in western Nigeria. I. Crop yield and soil physical properties. Soil Till. Res. 42 145-160. [ Links ]
LEBEL T, DELCALUX F, LE BARBÉ L and POLCHER J (2000) From GCM scales to hydrological scales: rainfall variability in West Africa. Stoch. Env. Res. Risk Ass. 14 275-295. [ Links ]
LOVE D, TWOMLOW S, MUPANGWA W, VAN DER ZAAG P and GUMBO B (2006) Implementing the millennium development food security goals challenges of the southern African context. Phys. Chem. Earth 31 731-737. [ Links ]
LOVE D, UHLENBROOK S, CORZO-PEREZ G and TWOMLOW S (in press) Rainfall-interception-evaporation-runoff relationships in a semi-arid meso-catchment, northern Limpopo Basin, Zimbabwe. Hydrol. Sci. J. 55 (5). [ Links ]
MAGOMBEYI MS and TAIGBENU AE (2008) Crop yield risk analysis and mitigation of smallholder farmers at quaternary catchment level: Case study of B72A in Olifants river basin, South Africa. Phys. Chem. Earth 33 744-756. [ Links ]
MAKARAU A and JURY MR (1998) Predictability of Zimbabwe summer rainfall. Int. J. Clim. 17 1421-1432. [ Links ]
MAPEDZA E, WRIGHT J and FAWCETT R (2003) An investigation of land cover change in Mafungautsi Forest Zimbabwe using GIS and participatory mapping. Appl. Geogr. 23 1-21. [ Links ]
MASON SJ, JURY MR and TYSON PD (1997) Climatic variability and change over southern Africa: a reflection on underlying processes. Prog. Phys. Geogr. 21 23-50. [ Links ]
MAZVIMAVI D (2004) Estimation of Flow Characteristics of Ungauged Catchments: Case Study in Zimbabwe. Ph.D. dissertation, Wageningen University and International Institute for Geo-Information Science and Earth Observation, Delft. 188 pp. [ Links ]
MAZVIMAVI D (2008) Investigating possible changes of extreme annual rainfall in Zimbabwe. Hydrol. Earth Syst. Sci. Disc. 5 1765-1785. [ Links ]
MILLY PCD, DUNNE KA and VECCHIA AV (2005) Global pattern of trends in streamflow and water availability in a changing climate. Nat. 438 347-350. [ Links ]
MILLY PCD, BETANCOURT J, FALKENMARK M, HIRSCH RM, KUNDZEWICZ ZW, LETTENMAIER DS and STOUFFER RJ (2008) Stationarity is dead: whither water management? Sci. 319 573-574. [ Links ]
MUPANGWA W, LOVE D and TWOMLOW S (2006) Soil-water conservation and other rainwater harvesting strategies in the semi-arid Mzingwane Catchments Limpopo Basin Zimbabwe. Phys. Chem. Earth 31 893-900. [ Links ]
MWENGE KAHINDA J-M, ROCKSTRÖM J, TAIGBENU AE and DIMES J (2007) Rainwater harvesting to enhance water productivity of rainfed agriculture in the semi-arid Zimbabwe. Phys. Chem. Earth 32 1068-1073. [ Links ]
MYERS JL and WELL AD (2002). Research Design and Statistical Analysis ( 2nd edn.). Lawrence Erlbaum, Philadelphia, 760 pp. [ Links ]
NEW M, HEIWITSON, B STEPHENSON DB, TSIGA A, KRUGER A, MANHIQUE A, GOMEZ B, COELHO CAS, MASISI DM, KULULANGA E, MBAMBALALA E, ADESINA F, SALEH H, KANYANGA J, ADOSI J, BULANE L, FORTUNATA L, MDOKA ML and LAJOIE R (2006) Evidence of trends in daily climate extremes over southern and west Africa. J. Geophys. Res. 111 D14102, DOI:10.1029/2005JD006289 [ Links ]
NYABEZE WR (2004) Estimating and interpreting hydrological drought indices using a selected catchment in Zimbabwe. Phys. Chem. Earth 29 1173-1180. [ Links ]
PETTITT AN (1979) A non-parametric approach to the change-point problem. Appl. Stat. 28 126-135. [ Links ]
RICHARDSON CJ (2007) How much did droughts matter? Linking rainfall and GDP growth in Zimbabwe. Afr. Aff. 106 463-478. [ Links ]
SAMAKANDE I, SENZANJE A and MANZUNGU E (2004) Sustainable water management in smallholder irrigation schemes: understanding the impact of field water management on maize productivity on two irrigation schemes in Zimbabwe. Phys. Chem. Earth 29 1075-1081. [ Links ]
SHAO Q, LI Z and XU Z (2009) Trend detection in hydrological time series by segment regression with application to Shiyang River Basin. Stoch. Environ. Res. Risk. Assess. DOI: 10.1007/s00477-009-0312-4. [ Links ]
SHUMBA EM and MAPOSA R (1996) An evaluation of the performance of six smallholder irrigation schemes in Zimbabwe. Irrig. Drain. Sys. 10 355-366. [ Links ]
STIGE LC, STAVE J, CHAN K-S, CIANNELLI L, PETTORELLI N, GLANTZ M, HERREN HR and STENSETH NC (2006) The effect of climate variation on agro-pastoral production in Africa. P. Natl. Acad. Sci. USA 103 3049-3053. [ Links ]
SURVEYOR GENERAL OF ZIMBABWE and FORESTRY COMMISSION (1996) 1:250000 Vegetation Map Series. Surveyor General of Zimbabwe, Harare, Zimbabwe. [ Links ]
THIRTLE C, ATKINS J, BOTTOMLEY P, GONESE N, GOVERE, HJ and KHATRI Y (1993) Agricultural productivity in Zimbabwe 1970-90. Econ. J. 103 474-480. [ Links ]
TILAHUN K (2006) The characterisation of rainfall in the arid and semi-arid regions of Ethiopia. Water SA 32 429-436. [ Links ]
TRENBERTH KE, JONES PD, AMBENJE P, BOJARIU R, EASTERLING D, KLEIN TANK A, PARKER D, RAHIMZADEH F, RENWICK JA, RUSTICUCCI M, SODEN B and ZHAI P (2007) Observations: Surface and Atmospheric Climate Change. In: Solomon SD, Qin M, Manning Z, Chen M, Marquis KB, Averyt M, Tignor and Miller HL (eds.) Climate Change 2007: The Physical Science Basis. Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge University Press, Cambridge. 235-336. [ Links ]
TU M, DE LAAT PJM, HALL MJ and DE WIT MJM (2005) Precipitation variability in the Meuse basin in relation to atmospheric circulation. Water Sci Technol. 51 5-14. [ Links ]
TWOMLOW S and BRUNEAU P (2000). Semi-arid soil water regimes in Zimbabwe. Geoderma 95 33-51. [ Links ]
TWOMLOW SJ, STEYN JT and DU PREEZ CC (2006) Dryland farming in southern Africa. In: Peterson GA, Ungwer PW and Payne WA (eds.) Dryland Agriculture (2nd edn.), American Society of Agronomy, Crop Science Society of America and Soil Science Society of America, Madison. 769-836. [ Links ]
TYSON PD (1986) Climatic Change and Variability in Southern Africa. Oxford University Press, Oxford. 220 pp. [ Links ]
UNGANAI LS and MASON SJ (2002) Long-range predictability of Zimbabwe summer rainfall. Int. J. Clim. 22 1091-1103. [ Links ]
VÖRÖSMARTY CJ, GREEN P, SALISBURY J and LAMMERS RB (2000) Global water resources: vulnerability from climate change and population growth. Sci. 289 284-288. [ Links ]
WOLTER K (2007) Multivariate ENSO Index National Oceanic and Atmospheric Administration URL: http://wwwcdcnoaagov/people/klauswolter/MEI/tablehtml (Accessed 15 December 2008). [ Links ]
WOLTER K and TIMLIN MS (1998) Measuring the strength of ENSO events how does 1997/98 rank? Weather 53 315-324. [ Links ]
YU B and NEIL DT (1993) Long-term variations in regional rainfall in the south-west of western Australia and the difference between average and high intensity rainfalls. Int. J. Clim. 13 77-88. [ Links ]
ZHANG X, ZHANG L, ZHAO J, RUSTOMJI P and HAIRSINE P (2008) Responses of streamflow to changes in climate and land use/cover in the Loess Plateau, China. Water Resour. Res. 44 W00A07, DOI:10.1029/2007WR006711 . [ Links ]
ZINYAMA L and WHITLOW R (2004) Changing patterns of population distribution in Zimbabwe. Geol. J. 13 365-384. [ Links ]
Received 7 January 2009; accepted in revised form 23 February 2010.