Spatiotemporal Variability of Seasonal Snow Cover over 25 Years in the Romanian Carpathians: Insights from a MODIS CGF-Based Approach

Understanding long-term snow cover dynamics is essential in mountain regions with limited meteorological or in situ observations. This study examines seasonal snow cover evolution across the Romanian Carpathians (2000–2025) using daily MODIS/Terra MOD10A1 Cloud-Gap-Filled data at 500 m resolution. Snow-covered pixels were identified using an NDSI ≥ 40 threshold, and snow cover duration (SCD), snow onset date (SOD), and snow end date (SED) were analyzed in relation to elevation and aspect from the FABDEM, complemented by snow-covered area (SCA) and snowline elevation (SLE) metrics. Across the entire range, the snow season shortens mainly due to later onset (+0.28 days/year) and earlier melt (−0.78 days/year), resulting in an SCD decrease of −1.14 days/year. High-elevation (>2000 m) areas show only small changes (SCD: −0.13 days/year; SOD: +0.46 days/year; SED: +0.32 days/year), while the strongest reductions occur at low and mid elevations, where snow persistence is most sensitive to warming; consistent declines in seasonal SCA and a pronounced monthly SLE cycle further document the spatial expression of this variability. Uncertainty was assessed by comparison with station-based snow cover duration (n = 230 station-years), indicating strong agreement (r = 0.95) with a modest negative bias (median: −8 days) and a mean absolute error (MAE) of 16.7 days. Climate correlations highlight air temperature as the dominant covariate of interannual snow-phenology variability, whereas precipitation associations are weaker. Overall, these shifts in snow phenology highlight increasing instability of the Carpathian snow regime and emphasize the value of long-term MODIS observations for tracking cryospheric change in a warming southeastern European mountain system.