Global Resilience of Tropical Forest and Savanna to Critical Transitions

See allHide authors and affiliations

Science  14 Oct 2011:
Vol. 334, Issue 6053, pp. 232-235
DOI: 10.1126/science.1210657


It has been suggested that tropical forest and savanna could represent alternative stable states, implying critical transitions at tipping points in response to altered climate or other drivers. So far, evidence for this idea has remained elusive, and integrated climate models assume smooth vegetation responses. We analyzed data on the distribution of tree cover in Africa, Australia, and South America to reveal strong evidence for the existence of three distinct attractors: forest, savanna, and a treeless state. Empirical reconstruction of the basins of attraction indicates that the resilience of the states varies in a universal way with precipitation. These results allow the identification of regions where forest or savanna may most easily tip into an alternative state, and they pave the way to a new generation of coupled climate models.

Tree cover is one of the defining variables of landscapes, their ecological functioning, and their impact on climate. Despite insights into the effects of resource availability and disturbances on tree growth and survival (14), our understanding of the mechanisms determining global patterns of tree cover remains fragmented. A major question is whether tree cover will respond smoothly to climatic change and other stressors (5) or exhibit sharp transitions between contrasting stable states at tipping points (6). In some regions, forest, savanna, and treeless (barren or grassy) states have been suggested to represent alternative attractors (79). However, the case for multiple stable states of tree cover is largely based on models and on local observations of sharp transitions (69). Systematic studies of tree-cover distributions could help distinguish between hypotheses (1) but have been largely restricted to particular continents or biome types (46, 10, 11). To explore whether global patterns of tree abundance suggest gradual responses or, rather, alternative stable states, we analyzed the MODIS (Moderate-Resolution Imaging Spectroradiometer) set of remotely sensed estimates of tree cover in 1-km2 blocks in tropical and subtropical zones of Africa, Australia, and South America. We relate these patterns to precipitation (12), which is a major driver of past (13) and recent (14) shifts in the extension of tropical forests and savannas.

Taking the entire data set together, the frequency distribution of tree cover (Fig. 1A) is strikingly tri-modal (table S1), implying that there are three distinct underlying states: forest, savanna, and treeless. A closer look at tree distributions for different precipitation ranges (Fig. 1C) shows that although there is an overall increase of average tree cover with precipitation, the distinct character of the three states remains. Indeed, the characteristic tree cover of savanna (around 20%) and of forest (around 80%) remains remarkably constant over a wide range of rainfall levels. Thus, rather than a gradual increase in tree cover with precipitation, we see a shifting probability of being in either of the three distinct states. Defining cutoff levels of 5 and 60% tree cover for treeless and forest states, respectively (fig. S1), we can compute the probability of being in each of the states as a function of precipitation (Fig. 1B). Although continents differ widely in their rainfall range and forest abundance (Fig. 2), we find little difference between the continents in the way in which the probability of finding forest or savanna varies with precipitation (fig. S2), suggesting that the probability of finding the three alternative states varies in a fairly universal way with rainfall for tropical and subtropical regions. The most conspicuous deviation is a relative lack of forest in high-rainfall places in Australia, which is most likely explained by their relatively long dry season (15).

Fig. 1

Tree cover distribution in tropical and subtropical regions of Africa, Australia, and South America. (A) Tri-modal relative frequency distribution of tree cover (T) indicating three distinct underlying states: forest, savanna, and treeless, defined for further analyses as T ≥ 60%, 5% ≤ T < 60%, and T < 5%, respectively (12). (B) The probability of being in each of the states as a function of the mean annual precipitation. Bars represent the relative frequency of the three vegetation types in rainfall classes. Curves represent logistic regression models (table S3). (C) Frequency distributions of tree cover at different precipitation levels, illustrating the fact that distinct vegetation types persist despite the overall shift from treeless to forest with increasing precipitation. In (A) and (C), the percent of tree cover has been arcsine-transformed.

Fig. 2

Tree cover (untransformed) in 1-km2 grid cells as a function of the mean annual precipitation for (A) Africa, (B) Australia, (C) South America, and (D) intercontinental data sets [between 35°S and 15°N (12)]. Although the precipitation distribution and forest abundance vary between continents, the statistical relationships of tree cover to precipitation are quite similar (fig. S2).

The patterns suggest a double hysteresis of tree cover in response to rainfall (Fig. 3). This can be inferred in a formal way from the data, assuming the distribution of states to result from the interplay between stochastic processes and the dynamics of an underlying deterministic system (16). The diminishing frequency at which the three states occur toward the end of their precipitation range corresponds to a decrease in the size of the basin of attraction of the states toward the bifurcation points (B), where the stability ends through a collision with the unstable (dashed) equilibria. The resilience, defined as the capacity to recover from perturbations (17, 18), declines toward such tipping points. For instance, approaching the bifurcation point BF,S at which a forest-savanna transition is inevitable, the basin of attraction of the forest state shrinks while the attraction basin around the savanna state is growing. As a consequence, the critical loss of tree cover beyond which forest is expected to shift toward a savanna state becomes smaller in the vicinity of this bifurcation point.

Fig. 3

Relationship between the resilience of tropical forest, savanna, and treeless states and mean annual precipitation (in millimeters per year). (A) The tree cover data (percent, bottom plane) suggest a double catastrophe-fold. Stable states correspond to solid parts of the curve on the bottom plane and to minima in the stability landscapes. Unstable equilibria correspond to the dashed parts of the curve and to hilltops in the stability landscapes. At bifurcation points (B), stable equilibria disappear through collision with unstable equilibria. Resilience measured as the width of the basin of attraction around a stable state diminishes toward such bifurcation points. (B) Potential landscapes as computed directly from the data. Stable states (solid dots) are minima and the unstable equilibria (open dots) are maxima at a given precipitation level. A three-dimensional animation is available at (12).

The inferred shapes of the basins of attraction show how the resilience of the different states changes with precipitation. Another obvious indicator of resilience of a state for a given annual precipitation level is simply the fraction of sites that are in that state. Our logistic regression models predict that fraction as a function of rainfall (Fig. 1B). One straightforward application is the mapping of estimated resilience for current biomes. In this interpretation, a currently forested site, for instance, has a low resilience if, given its rainfall, the site would be much more likely to be in a savanna or treeless state. One can imagine that such forests most easily shift to a contrasting state with fewer trees in response to perturbations such as drought or logging. As an example of such a resilience map, consider the estimated forest resilience for the South American continent (Fig. 4). The low predicted resilience of forest in the central-eastern and southeastern Amazon basin coincides with human pressures in the arc of deforestation. These areas are precisely the regions with the highest risk of enhanced drought in future climate scenarios (19). For savanna, resilience decreases toward the drier end, where the chances of turning into a treeless state are larger, but also toward the wetter end, where a forested state seems a more likely alternative. At the latter sites it seems reasonable to assume that forests could be restored relatively easily. A full set of resilience maps for the three continents can be found in (12).

Fig. 4

Forest resilience for South America. (A) Resilience of remaining forest expressed as the probability of finding forest at the local mean annual precipitation level, computed with the global logistic regression model depicted in Fig. 2B. Forest with low resilience (yellow dots) is predicted to be most likely to turn into a savanna or treeless state. (B) Current distribution of tree density obtained from remote sensing (12). Resilience maps of forest, savanna, and treeless states for South America, Africa, and Australia can be found in (12).

To interpret these results in terms of vulnerability and management options it is important to understand the causes of patterns. Many studies have related tree cover to abiotic and biotic drivers (111) as well as social drivers (1922), all of which will probably affect the resilience of the different states. An important pattern is the conspicuous rarity of places with roughly 60 or 5% tree cover. Our interpretation is that these specific situations are unstable, and this can only be explained from positive feedbacks (17). Fire might play a role in the instability around 60% tree cover. It can promote the openness of savanna when grasses produce sufficient flammable fuel (3, 21, 23), but once tree cover becomes sufficiently dense, inhibition of grass growth and the resulting fires could lead to a self-propagating shift to a closed forest (3, 4, 9, 15). The instability at very low tree cover might be related to a positive feedback, promoting further tree colonization due to facilitative amelioration of the stressors that would prevent tree seedling establishment in an entirely grassy or barren state (1, 24).

The multistability suggested by our results has practical implications. For instance, deforestation to the unstable threshold of 60% tree cover might induce a self-propagating shift to an open savanna over a range of rainfall levels. On the other hand, exceptionally wet years related, for example, to strong El Niño–Southern Oscillation events could enhance the probability of a shift to higher tree cover (25), implying a potential window of opportunity for restoration efforts to help tip the balance (26). Insights into the particular ways in which drivers such as herbivores, fire, soil biota, and nutrients interact to affect tree dynamics (111, 23, 25, 26) will be essential to design practical guidelines for managing these complex systems.

Because tree cover can substantially enhance rainfall in some parts of the world (22, 2730), an interactive vegetation-climate system may sometimes have alternative stable states even if vegetation by itself does not (8, 22, 29, 30). Because our analysis refers to current rainfall distributions, it will not reveal such hysteresis resulting from vegetation-climate coupling. A corollary is that our analysis does not detect situations of low resilience related to a possible collapse of a rainfall pattern that depends on tree cover. However, coupled vegetation-climate models focusing on this issue have so far neglected the possibility of intrinsic vegetation hysteresis. Our results pave the way for a new generation of such models, combining local and regional nonlinearities. Constraining the character of the local hysteresis requires further work, because our current data combine sites that differ in aspects such as topography, rainfall variability, soil characteristics, and human pressures, thus enlarging the range of rainfall levels over which states are found.

Although more-detailed studies are essential to understand dynamics in different regions, our analysis shows that global patterns may be used to infer resilience. Such empirical approaches are essential because it is becoming clear that accurate mechanistic models to predict tipping points are currently beyond our reach (18), and determining the resilience of complex systems to critical transitions remains one of the most challenging problems in environmental science today (17, 18).

Supporting Online Material

SOM Text

Figs. S1 to S3

Tables S1 to S3

References (3133)

Movie S1

References and Notes

  1. See supporting material on Science Online.
  2. Acknowledgments: This research was partly funded by the European Research Council–Early Warning grant and Spinoza award received by M.S. The data reported in this paper are extracted as described in the supporting online material from the publicly available sites of MODIS ( and the Climatic Research Unit (
View Abstract

Navigate This Article