Rich dynamics in a predator–prey model with both noise and periodic force

A spatial version of the predator–prey model with Holling III functional response, which includes some important factors such as external periodic forces, noise, and diffusion processes is investigated. For the model only with diffusion, it exhibits spiral waves in the two-dimensional space. However, combined with noise, it has the feature of chaotic patterns. Moreover, the oscillations become more obvious when the noise intensity is increased. Furthermore, the spatially extended system with external periodic forces and noise exhibits a resonant pattern and frequency-locking phenomena. These results may help us to understand the effects arising from the undeniable susceptibility to random fluctuations in the real ecosystems.     

Effect of resource subsidies on predator–prey population dynamics: a mathematical model

The influence of a resource subsidy on predator–prey interactions is examined using a mathematical model. The model arises from the study of a biological system involving arctic foxes (predator), lemmings (prey), and seal carcasses (subsidy). In one version of the model, the predator, prey and subsidy all occur in the same location; in a second version, the predator moves between two patches, one containing only the prey and the other containing only the subsidy. Criteria for feasibility and stability of the different equilibrium states are studied both analytically and numerically. At small subsidy input rates, there is a minimum prey carrying capacity needed to support both predator and prey. At intermediate subsidy input rates, the predator and prey can always coexist. At high subsidy input rates, the prey cannot persist even at high carrying capacities. As predator movement increases, the dynamic stability of the predator–prey-subsidy interactions also increases.     

Analysis of a competitive prey–predator system with a prey refuge

Gauss's competitive exclusive principle states that two competing species having analogous environment cannot usually occupy the same space at a time but in order to exploit their common environment in a different manner, they can co-exist only when they are active in different times. On the other hand, several studies on predators in various natural and laboratory situations have shown that competitive coexistence can result from predation in a way by resisting any one prey species from becoming sufficiently abundant to outcompete other species such that the predator makes the coexistence possible. It has also been shown that the use of refuges by a fraction of the prey population exerts a stabilizing effect in the interacting population dynamics. Further, the field surveys in the Sundarban mangrove ecosystem reveal that two detritivorous fishes, viz. Liza parsia and Liza tade (prey population) coexist in nature with the presence of the predator fish population, viz. Lates calcarifer by using refuges.     

Asexual reproduction changes predator population dynamics in a life predator–prey system

Many organisms display oscillations in population size. Theory predicts that these fluctuations can be generated by predator–prey interactions, and empirical studies using life model systems, such as a rotifer-algae community consisting of Brachionus calyciflorus as predator and Chlorella vulgaris as prey, have been successfully used for studying such dynamics. B. calyciflorus is a cyclical parthenogen (CP) and clones often differ in their sexual propensity, that is, the degree to which they engage into sexual or asexual (clonal) reproduction. Since sexual propensities can affect growth rates and population sizes, we hypothesized that this might also affect population oscillations. Here, we studied the dynamical behaviour of B. calyciflorus clones representing either CPs (regularly inducing sex) or obligate parthenogens (OPs). We found that the amplitudes of population cycles to be increased in OPs at low nutrient levels. Several other population dynamic parameters seemed unaffected. This suggests that reproductive mode might be an important additional variable to be considered in future studies of population oscillations.     

Multiple periodic solutions of a delayed predator–prey model with non-monotonic functional response and stage structure

The paper studies a periodic and delayed predator–prey system with non-monotonic functional responses and stage structure. In the system, both the predator and prey are divided into immature individuals and mature individuals by two fixed ages. It is assumed that the immature predators cannot attack preys, and the case of the mature predators attacking the immature preys is also ignored. Based on Mawhin's coincidence degree, sufficient conditions are obtained for the existence of two positive periodic solutions of the system. An example is presented to illustrate the feasibility of the main results.     

The complex dynamics of a diffusive prey–predator model with an Allee effect in prey

This paper investigates complex dynamics of a predator–prey interaction model that incorporates: (a) an Allee effect in prey; (b) the Michaelis–Menten type functional response between prey and predator; and (c) diffusion in both prey and predator. We provide rigorous mathematical results of the proposed model including: (1) the stability of non-negative constant steady states; (2) sufficient conditions that lead to Hopf/Turing bifurcations; (3) a prior estimates of positive steady states; (4) the non-existence and existence of non-constant positive steady states when the model is under zero-flux boundary condition. We also perform completed analysis of the corresponding ODE model to obtain a better understanding on effects of diffusion on the stability. Our analytical results show that the small values of the ratio of the prey's diffusion rate to the predator's diffusion rate are more likely to destabilize the system, thus generate Hopf-bifurcation and Turing instability that can lead to different spatial patterns. Through numerical simulations, we observe that our model, with or without Allee effect, can exhibit extremely rich pattern formations that include but not limit to strips, spotted patterns, symmetric patterns. In addition, the strength of Allee effects also plays an important role in generating distinct spatial patterns.     

Turing patterns in a predator–prey model with seasonality

Many ecological systems show striking non-homogeneous population distributions. Diffusion-driven instabilities are commonly studied as mechanisms of pattern formation in many fields of biology but only rarely in ecology, in part because some of the conditions seem quite restrictive for ecological systems. Seasonal variation is ubiquitous in temperate ecosystems, yet its effect on pattern formation has not yet been explored. We formulate and analyze an impulsive reaction–diffusion system for a resource and its consumer in a two-season environment. While the resource grows throughout the ‘summer’ season, the consumer reproduces only once per year. We derive conditions for diffusion-driven instability in the system, and we show that pattern formation is possible with a Beddington–DeAngelis functional response. More importantly, we find that a low overwinter survival probability for the resource enhances the propensity for pattern formation: diffusion-driven instability occurs even when the diffusion rates of prey and predator are comparable (although not when they are equal).

     

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.     

Permanence of a predator–prey discrete system with Holling-IV functional response and distributed delays

A predator–prey discrete-time model with Holling-IV functional response and distributed delays is investigated in this paper. By using the comparison theorem of the difference equation and some analysis technique, some sufficient conditions are obtained for the permanence of the discrete predator–prey system. Two examples are given to illustrate the feasibility of the obtained result.     

On the dynamics of one-prey-n-predator impulsive reaction-diffusion predator–prey system with ratio-dependent functional response

In this paper, a one-prey-n-predator impulsive reaction-diffusion periodic predator–prey system with ratio-dependent functional response is investigated. On the basis of the upper and lower solution method and comparison theory of differential equation, sufficient conditions on the ultimate boundedness and permanence of the predator–prey system are established. By constructing an appropriate auxiliary function, the conditions for the existence of a unique globally stable positive periodic solution are also obtained. Examples and numerical simulations are presented to verify the feasibility of our results. A discussion is conducted at the end.     

Structural sensitivity and resilience in a predator–prey model with density-dependent mortality

Numerous formulations with the same mathematical properties can be relevant to model a biological process. Different formulations can predict different model dynamics like equilibrium vs. oscillations even if they are quantitatively close (structural sensitivity). The question we address in this paper is: does the choice of a formulation affect predictions on the number of stable states? We focus on a predator–prey model with predator competition that exhibits multiple stable states. A bifurcation analysis is realized with respect to prey carrying capacity and species body mass ratio within range of values found in food web models. Bifurcation diagrams built for two type-II functional responses are different in two ways. First, the kind of stable state (equilibrium vs. oscillations) is different for 26.0–49.4% of the parameter values, depending on the parameter space investigated. Using generalized modelling, we highlight the role of functional response slope in this difference. Secondly, the number of stable states is higher with Ivlev's functional response for 0.1–14.3% of the parameter values. These two changes interact to create different model predictions if a parameter value or a state variable is altered. In these two examples of disturbance, Holling's disc equation predicts a higher system resilience. Indeed, Ivlev's functional response predicts that disturbance may trap the system into an alternative stable state that can be escaped from only by a larger alteration (hysteresis phenomena). Two questions arise from this work: (i) how much complex ecological models can be affected by this sensitivity to model formulation? and (ii) how to deal with these uncertainties in model predictions?     

Computing the heteroclinic bifurcation curves in predator–prey systems with ratio-dependent functional response

Predator-prey models with Michaelis-Menten-Holling type ratio- dependent functional response exhibit very rich and complex dynamical behavior, such as the existence of degenerate equilibria, appearance of limit cycles and heteroclinic loops, and the coexistence of two attractive equilibria. In this paper, we study heteroclinic bifurcations of such a predator-prey model. We first calculate the higher order Melnikov functions by transforming the model into a Hamiltonian system and then provide an algorithm for computing higher order approximations of the heteroclinic bifurcation curves.     

Analyzing the spatial dynamics of a prey–predator lattice model with social behavior

A lattice prey–predator model is studied. Transition rules applied sequentially describe processes such as reproduction, predation, and death of predators. The movement of predators is governed by a local particle swarm optimization algorithm, which causes the formation of swarms of predators that propagate through the lattice. Starting with a single predator in a lattice fully covered by preys, we observe a wavefront of predators invading the zones dominated by preys; subsequent fronts arise during the transient phase, where a monotonic approach to a fixed point is present. After the transient phase the system enters an oscillatory regime, where the amplitude of oscillations appears to be bounded but is difficult to predict. We observe qualitative similar behavior even for larger lattices. An empirical approach is used to determine the effects of the movement of predators on the temporal dynamics of the system. Our results show that the algorithm used to model the movement of predators increases the proficiency of predators.     

Bistability and limit cycles in generalist predator–prey dynamics

Since generalist predators feed on a variety of prey species they tend to persist in an ecosystem even if one particular prey species is absent. Predation by generalist predators is typically characterized by a sigmoidal functional response, so that predation pressure for a given prey species is small when the density of that prey is low. Many mathematical models have included a sigmoidal functional response into predator–prey equations and found the dynamics to be more stable than for a Holling type II functional response. However, almost none of these models considers alternative food sources for the generalist predator. In particular, in these models, the generalist predator goes extinct in the absence of the one focal prey. We model the dynamics of a generalist predator with a sigmoidal functional response on one dynamic prey and fixed alternative food source. We find that the system can exhibit up to six steady states, bistability, limit cycles and several global bifurcations.     

Global analysis of a predator–prey model with variable predator search rate

We consider a modified Holling-type II predator–prey model, based on the premise that the search rate of predators is dependent on the prey density, rather than constant. A complete analysis of the global behavior of the model is presented, and shows that the model exhibits a dichotomy similar to the classical Holling-type II model: either the coexistence steady state is globally stable; or it is unstable, and then a unique, globally stable limit cycle exists. We discuss the similarities, but also important differences between our model and the Holling-type II model. The main differences are that: 1. The paradox of enrichment which always occurs in the Holling-type II model, does not always occur here, and 2. Even when the paradox of enrichment occurs, predators can adapt by lowering their search rate, and effectively stabilize the system.

     

The role of noise in a predator–prey model with Allee effect

The existence and implications of alternative stable states in ecological systems have been investigated extensively within deterministic models. However, it is known that natural systems are undeniably subject to random fluctuations, arising from either environmental variability or internal effects. Thus, in this paper, we study the role of noise on the pattern formation of a spatial predator–prey model with Allee effect. The obtained results show that the spatially extended system exhibits rich dynamic behavior. More specifically, the stationary pattern can be induced to be a stable target wave when the noise intensity is small. As the noise intensity is increased, patchy invasion emerges. These results indicate that the dynamic behavior of predator–prey models may be partly due to stochastic factors instead of deterministic factors, which may also help us to understand the effects arising from the undeniable susceptibility to random fluctuations of real ecosystems.     

Global analysis of a stoichiometric producer–grazer model with Holling type functional responses

Cells, the basic units of organisms, consist of multiple essential elements such as carbon, nitrogen, and phosphorus. The scarcity of any of these elements can strongly restrict cellular and organismal growth. During recent years, ecological models incorporating multiple elements have been rapidly developed in many studies, which form a new research field of mathematical and theoretical biology. Among these models, the one proposed by Loladze et al. (Bull Math Biol 62:1137–1162, 2000) is prominent and has been highly cited. However, the global analysis of this nonsmooth model has never been done. The aim of this paper is to provide the complete global analysis for the model with Holling type I functional response and perform a bifurcation analysis for the model with Holling type II functional response.     

A predator–prey refuge system: Evolutionary stability in ecological systems

A refuge model is developed for a single predator species and either one or two prey species where no predators are present in the prey refuge. An individual’s fitness depends on its strategy choice or ecotype (predators decide which prey species to pursue and prey decide what proportion of their time to spend in the refuge) as well as on the population sizes of all three species. It is shown that, when there is a single prey species with a refuge or two prey species with no refuge compete only indirectly (i.e. there is only apparent competition between prey species), that stable resident systems where all individuals in each species have the same ecotype cannot be destabilized by the introduction of mutant ecotypes that are initially selectively neutral. In game-theoretic terms, this means that stable monomorphic resident systems, with ecotypes given by a Nash equilibrium, are both ecologically and evolutionarily stable. However, we show that this is no longer the case when the two indirectly-competing prey species have a refuge. This illustrates theoretically that two ecological factors, that are separately stabilizing (apparent competition and refuge use), may have a combined destabilizing effect from the evolutionary perspective. These results generalize the concept of an evolutionarily stable strategy (ESS) to models in evolutionary ecology. Several biological examples of predator–prey systems are discussed from this perspective.