A snow cover climatology for the Pyrenees from MODIS snow products

The seasonal snow in the Pyrenees is critical for hydropower production, crop irrigation and tourism in France, Spain and Andorra. Complementary to in situ observations, satellite remote sensing is useful to monitor the effect of climate on the snow dynamics. The MODIS daily snow products (Terra/MOD10A1 and Aqua/MYD10A1) are widely used to generate snow cover climatologies, yet it is preferable to assess their accuracies prior to their use. Here, we use both in situ snow observations and remote sensing data to evaluate the MODIS snow products in the Pyrenees. First, we compare the MODIS products to in situ snow depth (SD) and snow water equivalent (SWE) measurements. We estimate the values of the SWE and SD best detection thresholds to 40 mm water equivalent (w.e.) and 150 mm, respectively, for both MOD10A1 and MYD10A1. κ coefficients are within 0.74 and 0.92 depending on the product and the variable for these thresholds. However, we also find a seasonal trend in the optimal SWE and SD thresholds, reflecting the hysteresis in the relationship between the depth of the snowpack (or SWE) and its extent within a MODIS pixel. Then, a set of Landsat images is used to validate MOD10A1 and MYD10A1 for 157 dates between 2002 and 2010. The resulting accuracies are 97% (κ = 0.85) for MOD10A1 and 96% (κ = 0.81) for MYD10A1, which indicates a good agreement between both data sets. The effect of vegetation on the results is analyzed by filtering the forested areas using a land cover map. As expected, the accuracies decrease over the forests but the agreement remains acceptable (MOD10A1: 96%, κ = 0.77; MYD10A1: 95%, κ = 0.67). We conclude that MODIS snow products have a sufficient accuracy for hydroclimate studies at the scale of the Pyrenees range. Using a gap-filling algorithm we generate a consistent snow cover climatology, which allows us to compute the mean monthly snow cover duration per elevation band and aspect classes. There is snow on the ground at least 50% of the time above 1600 m between December and April. We finally analyze the snow patterns for the atypical winter 2011–2012. Snow cover duration anomalies reveal a deficient snowpack on the Spanish side of the Pyrenees, which seems to have caused a drop in the national hydropower production.

[1]  Y. Heymann,et al.  CORINE Land Cover. Technical Guide , 1994 .

[2]  András Bárdossy,et al.  Cloud removal methodology from MODIS snow cover product , 2009 .

[3]  Jocelyn Chanussot,et al.  Improving MODIS Spatial Resolution for Snow Mapping Using Wavelet Fusion and ARSIS Concept , 2008, IEEE Geoscience and Remote Sensing Letters.

[4]  J. Dozier.  Spectral Signature of Alpine Snow Cover from the Landsat Thematic Mapper , 1989 .

[5]  Dennis P. Lettenmaier,et al.  Evaluation of the snow‐covered area data product from MODIS , 2003 .

[6]  Richard Fernandes,et al.  Validation of VEGETATION, MODIS, and GOES + SSM/I snow‐cover products over Canada based on surface snow depth observations , 2003 .

[7]  S. Sandmeier,et al.  Radiometric corrections of topographically induced effects on Landsat TM data in an alpine environment , 1993 .

[8]  Klaus Seidel,et al.  Toward Snowmelt Runoff Forecast Based on Multisensor Remote-Sensing Informnation , 1987, IEEE Transactions on Geoscience and Remote Sensing.

[9]  J. I. López-Moreno,et al.  Influence of snow accumulation and snowmelt on streamflow in the central Spanish Pyrenees / Influence de l’accumulation et de la fonte de la neige sur les écoulements dans les Pyrénées centrales espagnoles , 2004 .

[10]  A. Klein,et al.  Improving snow cover mapping in forests through the use of a canopy reflectance model , 1998 .

[11]  Thomas H. Painter,et al.  Assessment of methods for mapping snow cover from MODIS , 2011 .

[12]  Mark A. Templin,et al.  Validation of the MODIS snow product and cloud mask using student and NWS cooperative station observations in the Lower Great Lakes Region , 2006 .

[13]  S. Vicente‐Serrano,et al.  Atmospheric circulation influence on the interannual variability of snow pack in the Spanish Pyrenees during the second half of the 20th century , 2007 .

[14]  T. Yao,et al.  Evaluation of a cloud-gap-filled MODIS daily snow cover product over the Pacific Northwest USA , 2011 .

[15]  Emmanuel Paquet,et al.  Retour d'expérience et perspectives d'exploitation des Nivomètres à Rayonnement Cosmique d'EDF , 2006 .

[16]  Albert Rango,et al.  Operational snowmelt runoff forecasting in the Spanish Pyrenees using the snowmelt runoff model , 2002 .

[17]  M. Alvarez,et al.  Research Note: Evidence of Employment Niches in Tourism — An Intra-Regional Approach , 2013 .

[18]  M. Bossard,et al.  CORINE land cover technical guide - Addendum 2000 , 2000 .

[19]  Simon Gascoin,et al.  Introducing Hysteresis in Snow Depletion Curves to Improve the Water Budget of a Land Surface Model in an Alpine Catchment , 2014 .

[20]  Jacob Cohen.  A Coefficient of Agreement for Nominal Scales , 1960 .

[21]  M. Marchamalo,et al.  Flow regime patterns and their controlling factors in the Ebro basin (Spain) , 2010 .

[22]  Dorothy K. Hall,et al.  An approach to using snow areal depletion curves inferred from MODIS and its application to land surface modelling in Alaska , 2005 .

[23]  Dorothy K. Hall,et al.  Development and evaluation of a cloud-gap-filled MODIS daily snow-cover product , 2010 .

[24]  Rui A. P. Perdigão,et al.  Estimation of regional snowline elevation (RSLE) from MODIS images for seasonally snow covered mountain basins , 2014 .

[25]  G. Blöschl,et al.  Validation of MODIS snow cover images over Austria , 2006 .

[26]  Günter Blöschl,et al.  Spatio‐temporal combination of MODIS images – potential for snow cover mapping , 2008 .

[27]  D. Ruelland,et al.  Modeling the current and future capacity of water resources to meet water demands in the Ebro basin , 2013 .

[28]  Albert Rango,et al.  II. Snow hydrology processes and remote sensing , 1993 .

[29]  Jeff Dozier,et al.  Snow Mapping and Classification from Landsat Thematic Mapper Data , 1987, Annals of Glaciology.

[30]  D. Hall,et al.  Accuracy assessment of the MODIS snow products , 2007 .

[31]  P. Houser,et al.  Evaluation of the MODIS snow cover fraction product , 2014 .

[32]  Andrew G. Klein,et al.  Validation of daily MODIS snow cover maps of the Upper Rio Grande River Basin for the 2000–2001 snow year , 2003 .

[33]  B. Markham,et al.  Summary of Current Radiometric Calibration Coefficients for Landsat MSS, TM, ETM+, and EO-1 ALI Sensors , 2009 .

[34]  Leo Breiman,et al.  Classification and Regression Trees , 1984 .

[35]  V. V. Salomonsona,et al.  Estimating fractional snow cover from MODIS using the normalized difference snow index , 2004 .

[36]  Sensibilité aux variations climatiques d'un bassin à forts enjeux - Le cas de la Garonne , 2008 .

[37]  Günter Blöschl,et al.  MODIS snow cover mapping accuracy in a small mountain catchment – comparison between open and forest sites , 2012 .

[38]  D. Hall,et al.  Development of methods for mapping global snow cover using moderate resolution imaging spectroradiometer data , 1995 .

[39]  Albert Rango,et al.  An algorithm to address the MODIS bowtie effect , 2004 .

[40]  N. DiGirolamo,et al.  MODIS snow-cover products , 2002 .

[41]  Gérard Dedieu,et al.  A multi-temporal method for cloud detection, applied to FORMOSAT-2, VENµS, LANDSAT and SENTINEL-2 images , 2010 .

[42]  Christophe Latry,et al.  Validation of an automatic image orthorectification processing , 2004, IGARSS 2004. 2004 IEEE International Geoscience and Remote Sensing Symposium.