Abstract
Urban greenspace is associated with a wide variety of human health benefits, from improved mental health to reduced violent crime and infant birth outcomes. This study investigates long-term trends in the distribution of greenspace across census tracts in the 260 U.S. cities with populations greater than 100,000 residents based on economic privilege and deprivation. Satellite-acquired Normalized Difference Vegetation Index (NDVI) images were used as a measure of greenspace, and the Index of Concentration at the Extremes (ICE) was used as a measure spatial social polarizations of deprived and priveliged populations. From 2001 to 2019, the most privileged, highest ICE quintile census tracts had both the highest mean NDVI and experienced a significantly greater increase in greenspace than the least privileged tracts, indicating that the disparity in greenspace access is widening. Public greenspace initiatives need to focus on programs that are more equitable across the spectrum of economic strata to reduce urban health disparities and address concerns of environmental justice.
Similar content being viewed by others
Background
Vegetation in cities can take many forms, from private lawns to public parks to street trees and grassy boulevards. Past studies have shown that greenspace reduces exposure to harmful air pollutants (Hwang et al. 2011; Nowak et al. 2006), and even small parks provide important shade and heat sinks to combat the urban heat island effect (Aram et al. 2019; Park et al. 2017). Parks provide increased opportunities for physical activity, and research has suggested the power of greenspace in promoting social cohesion and directly reducing depression and stress (Beyer et al. 2014; James et al. 2015; Jennings and Bamkole 2019; Twohig-Bennett and Jones 2018; Wendelboe-Nelson et al. 2019). Other studies have suggested links between increased urban greenness and higher birth weights (Agay-Shay et al. 2014; Dadvand et al. 2012; McEachan et al. 2016), lower violent crime rates (Kuo and Sullivan 2001; Shepley et al. 2019; Venter et al. 2022), and improved long-term health outcomes.
However, a large body of scholarly work has also shown the vast, systemic inequities in the distribution of urban vegetation. Multiple studies have demonstrated strong, positive correlations between income and tree canopy cover in urban census tracts (Casey et al. 2017; Schwarz et al. 2015). Other research has examined the relationship between the “redlining” practices of the Home Owners’ Loan Corporation (HOLC) in the 1930s and the distribution of greenspace today (Nardone et al. 2021). HOLC lending classifications were a systemically racist practice that “redlined” neighborhoods to indicate poor loan quality, often targeting areas with a high minority or immigrant population. They are still associated with present-day levels of racial segregation, poverty, and income inequality (Mitchell et al. 2018). Historically redlined neighborhoods have reduced present-day greenspace (Nardone et al. 2021), experience significantly hotter temperatures (Hoffman et al. 2020), and increased air pollution from vehicle exhaust (Nardone et al. 2020). Recent work has shown that neighborhoods with a higher concentration of non-Hispanic Black residents and residents under the poverty line have significantly less total greenspace as measured by remotely-sensed satellite imagery (NDVI; see methods) than comparable White neighborhoods (Casey et al. 2017).
These findings present serious problems for notions of equitable urban public health and environmental justice in our cities. People living below the poverty line are already more likely to have chronic conditions such as asthma, diabetes, and heart disease (Osborn et al. 2013; Winkleby et al. 1992). A landmark study found a 15-year gap in life expectancy between the richest 1% and poorest 1% of Americans (Chetty et al. 2016). Access to greenspace is inextricably linked to these problems. Inequities in urban greenspace along lines of income and race perpetuate the poor air quality, excessive ambient noise, and increased susceptibility to summer heat waves that lower income and minority neighborhoods experience. The U.S. Environmental Protection Agency (EPA) defines environmental justice as the “fair treatment of all people regardless of race, color, national origin, or income, with respect to the development, implementation, and enforcement of environmental laws, regulations, and policies” so that all can enjoy “the same degree of protection from environmental and health hazards.” (US EPA 2015). Under this framework, we can understand inequities in urban greenspace as an important issue of environmental justice.
COVID-19 further highlighted the vast racial and socioeconomic health inequities throughout the U.S. and the rest of the world. Urban greenness—vegetation—is one of the most important supports for urban public health but is also one of the most privileged resources in an increasingly dense and gray city landscape. Throughout the COVID-19 pandemic, outdoor spaces became a haven for those wanting to leave their home, exercise, and socialize without putting themselves at increased risk for contracting COVID-19 (Noszczyk et al. 2022; Pouso et al. 2021). Use of hiking trails and public parks increased greatly, with a study in Oslo, Norway finding that it increased by 291% during lockdown compared to a three year average for the period of time (Venter et al. 2020). The refuge service that outdoor spaces played during COVID-19 lockdowns reinforced their critical position in the built city matrix.
With the recent availability of high-resolution, satellite-acquired greenspace data, researchers can answer critical questions about greenspace resource availability across wide geographic areas. One recent study (Casey et al. 2017) examined change in greenspace across major metropolitan areas in the U.S. using the Index of Concentration at the Extremes (ICE) as a measure of spatial social polarizations of deprived and privileged populations and satellite-acquired Normalized Difference Vegetation Index (NDVI) as a proxy for neighborhood greenspace. However, that study only examined greenspace change from 2001 to 2010. We need to be updated with the most recent data to determine if disparities are widening or shrinking. This study examines how the relationship between satellite-acquired NDVI and ICE at the census tract level has changed from 2001 to 2019. We can then understand how disparities in greenspace access have evolved over this time period. Our analysis defines “greenspace” to include all urban vegetation, which gives us a holistic picture of greenspace distribution across a city. While recent research has investigated the present-day distribution of urban greenspace across census tracts of different levels of spatial social polarization, this study provides more useful and interpretable findings because we fix the census tract boundaries of our study. This allows us to understand how greenspace has changed in small, geographically constant units in relation to the spatial social polarization of that unit. These findings will give important feedback to policymakers and city officials, making them better equipped to deal with the challenges cities will face with a changing climate.
Methods
Area of interest
We set our study area to include all U.S. incorporated cities with greater than 100,000 residents. City boundaries and census tracts were drawn according to the U.S. Census Bureau 2018 Incorporated Place and Census Tract boundaries to ensure the study was based on a consistent geographic area. All statistics (ICE, NDVI, demographics) come from these census tract boundaries. While past research has used metropolitan statistical areas with greater than 100,000 residents as the study area (Casey et al. 2017), these metropolitan areas often include the outlying suburban or rural areas around a city center. We would expect the relationship between greenspace and income in these outlying rural areas to be different than that of the urban core. For instance, vegetation might be plentiful but property values significantly lower because of the lengthy commuting time to the city center. Other studies have chosen to focus their analyses within city boundaries (Barbosa et al. 2007; Gerrish and Watkins 2018) and take advantage of the dense building and population structures that make greenspace so limited. Therefore, to obtain a more accurate picture of urban greenspace distribution and change we focused our analysis inside major city boundaries. We hypothesize that this will allow for greater differentiation between urban neighborhoods.
U.S. Census demographic data
Demographic data was obtained from the American Community 5-year Survey (ACS) for 2010 and 2019 using the data.census.gov API. Key statistics include population density, and household income quintile thresholds. Controlling for population density would ensure that higher income areas aren’t more green only because there are less buildings and people living there.
We then measure the economic privilege of a census tract using Massey’s Index of Concentration at the Extremes (ICE) (Massey et al. 2001). While other economic indicators, such as the Gini Coefficient and median household income, could also be used to measure neighborhood social polarization, we choose ICE because of its well-established applications and justifications in public health monitoring. It captures extremes of privilege and deprivation, allowing us to explicitly understand the relationship between spatial social polarization and greenspace access (Krieger et al., 2016). We calculate it using ACS data. The ICE for Census Tract i in city j is calculated with the equation (Massey et al. 2001):
where Aij is the number of households in census tract i that have income >80th percentile for city j. Similarly, Pij is the number of households in census tract i that have income under the 20th percentile threshold for city j. Tij is the total number of households in census tract i. ICE ranges from −1 to 1, with values closer to 1 suggesting greater privilege in an area. An ICE value for a census tract can carry more information than a median household income statistic as it helps control for disparities in cost of living across cities. We used ICE measured from 2010 geographic boundaries and demographic data for our analysis. Because census tract boundaries are redrawn every ten years, ICE measured from 2001 data would not be compatible with ICE measured from 2010 to 2019 data. 2010 ICE was chosen since it was the temporal midpoint of the study.
NDVI satellite image acquisition
MODIS NDVI Images at 250 m resolution were obtained from the U.S. Geographical Survey (USGS) AρρEars portal in GeoTIFF format (US Geographical Survey 2022). MODIS satellites capture NDVI images every 16 days at 250m2 resolution, and have been used frequently in studies of urban vegetation (Casey et al. 2017; Gascon et al. 2016; Nardone et al. 2020). Since the U.S. NDVI tends to peak in the month of July, we used the first estimate in that month. Analyzing data at peak NDVI levels allows for more noticeable differences between census tracts. As a sensitivity analysis, an annual average NDVI was computed by obtaining MODIS images and averaging NDVI values from four data points (winter: January 1st; spring: April 7th; summer: July 12th; fall: September 30th). To gain insight into how urban greenspace has changed, we collected data on these dates in 2001—the earliest year NDVI data was available—2010, and 2019. 2019 was the end of our selected study period, and 2010 was a logical midpoint of this window.
After loading the GeoTIFF raster file into R Version 4.0.4 (R Core Team 2022), we used the exactextractr (Baston 2020) package to perform zonal statistics on census tract shapefiles. exactextractr was used to calculate the mean NDVI value per census tract. Raster cells that lie on a census tract boundary are weighted by the fraction of the cell that is within the census tract in question. All other cells are given a weight of “1” or “0” depending on if they lie entirely inside or outside the census tract. Less than 10 percent of census tracts covered less than one NDVI data cell.
LANDSAT, MODIS, and Sentinel-2 satellite systems all offer images that the Normalized Difference Vegetation Index (NDVI) can be extracted from. NDVI is the go-to metric for remotely sensed vegetation data, and is used across agricultural, forestry, and urban planning applications. We select MODIS as our satellite image system of choice for consistency with prior literature (Casey et al. 2017), and because the measurements are more reliable. The daily measurements of MODIS images, which the USGS combines into 16-day averages, are more robust to the effects of cloud cover than LANDSAT images which are measured less frequently but at higher resolution.
The heterogeneity of cities across the United States proved to be a challenge in analyzing NDVI trends. Cities existed in deserts, marine forests, arid plains and many other ecoregions. It is possible that cities in desert climates would have lower overall changes in NDVI than cities in wet, marine forests where vegetation is easier to maintain. While some of these city-specific effects are already addressed by our linear mixed-effects model, we include Omernik’s Level I Ecoregions to account for the ecological differences of each city in our study (Omernik 1987). Level I ecoregions divide the continental U.S. into 12 distinct ecoregions, from the “Great Plains” to the “North American Deserts” and “Wet Tropical Forests” and have been used in prior greenspace research (Casey et al. 2017).
Statistical analyses
We fit Linear Mixed-Effects (LME) models using the lme4 package in R Version 4.0.4 (Bates et al. 2015). LME models are highly effective when analyzing correlated measurements from grouped data. Census tracts within each city are likely to be correlated, since they share an ecological environment, municipal policies, and cultural values about greenspace. Each city spans hundreds of census tracts, and each census tract contains three separate years of NDVI measurements. We built a single-level LME model, fitting random intercepts at the city level to allow for varying levels of baseline average greenspace. We examined percentage greenspace change from 2001 to 2019, with 2001 as our baseline measurement. The percentage change is based on the nine-year intervals of our collected data (NDVI data collected in 2001, 2010, 2019).
Results
This study examined the relationship between census tract Index of Concentration at the Extremes (ICE) and the change in the Normalized Difference Vegetation Index (NDVI) from 2001 to 2019. The study area included 23,998 census tracts across 274 U.S. cities. Table 1 displays the results of Spearman’s correlation coefficient tests for the dependent variables used in our analysis.
Census tracts in the highest ICE quintile had a 7.85% increase in greenspace from 2001–2019 (Table 2). Compared with census tracts in the highest quintile (A), all census tracts had lower percentage gains in greenspace. However, the lowest ICE census tracts did comparatively better (6.23% increase) than the D and C quintile census tracts (5.40% and 5.89% increase respectively). Overall, the lowest (E) ICE quintile census tracts had greenspace growth 20.6% slower than the A quintile and the D quintile census tracts had greenspace growth 31.2% slower.
This trend of greenspace growth follows a distinctive “J” shape illustrated in Fig. 1, where the neighborhoods with the least gains in greenspace fall within the lower-middle ICE quintile. Scatterplots of the relationship between ICE and greenspace growth in the U.S.’ two largest cities by population, New York City and Los Angeles, further demonstrates this trend.
We fit two other models to describe the relationship between greenspace growth and neighborhood ICE. Table 3 displays the results from a model where neighborhood 2010 ICE was fit as a continuous variable. A one unit increase in 2010 census tract ICE is associated with a 1.98% increase in greenspace growth. For example, the typical census tract with ICE = 0 had a 6.36% increase in mid-July NDVI from 2001 to 2019, whereas a census tract with ICE = 1 would have an NDVI increase of 8.34%.
Table S1 displays the results of dichotomizing census tract ICE into two categories: census tracts in the highest ICE quintile (i.e., quintile “A” from Table 2)—and census tracts in the lower four B-E quintiles. The average census tract in the highest ICE quintile had a 7.85% increase in NDVI from 2001 to 2019, whereas the average census tract in the lower four ICE quintiles had a 6.01% increase in July NDVI. Figure 2 depicts the difference in the distributions of NDVI change between these two groups.
(a) Change in NDVI measured from 2001 to 2010 (2001 and 2010 satellite imagery). (b) Change in NDVI measured from 2010 to 2019 (2010 and 2019 satellite imagery). (c) Change in NDVI measured from 2001 to 2019 (2001 and 2019 satellite imagery). All ICE quintiles were measured from 2010 census tract boundaries and 2010 American Community Survey income statistics. Purple bars represent overlap between the top ICE quintile and the lower four ICE quintiles. Blue bars represent values of NDVI change where the lower four ICE quintiles have greater density, and red bars represent values where the top ICE quintile has greater density.
The relationship of ICE with change in NDVI was statistically significantly different from zero for all models, with a threshold of significance α <0.05. Surprisingly, the population density of a census tract was not significantly associated with the 20-year change in greenspace when controlling for ICE and ecoregion in any model.
Figure 3 presents a spatial representation of these findings, demonstrating the strong correlation between neighborhood 2010 ICE, 2010 NDVI, and net 2001–2019 change in NDVI across the three largest cities in the U.S. (New York City, Los Angeles, and Chicago). While only a small subset of the 260 cities included in our full analysis, we chose these cities because of their size and because we found them to be representative of the patterns in the entire study population.
(a) Chicago, (b) New York, (c) Los Angeles. Darker colors represent higher NDVI (green), higher greenspace growth (blue), and higher ICE (purple). Correlations in color shade across the three panels, especially noticeable in Los Angeles and Chicago, suggest a persistence and growth of the socioeconomic gradient of urban greenspace.
Discussion
There exist large disparities in health outcomes across different socioeconomic classes in the United States. Neighborhood greenspace is one of the many factors that can influence the health outcomes of city inhabitants. It is associated with improved air quality, better self-reported mental health, higher physical activity levels, and decreased rates of violent crime along with other impacts (Kuo and Sullivan 2001; Nowak et al. 2006; Pouso et al. 2021; Shepley et al. 2019; Wendelboe-Nelson et al. 2019). However historically redlined areas, as well as poorer areas more broadly, tend to have less greenspace than their more privileged neighbors (Casey et al. 2017; Nardone et al. 2021). Greenspace is one of the many resources attached to a higher social position in the United States. This leads to better health outcomes and a perpetuation of the socioeconomic health gradient. Therefore, a reduction in the greenspace gap between rich and poor neighborhoods could help reduce health disparities. In this study we demonstrate the persistence of disparities in urban greenspace over the last twenty years, 2001–2019. Furthermore, we uniquely suggest that these disparities may also be growing larger, even after considering seasonal changes.
This analysis demonstrated a significant, non-linear trend in the evolution of city greenspace over the course of an eighteen-year window. Most clearly, the top 20% of urban census tracts ranked by ICE Income quintile experienced the largest growth in their greenspace. This trend persists across a variety of different methods of looking at greenspace growth. It is somewhat surprising that the most privileged census tracts, which already had the most greenspace to begin with, also had the largest growth in greenspace. Since we accounted for Omernik ecoregion in our models and constructed models based on a census tract’s yearly average NDVI, it is unlikely that this trend was purely caused by large-scale geographic factors. The results from this study suggest that increased neighborhood privilege is associated with growth in that neighborhood’s greenspace, a pattern possibly explained through increased rates of park construction and tree planting programs.
These findings point to a consistent association between high-privileged neighborhoods and increased greenspace growth. While they cannot prove a causal relationship, multiple sociological mechanisms can support them. It has been widely established that higher ICE neighborhoods have significantly more greenspace than lower ICE neighborhoods even when controlling for housing tenure and population density. Areas with increased greenspace have significantly higher property values than those with less neighborhood greenspace, potentially excluding lower-income individuals from entering those areas (Conway et al. 2010). Higher-income homeowners or landlords may be more able to shoulder the non-trivial costs of purchasing and maintaining trees for private property (McPherson 1994). The correlation between income and educational attainment has been long established (Rowntree 1902), and so it is reasonable to expect that high-ICE neighborhoods will have reasonably high educational attainment. We also confirmed this with our own sensitivity analysis (Table S2). This gives these neighborhoods a tool that may increase their effectiveness in advocating with city government for more trees or other greenspace initiatives.
While the EPA and other governmental bodies can’t regulate all aspects of urban greenspace—for instance, private lawns contribute to the greenspace of a neighborhood but are maintained by private homeowners—it can promote tree planting initiatives, encourage the development of new public parks, and provide funding for green road verges. It is not surprising that higher income neighborhoods have greater availability of urban greenspace, given the increased control that wealthy homeowners have over the vegetation on their property. Because of this, we argue that cities should invest more resources into increasing greenspace in lower and lower-middle class neighborhoods. This will help account for the disproportionate growth in greenspace in higher-privileged census tracts. It is possible that to achieve true greenspace equity, people in the lowest ICE quintile should have a higher amount of public green space compared to people in the highest ICE quintile to counteract the inequity in the quality of privately owned green spaces.
While the highest ICE quintile neighborhoods had the largest increase in NDVI, areas in the bottom 20% also experienced significant improvements to their greenspace over the study period. Most notably, they tended to do better than census tracts in the “middle class”—neighborhoods in the 20th to 60th percentile. This trend can be visualized in Fig. 1, and deviates from the idea that increased neighborhood privilege will always lead to greater greenspace growth.
A census tract with multiple parks, minimal impervious surface area, and plentiful street trees may find it difficult to continue improving their green space access. Essential roads, buildings, and parking lots place a cap on the amount of greenspace a neighborhood can have. A neighborhood that lacks substantial greenspace may find it easier to make the same magnitude of greenspace change than a census tract with already plentiful parks and trees. Street trees can be planted, and vacant lots can even be converted into small parks. Along this line of reasoning, the presence of city greenspace initiatives, especially those targeting minority, lower income neighborhoods, may explain part of the increase among the “E” ranked census tracts.
However, the census tracts of the lower-middle class are lagging behind in greenspace growth. The 20-40th percentile ICE census tracts performed significantly worse than any other quintile. Despite possessing dramatically less greenspace on average than the most privileged census tracts, these tracts experienced 33% less of an NDVI increase. The lower performance of these neighborhoods follows other research suggesting the decline of the U.S. urban middle class through increased social polarization (Foster and Wolfson 2010) and white flight to the suburbs (Kye 2018). These middle-class census tracts may be neglected as cities focus on improving greenspace and social services for their least fortunate populations
Inequities in access to green space can be regarded as an example of how the Fundamental Cause Theory hypothesizes how a higher social position is linked with an array of resources and better outcomes even when health outcomes or desirable resources change. This accumulation of risk factors forms the foundation for Link and Phelan’s conception of Fundamental Cause Theory (FCT). Now a cornerstone of sociological theory, FCT states that social position should be considered a “fundamental” cause of health outcomes because it remains associated with health even when the intervening mechanisms change (Link and Phelan 1995). As long as there is a disparity in the access to resources that avoid disease, Fundamental Cause Theory states that a socioeconomic health gradient will exist. Policies that have a moderate effect on multiple diseases, or those that explicitly address inequality in fundamental resources, will have a more positive long-run impact than those that target intervening resources. Given urban greenspace’s wide-reaching impacts on human health and well-being, we posit that it is a “fundamental resource” worth addressing.
Urban gentrification can be a consequence of urban renewal and greenspace improvements (Jelks et al. 2021). However, we did not find any evidence of gentrification in this study. If gentrification were a significant factor in the observed greenspace trends, we might observe that a census tract’s change in NDVI is positively correlated with a change in ICE. However, we cannot directly observe this using ACS data since census tracts are redrawn every 10 years. The redrawing of census tracts in 2010 makes it extremely difficult to accurately compute a census tract’s change in ICE between 2001 and 2019. However, this calculation is possible for the years 2010 - 2019 for which ACS data is available and census tracts remain constant. We modeled the relationship between a census tract’s change in NDVI and change in ICE and found a small negative association (Table S3). This suggests that the effects of gentrification from greenspace are likely minimal if they exist.
Study findings in context
In a similar study assessing the distribution of greenspace across U.S. urban areas, Casey et al. examined the change in U.S. greenspace 2001-2010 (Casey et al. 2017). Using mid-July and yearly averaged NDVI as their objective measures, they focused on census tracts in U.S. Metropolitan Statistical Areas (MSAs) with a population above 100,000 residents. Their study focused on race-based disparities, finding that census tracts with a higher percentage of white residents tended to have significantly more greenspace in the first place and experienced a small, statistically significant increase in greenspace over the study period. On the other hand, tracts with a higher percentage of Hispanic residents had less greenspace to begin with and experienced a small decrease in greenspace from 2001 to 2010. However, their models found comparatively weaker effects when comparing percentages of Black/African-American or Asian/Pacific Islander residents in a census tract.
Our research refines the study area to focus on more densely populated urban zones, examining all tracts inside city boundaries with greater than 100,000 residents. While MSAs are more standardized than city boundaries, they also tend to overestimate the size of an urban area due to their definition on a county level (Oliveira et al. 2014; Ribeiro et al. 2019). Focusing on city boundaries allows for an examination of more ‘true’ urban areas, in the heart of the densely populated core, whereas MSAs would include census tracts from adjacent rural areas with potentially different relationships between greenspace and income.
Furthermore, a doubling of the study period to eighteen years allows for a more accurate estimation of city evolution. Improvements to city greenspace are unlikely to happen rapidly, with the construction of new parks a large-scale, expensive process that spans multiple years. For example, an initiative to increase urban greenspace in California, known as the “Urban Greening Program”, was signed into law in 2016. The program grew out of a bill originally signed in 2006, and grant recipients included projects that had been seeking funds for many years before that. The grant report for the Del Amo Community Park, in Carson, CA states that the park’s construction had been a “community priority for over two decades” and references multiple delays that blocked its road to funding and completion (Coastal Conservancy 2017). While the Del Amo Park may have experienced a longer timeline than other city greenery projects, its delays and complications reflect a reality of government funding where grants are not continuously or reliably available. Examining city greenspace change over ten or fewer years—as done in prior research (Casey et al. 2017)—allows for a snapshot of city evolution but not the full process where ideas are sparked, grants are applied for, and projects come to fruition.
Compounding the gradual process of park construction, the effects and benefits of new trees may not be noticeable until 15-20 years after planting. This is especially the case with city tree planting programs, where trees are often planted at a young age. Assuming a 20-year maturity window, a 20-year study period would allow trees planted in the early 2000s to have reached a state of more significant shade and maturity by the end of the study period. The urban tree growth model UrbTree also used a 20-year tree growth window as the foundation of its model (Kramer and Oldengarm 2010). This more long-term view of greenspace evolution allows for a greater chance of capturing the impact of efforts to improve urban greenspace. If a city initiated a tree planting program to combat a lack of urban greenery in the early 2000s, the data used in this study could detect a higher NDVI in the 2019 datapoint due to the growth of more mature trees.
Study strengths and limitations
This research presents an in-depth, holistic examination of greenspace growth across U.S. urban areas. By refining the study area to focus on urban census tracts and expanding the study period to nineteen years, this paper provides a needed update on U.S. greenspace equity. Nonetheless, while able to identify an association between privilege and a larger improvement in neighborhood greenspace, this study cannot make any causal claims about the fundamental drivers of greenspace growth. Greenspace may not depend on privilege and instead have a gentrifying effect, attracting people of higher income levels to live in a neighborhood. Municipal and state policymaking, as well as neighborhood social capital and civic capacity, may also impact the prioritization of greenspace development.
The choice of NDVI as a greenspace metric allows for a relatively straightforward data-collection process, but is unable to capture the quality of the greenspace in a neighborhood. Nardone et al. posits that NDVI measurements in cities with arid climates may not accurately reflect the proximity to natural environments. Past literature also suggests that the quality of greenspace can influence the strength of its positive influence on communities (Wood et al. 2018; Zhang et al. 2017). Quality metrics for greenspace include measures of biodiversity and availability of parks. Combining satellite-acquired data on greenspace prevalence with ground-level measurements on quality of greenspace maintenance is a challenging but important area for future research.
NDVI is also prone to measurement error, especially when significant cloud cover obscures an area of interest. While selecting four data points over the course of a year (January 1st, April 7th, July 12th, and September 30th) smooths some of this measurement error, data collection remains a potentially significant form of bias.
Furthermore, the definition of the study area using municipal boundaries will inevitably exclude some surrounding urban neighborhoods that we might wish to count. Municipal boundaries are often decided politically, rather than reflecting meaningful geographical differences. While discussed more extensively in the methods section, we made this choice in order to conservatively exclude more “rural” census tracts. Future research could expand this study area to more surrounding urban areas, such as the City Clustering Algorithm (Oliveira et al. 2014).
Future directions
Studies tracking greenspace growth serve as important accountability measures for cities as they attempt to create a more equitable urban landscape. Research efforts are needed to continue monitoring urban greenspace development, as well as developing better metrics to assess an area’s greenspace quality. Investigations into the sociological explanations behind observed greenspace trends may help city leaders most effectively target limited resources.
The U.S. is an appealing area to study urban greenspace because of its robust population statistics down to the census tract level. However, greenspace inequities extend far beyond the U.S.(Hoffimann et al. 2017; Hwang et al. 2020; Schüle et al. 2019) Further research is needed to understand the dynamics of greenspace growth outside the United States, particularly in fast-growing megacities such as Jakarta, Indonesia and Mumbai, India.
Conclusion
Greenspace distribution remains highly stratified, and these inequities have expanded over the last two decades. However, the non-linearities of this greenspace growth suggest complicated underlying dynamics. The study provides further evidence of the struggles of the American lower-middle class, where the lag in greenspace improvements finds parallels in larger macroeconomic trends.
Such studies are important for urban planners and city government officials in creating more equitable urban landscapes that align with the ideas of environmental justice. By ensuring that neighborhoods have equal access to greenery, and equal resources to improve it, cities can weaken the fundamental link between social position and health. These measures will only become more important as cities become increasingly crowded and hot with the effects of urbanization and global warming. But harnessing the many positive impacts of greenspace can create stronger, healthier neighborhoods for all.
Data availability
NDVI data can be found through this link (https://appeears.earthdatacloud.nasa.gov/), and demographic statistics from the U.S. Census API at data.census.gov. The cleaned data that support the findings of this study, as well as data analysis scripts, are available from the corresponding author upon reasonable request.
References
Agay-Shay K, Peled A, Crespo AV, Peretz C, Amitai Y, Linn S, Friger M, Nieuwenhuijsen MJ (2014) Green spaces and adverse pregnancy outcomes. Occup Environ Med 71:562–569. https://doi.org/10.1136/oemed-2013-101961
Aram F, Higueras García E, Solgi E, Mansournia S (2019) Urban green space cooling effect in cities. Heliyon 5:e01339. https://doi.org/10.1016/j.heliyon.2019.e01339
Barbosa O, Tratalos JA, Armsworth PR, Davies RG, Fuller RA, Johnson P, Gaston KJ (2007) Who benefits from access to green space? A case study from Sheffield, UK. Landsc Urban Plan 83:187–195. https://doi.org/10.1016/j.landurbplan.2007.04.004
Bates D, Mächler M, Bolker B, Walker S (2015) Fitting linear mixed-effects models using lme4. J Stat Softw 67(1):1–48. https://doi.org/10.18637/jss.v067.i01
Beyer KMM, Kaltenbach A, Szabo A, Bogar S, Nieto FJ, Malecki KM (2014) Exposure to neighborhood green space and mental health: evidence from the survey of the health of Wisconsin. Int J Environ Res Public Health 11:3453–3472. https://doi.org/10.3390/ijerph110303453
Casey JA, James P, Cushing L, Jesdale BM, Morello-Frosch R (2017) Race, ethnicity, income concentration and 10-year change in urban greenness in the United States. Int J Environ Res Public Health 14:1546. https://doi.org/10.3390/ijerph14121546
Chetty R, Stepner M, Abraham S, Lin S, Scuderi B, Turner N, Bergeron A, Cutler D (2016) The association between income and life expectancy in the United States, 2001–2014. JAMA 315:1750–1766. https://doi.org/10.1001/jama.2016.4226
Coastal Conservancy, 2017. Del Amo Neighborhood Park
Conway D, Li CQ, Wolch J, Kahle C, Jerrett M (2010) A spatial autocorrelation approach for examining the effects of urban greenspace on residential property values. J Real Estate Finance Econ 41:150–169. https://doi.org/10.1007/s11146-008-9159-6
Dadvand P, de Nazelle A, Figueras F, Basagaña X, Su J, Amoly E, Jerrett M, Vrijheid M, Sunyer J, Nieuwenhuijsen MJ (2012) Green space, health inequality and pregnancy. Environ Int 40:110–115. https://doi.org/10.1016/j.envint.2011.07.004
Baston D (2020) exactextractr: Fast Extraction from Raster Datasets using Polygons
Foster JE, Wolfson MC (2010) Polarization and the decline of the middle class: Canada and the U.S. J Econ Inequal 8:247–273. https://doi.org/10.1007/s10888-009-9122-7
Gascon M, Cirach M, Martínez D, Dadvand P, Valentín A, Plasència A, Nieuwenhuijsen MJ (2016) Normalized difference vegetation index (NDVI) as a marker of surrounding greenness in epidemiological studies: The case of Barcelona city. Urban For Urban Green, Special Section: Power in urban social-ecological systems: Processes and practices of governance and marginalization 19:88–94. https://doi.org/10.1016/j.ufug.2016.07.001
Gerrish E, Watkins SL (2018) The relationship between urban forests and income: A meta-analysis. Landsc Urban Plan 170:293–308. https://doi.org/10.1016/j.landurbplan.2017.09.005
Hoffimann E, Barros H, Ribeiro AI (2017) Socioeconomic Inequalities in Green Space Quality and Accessibility—Evidence from a Southern European City. Int J Environ Res Public Health 14(8):916. https://doi.org/10.3390/ijerph14080916
Hoffman JS, Shandas V, Pendleton N (2020) The effects of historical housing policies on resident exposure to intra-urban heat: A study of 108 US Urban Areas. Climate 8:12. https://doi.org/10.3390/cli8010012
Hwang H-J, Yook S-J, Ahn K-H (2011) Experimental investigation of submicron and ultrafine soot particle removal by tree leaves. Atmos Environ, Modeling of Air Quality Impacts, Forecasting and Interactions with Climate 45:6987–6994. https://doi.org/10.1016/j.atmosenv.2011.09.019
Hwang YH, Nasution IK, Amonkar D, Hahs A (2020) Urban green space distribution related to land values in fast-growing megacities, Mumbai and Jakarta–unexploited opportunities to increase access to greenery for the poor. Sustainability 12:4982. https://doi.org/10.3390/su12124982
James P, Banay RF, Hart JE, Laden F (2015) A review of the health benefits of greenness. Curr Epidemiol Rep 2:131–142. https://doi.org/10.1007/s40471-015-0043-7
Jelks NO, Jennings V, Rigolon A (2021) Green gentrification and health: A scoping review. Int J Environ Res Public Health 18:907. https://doi.org/10.3390/ijerph18030907
Jennings V, Bamkole O (2019) The Relationship between Social Cohesion and Urban Green Space: An Avenue for Health Promotion. Int J Environ Res Public Health 16:452. https://doi.org/10.3390/ijerph16030452
Kramer H, Oldengarm J (2010) UrbTree: a tree growth model for the urban environment. Carbon Balance Manage
Krieger N, Singh N, Waterman P. D (2016) Metrics for monitoring cancer inequities: residential segregation, the Index of Concentration at the Extremes (ICE), and breast cancer estrogen receptor status (USA, 1992-2012). Cancer Causes & Control 27:1139–1151
Kuo FE, Sullivan WC (2001) Environment and Crime in the Inner City: Does Vegetation Reduce Crime? Environ Behav 33:343–367. https://doi.org/10.1177/0013916501333002
Kye SH (2018) The persistence of white flight in middle-class suburbia. Soc Sci Res 72:38–52. https://doi.org/10.1016/j.ssresearch.2018.02.005
Link BG, Phelan J (1995). Social Conditions As Fundamental Causes of Disease. J Health Soc Behav 80–94. https://doi.org/10.2307/2626958
Massey D, Booth A, Crouter AC (2001) The Prodigal Paradigm Returns: Ecology Comes Back to Sociology. Does It Take Village Community Eff Child Adolesc Fam 41–48
McEachan RRC, Prady SL, Smith G, Fairley L, Cabieses B, Gidlow C, Wright J, Dadvand P, van Gent D, Nieuwenhuijsen MJ (2016) The association between green space and depressive symptoms in pregnant women: moderating roles of socioeconomic status and physical activity. J Epidemiol Community Health 70:253–259. https://doi.org/10.1136/jech-2015-205954
McPherson EG (1994) Benefits and costs of tree planting and care in Chicago. Gen Tech Rep, NE USA
Mitchell B, Franco J, Richardson J (2018) HOLC “redlining” maps: The persistent structure of segregation and economic inequality. https://doi.org/10.13140/RG.2.2.21841.48486
Nardone A, Casey JA, Morello-Frosch R, Mujahid M, Balmes JR, Thakur N (2020) Associations between historical residential redlining and current age-adjusted rates of emergency department visits due to asthma across eight cities in California: an ecological study. Lancet Planet Health 4:e24–e31. https://doi.org/10.1016/S2542-5196(19)30241-4
Nardone A, Rudolph KE, Morello-Frosch R, Casey JA (2021) Redlines and Greenspace: The Relationship between Historical Redlining and 2010 Greenspace across the United States. Environ Health Perspect 129:17006. https://doi.org/10.1289/EHP7495
Noszczyk T, Gorzelany J, Kukulska-Kozieł A, Hernik J (2022) The impact of the COVID-19 pandemic on the importance of urban green spaces to the public. Land Use Policy 113:105925. https://doi.org/10.1016/j.landusepol.2021.105925
Nowak DJ, Crane DE, Stevens JC (2006) Air pollution removal by urban trees and shrubs in the United States. Urban For Urban Green 4:115–123. https://doi.org/10.1016/j.ufug.2006.01.007
Oliveira EA, Andrade JS, Makse HA (2014) Large cities are less green. Sci Rep 4:4235. https://doi.org/10.1038/srep04235
Omernik JM (1987) Ecoregions of the Conterminous United States. Ann Assoc Am Geogr 77:118–125. https://doi.org/10.1111/j.1467-8306.1987.tb00149.x
Osborn CY, Groot M, de, Wagner JA (2013) Racial and ethnic disparities in diabetes complications in the Northeastern United States: The role of socioeconomic status. J Natl Med Assoc 105:51–58. https://doi.org/10.1016/S0027-9684(15)30085-7
Park J, Kim J-H, Lee DK, Park CY, Jeong SG (2017) The influence of small green space type and structure at the street level on urban heat island mitigation. Urban For Urban Green 21:203–212. https://doi.org/10.1016/j.ufug.2016.12.005
Pouso S, Borja Á, Fleming LE, Gómez-Baggethun E, White MP, Uyarra MC (2021) Contact with blue-green spaces during the COVID-19 pandemic lockdown beneficial for mental health. Sci Total Environ 756:143984. https://doi.org/10.1016/j.scitotenv.2020.143984
R Core Team (2022) R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing, Vienna, Austria
Ribeiro HV, Rybski D, Kropp JP (2019) Effects of changing population or density on urban carbon dioxide emissions. Nat Commun 10:3204. https://doi.org/10.1038/s41467-019-11184-y
Rowntree SB (1902) Poverty: A study of town life. Econ J 12:56–62. https://doi.org/10.2307/2957025
Schüle SA, Hilz LK, Dreger S, Bolte G (2019) Social inequalities in environmental resources of green and blue spaces: A review of evidence in the WHO European Region. Int J Environ Res Public Health 16:1216. https://doi.org/10.3390/ijerph16071216
Schwarz K, Fragkias M, Boone CG, Zhou W, McHale M, Grove JM, O’Neil-Dunne J, McFadden JP, Buckley GL, Childers D, Ogden L, Pincetl S, Pataki D, Whitmer A, Cadenasso ML (2015) Trees grow on money: Urban tree canopy cover and environmental justice. PLOS ONE 10:e0122051. https://doi.org/10.1371/journal.pone.0122051
Shepley M, Sachs N, Sadatsafavi H, Fournier C, Peditto K (2019) The impact of green space on violent crime in urban environments: An evidence synthesis. Int J Environ Res Public Health 16:5119. https://doi.org/10.3390/ijerph16245119
Twohig-Bennett C, Jones A (2018) The health benefits of the great outdoors: A systematic review and meta-analysis of greenspace exposure and health outcomes. Environ Res 166:628–637. https://doi.org/10.1016/j.envres.2018.06.030
US EPA (2015) Learn About Environmental Justice [WWW Document]. URL https://www.epa.gov/environmentaljustice/learn-about-environmental-justice (accessed 4.28.22)
US Geographical Survey (2022) AρρEEARS [WWW Document]. URL https://appeears.earthdatacloud.nasa.gov/ (accessed 9.29.22)
Venter ZS, Barton DN, Gundersen V, Figari H, Nowell M (2020) Urban nature in a time of crisis: recreational use of green space increases during the COVID-19 outbreak in Oslo, Norway. Environ Res Lett 15:104075. https://doi.org/10.1088/1748-9326/abb396
Venter ZS, Shackleton C, Faull A, Lancaster L, Breetzke G, Edelstein I (2022) Is green space associated with reduced crime? A national-scale study from the Global South. Sci Total Environ 825:154005. https://doi.org/10.1016/j.scitotenv.2022.154005
Wendelboe-Nelson C, Kelly S, Kennedy M, Cherrie JW (2019) A scoping review mapping research on green space and associated mental health benefits. Int J Environ Res Public Health 16:2081. https://doi.org/10.3390/ijerph16122081
Winkleby MA, Jatulis DE, Frank E, Fortmann SP (1992) Socioeconomic status and health: how education, income, and occupation contribute to risk factors for cardiovascular disease. Am J Public Health 82:816–820. https://doi.org/10.2105/AJPH.82.6.816
Wood E, Harsant A, Dallimer M, Cronin de Chavez A, McEachan RRC, Hassall C (2018) Not all green space is created equal: Biodiversity predicts psychological restorative benefits from urban Green Space. Front Psychol 9:2320
Zhang Y, Van den Berg AE, Van Dijk T, Weitkamp G (2017) Quality over quantity: Contribution of urban green space to neighborhood satisfaction. Int J Environ Res Public Health 14:535. https://doi.org/10.3390/ijerph14050535
Funding
Research reported in this publication was supported by the National Institute on Minority Health and Health Disparities [R00MD012615 (TTN), R01MD015716 (TTN), R01MD016037 (QCN], the National Library of Medicine [R01LM012849 (QCN)]. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.
Author information
Authors and Affiliations
Contributions
JCK: Formulated study design, executed all data analysis, drafted manuscript. TTN: Provided essential interpretation and feedback on study design, contributed to writing and drafting of the manuscript. Acquired Funding. QCN: Provided essential interpretation and feedback on study design, contributed to writing and drafting of the manuscript. YH: Senior author. Conceived of the presented idea. Supervised and advised study design, data analysis, manuscript drafting. Acquired funding. All authors discussed the results and contributed to the final manuscript.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Ethical approval
Statements of ethical approval and informed consent are not relevant for this research as it does not involve human subjects or non-publicly available data.
Informed consent
No human participants or sensitive data was involved in this study.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as 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 images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Kitch, J.C., Nguyen, T.T., Nguyen, Q.C. et al. Changes in the relationship between Index of Concentration at the Extremes and U.S. urban greenspace: a longitudinal analysis from 2001–2019. Humanit Soc Sci Commun 10, 862 (2023). https://doi.org/10.1057/s41599-023-02115-w
Received:
Accepted:
Published:
DOI: https://doi.org/10.1057/s41599-023-02115-w