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 Keywords and phrases. Coronavirus disease 2019 (COVID- 19) , deconfinement strategy, mathematical modeling, delayed stochastic differential equations (DSDEs), extinction.
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:
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; I s 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; I a 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; F b , F g and F c 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
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 F b , F g and F c , by the rates γ b , γ g and γ c , respectively. Then, γ b + γ g + γ c = 1. The mean recovery period of these forms are denoted by 1/r b , 1/r g and 1/r c , 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 I s and the three forms F b , F g and F c , respectively. Here, B 1 (t), B 2 (t) and B 3 (t) are independent standard Brownian motions defined on a complete probability space (Ω, F , P) with a filtration {F t } t≥0 and satisfying the usual conditions, that is, they are increasing and right continuous while F 0 contains all P-null sets and σ i represents the intensity of B i , 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) by perturbing the most sensitive parameters: β, α, δ, and ρ.
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 I s (t − τ 2 ) terms are the same as σ 3 dB 3 (t) although they are premultiplied by different constants. Indeed, the portion of diagnosed symptomatic infected population is completely distributed into the three forms F b , F g , and F c , by the rates γ b , γ g and γ c , respectively. Then,
In addition, we assume that the parameter α, which measures the efficiency of public health administration for hospitalization, undergoes random fluctuations. Remark 1.5. Temporarily asymptomatic individuals are included in the class I s of symptomatic, while individuals in I a , 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) satisfy:
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). A sufficient condition for the extinction is established 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.
Let us denote R 9
We begin by proving the following result.
Theorem 1. For any initial value satisfying condition (2), there is a unique solution
to the D-COVID-19 model (1) that remains in R 9 + with probability one. Proof. Since the coefficients of the Stochastic Differential Equations with several delays (1) are locally Lipschitz continuous, it follows from [12] that for any square integrable initial value x(0) ∈ R 9 + , 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 k 0 > 0 be sufficiently large for 1 k 0 < x(0) < k 0 . For each integer k ≥ k 0 , we define the stopping time τ k := inf t ∈ [0, τ e )/x i (t) / ∈ 1 k , k for some i = 1, 2, 3 , where inf ∅ = ∞. It is evident that τ k ≤ τ e . Let T > 0, and define the twice differentiable function V on R 3 + → R + as follows:
By Itô's formula, for any 0 ≤ t ≤ τ k ∧ T and k ≥ 1 we have
where L is the differential operator of function V :
Thus,
By applying the elementary inequality 2ab ≤ a 2 + b 2 , we can easily increase the right-hand side of the previous inequality to obtain that
where D is an adequate selected positive constant. By integrating both sides of the equality dV (x(t)) = LV (x(t))dt + σ(x(t))dB t between t 0 and t ∧ τ k and acting the expectation, which eliminates the martingale part, we get that
EV (x s ))ds.
Gronwall's inequality implies that
For ω ∈ {τ k ≤ T }, x i (τ k ) equals k or 1 k for some i = 1, 2, 3. Hence,
It follows that
Letting k → ∞, we get P (τ e ≤ T ) = 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.
In this section, we obtain a sufficient condition for the extinction of the disease.
) be a solution of the D-COVID-19 model (1) with positive initial value defined in (2) . Assume that
.
Then,
Namely, I s (t) tends to zero exponentially a.s., that is, the disease dies out with probability 1.
To simplify, we set
We then get
Hence,
where
We have
Then,
From the large number theorem for martingales [5] , we deduce that
We also have
Thus,
We conclude that if
The proof is complete.
In this section, we simulate the forecasts of the D-COVID-19 model (1), relating the deconfinement strategy adopted by Moroccan authorities with two scenarios. We assume u defined as follows: where u i ∈ (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 . [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) are summarized in Table 2 . 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) 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) 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) 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) 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 5, 6 , and 7, 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 (Figure 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 (Figures 6 and 7) . It should be mentioned that all our plots were obtained using the Matlab numerical computing environment by discretizing system (1) by means of the higher order method of Milstein presented in [7] and used in [11] .
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 Theorem 1) . Our result has shown a possible extinction of the disease when σ 2 1 is greater than a threshold parameter (see Theorem 2) . 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 Figure 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 Figure 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 ( Figure 5 ), the eradication of the disease from the population comes early compared to the cases when σ 2 = 0.10 ( Figure 6 ) and σ 2 = 0.15 ( Figure 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) 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.