Three glacier–rock glacier transitional landforms in the central Andes of Chile are investigated over the last decades in order to highlight and question the significance of their landscape and flow dynamics. Historical (1955–2000) aerial photos and contemporary (> 2000) Geoeye satellite images were used together with common processing operations, including imagery orthorectification, digital elevation model generation, and image feature tracking. At each site, the rock glacier morphology area, thermokarst area, elevation changes, and horizontal surface displacements were mapped. The evolution of the landforms over the study period is remarkable, with rapid landscape changes, particularly an expansion of rock glacier morphology areas. Elevation changes were heterogeneous, especially in debris-covered glacier areas with large heaving or lowering up to more than ±1 m yr⁻¹. The use of image feature tracking highlighted spatially coherent flow vector patterns over rock glacier areas and, at two of the three sites, their expansion over the studied period; debris-covered glacier areas are characterized by a lack of movement detection and/or chaotic displacement patterns reflecting thermokarst degradation; mean landform displacement speeds ranged between 0.50 and 1.10 m yr⁻¹ and exhibited a decreasing trend over the studied period. One important highlight of this study is that, especially in persisting cold conditions, rock glaciers can develop upward at the expense of debris-covered glaciers. Two of the studied landforms initially (prior to the study period) developed from an alternation between glacial advances and rock glacier development phases. The other landform is a small debris-covered glacier having evolved into a rock glacier over the last half-century. Based on these results it is proposed that morphological and dynamical interactions between glaciers and permafrost and their resulting hybrid landscapes may enhance the resilience of the mountain cryosphere against climate change.

Streamflow prediction is a key variable for water resources management. It becomes more important in semi-arid regions such as the Tensift river basin in Morocco, where water resources are facing a severe drought and the demand is continuously increasing. The present analysis focuses on evaluating Machine Learning techniques, namely support vector regression (SVR) and Random Forest (RF) against the multiple linear regression (MLR) for daily streamflow forecasting in the mountainous sub-basin of Rheraya between 2003 and 2016. The results show that SVR performed best, followed by RF and MLR. In measurable terms and regarding mean performance, SVR exhibited the higher Nash–Sutcliffe efficiency score (NSE = 0.59) and a lower root mean squared error (RMSE = 1.18 $$\text {m}^3\,\text {s}^{-1}$$m3s-1 ) compared to RF (NSE = 0.53, RMSE = 1.18 $$\text {m}^3\,\text {s}^{-1}$$m3s-1 ) and MLR (NSE = 0.54, RMSE = 1.01 $$\text {m}^3\,\text {s}^{-1}$$m3s-1 ). Furthermore,the available time series was too short to properly capture the full range of streamflow variability, which reduced the prediction performance outside of the calibration conditions. These findings suggest that ML algorithms, particularly SVR, can provide accurate streamflow estimation useful for water resources management when trained on a representative period. The results highlight the capacity of Machine Learning algorithms, specifically SVR, to augment streamflow prediction for enhanced water resource management in arid regions.

In the semiarid Andes of Chile, farmers and industry in the cordillera lowlands depend on water from snowmelt, as annual rainfall is insufficient to meet their needs. Despite the importance of snow cover for water resources in this region, understanding of snow depth distribution and snow mass balance is limited. Whilst the effect of wind on snow cover pattern distribution has been assessed, the relative importance of melt versus sublimation has only been studied at the point scale over one catchment. Analyzing relative ablation rates and evaluating uncertainties are critical for understanding snow depth sensitivity to variations in climate and simulating the evolution of the snowpack over a larger area and over time. Using a distributed snowpack model (SnowModel), this study aims to simulate melt and sublimation rates over the instrumented watershed of La Laguna (513 km², 3150–5630 m a.s.l., 30^∘ S, 70^∘ W), during two hydrologically contrasting years (i.e., dry vs. wet). The model is calibrated and forced with meteorological data from nine Automatic Weather Stations (AWSs) located in the watershed and atmospheric simulation outputs from the Weather Research and Forecasting (WRF) model. Results of simulations indicate first a large uncertainty in sublimation-to-melt ratios depending on the forcing as the WRF data have a cold bias and overestimate precipitation in this region. These input differences cause a doubling of the sublimation-to-melt ratio using WRF forcing inputs compared to AWS. Therefore, the use of WRF model output in such environments must be carefully adjusted so as to reduce errors caused by inherent bias in the model data. For both input datasets, the simulations indicate a similar sublimation fraction for both study years, but ratios of sublimation to melt vary with elevation as melt rates decrease with elevation due to decreasing temperatures. Finally results indicate that snow persistence during the spring period decreases the ratio of sublimation due to higher melt rates.

TanDEM-X digital elevation model (DEM) is a global DEM released by the German Aerospace Center (DLR) at outstanding resolution of 12 m. However, the procedure for its creation involves the combination of several DEMs from acquisitions spread between 2011 and 2014, which casts doubt on its value for precise glaciological change detection studies. In this work we present TanDEM-X DEM as a high-quality product ready for use in glaciological studies. We compare it to Aerial Laser Scanning (ALS)-based dataset from April 2013 (1 m), used as the ground-truth reference, and Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) V003 DEM and SRTM v3 DEM (both 30 m), serving as representations of past glacier states. We use a method of sub-pixel coregistration of DEMs by Nuth and Kääb (2011) to determine the geometric accuracy of the products. In addition, we propose a slope-aspect heatmap-based workflow to remove the errors resulting from radar shadowing over steep terrain. Elevation difference maps obtained by subtraction of DEMs are analyzed to obtain accuracy assessments and glacier mass balance reconstructions. The vertical accuracy (± standard deviation) of TanDEM-X DEM over non-glacierized area is very good at 0.02 ± 3.48 m. Nevertheless, steep areas introduce large errors and their filtering is required for reliable results. The 30 m version of TanDEM-X DEM performs worse than the finer product, but its accuracy, −0.08 ± 7.57 m, is better than that of SRTM and ASTER. The ASTER DEM contains errors, possibly resulting from imperfect DEM creation from stereopairs over uniform ice surface. Universidad Glacier has been losing mass at a rate of −0.44 ± 0.08 m of water equivalent per year between 2000 and 2013. This value is in general agreement with previously reported mass balance estimated with the glaciological method for 2012–2014.

Accurate knowledge of snow depth distributions in forested regions is crucial for applications in hydrology and ecology. In such a context, understanding and assessing the effect of vegetation and topographic conditions on snow depth variability is required. In this study, the spatial distribution of snow depth in two agro-forested sites and one coniferous site in eastern Canada was analyzed for topographic and vegetation effects on snow accumulation. Spatially distributed snow depths were derived by unmanned aerial vehicle light detection and ranging (UAV lidar) surveys conducted in 2019 and 2020. Distinct patterns of snow accumulation and erosion in open areas (fields) versus adjacent forested areas were observed in lidar-derived snow depth maps at all sites. Omnidirectional semi-variogram analysis of snow depths showed the existence of a scale break distance of less than 10 m in the forested area at all three sites, whereas open areas showed comparatively larger scale break distances (i.e., 11–14 m). The effect of vegetation and topographic variables on the spatial variability in snow depths at each site was investigated with random forest models. Results show that the underlying topography and the wind redistribution of snow along forest edges govern the snow depth variability at agro-forested sites, while forest structure variability dominates snow depth variability in the coniferous environment. These results highlight the importance of including and better representing these processes in physically based models for accurate estimates of snowpack dynamics.