Next Article in Journal
Flexible-Clustering Based on Application Priority to Improve IoMT Efficiency and Dependability
Previous Article in Journal
Evaluation of Soil Quality in a Composite Pecan Orchard Agroforestry System Based on the Smallest Data Set
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Exploring the Variations and Influencing Factors of Land Surface Temperature in the Urban Agglomeration on the Northern Slope of the Tianshan Mountains

1
School of Geography and Tourism, Xinjiang Normal University, Urumqi 830054, China
2
Research Centre for Urban Development of Silk Road Economic Belt, Xinjiang Normal University, Urumqi 830054, China
3
Xinjiang Key Laboratory of Lake Environment and Resources in Arid Zone, Urumqi 830054, China
4
Key Laboratory of Regional Sustainable Development Modeling, Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences, Beijing 100101, China
*
Author to whom correspondence should be addressed.
Sustainability 2022, 14(17), 10663; https://0-doi-org.brum.beds.ac.uk/10.3390/su141710663
Submission received: 20 July 2022 / Revised: 22 August 2022 / Accepted: 24 August 2022 / Published: 26 August 2022

Abstract

:
Changes in land surface temperature (LST) can have serious impacts on the water cycle and ecological environment evolution, which in turn threaten the sustainability of ecosystems. The urban agglomeration on the northern slopes of the Tianshan Mountains (UANSTM) is located in the arid and semi-arid regions of northwest China, with an extremely fragile ecological environment and sensitive to climate change. However, studies on the LST of the UANSTM have not received much attention. Therefore, this study explored the spatial distribution pattern, fluctuation characteristics, and influencing factors of the LST of the UANSTM from 2005 to 2021 based on MODIS time series LST data and the geo-detector model with optimal parameters. The results show that the UANSTM is dominated by medium- and high-temperature classes, with high- and extremely high-temperature classes clustered in Turpan City. The daytime and nighttime LST patterns are significantly different, with a typical “daytime cold island and nighttime heat island” feature in the oasis region. During 2005–2021, LST fluctuated greatly in the northwestern part of the UANSTM, with LST showing an increasing trend during both daytime and nighttime, and the warming rate was more intense during daytime than nighttime. The increasing trend of LST in Urumqi, Changji Hui Autonomous Prefecture, Shihezi, and Wujiaqu is very significant and will remain consistent in the future. Precipitation, DEM, and AOD are the most important influencing factors of LST in the UANSTM, where precipitation and DEM are negatively correlated with LST, and AOD is positively correlated with LST. Land cover factors (LULC, NDVI,, and NDBSI) are the next most influential, and socioeconomic factors (NTL, GDP, and POP) are the least influential. The results of this study can provide a scientific reference for the conservation and sustainable development of the ecological environment of the UANSTM.

1. Introduction

Land surface temperature (LST) is a key factor affecting the surface–atmosphere energy exchange and water cycle, and it is also an important manifestation of the surface energy balance [1,2,3]. With the rapid advancement of global urbanization, a large number of impervious surfaces dominated by cement and asphalt have destroyed the original natural landscape and caused changes in surface parameters such as surface albedo and emissivity, which make the LST changes in urban areas more obvious [4,5,6,7]. The change of urban LST will not only reduce the comfort of people’s living environment and endanger human health, it will also have a serious impact on the local climate and eco-environment evolution of the city [8,9]. The urban agglomeration is the space carrier for the development of regional economy, industry, and urbanization to a certain advanced level, which has become the most prominent manifestation of global urbanization [10,11]. At present, many large urban agglomerations have been formed in China, such as Beijing-Tianjin-Hebei [12], Yangtze River Delta [13], Pearl River Delta [14], and Chengdu-Chongqing [15] urban agglomerations. Therefore, it is very necessary to understand the spatial distribution pattern, changing trend, and influencing factors of LST in urban agglomerations for protecting the eco-environment and sustainable economic development.
In recent years, a large number of scholars have found fruitful research results in the spatial-temporal distribution pattern [16,17] and influencing factors [18,19,20,21] of LST, and many spatial statistical models [22,23] and landscape ecology theories [24,25,26] have also been introduced into the study of the evolution and influencing mechanism of LST. Previous studies have indicated that land cover is considered a key factor affecting LST [27]. The impervious surface will increase the sensible heat flux of the surface and increase the heat radiation energy outward from the surface, which will lead to the increase of LST [28]. Water and green space are the natural radiators of the city. Vegetation can effectively reduce the LST through transpiration and shading [29]. Water has a high specific heat capacity and can remove a large amount of heat through evapotranspiration, so it can keep a lower temperature [30]. The normalized difference vegetation index (NDVI), normalized difference water index (NDWI), and normalized differential build-up and bare soil index (NDBSI) can accurately characterize land cover information, and they have shown a good linear relationship with LST [31,32]. In addition, some climatic, socioeconomic, and topographic factors also have an assignable impact on LST [33,34]. For example, DEM and gradient have a significant negative linear correlation with LST [35]. Population density (POP), gross domestic product (GDP), nighttime lights (NTL), aerosol optical depth (AOD), and precipitation also have strong correlations with LST [31,36,37,38].
In previous studies, Landsat series data have often been used to analyze changes in the spatial distribution pattern of LST based on one or several sensor transit moments. Although Landsat data have a high spatial resolution, the temporal resolution is low and often limited by the cloud coverage conditions at the moment of sensor transit, making it difficult to provide continuous and stable observations of large-scale LST. In studies of influencing factors, simple linear regression models [23,39,40] are also mostly used as a research method to explore the relationship between single or multiple influencing factors and LST. However, the above methods cannot fully reveal the nonlinear relationship between each influencing factor and LST, nor can they analyze the effect on LST when different influencing factors act together. In addition, the selection of influencing factors mostly prefers the land cover factors and lacks consideration of climatic, socioeconomic, and topographical factors. Therefore, in this study, MODIS LST data with high temporal resolution were selected to analyze the spatial-temporal change trend of LST, and the geo-detector model (GDM) [41] was applied to explore the effects of various influencing factors on LST. The GDM can quantitatively reveal the linear and nonlinear relationships between each influencing factor and LST, and measure the impact on LST when different influencing factors act together [42].
This study was conducted in the urban agglomeration on the northern slope of Tianshan Mountains (UANSTM). The UANSTM is the most urbanized and most densely populated and industrialized area in Xinjiang, China [43]. With one-fifth of Xinjiang’s population and half of its GDP, it has an irreplaceable role in Xinjiang [44]. However, since the 1980s, urbanization has led to the deterioration of the eco-environment, and a large amount of ecological land has been converted into construction land, and the desert-oasis ecotone has continued to shrink [45,46]. In addition, the current studies on LST of urban agglomeration are still focused on coastal and developed regions, while very little attention has been paid to the study of urban agglomerations in inland arid and undeveloped regions, which is an important reason why the UANSTM has been chosen as the study area.
Therefore, to address the above issues, the main objective of this study is to analyze the spatial distribution pattern, fluctuation characteristics, and change trends of LST in the UANSTM based on MODIS LST data with high temporal resolution, and to analyze the effects of multiple influencing factors on LST using the geo-detector model. This study will help planners to provide decision-making and support for improving the quality of human habitat and surface thermal environment issues.

2. Study Area and Data Source

2.1. Study Area

The UANSTM is an emerging urban agglomeration in the inland arid region of northwest China [47]. It is located at the northern foot of the Tianshan Mountains and the southern edge of the Junggar Basin in Xinjiang, and it includes nine cities: Urumqi, Wujiaqu, Changji Hui Autonomous Prefecture, Turpan, Shihezi, Karamay, Kuitun, Usu, and Shawan (Figure 1). The UANSTM is a typical oasis urban agglomeration with scarce water resources, fragile eco-environment, and sensitivity to climate change [48]. Compared with urban agglomerations in developed countries and coastal China, the UANSTM is particularly uncompetitive in terms of urbanization, resources, and eco-environment carrying capacity [44,49]. The Chinese government has attached great importance to the eco-environment issues in recent years. To enhance the competitiveness of the UANSTM and promote its sustainable development transformation, it is necessary to search for the eco-environment problems in various aspects. LST, as an important parameter in the eco-environment, needs urgent attention. Therefore, under the influence of intense human activities, it is of very high application value to understand the changing characteristics and influencing factors of LST in the UANSTM.

2.2. Data Collection and Preprocessing

2.2.1. LST Data

The LST data for this study used all MODIS Aqua LST products MYD11A2 (version 6) covering the UANSTM from 2005 to 2021, which required track numbers “h23v04” and “h24v04” for mosaicking, with a total data volume of 1564 remote sensing images (https://ladsweb.modaps.eosdis.nasa.gov/, accessed on 15 April 2022). Each image contains “LST_Day_1km” and “LST_Night_1km” bands, which are daytime and nighttime LST data, respectively. The MYD11A2 product is calculated from the thermal infrared bands of MODIS 31 and 32 channels using a split-window algorithm with a spatial resolution of 1 km. This product is a composite LST product obtained by averaging 8 days of clear-sky LST, which can effectively reduce the effect of cloudiness, and the product has been validated with a series of accuracies and has been widely used in regional- or global-scale LST studies [50,51,52,53]. We calculated the average daytime and nighttime LST images for each year from 2005 to 2021 for the UANSTM and counted some details of the images, as shown in Table 1.

2.2.2. Influencing Factors’ Data

Based on previous studies and data availability, we considered four types of influencing factors: land cover, climate, socio-economic, and topographic factors, and finally selected 10 influencing factors, including: land use/land cover (LULC), NDVI, NDBSI, precipitation, AOD, NTL, GDP, POP, DEM, and gradient (Table 2). The LULC data are a 30 m annual Chinese land cover dataset from 1990 to 2019 produced by Yang et al. [54] (https://zenodo.org/, accessed on 19 April 2022). We reclassified it into six types, including: cropland, woodland, water area, grassland, construction land, and unused land. The NDVI, AOD, and NDBSI data come from the MODIS vegetation index product MYD13A2, terrestrial aerosol optical depth product MCD19A2, and surface reflectance product MYD09A1, respectively (https://ladsweb.modaps.eosdis.nasa.gov/, accessed on 15 April 2022). The calculation method of NDBSI is detailed in [47]. The GDP and POP data were from the Resource and Environmental Science and Data Center (https://www.resdc.cn/, accessed on 15 April 2022). The NTL data were from the NPP-VIIRS NTL monthly product produced by the Colorado School of Mines (https://eogdata.mines.edu/products/vnl/, accessed on 15 April 2022). The precipitation data were obtained from the monthly precipitation raster data of the National Tibetan Plateau Data Center (http://data.tpdc.ac.cn/, accessed on 16 April 2022). The DEM and gradient data were obtained from the Geospatial Data Cloud (http://www.gscloud.cn/, accessed on 16 April 2022).
We resampled all influencing factor data to a 1 km resolution to ensure that all data had the same spatial resolution. We calculated annual average rasters for NDVI, NDBSI, precipitation, AOD, and NTL to ensure that all data were annual averages.

3. Methodology

In this study, the spatial distribution pattern of LST was first analyzed by grading the LST using the standard deviation method. Next, the absolute variability, Sen’s slope analysis, Mann–Kendall trend test, and Hurst index were used to explore the fluctuation characteristics and change trends of LST. Finally, the effects of land cover factors, climatic factors, socioeconomic factors, and topographic factors on LST were analyzed based on the geo-detector model (the overall research flow is shown in Figure 2).

3.1. Land Surface Temperature Classification

To accurately describe the spatial pattern differences of LST, the standard deviation method was used to divide LST into five classes: extremely high temperature (EHT), high temperature (HT), medium temperature (MT), low temperature (LT), and extremely low temperature (ELT) (Table 3). The standard deviation method can well-characterize the concentration and volatility of LST by combining the mean and different standard deviation multiples [40].

3.2. Absolute Variability (AV)

The AV reflects the absolute deviation between the data and the mean over a certain time period, which means that it can reflect the degree of fluctuation of the data over a certain time period [52]. In this study, AV was used to characterize the fluctuations of LST in the daytime and nighttime between 2005 and 2021. The expression for AV is as follows:
A V = i = 1 n x i x m e a n n
where xi is the LST observation in year i, xmean is the average LST from 2005 to 2021, and n is the number of years (n = 17). The AV in Equation (1) is proportional to the fluctuation strength.

3.3. Trend Analysis Methods

3.3.1. Sen’s Slope Analysis

Sen’s slope analysis is a method used to estimate the changing trend of the data over a time series, which can effectively avoid the effect of missing values and outliers. Its fundamental principle is to first calculate the slope between all adjacent data in the time series data, and then take the median value of the slope as the value of the changing trend [55]. The basic equation is written as:
s l o p e = m e d i a n L S T j L S T i j i   ,     j > i
where median is the median function, and LSTj and LSTi are the observed values at moments j and i in the time series. In Equation (2), when the slope < 0, it indicates a decreasing trend, and when the slope > 0, it indicates an increasing trend.

3.3.2. Mann–Kendall Trend Test

The Mann–Kendall trend test is a nonparametric test method that is often used in conjunction with Sen’s slope analysis to determine the significance of a changing trend [56,57]. As with Sen’s slope analysis, the Mann–Kendall trend test is insensitive to missing and outlier data. Its normal distribution statistic, S, is calculated as follows:
S = i = 1 n 1 j = i + 1 n s g n L S T j L S T i
where LSTi and LSTj are the LST observations in year i and j, respectively, and sgn in Equation (3) is a symbolic function whose expression is presented in Equation (4):
s g n L S T j L S T i = 1 ,     L S T j L S T i > 0 0 ,     L S T j L S T i = 0 1 ,     L S T j L S T i < 0
From Equation (4), it can be seen that the value of the sgn function varies between −1, 0, and 1 based on (LSTjLSTi). The variance, Var(S), of S can be calculated based on Equation (5) as follows:
V a r S = n n 1 2 n + 5 18
Based on S and Var(S) in Equations (3) and (5), the standard orthogonal statistical distribution, Z, of the Mann–Kendall trend test can be calculated as shown in Equation (6):
Z = S 1 V a r S ,         S > 0                 0             ,     S = 0 S + 1 V a r S ,       S < 0
Referring to previous studies, this study uses 95% and 99% confidence intervals, which means that the changing trend is very significant when Z ≥ 2.58 and Z ≤ 2.58, significant when 1.96 ≤ Z < 2.58 and −2.58 < Z ≤ −1.96, and otherwise the trend is not significant. We superimposed the results of Sen’s slope analysis and the Mann–Kendall trend test to classify the changing trends of LST into six types (Table 4).

3.3.3. Hurst Index

The Hurst index can effectively describe the persistence of changes in time series data. It reflects the interrelationship between before and after changing trends, that is, the past changing trends can affect the current changing trends, and the current changing trends can affect the future changing trends [58,59]. Therefore, the Hurst index can be used to analyze the future trend of LST. The basic principle of the Hurst index is as follows:
Define a time series as: {LST(t), t = 1, 2, … n}, and for any integer α ≥ 1, the average series is shown in Equation (7):
L S T ¯ α = 1 α t = 1 α L S T t         α = 1 , 2 , , n
The cumulative deviation is calculated as shown in Equation (8):
X t , α = t = 1 t L S T t L S T ¯ α         1   t α  
Define the range of extreme difference, R(α), as in Equation (9):
R α = m a x 1 t τ X t , α m i n 1 t τ X t , α         α = 1 , 2 , , n
Calculate the standard deviation, S(α), as shown in Equation (10):
S α = 1 α t = 1 α L S T t L S T ¯ α 2 1 / 2         α = 1 , 2 , , n  
Finally, the Hurst index is calculated by Equations (11) and (12):
R α S α = c α H
log R S n = a + H × log n
where H is the Hurst index, c is the relation constant, and a is the intercept. The H value can be obtained by fitting Equation (12) using the least squares method. When H > 0.5, it indicates that the future changing trend of LST is most likely to be the same as the present, when H < 0.5, it indicates that the future changing trend of LST may be opposite to the present, and when H = 0.5, it indicates that the future changing trend of LST is random. In this study, H < 0.5 is defined as the future changing trend is unknown. To understand the future changing trend of LST, we superimpose the Hurst index and the LST changing trend types to obtain the future changing trend types of LST (Table 4).

3.4. Geo-Detector Model Based on Optimal Parameters

The GDM is a spatial statistical method used to explore influencing factors based on spatial differentiation theory. The GDM measures the degree of correlation between independent variables and dependent variables by measuring the degree of spatial consistency between them, which can effectively avoid the multicollinearity problem between dependent variables [41]. This study used four modules of GDM: factor detector, interaction detector, risk detector, and ecological detector, to analyze the effect of each influencing factor on LST. Due to the availability of data, this study only used 2019 as an example to analyze the influencing factors.
When using GDM, continuous variables must be discretized, and the discrete method can have a large impact on the results. The commonly used data discretization methods are natural breaks, equal interval, quantile, standard deviation, and geometric interval. However, most of the previous studies rely on personal experience in choosing the discretization methods, which is rather subjective, while the GD package of R studio can select the optimal method based on q values to obtain the optimal parameters for GDM. Therefore, in this study, firstly, 2 km grid points were established according to the study area to extract the values of LST and each influence factor, and then the optimal discretization of continuous type variables and calculation of GDM results were performed by GD package (https://cran.r-project.org/web/packages/GD/, accessed on 20 April 2022).

3.4.1. Factor Detector

The factor detector measures the impact of the independent variable on the dependent variable by detecting the extent to which the dependent variable explains the spatial heterogeneity of the independent variable, and this degree of impact is quantified using the q value. The expression is as follows:
q = 1 h = 1 L N h σ h 2 N σ 2
where the range of q is [0,1]. A larger q value indicates that the influencing factor has a greater impact on the LST. h = 1, 2, …, L is the classification or partition of the independent or dependent variable, Nh and N are the number of cells in layer h and the whole region, respectively, and σ h 2 and σ2 are the variance of the dependent variable in layer h and the whole region, respectively.

3.4.2. Interaction Detector

The interaction detector identifies interactions between different influencing factors by assessing the magnitude of the explanatory power of the LST when influencing factors X1 and X2 act together. The identification method is to first calculate the q values of the two independent variables X1 and X2 on LST: q(X1) and q(X2), and the q values of the two independent variables on LST when they interact: q(X1X2), and then compare q(X1), q(X2), and q(X1X2) to classify the different interaction types (Table 5).

3.4.3. Risk Detector

The risk detector is used to judge whether the mean values of dependent variables in different independent variable intervals are significantly different, and the t statistic is used to judge:
t = M ¯ h = 1 M ¯ h = 2 V a r M ¯ h = 1 n h = 1 + V a r M ¯ h = 2 n h = 2 1 2
where M ¯ h is the average value in the h interval, nh is the number of samples in the interval h, and Var() is the variance function.

3.4.4. Ecological Detector

The ecological detector is used to compare whether the influence of the two influencing factors X1 and X2 on the spatial distribution of the dependent variable is significantly different, measured by the F statistic:
F = N X 1 N X 2 1 S S W X 1 N X 2 N X 1 1 S S W X 2
where NX1 and NX2 denote the sample sizes of influences X1 and X2, and SSWX1 and SSWX2 denote the sum of variances within the X1 and X2 layers, respectively.

4. Results

4.1. Spatial Distribution Pattern of LST

To avoid chance, we calculated the average values of daytime and nighttime LST from 2005 to 2021 for the urban agglomeration on the northern slopes of Tianshan to analyze the spatial distribution pattern of LST during the daytime and nighttime (Figure 3). From the results, it is clear that the daytime LST ranged from −12.61 to 42.64 °C. The area proportions of HT and EHT were 33.28% and 2.23%, respectively, both mainly clustered in Turpan city. The area proportions of LT and ELT were 18.38% and 8.64%, respectively, mainly distributed in the high-altitude areas of the Tianshan Mountains in the southwestern edge and central part of the study area, as well as in the oasis area in the central part. The LST at night ranged from −25.58 to 10.61 °C. The proportions of HT and EHT areas were 20.43% and 5.98%, respectively, and both were equally clustered within Turpan city. The area of EHT almost doubled at night compared to daytime, while the area of HT significantly decreased. The area proportions of LT and ELT at night were 11.01% and 7.30%, respectively. During the nighttime, the area and spatial distribution of the ELT were almost the same as during the daytime. However, the area of LT was significantly reduced, mainly because the large area of LT in the oasis region in the central part of the study area changed to MT and HT at night. In general, the LST classes are dominated by MT and HT in the UANSTM, both during daytime and nighttime. The high-temperature classes (EHT and HT) were always clustered in Turpan city. In the oasis area of the study area, there was a significant difference in LST classes between daytime and nighttime.

4.2. The Fluctuation Characteristics of LST

We analyzed the fluctuation characteristics of daytime and nighttime LST of the UANSTM from 2005 to 2021 by absolute variability (AV), and higher values of AV represent more intense fluctuations of LST. From Figure 4, the AV ranged from 0.27 to 5.18 °C with a mean value of 1.09 °C during the daytime and from 0.26 to 2.89 °C with a mean value of 0.59 °C during the nighttime, indicating that LST fluctuations are more intense during daytime than nighttime. The daytime high AV areas are mainly clustered in Karamay, Shawan, and Changji Hui Autonomous Prefecture in the northwestern part of the study area, especially in Karamay, indicating strong LST fluctuations in these areas during the daytime. The nighttime high AV areas are mainly clustered in Usu, Karamay, Shawan, and parts of Changji Hui Autonomous Prefecture and Urumqi, especially in Usu, indicating that the LST fluctuations are stronger in these areas during the nighttime.

4.3. The Changing Trend of LST

4.3.1. The Changing Trend of LST from 2005 to 2019

The results of the Sen’s slope analysis (Figure 5) showed that the trend of LST during the daytime in the UANSTM ranged from −0.65 to 0.96 °C-yr−1 during 2005–2021, with a mean value of 0.05 °C-yr−1, and the proportion of the area with an increasing trend was 76.83%, while the proportion of the area with a decreasing trend was only 23.17%. The nighttime LST trend ranged from −0.22 to 0.49 °C-yr−1 with a mean value of 0.03 °C-yr−1, and the proportion of the area with an increasing trend was 86.59%, while only 13.41% of the area had a decreasing trend. In general, the UANSTM showed an increasing trend during both daytime and nighttime, and the increasing trend was stronger during daytime than nighttime.
The results of the Sen’s slope analysis and the Mann–Kendall trend test were superimposed, and the trends were classified into the following six levels according to Table 4: “Very significantly decreased”, “Significantly decreased”, “Not significantly decreased”, “Not significantly increased”, “Significantly increased”, and “Very significantly increased”. From the results (Figure 6), the trend of LST in most of the area of the UANSTM was insignificant (“Not significantly decreased” and “Not significantly increased”) in both daytime and nighttime, and only a small portion of the area had a significant or very significant trend (“Very significantly decreased”, “Significantly decreased”, “Very significantly increased”, and “Significantly increased”). During the daytime, the proportion of the area with an insignificant trend was 84.29%, while the proportion of the very significant and significant area was only 14.25%. Among them, the proportions of “Significantly increased” and “Very significantly increased” were 9.55% and 5.39%, respectively, mainly in Urumqi City, Wujiaqu City, Changji Hui Autonomous Prefecture, and Shihezi City, and the proportions of “Significantly decreased” and “Very significantly decreased” were 0.22% and 0.55%, mainly in Kuitun City. During the nighttime, the proportion of areas with insignificant trends was 93.43%, and only 6.57% of the areas were very significant and significant. The proportions of “Very significantly increased” and “Significantly increased” were 5.19% and 1.18%, respectively, with a discrete distribution and small clusters in Urumqi, the eastern part of Changji Hui Autonomous Prefecture, and Turpan City. The proportions of “Very significantly decreased” and “Significantly decreased” areas did not exceed 0.1%, and there were no obvious clustering areas.

4.3.2. The Future Changing Trend of LST

From the results of the Hurst index (Figure 7), during the daytime, the percentage of the area with sustained future trend (Hurst > 0.5) was 50.44%, and the percentage of the area with an unknown future trend (Hurst ≤ 0.5) was 49.56%. During the nighttime, the proportion of the area with a sustained future trend was 66.85%, and the proportion of the area with an unknown future trend was only 33.15%.
The Hurst index was coupled with the trend level, and the future trend of LST was classified into seven levels based on Table 4: “Sustained very significant decrease”, “Sustained significant decrease”, “Sustained no significant decrease”, “Sustained very significant increase”, “Sustained significant increase”, and “Sustained no significant increase”. During the daytime, 13.49% of the area showed a sustained decrease in the future trend of LST (“Sustained very significant decrease”, “Sustained significant decrease”, and “Sustained no significant decrease”), and 36.95% of the area showed a sustained increase (“Sustained very significant increase”, “Sustained significant increase”, and “Sustained no significant increase”). Among them, “Sustained very significant increase” and “Sustained significant increase” are mainly distributed in Urumqi, Wujiaqu, Changji Hui Autonomous Prefecture, and Shihezi, which means that the LST in these regions will continue to significantly and consistently increase in the future. During the nighttime, the proportion of the area that sustained an increase was 62.33%, while the proportion of the area that sustained a decrease was only 2.73%. The proportion of the area that “Sustained no significant increase” was 56.02%, which means that in the future, although the LST in most of the areas at night has a sustained increasing trend, the increasing trend is not significant (Figure 8).

4.4. Analysis of the Influencing Factors of LST

4.4.1. Discretization of Continuous Variables

When using the geo-detector model for influencing factor exploration, the independent variables are required to be type values. In this study, all variables except LULC are continuous type variables, thus requiring discretization of NDVI, NDBSI, precipitation, AOD, GDP, POP, DEM, and gradient data. We first defined five discretization methods, including: natural breaks, equal interval, quantile, standard deviation, and geometric interval methods. Then, we selected the optimal discretization method and the number of discretization intervals for daytime and nighttime based on the maximum q value by using the optidisc (optimal discretization for continuous variables and visualization) function in the GD package of R. Table 6 shows the selected discretization methods and the number of discretization intervals for each influence factor during the daytime and nighttime.

4.4.2. Impact of Each Influencing Factor on LST in the UANSTM

The factor detector results show the degree of impact of each influencing factor on LST. From the results (Figure 9), the impacts of each influence factor on LST differed between daytime and nighttime, and the impacts were greater in the daytime than in the nighttime.
In terms of the influencing factors, regardless of daytime or nighttime, precipitation had the greatest influence on LST, followed by DEM and AOD, and the q values of these three influencing factors were above 0.5 in the daytime and above 0.4 in the nighttime, which were obviously higher than other influencing factors. This indicates that precipitation, DEM, and AOD were the main factors influencing LST in the UANSTM. In terms of the types of influencing factors, regardless of daytime or nighttime, climatic factors (precipitation and AOD) had the greatest influence on LST, with average q values of 0.6476 and 0.5074 for daytime and nighttime, respectively, followed by topographic factors (DEM and gradient) and land cover factors (LULC, NDVI, and NDBSI). The socioeconomic factors (NTL, GDP, and POP) had the least effect on LST, with mean q values of 0.1747 and 0.0838 for daytime and nighttime, respectively.

4.4.3. Interaction between Influencing Factors

In this study, the interaction detector was used to detect the combined effect of two influencing factors on LST and the interaction relationship (Figure 10). In both the daytime and the nighttime, the influencing factors only showed two types of bivariate enhancement and nonlinear enhancement, there were no independent influencing factors, and the combined impact of the two influencing factors was greater than the impact of a single influencing factor on LST.
In the daytime, only the interaction relationships of NTL ∩ DEM and NTL ∩ gradient showed nonlinear enhancement, while the rest showed bivariate enhancement. The interaction of precipitation with other influencing factors was very strong, with q values above 0.7. The largest q value for precipitation ∩ LULC (q = 0.8325) indicated that the interaction between precipitation and LULC had the greatest effect on daytime LST in the UANSTM. In the nighttime, the nonlinear enhancement relationship only occurred between socioeconomic factors (NTL, GDP, and POP) or between socioeconomic factors and other influencing factors. The interaction effect of precipitation with other influencing factors remained very strong. The interaction q value between precipitation and DEM was the largest (q = 0.7481), which means that the interaction between precipitation and DEM had the largest effect on nighttime LST in the UANSTM.

4.4.4. Risk Interval Analysis of Each Influencing Factor

In this study, we used the risk detector module to analyze the risk intervals of each influencing factor during daytime and nighttime in the UANSTM. The risk detector can count the average LST within each discretized interval for each influence factor, and we defined the interval with the highest average LST as the risk interval. From the results of the risk detector (Figure 11), precipitation and DEM were negatively correlated with LST, while AOD was positively correlated with LST, both during the daytime and the nighttime. The risk interval of LST for precipitation, DEM, and AOD during daytime occurred in the intervals of 97–664 mm for precipitation, −227–−48.5 m for DEM, and 0.285–4 for AOD. The nighttime risk interval was within the intervals of 97–588 mm for precipitation, −227–163 m for DEM, and 0.285–4 for AOD. From the LULC perspective, the LULC daytime risk interval occurred in the unused land type because the unused land is basically desert and the daytime temperature is very high. The nighttime risk interval occurred in the construction land type because the desert cools quickly at night, while the temperature is higher in the construction land due to various human activities. The mean LST fluctuated with NDVI, NDBSI, and gradient. The NDVI and gradient daytime and nighttime risk intervals were the same, both occurring in the 0.01–0.06 interval for NDVI and the 2.43–7.75° interval for gradient. NDBSI daytime and nighttime risk intervals differed little, occurring in the intervals 0.612–0.647 and 0.607–0.648, respectively.

4.4.5. Differences in the Impact of Influencing Factors on LST

The ecological detector can detect whether the influence of each influencing factor on the spatial distribution of LST is significantly different. The results showed (Figure 12) that, regardless of daytime or nighttime, most of the influencing factors had significant differences in their effects on LST, and only a few influencing factors had no significant differences. Among them, the effects of precipitation and DEM were not significantly different. That is to say, the influence mechanisms of precipitation and DEM on LST were similar. From the spatial distribution, precipitation and DEM had a high consistency in spatial distribution, where the higher the DEM, the higher the rainfall, and from the results of the risk detector (Figure 11), it can be seen that there was a negative linear relationship between precipitation and DEM with LST, where the higher the precipitation and DEM, the lower the LST. Therefore, there was no significant difference between the effects of precipitation and DEM on LST. The socioeconomic factors were very closely related to each other, and their spatial distribution was highly consistent, with high-value regions of NTL indicating that their GDP and POP were also high. Therefore, there was no significant difference between their effects on LST.

5. Discussion

5.1. Spatial Distribution Patterns and Changing Trend of LST

Turpan is an extremely hot region in China with a mean interannual LST maximum above 75 °C [60,61]. This is consistent with our findings, where we found that Turpan is the main concentration of HT and EHT in both the daytime and nighttime. The main reason is that most of Turpan is a desert, which leads to direct solar radiation and fast warming in the daytime. In addition, Turpan is a typical graben basin, resulting in poor air mobility and slow heat dissipation. The cold island phenomenon in oases was first proposed in [62]. The UANSTM is a typical oasis urban agglomeration, and both the spatial distribution pattern of LST and the risk detector results of LULC indicate a cold island phenomenon during the daytime and a heat island phenomenon during the nighttime in the oasis region. This is consistent with the results of [63] based on MODIS LST data and meteorological station information for the Tarim Basin. However, some studies have shown that in some developed urban agglomerations in China, such as the Beijing-Tianjin-Hebei urban agglomeration [51] and the Pearl River Delta urban agglomeration [33], they exhibit a strong heat island phenomenon in both the daytime and the nighttime. The main reason for this difference is that the periphery of the oasis is mostly desert, bare land, and low-coverage grassland, which have low specific heat capacity and warm up fast by direct solar radiation in the daytime. The interior of the oasis has more artificial blue and green landscapes, which have a high specific heat capacity and warm up slowly. Therefore, in the daytime, the temperature inside the oasis is lower than that at the periphery of the oasis, while the opposite is true in the nighttime [64]. In addition, in arid regions, agriculture is dominated by irrigated farming, and some studies have shown that irrigation can enhance the cold island phenomenon in oases [65,66].
Since the 1960s, the biggest feature of global climate change has been large-scale warming [67]. However, intermittent global warming has also been found, which could also explain the fluctuation of the interannual mean LST to some extent [68]. In this study of change trends, we found that LST exhibited an increasing trend in both the daytime and the nighttime, with an average increase of 0.05 °C per year in the daytime and 0.03 °C per year in the nighttime. This is consistent with the change trend of LST in the whole of China. The authors of [53] explored the change trend of LST in China from 2003 to 2019 and found that China’s LST showed an increasing trend in both the daytime and the nighttime, and the daytime was stronger than the nighttime. It is noteworthy that there was a significant increasing trend of daytime LST in Urumqi, Wujiaqu, Changji Hui Autonomous Prefecture, and Shihezi during 2005–2021, and the same trend will remain in the future. This may be due to the following two reasons: (1) The “Urumqi-Changeji-Wujiaqu” metropolitan area and Shihezi have developed rapidly in recent years, and strong human activities and intense land use changes have led to the increase in LST. (2) The “Urumqi-Changeji-Wujiaqu” metropolitan area and Shihezi are the key development areas in the development plan of the UANSTM and will continue to lead the development of the Tianshan north slope urban agglomeration in the future [44].

5.2. Influencing Factors of LST

Based on the results of the GDM, we found that precipitation has the greatest impact on the LST of the UANSTM. Water shortage has been the most central problem limiting the development of the Xinjiang region [43]. Water is a natural heat sink due to its high specific heat capacity and slow warming, and precipitation can directly carry away a large amount of heat from the surface and increase surface runoff. There was also a significant correlation between precipitation and vegetation cover, so precipitation can also indirectly affect LST by influencing vegetation [69,70]. We found that AOD also has a strong impact on the LST of the UANSTM. Aerosol particles can influence the absorption of solar short-wave radiation at the surface and the radiation exchange with the atmosphere, which leads to changes in LST [71]. The Gurbantunggut desert to the north of the UANSTM is one of the major global sources of sand and dust aerosol emissions, which, together with the large amount of anthropogenic aerosol particle emissions, has led to serious aerosol particle pollution in the UANSTM [72,73]. The influence of DEM on LST is second only to precipitation, which may be related to the large difference in altitude between the UANSTM. The difference between the highest and lowest elevations in the UANSTM is above 5000 m. In the high-altitude region of the Tianshan Mountains, which is covered with snow all year round, LST is extremely low, while Turpan, which has the lowest altitude, has the highest LST in China. We also found that socioeconomic factors (NTL, GDP, and POP) have very little impact on LST in the UANSTM. This is contrary to the findings of [33] for the PRD urban agglomeration, who found that NTL, GDP, and POP have a very significant effect on LST. The reason for this discrepancy may be that the UANSTM, as an urban agglomeration in the northwestern frontier of China, has a wide spatial extent but a low urbanization rate, a relatively backward economic level, and a sparse population compared with other developed urban agglomerations, resulting in a smaller impact on LST.

5.3. Research Limitations and Future Work

Any geographical element has scale effects, and LST is no exception [74]. This means that there may be significant differences in study results at different temporal and spatial scales. This study deeply analyzed the spatial distribution pattern, changing trend, and various influencing factors of LST in the UANSTM, and provided a theoretical basis for the prevention and management of urban thermal environment problems, but there are still some limitations. The first is that this study analyzed the trend of LST only from the interannual scale, and in the future, we need to consider the variation on the seasonal scale, which will be more instructive for agricultural production [75]. Second, although the MODIS data have high temporal resolution, the spatial resolution is poor, which makes it difficult to reflect the details of changes when studying the LST within cities. Therefore, combining high spatial resolution data with downscaling fusion of MODIS data to obtain LST data with high spatial and temporal resolution will be an important direction for future urban surface thermal environment studies [76,77]. Third, although we considered 10 influencing factors in terms of climate, land cover, socioeconomics, and topography, LST is influenced by multiple factors, and more influencing factors need to be considered in future studies. Finally, we only used the geo-detector model to explore the influencing factors, although it possesses higher accuracy than the ordinary linear model. However, other methods may work better, such as random forest regression, logistic regression, and geographically weighted regression models [78,79,80].

6. Conclusions

This study analyzed the spatial distribution pattern of LST based on MODIS LST data, with the UANSTM as the study area. Absolute variability, Sen’s slope analysis, and the Mann–Kendall trend test were used to analyze the fluctuation characteristics and change trends of LST from 2005 to 2021, and the future change trends of LST were explored with the Hurst index. Finally, the influence of 10 factors, such as LULC, NDVI, NDBSI, and precipitation, on LST was analyzed by using the geo-detector model based on the optimal parameters. The detailed findings of this study are as follows:
1.
In terms of the spatial distribution pattern, both daytime and nighttime, the LST classes of the UANSTM are dominated by MT and HT, with EHT and HT clustered in Turpan city. As far as the diurnal differences are concerned, in the oasis region, the cold island feature is observed during the daytime and the heat island feature is observed at night.
2.
During 2005–2021, the northwestern part of the urban agglomeration on the northern slopes of the Tianshan Mountains fluctuated greatly, especially in Karamay and Usu. The LST showed an increasing trend in both the daytime and the nighttime, with an increase rate of 0.05 and 0.03 °C/yr−1, respectively. The increasing trend of LST in Urumqi, Changji Hui Autonomous Prefecture, Shihezi, and Wujiaqu was very significant during the daytime and still showed a significant increasing trend in the future, which needs to be investigated.
3.
The climatic and topographic factors of precipitation, DEM, and AOD are the main factors affecting the LST of the UANSTM, and they all have q values above 0.5 during the daytime and above 0.4 during the nighttime. The effects of land cover factors (LULC, NDVI, and NDBSI) was the second most important, and socioeconomic factors (NTL, GDP, and POP) had the least influence on LST. The interactions between the influencing factors all showed enhancement types (nonlinear enhancement and bivariate enhancement), and one of the two influencing factors that showed the nonlinear enhancement type must be a socioeconomic factor.
4.
Precipitation and DEM showed a negative linear correlation with LST, while AOD showed a positive linear correlation with LST. The mean LST values for each interval of precipitation and DEM decreased as the values within the interval increased, and the risk interval for precipitation and DEM occurred in the lowest value interval, while the mean LST values for each interval of AOD increased as the values within the interval increased, and the risk interval for AOD occurred in the highest value interval.

Author Contributions

Writing—original draft preparation, H.L.; funding acquisition, A.K.; software, H.M.; methodology, Y.Z.; data curation, X.Z.; visualization, B.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Third Xinjiang Scientific Expedition Program, grant number 2021xjkk0905, and the Postgraduate Research and Innovation Project of Xinjiang Uygur Autonomous Region, grant number XJ2022G209.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available within the article.

Acknowledgments

We thank the three anonymous reviewers for their constructive comments and suggestions that have helped to improve the original manuscript. Thanks also to the editorial staff. Many thanks to the team at the Oasis Urban Remote Sensing Laboratory in the Arid Zone for their hard work.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sobrino, J.A.; Jiménez-Muñoz, J. Minimum configuration of thermal infrared bands for land surface temperature and emissivity estimation in the context of potential future missions. Remote Sens. Environ. 2014, 148, 158–167. [Google Scholar] [CrossRef]
  2. Li, Z.; Tang, B.; Wu, H.; Ren, H.; Yan, G.; Wan, Z.; Trigo, I.F.; Sobrino, J.A. Satellite-derived land surface temperature: Current status and perspectives. Remote Sens. Environ. 2013, 131, 14–37. [Google Scholar] [CrossRef]
  3. Duan, S.; Li, Z.; Wu, H.; Leng, P.; Gao, M.; Wang, C. Radiance-based validation of land surface temperature products derived from Collection 6 MODIS thermal infrared data. Int. J. Appl. Earth Obs. Geoinf. 2018, 70, 84–92. [Google Scholar] [CrossRef]
  4. Grimm, N.B.; Faeth, S.H.; Golubiewski, N.E.; Redman, C.L.; Wu, J.; Bai, X.; Briggs, J.M. Global change and the ecology of cities. Science 2008, 319, 756–760. [Google Scholar] [CrossRef]
  5. Liu, G.; Zhang, Q.; Li, G.; Doronzo, D.M. Response of land cover types to land surface temperature derived from Landsat-5 TM in Nanjing Metropolitan Region, China. Environ. Earth Sci. 2016, 75, 1386. [Google Scholar] [CrossRef]
  6. Shen, H.; Huang, L.; Zhang, L.; Wu, P.; Zeng, C. Long-term and fine-scale satellite monitoring of the urban heat island effect by the fusion of multi-temporal and multi-sensor remote sensed data: A 26-year case study of the city of Wuhan in China. Remote Sens. Environ. 2016, 172, 109–125. [Google Scholar] [CrossRef]
  7. Chang, Y.; Xiao, J.; Li, X.; Middel, A.; Zhang, Y.; Gu, Z.; Wu, Y.; He, S. Exploring diurnal thermal variations in urban local climate zones with ECOSTRESS land surface temperature data. Remote Sens. Environ. 2021, 263, 112544. [Google Scholar] [CrossRef]
  8. Cao, Q.; Yu, D.; Georgescu, M.; Wu, J.; Wang, W. Impacts of future urban expansion on summer climate and heat-related human health in eastern China. Environ. Int. 2018, 112, 134–146. [Google Scholar] [CrossRef]
  9. Wilson, J.S.; Clay, M.; Martin, E.; Stuckey, D.; Vedder-Risch, K. Evaluating environmental influences of zoning in urban ecosystems with remote sensing. Remote Sens. Environ. 2003, 86, 303–321. [Google Scholar] [CrossRef]
  10. Fang, C.; Yu, D. Urban agglomeration: An evolving concept of an emerging phenomenon. Landsc. Urban Plann. 2017, 162, 126–136. [Google Scholar] [CrossRef]
  11. Fang, C. Important progress and future direction of studies on China’s urban agglomerations. J. Geogr. Sci. 2015, 25, 1003–1024. [Google Scholar] [CrossRef]
  12. Shen, Z.; Xu, X.; Xu, S.; Sun, D. A comparative study of land development patterns and regional thermal environments (RTEs) in typical urban agglomerations of China and America: A case study of Beijing-Tianjin-Hebei (BTH) and Boswash. Sci. Total Environ. 2022, 803, 149735. [Google Scholar] [CrossRef]
  13. Li, W.; Han, C.; Li, W.; Zhou, W.; Han, L. Multi-scale effects of urban agglomeration on thermal environment: A case of the Yangtze River Delta Megaregion, China. Sci. Total Environ. 2020, 713, 136556. [Google Scholar] [CrossRef]
  14. Yu, Z.; Yao, Y.; Yang, G.; Wang, X.; Vejre, H. Spatiotemporal patterns and characteristics of remotely sensed region heat islands during the rapid urbanization (1995–2015) of Southern China. Sci. Total Environ. 2019, 674, 242–254. [Google Scholar] [CrossRef]
  15. Liu, X.; Ming, Y.; Liu, Y.; Yue, W.; Han, G. Influences of landform and urban form factors on urban heat island: Comparative case study between Chengdu and Chongqing. Sci. Total Environ. 2022, 820, 153395. [Google Scholar] [CrossRef]
  16. Imhoff, M.L.; Zhang, P.; Wolfe, R.E.; Bounoua, L. Remote sensing of the urban heat island effect across biomes in the continental USA. Remote Sens. Environ. 2010, 114, 504–513. [Google Scholar] [CrossRef]
  17. Chen, L.; Wang, X.; Cai, X.; Yang, C.; Lu, X. Seasonal variations of daytime land surface temperature and their underlying drivers over Wuhan, China. Remote Sens. 2021, 13, 323. [Google Scholar] [CrossRef]
  18. Gao, Y.; Zhao, J.; Han, L. Exploring the spatial heterogeneity of urban heat island effect and its relationship to block morphology with the geographically weighted regression model. Sustain. Cities 2022, 76, 103431. [Google Scholar] [CrossRef]
  19. El Kenawy, A.M.; Hereher, M.; Robaa, S.M.; McCabe, M.F.; Lopez-Moreno, J.I.; Domínguez-Castro, F.; Gaber, I.M.; Al-Awadhi, T.; Al-Buloshi, A.; Al Nasiri, N. Nocturnal surface urban heat island over Greater Cairo: Spatial morphology, temporal trends and links to land-atmosphere influences. Remote Sens. 2020, 12, 3889. [Google Scholar] [CrossRef]
  20. Du, J.; Xiang, X.; Zhao, B.; Zhou, H. Impact of urban expansion on land surface temperature in Fuzhou, China using Landsat imagery. Sustain. Cities 2020, 61, 102346. [Google Scholar] [CrossRef]
  21. Niu, L.; Zhang, Z.; Peng, Z.; Liang, Y.; Liu, M.; Jiang, Y.; Wei, J.; Tang, R. Identifying surface urban heat island drivers and their spatial heterogeneity in China’s 281 cities: An empirical study based on multiscale geographically weighted regression. Remote Sens. 2021, 13, 4428. [Google Scholar] [CrossRef]
  22. Ivajnšič, D.; Kaligarič, M.; Žiberna, I. Geographically weighted regression of the urban heat island of a small city. Appl. Geogr. 2014, 53, 341–353. [Google Scholar] [CrossRef]
  23. Wu, Y.; Hou, H.; Wang, R.; Murayama, Y.; Wang, L.; Hu, T. Effects of landscape patterns on the morphological evolution of surface urban heat island in Hangzhou during 2000–2020. Sustain. Cities Soc. 2022, 79, 103717. [Google Scholar] [CrossRef]
  24. Ye, H.; Li, Z.; Zhang, N.; Leng, X.; Meng, D.; Zheng, J.; Li, Y. Variations in the effects of landscape patterns on the urban thermal environment during rapid urbanization (1990–2020) in megacities. Remote Sens. 2021, 13, 3415. [Google Scholar] [CrossRef]
  25. Lu, L.; Weng, Q.; Xiao, D.; Guo, H.; Li, Q.; Hui, W. Spatiotemporal variation of surface urban heat islands in relation to land cover composition and configuration: A multi-scale case study of Xi’an, China. Remote Sens. 2020, 12, 2713. [Google Scholar] [CrossRef]
  26. Gao, J.; Gong, J.; Yang, J.; Li, J.; Li, S. Measuring Spatial Connectivity between patches of the heat source and sink (SCSS): A new index to quantify the heterogeneity impacts of landscape patterns on land surface temperature. Landsc. Urban Plann. 2022, 217, 104260. [Google Scholar] [CrossRef]
  27. Zhou, W.; Qian, Y.; Li, X.; Li, W.; Han, L. Relationships between land cover and the surface urban heat island: Seasonal variability and effects of spatial and thematic resolution of land cover data on predicting land surface temperatures. Landsc. Ecol. 2014, 29, 153–167. [Google Scholar] [CrossRef]
  28. Kuang, W.; Yang, T.; Liu, A.; Zhang, C.; Lu, D.; Chi, W. An EcoCity model for regulating urban land cover structure and thermal environment: Taking Beijing as an example. Sci. China Earth Sci. 2017, 60, 1098–1109. [Google Scholar] [CrossRef]
  29. Farhadi, H.; Faizi, M.; Sanaieian, H. Mitigating the urban heat island in a residential area in Tehran: Investigating the role of vegetation, materials, and orientation of buildings. Sustain. Cities 2019, 46, 101448. [Google Scholar] [CrossRef]
  30. Wang, X.; Meng, Q.; Zhang, L.; Hu, D. Evaluation of urban green space in terms of thermal environmental benefits using geographical detector analysis. Int. J. Appl. Earth Obs. Geoinf. 2021, 105, 102610. [Google Scholar] [CrossRef]
  31. Peng, J.; Jia, J.; Liu, Y.; Li, H.; Wu, J. Seasonal contrast of the dominant factors for spatial distribution of land surface temperature in urban areas. Remote Sens. Environ. 2018, 215, 255–267. [Google Scholar] [CrossRef]
  32. Mushore, T.D.; Mutanga, O.; Odindi, J. Estimating urban LST using multiple remotely sensed spectral indices and elevation retrievals. Sustain. Cities Soc. 2022, 78, 103623. [Google Scholar] [CrossRef]
  33. Liu, W.; Meng, Q.; Allam, M.; Zhang, L.; Hu, D.; Menenti, M. Driving factors of land surface temperature in urban agglomerations: A case study in the pearl river delta, china. Remote Sens. 2021, 13, 2858. [Google Scholar] [CrossRef]
  34. Taripanah, F.; Ranjbar, A. Quantitative analysis of spatial distribution of land surface temperature (LST) in relation Ecohydrological, terrain and socio-economic factors based on Landsat data in mountainous area. Adv. Space Res. 2021, 68, 3622–3640. [Google Scholar] [CrossRef]
  35. Ma, Y.; Zhang, S.; Yang, K.; Li, M. Influence of spatiotemporal pattern changes of impervious surface of urban megaregion on thermal environment: A case study of the Guangdong–Hong Kong–Macao Greater Bay Area of China. Ecol. Indic. 2021, 121, 107106. [Google Scholar] [CrossRef]
  36. Yu, B.; Tang, M.; Wu, Q.; Yang, C.; Deng, S.; Shi, K.; Peng, C.; Wu, J.; Chen, Z. Urban built-up area extraction from log-transformed NPP-VIIRS nighttime light composite data. IEEE Geosci. Remote Sens. Lett. 2018, 15, 1279–1283. [Google Scholar] [CrossRef]
  37. Wang, Z.; Sun, D.; Hu, C.; Wang, Y.; Zhang, J. Seasonal contrast and interactive effects of potential drivers on land surface temperature in the Sichuan basin, China. Remote Sens. 2022, 14, 1292. [Google Scholar] [CrossRef]
  38. Cao, C.; Lee, X.; Liu, S.; Schultz, N.; Xiao, W.; Zhang, M.; Zhao, L. Urban heat islands in China enhanced by haze pollution. Nat. Commun. 2016, 7, 12509. [Google Scholar] [CrossRef]
  39. Huang, X.; Wang, Y. Investigating the effects of 3D urban morphology on the surface urban heat island effect in urban functional zones by using high-resolution remote sensing data: A case study of Wuhan, Central China. ISPRS J. Photogramm. Remote Sens. 2019, 152, 119–131. [Google Scholar] [CrossRef]
  40. Xiong, Y.; Zhang, F. Effect of human settlements on urban thermal environment and factor analysis based on multi-source data: A case study of Changsha city. J. Geogr. Sci. 2021, 31, 819–838. [Google Scholar] [CrossRef]
  41. Wang, J.F.; Li, X.H.; Christakos, G.; Liao, Y.L.; Zhang, T.; Gu, X.; Zheng, X.Y. Geographical detectors-based health risk assessment and its application in the neural tube defects study of the Heshun Region, China. Int. J. Geogr. Inf Sci. 2010, 24, 107–127. [Google Scholar] [CrossRef]
  42. Dong, Y.; Yin, D.; Li, X.; Huang, J.; Su, W.; Li, X.; Wang, H. Spatial–temporal evolution of vegetation NDVI in association with climatic, environmental and anthropogenic factors in the loess plateau, China during 2000–2015: Quantitative analysis based on geographical detector model. Remote Sens. 2021, 13, 4380. [Google Scholar] [CrossRef]
  43. Gao, Q.; Fang, C.; Liu, H.; Zhang, L. Conjugate evaluation of sustainable carrying capacity of urban agglomeration and multi-scenario policy regulation. Sci. Total Environ. 2021, 785, 147373. [Google Scholar] [CrossRef]
  44. Fang, C.; Gao, Q.; Zhang, X.; Cheng, W. Spatiotemporal characteristics of the expansion of an urban agglomeration and its effect on the eco-environment: Case study on the northern slope of the Tianshan Mountains. Sci. China Earth Sci. 2019, 62, 1461–1472. [Google Scholar] [CrossRef]
  45. Wei, Y.; Shi, Z.; Biswas, A.; Yang, S.; Ding, J.; Wang, F. Updated information on soil salinity in a typical oasis agroecosystem and desert-oasis ecotone: Case study conducted along the Tarim River, China. Sci. Total Environ. 2020, 716, 135387. [Google Scholar] [CrossRef]
  46. Wang, H.; Cai, Y.; Yang, Q.; Gong, Y.; Lv, G. Factors that alter the relative importance of abiotic and biotic drivers on the fertile island in a desert-oasis ecotone. Sci. Total Environ. 2019, 697, 134096. [Google Scholar] [CrossRef]
  47. Yan, Y.; Chai, Z.; Zibibula, S.; Yang, S. The temporal and spatial changes of the ecological environment quality of the urban agglomeration on the northern slope of Tianshan Mountain and the influencing factors. Ecol. Indic. 2021, 133, 108380. [Google Scholar]
  48. Li, C.; Wang, R.; Cui, X.; Wu, F.; Yan, Y.; Peng, Q.; Qian, Z.; Xu, Y. Responses of vegetation spring phenology to climatic factors in Xinjiang, China. Ecol. Indic. 2021, 124, 107286. [Google Scholar] [CrossRef]
  49. Ouyang, X.; Tang, L.; Wei, X.; Li, Y. Spatial interaction between urbanization and ecosystem services in Chinese urban agglomerations. Land Use Policy 2021, 109, 105587. [Google Scholar] [CrossRef]
  50. Wan, Z. New refinements and validation of the collection-6 MODIS land-surface temperature/emissivity product. Remote Sens. Environ. 2014, 140, 36–45. [Google Scholar] [CrossRef]
  51. Qiao, Z.; Tian, G.; Xiao, L. Diurnal and seasonal impacts of urbanization on the urban thermal environment: A case study of Beijing using MODIS data. ISPRS J. Photogramm. Remote Sens. 2013, 85, 93–101. [Google Scholar] [CrossRef]
  52. Wei, B.; Bao, Y.; Yu, S.; Yin, S.; Zhang, Y. Analysis of land surface temperature variation based on MODIS data a case study of the agricultural pastural ecotone of northern China. Int. J. Appl. Earth Obs. Geoinf. 2021, 100, 102342. [Google Scholar] [CrossRef]
  53. Song, Z.; Yang, H.; Huang, X.; Yu, W.; Huang, J.; Ma, M. The spatiotemporal pattern and influencing factors of land surface temperature change in China from 2003 to 2019. Int. J. Appl. Earth Obs. Geoinf. 2021, 104, 102537. [Google Scholar] [CrossRef]
  54. Yang, J.; Huang, X. The 30 m annual land cover dataset and its dynamics in China from 1990 to 2019. Earth Syst. Sci. Data 2021, 13, 3907–3925. [Google Scholar] [CrossRef]
  55. Yu, H.; Bian, Z.; Mu, S.; Yuan, J.; Chen, F. Effects of climate change on land cover change and vegetation dynamics in Xinjiang, China. Int. J. Environ. Res. Public Health 2020, 17, 4865. [Google Scholar] [CrossRef]
  56. Mahmood, R.; Jia, S. Assessment of hydro-climatic trends and causes of dramatically declining stream flow to Lake Chad, Africa, using a hydrological approach. Sci. Total Environ. 2019, 675, 122–140. [Google Scholar] [CrossRef]
  57. Mohammad, P.; Goswami, A. Temperature and precipitation trend over 139 major Indian cities: An assessment over a century. Model. Earth Syst. Environ. 2019, 5, 1481–1493. [Google Scholar] [CrossRef]
  58. Hurst, H.E. Long-term storage capacity of reservoirs. Trans. Am. Soc. Civ. Eng. 1951, 116, 770–799. [Google Scholar] [CrossRef]
  59. Mandelbrot, B.B.; Wallis, J.R. Robustness of the rescaled range R/S in the measurement of noncyclic long run statistical dependence. Resour. Res. 1969, 5, 967–988. [Google Scholar] [CrossRef]
  60. Yang, L.; Fu, R.; He, W.; He, Q.; Liu, Y. Adaptive thermal comfort and climate responsive building design strategies in dry–hot and dry–cold areas: Case study in Turpan, China. Energy Build. 2020, 209, 109678. [Google Scholar] [CrossRef]
  61. Pei, H.; Fang, S.; Lin, L.; Qin, Z.; Wang, X. Methods and applications for ecological vulnerability evaluation in a hyper-arid oasis: A case study of the Turpan Oasis, China. Environ. Earth Sci. 2015, 74, 1449–1461. [Google Scholar] [CrossRef] [Green Version]
  62. Oke, T.R. Boundary Layer Climates; Routledge: Abingdon-on-Thames, UK, 2002. [Google Scholar]
  63. Hao, X.; Li, W.; Deng, H. The oasis effect and summer temperature rise in arid regions-case study in Tarim Basin. Sci. Rep. 2016, 6, 35418. [Google Scholar] [CrossRef] [PubMed]
  64. Hao, X.; Li, W. Oasis cold island effect and its influence on air temperature: A case study of Tarim Basin, Northwest China. J. Arid Land 2016, 8, 172–183. [Google Scholar] [CrossRef]
  65. Han, S.; Yang, Z. Cooling effect of agricultural irrigation over Xinjiang, Northwest China from 1959 to 2006. Environ. Res. Lett. 2013, 8, 024039. [Google Scholar] [CrossRef]
  66. Zhu, X.; Liang, S.; Pan, Y. Observational evidence of the cooling effect of agricultural irrigation in Jilin, China. Clim. Change 2012, 114, 799–811. [Google Scholar] [CrossRef]
  67. Wen, K.; Ren, G.; Li, J.; Zhang, A.; Ren, Y.; Sun, X.; Zhou, Y. Recent surface air temperature change over mainland China based on an urbanization-bias adjusted dataset. J. Clim. 2019, 32, 2691–2705. [Google Scholar] [CrossRef]
  68. Solomon, S.; Rosenlof, K.H.; Portmann, R.W.; Daniel, J.S.; Davis, S.M.; Sanford, T.J.; Plattner, G.-K. Contributions of stratospheric water vapor to decadal changes in the rate of global warming. Science 2010, 327, 1219–1223. [Google Scholar] [CrossRef] [PubMed]
  69. Zhang, Y.; He, Y.; Li, Y.; Jia, L. Spatiotemporal variation and driving forces of NDVI from 1982 to 2015 in the Qinba Mountains, China. Environ. Sci. Pollut. Res. 2022, 29, 52277–52288. [Google Scholar] [CrossRef]
  70. Huang, C.; Yang, Q.; Zhang, H. Temporal and spatial variation of ndvi and its driving factors in qinling mountain. Water 2021, 13, 3154. [Google Scholar] [CrossRef]
  71. Ma, X.; Huang, G.; Cao, J. The significant roles of anthropogenic aerosols on surface temperature under carbon neutrality. Sci. Bull. 2022, 67, 470–473. [Google Scholar] [CrossRef]
  72. Choobari, O.A.; Zawar-Reza, P.; Sturman, A. The global distribution of mineral dust and its impacts on the climate system: A review. Atmos. Res. 2014, 138, 152–165. [Google Scholar] [CrossRef]
  73. Liu, J.; Ding, J.; Rexiding, M.; Li, X.; Zhang, J.; Ran, S.; Bao, Q.; Ge, X. Characteristics of dust aerosols and identification of dust sources in Xinjiang, China. Atmos. Environ. 2021, 262, 118651. [Google Scholar] [CrossRef]
  74. Lan, Y.; Zhan, Q.J.B. How do urban buildings impact summer air temperature? The effects of building configurations in space and time. Build. Environ. 2017, 125, 88–98. [Google Scholar] [CrossRef]
  75. Zhou, D.; Zhang, L.; Li, D.; Huang, D.; Zhu, C. Climate–vegetation control on the diurnal and seasonal variations of surface urban heat islands in China. Environ. Res. Lett. 2016, 11, 074009. [Google Scholar] [CrossRef]
  76. Yin, Z.; Wu, P.; Foody, G.M.; Wu, Y.; Liu, Z.; Du, Y.; Ling, F. Spatiotemporal fusion of land surface temperature based on a convolutional neural network. IEEE Trans. Geosci. Remote Sens. 2020, 59, 1808–1822. [Google Scholar] [CrossRef]
  77. Hutengs, C.; Vohland, M. Downscaling land surface temperatures at regional scales with random forest regression. Remote Sens. Environ. 2016, 178, 127–141. [Google Scholar] [CrossRef]
  78. Zhang, Y.; Liu, J.; Wen, Z. Predicting surface urban heat island in Meihekou city, China: A combination method of monte carlo and random forest. Chin. Geogr. Sci. 2021, 31, 659–670. [Google Scholar] [CrossRef]
  79. Chen, D.; Zhang, F.; Zhang, M.; Meng, Q.; Jim, C.Y.; Shi, J.; Tan, M.L.; Ma, X. Landscape and vegetation traits of urban green space can predict local surface temperature. Sci. Total Environ. 2022, 825, 154006. [Google Scholar] [CrossRef]
  80. Zhao, D.; Wang, J.; Zhao, X.; Triantafilis, J. Clay content mapping and uncertainty estimation using weighted model averaging. Catena 2022, 209, 105791. [Google Scholar] [CrossRef]
Figure 1. Location of the urban agglomeration on the northern slope of Tianshan Mountains.
Figure 1. Location of the urban agglomeration on the northern slope of Tianshan Mountains.
Sustainability 14 10663 g001
Figure 2. The research process of this study.
Figure 2. The research process of this study.
Sustainability 14 10663 g002
Figure 3. Spatial distribution of LST classes in the UANSTM from 2005 to 2021. Note: the numbers (1)–(9) indicate the city labels, and the specific names of the cities can be found in Figure 1.
Figure 3. Spatial distribution of LST classes in the UANSTM from 2005 to 2021. Note: the numbers (1)–(9) indicate the city labels, and the specific names of the cities can be found in Figure 1.
Sustainability 14 10663 g003
Figure 4. The spatial distribution of AV from 2005 to 2021. Note: the numbers (1)–(9) indicate the city labels, and the specific names of the cities can be found in Figure 1.
Figure 4. The spatial distribution of AV from 2005 to 2021. Note: the numbers (1)–(9) indicate the city labels, and the specific names of the cities can be found in Figure 1.
Sustainability 14 10663 g004
Figure 5. Spatial distribution of slope value from 2005 to 2021. Note: the numbers (1)–(9) indicate the city labels, and the specific names of the cities can be found in Figure 1.
Figure 5. Spatial distribution of slope value from 2005 to 2021. Note: the numbers (1)–(9) indicate the city labels, and the specific names of the cities can be found in Figure 1.
Sustainability 14 10663 g005
Figure 6. The changing trend levels of LST from 2005 to 2021. Note: the numbers (1)–(9) indicate the city labels, and the specific names of the cities can be found in Figure 1.
Figure 6. The changing trend levels of LST from 2005 to 2021. Note: the numbers (1)–(9) indicate the city labels, and the specific names of the cities can be found in Figure 1.
Sustainability 14 10663 g006
Figure 7. Spatial distribution of the Hurst index. Note: the numbers (1)–(9) indicate the city labels, and the specific names of the cities can be found in Figure 1.
Figure 7. Spatial distribution of the Hurst index. Note: the numbers (1)–(9) indicate the city labels, and the specific names of the cities can be found in Figure 1.
Sustainability 14 10663 g007
Figure 8. The future changing trend levels of LST. Note: the numbers (1)–(9) indicate the city labels, and the specific names of the cities can be found in Figure 1.
Figure 8. The future changing trend levels of LST. Note: the numbers (1)–(9) indicate the city labels, and the specific names of the cities can be found in Figure 1.
Sustainability 14 10663 g008
Figure 9. The q value of each influencing factor in the daytime and the nighttime.
Figure 9. The q value of each influencing factor in the daytime and the nighttime.
Sustainability 14 10663 g009
Figure 10. Interaction of each influencing factor in the daytime and the nighttime. Note: the dotted box represents nonlinear enhancement, and the rest are bivariate enhancement.
Figure 10. Interaction of each influencing factor in the daytime and the nighttime. Note: the dotted box represents nonlinear enhancement, and the rest are bivariate enhancement.
Sustainability 14 10663 g010
Figure 11. The risk interval or category of each influencing factor in the daytime and the nighttime. Note: the red indicates the highest average LST in the interval and the blue indicates the lowest.
Figure 11. The risk interval or category of each influencing factor in the daytime and the nighttime. Note: the red indicates the highest average LST in the interval and the blue indicates the lowest.
Sustainability 14 10663 g011
Figure 12. Differential judgment of the influence of each influencing factor on LST. Note: the yellow Y marks indicate significant differences in the effect on LST, while green N marks indicate no significant differences.
Figure 12. Differential judgment of the influence of each influencing factor on LST. Note: the yellow Y marks indicate significant differences in the effect on LST, while green N marks indicate no significant differences.
Sustainability 14 10663 g012
Table 1. Details of the annual mean LST images of the UANSTM from 2005 to 2021.
Table 1. Details of the annual mean LST images of the UANSTM from 2005 to 2021.
TimeMaximumMinimumAverageStandard DeviationTimeMaximumMinimumAverageStandard Deviation
2005 daytime46.50−13.1023.258.822013 nighttime10.57−25.12−0.254.09
2005 nighttime12.64−27.38−1.284.342014 daytime41.36−13.5223.148.38
2006 daytime46.19−12.3524.138.402014 nighttime9.59−25.68−1.414.18
2006 nighttime12.85−26.75−0.454.232015 daytime43.60−12.2524.318.67
2007 daytime41.75−12.8023.808.482015 nighttime10.92−22.720.194.16
2007 nighttime10.99−27.53−0.614.202016 daytime44.01−12.1423.598.79
2008 daytime42.70−12.8324.668.462016 nighttime12.72−28.190.034.21
2008 nighttime10.53−29.97−0.724.232017 daytime43.37−12.0824.298.77
2009 daytime41.92−13.0024.208.672017 nighttime13.84−28.11−0.044.38
2009 nighttime10.17−26.78−1.194.162018 daytime41.88−13.3823.038.55
2010 daytime42.46−13.2422.308.572018 nighttime9.87−29.60−1.674.07
2010 nighttime9.77−27.17−1.314.092019 daytime42.60−13.8124.098.48
2011 daytime44.43−12.8723.158.742019 nightime11.10−25.49−0.384.26
2011 nighttime9.86−27.21−1.384.202020 daytime43.61−12.3024.458.41
2012 daytime44.03−14.0923.078.872020 nighttime11.04−27.48−0.964.22
2012 nighttime10.37−28.24−1.764.302021 daytime43.25−12.7324.638.65
2013 daytime46.35−12.5424.638.602021 nighttime11.09−24.54−0.754.18
Note: the unit in Table 1 is °C.
Table 2. The details of influencing factors’ data.
Table 2. The details of influencing factors’ data.
Types of Influencing FactorsFactorsInitial Resolution/Resampling ResolutionTimeAccess
Land cover factorsLULC1000 m/1000 m2019https://zenodo.org/, accessed on 19 April 2022
NDVI1000 m/1000 m2019https://ladsweb.modaps.eosdis.nasa.gov/, accessed on 15 April 2022
NDBSI500 m/1000 m2019https://ladsweb.modaps.eosdis.nasa.gov/, accessed on 15 April 2022
Climatic factorsPrecipitation1000 m/1000 m2019http://data.tpdc.ac.cn/, accessed on 16 April 2022
AOD1000 m/1000 m2019https://ladsweb.modaps.eosdis.nasa.gov/, accessed on 15 April 2022
Socioeconomic factorsNTL500 m/1000 m2019https://eogdata.mines.edu/products/vnl/, accessed on 15 April 2022
GDP1000 m/1000 m2019https://www.resdc.cn/, accessed on 15 April 2022
POP1000 m/1000 m2019https://www.resdc.cn/, accessed on 15 April 2022
Topographical factorsDEM30 m/1000 m-http://www.gscloud.cn/, accessed on 16 April 2022
Gradient30 m/1000 m-http://www.gscloud.cn/, accessed on 16 April 2022
Table 3. The classification standard of land surface temperature.
Table 3. The classification standard of land surface temperature.
LST ClassesLST Range
EHTT > μ + 1.5 std
HTμ + 0.5 std < Tμ + 1.5 std
MTμ − 0.5 stdT < μ + 0.5 std
LTμ − 1.5 stdT < μ − 0.5 std
ELTT < μ − 1.5 std
Note: T is the LST value, μ is the mean, and std is the standard deviation.
Table 4. Condition and level definition of change trend.
Table 4. Condition and level definition of change trend.
Slope and Z ValueChanging TrendsSlope, Z Value and Hurst IndexFuture Changing Trends
Slope < 0, Z ≤ −2.58Very significantly decreasedSlope < 0, Z ≤ −2.58, H > 0.5Sustained very significant decrease
Slope < 0, −2.58 < Z ≤ −1.96Significantly decreasedSlope < 0, −2.58 < Z ≤ −1.96, H > 0.5Sustained significant decrease
Slope < 0, −1.96 < Z < 1.96Not significantly decreasedSlope < 0, −1.96 < Z < 1.96, H > 0.5Sustained no significant decrease
Slope > 0, 1.96 < Z < 1.96Not significantly increasedSlope > 0, 1.96 < Z < 1.96, H > 0.5Sustained no significant increase
Slope > 0, 1.96 ≤ Z < 2.58Significantly increasedSlope > 0, 1.96 ≤ Z < 2.58, H > 0.5Sustained significant increase
Slope > 0, Z ≥ 2.58Very significantly increasedSlope > 0, Z ≥ 2.58, H > 0.5Sustained very significant increase
H ≤ 0.5Unknown
Table 5. Types of interactions and judgment criteria.
Table 5. Types of interactions and judgment criteria.
Interaction TypesJudgment Standard
Weaken, nonlinearq(X1∩X2) < Min(q(X1), q(X2))
Weaken, univariateMin(q(X1), q(X2)) < q(X1∩X2) < Max(q(X1), q(X2))
Independentq(X1∩X2) = q(X1) + q(X2)
Enhance bivariateq(X1∩X2) > Max(q(X1), q(X2))
Enhance, nonlinearq(X1∩X2) > q(X1) + q(X2)
Note: Min(q(X1), q(X2)) is q(X1), q(X2), both to take the minimum. Max(q(X1), q(X2)) is q(X1), q(X2), both to take the maximum. q(X1) + q(X2) is q(X1), q(X2), both to sum.
Table 6. The optimal discretization method and number of discretization intervals for continuous variables.
Table 6. The optimal discretization method and number of discretization intervals for continuous variables.
Continuous VariablesDaytimeNighttime
MethodNumber of IntervalsMethodNumber of Intervals
NDVIstandard deviation7standard deviation7
NDBSInatural breaks7standard deviation5
Precipitationnatural breaks6natural breaks6
AODstandard deviation7standard deviation7
NTLquantile7quantile7
GDPquantile6quantile7
POPquantile7quantile7
DEMstandard deviation7standard deviation5
Gradientstandard deviation7standard deviation7
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Liang, H.; Kasimu, A.; Ma, H.; Zhao, Y.; Zhang, X.; Wei, B. Exploring the Variations and Influencing Factors of Land Surface Temperature in the Urban Agglomeration on the Northern Slope of the Tianshan Mountains. Sustainability 2022, 14, 10663. https://0-doi-org.brum.beds.ac.uk/10.3390/su141710663

AMA Style

Liang H, Kasimu A, Ma H, Zhao Y, Zhang X, Wei B. Exploring the Variations and Influencing Factors of Land Surface Temperature in the Urban Agglomeration on the Northern Slope of the Tianshan Mountains. Sustainability. 2022; 14(17):10663. https://0-doi-org.brum.beds.ac.uk/10.3390/su141710663

Chicago/Turabian Style

Liang, Hongwu, Alimujiang Kasimu, Haitao Ma, Yongyu Zhao, Xueling Zhang, and Bohao Wei. 2022. "Exploring the Variations and Influencing Factors of Land Surface Temperature in the Urban Agglomeration on the Northern Slope of the Tianshan Mountains" Sustainability 14, no. 17: 10663. https://0-doi-org.brum.beds.ac.uk/10.3390/su141710663

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop