Study on the influencing factors and the spatiotemporal heterogeneity of Urban Heat Island effect in Nanchang City of China

ABSTRACT Although some researchers have explored spatiotemporal characteristics of the land surface temperature (LST) and analyzed how the influencing factors impact the LST, the comprehensive influencing factors and specific spatial relationships that explore the Urban Heat Island (UHI) effect from long time scale have rarely received widespread attention. Taking Nanchang as an example, this paper first uses Landsat data to retrieve the LST in years 2000, 2010, and 2019. Then, the Geographically Weighted Regression (GWR) model is utilized to investigate the various influencing factors. The results showed as follows: (1) The UHI effect was intensified in the past 20 years in Nanchang City. From 2000 to 2019, the “heat island” increased from 63.06% to 64.07%, while the “green island” decreased from 36.94% to 35.93%. (2) The spatiotemporal anisotropy of the UHI effect was significant. (3) The GWR model could better reveal the spatial heterogeneity of the factors affecting the UHI effects. The relationship between Soil Brightness Index (NDSI), Normalized Difference Built-up Index (NDBI), and Night-lighting (NL) and the relative UHI effect intensity had significant spatial non-stationarity. This study can provide some references from reasonable urban planning to alleviate the urban heat island.


Introduction
Urbanization is a complex and dynamic process generated by the interaction between socioeconomic and environmental factors (Mu et al. 2016;Zhang et al. 2020;Zhou and Chen 2018). This process not only results in significant changes in the characteristics of the underlying surface of the city, but also produces a series of disadvantageous impacts on the urban environment, especially the urban thermal environment (Bai et al. 2017). The Urban Heat Island (UHI) effect, which refers to the phenomenon of urban areas that exhibits typically warmer than surrounding rural areas (Levermore et al. 2018), is one of the core representatives of the urban thermal environment (Nwakaire et al. 2020). Multiple studies have documented that the UHI effect does negatively impact on natural and socioeconomic aspects of cities (Li et al. 2019 ;Li et al. 2020;Zhou et al. 2011), e.g., air pollution (Sabrin et al. 2020;Wei et al. 2021), climatic warming (Emmanuel and Krüger 2012), and health risks for urban residents (Tan et al. 2010). Therefore, the study on the spatial pattern of the UHI effect and its influencing factors has become a research hotspot in various fields, such as climatology, ecology, and urban planning (Chen et al. 2021;Lu et al. 2021).
Traditionally, UHI surveys have focused on air UHI, where the air temperature is measured in situ at fixed isolated stations (e.g., meteorological stations, traverses of vehicles) and more special platforms (e.g., tall towers, radiosondes, and aircraft) (Zhao et al. 2018). Due to the restriction of the scarcity and inhomogeneity distribution of monitoring stations, the air UHIs measured by monitoring stations usually cannot provide enough effective spatial and temporal coverage for urban thermal environment research (Dewan et al. 2021;Naserikia et al. 2019). In contrast, the surface UHI (SUHI) is primarily measured by the Land Surface Temperature (LST) obtained from Remote Sensing (RS) data of thermal infrared (TIR), indicating the radiative temperature difference between urban and suburban surfaces (Chen et al. 2021). TIR remote sensing has the characteristics of space-time coverage, high accuracy, and low cost, which can greatly improve the ability to study urban thermal environment (SINGH 1984;Ulpiani 2021). Therefore, it is helpful to identify the spatiotemporal pattern of the UHI effect based on SUHI data (Athukorala and Murayama, 2021;Pathak et al. 2021), and further explore its influencing factors (Monteiro et al. 2021;Siqi and Yuhong 2020), which provides a scientific basis for better management of the urban thermal environment in metropolitan cities (Emmanuel and Krüger 2012;Estoque and Murayama 2017;Zhao et al., 2018;Zhi et al. 2020).
It has been proved that many factors can affect the UHI effect. These factors include urban morphology (two-dimensional and three-dimensional) (Guo et al. 2020;Qiao et al. 2020;Takebayashi and Senoo 2018;Yang et al. 2021), land cover (e.g., water bodies, and green plants) (Bao et al. 2016;Siqi and Yuhong 2020), RS spectral indices (e.g., Normalized difference vegetation index (NDVI), Normalized Difference Built-up Index (NDBI), Soil Brightness Index (NDSI)) (Estoque and Murayama 2015;Guha, Govil, and Besoya 2020;Pathak et al. 2021;Sun 2016;Weng, Lu, and Schubring 2004), topography Kashki et al. 2021), and spatialized socioeconomic indicators, e.g., night-lighting (Kurniati and Nitivattananon 2016;Zhang, Estoque, and Murayama 2017;Zhang and Li 2018). However, most of the current studies only focused on a single factor rather than from a comprehensive perspective of physical and socioeconomic factors' points of view (Kashki et al. 2021;Pathak et al. 2021). Therefore, it is of great scientific significance to use an integrated model with comprehensive influencing factor to clarify the UHI effect intensity or the LST.
Traditional regression models assume that the variables are spatially homogeneous, the spatial non-stationarity of spatial variables in geoscience is often ignored (Deilami, Kamruzzaman, and Liu 2018;Siqi and Yuhong 2020). However, there really exists obvious spatial differences in spatial variables (physical, ecological, and socioeconomic) in cities. Spatial heterogeneity is composed of the spatial differentiation of physics, ecology, and socioeconomic in urban areas , which leads to the spatial non-stationarity of the formation mechanism of the UHI effect . Thus, in order to better understand the spatial differences, the Geographically Weighted Regression (GWR) model was borrowed to simulate spatial heterogeneity (Zhao et al. 2018;Zhi et al. 2020). Luo et al. (Luo and Peng 2016) used the GWR model to quantify the relationships of the LST and impact factors at different spatial resolutions in Chongqing City, China, and indicated that GWR provides a better fit and is more localized in the formation. Kashki et al. (Kashki et al. 2021) utilized the GWR model to test the spatial non-stationarity of the relationships between geographical parameters and the LST in Shiraz City, Iran, and found that GWR has better performance in showing the spatial distribution of LSTs. Therefore, GWR can better reveal the spatial variation in the relationship between LST or UHI effect and influencing factors from a local perspective, which is of great significance to rational urban planning. Current research has mainly focused on large cities and megacities with rapid economic development (Dewan et al. 2021;Yang, Huang, and Tang 2019;Zeng et al. 2009). However, due to factors, such as geography, history and policy, the research on medium and small cities is seriously insufficient (Li et al. 2019). In addition, with the rapid development of social economy, continuous urbanization, population growth, and inappropriate urban planning have gradually resulted in increasing UHI effect in medium and small cities . For example, Liao et al. (Liao et al. 2022) analyzed the combined effect the abnormal heat island and UHI in Guiyang from multiple time series, and found that the impact of UHI intensified with rapid urbanization. Therefore, it is significant and necessary to explore the spatiotemporal distribution of the UHI effect from multiple time series and analyze the influencing factors in order to alleviate the thermal environment problems in medium and small cities.
Taking Nanchang City, the capital city of Jiangxi Province with rapid urbanization in the middle reaches of the Yangtze River in China as an example, this paper first retrieved the LST through a single-window algorithm Murayama 2015, 2017) under the support of Landsat RS images (Landsat 5 TM and Landsat 7 ETM+) in year 2000, 2010, and 2019. Then the Geographic Information System (GIS), RS, and GWR models were utilized to explore the spatiotemporal characteristics of the UHI effect and influencing factors in three typical periods. The objectives of this study are as follows: (1) Deriving and analyzing the UHI effect spatiotemporal distribution characteristics in Nanchang from satellite images. (2) Identifying spatiotemporal anisotropy of the UHI effect in some directions. (3) Investigating the various influencing factors using the GWR model and discovering the spatial heterogeneity of influencing factors.
The results of this study are expected to clarify the temporal and spatial variations of the UHI effect in Nanchang and influencing factors behind the variations, thus providing some references for rational urban planning for sake of thermal environment regulation and livable environment construction.

Research area
Nanchang city is located between 28°10'~29°11' N and 115°27'~116°35ʹE on the southwest bank of Poyang Lake, the largest freshwater lake in China. The Ganjiang River passes through the city. It is a relatively plain area with a subtropical monsoon climate characterized by abundant heat, rainfall, and light. The extreme maximum temperature is 40.9°C. It is one of the hottest cities in China, wellknown as the "four furnaces" of the country. It experienced rapid urbanization in the past 40 years, and the population has increased rapidly from 3.0682 million in 1978 to more than 5.606 million in 2019 (http://tjj.nc.gov. cn/ncstjj/tjnj/nav_list.shtml). The population growth and city expansion have brought about a series of environmental problems, one of which is the UHI effect has been exacerbated, presenting substantial adverse threats to both residential and environmental health in and around cities (Nwakaire et al. 2020;Zhang, Estoque, and Murayama 2017).
As the earliest built-up subway, Line 1 runs through the east and west areas of urban development. The midpoint of it passing through the Ganjiang River happens to be the city center. The area with a radius of 25 km covers the city's concentrated development area (e.g., Donghu District, Xihu District, Qingyunpu District, Qingshanhu District, and Honggutan District). Therefore, this paper takes the midpoint of subway line 1 with a radius of 25 km as the research area ( Figure 1).

Data source and preprocessing
The data used in this study were downloaded from the open dataset (e.g., https://earthexplorer.usgs.gov/, https://earthdata.nasa.gov/). They included the image data of Landsat with the path of 121 and row of 40 on 2 November 2000 (Landsat 5 TM), 6 November 2010, and 15 November 2019 (Landsat 7 ETM+), which are clear and visible with minimal cloud cover less than 5.00%, and thus meet the requirements of this research. Such data preprocessing as radiometric calibration, atmospheric correction, geometric correction, masking, and other treatments on image data from three time periods were conducted by ENVI 5.3 and ArcGIS10.2. LSTs were retrieved utilizing the single-window algorithm for better discovering the spatial pattern characteristics of the UHI effect on the study area Murayama 2015, 2017). The vegetation index, built-up index, and bare land (soil) indexes were calculated to represent the land cover information for exploring the response relationship between the land cover and the UHI effect (Guo et al. 2020;Luo and Peng 2016).
In addition, Night-lighting data revealing the relationship between various relevant socioeconomic information and the UHI effect (Zhang and Li 2018) were obtained from the Earth Observation Group (https://payneinstitute.mines.edu/). The Digital Elevation Model (DEM) data were obtained from the Geospatial Data Cloud (http://www.gscloud.cn/) and the United States Geological Survey (https://earthex plorer.usgs.gov/). The elevation and slope were extracted by the DEM to display the influence of topography on the UHI effect (Kashki et al. 2021). Administrative boundary data were acquired from the Nanchang Natural Resources Bureau. UTM projection with the WGS84 coordinate system was applied to all spatial data. The resolution was unified to that of image data, which is 30 m × 30 m.

Methods
Firstly, the LST was retrieved through the singlewindow algorithm and the relative UHI effect intensity was classified. Then, four directions profile maps, namely Along Ganjiang River to Poyang Lake, Honggutan district to Yaohu Forest Park, Meiling National Forest Park to Nanchang County, and Meiling National Forest Park to Poyang Lake were drawn using the profile analysis method, and the spatiotemporal anisotropy of the UHI effect was explored. Finally, 1152 grid samples, nearly 10% of the total, were randomly selected for obtaining influencing factors values and the relative UHI effect intensity values. Taking the various influencing factors as the independent variable and the relative UHI effect intensity as the dependent variable, the GWR models were established ( Figure 2). The main methods used in this study are as follows.

Classification of UHI effect
The relative UHI effect intensity was employed to classify the UHI effect based on existing researches (Bao et al. 2016). The larger the value is, the more obvious the UHI effect is. The equation is as follows: Where, UHIEI is the relative UHI effect intensity, which is a dimensionless value that reflects the strength of the UHI effect; T i is the LST of the ith point in the study area; T mean is the mean LST of the study area.

Spatiotemporal anisotropy of the UHI effect
Based on the Third Law of Geography, the similarities and differences in spatial localization or stratification of different areas correspond to the similarities and differences of target variables in these areas (Zhi et al. 2020).
Profile map was applied to determine the temporal and spatial anisotropy of the UHI effect in this study since it can intuitively describe the spatiotemporal differentiation of the urban thermal environment (Zhao et al., 2018). According to the terrain gradient and urban expansion directions of Nanchang City, the multi-directional profile map on the UHIEI in three periods including Along Ganjiang River to Poyang Lake, Honggutan district to Yaohu Forest Park, Meiling National Forest Park to Nanchang County, and Meiling National Forest Park to Poyang Lake were made. Firstly, the mean UHIEIvalue was calculated from a set of polygon grids captured on UHIEI maps in year 2000, 2010, and 2019. Secondly, the starting location of the profile analysis was located, which was called kilometer 0. Finally, the profile map was drawn based on its own distance interval of 3 km and the UHIEI-value ( Figure 4).

Influencing factors of the UHI effect
The UHI effect is strongly affected by combined factors, such as DEM, NDVI, NDBI (Deilami, Kamruzzaman, and Liu 2018;Guo et al. 2020). NDVI, NDBI, and NDSI are usually used to represent vegetation, built-up land, and bare land information on urban underlying surfaces, and then DEM and slope are used to characterize topographic information to explore the impact on land covers and topography on UHI effect (Ulpiani 2021). In addition, Night-lighting is an important indicator to reflect the heat release of human activities (Wei et al. 2021). Totally six impact indexes were selected as the UHIEI explanatory variables in this study. The details are as follows (Table 1).
(1) Topography factors. Slope is also a critical topographic factor affecting the LST (Kashki et al. 2021). When the elevation increases by 10 m, the temperature decreases by about 0.064°C. (2) Vegetation factors. Vegetation is expected to strongly affect land surface temperature (Sun 2016). The NDVI was applied to express vegetation. Its expression is as follows: (3) Building factors. The built-up or impervious surfaces were expressed using a sensitive indicator of the NDBI in this work. When the negative NDBI value usually represents the water body, close to zero value represents the vegetation, and a positive   (Zhang, Estoque, and Murayama 2017). The equation is as follows: (4) Bare land (soil) factors. The NDSI, can better distinguish the construction land from the bare land (Luo and Peng 2016) was used in this study. The equation is as follows: (5) Socioeconomic factors. Night light data have been widely used in urban research, such as analyzing urban spatial sprawl trends, mapping urban areas, and estimating the urban thermal environment (Guo et al. 2020;Zhang and Li 2018). The NL obtained by Inverse Distance Kriging Interpolation was applied to explain the relationship with the UHIEI.
In the above formula, NIR is the near-infrared band (band 4 of Landsat 5 TM and Landsat 7 ETM+), Red is the red band (band 3 of Landsat 5 TM and Landsat 7 ETM+), Green is the green band (band 2 of Landsat 5 TM and Landsat 7 ETM+), Blue is the blue band (band 1 of Landsat 5 TM and Landsat 7 ETM+), and MIR represents the midinfrared band (band 5 of Landsat 5 TM and Landsat 7 ETM+).

Regression analysis
Since the GWR model can reflect the spatial non-stationarity through the variation of the spatial position of the regression coefficients (Luo and Peng 2016), it is often used to analyze the spatial heterogeneity of explanatory variables. Based on the former research (Zhi et al. 2020) and actual conditions in the study area, the study area was divided into 420 m × 420 m grid cells. 1152 grid samples, approximately 10% of the total, were randomly selected to obtain the mean value of influencing factors and the UHIEI. In order to ensure the feasibility and improve the accuracy of the GWR model, the traditional ordinary linear least squares (OLS) model was used to conduct a preliminary test on the influencing factors. The variance inflation factor (VIF) tested was implemented to ensure there exists no multicollinearity between variables. The GWR model is represented as follows: where y i is the observed value(UHIEI); (μ i , ν i ) are the geographic coordinates of data sampling point i; β 0 (μ i , ν i ) is the regression attempt for sample i; β k (μ i , ν i ) is the regression parameter for independent variable k at data sampling grid i; k is the number of independent variables; x ik is the value of the independent variable x k (i.e., DEM, Slope, NDVI, NDBI, NDSI, NL) at point i; ε i is the random error term.

Spatiotemporal distribution characteristics of UHI effect
According to previous research (Bao et al. 2016), the UHIEI can be divided into six levels, i.e., Green Island Effect (GIE), Weak Green Island Effect (WGIE), Medium Heat Island Effect (MHIE), Relatively Strong Heat Island Effect (RSHIE), Strong Heat Island Effect (SHIE) and Extremely Strong Heat Island Effect (ESHIE) ( Table 2). As can be seen from Figure 3, spatially, the ESHIE was mainly distributed in the main urban area with dense population and increasing built-up land. The RSHIE was mainly located in the areas bordering the suburbs and urban areas. The GIE was mainly distributed in such nature conversation land as the Meiling National Forest Park, Yaohu Forest Park, and rivers and lakes' water surface.
From 2000 to 2019, we found that the "heat island" area of Nanchang City had increased from 1238.18 km 2 to 1258.02 km 2 , while the "Green Island" area had decreased from 725.31 km 2 to 705.47 km 2 (Table 3). In general, except for SHIE, the proportions of the other effect have undergone significant changes between 2000 and 2019. The area of the ESHIE, RSHIE, and GIE increased by 27.50, 31.04, and 58.39 km 2 , while that of the MHIE and WGIE decreased by 25.17 and 78.23 km 2 .

Spatiotemporal anisotropy of the UHI effect
From the profile map of the UHIEI (Figure 4), we found that the UHIEI varied from −0.4 to 0.4 with increasing distance in three periods, implying that there was spatiotemporal anisotropy of the UHI effect in Nanchang. Overall, the temporal and spatial anisotropy of the UHI effect was affected by various factors. When the UHIEI was less than 0, its area proportion was expanded along with increasing distance, indicating that the Ganjiang River could alleviate the UHI effect obviously (Figure 4(a)). Furthermore, the UHIEI was more than 0.5, and its value changed significantly as the distance increases (Figure 4(b,c)), suggesting that the main urban area with dense building, population density, and industrial layout caused the strengthening of the UHI effect. As the distance increased to 21 km, the UHIEI has been less than zero, implying that high altitude and vegetation relieved the UHI effect (Figure 4(d)).

Influencing factors identification
The regression parameters in year 2000, 2010, and 2019 established by OLS are shown in Table 4. The maximum value of the variance inflation factor (VIF max ) was less than 10, so there was no collinearity problem. So GWR was further applied to identify the influencing factors. It could be seen from Table 5 that the Adjusted R-squared were 0.811, 0.713, and 0.773, respectively, indicating that the GWR models can satisfy the accuracy requirement. From the spatial perspective, the local goodness-offit (Local R 2 ) in the center of the study area for the three periods was significantly lower than that in the northwest, southeast, and east directions ( Figure 5). The local goodness-of-fit (Local R 2 ) showed a decreasing first and then increasing later on the tendency. This implied that it gradually decreased from the suburbs to the city center. Its radiation shape invariably kept the maximum of three typical centers and shifted and changed over time. This indicated that the UHIEI was  mainly affected by the topography, vegetation, and bare land (soil) at first, and then by the urbanization process caused by the urban development and the built-up area expansion. Therefore, to effectively alleviate the future UHI effect of Nanchang City, we should scientifically and reasonably consider the layout of urban elements, such as expanding the construction of parks, strengthening the protection of nature reserves, and increasing the area of public green space.

Spatiotemporal heterogeneity of influencing factors
The regression coefficients for the various influencing factors in grid samples were obtained using the GWR models, as shown in Figures 6 and 7. These coefficients are applied to statistical calculations to obtain Mean, Maximum (Max), Minimum (Min), Standard Deviation (SD), and Median (Table 6). It can be seen from Table 6 that in 2000 and 2010, the overall NDSI regression coefficient was significantly larger than other variables, followed by NDBI and NDVI, while in 2019, the mean NDBI regression coefficient was significantly larger than other variables, followed by NDSI and NDVI. However, except for the NDBI is in 2019, the other coefficients contained more positive and negative values, such as DEM, NL, and highly correlated with the UHIEI, implying that influencing factors were higher spatial heterogeneity. As shown in Table 6 and Figure 6, compared to NDVI, the DEM and Slope were weaker correlated with the UHIEI. As for the spatial distribution of the DEM regression coefficient, the minimum value was mainly located in the Meiling mountains in the northwest of Nanchang City, gradually raising from the northwest to the southeast. For the spatial distribution of the slope regression coefficient, the maximum value of the slope was mainly distributed in the main urban area of the city, and the minimum value was located in the suburbs of the city. It decreased along both sides of the Ganjiang River to the suburbs over time. The spatial distribution of the NDVI regression coefficient changed significantly. In 2000 and 2010, the minimum value of NDVI was distributed in the northwest of Nanchang City with the higher forest cover area, while in 2019, the minimum value of NDVI was distributed in the South and northeast of Nanchang City with water body, suggesting that the vegetation can alleviate the UHI effect more than water body does. Table 6 and Figure 7 show that the regression coefficients of NDBI and NDSI were the largest among all variables. NDBI and NDSI had the greatest influence on the UHI effect when other variables were not considered. For NDBI, the largest area was mainly distributed in the northwest (Anyi county) and the southeast (Xihu District and Qinyunpu District). In contrast, the lowest area was mainly distributed in the southwest and the northeast. The spatial distribution of the NDSI regression coefficient was basically opposite to that of the NDBI. This is because most of the farmland was left fallow and exposed in November, leading to the intensification of UHI effect. As far as the spatial distribution was concerned, the spatial distribution of the minimum value of NL  regression coefficient was moving towards the central urban area and increased from the central urban area to the suburbs. It had little influence on the UHI effect.

Spatial regression model to identify influencing factors
It is imperative to choose an effective regression model to analyze the spatial variation of influencing factors. Traditional regressions are the primary tools for investigating the influencing factors of the UHI effect (Deilami, Kamruzzaman, and Liu 2018). However, most of the existing studies indicated that traditional regression model hides important details in the spatial variation and makes it impossible to capture the spatial dependence of the data (Guo et al. 2020). The assumption that evaluation coefficients without considering the spatial relationship between variables in the traditional regression model is unreasonable in this study because the UHI effect is proved to exhibit significant spatial heterogeneity (Zhao et al., 2018). Compared with traditional regressions, GWR can effectively reveal the spatial variation of the relationship between various factors and the UHI effect (Guo et al. 2020; Kashki  Zhi et al. 2020). Additionally, it further proved that the GWR model had the inherent potential to better understand the relationship between the phenomenon of the UHI effect and associated biophysical factors in China (Siqi and Yuhong 2020;Zhi et al. 2020).

Influencing factors of the UHI effect
The underlying surface of the city is a more complex ecosystem, which is produced by the joint action of nature, economy, and society. Researchers studying the UHI effect influencing factors have generally paid   attention to meteorological characteristics, landscape features, RS spectral information, land use types, urban morphology, topographic features, and socioeconomic characteristics (Kashki et al. 2021;Kurniati and Nitivattananon 2016;Zhi et al. 2020), and analyzed the impact of green space, waterbodies, and impervious surface on urban thermal environment (Bao et al. 2016;Du et al. 2016;Guha, Govil, and Besoya 2020;Weng, Lu, and Schubring 2004). The spatial characteristics of the UHI effect in the daytime, nighttime, seasonal, annual (Chen et al. 2021;Miles et al. 2017) and different spatial scales (Luo and Peng 2016) were also studied recently. Additionally, studies have begun to investigate the impact of socio-economic factors (Sabrin et al. 2020;Zhang, Estoque, and Murayama 2017) (e.g., air pollutants, gross domestic products (GDP), and population density) and two-dimensional and three-dimensional urban morphology (Qiao et al. 2020; Takebayashi and Senoo 2018) (e.g., building height, and building density) on the UHI effect. However, existing studies mostly focused single influencing factor and were lack of comprehensive analysis from the integration of natural and human perspectives (Deilami, Kamruzzaman, and Liu 2018;Dewan et al. 2021;Nwakaire et al. 2020;Ren et al. 2022;Xie et al. 2022;Zhao et al. 2021). The researches on the UHI effect are an inter-discipline topic which involves many disciplines such as geography, meteorology, urban planning, and geostatistics, but the traditional research cannot effectively integrate and make full use of the independent knowledge. In this work, GWR was adopted to solve the above mentioned problems. It can combine topography, vegetation, impervious surfaces, bare land (soil), and socioeconomic of the relative disciplines to reflect the UHI effect and determine the spatial distributions of influencing factors.
There are some limitations on this study needed to be further discussed and solved. First, the lack of threedimensional information can affect the results. In recent years, many studies have proved that the three-dimensional play a more important role than two-dimensional in shaping the urban thermal environment at different urban geometries (Berger et al. 2017;Sun et al. 2020). This is because the three-dimensional affect the surface energy balance and atmospheric transmission at local and urban scales, thereby changing the thermal environment inside the city, such as the intensification of the UHI effect (Alavipanah et al. 2018;Chen et al. 2022;Zheng et al. 2020). Therefore, three-dimensional information should be introduced in the following research to explore the impact on the UHI effect. Second, 420 m × 420 m spatial resolution was used for this study without considering the spatial scale effect. However, the scale effect could be an important factor affecting the precision of this study (Ren et al. 2022;Xie et al. 2022). Multiple scales are suggested to use for a better understanding of the UHI effect and influencing factors (Luo and Peng 2016;Miles et al. 2017). Third, the spatiotemporal distribution of the UHI effect in this study was obtained only in November. However, there are obvious seasonal and diurnal differences in the UHI effect (Chen et al. 2021). Therefore, the seasonal and diurnal spatiotemporal distributions and their driving factors should be further focused on for the study in the future.

Conclusion
Taking the midpoint of subway Line 1 with a radius of 25 km in Nanchang City as a case, the UHI effect of spatiotemporal characteristics in year 2000, 2010, and 2019 were analyzed. Furthermore, the influencing factors were explored by the GWR models. The main conclusions are as follows.
(1) The UHI effect of Nanchang City had increased from 2000 to 2019. Specifically, from 2000 to 2019, the "Heat Island" increased from 63.06% to 64.07%, while the "Green Island" decreased from 36.94% to 35.93%. The "Heat island" was mainly distributed in the main urban area with flat terrain and dense population, whereas the "Green Island" was concentrated in the northwestern region dominated by Meiling National Forest Park.
(2) In 2000, 2010, and 2019, the UHI effect of Nanchang City had significant temporal and spatial anisotropy. The UHIEI varied from −0.4 to 0.4 in the profile from Meiling National Forest Park to Poyang Lake, mainly affected by the topography, urban expansion, and land cover type. (3) The UHI effect was affected by many factors associated with underlying surface characteristics and had obvious heterogeneous characteristics. On the basis of identifying effective factors by OLS, the relationship between multi-temporal UHIEI and effective factors was analyzed by GWR. The spatial heterogeneity of regression coefficients of various factors can be better explained by the GWR model. The local regression coefficients of NDBI and NL show spatially significant differences, suggesting that there existed significant spatial non-stationarity. From the spatial distribution of Local-R 2 , it can be seen that with the expansion of the city, the fitting effect of the main urban area of the city shows a weakening trend as a whole. However, the overall radiation ranged showed an increasing trend, indicating that urban expansion played an important role in the UHI effect. Therefore, this study investigated the impact of urban underlying surface characteristics on the local thermal environment, and revealed the impact of influencing factors on the surrounding thermal environment by using the GWR. The results will help us better understand the influencing factors behind changes in the urban thermal environment, and provide support for future urban planning and construction to achieve sustainable urban development. Furthermore, in order to effectively improve the urban thermal environment, promotion of "zero carbon" buildings, high-density development of urban buildings control, reasonable urban planning for optimizing the layout of urban buildings and shaping a rational the urban green space are suggested to fulfil this task.

Disclosure statement
No potential conflict of interest was reported by the author(s).