Social distancing mediated generalized model to predict epidemic spread of COVID-19.

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.

Introduction

The prevention from coronavirus (COVID-19), all over the world, has become an utmost challenge for mankind [1-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-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-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-26] and recent pandemics, like COVID-19 [27-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-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.

Social Distancing Mediated Growth Model.

The phenomenological models mainly consist of generalized growth (GG) model , which at yields into simple exponential growth (EG) model, generalized Richard growth (GRG) model , which at converts to generalized logistic growth (GLG) model and further at yields in classical logistic growth (CLM) model [42, 43]. Here, () corresponds to the time derivative and is the cumulative number of confirmed cases in any network at time t. accounts for the exponential profile for growth trajectory, which will be constant at , sub-exponential at and exponential at . 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 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, 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 as social distancing parameter, the value of which can be ranged from to 1, as described in Table 1. Here corresponds to the zero distance between nodes leading to maximum spread scenario. 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. 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 . Finally, corresponds to the quarantine situation and there is no chance of social interaction. Both 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 means, we by ourselves share to all end-users and for means, we immediately shutdown the whole internet.

Table 1

Social distancing parameter categorization

Stages	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.
Forced confinement (Lockdown)	0>L_D≥-1	0>L_D>-1	Forced Confinement (lockdown) scenario where the social interactions are prohibited. The strength of forced confinement corresponds to the variation of L_D from 0 and -1.
		L_D=-1	Extreme lockdown (quarantine) scenario where the network nodes are so distant from one another that infection cannot spread.

The categorization of 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

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

Conventional phenomenological models and modified growth 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-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 shows almost perfect regressive growth patterns with respect to COVID-19 data of China, which, at slightly higher social distancing parameter (blue dashed curve), shows notably delayed growth to the saturated state.

Fig. 1

Conventional models in comparison with modified Richard model. a The simulated of cumulative growth versus time for comparison between EG (magenta), GG (green), GLG (purple), GRG (orange) and MRG (blue) models with respect to the data of real confirmed COVID-19 cases in China (CN), red curve. It can be seen that the MRG model shows more accurate fit to the real data at given parameters. b The comparison between cumulative growth modified models (MEG (magenta), MGG (green), MLG (brown) and MRG (blue)) with respect to the real COVID-19 data of China, red curve. Here and are for solid and dashed curves, respectively. It is illustrated that, even with modified daily growth rate r(t), the MRG model shows more compelling results with respect to COVID-19 growth of China. The remaining parameters are with , , and (with respect of COVID-19 data of China)

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. 1b. The interesting thing is that not only delays the growth of MRG model, but it also reduces the speed of aggression in all modified models, which further supports our concept of social distancing induced delay in phenomenological model.

The both and are crucial for the calibration of modified growth models, as illustrated in Fig. 2, where the detailed dynamics of MRG, as a function of and , versus time are shown. In a high social distancing scenario (i.e. ), the cumulative number (or total infections) moves from rapid growth to delayed, or slow, growth trends with higher to lower initial growth rate , respectively (see Fig. 2a), which is obvious. However, in lockdown regime (i.e. ), the will decrease from initial state, which is , to the zero infection scenario, as can be seen in Fig. 2b. The interesting thing is that it will decay faster at higher as compared to low values of . It is because the 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 , under influence of , can be seen in Fig. 2c and d, with and , respectively. The transforms from fast growth regime to fast decay regime by changing from to , which get more rapid with increase in , see Fig. 2d. It reveals that in growth regime (i.e. ), will reach sooner to the maximum network capacity K at higher values of and in decay regime (i.e. ), will decay sooner to zero infection state from the initial number of infections, at lower values of . At , will not be modified and will remain equivalent to the initial number of infections .

Fig. 2

The dynamics of MRG model with respect to social distancing parameter and initial daily growth rate . a and b are versus time for various (from 1 to 200), at and , respectively, where the total network size and initial number of confirmed patients . It is illustrated that the cumulative growth , in both social distancing and forced confinement domain , crucially depends on the initial daily rate . The higher initial rate causes faster growth (decay) of in () regimes. c) and (d are versus time for ranging from to 1, at and , respectively, where and . It further strengthens the argument of (a) and (b), regarding the effects of initial rate on in presence of . Remaining parameters used in numerical calculations are same as in Fig. 1

Robust influence of social distancing on cumulative growth of COVID-19.

The influence of social distancing is illustrated in Fig. 3, where the multiple outcomes of and r(t) are predicted with various calibrations for , with respect to the real COVID-19 data of China, USA and India. With a small increase in , grows rapidly and by strengthening the ratio of growth decreases, as illustrated in Fig. reffig3a, b and c. The blue curves (at for China, for USA, and 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 , the cumulative growth gets significantly enhanced (rapid) or damped (delayed), respectively, at a fixed arbitrary network size K. These modifications in 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 , providing the evidence for the reliability of our model.

Fig. 3

Social distancing induced robust decrease in COVID-19 cumulative growth. The MRG model versus various , for growth , (a), (b) and (c), and growth rate r(t), (d), (e), and (f), with respect to the real COVID-19 data of confirmed cases in China (CN), USA and India (IND), respectively. The real data of COVID-19 is represented by the red curves. Magenta, brown, blue, orange, and purple curves correspond to and 0.02 in (a) and (d), and in (b) and (e), and and 0.009 in (c) and (f), respectively. and for (a) and (d), and for (b) and (e), and and for (c) and (f). It can be observed that at MRG model almost exactly follows COVID-19 growth in China at particular initial rate . But in the case of COVID-19 in India and USA, the MRG model, so-far, is showing fit to real at lower values of but initial rates are very high, as compared to the case of China, yielding in rapid growth of COVID-19. Thus to deal with the rapid growth of COVID-19, especially in the case of USA, a more effective policies for social distancing are needed. The other parameters used in calculations are same as illustrated in Fig. 1

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 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 . With higher and K, a small change in will notably affect the growth trajectory, but for smaller and K, a strong change in will be needed in order to produce similar influences. In this way, our model provides more flexibility to map and predict real-time trajectories, see supplementary materials [44] for details. Not only this, but also by observing different trends for , with respect to , one can predict the future growth patterns with change in current social distancing scenarios. For example, a slight increase in to and to 0.031, will increase the daily rate to and to for USA and India, respectively. Similarly, a small decrease in yields in significant decrease in r(t), as illustrated in Fig. 3d, e and f.

Predicting outcomes of COVID-19: Daily growth rate in forced confinement scenarios

The rapid diffusion in the growth spectrum arrives in lockdown regime, , as illustrated in Fig. 4, where we used MRG model to predict spectrum of outcomes for r(t) versus various , by considering the peak value as an initial state. In the case of China, peak value of r(t) (that we consider) occurs around February. At that peak rate, which is around , the MRG model, at , 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. 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 was varied from to , the COVID-19 growth tends could have extended from February to April, even in the lockdown domain. Further increase in , in regime, could have led the country to worse scenarios, as illustrated in Fig. 4a.

Fig. 4

Predicting the outcomes of COVID-19 daily growth rate in forced confinement domain. The growth rate r(t) versus time with possible outcomes of MRG model from the state of maximum daily growth rate of COVID-19 confirm cases in China (a), USA (b), and India (c) in social distancing and forced confinement domain . Here with for (a), with for (b) and with for (c). The red curves correspond to real COVID-19 data. Black, magenta, brown, green, blue, orange, and purple curves correspond to and in (a), and in (b), and and in (c), respectively. It is observed that the increase in towards zero, or decrease in forced confinement policies, leads to the more rapid growth of daily COVID-19 cases, while decrease in , 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

In the case of USA, so-far, the maximum is occurring around July, which can decay to zero around 15th August if it follows the decay trajectory at . But it can be extended to 1st November with an increase in to , or it can continue to increase in less social distancing regime, i.e. , as can be seen in Fig. 4b. Similarly, for India, where the growth rate r(t) is still increasing, the maximum till 8th August can tend to decay with , and can continue rapid increase with , in spite implementation of social distancing. In lockdown regime, at , the infection growth rate can be zero around 20th September, but with a slight relaxation in lockdown with a small increase in can prolong the spread till 22th December, see Fig. 4c. Further increase in can yield in rapid and frightening spread of COVID-19. Thus, to stop the spread, policies with higher social distancing factor are needed.

Conclusion

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 (), the cumulative growth as well as growth rate appears to be delayed with , but it will still be increasing. However, in strong distancing or in lockdown regime (), 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.