In the USA, poor diet quality constituted the leading risk factor for mortality(Reference Mokdad, Ballestros and Echko1) with wide and persistent disparities in nutrition and related health outcomes(Reference Mensah, Brown and Pratt2). National surveys have identified differences across racial/ethnic groups in meeting fruit and vegetable intake as well as specific dietary components including energy intake, dietary fibre, added sugar and Na(Reference O’neil, Nicklas and Keast3). Gaps in diet quality also exist between lower-income and higher-income individuals and have worsened over time(Reference Wang, Leung and Li4). These trends parallel disparities in risk factors for and rates of chronic diseases including CVD and cancer; during recent decades, disparities between low- and high-socio-economic status groups and between racial/ethnic groups have either not improved or worsened(Reference Beckman, Herrin and Nasir5–Reference Withrow, Berrington de González and Spillane10). Between 1980 and 2014, these differences translated into large and growing inequalities in life expectancy among US counties, with variation explained by differences in socio-economic and race/ethnicity factors or behavioural and metabolic risk factors(Reference Dwyer-Lindgren, Bertozzi-Villa and Stubbs11).
Research in recent decades has shed light on non-individual factors that create and reinforce nutrition-related disparities(Reference Northridge, Sclar and Biswas12). Theory and evidence have pointed to supportive food environments, including those with close, easy and affordable access to healthy food options, as critical for nutrition outcomes(Reference Traverso-Yepez and Hunter13,Reference Caspi, Sorensen and Subramanian14) . Similar work has highlighted neighbourhood obstacles to healthy diets, including food deserts (i.e. neighbourhoods with limited access to supermarkets or other sources of healthy and affordable food(15)) and food swamps (i.e. areas with high density of establishments selling high-calorie fast food and junk food, relative to healthier options(Reference Cooksey-Stowers, Schwartz and Brownell16)). Notably, multiple studies have demonstrated patterning of retail food sources by neighbourhood socio-economic and racial/ethnic composition. Substantial evidence shows increases in fast food or other unhealthy options or absence of healthy stores and supermarkets for areas with lower socio-economic status and higher proportions of non-White residents(Reference Lamichhane, Warren and Puett17–Reference Raja, Ma and Yadav22). Combined, these point to underlying historic and persistent structural inequalities across neighbourhoods that show up in many dimensions of healthy amenities including retail food access. As a result, policies aimed to address retail food access including financing for fresh food options, shifting items available in smaller stores and increasing transportation options(Reference Chrisinger23–Reference Ulmer, Rathert and Rose26) target these high-need areas.
Almost no work has analysed trends in retail food environments(Reference Sanchez-Vaznaugh, Weverka and Matsuzaki19,Reference Berger, Kaufman and Bader27–Reference Widener and Shannon29) , with little examination at the national level. A national study across urban areas in New Zealand found a decrease in distance and time to both fast-food outlets and supermarkets, with the biggest decrease in distance for supermarkets seen in the most deprived areas(Reference Hobbs, Mackenbach and Wiki28). However, these patterns may not translate to other settings, including the USA. A study of fast-food and full-service restaurants across the USA found higher proportions of fast-food restaurants out of total restaurants in low-income and predominantly non-Hispanic Black neighbourhoods(Reference Powell, Chaloupka and Bao20). Yet, this work was cross-sectional and limited to unhealthy food outlets. We found no national research on temporal trends in the equitable distribution of food outlet resources; previous work on disparities in food access has centred around single or multi-city datasets(Reference Lamichhane, Warren and Puett17–Reference Powell, Chaloupka and Bao20).
Therefore, we aim to explore national disparities in retail food environment by neighbourhood composition (racial/ethnic and income). Using a national, longitudinal dataset spanning 25 years (1990–2014), we describe changes to disparities in built food environments (retail outlets) across time and space. We hypothesised that national patterns of food access will match single- and multi-city studies (greater access for neighbourhoods with higher income and non-Hispanic White residents). We further hypothesised that retail food environment disparities would mirror increased disparities in nutrition and health outcomes and would differ across geographic region and urbanicity.
Methods
Study sample
We compiled retail environment, demographic and socio-economic data between 1990 and 2014 for all census tracts (‘tracts’) in the contiguous USA (n 72 538). Tracts were excluded if they contained no land area (i.e. water tracts) or were non-residential (i.e. < 50 residents), leaving 71 547 tracts. For consistency despite boundary changes, all measures were assigned to 2010 census geographies. Sensitivity analyses examined associations only for tracts within Metropolitan Statistical Areas (MSA), and urban areas classified using Rural Urban Commuting Area (RUCA) codes (1; 1 or 2).
National food environment data and variables
Longitudinal food environment retail business data came from the National Establishment Time Series (NETS) database, licensed from Walls & Associates (Denver, CO). Detailed methods on the cleaning and classification of the NETS data can be found elsewhere(Reference Hirsch, Moore and Cahill30). Briefly, NETS represents a census of all US businesses and their industrial classification using Standard Industrial Classification (SIC) codes. It is considered one of the most comprehensive databases of US establishments available, with prior studies leveraging NETS to examine food environment(Reference Berger, Kaufman and Bader27,Reference Kaufman, Sheehan and Rundle31) . From NETS, records were categorised first using SIC codes (online Supplementary Table S1). To fully capture national chain stores, we searched for company or trade names that were on the Nielsen (New York) TDLinx list, Technomics/R&I chain name list, or trade channel matching corresponding categories.
The category ‘All Food Stores’ included any location where food can be purchased for off-premise consumption. ‘Healthy Food Stores’ (HF), a subcategory of All Food Stores, was composed of food stores offering a larger selection of healthy food items. We used stores with consistent food offerings across time and space. HF contained the subcategories ‘All Supermarkets’ and ‘Fruit and Vegetable Markets’ where supermarkets represent stores offering a wide variety of food and household items that typically have a wider selection than traditional grocery stores and were defined using SIC codes, sales amounts, employee counts and chain names. We defined ‘Unhealthy Food Sources’ (UF) as both food stores and restaurants (distinguished from stores as selling ready-to-eat items, including those for on-premise consumption) offering primarily unhealthy food items or limited supplies of healthy foods. This included ‘Convenience Stores’, ‘Fast Food Restaurants’, and ‘Bakeries, Candy, and Ice Cream Stores’.
NETS addresses were geocoded for each year a business was open between 1990 and 2014. We calculated annual counts of food establishments per tract and densities by dividing counts by land area (km2 with water bodies, from census TIGER files, removed).
Sociodemographic data
The Longitudinal Tract Database product (LTDB)(Reference Logan, Xu and Stults32) was used to estimate population density, percent of population in poverty, unemployed, aged 60+ years, aged 75+ years, with high school or college educations, who own their home, recent immigrant (within 10 years), foreign born, and persons who are not proficient in speaking English.
LTDB data were also used for our exposures of interest: racial/ethnic composition (based on percent within different race/ethnic categories) and inflation-adjusted median household income. For consistency with previous work(Reference Moore and Diez Roux33,34) , we defined a tracts as predominantly one racial/ethnic group if greater than 60 % of residents were of that particular racial/ethnic group (i.e. predominantly non-Hispanic White; predominantly non-Hispanic Black; predominantly Hispanic; predominantly Asian). Tracts that did not fall into any of these categories were classified as racially mixed areas. We also created longitudinal racial composition metrics defined using categories from previous work(34). We used inflation-adjusted income tertiles across all years (1990, 2000 and 2010) to examine socio-economic disparities (low – less than or equal to $43 135, middle – between $43 135 and less than or equal to $61 883, and high – greater than $61 883).
Statistical analyses and mapping
We calculated means, standard deviations, and frequencies for food densities and sociodemographic characteristics within each year and over time (1990–2014). We mapped densities of retail food environment variables in each year and for change between 1990 and 2014.
Tract and food environment characteristics were compared across categories of racial/ethnic composition. Because food environment variables are counts (number of stores within each tract), we used generalised Poisson regression to examine associations of racial/ethnic composition and income tertiles with number and trajectory of establishments. The count of each category was modelled as a function of time (in 10-year increments), racial/ethnic composition (time-varying), population density (time-varying z-score using mean and sd of 1990–2010 pooled data) and land area as an offset. Models estimated racial differences at baseline (1990), at 25 years (2014) and in linear time trend across the years (1990–2014). Analyses were clustered by tract to account for repeated measures. Exponentiated coefficients for the time trend (modelled yearly) can be interpreted as mean or expected change in stores per land area over 10 years, adjusted for population density. We did not isolate the effects of racial/ethnic composition from socio-economic composition due to strong correlations between the two; instead, we show similarly constructed models for both racial/ethnic composition and income categories.
We performed all analyses in ArcGIS (10.8) and SAS (9.4). Sensitivity analyses examined associations within restricted geographies (see above).
Results
Sample characteristics (1990–2014)
We note trends of increasing population density, older adult population and percent with at least a 4-year college education (Table 1). Tracts became more diverse, with increases in percent foreign born, recent immigrants, not proficient in English and non-White races/ethnicities. This pattern was also noted in racial composition categories; while a majority of tracts were predominantly non-Hispanic White in 1990 (79·5 %), far fewer were predominantly non-Hispanic White by 2010 (63·9 %), with 11 751 (16·4 %) changing from predominantly non-Hispanic White between 1990 and 2010 (online Supplementary Table S2). In all years, tracts that were predominantly non-Hispanic Black, Hispanic, and Asian experienced higher population density, poverty, and unemployment as well as lower educational attainment, household income, and owner-occupied housing (online Supplementary Table S3).
* Census tracts were excluded if they contained no land area (i.e. were water tracts) or had less than fifty residents.
† Harmonised data for the 1990, 2000 and 2010 population census, published in the Longitudinal Tract Database product (LTDB)(Reference Logan, Xu and Stults32), were used to estimate population density. LTDB data from the 1990 and 2000 population census and the 2008–2012 American Community Survey were used to estimate all other values.
‡ We defined a tract as predominantly one racial/ethnic group if greater than 60 % of residents were any particular racial/ethnic group. Tracts that did not fall into any of these categories were classified as racially mixed areas.
Change in retail food environment (1990–2014)
The density of All Food Stores increased between 1990 (mean 2·35, sd 7·07 stores/km2) and 2014 (mean 3·64, sd 11·41 stores/km2) (online Supplementary Table S4). Across all years, there were higher densities of UF than HF (Fig. 1). Notably, UF density more than doubled over the observed time period (mean 1·93 stores/km2 in 1990 v. mean 4·27 stores/km2 in 2014; average change of 2·35 stores/km2 between the two times), while HF density remained relatively flat (mean 0·25 stores/km2 in 1990 v. mean 0·27 stores/km2 in 2014; average change of 0·02 stores/km2 between the two times). These patterns were retained within most subcategories, with UF (i.e. fast food, convenience stores, bakeries, candy and ice cream) increasing, while supermarkets remained flat or slightly decreased (average change of -0·01 stores/km2 between the two time periods). Fruit and vegetable markets were an exception, increasing from an average of 0·06 stores/km2 in 1990 to 0·09 stores/km2 in 2014 (average increase 0·02 stores/km2 between the two times). Examining change by decade, we see slightly larger changes between 2000 and 2010 than 1990 and 2000 (online Supplementary Table S4). Graphing total counts for these retail categories within the USA reveals a steeper rise in UF (red line) compared with the relatively flat HF category (blue line) and a slight dip in all categories around 2009–2011.
When displaying annual densities geographically (not shown), higher densities of all food business categories were generally located in dense population centres. Maps of change between 1990 and 2014 similarly show larger increases in urban areas with most of the larger, rural tracts experiencing relatively little change in either direction (Fig. 2(a)–(c)). As shown in the inlaid maps of the four largest cities, most larger density changes (both increases and decreases) occurred in urban tracts. Additionally, maps illustrate larger increases in unhealthy food (Fig. 2(b) green) compared with the relative stability of healthy food (Fig. 2(c) white). Patterning of change by region of the USA illustrated more stability (lack of change) in the mountain west.
Retail food environment measures by racial/ethnic composition and socio-economic status
In both 1990 and 2014, tracts that were predominantly non-Hispanic White in 1990 had lower food store density than tracts that were predominantly non-Hispanic Black, Hispanic, Asian or mixed (Table 2). Predominantly non-Hispanic Black, Hispanic, Asian or mixed tracts in 1990 also experienced greater increases in food store density between 1990 and 2014. However, this was primarily driven by UF (increases for all tracts). HF densities declined, on average, between 1990 and 2014 for tracts that were predominantly non-Hispanic Black and predominantly Asian in 1990, with all tracts experiencing declines in supermarket density except predominantly non-Hispanic White tracts (which experienced no change).
* Food environment metrics derived from National Establishment Time Series (NETS) data 1990–2014 annually. Categories are not mutually exclusive. All Food Stores includes the other store categories (not the retail for consumption on site). Healthy Food Stores includes all Supermarkets and Fruit and Vegetable markets. All Unhealthy Food sources includes Fast Food, Convenience Stores, and the Bakery categories.
† Harmonised data for the 1990 census, 2000 census and the 2008–2012 American Community Survey, published in the Longitudinal Tract Database product (LTDB)(Reference Logan, Xu and Stults32), were used to estimate all percentage of each race/ethnicity within each tract. We defined a tract as predominantly one racial/ethnic group if greater than 60 % of residents were any particular racial/ethnic group. Tracts that did not fall into any of these categories were classified as racially mixed areas.
‡ Census tracts were excluded if they contained no land area (i.e. were water tracts) or had less than fifty residents.
From multi-level Poisson models, in 1990, predominantly non-Hispanic Black, Hispanic, Asian and mixed racial composition tracts had more food stores than predominantly non-Hispanic White tracts (Table 3). These differences persisted and slightly increased, by 2014 for all tracts aside from predominantly Asian tracts (where the disparity persisted but narrowed). Estimated time trends 1990–2014 showed that almost all tracts, except predominantly Asian, experienced increases in All Food Stores. However, only predominantly non-Hispanic Black, predominantly Hispanic and mixed racial composition tracts experienced increases in HF (1·06 CI 1·01, 1·11; 1·30, CI 1·17, 1·43; and 1·40, CI 1·33, 1·46, respectively). Predominantly Asian tracts experienced a decrease (0·71, CI 0·58, 0·88), and predominantly non-Hispanic White tracts experienced no statistically significant time trend (1·00, CI 0·99, 1·01). All tracts experienced increases in UF, although in predominantly Asian tracts this trend failed to reach statistical significance. When interpreting estimates and CI for predominantly Asian tracts, small sample sizes (n 30 in 1990, n 118 in 2000, and n 250 in 2010) should be considered.
* Food environment metrics derived from National Establishment Time Series (NETS) data 1990–2014 annually. Categories are not mutually exclusive. All Food Stores includes the other store categories (not the retail for consumption on site). Healthy Food Stores includes all Supermarkets and Fruit and Vegetable Markets. All Unhealthy Food sources includes Fast Food, Convenience Stores and the Bakery categories.
† Harmonised data for the 1990 census, 2000 census and the 2008–2012 American Community Survey, published in the Longitudinal Tract Database product (LTDB), were used to estimate all percentage of each race/ethnicity within each tract. We defined a tract as predominantly one racial/ethnic group if greater than 60 % of residents were any particular racial/ethnic group. Tracts that did not fall into any of these categories were classified as racially mixed areas.
‡ Census tracts were excluded if they contained no land area (i.e. were water tracts) or had less than fifty residents.
§ Estimates derived from generalised Poisson regression models modelled as a function of time (in years), the census tract racial/ethnic composition, the population density and the land area of the tract as an offset.
In 1990, compared with high-income tracts, low-income tracts had approximately 81 % fewer food stores (all types) per land area, controlling for population (0·19, CI 0·16, 0·21) and middle income had 72 % fewer (0·28, CI 0·25, 0·31) (Table 4). Although the gradient remained across all time periods, these disparities narrowed by 2014 (0·45 CI 0·40, 0·51 and 0·33, CI 0·29, 0·38, for low and middle income, respectively). This was driven by increasing time trends for low- and middle-income tracts and a stable time trend for high-income tracts. HF decreased for middle- and high-income tracts but increased for low-income tracts. While all tracts experienced increasing time trends for UF, the increases were greater for low- and middle-income than high-income tracts.
* Food environment metrics derived from National Establishment Time Series (NETS) data 1990–2014 annually. Categories are not mutually exclusive. All Food Stores includes the other store categories (not the retail for consumption on site). Healthy Food Stores includes all Supermarkets and Fruit and Vegetable Markets. All Unhealthy Food Sources includes Fast Food, Convenience Stores and the Bakery categories.
† Harmonised data for the 1990 census, 2000 census and the 2008–2012 American Community Survey, published in the Longitudinal Tract Database product (LTDB), were used to estimate inflation-adjusted, median household income within each tract. Categories represent tertiles.
‡ Census tracts were excluded if they contained no land area (i.e. were water tracts) or had less than fifty residents.
§ Estimates derived from generalised Poisson regression models modelled as a function of time (in years), the census tract racial/ethnic composition, the population density and the land area of the tract as an offset.
Sensitivity analyses showed similar results for urban areas (measured using MSA or RUCA) but with larger increases over time (online Supplementary Tables S5–S10).
Discussion
During the period 1990–2014, the prevalence of food stores increased across the USA. Importantly, at every time point, there were higher densities of unhealthy than healthy food retail and steep increases in availability of UF, while HF remained stable. Contrary to expectations, predominantly non-Hispanic Black, Hispanic, Asian and mixed tracts (unadjusted for income) had higher food store densities in 1990 and maintained a denser food environment through 2014 than predominantly non-Hispanic White tracts. In estimated time trends, we see that almost all tracts experienced increases in food stores but that only predominantly non-Hispanic Black, Hispanic and mixed racial composition tracts experienced increases in HF. Across all time periods, low- and medium-income tracts had lower food store density, but the disparity narrowed between 1990 and 2014 (unadjusted for race). Low-income tracts even experienced increases in HF density. Combined, our findings illustrate a national food landscape with both persistent and shifting resources in neighbourhoods with different racial/ethnic and socio-economic compositions.
Reconciling our contrary findings on racial disparities with prior, geographically restricted literature(Reference Walker, Keane and Burke21,Reference Raja, Ma and Yadav22,Reference Moore and Diez Roux33) is challenging. One explanation is a potential conflation of urbanicity and racial/ethnic composition. Specifically, that disparities observed prior in dense, urban settings in the northeast region may not be consistent for patterns across the entire country. This is especially true, since a large portion of census tracts in the USA are rural and predominantly non-Hispanic White. Indeed, our findings match patterns of walkability nationally. In a national study of five walkability indicators in almost 65 000 tracts, neighbourhoods with higher proportion non-Hispanic Black were more walkable (i.e. shorter block length, greater street node density, more developed land use and higher density of street segments)(Reference King and Clarke35). These walkability metrics represent more intensely developed neighbourhoods, which may have higher density of establishments, including all types of food stores. Due to historic and present-day discriminatory housing policies and practices, predominantly non-Hispanic Black, Hispanic, Asian and mixed neighbourhoods exist primarily in intensely developed, urban areas. It is important to consider rural, predominantly non-Hispanic White areas as high need for food access interventions. This rural–urban racial pattern is not isolated to our work; a study of change in fast-food establishments near public schools in California also found differences by urban and rural localities(Reference Sanchez-Vaznaugh, Weverka and Matsuzaki19). Additionally, results from South Carolina (a largely rural state) show lower minority populations in places with supermarkets and fast-food outlets(Reference Lamichhane, Warren and Puett17). Interestingly, when we restricted to the year and limited urban geographies of a prior study with similar objectives(Reference Moore and Diez Roux33), we find similar patterns of racial/ethnic disparities (less stores for non-White neighbourhoods) found in that paper (not shown). This implies that differences in disparity patterns are not the result of data differences and instead may be time and geographic scale-dependent. Nonetheless, when we restricted analyses to urban areas (defined in multiple ways), the resulting pattern of results was consistent with our national results. This indicates there may be a complicated story across geographies, urbanicity and race. Future work could explore spatial patterns of disparities using geographic weighted regression(Reference Chi, Grigsby-Toussaint and Bradford36) or stratification by more granular urbanicity and population density measures.
Notably, non-White neighbourhoods had more of both UF and HF. As posited by others, spatial co-occurrence and clustering of establishments may account for some racial/ethnic patterning(Reference Lamichhane, Warren and Puett17,Reference Yang, Wang and Qiu37) . The nuance of results for healthy and unhealthy food establishments shed light on some of the disparities seen previously, as intensely developed spaces (often with more racial/ethnic diversity) had higher UF and larger increases over time. This is consistent with research on ‘food swamps’, fast-food restaurants, or other unhealthy food options(Reference Fraser, Edwards and Cade18,Reference Powell, Chaloupka and Bao20,Reference Yang, Wang and Qiu37,Reference Goodman, Thomson and Landry38) . In particular, our results align with the national, cross-sectional evidence showing greater proportion fast-food restaurants in low-income and predominantly black neighbourhoods(Reference Powell, Chaloupka and Bao20). Our results for predominantly Hispanic and Asian tracts suggest that neighbourhoods may be ethnic enclaves with population-specific specialty stores including food stores with culturally appropriate options(Reference Osypuk, Roux and Hadley39).
The negative gradient we observed between income and food stores is consistent with many cross-sectional studies(Reference Fraser, Edwards and Cade18,Reference Powell, Chaloupka and Bao20,Reference Walker, Keane and Burke21,Reference Moore and Diez Roux33) as well as one, state-level longitudinal analyses of fast-food restaurants(Reference Sanchez-Vaznaugh, Weverka and Matsuzaki19). Most interestingly, we found a narrowing of income disparities across the 25 years included in this work. This corroborates the national study in New Zealand that found a larger decrease in distance for to supermarkets within over time for the most deprived areas(Reference Hobbs, Mackenbach and Wiki28). During this period, several national- and local-level policies attempted to address disparities in food environments(Reference Chrisinger23–Reference Ulmer, Rathert and Rose26). Specifically, the Obama administration (2009–2017) launched the Healthy Food Financing Initiative, modelled after several similar local models(Reference Giang, Karpyn and Laurison24,Reference Holzman25) . Notably, in 2011, USDA created their Food Desert Locator (now the Food Access Research Atlas) which combined information on income and access to a supermarket or large grocery store(15). This tool helped governments to channel resources towards increasing healthy food in high-need areas. Our work cannot provide causal evidence around the impact of these policies, and it is unlikely that the number of stores generated by these policies was enough to account for the changes we observed. However, it is encouraging to see the income gap narrowing. Nonetheless, income disparities remain even in the later time periods and UF increased for tracts of all income levels, reinforcing evidence regarding the rise of ‘food swamps’ and increase in unhealthy food options(Reference Lamichhane, Warren and Puett17–Reference Powell, Chaloupka and Bao20,Reference Yang, Wang and Qiu37) . It is also possible that tertiles did not entirely isolate the lowest resourced tracts where food retail access may have remained unchanged; additional work may benefit from examining income across additional categories (e.g. deciles) to see if patterns remain consistent.
Food retail environment disparities partially parallel worsening gaps in diet quality and chronic diseases. Increasing UF, especially in low-income and predominantly non-Hispanic Black, Hispanic, and mixed neighbourhoods, follow worsened disparities in cardiovascular risk and cancer over recent decades(Reference Beckman, Herrin and Nasir5–Reference Withrow, Berrington de González and Spillane10). Given our findings on healthy food advantage for these neighbourhoods, patterns of diet quality and health may be more driven by access to unhealthy options than access to healthy options(Reference Cooksey-Stowers, Schwartz and Brownell16,Reference Yang, Wang and Qiu37,Reference Goodman, Thomson and Landry38) . Translating food environment disparities to health highlights other important factors, including transportation, food pricing and cultural acceptability(Reference Andreyeva, Long and Brownell40,Reference James41) .
We expand existing literature in several significant ways. To our knowledge, this is the most comprehensive examination of built food environments in the USA; we displayed 25-year patterns for many types of food stores, across the entire contiguous USA. We limited error and bias by leveraging multiple linked databases, years of geoprocessing and comprehensive category refinement. We explicitly demonstrated that disparities in retail food environments exist nationally, rather than within a restricted set of geographies. Further, these patterns may be different than results from smaller areas. We also illustrated ways in which food retail changed longitudinally, including differences across geography, racial composition and income. This demonstrates that disparities in retail food environment are dynamic over time and modifiable. Finally, our multi-level modelling strategy accounted for the outcome variable (counts) and correlations within tracts over time.
Despite this paper’s strengths, research gaps remain. Large spatial and temporal scales required administrative datasets, which may not accurately represent on-the-ground reality(Reference Forsyth, Lytle and Van Riper42). Indeed, previous work has shown differences in commercial secondary data sources and stores present during field audits(Reference Fleischhacker, Evenson and Sharkey43–Reference Paquet, Daniel and Kestens46). Nonetheless, ground-truthing stores historically and nationally is logistically impossible. Future work could utilise crowdsourced data or commercial imagery (i.e. Google Street View) to capture and archive a prospective database of these detailed data(Reference Nesse and Airt47). Stores were also classified based on assumed purchasing patterns. Yet, we recognise that individuals may purchase healthy foods at locations classified as ‘unhealthy’ (e.g. convenience stores) or purchase unhealthy foods at ‘healthy’ locations (e.g. supermarkets). Importantly, food environments include additional elements beyond retail or business establishments, including marketing, labelling, prices and nutritional quality of the food supply(Reference Downs, Ahmed and Fanzo48). We were unable to account for purchasing patterns, store size, pricing information, or accessibility of culturally appropriate and acceptable foods(Reference Andreyeva, Long and Brownell40,Reference Peñalvo, Cudhea and Micha49,Reference Krukowski, West and Harvey-Berino50) . As we uncovered here, a complicated story exists for urbanicity, racial composition and food environment subtypes. Additional efforts should leverage advanced modelling techniques, such as latent growth curves(Reference Berger, Kaufman and Bader27), to elucidate neighbourhood typologies longitudinally and nationally. Tools such as geographically weighted regression(Reference Chi, Grigsby-Toussaint and Bradford36) or other advanced spatial methods(Reference Jin and Lu51) could identify whether magnitudes of disparities are patterned spatially. Related, racial and ethnic compositions represent a limited view of the complexities of race as a social construct; we were constrained by census-defined categories and by which racial/ethnic categories had large enough populations, densely located, to create ‘predominant’ tracts. Consideration of nuances around intra-group differences and variations in neighbourhood racial composition are important to ensure research encapsulates the complexity of race, ethnicity and racism within the USA. Similarly, median household incomes did not account for household size which may mask differences in purchasing power relevant to food acquisition. Most critically, due to collinearity, we were unable to disentangle racial composition from income distribution. Therefore, all results should be interpreted as the disparity for each component without adjustment for the other. Future work may benefit from the creation of combined categories of income and racial composition (for racial compositions that have enough tracts) to examine potential effect modification. Finally, while we discuss relationships between policies, food environment disparities and nutritional outcome disparities, we caution against discussing causal pathways linking these three outcomes. Future work should explicitly examine these relationships to best understand how policy shifts can impact food environments and subsequent nutrition.
Conclusion
This research is the first to quantify and display food environment resources nationally across time and to illustrate disparities in change in retail food environments within the USA. Results paint a complicated picture with nuance by racial/ethnic group and type of built food environment measure. Combined, the persistent (and growing) racial/ethnic and socio-economic disparities in unhealthy food access and nutrition-related health outcomes suggest that restricting placement of unhealthy food sources may be more important than support for new, healthy food store locations.
Acknowledgements: The authors thank the Urban Health Collaborative at Drexel University and the Built Environment and Health Research Group at Columbia University for their advice, input and support. Authorship: JH, A.R. and GL conceived of, designed and guided this research. J.Q. and S.M. performed all spatial analyses, overseen by JH, A.R., and GL. Complex methodology for classifying business data was created by GS, A.R., N.B., KA and JH. All analyses were carried out by Y.Z. and KA. Interpretation of results was done by all authors, drafted by JA and KM, and edited by all authors. Ethics of human subject participation: This research does not include human subjects.
Financial support:
This work was supported by the National Institute of Aging (grants 1R01AG049970, 3R01AG049970-04S1 and R01AG072634) and Commonwealth Universal Research Enhancement (C.U.R.E) programme funded by the Pennsylvania Department of Health – 2015 Formula award – SAP #4100072543.
Conflict of interest:
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.
Supplementary material
For supplementary material/s referred to in this article, please visit https://doi.org/10.1017/S1368980023000058