The Scientific World Journal

Volume 2014 (2014), Article ID 195864, 9 pages

http://dx.doi.org/10.1155/2014/195864

Research Article

Bifurcation Analysis in Models for Vector-Borne Diseases with Logistic Growth

Guihua Li and Zhen Jin

Department of Mathematics, North University of China, Taiyuan, Shanxi 030051, China

Received 21 November 2013; Accepted 10 February 2014; Published 25 March 2014

Academic Editors: H.-F. Huo and P. G. L. Leach

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.

Abstract.

We establish and study vector-borne 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 disease-free 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, saddle-node bifurcation, Hopf bifurcation, Bogdanov-Takens 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 vector-borne disease have been in many literatures. The earliest studies originated from the Ross-Macdonald 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 vector-borne diseases are in Table 2. For malaria, dengue fever, and West Nile virus, the three typical vector-borne diseases, the dynamics of most of the related models can be characterized by the reproduction number, but some compartmental models for these three mosquito-borne diseases undergo a backward bifurcation, an important feature of mosquito-borne diseases.

Table 1: The birth and incidence rates for different vector-borne disease models.

Disease
Literature
Birth
Incidence rate

Host
Vector

[3]

Malaria
[5, 6]

[7]

[8]

Dengue fever
[9]

[10]

[13]

WNV
[14]

[15, 18]

Table 2: Parameters involved in the models for the dynamics of vector-borne diseases.

Parameters in the models
Vector
Host

Per capita birth rate

The fraction of the infective immigrants

Number of times one mosquito would want to bite humans per unit time

The maximum number of mosquito bites a human can have per unit time

Recruitment rate per unit of time

Per capita birth rate

Vertical transmission rate in new birth of vectors

Natural death rate

Disease-induced death rate

Recovery rate

Biting rate (the average number of bites per mosquito per day)

Transmission probability (from infected vectors to host)

Transmission probability (from infected host to vectors)

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 self-limiting 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 four-dimensional 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 disease-free 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 disease-free 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 host-vector 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 vector-borne 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 disease-free 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 .

Figure 1: Bifurcation curves in plane when , , , , , , , , , , , and .

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 disease-free 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 disease-free one; then three ones appear and the larger two ones disappear and the smaller one lasts.

Figure 2: Bifurcation curves in plane when , , , , , , , , , , , and .

Figure 3: Bifurcation curves in plane when , , , , , , , , , , , and .

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 disease-free 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 disease-free 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 disease-free 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).

Figure 4: Bifurcation curves in plane when , , , , , , , , , , , and .

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 Bogdanov-Takens bifurcation points at , and , . At the cusp bifurcation point two branches of saddle-node bifurcation curve meet tangentially, forming a semicubic parabola. For nearby parameter values, the system has three equilibria which collide and disappear pairwise via the saddle-node bifurcations. For the Bogdanov-Takens (BT) bifurcation nearby parameter values, the system has two equilibria which collide and disappear via a saddle-node bifurcation. In the end, take any one of Bogdanov-Takens bifurcation points as initial point and draw Hopf bifurcation curve (the green curve of Figure 4). It shows that there is a Zero-Neutral Saddle (ZH) at and . Thus, we can obtain that all of points being saddle between two Bogdanov-Takens bifurcation points and Hopf points locate on both sides of Bogdanov-Takens 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.

Figure 5: Trajectory picture with time when , , , , , , , , , , , and .

Figure 6: Trajectory picture with time when , , , , , , , , , , , and .

Figure 7: Trajectory picture with time when , , , , , , , , , , , and .

Figure 8: Trajectory picture with time when , , , , , , , , , , , and .

Figure 9: Trajectory picture with time when , , , , , , , , , , , and .

4. Discussion

In previous literatures, the dynamic behaviors of models on vector-borne 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, saddle-node bifurcation, Hopf bifurcation, Bogdanov-Takens 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 (2013-087).

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. 7-8, pp. 747–763, 2000.

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.

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.

J. Tumwiine, J. Y. T. Mugisha, and L. S. Luboobi, “A host-vector model for malaria with infective immigrants,” Journal of Mathematical Analysis and Applications, vol. 361, no. 1, pp. 139–149, 2010.

L. Esteva and C. Vargas, “Analysis of a dengue disease transmission model,” Mathematical Biosciences, vol. 150, no. 2, pp. 131–151, 1998.

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.

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.

D. J. Gubler, “Dengue,” in The Arboviruses: Epidemiology and Ecology, vol. 2, pp. 223–261, CRC Press, New York, NY, USA, 1986.

Z. Feng and J. X. Velasco-Hernández, “Competitive exclusion in a vector-host model for the dengue fever,” Journal of Mathematical Biology, vol. 35, no. 5, pp. 523–544, 1997.

M. J. Wonham, T. de-Camino-Beck, 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.

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.

G. Cruz-Pacheco, L. Esteva, J. A. Montaño-Hirose, and C. Vargas, “Modelling the dynamics of West Nile virus,” Bulletin of Mathematical Biology, vol. 67, no. 6, pp. 1157–1172, 2005.

M. J. Wonham, M. A. Lewis, J. Rencławowicz, and P. van den Driessche, “Transmission assumptions generate conflicting predictions in host-vector disease models: a case study in West Nile virus,” Ecology Letters, vol. 9, no. 6, pp. 706–725, 2006.

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.

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.

O. Diekmann, J. A. Heesterbeek, and J. A. Metz, “On the definition and the computation of the basic reproduction ratio in models for infectious diseases in heterogeneous populations,” Journal of Mathematical Biology, vol. 28, no. 4, pp. 365–382, 1990.

P. van den Driessche and J. Watmough, “Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission,” Mathematical Biosciences, vol. 180, pp. 29–48, 2002.