Prepared by the Colorado COVID-19 Modeling Group
Colorado School of Public Health: Andrea Buchwald, Elizabeth Carlton, Debashis Ghosh, Irina Kasarskis, Jonathan Samet, Laura Timm, Emily Wu; University of Colorado School of Medicine: Kathryn Colborn; University of Colorado-Boulder Department of Applied Mathematics: Sabina Altus, David Bortz; University of Colorado-Denver: jimi adams; Colorado State University: Jude Bayham

For Contact:

Summary

Key messages in this report are:

Introduction

This report provides the results of epidemic models for regions of Colorado, using methods similar to that used for the state-level model. Estimates are presented for the 11 Local Public Health Agencies (LPHAs) regions in the state and for 8 selected counties with populations that are sufficiently large to allow for county-level estimates. The model results are subject to greater uncertainty than those for the entire state because there are fewer hospitalizations and cases in each region than in the state as a whole. Estimates are most uncertain for the regions with the smallest population size. We use the model as well as COVID-19 hospital and case data to generate three measures for each region. These measures can be used to gauge the current state of SARS-CoV-2 in each region.

Table 1. The estimated effective reproductive number, prevalence of infections and the percent of the population recovered to date by region. These metrics are estimated using hospitalization data from the Colorado COVID Patient Hospitalization Surveillance (COPHS) through 01/08/2021. In regions with smaller populations, reported cases are also used to generate these estimates.
Are infections increasing or decreasing?
How many people are infectious?
How many people have been infected to date?
Re Infections are… Prevalence per 100,000 People infectious Cumulative Infections to Date Proportion of population infected to date
LPHA Regions
Central 1.0 Flat 827 1 in 121 114,000 14.1
Central Mountains 0.8 Decreasing 532 1 in 188 18,400 10.1
East Central 1.2 Increasing 2,768 1 in 36 18,700 43.5
Metro 0.9 Decreasing 806 1 in 124 686,000 20.8
Northeast 0.5 Decreasing 744 1 in 134 186,000 24.3
Northwest 1.2 Increasing 763 1 in 131 24,200 11.9
San Luis Valley 0.6 Decreasing 1,678 1 in 60 9,790 21.1
South Central 0.7 Decreasing 1,925 1 in 52 71,700 29.5
Southeast 0.8 Decreasing 2,451 1 in 41 14,400 30.7
Southwest 0.9 Decreasing 746 1 in 134 10,500 10.3
West Central Partnership 0.4 Decreasing 316 1 in 316 7,850 7.4
Eight select counties
Adams 1.2 Increasing 1,280 1 in 78 179,000 33.9
Arapahoe 0.9 Decreasing 1,065 1 in 94 162,000 24.4
Boulder 0.7 Decreasing 637 1 in 157 39,300 11.9
Broomfield 0.9 Decreasing 367 1 in 273 7,690 10.6
Denver 0.9 Decreasing 1,216 1 in 82 232,000 31.4
Douglas 0.7 Decreasing 516 1 in 194 37,900 10.7
El Paso 0.9 Decreasing 831 1 in 120 115,000 15.5
Jefferson plus 0.3 Decreasing 185 1 in 540 60,600 10.1
Due to the small population sizes of Gilpin and Clear Creek counties, these counties are combined with Jefferson County. Jefferson County comprises 97% of the population in the Jefferson plus county cluster.
Due to lags between infection and hospitalization, the estimated effective reproductive number (Re) reflects the spread of infections approximately two weeks prior to the data of the last observed hospitalization.

Effective Reproductive Number

The figure below shows the estimated effective reproductive number for each region since March.

The effective reproduction number (Re) is a measure of how rapidly infections are spreading or declining in a region at a given point in time. When the effective reproductive number is below 1, infections are decreasing. When the effective reproductive number is above 1, infections are increasing.

The effective reproductive number is estimated using our age-structured SEIR model fit to hospitalization data. In the four LPHA regions with smaller populations, reported SARS-CoV-2 case data are also used (San Luis Valley, Southeast, Southwest, and West Central Partnership). Because we base our parameter estimates primarily on COVID-19 hospitalization data, and hospitalizations today generally reflect infections occurring approximately 13 days prior, our most recent estimates of the effective reproductive number likely reflect the spread of infections occurring on approximately 12/26/2020.

Figure 1. The estimated effective reproductive number (Re) over time in the 11 LPHA regions in Colorado, and 8 selected counties and county clusters. Estimates shown using COVID-19 hospitalization data through 01/08/2021.

Infection prevalence

Infection prevalence provides an estimate of the proportion of the population that is currently infected with SARS-CoV-2 and capable of spreading infections. At higher levels of infection prevalence, individuals are more likely to encounter infectious individuals among their contacts. Because many people experience no symptoms or mild symptoms of COVID-19, many infections are not identified by surveillance systems. The estimates we present here are intended to provide an approximation of all infections, including those not detected by the Colorado Electronic Disease Reporting System (CEDRS).

The figure below shows the estimated infection prevalence per 100,000 individuals for each region. These are estimated from SEIR models fit separately to each area’s reported data.

Figure 2. Estimated prevalence per 100,000 population for each of the 11 LPHA regions (top), plus the 8 selected counties and county clusters (bottom). All prevalence values over 1,000 per 100,000 are shown in dark red. Prevalence values estimated up to 01/08/2021.

The percent of the population recovered from infections to date

As more people develop immunity, due to vaccination or prior infection, the spread of infections slows because infectious individuals are less likely to encounter individuals that are not immune. At present, immunity to SARS-CoV-2 is incompletely understood and a vaccine is not yet available.

The figure below shows model-generated estimates of the percent of the population that has been infected and is now recovered to date for each region. This provides an estimate of the percent of the population that may be immune, although we still do not know how long immunity lasts after an infection. As a vaccine becomes available and our understanding of SARS-CoV-2 immunity changes, these estimates will be updated.

Figure 3. Estimated proportion of the population recovered to date for each of the 11 LPHA regions in Colorado (top) and each of the 8 selected counties and county clusters (bottom). Exposed proportion values estimated up to 01/08/2021. Black dashed line indicates mean of Colorado (top) and selected counties (bottom).

COVID-19 hospitalizations

The figures below show the daily number of individuals hospitalized with COVID-19 from each region. Hospitalization data are from the COVID Patient Hospitalization Surveillance (COPHS) maintained by the Colorado Department of Public Health and the Environment (CDPHE). Each COVID-19 patient is assigned to a region based on their home zip code. COVID-19 hospitalizations are shown per 100,000 population to allow comparability across regions.

COVID-19 hospitalizations are a sensitive measure of SARS-CoV-2 transmission and are an important indicator of the severity of infections in a region. While many SARS-CoV-2 infections are not captured by state surveillance systems, we expect that almost all COVID-19 hospitalizations are identified.

Figure 4. The daily number of people hospitalized with COVID-19 per capita for the 11 LPHA regions and the 8 selected counties and county clusters in Colorado over the past 12 weeks. Hospitalization data are from the COPHS hospital census data through 01/08/2021.

Figure 5. The daily number of people hospitalized with COVID-19 hospitalizations per capita for the 11 LPHA regions and the 8 selected counties and county clusters in Colorado since the first case was reported in March 2020. Hospitalization data are from the COPHS hospital census data through 01/08/2021.

Technical Summary

Within the State of Colorado, as for other states, the spread of SARS-CoV-2 varies across the regions of the state, differing between urban and rural locations and resort and non-resort areas, for example. The modeling carried out by the Colorado COVID-19 Modeling Group was initiated at the state level providing a picture that does not give detail at the regional or county-level. Such detail is needed for public health planning and action, and many local public health agencies have asked for model simulations for their jurisdictions.

To provide the needed detail, the Modeling Group has developed regional -level models, using approaches similar to those for the overall Colorado model. The models are generated for the 11 local public health agency (LPHA) regions, which cover the full state and for selected counties. County-level estimation is possible for selected counties due to the large population size. In providing the results for regions within Colorado, the model results are subject to greater uncertainty than those for the entire state because there are fewer hospitalizations and cases in each region than in the state as a whole. Estimates are most uncertain for the regions with the smallest population size. Additional details about the methods used to generate regional estimates are described in the Technical Summary, below.

Methods. We use data on COVID-19 hospitalizations and reported cases, and a mathematical model of SARS-CoV-2 transmission to estimate the current state of infections in each region. The approach follows that used for the state-wide model, adapted for the smaller population sizes of the LPHA regions and selected counties. A full description of the model and estimation approaches are provided in the documentation, available here. Prior modeling reports and documentation can be found here. Briefly, the model is a deterministic age-structured susceptible, exposed, infected, recovered (SEIR) model. It has been parameterized to Colorado-specific data whenever possible - for example, the length of time a COVID-19 patient spends in the hospital varies by age and is based on data provided by Colorado hospitals. We use model fitting approaches to estimate the level of transmission control for each two-week period of the outbreak. Transmission control is estimated by fitting the model to hospitalization data for each region. In the four regions where population size is small (San Luis Valley, Southeast, Southwest, and West Central Partnership), we fit the model to reported case data, using information from the state-level model to infer the proportion of infections detected by state surveillance systems.

Data Sources: COVID-19 hospitalizations are obtained from COPHS reported through 01/12/2021.Due to lags in reporting, these data are cleaved 4 days prior such that hospitalizations through 01/08/2021 are used in modeling and shown in this report.

Case data are based on CEDRS through 01/13/2021. Due to lags in reporting, these data are cleaved ten days prior on 01/03/2021. Both data sources are provided by CDPHE.

Limitations: These estimates are based on the available data and the assumptions in the SEIR model. Assumptions in the model are based on Colorado data when available and the current scientific understanding of SARS-CoV-2, which is evolving rapidly. Estimated prevalence and the percent of the population recovered are sensitive to model assumptions, which include: the probability an infected individual will be symptomatic and require hospital care and the estimated length of hospital stay, which varies over time and by age. Estimates for the smallest regions are subject to greater uncertainty than larger regions and may experience greater variation from week to week due to the limited data available in these regions.

Region Population Counties
LPHA Regions
Central 810,420 Chaffee, El Paso, Lake, Park, Teller
Central Mountains 182,689 Eagle, Garfield, Grand, Pitkin, Summit
East Central 43,032 Cheyenne, Elbert, Kit Carson, Lincoln
Metro 3,291,794 Adams, Arapahoe, Boulder, Broomfield, Clear Creek, Denver, Douglas, Gilpin, Jefferson*
Northeast 765,265 Larimer, Logan, Morgan, Phillips, Sedgwick, Washington, Weld, Yuma
Northwest 203,301 Jackson, Mesa, Moffat, Rio Blanco, Routt
San Luis Valley 46,472 Alamosa, Conejos, Costilla, Rio Grande, Saguache
South Central 243,196 Custer, Fremont, Huerfano, Las Animas, Pueblo
Southeast 46,938 Baca, Bent, Crowley, Kiowa, Otero, Prowers
Southwest 102,154 Archuleta, Dolores, La Plata, Montezuma, San Juan
West Central Partnership 106,839 Delta, Gunnison, Hinsdale, Mineral, Montrose, Ouray, San Miguel
Eight select counties
Adams 528,857 Adams
Arapahoe 664,988 Arapahoe
Boulder 330,978 Boulder
Broomfield 72,827 Broomfield
Denver 737,854 Denver
Douglas 354,331 Douglas
El Paso 737,354 El Paso
Jefferson plus 601,959 Clear Creek, Gilpin, Jefferson

Table 2: Description of the regions used in this report, including the 11 LPHA Regions. Population estimates are based on 2020 US Census Projections provided by the Colorado Demography Office.

*Clear Creek, Gilpin Counties and Jefferson counties are modeled as a single unit due to the small population of Clear Creek (9,379 residents) and Gilpin (5,924 residents) – populations too small to allow for stable estimation.

Figure 6. Map showing the 11 LPHA regions and 8 selected counties and county clusters for which estimates were generated.

Figure 7. Comparison of daily active COVID-19 hospitalizations by data source.

*EMResource recently began including VA hospitals and other institutions that may not be captured in current COPHS data sources, potentially explaining the under estimation of total hospitalizations from COPHS.

*The COPHS data were updated in late December to include records that had previously been unreported. At present, the COPHS dataset now estimates more hospitalizations on a given day than EMResources. Because both datasources show the same trend (e.g., both increase and decrease together), this should not impact our estimates of the effective reproductive number. However, the addition of previously unreported records will lead estimates of infection prevalence to be higher in some regions on a given day than they were previously. We have updated the estimates of prevalence since March and these can be viewed in the region and county-specific estimates, allowing apples to apples comparison of changes in prevalence over time.