key: cord-0897916-t9num9n9
authors: Bhadauria, Archana Singh; Pathak, Rachana; Chaudhary, Manisha
title: A SIQ mathematical model on COVID-19 investigating the lockdown effect
date: 2021-01-07
journal: Infect Dis Model
DOI: 10.1016/j.idm.2020.12.010
sha: 8df73e9a1f006937e2412f150dcf91b477a255b2
doc_id: 897916
cord_uid: t9num9n9

This research paper aims at studying the impact of lockdown on the dynamics of novel Corona Virus Disease (COVID-19) emerged in Wuhan city of China in December 2019. Perceiving the pandemic situation throughout the world, Government of India restricted international passenger traffic through land check post (Liang, 2020) and imposed complete lockdown in the country on 24 March 2020. To study the impact of lockdown on disease dynamics we consider a three-dimensional mathematical model using nonlinear ordinary differential equations. The proposed model has been studied using stability theory of nonlinear ordinary differential equations. Basic reproduction ratio is computed and significant parameters responsible to keep basic reproduction ratio less than one are identified. The study reveals that disease vanishes from the system only if complete lockdown is imposed otherwise disease will always persist in the population. However, disease can be kept under control by implementing contact tracing and quarantine measures as well along with lockdown if lockdown is imposed partially.

The novel Corona Virus Disease (COVID- 19) was first manifested in Wuhan, China in late 2019 that has travelled across different countries and spread at quite a fast pace by human to human infection [4, 6, 19] . Within a very short span of time the disease has spread in major part of the world leading to large number of deaths. Migration of individuals from China to other parts of world played a critical role in the spread of infectious disease COVID-19. China could control the disease by introducing lockdown in Wuhan since, it was only option available due to the absence of targeted drug and vaccine against the virus. This pandemic led the Government of India to impose complete lockdown in the country on March 24, 2020 and is still in lockdown stage as of May 17, 2020 . Literatures report that lockdown measure is significant because it has been able to reduce the transmission of COVID-19 by as much as 84 percent as in France [5] . Principle feature of lockdown is reduction of immigration of population to nearly zero. It is quite evident from the past that immigration of population plays a significant role in the spread of infectious diseases like flu. It has been reported that immigration has not only played roles in introducing infections into native populations but also in changing the incidence of infections, and in increasing the prospects of local transmission [27] . Immigration of infective population from COVID-19 region has posed risk to the local population and increased the geographic range of the disease at such a fast pace that it spread in 28 countries by 21 February 2020 [28] . This article focuses on the impact of lockdown as preventive measure through mathematical modeling approach. Modeling the effect of lockdown to prevent the spread of disease is a complex task however we have modelled it by introducing the immigration term in the model and comparing it with the model without immigration term. Contact tracing, immigration of individual has been taken as one of the most important features in studying the disease dynamics by some of the mathematical modellers given in [12, 14, 18, 30, 31, 32] . We have introduced immigration of susceptible individuals along with a fraction of infective individuals directly in the infective compartment and compared it with the model without immigration. In addition, we have also considered contact tracing and quarantine as other preventive measures in our mathematical model. Quarantine is one of the live options during disease outbreaks. Fraser et al in [15] , shows that the quarantine is an effective measure if asymptomatic infection are responsible for a large number of infections as seen in COVID-19 cases. The author further concluded in his study on SARS that its outbreak could be controlled effectively by sufficient isolation of symptomatic individuals. On the other hand, Chen et al. [10] reported that even with less efficient contact tracing an outbreak of SARS with effective isolation of symptomatic patients is sufficient to control SARS outbreak. Thus, objective of this paper is to study the combined effect of lockdown, contact tracing and quarantine on the COVID-19 dynamics. These measures are important since there is no specific antiviral treatment up to date hence prevention is the only way to keep the disease under control. Mathematical modelling can help J o u r n a l P r e -p r o o f by evaluating the effectiveness of these preventive measures in the early phase of epidemics of new emerging outbreak in the future [16] .

And to the best of our knowledge combined effect of all these measures, lockdown, contact tracing and quarantine has not been considered in any of the previous mathematical model on COVID-19. Our research paper is organised as follows: In section 2 mathematical model is discussed, section 3 represents boundedness of the system. In section 4 model dynamics without lockdown effect is studied followed by section 5 manifesting the lockdown effect in the system. In section 6, numerical analysis is performed and paper is concluded in result and discussion section 7.

We consider a three-dimensional model with susceptible population S(t), infective population I(t) and quarantined population Q(t). It is assumed that the population is homogeneously mixed and disease spreads via the direct contact between susceptible and infective individuals as well by immigration of individuals in susceptible class and in infective class. Thus, we have an SIQ model with immigration. Let N (t) be the total population at any time t in the region under consideration, which is the sum of three subpopulation. Susceptible population is entering in the system at the constant rate A and by immigration (1−θ)m. Whereas fraction θm is entering directly to the infective class with the assumption that a fraction of immigrants are infective are carrier of virus. We have studied the impact of lockdown in the model by taking immigration term as zero because immigration is restricted to a large extent during lockdown period. Moreover, by contact tracing, some of the infective individuals are quarantined and only fraction of the infective's are entering into the class of infective and spreading infection among the susceptible. We model here the average effect of contact-tracing based strategy [2] . Let k be the rate of contact tracing of individuals so that a fraction (1 − k) are responsible for spreading infection and rest are quarantined. We have assumed Holling type-II functional response to model the interaction between susceptible and infective populations, where α and η are the infection catching and handling time of the virus from the infective individual to the susceptible individuals. µ is the natural mortality rate of population in each compartment, δ 1 and δ 2 are the recovery rate of the infective and quarantined population whereas α 1 and α 2 are the disease related death rates of the infective and quarantined compartments respectively. Infected individuals are assumed to be quarantined at the rate σ. Thus, with these assumptions we have formulated a three dimensional epidemiological model taking susceptible S(t), infective I(t) and quarantined Q(t) population as system variables. A schematic flow figure 1 is represented to understand the dynamics of COVID-19.

The model equations with immigration are given by

J o u r n a l P r e -p r o o f 

with initial conditions S(0) > 0,

Description of parameters is defined in table 1 rate of recovery of infective population δ 2 rate of recovery of quarantined population σ transition rate from infective to quarantined population α 1 disease related death rate of infective population α 2 disease related death rate of quarantined population

In the following lemma, we state the bounds of system variables which is further used to prove the analytic results. Boundedness of the system may be interpreted as a natural restriction to indefinite growth of infective population because of various constraints either due to natural conditions or due to preventive habits acquired by the population to keep themselves safe from the disease. We will show that the solutions of system are bounded to ensure the biological validity of the model. 

According to comparison principle, it follows that N m = A+m µ . This completes the proof.

When m = 0 then the model (1)-(3) has no disease-free equilibrium point, only pandemic equilibrium exists. Thus, when immigration of population is allowed in the system and a fraction of them are infective; disease cannot be vanished to zero, although attempts may be made to control the disease.

In a dynamical system an equilibrium point is a state of the system that does not change with time. Thus, if the system starts at an equilibrium point, the state will remain at the equilibrium forever. We are interested in finding the equilibrium state of system when it is free from disease that is, disease free equilibrium point and the equilibrium state of persistence of disease, the endemic/pandemic equilibrium point. 

and I m is the root of the quadratic equation (4) 

which guarantees that unique positive equilibrium point E m exist iff R 0 > 1.

J o u r n a l P r e -p r o o f

Local stability analysis of equilibrium points gives an insight into the impact of small perturbations in the equilibrium state of the system. To study the stability of pandemic equilibrium point, we linearize the system about the arbitrary equilibrium points E m (S m , I m , Q m ) and obtain the corresponding Jacobian matrix as given by

Proof: From the Jacobean matrix J E m we observe that the eigenvalues of E m (S m , I m , Q m ) are given by λ 1 = −(µ + δ 2 + α 2 ) < 0 and other two eigenvalues λ 2 and λ 2 are the roots of the quadratic equation

Thus, the equilibrium E m will be locally asymptotically stable if (b 11 +b 22 ) < 0, which is obvious from the expression of b 11 and b 22 . Since both of them are negative if and only if pandemic equilibrium state of infective population is positive and equation (4) states that infective population is positive if and only if R 0 > 1. Thus, we may conclude that the pandemic equilibrium point is locally asymptotically stable if R 0 > 1. Biologically it may be interpreted by the statement that the system will return back to the pandemic equilibrium point and stays close to it when number of secondary infections produced by a single infective during his complete infectious period is more than one.

Here in this section we found persistence of the disease. The word 'Persistence' stands for the survival of each population in future time. We study persistence of the system to determine the conditions under which a disease persists in the system. Persistence of the model (1)-(3) is proved using the following Lemma 

S min is always positive. From equation (2), we have

Again, using system boundedness and comparison principle

Again using the comparison principle,

Q min remains always positive. This completes the proof of the Lemma.

To study the impact of lockdown effect, we study the model (1)-(3) for m = 0, that is when immigration of the population in the system is strictly prohibited. In this case, the model reduces to the following system of equations (5)-(7)

with initial conditions S(0) > 0,

J o u r n a l P r e -p r o o f

Basic Reproduction Number [7, 9, 13, 25] is a threshold number that determines number of secondary infections produced by an infective individual during his or her whole infectious period in a population in which every individual is susceptible. We determine basic reproduction number using next generation matrix approach [13] . The decomposition of the model into components R 1 and R 2 leads to the system (5)-(7)of the form

Since, the infected compartments are I and Q, at the disease free equilibrium point, we defineR

Note that,R 1 is non negative andR 2 is a non singular matrix withR 2 −1 as non-negative andR 1R2 −1 as non negative matrix;R 1R2 −1 is the next generation matrix. In this case if . Hence, Basic reproduction number is given by

.

From the mathematical expression of basic reproduction ratio we observe that it is directly proportional to the susceptible population at disease free state A µ , infectious period 1 (µ+δ1+σ+α1) and transmission rate of infection β. Along with these terms it is inversely proportional to rate of contact tracing k and half saturation constant α. Thus, spread of disease can be controlled by keeping a check on these parameter values.

Based on each of the parameters of R 0 , sensitivity analysis is performed to check the sensitivity of the basic reproduction number [3, 17] . In order to reduce the effect of corona virus, it is necessary to control the parameter values to make R 0 < 1. We are therefore interested in finding the rate of change of R 0 as the parameter values are changed. The rate of change of R 0 for a change in value of parameter h can be estimated from SI[h] = h R0 × ∂R0 ∂h . The normalized sensitivity indices of the reproduction number with respect to parameter k, α, µ, δ 1 , σ, α 1 are given by

from the above calculations it is obvious that R 0 is mainly sensitive to alters in β. The value of R 0 will be enhanced if we raise the value of β and rest of indices are negative i.e., the value of R 0 decreases as k, α, δ 1 , σ, α 1 and µ increase. Thus, basic reproduction number increases proportionally with the increase in rate of transmission of infection and decreases with the rate of contact J o u r n a l P r e -p r o o f tracing, half saturation constant, recovery rate of infective population, rate of transition from infective to quarantined population, disease related death rate of infective population and natural death rate of each population.

System (5)-(7) has two equilibrium points, one is the Disease free equilibrium point and another is pandemic equilibrium point.

Theorem: The SIQ model (5)-(7) has a unique Disease free equilibrium E 0 A µ , 0, 0 for all parameter values. As well as the system under consideration has also a unique pandemic equilibriumÊ(Ŝ,Î,Q) under the condition R 0 > 1. Proof: Pandemic equilibrium point is given byÊ(Ŝ,Î,Q), using equation (5) and (7) we haveŜ = A(α+ηÎ) µα+(β+µη)Î , andQ = βAÎ (µ+δ2+α2)[µα+(β+µη)Î] + σÎ µ+δ2+α2 respectively. Next putting the value ofŜ in equation (6), we getÎ 2 − (R0−1)µαÎ β+µη = 0. From above equation, we getÎ = 0 orÎ = (R0−1)µα β+µη . The equilibriumÊ exists if R 0 > 1.

To study the stability of the different equilibrium point the corresponding Jacobian matrix is given as: 

Proposition: For the system of equation (5)-(7), the disease free steady state E 0 is locally asymptotically stable if R 0 < 1. Proof: Eigenvalues corresponding to E 0 are given by

and λ 3 = (R 0 − 1)(µ + δ 1 + σ + α 1 ).

Thus the equilibrium E 0 will be locally asymptotically stable if R 0 < 1. Thus, we infer that if basic reproduction number is less than one, disease free equilibrium point is locally asymptotically stable. It implies that when number of secondary infections are less than one, a small perturbation in the disease free equilibrium state of the system always tends to move towards the disease free state and hence disease vanishes from the system.

Proposition: For the system of equation (5)-(7), the disease free steady state E 0 is globally asymptotically stable in the region R 2 + of I -Q plane if R 0 < 1. Proof: For the system of equation (5)-(7), the disease free steady state E 0 is globally asymptotically stable as

Let us consider g(I, Q) = 1 IQ and L(I, Q) = ∂ ∂I (gf 1 )+ ∂ ∂Q (gf 2 ). Since g(I, Q) > 0, for all I, Q > 0. So we have

Clearly, L(I,Q) does not change it's sign and is not zero identically in positive Quadrant of I − Q plane. Thus, by Bendixon-Dulac criteria, the system has no limit cycle in the positive quadrant of I − Q plane. Hence, a disease free equilibrium is always locally asymptotically stable when R 0 < 1, and is globally stable in (I − Q) plane of the region R 2 + if R 0 < 1. This theorem can be biologically interpreted as, when basic reproduction number is less than one system always tends towards the disease free equilibrium state no matter how much large the perturbation is or from whatever initial state we start, system will approach disease free equilibrium state and disease vanishes from the system.

Here, in this section we found the local stability of the pandemic equilibrium pointÊ(Ŝ,Î,Q), where

We assume β = β * , as the bifurcation parameter, predominantly as it has been explained in [11] that when R 0 = 1, we have

Now, the Jacobian of the system at disease free equilibrium E 0 when β = β * is as follows:

the eigenvalues of the matrix are given by (0, −µ, −α 2 − δ 2 − µ) T .

We observe that the matrix has simple zero eigenvalues, and the other eigenvalues are negative. We are now at the stage to apply center manifold theory [26] to analyze the dynamics of system (5)- (7) . Corresponding to zero eigenvalue the right eigenvector w = (w 1 , w 2 , w 3 ) T of the matrix is given by

Substitute the values of all the second-order derivatives calculated at disease free equilibrium as well as β = β * , we obtaina = 2(1−k)β α w 1 w 2 v 2 < 0 using the value of w 1 , w 2 and v 2 = 1. b = v2w2(1−k)A µα > 0. As a < 0 and b > 0 at β = β * , hence from [8, 9] a transcritical bifurcation takes place at R 0 = 1 as well as the inimitable pandemic equilibrium is locally asymptotically stable due to R 0 > 1. Hence, the proof. Thus, we may conclude that for a particular J o u r n a l P r e -p r o o f value of rate of transmission of infection, basic reproduction number becomes one and at this point equilibrium state changes its nature, when basic reproduction number is less than one, pandemic equilibrium is unstable and when it is greater than one pandemic equilibrium is stable subject to a small perturbation from the pandemic equilibrium state. It implies that disease spreads among the population when basic reproduction number is more than one.

In this section we have discussed the quantitative behavior of the model during lockdown in India. With and without immigration of individuals in the country we justify the analytical findings and perform numerical simulation using MATLAB and Mathematica softwares. Table 2 , describes the parameter used to perform numerical simulation. Most of the parameters have been taken from the previous literatures and some of them are assumed. 1.78 × 10 −5 per day [29] We have calculated the equilibrium points of the system without immigration when complete lockdown was imposed by the government. In this case system possess two equilibrium points, one is disease free equilibrium point (1.282750 × 10 9 , 0, 0) and the other is the pandemic equilibrium point (1.28237 × 10 9 , 14.03000, 35.59000). Basic reproduction ratio R 0 is computed and is found to be 3.80. All the conditions of local and global asymptotic stability are verified by the above set of parameters. However, when lockdown is not imposed and immigration of individuals is allowed in the system then our model possess only one equilibrium point given by (1.32976 × 10 9 , 22745.39525, 33564.48733). Above set of parameters satisfy the conditions of local stability of the equilibrium point and persistence. In addition, to study the impact of various parameters on the infective population we plotted time series graphs of infective population. Figure 2 shows the variation of infective population with time for different number of immigrants. It is observed that infective population rise significantly with the increase in number of immigrants. Hence it is imperative to impose complete/partial lockdown in the country to decrease the number of infective population. In the absence of lockdown, the spread of COVID-19 cases will shoot up. From figure 3 , we determine the impact of transition rate of infective individuals to quarantine compartment. It is observed that as transition rate increases number of infective individuals decrease. It implies that the rate σ at which infective individuals are quarantined plays a significant role in controlling the number of infective population and hence to control the disease. To determine the impact of various parameters on the threshold parameter Basic Reproduction number, we plot figures depicting variation of R 0 with different parameters. In figures 6a, 6b, 6c and 6d we plot the variation of R 0 with rate of contact tracing k, transmission rate of infection β, recruitment rate A, half saturation constant α respectively.

We observe that with the increase in rate of contact tracing k to 0.9, basic reproduction number reduces below one. Hence by the practice of contact tracing disease can be controlled. Further with the increase in β, basic reproduction number continues to rise linearly. This implies that if transmission rate of infection will rise continuously, we will not be able to control the disease and it will spread in the population. Further R 0 rises with the increase in recruitment rate A. Reproduction ratio R 0 decreases with the increase in half saturation constant, α. In figure 6e, 6f, 6g and 6h display the variation of reproduction number with natural mortality rate µ, recovery rate of infective's δ 1 , transition rate of infective's to quarantined class σ and disease related death rate of infective individuals α 1 . It is observed that R 0 decreases with the increase in these parameter values. Thus, from this study we confer that basic reproduction ratio can be reduced and maintained below one by controlling the parameters k, α, µ, β, σ, δ 1 and α 1 .

This paper represents a nonlinear mathematical model to study the impact of lockdown on the spread of COVID-19 with significant preventive measures adopted by the Indian government to control the spread of disease. Study of the proposed model is divided into two phases; complete lockdown state and partial lockdown or unlocked state. Qualitative analysis of the model in both the cases is performed and basic reproduction ratio is determined through next generation matrix approach. Nature of equilibrium points is found to be tion number is less than one and pandemic equilibrium exists only if basic reproduction number is greater than one. However, in the case of unlocked or partial lockdown state, diseases free equilibrium does not exist and disease will persist in the system. Thus, lockdown is most important measure to be imposed strictly until infection dies out from the country to reduce the infection to zero. Local stability of the pandemic equilibrium and conditions of persistence are determined in case of unlock or partial lockdown state. Further, to justify the analytical findings, numerical simulation is performed and graphs are plotted. It is observed that infective population rise significantly with the increase in number of immigrants in case of partial lockdown state. However, this situation can be kept under control by contact tracing and quarantine, since numerical simulation confirms the decrease in infective population with the increase in their rate. Hence, contact tracing and quarantine play an important role is keeping the disease under control in case of unlock or partial lockdown state. But in order to completely eradicate the disease, lockdown is most important measure to be followed by the government. In addition, we have identified the significant parameters by which basic reproduction number can be reduced and maintained below one. Contact tracing and half saturation constant play an important role in keeping the basic reproduction number below one along with the transmission rate of infection, susceptible population in disease free state and infectious period of an infective individual. Sensitivity analysis of the basic reproduction number is carried out and sensitivity index of each parameter is determined. It is observed from the analysis that transmission rate of infection is most sensitive factor in basic reproduction number and acts as transcritical bifurcation parameter. Our study reveals that complete lockdown is mandatory for reducing the infection to zero and control in immigration of population in the system is an outstanding step to control the spread of the disease. Still, the complete lockdown is not the ultimate solution since it may cause financial crisis in the near future and hence it is not possible to implement complete lockdown in highly dense country like India for a long tenure. Therefore, keeping an eye on the immigration of infective population in the country, government must emphasize on other preventive measures too like increase in rate of contact tracing and quarantining the confirmed cases of COVID-19 on a regular basis until infection dies out.

Population biology of infectious diseases, Part I

Preventive treatment of tuberculosis through contact tracing

Forward and adjoint sensitivity analysis: with applications in Dynamical Systems

Predictive Mathematical Modelling, medrxiv

Role of immigrants and migrants in emerging infectious diseases

Coro-naTracker: world-wide COVID-19 outbreak data analysis and prediction

Application of center manifold Theory

On the computation of ro and its role on, Mathematical approaches for emerging and reemerging infectious diseases: an introduction

Predictive models of control strategies involved in containing indoor airborne infections

Preliminary epidemiologic assessment of mers-cov outbreak in south korea

When is quarantine a useful control strategy for emerging infectious diseases

Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission

Mathematical modeling of COVID-19 transmission dynamics with a case study of Wuhan

Factors that make an infectious disease outbreak controllable

Prediction and analysis of coronavirus disease

Global stability and sensitivity analysis of SEIQR worm virus propagation model with quarantined state in mobile internet

Epidemic Models of Contact Tracing: Systematic Review of Transmission Studies of Severe Acute Respiratory Syndrome and Middle East Respiratory Syndrome, Computational and structural biotechnology journal

Early Transmission Dynamics in Wuhan, China, of Novel Coronavirus-Infected Pneumonia

Mathematical model of infection kinetics and its analysis for COVID-19

A model based study on the dynamics of COVID-19: Prediction and control

COVID-19 and tuberculosis: A mathematical model based forecasting in Delhi

Time series modelling to forecast the confirmed and recovered cases of COVID-19

From heroin epidemics to methamphetamine epiepidemics: Modelling substance abuse in a South African province

Mathematical Analysis of a COVID-19 Epidemic Model by using Data Driven Epidemiological Parameters of Dseases Spread in India, medRxiv

First cases of coronavirus disease 2019 (COVID-19) in the WHO European Region

Epidemic Models with Infective Immigrants and Vaccination

Modeling containing covid-19 infection: A conceptual model

Population dispersal and disease spread

An epidemic model with population dispersal and infection period

SIV epidemic models with age of infection

J o u r n a l P r e -p r o o f

☒ 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.☐The authors declare the following financial interests/personal relationships which may be considered as potential competing interests: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.