Stochastic SIRC epidemic model with time-delay for COVID-19

Environmental factors, such as humidity, precipitation, and temperature, have significant impacts on the spread of the new strain coronavirus COVID-19 to humans. In this paper, we use a stochastic epidemic SIRC model, with cross-immune class and time-delay in transmission terms, for the spread of COVID-19. We analyze the model and prove the existence and uniqueness of positive global solution. We deduce the basic reproduction number R0s\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}${\mathcal{R}}_{0}^{s}$\end{document} for the stochastic model which is smaller than R0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}${\mathcal{R}}_{0}$\end{document} of the corresponding deterministic model. Sufficient conditions that guarantee the existence of a unique ergodic stationary distribution, using the stochastic Lyapunov function, and conditions for the extinction of the disease are obtained. Our findings show that white noise plays an important part in controlling the spread of the disease; When the white noise is relatively large, the infectious diseases will become extinct; Re-infection and periodic outbreaks can occur due to the existence of feedback time-delay (or memory) in the transmission terms.


Introduction
The ongoing pandemic coronavirus disease  has become a worldwide emergency. This infectious disease is spreading fast, endangering a large number of people's health, and thus immediate actions and intensive studies are needed to control the disease in communities [1]. COVID-19 is the seventh member of the coronavirus (CoV) family, such as MERS-CoV and SARS-CoV [2]. Although SARS-CoV was more deadly, it was much less infectious than COVID- 19. There have been no outbreaks of SARS anywhere in the world since 2003. The symptoms of the COVID-19 infection include cough, fever, tiredness, diarrhea, and shortness of breath. Mostly in severe cases, COVID-19 causes pneumonia and death [3]. The primary studies show that the incubation period of COVID-19 is between 3-14 days or longer [4]. Additionally, the average of basic reproduction number R 0 for COVID-19 is about 2-2.8. The disease may still be infectious in the latent infection period. Studies to date suggest that the virus is very serious and spreads fast from person to person through close contact and respiratory droplets rather than through the air [4]. Table 1 shows the incubation period of several common infectious diseases. 1-10 days [31] Mathematical modeling of the infectious diseases has an important role in the epidemiological aspect of disease control [5][6][7][8]. Several epidemic models, with various characteristics, have been described and investigated in the literature. Most of these models are based on susceptible-infected-removed (SIR) model. Casagrandi et al. [9] introduced a SIRC model to describe the dynamical behaviors of influenza A by inserting a new compartment, namely cross-immunity (C) component a of people who have been recovered after being infected by different strains of the same viral subtype in previous years. The component C describes an intermediate state between the susceptible S and the recovered R. Rihan et al. [10] investigated the qualitative behaviors of a fractional-order SIRC model for salmonella bacterial infection. Recently in [11], the authors provided a deterministic SEIR epidemic model of fractional order to describe the dynamics of COVID- 19. In other descriptions, quarantine state (Q) may be included in the presence of subjects, such as SIRQ models [12].
In fact, stochastic perturbation factors, such as precipitation, absolute humidity, and temperature, have a significant impact on the infection force of all types of virus diseases to humans. Taking this into consideration enables us to present randomness into deterministic biological models to expose the environmental variability effect, whether it is environmental fluctuations in parameters or random noise in the differential systems [13][14][15][16][17]. Moreover, stochastic models give an extra degree of freedom and realism in comparison with their corresponding deterministic models. Stochastic population dynamics perturbed by white noise (or Brownian motion) has been studied extensively by many authors [18][19][20]. It has been investigated in [21] that an environmental Brownian noise can suppress explosions in population dynamics. Yuan et al. [22] discussed the results of stochastic viral infection, immune response dynamics and analyzed the human immunodeficiency virus infection. In [23], the author investigated the existence results of ergodic distribution for stochastic hepatitis B virus model based on Lyapunov function. In [24], the authors explored the dynamics of SIR epidemic model with environmental fluctuations. Additionally, they calculated a threshold parameter to demonstrate the persistence and extinction of the disease. Recently, Lakshmi et al. [25] identified some environmental factors such as geographic location of the countries, the upcoming climate, atmospheric temperature, humidity, sociobiological factors, etc., that influence the global spread of COVID-19.
Up-to-date studies have reported that there are many COVID-19 carriers who are not suffering from the disease. This may be due to cross-immunity of other virus survivors, people who have been recovered from the virus, such as other stains of coronavirus, H1N1, or influenza A. It has been reported in [2] that "T-cells that target SARS-CoV2, the virus that causes COVID-19, in the blood of people who had recovered from a coronavirus infection." Accordingly, in the present paper, we investigate an SIRC epidemic model of crossimmune class for the dynamics of COVID-19 infection among groups. We include timedelay in the transmission terms to represent the incubation period of the virus (the time between infection and symptom onset). We also incorporate white noise type of perturbations to reveal the effect of environmental fluctuations and variability in parameters. Based on the existing literature, this is the first work dealing with the persistence and extinction of a stochastic epidemic model for the COVID-19 infection. We investigate the impact of small and large values of white noise in the persistence and extinction of the disease. We also derive the existence results of stationary distribution and extinction of the disease using a novel combination of stochastic Lyapunov functional. This paper is presented as follows: We provide a stochastic SIRC model with time-delay in Sect. 2. In Sect. 3, we study the existence and uniqueness of a global positive solution for the stochastic delayed SIRC model. In Sect. 4, a stationary distribution and extinction analysis of the underlying model are investigated. Some virtual numerical examples are presented in Sect. 5. Finally, concluding remarks are given in Sect. 6.

Stochastic SIRC epidemic model
For the spread of the COVID-19 disease in humans, we classify the population into four categories: S(t), I(t), R(t), and C(t) are the proportion of susceptible, infected, recovered, and cross-immune ones at time t, respectively. Let N(t) = S(t)+I(t)+R(t)+C(t) be the total population. At this stage, we believe that a SIRC model efficiently describes the mechanism for the spreading of the COVID-19 virus. The classical SIRC model [9,32] takes the formṠ We incorporate a discrete time-delay τ into the SIRC model to represent the incubation period which is about 3-14 days [4]. All the parameters appearing in the model are nonnegative, see Table 2. In the absence of cross-immunity i.e. (1σ = 0), the SIRC model curtails to the SIRS model since the two individuals S and C become immunologically indistinguishable. Figure 1 shows the scheme of SIRC model. Mortality rate in every compartment assumed to be equal to the rate of newborn in the population [9] β Rate at which the cross-immune population becomes susceptible again  Time-delay τ > 0 is incorporated in the transmission terms to represent the incubation period of the viral infection, the time between infection and symptom onset. The current studies show that the average/median of incubation period of early confirmed cases of COVID-19 is about 5.5 days, which is similar to SARS-CoV. Presence of time-delay in the model may cause periodic solutions many times for different time-delay values τ [33].
Model (1) has a disease-free equilibrium E 0 = [1, 0, 0, 0] and an endemic equilibrium , and I * is a root of quadratic equation pI 2 + qI + r = 0, where Here R 0 = ξ η+α is known as the basic reproduction number of the deterministic model. In fact, there is an increasing indication that superior consistency with some phenomena can be contributed if the effects of environmental noises in the system are taken into account [34]. Epidemic model (1) assumes that the observed dynamics are driven exclusively by internal deterministic cases. Ignoring environmental variability in the modeling may affect the dynamics of the model and transmission of the disease. Accordingly, there is a need to extend the deterministic systems described by differential equations into stochastic differential equations (SDEs), where related parameters are modeled as suitable stochastic processes, added to the driving system equations.
From the mathematical and biological point of view, there are some assumptions to incorporate stochastic perturbations into the epidemiological model, such as Markov chain process, parameter perturbations, white noise type, etc. Here, we incorporate white noise type perturbation into model (1), which is proportional to the S, I, R, C classes, so that where W 1 (t), W 2 (t), W 3 (t), and W 4 (t) stand for the independent Brownian motions. ν 2 1 , ν 2 2 , ν 2 3 , and ν 2 4 represent the intensity of the environmental white noises, ν i > 0 (i = 1, 2, 3, 4) subject to the following initial conditions: such that C is a family of Lebesgue integrable functions from [-τ , 0] into R 4 + .

Existence and uniqueness of the positive solution
To investigate the dynamical characteristics of SDDEs (2), the first consideration is to verify if system (2) has a unique global positive solution. As the coefficients of system (2) satisfy the local Lipschitz condition together with the linear growth condition [35,36], there exists a unique local solution. Now, we need to prove that the solution is positive and global using the Lyapunov analysis method [36].

Theorem 1 System
(2) has a unique positive solution (S(t), I(t), R(t), C(t)) on t ≥ -τ , and the solution will remain in R 4 + for the given initial condition (3) with probability one. (3), as the coefficients of system (2) satisfy the local Lipschitz condition, so system (2) has a unique local solution (S(t),

Proof 1 For any initial value
where τ e represents the explosion time [36]. Our aim is to show that this solution is global i.e. τ e = ∞ a.s. Assume n 0 ≥ 1 to be sufficiently large such that S(θ ), I(θ ), R(θ ), and C(θ ) (θ ∈ [-τ , 0]) are lying in the interval [ 1 n 0 , n 0 ]. For each n ≥ n 0 , n ∈ N, define the stopping time we fix inf φ = ∞ (φ is the empty set). Apparently, τ n is increasing as n → ∞. Assume τ ∞ = lim n→∞ τ n , then τ ∞ ≤ τ e a.s. Therefore, we need to show that τ ∞ = ∞ a.s., then τ e = ∞ a.s. and (S(t), I(t), R(t), C(t)) ∈ R 4 + a.s. for all t ≥ -τ . If it is erroneous, there is a pair ∈ (0, 1) and T > 0 such that P{τ ∞ ≤ T} > . Then there is an integer n 1 ≥ n 0 such that We define a C 2 -function V : R 4 + → R + as follows: where κ > 0 is a constant to be determined. By Ito's formula, we can obtain Let κ = α-ξ ξ , then we have where M > 0 is a constant which is independent of S(t), I(t), R(t), and C(t). Therefore, Integrating (6) from 0 to τ n ∧ T = min{τ n , T} and then taking the expectation E on both sides, we have Let n = {τ n ≤ T}, for n ≥ n 1 and in view of (4), we obtain P( n ) ≥ such that, for every ω ∈ n , there is at least one of S(τ n , ω), I(τ n , ω), R(τ n , ω), or C(τ n , ω) equaling either n or 1 n , and then we obtain According to (7), we get where 1 n represents the indicator function of n . Letting n → ∞ yields which leads to a contradiction. It can be concluded that τ ∞ = ∞ a.s., which proves the theorem.

Existence of ergodic stationary distribution
Herein, we construct a suitable stochastic Lyapunov function to study the existence of a unique ergodic stationary distribution of the positive solutions to system (2). First, let us assume that X(t) is a regular time-homogenous Markov process in R d illustrated by the SDDE The diffusion matrix of the process X(t) is

Lemma 1 ([37]) The Markov process X(t) has a unique ergodic stationary distribution π(·)
if there exists a bounded domain U ⊂ R d with regular boundary , and Define the reproduction number of the stochastic model as follows: Theorem 2 Assume that R s 0 > 1 and η - Proof 2 First, we need to validate conditions (i) and (ii) of Lemma 1. To prove condition (i), the diffusion matrix of model (2) is described as follows: Then the matrix is positive definite for any compact subset of R 4 + , then condition (i) of Lemma 1 is satisfied.

Extinction
In order to show the extinction of the disease, we go through the following lemmas.  , then the solution of (2) satisfies the following: 2 )} < 0 and lim t→∞ S = 1 a.s.
Since R s 0 < 1, therefore, by taking the superior limit of both sides of (34), we have which implies that lim t→∞ I(t) = 0, lim t→∞ R(t) = 0, lim t→∞ C(t) = 0 a.s., which confirms that the disease I can die out with probability one. It is easy, by using (33) and (35), to show that lim t→∞ S = 1 a.s.

Numerical simulations and discussions
Numerical simulations are given to validate our theoretical results through Euler-Maruyama method for SDDEs reported in [40,41] to numerically solve SDDEs (2). The discretization transformation takes the form The independent Gaussian random variables denoted as ζ i,j (i = 1, 2, 3, 4), which follow the distribution N(0, 1), the time-delay is defined as τ = m t, m is an integer, and the step size t. Let ν i > 0, (i = 1, 2, 3, 4) be the white noise values.  = -0.0195 < 0. In this case, the conditions of Theorem 2 are not satisfied. From Fig. 3, we can clearly find that the disease goes to extinction. In Fig. 4 time-delay is increased to τ = 2.5, with white noises ν 1 = 0.01, ν 2 = 0.2, ν 3 = 0.02, ν 4 = 0.03, other parameter values are the same as in Fig. 3. Therefore R s 0 < 1 and η -      Thus, according to Theorem 3, the disease dies out exponentially with probability one; see Fig. 7. Therefore, the smaller values of white noise ensure the existence of unique stationary distribution, which gives the persistence of the disease; while larger values of white noise can lead to disease extinction.

Conclusion
In this work, we provided a stochastic SIRC epidemic model with time-delay for the new strain coronavirus COVID-19. The stochastic components, due to environmental variability, are incorporated in the model as Gaussian white noise. We established some sufficient conditions for persistence and extinction in the mean of the disease. The model has a unique stationary distribution which is ergodic if the intensity of white noise is small. Introduction of noise in the deterministic SIRC model modifies the basic reproductive number R 0 giving rise to a new threshold quantity R s 0 . It has been proved that the disease dies out if R s 0 < 1 < R 0 . On the other hand, if R s 0 > 1 and R 0 > 1, the disease persists with both models, but with different behaviors. In other words, extinction of the infection possibly occurs when R s 0 < 1 < R 0 and the intensity of white noise is large. This would not happen in the deterministic models. The potential of using stochastic SIRC model for COVID-19 is to consider the environmental fluctuation that all affects the spread of the virus. The periodicity of outbreaks is possible due to the presence of time-delay (memory) in the transmission terms.
The authors believe that the stochastic SIRC model is an attempt to understand epidemiological characteristics of COVID-19. The model provides new insights into epidemiological situations when the environmental noise (perturbations) and cross-immunity are considered in the COVID-19 epidemic models. The combination of white noise and timedelay, in the epidemic model, has a considerable impact on the persistence and extinction of the infection and enriches the dynamics of the model. This work can be extended to include control variables for vaccination, treatment, and/or quarantine actions. A more sophisticated model is also required to investigate the dynamics of COVID-19 with immune system in cells level [42]. Fractional derivatives can also be included in the model to consider long-run memory [43,44].