Abstract

In the present paper, I study a prey-predator model with multiple time delays where the predator population is regarded as generalist. For this regard, I consider a Holling-Tanner prey-predator system where a constant time delay is incorporated in the logistic growth of the prey to represent a delayed density dependent feedback mechanism and the second time delay is considered to account for the length of the gestation period of the predator. Predator’s interference in predator-prey relationship provides better descriptions of predator's feeding over a range of prey-predator abundances, so the predator's functional response here is considered to be Type II ratio-dependent. In accordance with previous studies, it is observed that delay destabilizes the system, in general, and stability loss occurs via Hopf bifurcation. There exist critical values of delay parameters below which the coexistence equilibrium is stable and above which it is unstable. Hopf bifurcation occurs when delay parameters cross their critical values. When delay parameters are large enough than their critical values, the system exhibits chaotic behavior and this abnormal behavior may be controlled by refuge. Numerical computation is also performed to validate different theoretical results. Lyapunov exponent, recurrence plot, and power spectral density confirm the chaotic dynamical behaviors.

1. Introduction

In ecology, predation describes a biological interaction where a predator (an organism that is hunting) feeds on its prey (the organism that is attacked). Predators may or may not kill their prey prior to feeding on them, but the act of predation often results in the death of its prey and the eventual absorption of the prey’s tissue through consumption. Predators can have profound impacts on the dynamics of their prey that depend on how predator consumption is affected by prey density (the predator’s functional response). A generalist species is able to thrive in a wide variety of environmental conditions and can make use of a variety of different resources (e.g., a heterotroph with a varied diet). Omnivores are usually generalists. Herbivores are often specialists, but those that eat a variety of plants may be considered generalists. Consumption by a generalist predator is expected to depend on the densities of all its major prey species (its multispecies functional response, or MSFR), but most studies of generalists have focused on their functional response to only one prey species [14]. Recently, many authors have explored the dynamics of a class of the semi-ratio-dependent [58]/Holling-Tanner [912]/Leslie-Grower [13, 14] predator-prey systems with functional responses as follows: where and are, respectively, the prey and predator densities at time . The prey population is assumed to grow logistically to its carrying capacity with intrinsic growth rate in absence of predator. is the maximal per capita prey consumption rate, and is the amount of prey necessary for the relative biomass growth rate of the predator to be half its maximum; the predator consumes the prey according to the functional response and grows logistically with growth rate and carrying capacity proportional to the population size of the prey (or prey abundance). is the measure of the food quality that the prey provides, which is converted to predator birth. For more background and biological adjustments of system (1), one can see [7, 15, 16] and the references cited therein.

Predator’s functional response is considered as Type II ratio-dependent [1725] because a ratio-dependent predator-prey model does not show the so called paradox of enrichment [17, 26, 27] and biological control paradox [28]. The simplification of the food web structure due to predator interference allows qualitative predictions concerning the response of a food web to an external perturbation [25, 2932]. Ginzburg and Akcakaya [33] demonstrated a positive relationship between the abundances of all trophic levels and concluded that the ratio-dependent models provide a better representation of predator-prey interactions than prey dependent models. Predator’s functional response, defined as the amount of prey catch per predator per unit of time, is affected by the structure of prey habitat and predator’s hunting ability [34, 35]. Anderson [36] and Johnson [37] showed that the dynamics of a local population largely depend upon attributes of the local habitats. In addition, a spatial refuge protects a constant proportion of prey from predation. Mite predator-prey interactions often exhibit spatial refugia which afford the prey some degree of protection from predation and reduce the chance of extinction due to predation. Maynard Smith [38] shows that the presence of a constant proportion refuge does not alter the dynamical stability of the neutrally stable Lotka-Volterra model, while a constant number refuge of any size replaces the neutrally stable behaviour with a stable equilibrium. Hassel [39] shows that adding a large refuge to a model, which in the absence of a refuge exhibits divergent oscillations, replaces the oscillatory behaviour with a stable equilibrium. These mathematical models and a number of experiments indicate that refugia have a stabilizing effect on predator-prey interactions, but, as Taylor [40] has mentioned, it would be an over simplification to assume this is always the case. Kar [41] and Jana [42] proposed and studied a predator-prey system where prey population is subjected to refuge from their predator population. These papers lead a model by incorporating a refuge protecting of the prey, where is the density of prey population and is constant measures of the degree or strength of prey refuge. This leaves of the prey available to the predator, so the Holling-Tanner predator-prey model with ratio-dependent functional response with prey refuge is given by

Delay is frequently used in a predator-prey model to make the model biologically more realistic. Recently, researchers are using more than one delay to study the effect of past history of the system populations [4347]. In the second phase of this study, I consider two delays in the model system (2). One discrete delay is considered in the specific growth rate of prey to incorporate the effect of density dependence feedback mechanism which takes units of time to respond to changes in the prey population [48]. The second delay is incorporated in the negative feedback of the predator density [4952]. I thus obtain the following multidelayed predator-prey model: The model system (3) has to be studied with the following initial conditions: The organization of the paper is as follows. Section 2 deals with the analysis of the model system (3). Numerical studies are given in Section 3, and, finally, a brief discussion is presented in Section 4.

2. Mathematical Results of System (3)

The system (3) has only one interior equilibrium point given by , where and . The equilibrium point is biologically feasible if . Linearizing the system (3) at , I get where The corresponding characteristic equation is given by where I now discuss the following cases.

Case 1 (). In this case, (7) becomes All roots of (9) will have negative real parts if and only if

(H1)

Theorem 1. The interior equilibrium of the system (3) exists and becomes locally asymptotically stable in absence of delays if and hold.

Case 2 (). In this case, the characteristic equation (7) becomes Let be a root of (11). Then I have This leads to It follows that (12) has no positive roots if the following conditions are satisfied:
(H2)
Hence, all roots of (13) will have negative real parts when if conditions of Theorem 1 and (H2) are satisfied.
Let
(H3) If Theorem 1 and (H3) hold, then (13) has a unique positive root . Substituting into (12), I have where is the unique positive root of (13).
Let
(H4) If (H1) and (H4) hold, then (13) has two positive roots and . Substituting into (12), we obtain If is a root of (11) satisfying , resp.) and (or , resp.), I obtain by hypothesis (H4).
Similarly, I can show that

Theorem 2. Assume that and conditions of Theorem 1 are satisfied, then the equilibrium is locally asymptotically stable for and unstable for . Furthermore, the system (3) undergoes a Hopf bifurcation at when , where

Case 3 (). In this case, I allow a gestation time period for the predator and also a constant time delay affecting the density dependent growth rate of the prey. I fix at some value from its stability range and regard as a free parameter. I also assume that the model parameters are such, that hypothesis (H1) holds. Let with be a root of (7). Then, where Note that is always negative. I define Then it is easy to verify that and , as . Thus, (22) has finitely many positive roots, denoted by . For every fixed , where , there exists a sequence , where with such that (22) holds. Let . When , (7) has a pair of pure imaginary roots for . In the following, I assume that
(H5) Therefore, by the general Hopf bifurcation theorem of functional differential equations, see, for example, [53], I obtain the following result considering the change of stability of system (3).

Theorem 3. Suppose that the parameters in system (3) are such that hypotheses and ( hold true, and . Then the coexistence equilibrium is locally asymptotically stable when , and it is unstable when . Moreover, Hopf bifurcation occurs when .

Case 4 (). For this choice of the delay parameters, I summarize my results in the following theorem. The proof follows similar arguments as the stability Theorem 2. in Case 2 above.

Theorem 4. Assume that and the conditions in hypothesis hold true. Then the equilibrium is locally asymptotically stable for and unstable for . Furthermore, the system (3) undergoes Hopf bifurcation when , where and is the unique positive root of the polynomial as follows:

Case 5 ( and is within its stability range). This choice of parameter regime corresponds to the biologically interesting case, when the predator individuals have a fixed gestation period and the growth of the prey population is affected by delayed density dependent mechanisms. The proof of the main result follows similar lines to that in Case 3 and I only summarize the stability results in the following theorem.

Theorem 5. Assume that the model parameters are such that the hypotheses in hold true, and . Then the coexistence equilibrium is asymptotically stable when , and it is unstable when . Hopf bifurcation occurs when , where and with

3. Numerical Computations

I perform numerical computations to observe various dynamics of the coexistence equilibrium for both the systems (2) and (3). I consider the fixed parameter values as , and . With these parameter values . Here, and , and initial value is considered as for each simulation. First I consider the nondelayed system (2). In this case, the system is asymptotically stable (Figure 1).

I now study the dynamic behavior of the delay-induced system (3). When (Case 2), one can compute from Theorem 2 that , , and . Therefore, the coexisting equilibrium is asymptotically stable for (Figure 2(b)) and unstable for (Figure 2(c)). When , the system (3) undergoes a Hopf bifurcation at . System behaviors in Case 2 can be demonstrated more prominently if I plot the bifurcation diagram in the three-dimensional space . Figure 2(a) shows that, when , the coexistence equilibrium is stable for but the instability sets in when .

For Case 3, I take any value of from its stability range , say that , and consider as a parameter. One can compute the value of as and the corresponding critical value of as . Thus, for fixed stable value of , the system (3) exhibits stable behavior around for and unstable oscillatory behavior for (Figure 3). A Hopf bifurcation occurs when takes the critical value .

In Case 4, when , one can compute from Theorem 4, and . One can also verify that and . Therefore, the coexisting equilibrium is asymptotically stable for , unstable for , and when , the system (3) undergoes a Hopf bifurcation at . These results can be observed from the bifurcation diagram Figure 4.

In the last case (Case 5), I take any value of from its stability range , say that , and consider as a parameter. As before, I calculate and the corresponding critical value . The bifurcation diagram Figure 5 demonstrates that the system (2) is locally asymptotically stable around when and unstable when . These results can be observed from the bifurcation diagram Figure 5.

3.1. Chaotic Dynamics

I have analytically studied the local behavior of the system (3) when the delay parameters are within or slightly above their critical values. One interesting topic in the delay-induced system is to study the dynamical behavior of the system when the delay parameters are far away from their critical values, or they assume large values. To observe the system dynamics, I consider and , both are beyond their stability range, and then the system shows chaotic behavior (Figure 6).

In order to characterize the irregular behavior, I perform the standard numerical diagnostics, for example, Lyapunov exponent, sensitivity of the solutions to initial condition, power spectral density, and recurrence plot. Figure 7(a) shows the Lyapunov exponent of with respect to the delay in a smaller range [0.45, 0.539] for fixed . Positive Lyapunov exponent indicates that the system exhibits chaotic behavior [42, 54, 55]. Sensitivity of the system trajectories to the initial values is shown in Figure 7(b). To show this, I compute the error , where and . This phenomena ensures the chaotic nature of the system [42, 54, 56]. The power spectra of the predator population are presented in Figure 7(c). The irregular broad peaks of this figure are indicative of chaos and randomness [42, 54, 57]. The recurrence plot of the system for the parameter values as in Figure 6 is represented in Figure 7(d). The random points on the time-time plane ensure that more or less the same values of phase trajectories take place without any rhythm, indicating the chaotic nature of the system [42, 54, 58].

3.2. Biological Control of Chaotic Dynamics

Many ecological situations show that their unstable or abnormal behavior may be controlled by any external perturbation [14, 21] or habitat structure [34, 3640]. Recently Jana [42] shows that chaotic behavior of a discrete predator-prey system at low level of prey refuge must be controlled to its stable coexistence by increasing the degree of prey refuge. Parameter set of Figure 6 depicts the chaotic behavior of the system (3). Standing at this situation, if I increase the intensity of degree of prey refuge (from to ), then again the system (3) settles down to its stable coexistence (Figure 8).

4. Discussion

Effect of prey refuge has an extended impact on the predator-prey models. However, both field and laboratory experiments confirm that intensity of prey refuge reduces predation rates by decreasing encounter rates between predator and prey. On the other hand, a predator-prey model becomes more realistic in presence of different delays which are unavoidable elements in physiological and ecological processes. In this paper, I have studied a multidelayed predator-prey model where the prey species is subject to partial refuge from predator population. Also I have considered here that the predator is a generalist type. A time delay is considered in the growth rate of the prey to represent density dependent feedback mechanism, and a second delay is introduced to account for the gestation time of the predator. The objective was to study the role of the effects of the interaction of multiple constant time delays and the prey refuge on the dynamical behavior of a predator-prey system. To have a better understanding of the complexities in natural systems, I can construct larger systems of differential equations containing more parameters. A second approach that is gaining prominence is the inclusion of time delay terms in the differential equations. Complex dynamical behaviour arises as a consequence of time delay in a biological system (with significant time delay) which may exhibit limited cycle oscillation and chaos. In case of the multidelayed system, I obtained sufficient conditions in terms of the system parameters for the stability of the coexistence equilibrium. I observed that the coexistence equilibrium is locally asymptotically stable if the delay parameters are lower than some critical values and it is unstable when the parameter values exceed these critical values. One important question in the multidelayed predator-prey system is to study the behavior of the system when the delay parameter is significantly large and hence it is far away from its critical value. As the estimated length of delays to preserve stability and the critical length of time delays for Hopf bifurcation are dependent on the system’s parameters, it is possible to impose some control, which will prevent the possible abnormal oscillation in the population density. To explore the behavior of the system in this case, I performed extensive numerical simulations. My simulation results show that a refuged predator-prey system with multiple delays may exhibit different interesting (e.g., chaotic) behavior when the delay is large enough. Finally I showed that the chaotic characteristic of the multidelayed system was controlled by increasing the intensity of the prey refuge.

Conflict of Interests

The author declares that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

This research is supported by UGC (Dr. D. S. Kothari Postdoctoral Fellowship), India; no. F.4-2/2006(BSR)/13-1004/2013(BSR).