The turning point and end of an expanding epidemic cannot be precisely forecast

Significance Susceptible–infected–removed (SIR) models and their extensions are widely used to describe the dynamics of infection spreading. Certain generic features of epidemics are well-illustrated by these models, which can be remarkably good at reproducing empirical data through suitably chosen parameters. However, this does not assure a good job anticipating the forthcoming stages of the process. To illustrate this point, we accurately describe the propagation of COVID-19 in Spain using one such model and show that predictions for its subsequent evolution are disparate, even contradictory. The future of ongoing epidemics is so sensitive to parameter values that predictions are only meaningful within a narrow time window and in probabilistic terms, much as what we are used to in weather forecasts. Epidemic spread is characterized by exponentially growing dynamics, which are intrinsically unpredictable. The time at which the growth in the number of infected individuals halts and starts decreasing cannot be calculated with certainty before the turning point is actually attained; neither can the end of the epidemic after the turning point. A susceptible–infected–removed (SIR) model with confinement (SCIR) illustrates how lockdown measures inhibit infection spread only above a threshold that we calculate. The existence of that threshold has major effects in predictability: A Bayesian fit to the COVID-19 pandemic in Spain shows that a slowdown in the number of newly infected individuals during the expansion phase allows one to infer neither the precise position of the maximum nor whether the measures taken will bring the propagation to the inhibition regime. There is a short horizon for reliable prediction, followed by a dispersion of the possible trajectories that grows extremely fast. The impossibility to predict in the midterm is not due to wrong or incomplete data, since it persists in error-free, synthetically produced datasets and does not necessarily improve by using larger datasets. Our study warns against precise forecasts of the evolution of epidemics based on mean-field, effective, or phenomenological models and supports that only probabilities of different outcomes can be confidently given.

[1]  M. R. Ferrández,et al.  Mathematical modeling of the spread of the coronavirus disease 2019 (COVID-19) taking into account the undetected infections. The case of China , 2020, Communications in Nonlinear Science and Numerical Simulation.

[2]  Herbert W. Hethcote,et al.  The Mathematics of Infectious Diseases , 2000, SIAM Rev..

[3]  T. de-Camino-Beck A modified SEIR Model with Confinement and Lockdown of COVID-19 for Costa Rica , 2020, medRxiv.

[4]  D. Sornette,et al.  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 , 2020, Nonlinear dynamics.

[6]  W. O. Kermack,et al.  A contribution to the mathematical theory of epidemics , 1927 .

[7]  David L. Smith,et al.  Mapping global variation in human mobility , 2020, Nature Human Behaviour.

[8]  H. Hethcote,et al.  Effects of quarantine in six endemic models for infectious diseases. , 2002, Mathematical biosciences.

[9]  Michael Y. Li,et al.  Why is it difficult to accurately predict the COVID-19 epidemic? , 2020, Infectious Disease Modelling.

[10]  H. Thieme,et al.  Recurrent outbreaks of childhood diseases revisited: the impact of isolation. , 1995, Mathematical biosciences.

[11]  Ruiyun Li,et al.  Substantial undocumented infection facilitates the rapid dissemination of novel coronavirus (SARS-CoV-2) , 2020, Science.

[12]  C. Beauchemin,et al.  Duration of SHIV production by infected cells is not exponentially distributed: Implications for estimates of infection parameters and antiviral efficacy , 2017, Scientific Reports.

[13]  Liangrong Peng,et al.  Epidemic analysis of COVID-19 in China by dynamical modeling , 2020, medRxiv.

[14]  Sergio Correia,et al.  Pandemics Depress the Economy, Public Health Interventions Do Not: Evidence from the 1918 Flu , 2020, The Journal of Economic History.

[15]  S. Bhatt,et al.  Report 13: Estimating the number of infections and the impact of non-pharmaceutical interventions on COVID-19 in 11 European countries , 2020 .

[16]  A. Flahault,et al.  Seroprevalence of anti-SARS-CoV-2 IgG antibodies in Geneva, Switzerland (SEROCoV-POP): a population-based study , 2020, The Lancet.

[17]  D. Brockmann,et al.  Effective containment explains subexponential growth in recent confirmed COVID-19 cases in China , 2020, Science.

[18]  Johannes Zierenberg,et al.  Inferring change points in the spread of COVID-19 reveals the effectiveness of interventions , 2020, Science.

[19]  David Soriano-Panos,et al.  Derivation of the effective reproduction number R for COVID-19 in relation to mobility restrictions and confinement , 2020, medRxiv.

[20]  J. Gómez-Gardeñes,et al.  A mathematical model for the spatiotemporal epidemic spreading of COVID19 , 2020, medRxiv.

[21]  X. Rodó,et al.  A modified SEIR model to predict the COVID-19 outbreak in Spain and Italy: Simulating control scenarios and multi-scale epidemics , 2020, Results in Physics.

[22]  J. Banga,et al.  Structural identifiability and observability of compartmental models of the COVID-19 pandemic , 2020, Annual Reviews in Control.

[23]  Martyn Plummer,et al.  JAGS: A program for analysis of Bayesian graphical models using Gibbs sampling , 2003 .

[24]  D. Sornette,et al.  Generalized logistic growth modeling of the COVID-19 outbreak in 29 provinces in China and in the rest of the world , 2020, medRxiv.

[25]  Timoteo Carletti,et al.  COVID-19: The unreasonable effectiveness of simple models , 2020, Chaos, Solitons & Fractals: X.

[26]  Ernesto Estrada COVID-19 and SARS-CoV-2. Modeling the present, looking at the future , 2020, Physics Reports.

[27]  M. Levitt,et al.  Predicting the Trajectory of Any COVID19 Epidemic From the Best Straight Line , 2020, medRxiv.

[28]  Carl A. B. Pearson,et al.  The effect of control strategies to reduce social mixing on outcomes of the COVID-19 epidemic in Wuhan, China: a modelling study , 2020, The Lancet Public Health.

[29]  A. Rinaldo,et al.  Spread and dynamics of the COVID-19 epidemic in Italy: Effects of emergency containment measures , 2020, Proceedings of the National Academy of Sciences.

[30]  E. Lorenz Deterministic nonperiodic flow , 1963 .

[31]  S. Leeder,et al.  A population based study , 1993, The Medical journal of Australia.

[32]  E. Lorenz,et al.  Predictability: Does the Flap of a Butterfly’s Wings in Brazil Set off a Tornado in Texas? , 2013 .

[33]  Ruiyun Li,et al.  Substantial undocumented infection facilitates the rapid dissemination of novel coronavirus (COVID-19) , 2020, medRxiv.

[34]  Tangchun Wu,et al.  Reconstruction of the full transmission dynamics of COVID-19 in Wuhan , 2020, Nature.

[35]  Sergei Maslov,et al.  Modeling COVID-19 dynamics in Illinois under non-pharmaceutical interventions , 2020, medRxiv.

[36]  M. Hernán,et al.  Prevalence of SARS-CoV-2 in Spain (ENE-COVID): a nationwide, population-based seroepidemiological study , 2020, The Lancet.

[37]  Marta Sales-Pardo,et al.  A Bayesian machine scientist to aid in the solution of challenging scientific problems , 2020, Science Advances.