ISRN Biomathematics

Volume 2013, Article ID 637640, 12 pages

http://dx.doi.org/10.1155/2013/637640

## Global Dynamics of an Exploited Prey-Predator Model with Constant Prey Refuge

^{1}Department of Mathematics, Sree Chaitanya College, Habra, North 24 Parganas 743268, India^{2}Department of Mathematics, Bengal Engineering and Science University, Shibpur, Howrah 711103, India

Received 27 May 2013; Accepted 8 July 2013

Academic Editors: M. T. Figge, M. José, A. A. Polezhaev, and J. H. Wu

Copyright © 2013 Uttam Das et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

This paper describes a prey-predator model with Holling type II functional response incorporating constant prey refuge and harvesting to both prey and predator species. We have analyzed the boundedness of the system and existence of all possible feasible equilibria and discussed local as well as global stabilities at interior equilibrium of the system. The occurrence of Hopf bifurcation of the system is examined, and it was observed that the bifurcation is either supercritical or subcritical. Influences of prey refuge and harvesting efforts are also discussed. Some numerical simulations are carried out for the validity of theoretical results.

#### 1. Introduction

Prey-predator models are of great interest to researchers in mathematics and ecology because they deal with environmental problems such as community’s morbidity and how to control it and optimal harvest policy to sustain a community. In the physical sciences, generic models can be constructed to explain a variety of phenomena. However, in the life sciences, a model only describes a particular situation. So a variety of models are needed due to the complexity of the ecosystem. Theoretical and numerical studies of these models are able to give us an understanding of the interactions that are taking place.

While investigating biological phenomena, there are many factors which affect dynamical properties of biological and mathematical models. One of the familiar nonlinear factors is the functional response. In population dynamics, a functional response of the predator to the prey density refers to the change in the density of prey attached per unit time per predator as the prey density changes [1]. Holling [2] suggested three different kinds of functional response for different kinds of species to model the phenomena of predation, which made the classical Lotka-Volterra system more realistic. The other factor which affects dynamical properties is harvesting. The subject of harvesting in predator-prey systems has been of interest to economists, ecologists, and natural resource managers for some time now. There are basically three types of harvesting reported in the literature: (i) constant harvesting where a constant number of individuals are harvested per unit time, (ii) proportional harvesting where the number of individuals harvested per unit time is proportional to the current population, and (iii) nonlinear harvesting. Most research has focused attention on optimal exploitation, guided entirely by profits from harvesting. First of all, depending on the nature of the applied harvesting strategy, the long-run stationary density of population may be significantly smaller than the long-run stationary density of a population in the absence of harvesting. Therefore, while a population can, in the absence of harvesting, be free of extinction risk, harvesting can lead to the incorporation of a positive extinction probability and therefore, to potential extinction in finite time. Secondly, if a population is subjected to a positive extinction rate, then harvesting can drive the population density to a dangerously low level at which extinction becomes a sure thing no matter how the harvester affects the population afterwards.

Effects of harvesting on various types of prey-predator models have been considered by many researchers. Mathematical modeling with harvesting renewable resources started with the studies of Clark [3, 4]. Xiao et al. [5] have investigated the dynamics of a system with constant harvesting on the predator, whereas Leard et al. [6] have studied the dynamics of ratio-dependent models that include nonconstant harvesting on the prey. Xia et al. [7] studied the combined effects of harvesting and time delay. Recently, Pahari and Kar [8] analyzed a prey-predator fishery model with harvesting where they applied a tax to regulate the fishery.

The existence of prey refuges can clearly have important effects on the coexistence of predators and their prey. According to Taylor [9], the different kinds of refuges can be arranged into three types: (a) those which provide permanent spatial protection for a small subset of the prey population, (b) those which provide temporary spatial protection, and (c) those which provide a temporal refuge in numbers, that is, decrease the risk of predation by increasing the abundance of vulnerable prey. It is also observed that refuge has a stabilizing effect on the equilibrium for a simple Lotka-Volterra model. The problem of predator-prey interactions under a prey refuge has been studied by some authors. McNair [10] showed that several kinds of refuges could exert a locally destabilizing effort and create stable, large-amplitude oscillations which would damp out if no refuge was present. McNair [11] obtained that a prey refuge with legitimate entry existing dynamics was quite capable of amplifying rather than damping predator-prey oscillation. Kar [12] proposed a predator-prey model incorporating a prey refuge and independent harvesting on either species. He showed that, using the harvesting efforts as control, it is possible to break the cyclic behavior of the system. Wang et al. [13] considered a prey-predator system where individuals from a prey fish population could hide in holes where predators were unable to enter. Ji and Wu [14] considered a predator-prey model incorporating a constant prey refuge and a constant-rate prey harvesting and showed that the system is controlled by using constant harvesting or prey refuge. Prepredator interactions with prey refuge also may be found in the work of Huang et al. [15], Wang and Pan [16], Kar et al. [17], and so forth.

The objective of this paper is to study the combined effects of harvesting and constant prey refuge on the dynamics of predator-prey model. We also examined the existence of a Hopf bifurcation and determined the conditions for which the bifurcation is either supercritical or subcritical.

González-Olivares and Ramos-Jiliberto [18] studied the following prey-predator model with constant number of prey using refuges as shown in the Figure 1: with , , where and denote the prey and predator populations, respectively, at any time and and , respectively, represent the intrinsic growth rate and carrying capacity of the prey. is a constant, which represents number of prey which seeks refuge from predation. is the death rate of the predator and is the conversion factor. The term denotes the functional response of the predator, which is termed Holling type II response function (see [2]). They examined the local stability of equilibria and existence of limit cycle. Chen et al. [19] considered the same model and examined also global stability and uniqueness of limit cycle. Also Ji and Wu [14] considered the prey-predator model with constant prey refuge and a constant-rate prey harvesting , and studied instability and global stability of the equilibria and uniqueness of limit cycle and showed also the influence of constant prey refuge and constant-rate prey harvesting.

Motivated by the paper of González-Olivares and Ramos-Jiliberto [18], we considered the following prey-predator system: where and denote the harvesting efforts for the prey and predator, respectively. and represent the catch of the prey and predator population, where and represent the catchability coefficients, respectively. , , , , , , , , and are positive constants.

Our paper is organized in the following way. In Section 2, we have discussed the boundedness of the solutions of system (3). All possible equilibria of the system (8) and the stability criterion at those equilibria and uniqueness of limit cycles at interior equilibrium are discussed in Section 3. In Section 4, we have discussed Hopf bifurcation at interior equilibrium. Influence of refuge parameter and harvesting efforts and are discussed in Section 5. Numerical simulations are given in Section 6. A brief concluding remark is given in Section 7.

#### 2. Boundedness of the System

Boundedness of a model guarantees its validity. The following theorem establishes the uniform boundedness of the model system (3).

Theorem 1. *All solutions of the system (3) which start in are uniformly bounded.*

*Proof. *Let be any solution of the system with positive initial conditions and .

Now we define the function .

Therefore, time derivative gives
Now for each , we have
We have .

Let us choose ; then the right-hand side is positive. As we assume that both and are bounded, the right-hand side is bounded for all .

Thus we choose a such that .

Applying the theory of differential inequality [20], we obtain
Thus, all solutions of the system (3) that start in are confined to the region , where
This completes the theorem.

Now for simplicity, let us introduce ; then the system (3) of equations changes to (still denote )

#### 3. 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 (8) as follows:(i)equilibrium in the absence of predator , where exists provided ,(ii)the interior equilibrium , where If holds and any of the following conditions hold:(H1), , , and ,(H2), , and ,(H3), , and ,where and , then the interior equilibrium exists.

For the stability analysis of the equilibria, let the variational matrix of the system (8) at an arbitrary point be .

At , the eigenvalues of the variational matrix are and .

Now if and , then both eigenvalues are negative and hence is locally asymptotically stable.

To determine the local stability character of the interior equilibrium , we compute the variational matrix at .

The characteristic equation of the variational matrix is given by where Routh-Hurwitz criterion states that all roots of the characteristic equation (10) have negative real parts if .

Therefore, interior equilibrium will be asymptotically stable if , and unstable if .

Theorem 2. *If exists with then is locally asymptotically stable.*

##### 3.1. Global Stability for the Interior Equilibrium

Theorem 3. *If , then the system (8) will be globally asymptotically stable around the interior equilibrium .*

*Proof. *To show the global stability of system (8), we define a Lyapunov function as follows:
where and are positive constants to be determined in the subsequent steps.

The time derivative along the trajectories of (8) is
where .

Choosing and , we have
Thus if , then . This completes the proof.

In Figure 2, we show that may be positive for some positive value of by using the numerical value of the parameters , , , , , , , , , , and .

##### 3.2. Uniqueness of Limit Cycles

It is known that for prey-predator 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 (8) in the form where , and .

Now we consider the following theorem (see [21]) regarding uniqueness of limit cycles of the previous system.

Theorem 4. *Suppose for system (15) that
**
in and . Then the previous system has exactly one limit cycle which is globally asymptotically stable with respect to the set
*

Following Theorem 4, we may state that when , system (8) has unique globally stable limit cycle. Thus we see that when the system is unstable, there exists a unique globally stable limit cycle.

#### 4. Hopf Bifurcation

To discuss Hopf bifurcation of the system (8) we take the help of the paper [22].

Let us now consider system (8) in the form where , and .

Let be the interior equilibrium of the system (18), where In this section, we examined the Hopf bifurcation occurring at taking as a bifurcation parameter.

Variational matrix of the system (18) at is where , , and .

The characteristic equation of the variational matrix is given by , where and .

Now if and , then is locally asymptotically stable. On the other hand, if and then , is unstable. Also and the Jacobian matrix has a pair of imaginary eigenvalues . Let be the roots of when is near ; then is a root of (as ; that is , satisfies the equation and also .

Therefore, as .

The system (18) undergoes a Hopf bifurcation at when . By further analysis, we determined that the bifurcation is either supercritical or subcritical by the first Lyapunov coefficient (see [23–25]) The computation of is technical, and the detailed calculations are given in the appendix and only the results are stated in the following theorem.

Theorem 5. *The system (18) undergoes a Hopf bifurcation at ; the Hopf bifurcation is supercritical and backward (subcritical and forward, resp.) if , where is defined in (23).*

#### 5. Influence of the Parameter and Harvesting Efforts and

Interior equilibrium of the system (8) is , where Notice that the scaling from system (3) to system (8) is . Let the interior equilibrium of the system (3) be , where The model (3) is a one-prey-one-predator model with constant prey refuge. Now we discuss the influence of the prey refuge on the model dynamics.

Now .

Therefore, is a strictly increasing function of .

That is, increasing the amount of prey refuge can increase prey population.

To check the effect of on predator, differentiating with respect to , we get , where .

*Case 1. *If , then is a monotonically decreasing function of .

That is, if the value of prey refuge gradually increases above threshold value , then the predator population gradually decreases.

*Case 2. *If , then is a monotonically increasing function of .

That is, if the value of prey refuge gradually decreases below threshold value , then the predator population gradually increases.

*Case 3. *If , then has a maximum value.

To construct Table 1 and Figure 3, we take , , , , , , , , , and in appropriate units. From the analysis shown we see that increasing the amount of prey refuge can increase prey population and that increasing the amount of prey refuge can increase the density of predator species, and this happened due to predator species still having enough food for predation with being small, but if the prey refuge is larger than a threshold, that is, as the prey refuge becomes large enough, then the increasing amount of prey refuge can decrease predator species, and this happened due to the loss of food for predator species.

*Influence of the Harvesting Efforts E*_{1}* and E*_{2}*.* Interior equilibrium of the system (8) is , where
At first let and be fixed.

It is observed that is independent of the parameter , whereas the value of depends on . Therefore has no effect on the interior equilibrium level on .

To check the effect of on predator, differentiating with respect to , we get

This shows that is a strictly decreasing function of .

To construct Table 2 and Figure 4, we take , , , , , , , , , and in appropriate units. From the analysis shown we see that as the harvesting effort increases, the prey species remain unchanged but predator species decrease and this happens due to loss of food for predator species and goes to extinction when is large.

Now let and be fixed; then we have Therefore is a strictly increasing function of .

Also Since , 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 decrease predator population.

To construct Table 3 and Figure 5, we take , , , , , , , , , and in appropriate units. From the analysis shown 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 will decrease predator species and go to extinction of the predator species when is large.

#### 6. 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. Let us consider the parameters of the system as , , , , , , , , , and in appropriate units. For these value of parameters, we get the critical value of as . Thus it is easy to verify that for this set of parameters, the system (8) is locally asymptotically stable around its interior equilibrium for and is unstable for . Thus for the system (8) undergoes a Hopf bifurcation. Now for , we have interior equilibrium which is asymptotically stable (see Figure 6), but , and the interior equilibrium is unstable (see Figure 7). Thus taking as a control parameter, it is possible to drive the prey-predator system to require equilibrium and to prevent the cycle behaviour of the system. From Figures 8, 9, 10, and 11, we see that and may also be used as controls for the system (8); Hopf bifurcation occurs when (here , ) and (here , ). For the previous values of parameters and , we obtained one value of say and . Thus we may conclude that the Hopf bifurcation around the interior equilibrium is supercritical and backward.

#### 7. Concluding Remarks

This paper deals with a prey-predator model with Holling type II functional response incorporating a constant prey refuge and independent harvesting in either species. Oscillatory behavior and existence of limit cycles in harvested prey-predator system are common in nature. It is noted that constant prey refuge plays an important role in the dynamics of the proposed model system. It is also observed from the obtained results that constant prey refuge can cause an unstable equilibrium to become stable, and even a simple Hopf bifurcation occurred when the parameter passes through its critical value. There exists a threshold value of , such that for the prey refuge smaller than this threshold, increasing the amount of prey refuge can increase the predator population, and if the prey refuge is larger than the threshold, increasing the amount of prey refuge can decrease the predator population. We have proved that exactly one stable limit cycle occurs when the positive equilibrium is unstable. We also determined the critical value of at which Hopf bifurcation occurs and observed that the bifurcation is supercritical and backward. It was also found that it is possible to control the system in such a way that the system approaches a required state, using the efforts and as controls.

Our analytical results and numerical simulation also indicate that dynamic behavior of the model not only depends on the prey refuge parameter but also depends on harvesting efforts and . Hence it is possible to control the system in such a way that the system approaches a required state, using the harvesting efforts and or prey refuge as controls.

In our model, we have considered the catch-rate function based on catch-per-unit-effort hypothesis; that is, ( and denote the prey or predator population and harvesting effort, resp.). But this type of catch-rate function embodies some defects in that (i) it assumes random search for fish, (ii) it assumes equal likelihood of being captured for every fish, (iii) there is unbounded linear increase of with respect to for fixed , and (iv) there is unbounded linear increase with respect to for a fixed . These unrealistic features can largely be removed by adopting the alternative functional form , where and are positive constants, but we leave it for our future research work. The entire study of the paper is mainly based on the deterministic framework. On the other hand, it will be more realistic if it is possible to consider the model system in the stochastic environment due to some ecological fluctuations and other factors. Thus, a future research problem would be considered in stochastic environment.

#### Appendix

#### Detailed Calculation of Formula (23)

First we translate the equilibrium to the origin by translation , . (Still denote and by and , resp.). Thus the system (18) becomes We write the system (A.1) as follows: where and .

Now compute the Taylor expansion of related functions where Then the system (A.1) becomes where Now we define a matrix where and ; then By using the linear transformation we have where Then system (A.5) becomes where and are defined in (21), and where Rewrite the system (A.12) in a polar coordinate form as: Then the Taylor expansion of (A.15) at yields In order to determine the stability of the periodic solution, we need to calculate the sign of the coefficient , which is given by where all partial derivatives are evaluated at the bifurcation point; that is, .

Since is linear in both and , we have that when, , , also .

At , we have by simple calculation that Now since .

#### References

- S. Zhang, L. Dong, and L. Chen, “The study of predator-prey system with defensive ability of prey and impulsive perturbations on the predator,”
*Chaos, Solitons and Fractals*, vol. 23, no. 2, pp. 631–643, 2005. View at Publisher · View at Google Scholar · View at Scopus - C. S. Holling, “The functional response of predator to prey density and its role mimicry and population regulation,”
*Memoirs of the Entomological Society of Canada*, vol. 45, pp. 3–60, 1965. View at Google Scholar - C. W. Clark,
*Bioeconomic Modelling and Fisheries Management*, Wiley, New York, NY, USA, 1985. - C. W. Clark,
*Mathematical Bioeconomics: The Optimal Management of Renewable Resource*, John Wiley and Sons, New York, NY, USA, 2nd edition, 1990. - D. Xiao, W. Li, and M. Han, “Dynamics in a ratio-dependent predator-prey model with predator harvesting,”
*Journal of Mathematical Analysis and Applications*, vol. 324, no. 1, pp. 14–29, 2006. View at Publisher · View at Google Scholar · View at Scopus - B. Leard, C. Lewis, and J. Rebaza, “Dynamics of ratio-dependent predator prey models with nonconstant harvesting,”
*Discrete and Continuous Dynamical Systems, Series S*, vol. 1, pp. 303–315, 2008. View at Google Scholar - J. Xia, Z. Liu, R. Yuan, and S. Ruan, “The effects of harvesting and time delay on predator-prey systems with Holling type II functional response,”
*SIAM Journal on Applied Mathematics*, vol. 70, no. 4, pp. 1178–1200, 2009. View at Google Scholar - U. K. Pahari and T. K. Kar, “Conservation of a resource based fishery through optimal taxation,”
*Nonlinear Dynamics*, vol. 72, pp. 591–603, 2013. View at Google Scholar - R. J. Taylor,
*Predation*, Chapman and Hall, New York, NY, USA, 1984. - J. N. McNair, “The effects of refuges on predator-prey interactions: a reconsideration,”
*Theoretical Population Biology*, vol. 29, no. 1, pp. 38–63, 1986. View at Google Scholar · View at Scopus - J. N. McNair, “Stability effects of prey refuges with entry-exit dynamics,”
*Journal of Theoretical Biology*, vol. 125, no. 4, pp. 449–464, 1987. View at Google Scholar · View at Scopus - T. K. Kar, “Modelling and analysis of a harvested prey-predator system incorporating a prey refuge,”
*Journal of Computational and Applied Mathematics*, vol. 185, no. 1, pp. 19–33, 2006. View at Publisher · View at Google Scholar · View at Scopus - H. Wang, W. Morrison, A. Singh, and H. Weiss, “Modeling inverted biomass pyramids and refuges in ecosystems,”
*Ecological Modelling*, vol. 220, no. 11, pp. 1376–1382, 2009. View at Publisher · View at Google Scholar · View at Scopus - L. Ji and C. Wu, “Qualitative analysis of a predator-prey model with constant-rate prey harvesting incorporating a constant prey refuge,”
*Nonlinear Analysis: Real World Applications*, vol. 11, no. 4, pp. 2285–2295, 2010. View at Publisher · View at Google Scholar · View at Scopus - Y. Huang, F. Chen, and L. Zhong, “Stability analysis of a prey-predator model with holling type III response function incorporating a prey refuge,”
*Applied Mathematics and Computation*, vol. 182, no. 1, pp. 672–683, 2006. View at Publisher · View at Google Scholar · View at Scopus - J. Wang and L. Pan, “Qualitative analysis of a harvested predator-prey system with Holling-type III functional response incorporating a prey refuge,”
*Advances in Difference Equations*, vol. 96, pp. 1–14, 2012. View at Google Scholar - T. K. Kar, A. Ghorai, and S. Jana, “Dynamics consequences of prey refuges in a two predator one prey system,”
*Journal of Biological Systems*, vol. 21, no. 2, Article ID 1350013, 28 pages, 2013. View at Google Scholar - E. González-Olivares and R. Ramos-Jiliberto, “Dynamic consequences of prey refuges in a simple model system: more prey, fewer predators and enhanced stability,”
*Ecological Modelling*, vol. 166, no. 1-2, pp. 135–146, 2003. View at Publisher · View at Google Scholar · View at Scopus - L. Chen, F. Chen, and L. Chen, “Qualitative analysis of a predator-prey model with Holling type II functional response incorporating a constant prey refuge,”
*Nonlinear Analysis: Real World Applications*, vol. 11, no. 1, pp. 246–252, 2010. View at Publisher · View at Google Scholar · View at Scopus - G. Birkoff and G. C. Rota,
*Ordinary Differential Equations*, Ginn, Cambridge, UK, 1982. - Y. Kuang and H. I. Freedman, “Uniqueness of limit cycles in Gause-type models of predator-prey systems,”
*Mathematical Biosciences*, vol. 88, no. 1, pp. 67–84, 1988. View at Google Scholar · View at Scopus - J. Wang, J. Shi, and J. Wei, “Predator-prey system with strong Allee effect in prey,”
*Journal of Mathematical Biology*, vol. 62, no. 3, pp. 291–331, 2011. View at Publisher · View at Google Scholar · View at Scopus - S. Wiggins,
*Introduction to Applied Nonlinear Dynamical Systems Chaos*, Texts in Applied Mathematics, Springer, New York, NY, USA, 2nd edition, 2003. - B. D. Hassard, N. D. Kazarinoff, and Y. Wan,
*Theory and Applications of Hopf Bifurcation*, vol. 41 of*London Mathematical Society Lecture Note Series*, Cambridge University Press, Cambridge, UK, 1981. - Y. A. Kuznetsov,
*Elements of Applied Bifurcation Theory*, vol. 112 of*Applied Mathematical Sciences*, Springer, New York, NY, USA, 2004.