Seasonal Contrast and Interactive Effects of Potential Drivers on Land Surface Temperature in the Sichuan Basin, China

: Little is known about the seasonal heterogeneity of land surface temperature (LST) and the interaction relationship between potential drivers in Sichuan Basin, China. In this study, based on exploring the spatial heterogeneity of LST in Sichuan Basin, China, multi-source remote sensing data as potential drivers were selected and a Geo-detector model was applied to analyze the main drivers and the interactive relationship between drivers on LST during different seasons. The results showed that the high-temperature areas in Sichuan Basin in different seasons all appeared in the cities near the high mountains on the edge of the basin. This phenomenon was summarized as “sinking heat island” by us. From the driving factors, the biophysical parameters (DEM, SLOPE and NDVI) had the greatest impact on LST in each season, reaching the peak in the transition season. The climate parameters (WIND, HUM, PRE and TEM) and socioeconomic parameters (LIGHT, POP and ROAD) also had a certain impact on LST. The inﬂuence of a single landscape parameter (SHDI, PD, LPI, ED and LSI) on LST is limited. From the effect of factor interaction on LST, the interaction of biophysical parameters, climatic parameters and landscape parameters from summer to the transitional season was strengthened obviously, and it showed a downward trend in the winter; in contrast, the socioeconomic parameters showed the opposite characteristics, indicating that the interaction between human activities and other factors affected LST more obviously in the winter. The results of this study are not only valuable for understanding the spatial features of LST but also important for formulating mitigation strategies and sustainable development of urban heat island in Sichuan Basin.


Introduction
The increase of population and economic activities worldwide have been promoted by urbanization and industrialization, and a series of environmental problems have been caused by the transformation of these activities on the natural environment, among which the urban heat island (UHI) has attracted much attention [1]. UHI generally refers to the phenomenon that the temperature in urban areas is higher than that in the surrounding suburbs or rural areas [2,3]. UHI has a profound impact on regional climate change (e.g., air pollution and extreme climate) [4,5], residents' health [6], labor productivity [7], ecosystem services [8] and energy consumption [9], which has attracted the attention of researchers in geography, urban planning, ecology and other fields [10,11]. Understanding the potential drivers of LST is essential to reduce the UHI effect, improve the quality of life of urban residents and promote regional sustainable development [12].
At present, the research foundation of quantitative UHI has been accumulated. Traditionally, urban thermal environments can be identified by comparing urban and suburban temperature differences using atmospheric temperature data observed by ground weather stations. However, the spatial representativeness of meteorological stations and the uncertainty of weather limit its application. In recent years, with the rapid development of Earth observation technology, multisource remote sensing data has been widely used in UHI research [13,14]. Remote sensing technology can provide effective urban environmental monitoring means at different spatial scales. Examples include Landsat Thematic Map-per/Enhanced Thematic Mapper+ (TM/ETM+), Moderate-resolution Imaging Spectroradiometer (Modis), Advanced Spaceborne Thermal Emission and radiation (Aster) and Advanced Very-High-Resolution Radiometer (AVHRR) satellite data [15][16][17]. The development of UHI research has been greatly promoted by the characteristics of high-resolution, multitemporal and multiscale of remote sensing data [18,19].
Although remote sensing data can help researchers depict the spatial and temporal patterns of LST, the purpose of exploring the law is to find solutions. The interaction between LST and influencing factors leads to the spatial heterogeneity of LST. Exploring this relationship can provide insights for reducing UHI phenomenon [20,21]. In general, the factors affecting LST mainly include four aspects: landscape indicators, socioeconomic indicators, climatic conditions and biophysical parameters [22]. Researchers have used a variety of mathematical analysis models to explore the relationship between UHI and spatiotemporal variables [23]. Gao et al. explored the spatial response of urban morphology factors to UHI using a geographically weighted regression model [20]. Based on a Pearson correlation analysis and principal component analysis (PCA), Sun et al. identified the contribution of various driving factors to Hangzhou UHI [24]. Taking Dalian as an example, Guo et al. used the spatial regression model (spatial lag model and spatial error model) to analyze the impact of the urban form on UHI at the community level [25]. These studies have provided valuable reference and basis for exploring the response of UHI to multiple driving factors and formulating mitigation strategies for UHI. However, the spatial characteristics of many influencing factors (e.g., evapotranspiration, vegetation and sky albedo) and their influence on LST changes with the seasons are ignored by many studies [26][27][28]. Analyzing the spatial characteristics of LST in different seasons and exploring the influence of multiple driving factors on LST are the basis for formulating mitigation measures of reducing the UHI effect. At the same time, most of the existing studies have mainly focused on the influence of a single factor and the determination of a nonlinear spatial relationship among multiple factors, but little is known about their interactions in different seasons [29,30].
The central topography of the basin is low, and the surrounding topography is high. The closed topography is not conducive to ventilation and heat dissipation, and the UHI effect is easy to form [31,32]. The Sichuan Basin is one of the most densely populated areas in China, with complex natural and social conditions, and some studies have focused on UHI of a typical city in the basin. For example, the average number of calm days in Chengdu is 240 days a year, and the UHI is becoming more intense during these years [33]. Chongqing, another big city in the Sichuan Basin, has the typical characteristics of a mountain city. The urban spatial form of Chongqing has a significant impact on the urban thermal environment and air quality [34,35]. However, the spatial pattern of different seasons of LST in the Sichuan Basin and the interactive relationship between influencing factors are still unclear. This study is dedicated to bridging this gap. In 2021, the outline of the construction plan for the Two-city Economic Circle in the Chengdu-chongqing region was issued by the Chinese government. The Sichuan Basin is required to take into account the multiple needs of economy, ecology, safety and health. Exploring the mitigation strategies of LST in the Sichuan Basin is of great significance to the regional high-quality and sustainable development.

Study Area
Sichuan Basin is located in the south central of the Asian continent (27 • 50 -33 • 03 N; 102 • 97 -110 • 11 E), the eastern part of the Qinghai-Tibet Plateau ( Figure 1). According to the classification data of climate zones in China, Sichuan Basin is mainly located in the middle subtropical zone and the north subtropical zone. The climate composition of Sichuan Basin is complex, and there is often heavy rainfall caused by a southwest vortex, valley wind caused by hilly topography and airflow disturbance in the Eastern Tibetan Plateau [23]. At the same time, Sichuan Basin is also one of the most populous regions in China and the world, and it is an important development pole in Western China. The Sichuan Basin includes 15 cities such as Chengdu and 29 districts and counties in Chongqing, with a total area of about 185,000 km 2 . According to the official data of the National Bureau of Statistics, the permanent population of Sichuan Basin urban agglomeration was 96 million, and the regional GDP was nearly $0.98 trillion in 2019. However, under the trend of urban development and construction in recent years, the phenomenon of extreme high temperature in the region has become increasingly prominent. According to the official data of the National Meteorological Administration, when the Sichuan Basin enters the third volt, the daily maximum temperature will reach 36-38 • C, of which the eastern and southern regions can reach more than 38 • C.

Study Area
Sichuan Basin is located in the south central of the Asian continent (27°50′-33°03′N; 102°97′-110°11′E), the eastern part of the Qinghai-Tibet Plateau( Figure 1). According to the classification data of climate zones in China, Sichuan Basin is mainly located in the middle subtropical zone and the north subtropical zone. The climate composition of Sichuan Basin is complex, and there is often heavy rainfall caused by a southwest vortex, valley wind caused by hilly topography and airflow disturbance in the Eastern Tibetan Plateau [23]. At the same time, Sichuan Basin is also one of the most populous regions in China and the world, and it is an important development pole in Western China. The Sichuan Basin includes 15 cities such as Chengdu and 29 districts and counties in Chongqing, with a total area of about 185,000 km 2 . According to the official data of the National Bureau of Statistics, the permanent population of Sichuan Basin urban agglomeration was 96 million, and the regional GDP was nearly $0.98 trillion in 2019. However, under the trend of urban development and construction in recent years, the phenomenon of extreme high temperature in the region has become increasingly prominent. According to the official data of the National Meteorological Administration, when the Sichuan Basin enters the third volt, the daily maximum temperature will reach 36-38 °C, of which the eastern and southern regions can reach more than 38 °C.

Data Source and Processing
Before collecting open access data, we reviewed all relevant service terms of the platform used to ensure that our studies were fully compliant with their agreements. As shown in Table 1, fifteen driving factors from four perspectives were selected: biophysical

Data Source and Processing
Before collecting open access data, we reviewed all relevant service terms of the platform used to ensure that our studies were fully compliant with their agreements. As shown  Table 1, fifteen driving factors from four perspectives were selected: biophysical parameters, climatic parameters, socioeconomic parameters and landscape parameters. The DEM, slope and NDVI are used to characterize the biophysical parameters [19,24]. The common climatic parameters include wind speed, humidity, precipitation and temperature [14,28]. Night light, population density and road network are commonly used to characterize the socioeconomic parameters [30,32]. The Shannon diversity index, patch density, largest patch index, edge density and landscape shape index are five commonly used landscape indicators as variables to explain landscape parameters [17,29]. LST: MODIS remote sensing images covering the Sichuan Basin in 2020 were selected for the study whose data are from LAADSDAAC (https://ladsweb.modaps.eosdis.nasa. gov, accessed on 23 December 2021). MODIS is an important sensor mounted on Terra and Aqua satellites. The spatial resolution of the data product is 250-1000 m, which can accurately reflect various information characteristics of the surface and atmosphere. Compared with other remote sensing images, MODIS adopts 8-day composite LST products, which can effectively reduce the impact of cloud cover compared to daily LST products. A large number of studies have shown the accuracy of MODIS LST data [36,37].
Sichuan Basin has obvious subtropical monsoon humid climate characteristics. Referring to Reference [38], the seasons of Sichuan Basin are divided into three types: summer, transition season and winter. As shown in Figure 2, the study area is mainly located in the middle subtropical zone, and very few regions are located in the north subtropic zone. Therefore, we need to consider the temperature of the two climate zones to divide the seasons. In general, the temperature and precipitation of the meteorological stations in the two climatic regions are similar. The average temperature from June to September is relatively high, and the distribution is relatively stable. The fluctuation range of the temperature in these four months is small, so these four months are designated as the summer [39]. Compared with the summer, the average temperature from December to March is lower, especially in February. Therefore, these four months are defined as the winter, while the average temperature from April to May and October to November are highly consistent, and they are often the transitional stage between summer and winter [40]. Therefore, these four months are defined as the transition season. Ⅱ mid-temperate zone, . south tempe Ⅲ rate zone, . Ⅳ plateau climate region, . north subtropical zone, . Ⅴ Ⅴ middle subtropical zone, .
Ⅺ south tropic zone and . Ⅻ the South China Sea islands).
Potential drivers: DEM data and NDVI data come from a geospatial data cloud platform (http://www.gscloud.cn/, accessed on 24 December 2021). Slope data is obtained from elevation data through ArcGIS 10.6 from the slope tool. Monthly average wind speed, monthly average humidity and monthly mean temperature data are from the National Earth System Science Data Center (http://www.geodata.cn/, accessed on 24 December 2021). The data of road density comes from the nonprofit map service platform open street map (http://www.openstreetmap.org/, accessed on 26 December 2021) through ArcGIS 10.6 preprocessing with the kernel density tool. The population density data were from Worldpop (https://www.worldpop.org/, accessed on 25 December 2021), and the night light data were from the world's first professional night light remote sensing satellite "Luojia-1" (http://59.175.109.173:8888/app/login.html, accessed on 25 December 2021). Land use data were from Globeland 30 (http://www.globallandcover.com, accessed on 26 December 2021). Globeland 30 is the first global geographic information public product provided by China to the United Nations. The results after processing with different parameters were shown in Figure 3. Potential drivers: DEM data and NDVI data come from a geospatial data cloud platform (http://www.gscloud.cn/, accessed on 24 December 2021). Slope data is obtained from elevation data through ArcGIS 10.6 from the slope tool. Monthly average wind speed, monthly average humidity and monthly mean temperature data are from the National Earth System Science Data Center

Landscape Pattern Index
The landscape pattern index can summarize the information of the landscape pattern [Error! Reference source not found.]. Five representative landscape pattern indices were selected by this study to reflect the characteristics of the landscape structure and spatial configuration( Table 2). The landscape pattern index was calculated by Fragstats4.2 software.

Landscape Pattern Index Calculation Formula Units
Ecological Meaning

Percent
It is helpful to determine the dominant types of landscape Reflect the shape complexity of the overall landscape

Standard Deviation Ellipse
A standard deviation ellipse is one of the classical methods to study the directional characteristics of the spatial distribution of geographical elements. The spatial distribution situation and evolution characteristics of the object are described quantitatively by a standard deviation ellipse with the basic parameters of the central point, azimuth angle, long semi-axis and short semi-axis, and the position change and moving direction of the center are identified; the calculation is as follows: where and are the latitude and longitude coordinates of the administrative center of the research unit, respectively, and is the value of the locale attribute of the study

Landscape Pattern Index
The landscape pattern index can summarize the information of the landscape pattern [38]. Five representative landscape pattern indices were selected by this study to reflect the characteristics of the landscape structure and spatial configuration ( Table 2). The landscape pattern index was calculated by Fragstats 4.2 software.

Standard Deviation Ellipse
A standard deviation ellipse is one of the classical methods to study the directional characteristics of the spatial distribution of geographical elements. The spatial distribution situation and evolution characteristics of the object are described quantitatively by a standard deviation ellipse with the basic parameters of the central point, azimuth angle, long semi-axis and short semi-axis, and the position change and moving direction of the center are identified; the calculation is as follows: where a i and b i are the latitude and longitude coordinates of the administrative center of the research unit, respectively, and X i is the value of the locale attribute of the study area, namely the heat island effect of the unit i. (A, B) is the study of the LST standard deviation elliptical central point coordinates.

Correlation Analysis
We began with the biophysical parameters, climatic parameters, socioeconomic parameters and landscape parameters influencing the distribution of the urban LST and used R language and ArcGIS 10.5 to analyze the correlation between them. We chose 100,000 points randomly throughout the sampling tool and then obtained their attribute values. The data were then imported into R language to carry out a Pearson correlation analysis and obtain the Pearson correlation coefficient between each index and the LST.

Geo-Detector
Geographic detector (Geo-detector) is a spatial analysis method based on a spatial analysis of variance. This method does not need a linear hypothesis, and for collinearity immunity, it can effectively identify the spatial differentiation of geospatial elements and reveal the driving factors behind it. In recent years, it has been applied to the research of UHI [22]. In this paper, LST in different seasons in the study area is taken as the dependent variable, and the factors affecting LST are independent variables. The "factor detector" and "interaction detector" in the geographic detector are selected to analyze the interpretation degree of each driving factor and factor combination on the spatial heterogeneity of LST in the study area.
(1) Factor detector module: Factor detection detects the impact of each influencing factor on LST spatial differentiation. The higher the Q value, the greater the impact on LST. The factor detection formula is as follows: where i represents the stratification of independent or dependent variables, N i represents the number of elements in layer i, N represents the number of elements in the whole domain, σ 2 i denotes the variance of dependent variables in layer i, σ 2 denotes the variance of dependent variables in the whole domain, SSW is the sum of the variances of the different layers and SST is the total variance of the whole region.
(2) Interactive detector module: Interaction detection is used to assess the degree of influence of different combination of driving factors on dependent variables. There are 5 types of interactions between two factors: if the q value of the two-factor interaction is less than any single-factor q value, (q(A) ∩ q(B) < MIN[q(A), q(B)]), it is nonlinearweakening. If the two-factor interaction q value between two single-factor q values is (MIN[q(A), q(B)] < q(A) ∩ q(B) < MAX[q(A), q(B)]), it is single-factor nonlinearweakening. If the q value of a two-factor interaction is greater than any single-factor q values, (q(A) ∩ q(B) > MAX[q(A), q(B)]), it is a two-factor enhancement. If the q value of the two-factor interaction is equal to the sum of two single-factor q values, (q(A) ∩ q(B) = q(A) + q(B)), it is independent. If the q value of a two-factor interaction is greater than the sum of two single-factor q values, (q(A) ∩ q(B) > q(A) + q()), it is a nonlinear enhancement.

Spatial Pattern of LST in Different Seasons
The surface temperatures of the Sichuan Basin in the summer, transition season and winter 2020 were obtained by supplementing MOD11A2 images of adjacent dates ( Figure 5). In the summer, represented by June, July, August and September, the monthly mean surface temperature is 29.25 • C, the maximum is 42.61 • C, the minimum is 3.67 • C and the standard deviation is 3.45 • . In the transition season represented by April, May, October and November, the monthly mean surface temperature is 19.82 • C, the maximum is 42.11 • C, the minimum is 2.01 • C and the standard deviation is 3.20 • . In the winter, represented by December, January, February and March, the monthly mean surface temperature is 9.99 • C, the maximum is 19.15 • C, the minimum is −15.29 • C and the standard deviation is 3.96 • .
Both the extreme hot and cold regions of LST in the Sichuan Basin throughout the year are located at the basin margins, with relatively stable temperatures at the center of the basin. Figure 6 and Table 3 show the direction of heat island development in the Sichuan Basin. The greater the oblateness is, the more directional the data. In 2020, the main axis of the spatial development of the heat island remained in the northwest-southeast direction, but the seasonal variations were different. The gravity center of the heat island shifted 67.64 km to the northwest from the summer to transitional season, and the deflection angle was 12.43 • . From the transitional season to the winter, the gravity center of the heat island shifted 16.97 km to the southeast, with a deflection angle of 2.08 • . On the whole, the evolution of the Sichuan Basin heat island is closely related to the spatial development pattern of Chengdu and Chongqing.

Correlation Analysis between LST and Influencing Factors in Different Seasons
In order to determine the relevance between LST and the biophysical parameters, climatic parameters, socioeconomic parameters and landscape parameters in different seasons of the Sichuan Basin, a Pearson correlation analysis was used to obtain the positive and negative relationships between the potential driving factors and LST in different seasons. The biophysical parameters and landscape parameters showed significant negative correlations with LST, indicating that the greater the above factors, the lower the corresponding LST. The results showed that the DEM, SLOPE, NDVI, SHDI, PD, LPI, ED and LSI have obvious effects on reducing LST, and all pass the significance test, which is consistent with the results of most previous studies. There is a significant positive correlation between the climatic and socioeconomic parameters and LST, indicating that WIND, HUM, TEM, LIGHT, POP and ROAD are effective in reducing LST; however, LST is influenced by precipitation, which showed a positive correlation in the summer and negative correlation in the winter and transition season (Figure 7).

Single Factor Attribution of LST in Different Seasons
Since the factors of the Geo-detector must be categories rather than values, 15 influencing factors such as biophysical parameters, climate parameters, socioeconomic parameters and landscape parameters are divided into five categories by the natural breakpoint method, so as to detect the main driving factors affecting the surface temperature of the Sichuan Basin. As shown in Figure 8, factor detection is used to analyze the impact of 15 driving factors on the surface temperature intensity in each season (q value). The results showed that the influencing factors passed the significance test. In the summer, the 15 potential drivers explained the intensity of the surface temperature from large to small: DEM (0. . Slope data had the greatest impact on the intensity of the surface temperature in the winter, followed by TEM and DEM. Compared with the landscape pattern parameters (q mean 0.0285), the biophysical parameters (q mean 0.6134) and climatic parameters (q mean 0.4572) had a greater impact on the intensity of the surface temperature in the winter.
(0.0256). The SLOPE data had the greatest impact on the surface temperature intensity the transition season, followed by TEM and DEM. The biophysical parameters (q mea 0.6581) and climatic parameters (q means 0.4935) had a great influence on the intensity the surface temperature in the transition season. In the winter, the explanations of 15 p tential driving factors for surface temperature intensity were: SLOPE (0.7760) > TE (0.7739) > DEM (0.6858) > HUM (0.5717) > POP (0.4108) > WIND (0.3911) > NDVI (0.37 > LIGHT (0.2842) > ROAD (0.2816) > PRE (0.0922) > ED (0.0378) > LSI (0.0353) > LPI (0.02 > SHDI (0.0227) > PD (0.0223). Slope data had the greatest impact on the intensity of surface temperature in the winter, followed by TEM and DEM. Compared with the lan scape pattern parameters (q mean 0.0285), the biophysical parameters (q mean 0.6134) a climatic parameters (q mean 0.4572) had a greater impact on the intensity of the surf temperature in the winter.

Interaction of Influencing Factors of LST in Different Seasons
The interaction between the driving factors of the surface temperature intensity in cities was tested by interactive detection. The results showed that the interaction q va of each pair of factors was greater than that of the two factors, and the interaction re

Interaction of Influencing Factors of LST in Different Seasons
The interaction between the driving factors of the surface temperature intensity in 15 cities was tested by interactive detection. The results showed that the interaction q value of each pair of factors was greater than that of the two factors, and the interaction relationship between the driving factors of the surface temperature intensity was synergistically enhanced. The interactive detection results of the summer surface temperature intensity showed that ( Figure 9) the maximum q (DEM ∩ ROAD) was 0.86, indicating that the interaction between DEM and ROAD had the greatest impact on the summer surface temperature intensity. Meanwhile, q (DEM ∩ NDVI) and q (LIGHT ∩ POP) both exceeded 0.85. The interactive detection results of the surface temperature intensity in the transition season showed that the maximum q (SLOPE ∩ TEM) was 0.91, indicating that the interaction between SLOPE and TEM had the greatest impact on the surface temperature intensity in the transition season, followed by q (DEM ∩ NDVI), q (SLOPE ∩ HUM), q (TEM ∩ SHDI) and q (TEM ∩ PD), which were all 0.89. According to the interactive detection results of the winter surface temperature intensity, the maximum values of q (SLOPE ∩ HUM) and q (TEM ∩ PD) were 0.85, indicating that the interactions between SLOPE and HUM and TEM and PD had the greatest impact on the winter surface temperature intensity, followed by q (TEM ∩ NDVI) and q (SLOPE ∩ SHDI), which were 0.84.

Spatial Pattern of LST
The study area is a large area with typical basin topography, including t

Spatial Pattern of LST
The study area is a large area with typical basin topography, including the Tibetan Plateau, Yunnan-Guizhou Plateau, Qin Mountains, Daba Mountains, Hengduan Mountains and other large mountains. The basin bottom can be roughly divided into three parts: parallel mountains and valleys in Eastern Sichuan, hills in Central Sichuan and the Chengdu Plain in Western Sichuan. The results showed that the heat island effect in the Sichuan Basin has obvious spatial differentiation. The relatively humid and cool microclimate environment formed is called cold island. The cold island mainly appears in the mountain area around the basin and the hilly area at the bottom of the basin. The heat island is mainly distributed in the parallel ridge valleys in Eastern Sichuan and the Chengdu Plain in Western Sichuan. In contrast to the overall spatial pattern of LST with "high in the middle and low periphery" in most regions [24,41], the LST in the Sichuan Basin presents the characteristics of "low in the outermost, high in the secondary periphery and low in the middle". From the LST of specific cities, the high-temperature areas are mainly concentrated in most areas of Chengdu and Chongqing; the southeast of Deyang and Mianyang; the north of Luzhou, Yibin and Leshan and the south of Dazhou and Nanchong. Here, there are dense populations, developed industry and agriculture, a large number of impervious water surfaces and cultivated land increasing rapidly, and the terrain of the basin does not easily dissipate heat, forming a large high-temperature area [13,21,42]. The low-value regions are mainly distributed in the central parts of the basin, such as Ziyang, Suining, Neijiang and Zigong, where the vegetation cover is dense; however, the area with flourishing vegetation usually has stronger evapotranspiration, which is characterized by a higher latent heat flux and lower sensible heat flux, thus forming a significant low-temperature area [43][44][45]. At the same time, the long axis and the short axis of the standard deviation ellipse are used to determine the spatial distribution trend of the high-temperature region, which reflects the spatial variation of the direction of the long axis in the northwest and southeast of the high-temperature region in the study area. In particular, from the summer to the transition season, the major axis of the ellipse decreases, the minor axis increases and the center of gravity shifts to the northwest, indicating that the east-west direction grows faster than the north-south direction. From the transitional season to the winter, the major axis and the minor axis of the ellipse increase, which shows that the high-temperature region of the outer edge of the ellipse grows faster than that of the central one.

Influencing Factors of LST
The current research has discussed the LST distribution pattern of plain cities. Generally speaking, the heat island effect of plain urban agglomeration is mainly affected by the socioeconomic environment [46][47][48]. The heat island effect basically shows a gradual weakening trend from the central urban area to the suburbs. However, this situation is not necessarily in a closed terrain (such as a basin or valley). Through the research, we found that the city with the hottest surface temperature in the Sichuan Basin is located at the edge of the Sichuan Basin, and the influence of the biophysical factors was greater than that of the other factors. This may be because the thermal air flow of the Sichuan Basin encounters the high mountains at the edge of the basin when it diffuses to the edge, which makes the thermal air flow sink, resulting in the temperature of the cities at the edge of the basin being higher than other regions; this "sinking heat island" has also been discussed in other studies. Most of the landscape parameters have negative correlations with LST, indicating that the landscape pattern can slow down the heat island, which has also been reflected in the research of others [17,40]. Most climatic parameters are positively correlated with LST, indicating that wind speed, humidity and temperature will enhance the heat island [46][47][48][49]. Precipitation shows significant seasonal variation characteristics on LST [50], showing a positive correlation in the summer and a negative correlation in the transition season and winter, which is different from the previous research that the total precipitation has a negative correlation with the annual LST [31,47]. It may be because the summer precipitation in the Sichuan Basin has a good negative correlation with the heating field of the Qinghai-Tibet Plateau, and the heating field of the Qinghai-Tibet Plateau increased in May. This leads to the decrease of precipitation in the south of the Sichuan Basin, which leads to a positive correlation between summer precipitation and LST in the Sichuan Basin [51]. The role of the socioeconomic activities has significantly enhanced the heat island, because urban expansion and population activities can indirectly increase the thermal energy intensity of LST by reducing vegetation and increasing road, building and energy consumption [52,53].
There are significant seasonal differences in the explanatory powers of the different driving factors on LST [33,54]. Specifically, in the hot summer, the ability of the socioeconomic parameters (LIGHT, POP and ROAD) to affect the LST intensity is stronger than that in the transition season and winter. Among them, POP with a q value of 0.4428 has a greater comprehensive influence than ROAD and LIGHT with a figure of almost 0.318. This shows that POP can fully explain the LST effect, and it can play a certain role in reducing UHI by adjusting the scale and structure of POP [55].
However, the ability of the biophysical parameters (DEM, SLOPE and NDVI) to affect LST is mostly stronger in the transition season, and the interpretation ability of the NDVI to LST is stronger in the winter, which is different from the previous research results, with the highest interpretation of the NDVI in the summer [17], which may be related to the unique vegetation coverage scale and structure formed by the terrain and climate environment of urban mountains in the Sichuan Basin agglomeration. Among them, the comprehensive impact of the SLOPE capability (q mean 0.7342) > comprehensive influence capability of DEM (q mean 0.7076) > comprehensive capability of NDVI (q mean 0.3478).
The comprehensive influence of the climate factor parameters on the LST intensity increases from the summer to transition season and then gradually decreases to the winter. Specifically, the changes of humidity and temperature are consistent with the comprehensive influence degree, but the influence ability of the WIND and PRE factors on LST decreases with the decrease of the seasonal temperature, comprehensive influence ability of TEM (q mean 0.7336) > comprehensive influence ability of HUM (q mean 0.5646) > comprehensive ability of WIND (q mean 0.4138) > comprehensive ability of PRE (q mean value 0.1359). It shows that the spatial variation of LST can be fully explained by TEM, HUM, WIND and PRE. The temperature is consistent with LST, which may be due to the long-wave radiation released from the surface that has an impact on the temperature [56]. High WIND can alleviate UHI by accelerating the temperature exchange [31,43,57]. Rainfall is usually accompanied by clouds blocking solar radiation, thus reducing LST, which has also shown in other studies that open rainfall can change soil moisture, which may affect plants and alleviate LST [48].
The overall trend of the influence of the landscape factor parameters on LST increased from the summer to transition season and decreased to the lowest level in the winter. As mentioned by previous scholars, the impact of the landscape factor parameters on LST was extremely limited [58][59][60], in which the LSI (q mean 0.0175) > ED (q mean 0.0170) > PD (q mean 0.0084) > LPI (q mean 0.0065) > SHDI (q mean 0.0064).
At the same time, the interaction effects of each driving factor on LST also have obvious seasonal differences. The biophysical parameters, climatic parameters and landscape parameters show an obvious interaction-strengthening trend from the summer to transition season, while, with the arrival of winter, they show an obvious downward trend, except the NDVI, LPI and other factors. The socioeconomic parameters show the opposite characteristics; the interaction from the summer to transition season shows a significant downward trend, and the interaction from the transition season to winter shows a significant strengthening trend. It shows that, with the decrease of the seasonal temperature, the impact of human comprehensive activities on LST is more and more significant. The above results are consistent with previous studies [30].

Research Limitations
This study still has some limitations and uncertainties. Firstly, we mainly used Geodetector to identify the spatial impact of different influencing factors on LST. Combined with other methods, an evaluation effect may be more reasonable, such as a questionnaire survey, random forest model, geographic weighted regression model and so on [61][62][63].
Secondly, typical representative potential driving factors from four aspects were mainly selected: biophysical parameters, climatic parameters, socioeconomic parameters and landscape parameters, but the factors affecting LST are varied in different fields. In addition to the possible driving factors considered in this study, other factors may also make non-negligible contributions to LST, such as haze pollution [51], sea-land breeze circulation [64], cloud amount [65], data spatial resolution [66] and other factors. Therefore, the interactions between other factors and LST need to be considered in order to reveal the driving mechanism more comprehensively.
Finally, there was uncertainty of the results due to the data collection time. The diurnal relationship between the valley wind and UHI should be studied on different scales. When the upward heat flow is blown to the topographic edge of the basin, it will be blocked by the mountains and lead to subsidence, so as to increase the heat island temperature in the topographic edge area and form the phenomenon of "sinking heat island". However, the interaction between the daytime "valley wind" and urban environment should be paid attention to, especially the day-night interactions between them. Since the peak can receive more heat during the day while the basin cannot, the temperature of the basin should be lower [67]. At night, however, the mountaintop radiates more heat, making the basin warmer [68].

Conclusions
Rapid urbanization has led to an obvious rise in LST all over the world, which has caused great harm to human and environmental health. How to effectively reduce LST has become an important proposition of urban ecology. The results show that the joint actions of different influencing factors form the spatial pattern of LST in different periods; the biophysical parameters have the greatest impact on the annual LST (the annual q mean is 0.5965), and the annual maximum is in the transition season, followed by the climate parameters (the annual q mean is 0.4620) and socioeconomic parameters (the annual q mean value is 0.3599), while the impact of a single landscape parameter on LST is weak, but interactions with other factors have a significant improvement effect. At the same time, there were significant seasonal differences in the interactions of the LST driving factors. The biophysical coefficients, climatic parameters and landscape parameters showed a significant strengthening effect from the summer to transition season, and the interaction increased with the arrival of winter. The interaction between human activities and other factors had a more significant impact on the winter LST. Therefore, rational use of the biophysical parameters and climatic parameters of the Sichuan Basin, especially DEM, SLOPE and TEM, and paying attention to the comprehensive impact of human activities and other factors in the winter are effective measures to alleviate UHI. The results of this study will help urban planners and policymakers make scientific and reasonable decisions and promote the future sustainable development of the Sichuan Basin. Data Availability Statement: All data for this paper are properly cited and referred to in the reference list.