Retrieval of land surface temperature from FY3D MERSI-II based on re-fitting Split Window Algorithm

ABSTRACT Medium Resolution Spectral Imager II (MERSI-II) is one of the core sensors mounted on the FengYun-3D (FY3D) satellite. Two adjacent 250 m long-wave thermal infrared (TIR) channels provide a considerable opportunity for retrieving Land Surface Temperature (LST) with high spatiotemporal resolution. In this paper, Thermodynamic Initial Guess Retrieval (TIGR) dataset and MODTRAN 4.0 model were used to re-fit the parameters of the Split-Window (SW) algorithm suitable for MERSI-II TIR channels, and then the daily 250 m resolution MERSI-II LST product was retrieved. The Radiance-based (R-based) method results showed that the bias value between simulated by MODTRAN4.0 and the input is 0.16 K, and the MAE value is 0.38 K. Inter-comparison method results showed that the MERSI-II LST and MODIS LST products were consistent in spatial distribution, but there were certain differences between MODIS LST and MERSI-II LST at different land cover types. T-based method results showed that R values between the site-observed LST and MERSI-II LST retrieved by SW algorithm exceeded 0.92, the bias value was between 3.6 K and 4.4 K, and the MAE value was between 2.6 K and 4.5 K. The above results indicating that the SW algorithm proposed in this study has good accuracy and applicability.


Introduction
Land Surface Temperature (LST) is an important physical parameter that reflects the processes of earth's surface change (Li et al., 2013 and 30 earth-atmosphere interactions (Anderson et al., 2008). It is widely used in research fields, such as the energy balance of earth's surface (Duan et al., 2014), soil moisture change (Sun et al., 2013(Sun et al., , 2012, and climate change (Sun et al., 2017). 35 Given the significant heterogeneity of land surface state parameters that affect LST, such as albedo, soil physical and thermal properties, and vegetation (Y. B. Liu et al., 2006;Neteler, 2010), LST changes rapidly in the spatio-temporal scale and 40 ground observation can hardly represent the spatio-temporal distributions of LST on the regional or global scale Prata et al., 1995). Instead, remote sensing technology offers the only possibility for measuring LST over the entire globe 45 with sufficiently high temporal resolution (Li et al., 2013;. The retrieval of the LST from remotely sensed Thermal Infrared (TIR) data has attracted much attention, and its history dates to the 1970s (McMillin, 50 1975). Over the past several decades, the algorithm of LST retrieval has been significantly improved . Many algorithms have been proposed to cater to the characteristics of various sensors onboard different satellites . Li et al. (2013) 55 summarized the retrieval algorithm of LST from remotely sensed TIR data into two groups: the first group of the LST retrieval algorithm with known land surface emissivity (LSE), which requires LSE as a priori knowledge, such as Single-channel method 60 (Jimnez- Muoz et al., 2009;Jimnez-Muoz & Sobrino, 2003, 2010Qin et al., 2001aQin et al., , 2001bZhou et al., 2010), Mulit-channel method (Becker & Li, 1990;Duan et al., 2020Duan et al., , 2017Duan et al., , 2014McMillin, 1975;Tang et al., 2008;Wan & Dozier, 1996) and Mulit-angle method (Prata, 65 1993;Sobrino et al., 2004Sobrino et al., , 1996; the second group of the LST retrieval algorithm with unknown LSE, which does not require LSE as a prior knowledge, including: (i) Stepwise retrieval methods (Ermida et al., 2020;Lan et al., 2021;Yin et al., 2020). The representatives of 70 these methods are the classification-based emissivity method (Peres & DaCamara, 2004;Sun & Pinker, 2003), the normalized difference vegetation index (NDVI)-based method (Sobrino et al., 2004(Sobrino et al., , 2008, and the temperature-independent spectral-indices 75 method (Dash et al., 2005;G.M. Jiang et al., 2006); (ii) Simultaneous LST and LSE retrieval methods with known atmospheric information. These methods can be roughly grouped into two categories: the multitemporal methods and multi (hyper)-spectral methods 80 Wan & Li, 2008);and (iii) Simultaneous retrieval of LST, LSE and atmospheric profiles. The representatives of these methods are the artificial neural network (ANN) method (Mao et al., 2011(Mao et al., , 2007(Mao et al., , 2008N. Wang et al., 2010) and the two-85 step physical retrieval method (X. L. Ma et al., 2002).
An appropriate algorithm is the key to retrieve LST with good accuracy. The SW algorithm is one of the most widely used methods at present (Hulley et al., 2019(Hulley et al., , 2021Zhu et al., 2021). It reduces 90 sensitivity to atmospheric parameters, and has the advantages of only a few input parameters, and can also maintain high retrieval accuracy. The SW algorithm is widely used in many sensors with multiple TIR bands, such as MODIS and AVHRR (Chen 95 et al., 2013;Wan & Dozier, 1996;Wan et al., 2002Wan et al., , 2004 etc. Similarly, for the FengYun-3 series satellites, the SW algorithm also shows good applicability . Jiang et al. (2011), Jiang et al. (2015 used the SW algo-100 rithm to retrieve high-precision LST products from VIRR TIR bands of FY3A and FY3B satellites respectively. MERSI-II is an integration and improvement of two imaging instruments (i.e. MERSI-1 and VIRR) 105 on the previous FY-3 (FY-3A/B/C) satellite, the SW algorithm is also suitable for MERSI-II TIR channels to retrieve good accuracy LST product H. Wang et al., 2019;Yan et al., 2021). To accurately obtain the MERSI-II LST product, this 110 study linearly simplified the Planck function based on the MERSI-II TIR channels characteristics to refit the SW algorithm coefficients suitable for MERSI-II, and used the TIGR dataset and MODTRAN4.0 model to fit the relationship between atmospheric 115 transmittance and Water Vapor Content (WVC). The objective is to obtain the atmospheric transmittance corresponding to the MERSI-II TIR channels by using WVC, and then retrieve MERSI-II LST. Lastly, atmospheric profile data, MYD11A1 data and the 120 four-component radiation data were used to evaluate the accuracy of the MERSI-II LST product using atmospheric simulate method, inter-comparison method and T-based method respectively.

Study area and data
125

Study area
The study area is in Northwest China as shown in Figure 1a, covering the coordinate range between 96°E -105.6°E, 36.45°N -43°N, including the part of Qinghai, Gansu, Ningxia, Inner Mongolia provinces 130 and part of Mongolia territory. The elevation in the study area ranges from 718 m to 5767 m, the southwest of the study area is part of the Tibetan Plateau, which is considered as the third pole of the earth, with an elevation from 3000 m to 5767 m. The black dotted 135 box shows the extent of the second largest inland river basin in China, the Heihe River Basin (HRB). From the upstream to the downstream, the HRB includes multiple natural landscapes such as ice and snow, frozen soil, forest, grassland, river, lake, oasis, desert, 140 and Gobi. The unique natural geographical environment makes HRB an important region for the study of inland river basin. At present, many scientific research and field experiments have been carried out in HRB, and accumulated abundant and valuable 145 resources, which has laid a good scientific data foundation for basin ecological research (S. M. Liu et al., 2016;Xu et al., 2017). One objective is to evaluate remote sensing models, algorithms, and products through purposeful validation experiment (Duan 150 et al., 2019;Jiang et al., 2015).

FY3D MERSI-II data
FY-3D is loaded with ten sets of advanced remote sensing sensors. MERSI-II is one of the core sensors mounted on the FY-3D satellite. It is the world's first 155 sensor that can provide global daily 250 m resolution TIR data, and can seamlessly acquire true color images with a resolution of 250 m globally every day (H. Wang et al., 2019;Yan et al., 2021). MERSI-II TIR channels have been upgraded from three infrared 160 windows of FY3A/B/C VIRR to 5. Table 1 shows the technical parameters of the MERSI-II TIR channels. Two adjacent long-wave TIR channels (24th band: 10.30-11.30 μm; 25th band: 11.50-12.50 μm) have similar interval setting with the 31st/32nd band of 165 MODIS (31st band: 10.78-11.28 μm; 32nd band: 11.77-12.27 μm), which provides better opportunity to retrieve LST based on SW algorithm (H. Wang et al., 2019;Zheng et al., 2020). Figure 2 shows a comparison diagram of the spectral response func-170 tions of the 24th and 25th bands of MERSI-II and the 31st and 32nd bands of MODIS. The comparison showed certain differences in the spectral response functions of the two TIR channels of MERSI-II and MODIS. Hence, it is inappropriate to directly use the 175 SW algorithm parameters by Qin et al. (2001b) to retrieve MERSI-II LST. In this study, we downloaded daily FY-3D MERSI-II level 1 (L1) multiband product data with a resolution of 250 m throughout the study area in 2019 to retrieve LST product.

AQUA MODIS data
MODIS is an important sensor mounted on two sunsynchronous polar-orbiting satellites, namely, TERRA and AQUA. Its product data are mainly used for longterm global observation of land surface, biosphere, 185 solid earth, atmosphere and ocean (S. Q. Yang et al., 2020;Sun et al., 2021). This study selected the LST product data from AQUA satellite (MYD11A1 V006) close to the overpass time of the FY3D satellite as the reference data, and used the inter-comparison method 190 to evaluate the MERSI-II LST retrieval accuracy. MYD11A1 V006 is the sixth-generation MODIS LST product data retrieved by Wan (2014) using the generalized SW algorithm (Wan & Dozier, 1996) and day-night algorithm (Wan & Li, 1997), with  195 spatial and temporal resolutions of 1 km and 1 d, respectively. Wan (2013) used the R-based method to verify MODIS LST using data from five bare soil sites in North Africa. The results showed that the average error of LST of 10 sets of data in 12 sets of verification 200 data was within 0.6 K, the average error of 1 set of data was below 0.8 K, and the average error of LST of another group of data was less than 1.9 K. Duan et al. (2019) used the T-based method to verify the MODIS LST data using homogeneous site measure-205 ment data. The results showed that the RMSE value between MODIS LST and measured data during the day at most sites was relatively large ( > 2K). Moreover, the RMSE value between the night MODIS LST and measured data was less than 2 K.  (Song et al., 2020. Table 2 shows the information of the nine observing stations. The measured LST was observed by the CNR1/CNR4 225 four component radiometer set-up in the automatic weather station, with averaged values provided every 10 minutes and extracted for validation at FY3D overpass time. Eq.1was used to calculate the LST. (1) where R up is the surface upward long-wave radiation, 230 R down is the atmospheric downward long-wave radiation at the surface, ε b is the surface broadband emissivity, and its calculation method is discussed in Section 3.4, and σ is the Stan-Boltzmann constant (5:67 Â 10 À8 Wm À2 K À4 ).
235 The excellent LST verification site should be a homogeneous surface with small spatial heterogeneity from point scale to kilometer scale (Coll et al., 2012). Table 3 shows the standard deviation (STD) of the extracted Landsat8 LST in the 3 Â 3 240 and 10 Â 10 windows with each site as the center. It is used to evaluate the spatial heterogeneity of the surface of the nine ground stations on the MERSI-II 250 m and MODIS 1000 m scales. Except for the YK site, the STD values of 245 Landsat8 LST in the regions of and 10 Â 10 at other sites were all below 1.06 K, it indicates that the surface spatial heterogeneity of HMZ, HEL, ZY, SDM, HRS, HZZ, and AR was small at scales of 250 m and 1000 m. Spatial heterogeneity of the 250 DSL site was small at a scale of 250 m, but its spatial heterogeneity was relatively poor at a scale of 1000 m. The surface spatial heterogeneity of YK was poor at scales of 250 m and 1000 m. Hence, YK is not suitable for evaluating the remotely 255 sensed LST product. By comparing the STD value of LST in the 3 Â 3 and 10 Â 10 windows, the STD value of LST in the 3 Â 3 window at each site was less than that in the 10 Â 10 window, thereby ideally reflecting the advantage of the LST products 260 with 250 m resolution in small-scale surface environmental monitoring.

TIGR atmospheric profile data
The TIGR atmospheric profile database is composed of 2311 typical atmospheric profiles, which 265 is widely used for the parameter simulation of LST and SST algorithms. TIGR atmospheric profile can be divided into five types: tropical atmospheric, mid-latitude summer, mid-latitude winter, polar summer, and polar winter profiles. The pressure 270 value of the TIGR atmospheric profile is set from 1013hpa to 0.05hpa with 40 pressure layers. Each layer contains elevation, humidity, temperature, and ozone information.
Given that only clear sky is considered when 275 retrieving LST, the cloud profile in the TIGR profile should be eliminated, and the criterion for determining cloud is that the relative humidity of any layer of the profile is above 90%, or that of two consecutive layers is over 85% (Galve et al., 2008;Tang et al., 2008). 280 This study combined the determination conditions and characteristics of the study area, and selected 1132 clear sky atmospheric profile data from the midlatitude summer profile as atmospheric input parameters.

Research methods
The SW Algorithm for FY3D MERSI-II LST is the main research target of thermal infrared remote sensing application. Accurately retrieving LST is significant to the study of earth-atmosphere energy 290 balance and ecosystems. Qin et al. (2001b) proposed the SW algorithm formula as follows: where T s is land surface temperature (K), T i and T j are the brightness temperatures (BT) of two thermal infrared bands with an atmospheric window in the range of 295 10-13 μm (one is near 11 μm and the other is near 12 μm), and A 0 , A 1 , and A 2 are the parameters of the SW algorithm. ccording to radiative transfer theory, radiant energy received by satellites under clear sky conditions 300 mainly comes from the earth's surface and atmosphere, and is represented as follows: where B λ T λ ð Þ is the radiance received by the sensor with a wavelength of λ, B λ T s ð Þ is the radiance of the black-body when the surface temperature is T s , ε λ is 305 the surface emissivity at the wavelength of λ, τ λ indicates atmospheric transmittance in the direction from the target to the sensor, R a" λ is the upwelling path radiances, and R a# λ is the atmospheric downwelling path radiances. Qin et al. (2001b) used Eq.4 to approx-310 imate R a" i θ ð Þ and R a# i , and held that using T a instead of T a# did not have considerable effect on the calculation result. Thus, Eq.3 can be simplified to Eq.5: To continue solving Eq.5, this study realizes the linearization of the Planck function, and constructed 315 a close linear relationship between MERSI-II TIR radiance and LST. Figure 3 shows the relationship between the radiation and LST for MERSI-II 24th and 25th bands in the LST range of 273 K to 322 K. The linear fitting results (Eq.6a and Eq.6b) showed a strong linear 320 relationship between the radiance and LST for MERSI-II TIR, and fitting accuracies R 2 were 0.9969 and 0.9978, respectively.
By incorporating the linear fitting results into Eq.6, the LST retrieval Eq.7 for the 24th and 25th bands of 325 MERSI-II can be derived as follows: The algorithm parameters are shown in Eq.8: B 24 ¼ 0:1417T 24 þ 32:717ε 24 τ 24 À 32:717: B 25 ¼ 0:1193T 25 þ 26:744ε 25 τ 25 À 26:744: (8d) Eq.8 shows that just the surface emissivity ε λ and atmospheric transmittance τ λ of the TIR channels of MERSI-II are needed to retrieve MERSI-II LST. 330 Through linearization of the Planck function, while ensuring high LST retrieval accuracy, the unknown parameters of the original LST retrieval were reduced from three to two, thereby substantially reducing the complexity of the LST retrieval.

BT of the MERSI-II TIR bands
According to the calibration coefficient, the original DN value of the TIR channels can be converted into the radiation intensity value (L λ ). Thereafter, the equivalent black-body brightness temperature value 340 (T eff ) can be calculated using the Planck equation.
Lastly, the band brightness temperature correction parameter was used to convert T eff into the band blackbody brightness temperature.
In Eq.9, gain is the gain value, offset is the offset, and 345 DN is the gray value of the image. In Eq.10, c 2 and c 1 are constants, c 2 ¼ 2hc 2 =λ 5 , and c 1 ¼ hc= Kλ ð Þ. In Eq.11, A and B are band correction coefficients derived from the attribute variables in the MERSI-II L1 file.

Estimation of surface emissivity
Land surface emissivity (LSE) is one of the important input parameters in the process of LST retrieval. The error of LSE has a significant impact on the accuracy of LST. The retrieval algorithm of LSE can be roughly 355 divided into three categories: Semi-Empirical Methods (SEMs), Multi-channel Temperature Emissivity Separation (TES), and Physically Based Methods (PBMs). SEM is the main algorithm to retrieve LSE in the process of LST retrieval (T. Song et al., 2015).
360 This study is used the method of Qin et al. (2004) to estimate the emissivity of the natural surface.
Water has a high emissivity in the thermal infrared band, which is markedly like the black-body, and LSE of the water is set to 0.995. The ε building and ε surface of 365 the two mixed pixels of urban and natural surfaces can be calculated using Eq.12 and Eq.13, respectively: where P v is the vegetation coverage, in which the calculation formula is shown in Eq.12; NDVI v and NDVI s represent the NDVI values of pure vegetation 370 (0.75) and pure bare soil (0.05) pixels, respectively; R v , R m , and R s are the temperature ratios of vegetation, buildings, and bare soil respectively, the calculation formula of which is shown in Eq.14; ε v , ε m , and ε s are the LSE of vegetation, buildings, and bare soil, respec-375 tively, in which the values of the 24th band of MERSI-II are 0.9826, 0.965, and 0.974 respectively, and the values of the 25th band of MERSI-II are 0.987, 0.975, and 0.979 respectively (H. Wang et al., 2019); and dε is the specific emissivity correction item, and its value is 380 related to the surface coverage, when P v 0:5, dε ¼ 0:0038P v ; when P v > 0:5, dε ¼ 0:0038ð1 À P v Þ; when P v ¼ 0:5, dε ¼ 0:0019; and when P v ¼ 0 or 1, dε ¼ 0 ( R Estimation of atmospheric transmittance 385 Atmospheric transmittance (τ) is also one of the important input parameters of the SW algorithm (Qin et al., 2001b;Wan et al., 2004). Given that τ is difficult to obtain directly in real time (Wan, 2008), indirect estimation of τ by establishing the relation-390 ship between WVC and τ becomes a feasible method (Mao & Qin, 2004;Mao et al., 2005). This study used 1132 TIGR clear sky atmospheric profile data as input parameters, and used MODTRAN4.0 to fit the relationship between WVC 395 and τ. Figure 4 shows the atmospheric transmittance and WVC scatter plot of the 24 th and 25 th bands of MERSI-II. The WVC values of the 1132 clear sky atmospheric profile data range from 0.06 g=cm 2 to 6.54 g=cm 2 , atmospheric transmittance fitted by 400 MODTRAN 4.0 was between 0.4 and 0.9, and atmospheric transmittance of the 24th and 25th bands of MERSI-II gradually decreased with an increase in WVC. Table 4 shows the fitting results of atmospheric transmittance and WVC of the 24th and 25th bands. The 405 accuracies of the three fitting methods have small differences. The 24th band has the lowest linear fitting accuracy (R 2 ¼ 0.974), and the cubic polynomial fitting has the highest accuracy (R 2 ¼ 0.977). R 2 of the three fitting methods of the 25th band was 0.989. Through compara-410 tive analysis, this research used cubic polynomial to fit the relationship between atmospheric transmittance in two thermal infrared bands of MERSI-II and WVC.
The 18th band in the MERSI-II sensor is the water vapor absorption zone, and its center wavelength is 415 0.940 μm. The 15th band (center wavelength is 0.865 μm) and 19th band (center wavelength is 1.03 μm) are the two adjacent atmospheric window areas. Eq.16 and Eq.17 are used to retrieve atmospheric transmittance and WVC, respectively. (Hu et al., 2011; 420 H. Q. Yang et al., 2012;X. Wang et al., 2012).

Results and analysis
MERSI-II LST accuracy verification using the 430 R-based method R-based method is a typical forward physical model, it simulates the TOA brightness temperature (BT) by inputting the initial parameters such as the satellitederived LST, atmospheric profile data and the emis-435 sivity into the atmospheric radiative transfer model ( Figure 5). During the simulation process, the input of the LST is continuously adjusted until the simulated TOA BT is consistent with the BT observed by the satellite sensor, and then the adjusted LST is consid-440 ered as the real LST. Finally, the adjusted LST is compared with the LST retrieved by remote sensing to evaluate the error of the satellite-derived LST product.
On the assumption that the LSE is between 0.91 445 and 1.0 to include most land surface types, the atmospheric profile screened out in Section 2.5 and relevant information of MERSI-II TIR channels were applied to evaluate the absolute error of the re-fitting SW algorithm by using the atmospheric 450 simulation model. The accuracy of the SW algorithm was characterized by bias, Mean Absolute Error (MAE) and Root Mean Square Error (RMSE) respectively (Eq.18, Eq.19 and Eq.20). Figure 6 shows the scatter plots of observed LST 455 versus the LST retrieved by the atmospheric simulation model, which shows that the R value between LST s and LST t is close to 1.0, the Bias value is 0.16 K, and the MAE value is 0.38 K, indicating that the re-fitting SW algorithm proposed in this 460 study has good accuracy and applicability.
(20) Table 5 shows the absolute error of SW algorithm with different values and WVC values, which indicated that SW algorithm maintains good accuracy when WVC and LST t change. The bias values between 465 LST s and LST t is range from −0.179 K to 0.364 K, the MAE values range from 0.28 K to 0.487 K, and the RMSE values range from 0.332 K to 0.552 K. In addition, the absolute error of SW algorithm increases with the increase of LST t and WVC.

MERSI-II LST accuracy verification using the inter-comparison method
The inter-comparison method compares LST retrieved by the satellite with another remotely sensed LST product with good accuracy. This method does 475 not require any ground measurement data and is a simple and feasible LST evaluation method (Coll et al., 2009;Gao et al., 2018). Figure 7 shows the spatial distribution of the MERSI-II LST(7a) and MODIS LST (7b) product on 21 October 2019. Moreover, Figure 8 480 indicates that the spatial distribution of MERSI-II LST on the north side of this study area was consistent with that of the MODIS LST products. The LST distribution map clearly distinguishes the Alxa Badain Jaran Desert. However, there were numerous low-value pix-485 els in the MERSI-II LST product on the south side of the study area. These pixels were identified as cloud pixels in the MODIS LST product. The reason is that the cloud recognition based on the CDI index in this paper failed to identify and eliminate all cloud pixels in 490 LST products. Figure 8 a and b shows the histogram and scatter plot of the differences between MERSI-II LST and MODIS LST, respectively. Figure 8a shows that the difference (ΔT) between MODIS LST and MERSI-II LST on 21 October 2019 was from −5 K to 495 7.5 K, and ΔT was centrally distributed between 2.5 K and 5.0 K. Figure 8b shows that the bias values between MERSI-II LST and MODIS LST were 3.57 K and the RMSE value was 3.89 K, indicating that the   MODIS LST values of many of the pixels were signifi-500 cantly higher than those of MERSI-II LST. Figure 9 shows the LST value changes of MODIS and MERSI-II in (a) bare soil, (b) water, (c) Forest, (d) crop, and (e) grassland with time. MODIS LST and MERSI-II LST values have a highly consistent trend 505 with time on different land cover types. Comparing Figure 9a and Figure 9e indicates that the value of ΔT of MODIS and MERSI-II in the water was the smallest, the bias value was1.15 K, followed by grassland, in which the bias value was 1.44 K. The bias value of 510 crops was 2.16 K, which is close to grassland. The bias value of Forest was 2.88 K. Meanwhile, the ΔT in time series in bare soil was the largest, with a bias value of 4.84 K between LST of the two. The results show that the value of MERSI-II LST is lower than that of 515 MODIS LST, and the ΔT between two estimated LSTs is relatively small for homogeneous land surface, such as Water, Grassland and Cropland, while the ΔT is large for heterogeneous land surface, particularly barren or sparsely vegetated land cover. 520 The inter-comparison method is usually used to evaluate the accuracy of LST product when there is a lack of site observation data, which reduces the cost of LST authenticity inspection. However, the uncertainty of the reference data itself, the difference of the 525 observation angle, observation time, or the spatial resolution between the reference data and the comparison data make the accuracy evaluation results is relative rather than absolute Zhu et al., 2021).

MERSI-II LST accuracy verification using the T-based method
The T-based verification method is used to directly compare LST obtained by ground observation with LST retrieved by remote sensing through temporal-535 spatial-angle matching. This method is the most direct LST verification method (Duan et al., 2019). Figure 10 shows the scatter plots of the MERSI-II LST retrieved by SW and GSW algorithm (respectively called LST s and LST g ) and LST measured by (a) SDM, (b) ZY, (c) AR, 540 (d) HEL, (e) HRS, (f) HZZ, (g) DSL, and (h) HMZ (LST st t). The scatter plots show that the R value between MERSI-II LST and LST st at each site exceeded 0.92, and the R value passed the 0.01-level two-sided significance test, indicating a significant correlation 545 between LST measured at the site and MERSI-II LST s .
Bias value between the LST st and MERSI-II LST s was between −3.6 K and 4.4 K, and the MAE value was between 2.6 K and 4.54 K. The R value between MERSI-II LST g and LST st at each site exceeded 0.92, except for 550 HZZ site. The bias between MERSI-II LST g and LST st was range from −6.25 K to 11.67 K. The high R value indicates that the time distribution of the MERSI-II LST product retrieved by the SW algorithm is consistent with the LST observed by the station. The relatively 555 large MAE value also indicates a large error in the MERSI-II LST g product. Table 6 shows the verification results of the MERSI-II LST retrieval accuracy under clear sky conditions. Except for March 2nd, the R values between MERSI-II 560 LST and LST measured at the HiWATER site on the other days were between 0.93 and 0.99, the absolute value of bias was between 0.06 K and 1.47 K, and the MAE value was between 1.0 K and 2.66 K. Higher R value, lower Bias and MAE indicate that MERSI-II 565 LST has higher accuracy under clear sky conditions, relatively, and can accurately characterize the thermal environment of the study area at a certain time. Figure  11 shows the overall accuracy assessment results of the eighth MERSI-II LST. R value between the eighth 570 MERSI-II LST and site-measured LST was 0.99, bias value was −0.23 K, and MAE value was 2.03 K.
Comparing Figures 10 and 11 shows that MAE value between LST measured by the site and MERSI-II LST is significantly reduced after removing some cloud pixels 575 that account for a relatively large amount of data. This result indicates that cloud shadow and thin cloud pixels have immense impact on the accuracy of LST retrieval evaluated using the T-based method.

Discussion
580 Apart from the absolute error of the algorithm, LST retrieval error is also affected by numerous factors, such as WVC, τ, and ε. This research used Eq.20 to evaluate the influence of the uncertainty of each parameter (Δx) on LST retrieval error: 585 where ΔT S is the error of LST retrieval, x is the parameter affecting the accuracy of LST retrieval, Δx is the uncertainty of the parameter, T S x þ Δx ð Þ and T S x ð Þ represent LST value inverted when the parameter takes the values of x þ Δx and x, respectively.

590
Influence of WVC on the accuracy of LST Table 7 shows LST error obtained by changing WVC uncertainty (ΔWVC). In general, the average value of ΔT S increased with an increase in ΔWVC, and the average value of ΔT S was between 0 and 0.58 K. However, 595 ΔT S does not follow this rule in all cases. When WVC was 3 g=cm 2 , ΔT S increased with a negative growth of ΔWVC. When WVC was 4.2 g=cm 2 , the value of (0.02 K) was smaller than that of ΔT S (0.11 K) when the WVC was 3.9 g=cm 2 . Figure 12 shows the influence 600 of ΔWVC on the accuracy of LST and MERSI-II's 24th and 25th band atmospheric transmittances. When ΔWVC was between −90% and 0, ΔT S , Δτ 24 , and Δτ 25 increased with an increase in ΔWVC. When was between 0 and 2.0 K, Δτ 24 and Δτ 25 were below 0.4. When 605 ΔWVC was between 0 and 90%, ΔT S did not increase with an increase in ΔWVC. When ΔT S was below 0.5, Δτ 24 and Δτ 25 were under 0.5. Kaufman and Gao (1992) confirmed that under cloud-free conditions, error of WVC retrieved by MODIS was between −13% and 610 13%. Given that MERSI-II and MODIS have similar TIR band settings, the sensitivity analysis of MERSI-II WVC can refer to the sensitivity analysis of MODIS (H. Wang et al., 2019). The preceding results indicate that LST and τ are not sensitive to WVC. Given that τ was 615 calculated based on the WVC using three polynomial fitting calculations, this situation indirectly indicates that LST inverted by the SW algorithm was not sensitive to τ.

620
The uncertainty of surface emissivity (Δε) is another key parameter of the SW algorithm. Figure 13 shows the effect of Δε on the accuracy of LST retrieval. When Δε 24 and Δε 25 changed equally between −0.0015 and 0.0015, LST retrieval error was between 0 and 0.68 K Note: ** indicates that there is a significant bilateral correlation at the level of 0.01 Figure 11. Overall accuracy evaluation results of MERSI-II LST products.  Figure 14a); when ε 25 remained unchanged (ε 25 ¼ 0:954) and Δε 24 changed from −0.0015 to 0.015, LST retrieval error was between 0 and 0.5 K (Figure 14b). When ε 24 was unchanged (ε 24 ¼ 0:944) and Δε 25 was between −0.0015 and 0.015, LST retrieval 630 error was between 0 and 1.2 K (Figure 14c). The preceding results showed that the LST retrieval error increased linearly with an increase in the absolute value of Δε. Moreover, the change of Δε 25 had a significantly greater impact on the accuracy of LST 635 retrieval than Δε 24 .
The numerical value of surface emissivity also affects the accuracy of LST retrieval. Figure 14 shows the effect of Δε on the accuracy of LST retrieval under different surface emissivity in the 24th and 25th bands 640 of MERSI-II. Under the condition that ε 25 ¼ 0:954 was unchanged, ε 24 changed from 0.944 to 1.0, when Δε 24 ¼ 0:001, LST retrieval error increased from  0.026 K to 0.035 K; when Δε 24 ¼ 0:005, LST retrieval error increased from 0.131 K to 0.177 K; and when 645 Δε 24 ¼ 0:01, LST retrieval error increased from 0.267 K to 0.359 K. Under the condition that ε 24 ¼ 0:944 was unchanged, ε 25 changed from 0.954 to 1.0; when Δε 25 ¼ 0:001, LST retrieval error decreased from 0.099 K to 0.071 K; when, LST retrieval error 650 decreased from 0.486 K to 0.350 K; and when Δε 25 ¼ 0:01, LST retrieval error decreased from 0.952 K to 0.688 K. The preceding results show that under the same uncertainty, the greater the value of ε 24 , the greater the influence of Δε 24 on LST retrieval 655 error; and the greater the value of ε 25 , the smaller the influence of Δε 25 on LST retrieval error.
Influence of cloud and snow pixels on LST Figure 8a shows there were still numerous outlier pixels in the MERSI-II LST product, and these pixel 660 values were significantly lower than LST value of the clear sky pixels. Figure 8b shows the MODIS LST product after effective cloud pixel recognition and removal. Comparing Figure 8 a and b, the low value areas of the MERSI-II LST product type correspond to  665 the cloud pixels in MODIS LST. The acquisition time of the two LST products was similar, which completely shows that the abnormal value of MERSI-II LST was the pixel affected by cloud shadow and thin cloud. In addition, snow cover is also one of the factors that 670 cannot be ignored to lead to abnormal value of MERSI-II LST pixels. By comparing Figure 8a and Figure 15, the differentiation constant pixels of MERSI-II LST product just corresponds to the snow pixel of MODIS snow cover product.

Conclusions
This study linearly simplified the Planck function based on the MERSI-II TIR band characteristics to re-fit the SW algorithm coefficients suitable for MERSI-II, and used the TIGR2000 data and 680 MODTRAN4.0 model to fit the relationship between atmospheric transmittance and WVC. The objective is to obtain the atmospheric transmittance corresponding to the MERSI-II TIR band and realize the MERSI-II LST retrieval. Lastly, atmospheric profile 685 data, MYD11A1 data and the four-component radiation data were used to evaluate the accuracy of the MERSI-II LST products using atmospheric simulate method, inter-comparison Method and T-based Method respectively. 690 The atmospheric simulate method shows that the R value between LST s and LST t is close to 1.0, the Bias value is 0.16 K, and the MAE value is 0.38 K, indicating that the re-fitting SW algorithm proposed in this study has good accuracy and applicability. The inter-695 comparison results show that the spatial distribution of MERSI-II LST and MODIS LST was relatively consistent, the difference LST (ΔT) between MODIS and MERSI-II was between −5 K and 7.5 K, and ΔT was mainly distributed between 2.5 K and 5.0 K. Bias 700 value was3.57 K, and RMSE value was 3.89 K.
T-based method results show the R value between MERSI-II LST s and the LST st at each site exceeded 0.92, and R value passed the 0.01 level two-sided significance test. Bias value was between 3.6 K and 705 4.4 K, and MAE value was between 2.6 K and 4.54 K. the R value between MERSI-II LST g and the LST st at each site exceeded 0.92, except for HZZ site, and the bias value was between −6.25 K and 11.67 K. Higher R value, lower bias, and RMSE indicate that 710 MERSI-II LST s has higher accuracy at all stages and, relatively, can accurately characterize the thermal environment of the study area at a certain time. In comparison, the accuracy of LST g is lower than that of LST s .

Disclosure statement
No potential conflict of interest was reported by the author(s).