An Assessment of Existing Methodologies to Retrieve Snow Cover Fraction from MODIS Data

The characterization of snow extent is critical for a wide range of applications. Since 1966, snow maps at different spatial resolutions have been produced using various satellite sensor images. Nowadays, the most widely used products are likely those derived from Moderate-Resolution Imaging Spectroradiometer (MODIS) data, which cover the whole Earth at a near-daily frequency. There are a variety of snow mapping methods for MODIS data, based on different methodologies and applied at different spatial resolutions. Up to now, all these products have been tested and evaluated separately. This study aims to compare the methods currently available for retrieving snow from MODIS data. The focus is on fractional snow cover, which represents the snow cover area at the subpixel level. We examine the two main approaches available for generating such products from MODIS data; namely, linear regression of the Normalized Difference Snow Index (NDSI) and spectral unmixing (SU). These two approaches have resulted in several methods, such as MOD10A1 (the NSIDC MODIS snow product) for NDSI regression, and MODImLAB for SU. The assessment of these approaches was carried out using higher resolution binary snow maps (i.e., showing the presence or absence of snow) at spatial resolutions of 10, 20, and 30 m, produced by SPOT 4, SPOT 5, and LANDSAT-8, respectively. Three areas were selected in order to provide landscape diversity: the French Alps (117 dates), the Pyrenees (30 dates), and the Moroccan Atlas (24 dates). This study investigates the impact of reference maps on accuracy assessments, and it is suggested that NDSI-based high spatial resolution reference maps advantage NDSI medium-resolution snow maps. For MODIS snow maps, the results show that applying an NDSI approach to accurate surface reflectance corrected for topographic and atmospheric effects generally outperforms other methods for the global retrieval of snow cover area. The improvements to the newer version of MOD10A1 (Collection 6) compared to the older version (Collection 5) are significant. Products based on SU provide a good alternative and more accurate retrieval of the snow fraction where wider ranges of land covers are concerned. The fusion process and its resulting 250 m spatial resolution product improve snow line retrieval. False detection in mixed pixels, probably due to the spectral variability associated with the various materials in the spectral mixture, has been identified as an area that will require improvement.