Discrete Dynamics in Nature and Society

Volume 2014 (2014), Article ID 856475, 6 pages

http://dx.doi.org/10.1155/2014/856475

Research Article

Canard Limit Cycle of the Holling-Tanner Model

Chongwu Zheng,¹ Fengqin Zhang,¹ and Jianquan Li²

¹Department of Mathematics, Yuncheng University, Yuncheng, Shanxi 044000, China

²Science College, Air Force Engineering University, Xi’an 710051, China

Received 11 April 2014; Accepted 2 July 2014; Published 16 July 2014

Academic Editor: Delfim F. M. Torres

Copyright © 2014 Chongwu Zheng 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.

By using the singular perturbation theory on canard cycles, we investigate the canard phenomenon for the Holling-Tanner model with the intrinsic growth rate of the predator small enough. The obtained result shows that there may be at most one canard limit cycle, and the range of small parameters is estimated. The phenomenon of outbreak is explained.

1. Introduction

The predator-prey system

is referred to as the Holling-Tanner model [1]. In (1), the variables and denote the numbers of preys and predators at time , respectively. is the intrinsic growth rate of the prey, is the carrying capacity of the prey, and function , known as the functional response of Holling type II, is the rate at which predators consume the prey, where is the maximum number of a predator and is the half saturation value corresponding to the number of preys necessary to achieve one half of the maximum value. It is assumed that the predator grows logistically with intrinsic growth rate and carrying capacity taking on the role of a prey-dependent carrying capacity for the predator, where is a measure of the quality of the prey as food for the predator. From the biological meaning all the parameters of (1) are positive.

It is easy to see that the set is positively invariant for (1), and there is always a unique positive equilibrium. The local stability of the positive equilibrium was analyzed in [2, 3]. In [4] Hsu and Huang considered the question of global stability of the equilibrium and obtained conditions under which its local stability implies the global stability. In 2003 Braza [5] investigated the existence of at least one limit cycle of (1) in a Hopf bifurcation. Gasull et al. [6] showed that (1) can exhibit two limit cycles by constructing an example. The bifurcation diagram of limit cycles of (1) was described to show the regions of parameters where two limit cycles or a semistable limit cycle appears [7]. Applying Andronov-Hopf bifurcation Theorem Hsu and Hwang proved that the Hopf bifurcation is subcritical and thus (1) may have at least two limit cycles for some parameter ranges and explained the phenomenon of outbreaking in [8] obtained by numerical simulations [9].

In the references mentioned above, it was assumed that all the parameters are positive constants, but there is not much consideration for the system with small parameters. It is interesting and important to analyze the bifurcation of limit periodic set of the system in the case when is very small; that is, the intrinsic growth rate of the predator remains at a very low level. For this case when plotting the trajectories of system (1) with the intrinsic growth rate small enough by numerical simulations, it would take much longer time when the trajectory passes by the vertical -axis. This behaves like the usual relaxation oscillation phenomena observed in slow-fast systems. By using the singular perturbation theory [10] Li and Zhu [11] studied the canard limit cycles for predator-prey system with functional response of four Holling types and developed criteria for the existence and multiplicity of canard slow-fast cycle (i.e., limit periodic set) in general predator-prey systems.

Although some references have considered the existence of two limit cycles for (1) [6, 7, 9], it is still interesting to discuss the existence of canard slow-fast cycle, since the canard slow-fast cycle consists of a part of the fast orbit and a part of the slow manifold and is the limit of a family of periodic orbits as the value of a positive parameter tends to zero. In this paper we mainly consider the existence of canard slow-fast cycle of (1) when is very small. Our result shows that there may be at most one canard limit cycle, and the range of small parameters is estimated. In next section we first reduce the number of parameters to three ones by some transformation for (1) and introduce some concepts of canard slow-fast cycle. In Section 3 we provide the complete inferences for the main result of this paper. Finally we give discussion according to the results of numerical simulation.

2. Simplification of the Model and Canard Limit Cycle

For the case where is very small, let

where is positive and may be negative, and we will choose and small enough to make in later discussion; then system (1) becomes

For system (3), making the scaling

and then dropping their bars yield

where

Since and , letting gives the equivalent system of (5):

where

If we introduce slow time , then system (7) is transformed to

By lettng for (7) and (9), we obtain the layer equations

and the reduced equations

respectively.

In this paper we will always assume that . The slow manifold has a unique extreme point at

and it is a simple maximum point; see Figure 1(a). Besides, is globally concave for . The fast movement of the layer equations (10) is shown in Figure 1(b), where the orbits outside are parallel lines with different constant , and the curve consists of singular points; hence it is also called the critical curve. Note that the points on left to are normally hyperbolic repelling while the points on right to are normally hyperbolic attracting.

(a)

(b)

Figure 1: The limit periodic cycle and the slow manifold .

To study the slow movement on , from the reduced equations (11) we have

and hence along the slow movement is controlled by

In order to show the existence of canard slow-fast cycle, from Figure 1(b) it is necessary to illustrate that the right hand side of (14) should be negative near the point .

In the next section, we will make some transformations such that the function in system (7) becomes , and for the right hand side of (14) is negative in an interval of , containing ; then we may construct a slow-fast cycle (i.e., limit periodic set) by combining a part of the fast orbit , where , and a part of the slow manifold , shown in Figure 1(a). As we mentioned above, on the left and right sides of the point the points on the slow manifold are normally hyperbolic repelling and attracting, respectively, the point is called a canard point, and the slow-fast cycle is called a canard type, or simply called a canard slow-fast cycle; see [10, 12].

3. Existence and Uniqueness of Canard Slow-Fast Cycles

To study the local behavior of system (7) near , we make the change of coordinates

For near system (7) becomes

where

Notice that as (i.e., ); then for some and making small enough. For such and , making one more change of coordinates

equation (16) becomes

In addition, we choose and such that

and then (19) can be rewritten as

where

From (17) and (20) we have

Furthermore, implies that the change of coordinates (18) is available.

According to (3.12) and (3.13) of [10], from (24) we have

and then

Therefore, the canard cycle of (7) is asymptotically stable if it exists.

On the other hand, for (14) we have

for and . The fact that for and implies that system (14) along the slow manifold is well defined near and that the point is a canard one.

Further, by Theorems 3.3 of [10], implies that a -family of canard cycles bifurcates from the limit periodic set , for each , and for some small and . Moreover, satisfies

where is a constant and

The inequality (28) implies that, when there is a family of canard cycles bifurcating from the canard slow-fast cycle, in (7) depends on , and as .

Further, direct computation shows

According to the symmetry of the slow manifold about , for (i.e., ) we suppose that the interval of -coordinate of the slow-fast cycle is , where and ; see Figure 1(a). By straightforwardly calculating, we have

where

Note that

for and ; then implies that for ; that is, for . Hence, it follows form the results in [12] that (7) has at most one canard slow-fast cycle.

Summarizing the above inference, we obtain the following result.

Theorem 1. When , (7) has at most one canard slow-fast cycle, which is asymptotically stable for and satisfying (28).

4. Numerical Simulation and Discussion

We first plot the trajectories of (5) and the corresponding curves of and as , , , and (see Figure 2), in which Figures 2(b) and 2(c) illustrate the existence of the fast stages of and .

(a)

(b)

(c)

Figure 2: The trajectories of (5) with the fast stage and the slow stage and the corresponding curves of and , where , , , and .

Note that the vertical and horizontal isoclines of (5) divide the first quadrant of the -plane into four regions:

And from (2) and (6) it follows that , , and , where , , , and are the parameters in (5); then, corresponding to (1), the value of parameter is very small relative to other parameters as is very small.

In order to more clearly show the fast stages, we locally enlarge the two curves in Figures 2(b) and 2(c) into Figures 3(a) and 3(b), respectively. From Figures 3(a) and 3(b), we can know that the fast stages of occur in region , and those of occur in region . Biologically, and represent the numbers of preys and predators at time , respectively; then, when the intrinsic growth rate of the predator is very small and there exists the periodic solution for (1), according to signs of derivatives of and with respect to time , the number of the preys increases quickly from the number small enough as , and the number of the predators decreases quickly as . The quick increase of the prey from the small value may explain the phenomenon of outbreak in [8].

(a)

(b)

Figure 3: The curves of and obtained by locally enlarging the curves of and in Figure 2, respectively.

Conflict of Interests

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

Acknowledgment

This work is supported by the National Natural Science Fund of China (nos. 11371369 and 11371313).

References

J. T. Tanner, “The stability and the intrinsic growth rates of prey and predator populations,” Ecology, vol. 5, pp. 855–867, 1975.

J. D. Murray, Mathematical Biology, Springer, Berlin, Germany, 1989.

R. M. May, Stability and Complexity in Model Ecosystems, Princeton University Press, Princeton, NJ, USA, 1974.

S. B. Hsu and T. W. Huang, “Global stability for a class of predator-prey systems,” SIAM Journal on Applied Mathematics, vol. 55, no. 3, pp. 763–783, 1995.

P. A. Braza, “The bifurcation structure of the Holling-Tanner model for predator-prey interactions using two-timing,” SIAM Journal on Applied Mathematics, vol. 63, no. 3, pp. 889–904, 2003.

A. Gasull, R. E. Kooij, and J. Torregrosa, “Limit cycles in the Holling-Tanner model,” Publicacions Matemàtiques, vol. 41, pp. 149–167, 1997.

E. Sáez and E. González, “Dynamics of a predator-prey model,” SIAM Journal on Applied Mathematics, vol. 59, no. 5, pp. 1867–1878, 1999.

D. J. Wollkind, J. B. Collings, and J. A. Logan, “Metastability in a temperature-dependent model system for predator-prey mite outbreak interactions on fruit trees,” Bulletin of Mathematical Biology, vol. 50, no. 4, pp. 379–409, 1988.

S. Hsu and T. Hwang, “Hopf bifurcation analysis for a predator-prey system of Holling and Leslie type,” Taiwanese Journal of Mathematics, vol. 3, no. 1, pp. 35–53, 1999.

M. Krupa and P. Szmolyan, “Relaxation oscillation and canard explosion,” Journal of Differential Equations, vol. 174, no. 2, pp. 312–368, 2001.

C. Li and H. Zhu, “Canard cycles for predator-prey systems with Holling types of functional response,” Journal of Differential Equations, vol. 254, no. 2, pp. 879–910, 2013.

P. De Maesschalck and F. Dumortier, “Canard cycles in the presence of slow dynamics with singularities,” Proceedings of the Royal Society of Edinburgh A. Mathematics, vol. 138, no. 2, pp. 265–299, 2008.