Changes in Andes snow cover from MODIS data, 2000–2016

The Andes span a length of 7000 km and are important for sustaining regional water supplies. Snow variability across this region has not been studied in detail due to sparse and unevenly distributed instrumental climate data. We calculated snow persistence (SP) as the fraction of time with snow cover for each year between 2000 and 2016 from Moderate Resolution Imaging Spectroradiometer (MODIS) satellite sensors (500 m, 8-day maximum snow cover extent). This analysis is conducted between 8 and 36∘ S due to high frequency of cloud (> 30 % of the time) south and north of this range. We ran Mann–Kendall and Theil–Sens analyses to identify areas with significant changes in SP and snowline (the line at lower elevation where SP = 20 %). We evaluated how these trends relate to temperature and precipitation from Modern-Era Retrospective Analysis for Research and Applications-2 (MERRA2) and University of Delaware datasets and climate indices as El Niño–Southern Oscillation (ENSO), Southern Annular Mode (SAM), and Pacific Decadal Oscillation (PDO). Areas north of 29∘ S have limited snow cover, and few trends in snow persistence were detected. A large area (34 370 km2) with persistent snow cover between 29 and 36∘ S experienced a significant loss of snow cover (2–5 fewer days of snow year−1). Snow loss was more pronounced (62 % of the area with significant trends) on the east side of the Andes. We also found a significant increase in the elevation of the snowline at 10–30 m year−1 south of 29–30∘ S. Decreasing SP correlates with decreasing precipitation and increasing temperature, and the magnitudes of these correlations vary with latitude and elevation. ENSO climate indices better predicted SP conditions north of 31∘ S, whereas the SAM better predicted SP south of 31∘ S.