 Research note
 Open Access
 Published:
SEIR model for COVID19 dynamics incorporating the environment and social distancing
BMC Research Notes volume 13, Article number: 352 (2020)
Abstract
Objective
Coronavirus disease 2019 (COVID19) is a pandemic respiratory illness spreading from persontoperson caused by a novel coronavirus and poses a serious public health risk. The goal of this study was to apply a modified susceptibleexposedinfectiousrecovered (SEIR) compartmental mathematical model for prediction of COVID19 epidemic dynamics incorporating pathogen in the environment and interventions. The next generation matrix approach was used to determine the basic reproduction number \(R_0\). The model equations are solved numerically using fourth and fifth order Runge–Kutta methods.
Results
We found an \(R_0\) 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 (SARSCoV2) 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 (COVID19) a global pandemic [2].
The objective of this study was to develop a modified SEIR compartmental mathematical model for prediction of COVID19 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 (MERSCoV) 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 BatsHosts–ReservoirPeople 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 nonlinear interactions between humans and the environment; and (3) thorough simulation studies were performed.
Main text
Methods
In the study, a mathematical model of the spread and transmission of SARSCoV2 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 \( N(t)=S(t)+E(t)+I_A(t)+I_S(t)+R(t).\)
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 COVID19 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 COVID19 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 COVID19 transmission model are given in Table 1.
The model culminates to a sixdimensional system of ordinary differential equations as follows.
with the initial conditions: \(S(0)>0, E(0)>0, I_A>0, I_S>0, R(0)=0, P(0)>0\)
The human population is born into the susceptible population at a rate b. The terms \(\beta _1SP\) and \(\beta _2S\left( I_A+I_S\right) \) 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 \(\frac{\beta _1SP}{1+\alpha _1P}\) and \(\frac{\beta _2S\left( I_A+I_S\right) }{1+\alpha _2\left( I_A+I_S\right) }\) respectively, where the interaction proportions \(\alpha _1\) and \(\alpha _2\) denotes reciprocal of the frequency with which susceptible individuals gets infected with COVID19 from the environment and from infectious individuals, respectively.
Equilibria and basic reproduction number of the SEIRP model
The relevant equilibrium points are obtained by solving the equations in (1) when the left hand side is equated to zero.
Existence of diseasefreeequilibrium 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 \(\left( \frac{b}{\mu },0,0,0,0,0\right) \).
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=\left( E,I_{A},I_{S},P\right) ^{T}\) then the model can be written as \(\frac{dx}{dt}=F\left( x\right) V\left( x\right) \), where
Evaluating the derivatives of F and V at the diseasefree equilibrium point, obtained above, yields \(\mathbf{FV} ^{1}\) as below:
where \(C_{1}=\psi +\mu +\omega \), \(C_{2}=\mu +\sigma +\gamma _{S}\) and \(C_{3}=\mu +\sigma +\gamma _{A}\). The reproduction number, \(R_0\), is the spectral radius of the product \(\mathbf{FV} ^{1}\) which is given by;
Denoting the basic reproduction numbers for human as \(R_{0}^{h}\) and for pathogens as \(R_{0}^{p}\), 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) using fourth and fifth order Runge–Kutta methods which are implemented via the ode45 function in MATLAB. The initial values used are \(S(0)=93000,E(0)=1000,I_A(0)=50,I_A(0)=50,I_S(0)=50,R(0)=0,P(0)=500\). Figure 2a depicts the change in the populations as time increases from 0 to 90 days. 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 COVID19 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 \(\alpha _1\) and \(\alpha _2\) to be reciprocal of the frequency with which individuals acquires the COVID19 from the environment and from infected individuals, respectively. In Fig. 2 the model shows that when \(\alpha _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 \(\alpha _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 \(\alpha _2=0.05\) since many people were infected quite rapidly, see Fig. 2d–f for duration \(020\) 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 \(\alpha _2=0.05\), as compared to when \(\alpha _1=0.05\) and \(\alpha _2=0.1\).
Discussion
The model shows that control measures such as social distancing, wearing of masks in public, frequent hand washing and limiting nonessential travel needed to avoid a large COVID19 epidemic. There is a growing concern that this disease could continue to ravage the human population globally since many aspects of the COVID19 are yet to be discovered, which also poses a challenge to the longterm 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.
Availability of data and materials
The data used to support findings of this study is available from the corresponding author upon request.
Abbreviations
 COVID19:

Coronavirus disease 2019
 MERSCoV:

Middle East respiratory syndrome coronavirus
 SARS:

Severe acute respiratory syndrome
 SARSCoV2:

Severe acute respiratory syndrome coronavirus 2
 SEIR:

Susceptible exposed infectious recovered
 WHO:

World Health Organization
References
 1.
Chen N, Zhou M, Dong X, Qu J, Gong F, Han Y, Qiu Y, Wang J, Liu Y, Wei Y, et al. Epidemiological and clinical characteristics of 99 cases of 2019 novel coronavirus pneumonia in Wuhan, China: a descriptive study. Lancet. 2020;395(10223):507–13.
 2.
WHO. Coronavirus disease 2019 (COVID19): situation report, 51. 2020.
 3.
Hsieh YH. 2015 Middle East respiratory syndrome coronavirus (MERSCoV) nosocomial outbreak in south korea: insights from modeling. PeerJ. 2015;3:1505.
 4.
Kim Y, Lee S, Chu C, Choe S, Hong S, Shin Y. The characteristics of middle eastern respiratory syndrome coronavirus transmission dynamics in South Korea. Osong Public Health Res Perspect. 2016;7(1):49–55.
 5.
Chen TM, Rui J, Wang QP, Zhao ZY, Cui JA, Yin L. A mathematical model for simulating the phasebased transmissibility of a novel coronavirus. Infect Dis Poverty. 2020;9(1):1–8.
 6.
Sookaromdee P, Wiwanitkit V, et al. Imported cases of 2019novel coronavirus (2019ncov) infections in Thailand: mathematical modelling of the outbreak. Asian Pac J Trop Med. 2020;13(3):139–40.
Acknowledgements
The authors appreciate ample time given by their respective universities towards this manuscript.
Funding
This study has not received any funding.
Author information
Affiliations
Contributions
SM was involved in writing, reviewing and approving the final manuscript. MM and VO were involved in result production and analyzing the results. DG and RM were involved in drafting the manuscript. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publish
Not applicable.
Competing interests
The authors declare that they have no competing interests.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Mwalili, S., Kimathi, M., Ojiambo, V. et al. SEIR model for COVID19 dynamics incorporating the environment and social distancing. BMC Res Notes 13, 352 (2020). https://doi.org/10.1186/s13104020051921
Received:
Accepted:
Published:
Keywords
 SEIR model
 COVID19 dynamics
 Social distancing
 Mathematical model
 Basic reproduction number
 Runge–Kutta method