1. Introduction
Urban areas have become increasingly important centres of economic, social, and cultural activity, attracting people from all walks of life. However, rapid urbanisation has also led to significant land use and cover changes, forming urban heat islands (UHIs). UHIs are characterised by higher temperatures in urban areas compared to their surrounding rural areas due to the absorption of solar radiation by urban surfaces such as buildings, pavements, and roads [
1]. Urban heat islands (UHIs) are a significant phenomenon that has received increasing attention, especially as a footprint of land use change. As cities undergo rapid development and land use change, the impact on UHI intensity becomes a critical concern for sustainable urban planning, denoting the effect (footprint) of those developments/land uses. Several studies have shown that rapid urbanisation-led land use can significantly impact UHI intensity worldwide [
2,
3,
4]. This research also intended to analyse land use changes over time, tracking the expansion of impervious surfaces to map the areas experiencing significant UHIs. However, this research focused more on the intensity of the effect and the analysis of other parameters that might affect the footprint.
UHI intensity is found higher in areas with a higher proportion of built-up land use than areas with more green space [
5,
6,
7,
8,
9,
10]. Thus, a city’s land use change trajectories might also indicate the UHI trajectories [
11,
12,
13]. However, the UHI effect might also vary seasonally for the same land use [
14], signifying other factors that might influence the Land Surface Temperature (LST). The Normalised Difference Vegetation Index (NDVI) represents the extent of green space in land use, and numerous studies have investigated the relationship between NDVI and LST, two important indicators of UHI. Research has examined the relationship between NDVI and LST in the rapidly urbanising area of China and found that vegetation can mitigate the urban heat island effect [
15,
16,
17] via evapotranspiration. Lack of green spaces and choice of construction materials were found to impact LST in semi-arid regions significantly [
18]. Urban vegetation also has a cooling effect on the built environment in subtropical desert cities, showing a strong inversely proportional relation between the oasis effect and NDVI [
19]. Shades from buildings and big trees are found to have a cooling effect on paved surfaces, especially during summer, reducing the LST in those areas. However, trees were found to have a distinct effect when considered an option to mitigate urban heat [
20]. The density of the vegetation is very important when considering reducing UHI effects through urban arboriculture. The research found that dense vegetation and large green spaces, having higher canopy areas, are more effective in reducing LST, than grasslands or scattered vegetation due to the effective reduction of temperature with evapotranspiration [
21,
22].
UHIs can cause a range of negative impacts, including increased energy consumption for cooling buildings, negative impacts on human health [
23,
24], and disruption of urban ecology [
25]. For example, research has shown that UHIs can increase the frequency and severity of heat waves, leading to increased mortality rates among vulnerable populations such as the elderly and children [
26,
27,
28]. UHIs can also disrupt urban ecosystems, with some plant and animal species being particularly sensitive to high temperatures [
29,
30]. In addition, UHIs contribute to climate change through increased energy consumption and greenhouse gas emissions [
31,
32,
33]. Thus, UHIs are a significant issue for urban areas globally, with numerous negative impacts on human health, urban ecology, and the environment. Addressing UHIs requires various strategies that incorporate urban design, planning, and environmental management. Understanding the causes and impacts of UHIs is essential for developing effective policies and strategies that promote sustainable urban development and enhance the resilience of urban areas to climate change.
It is evident globally that vegetation land cover has a strong connection with the UHI, and, therefore, this research focused on a spatial analysis of the relationship between urban land use characteristics (in the form of NDVI) and UHIs for a developing city of South Australia (City of Salisbury). Such a case study area can more clearly portray the developments over the years. The research would be useful in identifying the extent of planning required (regarding vegetation) to mitigate urban heat from the Australian arid climate perspective. To aid the planning strategies, the research further aimed to examine the potential mitigating effects of green spaces and other urban design features on UHI intensity. Research has been conducted in Australia assessing the relations between historical land cover changes and surface temperatures [
34,
35,
36]. However, research hasn’t been conducted that can integrate the planning to combat the local scale climate. Thus, this research provides useful correlations, not only on the effect of land use change but also on the mitigation options and the consequences.
Remote sensing data and geographic information system (GIS) techniques were used to identify the spatial distribution of NDVI in the case study area and analyse its impact on the formation of UHIs. Landsat images have been increasingly used to study urban heat islands (UHIs) due to their ability to capture large-scale land use and land cover changes over time. Chen et al. investigated the spatiotemporal dynamics of UHIs over 17 years (2000-2017) using Landsat images in Wuhan, China [
37]. Similarly, several research used Landsat images for analysing the relationships between the expansion of built-up areas and UHIs in Tehran (Iran), Hong Kong and Istanbul (Turkey) cities, justifying the wider acceptability of this approach, and was adopted in this research as well [
38,
39,
40,
41]. Further to the remote sensing analysis, meteorological data was analysed to understand better the relations, including the factors influencing the relations.
2. Materials and Methods
2.1. Study Area
Salisbury is a city located in the northern suburbs of Adelaide, South Australia (around 25 km north of Adelaide city). It covers an area of approximately 160 square kilometres and has a population of around 147,000 people [
42].
The city of Salisbury is known for its cultural diversity and vibrant community. It is home to a wide range of amenities and attractions, including parks, shopping centres, restaurants, and entertainment venues (
Figure 1 for details). A significant portion of the City of Salisbury covers a conservation area, with more than 50 wetlands, open space and recreation facilities.
Over the years, Salisbury has also become home to several schools and educational institutions, including the University of South Australia’s Mawson Lakes Campus and a campus of TAFE SA. The well-connected transport network to Adelaide via public transport, including bus and train services, has made the city a popular destination for living.
2.2. Data Collection
To investigate urban heat islands effects, three way analysis was conducted by analysing the land use and land cover changes through satellite maps, assessing the significance of meteorological conditions, and investigating the magnitude of land cover effects through aerial photographs.
- a.
Satellite maps
Landsat satellite images were collected for 23 years (1997 — 2020) for land use and LST analysis. Because of the high-altitude acquisition of the images, clouds can interfere and compromise the quality of the maps. Therefore, images were collected for the summer season (December to February) so that cloud free images could be obtained. Also, as the aim is to analyse LST, therefore, it is more logical to collect images for rain free summer days rather than cloudy and rainy winter days in Australia.
Landsat Level-1 data with precision and terrain correction (L1TP) images were collected, which provides radiometric and geodetic accuracy by incorporating ground control points while employing a Digital Elevation Model (DEM) for topographic displacement [
45]. To cover land use from 1997, Landsat 5 (covering the year 1997), Landsat 7 (covering 1999 — 2013) and Landsat 8 (covering 2013—2020) L1TP images were collected from the Landsat image data access portal of U.S. geological survey (USGS) [
46].
Considering the criteria for data collection and from available L1TP images, six years were selected (covering 23 years) for the time series LST analysis, which are 1997 (January), 2004 (February), 2007 (February), 2010 (January), 2016 (December) and 2020 (January).
Band 4 (Red), Band 5 (Near Infrared, NIR) and Band 10 (Thermal Infrared, TIRS) Landsat 8 images are specifically required for NDVI and LST analysis (demonstrated in Figure 5 in detail). The sample of Red, NIR and TIRS images are shown below in
Figure 2,
Figure 3 and
Figure 4. In the case of Landsat 4 to 7 images, Band 3, Band 4 and Band 6 images are required to represent Red, NIR and TIRS.
- b.
Meteorological data
Historical meteorological data for the City of Salisbury was collected from the Australian Bureau of Meteorology Department (BOM) [
47]. These data facilitated a comparison between the LST and meteorological parameters (e.g. rainfall, solar exposure, air temperature etc.), denoting if the LST changes are linked to meteorological parameters.
- c.
Aerial photographs
In order to analyse the magnitude of land use impact on LST with greater precision, high-resolution aerial photographs for the study area were obtained from the government of South Australia’s data portal [
48]. A single aerial photograph was not used for analysing the degree of land use impact on LST, rather the aim was to support the results (land use/cover vs LST) obtained from Landsat data analysis with a high-resolution image for a particular year. These aerial photographs were captured with an aerial thermal sensor during the daytime of 23 March 2018 from ~3000 m altitude, with a 2m resolution. LiDAR (Light Detection and Ranging) remote sensing method was used for the photographs, especially for analysing urban heat islands of the Adelaide region.
2.3. Data Analysis
Land surface temperature (LST) and the normalised difference vegetation index (NDVI) can both be derived from Landsat imagery. LST is estimated using thermal infrared bands (Landsat 8 band 10), which detect the amount of radiation emitted by the Earth’s surface. NDVI, on the other hand, is derived from the ratio of near-infrared and red reflectance bands (usually Landsat 8 bands 4 and 5).
Several steps were typically involved in deriving both LST and NDVI from a Landsat image. The first step was downloading the specific Landsat image bands. The images were then analysed through a series of formulas provided by USGS [
49,
50,
51] and using the ArcMap raster calculation tool.
The estimated LST and NDVI spatial variations were further analysed to assess the trends and correlations between LST and other factors. First, 700 spatial points were randomly selected using the ArcMap from the City of Salisbury (
Figure 6) to do the analysis. Then, the average of all points was compared with the meteorological parameters. NDVI impact on LST was analysed for each of the points.
Figure 5.
Schematic diagram of LST calculation from Landsat satellite image. References for calculation steps:
Steps (a)(b) [
49,
52,
53,
54];
Step (c) [
50,
51,
52,
53,
54];
Steps (d)(e) [
51,
52,
53,
54];
Step (f) [
51,
52,
54]; * Band 10 wavelength range 10.6—11.19 µm [
55], considered 10.8 for the calculation [
52].
Figure 5.
Schematic diagram of LST calculation from Landsat satellite image. References for calculation steps:
Steps (a)(b) [
49,
52,
53,
54];
Step (c) [
50,
51,
52,
53,
54];
Steps (d)(e) [
51,
52,
53,
54];
Step (f) [
51,
52,
54]; * Band 10 wavelength range 10.6—11.19 µm [
55], considered 10.8 for the calculation [
52].
Figure 6.
Randomised sampling points in the study area (City of Salisbury).
Figure 6.
Randomised sampling points in the study area (City of Salisbury).
3. Results
As land use has been changing over the years in the study area, a changed vegetation pattern is also expected in the area. Analysis of the Normalised Difference Vegetation Index (NDVI) from Landsat images can portray those changes. Although Landsat images have a resolution of 60m (every pixel covers 60m x 60 m area), the analysis therefore cannot demonstrate every detail of the changes but can show the bigger picture in general.
Figure 7 shows the analysis of NDVI for the study area since 1997, indicating a major shift from baren areas to bare soil, manifesting the development.
The NDV index in
Figure 7 indicates the expansion of shrubs or grasslands over the study area, which is logically aligned with the development of urban grasslands (lawns, domestic gardens, vacant land, parks, etc.) with new developments. Therefore, the development of urban grasslands can be indicative of the increase in concrete structures as well along those areas. Increased bare soil land area compared to barren areas might support this assumption as well. A significant decrease in water bodies has been observed over the years, relating to a noticeable decrease in vegetation density. While moderate vegetation/canopy development was evident in some places over the years, especially roadside arboriculture, new developments (residential/commercial) also caused reduced vegetation/canopy density at many places. As the Landsat image has 60-m resolution, thus any detected change in vegetation through analysis correlates with a change over a large area, which is significant. This research intended to cross-check whether the NDVI analysis matches with the land use/land cover map, and found a compelling match, as shown in
Figure 8. The comparison was done between 2020 Landsat NDVI analysis map with 2023 satellite land cover image as the 2020 land cover satellite image was not available and assuming the land cover pattern remains the same for 2020 and 2023. The approach was considered legitimate because the aim was to validate whether NDVI analysis done with Landsat images indicated a similar land cover pattern as is found in the satellite images.
Considering the land use and land cover changes over the years, the land surface temperature (LST) was analysed for the study area, depicting a consistent increase in the LST (
Figure 9). The surface temperature was found to increase up to 5
0C in 23 years (1997 to 2003) at any location. However, 40-60 percentile of LST over the city, which indicates the average LST range, spanned between 32—33
0C to 39—39.6
0C (
Figure 10a). If the top LST ranges are considered (more than 90 percentile of LST over the city), the range increased from 36—40
0C to 41—44.7
0C (
Figure 10b). In both the cases, LST increased about 5—7
0C, indicating a consistent increase in LST throughout the study area, which might impact the liveability of the city as well as the ecological vulnerability.
The increase in LST can be correlated with climatic factors; LST was found to be increased with decreasing rainfall and NDVI [
57]. As LST is the radiative temperature of land resulting from solar radiation [
58], a correlation may exist between solar radiation data and increasing LST in the study area. Therefore, in this research, correlation analysis was conducted between LST and rainfall, air temperature and solar radiation data obtained from the local meteorological station for the selected years. As the rainfall duration might affect the vegetation health, therefore, 3-months rainfall data (from the analysed satellite data) was considered to assess the correlation between rainfall and LST. The results are provided in
Figure 11.
Results show that Landsat analysed LST is not statistically significantly correlated with the maximum air temperature (
0C) of the analysis day and has a weak association (correlation coefficient 0.293) (
Figure 11a). This means air temperature (or ambient temperature) had very little and statistically insignificant influence on the LST for the selected year. Rainfall showed a moderate level of negative correlation with the LST (correlation coefficient -0.6), which is expected as the rainfall might increase the vegetation to some extent, yet the association is not statistically significant (
Figure 11b). That is, the level of change in rainfall during summertime couldn’t significantly influence the change in LST over the study area. The correlation between solar exposure and LST was found to be the weakest (coefficient 0.082) and statically insignificant, meaning solar exposure barely had any influence on the LST (
Figure 11c). A more detailed analysis for each year over a period of time might provide a more comprehensive correlation between LST and meteorological parameters. However, the current analysis depicts a lesser influence of the meteorological conditions on the LST. Therefore, the impact of land covers, especially the correlation of LST with NDVI was analysed for each concerned year, as presented in
Figure 12 and
Figure 13.
Figure 12 illustrates the analysis of LST versus NDVI, based on a random sample of 700 points within the study region. The findings indicate that regions without vegetation or with bare ground (NDVI values between -0.1 to 0.2) consistently experience elevated levels of LST and lack discernible differences from one another. The study consistently found that LST is low in areas with water bodies (NDVI ≤ -0.1) and regions with vegetation (NDVI ≥ 0.3). Compared to non-vegetative regions, areas with grasslands or shrubs demonstrated comparatively lower LST levels. The results imply that surface evapotranspiration and evaporation are significant factors that can effectively lower surface temperature.
The research aimed to examine the statistical significance of NDVI’s correlation with LST across various years (only for land areas). As per the findings depicted in
Figure 13, the results demonstrate that the correlation was statistically significant every year (moderate level of correlation in all the cases). This indicates that there is a significant decrease in LST with an increase in greenery in the study region. Further analysis shows that moderate vegetation or areas with widespread canopies can significantly lower the level of LST in urban regions compared to grasslands or shrub areas (
Figure 14).
Canopies reduce the amount of incoming solar radiation that reaches the ground, leading to lesser heat absorption by the ground surface. On the other hand, shrubs have an open canopy that permits more solar radiation to penetrate the ground, resulting in higher temperatures. Furthermore, canopies can enhance the evapotranspiration rate, which refers to the collective loss of water through plant transpiration and soil evaporation. This mechanism contributes to the cooling of the surrounding environment, leading to a decrease in LST. The findings from the analysis corroborated the theoretical justification, thus lending legitimacy to the results obtained in this research.
Figure 15 provides additional spatial analysis that illustrates the impact of different land covers on LST. The results reveal that moderate vegetation or canopy cover can result in a temperature difference of at least 4
0C compared to areas with no vegetation (Box 1 in
Figure 15). Non-vegetative land areas exhibit the highest temperature profile (Box 2). Depending on the type of land use, vegetation can lead to a temperature reduction ranging from 2 to 4
0C (Box 3). Furthermore, incorporating arboriculture along roads can potentially decrease the temperature by at least 1.5
0C compared to paved surfaces. Finally, dense canopies can effectively lower the land surface temperature even further (Box 4).
While the Landsat image analysis demonstrated the significant impact of canopy areas on the LST and showed the changes over the years, the magnitude of impact can be analysed more precisely with the aerial photographs (LiDAR image). This is because of the higher resolution (2 m) of LiDAR compared to Landsat images (30 m).
Figure 16 shows the LST analysis with LiDAR and Landsat images and depicts a similar trend for the impact of vegetation on LST. However, the segmentation of LST in the LiDAR image is more detailed and the impact of each vegetation is accounted. Nevertheless, the comparison shows the acceptability of Landsat image analysis to predict the LST changes, especially for a wide region and for historical changes.
LiDAR image analysis shows that the vegetation can significantly reduce the LST by around 8—13
0C, depending on the surfaces (
Figure 16). This is based on an image captured in March 2018 for the greater Adelaide region. Analysis shows that wherever the canopies are present, the LST remains around 33—34
0C, while wider roads and structures are found to increase LST above 40
0C. Road median strip arboriculture (with canopy >3m) is found to reduce LST in the middle of the wider roads, and so is found on narrow roads with dense arboriculture at the sides. Similarly, large structures are found to increase LST significantly.
Given the inevitability of urban development, the analysis recommends increasing the amount of greenery in order to mitigate the Land Surface Temperature (LST) in urban areas. Urban arboriculture on both sides of roads might increase the density of greenery and thus would reduce the LST further.