The emergent interactions that govern biodiversity change

Biodiversity in changing environments: An external-driver internal-topology framework to guide intervention

Accompanying the climate crisis is the more enigmatic biodiversity crisis. Rapid reorganization of biodiversity due to global environmental change has defied prediction and tested the basic tenets of conservation and restoration. Conceptual and practical innovation is needed to support decision making in the face of these unprecedented shifts. Critical questions include: How can we generalize biodiversity change at the community level? When are systems able to reorganize and maintain integrity, and when does abiotic change result in collapse or restructuring? How does this understanding provide a template to guide when and how to intervene in conservation and restoration? To this end, we frame changes in community organization as the modulation of external abiotic drivers on the internal topology of species interactions, using plant-plant interactions in terrestrial communities as a starting point. We then explore how this framing can help translate available data on species abundance and trait distributions to corresponding decisions in management. Given the expectation that community response and reorganization are highly complex, the external-driver internal-topology (EDIT) framework offers a way to capture general patterns of biodiversity that can help guide resilience and adaptation in changing environments.

Biodiversity modeling advances will improve predictions of nature's contributions to people

Accurate predictions of ecosystem functions and nature's contributions to people (NCP) are needed to prioritize environmental protection and restoration in the Anthropocene. However, our ability to predict NCP is undermined by approaches that rely on biophysical variables and ignore those describing biodiversity, which have strong links to NCP. To foster predictive mapping of NCP, we should harness the latest methods in biodiversity modeling. This field advances rapidly, and new techniques with promising applications for predicting NCP are still underutilized. Here, we argue that employing recent advances in biodiversity modeling can enhance the accuracy and scope of NCP maps and predictions. This enhancement will contribute significantly to the achievement of global objectives to preserve NCP, for both the present and an unpredictable future.

Seasonality and longer-term development generate temporal dynamics in the Populus microbiome

Temporal variation in community composition is central to our understanding of the assembly and functioning of microbial communities, yet the controls over temporal dynamics for microbiomes of long-lived plants, such as trees, remain unclear. Temporal variation in tree microbiomes could arise primarily from seasonal (i.e., intra-annual) fluctuations in community composition or from longer-term changes across years as host plants age. To test these alternatives, we experimentally isolated temporal variation in plant microbiome composition using a common garden and clonally propagated plants, and we used amplicon sequencing to characterize bacterial/archaeal and fungal communities in the leaf endosphere, root endosphere, and rhizosphere of two Populus spp. over four seasons across two consecutive years. Microbial community composition differed among seasons and years (which accounted for up to 21% of the variation in microbial community composition) and was correlated with seasonal dissimilarity in climatic conditions. However, microbial community dissimilarity was also positively correlated with time, reflecting longer-term compositional shifts as host trees aged. Together, our findings demonstrate that temporal patterns in tree microbiomes arise from both seasonal fluctuations and longer-term changes, which interact to generate unique seasonal patterns each year. In addition to shedding light on two important controls over the assembly of plant microbiomes, our results also suggest future studies of tree microbiomes should account for background temporal dynamics when testing the drivers of spatial patterns in microbial community composition and temporal responses of plant microbiomes to environmental change.IMPORTANCEMicrobiomes are integral to the health of host plants, but we have a limited understanding of the factors that control how the composition of plant microbiomes changes over time. Especially little is known about the microbiome of long-lived trees, relative to annual and non-woody plants. We tested how tree microbiomes changed between seasons and years in poplar (genus Populus), which are widespread and ecologically important tree species that also serve as important biofuel feedstocks. We found the composition of bacterial, archaeal, and fungal communities differed among seasons, but these seasonal differences depended on year. This dependence was driven by longer-term changes in microbial composition as host trees developed across consecutive years. Our findings suggest that temporal variation in tree microbiomes is driven by both seasonal fluctuations and longer-term (i.e., multiyear) development.

Joint species distribution modeling reveals a changing prey landscape for North Pacific right whales on the Bering shelf

The eastern North Pacific right whale (NPRW) is the most endangered population of whale and has been observed north of its core feeding ground in recent years with low sea ice extent. Sea ice and water temperature are important drivers for zooplankton dynamics within the whale's core feeding ground in the southeastern Bering Sea, seasonally forming stable fronts along the shelf that give rise to distinct zooplankton communities. A northward shift in NPRW distribution driven by changing distribution of prey resources could put this species at increased risk of entanglement and vessel strikes. We modeled the abundance of NPRW prey, Calanus glacialis, Neocalanus, and Thysanoessa species, using a dynamic biophysical food web model of nine zooplankton guilds in the Bering shelf zooplankton community during a period of warming (2006-2016). This model is unique among prior zooplankton studies from the region in that it includes density dependence, thereby allowing us to ask whether species interactions influence zooplankton dynamics. Modeling confirmed the importance of sea ice and ocean temperature to zooplankton dynamics in the region. Density-independent growth drove community dynamics, while dependent factors were comparatively minimal. Overall, Calanus responded to environment terms, with the strength and direction of response driven by copepodite stage. Neocalanus and Thysanoessa responses were weaker, likely due to their primary occurrence on the outer shelf. We also modeled the steady-state (equilibrium) abundance of Calanus in conditions with and without wind gusts to test whether advection of outer shelf species might disrupt the steady-state dynamics of Calanus abundance; the results did not support disruption. Given the annual fall sampling design, we interpret our results as follows: low-ice-extent winters induced stronger spring winds and weakened fronts on the shelf, thereby advecting some outer shelf species into the study region; increased development rates in these warm conditions influenced the proportion of C. glacialis copepodite stages over the season. Residual correlation suggests missing drivers, possibly predators, and phytoplankton bloom composition. Given the continued loss of sea ice in the region and projected continued warming, our findings suggest that C. glacialis will move northward, and thus, whales may move northward to continue targeting them.

Mechanistic modeling of climate effects on redistribution and population growth in a community of fish species

Understanding community responses to climate is critical for anticipating the future impacts of global change. However, despite increased research efforts in this field, models that explicitly include important biological mechanisms are lacking. Quantifying the potential impacts of climate change on species is complicated by the fact that the effects of climate variation may manifest at several points in the biological process. To this end, we extend a dynamic mechanistic model that combines population dynamics, such as species interactions, with species redistribution by allowing climate to affect both processes. We examine their relative contributions in an application to the changing biomass of a community of eight species in the Gulf of Maine using over 30 years of fisheries data from the Northeast Fishery Science Center. Our model suggests that the mechanisms driving biomass trends vary across space, time, and species. Phase space plots demonstrate that failing to account for the dynamic nature of the environmental and biologic system can yield theoretical estimates of population abundances that are not observed in empirical data. The stock assessments used by fisheries managers to set fishing targets and allocate quotas often ignore environmental effects. At the same time, research examining the effects of climate change on fish has largely focused on redistribution. Frameworks that combine multiple biological reactions to climate change are particularly necessary for marine researchers. This work is just one approach to modeling the complexity of natural systems and highlights the need to incorporate multiple and possibly interacting biological processes in future models.

Which demographic processes control competitive equilibria? Bayesian calibration of a size-structured forest population model

In forest communities, light competition is a key process for community assembly. Species' differences in seedling and sapling tolerance to shade cast by overstory trees is thought to determine species composition at late-successional stages. Most forests are distant from these late-successional equilibria, impeding a formal evaluation of their potential species composition. To extrapolate competitive equilibria from short-term data, we therefore introduce the JAB model, a parsimonious dynamic model with interacting size-structured populations, which focuses on sapling demography including the tolerance to overstory competition. We apply the JAB model to a two-"species" system from temperate European forests, that is, the shade-tolerant species Fagus sylvatica L. and the group of all other competing species. Using Bayesian calibration with prior information from external Slovakian national forest inventory (NFI) data, we fit the JAB model to short time series from the German NFI. We use the posterior estimates of demographic rates to extrapolate that F. sylvatica will be the predominant species in 94% of the competitive equilibria, despite only predominating in 24% of the initial states. We further simulate counterfactual equilibria with parameters switched between species to assess the role of different demographic processes for competitive equilibria. These simulations confirm the hypothesis that the higher shade tolerance of F. sylvatica saplings is key for its long-term predominance. Our results highlight the importance of demographic differences in early life stages for tree species assembly in forest communities.

Species traits and observer behaviors that bias data assimilation and how to accommodate them

Datasets that monitor biodiversity capture information differently depending on their design, which influences observer behavior and can lead to biases across observations and species. Combining different datasets can improve our ability to identify and understand threats to biodiversity, but this requires an understanding of the observation bias in each. Two datasets widely used to monitor bird populations exemplify these general concerns: eBird is a citizen science project with high spatiotemporal resolution but variation in distribution, effort, and observers, whereas the Breeding Bird Survey (BBS) is a structured survey of specific locations over time. Analyses using these two datasets can identify contradictory population trends. To understand these discrepancies and facilitate data fusion, we quantify species-level reporting differences across eBird and the BBS in three regions across the United States by jointly modeling bird abundances using data from both datasets. First, we fit a joint Species Distribution Model that accounts for environmental conditions and effort to identify reporting differences across the datasets. We then examine how these differences in reporting are related to species traits. Finally, we analyze species reported to one dataset but not the other and determine whether traits differ between reported and unreported species. We find that most species are reported more in the BBS than eBird. Specifically, we find that compared to eBird, BBS observers tend to report higher counts of common species and species that are usually detected by sound. We also find that species associated with water are reported less in the BBS. Species typically identified by sound are reported more at sunrise than later in the morning. Our results quantify reporting differences in eBird and the BBS to enhance our understanding of how each captures information and how they should be used. The reporting rates we identify can also be incorporated into observation models through detectability or effort to improve analyses across species and datasets. The method demonstrated here can be used to compare reporting rates across any two or more datasets to examine biases.

Pathways to global-change effects on biodiversity: new opportunities for dynamically forecasting demography and species interactions

In structured populations, persistence under environmental change may be particularly threatened when abiotic factors simultaneously negatively affect survival and reproduction of several life cycle stages, as opposed to a single stage. Such effects can then be exacerbated when species interactions generate reciprocal feedbacks between the demographic rates of the different species. Despite the importance of such demographic feedbacks, forecasts that account for them are limited as individual-based data on interacting species are perceived to be essential for such mechanistic forecasting-but are rarely available. Here, we first review the current shortcomings in assessing demographic feedbacks in population and community dynamics. We then present an overview of advances in statistical tools that provide an opportunity to leverage population-level data on abundances of multiple species to infer stage-specific demography. Lastly, we showcase a state-of-the-art Bayesian method to infer and project stage-specific survival and reproduction for several interacting species in a Mediterranean shrub community. This case study shows that climate change threatens populations most strongly by changing the interaction effects of conspecific and heterospecific neighbours on both juvenile and adult survival. Thus, the repurposing of multi-species abundance data for mechanistic forecasting can substantially improve our understanding of emerging threats on biodiversity.

Integrating gradient with scale in ecological and evolutionary studies

Gradient and scale are two key concepts in ecology and evolution that are closely related but inherently distinct. While scale commonly refers to the dimensional space of a specific ecological/evolutionary (eco-evo) issue, gradient measures the range of a given variable. Gradient and scale can jointly and interactively influence eco-evo patterns. Extensive previous research investigated how changing scales may affect the observation and interpretation of eco-evo patterns; however, relatively little attention has been paid to the role of changing gradients. Here, synthesizing recent research progress, we suggest that the role of scale in the emergence of ecological patterns should be evaluated in conjunction with considering the underlying environmental gradients. This is important because, in most studies, the range of the gradient is often part of its full potential range. The difference between sampled (partial) versus potential (full) environmental gradients may profoundly impact observed eco-evo patterns and alter scale-gradient relationships. Based on observations from both field and experimental studies, we illustrate the underlying features of gradients and how they may affect observed patterns, along with the linkages of these features to scales. Since sampled gradients often do not cover their full potential ranges, we discuss how the breadth and the starting and ending positions of key gradients may affect research design and data interpretation. We then outline potential approaches and related perspectives to better integrate gradient with scale in future studies.

Dynamic Species Distribution Modeling Reveals the Pivotal Role of Human-Mediated Long-Distance Dispersal in Plant Invasion

Plant invasions generate massive ecological and economic costs worldwide. Predicting their spatial dynamics is crucial to the design of effective management strategies and the prevention of invasions. Earlier studies highlighted the crucial role of long-distance dispersal in explaining the speed of many invasions. In addition, invasion speed depends highly on the duration of its lag phase, which may depend on the scaling of fecundity with age, especially for woody plants, even though empirical proof is still rare. Bayesian dynamic species distribution models enable the fitting of process-based models to partial and heterogeneous observations using a state-space modeling approach, thus offering a tool to test such hypotheses on past invasions over large spatial scales. We use such a model to explore the roles of long-distance dispersal and age-structured fecundity in the transient invasion dynamics of Plectranthus barbatus, a woody plant invader in South Africa. Our lattice-based model accounts for both short and human-mediated long-distance dispersal, as well as age-structured fecundity. We fitted our model on opportunistic occurrences, accounting for the spatio-temporal variations of the sampling effort and the variable detection rates across datasets. The Bayesian framework enables us to integrate a priori knowledge on demographic parameters and control identifiability issues. The model revealed a massive wave of spatial spread driven by human-mediated long-distance dispersal during the first decade and a subsequent drastic population growth, leading to a global equilibrium in the mid-1990s. Without long-distance dispersal, the maximum population would have been equivalent to 30% of the current equilibrium population. We further identified the reproductive maturity at three years old, which contributed to the lag phase before the final wave of population growth. Our results highlighted the importance of the early eradication of weedy horticultural alien plants around urban areas to hamper and delay the invasive spread.

Global change re-structures alpine plant communities through interacting abiotic and biotic effects

Global change is altering patterns of community assembly, with net outcomes dependent on species' responses to the abiotic environment, both directly and mediated through biotic interactions. Here, we assess alpine plant community responses in a 15-year factorial nitrogen addition, warming and snow manipulation experiment. We used a dynamic competition model to estimate the density-dependent and -independent processes underlying changes in species-group abundances over time. Density-dependent shifts in competitive interactions drove long-term changes in abundance of species-groups under global change while counteracting environmental drivers limited the growth response of the dominant species through density-independent mechanisms. Furthermore, competitive interactions shifted with the environment, primarily with nitrogen and drove non-linear abundance responses across environmental gradients. Our results highlight that global change can either reshuffle species hierarchies or further favour already-dominant species; predicting which outcome will occur requires incorporating both density-dependent and -independent mechanisms and how they interact across multiple global change factors.

Climate variability and aridity modulate the role of leaf shelters for arthropods: A global experiment

Current climate change is disrupting biotic interactions and eroding biodiversity worldwide. However, species sensitive to aridity, high temperatures, and climate variability might find shelter in microclimatic refuges, such as leaf rolls built by arthropods. To explore how the importance of leaf shelters for terrestrial arthropods changes with latitude, elevation, and climate, we conducted a distributed experiment comparing arthropods in leaf rolls versus control leaves across 52 sites along an 11,790 km latitudinal gradient. We then probed the impact of short- versus long-term climatic impacts on roll use, by comparing the relative impact of conditions during the experiment versus average, baseline conditions at the site. Leaf shelters supported larger organisms and higher arthropod biomass and species diversity than non-rolled control leaves. However, the magnitude of the leaf rolls' effect differed between long- and short-term climate conditions, metrics (species richness, biomass, and body size), and trophic groups (predators vs. herbivores). The effect of leaf rolls on predator richness was influenced only by baseline climate, increasing in magnitude in regions experiencing increased long-term aridity, regardless of latitude, elevation, and weather during the experiment. This suggests that shelter use by predators may be innate, and thus, driven by natural selection. In contrast, the effect of leaf rolls on predator biomass and predator body size decreased with increasing temperature, and increased with increasing precipitation, respectively, during the experiment. The magnitude of shelter usage by herbivores increased with the abundance of predators and decreased with increasing temperature during the experiment. Taken together, these results highlight that leaf roll use may have both proximal and ultimate causes. Projected increases in climate variability and aridity are, therefore, likely to increase the importance of biotic refugia in mitigating the effects of climate change on species persistence.

Disentangling key species interactions in diverse and heterogeneous communities: A Bayesian sparse modelling approach

Modelling species interactions in diverse communities traditionally requires a prohibitively large number of species‐interaction coefficients, especially when considering environmental dependence of parameters. We implemented Bayesian variable selection via sparsity‐inducing priors on non‐linear species abundance models to determine which species interactions should be retained and which can be represented as an average heterospecific interaction term, reducing the number of model parameters. We evaluated model performance using simulated communities, computing out‐of‐sample predictive accuracy and parameter recovery across different input sample sizes. We applied our method to a diverse empirical community, allowing us to disentangle the direct role of environmental gradients on species’ intrinsic growth rates from indirect effects via competitive interactions. We also identified a few neighbouring species from the diverse community that had non‐generic interactions with our focal species. This sparse modelling approach facilitates exploration of species interactions in diverse communities while maintaining a manageable number of parameters.

Jointly modeling marine species to inform the effects of environmental change on an ecological community in the Northwest Atlantic

Single species distribution models (SSDMs) are typically used to understand and predict the distribution and abundance of marine fish by fitting distribution models for each species independently to a combination of abiotic environmental variables. However, species abundances and distributions are influenced by abiotic environmental preferences as well as biotic dependencies such as interspecific competition and predation. When species interact, a joint species distribution model (JSDM) will allow for valid inference of environmental effects. We built a joint species distribution model of marine fish and invertebrates of the Northeast US Continental Shelf, providing evidence on species relationships with the environment as well as the likelihood of species to covary. Predictive performance is similar to SSDMs but the Bayesian joint modeling approach provides two main advantages over single species modeling: (1) the JSDM directly estimates the significance of environmental effects; and (2) predicted species richness accounts for species dependencies. An additional value of JSDMs is that the conditional prediction of species distributions can use not only the environmental associations of species, but also the presence and abundance of other species when forecasting future climatic associations.

Quantum inspired community detection for analysis of biodiversity change driven by land-use conversion and climate change

Community detection remains little explored in the analysis of biodiversity change. The challenges linked with global biodiversity change have also multiplied manifold in the past few decades. Moreover, most studies concerning biodiversity change lack the quantitative treatment central to species distribution modeling. Empirical analysis of species distribution and abundance is thus integral to the study of biodiversity loss and biodiversity alterations. Community detection is therefore expected to efficiently model the topological aspect of biodiversity change driven by land-use conversion and climate change; given that it has already proven superior for diverse problems in the domain of social network analysis and subgroup discovery in complex systems. Thus, quantum inspired community detection is proposed as a novel technique to predict biodiversity change considering tiger population in eighteen states of India; leading to benchmarking of two novel datasets. Elements of land-use conversion and climate change are explored to design these datasets viz.-Landscape based distribution and Number of tiger reserves based distribution respectively; for predicting regions expected to maximize Tiger population growth. Furthermore, validation of the proposed framework on the said datasets is performed using standard community detection metrics like-Modularity, Normalized Mutual Information (NMI), Adjusted Rand Index (ARI), Degree distribution, Degree centrality and Edge-betweenness centrality. Quantum inspired community detection has also been successful in demonstrating an association between biodiversity change, land-use conversion and climate change; validated statistically by Pearson's correlation coefficient and p value test. Finally, modularity distribution based on parameter tuning establishes the superiority of the second dataset based on the number of Tiger reserves-in predicting regions maximizing Tiger population growth fostering species distribution and abundance; apart from scripting a stronger correlation of biodiversity change with land-use conversion.

The persistence of ancient settlements and urban sustainability

We propose a dedicated research effort on the determinants of settlement persistence in the ancient world, with the potential to significantly advance the scientific understanding of urban sustainability today. Settlements (cities, towns, villages) are locations with two key attributes: They frame human interactions and activities in space, and they are where people dwell or live. Sustainability, in this case, focuses on the capacity of structures and functions of a settlement system (geography, demography, institutions) to provide for continuity of safe habitation. The 7,000-y-old experience of urbanism, as revealed by archaeology and history, includes many instances of settlements and settlement systems enduring, adapting to, or generating environmental, institutional, and technological changes. The field of urban sustainability lacks a firm scientific foundation for understanding the long durée, relying instead on narratives of collapse informed by limited case studies. We argue for the development of a new interdisciplinary research effort to establish scientific understanding of settlement and settlement system persistence. Such an effort would build upon the many fields that study human settlements to develop new theories and databases from the extensive documentation of ancient and premodern urban systems. A scientific foundation will generate novel insights to advance the field of urban sustainability.

Clustering Species With Residual Covariance Matrix in Joint Species Distribution Models

Modeling species distributions over space and time is one of the major research topics in both ecology and conservation biology. Joint Species Distribution models (JSDMs) have recently been introduced as a tool to better model community data, by inferring a residual covariance matrix between species, after accounting for species' response to the environment. However, these models are computationally demanding, even when latent factors, a common tool for dimension reduction, are used. To address this issue, Taylor-Rodriguez et al. (2017) proposed to use a Dirichlet process, a Bayesian nonparametric prior, to further reduce model dimension by clustering species in the residual covariance matrix. Here, we built on this approach to include a prior knowledge on the potential number of clusters, and instead used a Pitman–Yor process to address some critical limitations of the Dirichlet process. We therefore propose a framework that includes prior knowledge in the residual covariance matrix, providing a tool to analyze clusters of species that share the same residual associations with respect to other species. We applied our methodology to a case study of plant communities in a protected area of the French Alps (the Bauges Regional Park), and demonstrated that our extensions improve dimension reduction and reveal additional information from the residual covariance matrix, notably showing how the estimated clusters are compatible with plant traits, endorsing their importance in shaping communities.

Continent-wide tree fecundity driven by indirect climate effects

Indirect climate effects on tree fecundity that come through variation in size and growth (climate-condition interactions) are not currently part of models used to predict future forests. Trends in species abundances predicted from meta-analyses and species distribution models will be misleading if they depend on the conditions of individuals. Here we find from a synthesis of tree species in North America that climate-condition interactions dominate responses through two pathways, i) effects of growth that depend on climate, and ii) effects of climate that depend on tree size. Because tree fecundity first increases and then declines with size, climate change that stimulates growth promotes a shift of small trees to more fecund sizes, but the opposite can be true for large sizes. Change the depresses growth also affects fecundity. We find a biogeographic divide, with these interactions reducing fecundity in the West and increasing it in the East. Continental-scale responses of these forests are thus driven largely by indirect effects, recommending management for climate change that considers multiple demographic rates.

What processes must we understand to forecast regional-scale population dynamics?

An urgent challenge facing biologists is predicting the regional-scale population dynamics of species facing environmental change. Biologists suggest that we must move beyond predictions based on phenomenological models and instead base predictions on underlying processes. For example, population biologists, evolutionary biologists, community ecologists and ecophysiologists all argue that the respective processes they study are essential. Must our models include processes from all of these fields? We argue that answering this critical question is ultimately an empirical exercise requiring a substantial amount of data that have not been integrated for any system to date. To motivate and facilitate the necessary data collection and integration, we first review the potential importance of each mechanism for skilful prediction. We then develop a conceptual framework based on reaction norms, and propose a hierarchical Bayesian statistical framework to integrate processes affecting reaction norms at different scales. The ambitious research programme we advocate is rapidly becoming feasible due to novel collaborations, datasets and analytical tools.

The importance of genomic variation for biodiversity, ecosystems and people

The 2019 United Nations Global assessment report on biodiversity and ecosystem services estimated that approximately 1 million species are at risk of extinction. This primarily human-driven loss of biodiversity has unprecedented negative consequences for ecosystems and people. Classic and emerging approaches in genetics and genomics have the potential to dramatically improve these outcomes. In particular, the study of interactions among genetic loci within and between species will play a critical role in understanding the adaptive potential of species and communities, and hence their direct and indirect effects on biodiversity, ecosystems and people. We explore these population and community genomic contexts in the hope of finding solutions for maintaining and improving ecosystem services and nature's contributions to people.