Report

The Importance of Demographic Niches to Tree Diversity

See allHide authors and affiliations

Science  07 Jul 2006:
Vol. 313, Issue 5783, pp. 98-101
DOI: 10.1126/science.1124712

Abstract

Most ecological hypotheses about species coexistence hinge on species differences, but quantifying trait differences across species in diverse communities is often unfeasible. We examined the variation of demographic traits using a global tropical forest data set covering 4500 species in 10 large-scale tree inventories. With a hierarchical Bayesian approach, we quantified the distribution of mortality and growth rates of all tree species at each site. This allowed us to test the prediction that demographic differences facilitate species richness, as suggested by the theory that a tradeoff between high growth and high survival allows species to coexist. Contrary to the prediction, the most diverse forests had the least demographic variation. Although demographic differences may foster coexistence, they do not explain any of the 16-fold variation in tree species richness observed across the tropics.

Comparative studies of tree demography typically consider the entire community as a unit, ignoring species differences (1), simply because most tree inventories include small samples of many species (2, 3). Comparative studies show that tropical forests typically have higher turnover than do temperate forests (4) and that higher tree turnover associates with higher tree diversity (5). These studies cannot, however, test ecological hypotheses about diversity, coexistence, and demography (610).

A tradeoff between rapid growth and long life span permits species coexistence and can foster diversity: Species reproducing early in life persist despite poor competitive ability by growing rapidly on disturbed sites where resources are abundant. Long-lived species coexist by outliving the weedy invaders, persisting where resources are scarce. This is a familiar and widely known tradeoff in plant and animal communities (911) called the successional-niche hypothesis (7, 12). At a deterministic equilibrium, an indefinite number of species can coexist by this mechanism, each differing from all others along a continuum from short life span (with high growth) to long life span (and low growth). With stochastic demography, however, there is limiting similarity and the equilibrium species richness is finite (11, 13). This hypothesis is widely quoted as an explanation for tropical forest diversity (1416). Here, we ask whether species differences along a demographic axis explain why some tropical forests have many more species than others.

If demographic niches are a key force controlling forest diversity, then more diverse forests have more demographic niches. More niches could come about either by spreading demographic rates over a wider range or packing more in the same range. Here, we focus on the first prediction: Tropical forests gain diversity by having a wider range of demographic niches, as reflected by the range of mortality and growth rates across species.

We provide a direct test by quantifying mortality and growth of 4500 tree species in 10 different forests in America, Asia, and Africa (17). The 10 sites form a large-scale observation program, spanning a wide range of environmental conditions, designed to provide species-specific information for little-known tropical trees (18). At each site, a 20- to 52-ha tree census was set out in extensive, largely undisturbed forest (table S1). Species richness within the census plots differed by 16-fold, from 73 species per 50 ha in a dry forest at Mudumalai, India, to 1167 species per 50 ha in a wet dipterocarp forest in Sarawak, Malaysia (19).

Past studies on the demography of individual tree species were based on direct measures of rate constants. These excluded many rare species because their rate estimates are subject to high error (20, 21). To overcome this limitation, we did not simply record species' rates of mortality and growth; instead, we quantified the distribution of demographic rates across the entire community. A hierarchical Bayesian approach accomplishes this with explicit probability models covering both the observations of individual trees within species and the variation among species; all species, including rare ones, are included. For mortality, within-species distributions were modeled with the binomial distribution; for growth, we chose the log-normal based on the tendency for individual growth rates within a species to be highly skewed to the right. By separating within-species variation, the hierarchical model allows focus on the question of how species differ (10, 22, 23).

At the community level, we had to describe the variation in species' demographic rates across species, and again, we used the log-normal to account for the skew to the right. Histograms of mortality rate m and growth rate g (24) are fitted well by the log-normal when rare species are excluded (Fig. 1). The log-normal requires two parameters, μ and σ, the mean and standard deviation of the natural logarithm of m or g, respectively. We were able to estimate values of μ and σ that best describe a community's demography with the use of the Gibbs sampler, simultaneously producing for every species estimates of mortality and growth rates that are adjusted for abundance. That is, for abundant species, the estimate is barely different from the observed rate, but for rare species, it is guided by the community-wide pattern (25).

Fig. 1.

Distribution of sapling demographic rates of all species in the Lambir plot. (A) Annual mortality, m, for all individuals with dbh = 10 to 99 mm. Filled bars show the histogram of observed mortality rates for the 746 species with ≥50 individuals; open bars add the 415 species with <50 individuals. The open bar at m = 0 extends off the graph (162 species had no mortality; 6 of these had <50 individuals). The horizontal axis is m, expressed as a percent. The solid line is the fitted log-normal, based on all 1161 species. The dashed vertical line indicates the mean of the logarithm of the fitted distribution (parameter μ, Table 1), which is very close to the median. (B) Annual growth, g, for individuals 10 to 49 mm in diameter. Filled bars are the histogram for 995 species with ≥10 individuals; open bars for the remaining 154 species. The solid curve and dashed line are the fitted log-normal and the mean of the logarithm, respectively, as in (A). Both histograms are curtailed at 8% to accentuate details where most of the species fall. The number of species above 8% is indicated by arrows.

Fitted log-normal distributions for the Lambir forest in Malaysia are plotted with observed histograms of sapling mortality and growth (Fig. 1). The fit is close for more common species (filled bars), demonstrating that the large number of zeroes in the mortality histogram are sampling artifacts in rare species. Growth rates are also spread by rare species, though not as conspicuously. Fitted distributions for all the forests can be compared graphically (Fig. 2) or with estimated 2.5 and 97.5 percentiles (Table 1). In the supporting online material, tables of mortality and growth rates for all species are provided (database S1).

Fig. 2.

Comparing the fitted distributions of sapling demography in four forests. (A) Annual mortality rate, m. (B) Annual growth rate, g. The lower end of the growth distribution in saplings is limited by measurement accuracy (30).

Table 1.

Variation in sapling mortality and growth rates across species in tropical forests. For mortality, all individuals with dbh = 10 to 99 mm were included; for growth, those with dbh = 10 to 49 mm were included. Species number refers to those with at least one 10- to 99-mm sapling alive at the outset of a given census interval. Under mortality are percentiles of the distribution of mortality rate parameters (m) across species: the median plus lower and upper percentiles (2.5 and 97.5) of the fitted log-normal. Similarly, under growth are percentiles for the distribution of growth rates (g) across species. Rates are expressed as percentages (100m or 100g)—that is, 5 = 5% = 0.05. For each of the percentiles, confidence limits are given, based on the Gibbs sampler (25). Information about the sites is presented in table S1.

SiteYearsNo. of speciesAnnual mortality (%)Relative growth (%)
MedianLowerUpperMedianLowerUpper
BCI 82-85 284 3.14 ± 0.46 0.38 ± 0.10 26.0 ± 6.9 2.84 ± 0.16 1.35 ± 0.14 6.0 ± 0.7
BCI 85-90 282 2.56 ± 0.37 0.31 ± 0.08 21.5 ± 5.5 2.41 ± 0.18 0.85 ± 0.12 6.8 ± 1.0
BCI 90-95 282 2.85 ± 0.43 0.32 ± 0.09 25.3 ± 6.7 2.15 ± 0.13 0.89 ± 0.09 5.2 ± 0.7
BCI 95-00 285 3.35 ± 0.42 0.48 ± 0.11 23.3 ± 5.7 1.97 ± 0.12 0.81 ± 0.09 4.8 ± 0.6
BCI 00-05 285 2.91 ± 0.41 0.40 ± 0.10 21.4 ± 5.5 2.10 ± 0.16 0.73 ± 0.10 6.1 ± 0.9
Yasuni 96-03 1077 1.55 ± 0.10 0.31 ± 0.04 7.7 ± 0.9 1.67 ± 0.04 0.83 ± 0.04 3.4 ± 0.2
La Planada 97-03 218 3.22 ± 0.47 0.45 ± 0.13 22.9 ± 5.9 2.30 ± 0.17 0.93 ± 0.13 5.7 ± 0.8
Pasoh 87-90 802 1.04 ± 0.06 0.36 ± 0.04 3.0 ± 0.3 2.25 ± 0.05 1.38 ± 0.05 3.7 ± 0.2
Pasoh 90-95 801 1.35 ± 0.07 0.42 ± 0.05 4.3 ± 0.4 1.59 ± 0.03 1.02 ± 0.04 2.5 ± 0.1
Pasoh 95-00 804 1.69 ± 0.09 0.47 ± 0.06 6.0 ± 0.6 1.55 ± 0.03 1.07 ± 0.03 2.3 ± 0.1
Lambir 92-97 1161 1.32 ± 0.07 0.34 ± 0.03 5.2 ± 0.4 1.57 ± 0.03 0.96 ± 0.03 2.6 ± 0.1
HKK 93-99 256 4.11 ± 0.57 0.70 ± 0.19 24.1 ± 6.6 4.83 ± 0.45 1.53 ± 0.26 15.2 ± 2.7
Mudumalai 88-92 56 13.06 ± 3.48 2.65 ± 1.42 64.4 ± 36 7.87 ± 1.26 4.65 ± 1.52 13.3 ± 3.8
Mudumalai 92-96 52 17.06 ± 6.43 2.35 ± 1.51 124 ± 113 6.35 ± 1.67 2.57 ± 1.38 15.7 ± 11
Mudumalai 96-00 39 7.96 ± 2.70 1.73 ± 1.13 36.6 ± 25 5.71 ± 1.53 2.42 ± 1.69 13.4 ± 9.8
Sinharaja 95-01 205 1.35 ± 0.17 0.30 ± 0.07 6.0 ± 1.3 1.38 ± 0.07 0.75 ± 0.07 2.5 ± 0.2
Edoro 94-00 368 1.43 ± 0.20 0.21 ± 0.06 9.6 ± 2.5 1.41 ± 0.09 0.57 ± 0.07 3.5 ± 0.5
Lenda 94-00 346 1.26 ± 0.19 0.18 ± 0.06 8.8 ± 2.2 1.06 ± 0.04 0.66 ± 0.05 1.7 ± 0.1

Most of the forests were dominated by species with sapling mortality rates near 1% year–1 (Fig. 2). Even the high-mortality forests, such as Barro Colorado Island (BCI) and Hua Khae Khaeng National Park (HKK), had modes close to 1% year–1 and low rates around 0.4% year–1. The main feature separating these high-rate forests from the low-rate sites (such as Pasoh in Malaysia) is the long tail reaching 20% year–1 mortality; at Pasoh, nearly all species had mortality rates of <3% year–1. Thus, forests fell broadly in two groups: BCI, HKK, and La Planada had upper sapling mortality rates above 20%, whereas Sinharaja, Lambir, Pasoh, and Yasuni had upper rates below 8%. The Congo sites had exceedingly low mortality stretching to a modest 10% at the upper end.

Distributions of growth rates were similar to distributions of mortality, but growth was about half as variable across species (Fig. 2). Conspicuously, sites with less variation in mortality also had less variation in growth (Fig. 3). These patterns held for larger trees, although mortality and growth rates were lower (table S2 and fig. S1).

Fig. 3.

Range of sapling demographic rates for tree species within a community versus the number of species at the site. The range is the logarithm of the ratio between the 97.5 and 2.5 percentiles of the fitted distributions (Table 1). The range for mortality is given by filled circles; the range for growth is given by open triangles. Sites can be identified by the number of species; for example, Lambir is the most diverse and farthest to the right. Multiple censuses at BCI, Pasoh, and Mudumalai are included, and in each case, fall in a tight group.

Examples from individual species help illustrate. At BCI, a fast-growing understory treelet, Palicourea guianensis, had a population of 376 saplings in 1982, and every single one had died by 2005 (nevertheless, the population grew to 851). Although Palicourea's mortality rate is infinite by direct calculation, the Gibbs sampler produces an estimate of 33% year–1. Alloplectus schultzei, a small, weedy treelet at La Planada, also suffered 33% year–1 mortality, losing 284 of 335 individuals over 6 years. In contrast, of 1161 species at Lambir, none had mortality of 30% year–1, and only two had rates above 20%; at Pasoh between 1987 and 2000, the very highest Gibbs-corrected mortality rate among 802 species was 14% year–1, in Melastoma malabathricum.

At the other end of the distribution, Cupania sylvatica, a midsized tree of the BCI understory, lost only 10 of 1102 individuals between 1990 and 1995 (0.23% year–1), and Carapa guianensis at La Planada, a large and valuable timber tree, lost only 11 of 894 (0.28% year–1). In three census intervals at Pasoh, the lowest mortality rate recorded was 0.34% year–1, in Cynometra malaccensis.

The Mudumalai forest stood out. Saplings had greatly elevated mortality and growth, with rates stretching much higher than any other site. Between 1988 and 1992, every species at Mudumalai had sapling growth of >6% year–1, and only BCI and HKK had many rates this high. At Lenda and Sinharaja, no species grew by 6% year–1. For larger trees (≥100 mm diameter), however, Mudumalai was in line with other forests, having modes of mortality and growth near 1% year–1 (table S2). Indeed, trees at Mudumalai had among the lowest rates as well as the highest: Anogeissus latifolia had 116 deaths out of 2179 trees from 1988 to 2000, whereas Kydia calycina had 1272 of 1328 trees die over the same interval, many because of elephant herbivory (26). Their rates differ by 50-fold: 27% year–1 in Kydia compared with 0.47% year–1 in Anogeissus.

Three of the sites with long tails of elevated mortality and growth—BCI, HKK, and Mudumalai—have intense annual dry seasons (table S1). Mudumalai and HKK also burn in some years (26) (other plots do not burn and none suffer large-scale wind damage). It was not surprising that annual drought elevated mortality. Many species at these sites, however, had exceedingly low rates of mortality and evidently did not suffer much from drought. Conversely, forests lacking the tail of high growth and mortality had no or modest annual dry seasons, including the three forests dominated by Dipterocarpaceae (Sinharaja, Pasoh, and Lambir). Seasonality, however, was not the only factor predicting high variation in demography; the ever-wet cloud forest at La Planada, Colombia, had a wide spread of growth and mortality, comparable to the seasonally dry sites.

Mudumalai and HKK have relatively open canopies compared with all of the other sites, a typical feature of dry forests, and many saplings at Mudumalai are sprouts from large root systems. These are likely reasons for elevated sapling growth at the two sites. In contrast, both BCI and La Planada have dense canopies and dark understories, so canopy openness does not obviously account for the high-growth species at those two sites.

Contrary to the prediction that demographic variability begets species richness, diverse forests had the least variation in demography (Fig. 3 and fig. S1). If anything, the most diverse forests had the fewest demographic niches. At Lambir, high species richness coupled with a low diversity of demographic rates meant that 135 tree species coexisting in close proximity had sapling mortality rates in a narrow window from 0.8 to 1.0% year–1.

We do not question that demographic variability plays some role in species coexistence. In American forests, the familiar genus Cecropia is found exclusively in small forest clearings (or outside the forest), where it rapidly colonizes and rapidly dies. The upper end of sapling mortality and growth distributions in America is set by gap specialists: C. obtusifolia at BCI (12% year–1 mortality, 14% year–1 growth), C. sciadophylla at Yasuni (5.0% year–1 mortality, 6.3% year–1 growth), and C. monostachya at La Planada (8.8% year–1 mortality, 8.2% year–1 growth). Diverse Southeast Asian forests lacked species with such high rates (27).

The most diverse tropical forests are the least diverse demographically. It remains plausible that demographic niches are packed more tightly in some forests than others, but this seems unlikely, because packing should depend only on population size and turnover, which do not vary much. Moreover, the successional-niche hypothesis is not favored by the strong peak in demographic rates near 1 to 2% year–1; if demographic niches were crucial, then rates ought to be spread evenly over the entire range (28). Instead, the similarity in demography of many species suggests trait convergence (29). We believe that broad diversity differences are due to the source pool of different biogeographic regions, and that demographic differences play a minor role in species coexistence.

Supporting Online Material

www.sciencemag.org/cgi/content/full/1124712/DC1

Materials and Methods

Fig. S1

Tables S1 to S4

References

Database S1

Computer Codes

References and Notes

View Abstract

Navigate This Article