Commute Times, Food Retail Gaps, and Body Mass Index in North Carolina CountiesPosted on by
Volume 7: No. 5, September 2010
Commute Times, Food Retail Gaps, and Body Mass Index in North Carolina Counties
Stephanie B. Jilcott, PhD; Haiyong Liu, PhD; Justin B. Moore, PhD, MS; Jeffrey W. Bethel, PhD; James Wilson, PhD; Alice S. Ammerman, DrPH, RD
Suggested citation for this article:
Jilcott SB, Liu H, Moore JB, Bethel JW, Wilson J, Ammerman AS. Commute times, food retail gaps, and body mass index in North Carolina counties. Prev Chronic Dis 2010;7(5). http://www.cdc.gov/pcd/issues/2010/sep/09_0208.htm. Accessed [date].
The prevalence of obesity is higher in rural than in urban areas of the United States, for reasons that are not well understood. We examined correlations between percentage of rural residents, commute times, food retail gap per capita, and body mass index (BMI) among North Carolina residents.
We used 2000 census data to determine each county’s percentage of rural residents and 1990 and 2000 census data to obtain mean county-level commute times. We obtained county-level food retail gap per capita, defined as the difference between county-level food demand and county-level food sales in 2008, from the North Carolina Department of Commerce, and BMI data from the 2007 North Carolina Behavioral Risk Factor Surveillance System. To examine county-level associations between BMI and percentage of rural residents, commute times, and food retail gap per capita, we used Pearson correlation coefficients. To examine cross-sectional associations between individual-level BMI (n = 9,375) and county-level commute times and food retail gap per capita, we used multilevel regression models.
The percentage of rural residents was positively correlated with commute times, food retail gaps, and county-level BMI. Individual-level BMI was positively associated with county-level commute times and food retail gaps.
Longer commute times and greater retail gaps may contribute to the rural obesity disparity. Future research should examine these relationships longitudinally and should test community-level obesity prevention strategies.
In the United States, the prevalence of obesity is higher in rural than in urban populations (1-5). Area-level factors that contribute to this disparity are not well understood, but one underlying mechanism may be the food environment. Obesity prevalence is lower in census tracts containing a supermarket (6), and rural areas have few supermarkets, which generally have a healthier mix of low-cost food items compared with local convenience stores (7). Accessibility to healthy food is also difficult in rural areas because convenience stores are more common than supermarkets (8-10).
Rural residents may regularly travel to urban areas in neighboring counties to shop for food because of convenience along the route to work, better prices, wider selection, or one-stop shopping offered at discount “supercenters” (eg, Walmart) (11,12). This pattern of food shopping among rural residents may create a retail shortfall or “gap” for food venues in rural areas, causing rural food venues to have a decreased share of the market. A large food retail gap may exacerbate rural food deserts (13), or areas where residents have limited access to affordable, healthy food (14), when smaller food venues in underserved areas close as business is lost to nearby discount supercenters (13,15). Rural residents’ prolonged travel time to larger supermarkets or supercenters not only increases the retail gap in the rural county but decreases the frequency of food shopping. In turn, diet quality may decrease as rural residents purchase less fresh produce and more processed foods (16,17).
Another hypothesized mechanism underlying the rural-urban obesity disparity is that rural residents may spend more time traveling to work or to obtain goods and services than do their urban counterparts. Obesity is associated with urban sprawl (18-20), time spent in cars (21), and vehicle miles traveled per day (22). One Los Angeles-based study found that distance traveled to the nearest supermarket was positively associated with higher body mass index (BMI) (23). To our knowledge, no studies have examined the associations between distance to food shopping location, commute times, and BMI among rural and urban residents.
To better understand associations between area-level factors and obesity, we conducted ecologic analyses of associations between the percentage of rural-dwelling residents, commute times, food retail gaps, and BMI for all 100 North Carolina counties. We hypothesized that 1) the percentage of rural residents per county is positively correlated with commute time and food retail gap per capita, 2) county-level commute time is positively correlated with food retail gap per capita, and 3) both commute time and food retail gap per capita are positively correlated with county-level mean BMI. In separate individual-level, contextual analyses, we examined individual-level BMI as the dependent variable and county-level commute time and food retail gap per capita as independent variables. We hypothesized that longer commute times and greater food retail gaps per capita would be positively correlated with individual-level BMI.
Percentage of rural residents
We calculated the percentage of rural residents for all North Carolina counties by dividing the number of county residents who lived in a rural area according to 2000 census criteria (24) by the county population. The percentage of rural residents ranged from 4% to 100%.
We generated reports for county-level commute times for 1990 and 2000 from US census data from the North Carolina Department of Commerce Economic Development Intelligence System. Census data were derived from answers to the census long-form questionnaire. Respondents who worked outside the home estimated the number of minutes it took to get from home to work each day, and commute time was derived by dividing the total number of minutes by the number of workers aged 16 years or older who did not work at home. We examined associations by using the 1990 and 2000 commute times and the difference in commute times between 1990 and 2000. The difference in 1990 and 2000 commute times describes broad shifts in county-level commuting over 10 years.
Food retail gap
We defined the food retail gap as the difference between county-level demand for food and county-level sales of food. We obtained the food retail gap for each North Carolina county from the North Carolina Department of Commerce Economic Development Intelligence System. The Environmental Systems Research Institute (ESRI) calculated retail gaps by subtracting county-level retail sales (supply) of products for a particular industry category in 2008 from county-level demand for products in that industry category in 2008. ESRI estimated demand using data on consumer expenditures from the Bureau of Labor Statistics and InfoUSA, a commercial database marketing system.
ESRI calculates the food retail gap for North American Industry Classification System (NAICS) codes 445 (representing food and beverage stores) and 722 (representing the food services and drinking places) separately. For these analyses, we used food retail gaps calculated from individual and combined NAICS codes 445 and 722. Venues included in the food and beverage stores subsector (NAICS code 445) sell food and beverages from fixed point-of-sale locations, such as supermarkets, grocery stores, convenience stores, meat markets, produce markets, and specialty food stores. Venues included in the food services and drinking places subsector (NAICS code 722) prepare meals, snacks, and beverages to customer order for consumption on and off the premises, such as full-service restaurants, limited-service eating places (fast-food restaurants), special food services, and drinking places. To control for population density, we calculated the food retail gap per capita by dividing the ESRI-estimated food retail gap by the 2007 county population estimate provided by the US census. A negative retail gap indicated that county-level sales were greater than county-level demand; a positive retail gap indicated that county demand was greater than county sales. For example, if County X has 1 chain supermarket and neighboring County Y has a large discount supercenter, residents of County X may begin grocery shopping at the supercenter, creating a positive food retail gap in County X and a negative food retail gap in County Y as residents’ food dollars are spent in the neighboring county. This could result in closing of the 1 chain supermarket in County X, making travel to the discount supercenter a necessity for obtaining groceries.
Body mass index
We estimated county-level mean BMI using self-reported height and weight for respondents to the North Carolina Behavioral Risk Factor Surveillance System (BRFSS); responses were aggregated over 5 years (2003-2007). The 5-year aggregate provided an adequate number of responses for reliable estimates for counties with low population densities (single-year estimates for rural counties are unstable). We calculated mean weighted BMI using SUDAAN version 10.1 (Research Triangle Institute, Research Triangle Park, North Carolina), which accounts for BRFSS oversampling of minorities. The mean (standard deviation) county-level BMI was 27.7 (0.85) kg/m2. The median (interquartile range) was 27.6 (25.9-30.1) kg/m2.
We conducted individual-level, contextual analyses using data from the 2007 North Carolina BRFSS for respondents aged 18 to 65 years with valid county identifiers. Because of confidentiality concerns, BRFSS does not provide county identifiers for residents of counties with fewer than 50 respondents. We excluded those counties. The individual-level sample consisted of 9,375 respondents from 64 counties. The mean population of the 64 counties included was 123,968, and the mean population of the 36 counties excluded was 25,322. The individual-level mean (SD) BMI was 28.1 (6.4) kg/m2.
County-level census data
To control for economic interdependence of adjacent counties, we examined the Rural to Urban Continuum Codes (RUCC) as a covariate. The RUCC is a 9-level ordinal scale used by the Economic Research Service to classify counties according to adjacency to metropolitan areas (24). We included a diversity index as a potential covariate in analyses because of associations between racial/ethnic mix and availability of food venues (eg, supermarkets , fast-food restaurants ) and to account for North Carolina counties’ varied race/ethnicity distributions (27). The diversity index represents the percentage of times 2 randomly selected people in each county would differ by race/ethnicity (27). The index is calculated by squaring the proportions of residents in each racial/ethnic group, summing the squares, and subtracting the result from 1. We determined both the county-level diversity index and the percentage of residents who lived below the poverty level using 2000 census data. We calculated the percentage of residents who lived below the poverty level by dividing the number of residents below the poverty level in 1999 by the estimated 1999 county population. North Carolina is divided into 3 regions (Coastal Plain/Eastern, Appalachian Mountain/Western, and Piedmont Plateau) with distinct demographic and socioeconomic characteristics. Thus, we also examined the variable “region” as a potential covariate.
For county-level ecological analyses, we used SAS version 9.2 (SAS Institute, Inc, Cary, North Carolina) to calculate correlation coefficients for percentage of rural residents, food retail gap per capita, commute time, and BMI for all 100 North Carolina counties. We used backward selection to construct linear regression models to examine the associations among county-level independent variables of commute times and food retail gap per capita, using county-level mean BMI as the dependent variable. Percentage of rural residents, diversity index, percentage below poverty, and region were potential covariates and were eliminated from the model in successive steps if the P value for the parameter estimate was .05 or higher. We examined the potential multicollinearity among covariates by computing their corresponding tolerance values. The tolerance is the proportion of variance in a given independent variable that is not explained by all of the other covariates; we found a tolerance value for all of greater than 0.1, which has been widely used as the threshold for multicollinearity in linear regressions (28).
For individual-level, contextual analyses, we constructed multilevel linear regression models; the dependent variable was individual-level BMI from 2007 BRFSS respondents (n = 9,375). County-level independent variables were food retail gap per capita, commute time in 2000, and difference in commute times between 2000 and 1990. Sex, age, race/ethnicity, and education were individual-level covariates, and the RUCC was a county-level covariate. Region was added as a third level.
Multilevel regression analyses allowed us to assess associations between individual-level BMI and area-level factors, accounting for the fact that people who reside in the same county are not independent observations (29). We examined the association between individual-level BMI and the 5 county-level variables of interest (commute time in 2000, difference between 1990 and 2000 commute times, retail gap per capita for NAICS code 445 [food and beverage stores], retail gap per capita for NAICS code 722 [restaurants and drinking places], and combined retail gap per capita) in separate models. The first 3 models to examine the association between BMI and county-level variables of interest were 2-level random intercept models. Model 4 included additional regional dummy variables to account for fixed effects from region. We used SAS version 9.2 for individual-level, contextual analyses, with estimates weighted to adjust for BRFSS oversampling.
Summary statistics for the individual-level data among 2007 BRFSS respondents by region are reported in Table 1.
Percentage of rural residents was significantly correlated with both the commute times in 1990 and 2000 and the difference in commute times between the 2 years, food retail gap per capita for restaurants and drinking places, overall food retail gap per capita, and BMI (Table 2).
We found significant positive correlations between commute time and retail gap per capita (Table 3). There were significant positive correlations between total food retail gap per capita and BMI and between the difference in commute times from 1990 to 2000 and BMI.
In linear regression analyses adjusted for county-level diversity index and the percentage of residents below the poverty level, a positive association was found between commute time in 2000 and BMI (parameter estimate, 5.24; standard error, 1.86; P = .006). We also found a significant positive association between food retail gap per capita and BMI when controlling for region and population percentage below poverty (parameter estimate, 0.024; standard error, 0.006; P < .001).
In linear regression models with county-level mean BMI as the dependent variable and difference in commute times from 1990 to 2000 and retail gap per capita as independent variables, the most parsimonious model included the covariates population percentage below poverty and regional fixed effects and explained 43% of variance in county-level BMI. When 2000 commute time and food retail gap per capita were included as independent variables, controlling for diversity index and percentage below the poverty level, the model explained 40% of variance in county-level BMI.
The point estimates for each of the county-level variables of interest (commute time and retail gap per capita) are presented for 4 model specifications (Table 4). In Model 1, we did not include any additional covariates. Individual covariates were added in models 2 and 3. In model 4, regional fixed effects were added. All 5 measures of county-level commute time and food retail gap per capita were positively associated with individual-level BMI. These effects were significant in the unadjusted model (model 1), and the significance remained when individual-level and regional covariates were included in models 2, 3, and 4, with the exception of average commute time increase in model 4. When 2000 commute time and retail gap per capita were both included in the same model with individual-level and regional covariates, the parameter estimates for the county-level variables of interest were no longer significant.
Our results demonstrate a positive correlation between percentage of rural residents and 1) commute times and 2) food retail gap per capita, suggesting that counties with a higher percentage of rural residents have longer commute times and greater retail shortfalls, and thus residents may generally spend food dollars outside their county of residence. Previous studies have found positive associations between BMI and travel distance to grocery stores (23) and time spent in cars (21,22).
We found significant cross-sectional correlations between individual-level and county-level BMI and 1) commute times and 2) food retail gap per capita, but significance did not remain when both were included in the individual-level model. This attenuation could be due to model over-adjustment if commute time and retail gap are both on the causal pathway explaining the relationship between rural residence and BMI.
These analyses support strategies presented in Recommended Community Strategies and Measurements to Prevent Obesity in the United States (30) to improve geographic availability of supermarkets in underserved areas and provide incentives to food retailers to offer healthier food and beverage choices in underserved areas. If implemented, these strategies would decrease travel times necessary for accessing healthy, affordable foods among low-income and rural residents. When combined with health education efforts and mass media campaigns encouraging healthy food choices, more accessible and affordable healthy foods may lead to healthier food consumption patterns and to lower obesity prevalence in these groups.
In a qualitative study of rural Georgia adults, participants identified several barriers to obtaining healthy foods, including poor selection, limited time, fuel prices, and the distance (15-45 miles) to larger communities with bigger stores and better selection (31). Another study found that longer distance traveled to the primary grocery store was associated with higher BMI (23). This previous work, taken together with our results, supports the notion that rural residents who travel farther to shop for food may purchase less healthful food. However, we did not measure the distance to the locations where people shopped and assumed that a positive food retail gap indicated a general trend for rural residents to shop for food outside their county of residence. Future work should assess the relationship between commute times and the locations where they purchase food. Future work should also include mediational analyses to examine the relationships between commute time, food shopping frequency and location, diet quality, and BMI.
This study has several limitations. Foremost is the ecological design, which used several different data sources. The inconsistent timing of data collection for commute times (1990, 2000), food retail gaps (2008), and BMI (2003-2007) is an additional limitation. However, we used the most recent data available, and average commute time is a proxy for distance between place of employment and residence (32). A related limitation is the exclusion of people in the 36 counties where BRFSS did not provide county-level identifiers, pointing to the need for more work to examine rural populations. An additional caveat is that we used self-reported height and weight from BRFSS to calculate BMI, potentially biasing results toward the null if hypothesized relationships between commute times, food retail gaps, and BMI truly exist, because of potential underestimation of weight status. The use of a commercial business database (InfoUSA) to obtain sales data is also a limitation, because such databases may contain errors (33). Finally, in these analyses, we assumed commute time referred to time spent driving. Some people may walk or bike to work instead of drive; however, few Americans actively commute (34).
This study is the first to examine correlations between commute times, food retail gap per capita, and mean BMI in counties in North Carolina. We present an approach to studying the association between BMI and variables related to the built and economic environments, providing support for the notion that economic and built environment factors are related to obesity.
We thank Tammy Lister, Economic Development Intelligence System, North Carolina Department of Commerce, and Catherine Spisszak, ESRI, for their helpful comments and guidance. We thank Tom Hoerger from the RTI-UNC Center of Excellence in Health Promotion Economics for his guidance. We thank Satomi Imai, Jared McGuirt, Scott Ickes, and Lloyd Novick. The first author was supported by a faculty startup award from East Carolina University Division of Research and Graduate Studies during the time this work was completed.
Corresponding Author: Stephanie B. Jilcott, PhD, Department of Public Health, East Carolina University, 600 Moye Blvd, MS 660, Greenville, NC 27834. Telephone: 252-744-4034. E-mail: firstname.lastname@example.org.
Author Affiliations: Haiyong Liu, Justin B. Moore, Jeffrey W. Bethel, East Carolina University, Greenville, North Carolina; James Wilson, Northern Illinois University, DeKalb, Illinois; Alice S. Ammerman, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina.
- Lee SK, Sobal J, Frongillo EA, Olson CM, Wolfe WS. Parity and body weight in the United States: differences by race and size of place of residence. Obes Res 2005;13(7):1263-9.
- Dalbey M. Implementing smart growth strategies in rural America: development patterns that support public health goals. J Public Health Manag Pract 2008;14(3):238-43.
- McMurray RG, Harrell JS, Bangdiwala SI, Deng S. Cardiovascular disease risk factors and obesity of rural and urban elementary school children. J Rural Health 1999;15(4):365-74.
- Blankenau J. Nutrition, physical activity, and obesity in rural America: a series examining health care issues in rural America. Wayne State College Center for Rural Affairs; 2009. http://files.cfra.org/pdf/Nutrition-Physical-Activity-and-Obesity-in-Rural-America.pdf. Accessed April 26, 2010.
- Jackson JE, Doescher MP, Jerant AF, Hart LG. A national study of obesity prevalence and trends by type of rural county. J Rural Health 2005;21(2):140-8.
- Morland K, Diez Roux AV, Wing S. Supermarkets, other food stores, and obesity: the Atherosclerosis Risk in Communities study. Am J Prev Med 2006;30(4):333-9.
- Bustillos B, Sharkey JR, Anding J, McIntosh A. Availability of more healthful food alternatives in traditional, convenience, and nontraditional types of food stores in two rural Texas counties. J Am Diet Assoc 2009;109(5):883-9.
- Kaufman P. Rural poor have less access to supermarkets, large grocery stores. Rural Development Perspectives 1999;13(3):19-26.
- Liese AD, Weis KE, Pluto D, Smith E, Lawson A. Food store types, availability, and cost of foods in a rural environment. J Am Diet Assoc 2007;107(11):1916-23.
- Sharkey JR, Horel S. Neighborhood socioeconomic deprivation and minority composition are associated with better potential spatial access to the ground-truthed food environment in a large rural area. J Nutr 2008;138(3):620-7.
- Jilcott SB, Laraia BL, Evenson KR, Ammerman AS. Perceptions of the community food environment and related influences on food choice: a qualitative analysis. Women Health 2009;49(2-3):164-80.
- Robinson NC, Caraher M, Lang T. Access to shops: the views of low-income shoppers. Health Education J 2000;59(2):121-36.
- Guy C, Clarke G, Eyre H. Food retail change and the growth of food deserts: a case study of Cardiff. International Journal of Retail and Distribution Management 2004;32(2):72-88.
- Furey S, Strugnell C, McIlveen H. An investigation of the potential existence of ‘food deserts’ in rural and urban areas of Northern Ireland. Agriculture and Human Values 2001;18(4):447-57.
- Hallsworth A, Worthington S. Local resistance to larger retailers: the example of market towns and the food superstore in the UK. International Journal of Retail and Distribution Management 2000;28(4-5):207-16.
- Yoo S, Baranowski T, Missaghian M, Baranowski J, Cullen K, Fisher JO, et al. Food-purchasing patterns for home: a grocery store-intercept survey. Public Health Nutr 2006;9(3):384-93.
- Bhargava A. Socio-economic and behavioural factors are predictors of food use in the National Food Stamp Program Survey. Br J Nutr 2004;92(3):497-506.
- Lopez R. Urban sprawl and risk for being overweight or obese. Am J Public Health 2004;94(9):1574-9.
- Lee IM, Ewing R, Sesso HD. The built environment and physical activity levels: the Harvard Alumni Health Study. Am J Prev Med 2009;37(4):293-8.
- Pendola R, Gen S. BMI, auto use, and the urban environment in San Francisco. Health Place 2007;13(2):551-6.
- Frank LD, Andresen MA, Schmid TL. Obesity relationships with community design, physical activity, and time spent in cars. Am J Prev Med 2004;27(2):87-96.
- Lopez-Zetina J, Lee H, Friis R. The link between obesity and the built environment. Evidence from an ecological analysis of obesity and vehicle miles of travel in California. Health Place 2006;12(4):656-64.
- Inagami S, Cohen DA, Finch BK, Asch SM. You are where you shop: grocery store locations, weight, and neighborhoods. Am J Prev Med 2006;31(1):10-17.
- Hall SA, Kaufman JS, Ricketts TC. Defining urban and rural areas in US epidemiologic studies. J Urban Health 2006;83(2):162-75.
- Morland K, Wing S, Diez Roux A, Poole C. Neighborhood characteristics associated with the location of food stores and food service places. Am J Prev Med 2002;22(2):23-9.
- Block JP, Scribner RA, DeSalvo KB. Fast food, race/ethnicity, and income: a geographic analysis. Am J Prev Med 2004;27(3):211-7.
- Hero R, Tolbert C. A racial/ethnic diversity interpretation of politics and policy in the states of the U.S. American Journal of Political Science 1996;40(3):851-71.
- Lin FQ. Solving multicollinearity in the process of fitting regression model using the nested estimate procedure. Quality and Quantity 2008;42:417-26.
- Leyland AH, Groenewegen PP. Multilevel modelling and public health policy. Scand J Public Health 2003;31(4):267-74.
- Centers for Disease Control and Prevention. Recommended community strategies and measurements to prevent obesity in the United States. MMWR Recomm Rep 2009;58(RR-7):1-26.
- Kegler MC, Escoffery C, Alcantara I, Ballard D, Glanz K. A qualitative examination of home and neighborhood environments for obesity prevention in rural adults. Int J Behav Nutr Phys Act 2008;5:65.
- Guiliano G, Small KA. Is the journey to work explained by urban structure? Urban Studies 1993;30(9):1485-500.
- Jilcott SB, McGuirt J, Imai S, Evenson KR. Measuring the retail food environment in rural and urban North Carolina counties. J Public Health Manag Pract 2010. Forthcoming.
- French SA, Story M, Jeffrey RW. Environmental influences on eating and physical activity. Annu Rev Public Health 2001;22:309-35.