Dynamics of a Stage Structure Pest Control Model with Impulsive Effects at Different Fixed Time
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.
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 : 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 ). 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.4 (see ). 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
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.
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.
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.
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.
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
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
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