Research Article | Open Access

# Dynamics of a Stage Structure Pest Control Model with Impulsive Effects at Different Fixed Time

**Academic Editor:**Leonid Berezansky

#### Abstract

Many existing pest control models, which control pests by releasing natural enemies, neglect the effect that natural enemies may get killed. From this point of view, we formulate a pest control model with stage structure for the pest with constant maturation time delay (through-stage time delay) and periodic releasing natural enemies and natural enemies killed at different fixed time and perform a systematic mathematical and ecological study. By using the comparison theorem and analysis method, we obtain the conditions for the global attractivity of the pest-eradication periodic solution and permanence of the system. We also present a pest management strategy in which the pest population is kept under the economic threshold level (ETL) when the pest population is uniformly permanent. We show that maturation time delay, impulsive releasing, and killing natural enemies can bring great effects on the dynamics of the system. Numerical simulations confirm our theoretical results.

#### 1. Introduction

The warfare between man and pests has been on for thousands of years, and pest outbreaks often cause serious ecological and economic problems. With the development of society, human beings have come up with numerous methods to control pests, for instance, biological, cultural, physical, mechanical, and chemical tools in a way that minimize economics. Among those methods, releasing natural enemies is an effective and basic one, and it cannot pollute the environment. Morever, it is relatively simple to implement.

The pest control models by releasing natural enemies have been studied by many researchers [1–3]. All of them have invariably assumed that the releasing enemies can all have an effect on the models, whereas it is often the case that the natural enemies may get killed. For example, releasing frogs may prevent locusts. But in many places in China, especially in some southern areas, it is customary to eat frogs, and in markets the price of a frog is 15 yuan; so the number of the released frogs is reduced in a period as a result of their getting killed by people. According to the above biological background, we formulate a pest control model with periodic impulsive releasing natural enemies and natural enemies killed at different fixed time. We assume that the pest has a stage structure with constant maturation time delay (through-stage time delay). Many stage structure models with time delay were extensively studied (see [4–9]). In practice, from the principle of ecosystem balance, we need only to control the pest population under the economic threshold level (ETL) and not to eradicate natural enemy totally and hopes pest population and natural enemy population can coexist when the pests do not bring about immense economic losses. Then the questions that arise here are the following: how do the impulse period and time delay affect the extinction of the pest and permanence of the system? How many natural enemies should we release to control pests? How to keep the pests under ETL?

For these purposes, in Section 2 we suggest a delay impulsive differential equation to model the process of impulsive releasing natural enemies and natural enemies killed at different fixed time and introduce lemmas which will be used in this paper. Impulsive differential equations are found in almost every domain of applied science and have been studied in many investigations [10–17], and it can describe population dynamic models, since many life phenomena and human exploitation are almost impulsive in the natural world, and impulsive delay differential equations are almost analyzed in theory (see [18–20]). Time delay and impulse are introduced into pest control models with stage structure, which greatly enrich biologic background, but the system becomes nonautonomous and quite complicated, which causes great difficulties for us to study the model. In Sections 3 and 4, the conditions for the global attractivity of the pest-eradication periodic solution and permanence of the system are obtained. We give a brief conclusion of our results in the last section. Numerical simulations are presented to illustrate our theoretical results.

#### 2. Model Formulation and Auxiliary Lemmas

A model of single pest population growth incorporating stage structure as a reasonable generalization of the logistic model was derived as follows in [21]: where and represent the immature and mature pest population densities, respectively. This model assumes an average age to maturity which appears as a constant time delay reflecting a delayed birth of immature and a reduced survival of immature to their maturity. is the coefficient of birth rate of the mature population in this environment; , called maturation time delay, is the time to maturity; is the coefficient of death rate of the immature population; is the mature pest death and overcrowding rate; the term represents the immature pests which are born at the time (i.e., ) and survive at the time (with the immature pest death rate ) and therefore, represents the transformation of immatures to matures.

The basic model that we consider in this paper is the one based on the ideas that the natural enemies may be killed when they are released to control pest population; then we formulate the following two-stage pest control model with stage structure for pest and pulse releasing of natural enemies killed at different fixed time: where , , and represent the density of the immature pest, mature pest, and natural enemy population at time , respectively; , and ; is the natural enemy death's rate; is the natural enemy releasing amount at every impulsive period ; is the killed rate of natural enemies at every impulsive period . The meanings of other parameters are the same as those of Model (2.1).

The first equation in system (2.2) may be rewritten as follows: and ; that is can be linear expression by , and does not appear in the second and third equations in system (2.2), therefore in the rest of this paper, we will study the subsystem of (2.2) as follows: The initial conditions for (2.4) are From the biological point of view, we only consider system (2.4) in the biological meaning region:

Lemma 2.1 (see [22]). * Consider the following equation:
**
where , , , are all positive constants, , for ; one has the following.*(1)*If , then .*(2)*If , then *

Lemma 2.2 .. *For each solution of system (2.2) with being large enough, one has , where and .*

*Proof. *Let , . If and , we have
If , ; if , . We have
So is uniformly ultimately bounded. Hence, by the definition of , we obtain that each positive solution of system (2.2) is uniformly ultimately bounded. The proof is complete.

Lemma 2.3 (see [23]). *Consider the following system:
**
where are all positive constants; then, system (2.9) has a globally asymptotic stable positive periodic solution:
*

Lemma 2.4 (see [24]). *Let and . Assume that
**
where is continuous in and for exits, and is nondecreasing. Let be the maximal solution of the scalar impulsive differential equation:
**
exiting on . Then implies that , where is any solution of (2.11).*

#### 3. Global Attractivity of the “Pest-Eradication” Periodic Solution

Define

Theorem 3.1 .. *If holds, the population-extinction periodic solution of (2.4) is globally attractive.*

*Proof. *From the second equation in system (2.4), we have ; then consider the following system:
According to Lemma 2.3, the globally asymptotic stable positive periodic solution of system (3.1) is
According to Lemma 2.4, for any , there exits a such that
From the first equation of system (2.4), we can get
Then consider the following comparison equation:
Since , according to Lemma 2.1, we obtain ; then by the comparison theorem in differential equations, we get Without loss of generality, we assume that
Consider the following system:
The globally asymptotic stable positive periodic solution of system (3.7) is
By (3.3) and Lemma 2.4, for any , there exits a , when
Let , we can get
for being large enough, which implies This completes the proof.

#### 4. Permanence

Define

*Definition 4.1 .. **System (2.4) is said to be permanent if there are constants ** and a finite time ** such that for every positive solution ** with initial conditions (2.5) satisfies ** for all **. Here ** may depend on the initial condition (2.5).*

Theorem 4.2 .. *If holds, there exists a positive constant such that each positive solution of (2.4) satisfies with being large enough.*

*Proof. *The first equation of system (2.4) may be rewritten as follows:

Define
Calculating the derivative of along the solution to (2.4) gives
Since , we can choose . Let be small enough such that and
where For any positive constant , we claim that the inequality cannot hold for all . Otherwise, there is a positive constant such that for all . From the second and the fourth equations in system (2.4), we have
According to Lemma 2.4, there exists such that for ,
From (4.3) and (4.6) we have
Let we show that for all Otherwise, there exists a nonnegative constant , such that for , and

Thus, from the first equation of system (2.4) and (4.5), we easily see that
This is a contradiction. So we obtain that for all From (4.7), 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. Othervise, is oscillatory about . Let
In the following, we will show that for being large enough. There exist two positive constants such that
Since is continuous and bounded and is not effected by impulses, we conclude that is uniformly continuous. Then there exists a constant , such that for all

If , our aim is obtained.

If from the first equation of (2.4) we have that
Then we have
It is clear that for

If , by the first equation of (2.4), then we have that for Thus, proceeding exactly as the proof for above claim, we can obtain the inequality cannot hold for all so the same arguments can be continued for for and we can get for . Since the interval is arbitrarily chosen (we only need to be large), we get that for being large enough. In view of our arguments above, the choice of is independent of the positive solution of (2.4) which satisfies that for being sufficiently large . This completes the proof.

Theorem 4.3 .. *If , then system (2.4) is permanent.*

*Proof. *Suppose that is any solution of system (2.4) with initial condition (2.5). Let , by (3.3) we know , and according to Theorem 4.2, there exist positive constants , such that . Set
Then, is a bounded compact region which has positive distance from coordinate axes. By Theorem 4.2, one obtains that every solution to system (2.4) with the initial condition (2.5) eventually enters and remains in the region . The proof is completed.

The immature pest is hardly any harmful for the crop; so we just consider the effect of the mature pest. From the standpoint of ecological balance and saving resource, we only need to maintain the mature pest population under the economic threshold level (ETL) and not to eradicate the pest totally; then we have the following theorem.

Theorem 4.4 .. *Under the condition of Theorem 4.3, if
**
then the mature pest population is eventually under the economic threshold level *

*Proof. *From (3.3) we know for any that there exits a such that
From the first equation of system (2.4), we can get
Then consider the following comparison equation:
If holds, we can get ;then according to Lemma(2.1), we obtain
According to Lemma 2.4, we have as . The proof is completed.

#### 5. Discussion

In this paper, we discuss a pest control model with stage structure for the pest with constant maturation time delay (through-stage time delay) and periodic releasing natural enemies and natural enemies killed at different fixed time. From Theorems 3.1, 4.2, and 4.3, we can observe that the extinction and permanence of the population are very much dependent on . If is too large, from the condition of Theorem 3.1 we know that the population will be extinct. Although in theory pests can be eliminated completely, in fact it is hard to implement. In the natural world, the immature pest does not have any effect on the crop; so we only consider the mature pest in this paper. Even though pest is extinct, the food chain between pest and natural enemy is broken, that would be another disaster. Therefore we only need to control the mature pest population under (ETL) and not to eradicate natural enemy totally and hope that pest population and natural enemy population can coexist when the pests do not bring about immense economic losses.

To verify the theoretical results obtained in this paper, in the following we will give some numerical simulations and take ; ; ; ; ; ; ; ; ; (see Figure 1); here we can compute ; from Theorem 4.3 we know that system (2.4) is permanent. If we increase the natural enemy input amount to () or increase the maturation time delay (), other parameters are the same with those in Figure 1, and the pest will be extinct (see Figure 2 and Figure 3). If we decrease the killed rate of natural enemies to , then , the pest also will be extinct (see Figure 4). It shows that if the maturation time delay is too long, or the natural enemy releasing amount is too large, or fewer natural enemies are killed, the permanence of the system disappears and the pest population dies out. This implies that pulse releasing natural enemies, the maturation time delay and natural enemies killed bring great effects on the dynamics behaviors of the model. Suppose that the economic threshold level, , is 0.3, then Figure 5 shows that if we choose an appropriate pulse releasing natural enemy period (with ), we can control the pest population under the ETL and not eradicate natural enemy totally. This gives us some reasonable suggestions for pest management.

**(a)**

**(b)**

**(c)**

**(a)**

**(b)**

**(c)**

**(a)**

**(b)**

**(c)**

**(a)**

**(b)**

**(c)**

**(a)**

**(b)**

**(c)**

#### Acknowledgments

This work is supported by the National Natural Science Foundation of China (10671001 and 10971001) and the Science and Research Project Foundation of Liaoning Province Education Department.

#### References

- B. Liu, Y. Zhang, and L. Chen, “The dynamical behaviors of a Lotka-Volterra predator-prey model concerning integrated pest management,”
*Nonlinear Analysis: Real World Applications*, vol. 6, no. 2, pp. 227–243, 2005. View at: Publisher Site | Google Scholar | MathSciNet - H. Zhang, L. Chen, and J. J. Nieto, “A delayed epidemic model with stage-structure and pulses for pest management strategy,”
*Nonlinear Analysis: Real World Applications*, vol. 9, no. 4, pp. 1714–1726, 2008. View at: Publisher Site | Google Scholar | Zentralblatt MATH - 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. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - S. Tang and L. Chen, “Multiple attractors in stage-structured population models with birth pulses,”
*Bulletin of Mathematical Biology*, vol. 65, no. 3, pp. 479–495, 2003. View at: Publisher Site | Google Scholar - X.-Y. Song and L.-S. Chen, “Optimal harvesting and stability for a predator-prey system with stage structure,”
*Acta Mathematicae Applicatae Sinica*, vol. 18, no. 3, pp. 423–430, 2002. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - X. Song and H. Guo, “Global stability of a stage-structured predator-prey system,”
*International Journal of Biomathematics*, vol. 1, no. 3, pp. 313–326, 2008. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - G. H. Zhu, X. Z. Meng, and L. S. Chen, “The dynamics of a mutual interference age structured predator-prey model with time delay and impulsive perturbations on predators,”
*Applied Mathematics and Computation*, vol. 216, no. 1, pp. 308–316, 2010. View at: Google Scholar - S. Gao, L. Chen, and Z. Teng, “Pulse vaccination of an SEIR epidemic model with time delay,”
*Nonlinear Analysis: Real World Applications*, vol. 9, no. 2, pp. 599–607, 2008. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - X. Meng and L. Chen, “A stage-structured SI eco-epidemiological model with time delay and impulsive controlling,”
*Journal of Systems Science and Complexity*, vol. 21, no. 3, pp. 427–440, 2008. View at: Publisher Site | Google Scholar | Zentralblatt MATH - S. Sun and L. Chen, “Permanence and complexity of the eco-epidemiological model with impulsive perturbation,”
*International Journal of Biomathematics*, vol. 1, no. 2, pp. 121–132, 2008. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - 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. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - J. Jiao and L. Chen, “Global attractivity of a stage-structure variable coefficients predator-prey system with time delay and impulsive perturbations on predators,”
*International Journal of Biomathematics*, vol. 1, no. 2, pp. 197–208, 2008. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - 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. View at: Publisher Site | Google Scholar | MathSciNet - S. Gao and L. Chen, “The effect of seasonal harvesting on a stage-structured discrete model with birth pulses,”
*International Journal of Bifurcation and Chaos*, vol. 16, no. 9, pp. 2575–2586, 2006. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - J. Hui and L.-S. Chen, “Dynamic complexities in ratio-dependent predator-prey ecosystem models with birth pulse and pesticide pulse,”
*International Journal of Bifurcation and Chaos*, vol. 14, no. 8, pp. 2893–2903, 2004. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - K. Liu and L. Chen, “On a periodic time-dependent model of population dynamics with stage structure and impulsive effects,”
*Discrete Dynamics in Nature and Society*, vol. 2008, Article ID 389727, 15 pages, 2008. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - 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. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - J. Yan, “Stability for impulsive delay differential equations,”
*Nonlinear Analysis: Theory, Methods & Applications*, vol. 63, no. 1, pp. 66–80, 2005. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - L. Berezansky and E. Braverman, “Linearized oscillation theory for a nonlinear delay impulsive equation,”
*Journal of Computational and Applied Mathematics*, vol. 161, no. 2, pp. 477–495, 2003. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - X. Liu and G. Ballinger, “Boundedness for impulsive delay differential equations and applications to population growth models,”
*Nonlinear Analysis: Theory, Methods & Applications*, vol. 53, no. 7-8, pp. 1041–1062, 2003. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - W. G. Aiello and H. I. Freedman, “A time-delay model of single-species growth with stage structure,”
*Mathematical Biosciences*, vol. 101, no. 2, pp. 139–153, 1990. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - X. Song and L. Chen, “Optimal harvesting and stability for a two-species competitive system with stage structure,”
*Mathematical Biosciences*, vol. 170, no. 2, pp. 173–186, 2001. View at: Publisher Site | Google Scholar | MathSciNet - B. Liu, Z. Teng, and L. Chen, “Analysis of a predator-prey model with Holling II functional response concerning impulsive control strategy,”
*Journal of Computational and Applied Mathematics*, vol. 193, no. 1, pp. 347–362, 2006. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - D. Bainov and P. Simeonov,
*Impulsive Differential Equations: Periodic Solutions and Applications*, vol. 66 of*Pitman Monographs and Surveys in Pure and Applied Mathematics*, Longman Scientific & Technical, Harlow, UK, 1993. View at: MathSciNet

#### Copyright

Copyright © 2009 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.