Introduction

Concerns regarding the future viability of water resources in the Canadian Prairies (Reference Demuth, Keller, Demuth, Munro and YoungDemuth and Keller, 1997; Environment Canada, 2004) demonstrate that an improved monitoring network and understanding of processes regulating meltwater generation on alpine glaciers are needed. Glaciohydrological studies in the Canadian Rockies have largely been confined to Peyto Glacier (e.g. Reference YoungYoung, 1981; Reference MunroMunro, 1991; Reference Demuth and PietroniroDemuth and Pietroniro, 1999), with scattered results from other sites (Rae Glacier: Reference Smith, Guy and BarnardSmith, 1995). Thus, we present initial results and interpretations from a glaciohydrological study on the eastern slopes of the Canadian Rockies with an anticipated long-term monitoring program.

This research aims to (1) quantify relationships between hydrometeorological and glaciological variables and (2) identify surface processes affecting glacier discharge from an alpine glacier in the Canadian Rockies. The relationships between meteorological parameters and meltwater generation have been intensively studied (e.g. Reference BraithwaiteBraithwaite, 1995; Reference Gordon, Sharp, Hubbard, Smart, Ketterling and WillisHannah and McGregor, 1997; Reference Hannah and McGregorHannah and others, 2000; Reference Singh, Ramasastri, Kumar and ManoharSingh and others, 2000; Reference HodgkinsHodgkins, 2001; Reference Willis, Arnold and BrockWillis and others, 2002). We highlight the work of Reference Singh, Ramasastri, Kumar and ManoharSingh and others (2000) in particular, due to its use of lag correlation in evaluating relationships between discharge and meteorological parameters. Similar to Reference Jensen and LangJensen and Lang (1973), Reference Singh, Ramasastri, Kumar and ManoharSingh and others (2000) divided the discharge record into monthly intervals and found that lag times between air temperature and discharge varied throughout the season. One focus of this research is to examine the seasonal evolution of lag correlation relationships, and identify physical processes responsible for the observed changes.

Proglacial discharge from alpine glaciers generally exhibits pronounced increases in daily discharge maximum and range, with superimposed effects from precipitation events and freezing/snowfall occurrences (Reference Meier and TangbornMeier and Tang-born, 1961; Reference Smith, Guy and BarnardSmith, 1995; Reference FountainFountain, 1996; Reference Gordon, Sharp, Hubbard, Smart, Ketterling and WillisHannah and others, 2000; Reference Willis, Arnold and BrockWillis and others, 2002; Reference Verbunt, Gurtz, Jasper, Lang, Warmerdam and ZappaVerbunt and others, 2003). A critical issue in this regard is the role of the accumulation-area ratio (AAR) in the regulation of daily discharge range and totals. The AAR is the fraction of snow-covered glacier surface (0–1); changes in the AAR will produce significant changes in glacier-wide radiation reflectance and absorption, with consequences for runoff (Reference Smith, Guy and BarnardMoore and Demuth, 2001). Snow and firn also act as aquifers on a glacier (e.g. Reference StenborgStenborg, 1970), thereby serving to dampen discharge peaks and reduce the amplitude of diurnal flow variability. A decreasing AAR hence supports the rapid hydrological response of ice surfaces which do not store substantial amounts of meltwater, promoting increased diurnal variability. Model results at Haut Glacier d’Arolla, Switzerland, demonstrate that a reduction in AAR leads to a decreased lag between peak melt and peak discharge, and an increase in maximum discharge (Reference Willis, Arnold and BrockWillis and others, 2002). Two processes essential to evaluating seasonal glacier discharge are therefore embedded in the seasonal evolution of the AAR: first, it represents the change in snow/firn/ice partitioning of the surficial drainage network, and secondly, it constrains the average glacier-wide albedo and total absorbed radiation over the glacier surface. This study assesses the effect of snowline retreat, expressed through a modelled AAR, on the measured diurnal discharge of Haig Glacier, and its impact on linear multivariate regression analyses.

Field Site and Data Sources

Haig Glacier (50˚43′ N, 115˚18′ W), located approximately 100 km southwest of Calgary, Alberta, is the largest arm of a small icefield in the Front Ranges of the Canadian Rockies (Fig. 1). Mass-balance and meteorological measurements began in 2001 on the southeast arm of Haig Glacier, which has an area of 3.0 km² and ranges in elevation from 2800 m at the continental divide to 2450m at the toe. The main accumulation zone is a broad, gently sloping plateau located on the Continental Divide, which can receive upwards of 1800 mm snow water equivalent (SWE) annually. As a result of the low slope angles and modest relief of the glacier, slight variations in the equilibrium-line altitude (ELA) produce major changes in the accumulation area.

Fig. 1. Study area showing location of Haig Glacier (inset), automatic weather stations (squares), ablation-stake locations (solid dots), mass-balance density pits (open circles) and stream gauge (×). Bedrock contours in ma.s.l.

In 2002 and 2003, proglacial stream discharge was measured at a site 400 m downstream from the glacier toe (Fig. 1), where the single proglacial stream is channelized by bedrock. Stage height of the stream was measured with a Campbell Scientific SR-50 ultrasonic depth gauge at 10 min intervals, and glacial discharge records were calculated from rating curves created for each season from simultaneous stage and discharge measurements. Climate data were captured by two automatic weather stations (AWSs) recording 30 min averages of 10 s measurements of temperature, relative humidity, incoming and outgoing solar radiation (spectral range 0.3–2.8 μm), atmospheric pressure, wind speed and direction, and snow depth. The upper AWS is located along the longitudinal ablation monitoring line, near the midpoint of the glacier at 2640m (Fig. 1). The lower AWS is located below the toe of the glacier, roughly midway between the stream gauging site and the glacier terminus. Three weeks of data are missing from the upper AWS in 2003, as forest fires in the area forced an early field evacuation. Missing values at the upper AWS were replaced using regressions with the highly correlated lower-AWS data. Average surface albedo and daily average AAR values are obtained from the melt model developed for Haig Glacier by Reference AnslowAnslow (2004) (see Appendix).

Methods

Mean daily discharge (Q_i ), daily average temperature (T), total daily insolation (K), average daily detrended temperature (D _t), and standardized discharge range (Q _sr) are calculated for daily periods from 0000 to 2330 h, with

(1)

where Q _p is daily peak flow and Q _b is daily base flow (Reference Hannah and McGregorHannah and others, 2000).

A measure of radiation-detrended temperature was constructed by removing that part of the temperature signal that was linearly correlated with incoming solar radiation (Anslow, unpublished information). This provides a time series of ‘residual temperature’ which physically represents the temperature variability (hence, energy that is available for melt) associated with sensible and longwave heat fluxes. The resultant radiation-detrended temperature series is also independent of incoming solar radiation, avoiding statistical problems that arise due to correlation of solar radiation and actual temperature. It thus provides an estimate of the warming or cooling effects associated with boundary layer inversions, synoptic weather systems (i.e. the sensible-heat flux associated with different air masses), local winds (katabatic or valley winds that affect Haig Glacier) and the infrared radiative balance (e.g. cloudy vs clear-sky conditions).

We derive time series of daily average AAR and daily average surface albedo for the 2002 and 2003 ablation seasons using a blend of observations and model simulations. Hourly snow- and ice-melt values are simulated for the two seasons using a radiation-detrended temperature-indexed melt model (described in the Appendix and in greater detail by Reference AnslowAnslow (2004)), with daily totals calculated between 0000 and 2330 h during the ablation season.

To investigate seasonal variability of glacial meltwater storage and release, the transport efficiency (TE) of the glacier hydrologic system was calculated using total modelled glacier-wide daily melt (M _t) and observed total daily discharge (Q _t), both in units of m³:

(2)

Groundwater inputs and non-glacial snowmelt inputs to the proglacial stream were assumed to be negligible; the proglacial stream flows through a bedrock-incised channel with thin morainal debris, and Haig Glacier is surrounded by steep topography which is typically snow-free during the ablation season. Daily values of TE were examined for the 2003 ablation season, where stream-discharge data were sufficient to permit evaluation of melt–runoff relationships over most of the ablation season.

Statistical analysis

Analysis of the meteorological and hydrological datasets follows the two main purposes of the paper: lag autocorrelation is used to establish the varying strength of time-dependent relationships between radiation, discharge and temperature (Reference Singh, Ramasastri, Kumar and ManoharSingh and others, 2000); correlations between

daily discharge, standardized discharge range, temperature, insolation, albedo and AAR are investigated, and multivariate regression equations are derived to anticipate discharge forecasting from Haig Glacier.

Lag cross-correlation (r) establishes temporal relationships between two time series of data (u _t and y _t) through examination of sample cross-variance function c_uy and sample variances c_uu (0)c_yy (0) as a function of lag (k):

(3)

From the ablation season datasets, half-hour averages of temperature, radiation and discharge were divided into daily windows and 5 day windows to examine the seasonal evolution of the relationships between variables, and identify meteorological events which produce distinct deviations from the expected hydrologic response. Lag periods of up to 8 hours were investigated for both the daily and 5 day windows.

From the time series of mean daily discharge (Q_i ), mean previous day’s discharge (Q_{i -1}), mean daily temperature (T_i ), total daily insolation (K_i ), radiation-detrended temperature (DT) and modelled mean daily AAR, we conduct a suite of linear multivariate regression analyses. Normality of variables entered into a regression analysis was formally examined with the Kolmogorov–Smironov test (Reference Neter, Wasserman, Kutner and NachtseimNeter and others, 1996). All variables in 2002 passed tests for normality, but in 2003 this was only the case for the detrended temperature dataset.

Results and Discussion

Figure 2 displays summaries of the AWS datasets: snow depth, and incoming and outgoing radiation, albedo, and 5 day running means of temperature recorded between May 2002 and September 2003. Mean summer (June–August) temperatures were 5.5˚C in 2002 and 6.7˚C in 2003, and the 2002 melt season was punctuated by several episodes of below-freezing temperatures and minor snowfall events, whereas 5 day running means between 3 June and 9 September 2003 remained above 0˚C. As the snowline retreats past the upper AWS, rapid surface transitions are evident in the albedo record. A scatter plot of incoming and outgoing radiation between 1000 and 1600 h local time for the entire period (Fig. 3) further reveals distinct albedos for the three main surface types (fresh snow, aged snow/firn, and bare ice). Points which lie above the theoretical 1 : 1 line, i.e. values in the time series where reflected surface radiation is greater than incoming solar radiation, result from snowfall episodes where the sky-facing pyranometer is covered by snow.

Fig. 2. Summary of upper-AWS meteorological data: (a) daily average surface height (relative to end of summer 2002); (b) 5 day running mean temperature; (c) average daily incoming (dark curve) and outgoing (light curve) shortwave radiation; and (d) 5 day running mean of average daily albedo (calculated between 1000 and 1600 h).

Fig. 3. Scatter plot of half-hour incoming and outgoing solar radiation for points between 1000 and 1600 h, May 2002V September 2003.

Glacial discharge records, calculated from rating curves created for each season from simultaneous stage and discharge measurements, are shown in Figure 4. Logistical difficulties prevented the full collection of discharge for both seasons; the 2003 discharge record is significantly longer than that of 2002, with most of the ablation season captured in 2003. Differences are evident in the available data from each year, with 2003 demonstrating peak discharges of a much higher magnitude than 2002. Average baseflows calculated for the overlap period in each season show higher baseflows in 2002 (0.81m³ s^–1) than in 2003 (0.53m³ s^–1). In 2003 the dependent variable Q_i does not exhibit a normal distribution when all cases are considered. Discharge data for 2003 (Fig. 3) reveal a prominent shift in discharge mode near the end of the season. When the final 19 observed cases are excluded from the discharge dataset, the distribution passed formal normality testing, and significance levels attained in the subsequent trimmed regression analyses are valid. However, we also ran regression analyses with the inclusion of these 19 cases, and present the results in the following subsection for comparison. A significant question here is: does removal of these points increase the statistical strength of the regression model, or should the model be able to specify shut-off periods as well, despite the non-normality of the discharge data? Modelled AAR in 2003 likewise assumes a non-normal distribution; this, however, is a product of the hypsometry of the glacier itself, which leads to rapid changes in AAR when the snowline retreats over low-sloping sections on the upper glacier (Fig. 1). In 2002 this retreat only begins to occur by the end of the ablation season, and consequently the distribution of AAR assumes normality in 2002. As a key point of this investigation, however, we acknowledge the limitation of using non-normally distributed AAR data with regard to establishing significance levels on regression equations.

Fig. 4. 30 min averages of discharge from Haig Glacier proglacial stream (top) and flow indices Q _b, Q _sr, Q _i (bottom) for the 2002 and 2003 ablation seasons.

Glacier mass-balance values derived from ablation-stake surveys during the ablation season and from snow-depth and -density surveys during the accumulation season (cf. Fig. 1) are also substantially different for the two seasons (Table 1), 2002/03 being one of the most negative mass-balance years on record in the Canadian Rockies.

Table 1. Haig Glacier mass-balance values (mm.w.e.) for 2002 and 2003 mass-balance years

Detrended temperature values calculated from the regression of radiation and temperature are presented in Figure 5. The residual values (bottom panel) are entered in the multivariate regression analysis.

Fig. 5. Radiation, observed temperatures, predicted temperatures and residual values for the 2003 ablation season. A small amount of correlation between the residuals and the radiation series still exists (Table 3).

Cross-correlation analysis

Table 2 presents results from lag cross-correlation analysis of 5 day periods of average temperature, radiation and discharge for 2002 and 2003 at half-hour observation periods. Negative lag times indicate that peak discharge leads either temperature or radiation, thereby invalidating a causal relationship. Early-season discharge is more strongly correlated with temperature, yet both seasons evidence a shift to radiation-dominated discharge regimes by mid-summer. Snowfall events in 2002 (Julian days (JD) 244–248) and 2003 (JD 254–258) have a remarkable effect, causing the radiation–discharge correlation to collapse.

Table 2. Lag cross-correlation r values for discharge vs radiation and temperature in 2002 and 2003 ablation seasons

Notes: A negative lag value indicates that peak discharge occurs before radiation or temperature, thereby invalidating a causal relationship. Of temperature– and radiation–discharge correlations, the strongest periodic correlation is placed in bold to emphasize seasonal evolution of relationships. All correlations are significant at p = 0.05, n = 239.

Radiation–discharge lag times in 2003 decrease to 2.5 hours by JD 219–223, yet increase to 5 hours by JD 239–243, with no precipitation recorded during this period. A subsequent increase in lag time (JD 239–243) is the consequence of either an evolving glacier surface or a switch in subglacial drainage from a fast (channelized) to a slow (distributed) drainage network (Reference Nienow, Sharp and WillisNienow and others, 1996; Reference Gordon, Sharp, Hubbard, Smart, Ketterling and WillisGordon and others, 1998; Reference Willis, Arnold and BrockWillis and others, 2002). A crevasse zone in the bottom third of Haig Glacier offers numerous pathways for the quick transfer of melt from the surface to the subglacial drainage system. As the snowline retreats above this zone, melt is either distributed into up-glacier moulins with a longer subglacial routing, or travels along the surface to the lower crevasse zone, effectively increasing the lag time between peak radiation and discharge. Moulin development has been observed late in the ablation season, perhaps explaining the subsequent decrease in lag times that occurs on JD 249–253.

These observations from 5 day lag periods are borne out by an examination of daily lag cross-correlations (Fig. 6). The short duration of the 2002 discharge record prevents in-depth interpretation, yet compared with the 2003 daily cross-correlation record, the lack of a structured dominance of either radiation or temperature in 2002 presents the idea that AAR is an influential factor in discharge responses. In the 2003 record, the switch from a variable lag pattern to one of structured radiation-dominated discharge with a lag of approximately 2.5–3 hours occurs early in the record, when AAR values are still high. This is interpreted as either the establishment of a surface–bed connection, as crevasses in the lower part of the glacier are now open and draining snowpack runoff to the glacier bed (cf. Reference FountainFountain, 1996), or possibly the switch from a distributed subglacial network to one that is channelized. As the season progresses, fluctuations in daily lag pattern are consistent with varying incoming solar radiation inputs and a decreasing AAR, which gradually shortens the observed daily lag between peak radiation and peak discharge. The snowfall observed on JD 251 (AAR = 1) coincides with a return to the variable lag pattern, and thus a return to multi-day or variable discharge response. In 2002 the radiation-dominated discharge response is not observed, as a deep winter snowpack delayed the spring melt, and summer snowfall events disrupted the development of the supraglacial drainage system which was observed in 2003.

Single-variate correlation analysis

Results from univariate correlation analysis of Q_i , Q_{i ^– 1}, T_i , K_i AAR, Q _sr and DT are presented in Table 3. Strongest seasonal correlations are observed between Q_i and Q_i ^– ₁, with strong correlations also existing between daily average temperatures and discharge. Detrended temperature values offer an even stronger correlation with discharge than temperature itself, and so are preferentially employed in the multivariate regression analysis. Daily radiation–discharge correlations are comparatively lower than temperature–discharge correlations. In 2002, the greater winter mass balance (Table 1), lower summer temperatures and multiple summer snowfalls (Fig. 2) led to a low snowline position through much of the ablation season (cf. Fig. 6). The resultant high albedo reduced the efficiency of radiation in the daily production of meltwater, while higher baseflows resulted from meltwater storage and refreezing in the snow/firn aquifer (Reference Willis, Arnold and BrockWillis and others, 2002) that covered a large area of the glacier through the ablation season. Reference FountainFountain (1996) suggested that an observed seasonal increase in diurnal streamflow variation is the product of faster flow routing at the glacier surface, rather than an enlargement and increased efficiency of the subglacial hydraulic system (cf. Reference Nienow, Sharp and WillisNienow and others, 1996; Reference Willis, Arnold and BrockWillis and others, 2002). Low AAR values that persist through the 2003 season contribute to a much stronger AAR–S _qr correlation, and the timing of an increase in Q _sr coincides with rapidly decreasing AAR and albedo values. Correlations between AAR and standardized discharge range are negative and significant in both 2002 and 2003 (–0.39 and –0.67, respectively), suggesting that the surficial evolution is a strong contributing factor to the observed diurnal discharge.

Fig. 6. Daily lag/lead for temperature, radiation and discharge, 2002 and 2003 ablation seasons, with plotted glacier average albedo and AAR. Shaded areas indicate that discharge leads the variable in question (hence the relationship is null).

Table 3. 2002 and 2003 ablation season correlation matrix for

average discharge (Qi), average daily temperature (T), total daily insolation (K _t), AAR, daily detrended temperature (DT), daily standardized discharge range (S _qr) and previous day’s average discharge (Qi _–1)

Linear multivariate regression analysis

A suite of linear multivariate regression analyses for both ablation seasons is presented in Table 4. Regressions are generally more robust in 2003, though caution is warranted due to the non-normality of the AAR dataset. Below are listed several key observations regarding the multivariate regression analyses:

Table 4. Multivariate regression summary

• 1. The inclusion of all cases in 2003 regression analyses improves the regression results substantially. We suggest that as a regressive model, the performance is optimized when the model can accurately specify melt shut-down points through temperature and radiation data. An exploration of non-parametric regression methods is suggested as a way of surmounting this obstacle.

• 2. T and DT regressions against Qi produce inconclusive results, with a stronger regression from the detrended temperature values in 2003, and weaker regression in 2002. Note, however, that the use of DT values increases the Durbin–Watson score in both seasons. While the improved scores still demonstrate a high degree of autocorrelation, with the addition of multiple predictor variables the temperature data continue to cause an autocorrelation issue in the regression results. The use of detrended temperature data eliminates autocorrelation in these multivariate analyses, and is furthermore a normally distributed dataset in both seasons.

• 3. Using the previous day’s discharge, all regressions are improved substantially and autocorrelation is eliminated. Models without this term still perform robustly, and additionally offer the possibility of predicting glacier discharge using remotely sensed or commonly measured data (AAR, DT and K). Preferred models, according to criteria of normality and robustness, are shown in Table 4.

• 4. Table 5 shows a comparison of coefficient regression equations obtained in our study with those from other sites. Though it is not our preferred model, the strength of the relationships remains high, and coefficients are actually similar to those obtained at Mikkiglaciaren, Sweden (Reference StenborgStenborg, 1970). Preferred model coefficients are presented in Table 6, and interseasonal comparisons show that coefficients remain relatively stable, though the importance of AAR is emphasized in 2003 models.

Table 5. Comparison of ablation season multivariate regression coefficients, with regression coefficient (R ²) and standard error (SE)

Table 6. Summary of preferred regression model coefficients for constant (b) detrended temperature (bDT), radiation (bK), AAR (bAAR) and previous day’s mean discharge (bQi-1 ). For comparison between ablation seasons, the same regression models for 2002 and 2003 are shown sequentially

Transport efficiency analysis

The TE (Fig. 7), which describes the percentage of daily modelled melt that is captured by daily discharge measurement (Equation (2)), exhibits a long-term increase through the ablation season, and again this is matched by a decrease in the AAR, further suggesting the impact of snow/ice/firn partitioning on diurnal flow variability from an alpine glacier (Reference Willis, Arnold and BrockWillis and others, 2002). The TE record becomes dominated by values greater than 100% around JD 230, roughly the timing of the observed disappearance of a saturated ‘slush’ zone from mid-glacier. Precipitation events are superimposed on the overall increasing trend, and can lead to TE > 100% in some cases (e.g. JD 190, 195 and 237). Most notable in the TE record is the ‘purging’ of englacial and subglacial melt storage, as a snowfall event and negative temperatures (JD 251) effectively shut off modelled melt, resulting in high transport efficiencies as actual discharge drops off at a slower pace. Reference StenborgStenborg (1970) also observed the delayed flow in the early part of the ablation season, but saw evidence for a major evacuation of stored water in mid-season, with a switch to closely matching observed discharge and modelled melt by the end of the ablation season.

Fig. 7. (a) Total daily modelled melt and total daily discharge; (b) transport efficiency; (c) daily precipitation total; and (d) AAR.

Notes: A standard error (SE) of 20 000m³ is approximately 10% of average daily flow. Dependent variable is mean daily discharge (Q_i ). Predictors: T = mean daily temperature, K_i = total daily incoming solar radiation, AAR = accumulation-area ratio, DT = detrended temperature residual, Q_{i– 1} = previous day’s total discharge. Trim dataset in 2003 is all cases minus the final 19 cases of the discharge record. D-W is the Durbin–Watson statistic, with decision of autocorrelation based on number of cases and variables in analysis. Collinearity decision is based on condition indices of eigenvalues (not shown). All variables assume normal distributions in 2002; in 2003 only DT and Q_{i (trim)} are normal.

Available evidence from Haig Glacier points towards the evolution of the supraglacial drainage system as the leading control on the observed variability in proglacial discharge values. To corroborate this statement, the subglacial hydrology would need to be evaluated via dye-tracing experiments (e.g. Reference Nienow, Sharp and WillisNienow and others, 1996) or borehole observations (e.g. Reference FountainFountain, 1996; Reference Gordon, Sharp, Hubbard, Smart, Ketterling and WillisGordon and others, 1998).

Conclusions

We present glaciological and hydrometeorological relationships determined through correlations, multivariate regressions and lag cross-correlations of total daily discharge (Q_i ) vs total daily incoming solar radiation (K_i ), average daily temperature (T), average daily temperature residual (DT) and AAR. We also explore the idea of transport efficiency. Evolution of the ablation seasons shows a shift from temperature-dominated to radiation-dominated discharge as the AAR decreases, leading to (i) a decrease in snow and firn cover and associated water storage and (ii) increased absorption of solar radiation. Standardized daily discharge range is strongly correlated to the AAR, as increased ice exposure signifies increased efficiency of the glacier surface hydrology.

Previous day’s discharge, detrended temperature and AAR present the strongest seasonal correlations for discharge, and linear multiple regression models offer insight into hydrometeorological relationships. Detrended temperature values decrease autocorrelation in regression analysis, and present normal distributions where temperature datasets might not. The strength of the relationships between discharge, temperature and AAR suggest portability of the discharge regression model to other alpine sites, though these relationships must depend in some part on the crevasse/moulin configurations of a particular glacier and the hypsometry, which will affect AAR values as the snowline retreats.

In conclusion, the results of this study may provide direction for current melt modelling efforts. The use of a detrended temperature series offers a way of eliminating the collinearity of radiation and temperature for multivariate regression models while simultaneously reducing the autocorrelation of a regressed time series. Previous day’s discharge is a strong predictor variable for glacial discharge

and the autocorrelation of the discharge time series is indicative of the non-stationarity of the discharge system. We use a modelled AAR variable and radiation-detrended temperature to provide a robust independent variable for discharge predictions. While AAR provides an indirect measure of average glacier albedo, it contains additional information that pertains to delays in runoff associated with the supraglacial snow and firn aquifer. AAR therefore represents a bulk parameter indirectly measuring the glacier-scale supraglacial drainage efficiency.

Most features of the daily and seasonal evolution of runoff at Haig Glacier can be attributed to evolution of the supraglacial drainage system and surface energy balance through the ablation season. The extent of the winter snow-pack and the frequency of summer snowfalls proved very important to the drainage-system evolution and the total runoff from this glacier in the 2002 and 2003 ablation seasons. These effects can be broadly captured in runoff modelling through inclusion of AAR and average glacier albedo indices (via K) in a multivariate regression, accelerating efforts to develop improved hydrological forecasts and runoff models from glaciers of the eastern slopes of the Canadian Rockies. Meltwater runoff from glacierized catchments in the region is an important component of the hydrological budget in the drought-prone southern Canadian prairies.