Estimation of All-Weather 1 km MODIS Land Surface Temperature for Humid Summer Days

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.

All-sky 1 km MODIS Land Surface Temperature Reconstruction Considering Cloud Effects Based on Machine Learning

A high spatio-temporal resolution land surface temperature (LST) is necessary for various research fields because LST plays a crucial role in the energy exchange between the atmosphere and the ground surface. The moderate-resolution imaging spectroradiometer (MODIS) LST has been widely used, but it is not available under cloudy conditions. This study proposed a novel approach for reconstructing all-sky 1 km MODIS LST in South Korea during the summer seasons using various data sources, considering the cloud effects on LST. In South Korea, a Local Data Assimilation and Prediction System (LDAPS) with a relatively high spatial resolution of 1.5 km has been operated since 2013. The LDAPS model’s analysis data, binary MODIS cloud cover, and auxiliary data were used as input variables, while MODIS LST and cloudy-sky in situ LST were used together as target variables based on the light gradient boosting machine (LightGBM) approach. As a result of spatial five-fold cross-validation using MODIS LST, the proposed model had a coefficient of determination (R2) of 0.89–0.91 with a root mean square error (RMSE) of 1.11–1.39 °C during the daytime, and an R2 of 0.96–0.97 with an RMSE of 0.59–0.60 °C at nighttime. In addition, the reconstructed LST under the cloud was evaluated using leave-one-station-out cross-validation (LOSOCV) using 22 weather stations. From the LOSOCV results under cloudy conditions, the proposed LightGBM model had an R2 of 0.55–0.63 with an RMSE of 2.41–3.00 °C during the daytime, and an R2 of 0.70–0.74 with an RMSE of 1.31–1.36 °C at nighttime. These results indicated that the reconstructed LST has higher accuracy than the LDAPS model. This study also demonstrated that cloud cover information improved the cloudy-sky LST estimation accuracy by adequately reflecting the heterogeneity of the relationship between LST and input variables under clear and cloudy skies. The reconstructed all-sky LST can be used in a variety of research applications including weather monitoring and forecasting.

Global spatiotemporally continuous MODIS land surface temperature dataset

Land surface temperature (LST) plays a critical role in land surface processes. However, as one of the effective means for obtaining global LST observations, remote sensing observations are inherently affected by cloud cover, resulting in varying degrees of missing data in satellite-derived LST products. Here, we propose a solution. First, the data interpolating empirical orthogonal functions (DINEOF) method is used to reconstruct invalid LSTs in cloud-contaminated areas into ideal, clear-sky LSTs. Then, a cumulative distribution function (CDF) matching-based method is developed to correct the ideal, clear-sky LSTs to the real LSTs. Experimental results prove that this method can effectively reconstruct missing LST data and guarantee acceptable accuracy in most regions of the world, with RMSEs of 1–2 K and R values of 0.820–0.996 under ideal, clear-sky conditions and RMSEs of 4–7 K and R values of 0.811–0.933 under all weather conditions. Finally, a spatiotemporally continuous MODIS LST dataset at 0.05° latitude/longitude grids is produced based on the above method.

Measurement(s)	land surface temperature
Technology Type(s)	satellite imaging
Sample Characteristic - Environment	planetary surface
Sample Characteristic - Location	global

Reconstruction of All-Weather Daytime and Nighttime MODIS Aqua-Terra Land Surface Temperature Products Using an XGBoost Approach

Generating spatiotemporally continuous land surface temperature (LST) data is in great demand for hydrology, meteorology, ecology, environmental studies, etc. However, the thermal infrared (TIR)-based LST measurements are prone to cloud contamination with missing pixels. To repair the missing pixels, a new XGBoost-based linking approach for reconstructing daytime and nighttime Moderate Resolution Imaging Spectroradiometer (MODIS) LST measurements was introduced. The instantaneous solar radiation and two soil-related predictors from China Data Assimilation System (CLDAS) 0.0625°/1-h data were selected as the linking variables to depict the relationship with instantaneous MODIS LST data. Other land surface properties, including two vegetation indices, the water index, the surface albedo, and topographic parameters, were also used as the predictor variables. The XGBoost method was used to fit an LST linking model by the training datasets from clear-sky pixels and was then applied to the MODIS Aqua-Terra LSTs during summer time (June to August) in 2017 and 2018 across China. The recovered LST data was further rectified with the Savitzky–Golay (SG) filtering method. The results showed the distribution of the reconstructed LSTs present a reasonable pattern for different land-cover types and topography. The evaluation results using in situ longwave radiation measurements showed the RMSE varies from 3.91 K to 5.53 K for the cloud-free pixels and from 4.42 K to 4.97 K for the cloud-covered pixels. In addition, the reconstructed LST products correlated well with CLDAS LST data with similar LST spatial patterns. The variable importance analysis revealed that the two soil-related predictors and the elevation variable are key parameters due to their great contribution to the XGBoost model performance.

Retrieval of All-Weather 1 km Land Surface Temperature from Combined MODIS and AMSR2 Data over the Tibetan Plateau

Land surface temperature (LST) is one of the most valuable variables for applications relating to hydrological processes, drought monitoring and climate change. LST from satellite data provides consistent estimates over large scales but is only available for cloud-free pixels, greatly limiting applications over frequently cloud-covered regions. With this study, we propose a method for estimating all-weather 1 km LST by combining passive microwave and thermal infrared data. The product is based on clear-sky LST retrieved from Moderate-resolution Imaging Spectroradiometer (MODIS) thermal infrared measurements complemented by LST estimated from the Advanced Microwave Scanning Radiometer Version 2 (AMSR2) brightness temperature to fill gaps caused by clouds. Terrain, vegetation conditions, and AMSR2 multiband information were selected as the auxiliary variables. The random forest algorithm was used to establish the non-linear relationship between the auxiliary variables and LST over the Tibetan Plateau. To assess the error of this method, we performed a validation experiment using clear-sky MODIS LST and in situ measurements. The estimated all-weather LST approximated MODIS LST with an acceptable error, with a coefficient of correlation (r) between 0.87 and 0.99 and a root mean square error (RMSE) between 2.24 K and 5.35 K during the day. At night-time, r was between 0.89 and 0.99 and the RMSE was between 1.02 K and 3.39 K. The error between the estimated LST and in situ LST was also found to be acceptable, with the RMSE for cloudy pixels between 5.15 K and 6.99 K. This method reveals a significant potential to derive all-weather 1 km LST using AMSR2 and MODIS data at a regional and global scale, which will be explored in the future.

An Improved Cloud Gap-Filling Method for Longwave Infrared Land Surface Temperatures through Introducing Passive Microwave Techniques

Satellite-derived land surface temperature (LST) data are most commonly observed in the longwave infrared (LWIR) spectral region. However, such data suffer frequent gaps in coverage caused by cloud cover. Filling these ‘cloud gaps’ usually relies on statistical re-constructions using proximal clear sky LST pixels, whilst this is often a poor surrogate for shadowed LSTs insulated under cloud. Another solution is to rely on passive microwave (PM) LST data that are largely unimpeded by cloud cover impacts, the quality of which, however, is limited by the very coarse spatial resolution typical of PM signals. Here, we combine aspects of these two approaches to fill cloud gaps in the LWIR-derived LST record, using Kenya (East Africa) as our study area. The proposed “cloud gap-filling” approach increases the coverage of daily Aqua MODIS LST data over Kenya from <50% to >90%. Evaluations were made against the in situ and SEVIRI-derived LST data respectively, revealing root mean square errors (RMSEs) of 2.6 K and 3.6 K for the proposed method by mid-day, compared with RMSEs of 4.3 K and 6.7 K for the conventional proximal-pixel-based statistical re-construction method. We also find that such accuracy improvements become increasingly apparent when the total cloud cover residence time increases in the morning-to-noon time frame. At mid-night, cloud gap-filling performance is also better for the proposed method, though the RMSE improvement is far smaller (<0.3 K) than in the mid-day period. The results indicate that our proposed two-step cloud gap-filling method can improve upon performances achieved by conventional methods for cloud gap-filling and has the potential to be scaled up to provide data at continental or global scales as it does not rely on locality-specific knowledge or datasets.

A Review of Reconstructing Remotely Sensed Land Surface Temperature under Cloudy Conditions

Land surface temperature (LST) is an important environmental parameter in climate change, urban heat islands, drought, public health, and other fields. Thermal infrared (TIR) remote sensing is the main method used to obtain LST information over large spatial scales. However, cloud cover results in many data gaps in remotely sensed LST datasets, greatly limiting their practical applications. Many studies have sought to fill these data gaps and reconstruct cloud-free LST datasets over the last few decades. This paper reviews the progress of LST reconstruction research. A bibliometric analysis is conducted to provide a brief overview of the papers published in this field. The existing reconstruction algorithms can be grouped into five categories: spatial gap-filling methods, temporal gap-filling methods, spatiotemporal gap-filling methods, multi-source fusion-based gap-filling methods, and surface energy balance-based gap-filling methods. The principles, advantages, and limitations of these methods are described and discussed. The applications of these methods are also outlined. In addition, the validation of filled LST values’ cloudy pixels is an important concern in LST reconstruction. The different validation methods applied for reconstructed LST datasets are also reviewed herein. Finally, prospects for future developments in LST reconstruction are provided.

Modified Linear Scaling and Quantile Mapping Mean Bias Correction of MODIS Land Surface Temperature for Surface Air Temperature Estimation for the Lowland Areas of Peninsular Malaysia

MODIS land surface temperature data (MODIS Ts) products are quantified from the earth surface’s reflected thermal infrared signal via sensors onboard the Terra and Aqua satellites. MODIS Ts products are a great value to many environmental applications but often subject to discrepancies when compared to the air temperature (Ta) data that represent the temperature measured at 2 m above the ground surface. Although they are different in their nature, the relationship between Ts and Ta has been established by many researchers. Further validation and correction on the relationship between these two has enabled the estimation of Ta from MODIS Ts products in order to overcome the limitation of Ta that can only provide data in a point form with a very limited area coverage. Therefore, this study was conducted with the objective to assess the accuracy of MODIS Ts products, i.e., MOD11A1, MOD11A2, MYD11A1, and MYD11A2 against Ta and to identify the performance of a modified Linear Scaling using a constant and monthly correction factor (LS-MBC), and Quantile Mapping Mean Bias Correction (QM-MBC) methods for lowland area of Peninsular Malaysia. Furthermore, the correction factor (CF) values for each MBC were adjusted according to the condition set depending on the different bias levels. Then, the performance of the pre- and post-MBC correction for by stations and regions analysis were evaluated through root mean square error (RMSE), percentage bias (PBIAS), mean absolute error (MAE), and correlation coefficient (r). The region dataset is obtained by stacking the air temperature (Ta_r) and surface temperature (Ts_r) data corresponding to the number of stations within the identified regions. The assessment of pre-MBC data for both 36 stations and 5 regions demonstrated poor correspondence with high average errors and percentage biases, i.e., RMSE = 3.33–5.42 °C, PBIAS = 1.36–12.07%, MAE = 2.88–4.89 °C, and r = 0.16–0.29. The application of the MBCs has successfully reduced the errors and bias percentages, and slightly increased the r values for all MODIS Ts products. All post-MBC depicted good average accuracies (RMSE and MAE < 3 °C and PBIAS between ±5%) and r between 0.18 and 0.31. In detail, for the station analysis, the LS-MBC using monthly CF recorded better performance than the LS-MBC using constant CF or the QM-MBC. For the regional study, the QM-MBC outperformed the others. This study illustrated that the proposed LS-MBC, in spite of its simplicity, managed to perform well in reducing the error and bias terms of MODIS Ts as much as the performance of the more complex QM-MBC method.

A Two-Step Integrated MLP-GTWR Method to Estimate 1 km Land Surface Temperature with Complete Spatial Coverage in Humid, Cloudy Regions

There is an increasing demand for a land surface temperature (LST) dataset with both fine spatial and temporal resolutions due to the key role of LST in the Earth’s land–atmosphere system. Currently, the technique most commonly used to meet the demand is thermal infrared (TIR) remote sensing. However, cloud contamination interferes with TIR transmission through the atmosphere, limiting the potential of space-borne TIR sensors to provide the LST with complete spatio-temporal coverage. To solve this problem, we developed a two-step integrated method to: (i) estimate the 10-km LST with a high spatial coverage from passive microwave (PMW) data using the multilayer perceptron (MLP) model; and (ii) downscale the LST to 1 km and fill the gaps based on the geographically and temporally weighted regression (GTWR) model. Finally, the 1-km all-weather LST for cloudy pixels was fused with Aqua MODIS clear-sky LST via bias correction. This method was applied to produce the all-weather LST products for both daytime and nighttime during the years 2013–2018 in South China. The evaluations showed that the accuracy of the reproduced LST on cloudy days was comparable to that of the MODIS LST in terms of mean absolute error (2.29–2.65 K), root mean square error (2.92–3.25 K), and coefficients of determination (0.82–0.92) against the in situ measurements at four flux stations and ten automatic meteorological stations with various land cover types. The spatial and temporal analysis showed that the MLP-GTWR LST were highly consistent with the MODIS, in situ, and ERA5-Land LST, with the satisfactory ability to present the LST pattern under cloudy conditions. In addition, the MLP-GTWR method outperformed a gap-filling method and another TIR-PMW integrated method due to the local strategy in MLP and the consideration of temporal non-stationarity relationship in GTWR. Therefore, the test of the developed method in the frequently cloudy South China indicates the efficient potential for further application to other humid regions to generate the LST under cloudy condition.

The Gradient Effect on the Relationship between the Underlying Factor and Land Surface Temperature in Large Urbanized Region

Although research relating to the urban heat island (UHI) phenomenon has been significantly increasing in recent years, there is still a lack of a continuous and clear recognition of the potential gradient effect on the UHI—landscape relationship within large urbanized regions. In this study, we chose the Beijing-Tianjin-Hebei (BTH) region, which is a large scaled urban agglomeration in China, as the case study area. We examined the causal relationship between the LST variation and underlying surface characteristics using multi-temporal land cover and summer average land surface temperature (LST) data as the analyzed variables. This study then further discussed the modeling performance when quantifying their relationship from a spatial gradient perspective (the grid size ranged from 6 to 24 km), by comparing the ordinary least squares (OLS) and geographically weighted regression (GWR) methods. The results indicate that: (1) both the OLS and GWR analysis confirmed that the composition of built-up land contributes as an essential factor that is responsible for the UHI phenomenon in a large urban agglomeration region; (2) for the OLS, the modeled relationship between the LST and its drive factor showed a significant spatial gradient effect, changing with different spatial analysis grids; and, (3) in contrast, using the GWR model revealed a considerably robust and better performance for accommodating the spatial non-stationarity with a lower scale dependence than that of the OLS model. This study highlights the significant spatial heterogeneity that is related to the UHI effect in large-extent urban agglomeration areas, and it suggests that the potential gradient effect and uncertainty induced by different spatial scale and methodology usage should be considered when modeling the UHI effect with urbanization. This would supplement current UHI study and be beneficial for deepening the cognition and enlightenment of landscape planning for UHI regulation.

The Integration of Remote Sensing and Field Surveys to Detect Ecologically Damaged Areas for Restoration in South Korea

Ecological damage refers to the reduction in the value of the environment due to human activities such as development. The intensity of ecosystem damage is worsening worldwide. Although the importance of restoration projects to reduce ecosystem damage is increasing, they are difficult to carry out, owing to the absence of data and monitoring of damaged areas. In this study, ecologically damaged areas for restoration in South Korea were detected using remote sensing and field surveys. For the analysis, national standardized vector datasets and Google Earth images were used; field surveys were conducted from 2018 to early 2020. Our results showed that 62% of the ecological damage that occurred in South Korea existed in forest ecosystems; the damaged areas were mostly smaller than 50,000 m2. Additionally, most of the causes and types of damage due to human activities such as development were soil erosion related. The results also suggest the importance of obtaining monitoring data on ecologically damaged areas and the importance of establishing an appropriate restoration plan using this data.

Land Surface Temperature Retrieval from Passive Microwave Satellite Observations: State-of-the-Art and Future Directions

Land surface temperature (LST) is an important variable in the physics of land–surface processes controlling the heat and water fluxes over the interface between the Earth’s surface and the atmosphere. Space-borne remote sensing provides the only feasible way for acquiring high-precision LST at temporal and spatial domain over the entire globe. Passive microwave (PMW) satellite observations have the capability to penetrate through clouds and can provide data under both clear and cloud conditions. Nonetheless, compared with thermal infrared data, PMW data suffer from lower spatial resolution and LST retrieval accuracy. Various methods for estimating LST from PMW satellite observations were proposed in the past few decades. This paper provides an extensive overview of these methods. We first present the theoretical basis for retrieving LST from PMW observations and then review the existing LST retrieval methods. These methods are mainly categorized into four types, i.e., empirical methods, semi-empirical methods, physically-based methods, and neural network methods. Advantages, limitations, and assumptions associated with each method are discussed. Prospects for future development to improve the performance of LST retrieval methods from PMW satellite observations are also recommended.