Similar non-random processes maintain diversity in two tropical rainforests.

Quadrat-based analysis of two rainforest plots of area 50 ha, one in Panama (Barro Colorado Island, BCI) and the other in Malaysia (Pasoh), shows that in both plots recruitment is in general negatively correlated with both numbers and biomass of adult trees of the same species in the same quadrat. At BCI, this effect is not significantly influenced by treefall gaps. In both plots, recruitment of individual species is negatively correlated with the numbers of trees of all species in the quadrats, but not with overall biomass. These observations suggest, but do not prove, widespread frequency-dependent effects produced by pathogens and seed-predators that act most effectively in quadrats crowded with trees. Within-species correlations of mortality with numbers or biomass are not found in either plot, indicating that most frequency-dependent mortality takes place before the trees reach 1 cm in diameter. Stochastic effects caused by BCI's more rapid tree turnover may contribute to a larger variance in diversity from quadrat to quadrat at BCI, although they are not sufficient to explain why BCI has fewer than half as many tree species as Pasoh. Finally, in both plots quadrats with low diversity show a significant increase in diversity over time, and this increase is stronger at BCI. This process, like the frequency-dependence, will tend to maintain diversity over time. In general, these non-random forces that should lead to the maintenance of diversity are slightly stronger at BCI, even though the BCI plot is less diverse than the Pasoh plot.     

Spatial scale resolves the niche versus neutral theory debate

The 50‐ha long‐term forest plot on Barro Colorado Island in Panama was ‘ground zero’ for the development of ecology's ‘neutral theory’ and comparisons with its ‘niche theory’ counterpart. In this issue, Garzon‐Lopez and colleagues used tree distributions at this site to recast the unresolved (and unresolvable) debate to show that observational scale drives the perception of which processes predominate.     

Stochastically driven adult–recruit associations of tree species on Barro Colorado Island

The spatial placement of recruits around adult conspecifics represents the accumulated outcome of several pattern-forming processes and mechanisms such as primary and secondary seed dispersal, habitat associations or Janzen–Connell effects. Studying the adult–recruit relationship should therefore allow the derivation of specific hypotheses on the processes shaping population and community dynamics. We analysed adult–recruit associations for 65 tree species taken from six censuses of the 50 ha neotropical forest plot on Barro Colorado Island (BCI), Panama. We used point pattern analysis to test, at a range of neighbourhood scales, for spatial independence between recruits and adults, to assess the strength and type of departure from independence, and its relationship with species properties. Positive associations expected to prevail due to dispersal limitation occurred only in 16% of all cases; instead a majority of species showed spatial independence (≈73%). Independence described the placement of recruits around conspecific adults in good approximation, although we found weak and noisy signals of species properties related to seed dispersal. We hypothesize that spatial mechanisms with strong stochastic components such as animal seed dispersal overpower the pattern-forming effects of dispersal limitation, density dependence and habitat association, or that some of the pattern-forming processes cancel out each other.     

Intensive research activity alters short-term seedling dynamics in a tropical forest

Researchers can have unintentional, yet significant effects on their study systems. We tested for the effects of an intensive tree census on seedling dynamics in a 50-ha permanent forest plot on Barro Colorado Island, Panama. At the community level, and for different shade-tolerance guilds, we found no significant differences in seedling recruitment or survival inside compared to controls outside the plot. However, among growth forms, canopy trees and lianas exhibited significantly lower seedling survival inside the plot. Results suggest that intense researcher activity impacts short-term vegetation dynamics, but effects do not accumulate over time.     

Moving beyond abundance distributions: neutral theory and spatial patterns in a tropical forest

Assessing the relative importance of different processes that determine the spatial distribution of species and the dynamics in highly diverse plant communities remains a challenging question in ecology. Previous modelling approaches often focused on single aggregated forest diversity patterns that convey limited information on the underlying dynamic processes. Here, we use recent advances in inference for stochastic simulation models to evaluate the ability of a spatially explicit and spatially continuous neutral model to quantitatively predict six spatial and non-spatial patterns observed at the 50 ha tropical forest plot on Barro Colorado Island, Panama. The patterns capture different aspects of forest dynamics and biodiversity structure, such as annual mortality rate, species richness, species abundance distribution, beta-diversity and the species–area relationship (SAR). The model correctly predicted each pattern independently and up to five patterns simultaneously. However, the model was unable to match the SAR and beta-diversity simultaneously. Our study moves previous theory towards a dynamic spatial theory of biodiversity and demonstrates the value of spatial data to identify ecological processes. This opens up new avenues to evaluate the consequences of additional process for community assembly and dynamics.     

Comparing Allee effect-based and dispersal-based neutral models for species abundance distribution patterns

Dispersal is important for biodiversity maintenance in both neutral and niche theories. However, little is known about the potential role of Allee effect at the community level. In the present study, we developed neutral models for quantifying the separate and joint influences of the Allee effect and dispersal process, respectively, on species abundance distribution (SAD) patterns. Tree census data from Barro Colorado Island (BCI), Panama were used as the case to compare different neutral SAD models. Results showed that Allee effects were not detected in the BCI tree SAD curve. By contrast, the neutral models with the incorporation of dispersal process (including both immigration and emigration) can remarkably improve the fitting power of neutral models on the BCI tree SAD curve. Finally, even though the influence is not detectable, the Allee effect-based SAD models still might be alternative SAD models for model comparison and null hypothesis testing.     

On plotting species abundance distributions

1. There has been a revival of interest in species abundance distribution (SAD) models, stimulated by the claim that the log-normal distribution gave an underestimate of the observed numbers of rare species in species-rich assemblages. This led to the development of the neutral Zero Sum Multinomial distribution (ZSM) to better fit the observed data. 2. Yet plots of SADs, purportedly of the same data, showed differences in frequencies of species and of statistical fits to the ZSM and log-normal models due to the use of different binning methods. 3. We plot six different binning methods for the Barro Colorado Island (BCI) tropical tree data. The appearances of the curves are very different for the different binning methods. Consequently, the fits to different models may vary depending on the binning system used. 4. There is no agreed binning method for SAD plots. Our analysis suggests that a simple doubling of the number of individuals per species in each bin is perhaps the most practical one for illustrative purposes. Alternatively rank-abundance plots should be used. 5. For fitting and testing models exact methods have been developed and application of these does not require binning of data. Errors are introduced unnecessarily if data are binned before testing goodness-of-fit to models.     

Effect of tree host species on fungal community composition in a tropical rain forest in Panama

Abstract.  To determine whether host species influence the composition of fungal communities, the ascomycetes and basidiomycetes present on three tree species ( Prioria copaifera (Fabaceae), Quararibea asterolepis (Bombacaceae), and Trichilia tuberculata (Meliaceae)) were sampled on the 50-ha Forest Dynamics Project plot in lowland moist tropical forest on Barro Colorado Island, Panama. The most abundant fungal morphotypes of both ascomycetes and basidiomycetes were generalists found on all three hosts, but detrended correspondence analysis revealed distinct differences in fungal community composition among host trees. These differences among hosts were constant across census years. Randomization tests revealed that there were significantly fewer host-generalist fungi than expected for ascomycetes but not for basidiomycetes. These results indicate that host composition plays a role in structuring both ascomycete and basidiomycete fungal communities, but that the most successful fungal morphotypes are capable of colonizing multiple host species.     

Age structure in neutral theory resolves inconsistencies related to reproductive-size threshold

Aims Neutral theory consists of a suite of models that assume ecological equivalence among individual organisms. They have been most commonly applied to tropical forest tree communities either as null models or as approximations. Neutral models typically only include reproductive adults; therefore, fitting to empirical tree community data requires defining a reproductive-size threshold, which for trees is usually set arbitrarily to a diameter at breast height (DBH) of 100 mm. The inevitable exclusion of some reproductive adults and inclusion of some saplings cause a non-random sampling bias in neutral model fits. Here, we investigate this problem and resolve it by introducing simple age structure into a neutral model.Methods We compared the performance and sensitivity of DBH threshold of three variants of a spatially explicit neutral model: the traditional model, a model incorporating random sampling and a model with two distinct age classes—reproductive adults and saplings. In the age-structured model, saplings are offspring from adults that disperse according to a Gaussian dispersal kernel around the adults. The only extra parameter is the ratio of adults to saplings, which is not a free parameter but directly measurable. We used species–area relationships (SARs) to explore the predicted effect of saplings on the species richness at different scales in our model. We then evaluated the three model variations to find the parameters required to maintain the observed level of species richness in the 50-ha plot on Barro Colorado Island (BCI). We repeated our analysis filtering the data at different minimum tree-size thresholds in order to find the effect this threshold has on our results. Lastly, we used empirical species–individual relationships (SIRs) to test the pre-existing hypothesis that environmental filtering is the primary cause of differences between the assemblage of saplings and that of adults on BCI.Important findings Our age-structured neutral model was characterized by SARs that were insensitive to the presence of saplings at large scales and highly sensitive to them at small scales. Both models without age structure were highly sensitive to the DBH threshold chosen in a way that could not be explained based on random samplings alone. The age-structured neutral model, which allowed for non-random sampling based on life stage, was consistent with species richness observations. Our analysis of empirical SIRs did not support environmental filtering as a dominant force, but it did show evidence for other differences between age classes. Age can now be easily incorporated into future studies of neutral models whenever there is a concern that a sample is not entirely composed of reproductive adult individuals. More generally, we suggest that modeling studies using tree data subject to a minimum size threshold should consider the sensitivity of their results to that threshold.     

Eight microsatellite markers for the neotropical tree Luehea seemannii (Tiliaceae)

We isolated eight polymorphic microsatellites from the neotropical tree Luehea seemannii for gene flow and genetic structure studies. We used a streptavidin subtractive enrichment technique to develop a library of CA/GT repeats. Eight loci were screened for diversity from 96 individuals from Barro Colorado Island (BCI) and neighbouring Gigante peninsula in Panama. Luehea seemannii shows moderate levels of genetic diversity within these two populations. Allelic richness ranged from four to nine alleles and averaged 6.44 alleles per locus. Average expected heterozygosity was 0.65 on BCI and 0.60 on Gigante. Results are compared to microsatellite data from another wind‐dispersed gap colonizing species common in Panama.     

Liana Abundance,Diversity, and Distribution on Barro Colorado Island,Panama

Lianas are a key component of tropical forests; however, most surveys are too small to accurately quantify liana community composition, diversity, abundance, and spatial distribution – critical components for measuring the contribution of lianas to forest processes. In 2007, we tagged, mapped, measured the diameter, and identified all lianas ≥1 cm rooted in a 50-ha plot on Barro Colorado Island, Panama (BCI). We calculated liana density, basal area, and species richness for both independently rooted lianas and all rooted liana stems (genets plus clones). We compared spatial aggregation patterns of liana and tree species, and among liana species that varied in the amount of clonal reproduction. We also tested whether liana and tree densities have increased on BCI compared to surveys conducted 30-years earlier. This study represents the most comprehensive spatially contiguous sampling of lianas ever conducted and, over the 50 ha area, we found 67,447 rooted liana stems comprising 162 species. Rooted lianas composed nearly 25% of the woody stems (trees and lianas), 35% of woody species richness, and 3% of woody basal area. Lianas were spatially aggregated within the 50-ha plot and the liana species with the highest proportion of clonal stems more spatially aggregated than the least clonal species, possibly indicating clonal stem recruitment following canopy disturbance. Over the past 30 years, liana density increased by 75% for stems ≥1 cm diameter and nearly 140% for stems ≥5 cm diameter, while tree density on BCI decreased 11.5%; a finding consistent with other neotropical forests. Our data confirm that lianas contribute substantially to tropical forest stem density and diversity, they have highly clumped distributions that appear to be driven by clonal stem recruitment into treefall gaps, and they are increasing relative to trees, thus indicating that lianas will play a greater role in the future dynamics of BCI and other neotropical forests.     

Can neutral theory predict the responses of Amazonian tree communities to forest fragmentation?

We use Hubbell's neutral theory to predict the impact of habitat fragmentation on Amazonian tree communities. For forest fragments isolated for about two decades, we generate neutral predictions for local species extinction, changes in species composition within fragments, and increases in the probability that any two trees within a fragment are conspecific. We tested these predictions using fragment and intact forest data from the Biological Dynamics of Forest Fragments Project in central Amazonia. To simulate complete demographic isolation, we excluded immigrants--species absent from a fragment or intact forest plot in the initial census but present in its last census--from our tests. The neutral theory correctly predicted the rate of species extinction from different plots as a function of the diversity and mortality rate of trees in each plot. However, the rate of change in species composition was much faster than predicted in fragments, indicating that different tree species respond differently to environmental changes. This violates the key assumption of neutral theory. When immigrants were included in our calculations, they increased the disparity between predicted and observed changes in fragments. Overall, neutral theory accurately predicted the pace of local extinctions in fragments but consistently underestimated changes in species composition.     

How likely is speciation in neutral ecology?

Patterns of biodiversity predicted by the neutral theory rely on a simple phenomenological model of speciation. To further investigate the effect of speciation on neutral biodiversity, we analyze a spatially explicit neutral model based on population genetics. We define the metacommunity as a system of populations exchanging migrants, and we use this framework to introduce speciation with little or no gene flow (allopatric and parapatric speciation). We find that with realistic mutation rates, our metacommunity model driven by neutral processes cannot support more than a few species. Adding natural selection in the population genetics of speciation increases the number of species in the metacommunity, but the level of diversity found in the Barro Colorado Island is difficult to reach.     

When Can Species Abundance Data Reveal Non-neutrality?

Species abundance distributions (SAD) are probably ecology’s most well-known empirical pattern, and over the last decades many models have been proposed to explain their shape. There is no consensus over which model is correct, because the degree to which different processes can be discerned from SAD patterns has not yet been rigorously quantified. We present a power calculation to quantify our ability to detect deviations from neutrality using species abundance data. We study non-neutral stochastic community models, and show that the presence of non-neutral processes is detectable if sample size is large enough and/or the amplitude of the effect is strong enough. Our framework can be used for any candidate community model that can be simulated on a computer, and determines both the sampling effort required to distinguish between alternative processes, and a range for the strength of non-neutral processes in communities whose patterns are statistically consistent with neutral theory. We find that even data sets of the scale of the 50 Ha forest plot on Barro Colorado Island, Panama, are unlikely to be large enough to detect deviations from neutrality caused by competitive interactions alone, though the presence of multiple non-neutral processes with contrasting effects on abundance distributions may be detectable.     

Size-Class Effect Contributes to Tree Species Assembly through Influencing Dispersal in Tropical Forests

We have investigated the processes of community assembly using size classes of trees. Specifically our work examined (1) whether point process models incorporating an effect of size-class produce more realistic summary outcomes than do models without this effect; (2) which of three selected models incorporating, respectively environmental effects, dispersal and the joint-effect of both of these, is most useful in explaining species-area relationships (SARs) and point dispersion patterns. For this evaluation we used tree species data from the 50-ha forest dynamics plot in Barro Colorado Island, Panama and the comparable 20 ha plot at Bubeng, Southwest China. Our results demonstrated that incorporating an size-class effect dramatically improved the SAR estimation at both the plots when the dispersal only model was used. The joint effect model produced similar improvement but only for the 50-ha plot in Panama. The point patterns results were not improved by incorporation of size-class effects using any of the three models. Our results indicate that dispersal is likely to be a key process determining both SARs and point patterns. The environment-only model and joint-effects model were effective at the species level and the community level, respectively. We conclude that it is critical to use multiple summary characteristics when modelling spatial patterns at the species and community levels if a comprehensive understanding of the ecological processes that shape species’ distributions is sought; without this results may have inherent biases. By influencing dispersal, the effect of size-class contributes to species assembly and enhances our understanding of species coexistence.     

Species-area curves, diversity indices, and species abundance distributions: a multifractal analysis

Although fractals have been applied in ecology for some time, multifractals have, in contrast, received little attention. In this article, we apply multifractals to the species-area relationship and species abundance distributions. We highlight two results: first, species abundance distributions collected at different spatial scales may collapse into a single curve after appropriate renormalization, and second, the power-law form of the species-area relationship and the Shannon, Simpson, and Berger-Parker diversity indices belong to a family of equations relating the species number, species abundance, and area through the moments of the species abundance-probability density function. Explicit formulas for these diversity indices, as a function of area, are derived. Methods to obtain the multifractal spectra from a data set are discussed, and an example is shown with data on tree and shrub species collected in a 50-ha plot on Barro Colorado Island, Panama. Finally, we discuss the implications of the multifractal formalism to the relationship between species range and abundance and the relation between the shape of the species abundance distribution and area.     

The impact of neutrality, niche differentiation and species input on diversity and abundance distributions

We present a spatially-explicit generalization of Hubbell's model of community dynamics in which the assumption of neutrality is relaxed by incorporating dispersal limitation and habitat preference. In simulations, diversity and species abundances were governed by the rate at which new species were introduced (usually called 'speciation') and nearly unaffected by dispersal limitation and habitat preference. Of course, in the absence of species input, diversity is maintained solely by niche differences. We conclude that the success of the neutral model in predicting the abundance distribution has nothing to do with neutrality, but rather with the species-introduction process: when new species enter a community regularly as singletons, the typical J-shaped abundance distribution, with a long tail of rare species, is always observed, whether species differ in habitat preferences or not. We suggest that many communities are indeed driven by the introduction process, accounting for high diversity and rarity, and that species differences may be largely irrelevant for either.     

The variation of tree beta diversity across a global network of forest plots

Aims With the aim of understanding why some of the world's forests exhibit higher tree beta diversity values than others, we asked: (1) what is the contribution of environmentally related variation versus pure spatial and local stochastic variation to tree beta diversity assessed at the forest plot scale; (2) at what resolution are these beta‐diversity components more apparent; and (3) what determines the variation in tree beta diversity observed across regions/continents? Location World‐wide. Methods We compiled an unprecedented data set of 10 large‐scale stem‐mapping forest plots differing in latitude, tree species richness and topographic variability. We assessed the tree beta diversity found within each forest plot separately. The non‐directional variation in tree species composition among cells of the plot was our measure of beta diversity. We compared the beta diversity of each plot with the value expected under a null model. We also apportioned the beta diversity into four components: pure topographic, spatially structured topographic, pure spatial and unexplained. We used linear mixed models to interpret the variation of beta diversity values across the plots. Results Total tree beta diversity within a forest plot decreased with increasing cell size, and increased with tree species richness and the amount of topographic variability of the plot. The topography‐related component of beta diversity was correlated with the amount of topographic variability but was unrelated to its species richness. The unexplained variation was correlated with the beta diversity expected under the null model and with species richness. Main conclusions Because different components of beta diversity have different determinants, comparisons of tree beta diversity across regions should quantify not only overall variation in species composition but also its components. Global‐scale patterns in tree beta diversity are largely coupled with changes in gamma richness due to the relationship between the latter and the variation generated by local stochastic assembly processes.     

Soil phosphorus heterogeneity promotes tree species diversity and phylogenetic clustering in a tropical seasonal rainforest

The niche theory predicts that environmental heterogeneity and species diversity are positively correlated in tropical forests, whereas the neutral theory suggests that stochastic processes are more important in determining species diversity. This study sought to investigate the effects of soil nutrient (nitrogen and phosphorus) heterogeneity on tree species diversity in the Xishuangbanna tropical seasonal rainforest in southwestern China. Thirty‐nine plots of 400 m² (20 × 20 m) were randomly located in the Xishuangbanna tropical seasonal rainforest. Within each plot, soil nutrient (nitrogen and phosphorus) availability and heterogeneity, tree species diversity, and community phylogenetic structure were measured. Soil phosphorus heterogeneity and tree species diversity in each plot were positively correlated, while phosphorus availability and tree species diversity were not. The trees in plots with low soil phosphorus heterogeneity were phylogenetically overdispersed, while the phylogenetic structure of trees within the plots became clustered as heterogeneity increased. Neither nitrogen availability nor its heterogeneity was correlated to tree species diversity or the phylogenetic structure of trees within the plots. The interspecific competition in the forest plots with low soil phosphorus heterogeneity could lead to an overdispersed community. However, as heterogeneity increase, more closely related species may be able to coexist together and lead to a clustered community. Our results indicate that soil phosphorus heterogeneity significantly affects tree diversity in the Xishuangbanna tropical seasonal rainforest, suggesting that deterministic processes are dominant in this tropical forest assembly.     

The Bushmeat Harvest Alters Seedling Banks by Favoring Lianas, Large Seeds, and Seeds Dispersed by Bats, Birds, and Wind

We evaluated predictions that hunters favor lianas, large seeds, and seeds dispersed by bats, small birds, and mechanical means for seedling banks in central Panama. We censused 3201 trees in 20 1-ha plots and 38,250 seedlings in the central 64 m² of each plot. We found significant differences in the species composition of the seedling bank between nine protected sites in the Barro Colorado Nature Monument and 11 hunted sites in the contiguous Parque Nacional Soberanía. Lianas, species with large seeds, and species with seeds dispersed by bats, small birds, and mechanical means were all overrepresented at hunted sites. The latter two findings could also be evaluated relative to the species composition of reproductively mature adults for canopy trees. The tree species present in the seedling bank had significantly heavier seeds than the tree species present as adults at hunted sites but not at protected sites. The representation of seed dispersal modes among the species present in the seedling bank did not reflect pre-existing differences in the local species composition of adults. We hypothesize that hunting large seed predators favors large seeds by reducing predation and increasing survival. We also hypothesize that the harvest of large birds and mammals that disperse many seeds favors other species whose seeds are dispersed by bats, small birds, and mechanical means. This process also favors lianas because the seeds of disproportionate numbers of liana species are dispersed by wind.