key: cord-0685482-5e1qwebw
authors: Yasir, Kashif Ammar; Liu, Wu-Ming
title: Social distancing mediated generalized model to predict epidemic spread of COVID-19
date: 2021-04-11
journal: Nonlinear Dyn
DOI: 10.1007/s11071-021-06424-0
sha: c39fb12d18815b2cef9d6d87261da0246166a282
doc_id: 685482
cord_uid: 5e1qwebw

The extensive proliferation of recent coronavirus (COVID-19), all over the world, is the outcome of social interactions through massive transportation, gatherings and population growth. To disrupt the widespread of COVID-19, a mechanism for social distancing is indispensable. Also, to predict the effectiveness and quantity of social distancing for a particular social network, with a certain contagion, a generalized model is needed. In this manuscript, we propose a social distancing mediated generalized model to predict the pandemic spread of COVID-19. By considering growth rate as a temporal harmonic function damped with social distancing in generalized Richard model and by using the data of confirmed COVID-19 cases in China, USA and India, we find that, with time, the cumulative spread grows more rapidly due to weak social distancing as compared to the stronger social distancing, where it is explicitly decreasing. Furthermore, we predict the possible outcomes with various social distancing scenarios by considering highest growth rate as an initial state, and illustrate that the increase in social distancing tremendously decreases growth rate, even it tends to reach zero in lockdown regimes. Our findings not only provide epidemic growth scenarios as a function of social distancing but also provide a modified growth model to predict controlled information flow in any network.

The prevention from coronavirus , all over the world, has become an utmost challenge for mankind [1] [2] [3] [4] [5] [6] . It has cost hundreds of thousands of lives and affected millions, and still growing with its ruthless nature [7] . In 14th century, the bubonic Black Death took nine years to proliferate across Europe because it was restricted by the average 1.5 km daily travel dispersion of that time [8] . However, the recent COVID-19 has no such limitations and it just took days to spread all over the world. The reason behind such frightening difference is in the increase in social interactions stimulated by massive transportation networks, gatherings and overall population growth, which led to an aggressive spread of COVID-19 [6, 9, 10] . The only way to minimize social interactions is to implement social distancing mechanism, or policies, and, more importantly, closing or decreasing mobility in transportation networks-especially in social networks where the contagion factor is higher [11] [12] [13] [14] .

Over the years, the research on network epidemiology [15] has enabled us to analyze the growth of pandemic processes by developing a wide range of numerical methods and models [16] [17] [18] [19] . These models notably include SIR [20] , SEIR [21] , and neural network-based control models [22] , which are validated for the analysis and short-term prediction of various previous [23] [24] [25] [26] and recent pandemics, like COVID-19 [27] [28] [29] [30] [31] [32] . However, the mechanistic models are quite fragile because of the dependence on underlying spatio-temporal microscopic processes, which are often missing in ongoing metapopulation pandemics [33, 34] .

Meanwhile, the phenomenological processes-which are more robust in calibration with empirical data [35] -like exponential, generalized exponential, logistic and Richard growth models have been validated to provide the dynamical analysis of the pandemic, like Black Death, Ebola virus, influenza, Zika virus, SARS, MERS [36] [37] [38] [39] , even for COVID-19 [40, 41] . Among them, generalized logistic and Richard models, in spite of crucial dependence on enough empirical data, provide more reliable analysis than exponential models, because of the dynamical saturation at the later stages of pandemic [42, 43] . Although, the phenomenological models have been sufficiently tested, but there is an absence of a standardized model to incorporate the element of social distancing with respect to the growth rate. Thus, to successfully design, apply and monitor social distancing mechanism, an effective as well as generalized mathematical model is needed to predict the usefulness of the social distancing for local (national) or non-local (global) in less or more contagious social networks.

In this manuscript, we proposed a social distancing induced generalized model to predict the growth of COVID-19. By modeling daily growth rate as a time dependent harmonic function, which is damped with the distance between each entity and its surrounding entities in a social network, we study the effects of social distancing (or lockdown) on the cumulative growth of the spread. We illustrate the comparison between state of the art growth mathematical models with and without modified growth rate and find that the modified Richard model, with harmonically modified growth rate, provides the most reliable growth predictions. By using modified Richard model, in calibration with empirical data of China, USA and India, we show that the daily growth rate as well as the total cumulative number of confirmed cases are (or could) increasing at weak social distancing factor in each network. However, the higher values of social distancing are significantly reducing the growth rate, even can be reduced to zero in lockdown regime.

The phenomenological models mainly consist of generalized growth (GG) modelṄ c (t) = r N c (t) p , which at p = 1 yields into simple exponential growth (EG) model, generalized Richard growth (GRG) model

, which at α = 1 converts to generalized logistic growth (GLG) model and further at p = 1 yields in classical logistic growth (CLM) model [42, 43] . Here, (˙) corresponds to the time derivative and N c (t) is the cumulative number of confirmed cases in any network at time t. p ∈ [0, 1] accounts for the exponential profile for growth trajectory, which will be constant at p = 0, sub-exponential at 0 < p < 1 and exponential at p = 1. K accounts for the total asymptotic capacity, or total number of nodes in any social network. α is a scaling parameter and accounts for any asymmetric deviation of logistic trajectories. In all these models, the most crucial parameter is r , which defines the growth rate or the number of infections per day.

The growth rate r , or the quantity of spread (flow) per unit time, is the main ingredient to measure the cumulative amount of spread in any network, irrespective of mechanistic approach or phenomenological methods. In all of the previous models that are designed to measure the collective magnitude of the spread, the growth rate is used as a constant parameter. But, in actual, growth rate is the factor that varies with time depending upon the underlying microscopic network parameters, especially for meta-populated contagious network where epidemic spread is needed to be measured with respect to social distancing. Further, in accordance with the empirical data, daily growth rate r mostly follows time dependent nonlinear trends rather than linear trajectories. Therefore, in order to incorporate time (daily) fluctuation in growth rate, we consider r as a Parametric domains for L D Subdomains Policy description Self-induced social distancing +1 ≥ L D ≥ 0 L D = +1 Zero distance between network nodes leading to maximum spread scenario. +1 > L D ≥ 0 Self-induced social distancing protocols exist but there is no force confinement. The level of adaptation of protocols by individuals corresponds to L D variation between +1 and 0. The categorization of L D is based on the potential social distancing and forced confinement factors, like social gatherings, social interactions in transportation and educational institutes. It will be based on the policies adopted by the governments time dependent harmonic function, which is damped or perturbed over the distance between any node and its neighboring node(s),

where, K is the total number of nodes (entities) in the network. r (t) is total growth rate of the system defined over the transmission rate of each entity of network to its neighboring entities

Similarly, L D is the effective distance between the network nodes, which is defined over the distance (or the strength of social interactions) of each entity in the network with its surrounding entities

We named L D ∈ [−1, 1] as social distancing parameter, the value of which can be ranged from −1 to 1, as described in Table 1 . Here L D = 1 corresponds to the zero distance between nodes leading to maximum spread scenario. 1 > L D > 0 means that there is social distancing but still the infection can be transferred. It can be related to a situation where protocols for social distancing are present but there is no lockdown. 0 > L D > −1 represents a scenario of lockdown, means there is a considerable distance between each entity. In this case, the spread rate, instead of increasing, will decrease to zero from the initial infections state N c (0) = N 0 . Finally, L D = −1 corresponds to the quarantine situation and there is no chance of social interaction. Both L D = ±1 are ideal scenarios and in analogy, one can consider that there is a video which is going to be viral on internet and we want to stop it. Here L D = 1 means, we by ourselves share to all end-users and for L D = −1 means, we immediately shutdown the whole internet.

We consider the particular solution for second-order harmonic differential equation for r (t), which yield in a hyperbolic function r (t) = r 0 sinh(L D t), where r (0) = r 0 is the initial growth rate at t = 0. With perturbed r (t), the modified generalized growth (MGG) model will be read as N c (t) = r 0 sinh(L D t)N c (t) p , which, similarly at p = 1 yields in modified exponential growth (MEG) model. The modified Richard growth (MRG) model can be expressed as,

which again at α = 1 yields in modified logistic growth (MLG) model.

The comparison between EG, GG, GLG, GRG and MRG models is illustrated in Fig. 1a and comparison between modified models (MEG, MGG, MLG and MRG) is presented in Fig. 1b , by using data of cumulative COVID-19 growth in China as an example. (Note: The real COVID-19 data used in the manuscript was adopted from Refs. [1] [2] [3] [4] [5] , till August 8, 2020.) It can be clearly seen that EG model follows a rapid exponential growth (magenta curve), as expected, and GG model follows a (sort-of) sub-exponential growth (green curve) because of p, see Fig. 1a . On the other hand, both GLG (purple curve) and GRG (orange curve) models follow delayed cumulative growth patterns, but still they are incomparable with the actual growth in China at given parameters. However, MRG model (blue curve), at L D = 0.12 shows almost perfect regressive growth patterns with respect to COVID-19 data of China, which, at slightly higher social distancing parameter L D = 0.08 (blue dashed curve), shows notably delayed growth to the saturated state. Similarly, in the comparison between modified models (MEG, MGG, MLG and MRG), one can clearly note that the MRG model still provides a more reliable fit to the real data at given parameters, see Fig. 1 and r 0 , versus time are shown. In a high social distancing scenario (i.e. L D = 0.01), the cumulative number N c (t) (or total infections) moves from rapid growth to delayed, or slow, growth trends with higher to lower initial growth rate 1 ≤ r 0 ≤ 200, respectively (see Fig.  2a ), which is obvious. However, in lockdown regime (i.e. L D = −0.001), the N c (t) will decrease from initial state, which is N 0 = 500, to the zero infection scenario, as can be seen in Fig. 2b . The interesting thing is that it will decay faster at higher r 0 as compared to low values of r 0 . It is because the N c (t) will follow same initial rate to decay as it is following in growth scenario. It can be understood through an example that if we are driving a car at high speed and we suddenly apply breaks, then we will feel more deceleration as compared to the case when we apply breaks at a lower speed. More detailed behavior of N c (t), under influence of L D , can be seen in Fig. 2c and d, with r 0 = 20 and r 0 = 60, respectively. The N c (t) transforms from fast growth regime to fast decay regime by changing L D from +1 to −1, which get more rapid with increase in r 0 , see Fig. 2d . It reveals that in growth regime (i.e. 1 ≥ L D > 0), N c (t) will reach sooner to the maxi-mum network capacity K at higher values of L D and in decay regime (i.e. 0 > L D ≥ −1), N c (t) will decay sooner to zero infection state from the initial number of infections, at lower values of L D . At L D = 0, N c (t) will not be modified and will remain equivalent to the initial number of infections N 0 .

The influence of social distancing L D is illustrated in Fig. 3 , where the multiple outcomes of N c (t) and r (t) are predicted with various calibrations for L D , with respect to the real COVID-19 data of China, USA and India. With a small increase in L D , N c (t) grows rapidly and by strengthening L D the ratio of growth N c (t) decreases, as illustrated in Fig. reffig3a , b and c. The blue curves (at L D = 0.12 for China, L D = 0.61 × 10 −3 for USA, and L D = 0.022 for India) indicate the growth of cumulative number of infections matching with the real data. But when we slightly increase or decrease social distancing parameter L D , the cumulative growth N c (t) gets significantly enhanced (rapid) or damped (delayed), respectively, at a fixed arbitrary network size K . These modifications in N c (t) can also be validated from the spectrum of r (t), shown in Figs. 3d, e and f, where blue curves provide almost a perfect fit to the real growth rate of infections. By changing conditions for social distancing, the r (t) will either reach rapidly or slowly to the maximum network size with increase or decrease in L D , providing the evidence for the reliability of our model. It should be noted that the narrow and wide time spans, or the so-called standard deviation, correspond to the faster and slower approach to the maximum network size, which we considered fix and arbitrary (except for China because we have the state of total infections) in order to discuss rapid growth patterns. But it can be changed with the real social/sub-social network sizes. Further, one can question the notable variations in the quantitative values of L D for network of China, USA, and India. It is not only because of the different network sizes K , but it is also because of crucial role of r 0 . With higher r 0 and K , a small change in L D will notably affect the growth trajectory, but for smaller r 0 and K , a strong change in L D will be needed in order to produce similar influences. In this way, our model provides more flexibil-ity to map and predict real-time trajectories, see supplementary materials [44] for details. Not only this, but also by observing different trends for N c (t), with respect to L D , one can predict the future growth patterns with change in current social distancing scenarios. For example, a slight increase in L D = 0.61 × 10 −3 to 1.5 × 10 −3 and L D = 0.022 to 0.031, will increase the daily rate r (t) ≈ 80 × 10 3 /day to 120 × 10 3 /day and r (t) ≈ 100 × 10 3 /day to 130 × 10 3 /day for USA and India, respectively. Similarly, a small decrease in L D yields in significant decrease in r (t), as illustrated in Fig. 3d , e and f.

The rapid diffusion in the growth spectrum arrives in lockdown regime, L D < 0, as illustrated in Fig.  4 , where we used MRG model to predict spectrum of outcomes for r (t) versus various L D ∈ [−1, 1], by considering the peak value as an initial state. In the case of China, peak value of r (t) (that we consider) occurs around 5 th February. At that peak rate, which is around r (t) = 4000/day, the MRG model, at L D = −0.007, exactly follows the decay trend of COVID-19 in China, which illustrates the effectiveness of strict lockdown policies of China. Although, ideal scenarios, i.e. L D = ±1 are realistically impossible, but it can be seen that the scenario for COVID-19 dynamics in China could have been worse (or better) with relaxation (or increase) in lockdown and social distancing policies. For example, if L D was varied from −0.007 to −0.001, the COVID-19 growth tends could have extended from ≈ 25 th February to ≈ 5 th April, even in the lockdown domain. Further increase in L D , in L D > 0 regime, could have led the country to worse scenarios, as illustrated in Fig. 4a .

In the case of USA, so-far, the maximum r (t) ≈ 80000/day is occurring around ≈ 20th July, which can decay to zero around 15th August if it follows the decay trajectory at L D = −0.001. But it can be extended to 1st November with an increase in L D to −0.0001, or it can continue to increase in less social distancing regime, i.e. L D > 0, as can be seen in Fig.  4b . Similarly, for India, where the growth rate r (t) is still increasing, the maximum r (t) ≈ 650000/day till 8th August can tend to decay with L D < 0, and can continue rapid increase with L D > 0, in spite imple- It is observed that the increase in L D towards zero, or decrease in forced confinement policies, leads to the more rapid growth of daily COVID-19 cases, while decrease in L D , or more strict force confinement, leads to the robust decay in daily COVID-19 cases. Remaining parameters used in numerical calculation are same as in Fig. 1 mentation of social distancing. In lockdown regime, at L D = −0.05, the infection growth rate can be zero around 20th September, but with a slight relaxation in lockdown with a small increase in L D = −0.007 can prolong the spread till 22th December, see Fig. 4c . Further increase in L D can yield in rapid and frightening spread of COVID-19. Thus, to stop the spread, policies with higher social distancing factor L D are needed.

To conclude, we propose a social distancing mediated generalized model to predict the spread of COVID-19, with respect to social distancing and lockdown regimes. By considering growth rate as a time variant function of social distancing in generalized Richard model, we illustrate different possible outcomes of cumulative growth (total number of infections) and growth rate (infection per day) in different social distancing and lockdown scenarios. In weak distancing regime (1 > L D > 0), the cumulative growth as well as growth rate appears to be delayed with L D > 0, but it will still be increasing. However, in strong distancing or in lockdown regime (0 > L D > −1), the infection spread will not only stop but will also decay from its initial state with time. Further, by using our model, we predict different future scenarios (especially for the COVID-19 spread in USA and India) with highest growth rate as initial state and show that, in absence of strict lockdown policies, the COVID-19 pandemic spread will be worst and prolonged. Thus, in order to contain the spread, in any contagion network, strong social distancing polices are necessary. Furthermore, our model is not only effective in predicting the influence of social distancing for epidemic spreads, but it can also be validated for controlling information flow in any network. 

The authors declare no competing financial interests.

Numerical methods and codes In order to govern temporal evolution of modified Richard growth model, in calibration with the empirical data, we use the NDSolve function of Mathematica 10 programming application (https://www.wolfram.com /mathematica/) to solve equation (2) by considering given parameters in the manuscript. After generating data, we use graphical application OriginPro 9.0 (https://www.originlab. com/) to illustrate and compare computed results with the empirical data.

Appendix: The influence of r 0 on the cumulative growth N c (t) and growth rate r(t) Figure 5 illustrates the cumulative growth N c (t) and growth rate r ( t) trends for various network sizes K with (a) (b) Fig. 5 (a) N c (t) and (b) r ( t) as function of time for various initial growth rates r 0 and total network size K , for fix social distancing L D = 0.023. Red curves correspond to real COVID-19 data of India (IND), and Magenta, brown, blue, orange, and purple curves corresponds to r 0 = 9 with K = 80×10 6 , r 0 = 10 with K = 20 × 10 6 , r 0 = 12 with K = 10 × 10 6 , r 0 = 15 with K = 6 × 10 6 and r 0 = 16 with K = 4 × 10 6 , respectively. Remaining parameters are same as in Fig. 1 of main text different initial growth rates r 0 , at fixed social distancing scenario L D = 0.023, by considering COVID-19 data of India. It can be seen that, with higher initial growth rate r 0 = 16, the network size gets shrunk to K = 4 × 10 6 , and for lower rate r 0 = 16, the total network size increases to K = 80 × 10 6 , as can be seen in Fig. 5a . It is because the growth trajectory is fixed by the applied social distancing parameter L D , means social distancing is using same rate to stop the spread as it was initial state r 0 , as can also be seen in Fig. 2 of main text. Therefore, higher r 0 will be less effective as compared to the lower rates. It can be interpreted that the social distancing L D is more effective for higher initial rates in comparison to lower rates. Similarly, growth rate r (t) takes more time to grow and decay at lower r 0 as compared to the higher initial rate r 0 , as can be seen Fig. 5b .

COVID-19) pandemic, World Health Organization (WHO), United Nations

Wikipedia: Coronavirus disease 2019

United States

Chinese Center for Disease Control and Prevention (CCDC)

The effect of travel restrictions on the spread of the 2019 novel coronavirus (COVID-19) outbreak

COVID-19), dashboard, World Health Organization (WHO), United Nations

The bubonic plague

How will country-based mitigation measures influence the course of the COVID-19 epidemic?

Effect of the social distancing measures on the spread of COVID-19 in 10 highly infected countries

Lockdown, one, two, none, or smart Modeling containing covid-19 infection

Quantifying the impact of physical distance measures on the transmission of COVID-19 in the UK

Local decision making for implementing social distancing in response to outbreaks

Impacts of social and economic factors on the transmission of coronavirus disease 2019 (COVID-19) in China

Epidemic processes in complex networks

Epidemics and percolation in small-world networks

Mathematical models for COVID-19: applications, limitations, and potentials, Public Health Emerg

Tight lower bound for percolation threshold on an infinite graph

Percolation on sparse networks

Individual vaccination as nash equilibrium in a SIR model with application to the 2009-2010 Influenza A (H1N1) epidemic in France

Statistical inference in a stochastic epidemic SEIR model with control intervention: Ebola as a case study

Neural Network aided quarantine control model estimation of global Covid-19 spread

Modeling the SARS epidemic

Effect of social group dynamics on contagion

Reconstructing the 2003/2004 H3N2 influenza epidemic in Switzerland with a spatially explicit, individual-based model

Global disease spread: statistics and estimation of arrival times

Modelling transmission and control of the COVID-19 pandemic in Australia

The effect of control strategies to reduce social mixing on outcomes of the COVID-19 epidemic in Wuhan, China: a modelling study

Risk assessment of novel coronavirus COVID-19 outbreaks outside

COVID-19 outbreak on the diamond princess cruise ship: estimating the epidemic potential and effectiveness of public health countermeasures

A SIR model assumption for the spread of COVID-19 in different communities

Analysis and Prediction of COVID-19 Pandemic in Pakistan using Time-dependent SIR Model

Predictability of catastrophic events: Material rupture, earthquakes, turbulence, financial crashes, and human birth

Computational irreducibility and the predictability of complex physical systems

A flexible growth function for empirical use

Using Phenomenological Models to Characterize Transmissibility and Forecast Patterns and Final Burden of Zika Epidemics

Fitting dynamic models to epidemic outbreaks with quantified uncertainty: A Primer for parameter uncertainty, identifiability, and forecasts

A generalizedgrowth model to characterize the early ascending phase of infectious disease outbreaks

A novel sub-epidemic modeling framework for short-term forecasting epidemic waves

Short-term Forecasts of the COVID-19 Epidemic in Guangdong and Zhejiang

COVID-19 and Italy: what next?

Generalized logistic growth modeling of the COVID-19 outbreak: comparing the dynamics in the 29 provinces in China and in the rest of the world

COVID-19 Pandemic in Pakistan: Stages and Recommendations, (preprint) medRxiv

See Supplemental Material at, web link, for the correlation between social distancing factor L D and initial growth rate r 0

Acknowledgements This work was supported by the National