Drug overdose deaths are now the leading cause of mortality for Americans under 50 (
1). The age-adjusted rate of drug overdose deaths continues to rise, increasing by 9.6% from 2016 to 2017. Furthermore, 67.8% of the 70,237 drug overdose deaths in 2017 were attributable to prescription and illicit opioids (
2).
Rates of opioid use disorder have also increased; 2.1 million Americans met diagnostic criteria for the disorder in 2017 (
3). The medications buprenorphine, methadone, and naltrexone—in conjunction with behavioral therapy—are considered the gold standard in treating opioid use disorder (
4–
9). However, relatively few patients receive these medications, and rates of medication availability in the specialty substance use disorder treatment system remain low (
10). In 2017, only 30% of all specialty treatment programs offered either methadone or buprenorphine (
11).
Improving the accessibility and quality of opioid use disorder treatment is imperative, especially for the Medicaid population, which has high rates of both substance use and overdose deaths. However, there is little understanding of the treatment landscape, especially treatment at the local level. Previous research is limited by a lack of recent data on the availability of treatment programs that include geographic location at a local level. Older studies found that approximately 40% of counties in 2009 lacked a treatment program that accepted Medicaid (
12) and that, in 2012, many counties with high rates of opioid misuse lacked opioid treatment programs accepting Medicaid (
13). Recent state-level studies show increases in the percentage of programs accepting Medicaid (
14) and Medicaid-reimbursed treatment admissions (
15).
The geographic distribution of specialty treatment programs—including those that accept Medicaid insurance, offer medication, or both—remains unclear. This is important because living far from the nearest program can present substantial structural barriers to treatment. Barriers may be acute for employed individuals with inflexible or unpredictable work hours, those with childcare responsibilities, those with limited or unreliable transportation options, and those with disabilities that make travel difficult.
We contribute to the study of this question by analyzing a novel data set from the Substance Abuse and Mental Health Services Administration (SAMHSA) that provides the physical location of and services provided by all specialty substance use treatment programs in the United States. We stratified counties by opioid problem severity and provide a longitudinal analysis of county-level change in access to substance use disorder treatment programs from 2009 to 2017. We specifically examined changes—according to local opioid problem severity—in the availability of and distance to all treatment programs; programs offering methadone, buprenorphine, or both; programs that accept Medicaid; and programs that accept Medicaid and offer medication.
Methods
Data
We utilized data from the National Directory of Drug and Alcohol Abuse Treatment Facilities published in a PDF format annually by SAMHSA using data from the National Survey of Substance Abuse Treatment Services (NSSATS). NSSATS is an annual survey of all known substance use disorder treatment programs in the United States, with an average response rate of 92.4%. As such, it represents the census of programs nationwide (
16). We created a data set from the PDFs across the 2009–2017 directories, which includes the physical address of each location where treatment is provided, payment sources accepted, and medications offered. We retained data for all 50 states and the District of Columbia, which comprise 103,090 program-year observations.
The unit of analysis was the county-year. We assigned treatment programs to counties using a zip code–to–county crosswalk (
17). If a zip code spanned more than one county, the program was assigned to the county with the greatest percentage of business activity for that zip code. We aggregated the number of treatment programs to the county level and adjusted the totals to represent the number of programs per 100,000 county residents. Outcome variables were the number of programs overall, programs offering medication (i.e., methadone, buprenorphine, or both), programs accepting Medicaid insurance as payment for services, and programs both offering medication and accepting Medicaid. The final sample had 28,270 county-year observations.
We measured distance as the miles between county centroids (
18). In counties with at least one treatment program with the relevant characteristics (e.g., offers medication), we measured miles to the nearest county as zero. In counties with no such treatment programs, the distance to the closest county with a program was used as the distance to treatment.
For our principal independent variable of interest—county opioid problem severity—we divided counties into terciles of low, moderate, or high severity on the basis of three county-year variables: drug-related mortality rate (
19), opioid prescription rate (
20), and drug-related arrest rate (
21). We imputed missing values using the average from adjoining counties (for deaths and prescribing rates) or the average of the surrounding years’ values (for arrests).
We conducted a factor analysis using the annual county-level prescribing, mortality, and arrest rates. To account for large year-to-year shifts in values due to the infrequency of events such as mortality in small areas, we Winsorized values at the first and 99th percentiles. (The variables loaded on a single factor, with an eigenvalue of .718.) We separated all counties into thirds (low, moderate, or high) on the basis of their factor scores and lagged the severity variable by 1 year.
County-level demographic data on population, race-ethnicity, education levels, urbanicity, median income, and poverty level were obtained from the U.S. Census Bureau’s American Community Survey. Other county-year control variables included the unemployment rate (
22) and physicians per 1,000 population (
23). We also controlled for state policies, including Medicaid expansion (
24), prescription drug monitoring programs (both voluntary and mandatory access programs) (
25–
27), substance use treatment parity laws (
28–
30), pain clinic regulations (
31), and laws making naloxone available without a prescription (i.e., standing orders) (
32,
33). Finally, we controlled for per capita federal substance abuse prevention and treatment block grant dollars (2017 $) allocated to the states, including 2017 funding from the opioid state-targeted response grant.
Statistical Analyses
We present two sets of models estimating the relationship between opioid problem severity and treatment availability. The first set has dependent variables measuring the number of programs per 100,000 county residents for all treatment programs, programs offering medication, programs accepting Medicaid, and programs offering medication and accepting Medicaid. The second set of models estimates the number of miles to the nearest county with a treatment program (defined by the same criteria).
To address the zero-inflated distribution of our data, we employed a two-part model (TPM). In the first part, we used logistic regression to regress an indicator variable (having at least one such program) on our severity measure and control variables. (The first part of the TPM was the same for both outcomes.) A county’s inclusion in the second part of the TPM was conditional on having at least one program (for the outcome measuring number of programs) or on having no program (for the distance-to-treatment outcome). We used ordinary least squares (OLS) to regress the natural log of the number of programs or miles to the nearest program on the covariates. Log transformation of the dependent variables addressed the right-skewness in the values and allowed interpretation of the resulting coefficients as the percentage change in the dependent variable.
All models were conducted in Stata 15, with standard errors clustered at the county level. Each model included the county opioid severity level, the county- and state-level control variables described above, and year fixed effects to track changes in treatment availability over time. Because the analysis was focused on treatment programs, this study was not subject to institutional review board approval.
Results
Table 1 presents summary statistics for the full sample and by opioid severity level. (All statistics are weighted by county population to make them representative; unadjusted statistics are available in the
online supplement.) On average, counties had 3.44 programs offering treatment services per 100,000 people. Divided by severity, county averages ranged between 3.10 and 3.76 programs. The mean distance to the nearest county with a treatment program for the full sample was 1.5 miles. Distance to programs offering methadone, buprenorphine, or both was significantly higher (8.09 miles). Treatment programs that both offered medication and accepted Medicaid were the farthest, at 13.85 miles.
Figure 1 displays the average number of programs in a county per 100,000 population over time, by type of program and opioid problem severity. (Averages are weighted by county population.) The graphs clearly show that compared with counties with low-severity problems, there was more treatment access in counties with moderate-severity problems and greater access still in counties with high-severity problems. In addition, availability increased over time and increased more rapidly in the counties with higher-severity problems.
Figure 2 displays the average travel distance in miles to the nearest program in 2009 and 2017 (also weighted by county population). Because only counties that did not have such a program in 2009 are included, the graph represents counties with the lowest treatment accessibility in 2009. For each program type, there were substantial declines in the distances to treatment over the course of the study, with distances typically dropping by half or more.
Table 2 presents the estimated marginal effects from the first part of our TPM, which used logistic regression to predict changes in the probability of a county having at least one treatment program. (Full model results for all specifications are available in the
online supplement.) Across all outcomes, we found statistically significant and positive associations between the severity of the opioid problem and treatment program availability. Having a moderate-severity problem was associated with a 3–percentage point increase in the probability of having any treatment program (p<0.01). Moderate severity was associated with a 4–percentage point increase in the probability of having a treatment program offering methadone, buprenorphine, or both (p<0.01); a 6–percentage point increase in the probability of having a program that accepts Medicaid (p<0.01); and a 4–percentage point increase in the probability of having a program that offers medication and accepts Medicaid (p<0.01).
The direction and significance of changes in the availability of a treatment program was the same in counties with high-severity problems as in counties with moderate-severity problems. High severity was associated with a 3–percentage point increase in the probability of having any treatment program compared with low-severity counties (p<0.05). High opioid problem severity was also associated with a 7–percentage point increase in the probability of having a program that offers medication (p<0.01), a 5–percentage point increase in the probability of having a program that accepts Medicaid (p<0.01), and a 7–percentage point increase in the probability of having a program that offers medication and accepts Medicaid (p<0.01).
Table 3 presents OLS coefficients estimating the association between opioid severity and the log number of programs per 100,000 population among the subset of counties with at least one program. Differences between counties with low- and high-severity problems were not statistically significant for programs offering medication or for programs offering medication and accepting Medicaid. A significant and negative association was found between opioid problem severity and the log number of programs overall (b=–0.07, p<0.05) and programs offering Medicaid (b=–0.13, p<0.01). Coefficients estimating the association between moderate-severity problems and number of programs were negative and significant across program types. Between 2009 and 2017, there was also a decrease (–0.34, p<0.01) in the number of programs but an increase (0.25, p<0.05) in the number of programs offering medication and accepting Medicaid.
Table 3 displays the OLS coefficients from the second-part models estimating changes in the log distance to treatment (conditional on a county having no program). For all program types, there were no statistically significant differences between counties with low-, moderate-, and high-severity problems. All coefficients comparing distances to programs in 2017 versus 2009 were negative and significant (p<0.01)—all programs, –0.08; programs offering medication, –0.39; programs accepting Medicaid, –0.21; and programs offering medication and accepting Medicaid, –0.57.
Discussion
In this article, we examined the relationship of local opioid problem severity to changes in the geographic distribution of substance use disorder treatment programs over time. To our knowledge, this is the only study providing a longitudinal analysis of the substance use disorder treatment system using recent data at the local level. Furthermore, we know of no other work measuring geographic proximity to treatment. Our findings demonstrate that communities with the most acute need are more likely to have a program available in county, indicating a reduction in the distance required to travel to treatment. After controlling for year, socioeconomic and demographic characteristics, and state policy actions, the analyses showed that high- and moderate-severity problems were associated with modest increases (6.6%, or a 3–percentage point increase from 45.3% probability in low-severity counties) in the probability of having any treatment program, compared with low-severity status.
However, the likelihood of having a higher-quality, evidence-based program was substantially better for counties with higher-severity problems. Methadone and buprenorphine are considered the gold standard for achieving opioid use disorder recovery. Counties with moderate-severity problems were 34.5% (a 4–percentage point increase over the 11.6% average for low-severity counties) more likely than counties with low-severity problems to have at least one program offering these medications. The gains were larger for high-severity counties—a 60.3% increase (a 7–percentage point increase over the average for low-severity counties). (The findings hold for both methadone and buprenorphine when estimated separately [see online supplement].) This improvement represents closer proximity to treatment, given that residents of counties with these programs will not necessarily have to travel to other counties to access care. Closer access is important because it decreases the burden on vulnerable patients, especially those with fewer transportation options or with time constraints related to jobs or caregiving responsibilities.
These findings also reveal significant and meaningful evolution in the treatment landscape over time. The probability of a county having a program that offers medication increased by 14.4% between 2009 and 2017. Moreover, among counties without treatment programs, the distance to the nearest county with a program decreased across all program types between 2009 and 2017. There was a decrease of almost 40% in the distance to treatment for programs offering medication and of almost 60% in the distance to programs offering medication and accepting Medicaid. Strikingly, we did not see significant changes in the likelihood of a county having a program overall between 2009 and 2017, suggesting that existing programs had adopted medications or had been replaced by new programs that did.
These findings are consistent with studies that examine changes in treatment access by measuring the supply of physicians with waivers to prescribe buprenorphine. There is a positive correlation between the number of treatment programs and waivered physicians in a state (
34). The supply of waivered physicians is responsive to certain state and local conditions, including local problem severity (measured as mortality) (
35), and the number of counties with a shortage of waivered physicians fell sharply between 2002 and 2011 (
36).
Access to treatment for Medicaid beneficiaries was better in both counties with high-severity opioid problems and counties with moderate-severity problems compared with counties with low-severity problems. High-severity status was associated with a 13.1% increase in the likelihood of having a program that accepts Medicaid reimbursement (up 5 percentage points from 38.2% probability in low-severity counties) and a 76.9% increase in the likelihood of having a program offering medication and accepting Medicaid (a 7–percentage point increase from a 9.1% probability for low-severity counties). By 2017, counties were 25.3% more likely to have a program accepting Medicaid (a 13–percentage point increase from a baseline of 51.3% in 2009), and the travel distance to such programs for people in counties without one decreased by 21%.
Figure 1 illustrates that treatment availability accelerated around the time of Medicaid expansion through the Affordable Care Act in 2014. Measuring the effect of expansion on local changes in treatment access is a promising area of future research. Indeed, state-level analyses suggest that expansion led to an increase in the availability of programs that accept Medicaid (
24) and in Medicaid-reimbursed buprenorphine prescriptions (
37). Because Medicaid beneficiaries account for a heavy share of opioid use disorder diagnoses, expansion may be an essential aspect of providing high-quality, affordable treatment options.
Although the quality and affordability of opioid use disorder treatment has improved markedly over time, demonstrating responsiveness to the scope of the opioid problem, the full picture is not entirely positive. We observed large and significant improvements in the likelihood of having any treatment program in counties with a high-severity opioid problem, but our findings did not show declines in the distance to treatment in counties with high-severity problems and no program of their own nor did they show increases in the number of programs in counties with high-severity problems and at least one program. (In fact, we showed that fewer programs accepted Medicaid in counties with high-severity problems.)
These results suggest a treatment system that has broadened in geographic scope but not deepened its capacity to fully meet treatment needs. Other studies have documented gaps between rates of opioid abuse and treatment capacity (
38–
40). Most patients with opioid use disorder are not actively in treatment (
41). A recent cross-sectional review of 2016 data suggests the continued existence of areas with both high opioid-related death rates and low treatment accessibility (
42). Insufficient treatment access can arise from a variety of factors—affordability concerns, regulatory hurdles, stigma, and geographic barriers—and a multifaceted approach will be required in order to align treatment availability with the scale of the opioid problem (
43).
This study had several limitations of note. First, there is no publicly available county-level measure of opioid use disorder. We created a proxy measure using county opioid prescribing rates and drug-related death and arrest rates. Second, this study did not examine access to treatment outside the specialty addiction treatment system (e.g., buprenorphine and naltrexone prescribed by office-based providers or federally qualified health centers); therefore, treatment access is likely to be even better than reported. Longitudinal data on office-based providers were not available. Third, we did not examine access to oral or injectable naltrexone, two additional opioid treatment medications approved by the U.S. Food and Drug Administration. NSSATS directories do not distinguish between naltrexone used for the treatment of alcohol use disorder and opioid use disorder. Fourth, measuring distance by county centroids may underestimate or overestimate actual travel distance to treatment if centroids do not align with population concentrations.
Conclusions
This study is the first to provide a comprehensive view of recent changes in the treatment landscape for substance use disorders at the local level. Although our findings provide strong evidence that the substance use disorder treatment system is responsive to local opioid problem severity, these results should not be interpreted to suggest that there is sufficient treatment access, especially in communities hardest hit by the opioid crisis. Across the United States, people report problems locating the type of treatment they need and can afford. The observed adaptability of treatment systems, however, implies that programs may respond to incentives to provide higher-quality treatment where it is most needed. Policy makers should deploy an array of “carrots and sticks” to promote growth and improvements in the treatment system. Such efforts may entail higher and less restrictive Medicaid and Medicare reimbursement for medications, reduced reimbursement to programs that do not offer medication, Medicaid expansion in the remaining nonexpansion states, and other efforts to encourage adoption of medication by treatment programs.
Acknowledgments
The authors thank the following individuals for their insightful suggestions and feedback on this project: Jason M. Hockenberry, Janet R. Cummings, and Adam S. Wilk.