Seasonal and diurnal surface urban heat islands in China: an investigation of driving factors with three-dimensional urban morphological parameters

ABSTRACT Surface urban heat island (SUHI) can considerably influence the urban environment and the quality of life. It is vital to examine how underlying surface properties impact seasonal and diurnal SUHIs. However, the influence of three-dimensional (3D) urban morphological parameters (UMPs) on SUHIs has not been thoroughly studied under varying climatic settings. To fill this knowledge gap, the present study investigated seasonal and diurnal changes in SUHI intensities (ΔT) in 208 cities in China from 2014 to 2016 using moderate-resolution imaging spectroradiometer (MODIS) land surface temperature products. In addition, the influence of potential factors in urban surface energy balance, including two-dimensional (2D) and 3D UMPs, socio-economic indices, urban greening, and surface albedo, on seasonal and diurnal ΔT were assessed under different climatic settings and with different city sizes using the method of Geographically Weighted Regression (GWR). Results show that negative summer daytime ΔT was observed in some cities under dry climates. Generally, in summer, the ΔT during daytime was higher than at nighttime. The 3D UMPs (i.e. building height and volume) yielded more decisive influences on ΔT than 2D UMP (i.e. building coverage). This is particularly true for the summer diurnal cycle and under dry climatic settings. Building height was found to be negatively correlated with surface temperatures, while building volume was positively correlated. Additionally, the 3D UMPs yielded more influences on winter ΔT than summer ΔT. The capability of vegetation to regulate ΔT was more potent in dry climates than in wet climates and in small cities than in large cities. Varying climates and city sizes can modify the significance of the 2D and 3D UMPs on the urban surface energy balance, suggesting that urban thermal mitigation should consider climate background and population size.


Introduction
Over the past few decades, rapid urbanization creates complex urban systems with large populations, commercial and industrial facilities, forming distinctive urban climates, and resulting in various environmental problems (Imhoff et al. 2010;Zhou et al. 2014;Ward et al. 2016). Urban heat island (UHI), a phenomenon that surface or air temperatures in urban areas are higher than those in rural settings, gets a lot of attention from various communities (Bonafoni, Baldinelli, and Verducci 2017;Manoli et al. 2019;Gawuc et al. 2020;Chakraborty, Sarangi, and Lee 2021). UHI effects can yield positive and negative influences on urban ecosystems according to different climatic backgrounds. On one hand, UHI can increase water demand and heat-related risks in some warm regions (Santamouris 2014;Wang, Guo, and Han 2021). On the other hand, it may favorably reduce energy demand in some cold regions (Roxon, Ulm, and Pellenq 2020;Yang, Yan, and Zhang 2020;Karimi et al. 2021). Therefore, it is crucial to investigate the impacts of UHI to help sustainable urban developments.
Generally, UHI can be divided into three different categories according to the layer characteristics of urban environments, i.e. surface urban heat island (SUHI), canopy urban heat island (CUHI) and boundary urban heat island (BUHI) (Meng et al. 2018;Hu et al. 2019;Gawuc et al. 2020;Wu et al. 2020;Yao et al. 2021a). A reliable investigation of BUHI is difficult due to sensor measurement uncertainty by the aircraft or tower (Venter, Chakraborty, and Lee 2021). CUHI, a UHI type closely related to human health, relies on meteorological in-situ observation, but it usually features insufficient spatial details (Venter, Chakraborty, and Lee 2021;Hu et al. 2019). SUHI has been widely utilized to examine UHI phenomena because of its wide spatial coverages, good time synchronization, and continuous improvement of land surface temperature (LST) qualities (Zhou et al. 2019;Sobrino and Irakulis 2020;Stewart et al. 2021). Given that the LST can be regarded as a critical indicator concerning surface radiation and energy exchanges, studies of the spatiotemporal patterns of SUHIs and their driving factors are valuable (Weng 2009;Wang et al. 2015;Chakraborty and Lee 2019;Yao et al. 2021b;Kelly Turner et al. 2022).
Previous studies have demonstrated that varying seasonal and diurnal cycles of SUHIs were observed (Mathew et al. 2018;Ma et al. 2021;Karimi et al. 2021;Chang et al. 2021;Mohammad et al. 2021). Zhou et al. (2014) investigated the seasonal and diurnal SUHI intensity (ΔT) in China's 32 cities using Moderateresolution Imaging Spectroradiometer (MODIS) land surface temperature (LST) data from 2003 to 2011. They found higher ΔT in summer than in winter during the daytime, but there was an opposite observation at night for the majority of cities. In addition, previous studies have demonstrated that SUHI variations are influenced mainly by climatic factors . Climate backgrounds are essential to explore SUHI dynamics because the crucial factors in terms of surface energy balance, including precipitations, latitudes, and background temperatures, often change over a large spatial scale. Those factors exert significant influences on vegetation activity (Kalisa et al. 2019), soil moisture (Manoli et al. 2020), and surface albedo . For example, Hu et al. (2019) analyzed the ΔT using MODIS LST data in China across different climates. A considerable difference of seasonal and diurnal ΔT was observed between dry and wet climatic settings and SUHI effects were more evident in dry climates. These studies also imply that the composition of dominant factors with respect to SUHIs tends to be different under different climatic settings (Manoli et al. 2019). Meanwhile, previous scholars mainly focused on the spatiotemporal variations of SUHIs in large cities. However, SUHIs' variations and their influencing factors in medium and small cities were not fully understood. Therefore, it is necessary to investigate the drivers of SUHI in different climates and city sizes to help the heat regulation for such cities that experience rapid urbanization and global climate changes.
Studies have been widely performed that expose the relationships between SUHIs and twodimensional (2D) urban morphological parameters (UMPs), e.g. vegetation coverage and building coverage. For example, Li, Zha, and Zhang (2020b) employed the biophysical and socio-economic factors to investigate the driving factors of SUHI using MODIS LST data across 419 global cities. They found that vegetation activity had a strong cooling effect on SUHIs in summer daytime in Asia, and surface albedos yielded a significant cooling effect on SUHI in winter nighttime in East Asia. Niu et al. (2021) explored spatial variations and driving factors of SUHIs in 281 cities in China for multiple years using MODIS LST data and the multiscale geographically weighted regression model. They found that normalized difference vegetation index (NDVI) and nighttime light (NTL) intensity exert significant influence on SUHI dynamics. These studies have thoroughly investigated the impact of 2D UMPs on SUHI at a large scale. However, few studies had examined the influence of threedimensional (3D) UMPs on SUHI dynamics due to the difficulty in obtaining 3D building data at a large scale. With the rapid development of photogrammetry and remote sensing technology, 3D UMPs (e.g. building height and volume) have become more accessible. Scholars gradually employed 3D UMPs to investigate the exchanges of surface latent and sensible heats and SUHI variations (Berger et al. 2017;Duan et al. 2019;Sun et al. 2020;Yang et al. 2021).
For instance, Huang and Wang (2019) studied the relationships between 2D/3D UMPs (e.g. orientation variance, mean height, and height variance) and LSTs in Wuhan city using high-resolution remote sensing data. They found that the influence of 3D urban morphology on LST is complicated and contextdependent. Li et al. (2021) examined the impact of 3D building forms (e.g. building density, height, and shadow) on seasonal LSTs by using the random forest regression method in Wuhan, China. They observed the heating and cooling effects of building density and shadow on LSTs. Chen et al. (2021) examined diurnal cycle differences in the impacts of 2D and 3D UMPs on urban LSTs in Nanjing using advanced spaceborne thermal emission and reflection radiometer (ASTER) data, and found that 2D UMPs yielded a more decisive influence on urban thermal environments than 3D UMPs during the daytime. However, a stronger effect of 3D UMPs on urban thermal environments than 2D UMPs was observed at night. Although many studies have investigated the effect of 3D UMPs on SUHI dynamics and explained the influence mechanism, how, until now, 3D UMPs impact SUHI dynamics at different climatic settings and city sizes was poorly understood.
Various regression models are widely used to explore the driving factors regarding SUHIs. A limitation of classic global regression models (i.e. Ordinary Least Square Regression model) is that they are unable to reveal the spatial heterogeneity of driving factors (Ivajnšič, Kaligarič, and Žiberna 2014;Zhao et al. 2018;Niu et al. 2021). While, at large scales, the interpretation of variables may yield significant spatial non-stationarity, local regression models could better address this non-stationarity of influencing factors (Zhao et al. 2018;Li, Zha, and Zhang 2020b;Niu et al. 2021). Zhao et al. (2018) investigated the influence of underlying biophysical attributes on the Surface Urban Heat Island (SUHI) phenomenon using the global regression and geographically weighted regression (GWR) in Austin and San Antonio, Texas. They found that GWR had a higher explanatory power of the influence regarding underlying factors than the global regression. Li, Zha, and Zhang (2020b) established the GWR model to assess the relationships between SUHI and driving factors with respect to surface energy balance in 419 global cities. They found that, compared with ordinary least square (OLS) and stepwise multiple linear regression (SMLR) models, higher determination coefficients (R 2 ) were found in GWR models than in others. Therefore, using local regression models to explore the spatial heterogeneity of driving factors could be a considerable method.
Based on the above discussion, this study chose 208 cities in China to explore the dominant factors of SUHI in these cities by using GWR models. Attempts were made to address the following questions: (1) How does seasonal and diurnal ΔT change under varying climatic settings and city sizes? (2) How do factors, especially 3D UMPs, control the urban surface energy balance and impact ΔT dynamics? And (3) Can climate backgrounds and city sizes regulate the composition of dominant factors, and if so, how? 2 Study area and data

Study area
We chose 208 cities in China for an analysis of dominant factors regarding seasonal and diurnal ΔT (Figure 1(a)). It is noteworthy that the analytical procedure excluded cities with urban area of less than 50 km 2 . Moreover, to preserve a comparable number of cities for different sizes and regions, cities in neighboring locations with others were removed from the analysis. Especially, the threshold value of 800 mm yr −1 precipitation was utilized to determine dry and wet climates for the following two reasons. Firstly, the threshold value of 800 mm·yr −1 precipitation is highly consistent with the Qinling Mountains-Huaihe River dividing line, distinguishing the northern and southern parts of China (Zheng 2008;Yang et al. 2017). Secondly, it is close to the 0°C isotherm line in January as well as the line that divides subtropical and temperate monsoon climates. Based on the division line of dry and wet climates, the cities selected were divided into 100 and 108 cities in dry and wet climates, respectively. Following Fang (2014), the cities selected could be further divided into large (population > 1 million), medium (500 thousand < population < 1 million), and small cities (population < 500 thousand) at prefecture-level. The population data was obtained from China Statistical Yearbook 2016 (NBSC (National Bureau of Statistics of China) 2016). The numbers of the large, medium, and small cities are 93, 72, and 43, respectively. It is hard to guarantee exactly equal population distribution among cities, i.e. few small cities can be selected at the prefecture-level.

Data sets
Here, we collected MODIS LST products from the National Aeronautics and Space Administration to calculate ΔT. Specifically, we used the MODIS/Aqua version 6 product (i.e. MYD11A2) and MODIS/Terra version 6 product (i.e. MOD11A2) during the period of 2014-2016. Daytime LSTs could be obtained by composing two observations during the daytime (i.e. local time at 1:30 P.M. and 10:30 A.M.). Similarly, nighttime LSTs were composed of two nighttime observations (i.e. local time at 1:30 A.M. and 10:30 P.M.) ( Figure S1). The data are 8-day composite LST product and feature fine temporal and spatial coverage (1-km spatial resolution), which is suitable for seasonal and diurnal ΔT exploration at a continental scale (Wan 2014). According to Peng et al. (2012), the LST data features high accuracy and a low root mean square error compared with in-situ observation. Table 1 shows the driving factors of seasonal ΔT in this study and their distribution can be found in Figure 2. Two-dimensional UMPs, i.e. building coverage (BC), and 3D UMPs, i.e. mean building height (MBH) and mean building volume (MBV), were acquired from (Li et al. 2020a). The data includes building footprint, height, and volume with a spatial resolution of 1-km in 2015. It was retrieved by using the random forest algorithm combining spatial data (e.g. Landsat 8 spectral bands and Sentinel-SAR images) and reference data (e.g. 3D buildings collection from commercial map provider and street view map). Compared with the observation data, the estimated R 2 values of the building footprint, height, and volume are 0.90, 0.81, and 0.88, respectively (Li et al. 2020a). In addition to UMPs, we also collected socioeconomic indices, vegetation, and surface albedo to explore SUHI dynamics. Population density (PD) and nighttime light (NTL) were utilized to represent socioeconomic influences. Normalized difference vegetation index (NDVI) and normalized difference water index (NDWI) were chosen to reflect vegetation activities and water content. We selected white sky albedo (WSA) covering shortwave spectrum to represent surface albedo (Peng et al. 2012). Mean annual precipitation was obtained from Global Precipitation Measurement (GPM, Huffman et al. 2019) and used for dry-wet climate division. All MODIS data were averaged over three years (2014)(2015)(2016) and resampled to a spatial resolution of 1-km to match the LST data. Figure 3 provides the workflow of the study. We tried to explore human-environment interactions in SUHIs and selected the crucial factors associated with urban surface energy balance as of the main drivers, e.g. socio-economic indices, 2D and 3D UMPs, vegetation, and surface albedos. In particular, regional climates and city sizes were chosen as regulating factors to reveal the consequences of human activities and climate settings on seasonal and diurnal SUHIs. As shown in Figure 3, the main steps include the generation of urban and rural boundaries, ΔT calculation, and statistical analyses. The following sections detail the main steps.

Generation of urban and rural boundaries
In this study, the impervious surface coverage of 25% was selected as the threshold to generate urban boundaries (data from Kuang et al. 2020). The reasons are as follows: first, the 25% impervious surface coverage provides a threshold to separate residential areas of medium to high intensity from those of low intensity based on a study in the US (Lu and Weng 2006); Secondly, the 25% impervious surface coverage closely matches the urban area of land use remote sensing monitoring data in 2015 in China from Resource and Environmental Science and Data Center (RESDC) (https://www.resdc.cn/data.aspx? DATAID=184) (Figures 1(b-d)). Rural regions were defined as buffer zones around urban regions (Fu and Weng 2018; Niu et al. 2020). Following Zhou  ( et al. (2015), we tested the impact of rural buffer areas on diurnal ΔT calculation. As shown in Figure S2, stable daytime and nighttime ΔT were observed when rural buffers were designed as 3 and 4 times areas of corresponding urban regions. Therefore, 3 and 4 times areas of urban regions were chosen as rural buffer areas during daytime and nighttime ΔT calculation. Additionally, water bodies, snow, and ice were removed from urban and rural regions to avoid their impacts on ΔT calculation. Moreover, following Imhoff et al. (2010), pixels with elevations exceeding ± 50 m of the median value of elevations in urban regions were also removed to avoid the impacts of elevations. Figures 1(b-d) provide urban and rural boundaries in Beijing, Kaifeng, and Jingdezhen.

SUHI calculation
We defined diurnal ΔT as the urban-rural differences in LSTs and explored the influences of critical factors regarding urban surface energy balances on summer (from June to August) and winter (from December to February) ΔT (Fu and Weng 2018). The three-year (2014-2016) average LSTs in summer and winter were calculated for both daytime and nighttime. The formula is as follows: where ΔT denotes surface urban heat island intensity and the ΔT in summer daytime, summer nighttime, winter daytime and winter nighttime were denoted as ΔT sd , ΔT sn , ΔT wd and ΔT wn , respectively; T u and T r are the average LSTs in urban and rural areas, respectively.

Statistical analyses
The urban-rural differences in driving factors, e.g. ΔBC and ΔMBH, were calculated with the same urban and buffer area as ΔT in each city and used as independent variables. In addition, the mean-subtraction method was used before analysis to eliminate the large differences between independent variables (Cai, Huang, and Song 2017): where i denotes the i-th city, x � i is the meansubtracted value of i-th city for each variable, μ is the average value of each variable, σ is the standard deviation of each variable.
We employed GWR models to explore the relationships between driving factors and ΔT. In particular, the SUHI intensities were used as the dependent variable and the independent variables were critical factors involved in surface energy exchanges ( Table 1). The GWR model, a local regression model, has the advantage over global regression (e.g. OLS model) in reflecting the spatial distribution characteristics of relationships between variables (Xu, Zhang, and Li 2019;Niu et al. 2021). We further used the GRW model to reveal SUHI variations: where y i , x ij , and ε are the dependent variable (i.e. ΔT), the j-th independent variable (i.e. driving factors), and the random error at the city i, respectively. m and k and indicate the number of city and independent variables, respectively. β j μ i ; # i ð Þ is the local regression parameter to be calculated at the location μ i ; # i ð Þ, which represents the longitude and latitude coordinates of the point i. β 0 μ i ; # i ð Þ is the geographically varying intercept. The weights of independent variables in GWR models were estimated by the function of the distance from point i. Generally, nearby locations will be assigned higher weights than distant ones (Dziauddin 2019). Thus, the weights of independent variables change in a specific space, namely, the spatial kernel. The extent of the spatial kernel is also called bandwidth. There are two methods to construct a spatial kernel: fixed distance and adaptive method. The fixed distance method uses a fixed bandwidth to solve each local regression. In contrast, the adaptive method selects a specified number of neighbors to construct a spatial kernel. The denser the spatial distribution of regression location (i.e. the selected cities) is, the fewer the neighbor is. Here, we used the fixed kernel type to build GWR models due to the even and dense distribution of selected cities. To diagnose the performances of GWR models, we determined the coefficient of determination (R 2 ), Akaike Information Criterion (AICc) and standardized residual of models. In addition, the spatial autocorrelation of residuals was tested by Moran's Index (MI).

Seasonal and diurnal changes in SUHIs
Figure 4(a) shows that positive ΔT was observed in most cities during summer daytime (> 95%). About 30% of the cities selected yielded large ΔT sd (higher than 3°C). In addition, some cities associated with negative ΔT sd were observed in China's northwestern part, such as Jiayuguan (−2.29°C) and Hami (−0.57°C). Figure 4(b) shows that most cities had positive ΔT during summer nighttime (> 99%). Medium ΔT sn (1 < ΔT sn < 2°C) was observed in 67% of cities. Only 17 cities, mainly in dry climates, yielded large ΔT sn . It is noteworthy that, in summer, the daytime ΔT was higher than that at night (Figures 4(a,b), which may be due to the underlying surfaces absorbing solar radiations and producing more anthropogenic heat (e.g. air conditioning cooling) during the daytime.
In winter, about 39% of cities yielded negative daytime ΔT. Those cities are mainly located in dry climates, such as Lhasa (−1.90°C), Tongliao (−1.75°C), and Panjin (−0.57°C). However, in wet climates, most cities produced positive ΔT wd (Figure 4(c)). At night, large ΔT wn was found in dry climates. Additionally, in wet climates, a similar ΔT distribution at night was found compared with that in daytime (Figure 4(d)). Interestingly, nighttime ΔT was observed higher than daytime ΔT in winter (especially in dry climates), and this differs in summer diurnal cycles (Figures 4 (c,d)). In winter nighttime, cities in dry climates would produce more anthropogenic heat (i.e. heating in winter) than those in wet climates. Figures 4(e,f) show the average ΔT under varying city sizes in dry and wet climates. First, the city size can significantly impact the average ΔT in seasonal and diurnal cycles. Usually, a larger city size produces a higher ΔT. Notably, average summer daytime ΔT in wet climates was higher than that in dry climates. The possible reason is that, in dry climates, rural areas often have sparse vegetation, leading to limited heat mitigation of vegetation. In contrast, high-rise buildings in urban areas can enhance the efficiency of land-air convection, e.g. increasing shading of ground surfaces and accelerating air convection. All of these lead to low summer daytime SUHIs in dry climates. Meanwhile, average nighttime ΔT in wet climates was lower than that in dry climates, which could be explained by the warming effect of dense vegetation (compared with soils in dry climates) on rural LSTs in wet climates. It indicates that the influence of driving factors on ΔT may vary by different climates and by varying seasonal and diurnal cycles.

Diagnostics of GWR models
We inspected the correlations between the chosen variables using the variance inflation factors (VIFs, revealed by the OLS model) and correlation matrix to determine the multicollinearity level in our models. The VIFs are displayed in Table S1, which shows all VIFs of variables were below 10, indicating that no redundant variables were observed in the regression. Table S2 shows the correlation matrix concerning variables involved in analyses. As shown, most of the correlation values were below 0.5, indicating the chosen variables did not suffer multicollinearity. Table 2 shows the diagnostics of GWR models. As shown, seasonal and diurnal variations in adjusted determination coefficient (R 2 ) were observed, indicating that influencing factors used in the study had different total variance explained in seasonal and diurnal cycles. Overall, low Akaike information criterion (AICc) and root mean squared error (RMSE) values in GWR models were observed. It suggests the excellent performance of the models selected in revealing ΔT spatial variations, especially in summer daytime (AICc: 324.23; RMSE: 0.43) and winter nighttime (AICc: 358.46; RMSE: 0.43). In addition, the Moran's I (MI) of standardized residuals was close to 0 in GWR models (p-value < 0.05), suggesting that the residuals feature no evident spatial autocorrelation. Table 3 shows the mean absolute values of the T-statistic tests regarding driving factors using GWR models. As shown, the mean absolute values of the T-statistic test concerning ΔBC, ΔMBH, ΔMBV, and ΔWSA were higher than the threshold value of 1.96 in different periods, indicating these factors were statistically significant at the level of 0.05 for seasonal and diurnal ΔT. Notably, ΔNDVI had the highest mean absolute value of the T-statistic test in summer daytime (4.52), indicating ΔNDVI can significantly influence the spatial variations in ΔT sd . In addition, ΔPD was statistically significant for ΔT sd and ΔT wn (i.e. the mean absolute value > 1.96), ΔNTL was statistically significant for ΔT wn and ΔNDWI was statistically significant for ΔT sd and ΔT sn . Figure 5 shows the local R 2 values of each city in GWR models under varying seasonal and diurnal settings. As shown, the performance of the GWR model in summer daytime (the local R 2 values changed from 0.4 to 0.9) was better than that in other periods, indicating the driving factors used can better explain ΔT sd spatial variations ( Figure 5(a)). Low R 2 values were observed in China's central parts in summer nighttime, suggesting more influencing factors should be considered for ΔT sn explanations in these regions ( Figure 5(b)). Similarly, in winter daytime, Beijing, Tianjin, and Tangshan regions were associated with low local R 2 values ( Figure 5(c)). Notably, significant differences in local R 2 values of ΔT wn were observed between dry and wet climates, indicating that, in winter nighttime, climatic division can influence the performance of crucial factors regarding urban surface energy balance to control SUHI intensities ( Figure 5(d)). Figure 6 shows the spatial distribution of the regression coefficient (β) values concerning driving factors in summer daytime. As shown, ΔPD was positively correlated with ΔT sd , and a high β value in dry climates (β value = 0.27 ± 0.08, Table 4) was observed. However, it was found that the ΔPD influence was non-significant (p-value > 0.05) in most cities in wet climatic settings (Figure 6(a)). Additionally, slightly positive correlations between ΔNTL and ΔT sd were observed in 26.92% of municipalities, mainly located in eastern China with rapid economic developments and dense human activities (Figure 6(b)). Cities in this region generate more anthropogenic heats, leading to the rising ΔT sd . These findings suggest high human activity intensities can increase summer daytime ΔT, but the increase varied by regional climates and spatial variations.

Relationships between seasonal and diurnal ΔT and driving factors
Favorably positive links between ΔBC and ΔT sd were observed in China's northeastern parts (the highest β value = 0.40, Figure 6(c)), suggesting the warming effects of building coverages in summer daytime. As shown in Figure 6(d), ΔMBH was highly negatively correlated with ΔT sd , and more substantial influences of ΔMBH in dry climates (β = −0.37 ± 0.12, Table 4) was observed than those in wet climates (β = −0.18 ± 0.11, Table 4). The potential reason is that buildings cast shadows and accelerate air convection can reduce surface temperatures in urban areas. It might be more pronounced in dry climates because of the sparse vegetation and low solar altitude angle. As shown in Figure 6(e), positive links between ΔMBH and ΔT sd were observed and the links were stronger in dry climates (β = 0.27 ± 0.16, Table 4) than those in wet climates (β = 0.13 ± 0.12, Table 4).
ΔNDWI was slightly negatively correlated with ΔT sd (β = −0.13 ± 0.06, Figure 6(f) and Table 4), suggesting vegetation moisture contents can mitigate surface warming to some extents. Additionally, a strong negative correlation between ΔNDVI and ΔT sd was observed, especially in dry environments (minimum β value = −0.99) (Figure 6(g)), indicating vegetation coverage is a crucial cooling factor in summer daytime. It was found ΔWSA was positively correlated with ΔT sd (Figure 6(h)). Generally, higher urban-rural differences in albedo should reduce the energy absorbed by urban areas, leading to the reduction of SUHIs. However, low vegetation covers can increase albedo to some extent since vegetation components are "darker" than most built-up structures. This may lead to a positive correlation between ΔWSA and ΔT sd .   (Figure 7(a)). Additionally, ΔNTL was positively correlated with ΔT sn in China's southeastern parts, which was consistent with that in daytime (Figure 7(b)). These findings suggest human activities yielded a similar impact on ΔT in summer diurnal cycles.
It is noteworthy that slight positive correlations between ΔNDWI and ΔT sn were observed in dry climates (β = 0.10 ± 0.16, Table 4), which differs in summer daytime (β = −0.17 ± 0.08, Table 4) (Figure 7(f)). It suggests that distinguishable roles of ΔNDWI are played in surface energy balance during summer diurnal cycles. Generally, ΔNDVI significantly mitigated ΔT during summer diurnal cycles (Figure 7 (g)); however, the cooling effects of ΔNDVI on underlying daytime surfaces (β = −0.42 ± 0.10, Table 4) were observed higher than those on nighttime surfaces (β = −0.27 ± 0.20, Table 4). Notably, ΔWSA yielded strongly negative influences on ΔT sn in the cities selected (β = −0.24 ± 0.17), suggesting a crucial cooling effect of ΔWSA on ΔT sn (Figure 7(h) and Table 4). The possible reason is that a lower albedo leads to higher ground heat storage, and this heat is released at night, increasing the nighttime SUHI.  Figure 8 shows the β spatial distribution of different influencing factors in winter daytime. As shown in Figure 8(a), ΔPD was highly correlated with ΔT wd (β = 0.20 ± 0.11, Table 4). Additionally, the warming effects of ΔNTL in winter daytime (β = 0.20 ± 0.07, Table 4) were higher than those in summer daytime (β = 0.09 ± 0.08, Table 4) (Figure 8(b)). These findings suggest that, as typical heating factors, human activities in winter can exert more decisive influences on ΔT than those in summer. As to 2D UMP, the positive correlation between ΔBC and ΔT wd in dry climatic settings (the highest β value = 0.25) was observed more substantial than that in wet environments (β = 0.13 ± 0.04, Figure 8(c) and Table 4). In addition, the correlation in winter daytime (β = 0.15 ± 0.05, Table 4) was found lower than that in summer daytime (β = 0.28 ± 0.04, Table 4), suggesting that building coverage yields more warming effects on Earth's skins in summer than those in winter. As to 3D UMPs, negative correlations between ΔMBH and ΔT wd in high-latitude regions (i.e. dry climates, β = −0.45 ± 0.11) were observed higher than those in low-latitude regions (i.e. wet climates, β = −0.14 ± 0.09) (Figure 8(d)). The possible reason is that cities in high-latitude regions can cast large building shadows and thus reduce surface temperatures in urban areas. As shown in Figure 8(e), ΔMBV was highly positively correlated with ΔT wd , and the correlation in dry climates was found higher than that in wet climates (Table 4). It was found that ΔNDWI was negatively correlated with ΔT wd (β = −0.14 ± 0.12, Figure 8(f) and Table 4). Additionally, it was highly negatively correlated with ΔT wd in the northeastern parts of China (the minimum β value = −0.46), which could partly be explained by the vegetation created by special water and energy conditions in high latitudes (Sun et al. 2013). The negative correlation between ΔNDVI and ΔT wd in winter daytime (β = −0.14 ± 0.12) was lower than in summer daytime (β = −0.42 ± 0.10, Figure 8(g) and Table 4). Notably, ΔWSA showed a weak negative correlation with ΔT wd in the cities selected (β = −0.09 ± 0.20), which differs from that in summer daytime (Figure 8(h) and Table 4). Figure 9 provides the β spatial distribution of driving factors revealed by GWR models in winter nighttime. As shown in Figure 9(a), a weakly positive correlation between ΔPD and ΔT wn was observed (β = 0.16 ± 0.10, Table 4), and the correlation in winter nighttime was lower than that in other periods, e.g. summer nighttime. It was found that highly positive correlations between ΔNTL and ΔT wn in coastal cities in China's northern parts (the highest β value = 0.36, Figure 9(b)), which differs from those in winter daytime.
ΔBC showed slight positive correlations with ΔT wn (β = 0.12 ± 0.11, Table 4 and Figure 9(c)), lower than those in other periods. It suggests that the warming effects of building coverage are not evident compared with those in different periods, e.g. winter daytime. As to 3D UMPs, the links between 3D UMPs and ΔT wn were similar to those in daytime, i.e. 3D UMPs in dry climates yielded higher impacts on ΔT than those in wet climates (Figures 9(d,e)). It indicates that the influence of 3D UMPs on surface energy balances exist in diurnal inertia to some extent. Notably, positive correlations between ΔNDWI and ΔT wn were observed in dry climates (β = 0.22 ± 0.10, Figures 9(f) and Table 4), suggesting a warming effect of vegetation moisture on ΔT wn . Still, ΔNDVI and ΔWSA were negatively correlated with ΔT wn (Figures 9(g,h)), presenting their cooling effects in winter nighttime.

Influence of UMPs on seasonal and diurnal SUHIs: 2D vs. 3D
In this study, it was found that 3D UMPs (i.e. ΔMBH and ΔMBV) yielded more decisive influences on ΔT than 2D UMP (i.e. ΔBC); this is particularly true in summer daynight cycles and under a dry climatic setting. Our finding is in line with an investigation of the influence of urban structures (including 2D and 3D) on LSTs using remotely sensed observations by Ezimand, Azadbakht, and Aghighi (2021). Thus, it is essential to improve 3D UMPs to perform SUHI mitigations (He et al. 2021). Additionally, ΔBC and ΔMBV exerted significantly warming effects on Earth's skins. The possible reason is that buildings featured as impervious surfaces can absorb shortwave solar energy during the daytime and release it into atmospheres during the nighttime, leading to warming effects on underlying surfaces during diurnal cycles (Manoli et al. 2019;Hong et al. 2020;Morabito et al. 2021;Taurines et al. 2021). Meanwhile, building energy consumption increases anthropogenic heat emissions, resulting in SUHI rising (Santamouris 2014;Yu et al. 2021;Yang and Li 2013;Ye et al. 2021).
Seasonal variations in the impacts of 2D and 3D UMPs on ΔT were observed. The warming effects of ΔBC in summer were more substantial than those in winter. The possible reason is that buildings absorb higher solar radiation and create fewer building shadows in summer than in winter (Khamchiangta and Dhakal, 2019). However, the cooling effect of ΔMBH in winter was higher than in summer, partly because of the lower solar elevation angle in winter (Huang and Yang, 2019;He et al. 2021;Li et al. 2021). Beyond this, vegetation evapotranspiration becomes weak in winter; thus, the cooling effects of the land-air convection caused by building roughness determine the SUHI mitigation in winter (Berger et al. 2017;Niu et al. 2021). For the diurnal cycle of SUHIs, we found that the warming effect of 2D UMP on ΔT exhibits consistency during day and night observations, which is consistent with Zhou et al. (2014). However, significantly various performances of ΔMBH and ΔMBV on ΔT were found during diurnal cycles. ΔMBH can mitigate ΔT by casting shadows, which has been proved by Li et al. (2021). Therefore, the giant building shadows on surfaces can explain the more substantial cooling effect of the ΔMBH in daytime. In addition, the notable warming effect of ΔMBV on nighttime ΔT indicates that buildings would release more anthropogenic heat at night, which could generate heat upward to the atmosphere and then contribute to surface warming (Hu et al. 2019;Hong et al. 2020;Yao et al. 2021b).

The modification of climatic factors
Using GWR models, we revealed the spatial variations in the influences of dominant factors regarding ΔT under dry-wet climatic settings. It was found that ΔBC yielded a low warming effect on winter ΔT in wet environments (Figures 9(c)). The possible reason is that, in winter, building areas feature lower heat storage abilities than their rural settings with wet soil environments He et al. 2021). Additionally, the cooling effects of ΔMBH in the dry climate were higher than those in the wet climate ( Figures 6-Figures 9(d)). Previous studies have demonstrated that tall buildings could increase the land-air convection efficiency to help SUHI mitigations (Berger et al. 2017;Huang and Wang 2019;Yu et al. 2020;Budhiraja, Gawuc, and Agrawal 2019;Sun et al. 2020). In particular, in the dry climate, the cooling effects of low and sparse vegetation are limited. Thus, building roughness plays a critical role in surface cooling by increasing land-air convection efficiency. We also found that ΔMBV yielded a higher warming effect on winter ΔT in dry climates than that in wet climates (Figures 8(e) and 9(e)). It is because high latitude cities in China during heating periods could produce more anthropogenic heat emissions and thus increase winter SUHI. The cooling effect of ΔNDVI in dry climates (β = −0.60 ± 0.08) was higher than that in wet climates (β = −0.25 ± 0.05) in summer daytime (Figure 6(h) and Table 4). This finding is in line with (Manoli et al. 2019). The possible reason is that the cooling effect of vegetation evapotranspiration is more considerable for regions with low average precipitations and soil water contents (Buyantuyev and Wu 2012;Imhoff et al. 2010;Wang et al. 2015). Additionally, ΔWSA in summer daytime yielded higher warming effects on ΔT in wet climates than in dry climates (Figures 6h). The reasons can be attributed to the fact that a high WSA is often associated with higher building coverage, contributing to rising surface temperature (Bonafoni, Baldinelli, and Verducci 2017;Zhou et al. 2019;Sekertekin et al., 2021). Cities in wet climates are often surrounded by evergreen forests and croplands, which enlarge ΔWSA and thus significantly increase ΔT. Moreover, ΔWSA yielded cooling effects on winter ΔT in dry climates. It can be explained by the high brightness of bare soils featuring high albedo in the dry climate, which could help surface cooling (Haashemi et al. 2016).

Influence of city size on SUHIs
Generally speaking, in dry climates, ΔT in large and medium cities was observed higher than in small cities (Figures 4(e,f)). Large and medium cities produce more anthropogenic heat than small cities; additionally, cities in dry climates feature high thermal inertia, increasing surface temperatures by storing latent heat (Clinton and Gong 2013). In wet climates, small cities are usually surrounded by tall and dense vegetation, leading to a high evaporation rate and LST mitigations in rural areas (Ward et al. 2016). Meanwhile, building heights in small cities are often lower than those in large and medium cities, resulting in low land-air convection efficiencies in urban areas (Berger et al. 2017;Manoli et al. 2019).
In addition, the regulation of city sizes on ΔT was tested in this study. Table 5 provides the average and standard deviation of β values in different city sizes. As shown in Table 5, the influences of 2D UMP (e.g. ΔBC), 3D UMPs (e.g. ΔMBH and ΔMBV), and human activities (e.g. ΔPD) on seasonal and diurnal ΔT in large cities were more substantial than those in medium and small cities.
Additionally, it was found that ΔNDVI can help SUHI mitigation more in small cities. This finding is important since it suggests that the regulation of UMPs and population in large cities yields a better effect than in small cities (Yao et al. 2021a). In contrast, urban greening could help surface cooling in small cities more than in large cities.

Limitations and future research directions
We noted several limitations in this study. Firstly, we investigated the influences of 2D and 3D UMPs on seasonal and diurnal ΔT at a 1-km spatial scale to match the resolution of MODIS LST and other remotely sensed products that were used as driving parameters. Our study has considerable limitations in revealing 2D and 3D urban forms at the city-block scale, leaving uncertain driving analysis results (Bhatta, Saraswati, and Bandyopadhyay 2010). For instance, ΔMBH and its impacts on ΔT may underestimate in large cities due to the coarse resolution and averaging process in urban areas. Therefore, studies exploring the influences of 2D and 3D UMPs on ΔT with more satisfactory spatial resolution over a large scale are required.
Secondly, several advantages of GWR models are worthy of declaring and can be summarized as follows. Firstly, GWR models could better reveal the spatial heterogeneity of driving factors than global regression models. Secondly, GWR models can reflect the distinctive roles of driving factors under varying climate settings and across different sizes (Tables 4  and 5), e.g. we found that ΔMBH yielded a higher cooling effect on ΔT sd in dry climates than in wet climates. However, the R 2 values of GWR models were lower than 0.5 in several cities ( Figure 5), implying that current variables cannot fully explain ΔT in some regions. Factors, including landscape pattern, wind speed, solar radiation, and analyses at different spatial scales (e.g. city blocks), should be investigated in future studies (He et al. 2021;Li et al. 2021). Notably, in related research, the analytic results regarding the influence of critical factors on SUHIs should be carefully verified. Otherwise, we may reach misleading conclusions. For example, our results showed a positive correlation between ΔWSA and summer daytime SUHI, which makes no sense since higher urban-rural differences in albedo should reduce daytime SUHI by decreasing underlying surfaces' energy absorbed (Chakraborty, Sarangi, and Lee 2021). The possible reason for this statistical bias is probably the relation between albedo and vegetation.
Thirdly, the coarse division of urban-rural boundaries may lead to some uncertainty due to the heterogeneous surfaces in urban and rural areas (Stewart and Oke 2012). Combining the surface structure, cover, and human activities with local climate zones (LCZs) can provide a standard framework for reporting and comparing air/surface temperature observations .

Conclusions
In this study, the spatiotemporal patterns of seasonal and diurnal ΔT over 208 cities in China were investigated under varying climates and city sizes. We explored how key factors controlling the urban surface energy balance, i.e. socio-economic indices, 2D and 3D UMPs, vegetation, and surface albedos, impact seasonal and diurnal ΔT from 2014 to 2016. Regional climates and city sizes were further employed to reveal the consequences of human activities and climate settings on seasonal and diurnal cycles of ΔT. The main conclusions can be drawn as follows. • The 3D UMPs (i.e. ΔMBH and ΔMBV) yielded more decisive influences on ΔT than the 2D UMP (i.e. ΔBC); this is particularly true in the dry climate. The cooling effects of ΔMBH on ΔT in dry climates were higher than in wet climates. The warming effects of ΔMBV on ΔT wn in dry climates (β = 0.40 ± 0.12) were significantly higher than those in wet climates (β = 0.20 ± 0.08). • The cooling effect of ΔNDVI on ΔT sn in dry climates (β = −0.32 ± 0.15) was higher than that in wet climates (β = −0.22 ± 0.10). ΔWSA was positively (negatively) correlated with ΔT sd (ΔT sn ). The socio-economic indices, i.e. ΔPD and ΔNTL, had warming effects on seasonal and diurnal ΔT. Population density was an important heating factor in the northern part of China. In contrast, night light intensity mainly increased ΔT in the southeastern part of China. • Varying control patterns of dominant factors were observed in different climates and city sizes. The regulation of UMPs to control urban thermal environments in large cities yielded a better effect than in small cities. In addition, urban greening could help surface cooling in small cities more than in large cities (in dry climates more than in wet climates).
The study contributed to UHI studies in the following aspects. Firstly, the study provides a clear picture of the influence of various driving factors, especially 3D UMPs, on seasonal and diurnal ΔT at a large scale. Secondly, how human activities and climatic settings have interacted to modify the influence of these factors on the urban surface energy balance was determined. The finding suggests that urban thermal regulations should consider climate backgrounds and population sizes. Finally, we employed GWR models to reveal the varying influences of driving factors under different climatic settings and city sizes.

Highlights (for review)
• Driving factors of seasonal and diurnal surface urban heat islands were investigated under varying climates by geographically weighted regression models • Urban-rural difference in building height and volume yielded more decisive influence on surface urban heat islands than that in building coverage