Chaotic dynamics of a food web in a chemostat

We analyze a mathematical model of a simple food web consisting of one predator and two prey populations in a chemostat. Monod's model is employed for the dependence of the specific growth rates of the two prey populations on the concentration of the rate-limiting substrate and a generalization of Monod's model for the dependence of the specific growth rate of the predator on the concentrations of the prey populations. We use numerical bifurcation techniques to determine the effect of the operating conditions of the chemostat on the dynamics of the system and construct its operating diagram. Chaotic behavior resulting from successive period doublings is observed. Multistability phenomena of coexistence of steady and periodic states at the same operating conditions are also found.     

Revealing the role of predator interference in a predator–prey system with disease in prey population

Predation on a species subjected to an infectious disease can affect both the infection level and the population dynamics. There is an ongoing debate about the act of managing disease in natural populations through predation. Recent theoretical and empirical evidence shows that predation on infected populations can have both positive and negative influences on disease in prey populations. Here, we present a predator–prey system where the prey population is subjected to an infectious disease to explore the impact of predator on disease dynamics. Specifically, we investigate how the interference among predators affects the dynamics and structure of the predator–prey community. We perform a detailed numerical bifurcation analysis and find an unusually large variety of complex dynamics, such as, bistability, torus and chaos, in the presence of predators. We show that, depending on the strength of interference among predators, predators enhance or control disease outbreaks and population persistence. Moreover, the presence of multistable regimes makes the system very sensitive to perturbations and facilitates a number of regime shifts. Since, the habitat structure and the choice of predators deeply influence the interference among predators, thus before applying predators to control disease in prey populations or applying predator control strategy for wildlife management, it is essential to carefully investigate how these predators interact with each other in that specific habitat; otherwise it may lead to ecological disaster.     

Decomposing Predation: Testing for Parameters that Correlate with Predatory Performance by a Social Bacterium

Predator–prey interactions presumably play major roles in shaping the composition and dynamics of microbial communities. However, little is understood about the population biology of such interactions or how predation-related parameters vary or correlate across prey environments. Myxococcus xanthus is a motile soil bacterium that feeds on a broad range of other soil microbes that vary greatly in the degree to which they support M. xanthus growth. In order to decompose predator–prey interactions at the population level, we quantified five predation-related parameters during M. xanthus growth on nine phylogenetically diverse bacterial prey species. The horizontal expansion rate of swarming predator colonies fueled by prey lawns served as our measure of overall predatory performance, as it incorporates both the searching (motility) and handling (killing and consumption of prey) components of predation. Four other parameters—predator population growth rate, maximum predator yield, maximum prey kill, and overall rate of prey death—were measured from homogeneously mixed predator–prey lawns from which predator populations were not allowed to expand horizontally by swarming motility. All prey species fueled predator population growth. For some prey, predator-specific prey death was detected contemporaneously with predator population growth, whereas killing of other prey species was detected only after cessation of predator growth. All four of the alternative parameters were found to correlate significantly with predator swarm expansion rate to varying degrees, suggesting causal interrelationships among these diverse predation measures. More broadly, our results highlight the importance of examining multiple parameters for thoroughly understanding the population biology of microbial predation.     

Effects of a disease affecting a predator on the dynamics of a predator-prey system

We study the effects of a disease affecting a predator on the dynamics of a predator-prey system. We couple an SIRS model applied to the predator population, to a Lotka-Volterra model. The SIRS model describes the spread of the disease in a predator population subdivided into susceptible, infected and removed individuals. The Lotka-Volterra model describes the predator-prey interactions. We consider two time scales, a fast one for the disease and a comparatively slow one for predator-prey interactions and for predator mortality. We use the classical “aggregation method” in order to obtain a reduced equivalent model. We show that there are two possible asymptotic behaviors: either the predator population dies out and the prey tends to its carrying capacity, or the predator and prey coexist. In this latter case, the predator population tends either to a “disease-free” or to a “disease-endemic” state. Moreover, the total predator density in the disease-endemic state is greater than the predator density in the “disease-free” equilibrium (DFE).     

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.     

Mechanisms of coexistence in diverse herbivore–carnivore assemblages: demographic,temporal and spatial heterogeneities affecting prey vulnerability

Simple models coupling the dynamics of single predators to single prey populations tend to generate oscillatory dynamics of both predator and prey, or extirpation of the prey followed by that of the predator. In reality, such oscillatory dynamics may be counteracted by prey refugia or by opportunities for prey switching by the predator in multi‐prey assemblages. How these mechanisms operate depends on relative prey vulnerability, a factor ignored in simple interactive models. I outline how compositional, temporal, demographic and spatial heterogeneities help explain the contrasting effects of top predators on large herbivore abundance and population dynamics in species‐rich African savanna ecosystems compared with less species‐diverse northern temperate or subarctic ecosystems. Demographically, mortality inflicted by predation depends on the relative size and life history stage of the prey. Because all animals eventually die and are consumed by various carnivores, the additive component of the mortality inflicted is somewhat less than the predation rate. Prey vulnerability varies annually and seasonally, and between day and night. Spatial variation in the risk of predation depends on vegetation cover as well as on the availability of food resources. During times of food shortage, herbivores become prompted to occupy more risky habitats retaining more food. Predator concentrations dependent on the abundance of primary prey species may restrict the occurrence of other potential prey species less resistant to predation. The presence of multiple herbivore species of similar size in African savannas allows the top predator, the lion, to shift its prey selection flexibly dependent on changing prey vulnerability. Hence top–down and bottom–up influences on herbivore populations are intrinsically entangled. Models coupling the population dynamics of predators and prey need to accommodate the changing influences of prey demography, temporal variation in environmental conditions, and spatial variation in the relative vulnerability of alternative prey species to predation. Synthesis While re‐established predators have had major impacts on prey populations in northern temperate regions, multiple large herbivore species typically coexist along with diverse carnivores in African savanna ecosystems. In order to explain these contrasting outcomes, certain functional heterogeneities must be recognised, including relative vulnerability of alternative prey, temporal variation in the risk of predation, demographic differences in susceptibility to predation, and spatial contrasts in exposure to predation. Food shortfalls prompt herbivores to exploit more risky habitats, meaning that top–down and bottom–up influences on prey populations are intrinsically entangled. Models coupling the interactive dynamics of predator and prey populations need to incorporate these varying influences on relative prey vulnerability.     

食饵庇护所对斑块环境下Leslie-Gower捕食系统的影响

建立了一个显式含有空间庇护所的两斑块Leslie-Gower捕食者-食饵系统。假设只有食饵种群在斑块间以常数迁移率迁移,且在每个斑块上食饵间的迁移比局部捕食者-食饵相互作用发生的时间尺度要快。利用两个时间尺度,可以构建用来描述所有斑块总的食饵和捕食者密度的综合系统。数学分析表明,在一定条件下,存在唯一的正平衡点,并且此平衡点全局稳定。进一步,捕食者的数量随着食饵庇护所数量增加而降低;在一定条件下,食饵的数量随着食饵庇护所数量增加先增加后降低,在足够强的庇护所强度下,两物种出现灭绝。对比以往研究,利用显式含有和隐含空间庇护所的数学模型所得结论不一致,这意味着在研究庇护所对捕食系统种群动态影响时,空间结构可能起着重要作用。     

The Lotka-Volterra predator-prey model with foraging-predation risk trade-offs

This article studies the effects of adaptive changes in predator and/or prey activities on the Lotka-Volterra predator-prey population dynamics. The model assumes the classical foraging-predation risk trade-offs: increased activity increases population growth rate, but it also increases mortality rate. The model considers three scenarios: prey only are adaptive, predators only are adaptive, and both species are adaptive. Under all these scenarios, the neutral stability of the classical Lotka-Volterra model is partially lost because the amplitude of maximum oscillation in species numbers is bounded, and the bound is independent of the initial population numbers. Moreover, if both prey and predators behave adaptively, the neutral stability can be completely lost, and a globally stable equilibrium would appear. This is because prey and/or predator switching leads to a piecewise constant prey (predator) isocline with a vertical (horizontal) part that limits the amplitude of oscillations in prey and predator numbers, exactly as suggested by Rosenzweig and MacArthur in their seminal work on graphical stability analysis of predator-prey systems. Prey and predator activities in a long-term run are calculated explicitly. This article shows that predictions based on short-term behavioral experiments may not correspond to long-term predictions when population dynamics are considered.     

Consequences of symbiosis for food web dynamics

Basic Lotka-Volterra type models in which mutualism (a type of symbiosis where the two populations benefit both) is taken into account, may give unbounded solutions. We exclude such behaviour using explicit mass balances and study the consequences of symbiosis for the long-term dynamic behaviour of a three species system, two prey and one predator species in the chemostat. We compose a theoretical food web where a predator feeds on two prey species that have a symbiotic relationships. In addition to a species-specific resource, the two prey populations consume the products of the partner population as well. In turn, a common predator forages on these prey populations. The temporal change in the biomass and the nutrient densities in the reactor is described by ordinary differential equations (ODE). Since products are recycled, the dynamics of these abiotic materials must be taken into account as well, and they are described by odes in a similar way as the abiotic nutrients. We use numerical bifurcation analysis to assess the long-term dynamic behaviour for varying degrees of symbiosis. Attractors can be equilibria, limit cycles and chaotic attractors depending on the control parameters of the chemostat reactor. These control parameters that can be experimentally manipulated are the nutrient density of the inflow medium and the dilution rate. Bifurcation diagrams for the three species web with a facultative symbiotic association between the two prey populations, are similar to that of a bi-trophic food chain; nutrient enrichment leads to oscillatory behaviour. Predation combined with obligatory symbiotic prey-interactions has a stabilizing effect, that is, there is stable coexistence in a larger part of the parameter space than for a bi-trophic food chain. However, combined with a large growth rate of the predator, the food web can persist only in a relatively small region of the parameter space. Then, two zero-pair bifurcation points are the organizing centers. In each of these points, in addition to a tangent, transcritical and Hopf bifurcation a global heteroclinic bifurcation is emanating. This heteroclinic cycle connects two saddle equilibria where the predator is absent. Under parameter variation the period of the stable limit cycle goes to infinity and the cycle tends to the heteroclinic cycle. At this global bifurcation point this cycle breaks and the boundary of the basin of attraction disappears abruptly because the separatrix disappears together with the cycle. As a result, it becomes possible that a stable two-nutrient–two-prey population system becomes unstable by invasion of a predator and eventually the predator goes extinct together with the two prey populations, that is, the complete food web is destroyed. This is a form of over-exploitation by the predator population of the two symbiotic prey populations. When obligatory symbiotic prey-interactions are modelled with Liebigs minimum law, where growth is limited by the most limiting resource, more complicated types of bifurcations are found. This results from the fact that the Jacobian matrix changes discontinuously with respect to a varying parameter when another resource becomes most limiting.Revised version: 21 July 2003     

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.     

High reproductive rates result in high predation risks: a mechanism promoting the coexistence of competing prey in spatially structured populations

I tested the hypothesis that spatial structure provides a trade-off between reproduction and predation risk and thereby facilitates predator-mediated coexistence of competing prey species. I compared a cellular automata model to a mean-field model of two prey species and their common predator. In the mean-field model, the prey species with the higher reproductive rate (the superior competitor) always outcompeted the other species (the inferior competitor), both in the presence of and the absence of the predator. In the cellular automata model, both prey species, which differed only in their reproductive rates, coexisted for a long time in the presence of their common predator at intermediate levels of predation. At low predation rates, the superior competitor dominated, while high predation rates favored the inferior competitor. This discrepancy in the results of the different models was due to a trade-off that spontaneously emerged in spatially structured populations; that is, the more clustered distribution of the superior competitor made it more susceptible to predation. In addition, coexistence of competing prey species declined with increasing dispersal ranges of either prey or predator, which suggests that the trade-off that results from spatial structure becomes less important as either prey or predator disperse over a broader range.     

The role of prey taxis in biological control: a spatial theoretical model

We study a reaction-diffusion-advection model for the dynamics of populations under biological control. A control agent is assumed to be a predator species that has the ability to perceive the heterogeneity of pest distribution. The advection term represents the predator density movement according to a basic prey taxis assumption: acceleration of predators is proportional to the prey density gradient. The prey population reproduces logistically, and the local population interactions follow the Holling Type II trophic function. On the scale of the population, our spatially explicit approach subdivides the predation process into random movement represented by diffusion, directed movement described by prey taxis, local prey encounters, and consumption modeled by the trophic function. Thus, our model allows studying the effects of large-scale predator spatial activity on population dynamics. We show under which conditions spatial patterns are generated by prey taxis and how this affects the predator ability to maintain the pest population below some economic threshold. In particular, intermediate taxis activity can stabilize predator-pest populations at a very low level of pest density, ensuring successful biological control. However, very intensive prey taxis destroys the stability, leading to chaotic dynamics with pronounced outbreaks of pest density.     

An interfacial approach to regional segregation of two competing species mediated by a predator

We consider the problem of coexistence of two competing species mediated by the presence of a predator. We employ a reaction-diffusion model equation with Lotka-Volterra interaction, and speculate that the possibility of coexistence is,enhanced by differences in the diffusion rates of the prey and their predator. In the limit where the diffusion rate of the prey tends to zero, a new equation is derived and the dynamics of spatial segregation is discussed by means of the interfacial dynamics approach. Also, we show that spatial segregation permits periodic and chaotic dynamics for certain parameter ranges.     

Antagonistic species interaction drives selection for sex in a predator–prey system

The evolutionary maintenance of sexual reproduction has long challenged biologists as the majority of species reproduce sexually despite inherent costs. Providing a general explanation for the evolutionary success of sex has thus proven difficult and resulted in numerous hypotheses. A leading hypothesis suggests that antagonistic species interaction can generate conditions selecting for increased sex due to the production of rare or novel genotypes that are beneficial for rapid adaptation to recurrent environmental change brought on by antagonism. To test this ecology‐based hypothesis, we conducted experimental evolution in a predator (rotifer)–prey (algal) system by using continuous cultures to track predator–prey dynamics and in situ rates of sex in the prey over time and within replicated experimental populations. Overall, we found that predator‐mediated fluctuating selection for competitive versus defended prey resulted in higher rates of genetic mixing in the prey. More specifically, our results showed that fluctuating population sizes of predator and prey, coupled with a trade‐off in the prey, drove the sort of recurrent environmental change that could provide a benefit to sex in the prey, despite inherent costs. We end with a discussion of potential population genetic mechanisms underlying increased selection for sex in this system, based on our application of a general theoretical framework for measuring the effects of sex over time, and interpreting how these effects can lead to inferences about the conditions selecting for or against sexual reproduction in a system with antagonistic species interaction.     

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.     

Pulsed-resource dynamics constrain the evolution of predator-prey interactions

Although temporal variability in the physical environment plays a major role in population fluctuations, little is known about how it drives the ecological and evolutionary dynamics of species interactions. We studied experimentally how extrinsic resource pulses affect evolutionary and ecological dynamics between the prey bacterium Serratia marcescens and the predatory protozoan Tetrahymena thermophila. Predation increased the frequency of defensive, nonpigmented prey types, which bore competitive costs in terms of reduced maximum growth rate, most in a constant-resource environment. Furthermore, the predator densities of the pulsed-resource environment regularly fluctuated above and below the mean predator densities of the constant environment. These results suggest that selection favored fast-growing competitor prey types over defensive but slower-growing prey types more often in the pulsed-resource environment (abundance of resources and low predation risk). As a result, the selection for prey defense fluctuated more in the pulsed-resource environment, leading to a weaker mean response in prey defense. At the ecological level, the evolution of prey defense weakened the relative strength of top-down regulation on prey community. This was more evident in the constant-resource environment, whereas the slow emergence of defensive prey types gradually decreased the amplitude of predator peaks in the pulsed-resource environment. Our study suggests that rapid evolution plays a smaller role in the ecological dynamics of communities dominated by resource pulses.     

Mathematical models for continuous culture growth dynamics of mixed populations subsisting on a heterogeneous resource base. II. Predation and trophic structure

Models are presented for the joint dynamics of predators and prey, maintained in continuous flow chemostat culture. The predators are visualized as subsisting on one or more prey organisms, which in turn are visualized as subsisting on one or more substrate resources supplied by the investigator. The dynamic equations are translated into an analogous Lotka-Volterra predation model, and the criteria for the existence and stability of various equilibria are indicated. Denoting the number of different predator organisms as N_H, the number of different prey organisms by N_I and the number of different substrates as N_J, it is shown that the joint coexistence of all components requires 0 ? N_I ? N_H ? N_J. The model is extended to more complex situations by including additional trophic layers and by allowing trophic layer “leap-frogging.” The model may always be translated into an approximately quadratic differential equation of the Lotka-Volterra type. The α- and β-coefficients of these latter are really variables, and become quite complex for some of the multi-layered models.     

A comparison of foraging strategies in a patchy environment.

In this paper we compare foraging strategies that might be used by predators seeking prey in a patchy environment. The strategies differ in the extent to which predators aggregate in response to prey density. The approach to the comparison is suggested by the idea of evolutionarily stable strategies. A strategy is said to be evolutionarily stable if it cannot be invaded by another strategy. Thus we examine scenarios where a small number of individuals using one strategy are introduced into a situation where a large number of individuals using the other strategy are already present. However, our foraging models do not explicitly incorporate predator population dynamics, so we use net energy uptake as a surrogate for reproductive fitness. In cases where all of the patches visited by predators sustain prey populations, we find that for any pair of strategies one of them will have a higher net energy uptake than the other whether it is the resident or the introduced strain. However, which one is higher will typically depend on the total predator population, which is determined by the resident strain. If the predators leave prey densities high, the more aggregative strain will have the advantage. If the predators reduce prey densities to low levels the less aggregative strain will have the advantage. In cases where one strain of predators aggregates in response to prey density and the other does not, then there might be patches which do not contain prey but do contain (non-aggregating) predators. In those cases, there is the possibility that whichever strategy is used by the introduced strain will yield a higher energy uptake than that used by the resident strain. This suggests that if some patches are empty of prey then aggregative and non-aggregative strategies may be able to coexist.     

Predation across spatial scales in heterogeneous environments

A hierarchy of scales is introduced to the spatially heterogeneous Lotka-Volterra predator-prey diffusion model, and its effects on the model's spatial and temporal behavior are studied. When predators move on a large scale relative to prey, local coupling of the predator-prey interaction is replaced by global coupling. Prey with low dispersal ability become narrowly confined to the most productive habitats, strongly amplifying the underlying spatial pattern of the environment. As prey diffusion rate increases, the prey distribution spreads out and predator abundance declines. The model retains neutrally stable Lotka-Volterra temporal dynamics: different scales of predator and prey dispersal do not stabilize the interaction. The model predicts that, for prey populations that are limited by widely ranging predators, species with low dispersal ability should be restricted to discrete high density patches, and those with greater mobility should be more uniformly distributed at lower density.