Land Surface Temperature Products Research Articles

Correcting an Off-Nadir to a Nadir Land Surface Temperature Using a Multitemporal Thermal Infrared Kernel-Driven Model during Daytime

Land surface temperature (LST) is a fundamental parameter in global climate, environmental, and geophysical studies. Remote sensing is an essential approach for obtaining large-scale and frequently updated LST data. However, due to the wide field of view of remote sensing sensors, the observed LST with diverse view geometries suffers from inconsistency caused by the thermal radiation directionality (TRD) effect, which results in LST products being incomparable, especially during daytime. To address this issue and correct current off-nadir LSTs to nadir LSTs, a semi-physical time-evolved kernel-driven model (TEKDM) is proposed, which depicts multitemporal TRD patterns during the daytime. In addition, we employ a Bayesian optimization method to calibrate seven unknown parameters in the TEKDM. Validation results using the U.S. Climate Reference Network (USCRN) sites show that the RMSE (MBE) for GOES-16 and MODIS off-nadir LST products is reduced from 3.29 K (−2.0 K) to 2.34 K (−0.02 K), with an RMSE reduction of 0.95 K (29%) and a significant reduction in systematic bias. Moreover, the proposed method successfully eliminates the angular and temporal dependence of the LST difference between the satellite off-nadir LST and in situ nadir LST. In summary, this study presents a feasible approach for estimating the high-accuracy nadir LST, which can enhance the applicability of LST products in various domains.

Reconstruction of Hourly FY-4A AGRI Land Surface Temperature under Cloud-Covered Conditions Using a Hybrid Method Combining Spatial and Temporal Information

Land Surface Temperature (LST) products obtained by thermal infrared (TIR) remote sensing contain considerable blank areas due to the frequent occurrence of cloud coverage. The studies on the all-time reconstruction of the cloud-covered LST of geostationary meteorological satellite LST products are relatively few. To accurately fill the blank area, a hybrid method for reconstructing hourly FY-4A AGRI LST under cloud-covered conditions was proposed using a random forest (RF) regression algorithm and Savitzky-Golay (S-G) filtering. The ERA5-Land surface cumulative net radiation flux (SNR) reanalysis data was first introduced to represent the change in surface energy arising from cloud coverage. The RF regression method was used to estimate the LST correlation model based on clear-sky LST and the corresponding predictor variables, including the normalized difference vegetation index (NDVI), the normalized difference water index (NDWI), surface elevation and slope. The fitted model was then applied to reconstruct the cloud-covered LST. The S–G filtering method was used to smooth the outliers of reconstructed LST in the temporal dimension. The accuracy evaluation was performed using the measured LST of the representative meteorological stations after scale correction. The coefficients of determination derived with the reference LST were all above 0.73 on the three examined days, with a bias of −1.13–0.39 K, mean absolute errors (MAE) of 1.46–2.4 K, and root mean square errors (RMSE) of 1.77–3.2 K. These results indicate that the proposed method has strong potential for accurately restoring the spatial and temporal continuity of LST and can provide a solution for the production and research of gap-free LST products with high temporal resolution.

Solar zenith angle-based calibration of Himawari-8 land surface temperature for correcting diurnal retrieval error characteristics

The geostationary Himawari-8 satellite offers a unique opportunity to monitor sub-daily thermal dynamics over Asia and Oceania, and several operational land surface temperature (LST) retrieval algorithms have been developed for this purpose. However, studies have reported inconsistency between LST data obtained from geostationary and polar-orbiting platforms, particularly for daytime LST, which can arise from variations in viewing geometries and inherent differences in sensor types and LST algorithms. Despite this, previous research has primarily focused on analysing the directionality of LST without thoroughly exploring systematic differences between platforms. Hence, we presented a Solar Zenith Angle-based Calibration (SZAC) method to harmonise the daytime component of a split-window retrieved Himawari-8 LST (referred to here as the baseline) with the MODerate-resolution Imaging Spectroradiometer (MODIS) LST. SZAC describes the spatial heterogeneity and magnitude of diurnal LST discrepancies from different platforms, which is anticipated to complement typical directionality analyses. We evaluated the harmonised LST data, referred to as the Australian National University LST with SZAC (ANUSZAC), against MODIS LST and the Visible Infrared Imaging Radiometer Suite (VIIRS) LST, as well as in-situ LST from the OzFlux network. Two peer Himawari-8 LST products from Chiba University and the Copernicus Global Land Service were also collected for comparison. The median daytime bias of ANUSZAC LST against Terra-MODIS LST, Aqua-MODIS LST and VIIRS LST was 1.52 K, 0.98 K and −0.63 K, respectively, which demonstrated improved performance compared to baseline (5.37 K, 4.85 K and 3.02 K, respectively) and Chiba LST (3.71 K, 2.90 K and 1.07 K, respectively). All four Himawari-8 LST products showed comparable performance of unbiased root mean squared error (ubRMSE), ranging from 2.47 to 3.07 K, compared to LST from polar-orbiting platforms. In the evaluation against in-situ LST, the mean values of bias (ubRMSE) of baseline, Chiba, Copernicus and ANUSZAC LST during daytime were 4.23 K (3.74 K), 2.16 K (3.62 K), 1.73 K (3.31 K) and 1.41 K (3.24 K), respectively, based on 171,289 hourly samples from 20 OzFlux sites across Australia between 01/Jan/2016 and 31/Dec/2020. In summary, the SZAC method offers a promising approach to enhance the reliability of geostationary LST retrievals by incorporating the spatiotemporal characteristics observed by accurate polar-orbiting LST data. SZAC can also reduce LST angular effects due to its capability of quantifying spatial inhomogeneity of surface heat dynamics. Furthermore, it is possible to implement SZAC using LST data acquired by geostationary satellites in other regions, e.g., Europe, Africa and Americas. This could improve our understanding of the error characteristics of coincident imageries from geostationary and polar-orbiting platforms, allowing for targeted refinements and global harmonisations to further enhance applicability.

The contrasting trend of global urbanization-induced impacts on day and night land surface temperature from a time-series perspective

Increasing impervious surface area (ISA) is a major contributor to the rising land surface temperature (LST). Understanding the warming effect of ISA quantitively is crucial to alleviate the urban heat island effect. However, the quantitative analyses of urbanization-induced LST change globally remain to be determined over past decades. Here, we first characterized the response of LST to ISA (i.e., δLST) in 369 global cities using Terra MODIS LST product (1km) and the global artificial impervious area dataset (30m). We then investigated the spatiotemporal changes of δLST and its associated driving factors, including enhanced vegetation index (1km), albedo (500m), population (100m), air temperature (0.1°), and precipitation (0.1°). The analysis was conducted using linear regression models from 2000 to 2010 at the pixel level. The results illustrate that (1) the δLST in urban areas shows a negative trend during daytime (–0.2971°C/% and –0.0870°C/%) and a negative trend at night (0.0938°C/% and 0.0048°C/%) in arid and temperate regions, respectively. While in tropical regions, the δLST remains positive throughout the day, but it exhibits an opposite trend in snow regions. (2) The magnitude of nighttime δLST is stronger compared to daytime globally, which is more significant in arid regions with the largest ΔδLST values (i.e., diurnal variations) in urban (0.3909°C/%) and rural (0.3262°C/%) areas. Besides, the impervious surface has a distinct impact on LST in summer and winter, particularly at night. (3) In temperate regions, δALB (negative) and δAT (positive) are dominant factors in regulating δLST during daytime and nighttime. However, in arid regions, the daytime and nighttime δLSTs are negatively correlated by urban population and vegetation, respectively. These findings provide a quantitative understanding on long-term LST variations and driving factors behind at the pixel level, suggesting that urban greening and increasing surface albedo are effective strategies to mitigate the urbanization-induced surface warming.

A satellite view of the exceptionally warm summer of 2022 over Europe

Abstract. Summer heatwaves are becoming increasingly dangerous over Europe, and their close monitoring is essential for human activities. Typically, they are monitored using the 2 m temperature from meteorological weather stations or reanalysis datasets. In this study, the 2022 extremely warm summer over Europe is analysed using satellite land surface temperature (LST), specifically the LSA SAF (Land Surface Analysis Satellite Application Facility) all-sky LST product (available from 2004 onwards). Since climate applications of LST are still poorly explored, heatwave diagnostics derived from satellite observations are compared with those derived using ERA5/ERA5-Land reanalysis data. Results highlight the exceptionality of 2022 in different metrics such as the mean LST anomaly, area under extreme heat conditions, number of hot days and heatwave magnitude index. In all metrics, 2022 ranked first when compared with the remaining years. Compared to 2018 (next in all rankings), 2022 exceeded its LST anomaly by 0.7 °C and each pixel had on average 7 more hot days. Satellite LST complements reanalysis diagnostics, as higher LST anomalies occur over areas under severe drought, indicating a higher control and amplification of the heatwave by surface processes and vegetation stress. These cross-cutting diagnostics increase the confidence across satellite data records and reanalyses, fostering their usage in climate applications.

Validation of Remotely Sensed Land Surface Temperature at Lake Baikal’s Surroundings Using In Situ Observations

The accuracy of Land Surface Temperature (LST) products retrieved from satellite data in mountainous-coastal areas is not well understood. This study presents an analysis of the spatial and temporal variability of the differences between the LST and in situ observed air and surface temperatures (ISTs) for the southeastern slope of Lake Baikal’s surroundings. The IST was measured at 12 ground observation sites located on the southeastern macro-slope of the Primorskiy Ridge (Baikal, Russia) within an elevation range of 460–1656 m above sea level from 2009 to 2021. LST was estimated using 617 Landsat (7 and 8) images from 2009–2021, taking into account brightness temperature, surface emissivity and vegetation cover fraction. The comparison of the LST from satellite data and the IST from ground observation showed relatively high differences, which varied depending on the season and site type. A neural network was suggested and calibrated to improve the LST data. The corrected remote-sensed temperature was found to reproduce the IST very well, with mean differences of about 0.03 °C and linear correlation coefficients of 0.98 and 0.95 for the air and surface IST.

Filling gaps in cloudy Landsat LST product by spatial-temporal fusion of multi-scale data

Land surface temperature (LST) is an important factor in studies of surface energy fluxes between the Earth's surface and atmosphere. The Landsat LST product has been applied widely due to its fine spatial resolution and high data quality. Frequent cloud coverage, however, results in different degrees of gaps in the Landsat LST images, restricting greatly their application in practical cases requiring continuous LST time-series. Gap filling is an important technique used to reconstruct the missing data in the Landsat LST images. Using only temporally neighboring data from the same sensor, however, current gap filling methods struggle to deal with large spatial gaps and large temporal changes in LST. To solve these challenging issues, in this paper, multi-scale data (i.e., the 2 km, all-sky, hourly GOES LST product) were proposed for gap filling of Landsat LST images, and a spatial-temporal random forest (STRF) method was developed. The GOES LST product can provide effective information on LST in the gaps of the target Landsat LST data, with great potential to decreasing the uncertainty in gap filling when large LST change exists or the gap size is large. Through experiments in six regions of the USA, the proposed STRF method was shown to be more accurate than three typical gap filling and spatio-temporal fusion methods, and particularly advantageous for cases involving large temporal changes or large spatial gaps in LST. Specifically, STRF produces an average root mean square error (RMSE) of 1.6444 K and average correlation coefficient (CC) of 0.8653 for the six regions, which are 0.1133 K smaller and 0.0176 larger than the most competitive benchmark method (i.e., the spatial random forest (SRF) method). Conceptually, by spatial-temporal fusion of multi-scale data in gap filling, STRF provides a new paradigm for both spatial-temporal fusion and gap filling. Furthermore, the proposed STRF has the potential to be applied to various gap filling tasks in other situations.

The impact of clear-sky biases of land surface temperature on monthly evapotranspiration estimation

Remotely sensed land surface temperature (LST) is critical for retrieving evapotranspiration (ET). However, due to cloud contamination, LST is often limited to clear-sky conditions and the differences between clear-sky and all-sky LST will lead to clear-sky biases of LST. Consequently, the accuracy of ET varies drastically under different weather conditions. To evaluate the impact of clear-sky biases on ET estimation, the monthly ET under clear-sky and all-sky conditions was estimated using a nonparametric method at several sites with different humid conditions and vegetation coverage in 2014. 35.7 % of the clear-sky LST values were acquired in the arid region (Heihe River Basin) with different vegetation coverage. As the climate became more humid, only 14.2 % of the clear-sky LST data were available (Poyang Lake Basin). The clear-sky biases of LST variation at the sites resulted in a significant reduction in the accuracy of monthly ET, with an increase in the relative error (RE) of approximately 16.6 %. The impact of clouds can be reduced by at most half by the introduction of all-sky LST products, with a significant decline (6.3 % ∼ 34.2 % of ET) in the error contribution of LST to monthly ET. The variation in vegetation coverage of land cover and the change in humidity of the climate significantly affected the influence of the clear-sky biases of LST on the ET estimation and the error contribution to ET. Replacing all-sky LST with clear-sky LST in areas with less vegetation coverage exacerbates the underestimation of ET estimation and strengthens the error contribution of LST with the decline in RE (error contribution of LST) from −7.7 % (−6.3 %) in densely vegetated areas to −29.7 % (−34.2 %) in non-vegetated areas, especially in the summer. Similarly, clear-sky biases of LST tended to have a relatively significant impact on the error contributions of LST to monthly ET estimation when the climate became humid, with a relatively significant enlarging from −6.3 % in the arid area to −10.0 % in the humid area, especially in the rainy season in the humid area (−61.3 % of ET). This study contributes to the understanding of the relationship between clear-sky biases and vegetation coverage/humidity. Quantitative analysis of the impact of clear-sky biases of LST on monthly ET estimation will be helpful for improving the accuracy of ET retrieval in the future.

Estimating carbon emissions from thermal power plants based on thermal characteristics

Carbon emissions from thermal power plants (TPPs) are a significant source of anthropogenic carbon and a priority for carbon reduction efforts. The accurate estimating of carbon emissions from TPPs is essential research in the field of carbon reduction. In this study, we take TPPs located in the United States as examples to discuss the relationship between the thermal characteristics of the power plant and its annual carbon emissions. Our data consists of land surface temperature (LST) product generated from Landsat 8 thermal data, LST product from ERA5-Land reanalysis dataset, and carbon emission data reported by the U.S. Energy Information Administration (EIA). Based on pertinent research of urban heat islands, we propose the concept of thermal characteristic index (TCI), which can provide a quantitative description of power plant thermal characteristics. Furthermore, we convert TCI to annual thermal characteristic index (A-TCI) by annual end-point interpolation and Simpson integration of the TCI series to match the annual carbon emission data and then establish a TCI-based annual carbon emission estimation model (TACEE-Model). Using the TACEE-Model, we predict the annual carbon emission values and evaluate the model’s performance. Experimental results indicate that the TACEE-Model exhibits great predictive capability, with the mean absolute percentage error (MAPE) of 0.0325, and the average prediction percentage error of 3.24%. The application of annual end-point interpolation and Simpson integration can significantly enhance the model’s prediction accuracy, resulting in a 51.42% reduction in MAPE. Overall, this research presents a new method for verifying annual carbon emissions from TPPs, which can contribute to the supervisory oversight of carbon emission reduction from TPPs.

Multi-scale remote sensing retrieval model for land surface emissivity and its impact on land surface temperature estimation

ABSTRACT As a fundamental parameter of the surface radiation process, land surface emissivity (LSE) is the most direct factor affecting the land surface temperature (LST) retrieval accuracy. Based on the measured data from 2017 to 2023, the characteristics of LSE changes were analysed, different timescale LSE retrieval models were developed, and the impact of LSE on pixel scale observational LST acquisition was evaluated by using unmanned aerial vehicle (UAV) remote sensing images. The results show that with changes in atmospheric environmental conditions, the LSE changes are relatively less pronounced at 3.99 µm and 4.02 µm in the mid-infrared band and at 10.8 µm and 11.8 µm in the thermal infrared band than in the other bands. At the timescale above day, the LSE is significant correlation with the normalized differential vegetation index (NDVI). Vegetation is the main parameter that dominates the LSE change under heterogeneous underlying surfaces. On a daily timescale, LSE exhibits typical diurnal variation characteristics, and the amplitude of diurnal variation is influenced mainly by shallow soil water and heat conditions. The developed LSE model based on Fengyun-4A (FY-4A) satellite remote-sensing data can effectively reproduce the daily variation characteristics of LST and avoid the ‘jumping’ phenomenon that may occur in the NDVI threshold method. The root-mean-square error (RMSE) between the retrieved LST and in situ data is reduced to 3.78°C, and the retrieval accuracy is better than that of the published LST product. When using UAVs to obtain pixel-scale LSTs, dry bare soil is more sensitive to LSE. When the LSE increases from 0.95 to 1.0, the average LST obtained by the UAV deviates by approximately 1.0°C. This study preliminarily provides some spatiotemporal variation patterns of LSE and lays a foundation for the later-modified LST retrieval algorithms.

Geothermal target detection integrating multi-source and multi-temporal thermal infrared data

Thermal infrared remote sensing (TIRS) technology is intriguing for geothermal anomaly detection due to its time efficiency and cost-effectiveness. Land surface temperature (LST) derived from TIRS indicates potential geothermal anomalies. However, LST on different natural features, exhibiting relatively cold/hot anomalies in daytime and nighttime, affected the identification of abnormal areas. The aim of this study was to highlight geothermal anomalous areas by integrating multi-view daytime and nighttime LST data. Specifically, the winter LST time series of drillings on 2013–2023 was processed based on Landsat-8 TIRS in diurnal scenarios. An information aggregation classification method about Dempster-Shafer evidence theory was presented to target recognition at nighttime by the ASTER LST product. This day-night information fusion analysis effectively highlighted potential areas with geothermal anomalies. The results demonstrated that the geothermal anomalies were mainly distributed along the Ruili-Longling Fault and Wanding Fault in a Northeast (NE) – Southwest (SW) direction, and the NE side was more significant than the SW side. These findings suggested that the east side was closer to the magma heat source, and the presumed magma heat source location aligned closely with the detected geothermal anomalies. Geological data were employed for geological interpretation of the LST anomaly area, eliminating non-geothermal influences. Along with drilling data verification, this method was successfully used to identify geothermal anomaly areas in Ruili City, Yunnan Province. Overall, this study provided valuable insights into the detection of geothermal anomalies through TIRS, contributing to the successful development of geothermal resources.

A Spatial Downscaling Method for Solar-Induced Chlorophyll Fluorescence Product Using Random Forest Regression and Drought Monitoring in Henan Province

Drought is a frequent global phenomenon. Solar-induced chlorophyll fluorescence (SIF), an electromagnetic signal, has been proven to be an efficient tool for monitoring and assessing gross primary productivity (GPP) and drought. To address the issue of the sparse resolution of satellite-based SIF, researchers have developed different downscaling algorithms. Recently, the most frequently used SIF products had a spatial resolution of 0.05 degrees. However, these spatial resolution SIF data are not conducive to regional agricultural drought monitoring. In this study, we utilized the global ‘OCO-2’ solar-induced fluorescence (GOSIF) products along with normalized difference vegetation index (NDVI) and land surface temperature (LST) products. With the powerful advantages offered by Google Earth Engine (GEE), we could conveniently acquire the necessary data. Additionally, employing the random forest (RF) method, we successfully acquired downscaled SIF data at an enhanced spatial resolution of 1 km. Using those downscaled SIF results with 1 km resolution, an SIF anomaly index was established and calculated to monitor drought. Results showed that the RF-based downscaled SIF result followed the same trend as the GOSIF value. Subsequently, correlation coefficients between SIF and GPP were calculated. The downscaled SIF demonstrated a higher correlation with GPP from MODIS compared to 0.05-degree GOSIF, with coefficients of 0.74 and 0.68 in May 2018, respectively. Moreover, the SIF anomaly index showed positive correlations with crop yield; the correlation coefficients were 0.93 for wheat and 0.89 for maize. The drought index had a negative correlation with areas affected by drought, with a correlation coefficient of −0.58. Finally, the SIF anomaly index was used to monitor drought from 2001 to 2020 in Henan Province. The 1 km SIF results obtained through the RF-based downscaled method were deemed reliable, thereby establishing the suitability of the SIF anomaly index for drought monitoring at a regional scale.

Quantifying of surface urban heat island intensity in Isfahan metropolis using MODIS\Terra\LST data

Heat island characteristics depend on the background climate of the site where the city is located. Therefore, an index was defined for the Isfahan metropolitan area to quantify the surface urban heat island intensity. This new index is based on the representative pixels of urban and non-urban areas. For this purpose, MODIS land cover type product (MCD12Q1) data were used to distinguish between urban and non-urban areas. Also, data from the MODIS/Terra land surface temperature product (MOD11A1) from 2000 to 2018 were utilized for daytime and nighttime to study the surface heat island intensity. Then, the representative pixels of urban and non-urban areas were identified using the spatial correlation method, and the heat island index was calculated for the metropolitan area of Isfahan. The study showed that the frequency distribution of the nighttime heat island index follows a normal distribution and is often 3.5 to 4º K above the temperature of the surrounding areas of the city. The 365-day floating mean of the surface urban heat island reveals that this index has increased in recent years. The research of temporal behavior showed that the intensity of the surface urban heat island reaches its maximum in January and becomes weaker in summer, while the survey of spatial behavior showed that the core of the surface urban heat island extends towards downtown areas, where the oldest part of the city is located.

Spatiotemporal variation of land surface temperature and its driving factors in Xinjiang, China

Land surface temperature (LST) directly affects the energy balance of terrestrial surface systems and impacts regional resources, ecosystem evolution, and ecosystem structures. Xinjiang Uygur Autonomous Region is located at the arid Northwest China and is extremely sensitive to climate change. There is an urgent need to understand the distribution patterns of LST in this area and quantitatively measure the nature and intensity of the impacts of the major driving factors from a spatial perspective, as well as elucidate the formation mechanisms. In this study, we used the MOD11C3 LST product developed on the basis of Moderate Resolution Imaging Spectroradiometer (MODIS) to conduct regression analysis and determine the spatiotemporal variation and differentiation pattern of LST in Xinjiang from 2000 to 2020. We analyzed the driving mechanisms of spatial heterogeneity of LST in Xinjiang and the six geomorphic zones (the Altay Mountains, Junggar Basin, Tianshan Mountains, Tarim Basin, Turpan-Hami (Tuha) Basin, and Pakakuna Mountain Group) using geographical detector (Geodetector) and geographically weighted regression (GWR) models. The warming rate of LST in Xinjiang during the study period was 0.24°C/10a, and the spatial distribution pattern of LST had obvious topographic imprints, with 87.20% of the warming zone located in the Gobi desert and areas with frequent human activities, and the cooling zone mainly located in the mountainous areas. The seasonal LST in Xinjiang was at a cooling rate of 0.09°C/10a in autumn, and showed a warming trend in other seasons. Digital elevation model (DEM), latitude, wind speed, precipitation, normalized difference vegetation index (NDVI), and sunshine duration in the single-factor and interactive detections were the key factors driving the LST changes. The direction and intensity of each major driving factor on the spatial variations of LST in the study area were heterogeneous. The negative feedback effect of DEM on the spatial differentiation of LST was the strongest. Lower latitudes, lower vegetation coverage, lower levels of precipitation, and longer sunshine duration increased LST. Unused land was the main heat source landscape, water body was the most important heat sink landscape, grassland and forest land were the land use and land cover (LULC) types with the most prominent heat sink effect, and there were significant differences in different geomorphic zones due to the influences of their vegetation types, climatic conditions, soil types, and human activities. The findings will help to facilitate sustainable climate change management, analyze local climate and environmental patterns, and improve land management strategies in Xinjiang and other arid areas.

The anisotropy of MODIS LST in urban areas: A perspective from different time scales using model simulations

Remote sensing is one of the effective means to obtain urban land surface temperature (LST), but the observed temperature varies with sensor viewing angle due to urban thermal anisotropy (UTA) and biased sensor viewing angle. The anisotropy of satellite-based LST products (e.g., MODIS LST) varies at different time scales. Previous researches focus on the anisotropy of MODIS LST at the seasonal scale or at other time scales (e.g., hourly) but only for individual cities. The characteristics of anisotropy at different time scales and for cities at different latitudes are still unknown. In this study, we focus on analyzing the characteristics of anisotropy at hourly, daily, seasonal, and annual time scales separately for 80 global cities using the simulation data. The simulation data were generated by coupling a microscale urban surface temperature model and a thermal remote sensing model. In addition, the impact factors (i.e., sensor factor, meteorological factor, and surface factor) at different time scales are discussed and validated using satellite observations. The results show that the daytime instantaneous anisotropy varies from −4.2 K to 1.1 K at the hourly scale. For each city, the smaller anisotropy is captured at the local time 12:00 more easily for both Terra and Aqua overpass during daytime. This is because the magnitude of urban thermal anisotropy is more influenced by the sensor factor (i.e., sensor zenith angle), which varies with the overpass time and the relative position of the city in the track. However, if the seasonally or yearly aggregated MODIS LST is utilized to study mean surface temperature, the intensity of angular effect is much smaller with the value ranging from 0.1 K to 3.3 K at the seasonal scale (from 0.2 K to 2.1 K at the yearly scale) because the variation of sensor zenith angle is smoothed. At seasonal and yearly scales, the strength of the angular effect is greater in the morning during daytime, in summer among all seasons, and in the cities with lower latitudes, which is due to larger solar radiation. Among all meteorological factors, wind speed also has a larger impact on the magnitude of surface thermal anisotropy besides the downward surface shortwave radiation. Our study provides a reference for the correction of the angular effect in the MODIS LST products.

Downscaling Land Surface Temperature Derived from Microwave Observations with the Super-Resolution Reconstruction Method: A Case Study in the CONUS

Optical sensors cannot penetrate clouds and can cause serious missing data problems in optical-based Land Surface Temperature (LST) products. Under cloudy conditions, microwave observations are usually utilized to derive the land surface temperature. However, microwave sensors usually have coarse spatial resolutions. High-Resolution (HR) LST data products are usually desired for many applications. Instead of developing and launching new high-resolution satellite sensors for LST observations, a more economical and practical way is to develop proper methodologies to derive high-resolution LSTs from available Low-Resolution (LR) datasets. This study explores different algorithms to downscale low-resolution LST data to a high resolution. The existing regression-based downscaling methods usually require simultaneous observations and ancillary data. The Super-Resolution Reconstruction (SRR) method developed for traditional image enhancement can be applicable to high-resolution LST generation. For the first time, we adapted the SRR method for LST data. We specifically built a unique database of LSTs for the example-based SRR method. After deriving the LST data from the coarse-resolution passive microwave observations, the AMSR-E at 25 km and/or AMSR-2 at 10 km, we developed an algorithm to downscale them to a 1 km spatial resolution with the SRR method. The SRR downscaling algorithm can be implemented to obtain high-resolution LSTs without auxiliary data or any concurrent observations. The high-resolution LSTs are validated and evaluated with the ground measurements from the Surface Radiation (SURFRAD) Budget Network. The results demonstrate that the downscaled microwave LSTs have a high correlation coefficient of over 0.92, a small bias of less than 0.5 K, but a large Root Mean Square Error (RMSE) of about 4 K, which is similar to the original microwave LST, so the errors in the downscaled LST could have been inherited from the original microwave LSTs. The validation results also indicate that the example-based method shows a better performance than the self-similarity-based algorithm.

Diurnal variation in the urban thermal environment and its relationship to human activities in China: a Tencent location-based service geographic big data perspective.

The main factor of the formation and deterioration in China's urban thermal environment is human activity, which is difficult to describe and measure. A new perspective on the effect of human activity on the urban thermal environment can be obtained by examining the interaction between location-based service (LBS) data and the urban thermal environment in China. However, relevant research is still limited. In this study, we used Tencent LBS data, Terra/Aqua MODIS land surface temperature (LST) data, and land use data to investigate the relationship between LBS and the urban thermal environment, specifically the LST and surface urban heat island intensity (SUHII) across China and its provinces. Our results showed that (1) in summer, the heat island effect was an issue in 94% of the urban areas in China, which was worse during the day. The high- and low-value periods of LBS data on a given day coincided with the acquisition times of MODIS LST products during the day and at night, respectively. (2) During both the day and at night, there was a significant connection between LBS data and the urban thermal environment in China. The highest correlation coefficient (r) between LBS data and the LST could reach 0.55 (p < 0.01) at the provincial level, and the highest correlation coefficient (r) between LBS data and the SUHII could reach 0.78 (p < 0.01) at the provincial level. (3) The urban thermal environment diurnal difference and LBS data exhibited a significant relationship. The ΔLBS diurnal differences were significantly positively related to the SUHII diurnal differences in China. The overall study findings revealed that LBS data constitute an important parameter to represent the human activity intensity when investigating the formation of the urban thermal environment in China.

High-resolution (1 km) all-sky net radiation over Europe enabled by the merging of land surface temperature retrievals from geostationary and polar-orbiting satellites

Abstract. Surface net radiation (SNR) is a vital input for many land surface and hydrological models. However, most of the current remote sensing datasets of SNR come mostly at coarse resolutions or have large gaps due to cloud cover that hinder their use as input in models. Here, we present a downscaled and continuous daily SNR product across Europe for 2018–2019. Long-wave outgoing radiation is computed from a merged land surface temperature (LST) product in combination with Meteosat Second Generation emissivity data. The merged LST product is based on all-sky LST retrievals from the Spinning Enhanced Visible and InfraRed Imager (SEVIRI) onboard the geostationary Meteosat Second Generation (MSG) satellite and clear-sky LST retrievals from the Sea and Land Surface Temperature Radiometer (SLSTR) onboard the polar-orbiting Sentinel-3A satellite. This approach makes use of the medium spatial (approx. 5–7 km) but high temporal (30 min) resolution, gap-free data from MSG along with the low temporal (2–3 d) but high spatial (1 km) resolution of the Sentinel-3 LST retrievals. The resulting 1 km and daily LST dataset is based on an hourly merging of both datasets through bias correction and Kalman filter assimilation. Short-wave outgoing radiation is computed from the incoming short-wave radiation from MSG and the downscaled albedo using 1 km PROBA-V data. MSG incoming short-wave and long-wave radiation and the outgoing radiation components at 1 km spatial resolution are used together to compute the final daily SNR dataset in a consistent manner. Validation results indicate an improvement of the mean squared error by ca. 7 % with an increase in spatial detail compared to the original MSG product. The resulting pan-European SNR dataset, as well as the merged LST product, can be used for hydrological modelling and as input to models dedicated to estimating evaporation and surface turbulent heat fluxes and will be regularly updated in the future. The datasets can be downloaded from https://doi.org/10.5281/zenodo.8332222 (Rains, 2023a) and https://doi.org/10.5281/zenodo.8332128 (Rains, 2023b).

Estimation of Land Surface Temperature from Chinese ZY1-02E IRS Data

The role of land surface temperature (LST) is of the utmost importance in multiple academic disciplines, such as climatology, hydrology, ecology, and meteorology. To date, many methods have been proposed to estimate LST from satellite thermal infrared data. The single-channel (SC) algorithm can provide an accurate result in retrieving LST based on prior knowledge of known land surface emissivity (LSE). The SC algorithm is extensively employed for retrieving LST from Landsat series data due to its simplicity and its reliance on just one thermal infrared channel. The Thermal Infrared Sensor (IRS) on the Chinese ZY1-02E satellite is a pivotal instrument employed for gathering thermal infrared (TIR) data of land surfaces. The objective of this research is to evaluate the feasibility of a single-channel approach based on water vapor scaling (WVS) for deriving LST from ZY1-02E IRS data because of its wide spectrum range, i.e., 7~12 μm, which is affected strongly by both atmospheric water vapor and ozone. Three study areas, namely the Baotou, Heihe River Basin, and Yantai Sea sites, were selected as validation sites to evaluate the LST inversion accuracy. This evaluation was also conducted via cross-comparison between the retrieved LST and MODIS LST products. The results revealed that the WVS-based method exhibited an average bias of 0.63 K and an RMSE of 1.62 K compared to the in situ LSTs. The WVS-based method demonstrated reasonable accuracy through cross-validation with the MODIS LST product, with an average bias of 0.77 K and an RMSE of 2.0 K. These findings indicate that the WVS-based method is effective in estimating LST from ZY1-02E IRS data.

TRIMS LST: a daily 1 km all-weather land surface temperature dataset for China's landmass and surrounding areas (2000–2022)

Abstract. Land surface temperature (LST) is a key variable within Earth's climate system and a necessary input parameter required by numerous land–atmosphere models. It can be directly retrieved from satellite thermal infrared (TIR) observations, which contain many invalid pixels mainly caused by cloud contamination. To investigate the spatial and temporal variations in LST in China, long-term, high-quality, and spatiotemporally continuous LST datasets (i.e., all-weather LST) are urgently needed. Fusing satellite TIR LST and reanalysis datasets is a viable route to obtain long time-series all-weather LSTs. Among satellite TIR LSTs, the MODIS LST is the most commonly used, and a few corresponding all-weather LST products have been reported recently. However, the publicly reported all-weather LSTs were not available during the temporal gaps of MODIS between 2000 and 2002. In this study, we generated a daily (four observations per day) 1 km all-weather LST dataset for China's landmass and surrounding areas, the Thermal and Reanalysis Integrating Moderate-resolution Spatial-seamless (TRIMS) LST, which begins on the first day of the new millennium (1 January 2000). We used the enhanced reanalysis and thermal infrared remote sensing merging (E-RTM) method to generate the TRIMS LST dataset with the temporal gaps being filled, which had not been achieved by the original RTM method. Specifically, we developed two novel approaches, i.e., the random-forest-based spatiotemporal merging (RFSTM) approach and the time-sequential LST-based reconstruction (TSETR) approach, respectively, to produce Terra/MODIS-based and Aqua/MODIS-based TRIMS LSTs during the temporal gaps. We also conducted a thorough evaluation of the TRIMS LST. A comparison with the Global Land Data Assimilation System (GLDAS) and ERA5-Land LST demonstrates that the TRIMS LST has similar spatial patterns but a higher image quality, more spatial details, and no evident spatial discontinuities. The results outside the temporal gap show consistent comparisons of the TRIMS LST with the MODIS LST and the Advanced Along-Track Scanning Radiometer (AATSR) LST, with a mean bias deviation (MBD) of 0.09/0.37 K and a standard deviation of bias (SD) of 1.45/1.55 K. Validation based on the in situ LST at 19 ground sites indicates that the TRIMS LST has a mean bias error (MBE) ranging from −2.26 to 1.73 K and a root mean square error (RMSE) ranging from 0.80 to 3.68 K. There is no significant difference between the clear-sky and cloudy conditions. For the temporal gap, it is observed that RFSTM and TSETR perform similarly to the original RTM method. Additionally, the differences between Aqua and Terra remain stable throughout the temporal gap. The TRIMS LST has already been used by scientific communities in various applications such as soil moisture downscaling, evapotranspiration estimation, and urban heat island modeling. The TRIMS LST is freely and conveniently available at https://doi.org/10.11888/Meteoro.tpdc.271252 (Zhou et al., 2021).