Land surface temperature (LST) has a critical impact on the energy balance of land surface processes and ecosystem stability. Meanwhile, LST is controlled by multiple factors at the surface, resulting in heterogeneity of its spatial distribution. To understand the drivers of LST spatial heterogeneity and their contributions, the effects of air temperature, normalized difference vegetation index (NDVI), soil moisture, net surface radiation, precipitation, aerosol optical depth (AOD), evapotranspiration, water vapor, digital elevation model (DEM), climate type, and land cover type on LST spatial heterogeneity was analyzed in this study with GeoDetector. The results showed that the explanatory ability of air temperature to impact the spatial heterogeneity of LST was the largest in each year with a mean value of 0.74, followed by water vapor with a mean value of 0.7, and the driving effect of the factors on LST showed an increasing trend year by year. However, the land cover type did not have an effect on the spatial heterogeneity of LST for the univariate analysis in this study. In addition, the interaction analysis indicated that the spatial distribution of LST was jointly driven by all the driving factors. Among them, air temperature had the strongest interaction with other factors, with the strength of the effect in the range of 0.73–0.8. In terms of the highly sensitive area of LST for each driver, AOD has the largest driving area, accounting for 15.8% of the total area, followed by WV, TA, and ET at about 11%, and the remaining variables are less than 10%. During the study period, the area of the highly sensitive region of LST for each factor showed an overall decreasing trend, indicating that the influence of the driving factors on LST will be stronger and more concentrated. Generally, this study provides meaningful understanding of the spatial heterogeneity of LST since 2003 and provides a scientific reference for coping with climate change, analyzing surface environmental patterns, and protecting ecological environment.

The relationship between land surface temperature (LST) and environmental factors is complex and nonlinear. To determine these relationships for China, this study analyzed the driving effects of air temperature, vegetation index, soil moisture, net surface radiation, precipitation, aerosols, evapotranspiration, and water vapor on LST based on remote-sensing and reanalysis data from 2003–2018, using a convergent cross-mapping method. During the study period, air temperature and net surface radiation were the dominant drivers of LST with a cross-mapping skill above 0.9. Vegetation index and evapotranspiration were the secondary drivers of LST with a cross-mapping skill that was higher than 0.5. Except for air temperature and net surface radiation, the direction and strength of the effects of the driving factors on LST were related to the climate type. The effects of air temperature and net radiation on LST diminished from north to south, indicating that LST was more sensitive to air temperature and net radiation in energy-limited regions. However, the effects of vegetation index and evapotranspiration on LST varied significantly across climate zones; that is, positive effects were mostly in non-monsoonal zones and negative effects were primarily in monsoonal zones. Our results quantified the driving role of environmental factors on LST and provided a comprehensive understanding of LST dynamics.

The difference (DIF) between land surface temperature (Ts) and near surface air temperature (Ta) is the key indicator of the energy budget of the land surface, which has a more complex process than the individual Ts or Ta. However, the spatiotemporal variations and influencing factors of DIF remain incomplete. The contribution of vegetation and soil moisture (SM) as key driving factors to DIF is not yet clear. Here, we analyzed the spatiotemporal variation patterns of DIF in China from 2011 to 2023 using in situ Ts and Ta data. A convergent cross-mapping method was employed to explore the causal relationship between SM, NDVI and DIF, and subsequently calculated the contribution of NDVI and SM variations to DIF under different climatic backgrounds. The results indicate that during the study period, DIF values were all above 0 °C and showed a significant increasing trend with a national mean slope of 0.02 °C/a. In general, vegetation and SM have a driving effect on DIF, with vegetation contributing more to DIF (0.11) than SM (0.08) under different surface properties. The background values of SM and temperature have a significant effect on the spatial and temporal distribution of DIF, as well as the correlation of vegetation and soil moisture to DIF. The study outcomes contribute to a better understanding of the coupling relationship between the land surface and atmosphere, which are also crucial for addressing climate change and ecological environmental management.