Table of Contents Author Guidelines Submit a Manuscript
Journal of Applied Mathematics
Volume 2014 (2014), Article ID 324767, 16 pages
Research Article

Mathematical Analysis of a Cholera Model with Vaccination

1School of Science, Beijing University of Civil Engineering and Architecture, Beijing 100044, China
2College of Mathematics and Information Science, Xinyang Normal University, Xinyang, Henan 464000, China

Received 14 March 2013; Revised 18 May 2013; Accepted 26 September 2013; Published 13 February 2014

Academic Editor: Tin-Tai Chow

Copyright © 2014 Jing'an Cui et al. 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.


We consider a SVR-B cholera model with imperfect vaccination. By analyzing the corresponding characteristic equations, the local stability of a disease-free equilibrium and an endemic equilibrium is established. We calculate the certain threshold known as the control reproduction number . If , we obtain sufficient conditions for the global asymptotic stability of the disease-free equilibrium; the diseases will be eliminated from the community. By comparison of arguments, it is proved that if , the disease persists and the unique endemic equilibrium is globally asymptotically stable, which is obtained by the second compound matrix techniques and autonomous convergence theorems. We perform sensitivity analysis of on the parameters in order to determine their relative importance to disease transmission and show that an imperfect vaccine is always beneficial in reducing disease spread within the community.

1. Introduction

Cholera is an acute intestinal infection caused by the ingestion of food or water contaminated with the bacterium Vibrio cholera. Among the 200 serogroups of Vibrio cholera, it is only Vibrio cholera o1 and o139 that are known to be the cause of the cholera disease [1]. The etiological agent, Vibrio cholera o1 (and more recently vibrio cholera o139), passes through and survives the gastric acid barrier of the stomach and then penetrates the mucus lining that coats the intestinal epithelial [2]. Once they colonise the intestinal gut, then produce enterotoxin (which stimulates water and electrolyte secretion by the endothelial cells of the small intestine) that leads to copious, painless, and watery diarrhoea that can quickly lead to severe dehydration and death if treatment is not promptly given. Vomiting also occurs in most patients. In human volunteer studies, the infection was determined to be 102–103 [3]. Cholera can either be transmitted through interaction between humans (i.e., fecal-oral) or through interaction between humans and their environment (i.e., ingestion of contaminated water and food from the environment). To come on urgent, transmission fast, sweeping range widely are the characteristics of cholera which is one of the international quarantine infectious diseases as stipulated by the International Health Regulations (IHR), as well as one of a class of infectious diseases as stipulated by law for infectious diseases prevention and control of China.

Globally, cholera incidence has increased steadily since 2005 with cholera outbreaks affecting several continents (see Figure 1). Cholera continues to pose a serious public health problem among developing world populations which have no access to adequate water and sanitation resources. In 2011, 32% of cases were reported from Africa whereas between 2001 and 2009, 93% to 98% of total cases worldwide were reported from that continent [4]. In 2011, 61.2% of cases were reported from Americas where a large outbreak that started in Haiti at the end of October 2010 also affected the Dominican Republic. The outbreak was still ongoing at the end of 2011 with 522335 cases including 7001 deaths that were reported by 25 December in Haiti [4]. So, cholera remains a global threat and is one of the key indicators of social development. The history and reality have warned: we are facing the serious threat of cholera, the importance for the study of cholera’s pathogenesis, regular transmission and prevention and control strategy have become increasingly prominent, which has also become a major problem that needs to be solved. Up till now, a number of mathematical models have been used to study the transmission dynamics of cholera. Capasso and Serio [5] introduced an incidence rate in the form of (with human-to-human transmission model only) in 1973 [6]. Codeço [7] proposed an incidence form of (with environment-to-human transmission model only) in 2001 which, in the first time, explicitly incorporated the pathogen concentration into cholera modeling. Mukandavire et al. [8] included both transmission pathways in the form of . In 2012, Liao and Wang [9] generalised Codeco’s model [7], incorporating the theory of Volterra-Lyapunov stable matrices into the classic method of Lyapunov functions; they studied the global dynamics of the mathematical model.

Figure 1: In 1989–2011, cholera cases were reported to World Health Organization by year and by continent [4].

Vaccination is a major factor in the resurgence and epidemic outbreaks of some infectious diseases. Since the pioneering work of Edward Jenner on smallpox [10], vaccination has been a commonly used method for diseases control [1113] and works by reducing the number of susceptible individuals in a population. In modern times, vaccination has a large impact on the incidence and persistence of children infections, such as measles and whooping cough [14]. Hethcote [11] investigated a pertussis infectious model and showed that the vaccination can make the infection undergo fluctuation. Although vaccination offers a very powerful tool for disease control, generally, vaccines are not 100% affective and sometimes they only provide limited immunity due to the natural waning of immunity in the host or antigenic variation in the pathogen [15].

Here, we develop a cholera model with an additional equation for the vaccinated individuals in the population. Since Koch found Vibrio cholera in 1883, the research for cholera vaccine had been going on for over one hundred years. People have developed a variety of vaccines. However, these vaccines were parenteral, which have short effective protection and big side effects. In 1973, the World Health Organization canceled the vaccine inoculation which attracted a major concern to oral vaccines. At present, there are three kinds of oral vaccines (i.e., WC/BS vaccine, WC/rBS vaccine, and CVD103-HgR vaccine) have been proved to be safe, effective, and immunogenic, which were approved to apply in some countries [16].

In this paper, according to the natural history of cholera, we improve the model of [9] in the following two aspects. Firstly, if the cholera persists for a long time, it will cause death [17], especially in the area where water and sanitation resources are not adequate [4]; a parameter is added to describe the rate of disease-related death. Secondly, we propose a proportion of the vaccination in susceptible individuals as shown in the following differential equations:

The flow diagram of the model is depicted in Figure 2. Since the first three and last equations in (1) are independent of the variable , it suffices to consider the following reduced model:

Figure 2: Progression of infection from susceptible () and vaccinated () individuals through the infected () and recovered () compartments for the combined human-environment epidemiological model with an environmental component.

Here, , , , and refer to the susceptible individuals, infected individuals, vaccinated individuals, and recovered individuals, respectively. The pathogen population at time is given by . The parameter denotes the natural human birth and death rate, denotes the rate of recovery from the disease, represents the rate of human contribution to the growth of the pathogen, and represents the death rate of the pathogen in the environment. The coefficients and represent the contact rates for the human-environment and human-human interactions, respectively. Constants and adjust the appropriate form of the incidence which determines the rate of new infection. If , the corresponding incidence is reduced to the standard bilinear form based on the mass action law, which is most common in epidemiological models. If , then the corresponding incidence represents a consequence of saturation effects: when the infected number is high, the incidence rate will respond more slowly than linearly to the increase in . Similar meanings stand for . The rate at which the susceptible population is vaccinated is , and the rate at which the vaccine wears off is . All parameters are assumed nonnegative.

The organization of this paper is as follows: the positivity and boundedness of solutions are obtained in Section 2. In Section 3, we obtain the existence of the endemic equilibrium. We get the local and global stability of the disease-free equilibrium in Section 4. In Section 5, we present the persistence of the system. In Section 6, we show the local and global stability of the endemic equilibrium. We analyze the sensitivity of on the parameters, and we present the numerical simulation in Section 7. The paper ends with a conclusion in Section 8.

2. Positivity and Boundedness of Solutions

In the following, we show that the solutions of system (2) are positive with the nonnegative initial conditions.

Theorem 1. The solutions of model (2) are nonnegative for all with the non-negative initial conditions.

Proof. System (2) can be put into the matrix form where and is given by We have Therefore, Due to Lemma 2 in [18], any solution of system (2) is such that for all . This completes the proof of Theorem 1.

Theorem 2. All solutions of model (2) are bounded.

Proof. System (2) is split into two parts, the human population (i.e., , , and ) and pathogen population (i.e., ). It follows from the first three equations of system (2) that then it follows that . From the first equation, we can get Thus , as . It is easy to obtain Thus , as . From the last equation, we can obtain Hence, , when . Therefore, all solutions of model (2) are bounded.

From above discussion, we can see that the feasible region of human population for system (2) is and the feasible region of pathogen population for system (2) is

Define . Let denote the interior of . It is easy to verify that the region is a positively invariant region (i.e., the solutions with initial conditions in remain in ) with respect to system (2). Hence, we will consider the global stability of (2) in region .

3. Equilibria

In this section, we investigate the existence of equilibria of system (2). Solving the right hand side of model system (2) by equating it to zero, we obtain the following biologically relevant equilibria.

It is easy to see that model (2) always has a disease-free equilibrium (the absence of infection, i.e., ) , where and .

In the following, we will discuss the existence and uniqueness of the endemic equilibrium. The components of the endemic equilibrium satisfy

Substituting (13a), (13b), and (13d) into (13c), we obtain a single equation for : After dropping the solution , we obtain where Note that represents a straight line with a negative slope and a vertical intercept . Meanwhile, we have We see that is increasing for , and , where is the control reproduction number of infection. When . Hence, there is one and only one intersection between the curves of and ; that is, there is a unique solution to the equation . Consequently, , , and are uniquely determined by .

Theorem 3. System (2) has a unique endemic equilibrium when and no positive endemic equilibrium when .

4. Stability of Disease-Free Equilibrium

Now, we will discuss the local and global stability of the disease-free equilibrium.

Theorem 4. The disease-free equilibrium is locally asymptotically stable for and unstable for .

Proof. The Jacobian matrix of system (2) at is The characteristic polynomial of the matrix is given by where If , then Further, So , , , , , and (see Appendix A). Thus, using the Routh-Hurwitz criterion, all eigenvalues of have negative real part; is local asymptotically stable for system (2). If , then and we show that has at least one eigenvalue with nonnegative real parts. Consequently, is not asymptotically stable.

Theorem 5. When , the disease-free equilibrium is globally asymptotically stable.

We will prove the global stability of the disease-free equilibrium using Lemma 6.

Lemma 6 (see [19]). If a model system can be written in the form where denotes (its components) the number of uninfected individuals and denotes (its components) the number of infected individuals including latent and so forth, denotes the disease-free equilibrium of the system.
Assume that(H1)for , is globally asymptotically stable;(H2), for , where the Jacobian matrix is an Metzler matrix (the off-diagonal elements of are nonnegative) and is the region where the model makes biological sense. Then the fixed point is a globally asymptotically stable equilibrium of cholera model system (2) provided that .
We begin by showing condition (H1) as For the equilibrium , the system reduces to The characteristic polynomial of the system is given by There are two negative characteristic foots: . Hence, is always globally asymptotically stable.
Next, applying Lemma 6 to the cholera model system (2) gives So is a Metzler matrix. Meanwhile, we find . Hence, the disease-free equilibrium is globally asymptotically stable.

5. Persistence

Persistence is an important property of dynamical systems and of the systems in ecology, epidemics, and so forth, they are modeling. Biologically, persistence means the survival of all populations in future time. Mathematically, persistence of a system means that strictly positive solutions do not have any omega limit points on the boundary of the nonnegative cone [20]. In this section, we will present the persistence of system (2).

For various definitions of persistence [21, 22], we utilize the definitions of persistence developed by Freedman et al. [23]. System (2) can be defined to be uniformly persistent if for some for all initial points in .

A uniform persistence result given in [23] requires the following hypothesis (H) to be satisfied.

We denote that is a closed positively invariant subset of on which a continuous flow is defined and is the maximal invariant set of on . Suppose is a closed invariant set and there exists a cover of , where is a nonempty index set; , , and are pairwise disjoint closed invariant sets. Furthermore, we propose the following hypothesis.

Hypothesis (H):(a)all are isolated invariant sets of the flow ;(b) is acyclic; that is, any finite subset of does not form a cycle [24];(c)any compact subset of contains, at most, finitely many sets of .

Lemma 7 (see [24]). Let be a closed positively invariant subset of on which a continuous flow is defined. Suppose there is a constant such that is point dissipative on and the assumption (H) holds. Then the flow is uniformly persistent if and only if for any , where .

Now, we can obtain the following result.

Theorem 8. System (2) is uniformly persistent in if .

Proof. Suppose . We show that system (2) satisfies all the conditions of Lemma 7. Choose and . The vector field of system (2) is transversal to the boundary of on its faces except the S-axis and V-axis, which are invariant with respect to system (2) and on the S-axis and V-axis the equations for and are and , which implies that and as . Therefore, is the only -limit point on the boundary of . As the maximal invariant set on the boundary of is the singleton and is isolated when , thus the hypothesis (H) holds for system (2). The flow induced by is point dissipative by the positive invariance of . Because of , where is the omega limit set of , when , we have that is contained in the set and for , . Therefore, the uniform persistence of system (2) is equivalent to being unstable, and the theorem is proved.

Remark 9. Theorems 3 and 8 show that is a threshold parameter for the model; that is, when , its epidemiological implication is that the infected fraction of the population vanishes, so the cholera dies out; when , the disease is endemic and the infected fraction remains above a certain positive level for sufficiently large time.

6. Stability of the Endemic Equilibrium

6.1. Local Stability of the Endemic Equilibrium

Now we consider the case with . The stability of the endemic equilibrium is established as follows.

Theorem 10. If , is locally asymptotically stable.

Proof. Let The Jacobian matrix at is The characteristic polynomial of the matrix is given by where Based on (13c) and (13d), we have . It is then easy to observe that Further, So , , , , , and (see Appendix B). Using the well-known Routh-Hurwitz criterion, the proof is thus complete.

6.2. Global Stability of the Endemic Equilibrium

Theorem 11. When , is globally asymptotically stable in if and .

This approach to global dynamics is developed in the papers of Smith [25] and Li and Muldowney [2628]. Let , where ,,, and represent the right-hand sides of system (2), respectively. Furthermore, let . Then, the Jacobian matrix for system (2) is

The second additive compound matrix (see Appendix C for details) of iswhere

Let Then we havewhere In (41), is the directional derivative of in the direction of the vector field in system (2).

Here, we will use the theorem in [28] to give a sufficient condition on the parameters, which when satisfied, implies that the endemic equilibrium is globally asymptotically stable.

Lemma 12 (see [28]). If is a compact absorbing subset in the interior of and there exists and a Lozinskiǐ measure for all , then every omega point of system (2) in is an equilibrium in .

For , the disease-free equilibrium is repelling towards the interior. In fact, there is a compact absorbing set in int which attracts all orbits that intersect int. This gives the following results.

Corollary 13. If and there exists a Lozinskiǐ measure such that for all , then each orbit of system (2) which intersects limits to the endemic equilibrium.

For a norm on , the Lozinskiǐ measure associated with can be evaluated for an matrix as where is the right-hand derivative [29]. Hence, if we can find a norm on for which the associated Lozinskiǐ measure satisfied for all then the endemic equilibrium is globally asymptotically stable for .

We now define a norm on [30] for which the definition varies from one orthant to another. Let where , with components and

Theorem 14. Assume that Then there exits , such that for all .

Proof. We should show that

Case 1 ( and ). In this case, .

Subcase 1 (). Then and . Taking the right-hand derivative of , we obtain Since , we have Hence,
It is easy to see that (50) also holds for and when , which can be obtained by linearity.

Subcase 2 (). Then and .
Since , taking the right-hand derivative of , we obtain Therefore, It is easy to see that (52) also holds for and when , which can be obtained by linearity.

Case 2 ( and ). In this case, .

Subcase 3 (). Then and . Taking the right-hand derivative of , we obtain Therefore, It is easy to see that (54) also holds for and when , which can be obtained by linearity.

Subcase 4 (). Then and . Taking the right-hand derivative of , since , we obtain Therefore,

It is easy to see that (56) also holds for and when , which can be obtained by linearity.

Case 3 ( and ). In this case, .

Subcase 5 (). Then and . Taking the right-hand derivative of , we obtain Recalling that yields It is easy to see that (58) also holds for and when , which can be obtained by linearity.

Subcase 6 (). Then and . Taking the right-hand derivative of , we obtain Therefore, It is easy to see that (60) also holds for and when , which can be obtained by linearity.

Combing the results of the six cases presented here in (50)–(60), as well as the remaining 10 cases, we obtain the result when . Hence, . Therefore, in int  if and . Thus, Theorem 14 is satisfied. From Corollary 13, we can obtain that all solutions intersect the interior of and the endemic equilibrium is globally asymptotically stable, completing the proof of Theorem 11.

7. Sensitivity Analysis of

To facilitate the interpretation of the sensitivity of , we now present some numerical simulations by using the set of parameters values in Table 1.

Table 1: Estimation of parameters.

Now, we regard the vaccinated rate and the waning rate as the control parameters, while the other parameters are fixed. From Figures 3 and 4, the effects of various parameters, that is, and on the control reproduction number have been shown. It is noted that as the parameter increases, decreases and as decreases, decreases. In fact, we can obtain the critical values of and that reduce to 1: In Figure 3, we select , , , , corresponding , ,,, respectively. We can see that when the wanning rate has a greater value, then there is no such that . Similarly, in Figure 4, we select , , , , , corresponding , , , , , respectively. We can see that when is smaller, then there is no such that . Thus, the control reproduction number cannot reduce below unity only by increasing or decreasing . The critical values and are important in regulating the infection magnitude. In order to reduce to 1, a greater vaccinated rate than and a smaller wanning rate than have to be achieved simultaneously. We will deduce below 1 by using both and at the same time, which can control cholera. Otherwise, the cholera persists (see Figure 5).

Figure 3: The contour diagram of the control reproduction number with when has some fixed value.
Figure 4: The contour diagram of the control reproduction number with when has some fixed value.
Figure 5: The contour diagram of the control reproduction number with , variables. All the other parameter values are the same as those in Figure 3.

8. Conclusion

In this paper, we have conducted global stability analysis of a SVIR-B cholera model. Based on the imperfective vaccine, with the environment component incorporated and multiple transmission pathways coupled, the cholera models distinguish themselves from regular SIR epidemiological models. The mathematical results show that the control reproduction number satisfies a threshold property with threshold value 1. When , it has been proved that the disease-free equilibrium is globally asymptotically stable under some sufficient conditions. And, when , the unique endemic equilibrium is globally asymptotically stable. This shows that cholera can be eliminated from the community if the imperfect vaccine brings to a value less than unity.

Now, if we consider that the case there is no vaccination. Model (1) can be written as follows:

According to Theorem 2 in [31], the basic reproduction number of model is . We can express the control reproduction number as . Note that with equality only if . That is, despite being imperfect, the vaccine will always reduce the reproduction number of the disease, and in the absence of the vaccination, the disease transmission will be high. Further simulations in Section 7 also show that the vaccinated rate () and the wanning rate () play equal important roles in reducing the control reproduction number. The disease can be controlled if and only if the reproduction number is reduced to values less than unity if the vaccinated rate () exceed the threshold and the wanning rate () less than the threshold simultaneously. Therefore, vaccination is a good way to control cholera.

However, there are inherent disadvantages towards the vaccination modeling. For cholera with incubation period, it is hard to rapidly identify those with ambiguous symptoms [4]. Moreover, the vaccination does not always work well due to the limitations of medical development level and financial budget (some vaccines are very expensive and some portions of people cannot be covered) [32]. Nonetheless, [33] indicated that during cholera outbreaks (periodic control) vaccination campaigns can be a good strategy to control cholera epidemics. Besides, they pointed out that vaccination and improvement in the sanitation system and food/personal hygiene are the most efficient control strategies to prevent cholera transmission and outbreaks. Hence, incorporating some other control strategies, we may consider the more realistic ordinary differential equation model. The theoretical study of cholera models has been in progress and is an exciting area of future research.


A. The Proof of

The proof of =  +  +  + +