Study of isotopic seasonality to assess the water source of proglacial stream in Chhota Shigri Glaciated Basin,Western Himalaya

The impact of surface melt patterns and the Indian summer monsoon (ISM) is examined on the varying contributions of end member (snow, glacier ice, and rain) to proglacial streamflow during the ablation period (June–October) in the Chhota Shigri glaciated basin, Western Himalaya. Isotopic seasonality observed in the catchment precipitation was generally reflected in surface runoff (supraglacial melt and proglacial stream) and shows a shift in major water source during the melt season. Isotopically correlated (δ¹⁸O–δD) high deuterium intercept in the surface runoff suggests that westerly precipitation acts as the dominant source, augmenting the other snow- and ice-melt sources in the region. The endmember contributions to the proglacial stream were quantified using a three-component mixing. Overall, glacier ice melt is the major source of proglacial discharge. Snowmelt is the predominant source during the early ablation season (June) and the peak ISM period (August and September), whereas ice melt reaches a maximum in the peak melt period (July). The monthly contribution of rain is on the lower side and shows a steady rise and decline with onset and retreat of the monsoon. These results are persistent with the surface melt pattern observed in Chhota Shigri glacier, Upper Chandra basin. Moreover, the role of the ISM in Chhota Shigri glacier is unvarying to that observed in other glacierized catchments of Upper Ganga basin. Thus, this study augments the significant role of the ISM in glacier mass balance up to the boundary of the central-western Himalayan glaciated region.     

Effect of dust event timing on glacier runoff: sensitivity analysis for a Tibetan glacier

The impact of the timing of dust deposition on glacier runoff was evaluated using a glacier mass‐balance model with a newly improved scheme to track a dusted layer in a snow layer of a glacier. The lowering of surface albedo due to the dusted layer appearing leads to a drastic increase of glacier runoff even under the same meteorological conditions. Calculations of seasonal sensitivity, the relationship between dusted date and resulting runoff, have shown that dust deposition during a melting season might cause a drastic mass outflow from a glacier through changing the surface albedo during the melting season. Copyright © 2006 John Wiley & Sons, Ltd.     

Processes at the margins of supraglacial debris cover: Quantifying dirty ice ablation and debris redistribution

Current glacier ablation models have difficulty simulating the high-melt transition zone between clean and debris-covered ice. In this zone, thin debris cover is thought to increase ablation compared to clean ice, but often this cover is patchy rather than continuous. There is a need to understand ablation and debris dynamics in this transition zone to improve the accuracy of ablation models and the predictions of future debris cover extent. To quantify the ablation of partially debris-covered ice (or ‘dirty ice’), a high-resolution, spatially continuous ablation map was created from repeat unmanned aerial systems surveys, corrected for glacier flow in a novel way using on-glacier ablation stakes. Surprisingly, ablation is similar (range ~ 5 mm w.e. per day) across a wide range of percentage debris covers (~ 30–80%) due to the opposing effects of a positive correlation between percentage debris cover and clast size, countered by a negative correlation with albedo. Once debris cover becomes continuous, ablation is significantly reduced (by 61.6% compared to a partial debris cover), and there is some evidence that the cleanest ice (<~ 15% debris cover) has a lower ablation than dirty ice (by 3.7%). High-resolution feature tracking of clast movement revealed a strong modal clast velocity where debris was continuous, indicating that debris moves by creep down moraine slopes, in turn promoting debris cover growth at the slope toe. However, not all slope margins gain debris due to the removal of clasts by supraglacial streams. Clast velocities in the dirty ice area were twice as fast as clasts within the continuously debris-covered area, as clasts moved by sliding off their boulder tables. These new quantitative insights into the interplay between debris cover characteristics and ablation can be used to improve the treatment of dirty ice in ablation models, in turn improving estimates of glacial meltwater production. © 2020 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd     

Quantitative evaluation of different hydrological modelling approaches in a partly glacierized Swiss watershed

Mountain water resources management often requires hydrological models that need to handle both snow and ice melt. In this study, we compared two different model types for a partly glacierized watershed in central Switzerland: (1) an energy‐balance model primarily designed for snow simulations; and (2) a temperature‐index model developed for glacier simulations. The models were forced with data extrapolated from long‐term measurement records to mimic the typical input data situation for climate change assessments. By using different methods to distribute precipitation, we also assessed how various snow cover patterns influenced the modelled runoff. The energy‐balance model provided accurate discharge estimations during periods dominated by snow melt, but dropped in performance during the glacier ablation season. The glacier melt rates were sensitive to the modelled snow cover patterns and to the parameterization of turbulent heat fluxes. In contrast, the temperature‐index model poorly reproduced snow melt runoff, but provided accurate discharge estimations during the periods dominated by glacier ablation, almost independently of the method used to distribute precipitation. Apparently, the calibration of this model compensated for the inaccurate precipitation input with biased parameters. Our results show that accurate estimates of snow cover patterns are needed either to correctly constrain the melt parameters of the temperature‐index model or to ensure appropriate glacier surface albedos required by the energy‐balance model. Thus, particularly when only distant meteorological stations are available, carefully selected input data and efficient extrapolation methods of meteorological variables improve the reliability of runoff simulations in high alpine watersheds. Copyright © 2011 John Wiley & Sons, Ltd.     

Implications of decadal to century scale glacio‐hydrological change for water resources of the Hood River basin,OR, USA

In glacier‐fed rivers, melting of glacier ice sustains streamflow during the driest times of the year, especially during drought years. Anthropogenic and ecologic systems that rely on this glacial buffering of low flows are vulnerable to glacier recession as temperatures rise. We demonstrate the evolution of glacier melt contribution in watershed hydrology over the course of a 184‐year period from 1916 to 2099 through the application of a coupled hydrological and glacier dynamics model to the Hood River basin in Northwest Oregon, USA. We performed continuous simulations of glaciological processes (mass accumulation and ablation, lateral flow of ice and heat conduction through supra‐glacial debris), which are directly linked with seasonal snow dynamics as well as other key hydrologic processes (e.g. evapotranspiration and subsurface flow). Our simulations show that historically, the contribution of glacier melt to basin water supply was up to 79% at upland water management locations. We also show that supraglacial debris cover on the Hood River glaciers modulates the rate of glacier recession and progression of dry season flow at upland stream locations with debris‐covered glaciers. Our model results indicate that dry season (July to September) discharge sourced from glacier melt started to decline early in the 21st century following glacier recession that started early in the 20th century. Changes in climate over the course of the current century will lead to 14–63% (18–78%) reductions in dry season discharge across the basin for IPCC emission pathway RCP4.5 (RCP8.5). The largest losses will be at upland drainage locations of water diversions that were dominated historically by glacier melt and seasonal snowmelt. The contribution of glacier melt varies greatly not only in space but also in time. It displays a strong decadal scale fluctuations that are super‐imposed on the effects of a long‐term climatic warming trend. This decadal variability results in reversals in trends in glacier melt, which underscore the importance of long‐time series of glacio‐hydrologic analyses for evaluating the hydrological response to glacier recession. Copyright © 2016 John Wiley & Sons, Ltd.     

Supraglacial weathering crust dynamics inferred from cryoconite hole hydrology

Water levels in cryoconite holes were monitored at high resolution over a 3‐week period on Austre Brøggerbreen (Svalbard). These data were combined with melt and energy balance modelling, providing insights into the evolution of the glacier's near‐surface hydrology and confirming that the hydrology of the near‐surface, porous ice known as the ‘weathering crust’ is dynamic and analogous to a shallow‐perched aquifer. A positive correlation between radiative forcing of melt and drainage efficiency was found within the weathering crust. This likely resulted from diurnal contraction and dilation of interstitial pore spaces driven by variations in radiative and turbulent fluxes in the surface energy balance, occasionally causing ‘sudden drainage events’. A linear decrease in water levels in cryoconite holes was also observed and attributed to cumulative increases in near‐surface ice porosity over the measurement period. The transport of particulate matter and microbes between cryoconite holes through the porous weathering crust is shown to be dependent upon weathering crust hydraulics and particle size. Cryoconite holes therefore yield an indication of the hydrological dynamics of the weathering crust and provide long‐term storage loci for cryoconite at the glacier surface. This study highlights the importance of the weathering crust as a crucial component of the hydrology, ecology and biogeochemistry of the glacier ecosystem and glacierized regions and demonstrates the utility of cryoconite holes as natural piezometers on glacier surfaces. 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.     

Seasonal changes in the morphology of the subglacial drainage system,Haut Glacier d'Arolla,Switzerland

Dye tracing techniques were used to investigate the glacier-wide pattern of change in the englacial/subglacial drainage system of Haut Glacier d'Arolla during the ablation seasons of 1990 and 1991. Analysis of breakthrough curve characteristics indicate that over the course of a melt season, a system of major channels developed by headward growth at the expense of a hydraulically inefficient distributed system. By the end of the melt season, this channel system extended at least 3·3 km from the snout of the 4 km long glacier and drained the bulk of supraglacially derived meltwater passing through the glacier. The upper limit of the channel system closely followed the retreating snowline up-glacier. Rates of headward channel growth reached c. 65 m d⁻¹, although these rates decreased in the upper 1 km of the glacier where snowline retreat exposed a patchy firn aquifer. It appears that the removal of snow (with its high albedo and significant water storage capacity) from the glacier surface resulted in a dramatic increase in the volume of runoff into moulins, and in the peakedness of daily runoff cycles. This induced transient high water pressures within the distributed drainage system, which caused it to evolve rapidly into a channelised system. It is therefore likely that, at a local scale, channel growth occurred down-glacier from moulins, and that the overall up-glacier-directed pattern of channel formation was caused by the retreating snowline exposing new moulins and crevasses to inputs of ice-derived meltwater. Damping of diurnal melt inputs by storage in the firm aquifer accounts for the slowing of channel growth in the upper glacier. © 1998 John Wiley & Sons, Ltd.     

Long-term shifts in feedbacks among glacier surface change,melt generation and runoff,McMurdo Dry Valleys,Antarctica

Glaciers of the McMurdo dry valleys (MDVs) Antarctica are the main source of streamflow in this polar desert. Because summer air temperatures hover near 0°C small changes in the energy balance strongly affect meltwater generation. Here we demonstrate that increased surface roughness, which alters the turbulent transfer of energy between the ice surface and atmosphere, yields a detectable increase in meltwater runoff. At low elevations on the glaciers, basin-like topography became significantly rougher over 13 years between repeat lidar surveys, yielding greater melt. In contrast, the smoother ice at higher elevation exhibited no detectable change in roughness. We pose a conceptual model of the cycle of glacier surface change as a result of climate forcing whereby glacier surfaces transition from being dominated by sublimation to becoming increasingly melt-dominated, which is reversible under prolonged cool periods. This research advances our understanding of warm season effects on polar glaciers.     

Subglacial sediment production and snout marginal ice uplift during the late ablation season of a temperate valley glacier

Sediment export from glaciated basins involves complex interactions between ice flow, basal erosion and sediment transfer in subglacial and proglacial streams. In particular, we know very little about the processes associated with sediment transfer by subglacial streams. The Haut Glacier d'Arolla (VS, Switzerland) was investigated during the summer melt season of 2015. LiDAR survey revealed positive surface changes in the ablation zone, indicating glacier uplift, at the end of the morning during the period of peak ablation. Instream measures of sediment transport showed that suspended load and bedload responded differently to diurnal flow variability. Suspended load depended on the availability of fine material whereas bedload depended mainly on the competence of the flow. Interpretation of these results allowed development of a conceptual model of subglacial sediment transport dynamics. It is based upon the mechanisms of clogging (deposition) and flushing (transport/erosion) in sub-glacial channels as forced by diurnal flow variability. Through the melt season, the glacier hydrological response evolves from being buffered by glacier snow cover with a poorly developed subglacial drainage system to being dominated by more rapid ice melt with a more hydraulically efficient subglacial channel system. The resultant changes in the shape of diurnal discharge hydrographs, and notably higher peak flows and lower base flows, causes sediment transport to become discontinuous, with overnight clogging and late morning flushing of subglacial channels. Overnight clogging may be sufficient to reduce subglacial channel size, creating temporarily pressurized flow and lateral transfer of water away from the subglacial channels, leading to the late morning glacier surface uplift. However, without further data, we cannot exclude other hypotheses for the uplift. © 2018 John Wiley & Sons, Ltd.     

Modification and testing of a one‐dimensional energy and mass balance model for supraglacial snowpacks

A one‐dimensional energy and mass balance snow model (SNTHERM) has been modified for use with supraglacial snowpacks and applied to a point on Haut Glacier d'Arolla, Switzerland. It has been adapted to incorporate the underlying glacier ice and a site‐specific, empirically derived albedo routine. Model performance was tested against continuous measurements of snow depth and meltwater outflow from the base of the snowpack, and intermittent measurements of surface albedo and snowpack density profiles collected during the 1993 and 2000 melt seasons. Snow and ice ablation was simulated accurately. The timing of the daily pattern of meltwater outflow was well reproduced, although magnitudes were generally underestimated, possibly indicating preferential flow into the snowpack lysimeter. The model was used to assess the quantity of meltwater stored temporally within the unsaturated snowpack and meltwater percolation rates, which were found to be in agreement with dye tracer experiments undertaken on this glacier. As with other energy balance studies on alpine valley glaciers, the energy available for melt was dominated by net radiation (64%), with a sizable contribution from sensible heat flux (36%) and with a negligible latent heat flux overall, although there was more complex temporal variation on diurnal timescales. A basic sensitivity analysis indicated that melt rates were most sensitive to radiation, air temperature and snowpack density, indicating the need to accurately extrapolate/interpolate these variables when developing a spatially distributed framework for this model. Copyright © 2008 John Wiley & Sons, Ltd.     

Modelling suspended sediment concentration using artificial neural networks for Gangotri glacier

The dynamics of suspended sediment involves inherent non‐linearity and complexity because of existence of both spatial variability of the basin characteristics and temporal climatic patterns. This complexity, therefore, leads to inaccurate prediction by the conventional sediment rating curve (SRC) and other empirical methods. Over past few decades, artificial neural networks (ANNs) have emerged as one of the advanced modelling techniques capable of addressing inherent non‐linearity in the hydrological processes. In the present study, feed‐forward back propagation (FFBP) algorithm of ANNs is used to model stage–discharge–suspended sediment relationship for ablation season (May–September) for melt runoff released from Gangotri glacier, one of the largest glaciers in Himalaya. The simulations have been carried out on primary data of suspended sediment concentration (SSC) discharge and stage for ablation season of 11‐year period (1999–2009). Combinations of different input vectors (viz. stage, discharge and SSC) for present and previous days are considered for development of the ANN models and examining the effects of input vectors. Further, based on model performance indices for training and testing phase, a suitable modelling approach with appropriate model input structure is suggested. The conventional SRC method is also used for modelling discharge–sediment relationship and performance of developed models is evaluated by statistical indices, namely; root mean square error (RMSE), mean absolute error (MAE) and coefficient of determination (R²). Statistically, the performance of ANN‐based models is found to be superior as compared to SRC method in terms of the selected performance indices in simulating the daily SSC. The study reveals suitability of ANN approach for simulation and estimation of daily SSC in glacier melt runoff and, therefore, opens new avenues of research for application of hybrid soft computing models in glacier hydrology. Copyright © 2015 John Wiley & Sons, Ltd.     

The formation of supraglacial debris covers by primary dispersal from transverse englacial debris bands

Glaciological controls on debris cover formation are investigated from the perspective of primary dispersal of supraglacial debris across a melting ice surface. This involves the migration of angled debris septa outcrops across a melting, thinning glacier ablation zone. Three measures of a glacier's ability to evacuate supraglacial debris are outlined: (1) a concentration factor describing the focusing of englacial debris into specific supraglacial mass loads; (2) the rate of migration of a septum outcrop relative to the local ice surface; and (3) a downstream velocity differential between a slower septum outcrop and the faster ice surface velocity. Measures (1) and (2) are inversely related, while measure (3) increases down‐glacier to explain why slow‐moving, thinning ice rapidly becomes debris covered. Data from Glacier d'Estelette (Italian Alps) are used to illustrate these processes, and to explore the potential for debris cover formation and growth in different glaciological environments. The transition from a ‘clean’, transport‐dominated to a debris‐covered ablation‐dominated glacier is explained by the melting out of more closely‐spaced debris septa, in combination with the geometric interactions of angled septa and ice surface in a field of reducing flow and increasing ablation. The growth and shrinkage of debris covers are most sensitive to glaciological changes at glaciers with gently‐dipping debris‐bearing foliation, but less sensitive at high‐compression glaciers whose termini are constrained by moraine dams and other forms of obstruction. These findings show that a variety of debris‐covered glacier types will show a spectrum of response characteristics to negative mass balance. Copyright © 2013 John Wiley & Sons, Ltd.     

A spreadsheet‐based (Microsoft Excel) point surface energy balance model for glacier and snow melt studies

This paper describes a point surface energy balance model which runs within the Microsoft Excel spreadsheet package. The study incorporates a large amount of previous energy balance work and presents it in a useable form. The core model calculates the net shortwave and longwave radiation fluxes, the turbulent sensible and latent heat fluxes and the surface melt rate at a point on a melting ice or snow surface, from hourly inputs of incoming shortwave radiation, vapour pressure, air temperature and wind speed data. The latitude, longitude, slope angle, aspect, elevation, local temperature lapse rate, albedo and aerodynamic roughness of the study site, and the elevation of the meteorological station, can all be specified in the model. An output file containing the hourly and daily rates, and the totals of the energy fluxes is generated. The main advantages of the model are: first, that it requires only a PC or laptop computer running standard Microsoft Windows software, enabling it to be used at a desktop or in the field; and second, that it can be adapted quickly to different sites, meteorological data formats and other application requirements. Model calculations are compared with measured surface melt rates at five points on Haut Glacier d'Arolla, Switzerland, over a 115 day ablation period. Allowing for differences in shading between the meteorological station and the glacier, the root mean square error of the calculated melt rates is 2·0 mm day⁻¹ water equivalent melt (mean error +1·2 mm day⁻¹), for measured melt rates in the range 23 to 42 mm day⁻¹ water equivalent melt. Copyright © 2000 John Wiley & Sons, Ltd.     

From precipitation to runoff: stable isotopic fractionation effect of glacier melting on a catchment scale

Stable isotope variability and fractionation associated with transformation of precipitation/accumulation to firn to glacial river water is critical in a variety of climatic, hydrological and paleoenvironmental studies. This paper documents the modification of stable isotopes in water from precipitation to glacier runoff in an alpine catchment located in the central Tibetan Plateau. Isotopic changes are observed by sampling firnpack profiles, glacier surface snow/ice, meltwater on the glacier surface and catchment river water at different times during a melt season. Results show the isotopic fractionation effects associated with glacier melt processes. The slope of the δD‐δ¹⁸O regression line and the deuterium excess values decreased from the initial precipitation to the melt‐impacted firnpack (slope from 9.3 to 8.5 and average d‐excess from 13.4‰ to 7.4‰). The slope of the δD‐δ¹⁸O line further decreased to 7.6 for the glacier runoff water. The glacier surface snow/ice from different locations, which produces the main runoff, had the same δD‐δ¹⁸O line slope but lower deuterium excess (by 3.9‰) compared to values observed in the firnpack profile during the melt season. The δD‐δ¹⁸O regression line for the river water exhibited a lower slope compared to the surface snow/ice samples, although they were closely located on the δD‐δ¹⁸O plot. Isotope values for the river and glacier surface meltwater showed little scatter around the δD‐δ¹⁸O regression line, although the samples were from different glaciers and were collected on different days. Results indicate a high consistency of isotopic fractionation in the δD‐δ¹⁸O relationships, as well as a general consistency and temporal covariation of meltwater isotope values at the catchment scale. Copyright © 2013 John Wiley & Sons, Ltd.     

Characteristics and climatic sensitivities of runoff from a cold‐type glacier on the Tibetan Plateau

Model calculations are made in order to understand the characteristics and response to climate change of runoff from a cold glacier on the Tibetan Plateau. Some 20% of meltwater is preserved at the snow–ice boundary due to refreezing, since the glaciers in mid to northern Tibet are sufficiently cooled during the previous winter. Sensitivity to alterations in meteorological parameters has revealed that a change in air temperature would cause not only an increase in melting by sensible heat, but also a drastic increase in melting due to lowering of the albedo, since some of the snowfall changes to rainfall. In addition, it was suggested that a decrease in precipitation would cause a lowering of the surface albedo, with a resulting increase in the contribution of glacier runoff to the total runoff of river water. This study shows the first quantitative evaluation of the above effects, though they have been suggested qualitatively. The seasonal sensitivity of glacier runoff was examined by changing the dates given for a meteorological perturbation for a period of only 5 days. It was revealed that changes in both air temperature and precipitation during the melting season strongly affected glacier runoff by changing the surface albedo, though these perturbations only slightly altered the annual averages. Copyright © 2006 John Wiley & Sons, Ltd.     

Evolution of debris cover on glaciers of the Eastern Alps,Austria, between 1996 and 2015

Debris cover on glaciers is an important component of glacial systems as it influences climate–glacier dynamics and thus the lifespan of glaciers. Increasing air temperatures, permafrost thaw and rock faces freshly exposed by glacier downwasting in accumulation zones result in increased rockfall activity and debris input. In the ablation zone, negative mass balances result in an enhanced melt-out of englacial debris. Glacier debris cover thus represents a clear signal of climate warming in mountain areas. To assess the temporal development of debris on glaciers of the Eastern Alps, Austria, we mapped debris cover on 255 glaciers using Landsat data at three time steps. We applied a ratio-based threshold classification technique and analysed glacier catchment characteristics to understand debris sources better. Across the Austrian Alps, debris cover increased by more than 10% between 1996 and 2015 while glaciers retreated in response to climate warming. Debris cover distribution shows significant regional variability, with some mountain ranges being characterised by mean debris cover on glaciers of up to 75%. We also observed a general rise of the mean elevation of debris cover on glaciers in Austria. The debris cover distribution and dynamics are highly variable due to topographic, lithological and structural settings that determine the amount of debris delivered to and stored in the glacier system. Despite strong variation in debris cover, all glaciers investigated melted at increasing rates. We conclude that the retarding effects of debris cover on the mass balance and melt rate of Austrian glaciers is strongly subdued compared with other mountain areas. The study indicates that, if this trend continues, many glaciers in Austria may become fully debris covered. However, since debris cover seems to have little impact on melt rates, this would not lead to prolonged existence of debris-covered ice compared with clean ice glaciers.     

Effect of snowpack removal on energy balance,melt and runoff in a small supraglacial catchment

Modelling melt and runoff from snow‐ and ice‐covered catchments is important for water resource and hazard management and for the scientific study of glacier hydrology, dynamics and hydrochemistry. In this paper, a distributed, physically based model is used to determine the effects of the up‐glacier retreat of the snowline on spatial and temporal patterns of melt and water routing across a small (0·11 km²) supraglacial catchment on Haut Glacier d'Arolla, Switzerland. The melt model uses energy‐balance theory and accounts for the effects of slope angle, slope aspect and shading on the net radiation fluxes, and the effects of atmospheric stability on the turbulent fluxes. The water routing model uses simplified snow and open‐channel hydrology theory and accounts for the delaying effects of vertical and horizontal water flow through snow and across ice. The performance of the melt model is tested against hourly measurements of ablation in the catchment. Calculated and measured ablation rates show a high correlation (r² = 0·74) but some minor systematic discrepancies in the short term (hours). These probably result from the freezing of surface water at night, the melting of the frozen layer in the morning, and subsurface melting during the afternoon. The performance of the coupled melt/routing model is tested against hourly discharge variations measured in the supraglacial stream at the catchment outlet. Calculated and measured runoff variations show a high correlation (r² = 0·62). Five periods of anomalously high measured discharge that were not predicted by the model were associated with moulin overflow events. The radiation and turbulent fluxes contribute c. 86% and c. 14% of the total melt energy respectively. These proportions do not change significantly as the surface turns from snow to ice, because increases in the outgoing shortwave radiation flux (owing to lower albedo) happen to be accompanied by decreases in the incoming shortwave radiation flux (owing to lower solar incidence angles) and increases in the turbulent fluxes (owing to higher air temperatures and vapour pressures). Model sensitivity experiments reveal that the net effect of snow pack removal is to increase daily mean discharges by c. 50%, increase daily maximum discharges by >300%, decrease daily minimum discharges by c. 100%, increase daily discharge amplitudes by >1000%, and decrease the lag between peak melt rates and peak discharges from c. 3 h to c. 50 min. These changes have important implications for the development of subglacial drainage systems. Copyright © 2002 John Wiley & Sons, Ltd.