Cross-scale integration of knowledge for predicting species ranges: a metamodeling framework

Combining observational and experimental data to estimate environmental and species drivers of fungal metacommunity dynamics

Understanding the distribution and dynamics of species is central to ecology and important for managing biodiversity. The distributions of species in metacommunities are determined by many factors, including environmental conditions and interactions between species. Yet, it is difficult to quantify the effect of species interactions on metacommunity dynamics from observational data. We present an approach to estimate the importance of species interactions that combines data from two observational presence-absence inventories (providing colonization-extinction data) with data from species interaction experiments (providing informative prior distributions in the Bayesian framework). We further illustrate the approach on wood-decay fungi that interact within a downed log through competition for resources and space, and facilitate the succession of other species by decomposing the wood. Specifically, we estimated the relative importance of species interactions by examining how the presence of a species influenced the colonization and extinction probability of other species. Temporal data on fruit body occurrence of 12 species inventoried twice were jointly analyzed with experimental data from two laboratory experiments that aimed to estimate competitive interactions. Both environmental variables and species interactions affected colonization and extinction dynamics. Late-successional fungi had more colonization interactions with predecessor species than early-successional species. We identified several species interactions, and the presence of certain species changed the probability that later-successional species colonized by -81% to 512%. The presence of certain species increased the probability that other species went extinct from a log by 14%-61%. Including the informative priors from experimental data added two colonization interactions and one extinction interaction for which the observational field data was inconclusive. However, most species had no detectable interactions, either because they did not interact or because of low species occupancy, meaning data limitation. We show how temporal presence-absence data can be combined with experimental data to identify which species influence the colonization-extinction dynamics of others. Accounting for species interactions in metacommunity models, in addition to environmental drivers, is important because interactions can have cascading effects on other species.

Realised Thermal Niches in Marine Ectotherms Are Shaped by Ontogeny and Trophic Interactions

Understanding the response of marine organisms to temperature is crucial for predicting climate change impacts. Fundamental physiological thermal performance curves (TPCs), determined under controlled conditions, are commonly used to project future species spatial distributions or physiological performances. Yet, real-world performances may deviate due to extrinsic factors covarying with temperature (food, oxygen, etc.). Using a bioenergetic marine ecosystem model, we evaluate the differences between fundamental and realised TPCs for fish species with contrasted ecology and thermal preferences. Food limitation is the primary cause of differences, decreasing throughout ontogeny and across trophic levels due to spatio-temporal variability of low-trophic level prey availability with temperature. Deoxygenation has moderate impact, despite increasing during ontogeny. This highlights the lower sensitivity of early life stages to hypoxia, which is mechanistically explained by lower mass-specific ingestion at older stages. Understanding the emergence of realised thermal niches offers crucial insights to better determine population's persistence under climate warming.

Incorporating physiological knowledge into correlative species distribution models minimizes bias introduced by the choice of calibration area

Correlative species distribution models (SDMs) are important tools to estimate species' geographic distribution across space and time, but their reliability heavily relies on the availability and quality of occurrence data. Estimations can be biased when occurrences do not fully represent the environmental requirement of a species. We tested to what extent species' physiological knowledge might influence SDM estimations. Focusing on the Japanese sea cucumber Apostichopus japonicus within the coastal ocean of East Asia, we compiled a comprehensive dataset of occurrence records. We then explored the importance of incorporating physiological knowledge into SDMs by calibrating two types of correlative SDMs: a naïve model that solely depends on environmental correlates, and a physiologically informed model that further incorporates physiological information as priors. We further tested the models' sensitivity to calibration area choices by fitting them with different buffered areas around known presences. Compared with naïve models, the physiologically informed models successfully captured the negative influence of high temperature on A. japonicus and were less sensitive to the choice of calibration area. The naïve models resulted in more optimistic prediction of the changes of potential distributions under climate change (i.e., larger range expansion and less contraction) than the physiologically informed models. Our findings highlight benefits from incorporating physiological information into correlative SDMs, namely mitigating the uncertainties associated with the choice of calibration area. Given these promising features, we encourage future SDM studies to consider species physiological information where available.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-024-00226-0.

Spatial ecology of moose in Sweden: Combined Sr-O-C isotope analyses of bone and antler

The study of spatial (paleo)ecology in mammals is critical to understand how animals adapt to and exploit their environment. In this work we analysed the 87Sr/86Sr, δ18O and δ13C isotope composition of 65 moose bone and antler samples from Sweden from wild-shot individuals dated between 1800 and 1994 to study moose mobility and feeding behaviour for (paleo)ecological applications. Sr data were compared with isoscapes of the Scandinavian region, built ad-hoc during this study, to understand how moose utilise the landscape in Northern Europe. The 87Sr/86Sr isoscape was developed using a machine-learning approach with external geo-environmental predictors and literature data. Similarly, a δ18O isoscape, obtained from average annual precipitation δ18O values, was employed to highlight differences in the isotope composition of the local environment vs. bone/antler. Overall, 82% of the moose samples were compatible with the likely local isotope composition (n = 53), suggesting that they were shot not far from their year-round dwelling area. 'Local' samples were used to calibrate the two isoscapes, to improve the prediction of provenance for the presumably 'non-local' individuals. For the latter (n = 12, of which two are antlers and ten are bones), the probability of geographic origin was estimated using a Bayesian approach by combining the two isoscapes. Interestingly, two of these samples (one antler and one bone) seem to come from areas more than 250 km away from the place where the animals were hunted, indicating a possible remarkable intra-annual mobility. Finally, the δ13C data were compared with the forest cover of Sweden and ultimately used to understand the dietary preference of moose. We interpreted a difference in δ13C values of antlers (13C-enriched) and bones (13C-depleted) as a joint effect of seasonal variations in moose diet and, possibly, physiological stresses during winter-time, i.e., increased consumption of endogenous 13C-depleted lipids.

Potential sources of time lags in calibrating species distribution models

The Anthropocene is characterized by a rapid pace of environmental change and is causing a multitude of biotic responses, including those that affect the spatial distribution of species. Lagged responses are frequent and species distributions and assemblages are consequently pushed into a disequilibrium state. How the characteristics of environmental change-for example, gradual 'press' disturbances such as rising temperatures due to climate change versus infrequent 'pulse' disturbances such as extreme events-affect the magnitude of responses and the relaxation times of biota has been insufficiently explored. It is also not well understood how widely used approaches to assess or project the responses of species to changing environmental conditions can deal with time lags. It, therefore, remains unclear to what extent time lags in species distributions are accounted for in biodiversity assessments, scenarios and models; this has ramifications for policymaking and conservation science alike. This perspective piece reflects on lagged species responses to environmental change and discusses the potential consequences for species distribution models (SDMs), the tools of choice in biodiversity modelling. We suggest ways to better account for time lags in calibrating these models and to reduce their leverage effects in projections for improved biodiversity science and policy.

Uniting Experiments and Big Data to advance ecology and conservation

Many ecologists increasingly advocate for research frameworks centered on the use of 'big data' to address anthropogenic impacts on ecosystems. Yet, experiments are often considered essential for identifying mechanisms and informing conservation interventions. We highlight the complementarity of these research frameworks and expose largely untapped opportunities for combining them to speed advancements in ecology and conservation. With nascent but increasing application of model integration, we argue that there is an urgent need to unite experimental and big data frameworks throughout the scientific process. Such an integrated framework offers potential for capitalizing on the benefits of both frameworks to gain rapid and reliable answers to ecological challenges.

Integrating plant physiology into simulation of fire behavior and effects

Wildfires are a global crisis, but current fire models fail to capture vegetation response to changing climate. With drought and elevated temperature increasing the importance of vegetation dynamics to fire behavior, and the advent of next generation models capable of capturing increasingly complex physical processes, we provide a renewed focus on representation of woody vegetation in fire models. Currently, the most advanced representations of fire behavior and biophysical fire effects are found in distinct classes of fine-scale models and do not capture variation in live fuel (i.e. living plant) properties. We demonstrate that plant water and carbon dynamics, which influence combustion and heat transfer into the plant and often dictate plant survival, provide the mechanistic linkage between fire behavior and effects. Our conceptual framework linking remotely sensed estimates of plant water and carbon to fine-scale models of fire behavior and effects could be a critical first step toward improving the fidelity of the coarse scale models that are now relied upon for global fire forecasting. This process-based approach will be essential to capturing the influence of physiological responses to drought and warming on live fuel conditions, strengthening the science needed to guide fire managers in an uncertain future.

Forecasting the future of life in Antarctica

Antarctic ecosystems are under increasing anthropogenic pressure, but efforts to predict the responses of Antarctic biodiversity to environmental change are hindered by considerable data challenges. Here, we illustrate how novel data capture technologies provide exciting opportunities to sample Antarctic biodiversity at wider spatiotemporal scales. Data integration frameworks, such as point process and hierarchical models, can mitigate weaknesses in individual data sets, improving confidence in their predictions. Increasing process knowledge in models is imperative to achieving improved forecasts of Antarctic biodiversity, which can be attained for data-limited species using hybrid modelling frameworks. Leveraging these state-of-the-art tools will help to overcome many of the data scarcity challenges presented by the remoteness of Antarctica, enabling more robust forecasts both near- and long-term.

Integrating functional traits into correlative species distribution models to investigate the vulnerability of marine human activities to climate change

Climate change and particularly warming are significantly impacting marine ecosystems and the services they provided. Temperature, as the main factor driving all biological processes, may influence ectotherms metabolism, thermal tolerance limits and distribution species patterns. The joining action of climate change and local stressors (including the increasing human marine use) may facilitate the spread of non-indigenous and native outbreak forming species, leading to associated economic consequences for marine coastal economies. Marine aquaculture is one among the most economic anthropogenic activities threatened by multiple stressors and in turn, by increasing hard artificial substrates at sea would facilitate the expansion of these problematic organisms and face negative consequences regarding facilities management and farmed organisms' welfare. Species Distribution Models (SDMs) are considered powerful tools for forecasting the future occurrences and distributions of problematic species used to preventively aware stakeholders. In the current study, we propose the use of combined correlative SDMs and mechanistic models, based on individual thermal performance curve models calculated through non-linear least squares regression and Bayesian statistics (functional-SDM), as an ecological relevant tool to increase our ability to investigate the potential indirect effect of climate change on the distributions of harmful species for human activities at sea, taking aquaculture as a food productive example and the benthic cnidarian Pennaria disticha (one of the most pernicious fouling species in aquaculture) as model species. Our combined approach was able to improve the prediction ability of both mechanistic and correlative models to get more ecologically informed "whole" niche of the studied species. Incorporating the mechanistic links between the organisms' functional traits and their environments into SDMs through the use of a Bayesian functional-SDM approach would be a useful and reliable tool in early warning ecological systems, risk assessment and management actions focused on important economic activities and natural ecosystems conservation.

Understanding and Modeling Forest Disturbance Interactions at the Landscape Level

Disturbances, both natural and anthropogenic, affect the configuration, composition, and function of forested ecosystems. Complex system behaviors emerge from the interactions between disturbance regimes, the vegetation response to those disturbances, and their interplay with multiple drivers (climate, topography, land use, etc.) across spatial and temporal scales. Here, we summarize conceptual advances and empirical approaches to disturbance interaction investigation, and used those insights to evaluate and categorize 146 landscape modeling studies emerging from a systematic review of the literature published since 2010. Recent conceptual advances include formal disaggregation of disturbances into their constituent components, embedding disturbance processes into system dynamics, and clarifying terminology for interaction factors, types, and ecosystem responses. Empirical studies investigating disturbance interactions now span a wide range of approaches, including (most recently) advanced statistical methods applied to an expanding set of spatial and temporal datasets. Concurrent development in spatially-explicit landscape models, informed by these empirical insights, integrate the interactions among natural and anthropogenic disturbances by coupling these processes to account for disturbance stochasticity, disturbance within and across scales, and non-linear landscape responses to climate change. Still, trade-offs between model elegance and complexity remain. We developed an index for the degree of process integration (i.e., balance of static vs. dynamic components) within a given disturbance agent and applied it to the studies from our systematic review. Contemporary model applications in this line of research have applied a wide range process integration, depending on the specific question, but also limited in part by data and knowledge. Non-linear “threshold” behavior and cross-scaled interactions remain a frontier in temperate, boreal, and alpine regions of North America and Europe, while even simplistic studies are lacking from other regions of the globe (e.g., subtropical and tropical biomes). Understanding and planning for uncertainty in system behavior—including disturbance interactions—is paramount at a time of accelerated anthropogenic change. While progress in landscape modeling studies in this area is evident, work remains to increase model transparency and confidence, especially for understudied regions and processes. Moving forward, a multi-dimensional approach is recommended to address the uncertainties of complex human-ecological dynamics.

Ecological Dynamics: Integrating Empirical, Statistical, and Analytical Methods

Understanding ecological processes and predicting long-term dynamics are ongoing challenges in ecology. To address these challenges, we suggest an approach combining mathematical analyses and Bayesian hierarchical statistical modeling with diverse data sources. Novel mathematical analysis of ecological dynamics permits a process-based understanding of conditions under which systems approach equilibrium, experience large oscillations, or persist in transient states. This understanding is improved by combining ecological models with empirical observations from a variety of sources. Bayesian hierarchical models explicitly couple process-based models and data, yielding probabilistic quantification of model parameters, system characteristics, and associated uncertainties. We outline relevant tools from dynamical analysis and hierarchical modeling and argue for their integration, demonstrating the value of this synthetic approach through a simple predator-prey example.

Integrating intraseasonal grassland dynamics in cross-scale distribution modeling to support waterbird recovery plans

Despite much discussion about the utility of remote sensing for effective conservation, the inclusion of these technologies in species recovery plans remains largely anecdotal. We developed a modeling approach for the integration of local, spatially measured ecosystem functional dynamics into a species distribution modeling (SDM) framework in which other ecologically relevant factors are modeled separately at broad scales. To illustrate the approach, we incorporated intraseasonal water-vegetation dynamics into a cross-scale SDM for the Common Snipe (Gallinago gallinago), which is highly dependent on water and vegetation dynamics. The Common Snipe is an Iberian grassland waterbird characteristic of European agricultural meadows and a member of one of the most threatened bird guilds. The intraseasonal dynamics of water content of vegetation were measured using the standard deviation of the normalized difference water index time series computed from bimonthly images of the Sentinel-2 satellite. The recovery plan for the Common Snipe in Galicia (northwestern Iberian Peninsula) provided an opportunity to apply our modeling framework. Model accuracy in predicting the species' distribution at a regional scale (resulting from integration of downscaled climate projections with regional habitat-topographic suitability models) was very high (area under the curve [AUC] of 0.981 and Boyce's index of 0.971). Local water-vegetation dynamic models, based exclusively on Sentinel-2 imagery, were good predictors (AUC of 0.849 and Boyce's index of 0.976). The predictive power improved (AUC of 0.92 and Boyce's index of 0.98) when local model predictions were restricted to areas identified by the continental and regional models as priorities for conservation. Our models also performed well (AUC of 0.90 and Boyce's index of 0.93) when projected to updated water-vegetation conditions. Our modeling framework enabled incorporation of key ecosystem processes closely related to water and carbon cycles while accounting for other factors ecologically relevant to endangered grassland waterbirds across different scales, allowed identification of priority areas for conservation, and provided an opportunity for cost-effective recovery planning by monitoring management effectiveness from space.

A practical guide for combining data to model species distributions

Understanding and accurately modeling species distributions lies at the heart of many problems in ecology, evolution, and conservation. Multiple sources of data are increasingly available for modeling species distributions, such as data from citizen science programs, atlases, museums, and planned surveys. Yet reliably combining data sources can be challenging because data sources can vary considerably in their design, gradients covered, and potential sampling biases. We review, synthesize, and illustrate recent developments in combining multiple sources of data for species distribution modeling. We identify five ways in which multiple sources of data are typically combined for modeling species distributions. These approaches vary in their ability to accommodate sampling design, bias, and uncertainty when quantifying environmental relationships in species distribution models. Many of the challenges for combining data are solved through the prudent use of integrated species distribution models: models that simultaneously combine different data sources on species locations to quantify environmental relationships for explaining species distribution. We illustrate these approaches using planned survey data on 24 species of birds coupled with opportunistically collected eBird data in the southeastern United States. This example illustrates some of the benefits of data integration, such as increased precision in environmental relationships, greater predictive accuracy, and accounting for sample bias. Yet it also illustrates challenges of combining data sources with vastly different sampling methodologies and amounts of data. We provide one solution to this challenge through the use of weighted joint likelihoods. Weighted joint likelihoods provide a means to emphasize data sources based on different criteria (e.g., sample size), and we find that weighting improves predictions for all species considered. We conclude by providing practical guidance on combining multiple sources of data for modeling species distributions.

Biodiversity Models: What If Unsaturation Is the Rule?

Improving biodiversity predictions is essential if we are to meet the challenges posed by global change. As knowledge is key to feed models, we need to evaluate how debated theory can affect models. An important ongoing debate is whether environmental constraints limit the number of species that can coexist in a community (saturation), with recent findings suggesting that species richness in many communities might be unsaturated. Here, we propose that biodiversity models could address this issue by accounting for a duality: considering communities as unsaturated but where species composition is constrained by different scale-dependent biodiversity drivers. We identify a variety of promising advances for incorporating this duality into commonly applied biodiversity modelling approaches and improving their spatial predictions.

The majority of biodiversity modelling approaches do not explicitly address the question of saturation.

Theoretical and methodological implications of saturation or unsaturation in biodiversity modelling.

Addressing saturation or unsaturation is vital to produce more reliable conservation strategies.

Integrative community modelling frameworks may be the way forward.

Towards Process-based Range Modeling of Many Species

Understanding and forecasting species' geographic distributions in the face of global change is a central priority in biodiversity science. The existing view is that one must choose between correlative models for many species versus process-based models for few species. We suggest that opportunities exist to produce process-based range models for many species, by using hierarchical and inverse modeling to borrow strength across species, fill data gaps, fuse diverse data sets, and model across biological and spatial scales. We review the statistical ecology and population and range modeling literature, illustrating these modeling strategies in action. A variety of large, coordinated ecological datasets that can feed into these modeling solutions already exist, and we highlight organisms that seem ripe for the challenge.

Dominant forest tree species are potentially vulnerable to climate change over large portions of their range even at high latitudes

Projecting suitable conditions for a species as a function of future climate provides a reasonable, although admittedly imperfect, spatially explicit estimate of species vulnerability associated with climate change. Projections emphasizing range shifts at continental scale, however, can mask contrasting patterns at local or regional scale where management and policy decisions are made. Moreover, models usually show potential for areas to become climatically unsuitable, remain suitable, or become suitable for a particular species with climate change, but each of these outcomes raises markedly different ecological and management issues. Managing forest decline at sites where climatic stress is projected to increase is likely to be the most immediate challenge resulting from climate change. Here we assess habitat suitability with climate change for five dominant tree species of eastern North American forests, focusing on areas of greatest vulnerability (loss of suitability in the baseline range) in Quebec (Canada) rather than opportunities (increase in suitability). Results show that these species are at risk of maladaptation over a remarkably large proportion of their baseline range. Depending on species, 5–21% of currently climatically suitable habitats are projected to be at risk of becoming unsuitable. This suggests that species that have traditionally defined whole regional vegetation assemblages could become less adapted to these regions, with significant impact on ecosystems and forest economy. In spite of their well-recognised limitations and the uncertainty that remains, regionally-explicit risk assessment approaches remain one of the best options to convey that message and the need for climate policies and forest management adaptation strategies.