|Year : 2021 | Volume
| Issue : 1 | Page : 175-181
Modelling the spread of SARS-CoV-2 pandemic - Impact of lockdowns & interventions
Manindra Agrawal1, Madhuri Kanitkar2, M Vidyasagar3
1 Department of Computer Science & Engineering, Indian Institute of Technology Kanpur, Kanpur, Uttar Pradesh, India
2 Deputy Chief Integrated Defence Staff (Medical), HQ Integrated Defense Staff, Ministry of Defence, Government of India, New Delhi, India
3 Department of Artificial Intelligence, Indian Institute of Technology Hyderabad, Hyderabad, Telangana, India
|Date of Submission||28-Sep-2020|
|Date of Web Publication||26-Mar-2021|
Dr. Madhuri Kanitkar
Deputy Chief Integrated Defence Staff (Medical), HQ Integrated Defense Staff, Ministry of Defence, Government of India, New Delhi 110 011
Source of Support: None, Conflict of Interest: None
| Abstract|| |
Background & objectives: To handle the current COVID-19 pandemic in India, multiple strategies have been applied and implemented to slow down the virus transmission. These included clinical management of active cases, rapid development of treatment strategies, vaccines computational modelling and statistical tools to name a few. This article presents a mathematical model for a time series prediction and analyzes the impact of the lockdown.
Methods: Several existing mathematical models were not able to account for asymptomatic patients, with limited testing capability at onset and no data on serosurveillance. In this study, a new model was used which was developed on lines of susceptible-asymptomatic-infected-recovered (SAIR) to assess the impact of the lockdown and make predictions on its future course. Four parameters were used, namely β, γ, η and ε. β measures the likelihood of the susceptible person getting infected, and γ denotes recovery rate of patients. The ratio β/γ is denoted by R0 (basic reproduction number).
Results: The disease spread was reduced due to initial lockdown. An increase in γ reflects healthcare and hospital services, medications and protocols put in place. In Delhi, the predictions from the model were corroborated with July and September serosurveys, which showed antibodies in 23.5 and 33 per cent population, respectively.
Interpretation & conclusions: The SAIR model has helped understand the disease better. If the model is correct, we may have reached herd immunity with about 380 million people already infected. However, personal protective measures remain crucial. If there was no lockdown, the number of active infections would have peaked at close to 14.7 million, resulted in more than 2.6 million deaths, and the peak would have arrived by June 2020. The number of deaths with the current trends may be less than 0.2 million.
Keywords: Asymptomatic cases - India - lockdown - modelling - resources
|How to cite this article:|
Agrawal M, Kanitkar M, Vidyasagar M. Modelling the spread of SARS-CoV-2 pandemic - Impact of lockdowns & interventions. Indian J Med Res 2021;153:175-81
|How to cite this URL:|
Agrawal M, Kanitkar M, Vidyasagar M. Modelling the spread of SARS-CoV-2 pandemic - Impact of lockdowns & interventions. Indian J Med Res [serial online] 2021 [cited 2021 Apr 18];153:175-81. Available from: https://www.ijmr.org.in/text.asp?2021/153/1/175/298895
COVID-19 pandemic has affected India as much as it has the rest of the world and the updated data for India are available in the public domain,,. India, in spite of being the second most populated country in the world with less than ideal healthcare infrastructure, has been fortunate to have had a fairly low case-fatality rate until now,. The role of various measures to reduce transmission is now better defined,,,,. The exact reasons for the outcome in India are possibly multidimensional, but the various steps implemented by the Policy makers to slow down transmission of the virus have also been important,,. Besides aiding in flattening the curve, these have given the scientific community an opportunity to better understand the vector host environment dynamics,. These actions have also provided time for the existing healthcare organizations to plan, activate and reallocate resources to help prepare for the surge. There was an urgent need to be ready to cater to a large number of patients over a short period of time. The time gained through various government policies has also helped nurture individual innovators and industry partnerships to achieve self-sufficiency in resources, such as the personal protective equipment, sanitizers and hospital beds, at a time when the world grappled with the crises,. We have presented a mathematical model for a time series prediction and to analyze the impact of a lockdown in the country.
| Material & Methods|| |
A consultative committee was constituted by the Department of Science and Technology under the Ministry of Science and Technology, Government of India, to develop a supermodel consisting of mathematical predictions as related to the COVID-19 pandemic in India. During the deliberations of the committee, there were extensive discussions and literature review of the evolving pandemic and experience from other countries,,,,. Several mathematical models submitted to the committee, for the spread of a pandemic were analyzed in detail and the gaps identified, namely the inability of the standard susceptible-infected-recovered (SIR)/susceptible-exposed-infected-recovered (SEIR) models to account for asymptomatic patients, which were a novel feature of COVID-19. There is a paucity of data on the behaviour of the virus among Indian population,,. Limited testing capability in India at the time of onset of the pandemic,, non-availability of standardized tests for serosurveillance and non-availability of data on asymptomatic cases were other limitations. In the current COVID-19 pandemic, a large fraction of population showed little or no symptoms. The standard mathematical models could not explain the role of asymptomatic cases in unfolding of the pandemics.
A new susceptible-asymptomatic-infected-recovered (SAIR) model: Based on an in-depth analysis, a mathematical model was proposed, which built on an earlier model,. Before this, another mathematical model provided estimates of the inventories of vital equipment required to cope with the pandemic. This model developed on lines of susceptible-asymptomatic-infected-recovered (SAIR) was used to assess the magnitude of the impact of the lockdown on the disease burden. In addition, predictions on the future course of the pandemic were also made.
The standard model for pandemic dynamics called SIR classifies the population in three categories: susceptible (S), infected (I) and removed (R). A subcategory of removed is deceased (D). To differentiate between asymptomatic and symptomatic patients, a new category of asymptomatic (A) was introduced with the population in category I of SIR model divided into
A (asymptomatic patients) and I (symptomatic patients). Further, it has been assumed that whether a person is categorized in A or I upon infection solely depends on the physiology of the person (immunity level, genetic disposition, comorbidities, etc.). This splits category S into two: SA (asymptomatic upon getting infected) and SI (infection and manifesting disease upon getting infected).
Parameter estimation was based on four parameters in this model, namely β, γ, η and ε. In brief, β measures the likelihood that contacts between a susceptible person on the one hand, and an asymptomatic or an infected person on the other hand, leads to the susceptible person getting infected. The parameter γ corresponds to the rate at which asymptomatic or infected patients recover. The parameter η measures the rate at which infected patients die (the presumption is that all asymptomatic patients recover). Finally, the parameter ε measures the ratio between infected symptomatic and asymptomatic patient populations.
| Results|| |
The ratio β/γ is known as the basic reproduction ratio and is denoted by R0. If R0 <1, then the pandemic does not take off because it means that asymptomatic/infected patients recover faster than they can infect the susceptible population. However, if R0 >1, then the number of infections increases at first before decreasing to a steady state value. Therefore, to make predictions using the model, it is essential to form reliable estimates of these parameters based on the available data.
For data, the model used three daily time series reported: active infections, cumulative recovered and cumulative deaths. We took these three as the daily values of I, RI and D, respectively. This may not be exactly true since many asymptomatic patients are also reported and some deaths are not reported. Four parameters were used in this model, namely β, γ, η, ε and these four parameters change with time. The parameter β is a function of interactions between people and so, lockdowns directly impact it by reducing it. The parameters γ denotes recovery rate and η the deaths which depend on the medical care available. Improvements in medical care increases γ and reduces η. Parameter ε is the ratio of infected with disease and asymptomatic cases.
In principle, the parameter ε should not change as the ratio of infected and asymptomatic remains constant in the model. However, its measurement is not exact for two reasons. First, initially, the infection is limited to a small section of population in the region under study, and it takes time for it to spread to the entire region. Since we do not know the exact size of initial infected population, we get a significant underestimate of ε initially as we use total population of the region to compute it. The measured value slowly increases to eventually stabilize once the entire region is covered - this makes ε also dependent on lockdown. Second, measured values of active infections and cumulative removed (recovered plus deceased) also include parts of A and RA, respectively. Hence, actual values of I and RI are smaller than calculated by parameter estimation. This causes the value of ε to be overestimated. Adjustments can be made to get a better estimate of ε by estimating how the infection is spreading in the region and recording the fraction of asymptomatic cases in reported positive ones.
The details of these equations and calculations as developed by the authors and part of the mandate for the Committee are as follows:
Susceptible (S) is the uninfected fraction of population. It is divided in two: SH and SL fractions that transition to I and A, respectively, on infection. Asymptomatic (A) is the fraction of population who is infected but goes undetected. Infected (I) is the fraction of population who is infected and detected. Removed (R) is the fraction of population who had infection but no longer infected (either through recovery or by demise). It is also divided in two: I and A transition to RI and RA, respectively. The denominator dt is over a continuous timeframe.
Pandemic progression in India: The estimated parameters for the five phases are depicted in [Table 1]. The mathematical model and predictions are further depicted in [Figure 1], which depicts that the actual data suggest a peak on September 17, 2020. The model overestimates the actual growth by around 1.5 per cent and the peak arrives four days later. Cumulative deaths are predicted to be around 0.2 million. The β and γ values for the last phase imply R0 value around 1.39, and the herd immunity for this value of R0 at around 28 per cent of population in A, I and R categories. This translates to around 3.9 million population infected or with antibodies at the peak. On September 17, the numbers in I and RI categories were around 5.2 million. Using the 1/ε value (=67) of phase 6, the model predicts total population with infection or antibodies to be around 3.5 million.
|Figure 1: Progression of the pandemic in India. Source: Ref. 1 for actual data.|
Click here to view
Analysis of disease progression in Delhi: The estimated parameters for the phases 1-6 for Delhi are depicted in [Table 2]. For Delhi, first phase starts from April 1, as data were available in the open domain only from March and a window of 14 days is required for calculations by the model. The Delhi plot of projections are of special interest since the data from two serosurveys can be used to validate the calculations made by this model. The mathematical model and predictions are further depicted in [Figure 2], which depicts that the projected trajectory follows the actual line closely, except that the second wave starts about a week early, first peak is about 20 per cent higher, and second is about 10 per cent higher. Cumulative deaths are to be around 6500. July serosurvey showed antibodies in 23.5 per cent of population, and RI at the start of survey was around 55000. The model predicts R (=1/ɛ×RI) to be around 5.4 million (1/ɛ=104 for phase 3 assuming the population of Delhi to be around 20 million. This is close to 27.5 per cent of population which is 4 per cent higher than actual data. September serosurvey showed antibodies in 33 per cent of population and RI at the start of the survey was around 0.16 million. The model predicts R (=1/ɛ×RI) to be around 6.5 million (1/ɛ= 41 for phase 5). This amounts to 33.7 per cent which is 0.76 per cent higher than actual data. As evident by the calculations in this model, there was a good association between the predicted numbers and the actual data obtained in the serosurvey.
|Figure 2: Progression of the infection in Delhi. Source: Ref. 1 for actual data.|
Click here to view
The effect of lockdown and the spread in active infections and hypothetical effects of timing of lockdown are depicted in [Figure 3]. It revealed that if there was no lockdown, the number of active infections would have peaked at 14+ million and the peak would have arrived by mid-May. This would have resulted in overwhelming our hospitals. The impact of delay in lockdown is starkly reflected in projected mortality and the same is depicted in [Figure 4]. If there was no lockdown, it would have resulted in more than two million deaths. The number of deaths with current trends is projected to be less than 0.2 million.
| Discussion|| |
The presence of asymptomatic patients, who can infect the susceptible population when they come in contact with, changes the dynamics of disease. Among the total population who test positive, only the symptomatic and sick patients (which is a small fraction), require major medical intervention. Thus, the estimates for medical inventory, including oxygen, ventilators and intensive care unit beds, needs to be planned based on projections of this number and not the total infected persons. At the same time, asymptomatic patients need to be advised isolation so that they do not spread the infection. In the present model, time period of the disease progression in the country was divided in six phases. These were obtained by identifying the dates on which the value of ε increased significantly: phase 1: period up to April 5; phase 2: period from April 5 to 30; phase 3: period from May 1 to June 15; phase 4: period from June 16 to July 15; phase 5: period from July 16 to August 15; and phase 6: period from August 15 to date.
These phases correspond to different levels of lockdown introduced in the country with some delay. For each of these phases, parameters were estimated for India as a whole and specifically for the city of Delhi (due to lack of early data, for Delhi, the phases began from April 1 and there were slight changes in phase boundaries).
The progression of the pandemic in India revealed that there was a significant reduction in β due to initial lockdown. It was more or less stabilized despite relaxations in lockdown. This possibly reflects the protective measures being ensured by the susceptible population all over the country such as use of mask and social distancing. An increase in γ was observed following initial lockdown. This reflects healthcare and hospital services, medications and protocols put in place. η remained stationary in the first three phases and then started reducing. This possibly reflects improved medical interventions and increasing number of asymptomatic patients being diagnosed and counted towards total cases.
The progression of the disease in Delhi revealed that the value of 1/ɛ dropped and more so from phase 1 to 3. This suggested that the virus did not initially spread over the entire population uniformly. The value of β increased significantly from phase 2 to 3 and again from phase 4 to 5. This reflected relaxation of lockdown leading to increased human interactions and increasing cases. γ increased significantly from phase 2 to 3 reflecting the major improvements in medical care in June. Decreasing η reflected improvements in medical care provided, thus reducing mortality.
The pandemic progression actual data showed a peak on September 17. Our model predicted a slightly higher peak (by ~1.5%) four days later. Cumulative deaths are predicted to be slightly less than 0.2 million. The β and γ values for phase 6 imply R0 value around 1.39 and the herd immunity for this value of R0 is at around 28 per cent of population in A, I and R categories. This translates to around 3.85 million population with infection or with antibodies at peak. On September 17, the numbers in I and RI categories were around 5.2 million. Using the 1/ε value (=67) of phase 6, the model predicts total population with either infection or antibodies to be around 350 million. If the model is correct and there is no significant increase in R0 value, we perhaps have reached herd immunity. However, to ensure that R0 stays small, it is important to continue the use of personal protective measures such as mask, hand hygiene and physical distancing.
The mathematical model used in this study was developed to analyze the spread of the disease under three different hypothetical scenarios, namely no lockdown, lockdown starting April 1, 2020 and lockdown starting May 1, 2020. To simulate the trajectories, the assumptions made are that no lockdown on March 26 allows migrants and others to travel back home. Due to this, a small number of infected people reach many parts of the country. We incorporated this by reducing the value of 1/ε, the parameter that captures reach of the disease to different regions, to phase 3 value (1108) on April 5 instead of phase 2 value (20,850). There appeared a sharp reduction; however, even with strict lockdown, this reduction happened within a month. Hence, it is reasonable to assume that with no lockdown and significant travel all over, the reduction would have happened much faster.
For no lockdown, the value of 1/ε continues to reduce sharply, the value of 1/β remains as in phase 1 and γ and η improve slowly. Lockdown from April 1 changes the value of 1/β to phase 2 value. However, the value of γ improves slowly. This is because, with rapid spread of infection, hospitals start getting overwhelmed (as would be evident from plots below) and so are unable to ramp up and put in place a proper setup. Lockdown from May 1 reduces values of 1/ε further and values of 1/β; γ and η improve but slowly.
If there was no lockdown, the number of active infections would have peaked at 14+ million and the peak would have arrived by mid-May. This would have resulted in overwhelming our hospitals and widespread panic. There was little qualitative difference between two lockdown timings of April 1 and May 1, 2020. These would have resulted in a peak between 0 and 5 million active infections by mid-June. If there was no lockdown, it would have resulted in more than
2 million deaths. The two lockdowns (April 1 and May 1, 2020) would have resulted in between 0.5-1 million deaths. The number of deaths with current trends is projected to be less than 0.2 million.
The major limitation of our model was the non-availability of accurate data. While our modelling efforts have tried to match the available data, the utility of our projections is limited by the accuracy and reliability of the available data. Hence, more accurate and reliable input data to the model would result in more reliable projections.
In conclusions the new model is useful to calculate the possible disease burden in a realistic manner and may help plan for the resources. The impact of lockdown and interventions undertaken in a timely manner has been highlighted. The continuing importance of interventions such as use of mask, hand hygiene and physical distance is reinforced.
Financial support & sponsorship: None.
Conflicts of Interest: None.
| References|| |
Ministry of Health & Family Welfare, Government of India. Available from: https://www.mohfw.gov.in
, accessed on September 25, 2020.
Changotra R, Rajput H, Rajput P, Gautam S, Arora AS. Largest democracy in the world crippled by COVID-19: Current perspective and experience from India. Environ Dev Sustain
Weinberger DM, Cohen T, Crawford FW, Mostashari F, Olson D, Pitzer VE, et al
. Estimating the early death toll of COVID-19 in the United States. bioRxiv
2020. doi: 10.1101/2020.04.15.20066431.
Block P, Hoffman M, Raabe IJ, Dowd JB, Rahal C, Kashyap R, et al
. Social network-based distancing strategies to flatten the COVID-19 curve in a post-lockdown world. Nat Hum Behav
Leng T, White C, Hilton J, Kucharski A, Pellis L, Stage H, et al
. The effectiveness of social bubbles as part of a Covid-19 lockdown exit strategy, a modelling study. Wellcome Open Res
Mitze T, Kosfeld R, Rode J, Wälde K. Face masks considerably reduce COVID-19 cases in Germany. medRxiv
2020. doi: 10.1101/2020.06.21.20128181.
Balachandar S, Zaleski S, Soldati A, Ahmadi G, Bourouiba L. Host-to-host airborne transmission as a multiphase flow problem for science-based social distance guidelines. Int J Multiphas Flow
Cheng VC, Wong SC, Chan VW, So SY, Chen JH, Yip CC, et al
. Air and environmental sampling for SARS-CoV-2 around hospitalized patients with coronavirus disease 2019 (COVID-19). Infect Control Hosp Epidemiol
Chakraborty T, Ghosh I. Real-time forecasts and risk assessment of novel coronavirus (COVID-19) cases: A data-driven analysis. Chaos Solitons Fractals
Mahajan P, Kaushal J. Epidemic trend of COVID-19 transmission in India during lockdown-1 phase. J Community Health
Guan D, Wang D, Hallegatte S, Davis SJ, Huo J, Li S, et al
. Global supply-chain effects of COVID-19 control measures. Nat Hum Behav
Bateman IJ, Dannenberg A, Elliott R, Finus M, Koundouri P, Millock K, et al
. Perspectives on the Economics of the Environment in the Shadow of Coronavirus. Environ Resour Econ (Dordr)
Britton T, Ball F, Trapman P. A mathematical model reveals the influence of population heterogeneity on herd immunity to SARS-CoV-2. Science 2020; 369
Lee S, Kim T, Lee E, Lee C, Kim H, Rhee H, et al
. Clinical course and molecular viral shedding among asymptomatic and symptomatic patients with SARS-CoV-2 infection in a community treatment center in the Republic of Korea. JAMA Intern Med
2020. doi: 10.1001/jamainternmed.2020.3862.
Gomes MGM, Corder RM, King JG, Langwig KE, Souto-Maior C, Carneiro J, et al
. Individual variation in susceptibility or exposure to SARS-CoV-2 lowers the herd immunity threshold. medRxiv
2020. doi: 10.1101/2020.04.27.20081893.
McAloon C, Collins Á Hunt K, Barber A, Byrne AW, Butler F, et al
. Incubation period of COVID-19: A rapid systematic review and meta-analysis of observational research. BMJ Open
Rogers TN, Rogers CR, VanSant-Webb E, Gu LY, Yan B, Qeadan F. Racial disparities in COVID-19 mortality among essential workers in the United States. World Med Health Policy
2020. doi: 10.1002/wmh3.358.
Swaminathan A, Subramanian SV. Reflections on designing population surveys for COVID-19 infection and prevalence. Geroscience
2020. doi: 10.1007/s11357-020-00253-6.
Asnumali S, Kaushal S, Kumar A, Prakash M, Vidyasagar M. Modelling a pandemic with asymptomatic patients, impact of lockdown and herd immunity, with applications to COVID-19. Annu Rev Control
2020. doi: https://doi.org/10.1016/j.arcontrol.2020.10.003.
Robinson M, Stilianakis NI. A model for the emergence of drug resistance in the presence of asymptomatic infections. Math Biosci
[Figure 1], [Figure 2], [Figure 3], [Figure 4]
[Table 1], [Table 2]