Stochastic extinction of epidemics: how long does would it take for SARS-Cov-2 to die out without herd immunity?

Abstract

Worldwide, we are currently in an unprecedented situation with regard to the SARS-Cov-2 epidemic, where countries are using isolation and lock-down measures to control the spread of infection. This is a scenario generally not much anticipated by previous theory, and in particular, there has been little attention paid to the question of extinction as a means to eradicate the virus; the prevailing view appears to be that this is unfeasible without a vaccine. We use a simple well-mixed stochastic SIR model as a basis for our considerations, and calculate a new result, using branching process theory, for the distribution of times to extinction. Surprisingly, the distribution is an extreme value distribution of the Gumbel type, and we show that the key parameter determining its mean and standard deviation is the expected rate of decline ρ_e = γ(1 −R_e) of infections, where γ is the rate of recovery from infection and R_e is the usual effective reproductive number. The result also reveals a critical threshold number of infected I^† =1/(1 − R_e), below which stochastic forces dominate and need be considered for accurate predictions. As this theory ignores migration between populations, we compare against a realistic spatial epidemic simulator and simple stochastic simulations of sub-divided populations with global migration, to find very comparable results to our simple predictions; in particular, we find global migration has the effect of a simple upwards rescaling of R_e with the same Gumbel extinction time distribution we derive from our non-spatial model. Within the UK, using recent estimates of ≈ 37, 000 infected and R_e ≈ 0.9, this model predicts a mean extinction time of 616 ± 90 days or approximately ∼ 2 years, but could be as short as 123 ± 15 days, or roughly 4 months for R_e =0.4. Globally, the theory predicts extinction in less than 200 days, if the reproductive number is restricted to R_e < 0.5. Overall, these results highlight the extreme sensitivity of extinction times when R_e approaches 1 and the necessity of reducing the effective reproductive number significantly (R_e ≲ 0.5) for relatively rapid extinction of an epidemic or pandemic.