# Chemical Control for Host-Parasitoid Model within the Parasitism Season and Its Complex Dynamics

Authors: Tao Wang; Youtang Zhang
Journal: Discrete Dynamics in Nature and Society
 (2016)

---

## Abstract

In the present paper, we develop a host-parasitoid model with Holling type II functional response function and chemical control, which can be applied at any time of each parasitism season or pest generation, and focus on addressing the importance of the timing of application pesticide during the parasitism season or pest generation in successful pest control. Firstly, the existence and stability of both the host and parasitoid populations extinction equilibrium and parasitoid-free equilibrium have been investigated. Secondly, the effects of key parameters on the threshold conditions have been discussed in more detail, which shows the importance of pesticide application times on the pest control. Thirdly, the complex dynamics including multiple attractors coexistence, chaotic behavior, and initial sensitivity have been studied by using numerical bifurcation analyses. Finally, the uncertainty and sensitivity of all the parameters on the solutions of both the host and parasitoid populations are investigated, which can help us to determine the key parameters in designing the pest control strategy. The present research can help us to further understand the importance of timings of pesticide application in the pest control and to improve the classical chemical control and to make management decisions.


---

## Body

## 1. Introduction

Integrated pest management (IPM) is the selection and application of pest control actions that can ensure favourable economic and ecological consequences [1–3], which has been successfully applied to agricultural pest management situations. IPM employs a variety of tactics including cultural controls, biological controls, and chemical controls or pesticides, while biological pest control is one of the most common measures applied in IPM through the control and management of natural predators and parasites. For example, mosquitoes are often controlled by puttingBacillus thuringiensis (Bt) ssp.israelensis, a bacterium that infects and kills mosquito larvae, in local water sources. The aim of biological pest control is to eliminate a pest which can minimize harm to the ecological balance of the environment in its present form.Pesticide application or chemical control is another important component of IPM measures, which refers to the practical way in which pesticides including herbicides, fungicides, and insecticides are sprayed to the pest population. Public concern about the application of pesticides has highlighted the need to make this process as efficiently as possible which could minimize their release into the environment and human exposure [1–3]. Although pesticide and insecticide applications can cause a number of problems, in most cropping systems they are still the principle means of controlling pests once the economic threshold that defines the lowest population density at which the control actions should be applied has been reached. In practice, pesticides can be relatively cheap, are easy to apply and fast-acting, and in most instances can be relied on to control the pests [1–3].Undoubtedly, mathematical model is one of key tools to help us to understand those factors in the IPM strategies and the pest control. Recently, the continuous predator-prey models concerning IPM strategy have been developed and investigated [1, 4–10]. In particular, several factors including pest natural enemy ratios, starting densities, timings of natural enemy releases, dosages and timings of insecticide applications, and instantaneous killing rates of pesticides on both pests and natural enemies have been addressed by Tang et al. [10]. Moreover, the importance of timings of pesticide sprays and natural enemy releases has been studied through the stability threshold condition for a pest eradication periodic solution.However, the discrete or nonoverlapping generation is a common feature among the host and parasitoid populations [11–15]. If so the use of continuous-time models to describe the interaction between the pest population and its natural enemy population becomes questionable. Those show that the discrete host-parasitoid models are much more realistic when the populations have discrete and synchronized generations [16–21]. The complex dynamics of the host-parasitoid model with Holling II functional response function have been investigated by Tang and Chen [19]. The classical Nicholson-Bailey model for a two species’ host-parasitoid system with discrete generations and IPM strategies has been studied by Tang et al. [7]. Note that the pulse IPM control strategies assumed to be applied at the end of each periodic number of generations and some important issues concerning IPM have been addressed.According to the facts and main results of [7, 10], we know that the timings of pesticide application and releasing natural enemies play a key role in successful pest control. Note that there are several different pesticide spraying methods based on the pest growth generations and the parasitism season: (1) the pesticides are applied at the beginning of each generation; (2) the pesticides are sprayed at the end of each generation; and (3) the pesticides can be applied at any time point within the pest growth generation and the parasitism season. Therefore, the questions are whether there exists an optimal pesticide application time at which the number of the pest populations can be minimized, even eradicated, and how the pesticide application time and efficiency affect the successful pest control and complex outbreak patterns.To address those questions, it follows from the basic modelling methods proposed in literatures [22–24] that we assume that the chemical control tactic is applied instantaneously within the generation at any time. The main purpose is to extend the Holling II host-parasitoid model with a fraction of survival rate of parasitoid from one generation to next [25, 26] and then investigate the effects of timings of pesticide application and killing rate on this model. In particular, the existence and stability of both the host and parasitoid populations extinction equilibrium and parasitoid-free equilibrium have been investigated. The interesting results indicate that the different pesticide application times could result in significantly different size of the host population and consequently influence the pest control. Moreover, the effects of all important parameters including timings of pesticide application on the threshold conditions have been discussed in more detail. Further, the complex dynamics including multiple attractors coexistence, quasi-periodic windows, chaotic behavior, and initial sensitivity have been studied by using numerical bifurcation analyses. Finally, the uncertainty and sensitivity of all the parameters on the solutions of both the host and parasitoid populations have been investigated, which can help us to determine the key parameters in designing the successful pest control strategy. The present research can help us to further understand the importance of pesticide application times in the pest control and to improve the classical chemical control and to make management decisions.


## 2. The Host-Parasitoid Model with Chemical Control

### 2.1. Host Growth Model with Chemical Control

The based discrete map for host or pest population can be described as follows:(1)Hn+1=fHn=HnFHn,where Hn∈[0,∞) is the pest population size at generation n∈N and f,F:[0,∞)→R are, respectively, the population production and the per capita production. A common example is the Rocker model fH=Hexp⁡(r(1-H/K)) with r>0 and K>0.It follows from the basic modelling methods proposed in literatures [22–24] that we assume that the chemical control tactic is applied instantaneously within the generation [n,n+1]; that is, there exists a positive constant θ with 0≤θ≤1 such that the pesticide is sprayed at time point n+θ. Further, we assume that a proportional number (denoted by q) of the pest populations has been killed, and thus after n+θ, the growth of pest population and its production depends on (1-q)Hn. Therefore, model (1) with chemical control can be formulated as(2)Hn+1=1-qHnθFHn+1-θF1-qHn.According tof(Hn)=HnF(Hn), we have(3)Hn+1=θ1-qfHn+1-θf1-qHn.


### 2.2. Host-Parasitoid Model with the Effect of Overwintering Parasitoid and Chemical Control

Involving the parasitoid population into (1), the host-parasitoid system with discrete generations can be written [16, 17] as the generalized model(4)Hn+1=1-qHnθfHn+1-θf1-qHn·gθ+1-θ1-qHn,Pn,Pn+1=γ1-qHn1-gθ+1-θ1-qHn,Pn+δPn,where Hn and Pn represent the host and parasitoid population abundance in generation n, respectively, f(H) denotes the per capita net rate of the increase of the host population in the absence of parasitoid population, and g(H,P) is the proportion of host individuals that escape attack by the parasitoid. Note that P is generally interpreted as the abundance of adult parasitoid females and N as the abundance of host adults and thus γ represents the number of parasitoid eggs laid per host, the survival of the parasitoid in the attacked hosts, and the sex ratio of the emerging parasitoid adults, where δ∈[0,1) denotes the density-independent survival of parasitoid propagules at generation n. A fraction of pathogen can survive from one generation to the next. This discrete-generation framework characterizes a perfectly synchronized parasitoid interacting with a host that has distinct generations, which is frequent in host-parasitoid systems in temperate regions of the world and even some from more tropical regions when parasitism causes generation cycles within the overlapping generations of some multivoltine hosts [16, 17].In particular, ifθ=0 then model (4) can be changed as the following impulsive difference equation:(5)Hn+1=HnfHn·gHn,PnPn+1=γHn1-gHn,Pn+δPn,n=0,1,2,…,Hn+=1-qHn,Pn+=Pn,n=0,1,2,…,which indicates that the chemical control tactics have been applied at the beginning of each generation or parasitism season. For more details of impulsive difference equations, please see [7].Ifθ=1 then model (4) can be changed as the following impulsive difference equation:(6)Hn+1=HnfHn·gHn,PnPn+1=γHn1-gHn,Pn+δPn,n=0,1,2,…,Hn+1+=1-qHn+1,Pn+1+=Pn+1,n=0,1,2,…,which shows that the chemical control tactics have been applied at the ending of each generation or parasitism season.Holling (1959) was the first ecologist to investigate the interaction between host and parasitoid populations and discussed the functional relationship in depth through the experiments of the shrews and deer mice feeding on sawfly cocoons. Holling (1959) presents three different functional response classes, for which the following are nominal forms: type I (linear then constant), type II (decelerating rise to an upper asymptote), and type III (sigmoidal) [16, 17, 25, 26]. In the present work, we assume that the pest population follows the Ricker model in the absence of parasitoid population and the Holling II functional response function for parasitism; that is, we have(7)fHn=exp⁡r1-HnK,gHn,Pn=exp⁡-αTPn1+αThHn.Therefore, based on the above special choices we will focus on the existence and stability of the boundary equilibrium which concerns the outbreak of the host population and then address how the timing of spraying pesticides affects the values of this boundary equilibrium. Moreover, the completely numerical bifurcation analyses and sensitivity analyses have been carried out to show the effects of interesting parameters including instant killing rateq and pesticide spraying time point θ on the successful pest control and outbreak.


## 2.1. Host Growth Model with Chemical Control

The based discrete map for host or pest population can be described as follows:(1)Hn+1=fHn=HnFHn,where Hn∈[0,∞) is the pest population size at generation n∈N and f,F:[0,∞)→R are, respectively, the population production and the per capita production. A common example is the Rocker model fH=Hexp⁡(r(1-H/K)) with r>0 and K>0.It follows from the basic modelling methods proposed in literatures [22–24] that we assume that the chemical control tactic is applied instantaneously within the generation [n,n+1]; that is, there exists a positive constant θ with 0≤θ≤1 such that the pesticide is sprayed at time point n+θ. Further, we assume that a proportional number (denoted by q) of the pest populations has been killed, and thus after n+θ, the growth of pest population and its production depends on (1-q)Hn. Therefore, model (1) with chemical control can be formulated as(2)Hn+1=1-qHnθFHn+1-θF1-qHn.According tof(Hn)=HnF(Hn), we have(3)Hn+1=θ1-qfHn+1-θf1-qHn.


## 2.2. Host-Parasitoid Model with the Effect of Overwintering Parasitoid and Chemical Control

Involving the parasitoid population into (1), the host-parasitoid system with discrete generations can be written [16, 17] as the generalized model(4)Hn+1=1-qHnθfHn+1-θf1-qHn·gθ+1-θ1-qHn,Pn,Pn+1=γ1-qHn1-gθ+1-θ1-qHn,Pn+δPn,where Hn and Pn represent the host and parasitoid population abundance in generation n, respectively, f(H) denotes the per capita net rate of the increase of the host population in the absence of parasitoid population, and g(H,P) is the proportion of host individuals that escape attack by the parasitoid. Note that P is generally interpreted as the abundance of adult parasitoid females and N as the abundance of host adults and thus γ represents the number of parasitoid eggs laid per host, the survival of the parasitoid in the attacked hosts, and the sex ratio of the emerging parasitoid adults, where δ∈[0,1) denotes the density-independent survival of parasitoid propagules at generation n. A fraction of pathogen can survive from one generation to the next. This discrete-generation framework characterizes a perfectly synchronized parasitoid interacting with a host that has distinct generations, which is frequent in host-parasitoid systems in temperate regions of the world and even some from more tropical regions when parasitism causes generation cycles within the overlapping generations of some multivoltine hosts [16, 17].In particular, ifθ=0 then model (4) can be changed as the following impulsive difference equation:(5)Hn+1=HnfHn·gHn,PnPn+1=γHn1-gHn,Pn+δPn,n=0,1,2,…,Hn+=1-qHn,Pn+=Pn,n=0,1,2,…,which indicates that the chemical control tactics have been applied at the beginning of each generation or parasitism season. For more details of impulsive difference equations, please see [7].Ifθ=1 then model (4) can be changed as the following impulsive difference equation:(6)Hn+1=HnfHn·gHn,PnPn+1=γHn1-gHn,Pn+δPn,n=0,1,2,…,Hn+1+=1-qHn+1,Pn+1+=Pn+1,n=0,1,2,…,which shows that the chemical control tactics have been applied at the ending of each generation or parasitism season.Holling (1959) was the first ecologist to investigate the interaction between host and parasitoid populations and discussed the functional relationship in depth through the experiments of the shrews and deer mice feeding on sawfly cocoons. Holling (1959) presents three different functional response classes, for which the following are nominal forms: type I (linear then constant), type II (decelerating rise to an upper asymptote), and type III (sigmoidal) [16, 17, 25, 26]. In the present work, we assume that the pest population follows the Ricker model in the absence of parasitoid population and the Holling II functional response function for parasitism; that is, we have(7)fHn=exp⁡r1-HnK,gHn,Pn=exp⁡-αTPn1+αThHn.Therefore, based on the above special choices we will focus on the existence and stability of the boundary equilibrium which concerns the outbreak of the host population and then address how the timing of spraying pesticides affects the values of this boundary equilibrium. Moreover, the completely numerical bifurcation analyses and sensitivity analyses have been carried out to show the effects of interesting parameters including instant killing rateq and pesticide spraying time point θ on the successful pest control and outbreak.


## 3. Existence and Stability of Boundary Equilibria

It is easy to see that model (4) has a zero equilibrium (0,0) and its stability can be determined by the following Jacobian matrix:(8)J0,0=1-qθf0,0+1-θf0,0g0,001-q1-g0,0δ.It follows from (7) that the stability of (0,0) can be determined by the eigenvalues of following Jacobian matrix:(9)J0,0=1-qθer+1-θer00δ;that is, (0,0) is stable when the eigenvalues of the above matrix are less than one in magnitude. According to 0≤δ<1 we can see that the stability of (0,0) is determined by the value of (1-q)er. This indicates that if (1-q)er<1, then (0,0) is stable. Moreover, we note that the stability of (0,0) only depends on the intrinsic growth of the host population and the instant killing rate q, which means that both the host and parasitoid populations could die out if the pesticide is effective enough no matter when we spray the pesticides.Although the timing of pesticide applications (i.e.,θ) does not affect the stability of (0,0), from Figure 1 we note that it can significantly influence the behavior of solutions of model (4). For example, letting the parameter θ vary and fixing for all others as those in Figure 1 at which (1-q)er<1 holds true and further letting all solutions start from the same initial value, we can see that the host population tends to zero more quickly for small values of θ. This confirms that we should apply the pesticides at the ending of parasitism season or population generation, which could be best for the host control in practice.Figure 1
Illustration of the effects of pesticide application timing on the extinction of the host population. The parameter values are as follows:T=10,  Th=0.5,  K=50,  δ=0.3,  q=0.94,  γ=1,  α=0.007, and r=2.8.Furthermore, model (4) may have a boundary equilibrium (x∗,0), that is, parasitoid-free equilibrium, if x∗ satisfies the following equation: (10)1-qθfx∗+1-θf1-qx∗=1;that is, we have (11)1-qθexp⁡r1-x∗K+1-θexp⁡r1-1-qx∗K=1.Note that it is not easy to find the sufficient conditions of the existence ofx∗ for the above equation, so we first consider the following two special cases: (a) θ=0 and (b) θ=1. The analytical formula for x∗ can be easily obtained as follows:(i)
x1∗=K/rr-ln1/1-q for θ=1.(ii)
x0∗=K/r(1-q)r-ln1/1-q for θ=0.Note that, for the existence ofx∗ in the above two special cases, we need (12)r-ln⁡11-q>0;that is, x∗ is positive provided (1-q)er>1, which means that x∗ exists if (0,0) is unstable in those special cases. Obviously, we have x0∗>x1∗.It follows from(13)1-qθexp⁡r1-x∗K+1-θexp⁡r1-1-qx∗K=1that (14)er1-qθexp⁡-rx∗K+1-θexp⁡-r1-qx∗K=1,from which we can see that the necessary condition for existence of xθ∗ is er(1-q)>1.For stability of the boundary equilibrium(x∗,0), the two eigenvalues of Jacobian matrix are as follows: (15)λ1=1-qθexp⁡r1-x∗K+1-θexp⁡r1-1-qx∗K-1-qx∗Kθrexp⁡r1-x∗K+1-θr1-qexp⁡r1-1-qx∗K=1-1-qx∗Kθrexp⁡r1-x∗K+1-θr1-qexp⁡r1-1-qx∗K,λ2=1-qaTx∗1+aThθ+1-θ1-qx∗+δ.Obviously, we have λ1<1 and λ2>0.What we are interested in is how the parameters includingθ and q affect the values of x∗ and its stability. To address this, we fix the parameter values as those in Figure 2 and plot x∗ against the parameter θ for given q (denoted by xθ∗). It follows from Figure 2 that the values of x∗ are significantly changed as θ varies for large q (as shown in Figure 2 for q=0.9), and xθ∗ is a monotonically decreasing function of θ. All those indicate that carefully selecting the pesticides and designing the application time are quite important in pest control.Figure 2
The effects of instant killing rateq and pesticide application time θ on the values of xθ∗. The parameter values are as follows: T=10,  Th=0.5,  K=50,  δ=0.3,  q=0.1,  γ=1,  α=0.007, and  r=2.8.Therefore, we considerx∗ as the function of both parameters θ and q (denoted by xq,θ∗) and fix all other parameters as those in Figure 3(a). Meanwhile, the two eigenvalues λ1 and λ2 are also plotted in Figure 3(b). It follows from Figure 3 that the values of xq,θ∗ and their stabilities can be strongly influenced by the instant killing rate q and pesticide application time θ. Note that beyond the stable regions of xq,θ∗ the dynamical behavior of model (4) could be very complex. Therefore, we will address the complexity of dynamics by employing the bifurcation analysis in the following.Figure 3
The effects of instant killing rateq and pesticide application time θ on the values of xq,θ∗ and its stability. The parameter values are as follows: T=10,  Th=0.5,  K=50,  δ=0.3,  q=0.1,  γ=1,  α=0.007, and  r=2.8.
(a)
(b)


## 4. Bifurcation Analysis and Complex Dynamics

It follows from [19] that the rich dynamic behavior for model (4) can have multiple attractors coexistence and complex bifurcations. However, what we are interested in is how the key parameters of model (4) affect the successful pest control and biological implications. In particular, we will focus on the effects of the intrinsic growth rate, instant killing rate, timing of pesticide applications, and the initial densities of host and parasitoid populations on successful biological control in this section by carrying bifurcation analyses.Firstly, we show the bifurcation diagrams with respect to the intrinsic growth rate for different instant killing rate in Figure4. Figure 4(a) is similar to the bifurcation diagram shown in [19], from which we can see that model (4) has quite complex dynamics including chaotic windows, periodic windows, and quasi-periodic windows as parameter r varies. Moreover, the different attractors can coexist for certain range of parameter values of r. When we slightly increase parameter q from 0 to 0.05, we find that the bifurcation diagram is significantly changed. In particular, the multiple attractors of the host and parasitoid populations can coexist with a wide range of parameter r, as shown in Figure 4(b).Figure 4
Bifurcation diagrams with respect to the intrinsic growth rater with different q values. The parameter values are as follows: T=10,  Th=0.5,  K=50,  δ=0.3,  γ=1,  α=0.01,  θ=0,  q=0 in (a), and q=0.05 in (b).
(a)
(b)To show the effects of changing the parameter values on different attractors in more detail, we plot different attractors in Figure5, from which we can see that the instant killing rate can significantly influence the final state of host-parasitoid populations, as shown in Figures 5(a) and 5(b). Similarly, the intrinsic growth rate can also affect the final state of host-parasitoid populations, as shown in Figures 5(c) and 5(d).Figure 5
Attractors of model (4) for different parameter values r and q. The parameter values are as follows: T=10,  Th=0.5,  K=50,  δ=0.3,  γ=1,  α=0.01,  θ=0, and  q=0 in (a-b), q=0.05 in (c-d), r=2.3 in (a, c), and r=2.84 in (b, d).
(a)
q
=
0
(b)
q
=
0
(c)
q
=
0.05
(d)
q
=
0.05Thus, it follows from Figure4 that different host-parasitoid initial densities or host-parasitoid ratios can significantly affect the dynamical behavior of systems (4). For example, the solutions of model (4) starting from (5,5) and (2,1) will tend to two different attractors for fixed parameter values, as shown in Figure 6. Similar results can be obtained for the fixed parameter values in Figure 7. Those further clarify that the initial values play a key role in determining the final states of host-parasitoid populations.Figure 6
Two attractors’ coexistence. The parameter values are as follows:T=10,  Th=0.5,  K=50,  δ=0.3,  γ=1,  α=0.01,  θ=0,  q=0, and r=2.99. The initial values are (5,5) for (a-b) and (1,2) for (c-d).
(a)
(b)
(c)
(d)Figure 7
Two attractors’ coexistence. The parameter values are as follows:T=10,  Th=0.5,  K=50,  δ=0.3,  γ=1,  α=0.01,  θ=0,  q=0.05, and r=3.2. The initial values are (5,5) for (a-b) and (2,1) for (c-d).
(a)
(b)
(c)
(d)In order to address in more detail how the initial values affect the final states of host and parasitoid populations, we provide the basins of attraction of two attractors shown in Figures6 and 7, respectively, for a quite wide range of initial values of host and parasitoid population. The basins of attraction shown in Figures 8 and 9 indicate that the patterns of basins are quite complex and the final states of both host and parasitoid population can be influenced by slightly changing the initial values, which could result in complexity during the pest control due to random perturbation.Figure 8
The basins of attraction of model (4) for two attractors shown in Figure 6. The parameter values are as follows: T=10,  Th=0.5,  K=50,  δ=0.3,  γ=1,  α=0.01,  θ=0,  q=0, and r=2.99.Figure 9
The basins of attraction of model (4) for two attractors shown in Figure 7. The parameter values are as follows: T=10,  Th=0.5,  K=50,  δ=0.3,  γ=1,  α=0.01,  θ=0,  q=0.05, and r=3.2.The bifurcation diagrams with respect to survival rateδ given in Figure 10 show that the survival rate of the parasitoid population can greatly affect the dynamic behavior of model (4). Moreover, the periodicity could be strengthened as parameter δ increases and complexity can be reduced, which means that the larger the survival rate, the stronger the stability of model (4).Figure 10
Bifurcation diagram with respect to parameterθ. The parameter values are as follows: T=10,  Th=0.5,  K=50,  γ=1,  α=0.015,  θ=1,  q=0.2, and r=2.6.
(a)
(b)Furthermore, the bifurcation diagram with respect to the searching rateα shows that model (4) exists as a boundary equilibrium for small α (i.e., parasitoid-free equilibrium), and then model (4) has an interior equilibrium as α increases and exceeds some critical values (as shown in Figure 11). In this case, both host and parasitoid populations can coexist and stabilize at the interior equilibrium. The quasi-periodic solution appears once we further increase the parameter α, which could exist for a wide region of parameters.Figure 11
Bifurcation diagram with respect to parameterα. The parameter values are as follows: T=10,  Th=0.5,  K=50,  γ=1,  δ=0.5,  θ=1,  q=0.2, and r=2.8.
(a)
(b)


## 5. Sensitivity Analysis and Key Parameters

In this section, what we would like to address is how to determine the key parameters which can significantly affect the dynamics and consequently affect the pest control. Therefore, we employ the uncertainty and sensitivity analyses methods in this section. We first consider the simple case before we do this.The results shown in Figure12 based on the baseline parameter values in Figure 11 indicate that slightly changing the parameter values α and θ can strongly influence the stability of both host and parasitoid populations. Those can be confirmed by comparing Figures 12(a) and 12(b) with Figures 12(c) and 12(d), from which we can see that the slight increase of the value of α can result in strong oscillation of solutions for model (4).Figure 12
Sensitivity analyses of both parameter values ofα and θ. The parameter values are as follows: T=10,  Th=0.5,  K=50,  γ=1,  δ=0.5,  θ=1,  q=0.2, and r=2.8.
(a)
α
=
0.002
(b)
(c)
α
=
0.009
(d)To address those in more detail, in the following the sensitivity analysis was performed by evaluating the PRCCs [27, 28] for various input parameters against the solutions of both the host and parasitoid populations where we have calculated the values at 1000th generations, and then the most significant parameters (such as killing efficiency rate q and pesticide application time θ) were determined. A type of stratified Monte Carlo sampling method, named as LHS method, first proposed by Marino et al. [28] has been employed here, which has been applied to deterministic mathematical models [27]. PRCC can measure the influence of uncertainty in estimating the values of the input parameter on the imprecision in predicting the value of the output variable [27, 28]. Therefore, we performed uncertainty and sensitivity analyses for all parameters in model (4) using LHS with 1000 samples. A norm distribution function was used and tested for significant PRCCs for all parameters with mean values (i.e., the baseline values) as shown in Figures 13 and 14.Figure 13
Sensitivity analyses based on the PRCC values for the host population at 1000th generation. The baseline parameter values are as follows:T=10,  Th=0.5,  K=50,  α=0.005,  γ=1,  δ=0.5,  θ=1,  q=0.2, and r=2.8.Figure 14
Sensitivity analyses based on the PRCC values for the parasitoid population at 1000th generation. The baseline parameter values are as follows:T=10,  Th=0.5,  K=50,  α=0.005,  γ=1,  δ=0.5,  θ=1,  q=0.2, and r=2.8.PRCC results and PRCC scatter plots of all the parameters involving in model (4) have been provided in Figures 13 and 14 for both the host and parasitoid populations, respectively. In both figures the first row shows the PRCC results with sample size 1000 and all parameters were varied simultaneously. The second and third rows provide the PRCC scatter plots with sampling size 1000 and significant p<0.01. We considered absolute values of PRCC greater than 0.4 as indicating an important correlation between input parameters and output variables, values between 0.2 and 0.4 as moderate correlations, and values between 0 and 0.2 as not significantly different from zero. The positive sign of their PRCCs indicates that if the parameters are increased, the value of the solution at that generation increases (and vice versa). The negative sign suggests that if increased, the value of the solution at that generation decreases (and vice versa). All those show that the signs of PRCC values related to all parameters can help us to design the pest control strategies based on the sensitivity of each parameter with respect to the solution of model (4).It follows from Figures13 and 14 that the parameters K, r, and q are responsible for increasing the values of the host solution at 1000th generation. The PRCC value for killing rate q is a paradox result for pest control [23], because the larger q value is, the larger host population size we have. The main reason is that we apply the chemical control strategy at the end of each season or generation; that is, θ=1 or close to one. The parameters T,α, and δ are responsible for decreasing the values of the host solution at 1000th generation, so increasing all those parameters are beneficial for pest control. The most significant control parameters for the host population are T, K, α, and δ (see Figure 13).The sign of PRCC value of the parameterq related to the parasitoid population is negative and is important correlation to the size of the parasitoid population, which shows that increasing the killing rate for the host population may result in reducing the number of parasitoid populations, and in turn the host could outbreak again. The effects of all the other parameters on the parasitoid population have been provided in Figure 14.


## 6. Biological Conclusions

In this work we have proposed a discrete host-parasitoid model with pulse chemical control which can be implemented within each pest generation or parasitism season. This allows us to address the effects of timings of pesticide application and instant killing rate on the successful pest control and on the complex dynamics of proposed model.Therefore, based on the Holling II host-parasitoid model with a fraction of survival rate of parasitoid from one generation to next, we focus on the effects of timings of pesticide application and killing rate on the host or pest control and on the complex dynamics. To address those, we first discuss the existence and stability of both the host and parasitoid populations extinction equilibrium and parasitoid-free equilibrium, and then the threshold conditions which guarantee the local stabilities have been provided. The threshold condition for the stability of(0,0) equilibrium indicates that it only depends on the intrinsic growth of the host population and the instant killing rate q. This shows that both the host and parasitoid populations could be extinct simultaneously if the pesticide is effective enough no matter when we spray the pesticides. But the timing of pesticide applications can significantly influence the asymptotic behavior of solutions of model (4), as shown in Figure 1.The interesting results shown in Figures2 and 3 indicate that the different timings of pesticide applications and instant killing rate could result in significant different size of the host population and consequently influences the pest control. All those results confirm the importance of timings of pesticide application when the density dependent is involved into the model. Further, the complex dynamics including multiple attractors coexistence, quasi-periodic windows, chaotic behavior, and initial sensitivity have been studied by using numerical bifurcation analyses. Finally, the uncertainty and sensitivity of all the parameters on the solutions of both the host and parasitoid populations have been investigated, which can help us to determine the key parameters in designing the successful pest control strategy.This work focused entirely on the host-parasitoid model with chemical control, which means that we only focus on the single control tactic and study its effects on the pest control. Therefore, how can we formulate the host-parasitoid model with IPM strategies? And in particular how can we involve the biological control measure into model (4) with constant releasing rate? All those questions will be considered and studied in the near future, and our findings will be reported elsewhere.



---

*Source: 3989625-2016-03-17.xml*