• Views 356
• Citations 0
• ePub 15
• PDF 223
`Journal of Applied MathematicsVolume 2013 (2013), Article ID 678762, 8 pageshttp://dx.doi.org/10.1155/2013/678762`
Research Article

## Dynamics of a Stage Structured Pest Control Model in a Polluted Environment with Pulse Pollution Input

1Department of Mathematics, Anshan Normal University, Anshan, Liaoning 114007, China
2Department of Mathematics, Liaoning Normal University, Dalian, Liaoning 116029, China

Received 18 December 2012; Revised 3 August 2013; Accepted 23 August 2013

Copyright © 2013 Bing Liu 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 pollution model and impulsive delay differential equation, we formulate a pest control model with stage structure for natural enemy in a polluted environment by introducing a constant periodic pollutant input and killing pest at different fixed moments and investigate the dynamics of such a system. We assume only that the natural enemies are affected by pollution, and we choose the method to kill the pest without harming natural enemies. Sufficient conditions for global attractivity of the natural enemy-extinction periodic solution and permanence of the system are obtained. Numerical simulations are presented to confirm our theoretical results.

#### 1. Introduction

Nowadays, the problem of the world’s environmental pollution is serious, which has a frustrating effect on the ecosystem damage in the direct or indirect ways. Pollution leads to the living environmental change and gene mutation. It results in not only birth defects and deformities but also population variability, which decreases the number of the population in the nature and even makes them extinct. In order to assess the risk of the populations exposed to a polluted environment, in recent years, mathematical models concerning this topic have been studied extensively including continuous pollution input and impulsive pollution input [111].

As we all know, the predator-prey system can be used to model the process of controlling the pests by spraying pesticides, as well as relying on their natural enemies. However, in a polluted environment, some natural enemies are affected by pollution seriously and pests almost are not affected. For example, frogs are the natural enemies of beetles, locusts, and mole cricket, but some chemical plants discard waste products into rivers for their convenience, which cause severe water contamination, seriously injures frog’s reproductive system, and significantly decreases their fertility. Moreover, water pollution also causes large quantities of the fertilized eggs and tadpoles to die, resulting in the decrease of frogs. It is shown in a Sweden’s new study that male tadpoles can eventually grow into female frogs only in the environment similar to the nature but full of pollutants with estrogen. However, some male frogs have ovaries but no fallopian tubes, and they finally turn into lifelong infertile frogs, which are called “Yin and Yang frog”, and nearly one-third of the world’s frog species may be extinct because of the environmental pollution. People must control the period and quantity of emission of pollution to prevent natural enemy from extinction. In addition, too much pesticide spraying will reduce pests significantly; meanwhile, it also causes serious environmental pollution. Therefore, when controlling pests, we had better choose the method to kill the pests without polluting the environment and harming natural enemies at regular intervals.

The predator-prey models with stage structure for the predator were introduced or investigated by Hastings and Wang [1214]. Since the immature predator takes (which is called maturation time delay) units of time to mature, the death toll during the juvenile period should be considered, and time delays have important biological meanings in stage structured models. Recently, many models with time delay were extensively studied [1522].

According to the above biological background, in this paper, we suggest an impulsive predator-prey pollution model with stage structured for predator by introducing a constant periodic pollutant input and proportional killing pest at different fixed moments to model the process of pest control and polluted environment. Recently, there has been quite a lot of literatures on the applications of impulsive differential equations on population [1, 2, 8, 10, 11, 2031]. To our knowledge, there have been no results on this topic in the literature. The questions that arise here are as follows: how do we control the emission of pollution to prevent the extinction of natural enemies? Under what condition can the system be permanent? How can we control pests effectively?

The organization of this paper is as follows. In the next section, we formulate our model and give several lemmas which are useful for our main results. In Section 3 and Section 4, the sufficient conditions for the global attractivity of the “natural enemy-extinction” periodic solution and permanence of the system are obtained. We give a brief discussion of our results in Section 5. Numerical simulations are presented to illustrate our theoretical results.

#### 2. Model Formulation and Preliminaries

In this paper, we assume only that the natural enemies are affected by pollution and we choose the method to kill the pest without harming natural enemies. Then a pest control model with stage structure for natural enemy in a polluted environment by introducing a constant periodic pollutant input and killing pests at different fixed moment is formulated as follows: where ,  ,  ,  ,  . ,  , and represent the densities of prey (pest), immature, and mature predator (natural enemy) at time , respectively; , represent the concentration of pollution in the environment and organism at time , respectively; is intrinsic growth rate of the pests in the absence of natural enemies; is the pest capacity of environment; is the predation rate of natural enemy and represents the conversion rate at which ingested pest in excess of what is needed for maintenance is translated into natural enemy increase; is the saturation which represents that a certain amount of natural enemies can prey on a limited amount of pests, though the pests are numerous; and are the death rate of immature and mature natural enemies, respectively; in addition, we assume that juveniles suffer a mortality rate of (the through-stage death rate) and take units of time to mature. and are the dose-response parameters of species to the pollution in the immature and mature natural enemies, respectively; the exogenous quantity of impulsive input of pollutant into the environment at time is represented by ; () represents a proportional decrease of pest because of being harvested at time . The other parameters can be seen in [1].

The initial conditions for (1) are

Note that the variable does not appear in the first, third, forth, and fifth equations of system (1); hence, we only need to consider the subsystem of (1) as follows:

Lemma 1 (see [32]). Consider the following delay differential equation: where , , and are all positive constants and for ; we have(1)if , then ;(2)if , then .

Lemma 2 (see [1]). Consider the following subsystem of (1): Then system (5) has a unique positive -periodic solution , which is globally asymptotically stable, where and .

Lemma 3. There exists a constant such that , , , , and .

Proof. Define . Since , , in addition, , , thus for large enough. Define ; then, for , ; we have where . Consider the following impulse differential inequalities: We have Hence so is uniformly ultimately bounded. Therefore, by the definition of , system (1) is uniformly ultimately bounded. The proof is completed.

Lemma 4 (see [29]). If holds, system has a unique positive globally asymptotically stable periodic solution , .

Therefore, if holds, the system (3) has a natural enemy-extinction periodic solution . In this paper, we assume that always holds.

Remark 5 (see [1]). and are the concentration of pollution. To assure and , it is necessary that .

Remark 6 (see [1]). According to the biological significance, we assume .

#### 3. Global Attractivity of the “Natural Enemy-Extinction” Periodic Solution

In this section, we discuss under what condition the natural enemies will go extinct.

Denote where

Theorem 7. If , then the “mature natural enemy-extinction” periodic solution of system (3) is globally attractive.

Proof. Since , we have
By Lemma 2, for sufficiently small enough , there exists a positive constant such that holds for .
Note that Then we consider the following comparison system: According to Lemma 4, we know that is a unique globally asymptotically stable positive -periodic solution of system (17).
By using comparison theorem of impulsive differential equation, there exist a positive integer and a sufficiently small positive constant such that for all , holds. From (15), (19), and the second equation of (3), we obtain that for , holds.
Consider the following comparison equation:
By inequality (14), we have that holds; then, according to Lemma 1, we obtain that
By the comparison theorem of delay differential equation, we have .
Without loss of generality, we may assume that ( is sufficiently small positive constant such that ) for all ; by the first equation of system (3), we have
Consider the following comparison equation:
By Lemma 4, is a unique globally asymptotically stable positive -periodic solution of system (24). By using comparison theorem of impulsive differential equation, for above and large enough, we have
It follows from (19) and (26) that holds for large enough. Let ; we can get , so holds for large enough, which implies as . According to Lemma 2, , as . This completes the proof.

#### 4. Permanence

Definition 8. System (1) is said to be permanent if there are positive constants , , and a finite time such that for all solutions with initial conditions (2), , , , holds for all , .
Denote where

Theorem 9. If , there exists a positive constant such that for any solution of system (3).

Proof. Since , we can choose positive constants , , and such that and hold, where The second equation of system (3) can be written as Define Calculating the derivative of along the solution of (3), we have We claim that the inequality cannot hold for all . Otherwise, there is a positive constant such that for all . From the first equation of system (3), we have Consider the following comparison system: Then is a unique globally asymptotically stable positive -periodic solution of system (37). By using comparison theorem of impulsive differential equation, for , there exists a such that for , holds.
By Lemma 2, for , there exists a such that for , Let , and from (39) and (40), we have Let We will show that for all . Otherwise, there exists a nonnegative constant such that for , and . Thus from the second equation of (3), (31), and (41), we easily see that which is a contradiction. Hence we get that for all . Then we have which implies as . This is a contradiction to . Therefore, for any positive constant , the inequality cannot hold for all . If holds true for all large enough, then our aim is obtained; otherwise, is oscillatory about . Let In the following, we will show that . There exist two positive constants and such that When is large enough, the inequality holds true for . Since is continuous and bounded and not affected by impulses, we conclude that is uniformly continuous. Hence there exists a constant (, and is independent of the choice of ) such that for all .If , our aim is obtained.If , from the second equation of (3), we have that Then we have for . It is clear that for .
If , by the second equation of (3), then we have that for . The same arguments can be continued, and we can obtain for . Since the interval is arbitrarily chosen, we get that for t large enough. In view of our arguments above, the choice of is independent of the positive solution of (3). This completes the proof.

Theorem 10. If , the system (1) is permanent.

Proof. Suppose that is any positive solution of system (1) with initial conditions (2). By (39), we have for large enough. By Theorem 9, we have for large enough. From the second equation of system (1), we obtain So, Thus for large enough. By Lemma 2, we know for a sufficiently small positive , . Then from (15), Lemma 3, and Definition 8, we have that system (1) is permanent. The proof is completed.

#### 5. Discussion

In this paper, we discuss a pest control model with stage structure for natural enemy in a polluted environment by introducing a constant periodic pollutant input and killing pest at different fixed moments. From Theorems 7, 9, and 10, we can observe that the extinction and permanence of the population are very much dependent on , , and .

To verify the theoretical results obtained in this paper, in the following we will give some numerical simulations and take ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;  ;   (see Figure 1), and here we can compute , and from Theorem 10 we know the system (1) is permanent. If we decrease the period of pulsing () or increase the pollution input amount to (), and other parameters are the same with those in Figure 1, the natural enemy will be extinct (see Figures 2 and 3). If we increase the harvesting rate of pests to , and other parameters are the same with those in Figure 1, then , and the natural enemies will also be extinct (see Figure 4). Our results indicate that if impulsive period is short or or is too large, the natural enemy will go extinct, but we wish to protect natural enemy from extinction, so we should harvest the pests reasonably and control the period and quantity of emission of pollution into the environment efficiently. This offers us some reasonable suggestions for pest management.

Figure 1: Time series of system (1) with parameters ,  ,  ,  ,  ,  ,  ,  ,  ,  ,  ,  ,  ,  ,  ,  , and .
Figure 2: Time series of system (1) with parameters ,  ,  ,  ,  ,  ,  ,  ,  ,  ,  ,  ,  ,  ,  ,  , and .
Figure 3: Time series of system (1) with parameters ,  ,  ,  ,  ,  ,  ,  ,  ,  ,  ,  ,  ,  ,  ,  ,  and .
Figure 4: Time series of system (1) with parameters ,  ,  ,  ,  ,  ,  ,  ,  ,  ,  ,  ,  ,  ,  ,  ,  and .

#### Acknowledgments

This work is supported by National Natural Science Foundation of China (10971001, 11371030) and the Excellent Talents Support Project of Universities and Colleges in Liaoning.

#### References

1. B. Liu, L. S. Chen, and Y. J. Zhang, “The effects of impulsive toxicant input on a population in a polluted environment,” Journal of Biological Systems, vol. 11, no. 3, pp. 265–274, 2003.
2. B. Liu, L. Zhang, and Q. Zhang, “The effects of a single stage-structured population model with impulsive toxin input and time delays in a polluted environment,” Applicable Analysis, vol. 88, no. 8, pp. 1143–1155, 2009.
3. M. Liu and K. Wang, “Survival analysis of stochastic single-species population models in polluted environments,” Ecological Modelling, vol. 220, no. 9-10, pp. 1347–1357, 2009.
4. Z. Ma, Z. Luo, Z. Jin, and H. Zhao, “The threshold of survival for predator-prey Volterra system of three species in a polluted environment,” Systems Science and Mathematical Sciences, vol. 8, no. 4, pp. 373–382, 1995.
5. M. Zhien, Z. Wengang, and L. Zhixue, “The thresholds of survival for an $n$-dimensional food chain model in a polluted environment,” Journal of Mathematical Analysis and Applications, vol. 210, no. 2, pp. 440–458, 1997.
6. M. Debasis, “Persistence and global stability of a population in a polluted environment with delay,” Journal of Biological Systems, vol. 10, no. 3, pp. 225–232, 2002.
7. A. K. Pal and G. P. Samanta, “A single species population in a polluted environment,” International Journal of Biomathematics, vol. 3, no. 2, pp. 187–204, 2010.
8. F. Tao and B. Liu, “Dynamic behaviors of a single-species population model with birth pulses in a polluted environment,” The Rocky Mountain Journal of Mathematics, vol. 38, no. 5, pp. 1663–1684, 2008.
9. Y. Xiao and L. Chen, “How do the spatial structure and time delay affect the persistence of a polluted species,” Applicable Analysis, vol. 82, no. 3, pp. 253–267, 2003.
10. X. Yang, Z. Jin, and Y. Xue, “Weak average persistence and extinction of a predator-prey system in a polluted environment with impulsive toxicant input,” Chaos, Solitons & Fractals, vol. 31, no. 3, pp. 726–735, 2007.
11. Z. Zhao, L. Chen, and X. Song, “Extinction and permanence of chemostat model with pulsed input in a polluted environment,” Communications in Nonlinear Science and Numerical Simulation, vol. 14, no. 4, pp. 1737–1745, 2009.
12. A. Hastings, “Age-dependent predation is not a simple process. I. Continuous time models,” Theoretical Population Biology, vol. 23, no. 3, pp. 347–362, 1983.
13. A. Hastings, “Delays in recruitment at different trophic levels: effects on stability,” Journal of Mathematical Biology, vol. 21, no. 1, pp. 35–44, 1984.
14. W. Wang and L. Chen, “A predator-prey system with stage-structure for predator,” Computers & Mathematics with Applications, vol. 33, no. 8, pp. 83–91, 1997.
15. S. Gao, L. Chen, and Z. Teng, “Hopf bifurcation and global stability for a delayed predator-prey system with stage structure for predator,” Applied Mathematics and Computation, vol. 202, no. 2, pp. 721–729, 2008.
16. S. A. Gourley and Y. Kuang, “A stage structured predator-prey model and its dependence on maturation delay and death rate,” Journal of Mathematical Biology, vol. 49, no. 2, pp. 188–200, 2004.
17. H. Zhang and L. Chen, “A model for two species with stage structure and feedback control,” International Journal of Biomathematics, vol. 1, no. 3, pp. 267–286, 2008.
18. M. De la Sen, R. P. Agarwal, A. Ibeas, and S. Alonso-Quesada, “On a generalized time-varying SEIR epidemic model with mixed point and distributed time-varying delays and combined regular and impulsive vaccination controls,” Advances in Difference Equations, vol. 2010, Article ID 281612, 42 pages, 2010.
19. M. de la Sen, “On some structures of stabilizing control laws for linear and time-invariant systems with bounded point delays and unmeasurable states,” International Journal of Control, vol. 59, no. 2, pp. 529–541, 1994.
20. X. Song, M. Hao, and X. Meng, “A stage-structured predator-prey model with disturbing pulse and time delays,” Applied Mathematical Modelling, vol. 33, no. 1, pp. 211–223, 2009.
21. X. Meng and L. Chen, “Permanence and global stability in an impulsive Lotka-Volterra $n$-species competitive system with both discrete delays and continuous delays,” International Journal of Biomathematics, vol. 1, no. 2, pp. 179–196, 2008.
22. X. Meng, L. Chen, and H. Cheng, “Two profitless delays for the SEIRS epidemic disease model with nonlinear incidence and pulse vaccination,” Applied Mathematics and Computation, vol. 186, no. 1, pp. 516–529, 2007.
23. B. Liu, Q. Zhang, and Y. Gao, “The dynamics of pest control pollution model with age structure and time delay,” Applied Mathematics and Computation, vol. 216, no. 10, pp. 2814–2823, 2010.
24. Z. Liu and R. Tan, “Impulsive harvesting and stocking in a Monod-Haldane functional response predator-prey system,” Chaos, Solitons and Fractals, vol. 34, no. 2, pp. 454–464, 2007.
25. J. Hui and D. Zhu, “Dynamic complexities for prey-dependent consumption integrated pest management models with impulsive effects,” Chaos, Solitons and Fractals, vol. 29, no. 1, pp. 233–251, 2006.
26. Z. Jin, M. Haque, and Q. Liu, “Pulse vaccination in the periodic infection rate SIR epidemic model,” International Journal of Biomathematics, vol. 1, no. 4, pp. 409–432, 2008.
27. B. Liu, Z. D. Teng, and L. S. Chen, “The effect of impulsive spraying pesticide on stagestructured population models with birth pulse,” Journal of Biological Systems, vol. 13, no. 1, pp. 31–44, 2006.
28. B. Liu, Y. Zhang, and L. Chen, “Dynamic complexities in a Lotka-Volterra predator-prey model concerning impulsive control strategy,” International Journal of Bifurcation and Chaos in Applied Sciences and Engineering, vol. 15, no. 2, pp. 517–531, 2005.
29. X. Meng and L. Chen, “The dynamics of a new SIR epidemic model concerning pulse vaccination strategy,” Applied Mathematics and Computation, vol. 197, no. 2, pp. 582–597, 2008.
30. S. Tang and L. Chen, “Density-dependent birth rate, birth pulses and their population dynamic consequences,” Journal of Mathematical Biology, vol. 44, no. 2, pp. 185–199, 2002.
31. H. Zhang, P. Georgescu, and L. Chen, “An impulsive predator-prey system with Beddington-DeAngelis functional response and time delay,” International Journal of Biomathematics, vol. 1, no. 1, pp. 1–17, 2008.
32. Y. Kuang, Delay Differential Equations with Applications in Population Dynamics, vol. 191 of Mathematics in Science and Engineering, Academic Press, San Diego, Calif, USA, 1993.