Latitudinal patterns in stabilizing density dependence of forest communities

Numerous studies have shown reduced performance in plants that are surrounded by neighbours of the same species1,2, a phenomenon known as conspecific negative density dependence (CNDD)3. A long-held ecological hypothesis posits that CNDD is more pronounced in tropical than in temperate forests4,5, which increases community stabilization, species coexistence and the diversity of local tree species6,7. Previous analyses supporting such a latitudinal gradient in CNDD8,9 have suffered from methodological limitations related to the use of static data10-12. Here we present a comprehensive assessment of latitudinal CNDD patterns using dynamic mortality data to estimate species-site-specific CNDD across 23 sites. Averaged across species, we found that stabilizing CNDD was present at all except one site, but that average stabilizing CNDD was not stronger toward the tropics. However, in tropical tree communities, rare and intermediate abundant species experienced stronger stabilizing CNDD than did common species. This pattern was absent in temperate forests, which suggests that CNDD influences species abundances more strongly in tropical forests than it does in temperate ones13. We also found that interspecific variation in CNDD, which might attenuate its stabilizing effect on species diversity14,15, was high but not significantly different across latitudes. Although the consequences of these patterns for latitudinal diversity gradients are difficult to evaluate, we speculate that a more effective regulation of population abundances could translate into greater stabilization of tropical tree communities and thus contribute to the high local diversity of tropical forests.

Tropical tree ectomycorrhiza are distributed independently of soil nutrients

Mycorrhizae, a form of plant-fungal symbioses, mediate vegetation impacts on ecosystem functioning. Climatic effects on decomposition and soil quality are suggested to drive mycorrhizal distributions, with arbuscular mycorrhizal plants prevailing in low-latitude/high-soil-quality areas and ectomycorrhizal (EcM) plants in high-latitude/low-soil-quality areas. However, these generalizations, based on coarse-resolution data, obscure finer-scale variations and result in high uncertainties in the predicted distributions of mycorrhizal types and their drivers. Using data from 31 lowland tropical forests, both at a coarse scale (mean-plot-level data) and fine scale (20 × 20 metres from a subset of 16 sites), we demonstrate that the distribution and abundance of EcM-associated trees are independent of soil quality. Resource exchange differences among mycorrhizal partners, stemming from diverse evolutionary origins of mycorrhizal fungi, may decouple soil fertility from the advantage provided by mycorrhizal associations. Additionally, distinct historical biogeographies and diversification patterns have led to differences in forest composition and nutrient-acquisition strategies across three major tropical regions. Notably, Africa and Asia's lowland tropical forests have abundant EcM trees, whereas they are relatively scarce in lowland neotropical forests. A greater understanding of the functional biology of mycorrhizal symbiosis is required, especially in the lowland tropics, to overcome biases from assuming similarity to temperate and boreal regions.

Mycorrhizal feedbacks influence global forest structure and diversity

One mechanism proposed to explain high species diversity in tropical systems is strong negative conspecific density dependence (CDD), which reduces recruitment of juveniles in proximity to conspecific adult plants. Although evidence shows that plant-specific soil pathogens can drive negative CDD, trees also form key mutualisms with mycorrhizal fungi, which may counteract these effects. Across 43 large-scale forest plots worldwide, we tested whether ectomycorrhizal tree species exhibit weaker negative CDD than arbuscular mycorrhizal tree species. We further tested for conmycorrhizal density dependence (CMDD) to test for benefit from shared mutualists. We found that the strength of CDD varies systematically with mycorrhizal type, with ectomycorrhizal tree species exhibiting higher sapling densities with increasing adult densities than arbuscular mycorrhizal tree species. Moreover, we found evidence of positive CMDD for tree species of both mycorrhizal types. Collectively, these findings indicate that mycorrhizal interactions likely play a foundational role in global forest diversity patterns and structure.

Long-term seedling and small sapling census data from the Barro Colorado Island 50 ha Forest Dynamics Plot, Panama

Tropical forests are well known for their high woody plant diversity. Processes occurring at early life stages are thought to play a critical role in maintaining this high diversity and shaping the composition of tropical tree communities. To evaluate hypothesized mechanisms promoting tropical tree species coexistence and influencing composition, we initiated a census of woody seedlings and small saplings in the permanent 50 ha Forest Dynamics Plot (FDP) on Barro Colorado Island (BCI), Panama. Situated in old-growth, lowland tropical moist forest, the BCI FDP was originally established in 1980 to monitor trees and shrubs ≥1 cm diameter at 1.3 m above ground (dbh) at ca. 5-year intervals. However, critical data on the dynamics occurring at earlier life stages were initially lacking. Therefore, in 2001 we established a 1-m2 seedling plot in the center of every 5 × 5 m section of the BCI FDP. All freestanding woody individuals ≥20 cm tall and <1 cm dbh (hereafter referred to as seedlings) were tagged, mapped, measured, and identified to species in 19,313 1-m2 seedling plots. Because seedling dynamics are rapid, we censused these seedling plots every 1-2 years. Here, we present data from the 14 censuses of these seedling plots conducted between the initial census in 2001 to the most recent census, in 2018. This data set includes nearly 1 M observations of ~185,000 individuals of >400 tree, shrub, and liana species. These data will permit spatially-explicit analyses of seedling distributions, recruitment, growth, and survival for hundreds of woody plant species. In addition, the data presented here can be linked to openly-available, long-term data on the dynamics of trees and shrubs ≥1 cm dbh in the BCI FDP, as well as existing data sets from the site on climate, canopy structure, phylogenetic relatedness, functional traits, soil nutrients, and topography. This data set can be freely used for non-commercial purposes; we request that users of these data cite this data paper in all publications resulting from the use of this data set.

RecruitNet: A global database of plant recruitment networks

Plant recruitment interactions (i.e., what recruits under what) shape the composition, diversity, and structure of plant communities. Despite the huge body of knowledge on the mechanisms underlying recruitment interactions among species, we still know little about the structure of the recruitment networks emerging in ecological communities. Modeling and analyzing the community-level structure of plant recruitment interactions as a complex network can provide relevant information on ecological and evolutionary processes acting both at the species and ecosystem levels. We report a data set containing 143 plant recruitment networks in 23 countries across five continents, including temperate and tropical ecosystems. Each network identifies the species under which another species recruits. All networks report the number of recruits (i.e., individuals) per species. The data set includes >850,000 recruiting individuals involved in 118,411 paired interactions among 3318 vascular plant species across the globe. The cover of canopy species and open ground is also provided. Three sampling protocols were used: (1) The Recruitment Network (RN) protocol (106 networks) focuses on interactions among established plants ("canopy species") and plants in their early stages of recruitment ("recruit species"). A series of plots was delimited within a locality, and all the individuals recruiting and their canopy species were identified; (2) The paired Canopy-Open (pCO) protocol (26 networks) consists in locating a potential canopy plant and identifying recruiting individuals under the canopy and in a nearby open space of the same area; (3) The Georeferenced plot (GP) protocol (11 networks) consists in using information from georeferenced individual plants in large plots to infer canopy-recruit interactions. Some networks incorporate data for both herbs and woody species, whereas others focus exclusively on woody species. The location of each study site, geographical coordinates, country, locality, responsible author, sampling dates, sampling method, and life habits of both canopy and recruit species are provided. This database will allow researchers to test ecological, biogeographical, and evolutionary hypotheses related to plant recruitment interactions. There are no copyright restrictions on the data set; please cite this data paper when using these data in publications.

Distribution of biomass dynamics in relation to tree size in forests across the world

Tree size shapes forest carbon dynamics and determines how trees interact with their environment, including a changing climate. Here, we conduct the first global analysis of among-site differences in how aboveground biomass stocks and fluxes are distributed with tree size. We analyzed repeat tree censuses from 25 large-scale (4-52 ha) forest plots spanning a broad climatic range over five continents to characterize how aboveground biomass, woody productivity, and woody mortality vary with tree diameter. We examined how the median, dispersion, and skewness of these size-related distributions vary with mean annual temperature and precipitation. In warmer forests, aboveground biomass, woody productivity, and woody mortality were more broadly distributed with respect to tree size. In warmer and wetter forests, aboveground biomass and woody productivity were more right skewed, with a long tail towards large trees. Small trees (1-10 cm diameter) contributed more to productivity and mortality than to biomass, highlighting the importance of including these trees in analyses of forest dynamics. Our findings provide an improved characterization of climate-driven forest differences in the size structure of aboveground biomass and dynamics of that biomass, as well as refined benchmarks for capturing climate influences in vegetation demographic models.

Demographic composition, not demographic diversity, predicts biomass and turnover across temperate and tropical forests

The growth and survival of individual trees determine the physical structure of a forest with important consequences for forest function. However, given the diversity of tree species and forest biomes, quantifying the multitude of demographic strategies within and across forests and the way that they translate into forest structure and function remains a significant challenge. Here, we quantify the demographic rates of 1961 tree species from temperate and tropical forests and evaluate how demographic diversity (DD) and demographic composition (DC) differ across forests, and how these differences in demography relate to species richness, aboveground biomass (AGB), and carbon residence time. We find wide variation in DD and DC across forest plots, patterns that are not explained by species richness or climate variables alone. There is no evidence that DD has an effect on either AGB or carbon residence time. Rather, the DC of forests, specifically the relative abundance of large statured species, predicted both biomass and carbon residence time. Our results demonstrate the distinct DCs of globally distributed forests, reflecting biogeography, recent history, and current plot conditions. Linking the DC of forests to resilience or vulnerability to climate change, will improve the precision and accuracy of predictions of future forest composition, structure, and function.

The relative abundance of languages: Neutral and non-neutral dynamics

Credible estimates suggest that a large number of the nearly 7000 languages in the world could go extinct this century, a prospect with profound cultural, socioeconomic, and political ramifications. Despite its importance, we still have little predictive theory for language dynamics and richness. Critical to the language extinction problem, however, is to understand the dynamics of the number of speakers of languages, the dynamics of language abundance distributions (LADs). Many regional LADs are very similar to the bell-shaped distributions of relative species abundance predicted by neutral theory in ecology. Using the tenets of neutral theory, here we show that LADs can be understood as an equilibrium or disequilibrium between stochastic rates of origination and extinction of languages. However, neutral theory does not fit some regional LADs, which can be explained if the number of speakers has grown systematically faster in some languages than others, due to cultural factors and other non-neutral processes. Only the LADs of Australia and the United States, deviate from a bell-shaped pattern. These deviations are due to the documented higher, non-equilibrium extinction rates of low-abundance languages in these countries.

Local canopy disturbance as an explanation for long-term increases in liana abundance

Canopy disturbance explains liana abundance and distribution within tropical forests and thus may also explain the widespread pattern of increasing liana abundance; however, this hypothesis remains untested. We used a 10-year study (2007-2017) of 117,100 rooted lianas in an old-growth Panamanian forest to test whether local canopy disturbance explains increasing liana abundance. We found that liana density increased 29.2% and basal area 12.5%. The vast majority of these increases were associated with clonal stem proliferation following canopy disturbance, particularly in liana-dense, low-canopy gaps, which had far greater liana increases than did undisturbed forest. Lianas may be ecological niche constructors, arresting tree regeneration in gaps and thus creating a high-light environment that favours sustained liana proliferation. Our findings demonstrate that liana abundance is increasing rapidly and their ability to proliferate via copious clonal stem production in canopy gaps explains much of their increase in this and possibly other tropical forests.

Arbuscular mycorrhizal trees influence the latitudinal beta-diversity gradient of tree communities in forests worldwide

Arbuscular mycorrhizal (AM) and ectomycorrhizal (EcM) associations are critical for host-tree performance. However, how mycorrhizal associations correlate with the latitudinal tree beta-diversity remains untested. Using a global dataset of 45 forest plots representing 2,804,270 trees across 3840 species, we test how AM and EcM trees contribute to total beta-diversity and its components (turnover and nestedness) of all trees. We find AM rather than EcM trees predominantly contribute to decreasing total beta-diversity and turnover and increasing nestedness with increasing latitude, probably because wide distributions of EcM trees do not generate strong compositional differences among localities. Environmental variables, especially temperature and precipitation, are strongly correlated with beta-diversity patterns for both AM trees and all trees rather than EcM trees. Results support our hypotheses that latitudinal beta-diversity patterns and environmental effects on these patterns are highly dependent on mycorrhizal types. Our findings highlight the importance of AM-dominated forests for conserving global forest biodiversity.

Genome assemblies for two Neotropical trees: Jacaranda copaia and Handroanthus guayacan

The lack of genomic resources for tropical canopy trees is impeding several research avenues in tropical forest biology. We present genome assemblies for two Neotropical hardwood species, Jacaranda copaia and Handroanthus (formerly Tabebuia) guayacan, that are model systems for research on tropical tree demography and flowering phenology. For each species, we combined Illumina short-read data with in vitro proximity-ligation (Chicago) libraries to generate an assembly. For Jacaranda copaia, we obtained 104X physical coverage and produced an assembly with N50/N90 scaffold lengths of 1.020/0.277 Mbp. For H. guayacan, we obtained 129X coverage and produced an assembly with N50/N90 scaffold lengths of 0.795/0.165 Mbp. J. copaia and H. guayacan assemblies contained 95.8% and 87.9% of benchmarking orthologs, although they constituted only 77.1% and 66.7% of the estimated genome sizes of 799 and 512 Mbp, respectively. These differences were potentially due to high repetitive sequence content (>59.31% and 45.59%) and high heterozygosity (0.5% and 0.8%) in each species. Finally, we compared each new assembly to a previously sequenced genome for Handroanthus impetiginosus using whole-genome alignment. This analysis indicated extensive gene duplication in H. impetiginosus since its divergence from H. guayacan.

The interspecific growth-mortality trade-off is not a general framework for tropical forest community structure

Resource allocation within trees is a zero-sum game. Unavoidable trade-offs dictate that allocation to growth-promoting functions curtails other functions, generating a gradient of investment in growth versus survival along which tree species align, known as the interspecific growth-mortality trade-off. This paradigm is widely accepted but not well established. Using demographic data for 1,111 tree species across ten tropical forests, we tested the generality of the growth-mortality trade-off and evaluated its underlying drivers using two species-specific parameters describing resource allocation strategies: tolerance of resource limitation and responsiveness of allocation to resource access. Globally, a canonical growth-mortality trade-off emerged, but the trade-off was strongly observed only in less disturbance-prone forests, which contained diverse resource allocation strategies. Only half of disturbance-prone forests, which lacked tolerant species, exhibited the trade-off. Supported by a theoretical model, our findings raise questions about whether the growth-mortality trade-off is a universally applicable organizing framework for understanding tropical forest community structure.

Demographic trade-offs predict tropical forest dynamics

Understanding tropical forest dynamics and planning for their sustainable management require efficient, yet accurate, predictions of the joint dynamics of hundreds of tree species. With increasing information on tropical tree life histories, our predictive understanding is no longer limited by species data but by the ability of existing models to make use of it. Using a demographic forest model, we show that the basal area and compositional changes during forest succession in a neotropical forest can be accurately predicted by representing tropical tree diversity (hundreds of species) with only five functional groups spanning two essential trade-offs-the growth-survival and stature-recruitment trade-offs. This data-driven modeling framework substantially improves our ability to predict consequences of anthropogenic impacts on tropical forests.

Testing for changes in biomass dynamics in large-scale forest datasets

Tropical forest responses to climate and atmospheric change are critical to the future of the global carbon budget. Recent studies have reported increases in estimated above-ground biomass (EAGB) stocks, productivity, and mortality in old-growth tropical forests. These increases could reflect a shift in forest functioning due to global change and/or long-lasting recovery from past disturbance. We introduce a novel approach to disentangle the relative contributions of these mechanisms by decomposing changes in whole-plot biomass fluxes into contributions from changes in the distribution of gap-successional stages and changes in fluxes for a given stage. Using 30 years of forest dynamic data at Barro Colorado Island, Panama, we investigated temporal variation in EAGB fluxes as a function of initial EAGB (EAGB_i ) in 10 × 10 m quadrats. Productivity and mortality fluxes both increased strongly with initial quadrat EAGB. The distribution of EAGB (and thus EAGB_i ) across quadrats hardly varied over 30 years (and seven censuses). EAGB fluxes as a function of EAGB_i varied largely and significantly among census intervals, with notably higher productivity in 1985-1990 associated with recovery from the 1982-1983 El Niño event. Variation in whole-plot fluxes among census intervals was explained overwhelmingly by variation in fluxes as a function of EAGB_i , with essentially no contribution from changes in EAGB_i distributions. The high observed temporal variation in productivity and mortality suggests that this forest is very sensitive to climate variability. There was no consistent long-term trend in productivity, mortality, or biomass in this forest over 30 years, although the temporal variability in productivity and mortality was so strong that it could well mask a substantial trend. Accurate prediction of future tropical forest carbon budgets will require accounting for disturbance-recovery dynamics and understanding temporal variability in productivity and mortality.

Lightning is a major cause of large tree mortality in a lowland neotropical forest

The mortality rates of large trees are critical to determining carbon stocks in tropical forests, but the mechanisms of tropical tree mortality remain poorly understood. Lightning strikes thousands of tropical trees every day, but is commonly assumed to be a minor agent of tree mortality in most tropical forests. We use the first systematic quantification of lightning-caused mortality to show that lightning is a major cause of death for the largest trees in an old-growth lowland forest in Panama. A novel lightning strike location system together with field surveys of strike sites revealed that, on average, each strike directly kills 3.5 trees (> 10 cm diameter) and damages 11.4 more. Given lightning frequency data from the Earth Networks Total Lightning Network and historical total tree mortality rates for this site, we conclude that lightning accounts for 40.5% of the mortality of large trees (> 60 cm diameter) in the short term and probably contributes to an additional 9.0% of large tree deaths over the long term. Any changes in cloud-to-ground lightning frequency due to climatic change will alter tree mortality rates; projected 25-50% increases in lightning frequency would increase large tree mortality rates in this forest by 9-18%. The results of this study indicate that lightning plays a critical and previously underestimated role in tropical forest dynamics and carbon cycling.

Pervasive decreases in living vegetation carbon turnover time across forest climate zones

Forests play a major role in the global carbon cycle. Previous studies on the capacity of forests to sequester atmospheric CO2 have mostly focused on carbon uptake, but the roles of carbon turnover time and its spatiotemporal changes remain poorly understood. Here, we used long-term inventory data (1955 to 2018) from 695 mature forest plots to quantify temporal trends in living vegetation carbon turnover time across tropical, temperate, and cold climate zones, and compared plot data to 8 Earth system models (ESMs). Long-term plots consistently showed decreases in living vegetation carbon turnover time, likely driven by increased tree mortality across all major climate zones. Changes in living vegetation carbon turnover time were negatively correlated with CO2 enrichment in both forest plot data and ESM simulations. However, plot-based correlations between living vegetation carbon turnover time and climate drivers such as precipitation and temperature diverged from those of ESM simulations. Our analyses suggest that forest carbon sinks are likely to be constrained by a decrease in living vegetation carbon turnover time, and accurate projections of forest carbon sink dynamics will require an improved representation of tree mortality processes and their sensitivity to climate in ESMs.

Seed-to-seedling transitions exhibit distance-dependent mortality but no strong spacing effects in a Neotropical forest

Patterns of seed dispersal and seed mortality influence the spatial structure of plant communities and the local coexistence of competing species. Most seeds are dispersed in proximity to the parent tree, where mortality is also expected to be the highest, because of competition with siblings or the attraction of natural enemies. Whereas distance-dependent mortality in the seed-to-seedling transition was often observed in tropical forests, few studies have attempted to estimate the shape of the survival-distance curves, which determines whether the peak of seedling establishment occurs away from the parent tree (Janzen-Connell pattern) or if the peak attenuates but remains at the parent location (Hubbell pattern). In this study, we inferred the probability density of seed dispersal and two stages of seedling establishment (new recruits, and seedlings 20 cm or taller) with distance for 24 tree species present in the 50-ha Forest Dynamics Plot of Barro Colorado Island, Panama. Using data from seed traps, seedling survey quadrats, and tree-census records spanning the 1988-2014 period, we fit hierarchical Bayesian models including parameters for tree fecundity, the shape of the dispersal kernel, and overdispersion of seed or seedling counts. We combined predictions from multiple dispersal kernels to obtain more robust inferences. We find that Hubbell patterns are the most common and Janzen-Connell patterns are very rare among those species; that distance-dependent mortality may be stronger in the seed stage, in the early recruit stage, or comparable in both; and that species with larger seeds experience less overall mortality and less distance-dependent mortality. Finally, we describe how this modeling approach could be extended at a community scale to include less abundant species.

Adapterama II: universal amplicon sequencing on Illumina platforms (TaggiMatrix)

Next-generation sequencing (NGS) of amplicons is used in a wide variety of contexts. In many cases, NGS amplicon sequencing remains overly expensive and inflexible, with library preparation strategies relying upon the fusion of locus-specific primers to full-length adapter sequences with a single identifying sequence or ligating adapters onto PCR products. In Adapterama I, we presented universal stubs and primers to produce thousands of unique index combinations and a modifiable system for incorporating them into Illumina libraries. Here, we describe multiple ways to use the Adapterama system and other approaches for amplicon sequencing on Illumina instruments. In the variant we use most frequently for large-scale projects, we fuse partial adapter sequences (TruSeq or Nextera) onto the 5' end of locus-specific PCR primers with variable-length tag sequences between the adapter and locus-specific sequences. These fusion primers can be used combinatorially to amplify samples within a 96-well plate (8 forward primers + 12 reverse primers yield 8 × 12 = 96 combinations), and the resulting amplicons can be pooled. The initial PCR products then serve as template for a second round of PCR with dual-indexed iTru or iNext primers (also used combinatorially) to make full-length libraries. The resulting quadruple-indexed amplicons have diversity at most base positions and can be pooled with any standard Illumina library for sequencing. The number of sequencing reads from the amplicon pools can be adjusted, facilitating deep sequencing when required or reducing sequencing costs per sample to an economically trivial amount when deep coverage is not needed. We demonstrate the utility and versatility of our approaches with results from six projects using different implementations of our protocols. Thus, we show that these methods facilitate amplicon library construction for Illumina instruments at reduced cost with increased flexibility. A simple web page to design fusion primers compatible with iTru primers is available at: http://baddna.uga.edu/tools-taggi.html. A fast and easy to use program to demultiplex amplicon pools with internal indexes is available at: https://github.com/lefeverde/Mr_Demuxy.

Direct and indirect effects of climate on richness drive the latitudinal diversity gradient in forest trees

Climate is widely recognised as an important determinant of the latitudinal diversity gradient. However, most existing studies make no distinction between direct and indirect effects of climate, which substantially hinders our understanding of how climate constrains biodiversity globally. Using data from 35 large forest plots, we test hypothesised relationships amongst climate, topography, forest structural attributes (stem abundance, tree size variation and stand basal area) and tree species richness to better understand drivers of latitudinal tree diversity patterns. Climate influences tree richness both directly, with more species in warm, moist, aseasonal climates and indirectly, with more species at higher stem abundance. These results imply direct limitation of species diversity by climatic stress and more rapid (co-)evolution and narrower niche partitioning in warm climates. They also support the idea that increased numbers of individuals associated with high primary productivity are partitioned to support a greater number of species.

Effects of biotic interactions on tropical tree performance depend on abiotic conditions

Predicting biotic responses to environmental change requires understanding the joint effects of abiotic conditions and biotic interactions on community dynamics. One major challenge is to separate the potentially confounding effects of abiotic environmental variation and local biotic interactions on individual performance. The stress gradient hypothesis (SGH) addresses this issue directly by predicting that the effects of biotic interactions on performance become more positive as the abiotic environment becomes more stressful. It is unclear, however, how the predictions of the SGH apply to plants of differing functional strategies in diverse communities. We asked (1) how the effect of crowding on performance (growth and survival) of trees varies across a precipitation gradient, and (2) how functional strategies (as measured by two key traits: wood density and leaf mass per area, LMA) mediate average demographic rates and responses to crowding across the gradient. We built trait-based neighborhood models of growth and survival across a regional precipitation gradient where increasing precipitation is associated with reduced abiotic stress. In total, our dataset comprised ~170,000 individual trees belonging to 252 species. The effect of crowding on tree performance varied across the gradient; crowding negatively affected growth across plots and positively affected survival in the wettest plot. Functional traits mediated average demographic rates across the gradient, but we did not find clear evidence that the strength of these responses depends on species' traits. Our study lends support to the SGH and demonstrates how a trait-based perspective can advance these concepts by linking the diversity of species interactions with functional variation across abiotic gradients.

Contrasting outcomes of species- and community-level analyses of the temporal consistency of functional composition

Multiple anthropogenic drivers affect every natural community, and there is broad interest in using functional traits to understand and predict the consequences for future biodiversity. There is, however, no consensus regarding the choice of analytical methods. We contrast species- and community-level analyses of change in the functional composition for four traits related to drought tolerance using three decades of repeat censuses of trees in the 50-ha Forest Dynamics Plot on Barro Colorado Island, Panama. Community trait distributions shifted significantly through time, which may indicate a shift toward more drought tolerant species. However, at the species level, changes in abundance were unrelated to trait values. To reconcile these seemingly contrasting results, we evaluated species-specific contributions to the directional shifts observed at the community level. Abundance changes of just one to six of 312 species were responsible for the community-level shifts observed for each trait. Our results demonstrate that directional changes in community-level functional composition can result from idiosyncratic change in a few species rather than widespread community-wide changes associated with functional traits. Future analyses of directional change in natural communities should combine community-, species-, and possibly individual-level analyses to uncover relationships with function that can improve understanding and enable prediction.

Community-level species' correlated distribution can be scale-independent and related to the evenness of abundance

The spatial distribution of species is not random; instead, individuals tend to gather, resulting in a non-random pattern. Previous studies used the independent negative binomial distribution (NBD) to model the distributional aggregation of a single species, in which the independence of the distribution of individuals of a species in different quadrats had been assumed. This way of analyzing aggregation will result in the scale-dependent estimation of the aggregation or shape parameter. However, because non-random (and therefore non-independent) distribution of individuals of a species in a finite area can be caused by either correlated or clumped distribution of individuals of a species between neighboring sites, an alternative model would assume that the distribution of individuals of a species over different sampling areas is multinomial. Here, we showed that, by assuming that regional species abundance followed a NBD while using a multinomial distribution to assign individuals of species in different non-overlapped sampling quadrats that are from a partition of the entire region (quantifying positive correlation or synchrony), the estimation of the shape parameter in this probabilistic model, which is the negative multinomial distribution (NMD), was scale-invariant (i.e., the estimated shape parameter is identical across different partitions of the study region). Accordingly, the estimation of the shape parameter was related to regional species distribution alone. This implied that, the shape parameter at the community level, using the NMD model, reflected the evenness of interspecific abundance. As a comparison, if the distribution of individuals of a single species followed independent NBDs as studied previously, the shape parameter would measure the evenness of intraspecific abundance (quantifying single-species' distributional aggregation). Moreover, our study highlighted the necessity for adjusting the model for the effects of unsampled species when studying community-level distributional patterns. Collectively, as long as a target area is partitioned into non-overlapping quadrats (no matter how their sizes vary), the proposed NMD model in this study, along with the independent NBDs model, can be jointly formulated as a framework to reconcile the scale-dependent debate on the shape parameter, unifying the relationship between inter- or intraspecific abundance and distributional patterns.

Resolving the paradox of clumped seed dispersal: positive density and distance dependence in a bat-dispersed species

One of the hypothesized benefits of seed dispersal is to escape density- and distance-responsive, host-specific, natural enemies near maternal plants where conspecific seed and seedling densities are high. Such high conspecific neighbor densities typically result in lower offspring growth and survival (i.e., negative density-dependent effects), yet many dispersal modes result in clumped seed distributions. New World leaf-nosed bats transport fruits to their feeding roosts and deposit seeds, thereby creating high-density seed/seedling patches beneath feeding roosts in heterospecific trees away from maternal trees, which seemingly nullifies a key benefit of seed dispersal. Such dispersal may still be adaptive if negative density-dependent effects are reduced under feeding roosts or if the benefit of being dispersed away from maternal trees outweighs negative effects of conspecific seed/seedling density below roosts. We mapped the entire post-germination population of a bat-dispersed tree species Calophyllum longifolium (Calophyllaceae) in a 50-ha plot on Barro Colorado Island, Panama in each of three successive years. We tested two hypotheses: (1) distance-dependent effects are stronger than density-dependent effects on seedling performance because seedlings far from conspecific adults are more likely to escape natural enemies even when at high densities and (2) negative density-dependent effects will be reduced far from vs. near conspecific adults. Density and distance were naturally decoupled, as expected. However, in contrast to our expectation, we found positive density effects on seedling survival and density-dependent effects did not differ with distance from conspecific adults. Both density and distance had positive effects on seedling survival when considered together, while only year had a significant effect on seedling growth. Thus, both being dispersed under bat feeding roosts and escaping the vicinity of conspecific adults were beneficial for C. longifolium seedling survival, supporting the directed dispersal and escape hypotheses, respectively. Despite resulting in high densities of conspecific seedlings, favorable habitat under bat feeding roosts and lack of negative density-dependent effects appear to provide evolutionary advantages in C. longifolium.

Role of tree size in moist tropical forest carbon cycling and water deficit responses

Drought disproportionately affects larger trees in tropical forests, but implications for forest composition and carbon (C) cycling in relation to dry season intensity remain poorly understood. In order to characterize how C cycling is shaped by tree size and drought adaptations and how these patterns relate to spatial and temporal variation in water deficit, we analyze data from three forest dynamics plots spanning a moisture gradient in Panama that have experienced El Niño droughts. At all sites, aboveground C cycle contributions peaked below 50-cm stem diameter, with stems ≥ 50 cm accounting for on average 59% of live aboveground biomass, 45% of woody productivity and 49% of woody mortality. The dominance of drought-avoidance strategies increased interactively with stem diameter and dry season intensity. Although size-related C cycle contributions did not vary systematically across the moisture gradient under nondrought conditions, woody mortality of larger trees was disproportionately elevated under El Niño drought stress. Thus, large (> 50 cm) stems, which strongly mediate but do not necessarily dominate C cycling, have drought adaptations that compensate for their more challenging hydraulic environment, particularly in drier climates. However, these adaptations do not fully buffer the effects of severe drought, and increased large tree mortality dominates ecosystem-level drought responses.

Variation in hydroclimate sustains tropical forest biomass and promotes functional diversity

The fate of tropical forests under climate change is unclear as a result, in part, of the uncertainty in projected changes in precipitation and in the ability of vegetation models to capture the effects of drought-induced mortality on aboveground biomass (AGB). We evaluated the ability of a terrestrial biosphere model with demography and hydrodynamics (Ecosystem Demography, ED2-hydro) to simulate AGB and mortality of four tropical tree plant functional types (PFTs) that operate along light- and water-use axes. Model predictions were compared with observations of canopy trees at Barro Colorado Island (BCI), Panama. We then assessed the implications of eight hypothetical precipitation scenarios, including increased annual precipitation, reduced inter-annual variation, El Niño-related droughts and drier wet or dry seasons, on AGB and functional diversity of the model forest. When forced with observed meteorology, ED2-hydro predictions capture multiple BCI benchmarks. ED2-hydro predicts that AGB will be sustained under lower rainfall via shifts in the functional composition of the forest, except under the drier dry-season scenario. These results support the hypothesis that inter-annual variation in mean and seasonal precipitation promotes the coexistence of functionally diverse PFTs because of the relative differences in mortality rates. If the hydroclimate becomes chronically drier or wetter, functional evenness related to drought tolerance may decline.

BioTIME: A database of biodiversity time series for the Anthropocene

MOTIVATION

The BioTIME database contains raw data on species identities and abundances in ecological assemblages through time. These data enable users to calculate temporal trends in biodiversity within and amongst assemblages using a broad range of metrics. BioTIME is being developed as a community-led open-source database of biodiversity time series. Our goal is to accelerate and facilitate quantitative analysis of temporal patterns of biodiversity in the Anthropocene.

MAIN TYPES OF VARIABLES INCLUDED

The database contains 8,777,413 species abundance records, from assemblages consistently sampled for a minimum of 2 years, which need not necessarily be consecutive. In addition, the database contains metadata relating to sampling methodology and contextual information about each record.

SPATIAL LOCATION AND GRAIN

BioTIME is a global database of 547,161 unique sampling locations spanning the marine, freshwater and terrestrial realms. Grain size varies across datasets from 0.0000000158 km2 (158 cm2) to 100 km2 (1,000,000,000,000 cm2).

TIME PERIOD AND GRAIN

BioTIME records span from 1874 to 2016. The minimal temporal grain across all datasets in BioTIME is a year.

MAJOR TAXA AND LEVEL OF MEASUREMENT

BioTIME includes data from 44,440 species across the plant and animal kingdoms, ranging from plants, plankton and terrestrial invertebrates to small and large vertebrates.

SOFTWARE FORMAT

.csv and .SQL.

Abiotic niche partitioning and negative density dependence drive tree seedling survival in a tropical forest

In tropical tree communities, processes occurring during early life stages play a critical role in shaping forest composition and diversity through differences in species' performance. Predicting the future of tropical forests depends on a solid understanding of the drivers of seedling survival. At the same time, factors determining spatial and temporal patterns of seedling survival can play a large role in permitting species coexistence in diverse communities. Using long-term data on the survival of more than 45 000 seedlings of 238 species in a Neotropical forest, we assessed the relative importance of key abiotic and biotic neighbourhood variables thought to influence individual seedling survival and tested whether species vary significantly in their responses to these variables, consistent with niche differences. At the community level, seedling survival was significantly correlated with plant size, topographic habitat, neighbourhood densities of conspecific seedlings, conspecific and heterospecific trees and annual variation in water availability, in descending order of effect size. Additionally, we found significant variation among species in their sensitivity to light and water availability, as well as in their survival within different topographic habitats, indicating the potential for niche differentiation among species that could allow for species coexistence.

Response to Comment on "Plant diversity increases with the strength of negative density dependence at the global scale"

Hülsmann and Hartig suggest that ecological mechanisms other than specialized natural enemies or intraspecific competition contribute to our estimates of conspecific negative density dependence (CNDD). To address their concern, we show that our results are not the result of a methodological artifact and present a null-model analysis that demonstrates that our original findings-(i) stronger CNDD at tropical relative to temperate latitudes and (ii) a latitudinal shift in the relationship between CNDD and species abundance-persist even after controlling for other processes that might influence spatial relationships between adults and recruits.

Plant diversity increases with the strength of negative density dependence at the global scale

Theory predicts that higher biodiversity in the tropics is maintained by specialized interactions among plants and their natural enemies that result in conspecific negative density dependence (CNDD). By using more than 3000 species and nearly 2.4 million trees across 24 forest plots worldwide, we show that global patterns in tree species diversity reflect not only stronger CNDD at tropical versus temperate latitudes but also a latitudinal shift in the relationship between CNDD and species abundance. CNDD was stronger for rare species at tropical versus temperate latitudes, potentially causing the persistence of greater numbers of rare species in the tropics. Our study reveals fundamental differences in the nature of local-scale biotic interactions that contribute to the maintenance of species diversity across temperate and tropical communities.

Adult mortality in a low-density tree population using high-resolution remote sensing

We developed a statistical framework to quantify mortality rates in canopy trees observed using time series from high-resolution remote sensing. By timing the acquisition of remote sensing data with synchronous annual flowering in the canopy tree species Handroanthus guayacan, we made 2,596 unique detections of 1,006 individual adult trees within 18,883 observation attempts on Barro Colorado Island, Panama (BCI) during an 11-yr period. There were 1,057 observation attempts that resulted in missing data due to cloud cover or incomplete spatial coverage. Using the fraction of 123 individuals from an independent field sample that were detected by satellite data (109 individuals, 88.6%), we estimate that the adult population for this species on BCI was 1,135 individuals. We used a Bayesian state-space model that explicitly accounted for the probability of tree detection and missing observations to compute an annual adult mortality rate of 0.2%·yr^-1 (SE = 0.1, 95% CI = 0.06-0.45). An independent estimate of the adult mortality rate from 260 field-checked trees closely matched the landscape-scale estimate (0.33%·yr^-1 , SE = 0.16, 95% CI = 0.12-0.74). Our proof-of-concept study shows that one can remotely estimate adult mortality rates for canopy tree species precisely in the presence of variable detection and missing observations.

The use of soil pollen to determine the sex of overhead individuals of a temperate dioecious shrub

PREMISE OF THE STUDY

In dioecious species, determining the sex of individual plants from one-time phenological observations is rarely feasible when some individuals capable of reproducing are not flowering or fruiting at the time of observation. Currently, sexing those individuals requires long-term phenological data on individuals and populations, but such data are rarely available or feasible to collect. We tested the hypothesis that differences in soil pollen concentrations beneath the crowns of female and male plants would exist and be sufficient to reliably determine the sex of the individual plant overhead in a dioecious species. We predicted that soil pollen concentrations beneath male plants would be significantly higher than beneath female plants because only males produce pollen and pollen should accumulate in the soil underneath the male plants over repeated flowering events.

METHODS

We collected samples from surface soil under both sexes of the insect-pollinated dioecious shrub, Aucuba japonica (Garryaceae).

KEY RESULTS

Pollen grains were present in surface soil in both Oe and A horizons, and mean pollen concentration under males was significantly higher than under females. Pollen concentrations beneath males were positively correlated with male plant height, potentially reflecting greater pollen production by larger individuals.

CONCLUSIONS

Considering the small plant size and relatively low pollen production of A. japonica, this method may hold promise for sexing other dioecious species in the absence of direct phenological data. Our phenology-free and relatively low-cost method for sexing dioecious plants may be especially useful in tropical forests where many species are dioecious.

Surviving in a Cosexual World: A Cost-Benefit Analysis of Dioecy in Tropical Trees

Dioecy has a demographic disadvantage compared with hermaphroditism: only about half of reproductive adults produce seeds. Dioecious species must therefore have fitness advantages to compensate for this cost through increased survival, growth, and/or reproduction. We used a full life cycle approach to quantify the demographic costs and benefits associated with dioecy while controlling for demographic differences between dioecious and hermaphroditic species related to other functional traits. The advantage of this novel approach is that we can focus on the effect of breeding system across a diverse tree community. We built a composite integral projection model for hermaphroditic and dioecious tree populations from Barro Colorado Island, Panama, using long-term demographic and newly collected reproductive data. Integration of all costs and benefits showed that compensation was realized through increased seed production, resulting in no net costs of dioecy. Compensation was also facilitated by the low contribution of reproduction to population growth. Estimated positive effects of dioecy on tree growth and survival were small and insignificant for population growth rates. Our model revealed that, for long-lived organisms, the cost of having males is smaller than generally expected. Hence, little compensation is required for dioecious species to maintain population growth rates similar to those of hermaphroditic species.

Use of sonic tomography to detect and quantify wood decay in living trees

PREMISE OF THE STUDY

Field methodology and image analysis protocols using acoustic tomography were developed and evaluated as a tool to estimate the amount of internal decay and damage of living trees, with special attention to tropical rainforest trees with irregular trunk shapes.

METHODS AND RESULTS

Living trunks of a diversity of tree species in tropical rainforests in the Republic of Panama were scanned using an Argus Electronic PiCUS 3 Sonic Tomograph and evaluated for the amount and patterns of internal decay. A protocol using ImageJ analysis software was used to quantify the proportions of intact and compromised wood. The protocols provide replicable estimates of internal decay and cavities for trees of varying shapes, wood density, and bark thickness.

CONCLUSIONS

Sonic tomography, coupled with image analysis, provides an efficient, noninvasive approach to evaluate decay patterns and structural integrity of even irregularly shaped living trees.

Positive effects of neighborhood complementarity on tree growth in a Neotropical forest

Numerous grassland experiments have found evidence for a complementarity effect, an increase in productivity with higher plant species richness due to niche partitioning. However, empirical tests of complementarity in natural forests are rare. We conducted a spatially explicit analysis of 518 433 growth records for 274 species from a 50-ha tropical forest plot to test neighborhood complementarity, the idea that a tree grows faster when it is surrounded by more dissimilar neighbors. We found evidence for complementarity: focal tree growth rates increased by 39.8% and 34.2% with a doubling of neighborhood multi-trait dissimilarity and phylogenetic dissimilarity, respectively. Dissimilarity from neighbors in maximum height had the most important effect on tree growth among the six traits examined, and indeed, its effect trended much larger than that of the multitrait dissimilarity index. Neighborhood complementarity effects were strongest for light-demanding species, and decreased in importance with increasing shade tolerance of the focal individuals. Simulations demonstrated that the observed neighborhood complementarities were sufficient to produce positive stand-level biodiversity-productivity relationships. We conclude that neighborhood complementarity is important for productivity in this tropical forest, and that scaling down to individual-level processes can advance our understanding of the mechanisms underlying stand-level biodiversity-productivity relationships.

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.

Relating belowground microbial composition to the taxonomic, phylogenetic, and functional trait distributions of trees in a tropical forest

The complexities of the relationships between plant and soil microbial communities remain unresolved. We determined the associations between plant aboveground and belowground (root) distributions and the communities of soil fungi and bacteria found across a diverse tropical forest plot. Soil microbial community composition was correlated with the taxonomic and phylogenetic structure of the aboveground plant assemblages even after controlling for differences in soil characteristics, but these relationships were stronger for fungi than for bacteria. In contrast to expectations, the species composition of roots in our soil core samples was a poor predictor of microbial community composition perhaps due to the patchy, ephemeral, and highly overlapping nature of fine root distributions. Our ability to predict soil microbial composition was not improved by incorporating information on plant functional traits suggesting that the most commonly measured plant traits are not particularly useful for predicting the plot-level variability in belowground microbial communities.

CTFS-ForestGEO: a worldwide network monitoring forests in an era of global change

Global change is impacting forests worldwide, threatening biodiversity and ecosystem services including climate regulation. Understanding how forests respond is critical to forest conservation and climate protection. This review describes an international network of 59 long-term forest dynamics research sites (CTFS-ForestGEO) useful for characterizing forest responses to global change. Within very large plots (median size 25 ha), all stems ≥ 1 cm diameter are identified to species, mapped, and regularly recensused according to standardized protocols. CTFS-ForestGEO spans 25 °S-61 °N latitude, is generally representative of the range of bioclimatic, edaphic, and topographic conditions experienced by forests worldwide, and is the only forest monitoring network that applies a standardized protocol to each of the world's major forest biomes. Supplementary standardized measurements at subsets of the sites provide additional information on plants, animals, and ecosystem and environmental variables. CTFS-ForestGEO sites are experiencing multifaceted anthropogenic global change pressures including warming (average 0.61 °C), changes in precipitation (up to ± 30% change), atmospheric deposition of nitrogen and sulfur compounds (up to 3.8 g N m(-2) yr(-1) and 3.1 g S m(-2) yr(-1)), and forest fragmentation in the surrounding landscape (up to 88% reduced tree cover within 5 km). The broad suite of measurements made at CTFS-ForestGEO sites makes it possible to investigate the complex ways in which global change is impacting forest dynamics. Ongoing research across the CTFS-ForestGEO network is yielding insights into how and why the forests are changing, and continued monitoring will provide vital contributions to understanding worldwide forest diversity and dynamics in an era of global change.

Trade-offs in seedling growth and survival within and across tropical forest microhabitats

For niche differences to maintain coexistence of sympatric species, each species must grow and/or survive better than each of the others in at least one set of conditions (i.e., performance trade-offs). However, the extent of niche differentiation in tropical forests remains highly debated. We present the first test of performance trade-offs for wild seedlings in a tropical forest. We measured seedling relative growth rate (RGR) and survival of four common native woody species across 18 light, substrate, and topography microhabitats over 2.5 years within Hawaiian montane wet forest, an ideal location due to its low species diversity and strong species habitat associations. All six species pairs exhibited significant performance trade-offs across microhabitats and for RGR versus survival within microhabitats. We also found some evidence of performance equivalence, with species pairs having similar performance in 26% of comparisons across microhabitats. Across species, survival under low light was generally positively associated with RGR under high light. When averaged over all species, topography (slope, aspect, and elevation) explained most of the variation in RGR attributable to microhabitat variables (51-53%) followed by substrate type (35-37%) and light (11-12%). However, the relative effects of microhabitat differed among species and RGR metric (i.e., RGR for height, biomass, or leaf area). These findings indicate that performance trade-offs among species during regeneration are common in low-diversity tropical forest, although other mechanisms may better explain the coexistence of species with small performance differences.

Habitat filtering across tree life stages in tropical forest communities

Tropical tree communities are shaped by local-scale habitat heterogeneity in the form of topographic and edaphic variation, but the life-history stage at which habitat associations develop remains poorly understood. This is due, in part, to the fact that previous studies have not accounted for the widely disparate sample sizes (number of stems) that result when trees are divided into size classes. We demonstrate that the observed habitat structuring of a community is directly related to the number of individuals in the community. We then compare the relative importance of habitat heterogeneity to tree community structure for saplings, juveniles and adult trees within seven large (24-50 ha) tropical forest dynamics plots while controlling for sample size. Changes in habitat structuring through tree life stages were small and inconsistent among life stages and study sites. Where found, these differences were an order of magnitude smaller than the findings of previous studies that did not control for sample size. Moreover, community structure and composition were very similar among tree sub-communities of different life stages. We conclude that the structure of these tropical tree communities is established by the time trees are large enough to be included in the census (1 cm diameter at breast height), which indicates that habitat filtering occurs during earlier life stages.

Tropical rain forest conservation and the twin challenges of diversity and rarity

Data from a global network of large, permanent plots in lowland tropical forests demonstrate (1) that the phenomenon of tropical tree rarity is real and (2) that almost all the species diversity in such forests is due to rare species. Theoretical and empirically based reasoning suggests that many of these rare species are not as geographically widespread as previously thought. These findings suggest that successful strategies for conserving global tree diversity in lowland tropical forests must pay much more attention to the biogeography of rarity, as well as to the impact of climate change on the distribution and abundance of rare species. Because the biogeography of many tropical tree species is poorly known, a high priority should be given to documenting the distribution and abundance of rare tropical tree species, particularly in Amazonia, the largest remaining tropical forested region in the world.

A taxonomic comparison of local habitat niches of tropical trees

The integration of ecology and evolutionary biology requires an understanding of the evolutionary lability in species' ecological niches. For tropical trees, specialization for particular soil resource and topographic conditions is an important part of the habitat niche, influencing the distributions of individual species and overall tree community structure at the local scale. However, little is known about how these habitat niches are related to the evolutionary history of species. We assessed the relationship between taxonomic rank and tree species' soil resource and topographic niches in eight large (24-50 ha) tropical forest dynamics plots. Niche overlap values, indicating the similarity of two species' distributions along soil or topographic axes, were calculated for all pairwise combinations of co-occurring tree species at each study site. Congeneric species pairs often showed greater niche overlap (i.e., more similar niches) than non-congeneric pairs along both soil and topographic axes, though significant effects were found for only five sites based on Mantel tests. No evidence for taxonomic effects was found at the family level. Our results indicate that local habitat niches of trees exhibit varying degrees of phylogenetic signal at different sites, which may have important ramifications for the phylogenetic structure of these communities.

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.