Modeling the health impact of water and sanitation service deficits on waterborne disease transmission

Cholera is a waterborne disease that continues to pose serious public health problems in many developing countries. Increasing water and sanitation coverage is a goal for local authorities in these countries, as it can eliminate one of the root causes of cholera transmission. The SIWDR (susceptible–infected–water–dumpsite– recovered) model is proposed here to evaluate the effects of the improved coverage of water and sanitation services in a community at risk of a cholera outbreak. This paper provides a mathematical study of the dynamics of the water and sanitation (WatSan) deficits and their public health impact in a community. The theoretical analysis of the SIWDR model gave a certain threshold value (known as the basic reproductive number and denotedR0) to stop the transmission of cholera. It was found that the disease-free equilibrium was globally asymptotically stable whenever R0 ≤ 1. The unique endemic equilibrium was globally asymptotically stable wheneverR0 > 1. Sensitivity analysis was performed to determine the relative importance of model parameters to disease transmission and prevention. The numerical simulation results, using realistic parameter values in describing cholera transmission in Haiti, showed that improving the drinking water supply, wastewater and sewage treatment, and solid waste disposal services would be effective strategies for controlling the transmission pathways of this waterborne disease.


Introduction
The most common waterborne diseases that pose a major public health risk in communities and continue to be leading causes for illnesses and deaths worldwide, especially in developing countries, include diarrhea, cholera, giardia, dysentery, Escherichia coli, salmonella, shigella, typhoid, and hepatitis A [1]. Waterborne diseases are caused by a lack of safe drinking water and proper sanitation, which is a major cause [2]. According to the World Health Organization (WHO), cholera is an infectious disease characterized by severe watery diarrhea and is caused by eating food or drinking water contaminated with the bacterium Vibrio cholerae. Diarrhea kills an estimated 829,000 people globally each year as a result of contaminated drinking water, sanitation, and hand hygiene [3]. In addition, morbidity and mortality are worse in developing countries where communities lack the resources to provide sustained access to adequate water and sanitation services [4]. For instance, in 2010, Haiti suffered from a cholera outbreak with 665,000 cases and 8,183 deaths [5]. A cholera outbreak also occurred in Yemen in 2016 with more than 1.2 million cases and 2,500 deaths, of which 30% were children under the age of five [6]. In many low income countries, inadequate investment in water and sanitation infrastructure has perpetuated the problem. The inadequate investment resulted in illnesses among citizens. Loss of work productivity due to illnesses lowered individual income; hence, there was a lower national income. Consequently, low income countries had the least amount of money to invest in water and sanitation infrastructure that is necessary to combat the burden of waterborne diseases [7,8]. Furthermore, open dumpsites in developing countries play an important role in the transmission of diseases [9,10]. Open dumpsites can be an important and large reservoir of pathogens in these communities. Pathogens in open dumpsites can get into water sources through leachate, which contaminates unprotected groundwater sources, or through precipitation and runoff into surface water [11][12][13]. For example, open dumpsites in Yemen are located close to communities. This situation increases the health risks in these communities, especially from clinical waste, which has been demonstrated to contribute to cholera outbreaks in Yemen [14].
Since the connections between unclean water, sanitation, and health are well established [7,[15][16][17], combating waterborne diseases requires multiple strategies. These strategies include increasing access to safe drinking water and improved sanitation and cutting disease transmission pathways by having good hygiene practices [18]. According to the WHO, 94% of diarrhea cases are preventable by healthy environments, which include increasing access to clean water and sanitation services, and improved hygiene [19]. For sanitation services [7,20], systematic reviews concluded that improved sanitation can reduce diarrhea cases by 32%-37%. This indicates that providing adequate safe drinking water and sanitation services reduces the incidence of acute infectious diarrhea [21][22][23]. Moreover, a study in Brazil and El Salvador found that improving sewerage coverage led to a significant reduction in diarrhea prevalence [24][25][26][27]. In this paper, therefore, the impact of improving water and sanitation services and the management of open dumpsites are explored as mechanisms to control waterborne disease transmission.
Mathematical models have been used to understand the transmission dynamics of waterborne disease, to improve control and prevention strategies, and to inform public health impacts [28][29][30][31][32]. Tien and Earn [28] presented the susceptible-infected-recovered (SIR) model with pathogen concentration, which was called the susceptible-infected-waterrecovered (SIWR) model when applied to waterborne disease transmission. Their study showed that applying the SIR model rather than using the SIWR model underestimates the basic reproductive number and overestimates the infectious period when the pathogen decay rate is slow in the water compartment. Eisenberg et al. [29] showed that parameters in the SIWR model were globally identifiable. They estimated this by using the data from a cholera outbreak in Angola as a case study. Also, their study concluded that including environmental data was important when analyzing waterborne diseases. The SIWR model was further extended by having multiple patches within a common water resource. Their model assumed that the human-to-human transmission pathway occurs within a patch and the environment-to-human pathway happens across multiple patches with a shared water resource. Their study also explored the effects of heterogeneity on the final size of an outbreak and the efficiency of interventions [30]. This type of spatially explicit model [31][32][33] was also used to investigate waterborne disease epidemics that account for local communities of susceptible and infected individuals in a spatially explicit arrangement of nodes linked by networks having different topologies. A network epidemic model of waterborne disease was studied by Wang and Cao [34] and Collins and Govinder [35]. These models considered the heterogeneity of different water sources that were contaminated. The optimal control of waterborne disease having multiple patches was also studied. Tuite et al. [36] proposed a gravity model to predict the sequence and timing of regional cholera epidemics in Haiti (2010) and suggested that adaptive strategies for vaccination may provide a modest reduction in morbidity and mortality in the economically challenged country. Collins et al. [37] studied the effect of heterogeneity in transmission due to socioeconomic status on the dynamics of waterborne diseases. Even though all the aforementioned studies have contributed immensely to improving the understanding of waterborne disease dynamics, other pathogen reservoirs that cause waterborne disease outbreaks have not been studied. In this paper, we address this research gap by formulating an appropriate mathematical epidemiological model that incorporates open dumpsites.
Another important issue of interest is the study of globally asymptotic properties of epidemiological models. This can help design effective control strategies intended to permanently reduce pathogen spread, or even break the chain of disease transmission. Understanding transient phenomena in potential epidemiological scenarios can also be important. It can help to prevent transitory epidemics that are triggered by external perturbations to a system, in which endemic transmission is not possible. That is, the disease-free equilibrium (DFE) is globally asymptotically stable [38]. This results in the threshold value of epidemiological models, which is called the basic reproductive number and denoted by R 0 . The study conducted in [28] showed that estimates of R 0 and the infectious period can be sensitive to the relative contributions of the different transmission pathways and the pathogen's lifetime in water. Further studies [2,28,37,39,40] and references (therein) confirmed that the sensitivity of R 0 is relevant for dynamical studies of waterborne diseases because the basic life history and epidemiological traits remain uncertain for many waterborne pathogens.
Based on the above discussion, this paper proposes a compartmental mathematical model to study the health impact of water and sanitation service deficits on the spread of waterborne diseases. The paper is organized as follows. In Sect. 2, the model proposed by Tien and Earn [28] is modified by incorporating water and sanitation services together with open dumpsites. In Sect. 3, the proposed model is analyzed to gain insights into the mechanisms of waterborne disease transmission and the control measures for waterborne diseases. The existence and the number of equilibria are established. The threshold value of the proposed model is derived by using the next generation method. The global stabilities of equilibria are analyzed by constructing the suitable Lyapunov functions and using LaSalle's invariance principle. In Sect. 4, numerical simulations are illustrated to determine the appropriate parameter values used for the proposed model and to determine the parameters that are sensitive to disease prevalence and the impact of water and sanitation (WatSan) services. Finally, conclusions are given in Sect. 5.

System modeling
Individuals can be infected with waterborne disease if they come near to an open dumpsite and contact waste directly [13]. This is because pathogens in open dumpsites can get into water sources through leachate and contaminate unprotected groundwater sources, or through precipitation that washes into surface water [11,12]. Improving WatSan services is an effective means of reducing the number of infected individuals and the number of people who die each year as a result of diarrhea [41]. An increase in WatSan services would reduce the incidence of WatSan-related diseases, such as enteric diarrheal disease (EDD). Therefore, the system modeling included open dumpsites as reservoirs of pathogens and WatSan services, namely, the drinking water supply (DWS), wastewater and sewage treatment (WST), and municipal solid waste management (MSW). To investigate the impact caused by WatSan deficits on individual health and the state of the environment, we modified the model presented by Tien and Earn [28] with the following assumptions. We assumed that disease transmission occurred through two pathways, namely, a direct transmission pathway (human-to-human) and an indirect transmission pathway (environment-to-human). The indirect transmission pathway was formulated based on the SIR epidemic model for monitoring the dynamics of the sub-populations of susceptible individuals (S(t)), infectious individuals (I(t)), and recovered individuals (R(t)), so that the total population at time t is N(t) = S(t) + I(t) + R(t). The indirect transmission pathway consisted of the SIR epidemic model together with compartments W and D. W and D, respectively, represent pathogen concentration in water reservoirs and open dumpsites that are potential reservoirs of pathogens. Compartment D represents open dumpsites and other pathogen reservoirs that are not water reservoirs. Furthermore, to improve the WatSan services, let 0 ≤ r 1 ≤ 1, 0 ≤ r 2 ≤ 1, and 0 ≤ r 3 ≤ 1 represent the ratio of people who receive DWS, WST, and MSW services, respectively, in a community. A flow diagram of the formulated model, the susceptible-infected-water-dumpsite-recovered (SIWDR) model, is depicted in Fig. 1. The dynamics of each sub-population in Fig. 1 are described next.
The number of susceptible individuals increases at the recruitment rate of μN , where μ is the birth rate (intrinsic growth rate). The birth rate is assumed to be equal to the natural death rate so that the population size N remains constant. Thus, the susceptible    The pathogen shedding rates from infected individuals into water and dumpsite compartments are given by (1r 2 )α 1 I and (1r 2 )α 2 I, respectively. The pathogen shedding rate from the dumpsite compartment into the water compartment is given by (1r 3 )ψD.
The parameters ξ W and ξ D give the decay rate of pathogens in the water and dumpsite compartments, respectively.
Therefore, based on the above descriptions and Fig. 1, the SIWDR model for the pathways of infection in waterborne diseases is given by the following deterministic system of nonlinear differential equations: where the descriptions of the model parameters and sample units in system (2.1) are given in Tables 1-2. Since the above model monitors human and pathogen populations, all associated parameters and state variables have nonnegative values.

Basic properties
Adding the first, second, and fifth equations of system (2.1) gives dN/dt = 0. Thus, the total human population N is constant. For convenience of analysis, we rewrite system (2.1) in terms of new dimensionless quantities by setting the state variables: so that The model parameters are replaced by Using the expression above, system (2.1) is reduced to where r = 1si. This allows us to solve system (2.1) by studying system (3.2). Clearly, lim t→∞ sup(s + i) ≤ 1 implies lim t→∞ sup(i(t)) ≤ 1. It follows from the fourth equation of (3.2) that Thus, from the third equation of system (3.2), we have It is easy to show that .
Hence, the feasible region Thus, the existence, uniqueness, and continuation results for model (3.2) hold in . All solutions that start in remain in for all t ≥ 0 [42]. Hence, model (3.2) is mathematically and epidemiologically well posed, and it is sufficient to consider the dynamics of the flow generated by model (3.2) in .

Local stability
In the absence of infection and with adequate water and sanitation services (i = 0, w = 0 and d = 0), system (3.2) has a unique disease-free equilibrium (DFE) denoted by According to the next generation matrix and using the notation in [43], system (3.2) is rewritten as with Q 1 =r 1 β W sw +r 3 β D sd + β I si + μs and Q 2 =r 3 α 2 ξ W ψ, respectively. Clearly, x 0 = (0, 0, 0, 1) is a disease-free equilibrium of system (3.4), which is identical to ε 0 of system (3.2). It is easy to show that system (3.4) satisfies assumptions (A1)-(A5) in [43]. Further, the derivatives DF (x 0 ) and DV(x 0 ) are partitioned as where the matrices F (for the new infection terms) and V (for the transition terms) are given, respectively, by According to Theorem 1 in [43], the following result is established.
The local stability that resulted in Lemma 3.1 implies that for R 0 < 1 the disease dies out and for R 0 > 1 the infected population is not reduced to zero. Thus, the quality that R 0 > 1 is the threshold value of model (3.2). Furthermore, the result in Lemma 3.1 verifies that if the initial sizes of the sub-groups of the model are within the basin of attraction of ε 0 , the total number of infected people in the population can be reduced to zero whenever R 0 < 1. Moreover, it was found that the value of R 0 depends on the ratios of people who receive the DWS, WST, and MSW services: r 1 , r 2 , and r 3 , respectively. This theoretical determination of conditions of these ratios can make R 0 independent of the initial sizes of the five state variables. This is of great public health interest, which is later illustrated by analyzing the global stability.

Global stability of DFE
We claim the following result.

Proof We introduce a new function
. The derivative of L with the solutions of (3.2) is The Lyapunov-LaSalle theorem [44] implies that all paths in approach the largest positively invariant subset Thus, we easily obtain lim t→∞ w(t) = 0. Hence, all solution paths in approach the DFE ε 0 . This completes the proof of this theorem.

Local stability
Linearizing model (3.2) around its unique EE is analyzed by the use of the center manifold theory [45] as described in Theorem 4.1 [46]. The following theorem, thus, is established.
Proof The Jacobian of system (3.2) evaluated at ε 0 is given by Suppose, that β I = β * I is a bifurcation parameter. Solving R 0 = 1 for β I gives It was found that system (3.2) with β I = β * I has at least one non-hyperbolic equilibrium point. Thus, the centre manifold theory [45] as described in Theorem 4.1 [46] is used to analyze the LAS of EE ε 1 of the system (3.2) with β I = β * I . The right eigenvector of J(β * I ) with β I = β * I is given by w = [w 1 , w 2 , w 3 , w 4 , w 5 ] T , where , and v 2 > 0.
For system (3.2), the associated non-zero partial derivatives of the function f i (s, i, w, d, r), i = 1, . . . , 5 on the right-hand side of each equation are given by for variables x 1 = s, x 2 = i, x 3 = w, and x 4 = d.
Using the above expressions, parameters a and b as given in Theorem 4.1 [46] are computed.
This theorem indicates that if R 0 > 1, the EE ε 1 is locally asymptotically stable (that is, the disease will persist in the community if the initial sizes of the four-state variables are within the vicinity of ε 1 ). If the EE ε 1 is globally asymptotically stable, then any trajectories tend toward the attractor of ε 1 , of system (3.2), regardless of the initial conditions. This means that the attracting basin of trajectories in the dynamical system is the state space [47]. In a biological sense, this property gives the theoretical determination of conditions that can cause R 0 > 1, irrespective of the initial sizes of the five state variables that increase the risk of a disease outbreak. The Lyapunov function is constructed and LaSalle's invariance principle, which is the criterion for the asymptotic stability of an autonomous dynamical system [48][49][50][51], is applied to prove the global stability of ε 1 . The following theorem, therefore, is established.
Proof Consider the Lyapunov function where The derivative of V (t) and the solutions of system (3.2) are given by At steady state when R 0 > 1, and from the first four equations of system (3.2), we have μ =r 1 β W s * w * +r 3 β D s * d * + β I s * i * + μs * , Substituting expressions (3.9) into (3.8) and rearranging give It follows from the arithmetic mean is greater than or equal to the geometric mean that the following inequalities hold: Using the inequalities in (3.10) and that all the model parameters are nonnegative, it is clear that dV /dt ≤ 0 whenever R 0 > 1. Furthermore, dV /dt = 0 if and only if s = s * , i = i * , d = d * , and w = w * . Hence, V (t) is a Lyapunov function in . By LaSalle's invariance principle [52], it is found that s → s * , i → i * , d → d * , and w → w * as t → ∞. Hence, the EE (ε 1 ) of model (3.2) is GAS. That is, every solution of model (3.2) with initial condition in approaches ε 1 as t → ∞ whenever R 0 > 1.

Estimation values of model parameters
From the outbreak size of cholera in Haiti in 2010, the death toll was approximately 9,700, and there were over 800,000 suspected cases [53]. Model (3.2) used this cholera outbreak as a case study. The model was simulated with N = 2,200,000, which is the approximate population in the capital of Haiti in 2010 [54]. The parameter values are given in Table 3 for  2) and the cumulative real data from Haiti, 2010 [65]. Red points denote real data and the blue line denotes predicted data predicting the number of infected individuals. The results are then used to compute the cumulative number of predicted infected individuals by solving the differential equation where k = 0.99 is the rate of progression in the monthly infectious state. With all values given in Table 3, the basic reproductive number is R 0 = 1.2406, which is within the range (1.06, 2.78) [36,66]. This value indicates that cholera persists in the community and can result in a disease outbreak, as guaranteed by Theorem 3.4. As shown in Fig. 2, the cumulative predicted data are compared with the cumulative real cholera data in Haiti 2010, which are reported by the Pan American Health Organization (PAHO) [65]. Figure 2 shows that the cumulative predicted data produced by model (3.2) are close to the real data generated by using a statistical measure, the coefficient of determination, R 2 = 0.9946. This study result verifies that model (3.2) can be used to investigate the cholera transmission in Haiti (2010). Moreover, the values of model parameters given in Table 3 are appropri-ate for model (3.2). These values are used to investigate the influential factors that can mitigate the 2010 cholera outbreak in Haiti in the next sections.

Sensitivity analysis of the basic reproductive number
Sensitivity analysis [67] is a tool for determining the relative importance of the different factors responsible for disease transmission. The initial disease transmission is directly related to the basic reproductive number R 0 . We computed the sensitivity indices of R 0 in a model with parameter p, defined by using partial derivatives Table 4 shows the order of sensitivity indices of R 0 in relation to thirteen model parameters, evaluated at the baseline values given in Table 3. The values are arranged from the most to the least sensitive parameters. A negative value means that increasing a model parameter value decreases the value of R 0 , resulting in a decrease in the infected individual density. In contrast, a positive value means that increasing a model parameter value increases the value of R 0 , resulting in an increase in the infected individual density. The model parameters (β W , α 1 , α 2 , ψ, β D , β I ) affect the increase in the initial transmission, and β W is the most sensitive parameter. In contrast, the model parameters (r 1 , γ , r 2 , r 3 , ξ W , ξ D , μ) affect the decrease in the initial disease transmission, and r 1 is the most sensitive parameter. Regarding the coverage of DWS, WST, and MSW services, the coverage of DWS service (r 1 ) is the most sensitive parameter, followed by the coverage of WST service (r 2 ) and the coverage of MSW service (r 3 ). This result indicates that improving the coverage of the DWS service is the most effective intervention to control a waterborne disease outbreak.
Although an increase in the coverage of DWS is the most effective strategy, choosing only one strategy to control a waterborne outbreak might not be an effective strategy. Next, the results of combining two of these model parameters on the value of R 0 are explored using the contour plots of R 0 (see Fig. 3). The results show that an increase in the coverage r 1 , r 2 , and r 3 services decreases the value of R 0 to be less than unity. However, it is observed that if the coverage of r 1 is 0.65, the coverage of r 2 should be increased by at least 0.38 (compared to the baseline value r 2 = 0.225) to decrease R 0 to be less than 1 (see Table 4 Sensitivity indices of R 0 to model parameters for model (3.2), evaluated at the baseline parameter values given in Table 3 Model parameters Sensitivity indices ϒ  Table 2 Fig . 3(a)). In this case, an increase in the coverage of the r 3 service has the least impact in decreasing the value of R 0 (see Figs. 3(b) and 3(c)). In summary, this study suggests that the WatSan services coverage would be a crucial factor to mitigate the health impact of the 2010 cholera outbreak in Haiti. Generally, coverages of DWS, WST, and MSW services can decrease the initial disease transmission and interrupt the transmission of waterborne diseases.

Impacts of WatSan services
The 2013-2022 plan to combat cholera in Haiti has been put in place [68]. The plan is to spend approximately $1.6 billion U.S. dollars to improve WatSan service coverage. The budget allocation for improving the coverage of DWS, WST, and MSW services was set to be 49.5%, 28.1%, and 22.4%, respectively. This budget allocation satisfies the result obtained in Sect. 4.2. Our research results show that the coverage of DWS, WST, and MSW services significantly impacts the control of waterborne diseases.
In 2010, the government proposed that every citizen must have improved water supplies and sanitation services. However, only 1.5% of water supplies and only 12.5% of sanitation facilities improved from 2010 to 2017. According to the 2017 UNICEF report [69], the proportion of the Haitian population with improved water supplies and sanitation services was only 65% and 35% of the total population, respectively. In Haiti, the progress of WatSan service improvement was not done according to the 10-year plan proposed by the Haitian government [68]. Therefore, insufficient WatSan service provisions were likely the cause of the persistent cholera outbreaks in Haiti.
In this section, the impacts of improving WatSan services are investigated as strategies to interrupt the transmission pathways of a cholera outbreak. Model (3.2) was simulated to determine the fraction of infected individuals and the fraction of pathogens in water reservoirs and dumpsites with varying percentages of coverage of the DWS service (r 1 ), WST service (r 2 ), and MSW service (r 3 ). The other parameters used in the simulation are given in Table 2. The dynamic behavior of the fraction of infected individuals and the fraction of pathogens in water reservoirs and dumpsites is shown in Figs. 4-6. From the 2010 cholera outbreak data in Haiti, we set the starting values of r 1 , r 2 , and r 3 to be 65%, 35%, and 10 %, respectively. The fraction of infected individuals (Figs. 4-6(top)) uses fixed values of two out of three of the other variables (r 1 = 65%, r 2 = 35%, or r 3 = 10%). The lines in each graph show that a change in the coverage of DWS service (r 1 ) has the greatest impact in reducing the fraction of infected individuals. In Figs. 4-6(middle), the lines in each graph show that a change in the coverage of WST service (r 2 ) has the greatest impact in reducing the fraction of pathogens in water reservoirs. In Figs. 4-6(bottom), the lines in each graph show that a change in the coverage of WST service (r 2 ) has the greatest impact in reducing the fraction of pathogens in dumpsites. Figure 4 shows that the fraction of infected individuals, the fraction of pathogens in water reservoirs, and the fraction of pathogens in dumpsites decrease with an improvement in DWS service coverage. However, with a coverage of 65%, the outbreak seems to maintain a steady state of the fraction of infected individuals (noted by the blue horizontal line in Fig. 4(top)). Relatedly, the fraction of pathogens in water reservoirs and the fraction of pathogens in dumpsites (in this case) also seem to stay constant (noted by the blue horizontal line in Figs. 4(middle) and 4(bottom)). Note that with varying r 1 , the peaks in Fig. 4(bottom) occur approximately at the same position. This is because changing the coverage of DWS service (r 1 ) does not Figure 4 Dynamic behavior of the fraction of infected individuals (top), fraction of pathogens in water reservoirs (middle), and fraction of pathogens in dumpsites (bottom), with r 2 = 35%, r 3 = 10%, and varying r 1 from 65%, 70%, 80%, 90%, and 100% Figure 5 Dynamic behavior of the fraction of infected individuals (top), fraction of pathogens in water reservoirs (middle), and fraction of pathogens in dumpsites (bottom), with r 1 = 65%, r 3 = 10%, and varying r 2 from 35%, 40%, 60%, 80%, and 100% Figure 6 Dynamic behavior of the fraction of infected individuals (top), fraction of pathogens in water reservoirs (middle), and fraction of pathogens in dumpsites (bottom), with r 1 = 65%, r 2 = 35%, and varying r 3 from 10%, 20%, 40%, 60%, 80%, and 100% directly affect the fraction of pathogens in dumpsites. Similarly, Figs. 5 and 6 show that at starting values of r 2 = 35% and r 3 = 10%, respectively, the outbreak seems to maintain a steady state of the fraction of infected individuals, fraction of pathogens in water reservoirs, and fraction of pathogens in dumpsites. Figure 7 displays the dynamic behavior of the fraction of infected individuals and the fraction of pathogens in water reservoirs and dumpsites with varying the coverage of two or three services at a time. Looking at the extreme case where (r 1 , r 2 , r 3 ) = (80, 70, 50), the fraction of infected individuals (Fig. 7(top)) and the fraction of pathogens in water reservoirs (Fig. 7(middle)) are at the lowest. However, this does not reduce the fraction of pathogens in dumpsites (Fig. 7(bottom)) to be at the lowest. The cases of (r 1 , r 2 , r 3 ) = (65, 70, 10) and (80, 70, 10) have lower peaks than those of (r 1 , r 2 , r 3 ) = (80, 70, 50). This shows that an increase in the coverage of MSW service (r 3 ) or an improvement in waste or garbage management does not reduce the fraction of pathogens in dumpsites within the first 150 months. The fraction of pathogens will eventually be reduced to zero, but over a long time. In other words, better MSW management means that garbage or the pathogens in the garbage are less likely to get into water reservoirs, and a disease is less likely to be transmitted to individuals via water. This improvement cuts the transmission pathway where pathogens in dumpsites can get into water, but not the fraction of pathogens in dumpsites.
Consequently, an increase in the coverage of WST service (r 2 ) has a greater impact (short-term) on the fraction of pathogens in dumpsites than an increase in the coverage of MSW service (r 3 ). WST has a direct impact on the fraction of pathogens in dumpsites. In many developing countries, wastewater and dumpsites are in direct contact with each other. Thus, wastewater treatment reduces the fraction of pathogens in dumpsites.
The results can be explained using Fig. 1. An increase in the coverage of DWS service (r 1 ) also affects the fraction of pathogens in dumpsites. Figure 1 shows that a change in r 1 affects the fraction of infected individuals, which affects the quantity of infected garbage/waste. Also, Fig. 1 shows that a change in the coverage of MSW service (r 3 ) has the least impact on the fraction of pathogens in dumpsites, as it does not affect compartment D directly. Better solid waste management means that pathogens in garbage are less likely to get into compartment W . This leads to fewer pathogens in water reservoirs, which yields a smaller fraction of infected individuals and less infected garbage. Figure 7(bottom) also shows that the fraction of pathogens in dumpsites is at the lowest when (r 1 , r 2 , r 3 ) = (80, 70, 10). This implies that an increase in the coverage of DWS service (r 1 ) and WST service (r 2 ) affects the fraction of pathogens in dumpsites. Overall, an increase in variable r 1 or r 2 alone has less impact on the fraction of infected individuals. With increased coverage of DWS and WST services, the 10% coverage of MSW service (r 3 ) is enough to improve the situation of a cholera outbreak in Haiti according to Fig. 7(bottom) for approximately the first 100 months. This is an important implication.
In Fig. 7(bottom), the cases of (r 1 , r 2 , r 3 ) = (65, 70, 10) and (80, 70, 10) result in the lowest peaks. This implies that a change in the coverage of WST service (r 2 ) may have the greatest impact on the fraction of pathogens in dumpsites. An increase in the coverage of DWS service (r 1 ) delineates these two lines (see 50 to 150 months). For (r 1 , r 2 , r 3 ) = (65,70,50) and (80,70,50) where r 2 = 70, the peaks are the second lowest. For (r 1 , r 2 , r 3 ) = (65,35,10) and (80,35,10), where r 2 = 35, the peaks occur approximately at the same position and at the highest. That is, it takes longer to reduce the fraction of pathogens in dumpsites Figure 7 Effects of improving DWS, WST, and MSW services coverage (r 1 , r 2 , and r 3 ) on the fraction of infected individuals and the fraction of pathogens in water reservoirs and dumpsites that are predicted by model (3.2) that have less improvement in DWS. In summary, the coverage of WST service (r 2 ) plays a significant role in reducing the fraction of pathogens in dumpsites. The same pattern happens for any of the two lines where r 2 and r 3 are the same with only a change in r 1 . Peaks at any of the two lines occur at the same position. However, the decline in the fraction of pathogens (afterward) depends on the r 1 value. Figure 7(bottom) illustrates that when r 1 and r 2 are the same, with only a change in r 1 , there are intersections between any two selected observed lines. The coverage of MSW service (r 3 ) has a long-term impact. It reduces the fraction of pathogens significantly after approximately 170 months. In contrast, the coverage of DWS service (r 1 ) changes the fraction of pathogens significantly from about 25 to 100 months (the time after where the peak occurs). Consequently, decision makers must consider short-term and long-term goals depending on the budget at that time. The timing of eliminating a disease may be of concern. Figure 8 confirms that there is a threshold value. When fixing the value of r 3 , the coverage of DWS service (r 1 ) and the coverage of WST service (r 2 ) must be maintained. At the line where (r 1 , r 2 , r 3 ) = (60,30,10), the fraction of infected individuals, the fraction of pathogens in water reservoirs, and the fraction of pathogens in dumpsites are out of control.
This study suggests that controlling and preventing waterborne diseases could be achieved by improving the service coverage of WST and MSW, to be sufficient, to protect water reservoirs and dumpsites from contamination. However, improving the service coverage of the DWS service could have the most impact among all of the WatSan services. Therefore, the WatSan service impact is explored as follows.
The results obtained in Table 5 are used as the projection of the fraction of the infected population for the next five years to cope with the 10-year plan which ends in 2022. The current data of coverage of WatSan services in 2017 from UNICEF [69] are used as the current situation (scenario 1). In scenario 1, the coverage of DWS, WST, and MSW services is 65%, 35%, and 10%, respectively. These percentages correspond to the values of r 1 = 0.65, r 2 = 0.35, and r 3 = 0.10, respectively. These values, with other parameter values given in Table 4, lead to R 0 = 1.0073. In addition, the percentage reduction in the infected population is calculated by comparing the fraction of the infected population to the baseline, computed from the parameter values in Table 4. There is a 66.0674% reduction in the fraction of the infected population for scenario 1. By fixing all parameter values in scenario 1 except the coverage of each WatSan service, the threshold service coverage for each WatSan service that makes R 0 = 1 is 65.28% of the service coverage of DWS (r 1 ), or 35.49% of the service coverage of WST (r 2 ), or 13.95% of the service coverage of MSW (r 3 ). A little improvement in each WatSan service while maintaining the other services would alleviate the disease outbreak. The experiments are extended by increasing or decreasing the percentage of r 1 , r 2 , and r 3 by 5% from scenario 1, as shown in Table 5. Scenarios 2-8 decrease r 1 , r 2 , and r 3 , while scenarios 9-15 increase r 1 , r 2 , and r 3 by 5% from scenario 1. It is found that decreasing r 1 , r 2 , or r 3 from scenario 1 while maintaining the other services could result in a reduction in the infected population that is worse than scenario 1. R 0 in scenarios 2-8 are higher than R 0 in scenario 1. In scenarios 9-15, increasing r 1 , r 2 , or r 3 from scenario 1 while maintaining the other services could result in a percentage reduction in the infected population that is better than scenario 1. In addition, R 0 in scenarios 9-15 is less than R 0 in scenario 1. The impact of decreasing r 1 , r 2 , or r 3 is greater  than the impact of increasing r 1 , r 2 , or r 3 when comparing the percentage reduction in the infected population for scenarios 2-8 and 9-15. The result in scenario 5 shows that decreasing r 1 and r 2 decreases the reduction in the infected population significantly. Scenario 12 demonstrates that increasing r 1 and r 2 increases the reduction in the infected population remarkably. However, the impact on the reduction in the fraction of the infected population is greater for scenario 5 (from 66.0674 % to 8.5914%) than for scenario 12 (from 66.0674 % to 80.0980%). Thus, it is concluded that maintaining the service coverage of r 1 and r 2 is important to prevent the transmission of waterborne diseases. Figure 8 shows the impact of improving and reducing r 1 and r 2 . Reducing both services could make water reservoirs and dumpsites more contaminated. This impacts the fraction of the infected population. Our study shows that maintaining the service coverage of MSW (r 3 ) while improving the service coverage of DWS and WST (r 1 and r 2 ) could have the most health impact on a community. Thus, our study suggests that shifting the budget away from improving the service coverage of MSW to just maintaining this service is a desirable policy. In addition, the shifted budget could be used for improving the service coverage of DWS and WST instead. According to the 2013-2022 plan to combat cholera in Haiti, the budget allocation for improving the service coverage of DWS, WST, and MSW is 49.5%, 28.1%, and 22.4%, respectively. This allocation is suitable for fighting cholera. However, to efficiently eliminate cholera from Haiti, our study recommends that the allocated budget of 22.4% to improve MSW should be reduced. This reduction in the budget should be used to improve the service coverage of DWS and WST by emphasizing the improvement of DWS.
As described earlier, proposed model (2.1) is suitable when there are other pathogen reservoirs that are important in the disease transmission process. This study extends the Tien and Earn (2010) model that takes into account only one pathogen reservoir, i.e., a water reservoir [28]. This study finds that environmental factors play an important role in controlling and preventing disease outbreaks. Considering the water compartment (W ) and dumpsite compartment (D) is more effective than considering only compartment W if dumpsites in a particular community are contaminated. Additionally, our mathematical model quantifies the qualitative data that many researchers [41] have agreed upon. Wat-San services do indeed prevent and control waterborne diseases. According to this paper, WatSan services should not be given equal priority.

Conclusions
This paper presents the modified SIR model by incorporating access to WatSan services and the prevalence of open dumpsites as pathogen enablers. The proposed model is rigorously analyzed to gain insights into its dynamic behavior. Then, it is used to investigate the health impact of WatSan service deficits on waterborne disease transmission. Using cholera as our case study, our analysis yielded some useful new insights into the dynamics and control of waterborne diseases. The basic reproductive number (R 0 ) of the model is derived and given in Equation (3.6) by using the next generation method. It was found that the disease-free equilibrium is globally asymptotically stable whenever R 0 ≤ 1 as guaranteed by Theorem 3.2. When R 0 > 1, the endemic equilibrium (EE) of the proposed model exists and is globally asymptotically stable, which is consistent with Theorem 3.4. These theoretical studies reveal that the disease dies out when R 0 ≤ 1 and the disease becomes endemic when R 0 > 1. Therefore, control strategies that set R 0 ≤ 1 would generate public health benefits in making strategy decisions after infectious diseases appear. The appropriate parameter values of the model were obtained by comparing the predicted data with real data from the 2010 cholera outbreak in Haiti. These parameter values were used to investigate the most sensitive model parameters in controlling and preventing cholera outbreaks. It was found that β W is the most sensitive parameter to the initial transmission, and r 1 is the most sensitive parameter to the initial decrease in disease transmission. The coverage of DWS services is the most sensitive parameter to the initial decrease in disease transmission in relation to the coverage of WST and MSW services. In other words, this study verifies that improving the coverage of DWS services is the most effective intervention to control a cholera outbreak. Analyzing the contour plots of surface R 0 as a function of the service coverages (r 1 , r 2 , and r 3 ) demonstrates that improving all the service coverages (DWS, WST, and MSW) would certainly decrease the initial disease transmission. Moreover, analyses of single and double control strategies showed that control strategies that involve the combined improvement of multiple water and sanitation services are better at controlling the initial disease transmission and reducing the overall infections than single strategies.
Therefore, this study suggests that to control the spread of infectious waterborne diseases like cholera, the first priority should be to improve the access to safe drinking water supplies. As funding and other resources permit, the effectiveness of control measures would subsequently be increased by adding improved treatment of contaminated water sources and improved management of open dumpsites to the improved drinking water supply. To some degree these results are intuitive. Certainly, to improve access to safe drinking water supplies, the treatment of contaminated water sources would be required. By similar reasoning, improved management of open dumpsites could reduce the extent of leachate contamination of groundwater and control harmful runoff into surface water bodies that serve as sources for drinking water supplies. Thus, prioritizing the access to safe drinking water supplies is likely to require improved treatment of contaminated water bodies and improved management of community dumpsites as effective strategies for cutting the transmission pathways of waterborne diseases like cholera. This is an important implication and a contribution of this study.
The comparative results of the effects of each service are illustrated in Table 5 for the case of enteric diarrheal disease (EDD). In this case, we found that improving DWS and WST services yields a higher reduction in EDD than improvements in the MSW service. This result suggests that WatSan service managers should prioritize investment in DWS and WST services over the MSW service in a strategy to control the incidence of EDD from waterborne sources. The budget for infrastructure services may be limited, as it is in many low-and middle-income countries. Therefore, this strategy could be useful in budget allocation for improving WatSan services to manage outbreaks of waterborne diseases. Furthermore, model (2.1) is formulated by incorporating the state variables of WatSan services and the incidence of EDD. Thus, if reliable data on the parameters of a particular disease and the unit costs of WatSan services are available, model (2.1) can be applied to design optimization algorithms for budget allocation in improving WatSan services, to minimize the size of the infected population during disease outbreaks. This is the focus of our future work.