Palms, peccaries and perturbations: widespread effects of small-scale disturbance in tropical forests
© Queenborough et al; licensee BioMed Central Ltd. 2012
Received: 19 October 2011
Accepted: 19 March 2012
Published: 19 March 2012
Disturbance is an important process structuring ecosystems worldwide and has long been thought to be a significant driver of diversity and dynamics. In forests, most studies of disturbance have focused on large-scale disturbance such as hurricanes or tree-falls. However, smaller sub-canopy disturbances could also have significant impacts on community structure. One such sub-canopy disturbance in tropical forests is abscising leaves of large arborescent palm (Arececeae) trees. These leaves can weigh up to 15 kg and cause physical damage and mortality to juvenile plants. Previous studies examining this question suffered from the use of static data at small spatial scales. Here we use data from a large permanent forest plot combined with dynamic data on the survival and growth of > 66,000 individuals over a seven-year period to address whether falling palm fronds do impact neighboring seedling and sapling communities, or whether there is an interaction between the palms and peccaries rooting for fallen palm fruit in the same area as falling leaves. We tested the wider generalisation of these hypotheses by comparing seedling and sapling survival under fruiting and non-fruiting trees in another family, the Myristicaceae.
We found a spatially-restricted but significant effect of large arborescent fruiting palms on the spatial structure, population dynamics and species diversity of neighbouring sapling and seedling communities. However, these effects were not found around slightly smaller non-fruiting palm trees, suggesting it is seed predators such as peccaries rather than falling leaves that impact on the communities around palm trees. Conversely, this hypothesis was not supported in data from other edible species, such as those in the family Myristicaceae.
Given the abundance of arborescent palm trees in Amazonian forests, it is reasonable to conclude that their presence does have a significant, if spatially-restricted, impact on juvenile plants, most likely on the survival and growth of seedlings and saplings damaged by foraging peccaries. Given the abundance of fruit produced by each palm, the widespread effects of these small-scale disturbances appear, over long time-scales, to cause directional changes in community structure at larger scales.
Disturbances are an important process structuring forests worldwide and have long been considered as significant drivers of dynamics and diversity [1–3]. In tropical forests, disturbances such as hurricanes and tree-falls from lightning or wind events create a mosaic of forest patches of different microhabitats at varying stages of succession [4–6] superimposed upon background topographical and soil variation. If niche partitioning permits coexistence along the axis of tolerance to disturbance, different species should be selected by such disturbances, and we would expect to find predictable suites of species associated with different disturbance regimes . This is indeed the case. Because of a trade-off between growth and survival, fast-growing pioneer species occur predictably in tree-fall gaps. Conversely, slow-growing shade-tolerant species survive well in closed canopy forest [8–10]. Furthermore, resprouting of trees is high in areas affected by hurricanes . These large-scale disturbances are obvious in the forest and are important mechanisms of species coexistence and community structure [4, 12]. However, the contribution of gaps and their associated pioneer flora to species coexistence is likely to be limited, because full tree-fall gaps are relatively infrequent in lowland rain forest  and most pioneer species occur at low abundance . Are there other kinds of disturbance in tropical forests affecting far wider areas of forest that influence species composition and community dynamics?
Several authors have suggested that smaller disturbances such as branch-falls have a greater and more widespread effect on forest dynamics and composition than might be apparent at first glance [13–15]. In the absence of large canopy-opening disturbances, mortality in sub-canopy layers can increase light availability, free up resources and lead to faster sapling growth [16, 17]. Furthermore, these sub-canopy gaps occupy large proportions of forest . Sub-canopy openings occur from a variety of processes, including the death of sub-canopy trees, large branch falls, or in-filling of canopy gaps by trees in middle forest strata . In addition to freeing resources, these events may cause disturbances to younger life stages, such as branch-falls affecting the survival and growth of seedlings and saplings [18–21]. Falling branches are heavy and can kill or physically damage small plants . If these branch-falls were prevalent enough of a selective force, species intolerant of this damage could be selected against in areas of high branch-fall. Peters et al.  extend this argument to apply to arborescent palm leaves, which are heavy enough to physically damage anything underneath (> 15 kg dry mass and several metres in length ). Given the ubiquity of large palms in western Amazonian forests (densities often exceeding 60 palms > 10 cm dbh ha-1[24–29]), these authors suggested that falling leaves could select for species adapted to withstand this disturbance impact, such as those with storage organs or cotyledons that could easily resprout following damage. Indeed, Peters et al.  found differences in sapling (0.5-2.5 m in height) communities in areas close to large arborescent palms compared to those with no influence of palm leaf fall. Although theirs was only a small study over 2.25 ha, it would seem that falling palm leaves could have a large effect on the species composition of forest communities.
However, falling branches and leaf fronds are not the only widespread small-scale disturbances in forests. Large frugivores such as tapirs and peccaries are known to cause heavy localised disturbance, dramatically increase mortality of seedlings and saplings [30–33], and consume and disperse seeds [34, 35]. Peccaries also consume palm fruit  and in a repeat study of Peters et al. , Beck et al.  showed that sapling communities were impacted by the presence of palms in forest containing peccaries, but not in forest where they had been hunted out. Beck et al.  attributed the relationship between palms and saplings to peccaries rooting around palms and not to falling palm fronds. However, both studies [15, 36] suffered from using only static data, examining distributions of palms relative to other species without the ability to test how these distributions develop through differential growth or mortality.
Whatever the mechanism (falling fronds or rooting peccaries), the disturbance around palms should cause differences in seedling and sapling density, growth, mortality, and/or re-sprouting close to palms compared to areas further away. If falling fronds are the cause, there should be similar effects in the neighbourhoods of large, fruiting palms as well as large, but non-fruiting palms. Alternatively, if the disturbance is caused by peccaries, we would not expect an effect by non-fruiting palms. In addition, although palms are a common food source, many other species of large-fruited tree are also consumed by peccaries . Therefore, we would also expect an effect on the seedling and sapling communities around other large fruit-bearing tree species similarly eaten by peccaries. In this latter case, the effect of these small-scale disturbances may be even more widespread than Peters et al.  envisaged.
We should also consider that palms may have positive effects on some species. Evidence gathered by Connell et al.  showed that smaller sub-canopy disturbances led to faster growth of understorey saplings, presumably because of increased light penetration or freed-up resources from the death of a sub-canopy individual. These sub-canopy gaps (without complete opening of the forest canopy) were very common in the forest (covering up to about 48% of forest). Falling palm fronds likely create such sub-canopy gaps and seedlings that survive palm or peccary damage may benefit from having (i) fewer competitors (seen in the lower local densities of seedlings and saplings), (ii) more light (because palms have a less dense canopy than other trees and because the falling fronds knock other branches down) and, (iii) more soil resources (for the same reasons). Ultimately, this increased availability of resources adjacent to palms could lead to higher growth or survival of saplings of some species.
In this study we use spatially-explicit data from long-term monitoring of seedlings, saplings in the neighbourhoods of large arborescent Arecaceae (mostly Iriartea deltoidea) and reproductive Myristicaceae trees in a large forest dynamics plot in western Amazonia to examine the potential structuring effect of falling leaves and branches on tropical forest communities. We aim to test (i) whether large palms and other fruiting trees appear to impact seedling and sapling dynamics and distributions and (ii) whether any difference is due to the trees themselves or an interaction with large frugivores. We extend previous studies of these phenomena in two key areas. First, we use dynamic data to understand the processes that lead to the patterns observed in static data in previous studies. Second, we consider that small-scale disturbances and/or interactions with frugivores must also apply to other fruiting trees. In this case, we examined seedlings and saplings in the neighbourhood of fruiting adults in a common dioecious family at our site, the Myristicaceae, fruits of which are also often consumed by peccaries . The Myristicaceae have been studied here for 10 years, and detailed information is available on the sex of each reproductive tree. Comparing the sapling communities around male (non-fruiting) and female (fruiting) trees allows a rigorous comparison of fruit-eating by peccaries, and a counterpoint to understanding the impacts of disturbances associated with arborescent palms.
Abundances of the four arborescent palm species on the Yasuní FDP in the first census
DBH class (cm)
Iriartea deltoidea Ruiz & Pav.
Astrocaryum chambira Burret
Oenocarpus bataua Mart.
Socratea exorrhiza (Mart.) H. Wendl.
Spatial patterns of distribution and density
Spatial patterns of community structure
Abundances of tree species with significantly high or low relative abundances in areas influenced by large arborescent palms on the Yasuní FDP in the first census
Abundance in quadrats:
Proportion in non-palm
Dynamic demographic effects of large trees on seedlings
Dynamic demographic effects of large trees on saplings
When we looked at the effects of mature Myristicaceae trees on sapling dynamics, we found fewer significant effects. Sapling survival near to male but far from female Myrsticaceae trees was significantly higher compared to saplings far from both sexes (Figure 5A, red squares). However, these same saplings also had a significantly higher probability of negative growth (Figure 5B, red squares). There was no effect of Myristicaceae trees on resprouting rates (Figure 5C).
Using a combination of spatially-explicit static and population dynamic data, we conducted multiple tests to examine hypotheses that falling palm fronds and/or peccaries are major structuring forces of lowland Neotropical rain forests [15, 36]. We found a restricted but significant direct effect of large arborescent mature fruiting palms on the spatial structure and species diversity of neighbouring sapling and seedling communities. There was also a strong and consistent larger-scale environmental effect of large palms being located in areas of lower neighbourhood density and diversity of saplings. Further, these patterns were corroborated by process in that we also found significant differences in the population dynamics of saplings (reduced growth and increased resprouting) and seedlings (lower survival) in the vicinity of mature palms. Evidence such as this was taken by previous authors to implicate falling palm fronds as major causal agents of mortality and selection in tropical forests . However, when we attempted to exclude the confounding effects of frugivores by analysing still-large but immature (ie. non-fruiting) palm trees, we found no significant spatial signature, although there is some evidence of reduced growth of neighbouring saplings. This would suggest that large frugivores such as peccaries, and not falling fronds, are the more important cause of physical damage and mortality of seedlings and saplings in the vicinity of mature palms as they root around for fruit , because we would expect falling fronds to have an equal effect around mature and immature palm trees.
We then tested these potential effects of peccaries by examining the spatial structure and population dynamics of seedlings and saplings around another important food source, the Myristicaceae. Most species of this family are large canopy trees  and produce abundant crops of fruit that are eaten by peccaries  among other frugivores [37, 38]. Unlike palms though, these trees do not have large heavy leaves that could damage seedlings and saplings beneath the canopy. Furthermore, the family is dioecious and so we could compare fruiting female trees with non-fruiting male trees. However, we found no convincing evidence that the seedling and sapling communities near female Myristicaceae were structured at all, and the significant population dynamic results we did find were in the opposite direction to that hypothesized.
Why did we find these conflicting results? One possibility is that peccaries are not large consumers of Myristicaceae fruit. It is known that peccaries consume a wide variety of species in tropical forests . However, the abundance of each species in peccary diets is less well known . Peccaries may focus largely on palms and eat other fruit as and when they happen upon it [32, 39]. Arboreal frugivores such as toucans and primates will also likely have a large effect on the availability of fruit on the ground. Further data on peccary and other frugivore diets are required. In addition, palm seeds likely remain attractive to peccaries (if not viable) long after abscission because of their hard seed coat. Myristicaceae seeds are not protected in this way and quickly decompose.
Alternatively, there may be a difference in the pattern of abscission of palm and Myrsticaceae fruit and the seed shadows these generate. Palm fruits are held in a single infructescence at the top of the tree. These fruit then fall in rapid succession, creating a very dense carpet of fruit around the tree with a short dispersal tail . Animals rooting around in this small area will likely have a big impact on the juvenile plants growing there . Myristicaceae trees, on the other hand, have a much larger canopy, the fruit carpet is thinner and more spread out. Animals foraging for fruit here will likely have a lower impact on juvenile plants.
A final possibility is our assumption that the fronds of mature and immature palms are equivalent does not hold, and that immature palm trees have lighter and/or smaller leaves that cause less damage to seedlings and saplings when they fall. We were not able to compare leaves from both types of palm trees to verify this, although Iriartea leaves do change morphology throughout their life cycle [23, 41]. Furthermore, leaves falling from shorter immature palms will have lower velocity and might do less damage than leaves falling from the greater height of mature palms. Finally, the size cut-offs determined by Iriartea for mature and immature individuals might be different for other species. However, these species are much less abundant than Iriartea and are unlikely to significantly influence our results.
Implications for forest diversity and community structure
Given the difficulties of reliably measuring seedling growth, it was not possible to repeat all three tests of dynamic data in both seedlings and saplings. Growth and resprouting were analysed solely in saplings. However, it seems that significant mortality of seedlings occurs in the vicinity of palm trees, but that saplings can generally withstand physical damage caused by peccaries and/or palm fronds, and rather suffer reduced growth and can potentially resprout. This difference in mortality implies that different mechanisms are at work at different life history stages and that if seedlings can grow large enough, their fate is more secure . The differences in spatial structure and sapling density around palms are therefore more likely to be formed at the seedling stage than the sapling stage. The small-scale repulsion of saplings (Figure 1B) may be a consequence of the immediate negative effect of large palms on seedlings (falling palm fronds cause increased mortality and stem breakage in neighboring saplings) and the large-scale effects (Figure 1A) could be a long-term cumulative population dynamics effect that results from many generations of palms having a negative impact on nearby seedlings.
The impacts of palms on the forest structure are of smaller spatial scale than that hypthesized by Peters et al. , who suggested that the effect of palms extended out to 7 m and that at least a quarter of lowland forest could be impacted by falling palms fronds per year. We found significant negative effects of palm trees out to only 3 m (in terms of over-dispersion of saplings [spatial pattern]) and 5 m (in terms of species repelled [ISAR]). Given the spatial pattern and density of palm trees in our plot, the area of forest predicted to be directly affected by this process is therefore reduced to 17%, although we noted above that these effects could multiply over generations leading to large areas of lower sapling density.
The negative impacts identified above appear limited to large palms and do not involve other fruiting trees such as the Myristicaceae. We would conclude therefore that this process is not a disturbance that affects large swathes of forest (or everywhere there is a tree that produces large peccary-consumed fruit). If the effect is limited to the environment of palms, this creates heterogeneity in the forest, which may lead to a wide variety of regeneration niches for species that can withstand the effects of falling leaves, fruit and rooting frugivores . Once established, however, these individuals may have better access to light and soil resources, given the diffuse nature of palm canopies and the large amounts of decaying fruit and animal faeces that likely accumulate around the base of palm trees. The stilt-roots of many arborescent palms probably increase nutrient availability by trapping leaves and organic matter.
This kind of widespread small-scale disturbance creates heterogeneity in regeneration niches in both space and time. Fruiting and frugivory by peccaries is episodic, and palms generally mature one infructescence at a time. Furthermore, in the aseasonal forest of Yasuní, Iriartea fruit production is not confined to a few months of the year (Nancy Garwood, unpublished data). Confoundingly, palm fronds are also produced regularly throughout the year. Beck et al.  found no difference in leaf-fall rates between two widely separated forests in Peru. In seasonal sites such as Yasuní, therefore, where seedling germination takes place throughout the year, some individuals may escape damage simply by growing enough before the next frond or infructescence falls. In more seasonal forests where germination is often restricted to short periods, individuals are of similar ages and/or vulnerabilities and so may suffer greater impact. A further test of our hypotheses could take place in seasonal forest and examine if seedlings and saplings suffer greater mortality during the fruiting season (therefore more likely caused by peccaries) or during the non-fruiting season (more likely leaf-fall).
The ubiquitous nature of the palm family in Neotropical forests implies that given some effect, palm trees could be a powerful force structuring these forests . Several species of tree dominate large areas of lowland forests in western Amazonia, and arborescent palms such as Iriartea are among this 'oligarch' group [24, 26, 28]. Widespread small-scale mosaics created by these palms may affect the composition of advanced regeneration for when a larger canopy gap does open. However, we note that peccaries are among the first large mammals to be hunted out of forest near to human habitation. Hunting will likely have dramatic effects on forest structure .
If reduced survival and growth around palms is predictable and selective, what can we hypothesize about species that occur in the neighbourhood of palms? Species that are better able to survive damage by falling fronds or rooting peccaries would be more likely to survive in areas close to palm trees. From our community-level analyses we cannot identify traits of individual species. However, species with underground storage organs or cotyledons might find an advantage in germinating near to a large palm because of their increased resource availability allowing resprouting following damage . We identified a number of species that had significantly higher or lower abundance around palms - identification of key functional traits of these species implicated by this pattern would be useful. An advantage might even apply to the seedlings of arborescent palm trees themselves, if they can survive predation and falling fronds. Iriartea is one of the most widespread and abundant tree species in western Amazonia. If Iriartea seedlings can survive better in the adult zone of influence than the majority of non-palm species, this may go some way to explaining how such abundant populations can be maintained. Future work should therefore look for more evidence of selection potential of the process that Peters et al. , Beck et al.  and we have identified. What species are selected for? Does species diversity change predictably from seedlings to saplings? Are there species that appear more frequently together in proximity to palms than in other parts of the forest? Furthermore, a broader-scale comparison with other forests may be helpful. For example, the forests of SE Asia are not dominated by palms as are Neotropical ones, and yet peccaries have a significant impact on seedling and sapling communities .
Palms themselves are not immune to the effects of branch falls  and suffer negative density dependence , caused by falling fronds [22, 46–48], peccaries [32, 34], and pathogenic fungi . Although palm trees significantly impact non-palms, the combination of falling fronds and peccary seed predation and fungal attack inhibit population growth such that palms may currently be at maximum density in the forest.
On a much larger scale then previously tested, we examined hypotheses that peccaries and/or falling palm fronds are major structuring forces of lowland Neotropical rain forests. We found a restricted but significant direct effect of large arborescent mature fruiting palms on the spatial structure, population dynamics and species diversity of neighbouring sapling and seedling communities. These effects were not found in smaller non-fruiting palm trees. However, this evidence of a primary effect of peccaries on the seedling and sapling communities around palms was not supported in data from other edible species in the Myristicaceae.
Given the abundance of arborescent palm trees in Amazonian forests, it is reasonable to conclude that their presence does have a significant, if spatially-restricted, impact on juvenile plants. The widespread effects of these small-scale disturbances appear, over long time-scales, to cause directional changes in community structure. However, the population growth of large arborescent palms themselves is very likely limited by the very same two factors that impact their surrounding juvenile plant communities.
Tree and sapling data
A 25-ha permanent forest dynamics plot (FDP; http://www.ctfs.si.edu) is located inside Yasuní National Park (0°41S, 76°24W; ), a still largely wilderness, tropical lowland aseasonal rain forest in eastern Ecuador . Mean annual rainfall is approximately 2800 mm and total monthly rainfall is almost never ≤ 100 mm. Mean monthly temperature is 25-27°C . The FDP ranges from 216 to 248 m a.s.l: it includes two ridges and an intervening valley that floods for brief periods. From 1995-1999, all freestanding stems ≥ 1 cm diameter at breast height (dbh at 1.3 m), excluding lianas, were tagged, mapped and identified . In the first census, 1104 morph species were recorded, comprising 151,230 individual trees . The FDP was recensused in 2002-2003 .
There are four species of arborescent palm in the Yasuní FDP: Astrocaryum chambira Burret, Iriartea deltoidea Ruiz & Pav., Oenocarpus bataua Mart., and Socratea exorrhiza (Mart.) H. Wendl. The majority of individuals are Iriartea (Table 1). Although the relationship between stem diameter and height in palms is not as robust as that in dicotyledon trees, Iriartea generally reaches maturity around 15 cm dbh . Because the height of palms is not measured on the FDP, we took this diameter as the cut-off between mature and immature large palms. In order to examine the interaction between frugivory and leaf-falls, we compared mature fruiting palms (> 15 cm dbh) to large palms without fruit, being those 10-15 cm dbh.
The size at which tropical trees mature is difficult to examine and often unknown for many species and can depend on a number of factors such as light availability, nutrient and water resources. However, the Myristicaceae at Yasuní have been studied intensively for the previous 10 years, and trees have been visited multiple times to look for signs of reproduction. The family is dioecious and all reproductive (ie. mature) male and female trees have been identified and sexed from flowers and/or fruit [54–58]). The environment below a male and female Myristicaceae adult should be qualitatively similar, allowing the separation of the effect of fruit crop and associated frugivores from the effects of light availability or habitat quality on the performance of seedlings and saplings.
Within the 25-ha FDP, 339 1 m2 seedling plots were established along trails in 2002 . They were arranged around seed traps in sets of three as part of a long-term community phenology study. Survival, growth and recruitment are monitored every year and we use data from 2003 to 2008.
We conducted two sets of analyses. First we looked at the static spatial patterns of saplings and trees in the FDP. Second we examined the demographic rates of seedlings and saplings to understand the causes of these patterns. In all analyses, we tested the previously identified potential drivers of spatial patterns and demography. These were (i) large mature arborescent palms > 15 cm dbh, (ii) smaller, immature arborescent palms 10-15 cm dbh, (iii) large mature female Myristicaceae trees, (iv) large mature male Myristicaceae trees. In comparisons examining palms, differential responses based on proximity to larger or smaller palms are likely due to the effects of fruiting, while similar responses could be attributed to large, falling fronds common to both groups. Comparisons between male and female Myristicaceae isolate the effects of a large fruit crop and may generalize impacts of peccaries, if any exist, of the large palms. Unless otherwise stated, all analyses were conducted with R .
Spatial patterns of distribution and density
In order to examine patterns of attraction and repulsion between two sets of spatially-explicit points, we used an extension of the O-ring statistic for bivariate patterns  to test, for example, the density of saplings (1-2 cm dbh) at various distances r around large palms (> 15 cm dbh). More formally, the O-ring statistic O 12 (r) for a bivariate pattern composed of type 1 and type 2 points yields the average density of type 2 point within rings of radius r and width dr around the points of pattern 1. The O 12 (r) is also called neighbourhood density function and is closely related to the pair correlation function g 12 (r), i.e., O 12 (r) = λ 2g12(r)[61, 62] where λ2 is the average density of type 2 points in the plot.
To find out if the distribution of saplings in the neighbourhood of palms or Myristicaceae trees is different from a distribution expected by chance for each of these analyses we initially tested two different null models of spatial distribution. The first null model retained the spatial distribution of saplings and assumed a random distribution of palms (i.e., a homogenous Poisson null model ). This null model therefore compares the observed neighbourhood density of saplings around palms with sapling densities in random locations. If the neighbourhood density of saplings around focal trees is smaller (or larger) than expected for randomly relocated focal trees we can infer a negative (or positive) impact of focal trees on saplings. However, we could observe the same pattern if (larger-scale) environmental heterogeneity acted such that the focal trees were located in areas of higher or lower overall sapling density. We therefore used a second null model that approximately factors out such larger-scale environmental variability and reveals the "pure" smaller-scale effects of focal trees on saplings (e.g. [63–65]). The second test also assumed a random distribution of palms, but the palms were only displaced within a given neighborhood of 30 m. The neighborhood of 30 m is somewhat larger than the expected range of direct interactions between focal trees (i.e., palms and adult Myristicaceae trees) and saplings (< 10 m based on ePeters2004). Technically, this null model was implemented as a heterogeneous Poisson process with a non-parametric kernel estimate of the intensity function λ1(x) of the focal trees (for details see ).
We made these comparisons for (i): non-palm sapling distributions versus mature and immature palms and (ii) non-palm sapling distributions versus male and female Myristicaceae. All analyses were run for 199 simulations and 95% simulation envelopes were calculated. These analyses were conducted using the Programita software package .
Spatial patterns of community structure
While the previous analysis focused on overall sapling density we also carried out tests where we compared the species diversity of sapling communities that were influenced by palms to sapling communities that were not. For the first two tests we compared sapling (1-2 cm dbh) communities in 5 × 5 m quadrats that either contained ≥ 1 palm or were adjacent to a quadrat containing a palm to quadrats that were > 1 quadrat distant from a palm-quadrat. This accounted for the up to 10 m distance over which palm fronds fall. In the third test, we examined the individual species-area relationships between palms and saplings.
First, we tested whether sapling density, species richness and rarefied species richness per two individuals were lower in quadrats influenced by palms, using the 'vegan' R package . This provided a basic description of the community spatial structure with minimal consideration of species identity.
Second, we compared the community composition of saplings in quadrats using a multiple response permutation procedure (adonis), in 'vegan' . This method is analogous to a multivariate analysis of variance and provides a robust test of whether there is a significant difference between the species composition of palm and non-palm quadrats. Palm and non-palm quadrats may have equivalent species richness yet these species need not be identical. We also compared the abundance of every species in the plot in palm and non-palm quadrats using a proportion test, with the expected relative abundance set to 85% of individuals in non-palm quadrats (the proportion of all quadrats that were non-palm). Species with significantly higher or lower abundances than expected (P < 0.05) were identified.
ISAR is a statistic to analyse the spatial diversity structure in forest ecosystems and reconciles common species area relationships [67–69] and the individual perspective of point-pattern analysis [63, 70]. The ISAR allows for a subtle assessment of the effect of palms (or adult Myristicaceae) trees on the non-palm sapling diversity in their neighbourhood. If negative interactions of the focal species on saplings dominate, the target species would decrease the proportion of species in its neighborhood (i.e., a diversity repeller). Conversely, if the target species exerts positive interactions to non-palm saplings, the target species would accumulate and maintain an over-representative proportion of sapling diversity in its proximity (i.e., being a diversity accumulator). Finally, if positive and negative interactions are weak or even out, the focal species would behave neutrally.
Similar to the case of the neighbourhood density described above, environmental heterogeneity may cause repeller or accumulator effects in the absence of species interactions if the focal species is predominantly located in areas of low or high sapling diversity. We therefore contrasted the observed ISAR(r) function to that of the same null models as for the analyses of neighbourhood density. In the first null model we distributed individuals of palms randomly over the entire plot. The second test also assumed a random distribution of palms, but the palms were only displaced within a given neighborhood of 30 m to remove in the null model potential interactions with saplings and to guarantee that the displaced palms will be as well located in areas with low sapling density. The neighborhood of 30 m is somewhat larger than the expected range of direct interactions between focal trees. These analyses were conducted in Programita, for more details of ISAR, see Wiegand et al. .
We hypothesize that large palms are a diversity repeller because of the disturbance effects of falling fronds or rooting damage by palm fruit consumers in the vicinity of palms. Only a subset of species would have the traits necessary for resprouting and continued growth following damage from nearby palms.
Dynamic demographic effects of large trees on seedlings
To examine the effects of proximity to large palms and Myristicaceae on seedling survival, for each of the 339 seedling plots located in the FDP we calculated the total number of new recruits of all species between 2003 and 2007. We then calculated the number of each cohort surviving one year. We modelled the proportion of survivors in each plot as a function of location relative to mature and immature large trees with a generalised linear model with binomial errors as above.
Dynamic demographic effects of large trees on saplings
We examined the survival, growth and resprouting of saplings 1-2 cm dbh from the first to the second plot census as a function of whether the sapling was near to (within 10 m) or far from (> 10 m) focal trees (palm or Myristicaceae). First, we modelled the probability of survival to the second census of all saplings alive in the first census. Second, we modelled the probability of each sapling resprouting (new sprouts growing up and out from damaged or broken stems). Third, growth rate for each sapling from the first census to the second census was calculated as growth = (dbh1 dbh0)/(time1 time0) and divided into negative (1) or ≥ zero growth (0). We modelled the probability of a negative growth rate between censuses. In all three cases we used a generalised linear model with binomial errors. We modelled survival (or growth or resprouting, all binary response variables) as a function of proximity to the nearest focal tree (mature and immature palm or male and female Myristicaceae): survival ~ mature palms * immature palms. We coded each sapling as near or far (</> 10 m) from focal trees. We included the interaction term because saplings could be near to a mature palm but far from an immature palm, for example.
We thank two reviewers for their comments that improved the manuscript. The Forest Dynamics Plot of Yasuní National Park has been made possible through the generous support of the Pontifical Catholic University of Ecuador funds of donaciones del impuesto a la renta, the government of Ecuador, the US National Science Foundation, the Andrew W. Mellon Foundation, the Smithsonian Tropical Research Institute, and the University of Aarhus of Denmark. The Yasuní Forest Dynamics Plot is associated to the Center for Tropical Forest Science, a global network of large-scale demographic tree plots.
- Connell J: Diversity in tropical rain forests and coral reefs. Science. 1978, 199: 1302-1310. 10.1126/science.199.4335.1302. [[doi:10.1126/science.199.4335.1302]]View ArticlePubMedGoogle Scholar
- Sousa WP: The role of disturbance in natural communities. Annual Review of Ecology and Systematics. 1984, 15: 353-391. 10.1146/annurev.es.15.110184.002033. [http://arjournals.annualreviews.org/doi/abs/10.1146%2Fannurev.es.15.110184.002033]View ArticleGoogle Scholar
- The Ecology of Natural Disturbance and Patch Dynamics. Edited by: Pickett S, White P. 1987, Academic Press
- Denslow J: Tropical rainforest gaps and tree species diversity. Annual Review of Ecology and Systematics. 1987, 18: 431-451. 10.1146/annurev.es.18.110187.002243. [http://www.jstor.org/stable/2097139]View ArticleGoogle Scholar
- Condit R: Defining and mapping vegetation types in mega-diverse tropical forests. Trends in Ecology & Evolution. 1996, 11: 4-5. [http://www.ncbi.nlm.nih.gov/pubmed/21237740]View ArticleGoogle Scholar
- Hubbell SP, Foster RB, O'Brien S, Harms K, Condit R, Wechsler B, Wright S, de Lao SL: Light-gap disturbances, recruitment limitation, and tree diversity in a neotropical forest. Science. 1999, 283 (5401): 554-557. [http://www.sciencemag.org/content/283/5401/554.full]View ArticlePubMedGoogle Scholar
- Connell J, Slatyer R: Mechanisms of succession in natural communities and their role in community stability and organization. American Naturalist. 1977, 111 (982): 1119-1144. 10.1086/283241. [http://www.jstor.org/stable/2460259]View ArticleGoogle Scholar
- Swaine M, Whitmore T: On the definition of ecological species groups in tropical rain forests. Vegetatio. 1988, 75: 81-86. 10.1007/BF00044629.View ArticleGoogle Scholar
- Welden C, Hewett S, Hubbell SP, Foster R: Sapling survival, growth, and recruitment: relationship to canopy height in a Neotropical forest. Ecology. 1991, 72: 35-50. 10.2307/1938900. [http://www.esajournals.org/doi/abs/10.2307/1938900]View ArticleGoogle Scholar
- Wright SJ, Muller-Landau HC, Condit R, Hubbell SP: Gap-dependent recruitment, realized vital rates, and size distributions of tropical trees. Ecology. 2003, 84 (12): 3174-3185. 10.1890/02-0038. [http://www.esajournals.org/doi/abs/10.1890/02-0038]View ArticleGoogle Scholar
- Bellingham P, Sparrow A: Resprouting as a life history strategy in woody plant communities. Oikos. 2000, 89 (2): 409-416. 10.1034/j.1600-0706.2000.890224.x. [http://www.jstor.org/stable/3547338]View ArticleGoogle Scholar
- Wright SJ: Plant diversity in tropical forests: a review of mechanisms of species coexistence. Oecologia. 2002, 130: 1-14. [http://www.springerlink.com/index/LN6GP1XUKNV00391.pdf]View ArticleGoogle Scholar
- Connell J, Lowman M, Noble I: Subcanopy gaps in temperate and tropical forest. Aust J Ecol. 1997, 22: 163-168. 10.1111/j.1442-9993.1997.tb00655.x.View ArticleGoogle Scholar
- Wirth R, Weber B, Ryel RJ: Spatial and temporal variability of canopy structure in a tropical moist forest. Acta Oecologica. 2001, 22: 235-244. 10.1016/S1146-609X(01)01123-7.View ArticleGoogle Scholar
- Peters H, Pauw A, Silman M, Terborgh J: Falling palm fronds structure Amazonian rainforest sapling communities. Proceedings of the Royal Society of London. Series B: Biological Sciences. 2004, 271 (Suppl 5): S367-S369. [http://rspb.royalsocietypublishing.org/content/271/Suppl_5/S367.short]View ArticlePubMedPubMed CentralGoogle Scholar
- Rüger N, Huth A, Hubbell SP, Condit R: Response of recruitment to light availability across a tropical lowland rain forest community. Journal of Ecology. 2009, 97 (6): 1360-1368. 10.1111/j.1365-2745.2009.01552.x. [http://doi.wiley.com/10.1111/j.1365-2745.2009.01552.x]View ArticleGoogle Scholar
- Rüger N, Berger U, Hubbell SP, Vieilledent G, Condit R: Growth strategies of tropical tree species: disentangling light and size effects. PLoS ONE. 2011, 6 (9): e25330-10.1371/journal.pone.0025330. [http://dx.plos.org/10.1371/journal.pone.0025330]View ArticlePubMedPubMed CentralGoogle Scholar
- Gillman LN, Ogden J: Physical damage by litterfall to canopy tree seedlings in two temperate New Zealand forests. Journal of Vegetation Science. 2001, 12 (5): 671-676. 10.2307/3236907. [http://doi.wiley.com/10.2307/3236907]View ArticleGoogle Scholar
- Gillman L, Wright S, Ogden J: Response of forest tree seedlings to simulated litterfall damage. Plant Ecology. 2003, 169: 53-60. 10.1023/A:1026288306932. [http://www.springerlink.com/index/R528324220G50583.pdf]View ArticleGoogle Scholar
- Drake D, Pratt L: Seedling mortality in Hawaiian rain forest: the role of small-scale physical disturbance. Biotropica. 2001, 33 (2): 319-323. [http://onlinelibrary.wiley.com/doi/10.1111/j.1744-7429.2001.tb00182.x/abstract]View ArticleGoogle Scholar
- Alvarez-Clare S, Kitajima K: Susceptibility of tree seedlings to biotic and abiotic hazards in the understory of a moist tropical forest in Panama. Biotropica. 2009, 41: 47-56. 10.1111/j.1744-7429.2008.00442.x. [http://doi.wiley.com/10.1111/j.1744-7429.2008.00442.x]View ArticleGoogle Scholar
- Aide T: Limbfalls: a major cause of sapling mortality for tropical forest plants. Biotropica. 1987, 19: 284285-View ArticleGoogle Scholar
- Rich PM, Holbrook NM, Luttinger N: Leaf development and crown geometry of two Iriarteoid palms. Am J Bot. 1995, 82: 328-336. 10.2307/2445578.View ArticleGoogle Scholar
- Valencia R, Balslev H, Palacios W, Neill D, Josse C, Tirado M: Diversity and family composition of trees in different regions of Ecuador: a sample of 18 one-ha plots. Forest Biodiversity in North, Central and South America, and the Caribbean: Research and Monitoring. Man and Biosphere Series 21. Edited by: Dallmeier F, Komiskey J. 1998, The Parthenon Publishing Group, Paris, 569-584.Google Scholar
- Pitman NCA, Terborgh JW, Silman M, Nuñez VP, Neill DA, Cerón CE, Palacios WA, Aulestia M: Dominance and distribution of tree species in upper Amazonian terra firme forests. Ecology. 2001, 82 (8): 2101-2117. 10.1890/0012-9658(2001)082[2101:DADOTS]2.0.CO;2. [http://www.esajournals.org/doi/full/10.1890/0012-9658%282001%29082%5B2101%3ADADOTS%5D2.0.CO%3B2]View ArticleGoogle Scholar
- Pitman NCA, Terborgh JW, Silman MR, Nuñez VP, Neill DA, Ceron CE, Palacios WA, Aulestia M: A comparison of tree species diversity in two upper amazonian forests. Ecology. 2002, 83 (11): 3210-10.1890/0012-9658(2002)083[3210:ACOTSD]2.0.CO;2. [http://www.jstor.org/stable/3071854?origin=crossref]View ArticleGoogle Scholar
- Vormisto J, Svenning JC, Hall P: Diversity and dominance in palm (Arecaceae) communities in terra firme forests in the western Amazon basin. Journal of Ecology. 2004, 92 (4): 577-588. 10.1111/j.0022-0477.2004.00904.x.View ArticleGoogle Scholar
- Macía MJ, Svenning JC: Oligarchic dominance in western Amazonian plant communities. Journal of Tropical Ecology. 2005, 21 (06): 613-10.1017/S0266467405002579. [http://journals.cambridge.org/action/displayAbstract?fromPage=online&aid=345021]View ArticleGoogle Scholar
- Montufar R, Pintaud JC: Variation in species composition, abundance and microhabitat preferences among western Amazonian terra firme palm communities. Botanical Journal of the Linnean Society. 2006, 151: 127-140. 10.1111/j.1095-8339.2006.00528.x. [http://doi.wiley.com/10.1111/j.1095-8339.2006.00528.x]View ArticleGoogle Scholar
- Fragoso JMV: Tapir-generated seed shadows: scale-dependent patchiness in the Amazon rain forest. The Journal of Ecology. 1997, 85 (4): 519-10.2307/2960574. [http://www.jstor.org/stable/2960574?origin=crossref]View ArticleGoogle Scholar
- Ickes K, Dewalt SJ, Thomas SC: Resprouting of woody saplings following stem snap by wild pigs in a Malaysian rain forest. J Ecol. 2003, 91: 222-233. 10.1046/j.1365-2745.2003.00767.x.View ArticleGoogle Scholar
- Silman MR, Terborgh J, Kiltie R: Population regulation of a dominant rain forest tree by a major seed predator. Ecology. 2003, 84 (2): 431-438. 10.1890/0012-9658(2003)084[0431:PROADR]2.0.CO;2. [http://www.esajournals.org/doi/full/10.1890/0012-9658%282003%29084%5B0431%3APROADR%5D2.0.CO%3B2]View ArticleGoogle Scholar
- Beck H: Seed predation and dispersal by peccaries throughout the Neotropics and its consequence: a review and synthesis. Seed fate: predation, dispersal and seedling establishment. Edited by: Forget PM, Lambert J, Hulme P, Vander Wall S. 2005, CABI, 77-15.View ArticleGoogle Scholar
- Wyatt J, Silman MR: Distance-dependence in two Amazonian palms: effects of spatial and temporal variation in seed predator communities. Oecologia. 2004, 140: 26-35. 10.1007/s00442-004-1554-y. [http://www.springerlink.com/index/brvwa9wl18p0vh1j.pdf]View ArticlePubMedGoogle Scholar
- Fragoso JMV, Silvius KM, Correa JA: Long-distance seed dispersal by tapirs increases seed survival and aggregates tropical trees. Ecology. 2003, 84 (8): 1998-2006. 10.1890/01-0621. [http://www.esajournals.org/doi/abs/10.1890/01-0621]View ArticleGoogle Scholar
- Beck H: Synergistic impacts of ungulates and falling palm fronds on saplings in the Amazon. Journal of Tropical Ecology. 2007, 23 (05): 599-602. 10.1017/S0266467407004294. [http://journals.cambridge.org/action/displayAbstract?fromPage=online&aid=1304332]View ArticleGoogle Scholar
- Howe H: Annual variation in a neotropical seed-dispersal system. Tropical rain forests: ecology and management. Edited by: Sutton S, Whitmore T, Chadwick A. 1983, Blackwell, Oxford, UK, 211-22.Google Scholar
- Roosmalen M, Bardales P, Garcia O: Frutos da floresta Amazonica. Parte I: Myristicaceae. Acta Amazonica. 1996, 25: 209264-Google Scholar
- Beck H: A review of peccarypalm interactions and their ecological ramifications across the Neotropics. Journal of Mammalogy. 2006, 87 (3): 519-530. 10.1644/05-MAMM-A-174R1.1. [http://www.bioone.org/doi/abs/10.1644/05-MAMM-A-174R1.1]View ArticleGoogle Scholar
- Kiltie RA: Distribution of palm fruits on a rain forest floor: why White-Lipped Peccaries forage near objects. Biotropica. 1981, 13 (2): 141-145. 10.2307/2387716.View ArticleGoogle Scholar
- Terborgh J, Lopez L, Nuñez P, Rao M, Shahabuddin G, Orihuela G, Riveros M, Ascanio R, Adler GH, Lambert TD, Balbas L: Ecological meltdown in predator-free forest fragments. Science (New York, N.Y.). 2001, 294 (5548): 1923-6. 10.1126/science.1064397. [http://www.ncbi.nlm.nih.gov/pubmed/11729317]View ArticleGoogle Scholar
- Gilbert G, Harms K, Hamill D, Hubbell SP: Effects of seedling size, El Niño drought, seedling density, and distance to nearest conspecific adult on 6-year survival of Ocotea whitei seedlings in Panamá. Oecologia. 2001, 127 (4): 509-516. 10.1007/s004420000616. [http://www.springerlink.com/index/0CY7WTKN74RAVQEH.pdf]View ArticleGoogle Scholar
- Grubb PJ: The maintenance of species-richness in plant communities: the importance of the regeneration niche. Biological Reviews. 1977, 52: 107-145. [http://doi.wiley.com/10.1111/j.1469-185X.1977.tb01347.x]View ArticleGoogle Scholar
- Bullock S: Demography of an undergrowth palm in littoral Cameroon. Biotropica. 1980, 12 (4): 247-255. 10.2307/2387694. [http://www.jstor.org/stable/2387694]View ArticleGoogle Scholar
- Silva Matos DM, Freckleton RP, Watkinson AR: The role of density dependence in the population dynamics of a tropical palm. Ecology. 1999, 80 (8): 2635-2650. [http://www.esajournals.org/doi/abs/10.1890/0012-9658%281999%29080%5B2635%3ATRODDI%5D2.0.CO%3B2]View ArticleGoogle Scholar
- Vandermeer J: Notes on density dependence of Welfia georgii Wendl. ex Burrett (Palmae), a lowland rainforest species in Costa Rica. Brensia. 1977, 10/11: 9-15.Google Scholar
- Dransfield J: Growth forms of rain forest plants. Tropical Trees as Living Systems. Edited by: Tomlinson P, Zimmerman M. 1978, Cambridge University Press, Cambridge, 247-268.Google Scholar
- Sarukhan J: Studies on the demography of tropical trees. Tropical Trees as Living Systems. Edited by: Tomlinson P, Zimmerman M. 1978, Cambridge University Press, Cambridge, 247-268.Google Scholar
- Álvarez Loayza P, White JF, Torres MS, Balslev H, Kristiansen T, Svenning JC, Gil N: Light converts endosymbiotic fungus to pathogen, influencing seedling survival and niche-space filling of a common tropical tree, Iriartea deltoidea. PloS one. 2011, 6: e16386-10.1371/journal.pone.0016386. [http://www.plosone.org/article/info%253Adoi%252F10.1371%252Fjournal.pone.0016386]View ArticlePubMedPubMed CentralGoogle Scholar
- Valencia R, Foster RB, Villa G, Condit R, Svenning J, Hernandez C, Romoleroux K, Losos E, Magard E, Balslev H: Tree species distributions and local habitat variation in the Amazon: large forest plot in eastern Ecuador. J Ecol. 2004, 92: 214-229. 10.1111/j.0022-0477.2004.00876.x.View ArticleGoogle Scholar
- Bass M, Finer M, Jenkins C, Kreft H, Cisneros-Heredia D, McCracken S, Pitman NCA, English P, Swing K, Villa G, DiFiore A, Voigt CC, Kunz TH: Global conservation significance of Ecuador's Yasuní National Park. PloS one. 2010, 5: e8767-10.1371/journal.pone.0008767. [http://www.plosone.org/article/info:doi/10.1371/journal.pone.0008767]View ArticlePubMedPubMed CentralGoogle Scholar
- Valencia R, Condit R, Muller-Landau HC, Hernandez C, Navarrete H: Dissecting biomass dynamics in a large Amazonian forest plot. Journal of Tropical Ecology. 2009, 25: 473-482. 10.1017/S0266467409990095. [http://journals.cambridge.org/action/displayAbstract?fromPage=online&aid=5991576]View ArticleGoogle Scholar
- Homeier J, Breckle SW, Dalitz H, Leyers C, Ortiz R: Demography, spatial distribution, and growth of three aroborescent palm species in a tropical premontane rain forest in Costa Rica. Ecotropica. 2002, 8: 239-247.Google Scholar
- Queenborough SA, Burslem DFRP, Garwood NC, Valencia R: Determinants of biased sex ratios and inter-sex costs of reproduction in dioecious tropical forest trees. American journal of botany. 2007, 94: 67-78. 10.3732/ajb.94.1.67. [http://www.ncbi.nlm.nih.gov/pubmed/21642209]View ArticlePubMedGoogle Scholar
- Queenborough SA, Burslem DFRP, Garwood NC, Valencia R: Habitat niche partitioning by 16 species of Myristicaceae in Amazonian Ecuador. Plant Ecology. 2007, 192 (2): 193-207. 10.1007/s11258-007-9328-3. [http://www.springerlink.com/index/10.1007/s11258-007-9328-3]View ArticleGoogle Scholar
- Queenborough SA, Burslem DFRP, Garwood NC, Valencia R: Neighborhood and community interactions determine the spatial pattern of tropical tree seedling survival. Ecology. 2007, 88 (9): 2248-58. 10.1890/06-0737.1. [http://www.ncbi.nlm.nih.gov/pubmed/17918403]View ArticlePubMedGoogle Scholar
- Queenborough SA, Burslem DFRP, Garwood NC, Valencia R: Taxonomic scale-dependence of habitat niche partitioning and biotic neighbourhood on survival of tropical tree seedlings. Proceedings of the Royal Society B: Biological Sciences. 2009, 276 (1676): 4197-4205. 10.1098/rspb.2009.0921. [http://rspb.royalsocietypublishing.org/content/276/1676/4197.short]View ArticlePubMedPubMed CentralGoogle Scholar
- Queenborough SA, Mazer SJ, Vamosi SM, Garwood NC, Valencia R, Freckleton RP: Seed mass, abundance and breeding system among tropical forest species: do dioecious species exhibit compensatory reproduction or abundances?. Journal of Ecology. 2009, 97 (3): 555-566. 10.1111/j.1365-2745.2009.01485.x. [http://onlinelibrary.wiley.com/doi/10.1111/j.1365-2745.2009.01485.x/full]View ArticleGoogle Scholar
- Metz MR, Sousa WP, Valencia R: Widespread density-dependent seedling mortality promotes species coexistence in a highly diverse Amazonian rain forest. Ecology. 2010, 91 (12): 3675-85. 10.1890/08-2323.1. [http://www.ncbi.nlm.nih.gov/pubmed/21302838]View ArticlePubMedGoogle Scholar
- R Development Core Team: R: A Language and Environment for Statistical Computing. 2011, R Foundation for Statistical Computing, Vienna, Austria, [ISBN 3-900051-07-0], [http://www.R-project.org/]Google Scholar
- Wiegand T, A Moloney K: Rings, circles, and null-models for point pattern analysis in ecology. Oikos. 2004, 104 (2): 209-229. 10.1111/j.0030-1299.2004.12497.x. [http://doi.wiley.com/10.1111/j.0030-1299.2004.12497.x]View ArticleGoogle Scholar
- Perry G, Miller B, Enright N: A comparison of methods for the statistical analysis of spatial point patterns in plant ecology. Plant Ecol. 2006, 187: 59-82. 10.1007/s11258-006-9133-4.View ArticleGoogle Scholar
- Wiegand T, Gunatilleke S, Gunatilleke N: Species associations in a heterogeneous Sri Lankan dipterocarp forest. The American Naturalist. 2007, 170 (4): E77-95. 10.1086/521240. [http://www.ncbi.nlm.nih.gov/pubmed/17891727]View ArticlePubMedGoogle Scholar
- Wiegand T, Gunatilleke CVS, Gunatilleke IaUN, Huth A: How individual species structure diversity in tropical forests. Proceedings of the National Academy of Sciences of the United States of America. 2007, 104 (48): 19029-33. 10.1073/pnas.0705621104. [http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=2141902&tool=pmcentrez&rendertype=abstract]View ArticlePubMedPubMed CentralGoogle Scholar
- Wang X, Wiegand T, Hao Z, Li B, Ye J, Lin F: Species associations in an old-growth temperate forest in north-eastern China. Journal of Ecology. 2010, 98 (3): 674-686. 10.1111/j.1365-2745.2010.01644.x. [http://doi.wiley.com/10.1111/j.1365-2745.2010.01644.x]View ArticleGoogle Scholar
- Oksanen J, Blanchet FG, Kindt R, Legendre P, O' Hara RB, Simpson GL, Solymos P, Stevens MHH, Wagner H: vegan: Community Ecology Package. 2011, [R package version 1.17-11], [http://CRAN.R-project.org/package=vegan]Google Scholar
- He F, Legendre P: Species diversity patterns derived from species-area models. Ecology. 2002, 83: 1185-11988.Google Scholar
- Plotkin J, Potts M, Yu D, Bunyavejchewin S, Condit R, Foster R, Hubbell S, LaFrankie J, Manokaran N, Seng L, Sukumar R, Nowak MA, Ashton PS: Predicting species diversity in tropical forests. Proceedings of the National Academy of Sciences. 2000, 97 (20): 10850-10854. [http://www.pnas.org/content/97/20/10850.short]View ArticleGoogle Scholar
- Plotkin J, Potts M, Leslie N, Manokaran N, LaFrankie J, Ashton P: Species-area curves, spatial aggregation, and habitat specialization in tropical forests. J Theor Biol. 2000, 207: 81-99. 10.1006/jtbi.2000.2158.View ArticlePubMedGoogle Scholar
- Illian J, Penttinen A, Stoyan H, Stoyan D: Statistical Analysis and Modelling of Spatial Point Patterns. 2008, Chichester, UK: WileyGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.