Snow-Covered Area Retrieval from Himawari–8 AHI Imagery of the Tibetan Plateau

Daily snow-covered area retrieval using the imagery in solar reflective bands often encounters extensive data gaps caused by cloud obscuration. With the inception of geostationary satellites carrying advanced multispectral imagers of high temporal resolution, such as Japan’s geostationary weather satellite Himawari–8, considerable progress can now be made towards spatially-complete estimation of daily snow-covered area. We developed a dynamic snow index (normalized difference snow index for vegetation-free background and normalized difference forest–snow index for vegetation background) fractional snow cover estimation method using Himawari–8 Advanced Himawari Imager (AHI) observations of the Tibetan Plateau. This method estimates fractional snow cover with the pixel-by-pixel linear relationship of snow index observations acquired under snow-free and snow-covered conditions. To achieve reliable snow-covered area mapping with minimal cloud contamination, the daily fractional snow cover can be represented as the composite of the high temporal resolution fractional snow cover estimates during daytime. The comparison against reference fractional snow cover data from Landsat–8 Operational Land Imager (OLI) showed that the root–mean–square error (RMSE) of the Himawari–8 AHI fractional snow cover ranged from 0.07 to 0.16, and that the coefficient of determination (R2) reached 0.81–0.96. Results from the 2015/2016 and 2016/2017 winters indicated that the daily composite of Himawari–8 observations obtained a 14% cloud percentage over the Tibetan Plateau, which was less than the cloud percentage (27%) from the combination of Moderate Resolution Imaging Spectroradiometer (MODIS) onboard Terra and Aqua.