Bifurcation Analysis of a Delayed Predator-Prey Model with Holling Type III Functional Response and Predator Harvesting

This paper tries to highlight a delayed prey-predator model with Holling type III functional response and harvesting to predator species. In this context, we have discussed local stability of the equilibria, and the occurrence of Hopf bifurcation of the system is examined by considering the harvesting effort as bifurcation parameter along with the influences of harvesting effort of the system when time delay is zero. Direction of Hopf bifurcation and the stability of bifurcating periodic solutions are also studied by applying the normal form theory and the centermanifold theorem. Lastly some numerical simulations are carried out to draw for the validity of the theoretical results.


Introduction and Model Description
Differential equation models for interactions between species are one of the classical applications of mathematics to biology, dating back to the first half of this century.The development and use of analytical techniques and the growth of computer power have progressively improved our understanding of these types of models.
The study of population dynamics with harvesting is a subject of mathematical bioeconomics, which in turn related to the optimal management of renewable resources, Clark [1].Generally the concept of optimal resource management is based on the standard cost benefit criterion which maximizes present values of net economic revenues.This criterion is relevant to both private and public management decisions, although the specification of costs and benefits are not necessarily the same in both cases.
Regulation of exploitation of biological resources has become a problem of major concern nowadays in view of the dwindling resource stocks and the deteriorating environment.Exploitation reduces the biomass of the concerned species, exhibits oscillation, and even causes extinction of some other species.Rosenweig-MacAurtho model experiences oscillation under selective effort of Kar and Ghosh [2].Legovic et al. [3] have concluded that harvesting the prey species at maximum level causes the extinction of the predator species in traditional prey-predator system.Kar and Ghosh [4] and Ghosh and Kar [5] show that harvesting the prey species at maximum sustainable yield (MSY) level never causes the extinction of the predator species in both ratio-dependent and Holling-Tanner prey-predator systems.It is shown that harvesting the prey species at MSY level may or may not drive the predator population to extinction if intraspecific competition is present among the predator species, Kar and Ghosh [4].More recently, Ghosh and Kar [6] somehow managed an example to show that prey harvesting at MSY level may be sustainable in tritrophic food chain model.In conclusion, it is reasonable that harvesting the prey species to achieve MSY will cause the extinction of the predator if predation process follows the Holling type functional response.As MSY policy has a severe impact on other tropic levels in ecosystem, suitable management tools are required for smooth running of sustainable fishing.Creating marine protected areas, maximizing the economic rent, and taxation are some good examples to protect fish population from extinction, over exploitation, as well as to achieve desired stock levels for future generation.Kar and Matsuda [7] and Kar and Ghosh [8] have shown that marine protection can have positive impact from biological stock 2 Journal of Nonlinear Dynamics assessment but virtually reduces the economic rent of the fishermen.Recently, Kar and Das [9] have analyzed a preypredator model with harvesting where they introduce a tax to regulate the fishery.
It has been noticed over the years by many authors that time delays have come to play an important role in almost all branches of science, for example, ecology and biology.The time delay is considered in the population dynamics when the rate of change of the population is not only a function of the present population but also depends on the past population.Delay is frequently used in a predator-prey model to represent the biological process more accurately.
Wangersky and Cunningham [10] have considered a predator-prey system with time delay and studied the effect of time delay on the system.Martin and Ruan [11] have studied the combined effects of prey harvesting and delay on the dynamics of predator-prey systems and focused on three very well-studied delay predator-prey models.For first and third models they were shown that the time delay could cause not only instability and oscillations but also the switching of stabilities, while the prey harvesting changes only the equilibrium values but not the properties of solutions.But in the second model, the time delay induces instability and bifurcation but there is no switching of stabilities; however, increasing the prey harvesting level will help the system to regain its stability.This indicates that the prey harvesting has a stabilizing effect on the dynamics of the model.Chen and Changming [12] have considered a delayed model with stage structured for prey and investigated nonnegative equilibria, existence of Hopf bifurcation, the stability and direction of Hopf bifurcation by applying the normal from theory, and the centre manifold argument.Zhang et al. [13] have considered a delayed predator-prey model with constant rate prey harvesting and Holling type III functional response.They studied local stability of equilibria and existence of Hopf bifurcation at interior equilibrium.Direction of Hopf bifurcation and the stability of bifurcating periodic solutions are also studied by applying the normal form theory and the center manifold theorem.Also Kar [14], Kar and Matsuda [15], Kar and Pahari [16,17], Kar and Ghorai [18], Kuang [19], and some other authors have discussed delayed predator-prey system.
The purpose of the work is to illustrate the combined effects of harvesting and delay on the dynamics of predatorprey system.
In this paper, we consider the following prey-predator system: with initial conditions () = (), () = (), () ≥ 0, () ≥ 0,  ∈ [−, 0], (0) > 0, and (0) > 0.Here  and  denote the prey and predator populations, respectively, at time . and , respectively, represent the intrinsic growth rate and carrying capacity of the prey. is the maximum capture rate of the prey,  is the death rate of the predator, and  is the conversion factor.The term  2 /( +  2 ) denotes the functional response of the predator, which is termed as Holling type III response function (see [20]).We assume that the predator population is harvested according to catch-perunit effort (CPUE) hypothesis [1] which describes that catchper-unit effort is proportional to the stock level.Therefore, the harvest at time  is ℎ  = , where  is the catchability coefficient and  is the harvesting effort.The time delay  ≥ 0 is a constant based on the assumption that the change rate of the predators depends on the number of prey and of predators present at some previous time., , , , , , , , and  are positive constants.
For  = 0, system (1) becomes ( Our paper is organized in the following way: existence of equilibria of the system (2), the local stability criterion, and uniqueness of limit cycles at interior equilibrium are discussed in Section 2. Influence of harvesting effort  is discussed in Section 3. In presence of delay we studied stability and Hopf bifurcation in Section 4. Direction of Hopf bifurcation and the stability of bifurcating periodic solutions are also studied by applying the normal form theory and the center manifold theorem in Section 5. Numerical simulations are given in Section 6.A brief concluding remark is given in Section 7.

The Steady States and Their Stability
We now study the existence and nature of the steady states.Particularly we are interested in the interior equilibrium of the system.To begin with, we list all possible steady states of the system (2) as follows:  0 (0, 0),  1 (, 0), and the interior equilibrium  * ( * ,  * ), where If  > + holds, then the interior equilibrium  * ( * ,  * ) exists.
For the stability analysis of the equilibria, let the variational matrix of the system (2) at an arbitrary point (, ) be (, ).
The characteristic equation of the variational matrix ( * ,  * ) is given by where Routh-Hurwitz criterion states that all roots of the characteristic equation ( 4) have negative real parts if  1 > 0.
is locally asymptotically stable.

Uniqueness of Limit Cycles.
It is known that, for preypredator systems, existence and stability of a limit cycle are related to the existence and stability of a positive equilibrium.If the limit cycles do not exist, in this case the equilibrium is globally asymptotically stable.On the other hand if the positive equilibrium exists and is unstable, there must occur at least one limit cycle.
Let us consider system (2) in the form where () = (1 − (/)), ℎ() =  2 /( +  2 ), and () =  2 /( +  2 ).Now we consider the following Theorem 4.2 of Kuang and Freedman [21] regarding uniqueness of limit cycles of the above system.Theorem 2. Suppose for system (6) in 0 ≤  <  * and  * <  ≤ .Then the above system has exactly one limit cycles which is globally asymptotically stable with respect to the set Following Theorem 2, we may state that when  1 < 0, then system (6) has unique globally stable limit cycle.Thus we see that when the system is unstable, there exists unique globally stable limit cycle.

Bifurcation Analysis for the Parameter 𝐸.
The characteristic equation ( 4) has two purely imaginary roots if and only if  1 = 0 and  2 > 0 for some values of , say  * .Therefore, there is only one value of  * at which we have a Hopf bifurcation.Since  * and  * contain , it is not easy to find  explicitly.So we are trying to find  =  * numerically by solving the system of equations with suitable values of the parameters: Therefore the system undergoes a Hopf bifurcation at  =  * .The system is stable for  >  * and unstable  <  * .

Influence of the Harvesting Effort 𝐸
Let the interior equilibrium of the system (2) be  * ( * ,  * ), where Now, Therefore  * is strictly increasing function of .Also Since  * / > 0, therefore Thus as the harvesting effort  increases, the predator population increases when  is smaller than the threshold value.But if the harvesting effort  gradually increases above the threshold value, that is, as the harvesting effort  becomes large enough, then the increasing amount of the harvesting effort  can decreases predator population.
To construct Figure 1, we take  = 0.8,  = 1,  = 0.8,  = 0.9,  = 0.3,  = 18, and  = 0.5 in appropriate units.From the above analysis we see that, as the harvesting effort  increases, the prey and predator species increase.But if the harvesting effort  increases larger than a threshold value, that is, as the harvesting effort  becomes large enough, then increasing amount of the harvesting effort  can increase the prey species but decrease predator species and go to extinction of the predator species when  is large.

Stability and Hopf Bifurcation in the Presence of Delay
We will investigate the dynamics of delay system (1).Let ( * ,  * ) be the only interior equilibrium of the system (1), and let  =  −  * ,  =  −  * be the perturbed variables.After removing the nonlinear terms, we obtain the linear variational system, by using equilibrium conditions as where The characteristic equation of the linearized system ( 14) is the transcendental equation of the following form: where  = −( 1 +  5 ),  =  1  5 ,  = − 4 , and It is known that the steady state is asymptotically stable if all roots of the characteristic equation ( 16) have negative real parts.For a nonlinear delay equation, there are two types of stabilities: absolute stability (independent of the delay) and conditional stability (depending of the delay).
Case 2 ( ̸ = 0).We want to determine if the real part of some root increases to reach zero and eventually becomes positive as  varies.If  is a root of ( 16), then by separating the real and imaginary parts, we get cos  −  sin  = −.
From (20), it follows that, if then (20) does not have positive roots.Therefore, characteristic equation ( 16) does not have purely imaginary roots.Since ( 1 ) and ( 2 ) ensure that all roots of (17) have negative real parts by Rouche's theorem, it follows that the roots of (20) have negative real part too.This is summarized as follows.
Based on the discussion above, we can obtain the following results.Theorem 6.The sign of  2 determines the direction of Hopf bifurcation: if  2 > 0 ( 2 < 0), then the Hopf bifurcation is supercritical (subcritical);  2 determines the stability of the bifurcating periodic solution; the bifurcating periodic solution is stable (unstable) if  2 < 0 ( 2 > 0);  2 determines the period of the bifurcating periodic solution; the period increases (decreases) if  2 > 0 ( 2 < 0).

Numerical Simulation
As the problem is not a case study, the real-world data are not available for this model.We, therefore, take here some hypothetical data with the sole purpose of illustrating the results that we have established in the previous sections.
(i) Let us consider the parameters of the system (2) as  = 0.2,  = 1,  = 0.8,  = 0.5,  = 10,  = 2, and  = 0.9 in appropriate units.For this value of parameters we get the critical value of  as  * = 0.0949124.Thus it is easy to verify that for this set of parameters the system (2) is locally asymptotically stable around its interior equilibrium ( * ,  * ) for  >  * and is unstable for  <  * .Thus for  =  * = 0.0949124, the system (2) undergoes a Hopf bifurcation.Now for  = 0.1, we have interior equilibrium (5.916, 0.552) which is asymptotically stable (see Figures 2 and 3), but for  = 0.05, the interior equilibrium (2.236, 0.463) is unstable (see Figures 4 and 5).Thus taking  as control parameter, it is possible to drive the prey-predator system to require equilibrium and to prevent the cycle behaviour of the system.

Concluding Remarks
This paper deals with a delayed prey-predator model with Holling type III functional response and harvesting of predator species.Oscillatory behavior and existence of limit cycles in harvested predator-prey system are common in nature.In the absence of delay, we have proved that exactly one stable limit cycle occurs when positive equilibrium is unstable.Since harvesting is associated with economic interests, our analysis shows that the harvesting of predator species plays an important role in the shaping of the dynamical behavior of the system.From Figure 1, we see that harvesting effort affects the interior equilibrium of system (2), and, as the harvesting effort increases, the interior equilibrium may be vanished, which means that the population of predator may become extinct.It is also observed from the obtained results that harvesting effort can cause an unstable equilibrium to become stable and even a simple Hopf bifurcation occurred when the parameter  passes through its critical value.Taking delay as bifurcation parameters, we see that a Hopf bifurcation occurs, whenever delay increases a critical value.Furthermore, the direction of Hopf bifurcation and the stability of the bifurcating periodic solutions are determined by the normal form theory and center manifold theorem.

Figure 1 :
Figure 1: Variation of  * and  * with the change of the harvesting effort .