Multiple attractors and boundary crises in a tri-trophic food chain

Crisis, unstable dimension variability, and bifurcations in a system with high-dimensional phase space: coupled sine circle maps

The phenomenon of crisis in systems evolving in high-dimensional phase space can show unexpected and interesting features. We study this phenomenon in the context of a system of coupled sine circle maps. We establish that the origins of this crisis lie in a tangent bifurcation in high dimensions, and identify the routes that lead to the crisis. Interestingly, multiple routes to crisis are seen depending on the initial conditions of the system, due to the high dimensionality of the space in which the system evolves. The statistical behavior seen in the phase diagram of the system is also seen to change due to the dynamical phenomenon of crisis, which leads to transitions from nonspreading to spreading behavior across an infection line in the phase diagram. Unstable dimension variability is seen in the neighborhood of the infection line. We characterize this crisis and unstable dimension variability using dynamical characterizers, such as finite-time Lyapunov exponents and their distributions. The phase diagram also contains regimes of spatiotemporal intermittency and spatial intermittency, where the statistical quantities scale as power laws. We discuss the signatures of these regimes in the dynamic characterizers, and correlate them with the statistical characterizers and bifurcation behavior. We find that it is necessary to look at both types of correlators together to build up an accurate picture of the behavior of the system.

Alternative community compositional and dynamical states: the dual consequences of assembly history

1. Much work on ecological consequences of community assembly history has focused on the formation of history-induced alternative stable equilibria. We hypothesize that assembly history may affect not only community composition but also population dynamics, with assembled communities differing in species composition potentially residing in different dynamical states. 2. We provided an empirical test of the aforementioned hypothesis using a laboratory microcosm experiment that manipulated both the colonization order of three bacterivorous protist species in the presence of a protist predator and environmental productivity. 3. Both priority effects and random divergence emerged, resulting in two different community compositional states: one characterized by the dominance of one prey species and the other by the extinction of the same prey. While communities in the former state exhibited noncyclic dynamics, the majority of communities in the latter state exhibited cyclic dynamics driven by the interaction between another prey and the predator. 4. Temporal variability of total prey community biovolume consequently differed among communities with different histories. 5. Changing productivity altered priority effects on the structure and dynamics of communities experiencing only certain histories. 6. Our results support the dual (compositional and dynamical) consequences of assembly history and emphasize the importance of incorporating the dynamical view into the field of community assembly.

Ecological consequences of global bifurcations in some food chain models

Food chain models of ordinary differential equations (ode's) are often used in ecology to gain insight in the dynamics of populations of species, and the interactions of these species with each other and their environment. One powerful analysis technique is bifurcation analysis, focusing on the changes in long-term (asymptotic) behaviour under parameter variation. For the detection of local bifurcations there exists standardised software, but until quite recently most software did not include any capabilities for the detection and continuation of global bifurcations. We focus here on the occurrence of global bifurcations in four food chain models, and discuss the implications of their occurrence. In two stoichiometric models (one piecewise continuous, one smooth) there exists a homoclinic bifurcation, that results in the disappearance of a limit cycle attractor. Instead, a stable positive equilibrium becomes the global attractor. The models are also capable of bistability. In two three-dimensional models a Shil'nikov homoclinic bifurcation functions as the organising centre of chaos, while tangencies of homoclinic cycle-to-cycle connections 'cut' the chaotic attractors, which is associated with boundary crises. In one model this leads to extinction of the top predator, while in the other model hysteresis occurs. The types of ecological events occurring because of a global bifurcation will be categorized. Global bifurcations are always catastrophic, leading to the disappearance or merging of attractors. However, there is no 1-on-1 coupling between global bifurcation type and the possible ecological consequences. This only emphasizes the importance of including global bifurcations in the analysis of food chain models.

Efficient estimation of the robustness region of biological models with oscillatory behavior

Robustness is an essential feature of biological systems, and any mathematical model that describes such a system should reflect this feature. Especially, persistence of oscillatory behavior is an important issue. A benchmark model for this phenomenon is the Laub-Loomis model, a nonlinear model for cAMP oscillations in Dictyostelium discoideum. This model captures the most important features of biomolecular networks oscillating at constant frequencies. Nevertheless, the robustness of its oscillatory behavior is not yet fully understood. Given a system that exhibits oscillating behavior for some set of parameters, the central question of robustness is how far the parameters may be changed, such that the qualitative behavior does not change. The determination of such a “robustness region” in parameter space is an intricate task. If the number of parameters is high, it may be also time consuming. In the literature, several methods are proposed that partially tackle this problem. For example, some methods only detect particular bifurcations, or only find a relatively small box-shaped estimate for an irregularly shaped robustness region. Here, we present an approach that is much more general, and is especially designed to be efficient for systems with a large number of parameters. As an illustration, we apply the method first to a well understood low-dimensional system, the Rosenzweig-MacArthur model. This is a predator-prey model featuring satiation of the predator. It has only two parameters and its bifurcation diagram is available in the literature. We find a good agreement with the existing knowledge about this model. When we apply the new method to the high dimensional Laub-Loomis model, we obtain a much larger robustness region than reported earlier in the literature. This clearly demonstrates the power of our method. From the results, we conclude that the biological system underlying is much more robust than was realized until now.

Heteroclinic orbits indicate overexploitation in predator–prey systems with a strong Allee effect

Species establishment in a model system in a homogeneous environment can be dependent not only on the parameter setting, but also on the initial conditions of the system. For instance, predator invasion into an established prey population can fail and lead to system collapse, an event referred to as overexploitation. This phenomenon occurs in models with bistability properties, such as strong Allee effects. The Allee effect then prevents easy re-establishment of the prey species. In this paper, we deal with the bifurcation analyses of two previously published predator-prey models with strong Allee effects. We expand the analyses to include not only local, but also global bifurcations. We show the existence of a point-to-point heteroclinic cycle in these models, and discuss numerical techniques for continuation in parameter space. The continuation of such a cycle in two-parameter space forms the boundary of a region in parameter space where the system collapses after predator invasion, i.e. where overexploitation occurs. We argue that the detection and continuation of global bifurcations in these models are of vital importance for the understanding of the model dynamics.

Parametric dependence in model epidemics. I: contact-related parameters

One of the interesting properties of nonlinear dynamical systems is that arbitrarily small changes in parameter values can induce qualitative changes in behavior. The changes are called bifurcations, and they are typically visualized by plotting asymptotic dynamics against a parameter. In some cases, the resulting bifurcation diagram is unique: irrespective of initial conditions, the same dynamical sequence obtains. In other cases, initial conditions do matter, and there are coexisting sequences. Here we study an epidemiological model in which multiple bifurcation sequences yield to a single sequence in response to varying a second parameter. We call this simplification the emergence of unique parametric dependence (UPD) and discuss how it relates to the model's overall response to parameters. In so doing, we tie together a number of threads that have been developing since the mid-1980s. These include period-doubling; subharmonic resonance, attractor merging and subduction and the evolution of strange invariant sets. The present paper focuses on contact related parameters. A follow-up paper, to be published in this journal, will consider the effects of non-contact related parameters.

Parametric Analysis of a Predator–prey System Stabilized by a Top Predator

We present a complete parametric analysis of a predator-prey system influenced by a top predator. We study ecosystems with abundant nutrient supply for the prey where the prey multiplication can be considered as proportional to its density. The main questions we examine are the following: (1) Can the top predator stabilize such a system at low densities of prey? (2) What possible dynamic behaviors can occur? (3) Under which conditions can the top predation result in the system stabilization? We use a system of two nonlinear ordinary differential equations with the density of the top predator as a parameter. The model is investigated with methods of qualitative theory of ODEs and the theory of bifurcations. The existence of 12 qualitatively different types of dynamics and complex structure of the parametric space are demonstrated. Our studies of phase portraits and parametric diagrams show that a top predator can be an important factor leading to stabilization of the predator-prey system with abundant nutrient supply. Although the model here is applied to the plankton communities with fish (or carnivorous zooplankton) as the top trophic level, the general form of the equations allows applications of our results to other ecological systems.

Numerical bifurcation analysis of a tri-trophic food web with omnivory

We study the consequences of omnivory on the dynamic behaviour of a three species food web under chemostat conditions. The food web consists of a prey consuming a nutrient, a predator consuming a prey and an omnivore which preys on the predator and the prey. For each trophic level an ordinary differential equation describes the biomass density in the reactor. The hyperbolic functional response for single and multi prey species figures in the description of the trophic interactions. There are two limiting cases where the omnivore is a specialist; a food chain where the omnivore does not consume the prey and competition where the omnivore does not prey on the predator. We use bifurcation analysis to study the long-term dynamic behaviour for various degrees of omnivory. Attractors can be equilibria, limit cycles or chaotic behaviour depending on the control parameters of the chemostat. Often multiple attractor occur. In this paper we will discuss community assembly. That is, we analyze how the trophic structure of the food web evolves following invasion where a new invader is introduced one at the time. Generally, with an invasion, the invader settles itself and persists with all other species, however, the invader may also replace another species. We will show that the food web model has a global bifurcation, being a heteroclinic connection from a saddle equilibrium to a limit cycle of saddle type. This global bifurcation separates regions in the bifurcation diagram with different attractors to which the system evolves after invasion. To investigate the consequences of omnivory we will focus on invasion of the omnivore. This simplifies the analysis considerably, for the end-point of the assembly sequence is then unique. A weak interaction of the omnivore with the prey combined with a stronger interaction with the predator seems advantageous.

Quantitative aspects of metabolic organization: a discussion of concepts

Metabolic organization of individual organisms follows simple quantitative rules that can be understood from basic physical chemical principles. Dynamic energy budget (DEB) theory identifies these rules, which quantify how individuals acquire and use energy and nutrients. The theory provides constraints on the metabolic organization of subcellular processes. Together with rules for interaction between individuals, it also provides a basis to understand population and ecosystem dynamics. The theory, therefore, links various levels of biological organization. It applies to all species of organisms and offers explanations for body-size scaling relationships of natural history parameters that are otherwise difficult to understand. A considerable number of popular empirical models turn out to be special cases of the DEB model, or very close numerical approximations. Strong and weak homeostasis and the partitionability of reserve kinetics are cornerstones of the theory and essential for understanding the evolution of metabolic organization.