Next Article in Journal
Validity of the International Physical Activity Questionnaire Long Form for Assessing Physical Activity and Sedentary Behavior in Subjects with Chronic Stroke
Next Article in Special Issue
Neighborhood Characteristics and Cancer Survivorship: An Overview of the Current Literature on Neighborhood Landscapes and Cancer Care
Previous Article in Journal
Associations between Prenatal Exposure to Phthalates and Timing of Menarche and Growth and Adiposity into Adulthood: A Twenty-Years Birth Cohort Study
Previous Article in Special Issue
Neighborhood Disadvantage and Tobacco Retail Outlet and Vape Shop Outlet Rates
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJERPH
Volume 18
Issue 9
10.3390/ijerph18094728
ijerph-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Measuring Neighborhood Landscapes: Associations between a Neighborhood’s Landscape Characteristics and Colon Cancer Survival

by
Daniel Wiese
1,*,
Antoinette M. Stroup
2,3,
Aniruddha Maiti
4,
Gerald Harris
2,
Shannon M. Lynch
5,
Slobodan Vucetic
4,
Victor H. Gutierrez-Velez
1 and
Kevin A. Henry
1,5
1
Department of Geography and Urban Studies, Temple University, Philadelphia, PA 19122, USA
2
New Jersey Department of Health, New Jersey State Cancer Registry, Trenton, NJ 08625, USA
3
Rutgers School of Public Health, Rutgers Cancer Institute of New Jersey, New Brunswick, NJ 08901, USA
4
Department of Computer and Information Sciences, Temple University, Philadelphia, PA 19122, USA
5
Fox Chase Cancer Center, Division of Cancer Prevention and Control, Philadelphia, PA 19111, USA
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2021, 18(9), 4728; https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph18094728
Submission received: 25 March 2021 / Revised: 21 April 2021 / Accepted: 27 April 2021 / Published: 29 April 2021
(This article belongs to the Special Issue Impacts of the Neighborhood Environment on Cancer Prevention, Outcomes and Disparities)
Browse Figures
Versions Notes

Abstract

:
Landscape characteristics have been shown to influence health outcomes, but few studies have examined their relationship with cancer survival. We used data from the National Land Cover Database to examine associations between regional-stage colon cancer survival and 27 different landscape metrics. The study population included all adult New Jersey residents diagnosed between 2006 and 2011. Cases were followed until 31 December 2016 (N = 3949). Patient data were derived from the New Jersey State Cancer Registry and were linked to LexisNexis to obtain residential histories. Cox proportional hazard regression was used to estimate hazard ratios (HR) and 95% confidence intervals (CI95) for the different landscape metrics. An increasing proportion of high-intensity developed lands with 80–100% impervious surfaces per cell/pixel was significantly associated with the risk of colon cancer death (HR = 1.006; CI95 = 1.002–1.01) after controlling for neighborhood poverty and other individual-level factors. In contrast, an increase in the aggregation and connectivity of vegetation-dominated low-intensity developed lands with 20–<40% impervious surfaces per cell/pixel was significantly associated with the decrease in risk of death from colon cancer (HR = 0.996; CI95 = 0.992–0.999). Reducing impervious surfaces in residential areas may increase the aesthetic value and provide conditions more advantageous to a healthy lifestyle, such as walking. Further research is needed to understand how these landscape characteristics impact survival.

1. Introduction

Neighborhood characteristics can capture social, physical, and economic conditions of the environment in which a person lives [1]. According to Northridge et al. [2], social context and the built environment are two intermediate factors that influence health and well-being at the individual and population level. The social context of a neighborhood is most often defined through its socio-economic and demographic composition using census data (e.g., median income, % living below poverty) [3]. The built environment includes physical characteristics of neighborhoods (i.e., land covers) that are human-made or modified [4], such as buildings, roads, housing conditions, parks sidewalks or greenspace, that can potentially provide the setting for human activity [5].
Therapeutic landscape theory [6] can add depth to our understanding of how landscapes impact health and cancer outcomes [7]. Gesler’s therapeutic landscape theory posits that social and spatial factors are interconnected, and modification of place can have various effects on human health. Moreover, people may have experienced different feelings and healing effects in various places caused by perceptions of landscape aesthetics [6]. Some studies linked the aesthetical value of the landscapes, analyzed through photography of micro-landscapes, to emotional load and restoration [8,9,10,11]. Lee et al. [12] found that neighborhood satisfaction is related to landscape structure. Others suggest that configuration and composition of different types of land cover and land use classes could be influential on place perception and emotions because “perception, cognition, and evaluation are highly interrelated processes” [13], and perception of environmental landscape is important (among others) for movement and social interactions [14,15,16]. In addition, factors such as natural habitat fragmentation and impervious surfaces may also influence health outcomes considering the relation between urban design and mental health [17].
For numerous cancer types, studies report significant associations between neighborhood socio-economic conditions and cancer mortality and survival [18,19,20,21,22,23,24,25,26,27,28,29]. Generally, these studies have found that people living in low socio-economic-status neighborhoods have significantly higher mortality or shorter survival after diagnosis. The relationship is less consistent between the built environment and cancer mortality and survival. Several studies have examined these associations. For example, a study focused on greenspace availability and accessibility by James et al. [30] found that higher levels of green vegetation were associated with decreased cancer mortality. Several other studies on lung cancer-specific mortality did not find any significant associations [31,32,33,34]. For breast cancer, Keegan et al. [35] did not find any positive influence of park availability and survival, but did report an increased risk of breast cancer death in areas with higher traffic/road density, possibly through discouraging recreational-based physical activity.
Among the studies that have examined the relationship between the built environment and non-cancer health outcomes (e.g., cardiovascular and/or mental health illnesses), most have focused on greenspace. Mears et al. [36] summarized positive effects of the urban greenspace on health outcomes, including the reduction of traffic pollution and a reduction of heat-island effects, as well as an increase of emotional recreation effect and physical or social activities. In one study, Bratman et al. [37] examined evidence across the natural, social, and health sciences on the impacts of nature experience on mental health. They argue that the configuration and conception of greenspace is essential for human well-being. Kondo et al. [38] also concluded that navigating through urban greenspaces causes more positive emotions compared to the built urban environment. Additionally, the quality of the greenspace (e.g., well-maintained parks) has a positive influence for psychological well-being [39,40,41,42]. Moreover, urban gardens and well-maintained front yards may influence population health through the pathways of aesthetics and promote walking and outdoor physical activity [43]. Additionally, walkability, mix and type of businesses, and land use composition of a neighborhood may encourage healthy behavior and ‘facilitate integration of habits into a daily lifestyle” (p. 76) [5]. Given the importance of greenspace on health and a lack of studies focused specifically on cancer, more research is needed examining the relationship between the built environment and cancer outcomes in order to better understand these relationships and underlying pathways [44].
With the exception of a few studies [30,35,45], insufficient attention has been given to integration of both social and environmental landscape characteristics on cancer outcomes. Considering that cancer survival has shown to be associated with mental well-being [46,47], integrating landscape and built-environmental characteristics may be useful for quantifying neighborhood quality and help measure the effect of neighborhood land cover configuration and composition on cancer outcomes.
Estimating the aesthetical value of a landscape and finding the best measure for defining the built environment can be challenging because there are infinite ways to measure and operationalize the data in research [48] and further distil the essential components that impact health and health behaviors. Frequently, the built environment is integrated into modeling using census data such as housing density and quality [44]. However, recent technologies in image classification allow an estimation of neighborhood physical disorder (e.g., presence of abounded buildings, non-maintained roads) or greenspace availability and quality using remote sensing products and street view in health research [49,50,51]. Despite the increasing literature, however, when examining the role of greenspace on health outcomes, several gaps remain, including “standardization” of appropriate measures [38,52]. Therefore, it is important to develop clear concepts and metrics for quantifying and measuring the effects of a neighborhood’s landscape and the built environment on health outcomes to understand the underlying pathways.
Landscape metrics are a commonly used technique in landscape ecology—an environmental science approach of landscape characterization, evaluation, and design [53] to quantify landscape characteristics and features. This concept is frequently applied in urban planning, biodiversity, species richness, and conservation and infectious disease epidemiology studies [54,55]. Landscape metrics are also common in urban aesthetics evaluation [8,9]. To date, there have been no population-based cancer studies that integrate landscape metrics.
In this study, we examine associations between regional-stage colon cancer survival and several landscape metrics that quantify neighborhood’s built-environment using cancer surveillance data from the New Jersey State Cancer Registry, residential histories, and land cover and land use data from the National Land Cover Database.

2. Materials and Methods

2.1. Study Population

The NJSCR provided all colon cancer cases. The NJSCR is a population-based cancer registry that collects and maintain incidence data for the State of New Jersey. It was established in October 1978 and monitors cancer among the more than 8.9 million residents of New Jersey [56]. The study population includes all New Jersey residents, 18 years and older, with histologically confirmed, first primary regional stage colon cancer, as defined according to the International Classification of Diseases for Oncology, 3rd Edition (ICD-O3 C180–C189, C260; excluding histology codes 9050–9055, 9140, 9590–9992) [57] diagnosed between 1 January 2006 and 31 December 2011 (N = 4041). We restricted our analysis to regional stage colon cancers to simplify the interpretation of the results and reduce systematic variation in survival that could be attributable to stages at diagnosis. The study was reviewed and/or approved by Temple and Rutgers University Institutional Review Boards.
Individual-level factors included age at diagnosis, gender (male, female), race/ethnicity (Non-Hispanic (NH) White, NH Black, NH Asian/Pacific Islander (API), NH Other, and Hispanic (any race)), and vital status, including date of death and cause of death (if deceased) or date of last contact (if alive). Cases were followed until their deaths, relocation from the state of New Jersey, or until 31 December 2016. Deaths attributed to colon cancer were coded based on ICD-10 code C18 [57]. NJSCR regularly updates the database through linkages with state and national sources, including death data from the New Jersey Department of Health Office of Vital Statistics and Registry and the National Death Index, hospital discharge files, Centers for Medicare and Medicaid Services, Social Security Administration Services for Epidemiologic Researchers, and motor vehicle registration files [56].

2.2. Residential Histories

The NJSCR linked the study population to the commercial residential history database developed by LexisNexis, Inc. (Miamisburg, OH, USA; https://0-www-lexisnexis-com.brum.beds.ac.uk/en-us/products/public-records.page, accessed on 28 April 2021), similar to earlier studies [58,59,60,61]. The majority of regional stage colon cancer cases [n = 3949 (97.6%)] had residential information available for up to 20 of the most recent addresses between 1946 and 2018, with documented start and stop dates. A preprocessing technique described by Wiese et al. [24] was applied to establish a complete residential timeline for the time after diagnosis. All residential addresses were geocoded to the 2010 census tract (CT) boundaries using the North American Association of Central Cancer Registries (NAACCR) AGGIE Geocoder [62].

2.3. Socio-Economic Variables

The socio-economic variables included a widely used census tract poverty (CT-poverty) variable, defined as the proportion of population 18 years and older living below the Federal poverty level. Additionally, we included information on census tract measures of housing density and median year housing built. The required variables were obtained from publicly available U.S. Census and American Community Survey (ACS) data. U.S. Census 2010 and the ACS 5-year average data 2006–2010, 2007–2011, 2008–2012, 2009–2013, and 2010–2014 were used for residencies between 2006 and 2010. ACS 2011–2015, 2012–2016, and 2013–2017 were used for residencies between 2011 and 2016.
For each case, all corresponding CTs during the follow-up period were included, and every residential record received a corresponding value based on the earliest date of the residential appearance in the data set. If a patient remained at the same diagnosis CT over multiple years, we assigned annual values to capture changes within the neighborhood that could be caused by the gentrification.

2.4. Environmental Variables

The neighborhood built-environment factors were measured by using landscape metrics that captured characteristics such as land cover and land use composition and configuration of the neighborhood considering local spatial patterns.
The National Land Cover Database (NLCD) was used to extract land cover classes and proportion of imperviousness for the years 2006 (used for residential records 2006–2009), 2011 (used for residential records 2010–2014), and 2016 (used for residential records 2015–2016). The NLCD products are freely available raster files of 30 m spatial resolution that were classified using Landsat-based satellite imagery by the U.S. geological survey [63]. The land cover raster includes 16 categorical classes and the imperviousness raster, which has a continuous scale of proportion of impervious surfaces per pixel. To reduce the number of categories for further analysis, we reclassified the original NLCD raster using R package raster [64] into 7 classes (Forest, Grass, Shrubs, Developed Lands (Open, Low-, Medium-, and High-Intensity) (Table 1). Land use information was obtained from the New Jersey Geographic Information Network (njgin.nj.gov/njgin/edata/parcels/#!/, accessed on 28 April 2021), which includes publicly available information on parcel use for the state of New Jersey as of the year 2019. The original shapefile (spatial polygon) was rasterized (i.e., converted) using R package raster [64] based on the parcel use category and reclassified by keeping only industrial and commercial lands.
Previous studies have already examined the NLCD’s classification of developed lands to estimate greenspace availability in urban areas and have noted that open and low-intensity developed lands are suitable for identification of greenspaces and trees [34,65,66,67]. Typically, areas with predominantly large housing parcels, roads surrounded by greenspaces, urban housing with larger backyards, or isolated large roads would be classified as developed open land (NLCD class 21), assuming a larger proportion of greenspaces than of impervious/built surfaces (max. 20%) in the area. In contrast, a central business district, densely built inner-city housing, or a largely expanded shopping mall would be classified as high-intensity developed lands (NLCD class 24) with 80–100% imperviousness.
Landscape metrics on patch (square) and class levels (greenspaces) within a landscape (i.e., census tract) were calculated using R statistical software [68] (version 4.0.1), implementing packages landscapemetrics [69] and SDMTools [70]. Landscape metrics include more than 50 measures. Because most variables are highly correlated or difficult to interpret and compare between landscapes because of the open, non-fixed range scale [53,71], we developed a list of the influential landscape metrics, as defined by previous studies [36,39,42]. Table 2 summarizes all landscape metrics that were considered in the present study. Selected landscape metrics and the land cover classification were also mapped using QGIS v.3.10 (Figure 1).

2.5. Statistical Analysis

The survival time for every patient was calculated in months as the difference between the date of diagnosis and the date of death or date of last contact. Cases missing survival time (i.e., only ascertained through death certificates or autopsy) were excluded from this analysis. Patients were censored at the date of death if they died from causes other than colon cancer, the date the patient was lost to follow-up, at the end of the follow-up period (31 December 2016), or at the time of the relocation from the State of New Jersey, whichever occurred first. Additionally, we calculated time intervals to every CT location and assigned start and end dates based on residential histories. Every time interval received a corresponding socio-economic and environmental neighborhood value, as well as the vital status of the patient (1-dead, 0-alive).
We designed a process for variable selection and evaluation with minor modifications based on methodology from a previous neighborhood wide association study or NWAS [80]. All methodological steps are summarized in Figure 2. First, we developed a series of univariate, crude models to estimate the effect of each selected landscape metric and other neighborhood variables on the duration of survival time using Cox proportional hazard regression for time-varying covariates [81]. Cox proportional hazard regression is a widely used semi-parametric time-to-event modeling technique, where death is considered being the event. Cox proportional hazard regression allows incorporation of individual and area-level covariates. Additionally, it does not require the definition of a probability distribution in advance, and is suitable for time-varying covariates [82]. Variables were selected if they reached significance at p < 0.05.
Then we applied Spearman correlation analysis of the selected variables because many landscape metrics are highly correlated and redundant. After excluding all highly correlated variables (r2 > 0.7), the number of variables was reduced from 14 to 8 (Figure 3). The remaining variables included aggregation indices of all four intensity levels of developed lands, proportion of high intensity developed lands, forest contiguity index, Shannon diversity index, the patch richness density, and CT-poverty (Figure 1).
We then developed a set of models that included all individual-level variables (sex, age and race/ethnicity, sub-stage, and mover status) and each neighborhood variable or landscape metric. To estimate the risk of death from colon cancer by each individual and area-based variable, all coefficients were exponentiated and expressed as Hazard Ratios (HRs). For continuous variables, positive HRs indicate a positive association with the increase in risk of death. For categorical variables, the HRs are compared to the reference group. An HR = 1 indicates that the risk is similar across groups [83].
All models were run using R package survival [84] and survsim [81], and met the proportional hazard assumption based on the examination of Schoenfeld residuals using the cox.zph() function in the R survsim [81].

3. Results

3.1. Study Population

The study population included 3949 regional stage colon cancer cases. The spatially interpolated distribution of cases and total population density are presented in Figure S1. There were slightly fewer males (47.6%) than females (52.4%), and about three quarters (73.6%) were NH-White, 12.4% NH-Black, 8.1% Hispanic origin (any race), 3.6% NH-API, and 2.4% Other race. Approximately a third (27.5%) of all the patients died from the colon cancer by the end of follow-up, with a median survival of 66 months (range 1–139). During the follow-up period, 65.5% remained at their diagnosis CT, and 12.1% left New Jersey during the study period. Among those who moved, 18.5% only moved once, 12.3% moved twice, and 3.6% moved three or more times after cancer diagnosis. The average time spent at the CT at diagnosis was 7.5 years (Table 3).
The distribution of colon cancer patients by neighborhood/landscape characteristics at the time of diagnosis showed that the majority of patients (73.9%, n = 2919) lived in areas with a poverty level of less than 10%. Additionally, most patients (82.5%, n = 3259) were living in neighborhoods with a relatively low (≤30%) proportion of open developed lands (areas with large greenspace cover), while 89.8% (n = 3547) of patients living in neighborhoods where more than 30% of the total landscape was dominated by high-intensity developed lands (less than 20% greenspaces). Twenty-four percent (n = 939) of all colon cancer patients were residents in a neighborhood with no tree cover. The distribution of the study population for these and other landscape metrics are summarized in Figure 4.

3.2. Univariate Models

Figure 5 summarizes the model results for each variable. For every 10% increase in housing density, the risk of death increased by 5% (HR = 1.005; 95% CI = 1.002–1.009). For CT-poverty, cases living in CT with poverty levels in the range of 10–20% had a higher risk of death than those living in CTs with poverty levels in the range of 0–5% (HR = 1.23; 95% CI = 1.04–1.44).
Of the 27 landscape metrics, six were statistically significant. An increasing proportion of the high- and medium-intensity developed lands were positively associated with the risk of death (HR = 1.007; 95% CI = 1.003–1.01 and HR = 1.008; 95% CI = 1.005–1.011, respectively). Similarly, a positive association was found between an increasing aggregation index of high- and medium-intensity developed lands and risk of death (HR = 1.005; 95% CI = 1.001–1.009; HR = 1.01; 95% CI = 1.007–1.02, respectively). There was also a positive association between increasing patch richness density (HR = 1.02; 95% CI = 1.01–1.03) and elevated risk of death. Average imperviousness was also significant, indicating an approximate 7% risk increase for every 10% increase of the neighborhood’s imperviousness (HR = 1.007; 95% CI = 1.004–1.011).
Seven landscape metrics were significantly associated with a decrease in risk of colon cancer death, including aggregation of low-intensity developed lands (HR = 0.996; 95% CI = 0.991–0.999), proportion and aggregation of open developed lands (HR = 0.992; 95% CI = 0.988–0.997, HR = 0.996; 95% CI = 0.991–0.999, respectively), Shannon diversity index (HR = 0.79; 95% CI = 0.64–0.98), and proportion and contiguity (i.e., connectivity) of forest/trees (HR = 0.995; 95% CI = 0.992–0.998, HR = 0.997; 95% CI = 0.996–0.998, respectively).

3.3. Multivariate Models

In the last step, we developed ten multivariate models, adjusting for individual level variables (age and sub-stage, sex/gender, race/ethnicity, and mover status), CT-poverty, and each landscape metric. The individual level HRs are provided in Supplementary Table S1.
In multivariate models, the proportion of high-intensity developed lands had the strongest association with the colon cancer risk of death. For every 10% increase in the proportion, the risk of death increased by 6% (HR = 1.006; 95% CI = 1.002–1.01). Additionally, increasing aggregation index (i.e., more compact areas) of high-intensity developed lands was positively associated with the risk of colon cancer death (HR = 1.005; 95% CI = 1.001–1.009). Similar associations were found for the increasing aggregation index of medium intensity developed lands (HR = 1.009; 95% CI = 1.003–1.015). In contrast, a 6% decrease in risk of colon cancer death was estimated for every 10% increase in the proportion of open developed lands (HR = 0.994; 95% CI = 0.988–0.999). Increasing aggregation index of low-intensity developed lands (HR = 0.995; 95% CI = 0.99–0.999) and increasing forest contiguity index (HR = 0.998; 95% CI = 0.996–0.999) were all significantly negatively associated with the risk of colon cancer death. CT-poverty was no longer significant in multivariable models (Figure 6).

4. Discussion

To date, only a few studies have integrated landscape characteristics into cancer disparities research. Traditionally, socio-epidemiological neighborhood studies focus on examining associations between neighborhood factors and cancer survival using Census-based socio-economic status data. However, the physical environment [85], recreational activities [86], and urban design are important components in the selection of residence, especially among older adults [87]. While experiences of place are recognized and there have been several attempts to incorporate the role of the landscape characteristics into cancer research [35,44], often it is done through the utilization of road networks or based on parcel use, greenspace availability, and accessibility. Only few attempts were made using satellite imagery-based products [45,88,89].
We evaluated the relationship between the risk of death from colon cancer and several area-based landscape characteristics that attempt to describe the configuration and composition of the built environment after adjusting for individual-level factors and CT-poverty. This study adds to a body of literature on the effects of neighborhood landscape characteristics on cancer survival. Independent of CT-poverty and several individual demographic and prognostic factors, we found a significant relationship between the risk of death from colon cancer and the proportion and aggregation of high-intensity developed lands (i.e., areas dominated by buildings and roads). The risk of death increased as the proportion and aggregation of high-intensity developed lands increased. The reasons for this relationship may be attributed to a low prevalence of greenspaces in areas with more high-intensity developed lands. Less greenspace in these places could reduce access to recreational sites that promote physical activity. Furthermore, such places may evoke negative emotional feelings and psychological well-being that have been found to be associated with a lack of greenspace [37]. Creating landscapes that promote exercise and active transport such as biking and walking is particularly important to cancer patients. Many clinical trials have shown the positive effects that walking and physical activity can have on quality of life and survival time in cancer patients [90], and are known to reduce risk of colorectal cancer development [91,92].
The relationship of worse survival with increasing high-intensity developed lands independent of CT-poverty is likely related to what this metric is capturing. High-intensity developed lands are characterized by places that include less than 20% greenspace of total land area [63]. An increasing proportion and aggregation of high-intensity developed lands is generally characterized by very compact and large areas within neighborhoods, with few or no parks and where the majority of greenspace is a result of overgrown vegetation from abandoned lots. A study from Philadelphia (PA, USA) reported that abandoned buildings and lots were associated with negative health outcomes among residents because of dangerous physical and social environments and sanitation/garbage issues [93]. High-intensity developed lands are also characterized by large parking lots and multiple lane roads that can reduce the aesthetic value of a neighborhood. This can lead to a reduction in walkability. Several previous studies reported the importance of aesthetic value and quality of the neighborhood environment on mental health and physical activity [36,37,42,43,52,65,94]. Additionally, we cannot exclude that this finding could also result from unmeasured confounding factors such as environmental pollution in highly urbanized areas or food deserts.
The better survival in neighborhoods with a higher proportion of open developed lands is likely a result of the large and compact greenspaces found in these areas, which provide potential for recreational activities. Larger open developed lands also assume more contiguous areas with a large proportion of street tree canopy, urban gardens, or large backyards across the landscape (neighborhood). Having vegetation along the roads, green front and backyards could have positive emotional effects [36,43] essential for many cancer survivors. Kondo et al. [38] also conclude that navigating through urban greenspaces compared to built urban spaces leads to more positive emotions. Additionally, this could suggest that the availability of green infrastructure and neighborhood parks would increase walkability and exercising [43,95,96].
We also found that increasing the forest contiguity index was positively associated with a decrease in the risk of death. This suggests that lower shape complexity and lower interspersion of greenspaces (i.e., larger and better-connected greenspaces) may further decrease the risk of colon cancer death because of the availability of green corridors within cities. Large contiguous green space leads to more varied use and extended use for physical activity and provides an opportunity to strengthen social capital that improves survival. This aligns with an earlier study reporting a positive effect on general and mental health from having fewer, but larger patches/areas of greenspaces, rather than a high density of small patches [39]. This finding confirms our earlier hypothesis that lack of greenspaces in neighborhoods negatively influences colon cancer patients. As Mary Soderstrom [97] argues, increasing street greenery and number of greenspaces can make a difference and create a dense and pleasurable city at the same time.
While access to and availability of greenspace is recognized as an (predominantly positive) influential factor on health conditions, the reduction of impervious surfaces and spaces occupied by oversized residential-area roads and gigantic parking lots might be as important as the creation of green infrastructure through zoning or the re-use of abandoned lots. Unfortunately, this practice is rare in the U.S., but more attention must be given to the aesthetics of the neighborhoods. Therefore, we agree with Kondo et al. [38] that “urban planners and public health professionals need evidence of the impacts of specific therapeutic or place-based interventions to help address public health issues facing their constituents” (p. 22) and argue that integration of landscape metrics into health disparities research may provide the required evidence.
Remote sensing-based classification is a valuable tool for land cover analysis and can be customized depending on the research question and study area. Additionally, there is a growing amount of data on a high spatial resolution allowing fine classifications. However, the land cover classification schema developed by the NLCD is informative for analyzing urban and urbanized areas, and offers enormous opportunities for integration of land cover data into cancer research also on nationwide scale, including non-contiguous states and Puerto Rico [63]. Additionally, it is available for several time points, which allows temporal analysis and application of longitudinal study design like in our example. On the other hand, the calculation of various landscape metrics allows a straightforward integration of several measures of landscape characteristics. These could then be helpful for city planning and the establishment of specific cancer prevention and control strategies.
Moreover, while aesthetical value of a neighborhood is a very subjective measure and typically requires qualitative interviews or surveys, some landscape metrics offer an opportunity for quantification. The contiguity or aggregation indices used here may provide information about landscape configuration [53,71] and become an alternative to more complicated measures derived from Google street view, which are time-consuming and more expensive to process.
The present study has several limitations. While we found evidence for associations between the land cover configuration and risk of death from colon cancer, the results may not be generalizable to other states with different demographic, socio-economic, or landscape characteristics because the study population was limited to New Jersey only. Additionally, we did not have access to individual-level factors such as individual socio-economic status, general health conditions (e.g., obesity data) and behavior. Not only could these factors potentially confound the relationship between landscape characteristics and colon cancer survival, but accounting for individual-level factors may reduce the geographic variance explained by neighborhood and landscape characteristics [98].
Another limitation is that we restricted our study population to regional-stage colon cancer cases and followed patients for only 10 years after diagnosis. This was necessary to minimize extreme variations in survival and limit sources of variation in residential history measures. Therefore, analyzing colon cancer cases diagnosed in earlier (local) or later (metastatic) stages of the disease may result in a different conclusion.
Moreover, we utilized only residential histories collected from LexisNexis. We did not have access to self-reported information and could not validate and/or augment LexisNexis data. However, according to previous studies, the concordance between LexisNexis addresses and addresses collected from study participants (85–86%) is high [58,59]. In our study, we could only validate the residential location at the time of diagnosis between the LexisNexis and the NJSCR. The concordance rate was approximately 83%. Opening to a 6-month window before and after the diagnosis date substantially increased the concordance rate to 93%. Only 8% of the locations from LexisNexis cases did not match any locations reported by the NJSCR. Several factors such as incorrect links at LexisNexis, incorrect geocodes assigned by NJSCR for both registry and LexisNexis residential addresses, incorrect addresses reported to the registry by hospitals and other reporting facilities or geocodes assigned to addresses based on post office boxes could be the reason for address discordance. However, the extent of the bias in either direction would be minimal because of the low proportion of affected cases. In contrast to the residential histories from self-reported data, LexisNexis is an objective source of residential history data, and is not sensitive to potential recall bias.
Additionally, the application of landscape metrics in neighborhood research is not typical and is more common in natural landscapes for ecological analysis and modeling. In this study, we used census tract boundaries for defining neighborhoods or landscapes. The definition of a landscape as a census tract is a subject of modified area unit problem. The selection of the landscape boundaries is essential for calculating landscape metrics, and change in area size and shape will affect the values of multiple indices. Defining landscapes through a use of other administrative boundaries or grid system may result in different conclusions. However, census tract is a common unit of analysis in public health research and allows an uncomplicated merge of data from various sources.
Lastly, the spatial resolution of the land cover classification raster was 30 m and suggests that all features within an area of 900 m2 are generalized and defined as one class. The establishment of land cover classification with a fine spatial resolution would result in higher precision in the classification of the ground objects. However, 30 m spatial resolution is widely used in remote sensing discipline for land cover classification, and the utilized NLCD dataset is a well-known high-quality product.
The calculation of landscape metrics can be done on any spatial resolution, but values may vary with change in pixel size. More challenging is the selection of landscape metrics itself. There are many measures on various geographic levels. Thereby, most metrics are highly correlated and redundant. Previous research in landscape ecology suggests that a minimal number of metrics (e.g., number of patch types, mean edge/area ratio, contagion, average patch shape, fractal measurements) would be sufficient to quantify spatial heterogeneity [71,72,73]. However, in a public health context, the selection of landscape metrics should be done more carefully, selecting meaningful variables that can be easily translated to the policy makers and urban planners. Selection of other landscape metrics could result in different associations or cause complications in interpretation.

5. Conclusions

The associations between neighborhood socio-economic status and cancer survival are fairly established, where increased risk of death is often associated with high neighborhood poverty. However, the neighborhood environment is not limited to socio-demographic factors because buildings, roads, greenspaces, and other human-made objects dominate landscapes, especially in urban and urbanized areas. It is essential to understand the relationship between landscape characteristics and health outcomes to develop new policies in urban planning and design essential for population health, especially in urban and urbanized areas. Our results suggest that increasing proportions and connectivity of urban greenspaces may substantially decrease the risk of colon cancer death. This association did not change, even after adjusting for neighborhood poverty, which is reported to be associated with a lack of greenspaces in urban and suburban areas [99], and reflected in our correlation analysis. The integration of remote sensing-based products, the NLCD and the calculation of landscape metrics allow the exploration of undiscovered pathways between neighborhood characteristics and colon cancer survival and should be further evaluated in neighborhood studies with other cancer sites and outcomes such as stage at the diagnosis. Additionally, further research is needed to understand how these specific landscape characteristics impact survival, and evaluate opportunities for developing a socio-environmental deprivation index combing census-based variables and land cover metrics in order to identify neighborhoods in need of interventions. Moreover, future studies should include additional neighborhood variables, especially related to walkability, that could help to evaluate the association between neighborhood built environment and colon cancer survival.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/ijerph18094728/s1, Table S1: Hazard Rates of individual-level factors from the multivariate model. Figure S1: Distribution of population density by municipality and interpolated density of colon cancer patients’ residencies.

Author Contributions

Conceptualization, D.W.; methodology, D.W., V.H.G.-V., and K.A.H.; software, D.W. and A.M.; validation, K.A.H., V.H.G.-V., S.M.L., and A.M.S.; data curation, G.H.; writing—original draft preparation, D.W.; writing—review and editing, K.A.H., V.H.G.-V., A.M.S., S.M.L., and S.V.; visualization, D.W.; supervision, K.A.H., V.H.G.-V., and S.M.L.; project administration, K.A.H.; funding acquisition, K.A.H., S.V., and A.M.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Science Foundation (1560888). NJSCR data were collected using funding from NCI and the Surveillance, Epidemiology, and End Results (SEER) Program (HHSN261201300021I), and the CDC’s National Program of Cancer Registries (NPCR) (5U58DP003931-02/NU58DP006279-02-00), as well as the State of New Jersey and the Rutgers Cancer Institute of New Jersey.

Institutional Review Board Statement

The study was reviewed and/or approved by Temple and Rutgers University Institutional Review Boards.

Informed Consent Statement

Patient consent was waived per Institutional Review Board (IRB) policies and procedures governing minimal-risk, analysis-only studies using existing data sources. Authors received data through an honest-broker protocol and, therefore, did not have access to identifiable patient data.

Data Availability Statement

Due to the nature of this research, supporting patient data is not available.

Acknowledgments

We thank all members of the Remote Sensing and Sustainability Lab (RSENSUS, rsensus.org) at Temple University (Wilson Lara-Henao, Jeronimo Rodriguez, Matt Marcus, Victoria Sarmiento, Angela Mejia, Evelyn Gorey, and Alexis Saenz) for the technical support and recommendations.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Diez Roux, A.V. Investigating neighborhood and area effects on health. Am. J. Public Health 2001, 91, 1783–1789. [Google Scholar] [CrossRef]
  2. Northridge, M.E.; Sclar, E.D.; Biswas, P. Sorting out the connections between the built environment and health: A conceptual framework for navigating pathways and planning healthy cities. J. Urban Health 2003, 80, 556–568. [Google Scholar] [CrossRef] [Green Version]
  3. Link, B.G.; Phelan, J. Social Conditions As Fundamental Causes of Disease. J. Health Soc. Behav. 1995, 80–94. [Google Scholar] [CrossRef] [Green Version]
  4. Rapoport, A. The Mutual Interaction of People and Their Built Environment; De Gruyter Mouton: Berlin, Germany, 2011. [Google Scholar]
  5. Renalds, A.; Smith, T.H.; Hale, P.J. A systematic review of built environment and health. Fam. Community Health 2010, 33, 68–78. [Google Scholar] [CrossRef]
  6. Gesler, W.M. Therapeutic landscapes: Medical issues in light of the new cultural geography. Soc. Sci. Med. 1992, 34, 735–746. [Google Scholar] [CrossRef]
  7. Matsuoka, R.; Sullivan, W. Human Psychological Community Health. In The Routledge Handbook of Urban Ecology; Douglas, I., Goode, D., Houck, M.C., Maddox, D., Eds.; Taylor & Francis Group: London, UK, 2011. [Google Scholar]
  8. Dramstad, W.E.; Tveit, M.S.; Fjellstad, W.; Fry, G.L. Relationships between visual landscape preferences and map-based indicators of landscape structure. Landsc. Urban Plan. 2006, 78, 465–474. [Google Scholar] [CrossRef]
  9. Palmer, J.F. Using spatial metrics to predict scenic perception in a changing landscape: Dennis, Massachusetts. Landsc. Urban Plan. 2004, 69, 201–218. [Google Scholar] [CrossRef]
  10. Staats, H.; Kieviet, A.; Hartig, T. Where to recover from attentional fatigue: An expectancy-value analysis of environmental preference. J. Environ. Psychol. 2003, 23, 147–157. [Google Scholar] [CrossRef]
  11. van den Berg, A.E.; Koole, S.L.; van der Wulp, N.Y. Environmental preference and restoration: (How) are they related? J. Environ. Psychol. 2003, 23, 135–146. [Google Scholar] [CrossRef]
  12. Lee, S.-W.; Ellis, C.D.; Kweon, B.-S.; Hong, S.-K. Relationship between landscape structure and neighborhood satisfaction in urbanized areas. Landsc. Urban Plan. 2008, 85, 60–70. [Google Scholar] [CrossRef]
  13. Kaplan, S. Aesthetics, Affect, and Cognition:Environmental Preference from an Evolutionary Perspective. Environ. Behav. 1987, 19, 3–32. [Google Scholar] [CrossRef] [Green Version]
  14. Ittelson, W.H. Environmental Perception and Urban Experience. Environ. Behav. 1978, 10, 193–213. [Google Scholar] [CrossRef]
  15. Zube, E.H.; Sell, J.L.; Taylor, J.G. Landscape perception: Research, application and theory. Landsc. Plan. 1982, 9, 1–33. [Google Scholar] [CrossRef]
  16. Iverson Nassauer, J. Culture and changing landscape structure. Landsc. Ecol. 1995, 10, 229–237. [Google Scholar] [CrossRef]
  17. Jackson, L.E. The relationship of urban design to human health and condition. Landsc. Urban Plan. 2003, 64, 191–200. [Google Scholar] [CrossRef]
  18. Singh, G.K. Area Deprivation and Widening Inequalities in US Mortality, 1969–1998. Am. J. Public Health 2003, 93, 1137–1143. [Google Scholar] [CrossRef]
  19. Singh, G.K.; Jemal, A. Socioeconomic and Racial/Ethnic Disparities in Cancer Mortality, Incidence, and Survival in the United States, 1950–2014: Over Six Decades of Changing Patterns and Widening Inequalities. J. Environ. Public Health 2017, 2017, 19. [Google Scholar] [CrossRef]
  20. Singh, G.K.; Williams, S.D.; Siahpush, M.; Mulhollen, A. Socioeconomic, Rural-Urban, and Racial Inequalities in US Cancer Mortality: Part I-All Cancers and Lung Cancer and Part II-Colorectal, Prostate, Breast, and Cervical Cancers. J. Cancer Epidemiol. 2011, 2011, 107497. [Google Scholar] [CrossRef]
  21. Henry, K.A.; Niu, X.; Boscoe, F.P. Geographic disparities in colorectal cancer survival. Int. J. Health Geogr. 2009, 8, 48. [Google Scholar] [CrossRef] [Green Version]
  22. Henry, K.A.; Sherman, R.L.; McDonald, K.; Johnson, C.J.; Lin, G.; Stroup, A.M.; Boscoe, F.P. Associations of census-tract poverty with subsite-specific colorectal cancer incidence rates and stage of disease at diagnosis in the United States. J. Cancer Epidemiol. 2014, 2014, 823484. [Google Scholar] [CrossRef] [PubMed]
  23. Lian, M.; Schootman, M.; Doubeni, C.A.; Park, Y.; Major, J.M.; Stone, R.A.; Laiyemo, A.O.; Hollenbeck, A.R.; Graubard, B.I.; Schatzkin, A. Geographic variation in colorectal cancer survival and the role of small-area socioeconomic deprivation: A multilevel survival analysis of the NIH-AARP Diet and Health Study Cohort. Am. J. Epidemiol. 2011, 174, 828–838. [Google Scholar] [CrossRef] [PubMed]
  24. Wiese, D.; Stroup, A.M.; Maiti, A.; Harris, G.; Lynch, S.M.; Vucetic, S.; Henry, K.A. Socioeconomic Disparities in Colon Cancer Survival: Revisiting Neighborhood Poverty using Residential Histories. Epidemiology 2020, 31, 728–735. [Google Scholar] [CrossRef] [PubMed]
  25. Wiese, D.; Stroup, A.M.; Maiti, A.; Harris, G.; Lynch, S.M.; Vucetic, S.; Henry, K.A. Residential mobility and geospatial disparities in colon cancer survival. Cancer Epidemiol. Biomark. Prev. 2020, 29, 2119–2125. [Google Scholar] [CrossRef] [PubMed]
  26. Niu, X.; Pawlish, K.S.; Roche, L.M. Cancer survival disparities by race/ethnicity and socioeconomic status in New Jersey. J. Health Care Poor Underserved 2010, 21, 144–160. [Google Scholar] [CrossRef] [PubMed]
  27. Krieger, N.; Quesenberry, C.; Peng, T.; Horn-Ross, P.; Stewart, S.; Brown, S.; Swallen, K.; Guillermo, T.; Suh, D.; Alvarez-Martinez, L. Social class, race/ethnicity, and incidence of breast, cervix, colon, lung, and prostate cancer among Asian, Black, Hispanic, and White residents of the San Francisco Bay Area, 1988–1992 (United States). Cancer Cause Control 1999, 10, 525–537. [Google Scholar] [CrossRef]
  28. Wang, L.; Wilson, S.E.; Stewart, D.B.; Hollenbeak, C.S. Marital status and colon cancer outcomes in US Surveillance, Epidemiology and End Results registries: Does marriage affect cancer survival by gender and stage? Cancer Epidemiol. 2011, 35, 417–422. [Google Scholar] [CrossRef] [PubMed]
  29. Yang, C.-C.; Cheng, L.-C.; Lin, Y.-W.; Wang, S.-C.; Ke, T.-M.; Huang, C.-I.; Su, Y.-C.; Tai, M.-H. The impact of marital status on survival in patients with surgically treated colon cancer. Medicine 2019, 98, e14856. [Google Scholar] [CrossRef]
  30. James, P.; Hart, J.E.; Banay, R.F.; Laden, F. Exposure to greenness and mortality in a nationwide prospective cohort study of women. Environ. Health Perspect. 2016, 124, 1344–1352. [Google Scholar] [CrossRef] [Green Version]
  31. Mitchell, R.; Popham, F. Effect of exposure to natural environment on health inequalities: An observational population study. Lancet 2008, 372, 1655–1660. [Google Scholar] [CrossRef] [Green Version]
  32. Richardson, E.; Pearce, J.; Mitchell, R.; Day, P.; Kingham, S. The association between green space and cause-specific mortality in urban New Zealand: An ecological analysis of green space utility. BMC Public Health 2010, 10, 240. [Google Scholar] [CrossRef] [Green Version]
  33. Richardson, E.A.; Mitchell, R. Gender differences in relationships between urban green space and health in the United Kingdom. Soc. Sci. Med. 2010, 71, 568–575. [Google Scholar] [CrossRef] [Green Version]
  34. Richardson, E.A.; Mitchell, R.; Hartig, T.; De Vries, S.; Astell-Burt, T.; Frumkin, H. Green cities and health: A question of scale? J. Epidemiol. Community Health 2012, 66, 160–165. [Google Scholar] [CrossRef] [Green Version]
  35. Keegan, T.H.M.; Shariff-Marco, S.; Sangaramoorthy, M.; Koo, J.; Hertz, A.; Schupp, C.W.; Yang, J.; John, E.M.; Gomez, S.L. Neighborhood influences on recreational physical activity and survival after breast cancer. Cancer Cause Control 2014, 25, 1295–1308. [Google Scholar] [CrossRef] [Green Version]
  36. Mears, M.; Brindley, P.; Jorgensen, A.; Maheswaran, R. Population-level linkages between urban greenspace and health inequality: The case for using multiple indicators of neighbourhood greenspace. Health Place 2020, 62, 102284. [Google Scholar] [CrossRef] [PubMed]
  37. Bratman, G.N.; Anderson, C.B.; Berman, M.G.; Cochran, B.; De Vries, S.; Flanders, J.; Folke, C.; Frumkin, H.; Gross, J.J.; Hartig, T. Nature and mental health: An ecosystem service perspective. Sci. Adv. 2019, 5, eaax0903. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Kondo, M.C.; Fluehr, J.M.; McKeon, T.; Branas, C.C. Urban green space and its impact on human health. Int. J. Environ. Res. Public Health 2018, 15, 445. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Mears, M.; Brindley, P.; Jorgensen, A.; Ersoy, E.; Maheswaran, R. Greenspace spatial characteristics and human health in an urban environment: An epidemiological study using landscape metrics in Sheffield, UK. Ecol. Indic. 2019, 106, 105464. [Google Scholar] [CrossRef]
  40. Sugiyama, T.; Carver, A.; Koohsari, M.J.; Veitch, J. Advantages of public green spaces in enhancing population health. Landsc. Urban Plan. 2018, 178, 12–17. [Google Scholar] [CrossRef]
  41. Van Dillen, S.M.; de Vries, S.; Groenewegen, P.P.; Spreeuwenberg, P. Greenspace in urban neighbourhoods and residents’ health: Adding quality to quantity. J. Epidemiol. Community Health 2012, 66, e8. [Google Scholar] [CrossRef] [Green Version]
  42. Mears, M.; Brindley, P. Measuring urban greenspace distribution equity: The importance of appropriate methodological approaches. ISPRS Int. J. Geo-Inf. 2019, 8, 286. [Google Scholar] [CrossRef] [Green Version]
  43. Hunter, M.C.R.; Brown, D.G. Spatial contagion: Gardening along the street in residential neighborhoods. Landsc. Urban Plan. 2012, 105, 407–416. [Google Scholar] [CrossRef]
  44. Gomez, S.L.; Shariff-Marco, S.; DeRouen, M.; Keegan, T.H.M.; Yen, I.H.; Mujahid, M.; Satariano, W.A.; Glaser, S.L. The impact of neighborhood social and built environment factors across the cancer continuum: Current research, methodological considerations, and future directions. Cancer 2015, 121, 2314–2330. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Demoury, C.; Thierry, B.; Richard, H.; Sigler, B.; Kestens, Y.; Parent, M.-E. Residential greenness and risk of prostate cancer: A case-control study in Montreal, Canada. Environ. Int. 2017, 98, 129–136. [Google Scholar] [CrossRef]
  46. Ashing-Giwa, K.T.; Lim, J.-W.; Tang, J. Surviving cervical cancer: Does health-related quality of life influence survival? Gynecol. Oncol. 2010, 118, 35–42. [Google Scholar] [CrossRef]
  47. Cunningham, R.; Sarfati, D.; Stanley, J.; Peterson, D.; Collings, S. Cancer survival in the context of mental illness: A national cohort study. Gen. Hosp. Psychiatry 2015, 37, 501–506. [Google Scholar] [CrossRef] [Green Version]
  48. McCormack, G.R.; Cabaj, J.; Orpana, H.; Lukic, R.; Blackstaffe, A.; Goopy, S.; Hagel, B.; Keough, N.; Martinson, R.; Chapman, J. Evidence synthesis A scoping review on the relations between urban form and health: A focus on Canadian quantitative evidence. Health Promot. Chronic Dis. Prev. Can. Res. Policy Pract. 2019, 39, 187. [Google Scholar] [CrossRef]
  49. Plascak, J.J.; Rundle, A.G.; Babel, R.A.; Llanos, A.A.; LaBelle, C.M.; Stroup, A.M.; Mooney, S.J. Drop-and-spin virtual neighborhood auditing: Assessing built environment for linkage to health studies. Am. J. Prev. Med. 2020, 58, 152–160. [Google Scholar] [CrossRef]
  50. Rzotkiewicz, A.; Pearson, A.L.; Dougherty, B.V.; Shortridge, A.; Wilson, N. Systematic review of the use of Google Street View in health research: Major themes, strengths, weaknesses and possibilities for future research. Health Place 2018, 52, 240–246. [Google Scholar] [CrossRef] [PubMed]
  51. Maharana, A.; Nsoesie, E.O. Use of Deep Learning to Examine the Association of the Built Environment With Prevalence of Neighborhood Adult ObesityDeep Learning to Examine the Built Environment and Neighborhood Adult Obesity PrevalenceDeep Learning to Examine the Built Environment and Neighborhood Adult Obesity Prevalence. JAMA Netw. Open 2018, 1, e181535. [Google Scholar] [CrossRef] [PubMed]
  52. Wendelboe-Nelson, C.; Kelly, S.; Kennedy, M.; Cherrie, J.W. A scoping review mapping research on green space and associated mental health benefits. Int. J. Environ. Res. Public Health 2019, 16, 2081. [Google Scholar] [CrossRef] [Green Version]
  53. Turner, M.G.; Gardner, R.H.; O’neill, R.V.; O’Neill, R.V. Landscape Ecology in Theory and Practice; Springer: Berlin/Heidelberg, Germany, 2001; Volume 401. [Google Scholar]
  54. Uuemaa, E.; Antrop, M.; Roosaare, J.; Marja, R.; Mander, Ü. Landscape metrics and indices: An overview of their use in landscape research. Living Rev. Landsc. Res. 2009, 3, 1–28. [Google Scholar] [CrossRef]
  55. Uuemaa, E.; Mander, Ü.; Marja, R. Trends in the use of landscape spatial metrics as landscape indicators: A review. Ecol. Indic. 2013, 28, 100–106. [Google Scholar] [CrossRef]
  56. NJSCR. NJ State Cancer Registry. Available online: https://www.state.nj.us/health/ces/reporting-entities/njscr/ (accessed on 28 April 2021).
  57. Percy, C.; Holten, V.V.; Muir, C.S. International Classification of Diseases for Oncology; World Health Organization: Geneva, Switzerland, 1990. [Google Scholar]
  58. Hurley, S.; Hertz, A.; Nelson, D.O.; Layefsky, M.; Von Behren, J.; Bernstein, L.; Deapen, D.; Reynolds, P. Tracing a path to the past: Exploring the use of commercial credit reporting data to construct residential histories for epidemiologic studies of environmental exposures. Am. J. Epidemiol. 2017, 185, 238–246. [Google Scholar] [CrossRef] [Green Version]
  59. Wheeler, D.C.; Wang, A. Assessment of Residential History Generation Using a Public-Record Database. Int. J. Environ. Res. Public Health 2015, 12, 11670–11682. [Google Scholar] [CrossRef] [Green Version]
  60. Jacquez, G.M.; Slotnick, M.J.; Meliker, J.R.; AvRuskin, G.; Copeland, G.; Nriagu, J. Accuracy of commercially available residential histories for epidemiologic studies. Am. J. Epidemiol. 2010, 173, 236–243. [Google Scholar] [CrossRef] [Green Version]
  61. Stinchcomb, D.; Roeser, A. NCI/SEER Residential History Project Technical Report; Westat, Inc.: Rockville, MD, USA, 2016. [Google Scholar]
  62. Texas A &M University. NAACCR Geocoder Data Dictionary; NAACCR: Springfield, IL, USA, 2016; Available online: http://www.naaccr.org/wp-content/uploads/2016/11/NAACCR-Geocoder-Data-Dictionary.pdf (accessed on 28 April 2021).
  63. Homer, C.H.; Fry, J.A.; Barnes, C.A. The national land cover database. US Geol. Surv. Fact Sheet 2012, 3020, 1–4. [Google Scholar]
  64. Hijmans, R.J.; van Etten, J.; Cheng, M.; Mattiuzzi, M.; Sumner, J.A.; Greenberg, O.P.; Lamigueiro, A.; Bevan, E.B. Racine, and A. Shortridge. Package ‘Raster’; R Foundation for Statistical Computing: Vienna, Austria, 2016. [Google Scholar]
  65. Akpinar, A.; Barbosa-Leiker, C.; Brooks, K.R. Does green space matter? Exploring relationships between green space type and health indicators. Urban For. Urban Green. 2016, 20, 407–418. [Google Scholar] [CrossRef]
  66. Wu, J.; Jackson, L. Inverse relationship between urban green space and childhood autism in California elementary school districts. Environ. Int. 2017, 107, 140–146. [Google Scholar] [CrossRef]
  67. Callaghan, C.T.; Bino, G.; Major, R.E.; Martin, J.M.; Lyons, M.B.; Kingsford, R.T. Heterogeneous urban green areas are bird diversity hotspots: Insights using continental-scale citizen science data. Landsc. Ecol. 2019, 34, 1231–1246. [Google Scholar] [CrossRef]
  68. Team, R.C. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2013. [Google Scholar]
  69. Hesselbarth, M.H.; Sciaini, M.; With, K.A.; Wiegand, K.; Nowosad, J. Landscapemetrics: An open-source R tool to calculate landscape metrics. Ecography 2019, 42, 1648–1657. [Google Scholar] [CrossRef] [Green Version]
  70. Van Der Wal, J.; Falconi, L.; Januchowski, S.; Shoo, L.; Storlie, C.; Van Der Wal, M.J. Package ‘SDMTools’; R Foundation for Statistical Computing: Vienna, Austria, 2014. [Google Scholar]
  71. Riitters, K.H.; O’neill, R.; Hunsaker, C.; Wickham, J.D.; Yankee, D.; Timmins, S.; Jones, K.; Jackson, B. A factor analysis of landscape pattern and structure metrics. Landsc. Ecol. 1995, 10, 23–39. [Google Scholar] [CrossRef]
  72. Li, H.; Reynolds, J. On definition and quantification of heterogeneity. Oikos 1995, 73, 280–284. [Google Scholar] [CrossRef]
  73. Li, H.; Reynolds, J.F. A simulation experiment to quantify spatial heterogeneity in categorical maps. Ecology 1994, 75, 2446–2455. [Google Scholar] [CrossRef]
  74. Tsai, W.-L.; McHale, M.R.; Jennings, V.; Marquet, O.; Hipp, J.A.; Leung, Y.-F.; Floyd, M.F. Relationships between characteristics of urban green land cover and mental health in US metropolitan areas. Int. J. Environ. Res. Public Health 2018, 15, 340. [Google Scholar] [CrossRef] [Green Version]
  75. Greenfield, E.J.; Nowak, D.J.; Walton, J.T. Assessment of 2001 NLCD percent tree and impervious cover estimates. Photogramm. Eng. Remote Sens. 2009, 75, 1279–1286. [Google Scholar] [CrossRef]
  76. Nowak, D.J.; Greenfield, E.J. Evaluating the National Land Cover Database tree canopy and impervious cover estimates across the conterminous United States: A comparison with photo-interpreted estimates. Environ. Manag. 2010, 46, 378–390. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  77. Howden-Chapman, P. Housing standards: A glossary of housing and health. J. Epidemiol. Community Health 2004, 58, 162–168. [Google Scholar] [CrossRef] [Green Version]
  78. Shaw, M. Housing and public health. Annu. Rev. Public Health 2004, 25, 397–418. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  79. Rauh, V.A.; Landrigan, P.J.; Claudio, L. Housing and health: Intersection of poverty and environmental exposures. Ann. N. Y. Acad. Sci. 2008, 1136, 276–288. [Google Scholar] [CrossRef] [PubMed]
  80. Lynch, S.M.; Mitra, N.; Ross, M.; Newcomb, C.; Dailey, K.; Jackson, T.; Zeigler-Johnson, C.M.; Riethman, H.; Branas, C.C.; Rebbeck, T.R. A Neighborhood-Wide Association Study (NWAS): Example of prostate cancer aggressiveness. PLoS ONE 2017, 12, e0174548. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  81. Zhang, Z.; Reinikainen, J.; Adeleke, K.A.; Pieterse, M.E.; Groothuis-Oudshoorn, C.G.M. Time-varying covariates and coefficients in Cox regression models. Ann. Transl. Med. 2018, 6, 121. [Google Scholar] [CrossRef] [Green Version]
  82. Mills, M. Introducing Survival and Event History Analysis; SAGE: London, UK, 2010. [Google Scholar]
  83. Zwiener, I.; Blettner, M.; Hommel, G. Survival analysis: Part 15 of a series on evaluation of scientific publications. Dtsch. Arztebl. Int. 2011, 108, 163–169. [Google Scholar] [CrossRef] [PubMed]
  84. Therneau, T.M.; Lumley, T. Package ‘survival’. R Top. Doc. 2015, 128, 112. [Google Scholar]
  85. Rappaport, J. Moving to nice weather. Reg. Sci. Urban Econ. 2007, 37, 375–398. [Google Scholar] [CrossRef]
  86. Banzhaf, H.S.; Walsh, R.P. Do people vote with their feet? An empirical test of Tiebout. Am. Econ. Rev. 2008, 98, 843–863. [Google Scholar] [CrossRef] [Green Version]
  87. Litwak, E.; Longino, C.F., Jr. Migration patterns among the elderly: A developmental perspective. Gerontol. 1987, 27, 266–272. [Google Scholar] [CrossRef]
  88. Iyer, H.S.; Valeri, L.; James, P.; Chen, J.T.; Hart, J.E.; Laden, F.; Holmes, M.D.; Rebbeck, T.R. The contribution of residential greenness to mortality among men with prostate cancer: A registry-based cohort study of Black and White men. Environ. Epidemiol. 2020, 4, e087. [Google Scholar] [CrossRef] [Green Version]
  89. O’Callaghan-Gordo, C.; Pollan, M.; Molina de la Torre, A.J.; Dierssen-Sotos, T.; Ardanaz, E.; Fernández-Tardón, G.; Amiano, P.; Capelo, R.; Chirlaque, M.-D.; Nieuwenhuijsen, M. Green Spaces and Colon Cancer Risk: A Case-Control in Spain. In Proceedings of the ISEE Conference Abstracts, Ottawa, ON, Canada, 26–30 August 2018. [Google Scholar]
  90. Schmitz, K.H.; Stout, N.L.; Maitin-Shepard, M.; Campbell, A.; Schwartz, A.L.; Grimmett, C.; Meyerhardt, J.A.; Sokolof, J.M. Moving through cancer: Setting the agenda to make exercise standard in oncology practice. Cancer 2021, 127, 476–484. [Google Scholar] [CrossRef] [PubMed]
  91. Brown, J.C.; Winters-Stone, K.; Lee, A.; Schmitz, K.H. Cancer, physical activity, and exercise. Compr. Physiol. 2012, 2, 2775–2809. [Google Scholar] [CrossRef] [Green Version]
  92. Van Blarigan, E.L.; Meyerhardt, J.A. Role of physical activity and diet after colorectal cancer diagnosis. J. Clin. Oncol. 2015, 33, 1825–1834. [Google Scholar] [CrossRef] [Green Version]
  93. Garvin, E.; Branas, C.; Keddem, S.; Sellman, J.; Cannuscio, C. More Than Just An Eyesore: Local Insights And Solutions on Vacant Land And Urban Health. J. Urban Health 2013, 90, 412–426. [Google Scholar] [CrossRef] [Green Version]
  94. Chen, Y.; Stephens, M.; Jones, C.A. Does residents’ satisfaction with the neighbourhood environment relate to residents’ self-rated health? Evidence from Beijing. Int. J. Environ. Res. Public Health 2019, 16, 5051. [Google Scholar] [CrossRef] [Green Version]
  95. Zuniga, T.; Adriana, A. From Neighborhoods to Wellbeing and Conservation: Enhancing the Use of Greenspace through Walkability; The University of Arizona: Tucson, AZ, USA, 2015. [Google Scholar]
  96. Jennings, V.; Baptiste, A.K.; Jelks, O.; Skeete, R. Urban green space and the pursuit of health equity in parts of the United States. Int. J. Environ. Res. Public Health 2017, 14, 1432. [Google Scholar] [CrossRef] [Green Version]
  97. Soderstrom, M. The Walkable City: From Haussmann’s Boulevards to Jane Jacobs’ Streets and Beyond; Vehicule Press: Montreal, QC, Canada, 2008. [Google Scholar]
  98. Doubeni, C.A.; Major, J.M.; Laiyemo, A.O.; Schootman, M.; Zauber, A.G.; Hollenbeck, A.R.; Sinha, R.; Allison, J. Contribution of Behavioral Risk Factors and Obesity to Socioeconomic Differences in Colorectal Cancer Incidence. JNCI J. Natl. Cancer Inst. 2012, 104, 1353–1362. [Google Scholar] [CrossRef] [Green Version]
  99. Wen, M.; Zhang, X.; Harris, C.D.; Holt, J.B.; Croft, J.B. Spatial Disparities in the Distribution of Parks and Green Spaces in the USA. Ann. Behav. Med. 2013, 45, S18–S27. [Google Scholar] [CrossRef] [Green Version]
Ijerph 18 04728 g001 550
Figure 1. Land cover classification and landscape characteristics of New Jersey: (a) Reclassified land cover classification (1 = Water, 2 = Sand/Bare, 3 = Forest, 4 = Shrubs, 5 = Grass, 6 = Open, 7 = Low, 8 = Medium, 9 = High intensity developed lands); (b) Census tract poverty level by categories (0–5%, 5–10%, 10–20%, 20+%); (c) Shannon Diversity Index (SHDI) by quartiles low to high; (d) Proportion of open developed lands by quartiles low to high; (e) Proportion of high-intensity developed lands by quartiles, low to high; (f) Forest contiguity index by quartiles, low to high; (g) Aggregation index (AI) of open developed lands by quartiles, low to high; (h) Aggregation index (AI) of medium-intensity developed lands by quartiles, low to high; (i) Aggregation index (AI) of high-intensity developed lands by quartiles, low to high.
Figure 1. Land cover classification and landscape characteristics of New Jersey: (a) Reclassified land cover classification (1 = Water, 2 = Sand/Bare, 3 = Forest, 4 = Shrubs, 5 = Grass, 6 = Open, 7 = Low, 8 = Medium, 9 = High intensity developed lands); (b) Census tract poverty level by categories (0–5%, 5–10%, 10–20%, 20+%); (c) Shannon Diversity Index (SHDI) by quartiles low to high; (d) Proportion of open developed lands by quartiles low to high; (e) Proportion of high-intensity developed lands by quartiles, low to high; (f) Forest contiguity index by quartiles, low to high; (g) Aggregation index (AI) of open developed lands by quartiles, low to high; (h) Aggregation index (AI) of medium-intensity developed lands by quartiles, low to high; (i) Aggregation index (AI) of high-intensity developed lands by quartiles, low to high.
Ijerph 18 04728 g001
Ijerph 18 04728 g002 550
Figure 2. Workflow diagram.
Figure 2. Workflow diagram.
Ijerph 18 04728 g002
Ijerph 18 04728 g003 550
Figure 3. Correlation matrix before and after reduction of significant variables: (A) All significant variables based on the univariate models; (B) Reduced set of significant variables used in final multivariate models. (Note: Land cover classes are based on the reclassified NLCD raster: U21 = Open, U22 = Low, U23 = Medium, U24 = High intensity developed lands). AI = Aggregation Index, SHDI = Shannon Diversity Index, CI = Contiguity Index, PRD = Patch Richness Density.
Figure 3. Correlation matrix before and after reduction of significant variables: (A) All significant variables based on the univariate models; (B) Reduced set of significant variables used in final multivariate models. (Note: Land cover classes are based on the reclassified NLCD raster: U21 = Open, U22 = Low, U23 = Medium, U24 = High intensity developed lands). AI = Aggregation Index, SHDI = Shannon Diversity Index, CI = Contiguity Index, PRD = Patch Richness Density.
Ijerph 18 04728 g003
Ijerph 18 04728 g004 550
Figure 4. Distribution of colon cancer patients by neighborhood/landscape characteristics.
Figure 4. Distribution of colon cancer patients by neighborhood/landscape characteristics.
Ijerph 18 04728 g004
Ijerph 18 04728 g005 550
Figure 5. Hazard Ratios of neighborhood-level variables from the crude models, adjusted for each variable individually. Blue line indicates there is no difference in risk. Red bars indicate the 95% confidence interval of statistically significant variables. (Note: Significance levels <0.05 = *; <0.01 = **; <0.001 = *** Land cover classes are based on the reclassified NLCD raster. U21 = Open, U22 = Low, U23 = Medium, U24 = High intensity developed lands).
Figure 5. Hazard Ratios of neighborhood-level variables from the crude models, adjusted for each variable individually. Blue line indicates there is no difference in risk. Red bars indicate the 95% confidence interval of statistically significant variables. (Note: Significance levels <0.05 = *; <0.01 = **; <0.001 = *** Land cover classes are based on the reclassified NLCD raster. U21 = Open, U22 = Low, U23 = Medium, U24 = High intensity developed lands).
Ijerph 18 04728 g005
Ijerph 18 04728 g006 550
Figure 6. Hazard Ratios of neighborhood-level variables from the multivariate model adjusted for all individual-level factors and each neighborhood/landscape variable individually. Blue line indicates a zero-coefficient or no difference in risk. Red bars indicate the 95% confidence interval of statistically significant variables. (Note: Significance levels <0.05 = *; <0.01 = **; Land cover classes are based on the reclassified NLCD raster. U21 = Open, U22 = Low, U23 = Medium, U24 = High intensity developed lands).
Figure 6. Hazard Ratios of neighborhood-level variables from the multivariate model adjusted for all individual-level factors and each neighborhood/landscape variable individually. Blue line indicates a zero-coefficient or no difference in risk. Red bars indicate the 95% confidence interval of statistically significant variables. (Note: Significance levels <0.05 = *; <0.01 = **; Land cover classes are based on the reclassified NLCD raster. U21 = Open, U22 = Low, U23 = Medium, U24 = High intensity developed lands).
Ijerph 18 04728 g006
Table
Table 1. Overview of the land cover classes before and after the reclassification of the NLCD dataset.
Table 1. Overview of the land cover classes before and after the reclassification of the NLCD dataset.
The NLCD CodeReclassificationCommentary
21: Developed Open SpaceU21: Developed Open SpaceMax 20% imperviousness per cell/pixel
22: Developed Low IntensityU22: Developed Low Intensity20–49% imperviousness per cell/pixel
23: Developed Medium IntensityU23: Developed Medium Intensity50–79% imperviousness per cell/pixel
24: Developed High IntensityU24: Developed High Intensity80–100% imperviousness per cell/pixel
41: Deciduous ForestForest Dominated by tree canopy and includes any type of parks and squares
42: Evergreen Forest
43: Mixed Forest
52: Shrub/ScrubShrubsDominated by shrubs; present on empty housing parcels
71: Grasslands/HerbaceousGrasslandIn urban areas, may assume a low-quality green space
81: Pasture/Hay
82: Cultivated Crops
90: Woody WetlandForestWoody wetlands are common in southern New Jersey and have large proportions of deciduous trees
95: Emergent Herbaceous WetlandGrasslandHerbaceous (also grassy) wetlands are typical for many coastal regions.
Table
Table 2. Overview of the area-based and individual variables.
Table 2. Overview of the area-based and individual variables.
VariablesLand CoversDefinitionCommentary
Land Cover Class-Level Metrics
Class ProportionForest, Grass, Shrubs, Industrial, Developed Lands (Open, Low, Medium, High Intensities)Composition metric. Proportional coverage—% of the landscape covered by each type.Used by [72,73] and recommended by [36] for green and water space; [74] recommended for Forest, Shrubs and Grass
Aggregation Index (AI)Developed Lands (Open, Low, Medium, High Intensities), Industrial AreasConfiguration metric. Computed as an area-weighted mean class aggregation index, where each class is weighted by its proportional area in the landscape.Redundant with several other metrics of proportion, cohesion, and contiguity and may be a meaningful alternative [53]
Splitting IndexForest, Grass, ShrubsConfiguration metric. A large splitting index, results from land covers being split into many patches with an even size distribution.Correlated with the aggregation index.
Applied for green spaces only in relation to health outcomes [39].
Contiguity Index (CI)Developed Lands (Open, Low, Medium, High Intensities), Forest, GrassConfiguration metric. Large contiguous patches will result in larger contiguity index values.CI for green/tree land cover classes associated with health outcomes [39].
Landscape-Level Metrics
Shannon Diversity IndexBased on all Land Cover ClassesComposition metric. The more classes and the more equally distributed, the higher the index.Used for measuring the aesthetic value and diversity [8].
Associated with health outcomes [39].
Patch Richness Density (PRD)Based on all Land Cover ClassesNumber of patches per hectare. High values indicate high dispersionPRD for green areas and recreational lands associated with poor health [39].
Contagion IndexBased on all Land Cover ClassesComposition metric. High values indicate result from landscapes with a few large, contiguous patches and low dispersion and interspersion of patch types
Average Proportion of ImperviousnessCensus Tract Average based on NLCD dataset estimating imperviousness proportion per pixelComposition metric. Highly negatively correlated with Tree Canopy proportions but is more accurateHighly correlated (negative) with Tree Canopy Cover but more accurate [75,76]
Census-Based Variables
Poverty Level by CategoryPercentage of population 18 and older living below federal poverty level.Socio-economic statusAssociated with cancer outcomes including survival and mortality
Median Year Structures BuiltMedian that the areas residential buildings were constructed.Organized into categories.Housing age and conditions are associated with health outcomes [17,77] and poverty [78]
Housing DensityNumber of structures per area unit (acre)Continuous variable defined by census tractPotential intermediate factor in health outcomes [79]
Table
Table 3. Study population characteristics.
Table 3. Study population characteristics.
Overall (n = 3949)
Age
Mean (SD)65.8 (13.3)
Median [Min, Max]68.0 [21.0, 85.0]
Gender
Male1878 (47.6%)
Female2071 (52.4%)
Race/Ethnicity
NH-White2902 (73.5%)
NH-Black488 (12.3%)
Hispanic (any race)325 (8.2%)
NH-API141 (3.6%)
Other93 (2.4%)
Regional Stage Subcategory
Regional, direct extension only1339 (33.9%)
Regional, lymph nodes only1268 (32.1%)
Regional, both1342 (34.0%)
Vital Status
Censored2862 (72.5%)
Colon Cancer Death1087 (27.5%)
Survival Time (months)
Mean (SD)62.3 (38.0)
Median [Min, Max]66.0 [1.00, 139]
CT Changes (Type of “moves”)
CT at Date of Diagnosis Only2587 (65.5%)
Change in Residential CT within NJ885 (22.4%)
Change in Residential CT outside NJ477 (12.1%)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wiese, D.; Stroup, A.M.; Maiti, A.; Harris, G.; Lynch, S.M.; Vucetic, S.; Gutierrez-Velez, V.H.; Henry, K.A. Measuring Neighborhood Landscapes: Associations between a Neighborhood’s Landscape Characteristics and Colon Cancer Survival. Int. J. Environ. Res. Public Health 2021, 18, 4728. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph18094728

AMA Style

Wiese D, Stroup AM, Maiti A, Harris G, Lynch SM, Vucetic S, Gutierrez-Velez VH, Henry KA. Measuring Neighborhood Landscapes: Associations between a Neighborhood’s Landscape Characteristics and Colon Cancer Survival. International Journal of Environmental Research and Public Health. 2021; 18(9):4728. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph18094728

Chicago/Turabian Style

Wiese, Daniel, Antoinette M. Stroup, Aniruddha Maiti, Gerald Harris, Shannon M. Lynch, Slobodan Vucetic, Victor H. Gutierrez-Velez, and Kevin A. Henry. 2021. "Measuring Neighborhood Landscapes: Associations between a Neighborhood’s Landscape Characteristics and Colon Cancer Survival" International Journal of Environmental Research and Public Health 18, no. 9: 4728. https://0-doi-org.brum.beds.ac.uk/10.3390/ijerph18094728

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop