Search

A marginal snowpack is a recent concept that refers to snow cover that is relatively shallow and transient but still has important environmental and socioeconomic functions. This class of snow cover influences environments over large areas of the world and is very common in the Mediterranean climates, but also in the transition zone between persistent seasonal and ephemeral snowpacks all over the world. However, a globally-accepted definition to identify these environments has not been reached yet. Most of the modelling and monitoring tools for snow have been developed in cold and alpine areas with persistent snowpacks. The study and modelling of shallower snowpacks in milder environments are complicated by their high temporal variability, frequently patchy nature and typical coexistence with shrubs and forest. However, in recent years new remote sensing tools and more complex modelling approaches and specific parametrizations are facilitating their study. This paper provides a definition for marginal snowpacks, identifying them as those with: (i) seasonal peak depths that only occasionally surpass 1.5 m in favorable deposition zones such as wind drifts, ii) having typical seasonal durations ranging from 60 to 120 days, and (iii) that are predominantly in an isothermal state (on >75% of days) with high densification rates and several cycles of accumulation-ablation during each snow season. We review the main existing research challenges to improve scientific capabilities for their study and we discuss their hydrological and environmental relevance at different spatial scales.

Snow patterns in ice-free areas of Greenland play important roles in ecosystems. Within a changing climate, a comprehensive understanding of the snow responses to climate change is of interest to anticipate forthcoming dynamics in these areas. In this study, we analyze the future snowpack evolution of a polar maritime Arctic location, Qeqertarsuaq (Disko Island, Central-Western Greenland). A physically-based snow model (FSM2) is validated and forced with CMIP6 projections for SSP2–4.5 and SSP5–8.5 greenhouse gasses emission scenarios, using two models: CanESM5 and MIROC6. The future snowpack evolution is assessed through four key seasonal (October to May) snow climate indicators: snow depth, snow days, snowfall fraction and ablation rate. Comparison against the observed air temperature for the reference climate period demonstrates superior accuracies for MIROC6 SSP2.4–5, with anomalies at 19 %, compared to CanESM5 SSP5.8–5 (25 %) and CanESM5 SSP2.4–5 (78 %). In terms of precipitation, CanESM5 SSP2.4–5 and SSP2.4–5 exhibit smaller anomalies against the observed data (5 %) in contrast to MIROC6 SSP2.4–5 (15 %) and MIROC6 SSP2.8–5 (17 %). Results demonstrate distinct snowpack responses to climate change depending on the model and emission scenario. For CanESM5, seasonal snow depth anomalies with respect to the reference period range from – 38 % (SSP2–4.5, 2040–2050 period) to – 74 % (SSP5–8.5, 2090–2100 period). MIROC6 projects lower snowpack reductions, with a decrease ranging from – 38 % (SSP2–4.5, 2040–2050 period) to – 57 % (SSP5–8.5, 2090–2100 period). Similar reductions are anticipated for snowfall and snow days. Changes in the snowpack evolution are primarily driven by positive trends in downwelling longwave radiation and air temperature. The projected increase in precipitation by the mid to late 21st century will lead to more frequent rain-on-snow events, intensifying snowpack melting. These findings help enhance the comprehension of future snow dynamics in the

The Aneto glacier, although it may be considered a very small glacier (<0.5 km2), is the largest glacier in the Pyrenees. Its surface and thickness loss have been continuous in recent decades, and there have been signs of accelerated melting in recent years. In this study, thickness and surface losses of the Aneto glacier from 1981 to 2022 are investigated using historical aerial imagery, airborne lidar point clouds and unoccupied aerial vehicle (UAV) imagery. A ground-penetrating radar (GPR) survey conducted in 2020, combined with data from photogrammetric analyses, allowed us to reconstruct the current ice thickness and also the existing ice distribution in 1981 and 2011. Over the last 41 years, the total glacierised area has decreased by 64.7 %, and the ice thickness has decreased, on average, by 30.5 m. The mean remaining ice thickness in autumn 2022 was 11.9 m, as against the mean thickness of 32.9, 19.2 and 15.0 m reconstructed for 1981 and 2011 and observed in 2020, respectively. The results demonstrate the critical situation of the glacier, with an imminent segmentation into two smaller ice bodies and no evidence of an accumulation zone. We also found that the occurrence of an extremely hot and dry year, as observed in the 2021–2022 season, leads to a drastic degradation of the glacier, posing a high risk to the persistence of the Aneto glacier, a situation that could extend to the rest of the Pyrenean glaciers in a relatively short time.

This dataset includes very high spatial resolution observations at 1 m spatial resolution observations of the snow depth, the NDVI and the melt-out date (DOY of year) acquired with an Unmanned Aerial Vehicle at a sub-alpine site in the Pyrenees, the Izas Experimental Catchment. During two snow seasons (2019-2020 and 2020-2021), 14 NDVI and 17 snow depth distributions were acquired over 48ha. From the snow depth observations the melt-out dates have been derived. Also information on the main topographic variables (elevation, aspect and slope) is included, with same spatial resolution, in this dataset.

Accurate knowledge of the seasonal snow distribution is vital in several domains including ecology, water resources management, and tourism. Current spaceborne sensors provide a useful but incomplete description of the snowpack. Many studies suggest that the assimilation of remotely sensed products in physically based snowpack models is a promising path forward to estimate the spatial distribution of snow water equivalent (SWE). However, to date there is no standalone, open-source, community-driven project dedicated to snow data assimilation, which makes it difficult to compare existing algorithms and fragments development efforts. Here we introduce a new data assimilation toolbox, the Multiple Snow Data Assimilation System (MuSA), to help fill this gap. MuSA was developed to fuse remotely sensed information that is available at different timescales with the energy and mass balance Flexible Snow Model (FSM2). MuSA was designed to be user-friendly and scalable. It enables assimilation of different state variables such as the snow depth, SWE, snow surface temperature, binary or fractional snow-covered area, and snow albedo and could be easily upgraded to assimilate other variables such as liquid water content or snow density in the future. MuSA allows the joint assimilation of an arbitrary number of these variables, through the generation of an ensemble of FSM2 simulations. The characteristics of the ensemble (i.e., the number of particles and their prior covariance) may be controlled by the user, and it is generated by perturbing the meteorological forcing of FSM2. The observational variables may be assimilated using different algorithms including particle filters and smoothers as well as ensemble Kalman filters and smoothers along with their iterative variants. We demonstrate the wide capabilities of MuSA through two snow data assimilation experiments. First, 5 m resolution snow depth maps derived from drone surveys are assimilated in a distributed fashion in the Iz

This study extends knowledge of the evolution of glacier shrinkage in the Cocuy-Güican Mountains since the maximum glacier extent of the Little Ice Age (LIA). Mass balance data for the Ritacuba Glacier since 2009 were acquired and compared with available data for the Conejeras Glacier (Los Nevados National Park). This study also investigated the hydrological significance of Colombian glaciers, which is still largely unknown because the available information is very limited. Glaciers in the Cocuy-Güican Mountains covered 13.2 km2 in 2019 compared with 127.8 km2 during the maximum LIA, representing a shrinkage of 89.7%. Analysis of glacier cover observations made in 1955, 1994, 2010, and 2019 revealed that the rate of ice loss was greatest from 1994 to 2010 (0.59 km2 yr−1), and was then almost halved from 2010 to 2019 (0.34 km2 yr−1). This slowing of glacier retreat is consistent with the moderate negative mass balance measured for 2009–2019, and an accumulated loss of 1766 mm w.e. (mm water equivalent). The progressive confinement of glaciers to higher elevations, ice accumulation in topographic locations providing shelter from solar radiation, and an absence of recent marked climatic anomalies could explain why the Cocuy-Güican glaciers have temporally reached a near equilibrium state conditions. This contrasts with the Conejeras Glacier, where 47,000 mm w.e. has been lost in the same period. The available data on runoff and isotope tracers of streamflow and precipitation suggest that precipitation rather than glacier melt water exerts primary control over the hydrological variability at high elevation sites.

In mountain areas, the phenology and productivity of grassland are closely related to snow dynamics. However, the influence that snow melt timing has on grassland growing still needs further attention for a full understanding, particularly at high spatial resolution. Aiming to reduce this knowledge gap, this work exploits 1 m resolution snow depth and Normalized Difference Vegetation Index observations acquired with an Unmanned Aerial Vehicle at a sub-alpine site in the Pyrenees. During two snow seasons (2019–2020 and 2020–2021), 14 NDVI and 17 snow depth distributions were acquired over 48 ha. Despite the snow dynamics being different in the two seasons, the response of grasslands greening to snow melt-out exhibited a very similar pattern in both. The NDVI temporal evolution in areas with distinct melt-out dates reveals that sectors where the melt-out date occurs in late April or early May (optimum melt-out) reach the maximum vegetation productivity. Zones with an earlier or a later melt-out rarely reach peak NDVI values. The results obtained in this study area, suggest that knowledge about snow depth distribution is not needed to understand NDVI grassland dynamics. The analysis did not reveal a clear link between the spatial variability in snow duration and the diversity and richness of grassland communities within the study area.

The duration of the seasonal snowpack determines numerous aspects of the water cycle, ecology and the economy in cold and mountainous regions, and is a balance between the magnitude of accumulated snow and the rate of melt. The contribution of each component has not been well quantified under contrasting topography and climatological conditions although this may provide useful insights into how snow cover duration could respond to climate change. Here, we examined the contribution of the annual peak snow water equivalent (SWE) and the seasonal melt rate to define the duration of the snowpack over temperate mountains, using snow data for mountain areas with different climatological characteristics across the Iberian Peninsula. We used a daily snowpack database for the period 1980–2014 over Iberia to derive the seasonal peak SWE, melt rate and season snow cover duration. The influence of peak SWE and melt rates on seasonal snow cover duration was estimated using a stepwise linear model approach. The stepwise linear models showed high R-adjusted values (average R-adjusted = 0.7), without any clear dependence on the elevation or geographical location. In general, the peak SWE influenced the snow cover duration over all of the mountain areas analysed to a greater extent than the melt rates (89.1%, 89.2%, 81.6% 93.2% and 95.5% in the areas for the Cantabrian, Central, Iberian, Pyrenees and Sierra Nevada mountain ranges, respectively). At these colder sites, the melt season occurs mostly in the spring and tends to occur very fast. In contrast, the areas where the melt rates dominated snow cover duration were located systematically at lower elevations, due to the high interannual variability in the occurrence of annual peak SWE (in winter or early spring), yielding highly variable melt rates. However, in colder sites the melt season occurs mostly in spring and it is very fast in most of the years. The results highlight the control that the seasonal precipitation patterns, i

We performed high-resolution snow height measurements using drone-based lidar over two 0.5 km long areas near the community of Umiujaq, in Tasiapik valley. The North site is centered around (N:56.5684°; W:76.4895°, elevation 122 m) and the South site around (N:56.5594°; W:76.4816°, elevation 133 m). A snow-free survey was performed on 28-29 September 2017 and a snow survey on 25-26 April 2018. We produced 3 cm digital terrain and digital surface models (DTM and DSM) from the lidar cloud points for each site for both lidar campaigns.

The Aneto Glacier, is the largest glacier in the Pyrenees. Its shrinkage and wastage have been continuous in recent decades, and there are signs of accelerated melting in recent years. In this study, changes in the surface and ice thickness of the Aneto Glacier from 1981 to 2022 are investigated using historical aerial imagery, airborne LiDAR point clouds, and UAV imagery. A GPR survey conducted in 2020, combined with data from photogrammetric analyses, allowed us to reconstruct the current ice thickness and also the existing ice distribution in 1981 and 2011. Over the last 41 years, the total glaciated area has shrunk by 64.7% and the ice thickness has decreased, on average, by 30.5 m. The mean remaining ice thickness in autumn 2022 was 11.9 m, as against the mean thicknesses of 32.9 m, 19.2 m reconstructed for 1981 and 2011 and 15.0 m observed in 2020 respectively. The results demonstrate the critical situation of the glacier, with an imminent segmentation into two smaller ice bodies and no evidence of an accumulation zone. We also found that the occurrence of an extremely hot and dry year, as observed in the 2021–2022 season, leads to a drastic degradation of the glacier, posing a high risk to the persistence of the Aneto Glacier, a situation that could extend to the rest of the Pyrenean glaciers in a relatively short time.

Study region: Peñalara catchment, in the mountains of the Central System in Spain, Study focus: For the first time, we investigated the streamflow and streamwater isotopes during the snow accumulation and melting periods and over subsequent months in two snow seasons. The aim is to better understand the hydrological processes linked to snowmelt; to describe the temporal evolution and the interannual differences in isotopic streamwater; and to improve the understanding about the hydrological functioning of snowmelt water across the catchment. New hydrological insights for the region: The isotopic signal of the streamwater progressively became isotopically depleted from the beginning of the melt period until the snow cover depletion of the catchment. Higher snowfall led to depleted isotopic values in the stream compared to a year with low snowfall. The interannual variability of the isotopic signature of streamwater during snowmelt may represent a difficulty to establishing reference values to be used in mixing models for hydrograph separation. The streamflow isotopic values had very limited sub-daily variation and showed slow temporal changes, suggesting a central role for alpine aquifers in explaining the hydrological functioning of the catchment, pointing to piston flow as a key process in streamflow generation.

Expanding shrubs in the Arctic trap blowing snow, increasing snow height and accelerating permafrost warming. Topography also affects snow height as snow accumulates in hollows. The respective roles of topography and erect vegetation in snow accumulation were investigated using a UAV-borne lidar at two nearby contrasted sites in northern Quebec, Canada. The North site featured tall vegetation up to 2.5 m high, moderate snow height, and smooth topography. The South site featured lower vegetation, greater snow height, and rougher topography. There was little correlation between topography and vegetation height at both sites. Vegetation lower than snow height had very little effect on snow height. When vegetation protruded above the snow, snow height was well correlated with vegetation height. The topographic position index (TPI) was well correlated with snow height when it was not masked by the effect of protruding vegetation. The North site with taller vegetation therefore showed a good correlation between vegetation height and snow height, R2 = 0.37, versus R2 = 0.04 at the South site. Regarding topography, the reverse was observed between TPI and snow height, with R2 = 0.29 at the North site and R2 = 0.67 at the South site. The combination of vegetation height and TPI improved the prediction of snow height at the North site (R2 = 0.59) but not at the South site because vegetation height has little influence there. Vegetation was therefore the main factor determining snow height when it protruded above the snow. When it did not protrude, snow height was mostly determined by topography., [Significance Statement] Wind-induced snow drifting is a major snow redistribution process in the Arctic. Shrubs trap drifting snow, and drifting snow accumulates in hollows. Determining the respective roles of both these processes in snow accumulation is required to predict permafrost temperature and its emission of greenhouse gases, because thicker snow limits permafrost winter cooling. Using a UAV-borne lidar, we have determined snow height distribution over two contrasted sites in the Canadian low Arctic, with varied vegetation height and topography. When snow height exceeds vegetation height, topography is a good predictor of snow height, with negligible effect of buried vegetation. When vegetation protrudes above the snow, combining both topography and vegetation height is required for a good prediction of snow height.