BackgroundDrug overdose is one of the top leading causes of accidental death in the U.S., largely due to the opioid epidemic. Although the opioid epidemic is a nationwide issue, it has not affected the nation uniformly.MethodsWe combined multiple data sources, including emergency medical service response, American Community Survey data, and health facilities datasets to analyze distributions of heroin-related overdose incidents in Cincinnati, Ohio at the census block group level. The Ripley’s K function and the local Moran’s I statistics were performed to examine geographic variation patterns in heroin-related overdose incidents within the study area. Then, conditional cluster maps were plotted to examine a relationship between heroin-related incident rates and sociodemographic characteristics of areas as well as the resources for opioid use disorder treatment.ResultsThe global spatial analysis indicated that there was a clustered pattern of heroin-related overdose incident rates at every distance across the study area. The univariate local spatial analysis identified 7 hot spot clusters, 27 cold spot clusters, and 1 outlier cluster. Conditional cluster maps showed characteristics of neighborhoods with high heroin overdose rates, such as a higher crime rate, a high percentage of the male, a high poverty level, a lower education level, and a lower income level. The hot spots in the Southwest areas of Cincinnati had longer distances to opioid treatment programs and buprenorphine prescribing physicians than the median, while the hot spots in the South-Central areas of the city had shorter distances to those health resources.ConclusionsOur study showed that the opioid epidemic disproportionately affected Cincinnati. Multi-phased spatial clustering models based on various data sources can be useful to identify areas that require more policy attention and targeted interventions to alleviate high heroin-related overdose rates.