# The Effects of Resource Limitation on a Predator-Prey Model with Control Measures as Nonlinear Pulses

Authors: Wenjie Qin; Sanyi Tang; Robert A. Cheke
Journal: Mathematical Problems in Engineering
 (2014)

---

## Abstract

The dynamical behavior of a Holling II predator-prey model with control measures as nonlinear pulses is proposed and analyzed theoretically and numerically to understand how resource limitation affects pest population outbreaks. The threshold conditions for the stability of the pest-free periodic solution are given. Latin hypercube sampling/partial rank correlation coefficients are used to perform sensitivity analysis for the threshold concerning pest extinction to determine the significance of each parameter. Comparing this threshold value with that without resource limitation, our results indicate that it is essential to increase the pesticide’s efficacy against the pest and reduce its effectiveness against the natural enemy, while enhancing the efficiency of the natural enemies. Once the threshold value exceeds a critical level, both pest and its natural enemies populations can oscillate periodically. Further-more, when the pulse period and constant stocking number as a bifurcation parameter, the predator-prey model reveals complex dynamics. In addition, numerical results are presented to illustrate the feasibility of our main results.


---

## Body

## 1. Introduction

It is well known that pest outbreaks often cause serious ecological and economic problems, requiring complex control measures to reduce harm due to insect pests of agriculture and insect vectors of important plant, animal, and human diseases. Such measures include use of a variety of chemical pesticides, biological pesticides, and biological control.Biological control is the reduction of pest populations by other living organisms, often called natural enemies or beneficial species (see [1–3]). Virtually all pests have some natural enemies, and the key to successful pest control is to identify the pest and its main natural enemies and to release the beneficial insects early when pest levels are low. Chemical control relies mainly on the use of synthetic pesticides to suppress pests. Pesticides are useful because they quickly kill a significant portion of a pest population and they sometimes provide the only feasible method for preventing economic loss. However, pesticide pollution is also recognized as a major health hazard to people and to pest’s natural enemies. Also, overuse of a single control tactic is discouraged to avoid or delay the development of resistance by the pest to the control tactic, to minimize damage to nontarget organisms, and to preserve the quality of the environment. Therefore, it is natural to combine biological and chemical controls as components of integrated pest management (IPM).The concept of IPM was introduced in the late 1950s [4], was widely practised during the 1970s and 1980s [5], and is still often adopted. IPM emphasizes the importance of interactions between pests and its natural enemies and is a long-term management strategy to reduce pests to predetermined economic injury levels, with little cost and minimal effects on the environment. IPM has been shown to be more effective than the traditional methods, such as biological control or chemical control alone, both experimentally [6, 7] and theoretically [8, 9]. These results indicate that different pest control techniques should work together rather than against each other.In the last few decades, in order to consider the consequences of spraying pesticides and introducing additional predators into a natural predator-pest system, many authors have suggested impulsive differential equations to investigate the dynamics of pest control models. For the general theory of impulsive differential equations, see [10, 11]. Recently, many papers have been devoted to the analysis of mathematical models describing IPM strategies. See, for instance, [12–20] and so on.The emphasis of IPM is on control, not eradication of the pest, as the latter can be impossible, cost-ineffective, or potentially damaging to the environment. The above traditional predator-prey model with IPM has assumed that the fatality rate of pesticide applications with respect to the pest is constant, which implies that the agricultural resources such as pesticides, labor forces, equipments, and costs are very effective and sufficient for controlling pests. However, in reality, every community or country has an appropriate or limited capacity for pesticides, costs, and so forth. Meanwhile, the wide use of the pesticides will potentially damage the environment. A resource is limiting if changes in its availability affect the population equilibrium level [21], which is a function of individual survival and reproduction. Estimates of demography, therefore, can be linked with estimates of resource availability, as they change through time, to yield insights into the effects of resource limitation [22].Understanding resource limitation is critical to effective management and conservation of populations; however, resource limitation is difficult to quantify partly because it is a dynamic process [23]. In order to investigate the effect of limited resources on the outbreak of pest populations, a saturation phenomenon of the limited resources is considered. That is, we will study the dynamic behaviors of a predator-prey model with nonlinear pulse controls, which are suitable mathematical models to simulate processes with short-term perturbations during their development [11]. The main purposes of this paper are to construct a simple mathematical model including the features of periodic biological and chemical control for pest control to understand how limited resources affect pest outbreaks since the limited resources and the fatality rate of pesticide applications on the pest could depend on its population’s density. Also, in order to investigate the effect of the limited capacity of pesticides, a nonlinear continually differentiable function to characterize the saturation phenomenon of the limited resources is introduced. The model is based on the predator-prey model with generalized Holling II functional response.The organization of present paper is as follows. In the next section, a Holling II predator-prey model with nonlinear pulse is introduced. In Section3, by using the method of the differential inequality, qualitative analysis, and the discrete dynamical system determined by the stroboscopic map and Floquet theory, a set of sufficient conditions which guarantee the existence and stability of the pest-free periodic solution are obtained. Meanwhile, we explored the parameter space by performing an uncertainty analysis and sensitivity analysis using the Latin hypercube sampling (LHS) method and evaluating partial rank correlation coefficients (PRCCs) for various input parameters against threshold conditions and then the key factors which are most significantly related to the threshold conditions were determined (see [18, 24, 25]). Section 4 focuses on the sufficient conditions under which the positive periodic solution exists by using the bifurcation theorem. The paper ends with some interesting biological conclusions and numerical bifurcation analyses, which complement the theoretical findings.


## 2. Model Formulation

The predator-prey dynamics are modelled by two ordinary differential equations (see [26, 27]). Based on experiments, Holling [28] suggested three different kinds of functional responses for different kinds of species to model the phenomena of predation, which made the standard Lotka-Volterra model more realistic. So the basic model we consider is based on the following predator-prey model where the prey has logistic growth with no predators and the predator has Holling II functional response:
(1)
x
˙
=
a
x
-
b
x
2
-
c
x
y
1
+
ω
x
,
y
˙
=
μ
c
x
y
1
+
ω
x
-
d
y
,
where a, b, c, d, ω, μ are positive constants. x
(
t
), y
(
t
) are the densities of the prey and predator at time t, respectively.With the idea of impulsive perturbations, motivated by [29, 30], a predator-prey model concerning IPM strategies is proposed as follows:
(2)
x
˙
=
a
x
-
b
x
2
-
c
x
y
1
+
ω
x
,
y
˙
=
μ
c
x
y
1
+
ω
x
-
d
y
,
t
≠
n
T
,
n
∈
𝒩
,
x
(
t
+
)
=
(
1
-
p
1
)
x
,
y
(
t
+
)
=
(
1
-
p
2
)
y
+
σ
,
t
=
n
T
,
n
∈
𝒩
,
where 𝒩
=
{
0,1
,
2
,
…
} and T is the pulse period. 0
≤
p
1
<
1
(
0
≤
p
2
<
1
) represents the fraction of pests (predators) which die due to the pesticide at t
=
n
T, and σ
>
0 is the release amount of predator at t
=
n
T.It is assumed that the fatality rate of pesticide applications with respect to the pest (natural enemy) is a constantp
1
(
p
2
) after the spraying of pesticides every time for model (2). To take account of the resource limitation and saturation effects, we use the Hill function [31], and then model (2) can be rewritten as
(3)
x
˙
=
a
x
-
b
x
2
-
c
x
y
1
+
ω
x
,
y
˙
=
μ
c
x
y
1
+
ω
x
-
d
y
,
t
≠
n
T
,
n
∈
𝒩
,
x
(
t
+
)
=
(
1
-
p
1
max
⁡
x
(
t
)
x
(
t
)
+
θ
1
)
x
(
t
)
,
y
(
t
+
)
=
(
1
-
p
2
max
⁡
y
(
t
)
y
(
t
)
+
θ
2
)
y
(
t
)
+
σ
,
t
=
n
T
,
n
∈
𝒩
with nonnegative initial value (
x
0
,
y
0
), and p
1
max
⁡
(
p
2
max
⁡
) and θ
1
(
θ
2
) represent the maximal fatality rate and the half saturation constant for pest (natural enemy), respectively. That is, we can use a combination of biological (periodic releasing of natural enemies) and chemical (spraying pesticide) tactics to suppress the pest to a low level.


## 3. Mathematical Analysis of the Pest-Free Periodic Solution and Threshold Conditions

In this section, we first demonstrate the existence of the pest-free periodic solution of model (3) and focus on determining its global attractivity. Moreover, threshold values are derived, and LHS/PRCCs analysis techniques are employed to investigate the key control parameters which are most significantly related to threshold values.


### 3.1. Existence of the Pest-Free Periodic Solution

In order to illustrate the existence of the pest-free periodic solution, we first focus on the pest-free set:{
(
x
,
y
)
∈
ℛ
+
2
,
x
=
0
}, which is clearly invariant by model (3); within this set, model (3) becomes
(4)
y
˙
=
-
d
y
(
t
)
,
t
≠
n
T
,
y
(
t
+
)
=
(
1
-
p
2
max
⁡
y
(
t
)
y
(
t
)
+
θ
2
)
y
(
t
)
+
σ
,
t
=
n
T
,
y
(
0
+
)
=
y
0
.
Solving the above equation in interval (
n
T
,
(
n
+
1
)
T
], one yields
(5)
y
(
(
n
+
1
)
T
+
)
=
(
1
-
p
2
max
⁡
y
(
n
T
+
)
e
-
d
T
y
(
n
T
+
)
e
-
d
T
+
θ
2
)
y
(
n
T
+
)
e
-
d
T
+
σ
.
Denote y
(
n
T
+
)
=
y
n; then the above equation can be rewritten as the following difference equation:
(6)
y
n
+
1
=
(
1
-
p
2
max
⁡
)
e
-
2
d
T
y
n
2
+
θ
2
e
-
d
T
y
n
e
-
d
T
y
n
+
θ
2
+
σ
≐
F
(
y
n
)
,
which is the so-called stroboscopic map of model (4), and describes the relations of the number of natural enemies between any two successive pulse points. Consequently, the existence of the positive steady state of model (6) implies the existence of a positive periodic solution of model (4). Therefore, we first discuss the existence of a positive steady state of (6); that is, the derivative of F
(
y
n
) with respect to y
n yields
(7)
F
′
(
y
n
)
=
[
(
1
-
p
2
max
⁡
)
+
p
2
max
⁡
θ
2
2
(
e
-
d
T
y
n
+
θ
2
)
2
]
e
-
d
T
and it is easy to see that 0
<
F
′
(
y
n
)
<
1 holds true for all y
n.Next, we discuss the positive fixed point of the stroboscopic map (6), denoted by y
*, which satisfies the following equation:
(8)
A
1
y
*
2
+
B
1
y
*
+
C
1
=
0
,
where A
1
=
exp
⁡
(
-
d
T
)
[
1
-
(
1
-
p
2
max
⁡
)
exp
⁡
(
-
d
T
)
]
>
0, B
1
=
θ
2
-
σ
exp
⁡
(
-
d
T
)
-
θ
2
exp
⁡
(
-
d
T
), C
1
=
-
θ
2
σ
<
0.DenoteΔ
=
B
1
2
-
4
A
1
C
1, and by simple calculation, we have
(9)
Δ
=
(
θ
2
-
σ
e
-
d
T
-
θ
2
e
-
d
T
)
2
+
4
θ
2
σ
e
-
d
T
[
1
-
(
1
-
p
2
max
⁡
)
e
-
d
T
]
>
0
.
Obviously, (6) has a unique positive root; that is,
(10)
y
*
=
-
B
1
+
B
1
2
-
4
A
1
C
1
2
A
1
,
which is stable due to 0
<
F
′
(
y
*
)
<
1.According to the relations between a fixed point of the stroboscopic map (6) and the periodic solution of model (4), we conclude that the submodel (4) of (3) has a unique nontrivial positive periodic solution, denoted by y
p
(
t
) and
(11)
y
p
(
t
)
=
y
*
exp
⁡
[
-
d
(
t
-
n
T
)
]
for
n
T
<
t
≤
(
n
+
1
)
T
,
and it follows from 0
<
F
′
(
y
n
)
<
1 and [32] that we have the following theorem.Theorem 1.
Model (4) has a positive periodic solution y
p
(
t
) with periodic T and for every solution y
(
t
) of (4) such that |
y
(
t
)
-
y
p
(
t
)
|
→
0 as t
→
+
∞, where y
p
(
t
)
=
y
*
e
-
d
(
t
-
n
T
) and y
* is defined in (10).By using Theorem1, we conclude that the unique fixed point y
* of (6) is globally stable. Accordingly, the stability of the positive periodic solution is determined by the stability of the positive equilibria of the difference equation made by the sequence of impulsive points. So, the positive periodic solution y
p
(
t
) of model (4) is globally stable.Therefore, we obtain the general expression of an unique pest-free periodic solution of model (3) over the interval n
T
<
t
≤
(
n
+
1
)
T for all n
∈
𝒩 denoted by
(12)
(
x
p
(
t
)
,
y
p
(
t
)
)
=
(
0
,
y
*
e
-
d
(
t
-
n
T
)
)
.


### 3.2. Stability of the Pest-Free Periodic Solution

The threshold conditions which guarantee the global stability of the pest-free periodic solution (12) play a key role in pest control. To show this, we now concentrate on the stability of the pest-free solution (12) for model (3). We change model (3) variables to consider the deviations from the pest-free solution that are denoted by (
x
¯
(
t
)
,
y
¯
(
t
)
) so that
(13)
(
x
¯
(
t
)
,
y
¯
(
t
)
)
=
(
x
(
t
)
,
y
(
t
)
)
-
(
x
p
(
t
)
,
y
p
(
t
)
)
that yields
(14)
x
¯
˙
=
a
x
¯
-
b
x
¯
2
-
c
x
¯
1
+
ω
x
¯
[
y
¯
+
y
p
(
t
)
]
,
y
¯
˙
=
μ
c
x
¯
1
+
ω
x
¯
[
y
¯
+
y
p
(
t
)
]
-
d
y
¯
.
Assuming that (
x
¯
,
y
¯
) are small enough, we get the linear approximation of the deviation system around the periodic solution (
x
p
(
t
)
,
y
p
(
t
)
),
(15)
x
¯
˙
=
(
a
-
c
y
p
(
t
)
)
x
¯
,
y
¯
˙
=
μ
c
y
p
(
t
)
x
¯
-
d
y
¯
.Before stating Theorem2, for the sake of convenience, we set
(16)
R
0
1
=
a
d
T
c
y
*
(
1
-
e
-
d
T
)
,
R
0
2
=
c
d
T
M
s
c
y
*
(
1
-
e
-
d
T
)
-
d
ln
⁡
(
1
-
p
1
max
⁡
)
,
where M
s
=
(
a
ω
+
b
)
2
/
(
4
b
c
ω
).We are now in a position to present the sufficient conditions for the global stability of the pest-free periodic solution(
x
p
(
t
)
,
y
p
(
t
)
) of model (3) as follows.Theorem 2.
The pest-free periodic solution(
x
p
(
t
)
,
y
p
(
t
)
) of model (3) is locally stable in the first quadrant, provided that R
0
1
≤
1 and is globally stable if max
⁡
{
R
0
1
,
R
0
2
}
≤
1.The proof of Theorem2 is given in Appendix A.Remark 3.
According to the formulaeR
0
1 and R
0
2, we have
(17)
R
0
1
-
R
0
2
=
d
T
[
y
*
(
1
-
e
-
d
T
)
(
a
2
-
c
2
M
s
)
-
a
d
ln
⁡
(
1
-
p
1
max
⁡
)
]
c
y
*
(
1
-
e
-
d
T
)
[
a
y
*
(
1
-
e
-
d
T
)
-
d
ln
⁡
(
1
-
p
1
max
⁡
)
]
.
Therefore, if y
*
(
1
-
e
-
d
T
)
(
a
2
-
c
2
M
s
)
-
a
d
ln
⁡
(
1
-
p
1
max
⁡
)
>
0, then R
0
1
>
R
0
2, and consequently the pest-free periodic solution (
x
p
(
t
)
,
y
p
(
t
)
) of model (3) is globally stable if R
0
1
<
1, which shows that the local stability indicates the global stability. If y
*
(
1
-
e
-
d
T
)
(
a
2
-
c
2
M
s
)
-
a
d
ln
⁡
(
1
-
p
1
max
⁡
)
<
0, then R
0
1
<
R
0
2 (see Figure 1(a)). Thus, the pest-free periodic solution (
x
p
(
t
)
,
y
p
(
t
)
) of model (3) is globally stable if R
0
2
<
1 (see Figure 1(b)). If so, we cannot show that in this case the local stability indicates the global stability theoretically.The monotonicity ofR
0
1 and R
0
2 with respect to T. (a) R
0
1
>
R
2. All parameters are fixed as follows: a
=
0.5, b
=
0.8, c
=
1.6, ω
=
1, μ
=
0.4, d
=
0.37, σ
=
0.78, p
1
max
⁡
=
0.95, p
2
max
⁡
=
0.2, θ
1
=
1, θ
2
=
2. (b) R
0
1
<
R
2. All parameters are fixed as follows: a
=
0.1, b
=
0.8, c
=
2.5, ω
=
1, μ
=
0.4, d
=
0.37, σ
=
1, p
1
max
⁡
=
0.95, p
2
max
⁡
=
0.1, θ
1
=
1, and θ
2
=
2.
(a)
(b)


### 3.3. Threshold Conditions

Note that although the threshold valuesR
0
i depend on all parameters of model (3), the most interesting parameters here are the pulse period T and some parameters related to resource limitation such as the maximum fatality rates p
i
max
⁡
(
i
=
1,2
) and the half saturation constant θ
2, so it makes sense to know how T, p
i
max
⁡, θ
2 affect the threshold conditions (R
0
i
<
1; i
=
1,2) which guarantee the global stability of the pest-free periodic solution (12). Figure 1(a) indicates that R
0
i
(
i
=
1,2
) are monotonic increasing functions with respect to T
∈
(
0,20
]. Moreover, there are some critical values of T such that R
0
i
=
1, which indicate that there exists one maximum allowable impulsive period T
i
max
⁡ such that R
0
i
<
1 for all T
<
T
i
max
⁡ and i
=
1,2, respectively. Meanwhile, R
0
2 takes more time to reach the threshold value 1 than R
0
1 as a result of T
1
max
⁡
<
T
2
max
⁡. In more detail, the contour plots of Figures 2(a) and 2(c) show that the control strategies (decreasing p
2
max
⁡, increasing θ
2) can reduce the threshold value R
0
i
(
i
=
1,2
). However, increasing p
1
max
⁡ can reduce the threshold value R
0
2 from Figure 2(d). Therefore, in order to control the pest population, it is very important to increase the maximum fatality rate p
1
max
⁡ for the pest while decreasing the maximum fatality rate p
2
max
⁡ for the natural enemy, which means that it is essential to increase the pesticide’s efficacy and reduce its effectiveness against the natural enemy while enhancing the efficiency of the natural enemies. Also, it is important to carry out pulsed control actions more and more frequently under resource limitation according to Figure 2(b), and then we can eradicate the pest successfully.Contour plots ofR
0
1 and R
0
2: (a) plot of contours of R
0
1 versus p
2
max
⁡ and θ
2; (b) plot of contours of R
0
2 versus T and p
1
max
⁡; (c) plot of contours of R
0
2 versus p
2
max
⁡ and θ
2; (d) plot of contours of R
0
1 versus p
2
max
⁡ and p
1
max
⁡. All other parameters are fixed as follows: a
=
0.5, b
=
0.8, c
=
1.6, ω
=
1, μ
=
0.4, d
=
0.37, σ
=
0.78, p
1
max
⁡
=
0.95, p
2
max
⁡
=
0.2, θ
1
=
1, θ
2
=
2, T
=
6.5
(
A
), and T
=
12
(
C
-
D
).
(a)
(b)
(c)
(d)


### 3.4. Sensitivity to Variations

Previous analysis indicates thatR
0
i
(
i
=
1,2
) are significant threshold conditions affecting outbreaks of the pest population. Sensitivity analysis of the most significant parameters (such as intrinsic growth rate a, death rate c, the maximum fatality rate p
i
max
⁡, half saturation constant θ
2, constant stocking number σ, and pulse period T) was performed by evaluating the PRCCs for various input parameters against the threshold condition R
0
i
(
i
=
1,2
) with the LHS method (see [24, 25, 33]). PRCCs show how influential each parameter is on R
0
i
(
i
=
1,2
). Figure 3 shows the PRCCs results which illustrate the dependence of R
0
i
(
i
=
1,2
) on each parameter, and we considered |
PRCC
|
∈
(
0.4
,
∞
)
,
(
0.2,0.4
)
,
(
0,0.2
) as indicating correlations of very important, moderate, or little significance between input parameters and output variables, respectively. The positive sign of their PRCCs indicates that if the parameters are increased, the value of R
0
2 increases and vice versa for negative signs.Figure 3
Sensitivity analysis of thresholdsR
0
1 and R
0
2: the baseline parameter values are a
=
0.8, b
=
0.5, c
=
1.6, ω
=
0.9, μ
=
0.4, d
=
0.3, σ
=
0.5, p
1
max
⁡
=
0.95, p
2
max
⁡
=
0.2, θ
1
=
0.7, θ
2
=
2, and T
=
1. (a) PRCC results with R
0
1; (b) PRCC results with R
0
2.By using LHS with 3,000 samples, uncertainty and sensitivity analyses for all parameters in model (3) were determined. A uniform distribution function was used and tested for significant PRCCs for all parameters with wide ranges, making use of the data from Table 1, and the baseline values of all parameters are given in the figure legend of Figure 3. The PRCCs sensitivity analysis on R
0
1 and R
0
2 between all parameters is depicted in Figures 3(a) and 3(b), which shows that variations in the intrinsic growth rate a, the death rate c, the pulse period T, and the constant stocking rate σ have the four greatest effects on the outbreak of the pest with the first three increasing the R
0
i
(
i
=
1,2
) and the latter decreasing them. The half saturation constant θ
2 generates minor decreasing impact on two thresholds with similar amounts of influence in both figures. The maximum fatality rate p
2
max
⁡ has a minor increasing effect on R
0
i
(
i
=
1,2
). However, the other maximum fatality rate p
1
max
⁡ has a great decreasing impact on R
0
2; that is, p
1
max
⁡ plays an important role in the global stability, while it is without influence upon the local stability of (
x
p
(
t
)
,
y
p
(
t
)
) of model (3).Table 1
The main coefficient values of model (3). All other parameters are fixed as follows: b
=
0.5, c
=
1.6, ω
=
0.9, and μ
=
0.4.
Parameter
a
d
p
1
max
⁡
p
2
max
⁡
θ
2
σ
T
Definition
Intrinsic growth rate
Death rate
Maximum fatality rate
Maximum fatality rate
Half saturation constant
Constant stocking number
Pulse period
Range
0.1~10
0.01~1
0.01~1
0.01~1
0.01~10
0.01~10
0.5~20So far, the dynamics of the pest-free periodic solution (12) of model (3) with limited resource has been investigated completely. The results demonstrate that the pest population can be eradicated, provided that certain conditions are satisfied.


## 3.1. Existence of the Pest-Free Periodic Solution

In order to illustrate the existence of the pest-free periodic solution, we first focus on the pest-free set:{
(
x
,
y
)
∈
ℛ
+
2
,
x
=
0
}, which is clearly invariant by model (3); within this set, model (3) becomes
(4)
y
˙
=
-
d
y
(
t
)
,
t
≠
n
T
,
y
(
t
+
)
=
(
1
-
p
2
max
⁡
y
(
t
)
y
(
t
)
+
θ
2
)
y
(
t
)
+
σ
,
t
=
n
T
,
y
(
0
+
)
=
y
0
.
Solving the above equation in interval (
n
T
,
(
n
+
1
)
T
], one yields
(5)
y
(
(
n
+
1
)
T
+
)
=
(
1
-
p
2
max
⁡
y
(
n
T
+
)
e
-
d
T
y
(
n
T
+
)
e
-
d
T
+
θ
2
)
y
(
n
T
+
)
e
-
d
T
+
σ
.
Denote y
(
n
T
+
)
=
y
n; then the above equation can be rewritten as the following difference equation:
(6)
y
n
+
1
=
(
1
-
p
2
max
⁡
)
e
-
2
d
T
y
n
2
+
θ
2
e
-
d
T
y
n
e
-
d
T
y
n
+
θ
2
+
σ
≐
F
(
y
n
)
,
which is the so-called stroboscopic map of model (4), and describes the relations of the number of natural enemies between any two successive pulse points. Consequently, the existence of the positive steady state of model (6) implies the existence of a positive periodic solution of model (4). Therefore, we first discuss the existence of a positive steady state of (6); that is, the derivative of F
(
y
n
) with respect to y
n yields
(7)
F
′
(
y
n
)
=
[
(
1
-
p
2
max
⁡
)
+
p
2
max
⁡
θ
2
2
(
e
-
d
T
y
n
+
θ
2
)
2
]
e
-
d
T
and it is easy to see that 0
<
F
′
(
y
n
)
<
1 holds true for all y
n.Next, we discuss the positive fixed point of the stroboscopic map (6), denoted by y
*, which satisfies the following equation:
(8)
A
1
y
*
2
+
B
1
y
*
+
C
1
=
0
,
where A
1
=
exp
⁡
(
-
d
T
)
[
1
-
(
1
-
p
2
max
⁡
)
exp
⁡
(
-
d
T
)
]
>
0, B
1
=
θ
2
-
σ
exp
⁡
(
-
d
T
)
-
θ
2
exp
⁡
(
-
d
T
), C
1
=
-
θ
2
σ
<
0.DenoteΔ
=
B
1
2
-
4
A
1
C
1, and by simple calculation, we have
(9)
Δ
=
(
θ
2
-
σ
e
-
d
T
-
θ
2
e
-
d
T
)
2
+
4
θ
2
σ
e
-
d
T
[
1
-
(
1
-
p
2
max
⁡
)
e
-
d
T
]
>
0
.
Obviously, (6) has a unique positive root; that is,
(10)
y
*
=
-
B
1
+
B
1
2
-
4
A
1
C
1
2
A
1
,
which is stable due to 0
<
F
′
(
y
*
)
<
1.According to the relations between a fixed point of the stroboscopic map (6) and the periodic solution of model (4), we conclude that the submodel (4) of (3) has a unique nontrivial positive periodic solution, denoted by y
p
(
t
) and
(11)
y
p
(
t
)
=
y
*
exp
⁡
[
-
d
(
t
-
n
T
)
]
for
n
T
<
t
≤
(
n
+
1
)
T
,
and it follows from 0
<
F
′
(
y
n
)
<
1 and [32] that we have the following theorem.Theorem 1.
Model (4) has a positive periodic solution y
p
(
t
) with periodic T and for every solution y
(
t
) of (4) such that |
y
(
t
)
-
y
p
(
t
)
|
→
0 as t
→
+
∞, where y
p
(
t
)
=
y
*
e
-
d
(
t
-
n
T
) and y
* is defined in (10).By using Theorem1, we conclude that the unique fixed point y
* of (6) is globally stable. Accordingly, the stability of the positive periodic solution is determined by the stability of the positive equilibria of the difference equation made by the sequence of impulsive points. So, the positive periodic solution y
p
(
t
) of model (4) is globally stable.Therefore, we obtain the general expression of an unique pest-free periodic solution of model (3) over the interval n
T
<
t
≤
(
n
+
1
)
T for all n
∈
𝒩 denoted by
(12)
(
x
p
(
t
)
,
y
p
(
t
)
)
=
(
0
,
y
*
e
-
d
(
t
-
n
T
)
)
.


## 3.2. Stability of the Pest-Free Periodic Solution

The threshold conditions which guarantee the global stability of the pest-free periodic solution (12) play a key role in pest control. To show this, we now concentrate on the stability of the pest-free solution (12) for model (3). We change model (3) variables to consider the deviations from the pest-free solution that are denoted by (
x
¯
(
t
)
,
y
¯
(
t
)
) so that
(13)
(
x
¯
(
t
)
,
y
¯
(
t
)
)
=
(
x
(
t
)
,
y
(
t
)
)
-
(
x
p
(
t
)
,
y
p
(
t
)
)
that yields
(14)
x
¯
˙
=
a
x
¯
-
b
x
¯
2
-
c
x
¯
1
+
ω
x
¯
[
y
¯
+
y
p
(
t
)
]
,
y
¯
˙
=
μ
c
x
¯
1
+
ω
x
¯
[
y
¯
+
y
p
(
t
)
]
-
d
y
¯
.
Assuming that (
x
¯
,
y
¯
) are small enough, we get the linear approximation of the deviation system around the periodic solution (
x
p
(
t
)
,
y
p
(
t
)
),
(15)
x
¯
˙
=
(
a
-
c
y
p
(
t
)
)
x
¯
,
y
¯
˙
=
μ
c
y
p
(
t
)
x
¯
-
d
y
¯
.Before stating Theorem2, for the sake of convenience, we set
(16)
R
0
1
=
a
d
T
c
y
*
(
1
-
e
-
d
T
)
,
R
0
2
=
c
d
T
M
s
c
y
*
(
1
-
e
-
d
T
)
-
d
ln
⁡
(
1
-
p
1
max
⁡
)
,
where M
s
=
(
a
ω
+
b
)
2
/
(
4
b
c
ω
).We are now in a position to present the sufficient conditions for the global stability of the pest-free periodic solution(
x
p
(
t
)
,
y
p
(
t
)
) of model (3) as follows.Theorem 2.
The pest-free periodic solution(
x
p
(
t
)
,
y
p
(
t
)
) of model (3) is locally stable in the first quadrant, provided that R
0
1
≤
1 and is globally stable if max
⁡
{
R
0
1
,
R
0
2
}
≤
1.The proof of Theorem2 is given in Appendix A.Remark 3.
According to the formulaeR
0
1 and R
0
2, we have
(17)
R
0
1
-
R
0
2
=
d
T
[
y
*
(
1
-
e
-
d
T
)
(
a
2
-
c
2
M
s
)
-
a
d
ln
⁡
(
1
-
p
1
max
⁡
)
]
c
y
*
(
1
-
e
-
d
T
)
[
a
y
*
(
1
-
e
-
d
T
)
-
d
ln
⁡
(
1
-
p
1
max
⁡
)
]
.
Therefore, if y
*
(
1
-
e
-
d
T
)
(
a
2
-
c
2
M
s
)
-
a
d
ln
⁡
(
1
-
p
1
max
⁡
)
>
0, then R
0
1
>
R
0
2, and consequently the pest-free periodic solution (
x
p
(
t
)
,
y
p
(
t
)
) of model (3) is globally stable if R
0
1
<
1, which shows that the local stability indicates the global stability. If y
*
(
1
-
e
-
d
T
)
(
a
2
-
c
2
M
s
)
-
a
d
ln
⁡
(
1
-
p
1
max
⁡
)
<
0, then R
0
1
<
R
0
2 (see Figure 1(a)). Thus, the pest-free periodic solution (
x
p
(
t
)
,
y
p
(
t
)
) of model (3) is globally stable if R
0
2
<
1 (see Figure 1(b)). If so, we cannot show that in this case the local stability indicates the global stability theoretically.The monotonicity ofR
0
1 and R
0
2 with respect to T. (a) R
0
1
>
R
2. All parameters are fixed as follows: a
=
0.5, b
=
0.8, c
=
1.6, ω
=
1, μ
=
0.4, d
=
0.37, σ
=
0.78, p
1
max
⁡
=
0.95, p
2
max
⁡
=
0.2, θ
1
=
1, θ
2
=
2. (b) R
0
1
<
R
2. All parameters are fixed as follows: a
=
0.1, b
=
0.8, c
=
2.5, ω
=
1, μ
=
0.4, d
=
0.37, σ
=
1, p
1
max
⁡
=
0.95, p
2
max
⁡
=
0.1, θ
1
=
1, and θ
2
=
2.
(a)
(b)


## 3.3. Threshold Conditions

Note that although the threshold valuesR
0
i depend on all parameters of model (3), the most interesting parameters here are the pulse period T and some parameters related to resource limitation such as the maximum fatality rates p
i
max
⁡
(
i
=
1,2
) and the half saturation constant θ
2, so it makes sense to know how T, p
i
max
⁡, θ
2 affect the threshold conditions (R
0
i
<
1; i
=
1,2) which guarantee the global stability of the pest-free periodic solution (12). Figure 1(a) indicates that R
0
i
(
i
=
1,2
) are monotonic increasing functions with respect to T
∈
(
0,20
]. Moreover, there are some critical values of T such that R
0
i
=
1, which indicate that there exists one maximum allowable impulsive period T
i
max
⁡ such that R
0
i
<
1 for all T
<
T
i
max
⁡ and i
=
1,2, respectively. Meanwhile, R
0
2 takes more time to reach the threshold value 1 than R
0
1 as a result of T
1
max
⁡
<
T
2
max
⁡. In more detail, the contour plots of Figures 2(a) and 2(c) show that the control strategies (decreasing p
2
max
⁡, increasing θ
2) can reduce the threshold value R
0
i
(
i
=
1,2
). However, increasing p
1
max
⁡ can reduce the threshold value R
0
2 from Figure 2(d). Therefore, in order to control the pest population, it is very important to increase the maximum fatality rate p
1
max
⁡ for the pest while decreasing the maximum fatality rate p
2
max
⁡ for the natural enemy, which means that it is essential to increase the pesticide’s efficacy and reduce its effectiveness against the natural enemy while enhancing the efficiency of the natural enemies. Also, it is important to carry out pulsed control actions more and more frequently under resource limitation according to Figure 2(b), and then we can eradicate the pest successfully.Contour plots ofR
0
1 and R
0
2: (a) plot of contours of R
0
1 versus p
2
max
⁡ and θ
2; (b) plot of contours of R
0
2 versus T and p
1
max
⁡; (c) plot of contours of R
0
2 versus p
2
max
⁡ and θ
2; (d) plot of contours of R
0
1 versus p
2
max
⁡ and p
1
max
⁡. All other parameters are fixed as follows: a
=
0.5, b
=
0.8, c
=
1.6, ω
=
1, μ
=
0.4, d
=
0.37, σ
=
0.78, p
1
max
⁡
=
0.95, p
2
max
⁡
=
0.2, θ
1
=
1, θ
2
=
2, T
=
6.5
(
A
), and T
=
12
(
C
-
D
).
(a)
(b)
(c)
(d)


## 3.4. Sensitivity to Variations

Previous analysis indicates thatR
0
i
(
i
=
1,2
) are significant threshold conditions affecting outbreaks of the pest population. Sensitivity analysis of the most significant parameters (such as intrinsic growth rate a, death rate c, the maximum fatality rate p
i
max
⁡, half saturation constant θ
2, constant stocking number σ, and pulse period T) was performed by evaluating the PRCCs for various input parameters against the threshold condition R
0
i
(
i
=
1,2
) with the LHS method (see [24, 25, 33]). PRCCs show how influential each parameter is on R
0
i
(
i
=
1,2
). Figure 3 shows the PRCCs results which illustrate the dependence of R
0
i
(
i
=
1,2
) on each parameter, and we considered |
PRCC
|
∈
(
0.4
,
∞
)
,
(
0.2,0.4
)
,
(
0,0.2
) as indicating correlations of very important, moderate, or little significance between input parameters and output variables, respectively. The positive sign of their PRCCs indicates that if the parameters are increased, the value of R
0
2 increases and vice versa for negative signs.Figure 3
Sensitivity analysis of thresholdsR
0
1 and R
0
2: the baseline parameter values are a
=
0.8, b
=
0.5, c
=
1.6, ω
=
0.9, μ
=
0.4, d
=
0.3, σ
=
0.5, p
1
max
⁡
=
0.95, p
2
max
⁡
=
0.2, θ
1
=
0.7, θ
2
=
2, and T
=
1. (a) PRCC results with R
0
1; (b) PRCC results with R
0
2.By using LHS with 3,000 samples, uncertainty and sensitivity analyses for all parameters in model (3) were determined. A uniform distribution function was used and tested for significant PRCCs for all parameters with wide ranges, making use of the data from Table 1, and the baseline values of all parameters are given in the figure legend of Figure 3. The PRCCs sensitivity analysis on R
0
1 and R
0
2 between all parameters is depicted in Figures 3(a) and 3(b), which shows that variations in the intrinsic growth rate a, the death rate c, the pulse period T, and the constant stocking rate σ have the four greatest effects on the outbreak of the pest with the first three increasing the R
0
i
(
i
=
1,2
) and the latter decreasing them. The half saturation constant θ
2 generates minor decreasing impact on two thresholds with similar amounts of influence in both figures. The maximum fatality rate p
2
max
⁡ has a minor increasing effect on R
0
i
(
i
=
1,2
). However, the other maximum fatality rate p
1
max
⁡ has a great decreasing impact on R
0
2; that is, p
1
max
⁡ plays an important role in the global stability, while it is without influence upon the local stability of (
x
p
(
t
)
,
y
p
(
t
)
) of model (3).Table 1
The main coefficient values of model (3). All other parameters are fixed as follows: b
=
0.5, c
=
1.6, ω
=
0.9, and μ
=
0.4.
Parameter
a
d
p
1
max
⁡
p
2
max
⁡
θ
2
σ
T
Definition
Intrinsic growth rate
Death rate
Maximum fatality rate
Maximum fatality rate
Half saturation constant
Constant stocking number
Pulse period
Range
0.1~10
0.01~1
0.01~1
0.01~1
0.01~10
0.01~10
0.5~20So far, the dynamics of the pest-free periodic solution (12) of model (3) with limited resource has been investigated completely. The results demonstrate that the pest population can be eradicated, provided that certain conditions are satisfied.


## 4. The Existence of the Positive Periodic Solution

In this section, we adopt bifurcation theory [34] to analyze the existence of a positive periodic solution of model (3) near the pest-free periodic solution by setting the impulsive period T and constant stocking number σ as bifurcation parameters. According to the bifurcation theorem in [34], we first obtain the following result.Theorem 4.
The supercritical branch occurs at the pointT
0
>
0 satisfying R
0
1
(
T
0
)
=
1; that is, it will have a stable positive periodic solution when T
>
T
0 and is close to T
0, provided one of the following conditions of model (3) holds(C1)
M
1
<
0 and 1
-
A
3
A
4
<
0;
(C2)
M
1
>
0 and N
1
+
N
2
<
0,where
(18)
M
1
=
-
c
d
T
0
y
*
(
T
0
)
·
A
3
e
-
d
T
0
1
-
A
3
e
-
d
T
0
-
[
a
-
c
y
p
(
T
0
)
]
;
N
1
=
2
p
1
max
⁡
θ
1
-
1
·
e
2
∫
0
T
0
‍
[
a
-
c
y
p
(
ξ
)
]
d
ξ
-
A
5
;
N
2
=
-
2
c
T
0
·
1
-
A
3
A
4
1
-
A
3
e
-
d
T
0
·
e
∫
0
T
0
‍
[
a
-
c
y
p
(
ξ
)
]
d
ξ
;
A
3
=
1
-
p
2
max
⁡
+
p
2
max
⁡
θ
2
2
[
y
p
(
T
0
)
+
θ
2
]
2
;
A
4
=
μ
c
∫
0
T
0
‍
y
p
(
v
)
·
e
-
d
(
T
0
-
v
)
·
e
∫
0
v
‍
[
a
-
c
y
p
(
ξ
)
]
d
ξ
d
v
;
A
5
=
-
2
T
0
·
[
b
+
c
ω
y
p
(
v
)
]
·
e
∫
0
T
0
‍
[
a
-
c
y
p
(
ξ
)
]
d
ξ
-
c
∫
0
T
0
‍
{
e
∫
v
T
0
‍
[
a
-
c
y
p
(
ξ
)
]
d
ξ
·
e
∫
0
θ
‍
[
a
-
c
y
p
(
ξ
)
]
d
ξ
}
d
θ
·
∫
0
v
‍
{
μ
c
y
p
(
θ
)
·
e
-
d
(
v
-
θ
)
e
∫
0
θ
‍
[
a
-
c
y
p
(
ξ
)
]
d
ξ
·
e
∫
0
θ
‍
[
a
-
c
y
p
(
ξ
)
]
d
ξ
}
d
θ
}
d
v
.The proof of Theorem4 is given in Appendix B. Theorem 4 shows that there exists a positive periodic solution under some conditions, providing that the pest-free periodic solution becomes unstable. It follows from Theorem 4 that if T
>
T
0 and is closed to T
0, the periodic solution of model (3) is stable.Figure4(a) indicated that if the pulse period T is more than T
1
max
⁡, the pest-free solution becomes unstable and variable x begins to oscillate with a large amplitude that corresponds to periodic outbreaks of the pest population. If the pulse period is further increased, a sequence of period adding bifurcations interchanging with regions of chaos is observed. Figures 4 and 5 show that model (3) has more complex and interesting dynamic behaviors including periodic doubling bifurcation, chaotic solutions, quasiperiodic solutions, tangent bifurcation, multistability, chaos crisis, periodic windows, periodic halving bifurcation, and crises with the increasing of pulse period T and constant stocking number σ, respectively. For example, Figures 4(a) and 4(b) show that for model (3) there exists a quasiperiodic solution when T
=
8; see Figure 6 for full details.Bifurcation diagrams of model (3) with respect to bifurcation parameter T, where all parameter values are fixed as follows: a
=
8, b
=
5
,
c
=
1
,
ω
=
1
,
μ
=
0.95
,
d
=
0.2
,
σ
=
0.8
,
p
1
max
⁡
=
0.5
,
p
2
max
⁡
=
0.01
,
θ
1
=
0.2
,
θ
2
=
2, and (
x
0
,
y
0
)
=
(
2,1
). (a) Prey population x with T; (b) predator population y with T.
(a)
(b)Bifurcation diagrams of model (3) with respect to bifurcation parameter σ. (a) Prey population x with σ; (b) predator population y with σ. The other parameters are identical to those in Figure 4 and T
=
0.5.
(a)
(b)Figure 6
A typical quasiperiodic solution of model (3) over t
∈
[
17000,26000
]. (a) Phase portrait of x and y with t; (b) time series of x; (c) time series of y. The other parameters are identical to those in Figure 4 and T
=
8.Meanwhile, bifurcation analyses also indicate that multiple attractors can coexist for a wide range of parameters, for example, the two attractors with quite different pest amplitudes that can coexist whenσ
=
7.3 in Figure 5; see Figure 7 for more details.Two coexisting attractors of model (3) with σ
=
7.3. The initial values from left to right are (
x
0
,
y
0
)
=
(
2.5,1
)
,
(
2,5.1
), respectively. (a-b) Periodic attractor with period 12; (c-d) periodic attractor with period 16. The other parameters are identical to those in Figure 4 and T
=
6.
(a)
(b)
(c)
(d)Furthermore, we can clearly see that there are two examples of attractor coexistence of model (3) when σ
=
4.37 in Figure 5; when small changes are introduced in initial value, numerical simulations imply that this solution can switch to other attractors, that is, the switch-like transitions between the two attractors shown in Figure 8 with (
x
0
,
y
0
)
=
(
2,5
).Attractors switch-like behavior of model (3) with σ
=
4.37 and (
x
0
,
y
0
)
=
(
2,5
). (a) Time series of x; (b) time series of y. The other parameters are identical to those in Figure 7.
(a)
(b)In summary, numerical simulations show that with resource limitation model (3) has more and more complex dynamics. It is well known that the dosages of pesticide applied and numbers of natural enemies released are crucial for controlling pests, but when these key factors of controlling pests are limited, it makes it more difficult for pest control.


## 5. Discussion and Biological Conclusions

The main purpose of this paper was to understand the effect of resource limitation on outbreaks of a pest population, so we conducted a dynamical analysis of a Holling II predator-prey model which incorporates a nonlinear pulse. To this end, we introduced a nonlinear formp
i
(
ϕ
)
=
p
i
max
⁡
ϕ
/
(
ϕ
+
θ
i
) (i
=
1, ϕ
=
x; or i
=
2, ϕ
=
y) as fatality rate for pest or natural enemy. The nonlinear fatality rate taking density dependence into consideration gives a better portrayal of the reality, which can be used to obtain a more accurate evaluation of the pest-natural enemy system. However, it is worth noting that the mathematical model described here differs from the ones discussed by many other researchers since we have taken the nonlinear impulsive function into account. The nonlinear nature of the impulsive functions makes the dynamical behavior of solutions change dramatically and it gets more and more complicated, which means that it is a very difficult task to control pests under resource limitation. We discuss the control strategy for the pest through theoretical analysis and numerical simulations.The results indicate that pesticides and natural enemies as two approaches to control pests must be committed to a long-term strategy especially when being under resource limitation. Stability analysis for model (3) shows that R
0
i could be a global threshold in a sense that if all these control methods are adopted so that R
0
i
<
1
(
i
=
1,2
), then the strategies are effective enough to eradicate the pest population. We adopt LHS/PRCCs uncertainty and sensitivity analysis techniques to gain a better understanding of the impacts of parameters on the threshold when various parameters are changed within the ranges of values observed empirically (see Table 1 and Figure 3). The results illustrate that we should pay more attention to the impulsive period T, maximum fatality rate p
1
max
⁡ for the pest, and constant stocking number σ which maximizes the threshold values R
0
i to prevent outbreaks of the pest population. That is, it is essential to increase the pesticide efficacy and the efficiency of the natural enemy. Furthermore, when resources are limited, it is important to carry out pesticide and natural enemy programs more frequently than when resources are not limited. Our results also confirm that IPM is the optimal control strategy.Based on Theorem4, the results demonstrate that the dynamic behavior of model (3) may be dramatically affected by small changes in the value of initial densities of pest and natural enemy with resource limitation. By choosing impulsive perturbations σ and T as bifurcation parameters, we have obtained bifurcation diagrams (see Figures 4 and 5) including periodic doubling bifurcation, chaotic solutions, quasiperiodic solutions (see Figure 6 for details), tangent bifurcation, multistability (see Figure 7 for details), chaos crisis, periodic windows, periodic halving bifurcation, and crises, which can help us to further understand the application of nonlinear pulses in our model (for more information about bifurcation diagrams, refer to [15, 18, 35]). According to bifurcation diagrams, we can find that the routes to chaos are very complicated; that is, with resource limitation, there are several hidden factors that can adversely affect our control strategy. The increasing number of potential complexities predicted by the theory does not seem to make this task any easier. Nevertheless, identifying complicated, possibly chaotic, dynamics in IPM models may present a major challenge for controlling the pest populations in practice.In our study, we have investigated the dynamical behavior of a Holling II predator-prey model to understand how the limited resource affects pest outbreaks. To link the costs of developing and implementing controls to population dynamic modeling of disease epidemics, it is necessary to consider other resource limitation factors such as delayed responses and residual effects on pests of pesticides, as such nonlinear factors can affect the success of pest control. We leave these topics for further investigations.



---

*Source: 450935-2014-03-11.xml*