An international assessment of the COVID-19 pandemic using ensemble data assimilation

Evensen, G.; Amezcua, J.; Bocquet, M.; Carrassi, A.; Farchi, A.; Fowler, A.; Houtekamer, P.; Jones, C. K. R. T.; de Moraes, R.; Pulido, M.; Sampson, C.; Vossepoel, F.

This work shows how one can use iterative ensemble smoothers to effectively estimate parameters of an SEIR model with age-classes and compartments of sick, hospitalized, and dead. The data conditioned on are the daily numbers of accumulated deaths and the number of hospitalized. Also, it is possible to condition on the number of cases obtained from testing. We start from a wide prior distribution for the model parameters; then, the ensemble conditioning leads to a posterior ensemble of estimated parameters leading to model predictions in close agreement with the observations. The updated ensemble of model simulations have predictive capabilities and include uncertainty estimates. In particular, we estimate the effective reproductive number as a function of time, and we can assess the impact of different intervention measures. By starting from the updated set of model parameters, we can make accurate short-term predictions of the epidemic development given knowledge of the future effective reproductive number. Also, the model system allows for the computation of long-term scenarios of the epidemic under different assumptions. We have applied the model system on data sets from several countries with vastly different developments of the epidemic, and we can accurately model the development of the COVID-19 outbreak in these countries. We realize that more complex models, e.g., with regional compartments, may be desirable, and we suggest that the approach used here should be applicable also for these models.