Abstract

Stage-structured predator-prey models with disease in the prey are constructed. For the purpose of integrated pest management, two types of impulsive control strategies (impulsive release of infective prey and impulsive release of predator) are used. For Case  1, infective prey applications are more frequent than releases of predator (natural enemies). For Case  2, predator (natural enemies) releases are more frequent than infective prey applications. In both cases, we get the sufficient conditions for the global attractivity of the susceptible prey-eradication periodic solution. In addition, the persistence of the systems is also discussed. At last, the results are discussed and some possible future work is put forward.

1. Introduction

Pests, such as insects, mice and other animals, unwanted plants (weeds), fungi, microorganisms, and so forth, are living organisms that occur where they are not wanted or that cause damage to crops or other animals. How to minimize the loss caused by the injurious insects and injurious germ carrier to the important plants, animals, and human being is always the common problem concerned by the entomologists and society. Human has adopted some advanced and modern weapons such as chemical pesticides, biological pesticides, remote sensing and measuring, and so on to deal with pests, and some great achievements have been obtained [17].

The traditional chemical control only care about the current effect, but seldom take the influence on the ecosystem into consideration. And it caused many problems such as environment pollution, pest resistance to the pesticide and pest reemergence, and the like. In this regard, it has been observed that beneficial insects are often more susceptible to chemical pesticides than the target pests are. In the same time, the concentration of the pesticides in use tends to increase with time and usage, since many pests develop resistance to these chemicals. This kind of pest management strategy was considered by many authors [812]. At present more and more people are concerned about the effects of pesticide residues on human health and on the environment.

Compared to chemical treatment, nonchemical methods are safer to man and are generally effective for longer periods of time. One example of nonchemical pest control methods is biological treatment [1317], including microbial control with pathogens, as disease can be important natural controls of some pests. Insects, like humans and other animals, can be infected by disease-causing organisms such as bacteria, viruses, and fungi.

People also use natural enemy to control pest or regulate it to densities below the threshold for economic damage. Often with augmentation or release, the natural enemy is applied like a pesticide after the pest has reached or exceeded the economic threshold. There are many literatures concerning natural enemy for pest control [1825].

Many kinds of predator-prey models have been studied extensively [1, 2, 6, 20, 2628]. In the natural world there are many species whose individual members have a life history that takes them through two stages: immature and mature. In particular, we have in mind mammalian populations and some amphibious animals, which exhibit these two stages. In recent years, stage-structured models, with or without delays, have been studied by several authors [2633]. In addition, there are many control methods and results for complex dynamical network model [3438], from which we can learn for the proof of our main results.

Motivated by [1, 14, 24, 25, 33, 34], in this paper, we will consider predator-prey models with stage structure in the prey. The prey stands for the pest population and the predator stands for the natural enemy population. That is, we call the pest and natural enemy as prey and predator, respectively. Here, the pest population will be controlled by releasing natural enemy and infective pests together. The infective pests can be cultivated in the laboratory and the natural enemy can be migrated from other regions. Once the susceptible pest meets with the infective pest, there is a chance to be infected. The infective pests have more possibility of death due to the disease and have less damage to the crops and environment. In fact, there is such example: salt cedar leaf beetle is a pest, and it is hatched from eggs. We call the egg stage as immature pest, and mature pest after it is hatched. In view of its eggshell, pathogens may not be effective against pest eggs. That is, the disease only attacks the mature susceptible pest. Birds are the natural enemy of the beetle, and we call them predator.

The organization of this paper is as follows. In the next section, the main biological assumptions on which the models rely are formulated and the models are constructed. In Section 3, to prove our main results we give several definitions, notations, and lemmas. In Section 4, we analyze the first case and determine the sufficient conditions for the global attractivity of the susceptible pest-eradication periodic solution and permanence of the system (5). In Section 5, we analyze the second case by similar method and obtain the sufficient conditions for the global attractivity of the susceptible pest-eradication periodic solution and permanence of the system (6). In the last section, a brief discussion and some possible future work for pest management are provided.

2. Model Formulation

We assume that the life time of the prey population has two stages, immature stage and mature stage. represents the density of the immature prey (pest) population. , represent the densities of the susceptible and infective mature prey (pest) population, respectively. The predator population has only one life time stage. is the density of the predator (natural enemy) population. And the following assumptions hold. We suppose that the infective prey can neither produce offspring nor attack crops due to the disease, and only the susceptible prey can reproduce. At any time , birth into the immature prey population is proportional to the existing susceptible mature prey population with proportionality . The immature prey population will transfer to the mature prey class after its birth with a maturity period of . The term represents the immature prey that were born at time (i.e., ) and still survive at time (with the immature prey death rate ), and therefore represent the transformation from immature prey to mature prey. The immature prey population has the natural death rate . The parameters , are the death rate for the susceptible and infective mature prey, respectively. The predator population has the natural death rate . We also suppose the natural enemy only capture the susceptible mature prey population, since the immature and infective prey population are hidden in the sanctuary, and the predation functional response is type Holling II, is the conversion rate for predation. The mature prey is divided into two classes, susceptible and infective. The incidence rate is classic bilinear , and is the contact number per unit time for every infective prey with susceptible prey. We assume and are impulsive point series at which the infective prey and natural enemies are released, with the releasing amounts and , respectively. According to the above assumptions, we have two different cases as follows.

Case 1. The releases of infective pests are more frequent than releases of natural enemies.

We release the natural enemies impulsively with releasing amount at time . Assume holds for all , where is the period of releasing natural enemies. During the period, times of the infective pests are released. That is, there exists an integer , such that . For the convenience of calculation and research, we assume that . Accordingly, the mathematical model will be where we denote and throughout the paper.

Case 2. Natural enemies releases are more frequent than the release of the infective pests.

We release the infective pests impulsively with releasing amount . Assume that holds for all , where is the period of releasing the infective pests. During the period, times of natural enemy releases are applied. That is, there exists an integer , such that . For the simplicity of calculation and research, we assume that . Accordingly, the mathematical model will be For ecological reasons, we always assume that the initial values for system (1) and (2) satisfy where . For continuity of initial conditions, we require where represents the accumulated survivors of those prey members who were born between time and 0.

For the above two models, note that will not have impulsive perturbation, and it is not included in the other equations of system (1) and (2), so we only need to consider the following subsystems:

Accordingly, the initial condition of systems (5) and (6) becomes

3. Preliminary

We give some definitions and lemmas which will be useful for stating and proving our main results.

Definition 1. System (5) and (6) are said to be permanent if there are constants , (independent of the initial values) and a finite time such that for all solution with initial conditions (7), hold for all . Here may depend on the initial values.

Lemma 2 (see [39]). Consider the following equation: where are positive constants, and for , then one has the following.(i) If , then .(ii) If , then .

Lemma 3 (see [40]). Let and . Assume that where is continuous in and for each exists and is finite. is quasi-monotone nondecreasing in and is nondecreasing. Let be the maximal solution of the scalar impulsive differential equation defined on . Then implies that , , where is any solution of system (1) existing on .
Note that if one has some smoothness conditions of to guarantee the existence and uniqueness of solutions for (11), then is exactly the unique solution of (11).

Lemma 4 (see [1]). System has a positive periodic solution , and for every solution of system (12) with positive initial value , one has as , where for , and  .

Lemma 5. There exists a positive constant such that , , , for each solution of system (5) with positive initial values (7), where is large enough.

Proof. Define a function such that By simple computation, we see that when , , Since the right-hand side of the above inequality is a quadratic with negative quadratic coefficient, it is bounded from above for all . Hence there exists a positive constant , such that Obviously, we know
Consider system where ,  .
According to Lemmas 3 and 4, we derive Consequently, by the definition of we obtain that each solution of (5) with positive initial values is uniformly ultimately bounded above. This completes the proof.  

Remark 6. For system (6) with initial value (7), we have similar result as Lemma 5.

4. Analysis of System (5)

In this section, we determine the global attractive condition for the susceptible pest-eradication periodic solution and permanence of system (5).

If for all , then we get the following subsystem of system (5):

By Lemma 4, system (19) has a unique positive periodic solution as follows: where Thus, we get the presentation of the susceptible pest-eradication periodic solution of system (5) as . And we have the following results about this periodic solution.

Theorem 7. The susceptible pest-eradication periodic solution of system (5) is globally attractive provided that

Proof. Since (22) holds, we can choose a sufficiently small such that Let be any solution of system (5) with positive initial values. From system (5), we have According to Lemmas 3 and 4, for the above , there exists , such that hold for all .
Thus, for , we have Consider the following comparison equation: By Lemma 2, we have . And by Lemma 3, we have , for from which we get . Therefore, for any positive constant small enough , there exists an integer such that for all  .
From system (5), for all we have, From Lemmas 3 and 4, for any sufficiently small positive , there exists an integer such that for all , where are the solutions of the following comparison system: Since can be small enough, we have as . From (25) and (29), we deduce that as . Therefore, is globally attractive. This completes the proof.

Corollary 8. If , then the condition (22) of Theorem 7 becomes which means that if only natural enemies are released impulsively, then the release amount must be larger than to ensure the eradication of the susceptible pest.
If , then the condition (22) of Theorem 7 becomes which means that if only infective pests are released impulsively, then the release amount must be larger than to ensure the eradication of the susceptible pest.

Through Theorem 7 and Corollary 8, we can get the sufficient conditions for global attractivity of the susceptible pests-eradication periodic solution. That is, the susceptible pest population is eradicated totally. But in practice, from the view point of keeping ecosystem balance and preserving biological resources, it is not necessary to eradicate the susceptible pest population completely. In fact we hope the susceptible pests and natural enemies can coexist, and at the same time the susceptible pests do not cause immense economic loss. Thus, it is meaningful to study the permanence of system (5).

Theorem 9. The system (5) is permanent provided that

Proof. Since (34) holds, we can choose sufficiently small positive number ,   such that Suppose is a solution of system (5) with initial values . By Lemma 5, there exists a positive constant such that for all large enough. From (25), we know that are positive lower bound. Thus we only need to find such that for all large enough.
We claim that for any , it is impossible that for all . Suppose that the claim is not true. Then there is a such that for all .
Then, from system (5) we have From Lemmas 3 and 4, for the above selected in (35), there exists a time such that for all we have where are the solutions of the following comparison system: Note that the first equation of system (5) can be rewritten as In the following we define a Liapunov functional as Then the derivative of with respect to the solution of system (5) is From (37) and (41) we have Let .
We will show that for all . Otherwise there exists a nonnegative constant such that for , and . Thus from the first equation of system (5) and (35) we easily see that which is a contradiction. Hence we get that for all .
From (37) we have for all which implies as . This is a contradiction to . Therefore for any positive constant the inequality cannot hold for all  .
By the claim, we need to consider two cases.
Case  1.   for all large , then our aim is obtained.
Case  2.   oscillates about for all large . Denote In the following we will show that . Since oscillates about , there exist positive constants and , such that where is sufficiently large, and the inequalities (37) hold for . Since is continuous and bounded, and it is not effected by impulse. So there exists a constant (where and is independent of the choice of ) such that for all .
If , our aim is obtained.
If , from the first equation of system (5) we have, for then
It is clear that for .
If , then we have that for . Then, proceeding exactly as the proof for above claim, we see that for . Since this kind of interval is arbitrarily chosen (we only need to be large), we can get that for all large enough. In view of our arguments above, the choice of is independent of the positive solution of system (5) which satisfies that for sufficiently large . This completes the proof of the theorem.

Corollary 10. (1) If , then condition (34) of Theorem 9 becomes which means that if only natural enemies are released impulsively, and the release amount is less than , then system (5) is permanent. That is, the pest and natural enemy will coexist.
(2) If , then condition (34) of Theorem 9 becomes which means that if only infective pests are released impulsively, and the release amount is less than , then system (5) is permanent. That is, the pest and natural enemy will coexist.

5. Analysis of System (6)

In this section we will discuss system (6), the condition for the global attractive of the susceptible pest-eradication periodic solution, and the permanence of system (6) will be obtained.

If for all , then we get the following subsystem of system (6): By Lemma 4, we know that the system (51) has a unique positive periodic solution as follows:

Thus, system (6) has a unique nonnegative periodic solution , which is called as the susceptible pest-eradication periodic solution. Next, we will discuss the global attractivity of this periodic solution and the permanence of system (6).

Theorem 11. The susceptible pest-eradication periodic solution of system (6) is attractive provided that

Proof. Since (53) holds, we can choose a sufficiently small such that Let be any solution of system (6) with initial values (7). From system (6), we have From Lemmas 3 and 4, for the above selected , there exists , such that for all we have Thus, for , we have Consider the following comparison equation: By Lemma 2, we have and by Lemma 3, we have , for large , that is . Therefore, for any positive constant small enough (, there exists an integer such that for all .
Thus, from system (6), for all we have From Lemmas 3 and 4, for a sufficiently small , there exists an integer such that for all , whereis the solution of the following comparison system: Since can be small enough, we have as . From (56) and (60) we get as . Together with , we obtain that is globally attractive. This completes the proof.

Corollary 12. If , then condition (53) of Theorem 11 becomes which means that if only natural enemies are released impulsively, then the release amount must be larger than to ensure the eradication of the susceptible pest.
If , then condition (53) of Theorem 11 becomes which means that if only infective pests are released impulsively, then the release amount must be larger than to ensure the eradication of the susceptible pest.

Through Theorem 11 and Corollary 12, we can get the sufficient conditions for the susceptible pests-eradication periodic solution; that is, the susceptible pest population is eradicated totally. But in practice, from the view point of keeping ecosystem balance and preserving biological resources, it is not necessary to eradicate the susceptible pest population. In fact we hope the susceptible pests and enemies can coexist when the susceptible pests do not cause immense economic loss. Thus, it is necessary to discuss the permanence of system (6).

Theorem 13. System (6) is permanent provided that

Proof. Since (65) holds, we can choose sufficiently small positive number () such that Suppose is any solution of system (6) with initial values (7). By Lemma 4, there exists a positive constant such that for large enough. From (56), we know that are positively lower bounded.
Thus we only need to find such that for large enough. We claim that for any , it is impossible that for all . Suppose that the claim is not true. Then there is a such that for all . Then, from system (6) for all we have From Lemmas 3 and 4, for the above selected , there exists a time such that hold for all , where are the solution ofNote that the first equation of system (6) can be rewritten as In the following we define a Liapunov functional as Then the derivative of with respect to the solution of system (6) is From (68) and (72) we have Let .
We will show that for all . Otherwise there exists a nonnegative constant such that for , and . Thus from the first equation of system (6) and (66) we easily see that which is a contradiction.  Hence we get that for all . From (68), for all we have which implies as . This is a contradiction to . Therefore for any positive constant the inequality cannot hold for all .
By the claim, we need to consider two cases.
Case  1. for all large , then our aim is obtained.
Case  2. oscillates about for all large . Denote In the following we will show that . Since oscillates about , there exist positive constants and , such that where is sufficiently large, and the inequalities (68) hold for . Since is continuous and bounded, and it is not effected by impulse, there exists a constant (where and is independent of the choice of ) such that for all .
If , then our aim is obtained.
If , from the first equation of system (6) we have and then we have and it is clear that for .
 If , then we have that for . Then, proceeding exactly as the proof for above claim, we see that for . Since this kind of interval is arbitrarily chosen (we only need to be large), we can get that for all large enough. In view of our arguments above, the choice of is independent of the positive solution of system (6) which satisfies that for sufficiently large . This completes the proof of the theorem.

Corollary 14. If , then condition (65) of Theorem 13 becomes which means that if only natural enemies are released impulsively, and the release amount is less than , then system (6) is permanent. That is, the pest and natural enemy will coexist.
If , then the condition (65) of Theorem 13 becomes which means that if only infective pests are released impulsively, and the release amount is less than , then system (6) is permanent. That is, the pest and natural enemy will coexist.

6. Discussion

In this paper, a stage-structured predator-prey model with disease in the prey is considered. The prey stands for the pest population, and the predator stands for the natural enemy population. For the purpose of integrated pest management, two types of impulsive control strategies are used. In Case 1, infective prey applications are more frequent than releases of natural enemies. In Case 2, natural enemies releases are more frequent than releases of infected prey. In both cases, we analyzed the global attractivity of the susceptible pest-eradication periodic solution of the systems, and we also obtained the condition for the permanence of the systems.

In Section 4, the first case is analyzed. By Theorem 7, the sufficient condition for the global attraction of the susceptible pest-eradication periodic solution is obtained, which means that if the release amount of infective pest and natural enemy satisfy certain conditions, then the susceptible pest will be doomed. By the result of Theorem 9, the sufficient condition for the permanence of system (5) is also obtained, which means that the pest and the natural enemy will coexist for all time, if the release amounts of infective pest and natural enemy meet some critical values. Corollary 8 and 10 show that if only one control measure is taken, that is, either only infective pests are released or only natural enemies are released, then the release amount must satisfy certain conditions.

In Section 5, we use similar method to analyze the second case. In Theorem 11, the sufficient condition for the global attraction of the susceptible pest-eradication periodic solution is obtained, which means that if the release amounts of infective pest and natural enemy meet certain conditions, then the susceptible pest will be doomed. By the result of Theorem 13, the sufficient condition for the permanence of system (6) is also obtained, which means that the pest and the natural enemy will coexist for all time, if the release amounts of infective pest and natural enemy meet some critical values. Corollary 12 and 14 show that if only one measure is taken, that is, either only infective pests are released or only natural enemies are released, then the release amount must satisfy certain conditions.

From the above theoretical results obtained in Sections 4 and 5, we can choose different control strategies to control pest. We can do it by releasing infected pest more frequently than releasing natural enemy, or vice versa. However, in practice, every control measure needs certain cost, and we must think about the cost before adopting a measure. For example, the infective prey (pest) being cultivated in the laboratory and the natural enemy (predator) being migrated from other regions all need some cost. We can choose different pest management methods with different costs. In both models, different parameters mean different cost for control measures. Therefore, we have an interesting problem. In the view point of economy, under the premise of controlling the pest, which kind of control method is more suitable? This will be an optimal control problem, and we leave it as our future work.

Acknowledgments

The authors would like to thank the referees for their helpful suggestions, which improved the quality of this paper greatly. The first author is supported by Postdoctoral Science Foundation of China (no. 2011M501428). The second author is supported by NNSF of China (nos. 11171199, 81161120403).