Land Surface Temperature Dataset Research Articles

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).

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).

A Fast and Easy Way to Produce a 1-Km All-Weather Land Surface Temperature Dataset for China Utilizing More Ground-Based Data

A Fast and Easy Way to Produce a 1-Km All-Weather Land Surface Temperature Dataset for China Utilizing More Ground-Based Data

Spatial characterization of global heat waves using satellite-based land surface temperature

The warming world greatly suffers the increase of frequency, severity and duration of heat wave events, which would cause the significant societal and environmental damages and implications from local to global scales. In order to eliminate the limitations of site observations in accurate spatial extent identification and large scale monitoring, this study tried to investigate global heat waves in 2018 using spatial datasets of land surface temperature (LST) derived from AQUA TIR sensors. A 15-year daily maximum LST dataset was used to identify the hot days and heat wave events using the relative threshold method based on the probability density function analysis of LST sequences, then they were adopted in analyzing spatial patterns of duration, frequency and intensity of heat wave events using heat wave numbers (HWN), the frequency of the appeared heat wave days (HWF) and the amplitude of heat wave impacts (HWA) indices, respectively. Our results indicated that more heat wave events and longer heat wave duration happened in low and middle latitudes under the warmer land surface in 2018 than last decade, especially many heat wave events occurred in the environmentally fragile regions and the dense population area. Urban regions achieved the largest HWA by 45.5 ± 6.9℃. Both annual HWF and HWN highlighted seven hot spot regions influenced by heat waves, and most of them distributed around the latitude 30°N and the latitude 30°S. Extreme heat wave of hot spot regions in the Southern Hemisphere mainly happened in winter and spring in contrast to that in North Hemisphere appeared in summer and autumn. Dry climate could contribute to the occurrence of heat waves, and heat wave indices variations also implied the compound consequences between heat waves and drought. Our findings demonstrate that remote sensing datasets are capable of providing the continuous whole-Earth coverage of heat wave changes, and they would play more important roles in preventing or mitigating the impacts of extreme heat events on people and natural environment.

Assessment of the solar energy–agriculture–water nexus in the expanding solar energy industry of India: An initiative for sustainable resource management

The rapid expansion of solar energy in India poses challenges in the context of cropland encroachment and water-resource scarcity. To help address these challenges, we used open access satellite observations and GIS technologies over Earth Engine platform to devise the Solar Panel Index (SPI) for efficient detection of solar farms, achieving a Land use/Land cover classification accuracy of 89 %. We have used cropland, water availability, power grid and land surface temperature data-sets for generation of land suitability map. Sentinel-2 data-sets along with SPI were used to develop solar farm locations and their estimated installed capacity within the study area. The outcomes from the study depicted that over 40 % of all solar farms in the country are located on agricultural land, with the highest seen in Karnataka (73.55 %) followed by Tamil Nadu (68.81 %). Furthermore, high installed capacity coincides with low groundwater depths, exacerbating local water strain. Madhya Pradesh was found to have almost no power plant located within 5 km of surface water bodies. Our findings emphasize the significance of the solar-agriculture-water (SAW) nexus for sustainable development. Using high-resolution Land Suitability Map (LSM) as a decision-making criteria for allocating land for power generation reveals that Rajasthan holds the highest potential for solar energy installation. Our integrated approach considering the interplay between solar-energy production, agriculture, and water-resource management contributes to the sustainable growth of the solar-energy sector and regional development.

Drought Monitoring from Fengyun Satellite Series: A Comparative Analysis with Meteorological-Drought Composite Index (MCI)

Drought is a complex natural hazard that affects various regions of the world, causing significant economic and environmental losses. Accurate and timely monitoring and forecasting of drought conditions are essential for mitigating their impacts and enhancing resilience. Satellite-based drought indices have the advantage of providing spatially continuous and consistent information on drought severity and extent. A new drought product was developed from the thermal infrared observations of the Fengyun (FY) series of satellites. We proposed a data fusion algorithm to combine multiple FY satellites, including FY-2F, FY-2G, and FY-4A, to create a long time series of a land surface temperature (LST) data set without systematic bias. An FY drought index (FYDI) is then derived by coupling the long-term LST data set with the surface–atmospheric energy exchange model at 4 km spatial resolution over China from 2013 to present. The performance and reliability of the new FYDI product are evaluated in this study by comparing it with the Meteorological-drought Composite Index (MCI), one of the authoritative drought monitoring indices used in the Chinese meteorological services. The main objectives of this paper are: (1) to evaluate the performance of the FYDI in capturing the spatiotemporal patterns of drought events over China; (2) to quantitively analyze the consistency between the FYDI and MCI products; and (3) to explore the advantages and limitations of the FYDI for drought monitoring and assessment. The preliminary results show that the FYDI product has good agreement with the MCI, indicating that the FYDI can effectively identify the occurrence, duration, severity, and frequency of drought events over China. These two products have a strong correlation in terms of drought detection, with a correlation coefficient of approximately 0.7. The FYDI was found to be particularly effective in the regions where ground observation is scarce, with the capability of reflecting the spatial heterogeneity and variability of drought patterns more clearly. Overall, the FYDI can be a useful measure for operational drought monitoring and early warning, complementing the existing ground-based MCI drought indices.

Analysis of Long Time Series of Summer Surface Urban Heat Island under the Missing-Filled Satellite Data Scenario.

Surface urban heat islands (SUHIs) are mostly an urban ecological issue. There is a growing demand for the quantification of the SUHI effect, and for its optimization to mitigate the increasing possible hazards caused by SUHI. Satellite-derived land surface temperature (LST) is an important indicator for quantifying SUHIs with frequent coverage. Current LST data with high spatiotemporal resolution is still lacking due to no single satellite sensor that can resolve the trade-off between spatial and temporal resolutions and this greatly limits its applications. To address this issue, we propose a multiscale geographically weighted regression (MGWR) coupling the comprehensive, flexible, spatiotemporal data fusion (CFSDAF) method to generate a high-spatiotemporal-resolution LST dataset. We then analyzed the SUHI intensity (SUHII) in Chengdu City, a typical cloudy and rainy city in China, from 2002 to 2022. Finally, we selected thirteen potential driving factors of SUHIs and analyzed the relation between these thirteen influential drivers and SUHIIs. Results show that: (1) an MGWR outperforms classic methods for downscaling LST, namely geographically weighted regression (GWR) and thermal image sharpening (TsHARP); (2) compared to classic spatiotemporal fusion methods, our method produces more accurate predicted LST images (R2, RMSE, AAD values were in the range of 0.8103 to 0.9476, 1.0601 to 1.4974, 0.8455 to 1.3380); (3) the average summer daytime SUHII increased form 2.08 °C (suburban area as 50% of the urban area) and 2.32 °C (suburban area as 100% of the urban area) in 2002 to 4.93 °C and 5.07 °C, respectively, in 2022 over Chengdu City; and (4) the anthropogenic activity drivers have a higher relative influence on SUHII than other drivers. Therefore, anthropogenic activity driving factors should be considered with CO2 emissions and land use changes for urban planning to mitigate the SUHI effect.

Regional climate trend analyses for Aqua MODIS land surface temperatures

ABSTRACT Land surface temperature (LST) is the skin surface temperature of the Earth’s land surface and is an essential climate variable (ECV) that critically characterizes the Earth’s climate and surface energy processes. This is the first regional trend analysis for LST with uncertainties, using a stable LST climate data record suitable for climate trend analyses: the Aqua Moderate Imaging Spectroradiometer (MODIS) dataset (MYDCCI) produced in the European Space Agency (ESA) LST Climate Change Initiative (LST_cci) Project. The quality of the trends is verified here for Western Europe, for which consistent trends per decade of 0.5 K and 0.57 K across daytime and night-time data, respectively. Eight regions representative of a range of biomes and, including the Amazon, Western U.S.A., Greenland, the Sahel, Siberia, China, India, and Australia, have been analysed. Our most significant findings show substantial LST increases for Siberia and the Amazon, as well as evidence of recent increasing trends for the Western U.S.A. Siberia showed the most substantial daytime and night-time changes per decade of +0.87 K and +0.93 K, respectively. Furthermore, the Amazon and Western U.S.A. showed significant daytime LST trend increases of +0.82 K and +0.51 K, respectively. Western Europe and Western U.S.A. show the largest night-time change of +0.57 K and +0.54 K, respectively. India presents an atypical time series with a decreasing daytime and increasing night-time trend per decade of −0.62 K and +0.44 K, respectively, consistent with reported air temperatures. Overall the results show significant positive regional LST trends, against a background of both strong seasonal cycles and intermittent disruptive events, which demonstrate the importance of continued LST observations. The results also incentivize remote sensing science to derive ever more rigorous LST datasets and to continue to investigate emissivity and error correlations to allow improved trend analyses for smaller regions and individual locations.

Four-decade response of land surface temperature to urban expansion in Beijing

Urban expansion with increased impervious surface area (ISA) can affect the thermal regimes by altering the water and energy balance. Previous studies have indicated that ISA is the primary surrogate of rising land surface temperature (LST). However, urban-induced LST changes are not yet well quantified because of the temporal limitations of the LST datasets and the uncertainties from climate variabilities. In this study, we integrated multiple satellite-based products (vegetation coverage, albedo and radiation forcing) into a physically-based land surface model i.e., Variable Infiltration Capacity (VIC), to detect multi-temporal scale response of LST to the increasing ISA in Beijing. Heat transfer-related parameters in the model were updated along with urbanization processes. The integrated modeling presented high performance for simulating LST, as evaluated with ground-based observations and the MODerate-resolution Imaging Spectroradiometer (MODIS) LST product. We found that the average LST over Beijing increased from 10.3 °C in 1980–1990 to 11.2 °C in 2010–2020, with nearly 33% contributed by the urban growth. The urban-induced thermal effect was particularly strong in summer daytime and winter nighttime. The frequency of heat days (the day with the maximum LST over 40 °C) in a year for Beijing appeared to linearly correlate with the impervious surface fraction (ISF) when ISF > 25%. For the urbanizing area, moreover, the four-decade LST estimates indicated that the annual overheating duration was obviously extended with a rate of about 5 days per decade. We therefore conclude that the urban expansion in Beijing not only amplified heat stress but also altered heat phenology, which have implications for urban planning and the construction of heat mitigation facilities.

Determining the critical threshold of meteorological heat damage to tea plants based on MODIS LST products for tea planting areas in China

Tea (Camellia sinensis (L.) Kuntze) is a daily drink in modern life throughout the world. Global warming increases the exposure of tea plants to elevated heat conditions during the summer picking period. Therefore, it is a fundamental function of the meteorological service to determine the critical temperature threshold that triggers heat damage (HD) to tea plants. In this study, MODIS daytime land surface temperature (LST) products obtained from the Aqua satellite (MYD11A1), annual land cover type products (MCD12Q1), and actual HD records were jointly used to determine the critical LST threshold resulting in HD to tea planting areas in China. An LST dataset of HD and HD-free samples was constructed based on actual HD records. The box-plot method was then used to delineate the range of LST outliers in HD to tea plants. The cumulative frequency distribution (CFD) method was adopted to determine the critical LST threshold by searching for the maximum slope of CFD within the range of LST outliers based on 2901 actual HD samples from 42 agrometeorological stations in 11 provinces in China. Results showed that the critical LST threshold triggering HD to tea plants was 30.1 °C. Two different HD data sources (i.e., 68 counties in eight provinces from the Meteorological Disaster Management System (MDMS) and 37 counties in six cities from a published book) were used as HD validation samples. Validation accuracy analysis revealed that a critical LST threshold of 30.1 °C was reasonable, with an average accuracy of between 60% and 90% for the published-book dataset, and 75.7% in the northern Yangtze River area and 86.3% in the southern Yangtze River area for the MDMS dataset. This study confirmed that using remotely sensed LST products is an effective method for identifying temperature-related agrometeorological disasters owing to the method's ability to intuitively capture the surface temperature distribution of the disaster-affected object. The findings will be helpful for identifying the timing and location of HD to tea plants, and thereby potentially preventing HD to tea plants.

Spatiotemporal change patterns and driving factors of land surface temperature in the Yunnan-Kweichow Plateau from 2000 to 2020

In this study, the spatiotemporal change patterns and driving factors of land surface temperature (LST) on the Yunnan-Kweichow Plateau (YKP) during 2000–2020 are investigated by using the Thermal and Reanalysis Integrating Moderate-resolution Spatial-seamless (TRIMS) LST dataset provided by National Tibetan Plateau Data Center. The YKP LST spatiotemporal change patterns are revealed at annual, seasonal, monthly, and daily scales. Furthermore, seven driving factors such as air temperature, land cover types, normalized difference vegetation index, precipitation, solar radiation, elevation, and latitude are quantified the impacts on LST spatial heterogeneity at annual scale. The main findings are as follows: (1) Annual mean LST increases by 0.016 K/year. Annual mean daytime LST slightly decreases by 0.009 K/year. Annual mean nighttime LST significantly increases by 0.042 K/year. (2) The trend and seasonal components of the daily, daily mean daytime, and daily mean nighttime LST have five and four breakpoints respectively, indicating that the variation of LST is unstable during 2000–2020 on the YKP. (3) The LST lapse rates at nighttime are generally higher than those at daytime on the YKP at the annual, seasonal, and monthly scales. The LST maximum lapse rate is 0.59 K/100 m in summer nighttime, and the LST minimum lapse rate is 0.18 K/100 m in winter daytime. (4) The controlling effects of seven factors are generally stronger in the nighttime than those in the daytime. The factors of elevation and air temperature dominate the LST spatial distribution on the YKP, with a contribution rate of >70 %. In addition, the interactions among the seven factors are all enhancing the effects on the spatial distribution of annual mean LST, including bivariate enhancement and nonlinear enhancement. This study contributes to the mitigation and adaptation to climate change of LST in the plateau and plays a theoretical reference role in formulating corresponding policies for environmental protection.

Adapting an Existing Empirical Algorithm for Microwave Land Surface Temperature Retrieval in China for AMSR2 Data

To extend the time span of the microwave (MW) land surface temperature (LST) dataset in China, this study proposed an optimized empirical algorithm for Advanced Microwave Scanning Radiometer 2 (AMSR2) LST retrieval based on the algorithm for its predecessor, the AMSR-Earth Observing System (AMSR-E). A modified comprehensive classification system of environmental variables (CCSEV) that considered the impact of landform, landcover, atmospheric conditions, and solar radiation on the variation of LST was first constructed, and the LST for each class in the CCSEV was then retrieved through stepwise regression using the brightness temperature in different AMSR2 channels. The results indicate that the annual RMSE of the AMSR2 LST, compared to the reference Moderate Resolution Imaging Spectroradiometer (MODIS) LST from 2012 to 2020, varies between 3.26 K and 3.61 K in the daytime and 2.76 K and 2.96 K in the nighttime, respectively. The RMSE of the AMSR2 LST compared to the field measurements at the sites of the Beidahe river basin and Naqu regions varies between 4.16 K and 5.26 K in the daytime and 2.4 K and 5.17 K in the nighttime. The accuracy is relatively low in the warmer months and daytime due to the stronger solar radiation, and is also relatively low in western China due to the dominate highly fluctuating topography and barren and arid landcover. Generally, the accuracy of the AMSR2 LST is comparable with that of the AMSR-E LST retrieved by the predecessor algorithm, which facilitates coherent long-term applications using AMSR series LST datasets.

AGRO-TOURISM IMPACT ANALYSIS OF CLIMATE CHANGE USING GOOGLE EARTH ENGINE IN THE RAHOVEC WINE REGION OF KOSOVO

Abstract. Climate change is a global problem and has a major impact on the economy of the developing world. The use of open data and modeling provides a considerable resource for governments, municipalities, and NGOs to develop an understanding of how climate change will impact their communities. This paper discusses the statistical concepts used to model MODIS Land Surface Temperature dataset and calculation of water balance with interpretations of the results. The results of this paper will utilize open datasets, an analysis of the impacts of temperature changes on the developing tourist economy in the Republic of Kosovo, and the knowledge of the capacity for leveraging large geographic datasets for open climate change research. The data from this paper shows a high level of vulnerability to viticulture in Kosovo due to significant trends in water balance and land surface temperature trends within the Drini I Bardhë basin and the future impacts of this on the regions wine agro-tourism by understanding the development of climate change in the region.

Noninteger Dimension of Seasonal Land Surface Temperature (LST)

During the few last years, climate change, including global warming, which is attributed to human activities, and its long-term adverse effects on the planet’s functions have been identified as the most challenging discussion topics and have provoked significant concern and effort to find possible solutions. Since the warmth arising from the Earth’s landscapes affects the world’s weather and climate patterns, we decided to study the changes in Land Surface Temperature (LST) patterns in different seasons through nonlinear methods. Here, we particularly wanted to estimate the noninteger dimension and fractal structure of the Land Surface Temperature. For this study, the LST data were obtained during the daytime by a Moderate Resolution Imaging Spectroradiometer (MODIS) on NASA’s Terra satellite. Depending on the time of the year data were collected, temperatures changed in different ranges. Since equatorial regions remain warm, and Antarctica and Greenland remain cold, and also because altitude affects temperature, we selected Riley County in the US state of Kansas, which does not belong to any of these location types, and we observed the seasonal changes in temperature in this county. According to our fractal analysis, the fractal dimension may provide a complexity index to characterize different LST datasets. The multifractal analysis confirmed that the LST data may define a self-organizing system that produces fractal patterns in the structure of data. Thus, the LST data may not only have a wide range of fractal dimensions, but also they are fractal. The results of the present study show that the Land Surface Temperature (LST) belongs to the class of fractal processes with a noninteger dimension. Moreover, self-organized behavior governing the structure of LST data may provide an underlying principle that might be a general outcome of human activities and may shape the Earth’s surface temperature. We explicitly acknowledge the important role of fractal geometry when analyzing and tracing settlement patterns and urbanization dynamics at various scales toward purposeful planning in the development of human settlement patterns.

Evaluation of the Impact of Urban Forestry on Urban Heat Island

The purpose of this study is to statistically evaluate the impact of urban forestry on urban heat Island in Umuahia city, located in Eastern Nigeria using Land Surface Temperature (LST) time series. The LST datasets were extracted from the Moderate Resolution Imaging Spectroradiometer (MODIS, MYD11A2 V6) satellite using Python and JavaScript in Google Earth Engine (GEE) platform. Two microclimate environments were considered namely: The Densely and sparsely populated areas in three different locations each. Each microclimate location had 976 observations starting from the 18th of February 2000 to the 1st of May 2021. The Auto-regressive Moving Average (ARMA) was used to model and predict the time series future points of each location; accuracy was measured using Mean Absolute Percentage Error (MAPE) and Root Mean Square Error (RMSE). The Land surface temperature of each selected microclimate location were monitored and the results showed an average observed temperatures difference of 2.90C between the Densely Populated and Sparsely Populated areas. In conclusion, the study showed that inhabitants of densely populated areas experience higher temperature regime than the inhabitants living within the sparsely populated areas. This increase in temperature within the densely populated could be attributed to the gradual loss of vegetation cover, which is due to man-made activities such as land use, land cover changes, industrial activities, vehicular movements, and the heat generated from buildings within the city. Thus, a spatial development planning strategy, which controls the growth of the built area and add green spaces should be created for a friendly urban environment.

A global historical twice-daily (daytime and nighttime) land surface temperature dataset produced by Advanced Very High Resolution Radiometer observations from 1981 to 2021

Abstract. Land surface temperature (LST) is a key variable for monitoring and evaluating global long-term climate change. However, existing satellite-based twice-daily LST products only date back to 2000, which makes it difficult to obtain robust long-term temperature variations. In this study, we developed the first global historical twice-daily LST dataset (GT-LST), with a spatial resolution of 0.05∘, using Advanced Very High Resolution Radiometer (AVHRR) Level-1b Global Area Coverage (GAC) data from 1981 to 2021. The GT-LST product was generated using four main processes: (1) GAC data reading, calibration, and preprocessing using open-source Python libraries; (2) cloud detection using the AVHRR-Phase I algorithm; (3) land surface emissivity estimation using an improved method considering annual land cover changes; (4) LST retrieval based on a nonlinear generalized split-window algorithm. Validation with in situ measurements from Surface Radiation Budget (SURFRAD) sites and Baseline Surface Radiation Network sites showed that the overall root-mean-square errors (RMSEs) of GT-LST varied from 1.6 to 4.0 K, and nighttime LSTs were typically better than daytime LSTs. Intercomparison with the Moderate Resolution Imaging Spectroradiometer LST products (MYD11A1 and MYD21A1) revealed that the overall root-mean-square difference (RMSD) was approximately 3.0 K. Compared with MYD11A1 LST, GT-LST was overestimated, and relatively large RMSDs were obtained during the daytime, spring, and summer, whereas the significantly smaller positive bias was obtained between GT-LST and MYD21A1 LST. Furthermore, we compared our newly generated dataset with a global AVHRR daytime LST product at the selected measurements of SURFRAD sites (i.e., measurements of these two satellite datasets were valid), which revealed similar accuracies for the two datasets. However, GT-LST can additionally provide nighttime LST, which can be combined with daytime observations estimating relatively accurate monthly mean LST, with an RMSE of 2.7 K. Finally, we compared GT-LST with a regional twice-daily AVHRR LST product over continental Africa in different seasons, with RMSDs ranging from 2.1 to 4.3 K. Considering these advantages, the proposed dataset provides a better data source for a range of research applications. GT-LST is freely available at https://doi.org/10.5281/zenodo.7113080 (1981–2000) (Li et al., 2022a), https://doi.org/10.5281/zenodo.7134158 (2001–2005) (Li et al., 2022b), and https://doi.org/10.5281/zenodo.7813607 (2006–2021) (J. H. Li et al., 2023).

Spatio-temporal Analysis of Environmental Criticality: Planned Versus Unplanned Urbanization

The economic boom in the Indian cities is causing rapid urban growth, mostly in an unplanned fashion. The growth of urban built-up is primarily taking place by replacing vegetation and other low radiative surfaces, increasing the magnitude and spatial extent of heat concentrations. The resulting phenomenal increase in temperature, known as Urban Heat Island (UHI), raises environmental criticality. Satellite remote sensing provides a breakthrough for monitoring the spatiotemporal variations of UHI by estimating Land Surface Temperature (LST) and surface biophysical parameters. The objective of this work is to compare the changing pattern of LST that resulted from urban growth and associated biophysical characteristics in a planned city (Kalyani city) and an unplanned city (Barasat city) in West Bengal, India from 2005 to 2019. Using Landsat data, the study retrieved summer LST along with the prepared vegetation index (NDVI) and built-up index (NDBI) for the years 2005 and 2019. The Environmental Criticality Index (ECI) was calculated for the periods from LST, NDVI, and NDBI datasets. The long-term (1988-2019) LST has been derived using the cloud computation technique to analyze the trend. Over the years, though the average LST of Kalyani is relatively high from Barasat, a rapid increase in LST is noticed for Barasat city. Between 2005 – 2019, the rapid unplanned growth in Barasat city has not only increased the LST but also raised the concern for environmental criticality as compared to Kalyani City. The correlation of LST with NDVI and NDBI suggests that urban heating is significantly controlled by the surface characteristics that need to modify through proper planning for urban sustainability. This study may assist planners, administrators, and researchers in decision-making.

Large-Scale Estimation of Hourly Surface Air Temperature Based on Observations from the FY-4A Geostationary Satellite

Spatially continuous surface air temperature (SAT) is of great significance for various research areas in geospatial communities, and it can be reconstructed by the SAT estimation models that integrate accurate point measurements of SAT at ground sites with wall-to-wall datasets derived from remotely sensed observations of spaceborne instruments. As land surface temperature (LST) strongly correlates with SAT, estimation models are typically developed with LST as a primary input. Geostationary satellites are capable of observing the Earth’s surface across large-scale areas at very high frequencies. Compared to the substantial efforts to estimate SAT at daily or monthly scales using LST derived from MODIS, very limited studies have been performed to estimate SAT at high-temporal scales based on LST from geostationary satellites. Estimation models for hourly SAT based on the LST derived from FY-4A, the first geostationary satellite in China’s new-generation meteorological observation mission, were developed for the first time in this study. The models were fully cross-validated for a very large-scale region with diverse geographic settings using random forest, and specified differently to explore the influence of time and location variables on model performance. Overall predictive performance of the models is about 1.65–2.08 K for sample-based cross-validation, and 2.22–2.70 K for site-based cross-validation. Incorporating time or location variables into the hourly models significantly improves predictive performance, which is also confirmed by the analysis of predictive errors at temporal scales and across sites. The best-performing model with an average RMSE of 2.22 K was utilized for reconstructing maps of SAT for each hour. The hourly models developed in this study have general implications for future studies on large-scale estimating of hourly SAT based on geostationary LST datasets.

A stepwise framework for interpolating land surface temperature under cloudy conditions based on the solar-cloud-satellite geometry

Thermal infrared land surface temperature (LST) data from satellites often contain extensive missing values due to high cloudiness degree, which severely hinders their use in applications. Despite the many methods developed, common methods, such as fusion with microwave or reanalysis data and the surface energy budget approach, still remain subject to important limitations and uncertainties, such as dependency on coarse resolution data and difficulty in interpolation for large-scale missing LST data. In this study, we proposed a stepwise framework for estimating missing cloudy-sky LST values of Moderate Resolution Imaging Spectroradiometer (MODIS) from informative samples owing to the solar-cloud-satellite geometry (SCSG) effect, by which satellite imagery records the cloudy-sky LST values of a portion of pixels. We first estimated the clear-sky LST equivalents for all cloud-affected pixels via a similarity-based approach and then determined unknown LSTs for cloudy pixels by training a machine-learning model on cloudy-sky LST values observed owing to the SCSG effect. We demonstrated the utility of this approach by using MODIS/Aqua daytime LST data over Qinghai-Tibet Plateau (QTP) and validated the interpolation results against representative in-situ LST observations and two recently published all-weather LST datasets. When compared to the corresponding in-situ measurements, the interpolated cloudy-sky LST values showed satisfactory accuracy with a mean absolute error (MAE) value of 3.99 °C and a coefficient of determination (R2) value of 0.74, while the MODIS/Aqua clear-sky LST values led to an MAE value of 2.66 °C and an R2 value of 0.86. Compared to the two all-weather LST datasets, results of this study showed the highest accuracy over the data-gap-filled regions in terms of all quantitative performance metrics, more natural transition textures, and better representation of seasonal characteristics. The proposed framework has the advantage of relying on the MODIS family data and handling extensive missing data as well as triggers opportunities to leverage the SCSG effect to produce high-quality all-weather LST data.

Global hourly, 5 km, all-sky land surface temperature data from 2011 to 2021 based on integrating geostationary and polar-orbiting satellite data

Abstract. Land surface temperature (LST) plays a dominant role in the surface energy budget (SEB) and hydrological cycling. Thermal infrared (TIR) remote sensing is the primary method of estimating LST globally. However, cloud cover leaves numerous data gaps in satellite LST products, which seriously restricts their applications. Efforts have been made to produce gap-free LST products from polar-orbiting satellites (e.g., Terra and Aqua); however, satellite data from limited overpasses are not suitable for characterizing the diurnal temperature cycle (DTC), which is directly related to heat waves, plant water stress, and soil moisture. Considering the high temporal variability in LST and the importance of the DTC, we refined the SEB-based cloudy-sky LST recovery method by improving its feasibility and efficiency and produced a global hourly, 5 km, all-sky land surface temperature (GHA-LST) dataset from 2011 to 2021. The GHA-LST product was generated using TIR LST products from geostationary and polar-orbiting satellite data from the Copernicus Global Land Service (CGLS) and the Moderate Resolution Imaging Spectroradiometer (MODIS). Based on ground measurements at the 201 global sites from the Surface Radiation Budget (SURFRAD), Baseline Surface Radiation Network (BSRN), Fluxnet, AmeriFlux, Heihe River basin (HRB), and Tibetan Plateau (TP) networks, the overall root-mean-square error (RMSE) of the hourly GHA-LST product was 3.31 K, with a bias of −0.57 K and R2 of 0.95. Thus, this product was more accurate than the clear-sky CGLS and MODIS MYD21C1 LST samples. The RMSE value of the daily mean LST was 1.76 K. Validation results at individual sites indicate that the GHA-LST dataset has relatively larger RMSEs for high-elevation regions, which can be attributed to high surface heterogeneity and input data uncertainty. Temporal and spatial analyses suggested that GHA-LST has satisfactory spatiotemporal continuity and reasonable variation and matches the reference data well at hourly and daily scales. Furthermore, the regional comparison of GHA-LST with other gap-free hourly datasets (ERA5 and Global Land Data Assimilation System, GLDAS) demonstrated that GHA-LST can provide more spatial texture information. The monthly anomaly analysis suggests that GHA-LST couples well with global surface air temperature datasets and other LST datasets at daily mean and minimum temperature scales, whereas the maximum temperature and diurnal temperature range of LST and air temperature (AT) have different anomalous magnitudes. The GHA-LST dataset is the first global gap-free LST dataset at an hourly, 5 km scale with high accuracy, and it can be used to estimate global evapotranspiration, monitor extreme weather, and advance meteorological forecasting models. GHA-LST is freely available at https://doi.org/10.5281/zenodo.7487284 (Jia et al., 2022b) and http://glass.umd.edu/allsky_LST/GHA-LST (last access: 10 February 2023; Jia et al., 2022c).