Skip to main content
  • Research article
  • Open access
  • Published:

Modelling the impact of climate change on woody plant population dynamics in South African savanna



In Southern Africa savannas climate change has been proposed to alter rainfall, the most important environmental driver for woody plants. Woody plants are a major component of savanna vegetation determining rangeland condition and biodiversity. In this study we use a spatially explicit, stochastic computer model to assess the impact of climate change on the population dynamics of Grewia flava, a common, fleshy-fruited shrub species in the southern Kalahari. Understanding the population dynamics of Grewia flava is a crucial task, because it is widely involved in the shrub/bush encroachment process, a major concern for rangeland management due to its adverse effect on livestock carrying capacity and biodiversity.


For our study we consider four climate change scenarios that have been proposed for the southern Kalahari for the coming decades: (1) an increase in annual precipitation by 30–40%, (2) a decrease by 5–15%, (3) an increase in variation of extreme rainfall years by 10–20%, (4) and increase in temporal auto-correlation, i.e. increasing length and variation of periodic rainfall oscillations related to El Niño/La Niña phenomena. We evaluate the slope z of the time-shrub density relationship to quantify the population trend. For each climate change scenario we then compared the departure of z from typical stable population dynamics under current climatic conditions. Based on the simulation experiments we observed a positive population trend for scenario (1) and a negative trend for scenario (2). In terms of the projected rates of precipitation change for scenario (3) and (4) population dynamics were found to be relatively stable. However, for a larger increase in inter-annual variation or in temporal auto-correlation of rainfall population trends were negative, because favorable rainfall years had a limited positive impact due to the limited shrub carrying capacity.


We conclude that a possible increase in precipitation will strongly facilitate shrub encroachment threatening savanna rangeland conditions and regional biodiversity. Furthermore, the negative effects found for positive auto-correlated rainfall support current ecological theory stating that periodically fluctuating environments can reduce population viability because species suffer disproportionately from poor environmental conditions.


In order to assess biodiversity response under climate change, Hannah et al. [1] recently emphasized the need to apply simulation models operating on a regional scale. Moreover, species respond differently to climate change because of different adaptations to their environment [2]. As a consequence, single-species models with a regional focus are essential to fully understand the manifold impact of global climate change. However, even though recent simulation tools have occasionally been applied for climate-sensitive animal species (e.g. [37]), spatial plant population models are extremely scarce (see review in [8]).

In this study we show how climate change may affect the long-term population dynamics of the raisin bush, Grewia flava DC, a common, fleshy-fruited, woody plant species of South African savannas. Understanding the population dynamics of Grewia flava is a crucial task, because it is widely involved in the shrub/bush encroachment process (e.g. [9, 10]). Shrub encroachment, i.e. the increase in woody plant cover, is a major concern for conservation and savanna rangeland management, due to its adverse effect on livestock carrying capacity (e.g. [1113]) and biodiversity (e.g. [14]). In the context of global climate change, increase in woody plant cover has primarily been investigated in association with elevated CO2 (e.g. [15, 16]). For example, Bond et al. [16] suggested that higher rates of atmospheric CO2 will have a positive effect on the post-fire regrowth of woody plants resulting in an increase in woody plant cover. Unfortunately, appropriate data verifying this hypothesis are not yet available.

Here, we present a different approach. We used Spatial Grewia Model (SGM) (see [17, 18]), a stochastic, spatially-explicit computer model to evaluate the impact of precipitation pattern change on Grewia flava population dynamics. In semi-arid and arid savannas, rain is the most important environmental parameter governing crucial life history processes in woody plants [19]. Hence, climate change related shifts in precipitation pattern will potentially have severe consequences for woody plant population dynamics. As a regional focus for our study we considered the southern Kalahari, the near centre of Grewia flava's distribution. In this area, recent climatological studies proposed either a decrease in mean precipitation of 5–15% by the year 2050 [20], or an increase by up to 30–40% (e.g. [21]. Further studies suggest an increase in the frequency and variability of extreme rainfall events (e.g. [22]), as well as alternating phases with low and high rainfall, typical for Southern Africa [23]. The large divergence between the various precipitation scenarios raises the question how woody plants would react along this spectrum. Therefore, we have set up a systematic sensitivity study to explore Grewia flava population dynamics along precipitation gradients, so as to detect possible thresholds in the response of this species to continuous variations of a crucial environmental driver.


By implementing relevant ecological processes of Grewia flava population dynamics into the SGM (see Figure 1 and model description in method section), we developed a standard scenario, based on annual rainfall for the period 1940 – 2000. This standard scenario led to stable population dynamics (see Figure 2 and method section). To compare population trends between four climate change scenarios we then used simple linear regression to calculate the z-value, i.e. the slope of the year – shrub density relationship.

Figure 1
figure 1

Flow chart of SGM showing causal pathways of Grewia flava population dynamics. Bold arrows indicate processes where population parameters and variables are affected by the annual type of rainfall (for detailed model description see text).

Figure 2
figure 2

Typical SGM time series of the standard scenario with 500 annual time steps. Population trend is given as shrub density with simple linear regression y = 0.0032x + 15.73. Mean slope of population trend for 100 simulation runs was z = 0.000049x + 13.41 with standard deviation of +/-0.017. Further model output from top to bottom: number of high (black) and low rainfall (grey) years per decade, fire occurrence (indicated as a black bar), density of fire- and drought-killed shrubs and density of recruits. Decades with frequent high rainfall years are indicated by frequent fires and increased fire mortality. Decades with exceptionally low rainfall show low fire frequency and increased drought mortality. Peaks of recruitment occur mostly in years with high rainfall and absence of fire.

(1) Increase in precipitation

An increase in the frequency of high rainfall years resulted in an increase of the z-value (Figure 3 and Figure 4). For example, a 20% increase in high rainfall years yielded a mean z-value of 0.01 (Figure 4). With an increase of 50%, the maximum z-value was 0.02. This upper limit could be clearly assigned, because the carrying capacity of cell type T had been reached, i.e. no further recruits were able to establish in the sub-canopy of trees. Even though additional juveniles may still emerge in cell type M, they are prevented from reaching the open due to a low probability of P matrix . Additionally, fire mortality is higher in cell type M.

Figure 3
figure 3

Typical time series for climate change scenarios with z-value of the population trend: increase in precipitation (1); decrease in precipitation (2); increase in inter-annual variation of precipitation (3); increase in temporal auto-correlation of precipitation (4) for period length PL with 20 and 250 years. Columns indicate 20%, 60% and 100% variation, respectively (see parameter values in Table 1).

Figure 4
figure 4

Results of simulation experiments for climate change scenarios according to: increase in precipitation (1), decrease in precipitation (2) and increase in inter-annual variation of precipitation (3) given as mean population trend z for 100 simulation replicates. The mean population trend is given as z for 100 simulation replicates plotted against percent variation in the probability of occurrence of extreme rainfall years as compared to the standard scenario.

(2) Decrease in precipitation

A decrease in mean precipitation resulted in a negative population trend (scenario 2, Figure 3) with a lower limit of z at -0.027 (Figure 4). A 10% increase of P low , corresponding to a probability of occurrence of low rainfall years of P low = 0.165, yielded a mean z value of -0.008. This continued to z = -0.02 for a 40% increase of P low . Above a 60% increase of P low 30% of all replicate runs resulted an extinction of the Grewia flava population within the 500 year time frame (not shown, however compare Figure 3). Comparing the upper and lower limits of the z- value for scenario (1) and (2), negative effects associated with an increase in drought events were more effective than the positive effects related to higher recruitment rates associated with comparable rainfall excess (0.02 as opposed to -0.027 see Figure 4).

(3) Increase in inter-annual variation of precipitation

An increase in inter-annual variation between low and high rainfall years below 100% did not cause any significant departure from stable population dynamics. However, above a threshold value of 100%, corresponding to a reduction of P average from 0.72 to less than 0.44 (see Table 1), we observed a decrease in the mean z-value. For example, for an increase in variation of 150%, z was -0.02 (Figure 4). Here, rainfall almost exclusively occurred as either high or low (P low = 0.375, P high = 0.325). Considering the negative and positive effects of low and high rainfall years, negative effects, such as drought, outweighed positive effects associated with enforced recruitment in high rainfall years.

Table 1 Probability values for low (first value) and high rainfall years (second value) for climate change scenarios: increase in precipitation (1); decrease in precipitation (2); increase in inter-annual variation of precipitation (3); increase in temporal auto-correlation of precipitation (4) with 'good' phase (4a) and 'poor' phase (4b). Columns indicate 20%, 60% and 100% variation in the probability of occurrence of extreme rainfall years, respectively, as compared to the standard scenario (P low = 0.15, P average = 0.72, P high = 0.13).

(4) Increase in temporal auto-correlation of precipitation

Despite constant mean and inter-annual variation in precipitation, increasing period length (PL) and intra-cycle variation (ICV) led to negative population trends and a decrease in z-value (Figure 5). Here, we found varying ICV thresholds for the PL scenarios tested, i.e. an abrupt departure from stable population dynamics. For example, a 10-year alternating rainfall cycle with an increase of 80 % in ICV (corresponding to a 5-year 'poor' phase with P low = 0.27 / P high = 0.026 and a 'good' phase with P low = 0.03 / P high = 0.234) yielded a z-value of -0.003, whereas population dynamics were stable for an ICV of 70%. Further thresholds were found above 60% ICV for a PL value of 20 years, 50% for 50 years, 40% for 100 years and 15% for 250 years. In the latter case, a large proportion of the population became extinct below 50% ICV where z was close to its lower limit. It is noteworthy, that we found no significant difference between similar PL scenarios that began with different phase conditions.

Figure 5
figure 5

Results of simulation experiments for increase in temporal auto-correlation of precipitation (climate change scenario 4) according to a PL value of 10, 20, 50, 100 and 250 years. The mean population trend is given as z for 100 simulation replicates plotted against percent variation in the probability of occurrence of extreme rainfall years as compared to the standard scenario.


Based on our simulation results, we have demonstrated that climate change predictions of precipitation pattern for the southern Kalahari may significantly affect Grewia flava population dynamics. This enables us to estimate the possible consequences for land use management and biodiversity due to the keystone function of Grewia flava in the shrub/bush encroachment process.

What are the implications of climate change related population trends for regional biodiversity and rangeland condition?

(1) Increase in precipitation

A climate change related increase in precipitation is commonly rejected [20]). However, based on a 28-year atmospheric circulation model, Mason et al. [21] proposed an increase in mean annual rainfall of 30–40% for the southern Kalahari. Yet, due to the model's insensitivity to topography, the authors also emphasize uncertainties in their predictions. This scenario would result in a strong increase in shrub density as indicated by a positive mean population trend for the SGM simulation period (Figure 4). Here, fire and drought mortality rates are too low to compensate substantial increase in juvenile recruitment associated with high rainfall years. On rangelands with domestic livestock this process may be strongly enhanced: in an earlier model version it was shown that cattle feeding on the foliage of Grewia flava may disperse seeds into the open matrix vegetation through dung deposition and thus facilitate shrub encroachment [18]. In this case we propose that an increase in precipitation will result in a catastrophic deterioration of current rangeland conditions within several decades. As a further consequence, biodiversity will most likely decrease due to the homogenization of woody plant community structure [24]. An increase in structural homogeneity due to shrub encroachment has been shown to adversely affect species diversity in both plants and animals (e.g. [14, 25]).

(2) Decrease in precipitation

The majority of climatological studies propose a decrease in annual precipitation of 5–15% by the year 2050 [20]. For a 10% increase in the probability of occurrence of low rainfall years (at the cost of average years) we found a negative trend for Grewia flava population dynamics (Figure 3). High frequency of low rainfall years resulted in an increase of drought-related mortality and the reduction of shrub density. Increases in the probability of occurrence of low rainfall years of more than 40% lead occasionally to population extinction within the SGM time frame. However, we believe that under natural conditions population dynamics of Grewia flava may also stabilize at a lower level, as suggested by its distribution in arid parts of the southern Kalahari where droughts are more frequent (e.g. [26]). Even though precipitation decrease may mitigate the risk of local Grewia encroachment, a strong reduction in shrub density may lead to cascade effects in the food chain because the fruits are an important component in the diet of animals and humans in an environment where food resources are otherwise scarce [27].

(3) Increase in inter-annual variation of precipitation

The inter-annual variation in precipitation, i.e. frequency and magnitude of extreme rainfall events, has been predicted to increase with ongoing climate change within the next decades (e.g. [22, 28]). Rates may range between 10–20% for the southern Kalahari [21]. For this scenario SGM produced no significant departure form stable population dynamics (Figure 4). Thus, we propose that the predicted increase in inter-annual variability will have a low impact on the natural population dynamics and shrub density of Grewia flava. However, this may differ for rangelands with high cattle grazing: an earlier model version showed that increase in the probability of extreme rainfall years can increase shrub encroachment if additional seed dispersal by cattle into the matrix is considered [18].

An increase in variability of more than 100% revealed a significant threshold behavior with a negative Grewia flava population trend. Above this threshold frequent droughts kill the offspring and inhibit emergence. Ecological theory generally suggest that population viability decreases with increase in environmental variability or stochasticity [29], because favorable events have a limited positive impact due to the carrying capacity, whereas unfavorable events have a full negative impact [7]. However, to our knowledge, the general character of this relationship has not been studied in detail yet. Thus, the pattern described here provides new insights into the relationship between environmental variability and the corresponding population dynamics.

(4) Increase in temporal auto-correlation of precipitation

Tyson [23] showed evidence of periodic, non-random rainfall oscillations with a period of ca. 18–20 years in Southern Africa. This distinct pattern has recently been associated with El Niño/La Niña phenomena and may increase under climate change [30]. Since we did not simulate actual total rainfall we used values of ICV, i.e. intra-cycle variation and PL, i.e. period length to evaluate the impact of climate change related increase in positive auto-correlation of rainfall. In our study region, Kruger [30] found an oscillation with a PL value of 22 years between 1955 and 1991. The ICV values used in our simulation experiments refer to an increase in the probability of an extreme rainfall category and are thus difficult to compare with the cyclic variability of actual total rainfall. However, we compared the ICV values from our model with departures of annual rainfall from the inter-annual mean of a typical rainfall oscillation provided by Kruger [30]. Based on this, we infer that an ICV value above 60% resembles a realistic rainfall variation between favorable and unfavorable periods of a 22-year oscillation. Although we did not find a negative impact on Grewia flava population dynamics for this scenario, a slight increase in ICV combined with an increasing length of El Niño/La Niña phases might have a significant negative impact on the population viability of Grewia flava (see Figure 5). This will have profound consequences for other organisms which depend on, for example, fruit provision.

As a further corollary, our results confirm current ecological theory that positive auto-correlation with constant average and inter-annual variation of a driving environmental parameter can reduce population viability and lead to extinction (e.g. [17, 31, 32]). In terms of the internal model processes of SGM, this means that negative effects in drought years (i.e. drought mortality and lack of reproduction) outweigh positive effects in wet years (i.e. increased reproduction). This is noteworthy, as mean annual rainfall and inter-annual variation was kept constant. Until now, this has been primarily shown in solely theoretical studies with under-compensating population dynamics (e.g. [33]) or individual-based animal models (e.g. [7]). To our knowledge this is the first empirical-based, plant population model that showed evidence of population persistence being negatively influenced by positive auto-correlation through a periodically fluctuating environmental driver.

Parameter sensitivity and predictive power of SGM

In order to assess the predictive power of SGM it is important to discuss the inference of actual precipitation amounts on probability values. Threshold values for extreme rainfall years are based on expert knowledge and may vary depending on the plant species considered. Moreover, percentage changes in extreme rainfall years, as applied in our model, are not necessarily equivalent to percentages changes in mean annual rainfall. For example, an increase in 'high' rainfall years by 10% represents an increase in mean annual rainfall of approximately 10%, depending on actual precipitation amounts in each 'high' rainfall year. Thus, it is important to derive principal population trends from the SGM model results rather than absolute predictions. Another issue regards the temporal dimension: long-term SGM simulation periods deviate from short-term predictions of climate change. However, population trends applied for only a few decades may yield biased results due to the cyclicity of Grewia population dynamics (see [17]). Finally, we stress that potential deviations in the annual probability of drought mortality may modify the output of the model. For example, a sensitivity analysis showed that a 100% increase in drought mortality for adult shrubs decreased the model output by 10% (see [18]).


In this study we have shown that climate change may have severe and sometimes unexpected implications at a regional species level. A woody plant like Grewia flava is relatively long-lived, fire and drought-resistant. Intuitively, this would suggest a low climate change impact on Grewia flava population dynamics. However, we found that despite Grewia flava's capability to survive in a harsh environment, it may be strongly affected when rainfall decreases as predicted, or increases in periodical fluctuations. This may also include possible range shifts in regional distribution which have not been studied here. Based on the model and the climate change scenarios analysed, it would be inappropiate to forecasts changes in the geographical distribution of Grewia flava, because SGM parameters have been estimated and validated for populations in the Kimberley region of the southern Kalahari. Moreover, current geographical distribution suggests that buffer mechanisms may facilitate survival in more arid parts of the southern Kalahari (e.g. through higher variability in emergence as an adaption towards lower frequency of high rainfall years).

Decrease in rainfall may reduce the severity of Grewia flava encroachment in the southern Kalahari. However, increase in rainfall will likely enhance this process. A shift from typical open grassland with solitary trees towards wide spread, homogeneous Grewia flava thickets implies negative cascade effects for other species, resulting in the loss of biodiversity [24]. As a further consequence, local rangelands are likely to be reduced in their carrying capacity of domestic livestock with negative effects for economic sustainability. Moreover, a shift in precipitation pattern may have further consequences, such as the alteration of fire regimes, grass biomass production or woody plant carbon uptake.


A detailed description of the simulation model SGM and previous results have been presented elsewhere (see [17, 18]). SGM has been previously validated with empirical data and simulates population dynamics of Grewia flava under specific land use, fire and rain scenarios in southern Kalahari semiarid savannas. However, to facilitate a better understanding of the current results we will briefly describe the study species and relevant aspects of the simulation model.

The study species

In the open savannas of the southern Kalahari, Grewia flava typically grows beneath the dominant tree species Acacia erioloba [34], because bird-mediated seed dispersal predominately confines new establishments to woody plant microsites. However, occasionally large individuals may be found in the open grassland matrix at former tree sites, suggesting high longevity of Grewia flava. Under high cattle grazing a substantial proportion of seeds may be distributed into the open matrix vegetation, since cattle feed on the foliage and fleshy fruits of Grewia flava [18]. This may result in a substantial increase in Grewia cover, particularly around boreholes (e.g. [10]). Grewia flava has excellent resprouting capabilities after fire [35] and low drought mortality rates [36]. Size class distributions suggest a demographic bottleneck in early life stages due to low rates of emergence and high juvenile mortality [36].

General model structure

The computer model SGM (see [17, 18]) represents a grid-based approach iteratively simulating population dynamics of Grewia flava in annual time steps for a period of 500 years. Based on empirical demographic and spatial data from the Kimberley region of the southern Kalahari (see [36]), an initial population of 15 shrubs ha-1 was distributed on a 200 × 200 cell grid, with each cell representing 5 m × 5 m of savanna vegetation. The SGM grid is developed in two layers: a landscape and a population layer. In the first layer micro-site types of the savanna vegetation may change in the course of time. In the second layer SGM simulates population dynamics of Grewia flava. A cell type was classified as either tree (T) (i.e. occupied by Acacia erioloba) or matrix type (M) (i.e. grassland vegetation). However, it switched from T to M status when a tree died after a mean life span of 200 years. The converse occurred where a new tree establishes. Initial tree distribution and spatial recruitment pattern was random with a constant density of 5 trees ha-1 for the entire simulation period. In each time step, SGM simulated important life history stages, environmental conditions and the key ecological processes of Grewia flava, the majority of which are governed by annual rainfall patterns (see Figure 1). For example in the model, rainfall directly determines the likelihood of a fire, drought and fire mortality, fruit crop size, emergence and annual shrub growth rate.


Annual rainfall in the southern Kalahari is highly variable, with precipitation ranging between 200 and 700 mm yr-1. Rain falls almost exclusively during summer (November to April) with an inter-annual mean of 417 mm for the period 1940 – 2000 (Kimberley airport, South African Weather Service 2001, unpublished data). For this period we defined a threshold value of 150 mm below and above the long-term mean to classify years into 'low', 'average' and 'high' rainfall years (thresholds of 267 mm yr-1and 567 mm yr-1, respectively) (see [18]). Accordingly, frequency of extreme rainfall years resulted in an annual probability of P low = 0.15 for low and P high = 0.13 for high rainfall years (P average = 0.72). Similar classifications have been applied in other studies, e.g. in a spatial simulation model of Acacia raddiana in the Negev Desert [37].

Fruit production and seed dispersal

Fruit production rates are based on mean crop size of Grewia size classes and may vary depending on the annual rainfall (see [18]). Seed dispersal is mostly zoochorous and spatially aggregated with a large proportion of seeds deposited in woody plant microsites. It is a crucial factor for the population dynamics and long-term viability of Grewia flava [36] and represented by two parameters: probability of seeds removed from a shrub and deposited in cell type M, P matrix , and cell type T, P acacia . P matrix varies randomly per year between 0.1–0.01% and represents occasional seed distribution through, e.g. small mammals. P acacia refers to bird-mediated seed dispersal with a range of 1–5% in high rainfall years, 2.5–7.5% in average and 5–10% in low rainfall years. Based on estimates from empirical data [36], the assumption of relatively constant seed dispersal rates may be reasonable as the proportion of seeds removed is most likely to be higher in years of lower fruit set. Fruit production typically varies more than bird abundances, and unless birds switch in exactly compensatory fashion to other fruits to the degree that Grewia flava becomes less abundant, fruit removal may be higher in years of low population size.


Even though microsite types differ markedly in micro-environmental conditions, they do not differ in emergence rates of Grewia flava seeds [36] suggesting a similar emergence probability for cell type M and T. However, depending on annual rainfall conditions, the probability of emergence in the model may vary between 1–2% in average and 2–4% in high rainfall years (no emergence occurs in low rainfall years). As survival of Grewia flava seeds in the soil is very low, seeds that do not emerge are assumed to die.

Fire and drought

The occurrence and intensity of a fire depends on the amount of rainfall, since precipitation determines annual grass biomass production and thus fuel load [38]. In the SGM, fire was simulated probabilistically for the entire grid and may only occur in high and average rainfall years. Average frequency of 7.9 years in the model is supported by fire intervals reported from similar savannah types [39]. The impact of fire varies spatially and demographically: a seedling in the matrix vegetation has a 95% chance of being killed in a fire, compared with 4% for adult shrubs which largely resprout in the following year. For T cells, fire mortality probability is 0% for adults and 75% for seedlings. These model assumptions are realistic, since grass fuel accumulation and fire severity beneath trees is less than in tree inter-spaces (e.g. [40, 41]). In the SGM, the probability of drought mortality varies between life stages and the annual rainfall conditions (see also model description in [18]). Annual drought mortality for adults is restricted to low rainfall years with a probability of 3% for both cell types. This assumption is based on data from O'Connor [42] and Schurr [36]. For seedlings the annual drought mortality probability is 90% for average and 50% for high rainfall years (see [18]). All seedlings are assumed to die in years with low rainfall conditions.


In the model, each shrub was assigned a size class with an average canopy cover. For definition of shrub size classes we used the canopy volume to group Grewia flava individuals into categories of small (Grewia S ; <1 m3), medium-sized (Grewia M ; 1–10 m3) and large plants (Grewia L ; >10 m3) as well as seedlings (Grewia seedling ; temporary state after emergence, transforming to Grewia S with the following year). We assumed a maximum age for each size class member with MaxAge S = 5 years for small and MaxAge M = 25 years for medium-sized shrubs. These estimates are based on annual shoot growth rates and aerial photograph analysis of the study area (see [18]). In each annual time step an individual can accumulate a growth year and, if MaxAge is reached, attain the next size class with Ptransition(S)= 0.2 for small and Ptransition(M)= 0.1 for medium-sized shrubs. No growth occurs in low rainfall years.

Carrying capacity and competition

To incorporate intra-specific competition we used a simple causal approach that incorporates cell-based shrub cover. Therefore, we defined a carrying capacity of K = 25 m2 as maximum total cell shrub cover. If shrub cover exceeded K, individuals in the cell died. Density-dependent mortality was simulated annually by removing the smallest individuals first, i.e. in descending order of size class, until shrub cover <K (for further details see [18]). Inter-specific competition between Grewia flava seedlings and the grass layer were neglected: empirical tests showed that emergence rates were similar within and outside of the grassy vegetation matrix (see [36]). For the adult stage, inter-specific competition with other woody plants is of low importance because Grewia flava often occurs as a mono-dominant species, particularly on bush encroached rangelands.

Simulation experiments

The standard scenario of the model, based on rainfall for the period 1940 – 2000, led to stable population dynamics (Figure 2). As standard deviation was generally high we performed 100 replicate simulation runs (± 0.0173 for the standard scenario, see Figure 2). To determine and compare population trends between the climate change scenarios we then used a simple linear regression of years vs. shrub density to calculate the z-value, i.e. the slope of the year – shrub density relationship. This is a reasonable approach for analyzing time-abundance relationships (e.g. [43]). When population dynamics were stable for 500 years, i.e. z was 0, the number of new recruits more or less resembled the number of fire- and drought-killed shrubs (Figure 2). Significant recruitment events mostly occurred in high rainfall years without fire.

Based on the standard scenario we performed simulation experiments for each of the four different scenarios of climate change with variation in inter-annual mean of precipitation, coefficient of inter-annual variation, and temporal auto-correlation (see Table 1):

(1) increase in precipitation, i.e. stepwise increase of P high resulting in a higher probability of high rainfall years (increase in recruitment events). Probability of low rainfall years is constant whereas frequency of average years is reduced, accordingly.

(2) decrease in precipitation, i.e. stepwise increase of P low resulting in a higher probability of low rainfall years (increase in drought events),

(3) increase in inter-annual variation of precipitation, i.e. stepwise increase of P low and P high resulting in a higher probability of low and high rainfall years at the cost of average years (unaltered inter-annual mean),

(4) increase in temporal auto-correlation of rainfall with 'good' and 'poor' phases.

For the scenario (4) we varied two parameters:

(4a) introduction of rain cycles with increasing period length PL, i.e. longer phases of favorable and unfavorable rain conditions and,

(4b) increase in intra-cycle variation ICV, i.e. increasing variation between alternating favorable and unfavorable periods within one rain cycle.

For PL we applied a period length of 10, 20, 50, 100 and 250 years, respectively, with each period subdivided into a similarly long 'good' and 'poor' phase. For example, a 10-year period length was subdivided into a 5-year period of superior and a 5-year period of poor conditions. For an increase in ICV we alternately increased and decreased high and low rainfall probabilities in each period, respectively (see 4a and 4b in Table 1). For example, for an ICV value of 100% P low equaled 0.00 and P high 0.26 in a 'good' phase, whereas values were 0.30 and 0.00, respectively for a 'poor' phase. Through this procedure inter-annual mean and variation of rainfall probabilities were identical for each scenario and the default set of rainfall probabilities. Combined increase in PL and ICV resulted in an increase of positive auto-correlation. Near-decadal epochs of above- and below-normal rainfall have been identified for the period 1955–1991 [30] and may increase in ICV under predicted climate change.


  1. Hannah L, Midgley GF, Lovejoy T, Bond WJ, Bush M, Lovett JC, Scott D, Woodward FI: Conservation of biodiversity in a changing climate. Conservation Biology. 2002, 16: 264-268. 10.1046/j.1523-1739.2002.00465.x.

    Article  Google Scholar 

  2. Erasmus BFN, Van Jaarsveld AS, Chown SL, Kshatriya M, Wessels KJ: Vulnerability of South African animal taxa to climate change. Global Change Biology. 2002, 8: 679-693. 10.1046/j.1365-2486.2002.00502.x.

    Article  Google Scholar 

  3. Forchhammer MC, Stenseth NC, Post E, Langvatn R: Population dynamics of Norwegian red deer: density-dependence and climatic variation. Proceedings of the Royal Society of London Series B-Biological Sciences. 1998, 265: 341-350. 10.1098/rspb.1998.0301.

    Article  CAS  Google Scholar 

  4. Peterson JT, Kwak TJ: Modeling the effects of land use and climate change on riverine smallmouth bass. Ecological Applications. 1999, 9: 1391-1404.

    Article  Google Scholar 

  5. Saether BE, Tufto J, Engen S, Jerstad K, Rostad OW, Skatan JE: Population dynamical consequences of climate change for a small temperate songbird. Science. 2000, 287: 854-856. 10.1126/science.287.5454.854.

    Article  CAS  PubMed  Google Scholar 

  6. Wang GM, Hobbs NT, Singer FJ, Ojima DS, Lubow BC: Impacts of climate changes on elk population dynamics in Rocky Mountain National Park, Colorado, USA. Climatic Change. 2002, 54: 205-223. 10.1023/A:1015725103348.

    Article  Google Scholar 

  7. Wichmann MC, Jeltsch F, Dean WRJ, Moloney KA, Wissel C: Implications of climate change for the persistence of raptors in arid savanna. Oikos. 2003, 102: 186-202. 10.1034/j.1600-0706.2003.12044.x.

    Article  Google Scholar 

  8. Kickert RN, Tonella G, Simonov A, Krupa SV: Predictive modeling of effects under global change. Environmental Pollution. 1999, 100: 87-132. 10.1016/S0269-7491(99)00089-5.

    Article  CAS  PubMed  Google Scholar 

  9. Skarpe C: Shrub layer dynamics under different herbivore densities in an arid savanna, Botswana. Journal of Applied Ecology. 1990, 27: 873-885.

    Article  Google Scholar 

  10. Moleele NM, Perkins JS: Encroaching woody plant species and boreholes: is cattle density the main driving factor in the Olifants Drift communal grazing lands, southeastern Botswana. Journal of Arid Environments. 1998, 40: 245-253. 10.1006/jare.1998.0451.

    Article  Google Scholar 

  11. Scholes RJ, Walker BH: An African savanna: synthesis of the Nylsvley study. 1993, Cambridge University Press, 306-

    Chapter  Google Scholar 

  12. Jeltsch F, Milton SJ, Dean WRJ, van Rooyen N: Analysing shrub encroachment in the southern Kalahari: a grid-based modelling approach. Journal of Applied Ecology. 1997, 34: 1497-1508.

    Article  Google Scholar 

  13. Roques KG, O'Connor TG, Watkinson AR: Dynamics of shrub encroachment in an African savanna: relative influences of fire, herbivory, rainfall and density dependence. Journal of Applied Ecology. 2001, 38: 268-280. 10.1046/j.1365-2664.2001.00567.x.

    Article  Google Scholar 

  14. Meik JM, Jeo RM, Mendelson JR, Jenks KE: Effects of bush encroachment on an assemblage of diurnal lizard species in central Namibia. Biological Conservation. 2002, 106: 29-36. 10.1016/S0006-3207(01)00226-9.

    Article  Google Scholar 

  15. Polley HW: Implications of rising atmospheric carbon dioxide concentration for rangelands. Journal of Range Management. 1997, 50: 562-577.

    Article  Google Scholar 

  16. Bond WJ, Midgley JJ: A proposed CO2-controlled mechanism of woody plant invasion in grasslands and savannas. Global Change Biology. 2000, 6: 865-869. 10.1046/j.1365-2486.2000.00365.x.

    Article  Google Scholar 

  17. Tews J, Moloney K, Jeltsch F: Modelling seed dispersal in a variable environment: a case study of the fleshy-fruited savanna shrub Grewia flava. Ecological Modelling. 2004, 175: 65-76. 10.1016/j.ecolmodel.2003.10.008.

    Article  Google Scholar 

  18. Tews J, Schurr F, Jeltsch F: Seed dispersal by cattle may cause shrub encroachment of Grewia flava on southern Kalahari rangelands. Applied Vegetation Science. 2004, 7: 89-102.

    Google Scholar 

  19. Scholes RJ, Archer SR: Tree-grass interactions in savannas. Annual Review of Ecology and Systematics. 1997, 28: 517-544. 10.1146/annurev.ecolsys.28.1.517.

    Article  Google Scholar 

  20. IPCC: Climate change 2001: The regional impacts of climate change. 2001, []

    Google Scholar 

  21. Mason SJ, Joubert AM: Simulated changes in extreme rainfall over southern Africa. International Journal of Climatology. 1997, 17: 291-301. 10.1002/(SICI)1097-0088(19970315)17:3<291::AID-JOC120>3.3.CO;2-T.

    Article  Google Scholar 

  22. Katz RW, Brown BG: Extreme events in a changing climate – variability is more important than averages. Climatic Change. 1992, 21: 289-302. 10.1007/BF00139728.

    Article  Google Scholar 

  23. Tyson PD: Climatic Change and Variability in Southern Africa. 1987, Oxford University Press, Capetown, 220-

    Google Scholar 

  24. Tews J, Brose U, Grimm V, Tielbörger K, Wichmann MC, Schwager M, Jeltsch F: Animal species diversity driven by habitat heterogeneity/diversity: the importance of keystone structures. Journal of Biogeography. 2004, 31: 79-92.

    Article  Google Scholar 

  25. Milton SJ, Dean WRJ: South Africa's arid and semiarid rangelands: Why are they changing and can they be restored?. Environmental Monitoring and Assessment. 1995, 37: 245-264. 10.1007/BF00546893.

    Article  CAS  PubMed  Google Scholar 

  26. Jeltsch F, Moloney KA, Milton SJ: Detecting process from snapshot pattern: lessons from tree spacing in the southern Kalahari. Oikos. 1999, 85: 451-466.

    Article  Google Scholar 

  27. Skinner JD, Smithers RHN: The Mammals of the Southern African Subregion. 1990, University of Pretoria, South Africa, Pretoria

    Google Scholar 

  28. Gordon HB, Whetton PH, Pittock AB, Fowler AM, Haylock MR: Simulated changes in daily rainfall intensity due to the enhanced greenhouse-effect – implications for extreme rainfall events. Climate Dynamics. 1992, 8: 83-102. 10.1007/BF00209165.

    Article  Google Scholar 

  29. Menges ES: Population viability analyses in plants: challenges and opportunities. Trends in Ecology and Evolution. 2000, 15: 51-56. 10.1016/S0169-5347(99)01763-2.

    Article  PubMed  Google Scholar 

  30. Kruger AC: The influence of the decadal-scale variability of summer rainfall on the impact of El Niño and La Niña events in South Africa. International Journal of Climatology. 1999, 19: 59-68. 10.1002/(SICI)1097-0088(199901)19:1<59::AID-JOC347>3.3.CO;2-2.

    Article  Google Scholar 

  31. Ripa J, Lundberg P: Noise colour and the risk of population extinctions. Proceedings of the Royal Society of London Series B-Biological Sciences. 1996, 263: 1751-1753.

    Article  Google Scholar 

  32. Heino M: Noise colour, synchrony and extinctions in spatially structured populations. Oikos. 1998, 83: 368-375.

    Article  Google Scholar 

  33. Ripa J, Heino M: Linear analysis solves two puzzles in population dynamics: the route to extinction and extinction in coloured environments. Ecology Letters. 1999, 2: 219-222. 10.1046/j.1461-0248.1999.00073.x.

    Article  Google Scholar 

  34. Dean WRJ, Milton SJ, Jeltsch F: Large trees, fertile islands, and birds in an arid savanna. Journal of Arid Environments. 1999, 41: 61-78. 10.1006/jare.1998.0455.

    Article  Google Scholar 

  35. Skarpe C: Observations on two bushfires in the western Kalahari, Botswana. Acta Phytogeographica Suecica. 1980, 68: 131-140.

    Google Scholar 

  36. Schurr F: Grewia flava in the southern Kalahari – The population dynamics of a savanna shrub species. 2001, University of Jena, MSc thesis, University of Jena, Jena, Germany, 89-

    Google Scholar 

  37. Wiegand K, Jeltsch F, Ward D: Analysis of the population dynamics of desert-dwelling Acacia trees with a spatially-explicit computer simulation model. Ecological Modelling. 1999, 117: 203-224. 10.1016/S0304-3800(98)00199-9.

    Article  Google Scholar 

  38. Jeltsch F, Weber GE, Grimm V: Ecological buffering mechanisms in savannas: A unifying theory of long-term tree-grass coexistence. Plant Ecology. 2000, 161: 161-171. 10.1023/A:1026590806682.

    Article  Google Scholar 

  39. van Wilgen BW, Biggs HC, O'Regan SP, Mare N: A fire history of the savanna ecosystem in the Kruger National Park, South Africa, between 1941 and 1996. South African Journal of Science. 2000, 96: 167-178.

    Google Scholar 

  40. Hochberg ME, Menaut JC, Gignoux J: The influence of tree biology and fire in the spatial structure of the West African savannah. Journal of Ecology. 1994, 82: 217-226.

    Article  Google Scholar 

  41. Jeltsch F, Milton SJ, Dean WRJ, van Rooyen N: Tree spacing and coexistence in semiarid savannas. Journal of Ecology. 1996, 84: 583-595.

    Article  Google Scholar 

  42. O'Connor TG: Impact of sustained drought on a semi-arid Colophospermum mopane savanna. African Journal of Range and Forage Science. 1999, 15: 83-91.

    Article  Google Scholar 

  43. Thomas L: Monitoring long-term population change: Why are there so many analysis methods?. Ecology. 1996, 77: 49-58.

    Article  Google Scholar 

Download references


This work was funded by the German Ministry of Education and Research (BMBF) in the framework of BIOTA South Africa (01LC0024). I would like to thank E. Bell and N. Hanan for useful comments on an earlier version of this manuscript.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Jörg Tews.

Additional information

Author's contribution

JT developed the SGM simulation model, conceived the study and drafted the manuscript. FJ supervised the study and the manuscript writing. Both authors read and approved the final manuscript.

Authors’ original submitted files for images

Rights and permissions

Reprints and permissions

About this article

Cite this article

Tews, J., Jeltsch, F. Modelling the impact of climate change on woody plant population dynamics in South African savanna. BMC Ecol 4, 17 (2004).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: