Access the full text.

Sign up today, get DeepDyve free for 14 days.

Physics
, Volume 2020 (2012) – Aug 28, 2020

/lp/arxiv-cornell-university/spatial-interference-between-infectious-hotspots-epidemic-condensation-y2oV7gjoFz

- ISSN
- 0129-1831
- eISSN
- ARCH-3341
- DOI
- 10.1142/S0129183121500443
- Publisher site
- See Article on Publisher Site

December 24, 2020 1:39 manuscript_v2 Spatial interference between infectious hotspots: epidemic condensation and optimal windspeed Johannes Dieplinger Scuola Normale Superiore di Pisa, Piazza dei Cavalieri, Pisa, Italy johannes.dieplinger@t-online.de Sauro Succi Center for Life Nanosciences at La Sapienza, Italian Institute of Technology, Roma, Italy, Scuola Normale Superiore di Pisa, Piazza dei Cavalieri, Pisa, Italy, Physics Department, Harvard University, Cambridge, USA, sauro.succi@sns.it We discuss the eects of spatial interference between two infectious hotspots as a function of the mobility of individuals (wind speed) between the two and their relative degree of infectivity. As long as the upstream hotspot is less contagious than the downstream one, increasing the wind speed leads to a monotonic decrease of the infection peak in the downstream hotspot. Once the upstream hotspot becomes about between twice and ve times more infectious than the downstream one, an optimal wind speed emerges, whereby a local minimum peak intensity is attained in the downstream hotspot, along with a local maximum beyond which the bene cial eect of the wind is restored. Since this non-monotonic trend is reminiscent of the equation of state of non-ideal uids, we dub the above phenomena "epidemic condensation". When the relative infectivity of the upstream hotspot exceeds about a factor ve, the bene cial eect of the wind above the optimal speed is completely lost: any wind speed above the optimal one leads to a higher infection peak. It is also found that spatial correlation between the two hotspots decay much more slowly than their inverse distance. It is hoped that the above ndings may oer a qualitative clue for optimal con nement policies between dierent cities and urban agglomerates. Keywords : population dynamics, SIR model, interference, nonlinearity, pattern formation PACS Nos.: 87.23.Cc, 01.75.+m, 02.70.Bf 1. Introduction In early 2020 the world has been taken by a very aggressive global pandemic, the covid-19, which spread around the entire planet at unprecedented speed in mankind history. As we speak, the pandemic, originated in Wuhan, China, allegedly in Jan- uary 2020 has spread out over 100 countries, with over ten million contagion cases and over 500,000 casualties worldwide, as of end June 2020, giving rise to the trade-o between strong con nement measures to stave o devastating eects on 2,3 4{7 health systems and economic, social and psychological terms. A distinctive feature of the covid-19 pandemia is the heterogeneity of the viral arXiv:2012.12077v2 [physics.soc-ph] 23 Dec 2020 December 24, 2020 1:39 manuscript_v2 infection in space; in many countries a large fraction of the overall infection counts originated from very speci c hotspots, such as Lombardy in Italy and NYC in the USA. This strong inhomogeneity calls, among others, for a proper modelling of the mechanism by which the infection propagates in space and time. In this paper we address this issue by isolating a toy-problem, namely the inter- action between two infectious hotspots sitting at two separate locations in space. Special attention is paid to the spatial interference between the two hotspots, in particular the way that the presence of the rst aects the viral evolution in the second, depending on the mobility and infection rates in the two hotspots. To this purpose, spatial mobility is described by a simple advection-diusion SIR (ADSIR) model, in which diusion encodes small-distance mobility (say walking), while advection stands for mid-range mobility (say train or car driving). Rather than being an exact model of these real-world mobility schemes { for which network models would certainly be a more accurate choice { diusion and advection have the general purpose to realize two dierent mobility mechanisms with a dierent scaling in time and space. Even though human mobility proceeds by more complex 10{13 mechanisms than AD, typically encoded by mobility networks, the present AD- SIR model exposes nonetheless a number of interesting qualitative features related the spatio-temporal coupling between the two infectious hotspots. 2. Mathematical formulation We describe the covid-19 pandemic by means of a standard SIR model { which 7,11,15{20 is the starting point for numerous interesting models in epidemiology { coupled in space via an advection-diusion mechanism: @ s = r(Us + Drs) si (1) @ i = r(Ui + Dri) + si i (2) @ r = r(Ur + Drr) + i (3) where sir(x; y; t) is the population of Susceptible, Infected and Recovered indi- viduals at position x; y and time t, respectively. The coecients ; correspond to infection and recovery rate, respectively. In the above equations U is the wind speed, which we take aligned with the x-axis without loss of generality and D is the diusivity. The total number of individuals of species k = s; i; r at a given time t is thus given by integral over the entire domain of the corresponding densities: N (t) = n (x; y; t)dxdy; k = s; i; r; i = 1; 2 where the integral runs over the k;h k hotspot regions HS1 and HS2. The speed U will be compared to an important in- trinsic reference velocity: (s(t)=Nf )t 1 0 The infected population i(t) grows in HS1 as i(t) A e . When- ever the wind exceeds a critical speed U = w (s(t)=N f ) mitigation of c 1 epidemics is expected, due to the removal of infected individuals from the hotspot. December 24, 2020 1:39 manuscript_v2 Fig. 1. The simulation domain. Schematically shown is the simulation domain of length L. The boxes indicate the hotspots HS1 and HS2 with dierent contagion rate = f , where i i 0 0 is the basic contagion rate in the normal domain. The hotspots at distance d have width w in both directions. The homogeneous wind U is indicated by the black solid arrow and the red area indicates the small outbreak. For U=U > 1, we expect the rst hot-spot to become transparent and have little in uence on the second one. This of course largely depends on the a priori unknown parameters s(t)=N and the variable parameter f . For a reference, we take f = f = 50 and s(t)=N 0:5. 1 2 Hence, we de ne the reference speed U = w ( ). The main independent (dimensionless) parameters are then de ned as follows: The contagion rate is = 0:2 in the normal domain and = f in the hotspots 0 i i 0 i = 1; 2. The recovery rate is homogeneous = 0:15, such that the reproduction factor is R = f > 1. We x f = 50 while f ranges from 1 to 500, so that the i i 2 1 relative infectivity ratio f = f =f 2 (0:02; 10). The diusion coecient is D = 5. 1 2 The width of the hotspots is w = 10 and their distance is d = 50. The wind speed is measured in units of the reference speed, i.e. u = U=U . The grid spacing is 1 km and time is measured in days, corresponding to a reference speed U = 50 km=day and a diusivity D = 5 km =day. These are plausible scales for human mobility. We then study the solution of the ADSIR problem above as a function of the parameters u and f . In particular, we wish to assess under what conditions the presence of HS1 causes an increase of infections in HS2 in terms of both peak intensity and duration. 3. Simulation setup and results We set up two hotspots HS1 and HS2 of width w at position x = L=4 and x = 1 2 x + d + w respectively. The domain is a grid of size 64 1024. We place a small Gaussian outbreak at x = 0 and y = W=2 with a cuto at x = 10 in order to ensure that there are no initial infections in the hot-spots. The boundary conditions are chosen to be fully periodic. Fig. 1 shows the geometric set-up of the hot-spots in the domain. In Fig. 2 we show the typical evolution of the infected species in a single hotspot. The starting time of the hot-spot is de ned as the time t at which the infected s December 24, 2020 1:39 manuscript_v2 1.0 HS1 inf HS1 sus HS1 rec 0.8 0.6 0.4 0.2 0.0 0 100 200 300 400 500 600 time Fig. 2. Time evolution of the SIR populations in a single hotspot. The baseline parameters are: f = f = 5, = 0:5, U = 0, D = 0:05. The position of the hot-spot is at x = L=4. The red curve 1 2 shows the infected population in the hotspot, the black x's mark the start-up time, peak time and the decay time of the local epidemic, respectively. population reaches i(t ) = 0:002 N and the end time t at which the infected pop- s e ulation has dropped to i(t ) = 0:05 i . The time dierence t = t t can be e max e s de ned as the duration of the epidemic in this region. In Figure 3, we report a typical spatial interference pattern between HS1 on HS2. The gure clearly shows that the infected population generated in HS1 reaches up to HS2 and increases the local infection rate, thereby increasing the peak and possibly the duration as well. This is the typical scenario that HS2 policy makers endeavour to combat via lock-down measures. Peak and duration of the epidemics.{In Fig. 4, we summarize the eect of the wind and HS1 infectivity on the peak intensity of HS2. We measure the dimensionless peak value i =N in the second (downstream) max hotspot as a function of u and f , where N is the total number of individuals in the hotspot and i is the peak value of the infected population. max A few comments are in order. population December 24, 2020 1:39 manuscript_v2 t = t i1 t = t i2 t = t p1 t = t p2 HS 1 HS 2 100 150 200 250 300 350 400 450 500 Fig. 3. Spatial distribution of the infected population. The y-integrated population with respect to the x-position. The hotspots are marked with vertical lines. The parameters are f = f = 50, 1 2 = 0:2, = 0:15, and u = 0:28, d = 50, w = 10, D = 0:05. The snapshots are taken at the start-up time of both hotspots and the peak time of the second one. The development of a spatial interference between the two hotspots is clearly visible. The HS2 peak for the homogeneous case is 0:99 10 = 9:9, against an observed one of about 12, showing a 20 percent increase due to spatial interference from HS2. First, we see that the peak intensity is a fast decreasing function of the wind speed for all HS1 infection ratios well below the HS2 values. This is expected, since the infected in HS2 get replaced by less infected from HS1. However, upon increasing f in the vicinity and then above f , a shoulder appears 1 2 at intermediate wind speeds, indicating that a highly contagious mobile population from HS1 is capable of spoiling the bene cial eect of the wind. This is also a plausible result, since the infected removed by the wind in HS2 are quickly replaced by even more infected transmitted by HS1. This is the typical scenario dreaded by southern Italy towards the "stampede" from northern regions in the early stage of the Italian epidemics. Our simple model shows that such fears were indeed justi ed, an infectivity ratio f = 2 is already capable of producing a secondary peak in the curve and raising f only makes the situation worst, with the emergence of a whole range of wind speeds in which the peak intensity grows instead of decaying, almost reaching up to the dyN(x, y) December 24, 2020 1:39 manuscript_v2 f = 0.02 f = 0.1 0.8 f = 0.2 f = 0.4 f = 1.0 f = 1.4 0.6 f = 2.0 f = 3.0 f = 4.0 0.4 f = 7.0 f = 10.0 0.2 0.0 0.0 0.5 1.0 1.5 2.0 2.5 Fig. 4. Peak intensity in HS1 as a function of the wind speed at varying the infection rate in HS2 The simulation parameters are the same as in gure 3. We clearly observe the emergence of a non-monotonic wind speed regime in the range 1 < f < 5, followed by a loss of any bene cial wind eect above f 5. value of the windless case. Since this strongly reminds of the unstable region of a non-ideal equation of state, in which pressure goes down upon increasing density (condensation), we dub this eect "epidemic condensation". This is the main result of this paper, as it highlights the existence of an optimal wind speed u 0:5 which minimises the HS2 peak, and a second, higher, char- min acteristic speed u , beyond which the bene cial eects of the wind are restored. max By further increasing the relative infectivity of HS1, between ve and ten, no decay of the peak intensity at increasing wind speed above u is observed anymore in min the simulated window of the wind speed u, indicating that the presence of HS1 completely cancels any bene t of the wind speed above the optimal value u . min However, for a very large wind speed u ! 1 we expect again a decrease of the infected ratio, since for in nite wind speed the hot-spots become transparent again, and infectivity will have no impact. In Fig. 5 we report the duration of the epidemics as a function of u and f . A major peak is observed at low-wind, corresponding to the fact that the infected are convected away at very low rates. As expected, high infectivity goes with high peaks i /N max December 24, 2020 1:39 manuscript_v2 f= 0.02 f= 0.1 f= 0.2 f= 0.4 f= 1.0 f= 1.4 f= 2.0 f= 3.0 f= 4.0 f= 7.0 f= 10.0 2000 1000 0.00 0.02 0.04 0.06 0.08 0.10 0.12 0.0 0.5 1.0 1.5 2.0 2.5 Fig. 5. Duration of the epidemic in HS2. The curve shows a peak at very low wind speeds u 0:01, followed by a sequence of secondary peaks at higher speeds, all well below u . By and min large, wind speeds above u = 0:1 are consistently bene cial in shortening the epidemic duration. and short durations, the dreaded scenario for intensive care departments. As the wind speed increases, the local infected are eciently removed and the epidemic duration shortens. However, starting from comparatively low infectivity ratios, f = 0:2, further satellite peaks appear, indicating the existence of a sequence of wind speeds such that the duration grows back, if only mildly. This is again interpreted as a spatial interference eect, although we must caution that such measurement is very sensitive to small changes of the duration threshold, hence should be taken with great caution. 4. Qualitative scenario and discussion The ADSIR model presented in this paper focusses on the eects of spatial coupling, advection and diusion, on epidemic growth as dictated by local infection rates. It is well known that in the presence of random heterogeneities, such coupling can lead to highly nontrivial behaviour, such as the formation of striated infection highways. Here we take a simpler model problem, namely the eect of a primary hotspot (HS1) on the epidemic growth on a secondary hotspot (HS2) downstream HS1. In particular, we focus on the eect of a uniform "wind" at speed u, mimicking a t December 24, 2020 1:39 manuscript_v2 uniform human mobility across the two hotspots. In the absence of any wind, u = 0, and discounting diusion, the two hotspots evolve independently based on their corresponding infection rates. As soon as the wind is switched on, a bene cial eect is expected for both HSs because the wind sweeps infected individuals away into the "country side", where the chance to infect is much lower and healing can proceed nearly undisturbed. This is certainly true as soon as the wind speed exceeds the infection speed, namely the size of the hotspot divided by the typical infection timescale (reference speed), because, under such conditions, the wind blows susceptible individuals away before they have time to get signi cantly infected. So, the baseline expectation is that "wind is good", as it gives no time for infection to develop substantially. This is true for HS1, but not necessarily for HS2, which is exposed to the incoming ux of infected individuals from HS1. The quantitative question is whether, from the HS2 perspective, there exists an optimal wind speed which corresponds to a local minimum of the infection peak . In the following, we shall present evidence that the answer is in the armative. In particular, it is shown that as soon as HS1 is more infectious than HS2, the peak intensity in HS2 develops a much slower decay with the wind speed, and when HS1 is signi cantly more infectious than HS2, the HS2 peak increases at increasing wind speed, before it starts to decay again in the strong wind regime. In other words, the HS2 peak develops a non-monotonic dependence on the wind speed, with a local minimum, u at about half the reference speed and a local maximum u about min max twice as large. Such non-monotonic dependence bears an intriguing resemblance to a non-ideal equation of state, with the unstable branch in the wind speed region u u min u . Because of this close resemblance to equation of state of non-ideal gas, and max most notably to the unstable region where a density increase leads to a pressure decrease (condensation), we dub this eect epidemic condensation. We also monitor the duration of the epidemics as a function of the wind speed and infection rates. Note that while the peak intensity is the prime concern for health capacity issues, the duration bears directly on the mid-long term policies towards social and economic impact (many countries insisted on "curve attening" policies). Again, we nd that wind increase above a very low threshold is generally ben- e cial, although at increasing HS2 infectivity, the duration increases and shows repeated small-amplitude "sawtooth" oscillations. Such oscillations are yet another signature of spatial coupling, although their speci c nature remains to be fully ascertained. 5. Eect of the hotspot distance and the diusivity We also inspected the eect of the hotspot distance on epidemic condensation. To this purpose, we ran a series of simulations at dierent wind speeds and distances December 24, 2020 1:39 manuscript_v2 local minimum 0.80 local maximum power law fit 0.75 0.70 0.65 0.60 0.55 0.50 0.45 60 80 100 120 140 160 180 200 Fig. 6. The local minimum and maximum of the condensation curve. The infectivity ratio is set to be f = 4, all other parameters are chosen as in Fig. 4. Shown are the local minimum and maximum of the curve of the infected ratio in HS 2. The power-law t is performed with a small 0:32 exponent, d , hence the correlation eects decay slower than the distance. in the range 50 d 200, keeping a xed value f = 4. As expected, the local maximum observed in the condensation decreases with the distance and, less expectedly, so does the local minimum. Fig. 6 shows that both quantities decay according to an inverse power law d , with 1=3, indicating that the correlation between the two hotspots decays much more slowly than their inverse distance. To assess the eect of the diusivity, we computed the condensation curve { similar to Fig. 4 with xed f = 4, for dierent values of the diusion constants D. In Fig. 7 we observe a quantitative eect of the diusion parameter on the condensation curve, due to the fact that diusion smears out sharp spatial changes in population number, such as those observed at the hot-spot boundaries. Hence, the eect of increasing the diusivity D is similar to lowering the infec- tivity ratio f , as long as the diusivity remains suciently small enough, meaning by this that the Fisher speed U = Df remains well below the reference wind f 1 0 speed. In the simulations carried out here, U is always signi cantly smaller than U , hence these eects do not play any role. Whenever the condition U U is r f r violated, non trivial interference eects are expected, which may eventually lead to a revival of infectivity in the second hotspot. A detailed analysis of these eects i /N max December 24, 2020 1:39 manuscript_v2 0.9 0.8 0.7 0.6 0.5 D = 1 km /day D = 5 km /day 0.4 2 D = 10 km /day 0.0 0.5 1.0 1.5 2.0 2.5 Fig. 7. Eect of diusion. We change the diusivity and inspect its eect on the condensation curve for an intermediate high infectivity ratio f = 4. For small D, we observe no qualitative change of the curve. D = 5 km =day is the value we used throughout the rest of this paper. When the Fischer speed U = Df approaches the reference wind speed U , we expect additional 1 0 r interference eects which may lead a revival of infectivity in HS2. A detailed study of these phenomena is left to a future study. warrants a separate study on its own, hence it is deferred to future investigations. 6. Conclusions Summarizing, we have evidenced a non-monotonic relation between the wind speed and the peak intensity on the downstream hotspot as a function of the infectivity ratio with respect to the upstream one. Despite its drastic simpli cation of the mechanism of human mobility, it is hoped that the non-monotonic "constitutive relation" revealed by the present ADSIR model, as shown in Figures 4 and 5, may oer useful qualitative clues on the eects of spatial interference between infected hotspots. Acknowledgements SS kindly acknowledges funding from the European Research Council under the European Union's Horizon 2020 Framework Programme (No. FP/2014-2020)/ERC i /N max December 24, 2020 1:39 manuscript_v2 Grant Agreement No. 739964 (COPMAT). JD was supported by the ERASMUS program and Physics Advanced program of the Elite Network Bavaria (University of Regensburg, Germany). One of the authors would like to acknowledge useful discussions with Prof E. Marinari and G. Parisi in an early stage of the project. References 1. J. H. University, Covid-19 dashboard worldwide 2020. 2. B. Armocida, B. Formenti, S. Ussai, F. Palestra and E. Missoni, The Lancet Public Health 5, p. e253 (2020). 3. H. Legido-Quigley, J. T. Mateos-Garc a, V. R. Campos, M. Gea-S anchez, C. Muntaner and M. McKee, The Lancet Public Health 5, e251 (2020). 4. Q. Chen, M. Liang, Y. Li, J. Guo et al., The Lancet Psychiatry 7, e15 (2020). 5. L. Duan and G. Zhu, The Lancet Psychiatry 7, 300 (2020). 6. M. Nicola, Z. Alsa , C. Sohrabi, A. Kerwan et al., International Journal of Surgery 78, 185 (2020). 7. A. A. Toda, (2020) arXiv:2003.11221. 8. A. L. Lloyd and R. M. May, Journal of Theoretical Biology 179, 1 (1996). 9. K. Dietz, Journal of Mathematical Biology 8, 291 (1979). 10. F. Ball, D. Sirl and P. Trapman, Mathematical Biosciences 224, 53 (2010). 11. J. C. Lang, H. De Sterck, J. L. Kaiser and J. C. Miller, Journal of Complex Networks 6, 948 (2018). 12. M. U. Kraemer, C. H. Yang, B. Gutierrez, C. H. Wu et al., Science 368, 493 (2020). 13. M. Chinazzi, J. T. Davis, M. Ajelli, C. Gioannini et al., Science 368, 395 (2020). 14. W. O. Kermack and A. G. McKendrick, Proceedings of the Royal Society A: Mathe- matical, Physical and Engineering Sciences 115, 700 (1927). 15. S. Kaushal, A. S. Rajput, S. Bhattacharya, M. Vidyasagar et al., (2020) arXiv:2006.00045 . 16. E. Kaxiras, G. Neofotistos and E. Angelaki, Chaos, Solitons and Fractals 138, 1 (2020). 17. S. Pathak, A. Maiti and G. P. Samanta, Nonlinear Analysis: Modelling and Control 15, 71 (2010). 18. C. Ji, D. Jiang and N. Shi, Physica A: Statistical Mechanics and its Applications 390, 1747 (2011). 19. W. C. Roda, M. B. Varughese, D. Han and M. Y. Li, Infectious Disease Modelling 5, 271 (2020). 20. G. C. Cala ore, C. Novara and C. Possieri, (2020) arxiv:2003.14391. 21. T. Chotibut, D. R. Nelson and S. Succi, Physica A: Statistical Mechanics and its Applications 465, 500 (2017).

Physics – arXiv (Cornell University)

**Published: ** Aug 28, 2020

Loading...

You can share this free article with as many people as you like with the url below! We hope you enjoy this feature!

Read and print from thousands of top scholarly journals.

System error. Please try again!

Already have an account? Log in

Bookmark this article. You can see your Bookmarks on your DeepDyve Library.

To save an article, **log in** first, or **sign up** for a DeepDyve account if you don’t already have one.

Copy and paste the desired citation format or use the link below to download a file formatted for EndNote

Access the full text.

Sign up today, get DeepDyve free for 14 days.

All DeepDyve websites use cookies to improve your online experience. They were placed on your computer when you launched this website. You can change your cookie settings through your browser.