From theory to practice in pattern-oriented modelling: identifying and using empirical patterns in predictive models

Predicting emergent animal biodiversity patterns across multiple scales

Restoring biodiversity-based resilience and ecosystem multi-functionality needs to be informed by more accurate predictions of animal biodiversity responses to environmental change. Ecological models make a substantial contribution to this understanding, especially when they encode the biological mechanisms and processes that give rise to emergent patterns (population, community, ecosystem properties and dynamics). Here, a distinction between 'mechanistic' and 'process-based' ecological models is established to review existing approaches. Mechanistic and process-based ecological models have made key advances to understanding the structure, function and dynamics of animal biodiversity, but are typically designed to account for specific levels of biological organisation and spatiotemporal scales. Cross-scale ecological models, which predict emergent co-occurring biodiversity patterns at interacting scales of space, time and biological organisation, is a critical next step in predictive ecology. A way forward is to first capitalise on existing models to systematically evaluate the ability of scale-explicit mechanisms and processes to predict emergent patterns at alternative scales. Such model intercomparisons will reveal mechanism to process transitions across fine to broad scales, overcome approach-specific barriers to model realism or tractability and identify gaps which necessitate the development of new fundamental principles. Key challenges surrounding model complexity and uncertainty would need to be addressed, and while opportunities from big data can streamline the integration of multiple scale-explicit biodiversity patterns, ambitious cross-scale field studies are also needed. Crucially, overcoming cross-scale ecological modelling challenges would unite disparate fields of ecology with the common goal of improving the evidence-base to safeguard biodiversity and ecosystems under novel environmental change.

Climate change may shift metapopulations towards unstable source-sink dynamics in a fire-killed, serotinous shrub

Climate change, with warming and drying weather conditions, is reducing the growth, seed production, and survival of fire-adapted plants in fire-prone regions such as Mediterranean-type ecosystems. These effects of climate change on local plant demographics have recently been shown to reduce the persistence time of local populations of the fire-killed shrub Banksia hookeriana dramatically. In principle, extinctions of local populations may be partly compensated by recolonization events through long-distance dispersal mechanisms of seeds, such as post-fire wind and bird-mediated dispersal, facilitating persistence in spatially structured metapopulations. However, to what degree and under which assumptions metapopulation dynamics might compensate for the drastically increased local extinction risk remains to be explored. Given the long timespans involved and the complexity of interwoven local and regional processes, mechanistic, process-based models are one of the most suitable approaches to systematically explore the potential role of metapopulation dynamics and its underlying ecological assumptions for fire-prone ecosystems. Here we extend a recent mechanistic, process-based, spatially implicit population model for the well-studied fire-killed and serotinous shrub species B. hookeriana to a spatially explicit metapopulation model. We systematically tested the effects of different ecological processes and assumptions on metapopulation dynamics under past (1988-2002) and current (2003-2017) climatic conditions, including (i) effects of different spatio-temporal fires, (ii) effects of (likely) reduced intraspecific plant competition under current conditions and (iii) effects of variation in plant performance among and within patches. In general, metapopulation dynamics had the potential to increase the overall regional persistence of B. hookeriana. However, increased population persistence only occurred under specific optimistic assumptions. In both climate scenarios, the highest persistence occurred with larger fires and intermediate to long inter-fire intervals. The assumption of lower intraspecific plant competition caused by lower densities under current conditions alone was not sufficient to increase persistence significantly. To achieve long-term persistence (defined as >400 years) it was necessary to additionally consider empirically observed variation in plant performance among and within patches, that is, improved habitat quality in some large habitat patches (≥7) that could function as source patches and a higher survival rate and seed production for a subset of plants, specifically the top 25% of flower producers based on current climate conditions monitoring data. Our model results demonstrate that the impacts of ongoing climate change on plant demographics are so severe that even under optimistic assumptions, the existing metapopulation dynamics shift to an unstable source-sink dynamic state. Based on our findings, we recommend increased research efforts to understand the consequences of intraspecific trait variation on plant demographics, emphasizing the variation of individual traits both among and within populations. From a conservation perspective, we encourage fire and land managers to revise their prescribed fire plans, which are typically short interval, small fires, as they conflict with the ecologically appropriate spatio-temporal fire regime for B. hookeriana, and likely as well for many other fire-killed species.

It's about time: Feeding competition costs of sociality are affected more by temporal characteristics than spatial distribution

For most herbivorous animals, group-living appears to incur a high cost by intensifying feeding competition. These costs raise the question of how gregariousness (i.e., the tendency to aggregate) could have evolved to such an extent in taxa such as anthropoid primates and ungulates. When attempting to test the potential benefits and costs, previous foraging models demonstrated that group-living might be beneficial by lowering variance in intake, but that it reduces overall foraging success. However, these models did not fully account for the fact that gregariousness has multiple experiences and can vary in relation to ecological variables and foraging competition. Here, we present an agent-based model for testing how ecological variables impact the costs and benefits of gregariousness. In our simulations, primate-like agents forage on a variable resource landscape while maintaining spatial cohesion with conspecifics to varying degrees. The agents' energy intake rate, daily distance traveled, and variance in energy intake were recorded. Using Morris Elementary Effects sensitivity analysis, we tested the sensitivity of 10 model parameters, of which 2 controlled gregarious behavior and 8 controlled food resources, including multiple aspects of temporal and spatial heterogeneity. We found that, while gregariousness generally increased feeding competition, the costs of gregariousness were much lower when resources were less variable over time (i.e., when calorie extraction was slow and resource renewal was frequent). We also found that maintaining proximity to other agents resulted in lower variance in energy intake when resources were more variable over time. Thus, it appears that the costs and benefits of gregariousness are strongly influenced by the temporal characteristics of food resources, giving insight into the pressures that shaped the evolution of sociality and group living, including in our own lineage.

Extrapolating Metal (Cu, Ni, Zn) Toxicity from Individuals to Populations Across Daphnia Species Using Mechanistic Models: The Roles of Uncertainty Propagation and Combined Physiological Modes of Action

Mechanistic effect modeling is a promising tool to improve the ecological realism of environmental risk assessment. An open question for the mechanistic modeling of metal toxicity is whether the same physiological mode of action (PMoA) could be assumed for closely related species. The implications of various modeling choices, such as the use of parameter point estimates and assumption of simplistic toxicodynamic models, are largely unexplored. We conducted life-table experiments with Daphnia longispina, Daphnia magna, and Daphnia pulex exposed to the single metals Cu, Ni, and Zn, and calibrated toxicokinetic-toxicodynamic (TKTD) models based on dynamic energy budget theory. We developed TKTD models with single and combined PMoAs to compare their goodness-of-fit and predicted population-level sensitivity. We identified the PMoA reproduction efficiency as most probable in all species for Ni and Zn, but not for Cu, and found that combined-PMoA models predicted higher population-level sensitivity than single-PMoA models, which was related to the predicted individual-level sensitivity, rather than to mechanistic differences between models. Using point estimates of parameters, instead of sampling from the probability distributions of parameters, could also lead to differences in the predicted population-level sensitivity. According to model predictions, apical chronic endpoints (cumulative reproduction, survival) are conservative for single-metal population effects across metals and species. We conclude that the assumption of an identical PMoA for different species of Daphnia could be justified for Ni and Zn, but not for Cu. Single-PMoA models are more appropriate than combined-PMoA models from a model selection perspective, but propagation of the associated uncertainty should be considered. More accurate predictions of effects at low concentrations may nevertheless motivate the use of combined-PMoA models. Environ Toxicol Chem 2024;43:338-358. © 2023 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Mechanistic population models for ecological risk assessment and decision support: The importance of good conceptual model diagrams

The use of mechanistic population models as research and decision-support tools in ecology and ecological risk assessment (ERA) is increasing. This growth has been facilitated by advances in technology, allowing the simulation of more complex systems, as well as by standardized approaches for model development, documentation, and evaluation. Mechanistic population models are particularly useful for simulating complex systems, but the required model complexity can make them challenging to communicate. Conceptual diagrams that summarize key model elements, as well as elements that were considered but not included, can facilitate communication and understanding of models and increase their acceptance as decision-support tools. Currently, however, there are no consistent standards for creating or presenting conceptual model diagrams (CMDs), and both terminology and content vary widely. Here, we argue that greater consistency in CMD development and presentation is an important component of good modeling practice, and we provide recommendations, examples, and a free web app (pop-cmd.com) for achieving this for population models used for decision support in ERAs. Integr Environ Assess Manag 2024;00:1-9. © 2023 SETAC.

Ecology needs to overcome siloed modelling

Bahlburg et al. re-implemented eight growth models of Antarctic krill and showed that their predictions are all over the place. The authors discuss the reasons for this and how more coherence in modelling could be achieved through systematic model comparison and integration. For this, we need a common language.

Individual based models for the simulation of fish movement near barriers: Current work and future directions

River fragmentation is an increasing issue for water managers and conservationists. Barriers such as dams interfere with freshwater fish migration, leading to drastic population declines. While there are a range of widely implemented mitigation approaches, e.g. fish passes, such measures are often inefficient due to suboptimal operation and design. There is increasing need to be able to assess mitigation options prior to implementation. Individual based models (IBMs) are a promising option. IBMs can simulate the fine-scale movement of individual fish within a population as they attempt to find a fish pass, incorporating movement processes themselves. Moreover, IBMs have high transferability to other sites or conditions (e.g. changing mitigation, change in flow conditions), making them potentially valuable for freshwater fish conservation yet their application to the fine-scale movement of fish past barriers is still novel. Here, we present an overview of existing IBMs for fine-scale freshwater fish movement, with emphasis on study species and the parameters driving movement in the models. In this review, we focus on IBMs suitable for the simulation of fish tracks as they approach or pass a single barrier. The selected IBMs for modelling fine-scale freshwater fish movement largely focus on salmonids and cyprinid species. IBMs have many applications in the context of fish passage, such as testing different mitigation options or understanding processes behind movement. Existing IBMs include movement processes such as attraction and rejection behaviours, as reported in literature. Yet some factors affecting fish movement e.g. biotic interactions are not covered by existing IBMs. As the technology available for fine scale data collection continues to advance, such as increasing data linking fish behaviour to hydraulics, IBMs could become a more common tool in the design and implementation of fish bypass structures.

Evolutionary adaptation of trees and modelled future larch forest extent in Siberia

With changing climate, the boreal forest could potentially migrate north and become threatened by droughts in the south. However, whether larches, the dominant tree species in eastern Siberia, can adapt to novel situations is largely unknown but is crucial for predicting future population dynamics. Exploring variable traits and trait adaptation through inheritance in an individual-based model can improve our understanding and help future projections. We updated the individual-based spatially explicit vegetation model LAVESI (Larix Vegetation Simulator), used for forest predictions in Eastern Siberia, with trait value variation and incorporated inheritance of parental values to their offspring. Forcing the model with both past and future climate projections, we simulated two areas - the expanding northern treeline and a southerly area experiencing drought. While the specific trait of 'seed weight' regulates migration, the abstract 'drought resistance' protects stands. We show that trait variation with inheritance leads to an increase in migration rate (∼ 3% area increase until 2100). The drought resistance simulations show that, under increasing stress, including adaptive traits leads to larger surviving populations (17% of threatened under RCP 4.5 (Representative Concentration Pathway)). We show that with the increase expected under the RCP 8.5 scenario vast areas (80% of the extrapolated area) of larch forest are threatened and could disappear due to drought as adaptation plays only a minor role under strong warming. We conclude that variable traits facilitate the availability of variants under environmental changes. Inheritance allows populations to adapt to environments and promote successful traits, which leads to populations that can spread faster and be more resilient, provided the changes are not too drastic in both time and magnitude. We show that trait variation and inheritance contribute to more accurate models that can improve our understanding of responses of boreal forests to global change.

Calibrating spatiotemporal models of microbial communities to microscopy data: A review

Spatiotemporal models that account for heterogeneity within microbial communities rely on single-cell data for calibration and validation. Such data, commonly collected via microscopy and flow cytometry, have been made more accessible by recent advances in microfluidics platforms and data processing pipelines. However, validating models against such data poses significant challenges. Validation practices vary widely between modelling studies; systematic and rigorous methods have not been widely adopted. Similar challenges are faced by the (macrobial) ecology community, in which systematic calibration approaches are often employed to improve quantitative predictions from computational models. Here, we review single-cell observation techniques that are being applied to study microbial communities and the calibration strategies that are being employed for accompanying spatiotemporal models. To facilitate future calibration efforts, we have compiled a list of summary statistics relevant for quantifying spatiotemporal patterns in microbial communities. Finally, we highlight some recently developed techniques that hold promise for improved model calibration, including algorithmic guidance of summary statistic selection and machine learning approaches for efficient model simulation.

Temperature-dependent evolutionary speed shapes the evolution of biodiversity patterns across tetrapod radiations

Biodiversity varies predictably with environmental energy around the globe, but the underlaying mechanisms remain incompletely understood. The evolutionary speed hypothesis predicts that environmental kinetic energy shapes variation in speciation rates through temperature- or life history-dependent rates of evolution. To test whether variation in evolutionary speed can explain the relationship between energy and biodiversity in birds, mammals, amphibians, and reptiles, we simulated diversification over 65 million years of geological and climatic change with a spatially explicit eco-evolutionary simulation model. We modelled four distinct evolutionary scenarios in which speciation-completion rates were dependent on temperature (M1), life history (M2), temperature and life history (M3), or were independent of temperature and life-history (M0). To assess the agreement between simulated and empirical data, we performed model selection by fitting supervised machine learning models to multidimensional biodiversity patterns. We show that a model with temperature-dependent rates of speciation (M1) consistently had the strongest support. In contrast to statistical inferences, which showed no general relationships between temperature and speciation rates in tetrapods, we demonstrate how process-based modelling can disentangle the causes behind empirical biodiversity patterns. Our study highlights how environmental energy has played a fundamental role in the evolution of biodiversity over deep time.

A review of bioenergetic modelling for marine mammal populations

Bioenergetic models describe the processes through which animals acquire energy from resources in the environment and allocate it to different life history functions. They capture some of the fundamental mechanisms regulating individuals, populations and ecosystems and have thus been used in a wide variety of theoretical and applied contexts. Here, I review the development of bioenergetic models for marine mammals and their application to management and conservation. For these long-lived, wide-ranging species, bioenergetic approaches were initially used to assess the energy requirements and prey consumption of individuals and populations. Increasingly, models are developed to describe the dynamics of energy intake and allocation and predict how resulting body reserves, vital rates and population dynamics might change as external conditions vary. The building blocks required to develop such models include estimates of intake rate, maintenance costs, growth patterns, energy storage and the dynamics of gestation and lactation, as well as rules for prioritizing allocation. I describe how these components have been parameterized for marine mammals and highlight critical research gaps. Large variation exists among available analytical approaches, reflecting the large range of life histories, management needs and data availability across studies. Flexibility in modelling strategy has supported tailored applications to specific case studies but has resulted in limited generality. Despite the many empirical and theoretical uncertainties that remain, bioenergetic models can be used to predict individual and population responses to environmental change and other anthropogenic impacts, thus providing powerful tools to inform effective management and conservation.

Energy-mediated responses to changing prey size and distribution in marine top predator movements and population dynamics

Climate change is modifying the structure of marine ecosystems, including that of fish communities. Alterations in abiotic and biotic conditions can decrease fish size and change community spatial arrangement, ultimately impacting predator species which rely on these communities. To conserve predators and understand the drivers of observed changes in their population dynamics, we must advance our understanding of how shifting environmental conditions can impact populations by limiting food available to individuals. To investigate the impacts of changing fish size and spatial aggregation on a top predator population, we applied an existing agent-based model parameterized for harbour porpoises Phocoena phocoena which represents animal energetics and movements in high detail. We used this framework to quantify the impacts of shifting prey size and spatial aggregation on porpoise movement, space use, energetics and population dynamics. Simulated individuals were more likely to switch from area-restricted search to transit behaviour with increasing prey size, particularly when starving, due to elevated resource competition. In simulations with highly aggregated prey, higher prey encounter rates counteracted resource competition, resulting in no impacts of prey spatial aggregation on movement behaviour. Reduced energy intake with decreasing prey size and aggregation level caused population decline, with a 15% decrease in fish length resulting in total population collapse Increasing prey consumption rates by 42.8 ± 4.5% could offset population declines; however, this increase was 21.3 ± 12.7% higher than needed to account for changes in total energy availability alone. This suggests that animals in realistic seascapes require additional energy to locate smaller prey which should be considered when assessing the impacts of decreased energy availability. Changes in prey size and aggregation influenced movements and population dynamics of simulated harbour porpoises, revealing that climate-induced changes in prey structure, not only prey abundance, may threaten predator populations. We demonstrate how a population model with realistic animal movements and process-based energetics can be used to investigate population consequences of shifting food availability, such as those mediated by climate change, and provide a mechanistic explanation for how changes in prey structure can impact energetics, behaviour and ultimately viability of predator populations.