Use of Remote Sensing in Comprehending the Influence of Urban Landscape’s Composition and Configuration on Land Surface Temperature at Neighbourhood Scale

The spatial composition and configuration of land use land cover (LULC) in the urban landscape impact the land surface temperature (LST). In this study, we assessed such impacts at the neighbourhood level of the City of Edmonton. In doing so, we employed Landsat-8 Operational Land Imager (OLI) and Thermal Infrared Sensors (TIRS) satellite images to derive LULC and LST maps, respectively. We used three classification methods, such as ISODATA, random forest, and indices-based, for mapping LULC classes including built-up, water, and green. We obtained the highest overall accuracy of 98.53 and 97.90% with a kappa value of 0.96 and 0.92 in the indices-based method for the 2018 and 2015 LULC maps, respectively. Besides, we estimated the LST map from the brightness temperature using a single-channel algorithm. Our analysis showed that the highest contributors to LST were the industrial (303.51 K in 2018 and 295.99 K in 2015) and residential (303.47 K in 2018 and 296.56 K in 2015) neighbourhoods, and the lowest contributor was the riverine/creek (298.77 K in 2018 and 292.89 K in 2015) during the 2018 late summer and 2015 early spring seasons. We also found that the residential neighbourhoods exhibited higher LST in comparison with the industrial with the same LULC composition. The result was also supported by our surface albedo analysis, where industrial and residential neighbourhoods were giving higher and lower albedo values, respectively. This indicated that the rooftop materials played further role in impacting the LST. In addition, our spatial autocorrelation (local Moran’s I) and proximity (near distance) analyses revealed that the structural configurations would additionally play an important role in contributing to the LST in the neighbourhoods. For example, the cluster pattern with a small gap of minimum 2.4 m between structures in the residential neighbourhoods were showing higher LST in compared with the sparse pattern, with large gaps between structures in the industrial areas. The wide passages for wind flow through the large gaps would be responsible for cooling the LST in the industrial neighbourhoods. The outcomes of this study would help planners in planning and designing urban neighbourhoods, and policymakers and stakeholders in developing strategies to balance surface energy and mitigate local warming.


Introduction
Urbanization has become one of the most critical issues across the earth. In fact, about half of the world's population and 75% population of developed nations are living in urban cities [1]. The, world's urban population may increase from 4.22 billion in 2018 to 6.68 billion in 2050 due to better socioeconomic opportunities [2]. In order to accommodate the need of the growing population, it is common to convert forested land, agricultural land, wetland, water bodies, and open spaces into built-up/urban areas [3]. One of the major consequences of urbanization is the increase in local temperature as a result of a higher proportion of built-up/impervious surfaces in comparison to rural areas [4]. This is because built-up areas show higher sensible heat due to a lower amount of evaporative (i.e., water bodies) and transpirative (vegetation) surfaces [5]. Consequently, it is imperative to study the impact of urban land use and land cover (LULC) in terms of their composition (i.e., percentage coverage) and configuration (i.e., spatial arrangement) on land surface temperature (LST).
In order to examine the above-mentioned issue, it is worthwhile to explore the methods of generating LULC maps. One of the most commonly used methods is the traditional field survey; however, it is a very costly, time-consuming, and tedious process [6]. Alternatively, satellite-based remote sensing is becoming popular as these platforms provide continuous near real-time coverage of the earth surface that is readily available at a lower cost, with various spatial resolutions ranging from high, medium, to coarser ones. In general, high/finer spatial resolution images are useful to generate detailed LULC mapping of urban surfaces [7], but they are usually costly. Thus, the use of medium to coarse spatial resolution images is a viable option as they are available at little or no cost over a longer time-period, e.g., since the early 1980s. Consequently, it is common that researchers have considered using Landsat (medium resolution) and MODIS (coarse resolution) images more often for this purpose [8][9][10][11]. For example, Afrin et al. [8] applied ISODATA (iterative self-organizing data analysis) clustering (a form of unsupervised classification) to both Landsat and MODIS images over the Athabasca River watershed in Alberta, Canada. Akbar et al. [9] used normalized difference vegetation index (NDVI) thresholding concept (a form of indices-based classification) using Landsat images in Lahore, Pakistan. Abdullah et al. [11] used Landsat data with Random Forest (RF: a form of supervised classification) and gradient-boosting classification schemes over coastal region of Bangladesh. In this study, we opted to explore the above-mentioned classification types over an urban landscape.
Though it is known that both composition and configuration of LULC influence LST over space and time; however, researchers have made significant efforts to quantify the impact of composition. For example, Zhao et al. [12] analyzed the profile and concentric zonal relationship using Landsat images in Shenyang, China. For profile zonation, they created a circle with a radius of 14 km from the city center and divided the circle into eight equal zones for analyzing association between LULC composition and LST. They drew concentric circles at every 0.5 km intervals up to 9.5 km from the city center that resulted into 19 zones for further assessment. In another study, Ullah et al. [10] quantified the changes in LULC and its impact on LST over the period 1990-2017 using Landsat images for the Mansehra and Battagram districts in Northern Pakistan, in the lower Himalayan region. They conducted the analysis at regional scale upon considering the two districts as one entity. Furthermore, Huo et al. [13] explored LULC-LST relations between 1990-2015 using Landsat data in the Kuala Lumpur, Malaysia, where they evaluated such relations at city-scale. However, when researchers attempted to analyze LULC-LST relations at micro-level (e.g., dividing the city landscape into multiple segments), they observed that LST-values would vary for one segment to another even with the same LULC composition (e.g., 70% of built-up and 30% of green space). This might be attributed to the landscape configuration, in other words, the spatial arrangement of LULC features would be one of the major contributing factors [14][15][16][17][18].
In the context of LULC configuration and its impact on LST, researchers generally used two methods, i.e., land fragmentation metrics [3,7,17] and spatial autocorrelation [14,15,19] (see Table 1 for some example case studies). In case of fragmentation metrics when applied at landscape level, it would consider all the patch types concurrently. Thus, their application would often fail to quantify the impact of a specific LULC type on LST [14]. In order to address this issue, one of the solutions would be the use of spatial autocorrelation, e.g., local Moran's I [15,20]. Note that densely placed structures would show increased LST in comparison to the sparsely placed ones [14,15,20]. In the scope of this paper, our overall goal was to quantify the impact of both the composition and configuration of land surface on LST over the Canadian City of Edmonton. For analysing the relationship, we considered neighbourhood scale. Note that neighbourhood is the fundamental unit in the Canadian context that facilitates a mixture of land uses, including residential, commercial/industrial, and amenities/facilities focused around a neighbourhood node [21,22]. We therefore formulated a set of three specific objectives. Firstly, we delineated an urban LULC map at 30 m using Landsat 8 data with three major classes, viz. built-up, water body and green space with the aid of ISODATA Remote Sens. 2020, 12, 2508 4 of 17 clustering, RF, and indices-based classifiers. Then, we determined the best LULC map upon comparing them against visually identified LULC classes. Secondly, we derived an LST map at 30 m from Landsat 8 data and compared them against LULC classes at the neighbourhood level in quantifying their inter-relationships. Finally, we generated spatial dynamics of the rooflines of the structures at 2 m from a GIS database. We then calculated local Moran's I and minimum distance among structures at neighbourhood level and evaluated their relations with the Landsat 8-drived LST under similar LULC compositions (e.g., 0-10% water body and green space with 90-100% built-up). Note that the concept of applying the minimum distance and analysing at neighbourhood level is unique to the best of our knowledge.

General Description of the Study Area
We considered the corporate boundary of the City of Edmonton as our study area (see Figure 1), which is also the Capital of the Canadian Province of Alberta. It is the second largest city of the Province, situated between latitude 53 • 20 to 53 • 43 N and longitude 113 • 15 to 113 • 45 W with an area of 783 km 2 , i.e., about an average of 25 and 32 km in east-west and north-south directions, respectively. The city comprises 400 neighbourhoods, where the North Saskatchewan River passes through the heart of the city. The city's population has increased significantly over the 1960-2019 period, i.e., 269,314 to 972,223 as per the Edmonton Municipal Census [23]. The number of dwellings has increased about 19.8% over the 2008-2019 period [23]. The city landscape has high-rise buildings in the central zone, and high-density developments across the majority of the neighbourhoods. It also has riverine zones, parks and patches of green space, and agricultural lands in its fringe. In terms of climate, the city exhibits a daily average temperature of −10.4 • C in January (i.e., winter peak) and 17.7 • C in July (i.e., summer peak) over the period 1981-2010, and the annual precipitation of 455.7 mm, where 347.8 mm comes as rainfall and the rest as snowfall [24].

Data Requirements and Processing
We used three sets of data in this study, i.e., (i) Landsat-8 operational land imager (OLI) and thermal infrared sensors (TIRS) Level 2 products that were both atmospherically and geometrically corrected and provided by United States Geological Survey [25], (ii) Google Earth Pro images, and (iii) neighbourhood boundaries and rooflines shapefiles of the City of Edmonton. In particular to the Landsat-8 images, we found cloud-free images acquired on 07 September 2018 (late summer) and 10 May 2015 (early spring) with Universal Transverse Mercator projection system at Zone 12 Northern hemisphere (UTM Zone 12N) with World Geodetic System 1984 (WGS84) datum. Table 2 summarizes the details of Landsat-8 data and their usages.  In addition, we used Google Earth Pro to visually identify all LULC types for both training and validation purposes for executing the employed LULC classification methods. We obtained two shapefiles from the official geospatial data repository of the City of Edmonton [26], i.e., (i) 'neighbourhood boundaries' representing 400 neighbourhoods as of May 2020, for assessing the LST in comparison with the LULC classes at the neighbourhood level and (ii) 'rooflines' describing the In addition, we used Google Earth Pro to visually identify all LULC types for both training and validation purposes for executing the employed LULC classification methods. We obtained two shapefiles from the official geospatial data repository of the City of Edmonton [26], i.e., (i) 'neighbourhood boundaries' representing 400 neighbourhoods as of May 2020, for assessing the Remote Sens. 2020, 12, 2508 6 of 17 LST in comparison with the LULC classes at the neighbourhood level and (ii) 'rooflines' describing the roof outlines of the structures as of May 2017, for understanding their spatial arrangement and impact on LST. Both shapefiles were originally found in the geographic coordinate system (GCS) with North America Datum 1983 (NAD83), which we re-projected to the UTM Zone12N WGS84 projection system to conform with the Landsat data used in this study. Figure 2 illustrates the proposed methods in the form of a schematic diagram. Here, there were three major components: (i) generation of LULC maps from Landsat-8 OLI image using ISODATA clustering, RF, and indices-based classifiers and their comparisons; (ii) derivation of LST from Landsat-8 TIRS image, and establishment of relations between LST and LULC composition at the neighbourhood scale; (iii) analysis of the influence of the structural configurations on the LST at neighbourhood scale. Brief descriptions about these components are as follows. Figure 2 illustrates the proposed methods in the form of a schematic diagram. Here, there were three major components: (i) generation of LULC maps from Landsat-8 OLI image using ISODATA clustering, RF, and indices-based classifiers and their comparisons; (ii) derivation of LST from Landsat-8 TIRS image, and establishment of relations between LST and LULC composition at the neighbourhood scale; (iii) analysis of the influence of the structural configurations on the LST at neighbourhood scale. Brief descriptions about these components are as follows.

Generation of LULC Classifications and Accuracy Assessments
We implemented three different methods of ISODATA, RF, and indices-based classification on the Landsat-8 OLI image 2018. We intended to derive three major LULC classes including built-up, water, and green (see Table 3 for details). Table 3. Description of the LULC classes.

LULC Class Description
Built-up Developments in urban including residential areas, impervious surfaces (bitumen and concrete), and industrial areas. Water Open water bodies including the river, creeks, lakes, and pools (both natural and artificial).  For the ISODATA clustering method, we initially generated 50 classes with a convergence threshold of 0.995 and keeping an infinite number of iterations [27,28]. We further evaluated their

Generation of LULC Classifications and Accuracy Assessments
We implemented three different methods of ISODATA, RF, and indices-based classification on the Landsat-8 OLI image 2018. We intended to derive three major LULC classes including built-up, water, and green (see Table 3 for details). Table 3. Description of the LULC classes.

LULC Class Description
Built-up Developments in urban including residential areas, impervious surfaces (bitumen and concrete), and industrial areas. Water Open water bodies including the river, creeks, lakes, and pools (both natural and artificial).
Green Vegetation (including grasses, shrubs, and trees), agricultural lands (with and without crops), bare surfaces, and open spaces. For the ISODATA clustering method, we initially generated 50 classes with a convergence threshold of 0.995 and keeping an infinite number of iterations [27,28]. We further evaluated their class-specific spectral signatures to group into our intended three LULC classes (see Table 3). To aid this process, we consulted spectral bands of the Landsat image with various combinations, and Google Earth Pro images.
For the RF method, we generated training polygons on the Landsat image by observing feature clarity in the Google Earth Pro images. We then calculated and integrated the spectral attributes (i.e., the maximum, minimum, and mean surface reflectance values of the Landsat image) in each training polygon. Next, we derived a machine intellect data using the attributed training polygons, which was further used to generate a thematic image with our intended three LULC classes (see Table 3).
For the indices-based method, we calculated three indices, i.e., NDVI: a measure of greenness [29], modified normalized difference water index (MNDWI): a measure of surface water body [30], and normalized difference built-up index (NDBI): a measure of built-up in urban [31], using the following expressions where, ρ NIR , ρ Red , ρ Green , and ρ SWIR2 are the surface reflectance values of NIR, red, green and SWIR2 spectral bands, respectively. Next, we stacked these three index-based images, and identified samples of the intended LULC classes with the aid of Google Earth Pro images. This process led us to generate class-specific signatures, which were then used to derive the LULC map.
In assessing the accuracy of the LULC maps, we generated a set of reference polygons on the Landsat-8 OLI image with the help of Google Earth Pro images, which made up about 2% of the study area. Then, we calculated a confusion matrix by comparing the classes in each map with the corresponding reference polygons. We finally calculated the values of 'overall accuracy' and 'kappa statistics' for the three LULC maps to decide the best one. This best method was further applied on the Landsat-8 OLI 2015 image to derive an LULC map of 2015 and performed the accuracy assessment as well.

Derivation of LST and Its Relations with LULC Composition
We generated an LST image from the Landsat-8 TIRS image using single-channel algorithms [10,32,33] as follows where T B is the brightness temperature; λ is the center wavelength (10. In quantifying the relations between LST and LULC at the neighbourhood scale, we grouped the 400 neighbourhoods of the City of Edmonton into five major land-use categories, i.e., industrial, riverine/creek, freeway, agricultural, and residential. During the process, first we considered if land-use type was in the name attribute. Next, we visually identified from the spectral response, association, and pattern of the land use in the Landsat image. Table 4 shows a summary of the criteria we used for grouping the neighbourhoods in each category. Table 4. Criteria used for grouping the neighbourhoods into five different categories.

Industrial
Having 'industrial' word in the name attribute, large structural footprints, and occupied more than 50% of a neighbourhood.
Riverine/creek Having word 'river' or 'creek' in the attribute, recreational park, golf course, and located along river and creek with vegetation canopy. Freeway Roads with the words 'Anthony Henday' in the attribute. Agricultural More than 50% agricultural land, bare land, and open space. Residential Remaining neighbourhoods.
As we observed a visually clear distinction between the residential and industrial categories in terms of surface reflections, thus we considered to calculate 'surface albedo' using Equation (7) [34] to understand their impact on LST. Note that surface albedo is the fraction of incident sunlight reflected by the surface of the Earth, which is primarily governed by the properties of the land surfaces, i.e., LULC types.
For understanding the LULC composition, we calculated the percentage of the three LULC classes in each neighbourhood. Then, we subcategorized the industrial and residential neighbourhoods considering the total percentage of combined water and green classes in each neighbourhood, as shown in Table 5. Note that we did not subcategorized the remaining three neighbourhood categories (i.e., riverine/creek, freeway, and agricultural), because the percentages of combined water and green were over 50%. Next, we calculated both mean values of LST and albedo for each subcategory associated with both industrial and residential neighbourhood categories to quantify the LULC-LST relationships. Then, we performed two sets of linear regressions, and determined the coefficient of determination (r 2 )-values between subcategories and: (i) LST, and (ii) albedo, for each category of interest. Table 5. Criteria used for the subcategorization of the industrial and residential neighbourhoods.

Neighbourhood Subcategory Water and Green Combined (%) Industrial
Residential

Analysis of the Influence of Structural Configurations on the LST at Neighbourhood Scale
In this analysis, we considered to apply spatial autocorrelation and proximity analysis. In order to execute the spatial autocorrelation method, we converted the rooflines vector data to raster with a spatial resolution of 2 m to capture all gaps among structures. In the study area, the minimum distance Remote Sens. 2020, 12, 2508 9 of 17 between two structures would be 2.4 m in the study area according to the regulations [35]. Next, we calculated local Moran's I index, I i (d)-values of raster data with 3 × 3 window, as follows [36] where Z i is the attribute value of interest at location i, Z m is the mean of the attribute Z values, w ij (d) is a spatial weight matrix, and d (= 3 × 3) represents the size of the calculation window. The I i (d)-values range from −1 to +1 for highly dispersed to clustered patterns, respectively. For the proximity analysis, we applied 'near distance' calculations on the rooflines vector data, on the basis of finding the shortest separation between two structures. For analysing the influence of the structural configurations on the LST, we calculated zonal mean values of the both local Moran's I and near distance for each subcategory of both industrial and residential neighbourhoods. Next, we used these mean values to carry out two sets of linear regression analyses between LST and (i) local Moran's I index, (ii) near distance. In this case, we considered a set of three criteria: (i) use of residential and industrial neighbourhoods together, with the same percentages of combined water and green (e.g., I10 and R10 together); (ii) no use of I50+ and R50+ subcategories, because of their having more than 50% combined water and green; (iii) calculations of mean values using bins with a size of 0.03 and 5 m for local Moran's I index and near distance, respectively.

Evaluation of the LULC Maps and Dynamics of LST
We found the overall accuracies of 98.08%, 97.81%, and 98.53%; and kappa statistics of 0.95, 0.94, and 0.96 for the methods of ISODATA, RF, and indices-based respectively for 2018 (see Table 6). The best LULC map (i.e., indices-based) represented an area of 49.41%, 2.28%, and 48.32% for the built-up, water, and green classes, respectively (see Figure 3a). Furthermore, we estimated that the minimum, mean, and maximum LST were 286.08, 301.45, and 312.8 K, respectively. We observed the lower LST in and around the water and green LULC classes, and the higher LST in the built-up areas (see Figure 3b). For 2015, the overall accuracy and kappa statistics of the LULC map were found to be 97.90% and 0.92, respectively (see Table 6); and the estimated minimum, mean and maximum LST were 285.08, 296.08, and 302.88 K, respectively.

Relationships between LST and LULC Composition
The industrial, riverine/creek, freeway, agricultural, and residential categories consisted of 71, 27, 14, 25, and 263 neighbourhoods, respectively (see the Total column in Table 7; Figure 4a). Further, the number of subcategorized neighbourhoods of the industrial and residential are shown in the 'subtotal' columns of Table 7, and their spatial locations are shown in Figure 4b.     Our analysis revealed that the highest and lowest mean LST values were in the industrial and riverine/creek neighbourhood categories, respectively, where the values were very close to each other for the industrial and residential neighbourhoods (see 2nd and 3rd columns in Table 8). However, when the components of built-up, and combined green and water in the LULC composition were considered separately, we observed the highest and lowest mean LST values in the residential and riverine/creek categories, respectively (see 4th to 7th columns in Table 8). The regression analysis between the mean LST values and neighbourhood subcategories of: (a) I10 to I50+, and (b) R10 to R50+, showed strong relationships, i.e., r 2 -values of ~0.90 and 0.35 for industrial, and ~0.95 and 0.52 for the residential subcategories, respectively (see Figure 5). Note that we used I10/R10 (=1) to I50+/R50+(=6) for the neighbourhood subcategory in the regression analysis. In all cases, the mean LST values were gradually decreasing with the increasing percentage of combined green and water component in the LULC composition. However, residential subcategories were showing higher mean LST values in comparison with the same composition of the industrial Our analysis revealed that the highest and lowest mean LST values were in the industrial and riverine/creek neighbourhood categories, respectively, where the values were very close to each other for the industrial and residential neighbourhoods (see 2nd and 3rd columns in Table 8). However, when the components of built-up, and combined green and water in the LULC composition were considered separately, we observed the highest and lowest mean LST values in the residential and riverine/creek categories, respectively (see 4th to 7th columns in Table 8). The regression analysis between the mean LST values and neighbourhood subcategories of: (a) I10 to I50+, and (b) R10 to R50+, showed strong relationships, i.e., r 2 -values of~0.90 and 0.35 for industrial, and~0.95 and 0.52 for the residential subcategories, respectively (see Figure 5). Note that we used I10/R10 (=1) to I50+/R50+(=6) for the neighbourhood subcategory in the regression analysis. In all cases, the mean LST values were gradually decreasing with the increasing percentage of combined green and water component in the LULC composition. However, residential subcategories were showing higher mean LST values in comparison with the same composition of the industrial subcategories, which might be due to the differences in surface albedo characteristics as depicted in Figure 6a. In fact, we observed the higher albedo values (i.e., lower absorption of incident sunlight) in the industrial subcategories and considered this as one of the responsible factors for the relatively lower LST values when compared with the residential subcategories of neighbourhoods of both 2018 and 2015 (see Figure 6b).
Remote Sens. 2020, 12, x FOR PEER REVIEW  12 of 17 subcategories, which might be due to the differences in surface albedo characteristics as depicted in Figure 6a. In fact, we observed the higher albedo values (i.e., lower absorption of incident sunlight) in the industrial subcategories and considered this as one of the responsible factors for the relatively lower LST values when compared with the residential subcategories of neighbourhoods of both 2018 and 2015 (see Figure 6b).

Influence of Structural Configurations on LST
The  Table 9 for details). The negative slopes of the regressions (see 3 rd and 6 th columns in Table 9) indicated that the higher values of the both local Moran's I and near distance would contribute to the lessening of LST values. Table 9. Linear regression analysis between the mean LST values and mean values of (a) local Moran's I, and (b) near distance for the neighbourhood subcategories.   Table 9 for details). The negative slopes of the regressions (see 3 rd and 6 th columns in Table 9) indicated that the higher values of the both local Moran's I and near distance would contribute to the lessening of LST values. Table 9. Linear regression analysis between the mean LST values and mean values of (a) local Moran's I, and (b) near distance for the neighbourhood subcategories.

Year
Neighbourhood Subcategory

Discussion
Despite the similar overall accuracies (i.e., above 97%) and kappa statistics (i.e., above 94%) obtained from all the three classification methods in generating LULC maps, we preferred the indices-based LULC map for further analysis. This was due to the fact that this method would be simple and require less computational time in comparison to the RF and ISODATA clustering methods. In fact, our LULC map showed the green areas were in the neighbourhoods at the fringes of the city [37], which would be subject to continual development. Green areas were also found in neighbourhoods along the North Saskatchewan River, and other creeks, which were specifically developed for recreational parks and golf courses with no possibility of further developments [38]. Moreover, the central part of the city experienced the residential neighbourhoods' establishment at least 50 years ago [39], therefore they would have more grown tree/vegetation around the structures. Additionally, the industrial neighbourhoods at the fringes of the city had more underdeveloped land consisting of more green LULC type in comparison to the neighbourhoods in the city core [40]. In general, the LULC composition in the study area showed the cluster-development of built-up areas in the city core, which was also reported in other studies [41][42][43].
Our results for the five different categories of neighbourhoods in the study area identified that both industrial and residential neighbourhoods exhibited high mean LST values compared to the other three neighbourhoods (i.e., riverine/creek, agricultural, and freeway). The industrial and residential neighbourhoods comprised a high density of built-up with impervious (artificial) surface probably causing higher LST values. Similar findings were also reported in Phoenix [7,14] and Las Vegas [44] cities of the United States. We also found lower mean LST values in the riverine/creek and agricultural neighbourhoods, where riverine/creek was the lowest. The neighbourhoods in the riverine/creek category consisted of over 80% green and water LULC classes, which probably had a

Discussion
Despite the similar overall accuracies (i.e., above 97%) and kappa statistics (i.e., above 94%) obtained from all the three classification methods in generating LULC maps, we preferred the indices-based LULC map for further analysis. This was due to the fact that this method would be simple and require less computational time in comparison to the RF and ISODATA clustering methods. In fact, our LULC map showed the green areas were in the neighbourhoods at the fringes of the city [37], which would be subject to continual development. Green areas were also found in neighbourhoods along the North Saskatchewan River, and other creeks, which were specifically developed for recreational parks and golf courses with no possibility of further developments [38]. Moreover, the central part of the city experienced the residential neighbourhoods' establishment at least 50 years ago [39], therefore they would have more grown tree/vegetation around the structures. Additionally, the industrial neighbourhoods at the fringes of the city had more underdeveloped land consisting of more green LULC type in comparison to the neighbourhoods in the city core [40]. In general, the LULC composition in the study area showed the cluster-development of built-up areas in the city core, which was also reported in other studies [41][42][43].
Our results for the five different categories of neighbourhoods in the study area identified that both industrial and residential neighbourhoods exhibited high mean LST values compared to the other three neighbourhoods (i.e., riverine/creek, agricultural, and freeway). The industrial and residential neighbourhoods comprised a high density of built-up with impervious (artificial) surface probably causing higher LST values. Similar findings were also reported in Phoenix [7,14] and Las Vegas [44] cities of the United States. We also found lower mean LST values in the riverine/creek and agricultural neighbourhoods, where riverine/creek was the lowest. The neighbourhoods in the riverine/creek category consisted of over 80% green and water LULC classes, which probably had a strong cooling effect on the urban landscape [44], and thereby helped to lower the mean LST value. Additionally, agricultural neighbourhoods comprised agricultural lands with or without crops including some bare and exposed soils, which exhibited relatively higher LST in comparison to the riverine/creek ones. In fact, Ullah et al. [10] found that bare soil had higher LST values compared to water and vegetation, which might be similar to our riverine/creek category in terms of composition. Further, the freeway neighbourhood category comprised wider roads, grass, and bare lands; thus it showed a slightly higher mean LST value in comparison with riverine/creek and agricultural neighbourhoods that agreed quite well with the findings of Myint et al. [14].
Our analysis of the subcategories of both industrial and residential neighbourhoods revealed that the residential ones exhibited higher mean LST values with the same LULC composition as the industrial ones. In fact, our observations contradicted with other studies (e.g., [45,46]), because residential neighbourhoods ought to have lower LST compared to the industrial neighbourhoods. However, our albedo analysis, i.e., higher albedo values over the industrial subcategories to the residential ones, might be one of the potential explanations. Note that high albedo or the low absorption of incident sunlight was probably attributed to having rooftops with coated steels or metal sheets over the industrial neighbourhoods [47][48][49]. On the other hand, the rooftops/shingles in the residential neighbourhoods would be mainly made of asphalt with dark colour, so they might absorb more sunlight [49][50][51]. Consequently, they would exhibit higher LST values. In addition, the small distance of 2.4 m between structures in residential neighbourhoods probably facilitated densified clustered configuration. On the contrary, the dispersed configuration of industrial neighbourhoods with large open spaces and gaps between structures would create a large passage for wind flow that might help reducing the LST. This was also clear for local Moran's I (as implemented elsewhere [7,14]) and proximity (near distance) analyses. Both analyses showed that the densified clustered configuration of the residential neighborhoods caused the high LST, while the dispersed configuration of industrial caused lower LST.

Conclusions
In this study, we demonstrated the influence of LULC composition and configuration on LST at the 400 neighbourhoods of the City of Edmonton. To aid this process, we generated LULC maps using Landsat-8 OLI image, and found that the indices-based method of LULC classification performed the best. Besides, we estimated an LST map from the Landsat-8 TIRS image with single-channel algorithm and observed that the built-up and water LULC classes showed the highest and lowest LST values, respectively. We further analysed the impacts of LULC composition on the five neighbourhood categories. We found that the industrial and riverine/creek neighbourhoods had the highest and lowest mean LST values, respectively, although industrial and residential showed very close values with a slight difference, i.e., 0.04 K (2018) and 0.57 K (2015). On the other hand, while we compared this with the subcategories of having same LULC composition in the residential and industrial neighbourhoods, we observed that the residential subcategories showed higher mean LST values. Such a high LST in the residential would be because of using the dark-coloured roof-top materials in the study area, as confirmed by our surface albedo analysis. Therefore, we suggest adopting light-coloured roof-top materials with high albedo properties in the residential structures. Further, local Moran's I and proximity (near distance) analyses identified that structural configurations of cluster pattern would play a role in enhancing the LST.
Although our study explored the efficacy of Landsat data at the neighbourhood level, aerial photographs, or very high spatial resolution satellite image (e.g., less than 2 m) would explore further insight. Such images would facilitate generating a very high spatial resolution LULC composition to capture further details even between structures. This study would help policymakers, urban planners, and related stakeholders to have a clear understanding of the impacts of various landscape composition and structural configurations on surface energy balance in maintaining sustainable neighbourhoods in an urban setup. Despite the promising outcomes of this study, we strongly recommend evaluating the proposed method, and it should be investigated thoroughly before being applied in other cities.