key: cord-0681815-hs4vs2cu authors: Fredj, Haifa Ben; Chrif, Farouk title: Novel Corona virus Disease infection in Tunisia: Mathematical model and the impact of the quarantine strategy date: 2020-06-10 journal: Chaos Solitons Fractals DOI: 10.1016/j.chaos.2020.109969 sha: 66e28b5f2c84536bd0f3695c3a0dc3f88dd354ee doc_id: 681815 cord_uid: hs4vs2cu Corona virus disease (COVID-19) is an extremely serious infection with an extremely high death rate worldwide. In March, the disease was declared a ”global pandemic” by the World Health Organization (WHO). Until now, there is no known vaccine or drug, since the unknown things related to the disease are more important than our theoretical and empirical knowledge. However, mathematical modeling and the estimation of the basic number of reproductions can provide clarifications in order to determine the potential and severity of this epidemic and therefore provide essential information to identify the type of measures and interventions to be taken to control the intensity of the spread of the disease. Hence, in this paper, we propose a new deterministic compartmental model based on the clinical progression of the disease, the epidemiological state of the individuals and the intervention for the dynamics of COVID-19 infections. Our approach consists of seven phenotypes: the susceptible humans, exposed humans, infectious humans, the recovered humans, the quarantine population, there recovered-exposed and deceased population. We proved first through mathematical approach the positivity, boundness and existence of a solution to the considered model. We also studied the existence of the disease free equilibrium and corresponding stability. Our work shows, in particular, that the disease will decrease if the number of reproduction R(0) was less than one. Moreover, the impact of the quarantine strategies to reduce the spread of this disease is discussed. The theoretical results are validated by some numerical simulations of the system of the epidemic’s differential equations. It should be mentioned that, the error between the considered model and the official data curve is quite small. Since the beginning of the COVID-19 epidemic in Wuhan City on December 2019, the opinions of scientists, researchers and commentators contradict each other every day. On 7 January, the coronavirus disease (COVID-2019) which was named as a severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) by International Committee on Taxonomy of Viruses on 11 February, 2020, was identified as the causative virus by Chinese authorities [17] , and has become a pandemic especially by travelers [16] . This forced the World Health Organization (WHO) to consider the dramatic spread of the infection in March 2020 as a public health emergency of international concern. This epidemic is characterized by its rapid spread and its symptoms do not appear quickly. In particular, according to the WHO, the incubation period is from 2 days to 14 days [17] . In addition, there is no anti-viral treatment or vaccination officially approved for the management or prevention of this epidemic. In order to fight this outbreak, the public health decision and policy makers should decide and follow strategic and health-care management. Since the declaration of the first case of COVID-19 in Tunisia in early March 2020, thousands of screenings have been carried out in Tunisia. Following the first recommendations of the WHO Tunisia, as France, Italy or Spain, first opted for rare and targeted screenings on suspicious people, their entourage and certain sources of contamination; in particular for people arriving from abroad. Next, If a person is positive, the authorities try to trace all the people with whom he has been in contact and they are called to place themselves in self-containment without necessarily being tested. They can only be tested if they themselves have symptoms of COVID-19. It is clear that the low number of screenings, mainly on "suspect" cases or those presenting significant symptoms, does not give a precise idea of the number of people who could potentially be infected without knowing them. This gap between the day of infection and the day of diagnosis can have serious consequences on the spread of the epidemic. All this shows the complexity of the situation since the unknown things related to the epidemic are more important than the things we already know. These are the daily question that must be answered: how many people exactly recover from COVID-19? How many people are infected from COVID-19? How many people died from COVID-19? Recently, several mathematical models have been published in order to be able to study the dynamics and the evolution of this pandemic. One can refer to [1, 4, 6, 7, 8, 9, 13, 14, 18] and their references. It appears from clinical experiences and recent articles that knowing the data of infected people in the population would be very useful to have better models of when disease will peak and decline, and also when we can begin to let people go back to work. Also, knowing the real number of recovered people also could indicate how easily people can build immunity against the virus. Motivated by the above discussion our method comprises seven phenotypes; the susceptible humans, exposed humans, infectious humans, the recovered humans, the quarantine population, the recovered-exposed and the dead population in order to improve and adapt the susceptible-infected-recovered (SIR) model. In fact, we noticed that people previously infected and recovered from COVID-19 are generally excluded as susceptible individuals in the modeling, which would have an impact on predicting the number of cases that will occur in the near future. In our approach we chose to integrate them into the model as a recovered-exposed population and which we will note E r . In addition, at the individual level, if people can find out if they have been or slightly infected and cured, and if they show civility and respect the distancing, they can safely return to work once the general quarantine is off. Our main contributions in this paper are: -We gave a model of dynamical behavior of COVID-19 in Tunisia. -We studied the qualitative properties of our model. The reminder of the paper is organized as follows. In Section 2, we present our model. Section 3 is devoted to the mathematical analysis: Boundness, positivity, and the equilibrium point. the dynamics of exposed and infected population are also studied in this section. In Section 4, via Matab, we perform numerical simulations of three types of populations. Section 5 concludes the paper with some recommendations. The classical susceptible-infected-recovered (SIR) model in epidemiology [5] allows the determination of critical condition of disease development in the population irrespective of the total population size over a short period of time. The SIR is considered from the following simplest ODE system: For the sizes of the susceptible sub-population S, infected I, and recovered R. The term λ N IS describes the disease transmission rate due to the contacts between susceptible and infected individuals, βI characterizes the rate of recovery of infected people. It is assumed here that the recovered individuals do not return to susceptible class, that is, recovered individuals have immunity against the disease; they cannot become infected again and cannot infect susceptible either. Day by day, we see that the specifics of this virus require more complex models. Roughly speaking, it is necessary to question the duration of the incubation, or the presence of people who do not have symptoms after their infection, but who nevertheless participate in the spread of the virus and consider compartments for these different cases. In our work, we have used a 7-Phenotype model which takes into account people who are infected, exposed, recovered, those susceptible, the deceased, Quarantine but also recovered-exposed people who are not counted by the government. The interaction between the above sub-populations can be described by the compartmental diagram in Figure Obviously, the system (2.2) can be written as follows where with the initial conditions satisfying the following inequalities Given the non-negative initial conditions (3.2), then the solutions S(t), E(t), I(t), E r (t), R(t), D(t) and Q(t) are non-negative for all t ≥ 0 and bounded. Proof. First, let us prove that the solution of the system (3.1) is positive. Then, we have From the second case; S and E are not both negative as the same time. 3 rd case: Suppose that ∃t 1 > 0 such that 4 th case: Suppose that ∃t 1 > 0 such that We can write E as the following form witch is impossible. We can write E as the following form which is impossible. Then, we get E, I and S as positive functions. Therefore, we get E r , R, D and Q as increasing functions. Thus, we have E r , R, D and Q as positive functions. Now, we prove that the solution is bounded. From the positivity of the solution, we have the following inequality Then, by integrating We can deduce that S(t) ≤ S(0) and E(t) ≤ S(0) + E(0), for all positive real t. By the same way, we get E(t) + S(t) + I(t) ≤ (S(0) + E(0) + I(0))e −m 2 t ; where m 2 = min{λ 3 , α 2 , β 2 + β 2 }. Then, I(t) ≤ S(0) + E(0) + I(0) for t ≥ 0. 3.2 Disease free equilibrium, Reproduction number R 0 and extinction of infected population 8 For E r ,R, D and Q, using the recent inequalities, we have From integrating Thus, the proof is given. Firstly, the function X ∈ R 7 → AX + f (X) is continuous. Pose X = (S, E, I, E r , R, Q, D) and Y = (S , E , I , E r , R , Q , D ). From the bounded solution, we get Then, the Cauchy-Lipschitz is satisfied, and consquently, the model (3.1) possesses a unique solution [2] . At the disease-free state, there is no disease in the human population which implies E = I = R = D = Q = E r = 0. Thus, the disease-free equilibrium of the model (3.1) is given by N, 0, 0, 0, 0, 0, 0) . We can also verify whether or not the state-point is stable. In our case, according the positivity of our model, we have from the first equation of system (2.2) As a result, the state-point is unstable. Now, let us calculate the basic reproduction number R 0 which is the average number of secondary infections caused by an infectious individual during his or her entire period of infectiousness (Diekmann et al) [3] . The basic reproduction number is an important non-dimensional quantity in epidemiology as it sets the threshold in the study of a disease both for predicting its outbreak and for evaluating its control strategies. Using the next generation operator approach by van den Driessche and Watmough [15] , we have And we have the corresponding effective control reproduction number defined as in [3] R 0 (t) = This quantity provides us with a clear index to evaluate the control strategy for any time t. Theorem 3.3. If R 0 < 1 , then exposed and infected population will extinct. Proof. Let us define the Lyapunov functional The derivative of V is given by Thus, where V (0) = λ 2 I(0) + (β 1 + β 2 )E(0). Therefore, for R 0 < 1 we have This shows that the disease will be extinct if R 0 < 1. The effect of quarantine strategies 10 As most countries of the world, Tunisia fights vigorously against the spread of the COVID-19 epidemic with some disabilities such as the lack of medical equipment. Therefore, in order to brake the spread of this emerging epidemic the general quarantine is adopted. In addition, since the number of beds in intensive care in the hospital is reduced and insufficient, the infected population is not fully hospitalized except in cases with fairly severe symptoms. Besides, the rest of the population should respect total isolation. The natural question here is: What is the impact of this decision in other words this quarantine strategy? By equation one in the system (2.2), the susceptible population can be evaluate from the following expression It is clear that the parameter λ 3 reduces the number of people exposed to infection with the disease. From the evolution of exponential function, we can deduce that when the rate λ 3 is getting bigger, the susceptible population is getting smaller. When the λ 3 = 0, we have this diagram and the corresponding non-linear differential equations where the parameters were defined in the Table 1 . Clearly, by the first equation, one has which implies that the disease almost infected every susceptible person. This theoretical result which justifies the quarantine strategy will be consolidated by a numerical study in the following paragraph. Using data of Tunisia from March 14 th to April 8 th 2020, we can estimate the parameters'values. The results can be summarized as follows: parameters values And for stimulation, we chose the following initial condition: S(0) = N-E(0)-I(0), E(0) =200, I(0) =18, E r (0) = R(0) = D(0) = 0. In addition to the population plots in Figures 4.1 and 4 .2, we collected some meaningful quantitative information about the model parameters (Table 2 ) and the peak values for infected and exposed populations ( Table 3) . The results can be summarized as follows : Peak infected Table 3 : Infected and exposed peak values in Tunisia region. And we have R 0 < 1, therefore the disease will disappear. The human body needs around 10 days to get rid of the organism. However, cases where the body needed 21 days to get rid of symptoms are detected. So after the peak, it is logical that the rate of recovered population increases. Then, we can get the following Figures To make a better illustration of quarantine strategy, we tested different home quarantine rates (λ 3 ) in Figures (4.7.a-h) . The change in the λ 3 parameter generates a consequent variation in the rate of the infected population and the deceased persons (Figures 4.8.a-b) . The most spectacular is the case λ 3 = 0 which corresponds to Q = 0; or in other words, in the absence of a quarantine strategy. We can see that on the 20th day of the epidemic we get 3, 75x10 5 infected people (Figure 4.9 .a). This number is far above the capacity of health structures in Tunisia. Consequently, the choice of containment strategy has significantly reduced the number of deceased and infected persons. Based on the COVID-19 cases data collected from March 14 th to April 8 th , we estimated the basic reproduction numberR 0 and the prediction of the epidemic's evolution. According to the comparison between the official epidemic data based on our model (Figure 4 .1) one can see that the simulation curve model of the number of confirmed Infections I(t) and the number of dead people D(t) match the official data curve very well during the period of 30 days. On the one side, the case of this epidemic in Tunisia is not critical: the peak of infected cases less than 700 which was very comfortable than many Arabian countries (Emirate, Qatar, Morocco, Algeria, ...) and also Western countries (Portugal, France, Spania...). On the other side, the number of infected and dead due to this disease remains relatively low in Tunisia compared to what is happening elsewhere in the world. And the most remarkable thing is that this situation is better than our neighbors (Algeria, France, Italy). Health policies are probably there for something but is it enough to have a plausible theory? Based on the data at the start of the disease in March, we can deduce from the theoretical study and the numerical simulations that the modification of the parameter λ 3 generates a consequent variation in the rate of the infected population and deceased persons (Figures 4.8.a-b) . Naturally, the most spectacular case is the case λ 3 = 0 which corresponds to Q=0. In other words, in the absence of a quarantine strategy, we can see that on the 20th day of the epidemic, we get 3,75x10 5 infected people (Figure 4.9 .a). This number is much higher than the current capacity of health structures in Tunisia. Consequently, the choice of the containment strategy considerably reduced the number of dead and infected people. With the absence of drug and vaccination, we found that home quarantine can really limit the spread of the epidemic. However, this strategy cannot last for an infinite time, due to the economic consequences on the country. The indicator of infected cases is not the problem. The exposed cases is the most important indicator. We predicted the initial exposed population to be about 200. And after 50 days in Figure 4 .2, the indicator of infected cases shows less infected people. Accordingly, the government is considering the relaxation of the quarantine strategy. The proposed model is flexible and we believe it could be easily adapted to monitor various infected are as with different restriction policies. In this article, we presented a new model based on nonlinear differential equations allowing to model this COVID-19 in Tunisia. The aim of our work, at first, is to provide initial ideas and guidelines for a quantitative and qualitative study of our considered model. In particular, the positivity, boundness and the existence of a solution are established. The impact of the quarantine strategy is also investigated and discussed by two methods. It thus appears that in addition to a rigorous application of quarantine, accurately nearcasting the epidemic trend and projecting the peak time require real-time data and the knowledge about the implementation resources in order to construct the adequate policy and effective decisions of major public health interventions. We hope our work motivates new research to produce more elaborate and precise methods related to this model or to give significant improvements. Modeling and forecasting of epidemic spreading: The case of Covid-19 and beyond A, An introduction To Ordinary differential Equations On the definition and Computation of the basic reproduction ratio R 0 in the model of infectious disease in Heterogeneous populations Analysis and forecast of COVID-19 spreading in China The Mathematics of Infectious Diseases Modeling the Control of COVID-19: Impact of Policy Interventions and Meteorological Factors Early dynamics of transmission and control of COVID-19: a mathematical modelling study Lancet Infect Dis (2020) published online Mathematical modeling of COVID-19 transmission dynamics with a case study of Wuhan The effect of control strategies that reduce social mixing on outcomes of the COVID-19 epidemic in Wuhan, China, The Lancet Public Health Available online 25 An updated estimation of the risk of transmission of the novel coronavirus (2019-nCov) The effectiveness of quarantine and isolation determine the trend of the COVID-19 epidemics in the final phase of the current outbreak in China Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission Travellers give wings to novel coronavirus (2019-nCoV) Predicting turning point, duration and attack rate of COVID-19 outbreaks in major Western countries The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.