Spatiotemporal Analysis of Fentanyl-Associated Overdose Deaths in Chicago, IL, USA

Overdose deaths involving fentanyl represent a major public health crisis in the USA. This study investigates the spatiotemporal dynamics of fentanyl-involved deaths before, during, and after the COVID-19 pandemic and examines how sociodemographic factors influence these deaths across geographic regions. Using a retrospective ecological approach, we analyzed data on ZIP code-level fentanyl-related deaths in Cook County, IL, between 2018 and 2023, obtained from the Medical Examiner's Office and linked with sociodemographic data from the American Community Survey. We first mapped area-level death rates to assess their distribution and then conducted global and local clustering analyses to identify spatial autocorrelations and the locations of high- or low-death-rate areas. A geographically weighted Poisson regression (GWPR) model evaluated the associations between area-level fentanyl-related death rates and the area-level proportion of young adults, males, and individuals with at least a college degree, disability rate, and poverty rate. Spatial analyses found stronger spatial autocorrelations during (2020-2021) and after (2022-2023) the pandemic. Initially, high death rates were concentrated in the downtown area of Chicago, and they expanded to the surrounding areas during and after the pandemic. The GWPR model revealed that an increase in the area-level proportions of poverty, disability, and young adult residents increased the fentanyl-related death rates in most of the areas. Our findings highlight the urgent need to address the evolving dynamics of fentanyl-related overdoses through tailored public health interventions that account for the unique socioeconomic determinants of different regions. Importantly, a comprehensive approach to addressing differences in overdose death rates and their risk factors will be crucial to mitigating this public health crisis.