 Methodology
 Open Access
 Published:
Quantifying spatial accessibility in public health practice and research: an application to onpremise alcohol outlets, United States, 2013
International Journal of Health Geographicsvolume 17, Article number: 23 (2018)
Abstract
Objective
To assess spatial accessibility measures to onpremise alcohol outlets at census block, census tract, county, and state levels for the United States.
Methods
Using network analysis in a geographic information system, we computed distancebased measures (Euclidean distance, driving distance, and driving time) to onpremise alcohol outlets for the entire U.S. at the census block level. We then calculated spatial accessbased measures, specifically a populationweighted spatial accessibility index and populationweighted distances (Euclidean distance, driving distance, and driving time) to alcohol outlets at the census tract, county, and state levels. A multilevel modelbased sensitivity analysis was conducted to evaluate the associations between different onpremise alcohol outlet accessibility measures and excessive drinking outcomes.
Results
The national average populationweighted driving time to the nearest 7 onpremise alcohol outlets was 5.89 min, and the average populationweighted driving distance was 2.63 miles. At the state level, populationweighted driving times ranged from 1.67 min (DC) to 15.29 min (Arizona). Populationweighted driving distances ranged from 0.67 miles (DC) to 7.91 miles (Arkansas). At the county level, populationweighted driving times and distances exhibited significant geographic variations, and averages for both measures increased by the degree of county rurality. The populationweighted spatial accessibility indexes were highly correlated to respective populationweighted distance measures. Sensitivity analysis demonstrated that population weighted accessibility measures were more sensitive to excessive drinking outcomes than were population weighted distance measures.
Conclusions
These results can be used to assess the relationship between geographic access to onpremise alcohol outlets and health outcomes. This study demonstrates a flexible and robust method that can be applied or modified to quantify spatial accessibility to public resources such as healthy food stores, medical care providers, and parks and greenspaces, as well as, quantify spatial exposure to local adverse environments such as tobacco stores and fast food restaurants.
Background
Quantifying spatial accessibility in public health practice is essential for evaluating population exposure to local environments (e.g., alcohol and tobacco outlets or public parks) and population access to health care resources (e.g., primary care clinics or trauma hospitals). Three approaches have been reviewed by Zhang et al. [1] and summarized for application to the measurement of alcohol outlet density by CDC [2]. Table 1 summarizes the common metrics for quantifying spatial accessibility and basic relationships among them. Two commonly used approaches for quantifying spatial accessibility are distancebased and containerbased. The spatial interaction (or gravity) modelbased spatial accessibility index uses both distancebased metrics and the container concept. Populationweighted accessibility metrics are based on the spatial interaction modelbased spatial accessibility index, but at the same time account for the uneven population distribution within a study area. Populationweighted distance metrics use distancebased metrics and at the same time use the spatial interaction modelbased spatial accessibility index to construct differential probability access to destinations and further account for the uneven population distribution within a study area. The most complex populationweighted distance metrics aim to borrow the strengths and at the same time minimize the limits of classic distancebased and containerbased metrics, and integrate the power and flexibility of the spatial interaction modelbased approach and account for uneven population distributions.
Distancebased measures, such as distance to nearest outlets, are intuitive and relatively easy to generate, but often lack the ability to quantify and incorporate the potential geographic clustering effects of spatial destinations. The three most common distancebased metrics are listed in Table 1.
The containerbased approach is most commonly used to assess outlet density for a predefined area (container). Some common containerbased metrics are the number of outlets per square mile or per road mile, or within a predefined spanning distance (Table 1), but often there is an underlying assumption of equal accessibility within this predefined area. There could also be substantial boundary or edge effects, such that outlets near the predefined area boundaries could be either included or excluded in the density calculation. The spatial interaction model–based approach assumes the spatial association or interaction between spatial or geographic entities are proportional to their mass sizes and inversely proportional to their distances. It takes four basics steps to construct a spatial accessibility index, as detailed in Table 1.
The spatial interaction modelbased spatial accessibility index is often used by geographers because of its flexibility and robustness in quantifying spatial accessibility. However, interpretations of these indices are not intuitive for public health practitioners and communities. Populationweighted accessibility (PWA) metrics have similar interpretation challenges in public health practice.
The populationweighted distance (PWD) metrics, a form of spatial accessibility measure developed by Zhang et al. [1] accounts for spatial clustering of outlets, overcomes the unrealistic equal access assumption and potential edge effects of the containerbased approach, and also uses an intuitive form of a distancebased measure. In addition, PWD also accounts for uneven population distribution for the geographic area of interest. The flexibility of PWD measures makes them applicable to individual persons or households, as well as for a geographic area. The individualbased PWD measures could be applied in specific public health analyses with individual records. The arealevel PWD measures are based on the smallest census geographic units (census blocks); thus, PWD measures could be conveniently aggregated to any larger geographic levels as needed while avoiding the modifiable unit area problem [3]. Despite these potential advantages, PWD measures may be sensitive to the choice of the original distance measures upon which the resulting populationweighted measures are computed. The spatial interaction model to calculate populationweighted accessibility and accessibility index is described in Zhang et al. [1].
In the United States, the census block is the basic unit of census geographic hierarchy (https://www.census.gov/geo/reference/hierarchy.html). We used census blocklevel population as our demand area population. By using the population at finest level of census geography that is available in the United States, we minimized the geographic aggregation error for population data described as Source A type of error in Hillsman and Rhoda [4]. Spatial access metrics based on higher levels of census geography than census blocks (e.g. census tract, ZIP Code, or county), even by using population weighted centroids to estimate demand population locations, could still introduce substantial spatial aggregation errors or Source A type errors. In addition, the census blockbased spatial access metrics provide the flexibility to aggregate the metrics to any highlevel geography that could be linked to geocoded individual or aggregated health outcomes. For example, we could easily aggregate census blockbased spatial access metrics to higherlevel census geographic units, such as census tract, county, ZIP Code, school district, congressional district, or other administrative or customized geographic units, such as hospital service areas, hospital referral regions, and primary care service area. However, there were no studies that quantify spatial accessibility based on different census blocklevel distance measures (Euclidian distance, driving distance and driving time) and evaluate their correlations between populationweighted spatial accessibility index and populationweighted distance metrics.
In this study, we selected onpremise alcohol outlets because of significant public health interest in excessive alcohol use. Greater alcohol outlet density is associated with increased alcohol consumption and related harms [5]. Excessive alcohol consumption is responsible for 88,000 deaths annually in the United States [6] and accounted for $249 billion in economic costs in 2010 [7].
An alcohol outlet is defined as a place where alcohol may be legally sold to a buyer to consume there (e.g., onpremise outlets, such as bars or restaurants) or elsewhere (e.g., offpremise outlets, such as liquor stores) [5]. The onpremise alcohol outlet data that we obtained from HSIP GOLD is a point dataset that has been geocoded to the street address of the outlets. The pointlevel data eliminated a location error when geocoding data to ZIP Code area centroid. Regulating alcohol outlet density, which is typically defined as the number of alcohol outlets in a given area, is an effective strategy for the prevention of excessive alcohol consumption and related harms [8]. The three general approaches for quantifying spatial accessibility described previously can be applied to the measurement of alcohol outlets [2]. In the containerbased approach, the number of alcohol outlets in a given area is calculated. Different denominators can be used as follows: by population size (number of outlets per 1000 population); by area size (number of outlets per square mile); or by road mile length (number of outlets per road mile). Commonly, measurements of alcohol outlet density have used containerbased approaches, in which the number of outlets is divided by the population size of a particular area or by the land area itself [9]. However, these approaches could suffer from boundary and edge effects, especially for small geographic areas, and do not directly consider the spatial accessibility between alcohol outlets and the population [2]. It can result in overestimates for small areas with large numbers of alcohol outlets, or underestimates for small areas with small numbers or no alcohol outlets.
The distancebased approach calculates the distance to the nearest specified number of alcohol outlets. The distance can be calculated based on Euclidean distance, driving distance through a street network, and driving time through a street network with speed limits. The advantage to this method is that it is intuitive for pointtopoint units; however, it is less commonly used in alcohol outlet density studies. The limitation for the distancebased approach is that it ignores the possibility of multiple potential alcohol outlet destinations within an area. An alternative approach is to average the distances to a set of nearest (e.g. five or seven) alcohol outlets. However, this approach suffers from ignoring the unequal probability of access to any outlet. The populationweighted distance approach uses population as a weight to account for uneven population distributions. It incorporates unequal probability of access to nearby alcohol outlets, defines the choice set of alcohol outlets guided by human behavioral theory, and is able to account for alcohol outlet size (if data are available) and spatial clustering [9].
The objectives of this study are to use a comprehensive dataset of locations of onpremise alcohol outlets in the United States to create a populationweighted accessibility index and populationweighted distance metrics to quantify population access to onpremise alcohol outlets at the census block, census tract, county, and state levels for the entire U.S. Specifically, we use three different distance metrics: Euclidian (flight) distance, driving distance and driving time, to construct a spatial accessibility index to generate three populationweighted accessibility metrics and three more intuitive populationweighted distance metrics. We then compared and evaluated their similarities and correlations. This spatial accessibility modeling framework could be conveniently applied to any other countries or places with hierarchical census geography, such as, census output areas, super output areas, electoral wards, and electoral divisions in the UK; villages, street districts, towns, counties, cities, and provinces in China.
Methods
Data sources

1.
Onpremise alcohol outlets
We obtained geocoded data for onpremise alcohol outlets from the Homeland Security Infrastructure Program (HSIP) GOLD 2013 dataset [10]. It included 210,482 drinking establishments (HSIP terminology to represent onpremise alcohol outlets) in the 50 States and DC, with state averages ranging from 435 (DE) to 16,942 (TX) drinking establishments. The North American Industry Classification System (NAICS) is the standard used by Federal statistical agencies in classifying business establishments for the purpose of collecting, analyzing, and publishing statistical data related to the U.S. business economy. We used NAICS code 7224 and all subcategory NAICS codes of 7224 for Drinking Places (Alcoholic Beverages). The types of onpremise alcohol outlets in these NAICS codes include bars and lounges, drinking places, nightclubs, eating places, restaurants, hotels and motels, bowling centers, and recreation clubs. This dataset provided the exact street addresses of onpremise alcohol outlets.

2.
Geographic unit of analysis and population demographics
We used 2010 census blocks as our geographic unit of analysis. Census blocks are the smallest unit in the census geographic hierarchy, thus census blocklevel measures could be aggregated to any geographic levels of interest in a very flexible way and minimize the effects of the modifiable area unit problem. We included 6,207,027 census blocks with 2010 populations greater than zero in this study. The geometric centroid of the census block was used to calculate the distance to the location of onpremise alcohol outlets. The 2010 census population count was used as the weight to calculate the populationweighted spatial accessibility index, populationweighted driving distance, populationweighted driving time and populationweighted Euclidean distance at each geographic level.

3.
Street network data
We used Esri Data and Maps 10.2 U.S. and Canada Detailed Streets network dataset (Esri, Redlands, CA). These data provide distance and travel times for each street segment that was used in our network analyses.
Measures

1.
Distance metrics
We used ArcGIS 10.2.1 Network Analyst (Esri, Redlands, CA) to calculate the nearest 7 onpremise alcohol outlets for each populated census block. This tool produced a driving distance and a driving time for each origin (census block) and destination (onpremise alcohol outlet) pair. There were 2118 census blocks that did not have access to a street network, and thus were excluded in the driving distance and time calculations. We also calculated the Euclidian distance for each census block to its 7 nearest onpremise alcohol outlets using all census blocks in the U.S. We used SAS GEODIST function to calculate the Euclidian distance between census block centroids and alcohol outlets (SAS, Cary, North Carolina).

2.
Populationweighted spatial accessibility index
We used the spatial interaction model to calculate the populationweighted spatial accessibility index as described in the introduction. Because the size of each onpremise alcohol outlet is unknown, we used a weight of 1 to treat each outlet equally. The Census blocklevel spatial accessibility index was calculated as:
where i is ith census block, j is jth onpremise alcohol outlet and n = 7 was used for 7 nearest outlets in our study. The populationweighted spatial accessibility index for geographic units larger than census blocks was calculated as:
where N is the number of census blocks within a geographic area k (e.g. county), Pop_{ i } is the ith census block’s total population, and Pop_{ k } is the total population for the geographic area k.

3.
Populationweighted driving distance/time, Euclidian distance
We used distances and times from the distancebased measures along with blocklevel populations to calculate the populationweighted driving distance and populationweighted driving time for each census block. Because of the massive size of the data involved in the study, we stratified the U.S. census blocks into 17 smaller regions for computational purposes to perform the network analysis and then merged the results into one data file.
Populationweighted driving distance/time for census block i to 7 nearest onpremise alcohol outlets is calculated as:
where P_{ij} is the probability that a resident at a census block i will choose to visit an alcohol outlet j:
The probability was derived from Huff Model 1963 [11].
Populationweighted driving distance/time for area k is calculated as:
Using these methods, we calculated the populationweighted accessibility (PWA) index and populationweighted distance (PWD) distances (by driving distance, driving time, and Euclidian distance) to the 7 nearest onpremise alcohol outlets at the census block level. Then, we calculated the same spatial accessibility measures at U.S., state, county and census tractlevels respectively.
Data analysis
We compared the 6 spatial accessibility measures using Pearson correlation coefficients at the census block level. Specific measures that we analyzed—all at the county level—were driving distance, driving time, and Euclidean distance to the one nearest onpremise alcohol outlet; driving distance, driving time, and Euclidean distance to the nearest 7 onpremise alcohol outlets in A_{ij} = 1/(d _{ij} ^{^β} ) with a distance decay parameter β = 1.
We used countylevel spatial accessibility results linked with the 2013 National Center for Health Statistics (NCHS) urban–rural classification scheme for counties [12] and calculated spatial accessibility by urban and rural counties to determine whether urban/rural status was associated with spatial accessibility measures. The NCHS urban–rural classification for counties classified U.S. counties into the following 6 categories: (1) large central metro, defined as counties in a metropolitan statistical area (MSA) of at least 1 million residents that contain the entire population of the largest principal city of the MSA, or are completely contained within the largest principal city of the MSA, or contain at least 250,000 residents of any principal city in the MSA; (2) large fringe metro, defined as counties in an MSA of 1 million or more residents that does not qualify as a large central metro; (3) medium metro, defined as counties in an MSA of 250,000–999,999 residents; (4) small metro, defined as counties in an MSA of less than 250,000 residents; (5) micropolitan, defined as counties in a micropolitan statistical area; and (6) noncore (often called rural), defined as counties not in MSAs or micropolitan statistical areas. Categories 1–4 are metropolitan counties and categories 5–6 are nonmetropolitan counties.
Sensitivity analysis
We conducted a sensitivity analysis via linking population weighted spatial accessibility metrics to drinking outcomes from a nationwide health survey. We first calculated the county level population weighted spatial accessibility measures based on driving distance, driving time and Euclidian distance to the nearest 7 onpremise alcohol outlets. We then linked these county level alcohol outlet accessibility measures to individual binge drinking (4 or more drinks on an occasion for women, 5 or more drinks on an occasion for men) and heavy drinking outcomes (8 or more drinks per week for women, 15 or more drinks per week for men) from the restricted 2013 Behavioral Risk Factor Surveillance System (BRFSS) data with survey respondents’ county residence identifiers. With this linked dataset, we constructed a series of multilevel logistic models for both drinking outcomes and estimated the odd ratios associated with alcohol outlet access measures. These multilevel logistic models also included individual age, gender, race/ethnicity, education and income; countylevel poverty; and county and statelevel contextual effects (specified by countylevel and statelevel random effects). The adjusted ORs were used to detect the significant associations between alcohol outlet access and drinking outcomes and to measure their sensitivity to drinking outcomes. A OR larger than one means an increased risk for drinking behaviors, and a larger OR means more sensitive to drinking outcomes. SAS Proc GLIMMIX was used to implement these multilevel logistic models.
Results
In 2013, the average driving time to the nearest onpremise alcohol outlets (of the nearest 7 choices) for the overall U.S. population was 5.89 min and the average driving distance was 2.63 miles. For all census blocks, the median populationweighted driving time was 4.72 min; populationweighted driving distance was 1.99 miles; and populationweighted Euclidian distance was 1.36 miles (Table 2).
At the statelevel (Table 3), populationweighted driving time ranged from 1.67 (DC) to 15.29 min (AR), with a median of 6.67 min (LA); populationweighted driving distance ranged from 0.69 (DC) to 7.91 (AR) miles, with a median of 2.95 miles (LA). In general, states in the northeast, the upper Midwest, and the western U.S. had a driving distance of less than 2 miles to the nearest 7 onpremise alcohol outlets (Fig. 1). The same pattern was observed for driving time (not shown). At the county level, populationweighted driving distance (Fig. 2) and populationweighted driving time (not shown) varied greatly by geography. Shorter populationweighted distances were observed in counties along the northeast coast, around the Great Lakes, in Florida, portions of the Rocky Mountain states, along the west coast, and in large metropolitan areas. Figure 3 presents detailed populationweighted driving distances at the census tract level. All data for 6 measures at state, county and census tractlevel are included in Additional file 1.
The mean estimates of populationweighted driving time, driving distance, and Euclidian distance all increased from large central metro to noncore categories (Table 4). Populations living in medium metro counties had shorter median driving time/distance and Euclidian distance than those living in large fringe metro counties. For the minimum of countylevel populationweighted driving time, driving distance, and Euclidian distance, populations living in small metro counties had shorter driving time/distance and Euclidian distance than those living in medium metro counties. Noncore counties had shorter minimum distances than micropolitan counties, for both driving distance and Euclidian distance, but noncore counties had longer minimum driving time than micropolitan counties. For the maximum of countylevel populationweighted driving time, driving distance and Euclidian distance, driving time increased with increased rurality; for driving distance, small metro and micropolitan counties had shorter distances than medium metro counties; small metro counties had shorter Euclidian distances than medium metro counties.
Populationweighted distance measures based on driving time, driving distance, or Euclidian distance were strongly correlated with each other (p < 0.01) (Table 5). Additionally, populationweighted accessibility measures were also strongly correlated with each other (p < 0.01). However, populationweighted distance measures and accessibility measures were moderately correlated to each other (p < 0.01).
Table 6 presents the odds ratios (ORs) of excessive drinking for an interquartile range increase in populationweighted accessibility measures and populationweighted distance measures for onpremise alcohol outlets. The interquartile ORs are significantly larger than one, which means that an increased alcohol outlet access were associated with significantly increased risks for excessive drinking in the United States. In other words, the odds ratios represent the minimum risk increase when a person moves from one county with a lower access (first quartile) to alcohol outlets to another county with a higher access (fourth quartile) to alcohol outlets. Population weighted accessibility measures detected stronger associations (larger OR) with excessive drinking outcomes and were more sensitive to alcohol outcomes than population weighted distance measures. Among these spatial accessibility measures, driving distance/time and Euclidian distance had similar sensitivities to binge drinking/heavy drinking. Among the populationweighted distance measures, driving distance/timebased ones were much more sensitive to binge drinking/heavy drinking than one based on Euclidian distance.
Discussion
Our study has demonstrated a robust approach to quantify spatial accessibility of geographic entities from a population health perspective. The populationweighted bottomup approach provides great flexibility in generating spatial accessibility measures at any geographic level that could be linked with population health outcomes of interest. This is the first nationwide networkbased analysis of spatial accessibility to onpremise alcohol outlets at the level of the census block and their aggregated measures provide a more holistic and accurate picture of population spatial accessibility to onpremise alcohol outlets from the local to the national level.
Census blockbased populationweighted distance measures are very flexible and can be aggregated to any geographic level as needed (e.g., ZIP Code, neighborhood). The method accounts for uneven local population distributions, reduces the ecological bias in measuring alcohol outlet density, and is flexible in that it can incorporate more information when detailed data are available (e.g. size of the outlets, sales of the outlets). Populationweighted accessibility measures can be linked with other health outcomes at different geographic levels to study the relationship of alcohol access to related harms. Additional analysis is needed to evaluate its sensitivity to excessive drinking outcomes and related population health outcomes. Additional analysis is also needed to evaluate the difference between commercial datasets of licenses alcohol outlets with local alcohol license data.
In spatial interaction modeling and spatial choice modeling, how to characterize destination or spatial choice set is still challenging and has no certain answer. We used 7 onpremise alcohol outlets in our model based on Miller et al. and Saaty et al. [13, 14]. We give different probabilities for those 7 outlets based on the Huff model [4]. In this approach, the population in any given census block has a higher probability of accessing a nearer rather than farther onpremise alcohol outlet, and it also accounts for people who may not always chose the nearest outlet.
This study has several limitations. First, alcohol outlet size was treated as one unit (e.g., all alcohol outlets are equally weighted). Second, commercial datasets may not be updated frequently enough to reflect local businesses opening and closing, which may result in over or underestimations of spatial accessibility to onpremise alcohol outlets in some locations. Third, we were aware that people often purchase alcohol from offpremise alcohol outlets (stores) and consume the alcohol at home. Due to lack of data access, we could not include offpremise alcohol outlet in our analysis. This analysis is not meant to present the entire picture of access for both on and offpremise alcohol outlets. Further analysis may be conducted to include offpremise alcohol outlets in the model to obtain a more comprehensive measurement of the impact of alcohol outlet access, once such data becomes available. Finally, census block population counts are only available every 10 years when the decennial census is conducted. For the nondecennial years, census block population could be updated via small area population estimation techniques.
Populationweighted spatial accessibility index and populationweighted distance metrics are moderately correlated, which indicates that the populationweighed spatial accessibility index is conceptually different from populationweighted distance metrics and captures different aspects or dimensions of spatial accessibility to target destinations. Epidemiologically, we can expect that the populationweighted accessibility measures and populationweighted distance measures could have differential sensitivity to population health outcomes. Our sensitivity analysis between access to onpremise alcohol outlets and excessive drinking suggested that lessintuitive spatial accessibility measures are much more sensitive than those moreintuitive direct distance measures. It suggests that some conclusions in the current literature on alcohol outlet density and adverse health outcomes could underestimate the magnitude of associations between alcohol outlet density and harmful drinking behaviors. These moresensitive alcohol outlet density measures could be used in public health impact studies on alcohol outlets, provided data with the necessary spatial granularity are available. Alcohol use is a complicated human behavior and access to onpremise alcohol outlets is only one of many potential explanatory factors. It is also challenging to quantify people’s preference and selection of onpremise alcohol outlets. Our study method accounts for the probability that closer outlets get more access than more distant outlets, but also takes into consideration that people may not always go to the closest outlets as indicated by a recent USDA report (EIB138) [15].
Conclusions
These results can be used to assess the relationship between geographic access to onpremise alcohol outlets and health outcomes. Additionally, this study demonstrated a flexible and robust method that can be applied or modified to quantify spatial accessibility to public resources, such as stores that serve healthy food options, medical care providers, and parks and greenspaces, and spatial exposure to local adverse environments such as tobacco stores and fast food restaurants. This spatial accessibility modeling framework could be conveniently applied to any countries or places with hierarchical census geography. In modern society, census geographic hierarchy is the base for socioeconomic and population statistics. Our bottomup approach within a census geographic hierarchy could be widely applicable and substantially improve the accuracy and flexibility in quantifying population spatial access to public health resources or exposure to environmental hazards of interest.
Abbreviations
 CDC:

The Centers for Disease Control and Prevention
 PWA:

populationweighted accessibility
 PWD:

populationweighted distance
 HSIP:

The Homeland Security Infrastructure Program
 NAICS:

The North American Industry Classification System
 NCHS:

The National Center for Health Statistics
 MSA:

metropolitan statistical area
References
 1.
Zhang X, Lu H, Holt JB. Modeling spatial accessibility to parks: a national study. Int J Health Geogr. 2011;10:31.
 2.
Centers for Disease Control and Prevention. Guide for measuring alcohol outlet density. Atlanta: Centers for Disease Control and Prevention, US Department of Health and Human Services; 2017. https://www.cdc.gov/alcohol/researchinaction.html. Accessed 30 Aug 2017.
 3.
Openshaw S. The modifiable areal unit problem. Norwick: Geo Books; 1984.
 4.
Hillsman EL, Rhoda R. Errors in measuring distances from populations to service centers. Ann Reg Sci. 1978;12:74. https://doi.org/10.1007/BF01286124.
 5.
Campbell CA, Hahn RA, Elder R, Brewer R, Chattopadhyay S, Fielding J, Naimi TS, Toomey T, Lawrence B, Middleton JC, Task Force on Community Preventive Services. The effectiveness of limiting alcohol outlet density as a means of reducing excessive alcohol consumption and alcoholrelated harms. Am J Prev Med. 2009;37(6):556–69.
 6.
Stahre M, Roeber J, Kanny D, Brewer RD, Zhang X. Contribution of excessive alcohol consumption to deaths and years of potential life lost in the United States. Prev Chronic Dis. 2014;11:130293.
 7.
Sacks JJ, Gonzales KR, Bouchery EE, Tomedi LE, Brewer RD. 2010 National and state costs of excessive alcohol consumption. Am J Prev Med. 2015;49(5):e73–9.
 8.
Task Force on Community Services. Recommendations for reducing excessive alcohol consumption and alcoholrelated harms by limiting alcohol outlet density. Am J Prev Med. 2009;37(6):570–1.
 9.
Zhang X, Hatcher B, Clarkson L, Holt J, Bagchi S, Kanny D, Brewer RD. Changes in density of onpremises alcohol outlets and impact on violent crime, Atlanta, Georgia, 1997–2007. Prev Chronic Dis. 2015;12:140317.
 10.
Homeland Security Infrastructure Program (HSIP) https://www.dhs.gov/infrastructureinformationpartnerships. Accessed 20 June 2018.
 11.
David LH. A probabilistic analysis of shopping center trade areas. Land Econ. 1963;39(1):81–90.
 12.
Ingram DD, Franco SJ. 2013 NCHS urban–rural classification scheme for counties. Vital Health Stat 2. 2014;166:1–73. https://www.cdc.gov/nchs/data_access/urban_rural.htm. Accessed 30 June 2017.
 13.
Miller GA. The magical number seven, plus or minus: some limits on our capacity for processing information. Psychol Rev. 1956;63(2):81–97.
 14.
Saaty TL, Ozdemir MS. Why the magic number seven plus or minus two. Math Comput Model. 2003;38(3–4):233–44.
 15.
Ver Ploeg M, Mancino L, Todd JE, Clay DM, Scharadin B. Where do Americans usually shop for food and how do they travel to get there? Initial findings from the national household food acquisition and purchase survey. Economic Information Bulletin (EIB138), Economic Research Service, U.S. Department of Agriculture. March 2015.
Authors’ contributions
HL, XZ, JBH conceived the study and HL did all the network analysis and made all maps, XZ calculated PWD, PWA and statistical analysis. HL, XZ, JBH, DK and JBC wrote the draft and final version of the manuscript. All authors read and approved the final manuscript.
Authors’ information
The findings and conclusions in this report are those of the authors and do not necessarily represent the official position of the Centers for Disease Control and Prevention or Economic Research Service, U.S. Department of Agriculture.
Competing interests
The authors declare that they have no competing interests.
Availability of data and materials
The datasets (HISP and Esri Street Network) analyzed during the current study are not publicly available due to data license restriction. Census 2010 block level population can be download at www.census.gov. The datasets generated from this study are included in this published article and its additional file.
Consent for publication
Not applicable.
Ethics approval and consent to participate
Not applicable.
Funding
All authors are U.S. Government employee, therefore the work was founded by U.S. Government.
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Author information
Additional file
12942_2018_143_MOESM1_ESM.xlsx
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), 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 (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Received
Accepted
Published
DOI
Keywords
 Network analysis
 Spatial accessibility
 Alcohol outlet