Detecting Nonadditive Biotic Interactions and Assessing Their Biological Relevance among Temperate Rainforest Trees.

AbstractInteractions between and within abiotic and biotic processes generate nonadditive density-dependent effects on species performance that can vary in strength or direction across environments. If ignored, nonadditivities can lead to inaccurate predictions of species responses to environmental and compositional changes. While there are increasing empirical efforts to test the constancy of pairwise biotic interactions along environmental and compositional gradients, few assess both simultaneously. Using a nationwide forest inventory that spans broad ambient temperature and moisture gradients throughout New Zealand, we address this gap by analyzing the diameter growth of six focal tree species as a function of neighbor densities and climate, as well as neighbor × climate and neighbor × neighbor statistical interactions. The most complex model featuring all interaction terms had the highest predictive accuracy. Compared with climate variables, biotic interactions typically had stronger effects on diameter growth, especially when subjected to nonadditivities from local climatic conditions and the density of intermediary species. Furthermore, statistically strong (or weak) nonadditivities could be biologically irrelevant (or significant) depending on whether a species pair typically interacted under average or more extreme conditions. Our study highlights the importance of considering both the statistical potential and the biological relevance of nonadditive biotic interactions when assessing species performance under global change.

Functional trait trade-offs define plant population stability across different biomes.

Ecological theory posits that temporal stability patterns in plant populations are associated with differences in species' ecological strategies. However, empirical evidence is lacking about which traits, or trade-offs, underlie species stability, especially across different biomes. We compiled a worldwide collection of long-term permanent vegetation records (greater than 7000 plots from 78 datasets) from a large range of habitats which we combined with existing trait databases. We tested whether the observed inter-annual variability in species abundance (coefficient of variation) was related to multiple individual traits. We found that populations with greater leaf dry matter content and seed mass were more stable over time. Despite the variability explained by these traits being low, their effect was consistent across different datasets. Other traits played a significant, albeit weaker, role in species stability, and the inclusion of multi-variate axes or phylogeny did not substantially modify nor improve predictions. These results provide empirical evidence and highlight the relevance of specific ecological trade-offs, i.e. in different resource-use and dispersal strategies, for plant populations stability across multiple biomes. Further research is, however, necessary to integrate and evaluate the role of other specific traits, often not available in databases, and intraspecific trait variability in modulating species stability.

A function-based typology for Earth's ecosystems.

As the United Nations develops a post-2020 global biodiversity framework for the Convention on Biological Diversity, attention is focusing on how new goals and targets for ecosystem conservation might serve its vision of 'living in harmony with nature'1,2. Advancing dual imperatives to conserve biodiversity and sustain ecosystem services requires reliable and resilient generalizations and predictions about ecosystem responses to environmental change and management3. Ecosystems vary in their biota4, service provision5 and relative exposure to risks6, yet there is no globally consistent classification of ecosystems that reflects functional responses to change and management. This hampers progress on developing conservation targets and sustainability goals. Here we present the International Union for Conservation of Nature (IUCN) Global Ecosystem Typology, a conceptually robust, scalable, spatially explicit approach for generalizations and predictions about functions, biota, risks and management remedies across the entire biosphere. The outcome of a major cross-disciplinary collaboration, this novel framework places all of Earth's ecosystems into a unifying theoretical context to guide the transformation of ecosystem policy and management from global to local scales. This new information infrastructure will support knowledge transfer for ecosystem-specific management and restoration, globally standardized ecosystem risk assessments, natural capital accounting and progress on the post-2020 global biodiversity framework.

Global patterns of vascular plant alpha diversity.

Global patterns of regional (gamma) plant diversity are relatively well known, but whether these patterns hold for local communities, and the dependence on spatial grain, remain controversial. Using data on 170,272 georeferenced local plant assemblages, we created global maps of alpha diversity (local species richness) for vascular plants at three different spatial grains, for forests and non-forests. We show that alpha diversity is consistently high across grains in some regions (for example, Andean-Amazonian foothills), but regional 'scaling anomalies' (deviations from the positive correlation) exist elsewhere, particularly in Eurasian temperate forests with disproportionally higher fine-grained richness and many African tropical forests with disproportionally higher coarse-grained richness. The influence of different climatic, topographic and biogeographical variables on alpha diversity also varies across grains. Our multi-grain maps return a nuanced understanding of vascular plant biodiversity patterns that complements classic maps of biodiversity hotspots and will improve predictions of global change effects on biodiversity.

Co-limitation towards lower latitudes shapes global forest diversity gradients.

The latitudinal diversity gradient (LDG) is one of the most recognized global patterns of species richness exhibited across a wide range of taxa. Numerous hypotheses have been proposed in the past two centuries to explain LDG, but rigorous tests of the drivers of LDGs have been limited by a lack of high-quality global species richness data. Here we produce a high-resolution (0.025° × 0.025°) map of local tree species richness using a global forest inventory database with individual tree information and local biophysical characteristics from ~1.3 million sample plots. We then quantify drivers of local tree species richness patterns across latitudes. Generally, annual mean temperature was a dominant predictor of tree species richness, which is most consistent with the metabolic theory of biodiversity (MTB). However, MTB underestimated LDG in the tropics, where high species richness was also moderated by topographic, soil and anthropogenic factors operating at local scales. Given that local landscape variables operate synergistically with bioclimatic factors in shaping the global LDG pattern, we suggest that MTB be extended to account for co-limitation by subordinate drivers.

The adaptive challenge of extreme conditions shapes evolutionary diversity of plant assemblages at continental scales.

The tropical conservatism hypothesis (TCH) posits that the latitudinal gradient in biological diversity arises because most extant clades of animals and plants originated when tropical environments were more widespread and because the colonization of colder and more seasonal temperate environments is limited by the phylogenetically conserved environmental tolerances of these tropical clades. Recent studies have claimed support of the TCH, indicating that temperate plant diversity stems from a few more recently derived lineages that are nested within tropical clades, with the colonization of the temperate zone being associated with key adaptations to survive colder temperatures and regular freezing. Drought, however, is an additional physiological stress that could shape diversity gradients. Here, we evaluate patterns of evolutionary diversity in plant assemblages spanning the full extent of climatic gradients in North and South America. We find that in both hemispheres, extratropical dry biomes house the lowest evolutionary diversity, while tropical moist forests and many temperate mixed forests harbor the highest. Together, our results support a more nuanced view of the TCH, with environments that are radically different from the ancestral niche of angiosperms having limited, phylogenetically clustered diversity relative to environments that show lower levels of deviation from this niche. Thus, we argue that ongoing expansion of arid environments is likely to entail higher loss of evolutionary diversity not just in the wet tropics but in many extratropical moist regions as well.

Synchrony matters more than species richness in plant community stability at a global scale.

The stability of ecological communities is critical for the stable provisioning of ecosystem services, such as food and forage production, carbon sequestration, and soil fertility. Greater biodiversity is expected to enhance stability across years by decreasing synchrony among species, but the drivers of stability in nature remain poorly resolved. Our analysis of time series from 79 datasets across the world showed that stability was associated more strongly with the degree of synchrony among dominant species than with species richness. The relatively weak influence of species richness is consistent with theory predicting that the effect of richness on stability weakens when synchrony is higher than expected under random fluctuations, which was the case in most communities. Land management, nutrient addition, and climate change treatments had relatively weak and varying effects on stability, modifying how species richness, synchrony, and stability interact. Our results demonstrate the prevalence of biotic drivers on ecosystem stability, with the potential for environmental drivers to alter the intricate relationship among richness, synchrony, and stability.

Macroclimate and Topography Interact to Influence the Abundance of Divaricate Plants in New Zealand.

The abundance of the divaricate growth form in New Zealand has been interpreted as either (a) the response of an isolated flora to cool, dry, Plio-Pleistocene climates; or (b) a defense against large browsing birds (moa) that were hunted to extinction shortly after human arrival during the last millennium. We used patterns of divaricate plant abundance across present-day landscapes to test a novel synthetic hypothesis: that the divaricate form is of most value to plants on fertile soils that attract herbivores, on sites where climatic constraints prevent plants from quickly growing out of the browse zone. This hypothesis predicts that divaricate species should be most abundant on terraces (landforms that are both fertile and frost-prone) in regions that are cold and dry, and should be scarce across all topographic positions in the warmest (largely frost-free) regions. To address our hypothesis, we first tested the influence of topography on frost regimes and nutrient levels by measuring temperatures and soil total C, N, and P at four standard topographic positions at five localities differing widely in macroclimate. We then extracted a dataset of 236 surveys comprising 9,877 relevé plots from the New Zealand National Vegetation Survey databank. We calculated the proportion of arborescent species with a divaricate growth form and the proportion of total arborescentcover contributed by divaricates on each plot; we then fitted linear mixed-effect models predicting these response variables as functions of topographic position and climate. The number of frosts recorded averaged <1 yr-1 at the warmest of the five sites studied, to >60 yr-1 on all topographic positions at the coldest site. Terraces were subject to more frequent and harder frosts than any other topographic position. Topography had no significant influence on total N or C:N, but total P was higher on terraces and in gullies than on faces or ridges. Frost-free period was the dominant influence on both species representation and cover of divaricate plants throughout the country. The effect of topography was also significant, but weaker. The effect of frost-free period was stronger on sites with water deficits than on sites where precipitation exceeded evapotranspiration. Divaricates made their largest contributions on terraces in cold, dry regions; as predicted, they were scarce on all topographic positions on sites with frost-free periods >300 days. Our hypothesis was generally supported, although the effect of topography on divaricate abundance was not as strong as some previous studies led us to expect. Divaricates made their largest contributions to arborescent species richness and cover on sites where climatic restrictions on growth coincide with relatively high nutrient availability. The contemporary distribution of the divaricate form across New Zealand landscapes thus appears to be reasonably well explained by the hypothesized interaction of climate and fertility-mediated browsing, although experiments may provide more conclusive tests of this hypothesis.

The commonness of rarity: Global and future distribution of rarity across land plants.

A key feature of life's diversity is that some species are common but many more are rare. Nonetheless, at global scales, we do not know what fraction of biodiversity consists of rare species. Here, we present the largest compilation of global plant diversity to quantify the fraction of Earth's plant biodiversity that are rare. A large fraction, ~36.5% of Earth's ~435,000 plant species, are exceedingly rare. Sampling biases and prominent models, such as neutral theory and the k-niche model, cannot account for the observed prevalence of rarity. Our results indicate that (i) climatically more stable regions have harbored rare species and hence a large fraction of Earth's plant species via reduced extinction risk but that (ii) climate change and human land use are now disproportionately impacting rare species. Estimates of global species abundance distributions have important implications for risk assessments and conservation planning in this era of rapid global change.

Functional trait space and the latitudinal diversity gradient.

The processes causing the latitudinal gradient in species richness remain elusive. Ecological theories for the origin of biodiversity gradients, such as competitive exclusion, neutral dynamics, and environmental filtering, make predictions for how functional diversity should vary at the alpha (within local assemblages), beta (among assemblages), and gamma (regional pool) scales. We test these predictions by quantifying hypervolumes constructed from functional traits representing major axes of plant strategy variation (specific leaf area, plant height, and seed mass) in tree assemblages spanning the temperate and tropical New World. Alpha-scale trait volume decreases with absolute latitude and is often lower than sampling expectation, consistent with environmental filtering theory. Beta-scale overlap decays with geographic distance fastest in the temperate zone, again consistent with environmental filtering theory. In contrast, gamma-scale trait space shows a hump-shaped relationship with absolute latitude, consistent with no theory. Furthermore, the overall temperate trait hypervolume was larger than the overall tropical hypervolume, indicating that the temperate zone permits a wider range of trait combinations or that niche packing is stronger in the tropical zone. Although there are limitations in the data, our analyses suggest that multiple processes have shaped trait diversity in trees, reflecting no consistent support for any one theory.

Quantifying invasion resistance: the use of recruitment functions to control for propagule pressure.

Invasive species distributions tend to be biased towards some habitats compared to others due to the combined effects of habitat-specific resistance to invasion and non-uniform propagule pressure. These two factors may also interact, with habitat resistance varying as a function of propagule supply rate. Recruitment experiments, in which the number of individuals recruiting into a population is measured under different propagule supply rates, can help us understand these interactions and quantify habitat resistance to invasion while controlling for variation in propagule supply rate. Here, we constructed recruitment functions for the invasive herb Hieracium lepidulum by sowing seeds at five different densities into six different habitat types in New Zealand's Southern Alps repeated over two successive years, and monitored seedling recruitment and survival over a four year period. We fitted recruitment functions that allowed us to estimate the total number of safe sites available for plants to occupy, which we used as a measure of invasion resistance, and tested several hypotheses concerning how invasion resistance differed among habitats and over time. We found significant differences in levels of H. lepidulum recruitment among habitats, which did not match the species' current distribution in the landscape. Local biotic and abiotic characteristics helped explain some of the between-habitat variation, with vascular plant species richness, vascular plant cover, and light availability, all positively correlated with the number of safe sites for recruitment. Resistance also varied over time however, with cohorts sown in successive years showing different levels of recruitment in some habitats but not others. These results show that recruitment functions can be used to quantify habitat resistance to invasion and to identify potential mechanisms of invasion resistance.

Habitat area and climate stability determine geographical variation in plant species range sizes.

Despite being a fundamental aspect of biodiversity, little is known about what controls species range sizes. This is especially the case for hyperdiverse organisms such as plants. We use the largest botanical data set assembled to date to quantify geographical variation in range size for ~ 85 000 plant species across the New World. We assess prominent hypothesised range-size controls, finding that plant range sizes are codetermined by habitat area and long- and short-term climate stability. Strong short- and long-term climate instability in large parts of North America, including past glaciations, are associated with broad-ranged species. In contrast, small habitat areas and a stable climate characterise areas with high concentrations of small-ranged species in the Andes, Central America and the Brazilian Atlantic Rainforest region. The joint roles of area and climate stability strengthen concerns over the potential effects of future climate change and habitat loss on biodiversity.

Functional traits reveal processes driving natural afforestation at large spatial scales.

An understanding of the processes governing natural afforestation over large spatial scales is vital for enhancing forest carbon sequestration. Models of tree species occurrence probability in non-forest vegetation could potentially identify the primary variables determining natural afforestation. However, inferring processes governing afforestation using tree species occurrence is potentially problematic, since it is impossible to know whether observed occurrences are due to recruitment or persistence of existing trees following disturbance. Plant functional traits have the potential to reveal the processes by which key environmental and land cover variables influence afforestation. We used 10,061 survey plots to identify the primary environmental and land cover variables influencing tree occurrence probability in non-forest vegetation in New Zealand. We also examined how these variables influenced diversity of functional traits linked to plant ecological strategy and dispersal ability. Mean annual temperature was the most important environmental predictor of tree occurrence. Local woody cover and distance to forest were the most important land cover variables. Relationships between these variables and ecological strategy traits revealed a trade-off between ability to compete for light and colonize sites that were marginal for tree occurrence. Biotically dispersed species occurred less frequently with declining temperature and local woody cover, suggesting that abiotic stress limited their establishment and that biotic dispersal did not increase ability to colonize non-woody vegetation. Functional diversity for ecological strategy traits declined with declining temperature and woody cover and increasing distance to forest. Functional diversity for dispersal traits showed the opposite trend. This suggests that low temperatures and woody cover and high distance to forest may limit tree species establishment through filtering on ecological strategy traits, but not on dispersal traits. This study shows that 'snapshot' survey plot data, combined with functional trait data, may reveal the processes driving tree species establishment in non-forest vegetation over large spatial scales.

Status assessment of New Zealand's naturally uncommon ecosystems.

Globally, ecosystems are under increasing anthropogenic pressure; thus, many are at risk of elimination. This situation has led the International Union for Conservation of Nature (IUCN) to propose a quantitative approach to ecosystem-risk assessment. However, there is a need for their proposed criteria to be evaluated through practical examples spanning a diverse range of ecosystems and scales. We applied the IUCN's ecosystem red-list criteria, which are based on changes in extent of ecosystems and reductions in ecosystem processes, to New Zealand's 72 naturally uncommon ecosystems. We aimed to test the applicability of the proposed criteria to ecosystems that are naturally uncommon (i.e., those that would naturally occur over a small area in the absence of human activity) and to provide information on the probability of ecosystem elimination so that conservation priorities might be set. We also tested the hypothesis that naturally uncommon ecosystems classified as threatened on the basis of IUCN Red List criteria contain more threatened plant species than those classified as nonthreatened. We identified 18 critically endangered, 17 endangered, and 10 vulnerable ecosystems. We estimated that naturally uncommon ecosystems contained 145 (85%) of mainland New Zealand's taxonomically distinct nationally critical, nationally endangered, and nationally vulnerable plant species, 66 (46%) of which were endemic to naturally uncommon ecosystems. We estimated there was a greater number of threatened plant species (per unit area) in critically endangered ecosystems than in ecosystems classified as nonthreatened. With their high levels of endemism and rapid and relatively well-documented history of anthropogenic change, New Zealand's naturally uncommon ecosystems provide an excellent case-study for the ongoing development of international criteria for threatened ecosystems. We suggest that interactions and synergies among decline in area, decline in function, and the scale of application of the criteria be used to improve the IUCN criteria for threatened ecosystems.

Disturbance affects short-term facilitation, but not long-term saturation, of exotic plant invasion in New Zealand forest.

We investigate the spread of an exotic herb, Hieracium lepidulum, into a New Zealand Nothofagus forest with the aim of understanding how stand-development of tree populations, propagule pressure and invader persistence, affect invasion across the landscape and within communities. Using data repeatedly collected over 35 years, from 250 locations, we parametrize continuous-time Markov chain models and use these models to examine future projections of the invasion under a range of hypothetical scenarios. We found that the probability of invasion into a stand was relatively high following canopy disturbance and that local abundance of Hieracium was promoted by minor disturbances. However, model predictions extrapolated 45 years into the future show that neither the rate of landscape-level invasion, nor local population growth of Hieracium, was affected much by changing the frequency of canopy disturbance events. Instead, invasion levels were strongly affected by the ability of Hieracium to persist in the understorey following forest canopy closure, and by propagule supply from streams, forest edges and plants already established within the stand. Our results show that disturbance frequency has surprisingly little influence on the long-term trajectory of invasion, while invader persistence strongly determines invasion patterns.

Aboveground and belowground effects of single-tree removals in New Zealand rain forest.

There has been considerable recent interest in how human-induced species loss affects community and ecosystem properties. These effects are particularly apparent when a commercially valuable species is harvested from an ecosystem, such as occurs through single-tree harvesting or selective logging of desired timber species in natural forests. In New Zealand mixed-species rain forests, single-tree harvesting of the emergent gymnosperm Dacrydium cupressinum, or rimu, has been widespread. This harvesting has been contentious in part because of possible ecological impacts of Dacrydium removal on the remainder of the forest, but many of these effects remain unexplored. We identified an area where an unintended 40-year "removal experiment" had been set up that involved selective extraction of individual Dacrydium trees. We measured aboveground and belowground variables at set distances from both individual live trees and stumps of trees harvested 40 years ago. Live trees had effects both above and below ground by affecting diversity and cover of several components of the vegetation (usually negatively), promoting soil C sequestration, enhancing ratios of soil C:P and N:P, and affecting community structure of soil microflora. These effects extended to 8 m from the tree base and were likely caused by poor-quality litter and humus produced by the trees. Measurements for the stumps revealed strong legacy effects of prior presence of trees on some properties (e.g., cover by understory herbs and ferns, soil C sequestration, soil C:P and N:P ratios), but not others (e.g., soil fungal biomass, soil N concentration). These results suggest that the legacy of prior presence of Dacrydium may remain for several decades or centuries, and certainly well over 40 years. They also demonstrate that, while large Dacrydium individuals (and their removal) may have important effects in their immediate proximity, within a forest, these effects should only be important in localized patches containing high densities of large trees. Finally, this study emphasizes that deliberate extraction of a particular tree species from a forest can exert influences both above and below ground if the removed species has a different functional role than that of the other plant species present.

Context matters: matrix vegetation influences native and exotic species composition on habitat islands.

The extensive research on plant communities of natural-habitat islands has primarily focused on the "islands." The island analogy, however, potentially limits understanding of processes influencing composition on habitat islands because the nature of their matrix is overlooked. We determine how plant community structure of the surrounding matrix influences vegetation on volcanic outcrops in the modified landscape of Banks Peninsula, New Zealand. Our primary purpose is to address whether the matrix is more important for recently established exotic species than it is for well-established native species and whether such invasion by exotics has led to homogenization of the outcrop flora. To test this, we examined our data at three spatial scales: that of the entire outcrop flora, between individual outcrops and their immediate surrounding matrix, and between individual outcrop faces and the individual relevés of the immediate surrounding matrix. We found that 81% of the native flora and 90% of the exotic flora also occur in the matrix. This high level of species shared between the outcrop and matrix persists at the scale of individual outcrop faces (68% of the total flora of individual faces is shared with the matrix). We predicted that floras from different outcrops would vary in their distinctiveness from their immediate matrix. We found Bray-Curtis distance coefficient values to range from 0.26 to 0.64; these were even higher at the outcrop-face scale. Variability in outcrop distinctiveness relates primarily to the outcrop face properties of area, vegetation height, and soil depth, and matrix properties of vegetation structure and vegetation heterogeneity. The effect of the vegetation structure of the matrix is more pronounced on the exotic than on the native outcrop flora. The component of composition and structure of the matrix that was independent of outcrop properties and local environment accounts for similar levels of explainable variation in total and native composition (29-31%), but considerably more (40%) in composition of exotic species. Our results support our prediction that, as the surrounding matrix becomes more modified, invasion by exotics makes outcrop vegetation less distinct from its matrix.

Growth-size scaling relationships of woody plant species differ from predictions of the Metabolic Ecology Model.

The Metabolic Ecology Model predicts that tree diameter (D) growth (dD/dt) scales with D(1/3). Using data on diameter growth and height-diameter relationships for 56 and 40 woody species, respectively, from forests throughout New Zealand, we tested one prediction and two assumptions of this model: (i) the exponent of the growth-diameter scaling relationship equals 1/3 and is invariant among species and growth forms, (ii) small and large individuals are invariant in their exponents and (iii) tree height scales with D(2/3). We found virtually no support for any prediction or assumption: growth-diameter scaling exponents varied substantially among species and growth forms, correlated positively with species' maximum height, and shifted significantly with increasing individual size. Tree height did not scale invariantly with diameter. Based on a quantitative test, violation of these assumptions alone could not explain the model's poor fit to our data, possibly reflecting multiple, unsound assumptions, as well as unaccounted-for variation that should be incorporated.