- Research
- Open Access
- Published:
Dynamical modeling of the control of brown planthoppers by Beauveria bassiana and Cyrtorhinus lividipennis
Advances in Difference Equations volume 2019, Article number: 232 (2019)
Abstract
In this work, the control of brown planthopper, a major pest of rice, when a biological control agent (Cyrtorhinus lividipennis) and a pathogen (Beauveria bassiana) are utilized is investigated mathematically. An impulsive mathematical model accounted for the population densities of susceptible brown planthoppers (brown planthoppers that are susceptible to Beauveria bassiana), infected brown planthoppers (brown planthoppers that are infected by Beauveria bassiana), and Cyrtorhinus lividipennis (a natural enemy of brown planthopper) is developed. We analyze the model in terms of its stability and permanence so that we obtain the conditions that differentiate dynamic behaviors exhibited by the model. To illustrate our theoretical results, computer simulations are also presented.
Introduction
Rice is recognized as an important food crop for the population of the world [1]. To meet with the increased demand of rice due to the increased world’s population, the production of rice should be increased. Rice production has an important effect on the development of several countries’ economics. Any crises that reduce the production of this commodity can adversely affect these countries [2]. Approximately 25% of rice crop losses occur due to its pests [3]. One of the major pests of rice is brown planthopper (BPH). The infestation of BPH in a rice field can cause the damage known as hopperburn [4]. In Thailand, during the dry season of the year 2010, the outbreak of BPH caused the loss of approximately $52 million as reported in [5].
To avoid the side effects of insecticide such as biodiversity reduction and the decrease in population of natural enemies of BPH, biological control and pathogen are considered as the alternative methods for controlling the outbreak of BPH in the paddy field [6, 7].
There are many natural enemies that have been used as the biological control agent for controlling the outbreak of brown planthoppers such as Cyrtorhinus lividipennis and Lycosa pseudoannulata [8,9,10,11]. In this study, we focus on Cyrtorhinus lividipennis. The predatory activity of Cyrtorhinus lividipennis against BPH has been investigated widely and the results indicated that the Cyrtorhinus lividipennis’s preying on BPH’s eggs was an important cause of the decrease in BPH population [8, 12]. When the outbreak of BPH is severe and the use of Cyrtorhinus lividipennis alone might not be effective because the reproduction of Cyrtorhinus lividipennis is not rapid enough to control the outbreak, additional use of the pathogen Beauveria bassiana might be the appropriate way for controlling the outbreak of BPH.
Hence, we then investigate the effects of impulsive applications/releases of pathogen Beauveria bassiana and Cyrtorhinus lividipennis on the population dynamics of brown planthoppers in this paper. In the next section, we state a mathematical model that will be used to investigate the dynamic behaviors of the model when the pathogen Beauveria bassiana and Cyrtorhinus lividipennis are utilized.
Model development
Let \(x(t)\), \(y(t)\), and \(z(t)\) denote the population densities of susceptible brown planthoppers (brown planthoppers that are susceptible to Beauveria bassiana but not yet infected) at time t, infected brown planthoppers (brown planthoppers that have been infected by Beauveria bassiana) at time t, and Cyrtorhinus lividipennis (a natural enemy of brown planthopper) at time t, respectively. The following impulsive system is proposed to investigate the population dynamics of brown planthoppers when Cyrtorhinus lividipennis and the pathogen Beauveria bassiana are utilized:
with
where \(0 \leq \alpha \leq 1\). All parametric values of (1a)–(1f) are assumed to be positive. Once the rice plants in the field are not infested by BPH yet, the susceptible brown planthoppers have many rice plants to feed on. When brown planthoppers spread in the field as hopperburn, the number of rice plants available for the susceptible brown planthoppers to feed on is then decreased, and hence the logistic growth is assumed for the reproduction of susceptible BPH. The infection rate of susceptible BPH by Beauveria bassiana is assumed to vary directly to the encounters of the susceptible BPH and the infected BPH and hence the term \({a_{2}}xy\) is utilized. Since the consumption capability of a Cyrtorhinus lividipennis is limited, we then assume that the predation rate of susceptible BPH by Cyrtorhinus lividipennis is represented by the term \(\frac{{a_{3}}xz}{1+{k_{2}}x}\).
The impulsive period of application/release of the pathogen Beauveria bassiana and Cyrtorhinus lividipennis in the field is denoted by T, \(m\in {{Z}_{+}}\), \({{Z}_{+}}=\{1,2,3,\ldots\}\). \(\Delta x(t)=x({{t}^{+}})-x(t)\), \(\Delta y(t)=y({{t}^{+}})-y(t)\), and \(\Delta z(t)=z({{t}^{+}})-z(t)\). Once the pathogen Beauveria bassiana is applied, some of the population of the susceptible brown planthoppers will be infected. δ represents the fraction of the susceptible brown planthoppers that are infected by the pathogen Beauveria bassiana and becomes the infected brown planthoppers at time \(t=mT\), \(0<\delta <1\). However, there is no effect of the pathogen Beauveria bassiana on the population of Cyrtorhinus lividipennis. In addition, Cyrtorhinus lividipennis preys only on susceptible brown planthoppers, not on those infected by the pathogen Beauveria bassiana. γ represents the increase in the population of Cyrtorhinus lividipennis due to the release of Cyrtorhinus lividipennis at time \(t=mT\).
Model analysis
In what follows, we let \({{R}_{+}}\) denote \([0,\infty )\). \({R}_{+} ^{3}\) denotes \(\{Q\in {{R}^{3}}:Q=(x,y,z),\ x,y,z\in {R}_{+}\}\) and F denotes \(({{F}_{1}},{{F}_{2}}, {{F}_{3}})\). We also assume that the solution of (1a)–(1f) is piecewise continuous. To prove the main results, we first provide some lemmas which need the following definitions [13].
Definition 1
Let \(W:{{R}_{+}}\times R_{+}^{3}\to {{R}_{+}}\) satisfy the locally Lipschitzian condition in Q provided that W is continuous for \((t,Q)\in (mT,(m+1)T]\times R_{+}^{3}\) and \(\lim_{(t,S)\to (mT^{+},Q)} W(t,S)=W(mT^{+},Q)\) exists for each \(Q\in R_{+}^{3}\), \(m\in {{Z}_{+}}\). The upper right derivative of \(W(t,Q)\) with respect to the impulsive equations (1a)–(1f) is then defined as
for \((t,Q)\in (mT,(m+1)T]\times R_{+}^{3}\).
Note that \(\frac{dx}{dt}= 0\), \(\frac{dy}{dt}=0\), and \(\frac{dz}{dt}= 0\) whenever \(x(t)=0\), \(y(t)=0\), and \(z(t)=0\), respectively, for \(t\ne mT\), \(x(m{{T}^{+}})=(1-\delta )x(mT)\) and \(y(m{{T}^{+}})=y(mT)+\delta x(mT)\), \(0< \delta <1\), \(z(m{{T}^{+}})=z(mT)+ \gamma \). In addition, the smoothness properties of F also imply that the solution of the impulsive differential equations (1a)–(1f) exists and is unique [14].
Lemma 1
The solution of the impulsive differential equations (1a)–(1f), \(Q(t)=(x(t),y(t), z(t))\), is nonnegative for all \(t\ge 0\) if \(Q({{0}^{+}})\ge 0\). In addition, \(Q(t)\) is positive for all \(t\ge 0\) if \(Q({{0}^{+}})>0\).
Proof
We will prove by contradiction. Suppose that there exists \(\hat{t} \in (0,T]\) such that \(x(t)\geq 0\), \(y(t)\geq 0\), \(z(t)\geq 0\), \(x(\hat{t})=0\), \(\dot{x}(\hat{t})<0\), \(y(\hat{t})\geq 0\), and \(z(\hat{t})\geq 0\) for all \(t\in (0,\hat{t})\). From (1a), \(\dot{x}(\hat{t})=0\) which is a contradiction. Moreover, we also obtain
Therefore, \(x(t)\geq 0\) if \(x(0^{+})\geq 0\) and \(x(t)>0\) if \(x(0^{+})>0\) for \(t\in (0,T]\). The proof for \(y(t)\) and \(z(t)\) can be conducted in the same manner. □
Lemma 2
For sufficiently large t, the solution of the impulsive differential equations (1a)–(1f), \(Q(t)=(x(t),y(t),z(t))\), is bounded above.
Proof
We let \({M_{1}}=\frac{{a_{1}}{k_{1}}}{4}\), \(W(t,Q)=x+y+z\), and \(c = \min \{b_{1}, b_{2}, b_{3}\}\). It is obvious that
where \(Q_{1}=(x_{1},y_{1},z_{1})\), \(Q_{2}=(x_{2},y_{2},z_{2})\), and \(L=1\). Hence, \(W(t,Q)\) is locally Lipschitz in Q.
For \(t\ne mT\),
Hence \(D^{+}W\le -cW+{M_{1}}\).
For \(t=mT\),
Hence, for \(t\in (mT,(m+1)T]\), Lemma 2.2 of [13] implies that
That is, \(W(t)\) is uniformly ultimately bounded. Hence, when t is large enough, \(x(t)\), \(y(t)\), and \(z(t)\) are bounded above, which implies that the solution \(Q(t)=(x(t),y(t),z(t))\) of the impulsive differential equations (1a)–(1f) is bounded above as well. □
Next, let us consider the system of (1a)–(1f) when there is no brown planthopper (\(x=0\) and \(y=0\)):
We can see that the function
is a positive solution of equations (2)–(4) for \(t\in (mT,(m+1)T)\), \(m\in Z_{+}\), such that
Therefore,
is the solution of equations (2)–(4) for \(t\in (mT,(m+1)T)\), \(m\in Z_{+}\).
Lemma 3
Equations (2)–(4) have \(\tilde{z}(t)\) as a positive periodic solution. In addition, as \(t\to \infty \), \(z(t)\to \tilde{z}(t)\) for every solution \(z(t)\) of equations (2)–(4).
Therefore, at the vanishing of brown planthoppers, system (1a)–(1f) has a periodic solution
for \(t\in (mT,(m+1)T]\) and \(\tilde{z}(mT^{+})=\tilde{z}(0^{+})=\frac{ \gamma }{1-e^{-b_{3}T}}\), \(m\in {Z_{+}}\).
Theorem 1
Suppose that
and
then \((0,0,\tilde{z}(t))\), the solution of the impulsive differential equations (1a)–(1f), is locally asymptotically stable.
Proof
Consider a small perturbation of \((0,0,\tilde{z}(t))\):
Therefore, we have
where \({\varPhi }(t)\), which is the fundamental solution matrix, satisfies
with \({\varPhi }(0)=I\), the identity matrix. Hence,
We can see that the terms (*) and (**) will not be used in further calculation, then there is no need to obtain the exact expression for (*) and (**).
Linearization of (1d)–(1f) yields
The solution \((0,0,\tilde{z}(t))\) of the impulsive differential equations (1a)–(1f) is locally asymptotically stable, according to Floquet theory, if \(|\lambda _{1}|,|\lambda _{2}|,|\lambda _{3}|<1\) where \(\lambda _{1}\), \(\lambda _{2}\), \(\lambda _{3}\) are eigenvalues of
Here, the eigenvalues of P are
Since \(0<\delta <1\) and (5)–(6) hold, then all eigenvalues are positive and less than 1. Hence, the solution \((0,0,\tilde{z}(t))\) of the impulsive differential equations (1a)–(1f) is locally asymptotically stable. □
System permanence
Definition 2
If there are constants \(\bar{m}, \bar{M}>0\) and \(t_{0}>0\) for which every solution with positive initial condition \(x(0^{+})\), \(y(0^{+})\), and \(z(0^{+})\),
for all \(t>{t_{0}}\), system (1a)–(1f) is said to be permanent.
Theorem 2
Suppose that
and
where \(A\equiv a_{1}-b_{1}\) and \(B\equiv (\frac{a_{1}}{k_{1}}+a _{2}+a_{3} ) (\frac{a_{1}k_{1}}{4c}+\frac{\gamma e^{cT}}{e ^{cT}-1} )\), \(c=\min \{b_{1},b_{2},b_{3}\} \). The impulsive differential equations (1a)–(1f) are permanent if (7) and (8) hold.
Proof
Let \(Q(t)\) be a solution of the impulsive differential equations (1a)–(1f), \(Q(t)=(x(t),y(t),z(t))\), with \(x(0^{+}),y(0^{+}),z(0^{+})>0\). Then Lemma 2 guarantees that \(Q(t)\) is bounded above when t is large enough. Hence, a constant \(\bar{M}>0\) exists for which, when t is sufficiently large, we have \(x(t), y(t),z(t)\le \bar{M}\).
Consider (1c) and (1f), we know that
Consider the comparison system
We can see that, for \(t\in (mT,(m+1)T)\), \(m\in Z_{+}\),
is a positive solution of the comparison system (12)–(14) such that
Therefore, the solution of this comparison system (12)–(14) is
and \(r(t)\rightarrow \tilde{r}(t)\) as \(t\rightarrow \infty \).
According to [14], we obtain \(z(t)\geq r(t)\) by the comparison theorem and hence, when t is large enough,
Next, we will show that there exists a constant \({m_{2}}>0\) for which \(x(t)\geq {m_{2}}\) when t is large enough.
Step I. By contradiction, we will show for \(t\geq t_{1}\) that there is \(t_{1}>0\) for which \(x(t)\ge {m_{3}}\). For all \(t\geq 0\), we suppose that \(x(t)<{m_{3}}\).
Consider (1a) and (1d). Since Lemma 2 implies that there is \(t_{1}>0\) in which we can select \(\bar{M}= \frac{a_{1}k_{1}}{4c}+\frac{\gamma e^{cT}}{e^{cT}-1}\) where \(c=\min (b_{1},b_{2},b_{3})\) for which, for \(t\geq t_{1}\), we have \(x(t), y(t),z(t)\le \bar{M}\), then
for \(t\geq t_{1}\), with \(x(mT^{+})=(1-\delta )x(mT)\), \(t = mT\).
Let \(N\in {Z_{+}}\) and \(NT\ge {T_{1}}\). We integrate (18) over \((mT,(m+1)T], m\ge N\) and obtain
where \(\eta \equiv (1-\delta )\exp [a_{1}-b_{1}- (\frac{a _{1}}{k_{1}}+a_{2}+a_{3} )\bar{M} ]T\).
Consider
Since (7) and (8) hold, then \(\ln \eta > 0\), which implies that \(\eta > 1\). Therefore, \(x((m+k)T)\ge x(mT)\eta ^{k}\to \infty \) as \(k\to \infty \), and hence \(x(t)\) is not bounded above when t is large enough, and this is a contradiction. This means that there is \({t_{1}}>0\) for which \(x(t_{1})\ge {m_{3}}\).
Step II. For all \(t>{t_{1}}\), if \(x(t)\ge {m_{3}}\), the proof is then complete, or else there must be \(t' > t_{1}\) for which \(x(t')<{m_{3}}\). Next, we consider the following possible cases when we let \(t^{*}=\inf_{t>t_{1}}\{t:x(t)< m_{3}\}\).
Case 1. \(t^{*}=k_{1}T\) for some \(k_{1} \in Z_{+}\). This means that for \(t \in [t_{1},t^{*}]\), \(x(t)\geq m_{3}\) and
Next, \(k_{2}, k_{3} \in Z_{+}\) are chosen so that
where \(\eta _{1}=a_{1}-b_{1}- (\frac{a_{1}}{k_{1}}+a_{2}+a_{3} ) \bar{M} < 0\).
Now, we let \(T'=k_{2}T+k_{3}T\) and claim that \(t_{2} \in (t ^{\ast },t^{\ast }+T']\) must exist for which
or else, for \(t^{\ast }+k_{2}T \leq t \leq t^{\ast }+T'\), (18) holds and
For \(t \in [t^{\ast },t^{\ast }+k_{2}T]\),
By integrating (22) over \([t^{\ast },t^{\ast }+k_{2}T]\), we obtain
Therefore,
and a contradiction occurs. Thus, \(t_{2}\in (t^{\ast },t^{\ast }+T']\) must exist for which \(x(t_{2})>m_{3}\).
Letting \(\tilde{t}= \inf_{t>t^{\ast }}\{t:x(t)>m_{3}\}\). Since \(x(t)\) is left continuous and \(x(t^{+})=(1-\delta )x(t) \leq x(t)\) when \(t=mT\), we can conclude that \(x(t) \leq m_{3}\) whenever \(t \in (t^{ \ast },\tilde{t})\) and \(x(\tilde{t})=m_{3}\).
Next, we assume for \(t \in (t^{*},\tilde{t})\) that \(t\in (t ^{*},(l-1)T,t^{*}+lT]\) where \(l\in Z_{+}\) and \(l \leq k_{2}+k_{3}\).
Consider (21), we get
Therefore, we can conclude that \(x(t) \geq m'_{1}\) whenever \(t\in (t^{\ast },\tilde{t})\) and \(x(\tilde{t})\geq m_{3}\). With similar argument for \(t>\tilde{t}\), when t is sufficiently large, we obtain \(x(t)\geq m_{2}>0\).
Case 2. \(t^{*} \neq mT\) for all \(m \in Z_{+}\). For \(t \in (t _{1},t^{*})\), we obtain \(x(t) \geq m_{3}\) and \(x(t^{\ast })=m_{3}\). Suppose \(t^{*} \in (k'_{1}T,(k'_{1}+1)T)\) for some \(k'_{1}\in Z_{+}\).
Case 2.1. For all \(t \in (t^{*},(k'_{1}+1)T)\), \(x(t) \leq m _{3}\). We claim that \(t'_{2} \in [(n'_{1}+1)T,(n'_{1}+1)T+T']\) exists for which \(x(t'_{2})>m_{3}\). Otherwise, similar to Case 1, we get
For \(t \in (t^{*},(k'_{1}+1)T)\), (22) holds on \([t^{*},(k'_{1}+1+k _{2}+k_{3})T]\) and \(x(t) \leq m_{3}\), so that we obtain
Therefore,
and a contradiction occurs.
Next, letting \(\bar{t}= \inf_{t>t^{\ast }}\{t:x(t)>m_{3}\}\). Therefore, for \(t\in (t^{\ast },\bar{t})\),
and
Now we assume for \(t\in (t^{*},\bar{t})\) that \(t\in (k'_{1}T+(l'-1)T,k'_{1}T+l'T]\) where \(l' \in Z_{+}\) and \(l' \leq 1+k_{2}+k_{3}\). Thus, we obtain
So, for \(t \in (t^{\ast },\bar{t})\), \(x(t) \geq m_{1}\). For \(t>\bar{t}\), since \(x(\bar{t})\geq m_{3}\), the similar arguments can be applied. We thus get \(x(t) \geq m_{2} >0\) when t is sufficiently large.
Case 2.2. There is \(t \in (t^{*},(k_{1}+1)T)\) for which \(x(t)>m_{3}\). Therefore, we let \(\hat{t}= \inf_{t>t^{*}}\{t:x(t)>m _{3}\}\). Hence, for \(t\in (t^{\ast },\hat{t})\), we get \(x(t)\leq m _{3}\) and \(x(\hat{t})=m_{3}\).
Next, we integrate (22) on \((t^{*},\hat{t})\) where we note that (22) satisfies for \(t\in (t^{*},\hat{t})\) and obtain
With the fact that \(x(\hat{t})\geq m_{3}\), the above argument can be applied again for \(t>\hat{t}\). Thus, we obtain \(x(t) \geq m_{2} > 0\) for all \(t \geq t_{1}\).
Finally, we can also show that there exists a constant \({m_{4}}>0\) for which \(y(t)\geq {m_{4}}\) for sufficiently large t in the same manner as we showed above that \(x(t)\geq m_{2}\) for sufficiently large t by using the fact that \(x(t)\geq m_{2}\) for sufficiently large t, and hence the proof is omitted here.
Therefore, we can choose \(\bar{m}=\min \{m_{1},m_{2},m_{4}\}\) so that \(x(t),y(t),z(t)\geq \bar{m}>0\), and hence the system is permanent provided that (7) and (8) hold and the proof is complete. □
Numerical simulations
Computer simulations are presented in this section to illustrate our theoretical results in Sect. 3.
Here, two sets of parametric values of the model are chosen as examples to illustrate the theoretical results in which the parameters \(a_{3}\) and \(k_{2}\) are obtained from the literature [15].
A simulation result of the impulsive system of equations (1a)–(1f) with the parametric values \({a_{1}}=0.21\), \({a_{2}}=0.1\), \({a_{3}}=0.247\), \({b_{1}}=0.2\), \({b_{2}}=0.1\), \({b_{3}}=0.1\), \({k_{1}}=0.1\), \({k_{2}}=0.017\), \(\alpha =0.1\), \(\gamma =0.9\), \(\delta =0.2\), \(T=14\), \(x(0)=5\), \(y(0)=5\), and \(z(0)=5\), in which all the conditions in Theorem 1 are satisfied, are presented in Fig. 1. In Fig. 1, the time courses of the population densities of the susceptible BPH, infected BPH, and Cyrtorhinus lividipennis as well as the projection of the solution trajectory onto the \((x,z)\)-plane are presented. The solution trajectory in this case tends to the oscillating solution \((0,0,\tilde{z})\) for which both susceptible and infected BPH vanish as predicted in Theorem 1.
On the other hand, Fig. 2 shows a simulation result of the impulsive system of equations (1a)–(1f) with the parametric values \({a_{1}}=0.38\), \({a_{2}}=0.1\), \({a_{3}}=0.247\), \({b_{1}}=0.25\), \({b_{2}}=0.1\), \({b_{3}}=0.1\), \({k_{1}}=0.3\), \({k_{2}}=0.017\), \(\alpha =0.5\), \(\gamma =0.9\), \(\delta =0.2\), \(T=60\), \(x(0)=5\), \(y(0)=5\), and \(z(0)=5\) in which all the conditions in Theorem 2 are satisfied. The time courses of the population densities of the susceptible BPH, infected BPH, and Cyrtorhinus lividipennis as well as the projection of the solution trajectory onto the \((x,z)\)-plane are as shown in Fig. 2. The solution of the system in this case is permanent as predicted in Theorem 2. We can see that the population densities of both susceptible and infected BPH can be controlled within a certain level, while its natural enemy Cyrtorhinus lividipennis also survives in the paddy field.
Conclusion
We investigate the dynamic behaviors of the populations of susceptible and infected BPH when Cyrtorhinus lividipennis and the pathogen Beauveria bassiana are utilized to control the population of BPH in the paddy field mathematically. Here, the pathogen Beauveria bassiana does not have an effect on Cyrtorhinus lividipennis.
Brown plant hoppers (BPH) are rice’s insect pests. Therefore, the aim of this work is to obtain the conditions on the system parameters for which the populations of both susceptible and infected BPH tend to zero level as time passes. However, the cost for controlling BPH to zero level might be too high and the control of BPH levels \(x(t)\) and \(y(t)\) to lie within a certain range might be a better option in an economic point of view. Hence, we also provide the conditions for which the system is permanent and the levels of BPH can be controlled to lie within a certain range as well.
The examples of two sets of parameters are selected as examples to illustrate the theoretical results in which some parametric values are obtained from the literature [15]. Numerical simulations are as shown in Figs. 1 and 2. We can see that in Fig. 2, the population of BPH could be controlled below certain ranges, and we can also observe that the population of Cyrtorhinus lividipennis is approximately twice of the population of BPH. This result corresponds to what has been observed in the rice field in Thailand [16] that the population of BPH could be controlled to be lower than a certain level if the population of Cyrtorhinus lividipennis in the field is approximately twice of the population of BPH in the rice field.
In addition, the conditions in Theorem 1 and 2 depend on the duration T between the two consecutive applications of Beauveria bassiana and Cyrtorhinus lividipennis. Hence, the appropriate duration T could play an important role in controlling the population of BPH in the rice field. The current practices in Thailand that natural enemies of BPH or pathogens will be applied in the field whenever the spread of BPH in the rice field is detected would take some time before the population of BPH can be controlled. Hence, the applications of Beauveria bassiana and Cyrtorhinus lividipennis in the rice field with the appropriate duration T could lead to the more efficient control of BPH in the rice field.
References
- 1.
Kole, C.: Cereals and Millets, Genome mapping and molecular breeding in plants, vol. 1. Springer, Berlin (2006)
- 2.
Lichtfouse, E., Goyal, A.: Sustainable Agriculture Reviews, Sustainable Agriculture Reviews, vol. 16. Springer, Switzerland (2015)
- 3.
Oerke, E.C.: Crop losses to pests. J. Agric. Sci. 144, 13–43 (2006)
- 4.
Stout, M.J.: Chapter I: Host-plant resistance in pest management. In: Abrol, D.P. (ed.) Integrated Pest Management. Elsevier, Amsterdam (2014)
- 5.
Escalada, M.M., Luecha, M., Heong, K.L.: Chapter 10: Social impacts of planthopper outbreaks: case study from Thailand. In: Heong, K.L., Cheng, J.A., Escalada, M.M. (eds.) Rice Planthoppers: Ecology, Management, Socio Economics and Policy, pp. 209–226. Zhejiang University Press, Hangzhou; Springer, Dordrecht (2015). https://doi.org/10.1007/978-94-017-9535-7
- 6.
Huang, J., Qiao, F., Zhang, L., Rozelle, S.: Farm Pesticide, Rice Production and Human Health, Centre for Chinese Agricultural Policy. Chinese Academy of Sciences, Beijing (2000)
- 7.
Rola, A.C., Pingali, L.P.: Pesticides, Rice Productivity, and Farmers’ Health: An Economic Assessment. International Rice Research Institute, Baños Los, Philippines; World Resource Institute, Washington D.C. (1993)
- 8.
Katti, G., Pasalu, I.C., Padmakumari, A.P., Padmavathi, C., Jhansilakshmi, V., Krishnaiah, N.V., Bentur, J.S., Prasad, J.S., Rao, Y.K.: Biological control of insect pests of rice. Technical Bulletin, Directorate of Rice Research, Rajendranagar, Hyderabad, AP, India 22 22 (2007)
- 9.
Sigsgaard, L.: Early season natural control of the brown planthopper Nilaparvata lugens: the contribution and interaction of two spider species and a predatory bug. Bull. Entomol. Res. 97, 533–544 (2007)
- 10.
Khan, A.A., Misra, D.S.: Biology of woft spider Lycosa pseudoannulata Boesenberg and Strand (Araneae: Lycosidae). Environ. Ecol. 26, 796–799 (2008)
- 11.
Samiayyan, K.: Chapter 15: Spider-the generalist super predators in agro-ecosystem. In: Abrol, D.P. (ed.) Integrated Pest Management, 1st edn., pp. 283–310. Academic Press, Tamil nadu Agricultural University, Coimbatore (2014)
- 12.
Preetha, G., Stanley, J., Suresh, S., Samiyappan, R.: Risk assessment of insecticides used in rice on miridbug, Cyrtorhinus lividipennis Reuter, the important predator of brown planthopper, Nilaparvata lugens (Stal.). Chemosphere 80, 498–503 (2010)
- 13.
Liu, B., Zhi, Y., Chen, L.: The dynamics of a predator-prey model with Ivlev’s functional response concerning integrated pest management. Acta Math. Appl. Sin. 20(1), 133–146 (2004)
- 14.
Lakshmikantham, V., Bainov, D.D., Simeonov, P.S.: Theory of Impulsive Differential Equations. World Scientific, Singapore (1989)
- 15.
Jiang, X., Huang, Q., Ling, Y., Chen, Y., Xiao, G., Huang, S., Wu, B., Huang, F., Cai, J., Long, L.: Functional and numerical responses of Cyrtorhinus lividipennis to eggs of Nilaparvata ligens are not affected by genetically modified herbicide-tolerant rice. J. Integr. Agr. 14(10), 2019–2026 (2015)
- 16.
Rice Knowledge Bank, http://www.ricethailand.go.th/rkb/disease
Funding
We acknowledge the support by the Centre of Excellence in Mathematics, the Commission on Higher Education, Thailand.
Author information
Affiliations
Contributions
All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Competing interests
The author declares that they have no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
About this article
Cite this article
Rattanakul, C. Dynamical modeling of the control of brown planthoppers by Beauveria bassiana and Cyrtorhinus lividipennis. Adv Differ Equ 2019, 232 (2019). https://doi.org/10.1186/s13662-019-2171-8
Received:
Accepted:
Published:
Keywords
- Mathematical model
- Brown planthopper
- Cyrtorhinus lividipennis
- Beauveria bassiana