Abstract
Studies document inequitable tobacco retailer density by neighborhood sociodemographics, but these findings may not be robust to different density measures. Policies to reduce density may be less equitable depending on how the presence of store types differs by neighborhood characteristics. We built a 2018 list of probable tobacco retailers in the United States and calculated four measures of density for all census tracts (
Tobacco retailer availability is a geographic construct, representing the supply of tobacco retailers in a geographic area, and is often quantified through the use of tobacco retailer density measures (Kong et al., 2020; Valiente et al., 2020). In places with higher tobacco retailer density, individuals have greater smoking intentions, are more likely to start or continue smoking, smoke more cigarettes, and are less likely to quit smoking (Finan et al., 2019; Marsh et al., 2020; Valiente et al., 2020).
Several U.S.-based (Hyland et al., 2003; Peterson et al., 2005; Schneider et al., 2005; Siahpush et al., 2010) and international (Chaiton et al., 2013; Kite et al., 2014; Shortt et al., 2015) studies have documented higher tobacco retailer density in neighborhoods with a greater proportion of some shared demographic and social characteristics, potentially putting residents of these neighborhoods at a higher risk of smoking. However, few U.S. studies have assessed neighborhood inequities in tobacco retailer density at the national level. Using census tract sociodemographic characteristics from nearly all census tracts in 2000, Rodriguez et al. (2013, 2014) found that the proportion of Hispanic residents, Black residents, families living in poverty, and urbanicity were each uniquely and positively associated with the number of tobacco retailers per 1,000 people. However, no study of this scope has been completed with more recent data.
The first purpose of this study is to assess and compare associations of four common measures of tobacco retailer density with recent (2018) tract-level sociodemographic characteristics in the United States. There is substantial variation in how researchers and policymakers measure retailer density (Valiente et al., 2020) with little discussion on why specific measures were chosen. Some common density measures include the number of tobacco retailers per population, per land area, or per kilometers of roadway. There has also been little comparison across density measures both within and between studies (Mayers et al., 2012). As a result, conclusions about area-level sociodemographic inequities may not be robust to different measures, potentially under- or overestimating inequities. Additionally, understanding whether findings are similar across measures may be useful for past and future study comparison purposes.
Studies in select U.S. cities and states also suggest that the impact of policies that prohibit pharmacies from selling tobacco products may not be equitable (Craigmile et al., 2020; Giovenco et al., 2019; Hall et al., 2019; Kulbicki & Leslie, 2015). Additionally, some places have implemented policies that only allow specialty tobacco shops to sell tobacco products or certain types of tobacco products (CounterTobacco.org, 2020). While a study of U.S. vape shop density documented some area-level racial, ethnic, and socioeconomic inequities (Dai et al., 2017), it is not known if there are inequities in the availability of all specialty tobacco shops (e.g., vape shops
Method
Neighborhood Sociodemographic Characteristics
We conceptualized race and ethnicity as a social construct, resulting from racialization and discriminatory systems that create and sustain group-based hierarchies that advantage and disadvantage certain groups (Inwood & Yarbrough, 2010; Krieger et al., 2015; Pratto et al., 2006). Both individuals and places are racialized, resulting in spatially patterned racial inequities in the distribution of resources, health behaviors, and health outcomes (Inwood & Yarbrough, 2010; Lipsitz, 2007; Neely & Samura, 2011). To measure this construct, we used U.S. Census Bureau 2014–2018 American Community Survey population estimates, which measure the census tract-level percentage of the population in different racial and ethnic groups (non-Hispanic Black or African American [Black], Hispanic or Latino), based on self-report survey categories.
While most studies document inequities in retailer density by median household income or federal poverty level (FPL), very few studies have expanded beyond these economic measures. Vacant housing may be an indicator for residential or neighborhood stability: A lack of residential stability may reflect low attachment to a community, potentially leading to fewer opportunities for social cohesion and collective efficacy among its residents to empower a healthy community (Aiyer et al., 2015; Browning & Cagney, 2003; Kawachi & Berkman, 2000; Stockdale et al., 2007), such as organizing to prohibit tobacco retailers. In this study, we included two neighborhood economic variables: percentage of the tract-level population that was living below 150% of the FPL and vacant housing units.
We incorporated a measure of tract-level urbanicity as a control variable, using U.S. Department of Agriculture (2016) Rural-Urban Commuting Area Codes Documentation. Each census tract was categorized as urban, large rural city/town, or small and isolated rural town.
Tobacco Retailer Density
There is no national tobacco retailer licensing system in the United States, so we created a 2018 list of probable retailers that sold tobacco products intended for off-premise consumption, similar to previous studies (Kong et al., 2021; Ribisl et al., 2017; Rodriguez et al., 2013). In short, we used tobacco product sales data from the latest 2017 Economic U.S. Census to identify 11 North American Industry Classification System (NAICS) store type codes (e.g., Convenience Stores, Tobacco Stores, Pharmacies and Drug Stores) that account for approximately 99% of all retail tobacco product sales. Every NAICS code is associated with a subset of Standard Industrial Classification (SIC) codes, which describes the primary business activity of a retailer in more detail. For example, NAICS 453991 (Tobacco Stores) is associated with several SIC codes, including 599302 (Smoke Shops and Supplies) and 599306 (Electronic Cigarettes). Using these NAICS codes and ReferenceUSA, a database of business establishments that contains NAICS and associated SIC codes, retailer name, and geographic indicators for each retailer, we created a 2018 list of probable brick-and-mortar tobacco retailers. Specific retailer subtypes identified through SIC codes that seemed unlikely to sell tobacco products (e.g., craft galleries and dealers, marine services stations) were excluded from the sample. We further refined the list through a text search of store names to exclude those retailers known to not sell tobacco products (e.g., Target, Whole Foods, CVS Pharmacy).
Consistent with common measures in the literature, we operationalized tobacco retailer density in four ways: (1) total count of tobacco retailers, and tobacco retailers per (2) 1,000 people; (3) land area (square mile); and (4) 10 km of roadway. To calculate these measures, we used a spatial join in ArcMap 10.5 to assign each retailer to its respective census tract and then summed the total number of tobacco retailers within each tract. Publicly available roadway data were downloaded from the U.S. Census Bureau.
We calculated both the total number of tobacco-selling pharmacies and the number of tobacco shops in a tract. We then created two binary variables: One variable indicated whether a tract had at least one tobacco-selling pharmacy (=1), while the other variable indicated whether a tract had at least one tobacco shop (=1). As smaller, independent, and compounding pharmacies may be less likely to sell tobacco products (CounterTobacco.org, 2019), we included only those pharmacies that had at least 10 locations nationally. Similar to other work (Lee et al., 2018), we confirmed that the most frequently occurring pharmacy chain brands sold tobacco products (e.g., Walgreens, Duane Reade, Rite Aid, Thrifty White) and excluded those chain pharmacies known not to sell tobacco products (e.g., CVS, Medicine Shoppe). Using data collated by the American Nonsmoker’s Rights Foundation U.S. Tobacco Control Laws Database©, we identified places with policies that prohibit pharmacies from selling tobacco products and then excluded pharmacies that were located in those places. Finally, as current bans on tobacco sales in pharmacies also include stores (e.g., grocery, warehouse) that have a pharmacy counter (Craigmile et al., 2020), we included retailers that were identified as such (e.g., Kroger Pharmacy) and that met the above criteria.
Analytic Sample
In 2018, there were 72,377 census tracts with a population of at least one person and that had Rural-Urban Commuting Area and roadway data for all 50 states and the District of Columbia. As tracts are intended to range from 1,200 to 8,000 people and our per capita retailer density measure is per 1,000 people, we excluded those tracts with fewer than 1,000 people (
Analysis
To investigate associations between the four measures of tobacco retailer density and neighborhood sociodemographic characteristics, we fit unadjusted and adjusted multivariable linear regression models that controlled for the other tract-level variables described previously and area urbanicity. Because the majority of census tracts did not have a tobacco-selling pharmacy (72.3%) or tobacco shop (81.2%), and there was little variability in the count distribution for those that did, we fit unadjusted and adjusted logistic regression models to assess relationships between sociodemographic characteristics and two outcome variables: the presence of at least one (vs. none) tobacco-selling pharmacy and the presence of at least one (vs. none) tobacco shop.
All adjusted models also included a state fixed effect to account for potential state-level differences (e.g., historically tobacco growing state, tobacco retailer licensing laws). We did not find evidence of collinearity in multivariable models (average variance inflation factor was 1.34). To aid in interpretability, each sociodemographic variable was scaled to tens (e.g., 13% = 1.3) so that a 1-unit difference in a sociodemographic variable represents a 10–percentage point difference.
Results
Characteristics of the analytic sample are shown in Table 1. In 2018, there were an estimated 325,884 tobacco retailers. The average number of tobacco retailers in a tract was 4.6. Average retailer density per 1,000 people was 1.11; per square mile was 4.85; and per 10 km of roadway was 1.20. About 28% and 19% of tracts had at least one tobacco-selling pharmacy or tobacco shop, respectively. Comparing density measures, Pearson correlation coefficients were high for retailers per roadway and square mile (
Sociodemographic and Tobacco Retailer Availability Characteristics of Census Tract Neighborhoods, United States, 2018 (
Tobacco Retailer Density
For all measures of retailer density, unadjusted analyses indicated positive and statistically significant associations for tract-level composition of percentage of Black residents and residents living below 150% FPL (Table 2). Percentage of Hispanic or Latino composition was positive and significant for all density measures except per 1,000 people (Β = −0.01,
Unadjusted Analyses Testing Census Tract-Level Associations of Percentage Sociodemographics With Measures of Tobacco Retailer Density, United States, 2018 (
In multivariable adjusted models (Table 3), there were positive and significant associations between the percentage of Black residents and land area and roadway retailer density.
Adjusted Analyses Testing Census Tract-Level Associations of Percentage Sociodemographics With Measures of Tobacco Retailer Density, United States, 2018 (
When comparing two census tracts, one of which had a Black composition that was 10 percentage points higher than the other, we would expect to have 0.18 and 0.04 (
Associations for percentage of Hispanic or Latino residents were positive for all retailer density measures except per 1,000 people, which had a negative association (Β = −0.02,

Expected percentage difference in tobacco retailer density relative to average density by tract-level sociodemographic characteristics, United States, 2018 (
In a post hoc multivariable analysis, we included interaction terms between each sociodemographic variable and the three-level urbanicity variable. We report simple slope estimates for each urbanicity category and further tested whether associations significantly differed (
Pharmacies and Tobacco Shop Presence
Finally, we investigated associations between neighborhood sociodemographics and the presence of a pharmacy and the presence of a tobacco shop in a tract (Table 4). Adjusted results indicated that neighborhoods with a greater percentage of Black residents (adjusted odds ratio [a
Analyses Testing Tract-Level Associations of Percentage Sociodemographics With the Presence (vs. Absence) of a Tobacco-Selling Pharmacy or Tobacco Shop, United States, 2018 (
Discussion
In 2018, tobacco retailer density differed by neighborhood characteristics, but the extent and direction of these associations were sometimes sensitive to the density measure used. Neighborhoods with a greater percentage of residents living below 150% FPL were associated with higher tobacco retailer density in all models. In our results, however, the direction and significance of adjusted associations between retailer density and percentage of non-Hispanic Black residents, Hispanic or Latino residents, and vacant housing units were sensitive to the density measure operationalized.
Several U.S. studies have documented greater retailer density in neighborhoods with a higher percentage of Black residents (Hyland et al., 2003; Rodriguez et al., 2013). Our unadjusted results concur with these findings for all retailer density measures. However, in adjusted models, this association was negative for both total count of retailers and the number of retailers per 1,000 people. Although this finding may suggest that there is less tobacco availability in neighborhoods with a greater percentage of Black residents, adjusted results represent what would be observed in a counterfactual scenario where confounding sociodemographic variables in the model were held constant. Therefore, given that all unadjusted and some adjusted models indicated an inequity by percentage Black, greater retailer availability in these neighborhoods is still a concern in the real world. Similarly, a higher neighborhood percentage of Hispanic or Latino residents was associated with a greater total count of retailers, retailers per square mile, and retailers per 10 km of roadway in unadjusted models. In adjusted models, a higher neighborhood percentage of Hispanic or Latino residents was associated with fewer retailers per 1,000 people. The inverse adjusted associations between percentage Black and Hispanic or Latino residents and retailers per 1,000 people are consistent with a 2012 U.S. study limited to a sample of tracts within 97 counties (Lee et al., 2017). The tobacco industry has a long history of marketing tobacco products in the retail environment to Black and Latino and Hispanic individuals and communities (Iglesias-Rios & Parascandola, 2013; Kostygina et al., 2016; Yerger et al., 2007). Given associations of retailer availability and marketing with smoking behavior (Marsh et al., 2020; Robertson et al., 2015; Robertson et al., 2016; Valiente et al., 2020), greater tobacco retailer density in neighborhoods that have a higher percentage of these individuals is a public health concern.
We also found that for tracts with a 10–percentage point difference in vacant housing, the model-estimated density difference was 0.12 tobacco retailers per 1,000 people, almost 11% of the average density in the sample. This positive association of greater vacant housing and greater tobacco retailer density was consistent with two other studies (Fakunle et al., 2018; Lee et al., 2017). However, for two other measures of retailer density (total count, per km roadway), we documented inverse associations with vacant housing. Though this study is cross-sectional and we cannot make inferences about changes in retailer density or demographics over time, we posit one reason for these findings. Vacant housing may reflect a relatively low population count in a place that historically housed more people and may therefore be structurally equipped with enough roads to accommodate more people. In places with a high percentage of vacant dwellings, retailers may not have a desire to invest in neighborhoods where there may be little demand. Fewer retailers in neighborhoods with the same number of historical roads might produce a reduction in total count and roadway-based density measures over time. Yet a drop in retailers that also corresponds with a population decrease could actually increase population-based retailer density measures depending on the relative drop in each. Longitudinal research assessing changes over time is needed to further shed light onto these potential processes.
In this study, we use different measures of retailer density to help increase comparability of findings with past and future studies assessing inequities in tobacco retailer density. Taken together, our findings indicate that while there are neighborhood sociodemographic differences in the availability of tobacco retailers, the nature of these differences is dependent on the measure of retailer density used. This may be particularly critical for local and state jurisdictions that are currently using measures of retailer density in policies designed to reduce retailer inequities. For example, San Francisco, California, limits the number of tobacco retailers to 45 in each supervisorial district while Philadelphia, Pennsylvania, limits the number of tobacco retailers to one per 1,000 people per planning district (ChangeLab Solutions, 2019). Although statistical criteria can indicate which measure may best fit specific models, these analyses may not be generalizable to other data. More important, each measure may capture different aspects of the tobacco built and social environment. For example, per capita measures may reflect different levels of consumer demand for retailers, land area measures may also describe spatial proximity of consumers to retailers, and roadway measures may indicate the ease with which consumers can access retailers via existing infrastructure.
The four measures of retailer density used in this study are not exhaustive, and some studies have used other more computationally advanced measures, such as fixed and adaptive kernel density, which are too complex to be used in local policies. Additionally, proximity-based measures, such as the distance of tobacco retailers to one another, may be important for further understanding inequities in tobacco retailer availability. Which of these measures might most influence exposure to tobacco marketing, purchasing behavior, and tobacco use, however, requires greater research on how people interact within their activity spaces (Marsh et al., 2020; Valiente et al., 2020). Additionally, careful consideration should be taken when thinking about which measures might be most feasible to calculate or interpretable for communities and policymakers.
Given that policies regulating the sale of tobacco products in both pharmacies and tobacco shops have been implemented in several jurisdictions to reduce tobacco retailer density, we investigated whether neighborhood demographics were associated with the presence of tobacco-selling pharmacies or tobacco shops. We found that there was a lower odds of a tract having a tobacco-selling pharmacy in neighborhoods with a greater percentage of Black or Hispanic or Latino residents in adjusted models, suggesting that policies that restrict the sales of tobacco products in pharmacies may not decrease the number of tobacco retailers equitably across neighborhoods. These findings are consistent with other studies that have investigated the impact of a pharmacy ban on inequities in tobacco retailer density (Craigmile et al., 2020; Giovenco et al., 2019; Kulbicki & Leslie, 2015; Tucker-Seeley et al., 2016). While our study sample was limited to pharmacies that sell tobacco products, our findings are similar to those focused on pharmacy deserts, which find that pharmacies are less available in neighborhoods with a greater proportion of some racially and ethnically minoritized groups (Qato et al., 2014). Finally, the policy focus to date has been on restricting sales of tobacco products in pharmacies; however, one New Zealand modeling study considered the impact of only
Compared to White individuals who smoke, Black individuals who smoke are less likely to purchase cigarettes at tobacco shops (Groom et al., 2020), and our study found that the odds of a tract having at least one tobacco shop (vs. none) was lower in neighborhoods with a greater percentage of Black residents. Additionally, we found that the odds of a tract having a tobacco shop (vs. none) was higher in areas with a greater percentage of residents living below 150% FPL. Notably, the odds of having a tobacco shop was larger (a
Several considerations should be made when interpreting the results of this study. First, we created a probable list of tobacco retailers, and it is possible that our list contained retailers that do not sell tobacco, or there could be tobacco retailers missing. We have no reason to believe that this potential error is systematic. Second, as our analytic sample represents a near census of all tracts in the United States, we have high power to detect small associations, and caution should be taken when interpreting small effect sizes. On the other hand, because this is a near census of tracts, associations observed may be more likely to represent the “true” population parameter. Of importance is that we conceptualized the neighborhood as a census tract, but other neighborhood scales (e.g., block groups) may be appropriate. Additionally, there may be other important neighborhood sociodemographic characteristics that may contribute to the associations observed (e.g., commercial land use) that deserve future investigation. Finally, this study is cross-sectional, and temporality cannot be established. Regardless of temporality, tobacco products are not a health-promoting neighborhood commodity, and their availability and marketing could influence smoking behaviors.
In this national study, we use several common measures of tobacco retailer density to investigate associations with tract-level sociodemographic characteristics. While we document inequities in retailer density by area race, ethnicity, poverty, and vacant housing, these relationships were not consistent across all measures. Researchers and policymakers should consider how various measures of tobacco retailer density may capture different aspects of the tobacco retailer environment in their communities. Furthermore, attention to whether the availability of certain tobacco retailer types differs across and by neighborhood characteristics may be important when considering the varying impact of some tobacco retailer reduction policies. Overall, identifying the relationships between neighborhood sociodemographic characteristics and tobacco retailer availability may help communities better track place-based tobacco retailer inequities and design impactful pro-equity retailer-focused strategies.
Supplemental Material
sj-docx-1-heb-10.1177_10901981211008390 – Supplemental material for Neighborhood Inequities in Tobacco Retailer Density and the Presence of Tobacco-Selling Pharmacies and Tobacco Shops
Supplemental material, sj-docx-1-heb-10.1177_10901981211008390 for Neighborhood Inequities in Tobacco Retailer Density and the Presence of Tobacco-Selling Pharmacies and Tobacco Shops by Amanda Y. Kong, Paul L. Delamater, Nisha C. Gottfredson, Kurt M. Ribisl, Chris D. Baggett and Shelley D. Golden in Health Education & Behavior
Footnotes
Authors’ Note
The content is solely the responsibility of the authors and does not represent the official views of the National Institutes of Health.
Declaration of Conflicting Interests
The authors declared the following potential conflicts of interest with respect to the research, authorship, and/or publication of this article: KMR serves as an expert consultant in litigation against tobacco companies.
Funding
The authors disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was supported by the Paul C. Hardin Dissertation Completion Award from the Graduate School at the University of North Carolina at Chapel Hill and the National Cancer Institute of the National Institutes of Health (Award Numbers: F31CA239331, T32CA128582, P01CA225597). The content is solely the responsibility of the authors and does not represent the official views of the funders or the National Institutes of Health.
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
