Issue
Math. Model. Nat. Phenom.
Volume 15, 2020
Coronavirus: Scientific insights and societal aspects
Article Number 50
Number of page(s) 14
DOI https://doi.org/10.1051/mmnp/2020040
Published online 11 November 2020

© The authors. Published by EDP Sciences, 2020

Licence Creative CommonsThis is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1 Introduction

Coronavirus disease 2019 (COVID-19), reclassified as a pandemic by the World Health Organization (WHO) on March 11, 2020 [14], is an infectious disease caused by a new type of virus belonging to the coronaviruses family and recently named severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) [4]. All the countries affected by this disease have taken many preventive measures, including containment. The containment established by the Moroccan government and the public authorities at the right time made it possible to avoid the worst: according to the minister of Health, at least 6000 lives were saved thanks to the measures adopted to face the spread of this pandemic [13]. The resistance measures, regarded as necessary and urgent, cannot be sustainable.

Actually, the deconfinement is a new stage entered by the COVID-19 pandemic. Therefore, several countries strategically planned their deconfinement strategies. The extension of the state of emergency in Morocco until May 20, 2020 will no doubt have economic repercussions. If Morocco won the first round, or at least limited the consequences, especially in terms of limiting the pandemic and health management of the situation, the second seems difficult and complex. Indeed, it must not only be well thought out but also its axes of resistance have to be well-identified. In this context, all efforts should be focused on stabilizing the economy by intelligently relying on resources. Economic deconfinement is part of the solution and should be gradual and concerted. Indeed, it is absurd to think that the return to the normality is in the near months, because the unavailability of an effective vaccine implies that the virus will always be with us in the near future, which poses a risk for the population. This economic deconfinement should be prepared and accompanied by other related measures, in particular under health, security, education and social assistance. In this period of general crisis, the response must try to mitigate the impacts on priority sectors, such as agriculture, agrifood, transport and foreign trade, in relation to imports that are vital to the Moroccan economy. The challenge is to ensure resistance and a continuity of value creation while preventing a sector from being detached from the economic body. So to speak, priority must be given to vital sectors whose health directly affects all Moroccan activity, while protecting those bordering on chaos. According to the deconfinement strategy, which is applied by the Moroccan authorities, it is mandatory to study the occurrence of an eventual second wave and it’s magnitude.

Mathematical modeling through dynamical systems plays an important role to predict the evolution of COVID-19 transmission [9, 16]. However, while taking into account the deconfinement policies, the environmental effects and the social fluctuations should not be neglected in such a mathematical study in order to describe well the dynamics and consider an additional degree of realism [2, 3, 10, 18, 21]. For these reasons, we describe here the dynamics of the deconfinement strategy by a new D-COVID-19 model, governed by delayed stochastic differential equations (DSDE), as follows: { dS(t)=(ρC(t)δS(t)β(1u)S(t)Is(t)N)dtσ1(1u)S(t)Is(t)NdB1(t)+σ2(C(t)S(t))dB2(t),dC(t)=(δS(t)ρC(t))dt+σ2(S(t)C(t))dB2(t),dIs(t)=(βɛ(1u)S(tτ1)Is(tτ1)NαIs(t)(1α)(μs+ηs)Is(t))dt+σ1(ɛ(1u)S(tτ1)Is(tτ1)N)dB1(t)+σ3(μs+ηs1)Is(t)dB3(t), dIa(t)=(β(1ɛ)(1u)S(tτ1)Is(tτ1)NηaIa(t))dt+σ1(1ɛ)(1u)S(tτ1)Is(tτ1)NdB1(t),dFb(t)=(αγbIs(tτ2)(μb+rb)Fb(t))dt+σ3γbIs(tτ2)dB3(t),dFg(t)=(αγgIs(tτ2)(μg+rg)Fg(t))dt+σ3γgIs(tτ2)dB3(t),dFc(t)=(αγcIs(tτ2)(μc+rc)Fc(t))dt+σ3γcIs(tτ2)dB3(t), dR(t)=(ηs(1α)Is(tτ3)+ηaIa(tτ3)+rbFb(tτ4)+rgFg(tτ4)+rcFc(tτ4))dtσ3ηsIs(tτ3)dB3(t),dM(t)=(μs(1α)Is(tτ3)+μbFb(tτ4)+μgFg(tτ4)+μcFc(tτ4))dtσ3μsIs(tτ3)dB3(t),\begin{align*}\textrm{d}S(t)&=\left(\rho C(t)-\delta S(t)-\beta(1-u)\dfrac{S(t)I_{s}(t)}{N}\right)\textrm{d}t -\sigma_{1}(1-u)\dfrac{S(t)I_{s}(t)}{N}\textrm{d}B_{1}(t)+ \sigma_{2}(C(t)-S(t))\textrm{d}B_{2}(t),\nonumber\\ \textrm{d}C(t)&=\left(\delta S(t)- \rho C(t)\right)\textrm{d}t + \sigma_{2}(S(t)-C(t))\textrm{d}B_{2}(t),\nonumber \\ \textrm{d}I_{s}(t)&=\left(\beta\epsilon (1-u)\dfrac{S(t-\tau_1)I_{s}(t-\tau_1)}{N}-\alpha I_{s}(t)-(1-\alpha)(\mu_s+\eta_{s})I_s(t)\right)\textrm{d}t\nonumber \\ &\quad+\sigma_{1}\left(\epsilon (1-u)\dfrac{S(t-\tau_1)I_{s}(t-\tau_1)}{N}\right)\textrm{d}B_{1}(t) +\sigma_{3}(\mu_s+\eta_{s}-1)I_s(t)\textrm{d}B_{3}(t),\nonumber\\ \textrm{d}I_{a}(t)&=\left(\beta(1-\epsilon)(1-u)\dfrac{S(t-\tau_1)I_{s}(t-\tau_1)}{N} -\eta_{a}I_{a}(t)\right)\textrm{d}t +\sigma_{1}(1-\epsilon)(1-u)\dfrac{S(t-\tau_1)I_{s}(t-\tau_1)}{N}\textrm{d}B_{1}(t),\nonumber\\ \textrm{d}F_{b}(t)&=\bigg(\alpha\gamma_{b}I_{s}(t-\tau_2)-\big(\mu_b+r_b\big)F_{b}(t)\bigg)\textrm{d}t +\sigma_{3}\gamma_{b}I_{s}(t-\tau_2)\textrm{d}B_{3}(t),\nonumber\\ \textrm{d}F_{g}(t)&=\bigg(\alpha\gamma_{g}I_{s}(t-\tau_2)-\big(\mu_g+r_g\big)F_{g}(t)\bigg)\textrm{d}t +\sigma_{3}\gamma_{g}I_{s}(t-\tau_2)\textrm{d}B_{3}(t),\nonumber\\ \textrm{d}F_{c}(t)&=\bigg(\alpha\gamma_{c}I_{s}(t-\tau_2)-\big(\mu_c+r_c\big)F_{c}(t)\bigg)\textrm{d}t +\sigma_{3}\gamma_{c}I_{s}(t-\tau_2)\textrm{d}B_{3}(t),\nonumber\\ \textrm{d}R(t)&=\bigg(\eta_{s}(1-\alpha)I_s(t-\tau_3)+\eta_{a}I_{a}(t-\tau_3) +r_{b}F_{b}(t-\tau_4)+r_g F_{g}(t-\tau_4) +r_c F_{c}(t-\tau_4)\bigg)\textrm{d}t\nonumber \\ &\quad -\sigma_3\eta_{s}I_{s}(t-\tau_3)\textrm{d}B_3(t),\nonumber\\ \textrm{d}M(t)&=\bigg(\mu_s(1-\alpha)I_s(t-\tau_3)+\mu_bF_{b}(t-\tau_4) +\mu_gF_{g}(t-\tau_4)+\mu_cF_{c}(t-\tau_4)\bigg)\textrm{d}t -\sigma_{3}\mu_{s}I_{s}(t-\tau_3)\textrm{d}B_{3}(t),\nonumber \\ \end{align*}(1.1)

where S represents the susceptible sub-population, which is not infected and has not been infected before but is susceptible to develop the disease if exposed to the virus; C is the confined sub-population; Is is the symptomatic infected sub-population, which has not yet been treated, it transmits the disease, and outside of proper support it can progress to spontaneous recovery or death; Ia is the asymptomatic infected sub-population who is infected but does not transmit the disease, is not known by the health system and progresses spontaneously to recovery; Fb, Fg and Fc are the patients diagnosed, supported by the Moroccan health system and under quarantine, and subdivided into three categories: benign, severe, critical forms, respectively. Finally, R and M are the recovered and died classes, respectively. At each instant of time, the equation D(t):= μ s (1α) I s (t τ 3 )+ μ b F b (t τ 4 )+ μ g F g (t τ 4 )+ μ c F c (t τ 4 ) σ 3 μ s I s (t τ 3 ) Δ B 3 (t) Δt = ΔM(t) Δt \[ D(t):=\mu_s(1-\alpha)I_s(t-\tau_3)+\mu_bF_{b}(t-\tau_4) +\mu_gF_{g}(t-\tau_4)+\mu_cF_{c}(t-\tau_4) -\sigma_3\mu_{s}I_{s}(t-\tau_3)\dfrac{\UpDelta B_3(t)}{\UpDelta t} =\dfrac{\UpDelta M(t)}{\UpDelta t} \]

gives the number of the new dead due to disease. The parameter 1 − u represents the level of measures undertaken on the susceptible population while δ is the confinement rate and ρ represents the deconfinement rate. We adopt the bilinear incidence rate to describe the infection of the disease and use the parameter β to denote the transmission rate. It is reasonable to assume that the infected individuals are subdivided into individuals with symptoms and others without symptoms, for which we employ the parameter ɛ to denote the proportion for the symptomatic individuals and 1 − ɛ for the asymptomatic ones. The parameter α measures the efficiency of public health administration for hospitalization. Diagnosed symptomatic infected population is completely distributed into one of the three forms Fb, Fg and Fc, by the rates γb, γg and γc, respectively. Then, γb + γg + γc = 1. The mean recovery period of these forms are denoted by 1∕rb, 1∕rg and 1∕rc, respectively. The latter forms die also with the rates μb, μg and μc, respectively.Symptomatic infected population, which is not diagnosed, moves to the recovery compartment with a rate ηs or dies with a rate μs. On the other hand, asymptomatic infected population moves to the recovery compartment with a rate ηa. The time delays τ1 and τ2 denote the incubation period and the period of time needed before the charge by the health system, respectively. The time delays τ3 and τ4 denote the time required before the death of individuals coming from the compartments Is and the three forms Fb, Fg and Fc, respectively. Here, B1(t), B2 (t) and B3 (t) are independent standard Brownian motions defined on a complete probability space (Ω,F,) $(\UpOmega,\mathcal{F},\mathbb{P})$ with a filtration { F t } t0 $\{\mathcal{F}_{t}\}_{t\geq0}$ and satisfying the usual conditions, that is, they are increasing and right continuous while F 0 $\mathcal{F}_{0}$ contains all P-null sets and σi represents the intensity of Bi, i = 1, 2, 3. The schematic diagram of our extended model is illustrated in Figure 1.

Remark 1.1. The stochasticity is introduced in model (1.1) by perturbing the most sensitive parameters: β, α, δ, and ρ.

Remark 1.2. For the sake of simplicity, we have assumed that the parameters δ and ρ are perturbed with the same intensities, that is, we assume that Moroccan individuals possess the same behaviors and reactions towards the authorities instructions.

Remark 1.3. Note that the multipliers of Is(tτ2) terms are the same as σ3dB3(t) although they are premultiplied by different constants. Indeed, the portion of diagnosed symptomatic infected population is completely distributed into the three forms Fb, Fg, and Fc, by the rates γb, γg and γc, respectively. Then, γb + γg + γc = 1. In addition, we assume that the parameter α, which measures the efficiency of public health administration for hospitalization, undergoes random fluctuations.

Remark 1.4. For illustration and clarification purposes, let us suppose, as an example, that the disease progression is started from 15th of March and value of τ1 is 5. Then a susceptible individual, after contact with an infected one at instant t, becomes himself infected at instant t + τ1. Suddenly, the compartment of the infected is fed at the instant t by the susceptible infected at the instant tτ1. Therefore, in the considered situation, when the infection starts at March 15, the term βɛ(1 − u)S(T)Is(T)∕N of new infected is equal to zero on 12th, 13th and 14th March, due to the absence of the infection.

Remark 1.5. Temporarily asymptomatic individuals are included in the class Is of symptomatic, while individuals in Ia, who are permanently asymptomatic, will remain asymptomatic until recovery and will not spread the virus, a fact which has been recently confirmed by the World Health Organization.

For biological reasons, we assume that the initial conditions of system (1.1) satisfy: S(θ)= ϕ 1 (θ)0,C(θ)= ϕ 2 (θ)0, I a (θ)= ϕ 3 (θ)0, I s (θ)= ϕ 4 (θ)0, F b (θ)= ϕ 5 (θ)0, F g (θ)= ϕ 6 (θ)0, F c (θ)= ϕ 7 (θ)0,R(θ)= ϕ 8 (θ)0,M(θ)= ϕ 9 (θ)0, \begin{equation*}\begin{array}{ll} S(\theta) =\phi_{1}(\theta)\geq 0, \quad C(\theta) =\phi_{2}(\theta)\geq 0, \quad I_a(\theta)=\phi_{3}(\theta)\geq 0,\\ I_s(\theta)=\phi_{4}(\theta)\geq 0, \quad F_b(\theta) =\phi_{5}(\theta)\geq 0, \quad F_g(\theta)=\phi_{6}(\theta)\geq 0,\\ F_c(\theta)=\phi_{7}(\theta) \geq 0, \quad R(\theta)=\phi_{8}(\theta)\geq 0, \quad M(\theta)=\phi_{9}(\theta)\geq 0, \end{array} \end{equation*}(1.2)

where θ ∈ [−τ, 0] and τ = max{τ1, τ2 τ3, τ4}.

The rest of the paper is organized as follows. Section 2 deals with the existence and uniqueness of a positive global solution that ensures the well-posedness of the D-COVID-19 model (1.1). A sufficient condition for the extinction isestablished in Section 3. Then, some numerical scenarios, to assess the effectiveness of the adopted deconfinement strategy, are presented in Section 4. The paper ends up with Section 5 of conclusion.

thumbnail Figure 1

Schematic diagram of model (1.1).

2 Existence and uniqueness of a positive global solution

Let us denote + 9 :={( x 1 , x 2 , x 3 , x 4 , x 5 , x 6 , x 7 , x 8 , x 9 ) x i >0,i=1,2,,9} $\mathbb{R}^9_+ :=\{(x_1,x_2,x_3,x_4,x_5,x_6,x_7,x_8,x_9)\mid x_i>0, \, i=1,2,\ldots,9\}$. We begin by proving the following result.

Theorem 2.1. For any initial value satisfying condition (1.2), there is a unique solution

x(t)=(S(t),C(t), I s (t), I a (t), F b (t), F g (t), F c (t),R(t),M(t)) \[ x(t)=(S(t),C(t),I_s(t),I_a(t),F_b(t),F_g(t),F_c(t),R(t),M(t)) \]

to the D-COVID-19 model (1.1) that remains in + 9 $\mathbb{R}^9_+$ with probability one.

Proof. Since the coefficients of the Stochastic Differential Equations with several delays (1.1) are locally Lipschitz continuous, it follows from [12] that for any square integrable initial value x(0) + 9 $x(0)\in\mathbb{R}^9_&#x002B;$, which is independent of the considered standard Brownian motion B, there exists a unique local solution x(t) on t ∈ [0, τe), where τe is the explosion time. For showing that this solution is global, knowing that the linear growth condition is not verified, we need to prove that τe = . Let k0 > 0 be sufficiently large for 1 k 0 <x(0)< k 0 $\dfrac{1}{k_0}<x(0)<k_0$. For each integer kk0, we define the stopping time τ k :=inf{ t[0, τ e )/ x i (t)( 1 k ,k )for somei=1,2,3 } $\tau_k := \inf\left\lbrace t\in[0,\tau_e)/x_i(t)\notin \left(\dfrac{1}{k},k\right)\; \text{for some}\; i=1,2,3\right\rbrace$, where inf = .

It is evident that τkτe. Let T > 0, and define the twice differentiable function V on + 3 + $\mathbb{R}^3_&#x002B;\rightarrow \mathbb{R}^&#x002B;$ as follows: V(x):= ( x 1 + x 2 + x 3 ) 2 + 1 x 1 + 1 x 2 + 1 x 3 . \[ V(x):=(x_1&#x002B;x_2&#x002B;x_3)^2&#x002B;\frac{1}{x_1}&#x002B;\frac{1}{x_2}&#x002B;\frac{1}{x_3}. \]

By Itô’s formula, for any 0 ≤ tτkT and k ≥ 1 we have dV(x(t))=LV(x(t))dt+σ(x(t))d B t , \[ dV(x(t))=LV(x(t))\textrm{d}t&#x002B;\sigma(x(t))\textrm{d}B_t, \]

where L is the differential operator of function V : L V(x(t))=( 2(S(t)+C(t)+ I s (t)) 1 S 2 (t) )( ρC(t)δS(t)β(1u) S(t) I s (t) N )     + 1 2 ( 2+ 2 S 3 (t) )( ( σ 1 (1u) S(t) I s (t) N ) 2 + ( σ 2 (C(t)S(t))) 2 )     +( 2(S(t)+C(t)+ I s (t)) 1 C 2 (t) )(δS(t)ρC(t))+ 1 2 ( 2+ 2 C 3 (t) ) ( σ 2 (C(t)S(t))) 2     + 1 2 ( 2+ 2 C 3 (t) )( σ 2 (C(t)S(t)) ) 2     +( 2(S(t)+C(t)+ I s (t)) 1 I s 2 (t) )( βɛ(1u) S(t τ 1 ) I s (t τ 1 ) N α I s (t)(1α)( μ s + η s ) I s (t) )     + 1 2 ( 2+ 2 I s 3 (t) )( ( σ 3 I s (t)(1 μ s η s ) ) 2 + ( σ 1 ɛ(1u) S(t τ 1 ) I s (t τ 1 ) N ) 2 ). \begin{align*} \begin{split} L&V(x(t)) = \left(2(S(t)&#x002B;C(t)&#x002B;I_s(t))-\dfrac{1}{S^2(t)}\right) \left(\rho C(t)-\delta S(t)-\beta (1-u)\dfrac{S(t) I_s(t)}{N}\right)\\ &\qquad\qquad&#x002B; \dfrac{1}{2}\left(2&#x002B;\frac{2}{S^3(t)}\right) \left(\left( -\sigma_1(1-u)\dfrac{S(t) I_s(t)}{N}\right)^2 &#x002B;(\sigma_2(C(t)-S(t)))^2\right) \\ &\qquad\qquad&#x002B; \left(2(S(t)&#x002B;C(t)&#x002B;I_s(t))-\dfrac{1}{C^2(t)}\right) (\delta S(t) -\rho C(t))&#x002B;\dfrac{1}{2}\left(2&#x002B;\frac{2}{C^3(t)}\right) (-\sigma_2 (C(t)-S(t)))^2\\ &\qquad\qquad&#x002B; \dfrac{1}{2}\left(2&#x002B;\frac{2}{C^3(t)})(-\sigma_2(C(t)-S(t))\right)^2\\ &\qquad\qquad&#x002B; \left(2(S(t)&#x002B;C(t)&#x002B;I_s(t))-\frac{1}{I^2_s(t)}\right) \left( \beta \epsilon (1-u)\frac{S(t-\tau_1) I_s(t-\tau_1)}{N}-\alpha I_s(t) -(1-\alpha)(\mu_s&#x002B;\eta_s)I_s(t)\right) \\ &\qquad\qquad&#x002B; \dfrac{1}{2}\left(2&#x002B;\dfrac{2}{I^3_s(t)}\right) \left(\left( -\sigma_3 I_s(t)(1-\mu_s-\eta_s)\right) ^2&#x002B;\left(\sigma_1 \epsilon (1-u)\frac{S(t-\tau_1) I_s(t-\tau_1)}{N}\right) ^2\right). \end{split} \end{align*}

Thus, LV(x(t)) 2(S(t)+C(t)+ I s (t))ρC(t)+ δ S(t) + β(1u)S(t) I s (t) N S 2 (t) +( 1+ 1 S 3 (t) )( ( σ 1 (1u) S(t) I s (t) N ) 2 + ( σ 2 (C(t)S(t))) 2 ) +2(S(t)+C(t)+ I s (t))δS(t)+ ρ S(t) +( 1+ 1 C 3 (t) ) ( σ 2 (C(t)S(t))) 2 +2(S(t)+C(t)+ I s (t))βɛ(1u) S(t τ 1 ) I s (t τ 1 ) N + α I s (t) +(1α)( μ s + η s ) 1 I s (t) +( 1+ 1 I s 3 (t) )( ( σ 3 I s (t)(1 μ s η s ) 2 + ( σ 1 ɛ(1u) S(t τ 1 ) I s (t τ 1 ) N ) 2 ). \begin{align*} LV(x(t)) \leq & 2(S(t)&#x002B;C(t)&#x002B;I_s(t))\rho C(t)&#x002B;\frac{\delta}{S(t)} &#x002B;\frac{\beta (1-u)S(t) I_s(t)}{NS^{2}(t)}\\ &&#x002B; \left(1&#x002B;\dfrac{1}{S^3(t)}\right) \left(\left( \sigma_1(1-u)\dfrac{S(t) I_s(t)}{N}\right)^2 &#x002B; (\sigma_2(C(t)-S(t)))^2\right) \\ &&#x002B; 2(S(t)&#x002B;C(t)&#x002B;I_s(t))\delta S(t)&#x002B;\dfrac{\rho}{S(t)} &#x002B;\left(1&#x002B;\dfrac{1}{C^3(t)}\right) (\sigma_2(C(t)-S(t)))^2\\ &&#x002B; 2(S(t)&#x002B;C(t)&#x002B;I_s(t))\beta \epsilon (1-u)\frac{S(t-\tau_1) I_s(t-\tau_1)}{N}\\ &&#x002B; \frac{\alpha}{I_s(t)}&#x002B;(1-\alpha)(\mu_s&#x002B;\eta_s)\frac{1}{I_s(t)}\\ &&#x002B; \left(1&#x002B;\dfrac{1}{I^3_s(t)}\right) \left(\left( -\sigma_3 I_s(t)(1-\mu_s-\eta_s\right) ^2&#x002B;\left(\sigma_1 \epsilon (1-u)\dfrac{S(t-\tau_1) I_s(t-\tau_1)}{N}\right) ^2\right). \end{align*}

By applying the elementary inequality 2aba2 + b2, we can easily increase the right-hand side of the previous inequality to obtain that LV(x)D(1+V(x)), \[ LV(x)\leq D(1&#x002B;V(x)), \]

where D is an adequate selected positive constant. By integrating both sides of the equality dV(x(t))=LV(x(t))dt+σ(x(t))d B t \[ dV(x(t))=LV(x(t))\textrm{d}t&#x002B;\sigma(x(t))\textrm{d}B_t \]

between t0 and tτk and acting the expectation, which eliminates the martingale part, we get that E(V(x(t τ k )) =E(V( x 0 ))+E t 0 t τ k L V( x s ))ds E(V( x 0 ))+E t 0 t τ k D (1+V( x s ))ds E(V( x 0 ))+DT+ t 0 t τ k E V( x s ))ds. \begin{align*} E(V(x(t\wedge \tau_k)) &= E(V(x_0))&#x002B;E\int^{t\wedge \tau_k}_{t_0}LV(x_s))\textrm{d}s\\ &\leq E(V(x_0))&#x002B;E\int^{t\wedge \tau_k}_{t_0}D(1&#x002B;V(x_s))\textrm{d}s\\ &\leq E(V(x_0))&#x002B;DT&#x002B;\int^{t \wedge \tau_k}_{t_0}EV(x_s))\textrm{d}s. \end{align*}

Gronwall’s inequality implies that E(V(x(t τ k ))(EV( x 0 )+DT)exp(CT). \[ E(V(x(t\wedge \tau_k))\leq (EV(x_0)&#x002B;DT)\exp(CT). \]

For ω ∈{τkT}, xi (τk) equals k or 1 k $\dfrac{1}{k}$ for some i = 1, 2, 3. Hence, V( x i ( τ k ))( k 2 + 1 k )( 1 k 2 +k ). \[ V(x_i(\tau_k))\geq \left(k^2&#x002B;\frac{1}{k}\right) \wedge \left(\frac{1}{k^2}&#x002B;k\right). \]

It follows that (EV( x 0 )+DT)exp(CT) E( χ { τ k T} (ω)V( x τ k ) ) ( k 2 + 1 k )( 1 k 2 +k )P( τ k T ). \begin{align*} (EV(x_0)&#x002B;DT)\exp(CT) &\geq E\left(\chi_{\{\tau_k \leq T\}}(\omega) V(x_{\tau_k})\right)\\ &\geq \left(k^2&#x002B;\frac{1}{k}\right) \wedge \left(\frac{1}{k^2}&#x002B;k\right)P\left(\tau_k \leq T\right). \end{align*}

Letting k, we get P(τeT) = 0. Since T is arbitrary, we obtain P(τe = ) = 1. With the same technique, we also deduce that the rest of the variables of the system are positive on [0, ). This concludes the proof.

3 Extinction of the disease

In this section, we obtain a sufficient condition for the extinction of the disease.

Theorem 3.1. Let (S(t), C(t), Is(t), Ia(t), Fb(t), Fg(t), Fc(t), R(t), M(t)) be a solution of the D-COVID-19 model (1.1) with positive initial value defined in (1.2). Assume that

σ 1 2 > β 2 2(α+(1α)( μ s + η s )) . \[ \sigma_1^2> \dfrac{\beta^2}{2(\alpha&#x002B;(1-\alpha)(\mu_s&#x002B;\eta_s))}. \]

Then, limsup t ln I s (t) t <0. \[ \limsup_{t\to \infty} \ln\dfrac{I_s(t)}{t}<0. \]

Namely, Is(t) tends to zero exponentially a.s., that is, the disease dies out with probability 1.

Proof. Let dln I s (t) = [ 1 2 I s 2 (t) ( ( σ 1 βɛ(1u)S(t τ 1 ) I s (t τ 1 ) N ) 2 + ( σ 3 ( μ s + η s 1) I s (t) ) 2 ) 1 I s (t) ( βɛ(1u)S(t τ 1 ) I s (t τ 1 ) N α I s (t)(1α)( μ s + η s ) I s (t) ) 1 2 I s 2 (t) ( ( σ 1 βɛ(1u)S(t τ 1 ) I s (t τ 1 ) N ) 2 + ( σ 3 ( μ s + η s 1) I s (t) ) 2 ) ]dt + 1 I s (t) σ 1 βɛ(1u)S(t τ 1 ) I s (t τ 1 ) N d B 1 + 1 I s (t) σ 3 ( μ s + η s 1) I s (t)d B 3 . \begin{align*} d\ln I_s(t) &= \left[ \vphantom{\dfrac{1}{2I_s^{2}(t)}\left(\left(\sigma_1 \dfrac{\beta \epsilon (1-u)S(t-\tau_1)I_s(t-\tau_1)}{N} \right)^{2} &#x002B; \left(\sigma_3(\mu_s&#x002B;\eta_s-1)I_s(t)\right)^{2}\right)}\dfrac{1}{I_s(t)}\left(\dfrac{\beta \epsilon (1-u) S(t-\tau_1)I_s(t-\tau_1)}{N} - \alpha I_s(t) - (1-\alpha)(\mu_s&#x002B;\eta_s)I_s(t)\right) \right.\\ & \quad\left. -\dfrac{1}{2I_s^{2}(t)}\left(\left(\sigma_1 \dfrac{\beta \epsilon (1-u)S(t-\tau_1)I_s(t-\tau_1)}{N} \right)^{2} &#x002B; \left(\sigma_3(\mu_s&#x002B;\eta_s-1)I_s(t)\right)^{2}\right)\right] \textrm{d}t\\ &\quad&#x002B; \dfrac{1}{I_s(t)} \sigma_1 \dfrac{\beta \epsilon (1-u)S(t-\tau_1) I_s(t-\tau_1)}{N} \textrm{d}B_1 &#x002B; \dfrac{1}{I_s(t)} \sigma_3 (\mu_s&#x002B;\eta_s-1)I_s(t) \textrm{d}B_3. \end{align*}

To simplify, we set G(t) := ɛ(1u)S(t τ 1 ) I s (t τ 1 ) N , R 1 (t) := σ 1 I s (t) βɛ(1u)S(t τ 1 ) I s (t τ 1 ) N = β σ 1 I s (t) G, R 3 (t) := σ 3 I s (t) ( μ s + η s 1) I s (t), H :=α(1α)( μ s + η s ). \begin{align*} G(t)&:= \dfrac{ \epsilon (1-u)S(t-\tau_1)I_s(t-\tau_1)}{N},\\ R_{1}(t)&:= \dfrac{\sigma_1}{I_s(t)} \dfrac{\beta \epsilon (1-u) S(t-\tau_1)I_s(t-\tau_1)}{N} =\dfrac{\beta \sigma_1}{I_s(t)} G,\\ R_{3}(t)&:= \dfrac{\sigma_3}{I_s(t)} (\mu_s&#x002B;\eta_s-1)I_s(t),\\ H &:=-\alpha -(1-\alpha)(\mu_s&#x002B;\eta_s). \end{align*}

We then get dln I s (t) = βG(t) I s (t) H(t) 1 2 ( ( σ 1 G(t) I s (t) ) 2 + ( σ 3 ( μ s + η s 1) I s (t) ) 2 )+ R 1 (t)d B 1 + R 3 (t)d B 3 = σ 1 2 2 [ ( G(t) I s (t) ) 2 2β σ 1 2 G(t) I s (t) ]+H+ R 1 (t)d B 1 + R 3 (t)d B 3 = σ 1 2 2 [ ( G(t) I s (t) β σ 1 2 ) 2 β 2 σ 1 4 ]+H+ R 1 (t)d B 1 + R 3 (t)d B 3 β 2 2 σ 1 2 +H+ R 1 (t)d B 1 + R 3 (t)d B 3 . \begin{align*} \textrm{d}\ln I_s(t) &= \dfrac{\beta G(t)}{I_s(t)} - H(t) -\dfrac{1}{2} \left(\left( \dfrac{\sigma_1 G(t)}{I_s(t)} \right)^{2}&#x002B; \left(\sigma_3(\mu_s &#x002B;\eta_s-1)I_s(t)\right)^{2}\right)&#x002B; R_1(t) \textrm{d}B_{1} &#x002B; R_3(t) \textrm{d}B_{3}\\ &= -\dfrac{\sigma_1^{2}}{2}\left[ \left( \dfrac{G(t)}{I_s(t)}\right) ^{2}- \dfrac{2\beta}{\sigma_1^{2}} \dfrac{G(t)}{I_s(t)}\right] &#x002B; H &#x002B; R_1(t) \textrm{d}B_{1}&#x002B; R_3(t) \textrm{d}B_{3}\\ &= -\dfrac{\sigma_1^{2}}{2} \left[ \left(\dfrac{G(t)}{I_s(t)} - \dfrac{\beta}{\sigma_1^{2}} \right) ^{2} - \dfrac{\beta^{2}}{\sigma_1^4} \right] &#x002B; H &#x002B; R_1(t) \textrm{d}B_{1} &#x002B; R_3(t) \textrm{d}B_{3}\\ &\leq \,\dfrac{\beta^{2}}{2\sigma_1^{2}} &#x002B; H &#x002B; R_1(t) \textrm{d}B_{1} &#x002B; R_3(t) \textrm{d}B_{3}. \end{align*}

Hence, ln I s (t) t ln I s (0) t + β 2 2 σ 1 2 +H+ M 1 (t) t + M 3 (t) t , \begin{align*} \dfrac{\ln I_s(t)}{t} &\leq \dfrac{\ln I_s(0)}{t} &#x002B; \dfrac{\beta^{2}}{2\sigma_1^{2}} &#x002B; H &#x002B; \dfrac{M_1(t)}{t} &#x002B; \dfrac{M_3(t)}{t}, \end{align*}

where M 1 (t)= 0 t R 1 (s)d B 1 , M 3 (t)= 0 t R 3 (s)d B 3 . \[ M_1(t)= \int_{0}^{t}R_{1}(s)\textrm{d}B_{1}, \quad M_3(t)= \int_{0}^{t}R_{3}(s)\textrm{d}B_{3}. \]

We have M 1 , M 1 t = 0 t ( 1 I s (s) σ 1 G(s) ) 2 ds = 0 t σ 1 2 ɛ 2 (1u) 2 S (t τ 1 ) 2 I s (t τ 1 ) 2 N 2 I s 2 ds 0 t σ 1 2 ɛ 2 (1u) 2 N 4 N 2 1 I s 2 ds 0 t σ 1 2 ɛ 2 (1u) 2 N 2 ds. \begin{align*} \left\langle M_1,M_1\right\rangle_t &= \int^t_0\left(\frac{1}{I_s(s)}\sigma_1G(s)\right)^2\textrm{d}s\\ &= \int^t_0{\sigma_1}^2\epsilon^2(1-u)^2\dfrac{S(t-\tau_1)^2I_s(t-\tau_1)^2}{N^2I_s^2}\textrm{d}s\\ &\leq \int^t_0{\sigma_1}^2\epsilon^2(1-u)^2\frac{N^4}{N^2}\frac{1}{I^2_s}\textrm{d}s\\ &\leq \int^t_0{\sigma_1}^2\epsilon^2(1-u)^2N^2\textrm{d}s. \end{align*}

Then, limsup t M 1 , M 1 t t σ 1 2 ɛ 2 (1u) 2 N 2 <. \[ \underset{t\rightarrow\infty}{\limsup}\dfrac{\langle M_1,M_1\rangle_t}{t} \leq {\sigma_1}^2\epsilon^2(1-u)^2N^2<\infty. \]

From the large number theorem for martingales [5], we deduce that lim t M 1 (t) t =0. \[ \underset{t\rightarrow\infty}{\lim}\frac{M_1(t)}{t}=0. \]

We also have M 2 , M 2 t = 0 t ( 1 I s (s) σ 3 ( μ s + η s 1) I s (s) ) 2 ds = 0 t ( σ 3 2 ( μ s + η s 1) 2 ) ds =( σ 3 2 ( μ s + η s 1) 2 )t. \begin{align*} \langle M_2,M_2\rangle_t &= \int^t_0\left(\frac{1}{I_s(s)}\sigma_3(\mu_s&#x002B;\eta_s-1)I_s(s)\right)^2\textrm{d}s\\ &= \int^t_0\left(\sigma^2_3(\mu_s&#x002B;\eta_s-1)^2\right)\textrm{d}s\\ &= \left(\sigma^2_3(\mu_s&#x002B;\eta_s-1)^2\right)t. \end{align*}

Thus, limsup t M 2 , M 2 t t σ 3 2 ( μ s + η s 1)<. \[ \underset{t\rightarrow\infty}{\limsup}\dfrac{\langle M_2,M_2\rangle_t}{t} \leq \sigma^2_3(\mu_s&#x002B;\eta_s-1)<\infty. \]

We deduce that lim t M 2 (t) t =0. \[ \underset{t\rightarrow\infty}{\lim}\frac{M_2(t)}{t}=0. \]

Subsequently, limsup t ln I s (t) t β 2 2 σ 1 2 α(1α)( μ s + η s ). \[ \underset{t\rightarrow\infty}{\limsup}\dfrac{\ln I_s(t)}{t} \leq \frac{\beta^2}{2\sigma_1^2}-\alpha-(1-\alpha)(\mu_s&#x002B;\eta_s). \]

We conclude that if β 2 2 σ 1 2 α(1α)( μ s + η s )<0 $\dfrac{\beta^2}{2\sigma_1^2} -\alpha-(1-\alpha)(\mu_s&#x002B;\eta_s)<0$, then limI(t) t =0 $\underset{t\rightarrow\infty}{\lim I(t)}=0$. The proof is complete.

4 Results and discussion

In this section, we simulate the forecasts of the D-COVID-19 model (1.1), relating the deconfinement strategy adopted by Moroccan authorities with two scenarios. We assume u defined as follows: u={ u 0 , on[March 2,March10]; u 1 , on(March 10,March16]; u 2 , on(March 16,March20]; u 3 , on(March 20,April6]; u 4 , on(April 6,April25]; u 5 , fromApril 25 on; \[ u=\left\{ \begin{array}{ll} u_{0}, & \hbox{$ \text{on}\, [\text{March 2},\text{March 10]}$;} \\ u_{1}, & \hbox{$ \text{on}\, (\text{March 10},\text{March 16]}$;}\\ u_{2}, & \hbox{$ \text{on}\, (\text{March 16},\text{March 20]}$;}\\ u_{3}, & \hbox{$ \text{on}\, (\text{March 20},\text{April 6]}$;}\\ u_{4}, & \hbox{$ \text{on}\, (\text{April 6},\text{April 25]}$;}\\ u_{5}, & \hbox{$ \text{from}\, \text{April 25} \text{ on}$;}\\ \end{array} \right. \]

where ui ∈ (0, 1], for i = 0, 1, 2, 3, 4, 5, measures the effectiveness of applying the multiple preventive interventions imposed by the authorities and presented in Table 1.

COVID-19 is known as a highly contagious disease and its transmission rate, β, varies from country to country, according to the density of the country and movements of its population. Ozair et al. [17] assumed β to be [0.198 − 0.594] per day for Romania, and [0.097 − 0.291] per day for Pakistan. Further, Kuniya [10] estimated β as 0.26 (95%CI, 2.4 − 2.8). Observing the number of daily reported cases of COVID-19 in Morocco, we estimate β as 0.4517 (95%CI, 0.4484 − 0.455). After the infection, the patient remains in a latent period for 5.5 days [1, 20], in average, before becoming symptomatic and infectious or asymptomatic with a percentage that varies from 20.6% of infected population to 39.9% [15], while the time needed before his hospitalization is estimated to be 7.5 days [6, 8, 19]. All the parameter values chosen for the D-COVID-19 model (1.1) are summarized in Table 2.

Remark 4.1. From a biological point of view, the latency period is independent of the region or country under study, depending only on the structural nature of the SARS-CoV-2 coronavirus.

We consider that all measures and the adopted confinement strategy previously discussed are conserved. The evolution on the number of diagnosed infected positive individuals given by the D-COVID-19 model (1.1) versus the daily reported confirmed cases of COVID-19 in Morocco, from March 2 to May 6, is presented in Figure 2. We see that the curve generated by the D-COVID-19 model (1.1) follows the trend of the daily reported cases in Morocco. So, we confirm that the implemented measures taken by the authorities have an explicit impact on the propagation of the virus in the population since the curve of the D-COVID-19 model (1.1) has been flattening from April 17 and tends to go towards the extinction of the disease from May 05. In Figure 3, we see that Morocco has spent almost 40% of the total duration of the epidemic at May 11 and will reach extinction after four months, in average, from the start of the epidemic on March 2, 2020 (t = 0).

To prove the biological importance of delay parameters, we give the graphical results of Figure 4, which allow to compare the evolution of diagnosed positive cases with and without delays. We observe in Figure 4 a high impact of delays on the number of diagnosed positive cases. Indeed, the plot of model (1.1) without delays (τi = 0,  i = 1, 2, 3, 4) is very far from the clinical data. Thus, we conclude that delays play an important role in the study of the dynamical behavior of COVID-19 worldwide, especially in Morocco, and allows to better understand the reality.

In Figures 57, we consider the deconfinement of 30% of the population returning to work from May 20, and this proportion is immediately integrated into the susceptible population. Numerical simulations are presented for three possible scenarios. In the first, we consider that the whole population highly respects the majority of the measures announced by the authorities in relation with the deconfinement (Fig. 5). The second and third scenarios show the direct impact on the curves when the population moderately respects the measures with different levels, σ2 = 0.10 and σ2 = 0.15, respectively.With the last two scenarios we observe the growth in the final number of infected, deaths, severe and critical forms, which are the most important to monitor, since the health system should not be saturated. It is also important to note the appearance of a second significant peak and the fact that the time required for extinction becomes longer, which relates to the value of σ2 (Figs. 6 and 7).

It should be mentioned that all our plots were obtained using the Matlab numerical computing environment by discretizing system (1.1) by means of the higher order method of Milstein presented in [7] and used in [11].

Table 1

Summary of considered non-pharmaceutical interventions.

Table 2

Parameter values for the D-COVID-19 model (1.1).

thumbnail Figure 2

Evolution of COVID-19 confirmed cases in Morocco per day: curve predicted by our model (1.1) accordigly with Tables 1 and 2 versus real data.

thumbnail Figure 3

Evolution given by the D-COVID-19 model (1.1) without deconfinement (ρ = 0).

thumbnail Figure 4

Effect of delays on the diagnosed confirmed cases versus clinical data.

thumbnail Figure 5

Evolution of the D-COVID-19 model (1.1) with deconfinement (ρ = 0.3) from May 20, 2020 and high effectiveness of the measures (σ2 = 0.01).

5 Conclusion

In this work, we have proposed a delayed stochastic mathematical model to describe the dynamical spreading of COVID-19 in Morocco by considering all measures designed by authorities, such as confinement and deconfinement policies. More precisely, our model takes into account four types of delays: the first one is related to the incubation period, the second is the time needed to move from the symptomatic infected individuals to the three forms of diagnosed cases, the third is the time needed to move from the class of infected individuals to the recovered or dead class, while the last one is the time needed to pass from the three types of classes of individuals supported by the Moroccan health system, and under quarantine, to the recovered or dead compartments. Besides, to well describe reality, we have added a stochastic factor resulting from possible maladjustment of the population individuals to the measures.

To show that our model is mathematically and biologically well-posed, we have proved the global existence of a unique positive solution (see Thm. 2.1). Our result has shown a possible extinction of the disease when σ 1 2 $\sigma_1^{2}$ is greater thana threshold parameter (see Thm. 3.1).

In addition, numerical simulations have been performed to forecast the evolution of COVID-19. More precisely, we have shown that the evolution of our D-COVID-19 model follows the tendency of daily reported confirmed cases in Morocco (see Fig. 2). Further, if Moroccan people would maintain, strictly, their confinement policy, we observe that the disease dies out around four months from March 2, 2020 (see Fig. 3). On the other hand, in response to the decision of deconfinement represented by the liberation of the 30% of population, which took place at May 20, we simulate three scenarios corresponding to different values of the intensity σ2. When σ2 = 0.01 (Fig. 5), the eradication of the disease from the population comes early compared to the cases when σ2 = 0.10 (Fig. 6) and σ2 = 0.15 (Fig. 7). Additionally, the number of diagnosed confirmed cases mainly changes because of the value of this intensity and a small perturbation leads to relevant quantitative changes and significant variations on the time needed for extinction. Thus, we observe that the value of this perturbation has a high impact on the evolution of COVID-19, which means the Moroccan population has a big interest to respect the governmental measures announced May 20, 2020, in order to have a successful and good deconfinement strategy.

Here we have compared the predictions of the proposed model (1.1) with real data until middle of May 2020. We leave the comparison of the real data in Morocco till the end of 2020 to a future work, where we also plan to incorporate the predictions of the evolution of our COVID-19 model with respect to preventive Moroccan measures by regions and cities.

thumbnail Figure 6

Evolution of the D-COVID-19 model (1.1) with deconfinement (ρ = 0.3) from May 20, 2020 and moderate effectiveness of the measures (σ2 = 0.10).

thumbnail Figure 7

Evolution of the D-COVID-19 model (1.1) with deconfinement (ρ = 0.3) from May 20, 2020 and moderate effectiveness of the measures (σ2 = 0.15).

Acknowledgements

The authors are grateful to an Associate Editor and two anonymous referees, who kindly reviewed an earlier version of the manuscript and provided several valuable suggestions and comments.

References

  1. S.G. Baum, COVID-19 incubation period: An update, Journal Watch. Available from https://www.jwatch.org/na51083/2020/03/13/covid-19incubation-period-update (2020). [Google Scholar]
  2. M.S. Boudrioua and A. Boudrioua, Predicting the COVID-19 epidemic in Algeria using the SIR model. Preprint medRxiv 20079467v1 (2020). [Google Scholar]
  3. D. Fanelli and F. Piazza, Analysis and forecast of COVID-19 spreading in China, Italy and France. Chaos Solitons Fract. 134 (2020) 109761. [CrossRef] [Google Scholar]
  4. A.E. Gorbalenya, S.C. Baker, R.S. Baric, R.J. de Groot, C. Drosten, A.A. Gulyaeva, B.L. Haagmans, C. Lauber, A.M. Leontovich, B.W. Neuman, D. Penzar, S. Perlman, L.L.M. Poon, D.V. Samborskiy, I.A. Sidorov, I. Sola and J. Ziebuhr, The species severe acute respiratory syndrome related coronavirus: classifying 2019-nCoV and naming it SARS-CoV-2. Nat. Microbiol. (2020) 1–9. [PubMed] [Google Scholar]
  5. A. Grai, D. Greenhalgh, L. Hu, X. Mao and J. Pan, A stochastic differential equations SIS epidemic model. SIAM J. Appl. Math. 71 (2011) 876–902. [Google Scholar]
  6. Haut Conseil de la Santé Publique, Avis relatif aux recommandations thérapeutiques dans la prise en charge du COVID-19 (complémentaire A l’avis du 5mars 2020). Available from: https://sfar.org/avis-relatif-aux-recommandations-therapeutiques-dans-la-prise-en-charge-du-covid-19/ (2020). [Google Scholar]
  7. D.J. Higham, An algorithmic introduction to numerical simulation of stochastic differential equations. SIAM Rev. 43 (2001) 525–546. [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
  8. C. Huang, Y. Wang, X. Li, L. Ren, J. Zhao, Y. Hu, L. Zhang, G. Fan, J. Xu, X. Gu, Z. Cheng, T. Yu, J. Xia, Y. Wei, W. Wu, X. Xie, W. Yin, H. Li, M. Liu, Y. Xiao, H. Gao, L. Guo, J. Xie, G. Wang, R. Jiang, Z. Gao, Q. Jin, J. Wang and B. Cao, Clinical features of patients infectedwith 2019 novel coronavirus in Wuhan, China. Lancet 395 (2020) 497–506. [Google Scholar]
  9. B. Ivorra, M.R. Ferrandez, M. Vela-Perez and A.M. Ramos, Mathematical modeling of the spread of the coronavirusdisease 2019 (COVID-19) taking into account the undetected infections. The case of China. Commun. Nonlinear Sci. Numer. Simul. 88 (2020) 105303. [CrossRef] [PubMed] [Google Scholar]
  10. T. Kuniya, Prediction of the epidemic peak of coronavirus disease in Japan, 2020. J. Clin. Med. 9 (2020) 789. [Google Scholar]
  11. M. Mahrouf, K. Hattaf and N. Yousfi, Dynamics of a stochastic viral infection model with immune response. MMNP 12 (2017) 15–32. [Google Scholar]
  12. X. Mao, Stochastic differential equations and applications, second edition. Horwood Publishing Limited, Chichester, (2008). [CrossRef] [Google Scholar]
  13. Ministry of Health of Morocco, Department of Epidemiology and Disease Control. Available from: http://www.sante.gov.ma/Pages/Accueil.aspx (2020). [Google Scholar]
  14. Ministry of Health of Morocco, The official portal of Corona virus in Morocco. Available from: http://www.covidmaroc.ma/pages/Accueil.aspx (2020). [Google Scholar]
  15. K. Mizumoto, K. Kagaya, A. Zarebski and G. Chowell, Estimating the asymptomatic proportion of coronavirusdisease 2019 (COVID-19) cases on board the Diamond Princess cruise ship, Yokohama, Japan, 2020. Euro Surveill. 25 (2020) 2000180. [Google Scholar]
  16. F. Ndaïrou, I. Area, J.J. Nieto and D.F.M. Torres, Mathematical modeling of COVID-19 transmission dynamics with a case study of Wuhan. Chaos Solitons Fract. 135 (2020) 109846. [CrossRef] [Google Scholar]
  17. M. Ozair, T. Hussain, M. Hussain, A.U. Awan and D. Baleanu, Estimation of transmission potential and severity of COVID-19 in Romania and Pakistan. Preprint medRxiv 20088989v1 (2020). [Google Scholar]
  18. B. Tang, X. Wang, Q. Li, N. L. Bragazzi, S. Tang, Y. Xiao and J. Wu, Estimation of the transmission risk of the 2019-nCoV and its implication for public health interventions. J. Clin. Med. 9 (2020) 462. [Google Scholar]
  19. D. Wang, B. Hu, C. Hu, F. Zhu, X. Liu, J. Zhang, B. Wang, H. Xiang, Z. Cheng, Y. Xiong, Y. Zhao, Y. Li, X. Wang and Z. Peng, Clinical characteristics of 138 hospitalized patientswith 2019 novel coronavirus-infected pneumonia in Wuhan, China. JAMA 323 (2020) 1061–1069. [Google Scholar]
  20. WHO, Coronavirus disease 2019 (COVID-19) Situation Report 73. Available from: https://apps.who.int/iris/handle/10665/331686 (2020). [Google Scholar]
  21. J.T. Wu, K. Leung and G.M. Leung, Nowcasting and forecasting the potential domestic and international spread of the 2019-nCoV outbreak originating in Wuhan, China: a modelling study. The Lancet 395 (2020) 689–697. [CrossRef] [PubMed] [Google Scholar]

All Tables

Table 1

Summary of considered non-pharmaceutical interventions.

Table 2

Parameter values for the D-COVID-19 model (1.1).

All Figures

thumbnail Figure 1

Schematic diagram of model (1.1).

In the text
thumbnail Figure 2

Evolution of COVID-19 confirmed cases in Morocco per day: curve predicted by our model (1.1) accordigly with Tables 1 and 2 versus real data.

In the text
thumbnail Figure 3

Evolution given by the D-COVID-19 model (1.1) without deconfinement (ρ = 0).

In the text
thumbnail Figure 4

Effect of delays on the diagnosed confirmed cases versus clinical data.

In the text
thumbnail Figure 5

Evolution of the D-COVID-19 model (1.1) with deconfinement (ρ = 0.3) from May 20, 2020 and high effectiveness of the measures (σ2 = 0.01).

In the text
thumbnail Figure 6

Evolution of the D-COVID-19 model (1.1) with deconfinement (ρ = 0.3) from May 20, 2020 and moderate effectiveness of the measures (σ2 = 0.10).

In the text
thumbnail Figure 7

Evolution of the D-COVID-19 model (1.1) with deconfinement (ρ = 0.3) from May 20, 2020 and moderate effectiveness of the measures (σ2 = 0.15).

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.