
==== Front
PLoS One
PLoS One
plos
PLOS ONE
1932-6203
Public Library of Science San Francisco, CA USA

10.1371/journal.pone.0307711
PONE-D-24-01189
Research Article
Ecology and Environmental Sciences
Terrestrial Environments
Urban Environments
Engineering and Technology
Remote Sensing
Earth Sciences
Geography
Earth Sciences
Geography
Human Geography
Land Use
Social Sciences
Human Geography
Land Use
Computer and Information Sciences
Geoinformatics
Spatial Autocorrelation
Earth Sciences
Geography
Geoinformatics
Spatial Autocorrelation
Earth Sciences
Geography
Human Geography
Urban Geography
Urbanization
Social Sciences
Human Geography
Urban Geography
Urbanization
Earth Sciences
Geography
Human Geography
Urban Geography
Urban Areas
Social Sciences
Human Geography
Urban Geography
Urban Areas
Earth Sciences
Geography
Geographic Areas
Urban Areas
Earth Sciences
Geography
Human Geography
Urban Geography
Social Sciences
Human Geography
Urban Geography
Spatial and temporal dynamics of urban heat environment at the township scale: A case study in Jinan city, China
Spatial and temporal dynamics of urban heat environment at the township scale
https://orcid.org/0000-0002-3362-0631
Wang Dongchao Conceptualization Methodology Resources Validation Visualization Writing – original draft 1
Cao Jianfei Visualization Writing – original draft 1
https://orcid.org/0000-0001-8866-6728
Zhang Baolei Project administration Supervision Writing – original draft 1 *
Kong Kangning Writing – review & editing 2
Wang Run Writing – review & editing 2
1 College of Geography and Environment, Shandong Normal University, Jinan, Shandong, China
2 Shandong Provincial Territorial Spatial Ecological Restoration Center, Jinan, Shandong, China
Son Nguyen-Thanh Editor
National Central University, TAIWAN
Competing Interests: The authors have declared that no competing interests exist.

* E-mail: sdnuzbl@outlook.com
16 9 2024
2024
19 9 e030771110 1 2024
9 7 2024
© 2024 Wang et al
2024
Wang et al
https://creativecommons.org/licenses/by/4.0/ This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

The prolonged dependence on industrial development has accentuated the cumulative effects of pollutants. Simultaneously, influenced by land construction activities and green space depletion, the Urban Heat Island (UHI) effect in cities has intensified year by year, jeopardizing the foundation of sustainable urban development. Prudent urban spatial planning holds the potential to robustly ameliorate the persistent deterioration of the UHI phenomenon. This study selects Jinan City as a case study and employs spatial autocorrelation and spatial regression algorithms to explore the spatiotemporal evolution of urban-rural patterns at the township scale. The aim is to identify key factors driving the spatiotemporal differentiation of Land Surface Temperature (LST) from 2013 to 2022. The research reveals a trend of initially rising and subsequently falling LST in various townships, with low-temperature concentration areas in the southern mountainous region and the northern plain area. The "West-Central-East" main urban axis and the southeast Laiwu District exhibit high-temperature zones. Significant influences on LST are attributed to pollution levels, topographical factors, urbanization levels, and urban greenness. The global Moran’s Index for LST exceeds 0.7, indicating a strong positive spatial correlation. Cluster analysis results indicate High-High (HH) clustering in the central Shizhong District and Low-Low (LL) clustering in the northern Shanghe County. Multiscale Geographically Weighted Regression (MGWR) outperforms Geographically Weighted Regression (GWR) and Ordinary Linear Regression (OLR), providing a more accurate reflection of the regression relationships between variables. By investigating the spatiotemporal evolution of LST and its driving factors at the township scale, this study contributes insights for future urban planning and sustainable development.

Shandong Province Undergraduate Teaching Reform Project Z20220004 https://orcid.org/0000-0001-8866-6728
Zhang Baolei Jinan City-School Integration Project JNSX2023036 https://orcid.org/0000-0001-8866-6728
Zhang Baolei http://dx.doi.org/10.13039/501100001809 National Natural Science Foundation of China 42201308 https://orcid.org/0000-0001-8866-6728
Zhang Baolei Natural Science Foundation of Shandong Province, China ZR2021QD127, ZR2021ME203 https://orcid.org/0000-0001-8866-6728
Zhang Baolei B.Z. acknowledges financial support from Shandong Province Undergraduate Teaching Reform Project (Z20220004), Jinan City-School Integration Project (JNSX2023036), National Natural Science Foundation of China (42201308), and Natural Science Foundation of Shandong Province, China (ZR2021QD127, ZR2021ME203). Data AvailabilityThe authors acknowledge the LP DAAC for MODIS data, the EROS center for Landsat and DEM data, the ESRI for land use data, and the EC for GHSL data. The data used in this article are all public open source, and can be obtained from https://lpdaac.usgs.gov/products/mcd19a2v061/, https://www.usgs.gov/centers/eros/science/usgs-eros-archive-landsat-archives-landsat-8-9-olitirs-collection-2-level-2, https://www.usgs.gov/centers/eros/science/usgs-eros-archive-digital-elevation-global-multi-resolution-terrain-elevation, https://www.arcgis.com/home/item.html?id=cfcb7609de5f478eb7666240902d4d3d, and https://ghsl.jrc.ec.europa.eu/datasets.php.
Data Availability

The authors acknowledge the LP DAAC for MODIS data, the EROS center for Landsat and DEM data, the ESRI for land use data, and the EC for GHSL data. The data used in this article are all public open source, and can be obtained from https://lpdaac.usgs.gov/products/mcd19a2v061/, https://www.usgs.gov/centers/eros/science/usgs-eros-archive-landsat-archives-landsat-8-9-olitirs-collection-2-level-2, https://www.usgs.gov/centers/eros/science/usgs-eros-archive-digital-elevation-global-multi-resolution-terrain-elevation, https://www.arcgis.com/home/item.html?id=cfcb7609de5f478eb7666240902d4d3d, and https://ghsl.jrc.ec.europa.eu/datasets.php.
==== Body
pmcIntroduction

The intricate interplay of urban ecological environments and economic activities complicates the interpretation of urban spatial patterns, making them inherently complex and variable [1]. Urban patterns significantly alter pollution levels [2], natural environments [3], land use [4], and population distribution [5]. The rapid expansion of urban areas exacerbates negative resident sentiments [6] and triggers various issues related to transportation [7], environment [8], and development [9]. The escalating scale of economic activities within urban spaces results in a surge in urban population and impermeable surface area, while the proportion of vegetation and water bodies continuously diminishes [10]. These urban development challenges elevate urban Land Surface Temperature (LST), widening the temperature difference between urban and rural areas, thereby reinforcing the exposure risk of the Urban Heat Island (UHI) effect [11, 12]. As urban spatial patterns adjust with increasing urbanization, changes in urban characteristics (pollution levels, land use, architectural forms, and population size) significantly alter, indirectly promoting the increase in regional heat production and absorption capacity [13]. In the context of global warming, population expansion, and rapid urbanization, the UHI effect poses a substantial threat to current living standards and future sustainable development [14]. Studying the spatiotemporal patterns of the urban thermal environment and exploring the impact of urban spatial patterns on the thermal environment are crucial for alleviating the increasingly severe negative effects of the UHI.

LST, obtained through thermal radiation measurement techniques [15], is crucial for understanding the UHI effect [16]. Recent advancements in remote sensing technology have greatly expanded the perspectives and methods in the field of UHI. Scholars use remote sensing technology to identify urban impervious layers, establish regression models between impervious layers and LST, and further measure the intensity of the UHI effect [17]. Alternatively, they divide local climate zones [18] within cities according to different criteria, combine building landscape forms [19], and simulate the correlation between various urban indicators [20] and LST. Advanced machine learning techniques are also applied in the field of LST measurement [21], combined with disciplines such as urban climate, human geography, and urban planning, collectively promoting the research on the urban thermal environment [22].

The resolution and accuracy of remote sensing data are crucial for accurately capturing the spatiotemporal variations in the UHI effect [23]. However, current remote sensing technology still has room for improvement in terms of high resolution and high accuracy [24]. The continuous changes in factors such as urban structure, greenery rate, and land use during the process of urbanization have a significant impact on the formation and development of the UHI effect [25]. Therefore, comprehensive monitoring and analysis using multi-source, multi-temporal remote sensing data are necessary. However, current remote sensing technology faces technical challenges in obtaining long time series data and achieving large-area coverage [26]. While remote sensing technology provides information on LST, it needs to be combined with multidisciplinary data such as meteorology, geology, and ecology to comprehensively and deeply understand the mechanisms and impacts of the UHI effect [27]. In addition, collaborative efforts with urban planning and management are essential directions for future research [28]. In summary, although remote sensing technology has achieved some results in the field of the UHI effect and LST inversion, continuous improvements are needed in data resolution, spatiotemporal monitoring capabilities, and multi-source data fusion to better address environmental issues in the urbanization process and provide scientific support for urban sustainable development [29].

This study unlocks a perspective on urban spatial patterns at the township scale, employing spatial regression statistical models and spatial autocorrelation analysis tool to decipher complex urban spatial feature factors and delineate the spatiotemporal patterns and driving elements of the UHI effect. We quantify the spatiotemporal driving effect of township spatial patterns on urban LST from a new perspective, bridging the gap in the study of urban heat environment at the township scale. The main objectives of this paper are as follows: (1) preprocess spatiotemporal image sequences, use a mono-window algorithm to invert LST, and analyze its spatiotemporal patterns; (2) quantify the spatial dependence and spatial clustering characteristics of LST using spatial autocorrelation methods; (3) utilize geographically weighted spatial regression methods to explore the local spatial relationships between urban characteristics and LST.

Materials and methods

Study area

Jinan City is located in the central part of the Shandong Peninsula, with Mount Tai to the south and the Yellow River to the north (Fig 1). Geographically adjacent to the North China Plain, it experiences a typical temperate monsoon climate. Jinan comprises 12 districts and counties, namely, Shizhong (17 townships), Lixia (14 townships), Huaiyin (14 townships), Tianqiao (15 townships), Licheng (21 townships), Changqing (10 townships), Zhangqiu (20 townships), Jiyang (10 townships), Laiwu (15 townships), Gangcheng (5 townships), Pingyin (8 townships), and Shanghe (12 townships). As the capital city of Shandong Province, Jinan possesses substantial political, economic, industrial, and population resources. In recent years, with the continuous enhancement of the city’s overall strength, urban spatial patterns and functional zoning have undergone constant evolution and updates, significantly influencing the UHI effect.

10.1371/journal.pone.0307711.g001 Fig 1 Location and study area (Jinan City, China).

The expanding urban volume continually increases the demand for land and natural resources, leading to intensified pollutant emissions [30]. Townships are crucial administrative units for managing Chinese cities, capable of responding to various strategic urban development plans [31]. This study explores the spatiotemporal patterns of LST and its driving factors at the township scale, providing valuable references for future urban spatial planning.

Data

Data sources

The research dataset includes remote sensing images, digital elevation models, land use data, built-up area information, population data, and urbanization levels (Table 1). MODIS images provide aerosol thickness data for assessing urban air pollution. Landsat images provide surface reflectance data, and a mono-window algorithm is applied for LST inversion. Digital elevation models are obtained from the U.S. Geological Survey (USGS) Earth Resources Observation and Science (EROS) center, covering the entire globe. Land use data utilize the latest 10-meter resolution data released by ESRI (it did not include data prior to 2017, so the data of 2017 were used for 2013 and 2016 in the experiment). Built-up area, population, and urbanization level data are sourced from the Global Human Settlement Layer (GHSL) by the European Commission (EC).

10.1371/journal.pone.0307711.t001 Table 1 Data sources.

Data source	Data provider	Spatial resolution	Time	
MODIS (MCD19A2v061)	NASA LP DAAC at the USGS EROS Center [32]	1000 m	2013	
2016	
2019	
2022	
Landsat 8 OLI/TIRS	USGS EROS Center [33]	30 m	2013	
2016	
2019	
2022	
DEM	USGS EROS Center [34]	7.5 arc-second	2010	
ESRI 2020 Global Land Use Land Cover from Sentinel-2	Impact Observatory for Esri [35]	10 m	2017	
2019	
2022	
GHS built-up surface grid (GHS-BUILT-S)	Global Human Settlement Layer (GHSL), European Commission (EC), Joint Research Centre (JRC) [36]	100 m	2020	
GHS built-up height grid (GHS-BUILT-H)	GHSL, EC, JRC [37]	100 m	2018	
GHS population grid (GHS-POP)	GHSL, EC, JRC [38]	100 m	2020	
GHS urbanisation-degree grid (GHS-SMOD)	GHSL, EC, JRC [39]	1000 m	2020	

Explanatory variables

The relationship between urban thermal environments and urban spatial features is complex, influenced by changes in industry, economy, culture, and public health, all impacting the distribution of urban LST. To study the driving factors of the UHI effect, various types of factors need comprehensive consideration [40]. This study collects six types of factors: Environmental condition, Remote sensing indicators, Terrain, Land use, Building scale, and Socioeconomic factors, with a total of 16 variables (Table 2) [41–44]. Fig 2 illustrates the basic spatial distribution characteristics of explanatory variables at the township scale.

10.1371/journal.pone.0307711.g002 Fig 2 Spatial variations of explanatory variables related to LST.

10.1371/journal.pone.0307711.t002 Table 2 Variable explanations.

Type	Variable	Description	Formula	
Environmental condition	Aerosol Optical Depth (AOD)	Average of AOD per township	—	
Remote sensing indicators	Normalized Difference Vegetation Index (NDVI)	Average of NDVI per township	(Bnir ‐ Bred) / (Bnir + Bred)	
Normalized Difference Water Index (NDWI)	Average of NDWI per township	(Bgreen ‐ Bnir) / (Bgreen + Bnir)	
Normalized Difference Built-up Index (NDBI)	Average of NDBI per township	(Bmir ‐ Bnir) / (Bmir + Bnir)	
Built-up Index (BU)	Average of BU per township	NDBI ‐ NDVI	
Terrain	Elevation (HGT)	Average of elevation per township	—	
Slope	Average of slope per township	—	
Aspect	Average of aspect per township	—	
Land use	Proportion of water land (LU1)	Water proportion in each township	AREAw / AREAt	
Proportion of vegetation land (LU2)	Vegetation proportion in each township	AREAv / AREAt	
Proportion of cultivated land (LU3)	Cultivated land proportion in each township	AREAcl / AREAt	
Proportion of construction land (LU4)	Construction proportion in each township	AREAc / AREAt	
Building scale	Building height (BH)	Average of BH per township	—	
Building area (BA)	Average of BA per township	—	
Socioeconomic factors	Population count (PC)	Average of PC per township	—	
Urbanisation degree (UD)	Average of UD per township	—	

Prior to using these explanatory variables, a pre-detection process is necessary to address significant multicollinearity issues between variables. Table 3 lists the basic statistical information (Mean, Variance, and Standard deviation), correlations (Pearson correlation coefficient), global Moran’s Index, and variance inflation factors (VIF) for each variable. The results indicate that the VIF values for NDVI, NDWI, NDBI, and BU in remote sensing indicators are higher than 10.0, revealing significant multicollinearity. To correct this issue, a logarithmic transformation is applied to all 16 variables, eliminating interference between variables and ensuring the interpretability of regression results.

10.1371/journal.pone.0307711.t003 Table 3 Variable statistical description results.

Variable	Mean	Variance	Standard deviation	Pearson correlation coefficient	Global Moran’s Index	VIF	Logarithmized VIF	
AOD	0.421	0.002	0.049	0.430	0.121	9.535	7.275	
NDVI	0.422	0.017	0.131	-0.304	0.828	17.771	9.154	
NDWI	-0.448	0.015	0.123	0.290	0.911	13.451	8.385	
NDBI	-0.087	0.009	0.094	-0.031	0.367	13.525	8.426	
BU	-0.509	0.038	0.196	0.184	0.826	16.631	8.864	
HGT	166.918	27914.053	167.075	-0.383	0.273	8.127	7.203	
Slope	4.356	35.114	5.926	-0.368	0.216	9.835	8.965	
Aspect	172.126	12064.994	109.841	0.057	0.154	1.929	1.679	
LU1	0.040	0.000	0.000	0.017	0.220	4.431	2.171	
LU2	0.528	0.000	0.000	-0.360	0.248	8.609	6.092	
LU3	0.884	0.000	0.000	0.213	0.734	7.588	4.791	
LU4	0.542	0.000	0.000	0.306	0.805	9.166	8.491	
BH	0.477	2.107	1.451	0.192	0.875	7.362	6.395	
BS	0.132	0.000	0.000	0.282	0.840	9.951	8.519	
PC	8.566	1554.950	39.433	0.181	0.915	7.180	4.070	
UD	16.277	38.230	6.183	0.315	0.971	8.187	7.872	

Methods

This study begins at the township scale, using Landsat remote sensing images to invert LST and employing spatial autocorrelation and spatial regression methods to explore the spatiotemporal patterns of LST and its driving factors (Fig 3). First, the statistical mono-window algorithm is utilized to calculate the inversion products of LST based on the radiance brightness temperature and atmospheric parameters derived from satellite images. Then, several explanatory variables related to LST, such as aerosol optical thickness (AOD), digital elevation model (DEM), land use (LU), urbanization degree (UD), normalized difference vegetation index (NDVI), etc., are selected as independent variables for spatial regression analysis. Finally, global (OLR) and local (GWR, and MGWR) regression methods are employed to establish spatial relationship models between LST and explanatory variables, analyzing the spatiotemporal distribution characteristics and driving factors of LST.

10.1371/journal.pone.0307711.g003 Fig 3 Research technology flowchart.

LST inversion

The Landsat LST standard product often show striped chromatic aberration at mesoscale [45]. The problem is mainly caused by several factors, such as satellite optical system errors, image sensor errors, color processing errors, and the influence of ambient light [46]. It is necessary to repair the strip chromatic aberration present in the data when using it to eliminate the serious impact of data errors on the results of subsequent experiments [47].

The mono-window inversion algorithm for LST is a method used to extract LST information from remote sensing data [48]. Its core idea is to utilize the interaction of thermal radiation observation signals between the atmosphere and the land surface. By establishing a radiative transfer equation, the algorithm inverts the radiance brightness temperature of a single band or window, providing an estimate of LST. In recent years, many researchers have adopted this method for LST inversion to ensure data accuracy and stability [49–51].

In the mono-window algorithm, considerations typically include atmospheric transmittance, absorption, scattering, and emission factors to establish a complex radiative transfer model. By utilizing observed radiance brightness temperature data and mathematical optimization algorithms, model parameters are adjusted to make the model calculation results as consistent as possible with the observed data, obtaining the optimal estimate of LST. The specific algorithm is as follows: C=ετ (1)

D=(1−τ)[1+(1−ε)τ] (2)

Lλ=scale*DN+offset (3)

Ta=K2ln1+K1Lλ (4)

Ts=a(1−C−D)+Ta(b(1−C−D)+C+D)−DTbC−273.15 (5)

In the formula, ε represents the land surface emissivity; τ represents atmospheric transmittance; K1 and K2 are preset constants for Landsat-8, where K1 = 774.89, and K2 = 1321.08; a and b are fitting coefficients with values -67.355351 and 0.458606, respectively; C and D are intermediate variables; Lλ represents radiance intensity values, with gain parameter scale = 0.0003342, offset = 0.1, and digital number value DN representing pixel radiance; Ta represents brightness temperature, Tb represents atmospheric average temperature, and Ts is the target LST.

Spatial regression models

Geographically Weighted Regression (GWR) is a regression analysis method used for spatial data [52]. Unlike traditional global regression models, GWR allows the model’s parameters to vary spatially, better capturing heterogeneity and local correlations in geographic spatial data [53]. The regression model of GWR can be expressed as: Yi=βi0+∑k=1PβikXik+εi (6)

Where Yi is the observed value of the dependent variable at geographic location i; P corresponds to the number of independent variables; Xikis the observed value of independent variable k at geographic location i; βi0 and βik are the intercept and regression coefficient at geographic location i; εi is the error term.

The key point of GWR is the use of different regression coefficients for each geographic location. The spatial weight matrix reflects the variation of spatial dependence within the study area. This allows GWR to have the ability to detect spatial heterogeneity and improve the accuracy of regression fitting. GWR generally uses a specific distance decay function to determine the elements of the weight matrix, and the spatial bandwidth controls the decay rate. Samples closer in distance have larger weights, while those farther away have smaller weights.

Multiscale Geographically Weighted Regression (MGWR) is an improved spatial regression model that can be used to explore multiscale variations in spatial heterogeneity [54]. Compared to GWR, which uses a fixed global bandwidth, MGWR assigns an independent spatial bandwidth for each variable. Therefore, MGWR can flexibly capture spatial heterogeneity of different regression coefficients and more accurately reflect local relationships in spatial data [55]. Moreover, MGWR can consider both global and local information, providing a more flexible balance between global and local factors.

Spatial autocorrelation analysis

In urban thermal environment research, spatial autocorrelation is a method used to describe the spatial correlation of geographic data. Moran’s I is a commonly used spatial autocorrelation statistic, measuring the overall spatial correlation degree of geographical data [56]. It can be divided into global Moran’s Index and local Moran’s Index. The calculation formulas are as follows: GlobalMoran'sI=n∑i=1n∑j≠inWijxi−x¯xj−x¯∑i=1n∑j≠inWij∑i=1nxi−x¯2 (7)

LocalMoran'sI=nxi−x¯∑j=inWijxj−x¯∑i=inxi−x¯2 (8)

Where n is the number of geographic units; xi is the observed value of the ith geographic unit; x¯is the mean value of all geographic units; Wij is the spatial weight between geographic unit i and geographic unit j. The Moran’s Index ranges from -1 to 1, where positive values indicate positive correlation and negative values indicate negative correlation. Values close to 1 or -1 suggest significant spatial clustering patterns.

Local Indicators of Spatial Association (LISA) aggregation analysis is used to detect local spatial clustering patterns in geographic space. By calculating the Moran’s Index between each geographic unit and its neighboring units, geographic units with significant spatial clustering or dispersion patterns can be identified. LISA aggregation analysis produces four quadrants: High-High (HH), Low-Low (LL), High-Low (HL), and Low-High (LH). In the HH and LL quadrants, the observed values of geographic units exhibit clustering patterns, while in the HL and LH quadrants, they show dispersion patterns. This helps understand the local clustering characteristics in spatial distribution.

Analysis results

Spatiotemporal patterns of LST

Utilizing the mono-window inversion algorithm for LST, this study obtained spatiotemporal distribution maps of LST in Jinan for the years 2013–2022 (Fig 4). The daytime annual average LST in Jinan exhibited a continuous upward trend from 2013 to 2019, followed by a decline over the subsequent three years. Originating from the southwest in Pingyin County, the temperature rise extended in a northeastern direction, connecting Changqing District, Huaiyin District, Shizhong District, Tianqiao District, Licheng District, and Lixia District, ultimately spanning the "Western-Central-Eastern" main urban axis. This axis, along with the southeastern Laiwu District, characterized by flat terrain dominated by urban land use, industrial development, and high population density, demonstrated higher LST. In contrast, the southern mountainous regions and the northern plain areas, separated by the Yellow River from the city center, exhibited relatively lower temperatures, primarily featuring forests and farmland.

10.1371/journal.pone.0307711.g004 Fig 4 Spatiotemporal distribution characteristics of LST in 2013, 2016, 2019, and 2022.

This study conducted regional statistical analysis on LST raster data (Fig 4), yielding township scale results. In 2013, Tangye Street in Licheng District recorded the highest LST at 24.2°C, while Duozhuang Township in Zhangqiu District exhibited the lowest LST at 19.6°C. In 2016, Ancheng Township in Pingyin County and Xiying Township in Licheng District reported the highest LST at 25.4°C and the lowest LST at 18.1°C, respectively. By 2019, the peak LST shifted to Wande Township in Changqing District, reaching 26.8°C, while Xiying Township in Licheng District maintained the lowest LST at 21.0°C. In 2022, Jiefanglu Street in Lixia District revealed the highest LST at 23.6°C, while Wenyuan Street in Gangcheng District marked the lowest at 20.5°C. Throughout the nine years, Wande Township in Changqing District and Xiying Township in Licheng District documented the highest and lowest LSTs in Jinan, respectively.

Conducting a detailed spatial statistical analysis at the township level for each county in Jinan, Fig 5 illustrates the minimum, mean, and maximum LST for the four periods. Overall, a rising-then-falling temperature trend is observed. Lixia District, Shizhong District, and Pingyin County consistently exhibited higher average temperatures, while Shanghe County and Jiyang District displayed lower values. Examining the highest and lowest temperatures in different years, Changqing District and Pingyin County recorded the highest temperatures of 31.4°C and 28.5°C in 2019, whereas Licheng District and Gangcheng District experienced the lowest temperatures of 14.1°C and 17.3°C in 2016. Notably, Laiwu District, Zhangqiu District, Jiyang District, and Shanghe County demonstrated the most significant average temperature variations, with maximum increases of 3.1°C, 3.1°C, 3.0°C, and 2.8°C, respectively, over the nine-year period.

10.1371/journal.pone.0307711.g005 Fig 5 Statistical analysis of LST in counties of Jinan.

The numbers show the LST data; the larger the circle radius the higher the LST value; Min, Mean, and Max represent the minimum, mean, and maximum values for each county or district.

LST regression results and influencing factors

To visually compare the regression effects of different models on LST, the OLR, GWR, and MGWR models were successively applied to fit the temperature data. In this experiment, 16 explanatory variables from various townships in 2022 were selected as independent variables, with the inverted LST results serving as the dependent variable. Prior to the regression experiment, logarithmic transformations were applied to all variables to correct data distribution skewness, eliminate multicollinearity among variables, and enhance the explanatory power of the regression model.

Table 4 illustrates the regression analysis results for OLR, GWR, and MGWR models, with evaluation metrics including R2, adjusted R2, root mean square error (RMSE), and Akaike information criterion (AIC). Compared to OLR and GWR, the MGWR model achieved higher R2 and adjusted R2 values, along with lower RMSE and AIC levels. The respective values for MGWR were 0.872, 0.685, 0.521, and 329.840. Clearly, the MGWR model outperformed OLR and GWR in terms of fitting effectiveness, leading to its selection for explaining the regression relationships between variables.

10.1371/journal.pone.0307711.t004 Table 4 The regression analysis results for OLR, GWR, and MGWR models.

Type	Variable	OLR	GWR (Mean)	MGWR (Mean)	
	Constant	266.461	258.779	257.971	
Environmental condition	Log(AOD)	18.670	19.544	7.300	
Remote sensing indicators	Log(NDVI)	-579.959	-563.490	-549.528	
Log(NDWI)	10.399	10.073	3.060	
Log(NDBI)	599.257	586.002	580.586	
Log(BU)	-592.151	-577.388	-563.303	
Terrain	Log(HGT)	-0.001	-0.001	-0.009	
Log(Slope)	0.072	0.075	0.050	
Log(Aspect)	0.005	0.004	-0.005	
Land use	Log(LU1)	-0.660	-0.598	-2.431	
Log(LU2)	0.009	-0.033	-0.150	
Log(LU3)	0.097	0.107	0.389	
Log(LU4)	0.895	1.074	1.509	
Building scale	Log(BH)	0.038	0.063	0.064	
Log(BA)	-4.067	-5.087	-8.750	
Socioeconomic factors	Log(PC)	0.002	0.002	-0.018	
Log(UD)	0.059	0.067	0.006	
	R2	0.445	0.550	0.872	
	Adjusted R2 (R2 adj)	0.383	0.453	0.685	
	RMSE	0.774	0.651	0.521	
	AIC	408.192	399.205	329.840	

LST spatial autocorrelation analysis

Global and local spatial autocorrelation analyses were employed to explore the spatial heterogeneity of LST. Fig 6 presents Moran scatter plots illustrating the fitted relationship between remotely sensed LST and MGWR estimations for the four periods, establishing a foundation for quantifying the spatial differentiation of LST data. The results indicate that the global Moran indices for Jinan LST during 2013–2022 were 0.75, 0.71, 0.72, and 0.78, respectively, demonstrating a significant positive spatial correlation.

10.1371/journal.pone.0307711.g006 Fig 6 Moran scatter plots of LST in Jinan for various periods.

Fig 7 provides LISA distribution maps obtained through spatial clustering analysis for the four periods. These maps effectively explain the local spatial autocorrelation variations of LST at the township level. The spatial distribution patterns of LST during 2013–2022 generally remained similar. High-high correlation (HH) primarily occurred in the central urban areas, with occasional distribution in Lixia District, Zhangqiu District, Laiwu District, and Pingyin County. Low-low correlation (LL) was concentrated in Shanghe County, with scattered occurrences in Gangcheng District, Zhangqiu District, Laiwu District, and Jiyang District. Low-high correlation (LH) generally surrounded the high-high correlation (HH) areas, geographically close to the city center.

10.1371/journal.pone.0307711.g007 Fig 7 The spatial distribution clustering of LST in the years 2013, 2016, 2019, and 2022.

Discussion

Driving factors of urban heat island effect

Traditional approaches to uncovering drivers of LST often rely on a singular or a few factors [57], such as local climatic zones, urban landscapes, land use, visible light remote sensing indices, and impervious urban surfaces. However, the impact of the UHI effect extends beyond the city center and may affect remote rural areas or neighboring cities [58]. As the natural hub for residents’ production and life in the field, the normal functioning of urban systems depends on the coordinated interaction of various factors. Changes in air quality, topography, land use, socioeconomic conditions, and remote sensing indicators all contribute to LST in urban areas. Considering as many influencing factors as possible, the derived factors hold more scientific guidance significance.

Table 5 presents the diagnostic results of the correlation between LST and its influencing factors for four consecutive years. It is evident from the table that the environmental factor AOD consistently exhibits the strongest positive correlation with LST (r2013 = 0.455, r2016 = 0.400, r2019 = 0.386, r2022 = 0.430; P<0.001). Terrain factors such as HGT and Slope, as well as remote sensing indicators like NDVI, demonstrate a strong negative correlation with LST (r<-0.300; P<0.001). Significant positive correlations (r>0.300; P<0.001) are observed between land use type LU4, socioeconomic UD, and LST. However, building scale such as BH and BS show no significant positive correlation with LST (r<0.290; P<0.001). Based on the correlation analysis results of driving factors, it is determined that environmental conditions, topography, land use types, socioeconomic conditions, and remote sensing indicators are more critical than building scale.

10.1371/journal.pone.0307711.t005 Table 5 Correlation analysis between LST and variables in the years 2013, 2016, 2019, and 2022.

Variable	2013	2016	2019	2022	
r-value	P-value	r-value	P-value	r-value	P-value	r-value	P-value	
AOD	0.455	0.0000	0.400	0.0000	0.386	0.0001	0.430	0.0000	
NDVI	-0.300	0.0000	-0.304	0.0000	-0.306	0.0000	-0.304	0.0000	
NDWI	0.269	0.0000	0.195	0.0000	0.215	0.0000	0.290	0.0000	
NDBI	0.278	0.0000	-0.025	0.0000	0.156	0.0000	-0.031	0.0001	
BU	0.181	0.0000	0.112	0.0001	0.134	0.0000	0.184	0.0001	
HGT	-0.395	0.0001	-0.353	0.0000	-0.350	0.0000	-0.383	0.0000	
Slope	-0.359	0.0001	-0.326	0.0000	-0.364	0.0002	-0.368	0.0000	
Aspect	0.304	0.0000	0.036	0.0001	0.097	0.0001	0.057	0.0002	
LU1	-0.051	0.0002	-0.067	0.0002	-0.048	0.0000	0.017	0.0002	
LU2	-0.281	0.0002	-0.304	0.0000	-0.268	0.0001	-0.360	0.0000	
LU3	-0.073	0.0004	0.210	0.0000	0.185	0.0000	0.213	0.0000	
LU4	0.394	0.0000	0.381	0.0002	0.332	0.0000	0.306	0.0000	
BH	0.124	0.0000	0.127	0.0001	0.111	0.0000	0.192	0.0001	
BS	0.166	0.0000	0.183	0.0002	0.143	0.0000	0.282	0.0000	
PC	0.151	0.0000	0.113	0.0001	0.117	0.0000	0.181	0.0002	
UD	0.360	0.0000	0.310	0.0000	0.350	0.0001	0.315	0.0000	

This study conducted grouped statistical analysis of LST based on AOD, HGT, Slope, LU, UD, and NDVI. The results are presented in Fig 8 in the forms of cloud-rain, box, and scatter plots. Land use types include water bodies, vegetation, cropland, and built-up areas. Here, these four land types are amalgamated into subcategories of a single indicator, and their LSTs are uniformly subjected to statistical analysis. AOD is positively correlated with LST. Beyond an altitude of 120 meters, higher HGT is associated with lower LST. As the slope increases, LST decreases. Among land use types, water bodies exhibit the lowest LST, followed by vegetation, cropland, and built-up areas with the highest LST. Larger UD values are associated with higher LST, while NDVI shows a negative correlation with LST. These six factors drive LST variations in different ways, consequently influencing the UHI effect.

10.1371/journal.pone.0307711.g008 Fig 8 Regulatory effect of warming and cooling factors on LST.

Temporal resilience of urban sustainability

Temporal resilience of urban sustainability pertains to a city’s adaptive capacity over time to various shocks and changes, reflecting the degree of sustainable development in the face of challenges at different periods [59]. Emphasizing the dynamic, adaptive, and responsive aspects of urban sustainability over long-term evolution, temporal resilience underscores the need for cities to maintain stability and development across different periods. As urbanization accelerates, the UHI effect has become a decisive force in urban microclimates [60], progressively impinging on sustainable development. Through a correlation analysis of variables over time, sensitivity factors for LST can be identified. Fig 9 presented in the form of a segmented dumbbell plot, illustrates the changing strength of correlations between explanatory variables (AOD, NDVI, NDWI, NDBI, BU) and LST in 2013, 2016, 2019, 2020 and 2022.

10.1371/journal.pone.0307711.g009 Fig 9 Temporal changes in correlation strength between explanatory variables and LST.

Distinguished from the other four variables, the correlation dynamics of AOD appear unique, registering its lowest value in 2020 but surpassing the 2013 level in 2022. Plausibly, this divergence could be attributed to the significant global events during 2019–2022, particularly the outbreak of the COVID-19 pandemic at the end of 2019 [61]. Research has demonstrated substantial global impacts of the pandemic, including financial market volatility, production disruptions, and increased unemployment. The specific effects varied across countries, regions, and industries. During the pandemic, the reduced economic activity in Jinan likely led to decreased air pollution sources [62], contributing to an improvement in the self-purification ability of the environment. The phenomenon indirectly drives the urban LST fallback event from 2020 to 2022.

Limitations and future prospects

This study, relying on remote sensing data and other geographical data, employed spatial statistical methods and spatial autocorrelation analysis to explore the spatiotemporal characteristics and driving factors of LST at the township level. Nevertheless, there are several limitations and shortcomings in the content and methodology of the study. (1) The experiment’s acquired LST data represent daytime averages for the entire year, without considering temporal variations such as day-night differences, monthly variations, or seasonal changes. (2) The study exclusively addresses spatial scales at the township level, overlooking potential errors arising from spatial heterogeneity expression at different scales. (3) Due to limitations in data acquisition frequency, the availability of explanatory variables for all periods was constrained. The use of variables from a single period inadequately captures the real-time dynamics of temperature.

Future research should prioritize investigations into multiple time scales, multiple spatial scales, and real-time multi-source data. A comprehensive exploration of the spatiotemporal characteristics and driving factors of LST necessitates the combined application of global statistical methods, local statistical methods, and spatiotemporal statistical methods. A few points need to be noted when applying the research methodology of this study to other regions or scales. (1) The climatic type and intensity of human activities in the target area strongly influence the spatiotemporal distribution patterns of LST driving factors [63], and the degree of correlation will also vary geographically. (2) Different research scales yield different results. District-level analysis is suitable for large areas with sparse population densities [64], while Street-level analysis is more suitable for densely populated areas [65]. The choice of scale should be based on the specific research requirements. In conclusion, researchers need to consider data selection, preprocessing, and research objective setting to eliminate the negative impact of uncertainties on the feasibility of the study.

Conclusion

This study, employing the classical mono-window algorithm, inverted LST in Jinan city. Utilizing spatial regression statistical models and spatial autocorrelation analysis tool, it delineated the spatiotemporal patterns and driving factors of the UHI effect at the township level. The study’s main findings can be summarized as follows:

From 2013 to 2019, LST in various townships gradually increased, followed by a subsequent decline. The southern mountainous regions with dense vegetation and the northern plain areas primarily featuring farmland exhibited lower temperatures. The "Western-Central-Eastern" main urban axis, along with the southeastern Laiwu District characterized by urban land use, industrial development, and high population density, experienced higher LST.

During the period of 2013–2022, both global and local Moran indices for LST in Jinan were consistently high, affirming a significant positive spatial correlation. LISA results demonstrated that HH correlation mainly occurred in the city center, while LL correlation was concentrated in the northern regions. LH correlation typically surrounded the HH correlation areas, geographically close to the city center.

Factors driving LST are diverse, with environmental conditions exhibiting the strongest correlation, followed by topography and land use types. Subsequently, socioeconomic conditions and remote sensing indicators demonstrated correlation, while the correlation with building scale was weakest. In comparison to OLR and GWR models, the MGWR model, addressing the multiscale issue of regression coefficients, exhibited superior performance in local spatial regression.

In conclusion, this study provides valuable insights into the spatiotemporal dynamics and driving factors of LST in Jinan. The findings contribute to a more comprehensive understanding of the UHI effect, facilitating informed decision-making for sustainable urban development.
==== Refs
References

1 Jiang F , Ma J , Webster CJ , Chiaradia AJF , Zhou Y , Zhao Z , et al . Generative urban design: A systematic review on problem formulation, design generation, and decision-making. Prog Plan. 2023; 100795. doi: 10.1016/j.progress.2023.100795
2 Wu L , Chen C . Does pattern matter? Exploring the pathways and effects of urban green space on promoting life satisfaction through reducing air pollution. Urban For Urban Green. 2023;82 : 127890. doi: 10.1016/j.ufug.2023.127890
3 Goličnik Marušić B , Dremel M , Ravnikar Ž . A frame of understanding to better link nature-based solutions and urban planning. Environ Sci Policy. 2023;146 : 47–56. doi: 10.1016/j.envsci.2023.05.005
4 Gaur S , Singh R . A Comprehensive Review on Land Use/Land Cover (LULC) Change Modeling for Urban Development: Current Status and Future Prospects. Sustainability. 2023;15 : 903. doi: 10.3390/su15020903
5 Newsham N , Rowe F . Understanding trajectories of population decline across rural and urban Europe: A sequence analysis. Popul Space Place. 2023;29 : e2630. doi: 10.1002/psp.2630
6 He H , Sun R , Li J , Li W . Urban landscape and climate affect residents’ sentiments based on big data. Appl Geogr. 2023;152 : 102902. doi: 10.1016/j.apgeog.2023.102902
7 Khumvongsa K , Guo J , Theepharaksapan S , Shirakawa H , Tanikawa H . Uncovering urban transportation infrastructure expansion and sustainability challenge in Bangkok: Insights from a material stock perspective. J Ind Ecol. 2023;27 : 476–490. doi: 10.1111/jiec.13342
8 Zhang M , Kafy A-A , Xiao P , Han S , Zou S , Saha M , et al . Impact of urban expansion on land surface temperature and carbon emissions using machine learning algorithms in Wuhan, China. Urban Clim. 2023;47 : 101347. doi: 10.1016/j.uclim.2022.101347
9 Zhong C , Guo H , Swan I , Gao P , Yao Q , Li H . Evaluating trends, profits, and risks of global cities in recent urban expansion for advancing sustainable development. Habitat Int. 2023;138 : 102869. doi: 10.1016/j.habitatint.2023.102869
10 Zhou X , Wu B , Liu Y , Zhou Q , Cheng W . Synergistic effects of heat and carbon on sustainable urban development: Case study of the Wuhan Urban Agglomeration. J Clean Prod. 2023;425 : 138971. doi: 10.1016/j.jclepro.2023.138971
11 He T , Zhou R , Ma Q , Li C , Liu D , Fang X , et al . Quantifying the effects of urban development intensity on the surface urban heat island across building climate zones. Appl Geogr. 2023;158 : 103052. doi: 10.1016/j.apgeog.2023.103052
12 You M , Huang J , Guan C . Are New Towns Prone to Urban Heat Island Effect? Implications for Planning Form and Function. Sustain Cities Soc. 2023;99 : 104939. doi: 10.1016/j.scs.2023.104939
13 Dewa DD , Buchori I . Impacts of rapid urbanization on spatial dynamics of land use–based carbon emission and surface temperature changes in the Semarang Metropolitan Region, Indonesia. Environ Monit Assess. 2023;195 : 259. doi: 10.1007/s10661-022-10839-6 36595039
14 Liu X , Wang N , Li Z , Jia R , Qiao Z . Research on Time Series and Spatial Gradient of Urban Heat Island Expansion from the Perspective of Urban Renewal. IEEE J Sel Top Appl Earth Obs Remote Sens. 2023;16 : 8680–8688. doi: 10.1109/JSTARS.2023.3314970
15 Li Z-L , Wu H , Duan S-B , Zhao W , Ren H , Liu X , et al . Satellite Remote Sensing of Global Land Surface Temperature: Definition, Methods, Products, and Applications. Rev Geophys. 2023;61 : e2022RG000777. doi: 10.1029/2022RG000777
16 Wang C , Zhan W , Li L , Wang S , Wang C , Miao S , et al . Urban heat islands characterized by six thermal indicators. Build Environ. 2023;244 : 110820. doi: 10.1016/j.buildenv.2023.110820
17 Pan T , Li B , Ning L . Impervious Surface Area Patterns and Their Response to Land Surface Temperature Mechanism in Urban–Rural Regions of Qingdao, China. Remote Sens. 2023;15 : 4265. doi: 10.3390/rs15174265
18 Wang Z , Ishida Y , Mochida A . Effective Factors for Reducing Land Surface Temperature in Each Local Climate Zone Built Type in Tokyo and Shanghai. Remote Sens. 2023;15 : 3840. doi: 10.3390/rs15153840
19 Xu H , Li C , Hu Y , Li S , Kong R , Zhang Z . Quantifying the effects of 2D/3D urban landscape patterns on land surface temperature: A perspective from cities of different sizes. Build Environ. 2023;233 : 110085. doi: 10.1016/j.buildenv.2023.110085
20 Li K , Chen Y . Characterizing the indicator-based, day-and-night, and climate-based variations in response of surface urban heat island during heat wave across global 561 cities. Sustain Cities Soc. 2023;99 : 104877. doi: 10.1016/j.scs.2023.104877
21 Bindajam AA , Mallick J , Talukdar S , Shahfahad ,Shohan AAA , Rahman A . Modeling the spatiotemporal heterogeneity of land surface temperature and its relationship with land use land cover using geo-statistical techniques and machine learning algorithms. Environ Sci Pollut Res. 2023;30 : 106917–106935. doi: 10.1007/s11356-022-23211-5 36178650
22 Wang M , Lu H , Chen B , Sun W , Yang G . Fine-Scale Analysis of the Long-Term Urban Thermal Environment in Shanghai Using Google Earth Engine. Remote Sens. 2023;15 : 3732. doi: 10.3390/rs15153732
23 Leichtle T , Kühnl M , Droin A , Beck C , Hiete M , Taubenböck H . Quantifying urban heat exposure at fine scale ‐ modeling outdoor and indoor temperatures using citizen science and VHR remote sensing. Urban Clim. 2023;49 : 101522. doi: 10.1016/j.uclim.2023.101522
24 Reiners P , Sobrino J , Kuenzer C . Satellite-Derived Land Surface Temperature Dynamics in the Context of Global Change—A Review. Remote Sens. 2023;15 : 1857. doi: 10.3390/rs15071857
25 Deng X , Gao F , Liao S , Liu Y , Chen W . Spatiotemporal evolution patterns of urban heat island and its relationship with urbanization in Guangdong-Hong Kong-Macao greater bay area of China from 2000 to 2020. Ecol Indic. 2023;146 : 109817. doi: 10.1016/j.ecolind.2022.109817
26 Tahooni A , Kakroodi AA , Kiavarz M . Monitoring of land surface albedo and its impact on land surface temperature (LST) using time series of remote sensing data. Ecol Inform. 2023;75 : 102118. doi: 10.1016/j.ecoinf.2023.102118
27 Zhang Y , Yang Y , Pan X , Ding Y , Hu J , Dai Y . Multiinformation Fusion Network for Mapping Gapless All-Sky Land Surface Temperature Using Thermal Infrared and Reanalysis Data. IEEE Trans Geosci Remote Sens. 2023;61 : 1–15. doi: 10.1109/TGRS.2023.3269622
28 Al Shawabkeh R , AlHaddad M , Al-Fugara A , Al-Hawwari L , Al-Hawwari MI , Omoush A , et al . Modeling the impact of urban land cover features and changes on the land surface temperature (LST): The case of Jordan. Ain Shams Eng J. 2024;15 : 102359. doi: 10.1016/j.asej.2023.102359
29 Rajagopal P , Priya RS , Senthil R . A review of recent developments in the impact of environmental measures on urban heat island. Sustain Cities Soc. 2023;88 : 104279. doi: 10.1016/j.scs.2022.104279
30 Qu Y , Zhan L , Zhang Q , Si H , Jiang G . Towards sustainability: The impact of the multidimensional morphological evolution of urban land on carbon emissions. J Clean Prod. 2023;424 : 138888. doi: 10.1016/j.jclepro.2023.138888
31 Wang L. China’s new town movements since 1949: A state/space perspective. Prog Plan. 2022;155 : 100514. doi: 10.1016/j.progress.2020.100514
32 Lyapustin A , Wang Y . MODIS/Terra+Aqua Land Aerosol Optical Depth Daily L2G Global 1km SIN Grid V061. NASA EOSDIS Land Processes Distributed Active Archive Center; 2022. doi: 10.5067/MODIS/MCD19A2.061
33 Earth Resources Observation And Science (EROS) Center. Collection-2 Landsat 8–9 OLI (Operational Land Imager) and TIRS (Thermal Infrared Sensor) Level-2 Science Products. U.S. Geological Survey; 2013. doi:10.5066/P9OGBGM6
34 Earth Resources Observation And Science (EROS) Center. Global Multi-resolution Terrain Elevation Data 2010 (GMTED2010). U.S. Geological Survey; 2017. doi:10.5066/F7J38R2N
35 Karra K , Kontgis C , Statman-Weil Z , Mazzariello JC , Mathis M , Brumby SP . Global land use / land cover with Sentinel 2 and deep learning. 2021 IEEE International Geoscience and Remote Sensing Symposium IGARSS. Brussels, Belgium: IEEE; 2021. pp. 4704–4707. doi: 10.1109/IGARSS47720.2021.9553499
36 Pesaresi M. GHS-BUILT-S R2023A - GHS built-up surface grid, derived from Sentinel2 composite and Landsat, multitemporal (1975–2030). European Commission, Joint Research Centre (JRC); 2023. doi: 10.2905/9F06F36F-4B11-47EC-ABB0-4F8B7B1D72EA
37 Pesaresi M. GHS-BUILT-H R2023A - GHS building height, derived from AW3D30, SRTM30, and Sentinel2 composite (2018). European Commission, Joint Research Centre (JRC); 2023. doi:10.2905/85005901-3A49-48DD-9D19-6261354F56FE
38 Schiavina M , Freire S , MacManus K . GHS-POP R2023A - GHS population grid multitemporal (1975–2030). European Commission, Joint Research Centre (JRC); 2023. doi: 10.2905/2FF68A52-5B5B-4A22-8F40-C41DA8332CFE
39 Schiavina M , Melchiorri M , Pesaresi M . GHS-SMOD R2023A - GHS settlement layers, application of the Degree of Urbanisation methodology (stage I) to GHS-POP R2023A and GHS-BUILT-S R2023A, multitemporal (1975–2030). European Commission, Joint Research Centre (JRC); 2023. doi: 10.2905/A0DF7A6F-49DE-46EA-9BDE-563437A6E2BA
40 Liu K , Li X , Wang S , Li Y . Investigating the impacts of driving factors on urban heat islands in southern China from 2003 to 2015. J Clean Prod. 2020;254 : 120141. doi: 10.1016/j.jclepro.2020.120141
41 Taripanah F , Ranjbar A . Quantitative analysis of spatial distribution of land surface temperature (LST) in relation Ecohydrological, terrain and socio- economic factors based on Landsat data in mountainous area. Adv Space Res. 2021;68 : 3622–3640. doi: 10.1016/j.asr.2021.07.008
42 Hussain S , Mubeen M , Ahmad A , Majeed H , Qaisrani SA , Hammad HM , et al . Assessment of land use/land cover changes and its effect on land surface temperature using remote sensing techniques in Southern Punjab, Pakistan. Environ Sci Pollut Res. 2023;30 : 99202–99218. doi: 10.1007/s11356-022-21650-8 35768713
43 Warth G , Braun A , Assmann O , Fleckenstein K , Hochschild V . Prediction of Socio-Economic Indicators for Urban Planning Using VHR Satellite Imagery and Spatial Analysis. Remote Sens. 2020;12 : 1730. doi: 10.3390/rs12111730
44 Ezimand K , Azadbakht M , Aghighi H . Analyzing the effects of 2D and 3D urban structures on LST changes using remotely sensed data. Sustain Cities Soc. 2021;74 : 103216. doi: 10.1016/j.scs.2021.103216
45 Emery W , Camps A . Introduction to Satellite Remote Sensing: Atmosphere, Ocean, Land and Cryosphere Applications. Elsevier; 2017. Available: 10.1016/C2015-0-04517-8
46 Chen S , Ren H , Ye X , Dong J , Zheng Y . Geometry and adjacency effects in urban land surface temperature retrieval from high-spatial-resolution thermal infrared images. Remote Sens Environ. 2021;262 : 112518. doi: 10.1016/j.rse.2021.112518
47 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. doi: 10.1016/j.rse.2014.02.003
48 Wan J , Zhu M , Ding W . Accuracy Evaluation and Parameter Analysis of Land Surface Temperature Inversion Algorithm for Landsat-8 Data. Adv Meteorol. 2021;2021 : e9917145. doi: 10.1155/2021/9917145
49 Wang S , Luo Y , Li X , Yang K , Liu Q , Luo X , et al . Downscaling Land Surface Temperature Based on Non-Linear Geographically Weighted Regressive Model over Urban Areas. Remote Sens. 2021;13 : 1580. doi: 10.3390/rs13081580
50 Chao J , Zhao Z , Lai Z , Xu S , Liu J , Li Z , et al . Detecting geothermal anomalies using Landsat 8 thermal infrared remote sensing data in the Ruili Basin, Southwest China. Environ Sci Pollut Res. 2023;30 : 32065–32082. doi: 10.1007/s11356-022-24417-3 36462073
51 Al-Ruzouq R , Shanableh A , Khalil MA , Zeiada W , Hamad K , Abu Dabous S , et al . Spatial and Temporal Inversion of Land Surface Temperature along Coastal Cities in Arid Regions. Remote Sens. 2022;14 : 1893. doi: 10.3390/rs14081893
52 Brunsdon C , Fotheringham AS , Charlton ME . Geographically Weighted Regression: A Method for Exploring Spatial Nonstationarity. Geogr Anal. 1996;28 : 281–298. doi: 10.1111/j.1538-4632.1996.tb00936.x
53 Gao Y , Zhao J , Han L . Exploring the spatial heterogeneity of urban heat island effect and its relationship to block morphology with the geographically weighted regression model. Sustain Cities Soc. 2022;76 : 103431. doi: 10.1016/j.scs.2021.103431
54 Fotheringham AS , Yang W , Kang W . Multiscale Geographically Weighted Regression (MGWR). Ann Am Assoc Geogr. 2017;107 : 1247–1265. doi: 10.1080/24694452.2017.1352480
55 Yang L , Yu K , Ai J , Liu Y , Yang W , Liu J . Dominant Factors and Spatial Heterogeneity of Land Surface Temperatures in Urban Areas: A Case Study in Fuzhou, China. Remote Sens. 2022;14 : 1266. doi: 10.3390/rs14051266
56 Fan C , Wang Z . Spatiotemporal Characterization of Land Cover Impacts on Urban Warming: A Spatial Autocorrelation Approach. Remote Sens. 2020;12 : 1631. doi: 10.3390/rs12101631
57 Shi H , Xian G , Auch R , Gallo K , Zhou Q . Urban Heat Island and Its Regional Impacts Using Remotely Sensed Thermal Data—A Review of Recent Developments and Methodology. Land. 2021;10 : 867. doi: 10.3390/land10080867
58 Wang Z-H . Reconceptualizing urban heat island: Beyond the urban-rural dichotomy. Sustain Cities Soc. 2022;77 : 103581. doi: 10.1016/j.scs.2021.103581
59 Amirzadeh M , Sobhaninia S , Sharifi A . Urban resilience: A vague or an evolutionary concept? Sustain Cities Soc. 2022;81 : 103853. doi: 10.1016/j.scs.2022.103853
60 Bherwani H , Singh A , Kumar R . Assessment methods of urban microclimate and its parameters: A critical review to take the research from lab to land. Urban Clim. 2020;34 : 100690. doi: 10.1016/j.uclim.2020.100690
61 Acharya P , Barik G , Gayen BK , Bar S , Maiti A , Sarkar A , et al . Revisiting the levels of Aerosol Optical Depth in south-southeast Asia, Europe and USA amid the COVID-19 pandemic using satellite observations. Environ Res. 2021;193 : 110514. doi: 10.1016/j.envres.2020.110514 33245884
62 Nichol JE , Bilal M , Ali MA , Qiu Z . Air Pollution Scenario over China during COVID-19. Remote Sens. 2020;12 : 2100. doi: 10.3390/rs12132100
63 Geng X , Zhang D , Li C , Yuan Y , Yu Z , Wang X . Impacts of climatic zones on urban heat island: Spatiotemporal variations, trends, and drivers in China from 2001–2020. Sustain Cities Soc. 2023;89 : 104303. doi: 10.1016/j.scs.2022.104303
64 Stow DA , Weeks JR , Shih H , Coulter LL , Johnson H , Tsai Y-H , et al . Inter-regional pattern of urbanization in southern Ghana in the first decade of the new millennium. Appl Geogr. 2016;71 : 32–43. doi: 10.1016/j.apgeog.2016.04.006
65 Song J , Tong X , Wang L , Zhao C , Prishchepov AV . Monitoring finer-scale population density in urban functional zones: A remote sensing data fusion approach. Landsc Urban Plan. 2019;190 : 103580. doi: 10.1016/j.landurbplan.2019.05.011
