- About this Journal ·
- Abstracting and Indexing ·
- Aims and Scope ·
- Annual Issues ·
- Article Processing Charges ·
- Articles in Press ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents

Discrete Dynamics in Nature and Society

Volume 2012 (2012), Article ID 150359, 12 pages

http://dx.doi.org/10.1155/2012/150359

## Seasonally Perturbed Prey-Predator Ecological System with the Beddington-DeAngelis Functional Response

^{1}School of Mathematics and Information Science, Wenzhou University, Wenzhou, Zhejiang 325035, China^{2}School of Life and Environmental Science, Wenzhou University, Wenzhou, Zhejiang 325027, China

Received 1 October 2011; Accepted 2 November 2011

Academic Editor: Vladimir Gontar

Copyright © 2012 Hengguo Yu and Min Zhao. 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

On the basis of the theories and methods of ecology and ordinary differential equation, a seasonally perturbed prey-predator system with the Beddington-DeAngelis functional response is studied analytically and numerically. Mathematical theoretical works have been pursuing the investigation of uniformly persistent, which depicts the threshold expression of some critical parameters. Numerical analysis indicates that the seasonality has a strong effect on the dynamical complexity and species biomass using bifurcation diagrams and Poincaré sections. The results show that the seasonality in three different parameters can give rise to rich and complex dynamical behaviors. In addition, the largest Lyapunov exponents are computed. This computation further confirms the existence of chaotic behavior and the accuracy of numerical simulation. All these results are expected to be of use in the study of the dynamic complexity of ecosystems.

#### 1. Introduction

Population communities are embedded in periodically varying environments. Therefore, it is appropriate to identify the functional role that seasons play in the behavior of population communities [1]. This seasonal variation can cause changes in the dynamics of an ecological system. The study of ecological systems which are subjected to seasonal perturbations is important for both theoretical and experimental ecologists. This study should also take into account the intrinsic nature of environmental and seasonal perturbations [2]. The basic problem is to understand the relationship between the magnitude of the seasonal variations and the complexity of the system. Numbers of studies have been performed on the interactions between the seasons and the internal biological rhythms of simple prey-predator ecosystems [3–6]. These studies have shown that these interactions can have spectacular consequences, such as multiplicity of attractors, catastrophes, and chaos [7].

In recent decades, it has been demonstrated that complex dynamics can appear in continuous-time models with three or more species [8–12], and specifically that nonlinear dynamics, including cycles, quasicycles, and chaos, can occur in such biological systems. In this context, continuous-time and discrete-time systems [13–21] have been discussed extensively by a number of researchers. The Beddington-DeAngelis functional response is introduced by Beddington [22] and DeAngelis et al. [23], independently. The main difference of this functional response from other functional responses is that it contains an extra term presenting mutual interference by predators [24]. Although a direct link between the predators and preys cannot be established unless quantitative methods are used, the precious works clearly show that the amount of three species are often related, and a change in one species can cause a change in the others, especially predator. Thus, we apply Beddington-DeAngelis functional response to describe their relationship with sufficient accuracy in this paper.

This paper considers a seasonally perturbed prey-predator system with the Beddington-DeAngelis functional response, which can be described by the following differential equations:

where , and are the densities of one prey and two predators at time , respectively, are the cropping rates, denote the efficiency with which preys are converted by new predators, is the carrying capacity of prey , and is the value of the limiting prey. In other words, is the theoretical carrying capacity under ideal conditions if there is no wastage in preys, which is impossible in reality. expresses the efficiency of nutrient utilization by a species and has a value between zero and one. If the ratio approaches unity, the efficiency is high; a lower ratio indicates that population increase is quickly restricted by the limiting prey. are saturation constants, and are scaling factors expressing the impact of predator interference, and is the relative superiority of predator . Without loss of generality, it can be assumed that ; accordingly, a predator with density is an inferior predator, while a predator with density is a superior predator because the ratio increases with time. , , and are the average values of , , and , respectively, and is the angular frequency of the fluctuations caused by seasonality. The parameters represent the degree of seasonality, and , and are the magnitudes of the perturbations in , , and , respectively [25]. Finally, the parameter , where , can be interpreted as the difference in phase angle between the seasonality of the above three parameters. Clearly, when, there is synchronous variation in the intrinsic growth rate and death rate, while when , the variation is antisynchronous. In this paper, only three values of phase angles will be considered: , and .

#### 2. Mathematical Analysis

*Definition 2.1. *System (1.1) is said to be uniformly persistent if, for any positive solution of system (1.1), there exist positive constants , , , , , , and such that , , and for .

Lemma 2.2 (see [26]). *If , , and , then, for any small constant , there exists a positive constant such that , for .*

Lemma 2.3 (see [27]). *If , , , and , then for .*

Now we can obtain the threshold expression of some critical parameters under the condition of all species persistence.

Theorem 2.4. *There exist positive constants , , and for any positive solution , , and of system (1.1) with all sufficiently large.*

*Proof. *Assume that is an arbitrary positive solution of system (1.1), then the first equation of system (1.1) can yield
From Lemma 2.2, we get that, for any small positive constant , there exists a constant such that for .

Defining that , then

Let , then
Thus is bound. Set such that . According to Lemma 2.3, . Therefore is ultimately bounded, and it follows that each positive solution of system (1.1) is uniformly ultimately bounded. Hence, there exist two positive , and such that and , for . This completes the proof.

Theorem 2.5. *If , and , then system (1.1) is uniformly persistent.*

*Proof. *From Theorem 2.4, it is obvious that there exist positive constants, , , and such that , , for .

On the other hand, we get from system (1.1) that
Then it follows from Lemma 2.2 that if , there is a positive constant such that for .

From the second equation of system (1.1), we can obtain that . It is obvious to know that if , there are two positive and such that for .

From the third equation of system (1.1), we can obtain that . It is obvious to know that if , there are two positive constants and such that for . Let then we have , , and for . This completes the proof.

By using the theoretical analysis, we obtain the threshold expression of some critical parameters under the condition of all species persistence, which in turn provides a theoretical basis for the numerical simulation.

#### 3. Simulation Analysis and Results

To study complex dynamics of the prey-predator system with seasonality, the solution of the system (1.1) with initial conditions is obtained numerically for the biologically feasible range of parametric values. In order to investigate the interplay among seasonal perturbation and some critical parameters, the diagrams as a function of and are obtained for three different cases: , . All other parameters are fixed as follows:, , , , , , , , , , , , , , and .

Let us investigate the critical parameter. Figure 1 shows the bifurcation diagram and the corresponding largest Lyapunov exponents as a function of in the range for , , and . It is evident from Figure 1 that the system produces complex dynamics for when , and , such as chaotic band, period-doubling bifurcation, and chaotic band with wide or narrow periodic windows. Nonetheless, it is interesting to notice that the preybiomass level increases with increase of , but it can be found that when , the prey biomass level is relatively high. Comparison of Figures 1(a), 1(b), and 1(c) suggests that seasonality may still give a region of periodic solutions for . The evidence for the cascade of periodic doubling leading to chaos can be seen in Figure 1 for , , and . However, Figure 1(a) shows chaotic bands with a narrow periodic window, and Figures 1(b) and 1(c) show chaotic bands with wide or narrow periodic windows. It is obvious that the seasonal disturbance not only enhances the prey biomass level, but also changes the complex dynamics when .

To investigate more carefully the effect of synchronous variation in the intrinsic growth rate and death rate, three strange attractors are plotted for and three different values of the synchronous variation parameter: , , . Figure 2 shows that when is 0 or , system (1.1) has a chaotic attractor at , but when, system (1.1) has a periodic attractor at . It can be observed from Figure 2 that seasonality in the intrinsic growth rate and death rate can change the chaotic behavior of the system at various values of the synchronous variation.

For further analysis of the effects of seasonality, three chaotic attractors as well as a Poincaré section (the Poincaré section is a classical technique for analyzing dynamic systems [27, 28]) have been obtained for and three different values of synchronous variation: , , . In this paper, the plane () was chosen, and values of were plotted against . The resulting attractors and Poincaré section are shown in Figure 3. Figure 3(a) shows a chaotic attractor at , and Figure 3(b) shows a Poincaré section for the chaotic attractor shown in (a). Three wide strips can be seen in Figure 3(b). Figure 3(c) shows a chaotic attractor at , and Figure 3(d) shows a Poincaré section for the chaotic attractor shown in (c). One wide strip can be seen in Figure 3(d). Figure 3(e) shows a chaotic attractor at , and Figure 3(f) shows a Poincaré section for the chaotic attractor shown in (e). Four wide strips can be seen in Figure 3(f). By comparison of the three Poincaré sections, it is apparent that the three chaotic attractors have different structures in the three different cases , , . However, they might be better described as stranger chaotic attractors and fractals. These results show that seasonal perturbation can change chaotic attractor structure. In a word, seasonal perturbation is sufficient to change the chaotic attractor trajectory in system (1.1).

Moreover, convincing evidence for deterministic chaos has come from several recent experiments [29–34]. The results of chaos studies have confirmed the importance of detecting and exploring chaos. Here, the largest Lyapunov exponents are considered because they have proved to offer the most useful diagnostics for a chaotic system [35–40]. The Lyapunov exponents take into account the average exponential rates of divergence or convergence of nearby orbits in phase space [29]. For a chaotic attractor, the values of the largest Lyapunov exponent must be positive; if they are always negative, chaos cannot be observed, and a stable state or a period attractor will instead appear. In the bifurcation diagram for resource shown in Figure 1(b), the corresponding largest Lyapunov exponent can be calculated for the system described by (1.1), which is shown in Figure 1(d). It should be stressed that the result is consistent with Figure 1(b), which shows the accuracy and effectiveness of numerical simulation. Moreover, using the simulation of the largest Lyapunov exponents, the existence of chaotic behavior in system (1.1) can be further confirmed.

The other critical parameter to be investigated is . Figure 4 depicts the local long-term dynamics maxima for , , and the range and . As increases from 0 to 0.14, Figure 4(a) shows that the system has rich dynamics, including periodic windows, period doubling, chaos, and a chaotic crisis. However, only chaos, chaotic crisis, and periodic windows are observed in Figures 4(b) and 4(c). Figure 4(b) shows that when is slightly increased beyond 0.14, the size of the chaotic attractor abruptly changes, constituting a type of attractor crisis. As continues to increase, a cascade of period-doubling bifurcation occurs, the period-2 attractor changes to a period-4 attractor, then a period-8 attractor and period-doubling bifurcation finally leads system (1.1) into chaos. However, only chaos and periodic windows are observed in Figures 4(a) and 4(c). In the range , the solution is chaotic with intermittent periodic windows. Similar behavior is observed in Figures 4(a) and 4(c). Hence, it is obvious that the seasonality has a strong effect on the dynamical behaviors of system (1.1) with the synchronous variation varying.

In order to study the effect of seasonality, three chaotic attractors as well as a Poincaré section have been obtained for and three different values of the synchronous variation parameter: , , and the plane . The attractors and Poincaré sections are shown in Figure 5. Figure 5(a) shows a chaotic attractor at , and Figure 5(b) shows a Poincaré section for the chaotic attractor shown in (a). Two wide point regions can be seen in Figure 5(b). Figure 5(c) shows a chaotic attractor at , and Figure 5(d) shows a Poincaré section for the chaotic attractor shown in (c). Two point regions which are wider than in the Poincaré section shown in (b) can be seen in Figure 5(d). Figure 5(e) shows a chaotic attractor at, and Figure 5(f) shows a Poincaré section for the chaotic attractor shown in (e). Two point regions which are narrower than in the other two Poincaré sections can be seen in Figure 5(f). By comparison of the three Poincaré sections, it is apparent that the three chaotic attractors have different structures in the three different cases , , . However, they are all fractals. These results show that the angular frequency of the fluctuations caused by seasonality can change chaotic attractor structure with the synchronous variation varying.

Based on the above analysis, it can be seen that the seasonal disturbance can enhance the prey biomass level for, in which result is agreed with some results in reality. Further, it is also interesting to point out that the seasonality in three different parameters can come into rich and complex dynamical behaviors, but these dynamical behaviors are different. Moreover, the use of mathematical model with seasonality is considered to investigate some biological problems, and the numerical simulation provides an approximation of the real biological system behaviors; hence, these results can promote the study of ecological dynamics.

#### 4. Conclusions and Remarks

In this paper, a prey-predator model with Beddington-DeAngelis functional response and seasonal perturbation has been studied analytically and numerically. Mathematical theoretical works have been pursuing the investigation of uniformly persistent, which depicts the threshold expression of some critical parameters and in turn provides a theoretical basis for the numerical simulation. Numerical analysis indicates that the seasonality has a strong effect on the dynamical complexity and species biomass using bifurcation diagrams and Poincaré sections. Bifurcation diagrams show that the seasonality in three different parameters can give rise to rich and complex dynamical behaviors, including periodic, period-doubling, period-halving, chaotic crises, and chaos. By comparing the Poincaré sections of various chaotic attractors, it can be determined that these chaotic attractors have different structures with the seasonality in three different parameters. Using numerical simulation of the largest Lyapunov exponents, the existence of chaotic behavior in system (1.1) and the accuracy and effectiveness of numerical simulation can be further confirmed. All these results are expected to be of use in the study of the dynamic complexity of ecosystems.

#### Acknowledgments

This work was supported by the National Natural Science Foundation of China (Grant no. 30970305 and Grant no. 31170338) and also by the Zhejiang Provincial Natural Science Foundation of China (Grant no. Y505365).

#### References

- S. Gakkhar and R. K. Naji, “Chaos in seasonally perturbed ratio-dependent prey-predator system,”
*Chaos, Solitons and Fractals*, vol. 15, no. 1, pp. 107–118, 2003. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - R. K. Upadhyay and S. R. K. Iyengar, “Effect of seasonality on the dynamics of 2 and 3 species prey-predator systems,”
*Nonlinear Analysis*, vol. 6, no. 3, pp. 509–530, 2005. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - S. Rinaldi, S. Muratori, and Y. A. Kuznetsov, “Multiple attractors, catastrophes and chaos in seasonally perturbed predator-prey communities,”
*Bulletin of Mathematical Biology*, vol. 55, no. 1, pp. 15–35, 1993. View at Publisher · View at Google Scholar · View at Scopus - H. Baek, “An impulsive two-prey one-predator system with seasonal effects,”
*Discrete Dynamics in Nature and Society*, vol. 2009, Article ID 793732, 19 pages, 2009. View at Google Scholar · View at Zentralblatt MATH - A. Gragnani and S. Rinaldi, “A universal bifurcation diagram for seasonally perturbed predator-prey models,”
*Bulletin of Mathematical Biology*, vol. 57, no. 5, pp. 701–712, 1995. View at Google Scholar · View at Scopus - H. G. Yu, S. M. Zhong, R. P. Agarwal, and S. K. Sen, “Effect of seasonality on the dynamical behavior of an ecological system with impulsive control strategy,”
*Journal of the Franklin Institute*, vol. 348, no. 4, pp. 652–670, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - T. S. Parker and L. O. Chua,
*Practical Numerical Algorithms for Chaotic System*, Springer, New York, NY, USA, 1989. - W. M. Schaffer, “Order and chaos in ecological systems,”
*Ecology*, vol. 66, no. 1, pp. 93–106, 1985. View at Google Scholar · View at Scopus - R. K. Upadhyay and V. Rai, “Crisis-limited chaotic dynamics in ecological systems,”
*Chaos, Solitons and Fractals*, vol. 12, no. 2, pp. 205–218, 2001. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - S. J. Lv and M. Zhao, “The dynamic complexity of a three species food chain model,”
*Chaos, Solitons and Fractals*, vol. 37, no. 5, pp. 1469–1480, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - S. Gakkhar and R. K. Naji, “Order and chaos in a food web consisting of a predator and two independent preys,”
*Communications in Nonlinear Science and Numerical Simulation*, vol. 10, no. 2, pp. 105–120, 2005. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - F. Y. Wang, C. P. Hao, and L. S. Chen, “Bifurcation and chaos in a Monod-Haldene type food chain chemostat with pulsed input and washout,”
*Chaos, Solitons and Fractals*, vol. 32, no. 1, pp. 181–194, 2007. View at Publisher · View at Google Scholar - Y. K. Xue and X. F. Duan, “The dynamic complexity of a Holling type-IV predator-prey system with stage structure and double delays,”
*Discrete Dynamics in Nature and Society*, vol. 2011, Article ID 509871, 19 pages, 2011. View at Google Scholar · View at Zentralblatt MATH - M. Zhao and S. J. Lv, “Chaos in a three-species food chain model with a Beddington-DeAngelis functional response,”
*Chaos, Solitons and Fractals*, vol. 40, no. 5, pp. 2305–2316, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - S. Gakkhar and R. K. Naji, “Seasonally perturbed prey-predator system with predator-dependent functional response,”
*Chaos, Solitons and Fractals*, vol. 18, no. 5, pp. 1075–1083, 2003. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - T. W. Hwang, “Uniqueness of limit cycles of the predator-prey system with Beddington-DeAngelis functional response,”
*Journal of Mathematical Analysis and Applications*, vol. 290, no. 1, pp. 113–122, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - N. Zhang, F. D. Chen, Q. Q. Su, and T. Wu, “Dynamic behaviors of a harvesting Leslie-Gower predator-prey model,”
*Discrete Dynamics in Nature and Society*, vol. 2011, Article ID 473949, 14 pages, 2011. View at Google Scholar · View at Zentralblatt MATH - V. Křivan and J. Eisner, “The effect of the Holling type II functional response on apparent competition,”
*Theoretical Population Biology*, vol. 70, no. 4, pp. 421–430, 2006. View at Publisher · View at Google Scholar · View at Scopus - X. Li and W. Yang, “Permanence of a discrete predator-prey systems with Beddington-deAngelis functional response and feedback controls,”
*Discrete Dynamics in Nature and Society*, vol. 2008, Article ID 149267, 8 pages, 2008. View at Google Scholar - R. X. Wu and L. Li, “Permanence and global attractivity of discrete predator-prey system with hassell-varley type functional response,”
*Discrete Dynamics in Nature and Society*, vol. 2009, Article ID 323065, 17 pages, 2009. View at Google Scholar · View at Zentralblatt MATH - L. J. Chen, J. Y. Xu, and Z. Li, “Permanence and global attractivity of a delayed discrete predator-prey system with general holling-type functional response and feedback controls,”
*Discrete Dynamics in Nature and Society*, vol. 2008, Article ID 629620, 17 pages, 2008. View at Google Scholar · View at Zentralblatt MATH - J. R. Beddington, “Mutual interference between parasites or predator and its effect on searching efficiency,”
*Journal of Animal Ecology*, vol. 44, no. 1, pp. 331–340, 1975. View at Google Scholar - D. L. DeAngelis, R. A. Goldstein, and R. V. Neill, “A model for trophic interaction,”
*Ecology*, vol. 56, no. 4, pp. 881–892, 1975. View at Google Scholar - H. Baek, “Species extinction and permanence of an impulsively controlled two-prey one-predator system with seasonal effects,”
*BioSystems*, vol. 98, no. 1, pp. 7–18, 2009. View at Publisher · View at Google Scholar · View at Scopus - S. Gakkhar and R. K. Naji, “Seasonally perturbed prey-predator system with predator-dependent functional response,”
*Chaos, Solitons and Fractals*, vol. 18, no. 5, pp. 1075–1083, 2003. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - K. Wang, “Permanence and global asymptotical stability of a predatorprey model with mutual interference,”
*Nonlinear Analysis: Real World Applications*, vol. 12, no. 2, pp. 1062–1071, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - F. D. Chen, “On a nonlinear nonautonomous predator-prey model with diffusion and distributed delay,”
*Journal of Computational and Applied Mathematics*, vol. 180, no. 1, pp. 33–49, 2005. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - R. Barrio, “Sensitivity tools vs. Poincaré sections,”
*Chaos, Solitons and Fractals*, vol. 25, no. 3, pp. 711–726, 2005. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - C. Masoller, A. C. S. Schifino, and L. Romanelli, “Characterization of strange attractors of lorenz model of general circulation of the atmosphere,”
*Chaos, Solitons and Fractals*, vol. 6, pp. 357–366, 1995. View at Google Scholar · View at Scopus - M. Zhao and L. M. Zhang, “Permanence and chaos in a host-parasitoid model with prolonged diapause for the host,”
*Communications in Nonlinear Science and Numerical Simulation*, vol. 14, no. 12, pp. 4197–4203, 2009. View at Publisher · View at Google Scholar · View at Scopus - M. Hagmüller and G. Kubin, “Poincaré pitch marks,”
*Speech Communication*, vol. 48, no. 12, pp. 1650–1665, 2006. View at Publisher · View at Google Scholar · View at Scopus - H. G. Yu, M. Zhao, S. J. Lv, and L. L. Zhu, “Dynamic complexities in a parasitoid-host-parasitoid ecological model,”
*Chaos, Solitons and Fractals*, vol. 39, no. 1, pp. 39–48, 2009. View at Publisher · View at Google Scholar · View at Scopus - J. G. Sportt,
*Chaos and Time-Series Analysis*, Oxford University Press, New York, NY, USA, 2003. - M. T. Rosenstein, J. J. Collins, and C. J. de Luca, “A practical method for calculating largest Lyapunov exponents from small data sets,”
*Physica D*, vol. 65, no. 1-2, pp. 117–134, 1993. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - S. J. Lv and M. Zhao, “The dynamic complexity of a host-parasitoid model with a lower bound for the host,”
*Chaos, Solitons and Fractals*, vol. 36, no. 4, pp. 911–919, 2008. View at Publisher · View at Google Scholar - Q. Jia, “Hyperchaos generated from the Lorenz chaotic system and its control,”
*Physics Letters, Section A*, vol. 366, no. 3, pp. 217–222, 2007. View at Publisher · View at Google Scholar · View at Scopus - M. Zhao, L. M. Zhang, and J. Zhu, “Dynamics of a host-parasitoid model with prolonged diapause for parasitoid,”
*Communications in Nonlinear Science and Numerical Simulation*, vol. 16, no. 1, pp. 455–462, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH - M. Zhao, H. G. Yu, and J. Zhu, “Effects of a population floor on the persistence of chaos in a mutual interference host-parasitoid model,”
*Chaos, Solitons and Fractals*, vol. 42, no. 2, pp. 1245–1250, 2009. View at Publisher · View at Google Scholar · View at Scopus - F. Grond, H. H. Diebner, S. Sahle, A. Mathias, S. Fischer, and O. E. Rossler, “A robust, locally interpretable algorithm for Lyapunov exponents,”
*Chaos, Solitons and Fractals*, vol. 16, no. 5, pp. 841–852, 2003. View at Publisher · View at Google Scholar - L. L. Zhu and M. Zhao, “Dynamic complexity of a host-parasitoid ecological model with the Hassell growth function for the host,”
*Chaos, Solitons and Fractals*, vol. 39, no. 1, pp. 1259–1269, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH