SnowGalileo: A Pre-trained Earth Observation Transformer for Daily, 100 m Fractional Snow Cover Mapping

Abstract. Mountain snow is an important component of the cryosphere that directly affects downstream livelihoods. Accurate monitoring of mountain snow is crucial, yet remains challenging due to complex topography and frequent cloud cover. Although combining data from multiple Earth Observation (EO) satellites can improve spatial and temporal coverage, extracting Fractional Snow Cover (FSC) from sensors with different spatial and temporal resolutions remains difficult. AI-based Earth foundation models can process diverse sensor inputs and have been successfully applied to various remote sensing tasks; however, they often lack snow-specific design considerations. This paper presents SnowGalileo, a pre-trained transformer model designed to integrate EO data to map snow-covered areas. We evaluated SnowGalileo's potential for generating daily, gap-free FSC maps at 100m resolution in mountainous regions. SnowGalileo combines one week of satellite time-series data from multiple sensors, including Sentinel-1, Sentinel-2, Landsat, Sentinel-3, MODIS, and VIIRS, with topographic, land cover, and meteorological information to predict FSC for a given day. The model was pre-trained using masked autoencoding and fine-tuned using labeled data from various mountain ranges across the Northern Hemisphere. In addition to clear-sky conditions, SnowGalileo is evaluated under a wider variety of conditions than was previously possible, including cloud cover and lack of high-resolution (~10–30 m) satellite imagery. For the Canadian Rockies and the Swiss Alps, respectively, SnowGalileo achieves RMSEs of 0.099 and 0.124 on clear days, 0.144 and 0.209 on cloudy days, 0.185 and 0.262 on days without high-resolution imagery, and 0.200 and 0.299 on cloudy days without high-resolution imagery. SnowGalileo also consistently outperforms random forests, support vector regressors, and multi-layer perceptrons. While the current product is a proof of concept that has undergone limited validation across geographic regions, operates on 1km × 1km tiles rather than full maps, and has restricted capabilities in challenging conditions, this approach could ultimately enable the continuous, gap-free operational generation of FSC time series for mountain regions worldwide.