Mountain regions have experienced above-average warming in the 20th century and this trend is likely to continue. These accelerated temperature changes in alpine areas are causing reduced snowfall and changes in the timing of snowfall and melt. Snow is a critical component of alpine areas - it drives hibernation of animals, determines the length of the growing season for plants and the soil microbial composition. Thus, changes in snow patterns in mountain areas can have serious ecological consequences. Here we use 35 years of Landsat satellite images to study snow changes in the Mocho-Choshuenco Volcano in the Southern Andes of Chile. Landsat images have 30 m pixel resolution and a revisit period of 16 days. We calculated the total snow area in cloud-free Landsat scenes and the snow frequency per pixel, here called “snow persistence” for different periods and seasons. Permanent snow cover in summer was stable over a period of 30 years and decreased below 20 km2 from 2014 onward at middle elevations (1,530–2,000 m a.s.l.). This is confirmed by negative changes in snow persistence detected at the pixel level, concentrated in this altitudinal belt in summer and also in autumn. In winter and spring, negative changes in snow persistence are concentrated at lower elevations (1,200–1,530 m a.s.l.). Considering the snow persistence of the 1984–1990 period as a reference, the last period (2015–2019) is experiencing a −5.75 km2 reduction of permanent snow area (snow persistence > 95%) in summer, −8.75 km2 in autumn, −42.40 km2 in winter, and −18.23 km2 in spring. While permanent snow at the high elevational belt (>2,000 m a.s.l.) has not changed through the years, snow that used to be permanent in the middle elevational belt has become seasonal. In this study, we use a probabilistic snow persistence approach for identifying areas of snow reduction and potential changes in alpine vegetation. This approach permits a more efficient use of remote sensing data, increasing by three times the amount of usable scenes by including images with spatial gaps. Furthermore, we explore some ecological questions regarding alpine ecosystems that this method may help address in a global warming scenario.