Land Surface Temperature Research Articles

Green space-building integration for Urban Heat Island mitigation: Insights from Beijing's fifth ring road district

In this research, we delve into the complex arrangement of urban landscapes, where green spaces and buildings are not merely co-existing but are interwoven into a cohesive fabric that shapes the thermal environment. Our approach transcends the conventional methods of analysis, which typically isolate the roles of greenery or built environments. Instead, we adopt a synergistic perspective that recognizes the collective influence of these landscape constituents on the urban thermal pattern. Key insights are: (1) A linear decrease in average land surface temperature with increasing green space coverage is observed. However, substantial temperature variations (up to 8 °C) within the same coverage interval highlight the significant impact of built-up pattern on thermal conditions; (2) High Building Height and Floor Area Ratio, and low Building Coverage Ratio and Sky View Factor, are linked to cooler temperatures in areas with up to 50 % green space; (3) The study suggests that low-temperature areas can inform the adjustment of built-up patterns in high-temperature areas, offering a strategy for thermal environment optimization within specific green space coverage intervals. This research contributes insights into the integrated planning of green spaces and buildings, with implications for urban development and renewal initiatives aiming to enhance the urban thermal environment.

High resolution (1-km) surface soil moisture generation from SMAP SSM by considering its difference between freezing and thawing periods in the source region of the Yellow River

Soil moisture (SM) is a critical component of the land surface hydrological cycle, significantly impacting various sectors such as hydrology, meteorology, and agriculture. Accurate, high-resolution SM data are essential for effective flood forecasting, water resource management, and understanding the soil freeze-thaw processes in cold regions. This study aims to generate 1 km resolution liquid surface SM (SSM) data with a twice-daily update frequency by downscaling SMAP Level-4 SSM data using random forest (RF) and multiple linear regression (MLR) in the source region of the Yellow River (SRYR), by considering the differences in SM changes between freezing and thawing periods. To obtain the SSM data, 16 downscaling schemes of both RF and MLR were designed for each of the three scenarios. In each downscaling process, both land surface temperature (LST) and normalized difference vegetation index (NDVI) were utilized in MLR and RF models, alongside various combinations of additional variables such as albedo, elevation, leaf area index (LAI), soil texture. Results showed that during the freezing period, RF produced superior SSM estimates when supplemented with NDVI, LST, albedo, elevation, LAI, and soil texture. MLR was more effective during the thawing period when paired with NDVI, LST, elevation, LAI, and soil texture. During the freezing period, the downscaled SMAP SSM exhibited average R, RMSE, ubRMSE of 0.76, 0.029 m3·m-3, and 0.023 m3·m-3, respectively, when compared with in-situ measurements. During the thawing period, the average R, MAE, RMSE, and ubRMSE between the downscaled SMAP SSM and in-situ measurements were 0.52, 0.057 m3·m-3, 0.067 m3·m-3, and 0.054 m3·m-3, respectively, compared to 0.45, 0.070 m3·m-3, 0.083 m3·m-3, and 0.060 m3·m-3 for the original SMAP SSM. Thus, the research significantly enhances both the accuracy and spatial resolution of SMAP SSM estimations, underscoring its vital role in advancing hydrological studies within the SRYR.

Spatiotemporal dynamics of vegetation net primary productivity and its response to climate variability

Despite its significance, net primary productivity is rarely used as an indicator of climate change. The purpose of this study was to assess the spatiotemporal dynamics of vegetation net primary productivity and its response to climate variability in the Welmel watershed, southeastern Ethiopia. The investigation utilized data from the Moderate Resolution Imaging Spectroradiometer, and the Climate Hazards Group Infrared Precipitation with Stations accessed through the Google Earth Engine cloud computing platform. To evaluate the trend and response of net primary productivity to climatic factors, the Mann-Kendall and Theil Sen’s tests, the Pearson correlation coefficient, and stability analysis through the coefficient of variation of net primary productivity were utilized. The analysis of the temporal trend of net primary productivity revealed that the entire watershed, forest, wood, and cultivated land areas exhibited decreasing trends with net primary productivity values of -0.681, -4.361, -1.41, and − 3.25 g C m− 2/year, respectively. On the other hand, shrubs and grazing land displayed increasing trends of 1.15 and 2.04 g C m− 2/year, respectively. The wood and shrubland trends were statistically significant (p < 0.05). In the spatial analysis, an area of 64.6% showed a decreasing trend, whereas 33.78% displayed an increasing trend in net primary productivity values. The areas with stable and unstable variations in net primary productivity accounted for 26.47% and 8.02%, respectively. The net primary productivity of the land cover types and entire watershed were positively correlated with rainfall, however, only forests and woodlands were statistically significantly correlated (p < 0.05). The net primary productivity and land surface temperature were negatively correlated, except for forests and cultivated land, which were positively correlated. The results of this study provide essential information for managing ecosystems, supporting agricultural practices, and enhancing community resilience in a changing climate.

Delineating Urbanicity and Rurality: Impact on Environmental Exposure Assessment.

Environmental exposures and their health impacts can vary substantially between urban and rural areas. However, different methods for classifying these areas could lead to inconsistencies in environmental exposure and health studies, which are often overlooked. We constructed different urban/rural classification systems based on multiple population-based (e.g., total population, population density, and commuting) and built-environment-based (e.g., nighttime light intensity, building density, road density, distance to urban centers, point of interest density, and urban area coverage) indicators and various classification schemes. These classification systems were applied to Virginia and West Virginia, United States. We compared differences in urban/rural spatial patterns, demographic compositions, and exposures of particulate matter (PM2.5), greenspace, and land surface temperature using these urban/rural classification systems to understand their impacts on environmental exposure and health research. Our findings reveal clear differences in spatial patterns and demographic compositions across various systems. We also observed that different systems can lead to changes in the magnitude and direction of urban/rural disparities in environmental exposure assessment. Addressing the complexities in delineating urbanicity and rurality may include careful consideration of classification systems to reflect those aspects of urbanicity and rurality that are relevant to the research question or the use of multiple, complementary systems.

Identification of the Driving factors impacts of Land Surface Albedo over Iran: An analysis with the MODIS data

Albedo is a key parameter in climatic research and depends on environmental and climatic factors. Modeling these factors greatly contributes to understanding environmental variations. To this end, the data of Land Surface Albedo, Land Surface Temperature (LST), Vegetation, Snow, Elevation, Slope, and Aspect of the MODIS sensor from 1/1/2001 to 30/12/2021 with a 1000-m spatial resolution were used. After pre-processing, monthly, seasonal, and annual albedo modeling was performed using multiple linear regression (MLR) in the highlands of Iran. The results of monthly modeling revealed the salient direct role of snow on the albedo of Iran's highlands in all months, except for July, August, and September. In these months, due to the lack of snow coverage and the fruiting of agricultural lands and gardens, the inverse role of vegetation on albedo variations is determining. Seasonal examinations also showed that snow plays a significant role on the albedo of Iran's highlands in winter, spring, and fall; however, vegetation has a determining role in the summer. The annual results indicated that snow, vegetation, elevation, slope, LST, and aspect, respectively, are the factors affecting albedo in the highlands of Iran. Furthermore, the role of snow, LST, and aspect is positive, while the role of vegetation, elevation, and slope is negative on albedo.

Analysis of the Biennial Productivity of Arabica Coffee with Google Earth Engine in the Northeast Region of São Paulo, Brazil

Numerous challenges are associated with the classification of satellite images of coffee plantations. The spectral similarity with other types of land use, variations in altitude, topography, production system (shaded and sun), and the change in spectral signature throughout the phenological cycle are examples that affect the process. This research investigates the influence of biennial Arabica coffee productivity on the accuracy of Landsat-8 image classification. The Google Earth Engine (GEE) platform and the Random Forest algorithm were used to process the annual and biennial mosaics of the Média Mogiana Region, São Paulo (Brazil), from 2017 to 2023. The parameters evaluated were the general hits of the seven classes of land use and coffee errors of commission and omission. It was found that the seasonality of the plant and its development phases were fundamental in the quality of coffee classification. The use of biennial mosaics, with Landsat-8 images, Brightness, Greenness, Wetness, SRTM data (elevation, aspect, slope), and LST data (Land Surface Temperature) also contributed to improving the process, generating a classification accuracy of 88.8% and reducing coffee omission errors to 22%.

Cooling benefits from urban planting depend on local background canopy cover: A continental cross-city comparison

Numerous studies have shown that the cooling efficiency (CE) of urban trees varies by cities with different climate backgrounds, and recent research further indicated that there may be large within-city variations in CE. However, how such within-city variations differ across cities, and their relations to the local percent of urban tree canopy (Ptree) remain poorly understood. This study aims to fill this gap based on a comparison study across 118 cities with different biomes and climates in the continental USA. We used the tree canopy layer of the National Land Cover Dataset (NLCD) 2011 to measure urban tree canopy (UTC), and calculated land surface temperature (LST) based on Landsat thermal bands. We found: 1) CE had larger within-city and cross-city spatial variations in cities located in arid and semi-arid biomes. 2) CE was related to Ptree in nonlinear ways in >90 % of the study cities. In most cities (approximately 70 %), CE had an L-shaped relationship with Ptree, showing that CE first declined quickly with the increase of Ptree, but then gradually dropped in a slower way or became relatively stable after Ptree reached a certain threshold. 3) While there was no significant difference in the types of CE-Ptree relationship among biomes and climates, the threshold of Ptree in CE-Ptree nonlinear ways was smaller in arid cities. The results of this threshold linking cooling benefits and current UTC can serve as a useful tool to prioritize locations for urban planting to maximize cooling benefits.

An integrated approach for analysing land surface temperature, heatwave and air pollution for sustainable planning in Delhi and its environs, India

An integrated approach for analysing land surface temperature, heatwave and air pollution for sustainable planning in Delhi and its environs, India

Forecasting malaria dynamics based on causal relations between control interventions, climatic factors, and disease incidence in western Kenya.

Malaria remains one of the deadliest diseases worldwide, especially among young children in sub-Saharan Africa. Predictive models are necessary for effective planning and resource allocation; however, statistical models suffer from association pitfalls. In this study, we used empirical dynamic modelling (EDM) to investigate causal links between climatic factors and intervention coverage with malaria for short-term forecasting. Based on data spanning the period from 2008 to 2022, we used convergent cross-mapping (CCM) to identify suitable lags for climatic drivers and investigate their effects, interaction strength, and suitability ranges on malaria incidence. Monthly malaria cases were collected at St. Elizabeth Lwak Mission Hospital. Intervention coverage and population movement data were obtained from household surveys in Asembo, western Kenya. Daytime land surface temperature (LSTD), rainfall, relative humidity (RH), wind speed, solar radiation, crop cover, and surface water coverage were extracted from remote sensing sources. Short-term forecasting of malaria incidence was performed using state-space reconstruction. We observed causal links between climatic drivers, bed net use, and malaria incidence. LSTD lagged over the previous month; rainfall and RH lagged over the previous two months; and wind speed in the current month had the highest predictive skills. Increases in LSTD, wind speed, and bed net use negatively affected incidence, while increases in rainfall and humidity had positive effects. Interaction strengths were more pronounced at temperature, rainfall, RH, wind speed, and bed net coverage ranges of 30-35°C, 30-120 mm, 67-80%, 0.5-0.7 m/s, and above 90%, respectively. Temperature and rainfall exceeding 35°C and 180 mm, respectively, had a greater negative effect. We also observed good short-term predictive performance using the multivariable forecasting model (Pearson correlation coefficient = 0.85, root mean square error = 0.15). Our findings demonstrate the utility of CCM in establishing causal linkages between malaria incidence and both climatic and non-climatic drivers. By identifying these causal links and suitability ranges, we provide valuable information for modelling the impact of future climate scenarios.

Assessment of Wind Energy Potential and Optimal Site Selection for Wind Energy Plant Installations in Igdir/Turkey

Wind energy is an eco-friendly, renewable, domestic, and infinite resource. These factors render the construction of wind turbines appealing to nations, prompting numerous governments to implement incentives to augment their installed capacity of wind turbines. Alongside augmenting the installed capacity of wind turbines, identifying suitable locations for their installation is crucial for optimizing turbine performance. This study aims to evaluate potential sites for wind power plant installation via a GIS, a mapping technique. The Analytic Hierarchy Process (AHP) was employed to assess the locations, including both quantitative and qualitative aspects that significantly impact the wind farm suitability map. Utilizing the GIS methodology, all datasets were examined through height and raster transformations of land surface temperature, plant density index, air pressure, humidity, wind speed, air temperature, land cover, solar radiation, aspect, slope, and topographical characteristics, resulting in the creation of a wind farm map. The correlation between the five-year meteorological data and environmental parameters (wind direction, daily wind speed, daily maximum and minimum air temperatures, daily relative humidity, daily average air temperature, solar radiation duration, daily cloud cover, air humidity, and air pressure) influencing the wind power plant in Iğdır province, including Iğdır Airport, Karakoyunlu, Aralık, and Tuzluca districts, was analyzed. If wind energy towers are installed at 1 km intervals across an area of roughly 858,180 hectares in Igdir province, an estimated 858,180 GWh of wind energy can be generated. The GIS-derived wind power plant map indicates that the installation sites for wind power plants are located in regions susceptible to wind erosion.

Natural forest regeneration is projected to reduce local temperatures

Forest regeneration is a crucial strategy for mitigating and adapting to global warming. Yet its precise impact on local climate remains uncertain, a factor that complicates decision-making when it comes to prioritizing investments. Here, we developed global maps illustrating how natural forest regeneration influences key local climate drivers—land surface temperature (LST), albedo, and evapotranspiration—using models fitted at a 1-km spatial resolution with a random forest classifier. We found that natural forest regeneration can alter annual mean LST by 0.01 °C, −0.59 °C, −0.50 °C, and −2.03 °C in Boreal, Mediterranean, Temperate, and Tropical regions, respectively. These variations underscore the region-specific effects of forest regeneration. Importantly, natural forest regeneration reduces LST across 64% of 1 billion hectares and 75% of 148 million hectares of potentially restorable land under different scenarios. These findings improve understanding of how forest regeneration can help regulate local climate, supporting climate adaptation efforts.

Assessing Ecological Environment in Dar es Salaam City: A Settlement Surface Ecological Index (SSEI) Analysis

Urbanization dramatically alters landscapes and ecological conditions. This study examines the environmental conditions within Dar es Salaam, Tanzania, using the Settlement Surface Ecological Index (SSEI). The SSEI combines biophysical parameters from satellite imagery to assess ecological health, revealing disparities across diverse urban typologies. This study’s methodology demonstrates advancements in using the Calculate Composite Index tool for streamlined index creation and greater flexibility in urban analysis. Statistical analysis reveals crucial relationships between ecological conditions and variables like tree cover, impervious surfaces, vegetation moisture, and Land Surface Temperature (LST). Multivariate clustering highlights distinct areas of high and Low SSEI, emphasizing disparities in urban ecological health. This study demonstrates how composite indices like the SSEI can inform targeted urban planning, promote green infrastructure, and mitigate the negative impacts of urbanization, ensuring both equitable and sustainable development in the face of climate change.

Post-fire vegetation dynamic patterns and drivers in Greater Hinggan Mountains: Insights from long-term remote sensing data analysis

Fire has become a major disturbing factor in boreal forests, and giant forest disturbances play a vital role in regulating the climate under global warming. Therefore, it is essential to investigate the spatiotemporal patterns and main drivers of post-fire vegetation recovery for forest ecological research and post-fire recovery management. However, previous studies have focused on the post-fire forest change within the entire fire perimeter, lacking separate analysis and comparison of the burned zone (BZ) and unburned zone (UNBZ). Here, we propose the utilization of Moderate Resolution Imaging Spectroradiometer land cover type and vegetation index data to monitor vegetation dynamics and explore its drivers after the most serious forest fire in the history of P.R. China in the Greater Hinggan Mountains (GHM). The temporal and spatial patterns of vegetation recovery in the BZ/UNBZ in the GHM were analyzed using the Sen & Mann-Kendall method, Hurst index and coefficient of variation, and their driving mechanisms were explored using GeoDetector and geographically weighted regression. The results showed that there were significant differences in the spatial distribution and fluctuation of vegetation between the BZ and UNBZ, and that the BZ exhibited higher productivity and vigor. Vegetation recovery was influenced by different dominant factors and changed over time, in which land surface temperature and precipitation dominated all the time, whereas topographic relief and elevation had a more significant contribution to vegetation recovery in the BZ and UNBZ, respectively. This study provides a scientific basis for the protection and management of vegetation in disturbed forested areas, particularly after fires.

Exploring the Configurational Relationships between Urban Heat Island Patterns and the Built Environment: A Case Study of Beijing

The spatial heterogeneity of land surface temperature (LST) within cities is profoundly influenced by the built environment. Although significant progress has been made in the study of the urban thermal environment, there is still a lack of research on how the pattern and structural layout of the built environment affects the thermal environment. In this study, we take the Fifth Ring Road of Beijing as an example, invert the urban LST on the basis of multisource spatial data, characterize the built environment, and use k-means cluster analysis to investigate the main influencing factors of the LST in different functional areas and building patterns within the city, as well as the spatial relationship between the built environment and the urban LST. The results show the following: (1) The urban heat island (UHI) effect occurs to varying degrees over a large part of the study area, and these UHI areas are mainly concentrated in the southwestern part of the city, forming a large contiguous area between the second and fifth ring roads. (2) Class 1 is dominated by transport blocks, Class 3 is dominated by commercial blocks, and Class 5 is dominated by green space blocks, with a clustering index of 0.38. (3) The high-density, high-height class (HH-Class 2) has a greater number of blocks distributed in a ring shape around the periphery of the second ring road. The high-density, low-height class (HL-Class 2) has a relatively small number of blocks but a relatively large area, and the largest blocks are located in the western part of the study area. (4) In the HH and HL building patterns, extreme heat scenarios often occur; from the perspective of functional areas, the probability of extreme heat in the transport block is much higher than that of other functional areas, and except for the HH scenario, the green space functional area plays a very important role in reducing the temperature. This study explores the characteristics of the built environment that influence the urban LST from the perspective of different urban functional zones in cities to provide decision support for quantitative territorial spatial planning, optimization, and management.

Evaluating MPAS-A Performance for Mesoscale Simulation in a Tropical Region: A Case Study of Extreme Heat in Jakarta, Indonesia

The Model for Prediction Across Scales–Atmosphere (MPAS-A) has been widely used for larger scale simulations, but its performance in mesoscale, particularly in tropical regions, is less evaluated. This study aimed to assess MPAS-A in simulating extreme surface air temperature in Jakarta during the hot spells of October 2023 with eight different simulation setups. Several validation metrics were applied to near-surface meteorological variables, land surface temperature (LST), and vertical atmospheric profile. From the eight simulations, MPAS-A captured diurnal patterns of the near-surface variables well, except for wind direction. The model also performed well in LST simulations. Moreover, the biases in the vertical profiles varied with height and were sensitive to the initial/boundary conditions used. Simulations with modified terrestrial datasets showed higher LST and air temperatures over the sprawling urban areas. MPAS-A successfully simulated the extreme event, showing higher air temperatures in southern Jakarta (over 36 °C) compared to the northern part. Negative temperature advection by sea breeze helped lower air temperature in the northern area. This study highlights the role of sea breezes as natural cooling mechanisms in coastal cities. Additionally, MPAS-A is feasible for several applications for urban climate studies and climate projection, although further development is needed.

Analysis of radiative heat flux using ASTER thermal images: Climatological and volcanological factors on Java Island, Indonesia

This study focuses on analysing natural Radiative Heat Flux (RHF) anomalies to map out the heat distribution across the Java Island. Leveraging remote sensing techniques, we calculated natural RHF anomalies using Land Surface Temperature (LST) and Land Surface Emissivity (LSE) data obtained from Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) imagery. A key aspect of our approach was distinguishing between natural and anthropogenic heat sources by cross-referencing the LST Map with the Land Use Land Cover (LULC) map of Java Island. The study interprets natural RHF anomalies by examining regional trends in non-volcanic areas and local trends within volcanic regions, considering climatological and volcanological factors. Relation with climatological factors involves assessing soil moisture parameters from Soil Moisture Active Passive (SMAP) data, precipitation from monthly Global Precipitation Measurement (GPM) data, and classifications according to the Köppen-Geiger climate schema. Our regional analysis reveals high natural RHF anomalies in the northern regions of West Java, parts of Central Java, and most of East Java, attributed to low soil moisture and low precipitation in savanna and monsoon climates. On a more localised scale, RHF values are significantly high in volcanic areas, particularly around Central and East Java's volcanoes, such as Mt. Merapi, Mt. Slamet, Mt. Semeru, the Sidoarjo Mud Volcano, and Mt. Ijen. The Natural RHF anomalies at volcanoes in West Java were identified as not being high except at Mt Patuha. These areas exhibit average natural RHF anomalies ranging between 32.22 W/m2 and 115.13 W/m2, indicating strong and intense volcanic activity. The insights obtained from these findings explain the overall thermal characteristics of Java Island and highlight the presence of subsurface thermal zones associated with volcanic activity and geothermal potential.

Exploring Urban Heat Distribution via Intra- and Extra-Block Morphologies with Integrated Stacked Models

The spatial variability of land surface temperature (LST) is considerably affected by urban morphology. Previous research has focused separately on the thermal effects of urban morphology and the cooling effects of water bodies and urban parks. However, the combined influence of intra- and extra-block factors on LST has not been thoroughly examined. To bridge this research gap, we conducted an extensive analysis of 17 urban morphology factors in Hangzhou by employing a novel stacked ensemble approach. Results showed that the stacked ensemble models outperformed commonly used techniques, such as random forest and boosted regression trees. Extra-block factors, alongside building density, average building height, and vegetation coverage within blocks, predominantly influenced the LST distribution across all seasons. Building density was positively correlated with LST, with a maximum influence of 1.5 °C in spring, whereas building height was negatively correlated with it, with a maximum influence of 1.8 °C in winter. The cooling distance of the Qiantang River extends up to 2500 m into the urban blocks and has a maximum effect of 2 °C in summer. These insights deepen our comprehension of the interplay between LST and intra- and extra-block urban morphologies, thus offering valuable guidance for urban planners and policymakers.

Evaluating urban heat island to achieve sustainable development goals: A case study of Tiruchirappalli city, India

Tiruchirappalli City the fourth largest municipal corporation in Tamilnadu, India is witnessing accelerated urban growth. The city has exacerbated the microclimate indicating a significant increase in land surface temperature. The study aims to establish the differentiation of land surface temperature of city between 2001, 2011, and 2021 of summer months using Landsat images to examine the effects of urban heat island. The results shows that the temperature ranged from 31.23°C to 37.45°C between 2011 and 2021. The assessed land surface temperature (LST) is correlated and the result shows a positive correlation between the LST-normalised difference built-up index (NDBI) and a negative correlation between the LST-normalised difference vegetation index (NDVI) and a weak correlation between LST-BSI and LST-MNDWI of summer months examined. Further, the urban thermal field variance index (UTFVI) of the city is analysed. The thermal condition of the city during April has significantly changed from 5.14 % of area in 2001 to 12.19 % of area in 2021 (UTFVI >0.020, Strongest) of worst ecological valuation index. The selected city is experiencing major thermal conditions (UTFVI 0.010–0.015) of ecological evaluation index from bad to worst. This study will help mitigate the urban heat island effect and promote sustainable development for rapid urbanization.

The Characteristics of Precipitation with and without Bright Band in Summer Tibetan Plateau and Central-Eastern China

The bright band (BB) is an important symbol of the ice–water transition zone in stratiform precipitation, and the presence or absence of BB will lead to different microphysical processes. In this paper, the characteristics of BB and precipitation characteristics with and without BB in summer at Tibetan Plateau (TP) as well as Central-eastern China (CEC) are analyzed by using Global Precipitation Measurement (GPM) and the fifth generation ECMWF atmospheric reanalysis of the global climates (ERA5) datasets. The results show the freezing level height and BB height in TP are 0.5 km higher than those in CEC. With the increase in rain rate, the BB height decreases in TP but increases in CEC. The BB width becomes wider with the increase in maximum radar reflectivity. Secondly, the maximum reflectivity factor and particle diameter of stratiform precipitation with BB appear at 5 km, while the maximum reflectivity factor of stratiform precipitation without BB and convective precipitation appear near the ground. The particle diameter first decreases and then increases from the cloud top to the ground. Thirdly, the land surface temperature of convective precipitation is about 2.5 °C higher than stratiform precipitation with BB, indicating higher land surface temperatures are more likely to trigger convection. Lastly, BB can lead to a decrease in brightness temperature and an increase in polarized difference at 89 GHZ and 166 GHZ in CEC, likely due to the increasing ice particles in stratiform precipitation with BB.

Surface energy balance-based surface urban heat island decomposition at high resolution

Urban heat island (UHI) is among the most pronounced human impacts on Earth. To formulate locally adapted mitigation strategies, a comprehensive understanding of the influencing mechanisms of UHI at high resolution is imperative. Based on surface energy balance, we attributed surface UHI (SUHI) into five biophysical terms (surface radiation, anthropogenic heat, convection, evapotranspiration and heat storage term) using Sentinel-2 and Landsat-8 images in Beijing. The simulated SUHI intensity, derived by combining all five contribution terms, exhibited a good consistency but a higher spatial resolution, than SUHI intensity extracted from Landsat-8 land surface temperature product. SUHI intensity tended to decrease from the old city to outsides, attributed to the decrease of evapotranspiration, solar radiation and anthropogenic heat term. The convection and heat storage term play a positive role in reducing SUHI. Among urban morphological blocks, low-rise and high-density blocks had the strongest SUHI, with the evapotranspiration term contributing the most. The results highlighted the capacity of the urban surface to evaporate water in affecting Beijing SUHI. The proposed method provides one useful tool to analyze the drivers of SUHI from the aspect of heat formation, which can be potentially applied worldwide for large-scale comparisons of how urbanization affects UHI.