Network-based simulations of re-emergence and spread of COVID-19 in Aotearoa New Zealand

We simulate the late July/early August re-emergence and spread of COVID-19 in Aotearoa New Zealand. We use a stochastic, individual-based network model of all ≈5 million individuals in Aotearoa, and run simulations for a period of 30 days. Based on these simulations, we calculate: the expected time to detection of the first case after initial seed cases; the number of cases at the time of detection; the time until detection of a first case outside of Auckland; and how the overall number of cases increases without intervention. Our model includes interaction pathways, referred to as ‘contexts’ in the network, broken down into network ‘layers’ representing home, work, school, and community structure. Each simulation starts from initial (seed) cases corresponding to the first detected re-emergence cases in August 2020. We run 50 realisations of each simulation for 30 days — each simulation scenario corresponding to one of three different levels of transmission rate. To model the behaviour of individuals in the weeks prior to the August 11th re-emergence, we assume a moderate rate of people getting tested if mildly symptomatic. No contact tracing or intervention is present in this scenario, other than cases that test positive being isolated to their dwelling.

[1]  D. Kendall An Artificial Realization of a Simple “Birth-And-Death” Process , 1950 .

[2]  D. Gillespie A General Method for Numerically Simulating the Stochastic Time Evolution of Coupled Chemical Reactions , 1976 .

[3]  M. Lipsitch,et al.  How generation intervals shape the relationship between growth rates and reproductive numbers , 2007, Proceedings of the Royal Society B: Biological Sciences.

[4]  S. Blower,et al.  Theory versus Data: How to Calculate R0? , 2007, PloS one.

[5]  R. Mikolajczyk,et al.  Social Contacts and Mixing Patterns Relevant to the Spread of Infectious Diseases , 2008, PLoS medicine.

[6]  Lucas C. Parra,et al.  Origins of power-law degree distribution in the heterogeneity of human activity in social networks , 2013, Scientific Reports.

[7]  P. Boxall,et al.  International Social Survey Programme (ISSP) for New Zealand , 2015 .

[8]  Joel C. Miller,et al.  Mathematics of Epidemics on Networks: From Exact to Approximate Models , 2017 .

[9]  Dennis Andersson,et al.  A retrospective cohort study , 2018 .

[10]  Joel C. Miller,et al.  EoN (Epidemics on Networks): a fast, flexible Python package for simulation, analytic approximation, and analysis of epidemics on networks , 2019, J. Open Source Softw..

[11]  M. Plank,et al.  Suppression and Mitigation Strategies for Control of COVID-19 in New Zealand , 2020, medRxiv.

[12]  Michael J. Plank,et al.  A stochastic model for COVID-19 spread and the effects of Alert Level 4 in Aotearoa New Zealand , 2020, medRxiv.

[13]  L. Kucirka,et al.  Variation in False-Negative Rate of Reverse Transcriptase Polymerase Chain Reaction–Based SARS-CoV-2 Tests by Time Since Exposure , 2020, Annals of Internal Medicine.

[14]  Yongsheng Wu,et al.  Epidemiology and transmission of COVID-19 in 391 cases and 1286 of their close contacts in Shenzhen, China: a retrospective cohort study , 2020, The Lancet Infectious Diseases.

[15]  Lucie Abeler-Dörner,et al.  Quantifying SARS-CoV-2 transmission suggests epidemic control with digital contact tracing , 2020, Science.

[16]  Demival Vasques Filho,et al.  Transitivity and degree assortativity explained: The bipartite structure of social networks , 2019, Physical review. E.

[17]  Mikhail Prokopenko,et al.  Modelling transmission and control of the COVID-19 pandemic in Australia , 2020, Nature Communications.

[18]  Junling Ma Estimating epidemic exponential growth rate and basic reproduction number , 2020, Infectious Disease Modelling.

[19]  N. Linton,et al.  Serial interval of novel coronavirus (COVID-19) infections , 2020, International Journal of Infectious Diseases.