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.

Damage to living trees contributes to almost half of the biomass losses in tropical forests.

Accurate estimates of forest biomass stocks and fluxes are needed to quantify global carbon budgets and assess the response of forests to climate change. However, most forest inventories consider tree mortality as the only aboveground biomass (AGB) loss without accounting for losses via damage to living trees: branchfall, trunk breakage, and wood decay. Here, we use ~151,000 annual records of tree survival and structural completeness to compare AGB loss via damage to living trees to total AGB loss (mortality + damage) in seven tropical forests widely distributed across environmental conditions. We find that 42% (3.62 Mg ha-1 year-1 ; 95% confidence interval [CI] 2.36-5.25) of total AGB loss (8.72 Mg ha-1 year-1 ; CI 5.57-12.86) is due to damage to living trees. Total AGB loss was highly variable among forests, but these differences were mainly caused by site variability in damage-related AGB losses rather than by mortality-related AGB losses. We show that conventional forest inventories overestimate stand-level AGB stocks by 4% (1%-17% range across forests) because assume structurally complete trees, underestimate total AGB loss by 29% (6%-57% range across forests) due to overlooked damage-related AGB losses, and overestimate AGB loss via mortality by 22% (7%-80% range across forests) because of the assumption that trees are undamaged before dying. Our results indicate that forest carbon fluxes are higher than previously thought. Damage on living trees is an underappreciated component of the forest carbon cycle that is likely to become even more important as the frequency and severity of forest disturbances increase.

Individual tree damage dominates mortality risk factors across six tropical forests.

The relative importance of tree mortality risk factors remains unknown, especially in diverse tropical forests where species may vary widely in their responses to particular conditions. We present a new framework for quantifying the importance of mortality risk factors and apply it to compare 19 risks on 31 203 trees (1977 species) in 14 one-year periods in six tropical forests. We defined a condition as a risk factor for a species if it was associated with at least a doubling of mortality rate in univariate analyses. For each risk, we estimated prevalence (frequency), lethality (difference in mortality between trees with and without the risk) and impact ('excess mortality' associated with the risk, relative to stand-level mortality). The most impactful risk factors were light limitation and crown/trunk loss; the most prevalent were light limitation and small size; the most lethal were leaf damage and wounds. Modes of death (standing, broken and uprooted) had limited links with previous conditions and mortality risk factors. We provide the first ranking of importance of tree-level mortality risk factors in tropical forests. Future research should focus on the links between these risks, their climatic drivers and the physiological processes to enable mechanistic predictions of future tree mortality.

Hydraulically-vulnerable trees survive on deep-water access during droughts in a tropical forest.

Deep-water access is arguably the most effective, but under-studied, mechanism that plants employ to survive during drought. Vulnerability to embolism and hydraulic safety margins can predict mortality risk at given levels of dehydration, but deep-water access may delay plant dehydration. Here, we tested the role of deep-water access in enabling survival within a diverse tropical forest community in Panama using a novel data-model approach. We inversely estimated the effective rooting depth (ERD, as the average depth of water extraction), for 29 canopy species by linking diameter growth dynamics (1990-2015) to vapor pressure deficit, water potentials in the whole-soil column, and leaf hydraulic vulnerability curves. We validated ERD estimates against existing isotopic data of potential water-access depths. Across species, deeper ERD was associated with higher maximum stem hydraulic conductivity, greater vulnerability to xylem embolism, narrower safety margins, and lower mortality rates during extreme droughts over 35 years (1981-2015) among evergreen species. Species exposure to water stress declined with deeper ERD indicating that trees compensate for water stress-related mortality risk through deep-water access. The role of deep-water access in mitigating mortality of hydraulically-vulnerable trees has important implications for our predictive understanding of forest dynamics under current and future climates.

Phylogenetic classification of the world's tropical forests.

Knowledge about the biogeographic affinities of the world's tropical forests helps to better understand regional differences in forest structure, diversity, composition, and dynamics. Such understanding will enable anticipation of region-specific responses to global environmental change. Modern phylogenies, in combination with broad coverage of species inventory data, now allow for global biogeographic analyses that take species evolutionary distance into account. Here we present a classification of the world's tropical forests based on their phylogenetic similarity. We identify five principal floristic regions and their floristic relationships: (i) Indo-Pacific, (ii) Subtropical, (iii) African, (iv) American, and (v) Dry forests. Our results do not support the traditional neo- versus paleotropical forest division but instead separate the combined American and African forests from their Indo-Pacific counterparts. We also find indications for the existence of a global dry forest region, with representatives in America, Africa, Madagascar, and India. Additionally, a northern-hemisphere Subtropical forest region was identified with representatives in Asia and America, providing support for a link between Asian and American northern-hemisphere 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.

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.

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.

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.

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.

Asymmetric density dependence shapes species abundances in a tropical tree community.

The factors determining species commonness and rarity are poorly understood, particularly in highly diverse communities. Theory predicts that interactions with neighbors of the same (conspecific) and other (heterospecific) species can influence a species' relative abundance, but empirical tests are lacking. By using a hierarchical model of survival for more than 30,000 seedlings of 180 tropical tree species on Barro Colorado Island, Panama, we tested whether species' sensitivity to neighboring individuals relates to their relative abundance in the community. We found wide variation among species in the effect of conspecific, but not heterospecific, neighbors on survival, and we found a significant relationship between the strength of conspecific neighbor effects and species abundance. Specifically, rare species suffered more from the presence of conspecific neighbors than common species did, suggesting that conspecific density dependence shapes species abundances in diverse communities.

Error propagation and scaling for tropical forest biomass estimates.

The above-ground biomass (AGB) of tropical forests is a crucial variable for ecologists, biogeochemists, foresters and policymakers. Tree inventories are an efficient way of assessing forest carbon stocks and emissions to the atmosphere during deforestation. To make correct inferences about long-term changes in biomass stocks, it is essential to know the uncertainty associated with AGB estimates, yet this uncertainty is rarely evaluated carefully. Here, we quantify four types of uncertainty that could lead to statistical error in AGB estimates: (i) error due to tree measurement; (ii) error due to the choice of an allometric model relating AGB to other tree dimensions; (iii) sampling uncertainty, related to the size of the study plot; (iv) representativeness of a network of small plots across a vast forest landscape. In previous studies, these sources of error were reported but rarely integrated into a consistent framework. We estimate all four terms in a 50 hectare (ha, where 1 ha = 10(4) m2) plot on Barro Colorado Island, Panama, and in a network of 1 ha plots scattered across central Panama. We find that the most important source of error is currently related to the choice of the allometric model. More work should be devoted to improving the predictive power of allometric models for biomass.

An ecosystem report on the Panama Canal: monitoring the status of the forest communities and the watershed.

In 1996, the Smithsonian Tropical Research Institute and the Republic of Panama's Environmental Authority, with support from the United States Agency for International Development, undertook a comprehensive program to monitor the ecosystem of the Panama Canal watershed. The goals were to establish baseline indicators for the integrity of forest communities and rivers. Based on satellite image classification and ground surveys, the 2790 km2 watershed had 1570 km2 of forest in 1997, 1080 km2 of which was in national parks and nature monuments. Most of the 490 km2 of forest not currently in protected areas lies along the west bank of the Canal, and its management status after the year 2000 turnover of the Canal from the U.S. to Panama remains uncertain. In forest plots designed to monitor forest diversity and change, a total of 963 woody plant species were identified and mapped. We estimate there are a total of 850-1000 woody species in forests of the Canal corridor. Forests of the wetter upper reaches of the watershed are distinct in species composition from the Canal corridor, and have considerably higher diversity and many unknown species. These remote areas are extensively forested, poorly explored, and harbor an estimated 1400-2200 woody species. Vertebrate monitoring programs were also initiated, focusing on species threatened by hunting and forest fragmentation. Large mammals are heavily hunted in most forests of Canal corridor, and there was clear evidence that mammal density is greatly reduced in hunted areas and that this affects seed predation and dispersal. The human population of the watershed was 113 000 in 1990, and grew by nearly 4% per year from 1980 to 1990. Much of this growth was in a small region of the watershed on the outskirts of Panama City, but even rural areas, including villages near and within national parks, grew by 2% per year. There is no sewage treatment in the watershed, and many towns have no trash collection, thus streams near large towns are heavily polluted. Analyses of sediment loads in rivers throughout the watershed did not indicate that erosion has been increasing as a result of deforestation, rather, erosion seems to be driven largely by total rainfall and heavy rainfall events that cause landslides. Still, models suggest that large-scale deforestation would increase landslide frequency, and failure to detect increases in erosion could be due to the gradual deforestation rate and the short time period over which data are available. A study of runoff showed deforestation increased the amount of water from rainfall that passed directly into streams. As a result, dry season flow was reduced in a deforested catchment relative to a forested one. Currently, the Panama Canal watershed has extensive forest areas and streams relatively unaffected by humans. But impacts of hunting and pollution near towns are clear, and the burgeoning population will exacerbate these impacts in the next few decades. Changes in policies regarding forest protection and pollution control are necessary.

Beta-diversity in tropical forest trees.

The high alpha-diversity of tropical forests has been amply documented, but beta-diversity-how species composition changes with distance-has seldom been studied. We present quantitative estimates of beta-diversity for tropical trees by comparing species composition of plots in lowland terra firme forest in Panama, Ecuador, and Peru. We compare observations with predictions derived from a neutral model in which habitat is uniform and only dispersal and speciation influence species turnover. We find that beta-diversity is higher in Panama than in western Amazonia and that patterns in both areas are inconsistent with the neutral model. In Panama, habitat variation appears to increase species turnover relative to Amazonia, where unexpectedly low turnover over great distances suggests that population densities of some species are bounded by as yet unidentified processes. At intermediate scales in both regions, observations can be matched by theory, suggesting that dispersal limitation, with speciation, influences species turnover.