The MOD10A1 daily snow cover product offers a fine temporal resolution, which is particularly important for time-critical applications. Although the reliability of this product was mostly verified by ground measurements before use, the inadequate representativeness of point-scale observations may lead to errors in the assessment of MOD10A1. In this study, based on the Google Earth Engine, MOD10A1 was thoroughly evaluated by a total of 20,131 binary snow maps generated from 30 m Landsat imagery for the 2000-2016 snow seasons in three typical snow regions across China. The three typical snow regions included Northern Xinjiang, Northeast China and the Tibetan Plateau. In general, the accuracy of MOD10A1 under clear sky conditions is 86.5% based on ground observations from meteorological stations and 90.3% based on Landsat imagery. Compared with ground measurements, snow cover derived from Landsat images is in better agreement with MOD10A1 due to a better spatial match. The error analyses also indicate that 1) compared with land cover, topography has larger effects on the accuracy of MOD10A1 over the three snow regions with statistically significant negative relationships, particularly in Northern Xinjiang (R = -0.604, P < 0.05); and 2) MOD10A1 shows a relatively lower accuracy of 81.8% from 2006 to 2008 due to the decrease in snow depth during this period (R = 0.58, P < 0.05) based on the Mann-Kendall test.

Accurate cloud detection via satellites is important for cloud radiative forcing estimation and disaster weather monitoring. Current polar-orbiting satellite cloud observation are limited during early morning orbit and contain notable uncertainty due to dimness measurements in visible bands. FY-3E\MERSI-LL is the first early morning orbit satellite worldwide and can realize global cloud observation under early morning scenarios. In this study, a dynamic threshold cloud detection algorithm is proposed based on the FY-3E\MERSI-LL infrared channel, combined with auxiliary data such as sea surface temperature, land surface temperature, snow cover mask and terrain elevation. The algorithm can detect clouds against complex land surface background, but faces classification difficulties over some plateau, high-latitude and snow surface regions, especially during early morning observation periods. Compared to coincident Himawari-8 and GOES-16 cloud measurements in the Eastern and Western Hemispheres, respectively, our algorithm recognizes reasonable cloud distributions. Furthermore, Himawari-8 and GOES-16 cloud products are used for quantitative cloud algorithm evaluation. The results show that at low-middle latitudes (60°N-60°S), the average cloud and clear hit rates during the various seasons are 73.24% and 76.46%, respectively, the cloud leakage and false alarm rates are 14.46% and 8.15%, respectively, and the total accuracy (cloud and clear) is 77.33%. The algorithm performance is better over the ocean than over land. Ground site MPLCMASK products are also used to verify the FY-3E cloud results in middle- and high-latitude areas. This algorithm provides a cloud detection reference during early morning orbit based on infrared channels.

Accurate long-term snow-covered-area mapping is essential for climate change studies and water resource management. The NOAA AVHRR/2 provides a unique data source for long-term, large-spatial-scale monitoring of snow-covered areas at a daily scale. However, the value of AVHRR/2 in mapping snow-covered areas is limited, due to its lack of a shortwave infrared band for snow/cloud discrimination. We simulated the reflectance in the 3.75 µm mid-infrared band with a radiative transfer model and then developed three fractional-snow-cover retrieval algorithms for AVHRR/2 imagery at 1 km and 5 km resolutions. These algorithms are based on the multiple endmember spectral mixture analysis algorithm (MESMA), snow index (SI) algorithm, and non-snow/snow two endmember model (TEM) algorithm. Evaluation and comparison of these algorithms were performed using 313 scenarios that referenced snow-cover maps from Landsat-5/TM imagery at 30 m resolution. For all the evaluation data, the MESMA algorithm outperformed the other two algorithms, with an overall accuracy of 0.84 (0.85) and an RMSE of 0.23 (0.21) at the 1 km (5 km) scale. Regarding the effect of land cover type, we found that the three AVHRR/2 fractional-snow-cover retrieval algorithms have good accuracy in bare land, grassland, and Himalayan areas; however, the accuracy decreases in forest areas due to the shading of snow by the canopy. Regarding the topographic effect, the accuracy evaluation indices showed a decreasing and then increasing trend as the elevation increased. The accuracy was worst in the 4000–5000 m range, which was due to the severe snow fragmentation in the High Mountain Asia region; the early AVHRR/2 sensors could not effectively monitor the snow cover in this region. In this study, by increasing the number of bands of AVHRR/2 1 km data for fractional-snow-cover retrieval, a good foundation for subsequent long time series kilometre- resolution snow-cover monitoring has been laid.

As one of the major components in the hydrological system, seasonal snow cover in Northeast China has drawn much attention recently. Because of the coarse spatial resolution of the passive microwave (PMW), heterogeneity of snowpack, and forest cover, it is difficult for existing snow products to achieve high precision snow parameters [e.g., snow depth (SD) or snow water equivalent (SWE)] assessment and hydrological research in fine scale. In this study, a novel SD retrieval algorithm that considered both the spatiotemporal dynamic of snow characteristics and forest attenuation was developed by combining the Calibrated Enhanced-Resolution Brightness Temperature (TB) data and other auxiliary information and produced a fine resolution (i.e., 6.25 km $\times6.25$ km) and high accuracy SD data in Northeast China. Instead of complex physical models, the machine learning was used to untangle the nonlinear complex relationship between SD and the enhanced-resolution TB, forest fraction (FF), and snow characteristics. The verification results at ground weather stations showed that the retrieved SD by the proposed algorithm had high consistency with the observed SD, its RMSE, bias, and correlation coefficient ( $R$ ) of 6.32 cm, −0.23 cm, and 0.63, respectively. Compared with the existing SD products (WESTDC and AMSR2), the developed model greatly improved both spatial resolution and retrieval accuracy. In general, the fine-resolution SD inversion model achieved satisfactory accuracy and stability, and it will be used to generate a long-term SD dataset service for climate change and hydrological research in the future.

Since 2013, frequent haze pollution events in China have been attracting public attention, generating a demand to identify the haze areas using satellite observations. Many studies of haze recognition algorithms are based on observations from space-borne imagers, such as the Moderate Resolution Imaging Spectroradiometer (MODIS), the Visible Infrared Imaging Radiometer Suite (VIIRS) and the Advanced Himawari Imager (AHI). Since the haze pixels are frequently misidentified as clouds in the official cloud detection products, these algorithms mainly focus on recovering them from clouds. There are just a few studies that provide a more precise distinction between haze and clear pixels. The Medium Resolution Imaging Spectrometer II (MERSI-II), the imager aboard the FY-3D satellite, has similar bands to those of MODIS, hence, it appears to have equivalent application potential. This study proposes a novel MERSI haze mask (MHAM) algorithm to directly categorize haze pixels in addition to cloudy and clear ones. This algorithm is based on the fact that cloudy and clear pixels exhibit opposing visible channel reflectance and infrared channel brightness temperature characteristics, and clear pixels are relative brighter, and as well as this, there is a positive difference between their apparent reflectance values, at 0.865 μm and 1.64 μm, respectively, over bright surfaces. Compared with the Aqua/MODIS and MERSI-II official cloud detection products, these two datasets treat the dense aerosol loadings as certain clouds, possible clouds and possible clear pixels, and they treat distinguished light or moderate haze as possible clouds, possible clear pixels and certainly clear pixels, while the novel algorithm is capable of demonstrating the haze region’s boundary in a manner that is more substantially consistent with the true color image. Using the PM_2.5 (particle matter with a diameter that is less than 2.5 μm) data monitored by the national air quality monitoring stations as the test source, the results indicated that when the ground-based PM_2.5 ≥ 35 μg/cm³ is considered to be haze days, the samples with the recognition rate that is higher than 85% accounted for 72.22% of the total samples. When PM_2.5 ≥ 50 μg/cm³ is considered as haze days, 83.33% of the samples had an identification rate that was higher than 85%. A cross-comparison with similar research methods showed that the method proposed in this study had better sensitivity to bright surface clear and haze areas. This study will provide a haze mask for subsequent quantitative inversion of aerosol characteristics, and it will further exert the application benefits of MERSI-II instrument aboard on FY3D satellite.