Bifurcation structure of a chemostat model for an age-structured predator and its prey

We model a chemostat containing an age-structured predator and its prey using a linear function for the uptake of substrate by the prey and two different functional responses (linear and Monod) for the consumption of prey by the predator. Limit cycles (LCs) caused by the predator's age structure arise at Hopf bifurcations at low values of the chemostat dilution rate for both model cases. In addition, LCs caused by the predator–prey interaction arise for the case with the Monod functional response. At low dilution rates in the Monod case, the age structure causes cycling at lower values of the inflowing resource concentration and conversely prevents cycling at higher values of the inflowing resource concentration. The results shed light on a similar model by Fussmann et al. [G. Fussmann, S. Ellner, K. Shertzer, and N. Hairston, Crossing the Hopf bifurcation in a live predator–prey system, Science 290 (2000), pp. 1358–1360.], which correctly predicted conditions for the onset of cycling in a chemostat containing an age-structured rotifer population feeding on algal prey.     

Limit cycles in a chemostat model for a single species with age structure

We model an age-structured population feeding on an abiotic resource by combining the Gurtin-MacCamy [Math. Biosci. 43 (1979) 199] approach with a standard chemostat model. Limit cycles arise by Hopf bifurcations at low values of the chemostat dilution rate, even for simple maternity functions for which the original Gurtin-MacCamy model has no oscillatory solutions. We find the exact location in parameter space of the Hopf bifurcations for special cases of our model. The onset of cycling is largely independent of both the form of the resource uptake function and the shape of the maternity function.     

Effects of rapid prey evolution on predator–prey cycles

We study the qualitative properties of population cycles in a predator-prey system where genetic variability allows contemporary rapid evolution of the prey. Previous numerical studies have found that prey evolution in response to changing predation risk can have major quantitative and qualitative effects on predator-prey cycles, including: (1) large increases in cycle period, (2) changes in phase relations (so that predator and prey are cycling exactly out of phase, rather than the classical quarter-period phase lag), and (3) "cryptic" cycles in which total prey density remains nearly constant while predator density and prey traits cycle. Here we focus on a chemostat model motivated by our experimental system (Fussmann et al. in Science 290:1358-1360, 2000; Yoshida et al. in Proc roy Soc Lond B 424:303-306, 2003) with algae (prey) and rotifers (predators), in which the prey exhibit rapid evolution in their level of defense against predation. We show that the effects of rapid prey evolution are robust and general, and furthermore that they occur in a specific but biologically relevant region of parameter space: when traits that greatly reduce predation risk are relatively cheap (in terms of reductions in other fitness components), when there is coexistence between the two prey types and the predator, and when the interaction between predators and undefended prey alone would produce cycles. Because defense has been shown to be inexpensive, even cost-free, in a number of systems (Andersson et al. in Curr Opin Microbiol 2:489-493, 1999: Gagneux et al. in Science 312:1944-1946, 2006; Yoshida et al. in Proc Roy Soc Lond B 271:1947-1953, 2004), our discoveries may well be reproduced in other model systems, and in nature. Finally, some of our key results are extended to a general model in which functional forms for the predation rate and prey birth rate are not specified.     

A comparison of two predator-prey models with Holling's type I functional response

In this paper, we analyze a laissez-faire predator-prey model and a Leslie-type predator-prey model with type I functional responses. We study the stability of the equilibrium where the predator and prey coexist by both performing a linearized stability analysis and by constructing a Lyapunov function. For the Leslie-type model, we use a generalized Jacobian to determine how eigenvalues jump at the corner of the functional response. We show, numerically, that our two models can both possess two limit cycles that surround a stable equilibrium and that these cycles arise through global cyclic-fold bifurcations. The Leslie-type model may also exhibit super-critical and discontinuous Hopf bifurcations. We then present and analyze a new functional response, built around the arctangent, that smoothes the sharp corner in a type I functional response. For this new functional response, both models undergo Hopf, cyclic-fold, and Bautin bifurcations. We use our analyses to characterize predator-prey systems that may exhibit bistability.     

Individual base model of predator-prey system shows predator dominant dynamics

Individual base model of predator-prey system is constructed. Both predator and prey species have age structure and cohorts of early reproductive age have competitive advantage. The model has linear functional response in predation behavior and includes the effect of interference among predators and delay of population growth from resource intake, not by functional response but by calculation procedure. Each foraging action is calculated successively and surplus or scarce of acquired resources is interpreted into population size through individual birth and death. This model shows that biomass of prey killed by predator is dependent on demand of predator and that heterogeneity in predator population is essential in persistency and stability of predator-prey system. Heterogeneity of predator makes predator individuals of less competing ability die rapidly. Rapid death of weak individuals causes rapid decrease of total demand of predator and that makes enough room for survived predators. Therefore, the biomass of killed prey is dependent on predator's demand. As young or infant population of predator are the more vulnerable to shortage of prey, and when many of them cannot survive to reproductive age, they can stabilize the system by wasting excessive prey with only temporal numerical increase of predator population.     

Complex dynamics in a model microbial system

The forced double-Monod model (for a chemostat with a predator, a prey and periodically forced inflowing substrate) displays quasiperiodicity, phase locking, period doubling and chaotic dynamics. Stroboscopic sections reveal circle maps for the quasiperiodic regimes and noninvertible maps of the interval for the chaotic regimes. Criticality in the circle maps sets the stage for chaos in the model. This criticality may arise with an increase in the period or amplitude of forcing.     

Alternative stable states in communities with intraguild predation

Intraguild predation (IGP) is a widespread ecological phenomenon in which two consumers that share a resource also engage in a predator-prey interaction. Theory on IGP predicts the occurrence of alternative stable states, but empirical evidence of such states is scarce. This raises the question of whether alternative states are a rare phenomenon that is unlikely to be observed in nature. Here we analyze a model in which the resource exhibits logistic or chemostat dynamics and consumers have saturating (Type II) functional responses. We show that alternative states can arise for a wide range of biological scenarios and that environmental constraints can make their detection difficult. Our analysis identifies three possible combinations of alternative states: (i) IG prey or IG predator, (ii) coexistence or IG predator, and (iii) coexistence or IG prey. Bifurcation diagrams reveal that alternative states are possible over large regions of the parameter space. However, they can be limited to narrow ranges along the resource productivity axis, which may make it difficult to observe the occurrence of alternative states in communities with IGP. Microcosm experiments provide a promising avenue for detecting combinations of asymptotically stable states along a productivity gradient.     

Endogenous metabolism and the stability of microbial prey-predator systems

We formulate a variant of the "double Monod model" which takes explicit account of endogenous metabolism. Using parameter values appropriate to carbohydrate-limited substrate, bacterial prey, and protozoan predator, we study the stability of steady states under chemostat conditions. We conclude that the predator's endogenous metabolism may have a stabilizing effect at low dilution rates.     

Effect of predator density dependent dispersal of prey on stability of a predator-prey system

This work presents a predator-prey Lotka-Volterra model in a two patch environment. The model is a set of four ordinary differential equations that govern the prey and predator population densities on each patch. Predators disperse with constant migration rates, while prey dispersal is predator density-dependent. When the predator density is large, the dispersal of prey is more likely to occur. We assume that prey and predator dispersal is faster than the local predator-prey interaction on each patch. Thus, we take advantage of two time scales in order to reduce the complete model to a system of two equations governing the total prey and predator densities. The stability analysis of the aggregated model shows that a unique strictly positive equilibrium exists. This equilibrium may be stable or unstable. A Hopf bifurcation may occur, leading the equilibrium to be a centre. If the two patches are similar, the predator density dependent dispersal of prey has a stabilizing effect on the predator-prey system.     

Egg-eating age-structured predators in interaction with age-structured prey

In this paper, we consider a system of integrodifferential equations which models a predator-prey system with both species (predator and prey) age-structured and predators living only on the eggs of prey. The present model is a generalization of the model given in [20]. The existence, stability, and instability of nonnegative equilibria is studied assuming a general fecundity rate function for the prey. With a special choice of fecundity rate function for the predator it is shown here that a large maturation period m of the predator leads to stability. This seems to be contrary to the usual rule of thumb that increasing delays in growth rate responses cause instabilities.     

The dynamics of two diffusively coupled predator-prey populations

I analyze the dynamics of predator and prey populations living in two patches. Within a patch the prey grow logistically and the predators have a Holling type II functional response. The two patches are coupled through predator migration. The system can be interpreted as a simple predator-prey metapopulation or as a spatially explicit predator-prey system. Asynchronous local dynamics are presumed by metapopulation theory. The main question I address is when synchronous and when asynchronous dynamics arise. Contrary to biological intuition, for very small migration rates the oscillations always synchronize. For intermediate migration rates the synchronous oscillations are unstable and I found periodic, quasi-periodic, and intermittently chaotic attractors with asynchronous dynamics. For large predator migration rates, attractors in the form of equilibria or limit cycles exist in which one of the patches contains no prey. The dynamical behavior of the system is described using bifurcation diagrams. The model shows that spatial predator-prey populations can be regulated through the interplay of local dynamics and migration.     

Effects of Density Dependent Migrations on the Dynamics of a Predator Prey Model

We study the effects of density dependent migrations on the stability of a predator-prey model in a patchy environment which is composed with two sites connected by migration. The two patches are different. On the first patch, preys can find resource but can be captured by predators. The second patch is a refuge for the prey and thus predators do not have access to this patch. We assume a repulsive effect of predator on prey on the resource patch. Therefore, when the predator density is large on that patch, preys are more likely to leave it to return to the refuge. We consider two models. In the first model, preys leave the refuge to go to the resource patch at constant migration rates. In the second model, preys are assumed to be in competition for the resource and leave the refuge to the resource patch according to the prey density. We assume two different time scales, a fast time scale for migration and a slow time scale for population growth, mortality and predation. We take advantage of the two time scales to apply aggregation of variables methods and to obtain a reduced model governing the total prey and predator densities. In the case of the first model, we show that the repulsive effect of predator on prey has a stabilizing effect on the predator-prey community. In the case of the second model, we show that there exists a window for the prey proportion on the resource patch to ensure stability.     

Prey-taxis destabilizes homogeneous stationary state in spatial Gause–Kolmogorov-type model for predator–prey system

We consider a continuous taxis-diffusion-reaction system of partial-differential equations describing spatiotemporal dynamics of a predator–prey system. The local kinetics of the system is defined by general Gause–Kolmogorov-type model. The predator ability to pursue the prey is modelled by the Patlak–Keller–Segel taxis model, assuming that movement velocities of predators are proportional to the gradients of specific cues emitted by prey (e.g., odour, pheromones, exometabolites). The linear stability analysis of the model showed that the non-trivial homogeneous stationary regime of the model becomes unstable with respect to small heterogeneous perturbations with increase of prey-taxis activity; an Andronov–Hopf bifurcation occurs in the system when the taxis coefficient of predator exceeds its critical bifurcation value that exists for all admissible values of model parameters. These findings generalize earlier results obtained for particular cases of the Gause–Kolmogorov-type model assuming logistic reproduction of the prey population and the Holling types I and II functional responses of the predator population. Numerical simulations with theta-logistic growth of the prey population and the Ivlev functional response of predators illustrate and support results of the analytical study.     

PREY ADAPTATION AS A CAUSE OF PREDATOR-PREY CYCLES

We analyze simple models of predator-prey systems in which there is adaptive change in a trait of the prey that determines the rate at which it is captured by searching predators. Two models of adaptive change are explored: (1) change within a single reproducing prey population that has genetic variation for vulnerability to capture by the predator; and (2) direct competition between two independently reproducing prey populations that differ in their vulnerability. When an individual predator's consumption increases at a decreasing rate with prey availability, prey adaptation via either of these mechanisms may produce sustained cycles in both species' population densities and in the prey's mean trait value. Sufficiently rapid adaptive change (e.g., behavioral adaptation or evolution of traits with a large additive genetic variance), or sufficiently low predator birth and death rates will produce sustained cycles or chaos, even when the predator-prey dynamics with fixed prey capture rates would have been stable. Adaptive dynamics can also stabilize a system that would exhibit limit cycles if traits were fixed at their equilibrium values. When evolution fails to stabilize inherently unstable population interactions, selection decreases the prey's escape ability, which further destabilizes population dynamics. When the predator has a linear functional response, evolution of prey vulnerability always promotes stability. The relevance of these results to observed predator-prey cycles is discussed.     

Impact of fear effect on the growth of prey in a predator-prey interaction model

Several field data and experiments on a terrestrial vertebrates exhibited that the fear of predators would cause a substantial variability of prey demography. Fear for predator population enhances the survival probability of prey population, and it can greatly reduce the reproduction of prey population. Based on the experimental evidence, we proposed and analyzed a prey-predator system introducing the cost of fear into prey reproduction with Holling type-II functional response. We investigate all the biologically feasible equilibrium points, and their stability is analyzed in terms of the model parameters. Our mathematical analysis exhibits that for strong anti-predator responses can stabilize the prey-predator interactions by ignoring the existence of periodic behaviors. Our model system undergoes Hopf bifurcation by considering the birth rate r₀ as a bifurcation parameter. For larger prey birth rate, we investigate the transition to a stable coexisting equilibrium state, with oscillatory approach to this equilibrium state, indicating that the greatest characteristic eigenvalues are actually a pair of imaginary eigenvalues with real part negative, which is increasing for r₀. We obtained the conditions for the occurrence of Hopf bifurcation and conditions governing the direction of Hopf bifurcation, which imply that the prey birth rate will not only influence the occurrence of Hopf bifurcation but also alter the direction of Hopf bifurcation. We identify the parameter regions associated with the extinct equilibria, predator-free equilibria and coexisting equilibria with respect to prey birth rate, predator mortality rates. Fear can stabilize the predator-prey system at an interior steady state, where all the species can exists together, or it can create the oscillatory coexistence of all the populations. We performed some numerical simulations to investigate the relationship between the effects of fear and other biologically related parameters (including growth/decay rate of prey/predator), which exhibit the impact that fear can have in prey-predator system. Our numerical illustrations also demonstrate that the prey become less sensitive to perceive the risk of predation with increasing prey growth rate or increasing predators decay rate.     

Enrichment and ecosystem stability: effect of toxic food

Enrichment in resource availability theoretically destabilizes predator-prey dynamics (the paradox of enrichment). However, a minor change in the resource stoichiometry may make a prey toxic for the predator, and the presence of toxic prey affects the dynamics significantly. Here, theoretically we explore how, at increased carrying capacity, a toxic prey affects the oscillation or destabilization of predator-prey dynamics, and how its presence influences the growth of the predator as well as that of a palatable prey. Mathematical analysis determines the bounds on the food toxicity that allow the coexistence of a predator along with a palatable and a toxic prey. The overall results demonstrate that toxic food counteracts oscillation (destabilization) arising from enrichment of resource availability. Moreover, our results show that, at increased resource availability, toxic food that acts as a source of extra mortality may increase the abundance of the predator as well as that of the palatable prey.     

Impacts of Biotic Resource Enrichment on a Predator–Prey Population

The environmental carrying capacity is usually assumed to be fixed quantity in the classical predator–prey population growth models. However, this assumption is not realistic as the environment generally varies with time. In a bid for greater realism, functional forms of carrying capacities have been widely applied to describe varying environments. Modelling carrying capacity as a state variable serves as another approach to capture the dynamical behavior between population and its environment. The proposed modified predator–prey model is based on the ratio-dependent models that have been utilized in the study of food chains. Using a simple non-linear system, the proposed model can be linked to an intra-guild predation model in which predator and prey share the same resource. Distinct from other models, we formulate the carrying capacity proportional to a biotic resource and both predator and prey species can directly alter the amount of resource available by interacting with it. Bifurcation and numerical analyses are presented to illustrate the system’s dynamical behavior. Taking the enrichment parameter of the resource as the bifurcation parameter, a Hopf bifurcation is found for some parameter ranges, which generate solutions that posses limit cycle behavior.     

Testing for predator dependence in predator-prey dynamics: a non-parametric approach

The functional response is a key element in all predator-prey interactions. Although functional responses are traditionally modelled as being a function of prey density only, evidence is accumulating that predator density also has an important effect. However, much of the evidence comes from artificial experimental arenas under conditions not necessarily representative of the natural system, and neglecting the temporal dynamics of the organism (in particular the effects of prey depletion on the estimated functional response). Here we present a method that removes these limitations by reconstructing the functional response non-parametrically from predator-prey time-series data. This method is applied to data on a protozoan predator-prey interaction, and we obtain significant evidence of predator dependence in the functional response. A crucial element in this analysis is to include time-lags in the prey and predator reproduction rates, and we show that these delays improve the fit of the model significantly. Finally, we compare the non-parametrically reconstructed functional response to parametric forms, and suggest that a modified version of the Hassell-Varley predator interference model provides a simple and flexible function for theoretical investigation and applied modelling.     

Modeling and analysis of a predator-prey model with disease in the prey

A system of retarded functional differential equations is proposed as a predator-prey model with disease in the prey. Mathematical analyses of the model equations with regard to invariance of non-negativity, boundedness of solutions, nature of equilibria, permanence and global stability are analyzed. If the coefficient in conversing prey into predator k=k(0) is constant (independent of delay tau;, gestation period), we show that positive equilibrium is locally asymptotically stable when time delay tau; is suitable small, while a loss of stability by a Hopf bifurcation can occur as the delay increases. If k=k(0)e(-dtau;) (d is the death rate of predator), numerical simulation suggests that time delay has both destabilizing and stabilizing effects, that is, positive equilibrium, if it exists, will become stable again for large time delay. A concluding discussion is then presented.