Urban-rural variation in the socioeconomic determinants of opioid overdose

Prescription opioid overdose (POD) and heroin overdose (HOD) rates have quadrupled since 1999. Community-level socioeconomic characteristics are associated with opioid overdoses, but whether this varies by urbanicity is unknown.

In this serial cross-sectional study of zip codes in 17 states, 2002-2014 (n = 145,241 space-time units), we used hierarchical Bayesian Poisson space-time models to analyze the association between zip code-level socioeconomic features (poverty, unemployment, educational attainment, and income) and counts of POD or HOD hospital discharges. We tested multiplicative interactions between each socioeconomic feature and zip code urbanicity measured with Rural-Urban Commuting Area codes.

Percent in poverty and of adults with ≤ high school education were associated with higher POD rates (Rate Ratio [RR], 5% poverty: 1.07 [95% credible interval: 1.06-1.07]; 5% low education: 1.02 [1.02-1.03]), while median household income was associated with lower rates (RR, $10,000: 0.88 [0.87-0.89]). Urbanicity modified the association between socioeconomic features and HOD. Poverty and unemployment were associated with increased HOD in metropolitan areas (RR, 5% poverty: 1.12 [1.11-1.13]; 5% unemployment: 1.04 [1.02-1.05]), and median household income was associated with decreased HOD (RR, $10,000: 0.88 [0.87-0.90]). In rural areas, low educational attainment alone was associated with HOD (RR, 5%: 1.09 [1.02-1.16]).

Regardless of urbanicity, elevated rates of POD were found in more economically disadvantaged zip codes. Economic disadvantage played a larger role in HOD in urban than rural areas, suggesting rural HOD rates may have alternative drivers. Identifying social determinants of opioid overdoses is particularly important for creating effective population-level interventions.