Spatial and temporal trends in social vulnerability and COVID-19 incidence and death rates in the United States

Abstract Background Emerging evidence suggests that socially vulnerable communities are at higher risk for coronavirus disease 2019 (COVID-19) outbreaks in the United States. However, no prior studies have examined temporal trends and differential effects of social vulnerability on COVID-19 incidence and death rates. The purpose of this study was to examine temporal trends among counties with high and low social vulnerability and to quantify disparities in these trends over time. We hypothesized that highly vulnerable counties would have higher incidence and death rates compared to less vulnerable counties and that this disparity would widen as the pandemic progressed. Methods We conducted a retrospective longitudinal analysis examining COVID-19 incidence and death rates from March 1 to August 31, 2020 for each county in the US. We obtained daily COVID-19 incident case and death data from USAFacts and the Johns Hopkins Center for Systems Science and Engineering. We classified counties using the Social Vulnerability Index (SVI), a percentile-based measure from the Centers for Disease Control and Prevention in which higher scores represent more vulnerability. Using a Bayesian hierarchical negative binomial model, we estimated daily risk ratios (RRs) comparing counties in the first (lower) and fourth (upper) SVI quartiles. We adjusted for percentage of the county designated as rural, percentage in poor or fair health, percentage of adult smokers, county average daily fine particulate matter (PM2.5), percentage of primary care physicians per 100,000 residents, and the proportion tested for COVID-19 in the state. Results In unadjusted analyses, we found that for most of March 2020, counties in the upper SVI quartile had significantly fewer cases per 100,000 than lower SVI quartile counties. However, on March 30, we observed a “crossover effect” in which the RR became significantly greater than 1.00 (RR = 1.10, 95% PI: 1.03, 1.18), indicating that the most vulnerable counties had, on average, higher COVID-19 incidence rates compared to least vulnerable counties. Upper SVI quartile counties had higher death rates on average starting on March 30 (RR = 1.17, 95% PI: 1.01,1.36). The death rate RR achieved a maximum value on July 29 (RR = 3.22, 95% PI: 2.91, 3.58), indicating that most vulnerable counties had, on average, 3.22 times more deaths per million than the least vulnerable counties. However, by late August, the lower quartile started to catch up to the upper quartile. In adjusted models, the RRs were attenuated for both incidence cases and deaths, indicating that the adjustment variables partially explained the associations. We also found positive associations between COVID-19 cases and deaths and percentage of the county designated as rural, percentage of resident in fair or poor health, and average daily PM 2.5 . Conclusions Results indicate that the impact of COVID-19 is not static but can migrate from less vulnerable counties to more vulnerable counties over time. This highlights the importance of protecting vulnerable populations as the pandemic unfolds.