## Stability and Bifurcation Analysis of Differential Equations and its Applications

View this Special IssueResearch Article | Open Access

# Stability and Hopf Bifurcation in a Delayed Predator-Prey System with Herd Behavior

**Academic Editor:**Yonghui Xia

#### Abstract

A special predator-prey system is investigated in which the prey population exhibits herd behavior in order to provide a self-defense against predators, while the predator is intermediate and its population shows individualistic behavior. Considering the fact that there always exists a time delay in the conversion of the biomass of prey to that of predator in this system, we obtain a delayed predator-prey model with square root functional response and quadratic mortality. For this model, we mainly investigate the stability of positive equilibrium and the existence of Hopf bifurcation by choosing the time delay as a bifurcation parameter.

#### 1. Introduction

During the last few decades, there has been great interest in the construction and study of models for the population dynamics of predator-prey systems. The classical predator-prey model can be written in the generalized form (see [1]) where , denote the number or density of the prey and predator species at time , respectively, is the relative growth function of prey species in the absence of predator, is the relative mortality function of predator species in the absence of prey, is the functional response function of the predator to the prey which can be interpreted as the change in the number of prey species attached per unit time per predator as the prey number changes, and is the biomass conversion or consumption rate.

For system (1), the dynamic behavior is determined by the three crucial components: , , and . Chen et al. in [1] summarized some types of and which determine the species intrinsic population dynamics (primary production and mortality, resp.). For example, Functional response function also can be classified into many different types: Holling I–IV types [2, 3], Beddington-DeAngelis type [4, 5], Ivlev type [6], ratio-dependent type [7], and the modified forms of these types [8–18]. Over the years, the mathematical models with different , , and have been constructed by many researchers for describing different predator-prey systems; see [2–18].

Recently, one kind of predator-prey system with herd behavior is considered by some researchers [19–23]. In this kind of systems, such as herbivores on the large savanna and their large predators, the prey species join together in herds in order to provide a self-defense against predators, and the predator interacts with the prey along the outer corridor of the herd of prey [20, 21]. Therefore, it is more appropriate to model the response functions of prey in terms of the square root of the number of prey species for this class of systems. In [20], Braza considered a predator-prey system with square root functional response and compared the dynamics of this system with the dynamics of predator-prey system that use a typical Lotka-Volterra interaction term.

Moreover, there are some papers concerned with the effects of the form of mortality terms for the dynamics [23–26]. Brentnall et al. [24] pointed out that the mortality term of predator is usually described by “linear form” or “quadratic form” and the quadratic mortality is suited to intermediate predator (such as piscivore). Recently, Yuan et al. [23] studied a spatial predator-prey model with herd behavior and concluded that the Turing pattern is induced by quadratic mortality. Fulton et al. in [25] showed that, in the large and interlinked webs used in ecosystem models, model behavior is far more sensitive to the form of the grazing term than to that of the mortality terms that are close to the modeled food web.

On the other hand, for most of the natural ecosystems, every species does not respond instantaneously to changes in the environment or the interactions with other species within the community. Thus models with delay are much more realistic [27–31], and many researchers paid a lot of attention to the delayed predator-prey systems; see [32–40] and the references cited therein. For example, Xiao and Ruan [34] discussed a delayed predator-prey system with nonmonotonic functional response and showed that there is a B-T singularity for any time delay value and the system can exhibit Hopf bifurcation as the time delay passes through some critical values. Xu et al. in [36] proposed a delayed periodic L-V type predator-prey system with prey dispersal in two-patch environments and gave the sufficient conditions of the existence, uniqueness, and global stability of positive periodic solution.

Motivated by the above, in the present paper, we are intended to consider a predator-prey system in which the prey species exhibits herd behavior, the predator species is intermediate, and there exists a time delay in the convention of the biomass of prey to that of predator.

This paper is organized as follows. We first introduce our working system in the next section. In Section 3, by analyzing the characteristic equation, we discuss the local stability of a positive equilibrium and study the existence of Hopf bifurcations at the positive equilibrium. In Section 4, using the normal form theory and center manifold argument, the explicit formulae are derived to determine the direction of bifurcations and the properties of bifurcating periodic solutions. Several numerical simulations and a simple discussion are given in the last Section.

#### 2. The Working System

Braza [20] proposed the basic predator-prey system with logistic growth in the prey species, linear mortality in the predator species, and a square root functional response function. It is given by where , denote the number of the prey and predator species at time , respectively. The parameter is the growth rate of prey species, is its carrying capacity, is the death rate of predator species in the absence of prey, is the search efficiency of predator for prey, is the biomass conversion rate, and is predator's average handling time of prey.

Following References [23] and [24], we choose the quadratic mortality for predator species in system (3). If we further consider the delay fact, system (3) should be modified as the following form: where represents the quadratic mortality of predator species and is the time delay which means the growth rate of predator species to depend on the number of the prey species units of time earlier.

By the same way of [20, 23], we make some scaling and assume that the average handling time is zero. Then the working system is The initial conditions for system (5) take the form where , .

It is easy to verify that system (5) always has a trivial equilibrium and a boundary equilibrium for any feasible parameters. If we assume , then system (5) has a unique positive equilibrium , where

#### 3. Stability Analysis and Hopf Bifurcation

Due to biological interpretation of the system, we only consider the positive equilibrium. In this section, by choosing as the bifurcation parameter and analyzing the corresponding linearized system, we investigate the stability of the positive equilibrium and the effects of the time delay on the dynamics of system (5).

To study the local stability of the positive equilibrium , we first use the linear transformation , , for which system (5) can be written in the form where Thus, the positive equilibrium of system (5) is transformed into the zero equilibrium of system (8).

The linear equations corresponding to (8) are given by and its characteristic equation is To study the stability of the zero equilibrium and the Hopf bifurcation of system (8), we should analyze the distribution of the roots of characteristic equation (11). A simple calculation shows that ; thus is not a root of (11).

When there is no delay, that is, , characteristic equation (11) reduces to If , the two roots of (12) will have always negative real parts. Thus, we can obtain the following result.

Lemma 1. *Assume and ; then the two roots of characteristic equation (11) with have always negative real parts; that is, the zero equilibrium of system (8) with is locally asymptotically stable.*

Now for , let () be a root of (11); then satisfies the following equation: Separating the real and imaginary parts, we have which imply that Notice ; we can find that (15) has no positive roots if and has a positive root if , where Defining we have the following result.

Lemma 2. *Assume and ; then the following statements are true. *(1)*If , the roots of (11) have always negative real parts for all .*(2)*If , (11) has a pair of purely imaginary roots for .*

Next, under the condition of Lemma 2(2), we will investigate whether the transversality condition is satisfied. Denote by the root of (11) near satisfying , , . Differentiating the two sides of (11) with respect to , we get that is We can further get Then we have the following conclusion.

Lemma 3. *Assume , , and ; then the transversality conditions are satisfied; that is,
*

Summarizing the above three lemmas, we can obtain the following theorem on stability and Hopf bifurcation of system (5).

Theorem 4. *Assume and . For system (5), the following results are true.*(1)*If , the positive equilibrium is locally asymptotically stable for all .*(2)*If , then is locally asymptotically stable when and unstable when and system (5) undergoes a Hopf bifurcation at when , .*

#### 4. Direction and Stability of the Hopf Bifurcation

In the previous section, we have obtained the conditions under which a family of periodic solutions bifurcate from the positive equilibrium as crosses through the critical values . In this section, we will study the direction of Hopf bifurcation and the stability and period of the bifurcating periodic solutions by applying the normal form theory and the center manifold theorem from Hassard et al. [41]. Since the methods used are standard, we omit the detailed process and only give the main results. The readers can see [41] for more details on the derivation process.

Assume that system (5) undergoes Hopf bifurcations at for , and then are the corresponding purely imaginary roots of characteristic equation (11). Let We can choose such that is the eigenvector of operator corresponding to , is the eigenvector of operator corresponding to , and Thus, we can compute where Now, let which can be determined by the parameters in system (5). Defining we have the following result.

Theorem 5. *The direction of the Hopf bifurcation is determined by : if , then the Hopf bifurcation is supercritical (subcritical) and the bifurcating periodic solutions exist for ; determines the stability of the bifurcating periodic solutions: the bifurcating periodic solutions are stable (unstable) if ; and determines the period of the bifurcating periodic solutions: the period increases (decreases) if .*

#### 5. Numerical Simulations and Discussions

In this section, we present some numerical simulations to verify our theoretical results proved in previous sections by using MATLAB DDE solver. We simulate the system (5) with the initial value . In Figures 1–3 the initial point and the positive equilibrium are represented by green star and red star, respectively.

**(a)**

**(b)**

**(c)**

**(d)**

First, we choose parameters , . In this case, system (5) has only one positive equilibrium . Simple calculations show that , , which satisfy the conditions of Theorem 4(1). By Theorem 4(1), we expect that is stable for any delay value. We further choose four different delay values: , , , and . The computation simulations are depicted in Figure 1. Obviously, the computation simulations can support the result of Theorem 4(1).

Next, we choose parameters , . The positive equilibrium of system (5) is . It is easy to show that , , which satisfy the conditions of Theorem 4(2). From the formulae in the previous section, we can compute the values of , , , and as from which we conclude that is asymptotically stable for and the Hopf bifurcation of system (5) occurring at the critical value is supercritical and the bifurcating periodic solution exists when crosses to the right; also the bifurcating periodic solution is stable. By taking and , we can show the computation simulations in Figures 2 and 3, respectively. It is easy to see that the computation simulations support the results of Theorem 4(2) and Theorem 5.

In the present paper, we have considered a delayed predator-prey system in which the prey species exhibits herd behavior and the predator species with quadratic mortality. Our research shows that, for system (5), the positive equilibrium is always asymptotically stable under certain conditions (see Theorem 4(1)). In other words, under such conditions, the dynamics behavior of the corresponding ODE system is robust with respect to the time delay. But, under some other conditions (see Theorem 4(2)), will lose its stability and Hopf bifurcations occur when the delay passes through some critical values. In this case, the time delay has an important effect for the system dynamics.

Our results may enrich the dynamics in the predator-prey system and help us to better understand the interaction of predator with prey in a real ecosystem. Further studies are necessary to analyze the dynamics of more realistic but complex systems, such as delayed diffusive predator-prey system with herd behavior.

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgments

The authors wish to thank the handling editor and reviewers for their valuable comments and suggestions which lead to truly significant improvement of the paper. This research is supported by the National Natural Science Foundation of China (no. 11271260), Shanghai Leading Academic Discipline Project (no. XTKX2012), Innovation Program of Shanghai Municipal Education Commission (no. 13ZZ116), and the Innovation Fund Project For Graduate Student of Shanghai (no. JWCXSL1401).

#### References

- L. Chen, X. Song, and Z. Lu,
*Mathematical Ecological Models and Research Methods*, Sichuan Science and Technology Press, Chengdu, China, 2003. - C. S. Holling, “The functional response of predators to prey density and its role in mimicry and population regulation,”
*Memoirs of the Entomological Society of Canada*, vol. 45, supplement S45, pp. 5–60, 1965. View at: Publisher Site | Google Scholar - J. F. Andrews, “A mathematical model for the continuous culture of microorganisms utilizing inhibitory substrates,”
*Biotechnology and Bioengineering*, vol. 10, no. 6, pp. 707–723, 1968. View at: Publisher Site | Google Scholar - J. Beddington, “Mutual interference between parasites or predators and its effect on searching effciency,”
*Journal of Animal Ecology*, vol. 44, no. 1, pp. 331–340, 1975. View at: Publisher Site | Google Scholar - D. L. DeAngelis, R. A. Goldsten, and R. V. O'Neill, “A model for trophic interaction,”
*Ecology*, vol. 56, no. 4, pp. 881–892, 1975. View at: Publisher Site | Google Scholar - R. E. Kooij and A. Zegeling, “A predator-prey model with Ivlev's functional response,”
*Journal of Mathematical Analysis and Applications*, vol. 198, no. 2, pp. 473–489, 1996. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - R. Arditi and L. R. Ginzburg, “Coupling in predator-prey dynamics: ratio-dependence,”
*Journal of Theoretical Biology*, vol. 139, no. 3, pp. 311–326, 1989. View at: Google Scholar - Y. Kuang and E. Beretta, “Global qualitative analysis of a ratio-dependent predator-prey system,”
*Journal of Mathematical Biology*, vol. 36, no. 4, pp. 389–406, 1998. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - S. Ruan and D. Xiao, “Global analysis in a predator-prey system with nonmonotonic functional response,”
*SIAM Journal on Applied Mathematics*, vol. 61, no. 4, pp. 1445–1472, 2000. View at: Publisher Site | Google Scholar | MathSciNet - D. Xiao and H. Zhu, “Multiple focus and Hopf bifurcations in a predator-prey system with nonmonotonic functional response,”
*SIAM Journal on Applied Mathematics*, vol. 66, no. 3, pp. 802–819, 2006. View at: Publisher Site | Google Scholar | MathSciNet - H. Zhu, S. A. Campbell, and G. S. K. Wolkowicz, “Bifurcation analysis of a predator-prey system with nonmonotonic functional response,”
*SIAM Journal on Applied Mathematics*, vol. 63, no. 2, pp. 636–682, 2002. View at: Publisher Site | Google Scholar | MathSciNet - X.-A. Zhang, L. Chen, and A. U. Neumann, “The stage-structured predator-prey model and optimal harvesting policy,”
*Mathematical Biosciences*, vol. 168, no. 2, pp. 201–210, 2000. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - M. A. Aziz-Alaoui and M. D. Okiye, “Boundedness and global stability for a predator-prey model with modified Leslie-Gower and Holling-type II schemes,”
*Applied Mathematics Letters*, vol. 16, no. 7, pp. 1069–1075, 2003. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - A. F. Nindjin, M. A. Aziz-Alaoui, and M. Cadivel, “Analysis of a predator-prey model with modified Leslie-Gower and Holling-type II schemes with time delay,”
*Nonlinear Analysis: Real World Applications*, vol. 7, no. 5, pp. 1104–1118, 2006. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - L. Han, Z. Ma, and H. W. Hethcote, “Four predator prey models with infectious diseases,”
*Mathematical and Computer Modelling*, vol. 34, no. 7-8, pp. 849–858, 2001. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - R. Xu and Z. Ma, “Stability and Hopf bifurcation in a predator-prey model with stage structure for the predator,”
*Nonlinear Analysis: Real World Applications*, vol. 9, no. 4, pp. 1444–1460, 2008. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - M. Liu and K. Wang, “Global stability of stage-structured predator-prey models with Beddington-DeAngelis functional response,”
*Communications in Nonlinear Science and Numerical Simulation*, vol. 16, no. 9, pp. 3792–3797, 2011. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - C. Ji, D. Jiang, and X. Li, “Qualitative analysis of a stochastic ratio-dependent predatorprey system,”
*Journal of Computational and Applied Mathematics*, vol. 235, no. 5, pp. 1326–1341, 2011. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - V. Ajraldi, M. Pittavino, and E. Venturino, “Modeling herd behavior in population systems,”
*Nonlinear Analysis: Real World Applications*, vol. 12, no. 4, pp. 2319–2338, 2011. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - P. A. Braza, “Predator-prey dynamics with square root functional responses,”
*Nonlinear Analysis: Real World Applications*, vol. 13, no. 4, pp. 1837–1843, 2012. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - E. Venturino and S. Petrovskii, “Spatiotemporal behavior of a prey-predator system with a group defense for prey,”
*Ecological Complexity*, vol. 14, pp. 37–47, 2013. View at: Publisher Site | Google Scholar - S. Belvisi and E. Venturino, “An ecoepidemic model with diseased predators and prey group defense,”
*Simulation Modelling Practice and Theory*, vol. 34, pp. 144–155, 2013. View at: Publisher Site | Google Scholar - S. Yuan, C. Xu, and T. Zhang, “Spatial dynamics in a predator-prey model with herd behavior,”
*Chaos*, vol. 23, no. 3, Article ID 033102, 2013. View at: Publisher Site | Google Scholar - S. J. Brentnall, K. J. Richards, J. Brindley, and E. Murphy, “Plankton patchiness and its effect on larger-scale productivity,”
*Journal of Plankton Research*, vol. 25, no. 2, pp. 121–140, 2003. View at: Publisher Site | Google Scholar - E. A. Fulton, A. D. M. Smith, and C. R. Johnson, “Mortality and predation in ecosystem models: is it important how these are expressed?”
*Ecological Modelling*, vol. 169, no. 1, pp. 157–178, 2003. View at: Publisher Site | Google Scholar - M. Baurmann, T. Gross, and U. Feudel, “Instabilities in spatially extended predator-prey systems: spatio-temporal patterns in the neighborhood of Turing-Hopf bifurcations,”
*Journal of Theoretical Biology*, vol. 245, no. 2, pp. 220–229, 2007. View at: Publisher Site | Google Scholar | MathSciNet - N. MacDonald,
*Biological Delay Systems: Linear Stability Theory*, Cambridge University Press, Cambridge, UK, 1989. View at: MathSciNet - J. K. Hale and S. M. V. Lunel,
*Introduction to Functional-Differential Equations*, Springer, New York, NY, USA, 1993. View at: MathSciNet - Y. Kuang,
*Delay Differential Equations with Applications in Population Dynamics*, Academic Press, Boston, Mass, USA, 1993. View at: MathSciNet - J. Wu,
*Theory and Applications of Partial Functional-Differential Equations*, Springer, New York, NY, USA, 1996. View at: Publisher Site | MathSciNet - R. Xu, Q. Gan, and Z. Ma, “Stability and bifurcation analysis on a ratio-dependent predator-prey model with time delay,”
*Journal of Computational and Applied Mathematics*, vol. 230, no. 1, pp. 187–203, 2009. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - Y. Song and S. Yuan, “Bifurcation analysis in a predator-prey system with time delay,”
*Nonlinear Analysis: Real World Applications*, vol. 7, no. 2, pp. 265–284, 2006. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - Y. Song and J. Wei, “Local Hopf bifurcation and global periodic solutions in a delayed predator-prey system,”
*Journal of Mathematical Analysis and Applications*, vol. 301, no. 1, pp. 1–21, 2005. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - D. Xiao and S. Ruan, “Multiple bifurcations in a delayed predator-prey system with nonmonotonic functional response,”
*Journal of Differential Equations*, vol. 176, no. 2, pp. 494–510, 2001. View at: Publisher Site | Google Scholar | MathSciNet - S. Ruan, “Absolute stability, conditional stability and bifurcation in Kolmogorov-type predator-prey systems with discrete delays,”
*Quarterly of Applied Mathematics*, vol. 59, no. 1, pp. 159–173, 2001. View at: Google Scholar | Zentralblatt MATH | MathSciNet - R. Xu, M. A. J. Chaplain, and F. A. Davidson, “Periodic solutions for a delayed predator-prey model of prey dispersal in two-patch environments,”
*Nonlinear Analysis: Real World Applications*, vol. 5, no. 1, pp. 183–206, 2004. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - Z. Liu and R. Yuan, “Stability and bifurcation in a delayed predator-prey system with Beddington-DeAngelis functional response,”
*Journal of Mathematical Analysis and Applications*, vol. 296, no. 2, pp. 521–537, 2004. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - M. Fan and K. Wang, “Periodicity in a delayed ratio-dependent predator-prey system,”
*Journal of Mathematical Analysis and Applications*, vol. 262, no. 1, pp. 179–190, 2001. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - E. Beretta and Y. Kuang, “Global analyses in some delayed ratio-dependent predator-prey systems,”
*Nonlinear Analysis: Theory, Methods and Applications*, vol. 32, no. 3, pp. 381–408, 1998. View at: Publisher Site | Google Scholar | MathSciNet - T. Zhao, Y. Kuang, and H. L. Smith, “Global existence of periodic solutions in a class of delayed gause-type predator-prey systems,”
*Nonlinear Analysis: Theory, Methods and Applications*, vol. 28, no. 8, pp. 1373–1394, 1997. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - B. D. Hassard, N. D. Kazarinoff, and Y. H. Wan,
*Theory and Applications of Hopf Bifurcation*, Cambridge University Press, Cambridge, UK, 1981. View at: Zentralblatt MATH | MathSciNet

#### Copyright

Copyright © 2014 Chaoqun Xu and Sanling Yuan. 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.