Next Article in Journal
Drinking Water Consumption Patterns: An Exploration of Risk Perception and Governance in Two First Nations Communities
Next Article in Special Issue
Measurement of 30-Year Urban Expansion Using Spatial Entropy in Changwon and Gimhae, Korea
Previous Article in Journal
Design of Strategies for the Implementation and Management of a Complementary Monetary System Using the SWOT-AHP Methodology
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effects of Green Space Patterns on Urban Thermal Environment at Multiple Spatial–Temporal Scales

1
Institute of Remote Sensing and Earth Sciences, College of Science, Hangzhou Normal University, Hangzhou 311121, China
2
Zhejiang Provincial Key Laboratory of Urban Wetlands and Regional Change, Hangzhou 311121, China
3
College of Geomatics & Municipal Engineering, Zhejiang University of Water Resources and Electric Power, Hangzhou 310018, China
4
Department of Forestry and Natural Resources, Purdue University, West Lafayette, IN 47907, USA
*
Author to whom correspondence should be addressed.
Sustainability 2020, 12(17), 6850; https://0-doi-org.brum.beds.ac.uk/10.3390/su12176850
Submission received: 12 June 2020 / Revised: 23 July 2020 / Accepted: 28 July 2020 / Published: 24 August 2020

Abstract

:
Land use/land cover (LULC) pattern change due to human activity is one of the key components of regional and global climate change drivers. Urban green space plays a critical role in regulating urban thermal environment, and its cooling effect has received widespread attention in urban heat island (UHI) related studies. To fully understand the effects of the landscape pattern of an urban green space in regulating the urban thermal environment, it is necessary to further study the thermal effects of the landscape pattern of the urban green space and its characteristics under varied spatial–temporal scales. In this paper, we took the urban core area of Hangzhou City as the study area and analyzed the relationships between the landscape metrics of the urban green space and land surface temperature (LST) under varied spatial scales by using correlation analysis and redundancy analysis (RDA) methods. Multi-temporal Landsat 8 thermal infrared sensor data were used to retrieve the spatial and temporal dynamics of LSTs in four consecutive seasons, and the land use classification was interpreted using SPOT (Systeme Probatoire d’Observation de la Terre) satellite imagery. The results showed that landscape dominance metrics—e.g., percentage of landscape (PLAND) and largest patch index (LPI)—were the most influential factors on urban LST. The spatial configuration of urban landscape, as represented by the fragmentation and aggregation and connectedness, also showed significant effects on LST. Furthermore, landscape pattern metrics had varied spatial scale effects on LST. Specifically, the landscape dominance metrics of urban forest showed an increased influence on LST as the spatial scale increased, while for urban water, the trend was opposite. These findings might have some practical significance for urban planning about how to spatially arrange urban green space to alleviate UHI at local and regional scales.

1. Introduction

The land use/land cover (LULC) pattern change due to human activity is one of the key forces of regional and global climate change [1,2,3,4], and large-scale urbanization causing changes in land surface composition and structure during the past decades has been a main factor in LULC [5,6]. The expansion of urban land directly results in the relative reduction of natural vegetation and an increase of impervious surface. As a result, the changed evapotranspiration and thermal inertia [7,8,9] of the land surface unit, together with the heat produced by human activities during energy consumption [10,11], will increase the surface temperature of the urban area, and this is the so-called urban heat island (UHI) effect [12,13,14]. As the land surface temperature (LST) increases continuously in the context of global warming and urban areas continue to expand, especially in the developing countries and regions, there is an increased risk of environmental effects of UHI [15,16,17,18].
Urban green space is critical in regulating the urban thermal environment and improving the local microclimate of the urban area, and its cooling effect has received widespread attention in UHI-related studies [19,20,21]. Over the past decades, local-scale field observations [22,23], urban-scale remote sensing investigations [24,25,26], and model simulations [27,28,29] have been applied to study the relationship between urban green space and LST. At the single patch scale, the spatial characteristics—e.g., area, perimeter, and shape indices—of urban green space were found to be significantly correlated with LST [30,31,32]. At multi-patch scales—or, rather, the landscape scale—one of the common consensuses in UHI studies was that vegetation coverage or normalized vegetation index (NDVI) plays an important role in reducing LST [33,34]. The spatial layout of urban green patches and non-green patches is of great significance for the regulation of microclimate at urban or regional scales [35,36].
The relationship between landscape pattern and the process and function of a particular ecosystem under varied spatial scales is one of the key problems in landscape ecology [37,38]. The regulation of urban green space on the urban thermal environment by its spatial pattern has received increased attention in recent years. Numerous studies have tried to establish the relationship between landscape indices: e.g., percentage of landscape (PLAND), edge density (ED), path density (PD), landscape shape index (LSI), and aggregation index (AI) [39,40,41,42,43,44], and LST. It is acknowledged that the composition or configuration of green space could significantly affect LST. However, results from previous studies showed inconsistent or even contradictory relationships between LST and the spatial patterns of urban green space [45,46], and this suggested that the correlations between the landscape patterns of the urban green space and LST might be scale-dependent [47,48]. Meanwhile, previous studies have shown that there are seasonal variations in the relationships between LST and land cover features. Urban green spaces have been confirmed in forming cool islands and improving outdoor thermal comfort in summer both in field observation and remote sensing-related studies [21,49]. However, the growth conditions of green spaces differ significantly in different seasons, and the intensity of UHI has been found to change seasonally due to the high variability in vegetation growth [25,50]. Therefore, a comprehensive understanding of how the urban green spaces affect temperature distributions all year round is worth exploring. To fully understand the effects of the landscape pattern of an urban green space in regulating the urban thermal environment, it is necessary to further study the thermal effects of the landscape pattern of the urban green space and its characteristics under varied spatial–temporal scales.
In this study, we took the urban core area of Hangzhou City—one of the megacities in China—as the study area and explored the relationship between the landscape pattern of the urban green space and LST. We retrieved the LSTs in four consecutive seasons by Landsat 8 images and obtained the land use information from one SPOT images. We analyzed the relationships between the spatial patterns of the urban green space and LST with a correlation analysis and redundancy analysis (RDA). The objectives were to quantitatively address the following questions: (1) for typical urbanized areas, how does the green space affect the urban thermal environment across different seasons? (2) What are the most influential landscape pattern metrics correlated with urban LST? (3) Is there a spatial scale dependence between landscape metrics and urban LST?

2. Materials and Methods

2.1. Study Area

Hangzhou City, located South of the Yangtze River Delta (YRD) (Figure 1), has experienced rapid urbanization over the past decades [51]. Hangzhou is dominated by a monsoon climate, with an annual mean temperature ranging from 15.9 to 17.0 °C. The extreme maximum temperature was 39.8–42.9 °C, and the extreme minimum temperature was −7.1 to −15.0 °C. The annual mean precipitation of this city was 1378.5 mm, and the number of annual precipitation days were 134; the annual sunshine hours were 1471.2, and the average annual relative humidity was 71%. Hangzhou has experienced drastic urbanization over the past two decades. For example, the population has increased from 1.79 to 6.35 million, and the built-up area increased from 177.18 to 615.22 km2 during 2000–2018 in the urban district, respectively (Hangzhou Statistic Bureau, available at http://tjj.hangzhou.gov.cn/). Meanwhile, the contradiction between urban development and greenspace protection in Hangzhou is becoming more prominent [52,53]. We selected the urban core area bordered by the expressway around the city (about 800 km2) as the study area, which covered the entire downtown of Hangzhou City with the densest population.

2.2. Classification of Urban Land Use

Urban green spaces refer to all kinds of green infrastructures that use vegetated spaces and water bodies for the public well-being. Impervious surfaces are land surfaces, such as roads, parking lots, and building roofs, that repel rainwater and do not permit it to soak into the ground. It is noted that the study area is the core area of the city, so a large part of the land uses were built-up lands, mixed with urban forest composed by the mountainous area, wet land, parks, and trees along the roads. The cropland mainly distributed in the periphery of the urban core and its proportion was relatively small. We classified the land use in the study area into four Level I classes with 15 Level II sub-classes. The four Level I classes were Cropland (1), Forest (2), Water (3), and Built-up (4), and the corresponding 15 Level II sub-classes were Cropland (11), Street trees (21), Park trees (22), Mountain forest (23), Other forest (24), River (31), Lake (32), Pond (33), Residential (41), Industrial (42), Commercial (43), Administrative (44), Educational (45), Road (46), and Under construction (47). The first digit of the land classification code of the Level II class referred to its corresponding Level I class. The land use information in the study area was extracted from one cloudless SPOT 6 image obtained on 28 December, 2015. Imageries with higher spatial resolution could more accurately quantify the spatial pattern of urban land use, such as the smaller elements of greenspace, as well as changes in the configuration of larger patches [54]. The multispectral and panchromatic bands of the SPOT image were fused to get a 1.5 m resolution image using the Gram–Schmidt method. The SPOT image was radiometrically and atmospherically corrected, and the accuracy of the geometric correction was under 0.5 pixel. The pre-processed fused image was classified into 15 classes using the maximum likelihood classification method. The training sites of each Level II class were collected by field survey and visually identified polygons on the SPOT satellite imagery (R4, G3, and B2), and then each of the classes was assigned and recoded into four Level I classes, respectively (Figure 2). The post-classification was assisted by using a 1:10,000 land use map and visual interpretation under the aid of Google Earth satellite images. The classification accuracy was assessed by stratified random sampling with reference to ground-truth data collected in the field survey. The assessment results showed that the overall accuracy of land use classification was above 95% and the kappa coefficient was 0.93.

2.3. LST Retrieval

The Landsat satellite series with a thermal infrared band is the most common imagery source for LST retrieval. From the Landsat 4 Thematic Mapper to Landsat 8, several widely used algorithms have been developed for retrieving LST from the thermal channel of Landsat imageries: for example, the radiative transfer equation (RTE) and single channel algorithm (SCA) [55,56], the mono-window algorithm (MWA) [57,58], and the spit-window algorithm (SWA) [59].
LST observations acquired by remote sensing technologies have been used to assess the UHI, to develop land surface–atmosphere exchange models, and to analyze the relationship between temperature and land use and land cover in urban areas. The multi-temporal satellite data is effective and useful in diurnal, seasonal, and interannual urban land surface temperature variation studies [12,60]. In this study, we retrieved LSTs from four Landsat 8 imageries representing four consecutive seasons. The acquisition dates of the Landsat 8 images were 13 October, 2015, 16 December, 2015, 15 March, 2016, and 27 July, 2016, all at 10:31 a.m. local time. The 13 October and 16 December represented autumn and winter, the cold and dry seasons, respectively; the 15 March represented spring, the dry and warm season; and 27 July represented summer, the hot and humid season. The images were generally cloud free and the imaging dates were approximately in the middle of the seasons. These dates were no more than 8 months from the acquisition date of the SPOT image used for urban green space mapping, and it was reasonable to assume that there was no significant change in urban land uses [45]. To retrieve LST, Landsat 8 Thermal Infrared Sensor (TIRS) band 10 is more commonly used than band 11 due to its higher accuracy [61,62], and the RTE applied is:
L sensor = [ ε B ( T s ) + ( 1 ε ) L up ] τ λ + L down
where Lsensor is the radiance measured by the sensor, ε is the land surface emissivity, Ts is the land surface temperature according to B(Ts), which is the blackbody radiance given by the Planck’s law, τ is the total atmospheric transmissivity for channel λ between the surface and the sensor, and Lup and Ldown, are the downwelling and upwelling atmospheric radiance, respectively.
The atmospheric parameters τ, Lup and Ldown can be obtained by MODTRAN (MODerate resolution atmospheric TRANsmission) via the atmospheric correction parameter calculator provided by NASA (https://atmcorr.gsfc.nasa.gov/) [63]. The value of Lsensor can be directly converted from Landsat 8 TIRS band data by using the radiance rescaling factors provided in the metadata file of Landsat 8 data [62]. The value of ε was estimated by the NDVI thresholds method (NTM), which uses certain NDVI values to distinguish between soil pixels, full vegetation pixels, and mixed pixels composed of soil and vegetation [64]. The emissivities of water and built-up areas were adopted from the ASTER Global Emissivity Database 4.0 [58,65].
With the thermal radiance measured at sensor level, the atmospheric parameters obtained from the radiative transfer model and the land surface emissivity estimated using NTM, the Ts could be retrieved as follows [61]:
T s = K 2 ln ( K 1 L sensor L up τ ( 1 ε ) L down τ ε + 1 )
where K1 is 774.89 Wm−2sr−1μm−1, and K2 is 1321.08 K [66].

2.4. Landscape Pattern Metrics

The landscape metrics are quantitative representations of the spatial composition of a particular landscape pattern. In this study, we used four categories of landscape pattern metrics—i.e., dominance, shape complexity, aggregation and connectedness, and fragmentation [39,40,41,42,43,44]—to quantify the landscape pattern and explore the corresponding correlations with LST (Table 1 and Table S1 in the supplementary materials). The landscape pattern metrics were analyzed by the software FRAGSTATS 4.3.

2.5. Statistical Analysis

To explore the effects of urban green space patterns on LST, different grid scales were utilized. In the literature, the grid size generally ranges from 200 to 1000 m, or even 3000 m, according to the greenspace condition or the sampling density in a case-specific manner [40,42,68,69]. In this study, we repeatedly divided the study area into fishnet-like grids with different spatial scales (Figure 3)—i.e., 0.5 to 4.0 km every 500 m [44]—to study the relationships between the landscape metrics of the urban green space and LST under varied spatial scales. Since this study mainly concerned the influence of the landscape pattern of urban green space on LST, the sampling grids containing cropland parcels were excluded; in addition, the grids at the edge of the study area with blank areas of more than 20% of the total area of the specific grid were also excluded, including the grids with cloud or cloud shadow. As a result, the valid sampling grid numbers for the spatial scale ranging from 0.5–4.0 km every 500 m were 2000, 394, 156, 77, 41, 29, 23, and 19, respectively.
The Pearson correlation analysis and direct gradient analysis were used to analyze the relationships between urban LST and the landscape metrics of the urban green space. In environmental studies, one of the most common tasks is to identify the most influential environmental factors. Compared with a regression analysis, factor analysis, or principal analysis, the direct gradient analysis is special in its ability to quantify the unique variation belonging to specific explanatory variables [70,71,72]. In the direct gradient analysis, the detrended correspondence analysis (DCA) was used to detect gradient length in the data sets. In the case of a short gradient—i.e., the longest gradient length obtained with DCA was less than 3 standard deviations—the linear model with redundancy analysis (RDA) was applied; otherwise, the canonical correspondence analysis (CCA) was used. We applied DCA to LST and found that all gradient lengths were less than 3 standard deviations, so RDA was applied to explore the relationships among LST and landscape metrics: i.e., the species and environmental variables in the RDA context, respectively. Two important outputs could be obtained by RDA: i.e., the proportions (%) of the total variance of LST explained by each canonical axis, and the ordination diagrams, which explicitly reveal the relationships between landscape pattern metrics and LST. A positive correlation is expected when the arrows of two variables point to the same direction, and vice versa. The RDAs were performed using the CANOCO 4.5 program [73].
A comprehensive flowchart was designed to show the variables used in the statistical analysis and how they were combined with the spatial and temporal analysis (Figure 4). Firstly, the satellite imagery-based land use and LST data were retrieved from SPOT 6 and Landsat 8 images, respectively, and then the landscape pattern metrics and mean LST of each grid at different spatial scales were calculated under varied grid scales ranging from 0.5 to 4 km. Finally, the relationships between urban green space patterns and LST at multiple spatial–temporal scales were analyzed by using a correlation analysis and RDA.

3. Results

3.1. Spatial Pattern of LST across Different Seasons

There were remarkable seasonal and spatial variations in LST. Figure 5 shows the spatial distribution of LSTs in four seasons. The average LSTs were 24.93, 10.42, 25.71, and 34.58 °C in autumn (October), winter (December), spring (March), and summer (July), respectively (Table 2). The maximum standard deviation (SD) of LST, 4.23 °C, was in July, followed by 3.33, 2.37, and 1.32 °C in March, October, and December, respectively.
Table 3 shows the LST statistics of the different land uses in four seasons. The impervious surface generally had the highest mean LSTs in four seasons, except in December, 2015, when its mean LST was equal to the LST of urban water. The urban water had the lowest mean LSTs in all seasons except winter (December). The mean LSTs of urban forest were higher than those of urban water (except December), but lower than cropland. The largest LST variation was found in summer (July) when the impervious surface and urban forest had the largest temperature difference compared with the other seasons. In contrast, the mean LSTs of all the land uses were much closer in winter (December).

3.2. Relationships between Land Use Proportion and LST

The correlation analysis showed that the proportions of vegetation and water at grid level had significant negative correlations with LST on most spatiotemporal scales, except winter for water; for the impervious surface, significant positive correlations with LST could be found on all scales in each season (Figure 6). The absolute value of the correlation coefficients between forest proportion and LST was the largest in summer, followed by autumn, spring, and winter, and increased with the increase of grid scales. The water proportion was significantly negatively correlated with LST in spring, summer, and autumn, but significantly positively correlated with LST in winter at small scales: e.g., 0.5–2.0 km. The absolute values of the correlation coefficients decreased with the increase of grid scales. The positive correlation coefficients between impervious surface and LST increased with the increase of grid scales, and were larger in summer and autumn, followed by spring and winter.

3.3. Relationships between Urban Green Space Landscape Pattern and LST

In the Pearson correlation analysis, the correlation coefficients between land uses proportion (at grid scale) and LST increased with the increasing of the spatial scale (Figure 6). To simplify the analysis of the effects of landscape patterns on the spatial variation of urban LSTs, RDA was applied at 1, 2, and 4 km grid scales, respectively. A Monte Carlo permutation test (499 permutations) was used to determine the statistical validity of RDA, and the significant explanatory variables were identified. The explanations of landscape pattern metrics on LST at different spatial scales and the biplots of LST and significant landscape metrics according to RDA were shown in Table 4 and Figure 7, respectively.
At 1, 2, and 4 km grid scales, the 32.9%, 63.7%, and 99.8% variations in LST could be explained by the landscape patterns of urban forest, respectively; the corresponding percentages for urban water were 71.6%, 82.8%, and 99.3%, respectively; for the impervious surface, the percentages were 69.7%, 84.6%, and 98.8%, respectively. The landscape metrics of PLAND, LPI, area-weighted mean perimeter–area ratio (PAAM), landscape division index (DIVISION), PLAND, AI, cohesion index (COHESION), and area-weighted mean contiguity index (CONAM) showed relatively high explanatory capabilities at all three scales, in which PLAND, LPI, and DIVISION were the most influential explanatory variables. There were 4%, 8%, and 11% variations in LST at 1, 2, and 4 km grid scales, which could be explained by the landscape dominance metrics PLAND and LPI, respectively; for urban water, about 8% of variations in LST could be explained by PLAND and LPI at all three scales, but the explained variations slightly decreased as the spatial scale increased; for the impervious surface, all the explained variations in LST by PLAND and LPI were greater than 10% but decreased as the spatial scale increased. The metric DIVISION, which represented the landscape fragmentation of urban forest, urban water, and impervious surface, showed similar explanatory capabilities to LST variations, such as PLAND and LPI, at all three scales.
The landscape complexity metrics—except mean perimeter–area ratio (PAMN) and LSI—of urban water showed relative high explanatory capabilities on LST compared with those of urban forest and impervious surface, and the explained variations in LST increased as the spatial scale increased (Table 4). For example, the explained variations in LST by mean shape index (SHMN), area-weighted mean shape index (SHAM), mean fractal dimension index (FRMN), and area-weighted mean fractal dimension index (FRAM) of the urban water increased from 3.6%–6.7% at 1 km to 6.1%–8.1% and 7.3%–13.4% at 2 and 4 km spatial scales, respectively. For the aggregation and connectedness metrics, the explanatory capabilities of proportion of like adjacencies (PLADJ), AI, COHESION, and CONAM of the urban forest and impervious surface increased as the spatial scale increased; the explained variations in LST by COHESION of urban water showed an increasing trend as the spatial scale increased, while PLADJ, AI, and CONAM were the opposite.
Figure 7 shows the ordination diagrams of the RDA analysis of seasonal LST with respect to the landscape metrics of urban forest, urban water, and impervious surface at 1, 2, and 4 km grid scales, respectively. The first axis in the ordination diagrams represents changes of LST in three seasons (March, July, and October), and in one season (December) to a lesser extent. The arrows of two variables pointing to the same direction indicate a positive correlation, and the angle between two arrows is inversely proportional to the degree of their correlation. The LSTs in all seasons were negatively correlated with the landscape dominance, and the aggregation and connectedness metrics—e.g., the PLAND, LPI, PLADJ, AI, and clumpiness index (CLUMPY)—of urban forest at all scales, and were positively correlated with landscape fragmentation and complexity metrics: e.g., DIVISION and PAAM, at all scales. The correlations between LST and the landscape metrics of urban water were not consistent among seasons. In March, July, and October, urban water had similar effects on LST as those of urban forest; however, the correlations between LST in December and landscape dominance, aggregation and connectedness, fragmentation, and shape complexity—as expressed by the metrics of PLAND, LPI, PLADJ, AI, CLUMPY, DIVISION, and PAAM—of urban water were much smaller compared with the other seasons. The impervious surface showed opposite effects on LST compared with that of urban forest. The LSTs in all seasons were negatively correlated with the landscape shape complexity and fragmentation indices—e.g., PAAM and DIVISION—at all spatial scales, and negatively correlated with landscape dominance and aggregation and connectedness of the impervious surface: e.g., PLAND, LPI, and PLADJ.

4. Discussion

4.1. Seasonal Variability in the Relationships between Urban Green Space and LST

In this study, the mean LSTs of impervious surface in all seasons except winter (December) were higher than the other land uses (Table 3), and the correlation analysis further indicated that LST was positively correlated with the proportion of impervious surface in all seasons (Figure 6), while negatively correlated with that of urban forest; for urban water, the correlation with LST was negative except in winter. This was consistent with the general consensus that vegetation and water could alleviate UHI, while impervious surface was the opposite due to the different specific heat capacities [74,75,76,77,78]. Furthermore, we showed that the mean LSTs of different land uses differed mostly in July, when the mean LSTs of the impervious surface were 4.4 and 8.7 °C higher than the urban forest and urban water, respectively. The mean LSTs of different land uses were closest in December: for example, the LST difference between urban forest and impervious surface was only 1.3 °C; the LST differences between impervious surface and the other land uses in October and March were in between those of July and December (Table 3). This indicated that different land uses might have varied thermal environment effects.
Numerous studies have showed that there are seasonal differences in the cooling effect of urban green space—generally, the strongest is in summer while the weakest is in winter—and it might have a close relationship with the evapotranspiration of the urban forest [79,80]. This was consistent with the Pearson correlation analysis between LSTs and the proportion of urban forest on multiple spatial scales (Figure 6), in which the correlation coefficients in warmer seasons were much higher than those in winter. Previous studies also showed that the correlations between land uses and LST in winter were relatively weak compared with the other seasons [81,82]. The RAD analysis also showed that landscape pattern metrics could explain more variations in LST in July than December (Figure 7). This suggested that integrating urban forests with urban water, such as urban rivers, lakes, and ponds, may have great potential to improve human thermal comfort in urban areas. Additionally, although cropland is slightly weaker than that of urban water in regulating the thermal environment all year around, it should be strengthened in the future together with the optimal allocation and design of other urban green spaces as a sustainable food production system with economic and environmental benefits.

4.2. Influential Landscape Pattern Metrics on LST

The RDA analysis showed that landscape dominance metrics PLAND and LPI were the most influential factors on urban LST (Table 4). Generally, the explanatory capabilities of the landscape dominance of urban forest and water were higher than the other landscape pattern metrics, and this suggested that the cooling effects might be closely related with the enhanced evapotranspiration as a result of increased forest and water coverage at local scales. In the literature, the land use composition was found to be more influential on LST than its spatial configuration [41,43,83,84], and this was consistent with our findings. This means that increasing vegetation coverage and expanding the water area in urban areas are still the effective measures to increase surface cooling coverage.
However, we found that the spatial configuration of the urban landscape also showed significant effects on LST. Specifically, the explanatory capability of DIVISION, one metric of the landscape fragmentation, was almost comparable with that of LPI, and it meant that the more fragmented and dispersed the urban forest, including water area and contiguous sprawl of construction land, the higher the LST. Meanwhile, the PD and number of patches (NP) of urban forest showed significant positive correlations with LST, and this implied that the fragmentation of urban green space might be an important factor leading to the increase of LST; similar reports could be found in [41,42,43,85]. On the contrary, compared with DIVISION, the more the aggregation and connectedness of the urban forest and urban water, the better the cooling effect, while for the impervious surface, the trend was opposite.
In terms of implications for urban planning and land use management, it might still be necessary to increase the proportion of urban green spaces—just as with the common practices in urban green space planning—to relieve UHI. However, the focus should also be on optimizing the layout of urban green spaces to reduce LST due to the shortage of urban land [30,68,86]. The relationships of landscape fragmentation and aggregation and connectedness of urban green space to LST revealed in this study means that, under the premise that the area of urban green space remains unchanged, the larger and more concentrated and connected urban green spaces might be more effective in cooling land surface, rather than numerous small green space patches.

4.3. Scaling Effects of Urban Green Space on LST

The landscape pattern metrics could explain more than 60% variation in the urban thermal environment at all spatial scales except at the 1 km grid scale, in which 32.9% of variations in LST were explained (Table 4). This confirmed that the composition and configuration of urban landscape might be the key factors controlling the spatial variation of LST [68,69].
Meanwhile, the landscape pattern metrics showed varied spatial scale effects on LST. For example, the RDA analysis showed that the landscape dominance metrics (PLAND and LPI) of urban forest had an increased influence on LST as the spatial scale increased, while for urban water, the trend was opposite (Figure 7). The correlation analysis also revealed that the proportions of urban forest and urban water had opposite trends in correlation with LST as the spatial scale increased (Figure 6). Additionally, we found that a large part of landscape metrics—for example, the aggregation and connectedness ones of urban forest, urban water, and impervious surface, and the shape complexity ones of urban water—had increased effects on LST as the spatial scale increased. These findings might have some practical significance for urban planning about how to arrange urban green space to alleviate UHI. On a smaller block scale, for example, a reasonably increasing urban water area, if possible, and avoiding the continuous development of construction land might achieve a certain cooling effect, while on urban or regional scales, it might be more effective to alleviate UHI by optimizing the spatial pattern of urban forest, such as increasing the coverage and improving the agglomeration and connectivity of urban forest. Planning and managing urban green space and urban heat island mitigation needs to be approached holistically, taking into account diverse spatial–temporal dynamics [87,88].

5. Conclusions

There were remarkable seasonal and spatial variations in LST correlated with specific land uses. The impervious surface and urban water had the highest and lowest mean LSTs in all seasons except winter (December), respectively. The largest LST variation was found in summer (July) compared with the other seasons. In contrast, the mean LSTs of all the land uses were much closer in winter (December). In all seasons, the LSTs were positively correlated with the proportion of impervious surface, while negatively correlated with that of urban forest; for urban water, the correlation was negative, except in winter. The cooling effect of urban forest showed seasonal differences: generally, the strongest were in summer while the weakest were in winter.
Generally, the explanatory capabilities of landscape dominance metrics were higher than the other landscape pattern metrics. The spatial configuration of urban landscape, as represented by the fragmentation and aggregation and connectedness of the land use patches, also showed significant effects on LST. We suggested herein that it might still be necessary to increase the proportion of urban green spaces, just as with the common practices applied in urban green space planning, to relieve UHI; however, for cases where there is no more land allocated for urban green space, an optimized spatial arrangement of urban green space might also be an available means to counteract UHI.
The landscape pattern metrics had varied spatial scale effects on LST. Specifically, the landscape dominance metrics of urban forest showed that the influence on LST increased as the spatial scale increased, while for urban water, the trend was opposite. These findings might have some practical significance for urban planning about how to spatially arrange urban green space to alleviate UHI at local and regional scales.
Although apparently conclusive, there were limitations in interpreting the relationship between landscape pattern metrics and LST at multiple spatial scales. Numerous studies about the relationship between LST and urban green space are still case-specific. Further research should also consider the influential factors of urban morphology, such as the tree canopy cover, urban block form, and mountain topography, and comparison studies across different cities under varied climatic conditions will be of particular interest. Additionally, it would be helpful to establish a common set of landscape metrics to avoid redundancy and improve the comparability across similar studies.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2071-1050/12/17/6850/s1, Table S1: Landscape pattern metrics used in this study after McGarigal et al.

Author Contributions

Y.S. and X.S. conceived and designed the framework of the research structure. All authors contributed to the data analysis, interpretation, and manuscript writing. All authors have read and agreed to the published version of the manuscript.

Funding

This research is supported by National Natural Science Foundation of China (grant number: 41301484), Entrepreneurship and Innovation Project for High-level Overseas Returnees in Hangzhou City in 2019, and Science and Technology Project of the Water Resources Department of Zhejiang Province (grant number: RC1810).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Feddema, J.J. The Importance of Land-Cover Change in Simulating Future Climates. Science 2005, 310, 1674–1678. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Pielke, R.A.; Marland, G.; Betts, R.A.; Chase, T.N.; Running, S.W. The influence of land-use change and landscape dynamics on the climate system: Relevance to climate-change policy beyond the radiative effect of greenhouse gases. Philos. Trans. R. Soc. A 2002, 360, 1705–1719. [Google Scholar] [CrossRef] [PubMed]
  3. Kalnay, E.; Cai, M. Impact of urbanization and land-use change on climate. Nature 2003, 423, 528–531. [Google Scholar] [CrossRef] [PubMed]
  4. Meli, P.; Holl, K.D.; Benayas, J.M.R.; Jones, H.P.; Jones, P.C.; Montoya, D.; Mateos, D.M. A global review of past land use, climate, and active vs. passive restoration effects on forest recovery. PLoS ONE 2017, 12, e0171368. [Google Scholar]
  5. Deng, J.S.; Ke, W.; Yang, H.; Qi, J.G. Spatio-temporal dynamics and evolution of land use change and landscape pattern in response to rapid urbanization. Landsc. Urban Plan. 2009, 92, 187–198. [Google Scholar] [CrossRef]
  6. Dadashpoor, H.; Azizi, P.; Moghadasi, M. Land use change, urbanization, and change in landscape pattern in a metropolitan area. Sci. Total Environ. 2019, 655, 707–719. [Google Scholar] [CrossRef]
  7. Hart, M.A.; Sailor, D.J. Quantifying the influence of land-use and surface characteristics on spatial variability in the urban heat island. Theor. Appl. Climatol. 2009, 95, 397–406. [Google Scholar] [CrossRef]
  8. Taha, H. Urban climates and heat islands: Albedo, evapotranspiration, and anthropogenic heat. Energy Build. 1997, 25, 99–103. [Google Scholar] [CrossRef] [Green Version]
  9. Poë, S.; Stovin, V.; Berretta, C. Parameters influencing the regeneration of a green roof’s retention capacity via evapotranspiration. J. Hydrol. 2015, 523, 356–367. [Google Scholar] [CrossRef]
  10. Bokaie, M.; Zarkesh, M.K.; Arasteh, P.D.; Hosseini, A. Assessment of Urban Heat Island Based on the Relationship between Land Surface Temperature and Land Use/Land Cover in Tehran. Sustain. Cities Soc. 2016, 23, 94–104. [Google Scholar] [CrossRef]
  11. Ichinose, T.; Shimodozono, K.; Hanaki, K. Impact of anthropogenic heat on urban climate in Tokyo. Atmos. Environ. 1999, 33, 3897–3909. [Google Scholar] [CrossRef]
  12. Amiri, R.; Weng, Q.; Alimohammadi, A.; Alavipanah, S.K. Spatial–temporal dynamics of land surface temperature in relation to fractional vegetation cover and land use/cover in the Tabriz urban area, Iran. Remote Sens. Environ. 2009, 113, 2606–2617. [Google Scholar] [CrossRef]
  13. Memon, R.A.; Leung, D.Y.C.; Chunho, L. A review on the generation, determination and mitigation of Urban Heat Island. J. Environ. Sci. 2008, 20, 120–128. [Google Scholar]
  14. Zhou, B.; Rybski, D.; Kropp, J.P. The role of city size and urban form in the surface urban heat island. Sci. Rep. 2017, 7, 4791. [Google Scholar] [CrossRef] [PubMed]
  15. Chapman, S.; Watson, J.E.M.; Salazar, A.; Thatcher, M.; Mcalpine, C.A. The impact of urbanization and climate change on urban temperatures: A systematic review. Landsc. Ecol. 2017, 32, 1921–1935. [Google Scholar] [CrossRef]
  16. Jay, S.G. The Built Environment Induced Urban Heat Island Effect in Rapidly Urbanizing Arid Regions—A Sustainable Urban Engineering Complexity. Environ. Sci. 2004, 1, 321–349. [Google Scholar]
  17. Singh, P.; Kikon, N.; Verma, P. Impact of land use change and urbanization on urban heat island in Lucknow city, Central India. A remote sensing based estimate. Sustain. Cities Soc. 2017, 32, 100–114. [Google Scholar] [CrossRef]
  18. Yang, J.; Sun, J.; Ge, Q.; Li, X. Assessing the impacts of urbanization-associated green space on urban land surface temperature: A case study of Dalian, China. Urban For. Urban Green. 2017, 22, 1–10. [Google Scholar] [CrossRef]
  19. Gill, S.E.; Handley, J.F.; Ennos, A.R.; Pauleit, S. Adapting Cities for Climate Change: The Role of the Green Infrastructure. Built Environ. 2007, 33, 115–133. [Google Scholar] [CrossRef] [Green Version]
  20. Norton, B.A.; Coutts, A.M.; Livesley, S.J.; Harris, R.J.; Hunter, A.M.; Williams, N.S.G. Planning for cooler cities: A framework to prioritise green infrastructure to mitigate high temperatures in urban landscapes. Landsc. Urban Plan. 2015, 134, 127–138. [Google Scholar] [CrossRef]
  21. Oliveira, S.; Andrade, H.; Vaz, T. The cooling effect of green spaces as a contribution to the mitigation of urban heat: A case study in Lisbon. Build. Environ. 2011, 46, 2186–2194. [Google Scholar] [CrossRef]
  22. Nichol, J.E. High-Resolution Surface Temperature Patterns Related to Urban Morphology in a Tropical City: A Satellite-Based Study. J. Appl. Meteorol. 1996, 35, 135–146. [Google Scholar] [CrossRef] [Green Version]
  23. Qiu, G.Y.; Zou, Z.; Li, X.; Li, H.; Guo, Q.; Yan, C.; Tan, S. Experimental studies on the effects of green space and evapotranspiration on urban heat island in a subtropical megacity in China. Habitat Int. 2017, 68, 30–42. [Google Scholar] [CrossRef]
  24. Chen, X.L.; Zhao, H.M.; Li, P.X.; Yin, Z.Y. Remote sensing image-based analysis of the relationship between urban heat island and land use/cover changes. Remote Sens. Environ. 2006, 104, 133–146. [Google Scholar] [CrossRef]
  25. Fu, P.; Weng, Q. A time series analysis of urbanization induced land use and land cover change and its impact on land surface temperature with Landsat imagery. Remote Sens. Environ. 2016, 175, 205–214. [Google Scholar] [CrossRef]
  26. Tran, H.; Uchihama, D.; Ochi, S.; Yasuoka, Y. Assessment with satellite data of the urban heat island effects in Asian mega cities. Int. J. Appl. Earth Obs. Geoinf. 2006, 8, 34–48. [Google Scholar] [CrossRef]
  27. Kusaka, H.; Kimura, F. Coupling a Single-Layer Urban Canopy Model with a Simple Atmospheric Model: Impact on Urban Heat Island Simulation for an Idealized Case. J. Meteorol. Soc. Jpn. 2004, 82, 67–80. [Google Scholar] [CrossRef] [Green Version]
  28. Morini, E.; Touchaei, A.G.; Rossi, F.; Cotana, F.; Akbari, H. Evaluation of albedo enhancement to mitigate impacts of urban heat island in Rome (Italy) using WRF meteorological model. Urban Clim. 2018, 24, 551–566. [Google Scholar] [CrossRef]
  29. Vahmani, P.; Ban-Weiss, G.A. Impact of remotely sensed albedo and vegetation fraction on simulation of urban climate in WRF-urban canopy model: A case study of the urban heat island in Los Angeles. J. Geophys. Res. Atmos. 2016, 121, 1511–1531. [Google Scholar] [CrossRef]
  30. Yu, Z.; Guo, X.; Jørgensen, G.; Vejre, H. How can urban green spaces be planned for climate adaptation in subtropical cities? Ecol. Indic. 2017, 82, 152–162. [Google Scholar] [CrossRef]
  31. Doick, K.J.; Peace, A.; Hutchings, T.R. The role of one large greenspace in mitigating London’s nocturnal urban heat island. Sci. Total Environ. 2014, 493, 662–671. [Google Scholar] [CrossRef] [PubMed]
  32. Bowler, D.E.; Buyung-Ali, L.; Knight, T.M.; Pullin, A.S. Urban greening to cool towns and cities: A systematic review of the empirical evidence. Landsc. Urban Plan. 2010, 97, 147–155. [Google Scholar] [CrossRef]
  33. Senanayake, I.P.; Welivitiya, W.D.D.P.; Nadeeka, P.M. Remote sensing based analysis of urban heat islands with vegetation cover in Colombo city, Sri Lanka using Landsat-7 ETM+ data. Urban Clim. 2013, 5, 19–35. [Google Scholar] [CrossRef]
  34. Weng, Q.; Lu, D.; Schubring, J. Estimation of land surface temperature–vegetation abundance relationship for urban heat island studies. Remote Sens. Environ. 2004, 89, 467–483. [Google Scholar] [CrossRef]
  35. Solecki, W.D.; Rosenzweig, C.; Parshall, L.; Pope, G.; Clark, M.; Cox, J.; Wiencke, M. Mitigation of the heat island effect in urban New Jersey. Glob. Environ. Chang. Part B Environ. Hazards 2005, 6, 39–49. [Google Scholar] [CrossRef]
  36. Santamouris, M. Cooling the cities-A review of reflective and green roof mitigation technologies to fight heat island and improve comfort in urban environments. Sol. Energy 2014, 103, 682–703. [Google Scholar] [CrossRef]
  37. Turner, M.G. Landscape ecology: The effect of pattern on process. Annu. Rev. Ecol. Syst. 1989, 20, 171–197. [Google Scholar] [CrossRef]
  38. Ward, J.; Malard, F.; Tockner, K. Landscape ecology: A framework for integrating pattern and process in river corridors. Landsc. Ecol. 2002, 17, 35–45. [Google Scholar] [CrossRef]
  39. Asgarian, A.; Amiri, B.J.; Sakieh, Y. Assessing the effect of green cover spatial patterns on urban land surface temperature using landscape metrics approach. Urban Ecosyst. 2015, 18, 209–222. [Google Scholar] [CrossRef]
  40. Kong, F.; Yin, H.; James, P.; Hutyra, L.R.; He, H.S. Effects of spatial pattern of greenspace on urban cooling in a large metropolitan area of eastern China. Landsc. Urban Plan. 2014, 128, 35–47. [Google Scholar] [CrossRef]
  41. Li, X.; Zhou, W.; Ouyang, Z.; Xu, W.; Zheng, H. Spatial pattern of greenspace affects land surface temperature: Evidence from the heavily urbanized Beijing metropolitan area, China. Landsc. Ecol. 2012, 27, 887–898. [Google Scholar] [CrossRef]
  42. Maimaitiyiming, M.; Ghulam, A.; Tiyip, T.; Pla, F.; Latorre-Carmona, P.; Halik, U.; Sawut, M.; Caetano, M. Effects of green space spatial pattern on land surface temperature: Implications for sustainable urban planning and climate change adaptation. ISPRS J. Photogramm. Remote Sens. 2014, 89, 59–66. [Google Scholar] [CrossRef] [Green Version]
  43. Masoudi, M.; Tan, P.Y. Multi-year comparison of the effects of spatial pattern of urban green spaces on urban land surface temperature. Landsc. Urban Plan. 2019, 184, 44–58. [Google Scholar] [CrossRef]
  44. Yang, C.; He, X.; Wang, R.; Yan, F.; Yu, L.; Bu, K.; Yang, J.; Chang, L.; Zhang, S. The Effect of Urban Green Spaces on the Urban Thermal Environment and Its Seasonal Variations. Forests 2017, 8, 153. [Google Scholar] [CrossRef] [Green Version]
  45. Chen, A.; Yao, X.A.; Sun, R.; Chen, L. Effect of urban green patterns on surface urban cool islands and its seasonal variations. Urban For. Urban Green. 2014, 13, 646–654. [Google Scholar] [CrossRef]
  46. Cohen, P.; Potchter, O.; Matzarakis, A. Daily and seasonal climatic conditions of green urban open spaces in the Mediterranean climate and their impact on human comfort. Build. Environ. 2012, 51, 285–295. [Google Scholar] [CrossRef]
  47. Li, X.; Zhou, W.; Ouyang, Z. Relationship between land surface temperature and spatial pattern of greenspace: What are the effects of spatial resolution? Landsc. Urban Plan. 2013, 114, 1–8. [Google Scholar] [CrossRef]
  48. Liu, H.; Weng, Q. Scaling Effect on the Relationship between Landscape Pattern and Land Surface Temperature: A Case Study of Indianapolis, United States. Photogramm. Eng. Remote Sens. 2009, 74, 291–304. [Google Scholar] [CrossRef] [Green Version]
  49. Zhou, W.; Qian, Y.; Li, X.; Li, W.; Han, L. Relationships between land cover and the surface urban heat island: Seasonal variability and effects of spatial and thematic resolution of land cover data on predicting land surface temperatures. Landsc. Ecol. 2014, 29, 153–167. [Google Scholar] [CrossRef]
  50. Peng, J.; Jia, J.; Liu, Y.; Li, H.; Wu, J. Seasonal contrast of the dominant factors for spatial distribution of land surface temperature in urban areas. Remote Sens. Environ. 2018, 215, 255–267. [Google Scholar] [CrossRef]
  51. Li, S.; Lu, D.; Huang, J. Impacts of land-cover types on an urban heat island in Hangzhou, China. Int. J. Remote Sens. 2015, 36, 1584–1603. [Google Scholar]
  52. Byrne, J.A.; Lo, A.Y.; Jianjun, Y. Residents’ understanding of the role of green infrastructure for climate change adaptation in Hangzhou, China. Landsc. Urban Plan. 2015, 138, 132–143. [Google Scholar] [CrossRef] [Green Version]
  53. Sheng, L.; Tang, X.; You, H.; Gu, Q.; Hu, H. Comparison of the urban heat island intensity quantified by using air temperature and Landsat land surface temperature in Hangzhou, China. Ecol. Indic. 2017, 72, 738–746. [Google Scholar] [CrossRef]
  54. Qian, Y.; Zhou, W.; Yu, W.; Pickett, S.T. Quantifying spatiotemporal pattern of urban greenspace: New insights from high resolution data. Landsc. Ecol. 2015, 30, 1165–1173. [Google Scholar] [CrossRef]
  55. Sobrino, J.A.; Jimenez-Munoz, J.C.; Paolini, L. Land surface temperature retrieval from LANDSAT TM 5. Remote Sens. Environ. 2004, 90, 434–440. [Google Scholar] [CrossRef]
  56. Jimenez-Munoz, J.C.; Sobrino, J.A.; Skokovic, D.; Mattar, C.; Cristobal, J. Land Surface Temperature Retrieval Methods from Landsat-8 Thermal Infrared Sensor Data. IEEE Geosci. Remote Sens. Lett. 2014, 11, 1840–1843. [Google Scholar] [CrossRef]
  57. Qin, Z.; Karnieli, A.; Berliner, P. A mono-window algorithm for retrieving land surface temperature from Landsat TM data and its application to the Israel-Egypt border region. Int. J. Remote Sens. 2001, 22, 3719–3746. [Google Scholar] [CrossRef]
  58. Wang, F.; Qin, Z.; Song, C.; Tu, L.; Karnieli, A.; Zhao, S. An improved mono-window algorithm for land surface temperature retrieval from Landsat 8 thermal infrared sensor data. Remote Sens. 2015, 7, 4268–4289. [Google Scholar] [CrossRef] [Green Version]
  59. Rozenstein, O.; Qin, Z.; Derimian, Y.; Karnieli, A. Derivation of Land Surface Temperature for Landsat-8 TIRS Using a Split Window Algorithm. Sensors 2014, 14, 5768–5780. [Google Scholar] [CrossRef]
  60. Peng, S.; Piao, S.; Ciais, P.; Friedlingstein, P.; Ottle, C.; Breon, F.M.; Nan, H.; Zhou, L.; Myneni, R.B. Surface Urban Heat Island Across 419 Global Big Cities. Environ. Sci. Technol. 2012, 46, 696–703. [Google Scholar] [CrossRef]
  61. Yu, X.; Guo, X.; Wu, Z. Land surface temperature retrieval from Landsat 8 TIRS—Comparison between radiative transfer equation-based method, split window algorithm and single channel method. Remote Sens. 2014, 6, 9829–9852. [Google Scholar] [CrossRef] [Green Version]
  62. García-Santos, V.; Cuxart, J.; Martínez-Villagrasa, D.; Jiménez, M.; Simó, G. Comparison of Three Methods for Estimating Land Surface Temperature from Landsat 8-TIRS Sensor Data. Remote Sens. 2018, 10, 1450. [Google Scholar] [CrossRef] [Green Version]
  63. Deng, Y.; Wang, S.; Bai, X.; Tian, Y.; Wu, L.; Xiao, J.; Chen, F.; Qian, Q. Relationship among land surface temperature and LUCC, NDVI in typical karst area. Sci. Rep. 2018, 8, 641. [Google Scholar] [CrossRef] [PubMed]
  64. Sobrino, J.A.; Jiménez-Muñoz, J.C.; Sòria, G.; Romaguera, M.; Guanter, L.; Moreno, J.; Plaza, A.; Martínez, P. Land surface emissivity retrieval from different VNIR and TIR sensors. IEEE Trans. Geosci. Remote Sens. 2008, 46, 316–327. [Google Scholar] [CrossRef]
  65. Weng, Q.; Fu, P.; Gao, F. Generating daily land surface temperature at Landsat resolution by fusing Landsat and MODIS data. Remote Sens. Environ. 2014, 145, 55–67. [Google Scholar] [CrossRef]
  66. USGS. Landsat 8 (L8) Data Users Handbook. Available online: https://www.usgs.gov/land-resources/nli/landsat/landsat-8-data-users-handbook (accessed on 27 November 2019).
  67. McGarigal, K.; Cushman, S.A.; Ene, E. FRAGSTATS v4: Spatial Pattern Analysis Program for Categorical and Continuous Maps. Available online: http://www.umass.edu/landeco/research/fragstats/fragstats.html (accessed on 1 May 2019).
  68. Estoque, R.C.; Murayama, Y.; Myint, S.W. Effects of landscape composition and pattern on land surface temperature: An urban heat island study in the megacities of Southeast Asia. Sci. Total Environ. 2016, 577, 349–359. [Google Scholar] [CrossRef]
  69. Li, J.; Song, C.; Cao, L.; Zhu, F.; Meng, X.; Wu, J. Impacts of landscape structure on surface urban heat islands: A case study of Shanghai, China. Remote Sens. Environ. 2011, 115, 3249–3263. [Google Scholar] [CrossRef]
  70. Liu, Q. Variation partitioning by Partial Redundancy Analysis (RDA). Environmetrics 1997, 8, 75–85. [Google Scholar] [CrossRef]
  71. Ter Braak, C.J. The analysis of vegetation-environment relationships by canocical correspondence analysis. Vegetatio 1987, 69, 69–77. [Google Scholar] [CrossRef]
  72. Legendre, P.; Legendre, L. Canonical analysis. Dev. Environ. Model. 2012, 24, 625–710. [Google Scholar]
  73. ter Braak, C.J.; Smilauer, P. Canoco Reference Manual and User’s Guide: Software for Ordination, version 5.0; Microcomputer Power: Ithaca, NY, USA, 2012. [Google Scholar]
  74. Peng, J.; Xie, P.; Liu, Y.; Ma, J. Urban thermal environment dynamics and associated landscape pattern factors: A case study in the Beijing metropolitan region. Remote Sens. Environ. 2016, 173, 145–155. [Google Scholar] [CrossRef]
  75. Yuan, F.; Bauer, M.E. Comparison of impervious surface area and Normalized Difference Vegetation Index as indicators of surface Urban Heat Island effects in Landsat Imagery. Remote Sens. Environ. 2007, 106, 375–386. [Google Scholar] [CrossRef]
  76. Qiu, G.; Li, H.; Zhang, Q.; Chen, W.; Liang, X.; Li, X. Effects of evapotranspiration on mitigation of urban temperature by vegetation and urban agriculture. J. Integr. Agric. 2013, 12, 1307–1315. [Google Scholar] [CrossRef]
  77. Gunawardena, K.R.; Wells, M.J.; Kershaw, T. Utilising green and bluespace to mitigate urban heat island intensity. Sci. Total Environ. 2017, 584, 1040–1055. [Google Scholar] [CrossRef] [PubMed]
  78. Du, H.; Cai, W.; Xu, Y.; Wang, Z.; Wang, Y.; Cai, Y. Quantifying the cool island effects of urban green spaces using remote sensing Data. Urban For. Urban Green. 2017, 27, 24–31. [Google Scholar] [CrossRef]
  79. Hamada, S.; Ohta, T. Seasonal variations in the cooling effect of urban green areas on surrounding urban areas. Urban For. Urban Green. 2010, 9, 15–24. [Google Scholar] [CrossRef]
  80. Chun, B.; Guldmann, J.M. Impact of greening on the urban heat island: Seasonal variations and mitigation strategies. Comput. Environ. Urban Syst. 2018, 71, 167–176. [Google Scholar] [CrossRef]
  81. Greene, C.S.; Millward, A.A. Getting closure: The role of urban forest canopy density in moderating summer surface temperatures in a large city. Urban Ecosyst. 2017, 20, 141–156. [Google Scholar] [CrossRef]
  82. Zhang, Z.; Lv, Y.; Pan, H. Cooling and humidifying effect of plant communities in subtropical urban parks. Urban For. Urban Green. 2013, 12, 323–329. [Google Scholar] [CrossRef]
  83. Wong, N.H.; Yu, C. Study of green areas and urban heat island in a tropical city. Habitat Int. 2005, 29, 547–558. [Google Scholar] [CrossRef]
  84. Zhou, W.; Huang, G.; Cadenasso, M.L. Does spatial configuration matter? Understanding the effects of land cover pattern on land surface temperature in urban landscapes. Landsc. Urban Plan. 2011, 102, 54–63. [Google Scholar] [CrossRef]
  85. Yokohari, M.; Brown, R.D.; Kato, Y.; Moriyama, H. Effects of paddy fields on summertime air and surface temperatures in urban fringe areas of Tokyo, Japan. Landsc. Urban Plan. 1997, 38, 1–11. [Google Scholar] [CrossRef]
  86. Swanwick, C.; Dunnett, N.; Woolley, H. Perspectives on Urban Greenspace in Europe || Nature, Role and Value of Green Space in Towns and Cities: An Overview. Built Environ. 2003, 29, 94–106. [Google Scholar] [CrossRef]
  87. Fisher, B.; Turner, R.K.; Morling, P. Defining and Classifying Ecosystem Services For Decision Making. Ecol. Econ. 2009, 68, 643–653. [Google Scholar] [CrossRef] [Green Version]
  88. Demuzere, M.; Orru, K.; Heidrich, O.; Olazabal, E.; Geneletti, D.; Orru, H.; Bhave, A.G.; Mittal, N.; Feliu, E.; Faehnle, M. Mitigating and adapting to climate change: Multi-functional and multi-scale assessment of green urban infrastructure. J. Environ. Manag. 2014, 146, 107–115. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Location of the study area (RGB combination of Landsat 8 bands 7, 4, and 3).
Figure 1. Location of the study area (RGB combination of Landsat 8 bands 7, 4, and 3).
Sustainability 12 06850 g001
Figure 2. Land use classification schema.
Figure 2. Land use classification schema.
Sustainability 12 06850 g002
Figure 3. Analytical units for different grid scales (taking 1, 2, and 4 km grid scales as an example).
Figure 3. Analytical units for different grid scales (taking 1, 2, and 4 km grid scales as an example).
Sustainability 12 06850 g003
Figure 4. Flowchart of the research methods for investigating the relationship between landscape patterns and land surface temperature (LST).
Figure 4. Flowchart of the research methods for investigating the relationship between landscape patterns and land surface temperature (LST).
Sustainability 12 06850 g004
Figure 5. LSTs of (a) 13 October 2015, (b) 16 December 2015, (c) 15 March 2016, and (d) 27 July 2016. The blank areas in subplot (d) are clouds and cloud shadows.
Figure 5. LSTs of (a) 13 October 2015, (b) 16 December 2015, (c) 15 March 2016, and (d) 27 July 2016. The blank areas in subplot (d) are clouds and cloud shadows.
Sustainability 12 06850 g005
Figure 6. Correlation coefficients between LST and the proportion of different land uses at multiple temporal and spatial scales (500 m: n = 2000, 1000 m: n = 394, 1500 m: n = 156, 2000 m: n = 77, 2500 m: n = 41, 3000 m: n = 29, 3500 m: n = 23, 4000 m: n = 19). * and ** indicate p-value less than 0.05 and 0.01, respectively.
Figure 6. Correlation coefficients between LST and the proportion of different land uses at multiple temporal and spatial scales (500 m: n = 2000, 1000 m: n = 394, 1500 m: n = 156, 2000 m: n = 77, 2500 m: n = 41, 3000 m: n = 29, 3500 m: n = 23, 4000 m: n = 19). * and ** indicate p-value less than 0.05 and 0.01, respectively.
Sustainability 12 06850 g006
Figure 7. Biplots of LST (represented by blue lines) and landscape metrics (red lines) according to RDA. RDA was conducted separately under 1 km (ac), 2 km (df), and 4 km (gi) scales. The length of the arrow represents the standard deviation of variables in the sorting space, and the direction indicates the change of gradient direction. The cosine value of the angle between the arrow and axis represents the correlation. The explanatory landscape pattern variables were selected according to p-value < 0.05.
Figure 7. Biplots of LST (represented by blue lines) and landscape metrics (red lines) according to RDA. RDA was conducted separately under 1 km (ac), 2 km (df), and 4 km (gi) scales. The length of the arrow represents the standard deviation of variables in the sorting space, and the direction indicates the change of gradient direction. The cosine value of the angle between the arrow and axis represents the correlation. The explanatory landscape pattern variables were selected according to p-value < 0.05.
Sustainability 12 06850 g007
Table 1. Landscape pattern metrics used in this study after McGarigal et al. [67].
Table 1. Landscape pattern metrics used in this study after McGarigal et al. [67].
CategoryMetrics (Abbreviation)Unit
DominancePercentage of Landscape (PLAND)%
Largest Patch Index (LPI)%
Shape complexityMean Shape Index (SHMN)unitless
Area-Weighted Mean Shape Index (SHAM)unitless
Mean Fractal Dimension Index (FRMN)unitless
Area-Weighted Mean Fractal Dimension Index (FRAM)unitless
Mean Perimeter–Area Ratio (PAMN)unitless
Area-Weighted Mean Perimeter–Area Ratio (PAAM)unitless
Landscape Shape Index (LSI)unitless
FragmentationTotal Edge (TE)m
Edge Density (ED)m/ha
Number of Patches (NP)unitless
Patch Density (PD)n/km2
Landscape Division Index (DIVISION)%
Aggregation and ConnectednessClumpiness Index (CLUMPY)%
Proportion of Like Adjacencies (PLADJ)%
Interspersion and Juxtaposition Index (IJI)%
Aggregation Index (AI)%
Cohesion Index (COHESION)unitless
Mean Contiguity Index (CONMN)unitless
Area-Weighted Mean Contiguity Index (CONAM)unitless
Table 2. Statistics of LST in different seasons (°C).
Table 2. Statistics of LST in different seasons (°C).
DateMaxMinMeanSD
October 201544.799.2624.932.37
December 201520.81−8.2910.421.32
March 201645.7110.6525.713.33
July 201653.7721.8134.584.23
Table 3. Mean LSTs of different land use types (°C).
Table 3. Mean LSTs of different land use types (°C).
DateUrban ForestUrban WaterCroplandImpervious Surface
October 201523.4 (1.8)22.5 (1.6)24.3 (1.1)26.2 (2.0)
December 20159.5 (1.4)10.8 (0.8)10.2 (0.7)10.8 (1.3)
March 201625.0 (1.8)19.9 (4.3)26.1 (1.3)27.2 (2.1)
July 201632.6 (3.1)28.3 (4.3)33.6 (2.0)37.0 (2.7)
Table 4. Proportions of the explained variations in LST by landscape pattern metrics by using a redundancy analysis (RDA) at different spatial scales.
Table 4. Proportions of the explained variations in LST by landscape pattern metrics by using a redundancy analysis (RDA) at different spatial scales.
CategoryMetrics1 km2 km4 km
UF 1UW 2IS 3UFUWISUFUWIS
DominancePLAND3.7 ** 9.4 ** 16.7 ** 6.1 ** 8.6 ** 13.9 ** 10.4 ** 7.0 * 11.6 **
LPI4.6 ** 9.0 ** 15.6 ** 6.7 ** 7.9 ** 13.5 ** 11.3 ** 7.6 ** 11.5 **
Shape complexitySHMN0.5 ** 3.6 ** 0.3 1.0 6.1 ** 0.7 0.5 13.4 ** 0.7
SHAM0.4 ** 4.1 ** 0.1 0.5 6.2 ** 0.8 0.5 7.3 ** 3.2
FRMN0.3 **5.6 ** 1.4 0.6 6.7 ** 0.9 1.2 13.4 ** 1.2
FRAM0.2 6.7 ** 0.5 0.4 8.1 ** 0.2 0.5 11.1 ** 2.3
PAMN0.2 0.9 ** 0.1 0.3 0.2 0.0 1.0 0.4 1.4
PAAM2.0 ** 3.8 ** 2.7 ** 6.0 ** 5.5 ** 5.0 ** 8.6 ** 4.5 * 5.2 **
LSI1.0 ** 0.8 ** 0.1 1.3 1.3 * 0.3 2.8 2.2 0.4
FragmentationTE0.1 0.1 3.4 ** 0.1 0.1 3.5 ** 0.2 0.1 4.6 **
ED0.1 0.1 3.4 ** 0.1 0.1 3.8 ** 0.4 0.2 5.6 **
NP2.8 ** 0.1 0.1 3.6 ** 0.1 0.4 6.1 ** 0.1 0.1
PD2.8 ** 0.1 0.1 3.8 ** 0.1 0.3 6.1 ** 0.1 0.7
DIVISION4.2 ** 7.8 ** 13.5 ** 6.0 ** 6.8 ** 13.0 ** 9.6 ** 7.2 ** 10.1 **
Aggregation and ConnectednessCLUMPY0.6 ** 3.7 ** 0.5 1.2 4.8 ** 5.4 ** 2.8 4.9 * 8.0 **
PLADJ2.0 ** 3.8 ** 2.7 ** 6.0 ** 5.5 ** 5.0 ** 8.6 ** 4.5 * 5.2 **
IJI0.1 0.8 ** 0.8 0.2 0.6 3.1 ** 1.5 1.1 7.2 **
AI2.6 ** 4.7 ** 2.9 ** 5.6 ** 6.1 ** 4.5 ** 8.4 ** 5.9 * 4.6 **
COHESION2.3 ** 1.5 ** 2.2 7.5 ** 2.1 ** 5.8 ** 9.5 ** 3.4 8.7 **
CONMN0.2 1.0 ** 0.1 0.3 0.3 0.1 0.9 0.3 1.6
CONAM2.2 ** 3.8 ** 2.3 6.3 ** 5.6 ** 4.4 ** 8.9 ** 4.5 5.0 **
Total32.9 71.6 69.7 63.7 82.8 84.6 99.8 99.3 98.8
   — >10,    — 8–10,    — 6–8,    — 4–6,    — 2–4,    — <2. Note: * and ** indicate p-value less than 0.05 and 0.01, respectively. 1 UF: Urban forest, 2 UW: urban water, 3 IS: impervious surface.

Share and Cite

MDPI and ACS Style

Song, Y.; Song, X.; Shao, G. Effects of Green Space Patterns on Urban Thermal Environment at Multiple Spatial–Temporal Scales. Sustainability 2020, 12, 6850. https://0-doi-org.brum.beds.ac.uk/10.3390/su12176850

AMA Style

Song Y, Song X, Shao G. Effects of Green Space Patterns on Urban Thermal Environment at Multiple Spatial–Temporal Scales. Sustainability. 2020; 12(17):6850. https://0-doi-org.brum.beds.ac.uk/10.3390/su12176850

Chicago/Turabian Style

Song, Yu, Xiaodong Song, and Guofan Shao. 2020. "Effects of Green Space Patterns on Urban Thermal Environment at Multiple Spatial–Temporal Scales" Sustainability 12, no. 17: 6850. https://0-doi-org.brum.beds.ac.uk/10.3390/su12176850

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