Snow Cover Reconstruction in the Brunswick Peninsula, Patagonia, Derived from a Combination of the Spectral Fusion, Mixture Analysis, and Temporal Interpolation of MODIS Data

Several methods based on satellite data products are available to estimate snow cover properties, each one with its pros and cons. This work proposes and implements a novel methodology that integrates three main processes applied to MODIS satellite data for snow cover property reconstruction: (1) the increase in the spatial resolution of MODIS (MOD09) data to 250 m using a spectral fusion technique; (2) a new proposal of snow-cloud discrimination; (3) the daily spatio-temporal reconstruction of snow extent and its albedo signature using the endmembers extraction and spectral mixture analyses. The snow cover reconstruction method was applied to the Brunswick Peninsula, Chilean Patagonia, a low-elevation (<1500 m a.s.l.) mid-latitude area. The results show a 98% agreement between MODIS snow detection and ground-based snow measurements at the automatic weather station, Tres Morros (53.3174°S, 71.2790°W), with fractional snow cover values between 20% and 50%, showing a close relationship between snow and vegetation type. The number of snow days compiled from the MODIS data indicates a good performance (Pearson’s correlation of 0.9) compared with the number of skiing days at the Cerro Mirador ski center, Punta Arenas. Although the number of seasonal snow days showed a significant increasing trend of 0.54 days/year in the Brunswick Peninsula during the 2000–2020 period, a significant decrease of −4.64 days/year was detected in 2010–2020.

[1]  R. Rozzi,et al.  The extreme rainfall gradient of the Cape Horn Biosphere Reserve and its impact on forest bird richness , 2022, Biodiversity and Conservation.

[2]  J. Thepaut,et al.  Supplementary material to "ERA5-Land: A state-of-the-art global reanalysis dataset for land applications" , 2021, Earth System Science Data.

[3]  A. Rivera,et al.  Surface ablation and its drivers along a west–east transect of the Southern Patagonia Icefield , 2021, Journal of Glaciology.

[4]  D. Zheng,et al.  The ERA5-Land soil temperature bias in permafrost regions , 2020 .

[5]  Maogui Hu,et al.  atakrig: An R package for multivariate area-to-area and area-to-point kriging predictions , 2020, Comput. Geosci..

[6]  B. Davies The evolution of the Patagonian Ice Sheet from 35 ka to the present day (PATICE) , 2020, Earth-Science Reviews.

[7]  R. R. de Lima,et al.  Space–Time Kriging of Precipitation: Modeling the Large-Scale Variation with Model GAMLSS , 2019, Water.

[8]  L. Gallardo,et al.  Dynamical downscaling over the complex terrain of southwest South America: present climate conditions and added value analysis , 2019, Climate Dynamics.

[9]  P. Sirguey,et al.  Characterising spatio-temporal variability in seasonal snow cover at a regional scale from MODIS data: the Clutha Catchment, New Zealand , 2019, Hydrology and Earth System Sciences.

[10]  T. Sauter Revisiting extreme precipitation amounts over southern South America and implications for the Patagonian Icefields , 2019, Hydrology and Earth System Sciences.

[11]  Robert M. Graham,et al.  Comparison of ERA5 and ERA-Interim near-surface air temperature, snowfall and precipitation over Arctic sea ice: effects on sea ice thermodynamics and evolution , 2019, The Cryosphere.

[12]  Manfred Mudelsee,et al.  Trend analysis of climate time series: A review of methods , 2019, Earth-Science Reviews.

[13]  G. Casassa,et al.  Snow Cover Change as a Climate Indicator in Brunswick Peninsula, Patagonia , 2018, Front. Earth Sci..

[14]  Y. Hong,et al.  Developing a composite daily snow cover extent record over the Tibetan Plateau from 1981 to 2016 using multisource data , 2018, Remote Sensing of Environment.

[15]  R. Villalba,et al.  Climate and Nothofagus pumilio Establishment at Upper Treelines in the Patagonian Andes , 2018, Frontiers in Earth Science.

[16]  C. Mattar,et al.  Decrease in Snow Cover over the Aysén River Catchment in Patagonia, Chile , 2018 .

[17]  T. Sauter,et al.  A 17-year Record of Meteorological Observations Across the Gran Campo Nevado Ice Cap in Southern Patagonia, Chile, Related to Synoptic Weather Types and Climate Modes , 2018, Front. Earth Sci..

[18]  R. Fensholt,et al.  Snow cover and snow albedo changes in the central Andes of Chile and Argentina from daily MODIS observations (2000–2016) , 2018 .

[19]  Gerard B. M. Heuvelink,et al.  Geographically Weighted Area-to-Point Regression Kriging for Spatial Downscaling in Remote Sensing , 2018, Remote. Sens..

[20]  Pablo A. Mendoza,et al.  The CAMELS-CL dataset: catchment attributes and meteorology for large sample studies – Chile dataset , 2018, Hydrology and Earth System Sciences.

[21]  Michael Lehning,et al.  The European mountain cryosphere: a review of its current state, trends, and future challenges , 2018 .

[22]  G. Liston,et al.  The Andes Cordillera. Part I: snow distribution, properties, and trends (1979–2014) , 2017 .

[23]  Christine Pohl,et al.  Remote Sensing Image Fusion: A Practical Guide , 2016 .

[24]  A. Stehr,et al.  Snow cover dynamics in Andean watersheds of Chile (32.0–39.5° S) during the years 2000–2016 , 2016 .

[25]  D. Easterling,et al.  Detection and attribution of climate extremes in the observed record , 2016 .

[26]  Edward Hanna,et al.  Mass loss and imbalance of glaciers along the Andes Cordillera to the sub-Antarctic islands , 2015 .

[27]  D. Roberts,et al.  A multi-temporal spectral library approach for mapping vegetation species across spatial and temporal phenological gradients , 2015 .

[28]  P. Atkinson,et al.  Downscaling MODIS images with area-to-point regression kriging , 2015 .

[29]  Steven A. Margulis,et al.  Analysis of sub-pixel snow and ice extent over the extratropical Andes using spectral unmixing of historical Landsat imagery , 2014 .

[30]  Martin Charlton,et al.  GWmodel: An R Package for Exploring Spatial Heterogeneity Using Geographically Weighted Models , 2013, 1306.0413.

[31]  Marius Schaefer,et al.  Modeling past and future surface mass balance of the Northern Patagonia Icefield , 2013 .

[32]  Phaedon C. Kyriakidis,et al.  Fusion of MODIS Images Using Kriging With External Drift , 2013, IEEE Transactions on Geoscience and Remote Sensing.

[33]  M. Rojas,et al.  Large-Scale Control on the Patagonian Climate , 2013 .

[34]  Y. Arnaud,et al.  Linking glacier annual mass balance and glacier albedo retrieved from MODIS data , 2012 .

[35]  E. Pebesma spacetime: Spatio-Temporal Data in R , 2012 .

[36]  Feng Gao,et al.  A Data Mining Approach for Sharpening Thermal Satellite Imagery over Land , 2012, Remote. Sens..

[37]  Frank Lamy,et al.  A review of Glacial and Holocene paleoclimate records from southernmost Patagonia (49–55°S) , 2012 .

[38]  Shihyan Lee,et al.  A review of global satellite-derived snow products , 2012 .

[39]  Matthew F. McCabe,et al.  Satellite based observations for seasonal snow cover detection and characterisation in Australia , 2012 .

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

[41]  Pol Coppin,et al.  Endmember variability in Spectral Mixture Analysis: A review , 2011 .

[42]  Yuedong Wang,et al.  Smoothing Splines: Methods and Applications , 2011 .

[43]  V. Markgraf,et al.  Late and postglacial vegetation and fire history in Southern Patagonia and Tierra del Fuego , 2010 .

[44]  Nicholas C. Coops,et al.  Forest canopy effects on snow accumulation and ablation: an integrative review of empirical results. , 2010 .

[45]  Shihyan Lee,et al.  A comparison of optical-band based snow extent products during spring over North America , 2010 .

[46]  Mario Chica-Olmo,et al.  DSCOKRI: A library of computer programs for downscaling cokriging in support of remote sensing applications , 2010, Comput. Geosci..

[47]  B. Luckman,et al.  Spatio‐temporal rainfall patterns in Southern South America , 2009 .

[48]  M. Grosjean,et al.  Long-term multi-proxy climate reconstructions and dynamics in South America (LOTRED-SA): State of the art and perspectives , 2009 .

[49]  Thomas H. Painter,et al.  Retrieval of subpixel snow covered area, grain size, and albedo from MODIS , 2009 .

[50]  A. Zuur,et al.  Mixed Effects Models and Extensions in Ecology with R , 2009 .

[51]  John Turner,et al.  State of the Antarctic and Southern Ocean climate system , 2009 .

[52]  Y. Arnaud,et al.  Subpixel monitoring of the seasonal snow cover with MODIS at 250 m spatial resolution in the Southern Alps of New Zealand: Methodology and accuracy assessment , 2009 .

[53]  Thomas H. Painter,et al.  Time-space continuity of daily maps of fractional snow cover and albedo from MODIS , 2008 .

[54]  N. Molotch,et al.  Effects of vegetation on snow accumulation and ablation in a mid‐latitude sub‐alpine forest , 2008 .

[55]  Richard A. Frey,et al.  Cloud Detection with MODIS. Part I: Improvements in the MODIS Cloud Mask for Collection 5 , 2008 .

[56]  Yves Arnaud,et al.  Snow cover monitoring in the Northern Patagonia Icefield using MODIS satellite images (2000–2006) , 2008 .

[57]  S. Baronti,et al.  Multispectral and panchromatic data fusion assessment without reference , 2008 .

[58]  M. Dubey,et al.  Remote sensing of Greenland ice sheet using multispectral near‐infrared and visible radiances , 2007 .

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

[60]  David A. Seal,et al.  The Shuttle Radar Topography Mission , 2007 .

[61]  M. Canty Image Analysis, Classification, and Change Detection in Remote Sensing , 2006 .

[62]  Matthew Rodell,et al.  Updating a Land Surface Model with MODIS-Derived Snow Cover , 2004 .

[63]  Edzer J. Pebesma,et al.  Multivariable geostatistics in S: the gstat package , 2004, Comput. Geosci..

[64]  Thomas H. Painter,et al.  MULTISPECTRAL AND HYPERSPECTRAL REMOTE SENSING OF ALPINE SNOW PROPERTIES , 2004 .

[65]  G. Marshall Trends in the Southern Annular Mode from Observations and Reanalyses , 2003 .

[66]  L. Hughes Climate change and Australia: Trends, projections and impacts , 2003 .

[67]  Chris Chatfield,et al.  The Analysis of Time Series : An Introduction, Sixth Edition , 2003 .

[68]  Luciano Alparone,et al.  Image fusion—the ARSIS concept and some successful implementation schemes , 2003 .

[69]  G. Casassa,et al.  Weather Observations Across the Southern Andes at 53°S , 2003 .

[70]  Andrea Garzelli,et al.  Context-driven fusion of high spatial and spectral resolution images based on oversampled multiresolution analysis , 2002, IEEE Trans. Geosci. Remote. Sens..

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

[72]  A. Bovik,et al.  A universal image quality index , 2002, IEEE Signal Processing Letters.

[73]  Xavier Otazu,et al.  Multiresolution-based image fusion with additive wavelet decomposition , 1999, IEEE Trans. Geosci. Remote. Sens..

[74]  Thomas H. Painter,et al.  The Effect of Grain Size on Spectral Mixture Analysis of Snow-Covered Area from AVIRIS Data , 1998 .

[75]  Margaret E. Gardner,et al.  Mapping Chaparral in the Santa Monica Mountains Using Multiple Endmember Spectral Mixture Models , 1998 .

[76]  C. Wessman,et al.  Long-term studies of snow-vegetation interactions , 1993 .

[77]  P. Sen Estimates of the Regression Coefficient Based on Kendall's Tau , 1968 .

[78]  H. B. Mann Nonparametric Tests Against Trend , 1945 .

[79]  E. Zsoter,et al.  Author’s Responses to RC1’s comments on “Brief communication: Improving ERA5-Land soil temperature in permafrost regions using an optimized multi-layer snow scheme” , 2022 .

[80]  Kazutoshi Sato,et al.  Comparison of Arctic sea ice thickness and snow depth estimates from CFSR with in situ observations , 2017, Climate Dynamics.

[81]  R. Rozzi,et al.  La reserva de biósfera Cabo de Hornos: un desafío para la conservación de la biodiversidad e implementación del desarrollo sustentable en el extremo austral de América , 2007 .

[82]  P. Aceituno,et al.  RECENT TEMPERATURE VARIATIONS IN SOUTHERN SOUTH AMERICA , 1997 .

[83]  V. K. Shettigara,et al.  A generalized component substitution technique for spatial enhancement of multispectral images using , 1992 .