Next Article in Journal
Do Drinking Motives Mediate the Relationship between Neighborhood Characteristics and Alcohol Use among Adolescents?
Previous Article in Journal
Is It Possible to Change from a Linear to a Circular Economy? An Overview of Opportunities and Barriers for European Small and Medium-Sized Enterprise Companies
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spatiotemporal Contextual Uncertainties in Green Space Exposure Measures: Exploring a Time Series of the Normalized Difference Vegetation Indices

Department of Human Geography and Spatial Planning, Faculty of Geosciences, Utrecht University, Princetonlaan 8a, 3584 CB Utrecht, The Netherlands
Int. J. Environ. Res. Public Health 2019, 16(5), 852; https://doi.org/10.3390/ijerph16050852
Submission received: 10 February 2019 / Revised: 4 March 2019 / Accepted: 5 March 2019 / Published: 8 March 2019
(This article belongs to the Section Environmental Health)

Abstract

:
Environmental health studies on green space may be affected by contextual uncertainties originating from the temporality of environmental exposures and by how the spatial context is delimitated. The Normalized Difference Vegetation Index (NDVI) is frequently used as an outdoor green space metric capturing the chlorophyll content in the vegetation canopy. This study assessed (1) whether residential NDVI exposures vary over time, and (2) how these time series of NDVI scores vary across spatial context delimitations. Multi-temporal NDVI data for the period 2006–2017 for the Netherlands were obtained from the Moderate Resolution Imaging Spectroradiometer (MODIS) satellite platform. Annual NDVI exposures were determined across multiple buffer sizes (i.e., 300, 600, and 1000 m) centered on a random sample of 10,000 Dutch residential addresses. Besides the descriptive statistics, pairwise Wilcoxon tests and Fligner–Killeen tests were used to determine mean and variance differences in annual NDVI scores across buffer widths. Heat maps visualized the correlation matrices. Significance levels were adjusted for multiple hypotheses testing. The results indicated that annual NDVI metrics were significantly correlated but their magnitude varied notably between 0.60 to 0.97. Numerous mean and variance differences in annual NDVI exposures were significant. It seems that the disparate buffers (i.e., 300 and 1000 m) were less strongly correlated, possibly because variance heterogeneity is reduced in larger buffers. These results have been largely consistent over the years and have passed Monte Carlo-based sensitivity tests. In conclusion, besides assessing green space exposures along different buffer sizes, our findings suggest that green space–health studies should employ NDVI data that are well-aligned with epidemiological data. Even an annual temporal incompatibility may obscure or distort green space–health associations. Both strategies may diminish contextual uncertainties in environmental exposure assessments.

Graphical Abstract

1. Introduction

It is increasingly recognized that the human body responds to the environments to which it is exposed [1,2,3]. The assessment of health-threatening and health-promoting environmental exposures has received fresh impetus from progress in earth observation coupled with geographic information systems (GIS) leading to high-quality environmental information [4,5]. To date, satellite images represent the gold standard for capturing conditions of the natural environment (i.e., vegetation) in an objective manner [5,6,7].
One prominent environmental factor is green space—namely forests, parks, etc.—which is capable of reducing harm caused by exposure to physical and mental health stressors [2,8]. Studies found that being surrounded by outdoor residential green space has beneficial effects for diverse population groups and on multiple health outcomes, such as adiposity [9], birth outcomes [10], suicide [11], physical activity [12], children’s brain development and cognitive function [13], the elderly’s mortality [14], and depression [15]. Pathways connecting green space and health include the mitigation of harmful exposures (e.g., air pollution, noise), the promotion of physical activity, and the support of restoring capacities [8].
Vegetation metrics [7] derived from remote sensing images such as the Normalized Difference Vegetation Index (NDVI) [16,17] are widely used to describe the distribution and amount of green space in population-based epidemiological studies [14,18,19]. However, findings that relate NDVI to health outcomes are not always consistent. For example, while [13] reported beneficial effects of green space on brain development, another study found an insignificant direct effect of residential NDVI extracted from the Moderate Resolution Imaging Spectroradiometer (MODIS) on major depressive disorders and psychological distress [20].
While there are manifold reasons for these contradictory findings (e.g., different research designs, study site specific differences), contextual uncertainties also contribute, to a certain extent, to the heterogeneity of findings on green space and health outcomes [21]. First, the temporality of environmental exposures contributes to contextual uncertainties. Green space varies naturally over time, but ambiguities can also arise artificially as a result of applied data processing [22]. Despite such temporal variations in NDVI scores, some studies [20,23] defend the application of NDVI data that are misaligned with epidemiological data. However, empirical results supporting such practices are lacking. Since longitudinal studies incorporating people’s life courses and their exposure to (temporally weighted) time series of green space are receiving increasing attention [9,24], well-aligned time series of NDVI data are crucial to circumvent exposure misclassifications.
Spatial contextual uncertainties arise due to the diverse ways in which the spatial context that people belong to is defined. Several studies have shown that assigning exposures based on administrative areas (e.g., postal codes) leads to conceptual and methodological flaws (for details, see [3]). Refined green space assessments can be achieved by focusing only on the vicinity of an individual’s actual residential location as represented through GIS-based buffers. As the “true” health influencing spatial context is unknown, the ways to determine an adequate buffer size are still under debate. Consensus has been reached, however, that different buffers sizes need to be explored to validate the stability of health–NDVI associations [25,26].
Motivated by these challenges, this study addressed spatiotemporal contextual uncertainties in exposure to outdoor NDVI for a large sample of addresses in the Netherlands. The first objective was to examine whether NDVI exposure measures vary significantly over time. Such analyses empirically underpin common empirical practices whereby temporally misaligned indicators are used as substitutes in assessing health–green space correlations. The second objective was to examine how variations in the spatial context—here defined as buffers with different sizes centered on residential addresses—lead to temporally misaligned exposure measures.
The empirical findings contribute to a better understanding of contextual uncertainties about NDVI scores. As spatial and temporal contextual uncertainties jointly affect green space assessments, this study puts both scale issues and the temporal alignment of green space exposures centrally. The results may assist those carrying out future studies to decide whether misaligned vegetation indices could serve as substitutes. To the best of our knowledge, no previous study has addressed these elementary but important issues, which have the potential to result in biased green space exposure assessments.

2. Materials and Methods

2.1. Study Area

This analysis was conducted in the Netherlands as part of the NEEDS study (www.needs.sites.uu.nl). The Netherlands was selected as a case study due to its diverse geography, which ranges from rural areas with much green space in the north of the country, to one of Europe’s most populated and urbanized areas in the middle of the country (i.e., the Randstad area).

2.2. Address Data

Like previous studies [23,27], the analysis was carried out at an address level. Address locations stem from the Dutch cadaster (version September 2017), which contains 9.1 million addresses. We extracted those for residential function and mixed use. From these 8.7 million addresses, a random sample of 10,000 locations was drawn; each had the same selection probability. Such a sample size corresponds to epidemiological studies. Figure A1 in the Appendix shows the distribution of the sampled addresses.

2.3. Remote Sensing-Based NDVI Data

The widely implemented satellite-derived Normalized Difference Vegetation Index (NDVI) was used to measure the level of greenness [17]. Such a multispectral vegetation index reflects the fact that chlorophyll absorbs the red light (i.e., 0.6–0.7 µm) and the mesophyll leaf structure scatters the near-infrared (0.7–1.1 µm) [16,28,29]. The NDVI is composed of the near-infrared radiation minus visible radiation divided by near-infrared radiation plus visible radiation [17]. The NDVI ranges between −1 and +1, where negative values refer to non-biomass (e.g., clouds, water, and snow) and values around 0 represent rocks and bare soil. The higher the positive values, the higher the amounts of biomass. To reduce the likelihood of biased exposure determination, values <0 (representing water bodies) are filtered for outdoor green space assessments [23].
We focused on MODIS NDVI satellite-based measures due to their global coverage as pre-processed time series while being less computationally intensive than Landsat or other higher resolution data. Annual MODIS NDVI 16-day composites (June 25–July 10, with the product ID MOD13Q1) in the hierarchical data format were acquired for the Netherlands spanning the period 2006–2017 [7,29]. Within the 16-day acquisition period, the algorithm to pre-process the MODIS data selected pixels with the lowest cloud coverage, a low view angle, as well as having the highest NDVI values [29]. Data were obtained through the online MODIS data repository (https://modis.gsfc.nasa.gov/data/dataprod/mod13.php). As vegetation follows a seasonal cycle with the highest amount of green during the summer, we consistently used annual data for June/July. To match with other available data, the original spatial resolution of 250 meters was down-scaled to 100 meters using bilinear interpolation. After clipping and re-projecting the data to the national Dutch reference system, the data were entered into a multilayer image stack, and the recommended scaling factor of 0.0001 was applied [29]. Figure 1 illustrates an example of the NDVI time series for Houten and its surroundings.

2.4. GIS-Based Exposure Assessment

To represent the spatial environmental context, we utilized a buffer approach using GIS [25,27]. Circular buffers with radii of 300, 600, and 1000 m were centered on each sampled address location. Smaller buffer sizes reflect the immediate surroundings, while 1000 m buffers capture the wider surroundings. In addition, such circular buffers allow for the exploration of scale effects and are more suitable for long-term analyses because they are independent of auxiliary data (i.e., the road network), which also change over time but are rarely available. Next, we determined the average NDVI values for every address across the three buffer sizes for each year. These proxy measures for green space exposure are attached to each location.

2.5. Statistical Analyses

Methods included descriptive statistics describing the central tendency including the mean, standard deviation, etc., to summarize the residential NDVI scores. Individual NDVI trajectories for each address and year were displayed through time series plots. Non-parametric correlations after Spearman coupled with a heat map [30] were implemented to assess associations across the different buffer sizes as well as their correlations over time. Pairwise Wilcoxon rank sum tests were used to test whether NDVI scores differ statistically over the years and across different buffer sizes [31]. Non-parametric Fligner–Killeen tests, which are robust against non-normally distributed values, were conducted to test whether the variances across the annual NDVI scores and buffer sizes differ [32].
Since we conducted a large number of tests, the multiple testing problem arose, namely the probability of false positives increased (i.e., finding differences by chance). To handle multiple hypotheses testing in the applied statistics, the p-values were adjusted following Holm’s procedure [33]. Finally, because the sample characteristics could vary by chance, sensitivity tests were performed. A Monte Carlo procedure [34] was used to examine whether the initial sample mean varied significantly compared to 100 alternative random samples, each 1000 addresses large, drawn from all residential addresses (i.e., the initial 10,000 addresses were excluded). The annual mean and standard deviation (SD) of the resulting reference distribution was compared to the initial sample per year. The analyses were carried out in the R environment for statistical computing [35].

3. Results

Table 1 together with Figure 2 summarize the annual NDVI scores aggregated across the sample. The boxplots confirm what the measures of central tendency indicated, namely that the NDVI values fluctuate slightly between 2006 and 2017. For example, the NDVI values for 2010 seem to be lower compared to the other years. There seems to be a minor overall temporal trend of increasing NDVI scores since 2006. Regression lines added to the boxplots confirm this impression. Similar patterns were observed when cross-comparing the different buffer sizes.
Because such aggregated analyses hide individual temporal patterns, the individual time series of the average 16-days NDVI composite values per address location were also examined. Figure 3 displays the time series plots. While the overall plots do not indicate striking patterns, an individual investigation showed some striking temporal changes of the NDVI scores per location.
The core of the analyses was, as summarized in Figure 4, assessing the NDVI correlations over time and across buffer sizes. While all correlations were highly statistically significant (p < 0.05), even after adjusting for multiple testing, some patterns were recognizable. For example, the 1000 m buffers were less strongly correlated than the 600 and 300 m buffers. Further, it seemed that some years were less associated, although still statistically significant. The years 2007 and 2016 were, independent of the examined buffer size, less associated with the other years of NDVI data. This was in line with Figure 1, whereas deviations in the NDVI patterns were recognizable. Even annual differences of a single year may result in some correlation differences.
The results of the pairwise mean comparisons between years with corrections for multiple testing are given in Table 2. As indicated by the p-values, several years of our sampled NDVI values have different means, independent of which buffer size is considered. Moreover, the Fligner–Killeen tests to assess the homogeneity of variances rejected the null hypothesis of equal variances over the years and buffer sizes (all p < 0.001). Our final sensitivity tests, using 100 random samples of 1000 address samples and comparing mean NDVI scores, indicated no statistically significant differences. The initial mean NDVI score was well-situated in the simulated distribution of mean NDVI values plus/minus one standard deviation. Figure A2 in the Appendix illustrates the result for NDVI values for 2014 and the 600 m buffer.

4. Discussion

Accurately assessing the level of outdoor green space surrounding people’s residential homes is crucial for epidemiological studies. Due to its cost-effectivity and large coverage of high-quality time series, the NDVI is commonly used for regional or even nationwide studies [14,19]. To limit contextual uncertainties—a fundamental methodological problem [21]—it is important to have NDVI data that are temporally well-aligned with the epidemiological data, although it is not common practice [20,23]. Such data mismatch may translate into biased outcomes, that is, lead to an over- or underestimation of the NDVI effect size in regressions. This study shed light on the often ignored circumstance of temporally incompatible green space data by exploring annual time series of NDVI exposures derived from buffers centered on a large sample of residential addresses randomly distributed across the Netherlands.

4.1. Interpretation of the Results

The results showed that, as expected, individual years of NDVI scores were significantly correlated; however, a few years were less strongly correlated than others. These differences were further confirmed statistically through mean and variance differences in NDVI exposures. Our findings question the practice of having a temporally poorly-aligned linkage of green space and epidemiological data. For example, Rugel et al. [20] used survey data for 2012 while the NDVI was for 2014–2016. At least partly due to temporally incompatible data, an insignificant green space effect on mental health outcomes was reported by that study. It could be that the level of green space changed, which translated into less accurate exposure assessments. However, the NDVI differences could also emerge due to mixed pixels, misregistration of the satellite location, etc. (see [16] for a discussion).
There is no consensus on the buffer size to delineate the residential context [25,26]. Our tested buffer widths followed the literature, reflecting the closer and extended residential environment. The results show that the NDVIs for different buffer widths over time are highly positively correlated, ranging from 0.60 to 0.97. A similar range of correlations for MODIS NDVIs was reported elsewhere [26]. Furthermore, there seems to be a tendency for the disparate buffers to be less correlated, which is in concordance with a study in New York City [26]. An explanation could be that smaller buffers have a pronounced variance while NDVI differences average out in larger buffers. Roughly, these results are consistent over the time series of NDVI values. Utilizing a Monte Carlo approach, these results passed sensitivity tests. It can be ruled out that the initial sample represented an extreme case. Given these results, it is recommended not to use temporally incompatible NDVI and epidemiological survey data. To reduce temporal uncertainties, the date should be aligned as closely as possible with the survey data. Without having conducted the obligatory tests, we speculate that discrepancies may be attenuated when cumulative green space exposures over time (or time-weighted ones), rather than purely cross-sectional metrics, are employed. Here, systematic investigations across different study designs, health outcomes, study sites, etc. require further attention.
Because MODIS NDVI maps are retrieved from daily images to generate 16-day composite measures, exact temporal matches are possible for most parts of the world, which is challenging for other popular sensor platforms such as Landsat. The promising temporal granularity of MODIS images is challenged by a moderate spatial resolution of 250 meters and a minimum mapping size of detectable vegetation. This may be problematic for urban areas where small-scale green spaces exist, which may lead to an over- or underestimation of the true green space exposure. Besides these technological constraints, it is debatable whether green space measures from downward-facing satellites (including all platforms) represent accurately how people on the ground perceive green space [36,37]. When comparing neighborhood-based NDVI scores and street view-based green space for Beijing (China), no significant associations were found across both metrics [37]. Such overhead-view assessments only capture the available vegetation; they do not provide insights into the quality, accessibility, and actual use of these areas [20,38].
The application of green space metrics is more complex than it seems at first sight due not only to different sensor platforms (e.g., MODIS, Landsat) but also to complications ensuing from the variety of measures proposed. While the NDVI correlates with expert ratings [39], the enhanced vegetation index (i.e., an alternative vegetation index), for example, is designed to overcome methodological limitations of the NDVI including the elimination of atmospheric and canopy background effects contaminating standard NDVI while having better sensitivity at higher vegetation levels [7,29] (for an application, see [20]). Others implemented classified land use information [15,27,40]—which is challenging because of a varying number of land use categories or sensitivity to the underlying image-based classification algorithm—or used measures such as the leaf area index [41] and the vegetation continuous fields tree cover [42,43]. However, our results are of more general interest because the MODIS NDVI was found to be highly correlated with alternative vegetation measures such as the leaf area index [41].

4.2. Study Limitations

There are numerous limitations to this analysis. Firstly, the present analytical study design is restricted. While imitating real-word ego-centric exposure assessments, the study did not evaluate the actual consequences of misaligned green space measures on physical and mental health outcomes across different settings (e.g., countries, cities and rural areas, spatial contexts). Studies that rule out incompatible green space data and that do not lead to biased regression coefficients are advised. Secondly, although down-scaling increased the number of cells per buffer while limiting the influence of edge effects, we cannot exclude that this procedure artificially over- or underestimated the average NDVI values. Furthermore, when interpreting the results, it is critical to stress that all NDVI values given, within the threshold buffer distance, contribute equally (i.e., mean of NDVI cells per buffer). A refinement could be distance-based weighting schemes because cells close to homes are more important than those further away. Thirdly, notwithstanding that NDVI was found to be fairly constant across satellite platforms [6], only MODIS NDVI measures were explored. Thus, the results should not be uncritically transferred to other platforms (e.g., RapidEye, Sentinel 2). Such high-resolution sensors support more precise inner-urban analyses. Fourthly, for the sake of simplicity, we assumed the residential home location to be the source of exposure, as did others [19,23]. This conventional conceptualization leaves aside exposures along people’s day-to-day mobility, which is relevant for more accurate exposure assessments [3,21]. Finally, as only a large, representative sample for the Netherlands was investigated, the generalization of findings could be questioned.

5. Conclusions

This paper addressed contextual uncertainties originating from both the temporality of environmental exposures and how the spatial context is delimitated. The empirical results show that annual NDVI scores derived from different buffer sizes differed significantly in their means and variances. Furthermore, within the range of 0.60 to 0.97, the magnitude of the correlations varied notably over the years and buffer sizes.
These findings are vital for exposure assessments that often enrich epidemiological data with temporally incompatible NDVI data, thus potentially obscuring green space–health associations. Even a misalignment of a single year may distort correlations. To mitigate contextual uncertainties, it is advised to integrate temporally well-aligned green space data. While the present study took a first step toward an empirical basis addressing spatiotemporal contextual uncertainties, future studies should systematically research the consequences for correlation studies.

Funding

This research is part of the NEEDS project (Dynamic Urban Environmental Exposures on Depression and Suicide, http://needs.sites.uu.nl/), which has received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation program (grant agreement No 714993).

Acknowledgments

I appreciate the reviewers’ critical and constructive comments on the manuscripts, which contributed to the quality of this publication. I thank the IJERPH for providing a waiver for the article processing charge.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Appendix A

Figure A1. Sampled address locations. The red dot in the center indicates the location shown in Figure 1.
Figure A1. Sampled address locations. The red dot in the center indicates the location shown in Figure 1.
Ijerph 16 00852 g0a1
Figure A2. Simulated reference distribution of mean NDVI values for 2014 using a 600 m buffer based on 100 samples each 1000 comprising addresses. The dotted line refers to the mean NDVI value of the initial sample.
Figure A2. Simulated reference distribution of mean NDVI values for 2014 using a 600 m buffer based on 100 samples each 1000 comprising addresses. The dotted line refers to the mean NDVI value of the initial sample.
Ijerph 16 00852 g0a2

References

  1. Frumkin, H. Environmental Health: From Global to Local; John Wiley & Sons: New York, NY, USA, 2016. [Google Scholar]
  2. Van den Bosch, M.; Sang, Å. Urban natural environments as nature-based solutions for improved public health—A systematic review of reviews. Environ. Res. 2017, 158, 373–384. [Google Scholar] [CrossRef] [PubMed]
  3. Helbich, M. Toward dynamic urban environmental exposure assessments in mental health research. Environ. Res. 2018, 161, 129–135. [Google Scholar] [CrossRef] [PubMed]
  4. Toth, C.; Jóźków, G. Remote sensing platforms and sensors: A survey. ISPRS J. Photogramm. Remote Sens. 2016, 115, 22–36. [Google Scholar] [CrossRef]
  5. Viana, J.; Santos, J.V.; Neiva, R.M.; Souza, J.; Duarte, L.; Teodoro, A.C.; Freitas, A. Remote Sensing in Human Health: A 10-Year Bibliometric Analysis. Remote Sens. 2017, 9, 1225. [Google Scholar] [CrossRef]
  6. Brown, M.E.; Pinzón, J.E.; Didan, K.; Morisette, J.T.; Tucker, C.J. Evaluation of the consistency of long-term NDVI time series derived from AVHRR, SPOT-vegetation, SeaWiFS, MODIS, and Landsat ETM+ sensors. IEEE Trans. Geosci. Remote Sens. 2006, 44, 1787–1793. [Google Scholar] [CrossRef]
  7. Huete, A.; Didan, K.; van Leeuwen, W.; Miura, T.; Glenn, E. MODIS vegetation indices. In Land Remote Sensing and Global Environmental Change; Springer: New York, NY, USA, 2010; pp. 579–602. [Google Scholar]
  8. Hartig, T.; Mitchell, R.; De Vries, S.; Frumkin, H. Nature and health. Annu. Rev. Public Health 2014, 35, 207–228. [Google Scholar] [CrossRef]
  9. Persson Åand Pyko, A.; Lind, T.; Bellander, T.; Östenson, C.-G.; Pershagen, G.; Eriksson, C.; Lõhmus, M. Urban residential greenness and adiposity: A cohort study in Stockholm County. Environ. Int. 2018, 121, 832–841. [Google Scholar] [CrossRef]
  10. Ebisu, K.; Holford, T.R.; Bell, M.L. Association between greenness, urbanicity, and birth weight. Sci. Total Environ. 2016, 542, 750–756. [Google Scholar] [CrossRef]
  11. Helbich, M.; de Beurs, D.; Kwan, M.-P.; O’Connor, R.C.; Groenewegen, P.P. Natural environments and suicide mortality in the Netherlands: A cross-sectional, ecological study. Lancet Planet. Health 2018, 2, e134–e139. [Google Scholar] [CrossRef]
  12. Chong, S.; Byun, R.; Mazumdar, S.; Bauman, A.; Jalaludin, B. Effects of Distant Green Space on Physical Activity in Sydney, Australia. J. Phys. Act. Health 2017, 14, 29–35. [Google Scholar] [CrossRef]
  13. Dadvand, P.; Basagaña, X.; Jerrett, M.; Nieuwenhuijsen, M.J.; Sunyer Deu, J.; Pujol, J.; Macià, D.; Martinez Vilavella, G.; Blanco Hinojo, L.; Mortamais, M.; et al. The association between lifelong greenspace exposure and 3-dimensional brain magnetic resonance imaging in Barcelona schoolchildren. Environ. Health Perspect. 2018, 126, 027012. [Google Scholar] [CrossRef]
  14. Ji, J.S.; Zhu, A.; Bai, C.; Wu, C.-D.; Yan, L.; Tang, S.; Zeng, Y.; James, P. Residential greenness and mortality in oldest-old women and men in China: A longitudinal cohort study. Lancet Planet. Health 2019, 3, e17–e25. [Google Scholar] [CrossRef]
  15. Helbich, M.; Klein, N.; Roberts, H.; Hagedoorn, P.; Groenewegen, P.P. More green space is related to less antidepressant prescription rates in the Netherlands: A Bayesian geoadditive quantile regression approach. Environ. Res. 2018, 166, 290–297. [Google Scholar] [CrossRef]
  16. Pettorelli, N.; Vik, J.O.; Mysterud, A.; Gaillard, J.-M.; Tucker, C.J.; Stenseth, N.C. Using the satellite-derived NDVI to assess ecological responses to environmental change. Trends Ecol. Evol. 2005, 20, 503–510. [Google Scholar] [CrossRef] [PubMed]
  17. Tucker, C.J. Red and photographic infrared linear combinations for monitoring vegetation. Remote Sens. Environ. 1979, 8, 127–150. [Google Scholar] [CrossRef] [Green Version]
  18. Balseviciene, B.; Sinkariova, L.; Grazuleviciene, R.; Andrusaityte, S.; Uzdanaviciute, I.; Dedele, A.; Nieuwenhuijsen, M.J. Impact of residential greenness on preschool children’s emotional and behavioral problems. Int. J. Environ. Res. Public Health 2014, 11, 6757–6770. [Google Scholar] [CrossRef] [PubMed]
  19. Song, H.; Lane, K.J.; Kim, H.; Kim, H.; Byun, G.; Le, M.; Choi, Y.; Park, C.R.; Lee, J.-T. Association between Urban Greenness and Depressive Symptoms: Evaluation of Greenness Using Various Indicators. Int. J. Environ. Res. Public Health 2019, 16, 173. [Google Scholar] [CrossRef]
  20. Rugel, E.J.; Carpiano, R.M.; Henderson, S.B.; Brauer, M. Exposure to natural space, sense of community belonging, and adverse mental health outcomes across an urban region. Environ. Res. 2019. [Google Scholar] [CrossRef] [PubMed]
  21. Kwan, M.-P. The Limits of the Neighborhood Effect: Contextual Uncertainties in Geographic, Environmental Health, and Social Science Research. Ann. Am. Assoc. Geogr. 2018, 108, 1482–1490. [Google Scholar] [CrossRef]
  22. Wenxia, G.; Huanfeng, S.; Liangpei, Z.; Wei, G. Normalization of NDVI from Different Sensor System using MODIS Products as Reference. In IOP Conference Series: Earth and Environmental Science; IOP Publishing: Bristol, UK, 2014; Volume 17, p. 12225. [Google Scholar]
  23. Gascon, M.; Sánchez-Benavides, G.; Dadvand, P.; Martinez, D.; Gramunt, N.; Gotsens, X.; Cirach, M.; Vert, C.; Molinuevo, J.L.; Crous-Bou, M.; et al. Long-term exposure to residential green and blue spaces and anxiety and depression in adults: A cross-sectional study. Environ. Res. 2018, 162, 231–239. [Google Scholar] [CrossRef]
  24. Pearce, J.; Shortt, N.; Rind, E.; Mitchell, R. Life course, green space and health: Incorporating place into life course epidemiology. Int. J. Environ. Res. Public Health 2016, 13, 331. [Google Scholar] [CrossRef] [PubMed]
  25. Frank, L.D.; Fox, E.H.; Ulmer, J.M.; Chapman, J.E.; Kershaw, S.E.; Sallis, J.F.; Conway, T.L.; Cerin, E.; Cain, K.L.; Adams, M.A.; et al. International comparison of observation-specific spatial buffers: Maximizing the ability to estimate physical activity. Int. J. Health Geogr. 2017, 16, 4. [Google Scholar] [CrossRef] [PubMed]
  26. Reid, C.E.; Kubzansky, L.D.; Li, J.; Shmool, J.L.; Clougherty, J.E. It’s not easy assessing greenness: A comparison of NDVI datasets and neighborhood types and their associations with self-rated health in New York City. Health Place 2018, 54, 92–101. [Google Scholar] [CrossRef] [PubMed]
  27. Boers, S.; Hagoort, K.; Scheepers, F.; Helbich, M. Does residential green and blue space promote recovery in psychotic disorders? A cross-sectional study in the province of Utrecht, the Netherlands. Int. J. Environ. Res. Public Health 2018, 15, 2195. [Google Scholar] [CrossRef]
  28. Myneni, R.B.; Hall, F.G.; Sellers, P.J.; Marshak, A.L. The interpretation of spectral vegetation indexes. IEEE Trans. Geosci. Remote Sens. 1995, 33, 481–486. [Google Scholar] [CrossRef]
  29. Didan, K.; Munoz, A.B.; Solano, R.; Huete, A. MODIS vegetation index user’s guide (MOD13 Series). 2015. Available online: https://vip.arizona.edu/documents/MODIS/MODIS_VI_UsersGuide_June_2015_C6.pdf (accessed on 10 February 2019).
  30. Zhao, S.; Guo, Y.; Sheng, Q.; Shyr, Y. Advanced heat map and clustering analysis using heatmap3. Biomed. Res. Int. 2014, 2014. [Google Scholar] [CrossRef]
  31. Kanji, G.K. 100 Statistical Tests; Sage: Thousand Oaks, CA, USA, 2006. [Google Scholar]
  32. Conover, W.J.; Johnson, M.E.; Johnson, M.M. A comparative study of tests for homogeneity of variances, with applications to the outer continental shelf bidding data. Technometrics 1981, 23, 351–361. [Google Scholar] [CrossRef]
  33. Holm, S. A simple sequentially rejective multiple test procedure. Scand. J. Stat. 1979, 65–70. [Google Scholar]
  34. Hope, A.C.A. A simplified Monte Carlo significance test procedure. J. R. Stat. Soc. Ser. B 1968, 582–598. [Google Scholar] [CrossRef]
  35. R Core Team. R: A Language and Environment for Statistical Computing; R Core Team: Vienna, Austria, 2018. [Google Scholar]
  36. Dong, R.; Zhang, Y.; Zhao, J. How Green Are the Streets Within the Sixth Ring Road of Beijing? An Analysis Based on Tencent Street View Pictures and the Green View Index. Int. J. Environ. Res. Public Health 2018, 15, 1367. [Google Scholar] [CrossRef]
  37. Helbich, M.; Yao, Y.; Liu, Y.; Zhang, J.; Liu, P.; Wang, R. Using deep learning to examine street view green and blue spaces and their associations with geriatric depression in Beijing, China. Environ. Int. 2019, 126, 107–117. [Google Scholar] [CrossRef] [PubMed]
  38. Ekkel, E.D.; de Vries, S. Nearby green space and human health: Evaluating accessibility metrics. Landsc. Urban Plan. 2017, 157, 214–220. [Google Scholar] [CrossRef]
  39. Rhew, I.C.; Vander Stoep, A.; Kearney, A.; Smith, N.L.; Dunbar, M.D. Validation of the normalized difference vegetation index as a measure of neighborhood greenness. Ann. Epidemiol. 2011, 21, 946–952. [Google Scholar] [CrossRef]
  40. Groenewegen, P.P.; Zock, J.-P.; Spreeuwenberg, P.; Helbich, M.; Hoek, G.; Ruijsbroek, A.; Strak, M.; Verheij, R.; Volker, B.; Waverijn, G.; et al. Neighbourhood social and physical environment and general practitioner assessed morbidity. Health Place 2018, 49, 68–84. [Google Scholar] [CrossRef] [PubMed]
  41. Chen, P.-Y.; Fedosejevs, G.; Tiscareno-Lopez, M.; Arnold, J.G. Assessment of MODIS-EVI, MODIS-NDVI and VEGETATION-NDVI composite data using agricultural measurements: An example at corn fields in western Mexico. Environ. Monit. Assess. 2006, 119, 69–82. [Google Scholar] [CrossRef] [PubMed]
  42. Dadvand, P.; Tischer, C.; Estarlich, M.; Llop, S.; Dalmau-Bueno, A.; López-Vicente, M.; Valentin, A.; de Keijzer, C.; Fernández-Somoano, A.; Lertxundi, N.; et al. Lifelong residential exposure to green space and attention: A population-based prospective study. Environ. Health Perspect. 2017, 125, 097016. [Google Scholar] [CrossRef]
  43. Sexton, J.O.; Song, X.-P.; Feng, M.; Noojipady, P.; Anand, A.; Huang, C.; Kim, D.-H.; Collins, K.M.; Channan, S.; DiMiceli, C.; et al. Global, 30-m resolution continuous fields of tree cover: Landsat-based rescaling of MODIS vegetation continuous fields with lidar-based estimates of error. Int. J. Digit. Earth 2013, 6, 427–448. [Google Scholar] [CrossRef] [Green Version]
1
© 2019 by the authors. Submitted for possible open access publication under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).
Figure 1. Example of the annual NDVI time series from 2006–2017 for the municipality of Houten, the Netherlands.
Figure 1. Example of the annual NDVI time series from 2006–2017 for the municipality of Houten, the Netherlands.
Ijerph 16 00852 g001
Figure 2. Boxplots of the NDVI scores per year across different buffer sizes. The black regression line refers to the trend over time.
Figure 2. Boxplots of the NDVI scores per year across different buffer sizes. The black regression line refers to the trend over time.
Ijerph 16 00852 g002
Figure 3. NDVI time series for the sampled addresses using different buffer radii. Each horizontal line refers to the NDVI trend of an address (N = 10,000).
Figure 3. NDVI time series for the sampled addresses using different buffer radii. Each horizontal line refers to the NDVI trend of an address (N = 10,000).
Ijerph 16 00852 g003
Figure 4. Heat map of the correlations for the NDVI time series for the 1000, 600, and 300 m buffers. All plotted correlation values were significant after adjusting for multiple testing.
Figure 4. Heat map of the correlations for the NDVI time series for the 1000, 600, and 300 m buffers. All plotted correlation values were significant after adjusting for multiple testing.
Ijerph 16 00852 g004
Table 1. Summary statistics for the residential NDVI scores (N = 10,000).
Table 1. Summary statistics for the residential NDVI scores (N = 10,000).
Buffer (m)YearMeanSDMedian25th Percentile75th PercentileMin.MaxSkewnessKurtosis
100020060.550.110.560.480.630.1970.84−0.1595−0.2427
20070.530.140.540.440.630.1190.88−0.3160−0.1868
20080.550.110.550.480.630.1800.85−0.2441−0.0829
20090.570.110.580.500.660.1850.85−0.2540−0.2093
20100.520.100.520.460.590.1700.82−0.1243−0.1466
20110.570.120.570.490.660.1900.86−0.1028−0.3292
20120.530.130.540.450.630.1330.84−0.3213−0.2992
20130.590.110.590.520.670.1980.87−0.2091−0.2376
20140.560.120.560.480.640.1630.85−0.2075−0.2629
20150.560.110.560.480.630.1820.86−0.1228−0.1823
20160.540.130.540.460.630.1580.89−0.1377−0.3493
20170.570.120.570.500.650.0880.88−0.30840.1690
60020060.540.110.540.470.610.1420.86−0.1232−0.0792
20070.520.140.530.420.620.0820.88−0.2706−0.2122
20080.540.110.540.470.610.1230.86−0.20670.1153
20090.560.110.560.490.640.1180.86−0.2223−0.0631
20100.510.100.510.450.580.0920.84−0.07080.0195
20110.560.120.550.480.630.1400.86−0.0320−0.1048
20120.520.130.530.430.620.1070.86−0.3209−0.2185
20130.570.110.570.500.650.1260.87−0.1449−0.0414
20140.550.120.550.470.630.1320.86−0.1682−0.0987
20150.540.110.540.470.620.1190.87−0.06270.0493
20160.530.130.530.440.620.1020.89−0.1643−0.2389
20170.560.120.560.490.630.0560.89−0.31250.4280
30020060.530.110.530.450.600.0860.87−0.09370.0042
20070.510.150.520.410.610.0240.88−0.2453−0.2379
20080.520.110.530.450.600.0550.86−0.18370.2192
20090.550.120.550.480.630.0870.87−0.21080.0401
20100.500.100.500.430.570.0620.85−0.02090.1092
20110.540.120.540.470.620.1050.870.02410.0100
20120.510.140.530.420.610.0690.86−0.3458−0.1663
20130.560.110.560.490.630.0740.87−0.09030.1060
20140.530.120.540.460.610.0820.87−0.15690.0586
20150.530.110.530.460.600.0720.88−0.00470.1354
20160.520.140.530.430.620.0450.93−0.2037−0.1898
20170.550.120.550.470.620.0500.89−0.31840.5647
Table 2. Significance levels for the Wilcoxon tests (p-values adjusted after Holm).
Table 2. Significance levels for the Wilcoxon tests (p-values adjusted after Holm).
Buffer (m)Year20062007200820092010201120122013201420152016
100020070
20080.0340
2009000
20100000
20110000.0670
201200.3310000
20130000000
20140.0060000000
20150.33100000000.331
20160000000.065000
20170000.02300.70800000
60020070
20080.1320
2009000
20100000
201100000
201200.1320000
20130000000
201400000000
20150.13200000000.064
20160.00100.1320000.015000
20170000.00800.13200000
30020070
20080.4370
2009000
20100000
201100000
201200.0580.027000
20130000000
2014000000.06700
20150.18100.006000000.004
20160.6800.680000.009000.027
20170000.0050000000

Share and Cite

MDPI and ACS Style

Helbich, M. Spatiotemporal Contextual Uncertainties in Green Space Exposure Measures: Exploring a Time Series of the Normalized Difference Vegetation Indices. Int. J. Environ. Res. Public Health 2019, 16, 852. https://doi.org/10.3390/ijerph16050852

AMA Style

Helbich M. Spatiotemporal Contextual Uncertainties in Green Space Exposure Measures: Exploring a Time Series of the Normalized Difference Vegetation Indices. International Journal of Environmental Research and Public Health. 2019; 16(5):852. https://doi.org/10.3390/ijerph16050852

Chicago/Turabian Style

Helbich, Marco. 2019. "Spatiotemporal Contextual Uncertainties in Green Space Exposure Measures: Exploring a Time Series of the Normalized Difference Vegetation Indices" International Journal of Environmental Research and Public Health 16, no. 5: 852. https://doi.org/10.3390/ijerph16050852

APA Style

Helbich, M. (2019). Spatiotemporal Contextual Uncertainties in Green Space Exposure Measures: Exploring a Time Series of the Normalized Difference Vegetation Indices. International Journal of Environmental Research and Public Health, 16(5), 852. https://doi.org/10.3390/ijerph16050852

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