We apply the process-based, distributed TOPKAPI-ETH glacio-hydrological model to a glacierized catchment (19% glacierized) in the semiarid Andes of central Chile. The semiarid Andes provides vital freshwater resources to valleys in Chile and Argentina, but only few glacio-hydrological modelling studies have been conducted, and its dominant hydrological processes remain poorly understood. The catchment contains two debris-free glaciers reaching down to 3900m asl (Bello and Yeso glaciers) and one debris-covered avalanche-fed glacier reaching to 3200m asl (Piramide Glacier). Our main objective is to compare the mass balance and runoff contributions of both glacier types under current climatic conditions. We use a unique dataset of field measurements collected over two ablation seasons combined with the distributed TOPKAPI-ETH model that includes physically oriented parameterizations of snow and ice ablation, gravitational distribution of snow, snow albedo evolution and the ablation of debris-covered ice. Model outputs indicate that while the mass balance of Bello and Yeso glaciers is mostly explained by temperature gradients, the Piramide Glacier mass balance is governed by debris thickness and avalanches and has a clear nonlinear profile with elevation as a result. Despite the thermal insulation effect of the debris cover, the mass balance and contribution to runoff from debris-free and debris-covered glaciers are similar in magnitude, mainly because of elevation differences. However, runoff contributions are distinct in time and seasonality with ice melt starting approximately four weeks earlier from the debris-covered glacier, what is of relevance for water resources management. At the catchment scale, snowmelt is the dominant contributor to runoff during both years. However, during the driest year of our simulations, ice melt contributes 42 +/- 8% and 67 +/- 6% of the annual and summer runoff, respectively. Sensitivity analyses show that runoff is most sensitive to temperature and precipitation gradients, melt factors and debris cover thicknessBecas Chile scholarship program FONDECYT Initiation into Investigation 11130484 FONDECYT 1121184 Chilean Water Authority (Direccion General de Aguas, DGA

No abstract

Studying subglacial drainage networks is important for understanding the potential relationship between channel dynamics and rapid glacier recession as well as the role of subglacial channels in subglacial sediment evacuation. In order to delineate the planform geometry of snout marginal subglacial channels, densely spaced ground-penetrating radar (GPR) measurements at a frequency of ~70 MHz were carried out over the snout marginal zones of two temperate glaciers in the southwestern Swiss Alps, the Haut Glacier d'Arolla and the Glacier d'Otemma. Three-dimensional (3-D) data processing and amplitude analysis of the GPR reflection along the glacier bed was used to map the channels. At the Haut Glacier d'Arolla, two relatively straight channels of several meters in width were identified. The positions of these channels correspond well with the locations of channel outlets at the glacier terminus, as well as with fractures appearing on the glacier surface one month after the GPR data acquisition. The latter are believed to represent the beginning of ice collapse above the subglacial channels. At the Glacier d'Otemma, a major subglacial conduit was detected with similar dimensions to those identified at the Haut Glacier d'Arolla, but greater sinuosity. The position of this channel was confirmed by drone-based imagery acquired after glacier margin collapse. Our results confirm that high-density 3-D GPR surveys can be used to map subglacial channels near temperate alpine glacier margins.

Negative glacier mass balance can also be driven by reduced snow accumulation. Less considered is basal or internal ablation. This can involve the collapse of subglacial channels in the snout marginal zone, driven by thinning ice combined with slow creep closure. After the collapse, ice is removed via the channel to the glacier outlet. This mechanism of glacier retreat was first described some time ago as "subglacial stoping" or "block caving" (Loewe, 1957;Paige, 1956).

<p>Alpine glacier retreat has increased markedly since the late 1980s and is commonly linked to the effects of rising air temperature on surface melt. Less considered are processes associated with glacier snout-marginal surface collapse. A survey of 22 retreating Swiss glaciers suggests that collapse events have increased in frequency since the early 2000s, driven by ice thinning and reductions in glacier-longitudinal ice flux.</p><p>Detailed measurement of a collapse event at one glacier with Uncrewed Aerial Vehicles and ablation stakes showed 0.02 m/day vertical surface deformation above a meandering main subglacial channel, the planform of which was mapped with Ground Penetrating Radar measurements. However, with low rates of longitudinal flux (<1.3 m/year), ice creep was insufficient to close the channel in the snout marginal zone. We hypothesize that an open channel maintains contact between subglacial ice and the atmosphere, allowing greater incursion of warm air up-glacier, thus enhancing melt from below. The associated meandering of subglacial channels at glacier snouts leads to surface collapse due to erosion and internal melt as well as removal of ice via fluvial processes.</p>