Traveling Wave Solutions for Epidemic Cholera Model with Disease-Related Death

Zhang, Tianran; Gou, Qingming

doi:https://doi.org/10.1155/2014/409730

Abstract

Based on Codeço’s cholera model (2001), an epidemic cholera model that incorporates the pathogen diffusion and disease-related death is proposed. The formula for minimal wave speed is given. To prove the existence of traveling wave solutions, an invariant cone is constructed by upper and lower solutions and Schauder’s fixed point theorem is applied. The nonexistence of traveling wave solutions is proved by two-sided Laplace transform. However, to apply two-sided Laplace transform, the prior estimate of exponential decrease of traveling wave solutions is needed. For this aim, a new method is proposed, which can be applied to reaction-diffusion systems consisting of more than three equations.

1. Introduction

In the past and at present, cholera has been a serious threat to human health, which is an acute, diarrheal illness caused by infection of the intestine with the bacterium Vibrio cholera. An estimated 3–5 million cases and over 100,000 deaths occur each year around the world [1]. The cholera bacterium is usually found in water or food sources that have been contaminated by feces from a person infected with cholera. Cholera is most likely to be found and spread in places with inadequate water treatment, poor sanitation, and inadequate hygiene. Therefore, cholera outbreaks have been occurring in developing countries—for example, Iraq (2007-2008), Guinea Bissau (2008), Zimbabwe (2008-2009), Haiti (2010), Democratic Republic of Congo (2011-2012), and Sierra Leone (2012) [2].

Many mathematical models were proposed to understand the propagation mechanism of cholera, the earlier one of which was established by Capasso and Paveri-Fontana [3] to study the 1973 cholera epidemic in the Mediterranean region as follows: where and denote the concentrations of the pathogen and the infective populations, respectively. In addition, Codeço [4] investigated the role of the aquatic pathogen in dynamics of cholera through the following susceptible-infective-pathogen mode: where is the susceptible individuals. In this model, human is divided into two groups: the susceptible and the infective. As pointed out in [4–8], bacterium Vibrio cholera can spread by direct human-to-human and indirect environment-to-human modes. To understand the complex dynamics of cholera, model (2) is extended by [6, 9–15], and so forth.

In all previous models, the influences of space distribution of human on the transmission of cholera are omitted. Cholera usually spreads in spacial wave [16]. Cholera bacteria live in rivers and interact with the plankton on the surface of the water [17]. When individuals drink contaminated water and are infected, then they will release cholera bacteria through excretion [18]. Capasso et al. [19–22] developed model (1) by incorporating the bacterium diffusion in a bounded area and studied the existence and stability of solutions. To deeply investigate the interaction of transmission modes and bacterium diffusion, Bertuzzo et al. [23, 24] incorporated patchy structure into model (2) and supposed that the pathogen in water could diffuse among these patches. Furthermore, Mari et al. [25] studied the influence of diffusion of both human and pathogen on cholera dynamics through a patchy model.

Infectious case usually is found firstly at some location and then spreads to other areas. Consequently, the most important question for cholera is as follows: what is the spreading speed of cholera? However, the above spacial models mainly focus on the stability of solutions, not the spreading speed. Traveling wave solution is an important tool used to study the spreading speed of infectious diseases [26–28]. Based on Capasso’s model (1), Zhao and Wang [29], Xu and Zhao [30], Jin and Zhao [31], and Hsu and Yang [32] studied the influences of pathogen diffusion on the spread speed of cholera.

In above diffusive cholera models, diffusion of aquatic pathogen is neglected. In this paper, we investigate the effects of the disease-related death and aquatic pathogen in cholera epidemic dynamics by developing model (2). Based on model (2) and ignoring natural birth and death, a general diffusive epidemic cholera model incorporating the disease-related death and aquatic pathogen dynamics can be formulated as the following reaction-diffusion system: where and denote the concentrations of susceptible and infected individuals, respectively, and is the concentration of the infectious agents. is the disease-related death rate, denotes the contribution of each infected person to the concentration of cholera, and is the net death rate of the vibrio. is the environment-to-human transmission incidence. Similar to [15], we assume that satisfies(A1), , , for .From hypothesis (A1), we have .

In this paper, we study the traveling wave solutions of model (3). The formula for minimal wave speed is given. To prove the existence of traveling wave solutions for , an invariant cone is constructed and Schauder’s fixed point theorem is introduced. Schauder’s fixed point theorem is applied widely to prove the existence of traveling wave solutions (e.g., [26, 33, 34]). However, unlike Wang and Wu [34], the cone in our paper is bounded. Motivated by [34–37], we prove the nonexistence of traveling wave solutions for by two-sided Laplace transform, which was firstly introduced to prove the nonexistence of traveling wave solutions by Carr and Chmaj [37] and then was applied by [34–36]. To apply two-sided Laplace transform, the exponential decrease of traveling wave solutions is needed, which is proved in [34] by analysis method. However, it cannot be applied to our model due to the nonlinearity of cholera incidence. Therefore, in this paper, a new method is proposed to get the exponential decrease of traveling wave solutions, which is inspired by the proof of Stable Manifold Theorem in [38]. In addition, our method can be applied to reaction-diffusion systems consisting of more than three equations.

This paper is organized as follows. Section 2 is focused on the existence of traveling wave solutions. Firstly, the existence of traveling wave solutions for original system is proved to be equivalent to that of a new simple system. Then, two pairs of upper and lower solutions are constructed to get an invariant cone and Schauder’s fixed point theorem is applied for new system. Section 3 is devoted to the nonexistence of traveling wave solutions. For this aim, a new method is proposed to show the exponential decrease of traveling wave solutions and two-sided Laplace transform is used.

2. Existence of Traveling Wave Solutions

For convenience in discussing the model, we introduce dimensionless variables and parameters. Setting we obtain where Obviously, also satisfies assumption (A1) with being replaced by a new constant .

A traveling wave solution of system (5) is a nonnegative nontrivial solution of the form satisfying boundary condition where .

Before giving the main theorem, we introduce the equation for minimal wave speed: where The main result of this section is given as follows.

Theorem 1. Suppose . Then, there exists a positive constant which is the only positive root of (9). For any , system (5) has a traveling wave solution satisfying boundary condition (8) such that is nonincreasing in . Furthermore, one has that

To study the existence of traveling wave solutions, using constant variation method in the second equation of (5) gives Then, system (5) changes into

Lemma 2. is a traveling wave solution of system (13) satisfying boundary condition if and only if is a traveling wave solution of system (5) satisfying boundary condition (8), where

Proof. Assume is a traveling wave solution of system (13) satisfying boundary condition (14). Obviously, is a solution of system (5). To prove the necessity, it is enough to show that . Consider where the third equality is due to L’Hopital principal. The sufficiency is clear and is omitted. The proof is completed.

From Lemma 2, we only need to study traveling wave solutions of (13) satisfying boundary condition (14). Substituting traveling profile into system (13) yields the following equations: where denotes the derivative with respect to .

In the following, we will use Schauder’s fixed point theorem to prove the existence of traveling wave solutions. To achieve this goal, we firstly linearize the second equation of (18) at and obtain Substituting into (19), we get the characteristic equation that is, where , , and . To investigate distribution of roots of (21), denote and introduce the following lemma [39].

Lemma 3. (a) If , (21) has one real root and two nonreal complex conjugate roots.
(b) If , (21) has a multiple root and all its roots are real.
(c) If , (21) has three distinct real roots.

Then, we get the following lemma about the distribution of eigenvalues.

Lemma 4. Assume . Then, there exists a constant which is the only positive root of (9) such that(a)if , (21) has a negative real root and two nonreal complex conjugate roots with positive real parts;(b)if , (21) has a negative real root and a positive real multiple root;(c)if , (21) has a negative real root and two different positive real roots;(d)one can assume and let be the two positive roots of (21); then if .

Proof. Obviously, , , and . By Descartes’ rule of signs, has only one positive root such that for and for . Direct calculations show that . Since and , Descartes’ rule of signs shows that (21) has only one negative real root and Routh-Hurwitz criterion indicates that (21) has roots with positive real parts. Then, the combination of Lemma 3 and above analysis completes the proof of (a)–(c). (d) is obviously true, since is a cubic polynomial.

In this section, we always suppose and unless other conditions are specified. Denote to be the two positive roots of (21) and define where , .

Lemma 5. The function satisfies inequality for any .

Proof. Firstly, assume and, therefore, . Since satisfies (19) and for any , we have where .
Secondly, set , which implies . We have that The proof is completed.

Lemma 6. For sufficiently small and sufficiently large, the function satisfies for any .

Proof. Let be sufficiently large to ensure . When , then and the lemma is obviously true. Now, suppose . Then, and where . Let . Since we can find sufficiently small and sufficiently large such that Thus, the proof is completed.

Lemma 7. Let . Then, for large enough, the function satisfies for any .

Proof. It is clear that if and only if , that if and only if , and that if and only if . Let . When , then , , and Lemma 7 holds.
In this paragraph, assume . Then, , , and . To prove this lemma, it is enough to show where . Since we only need to show Since , then Since , inequality (34) satisfies if The proof is completed.

To apply Schauder’s fixed point theorem, we will introduce a topology in . Let be a positive constant which will be specified in the following. For , define We will find traveling wave solutions in the following profile set: Obviously, is closed and convex in . Firstly, we change system (18) into the following form: where , Suppose are the two roots of equation . Furthermore, define by In the remainder of this paper, it is always assumed that .

Lemma 8. Consider .

Proof. Let ; that is , . Then, we need to prove that
First of all, we have From Lemma 6 and system (39), we get where the second inequality is due to . Then, Therefore, we have proved .
Now, we study . If , then , which implies that since , . Assume . From Lemma 7 and system (39), it is clear that which implies where the final inequality is due to and . In conclusion, for any .
Similarly, we can show and the proof is completed.

Lemma 9. Map is continuous with respect to the norm in .

Proof. For , we have where , and is between and . Therefore, we have Then, Thus, when , we get
When , it follows that Consequently, we have proved that for any , where That is, In conclusion, is continuous with respect to the norm in .
In addition, consider . Firstly, we have where . Therefore, If , it holds that If , we have Consequently, we conclude that where Thus, is continuous with respect to the norm in . The proof is completed.

Lemma 10. Map is compact with respect to the norm in .

Proof. For any , it is clear that Since , we have which implies that . Since , we get Then, which implies that. Consequently, and are bounded, which shows that is uniformly bounded and equicontinuous with respect to the norm in .
Furthermore, for any positive integer , define Obviously, for fixed , is uniformly bounded and equicontinuous with respect to the norm in , which implies that is a compact operator. Since we have that Similarly, we can prove that when . Thus, when . By Proposition 2.1 in Zeilder [40], we have that converges to in with respect to the norm . Consequently, is compact with respect to the norm . The proof is completed.

Proof of Theorem 1. Combination of Schauder’s fixed point theorem and Lemmas 8, 9, and 10 shows that there exists a nonnegative traveling wave solution such that when . By the definition of and , there is a such that and when . Therefore, if , we have that which implies that . It is clear that is nonincreasing in .
From Lemma 2 and system (18), it follows that By the first and second equation of (69), we have Since the traveling wave solution and is the fixed point of the operator , thus . The proof of Lemma 10 shows that is bounded. Using L’Hopital principal shows that Integrating the third equation of (69) from to yields Then, the boundedness of and implies , which shows that and .
Next, we prove , which is motivated by Wang and Wu [34]. Define It is clear that satisfies , and . Denote , which satisfies Due to , it follows that which implies that is nondecreasing in . Since and , we have for all . The proof is completed.

3. Nonexistence of Traveling Wave Solutions

In this section, we give the conditions on which system (5) has no traveling wave solutions.

Theorem 11. (I) Assume . Then, for any , system (5) has no nonnegative nontrivial traveling wave solutions satisfying boundary condition (8).
(II) Suppose . Then, for any , system (5) has no traveling wave solutions satisfying boundary condition (8).

To prove (I) by two-sided Laplace transform, the following lemma about the exponential decrease of traveling wave solutions is needed.

Lemma 12. Assume . If is a nonnegative nontrivial traveling wave solution of (5) satisfying boundary condition (8), there exists a positive constant such that

Proof. Substituting traveling wave profile into the second and third equations of (5), it follows that Setting , we have Furthermore, system (80) can be written as where and . Obviously, . Since and is bounded in , for any small constant , there exists such that for all .
Since is the solution of (79) if and only if is the solution of (79) for any and they satisfy the same boundary condition at , we can select large such that (84) holds with and being replaced by . Therefore, we suppose .
Calculations show that the eigenfunction of is defined by (21). Then, by Lemma 4, there exists constant matrix such that where and the eigenvalues , of matrix have positive real parts. Setting , the system (81) then has the form Obviously, and satisfies (84). Denote Then, , , and It is not difficult to see that we can choose positive constants and such that
It is clear that (86) is equivalent to where Since , then multiplying the first equality of (90) by and setting yield or Thus, (90) has the form where .
To study the properties of , we need to construct a functional sequence. Define for . Obviously, we have for any . Assume that the induction hypothesis holds for and . Then, using the condition (84) satisfied by the function , it follows from the induction hypothesis that for we have Setting yields It then follows by induction that (97) holds for all and . Thus, for any and , we have which implies that is a Cauchy sequence of continuous functions. It follows that uniformly for all . Taking the limit of both sides of (95), it follows from the uniform convergence that the continuous function satisfies the integral equation (94). Then, the uniqueness of solution of (94) implies . Consequently, (97) and show that for any . Since , it follows . Thus, and the proof is completed.

Proof of Theorem 11(I). Assume (5) has a traveling wave solution satisfying boundary condition (8). From the first equation of (18), it follows By Lemma 12, there exists a positive constant such that For , define two-sided Laplace transform as follows: Obviously, there exists and is increasing in such that satisfying or .
Furthermore, we have The second equality of (18) can be rewritten as where Denote , from which it follows that . Since satisfies hypothesis (A1), we have for any .
Taking two-sided Laplace transform of (109), we have If , it is clear that . Since Lemma 12, inequality (106), and imply . However, Lemma 4 shows that , which contradicts (111). Hence, we have . Furthermore, (111) can be rewritten as Since , it is impossible that the above equality holds, again a contradiction.

Proof of Theorem 11(II). Suppose is a traveling wave solution of system (5) satisfying boundary condition (8). Then, by the proof of Theorem 1, it is clear that . Furthermore, which is a contradiction. The proof is completed.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors are supported by the Fundamental Research Funds for the Central Universities (XDJK2012C042) and NSFC Grants (11201380).

References

Centers for Disease Control and Prevention, USA, http://www.cdc.gov/cholera/general/.
World Health Organization, http://www.who.int/csr/don/archive/disease/cholera/en/.
V. Capasso and S. L. Paveri-Fontana, “Mathematical model for the 1973 cholera epidemic in the European Mediterranean region,” Revue d'Epidemiologie et de Sante Publique, vol. 27, no. 2, pp. 121–132, 1979. View at: Google Scholar
C. T. Codeço, “Endemic and epidemic dynamics of cholera: the role of the aquatic reservoir,” BMC Infectious Diseases, vol. 1, article 1, 2001. View at: Publisher Site | Google Scholar
K. T. Goh, S. H. Teo, S. Lam, and M. K. Ling, “Person-to-person transmission of cholera in a psychiatric hospital,” Journal of Infection, vol. 20, no. 3, pp. 193–200, 1990. View at: Publisher Site | Google Scholar
J. H. Tien and D. J. D. Earn, “Multiple transmission pathways and disease dynamics in a waterborne pathogen model,” Bulletin of Mathematical Biology, vol. 72, no. 6, pp. 1506–1533, 2010. View at: Publisher Site | Google Scholar
A. A. Weil, A. I. Khan, F. Chowdhury et al., “Clinical outcomes in household contacts of patients with cholera in Bangladesh,” Clinical Infectious Diseases, vol. 49, no. 10, pp. 1473–1479, 2009. View at: Publisher Site | Google Scholar
J. R. Andrews and S. Basu, “Transmission dynamics and control of cholera in Haiti: an epidemic model,” The Lancet, vol. 377, no. 9773, pp. 1248–1255, 2011. View at: Publisher Site | Google Scholar
D. M. Hartley, J. G. Morris Jr., and D. L. Smith, “Hyperinfectivity: a critical element in the ability of V. cholerae to cause epidemics?” PLoS Medicine, vol. 3, no. 1, pp. 63–69, 2006. View at: Publisher Site | Google Scholar
M. A. Jensen, S. M. Faruque, J. J. Mekalanos, and B. R. Levin, “Modeling the role of bacteriophage in the control of cholera outbreaks,” Proceedings of the National Academy of Sciences of the United States of America, vol. 103, no. 12, pp. 4652–4657, 2006. View at: Publisher Site | Google Scholar
Z. Mukandavire, S. Liao, J. Wang, H. Gaff, D. L. Smith, and J. G. Morris Jr., “Estimating the reproductive numbers for the 2008-2009 cholera outbreaks in Zimbabwe,” Proceedings of the National Academy of Sciences of the United States of America, vol. 108, no. 21, pp. 8767–8772, 2011. View at: Publisher Site | Google Scholar
R. P. Sanches, C. P. Ferreira, and R. A. Kraenkel, “The role of immunity and seasonality in cholera epidemics,” Bulletin of Mathematical Biology, vol. 73, no. 12, pp. 2916–2931, 2011. View at: Publisher Site | Google Scholar
Z. Shuai and P. van den Driessche, “Global dynamics of cholera models with differential infectivity,” Mathematical Biosciences, vol. 234, no. 2, pp. 118–126, 2011. View at: Publisher Site | Google Scholar
Z. Shuai, J. H. Tien, and P. V. V. Driessche, “Cholera Models with hyperinfectivity and temporary immunity,” Bulletin of Mathematical Biology, vol. 74, no. 10, pp. 2423–2445, 2012. View at: Google Scholar
J. P. Tian and J. Wang, “Global stability for cholera epidemic models,” Mathematical Biosciences, vol. 232, no. 1, pp. 31–41, 2011. View at: Publisher Site | Google Scholar
A. Mutreja, D. W. Kim, N. R. Thomson et al., “Evidence for several waves of global transmission in the seventh cholera pandemic,” Nature, vol. 477, no. 7365, pp. 462–465, 2011. View at: Publisher Site | Google Scholar
M. Emch, C. Feldacker, M. S. Islamand, and M. Ali, “Global climate and infectious disease: the cholera paradigm,” Science, vol. 274, no. 5295, pp. 2025–2031, 1996. View at: Google Scholar
S. M. Faruque, M. J. Islam, Q. S. Ahmad et al., “Self-limiting nature of seasonal cholera epidemics: role of host-mediated amplification of phage,” Proceedings of the National Academy of Sciences of the United States of America, vol. 102, no. 17, pp. 6119–6124, 2005. View at: Publisher Site | Google Scholar
V. Capasso and K. Kunisch, “A reaction-diffiusion system arising in modelling man- environment diseases,” Quarterly of Applied Mathematics, vol. 46, no. 3, pp. 431–450, 1988. View at: Google Scholar
V. Capasso and L. Maddalena, “Convergence to equilibrium states for a reaction-diffusion system modelling the spatial spread of a class of bacterial and viral diseases,” Journal of Mathematical Biology, vol. 13, no. 2, pp. 173–184, 1981. View at: Google Scholar
V. Capasso and R. E. Wilson, “Analysis of a reaction-diffusion system modeling man-environment-man epidemics,” SIAM Journal on Applied Mathematics, vol. 57, no. 2, pp. 327–346, 1997. View at: Google Scholar
S. Aniţa and V. Capasso, “Stabilization of a reaction-diffusion system modelling a class of spatially structured epidemic systems via feedback control,” Nonlinear Analysis: Real World Applications, vol. 13, no. 2, pp. 725–735, 2012. View at: Publisher Site | Google Scholar
E. Bertuzzo, S. Azaele, A. Maritan, M. Gatto, I. Rodriguez-Iturbe, and A. Rinaldo, “On the space-time evolution of a cholera epidemic,” Water Resources Research, vol. 44, no. 1, Article ID W01424, 2008. View at: Publisher Site | Google Scholar
E. Bertuzzo, R. Casagrandi, M. Gatto, I. Rodriguez-Iturbe, and A. Rinaldo, “On spatially explicit models of cholera epidemics,” Journal of the Royal Society Interface, vol. 7, no. 43, pp. 321–333, 2010. View at: Publisher Site | Google Scholar
L. Mari, E. Bertuzzo, L. Righetto et al., “Modelling cholera epidemics: the role of waterways, human mobility and sanitation,” Journal of the Royal Society Interface, vol. 9, no. 67, pp. 376–388, 2012. View at: Publisher Site | Google Scholar
O. Diekmann, “Thresholds and travelling waves for the geographical spread of infection,” Journal of Mathematical Biology, vol. 6, no. 2, pp. 109–130, 1978. View at: Google Scholar
M. Lewis, J. Rencławowicz, and P. van den Driessche, “Traveling waves and spread rates for a west nile virus model,” Bulletin of Mathematical Biology, vol. 68, no. 1, pp. 3–23, 2006. View at: Publisher Site | Google Scholar
J. Radcliffe and L. Rass, “The asymptotic speed of propagation of the deterministic non-reducible n-type epidemic,” Journal of Mathematical Biology, vol. 23, no. 3, pp. 341–359, 1986. View at: Google Scholar
X.-Q. Zhao and W. Wang, “Fisher waves in an epidemic model,” Discrete and Continuous Dynamical Systems B, vol. 4, no. 4, pp. 1117–1128, 2004. View at: Google Scholar
D. Xu and X.-Q. Zhao, “Bistable waves in an epidemic model,” Journal of Dynamics and Differential Equations, vol. 16, no. 3, pp. 679–707, 2004. View at: Google Scholar
Y. Jin and X.-Q. Zhao, “Bistable waves for a class of cooperative reaction-diffiusion systems,” Journal of Biological Dynamics, vol. 2, no. 2, pp. 196–207, 2008. View at: Google Scholar
C.-H. Hsu and T.-S. Yang, “Existence, uniqueness, monotonicity and asymptotic behaviour of travelling waves for epidemic models,” Nonlinearity, vol. 26, no. 1, pp. 121–139, 2013. View at: Google Scholar
W.-T. Li, G. Lin, and S. Ruan, “Existence of travelling wave solutions in delayed reaction-diffusion systems with applications to diffusion-competition systems,” Nonlinearity, vol. 19, no. 6, pp. 1253–1273, 2006. View at: Publisher Site | Google Scholar
Z.-C. Wang and J. Wu, “Travelling waves of a diffiusive KermackCMcKendrick epidemic model with non-local delayed transmission,” Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences, vol. 466, no. 2113, pp. 237–261, 2010. View at: Google Scholar
Z.-C. Wang, W.-T. Li, and S. Ruan, “Traveling fronts in monostable equations with nonlocal delayed effects,” Journal of Dynamics and Differential Equations, vol. 20, no. 3, pp. 573–607, 2008. View at: Publisher Site | Google Scholar
Z.-C. Wang, W.-T. Li, and S. Ruan, “Entire solutions in bistable reaction-diffusion equations with nonlocal delayed nonlinearity,” Transactions of the American Mathematical Society, vol. 361, no. 4, pp. 2047–2084, 2009. View at: Publisher Site | Google Scholar
J. Carr and A. Chmaj, “Uniqueness of travelling waves for nonlocal monostable equations,” Proceedings of the American Mathematical Society, vol. 132, no. 8, pp. 2433–2439, 2004. View at: Publisher Site | Google Scholar
L. Perko, Diffierential Equations and Dynamical Systems, Springer, New York, NY, US, 2001.
R. S. Irving, Integers, Polynomials, and Rings, Springer, New York, NY, USA, 2004.
E. Zeilder, Nonlinear Functional Analysis and Its Applications I, Springer, New York, NY, USA, 1986.

Copyright

Copyright © 2014 Tianran Zhang and Qingming Gou. 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.

The Scientific World Journal

Dynamics of Nonlinear Systems