Abstract

An integrated pest management Gompertz model with interval impulsive control is put forward. Through pest density monitoring, an integrated control strategy is adopted; that is when pest density reaches an environmental damage level, pesticide is used as a control method; when pest density is lower than the damage level, predators as its natural enemy are released in case pest density is higher than the slightly harmful level and predator density is below its maintainable level. The analysis on the existence of order-1 or order-2 periodic orbit is carried out by the construction of Poincaré map of semicontinuous dynamical system. The stability and attractiveness of the periodic state are obtained by geometry approach, which ensures a certain robustness of control. The analytical results presented in the work are validated by numerical simulations for a specific model.

1. Introduction

Agricultural pests are harmful to the crops, and thus pest management plays an important role in agricultural sustainable development, which also becomes an interesting and significant topic in real life. The traditional and efficient method for pest control is to spray pesticide, which can quickly kill the pest. However, unrestrained use of persistent pesticide not only increases the incidence of pesticide-resistant pest varieties but also inflicts harmful effects on humans through the accumulation of hazardous chemicals in their food chain [1]. An alternative way is the biological control by launching predators or enhancing predators’ genes to improve the effectiveness of the pest control [2]. This can be achieved by mass production and periodic release of natural enemies of the pest. However, the cultivation of the natural enemy in laboratories is uneasy, and the cost is very high in general. Integrated pest management (IPM) is the comprehensive utilization of agricultural, biological, chemical, and physical methods to control pests with the goal of controlling the number of the pests under an economic threshold (ET) [35]. Compared to natural growth process, spraying pesticides (or releasing the natural enemies) can cause the density of pests (or the natural enemies) to change sharply in a short time, which results in the discontinuity of the system states.

In the farmland and forest, the state of pests is always monitored. The administrative authority makes decision (i.e., whether the control should be taken and what kind of control should be adopted) according to the monitored pest level [6]. Based on the adopted control action, a series of integrated pest management models have been built in the literature, for example, the periodic release of predators [79]; the periodic release of pests infected by a disease [10]; the periodic release of predators and infected pests [11, 12]; the periodic release of infected pests combined with periodic applications of pesticides [13]; the periodic release of predators and pests combined with periodic applications of pesticides [1416]; and state-dependent release of predators combined with applications of pesticides [35, 1722].

As far as the state-dependent impulsive control is concerned, the earlier works in pest management belong to Tang et al. [35], where the biological and chemical controls are assumed to be taken at the pest economic threshold (ET). Notice that, in real applications, the biological control and chemical control sometimes are adopted at different pest levels; some researchers began to investigate the pest control models by assuming that releasing natural enemies and spraying pesticide are taken at different pest thresholds [2328]. From theoretical and practical points of view, the pest control model with twice impulsive controls looks more reasonable than the model with one impulsive control, and it can be seen as an extension of the one-impulsive-control models.

The ideas of involving biological and chemical controls at different prey densities is interesting and also has practical significance. Biological control with an appropriate yield of release of the predator is adopted in advance, to extend the time for pest density increasing to the damage level (i.e., reduce the operating frequency of adopting chemical control or pesticide). But there exists a problem in modelling the real system; that is, the biological control is adopted when the pest density reaches the first control level, but for a higher pest density between the biological and chemical control levels, there is no control strategy adopted. This is obviously unreasonable. Since the biological control and chemical control are activated at different pest levels, a more reasonable model should also consider the control action when the pest density lies between the two levels; that is to say, when the pest density increases to the chemical level, the chemical control has to be carried out, which causes a certain proportion to pest and predator to be killed. When the pest density increases to or exceeds the biological control activation level but is lower than the chemical control level, the biological control is sufficient, which is also necessary when the predator density is lower than its maintainable level, while in case of low pest density (i.e., below the biological control excitation level) it is not necessity to take any control action. Motivated by this control strategy, the current work presents and studies a prey-predator system involving interval state impulsive control.

This paper is organized as follows. In Section 2, a pest control prey-predator model with interval impulsive control is put forward. In Section 3, the Poincaré map, successor function, and some basic definitions are given, followed by a detailed dynamics analysis in case of the chemical control strength. In Section 4, numerical simulations are carried out with a specific model to verify the theoretical results step by step. Finally, conclusions are presented in Section 5.

2. Model Formulation and Preliminaries

Let and denote the pest and its natural enemy densities at time . Motivated by the control strategy, the following pest management Gompertz model is proposed:where denotes the intrinsic growth rate of the prey in the absence of the predator, denotes the environment carrying capacity, denotes the predation rate by natural enemies, denotes the transformation rate and denotes the death rate of the predator, and is a positive constant. The control parameters and represent the effect of pesticide to prey and predator species. Since the pesticide is relative to the pest, the effect on the predator is limited; here it is assumed that and . The continuous functions and () are, respectively, the predator reference level and the yield release of the predator, where .

The trajectory of the continuous system (1) passing through is described by , where is the first integral of the continuous system (1). Define the implicit function determined by the low branch above as (or ). Denote

Definition 1 (Poincaré map). For any trajectory in system (1) starting from point on , there will be two cases: for , the trajectory first reaches the point on section ; then, it jumps from to the point on due to the impulsive effects and and then jumps from to on due to the impulsive effects and ; for , the trajectory first reaches the point on section ; then, it jumps from to on due to the impulsive effects and and then the trajectory starting from reaches the point on section then jumps from to on due to the impulsive effects and . In any case, is a function of , , , , and . The Poincaré map on is defined as

Definition 2 (successor function). The successor function on is defined bywhich is continuous on , where is the Poincaré map determined by (3).

Definition 3 (periodic solution). The solution of system (1) passing through is called periodic if .

3. Dynamic Analysis of System (1)

As illustrated in Figure 1, system (1) has three equilibria , , and , where and if and holds. The reference level is established by a statistics method in practice, which reflects a certain consistency of the system; that means the biological control taken at any level should take the same or similar effect in determining the system’s dynamics. Thus, in this study a trajectory of the solution of system (1) starting from the point on is selected as a predator reference. To make practical sense, the study is restricted in the left region of ; that is, .

3.1. Existence of Periodic Solution

The discussions will be divided into three cases according to the magnitude between and : the position of the trajectory , the magnitude of and , and the location of the point .

Case 1 ( is greater than or equal to ). Suppose the intersection point between the line and the trajectory staring from is . The trajectory starting from the point between and intersects with at and then jumps to the point on under the effect of pulse. If the point , then it will jump upward to the point .
Case  1.1 ( overlaps with ). Then, the curve constitutes an order-1 periodic orbit (see Figure 2(a)).
Case  1.2 ( is under ). Then, there is . Next, choose a point adjacent to sufficiently (i.e., ). The trajectory of system (1) starting from intersects with at then jumps to under impulsive effects and and next jumps to above under impulsive effects and . Thus, there is (see Figure 2(b)), which means that an order-1 periodic solution exists.
Case  1.3 ( is above ). In this case, there is . Then, choose a point between and , which is close to sufficiently; the trajectory starting from intersects with at and then jumps to on under impulsive effects and . If , then it jumps to under impulsive effects and , and must be under ; then, there is (see Figure 2(c)); that is, an order-1 periodic solution exists. Otherwise, let be the point on such that . If , then an order-1 periodic solution exists.

Case 2 ( is smaller than ). Suppose the trajectory starting from intersects with the line at ; then, it jumps to the point under the impulsive effects and . And then the trajectory starting from intersects with at point and then jumps to at the point under the effect of pulses and . According to the relative position of and , there are three cases to be considered.
Case  2.1 ( overlaps with ). Then, the curve constitutes an order-2 periodic orbit (see Figure 3(a)).
Case  2.2 ( is under ). Then, there is . Next, choose a point adjacent to (i.e., ), where is the phase point of . The trajectory of system (1) starting from intersects with at point and then jumps to the point under the impulsive effects and . The trajectory starting from intersects with at point and then jumps to a point under the impulsive effects and . Since is close to , then must be above ; then, there is (see Figure 3(b)), which means that an order-2 periodic solution exists.
Case  2.3 ( is above ). In this case, there is . Let be a point on such that the trajectory starting from intersects with the line at . If the trajectory starting from does not intersect with , then choose the point as . Obviously there is . If the trajectory starting from intersects with at and with , then choose the point from and with smaller label as . The trajectory starting from intersects with at below and then jumps to . The trajectory starting from intersects with at and then jumps to . If , there is (see Figure 3(c)), which means that an order-2 periodic solution exists.

3.2. Stability of the Periodic Solutions

Theorem 4. For the case , if the order-1 periodic solution exists, it is unique and orbitally asymptotically stable.

Proof.
(1) The Uniqueness of the Order-1 Periodic Solution. Assume system (1) admits two different order-1 periodic solutions and which start from and on with , respectively. Then, there is . Denote , where and with . Then, there is which means that ; that is, is a monotone decreasing function on . Then, it yields that that is, . Then, there is which leads to a contradiction. Thus, the order-one periodic solution is unique.
(2) Stability of the Order-1 Periodic Solution. Suppose the impulsive point of is ; obviously, the trajectory of system (1) starting from any point will arrive at the segment . Thus, it is only necessary to consider the tendency of the trajectory starting from . If , we should consider the tendency of the trajectory starting from and , respectively. For any , a sequence on can be obtained such that ; and for any , a sequence on can be obtained such that . From the proof of the uniqueness of the order-one periodic, it can be observed that the successor function satisfies the following properties: (i) if and only if ;(ii) for , ;(iii) for , .This means that is a monotone increasing sequence with an upper bound , and is a monotone decreasing sequence with a lower bound ; then, there exist and such that and when . But implies that . Thus, the unique order-1 periodic is orbitally asymptotically stable. Since the points and are arbitrary, the orbit asymptotical stability implies the global attractiveness.

Theorem 5. For the case , if the order-2 periodic solution exists, it is unique, orbitally asymptotically stable, and globally attractive.

Proof.
(1) The Uniqueness of the Order-2 Periodic Solution. Assume system (1) admits two different period-2 solutions and which start from and with , respectively. Then, there is and which leads to a contradiction. Thus, the order-2 periodic solution is unique.
(2) The Stability of the Order-2 Periodic Solution. Assume that and are the impulsive points of and under the pulses and , respectively. The trajectory meets with at the point in the opposite direction; is the impulsive point of . Obviously, the trajectory of system (1) starting from any point will arrive at the segment if or if . Thus, it is only necessary to consider the tendency of the trajectory starting from or . Denote the order-2 periodic by with on the segment or . Similarly, let us consider the tendency of the trajectory starting from or and , respectively. For any or , a sequence can be obtained such that ; and for any , a sequence on can be obtained such that . From the proof of the uniqueness of the order-2 periodic solution, it can be observed that the successor function satisfies the following properties: (i) if and only if ;(ii) for , or , ;(iii) for , .This means that is a monotone increasing sequence with an upper bound , and is a monotone decreasing sequence with a lower bound ; then, there exist and such that and when . But implies that . Thus, the unique order-2 periodic solution is orbitally asymptotically stable. Since the points and are arbitrary, the orbit asymptotical stability implies the global attractiveness.

4. Numerical Simulations

In this section, a specific example is given to verify the theoretical results obtained in the previous section by considering the change of the control parameters , , and . Let , , , , , and . With a simple calculation, the steady state of free system is . The chemical control is assumed to be taken at , and the biological control is assumed to be taken at . The predator reference level at is assumed to be about of steady predator density of the free system; that is, . The curve intersects with and at the points and , respectively. Then, there is and .

4.1. Verification for the Case

Let . Then, there is . The predator reference level at is about . The curve intersects with at the point ; then, there is . For the yield release , an order-1 periodic solution is formed, as illustrated in Figure 4.

When increases to 0.2, that is, , the order-1 periodic solution disappears, and an order-2 periodic solution is formed, as illustrated in Figure 5.

For the yield release , an order-3 periodic solution is formed (Figure 6).

4.2. Verification for the Case

Let and . Then, . For the release , an order-2 periodic solution is formed (Figure 7).

For the yield release , an order-3 periodic solution is formed (Figure 8).

For the yield release , an order-4 periodic solution is formed (Figure 9).

For the yield release , an order-5 periodic solution is formed (Figure 10).

5. Conclusion

In this work, a pest control prey-predator model was analyzed, where the biological control and chemical control are considered at different thresholds. At the early stage of the pest damage outbreaks, the biological control is adopted in case of the predator density in the environment lower than its maintainable level. Once the pest density reaches a critical level, the chemical control with a given strength is taken since higher pest density may cause a serious damage to environment. Different to the models in literature, the proposed model is more consistent with practice.

The theoretical analysis indicated that the yield of releases of the predator plays a key role in determining the existence of order-1 periodic orbit. The practical significance to studying the existence of order-1 or order-2 periodic solution lies in that it could provide a possibility to determine the frequency of using chemical pesticide and yield of releases of the predator, which makes the control a periodic one without real-time monitoring of the species while keeping the prey density below the damage level. The stability and attractiveness could ensure a certain robustness of control; that is, even though the species density is detected inaccurately or with a deviation, the system will be eventually stable at the periodic solution under the control action.

Competing Interests

The authors declare that there is no conflict of interests.

Acknowledgments

This work was supported in part by the National Natural Science Foundation of China (nos. 11401068 and 61473327) and the Liaoning Province Natural Science Foundation of China (no. 2014020133).