Investigation of the MODIS snow mapping algorithm during snowmelt in the northern boreal forest of Canada (2006)

AbstractA major limitation in mapping snow-covered area (SCA) when using optical satellite sensors is encountered in forested areas, particularly during snowmelt. Consequently, using land cover information to map snow in heterogeneous environments that include dense forest might provide a solution. This study investigated the National Aeronautics and Space Administration (NASA) Goddard Space Flight Center version 4 snow mapping algorithm developed for the moderate resolution imaging spectroradiometer (MODIS) Terra sensor. The functionality of this algorithm for mapping snow during the snowmelt period in the northern boreal forest of Manitoba was tested by comparing the results of this algorithm with those produced by two other algorithms that rely on land cover data instead of the normalized difference vegetation index (NDVI). Comparisons are also made between these three algorithms and the SCA provided by the MOD10A1 MODIS Terra product. For the 2001 snowmelt period, the results from the algorithms and the Terra product showed no noticeable difference in the amount of SCA predicted, but a difference of as much as 20% existed over the 2002 melt period, with the MODIS product predicting more snow than the algorithms. Atmospheric corrections and thresholds used to detect dark targets also had noticeable influences on the amount of SCA predicted by the algorithms. The location of areas mapped with snow also differed between the algorithms and the MOD10A1 product. The agreement in the location of SCA between that predicted by the algorithms and that provided by MOD10A1 differed by as much as 75% towards the end of the 2001 melt period, and even larger discrepancies were observed in 2002. Analysis indicated that the greatest differences occurred in marsh and coniferous areas, and this was consistent throughout the melt period. In addition, higher elevations caused the greatest differences in early stages of the melt period, but lower elevations contributed to the differences in later stages of melt. SCA data for both years showed an exponential correlation with accumulated degree-days for 2001, 2002, and 2003, but the relationship varied depending on the length of the melt period. Snow depletion behavior at different elevation ranges in the study area indicated that snow located at higher elevation ranges melted faster than snow located at lower elevation ranges for 2001, 2002, and 2003. This analysis is based on the MODIS Terra version 4 snow mapping algorithm (as opposed to the algorithm used for the MODIS Aqua sensor), and it is still the most up-to-date version as of early 2006.Cartographier la couverture nivale (CN) à partir d'images satellitaires est particulièrement problématique dans les secteurs couverts de forêts surtout lors de la période de la fonte des neiges. L'utilisation des types de couvertures terrestres pourrait aider à cartographier l'enneigement dans des zones hétérogènes incluant les secteurs forestiers. Ce papier examine l'applicabilité de la version 4 de l'algorithme de cartographie de la couverture nivale du centre NASA (« National Aeronautics and Space Administration ») Goddard Space Flight, version développée pour les données MODIS (« moderate resolution imaging spectroradiometer ») à bord du satellite Terra. Cette version de l'algorithme est testée pour estimer la couverture nivale de la forêt boréale au nord du Manitoba. Ses résultats sont comparés à ceux de deux autres algorithmes basés sur des données de couvertures des terres et non sur l'indice de végétation NDVI, ainsi qu'aux cartes de couverture nivale MOD10A1, produits de Terra-MODIS. Pour la période de fonte des neiges de 2001, les résultats provenant des algorithmes et les cartes MOD10A1 ne montrent pas de différences notables, mais il existe une différence d'au plus 20 % entre l'épaisseur du manteau neigeux estimé par MOD10A1 et les épaisseurs estimés par les trois algorithmes pour la période de fonte des neiges de 2002. Les corrections atmosphériques et les seuils utilisés pour détecter les objets noirs influencent grandement l'estimation de l'épaisseur du manteau neigeux par les algorithmes. De plus, la localisation des zones enneigées diffère entre les cartes MOD10A1 et les cartes produites par les algorithmes : il y a au plus 75 % de différence vers la fin de la période de fonte en 2001 et encore plus de différences sont observées en 2002. L'analyse des résultats montre que les plus grands écarts se situent dans les marais et les forêts de conifères. Ces écarts se retrouvent durant toute la période de la fonte des neiges. Il faut, de plus, noter que les plus grandes différences se situent en hautes altitudes au début de la période et en basses altitudes vers la fin de la période de fonte. Les données CN pour les deux années sont exponentiellement corrélées aux températures de surface pour 2001, 2002 et 2003, mais cette relation dépend de la longueur de la période de fonte. L'analyse montre aussi que le manteau neigeux dans la zone d'étude décroît plus rapidement en hautes altitudes qu'en basses altitudes. Cette analyse est basée sur la version 4 de l'algorithme de cartographie de la couverture nivale pour Terra-MODIS, version toujours la plus récente jusqu'à présent.