Land surface temperature (LST) is essentially considered to be one of the most important indicators used for assessment of the urban thermal environment. It is quite evident that land-use/land-cover (LULC) and landscape patterns have ecological implications at varying spatial scales, which in turn influence the distribution of habitat and material/energy fluxes in the landscape. This article attempts to quantitatively analyse the complex interrelationships between urban LST and LULC landscape patterns with the purpose of elucidating their relation to landscape processes. The study employed an integrated approach involving remote-sensing, geographic information system (GIS), and landscape ecology techniques on bi-temporal Landsat Thematic Mapper images of Southwestern Sydney metropolitan region and the surrounding fringe, taken at approximately the same time of the year in July 1993 and July 2006. First, the LULC categories and LST were extracted from the bi-temporal images. The LST distribution and changes and LST of the LULC categories were then quantitatively analysed using landscape metrics and LST zones. The results show that large differences in temperature existed in even a single LULC category, except for variations between different LULC categories. In each LST zone, the regressive function of LST with fractional vegetation cover (FVC) indicated a significant relationship between LST and FVC. Landscape metrics of LULC categories in each zone in relation to the other zones showed changing patterns between 1993 and 2006. This study also illustrates that a method integrating retrieval of LST and FVC from remote-sensing images combined with landscape metrics provides a novel and feasible way to describe the spatial distribution and temporal variation in urban thermal patterns and associated LULC conditions in a quantitative manner.

The urban thermal environment is impacted by changes in urban landscape patterns resulting from urban expansion and seasonal variation. In order to cope effectively with urban heat island (UHI) effects and improve the urban living environment and microclimate, an analysis of the heating effect of impervious surface areas (ISA) and the cooling effects of vegetation is needed. In this study, Landsat 8 data in four seasons were used to derive the percent ISA and fractional vegetation cover (FVC) by spectral unmixing and to retrieve the land surface temperature (LST) from the radiative transfer equation (RTE). The percent ISA and FVC were divided into four different categories based on ranges 0–25%, 25–50%, 50–75%, and 75–100%. The LST with percent ISA and FVC were used to calculate the surface heating rate (SHR) and surface cooling rate (SCR). Finally, in order to analyze the heating effect of ISA and the cooling effect of vegetation, the variations of LST with SHR and SCR were compared between different percent ISA and FVC categories in the four seasons. The results showed the following: (1) In summer, SHR decreases as percent ISA increases and SCR increases as FVC increases in the study area. (2) Unlike the dependence of LST on percent ISA and FVC, the trends of SHR/SCR as a function of percent ISA/FVC are more complex for different value ranges, especially in spring and autumn. (3) The SHR (heating capacity) decreases with increasing percent ISA in autumn. However, the SCR (cooling capacity) decreases with increasing FVC, except in summer. This study shows that our methodology to analyze the variation and change trends of SHR, SCR, and LST based on different ISA and FVC categories in different seasons can be used to interpret urban ISA and vegetation cover, as well as their heating and cooling effects on the urban thermal environment. This analytical method provides an important insight into analyzing the urban landscape patterns and thermal environment. It is also helpful for urban planning and mitigating UHI.

The urban heat island (UHI) effect is the phenomenon of increased surface temperatures in urban environments compared to their surroundings. It is linked to decreased vegetation cover, high proportions of artificial impervious surfaces, and high proportions of anthropogenic heat discharge. We evaluated the surface heat balance to clarify the contribution of anthropogenic heat discharges into the urban thermal environment. We used a heat balance model and satellite images (Landsat Thematic Mapper (TM) and Enhanced Thematic Mapper Plus (ETM+) images acquired in 1989 and 2001), together with meteorological station data to assess the urban thermal environment in the city of Fuzhou, China. The objective of this study was to estimate the anthropogenic heat discharge in the form of sensible heat flux in complex urban environments. In order to increase the accuracy of the anthropogenic heat flux analysis, the sub-pixel fractional vegetation cover (FVC) was calculated by linear spectral unmixing. The results were then used to estimate latent heat flux in urban areas and to separate anthropogenic heat discharge from heat radiation due to insolation. Spatial and temporal distributions of anthropogenic heat flux were analysed as a function of land-cover type, percentage of impervious surface area, and FVC. The accuracy of heat fluxes was assessed using the ratios of sensible heat flux (H), latent heat flux (L), and ground heat flux (G) to net radiation (R _n), which were compared to the results from other studies. It is apparent that the contribution of anthropogenic heat is smaller in suburban areas and larger in high-density urban areas. However, seasonal disparities of anthropogenic heat discharge are small, and the variance of anthropogenic heat discharge is influenced by urban expansion, land-cover change, and increasing energy consumption. The results suggest that anthropogenic heat release probably plays a significant role in the UHI effect, and must be considered in urban climate change adaptation strategies. Remote sensing can play a role in mapping the spatial and temporal patterns of UHIs and can differentiate the anthropogenic heat from the solar radiative fluxes. The findings presented here have important implications for urban development planning.

Estimating the distribution of impervious surfaces and vegetation is important for analysing urban landscapes and their thermal environment. The application of a crisp classification of land-cover types to analyse urban landscape patterns and land surface temperature (LST) in detail presents a challenge, mainly due to the complex characteristics of urban landscapes. In this article, sub-pixel percentage impervious surface areas (ISAs) and fractional vegetation cover (FVC) were extracted from bitemporal Thematic Mapper/Enhanced Thematic Mapper Plus (TM/ETM+) data by linear spectral mixture analysis (LSMA). Their accuracy was assessed with proportional area estimates of the impervious surface and vegetation extracted from high-resolution data. A range approach was used to classify percentage ISA into different categories by setting thresholds of fractional values and these were compared for their LST patterns. For each ISA category, FVC, LST, and percentage ISA were used to quantify the urban thermal characteristics of different developed areas in the city of Fuzhou, China. Urban LST scenarios in different seasons and ISA categories were simulated to analyse the seasonal variations and the impact of urban landscape pattern changes on the thermal environment. The results show that FVC and LST based on percentage ISA can be used to quantitatively analyse the process of urban expansion and its impacts on the spatial–temporal distribution patterns of the urban thermal environment. This analysis can support urban planning by providing knowledge on the climate adaptation potential of specific urban spatial patterns.

Analyzing the impacts of urban landscape patterns on the thermal environment has become one of the key research areas in addressing urban heat islands (UHIs) and improving the living environment. A case study was carried out in Fuzhou, Fujian Province of China, and bi-temporal Landsat imagery was selected to calculate land surface temperature (LST), percent impervious surface area (ISA), and fractional vegetation cover (FVC). The urban area was further divided into three concentric urban zones, ranging from the city center to the urban periphery, based on urban development densities. The spatial pattern of LST and its variance were analyzed and compared between different zones and different dates. The thermal environment contribution of different zones was also quantified to indicate the change in urban landscape patterns resulting from urban expansion in different zones. Furthermore, Geodetector was used to explore the single factors and interaction factors controlling the spatial patterns of LST in each zone. The results showed that (i) urban expansion primarily increased in Zone 2 and Zone 3, and the areal proportion of high and sub-high LST areas increased from 56.11% and 21.08% to 62.03% and 32.49% in Zone 2 and Zone 3, respectively, from 2004 to 2021; (ii) the heat effect contribution of Zones 2 and 3 reached from 75.16% in 2004 to 89.40% in 2021, indicating that the increase in ISA with >LST_mean was more pronounced in Zone 3 and Zone 2 during the period; (iii) the driving factors of LST spatial distribution were regionally different because of the different landscape patterns, and the explanatory power for the heterogeneity of LST in Zone 1 was weaker than in Zone 2 and Zone 3 in the study area; (iv) the interaction of different factors had a higher explanatory power in the spatial distribution of LST than a single factor in each zone because the distributions of land cover types are heterogeneous in urban areas. The results of this study can be used to improve urban planning for urban ecology and UHI mitigation.

Surface reflectance (SR) estimation is the most essential preprocessing step for multi-sensor remote sensing inversion of geophysical parameters. Therefore, accurate and stable atmospheric correction is particularly important, which is the premise and basis of the quantitative application of remote sensing. It can also be used to directly compare different images and sensors. The Landsat-8 Operational Land Imager (OLI) and Sentinel-2 Multi-Spectral Instrument (MSI) surface reflectance products are publicly available and demonstrate high accuracy. However, there is not enough validation using synchronous spectral measurements over China’s land surface. In this study, we utilized Moderate Resolution Imaging Spectroradiometer (MODIS) atmospheric products reconstructed by Categorical Boosting (CatBoost) and 30 m ASTER Global Digital Elevation Model (ASTER GDEM) data to adjust the relevant parameters to optimize the Second Simulation of Satellite Signal in the Solar Spectrum (6S) model. The accuracy of surface reflectance products obtained from the optimized 6S model was compared with that of the original 6S model and the most commonly used Fast Line-of-sight Atmospheric Analysis of Spectral Hypercubes (FLAASH) model. Surface reflectance products were validated and evaluated with synchronous in situ measurements from 16 sites located in five provinces of China: Fujian, Gansu, Jiangxi, Hunan, and Guangdong. Through the indirect and direct validation across two sensors and three methods, it provides evidence that the synchronous measurements have the higher and more reliable validation accuracy. The results of the validation indicated that, for Landsat-8 OLI and Sentinel-2 MSI SR products, the overall root mean square error (RMSE) calculated results of optimized 6S, original 6S and FLAASH across all spectral bands were 0.0295, 0.0378, 0.0345, and 0.0313, 0.0450, 0.0380, respectively. R² values reached 0.9513, 0.9254, 0.9316 and 0.9377, 0.8822, 0.9122 respectively. Compared with the original 6S model and FLAASH model, the mean percent absolute error (MPAE) of the optimized 6S model was reduced by 32.20% and 15.86% for Landsat-8 OLI, respectively. On the other, for the Sentinel-2 MSI SR product, the MPAE value was reduced by 33.56% and 33.32%. For the two kinds of data, the accuracy of each band was improved to varying extents by the optimized 6S model with the auxiliary data. These findings support the hypothesis that reliable auxiliary data are helpful in reducing the influence of the atmosphere on images and restoring reality as much as is feasible.