Skip directly to site content Skip directly to page options Skip directly to A-Z link Skip directly to A-Z link Skip directly to A-Z link
Volume 27, Number 9—September 2021
Research

Estimating the Impact of Statewide Policies to Reduce Spread of Severe Acute Respiratory Syndrome Coronavirus 2 in Real Time, Colorado, USA

Andrea G. Buchwald1, Jude Bayham, Jimi Adams, David Bortz, Kathryn Colborn, Olivia Zarella, Meghan Buran, Jonathan Samet, Debashis Ghosh, Rachel Herlihy, and Elizabeth J. CarltonComments to Author 
Author affiliations: Colorado School of Public Health, Aurora, Colorado, USA (A.G. Buchwald, O. Zarella, M. Buran, J. Samet, D. Ghosh, E.J. Carlton); Colorado State University, Fort Collins, Colorado, USA (J. Bayham); University of Colorado, Denver, Colorado, USA (J. Adams); University of Colorado, Boulder, Colorado, USA (D. Bortz); University of Colorado School of Medicine, Aurora (K. Colborn); Colorado Department of Public Health and Environment, Denver (R. Herlihy)

Cite This Article

Abstract

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic necessitated rapid local public health response, but studies examining the impact of social distancing policies on SARS-CoV-2 transmission have struggled to capture regional-level dynamics. We developed a susceptible-exposed-infected-recovered transmission model, parameterized to Colorado, USA‒specific data, to estimate the impact of coronavirus disease‒related policy measures on mobility and SARS-CoV-2 transmission in real time. During March‒June 2020, we estimated unknown parameter values and generated scenario-based projections of future clinical care needs. Early coronavirus disease policy measures, including a stay-at-home order, were accompanied by substantial decreases in mobility and reduced the effective reproductive number well below 1. When some restrictions were eased in late April, mobility increased to near baseline levels, but transmission remained low (effective reproductive number <1) through early June. Over time, our model parameters were adjusted to more closely reflect reality in Colorado, leading to modest changes in estimates of intervention effects and more conservative long-term projections.

Mathematical transmission models are useful tools for predicting the magnitude, duration, and severity of the severe acute respiratory coronavirus 2 (SARS-CoV-2) pandemic. However, widely used national-level models might not capture regional heterogeneity. The coronavirus disease (COVID-19) outbreak in Colorado, USA, has been the subject of numerous discrepant projections from the Institute for Health Metrics and Evaluation and other modeling groups (1), which might have structural and data source explanations, highlighting the need for ensuring that models are fit to local epidemiologic data (24).

We report on our experience using a locally tailored model to inform policy in Colorado. Social distancing policies, intended to decrease contact rates, have been cornerstone public health tools for pandemic control, and these strategies have been adopted to control SARS-CoV-2 globally (2,5). Until recently, evidence of the effects of social distancing has come primarily from studies of the consequences of school and transit closures on influenza transmission (3,4,6). Early evidence suggests that social distancing policies can suppress transmission of SARS-CoV-2 (7,8), and recent evidence suggests a strong correlation between mobility and transmission reduction (9). However, these studies largely focused on periods when social distancing policies were in place, leaving critical questions unanswered regarding how long populations will comply with such measures and what happens when policies are relaxed.

One of the defining characteristics of the COVID-19 pandemic is the need for rapid response in the face of imperfect and incomplete information. Mathematical models of infectious disease transmission can be used in real-time to estimate parameters, such as the effective reproductive number (Re) and the efficacy of current and future intervention measures, providing time-sensitive data to policy-makers (10). We describe development of such a model, in close collaboration with the Colorado Department of Health and Environment and the Governor’s office, to gauge the current and future effects of early policies to decrease social contacts and, later, the gradual relaxing of stay-at-home orders.

We developed a compartmental susceptible-exposed-infected-recovered (SEIR) model calibrated to statewide COVID-19 case and hospitalization data to estimate changes in the contact rate and the Re after emergence of SARS-CoV-2 and the implementation of statewide social distancing policies in Colorado. We supplemented model estimates with an analysis of mobility by using mobile-device location data. Estimates were generated in near real time, at multiple time-points, with a rapidly evolving understanding of SARS-CoV-2. At each time point, we generated projections of the possible course of the outbreak under future social distancing scenarios. Findings were regularly provided to key Colorado decision-makers. We present estimates generated at multiple time points to document how our model, estimates and projections evolved over time. Although our analysis is specific to Colorado, our experience highlights the need for locally calibrated transmission models to inform public health preparedness and policymaking, along with ongoing analyses of the impact of policies to slow the spread of SARS-CoV-2.

Methods

COVID-19 Timeline and Policies

The first SARS-CoV-2 case was reported in Colorado on March 5, 2020, and a rapid succession of policies to control transmission followed (Table 1). The Colorado governor formally declared a disaster emergency on March 11. During March 14‒April 16, a total of 35 executive orders were mandated to curb SARS-CoV-2 transmission, including school closures, reduction in workforce percentages, and shelter-in-place (stay-at-home) orders. In conjunction with state executive orders, the Colorado Department of Health and Environment issued orders closing restaurants, bars, and other congregate environments on March 17 and prohibiting gatherings of >10 persons on March 19. A state-wide stay-at-home order was in effect during March 26‒April 26. Transition to a less restrictive phase, safer-at-home, began on April 27, which enabled some businesses to reopen with restrictions. The metropolitan Denver counties, comprising ≈50% of the population of Colorado, were under extended stay-at-home orders until May 8.

Reported Case and Hospitalization Data

Hospitalization data are a robust indicator of transmission trends compared with reported case data because reported case data are sensitive to testing capacity. However, because COVID-19 hospitalization data were sparse early in the epidemic, we initially fit models to reported COVID-19 cases from the Colorado Electronic Disease Reporting System (CEDRS). We fit models to the daily number of symptom onsets to reflect a biologically meaningful process (report date can be sensitive to testing lags). Missing onset dates were imputed as report date minus 7 days, the median onset-to-report lag. In May, we began fitting models to the daily number of hospitalized COVID-19 patients because we suspected that reported cases captured a variable proportion of infections over time because of increases in testing capacity. Daily hospital census records were obtained from EMResource (https://emresource.juvare.com). Because EMResource appeared to underreport COVID-19 hospitalizations during March, we inferred COVID-19 hospitalizations by using CEDRS before April 8.

Model Description

We used a deterministic age-structured SEIR model with 3 age groups (<30, 30–59, and >60 years of age) to estimate key model parameters and project the number of COVID-19 hospitalizations. In the model, we assume a single virus introduction event occurring on January 24, a date extrapolated from the first reported cases in Colorado.

In the model, the probability that an infected person shows development of symptoms (13) and needs hospitalization or ICU care is age dependent (14). All persons have an equal probability of exposure and infection, regardless of age. Initially, we used published estimates (15) for the proportion of symptomatic case-patients requiring hospitalization and critical care. Starting in May, with sufficient hospitalization data, we generated Colorado-specific estimates of the probability of hospitalization and critical care among case-patients by using model-fitting approaches, which enabled us to better account for underlying health status and patterns of care in Colorado (16).

Figure 1

Emergence of COVID-19, Colorado, USA, 2020, showing policy-based responses (A) and definition of 4 distinct social distancing phases (B) corresponding with early closures (phase 1, March 17‒25); statewide stay-at-home (phase 2, March 26‒April 26), statewide partial transition to safer-at-home (phase 3, April 27–May 8); statewide safer-at-home (phase 4, May 5–June 3). Social distancing parameters were estimated at 4 points during March‒June by using model fitting procedures and reported case data (fits 1 and 2) and hospital census data (fits 3 and 4). In light of the 5.1 day mean incubation period, the ≈7-day lag between symptom onset and case report, and the ≈8-days between symptom onset and hospitalization based on state records, there are ≈12 and 13 day lags between infection and case report, and infection and hospitalization, respectively (gray boxes). Thus, at each model fit, we could estimate social distancing parameters reflecting transmission conditions up to 12 (fits 1 and 2) or 13 (fits 3 and 4) days before the fit date. Asterisks (*) indicate estimate generated for only part of the period. COVID-19, coronavirus disease; P, phase.

Figure 1. Emergence of COVID-19, Colorado, USA, 2020, showing policy-based responses (A) and definition of 4 distinct social distancing phases (B) corresponding with early closures (phase 1, March 17‒25); statewide stay-at-home (phase...

The model includes 3 types of transmission-reducing parameters: social distancing, mask wearing, and self-isolation of symptomatic persons. Social distancing was modeled as a reduction in the contact rate between susceptible and infectious persons by multiplying the transmission parameter, β, by (1 ‒ social distancing). We defined contacts as interactions that could enable spread of infections from an infected person to a susceptible person. The term social distancing is used to encompass all contact-reducing behaviors and policies, including working from home, school closures, maintaining physical distancing, socializing outdoors (vs. indoors), and increased hygiene. Social distancing was modeled in phases coinciding with major events and policy measures (Figure 1). Phase 1 (March 17–25) corresponds with mid-March policies and increasing public concern regarding COVID-19, phase 2 (March 26–April 26) corresponds with the state-wide stay-at-home order, phase 3 (April 27–May 8) is the period when half the state transitioned to safer-at-home, and phase 4 (May 9–June 3) is the period when safer-at-home was in effect statewide.

We added mask wearing to the model in May (fits 3 and 4) in response to increasing evidence that masks are effective for controlling transmission (17,18). We modeled the effect of mask wearing as a reduction in the spread of infections from asymptomatic and presymptomatic persons to nonhousehold contacts. More recent evidence suggests that masks might also protect the wearer, an added benefit not considered here (19). The effectiveness of mask wearing depends on the ability of the mask to trap infectious particles and the proportion of the population wearing masks. We assume in the model that masks made from household materials are ≈50% effective in trapping infectious particles when worn properly (17,18). Previous studies estimated that ≈23% of contacts occur at home (20). Because persons spent more time at home during the pandemic, we assumed that 67% of the contact of an individual is with nonhousehold members. We assumed that mask wearing was uniform across asymptomatic and presymptomatic persons and fit the proportion of the population wearing masks beginning on April 4, the date of the governor’s press conference advising persons living in Colorado to wear masks. Because some transmission might also occur by fomites, we modeled mask effectiveness as a net 27% reduction in infectiousness among asymptomatic persons wearing masks. In addition, we assume mask wearing decreases transmission by presymptomatic persons (21,22); this is modeled as a 3.4% reduction in infectiousness for symptomatic persons wearing masks (assuming that symptomatic persons are asymptomatic on 1 of 8 infectious days). This model does not account for potential reduction in infectiousness by symptomatic persons who are assumed to isolate (23).

We modeled self-isolation assuming that a proportion of symptomatic case-patients self-isolate 24 hours after the onset of symptoms, and that self-isolation reduces transmission by symptomatic persons to nonhousehold contacts. This assumption is modeled as a 59% reduction in contacts by symptomatic persons who self-isolate. Self-isolation begins in the model on March 5 and the proportion of symptomatic persons who self-isolate is fit to the data.

Estimating Social Distancing and Other Transmission-Reducing Interventions

We inferred the effect of social distancing and other interventions on transmission by using an algorithm-based optimization procedure at 4 different time points from April through June. We used the same approaches to estimate parameters that might vary regionally or for which there was considerable uncertainty in the literature (Table 2). We identified best-fitting parameter values by using a least-squares cost function that minimized difference between the model-estimated and observed number of reported SARS-CoV-2 cases in Colorado (fits 1 and 2) and the observed number of COVID-19 hospitalizations (fits 3 and 4). When fitting to cases (fits 1 and 2), it was necessary to also fit a parameter for the estimated probability that cases would be detected by state surveillance. We minimized the cost function by using a 2-stage fitting algorithm in R (27) and used the FME package (28) by first applying a pseudo-random optimization algorithm (29) to find a region of minimum difference between the model and the data. The second phase used least-squares optimization applying the Levenberg-Marquardt algorithm (30). We calculated 95% CIs by using a Markov Chain Monte Carlo simulation with an adaptive Metropolis algorithm with 1,500 iterations (28). This method obtains 95% CIs by sampling from a Gaussian distribution around the mean trajectory of the ordinary differential equation model. By the end of March, the differential sensitivity matrix was full rank, and thus all 6 parameters were identifiable for all datasets used (Appendix Figure 3).

Because of the median 7-day lag between symptom onset and reporting, on April 3 (fit 1), we included cases that had a symptom onset date through March 26 for model fitting, which enabled us to generate a preliminary estimate of phase 1 social distancing. On April 16 (fit 2), we included cases that had an onset date through April 9 for model fitting, which enabled us to refine estimates of social distancing during phase 1 and generate preliminary estimates of social distancing during phase 2. We generated a preliminary estimate of phase 3 social distancing on May 15 (fit 3) and then re-estimated on June 16 (fit 4), when social distancing during phase 4 was estimated. We estimated the proportion of the population wearing masks in fits 3 and 4 with the effectiveness assumptions defined and estimated the effective reproduction number (Re) from the model output. We additionally estimated the overall number of hospitalizations prevented as a result of decreasing contacts by comparing the fitted model on June 16 with a reference scenario assuming no reduction of the contact rate (social distancing = 0%).

Projections of COVID-19 Hospitalizations and ICU Need

We used the fitted parameters to generate projections of future hospital and critical care needs under different scenarios. Changes in social distancing were implemented beginning 2 weeks after the date of model fitting to account for lags in policy implementation. All other parameters were held fixed as estimated from the model. Preparing for and meeting ICU load was a critical decision-making issue.

Population Mobility

We used SafeGraph (https://www.safegraph.com) data to examine changes in mobility in Colorado from March through June. Specifically, we used an aggregated and anonymized measure of time away from home reported at the census block group. We calculated changes in mobility as a percent decrease in time away from home relative to pre-epidemic baseline: January 29‒February 15 (mean 2.3 hours).

We examined the relationship between mobility and transmission by calibrating the transmission parameter, β, conditional on time away from home at baseline. We then projected the model forward, enabling the parameter for the daily time away from home to change according to observed mobility data. These projections assume no other transmission reducing behavior (i.e., no mask wearing or self-isolation) to avoid conflating parametric assumptions with changes in observed mobility, nor changes to any other aspects of the model.

Results

Estimating Efficacy of Social Distancing and Other Transmission-Reducing Interventions

Figure 2

Observed (black bars) versus model-estimated (green line) number of reported coronavirus disease cases (panels A, B) and hospitalizations (panels C, D), Colorado, USA, 2020, based on models calibrated at 4 points in the early months of the epidemic. Model-based estimates were generated by using an age-structured susceptible-exposed-infected-recovered model, and best-fit parameter values were estimated based on observed data shown. Reported cases are shown by using symptom onset date or report date minus 7 days if onset date was missing, in accordance with onset to report lags for Colorado during this period.

Figure 2. Observed (black bars) versus model-estimated (green line) number of reported coronavirus disease cases (panels A, B) and hospitalizations (panels C, D), Colorado, USA, 2020, based on models calibrated at 4...

Figure 3

Estimated Re over time, Colorado, USA 2020, based on susceptible-exposed-infected-recovered models fit to data at 4 time points in the early months of the epidemic. The reproductive number was estimated from model output at the time of each fit. A) Fit 1 on April 3; B) fit 2 on April 16; C) fit 3 on May 15; D) fit 4 on June 16. Dashed lines indicate an Re of 1, below which the rate of new infections decreases and above which the rate of new infections increases. Re, effective reproductive number.

Figure 3. Estimated Reover time, Colorado, USA 2020, based on susceptible-exposed-infected-recovered models fit to data at 4 time points in the early months of the epidemic. The reproductive number was...

On April 3 (fit 1), we generated a preliminary estimate of social distancing during phase 1, which equated to a ≈45% decrease in the contact rate (Table 2; Figure 2), and Re decreased from 5.3 to 2.4 (Figure 3). Because of the ≈12-day lag between infection and case report, this preliminary estimate was through March 21. On April 16 (fit 2), with more complete case data, the updated estimate of social distancing during phase 1 was greater: a 65% decrease in the contact rate. At this point, we generated a preliminary estimate of the level of social distancing during the first 2 weeks of phase 2 (March 26‒April 4): 76%. Re was estimated to be 0.9. Incorporating increasing data and using COVID-19 hospitalizations instead of case reports, on May 15 (fit 3), we re-estimated phase 1 and phase 2 parameters, which indicated transmission reduction had been more moderate initially (social distancing = 52% for phase 1), and greater for phase 2 (social distancing = 80%). On May 15, which was 18 days after the end of the stay-at-home order, there was no evidence of an increase in hospitalizations or contact rate due to decreased restrictions, noting that because of the ≈13-day lag between infection and hospitalization, this estimate reflects transmission through May 2. A month later on June 16 (fit 4), a greater decrease in transmission was estimated for phase 3 (social distancing = 85%), and the estimated decrease in contact rates during safer-at-home (phase 4) was 90%. Estimated Re decreased to 0.6 during phase 4.

Estimated Reduction in Hospitalizations from Decreased Contacts

As of June 16, a total of 5,272 COVID-19 hospitalizations in Colorado had been reported to CEDRS, and EMResource data strongly suggested underreporting of COVID-19 hospitalizations to CEDRS during April and May. Using CEDRS and EMResource data, we found that the SEIR model estimated a cumulative of 5,344 COVID-19 hospitalizations in Colorado by Jun 16 (Figure 2, panel D). Assuming that no measures had been taken to alter individual behavior or risk perception (social distancing 0% throughout), we estimate that >173,000 hospitalizations would have occurred by that same date, suggesting that ≈97% of potential hospitalizations were avoided as a result of decreases in effective contacts.

Projecting Hospitalizations and ICU Need

Figure 4

Projected coronavirus disease hospitalizations, Colorado, USA, 2020, if current trajectory continued (black line) and for a range of social distancing scenarios (colored lines) generated by models calibrated at 4 time points during April‒June (fit 1: Apr 3; fit 2: April 16; fit 3: May 15; fit 4: June 16). Current trajectory was based on estimated parameters generated for each fit. Social distancing is modeled as a percent reduction in the contact rate (from baseline), and changes in social distancing are introduced 2 weeks after model fitting date. All other fitted parameters are held at the estimated values for each fit. Because peak hospitalization estimates from fit 1 were substantially higher than estimated for later fits, the y-axis is scaled to 50,000 as opposed to 25,000 for fits 2–4. Numbers in parentheses are current values.

Figure 4. Projected coronavirus disease hospitalizations, Colorado, USA, 2020, if current trajectory continued (black line) and for a range of social distancing scenarios (colored lines) generated by models calibrated at 4 time...

Figure 5

Projected coronavirus disease ICU needs, Colorado, USA, 2020, if current trajectory continues (black line) and for a range of social distancing scenarios (colored lines) generated by using models calibrated at 4 time points during April‒ June (fit 1: Apr 3; fit 2: April 16; fit 3: May 15; fit 4: June 16). Current trajectory was based on estimated parameters generated for each fit. Social distancing is modeled as a percent reduction in the contact rate from baseline, and changes in social distancing are introduced 2 weeks after model fitting date. All other fitted parameters are held at the estimated values for each fit. Because ICU estimates from fit 1 were substantially higher than for later fits, the y-axis is scaled to 20,000 as opposed to 10,000 for fits 2–4. Colorado estimated ICU capacity (1,800 beds) is indicated by the dashed gray horizontal line. ICU, intensive care unit. Numbers in parentheses are current values.

Figure 5. Projected coronavirus disease ICU needs, Colorado, USA, 2020, if current trajectory continues (black line) and for a range of social distancing scenarios (colored lines) generated by using models calibrated at...

We provide projected hospitalizations (Figure 4) and ICU need (Figure 5) that were generated from each of the 4 model fits. Fit 1, produced with the least available data, substantially overestimated hospitalizations and ICU need compared with later fits and predicted that ICU capacity would be exceeded even if 80% contact reduction was achieved. As data accumulated and transmission slowed in the state, the estimated peaks under all possible levels of social distancing decreased and shifted later in the year, and projections indicated contact reduction would need to remain at or above ≈70% to prevent exceeding ICU capacity.

Association between Changes in Mobility and Contact Rates

Figure 6

Changes in population mobility before and after emergence of coronavirus disease, Colorado, USA, 2020, and comparison between mobility and estimated social distancing. A) Changes in mobility measured by the number hours spent away from home per day (source: SafeGraph, https://www.safegraph.com). Gray line indicates daily observations, and black line indicates a smoothed line using locally estimated scatterplot smoothing in R (https://www.r-project.org). The ribbon at the bottom indicates the 4 social distancing phases. B) Comparison between susceptible-exposed-infected-recovered model-estimated social distancing (colored boxes) and reduction in mobility relative to the preintervention period, January 29–February 15 (colored lines). Colors correspond to the 4 social distancing phases. Reductions in mobility are calculated as percentage decreases in time away from home relative to pre-epidemic baseline. C) Observed hospitalization data (red) and the simulated hospitalizations based on time away from home relative to a baseline mean during January 29–February 15. In the simulation, it is assumed there is no self-isolation of symptomatic infectious, no mask wearing, and no other transmission reduction to highlight the role of the mobility data in the simulation. P, phase.

Figure 6. Changes in population mobility before and after emergence of coronavirus disease, Colorado, USA, 2020, and comparison between mobility and estimated social distancing. A) Changes in mobility measured by the number...

Residents of Colorado decreased activity outside the home starting in early March (Figure 6, panel A). The trends in mobility data suggest that, on average, the time spent away from home decreased by ≈60% from February to mid-April. Time away from home began to increase in late April, before the end of statewide stay-at-home orders, and increased steadily through June. Mobility metrics initially reinforced estimated social distancing levels, and percent reduction in time away from home coincided with estimated transmission reduction resulting from social distancing (Figure 6, panel B). However, increased time away from home in late April contrasted with estimated infectious contact reduction, which remained high through June. We compared hospitalizations simulated from mobility data to observed and observed a relatively strong association up until late April (Figure 6, panel C). After that, the modeled and observed trends decoupled, indicating that other behaviors or interventions not captured by mobility played a major role in transmission reduction.

Discussion

We used an SEIR transmission model, calibrated to local COVID-19 case and hospital data, to estimate the collective impact of individual behaviors and public policy measures in reducing COVID-19 transmission in Colorado during 2020, providing time-sensitive estimates of the pandemic trajectory to policy-makers to assist in decision-making. COVID-19 policies introduced during March and April were followed by substantial decreases in contact rates and suppression of the Re well below 1, in agreement with other studies of nonpharmaceutical interventions to decrease SARS-CoV2 transmission (3134). These values remained low after the stay-at-home order was lifted and mobility increased.

The continued suppression of transmission might be explained, in part, by transmission control policies that remained in effect and/or were implemented after stay-at-home ended: outbreak prevention and control strategies in long-term care facilities persisted, capacity limits were implemented for businesses and restaurants, and mass gatherings were banned. Moreover, the state reopened slowly during April and May. Given the typical lags between infection and hospitalization, the estimates do not capture the impact of reopening measures implemented during or after late May (e.g., restaurant openings).

Mobility, assessed by using mobile-device data, generally reflected state-level policy during March‒June but mirrored transmission only in the early months of the pandemic. Mobility decreased rapidly in March in concert with early transmission control policies and the statewide stay-at-home order, and mobility increased as social distancing policies were relaxed. Consistent with our findings, others have found that US population mobility was reactive to policy during March: greater perceived disease prevalence and governmental stay-at-home orders resulted in less mobility and social interaction (3537). However, in Colorado, lifting stay-at-home orders and concurrent increases in mobility do not appear to have led to increased transmission, indicating the limitations in using mobility data to predict transmission. These results warrant further investigation in other contexts to help clarify the utility of mobility data in SARS-CoV-2 forecasting, particularly during reopening phases.

Mobility data can be used to estimate when and where persons are congregating, a precondition for transmission, but do not sufficiently capture behaviors such as mask wearing, physical distancing, or moving activities outside. Those individual behaviors can play a critical role in spreading infections, and understanding what drives those behaviors can improve epidemic response. Public perception of and reaction to perceived risk is multifactorial, and, although clearly influenced by policy, in a time of heightened fear, local policy probably plays only a partial role in determining risk-reducing behavior. Rapid and frequent introduction of COVID-19‒related policy measures and public communication by government officials at the national, state, and local scales probably increased fear and public risk perception regarding COVID-19 transmission, and contributed to adoption of transmission-reducing behaviors before the start of and beyond the end of stay-at-home orders (38). Conversely, persons might perceive decreased risks and abandon risk-reducing behaviors when transmission control policies are relaxed, a phenomenon we suspect contributed to an increase in COVID-19 transmission in Colorado during July 2020 (39). Research on how risk perception and behavior fluctuates in relationship to the epidemic trajectory can improve communication and policy making.

Real-time estimation of contact reduction enabled us to respond to urgent requests to actively inform rapidly changing public health policy amidst a pandemic. In early stages, the urgent need was to flatten the curve (Figure 4, panels A, B; Figure 5, panels A, B). Policymakers used initial projections to support decision making on the timing and scope of proposed social distancing policies and to compare potential healthcare need and existing capacity under different scenarios. Once infections began to decrease, there was interest in the degree of increased social contact that could be tolerated as the economy reopened without leading to overwhelmed hospitals (Figure 4, panels C, D; Figure 5, panels C, D). Model estimates were used to evaluate the impact of past policies and to forecast the impact of future proposed interventions, including permutations of layered policies or interventions. Using locally derived estimates enabled policymakers to evaluate potential disease control strategies that were relevant to the current transmission trends in Colorado. For example, as the age distribution shifted within the epidemic in Colorado, estimates with contact rates that varied by age group were produced and used to evaluate policies or interventions targeted to specific age groups, such as return to school and alcohol last call policies.

A key challenge we faced was generating projections of hospital and critical care needs with limited data and rapidly evolving science. Early model estimates were imprecise because data were sparse and poor quality, leading projections to overestimate hospital needs. Estimates and, consequently, projections improved with greater data quantity and quality. Another challenge was the need to generate estimates under extreme time constraints. We designed a preliminary model in a matter of days and adapted it regularly to accommodate new data streams and science. The need for efficiency led to tradeoffs: for example, we did not account for age-specific contact rates (40). We present this material not as a finished work but to illustrate how models can be constructed and adapted in real time to inform critical policy questions.

The model findings were provided weekly to decision-makers in Colorado by written reports and briefings. These collaborative interactions provided an opportunity to review findings and define projection scenarios useful for decision-making. The rapidly developed, locally calibrated transmission model provided timely evidence of a moderate decrease in transmission in Colorado after an early transmission control policy began and a substantial decrease in the contact rate after stay-at-home orders that persisted after these policies were partially relaxed. Decreases in transmission mirrored changes in mobility through the end of stay-at-home, at which point mobility ceased to be a clear proxy for transmission. Locally calibrated models have local credibility and can be used to provide time-sensitive, tailored information to policymakers to assist their decision-making.

At the time of this study, Dr. Buchwald was a research associate at the Colorado School of Public Health, University of Colorado, Aurora, CO. She is currently a research associate at the University of Maryland School of Medicine, Baltimore, MD. Her primary research interests include transmission of asymptomatic infections and vector-borne diseases.

Top

Acknowledgment

This study was supported by the Colorado Department of Public Health and Environment (Contract PO:202000011320-204174) and the Colorado Health Foundation. A.G.B. and E.J.C. also received support from the National Institutes of Health (grant R01AI134673). R.H. was supported by grants from the Centers for Disease Control and Prevention Emerging Infections Program and the Centers for Disease Control and Prevention Epidemiology and Laboratory Capacity during this study.

Top

References

  1. Centers for Disease Control and Prevention. Observed and forecasted new and total reported COVID-19 deaths as of August 10, 2020 [cited 2021 Jun 2]. https://www.cdc.gov/coronavirus/2019-ncov/covid-data/forecasting-us.html#state-forecasts
  2. Bell  D, Nicoll  A, Fukuda  K, Horby  P, Monto  A, Hayden  F, et al.; World Health Organization Writing Group. Non-pharmaceutical interventions for pandemic influenza, national and community measures. Emerg Infect Dis. 2006;12:8894.PubMedGoogle Scholar
  3. Cauchemez  S, Valleron  A-J, Boëlle  P-Y, Flahault  A, Ferguson  NM. Estimating the impact of school closure on influenza transmission from Sentinel data. Nature. 2008;452:7504. DOIPubMedGoogle Scholar
  4. Ferguson  NM, Cummings  DA, Fraser  C, Cajka  JC, Cooley  PC, Burke  DS. Strategies for mitigating an influenza pandemic. Nature. 2006;442:44852. DOIPubMedGoogle Scholar
  5. Fineberg  HV. Rapid expert consultation on social distancing for the COVID-19 pandemic. Washington (DC): National Academies Press; 2020 Mar 19 [cited 2021 Jun 22]. https://www.nap.edu/read/25753/chapter/1
  6. Jackson  C, Vynnycky  E, Hawker  J, Olowokure  B, Mangtani  P. School closures and influenza: systematic review of epidemiological studies. BMJ Open. 2013;3:e002149. DOIPubMedGoogle Scholar
  7. Hsiang  S, Allen  D, Annan-Phan  S, Bell  K, Bolliger  I, Chong  T, et al. The effect of large-scale anti-contagion policies on the COVID-19 pandemic. Nature. 2020;584:2627. DOIPubMedGoogle Scholar
  8. Matrajt  L, Leung  T. Evaluating the effectiveness of social distancing interventions to delay or flatten the epidemic curve of coronavirus disease. Emerg Infect Dis. 2020;26:17408. DOIPubMedGoogle Scholar
  9. Klein  B, LaRock  T, McCabe  S, Torres  L, Friedland  L, Privitera  F, et al. Reshaping a nation: mobility, commuting, and contact patterns during the COVID-19 outbreak, 2020 [cited 2021 Jun 2]. https://www.mobs-lab.org/uploads/6/7/8/7/6787877/covid19mobility_report2.pdf
  10. Buchwald  AG, Adams  J, Bortz  DM, Carlton  EJ. Infectious disease transmission models to predict, evaluate, and improve understanding of SARS-COV-2 trajectory and interventions. Ann Am Thorac Soc. 2020;17:12046. DOIPubMedGoogle Scholar
  11. Colorado Governor Jared Polis. Executive orders, 2020 [cited 2021 Jun 2]. https://www.colorado.gov/governor/2020-executive-orders
  12. Colorado Department of Public Health and Environment. Public health and executive orders, 2020 [cited 2021 Jun 2]. https://covid19.colorado.gov/prepare-protect-yourself/prevent-the-spread/public-health-executive-orders
  13. Davies  NG, Klepac  P, Liu  Y, Prem  K, Jit  M, Eggo  RM; CMMID COVID-19 working group. Age-dependent effects in the transmission and control of COVID-19 epidemics. Nat Med. 2020;26:120511. DOIPubMedGoogle Scholar
  14. Verity  R, Okell  LC, Dorigatti  I, Winskill  P, Whittaker  C, Imai  N, et al. Estimates of the severity of coronavirus disease 2019: a model-based analysis. Lancet Infect Dis. 2020;20:66977. DOIPubMedGoogle Scholar
  15. Ferguson  N, Laydon  D, Nedjati-Gilani  G, Imai  N, Ainslie  K, Baguelin  M, et al. Report 9: impact of non-pharmaceutical interventions (NPIs) to reduce COVID-19 mortality and healthcare demand, 2020. [cited 2021 Jun 2]. https://www.imperial.ac.uk/media/imperial-college/medicine/sph/ide/gida-fellowships/Imperial-College-COVID19-NPI-modelling-16-03-2020.pdf
  16. Williamson  EJ, Walker  AJ, Bhaskaran  K, Bacon  S, Bates  C, Morton  CE, et al. Factors associated with COVID-19-related death using OpenSAFELY. Nature. 2020;584:4306. DOIPubMedGoogle Scholar
  17. Tian  L, Li  X, Qi  F, Tang  QY, Tang  V, Liu  J, et al. Harnessing peak transmission around symptom onset for non-pharmaceutical intervention and containment of the COVID-19 pandemic. Nat Commun. 2021;12:1147. DOIPubMedGoogle Scholar
  18. Howard  J, Huang  A, Li  Z, Tufekci  Z, Zdimal  V, van der Westhuizen  HM, et al. An evidence review of face masks against COVID-19. Proc Natl Acad Sci U S A. 2021;118:e2014564118. DOIPubMedGoogle Scholar
  19. Clapp  PW, Sickbert-Bennett  EE, Samet  JM, Berntsen  J, Zeman  KL, Anderson  DJ, et al.; US Centers for Disease Control and Prevention Epicenters Program. Evaluation of cloth masks and modified procedure masks as personal protective equipment for the public during the COVID-19 pandemic. JAMA Intern Med. 2021;181:4639. DOIPubMedGoogle Scholar
  20. Mossong  J, Hens  N, Jit  M, Beutels  P, Auranen  K, Mikolajczyk  R, et al. Social contacts and mixing patterns relevant to the spread of infectious diseases. PLoS Med. 2008;5:e74. DOIPubMedGoogle Scholar
  21. He  X, Lau  EHY, Wu  P, Deng  X, Wang  J, Hao  X, et al. Temporal dynamics in viral shedding and transmissibility of COVID-19. Nat Med. 2020;26:6725. DOIPubMedGoogle Scholar
  22. Huff  HV, Singh  A. Asymptomatic transmission during the COVID-19 pandemic and implication for public health strategies. Clin Infect Dis. 2020;71:27526. DOIPubMedGoogle Scholar
  23. Leung  NHL, Chu  DKW, Shiu  EYC, Chan  KH, McDevitt  JJ, Hau  BJP, et al. Respiratory virus shedding in exhaled breath and efficacy of face masks. Nat Med. 2020;26:67680. DOIPubMedGoogle Scholar
  24. Models of Infectious Disease Agent Study. MIDAS online COVID-19 portal, 2020 [cited 2021 Jun 2]. https://github.com/midas-network/COVID-19/tree/master/parameter_estimates/2019_novel_coronavirus
  25. Li  R, Pei  S, Chen  B, Song  Y, Zhang  T, Yang  W, et al. Substantial undocumented infection facilitates the rapid dissemination of novel coronavirus (SARS-CoV-2). Science. 2020;368:48993. DOIPubMedGoogle Scholar
  26. Zou  L, Ruan  F, Huang  M, Liang  L, Huang  H, Hong  Z, et al. SARS-CoV-2 viral load in upper respiratory specimens of infected patients. N Engl J Med. 2020;382:11779. DOIPubMedGoogle Scholar
  27. R Core Team. R: A language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2019.
  28. Soetaert  K, Petzoldt  T. Inverse modelling, sensitivity and Monte Carlo analysis in R using package FME. J Stat Softw. 2010;33:128. DOIGoogle Scholar
  29. Price  WL. A controlled random search procedure for global optimisation. Comput J. 1977;20:36770. DOIGoogle Scholar
  30. Moré  JJ. The Levenberg-Marquardt algorithm: implementation and theory. In: Watson GA, editor. Numerical analysis. Heidelberg (Germany): Springer; 1978. p. 105–16.
  31. Flaxman  S, Mishra  S, Gandy  A, Unwin  HJT, Mellan  TA, Coupland  H, et al.; Imperial College COVID-19 Response Team. Estimating the effects of non-pharmaceutical interventions on COVID-19 in Europe. Nature. 2020;584:25761. DOIPubMedGoogle Scholar
  32. Brauner  JM, Mindermann  S, Sharma  M, Johnston  D, Salvatier  J, Gavenciak  T, et al. Inferring the effectiveness of government interventions against COVID-19. Science. 2020.
  33. Courtemanche  C, Garuccio  J, Le  A, Pinkston  J, Yelowitz  A. Strong social distancing measures in the United States reduced the COVID-19 growth rate. Health Aff (Millwood). 2020;39:123746. DOIPubMedGoogle Scholar
  34. Hsiang  S, Allen  D, Annan-Phan  S, Bell  K, Bolliger  I, Chong  T, et al. The effect of large-scale anti-contagion policies on the COVID-19 pandemic. Nature. 2020;584:2627. DOIPubMedGoogle Scholar
  35. Wellenius  GA, Vispute  S, Espinosa  V, Fabrikant  A, Tsai  TC, Hennessy  J, et al. Impacts of state-level policies on social distancing in the United States using aggregated mobility data during the COVID-19 pandemic. Nat Commun. 2021;12:3118. DOIPubMedGoogle Scholar
  36. Engle  S, Stromme  J, Zhou  A. Staying at home: mobility effects of COVID-19. COVID Economics. 2020;4:86–102. https://cepr.org/sites/default/files/news/CovidEconomics4.pdf
  37. Abouk  R, Heydari  B. The immediate effect of COVID-19 policies on social-distancing behavior in the United States. Public Health Rep. 2021;136:24552. DOIPubMedGoogle Scholar
  38. Albright  E, Banacos  N, Birkland  T, Crow  D, DeLeo  R, Dickinson  KL, et al. COVID-19 technical report wave one, June 22, 2020 [cited 2021 June 2]. https://www.riskandsocialpolicy.org/our-work
  39. Yan  Y, Bayham  J, Richter  A, Fenichel  EP. Risk compensation and face mask mandates during the COVID-19 pandemic. Sci Rep. 2021;11:3174. DOIPubMedGoogle Scholar
  40. Mistry  D, Litvinova  M, Pastore Y Piontti  A, Chinazzi  M, Fumanelli  L, Gomes  MFC, et al. Inferring high-resolution human mixing patterns for disease modeling. Nat Commun. 2021;12:323. DOIPubMedGoogle Scholar

Top

Figures
Tables

Top

Cite This Article

DOI: 10.3201/eid2709.204167

Original Publication Date: June 30, 2021

1Current affiliation: University of Maryland School of Medicine, Baltimore, Maryland, USA.

Table of Contents – Volume 27, Number 9—September 2021

Comments

Please use the form below to submit correspondence to the authors or contact them at the following address:

Elizabeth J. Carlton, Department of Environmental and Occupational Health, Colorado School of Public Health, 17001 E 17th Pl, B119, Aurora, CO 80045, USA

Send To

10000 character(s) remaining.

Top

Page created: June 02, 2021
Page updated: August 17, 2021
Page reviewed: August 17, 2021
The conclusions, findings, and opinions expressed by authors contributing to this journal do not necessarily reflect the official position of the U.S. Department of Health and Human Services, the Public Health Service, the Centers for Disease Control and Prevention, or the authors' affiliated institutions. Use of trade names is for identification only and does not imply endorsement by any of the groups named above.
file_external