Research Article  Open Access
Bifurcation Analysis in Models for VectorBorne Diseases with Logistic Growth
Abstract
We establish and study vectorborne models with logistic and exponential growth of vector and host populations, respectively. We discuss and analyses the existence and stability of equilibria. The model has backward bifurcation and may have no, one, or two positive equilibria when the basic reproduction number is less than one and one, two, or three endemic equilibria when is greater than one under different conditions. Furthermore, we prove that the diseasefree equilibrium is stable if is less than 1, it is unstable otherwise. At last, by numerical simulation, we find rich dynamical behaviors in the model. By taking the natural death rate of host population as a bifurcation parameter, we find that the system may undergo a backward bifurcation, saddlenode bifurcation, Hopf bifurcation, BogdanovTakens bifurcation, and cusp bifurcation with the saturation parameter varying. The natural death rate of host population is a crucial parameter. If the natural death rate is higher, then the host population and the disease will die out. If it is smaller, then the host and vector population will coexist. If it is middle, the period solution will occur. Thus, with the parameter varying, the disease will spread, occur periodically, and finally become extinct.
1. Introduction
The dynamical modeling and studies of vectorborne disease have been in many literatures. The earliest studies originated from the RossMacdonald model for malaria [1] in 1911 and major extensions are described in Macdonaldâ€™s 1957 book [2]. For malaria, there aremany literature [3â€“7]. Research literatures on dengue fever have [8â€“12], and so forth. For West Nile virus (WNV), there are also many literatures, [13â€“18], and so forth. Here, we classify vector population as susceptible and infected and classify the host population as susceptible , infected , and recovered . In these literatures, birth function of susceptible host and vector population and incidence rate are taken as Table 1. Parameters involved in Table 1 for the dynamics of vectorborne diseases are in Table 2. For malaria, dengue fever, and West Nile virus, the three typical vectorborne diseases, the dynamics of most of the related models can be characterized by the reproduction number, but some compartmental models for these three mosquitoborne diseases undergo a backward bifurcation, an important feature of mosquitoborne diseases.

For compartmental models for malaria, dengue fever, and WNV, different forms of birth functions have been used; we summarize and find that the logistic birth function which describes the selflimiting growth of a host population was rarely considered. In the paper, we will build the following model and analyse the dynamic behaviors: In the model (1), the host and vector populations satisfy the following equations, respectively: It means that the total number of the vector populations is a constant and denotes . Furthermore, using the relations and , we lead to the following fourdimensional nonlinear system of ODEs:
One can verify that the positive cone is positive invariant.
The remainder of this paper is organized as follows. In Section 2, we consider the number and existence of equilibria. Section 3 studies local stability of the diseasefree equilibrium and gives some examples to explain the existence of endemic equilibria and bifurcation with the basic reproduction number changing. The paper concludes with a brief discussion in Section 4.
2. Existence of Equilibria
To study the existence of equilibria, let the right side of each of the four differential equations be equal to zero in system (3), obtaining the following equations:
From the above equations, one can verify that system (3) has a diseasefree equilibrium at . The linear stability of is governed by the basic reproduction number which can be found from the next generation matrix [19] for the system. Noting that the model has infected populations, namely, and , it follows that, using the notation of van den Driessche and Watmough [20], the matrices and , for the new infection terms and the remaining transfer terms, respectively, are given in partitioned form by where Here, is a nonnegative matrix of rank 2 and is nonsingular. It is easy to show that, at steady state, the spectral radius (dominant eigenvalue) of the nonnegative matrix , denoted by , is equal to ; hence, One can see that is the geometric mean of the number of new diseased generated in both the host and vector population, respectively. The first part is the number of host infections caused by one infected vector; the second part is the number of vector infections caused by one infected host. Discussion in detail on basic reproductive number of hostvector disease models see literature [16].
In order to find positive equilibria, it follows from the first and last two equations of (5), we have Substituting (10) into the second equation of (5), that is, if an endemic equilibrium exists, its coordinate satisfies where
should satisfy to ensure a positive equilibrium. Then we will consider the existence of positive roots for (11) from the following two cases.
2.1. Case 1: ,
In this case, system (3) for vectorborne diseases is reduced to For system (13), satisfies where
Now we consider the case . One can see that it is the number of roots in the interval that determines the number of positive equilibria for model (13). Note If , , has a unique root in . If , , then there may be two positive equilibria or no positive equilibrium. We have the following theorem.
Theorem 1. If , system (13) has a unique positive equilibrium .
Theorem 2. If , one has the following.(a)If , then system (13) has no endemic equilibrium.(b)If , then(i)when , system (13) has no endemic equilibrium;(ii)when , system (13) has a unique repeated endemic equilibrium ;(iii)when , system (13) has two endemic equilibria , , where .
2.2. Case 2: ,
For of system (3), If , and , . Then we will analyze the existence of positive roots of . The derivative of is Denote Then we have the following theorem.
Theorem 3. If , one has the following.(a)System (3) has no endemic equilibrium, if any one of the following conditions holds:(1) and ,(2), , and ,(3), , , and .(b)System (3) has two endemic equilibria, if , , , and .
Theorem 4. If , one has the following.(a)System (3) has a unique endemic equilibrium, if any one of the following conditions holds:(1) and ;(2), , ;(3), , , and ;(4), , , , ;(5), , , , , and .(b)System (3) has three endemic equilibria, if , , , , , and .
It is easy to prove Theorem 4, we will omit the proof. In the following part, we will give two examples to show changes of equilibria with the basic reproduction number by numerical simulation when and .
Example 5. Fix parameters , , , , , , , , , , and . Figure 1 shows that there is always a unique endemic equilibrium when . As increases, the equilibrium will be keep. If equals one, a new positive equilibrium is branched from the diseasefree equilibrium. If continue to be increased, the two positive equilibria may be repeated and then disappeared if exceeds a certain value, that is, there are two positive equilibria or no positive equilibrium when if .
Example 6. Take , , , , , , , , , and as same values of Example 5 and . From Figure 2, we find that the number of positive equilibria is from zero to two with increasing when . If is more than one, it is found that the two positive equilibria still be keep and there is a new endemic one to be branched from the diseasefree equilibrium. Now system (3) has three endemic equilibria. Continuing to increase , we find that the larger two equilibria will be repeated and then disappear and the smaller equilibrium always be keep. Let ; from Figure 3, it can be obtained that there is no positive equilibrium when . With increasing, a positive equilibrium is branched from the diseasefree one; then three ones appear and the larger two ones disappear and the smaller one lasts.
3. Stability and Bifurcation Analysis
Then we study the stability of the steady states of system (3) and possible bifurcations. Firstly, we discuss the stability of diseasefree equilibrium . At this equilibrium point, the Jacobian matrix is where The characteristic polynomial of (20) is It means that or
The roots of the two quadratic polynomial have negative real parts if and only if their coefficients are positive; that is, . Therefore, the diseasefree equilibrium is always locally asymptotically stable when , and it is unstable if .
For any positive endemic equilibrium , the Jacobian matrix of (3) becomes where Here, it is difficult to determine a sign for the real part of eigenvalue of , so we will give the dynamical behaviors for system (3) by numerical simulation.
Example 7. Fix parameters , , , , , , , , , , and . In Figure 1, we find that a unique positive equilibrium is unstable when . With increasing and exceeding one, the unstable positive equilibrium still be keep and a new positive one branching from the diseasefree one is stable. After is more than a threshold, the two equilibria will disappear, and all of trajectories tend to equilibrium , which means that the host population will be extinct If we fix and let change, we find that Hopf bifurcation maybe occurs from the blue curve of Figure 4 for system (13).
Example 8. Fix parameters , , , , , , , , , and . Then we draw the bifurcation figure with changed when , respectively. We use blue, black, and brown curves to represent the bifurcation curves when . In Figure 4, we firstly choose as bifurcation parameter and draw the bifurcation curve (blue (), black (), and brown curves ()) with . When , we find that is always more than no matter what the value of is, which means that system (3) has no or two positive equilibria. From the blue curve of Figure 4, we can see that there are a limit point (LP) at and a Hopf point (H) at whose first Lyapunov coefficient is . When , it is found that two limit points at and and a Hopf point at whose first Lyapunov coefficient is for system (3). This is consistent with the existence of the positive equilibrium. If we take , by numerical simulation, we see that there are two limit points at and , two neutral saddle points at and , and two Hopf points at (first Lyapunov coefficient ) and (first Lyapunov coefficient ). Then take a limit point as initial point and and as bifurcation parameters and obtain a fold curve (red curve in Figure 4). In the fold curve, there are a cusp bifurcation point at , , and two BogdanovTakens bifurcation points at , and , . At the cusp bifurcation point two branches of saddlenode bifurcation curve meet tangentially, forming a semicubic parabola. For nearby parameter values, the system has three equilibria which collide and disappear pairwise via the saddlenode bifurcations. For the BogdanovTakens (BT) bifurcation nearby parameter values, the system has two equilibria which collide and disappear via a saddlenode bifurcation. In the end, take any one of BogdanovTakens bifurcation points as initial point and draw Hopf bifurcation curve (the green curve of Figure 4). It shows that there is a ZeroNeutral Saddle (ZH) at and . Thus, we can obtain that all of points being saddle between two BogdanovTakens bifurcation points and Hopf points locate on both sides of BogdanovTakens bifurcation points from Hopf bifurcation curve of Figure 4.
In the following example, we will analyse how to change the state variables and with time .
Example 9. Take , , , ,, ,,, and as the same with the above example and . We see how to change for trajectories and with time when , and , respectively. From Figures 5 and 9, we can find that all of trajectories tend to an endemic equilibrium when is smaller, and all of trajectories tend to an equilibrium when is larger, that is, the host population will be distinct, and the vector population will be keep. If , respectively, from Figures 6, 7, and 8, we find that system (3) has period solution. Furthermore, from these figures, we can see that the numbers of the infective host population and the infective vector population are on the decline with the natural death rate increase, but the rate of the decline for host population is much faster than that for vector population. Finally, the host population decays to zero.
4. Discussion
In previous literatures, the dynamic behaviors of models on vectorborne disease may have backward bifurcation and are not complex. In this paper, we consider the behaviors of a model with standard incidence rate and the logistic birth function for host population and the exponential birth rate for vector population. The rich dynamic behaviors are found. We find that the complex behaviors will not occur if the birth function is not logistic. Taking the natural death rate of host population as a bifurcation parameter, we find that the system may undergo a backward bifurcation, saddlenode bifurcation, Hopf bifurcation, BogdanovTakens bifurcation, and cusp bifurcation with the saturation parameter varying. We find that the natural death rate still is a key factor. If the natural death rate is higher, then the host population and the disease will die out; if it is smaller, then the host and vector population will coexist; that is, the disease will spread. If it is middle, the period solution will occur. In addition, the existence of equilibria of our model is different from previous model. If the incidence rate function is not affected by the number of vector populations, we find that there is always a unique positive equilibrium when the basic reproduction number is less than one, and there are two endemic equilibria or no positive equilibrium when is more than one. If the incidence rate function depends on the number of vector populations, we obtained that there may be two positive equilibria or no positive equilibrium when and that there exist one or three endemic equilibria when where the unique endemic one is always unstable.
Conflict of Interests
The authors declare that they have no conflict of interests regarding the publication of this paper.
Acknowledgments
This work is supported by the National Science Foundation of China (11201434) and Shanxi Scholarship Council of China (2013087).
References
 R. Ross, The Prevention of Malaria, John Murray, London, UK, 1911.
 G. Macdonald, The Epidemiology and Control of Malaria, Oxford University Press, London, UK, 1957.
 G. A. Ngwa and W. S. Shu, â€śA mathematical model for endemic malaria with variable human and mosquito populations,â€ť Mathematical and Computer Modelling, vol. 32, no. 78, pp. 747â€“763, 2000. View at: Publisher Site  Google Scholar
 G. A. Ngwa, â€śModelling the dynamics of endemic malaria in growing populations,â€ť Discrete and Continuous Dynamical Systems B, vol. 4, no. 4, pp. 1173â€“1202, 2004. View at: Google Scholar
 N. Chitnis, Using mathematical models in controlling the spread of malaria [Ph.D. thesis], Program in Applied Mathematics, University of Arizona, Tucson, Ariz, USA, 2005.
 N. Chitnis, J. M. Cushing, and J. M. Hyman, â€śBifurcation analysis of a mathematical model for malaria transmission,â€ť SIAM Journal on Applied Mathematics, vol. 67, no. 1, pp. 24â€“45, 2006. View at: Publisher Site  Google Scholar
 J. Tumwiine, J. Y. T. Mugisha, and L. S. Luboobi, â€śA hostvector model for malaria with infective immigrants,â€ť Journal of Mathematical Analysis and Applications, vol. 361, no. 1, pp. 139â€“149, 2010. View at: Publisher Site  Google Scholar
 L. Esteva and C. Vargas, â€śAnalysis of a dengue disease transmission model,â€ť Mathematical Biosciences, vol. 150, no. 2, pp. 131â€“151, 1998. View at: Publisher Site  Google Scholar
 L. Esteva and C. Vargas, â€śA model for dengue disease with variable human population,â€ť Journal of Mathematical Biology, vol. 38, no. 3, pp. 220â€“240, 1999. View at: Google Scholar
 L. Esteva and C. Vargas, â€śInfluence of vertical and mechanical transmission on the dynamics of dengue disease,â€ť Mathematical Biosciences, vol. 167, no. 1, pp. 51â€“64, 2000. View at: Publisher Site  Google Scholar
 D. J. Gubler, â€śDengue,â€ť in The Arboviruses: Epidemiology and Ecology, vol. 2, pp. 223â€“261, CRC Press, New York, NY, USA, 1986. View at: Google Scholar
 Z. Feng and J. X. VelascoHernández, â€śCompetitive exclusion in a vectorhost model for the dengue fever,â€ť Journal of Mathematical Biology, vol. 35, no. 5, pp. 523â€“544, 1997. View at: Google Scholar
 M. J. Wonham, T. deCaminoBeck, and M. A. Lewis, â€śAn epidemiological model for West Nile virus: invasion analysis and control applications,â€ť Proceedings of the Royal Society B, vol. 271, no. 1538, pp. 501â€“507, 2004. View at: Publisher Site  Google Scholar
 C. Bowman, A. B. Gumel, P. van den Driessche, J. Wu, and H. Zhu, â€śA mathematical model for assessing control strategies against West Nile virus,â€ť Bulletin of Mathematical Biology, vol. 67, no. 5, pp. 1107â€“1133, 2005. View at: Publisher Site  Google Scholar
 G. CruzPacheco, L. Esteva, J. A. MontañoHirose, and C. Vargas, â€śModelling the dynamics of West Nile virus,â€ť Bulletin of Mathematical Biology, vol. 67, no. 6, pp. 1157â€“1172, 2005. View at: Publisher Site  Google Scholar
 M. J. Wonham, M. A. Lewis, J. Rencławowicz, and P. van den Driessche, â€śTransmission assumptions generate conflicting predictions in hostvector disease models: a case study in West Nile virus,â€ť Ecology Letters, vol. 9, no. 6, pp. 706â€“725, 2006. View at: Publisher Site  Google Scholar
 J. Jiang, Z. Qiu, J. Wu, and H. Zhu, â€śThreshold conditions for West Nile virus outbreaks,â€ť Bulletin of Mathematical Biology, vol. 71, no. 3, pp. 627â€“647, 2009. View at: Publisher Site  Google Scholar
 H. Wan and H. Zhu, â€śThe backward bifurcation in compartmental models for West Nile virus,â€ť Mathematical Biosciences, vol. 227, no. 1, pp. 20â€“28, 2010. View at: Publisher Site  Google Scholar
 O. Diekmann, J. A. Heesterbeek, and J. A. Metz, â€śOn the definition and the computation of the basic reproduction ratio ${R}_{0}$ in models for infectious diseases in heterogeneous populations,â€ť Journal of Mathematical Biology, vol. 28, no. 4, pp. 365â€“382, 1990. View at: Google Scholar
 P. van den Driessche and J. Watmough, â€śReproduction numbers and subthreshold endemic equilibria for compartmental models of disease transmission,â€ť Mathematical Biosciences, vol. 180, pp. 29â€“48, 2002. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2014 Guihua Li and Zhen Jin. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.