Borehole drainage and its implications for the investigation of glacier hydrology: experiences from Haut Glacier d'Arolla,Switzerland

Studies of glacier hydrology rely increasingly on measurements made in boreholes as a basis for reconstructing the character and behaviour of subglacial drainage systems. In temperate glaciers, in which boreholes remain open to the atmosphere following drilling, the interpretation of such data may be complicated by supraglacial or englacial water flows to and from boreholes. We report on a suite of techniques used to identify borehole water sources and to reconstruct patterns of water circulation within boreholes at Haut Glacier d'Arolla, Switzerland. Results are used to define a number of borehole ‘drainage’ types. Examples of each drainage type are presented, along with the manner in which they influence interpretations of borehole water‐levels, borehole water‐quality data, and borehole dye traces. The analysis indicates that a full understanding of possible borehole drainage modes is required for the correct interpretation of many borehole observations, and that those observations provide an accurate indication of subglacial conditions only under relatively restricted circumstances. Copyright © 2001 John Wiley & Sons, Ltd.     

Focal mechanisms of intraplate earthquakes in Bolivia,South America

Intraplate seismic activity in Bolivia is mainly located in the central region (16°–19°S, 63°–67°W) which includes the East Andean Cordillera and the Sub-Andean Sierras. At this region there is a bend in the trend of the main geological structures from NW-SE in the north to N-S in the south. Focal mechanisms have been calculated for 10 earthquakes of magnitudes 4.9–5.6, using first motionP-waves from long period instruments. Their solutions correspond to reverse faulting, some with a large component of strike-slip motion. Their solutions can be grouped into two types; one with pure reverse faulting on planes with azimuth NW-SE and the other with a large strike-slip component on planes with azimuths nearly N-S or WNW-ESE. The maximum stress axis (P-axis) is practically horizontal (dipping less than 5°) oriented in a mean N56°E direction. This orientation may be related with the direction of compression resulting from the collision of the Nazca plate against the western margin of the South American continent. Wave-form analysis of long-periodP-waves for one event restricts the focal depth to 8 km in the Sub-Andean region. Seismic moments and source dimensions determined from spectra of Rayleigh waves are in the range of 10¹⁶–10¹⁷Nm and 17–24 km, respectively. The Central Bolivia region can be considered as a zone of intraplate deformation situated between the Bolivian Altiplano and the Brazil shield.     

Drivers of peatland water table dynamics in the central Andes,Bolivia and Peru

Cushion plant dominated peatlands are key ecosystems in tropical alpine regions of the Andes in South America. The cushion plants have formed peat bodies over thousands of years that fill many valley bottoms, and the forage produced by the plants is critical for native and nonnative domesticated mammals. The sources and flow paths of water supporting these peatlands remain largely unknown. Some studies have suggested that glacier meltwater streams support some peatlands, and that the ongoing loss of glaciers and their meltwaters could lead to the loss or diminishment of peatlands. We analysed the hydrologic regime of 10 peatlands in four mountain regions of Bolivia and Peru using groundwater monitoring. Groundwater levels in peatlands were relatively stable and within 20 cm of the ground surface during the rainy season, and many sites had water tables 40–90 cm below the ground surface in the dry season. Topographic and groundwater elevations in the peatlands demonstrated that the water source of all 10 peatlands was hillslope groundwater flowing from lateral moraines, talus, colluvium, or bedrock aquifers into the peatlands. There was little to no input from streams, whether derived from glacier melt or other sources, and glacier melt could not have recharged the hillslope aquifers supporting peatlands. We measured the stable water isotopes in water samples taken during different seasons, distributed throughout the catchments, and the values are consistent with this interpretation. Our findings indicate that peatlands in the study region are recharged by hillslope groundwater discharge rather than stream water and may not be as vulnerable to glacial decline as other studies have indicated. However, both glaciers and peatlands are susceptible to changing thermal and precipitation regimes that could affect the persistence of peatlands.     

Impact of a rock avalanche on a moraine‐dammed proglacial lake: Laguna Safuna Alta,Cordillera Blanca,Peru

Moraines that dam proglacial lakes pose an increasing hazard to communities in the Andes and other mountain ranges. The moraines are prone to failure through collapse, overtopping by lake waters or the effect of displacement waves resulting from ice and rock avalanches. Resulting floods have led to the loss of thousands of lives in the Cordillera Blanca mountains of Peru alone in the last 100 years. On 22 April 2002 a rock avalanche occurred immediately to the south‐west of Laguna Safuna Alta, in the Cordillera Blanca. The geomorphic evidence for the nature, magnitude and consequences of this event was investigated in August 2002. Field mapping indicated that the avalanche deposited 8–20 × 10⁶ m³ of rock into the lake and onto the surface of the frontal region of Glaciar Pucajirca, which flows into the lake. Repeated bathymetric surveying indicated that ～5 × 10⁶ m³ of this material was deposited directly into the lake. The immediate effect of this event was to create a displacement wave that gained in height as it travelled along the lake basin, overtopping the impounding moraine at the lake's northern end. To achieve overtopping, the maximum wave height must have been greater than 100 m. This, and subsequent seiche waves, caused extensive erosion of both the proximal and distal faces of the impounding terminal moraine. Further deep gullying of the distal face of this moraine resulted from the supply of pressurized water to the face via a relief overflow tunnel constructed in 1978. Two‐dimensional, steady‐state analysis of the stability of the post‐avalanche moraine rampart indicates that its proximal face remains susceptible to major large‐scale rotational failure. Copyright © 2005 John Wiley & Sons, Ltd.     

Holocene calcrete crust deposits on the moraine of Batura Glacier, northern Pakistan

Abstract   Calcretes can be observed on the surface of old moraines around Batura Glacier in the upper Hunza Valley, Karakoram Mountains, Pakistan. They develop as a calcareous crust cementing small gravels under boulders. In order to understand the genesis of the calcrete crust, a variety of methods were employed: (i) study of mineralogy and geochemistry of a calcrete crust precipitated on the lateral moraine using X-ray diffractometer and electron probe microanalysis; (ii) analysis of solute chemistry of surface water and ice bodies around the Batura Glacier; and (iii) accelerator mass spectrometry ¹⁴C dating of the crust itself. The results indicate that the calcrete crust has definite laminated layers composed of a fine-grain and compact calcite layer, and a mineral fragment layer. The chemical composition of the calcite layer is approximately 60% CaO and 1% MgO. The mineral fragment layer consists of rounded grain materials up to 0.2 mm in diameter. It shows a graded bedding structure with fine grains of quartz, albite and muscovite. Meanwhile, as the Paleozoic Pasu limestone is distributed around the terminal of Batura Glacier, Ca cations dissolve in the melt water of the glacier. Accordingly, the calcrete crust is precipitated by decreases in CO₂ partial pressure from glacier ice and evaporation of the melt water, including high concentration of Ca²⁺ at ephemeral streams and small ponds stagnating between the moraine and glacial ice. On the basis of the AMS ¹⁴C age, the calcrete is considered to have formed approximately 8200 calibrated years bp under the Batura glacial stage.     

Water and sediment outbursts from advanced Franz Josef Glacier,New Zealand

The Franz Josef Glacier, Westland, New Zealand, has a history of catastrophic sediment‐laden outburst ?oods associated with extreme rainfall events when the glacier toe is advanced over its own sediments. Consideration of these events and inspection of recent sediment deposits suggest that there are three distinct modes of outburst. The ?rst is associated with fans fed by over?ow along the glacier margin. As the glacier has advanced across its own fore?eld gravels, it is inferred that the primary drainage conduit has developed a reach of negative slope. In high ?ows massive boulders can block the conduit, trapping lesser clasts. The resulting backup of water causes over?ows through marginal moulins, producing the fan type of deposit. The second type of outburst deposits massive imbricated boulders at a greater or lesser distance from the glacier portal. In this case, pressure buildup drives the blockage out of the portal where the boulders deposit. Smaller materials are generally carried away. The third type consists of very shallow ?ows, and produces massive gravel deposits of uncertain provenance. In this condition, the excess pressure in the conduit results in slight uplift of the glacier and widespread discharge of water and sediment below the glacier snout; gravels and smaller sediments are laid down in a massive deposit across the fore?eld. The massive, boulder‐veneered deposit from the December 1995 outburst is interpreted in the light of the above mechanisms as a hyperconcentrated ?ow deposit from hydraulic jacking, overlain by boulders emplaced by a subsequent conduit outburst. A possible association of outbursts with the present advanced position of the glacier is suggested. Copyright © 2003 John Wiley & Sons, Ltd.     

Revised palaeoclimatic significance of Mueller Glacier moraines,Southern Alps,New Zealand

A sequence of Late Holocene moraines on the foreland of the Mueller Glacier, Southern Alps, New Zealand, forms part of a local moraine‐age database used to establish a regional glacier chronology and subsequently to investigate potential intra‐hemispheric and global climate forcing mechanisms. We present new sedimentological and geomorphological evidence that a set of these moraine ridges, previously considered to represent individual advances, constitutes a single moraine complex (the ‘Mueller Memorial Moraine’) formed by supraglacial transport of a large volume of landslide debris to the glacier terminus. Because a moraine formed in this way is not necessarily associated with an advance triggered by a climate event, we question the palaeoclimatic significance of the Mueller Memorial Moraine, as well as that of other moraines in comparable settings. Our findings suggest that the mode of formation and glacio‐dynamical context of moraines whose ages contribute to existing palaeoclimate reconstructions need to be re‐examined in order to assess the reliability of these reconstructions. Copyright © 2015 John Wiley & Sons, Ltd.     

Heterogeneous water storage and thermal regime of supraglacial ponds on debris‐covered glaciers

The water storage and energy transfer roles of supraglacial ponds are poorly constrained, yet they are thought to be important components of debris‐covered glacier ablation budgets. We used an unmanned surface vessel (USV) to collect sonar depth measurements for 24 ponds to derive the first empirical relationship between their area and volume applicable to the size distribution of ponds commonly encountered on debris‐covered glaciers. Additionally, we instrumented nine ponds with thermistors and three with pressure transducers, characterizing their thermal regime and capturing three pond drainage events. The deepest and most irregularly‐shaped ponds were those associated with ice cliffs, which were connected to the surface or englacial hydrology network (maximum depth = 45.6 m), whereas hydrologically‐isolated ponds without ice cliffs were both more circular and shallower (maximum depth = 9.9 m). The englacial drainage of three ponds had the potential to melt ~100 ± 20 × 10³ kg to ~470 ± 90 × 10³ kg of glacier ice owing to the large volumes of stored water. Our observations of seasonal pond growth and drainage with their associated calculations of stored thermal energy have implications for glacier ice flow, the progressive enlargement and sudden collapse of englacial conduits, and the location of glacier ablation hot‐spots where ponds and ice cliffs interact. Additionally, the evolutionary trajectory of these ponds controls large proglacial lake formation in deglaciating environments. Copyright © 2017 John Wiley & Sons, Ltd.     

Impact of post-mixing chemical reactions on the major ion chemistry of bulk meltwaters draining the haut glacier d'arolla,valais, Switzerland

Until now, alpine glacial meltwaters have been assumed to consist of two components, dilute quickflow and concentrated delayed flow, the mixing of which has been regarded as chemically conservative for the major dissolved ions and electrical conductivity. Dye tracing results suggest that this two-component model adequately represents the sub-glacial hydrology of the Haut Glacier d'Arolla, Switzerland. However, laboratory dissolution experiments in which various concentrations of glacial rock flour are placed in dilute solutions show that this rock flour is highly reactive and suggest that bulk meltwaters may acquire significant amounts of solute through rapid chemical reactions with suspended sediment which occur after mixing of the two components. This view is supported by detailed analysis of variations in the hydrochemistry of meltwaters draining from the Haut Glacier d'Arolla over three diurnal cycles during the 1989 melt season. Variations in the composition of bulk meltwaters are controlled by two main factors: dilution of the delayed flow component by quickflow, and the extent of post-mixing reactions. The latter depends on the suspended sediment concentration in bulk meltwaters and on the duration of contact between these waters and suspended sediment. Seasonal changes in the magnitude of these factors result in changes in the character and causes of diurnal variations in meltwater chemistry. In June, these variations reflect discharge-related variations in residence time within a distributed subglacial drainage system; in July, when a channelized drainage system exists beneath the lower glacier, they primarily reflect the dilution of delayed flow by quickflow; in August, when suspended sediment concentrations are particularly high, they reflect varying degrees of solute acquisition by post-mixing reactions with suspended sediment that take place in arterial channels at the glacier bed.     

A conceptual model of supra‐glacial lake formation on debris‐covered glaciers based on GPR facies analysis

Supra‐glacial lakes and ponds can create hotspots of mass loss on debris‐covered glaciers. While much research has been directed at understanding lateral lake expansion, little is known about the rates or processes governing lake deepening. To a large degree, this knowledge gap persists due to sparse observations of lake beds. Here we report on the novel use of ground penetrating radar (GPR) surveys to simultaneously collect supra‐glacial lake bathymetry and bottom composition data from Spillway Lake (surface area of 2.4 × 10⁵ m²; volume of 9.5 × 10⁴ m³), which is located in the terminus region of the Ngozumpa Glacier in the Khumbu region of the Nepal Himalaya. We identified two GPR bottom signals corresponding to two sedimentary facies of (1) sub‐horizontal layered fine sediment drape and (2) coarse blocky diamict. We provide an understanding of the changes in subaqueous debris distribution that occur through stages of lake expansion by combining the GPR results with in situ observations of shoreline deposits matching the interpreted facies. From this, we present an updated conceptual model of supra‐glacial lake evolution, with the addition of data on the evolving debris environment, showing how dominant depositional processes can change as lakes evolve from perched lakes to multi‐basin base‐level lakes and finally onto large moraine‐dammed lakes. Throughout lake evolution, processes such as shoreline steepening, lakebed collapse into voids and conduit interception, subaerial and subaqueous calving and rapid areal expansion alter the spatial distribution and makeup of lakebed debris and sediments forcing a number of positive and negative feedbacks on lake expansion. Copyright © 2016 John Wiley & Sons, Ltd.     

Glacial conditions that contribute to the regeneration of Fountain Glacier proglacial icing,Bylot Island,Canada

Proglacial icings are one of the most common forms of extrusive ice found in the Canadian Arctic. However, the icing adjacent to Fountain Glacier, Bylot Island, is unique due to its annual cycle of growth and decay, and perennial existence without involving freezing point depression of water due to chemical characteristics. Its regeneration depends on the availability of subglacial water and on the balance between ice accretion and hydro‐thermal erosion. The storage and conduction of the glacial meltwater involved in the accretion of the icing were analyzed by conducting topographic and ground penetrating radar surveys in addition to the modelling of the subglacial drainage network and the thermal characteristics of the glacier base. The reflection power analysis of the geophysical data shows that some areas of the lower ablation zone have a high accumulation of liquid water, particularly beneath the centre part of the glacier along the main supraglacial stream. A dielectric permittivity model of the glacier – sediment interface suggests that a considerable portion of the glacier is warm based; allowing water to flow through unfrozen subglacial sediments towards the proglacial outwash plain. All these glacier‐related characteristics contribute to the annual regeneration of the proglacial icing and allow for portions of the icing to be perennial. Copyright © 2013 John Wiley & Sons, Ltd.     

Athabasca Glacier,Canada – a field example of subglacial ice and till erosion?

The Athabasca Glacier, resting on a rigid bed, provides an excellent example of subglacial ice and till erosion. The presence of a thin mobile till layer is shown by the presence of flutes, saturated till layer, push moraines and ploughed boulders. Cross‐cutting striations, v‐shaped striations and reversed stoss‐and‐lee clasts are indicative of clasts rotating within this layer. As the till moves it erodes the bedrock and clasts within it. A combination of erosion by ice and till produces stoss‐and‐lee‐clasts and generates striations on flutes and embedded clasts, as well as eroding the bedrock into a continuum of smoothed, rounded and streamlined forms. Copyright © 2005 John Wiley & Sons, Ltd.     

Hydro‐meteorological drivers and sources of suspended sediment flux in the pro‐glacial zone of the retreating Castle Creek Glacier,Cariboo Mountains,British Columbia,Canada

Glaciers are major agents of erosion that increase sediment load to the downstream fluvial system. The Castle Creek Glacier, British Columbia, Canada, has retreated ~1.0 km in the past 70 years. Suspended sediment concentration (SSC) and streamflow (Q) were monitored independently at five sites within its pro‐glacial zone over a 60 day period from July to September 2011, representing part of the ablation season. Meteorological data were collected from two automatic weather stations proximal to the glacier. The time‐series were divided into hydrologic days and the shape and magnitude of the SSC response to hydro‐meteorological conditions (‘cold and wet’, ‘hot and dry’, ‘warm and damp’, and ‘storm’) were categorized using principal component analysis (PCA) and cluster analysis (CA). Suspended sediment load (SSL) was computed and summarized for the categories. The distribution of monitoring sites and results of the multivariate statistical analyses describe the temporal and spatial variability of suspended sediment flux and the relative importance of glacial and para‐glacial sediment sources in the pro‐glacial zone. During the 2011 study period, ~ 60% of the total SSL was derived from the glacial stream and sediment deposits proximal to the terminus of the glacier; during ‘storm’ events, that contribution dropped to ~40% as the contribution from diffuse and point sources of sediment throughout the pro‐glacial zone and within the meltwater channels increased. While ‘storm’ events accounted for just 3% of the study period, SSL was ~600% higher than the average over the monitoring period, and ~20% of the total SSL was generated in that time. Determining how hydro‐meteorological conditions and sediment sources control sediment fluxes will assist attempts to predict how pro‐glacial zones respond to future climate changes. Copyright © 2015 John Wiley & Sons, Ltd.     

Basal sediment evacuation by subglacial meltwater: suspended sediment transport from Haut Glacier d'Arolla,Switzerland

Proglacial suspended sediment transport was monitored at Haut Glacier d'Arolla, Switzerland, during the 1998 melt season to investigate the mechanisms of basal sediment evacuation by subglacial meltwater. Sub‐seasonal changes in relationships between suspended sediment transport and discharge demonstrate that the structure and hydraulics of the subglacial drainage system critically influenced how basal sediment was accessed and entrained. Under hydraulically inefficient subglacial drainage at the start of the melt season, sediment availability was generally high but sediment transport increased relatively slowly with discharge. Later in the melt season, sediment transport increased more rapidly with discharge as subglacial meltwater became confined to a spatially limited network of channels following removal of the seasonal snowpack from the ablation area. Flow capacity is inferred to have increased more rapidly with discharge within subglacial channels because rapid changes in discharge during highly peaked diurnal runoff cycles are likely to have been accommodated largely by changes in flow velocity. Basal sediment availability declined during channelization but increased throughout the remainder of the monitored period, resulting in very efficient basal sediment evacuation over the peak of the melt season. Increased basal sediment availability during the summer appears to have been linked to high diurnal water pressure variation within subglacial channels inferred from the strong increase in flow velocity with discharge. Basal sediment availability therefore appears likely to have been increased by (1) enhanced local ice‐bed separation leading to extra‐channel flow excursions and[sol ]or (2) the deformation of basal sediment towards low‐pressure channels due to a strong diurnally reversing hydraulic gradient between channels and areas of hydraulically less‐efficient drainage. Copyright © 2005 John Wiley & Sons, Ltd.     

Trace element and major ion concentrations and dynamics in glacier snow and melt: Eliot Glacier,Oregon Cascades

We present concentrations of environmentally available (unfiltered acidified 2% v/v HNO₃) As, Cu, Cd, Pb, V, Sr, and major ions including Ca²⁺, Cl^?, and SO₄^2? in a July 2005 and a March 2006 shallow snow profile from the lower Eliot Glacier, Mount Hood, Oregon, and its proglacial stream, Eliot Creek. Low enrichment factors (EF) with respect to crustal averages suggests that in fresh March 2006 snow environmentally available elements are derived primarily from lithogenic sources. Soluble salts occurred in lower and less variable concentrations in July 2005 snow than March 2006. Conversely, environmentally available trace elements occurred in greater and more variable concentrations in July 2005 than March 2006 snow. Unlike major solutes, particulate‐associated trace elements are not readily eluted during the melt season. Additionally, elevated surface concentrations suggest that they are likely added throughout the year via dry deposition. In a 1‐h stream sampling, ratios of dissolved (<0·45 µm) V:Cl^?, Sr:Cl^?, and Cu:Cl^? are enriched in the Eliot Stream with respect to their environmentally available trace element to Cl^? ratios in Eliot Glacier snow, suggesting chemical weathering additions in the stream waters. Dissolved Pb:Cl^? is depleted in the Eliot Stream with respect to the ratio of environmentally available Pb to Cl^? in snow, corresponding to greater adsorption onto particles at greater pH values. Copyright © 2009 John Wiley & Sons, Ltd.     

Streamwater hydrograph separation in an alpine glacier area in the Qilian Mountains,northwestern China

ABSTRACT

Glacier-melt-induced changes in runoff are of concern in northwestern China where glacier runoff is a major source for irrigation, industries and ecosystems. Samples were collected in different water mediums such as precipitation, glacial ice/snowcover, meltwater, groundwater and streamwater for the analysis of stable isotopes and solute contents during the 2009 runoff season in the Laohugou Glacial Catchment. The multi-compare results of δ¹⁸O values showed that significant difference existed in different water mediums. Source waters of streamflow were determined using data of isotopic and geochemical tracers and a three-component hydrograph separation model. The results indicated that meltwater dominated (69.9 ± 2.7%) streamflow at the catchment. Precipitation and groundwater contributed 17.3 ± 2.3% and 12.8 ± 2.4% of the total discharge, respectively. According to the monthly hydrograph, the contribution of snow and glacier meltwater varied from 57.4% (September) to 79.1% (May), and that of precipitation varied from 0% (May) to 34.6% (September). At the same time, the monthly contribution of groundwater kept relatively steady, varying from 9.7% (June) to 20.9% (May) in the runoff season. Uncertainties for this separation were mainly caused by the variation of tracer concentrations. It is suggested that the end-member mixing analysis (EMMA) method can be used in the runoff separation in an alpine glacial catchment.

Editor Z.W. Kundzewicz; Associate editor Not assigned     

What are the climate controls on δD in precipitation in the Zongo Valley (Bolivia)? Implications for the Illimani ice core interpretation

Controversy has surrounded the interpretation of the water isotopic composition (δD or δ¹⁸O) in tropical and subtropical ice cores in South America. Although recent modeling studies using AGCM have provided useful constraints at interannual time scales, no direct calibration based on modern observations has been achieved. In the context of the recent ice core drilling at Nevado Illimani (16°39′S-67°47′W) in Bolivia, we examine the climatic controls on the modern isotopic composition of precipitation in the Zongo Valley, located on the northeast side of the Cordillera Real, at about 55 km from Nevado Illimani. Monthly precipitation samples were collected from September 1999 to August 2004 at various altitudes along this valley. First we examine the local and regional controls on the common δD signal measured along this valley. We show that (1) local temperature has definitely no control on δD variations, and (2) local rainout is a poor factor to explain δD variations. We thus seek regional controls upstream the Valley potentially affecting air masses distillation. Based on backtrajectory calculations and using satellite data (TRMM precipitation, NOAA OLR) and direct observations of precipitation (IAEA/GNIP), we show that moisture transport history and the degree of rainout upstream are more important factors explaining seasonal δD variations. Analysis of a 92-yr simulation from the ECHAM-4 model (T30 version) implemented with water stable isotopes confirms our observations at seasonal time scale and emphasize the role of air masses distillation upstream as a prominent factor controlling interannual δD variations. Lastly, we focus on the isotopic depletion along the valley when air masses are lifted up. Our results suggest that, if the temperature gradient between the base and the top of the Andes was higher by a few degrees during the Last Glacial Maximum (LGM), less than 10% of the glacial to interglacial isotopic variation recorded in the Illimani ice core could be accounted for by this temperature change. It implies that the rest of the variation would originate from wetter conditions along air masses trajectory during LGM.     

Flow prediction in high altitude ungauged catchments: A case study in the Italian Alps (Pantano Basin,Adamello Group)

We develop a simple model to evaluate the daily flow discharges in the ablation season for the 11 km² Pantano basin in the Retiche Italian Alps, based upon the data gathered during a three years field campaign. The Pantano basin embeds the Venerocolo debris covered and the Avio debris free glaciers, covering 2.14 km² in the Adamello Group, where the widest Italian glacier Adamello is located. First, degree-day models based upon air temperature are tuned to calculate snow and ice melt at daily scale. Glaciers’ meteorological data are collected from an automatic weather station (AWS), operating on the glacier during summer 2007. The melt factors in the debris covered areas of the glacier are estimated against debris thickness, using a data driven parameterization. The flow discharge from the catchment is estimated using semi distributed flow routing for the ablation seasons of four years, from 2006 to 2009. The predicted discharges are compared to those derived from inverse reservoir's routing at the Benedetto lake, catching the basin outflow. The proposed approach is valuable as a tool to investigate the hydrology of poorly gauged glacierized areas, including those with debris covered ice, widely diffused and yet poorly understood. Pending accurate parameterization the approach is usable for water resources evaluation and for long term assessment of the climate change impact on the glacierized areas within the Alps.