Landsat, MODIS, and VIIRS snow cover mapping algorithm performance as validated by airborne lidar datasets

Snow cover mapping algorithms utilizing multispectral satellite data at various spatial resolutions are available, each treating subpixel variation differently. Past evaluations of snow mapping accuracy typically relied on satellite data collected at a higher spatial resolution than the data in ques...

Full description

Bibliographic Details
Published in:The Cryosphere
Main Authors: T. Stillinger, K. Rittger, M. S. Raleigh, A. Michell, R. E. Davis, E. H. Bair
Format: Article in Journal/Newspaper
Language:English
Published: Copernicus Publications 2023
Subjects:
Online Access:https://doi.org/10.5194/tc-17-567-2023
https://doaj.org/article/5f70dc89f8e94b4da40ca56b5134f757
Description
Summary:Snow cover mapping algorithms utilizing multispectral satellite data at various spatial resolutions are available, each treating subpixel variation differently. Past evaluations of snow mapping accuracy typically relied on satellite data collected at a higher spatial resolution than the data in question. However, these optical data cannot characterize snow cover mapping performance under forest canopies or at the meter scale. Here, we use 3 m spatial resolution snow depth maps collected on 116 d by an aerial laser scanner to validate band ratio and spectral-mixture snow cover mapping algorithms. Such a comprehensive evaluation of sub-canopy snow mapping performance has not been undertaken previously. The following standard (produced operationally by an agency) products are evaluated: NASA gap-filled Moderate Resolution Imaging Spectroradiometer (MODIS) MOD10A1F, NASA gap-filled Visible Infrared Imaging Radiometer Suite (VIIRS) VNP10A1F, and United States Geological Survey (USGS) Landsat 8 Level-3 Fractional Snow Covered Area. Two spectral-unmixing approaches are also evaluated: Snow-Covered Area and Grain Size (SCAG) and Snow Property Inversion from Remote Sensing (SPIReS), both of which are gap-filled MODIS products and are also run on Landsat 8. We assess subpixel snow mapping performance while considering the fractional snow-covered area (fSCA), canopy cover, sensor zenith angle, and other variables within six global seasonal snow classes. Metrics are calculated at the pixel and basin scales, including the root-mean-square error (RMSE), bias, and F statistic (a detection measure). The newer MOD10A1F Version 61 and VNP10A1F Version 1 product biases ( − 7.1 %, −9.5 %) improve significantly when linear equations developed for older products are applied (2.8 %, −2.7 %) to convert band ratios to fSCA. The F statistics are unchanged (94.4 %, 93.1 %) and the VNP10A1F RMSE improves (18.6 % to 15.7 %), while the MOD10A1F RMSE worsens (12.7 % to 13.7 %). Consistent with previous studies, spectral-unmixing approaches ...