Social Determinants of Mortality of COVID-19 and Opioid Overdose in American Rural and Urban Counties.

OBJECTIVES: Both COVID-19 deaths and opioid overdose deaths continue to increase in the United States. Little is known about the characteristics of counties with high rates of mortality for both.
METHODS: We analyzed county-level data on COVID-19 mortality from January 1 to May 31, 2020, and on opioid overdose mortality during 2014-2018. The outcome variable, "high-risk county" was a binary indicator of high mortality rates (above 75% quartile) for both COVID-19 and opioid overdose. We conducted geospatial logistic regression models separately for urban and rural counties to identify social determinants of health associated with being a high-risk county.
RESULTS: After adjusting for other covariates, the overall mortality rate of COVID-19 is higher in counties with larger population size and a higher proportion of racial/ethnic minorities, although counties with high rates of opioid overdose mortality have lower proportions of racial/ethnic minorities, a higher proportion of females, and are more economically disadvantaged. Significant predictors of rural counties with high mortality rates for both COVID-19 and opioid overdose include higher proportions of Blacks (Adjusted odds ratio [aOR], 1.04; 95%CI, 1.01-1.07), American Indians and Alaska Natives (aOR, 1.07; 95%CI, 1.02-1.13), and two or more races (aOR, 1.34; 95%CI, 1.13-1.60). Additional predictors for high-risk urban counties include population density (aOR, 1.12; 95%CI, 1.04-1.22) and higher unemployment rates during the COVID-19 pandemic (aOR, 1.23; 95%CI, 1.07-1.41).
CONCLUSIONS: Rural counties with high proportions of racial/ethnic minorities and urban counties with high unemployment rates are at high mortality risk for COVID-19 and opioid overdose.

Although the United States is still in the midst of the COVID-19 pandemic, the nation's opioid crisis is on the rise again. After preliminary success to bring down opioid overdose deaths from 47,885 in 2017 to 47,096 in 2018, there were 50,312 in 2019, and actual numbers could be worse in 2020.[1,2] More than 40 states have reported increases in opioid-related deaths during the COVID-19 pandemic.[3] Both COVID-19 and opioid overdose are deadly. By the end of May 2020, more than 100,000 individuals in the U.S. died from COVID-19.[4] Approximately 450,000 people died from an overdose involving any opioid between 1999 and 2018.[5]

Thus far, there is no effective medication to treat COVID-19. Although there are approved medications for the treatment of opioid use disorder (MOUD), the availability of MOUD has been limited in many U.S. counties, particularly in rural areas. COVID-19 could worsen the opioid crisis as COVID-19 shelter-in-place orders have pushed individuals battling recovery into isolation and have reduced access to treatment.[6] Although federal and state policies have made quick adjustments to allow easier access to MOUD during this public health emergency, socio-economic conditions and healthcare infrastructure within counties are slow to change. Social determinants of health (SDOH), which are conditions in the environments that affect people's health risks and outcomes, provide a useful context to help characterize the factors that contribute to health disparities in communities.[7] Hence, the present study aims to identify SDOH for rural and urban U.S. counties with high mortality rates for COVID-19 and opioid overdose during the initial phase of the COVID-19 pandemic.

METHODS

Study Population

We analyzed characteristics associated with U.S. rural (n = 1,976) and urban (n = 1,166) counties with high mortality rates of COVID-19 and opioid overdose using geospatial cross-sectional analysis in 2020. We obtained county-level data including COVID-19 mortality from January 1 to May 31, 2020, and opioid overdose mortality during 2014–2018 from the Center for Disease Control and Prevention (Appendix definitions of primary outcomes, available online). Other county-level data were from several sources (Appendix Table 1, available online). Institutional review board approval was not required for this study, which used only publicly available data.

Measures

Separately examining rural and urban counties, we defined counties with rates of both COVID-19 and opioid overdose mortality in the above 75% quartile as “high risk” counties, and counties with rates of both below 75% quartile are “low risk” counties. The county-level SDOH includes demographics, education, income, unemployment rate, population size, population density, urban-rural location, and availability and access to healthcare professionals (Appendix Table 2, available online).

Statistical Analysis

First, we fit 4 generalized additive models with a thin-plate spline term of spatial coordinates[8] to identify county-level factors associated with COVID-19 and opioid overdose mortality rates of rural and urban U.S. counties. Adding county latitude and longitude into the model was to account for the spatial trend. Second, we modeled the probability of being a high-risk county as a function of all county-level covariates using generalized additive models with a logit link function and spatial coordinates among rural and urban counties, respectively. Significance was defined as P < 0.05. We carried out all analyses in R and SAS software.

RESULTS

Table 1 shows the results of four separate regression models investigating characteristics associated with mortality rates of COVID-19 and opioid overdose. Rural and urban counties with higher rates of COVID-19 mortality had larger population size, higher proportions of racial/ethnic minorities, and higher median household income. Rural and urban counties with higher rates of opioid overdose mortality had lower proportions of Blacks, lower proportions of males, and lower median household income.

TABLE 1

County-Level Factors Associated with Rates of COVID-19 and Opioid Overdose Deaths (per 100,000 individuals) in U.S. Rural Counties (1,976) and Urban Counties (1,166)

	Regression Coefficient (95% CI)†
	COVID-19		Opioid Overdose
	Rural	Urban	Rural	Urban
Population characteristics
Population (10,000s)	0.90 (0.33, 1.47) ∗∗	0.04 (0.01, 0.08) ∗	–0.02 (–0.20, 0.17)	0.00 (–0.01, 0.02)
Population density per square mile (100s)	–0.16 (–1.30, 0.99)	0.33 (0.28, 0.39) ∗∗∗	0.06 (–0.31, 0.43)	–0.01 (–0.02, 0.01)
Race/(Ethnicity), %
(Hispanic or Latinx)	–0.02 (–0.17, 0.12)	0.31 (0.04, 0.58) ∗	0.01 (–0.04, 0.06)	–0.06 (–0.14, 0.03)
Black	0.65 (0.55, 0.76) ∗∗∗	0.81 (0.65, 0.97) ∗∗∗	–0.09 (–0.12, –0.05)∗∗∗	–0.06 (–0.11, 0.00)∗
American Indian and Alaska Native	0.22 (0.07, 0.37) ∗∗	0.94 (0.43, 1.45) ∗∗∗	0.04 (–0.01, 0.09)	–0.13 (–0.29, 0.03)
Asian	0.58 (–0.28, 1.44)	1.05 (0.38, 1.73) ∗∗	–0.09 (–0.37, 0.19)	–0.13 (–0.34, 0.08)
Native Hawaiian and Other Pacific islanders	–3.34 (–7.78, 1.10)	3.49 (–2.12, 9.10)	–1.25 (–2.71, 0.22)	1.49 (–0.29, 3.28)
Other races	3.95 (–0.05, 7.95)	23.23 (15.10, 31.36) ∗∗∗	0.83 (–0.46, 2.11)	2.40 (–0.17, 4.96)
Two or more races	–0.15 (–0.87, 0.56)	–1.94 (–3.44, –0.45)∗	0.32 (0.09, 0.56) ∗∗	0.83 (0.35, 1.30) ∗∗∗
Male, %	–0.25 (–0.70, 0.21)	0.29 (–0.65, 1.23)	–0.20 (–0.35, –0.06)∗∗	–0.39 (–0.69, –0.09)∗∗
Age, %
25 to 64	–0.03 (–0.44, 0.38)	–1.30 (–1.94, –0.67)∗∗∗	0.33 (0.20, 0.46) ∗∗∗	0.42 (0.21, 0.62) ∗∗∗
≥ 65	0.18 (–0.13, 0.48)	0.69 (0.20, 1.18) ∗∗	0.06 (–0.04, 0.15)	0.01 (–0.14, 0.17)
Education, %
High school	–0.07 (–0.39, 0.26)	–0.34 (–1.09, 0.41)	0.02 (–0.08, 0.13)	0.13 (–0.11, 0.36)
Some college or higher	–0.25 (–0.52, 0.02)	–1.03 (–1.57, –0.50)∗∗∗	–0.05 (–0.14, 0.04)	0.04 (–0.13, 0.21)
Socioeconomic characteristics
Median household income ($1,000s)	0.23 (0.06, 0.39) ∗∗	0.56 (0.39, 0.73) ∗∗∗	–0.08 (–0.13, –0.03)∗∗	–0.08 (–0.14, –0.03)∗∗
Unemployed, %	0.01 (–0.25, 0.27)	0.88 (0.43, 1.32) ∗∗∗	0.03 (–0.06, 0.11)	0.17 (0.03, 0.31) ∗
Health care characteristics
Physicians per 10,000 population	–0.06 (–0.19, 0.08)	–0.16 (–0.27, –0.04)∗∗	–0.02 (–0.06, 0.02)	0.04 (0.00, 0.07) ∗
Healthcare professional shortage area score	0.07 (–0.16, 0.30)	0.00 (–0.26, 0.25)	0.14 (0.06, 0.21) ∗∗∗	0.04 (–0.04, 0.12)

CI = 95% confidence interval. Boldface indicates statistical significance.

We combined the data from multiple county-level databases. The total number of counties is 3,142. Generalized additive models were conducted using information from 3,034 counties; 108 counties were excluded because of missing data (Missing 2 values among population density, one among median household income, one among the percentage of unemployment, and 105 among healthcare professional shortage area score).

P < 0.05.

P < 0.01.

P < 0.001.

In Table 2, factors associated with being at greater risk for both high mortality rates of COVID-19 and opioid among rural counties were: a larger population (Adjusted odds ratio [aOR], 1.17; 95% CI, 1.03–1.33 per 10,000 increment); higher proportions of Blacks (aOR, 1.04; 95%CI, 1.01–1.07), American Indians and Alaska Natives (aOR, 1.07; 95%CI, 1.02–1.13), or two or more races (aOR, 1.34; 95%CI, 1.13–1.60); and a higher proportion of the population ≥ 65 years of age (aOR, 1.12; 95%CI, 1.02–1.12). Urban counties with higher population density were associated with greater odds of being high risk (aOR, 1.12; 95%CI, 1.04–1.22). Higher proportions of Native Hawaiians and Other Pacific islanders (aOR, 0.00; 95%CI, 0.00–0.10) or other races (aOR, 0.15; 95%CI, 0.03–0.90) were inversely associated with increased odds of being high-risk urban counties. A higher unemployment rate during the COVID-19 pandemic increased the odds of being a high-risk urban county (aOR, 1.23; 95%CI, 1.07–1.41).

TABLE 2

County-level Factors Associated with County with High Mortality Rates of both COVID-19 and Opioid Overdose among U.S. Rural and Urban Counties

	Odds Ratio (95% CI)†,‡
	Rural (n = 1,182)	Urban (n = 789)
Population characteristics
Population (10,000)	1.17 (1.03, 1.33) ∗	1.01 (1.00, 1.02)
Population density per square mile (100)	1.06 (0.83, 1.36)	1.12 (1.04, 1.22) ∗∗
Race/ Ethnicity, %
Hispanic or Latinx	1.02 (0.96, 1.09)	0.98 (0.87, 1.10)
Black	1.04 (1.01, 1.07) ∗	1.03 (0.97, 1.09)
American Indian and Alaskan Native	1.07 (1.02, 1.13) ∗∗	0.86 (0.40, 1.88)
Asian	0.54 (0.25, 1.14)	0.95 (0.73, 1.24)
Native Hawaiian and Other Pacific islanders	0.68 (0.07, 7.07)	0.00 (0.00, 0.10) ∗
Other races	1.12 (0.36, 3.50)	0.15 (0.03, 0.90) ∗
Two or more races	1.34 (1.13, 1.60) ∗∗∗	1.23 (0.71, 2.13)
Male, %	0.93 (0.81, 1.07)	0.78 (0.52, 1.19)
Age, %
25 to 64	1.12 (0.98, 1.29)	1.03 (0.84, 1.28)
≥ 65	1.12 (1.02, 1.23) ∗	1.09 (0.93, 1.27)
Education, %
High school	1.04 (0.94, 1.14)	0.85 (0.67, 1.06)
Some college or higher	1.02 (0.94, 1.12)	0.88 (0.74, 1.05)
Socioeconomic characteristics
Median household income ($1,000s)	0.99 (0.94, 1.04)	1.03 (0.98, 1.08)
Unemployed, %	0.95 (0.88, 1.02)	1.23 (1.07, 1.41) ∗∗
Health care characteristics
Physicians per 10,000 population	1.00 (0.97, 1.03)	0.99 (0.97, 1.02)
Healthcare professional shortage area score	1.03 (0.97, 1.10)	0.99 (0.92, 1.06)

CI = 95% confidence interval. Boldface indicates statistical significance.

High mortality of COVID-19 and opioid counties are defined as those with mortality rates of COVID-19 and opioid overdose both above 75% quartile (Rural = 97 vs Urban = 103).

The reference group of generalized additive models was counties with rates of COVID-19 and opioid overdose both below 75% quartile (Rural = 1,085 vs Urban = 686).

P < 0.05.

P < 0.01.

P < 0.001.

DISCUSSION

Our study findings revealed distinctive community-level factors associated with COVID-19 and opioid overdose mortality rates for rural and urban U.S. counties. Higher proportions of racial/ethnic minorities increased the odds of being “high risk” rural counties, whereas higher unemployment rates during the COVID-19 pandemic increased odds of being “high risk” urban counties.

Our analyses indicate the existence of population-level disparities and confirmed that racial/ethnic minorities are more impacted by COVID-19. Furthermore, having higher proportions of racial/ethnic minorities was associated with being a “high risk” rural county. Future studies should investigate underlying factors and ecological influences that may contribute to this disparity (eg, lack of health care access, or other indirect health determinants).[9] We observed a positive association between median household income and COVID-19 mortality rate. When we divided median household income into four groups, counties below 25% and above 75% of median household income had higher average rates of COVID-19 mortality compared with the remaining counties. The potential reason could be that those with low income are more affected by unprotected workplace and dense housing, whereas those with high income could be more affected by international and vacation travel in the early stage of the pandemic.

Previous studies have shown positive associations between unemployment and overdose.[10] Unemployment causes adverse changes in health behaviors (eg, substance use), which in turn leads to deterioration of health.[11] Alternatively, pre-existing health conditions may be a common cause of both unemployment and mortality,[12] which could be even more devastating among the poor. During the COVID-19 pandemic, U.S. unemployment reached 13.3% of the workforce in May 2020[13] and is expected to remain elevated. Efforts to mitigate COVID-19 related socioeconomic consequences and to reduce mortality should involve welfare and healthcare systems (eg, expansion of Medicaid coverage)[14] that helps prevent risk-taking behaviors, including substance use leading to overdose mortality.

Limitations

This study has several limitations. First, we used measures at the county level, which cannot identify within-county differences. Ecological study potentially suffers from the ecological fallacy, which limits the extension of findings to individuals within counties. Second, the cross-sectional study design limited the ability to make causal inferences. Given the continuing pandemic, factors influencing subsequent trends may change. Finally, residual confounding bias may exist because of unmeasured confounders.

CONCLUSIONS

The present study identified racial/ethnic minority and unemployment as upstream factors of health in rural and urban counties, indicating that those vulnerable populations were more severely affected by the COVID-19 pandemic and the opioid overdose crisis. Society needs to consider that the unintended consequences of strategies to reduce COVID-19 impacts may further exacerbate the opioid overdose crisis through increased social and economic marginalization.[15] The study findings could encourage policymakers and program leaders to take geographic needs into consideration when developing practical responses to combat these public health crises.