Mathematical Model of Schistosomiasis under Flood in Anhui Province

and Applied Analysis 3 basic reproduction numberR 0 provides a threshold condition for parasite extinction. Theorem 1. There exist at most two equilibria: (i) ifR 0 ≤ 1, system (2) has a disease free equilibriumE 0 = (A x /μ x , 0, A y /μ y , 0, 0); (ii) if R 0 > 1, system (2) has two equilibria, the disease free equilibrium E 0 and the unique endemic equilibrium


Introduction
As we know, schistosomiasis is a serious water-borne disease. It is not easy to control because of many reasons such as flood. Many reports have shown that flood leads to a serious outbreak of schistosomiasis [1][2][3]. During the flood period there are a lot of people that come into contact with contaminated water, which may lead to the fact that a lot of people are infected by schistosome [1][2][3]. In China, Anhui province often encounters floods; in particular in 1998 the flood was one of the most serious flood [1]. Based on the observation data from Tongcheng Schistosomiasis Control Station in Anhui province (Figure 1), we can see that the number of patients and the area of snails increase by a big margin after 1998 in Tongcheng city in Anhui province. Although people know the phenomenon that schistosomiasis will be serious after flood, people do not know the reason and there are only some live reports. Hence, it is necessary to investigate theoretically the effect of flood on the schistosomiasis transmission.
After flood the infected human by cercaria will have an incubation period to become an infectious human. In fact, it is about five weeks from the time of cercaria penetration through skins of human to the time when eggs are discharged [4]. Adult schistosomes in human are capable of producing eggs for a number of years [5]. This leads to breakout of schistosomiasis in many places after flood. For example, the catastrophic flood in 1998 brought a serious impact on the prevalence of schistosomiasis in Anhui province from 1998 to 2000 [1]. Furthermore, the data from Tongcheng Schistosomiasis Control Station ( Figure 1) and the report of Ge et al. [1] both show that schistosomiasis is more serious in three years after flood than in the flood year. This phenomenon is called the delayed effect of flood [1]. In this paper we want to investigate how flood affects the dynamical behavior of schistosomiasis.
Many schistosomiasis models have involved many aspects such as drug-resistant, age-structure, incubation period of snail, and chemotherapy [6][7][8][9][10]. Their results imply that many factors affect the transmission of schistosomiasis. However, there are few mathematical models considering the effect of flood in previous papers.
To study the delayed effect of flood, we use a time delay to reflect the incubation period in the infected human. We modified the model in [11]. Distribute human into susceptible ( ) and infectious ( ) and snails into susceptible ( ), preshedding ( ), and infectious ( ). The model in [11] where is the recruitment rate of human, is the death rate of human, is the disease-induced death rate of human, is the transmission rate from infectious snails to susceptible human, is the recruitment rate of snail host, is the death rate of snail host, is the disease-induced death rate of snail host, is the transmission rate from infectious human to susceptible snails, and is the translate rate from infected and preshedding snails to shedding snails. In the model, we have studied the stability of equilibria and preferable control strategies.
The goal of this paper is to study the impact of flood on the basic reproduction number and the dynamics of the schistosomiasis transmission. This paper is organized as follows. In Section 2 we establish a schistosomiasis model with a time delay and define the basic reproduction number 0 . The stability of the disease free equilibrium is obtained in Section 3. We devote Section 4 to the Hopf bifurcation analysis. Section 5 examines mathematical results by numerical simulations.

The Delayed Model
By incorporating a time delay in human, we have the following model: where is the incubation period in the infected human, that is, the time from cercaria penetration through skins to the time when eggs are discharged. Define the basic reproduction number according to biological meanings: These quantities have a clear biological interpretation. Consider the case when an infectious snail is introduced into a purely susceptible people population with size / . The size of susceptible people who become infectious people per unit time is ( / ). 1/( + ) is the mean infective period of the infectious people and − represents the survived rate of people during his infection. On the other hand, infectious people can infect ( / ) susceptible snails which should get through the latent time where the rate of transmission is and then infective period is 1/( + )( + ). Thus, 0 gives the total number of secondary infectious snails produced by a typical infected snail during its entire period of infectiousness in a completely susceptible population. The following section shows that the Abstract and Applied Analysis 3 basic reproduction number 0 provides a threshold condition for parasite extinction.

Stability Analysis of 0
In this section, we will analyze the stability of the disease free equilibrium 0 of the delayed model (2) in the two cases: 0 < 1 and 0 > 1.

Theorem 2.
The disease free equilibrium 0 of the system (2) is locally asymptotically stable if 0 < 1 and unstable if 0 > 1.
From [12,Lemma 3.31], if ≥ 0, > 0, ≥ 0, then (10) has no positive real roots. This implies that (7) does not have positive solution since 0 < 1. Therefore, (5) does not have purely imaginary roots. Consequently, the real parts of all eigenvalues of 0 are negative for all positive . This indicates that the disease free equilibrium 0 is locally asymptotically stable if 0 < 1.

Hopf Bifurcation Analysis
In this section, we turn to the study of the stability of the endemic equilibrium when 0 > 1. Notice that 0 > 1 is equivalent to The characteristic equation of is In the following, it can be shown that (12) does not have nonnegative real roots for > 0. Let Note that̃3 > 0,̃4 > 0,̃5 > 0 for all ≥ 0 and > 0. We rewrite (12) in the following form: It is easy to see that the left-hand side in (15) is positive while the right-hand side is negative for all ≥ 0. Then (12) does not have nonnegative real solutions. Now we consider whether or not (12) has purely imaginary solutions. Suppose = , > 0 for some > 0, is a root of (12). Then we have From (17), we obtain that is, Let = 2 again; we obtain Abstract and Applied Analysis 5 Because (20) is very complex, the roots cannot easily be found. However, we know there are positive roots in some conditions. For example, if 5 < 0, then (20) has at least a positive root, say 0 , and (19) has at least a positive root 0 = √ 0 . Consequently, the endemic equilibrium may lose stability and lead to oscillations because the time delay > 0. In this case, we will do bifurcation analysis by as bifurcation parameter in the following.
Theorem 3. If 0 > 1, 5 < 0 and 0 is the largest positive simple root of (19), a Hopf bifurcation occurs around the endemic equilibrium of the delayed model (2).

Numerical Simulations
Based on the observation data from the investigation of Tongcheng Schistosomiasis Control Station in Anhui province, we estimated transmission rates in our model. Also according to the previous papers [7-9, 11, 13], we choose the parameter values in Table 1. Thus, 0 > 1, * = 327, and 5 < 0 when = 0.1.
Note that the bifurcation parameter 0 = 3 at this time. We performed some simulations and obtained Figure 2. From Figure 2 we can see that Hopf bifurcation can occur when 0 = 3. This implies that schistosomiasis will break out in about three years after flood. It is also in accord with the investigation of Tongcheng Schistosomiasis Control Station. This phenomenon is also in accord with the report of the whole Anhui province [1]. From our theoretical results and

Discussion
In this paper, based on the observation data in Tongcheng Schistosomiasis Control Station in Anhui province we have modified our previous model by including a time delay that describes the incubation period of schistosome within infected human. We define the basic reproduction number 0 according to biological meanings and give the existence of the disease free equilibrium and the endemic equilibrium. We find that, if 0 < 1, then the disease free equilibrium is locally asymptotically stable. However, the stability of the unique endemic equilibrium may be changed under some condition even if the basic reproduction number is larger than one. The results imply that the time delay can destabilize the system and periodic solutions can arise by Hopf bifurcation.
Numerical simulations imply that schistosomiasis will break out in about three years after flood. Furthermore the observation data show that schistosomiasis will be the most serious in about the third year after flood. From Figure 1, we can see that the number of patients and snails did not greatly change in 1998 and 1999. However, in 2001 the number of patients became about 5 times that of 1998, and the area of snails became about two times that of 1998. In our simulations, there is a little difference. Our results are higher than the observation data. We think the reason may be that after flood the government dispatched a large number of manpower and material resources to control the spread of the disease. In summary, our theoretical results are in accord with the investigation of Tongcheng Schistosomiasis Control Station and the report of Anhui Province Institute of Schistosomiasis for the whole Anhui province. Hence, we can obtain the result that after flood the delayed effect of flood may be caused by the incubation period of schistosome in the definitive host. Furthermore, the period of outbreak is about three years after flood in Anhui province.