Research Article | Open Access
Jing Hui, Jian-Hua Pang, Dong-Rong Lin, "Mathematical Analysis of an Epidemic-Species Hybrid Dynamical System", Mathematical Problems in Engineering, vol. 2014, Article ID 810312, 11 pages, 2014. https://doi.org/10.1155/2014/810312
Mathematical Analysis of an Epidemic-Species Hybrid Dynamical System
We consider an epidemic-species hybrid dynamical system. The disease is spread among the prey only and the infected prey can reproduce virus. The predator only eats the infected prey. Mathematical analyses are given for the system with regard to the existence of equilibria, local stability, Hopf bifurcation, and the orbital stability of the Hopf bifurcating limit cycle. We further analyse the system under impulsive releasing of virus and predator.
Epidemic models and species models have received much attention from scientists, respectively. There are many literatures about them. Now we consider an epidemic-species hybrid dynamical system which is motivated by the integrated pest management (IPM). As we know, IPM is becoming more and more popular among farmers, researchers, and policy makers. IPM seeks to minimize reliance on pesticides by emphasizing the contribution of other control methods, including biological control, host-plant resistance breeding, and cultural techniques. Potentially, the use of viruses, fungi, and bacteria is one of the most effective biological methods for controlling pests. One reason for this optimism is that the generation time lapse of microbes is much shorter than that of an insect pest . The role of microbial pesticides in the integrated management of insect pests has been recently reviewed for agriculture [2–5], forestry [6, 7], and public health . In most cases no single microbial control agent will provide sustainable control of an insect pest or complex of pests. As components of an integrated approach, entomopathogens can provide significant and selective insect control. In the not too distant future we can envision a broader appreciation for the attributes of entomopathogens and expect to see synergistic combinations of microbial control agents with other technologies (in combination with semichemical technology, soft chemical pesticides, other natural enemies, resistant plants, chemigation, remote sensing, etc.) that will enhance the effectiveness and sustainability of integrated control strategies.
There is a vast amount of literature on the applications of microbial disease to suppress pests [9, 10]. But there are only a few papers on mathematical models of the dynamics of microbial disease in pest control [11, 12]. On the other hand, there are many mathematical models of infectious diseases of human beings [13–16]. All the earlier pest control models are normally two-dimensional consisting of susceptible and infected pests, and infectious disease models are at most three-dimensional with a susceptible, infected, and removable class; the dynamics of the diseases-causing organism is not considered in the model. However, presently the insect pathogens are used in two ways . In the first method, a small amount of pathogens is introduced in the pest population with the expectation that it will generate an epidemic and will be subsequently endemic. The success of the approach depends on the survival of the microbes which in turn depends on environmental factors, such as temperature, humidity, and crop conditions. In the second method, an insect pathogen is used like a nonresidual chemical insecticide. In this case, it is applied whenever a pest population is at an economically extreme level for pest damage and there is no expectation that the pathogen will survive for an appreciable period. In this paper, we examine the use of pathogens in a four-dimensional prey-predator model for the agricultural or forestry ecosystem with the assumption that viral disease spreads only among the insect pests. In nature, the infection by baculovirus begins when an insect eats virus particles on a plant, perhaps from a sprayed treatment. Virus infection causes the cell lysis in the host body and produces more virus particles (virus replication), until the cell and ultimately the insect die. Most baculoviruses cause the host insect to die in a way that maximizes the chance that other insects come in contact with the virus and become infected by eating the foliage that has been contaminated by virus-killed larvae. We further assume that the natural enemy (predator) in the system survives on the infected prey. This is due to the fact that the viral infection makes some behavioral changes and sublethal effects on host, which make them more vulnerable to predation by natural enemies than healthy hosts before death .
The dynamical behavior of the considered system is investigated from the point of view of stability and persistence. The model shows that infection can be sustained only above a threshold force of infection resulting from virus replication parameter . Still, on increasing the value of , the endemic equilibrium bifurcates towards a periodic solution. The stability of the limit cycles arising from Hopf bifurcation is analysed using Poore’s condition . Whereas, in consideration of the practice in pest control, the natural enemy and the pathogens are introduced discontinuously at some fixed moment, the impulsive releasing of the virus and natural enemy is added in the system. We prove that the pest extinction periodic solution exists and is globally asymptotically stable when the impulsive period is less than the critical value; otherwise, the system can be permanent. The arrangement of the paper is as follows. In Section 2 we describe the model formulation. Section 3 shows its equilibria. Section 4 discusses the local stability and Hopf bifurcation. In Section 5 we give the details about the orbital stability of the periodic orbit arising from Hopf bifurcation by Poore’s condition. In Section 6 we consider the corresponding impulsive system by releasing virus and natural enemy at some fixed moment. Lastly in Section 7 we end the paper with a concluding discussion.
2. The Model
The model contains three species, namely, pest, predator, and virus. The total pest population is divided into two subpopulation classes: the susceptible pest, denoted by , and the virus-infected pest, denoted by . Therefore, at any time , . We assume that only susceptible pest is capable of reproducing with logistic law; that is, the infected class of pest is removed by lysis before having the possibility of reproducing. However, they still contribute with to population growth towards the carrying capacity . is the intrinsic birth rate. Then the growth equation of susceptible pest is given by A susceptible pest becomes infected as the two mentioned methods ahead. The incidence is assumed to be the simple mass action incidence ; represents the effective per pest contact rate with viruses. Hence, the evolution equation of the susceptible pest is The infected pest has a latent period, the period between the instant of infection and that of lysis, during which the virus is reproduced inside the larval tissue. The lysis death rate gives a measure of such latency period , given by . The equation of infected class takes the form The lysis largely produces virus polyhedra or polyhedral inclusion bodies (PIB) on average PIB per insect; is called the virus replication parameter. The virus particles have natural mortality denoted by due to temperature changes, enzymatic attack, pH dependence, and so forth. The equation of the virus is given by We assume the natural enemy or predator consumes only the infected prey , the consuming means is the standard bilinear expression , and is the predation coefficient. We assume the natural enemies in the system remain unharmed from the virus. Then the equations of the infected pest and the natural enemy are where is the death rate of the predator, is the density restriction coefficient of predator, and is the conversion factor for the predator. Keeping these in view, the dynamics of the epidemic-species may be governed by the following autonomous system of differential equations: In the above model, all of the coefficients, , , , , , , , , , and , are positive real numbers.
If using the dimensionless time and the transformations , , , and , we have the dimensionless form of model (6). For convenience, we still write in place of ; the corresponding dimensionless form of model (6) is given as follows: where we denote , , , , and .
3. Equilibria of System (7)
Equilibria of model (7) are obtained by solving . It can be checked that model (7) has the following three boundary equilibria: the vanishing equilibrium ; the disease-free equilibrium ; and the predator-free equilibrium , where , , and ; when , the equilibrium is nonnegative. The positive equilibrium is obtained by solving the following equations: From the fourth equation of the above equations, we have ; substituting it into the first equation, then ; from the third equation of (8), ; substituting , , and into the second equation, we know the infected pest satisfies the following equation: , so . For the positivity of , we assume and . If , then the predator-free equilibrium and the positive equilibrium do not exist, which implies the virus replication parameter is too small to support the virus invasion in insect pest. When the virus replication parameter is larger than and very close to , then the predator-free equilibrium will collapse to the disease-free equilibrium . Thus it can be seen that the virus replication parameter plays an important role for the dynamics of system (7). In summary, we have the following theorem about the equilibria of system (7).
Theorem 1. System (7) always has two equilibria, namely, the vanishing equilibrium and the disease-free equilibrium ; if , there is one more predator-free equilibrium ; further, if the virus replication parameter satisfies and , then the positive equilibrium exists.
4. Local Stability and Hopf Bifurcation
The dynamical behaviour of equilibria can be studied by computing variational matrices corresponding to each equilibrium. The variational matrix of (7) is
For the vanishing equilibrium , correspondingly, variational matrix is Obviously, is unstable saddle point with three-dimensional stable manifold and one-dimensional unstable manifold.
For the disease-free equilibrium , the corresponding variational matrix is The characteristic equation is ; that is, . When , there are four real negative eigenvalues and consequently the disease-free equilibrium is locally asymptotically stable; when , corresponding to the disease-free equilibrium , there are a simple eigenvalue 0 and another three negative eigenvalues. The system (7) enters into a saddle-node bifurcation at ; when , there are three real negative eigenvalues and one positive eigenvalue so the disease-free equilibrium is unstable.
For the predator-free equilibrium , the corresponding variational matrix is The characteristic equation is , where , , and . The characteristic equation has clearly one negative real root, namely, , and another three roots are given by the equation . , , and can be looked as the function of ; since and the predator-free equilibrium exists if , . Now we apply Routh-Hurwitz criteria and the definition of a simple criterion of Hopf bifurcation  to show the local stability of and the Hopf bifurcation at . That is to say, for any , is locally asymptotically stable if and only if , , and . A Hopf bifurcation of the predator-free equilibrium occurs at if and only if , .
In fact, , , is obvious under the condition . For the function , , since and , based on the intermediate value theorem, there exists a so that . In addition, and at ; at the same time , so . Thus at which implies that the predator-free equilibrium is locally asymptotically stable in and unstable in ; further a Hopf bifurcation of the predator-free equilibrium occurs at ; for decreasing , it approaches a periodic solution, whose period is , where .
For the positive equilibrium , the corresponding variational matrix is For the local stability of positive equilibrium , we have the following theorem.
Theorem 2. The positive equilibrium is locally asymptotically stable if ; the stable manifold is .
Proof. We choose the Lyapunov function as follows: where , , is to be determined suitably. Obviously is positive definite. The derivative of along the solution of the equation , where , is as follows: The symmetric matrix corresponding to is the following:The positive equilibrium is locally asymptotically stable if is negative definite, which in turn follows if the symmetric matrix is negative definite; that is, the odd rank principal minor in order is negative and the even rank principal minor in order is positive, which in turn follows if(i);(ii);(iii);(iv)[the left expression of the inequality (iii)].We choose , , and ; then from the above expression, we know if , that is, , for the symmetric matrix , the odd rank principal minor in order is negative and the even rank principal minor in order is positive. That is to say, if , the positive equilibrium is locally asymptotically stable. Furthermore, to find out the basin of attraction, we choose the Lyapunov function as follows: Differentiating along the solution of the model and choosing , we can obtain the stable manifold . The proof is completed.
Now, we will find out the conditions for which the equilibrium enters into Hopf bifurcation. The characteristic equation for is where Since , , and can be written as , +, and , respectively, under the condition , it is not easy to know , . Let be the following continuously differentiable function: . A Hopf bifurcation of the equilibrium occurs at if and only if (i) , , where is purely imaginary at , and (ii) all other eigenvalues are of negative real parts.
The existence of can be obtained by solving . At , ; that is, , which deduces that the characteristic equation can be written as . So there exist a pair of purely imaginary , respectively; other eigenvalues , satisfy , . Denote ; then since , , and . If and are complex conjugates, then from we know ; if and are real roots, then based on the fact that ; furthermore, since , then , . In summary, for the characteristic equation, besides and , other eigenvalues are negative real parts. Now we show . In fact, let ; substituting it into the characteristic equation for the positive equilibrium and calculating the derivative on , the following expression can be given through some calculations and arrangement: where It deduces that Solving the above equations on , , we have At this point, the conditions of Hopf bifurcation for positive equilibrium are verified. In summary, we have the following theorem.
Theorem 3. There exists a ; the Hopf bifurcation for the positive equilibrium occurs at . The period of the bifurcating periodic orbits close to is , where .
In the next section we will consider the stability of bifurcating periodic orbits applying Poore’s condition .
5. Stability of the Bifurcating Periodic Solution
We apply Poore’s condition for verification of the orbital stability of the Hopf bifurcating limit cycle. For convenience, we list the preliminaries firstly. Let a real, -dimensional (), first-order system of autonomous differential equations be of the form where denotes a vector of -real parameters. There must firstly exist a combination of the parameters, say, ; a critical point , such that the variational matrix has exactly two, nonzero, purely imaginary eigenvalues, say, with ; and other eigenvalues with nonzero real parts. To vary one, some, or all of the parameters, an -dimensional vector function is introduced with the property that , and hence we confine our analysis to the following system of ODE: It follows from the definition of that and eigenvalues of are the same as those of . It is assumed that , where and is a domain in containing and . As there is no nonzero eigenvalue of the variational matrix, , and hence the implicit function theorem guarantees the existence of a critical point which is -times continuously differentiable in and satisfies for in a small neighborhood of . Using this definition of , a change of variables is introduced: This reduces the differential (25) to the following form: where and is defined as given in . Thus, the problem of periodic solutions of (25) is reduced to a perturbation problem in the small parameter of .
Now, the stability information of the bifurcating periodic orbits is contained in the following theorem of Poore, coupled with an algebraic expression, which completely reduces the determination of stability to an algebraic problem. By the assumptions of Theorem 2.3 , the differential equation in (27) is continuously differentiable in , where and , and in the function in a neighborhood of the periodic orbit. Thus, the existing periodic orbit will be asymptotically orbitally stable with asymptotic phase if of the characteristic multipliers of the variational equation have moduli less than one. The following theorem developed by Poore about the modulus of each of the characteristic multipliers reduces the condition of orbital stability of the Hopf bifurcating limit cycle to an algebraic expression as follows.
Theorem 4. Let satisfy the hypotheses in Theorem 2.3  and let and denote the left and right eigenvectors, respectively, for the eigenvalues of the variational matrix . If and are normalized in the sense that , then where , , , is the value of the two-dimensional vector , which occurs in , and denotes the complex conjugate of .
Written out in component form, the above expression reduces to where the repeated indices within each term imply a sum from 1 to and all the derivatives of are evaluated at the equilibrium . The sign of the real and imaginary parts of the right-hand side of expression in Theorem 4 is independent of the choice of in and the eigenvectors and so long as and . So positivity of the real part of the above expression in parenthesis really indicates the orbital stability of the periodic solution arising out of Hopf bifurcation.
For model (7) we considered, let the right four equations be , , , and , respectively. To verify that the conditions of Theorem 4 are satisfied, we will calculate the second- and the third-order derivatives of , , as follows: Other second-order and third-order derivatives are all zero. The inverse matrix of and the left and right eigenvectors and , respectively, in the sense that , of the variational matrix can be calculated. Based on the preliminary work, we can calculate expression (29) in parenthesis as follows: where Putting the values of , , , and and components of the matrix in terms of the parameters of the model, the positivity of the real part of expression (29) in parenthesis can be deduced. This in turn indicates the orbital stability of the limit cycle arising out of Hopf bifurcation.
6. Extinction and Permanence with Impulsive Effect
In the sections ahead, we focus our attention on model (6) or (7), in which the growth of the predator and the virus is continuous. But in fact, it is needed to throw in the predator or virus according to the pest control at some moments for better effects. So in this section, we develop model (6) into an impulsive system as follows and the presence of impulses gives the system a mixed nature, both continuous and discrete: where the predator and virus are thrown at , , is the period, and , are the amount of impulsive stocking on the predator and the virus; , are positive constants. Other parameters are the same as model (6). , , , and represent the densities of the susceptible pest, the infected pest by virus, the predator, and the virus, respectively: , , , and ; , , , and are left continuous at ; that is, , , , and . In this section, we mainly consider the global stability of pest-eradication solution and the uniform permanence of system (33); before introducing our main results, we give some preliminaries needed next.
For the system we have the following lemma.
Lemma 5. System (34) has a unique positive periodic solution with period if and only if the initial value , where , , and is globally asymptotically stable; in addition, every solution of system (34) can be written as , .
For the system we have the following lemma.
Lemma 6. System (35) has a unique positive periodic solution with period if and only if the initial value , where , , , and is globally asymptotically stable.
Lemma 7. There exists a constant such that , , , and for each solution of (33) with being large enough.
The solution of system (33) corresponding to is called pest-eradication solution whose explicit form can be obtained by Lemmas 5 and 6, where and , , with initial values , as Lemmas 5 and 6, respectively. Now we consider the stability of pest-eradication solution. The Jacobi matrix at is we can calculate the monodromy matrix , where is identity matrix and is just as follows:there is no need to calculate the exact forms of since they are not required for the following analyses. The Floquet multipliers are , , , and , respectively. Obviously, , so the stability of is decided by whether ; that is, when , the pest-eradication solution is locally stable; further we can prove that it is globally asymptotically stable.
Theorem 8. If , then the pest-eradication solution of system (33) is globally asymptotically stable.
The pest-eradication solution is globally asymptotically stable when ; that is to say, the pest population is eradicated totally. But in practice, we only need to control the pest population under the economic threshold level (ETL), instead of eradicating it totally, and hope the pest population, the predator (natural enemy), and the virus population can coexist when the pests do not bring about immense economic losses. As for this, we have the uniform permanence of system (33); see the next theorem.
Theorem 9. If , then system (33) is uniformly permanent.
The proofs of the lemmas and the theorems listed in this section are all omitted here since they are similar to Hui and Zhu .
7. Numerical Simulations
In this section, we use numerical simulations to illustrate the dynamical behaviors of system (7) and (33). We choose as bifurcation parameter; for system (7), the other parameter values are taken as , , , , , , and ; for system (33), the other parameter values are taken as , , , , , , , , , , , and .
Firstly, Figure 1 illustrates that the disease-free equilibrium is locally asymptotically stable when . It is calculated that () and with . Therefore, the disease-free equilibrium is unstable and is locally asymptotically stable when (see Figure 2). A Hopf bifurcation of the predator-free equilibrium occurs at due to and . In this case, there is a periodic solution for system (7) (see Figure 3). When , the positive equilibrium exists and it is locally asymptotically stable. Figure 4 illustrates the stability of the positive equilibrium .
Figure 5 is the bifurcation graph of system (7). The virus replication factor is chosen as the bifurcation parameter. When , there is a stable positive equilibrium ; when , the positive equilibrium is unstable and the Hopf bifurcation appears at ; the stable periodic solution appears. Figure 6 illustrates the bifurcation phenomenon of the impulsive system (33). When , there exists a stable periodic solution; when , the system (33) enters into the chaos state.
In this paper, from Sections 1 to 5 we focus our attention on an epidemic-species hybrid system without impulse influence. From our analyses, we found the virus replication factor plays a very important role for dynamics of system. The system always has two equilibria, namely, the vanishing equilibrium and the disease-free equilibrium . When the virus replication factor is larger than its natural death rate , that is, , the predator-free equilibrium appears; furthermore, the virus replication factor satisfies and ; the positive equilibrium exists. The vanishing equilibrium is a saddle point with three-dimensional stable and one-dimensional unstable manifold. Biologically, it means that in the absence of prey, that is, the insect pest, the predator and virus population becomes extinct. When , the disease-free equilibrium is locally asymptotically stable; it means that if the virus replication parameter or virulence of the applied virus on the system is less than its natural death rate , then only viruses obtained from lysis of the infected pest are not enough to maintain the viral infection in the pest population. When , the system passes through saddle-node bifurcation and for value of higher than , the system gives two branches of equilibria, namely, the disease-free equilibrium and the predator-free equilibrium , one hyperbolic equilibrium with three-dimensional stable manifold and one nonhyperbolic equilibrium, respectively; that is, for , becomes unstable and becomes feasible. Biologically, for spreading of infection into the system, the virus replication factor must be greater than . When , the predator-free equilibrium is feasible and stable. For higher value of , there exists a which deduces that there exists a such that at the system enters into a Hopf bifurcation. When the solution becomes unstable and gives a periodic orbit of small amplitude near . For the positive equilibrium , there exists a , the positive equilibrium is locally asymptotically stable if , and the Hopf bifurcation occurs at , bifurcating a periodic orbit which is orbitally stable applying Poore’s condition for verification. Considering, in practice, that the predator and the virus are thrown at some moments, lastly in Section 6, we think about the corresponding impulsive equation to model (6). The results have shown that when the impulsive releasing period , then the pest-eradication solution of system is globally asymptotically stable; otherwise the system is uniformly permanent. So the pest management can be realized by controlling the impulsive releasing period.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
The authors would like to thank the referee, Professor H. F. Huo, for his careful reading of the original paper and giving meaningful suggestions to improve the quality of the paper. At the same time, they thank Professor S. Y. Tang for his instructive advice and useful suggestions for the numerical simulation. This research was supported by Guangxi Natural Science Foundation 2012GXNSFAA053011 and Foundation of Guangxi Educational Committee (no. YB2014203).
- B. S. Goh, Management and Analysis of Biological Populations, Elsevier Science, Amsterdam, The Netherlands, 1980.
- L. A. Lacey and M. S. Goettel, “Current developments in microbial control of insect pests and prospects for the early 21st century,” Entomophaga, vol. 40, no. 1, pp. 3–27, 1995.