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.


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.

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 x j t and x t 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.α > 0 is the coefficient of birth rate of the mature population in this environment; τ, called maturation time delay, is the time to maturity; d > 0 is the coefficient of death rate of the immature population; β > 0 is the mature pest death and overcrowding rate; the term αe −dτ x t − τ represents the immature pests which are born at the time t − τ i.e., αx t − τ and survive at the time t with the immature pest death rate d 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: > 0 is the natural enemy death's rate; μ > 0 is the natural enemy releasing amount at every impulsive period nT ; 0 ≤ p ≤ 1 is the killed rate of natural enemies at every impulsive period n k −1 T .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 x j 0 0 −τ αe −ds x s ds; that is x j t can be linear expression by x t , and x j t 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: y t y t μ, t nT.

2.4
The initial conditions for 2.4 are From the biological point of view, we only consider system 2.4 in the biological meaning region: D {x t , y t | x t ≥ 0, y t ≥ 0}.
Lemma 2.1 see 22 .Consider the following equation: where a, b, c, τ are all positive constants, x t > 0, for −τ ≤ t ≤ 0; one has the following.

2.8
So V t is uniformly ultimately bounded.Hence, by the definition of V t , we obtain that each positive solution of system 2.2 is uniformly ultimately bounded.The proof is complete.where d 2 , p, μ, k are all positive constants; then, system 2.9 has a globally asymptotic stable positive periodic solution:
Proof.From the second equation in system 2.4 , we have dy t /dt ≥ −d 2 y t ; then consider the following system:

3.1
According to Lemma 2.3, the globally asymptotic stable positive periodic solution of system 3.1 is

3.2
According to Lemma 2.4, for any 1 > 0, there exits a n 1 such that From the first equation of system 2.4 , we can get Then consider the following comparison equation: Since ζ 1 < 1, according to Lemma 2.1, we obtain lim t → ∞ z t 0; then by the comparison theorem in differential equations, we get lim t → ∞ x t 0. Without loss of generality, we assume that 0 < x t < , ∀t ≥ 0.

3.6
Consider the following system:

3.7
The globally asymptotic stable positive periodic solution of system 3.7 is

3.8
By 3.3 and Lemma 2.4, for any 2 > 0, there exits a T * , when t > T * : 3.9 Let → 0, we can get for t being large enough, which implies lim t → ∞ y t y * t .This completes the proof.

Permanence
Define

4.5
According to Lemma 2.4, there exists T 1 ≥ t 0 τ such that for t ≥ T 1 , From 4.3 and 4.6 we have Let m 2 min t∈ T 1 ,T 1 τ x t , we show that x t ≥ m 2 for all t > T 1 .Otherwise, there exists a nonnegative constant T 2 , such that Thus, from the first equation of system 2.4 and 4.5 , we easily see that

4.8
This is a contradiction.So we obtain that x t ≥ m 2 , for all t > T 1 .From 4.7 , we have Therefore, for any positive constant t 0 , the inequality x t < m 1 cannot hold for all t ≥ t 0 .If x t ≥ m 1 holds true for all t large enough, then our aim is obtained.Othervise, x t is oscillatory about m 1 .Let

4.10
In the following, we will show that x t ≥ q for t being large enough.There exist two positive constants t, ω such that x t < m 1 for t < t < t ω.

4.11
Since x t is continuous and bounded and is not effected by impulses, we conclude that x t is uniformly continuous.Then there exists a constant T 3 , such that x t > m 1 /2 for all t ≤ t ≤ t T 3 .
If ω ≤ T 3 , our aim is obtained.If T 3 < ω ≤ τ, from the first equation of 2.4 we have that Then we have It is clear that x t ≥ q for t < t ≤ t ω.
If ω > τ, by the first equation of 2.4 , then we have that x t ≥ q for t < t ≤ t τ.Thus, proceeding exactly as the proof for above claim, we can obtain the inequality x t < m 1 cannot hold for all t > t τ, so the same arguments can be continued for for t τ < t ≤ t ω, and we can get x t ≥ q for t τ < t ≤ t ω.Since the interval t, t ω is arbitrarily chosen we only need t to be large , we get that x t ≥ q for t being large enough.In view of our arguments above, the choice of q is independent of the positive solution of 2.4 which satisfies that x t ≥ q for being sufficiently large t.This completes the proof.
Proof.Suppose that x t , y t is any solution of system 2.4 with initial condition 2.5 .Let q * μ 1−p e −d 2 T /1− 1−p e −d 2 T , by 3.3 we know y t ≥ q * , and according to Theorem 4.2, there exist positive constants q, such that x t ≥ q.Set Then, D 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 D. 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.Proof.From 3.3 we know for any 1 > 0 that there exits a n 1 such that From the first equation of system 2.

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 T, τ, μ, p.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 d 0.6; d 1 0.4; α 0.8; τ 0.8; d 2 0.5; β 0.8; λ 0.9; p 0.3; μ 0.2; T 0.2 see Figure 1 ; here we can compute ζ 2 1.253 > 1; from Theorem 4.3 we know that system 2.4 is permanent.If we increase the natural enemy input amount to μ 0.5 ζ 1 0.716 < 1 or increase the maturation time delay τ 2 ζ 1 0.872 < 1 , 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 p 0.1, then ζ 1 0.635 < 1, 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, E, is 0.3, then Figure 5 shows that if we choose an appropriate pulse releasing natural enemy period with T , we can control the pest population under the ETL and not eradicate natural enemy totally.This gives us some reasonable suggestions for pest management.
0 such that x t < m 1 for all t ≥ t 0 .From the second and the fourth equations in system 2.4 , we have ζ 2 1 − 1 − p e −d 2 T αe −dτ /βμe −d 2 T .Definition 4.1.System 2.4 is said to be permanent if there are constants l, L > 0 and a finite time T 0 such that for every positive solution x j t , x t , y t ∈ R 3 with initial conditions 2.5 satisfies l ≤ x j t ≤ L, l ≤ x t ≤ L, l ≤ y t ≤ L for all t ≥ T 0 .Here T 0 may depend on the initial condition 2.5 .Theorem 4.2.If ζ 2 > 1 holds, there exists a positive constant q such that each positive solution x t , y t of 2.4 satisfies x t ≥ q with t being large enough.Proof.The first equation of system 2.4 may be rewritten as follows: − d 2 − λβm 1 T / 1 − 1 − p e − d 2 −λβm 1 T − .For any positive constant t 0 > 0, we claim that the inequality x t < m 1 cannot hold for all t ≥ t 0 .Otherwise, there is a positive constant t