# Stochastic Periodic Solution of a Susceptible-Infective Epidemic Model in a Polluted Environment under Environmental Fluctuation

Authors: Yu Zhao; Jiangping Li; Xu Ma
Journal: Computational and Mathematical Methods in Medicine
 (2018)

---

## Abstract

It is well known that the pollution and environmental fluctuations may seriously affect the outbreak of infectious diseases (e.g., measles). Therefore, understanding the association between the periodic outbreak of an infectious disease and noise and pollution still needs further development. Here we consider a stochastic susceptible-infective (SI) epidemic model in a polluted environment, which incorporates both environmental fluctuations as well as pollution. First, the existence of the global positive solution is discussed. Thereafter, the sufficient conditions for the nontrivial stochastic periodic solution and the boundary periodic solution of disease extinction are derived, respectively. Numerical simulation is also conducted in order to support the theoretical results. Our study shows that (i) large intensity noise may help the control of periodic outbreak of infectious disease; (ii) pollution may significantly affect the peak level of infective population and cause adverse health effects on the exposed population. These results can help increase the understanding of periodic outbreak patterns of infectious diseases.


---

## Body

## 1. Introduction

In Northern China, coal fire-power industries and heating systems, as well as vehicle emissions, all conduce to air pollution (airborne fine particulate matterPM2.5, PM10, and SO2, etc.), which has threatened the survival of exposed human population and affected the transmission of infectious diseases [1, 2]. Numerous studies have provided cumulative evidence of the health effects of particulate air pollution on the spread of infectious diseases (e.g., measles) [3, 4]. Therefore, investigating the role of pollution on the outbreak of infectious diseases is one of the most interesting and meaningful issues in the recent past [5].Dynamic mathematical models have provided a deeper understanding of the transmission process of infectious diseases [6, 7]. There are many interesting results (see, e.g., [8–10]), which show that the simple susceptible-infective (SI) model can fit the transmission process of some diseases (measles, chicken pox, etc.) well. Thus, based on the interaction between the environment and the population (see Figure 1), we incorporate the environmental pollution into the SI epidemic model. Let S(t), I(t), C0(t), and Ce(t) denote the number of people in the susceptible population, the number of people in the infective population, and the concentration of pollution in the organism and in the environment at time t, respectively. The SI epidemic model in a polluted environment is as follows: (1)dStdt=γSt1-St+ItK-βStIt-r10C0tSt-ξSt,dItdt=βStIt-r20C0tIt-vIt,dC0tdt=αCet-g+mC0t,dCetdt=ut-hCet;all of the parameters are positive constants and the corresponding biological meanings are listed in Table 1. Liu et al. [8] obtained the sufficient conditions of the ultimate boundedness of solutions and the global asymptotical stability of the equilibria.Table 1
Biological meanings of the parameters in model (1).
Parameters Biological meanings Unit γ The intrinsic growth rate in absence of the toxicant t - 1 K Carrying capacity of the population in absence of the toxicant Person β Probability of infection - ξ The population natural death rate t - 1 v The death rate caused by disease t - 1 r10 The dose-response rate due to uptake of pollution for the susceptible t - 1 r20 The dose-response rate due to uptake of pollution for the infected t - 1 α The organisms net uptake rate of pollution from environment t - 1 g The egestion rate of pollution in the organism (metabolism) t - 1 m The depuration rate of pollution in the organism t - 1 h The environmental pollution loss rate due to natural degradation t - 1 u(t) The exogenous rate of pollutant input into the environment t - 1Figure 1
Flow chart of the interaction between environmental pollution and population, and the disease transmits from the susceptible subpopulation to infected subpopulation.In real situations, as was pointed by Britton et al. [11], the transmission of infectious diseases is inevitably disrupted by unpredictable environmental conditions (e.g., absolute humidity [12], temperature [13]) making it more appropriate to use a stochastic model for biological parameters. For example, Yang et al. [14] observed the nonlinear effects of temperature and relative humidity on the incidence of measles. Thus, it is reasonable to model the environmental fluctuation as a stochastic transmission coefficient [15, 16]. If the transmission parameter β in the model (1) is subjected to some random environmental effects (temperature, humidity, etc.), then it is natural to consider that the transmission rate β is replaced by a random variable: (2)β~=β+σςt,where ς(t) is the Gaussian white noise with mean zero and variance one, ς(t)=dB(t)/dt, and B(t) is a scalar Wiener Process defined in (Ω,F,{Ft}t≥0,P), which is a complete probability space with a filtration {Ft}t≥0 satisfying the usual conditions (i.e., it is right continuous and F0 contains all P-null sets). σ is the intensity of the white noise. Thus, we can incorporate the environmental white noises into model (3): (3)dSt=Stγ1-St+ItK-βIt-r10C0t-ξdt-σStItdBt,dIt=ItβSt-r20C0t-vdt+σStItdBt,dC0tdt=αCet-g+mC0t,dCetdt=ut-hCet.Moreover, outbreaks of infectious diseases always fluctuate over time and exhibit seasonal patterns of incidence [17]. Ferrari et al. [18] pointed out that outbreaks of measles in the tropics have more variable seasonal patterns driven by accumulation and decline of susceptible individuals. Regular oscillatory patterns of measles outbreaks in Baltimore (USA) with an average period of three years have also been reported [19]. To describe the seasonal effect in the model, many existing studies [20–23] assume that the system parameters are subjected to a periodic rhythm. Therefore, we can further consider the following nonautonomous stochastic SI epidemic model as follows: (4)dSt=Stγt1-St+ItKt-βtIt-r10tC0t-ξtdt-σtStItdBt,dIt=ItβtSt-r20tC0t-vtdt+σtStItdBt,dC0tdt=αCet-g+mC0t,dCetdt=ut-hCet,with initial data (5)S0=S0≥0,I0=I0≥0,0≤C00≤1,0≤Ce0≤1,where γ(t), K(t), β(t), r10(t), r20(t), ξ(t), σ(t) are all positive, bounded, continuous θ-periodic functions.Considering the periodic variation and pollution exposure of epidemic models and exploring the existence of stochastic periodic solutions are meaningful to predict and control the outbreaks of infectious diseases. Such analysis has benefited from the theoretical contributions about the nonautonomous stochastic system [24, 25]. We also see that there has been some research in this respect [20, 21, 26, 27]. For example, Jiang et al. [21] considered a stochastic nonautonomous competitive Lotka-Volterra model in a polluted environment and then derived sufficient criteria for the existence and global attractivity of a nontrivial positive periodic solution. Xie et al. [27] presented a stochastic hepatitis B virus infection model with logistic hepatocyte growth and showed that the model has at least one periodic solution. More related results can be found in [22, 28, 29]. To the best of our knowledge, there are few results about the periodic solution of a stochastic SI epidemic model in a polluted environment. Therefore, the main objective of this paper is to concentrate on the effects of pollution and environmental fluctuation on the existence of the positive periodic solution.The rest of this paper is organized as follows: in the next Section2, we present the underlying mathematical analysis: the existence of the global positive solution, the sufficient conditions for the nontrivial stochastic periodic solution, and the boundary periodic solution of disease extinction are derived. The subsequent Section 3 describes the numerical simulation, based on the case of measles, carried out to support the theoretical results. Finally, in the last Section 4, the conclusion is presented.


## 2. Mathematical Analysis

### 2.1. Preliminary

SinceC0(t), Ce(t) are the concentrations of the pollution and 0≤C0(t)<1 and 0≤Ce(t)<1 must be satisfied, we assume the following.Assumption 1 (0<α≤g+m, u(t)≤h).
Notice thatu(t) is a positive θ-periodic continuous function, so we can prove the following.Lemma 2.
For model (4), we have limt→∞C0t-C0∗t=0; here (6)C0∗t=∫tt+θeg+ms-tαCesdseg+mθ-1.The proof of this lemma is provided in Appendix.From now on, we will only consider the following system:(7)dSt=Stγt1-St+ItKt-βtIt-r10tC0∗t-ξtdt-σtStItdBt,dIt=ItβtSt-r20tC0∗t-vtdt+σtStItdBt.For a bounded function on0,∞, say, f(t), define (8)fu=supt∈0,∞⁡ft,fl=inft∈0,∞⁡ft,fθ=1θ∫0θfsds.Now, we shall give some definitions and Lemmas with respect to the periodic Markov process X(t) as the solution of stochastic system (9)Xt=Xt0+∫t0tbs,Xsds+∑r=1k∫t0tσrs,XsdBrs,x∈Rl.Definition 3 (see [24]).
A stochastic processξ(t)=ξ(t,ω)(-∞<t<∞) is said to be periodic with period θ, if for every finite sequence of numbers t1,t2,…,tn the joint distribution of the random variables ξ(t1+h),…,ξ(tn+h) is independent of h, where h=kθ, k=±1,±2,….Remark 4.
It follows from [24] that a stochastic Markov process X(t) is θ-periodic if and only if its transition probability function is θ-periodic and the function P0(t,A)=PXt∈A satisfies(10)P0s,A=∫RlP0s,dxPs,x,s+θ,A≡P0s+θ,A,for every A∈B, where B denotes the Borel σ-algebra in Rl.LetL be a linear operator defined by (11)L=∂∂t+∑i=1nbit,x∂∂xi+12∑i,j=1naij∂2∂xi∂xj,aij=∑r=1kσrit,xσrjt,x.Lemma 5 (see [24]).
Suppose that the coefficients of system (13) are θ-periodic in t and satisfy (12)bs,x-bs,y+∑r=1kσrs,x-σrs,x≤Cx-y,bs,x+∑r=1kσrs,x≤C1+x,in every cylinder I×U, where C is a constant; and suppose further that there exists a function V(t,x)∈Rl which is θ-periodic in t, satisfying (13)infx>R⁡Vt,x⟶∞,asR⟶∞,(14)LVt,x≤-1,outside  some  compact  set.Then, there exists a solution of system (9) which is a θ-periodic Markov process.Remark 6.
According to the proof of Lemma 2.1 in [24], condition (14) is only used to guarantee the existence and uniqueness of the solution of system (9).Next, we have the following theorem.Theorem 7.
For any given initial value (5), model (7) has a unique positive solution (S(t),I(t)) on t≥0, and the solution will remain in R+2 with probability one.The proof of this Theorem is provided in Appendix.


### 2.2. Existence of the Positive Stochastic Periodic Solution

In this section, we shall prove the existence of a positive stochastic periodic solution of models (4) and (7). Firstly, we define (15)λθ=1θ∫0θγs-r10s+r20sC0∗s-12γuKuγl2σ2s-ξs-vsds.Now, we obtain the following result regarding the existence of a positive periodic solution of model (7).Theorem 8.
Ifλθ>0, then model (7) at least has one positive θ-periodic solution (S∗(t),I∗(t)).Proof.
To prove the existence of a positiveθ-periodic solution of model (7), it follows from Lemma 2 and Remark 6 that we need to find a C2-function V(t,S,I) and a closed set Θ∈R+2 such that (13) and (14) hold. Firstly, we assume that γu+γu/Kl-βl-(r10C0∗+ξ)l>0 and define a nonnegative function as follows: (16)Vt,S,I=V1t,S+V2t,S,I+V3t,where (17)
V
1
t
,
S
=
-
M
ln
⁡
S
t
-
M
ln
⁡
I
t
,
V
2
t
,
S
,
I
=
1
2
S
t
+
I
t
2
+
M
S
t
+
I
t
,
V
3
t
=
M
ω
¯
t
,
M
=
2
λ
θ
max
⁡
1
,
sup
S
,
I
∈
R
+
2
⁡
-
γ
l
K
u
S
3
+
γ
u
-
r
10
C
0
∗
+
ξ
l
-
M
γ
l
K
u
+
σ
u
2
I
2
S
2
+
γ
u
-
r
20
C
0
∗
+
v
l
-
r
10
C
0
∗
+
ξ
l
-
M
γ
l
K
u
S
I
-
r
20
C
0
∗
+
v
l
I
2
+
M
γ
u
K
l
+
γ
u
-
β
l
-
r
10
C
0
∗
+
ξ
l
S
+
M
γ
l
K
u
S
3
+
γ
u
-
r
10
C
0
∗
+
ξ
l
-
M
γ
l
K
u
+
σ
u
2
I
2
γ
u
K
l
+
β
u
-
r
20
C
0
∗
+
v
l
I
. It is easy to see that Mλθ/2≥1 and ω¯(t) satisfies the following: (18)ω¯˙t=-λθ+γt-r10t+r20tC0∗t-12γuKuγl2σ2t-ξt-vt,ω¯0=0.
Integrating (18) from t to t+θ yields (19)ω¯t+θ-ω¯t=∫tt+θω¯˙sds=∫tt+θ-λθ+γs-r10s+r20sC0∗s-12γuKuγl2σ2s-ξs-vsds-∫0θ-λθ+γs-r10s+r20sC0∗s-12γuKuγl2σ2s-ξs-vsds=0.Thus, we can see that ω¯(t) is a θ-periodic function on 0,∞ and (20)liminfϵ→0,S,I∈R+2∖Θϵ⁡Vt,S,I=∞,where Θϵ={(S,I):(S,I)∈(ϵ,1/ϵ)×(ϵ,1/ϵ)}. Thus, V(t,S,I) is θ-periodic with respect to t.
Now, we have the requisite information to verify (14) in Lemma 5. Applying the Itô formula to V(t,S,I), one can obtain (21)LV1=-Mγt-γtktS+I-r10tC0∗t-βtI-ξt-σ2tI22+βtS-r20tC0∗t-vt-σ2tS22≤M-γt-r10t+r20tC0∗t-ξt-vt+γuKlS+γuKl+βuI+γuKuγl2σ2t-βlS=M-γt-r10t+r20tC0∗t-12γuKuγl2σ2t-ξt-vt+γuKlS+γuKl+βuI-βlS,LV2=S+IγtS1-S+IKt-r10tC0∗t+ξtS-r20tC0∗t+vtI+MγtS1-S+IKt-r10tC0∗t+ξtS-r20tC0∗t+vtI+σ2tS2I2≤γuS2-γlKuS3-γlKuS2I-r10C0∗+ξlS2-r20C0∗+vlSI-γlKuS2I-γlKuSI2-r10C0∗+ξlSI-r20C0∗+vlI2+γuSI+σu2S2I2+MγuS-MγlKuS2-MγlKuIS-Mr10C0∗+ξlS-Mr20C0∗+vlI≤-γlKuS3+γu-r10C0∗+ξl-MγlKu+σu2I2S2-r20C0∗+vlI2+γu-r20C0∗+vl-r10C0∗+ξl-MγlKuSI+Mγu-r10C0∗+ξlS-Mr20C0∗+vlI.From (18) and (21), we can get that (22)LV=-Mλθ+MγuKl-βlS+MγuKl+βuI-γlKuS3+γu-r10C0∗+ξl-MγlKu+σu2I2S2-r20C0∗+vlI2+γu-r20C0∗+vl-r10C0∗+ξl-MγlKuSI+Mγu-r10C0∗+ξlS-Mr20C0∗+vlI.Define a closed set (23)Θϵ=S,I∈R+2:ϵ≤S≤1ϵ,ϵ≤I≤1ϵ,where 0<ϵ<1 is a sufficiently small number such that (24)-Mλθ+MγuKl+γu-βl-r10C0∗+ξlϵ+C1≤-1,(25)-Mλθ-Mβl+r10C0∗+ξl1ϵ+C2≤-1,(26)-Mλθ-r20C0∗+vl1ϵ+C3≤-1,where Ci, i=1,2,3 are positive constants defined in (29), (32), and (34) later. Moreover, we denote (27)Θϵ1=S,I∈R+2:0<S<ϵ,Θϵ2=S,I∈R+2:0<I<ϵ,Θϵ3=S,I∈R+2:S>1ϵ,Θϵ4=S,I∈R+2:I>1ϵ.Then Θϵc=R+2∖Θϵ=Θϵ1∪Θϵ2∪Θϵ3∪Θϵ4. Next, we shall prove LV(t,S,I)≤-1 on [0,∞]×Θϵc.
Case  1.On Θϵ1, we have 0<S<ϵ: (28)LV=-Mλθ+MγuKl+γu-βl-r10C0∗+ξlS-γlKuS3+γu-r10C0∗+ξl-MγlKu+σu2I2S2-r20C0∗+vlI2+γu-r20C0∗+vl-r10C0∗+ξl-MγlKuSI+MγuKl+βu-r20C0∗+vlI≤-Mλθ+MγuKl-βl+γu-r10C0∗+ξlϵ+C1,where (29)C1=supS,I∈R+2⁡-γlKuS3+γu-r10C0∗+ξl-MγlKu+σu2I2S2-r20C0∗+vlI2+γu-r20C0∗+vl-r10C0∗+ξl-MγlKuSI-Mr20C0∗+vl-γuKl-βuI.Therefore, we can say that LV(t,S,I)≤-1 on [0,∞]×Θϵ1 in lieu of (24).
Case  2.On Θϵ2, we have 0<I<ϵ: (30)LV≤-Mλθ2+-Mλθ2+supS,I∈R+2-γlKuS3+γu-r10C0∗+ξl-MγlKu+σu2I2S2-r20C0∗+vlI2+γu-r20C0∗+vl-r10C0∗+ξl-MγlKuSI+MγuKl+γu-βl-r10C0∗+ξlS+MγuKl+βu-r20C0∗+vlI.It follows from the definition of M that LV≤-Mλθ/2≤-1, which implies LV(t,S,I)≤-1 on [0,∞]×Θϵ2.
Case  3.On Θϵ3, we have S>1/ϵ: (31)LV=-Mλθ-Mβl+r10C0∗+ξlS+-γlKuS3+γu-r10C0∗+ξl-MγlKu+σu2I2S2-r20C0∗+vlI2+γu-r20C0∗+vl-r10C0∗+ξl-MγlKuSI+MγuKl+βu-r20C0∗+vlI+MγuKl+γuS≤-Mλθ-Mβl+r10C0∗+ξl1ϵ+C2,where (32)C2=supS,I∈R+2⁡-γlKuS3+γu-r10C0∗+ξl-MγlKu+σu2I2S2-r20C0∗+vlI2+γu-r20C0∗+vl-r10C0∗+ξl-MγlKuSI+MγuKl+βu-r20C0∗+vlI+MγuKl+γuS.According to (25), one can get that LV(t,S,I)≤-1 on [0,∞]×Θϵ3.
Case  4.On Θϵ4, we have I>1/ϵ: (33)LV=-Mλθ-r20C0∗+vlI+-γlKuS3+γu-r10C0∗+ξl-MγlKu+σu2I2S2-r20C0∗+vlI2+γu-r20C0∗+vl-r10C0∗+ξl-MγlKuSI+MγuKl+βuI+MγuKl+γu-βl-r10C0∗+ξlS≤-Mλθ-r20C0∗+vl1ϵ+C3,where (34)C3=supS,I∈R+2⁡-γlKuS3+γu-r10C0∗+ξl-MγlKu+σu2I2S2-r20C0∗+vlI2+γu-r20C0∗+vl-r10C0∗+ξl-MγlKuSI+MγuKl+βuI+MγuKl+γu-βl-r10C0∗+ξlS.In view of (26), we can deduce that LV(t,S,I)≤-1 on [0,∞]×Θϵ4.
In summary, we can draw the conclusion that(35)LVt,S,I≤-1for  every0,∞×Θϵc.Thus, the condition (14) of Lemma 5 is satisfied. Consequently, model (7) at least has one positive stochastic θ-periodic solution.Remark 9.
Theorem8 implies that the intrinsic growth rate of population should overcome the extinction risks of infected disease and pollution in order to guarantee the survival of the population. In addition, the condition λθ means that the susceptible population evolution dynamic (γ,K), the dose-response rates (ri0, i=1,2), and intensity of noise (σ) play an important role in determining the periodic outbreak of infectious disease; that is, reducing the possibility of λθ>1 is beneficial to the control of the periodic outbreak of infectious diseases.Remark 10.
In contrast to the authors in [8] assume that the exogenous rate of pollutant input u(t) has constant limit(i.e., limt→∞u(t)=u∗) and derived the stability of equilibria, which ignored the periodicity of model. In this paper, we consider the limit periodic system with attractiveness (i.e., limt→∞C0t-C0∗t=0) and obtain the positive stochastic periodic solution, which extends the results in [8] to a stochastic nonautonomous situation.Combining Lemma2 and Theorem 8, we can have the following result on the existence of positive stochastic periodic solution of model (4).Theorem 11.
Under the conditions of Theorem8, model (4) at least has one positive stochastic θ-periodic solution (S∗(t),I∗(t),C0∗(t),Ce∗(t)).


### 2.3. The Boundary Periodic Solution of Disease Extinction

In this section, we shall obtain the sufficient conditions for disease extinction. Firstly, we define(36)λθ0=1θ∫0θβ2s2σ2s-r20sC0∗s-vsds.Next, we prove the following theorem.Theorem 12.
For model (7), if λθ0<0, then the disease I(t) goes to extinct almost surely.Proof.
DefineW(I)=ln⁡It for I∈[0,γuKu/γl]. Utilizing the Itô formula to model (7) yields (37)dWIt=LWItdt+σtStdBt,where (38)LWIt=βtSt-r20tC0∗t-vt-σ2t2S2t≤β2t2σ2t-r20tC0∗t-vt.Substituting (38) into (37) and integrating both sides of (37), one can deduce that (39)1t∫0tln⁡ItI0≤1t∫0tβ2s2σ2s-r20sC0∗s-vsds+Mtt,where M(t)=∫0tσ(s)S(s)dB(s) is a martingale with the following quadratic variation: (40)Mt,Mtt=∫0tσuSu2du≤sups≥0⁡σuγuKuγl2,a.s.According to the strong law of large numbers for martingales [32], we can get (41)limt→∞Mtt=0,a.s.Combining (39), (41), and the coefficients’ periodicity of model (7), we get (42)limt→∞1t∫0tln⁡ItI0≤limt→∞⁡1t∫0tβ2s2σ2s-r20sC0∗s-vsds+limt→∞⁡Mtt=1θ∫0θβ2s2σ2s-r20sC0∗s-vsds=λθ0<0,and hence limt→∞I(t)=0, a.s.Note the fact that whenlimt→∞I(t)=0, model (7) reduces to the following nonautonomous system: (43)dSt=Stat-btSt,where a(t)=γ(t)-r10(t)C0∗(t)-ξ(t), b(t)=γ(t)/K(t) are all θ-periodic functions. Define (44)λθS=1θ∫0θγs-r10sC0∗s-ξsds.Therefore, we have the following θ-periodic solution result of system (43).Lemma 13 (see Globalism [33]).
For system (43), if λθS>0 then it has a stable positive θ-periodic solution Sθ0(t) which satisfies (45)1Sθ0t=∫tt+θexp⁡∫tsaτdτbsdsexp⁡∫0θaτdτ-1,t≥0.In summary, we obtain the following.Theorem 14.
For model (7), if λθ0<0 and λθS>0, then it has a boundary periodic solution of disease extinction (Sθ0(t),0).Remark 15.
The transmission coefficientβ, intensity of noise σ, and pollution level r20C0∗ may determine the fate of the evolution of infected population.


## 2.1. Preliminary

SinceC0(t), Ce(t) are the concentrations of the pollution and 0≤C0(t)<1 and 0≤Ce(t)<1 must be satisfied, we assume the following.Assumption 1 (0<α≤g+m, u(t)≤h).
Notice thatu(t) is a positive θ-periodic continuous function, so we can prove the following.Lemma 2.
For model (4), we have limt→∞C0t-C0∗t=0; here (6)C0∗t=∫tt+θeg+ms-tαCesdseg+mθ-1.The proof of this lemma is provided in Appendix.From now on, we will only consider the following system:(7)dSt=Stγt1-St+ItKt-βtIt-r10tC0∗t-ξtdt-σtStItdBt,dIt=ItβtSt-r20tC0∗t-vtdt+σtStItdBt.For a bounded function on0,∞, say, f(t), define (8)fu=supt∈0,∞⁡ft,fl=inft∈0,∞⁡ft,fθ=1θ∫0θfsds.Now, we shall give some definitions and Lemmas with respect to the periodic Markov process X(t) as the solution of stochastic system (9)Xt=Xt0+∫t0tbs,Xsds+∑r=1k∫t0tσrs,XsdBrs,x∈Rl.Definition 3 (see [24]).
A stochastic processξ(t)=ξ(t,ω)(-∞<t<∞) is said to be periodic with period θ, if for every finite sequence of numbers t1,t2,…,tn the joint distribution of the random variables ξ(t1+h),…,ξ(tn+h) is independent of h, where h=kθ, k=±1,±2,….Remark 4.
It follows from [24] that a stochastic Markov process X(t) is θ-periodic if and only if its transition probability function is θ-periodic and the function P0(t,A)=PXt∈A satisfies(10)P0s,A=∫RlP0s,dxPs,x,s+θ,A≡P0s+θ,A,for every A∈B, where B denotes the Borel σ-algebra in Rl.LetL be a linear operator defined by (11)L=∂∂t+∑i=1nbit,x∂∂xi+12∑i,j=1naij∂2∂xi∂xj,aij=∑r=1kσrit,xσrjt,x.Lemma 5 (see [24]).
Suppose that the coefficients of system (13) are θ-periodic in t and satisfy (12)bs,x-bs,y+∑r=1kσrs,x-σrs,x≤Cx-y,bs,x+∑r=1kσrs,x≤C1+x,in every cylinder I×U, where C is a constant; and suppose further that there exists a function V(t,x)∈Rl which is θ-periodic in t, satisfying (13)infx>R⁡Vt,x⟶∞,asR⟶∞,(14)LVt,x≤-1,outside  some  compact  set.Then, there exists a solution of system (9) which is a θ-periodic Markov process.Remark 6.
According to the proof of Lemma 2.1 in [24], condition (14) is only used to guarantee the existence and uniqueness of the solution of system (9).Next, we have the following theorem.Theorem 7.
For any given initial value (5), model (7) has a unique positive solution (S(t),I(t)) on t≥0, and the solution will remain in R+2 with probability one.The proof of this Theorem is provided in Appendix.


## 2.2. Existence of the Positive Stochastic Periodic Solution

In this section, we shall prove the existence of a positive stochastic periodic solution of models (4) and (7). Firstly, we define (15)λθ=1θ∫0θγs-r10s+r20sC0∗s-12γuKuγl2σ2s-ξs-vsds.Now, we obtain the following result regarding the existence of a positive periodic solution of model (7).Theorem 8.
Ifλθ>0, then model (7) at least has one positive θ-periodic solution (S∗(t),I∗(t)).Proof.
To prove the existence of a positiveθ-periodic solution of model (7), it follows from Lemma 2 and Remark 6 that we need to find a C2-function V(t,S,I) and a closed set Θ∈R+2 such that (13) and (14) hold. Firstly, we assume that γu+γu/Kl-βl-(r10C0∗+ξ)l>0 and define a nonnegative function as follows: (16)Vt,S,I=V1t,S+V2t,S,I+V3t,where (17)
V
1
t
,
S
=
-
M
ln
⁡
S
t
-
M
ln
⁡
I
t
,
V
2
t
,
S
,
I
=
1
2
S
t
+
I
t
2
+
M
S
t
+
I
t
,
V
3
t
=
M
ω
¯
t
,
M
=
2
λ
θ
max
⁡
1
,
sup
S
,
I
∈
R
+
2
⁡
-
γ
l
K
u
S
3
+
γ
u
-
r
10
C
0
∗
+
ξ
l
-
M
γ
l
K
u
+
σ
u
2
I
2
S
2
+
γ
u
-
r
20
C
0
∗
+
v
l
-
r
10
C
0
∗
+
ξ
l
-
M
γ
l
K
u
S
I
-
r
20
C
0
∗
+
v
l
I
2
+
M
γ
u
K
l
+
γ
u
-
β
l
-
r
10
C
0
∗
+
ξ
l
S
+
M
γ
l
K
u
S
3
+
γ
u
-
r
10
C
0
∗
+
ξ
l
-
M
γ
l
K
u
+
σ
u
2
I
2
γ
u
K
l
+
β
u
-
r
20
C
0
∗
+
v
l
I
. It is easy to see that Mλθ/2≥1 and ω¯(t) satisfies the following: (18)ω¯˙t=-λθ+γt-r10t+r20tC0∗t-12γuKuγl2σ2t-ξt-vt,ω¯0=0.
Integrating (18) from t to t+θ yields (19)ω¯t+θ-ω¯t=∫tt+θω¯˙sds=∫tt+θ-λθ+γs-r10s+r20sC0∗s-12γuKuγl2σ2s-ξs-vsds-∫0θ-λθ+γs-r10s+r20sC0∗s-12γuKuγl2σ2s-ξs-vsds=0.Thus, we can see that ω¯(t) is a θ-periodic function on 0,∞ and (20)liminfϵ→0,S,I∈R+2∖Θϵ⁡Vt,S,I=∞,where Θϵ={(S,I):(S,I)∈(ϵ,1/ϵ)×(ϵ,1/ϵ)}. Thus, V(t,S,I) is θ-periodic with respect to t.
Now, we have the requisite information to verify (14) in Lemma 5. Applying the Itô formula to V(t,S,I), one can obtain (21)LV1=-Mγt-γtktS+I-r10tC0∗t-βtI-ξt-σ2tI22+βtS-r20tC0∗t-vt-σ2tS22≤M-γt-r10t+r20tC0∗t-ξt-vt+γuKlS+γuKl+βuI+γuKuγl2σ2t-βlS=M-γt-r10t+r20tC0∗t-12γuKuγl2σ2t-ξt-vt+γuKlS+γuKl+βuI-βlS,LV2=S+IγtS1-S+IKt-r10tC0∗t+ξtS-r20tC0∗t+vtI+MγtS1-S+IKt-r10tC0∗t+ξtS-r20tC0∗t+vtI+σ2tS2I2≤γuS2-γlKuS3-γlKuS2I-r10C0∗+ξlS2-r20C0∗+vlSI-γlKuS2I-γlKuSI2-r10C0∗+ξlSI-r20C0∗+vlI2+γuSI+σu2S2I2+MγuS-MγlKuS2-MγlKuIS-Mr10C0∗+ξlS-Mr20C0∗+vlI≤-γlKuS3+γu-r10C0∗+ξl-MγlKu+σu2I2S2-r20C0∗+vlI2+γu-r20C0∗+vl-r10C0∗+ξl-MγlKuSI+Mγu-r10C0∗+ξlS-Mr20C0∗+vlI.From (18) and (21), we can get that (22)LV=-Mλθ+MγuKl-βlS+MγuKl+βuI-γlKuS3+γu-r10C0∗+ξl-MγlKu+σu2I2S2-r20C0∗+vlI2+γu-r20C0∗+vl-r10C0∗+ξl-MγlKuSI+Mγu-r10C0∗+ξlS-Mr20C0∗+vlI.Define a closed set (23)Θϵ=S,I∈R+2:ϵ≤S≤1ϵ,ϵ≤I≤1ϵ,where 0<ϵ<1 is a sufficiently small number such that (24)-Mλθ+MγuKl+γu-βl-r10C0∗+ξlϵ+C1≤-1,(25)-Mλθ-Mβl+r10C0∗+ξl1ϵ+C2≤-1,(26)-Mλθ-r20C0∗+vl1ϵ+C3≤-1,where Ci, i=1,2,3 are positive constants defined in (29), (32), and (34) later. Moreover, we denote (27)Θϵ1=S,I∈R+2:0<S<ϵ,Θϵ2=S,I∈R+2:0<I<ϵ,Θϵ3=S,I∈R+2:S>1ϵ,Θϵ4=S,I∈R+2:I>1ϵ.Then Θϵc=R+2∖Θϵ=Θϵ1∪Θϵ2∪Θϵ3∪Θϵ4. Next, we shall prove LV(t,S,I)≤-1 on [0,∞]×Θϵc.
Case  1.On Θϵ1, we have 0<S<ϵ: (28)LV=-Mλθ+MγuKl+γu-βl-r10C0∗+ξlS-γlKuS3+γu-r10C0∗+ξl-MγlKu+σu2I2S2-r20C0∗+vlI2+γu-r20C0∗+vl-r10C0∗+ξl-MγlKuSI+MγuKl+βu-r20C0∗+vlI≤-Mλθ+MγuKl-βl+γu-r10C0∗+ξlϵ+C1,where (29)C1=supS,I∈R+2⁡-γlKuS3+γu-r10C0∗+ξl-MγlKu+σu2I2S2-r20C0∗+vlI2+γu-r20C0∗+vl-r10C0∗+ξl-MγlKuSI-Mr20C0∗+vl-γuKl-βuI.Therefore, we can say that LV(t,S,I)≤-1 on [0,∞]×Θϵ1 in lieu of (24).
Case  2.On Θϵ2, we have 0<I<ϵ: (30)LV≤-Mλθ2+-Mλθ2+supS,I∈R+2-γlKuS3+γu-r10C0∗+ξl-MγlKu+σu2I2S2-r20C0∗+vlI2+γu-r20C0∗+vl-r10C0∗+ξl-MγlKuSI+MγuKl+γu-βl-r10C0∗+ξlS+MγuKl+βu-r20C0∗+vlI.It follows from the definition of M that LV≤-Mλθ/2≤-1, which implies LV(t,S,I)≤-1 on [0,∞]×Θϵ2.
Case  3.On Θϵ3, we have S>1/ϵ: (31)LV=-Mλθ-Mβl+r10C0∗+ξlS+-γlKuS3+γu-r10C0∗+ξl-MγlKu+σu2I2S2-r20C0∗+vlI2+γu-r20C0∗+vl-r10C0∗+ξl-MγlKuSI+MγuKl+βu-r20C0∗+vlI+MγuKl+γuS≤-Mλθ-Mβl+r10C0∗+ξl1ϵ+C2,where (32)C2=supS,I∈R+2⁡-γlKuS3+γu-r10C0∗+ξl-MγlKu+σu2I2S2-r20C0∗+vlI2+γu-r20C0∗+vl-r10C0∗+ξl-MγlKuSI+MγuKl+βu-r20C0∗+vlI+MγuKl+γuS.According to (25), one can get that LV(t,S,I)≤-1 on [0,∞]×Θϵ3.
Case  4.On Θϵ4, we have I>1/ϵ: (33)LV=-Mλθ-r20C0∗+vlI+-γlKuS3+γu-r10C0∗+ξl-MγlKu+σu2I2S2-r20C0∗+vlI2+γu-r20C0∗+vl-r10C0∗+ξl-MγlKuSI+MγuKl+βuI+MγuKl+γu-βl-r10C0∗+ξlS≤-Mλθ-r20C0∗+vl1ϵ+C3,where (34)C3=supS,I∈R+2⁡-γlKuS3+γu-r10C0∗+ξl-MγlKu+σu2I2S2-r20C0∗+vlI2+γu-r20C0∗+vl-r10C0∗+ξl-MγlKuSI+MγuKl+βuI+MγuKl+γu-βl-r10C0∗+ξlS.In view of (26), we can deduce that LV(t,S,I)≤-1 on [0,∞]×Θϵ4.
In summary, we can draw the conclusion that(35)LVt,S,I≤-1for  every0,∞×Θϵc.Thus, the condition (14) of Lemma 5 is satisfied. Consequently, model (7) at least has one positive stochastic θ-periodic solution.Remark 9.
Theorem8 implies that the intrinsic growth rate of population should overcome the extinction risks of infected disease and pollution in order to guarantee the survival of the population. In addition, the condition λθ means that the susceptible population evolution dynamic (γ,K), the dose-response rates (ri0, i=1,2), and intensity of noise (σ) play an important role in determining the periodic outbreak of infectious disease; that is, reducing the possibility of λθ>1 is beneficial to the control of the periodic outbreak of infectious diseases.Remark 10.
In contrast to the authors in [8] assume that the exogenous rate of pollutant input u(t) has constant limit(i.e., limt→∞u(t)=u∗) and derived the stability of equilibria, which ignored the periodicity of model. In this paper, we consider the limit periodic system with attractiveness (i.e., limt→∞C0t-C0∗t=0) and obtain the positive stochastic periodic solution, which extends the results in [8] to a stochastic nonautonomous situation.Combining Lemma2 and Theorem 8, we can have the following result on the existence of positive stochastic periodic solution of model (4).Theorem 11.
Under the conditions of Theorem8, model (4) at least has one positive stochastic θ-periodic solution (S∗(t),I∗(t),C0∗(t),Ce∗(t)).


## 2.3. The Boundary Periodic Solution of Disease Extinction

In this section, we shall obtain the sufficient conditions for disease extinction. Firstly, we define(36)λθ0=1θ∫0θβ2s2σ2s-r20sC0∗s-vsds.Next, we prove the following theorem.Theorem 12.
For model (7), if λθ0<0, then the disease I(t) goes to extinct almost surely.Proof.
DefineW(I)=ln⁡It for I∈[0,γuKu/γl]. Utilizing the Itô formula to model (7) yields (37)dWIt=LWItdt+σtStdBt,where (38)LWIt=βtSt-r20tC0∗t-vt-σ2t2S2t≤β2t2σ2t-r20tC0∗t-vt.Substituting (38) into (37) and integrating both sides of (37), one can deduce that (39)1t∫0tln⁡ItI0≤1t∫0tβ2s2σ2s-r20sC0∗s-vsds+Mtt,where M(t)=∫0tσ(s)S(s)dB(s) is a martingale with the following quadratic variation: (40)Mt,Mtt=∫0tσuSu2du≤sups≥0⁡σuγuKuγl2,a.s.According to the strong law of large numbers for martingales [32], we can get (41)limt→∞Mtt=0,a.s.Combining (39), (41), and the coefficients’ periodicity of model (7), we get (42)limt→∞1t∫0tln⁡ItI0≤limt→∞⁡1t∫0tβ2s2σ2s-r20sC0∗s-vsds+limt→∞⁡Mtt=1θ∫0θβ2s2σ2s-r20sC0∗s-vsds=λθ0<0,and hence limt→∞I(t)=0, a.s.Note the fact that whenlimt→∞I(t)=0, model (7) reduces to the following nonautonomous system: (43)dSt=Stat-btSt,where a(t)=γ(t)-r10(t)C0∗(t)-ξ(t), b(t)=γ(t)/K(t) are all θ-periodic functions. Define (44)λθS=1θ∫0θγs-r10sC0∗s-ξsds.Therefore, we have the following θ-periodic solution result of system (43).Lemma 13 (see Globalism [33]).
For system (43), if λθS>0 then it has a stable positive θ-periodic solution Sθ0(t) which satisfies (45)1Sθ0t=∫tt+θexp⁡∫tsaτdτbsdsexp⁡∫0θaτdτ-1,t≥0.In summary, we obtain the following.Theorem 14.
For model (7), if λθ0<0 and λθS>0, then it has a boundary periodic solution of disease extinction (Sθ0(t),0).Remark 15.
The transmission coefficientβ, intensity of noise σ, and pollution level r20C0∗ may determine the fate of the evolution of infected population.


## 3. Numerical Simulation

In this section, we shall verify the above theoretical results and illustrate the effects of environmental fluctuation and pollution on the periodic outbreak of infectious disease. With help from MATLAB (Mathworks, Inc., Natick, MA, USA) and Milsteins higher order method [34], which is a powerful tool for solving stochastic differential equations, we consider the following discretized equation of model (7) at t=(k+1)Δt, k=0,1,…: (46)Sk+1=Sk+SkγkΔt1-Sk+IkKkΔt-βkΔtIk-r10kΔtC0∗kΔt-ξkΔtΔt-SkIkσkΔtkΔtξk+σ2kΔt2ξk2Δt-Δt,Ik+1=Ik+IkβkΔtSk-r20kΔtC0∗kΔt-vkΔtΔt+SkIkσkΔtkΔtξk+σ2kΔt2ξk2Δt-Δt,where ξk are the N(0,1)-distribution independent Gaussian random variables. Let us assume that (47)γt=3.5+0.5cos⁡t12,Kt=100+5cos⁡t12,ξt=0.34+0.05cos⁡t12,vt=0.05+0.01cos⁡t12,C0∗t=0.5+0.05cos⁡t12,σt=0.01+0.005cos⁡t12,βt=0.25+0.05cos⁡t12.Example 16.
To illustrate the effect of pollution on the periodic outbreaks of infectious disease, we look at the following two cases that differ with respect to the average pollution level.
Case (i). Assume that r10(t)=0.15+0.05cos(t/12),r20(t)=0.2+0.05cos(t/12). After some simple calculations, we can see that γu+γu/Kl-βl-(r10C0∗+ξ)l=1.8711>0 and (48)λθ=3.3871>0.Thus, it follows from Theorem8 that model (7) has at least one positive 24π-periodic solution. As shown in Figures 2(a) and 2(c), the probability density functions (PDFs) of (S(t),I(t)) of model (7) are nearly equal (from the shape of stationary distribution aspect) to each other in different periods, which supports the definition of a periodic Markov process (see Figures 2(b) and 2(d)). Therefore, the solution process of model (7) is a 24π-periodic Markov process. Additionally, in the absence of white noise, model (7) reduces to a deterministic system. Hence, we also plot the trajectories of the corresponding deterministic model (7) in Figures 2(a) and 2(c) (red lines). In summary, it can be observed that the sample trajectories of the stochastic model (7) have regular periodicity under small environmental fluctuations, and fluctuation happens around the periodic solution of the corresponding deterministic counterparts.Figure 2
(a), (c) The sample trajectory ofS(t) and I(t) of model (7) (blue lines) and their corresponding deterministic periodic solution (red lines), respectively. (b), (d) the probability density functions (PDFs) of S(t) and I(t) of model (7) in different periods, respectively. The initial value is (1,11) and the parameter values are used as in Example 16.
(a) (b) (c) (d)Case (ii). Assume that r10(t)=0.1+0.05cost/12, r20(t)=0.15+0.05cos(t/12). After some simple calculations, we can check that γu+γu/Kl-βl-(r10C0∗+ξ)l=1.9711>0 and λθ=3.4121>0. According to Theorem 8, model (7) has a stochastic 24π-periodic solution (see Figure 3(a), red lines). We can also observe from Figure 3(a) that due to the decrease in pollution, the peak level of infective population increases. Meanwhile, the corresponding PDF moves to the right position, which implies a higher number for population of I(t) (see Figure 2(b)). Thus, pollution may increase the peak level of infective population.Figure 3
(a) The sample trajectoryI(t) of model (7) corresponding to different average pollution level ri0C0∗, respectively. (b) the PDFs of I(t) of model (7) in the same period with different average pollution level, respectively.
(a) (b)Next, we shall check the existence of a boundary periodic solution.Example 17.
Let us assume thatr10(t)=0.15+0.05cost/12, r20(t)=0.2+0.05cost/12, β(t)=0.0015+0.0005cos(t/12); then we can calculate that (49)λθ0=-0.2388<0,λθS=3.011>0.It follows from Theorem 12 that there exists a boundary periodic solution of disease extinction (Sθ0(t),0) of model (7), which is consistent with the simulation results as shown in Figure 4.Figure 4
The sample trajectory of the boundary periodic solution(Sθ0(t),0) of model (7), with initial values (75,10).
(a) (b)Now, we are in a position to see the fit of model (7) for a real-world situation (the case of measles).Example 18.
Measles is a highly contagious airborne infectious disease caused by the measles virus, which spreads easily through coughing and sneezing of infected people. Major epidemics occur approximately every 2-3 years, causing an estimated 2.6 million deaths each year [31]. Recent research has showed that the incidence of measles is related to air pollution in China [1, 3, 35] and has provided cumulative evidence of the adverse health effects of particulate air pollution and dust. Thus, using measles as an example, we have the requisite information to fit a real-world situation such as the outbreak of measles by using model (7). The data source for the cases of measles was from the Chinese center for disease control and prevention (CCDCP) [36]. The parameters of the simulation are listed in Table 2.Table 2
The values of the parameters in model (46).
Parameters Mean value Source γ 0.01295 [30] β 0.00495 Estimated ξ 0.0067 [30] σ 0.05 [30] K 1.378 × 1 0 9 [30] v 0.00175 [31] r i 0, i=1,2 0.0015 Estimated C 0 ∗ 1.05 [1]It can be seen from Figure5 that our simulation, based on model (7), is a good fit compared to the data on the observed cases of measles in the period from Jan 2014 to Dec 2016. However, the results of fitting the model in the period from Jan 2017 to Dec 2017 are not good; this may be due to the following reasons: (1) by 2016, the government and prevention departments push to improve vaccine coverage may have resulted in the drop in the cases of measles [31]; (2) the Chinese government paid more attention to control the pollution in the environment and enforced strict emission standards, which reduced the negative effect on the population health [35]. Moreover, one interesting finding is that the decrease of susceptible S(t) is beneficial to the decreasing trend of infected cases I(t); that is, if we can reduce the number of the susceptible subpopulation, the cases of measles will show a decreasing trend in the future (see Figure 5, blue lines). Thus, strengthening the coverage of measles vaccination and the environmental quality improvement are advantageous in controlling the outbreak of infectious diseases (e.g., measles).Figure 5
Observed (red double cross pattern) of cases of measlesI(t) during Jan 2013 to Dec 2016 and model fitted and predicted cases of measles (blue line) from Jan 2013 to Dec 2019. The inner panel is the sample path of S(t).


## 4. Conclusion

Generally, humans are exposed to some kinds of infectious diseases because the diseases propagate through a polluted environment. Examples include measles spreading through air pollution, snail fever spreading through water pollution, and diarrhea spreading through food pollution. Understanding the transmission of an infectious disease is crucial to predict and prevent major outbreaks of an epidemic [17]. Thus, one of the fundamental questions for the dynamical models of infectious diseases is to find the conditions that identify whether an infectious disease will exhibit a periodic outbreak or not and determine the risk factors of pollution exposure for such an outbreak in the population. In this paper, we considered a stochastic SI epidemic model in a polluted environment and incorporated the effect of environmental fluctuations as well as pollution. First, we discussed the existence of the global positive solution to the model. The main result of this paper was to obtain the sufficient conditions of the nontrivial stochastic periodic solution (see Figure 2) and the boundary periodic solution of disease extinction (see Figure 4).Compared to existing research, the main breakthrough of this paper is that we incorporated both environmental white noise as well as pollution into an SI epidemic model, which described two kinds of common phenomena in the transmission process of infectious diseases and explored the effects of environmental fluctuations (noise and pollution) on the dynamical behaviors of an epidemic. The numerical simulations based on the cases of measles showed the following:(i)
The environmental noiseσ on β may play an important role in determining the epidemic pattern: (1) it follows from Theorems 8 and 12 that the large intensity of noise σ may adversely affect the existence of the stochastic periodic of model (7) and accelerate the extinction of infectious disease. Thus, the large intensity noise may help the control of periodic outbreak of infectious disease; (2) according to Figure 2, we can see that environmental fluctuations may be responsible for the variations in the seasonal outbreak pattern of a disease in a polluted environment.(ii)
The pollution level (ri0C0∗) plays an important role in susceptible populations, in that it may reduce the number of susceptible population due to the effect of pollution; that is, the pollution causes serious harm to the susceptible population. Therefore, the pollution level may have adverse health effects on the susceptible exposed population, which is also be supported by the measles data (see Figure 5; the decreasing susceptible S(t) may be responsible for the decreasing tendency of I(t).) Moreover, by comparing the peak level of infective population with different average pollution levels ri0C0∗ (see Figure 2), we can see that the peak level of infective population increases with the levels of pollution.(iii)
From an epidemiological viewpoint, our results may provide some theoretical evidence for controlling the infectious disease. For example, in the cases of measles, the strengthening coverage of the measles vaccination and environmental quality improvement are still effective prevention measures in a polluted environment. That is to say, the lesser the population falling within the scope of susceptible subpopulation, the less infected the patients. Therefore, the demographic characteristics of susceptible population may affect the periodic outbreaks of infectious disease. In addition, the pollution control is beneficial to population health, which is consistent with the environmental research results of hemorrhagic fever [37] and influenza [38].However, this study also has several limitations:(i)
A key assumption of our model is that the pollution affects population dynamics with a linear dose-response function. Since the complicated mechanism of interaction between the pollution and population is still unclear, the dose-response parameter estimation is difficult, since the effect of pollution in vivo is not measurable for human patients. Thus, our model can not accurately describe this interaction, and the numerical simulations do not yet use the polluted data (such as PM2.5 or PM10) to check the effect of the pollution. (ii)
The times series of pollution concentration presents significant variability [39]; however, our model has not included the variability of pollution. (iii)
Some other issues also need to be considered in future, for example, the age-structured modeling [40], the impulse pollution input [41, 42], or the population with partial pollution tolerance [43].



---

*Source: 7360685-2018-05-02.xml*