SEIR model for COVID-19 dynamics incorporating the environment and social distancing

Coronavirus disease 2019 (COVID-19) is a pandemic respiratory illness spreading from person-to-person caused by a novel coronavirus and poses a serious public health risk. The goal of this study was to apply a modified susceptible-exposed-infectious-recovered (SEIR) compartmental mathematical model for prediction of COVID-19 epidemic dynamics incorporating pathogen in the environment and interventions. The next generation matrix approach was used to determine the basic reproduction number R0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R_0$$\end{document}. The model equations are solved numerically using fourth and fifth order Runge–Kutta methods. We found an R0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R_0$$\end{document} of 2.03, implying that the pandemic will persist in the human population in the absence of strong control measures. Results after simulating various scenarios indicate that disregarding social distancing and hygiene measures can have devastating effects on the human population. The model shows that quarantine of contacts and isolation of cases can help halt the spread on novel coronavirus.


Introduction
Coronaviruses are a large family of viruses that are known to cause illness ranging from the common cold to more severe diseases such as severe acute respiratory syndrome (SARS). The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) was identified as the cause of a cluster of pneumonia cases in Wuhan [1], a city in the Hubei Province of China, at the end of 2019. It subsequently spread throughout China and elsewhere, becoming a global health emergency. In February 2020, the World Health Organization (WHO) designated the disease coronavirus disease 2019 (COVID-19) a global pandemic [2].
The objective of this study was to develop a modified SEIR compartmental mathematical model for prediction of COVID-19 epidemic dynamics considering different intervention scenarios which might give insights on the best interventions to reduce the epidemic risk.
Several authors have worked on mathematical modeling of the novel coronavirus. A mathematical model for Middle East respiratory syndrome coronavirus (MERS-CoV) transmission dynamics was used to estimate the transmission rates in two periods due to the implementation of intensive interventions [3][4][5][6] Further and related to this work, a Bats-Hosts-Reservoir-People transmission network model for simulating the potential transmission from the infection source to the human infection was developed [5]. This article, however, differs from [5] in the sense that (1) the compartmental models are different; (2) an additional compartment for the pathogens was included to allow for

Methods
In the study, a mathematical model of the spread and transmission of SARS-CoV-2 was formulated. We consider two interacting populations, the human population as hosts and the pathogens. The model subdivides the total human population size at time t denoted as N(t) into susceptible S(t), exposed E(t), asymptomatic infectious I A (t) , symptomatic infectious I S (t) and the recovered as R(t). The pathogen in the environment is denoted as P(t). Hence for the human population we have Studies have shown that the virus can be transmitted in two ways, namely: human to human and environment to human. The epidemic data indicates that both asymptomatic I A (t) and symptomatic I S (t) infected individuals spread the COVID-19 virus to susceptible persons S(t) with whom they are in close contact. In addition, when infected individuals sneeze or cough, without taking the necessary precautions, the virus spreads to the environment they are in. Since the pathogen P(t) is known to survive in the environment for some days, susceptible individuals S(t) in close contact to this environment are likely to get exposed to these pathogens, especially in the early days of the COVID-19 outbreak before hygiene protocols are implemented. In the process of disease spread, the susceptible individual first moves to the exposed population E(t) since the host has an incubation period [5]. The exposed individual moves to either asymptomatic I A (t) or symptomatic I S (t) infectious population. P(t) is the number or quantity of pathogens present during interaction of human beings at time t. The majority of infectious individuals recovers and move to the recovered human population R(t).
The compartmental model depicting the interaction between the human population, and the pathogens in the environment is shown in Fig. 1.
The parameters used in the COVID-19 transmission model are given in Table 1.
The model culminates to a six-dimensional system of ordinary differential equations as follows.
The human population is born into the susceptible population at a rate b. The terms β 1 SP and β 2 S(I A + I S ) describes the rate at which susceptible individuals S(t) gets infected by pathogens in the environment, P(t) and from infectious humans I A (t) and I S (t) respectively. Health experts and governments have been advising people, during this outbreak, to minimize contact with infectious individuals through social distancing. Therefore in our model we propose to have new infections occur in the form β 1 SP 1+α 1 P and β 2 S(I A +I S ) 1+α 2 (I A +I S ) respectively, where the interaction proportions α 1 and α 2 denotes reciprocal of the frequency with which susceptible individuals gets infected with COVID-19 from the environment and from infectious individuals, respectively.

Equilibria and basic reproduction number of the SEIR-P model
The relevant equilibrium points are obtained by solving the equations in (1) when the left hand side is equated to zero. (1)

Existence of disease-free-equilibrium point (DFE)
In this case I A = I S = P = 0 , which implies that E = 0 and R = 0 too. Hence we have: Therefore DFE is given by b µ , 0, 0, 0, 0, 0 .

The basic reproduction number
The basic reproduction number, usually denoted as R 0 defines the average number of secondary infections caused by an individual in an entirely susceptible population. This number indicates whether the infection will spread through the population or not. The next generation matrix approach is used to obtain R 0 . Let x = (E, I A , I S , P) T then the model can be written as Evaluating the derivatives of F and V at the disease-free equilibrium point, obtained above, yields FV −1 as below: The reproduction number, R 0 , is the spectral radius of the product FV −1 which is given by; Denoting the basic reproduction numbers for human as R h 0 and for pathogens as R p 0 , we make the following deductions: Therefore, Notice that the basic reproduction number R 0 consists of two parts, representing the two modes of transmission of the coronavirus.

Results
In this section, we approximate solutions to the model equations (1)  During the first 10 days, the number of susceptible humans declines rapidly as the number of exposed individuals increases rapidly due to contact with infected individuals ( I A and I S ) and also the virus in the environment (P). After the latency period, and without mitigating the epidemic, the number of infected individuals surges, surpassing the hospital bed capacity, set here as 8000. The infected individuals who exhibit mild or no symptoms I A are considered to be 30% of the total infected population. The model parameters used in this simulation study are shown in Table 1.
Since the symptomatic individuals I S are assumed to be more infectious than the asymptomatic I A , the transmission of COVID-19 through contacts in households, workplaces, schools, from foodstuffs, or during commute rises. This leads to a surge of the virus in environments such as workplace, school, foodstuffs, and public transport, see Fig. 2b, and consequently more cases of the coronavirus are confirmed, see Fig. 2a between 10 and 35 days. In this model we take the constants α 1 and α 2 to be reciprocal of the frequency with which individuals acquires the COVID-19 from the environment and from infected individuals, respectively. In Fig. 2 the model shows that when α 1 = 0.05 i.e. there is a high risk of getting infected by a contaminated environment, as compared to an infected individual, the number of exposed, asymptomatic and symptomatic individuals increases. However, when α 2 = 0.05 i.e. higher chances of getting infected by an individual, as compared to a contaminated environment. Moreover, in Fig. 2c the number of susceptible vanishes by the 23rd day for α 2 = 0.05 since many people were infected quite rapidly, see Fig. 2d-f for duration 0 − 20 days. Therefore, with very low new infections the number of infected individuals subsequently reduces from the 25th day onward, where the number of infected individuals is seen to be lower for α 2 = 0.05 , as compared to when α 1 = 0.05 and α 2 = 0.1.

Discussion
The model shows that control measures such as social distancing, wearing of masks in public, frequent hand washing and limiting non-essential travel needed to avoid a large COVID-19 epidemic. There is a growing concern that this disease could continue to ravage the human population globally since many aspects of the COVID-19 are yet to be discovered, which also poses a challenge to the long-term mathematical modeling of the disease.

Limitations
This model was designed to look at transmission dynamics so does not describe disease severity and death. Given we made assumptions of the parameters at onset of the pandemic, there is a possibility that the model may overestimate the pandemic at later period of time.