Skip to main content

Predictors of historical change in drug treatment coverage among people who inject drugs in 90 large metropolitan areas in the USA, 1993–2007



Adequate access to effective treatment and medication assisted therapies for opioid dependence has led to improved antiretroviral therapy adherence and decreases in morbidity among people who inject drugs (PWID), and can also address a broad range of social and public health problems. However, even with the success of syringe service programs and opioid substitution programs in European countries (and others) the US remains historically low in terms of coverage and access with regard to these programs. This manuscript investigates predictors of historical change in drug treatment coverage for PWID in 90 US metropolitan statistical areas (MSAs) during 1993–2007, a period in which, overall coverage did not change.


Drug treatment coverage was measured as the number of PWID in drug treatment, as calculated by treatment entry and census data, divided by numbers of PWID in each MSA. Variables suggested by the Theory of Community Action (i.e., need, resource availability, institutional opposition, organized support, and service symbiosis) were analyzed using mixed-effects multivariate models within dependent variables lagged in time to study predictors of later change in coverage.


Mean coverage was low in 1993 (6.7%; SD 3.7), and did not increase by 2007 (6.4%; SD 4.5). Multivariate results indicate that increases in baseline unemployment rate (β = 0.312; pseudo-p < 0.0002) predict significantly higher treatment coverage; baseline poverty rate (β = − 0.486; pseudo-p < 0.0001), and baseline size of public health and social work workforce (β = 0.425; pseudo-p < 0.0001) were predictors of later mean coverage levels, and baseline HIV prevalence among PWID predicted variation in treatment coverage trajectories over time (baseline HIV * Time: β = 0.039; pseudo-p < 0.001). Finally, increases in black/white poverty disparity from baseline predicted significantly higher treatment coverage in MSAs (β = 1.269; pseudo-p < 0.0001).


While harm reduction programs have historically been contested and difficult to implement in many US communities, and despite efforts to increase treatment coverage for PWID, coverage has not increased. Contrary to our hypothesis, epidemiologic need, seems not to be associated with change in treatment coverage over time. Resource availability and institutional opposition are important predictors of change over time in coverage. These findings suggest that new ways have to be found to increase drug treatment coverage in spite of economic changes and belt-tightening policy changes that will make this difficult.


A key pillar of public health planning is that the magnitude of a response needs to match the magnitude of a problem. In the United States (US) and its large metropolitan areas, however, despite repeated calls for expansion of drug treatment, treatment coverage for people who inject drugs did not increase overall during the period 1993–2007 and continues to fall far short of need [1,2,3,4].

A second pillar of public health is that the distribution of programs across geographic areas should reflect the geographic distribution of need for that program. However, research suggests that local need for a wide variety of different types of programs fails to predict local program presence or coverage. Friedman and colleagues, for example, have found that syringe service programs (SSPs) and drug treatment coverage for people who inject drugs (PWID) varies greatly across metropolitan statistical areas (MSAs), and that local need does not predict these variations [4,5,6,7]. Rather, the political influence of men who have sex with men (MSM) is associated with more program coverage for PWID, and government budget limitations (i.e., long term debt per capita) predict less coverage [4,5,6]. Need also does not predict the presence or coverage of other health and social service programs (e.g., programs against drunk drivers or smoking) [8,9,10,11,12,13,14,15,16]. Such research suggests that the presence and coverage of public and social service programs is influenced by local policy environment, and not by local need. In order to adequately address the elimination of new HIV transmissions and response to epidemiological need in low resource areas, community-based organizations need the ability to rapidly change and assemble new prevention services to meet the challenge of changing epidemiology, population demographics, and advances in technology, or policy/political imperatives.

Low treatment coverage for PWID may produce a high cost to society in terms of the spread of HIV, hepatitis B and C and other infectious diseases among injectors, their partners, and the broader community [17, 18]. Evidence-based drug treatment such as methadone maintenance therapy and buprenorphine can address a broad range of social and public health problem valued in communities affected by PWID [18, 19]. Adequate access to effective treatment and medication assisted therapies for opioid dependence has led to a decrease in HIV transmission, improved ART adherence and decrease in morbidity and mortality not only for opioid overdose but also HIV/AIDS related disease [18, 20,21,22,23,24]. Research is needed to address what policy and structural changes affect variations and changes in treatment coverage - and, in particular, what combinations of factors lead to increases in treatment coverage.

A previous paper showed that the magnitude of drug treatment coverage for people who inject drugs did not increase in large US metropolitan areas over the 15 years, 1993–2007 [1,2,3,4]. Here, we study whether the increases and decreases in coverage among various US metropolitan areas seemed to respond to the need in those metropolitan areas. Thus, this paper presents historic trends and predictors of change in drug treatment coverage for PWID in 90 US MSAs during 1993–2007. Our drug treatment sample for calculating treatment coverage includes clients enrolled in residential or ambulatory inpatient/outpatient care, detoxification services, and methadone maintenance therapy at publicly- and privately-funded substance abuse agencies receiving public funds. Coverage was measured as the number of PWID in drug treatment, calculated by using data from the Substance Abuse and Mental Health Service Administration, divided by numbers of PWID in each MSA.

Additionally, we present theory-based predictors of metropolitan treatment coverage rates as a function of program need, resource availability, institutional opposition, organized support, and service symbiosis, factors which may contribute to greater coverage of drug treatment programs and coverage. The present study extends our research on the predictors of drug treatment coverage for PWID to include longitudinal data. Understanding which metropolitan characteristics are related to changes in treatment coverage can assist public health policy planners, treatment providers and grassroots organizations in improving access to treatment and in facilitating its spread in areas of need.

Theoretical framework and selection of predictors

Here, we present a theoretical framework for predicting program presence using the “theory of community action” (TCA). Figure 1 displays the conceptual model of TCA. This framework utilizes concepts from urban studies [25,26,27] social movement theory [28,29,30] and diffusion of innovations theory [31,32,33,34]. It has been used to identify a variety of place characteristics that are likely to affect the extent to which a community carries out and sustains an action. Rosser & Horvath for example, found that successful rural HIV prevention was less likely in states with more religious and Evangelical Protestant adherents and more successful in states with more “gay community” infrastructure [35]; and providing condom distribution interventions or programs [20, 36], counseling and testing services [37, 38], or specialized programs to reduce opiate misuse [3,4,5,6,7, 39, 40] or treatment for people who are mentally ill [41, 42]. As previously conceived [3,4,5,6,7] this framework emphasizes five types of place characteristics as important to the distribution and implementation of institutionalized programs in cities or MSAs. We define each of these domains below:

Fig. 1
figure 1

Theory of Community Action Framework and Contextual Factors Predicting Program Presence

In this analysis, need refers to rates of epidemiologic factors related to drug use, injection drug use, and HIV (e.g., such as PWID AIDS cases per 10,000 population). The idea that need leads to response suggests that MSAs with greater epidemiologic need will respond to such need by providing more access to services.

Resource availability refers to the extent to which assets that exist locally are useful in initiating and implementing a program. Higher levels of local resource availability are hypothesized to be associated with higher rates of service provision. Local resources can be categorized as either general or specialized. General resources might support a wide range of activities [3,4,5,6,7] and may be positively associated with a number of social and health programs that may have an impact (either directly or indirectly) on programs for PWID. General resources may include economic prosperity in an MSA, the financial health of local governments, and the percent of the local population who are college educated [5, 6]. Specialized resources support the development of only a narrow range of programs, such as breast cancer programming, mental health programs, or health insurance for children [8, 12, 13]. For the purposes of this research, specialized resources will be defined as resources that support HIV and AIDS programs for PWID. These resources include the skilled labor force coming from medical and public health schools; hospitals with special HIV/AIDS units and staff; and emergency substance abuse services, including HIV/AIDS counselors. Levels of specialized resources are expected to affect key outcomes and the presence of a substance abuse treatment program. In addition, general resources such as a larger, more educated labor pool may increase the efficiency of service provision [5, 6, 37].

Organized and potential institutional opposition to programs (even presumed opposition) can hinder or prevent formation of a program. Strong organized opposition can kill a proposal or result in a program operating only for a short time. Drug treatment services for PWID in US communities remain controversial and face ongoing obstacles from law enforcement and local communities, including ‘not-in-my-backyard’ community opposition [37, 39, 43, 44]. As such, organized opposition is hypothesized by the present research to predict declines in, or lower levels of, drug treatment coverage. Institutional opposition to programs may include “legal repressiveness” as a competing strategy for maintaining social order [39, 43,44,45,46,47] by controlling drug use. Institutional opposition to programs may also include the mixing of drug policy and racial subordination by the enactment of harsher penalties for possession of drugs among ethnic communities [43,44,45,46,47].

Previous research has shown that racial economic and political structures may adversely affect the health of a variety of racial/ethnic groups as well as reducing the resources and services available to them, including local services for drug users [37, 39, 40, 43,44,45,46,47,48]. Racial structures such as greater inequality between racial groups, for example, might contribute to, or be reflective of, local environments that are more institutionally oppressive to some groups (i.e. racial minorities), and therefore less likely to provide sufficient resources and services that meet the needs of these groups. Thus, we hypothesize that structural conditions of racial/ethnic inequality may be associated with stronger opposition to drug treatment, and thus may be associated with less drug treatment coverage.

The strength of institutional opposition has often contributed to greater delays in implementation or total lack of implementation of HIV prevention programs for PWID and other public health intervention programs for drug users in U.S. cities and metropolitan areas [37,38,39,40, 43, 44, 47, 48]. Research by Shaw, for example (2006) [40] found that strong community opposition, including negative and stigmatizing attitudes toward drug users, was crucial in the failure to establish public harm reduction programs in Springfield, MA. Similar results suggesting that opposition can weaken or prevent harm reduction programs has been found by Tempalski et al., [6, 43, 44]; Downing et al., [37] and Des Jarlais et al., [47, 48]. Likewise, opposition has been identified as a hindrance to services and programs for the homeless, to housing programs for the mentally ill, to half- way houses for people living with AIDS, and to health-related facilities for persons with AIDS [14, 15, 49,50,51,52,53,54].

Organized and potential institutional support for programs. Theoretically converse to the notion of opposition, organized support for programs is hypothesized by the present research to predict increases in, or higher levels of, drug treatment coverage. Indicators of organized or potential support can come from the presence of outreach efforts, ACT UP chapters, local American Public Health Association units, or Planned Parenthood strength. Organized or potential support for a program can determine its presence, its size, and its longevity [5,6,7]. Such support can come from either of two kinds of sources: institutionalized sources like public health departments, or “mobilization from below” in the form of social movements. Institutional support from health officials, medical and public health researchers, church leaders, law enforcement, and/or school administrators can provide support and outreach for a program [4,5,6,7, 37, 39, 40, 43]. Institutional support can also come in the form of funding for public policies and programs. We can think of this as a general willingness of government to spend money on social services for the public good. Here, we hypothesize that drug treatment is generally a “liberal” response [4, 47, 48] to the problem of drug misuse, and that financial or other support from policy and other institutional sources for drug treatment programs therefore reflects a higher degree of “liberalism” than does a lack of such support. Mobilization from below: Social movements, ranging from the feminist health movement to AIDS activism, have helped to shape programs, and policies in the U.S. (such as treatment services, health care reform, and AIDS policy), and have also helped to shape public opinion about health-related issues (such as the de-stigmatization of groups, including PWID) [5, 6, 9, 14, 53,54,55,56]. As an example of the potential influence of social movements on programs and policies for PWID specifically, there has been support for the hypothesis that gay political influence and concern among the gay community regarding HIV/AIDS prevention (including the direct involvement of gay and lesbian activists in ACT UP) helped form and sustain harm reduction programs for PWID in the US [5,6,7].

Service symbiosis occurs when the increased or newly introduced presence of one type of service is associated with increases in other services; we have added this domain to the standard TCA model. For the present research, specifically, the presence of syringe exchange programs is hypothesized to be associated with higher levels of other drug treatment. This hypothesis is based on reports that syringe exchange programs provide many referrals to treatment [28, 47, 48, 55].

Material and methods

Unit of analysis and sample

The unit of analysis in this study is the MSA. The US Census Bureau and Office of Management and Budget define an MSA as a set of contiguous counties that include one or more central cities of at least 50,000 people that collectively form a single cohesive socioeconomic unit, defined by inter-county commuting patterns and socioeconomic integration [57]. The MSA was selected as the unit of analysis because data were readily available at this geographic level and because it is posited that MSAs are meaningful epidemiologic units with which to study injectors and services designated for them [3, 4]. Also, it is appropriate to include counties that comprise MSAs, as opposed to including only counties containing central cities, given that drug-related epidemics travel from central cities to their surrounding suburbs, as injectors often live in suburbs but buy drugs and perhaps receive drug-related social services in the central city [3, 4].

The sample of MSAs included in the present study was obtained by selecting all MSAs (N = 96) in the United States which had a population greater than 500,000 in 1993. Six of these MSAs (Gary, IN; Hartford, CT; New Haven-Bridgeport-Danbury, CT; Phoenix-Mesa, AZ; San Juan, PR; and Tucson, AZ) are missing from the present analyses because they did not report treatment data necessary to estimate treatment coverage, resulting in a sample of 90 MSAs for the present study.

Because this is a study of 90 MSAs with populations of 500,000 or more in 1993 that had data available on our key variables, our sample is a fully enumerated universe. This means there is no sampling error and that p-values are not meaningful as estimates of the probability of arriving at estimates based on chance introduced by selecting a sample. Nonetheless, we report statistical significance as a heuristic guide to the importance of variables in our equations. We compute them as if we had a random sample of MSAs, but report results as “pseudo-p-values” to guide our interpretation (as in previous articles: [19, 22]. Thus, for purpose of this paper, we used pseudo-p < 0.05 as a heuristic criterion.

Dependent variable: calculating drug treatment coverage

We calculated treatment coverage rates for each year from 1993 to 2007 (excluding years 1994, 1999, and 2001 due to data missingness) using information from two databases from the Substance Abuse and Mental Health Service Administration (SAMHSA) [58,59,60,61], and estimates of PWID from previous research [62]. Table 1 describes each database utilized to calculate treatment coverage rates. We define treatment coverage as the ratio of PWID in treatment to PWID in the MSA. Treatments included in our coverage estimates are residential or ambulatory inpatient/outpatient care, detoxification services and methadone maintenance therapy at publicly- and privately-funded drug treatment agencies receiving public funds. These are facilities licensed, certified, or otherwise approved by State treatment agencies to provide substance use treatment.

Table 1 Description of Data Sources Utilized to Calculate Drug Treatment Coverage Rates

Treatment coverage for PWID is estimated using TEDS and UFDS/N-SSATS. We use both TEDS and UFDS/N-SSATS to calculate treatment coverage to maximize validity and reliability of our estimates. Each of our data set differs in counts of drug treatment clients. TEDS counts each admission in a given year. Therefore, an individual admitted to treatment twice in a calendar year is counted as two admissions which inflates annual treatment entries, but only produces bias in the proportion of entrants who are PWID to the extent that such double-counting varies systematically by route of administration. In contrast, UFDS/N-SSATS is a one day census of treatment.

The following equation calculates drug treatment coverage rateFootnote 1:

$$ Ajt=\left( Djt\ast \left( Bjt/ Cjt\right)\right)/ Ejt\ast 100 $$


Ajt = treatment coverage rate for an MSA j in year t

Bjt = number of PWID entering drug treatment as reported by TEDS for an MSA j in year t

Cjt = number of PWID and number of non-injectors entering drug treatment as reported by TEDS for an MSA j in year t

Djt = number of drug users entering drug treatment reported by UFDS/N-SSATS for an MSA j in year t

Ejt = estimated number of PWID as estimated by Tempalski et al. 2013 [62] for an MSA j in year t.

First, the TEDS data series identifies the number and attributes of clients who enter substance use treatment programs that receive any state and federal funding. From TEDS, we calculated the proportion of treatment entrants who reported they injected drugs as a mode of administration. Our second SAMHSA data source comes from the annual census of drug treatment facilities originally referred to as the UFDS – but since renamed the N-SSATS. UFDS/N-SSATS data measure client characteristics and use of privately- and publicly-funded substance use treatment programs in the U.S. on October 1 for each year. However, UFDS/N-SSATS data were unavailable for 1992, 1994, 1999, and 2001. As a result of this limited availability, our coverage estimates were only created for years where data were available. Thus, our final drug treatment coverage estimates only provide data for 1993, 1995, 1996–1998, 2000, and 2002–2007.

Calculating number of PWID

Because estimation of the total numbers of injectors is discussed in detail elsewhere [62], it is described only briefly here. Tempalski and colleagues first estimated the number of PWID in the US each year from 1992 to 2007 and then apportioned these estimates to MSAs using multiplier methods. Four different types of data indicating drug injection were used to allocate national annual totals to MSAs, creating four distinct series of estimates of the number of injectors in each MSA. These estimates rely on using (1) HIV counseling and testing data from the Centers for Disease Control (CDC) [63]; (2) SAMSHA’s UFDS and TEDS data [58,59,60,61]; (3) CDC’s diagnoses of PWIDs with HIV/AIDS [63]; and (4) an estimate derived from published estimates of the number of injectors living in each MSA in 1992 [64] and in 1998 (3). Each series was smoothed over time using loess regression and the mean value of the four component estimates was taken as the best estimate of PWID for that MSA and year. In order to avoid circularity, the estimated numbers of PWID in the population used in this study modify the Tempalski estimates [62] so that they do not rely on data on the numbers of PWID in drug treatment from SAMSHA.

Independent variables

Data at the MSA-level were available on a range of variables measuring the theoretically supported domains of program need, resource availability, institutional opposition, organized support, and service symbiosis. Table 2 describes the statistical distribution of all of these independent variables across MSAs and within each theoretical domain. Additional file 1: Table S1 depicts the bivariate correlations among all variables.

Table 2 Statistical description of independent variables across MSAs & across all years for which outcome data was available

Analytic approach

We utilize a series of mixed-effects models [65] to examine trends in drug treatment coverage across the study period and to test all study hypotheses. This method used maximum likelihood estimation to assess the associations of interest while adjusting for variance shared within MSAs across time.


Where possible, we included time lags in our measurement strategy to ensure that our independent variables had time to affect treatment coverage. We measured treatment coverage for 1993–2007. Thus, we measured each independent variable before 1993. We chose three-year lags in most cases (e.g., demographic and economic variables collected in the 1990 US Census). Our inclusion of time lags also reflects the likelihood that many of these variables change slowly (and that therefore any change in the outcome which might occur as a response to or in accordance with changes in other setting characteristics could take years), and acknowledges the time required to create or change treatment programs. There are also three variables (presence of SSPs; Right-to-work State; and Number of types of “early” groups) which were only measured once, in 1993, and for which change over time was not assessed. To facilitate interpretation of intercepts and of the effects of predictors which interact with time, we centered independent variables at the first year for which we measured them.

Mixed-effects models

Growth curve models were utilized first in order to assess the nature of the relationship between time and treatment coverage. Linear, quadratic, and cubic functions for time were modeled in this “univariate” first step in order to assess the functional form of change in treatment coverage over time. Next, in order to select the most empirically relevant set of independent predictors of treatment coverage from the large number (relative to the number of MSAs) of theoretically-relevant potential independent variables, we developed a four-step process. Each stage in this four-step process utilized a logarithm-transformed version of the dependent variable to address the non-normal distribution of the treatment coverage variable.

Step 1: bivariate model section

First, we conducted bivariate analyses to determine which independent variables might be associated with treatment coverage. For each of the potential independent variables, separately, we used mixed-effects models [66] to assess the strength and nature of its relationship to treatment coverage. We developed three models appropriate for understanding potentially nuanced relationships varying in a curvilinear manner over time in a multilevel framework:

$$ \mathrm{Model}\kern0.28em 1\Big)\kern0.28em \hat{Y}=A+B+C\kern0.28em Time+ Tim{e}^2 $$
$$ \mathrm{Model}\;2\Big)\;\hat{Y}=A+B+C+ Tim e+ Tim{e}^2+B\ast Tim e+B\ast Tim{e}^2 $$


$$ \mathrm{Model}\;3\Big)\;\hat{Y}=A+B+C+ Tim e+ Tim{e}^2+C\ast Tim e+C\ast Tim{e}^2 $$

A = Intercept

Ŷ = predicted treatment coverage

B = baseline values of each potential independent variable

C = change in independent variable baseline at each stage.

We then compared each of these three nested models for each potential independent variable, and selected the model with the “best fit” for each construct based on Akaike’s Information Criterion (AIC) [67, 68]. The “best model” was chosen based on the following criterion: Model 1 was the default “best” model, unless the AIC for Model 2 or Model 3 was lower than that for Model 1 by at least 2, in which case the model with the lowest AIC was chosen. This criterion was utilized as a mechanism for ensuring that interactions were only included if they improved model fit by a meaningful margin.

Step 2: bivariate analyses

In the second step, we ran the “best fit” model for each construct using standardized variables (z-scores), and compared standardized coefficients from these “best fit” models for each potential independent variable to determine eligibility for entry into the next step of analysis (domain analysis). The somewhat standard use of pseudo-p-values as a criterion for bivariate selection of independent variables into multivariate models was not possible due to the nature of the “best fit” mixed models, which included multiple parameters that together contributed to the ability of each construct as a whole to meaningfully explain variation in treatment coverage. Standardized coefficients, however, are an appropriate measure of effect size [69] which can serve as indicators of the relative explanatory importance of each variable in predicting treatment coverage. Given our desire to consider, as a whole, the ability of the multiple parameters composing each construct to predict our outcome, we summed the standardized coefficients from the multiple parameters composing each construct. Considering Ferguson’s [69] recommendation of a 0.2 minimum effect size for “strength of association” measures, including standardized coefficients, we used the criterion that the sum of the absolute value of the standardized coefficients from the best fit model for the construct (not including the coefficients for the Intercept or for Time or Time2, but including the coefficients for B, for C, and for any applicable interactions) must be equal to or greater than 0.25. A lower criterion of 0.20 was applied to constructs for which only one or two coefficients were included in its “best fit” model, either due to Model 1 (with no interactions) being selected, or due to the construct being time-invariant and therefore not having change scores included in its model. Because standardized versions of all variables were used in these analyses, the comparison of their standardized coefficients from models predicting treatment coverage should serve as a reasonable comparison of their relative explanatory importance in understanding variation in treatment coverage.

Step 3: domain analyses

In the third step, we selected only independent variables which met the criterion applied to the bivariate analyses in Step 2, and ran a set of mixed-effect models, each of which included the “best fit” models for each of the eligible constructs in a specific theoretical domain. The same criteria used in Step 2 for bivariate analyses (>.25 or > .20 sum of model coefficients) were then applied to the results of these domain analyses to determine eligibility for inclusion of each construct in the final analytic model. This step allowed us to limit potential multicollinearity by identifying the most empirically important predictors of treatment coverage from a set of highly conceptually related variables and eliminating the rest.

Step 4: multi-domain analyses

In the next step, the “best fit” models from all domains which met the eligibility criterion in Step 3 were included into a single mixed-effects multi-domain model predicting treatment coverage, to estimate the relationships of each eligible independent variable to treatment coverage, net of the predictive influence of all other eligible independent variables. Finally, to achieve the most parsimonious multi-domain model, model parameters which did not meaningfully contribute to the multi-domain prediction of treatment coverage were identified for removal using the following process: model AIC was compared among versions of the multi-domain model which systematically and individually removed either a) interactions with Time2; b) non-“significant” interactions; or c) constructs for which neither the baseline nor change score were “significant” predictors of treatment coverage. The model with the fewest parameters and lowest AIC was selected. Constructs with a “significant” coefficient at either their baseline or change score were not removed from the original multi-domain model, and no parameters were removed for which removal resulted in a > 2.0 increase in AIC. This process allowed us to assess the importance to the model of including each interaction and each non- “significant” construct. All analyses were conducted using SAS software. Mixed-effects models were conducted using PROC MIXED in SAS [66].


Descriptive statistics

A comparison of the beginning and end points of the study period reveals little variation in treatment coverage over time. Coverage overall was very similar in 1993 and 2007 (Additional file 1: Table S2). Mean coverage was only 6.4% (SD = 4.5) in 2007, which was quite similar to the estimated 6.7% coverage in 1993 (SD = 3.7). Median treatment coverage was also quite similar in 1993 (5.6%) and in 2007 (5.2%) among injection drug users in 90 US MSAs.

Growth curve model

Mixed-effects growth curve models were utilized to examine the nature of average changes in treatment coverage between the 1993 and 2007 time points (linear, quadratic, and cubic). The model for a quadratic function of time contained “significant” coefficients for both linear time (β = 0.43; S.E. =0.10; pseudo-p < 0.0001) and quadratic time (β = − 0.03; S.E. = 0.01; pseudo-p < 0.0001), suggesting that there is a “significant” curvilinear trend for change in treatment coverage over time, averaging across MSAs. When modeled as a quadratic function of time (see Fig. 2), treatment coverage increases across MSAs, on average, from 1993 to 2000, and then decreases across MSAs, on average, from 2001 to 2007. The curvilinear increase in treatment coverage rising in the 90s and peaking in 2000 might have been due to the rise in nationwide opiate overdoses which may have driven some increase in treatment availability during this period [70, 71].

Fig. 2
figure 2

Quadratic Growth Curve for Treatment Coverage, 1993–2007

Bivariate and multivariate models

Table 3 shows the results of both the bivariate analyses and the domain-specific analyses. A logarithm transformation was used on the dependent variable for these and all subsequent models to address the non-normal distribution of the treatment coverage variable. Ten out of seventeen potential constructs met the eligibility criterion based on bivariate analysis to be included in the domain specific analysis. Of the ten constructs included in the domain analyses, nine met the criterion to be included in the multi- domain model.

Table 3 Standardized Predictors in Bivariate and Domain-Specific Analyses

Additional file 1: Table S3 shows the results of the first multi-domain model which includes all qualifying independent variables from across the domain analyses. Table 4 shows the results of the final multi-domain model, from which parameters were removed based on

Table 4 Results from Parsimonious Multi-Domain Mixed-Effects Model

the process of AIC comparisons described above. Findings from the final model in Table 4 are described below, followed by a brief comparison of coefficients from the two multi- domain models. In the final model both Time (β = 0.169; pseudo-p < 0.0001) and Time2 (β = − 0.008; pseudo-p = 0.001) were found to be statistically significant predictors of treatment coverage, net of the effects of all other independent variables in the model.

From the need domain, HIV prevalence among PWID was included in the final model. Higher HIV prevalence among PWID at baseline was found to predict significantly lower treatment coverage (β = − 0.360; pseudo-p = 0.0003), on average. The interaction of HIV at baseline with Time was also statistically significant (β = 0.039; pseudo-p = 0.001), indicating that baseline values of HIV significantly predict variation in the trajectory of treatment coverage over time.

From the resource availability domain, percent of population in poverty; unemployment rate; and community, public health, and social work workforce were included in the final model. Both higher unemployment rates at baseline (β = 0.312; pseudo-p = 0.0002) and increases from baseline in unemployment (β = 0.407; pseudo-p = 0.002) were found to predict significantly higher treatment coverage, on average. In addition to a significant interaction with Time, the interaction between change in unemployment and Time2 was statistically significant (β = 0.008; pseudo-p < 0.0001), indicating that change in unemployment rate from baseline significantly predicted variation in treatment coverage trajectories over time.

Although change in poverty rate from baseline was not found to be a significant predictor of treatment coverage (β = − 0.040; pseudo-p = 0.357), higher poverty rate at baseline was found to predict significantly lower levels of treatment coverage, on average (β = − 0.486; pseudo-p < 0.0001). Larger community, public health, and social work workforces at baseline were also found to predict significantly higher treatment coverage levels, on average (β = 0.425; pseudo-p < 0.0001).

From the institutional opposition domain, only the black/white poverty disparity construct was included in the final model. Although the level of disparity in poverty rates among Black and White households at baseline was not significantly predictive of variation in treatment coverage (β = − 0.138; pseudo-p = 0.058), on average, increases in black/white poverty disparity from baseline predicted significantly higher treatment coverage among MSAs, on average (β = 1.269; pseudo-p < 0.0001). In addition to its interaction with Time, the interaction of change in Black/White poverty disparities from baseline with Time2 was also statistically significant (β = 0.006; pseudo-p = 0.009), suggesting that change in Black/White poverty disparities significantly predicted variation in treatment coverage trajectories over time. None of the variables from either the institutional support or the service symbiosis domains met all criteria for inclusion in the final model.

A comparison of this final model to the models in S3 (which is the less parsimonious model including all variables that were originally eligible for multi-domain analysis based on the results of domain analyses) illustrates that the positive or negative valence of the coefficients for all independent variables except one (the non-significant coefficient for change in percent population in poverty) remained constant across these two models (and also across all tested versions of the multi-domain model), suggesting stability of our findings across models with various non-significant parameters (and other parameters not contributing to overall model fit) removed.


As previously reported by Tempalski and colleagues [1,2,3,4], treatment coverage for PWID in large US metropolitan areas is far below international standards. Some European Union countries, for example, maintain coverage levels of 65% or higher [72]. Although the overall level of treatment coverage for PWID in our 90 MSAs was similarly low at the beginning (6.7%), and end (6.4%) of the study period, it changed significantly during this time period, increasing on average until 2000, then decreasing to its original level.

Contrary to our hypothesis based on the TCA, epidemiologic need, as measured both by the prevalence of AIDS cases per 10,000 population and by HIV prevalence among PWID, seems not to be associated with change in treatment coverage for PWID over time. The present study does find that baseline need as measured by 1993 HIV prevalence rate predicted both level of treatment coverage and variation in trajectories of treatment coverage. These findings, that epidemiologic changes in HIV after 1993 did not correspond with related changes in treatment coverage therefore suggest that service systems are not adequately or efficiently attending to changes in need. They could suggest that need is not, in fact, engendering a direct service provision response at all, which would be consistent with the findings of our previous studies [4,5,6,7].

Such results have important public health implications given the current opiate overdose epidemic, and may indicate that the US should find new ways to allocate resources to drug treatment programs in order to allow responses to a changing need environment. Clearly, drug treatment programs are the basic tools with which public health agencies try to influence and reduce the harms associated with substance misuse. They are thus central to our efforts in helping those vulnerable to substance misuse lead healthy lives. Yet, the US today still struggles with implementing some of the most basic of services for those in need. Public policies need to address the broad individual, environmental, and societal factors that influence substance misuse and its consequences. For example, local politicians have power to shape the nature of care and funding for publicly-funded treatment programs. State licensing and financing policies can provide incentives to programs to offer the full continuum of care (i.e., residential, outpatient, continuing care, and recovery supports) including behavioral treatments and therapy for mental issues and opiate medications, such as buprenorphine. Service providers, harm reduction advocates and researchers can also work toward developing local levels of social and policy support for expanding continuum of care programs in areas experiencing program NIMBYism [39, 40, 43, 46].

As hypothesized, results imply that resource availability does seem to shape service provision of drug treatment. Our model indicates that general resources (both poverty and unemployment rate) are important economic indicators that affect drug treatment coverage (either directly or indirectly). Here, consistent with hypotheses, higher poverty at baseline was found to predict lower levels of treatment coverage.

However, higher baseline rates of unemployment and increases in unemployment over time were associated with higher levels of treatment coverage, and change in unemployment also significantly predicted variation in treatment coverage trajectories over time. This relationship between unemployment and treatment coverage would therefore appear to be driven by processes which are not accounted for by the TCA. One plausible explanation for this positive relationship between unemployment and treatment coverage may be that higher levels of unemployment result in more people entering treatment due to job placement programs that require treatment, or simply due to the availability of people who are not working to participate in treatment [73,74,75,76,77]. For example, research by Popovici and French (2013) [78] and Henkel (2011) [77] have found that both drinking and smoking patterns increase when the economy declines and unemployment rate increases. Both research suggests that the need for treatment services appear to be procyclical with economic turndown. As such, our finding suggests the need for more research and understanding on the effect of changes in unemployment and drug treatment coverage over time.

Additionally, we hypothesized that specialized resources might suggest a concentration of local assets that directly supports a narrow range of programs and/or support a movement for a program or expansion of a program. Previous research measured specialized resources useful for the community action being studied. Such actions have included SPP presence and HIV testing and treatment among PWID [4,5,6,7, 38,39,40], as well as a broad range of programs such as providing condom distribution interventions or programs [20, 21, 36], or specialized programs for treatment of people who are mentally ill [40, 41]. In the present study we found that, consistent with our hypotheses, having a larger community, public health, and social work workforce in MSAs predicted higher treatment coverage levels.

As such these findings strengthens the argument that specialized resources in support of treatment provision are needed especially in communities hit hard by the opioid epidemic. Assets that directly supports a narrow range of programs (i.e., opioid overdose prevention, naloxone distribution) and implementation strategies (i.e., coordinated multi-system & multi-sector public health response driven by community engagement) can take into account special needs and resources in local communities. More specifically, service organizations should articulate the aims of the proposed service in terms that fit in with the local community’s epidemiological needs. Thus, increasing specialized resources for drug treatment may be an issue of pushing for more local government funding for treatment services and developing local coalitions such to do so.

Variables within the institutional opposition domain, such as drug arrests and correction expenditures (i.e., organized and potential opposition to programs), which were hypothesized to predict declines in or lower levels of services, were not found to be associated with treatment coverage. However, evidence from our model suggests that racial structures within this domain (ratio of Black to White poverty) do significantly predict treatment coverage for PWID. Here, we found that in MSAs where Black to White poverty disparity has increased over time, treatment coverage on average is higher than it is in MSAs which experienced a decrease in poverty disparity, with change in disparity not only predicting level of treatment coverage, but also change in treatment coverage over time.

One possible explanation for this relationship could be that, as economic disparities between racial groups increase, and relative disadvantage increases among already disadvantaged groups, drug use may also increase, resulting in a higher number of individuals entering drug treatment [79, 80]. This finding may therefore have important implications both for future directions in epidemiological research which aims to understand factors predicting need, as well as for service systems which strive to meet changing demands in need [81].

Historically, these data report on outcomes from the early 1990s until 2007. This was a time of high need for effective programs related to hard drug use and injection. HIV prevalence and mortality among PWUD in 1993 was at an all time high, and remained so throughout most of that decade. Overdose deaths were continuing in an exponential growth curve dating back at least until 1980 [70, 71, 82]. By standard economic indicators, this was a period of relative prosperity; the Great Recession began only near the end of 2007. Politically, this was a period of War on Drugs but also a period when drug treatment was widely supported as an HIV prevention strategy and as a way to prevent overdoses. Further, when SSPs were proposed and/or established in various cities of the US, opponents often rallied around the watchword that treatment was what was needed and that syringe exchange was a diversion of resources from treatment [43, 44, 47, 48, 82,83,84] Even still, SSPs are banned in 15 states (see

On the one hand, despite the claims of critics of SSPs (i.e, public injecting, inappropriate disposal of used syringes, mortality of drug use and associated illegal activity) [43, 84,85,86,87,88], the presence of a syringe exchange was not associated with decreases in treatment coverage. Need had only a weak relationship with treatment coverage, and changes in need were not associated with changes in coverage. Resource availability indicators (baseline poverty rate and social work workforce) were associated with treatment coverage at baseline in directions suggesting that resource lack retards treatment coverage—yet, in spite of this, at a time of relative economic prosperity, overall treatment coverage in the US did not increase.

The US is today faced with a severe crisis of overdose mortality based primarily on opioid use. Opioid use, unlike most other forms of drug use, has reasonably effective forms of medically assisted treatment available. Federal efforts like the HEALing Community initiative (see have been funded to conduct research into how to provide more treatment for drug users—a critical issue well-deserving of research. Nonetheless, it should be noted that HEALing Community is a research project, and is only funded and design to cover approximately 2 % of US counties.

It is unclear whether funders and policy makers will respond to the results of HEALing Communities with the needed expansion and continuum of care of drug treatment and other services for people who use drugs. Our results in this paper, and in previous research by this team [1,2,3,4,5,6,7, 38, 43, 44] suggest that efforts to increase treatment to have public health scale impacts on overdose deaths face severe obstacles.


Certain data limitations must be taken into account when interpreting the findings from the present study. As discussed in the Methods section, the limitations of both the TEDS and the UFDS/N SSATS data sources used to calculate our estimates of treatment coverage may have resulted in some bias in our estimation of the outcome variable. In addition, change in our estimated numbers of PWID in treatment in an MSA might in part result from measurement error, specifically from change in which and how many treatment facilities in an MSA respond to SAMSHA surveys. The survey response rate increased from 87% in 1995 to 94.5% in 2007, producing a 7.5% increase in reported US client totals from 1995 to 2007. SAMSHA attempts to obtain responses from all known treatment facilities, but the survey is voluntary and no adjustments for facility non- response are made. As a result, the estimated changes in treatment coverage may partially reflect changes in SAMSHA survey methodology over time. UFDS/N-SSATS data were unavailable for 1994, 1999, and 2001. Consequently, our drug treatment estimates were only created for years where data were available. Thus, our final coverage estimates only include data for 1993, 1995, 1996–1998, 2000, and 2002–2007. Finally, PWID estimates beyond 2007 were not available for our coverage estimates. As such, our data and analyses do not extend beyond 2007 due to the lack of PWID population denominators. Although exact estimates are hard to come by, estimating the contribution of drug treatment availability in preventing opioid-related morbidity and mortality among PWID is key for public health and common sense health policy in reducing harms.


Programs are the basic tools with which public health agencies try to influence HIV transmission and disease progression. Despite this, however, relatively little is known about what determines their presence and reach. Much health policy discourse assumes that need for a program is associated with program presence or magnitude— however our previous studies of the determinants of drug treatment coverage and syringe exchange presence have found that need is not a predictor of these programs for PWID [4,5,6,7, 38]. The present study, similarly, has found that several key indicators of need were not related to 1993–2007 drug treatment coverage. While one 1993 indicator of need (HIV prevalence) was found by the present study to be related to differences in coverage, lack of a relationship between treatment coverage and change in need over time may still point to a lack of synergy between local need and service provision.

The present study also has found that several indicators of resource availability and economic conditions at the MSA level were related to variation in levels of treatment coverage from 1993 to 2007. These associations are especially important findings of the present study given the flux in economic conditions around the US and increased disparity among racial/ethnic groups. Specifically, a time of increased economic difficulty and instability since the early 2000’s has led to higher rates of unemployment, poverty, and evictions and foreclosures. This may have resulted in an increase in the level of economic disadvantage within neighborhoods.

Economic disadvantage has been associated with a variety of social problems, including income inequality, housing instability and crime, and to increased substance use prevalence rates [79,80,81, 89]. Economic changes appear to be important factors in predicting changes in treatment services [76, 90,91,92,93,94]. Given these empirical and theoretical links between treatment services and economic conditions, future research should consider the implications of the past and present financial conditions for treatment coverage.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.


  1. During our study period, SAMSHA eliminated questions from UFDS about the number of PWIDs in a treatment program. To correct for this, we multiplied the proportion of drug users who inject drugs (from TEDS) in each MSA and year by the total number of drug users in treatment as reported by UFDS/N-

    SSATS. Second, these data sets differ in what they count: TEDS counts each admission in a given year, so an individual who enters drug treatment twice or more in a year is counted as two or more independent cases. In contrast, UFDS/N-SSATS is a point-prevalence survey of treatment facilities and clients on a specific date. Thus, if PWIDs differ from non-PWIDs in the ratio of admissions to those remaining in treatment, our estimates will be biased. In addition, these data do not capture medication-assisted treatment (MAT) operating out of private medical offices. Lastly, we need to acknowledge much of any system response to the opioid epidemic would have taken place after 2007, and so our data would miss those more recent changes.



AIDS Coalition to Unleash Power


Acquired immune deficiency syndrome




Antiretroviral therapy




Centers for Disease Control


Hepatitis C virus


Human immunodeficiency virus






US metropolitan statistical areas


Men who have sex with Men


North Carolina


New Jersey


National Survey of Substance Abuse Treatment Services


New York






People who inject drugs


People who use drugs


Substance Abuse and Mental Health Service Administration


Syringe Service Programs


Theory of Community Action


Treatment Episode Data Set


Uniform Facility Data Set


United States




  1. Tempalski B, Cleland CM, Williams LD, Cooper HLF, Friedman SR. Change and variability in drug treatment coverage among people who inject drugs in 90 large metropolitan areas in the USA, 1993-2007. Subst Abuse Treat Prev Policy. 2018;13:1.

    Article  Google Scholar 

  2. Tempalski B, Cleland CM, Pouget ER, Chatterjee S, Friedman SR. Persistence of low drug treatment coverage for injection drug users in large US metropolitan areas. Subst Abuse Trea Prev Policy. 2010;5:3.

    Article  Google Scholar 

  3. Friedman SR, Tempalski B, Cooper H, Perlis T, Keem M, et al. Estimating numbers of injecting drug users in metropolitan areas for structural analyses of community vulnerability and for assessing relative degrees of service provision for injecting drug users. J Urban Health. 2004;81:377–400.

    PubMed  PubMed Central  Article  Google Scholar 

  4. Friedman SR, Tempalski B, Brady J, Friedman JJ, Cooper H, Flom PL, McGrath MM, Gostnell K, Des Jarlais DC. Predictors of the degree of drug treatment coverage for injection drug users in 94 metropolitan areas in the United States. Int J Drug Pol. 2007;18:475–85.

    Article  Google Scholar 

  5. Tempalski B, Friedman SR, Des Jarlais DC, McKnight C, Keem M, Friedman R. What predicts which metropolitan areas in the USA have syringe exchanges? Int J Drug Policy. 2003;14:417–24.

    Article  Google Scholar 

  6. Tempalski B, Flom PL, Friedman SR, Des Jarlais DC, Friedman JJ, McKnight C, Friedman R. Social and political factors predicting the presence of syringe exchange programs in 96 metropolitan areas in the United States. Am J Public Health. 2007;97:3.

    Article  Google Scholar 

  7. Tempalski B, Cooper H, Friedman SR, Des Jarlais DC, Brady J. Correlates of syringe coverage for heroin injection in 35 large metropolitan areas in the US in which heroin is the dominant injected drug. Int J Drug Policy. 2008;19:47–58.

    Article  Google Scholar 

  8. Sheppard VB, Hannon M, Barbosa CE, Pestian JP. From research to community action: an assessment of child and adolescent hospitalizations. Family Comm Health. 1999;22:1.

    Article  Google Scholar 

  9. Poland B, Boutilier M, Tobin S, Badgley R. The policy context for community development practice in public health: a Canadian case study. J Public Health Policy. 2000;21:1.

    Article  Google Scholar 

  10. Miller BA. Geography and social movements: comparing antinuclear activism in the Boston area. Minneapolis: University of Minnesota Press; 2000.

    Google Scholar 

  11. Nathanson C. Disease prevention as social change: toward a theory of public health. Popul Dev Rev. 1996;22:609.

    Article  Google Scholar 

  12. McCarthy JD, Wolfson M, Baker DP, Mosakowski E. The founding of social movement organizations: local citizens’ groups opposing drunken driving. In G. R. Carroll (Ed.), ecological models of organizations, pp. 71–82. Cambridge: Ballinger; 1988.

    Google Scholar 

  13. Brown P, Zavestoski S, McComiick S, Mayer B, Morello-Frosch R, Gasior R. Embodied health movements: new approaches to social movements in health. Sociol Health Illn. 2004;26:50–80.

    PubMed  Article  Google Scholar 

  14. DeMarco R, Johnsen C. Taking action in communities: women living with HIV/AIDS lead the way. J Comm Health Nursing. 2003;20:1.

    Article  Google Scholar 

  15. Chiotti QP, Joseph AE. Casey house: interpreting the location of a Toronto AIDS hospice. Soc Sci Med. 1995;41:131–40.

    CAS  PubMed  Article  Google Scholar 

  16. Barnett R, Barnett P. If you want to sit on your butts you'll get nothing!: community activism in response to threats of rural hospital closure in southern New Zealand. Health Place. 2003;9:2.

    Google Scholar 

  17. Bruggmann P, Grebely J. Prevention, treatment and care of hepatitis C virus infection among people who inject drugs. Int J Drug Policy. 2005;26:S22–6.

    Article  Google Scholar 

  18. Wilson DP, Donald B, Shattock AJ, Wilson D, Frazer-Hunt N. The cost-effectiveness of harm reduction. Int J Drug Policy. 2015;26:S5–11.

    PubMed  Article  Google Scholar 

  19. Friedman SR, West BS, Tempalski B, Morton CM, Cleland CM, Des Jarlais DC, Hall HI, Cooper HL. Do metropolitan HIV epidemic histories and programs for people who inject drugs and men who have sex with men predict AIDS incidence and mortality among heterosexuals? Ann Epidemiol. 2014;24:304.

    PubMed  PubMed Central  Article  Google Scholar 

  20. Abdul-Quader AS, Feelemyer J, Modi S, Stein ES, Briceno A, Semaan S, et al. Effectiveness of structural-level needle/syringe programs to reduce HCV and HIV infection among people who inject drugs: a systematic review. AIDS Behav. 2013;17:9.

    Article  Google Scholar 

  21. Arfken CL, Johanson C, di Menza S, Schuster CR. Expanding treatment capacity for opioid dependence with office-based treatment with buprenorphine: national surveys of physicians. J Subst Abus Treat. 2010;39:96–104.

    Article  Google Scholar 

  22. Friedman SR, West BS, Pouget ER, Hall HI, Cantrell J, Tempalski B, et al. Metropolitan social environments and pre-HAART/HAART era changes in mortality rates (per 10,000 adult residents) among injection drug users living with AIDS. PLoS One. 2013;8:2.

    Article  Google Scholar 

  23. Kresina TF, Lubran R. Improving public health through access to and utilization of medication assisted treatment. Int J Environ Res Public Health. 2011;8:4102.

    PubMed  PubMed Central  Article  Google Scholar 

  24. Metzger DS, Woody GE, O’Brian CP. Drug treatment as HIV prevention: A research update. J Acquir Immune Defic Syndr. 2010;55:S32.

    PubMed  PubMed Central  Article  Google Scholar 

  25. Smith N. Uneven development: nature, capital and the production of space. New York: Basil Blackwell; 1991.

    Google Scholar 

  26. Mitchell D. The right to the City: social justice and the fight for public space in America. New York: Guilford; 2003.

    Google Scholar 

  27. Blomley NK. Unsettling the city: urban land and the politics of property. New York: Routledge; 2004.

    Book  Google Scholar 

  28. Stoller NE. Lessons from the damned: queers, whores, and junkies respond to AIDS. New York, NY: Routledge; 1998.

    Google Scholar 

  29. Risse GB. Epidemics and history: ecological perspectives and social responses. In: Fee E, Fox DM, editors. AIDS: the burdens of history. Berkeley: University of California Press; 1988. p. 33–66.

    Google Scholar 

  30. Keith M, Pile S, editors. Geographies of resistance. London: Routledge; 1997.

    Google Scholar 

  31. Rogers E. Diffusion of innovations. New York: The Free Press; 2003.

    Google Scholar 

  32. Barker K. Diffusion of innovations: a world tour. J Commun. 2004;9:131–7.

    Google Scholar 

  33. Bero LA, Grilli R, Grimshaw JM, Harvey E, Oxman AD, Thomson MA. Closing the gap between research and practice: an overview of systematic reviews of interventions to promote the implementation of research findings. Br Med J. 1998;317:465–8.

    CAS  Article  Google Scholar 

  34. Des Jarlais DC, Sloboda Z, Friedman SR, Tempalski B, McKnight C, Braine N. Comparing the diffusion of D.A.R.E and syringe exchange programs. Am J Public Health. 2006;96(8):1354–8.

    PubMed  PubMed Central  Article  Google Scholar 

  35. Rosser BRS, Horvath KJ. Predictors of success in HIV prevention in rural America: a state level structural factor analysis of HIV prevention targeting men who have sex with men. AIDS & Behavior. 2008;12:159–68.

    Article  Google Scholar 

  36. Mahnaz RC, Crepaz N, Guenther-Gray C, Henny K, Liau A, Willis LA, Lyles CM. Efficacy of Structural-Level Condom Distribution Interventions: A Meta-Analysis of U.S. and International Studies, 1998–2007. AIDS Behav. 2011;15:7.

    Google Scholar 

  37. Downing M, Riess TH, Vernon K, Mulia N, Hollinquest M, McKnight C, et al. It takes a community: successful implementation of syringe exchange programs. AIDS Educ Prev. 2005;17:1.

    Article  Google Scholar 

  38. Tempalski B, Cooper HLF, Kelley ME, Linton SL, Semaan S, DiNenno E, Wejnert C, Broz D, Paz-Bailey G, National HIV Behavioral Surveillance Study Group. Identifying which place characteristics are associated with the odds of recent HIV testing in a large sample of people who inject drugs in 19 US metropolitan areas. AIDS Behav. 2018;23(2):318–35.

    Article  Google Scholar 

  39. Bluthenthal RN. Syringe exchange as a social movement: a case study of harm reduction in Oakland, California. Subst Use Misuse. 1998;33:5.

    Article  Google Scholar 

  40. Shaw SJ. Public citizens, marginalized communities: the struggle for syringe exchange in Springfield, Massachusetts. Med Anthropol. 2006;25:31.

    PubMed  Article  Google Scholar 

  41. Drainoni ML, Lee-Hood E, Tobia C, Bachman S, Andrew J, Maisels L. Cross-disability experiences of barriers of health-care access. Journal of Disability Policy Studies. 2006;17:2.

    Article  Google Scholar 

  42. Beatty PW, Hagglund KJ, Neri MT, Dhont KR, Clark MJ, Hilton S. Access to health care services among people with chronic or disabling conditions: patterns and predictors. Arch Phys Med Rehabil. 2003;84:1417.

    PubMed  Article  Google Scholar 

  43. Tempalski B, Friedman R, Keem M, Cooper H, Friedman SR. NIMBY localism and national inequitable exclusion alliances: the case of syringe exchange programs in the United States. Geoforum. 2007;38:1250.

    PubMed  PubMed Central  Article  Google Scholar 

  44. Tempalski B. Placing the dynamics of syringe exchange programs in the United States. Health & Place. 2007;13:2.

    Article  Google Scholar 

  45. Friedman SR, Pouget ER, Chatterjee S, Cleland CM, Tempalski B, Brady JE, Cooper HLF. Do drug arrests deter injection drug use? Am J Public Health. 2011;101:344.

    PubMed  PubMed Central  Article  Google Scholar 

  46. Cooper HLF, Des Jarlais DC, Tempalski B, Bossak B, Ross Z, Friedman SR. Drug-related arrest rates and spatial access to syringe exchange programs in New York City health districts: combined effects on the risk of injection-related infections among injectors. Health Place. 2012;18:2.

    Google Scholar 

  47. Des Jarlais DC, Paone D, Friedman SR, Peyser N, Newman RG. Public health then and now: regulating controversial programs for unpopular people: methadone maintenance and syringe-exchange programs. Am J Public Health. 1995;85:11.

    Google Scholar 

  48. Des Jarlais DC, Pinkerton S, Hagan H, Guardino V, Feelemyer J, Cooper HL, et al. 30 years on selected issues in the prevention of HIV among persons who inject drugs. Adv Prev Med. 2013.

    Article  Google Scholar 

  49. Dear M, Wilton R, Gaber SL, Takahashi LM. Seeing people differently: the sociospatial construction of disability. Environ Plann. 1997;15:455.

    Article  Google Scholar 

  50. Law R. Communities, citizens, and the perceived importance of AIDS- related services in West Hollywood, California. Health Place. 2003;9:7.

    PubMed  Article  Google Scholar 

  51. Law RM, Takahashi LM. HIV, AIDS and human services: exploring public attitudes in West Hollywood, California. Health Soc Care Comm. 2000;8(2):90–108.

    Article  Google Scholar 

  52. Lyon-Callo V. Making sense of NIMBY: poverty, power and community opposition to homeless shelters. City Community. 2001;13:183–209.

    Google Scholar 

  53. Petchey R, Williams J, Farnsworth B, Starkey K. A tale of two (low prevalence) cities: social movement organizations and the local policy response to HIV/AIDS. Soc Sci Med. 1998;47:9.

    Article  Google Scholar 

  54. Takahashi LM. Homelessness, AIDS, and stigmatization. Oxford: Clarendon Press; 1998.

    Google Scholar 

  55. Paone D, Clark J, Shi Q, Purchase D, Des Jarlais DC. Syringe exchange in the United States, 1996: A national profile. Am J Public Health. 1999;89:43.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  56. Hoffman B. Health care reform and social movements in the United States. Am J Public Health. 2003;93:1.

    Google Scholar 

  57. Office of Management and Budget. Standards for defining metropolitan and micropolitan statistical areas. In: Federal Register, vol. 65; 2000. p. 8228–82238.

    Google Scholar 

  58. US Department of Health and Human Services. Substance Abuse and Mental Health Services Administration, Office of Applied Studies. Treatment episode data set, 1992–2008 (concatenated). Ann Arbor: Inter-university Consortium for Political and Social Research; 2010.

    Google Scholar 

  59. US Department of Health and Human Services, Substance Abuse and Mental Health Services Administration, Office of Applied Studies. Uniform Facility Data Set (UFDS): 1993–1995–Data on Substance Abuse Treatment Facilities; 1998.

  60. US Department of Health and Human Services. Substance Abuse and Mental Health Services Administration, Office of Applied Studies. Uniform Facility Data Set (UFDS): 1997–1998–Data on substance abuse treatment facilities. [Computer file]. 2nd ICPSR version. Ann Arbor: Inter-university Consortium for Political and Social Research; 2004.

    Google Scholar 

  61. US Department of Health and Human Services. Substance Abuse and Mental Health Services Administration, Office of Applied Studies. National Survey of Substance Abuse Treatment Services (N-SSATS), 2000–2008. [Computer file]. 3rd ICPSR version. Ann Arbor: Inter-university Consortium for Political and Social Research; 2010.

    Google Scholar 

  62. Tempalski B, Pouget ER, Cleland CM, Brady JE, Cooper HLF, Hall HI, Lansky A, West BS, Friedman SR. Trends in the population prevalence of people who inject drugs in US metropolitan areas 1992-2007. PLoS One. 2013;8(6):e64789.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  63. Centers for Disease Control and Prevention. HIV Counseling and Testing in Publicly Funded Sites. Atlanta: Centers for Disease Control and Prevention; 2010.

    Google Scholar 

  64. Holmberg S. The estimated prevalence and incidence of HIV in 96 large US metropolitan areas. Am J Public Health. 1996;86:642.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  65. Singer JD, Willet JB. (Eds). Applied longitudinal data analysis. Modeling change and event occurrence. New York: Oxford University Press; 2003.

    Book  Google Scholar 

  66. SAS Institute Inc. 2002–SAS 9.2. Cary: SAS Institute Inc; 2009.

  67. Vrieze SI. Model selection and psychological theory: a discussion of the differences between the Akaike information criterion (AIC) and the Bayesian information criterion (BIC). Psychol Methods. 2012;17:228.

    PubMed  PubMed Central  Article  Google Scholar 

  68. Kuha J. AIC and BIC: comparisons of assumptions and performance. Sociol Methods Res. 2014;33:188.

    Article  Google Scholar 

  69. Ferguson CJ. An effect size primer: a guide for clinicians and researchers. Prof Psychol Res Pract. 1990;40:5.

    Google Scholar 

  70. Unick GJ, Rosenblum D, Mars S, Ciccarone D. Intertwined epidemics: National Demographic Trends in hospitalizations for heroin- and opioid-related overdoses, 1993–2009. PLoS One. 2013;8(2):e54496.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  71. Ciccarone D, Unick GJ, Kraus A. Impact of south American heroin on the US heroin market 1993–2004. Int J Drug Policy. 2009;20:392–401.

    PubMed  PubMed Central  Article  Google Scholar 

  72. European Monitoring Centre for Drugs and Drug Addiction Health and social responses to drug problems. a European guide. Luxembourg: Publications Office of the European Union; 2017.

    Google Scholar 

  73. Kemp PA, Neale J. Employability and problem drug users. Crit Soc Policy. 2005;25:28.

    Article  Google Scholar 

  74. Magura S. The role of work in substance dependency treatment: a preliminary overview. Subst Use Misuse. 2003;38:1865.

    PubMed  Article  Google Scholar 

  75. Magura S, Staines GL, Blankertz L, Madison E. The effectiveness of vocational services for substance users in treatment. Subst Use Misuse. 2004;39:2165.

    PubMed  Article  Google Scholar 

  76. McIntosh J, Bloor M, Robertson M. Drug treatment and the achievement of paid employment. Addict Res Theory. 2008;16:37.

    Article  Google Scholar 

  77. Henkel D. Unemployment and substance use: a review of the literature 1990-2010. Current Drug Abuse Reviews. 2011;4:1.

    Article  Google Scholar 

  78. Popovici I, French MT. Does Unemployment Lead to Greater Alcohol Consumption? Ind Relat (Berkeley). 2013;52(2):444–66. Epub 2013 Mar 18. PMID: 23543880; PMCID: PMC3609661.

    Article  Google Scholar 

  79. Winstanley EL, Steinwachs DM, Ensminger ME, Latkin CA, Stitzer ML, Olsen Y. The association of self-reported neighborhood disorganization and social capital with adolescent alcohol and drug use, dependence, and access to treatment. Drug Alcohol Depend. 2008;92:1–3.

    Article  Google Scholar 

  80. Freisthler B, LaScala EA, Gruenewald PJ, Treno AJ. An examination of drug activity: effects of neighborhood social organization on the development of drug distribution systems. Subst Use Misuse. 2005;40:5.

    Article  Google Scholar 

  81. Niewczyk PM, Lwebuga-Mukasa JS. Is poverty the Main factor contributing to health care disparities? an investigation of individual level factors contributing to health care disparities. J Equity Health. 2008;1:1.

    Google Scholar 

  82. Jalal H, Buchanich JM, Roberts MS, Balmert LC, Zhang K, Burke DS. Changing dynamics of the drug overdose epidemic in the United States from 1979 through 2016. Science. 2018.

    Article  CAS  PubMed  Google Scholar 

  83. Bastos FI, Strathdee SA. Evaluating effectiveness of syringe exchange programmes: current issues and future prospects. Soc Sci Med. 2000;51:12.

    Article  Google Scholar 

  84. Des Jarlais DC. Research, politics, and needle exchange. Am J Pub Health. 2000;90(9):1392–4.

    CAS  Article  Google Scholar 

  85. Fernandes RM, Cary M, Duarte G, et al. Effectiveness of needle and syringe programs in people who inject drugs – an overview of systematic reviews. BMC Public Health. 2017;17(1):309.

    PubMed  PubMed Central  Article  Google Scholar 

  86. Backes G, Rose VJ. Primary and secondary analysis of local elected officials’ decision to support or oppose pharmacy sale of syringes in California. J Urban Health. 2010;87(4):553–60.

    PubMed  PubMed Central  Article  Google Scholar 

  87. Moss AR. Epidemiology and the politics of needle exchange. Am J Public Health. 2000;90(9):1385–7.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  88. March JC, Oviedo-Joekes E, Romero M. Drugs and social exclusion in ten European cities. Eur Addict Res. 2006;12:33.

    PubMed  Article  Google Scholar 

  89. Link BG, Phelan JC. Fundamental sources of health inequalities. In: Policy Challenges in Modern Health Care. Piscataway: Rutgers University Press; 2005.

    Google Scholar 

  90. Costa Storti C, De Grauwe P, Sabbadash A, Montanari L. Unemployment and drug treatment. Int J Drug Policy. 2011;22:5.

    Article  Google Scholar 

  91. Fergusson DM, Horwood LJ, Woodward LJ. Unemployment and psychosocial adjustment in young adults: causation or selection? Soc Sci Med. 2001;53:305.

    CAS  PubMed  Article  Google Scholar 

  92. McCoy C, Comerford M, Metsch L. Employment among chronic drug users at baseline and 6-month follow-up. Subst Use Misuse. 2007;42:7.

    Google Scholar 

  93. Ruhm CJ. Good times make you sick. J Health Econ. 2003;24:637.

    Article  Google Scholar 

  94. Storti CC, De Grauwe P, Reuter P. Economic recession, drug use and public health. Int J Drug Policy. 2011;22:5.

    Article  Google Scholar 

Download references


The authors thank Dr. Charles M. Cleland from the New York University, Population Health New York, NY, and Chunki Fong from NDRI, Inc., New York, NY, USA for providing statistical guidance. We kindly acknowledge from the Department of Psychiatry at the University of California San Francisco, Professor’s Elinore F. McCance-Katz and James L. Sorensen; Professor Joseph R. Guydish from Philip R. Lee Institute for Health Policy Studies, University of California, San Francisco; Dr’s Darren Urada and Richard A. Rawson at the UCLA Integrated Substance Abuse Programs in Los Angeles, CA; and Dr. Alex H. Kral Director, Urban Health Program at RTI International, San Francisco, CA for their knowledge and reflection of drug treatment coverage for PWID in the State of California. The authors also thank Mr. Lucas Wiessing from the European Monitoring Centre for Drugs and Drug Addiction in Lisbon, Portugal for his thoughtful discussion and knowledge regarding drug treatment coverage in the European Union and elsewhere. We thank from the Substance Abuse and Mental Health Services Administration, Dr. Cathie E. Alderks for her assistance on SAMHSA drug treatment data issues.


Research reported in this publication was supported by the National Institute of Drug Abuse of the National Institutes of Health under Award Number R01 DA037568; Metropolitan Trajectories of HIV Epidemics and Responses in US Key Populations. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. As such, the content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.

Author information

Authors and Affiliations



Each author has contributed to the conception and design of the work, the acquisition of data or the analysis of the data in a manner substantial enough to take public responsibility for it. In addition, each author believes that the paper represents valid work and has reviewed the final version of the manuscript and approves consent for publication. Study concept and design: BT, SRF. Acquisition of data: BT, BSW. Analysis and interpretation of data: BT, LDW, SRF. Drafting of the manuscript: BT, LDW, BSW, HLFC, SRF, SB, UI. Obtaining Funding: HLFC, SRF. Study supervision: SRF. BT, BSW, LDW had full access to all of the data and takes responsibility for the integrity of the data and the accuracy of the data analysis.

Corresponding author

Correspondence to Barbara Tempalski.

Ethics declarations

Ethics approval and consent to participate

Not applicable

Consent for publication

Each author has contributed to the conception and design of the work, the acquisition of data or the analysis of the data in a manner substantial enough to take public responsibility for it. In addition, each author believes that the paper represents valid work and has reviewed the final version of the manuscript and approves consent for publication.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1: Table S1.

Bivariate correlation matrix among all variables. Table S2. Estimated drug treatment coverage rates. Table S3. Results from Full Multi-Domain Mixed-Effects Model by Domain.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Tempalski, B., Williams, L.D., West, B.S. et al. Predictors of historical change in drug treatment coverage among people who inject drugs in 90 large metropolitan areas in the USA, 1993–2007. Subst Abuse Treat Prev Policy 15, 3 (2020).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Injection drug use
  • Predictors
  • Drug treatment coverage
  • Longitudinal
  • Mixed-effects multivariate models
  • Metropolitan areas
  • Drug policy
  • Theory of community action