Soils and topography control natural disturbance rates and thereby forest structure in a lowland tropical landscape.

Tree mortality is a major control over tropical forest carbon stocks globally but the strength of associations between abiotic drivers and tree mortality within forested landscapes is poorly understood. Here, we used repeat drone photogrammetry across 1500 ha of forest in Central Panama over 5 years to quantify spatial variation in canopy disturbance rates and its predictors. We identified 11,153 canopy disturbances greater than 25 m2 in area, including treefalls, large branchfalls and standing dead trees, affecting 1.9% of area per year. Soil type, forest age and topography explained up to 46%-67% of disturbance rate variation at spatial grains of 58-64 ha, with higher rates in older forests, steeper slopes and local depressions. Furthermore, disturbance rates predicted the proportion of low canopy area across the landscape, and mean canopy height in old growth forests. Thus abiotic factors drive variation in disturbance rates and thereby forest structure at landscape scales.

Simulating environmentally-sensitive tree recruitment in vegetation demographic models.

Vegetation demographic models (VDMs) endeavor to predict how global forests will respond to climate change. This requires simulating which trees, if any, are able to recruit under changing environmental conditions. We present a new recruitment scheme for VDMs in which functional-type-specific recruitment rates are sensitive to light, soil moisture and the productivity of reproductive trees. We evaluate the scheme by predicting tree recruitment for four tropical tree functional types under varying meteorology and canopy structure at Barro Colorado Island, Panama. We compare predictions to those of a current VDM, quantitative observations and ecological expectations. We find that the scheme improves the magnitude and rank order of recruitment rates among functional types and captures recruitment limitations in response to variable understory light, soil moisture and precipitation regimes. Our results indicate that adopting this framework will improve VDM capacity to predict functional-type-specific tree recruitment in response to climate change, thereby improving predictions of future forest distribution, composition and function.

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.

Tree crown damage and its effects on forest carbon cycling in a tropical forest.

Crown damage can account for over 23% of canopy biomass turnover in tropical forests and is a strong predictor of tree mortality; yet, it is not typically represented in vegetation models. We incorporate crown damage into the Functionally Assembled Terrestrial Ecosystem Simulator (FATES), to evaluate how lags between damage and tree recovery or death alter demographic rates and patterns of carbon turnover. We represent crown damage as a reduction in a tree's crown area and leaf and branch biomass, and allow associated variation in the ratio of aboveground to belowground plant tissue. We compare simulations with crown damage to simulations with equivalent instant increases in mortality and benchmark results against data from Barro Colorado Island (BCI), Panama. In FATES, crown damage causes decreases in growth rates that match observations from BCI. Crown damage leads to increases in carbon starvation mortality in FATES, but only in configurations with high root respiration and decreases in carbon storage following damage. Crown damage also alters competitive dynamics, as plant functional types that can recover from crown damage outcompete those that cannot. This is a first exploration of the trade-off between the additional complexity of the novel crown damage module and improved predictive capabilities. At BCI, a tropical forest that does not experience high levels of disturbance, both the crown damage simulations and simulations with equivalent increases in mortality does a reasonable job of capturing observations. The crown damage module provides functionality for exploring dynamics in forests with more extreme disturbances such as cyclones and for capturing the synergistic effects of disturbances that overlap in space and time.

Liana optical traits increase tropical forest albedo and reduce ecosystem productivity.

Lianas are a key growth form in tropical forests. Their lack of self-supporting tissues and their vertical position on top of the canopy make them strong competitors of resources. A few pioneer studies have shown that liana optical traits differ on average from those of colocated trees. Those trait discrepancies were hypothesized to be responsible for the competitive advantage of lianas over trees. Yet, in the absence of reliable modelling tools, it is impossible to unravel their impact on the forest energy balance, light competition, and on the liana success in Neotropical forests. To bridge this gap, we performed a meta-analysis of the literature to gather all published liana leaf optical spectra, as well as all canopy spectra measured over different levels of liana infestation. We then used a Bayesian data assimilation framework applied to two radiative transfer models (RTMs) covering the leaf and canopy scales to derive tropical tree and liana trait distributions, which finally informed a full dynamic vegetation model. According to the RTMs inversion, lianas grew thinner, more horizontal leaves with lower pigment concentrations. Those traits made the lianas very efficient at light interception and significantly modified the forest energy balance and its carbon cycle. While forest albedo increased by 14% in the shortwave, light availability was reduced in the understorey (-30% of the PAR radiation) and soil temperature decreased by 0.5°C. Those liana-specific traits were also responsible for a significant reduction of tree (-19%) and ecosystem (-7%) gross primary productivity (GPP) while lianas benefited from them (their GPP increased by +27%). This study provides a novel mechanistic explanation to the increase in liana abundance, new evidence of the impact of lianas on forest functioning, and paves the way for the evaluation of the large-scale impacts of lianas on forest biogeochemical cycles.

Patterns and mechanisms of spatial variation in tropical forest productivity, woody residence time, and biomass.

Tropical forests vary widely in biomass carbon (C) stocks and fluxes even after controlling for forest age. A mechanistic understanding of this variation is critical to accurately predicting responses to global change. We review empirical studies of spatial variation in tropical forest biomass, productivity and woody residence time, focusing on mature forests. Woody productivity and biomass decrease from wet to dry forests and with elevation. Within lowland forests, productivity and biomass increase with temperature in wet forests, but decrease with temperature where water becomes limiting. Woody productivity increases with soil fertility, whereas residence time decreases, and biomass responses are variable, consistent with an overall unimodal relationship. Areas with higher disturbance rates and intensities have lower woody residence time and biomass. These environmental gradients all involve both direct effects of changing environments on forest C fluxes and shifts in functional composition - including changing abundances of lianas - that substantially mitigate or exacerbate direct effects. Biogeographic realms differ significantly and importantly in productivity and biomass, even after controlling for climate and biogeochemistry, further demonstrating the importance of plant species composition. Capturing these patterns in global vegetation models requires better mechanistic representation of water and nutrient limitation, plant compositional shifts and tree mortality.

Leaf turgor loss point shapes local and regional distributions of evergreen but not deciduous tropical trees.

The effects of climate change on tropical forests will depend on how diverse tropical tree species respond to drought. Current distributions of evergreen and deciduous tree species across local and regional moisture gradients reflect their ability to tolerate drought stress, and might be explained by functional traits. We measured leaf water potential at turgor loss (i.e. 'wilting point'; πtlp ), wood density (WD) and leaf mass per area (LMA) on 50 of the most abundant tree species in central Panama. We then tested their ability to explain distributions of evergreen and deciduous species within a 50 ha plot on Barro Colorado Island and across a 70 km rainfall gradient spanning the Isthmus of Panama. Among evergreen trees, species with lower πtlp were associated with drier habitats, with πtlp explaining 28% and 32% of habitat association on local and regional scales, respectively, greatly exceeding the predictive power of WD and LMA. In contrast, πtlp did not predict habitat associations among deciduous species. Across spatial scales, πtlp is a useful indicator of habitat preference for tropical tree species that retain their leaves during periods of water stress, and holds the potential to predict vegetation responses to climate change.

Global patterns of forest autotrophic carbon fluxes.

Carbon (C) fixation, allocation, and metabolism by trees set the basis for energy and material flows in forest ecosystems and define their interactions with Earth's changing climate. However, while many studies have considered variation in productivity with latitude and climate, we lack a cohesive synthesis on how forest carbon fluxes vary globally with respect to climate and one another. Here, we draw upon 1,319 records from the Global Forest Carbon Database, representing all major forest types and the nine most significant autotrophic carbon fluxes, to comprehensively review how annual C cycling in mature, undisturbed forests varies with latitude and climate on a global scale. Across all flux variables analyzed, rates of C cycling decreased continuously with absolute latitude-a finding that confirms multiple previous studies and contradicts the idea that net primary productivity of temperate forests rivals that of tropical forests. C flux variables generally displayed similar trends across latitude and multiple climate variables, with no differences in allocation detected at this global scale. Temperature variables in general, and mean annual temperature or temperature seasonality in particular, were the best single predictors of C flux, explaining 19%-71% of variation in the C fluxes analyzed. The effects of temperature were modified by moisture availability, with C flux reduced under hot and dry conditions and sometimes under very high precipitation. Annual C fluxes increased with growing season length and were also influenced by growing season climate. These findings clarify how forest C flux varies with latitude and climate on a global scale. In an era when forests will play a critical yet uncertain role in shaping Earth's rapidly changing climate, our synthesis provides a foundation for understanding global patterns in forest C cycling.

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.

Pantropical geography of lightning-caused disturbance and its implications for tropical forests.

Lightning is a major agent of disturbance, but its ecological effects in the tropics are unquantified. Here we used ground and satellite sensors to quantify the geography of lightning strikes in terrestrial tropical ecosystems, and to evaluate whether spatial variation in lightning frequency is associated with variation in tropical forest structure and dynamics. Between 2013 and 2018, tropical terrestrial ecosystems received an average of 100.4 million lightning strikes per year, and the frequency of strikes was spatially autocorrelated at local-to-continental scales. Lightning strikes were more frequent in forests, savannas, and urban areas than in grasslands, shrublands, and croplands. Higher lightning frequency was positively associated with woody biomass turnover and negatively associated with aboveground biomass and the density of large trees (trees/ha) in forests across Africa, Asia, and the Americas. Extrapolating from the only tropical forest study that comprehensively assessed tree damage and mortality from lightning strikes, we estimate that lightning directly damages c. 832 million trees in tropical forests annually, of which c. 194 million die. The similarly high lightning frequency in tropical savannas suggests that lightning also influences savanna tree mortality rates and ecosystem processes. These patterns indicate that lightning-caused disturbance plays a major and largely unappreciated role in pantropical ecosystem dynamics and global carbon cycling.

Allometric constraints and competition enable the simulation of size structure and carbon fluxes in a dynamic vegetation model of tropical forests (LM3PPA-TV).

Tropical forests are a key determinant of the functioning of the Earth system, but remain a major source of uncertainty in carbon cycle models and climate change projections. In this study, we present an updated land model (LM3PPA-TV) to improve the representation of tropical forest structure and dynamics in Earth system models (ESMs). The development and parameterization of LM3PPA-TV drew on extensive datasets on tropical tree traits and long-term field censuses from Barro Colorado Island (BCI), Panama. The model defines a new plant functional type (PFT) based on the characteristics of shade-tolerant, tropical tree species, implements a new growth allocation scheme based on realistic tree allometries, incorporates hydraulic constraints on biomass accumulation, and features a new compartment for tree branches and branch fall dynamics. Simulation experiments reproduced observed diurnal and seasonal patterns in stand-level carbon and water fluxes, as well as mean canopy and understory tree growth rates, tree size distributions, and stand-level biomass on BCI. Simulations at multiple sites captured considerable variation in biomass and size structure across the tropical forest biome, including observed responses to precipitation and temperature. Model experiments suggested a major role of water limitation in controlling geographic variation forest biomass and structure. However, the failure to simulate tropical forests under extreme conditions and the systematic underestimation of forest biomass in Paleotropical locations highlighted the need to incorporate variation in hydraulic traits and multiple PFTs that capture the distinct floristic composition across tropical domains. The continued pressure on tropical forests from global change demands models which are able to simulate alternative successional pathways and their pace to recovery. LM3PPA-TV provides a tool to investigate geographic variation in tropical forests and a benchmark to continue improving the representation of tropical forests dynamics and their carbon storage potential in ESMs.

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 (EAGBi ) in 10 × 10 m quadrats. Productivity and mortality fluxes both increased strongly with initial quadrat EAGB. The distribution of EAGB (and thus EAGBi ) across quadrats hardly varied over 30 years (and seven censuses). EAGB fluxes as a function of EAGBi 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 EAGBi , with essentially no contribution from changes in EAGBi 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.

Signs of stabilisation and stable coexistence.

Many empirical studies motivated by an interest in stable coexistence have quantified negative density dependence, negative frequency dependence, or negative plant-soil feedback, but the links between these empirical results and ecological theory are not straightforward. Here, we relate these analyses to theoretical conditions for stabilisation and stable coexistence in classical competition models. By stabilisation, we mean an excess of intraspecific competition relative to interspecific competition that inherently slows or even prevents competitive exclusion. We show that most, though not all, tests demonstrating negative density dependence, negative frequency dependence, and negative plant-soil feedback constitute sufficient conditions for stabilisation of two-species interactions if applied to data for per capita population growth rates of pairs of species, but none are necessary or sufficient conditions for stable coexistence of two species. Potential inferences are even more limited when communities involve more than two species, and when performance is measured at a single life stage or vital rate. We then discuss two approaches that enable stronger tests for stable coexistence-invasibility experiments and model parameterisation. The model parameterisation approach can be applied to typical density-dependence, frequency-dependence, and plant-soil feedback data sets, and generally enables better links with mechanisms and greater insights, as demonstrated by recent studies.

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.

Drivers and mechanisms of tree mortality in moist tropical forests.

Tree mortality rates appear to be increasing in moist tropical forests (MTFs) with significant carbon cycle consequences. Here, we review the state of knowledge regarding MTF tree mortality, create a conceptual framework with testable hypotheses regarding the drivers, mechanisms and interactions that may underlie increasing MTF mortality rates, and identify the next steps for improved understanding and reduced prediction. Increasing mortality rates are associated with rising temperature and vapor pressure deficit, liana abundance, drought, wind events, fire and, possibly, CO2 fertilization-induced increases in stand thinning or acceleration of trees reaching larger, more vulnerable heights. The majority of these mortality drivers may kill trees in part through carbon starvation and hydraulic failure. The relative importance of each driver is unknown. High species diversity may buffer MTFs against large-scale mortality events, but recent and expected trends in mortality drivers give reason for concern regarding increasing mortality within MTFs. Models of tropical tree mortality are advancing the representation of hydraulics, carbon and demography, but require more empirical knowledge regarding the most common drivers and their subsequent mechanisms. We outline critical datasets and model developments required to test hypotheses regarding the underlying causes of increasing MTF mortality rates, and improve prediction of future mortality under climate change.

Vegetation demographics in Earth System Models: A review of progress and priorities.

Numerous current efforts seek to improve the representation of ecosystem ecology and vegetation demographic processes within Earth System Models (ESMs). These developments are widely viewed as an important step in developing greater realism in predictions of future ecosystem states and fluxes. Increased realism, however, leads to increased model complexity, with new features raising a suite of ecological questions that require empirical constraints. Here, we review the developments that permit the representation of plant demographics in ESMs, and identify issues raised by these developments that highlight important gaps in ecological understanding. These issues inevitably translate into uncertainty in model projections but also allow models to be applied to new processes and questions concerning the dynamics of real-world ecosystems. We argue that stronger and more innovative connections to data, across the range of scales considered, are required to address these gaps in understanding. The development of first-generation land surface models as a unifying framework for ecophysiological understanding stimulated much research into plant physiological traits and gas exchange. Constraining predictions at ecologically relevant spatial and temporal scales will require a similar investment of effort and intensified inter-disciplinary communication.

Tree diversity in relation to maximum tree height: evidence for the harshness hypothesis of species diversity gradients.

Marks et al. (Ecology Letters, 19, 2016, 743) showed tree species richness correlates with maximum tree height, and interpret this as evidence that the environmental stressors that limit tree height also act as ecological filters on species richness. Here, we strengthen these arguments by further addressing the roles of environmental covariates and beta diversity.

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.

Drought-induced mortality patterns and rapid biomass recovery in a terra firme forest in the Colombian Amazon.

Extreme climatic events affecting the Amazon region are expected to become more frequent under ongoing climate change. In this study, we assessed the responses to the 2010 drought of over 14,000 trees ≥10 cm dbh in a 25 ha lowland forest plot in the Colombian Amazon and how these responses varied among topographically defined habitats, with tree size, and with species wood density. Tree mortality was significantly higher during the 2010-2013 period immediately after the drought than in 2007-2010. The post-drought increase in mortality was stronger for trees located in valleys (+243%) than for those located on slopes (+67%) and ridges (+57%). Tree-based generalized linear mixed models showed a significant negative effect of species wood density on mortality and no effect of tree size. Despite the elevated post-drought mortality, aboveground biomass increased from 2007 to 2013 by 1.62 Mg ha-1  yr-1 (95% CI 0.80-2.43 Mg ha-1  yr-1 ). Biomass change varied among habitats, with no significant increase on the slopes (1.05, 95% CI -0.76 to 2.85 Mg ha-1  yr-1 ), a significant increase in the valleys (1.33, 95% CI 0.37-2.34 Mg ha-1  yr-1 ), and a strong increase on the ridges (2.79, 95% CI 1.20-4.21 Mg ha-1  yr-1 ). These results indicate a high carbon resilience of this forest to the 2010 drought due to habitat-associated and interspecific heterogeneity in responses including directional changes in functional composition driven by enhanced performance of drought-tolerant species that inhabit the drier ridges.

Cascading effects of defaunation on the coexistence of two specialized insect seed predators.

Identification of the mechanisms enabling stable coexistence of species with similar resource requirements is a central challenge in ecology. Such coexistence can be facilitated by species at higher trophic levels through complex multi-trophic interactions, a mechanism that could be compromised by ongoing defaunation. We investigated cascading effects of defaunation on Pachymerus cardo and Speciomerus giganteus, the specialized insect seed predators of the Neotropical palm Attalea butyracea, testing the hypothesis that vertebrate frugivores and granivores facilitate their coexistence. Laboratory experiments showed that the two seed parasitoid species differed strongly in their reproductive ecology. Pachymerus produced many small eggs that it deposited exclusively on the fruit exocarp (exterior). Speciomerus produced few large eggs that it deposited exclusively on the endocarp, which is normally exposed only after a vertebrate handles the fruit. When eggs of the two species were deposited on the same fruit, Pachymerus triumphed only when it had a long head start, and the loser always succumbed to intraguild predation. We collected field data on the fates of 6569 Attalea seeds across sites in central Panama with contrasting degrees of defaunation and wide variation in the abundance of vertebrate frugivores and granivores. Speciomerus dominated where vertebrate communities were intact, whereas Pachymerus dominated in defaunated sites. Variation in the relative abundance of Speciomerus across all 84 sampling sites was strongly positively related to the proportion of seeds attacked by rodents, an indicator of local vertebrate abundance. SYNTHESIS: We show that two species of insect seed predators relying on the same host plant species are niche differentiated in their reproductive strategies such that one species has the advantage when fruits are handled promptly by vertebrates and the other when they are not. Defaunation disrupts this mediating influence of vertebrates and strongly favours one species at the expense of the other, providing a case study of the cascading effects of defaunation and its potential to disrupt coexistence of non-target species, including the hyperdiverse phytophagous insects of tropical forests.