Threshold Dynamics of the Stochastic SIRC Epidemic Model

In this paper, we discuss the dynamics of a stochastic SIRC epidemic model with infection rate affected by white noise. We prove that this stochastic model has a unique nonnegative solution globally. A threshold is identified. When the noise is small, the solution of the stochastic model converges to the disease-free equilibrium point of the deterministic model if 0 1 s R < , which means the basic reproductive number of the stochastic model. And if 0 1 s R > , the solution of the stochastic model fluctuates around the epidemic equilibrium of the deterministic model. When the noise is large, the disease tends to extinction. The results are illustrated by computer simulations.


Introduction
Mathematical models have become important instruments in the analysis and control of infectious diseases.Several models corresponding to epidemic with different characteristics have been proposed and discussed in literature [1] [2] [3] [4] [5].The fundamental assumption in epidemic models is that the population can be divided into distinct groups; the most common are: the susceptible (S) that are the subjects that may catch the disease; the infected (I) that are the subjects that are already infected and can spread the disease to susceptible individuals; the removed (R) that are the subjects that are immune for life.Therefore, these models are referred to as SIR models.SI, SIS, SIR models are the basic epidemic models.The differences of these models are that the population is divided into different distinct groups and the G. Y. Zhang relationships between groups in these models are different.The SI model is suitable for the disease without recurrent seizures.The SIS model can describe the patient getting sick repeatedly.The SIR model is used when the patient are immunized after infection.Other descriptions may include the presence of subjects in the quarantine state (Q) and in this situation SIRQ models are analyzed [6].The models are more and more complex with the development of epidemiology.
Recently the class of cross-immune individuals (C) in the population has been introduced [7]: it is an intermediate state between the fully susceptible state (S) and the fully protected one (R).The obtained SIRC model takes temporary partial immunity into account.The author gave a brief description of the deterministic SIRC model and calculated the basic regeneration number [8].In 2013, Iacoviello and Stasio studied the SIRC model more closely [9], and discussed the optimal control method in the case of the outbreak of the SIRC model based on the mathematical results.
In the study of stochastic epidemiological model, the authors give the proof of asymptotic stability of stochastic epidemic models and the theoretical results are verified by simulation experiments.The global uniqueness and asymptotic stability of the solution of the stochastic SIR model [10] were deduced in 2011.
Zhao studied the stochastic SIS model in 2013 [11], and the global uniqueness and asymptotic stability of the solution of the stochastic SIS model were deduced based on the basic regeneration number.Liu published a study of the stochastic SIRS model in 2017 [12]; the global uniqueness and asymptotic stability of the solution of a stochastic SIRS model are deduced under the condition that the basic regeneration number is greater than 1 or less than 1.
In the existing research, there is no literature to give the basic regeneration number of stochastic SIRC model and other related characteristics.Therefore, this paper mainly studies the stochastic SIRC model, and gives the basic regeneration number of the model.It is proved that when the basic regeneration number is less than 1, the disease-free equilibrium point is randomly asymptotically stable, and when the basic regeneration number is more than 1, the solution of the stochastic model revolves around the endemic equilibrium point of the corresponding deterministic model, and computer simulation is used to verify the correctness of the conclusions.
This paper is organized as follows.We give the model description in Section 2. In Section 3, we show that there is a unique positive solution of system (2.2).In Section 4, we investigate system (2.2) is exponential stability when the noise is large.In this case, the infective decays exponentially to zero.When the noise is small, we deduce the condition 0 1 s R < which will enable the disease to die out exponentially in Section 4 and the condition 0 1 s R > for the disease being per- sistent is given in Section 5.In Section 6, we show that 0 1 s R ≥ , the result of sto- chastic model is asymptotically stable and is oscillating around the endemic equilibrium of the deterministic model when t → ∞ .Simulations in Section 7 are reported to support the analytical results.

The SIRC Model
The SIRC model and its correlations between the variables are shown in Figure 1.The population can be divided into four parts.
The natural death rates of different compartments can be different for a long period.But in a short time, we can consider that the rates are the same.In this paper, we assumed that the natural death rate of S, I, R, C compartments are the same.
Here, ( ) A: a consistent input of new members into the population per unit time; μ: natural death rate of S, I, R, C compartments; β: transmission coefficient between compartments S and I; α: transmission coefficient between compartments I and R; ε: transmission coefficient between compartments R and C; γ: transmission coefficient between compartments C and S.
All parameter values are assumed to be nonnegative and 0 1 d < < .
Consider the following SIRC model: The threshold of system (2.1) is 1) always has the disease-free equilibrium: 0 ,0,0,0 If 0 1 R > , then 0 P is unstable and there is an epidemic equilibrium.If  , , , : 0, 0, 0, 0, In fact, epidemic models are inevitably affected by environmental white noise which is an important component in realism, because it can provide an additional degree of realism in comparison to their deterministic counterparts.Many stochastic models for epidemic population have been developed in [13]- [21].For example, stochastic models are able to take care of randomness of infectious contacts occurring in the latent and infectious periods [22].Dalal, Greenhalgh and Mao [20] have previously used the technique of parameter perturbation to examine the effect of environmental stochasticity in a model of AIDS and condom use.They found that the introduction of stochastic noise changes the basic reproduction number of the disease and can stabilize an otherwise unstable system.Nasell [23] formulated stochastic models to show that some stochastic models are a better approach to describe epidemics for a large range of realistic parameter values in comparison with their deterministic counterparts.
In this paper, we mainly discuss while β is affected by noise, ( ) where ( ) B t is standard Brownian motions with ( ) 0 0 B = , and with the in- tensity of white noise 2 0 σ > .The stochastic version corresponding to the de- terministic model (2.1) takes the following form:

Existence and Uniqueness of Positive Solution
In the biological dynamical system, the researcher is concerned about whether the positive solution of the system exists globally.Therefore, in this section, we show that the solution of system (2.2) is global and positive.Theorem 3.1.There is a unique solution 0 , 0 , 0 , 0 Proof.Since the coefficients of the equation are locally Lipschitz continuous for all given initial value 0 , 0 , 0 , 0 , where e τ is the explosion time (see [24]).To show that this solution is global, we need to show that e τ = ∞ , a.s.Let 0 0 k ≥ be sufficiently large so that where throughout this paper, we set inf ∅ = ∞ (as usual ∅ donates the empty set).According to the definition, k τ is increasing as k → ∞ .Set , , , S t I t R t C t + ∈ R , a.s.for all 0 t ≥ .In other words, to complete the proof all we need to show is that e τ = ∞ , a.s.If statement is false, then there exist a pair of constants 0 T > and ( ) Hence there is an integer , , , : 0, 0, 0, 0, , . .
is a positively invariant set of system (2.2) on * Γ , which is similar to Γ of sys- tem (2.1).From now on, we always assume that Define a C 2 -function W: The nonnegativity of this function can be seen form T > be arbitrary.Applying the Itô formula, we obtain ( ) Therefore, we get Formula (3.2) It then follows from (3.1) and ( The former formula is contradictory and we must therefore have τ ∞ = ∞ , a.s.
For convenience we introduce the notation:

Extinction
In this section, we investigate the conditions for the extinction of disease.Proof.An integration of system (2.2) yields, the Formula (4.1) is as follows Applying the Itô formula, we obtain Then Formula (4.2) is as follows If the condition (a) is satisfied, the following statement is tenable.
( ) If the conditions (b) are satisfied, Formula (4.5) can be rewritten in another way as follows

I t I t S dC S dC S r dC r B r
From Formula (3.1), we have So Formula (5.1) can be rewritten as Taking inferior limit of Formula (5.4), we have Formula (

Asymptotic Property around the Epidemic Equilibrium
Consider the basic reproductive number of stochastic model 0 s R , if 0 1 s R < , the disease becomes extinct and the solution of stochastic model converges to the disease-free equilibrium point of the deterministic model when t → ∞ .There- fore, the following part is to introduce the characteristic of stochastic model when the disease survives.Theorem 6.1.When 0 1 s R ≥ , the result of stochastic model is asymptotically stable and is oscillating around the endemic equilibrium of the deterministic model when t → ∞ .That is , Formula (6.2) can be rewritten like this

Conclusions
The t denotes the number of members who are susceptible to an infec- tion at time t. ( ) I t denotes the number of members who are infective to an in- fection at time t. ( ) R t denotes the number of members who are immune to an infection at time t. ( ) C t denotes the number of members who are cross-immune to an infection at time t.The parameters in the model are summarized in the following list:

and the solu- tion will remain in 4 +R
with probability 1, namely the solution of system (2.2) has the following property: And the disease-free equilibrium ( ) can be rewritten in the same way

Figure 4 .
Figure 4.The sample paths of model (2.2) with 0 1 R > and

1 sR < . And if 0 1 sR
solution of our stochastic SIRC model is affected by basic reproduction number and noise strength.When the noise is small, the solution of the stochastic model converges to the disease-free equilibrium point of the deterministic model if 0 > , the solution of the stochastic model fluctuates around the epidemic equilibrium of the deterministic model.When the noise is large, the disease tends to extinction.In Section 6, we get that the solution of the stochastic model fluctuates around the epidemic equilibrium of the deterministic model.Our next work is to prove this conclusion is valid by probability; that is ) (