EPI’s Covid-19 forecast model: Sars-CoV-2 spread in Florida

This model was produced in partnership between the authors listed below and the Center for Inference and Dynamics of Infectious Diseases (CIDID). For additional COVID-19 models by this team which cover the contiguous 48 U.S. states, please visit the Global Epidemic and Mobility Modeling (GLEAM) project site.

Authors

Kunpeng Mu {1}, Ana Pastore y Piontti {2}, Matteo Chinazzi {1}, Jessica T. Davis {1}, Xinyue Xiong {1}, and Alessandro Vespignani {1}, Ira M Longini {2}

  1. Laboratory for the Modeling of Biological and Socio-technical Systems, Northeastern University, Boston, MA USA
  2. Department of Biostatistics, College of Public Health and Health Professions, and Emerging Pathogens Institute, University of Florida, Gainesville, USA

Introduction

Motivated by the local spread of SARS-CoV-2 in the United States we study the impact of different non-pharmaceutical interventions to mitigate the outbreak by using modeling scenarios. In particular, we focus on the effects of social distancing  measures (school closure, “stay home” or “shelter in place” policies) that may mitigate or suppress the epidemic. Modeling suggests that “Stay Home” mitigation policies are effective if sustained for several weeks. We provide a specific model that takes into account the mitigation measures as implemented in the State of Florida and project infections and deaths in the next four weeks.

Disclaimer: There are very large uncertainties around the transmission of SARS-CoV-2, the  effectiveness of different policies and the extent to which the population is compliant to social distancing policies. The presented material is based on modeling scenario assumptions  informed by current knowledge of the disease, and subject to change as more data will be available. Future decisions on when and for how long to relax policies will need to be informed by ongoing surveillance. Additional modeling and data studies are required to assess the level and effectiveness of additional non-pharmaceuticals interventions required to lift current social distancing interventions.

Results

Key projections Results

  • Our model for Florida points to the week of April 5,2020 as the peak day for hospital and ICU beds needed.
  • Based on the last projections, a total of 516 COVID-19 deaths (range 317-745) are currently predicted through April 26th, 2020.

Figure 1A shows the daily incidence of deaths for an unmitigated scenario (no-interventions) and the scenario implementing social distancing. On average the peak for the deaths daily incidence is in the first two weeks of April, 2020. In Table 1 we report  the mean expected cumulative number of deaths and the 90%CI, in the case of a “stay home” policy extended during the month of April, 2020.

Figure 1B shows the daily incidence of new infections in Florida. This number is indicative of the epidemic trajectory in the state. However, it is not directly comparable with the number of confirmed cases because of the large fraction of subclinical cases and the variations in testing rate. The plot shows however that the number of infections is expected to decline starting in the second week of April. In fig. 1B) we report the cumulative attack rate (infections per 10,000 individuals in the population). In Table 1, we report the cumulative attack rate (percentage of infected population) during the month of April 2020. These estimates should be revisited as more data on clinical and subclinical infections by age will be available (see MedRxiv article, Age-dependent effects in the transmission and control of COVID-19 epidemics).

It is important to stress that the curves for new infections and deaths are not accounting for delay in reporting, so inflection points in real data may be delayed even more with respect to the start of interventions. Real data curves are also expected to be more noisy because of the clumping of cases or deaths from different dates that are reported on a single day.

In Fig. 1B we report the number of hospitalized patients with COVID-19. Figure 2B shows the number of ICU beds needed under the different assumptions. In an unmitigated scenario the state capacity could have been surpassed by the end of March. In the case of considering interventions, the number of ICU beds appear within capacity. However, possible heterogeneities across the state and potential clusters of cases in senior facilities might generate spikes in demand that the model cannot account for. The results on hospital and ICU occupancy is based on the assumption of a 30% overall rate of sub-clinical cases among infected people (see MedRxiv article, Age-dependent effects in the transmission and control of COVID-19 epidemics). Sensitivity analysis for different rates ranging from 20% to 50% is available on request.

Discussion

In summary the results show that aggressive social distancing  measures (“stay home” or “shelter in place” policies) may keep the epidemic under control and induce a decay of the epidemic activity only if sustained for several weeks. These results suggest that additional non-pharmaceuticals interventions at scale such as testing, contact tracing, early isolation of exposed contacts and asymptomatic individuals must be implemented before lifting stringent social distancing interventions. Additional modeling and data studies are required to assess the level and effectiveness of additional non-pharmaceuticals interventions required to lift current social distancing interventions.

A figure that shows two line graphs for the two scenarios: one for number of hospitalized individuals and one for number of ICU beds.
Figure 2. A) Number of hospitalized individuals and B) Number of ICU beds needed for the two scenarios we consider here: an “Unmitigated” scenario and a “Mitigated” scenario. The solid lines represent the mean values, while the shaded areas the 90% confidence interval. The dark blue line is the number of ICUs occupied by COVID cases by using a 70% rate clinical case over all infections (1.4 the model estimate).
A figure that shows two line graphs for the two scenarios: one for number of hospitalized individuals and one for number of ICU beds.
Figure 2. A) Number of hospitalized individuals and B) Number of ICU beds needed for the two scenarios we consider here: an “Unmitigated” scenario and a “Mitigated” scenario. The solid lines represent the mean values, while the shaded areas the 90% confidence interval. The dark blue line is the number of ICUs occupied by COVID cases by using a 70% rate clinical case over all infections (1.4 the model estimate).
This table showcases the cumulative attack rate (AR) and cumulative number of deaths in the unmitigated scenario versus the scenario with a stay at home order.
Table 1: Cumulative attack rate (AR) and cumulative number of deaths in the unmitigated scenario versus the scenario where we consider a stay at home order from April 3 until the end of April. The deaths estimate assumes the infectious fatality rate reported in [5].

Projection Model

We developed an age-structured compartmental model of SARS-CoV-2 transmission. This model accounts specifically for the different classes of infectious individuals (symptomatic, pre-symptomatic, and asymptomatic), as well as patients who are hospitalized and in need of intensive care units (ICU). Age mixing patterns specific for each state in the United States are derived from Mistry et al. [1]. To initialize the model in each state we consider an initialization provided by the Global Epidemic and Mobility model (GLEAM) used to study the worldwide spread of the COVID-19 outbreak from mainland China [2]. The model is calibrated using international exportations from China. The parameters defining some of the major characteristics of the model are reported in the Appendix 1.

It is important to stress that the following analysis is based on models that draw on previous research. The assumptions and caveats are discussed in a separate section. Furthermore, the model may change as new data become available, thus providing specific characteristic rates and times estimates.

We streamlined the interventions according to the following timeline:

  • Before March 17: Business as usual
  • March18: Schools are closed.
  • March 20: Social distancing is enforced with the closing of non essential businesses.
  • April 1st  – Governor issues a statewide “Stay Home” order to be effective on April 3 for 30 days, limiting all activity in Florida to essential services.

This is translated in three different modeling implementations of the social contacts in the population:

  • No intervention (NI): considers schools are open until the summer following a usual schedule and re-open in September. In addition, we consider that schools are closed over the weekends. In this way the NI scenario, even though it does not consider interventions, accounts for the reduction of the number of contacts due to weekend school closure that would normally happen.
  • School closure: schools are closed.
  • SC-WR-SD: school closure, work reduction, social distancing (regular interaction in households, no school, 50% reduced work interaction, 10% casual social contacts). This corresponds to a transmissibility reduction of about 60%.
  • “Stay Home” scenario: we assume a stay at home policy (in some states referred to as “shelter in place”) from March 16 through March 31st and then from April 16 onward for different periods of time, namely 2, 4, 8 and 12 weeks. This corresponds to a transmissibility reduction of about 70%.

We initialize our model according to the GLEAM outputs. Until March 17th, we consider there is no/minimal intervention in place. On March 18th according to the above timeline we start interventions ramping up on April 3rd to the “stay home” scenario  that stays until the end of April.

Note: we do not consider in the scenarios specific strategies for enhancing contact tracing, testing and early isolation of cases. We are also not considering potential changes to the virus transmissibility due to environmental factors, in particular seasonal drivers such as temperature and humidity.

The model is stochastic and age structured, considering the interaction of individuals in single years from 0 to 85+. The contact patterns account for the interaction of individuals in different settings: households, schools, workplaces, and the general community. This compartmentalization of settings allows us to implement the different interventions, controlling for the level of interaction in each one of them separately. In particular we can provide school closure, reduction of work-place contacts and causal community contacts.

Caveats and assumption

These results were obtained under several assumptions. The first is that we use modeling estimates for the effect of school closures, smart working and social distancing effects on the transmissibility of SARS-CoV-2.  The model is considering the ICU capacity at the state level, however relevant heterogeneity may lead to exceed this capacity in specific urban areas earlier than at the state level. All estimates do not consider the likely introduction of specific suppression or mitigation policies issued to lower the transmissibility in specific states that experience elevated epidemic activity. Finally, we are not considering superspreading events and differential transmissibility across age brackets. Even in the presence of these limitations we hope that the information provided here can be of help in assessing the risk of transmission of SARS-CoV-2 in the US.

References

  1. Mistry, D., Litvinova, M., Chinazzi, M., Fumanelli, L., Gomes, M.F., Haque, S.A., Liu, Q.H., Mu, K., Xiong, X., Halloran, M.E. and Longini Jr, I.M., Merler, S., Ajelli, M, and Vespignani, A., 2020. Inferring high-resolution human mixing patterns for disease modeling. arXiv preprint arXiv:2003.01214.
  2. Chinazzi, Matteo, Jessica T. Davis, Marco Ajelli, Corrado Gioannini, Maria Litvinova, Stefano Merler, Ana Pastore y Piontti et al. “The effect of travel restrictions on the spread of the 2019 novel coronavirus (COVID-19) outbreak.” Science (2020).
  3. The White House, 15 days to slow the spread. Available online: https://www.whitehouse.gov/wp-content/uploads/2020/03/03.16.20_coronavirus-guidance_8.5x11_315PM.pdf
  4. Data were obtained from the American Hospital Directory and are based on public records obtained from the US Centers for Medicare and Medicaid Services. For additional information see https://www.ahd.com/data_sources.html
  5. Verity, R. et al. “Estimates of the severity of coronavirus disease 2019: a model-based analysis”, Lancet Infectious diseases, https://doi.org/10.1016/S1473-3099(20)30243-7

Appendix 1

ParameterValue
R-naught (R0) 2.2-2.8
Incubation period 4 days
Infection fatality rate (%)0-4 years old: 0.1%
5-17 years old: 0.2166%
18-49 years old: 1.475%
50-64 years old: 12.33%
65+ years old: 22.0813%
Hospitalization ration (%)0-4 years old: 0.002%
5-17 years old: 0.00433%
18-49 years old: 0.0364%
50-64 years old: 1.1333%
65+ years old: 4.45625%
ICU among hospitalized0-4 years old: 5%
5-17 years old: 5%
18-49 years old: 5%
50-64 years old: 17.267%
65+ years old: 39.99%

Written by: DeLene Beeland