 Research note
 Open Access
 Published:
Modelling cholera transmission dynamics in the presence of limited resources
BMC Research Notesvolume 12, Article number: 475 (2019)
Abstract
Objectives
We study the transmission dynamics of cholera in the presence of limited resources, a common feature of the developing world. The model is used to gain insight into the impact of available resources of the health care system on the spread and control of the disease. A deterministic model that includes a nonlinear recovery rate is formulated and rigorously analyzed. Limited treatment is described by inclusion of a special treatment function. Center manifold theory is used to show that the model exhibits the phenomenon of backward bifurcation. Matlab has been used to carry out numerical simulations to support theoretical findings.
Results
The model analysis shows that the disease free steady state is locally stable when the threshold \({\mathcal {R}}_{0} < 1\). It is also shown that the model has multiple equilibria and the model exhibits the phenomenon of backward bifurcation whose implications to cholera infection are discussed. The results are useful for the public health planning in resource allocation for the control of cholera transmission.
Introduction
Cholera is an acute gastrointestinal infection and waterborne disease which is caused by the bacterium Vibrio Cholerae, V. cholerae \(\text {O}1\) or \(\text {O}139\). Vomiting and diarrhoea are its major characteristics and when patients are treated with delay, it can lead to severe dehydration and death within few hours. The disease has two modes of transmission: direct and the indirect transmission. Direct transmission (human–human) is very uncommon as compared to the indirect (environment–human) which occurs by ingesting contaminated food or water [1, 2]. An estimated 100,000–120,000 deaths are due to cholera every year in the world with only a small proportion being reported to World Health Organization (WHO) [3].
Cholera remains a significant threat to public health in the developing world, with cyclic outbreaks occurring twice per year in endemic areas [4]. For instance, more recently on the 6th of September 2018, a cholera outbreak in Harare was declared by the Ministry of Health and Child Care (MoHCC) of Zimbabwe [5]. As of 15 September 2018, 3621 cumulative suspected cases, including 71 confirmed cases, and 32 deaths had been reported (case fatality ratio: \(0.8 \%\)); of these, \(98\%\) (3564 cases) were reported from the densely populated capital Harare [5]. The City of Harare is facing a plethora of challenges, notably insufficient safe water supplies, frequent sewer pipe bursts, uncollected refuse and rampant illegal vending [6]. This has negatively impacted on public health in the city exposing residents to diarrhoeal disease outbreaks, an upsurge in typhoid fever cases and sporadic outbreaks of cholera [6]. As of the year 2018, the disease has also claimed more than 67 lives in Zambia and Malawi did report some cases of cholera in Lilongwe. Thus, the cholera tragedy continues to devastate disadvantaged countries and communities. For more information about cholera and its occurrence, we refer the reader to [7, 8].
Several mathematical models describing cholera dynamics have been proposed and analyzed; see for instance [8,9,10,11,12,13,14,15,16,17]. These models differ from each other in some aspects. In most of these models, the recovery rate is assumed to be a constant. However, in reality the recovery rate depends on time of recovering process, which can be related to the number of infectious individuals seeking treatment and the basic factor; availability of health resources to the public. The resources of the health system includes the number of health care workers (physicians, nurses, pharmacists, etc.), capacity of the hospital settings (number of hospital beds and medicines) and the effectiveness and efficiency of the treatment. In many developing countries, resources of treatment are extremely scarce. So this approximation cannot reflect the real cure rate.
Motivated by the recent cholera outbreak in Zimbabwe which has overwhelmed the resources of the health care system, we formulate a compartmental deterministic mathematical model with a suitable treatment function in order to study the impact of limited hospital resource capacity on Cholera disease. The number of available hospital beds per 10,000 (hospital bedpopulation ratio) is used by health planners as a method of estimating resource availability to the public [18, 19]. Due to the significance of hospital bedpopulation ratio (HBPR), we will formulate the recovery rate incorporating the impact of limited resource capacity of the health care system in terms of HBPR for this work.
The paper is arranged as follows; in “Main text” section, we formulate and establish the basic properties of the model. The model is analysed for stability in “Analysis of the model” section. In “Numerical simulations” section, we carry out some numerical simulations. Parameter estimation and numerical results are also presented in this section. The paper is concluded in “Conclusions” section.
Main text
The model
The cholera model classifies the human population at time t, denoted by N(t), into susceptible individuals S(t), cholera infected individuals I(t) and recovered individuals R(t) such that,
An additional compartment B(t), representing the concentration of vibrios in contaminated water has also been incorporated in the model. We use the model in [8] and incorporate aspects of limited resources as proposed in [20]. The model involves assumptions which are of critical importance and these are:

i.
The recovery rate depends on both the number of infectious individuals (I) and the hospital bedpopulation ratio (b).

ii.
All the infected individuals cannot recover unless they get treated in hospitals.

iii.
Recovered individuals are not permanently immune to the disease.
Susceptible individuals are recruited into the community either by birth or immigration at a rate \(\mu N\). Susceptible individuals can be infected either through humantohuman transmission or by ingesting environmental vibrios from contaminated aquatic reservoirs at the rates \(\beta _{1}I\) and \(\beta _{2}\dfrac{B}{B + k}\), respectively. The recovery rate of infected individuals is given by \(\gamma \). This recovery rate includes the hospital bedpopulation ratio, \(b>0\) and also depends on infected individuals I, such that, it is a function of both b and I. The recovery rate \(\gamma (b, I)\) is thus given as follows:
where \(\gamma _{1}\) is the maximum per capita recovery rate due to the sufficient health care resource and few infectious individuals as well as the inherent property of a specific disease, \(\gamma _{0}\) is the minimum per capita recovery rate due to the function of basic clinical resources. This recovery function was firstly used in [21]. The following assumptions for the recovery rate \(\gamma (b, I)\) are made:
 (H1):

\(\gamma (b, I) > 0\) for \(I \ge 0\), \(b > 0\), and \(\gamma (b, 0) = \gamma _{1} > 0\),
 (H2):

\(\frac{\partial \gamma (b, I)}{\partial I} < 0\), \(\displaystyle \lim _{I\rightarrow \infty } \gamma (b, I) = \gamma _{0} > 0\) and \(\displaystyle \lim _{I\rightarrow \ 0} \gamma (b, I) = \gamma (b,0) = \gamma _{1}\),
 (H3):

\(\frac{\partial \gamma (b, I)}{\partial b} > 0\), \(\displaystyle \lim _{b\rightarrow \infty } \gamma (b, I ) = \gamma _{1}\) and \(\displaystyle \lim _{b\rightarrow \ 0} \gamma (b, I ) = \gamma _{0}\).
We assume a constant size population with natural death rate given by \(\mu \). Cholerainfected individuals contribute to V. cholerae in the aquatic environment at rate \(\alpha \) and vibrios have a net death rate \(\delta \) in the environment. The differential equations for the Cholera model are;
We assume that all parameters are positive and the initial conditions of system (2) are given by: \(S(0) = S_{0} > 0,\) \(I(0) = I_{0} \ge 0,\) \(R(0) = R_{0} \ge 0,\) \(B(0) = B_{0} > 0.\)
Analysis of the model
Nondimensionalization of the model
Our system of equations has different dimensions with respect to the human population and V. cholerae. To make system (2) dimensionless, the following substitutions are made: \(S = sN\), \(I = iN\), \(R = rN\), \(B=xN\), \(k={\hat{k}}N\) and \(b={\hat{b}}N\) with \(s + i + r = 1\). The new system becomes:
Positivity of solutions
Since system (3) describes changes in the population of humans, it is considered mathematically and epidemiologically wellposed if it satisfies the positivity and boundedness conditions.
Lemma 1
Given that the initial conditions of system (3) are positive, the solutions s(t), i(t), r(t) and x(t) are nonnegative for all \(t> 0\).
Proof
Assume that \(t_{1} = \sup \{ t> 0: s> 0, i> 0, x > 0 \} \in (0, t]\). Thus, \(t_{1} > 0\). Let \(\lambda (t) = \beta _{1} i + \dfrac{\beta _{2} x}{x + {\hat{k}}}\), it follows from the first equation of (3) that,
so that,
Similarly, it can be shown that \(i(t) > 0\), \(r(t)>0\) and \(x(t) > 0\), for all time \(t > 0\). \(\square \)
Invariant region
Theorem 1
Let \(\left( s(t), i(t), r(t), x(t) \right) \) be the solution of system (3) with initial conditions \( ( s_{0}, i_{0}, r_{0}, x_{0} ) \). The compact set,
is positively invariant and attracts all solutions in \({\mathbb {R}}^{4}_{+} \).
Proof
We follow the proof given in [22]. Consider, \(W(t) = (W_{H}, W_{B}) = (s + i+r, x)\). The time derivative of W(t) is given by
This gives
From (4), we have \(\dfrac{dW}{dt} \le 0\) which implies that \(\Phi \) is a positively invariant set. We also note that by solving (4) we have;
where \(W_{H} (0)\) and \(W_{B} (0)\) are the initial conditions of \(W_{H} (t)\) and \(W_{B}(t)\) respectively. Thus, \( 0 \le \left( W_{H} (t), W_{B}(t) \right) \le \left( 1, \dfrac{\alpha }{\delta } \right) \) as \(t\rightarrow \infty \) and hence \(\Phi \) is an attractive set. \(\square \)
Disease free steady state and the basic reproduction number
System (3) has a disease free steady state given by
a scenario depicting an infectionfree state in the community or society. The basic reproduction number, \({\mathcal {R}}_{0}\), defined as the expected number of secondary cases produced by a single infectious individual in a completely susceptible population over the duration of its infectious period, is a threshold parameter that allows us to predict whether the disease will die out or persist [23]. Generally, \({\mathcal {R}}_{0} < 1\) means that the disease cannot invade the population and \({\mathcal {R}}_{0} > 1\) means that each infected individual produces more than one secondary infected individual.
Denote the basic reproduction number of system (3) by
Here, \({\mathcal {R}}_0\) is the sum of two subreproduction numbers representing the contributions of individuals in compartments i and x respectively.
Local stability of the diseasefree steady state
We now show that the diseasefree equilibrium point \({\mathcal {E}}_0\) is locally asymptotically stable whenever \({\mathcal {R}}_0<1\).
Theorem 2
The diseasefree equilibrium point \({\mathcal {E}}_0\) of system (3) is locally asymptotically stable if \({\mathcal {R}}_0< 1\) and is unstable if \({\mathcal {R}}_0>1\).
Proof
The Jacobian matrix of system (3) at \({\mathcal {E}}_0\) is given by
We determine the local stability of the diseasefree equilibrium by the following submatrix of \(J({\mathcal {E}}_0)\),
We note that all offdiagonal elements are positive, thus we now consider matrix \(J_1({\mathcal {E}}_0)\). We claim that \(J_1({\mathcal {E}}_0)\) is an Mmatrix. Multiplying matrix \(J_1({\mathcal {E}}_0)\) by the positive \(3\times 1\) matrix \(X_1=\left[ \mu \delta ,~\delta \gamma _1,~\alpha \mu \right] ^{T}\), we have
where \(X_2\) is a positive \(3\times 1\) matrix given by \(X_2=\left[ \mu \delta (\mu +\gamma _1),~0,~0\right] ^{T}\). Since \(J_1({\mathcal {E}}_0)\) is an Mmatrix, it follows that all eigenvalues of \(J_1({\mathcal {E}}_0)\) have negative real parts, which implies the local asymptotic stability of the diseasefree equilibrium if \({\mathcal {R}}_0<1\). Also, we show that the determinant of \(J_1({\mathcal {E}}_0)\) is given by
Thus, the matrix \(J_1({\mathcal {E}}_0)\) has eigenvalues with negative real parts if \({\mathcal {R}}_0<1\), which implies the stability of the diseasefree equilibrium. This completes the proof. \(\square \)
Endemic steady state
The endemic equilibrium of system (3) always satisfies
From the third and last equation of (6), we have that
Substituting the expression for \(x^*\) in (7) into the first equation of (6) we obtain
Substituting (7) and (8) into the second equation of (6) leads to the following fourth order polynomial equation
Solving (9) gives \(i^*=0\) which corresponds to the diseasefree equilibrium or
where
We can clearly note that, \(\nu _0>0\Leftrightarrow {\mathcal {R}}_0<1\) and \(\nu _0<0\Leftrightarrow {\mathcal {R}}_0>1\). The number of possible positive real roots of polynomial (10) are determined using the Descartes Rule of Signs. The various possibilities for the roots are shown in the presentation below.
\(\nu _{3} > 0\)  

\(\nu _{2} > 0\)  \(\nu _{2} < 0\)  
\(\nu _{1} > 0\)  \(\nu _{1} < 0\)  \(\nu _{1} > 0\)  \(\nu _{1} < 0\)  
\(\nu _{0} > 0\)  \(\nu _{0} < 0\)  \(\nu _{0} > 0\)  \(\nu _{0} < 0\)  \(\nu _{0} > 0\)  \(\nu _{0} < 0\)  \(\nu _{0} > 0\)  \(\nu _{0} < 0\)  
\((R_{0} < 1)\)  \((R_{0} > 1)\)  \((R_{0} < 1)\)  \((R_{0} > 1)\)  \((R_{0} < 1)\)  \((R_{0} > 1)\)  \((R_{0} < 1)\)  \((R_{0} > 1)\)  
\(i^{*}\)  0  1  2  1  2  3  3  1 
Existence of backward bifurcation
We establish conditions for the existence of backward bifurcation following Theorem 4.1 proven in [24]. We shall make the following change of variables: \(s=x_{1},~i=x_2,~r=x_3,~x=x_4\), so that \(\text{ N }=\displaystyle \sum\nolimits _{n=1}^{4}{x_n}\). We now use the vector notation \(X=(x_{1},x_{2},x_{3},x_{4})^{T}\). Then, system (3) can be written in the form \(\dfrac{dX}{dt}=F(t,x(t))=(f_{1},f_{2},f_{3},f_{4})^T\), where
Here \(\beta _2=\eta \beta _1\), with the following possibilities on the value of \(\eta \); \(\eta =1\), \(\eta \in (0,1)\) or \(\eta >1\). Let \(\beta _1\) be the bifurcation parameter, \({\mathcal {R}}_0=1\) corresponds to
The Jacobian matrix of system (3) at \({\mathcal {E}}_0\) when \(\beta _1=\beta ^*_1\) is given by
System (11), with \(\beta _1=\beta ^*_1\) has a simple eigenvalue, hence the center manifold theory can be used to analyze the dynamics of system (3) near \(\beta _1=\beta ^*_1\). It can be shown that \(J^*({\mathcal {E}}_0)\), has a right eigenvector given by \(w=(w_1,w_2,w_3,w_4)^{T}\), where
Further, the left eigenvector of \(J^*({\mathcal {E}}_0)\), associated with the zero eigenvalue at \(\beta _1=\beta ^*_1\) is given by \(v=(v_1,v_2,v_3,v_4)^{T}\), where
We compute a and b in order to apply Theorem 4.1 in [24]. For system (11), the associated nonzero partial derivatives of F at the diseasefree equilibrium are as follows:
It thus follows that
with
Note that if \({\hat{b}} <{\hat{b}}^*\) then \(\text{ a }>0\) and \(\text{ a }<0\) if \({\hat{b}}>{\hat{b}}^*\). Lastly,
We thus have the following result
Theorem 3
If \({\hat{b}} <{\hat{b}}^*\), then system (3) has a backward bifurcation at \({\mathcal {R}}_0=1\).
We observe from the results above that a backward bifurcation occurs at \({\mathcal {R}}_0 = 1\) if and only if \({\hat{b}}<{\hat{b}}^*\) is satisfied. From this, we can deduce that when the hospital bedpopulation ratio, \({\hat{b}}\) is below the critical threshold \({\hat{b}}^*\), then the number of hospital beds available to the population are below capacity and thereby lead to some individuals failing to access treatment. In such a case the prevalence of cholera infection remains high leading to a backward bifurcation, see Fig. 1. The existence of a backward bifurcation is also illustrated through numerical example by creating bifurcation diagram around \({\mathcal {R}}_0 =1\) (Fig. 1). To draw a bifurcation curve (the graph of \(i^*\) as a function of \({\mathcal {R}}_0\)), we fix \(\mu = 0.03;~\beta _1 = 0.2;~\beta _2 = 0.1;~{\hat{k}} = 0.9;~{\hat{b}} = 0.1;~\gamma _0 = 0.006;~\gamma _1 = 0.13;~\alpha = 0.282;~\delta = 0.5\). For this case we have that \({\hat{b}}^*=0.1676 >{\hat{b}}\). The solid lines denote stable states and the dotted lines denote unstable states.
Remark
When the model exhibits backward bifurcation, reducing \({\mathcal {R}}_0\) below unit is not sufficient to control the cholera epidemic.
Results and discussion
Numerical simulations
We perform some numerical simulations of system (3) to support our theoretical findings.
Estimation of parameters
Parameter values used for numerical simulations are given in Table 1.
Numerical results
Using the parameter values from Table 1, we obtain \(R_{ 0} < 1\) for the initial conditions \(s(0) = 0.80\), \( i(0) = 0.15\), \(r(0) = 0.05\), \( x(0) = 0.40\).
The limited resource parameter b, is varied in Fig. 2. It is shown that as b is increased, the infection population decrease. This means that increasing b, that is, when hospitalbeds are increased during an outbreak, there is a high chance that the disease will not persist.
Conclusions
A deterministic compartmental model with a nonlinear recovery rate was formulated to study and analyze the impact of available resources of the health care system on the transmission dynamics of Cholera. The recovery rate accounts for the number of available hospital beds per 10,000 population represented by the parameter \({\hat{b}}\) which is the critical index reflecting the resources of the health care system available to the public. Compared with previous cholera models, the work contained in this study is the first attempt to model the impact of limited resources of the health care system on the spread of cholera, with particular emphasis on the hospital beds.
It was shown that the disease free steady state is locally asymptotically stable whenever \({\mathcal {R}}_{0}<1\) and unstable otherwise. Inclusion of a non linear recovery rate has resulted in the existence of multiple endemic equilibria and the model exhibiting the phenomenon of backward bifurcation. The classical \({\mathcal {R}}_0\)threshold is not the key to control disease spread within a population. This was shown to result, in particular when the parameter \({\hat{b}}\) is low enough below \({\hat{b}}^*\). However, the cases of cholera infection decrease if there are a sufficient number of hospital beds, that is, when \({\hat{b}}>{\hat{b}}^*\). Therefore, in order to eradicate the disease in a community, effort must be targeted to increasing hospital resources.
Limitations
Like in any model development, the model is not without limitations. The model can be extended by inclusion of other control measures such as vaccination and disinfection.
Availability of data and materials
Estimation of parameters have been stated throughout the body of the paper and included in the reference section. The graphs were produced using the MATLAB software that is available from https://www.mathworks.com/products/matlab.html.
Abbreviations
 WHO:

World Health Organization
 MoHCC:

Ministry of Health and Child Care
 HBPR:

hospital bedpopulation ratio
References
 1.
Brachman PS, Abrutyn E. Bacterial infections of humans: epidemiology and control. New York: Springer; 2009.
 2.
CDC. Cholera–vibrio cholerae infection. Antlanta: Centers for Disease Control and Prevention; 2004.
 3.
Chaignat CL. 10 facts on cholera. Geneva: World Health Organization; 2014.
 4.
Jensen MA, Faruque SM, Mekalanos JJ, Levin BR. Modeling the role of bacteriophage in the control of cholera outbreaks. Proc Natl Acad Sci USA. 2006;103(12):4652–7.
 5.
World Health Organization. CholeraZimbabwe. Disease outbreak news. 2018. http://www.who.int/csr/don/20september2018cholerazimbabwe/en/.
 6.
The Herald Newspaper Zimbabwe. Cholera death toll hits 30. 2018. https://www.herald.co.zw/choleradeathtollhits30/.
 7.
de Magny GC, Murtugudde R, Sapiano MR, Nizam A, Brown CW, Busalacchi AJ, Yunus M, Nair GB, Gil AI, Lanata CF. Environmental signatures associated with cholera epidemics. Proc Natl Acad Sci. 2008;105(46):17676–81.
 8.
HoveMusekwa SD, Nyabadza F, Chiyaka C, Das P, Tripathi A, Mukandavire Z. Modelling and analysis of the effects of malnutrition in the spread of cholera. Math Comput Model. 2011;53(9):1583–95.
 9.
AlArydah M, Mwasa A, Tchuenche JM, Smith RJ. Modeling cholera disease with education and chlorination. J Biol Syst. 2013;21(04):1340007.
 10.
Bayleyegn YN. Mathematical analysis of a model of cholera transmission dynamics. Ph.D. thesis. Canada: University of Manitoba; 2009.
 11.
Codeço CT. Endemic and epidemic dynamics of cholera: the role of the aquatic reservoir. BMC Infect Dis. 2001;1(1):1.
 12.
Manfredi P, D’Onofrio A. Modeling the interplay between human behavior and the spread of infectious diseases. Berlin: Springer; 2013.
 13.
Mukandavire Z, Liao S, Wang J, Gaff H, Smith DL, Morris JG. Estimating the reproductive numbers for the 2008–2009 cholera outbreaks in Zimbabwe. Proc Natl Acad Sci. 2011;108(21):8767–72.
 14.
Mwasa A, Tchuenche JM. Mathematical analysis of a cholera model with public health intervention. Biosystems. 2011;105(3):190–200.
 15.
Neilan RLM, Schaefer E, Gaff H, Fister KR, Lenhart S. Modeling optimal intervention strategies for cholera. Bull Math Biol. 2010;72(8):2004–18.
 16.
Nelson EJ, Harris JB, Morris JG, Calderwood SB, Camilli A. Cholera transmission: the host, pathogen and bacteriophage dynamic. Nat Rev Microbiol. 2009;7(10):693–702.
 17.
Wang X, Gao D, Wang J. Influence of human behavior on cholera dynamics. Math Biosci. 2015;267:41–52.
 18.
Boaden R, Proudlove N, Wilson M. An exploratory study of bed management. J Manag Med. 1999;13(4):234–50.
 19.
World Health Organization. World health statistics. World Health Organization (2005–2011).
 20.
Abdelrazec A, Bélair J, Shan C, Zhu H. Modeling the spread and control of dengue with limited public health resources. Math Biosci. 2016;271:136–45.
 21.
Shan C, Zhu H. Bifurcations and complex dynamics of an sir model with the impact of the number of hospital beds. J Diff Equ. 2014;257(5):1662–88.
 22.
Mushanyu J, Nyabadza F, Muchatibaya G, Mafuta P, Nhawu G. Assessing the potential impact of limited public health resources on the spread and control of typhoid. J Math Biol. 2018;77:647. https://doi.org/10.1007/s0028501812199.
 23.
van den Driessche P, Watmough J. Reproduction numbers and subthreshold endemic equilibria for compartmental models of disease transmission. Math Biosci. 2002;180:29–48.
 24.
CastilloChavez C, Song B. Dynamical models of tuberculosis and their applications. Math Biosci Eng. 2004;1(2):361–404.
Acknowledgements
Farai Nyabadza acknowledges with gratitude the support from National Research Foundation and University of Johannesburg for the production of this manuscript. Jennifer Mawunyo Aduamah would like to thank the support from the African Institute for Mathematical Sciences (AIMS). Josiah Mushanyu acknowledge, with thanks, the support of the Department of Mathematics, University of Zimbabwe.
Funding
Not applicable.
Author information
Affiliations
Contributions
FN conceived of the study, and participated in model formulation and helped to draft the manuscript. JMA participated in formulating the model, and carried out the stability analysis of the model steady states and performed the numerical analysis. JM participated in the stability analysis of the model steady states, numerical analysis and drafting of the manuscript. All authors read and approved the final manuscript.
Corresponding author
Correspondence to Josiah Mushanyu.
Ethics declarations
Ethics approval and consent to participate
No ethical approval was required for this project as this is secondary research.
Consent to publish
Not applicable.
Competing interests
The authors declare 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. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Received
Accepted
Published
DOI
Keywords
 Cholera
 Nonlinear recovery rate
 Hospital bed
 Backward bifurcation
 Basic reproduction number