Does deterministic coexistence theory matter in a finite world?

Contemporary studies of species coexistence are underpinned by deterministic models that assume that competing species have continuous (i.e., noninteger) densities, live in infinitely large landscapes, and coexist over infinite time horizons. By contrast, in nature, species are composed of discrete individuals subject to demographic stochasticity and occur in habitats of finite size where extinctions occur in finite time. One consequence of these discrepancies is that metrics of species' coexistence derived from deterministic theory may be unreliable predictors of the duration of species coexistence in nature. These coexistence metrics include invasion growth rates and niche and fitness differences, which are now commonly applied in theoretical and empirical studies of species coexistence. In this study, we tested the efficacy of deterministic coexistence metrics on the duration of species coexistence in a finite world. We introduce new theoretical and computational methods to estimate coexistence times in stochastic counterparts of classic deterministic models of competition. Importantly, we parameterized this model using experimental field data for 90 pairwise combinations of 18 species of annual plants, allowing us to derive biologically informed estimates of coexistence times for a natural system. Strikingly, we found that for species expected to deterministically coexist, community sizes containing only 10 individuals had predicted coexistence times of more than 1000 years. We also found that invasion growth rates explained 60% of the variation in intrinsic coexistence times, reinforcing their general usefulness in studies of coexistence. However, only by integrating information on both invasion growth rates and species' equilibrium population sizes could most (>99%) of the variation in species coexistence times be explained. This integration was achieved with demographically uncoupled single-species models solely determined by the invasion growth rates and equilibrium population sizes. Moreover, because of a complex relationship between niche overlap/fitness differences and equilibrium population sizes, increasing niche overlap and increasing fitness differences did not always result in decreasing coexistence times, as deterministic theory would predict. Nevertheless, our results tend to support the informed use of deterministic theory for understanding the duration of species' coexistence while highlighting the need to incorporate information on species' equilibrium population sizes in addition to invasion growth rates

Trait-based tests of coexistence mechanisms

Abstract Recent functional trait studies have shown that trait differences may favour certain species (environmental filtering) while simultaneously preventing competitive exclusion (niche partitioning). However, phenomenological trait-dispersion analyses do not identify the mechanisms that generate niche partitioning, preventing trait-based prediction of future changes in biodiversity. We argue that such predictions require linking functional traits with recognised coexistence mechanisms involving spatial or temporal environmental heterogeneity, resource partitioning and natural enemies. We first demonstrate the limitations of phenomenological approaches using simulations, and then (1) propose trait-based tests of coexistence, (2) generate hypotheses about which plant functional traits are likely to interact with particular mechanisms and (3) review the literature for evidence for these hypotheses. Theory and data suggest that all four classes of coexistence mechanisms could act on functional trait variation, but some mechanisms will be stronger and more widespread than others. The highest priority for future research is studies of interactions between environmental heterogeneity and trait variation that measure environmental variables at within-community scales and quantify species' responses to the environment in the absence of competition. Evidence that similar trait-based coexistence mechanisms operate in many ecosystems would simplify biodiversity forecasting and represent a rare victory for generality over contingency in community ecology

Divergent drivers of leaf trait variation within species, among species, and among functional groups.

Understanding variation in leaf functional traits-including rates of photosynthesis and respiration and concentrations of nitrogen and phosphorus-is a fundamental challenge in plant ecophysiology. When expressed per unit leaf area, these traits typically increase with leaf mass per area (LMA) within species but are roughly independent of LMA across the global flora. LMA is determined by mass components with different biological functions, including photosynthetic mass that largely determines metabolic rates and contains most nitrogen and phosphorus, and structural mass that affects toughness and leaf lifespan (LL). A possible explanation for the contrasting trait relationships is that most LMA variation within species is associated with variation in photosynthetic mass, whereas most LMA variation across the global flora is associated with variation in structural mass. This hypothesis leads to the predictions that (i) gas exchange rates and nutrient concentrations per unit leaf area should increase strongly with LMA across species assemblages with low LL variance but should increase weakly with LMA across species assemblages with high LL variance and that (ii) controlling for LL variation should increase the strength of the above LMA relationships. We present analyses of intra- and interspecific trait variation from three tropical forest sites and interspecific analyses within functional groups in a global dataset that are consistent with the above predictions. Our analysis suggests that the qualitatively different trait relationships exhibited by different leaf assemblages can be understood by considering the degree to which photosynthetic and structural mass components contribute to LMA variation in a given assemblage

Flowering Date of Taxonomic Families Predicts Phenological Sensitivity to Temperature: Implications for Forecasting the Effects of Climate Change on Unstudied Taxa

Premise of the study: Numerous long-term studies in seasonal habitats have tracked interannual variation in fi rst fl owering date (FFD) in relation to climate, documenting the effect of warming on the FFD of many species. Despite these efforts, long-term phenological observations are still lacking for many species. If we could forecast responses based on taxonomic affi nity, however, then we could leverage existing data to predict the climate-related phenological shifts of many taxa not yet studied; Methods: We examined phenological time series of 1226 species occurrences (1031 unique species in 119 families) across seven sites in North America and England to determine whether family membership (or family mean FFD) predicts the sensitivity of FFD to standardized interannual changes in temperature and precipitation during seasonal periods before fl owering and whether families differ signifi cantly in the direction of their phenological shifts; Key results: Patterns observed among species within and across sites are mirrored among family means across sites; earlyfl owering families advance their FFD in response to warming more than late-fl owering families. By contrast, we found no consistent relationships among taxa between mean FFD and sensitivity to precipitation as measured here; Conclusions: Family membership can be used to identify taxa of high and low sensitivity to temperature within the seasonal, temperate zone plant communities analyzed here. The high sensitivity of early-fl owering families (and the absence of earlyfl owering families not sensitive to temperature) may refl ect plasticity in fl owering time, which may be adaptive in environments where early-season conditions are highly variable among years

Intraspecific leaf trait variability along a boreal-to-tropical community diversity gradient

Disentangling the mechanisms that shape community assembly across diversity gradients is a central matter in ecology. While many studies have explored community assembly through species average trait values, there is a growing understanding that intraspecific trait variation (ITV) can also play a critical role in species coexistence. Classic biodiversity theory hypothesizes that higher diversity at species-rich sites can arise from narrower niches relative to species-poor sites, which would be reflected in reduced ITV as species richness increases. To explore how ITV in woody plant communities changes with species richness, we compiled leaf trait data (leaf size and specific leaf area) in a total of 521 woody plant species from 21 forest communities that differed dramatically in species richness, ranging from boreal to tropical rainforests. At each forest, we assessed ITV as an estimate of species niche breadth and we quantified the degree of trait overlap among co-occurring species as a measure of species functional similarity. We found ITV was relatively invariant across the species richness gradient. In addition, we found that species functional similarity increased with diversity. Contrary to the expectation from classic biodiversity theory, our results rather suggest that neutral processes or equalizing mechanisms can be acting as potential drivers shaping community assembly in hyperdiverse forests.Leaf data used to come to this manuscript resulted from many different studies involving different funding sources: European Union Seventh Framework Programme (FP7/2007-2013) under grant agreement no. 265171; Spanish-funded project REMEDINAL3eCM (S2013/MAE-2719); Ecometas excellence network (CGL2014-53840-REDT); a FPU grant by the Spanish Government (AP2010-5600); a Marie Curie IEF fellowship (FP7-PEOPLE-2011-IEF. no. 302445);Peer Reviewe

Divergent drivers of leaf trait variation within species, among species, and among functional groups

Understanding variation in leaf functional traits—including rates of photosynthesis and respiration and concentrations of nitrogen and phosphorus—is a fundamental challenge in plant ecophysiology. When expressed per unit leaf area, these traits typically increase with leaf mass per area (LMA) within species but are roughly independent of LMA across the global flora. LMA is determined by mass components with different biological functions, including photosynthetic mass that largely determines metabolic rates and contains most nitrogen and phosphorus, and structural mass that affects toughness and leaf lifespan (LL). A possible explanation for the contrasting trait relationships is that most LMA variation within species is associated with variation in photosynthetic mass, whereas most LMA variation across the global flora is associated with variation in structural mass. This hypothesis leads to the predictions that (i) gas exchange rates and nutrient concentrations per unit leaf area should increase strongly with LMA across species assemblages with low LL variance but should increase weakly with LMA across species assemblages with high LL variance and that (ii) controlling for LL variation should increase the strength of the above LMA relationships. We present analyses of intra- and interspecific trait variation from three tropical forest sites and interspecific analyses within functional groups in a global dataset that are consistent with the above predictions. Our analysis suggests that the qualitatively different trait relationships exhibited by different leaf assemblages can be understood by considering the degree to which photosynthetic and structural mass components contribute to LMA variation in a given assemblage

Tree mortality across biomes is promoted by drought intensity, lower wood density and higher specific leaf area

Drought events are increasing globally, and reports of consequent forest mortality are widespread. However, due to a lack of a quantitative global synthesis, it is still not clear whether drought-induced mortality rates differ among global biomes and whether functional traits influence the risk of drought-induced mortality. To address these uncertainties, we performed a global meta-analysis of 58 studies of drought-induced forest mortality. Mortality rates were modelled as a function of drought, temperature, biomes, phylogenetic and functional groups, and functional traits. We identified a consistent global-scale response, where mortality increased with drought severity (log mortality (trees trees-1 year-1) increased 0.46 (95% CI=0.2-0.7) with one SPEI unit drought intensity). We found no significant differences in the magnitude of the response depending on forest biomes or between angiosperms and gymnosperms or evergreen and deciduous tree species. Functional traits explained some of the variation in drought responses between species (i.e. increased from 30 to 37% when wood density and specific leaf area were included). Tree species with denser wood and lower specific leaf area showed lower mortality responses. Our results illustrate the value of functional traits for understanding patterns of drought-induced tree mortality and suggest that mortality could become increasingly widespread in the future

Megafauna extinction, tree species range reduction, and carbon storage in Amazonian forests

During the Late Pleistocene and early Holocene 59 species of South American megafauna went extinct. Their extinction potentially triggered population declines of large-seeded tree species dispersed by the large-bodied frugivores with which they co-evolved, a theory first proposed by Janzen and Martin (1982). We tested this hypothesis using species range maps for 257 South American tree species, comparing 63 species thought to be primarily distributed by megafauna with 194 distributed by other animals. We found a highly significant (p 95% following disperser extinction. A numerical gap dynamic simulations suggests that over a 10 000 yr period following the disperser extinctions, the average convex hull range size of large-seeded tree species decreased by ∼ 31%, while the estimated decrease in population size was ∼ 54%, indicating a likely greater decrease in species population size than indicated by the empirical range patterns. Finally, we found a positive correlation between seed size and wood density of animal-dispersed tree species implying that the Late Pleistocene and early Holocene megafaunal extinctions reduced carbon content in the Amazon by ∼ 1.5 ± 0.7%. In conclusion, we 1) provide some empirical evidence that megafauna distributed fruit species have a smaller mean range size than wind, water or other animal-dispersed species, 2) demonstrate mathematically that such range reductions are expected from megafauna extinctions ca 12 000 yr ago, and 3) illustrate that these extinctions may have reduced the Amazon's carbon storage capacity

The ecological causes of functional distinctiveness in communities

Recent work has shown that evaluating functional trait distinctiveness, the average trait distance of a species to other species in a community offers promising insights into biodiversity dynamics and ecosystem functioning. However, the ecological mechanisms underlying the emergence and persistence of functionally distinct species are poorly understood. Here, we address the issue by considering a heterogeneous fitness landscape whereby functional dimensions encompass peaks representing trait combinations yielding positive population growth rates in a community. We identify four ecological cases contributing to the emergence and persistence of functionally distinct species. First, environmental heterogeneity or alternative phenotypic designs can drive positive population growth of functionally distinct species. Second, sink populations with negative population growth can deviate from local fitness peaks and be functionally distinct. Third, species found at the margin of the fitness landscape can persist but be functionally distinct. Fourth, biotic interactions (positive or negative) can dynamically alter the fitness landscape. We offer examples of these four cases and guidelines to distinguish between them. In addition to these deterministic processes, we explore how stochastic dispersal limitation can yield functional distinctiveness. Our framework offers a novel perspective on the relationship between fitness landscape heterogeneity and the functional composition of ecological assemblages