Generalized Split-window Algorithm Research Articles

Retrieval and validation of the dawn-dusk land surface temperature from FY-3E MERSI-LL

ABSTRACT Land surface temperature (LST) is an important physical parameter that reflects changes in surface processes. Remote sensing makes it possible to obtain global or regional LSTs with high spatial and temporal resolution. Feng-Yun 3E (FY-3E) is the world’s first civilian dawn-dusk orbiting meteorological satellite on board the Medium Resolution Spectral Imager – Low-Light (MERSI-LL) sensor capable of acquiring LST at 250 m spatial resolution at dawn and dusk. By complementing with existing polar-orbiting satellites, FY-3E can fill the satellite observation gap within the 6-hour assimilation window and further enhance the observation capability of LST. To retrieve accurate LST from MERSI-LL, the generalized split-window (GSW) method was developed. The GSW coefficients were derived from the simulation dataset produced by Moderate Resolution Transmittance Code 5.2 with an atmospheric profile database. As a key input to the GSW algorithm, the land surface emissivity was dynamically estimated by the Advanced Spaceborne Thermal Emission and Reflection Radiometer Global Emissivity Dataset in combination with the fractional vegetation cover. The water vapour content is calculated in real time by the split-window covariance/variance ratio method. The comparison between retrieved LSTs and in situ LSTs shows that the GSW algorithm can effectively retrieve the LSTs of vegetation cover surfaces. The root mean square error ranged from 0.79 to 1.79 K at dawn and from 1.32 to 1.92 K at dusk.

A Physical-Based Method for Pixel-by-Pixel Quantifying Uncertainty of Land Surface Temperature Retrieval From Satellite Thermal Infrared Data Using the Generalized Split-Window Algorithm

Land surface temperature (LST) is an important physical parameter at the interface between the Earth’s surface and the atmosphere. Accurately quantifying LST uncertainty is essential for the generation of a long-term and consistent LST Climate Data Record (CDR) or Earth System Data Record (ESDR) from either multiple sensors or algorithms. In this study, a physical-based method was proposed to quantify the uncertainty of LST retrieval from satellite thermal infrared data using the generalized split-window (GSW) algorithm. LST uncertainties were parameterized as a function of brightness temperature at the top of the atmosphere (TOA) and surface emissivity in two split-window channels, which are two key input parameters in the GSW algorithm, as well as their uncertainties. The performance of the parameterized uncertainty model was evaluated according to simulation dataset at six prescribed viewing zenith angles (VZAs) of 0°, 33.56°, 44.42°, 51.32°, 56.25°, and 60°, with a root mean squared error (RMSE) of 0.001 K. The coefficients of the parameterized uncertainty model at arbitrary VZA within a sensor’ field of view (FOV) can be obtained by linear interpolation of the coefficients at the six prescribed VZAs. Once the coefficients of the parameterized uncertainty model for each pixel are available, total LST uncertainties can be quantified on a pixel-by-pixel basis. As an example, the parameterized uncertainty model was applied to actual MODIS data for displaying the spatial distribution of LST uncertainties. The results indicate that the parameterized uncertainty model can well characterize the spatial variation in LST uncertainties over various land cover types.

A Thermal Infrared Land Surface Temperature Retrieval Algorithm for Thin Cirrus Skies Using Cirrus Optical Properties

To acquire daytime land surface temperature (LST) in thin cirrus cloudy skies, we have developed a three-channel LST retrieval algorithm based on a widely used two-channel LST retrieval algorithm for the clear-sky conditions. In this algorithm, the LST is expressed as a multiple linear function of MODIS channels 29, 31 and 32 with the coefficients of the linear function dependent on the cirrus optical depth (COD) and cirrus effective radius (R). The influences from land surface emissivities (LSEs) are also considered in this algorithm. The simulated dataset shows that the LST could be estimated using the proposed algorithm with the root mean squire error (RMSE) less than 2.2 K in thin cirrus cloudy skies (COD less than 0.7) when viewing zenith angle (VZA) equivalent to 0°. As VZA is equivalent to 60°, the maximum RMSE are 2.7 K. The widely used generalized split-window (GSW) algorithm proposed for clear-sky conditions are used in cirrus cloudy skies, and the RMSEs of GSW algorithm estimated LST are 16.89 K and 22.32 K for VZA =0° and VZA =60° respectively when COD is 0.7. It indicates that the proposed three-channel algorithm can significantly improve the LST retrieval accuracy using thermal infrared data in cirrus cloudy skies. To estimate the LST errors caused by the uncertainties of COD, R, LSE and instrument noise, a sensitivity analysis was performed. It shows that the accuracy of cirrus COD is more important for the retrieval of LST compared with other parameters. The maximum total LST errors, taking into account all the input parameters' uncertainty and algorithm error itself, are 3.8 K and 4.3 K when VZA =0° and VZA =60° respectively.

Influence of emissivity angular variation on land surface temperature retrieved using the generalized split-window algorithm

The angular variation of land surface emissivity (LSE) is rarely considered in the split-window algorithm for retrieving land surface temperature (LST), and this can cause large uncertainties in LST retrievals. To analyze the influence of angular LSE variation on LST retrievals, we built a look-up table (LUT) of directional emissivities from the MYD21A1 LST/LSE product in the Moderate Resolution Imaging Spectroradiometer (MODIS) split-window channels. The extracted directional emissivities were then input into the MODIS generalized split-window (GSW) algorithm to substitute for the classification-based emissivities. A simulation analysis was first conducted based on the LUT. Furthermore, the LST retrievals estimated from MODIS observations using the directional emissivities were compared with those estimated using the classification-based emissivities. In-situ measurements from the US SURFRAD and China’s HiWATER networks were used to evaluate LST retrievals obtained using the two different emissivities. The results showed that angular LSE variations in the split-window channels for vegetated surfaces were generally minor during the daytime, but more pronounced during the night-time (approximately 0.005 between nadir and 60°). For barren surfaces, the angular LSE variation in the ˜12 μm channel was negligible but reached approximately 0.01 in the ˜11 μm channel. In the simulation, the influence of angular LSE variation was minor for view-zenith angles (VZA) <40°, but pronounced for VZA >40° reaching approximately 1.0 and 0.7 K at VZA 65° for barren and vegetated surfaces, respectively. In the evaluation, the LST estimated using the directional emissivities showed a higher accuracy than those estimated using the classification-based emissivities, especially over barren surfaces where the improvement reached >1 K. We conclude that angular LSE variation cannot be ignored in LST estimation using the GSW algorithm when VZA is >40°, especially over barren surfaces. The accuracy of the GSW algorithm is improved pronouncedly by using the directional emissivities extracted from the MYD21 product.

Land Surface Temperature Retrieval From Landsat-8 Data With the Generalized Split-Window Algorithm

This paper addresses land surface temperature (LST) retrieval from Landsat-8 data using the generalized split-window (GSW) algorithm. First, radiative transfer modeling experiment is conducted using the moderate spectral resolution atmospheric transmittance algorithm and computer model fed with SeeBor V5.0 atmospheric profile database to build a data set of LST related to brightness temperatures in Thermal Infrared Sensor (TIRS) bands 10 and 11, land surface emissivities (LSEs), and total precipitable water (TPW). Then, the GSW algorithm’s coefficients are obtained through linear regression, in which the simulated data are grouped into several sub-ranges to improve the accuracy. The impacts of noise equivalent temperature difference and uncertainty of LSEs and TPW on derived LST are evaluated. Next, the TIRS channels 10 and 11 are inter-calibrated against the channels of infrared atmospheric sounding interferometer on board Metop-B. After that, LST is retrieved from the re-calibrated and clear sky Landsat-8 data using the GSW algorithm, where LSEs are estimated from the measurements of operational land imager on Landsat-8 by the modified normalized difference vegetation index (NDVI) based emissivity method, and TPW is extracted from the european centre for medium-range weather forecasts reanalysis data. Finally, the retrieved LST is cross-validated with the MOD11_L2 V6 product. The results show that the GSW algorithm can accurately retrieve LST from Landsat-8 data, and errors mainly come from the uncertainty of LSEs and TPW. Against the MOD11_L2 V6 product, the LST errors are −1.45 ± 0.80 K and −0.49 ± 0.78 K before and after the correction of LSEs and TPW, respectively.

Estimation of Land Surface Temperature From MODIS Data for the Atmosphere With Air Temperature Inversion Profile

Air temperature inversion (ATI), which usually occurs at the near surface boundary layer in the atmosphere, influences the thermal path atmospheric upwelling and downwelling radiances. Because it is difficult to determine the occurrence of temperature inversion from satellite data, the influence of ATI on the retrieval of land surface temperature (LST) was not considered in the development of LST retrieval algorithm. This paper aims to analyze and reduce the influence of ATI on LST retrieval in the generalized split-window (GSW) algorithm. The GSW coefficients are reparameterized by using the ATI profiles from the thermodynamic initial guess retrieval cloud-free database. Comparison of the root-mean-square errors calculated before and after using the reparameterized coefficients in the GSW algorithm from the numerical simulations showed that the LST retrieval accuracy could be improved by 0.71 K for viewing zenith angle equivalent to 60°. Intercomparisons using the moderate resolution imaging spectroradiometer products MOD11_L2/MYD11_L2 and in situ LST measured at the Hailar field site in northeastern Inner Mongolia, China, showed that the LST retrieval accuracy could be improved by 0.4 K using the reparamerization coefficients in the GSW algorithm when the atmosphere is occurred by ATI.

Meteosat Land Surface Temperature Climate Data Record: Achievable Accuracy and Potential Uncertainties

The European Organization for the Exploitation of Meteorological Satellites’ (EUMETSAT) Meteosat satellites provide the unique opportunity to compile a 30+ year land surface temperature (LST) climate data record. Since the Meteosat instrument on-board Meteosat 2–7 is equipped with a single thermal channel, single-channel LST retrieval algorithms are used to ensure consistency across Meteosat satellites. The present study compares the performance of two single-channel LST retrieval algorithms: (1) A physical radiative transfer-based mono-window (PMW); and (2) a statistical mono-window model (SMW). The performance of the single-channel algorithms is assessed using a database of synthetic radiances for a wide range of atmospheric profiles and surface variables. The two single-channel algorithms are evaluated against the commonly-used generalized split-window (GSW) model. The three algorithms are verified against more than 60,000 LST ground observations with dry to very moist atmospheres (total column water vapor (TCWV) 1–56 mm). Except for very moist atmospheres (TCWV &gt; 45 mm), results show that Meteosat single-channel retrievals match those of the GSW algorithm by 0.1–0.5 K. This study also outlines that it is possible to put realistic uncertainties on Meteosat single-channel LSTs, except for very moist atmospheres: simulated theoretical uncertainties are within 0.3–1.0 K of the in situ root mean square differences for TCWV &lt; 45 mm.

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

This paper uses the refined Generalized Split-Window (GSW) algorithm to derive the land surface temperature (LST) from the data acquired by the Visible and Infrared Radiometer on FengYun 3B (FY-3B/VIRR). The coefficients in the GSW algorithm corresponding to a series of overlapping ranges for the mean emissivity, the atmospheric Water Vapor Content (WVC), and the LST are derived using a statistical regression method from the numerical values simulated with an accurate atmospheric radiative transfer model MODTRAN 4 over a wide range of atmospheric and surface conditions. The GSW algorithm is applied to retrieve LST from FY-3B/VIRR data in an arid area in northwestern China. Three emissivity databases are used to evaluate the accuracy of different emissivity databases for LST retrieval, including the ASTER Global Emissivity Database (ASTER_GED) at a 1-km spatial resolution (AG1km), an average of twelve ASTER emissivity data in the 2012 summer and emissivity spectra extracted from spectral libraries. The LSTs retrieved from the three emissivity databases are evaluated with ground-measured LST at four barren surface sites from June 2012 to December 2013 collected during the HiWATER field campaign. The results indicate that using emissivity extracted from ASTER_GED can achieve the highest accuracy with an average bias of 1.26 and −0.04 K and an average root mean square error (RMSE) of 2.69 and 1.38 K for the four sites during daytime and nighttime, respectively. This result indicates that ASTER_GED is a useful emissivity database for generating global LST products from different thermal infrared data and that using FY-3B/VIRR data can produce reliable LST products for other research areas.

Daytime Land Surface Temperature Extraction from MODIS Thermal Infrared Data under Cirrus Clouds.

Simulated data showed that cirrus clouds could lead to a maximum land surface temperature (LST) retrieval error of 11.0 K when using the generalized split-window (GSW) algorithm with a cirrus optical depth (COD) at 0.55 μm of 0.4 and in nadir view. A correction term in the COD linear function was added to the GSW algorithm to extend the GSW algorithm to cirrus cloudy conditions. The COD was acquired by a look up table of the isolated cirrus bidirectional reflectance at 0.55 μm. Additionally, the slope k of the linear function was expressed as a multiple linear model of the top of the atmospheric brightness temperatures of MODIS channels 31–34 and as the difference between split-window channel emissivities. The simulated data showed that the LST error could be reduced from 11.0 to 2.2 K. The sensitivity analysis indicated that the total errors from all the uncertainties of input parameters, extension algorithm accuracy, and GSW algorithm accuracy were less than 2.5 K in nadir view. Finally, the Great Lakes surface water temperatures measured by buoys showed that the retrieval accuracy of the GSW algorithm was improved by at least 1.5 K using the proposed extension algorithm for cirrus skies.

Extension of the Generalized Split-Window Algorithm for Land Surface Temperature Retrieval to Atmospheres With Heavy Dust Aerosol Loading

It is worth noting that the influences of dust aerosol type and different aerosol loadings were not considered in the development of the generalized split-window (GSW) algorithm. However, numerical simulations showed that the influence of dust aerosol could lead to a maximum land surface temperature (LST) retrieval error of 5.12 K when the aerosol optical depth (AOD) in the atmosphere is 1.0 and viewing zenith angle (VZA) is 60°. This paper focuses on reducing the influence of dust aerosol on the LST retrieval error of the GSW algorithm. A linear function was developed to reduce such influence with respect to the AOD. The slope could be expressed as a function of the difference between the MODIS channel brightness temperatures T31 and T32 measured at the top of the atmosphere (TOA) and difference and mean of the two-channel emissivities, and the offset could be used as a constant value for each VZA. The results showed that the retrieval accuracy could be improved by approximately 4 K for AOD = 1.0 and VZA = 60°. Sensitivity analysis in terms of the uncertainties of the input parameters showed that the maximum LST retrieval error is 1.15 K for VZA = 0°. Some of the in situ measurements observed at the Yingke site in northwest China and Arvaikheer site in south Mongolia were used to test the proposed method, respectively. The results showed that the proposed method could improve the LST retrieval accuracy by at least 1 K for the GSW algorithm in atmospheres with heavy dust aerosol loading.

A generalized split-window algorithm for land surface temperature estimation from MSG-2/SEVIRI data

This paper aims to determine land surface temperature (LST) using data from a spinning enhanced visible and infrared imager (SEVIRI) on board Meteosat Second Generation 2 (MSG-2) by using the generalized split-window (GSW) algorithm. Coefficients in the GSW algorithm are pre-determined for several overlapping sub-ranges of the LST, land surface emissivity (LSE), and atmospheric water vapour content (WVC) using the data simulated with the atmospheric radiative transfer model MODTRAN 4.0 under various surface and atmospheric conditions for 11 view zenith angles (VZAs) ranging from 0° to 67°. The results show that the root mean square error (RMSE) varies with VZA and atmospheric WVC and that the RMSEs are within 1.0 K for the sub-ranges in which the VZA is less than 30° and the atmospheric WVC is less than 4.25 g cm−2. A sensitivity analysis of LSE uncertainty, atmospheric WVC uncertainty, and instrumental noise (NEΔT) is also performed, and the results demonstrate that LSE uncertainty can result in a larger LST error than other uncertainties and that the total error for the LST is approximately 1.21 and 1.45 K for dry atmosphere and 0.86 and 2.91 K for wet atmosphere at VZA = 0° and at VZA = 67°, respectively, if the uncertainty in the LSE is 1% and that in the WVC is 20%. The GSW algorithm is then applied to the MSG-2 – SEVIRI data with the LSE determined using the temperature-independent spectral indices method and the WVC either determined using the measurements in two split-window channels or interpolated temporally and spatially using European Centre for Medium Range Weather Forecasting (ECMWF) data. Finally, the SEVIRI LST derived in this paper (SEVIRI LST1) is evaluated through comparisons with the SEVIRI LST provided by the land surface analysis satellite applications facility (LSA SAF) (SEVIRI LST2) and the Moderate Resolution Imaging Spectroradiometer (MODIS) LST product (MOD11B1 LST product). The results show that more than 80% of the differences between SEVIRI LST1 and SEVIRI LST2 are within 2 K, and approximately 70% of the differences between SEVIRI LST1 and MODIS LST are within 4 K. Furthermore, compared to MODIS LST, for four specific areas with different land surfaces, our GSW algorithm overestimates the LST by up to 1.0 K for vegetated surfaces and by 1.3 K for bare soil.