Assessing FY-3E HIRAS-II Radiance Accuracy Using AHI and MERSI-LL

: The FY-3E/HIRAS-II (Hyperspectral Infrared Atmospheric Sounder-II), as an infrared hyperspectral instrument onboard the world’s ﬁrst early morning polar-orbiting satellite, plays a major role in improving the accuracy and timeliness of global numerical weather predictions. In order to assess its observation quality, the geometrically, temporally, and spatially matched scene homogeneous HIRAS-II hyperspectral observations were convolved to the channels corresponding to the Himawari-8/AHI (Advanced Himawari Imager) and FY-3E/MERSI-LL (Medium-Resolution Spectral Imager) imagers from 15 March to 21 April 2022, and their brightness temperature deviation characteristics were statistically calculated in this paper. The results show that the HIRAS-II in-orbit observed brightness temperatures are slightly warmer than the AHI observations in all the matched AHI channels (long wave infrared channel 8 to channel 16) with a mean brightness temperature bias less than 0.65 K. The bias of the atmospheric absorption channel is slightly larger than that of the window channel. A standard deviation less than 0.31 K and a correlation coefﬁcient higher than 0.98 in all channels means that the quality of the observation is satisfactory. The thresholds chosen for the colocation approximation factors (e.g., observation geometry angle, scene uniformity, observation azimuth, and observation time) for matching the HIRAS-II with AHI contribute little and negligible uncertainty to the bias assessment, so the difference between the two observed radiations is considered to be mainly from the systematic bias of the two-instrument measurement. Compared with MERSI-LL window channel 5, the observations of both instruments are very close, with a mean bias of 0.002 K and a standard deviation of 0.31 K. The mean brightness temperature bias (HIRAS-II minus MERSI-LL) of the MERSI-LL water vapor channel 4 is 0.66 K with a standard deviation of 0.22 K. The mean brightness temperature bias of channel 6 and channel 7 is 0.63 K (the standard deviation is 0.36 K) and 0.5 K (the standard deviation is 0.3 K), respectively. The biases of channel 4 are signiﬁcantly and positively correlated with the target scene temperature, and the biases of channel 6 and 7 show a U-shaped change with the increase in the scene temperature, and the biases are smallest (close to 0 K) when the scene temperature is between 250 K and 280 K. The statistical characteristics of the HIRAS-II–MERSI-LL difference vary minimally and almost constantly over a period of time, indicating that the performance of the HIRAS-II instrument is stable.


Introduction
Satellite-borne infrared hyperspectral atmospheric sounders can obtain global meteorological observations with high precision and high spectral resolution and have been frequently applied to retrieving atmospheric temperature and humidity profiles, data assimilation, and climate studies [1]. In addition, quality control of meteorological satellite observations is a pivotal step before using satellite data for assimilation and retrieval, and it is the fundamental basis for building long-term infrared hyperspectral datasets [2].
FY-3E, the world's first early morning polar-orbiting meteorological satellite, was successfully launched on 5 July 2021. It effectively fills a void for global satellite observations

Data Used in the Research
The HIRAS-II is an interferometric Fourier transform spectrometer carried in a polar orbit 836 km above the ground. HIRAS-II views the ground in the conventional mode through a cross-track rotary scan mirror that provides ±50.4 • ground coverage every 8 s. Each scan line observes 32 fields of regard (FORs), including 28 continuous Earth targets, 2 cold space targets, and 2 blackbody targets on the satellite. Each field of regard (FOR) includes a 3 × 3 field of view (FOV) with a spatial resolution of 14 km at the nadir. HIRAS-II covers the 3.92-15.38 µm infrared band with 3041 continuous channels at a spectral resolution of 0.625 cm −1 . The HIRAS-II Level 1 radiance observations from 15  April 2022 are used in this paper. The data can be found on the Chinese Feng Yun satellite remote sensing data service network (http://data.nsmc.org.cn accessed on 11 April 2022).
The Moderate-Resolution Spectral Imager-Low Light (MERSI-LL) is an important optical instrument onboard FY-3E with microlight and infrared detection capabilities. It is equipped with one visible channel operable with low-level illumination and six infrared channels. The spatial resolution of the two infrared split-window channels is 250 m, and the remaining channels are 1000 m. The MERSI-LL Level 1 radiance observations with a spatial resolution of 1000 m from 15 to 22 March 2022 are used in this paper and can be downloaded from http://data.nsmc.org.cn (accessed on 11 April 2022).
The AHI on the geostationary satellite Himawari-8 successfully launched in October 2014 and is one of the most advanced spaceborne imagers in the world. It has 16 observation channels (3 visible, 3 near-infrared, and 10 infrared), of which the spatial resolution of the infrared channel is 2 km and the temporal resolution is 10 min. Himawari-8/AHI radiation data obtained from the Japan Earth Observation Data Center (https://www.eorc.jaxa.jp/ ptree/index.html accessed on 11 April 2022) from 15 March to 21 April 2022 are analyzed.
The channel settings and performance of the thermal infrared band covered by HIRAS-II, AHI, and MESI-LL are shown in Table 1. The last row of the table (spectral coverage) specifies the central wavelength and the corresponding peak height of the weighting function (in parentheses) for each channel. The AHI has nine channels that can be completely spectrum matched with the HIRAS-II spectrum, of which channels 8, 9, and 10 are water vapor absorption channels; channels 13, 14, and 15 are window channels; and channels 11, 12, and 16 are SO 2 , O 3 , and CO 2 absorption channels, respectively. The weighting function peak heights of AHI channels 8, 9, 10, 12, and 16 are 300 hPa, 371 hPa, 532 hPa, 40 hPa, and 863 hPa, respectively. The weighting function heights of the remaining channels are almost near the surface. For MERSI-LL, only channels 4, 5, 6, and 7 can be completely spectrum matched. Channel 4 is a water vapor channel (the peak height of the weighting function is approximately 400 hPa), and channels 5, 6, and 7 are window channel with a central wavelength of 8.55 µm, 10.8 µm, and 12.0 µm, respectively. matching process is shown in Figure 1. The following will be a step-by-step detailed introduction.

Spectral Convolution
The most fundamental problem in the intercomparison between the spaceborne imager and sounder observation is radiance spectral matching. The two types of instruments with different spectral resolutions cannot be directly compared. To be compared with the imager instrument observations, the hyperspectral radiance must be convolved to match the spectral response function (SRF) of the broadband imager [14,15]. The computational formula to achieve spectral convolution is where ν is the wavenumber, L is the radiance calculated by convolution, R(ν) is the hyperspectral radiance at the corresponding wavenumber and S(ν) is the spectral response function of the imager. The simulated HIRAS-II brightness temperature spectrum using fast radiative transfer mode RTTOV under American standard atmospheric conditions, as well as the AHI and MERSI-LL spectral response functions (SRF), are given in Figure 2. The bottom horizontal coordinate is the wavelength, the upper horizontal coordinate is the wavenumber, the left vertical coordinate is the simulated brightness temperature of HIRAS-II, and the right vertical coordinate is the spectral response function. The solid line in the figure is the spectral response function of the matched AHI channels 8-16, and the dashed line denotes the spectral response function of MERSI channels 4-7.

Spectral Convolution
The most fundamental problem in the intercomparison between the spaceborne imager and sounder observation is radiance spectral matching. The two types of instruments with different spectral resolutions cannot be directly compared. To be compared with the imager instrument observations, the hyperspectral radiance must be convolved to match the spectral response function (SRF) of the broadband imager [14,15]. The computational formula to achieve spectral convolution is where ν is the wavenumber, L is the radiance calculated by convolution, R(ν) is the hyperspectral radiance at the corresponding wavenumber and S(ν) is the spectral response function of the imager. The simulated HIRAS-II brightness temperature spectrum using fast radiative transfer mode RTTOV under American standard atmospheric conditions, as well as the AHI and MERSI-LL spectral response functions (SRF), are given in Figure 2.
The bottom horizontal coordinate is the wavelength, the upper horizontal coordinate is the wavenumber, the left vertical coordinate is the simulated brightness temperature of HIRAS-II, and the right vertical coordinate is the spectral response function. The solid line in the figure is the spectral response function of the matched AHI channels 8-16, and the dashed line denotes the spectral response function of MERSI channels 4-7.

Observation Geometry and Temporal Matching
The most important step in the intercomparison between the imager and the sounder observation is to find the consistent (same) FOV [16]. The polar-orbiting satellite FY-3E passes the geostationary satellite Himawari's nadir 140.7 E at approximately 0830 or 2030 UTC daily, the observation time difference remained within 10 min. Since the measurements of the cross-track scanning instruments are sensitive to the scan angle (satellite zenith angle), the scan angle of both instruments in the matched field of view is specified to be less than 5 • to ensure that HIRAS-II and AHI observe the same scene. To Remote Sens. 2022, 14, 4309 5 of 13 minimize the difference in observation geometry and ensure that both instruments have similar observation geometry paths, the scan angle is further constrained [7,17] cos θ 1 cos θ 2 − 1 < 0.002 (2) where θ 1 and θ 2 are the scan angles of the geostationary and polar-orbiting satellites, respectively.

Observation Geometry and Temporal Matching
The most important step in the intercomparison between the imager and the sounder observation is to find the consistent (same) FOV [16]. The polar-orbiting satellite FY-3E passes the geostationary satellite Himawari's nadir 140.7 E at approximately 0830 or 2030 UTC daily, the observation time difference remained within 10 min. Since the measurements of the cross-track scanning instruments are sensitive to the scan angle (satellite zenith angle), the scan angle of both instruments in the matched field of view is specified to be less than 5° to ensure that HIRAS-II and AHI observe the same scene. To minimize the difference in observation geometry and ensure that both instruments have similar observation geometry paths, the scan angle is further constrained [7,17] cos θ 1 cos θ 2 − 1 <0.002 (2) where θ 1 and θ 2 are the scan angles of the geostationary and polar-orbiting satellites, respectively.

Spatial Matching and Uniformity Check
The spatial field of view collocation of the FY-3E/HIRAS-II and Himawari-8/AHI in the infrared channel is shown in Figure 3, where the grid is the AHI pixel, the large circle is a HIRAS-II field of regard (FOR) centered on the field of view (FOV) to be matched (target) and the small circle is the HIRAS-II FOV within its FOR. The spatial resolution of the HIRAS-II FOV (14 km at the nadir) is coarser than that of the AHI infrared channel (2 km). Approximately 7 × 7 AHI pixels (small squares in Figure 3) fall in a HIRAS-II FOV (the small circle in Figure 3). The average of all these AHI pixel measurements is taken as the imager's measurement in the matching FOV [18]. This requires that the observation target be relatively homogeneous. Since there are many AHI pixels collocated in the HIRAS-II FOV, some of which are clear sky and some are cloudy, at the same time the underlying surface of the field of view is not homogeneous, it is necessary to check the uniformity of each HIRAS-II FOV, which is carried out by using the ratio of the standard deviation of the matched AHI radiations to its mean value where Std fov denotes the standard deviation of the observed radiance from all the matched AHI pixels within each HIRAS-II FOV, and Mean fov is the mean radiance of all

Spatial Matching and Uniformity Check
The spatial field of view collocation of the FY-3E/HIRAS-II and Himawari-8/AHI in the infrared channel is shown in Figure 3, where the grid is the AHI pixel, the large circle is a HIRAS-II field of regard (FOR) centered on the field of view (FOV) to be matched (target) and the small circle is the HIRAS-II FOV within its FOR. The spatial resolution of the HIRAS-II FOV (14 km at the nadir) is coarser than that of the AHI infrared channel (2 km). Approximately 7 × 7 AHI pixels (small squares in Figure 3) fall in a HIRAS-II FOV (the small circle in Figure 3). The average of all these AHI pixel measurements is taken as the imager's measurement in the matching FOV [18]. This requires that the observation target be relatively homogeneous. Since there are many AHI pixels collocated in the HIRAS-II FOV, some of which are clear sky and some are cloudy, at the same time the underlying surface of the field of view is not homogeneous, it is necessary to check the uniformity of each HIRAS-II FOV, which is carried out by using the ratio of the standard deviation of the matched AHI radiations to its mean value where Std fov denotes the standard deviation of the observed radiance from all the matched AHI pixels within each HIRAS-II FOV, and Mean fov is the mean radiance of all these AHI pixels. The threshold value is set to 0.01. Only uniform scenes are selected for intercomparison to reduce the uncertainty introduced by the field of view averaging. Furthermore, a constraint of background environmental uniformity of the HIRAS-II FOV is needed to compensate for the minor error of the spatial collocation, as well as to reduce the uncertainty due to different azimuths. The 3 × 3 HIRAS-II FOV (the large circle in Figure 3) centered on the HIRAS-II FOV to be matched is considered the background environment area.
Std env /Mean env < 0.05 (4) where Std env denotes the standard deviation of all the matched AHI pixels' radiance within the background area, and Mean env denotes the mean radiance of these matched AHI pixels. in Figure 3) centered on the HIRAS-II FOV to be matched is considered the background environment area.
Std env /Mean env < 0.05 (4) where Std env denotes the standard deviation of all the matched AHI pixels' radiance within the background area, and Mean env denotes the mean radiance of these matched AHI pixels.

Statistical Calculation
For each collocated HIRAS-II-AHI FOV, the brightness temperature difference (BT diff=BT HIRAS − BT AHI ) and the standard deviation of the brightness temperature dif- ) are counted (n is the number of samples).
A total of 458 pairs of samples are obtained after collocating the HIRAS-II observations with the AHIs from 15 March to 21 April 2022 based on the above matching steps.

Matching with MERSI-LL on the Same Platform
When the HIRAS-II observations are compared with those from the AHI nadir on the geostationary satellite, the matched fields of view are concentrated near the tropical equator, and the dynamic range of the observed brightness temperature is narrow. The HIRAS-II observations will be compared with measurements from one imager that is carried on a polar-orbiting platform to evaluate its observation accuracy on a global scale. However, FY-3E is the first early morning polar-orbiting satellite with a significant observation time difference (even more than 8 h) from the established mid-morning or afternoon orbit satellites. Therefore, the MERSI-LL imager on the same FY-3E platform is chosen to perform the HIRAS-II calibration in this paper, and the MERSI-LL channels 4-7 are spectral matched. Since HIRAS-II and MERSI-LL are on the same polar-orbiting satellite platform and are almost observed simultaneously, the matching process of the two instruments is relatively simple. Only the nadir HIRAS-II FOV is matched to ensure the same observation scene. The spectral matching, spatial matching, and uniformity checking steps are the same as those for AHI. Finally, a total of 12,395 pairs of HIRAS-II and MERSI-LL observation samples were matched over 8 days from 15 to 22 March 2022.

Statistical Calculation
For each collocated HIRAS-II-AHI FOV, the brightness temperature difference (BT diff = BT HIRAS − BT AHI ) and the standard deviation of the brightness temperature A total of 458 pairs of samples are obtained after collocating the HIRAS-II observations with the AHIs from 15 March to 21 April 2022 based on the above matching steps.

Matching with MERSI-LL on the Same Platform
When the HIRAS-II observations are compared with those from the AHI nadir on the geostationary satellite, the matched fields of view are concentrated near the tropical equator, and the dynamic range of the observed brightness temperature is narrow. The HIRAS-II observations will be compared with measurements from one imager that is carried on a polar-orbiting platform to evaluate its observation accuracy on a global scale. However, FY-3E is the first early morning polar-orbiting satellite with a significant observation time difference (even more than 8 h) from the established mid-morning or afternoon orbit satellites. Therefore, the MERSI-LL imager on the same FY-3E platform is chosen to perform the HIRAS-II calibration in this paper, and the MERSI-LL channels 4-7 are spectral matched. Since HIRAS-II and MERSI-LL are on the same polar-orbiting satellite platform and are almost observed simultaneously, the matching process of the two instruments is relatively simple. Only the nadir HIRAS-II FOV is matched to ensure the same observation scene. The spectral matching, spatial matching, and uniformity checking steps are the same as those for AHI. Finally, a total of 12,395 pairs of HIRAS-II and MERSI-LL observation samples were matched over 8 days from 15 to 22 March 2022.

Comparison of HIRAS-II with AHI
There are nine Himawari-8/AHI channels (channels 8 to 16) that are spectral matched with FY-3E/HIRAS-II in the longwave infrared band. The scatter plots of the convolved HIRAS-II observed brightness temperature with AHI measurements in channel 8 to 16 are given sequentially in Figure 4. The horizontal coordinate is the HIRAS-II observed brightness temperature, the vertical axis is the AHI observation, the dashed line is the y = x line, and the solid line is the linear fitting result. Because the matched samples are concentrated near the nadir of the geostationary satellite, the dynamic range of the observed brightness temperature for each channel is narrow, and the value gradually decreases as the peak height of the weighting functions of atmospheric absorption channels increases. The observations of the two instruments are very close, with correlation coefficients higher than 0.98, and the fitting lines almost coincide with the y = x line in all the channels. The statistical bias and standard deviation for the HIRAS-II-AHI matching samples from 15 March to 21 April 2022 are listed in Table 2. Figure 4 and Table 2 show that all the channels have a slightly positive bias; namely, that the HIRAS-II convolved observations are slightly warmer than AHI with a maximum bias of 0.65 K (channel 9 in the water vapor wing), and the minimum is 0. 22 K (window channel 14). The standard deviation of all the channel biases ranges from 0.22 to 0.31 K with small values and little difference between the channels. Water vapor absorption channels 8-10 and ozone absorption channel 12 have relatively larger biases and small standard deviations, while the biases of the window channels (such as channels 14 and 15) are relatively small, and the standard deviations are slightly large. In addition, the closer the peak height of the weighting function is to the surface, the larger the standard deviation is. This is because the value range of the observed brightness temperature of the window channels is relatively larger than that of the absorption channels, so the dispersion is larger in the window channel.
given sequentially in Figure 4. The horizontal coordinate is the HIRAS-II observed brightness temperature, the vertical axis is the AHI observation, the dashed line is the y = x line, and the solid line is the linear fitting result. Because the matched samples are concentrated near the nadir of the geostationary satellite, the dynamic range of the observed brightness temperature for each channel is narrow, and the value gradually decreases as the peak height of the weighting functions of atmospheric absorption channels increases. The observations of the two instruments are very close, with correlation coefficients higher than 0.98, and the fitting lines almost coincide with the y = x line in all the channels. The statistical bias and standard deviation for the HIRAS-II-AHI matching samples from 15 March to 21 April 2022 are listed in Table 2. Figure 4 and Table 2 show that all the channels have a slightly positive bias; namely, that the HIRAS-II convolved observations are slightly warmer than AHI with a maximum bias of 0.65 K (channel 9 in the water vapor wing), and the minimum is 0. 22 K (window channel 14). The standard deviation of all the channel biases ranges from 0.22 to 0.31 K with small values and little difference between the channels. Water vapor absorption channels 8-10 and ozone absorption channel 12 have relatively larger biases and small standard deviations, while the biases of the window channels (such as channels 14 and 15) are relatively small, and the standard deviations are slightly large. In addition, the closer the peak height of the weighting function is to the surface, the larger the standard deviation is. This is because the value range of the observed brightness temperature of the window channels is relatively larger than that of the absorption channels, so the dispersion is larger in the window channel.   The ideal condition for the realization of the cross-calibration using the simultaneous nadir overpass (SNO) method is that the two instruments observe the same target at the same time through the same atmospheric path. However, the matching condition is ap- The ideal condition for the realization of the cross-calibration using the simultaneous nadir overpass (SNO) method is that the two instruments observe the same target at the same time through the same atmospheric path. However, the matching condition is appropriately approximated in the actual application to obtain enough samples, and the introduction of the matching threshold may bring uncertainty to the calibration evaluation. For the HIRAS-II and AHI cross-calibration, the observation inconsistency may be partly due to the random errors caused by the differences in observation geometry, scene  Figure 5. The larger the value of the horizontal coordinate in Figure 5a denotes the larger the difference in the observation geometric viewpoint. A larger value of the x-axis in Figure 5b indicates worse uniformity within the field of view. The larger x value of Figure 5c corresponds to the larger difference in the observation azimuth of the two instruments. The solid line shows the linear fit pattern of the brightness temperature differences with the approximation factors. The brightness temperature differences are randomly and uniformly distributed with these factors and do not increase with the decrease in various matching degrees. There is no obvious linear variation characteristic with these matching factors, indicating that the influence of various matching factor differences within their threshold on observation bias can be neglected and that these reasonable matching thresholds bring little uncertainty to the bias assessment. Figure 5e shows the distribution of the brightness temperature differences with the HIRAS-II observed brightness temperature. These is also no significant scene temperature-dependent bias. The results of the remaining channels are similar and omitted. After excluding the random errors caused by the matching factors mentioned above, it can be concluded that the HIRAS-II-AHI brightness temperature differences mainly represent the systematic observation bias of the two instruments.

Comparison of HIRAS-II with MERSI-LL
The observed biases and standard deviations were counted based on the 12,395 pairs of samples matched by HIRAS-II and MERSI from 15-22 March 2022. The distribution of the brightness temperature differences (HIRAS-II minus MERSI-LL) with the HIRAS-II observed scene temperature for channel 4-7 are shown in the left subplot of Figure 6 with the vertical coordinates representing the brightness temperature differences and the horizontal coordinates representing the observed scene temperature. The color distinguishes

Comparison of HIRAS-II with MERSI-LL
The observed biases and standard deviations were counted based on the 12,395 pairs of samples matched by HIRAS-II and MERSI from 15-22 March 2022. The distribution of the brightness temperature differences (HIRAS-II minus MERSI-LL) with the HIRAS-II observed scene temperature for channel 4-7 are shown in the left subplot of Figure 6 with the vertical coordinates representing the brightness temperature differences and the horizontal coordinates representing the observed scene temperature. The color distinguishes the scene uniformity, and the dashed line gives the mean of the biases. The probability density distributions of the brightness temperature differences are given in the right subplot with the horizontal coordinate as the sample probability density. Channel 4 of MERSI-LL is a water vapor absorption band with a central wavelength of 7.22 µm (peak height of the weighting function 400 hPa), and channel 5-7 are window channels with a central wavelength of 8.55 µm, 10.8 µm, and 12.0 µm, respectively. The dynamic range of the channel 4 target brightness temperature is between 220 and 280 K, and the HIRAS-II measurement is slightly higher than the MERSI-LL observations (the mean bias is 0.6643 K) with a standard deviation of 0.2229 K. The dynamic range of channel 5 is slightly larger at approximately 220-300 K, and the observed brightness temperatures of HIRAS-II and MERSI are close, with a mean bias of 0.0023 K and a standard deviation of 0.3135 K. The dynamic range of channel and 7 target brightness temperature are between 210 and 310 K, and the biases of channel 6 and 7 show a U-shaped change with the increase in the scene temperature, and the biases are smallest (close to 0 K) when the scene temperature is between 250 K and 280 K. Channel 4 has smaller brightness temperature differences at lower scene temperatures (i.e., high latitudes) and relatively larger brightness temperature differences at higher scene temperatures (i.e., low latitudes), especially when the bias increases to 1.2 K near the equator. Although the channel 5 bias takes values of approximately 0 K, the bias dispersion increases as the HIRAS-II observed scene temperature increases. Both channel 6 and 7 have relatively larger brightness temperature differences at lower scene temperatures and higher scene temperatures, and the maximum value is close to 1.75 K. From the right subplot, it can be seen that the probability density distributions of the brightness temperature bias for channels 4-7 all conform to the normal distribution.
It is noteworthy that the brightness temperature differences of the water vapor channel 4 in Figure 6 are obviously positively correlated with the target scene temperature, and the window channels 6 and 7 also have an obvious scene temperature-dependence, while window channel 5 shows no scene temperature-dependent bias. At the same time, AHI water vapor channels 9 and 10-whose spectral positions are close to MERSI-LL channel 4-also do not find bias scene-dependent characteristics. Since HIRAS-II and MERSI-LL are mounted on the same platform, the scene uniformity is the only factor that introduces matching uncertainty into the intercomparison. Figure 7 shows the scatter distribution of MERSI-LL channel 4 (a), channel 5 (b), channel 6 (c), and channel 7 (d) brightness temperature differences (HIRAS-II minus MERSI-LL) with scene uniformity. A larger value of the horizontal coordinate in Figure 7 indicates worse scene uniformity, and the solid line indicates the linear fitting result. The brightness temperature differences of channel 5-7 are uniformly distributed with the scene uniformity and do not have linear variation characteristics (in Figure 7b). However, the brightness temperature differences of channel 4 show an obvious linear trend with the scene uniformity, and the biases gradually decrease as the scene uniformity worsens (in Figure 7a). Combined with the scatter color of channel 4 in Figure 6, the scene uniformity is relatively poor (yellow) in the high latitudes with a low brightness temperature, and the scene uniformity is good in the low latitudes with a high brightness temperature. This is because the underlying surface in the field of view varies greatly in the polar region when the instrument is scanning with the same spatial resolution and swath, especially the Arctic has greater underlying surface variability due to the presence of different surface types (e.g., land, snow, ocean, glacier, etc.) with higher variability in absolute temperature. Theoretically, the bias is smaller when scene uniformity is better. However, Figure 6 shows that the scene uniformity gradually improves with the increasing scene temperature, while the bias increases instead. This indicates that the scene uniformity is not the cause of the scene temperature-dependent bias.

Conclusions
To assess HIRAS-II's on-orbit observation quality, the geometrically, temporally, and spatially matched scene homogeneous HIRAS-II hyperspectral observations were convolved to the longwave infrared channels corresponding to the Himawari-8/AHI and FY-3E/MERSI-LL from 15 March to 21 April 2022, and their brightness temperature deviation characteristics were statistically calculated in this paper. The matching samples of HIRAS-II and AHI are concentrated near the equator, and the spectral matching channels are longwave infrared channel 8 to channel 16 onboard the same polar orbiting satellite platform FY-3E. The matching samples of HIRAS-II and MERSI-LL are evenly distributed all over the world with spectral matching channels 4 to channel 7. The following conclusions can be made based on this analysis: 1. The HIRAS-II on-orbit observed brightness temperatures are slightly warmer than the AHI observations, with a small positive bias in all the matched channels. The brightness temperature scatters of both observations are distributed near the y = x line with a correlation coefficient higher than 0.98 in all channels. The biases of water The main factors that may cause spaceborne radiation imager calibration errors mainly include blackbody emissivity and spectral response function instrument nonlinearity. MERSI-LL channel 4 is located in the wing area of the water vapor absorption band, and a very small drift in the spectral response function can also lead to a temperaturedependent bias in the scene. However, to date, there have been no specific references about the design of the black bodies, the calibration system, and so on of these two instruments onboard FY-3E. These matters require further study in the future.
The day-to-day variations in the mean biases and standard deviations are counted based on the HIRAS-II-MERSI-LL matched samples. The daily mean biases of channel 4 range from 0.64 to 0.68 K, and the standard deviations are all approximately 0.2 K. The daily mean biases of channel 5 range from −0.02 to 0.02 K, and the standard deviations are approximately 0.3 K. The daily mean biases of channel 6 range from 0.62 K to 0.64 K, and the standard deviations are approximately 0.35 K. The daily mean biases of channel 7 range from 0.48 K to 0.51 K, and the standard deviations are approximately 0.3 K. The biases of the two instruments vary minimally and almost constantly over a period of time, indicating that the performance of the HIRAS-II instrument is stable.

Conclusions
To assess HIRAS-II's on-orbit observation quality, the geometrically, temporally, and spatially matched scene homogeneous HIRAS-II hyperspectral observations were convolved to the longwave infrared channels corresponding to the Himawari-8/AHI and FY-3E/MERSI-LL from 15 March to 21 April 2022, and their brightness temperature deviation characteristics were statistically calculated in this paper. The matching samples of HIRAS-II and AHI are concentrated near the equator, and the spectral matching channels are longwave infrared channel 8 to channel 16 onboard the same polar orbiting satellite platform FY-3E. The matching samples of HIRAS-II and MERSI-LL are evenly distributed all over the world with spectral matching channels 4 to channel 7. The following conclusions can be made based on this analysis:

1.
The HIRAS-II on-orbit observed brightness temperatures are slightly warmer than the AHI observations, with a small positive bias in all the matched channels. The brightness temperature scatters of both observations are distributed near the y = x line with a correlation coefficient higher than 0.98 in all channels. The biases of water vapor channels 8-10 and ozone absorption channel 12 are relatively large, with a maximum of 0.65 K (channel 9 in the water vapor wing), and the biases of the window channels are relatively small, with a minimum of 0.22 K (channel 14). The standard deviations for all channels are small (0.22-0.31 K) and there is little difference between the channels.

2.
The thresholds chosen for the colocation approximation factors (e.g., observation geometry angle, field of view uniformity, observation azimuth, and observation time) when matching the HIRAS-II with AHI contribute little and negligible uncertainty to the bias assessment, so the difference between the two observed radiations is considered to be mainly from the systematic bias of the two instrument measurements. 3.
Since HIRAS-II and MERSI-LL are mounted on the same platform, the scene uniformity is the only factor introducing matching uncertainty in the intercomparison. The mean brightness temperature bias (HIRAS-II minus MERSI-LL) of the MERSI-LL water vapor channel 4 is 0.66 K with a standard deviation of 0.22 K. To window channel 5, the observations of both instruments are very close, with a mean bias of 0.002 K and a standard deviation of 0.31 K. Both channel 6 and 7 have relatively larger brightness temperature differences at lower scene temperatures and higher scene temperatures, with a mean bias of 0.63 K (the standard deviation is 0.36 K) and 0.5 K (the standard deviation is 0.3 K), respectively. 4.
The biases of MERSI-LL channel 4 are obviously positively correlated with the target scene temperature. The biases of channel 6 and 7 show a U-shaped change with the increase in the scene temperature, and the biases are smallest (close to 0 K) when the scene temperature is between 250 K and 280 K. The statistical characteristics of the HIRAS-II-MERSI-LL difference vary minimally and almost constantly over a period of time, indicating that the performance of the HIRAS-II instrument is stable.
As a final note, we just found the phenomenon of bias distribution, which is not yet fully explained due to lack of relevant references. Therefore, we will use NWP data, double-difference method to further evaluate the accuracy of HIRAS-II in future studies.

Data Availability Statement:
The HIRAS-II and MERSI-LL Level 1 data can be obtained at the Chinese Feng Yun satellite remote sensing data service network (http://data.nsmc.org.cn accessed on 11 April 2022). Himawari-8/AHI radiation data obtained from the Japan Earth Observation Data Center (https://www.eorc.jaxa.jp/ptree/index.html accessed on 11 April 2022).