Next Article in Journal
Influence of Flight Altitude and Surface Characteristics on UAS-LiDAR Ground Height Estimate Accuracy in Juncus roemerianus Scheele-Dominated Marshes
Next Article in Special Issue
Primary Impact Evaluation of Surface Temperature Observations for Microwave Temperature Sounding Data Assimilation over Land
Previous Article in Journal
Exploring Spatial and Temporal Dynamics of Red Sea Air Quality through Multivariate Analysis, Trajectories, and Satellite Observations
Previous Article in Special Issue
Reconstruction of Land Surface Temperature Derived from FY-4A AGRI Data Based on Two-Point Machine Learning Method
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

1
Key Laboratory of Physical Electronics and Devices, Ministry of Education, Faculty of Electronic and Information Engineering, Xi’an Jiaotong University, Xi’an 710049, China
2
Land Satellite Remote Sensing Application Center, Ministry of Natural Resources, Beijing 100094, China
3
Key Laboratory of Digital Earth Sciences, Aerospace Information Research Institute, Chinese Academy of Sciences, Beijing 100094, China
4
International Research Center of Big Data for Sustainable Development Goals, Beijing 100094, China
5
University of Chinese Academy of Sciences, Beijing 100049, China
6
China Sanya Institute of South China Sea Geology, Guangzhou Marine Geological Survey, Sanya 572025, China
7
The Fourth Topographic Survey Team, ACCUR Ministry of Natural Resources, Harbin 150028, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2024, 16(2), 383; https://doi.org/10.3390/rs16020383
Submission received: 1 December 2023 / Revised: 11 January 2024 / Accepted: 15 January 2024 / Published: 18 January 2024
(This article belongs to the Special Issue Land Surface Temperature Estimation Using Remote Sensing II)

Abstract

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

1. Introduction

Land surface temperature (LST) is a crucial parameter in multiple fields of study, such as in climatology [1,2,3,4,5], ecology [6,7], and surface energy [8,9,10]. It is widely recognized as a significant Earth surface parameter [11] and is considered one of the ten essential climate variables in the land biosphere, as defined by the Global Climate Observing System [12,13].
Satellite observations offer a precise method for obtaining LST on a global scale due to their large spatial resolution and regular temporal revisiting [11]. In recent decades, various algorithms have been developed to retrieve LST from thermal infrared instruments carried on satellites, such as the single-channel (SC) algorithm [14,15,16], the split-window (SW) algorithm [17,18], multiple angle methods [17,19,20,21], and the temperature and emissivity separation (TES) algorithm [22,23,24,25,26]. A comprehensive overview of LST retrieval methods can be seen in previous research [11,13]. SW algorithm provides a method to reduce the atmospheric effect utilizing the differential absorption in two adjacent TIR channels [17,27]. The TES algorithm increases the number of equations by constructing the empirical relationship between spectral contrast and minimum emissivity and finally completes the solution of the LST inversion problem [22].
With the exception of single-channel methods, other methods require two or more thermal infrared channels. The single-channel algorithm [28] uses the radiance measured by the satellite sensor in a single channel and conducts atmospheric correction using the atmospheric radiative transfer model. LST is then retrieved from the at-sensor radiance based on the radiative transfer equation and prior knowledge of known land surface emissivity (LSE) [28,29]. Some previous research has proposed several efficient methods to estimate LST from satellite data based on the SC algorithm. Qin et al. [14] used near-surface air temperature and water vapor content instead of atmospheric profiles to estimate the LST from Landsat-5 Thematic Mapper channel 6 data. The single-channel algorithm can provide theoretically accurate LST retrieval if the LSE is known in advance. Jiménez-Muñoz et al. proposed an update and extension of the generalized single-channel algorithm which extended the application of the SC algorithm to the Landsat-4 TM sensor and enhanced the TM plus sensor onboard the Landsat-7. A series of research studies have indicated that the SC algorithm is extensively employed for retrieving LST from Landsat series data due to its simplicity and its reliance on just one thermal infrared channel [30].
The ZY1-02E satellite was successfully launched on 26 December 2021 from the Taiyuan Satellite Launch Centre, China, aboard a Long March-4 carrier rocket. Developed by the Fifth Academy of Aerospace Science and Technology Corporation, the ZY1-02E satellite is classified as a medium-resolution remote-sensing satellite. Its primary purpose is to conduct land resource surveying and monitoring, providing domestic data support for the investigation, monitoring, supervision, and production capacity monitoring of land resources. Additionally, it contributes to disaster reduction, environmental protection, housing and construction, transportation, agriculture, forestry, marine surveying, and mapping industries.
The ZY1-02E satellite carries a payload comprising three instruments—the Visible Near-Infrared Multispectral Camera (VNIC), the Advanced Hyperspectral Imager (AHSI), and the Thermal Infrared Sensor (IRS). The AHSI instrument features 166 bands, offering a resolution of 30 m and a width of 60 km. The VNIC sensors provide panchromatic imagery at 2.5 m resolution and multispectral data with eight bands at 10 m resolution covering a width of over 100 km. The IRS has a single thermal infrared band, providing a spatial resolution of 16 m and covering a width of 60 km.
Unlike the traditional thermal sensor, the ZY1-02E IRS sensor has a wider spectral coverage, i.e., 7~12 μm, which includes a stronger atmospheric water vapor absorption and ozone absorption band. The aim of this study was to evaluate the feasibility of temperature inversion using a single-channel algorithm based on WVS from single-channel ZY1-02E IRS data because of its wide spectrum range affected by both atmospheric water vapor and ozone. This paper is structured as follows: Section 2 introduces the used ZY1-02E IRS and VNIC data, atmospheric profiles, and in situ field measurements. Section 3 shows the WVS-based LST inversion method. Section 4 and Section 5 give the results and discussion. Section 6 outlines the conclusions.

2. Study Area and Data

2.1. Study Area

The global LST validation sites were summarized, and it can be found that the validation sites were selected as uniform surfaces, including water, vegetation, sand, agricultural fields, barren soil, etc., to evaluate the LST products, such as MODIS, SEVERI, VIRRS, ABI, SLSTR, ASTER, VIMS LST products. Among them, SURFRAD (Surface Radiation Budget) Network is currently the most widely used authenticity LST testing site, including grassland, rice, forest land, desert, etc. It can be seen that these verification sites are located in a large area of plains. This article is not aimed at the LST inversion over mountainous areas because the LST validation is very difficult due to the complex mountainous terrain and because the topography changes the surface structure, which leads to changes in the geometric relationship between the ground and the sensor, resulting in changes in the composition of the radiant energy received by the sensor. In addition, the shadow also varies from the sun’s location in mountainous areas. Thus, field campaigns including the commonly used surface types were conducted to validate the accuracy of the algorithms using ZY1-02E IRS data. For validation purposes, three field sites in Figure 1 were chosen: Baotou, Heihe, and the Yantai Sea. The selection of these three sites seen in Table 1 mainly considers regional differences, climate differences, and surface type differences, including desert, bare soil, water body, vegetation, wetland, which can provide a more comprehensive evaluation of temperature inversion.
The Baotou site [31,32] is located in Urad Front Banner, in Western Inner Mongolia, with coordinates of approximately 40.85°N latitude and 109.6°E longitude. The Baotou site encompasses two different target categories for validation. The first category consists of a crop target measuring 300 × 600 m, primarily comprising corn, which is ideal for validating LST. The second category comprises a desert target, measuring around 300 × 300 m. Figure 1 shows both the crop target and natural desert target at the Baotou site, with a separation distance of approximately 6 km between them. Water bodies serve as natural targets for thermal infrared calibration and validation due to their high thermal inertia and minimal spatiotemporal temperature variability. Thus, Ulansuhai Lake, located in the Urad Qianqi, was used to obtain the water temperature for validation. Additionally, bare soil and the Kubuqi desert surrounding Ulansuhai Lake were also selected for validation purposes. The Baotou site and Ulansuhai Lake have average ground elevations of 1290 m and 1021 m, respectively.
The Heihe River Basin (HRB) serves as the second site in this study. Situated in the arid region of Northwestern China, it is the second largest endorheic basin in China. The HRB is well suited for investigating land surface processes due to its diverse landscapes [33]. The hydrometeorological observatory consists of over 20 observation stations that cover the main land surfaces in the HRB. These surfaces include alpine meadows, forestlands, croplands, deserts, bare lands, and wetlands. For the purposes of validation, two stations, namely Zhangye wetland and desert, were chosen. These stations are located in the upstream region, midstream region, and downstream region in the HRB, respectively.
The third site, the Yantai Sea, is located in eastern China at a latitude of 37.51°N and a longitude of 121.47°E, which is close to the Yellow Sea and Bohai Sea, with an average elevation of 0 m. The Yantai Sea experiences a warm-temperature monsoonal continental climate, characterized by mild temperatures and abundant precipitation throughout the year. The Yantai Sea was selected as an important validation site.

2.2. ZY1-02E Data

The ZY1-02E satellite captures visible and infrared imagery to measure various parameters of the land, atmosphere, and oceans. As a cutting-edge, operational high-resolution imaging instrument, it belongs to a new generation of remote sensing technology. The ZY1-02E satellite enables the generation of various critical environmental products, such as snow and ice cover, clouds and city monitoring, sea surface temperature and LST, and vegetation and surface albedo. Table 2 provides a comprehensive summary of its detailed spectral characteristics in its 10 spectral bands, which range from 0.45 to 12.0 μm. The data were obtained from both the VNIC and IRS sensors. The VNIC sensor has nine bands with a spatial resolution of 2.5 m (Pan) and 10 m (multispectral) at nadir. Additionally, the IRS sensor includes one thermal infrared band with a 16 m spatial resolution at nadir, specifically designed for LST retrieval. The ZY1-02E VNIC and IRS spectral response functions are shown in Figure 2 and Figure 3.
Figure 3 shows the transmittance of atmospheric water vapor and ozone from 7 μm to 12 μm using MODTRAN (moderate resolution atmospheric transmission) US 1976 Standard atmosphere. The minimum of ozone transmittance ranging from 9 to 10 μm is nearly 0.2. Meanwhile, the IRS channel covers both continuous water vapor absorption and line absorption, especially 8~9 μm. Compared to the traditional Thermal Infrared Sensor used for LST estimation, the ZY1-02E IRS channel is strongly influenced by water vapor and ozone.

2.3. MODIS Data

The Moderate Resolution Imaging Spectroradiometer (MODIS) onboard the NASA Aqua and Terra Earth Observing System satellites is one of the key instruments mainly used to observe biological and physical processes around the world. MODIS provides per-pixel LST (https://modis.gsfc.nasa.gov/data/dataprod/mod11.php, accessed on 7 July 2022 to 29 March 2023) and sea surface temperature (SST, https://oceancolor.gsfc.nasa.gov/cgi/browse.pl?sub=level3&prm=SST, accessed on 7 July 2022 to 29 March 2023) products, covering from 5 min temporal to multi-day averages. In this study, two temperature products (i.e., MOD11_L2 and SST data) were used to validate the accuracy of the WVS-based method. The spatial resolution of both temperature products is 1 km. The MOD11_L2 LST product (https://doi.org/10.5067/MODIS/MOD11_L2.061, accessed on 7 July 2022 to 29 March 2023) was retrieved by using the generalized split-window algorithm [34]. A comparison was performed between the MOD11_L2 LST and ground measurements at several sites: Brookings, Audubon, Canaan Valley, and Black Hills [35]. Results show that the absolute biases and RMSEs between MOD11_L2 LST and ground measurements were less than 0.8 K and 1.7 K, respectively.

2.4. The ASTER Spectral Library

The Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) spectral library v2.0 [36], developed jointly by the Jet Propulsion Laboratory (JPL), Johns Hopkins University (JHU), and the United States Geological Survey (USGS), comprises a comprehensive collection of data. It encompasses over 2300 spectra covering a broad range of materials, such as vegetation, man-made materials, soil, minerals, rocks, water, ice, and snow. The wavelength range covered is from 0.4 µm to 15.4 µm, with specific wavelengths falling between 3.0 μm and 14.5 μm. From this library, 108 emissivity spectra were selected to represent various surface features, including 70 soil/mineral types, 28 vegetation types, and 10 man-made material types.

2.5. ERA5 Atmospheric Profiles

In order to determine the LST, atmospheric correction is essential. In this study, the European Centre for Medium-Range Weather Forecasts (ECMWF) ERA5 atmospheric profile dataset (https://cds.climate.copernicus.eu/, accessed on 7 July 2022 to 29 March 2023) was used to obtain the atmospheric parameters by using the radiative transfer code (MODTRAN 5.3) developed by AFRL/VSBT (The Air Force Research Laboratory (AFRL), established on 31 October 1997, is a dedicated scientific research organization of the United States Air Force) in collaboration with Spectral Sciences, Inc. (Spectral Sciences Incorporated (SSI), established on 1981, is an internationally recognized resource hub for expertise in spectroscopy, remote sensing and imaging, in Burlington, United States). ERA5, the fifth-generation ECMWF reanalysis, offers global climate and weather data spanning the past eight decades. The dataset begins from 1940, replacing the ERA-Interim reanalysis. By leveraging the ERA5 profiles, atmospheric parameters, such as atmospheric transmittance, atmospheric downwelling, and upwelling radiances, are derived.
ERA5 profiles provide geopotential height, temperature, relative humidity information at 37 pressure levels, and total column ozone on a 0.25° × 0.25° grid, with updates every hour in UTC time. ERA5 atmospheric profiles provide the pressure, geopotential height, air temperature, relative humidity, etc., including 0.25° of the spatial resolution, 37 pressure levels (1, 2, 3, 5, 7, 10, 20, 30, 50, 70, 100, 125, 150, 175, 200, 225, 250, 300, 350, 400, 450, 500, 550, 600, 650, 700, 750, 775, 800, 825, 850, 875, 900, 925, 950, 975, and 1000), 48 km height, and 1 h temporal resolution. It is worth noting that, due to the wide spectrum of ZY1-02E covering the ozone absorption band, the total column ozone is also applied in the algorithm process and simulation for WVS coefficients. The ERA5 profiles were used to calculate the atmospheric transmittance and upwelling and downwelling radiances with the aid of MODTRAN 5.3.

2.6. In Situ Data

At the Baotou site, SI-111 thermometers with a spectral range of 8~12 µm were employed to measure the temperature of various surfaces, including sand, crops, Ulansuhai Lake, bare soil, Kubuqi desert, etc. At the Baotou site, two thermometers were distributed in a sand area, bare soil, Kubuqi desert, and four thermometers were distributed in crop areas to capture the in situ field LSTs, and the temporal sampling interval was 2 s with an altitude of approximately 2.0 m. An additional SI-111 thermometer was used to observe the sky at an angle of 53° relative to the zenith and measure the atmospheric downwelling radiance to correct for the atmospheric effect. For Ulansuhai Lake, a 102 F Fourier transform infrared spectroradiometer (FTIR) with a spectrum ranging from 2 to 16 μm, a spectral resolution of 4 cm−1, and a field of view (FOV) of 4.8° and two thermometers were deployed on a boat to measure the water temperature. At the Heihe River Basin, two SI-111 thermometers were mounted at 6m, facing due south, with the probe facing straight down to measure the LST. Land emissivity measurements were performed using the 102F FTIR.
Two automatic water observation buoys in the Yantai Sea were equipped with an automatic weather station to measure wind speed, wind direction, atmospheric temperature, humidity, pressure, and rainfall. An SI-111 thermometer was used to measure the sea temperature.
The next equation was used to calculate the in situ LSTs from surface-leaving radiance and atmospheric downwelling radiance:
T s = B 1 ( L ( 1 ε ) L a t m ε )
where Ts represents LST, B denotes the Planck function, L is the radiance measured by the SI-111 radiometer, ε is the channel-effective land surface emissivity (LSE) specifically for the SI-111 thermometer, and Latm↓ is the atmospheric downwelling radiance calculated by the spectral response function of the SI-111 thermometer.
In this study, several natural scenes (refer to Figure 4) were selected to evaluate the accuracy of LST estimated from the ZY1-02E IRS data. Table 3 provides a summary of the main technical characteristics of the thermal instruments. Figure 4 illustrates the primary attributes of the measurement targets.

3. Methodology

3.1. WVS-Based LST Method

The theoretical basis of radiative transfer in the TIR spectral region (8–14 m) has been developed [11]. In the case of a cloud-free atmosphere at the local thermodynamic equilibrium, the radiative transfer equation (RTE) can approximately be expressed as:
L λ ( T ) = τ λ ε λ B λ ( T s ) + ( 1 - ε λ ) L a t m _ λ + L a t m _ λ
where T is the at-sensor brightness temperature (BT) and Lλ(T) is the at-sensor radiance at wavelength λ. τλ is the atmospheric transmittance at wavelength λ. ελ is the land surface emissivity. Bλ(Ts) is the Planck function of the wavelength λ. Ts is the LST. L a t m _ λ is the atmospheric upwelling radiance, and L a t m _ λ is the downwelling radiance. Obviously, the atmospheric parameters were estimated by using a radiative transfer code, such as MODTRAN, with atmospheric profiles and elevation data. Thus, the LST based on the RTE method is derived from the at-sensor radiance.
As is well known, the existing single-channel method for LST retrieval, the reanalysis of atmospheric profile datasets, such as the European Centre for Medium-Range Weather Forecasts (ECMWF) and National Centers for Environmental Prediction (NCEP), is generally preferred to correct the atmospheric effect [37,38]. However, the reanalysis of atmospheric profile data has a certain uncertainty, which may influence the accuracy of atmospheric transmittance and atmospheric upward/downward radiance, and then, the retrieval accuracy of LST will be affected. The WVS-based method proposed by Tonooka [39,40] is based on the standard atmospheric correction algorithm for ASTER/TIR. A significant enhancement of this method involves determining a water vapor scaling factor γ, which serves to mitigate the errors resulting from multiple sources. These errors include not only a bias error in the water vapor profile but also random errors in the water vapor profile, errors in the air temperature profile, and elevation errors since all of these errors contribute to the bias error in the water vapor profile. Once an appropriate γ is determined, the improved atmospheric parameters are calculated as:
τ ( θ , γ ) = τ ( θ , γ 1 ) ( γ β γ 2 β ) / ( γ 1 β γ 2 β ) τ ( θ , γ 2 ) ( γ 1 β γ β ) / ( γ 1 β γ 2 β )
L a t m _ λ ( θ , γ ) = L a t m _ λ ( θ , γ 1 ) 1 τ ( θ , γ ) 1 τ ( θ , γ 1 )
L a t m _ λ ( γ ) = a + b L a t m _ λ ( 0 , γ ) + c L a t m _ λ ( 0 , γ ) 2
γ = U U = ω ( z ) ω ( z )
where γ is the water vapor scaling factor. U and U are the total water vapor content of atmospheric parameters for water vapor profiles P and P’, respectively, where P’ is a water vapor profile scaled from a water vapor profile P by a factor of γ . ω z and ω z are water vapor amounts at an arbitrary height Z. γ 1 and γ 2 are different appropriate values (e.g., γ 1  = 1 and γ 2  = 0.7 in this paper).
γ = ln τ ( γ 2 ) γ 1 β τ ( γ 1 ) γ 2 β B ( T g ) L ( γ 1 ) 1 τ ( γ 1 ) B ( T sensor ) L ( γ 1 ) 1 τ ( γ 1 ) ( γ 1 β γ 2 β ) / ln τ ( γ 2 ) τ ( γ 1 ) 1 / β
where B ( T g ) is at-surface radiance, i.e., B ( T g ) = ε B ( T s ) + ( 1 ε ) L .
From Figure 3, it can be found that the IRS has been influenced by atmospheric water vapor and ozone. Therefore, in this study, the WVS coefficients were simulated by using the global atmospheric profile library, i.e., the Thermodynamic Initial Guess Retrieval (TIGR) database [30], including the band model parameter β and regression coefficients a , b , and c of atmospheric downward radiation. The TIGR atmospheric profile database contains 2311 representative atmospheric situations. Each profile records the atmospheric pressure, temperature, humidity, and ozone content. The atmospheric transmittance with different water vapor scaling factors γ (0.7, 0.9, and 1.0, respectively) was obtained using MODTRAN 5.3 with the band model parameter, and the β parameters were obtained using least square regressions. β is equal to 1.4072 calculated from a large amount of simulation data. The coefficients a , b , and c can be obtained from the simulated atmospheric downwelling and upwelling data, and their values are −0.3630, 2.2013, and −0.1080, respectively. Obviously, the calculation of atmospheric parameters takes into account the influence of water vapor and ozone, which provides more accurate atmospheric data for LST inversion.

3.2. Land Surface Emissivity Inversion

The estimation of land surface emissivity (LSE) relies on the fraction vegetation coverage method (FVC), which integrates two constant emissivity values representing the bare ground and full vegetation conditions of each pixel. The real-time emissivity is then adjusted according to the FVC. The classification of land surfaces into distinct types necessitates the adoption of varied processing approaches owing to their diverse thermal emission characteristics. The inland water and full vegetation emissivities are directly converted from the mean value determined according to the ASTER spectral library dataset, and the bare soil emissivity is calculated based on the relationship between emissivity and the reflectance of the red band (Tang et al., 2015 [41]), i.e., ε = a + b ρ r e d , where a and b are the channel-dependent regression coefficients and ρ r e d is the reflectance of the red band. For mixed pixels, the emissivity can be expressed as a linear relationship between the bare soil and vegetation fraction [42]:
ε = ε v f v + ε s ( 1 f v ) + d ε
where ε v and ε s are the emissivities of the vegetation and bare soil, respectively, and f v is the fractional vegetation coverage (FVC) [43]:
f v = N D V I N D V I min N D V I max N D V I min
where normalized difference vegetation index (NDVI) is acquired from the target pixel and N D V I max and N D V I min are calculated from the full vegetation and bare soil pixels in the whole image.
d ε = 4 ( 1 ε s ) ε v F f v ( 1 f v ) 2
where F is the shape factor depending on the vegetation structure information in which F is set up to be 0.55 [44].
In total, the emissivity can be written as Equation (11), where pixels with NDVI values below 0 are considered as water bodies (ε = εwater). The pixels with NDVI values below 0.2 are considered bare soil. NDVI values above 0.5 represent full vegetation pixels. Therefore, the pixels with NDVI values ranging from 0.2 to 0.5 are considered as a mixture of soil and vegetation [41,45].
ε = ε water , NDVI < 0 ε = a + b ρ r e d , 0 NDVI < 0.2 ε v f v + ε s ( 1 f v ) + d ε , 0.2 NDVI 0.5 ε v , NDVI > 0.5
where ε water is the emissivity of water, ε s is the emissivity of soil, and ε v is the emissivity of vegetation. ρ red is the land surface reflectance of the red band.

3.3. Atmospheric Parameter Inversion

To obtain the atmospheric parameters for each pixel of ZY1-02E IRS images, the interpolation of atmospheric profiles in time and space must be performed, specifically temporal linear interpolation of the atmospheric parameters in terms of the time of the ERA5 profiles (https://cds.climate.copernicus.eu/#!/home, accessed on 7 July 2022 to 29 March 2023) and the ZY1-02E acquisition time and spatial interpolation of the atmospheric parameters in terms of geographic latitudes and longitudes of the closest four grid points. In addition, there is a certain uncertainty in the reanalysis of atmospheric profile data, and the water vapor scaling factor γ should be calculated with simulation data using the WVS-based LST method. Then, the modified atmospheric parameters can be obtained.

3.4. WVS-Based LST Retrieval

DN values in ZY1-02E IRS band data should be converted to top-of-the-atmosphere (TOA) radiance using the calibration coefficients provided in the metadata file.
L = a D N + b
where L is TOA radiance, DN is a digital number of IRS data, and a and b are calibration coefficients.
In addition, the TOA radiance should be converted to TOA brightness temperature (BT). Considering the wide spectral range (7~12 μm) of ZY1-02E IRS, the look-up table (LUT) between radiance and temperature of ZY1-02E should be set up from 200 K to 400 K with a step of 0.1 K.
The LST retrieval process is depicted in Figure 5 as the following steps:
  • Land surface emissivity retrieval. Firstly, land surface reflectance is retrieved from the ZY1-02E VNIC data using the ERA 5 atmospheric profile from MODTRAN. Then, the NDVI and FVC can be obtained from land surface reflectance data. Subsequently, the land surface emissivity is accurately retrieved based on FVC and NDVI.
  • TOA brightness temperature retrieval: The DN values received by the ZY1-02E IRS data are converted into radiance data with the calibration parameters. Subsequently, the TOA radiance is converted into TOA BT using the look-up table between radiance and brightness temperature.
  • Adjusted atmospheric parameters: The atmospheric parameters in the ERA5 spatial scale can be calculated using ERA5 profile data from MODTRAN. Then, the atmospheric parameters in the ZY1-02E spatial scale can be estimated using temporal and spatial interpolation. Finally, the WVS method is used to adjust the ZY1-02E atmospheric parameters.
  • Based on the adjusted atmospheric parameters and LSE data, the ZY1-02E LST can be inversed using the thermal radiative transfer Equation (2).

4. Results

4.1. LST Results

Figure 6 shows the false color band composition (Band 7, 4, 3 | red, green, blue) from ZY1-02E VNIC images for the Heihe River Basin and Baotou, encompassing surfaces characterized by vegetation, rock, soil, and urban areas. Table 4 presents the number of images used for the validation of sites in each study area.
The LST was estimated from ZY1-02E IRS data using a WVS-based method. Based on the FVC and NDVI calculated using land surface reflectance and the emissivity for the Ulansuhai, Kubuqi desert, Baotou desert, Zhangye wetland, desert, and Yantai Sea sites were estimated. Emissivity measurements for bare soil were obtained on different days due to variations in local FVC and NDVI.
Some of the inversed LST images of three study areas derived from ZY1-02E IRS data are presented in Figure 7. We selected four specific days of ZY1-02E IRS images that cover the Baotou site, HRB, and Yantai Sea. The average UTC time of the ZY1-02E overpass at the Baotou site, HRB, and Yantai Sea is 03:50, 04:20, and 03:00, respectively. Figure 7 demonstrates that the temperature of Ulansuhai Lake during local noon is typically lower than other landscapes while bare soil and sand exhibit higher temperatures. The reflectance of water bodies exceeds that of other land cover types.

4.2. Validation

Two validation methods were selected, i.e., based on in situ data and based on cross-validation with the MODIS LST product, and employed to evaluate the LST inversion accuracy. A total of 25 scenes of ZY1-02E satellite images have been used for statistical analysis, covering the period from July 2022 to March 2023. The number of sampled points from the Baotou site, HRB, and Yantai Sea are 14, 6, and 5, respectively. In addition, cross-validation with MODIS LST products has also been conducted for further analysis.

4.2.1. Validation Based on In Situ Data

The LST derived from the WVS-based method was evaluated using comparisons with the ground measurements (GMs) LST obtained from the Baotou, HRB, and Yantai Sea sites. The ZY-1F satellite images, including those from July, September, October, November 2022 and February and March 2023, have some images for validation. The in situ averaged LST before and after 10 min against the ZY1-02E overpass time was selected as the validation data and compared with the inversed LST from these sites. During the validation of the LST retrieval, we employed the average LST value of 3 × 3 pixels surrounding the target site in the ZY1-02E IRS images. Figure 7 shows scatterplots depicting the inversed LST versus the in situ LST at the eight sites. The results in Figure 8a reveal that the average bias and root mean square error (RMSE) between the WVS-based method and in situ LSTs are 0.63 K and 1.62 K, respectively. This shows a high agreement between the WVS method and field measurement data. Additionally, the results in Figure 7b indicate that the LST inversion in the water body, desert, sand, reed wetland, and vegetation exhibits the closest proximity to the measured temperature data, with an average LST bias of 0.24, −0.17, −0.7, 0.27, and 0.27 K, respectively. Comparatively, the average LST bias for the bare soil is 1.7 K. And the RMSEs of these six land surface types (water body, bare soil, desert, sand, reed wetland, and vegetation) are 1.39 K, 1.77 K, 0.53 K, 0.99 K, 0.42 K, and 0.42 K, respectively.
In addition, according to each surface type, the statistical analysis shows that LST inversion in the water body, bare soil, desert, sand, reed wetland, and vegetation exhibit the closest proximity to the measured temperature data, with an average LST bias and RMSE of 0.24 K and 1.39 K, 1.7 K and 1.77 K, −0.17 K and 0.53 K, −0.7 K and 0.99 K, 0.27 K and 0.42 K, and 0.27 K and 0.42 K, respectively.

4.2.2. Cross-Validation Compared to MODIS LST and SST Products

As is well known, the MODIS LST product offers high accuracy, making the MODIS Terra daily LST product (MOD11_L2) with a spatial resolution of 1 km suitable for cross-validation. According to satellite overpass time, the time differences between the MOD11_L2 and ZY1-02E data range from 0 to 50 min for study areas A and B. To minimize the error caused by the difference due to overpass time, a temporal correction method was considered to correct the MODIS temperature product using the measured surface temperatures over the two satellites’ overpass times, i.e., LST = a × time + b. The coefficients a and b can be calculated using field measurement data. This is because the surface temperature shows a certain change rule over a short period of time. The ground measurement temperature data were plotted as a function of viewing time in Figure 9, and then, the time-corrected MODIS LST values were obtained with a linear function expression.
The results in Figure 10a indicate that the average bias and RMSE between the inversed LSTs from the WVS-based method and the MODIS LST/SST product are 0.77 K and 2.0 K, respectively. Figure 10b reveals that the LST biases between the WVS-based method and the MODIS LST product are consistent across five land surface types (water body, bare soil, desert, sand, and reed wetland), with average biases of 0.88 K, 1.5 K, −0.73 K, 1.3 K, and 3.9 K, respectively. And the RMSEs are 1.75 K, 1.51 K, 0.71 K, 1.86 K, and 4.05 K, respectively.
The data collection period for the ZY1-02E IRS scenes spanned from July 2022 to March 2023. For evaluating the accuracy of LST retrieval from ZY1-02E IRS data, Figure 11 presents the results that the absolute errors (AEs) between the WVS-based method and ground measurement range from 0.1 K to 3.5 K, with an average AE of 1.31 K. Moreover, the average absolute error between the WVS-based method and MODIS LST&SST product is 1.82 K.
To validate the reliability of the WVS-based LST method, we conducted calculations using the variance of the LST values for the target sites and their surrounding eight pixels in the retrieved LST images, with a spatial resolution of 48m for the variance statistics. The results in Figure 10 demonstrate that the average variance of LST at the target site was 0.134 K, indicating a minimal deviation in the retrieved LST for the target stations. These findings affirm that the proposed method can accurately retrieve LST from single-channel TIR data.

5. Discussion

In this study, the single-channel algorithm based on a water vapor scaling (WVS) method proposed by Tonooka [39] to correct the error caused by uncertain atmospheric profile data was used to estimate LST from ZY1-02E thermal infrared data. The ZY1-02E IRS data is more affected by the atmosphere due to its wide channel compared to Landsat 8, 9, and ASTER. Therefore, the LST retrieval method based on WVS considers the effects of the strong absorption of water vapor (7~8.2 μm) and ozone (~9.6 μm) during the model construction process.
In terms of the validation based on in situ data, the WVS-based method exhibits poorer performance on bare soil and Zhangye wetland sites, as illustrated in Figure 12. The challenge of accurately defining the emissivity of Zhangye wetland and Baotou bare soil sites results in significant LST deviation, leading to the reduced accuracy of LST inversion for these locations. To quantitatively describe the accuracy of LST inversion, the biases and RMSEs were calculated as the evaluation index [46,47]. The average LST biases between the WVS-based method and in situ LSTs for the Ulansuhai, bare soil, Kubuqi desert, Baotou sand, Zhangye wetland, desert, and Yantai Sea sites are as follows: 1.3 K, 1.4 K, −1.3 K, 1.3 K, 3.9 K, −0.69 K, and 0.65 K, respectively. Therefore, the WVS-based method demonstrates higher accuracy for uniform land cover, such as water bodies and deserts. The results show the effectiveness of LST inversed from ZY1-02E IRS data using the WVS-based method. Although the spatial resolution of Landsat 9’s Thermal Infrared Sensor 2 (TIRS-2) (Landsat 9 is the ninth satellite in the U.S. Landsat program, which carries a NASA-built thermal infrared sensor 2) and ECOTRESS LST product are higher compared to the MODIS LST product, by 100 m and 70 m, respectively, the transit dates for either Landsat or ECOTRESS are different with ZY1-02E. Cross-validation based on MODIS LST and SST products was performed. To ensure comprehensive coverage of the corresponding areas in the MODIS LST product, we analyzed the average LST value and LST variance per site in ZY1-02E IRS images using the area-weighted pixel aggregation algorithm [48]. The method used to correct the temporal difference between ZY1-02E and MODIS has been introduced in Section 4.2.2. Figure 12 illustrates that the average LST biases between the WVS-based method and MODIS are lower in Baotou sites compared to HRB sites and Yantai Sea. The vegetation area of the Baotou Crop site is less than the MODIS single pixel size of 1 × 1 km, so the retrieval LST of the Baotou Crop was not considered in the comparison with the MODIS LST product. The dissimilarity in spatial resolution accounts for another possible reason since the ZY1-02E has a resolution of 16m, while the MODIS temperature product has a resolution of 1KM. This discrepancy can result in inaccuracies when observing bare soil, mainly due to scaling effects within the area-averaged FOV [49]. The results of average LST bias between the WVS-based method and the MODIS product for each study site indicate that the homogeneity of land cover types plays a critical role in LST retrieval across various spatial scales. This is evidence that the lower LST biases can be found in water bodies, deserts, and vegetated areas compared to other mixed surface features.
This study employed two validation methods to assess the accuracy of retrieved LSTs. The results demonstrate a close approximation between the WVS-based method and the in situ LST, exhibiting an average bias of 0.63 K and RMSE of 1.62 K. However, this study’s focus was limited to six types of landscapes across three study areas. Future research should explore a wider range of landscape types to comprehensively assess the applicability of the WVS-based method. Additional and more data samples should also be collected for specific stations in future research.

6. Conclusions

In this study, a water vapor scaling (WVS-based) method was used to estimate the LST from ZY1-02E IRS data and evaluate its feasibility. The ERA5 atmospheric profile dataset and the radiative transfer code were used to correct atmospheric effects. A water vapor scaling factor was used to reduce the errors in atmospheric correction caused by various atmospheric factors. Additionally, the LSE was estimated based on NDVI and FVC data. We obtained ZY1-02E data for three study areas from July 2022 to March 2023. To ensure the accuracy and reliability of the data, we conducted a series of validation and screening procedures and ultimately determined the samples for comparison. We compared the LST retrieval from ZY1-02E IRS data with the in situ LST and MOD11_L2 products.
The estimated LST was validated by collecting in situ LST data from three study areas. For study area A (Baotou site), the WVS-based method showed an absolute bias of approximately 0.2 to 3.3 K compared to the in situ LST, with the RMSEs ranging from 0.42 to 2.3 K. In study area B (Heihe River Basin), the LST biases ranged from 0.1 to 2.9 K, with the RMSEs between 0.5 and 2.4 K. In study area C (Yantai Sea), the SST biases ranged from 0.1 to 1.6 K, with an RMSE of ~0.8 K. In addition, according to each surface type, the statistical analysis shows that LST inversion in the water body, bare soil, desert, sand, reed wetland, and vegetation exhibit the closest proximity to the measured temperature data, with an average LST bias and RMSE of 0.24 K and 1.39 K, 1.7 K and 1.77 K, −0.17 K and 0.53 K, −0.7 K and 0.99 K, 0.27 K and 0.42 K, and 0.27 K and 0.42 K, respectively. These results show that the WVS-based method can accurately retrieve LST from single-channel thermal infrared data, though there is a wide spectrum range, i.e., 7~12 μm, affected strongly by both atmospheric water vapor and ozone. However, further validation work is necessary to evaluate the method’s performance across different land cover types and geographical locations.

Author Contributions

Methodology, X.D., K.L., C.M. and Y.Q.; Validation, J.L. and F.W.; Investigation, J.Y.; Writing—review & editing, Q.Z., Y.L. and X.L.; Supervision, H.T. and J.C.; Project administration, T.W. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported in part by the Key Program of National Natural Science Foundation of China under Grant 42130108 and in part by the National Natural Science Foundation of China under Grant 42022045, International Cooperation in Science and Technology Innovation among Governments under Grant 2019YFE0127200, and Shan’xi Key Research and Development Program under Grant 2022ZDLSF06-09.

Data Availability Statement

Data are contained within the article.

Acknowledgments

The authors sincerely acknowledge the support of the China Centre for Resources Satellite Data and Application for providing the Yantai Sea data.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bechtel, B. A New Global Climatology of Annual Land Surface Temperature. Remote Sens. 2015, 7, 2850–2870. [Google Scholar] [CrossRef]
  2. Bright, R.M.; Davin, E.; O’Halloran, T.; Pongratz, J.; Zhao, K.; Cescatti, A. Local temperature response to land cover and management change driven by non-radiative processes. Nat. Clim. Change 2017, 7, 296–302. [Google Scholar] [CrossRef]
  3. Cai, M.; Ren, C.; Xu, Y.; Lau, K.K.-L.; Wang, R. Investigating the relationship between local climate zone and land surface temperature using an improved WUDAPT methodology—A case study of Yangtze River Delta, China. Urban Clim. 2018, 24, 485–502. [Google Scholar] [CrossRef]
  4. Keenan, T.F.; Riley, W.J. Greening of the land surface in the world’s cold regions consistent with recent warming. Nat. Clim. Change 2018, 8, 825–828. [Google Scholar] [CrossRef] [PubMed]
  5. Shen, Y.; Shen, H.; Cheng, Q.; Zhang, L. Generating Comparable and Fine-Scale Time Series of Summer Land Surface Temperature for Thermal Environment Monitoring. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2021, 14, 2136–2147. [Google Scholar] [CrossRef]
  6. Myneni, R.B.; Keeling, C.D.; Tucker, C.J.; Asrar, G.; Nemani, R.R. Increased plant growth in the northern high latitudes from 1981 to 1991. Nature 1997, 386, 698–702. [Google Scholar] [CrossRef]
  7. Peng, S.-S.; Piao, S.; Zeng, Z.; Ciais, P.; Zhou, L.; Li, L.Z.X.; Myneni, R.B.; Yin, Y.; Zeng, H. Afforestation in China cools local land surface temperature. Proc. Natl. Acad. Sci. USA 2014, 111, 2915–2919. [Google Scholar] [CrossRef]
  8. Tajfar, E.; Bateni, S.M.; Lakshmi, V.; Ek, M. Estimation of surface heat fluxes via variational assimilation of land surface temperature, air temperature and specific humidity into a coupled land surface-atmospheric boundary layer model. J. Hydrol. 2020, 583, 124577. [Google Scholar] [CrossRef]
  9. Xu, T.; Guo, Z.; Xia, Y.; Ferreira, V.G.; Liu, S.; Wang, K.; Yao, Y.; Zhang, X.; Zhao, C. Evaluation of twelve evapotranspiration products from machine learning, remote sensing and land surface models over conterminous United States. J. Hydrol. 2019, 578, 124105. [Google Scholar] [CrossRef]
  10. Li, K.; Qian, Y.; Wang, N.; Li, W.; Qiu, S.; Ma, L.; Li, C.; Sun, D.; Liu, Y.; Ni, L. A Four-Component Parameterized Directional Thermal Radiance Model for Row Canopies. IEEE Trans. Geosci. Remote Sens. 2022, 60, 1–15. [Google Scholar] [CrossRef]
  11. Li, Z.-L.; Tang, B.-H.; Wu, H.; Ren, H.; Yan, G.; Wan, Z.; Trigo, I.F.; Sobrino, J.A. Satellite-derived land surface temperature: Current status and perspectives. Remote Sens. Environ. 2013, 131, 14–37. [Google Scholar] [CrossRef]
  12. Hollmann, R.; Merchant, C.J.; Saunders, R.; Downy, C.; Buchwitz, M.; Cazenave, A.; Chuvieco, E.; Defourny, P.; de Leeuw, G.; Forsberg, R.; et al. The ESA Climate Change Initiative: Satellite Data Records for Essential Climate Variables. Bull. Am. Meteorol. Soc. 2013, 94, 1541–1552. [Google Scholar] [CrossRef]
  13. Li, Z.-L.; Wu, H.; Duan, S.-B.; Zhao, W.; Ren, H.; Liu, X.; Leng, P.; Tang, R.; Ye, X.; Zhu, J.; et al. Satellite Remote Sensing of Global Land Surface Temperature: Definition, Methods, Products, and Applications. Rev. Geophys. 2023, 61, e2022RG000777. [Google Scholar] [CrossRef]
  14. Qin, Z.; Karnieli, A.; Berliner, P. A mono-window algorithm for retrieving land surface temperature from Landsat TM data and its application to the Israel-Egypt border region. Int. J. Remote Sens. 2001, 22, 3719–3746. [Google Scholar] [CrossRef]
  15. Jiménez-Muñoz, J.C.; Sobrino, J.A. A generalized single-channel method for retrieving land surface temperature from remote sensing data. J. Geophys. Res. Atmos. 2003, 108. [Google Scholar] [CrossRef]
  16. Jimenez-Munoz, J.C.; Cristobal, J.; Sobrino, J.A.; Soria, G.; Ninyerola, M.; Pons, X. Revision of the Single-Channel Algorithm for Land Surface Temperature Retrieval From Landsat Thermal-Infrared Data. IEEE Trans. Geosci. Remote Sens. 2009, 47, 339–349. [Google Scholar] [CrossRef]
  17. McMillin, L.M. Estimation of Sea Surface Temperatures from Two Infrared Window Measurements with Different Absorption. J. Geophys. Research. Part C Ocean. 1975, 80, 5113–5117. [Google Scholar] [CrossRef]
  18. Atitar, M.; Sobrino, J.A. A Split-Window Algorithm for Estimating LST From Meteosat 9 Data: Test and Comparison with In Situ Data and MODIS LSTs. IEEE Geosci. Remote Sens. Lett. 2009, 6, 122–126. [Google Scholar] [CrossRef]
  19. Becker, F. The impact of spectral emissivity on the measurement of land surface temperature from a satellite. Int. J. Remote Sens. 1987, 8, 1509–1522. [Google Scholar] [CrossRef]
  20. Barton, I.J.; Zavody, A.M.; O’Brien, D.M.; Cutten, D.R.; Saunders, R.W.; Llewellyn-Jones, D.T. Theoretical algorithms for satellite-derived sea surface temperatures. J. Geophys. Res. 1989, 94, 3365–3375. [Google Scholar] [CrossRef]
  21. Caselles, V.; Coll, C.; Valor, E. Land surface emissivity and temperature determination in the whole HAPEX-Sahel area from AVHRR data. Int. J. Remote Sens. 1997, 18, 1009–1027. [Google Scholar] [CrossRef]
  22. Gillespie, A.; Rokugawa, S.; Matsunaga, T.; Cothern, J.S.; Hook, S.; Kahle, A.B. A temperature and emissivity separation algorithm for Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) images. IEEE Trans. Geosci. Remote Sens. 1998, 36, 1113–1126. [Google Scholar] [CrossRef]
  23. Gillespie, A.R.; Abbott, E.A.; Gilson, L.; Hulley, G.; Jiménez-Muñoz, J.-C.; Sobrino, J.A. Residual errors in ASTER temperature and emissivity standard products AST08 and AST05. Remote Sens. Environ. 2011, 115, 3681–3694. [Google Scholar] [CrossRef]
  24. Jiménez-Muñoz, J.C.; Sobrino, J.A.; Gillespie, A.; Sabol, D.; Gustafson, W.T. Improved land surface emissivities over agricultural areas using ASTER NDVI. Remote Sens. Environ. 2006, 103, 474–487. [Google Scholar] [CrossRef]
  25. Malakar, N.K.; Hulley, G.C. A water vapor scaling model for improved land surface temperature and emissivity separation of MODIS thermal infrared data. Remote Sens. Environ. 2016, 182, 252–264. [Google Scholar] [CrossRef]
  26. Ma, C.; Qian, Y.; Li, K.; Dou, X.; Shen, H.; Tang, H.; Qiu, S.; Zhang, L.; Jia, Y.; Ou-Yang, G. Temperature and Emissivity Retrieval From Hyperspectral Thermal Infrared Data Using Dictionary-Based Sparse Representation for Emissivity. IEEE Trans. Geosci. Remote Sens. 2023, 61, 1–16. [Google Scholar] [CrossRef]
  27. Becker, F.; Li, Z.-L. Towards a local split window method over land surfaces. Int. J. Remote Sens. 1990, 11, 369–393. [Google Scholar] [CrossRef]
  28. Hook, S.J.; Gabell, A.R.; Green, A.A.; Kealy, P.S. A comparison of techniques for extracting emissivity information from thermal infrared data for geologic studies. Remote Sens. Environ. 1992, 42, 123–135. [Google Scholar] [CrossRef]
  29. Mushkin, A.; Balick, L.K.; Gillespie, A.R. Extending surface temperature and emissivity retrieval to the mid-infrared (3–5 μm) using the Multispectral Thermal Imager (MTI). Remote Sens. Environ. 2005, 98, 141–151. [Google Scholar] [CrossRef]
  30. Wang, M.; Zhang, Z.; Hu, T.; Liu, X. A Practical Single-Channel Algorithm for Land Surface Temperature Retrieval: Application to Landsat Series Data. J. Geophys. Res. Atmos. 2019, 124, 299–316. [Google Scholar] [CrossRef]
  31. Ma, L.; Zhao, Y.; Woolliams, E.R.; Dai, C.; Wang, N.; Liu, Y.; Li, L.; Wang, X.; Gao, C.; Li, C.; et al. Uncertainty Analysis for RadCalNet Instrumented Test Sites Using the Baotou Sites BTCN and BSCN as Examples. Remote Sens. 2020, 12, 1696. [Google Scholar] [CrossRef]
  32. Liu, Y.; Ma, L.; Wang, N.; Qian, Y. Vicarious radiometric calibration/validation of Landsat-8 operational land imager using a ground reflected radiance-based approach with Baotou site in China. J. Appl. Remote Sens. 2017, 11, 044004. [Google Scholar] [CrossRef]
  33. Liu, S.; Li, X.; Xu, Z.; Che, T.; Xiao, Q.; Ma, M.; Liu, Q.; Jin, R.; Guo, J.; Wang, L.; et al. The Heihe Integrated Observatory Network: A Basin-Scale Land Surface Processes Observatory in China. Vadose Zone J. VZJ 2018, 17, 1–21. [Google Scholar] [CrossRef]
  34. Wan, Z.; Hook, S.; Hulley, G. MODIS/Terra Land Surface Temperature/Emissivity 5-Min L2 Swath 1 km V061. 2021. Available online: https://lpdaac.usgs.gov/products/mod11_l2v061/ (accessed on 8 July 2020).
  35. Wang, W.; Liang, S.; Meyers, T. Validating MODIS land surface temperature products using long-term nighttime ground measurements. Remote Sens. Environ. 2008, 112, 623–635. [Google Scholar] [CrossRef]
  36. Baldridge, A.M.; Hook, S.J.; Grove, C.I.; Rivera, G. The ASTER spectral library version 2.0. Remote Sens. Environ. 2009, 113, 711–715. [Google Scholar] [CrossRef]
  37. Yang, J.; Duan, S.-B.; Zhang, X.; Wu, P.; Huang, C.; Leng, P.; Gao, M. Evaluation of Seven Atmospheric Profiles from Reanalysis and Satellite-Derived Products: Implication for Single-Channel Land Surface Temperature Retrieval. Remote Sens. 2020, 12, 791. [Google Scholar] [CrossRef]
  38. Skoković, D.; Sobrino, J.A.; Jiménez-Muñoz, J.C. Vicarious Calibration of the Landsat 7 Thermal Infrared Band and LST Algorithm Validation of the ETM+ Instrument Using Three Global Atmospheric Profiles. IEEE Trans. Geosci. Remote Sens. 2017, 55, 1804–1811. [Google Scholar] [CrossRef]
  39. Tonooka, H. An atmospheric correction algorithm for thermal infrared multispectral data over land-a water-vapor scaling method. IEEE Trans. Geosci. Remote Sens. 2001, 39, 682–692. [Google Scholar] [CrossRef]
  40. Tonooka, H. Accurate atmospheric correction of ASTER thermal infrared imagery using the WVS method. IEEE Trans. Geosci. Remote Sens. 2005, 43, 2778–2792. [Google Scholar] [CrossRef]
  41. Tang, B.H.; Shao, K.; Li, Z.L.; Wu, H.; Tang, R. An improved NDVI-based threshold method for estimating land surface emissivity using MODIS satellite data. Int. J. Remote Sens. 2015, 36, 4864–4878. [Google Scholar] [CrossRef]
  42. Sobrino, J.A.; Jimenez-Munoz, J.C.; Paolini, L. Land surface temperature retrieval from LANDSAT TM 5. Remote Sens. Environ. 2004, 90, 434–440. [Google Scholar] [CrossRef]
  43. Shobairi, O.; Li, M. Dynamic Modelling of VFC from 2000 to 2010 Using NDVI and DMSP/OLS Time Series: A Study in Guangdong Province, China. J. Geogr. Inf. Syst. 2016, 8, 205–223. [Google Scholar] [CrossRef]
  44. Zhao, E.; Gao, C.; Han, Q.; Yao, Y.; Wang, Y.; Yu, C.; Yu, H. An Operational Land Surface Temperature Retrieval Methodology for Chinese Second-Generation Huanjing Disaster Monitoring Satellite Data. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2022, 15, 1283–1292. [Google Scholar] [CrossRef]
  45. Sobrino, J.A.; Raissouni, N.; Li, Z.-L. A Comparative Study of Land Surface Emissivity Retrieval from NOAA Data. Remote Sens. Environ. 2001, 75, 256–266. [Google Scholar] [CrossRef]
  46. Kumar, D.; Shekhar, S. Statistical analysis of land surface temperature–vegetation indexes relationship through thermal remote sensing. Ecotoxicol. Environ. Saf. 2015, 121, 39–44. [Google Scholar] [CrossRef] [PubMed]
  47. Kisi, O.; Kerem Cigizoglu, H. Comparison of different ANN techniques in river flow prediction. Civ. Eng. Environ. Syst. 2007, 24, 211–231. [Google Scholar] [CrossRef]
  48. Qian, Y.-G.; Li, Z.-L.; Nerry, F. Evaluation of land surface temperature and emissivities retrieved from MSG/SEVIRI data with MODIS land surface temperature and emissivity products. Int. J. Remote Sens. 2013, 34, 3140–3152. [Google Scholar] [CrossRef]
  49. Obata, K.; Miura, T.; Yoshioka, H. Analysis of the Scaling Effects in the Area-Averaged Fraction of Vegetation Cover Retrieved Using an NDVI-Isoline-Based Linear Mixture Model. Remote Sens. 2012, 4, 2156–2180. [Google Scholar] [CrossRef]
Figure 1. The study areas used for validation.
Figure 1. The study areas used for validation.
Remotesensing 16 00383 g001
Figure 2. Relative spectral response function (a) ZY1-02E VNIC; (b) ZY1-02E IRS.
Figure 2. Relative spectral response function (a) ZY1-02E VNIC; (b) ZY1-02E IRS.
Remotesensing 16 00383 g002
Figure 3. Total transmittance of H2O and ozone from 7 μm to 12 μm.
Figure 3. Total transmittance of H2O and ozone from 7 μm to 12 μm.
Remotesensing 16 00383 g003
Figure 4. Photos of targets in study area.
Figure 4. Photos of targets in study area.
Remotesensing 16 00383 g004
Figure 5. Flow chart of LST retrieval from ZY1-02E IRS data.
Figure 5. Flow chart of LST retrieval from ZY1-02E IRS data.
Remotesensing 16 00383 g005
Figure 6. False color images of sites. (a) Ulansuhai Lake on 13 July 2022; (b) Baotou site on 10 July 2022; (c) Yantai Sea on 7 September 2022; (d) HRB site on 22 November 2022.
Figure 6. False color images of sites. (a) Ulansuhai Lake on 13 July 2022; (b) Baotou site on 10 July 2022; (c) Yantai Sea on 7 September 2022; (d) HRB site on 22 November 2022.
Remotesensing 16 00383 g006
Figure 7. The retrieved LST images of Ulansuhai Lake (a), Baotou site (b), Yantai Sea (c), and HRB (d) from ZY1-02E IRS data based on WVS method.
Figure 7. The retrieved LST images of Ulansuhai Lake (a), Baotou site (b), Yantai Sea (c), and HRB (d) from ZY1-02E IRS data based on WVS method.
Remotesensing 16 00383 g007
Figure 8. (a) Scatterplots of the inversed LST using WVS-based method versus the in situ LST at six land surface types (water body, bare soil, desert, sand, reed wetland, and vegetation). (b) LST bias between WVS-based method and ground measurements. Different colors represent different surface types.
Figure 8. (a) Scatterplots of the inversed LST using WVS-based method versus the in situ LST at six land surface types (water body, bare soil, desert, sand, reed wetland, and vegetation). (b) LST bias between WVS-based method and ground measurements. Different colors represent different surface types.
Remotesensing 16 00383 g008
Figure 9. (a) LST change in desert site on 22 October 2022. (b) LST change in Kubuqi desert site on 7 July 2022.
Figure 9. (a) LST change in desert site on 22 October 2022. (b) LST change in Kubuqi desert site on 7 July 2022.
Remotesensing 16 00383 g009
Figure 10. (a) Scatterplots of the inversed Ts (LST) versus the MODIS LST product at five land surface types (water body, bare soil, desert, sand, reed wetland). (b) LST bias between WVS-based method and MOD11_L2.
Figure 10. (a) Scatterplots of the inversed Ts (LST) versus the MODIS LST product at five land surface types (water body, bare soil, desert, sand, reed wetland). (b) LST bias between WVS-based method and MOD11_L2.
Remotesensing 16 00383 g010
Figure 11. Evaluating accuracy of LST retrieval using WVS-based method. The boxplots are centered on the errors of retrieved LST using WVS-based method and other retrieval LST methods.
Figure 11. Evaluating accuracy of LST retrieval using WVS-based method. The boxplots are centered on the errors of retrieved LST using WVS-based method and other retrieval LST methods.
Remotesensing 16 00383 g011
Figure 12. The LST bias retrieved by WVS-based method and MODIS in all study areas.
Figure 12. The LST bias retrieved by WVS-based method and MODIS in all study areas.
Remotesensing 16 00383 g012
Table 1. Information of the seven stations selected in the study area for validation.
Table 1. Information of the seven stations selected in the study area for validation.
No.StationLand CoverStudy AreaLongitude (°)Latitude (°)Elevation (m)LSE
1UlansuhaiWater bodyBaotou108.7706E40.8476N9770.9869
2Bare soilBare soilBaotou108.8176E40.7978N9770.9257
3Kubuqi desertDesertBaotou108.6203E40.4551N9770.8836
4Baotou sandSandBaotou109.6187E40.8659N12960.9159
5Baotou CropVegetationBaotou109.5537E40.8708N12950.9718
6Zhangye wetlandReed wetlandHRB100.4464E38.9751N14600.9869
7DesertReaumuria desertHRB100.9872E42.1135N10540.8836
8Yantai SeaWater bodyYantai121.4653E37.5148N00.9869
Table 2. Spectral characteristics of ZY1-02E VNIC/IRS.
Table 2. Spectral characteristics of ZY1-02E VNIC/IRS.
BandsBands No.Spectral Range (µm)Resolution (m)NEDT/SNR
VNICPan0.45~0.902.5≥28 dB@sun altitude angle is 30° and surface reflectance is 0.03
≥48 dB@sun altitude angle is 70° and surface reflectance is 0.5
B10.45~0.5210
B20.52~0.5910
B30.63~0.6910
B40.77~0.8910
B50.40~0.4510
B60.59~0.62510
B70.705~0.74510
B80.860~1.04010
IRSB97~1216NEΔT ≤ 0.1 K@300K
Table 3. Main technical specifications for the thermal instruments.
Table 3. Main technical specifications for the thermal instruments.
InstrumentSpectral Range (µm)Operating Environment (°C)Accuracy ResolutionFOV (°)
SI-1118~14−55~80±0.2 K0.1 K44
KT-159.6~11.50~55±0.5 K0.06 K2
102 F2~1615~351 cm−14 cm−14.8
Table 4. The number of images used for validation of sites in each study area.
Table 4. The number of images used for validation of sites in each study area.
Study AreaSiteNumber of Images
BaotouUlansuhai3
BaotouBare soil3
BaotouKubuqi desert3
BaotouBaotou sand2
BaotouBaotou Crop3
HRBZhangye wetland4
HRBDesert3
YantaiYantai Sea5
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Dou, X.; Li, K.; Zhang, Q.; Ma, C.; Tang, H.; Liu, X.; Qian, Y.; Chen, J.; Li, J.; Li, Y.; et al. Estimation of Land Surface Temperature from Chinese ZY1-02E IRS Data. Remote Sens. 2024, 16, 383. https://doi.org/10.3390/rs16020383

AMA Style

Dou X, Li K, Zhang Q, Ma C, Tang H, Liu X, Qian Y, Chen J, Li J, Li Y, et al. Estimation of Land Surface Temperature from Chinese ZY1-02E IRS Data. Remote Sensing. 2024; 16(2):383. https://doi.org/10.3390/rs16020383

Chicago/Turabian Style

Dou, Xianhui, Kun Li, Qi Zhang, Chenyang Ma, Hongzhao Tang, Xining Liu, Yonggang Qian, Jun Chen, Jinglun Li, Yichao Li, and et al. 2024. "Estimation of Land Surface Temperature from Chinese ZY1-02E IRS Data" Remote Sensing 16, no. 2: 383. https://doi.org/10.3390/rs16020383

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop