Refine
Year of publication
Language
- English (80)
Is part of the Bibliography
- yes (80)
Keywords
- individual-based model (8)
- Individual-based model (7)
- modeling (4)
- BEEHAVE (3)
- Climate change (3)
- Risk assessment (3)
- Stage-based model (3)
- biodiversity (3)
- coexistence (3)
- foraging (3)
The BEEHAVE model simulates the population dynamics and foraging activity of a single honey bee colony (Apis mellifera) in great detail. Although it still makes numerous simplifying assumptions, it appears to capture a wide range of empirical observations.
It could, therefore, in principle, also be used as a tool in beekeeper education, as it allows the implementation and comparison of different management options.
Here, we focus on treatments aimed at controlling the mite Varroa destructor. However, since BEEHAVE was developed in the UK, mite treatment includes the use of a synthetic acaricide, which is not part of Good Beekeeping Practice in Germany.
A practice that consists of drone brood removal from April to June, treatment with formic acid in August/September, and treatment with oxalic acid in November/December. We implemented these measures, focusing on the timing, frequency, and spacing between drone brood removals.
The effect of drone brood removal and acid treatment, individually or in combination, on a mite-infested colony was examined. We quantify the efficacy of Varroa mite control as the reduction of mites in treated bee colonies compared to untreated bee colonies. We found that drone brood removal was very effective, reducing mites by 90% at the end of the first simulation year after the introduction of mites. This value was significantly higher than the 50-67% reduction expected by bee experts and confirmed by empirical studies.
However, literature reports varying percent reductions in mite numbers from 10 to 85% after drone brood removal. The discrepancy between model results, empirical data, and expert estimates indicate that these three sources should be reviewed and refined, as all are based on simplifying assumptions.
These results and the adaptation of BEEHAVE to the Good Beekeeping Practice are a decisive step forward for the future use of BEEHAVE in beekeeper education in Germany and anywhere where organic acids and drone brood removal are utilized.
Riverine ecosystems provide various ecosystem services. One of these services is the biological control of eutrophication by grazing macroinvertebrates.
However, riverine ecosystems are subject to numerous stressors that affect community structure, functions, and stability properties. To manage rivers in response to these stressors, a better understanding of the ecological functions underlying services is needed.
This requires consideration of local and regional processes, which requires a metacommunity approach that links local food webs through drift and dispersal. This takes into account long-distance interactions that can compensate for local effects of stressors.
Our modular model MASTIFF (Multiple Aquatic STressors In Flowing Food webs) is stage-structured, spatially explicit, and includes coupled food webs consisting of benthic resource-consumer interactions between biofilm and three competing macroinvertebrate functional types. River segments are unidirectionally connected through organismal drift and bidirectionally connected through dispersal. Climate and land use stressors along the river can be accounted for. Biocontrol of biofilm eutrophication is used as an exemplary functional indicator.
We present the model and the underlying considerations, and show in an exemplary application that explicit consideration of drift and dispersal is essential for understanding the spatiotemporal biocontrol of eutrophication.
The combination of drift and dispersal reduced eutrophication events. While dispersal events were linked to specific periods in the species life cycles and therefore had limited potential to control, drift was ubiquitous and thus responded more readily to changing habitat conditions.
This indicates that drift is an important factor for coping with stress situations.
Finally, we outline and discuss the potential and possibilities of MASTIFF as a tool for mechanistic, cross-scale analyses of multiple stressors to advance knowledge of riverine ecosystem functioning.
Both climate change and land use regimes affect the viability of populations, but they are often studied separately. Moreover, population viability analyses (PVAs) often ignore the effects of large environmental gradients and use temporal resolutions that are too coarse to take into account that different stages of a population's life cycle may be affected differently by climate change. Here, we present the High-resolution Large Environmental Gradient (HiLEG) model and apply it in a PVA with daily resolution based on daily climate projections for Northwest Germany. We used the large marsh grasshopper (LMG) as the target species and investigated (1) the effects of climate change on the viability and spatial distribution of the species, (2) the influence of the timing of grassland mowing on the species and (3) the interaction between the effects of climate change and grassland mowing. The stageand cohort-based model was run for the spatially differentiated environmental conditions temperature and soil moisture across the whole study region. We implemented three climate change scenarios and analyzed the population dynamics for four consecutive 20-year periods. Climate change alone would lead to an expansion of the regions suitable for the LMG, as warming accelerates development and due to reduced drought stress. However, in combination with land use, the timing of mowing was crucial, as this disturbance causes a high mortality rate in the aboveground life stages. Assuming the same date of mowing throughout the region, the impact on viability varied greatly between regions due to the different climate conditions. The regional negative effects of the mowing date can be divided into five phases: (1) In early spring, the populations were largely unaffected in all the regions; (2) between late spring and early summer, they were severely affected only in warm regions; (3) in summer, all the populations were severely affected so that they could hardly survive; (4) between late summer and early autumn, they were severely affected in cold regions; and (5) in autumn, the populations were equally affected across all regions. The duration and start of each phase differed slightly depending on the climate change scenario and simulation period, but overall, they showed the same pattern. Our model can be used to identify regions of concern and devise management recommendations. The model can be adapted to the life cycle of different target species, climate projections and disturbance regimes. We show with our adaption of the HiLEG model that high-resolution PVAs and applications on large environmental gradients can be reconciled to develop conservation strategies capable of dealing with multiple stressors.
Forage availability has been suggested as one driver of the observed decline in honey bees. However, little is known about the effects of its spatiotemporal variation on colony success. We present a modeling framework for assessing honey bee colony viability in cropping systems. Based on two real farmland structures, we developed a landscape generator to design cropping systems varying in crop species identity, diversity, and relative abundance. The landscape scenarios generated were evaluated using the existing honey bee colony model BEEHAVE, which links foraging to in-hive dynamics. We thereby explored how different cropping systems determine spatiotemporal forage availability and, in turn, honey bee colony viability (e.g., time to extinction, TTE) and resilience (indicated by, e.g., brood mortality). To assess overall colony viability, we developed metrics,P(H)andP(P,)which quantified how much nectar and pollen provided by a cropping system per year was converted into a colony's adult worker population. Both crop species identity and diversity determined the temporal continuity in nectar and pollen supply and thus colony viability. Overall farmland structure and relative crop abundance were less important, but details mattered. For monocultures and for four-crop species systems composed of cereals, oilseed rape, maize, and sunflower,P(H)andP(P)were below the viability threshold. Such cropping systems showed frequent, badly timed, and prolonged forage gaps leading to detrimental cascading effects on life stages and in-hive work force, which critically reduced colony resilience. Four-crop systems composed of rye-grass-dandelion pasture, trefoil-grass pasture, sunflower, and phacelia ensured continuous nectar and pollen supply resulting in TTE > 5 yr, andP(H)(269.5 kg) andP(P)(108 kg) being above viability thresholds for 5 yr. Overall, trefoil-grass pasture, oilseed rape, buckwheat, and phacelia improved the temporal continuity in forage supply and colony's viability. Our results are hypothetical as they are obtained from simplified landscape settings, but they nevertheless match empirical observations, in particular the viability threshold. Our framework can be used to assess the effects of cropping systems on honey bee viability and to develop land-use strategies that help maintain pollination services by avoiding prolonged and badly timed forage gaps.
Forage availability has been suggested as one driver of the observed decline in honey bees. However, little is known about the effects of its spatiotemporal variation on colony success. We present a modeling framework for assessing honey bee colony viability in cropping systems. Based on two real farmland structures, we developed a landscape generator to design cropping systems varying in crop species identity, diversity, and relative abundance. The landscape scenarios generated were evaluated using the existing honey bee colony model BEEHAVE, which links foraging to in-hive dynamics. We thereby explored how different cropping systems determine spatiotemporal forage availability and, in turn, honey bee colony viability (e.g., time to extinction, TTE) and resilience (indicated by, e.g., brood mortality). To assess overall colony viability, we developed metrics,P(H)andP(P,)which quantified how much nectar and pollen provided by a cropping system per year was converted into a colony's adult worker population. Both crop species identity and diversity determined the temporal continuity in nectar and pollen supply and thus colony viability. Overall farmland structure and relative crop abundance were less important, but details mattered. For monocultures and for four-crop species systems composed of cereals, oilseed rape, maize, and sunflower,P(H)andP(P)were below the viability threshold. Such cropping systems showed frequent, badly timed, and prolonged forage gaps leading to detrimental cascading effects on life stages and in-hive work force, which critically reduced colony resilience. Four-crop systems composed of rye-grass-dandelion pasture, trefoil-grass pasture, sunflower, and phacelia ensured continuous nectar and pollen supply resulting in TTE > 5 yr, andP(H)(269.5 kg) andP(P)(108 kg) being above viability thresholds for 5 yr. Overall, trefoil-grass pasture, oilseed rape, buckwheat, and phacelia improved the temporal continuity in forage supply and colony's viability. Our results are hypothetical as they are obtained from simplified landscape settings, but they nevertheless match empirical observations, in particular the viability threshold. Our framework can be used to assess the effects of cropping systems on honey bee viability and to develop land-use strategies that help maintain pollination services by avoiding prolonged and badly timed forage gaps.
In many species, dispersal is decisive for survival in a changing climate. Simulation models for population dynamics under climate change thus need to account for this factor. Moreover, large numbers of species inhabiting agricultural landscapes are subject to disturbances induced by human land use. We included dispersal in the HiLEG model that we previously developed to study the interaction between climate change and agricultural land use in single populations. Here, the model was parameterized for the large marsh grasshopper (LMG) in cultivated grasslands of North Germany to analyze (1) the species development and dispersal success depending on the severity of climate change in subregions, (2) the additional effect of grassland cover on dispersal success, and (3) the role of dispersal in compensating for detrimental grassland mowing. Our model simulated population dynamics in 60-year periods (2020-2079) on a fine temporal (daily) and high spatial (250 x 250 m(2)) scale in 107 subregions, altogether encompassing a range of different grassland cover, climate change projections, and mowing schedules. We show that climate change alone would allow the LMG to thrive and expand, while grassland cover played a minor role. Some mowing schedules that were harmful to the LMG nevertheless allowed the species to moderately expand its range. Especially under minor climate change, in many subregions dispersal allowed for mowing early in the year, which is economically beneficial for farmers. More severe climate change could facilitate LMG expansion to uninhabited regions but would require suitable mowing schedules along the path. These insights can be transferred to other species, given that the LMG is considered a representative of grassland communities. For more specific predictions on the dynamics of other species affected by climate change and land use, the publicly available HiLEG model can be easily adapted to the characteristics of their life cycle.
The pace-of-life syndrome (POLS) hypothesis posits that suites of traits are correlated along a slow-fast continuum owing to life history trade-offs. Despite widespread adoption, environmental conditions driving the emergence of POLS remain unclear. A recently proposed conceptual framework of POLS suggests that a slow-fast continuum should align to fluctuations in density-dependent selection. We tested three key predictions made by this framework with an ecoevolutionary agent-based population model. Selection acted on responsiveness (behavioral trait) to interpatch resource differences and the reproductive investment threshold (life history trait). Across environments with density fluctuations of different magnitudes, we observed the emergence of a common axis of trait covariation between and within populations (i.e., the evolution of a POLS). Slow-type (fast-type) populations with high (low) responsiveness and low (high) reproductive investment threshold were selected at high (low) population densities and less (more) intense and frequent density fluctuations. In support of the predictions, fast-type populations contained a higher degree of variation in traits and were associated with higher intrinsic reproductive rate (r(0)) and higher sensitivity to intraspecific competition (gamma), pointing to a universal trade-off. While our findings support that POLS aligns with density-dependent selection, we discuss possible mechanisms that may lead to alternative evolutionary pathways.
Resilience trinity
(2020)
Ensuring ecosystem resilience is an intuitive approach to safeguard the functioning of ecosystems and hence the future provisioning of ecosystem services (ES). However, resilience is a multi-faceted concept that is difficult to operationalize. Focusing on resilience mechanisms, such as diversity, network architectures or adaptive capacity, has recently been suggested as means to operationalize resilience. Still, the focus on mechanisms is not specific enough. We suggest a conceptual framework, resilience trinity, to facilitate management based on resilience mechanisms in three distinctive decision contexts and time-horizons: 1) reactive, when there is an imminent threat to ES resilience and a high pressure to act, 2) adjustive, when the threat is known in general but there is still time to adapt management and 3) provident, when time horizons are very long and the nature of the threats is uncertain, leading to a low willingness to act. Resilience has different interpretations and implications at these different time horizons, which also prevail in different disciplines. Social ecology, ecology and engineering are often implicitly focussing on provident, adjustive or reactive resilience, respectively, but these different notions of resilience and their corresponding social, ecological and economic tradeoffs need to be reconciled. Otherwise, we keep risking unintended consequences of reactive actions, or shying away from provident action because of uncertainties that cannot be reduced. The suggested trinity of time horizons and their decision contexts could help ensuring that longer-term management actions are not missed while urgent threats to ES are given priority.
Resilience trinity
(2020)
Ensuring ecosystem resilience is an intuitive approach to safeguard the functioning of ecosystems and hence the future provisioning of ecosystem services (ES). However, resilience is a multi-faceted concept that is difficult to operationalize. Focusing on resilience mechanisms, such as diversity, network architectures or adaptive capacity, has recently been suggested as means to operationalize resilience. Still, the focus on mechanisms is not specific enough. We suggest a conceptual framework, resilience trinity, to facilitate management based on resilience mechanisms in three distinctive decision contexts and time-horizons: 1) reactive, when there is an imminent threat to ES resilience and a high pressure to act, 2) adjustive, when the threat is known in general but there is still time to adapt management and 3) provident, when time horizons are very long and the nature of the threats is uncertain, leading to a low willingness to act. Resilience has different interpretations and implications at these different time horizons, which also prevail in different disciplines. Social ecology, ecology and engineering are often implicitly focussing on provident, adjustive or reactive resilience, respectively, but these different notions of resilience and their corresponding social, ecological and economic tradeoffs need to be reconciled. Otherwise, we keep risking unintended consequences of reactive actions, or shying away from provident action because of uncertainties that cannot be reduced. The suggested trinity of time horizons and their decision contexts could help ensuring that longer-term management actions are not missed while urgent threats to ES are given priority.
Building and changing a microbiome at will and maintaining it over hundreds of generations has so far proven challenging. Despite best efforts, complex microbiomes appear to be susceptible to large stochastic fluctuations. Current capabilities to assemble and control stable complex microbiomes are limited. Here, we propose a looped mass transfer design that stabilizes microbiomes over long periods of time. Five local microbiomes were continuously grown in parallel for over 114 generations and connected by a loop to a regional pool. Mass transfer rates were altered and microbiome dynamics were monitored using quantitative high-throughput flow cytometry and taxonomic sequencing of whole communities and sorted subcommunities. Increased mass transfer rates reduced local and temporal variation in microbiome assembly, did not affect functions, and overcame stochasticity, with all microbiomes exhibiting high constancy and increasing resistance. Mass transfer synchronized the structures of the five local microbiomes and nestedness of certain cell types was eminent. Mass transfer increased cell number and thus decreased net growth rates mu'. Subsets of cells that did not show net growth mu'SCx were rescued by the regional pool R and thus remained part of the microbiome. The loop in mass transfer ensured the survival of cells that would otherwise go extinct, even if they did not grow in all local microbiomes or grew more slowly than the actual dilution rate D would allow. The rescue effect, known from metacommunity theory, was the main stabilizing mechanism leading to synchrony and survival of subcommunities, despite differences in cell physiological properties, including growth rates.
Biodiversity loss is a result of interacting ecological and economic factors, and it must be addressed through an analysis of biodiversity conservation policies. Ecological-economic modelling is a helpful approach to this analysis, but it is also challenging since modellers often have a specific disciplinary background and tend to misrepresent either the ecological or economic aspects. Here, we introduce some of the most important concepts from both disciplines, and since the two modelling cultures also differ between the two disciplines, we present an integrated, consistent guide through all the steps of generic ecological-economic modelling, such as formulation of the research question, development of the conceptual model, model parametrisation and analysis, and interpretation of model results. Although we focus on generic models aimed at a general understanding of causes and remedies for biodiversity loss, the concepts and guidance provided here may also help in the modelling of more specific conservation problems. This guide is aimed at the intersection of three disciplines: ecology, economics and mathematical modelling, and addresses readers who have some knowledge in at least one of these disciplines and want to learn about the others to build and analyse generic ecological-economic models. Compared to textbooks, the guide focuses on the practice of modelling rather than lengthy explanations of theoretical concepts. We attempt to demonstrate that generic ecological-economic modelling does not require magical powers and instead is a manageable exercise.
Editorial
(2020)
Movement behavior is an essential element of fundamental ecological processes such as competition and predation. Although intraspecific trait variation (ITV) in movement behaviors is pervasive, its consequences for ecological community dynamics are still not fully understood. Using a newly developed individual-based model, we analyzed how given and constant ITVs in foraging movement affect differences in foraging efficiencies between species competing for common resources under various resource distributions. Further, we analyzed how the effect of ITV on emerging differences in competitive abilities ultimately affects species coexistence. The model is generic but mimics observed patterns of among-individual covariation between personality, movement and space use in ground-dwelling rodents. Interacting species differed in their mean behavioral types along a slow-fast continuum, integrating consistent individual variation in average behavioral expression and responsiveness (i.e. behavioral reaction norms). We found that ITV reduced interspecific differences in competitive abilities by 5-35% and thereby promoted coexistence via an equalizing mechanism. The emergent relationships between behavioral types and foraging efficiency are characteristic for specific environmental contexts of resource distribution and population density. As these relationships are asymmetric, species that were either 'too fast' or 'too slow' benefited differently from ITV. Thus, ITV in movement behavior has consequences for species coexistence but to predict its effect in a given system requires intimate knowledge on how variation in movement traits relates to fitness components along an environmental gradient.
The cultivation of energy crops leads to direct and indirect land use changes that impair the biodiversity of the agricultural landscape. In our study, we analyse the effects of mitigation measures on the European brown hare (Lepus europaeus), which is directly affected by ongoing land use change and has experienced widespread decline throughout Europe since the 1960s. Therefore, we developed a spatially explicit and individual-based ecological model to study the effects of different landscape configurations and compositions on hare population development. As an input, we used two 4 x 4 km large model landscapes, which were generated by a landscape generator based on real field sizes and crop proportions and differed in average field size and crop composition. The crops grown annually are evaluated in terms of forage suitability, breeding suitability and crop richness for the hare. In six mitigation scenarios, we investigated the effects of a 10 % increase in the following measures: (1) mixed silphie, (2) miscanthus, (3) grass-clover ley, (4) alfalfa, (5) set-aside, and (6) general crop richness. All mitigation measures had significant effects on hare population development. Compared to the base scenario, the relative change in hare abundance ranged from a factor of 0.56 in the grass-clover ley scenario to-0.16 in the miscanthus scenario. The mitigation measures of mixed silphie, grass-clover ley and increased crop richness led to distinct increases in hare abundance in both landscapes ( > 0.3). The results show that both landscape configuration and composition have a significant effect on hare population development, which responds particularly strongly to compositional changes. The increase in crop diversity, e.g., through the cultivation of alternative energy crops such as mixed silphie and grass-clover ley, proves to be beneficial for the brown hare.
Pattern-oriented modelling as a novel way to verify and validate functional-structural plant models
(2018)
Background and Aims Functional-structural plant (FSP) models have been widely used to understand the complex interactions between plant architecture and underlying developmental mechanisms. However, to obtain evidence that a model captures these mechanisms correctly, a clear distinction must be made between model outputs used for calibration and thus verification, and outputs used for validation. In pattern-oriented modelling (POM), multiple verification patterns are used as filters for rejecting unrealistic model structures and parameter combinations, while a second, independent set of patterns is used for validation. Key Results After calibration, our model simultaneously reproduced multiple observed architectural patterns. The model then successfully predicted, without further calibration, the validation patterns. The model supports the hypothesis that carbon allocation can be modelled as being dependent on current organ biomass and sink strength of each organ type, and also predicted the observed developmental timing of the leaf sink-source transition stage.
Terrestrial environmental systems are characterised by numerous feedback links between their different compartments. However, scientific research is organized into disciplines that focus on processes within the respective compartments rather than on interdisciplinary links. Major feedback mechanisms between compartments might therefore have been systematically overlooked so far. Without identifying these gaps, initiatives on future comprehensive environmental monitoring schemes and experimental platforms might fail. We performed a comprehensive overview of feedbacks between compartments currently represented in environmental sciences and explores to what degree missing links have already been acknowledged in the literature. We focused on process models as they can be regarded as repositories of scientific knowledge that compile findings of numerous single studies. In total, 118 simulation models from 23 model types were analysed. Missing processes linking different environmental compartments were identified based on a meta-review of 346 published reviews, model inter-comparison studies, and model descriptions. Eight disciplines of environmental sciences were considered and 396 linking processes were identified and ascribed to the physical, chemical or biological domain. There were significant differences between model types and scientific disciplines regarding implemented interdisciplinary links. The most wide-spread interdisciplinary links were between physical processes in meteorology, hydrology and soil science that drive or set the boundary conditions for other processes (e.g., ecological processes). In contrast, most chemical and biological processes were restricted to links within the same compartment. Integration of multiple environmental compartments and interdisciplinary knowledge was scarce in most model types. There was a strong bias of suggested future research foci and model extensions towards reinforcing existing interdisciplinary knowledge rather than to open up new interdisciplinary pathways. No clear pattern across disciplines exists with respect to suggested future research efforts. There is no evidence that environmental research would clearly converge towards more integrated approaches or towards an overarching environmental systems theory. (c) 2017 Elsevier B.V. All rights reserved.
Non-consumptive effects of predators within ecosystems can alter the behavior of individual prey species, and have cascading effects on other trophic levels. In this context, an understanding of non-consumptive predator effects on the whole prey community is crucial for predicting community structure and composition, hence biodiversity patterns. We used an individual-based, spatially-explicit modelling approach to investigate the consequences of landscapes of fear on prey community metrics. The model spans multiple hierarchical levels from individual home range formation based on food availability and perceived predation risk to consequences on prey community structure and composition. This mechanistic approach allowed us to explore how important factors such as refuge availability and foraging strategy under fear affect prey community metrics. Fear of predators affected prey space use, such as home range formation. These adaptations had broader consequences for the community leading to changes in community structure and composition. The strength of community responses to perceived predation risk was driven by refuge availability in the landscape and the foraging strategy of prey animals. Low refuge availability in the landscape strongly decreased diversity and total biomass of prey communities. Additionally, body mass distributions in prey communities facing high predation risk were shifted towards small prey animals. With increasing refuge availability the consequences of non-consumptive predator effects were reduced, diversity and total biomass of the prey community increased. Prey foraging strategies affected community composition. Under medium refuge availability, risk-averse prey communities consisted of many small animals while risk-taking prey communities showed a more even body mass distribution. Our findings reveal that non-consumptive predator effects can have important implications for prey community diversity and should therefore be considered in the context of conservation and nature management.
There is an increasing need for an assessment of the impacts of land use and land use change (LUCC). In this context, simulation models are valuable tools for investigating the impacts of stakeholder actions or policy decisions. Agricultural landscape generators (ALGs), which systematically and automatically generate realistic but simplified representations of land cover in agricultural landscapes, can provide the input for LUCC models. We reviewed existing ALGs in terms of their objectives, design and scope. We found eight ALGs that met our definition. They were based either on generic mathematical algorithms (pattern-based) or on representations of ecological or land use processes (process-based). Most ALGs integrate only a few landscape metrics, which limits the design of the landscape pattern and thus the range of applications. For example, only a few specific farming systems have been implemented. We conclude that existing ALGs contain useful approaches that can be used for specific purposes, but ideally generic modular ALGs are developed that can be used for a wide range of scenarios, regions and model types. We have compiled features of such generic ALGs and propose a possible software architecture. Considerable joint efforts are required to develop such generic ALGs, but the benefits in terms of a better understanding and development of more efficient agricultural policies would be high.
Intraspecific trait variation (ITV) is thought to play a significant role in community assembly, but the magnitude and direction of its influence are not well understood. Although it may be critical to better explain population persistence, species interactions, and therefore biodiversity patterns, manipulating ITV in experiments is challenging. We therefore incorporated ITV into a trait‐ and individual‐based model of grassland community assembly by adding variation to the plants’ functional traits, which then drive life‐history tradeoffs. Varying the amount of ITV in the simulation, we examine its influence on pairwise‐coexistence and then on the species diversity in communities of different initial sizes. We find that ITV increases the ability of the weakest species to invade most, but that this effect does not scale to the community level, where the primary effect of ITV is to increase the persistence and abundance of the competitively‐average species. Diversity of the initial community is also of critical importance in determining ITV's efficacy; above a threshold of interspecific diversity, ITV does not increase diversity further. For communities below this threshold, ITV mainly helps to increase diversity in those communities that would otherwise be low‐diversity. These findings suggest that ITV actively maintains diversity by helping the species on the margins of persistence, but mostly in habitats of relatively low alpha and beta diversity.
Ecosystems respond in various ways to disturbances. Quantifying ecological stability therefore requires inspecting multiple stability properties, such as resistance, recovery, persistence and invariability. Correlations among these properties can reduce the dimensionality of stability, simplifying the study of environmental effects on ecosystems. A key question is how the kind of disturbance affects these correlations. We here investigated the effect of three disturbance types (random, species-specific, local) applied at four intensity levels, on the dimensionality of stability at the population and community level. We used previously parameterized models that represent five natural communities, varying in species richness and the number of trophic levels. We found that disturbance type but not intensity affected the dimensionality of stability and only at the population level. The dimensionality of stability also varied greatly among species and communities. Therefore, studying stability cannot be simplified to using a single metric and multi-dimensional assessments are still to be recommended.
Environmental factors shape the spatial distribution and dynamics of populations. Understanding how these factors interact with movement behavior is critical for efficient conservation, in particular for migratory species. Adult female green sea turtles, Chelonia mydas, migrate between foraging and nesting sites that are generally separated by thousands of kilometers. As an emblematic endangered species, green turtles have been intensively studied, with a focus on nesting, migration, and foraging. Nevertheless, few attempts integrated these behaviors and their trade‐offs by considering the spatial configurations of foraging and nesting grounds as well as environmental heterogeneity like oceanic currents and food distribution. We developed an individual‐based model to investigate the impact of local environmental conditions on emerging migratory corridors and reproductive output and to thereby identify conservation priority sites. The model integrates movement, nesting, and foraging behavior. Despite being largely conceptual, the model captured realistic movement patterns which confirm field studies. The spatial distribution of migratory corridors and foraging hot spots was mostly constrained by features of the regional landscape, such as nesting site locations, distribution of feeding patches, and oceanic currents. These constraints also explained the mixing patterns in regional forager communities. By implementing alternative decision strategies of the turtles, we found that foraging site fidelity and nesting investment, two characteristics of green turtles' biology, are favorable strategies under unpredictable environmental conditions affecting their habitats. Based on our results, we propose specific guidelines for the regional conservation of green turtles as well as future research suggestions advancing spatial ecology of sea turtles. Being implemented in an easy to learn open‐source software, our model can coevolve with the collection and analysis of new data on energy budget and movement into a generic tool for sea turtle research and conservation. Our modeling approach could also be useful for supporting the conservation of other migratory marine animals.
Anthropogenic pressures increasingly alter natural systems. Therefore, understanding the resilience of agent-based complex systems such as ecosystems, i.e. their ability to absorb these pressures and sustain their functioning and services, is a major challenge. However, the mechanisms underlying resilience are still poorly understood. A main reason for this is the multidimensionality of both resilience, embracing the three fundamental stability properties recovery, resistance and persistence, and of the specific situations for which stability properties can be assessed. Agent-based models (ABM) complement empirical research which is, for logistic reasons, limited in coping with these multiple dimensions. Besides their ability to integrate multidimensionality through extensive manipulation in a fully controlled system, ABMs can capture the emergence of system resilience from individual interactions and feedbacks across different levels of organization. To assess the extent to which this potential of ABMs has already been exploited, we reviewed the state of the art in exploring resilience and its multidimensionality in ecological and socio-ecological systems with ABMs. We found that the potential of ABMs is not utilized in most models, as they typically focus on a single dimension of resilience by using variability as a proxy for persistence, and are limited to one reference state, disturbance type and scale. Moreover, only few studies explicitly test the ability of different mechanisms to support resilience. To overcome these limitations, we recommend to simultaneously assess multiple stability properties for different situations and under consideration of the mechanisms that are hypothesised to render a system resilient. This will help us to better exploit the potential of ABMs to understand and quantify resilience mechanisms, and hence support solving real-world problems related to the resilience of agent-based complex systems.
Simulation models that describe autonomous individual organisms (individual based models, IBM) or agents (agent-based models, ABM) have become a widely used tool, not only in ecology, but also in many other disciplines dealing with complex systems made up of autonomous entities. However, there is no standard protocol for describing such simulation models, which can make them difficult to understand and to duplicate. This paper presents a proposed standard protocol, ODD, for describing IBMs and ABMs, developed and tested by 28 modellers who cover a wide range of fields within ecology. This protocol consists of three blocks (Overview, Design concepts, and Details), which are subdivided into seven elements: Purpose, State variables and scales, Process overview and scheduling, Design concepts, Initialization, Input, and Submodels. We explain which aspects of a model should be described in each element, and we present an example to illustrate the protocol in use. In addition, 19 examples are available in an Online Appendix. We consider ODD as a first step for establishing a more detailed common format of the description of IBMs and ABMs. Once initiated, the protocol will hopefully evolve as it becomes used by a sufficiently large proportion of modellers. (c) 2006 Elsevier B.V. All rights reserved.
Polychlorinated biphenyls (PCBs) can cause endocrine disruption, cancer, immunosuppression, or reproductive failure in animals. We used an individual-based model to explore whether and how PCB-associated reproductive failure could affect the dynamics of a hypothetical polar bear (Ursus maritimus) population exposed to PCBs to the same degree as the East Greenland subpopulation. Dose-response data from experimental studies on a surrogate species, the mink (Mustela vision), were used in the absence of similar data for polar bears. Two alternative types of reproductive failure in relation to maternal sum-PCB concentrations were considered: increased abortion rate and increased cub mortality. We found that the quantitative impact of PCB-induced reproductive failure on population growth rate depended largely on the actual type of reproductive failure involved. Critical potencies of the dose-response relationship for decreasing the population growth rate were established for both modeled types of reproductive failure. Comparing the model predictions of the age-dependent trend of sum-PCBs concentrations in females with actual field measurements from East Greenland indicated that it was unlikely that PCB exposure caused a high incidence of abortions in the subpopulation. However, on the basis of this analysis, it could not be excluded that PCB exposure contributes to higher cub mortality. Our results highlight the necessity for further research on the possible influence of PCBs on polar bear reproduction regarding their physiological pathway. This includes determining the exact cause of reproductive failure, i.e., in utero exposure versus lactational exposure of offspring; the timing of offspring death; and establishing the most relevant reference metrics for the dose-response relationship.
Population viability analysis (PVA) models are used to estimate population extinction risk under different scenarios. Both simple and complex PVA models are developed and have their specific pros and cons; the question therefore arises whether we always use the most appropriate model type. Generally, the specific purpose of a model and the availability of data are listed as determining the choice of model type, but this has not been formally tested yet. We quantified the relative importance of model purpose and nine metrics of data availability and resolution for the choice of a PVA model type, while controlling for effects of the different life histories of the modelled species. We evaluated 37 model pairs: each consisting of a generally simpler, population-based model (PBM) and a more complex, individual-based model (IBM) developed for the same species. The choice of model type was primarily affected by the availability and resolution of demographic, dispersal and spatial data. Low-resolution data resulted in the development of less complex models. Model purpose did not affect the choice of the model type. We confirm the general assumption that poor data availability is the main reason for the wide use of simpler models, which may have limited predictive power for population responses to changing environmental conditions. Conservation biology is a crisis discipline where researchers learned to work with the data at hand. However, for threatened and poorly-known species, there is no short-cut when developing either a PBM or an IBM: investments to collect appropriately detailed data are required to ensure PVA models can assess extinction risk under complex environmental conditions. (C) 2015 Elsevier B.V. All rights reserved.
Robustness analysis: Deconstructing computational models for ecological theory and applications
(2016)
The design of computational models is path-dependent: the choices made in each step during model development constrain the choices that are available in the subsequent steps. The actual path of model development can be extremely different, even for the same system, because the path depends on the question addressed, the availability of data, and the consideration of specific expert knowledge, in addition to the experience, background, and modelling preferences of the modellers. Thus, insights from different models are practically impossible to integrate, which hinders the development of general theory. We therefore suggest augmenting the current culture of communicating models as working just fine with a culture of presenting analyses in which we try to break models, i.e., model mechanisms explaining certain observations break down. We refer to the systematic attempts to break a model as “robustness analysis” (RA). RA is the systematic deconstruction of a model by forcefully changing the model's parameters, structure, and representation of processes. We discuss the nature and elements of RA and provide brief examples. RA cannot be completely formalized into specific techniques and instead corresponds to detective work that is driven by general questions and specific hypotheses, with strong attention focused on unusual behaviours. Both individual modellers and ecological modelling in general will benefit from RA because RA helps with understanding models and identifying “robust theories”, which are general principles that are independent of the idiosyncrasies of specific models. Integrating the results of RAs from different models to address certain systems or questions will then provide a comprehensive overview of when certain mechanisms control system behaviour and when and why this control ceases. This approach can provide insights into the mechanisms that lead to regime shifts in actual ecological systems.
Climate change and land use management practices are major drivers of biodiversity in terrestrial ecosystems. To understand and predict resulting changes in community structures, individual-based and spatially explicit population models are a useful tool but require detailed data sets for each species. More generic approaches are thus needed. Here we present a trait-based functional type approach to model savanna birds. The aim of our model is to explore the response of different bird functional types to modifications in habitat structure. The functional types are characterized by different traits, in particular body mass, which is related to life-history traits (reproduction and mortality) and spatial scales (home range area and dispersal ability), as well as the use of vegetation structures for foraging and nesting, which is related to habitat quality and suitability. We tested the performance of the functional types in artificial landscapes varying in shrub:grass ratio and clumping intensity of shrub patches. We found that an increase in shrub encroachment and a decrease in habitat quality caused by land use mismanagement and climate change endangered all simulated bird functional types. The strength of this effect was related to the preferred habitat. Furthermore, larger-bodied insectivores and omnivores were more prone to extinction due to shrub encroachment compared to small-bodied species. Insectivorous and omnivorous birds were more sensitive to clumping intensity of shrubs whereas herbivorous and carnivorous birds were most affected by a decreasing amount of grass cover. From an applied point of view, our findings emphasize that policies such as woody plant removal and a reduction in livestock stocking rates to prevent shrub encroachment should prioritize the enlargement of existing grassland patches. Overall, our results show that the combination of an individual-based modelling approach with carefully defined functional types can provide a powerful tool for exploring biodiversity responses to environmental changes. Furthermore, the increasing accumulation of worldwide data sets on species’ core and soft traits (surrogates to determine core traits indirectly) on one side and the refinement of conceptual frameworks for animal functional types on the other side will further improve functional type approaches which consider the sensitivities of multiple species to climate change, habitat loss, and fragmentation.
Current rates of environmental change are exceeding the capacity of many populations to adapt to new conditions and thus avoid demographic collapse and ultimate extinction. In particular, cold-water freshwater fish species are predicted to experience strong selective pressure from climate change and a wide range of interacting anthropogenic stressors in the near future. To implement effective management and conservation measures, it is crucial to quantify the maximum rate of change that cold-water freshwater fish populations can withstand. Here, we present a spatially explicit eco-genetic individual-based model, inSTREAM-Gen, to predict the eco-evolutionary dynamics of stream-dwelling trout under anthropogenic environmental change. The model builds on a well-tested demographic model, which includes submodels of river dynamics, bioenergetics, and adaptive habitat selection, with a new genetic module that allows exploration of genetic and life-history adaptations to new environments. The genetic module models the transmission of two key traits, size at emergence and maturity size threshold. We parameterized the model for a brown trout (Salmo trutta L.) population at the warmest edge of its range to validate it and analyze its sensitivity to parameters under contrasting thermal profiles. To illustrate potential applications of the model, we analyzed the population's demographic and evolutionary dynamics under scenarios of (1) climate change-induced warming, and (2) warming plus flow reduction resulting from climate and land use change, compared to (3) a baseline of no environmental change. The model predicted severe declines in density and biomass under climate warming. These declines were lower than expected at range margins because of evolution towards smaller size at both emergence and maturation compared to the natural evolution under the baseline conditions. Despite stronger evolutionary responses, declining rates were substantially larger under the combined warming and flow reduction scenario, leading to a high probability of population extinction over contemporary time frames. Therefore, adaptive responses could not prevent extinction under high rates of environmental change. Our model demonstrates critical elements of next generation ecological modelling aiming at predictions in a changing world as it accounts for spatial and temporal resource heterogeneity, while merging individual behaviour and bioenergetics with microevolutionary adaptations.
Both dispersal and local demographic processes determine a population's distribution among habitats of varying quality, yet most theory, experiments, and field studies have focused on the former. We use a generic model to show how both processes contribute to a population's distribution, and how the relative importance of each mechanism depends on scale. In contrast to studies only considering habitat-dependent dispersal, we show that predictions of ideal free distribution (IFD) theory are relevant even at landscape scales, where the assumptions of IFD theory are violated. This is because scales that inhibit one process, promote the other's ability to drive populations to the IFD. Furthermore, because multiple processes can generate IFDs, the pattern alone does not specify a causal mechanism. This is important because populations with IFDs generated by dispersal or demography respond much differently to shifts in resource distributions.
The causes underlying the increased mortality of honeybee Apis mellifera colonies observed over the past decade remain unclear. Since so far the evidence for monocausal explanations is equivocal, involvement of multiple stressors is generally assumed. We here focus on various aspects of forage availability, which have received less attention than other stressors because it is virtually impossible to explore them empirically. We applied the colony model BEEHAVE, which links within-hive dynamics and foraging, to stylized landscape settings to explore how foraging distance, forage supply, and “forage gaps”, i.e. periods in which honeybees cannot find any nectar and pollen, affect colony resilience and the mechanisms behind. We found that colony extinction was mainly driven by foraging distance, but the timing of forage gaps had strongest effects on time to extinction. Sensitivity to forage gaps of 15 days was highest in June or July even if otherwise forage availability was sufficient to survive. Forage availability affected colonies via cascading effects on queen's egg-laying rate, reduction of new-emerging brood stages developing into adult workers, pollen debt, lack of workforce for nursing, and reduced foraging activity. Forage gaps in July led to reduction in egg-laying and increased mortality of brood stages at a time when the queen's seasonal egg-laying rate is at its maximum, leading to colony failure over time. Our results demonstrate that badly timed forage gaps interacting with poor overall forage supply reduce honeybee colony resilience. Existing regulation mechanisms which in principle enable colonies to cope with varying forage supply in a given landscape and year, such as a reduction in egg-laying, have only a certain capacity. Our results are hypothetical, as they are obtained from simplified landscape settings, but they are consistent with existing empirical knowledge. They offer ample opportunities for testing the predicted effects of forage stress in controlled experiments.
Animal personality may affect an animal’s mobility in a given landscape, influencing its propensity to take risks in an unknown environment. We investigated the mobility of translocated common voles in two corridor systems 60 m in length and differing in width (1 m and 3 m). Voles were behaviorally phenotyped in repeated open field and barrier tests. Observed behavioral traits were highly repeatable and described by a continuous personality score. Subsequently, animals were tracked via an automated very high frequency (VHF) telemetry radio tracking system to monitor their movement patterns in the corridor system. Although personality did not explain movement patterns, corridor width determined the amount of time spent in the habitat corridor. Voles in the narrow corridor system entered the corridor faster and spent less time in the corridor than animals in the wide corridor. Thus, landscape features seem to affect movement patterns more strongly than personality. Meanwhile, site characteristics, such as corridor width, could prove to be highly important when designing corridors for conservation, with narrow corridors facilitating faster movement through landscapes than wider corridors.
Animal personality may affect an animal’s mobility in a given landscape, influencing its propensity to take risks in an unknown environment. We investigated the mobility of translocated common voles in two corridor systems 60 m in length and differing in width (1 m and 3 m). Voles were behaviorally phenotyped in repeated open field and barrier tests. Observed behavioral traits were highly repeatable and described by a continuous personality score. Subsequently, animals were tracked via an automated very high frequency (VHF) telemetry radio tracking system to monitor their movement patterns in the corridor system. Although personality did not explain movement patterns, corridor width determined the amount of time spent in the habitat corridor. Voles in the narrow corridor system entered the corridor faster and spent less time in the corridor than animals in the wide corridor. Thus, landscape features seem to affect movement patterns more strongly than personality. Meanwhile, site characteristics, such as corridor width, could prove to be highly important when designing corridors for conservation, with narrow corridors facilitating faster movement through landscapes than wider corridors.
Give chance a chance
(2019)
A large part of biodiversity theory is driven by the basic question of what allows species to coexist in spite of a confined number of niches. A substantial theoretical background to this question is provided by modern coexistence theory (MCT), which rests on mathematical approaches of invasion analysis to categorize underlying mechanisms into factors that reduce either niche overlap (stabilizing mechanisms) or the average fitness differences of species (equalizing mechanisms). While MCT has inspired biodiversity theory in the search for these underlying mechanisms, we feel that the strong focus on coexistence causes a bias toward the most abundant species and neglects the plethora of species that are less abundant and often show high local turnover. Given the more stochastic nature of their occurrence, we advocate a complementary cross-level approach that links individuals, small populations, and communities and explicitly takes into account (1) a more complete inclusion of environmental and demographic stochasticity affecting small populations, (2) intraspecific trait variation and behavioral plasticity, and (3) local heterogeneities, interactions, and feedbacks. Focusing on mechanisms that drive the temporary coviability of species rather than infinite coexistence, we suggest a new approach that could be dubbed coviability analysis (CVA). From a modeling perspective, CVA builds on the merged approaches of individual-based modeling and population viability analysis but extends them to the community level. From an empirical viewpoint, CVA calls for a stronger integration of spatiotemporal data on variability and noise, changing drivers, and interactions at the level of individuals. The resulting large volumes of data from multiple sources could be strongly supported by novel techniques tailored to the discovery of complex patterns in high-dimensional data. By complementing MCT through a stronger focus on the coviability of less common species, this approach can help make modern biodiversity theory more comprehensive, predictive, and relevant for applications.
Give chance a chance
(2019)
A large part of biodiversity theory is driven by the basic question of what allows species to coexist in spite of a confined number of niches. A substantial theoretical background to this question is provided by modern coexistence theory (MCT), which rests on mathematical approaches of invasion analysis to categorize underlying mechanisms into factors that reduce either niche overlap (stabilizing mechanisms) or the average fitness differences of species (equalizing mechanisms). While MCT has inspired biodiversity theory in the search for these underlying mechanisms, we feel that the strong focus on coexistence causes a bias toward the most abundant species and neglects the plethora of species that are less abundant and often show high local turnover. Given the more stochastic nature of their occurrence, we advocate a complementary cross-level approach that links individuals, small populations, and communities and explicitly takes into account (1) a more complete inclusion of environmental and demographic stochasticity affecting small populations, (2) intraspecific trait variation and behavioral plasticity, and (3) local heterogeneities, interactions, and feedbacks. Focusing on mechanisms that drive the temporary coviability of species rather than infinite coexistence, we suggest a new approach that could be dubbed coviability analysis (CVA). From a modeling perspective, CVA builds on the merged approaches of individual-based modeling and population viability analysis but extends them to the community level. From an empirical viewpoint, CVA calls for a stronger integration of spatiotemporal data on variability and noise, changing drivers, and interactions at the level of individuals. The resulting large volumes of data from multiple sources could be strongly supported by novel techniques tailored to the discovery of complex patterns in high-dimensional data. By complementing MCT through a stronger focus on the coviability of less common species, this approach can help make modern biodiversity theory more comprehensive, predictive, and relevant for applications.
Ecologists urgently need a better ability to predict how environmental change affects biodiversity. We examine individual-based ecology (IBE), a research paradigm that promises better a predictive ability by using individual-based models (IBMs) to represent ecological dynamics as arising from how individuals interact with their environment and with each other. A key advantage of IBMs is that the basis for predictions-fitness maximization by individual organisms-is more general and reliable than the empirical relationships that other models depend on. Case studies illustrate the usefulness and predictive success of long-term IBE programs. The pioneering programs had three phases: conceptualization, implementation, and diversification. Continued validation of models runs throughout these phases. The breakthroughs that make IBE more productive include standards for describing and validating IBMs, improved and standardized theory for individual traits and behavior, software tools, and generalized instead of system-specific IBMs. We provide guidelines for pursuing IBE and a vision for future IBE research.
The Southern Ocean ecosystem is characterized by extreme seasonal changes in environmental factors such as day length, sea ice extent and food availability. The key species Antarctic krill (Euphausia superba) has evolved metabolic and behavioural seasonal rhythms to cope with these seasonal changes. We investigate the switch between a physiological less active and active period for adult krill, a rhythm which seems to be controlled by internal biological clocks. These biological clocks can be synchronized by environmental triggers such as day length and food availability. They have evolved for particular environmental regimes to synchronize predictable seasonal environmental changes with important life cycle functions of the species. In a changing environment the time when krill is metabolically active and the time of peak food availability may not overlap if krill's seasonal activity is solely determined by photoperiod (day length). This is especially true for the Atlantic sector of the Southern Ocean where the spatio-temporal ice cover dynamics are changing substantially with rising average temperatures. We developed an individual-based model for krill to explore the impact of photoperiod and food availability on the growth and demographics of krill. We simulated dynamics of local krill populations (with no movement of krill assumed) along a south-north gradient for different triggers of metabolic activity and different levels of food availability below the ice. We also observed the fate of larval krill which cannot switch to low metabolism and therefore are likely to overwinter under ice. Krill could only occupy the southern end of the gradient, where algae bloom only lasts for a short time, when alternative food supply under the ice was high and metabolic activity was triggered by photoperiod. The northern distribution was limited by lack of overwintering habitat for krill larvae due to short duration of sea ice cover even for high food content under the ice. The variability of the krill's length-frequency distributions varied for different triggers of metabolic activity, but did not depend on the sea ice extent. Our findings suggest a southward shift of krill populations due to reduction in the spatial sea ice extent, which is consistent with field observations. Overall, our results highlight the importance of the explicit consideration of spatio-temporal sea ice dynamics especially for larval krill together with temporal synchronization through internal clocks, triggered by environmental factors (photoperiod and food) in adult krill for the population modelling of krill. (C) 2015 Elsevier B.V. All rights reserved.
There are two major limitations to the potential of computational models in ecology for producing general insights: their design is path-dependent, reflecting different underlying questions, assumptions, and data, and there is too little robustness analysis exploring where the model mechanisms explaining certain observations break down. We here argue that both limitations could be overcome if modellers in ecology would more often replicate existing models, try to break the models, and explore modifications. Replication comprises the re-implementation of an existing model and the replication of its results. Breaking models means to identify under what conditions the mechanisms represented in a model can no longer explain observed phenomena. The benefits of replication include less effort being spent to enter the iterative stage of model development and having more time for systematic robustness analysis. A culture of replication would lead to increased credibility, coherence and efficiency of computational modelling and thereby facilitate theory development.
Population models in ecology are often not good at predictions, even if they are complex and seem to be realistic enough. The reason for this might be that Occam's razor, which is key for minimal models exploring ideas and concepts, has been too uncritically adopted for more realistic models of systems. This can tic models too closely to certain situations, thereby preventing them from predicting the response to new conditions. We therefore advocate a new kind of parsimony to improve the application of Occam's razor. This new parsimony balances two contrasting strategies for avoiding errors in modeling: avoiding inclusion of nonessential factors (false inclusions) and avoiding exclusion of sometimes-important factors (false exclusions). It involves a synthesis of traditional modeling and analysis, used to describe the essentials of mechanistic relationships, with elements that arc included in a model because they have been reported to be or can arguably be assumed to be important under certain conditions. The resulting models should be able to reflect how the internal organization of populations change and thereby generate representations of the novel behavior necessary for complex predictions, including regime shifts.
The importance of a careful choice of the appropriate scale for studying ecological phenomena has been stressed repeatedly. However, issues of spatial scale in metapopulation dynamics received much more attention compared to temporal scale. Moreover, multiple calls were made to carefully choose the appropriate model structure for Population Viability Analysis (PVA). We assessed the effect of using coarser resolution in time and model structure on population dynamics. For this purpose, we compared outcomes of two PVA models differing in their time step: daily individual-based model (dIBM) and yearly stage-based model (ySBM), loaded with empirical data on a well-known metapopulation of the butterfly Boloria eunomia. Both models included the same environmental drivers of population dynamics that were previously identified as being the most important for this species. Under temperature change scenarios, both models yielded the same qualitative scenario ranking, but they quite substantially differed quantitatively with dIBM being more pessimistic in absolute viability measures. We showed that these differences stemmed from inter-individual heterogeneity in dIBM allowing for phenological shifts of individual appearance. We conclude that a finer temporal resolution and an individual-based model structure allow capturing the essential mechanisms necessary to go beyond mere PVA scenario ranking. We encourage researchers to carefully chose the temporal resolution and structure of their model aiming at (1) depicting the processes important for (meta)population dynamics of the species and (2) implementing the environmental change scenarios expected for their study system in the future, using the temporal resolution at which such changes are predicted to operate.
BEEHAVE offers a valuable tool for researchers to design and focus field experiments, for regulators to explore the relative importance of stressors to devise management and policy advice and for beekeepers to understand and predict varroa dynamics and effects of management interventions. We expect that scientists and stakeholders will find a variety of applications for BEEHAVE, stimulating further model development and the possible inclusion of other stressors of potential importance to honeybee colony dynamics.
The potential of ecological models for supporting environmental decision making is increasingly acknowledged. However, it often remains unclear whether a model is realistic and reliable enough. Good practice for developing and testing ecological models has not yet been established. Therefore, TRACE, a general framework for documenting a model's rationale, design, and testing was recently suggested. Originally TRACE was aimed at documenting good modelling practice. However, the word 'documentation' does not convey TRACE's urgency. Therefore, we re-define TRACE as a tool for planning, performing, and documenting good modelling practice. TRACE documents should provide convincing evidence that a model was thoughtfully designed, correctly implemented, thoroughly tested, well understood, and appropriately used for its intended purpose. TRACE documents link the science underlying a model to its application, thereby also linking modellers and model users, for example stakeholders, decision makers, and developers of policies. We report on first experiences in producing TRACE documents. We found that the original idea underlying TRACE was valid, but to make its use more coherent and efficient, an update of its structure and more specific guidance for its use are needed. The updated TRACE format follows the recently developed framework of model 'evaludation': the entire process of establishing model quality and credibility throughout all stages of model development, analysis, and application. TRACE thus becomes a tool for planning, documenting, and assessing model evaludation, which includes understanding the rationale behind a model and its envisaged use. We introduce the new structure and revised terminology of TRACE and provide examples. (C) 2014 Elsevier B.V. All rights reserved.
Confusion about model validation is one of the main challenges in using ecological models for decision support, such as the regulation of pesticides. Decision makers need to know whether a model is a sufficiently good representation of its real counterpart and what criteria can be used to answer this question. Unclear terminology is one of the main obstacles to a good understanding of what model validation is, how it works, and what it can deliver. Therefore, we performed a literature review and derived a standard set of terms. 'Validation' was identified as a catch-all term, which is thus useless for any practical purpose. We introduce the term 'evaludation', a fusion of 'evaluation' and 'validation', to describe the entire process of assessing a model's quality and reliability. Considering the iterative nature of model development, the modelling cycle, we identified six essential elements of evaludation: (i) 'data evaluation' for scrutinising the quality of numerical and qualitative data used for model development and testing; (ii) 'conceptual model evaluation' for examining the simplifying assumptions underlying a model's design; (iii) 'implementation verification' for testing the model's implementation in equations and as a computer programme; (iv) 'model output verification' for comparing model output to data and patterns that guided model design and were possibly used for calibration; (v) 'model analysis' for exploring the model's sensitivity to changes in parameters and process formulations to make sure that the mechanistic basis of main behaviours of the model has been well understood; and (vi) 'model output corroboration' for comparing model output to new data and patterns that were not used for model development and parameterisation. Currently, most decision makers require 'validating' a model by testing its predictions with new experiments or data. Despite being desirable, this is neither sufficient nor necessary for a model to be useful for decision support. We believe that the proposed set of terms and its relation to the modelling cycle can help to make quality assessments and reality checks of ecological models more comprehensive and transparent. (C) 2013 Elsevier B.V. All rights reserved.
Current chemical risk assessment procedures may result in imprecise estimates of risk due to sometimes arbitrary simplifying assumptions. As a way to incorporate ecological complexity and improve risk estimates, mechanistic effect models have been recommended. However, effect modeling has not yet been extensively used for regulatory purposes, one of the main reasons being uncertainty about which model type to use to answer specific regulatory questions. We took an individual-based model (IBM), which was developed for risk assessment of soil invertebrates and includes avoidance of highly contaminated areas, and contrasted it with a simpler, more standardized model, based on the generic metapopulation matrix model RAMAS. In the latter the individuals within a sub-population are not treated as separate entities anymore and the spatial resolution is lower. We explored consequences of model aggregation in terms of assessing population-level effects for different spatial distributions of a toxic chemical. For homogeneous contamination of the soil, we found good agreement between the two models, whereas for heterogeneous contamination, at different concentrations and percentages of contaminated area, RAMAS results were alternatively similar to IBM results with and without avoidance, and different food levels. This inconsistency is explained on the basis of behavioral responses that are included in the IBM but not in RAMAS. Overall, RAMAS was less sensitive than the IBM in detecting population-level effects of different spatial patterns of exposure. We conclude that choosing the right model type for risk assessment of chemicals depends on whether or not population-level effects of small-scale heterogeneity in exposure need to be detected. We recommend that if in doubt, both model types should be used and compared. Describing both models following the same standard format, the ODD protocol, makes them equally transparent and understandable. The simpler model helps to build up trust for the more complex model and can be used for more homogeneous exposure patterns. The more complex model helps detecting and understanding the limitations of the simpler model and is needed to ensure ecological realism for more complex exposure scenarios. (C) 2013 Elsevier B.V. All rights reserved.
Current environmental risk assessment (ERA) of chemicals for aquatic invertebrates relies on standardized laboratory tests in which toxicity effects on individual survival, growth and reproduction are measured. Such tests determine the threshold concentration of a chemical below which no population-level effects are expected. How well this procedure captures effects on individuals and populations, however, remains an open question. Here we used mechanistic effect models, combining individual-level reproduction and survival models with an individual-based population model (IBM), to understand the individuals' responses and extrapolate them to the population level. We used a toxicant (Dispersogen A) for which adverse effects on laboratory populations were detected at the determined threshold concentration and thus challenged the conservatism of the current risk assessment method. Multiple toxicity effects on reproduction and survival were reported, in addition to effects on the F1 generation. We extrapolated commonly tested individual toxicity endpoints, reproduction and survival, to the population level using the IBM. Effects on reproduction were described via regression models. To select the most appropriate survival model, the IBM was run assuming either stochastic death (SD) or individual tolerance (IT). Simulations were run for different scenarios regarding the toxicant's effects: survival toxicity, reproductive toxicity, or survival and reproductive toxicity. As population-level endpoints, we used population size and structure and extinction risk. We found that survival represented as SD explained population dynamics better than IT. Integrating toxicity effects on both reproduction and survival yielded more accurate predictions of population effects than considering isolated effects. To fully capture population effects observed at high toxicant concentrations, toxicity effects transmitted to the F1 generation had to be integrated. Predicted extinction risk was highly sensitive to the assumptions about individual-level effects. Our results demonstrate that the endpoints used in current standard tests may not be sufficient for assessing the risk of adverse effects on populations. A combination of laboratory population experiments with mechanistic effect models is a powerful tool to better understand and predict effects on both individuals and populations. Mechanistic effect modelling thus holds great potential to improve the accuracy of ERA of chemicals in the future. (C) 2013 The Authors. Published by Elsevier B.V. All rights reserved.
The relative contribution of personal and social information to explain individual and collective behavior in different species and contexts is an open question in animal ecology. In particular, there is a major lack of studies combining theoretical and empirical approaches to test the relative relevance of different hypothesized individual behaviors to predict empirical collective patterns. We used an individual-based model to confront three hypotheses about the information transfer between social scavengers (Griffon Vultures, Gyps fulvus) when searching for carrion: (1) Vultures only use personal information during foraging ("nonsocial" hypothesis); (2) they create long chains of vultures by following both other vultures that are flying towards carcasses and vultures that are following other vultures that are flying towards carcasses ("chains of vultures" hypothesis); and (3) vultures are only attracted by other vultures that are sinking vertically to a carcass ("local enhancement" hypothesis). The chains of vultures hypothesis has been used in existing models, but never been confronted with field data. Testing is important, though, because these hypotheses could have different management implications. The model was parameterized to mimic the behavior and the densities of both Griffon Vultures and carcasses in a 10 000-km(2) study area in northeastern Spain. We compared the number of vultures attending simulated carcasses with those attending 25 continuously monitored experimental carcasses in the field. Social hypotheses outperformed the nonsocial hypothesis. The chains of vultures hypothesis overestimated the number of vultures feeding on carcasses; the local enhancement hypothesis fitted closely to the empirical data. Supported by our results, we discuss mechanistic and adaptive considerations that reveal that local enhancement may be the key social mechanism behind collective foraging in this and likely other avian scavengers and/or social birds. It also highlights the current need for more studies confronting alternative models of key behaviors with empirical patterns in order to understand how collective behavior emerges in animal societies.
Females may select a mate based on signalling traits that are believed to accurately correlate with heritable aspects of male quality. Anthropogenic actions, in particular chemicals released into the environment, are now disrupting the accuracy of mating signals to convey information about male quality. The long-term prediction for disrupted mating signals is most commonly loss of female preference. Yet, this prediction has rarely been tested using quantitative models. We use agent-based models to explore the effects of rapid disruption of mating signals. In our model, a gene determines survival. Males signal their level of genetic quality via a signal trait, which females use to select a mate. We allowed this system of sexual selection to become established, before introducing a disruption between the male signal trait and quality, which was similar in nature to that induced by exogenous chemicals. Finally, we assessed the capacity of the system to recover from this disruption. We found that within a relatively short time frame, disruption of mating signals led to a lasting loss of female preference. Decreases in mean viability at the population-level were also observed, because sexual-selection acting against newly arising deleterious mutations was relaxed. The ability of the population to recover from disrupted mating signals was strongly influenced by the mechanisms that promoted or maintained genetic diversity in traits under sexual selection. Our simple model demonstrates that environmental perturbations to the accuracy of male mating signals can result in a long-term loss of female preference for those signals within a few generations. What is more, the loss of this preference can have knock-on consequences for mean population fitness.