Parameter Sensitivity of Transient Community Dynamics

AbstractTransient dynamics have always intrigued ecologists, but current rapid environmental change (inducing transients even in previously undisturbed systems) has highlighted their importance more than ever. Here, I introduce a method for analyzing the sensitivity of transient ecological dynamics to parameter perturbations. The question the method answers is: how would the community dynamics have unfolded for some time horizon had the parameters been slightly different? I apply the method to three empirically parameterized models: competition between native forbs and exotic grasses in California, a host-parasitoid system, and an experimental chemostat predator-prey model. These applications showcase the ecological insights one can gain from models using transient sensitivity analysis. First, one can find parameters and their combinations whose perturbations disproportionately affect a system. Second, one can identify particular windows of time during which the predicted deviation from the unperturbed trajectories is especially large and utilize this information for management purposes. Third, there is an inverse relationship between transient and long-term sensitivities whenever the interacting populations are ecologically similar; paradoxically, the smaller the immediate response of the system, the more extreme its long-term response will be.

No sensitivity to functional forms in the Rosenzweig-MacArthur model with strong environmental stochasticity

The classic Rosenzweig-MacArthur predator-prey model has been shown to exhibit, like other coupled nonlinear ordinary differential equations (ODEs) from ecology, worrying sensitivity to model structure. This sensitivity manifests as markedly different community dynamics arising from saturating functional responses with nearly identical shapes but different mathematical expressions. Using a stochastic differential equation (SDE) version of the Rosenzweig-MacArthur model with the three functional responses considered by Fussmann & Blasius (2005), I show that such sensitivity seems to be solely a property of ODEs or stochastic systems with weak noise. SDEs with strong environmental noise have by contrast very similar fluctuations patterns, irrespective of the mathematical formula used. Although eigenvalues of linearized predator-prey models have been used as an argument for structural sensitivity, they can also be an argument against structural sensitivity. While the sign of the eigenvalues' real part is sensitive to model structure, its magnitude and the presence of imaginary parts are not, which suggests noise-driven oscillations for a broad range of carrying capacities. I then discuss multiple other ways to evaluate structural sensitivity in a stochastic setting, for predator-prey or other ecological systems.

Direct input of monitoring data into a mechanistic ecological model as a way to identify the phytoplankton growth-rate response to temperature variations

We present an approach (knowledge-and-data-driven, KDD, modeling) that allows us to get closer to understanding the processes that affect the dynamics of plankton communities. This approach, based on the use of time series obtained as a result of ecosystem monitoring, combines the key features of both the knowledge-driven modeling (mechanistic models) and data-driven (DD) modeling. Using a KDD model, we reveal the phytoplankton growth-rate fluctuations in the ecosystem of the Naroch Lakes and determine the degree of phase synchronization between fluctuations in the phytoplankton growth rate and temperature variations. More specifically, we estimate a numerical value of the phase locking index (PLI), which allows us to assess how temperature fluctuations affect the dynamics of phytoplankton growth rates. Since, within the framework of KDD modeling, we directly include the time series obtained as a result of field measurements in the model equations, the dynamics of the phytoplankton growth rate obtained from the KDD model reflect the behavior of the lake ecosystem as a whole, and PLI can be considered as a holistic parameter.

Bridging Theories for Ecosystem Stability Through Structural Sensitivity Analysis of Ecological Models in Equilibrium

Ecologists are challenged by the need to bridge and synthesize different approaches and theories to obtain a coherent understanding of ecosystems in a changing world. Both food web theory and regime shift theory shine light on mechanisms that confer stability to ecosystems, but from different angles. Empirical food web models are developed to analyze how equilibria in real multi-trophic ecosystems are shaped by species interactions, and often include linear functional response terms for simple estimation of interaction strengths from observations. Models of regime shifts focus on qualitative changes of equilibrium points in a slowly changing environment, and typically include non-linear functional response terms. Currently, it is unclear how the stability of an empirical food web model, expressed as the rate of system recovery after a small perturbation, relates to the vulnerability of the ecosystem to collapse. Here, we conduct structural sensitivity analyses of classical consumer-resource models in equilibrium along an environmental gradient. Specifically, we change non-proportional interaction terms into proportional ones, while maintaining the equilibrium biomass densities and material flux rates, to analyze how alternative model formulations shape the stability properties of the equilibria. The results reveal no consistent relationship between the stability of the original models and the proportionalized versions, even though they describe the same biomass values and material flows. We use these findings to critically discuss whether stability analysis of observed equilibria by empirical food web models can provide insight into regime shift dynamics, and highlight the challenge of bridging alternative modelling approaches in ecology and beyond.

Host adaptation to novel pathogen introduction: Predicting conditions that promote evolutionary rescue

Novel pathogen introduction can have drastic consequences for naive host populations, and outcomes can be difficult to predict. Evolutionary rescue (ER) provides a foundation for understanding whether hosts are driven to extinction or survive via adaptation. Currently, patterns of host population dynamics alongside evidence of adaptation are used to infer ER. However, the gap between established ER theory and complexity inherent in natural systems makes interpreting empirical patterns difficult because they can be confounded with ecological drivers of survival under current theory. To bridge this gap, we expand ER theory to include biological selective agents, such as pathogens. We find birth processes to be more important than previously theorised in determining ER potential. We employ a novel framework evaluating ER potential within natural systems and gain ability to identify system characteristics that make ER possible. Identifying these characteristics allows a shift from retrospective observation to a predictive mindset, and our findings suggest that ER occurrence may be more limited than previously thought. We use the plague system of Yersinia pestis infecting Cynomys ludovicianus (black-tailed prairie dogs) and Spermophilus beecheyi (California ground squirrels) as a case study.

Hidden layers of density dependence in consumer feeding rates

Functional responses relate a consumer's feeding rates to variation in its abiotic and biotic environment, providing insight into consumer behaviour and fitness, and underpinning population and food-web dynamics. Despite their broad relevance and long-standing history, we show here that the types of density dependence found in classic resource- and consumer-dependent functional-response models equate to strong and often untenable assumptions about the independence of processes underlying feeding rates. We first demonstrate mathematically how to quantify non-independence between feeding and consumer interference and between feeding on multiple resources. We then analyse two large collections of functional-response data sets to show that non-independence is pervasive and borne out in previously hidden forms of density dependence. Our results provide a new lens through which to view variation in consumer feeding rates and disentangle the biological underpinnings of species interactions in multi-species contexts.

Identifying the sources of structural sensitivity in partially specified biological models

Biological systems are characterised by a high degree of uncertainty and complexity, which implies that exact mathematical equations to describe biological processes cannot generally be justified. Moreover, models can exhibit sensitivity to the precise formulations of their component functions—a property known as structural sensitivity. Structural sensitivity can be revealed and quantified by considering partially specified models with uncertain functions, but this goes beyond well-established, parameter-based sensitivity analysis, and currently presents a mathematical challenge. Here we build upon previous work in this direction by addressing the crucial question of identifying the processes which act as the major sources of model uncertainty and those which are less influential. To achieve this goal, we introduce two related concepts: (1) the gradient of structural sensitivity, accounting for errors made in specifying unknown functions, and (2) the partial degree of sensitivity with respect to each function, a global measure of the uncertainty due to possible variation of the given function while the others are kept fixed. We propose an iterative framework of experiments and analysis to inform a heuristic reduction of structural sensitivity in a model. To demonstrate the framework introduced, we investigate the sources of structural sensitivity in a tritrophic food chain model.

Exploring the role of spatial and stoichiometric heterogeneity in the top-down control in eutrophic planktonic ecosystems

Understanding the impact of eutrophication on the dynamics of aquatic food webs, remains a long-term challenge in ecology. Mathematical models generally predict the destabilisation of such webs, under increasing eutrophication levels, with large oscillations of species densities that can result in their extinction. This is at odds with a number of ecological observations that show stable dynamics even for high nutrient loads. The apparent discrepancy between theory and observations is known as the Rosenzweig's 'paradox of enrichment' and various solutions to the problem have been proposed over the years. In this study, we explore the stabilisation of dynamics of a tri-trophic plankton model in a eutrophic environment which occurs as a result of interplay of space heterogeneity, ecological stoichiometry, and food taxis of predators. We build a variety of models of increasing complexity, to explore various scenarios of phytoplankton growth, zooplankton food-dependent vertical movement, and different stoichiometric limitations of zooplankton. We show that the synergy among the vertical gradient in phytoplankton growth, phytoplankton structuring in terms of their stoichiometric ratio, and food-dependent vertical movement of zooplankton, would result in a postponing of destabilisation of eutrophic systems as compared to a well-mixed system. Our approach reveals a high complexity of the bifurcation structure of the system when key model parameters, such as the degree of eutrophication and light shading, are varied. We find coexistence of limit cycles and stable equilibria and that the possibility of multiple attractors in the system can result in hysteresis phenomena when the nutrient load is manipulated. These results are relevant and should be taken into account in lake restoration programs.

Long transients in ecology: Theory and applications

This paper discusses the recent progress in understanding the properties of transient dynamics in complex ecological systems. Predicting long-term trends as well as sudden changes and regime shifts in ecosystems dynamics is a major issue for ecology as such changes often result in population collapse and extinctions. Analysis of population dynamics has traditionally been focused on their long-term, asymptotic behavior whilst largely disregarding the effect of transients. However, there is a growing understanding that in ecosystems the asymptotic behavior is rarely seen. A big new challenge for theoretical and empirical ecology is to understand the implications of long transients. It is believed that the identification of the corresponding mechanisms along with the knowledge of scaling laws of the transient's lifetime should substantially improve the quality of long-term forecasting and crisis anticipation. Although transient dynamics have received considerable attention in physical literature, research into ecological transients is in its infancy and systematic studies are lacking. This text aims to partially bridge this gap and facilitate further progress in quantitative analysis of long transients in ecology. By revisiting and critically examining a broad variety of mathematical models used in ecological applications as well as empirical facts, we reveal several main mechanisms leading to the emergence of long transients and hence lays the basis for a unifying theory.

The mechanistic basis for higher-order interactions and non-additivity in competitive communities

Motivated by both analytical tractability and empirical practicality, community ecologists have long treated the species pair as the fundamental unit of study. This notwithstanding, the challenge of understanding more complex systems has repeatedly generated interest in the role of so-called higher-order interactions (HOIs) imposed by species beyond the focal pair. Here we argue that HOIs - defined as non-additive effects of density on per capita growth - are best interpreted as emergent properties of phenomenological models (e.g. Lotka-Volterra competition) rather than as distinct 'ecological processes' in their own right. Using simulations of consumer-resource models, we explore the mechanisms and system properties that give rise to HOIs in observational data. We demonstrate that HOIs emerge under all but the most restrictive of assumptions, and that incorporating non-additivity into phenomenological models improves the quantitative and qualitative accuracy of model predictions. Notably, we also observe that HOIs derive primarily from mechanisms and system properties that apply equally to single-species or pairwise systems as they do to more diverse communities. Consequently, there exists a strong mandate for further recognition of non-additive effects in both theoretical and empirical research.

Community dynamics and sensitivity to model structure: towards a probabilistic view of process-based model predictions

Statistical inference and mechanistic, process-based modelling represent two philosophically different streams of research whose primary goal is to make predictions. Here, we merge elements from both approaches to keep the theoretical power of process-based models while also considering their predictive uncertainty using Bayesian statistics. In environmental and biological sciences, the predictive uncertainty of process-based models is usually reduced to parametric uncertainty. Here, we propose a practical approach to tackle the added issue of structural sensitivity, the sensitivity of predictions to the choice between quantitatively close and biologically plausible models. In contrast to earlier studies that presented alternative predictions based on alternative models, we propose a probabilistic view of these predictions that include the uncertainty in model construction and the parametric uncertainty of each model. As a proof of concept, we apply this approach to a predator-prey system described by the classical Rosenzweig-MacArthur model, and we observe that parametric sensitivity is regularly overcome by structural sensitivity. In addition to tackling theoretical questions about model sensitivity, the proposed approach can also be extended to make probabilistic predictions based on more complex models in an operational context. Both perspectives represent important steps towards providing better model predictions in biology, and beyond.

Temperature-dependent virus lifecycle choices may reveal and predict facets of the biology of opportunistic pathogenic bacteria

Melioidosis, a serious illness caused by Burkholderia pseudomallei, results in up to 40% fatality in infected patients. The pathogen is found in tropical water and soil. Recent findings demonstrated that bacterial numbers can be regulated by a novel clade of phages that are abundant in soil and water. These phages differentially infect their bacterial hosts causing lysis at high temperatures and lysogeny at lower temperatures. Thus seasonal and daily temperature variations would cause switches in phage-bacteria interactions. We developed mathematical models using realistic parameters to explore the impact of phages on B. pseudomallei populations in the surface water of rice fields over time and under seasonally changing environmental conditions. Historical records were used to provide UV radiation levels and temperature for two Thailand provinces. The models predict seasonal variation of phage-free bacterial numbers correlates with the higher risk of melioidosis acquisition during the “warm and wet” season. We find that enrichment of the environment may lead to irregular large amplitude pulses of bacterial numbers that could significantly increase the probability of disease acquisition. Our results suggest that the phages may regulate B. pseudomallei populations throughout the seasons, and these data can potentially help improve the melioidosis prevention efforts in Southeast Asia.

Is structural sensitivity a problem of oversimplified biological models? Insights from nested Dynamic Energy Budget models

Many current issues in ecology require predictions made by mathematical models, which are built on somewhat arbitrary choices. Their consequences are quantified by sensitivity analysis to quantify how changes in model parameters propagate into an uncertainty in model predictions. An extension called structural sensitivity analysis deals with changes in the mathematical description of complex processes like predation. Such processes are described at the population scale by a specific mathematical function taken among similar ones, a choice that can strongly drive model predictions. However, it has only been studied in simple theoretical models. Here, we ask whether structural sensitivity is a problem of oversimplified models. We found in predator-prey models describing chemostat experiments that these models are less structurally sensitive to the choice of a specific functional response if they include mass balance resource dynamics and individual maintenance. Neglecting these processes in an ecological model (for instance by using the well-known logistic growth equation) is not only an inappropriate description of the ecological system, but also a source of more uncertain predictions.

Does a better model yield a better argument? An info-gap analysis

Theories, models and computations underlie reasoned argumentation in many areas. The possibility of error in these arguments, though of low probability, may be highly significant when the argument is used in predicting the probability of rare high-consequence events. This implies that the choice of a theory, model or computational method for predicting rare high-consequence events must account for the probability of error in these components. However, error may result from lack of knowledge or surprises of various sorts, and predicting the probability of error is highly uncertain. We show that the putatively best, most innovative and sophisticated argument may not actually have the lowest probability of error. Innovative arguments may entail greater uncertainty than more standard but less sophisticated methods, creating an innovation dilemma in formulating the argument. We employ info-gap decision theory to characterize and support the resolution of this problem and present several examples.

Towards developing a general framework for modelling vertical migration in zooplankton

Diel vertical migration (DVM) of zooplankton is a widespread phenomenon in both oceans and lakes, and is generally considered to be the largest synchronized movement of biomass on Earth. Most existing mathematical models of DVM are based on the assumption that animals maximize a certain criterion such as the expected reproductive value, the venturous revenue, the ratio of energy gain/mortality or some predator avoidance function when choosing their instantaneous depth. The major shortcoming of this general point of view is that the predicted DVM may be strongly affected by a subjective choice of a particular optimization criterion. Here we argue that the optimal strategy of DVM can be unambiguously obtained as an outcome of selection in the underlying equations of genotype/traits frequency dynamics. Using this general paradigm, we explore the optimal strategy for the migration across different depths by zooplankton grazers throughout the day. To illustrate our ideas we consider four generic DVM models, each making different assumptions on the population dynamics of zooplankton, and demonstrate that in each model we need to maximize a particular functional to find the optimal strategy. Surprisingly, patterns of DVM obtained for different models greatly differ in terms of their parameters dependence. We then show that the infinite dimensional trait space of different zooplankton trajectories can be projected onto a low dimensional space of generalized parameters and the genotype evolution dynamics can be easily followed using this low-dimensional space. Using this space of generalized parameters we explore the influence of mutagenesis on evolution of DVM, and we show that strong mutagenesis allows the coexistence of an infinitely large number of strategies whereas for weak mutagenesis the selection results in the extinction of most strategies, with the surviving strategies all staying close to the optimal strategy in the corresponding mutagenesis-free system.

Quantifying uncertainty in partially specified biological models: how can optimal control theory help us?

Mathematical models in biology are highly simplified representations of a complex underlying reality and there is always a high degree of uncertainty with regards to model function specification. This uncertainty becomes critical for models in which the use of different functions fitting the same dataset can yield substantially different predictions-a property known as structural sensitivity. Thus, even if the model is purely deterministic, then the uncertainty in the model functions carries through into uncertainty in model predictions, and new frameworks are required to tackle this fundamental problem. Here, we consider a framework that uses partially specified models in which some functions are not represented by a specific form. The main idea is to project infinite dimensional function space into a low-dimensional space taking into account biological constraints. The key question of how to carry out this projection has so far remained a serious mathematical challenge and hindered the use of partially specified models. Here, we propose and demonstrate a potentially powerful technique to perform such a projection by using optimal control theory to construct functions with the specified global properties. This approach opens up the prospect of a flexible and easy to use method to fulfil uncertainty analysis of biological models.

Long-term transients and complex dynamics of a stage-structured population with time delay and the Allee effect

Traditionally, mathematical modeling in population ecology is mainly focused on asymptotic behavior of the model, i.e. as given by the system attractors. Recently, however, transient regimes and especially long-term transients have been recognized as playing a crucial role in the dynamics of ecosystems. In particular, long-term transients are a potential explanation of ecological regime shifts, when an apparently healthy population suddenly collapses and goes extinct. In this paper, we show that the interplay between delay in maturation and a strong Allee effect can result in long-term transients in a single species system. We first derive a simple 'conceptual' model of the population dynamics that incorporates both a strong Allee effect and maturation delay. Unlike much of the previous work, our approach is not empirical since our model is derived from basic principles. We show that the model exhibits a high complexity in its asymptotic dynamics including multi-periodic and chaotic attractors. We then show the existence of long-term transient dynamics in the system, when the population size oscillates for a long time between locally stable stationary states before it eventually settles either at the persistence equilibrium or goes extinct. The parametric space of the model is found to have a complex structure with the basins of attraction corresponding to the persistence and extinction states being of a complicated shape. This impedes the prediction of the eventual fate of the population, as a small variation in the maturation delay or the initial population size can either bring the population to extinction or ensure its persistence.

Bifurcation analysis of models with uncertain function specification: how should we proceed?

When we investigate the bifurcation structure of models of natural phenomena, we usually assume that all model functions are mathematically specified and that the only existing uncertainty is with respect to the parameters of these functions. In this case, we can split the parameter space into domains corresponding to qualitatively similar dynamics, separated by bifurcation hypersurfaces. On the other hand, in the biological sciences, the exact shape of the model functions is often unknown, and only some qualitative properties of the functions can be specified: mathematically, we can consider that the unknown functions belong to a specific class of functions. However, the use of two different functions belonging to the same class can result in qualitatively different dynamical behaviour in the model and different types of bifurcation. In the literature, the conventional way to avoid such ambiguity is to narrow the class of unknown functions, which allows us to keep patterns of dynamical behaviour consistent for varying functions. The main shortcoming of this approach is that the restrictions on the model functions are often given by cumbersome expressions and are strictly model-dependent: biologically, they are meaningless. In this paper, we suggest a new framework (based on the ODE paradigm) which allows us to investigate deterministic biological models in which the mathematical formulation of some functions is unspecified except for some generic qualitative properties. We demonstrate that in such models, the conventional idea of revealing a concrete bifurcation structure becomes irrelevant: we can only describe bifurcations with a certain probability. We then propose a method to define the probability of a bifurcation taking place when there is uncertainty in the parameterisation in our model. As an illustrative example, we consider a generic predator-prey model where the use of different parameterisations of the logistic-type prey growth function can result in different dynamics in terms of the type of the Hopf bifurcation through which the coexistence equilibrium loses stability. Using this system, we demonstrate a framework for evaluating the probability of having a supercritical or subcritical Hopf bifurcation.

Defining and detecting structural sensitivity in biological models: developing a new framework

When we construct mathematical models to represent biological systems, there is always uncertainty with regards to the model specification--whether with respect to the parameters or to the formulation of model functions. Sometimes choosing two different functions with close shapes in a model can result in substantially different model predictions: a phenomenon known in the literature as structural sensitivity, which is a significant obstacle to improving the predictive power of biological models. In this paper, we revisit the general definition of structural sensitivity, compare several more specific definitions and discuss their usefulness for the construction and analysis of biological models. Then we propose a general approach to reveal structural sensitivity with regards to certain system properties, which considers infinite-dimensional neighbourhoods of the model functions: a far more powerful technique than the conventional approach of varying parameters for a fixed functional form. In particular, we suggest a rigorous method to unearth sensitivity with respect to the local stability of systems' equilibrium points. We present a method for specifying the neighbourhood of a general unknown function with [Formula: see text] inflection points in terms of a finite number of local function properties, and provide a rigorous proof of its completeness. Using this powerful result, we implement our method to explore sensitivity in several well-known multicomponent ecological models and demonstrate the existence of structural sensitivity in these models. Finally, we argue that structural sensitivity is an important intrinsic property of biological models, and a direct consequence of the complexity of the underlying real systems.

Feeding on multiple sources: towards a universal parameterization of the functional response of a generalist predator allowing for switching

Understanding of complex trophic interactions in ecosystems requires correct descriptions of the rate at which predators consume a variety of different prey species. Field and laboratory data on multispecies communities are rarely sufficient and usually cannot provide an unambiguous test for the theory. As a result, the conventional way of constructing a multi-prey functional response is speculative, and often based on assumptions that are difficult to verify. Predator responses allowing for prey selectivity and active switching are thought to be more biologically relevant compared to the standard proportion-based consumption. However, here we argue that the functional responses with switching may not be applicable to communities with a broad spectrum of resource types. We formulate a set of general rules that a biologically sound parameterization of a predator functional response should satisfy, and show that all existing formulations for the multispecies response with prey selectivity and switching fail to do so. Finally, we propose a universal framework for parameterization of a multi-prey functional response by combining patterns of food selectivity and proportion-based feeding.