Abstract

To investigate the spreading speed of cholera, Codeço’s cholera model (2001) is developed by a reaction-diffusion model that incorporates both indirect environment-to-human and direct human-to-human transmissions and the pathogen diffusion. The two transmission incidences are supposed to be saturated with infective density and pathogen density. The basic reproduction number is defined and the formula for minimal wave speed is given. It is proved by shooting method that there exists a traveling wave solution with speed for cholera model if and only if .

1. Introduction

Cholera has been a serious threat to human health in the past and at present, 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 to spread in places with inadequate water treatment, poor sanitation, and inadequate hygiene. Therefore, cholera outbreaks have occurred 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].

To understand the propagation mechanism of cholera, many mathematical models were proposed, whose earlier one 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 model: where is the susceptible individuals. In this model, human is divided into two groups: the susceptible group and the infective group. As pointed out in [48], 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 [815] 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 spatial 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, they will release cholera bacteria through excretion [18]. Capasso et al. [1923] 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. [24, 25] incorporated patchy structure into model (2) and supposed that pathogen in water could diffuse among these patches. Furthermore, Mari et al. [26] studied the influence of diffusion of both human and pathogen on cholera dynamics through a patchy model.

Infectious case is usually found firstly at some location and then spreads to other areas. Consequently, the most important question for cholera is what the spreading speed of cholera is. However, the above spatial 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 [2729]. Based on Capasso’s model (1), Zhao and Wang [30], Xu and Zhao [31], Jin and Zhao [32], and Hsu and Yang [33] studied the influences of pathogen diffusion on the spread speed of cholera.

The studies of traveling wave solutions of Capasso’s model (1) incorporating pathogen diffusion provide insight into the spreading speed of cholera. However, some pieces of information are omitted, such as the interaction of direct human-to-human and indirect environment-to-human transmissions. In this paper, a reaction-diffusion model with pathogen diffusion and both transmission paths is proposed by developing Codeço’s model (2). Based on model (2) and ignoring the disease-related death, a general diffusive cholera model 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 total human population, stands for the natural birth and death rate, denotes the contribution of each infected person to the concentration of cholera, and is the net death rate of vibrio cholera. and are the human-to-human and environment-to-human transmission incidences, respectively. Similar to [10], we assume that and satisfy (A1), , ;(A2), , , , and is strictly monotonously increasing in .It is easy to conclude that , , and and are nonincreasing. Obviously, hypotheses (A1) and (A2) imply that the two transmission paths are saturated. In Tian and Wang [10], and have the following expressions: Obviously, as a special case, such selections satisfy (A1) and (A2).

Shooting method is very important in proving the existence of traveling wave solutions, which was proposed by Dunbar [34, 35] and was applied to many models (e.g., [3640]). In this paper, the existence of traveling wave solutions of system (3) will be proved by shooting method and the formula for minimal wave speed will be given.

This paper is organized as follows. In next section, the main theorem and the formula for minimal wave speed will be given. In Section 3, the nonexistence of the traveling wave solutions for is proved by geometric method. Section 4 is devoted to shooting arguments and the construction of Wazewski set. In Section 5, we prove the existence of traveling wave solutions for and then give the existence of traveling wave solution for by limit arguments. The final section is devoted to the simulations.

2. Main Results

For convenience, we introduce dimensionless variables and parameters. By setting model (3) has the form where and .

Denote , which is the basic reproduction number of (6). Then hypotheses (A1) and (A2) imply that system (6) has two nonnegative constant solutions and if and only if , where is the only one positive root of equation and . Biologically, corresponds to disease-free equilibrium and corresponds to endemic equilibrium. To study the spreading wave of cholera, it is assumed that holds in this paper; that is

A traveling wave solution of system (6) with speed is a nonnegative solution of the form

Substituting traveling profile into system (6) yields the following equations: where denotes . To investigate invasion question by cholera, we will study the positive solutions of (10) such that

Before giving the main theorem, we introduce the equation for minimal wave speed where

Theorem 1. There exists a constant which is the greatest positive root of (12). When , system (6) has a traveling wave solution satisfying boundary condition (11). When , system (6) has no traveling wave solutions satisfying boundary condition (11).

3. Nonexistence of Traveling Wave Solutions for

From (10), we have

Consequently, if , then

Hence, the traveling profile with boundary condition (11) must satisfy

Therefore, to study traveling wave solutions we assume (16) satisfies. Setting in system (10) and noticing (16), it follows

If , then by system (10). Therefore, we suppose for any ; that is, .

Obviously, system (17) has two equilibria and . A profile solution of (10) which satisfies boundary condition (11) corresponds to the positive solution of system (17) which satisfies where . Therefore, to study the solutions of (10), it is sufficient to study those of system (17) satisfying boundary condition (18).

Firstly, we investigate the dynamics near . Simple calculations show that the characteristic equation of the linearization of system (17) at is where

Because (19) has a negative real root, which is denoted by . Let and be the other two eigenvalues of (19) and suppose that . To investigate the distribution of roots of (19), denote and introduce the following lemma [41].

Lemma 2.    If , (19) has one real root and two nonreal complex conjugate roots.
   If , (19) has a multiple root and all its roots are real.
   If , (19) has three distinct real roots.
Direct calculations show that , where is defined by (12).

Lemma 3.   The real parts of and are positive.
  Assume . Then, there exists which is the only positive root of . When , , and are real. When , , and are complex and nonreal.
  Assume that . Then, there exist two positive constants which are all positive roots of . and are complex and nonreal if and only if . If , then ; if , then , where .
   if and only if or .

Proof. Suppose is the root of (19). Substituting into (19) and comparing real and imaginary parts show that and . Since , then . However, it is impossible that and by the expressions of and . Therefore, the real parts of and are not zero. Furthermore, since it is impossible that and , Routh-Hurwitz theorem implies that it is impossible that the real parts of both and are negative. Consequently, there are two cases: (i) and are complex conjugate roots with positive real parts; (ii) and are real and at least one is positive. However, Descartes’ rule of signs shows that the number of positive roots of (19) is zero or two. Thus, if case (ii) is true, both of and are real and positive. Therefore, (a) is proved.
In this paragraph, we consider the case . Firstly, suppose that . Obviously, and . By the expression of , we have since . Now, assume ; that is, . Then. , , , and . Then, if , Descartes' rule of signs shows that there exists which is the only positive root of , where for and for . Using Lemma 2 completes the proof of (b).
Suppose that in this paragraph and, thus, . Calculations show that and that has two roots and , where and . By letting and using trivial calculations, we get (see Figure 1)
Therefore, if , then . Since is the only minimum-value point of , and then for any and both of and are not real. Lemma 2 shows that . Thus, since and , there exist two positive roots for equation such that . Then, using (25) and Lemma 2 completes the proof of (c) and (d).

Direct calculations show that corresponding eigenvectors of eigenvalue are where . Since and thus

Then, we have the following lemma.

Lemma 4. If , there exist no traveling wave solutions which satisfy boundary condition (11).

Proof. Assume that and . Then, (b) of Lemma 3 implies that and are complex conjugate eigenvalues and there exits locally unstable manifold and locally stable manifold . If a solution of (17) tends to when , then it will be spiral on . By the structures of and , at some time , which shows that there exist no traveling wave solutions departing from .
Suppose that . If , (c) of Lemma 3 shows that and are complex conjugate eigenvalues and similar arguments to that of previous paragraph finish the proof. If , (c) of Lemma 3 shows that and are real; however, . If a solution of (17) tends to when , structures of and indicate that there is an such that . The proof is completed.

From Section 4 to Section 5.2, we suppose that , which implies .

4. Shooting Method and Wazewski Set

To prove the existence of traveling wave, shooting method developed by Dunbar [34] is used. Firstly, we give the shooting arguments.

Consider the differential equation where from to satisfies Lipschitz condition about . Let denote the unique solution of (29) with initial value . It is convenient to give the notations and . To describe the shooting method (or Wazewski theorem), some definitions are necessary.

Definition 5. (a) For , define immediate exit set of as
(b) For , let .
(c) Given , define exit time of by
Then, Wazewski theorem is formulated as follows.

Lemma 6 (see [34]). Suppose that (1)if and , then .(2)If , and , then there exists an open set about disjoint from .(3)If , is compact and intersects a trajectory of (29) only once.
Then, the mapping is a homeomorphism from to its image on .

A set satisfying conditions (1) and (2) of Lemma 6 is called a Wazewski set. In the following, we first construct the Wazewski set . Fundamental idea to construct a Wazewski set is that the characteristic vectors corresponding eigenvalues with positive real parts should be removed from and that those characteristic vectors corresponding eigenvalues with negative real parts should be included. Therefore, we set where

It is obvious that . Firstly, we give the construction of , which is described in Figure 2.

Lemma 7. The construction of is as follows: where .

Proof. It is enough to analyze the behavior of solution on . We only study and omit the proof of since the analysis of is similar to that of and is simpler. In the process of this proof, we use some notations to simplify the proof. Set
From hypotheses (A1) and (A2), we find that and are monotonously decreasing, is strictly monotonously decreasing for , and is the only positive root of . The set is classified into two cases according to variable . (a) Case . This case is classified as follows.(1) Case . Then and the solution of (17) will enter .(2) Case . Then The solution of (17) will enter .(3) Case . Then The solution of (17) will enter .(4) Case . Then and the solution of (17) will enter .(5) Case and . The solution of (17) will enter .(6) Case . Then , and . Therefore, the solution of (17) will enter .(b)  Case . This case is classified as follows.(1) Case . Then and the solution of (17) will enter .(2) Case . Then , . The solution of (17) will enter .(3) Case . Then , , and . The solution of (17) will enter .(4) Case . Then and . The solution of (17) will enter .(5) Case . In this case, is equilibrium and is constant.(6) Case and . Then and . The solution of (17) will enter .(7) Case . Then is equilibrium and is constant.
The proof is completed.

5. Existence of Traveling Wave Solution for

In this section, we prove the existence of traveling wave solution for . Firstly, we study the behaviors of solutions near .

5.1. Behaviors of Solutions Near

Lemma 8. Suppose is a solution of (17) satisfying initial conditions where . Then, for every , we have

Proof. From Lemma 3, we have . To finish the proof, it is sufficient to prove that the set is positively invariant. It is obvious that where
Suppose that . Then, and
The last inequality is given since . Suppose that . If , then
If , we have Consequently, the solution of system (17) departing from cannot intersect . If , then . Since is equilibrium, in summary, is positive invariant.

Since , stable manifold theorem implies that there exists a one-dimensional strong unstable manifold tangent to at such that the point on near can be expressed by Furthermore, there is a two-dimensional unstable manifold tangent to span at such that near can be expressed by

Lemma 9. Suppose that is a solution of (17) such that for small . Then, will leave and enter .

Proof. Obviously, satisfies initial condition (39) by the structure of , and Lemma 8 implies ( means that and , ) for every .
Furthermore, Lemma 8 shows that , implying . Since , it follows . Suppose that for every . Then for large since and are strictly monotonous increasing with respect to and , respectively. Thus, we have that , contradicting for any . Therefore, there exists such that . Without losing generality, let . Obviously, we have . If , thenwhich is a contradiction. Therefore, and . Then, the construction of shows that will leave and enter .

Let be a small circle on centered at . Then, points on can be expressed in terms of local coordinate by where , , and is chosen such that lies on with . Then, stable manifold theorem shows that when . Denote .

Lemma 10. There exists a such that and that for .

Proof. From (51), we have where , , and since . Therefore, and imply that . Obviously, for any . However,
Then, equality , together with the last of (55), reveals ; that is, . For , the first and second equalities of (55) imply that where since and .

Let

By Lemma 10, is an arc of circle, , and the solution of (17) with initial value being the endpoint will enter since . From Lemma 9, the solution of (17) with initial value being the endpoint will enter .

5.2. Traveling Wave Solution for

Lemma 11. Let be a solution of (17) such that . If for any , then for any , where and .

Proof. Set . Suppose the conclusion is false; that is, . Obviously, and where
In Figure 3, we find , , , , , , and .
Since , thus . If , we have because for and . However, which is a contradiction. Therefore, . If , thencontradicting . If , then since , contradicting . If , then which is a contradiction. In conclusion, . If , then and for any . Hence, for any , which implies that . From this contradiction we find . Because is a constant solution, we get . In summary, and . The proof is completed.

Lemma 12. There exists a point such that the solution of (17) with initial value being will stay in for any .

Proof. It is sufficient to prove . Suppose that . Firstly, we verify Conditions and of Lemma 6. Condition of Lemma 6 is valid since is closed.
Suppose , and . Then, and since . The structure of implies that , , and . By the proof of Lemma 11, we have that for . Therefore, and . Condition of Lemma 6 holds.
Lemma 6 shows that is homeomorphic to . Since and is disconnected, we have that is disconnected, contradicting the connection of . Thus, and the proof is completed.

Lemma 13. Let . Then, there exists a positive solution of (17) such that

Proof. By Lemma 12 there exists a point such that the solution of (17) with initial value being will stay in for any . Furthermore, Lemma 11 shows for any . Stable manifold theorem implies that for any and . Therefore, is a positive solution.
To complete the proof, it is sufficient to show that . By Lemma 11, we know that for any since remains in for all . Because , then the limit of exists; that is, and . Suppose that . The first equation of (17) shows that for large , which implies that there is an such that . This is a contradiction, and thus . From the first equation of (17), we have where is the only positive root of algebra equation At the same time, the third equation of (17) implies and or . It is impossible that due to the boundedness of . In conclusion, the limit exists and is finite. By [42], must be equilibrium. Since , then .

Noticing the relation of systems (17) and (10) completes the proof of Theorem 1 for case .

5.3. Traveling Wave Solution for

Firstly, suppose and let be the traveling wave solution of (17). Then, Lemma 11 implies that for all . From the proof of Lemma 13, we find . Therefore, for all , we have where

Let be a sequence such that for any and . Set and Then, for any .

Lemma 13 shows that there is a positive solution for system satisfying boundary condition (62) such that for any .

Lemma 14. Let . Then, there exists a traveling wave solution for system (6) satisfying boundary condition (11).

Proof. It is sufficient to prove that there exists a positive solution of (17) satisfying boundary condition (62).
Firstly, we show that sequences , , , , , and are uniformly bounded and equicontinuous. The idea of Lemma 11 in [34] is used. Obviously, , , and are uniformly bounded since for any . Because is the solution of (67), , , and are also uniformly bounded. Since where , then is equicontinuous. Similarly, and are also equicontinuous. By differentiating the equations of (67) and using the previous bounds, we can get that , , and are uniformly bounded, and hence , , and are equicontinuous.
The previous paragraph and Arzelà-Ascoli theorem imply that there exist subsequences, again denoted by , , and and functions , , and such that uniformly on compact subsets of , thus pointwise on . Same arguments imply that , , and are also uniformly convergent on compact subsets of and pointwise convergent on . Consequently, we get Since is the solution of (67), then is the solution of (17) for and , where is the closer of . Because system (67) is autonomous and satisfies boundary condition (62), we can assume that for all ; thus, . Then, similar to the proof of Lemma 13, we have that the solution satisfies boundary condition (62).

6. Simulations

In this section, we present some simulations to confirm the theoretical results. Set and assign numerical values to parameters as follows:

Obviously, such selection for and satisfies (A1) and (A2). Then, the traveling wave solution is described in Figure 4.

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 (Grants nos. XDJK2012C042 and SWU113048) and NSFC (Grant no. 11201380).