Evaluation of Land Surface Temperature Retrieval from FY-3B/VIRR Data in an Arid Area of Northwestern China

Land Surface Temperature Retrieval from Fengyun-3D Medium Resolution Spectral Imager II (FY-3D MERSI-II) Data with the Improved Two-Factor Split-Window Algorithm

Land surface temperature (LST) is an essential parameter widely used in environmental studies. The Medium Resolution Spectral Imager II (MERSI-II) boarded on the second generation Chinese polar-orbiting meteorological satellite, Fengyun-3D (FY-3D), provides a new opportunity for LST retrieval at a spatial resolution of 250 m that is higher than that of the already widely used Moderate Resolution Imaging Spectrometer (MODIS) LST data of 1000 m. However, there is no operational LST product from FY-3D MERSI-II data available for free access. Therefore, in this study, we developed an improved two-factor split-window algorithm (TFSWA) of LST retrieval from this data source as it has two thermal-infrared (TIR) bands. The essential coefficients of the TFSWA algorithm have been carefully and precisely estimated for the FY-3D MERSI-II TIR thermal bands. A new approach for estimating land surface emissivity has been developed using the ASTER Global Emissivity Database (ASTER GED) and the International Geosphere-Biosphere Program (IGBP) data. A model to estimate the atmospheric water vapor content (AWVC) from the three atmospheric water vapor absorption bands (bands 16, 17, and 18) has been developed as AWVC has been recognized as the most important factor determining the variation of AT. Using MODTRAN 5.2, the equations for the AT estimate from the retrieved AWVC were established. In addition, the AT of the pixels at the far edge of FY-3D MERSI-II data may be strongly affected by the increase of the optical path. Viewing zenith angle (VZA) correction equations were proposed in the study to correct this effect on AT estimation. Field data from four stations were applied to validate the improved TFSWA in the study. Cross-validation with MODIS LST (MYD11) was also conducted to evaluate the improved TFSWA. The cross-validation result indicates that the FY-3D MERSI-II LST from the improved TFSWA are comparable with MODIS LST while the correlation coefficients between FY-3D MERSI-II LST and MODIS LST over the Mid-East China region are in the range of 0.84~0.98 for different seasons and land cover types. Validation with 318 field LST samples indicates that the average MAE and R2 of the scenes at the four stations are about 1.97 K and 0.98, respectively. Thus, it could be concluded that the improved TFSWA developed in the study can be a good algorithm for LST retrieval from FY-3D MERSI-II data with acceptable accuracy.

Spatial Downscaling of Land Surface Temperature over Heterogeneous Regions Using Random Forest Regression Considering Spatial Features

Land surface temperature (LST) is one of the crucial parameters in the physical processes of the Earth. Acquiring LST images with high spatial and temporal resolutions is currently difficult because of the technical restriction of satellite thermal infrared sensors. Downscaling LST from coarse to fine spatial resolution is an effective means to alleviate this problem. A spatial random forest downscaling LST method (SRFD) was proposed in this study. Abundant predictor variables—including land surface reflection data, remote sensing spectral indexes, terrain factors, and land cover type data—were considered and applied for feature selection in SRFD. Moreover, the shortcoming of only focusing on information from point-to-point in previous statistics-based downscaling methods was supplemented by adding the spatial feature of LST. SRFD was applied to three different heterogeneous regions and compared with the results from three classical or excellent methods, including thermal image sharpening algorithm, multifactor geographically weighted regression, and random forest downscaling method. Results show that SRFD outperforms other methods in vision and statistics due to the benefits from the supplement of the LST spatial feature. Specifically, compared with RFD, the second-best method, the downscaling results of SRFD are 10% to 24% lower in root-mean-square error, 5% to 20% higher in the coefficient of determination, 11% to 25% lower in mean absolute error, and 4% to 17% higher in structural similarity index measure. Hence, we conclude that SRFD will be a promising LST downscaling method.

An Operational Split-Window Algorithm for Retrieving Land Surface Temperature from Geostationary Satellite Data: A Case Study on Himawari-8 AHI Data

An operational split-window (SW) algorithm was developed to retrieve high-temporal-resolution land surface temperature (LST) from global geostationary (GEO) satellite data. First, the MODTRAN 5.2 and SeeBor V5.0 atmospheric profiles were used to establish a simulation database to derive the SW algorithm coefficients for GEO satellites. Then, the dynamic land surface emissivities (LSEs) in the two SW bands were estimated using the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) Global Emissivity Dataset (GED), fractional vegetation cover (FVC), and snow cover products. Here, the proposed SW algorithm was applied to Himawari-8 Advanced Himawari Imager (AHI) observations. LST estimates were retrieved in January, April, July, and October 2016, and three validation methods were used to evaluate the LST retrievals, including the temperature-based (T-based) method, radiance-based (R-based) method, and intercomparison method. The in situ night-time observations from two Heihe Watershed Allied Telemetry Experimental Research (HiWATER) sites and four Terrestrial Ecosystem Research Network (TERN) OzFlux sites were used in the T-based validation, where a mean bias of −0.70 K and a mean root-mean-square error (RMSE) of 2.29 K were achieved. In the R-based validation, the biases were 0.14 and −0.13 K and RMSEs were 0.83 and 0.86 K for the daytime and nighttime, respectively, over four forest sites, four desert sites, and two inland water sites. Additionally, the AHI LST estimates were compared with the Collection 6 MYD11_L2 and MYD21_L2 LST products over southeastern China and the Australian continent, and the results indicated that the AHI LST was more consistent with the MYD21 LST and was generally higher than the MYD11 LST. The pronounced discrepancy between the AHI and MYD11 LST could be mainly caused by the differences in the emissivities used. We conclude that the developed SW algorithm is of high accuracy and shows promise in producing LST data with global coverage using observations from a constellation of GEO satellites.

Improvement of Split-Window Algorithm for Land Surface Temperature Retrieval from Sentinel-3A SLSTR Data Over Barren Surfaces Using ASTER GED Product

Land surface temperature (LST) is a key variable influencing the energy balance between the land surface and the atmosphere. In this work, a split-window algorithm was used to calculate LST from Sentinel-3A Sea and Land Surface Temperature Radiometer (SLSTR) thermal infrared data. The National Centers for Environmental Prediction (NCEP) reanalysis atmospheric profiles combined with the radiation transport model MODerate resolution atmospheric TRANsmission version 5.2 (MODTRAN 5.2) were utilized to obtain atmospheric water vapor content (WVC). The ASTER Global Emissivity Database Version 3 (ASTER GED v3) product was utilized to estimate surface emissivity in order to improve the accuracy of LST estimation over barren surfaces. Using a simulation database, the coefficients of the algorithm were fitted and the performance of the algorithm was evaluated. The root-mean-square error (RMSE) values of the differences between the estimated LST and the actual LST of the MODTRAN radiative transfer simulation at each WVC subrange of 0–6.5 g/cm2 were less than 1.0 K. To validate the retrieval accuracy, ground-based LST measurements were collected at two relatively homogeneous desert study sites in Dalad Banner and Wuhai, Inner Mongolia, China. The bias between the retrieved LST and the in situ LST was about 0.2 K and the RMSE was about 1.3 K at the Dalad Banner site, whereas they were approximately -0.4 and 1.0 K at the Wuhai site. As a reference, the retrieved LST was compared with the operational SLSTR LST product in this study. The bias between the SLSTR LST product and the in situ LST was approximately 1 K and the RMSE was approximately 2 K at the Dalad Banner site, whereas they were approximately 1.1 and 1.4 K at the Wuhai site. The results demonstrate that the split-window algorithm combined with improved emissivity estimation based on the ASTER GED product can distinctly obtain better accuracy of LST over barren surfaces.

A Split Window Algorithm for Retrieving Land Surface Temperature from FY-3D MERSI-2 Data

The thermal infrared (TIR) data from the Medium Resolution Spectral Imager II (MERSI-2) on the Chinese meteorological satellite FY-3D have high spatiotemporal resolution. Although the MERSI-2 land surface temperature (LST) products have good application prospects, there are some deviations in the TIR band radiance from MERSI-2. To accurately retrieve LSTs from MERSI-2, a method based on a cross-calibration model and split window (SW) algorithm is proposed. The method is divided into two parts: cross-calibration and LST retrieval. First, the MODTRAN program is used to simulate the radiation transfer process to obtain MERSI-2 and Moderate Resolution Imaging Spectroradiometer (MODIS) simulation data, establish a cross-calibration model, and then calculate the actual brightness temperature (BT) of the MERSI-2 image. Second, according to the characteristics of the near-infrared (NIR) bands, the atmospheric water vapor content (WVC) is retrieved, and the atmospheric transmittance is calculated. The land surface emissivity is estimated by the NDVI-based threshold method, which ensures that both parameters (transmittance and emissivity) can be acquired simultaneously. The validation shows the following: 1) The average accuracy of our algorithm is 0.42 K when using simulation data; 2) the relative error of our algorithm is 1.37 K when compared with the MODIS LST product (MYD11A1); 3) when compared with ground-measured data, the accuracy of our algorithm is 1.23 K. Sensitivity analysis shows that the SW algorithm is not sensitive to the two main parameters (WVC and emissivity), which also proves that the estimation of LST from MERSI-2 data is feasible. In general, our algorithm exhibits good accuracy and applicability, but it still requires further improvement.

Estimating Land Surface Temperature from Landsat-8 Data using the NOAA JPSS Enterprise Algorithm

Land surface temperature (LST) is one of the key parameters in hydrology, meteorology, and the surface energy balance. The National Oceanic and Atmospheric Administration (NOAA) Joint Polar Satellite System (JPSS) Enterprise algorithm is adapted to Landsat-8 data to obtain the estimate of LST. The coefficients of the Enterprise algorithm were obtained by linear regression using the analog data produced by comprehensive radiative transfer modeling. The performance of the Enterprise algorithm was first tested by simulation data and then validated by ground measurements. In addition, the accuracy of the Enterprise algorithm was compared to the generalized split-window algorithm and the split-window algorithm of Sobrino et al. (1996). The validation results indicate the Enterprise algorithm has a comparable accuracy to the other two split-window algorithms. The biases (root mean square errors) of the Enterprise algorithm were 1.38 (3.22), 1.01 (2.32), 1.99 (3.49), 2.53 (3.46), and −0.15 K (1.11 K) at the SURFRAD, HiWATER_A, HiWATER_B, HiWATER_C sites and BanGe site, respectively, whereas those values were 1.39 (3.20), 1.0 (2.30), 1.93 (3.48), 2.53 (3.35), and −0.35 K (1.16 K) for the generalized split-window algorithm, 1.45 (3.39), 1.08 (2.41), 2.16 (3.67), 2.52 (3.58), and 0.02 K (1.12 K) for the split-window algorithm of Sobrino, respectively. This study provides an alternative method to estimate LST from Landsat-8 data.