3.1. GPS documentation of ice motion

In the spring and early summer, we deployed five GPS units spaced evenly over 12.5 km of the ablation area center line, the topmost site being up-glacier from HCL (Fig. 1). Each Trimble Zephyr Geodetic antenna was attached to a triangular base fixed to a tripod of 3 m electrical conduit steam-drilled into the ice surface in the manner of Reference AndersonAnderson and others (2004). Data were recorded on Trimble 5700 receivers at 15 s intervals. Over the course of the summer, each GPS station needed to be re-drilled into the ice following significant ablation. During re-drills, local glacier velocities were extrapolated from those before and after the re-drill. During one of these reinstallations, on day 171, GPS3 was relocated 725 m to the west-northwest of its initial position, to place it more in line with the other receivers. It is this second position that is illustrated In Figure 1. A sixth GPS system, 1.5 km south of the glacier terminus, served as a base station for differential processing of GPS baselines. Table 1 describes aspects of each on-ice GPS site.

During post-processing, 24 hour files were broken into 1 hour increments and differentially, statically, processed using Trimble Geomatics Office and kinematically processed using the GIPSY software package from the NASA Jet Propulsion Laboratory (Reference GregoriusGregorius, 1996). Both processing routines revealed similar patterns of motion; Trimble Geomatics Office solutions are plotted in all figures. With a coordinate rotation, we separated the GPS displacements into along-mean-flow and across-mean-flow components. The standard deviation of horizontal, across-flow positions during periods with minimal high-frequency velocity fluctuations is ∼3 mm for the majority of our time series, half the hourly displacement for one of our slower-moving GPS stations. However, on day 191, the GPS base station began to lose power intermittently (usually at night), increasing the standard deviation of our horizontal position solutions to 60 mm. Therefore, in order to ease the interpretation of velocity fluctuations over periods of several hours and longer, we smoothed our entire horizontal position time series with the robust locally weighted scatter-plot smoothing algorithm, rloess, implemented with the software package MATLAB^™ (v. 7.9) (Reference Truffer, Motyka, Hekkers, Howat and KingTruffer and others, 2009). We used a 1 day smoothing window that preserved velocity fluctuations with periods greater than several hours but eliminated fluctuations with higher frequencies. The 1 day smoothing window has the potential to reduce the magnitude of the diurnal fluctuations by 25% when compared with a 0.5 day smoothing window. This smoothing algorithm allows us to describe the horizontal velocity of Kennicott Glacier between days 191 and 229, data we have not been able to present in previous work (Reference Bartholomaus, Anderson and AndersonBartholomaus and other, 2008).

Ice surface slope at the GPS receivers ranged from 0.026 to 0.037 m m⁻¹ (Table 1); if mean ice motion is surface-parallel, then any bed-parallel horizontal motion will result in lowering of the GPS antenna. We remove this vertical component of surface parallel motion from the vertical displacement records of each GPS receiver, and refer to the remainder of the vertical record as the ‘slope-corrected vertical displacement’ (Figs 2c and 3c). The vertical positions calculated without a base station (days 191–229) exhibit a clear bias towards lower elevations; thus, the smoothing algorithm could not be used for these vertical position data. We have shaded these positions more lightly when plotting them in Figure 2.

3.2. Measurements of glacier hydrology

Measurements relevant to glacier hydrology consisted of stage, water and air temperature, and water chemistry. The HCL and DFL stage sensors were installed ∼70 and ∼30 m above their respective lake bottoms. The depth of Erie Lake is unknown, but the stage sensor was likely placed on the order of 50 m above the lake bottom. On day 159, after DFL drained following its monotonic spring and early-summer filling, the DFL sensor was relocated to the bottom of the empty basin in anticipation of refilling during the HCL flood. A temperature sensor in Jumbo Lake recorded the presence or absence of water. Each of these sensors recorded data every 10 min.

A sonic ranging sensor was fixed to the west footbridge over the Kennicott River to measure stage every 15 min. The cross section of the Kennicott River did not change significantly between 1999/2000 field seasons (Reference Anderson, Longacre and KraalAnderson and others, 2003b) and the HCL jökulhlaup of 2006. After the 2006 jökulhlaup, the largest maximum discharge in 10 years, we resurveyed the river cross section. During the ∼24 hours of elevated discharge, the river scoured up to 3 m in the thalweg (deepest portion of the channel cross section). Thus, we calculated separate rating curves for the pre- and post-jökulhlaup portions of our record and transitioned smoothly between them with an arctangent function. Both rating curves are set by logarithmic velocity profiles (Reference JulienJulien, 1995, p.94) from which we calculate a depth-averaged velocity for various portions of the channel and sum the velocity–depth–width products to calculate the total discharge through the river channel. This method requires estimation of a channel bed roughness height; we use the 84th percentile boulder size, D ₈₄, observed during low flows to estimate the roughness height, z₀ = D ₈₄/10. These steps reproduce the rating curve developed during the 1999 and 2000 field campaigns by actual discharge measurements, and allow a smooth transition through the flood to a new rating curve that pertains to the remainder of the season.

Near the same footbridge where stage was measured, we also measured water electrical conductivity (abbreviated as conductivity) and temperature. Outlet water chemistry was further assessed with near-daily sampling at the location of the conductivity/temperature probe. We filtered our samples in the field (<0.45 μm) and later analyzed them for Cl⁻, NO₃ ⁻ and SO₄ ²⁻ concentrations.

Air temperature was recorded by radiation-shielded sensors 0.5–0.7 m above the ice surface at GPS1–GPS3 at 10 or 30 min intervals. We measured ablation at the support poles of each GPS station at irregular intervals ranging from 3 to 40 days.

4.1. Ice motion

Horizontal ice displacement reflects two types of motion: a steady, background component, and an unsteady component that varies on diurnal and seasonal timescales. Over our study area, the magnitude of the steady speed increases with distance from the glacier terminus (Table 1), whereas the magnitude of higher-frequency speed fluctuations is independent of position on the glacier. It is these fluctuations that we focus on.

Over the 61 (GPS5) to 91 day (GPS3) records, instantaneous (mean hourly) horizontal ice surface speeds ranged between 0.14 and 0.95 m d⁻¹ during non-jökulhlaup times, and are generally well correlated with each other (Fig. 2b). Average rates of surface horizontal motion over the entire period of observation range between 0.24 m d⁻¹ (GPS1, for 62.8 days) and 0.47 m d⁻¹ (GPS5, for 60.8 days).

Fig. 2.

Time series of data collected during the 2006 field season. Prior to day 188, vertical dashed lines identify two distinct glacier Modes, A and B. After day 188, Mode switches are indistinct. (a) Half-hour and 24 hour running average air temperature at GPS3. (b) Horizontal ice surface speeds calculated using smoothed displacements at each GPS receiver. Due to collapse of the GPS3 support platform, data for this receiver are not plotted after day 210. (c) Slope-corrected vertical displacement (displacement minus the surface-parallel trajectory) at each GPS receiver. During the failure of our base station (days 191–226), some vertical positions exhibit a clear negative bias. These positions are lightly colored. (d) Lake level record at DFL, HCL and Erie Lake. DFL stage is relative to the lake basin floor. HCL and Erie Lake stage are with reference to the stage at time of instrument deployment. Total HCL depth was ∼80 m. The depth of Erie Lake is not known. (e) Discharge, Cl⁻ concentration and adjusted electrical conductivity of the Kennicott River. We linearly adjusted the value of the electrical conductivity so that it plots on the same axis as the Cl⁻ concentration. Adjusted conductivity = conductivity/7.41–7.00.

Diurnal variations in the horizontal speed with amplitudes of up to 0.5 m d⁻¹ are present in portions of each of the five GPS records, although they are absent or quite weak at other times. The strength of the diurnal variation may be described by its diurnal velocity ratio: the range of the daily speed divided by the minimum speed for the day. This ratio is analogous to the percent speed-up sometimes used to describe the impact of surface water delivery on glacier flow (e.g. Reference Bartholomew, Nienow, Mair, Hubbard, King and SoleBartholomew and others, 2010). Prior to the jökulhlaup, the diurnal velocity ratio tended to vary in blocks of time with generally higher and lower values (Fig. 4b). The switches between time intervals with high and low diurnal velocity ratios occurred nearly synchronously across the instrumented reach of the glacier; we refer to these time periods as Mode A and Mode B (Reference Bartholomaus, Anderson and AndersonBartholomaus and others, 2008). Periods with strong diurnal velocity variations (e.g. diurnal velocity ratios greater than 0.4) are referred to as Mode A, whereas Mode B is characterized by smaller diurnal velocity ratios. The highest diurnal velocity ratios not coincident with the jökulhlaup immediately precede the drainage of HCL, from day 178 to day 186, when they reach values between 2.4 (at GPS2; Fig. 4b) and 1.6 (at GPS1; not shown). Typical diurnal velocity ratios during Mode B (e.g. between days 172 and 178) are <0.3. Following the jökulhlaup, the diurnal velocity ratio is more consistent and close to 0.3; they lack both the highest or lowest values present during the period before the jökulhlaup.

Uplift rate, following Reference Iken, Röthlisberger, Flotron and HaeberliIken and others (1983), is the rate of change of slope-corrected vertical displacement. Maximum uplift rate (slopes of vertical displacement time series; Figs 2c and 3c) coincided with maximum horizontal speed (Figs 2b and 3b); thus generally Mode A is associated with positive uplift rate, whereas Mode B is marked by zero or negative uplift rate. Until day 170, each GPS receiver recorded increases in slope-corrected vertical displacement. GPS3, the earliest installed receiver, rose (against its surface slope) at an average rate of 0.036 m d⁻¹ from day 136 until day 157, after which uplift rate decreased but remained positive. From day 170 until days 178–179, GPS2–GPS5 fell (i.e. negative uplift rate) (Mode B). The receivers showed renewed positive uplift rates from day 178–179 up to the drainage of HCL (Mode A). In contrast, GPS1 experienced nearly steady uplift rates for the entirety of its record, punctuated only by increased uplift rates during the drainage of HCL. After the drainage of HCL, GPS1 resumed the steady climb it had recorded prior to day 186. Several days after the jökulhlaup, on day 191, GPS2–GPS5 resumed moving roughly parallel to the ice surface for the remainder of the measurement period. Net slope-corrected vertical displacement over the measurement period ranged from −0.70 m (GPS5, over 60.2 days) to 2.00 m (GPS1, over 62.9 days).

Fig. 3.

Detail of observations during period immediately before and during the drainage of HCL. (a) Hydrology records showing HCL and DFL stage, Kennicott River discharge, and status of Jumbo Lake and a supraglacial ‘geyser’ observed near the Kennecott town site. The symbols for Jumbo Lake and the ‘geyser’ represent only the periods when they were full or active, respectively; no information exists about their stage or discharge. (b) Horizontal ice surface speeds at each GPS receiver. To allow for the resolution of higher-frequency changes in motion, the smoothing window used to filter the positions for the speed calculation of (b) is half the duration of that used for Figure 2b. (c) Slope-corrected vertical displacement of each GPS receiver, smoothed in the same manner as (b).

4.2. Kennicott River discharge record

As is typical of other glacier-fed rivers, the Kennicott River discharge fluctuated at several timescales (Figs 2e and 3a) (e.g. Reference Anderson, Fernald, Anderson and HumphreyAnderson and others, 1999). Discharge increased from a few m³ s⁻¹ when the winter ice broke out on day 130 to a summer mean of 150 m³ s⁻¹ by day 157. For the remainder of the summer, with the exception of the HCL jökulhlaup, discharge varied between ∼100 and 250 m³ s ⁻¹. Discharge also varied diurnally with an amplitude of 10 m³ s⁻¹, roughly half a day out of phase with air temperature on the glacier. The jökulhlaup observations are discussed in greater depth below, along with other related observations.

4.3. Water chemistry

Solute concentrations in the Kennicott River vary significantly over the period of observation. The concentration of Cl⁻ ranged from 0.26 to 16.03 mg L⁻¹, and tracks well with conductivity (Fig. 2e). Winter snowfall collected on day 134 from Root Glacier had 0.07 mg L⁻¹ Cl⁻, whereas Jumbo Creek carried 0.11 mg L⁻¹ Cl⁻, at base flow on day 229. This provides further evidence, as had been found before, that Cl⁻ in the Kennicott River comes not from snowfall and associated runoff, but from the glacier bed (Reference AndersonAnderson and others, 2003a,Reference Anderson, Longacre and Kraalb). Prior to the HCL jökulhlaup, Cl⁻ concentration in the river can be described as low or high (<3 mg L⁻¹ Cl⁻ or >5 mg L⁻¹ Cl⁻), and switched rapidly between these levels, roughly synchronized with the Mode A/B switches of the horizontal velocity record. After the jökulhlaup, Cl⁻ concentrations were less variable, at levels between the highs and lows measured earlier.

Electrical conductivity followed similar patterns until the jökulhlaup, when the conductivity sensor cable was severed and the sensor was lost. Low values of conductivity were within the range 50–70 μs cm⁻¹ and higher values ranged between 115 and 170 μs cm⁻¹ (Fig. 2e). Low conductivity and Cl⁻ concentrations coincide with Mode A; higher values coincide with Mode B.

4.4. Lake level records

The ice-dammed lakes along the margin of Kennicott Glacier offer an important tool for assessing glacier hydrology (Figs 2d and 3a). The timing of their annual spring–summer drainage allows one to track qualitatively the development and up-glacier propagation of a high-discharge, low-pressure conduit (Reference Anderson, Longacre and KraalAnderson and others, 2003b; Reference Kessler and AndersonKessler and Anderson, 2004). Once a lake empties, its basin is generally well connected to the subglacial hydrologic network and can therefore serve as a manometer during high subglacial water-pressure events (Reference Anderson, Longacre and KraalAnderson and others, 2003b). After filling through the spring, Jumbo Lake drained between days 148 and 151. The DFL pressure logger recorded a complicated series of fillings and drainings between day 143, when the level first fell significantly, and day 166, when the lake drained for the last time of the summer (excepting the HCL flood on day 186). HCL and Erie Lake began to drain on days 184 and 190, respectively; the draining of HCL resulted in the jökulhlaup. The Gates Lakes remained full throughout the observation period.

4.5. Weather records

Air temperatures at GPS1–GPS3 were quite similar and generally followed the same seasonal patterns. At GPS3, 90% of temperature readings were between 2°C and 11°C, although air temperature prior to day 160 occasionally fell below 0°C (Fig. 2a). Ablation rates varied between 0.038 and 0.087 m d⁻¹, with higher rates at lower elevations. The mean, midsummer ablation rate at GPS3 was 0.062 m d⁻¹.

Precipitation totaled 158 mm between days 132 and 230 at the National Weather Service McCarthy 3 SW weather station. The most intense storm brought 2.5 cm of rain during day 172–173. Neither this rain event, nor any other during the summer, had any discernible impact on glacier motion.

4.6. Response to the Hidden Creek Lake jökulhlaup

One of the primary draws to study ice motion and subglacial hydrology on Kennicott Glacier is the natural experiment performed annually when HCL drains 20–30 × 10⁶ m³ of water beneath the glacier over ∼3 days. On day 186 (5 July 2006), the first day of the HCL jökulhlaup, horizontal ice surface speeds reached maxima of 3.0 m d⁻¹, up to six times greater than mean speeds recorded during non-flood times (Fig. 3). At the time of HCL drainage initiation late on day 184, both Jumbo Lake and DFL basins were empty, with streams passing along their bottoms to enter tunnel mouths in Root Glacier. During the passage of jökulhlaup floodwaters through Kennicott Glacier, both of these lake basins filled for 24–30 hours before re-draining beneath the ice. In 2000 and 2006, discharge of the streams into the drained basins and chemistry of the water that filled the basin indicated that DFL filled with water emerging from the glacier during the jökulhlaup (Reference AndersonAnderson and others, 2003a).

On day 185, as the discharge from HCL was increasing but the Kennicott River had yet to flood significantly, each GPS receiver accelerated. The most pronounced responses were recorded below HCL, where in quick succession GPS3 speeded up, then GPS2 and finally GPS1. This wave of acceleration progressed down-glacier at an average speed of 1.3 km h⁻¹. Peak speed occurred in the same order: GPS3 first, followed by GPS2 and then GPS1. The responses of the two receivers at and above HCL (GPS4 and GPS5) were more muted, with only minor increases in speed above typical diurnal maxima at GPS5. Each GPS receiver attained a maximum speed of 1.2 (GPS5) to 3.0 m d⁻¹ (GPS1–GPS3). Speeds began to decline simultaneously around day 186.3, falling back towards pre-flood rates of 0.2–0.8 m d⁻¹.

In contrast to the horizontal speeds, for which the most significant jökulhlaup response was at GPS1–GPS3, GPS2–GPS4 experienced the most pronounced rapid increases in slope-corrected vertical displacement, up to a total of 0.55 m at GPS2 (Fig. 3c). There appear to be two phases to the jökulhlaup-related uplift history of GPS2–GPS4. The initial rate of uplift at day 185.5 was slow for each of these receivers. At day 186.1, the rate of uplift at GPS2 and GPS3 increased to 0.72 m d⁻¹, but declined slightly at GPS4. The timing of the maximum uplift varied slightly between stations (from day 186.5 to 186.8) and was attained earliest at up-glacier sites. The greatest horizontal speed coincided with the greatest uplift rate (Fig. 3). Initiation of collapse occurs first at GPS4 and progresses down-glacier. At GPS4 through GPS2, following the peak in uplift, roughly a day of high vertical relaxation rates (at roughly −0.50 m d⁻¹) is followed by significantly slower collapse until day 191.

The vertical response of GPS1 appears as a step of 0.20 m set within the otherwise steady rate of background uplift. This was accomplished by enhanced uplift rate of 0.27 m d⁻¹ that lasted from day 186.0 to 186.7, roughly coincident with the time of the greatest uplift rate at GPS2– GPS4. GPS5, 3 km up-glacier of HCL, displayed little or no vertical response that can be clearly attributed to the HCL flood. This most up-glacier receiver maintains a nearly steady negative uplift rate, falling 0.20 m between days 185 and 191.

The onset of horizontal deceleration towards baseline ice speeds coincided with the onset of the rapid increase in discharge of the Kennicott River. The flood peaked at 900 m³ s⁻¹ at day 186.9, before declining exponentially over ∼4 days back toward pre-flood discharges. The timing of down-glacier lake filling and draining is also closely associated with ice speed. The linear rising limb of DFL stage during the jökulhlaup nearly perfectly coincides with the rise of speed at the closest GPS receiver, GPS2. DFL stage and GPS2 speed peaked simultaneously. However, although ice at GPS2 began to decelerate shortly after attaining its maximum speed, DFL stage remained high, up to 45 m above its basin floor, until Kennicott River discharge began to decline from its maximum. Based upon our temperature logger record, Jumbo Lake began to fill 0.5 days after up-glacier DFL, and drained only 0.1 days after DFL drained.

As yet another proxy for water pressure in the glacier during the outburst flood, a supraglacial ‘geyser’, several meters in height above the ice surface, was observed between days 186.4 and 186.5 several hundred meters west of the town site of Kennecott (personal communication from N. Bay, 2006).

5.1. Conceptual model

Our conceptual model (Reference Bartholomaus, Anderson and AndersonBartholomaus and others, 2008) is similar to that used to explain observations at Haut Glacier d’Arolla (Reference Hubbard, Sharp, Willis, Nielsen and SmartHubbard and others, 1995; Reference Mair, Nienow, Sharp, Wohlleben and WillisMair and others, 2002) and Aletschgletscher, Switzerland (Reference Hock, Iken and WanglerHock and others, 1999). We focus on the daily-averaged rate of change of water volume stored within the glacier, 24 hour averaged ∂S/∂t, as it allows us to discern broader temporal patterns of water storage, beyond the daily fluctuations driven by diurnal melt variability. We interpret a positive 24 hour averaged ∂S/∂t as indicating that the subglacial hydrologic system is not of sufficient capacity to transmit the water delivered to it from the glacier surface or margins. Therefore, when the instantaneous ∂S/∂t is positive, the increasing stored water must be englacial (Mode A). Englacial water, in macroporous moulins and crevasses, exerts a pressure on the bed, reducing the effective ice pressure at the bed. The effective pressure affects glacier basal motion through an inverse power-law relationship (e.g. Reference BindschadlerBindschadler, 1983), which feeds back on the englacial water storage by opening up space at the glacier bed into which englacial water is drawn. This dynamic subglacial storage is often considered to be cavities in the lee of bedrock bumps or asperities, although basal crevasses are also possible (Reference Harper, Bradford, Humphrey and MeierbachtolHarper and others, 2010).

In this conceptual model, the state of the subglacial hydrologic system is of paramount importance in setting the basal motion response to changes in water supply: the glacier response to the input of water cannot be predicted without knowing the condition of the subglacial hydrologic system. In the present study, we use the 24 hour averaged ∂S/∂t as a proxy for this condition, and see good correspondence between positive 24 hour averaged ∂S/∂t and Mode A behavior (Fig. 4a). Only when the 24 hour averaged ∂S/∂t is positive will instantaneous increases in ∂S/∂t lead to increased basal motion and to surface motion above the background rate.

Fig. 4.

(a) Time series of normalized ∂S/∂t and non-normalized horizontal speed at GPS receivers, as in Figure 2b. Speeds greater than 1 m d⁻¹ are not shown. (b) Diurnal velocity ratio for each day in our time series at two representative GPS stations. GPS2 reaches a maximum diurnal velocity ratio of 8.6 on day 185. Also on day 185, GPS3 reaches a maximum ratio of 5.4. (c) Cl⁻ concentration and adjusted electrical conductivity, as in Figure 2e.

5.2. Model development

We explore the above conceptual model by recreating the systems and interactions outlined above in a numerical model. A sketch of our model is as follows, with the complete description of the model numerics following thereafter. A key feature of the model is the presence of two water reservoirs within the glacier. Water supplied to the glacier surface or sides fills a subglacial reservoir with dynamic capacity. Water beyond the volume of the subglacial reservoir fills an englacial reservoir that has static macroporosity and is well connected to the bed. The basal water pressure scales linearly with the volume of englacial water storage and this basal water pressure in turn controls the rate of basal motion. Any increases in basal motion can expand the volume of the subglacial reservoir, thereby drawing water down out of the englacial reservoir, exerting a negative feedback on basal motion. Here we envision a subglacial network of interconnected cavities, formed in the lee of bedrock bumps, which grow as a result of basal motion and collapse due to viscous creep. In reality, the exact mechanism by which water travels and is stored at the glacier bed is less important than the characteristic that increasing basal motion increases the available basal storage, a situation applicable to both hard- and soft-bedded glaciers.

To capture the essence of the velocity response to changing hydrology, we make several simplifying assumptions: that the hydrologic system is uniform across the glacier and that the water pressure in subglacial conduits does not significantly affect basal sliding. Other models (e.g. Reference ClarkeClarke, 2003; Reference Kessler and AndersonKessler and Anderson, 2004; Reference Pimentel, Flowers and SchoofPimentel and others, 2010; Reference SchoofSchoof, 2010) have successfully predicted the patterns of water pressure and ice motion observed at other glaciers without assuming uniformity. We, however, use mean water pressures in the glacial hydrologic system, aware that in reality they are spatially variable no doubt on several length scales. Although others have shown this is rarely the case (e.g. Reference Fudge, Harper, Humphrey and PfefferFudge and others, 2005), for most of our period of observation the instrumented reach of glacier operates in synchrony (as reflected in switches between Mode A and B states shown in Figure 2). We therefore treat the study reach of the glacier as a single, homogeneous slab, with no spatial variability in either the water balance or the control that water pressure exerts on basal motion. While we lose the opportunity to investigate the phasing among sites on the glacier, we gain the constraint of a water balance applicable to the entire instrumented reach.

Linked cavities are expected to cover much of a glacier’s bed. When supplied with water at a sufficient rate, pressurization of water in this linked cavity system can promote basal motion (e.g. Reference KambKamb, 1987; Reference AndersonAnderson and others, 2004). For this reason, our model focuses exclusively on cavity evolution, to the exclusion of conduits or water storage in any basal till. In our model, cavities grow by basal sliding and by melt, and collapse by viscous deformation of their ice roofs, which in turn is dependent on the effective pressure at the bed.

The model glacier is a rectangular prism with dimensions L × W × H chosen to approximate the instrumented reach of Kennicott Glacier, with sides of 15, 4 and 0.4 km respectively. The model glacier bed is given a uniform distribution of identically sized cavities, each reacting in the same way to the evolution of water storage. Following Reference Kessler and AndersonKessler and Anderson (2004), bedrock bumps are assumed to be h tall, w _c wide, and are spaced λ_x and λ_y in down-glacier and cross-glacier directions, respectively. Effectively, λ_x and λ_y define a unit cell of the bed (Fig. 5), containing one bedrock bump. One may define a dimensionless cavity-forming capacity of the glacier bed f = (hw _c/λ_x λ_y ). In the simulations to follow, our cavities are spaced evenly with λ_x = 10 m and λ_y = 10 m; each cavity has width w _c = 5 m, and height h = 1 m at the right-angled lip where ice decouples from the bump; hence f = 0.05.

Fig. 5.

Schematic of a representative portion of the model domain, in which the glacier bed is characterized by a set of bedrock bumps that serve to nucleate cavities. A unit cell of the glacier bed consists of one bump and its associated cavity, centered in a domain of width λ_y and length λ_x . Spacing between the bumps, the widths, and the heights of the bumps govern the cavity-generating capacity of the bed. Changes in total storage of water within the glacier system are governed by differences in the time series of water inputs through surface melt and through drainage of side-glacier lakes, Q _in, and of water outputs from the glacier, Q _out. Total water storage, S, is apportioned between subglacial and englacial reservoirs. Englacial water stored within the connected macroporosity of the glacier, φ, governs the height of the effective water table, z _w, which in turn controls water pressure, P _w, and hence the effective pressure, N, at the glacier bed. High water pressure promotes sliding, u _b (inset). Cavities also grow by melt of their roofs. Cavities shrink by viscous deformation of their ice roofs, the rate of collapse controlled by the effective pressure. A negative feedback exists in the system in that sliding enlarges cavity cross-sectional area, which draws down water from the englacial reservoir. This then reduces the water table height, z _w, and the associated water pressure, which reduces the sliding speed, u _b.

At each model time-step, subglacial storage, S _sub, calculated as the available cavity volume, is filled from the total volume of storage. The remainder of the total storage is englacial, S _en. We assume the glacier has a macroscale, effective porosity, φ, that is well connected to the subglacial cavity reservoir. In reality, crevasses and moulins are present on the surface of Kennicott Glacier in distinct bands. We are aware of at least two such bands, in the vicinity of GPS2 and GPS4. Given the connected porosity φ, englacial water rises above the glacier bed to a height z _w, similar to water rising in piezometers tapping a confined aquifer:

The effective pressure, N, is then the ice overburden pressure less the englacial water pressure:

where ρ _i and ρ _w are the densities of ice and water, respectively, H is the ice thickness and g is the acceleration due to gravity.

We assume that the glacier slides according to a power-law relationship (e.g. Reference BindschadlerBindschadler, 1983),

where u _b is the basal velocity, C and γ are parameters, τ _b is the basal shear stress and n is the flow-law exponent. Although there are limitations to this sliding law, we (as well as other authors: e.g. Reference Vieli, Funk and BlatterVieli and others, 2001; Reference SchoofSchoof, 2010) prefer the formulation in Equation (3) due to its simplicity and its continued success in describing observations across a number of different glaciers (e.g. Reference Iken and BindschadlerIken and Bindschadler, 1986; Reference Sugiyama and GudmundssonSugiyama and Gudmundsson, 2004). The sensitivity of this model to the three parameters γ, and φ is discussed in Section 5.5. Sensitivity to the cavity-forming capacity, f, is addressed in Section 6.2.

A representative cavity cross-sectional area, A _c, increases at a rate governed by the sliding speed and the melt rate of cavity walls, and declines at a rate governed by creep closure (Reference NyeNye, 1976):

where γ _R is Röthlisberger’s dimensionless constant (0.313) (Reference RöthlisbergerRöthlisberger, 1972; Reference Kessler and AndersonKessler and Anderson, 2004), L _f is the latent heat of fusion of ice, is the water outflow discharge per cavity in a transverse section of our model (Q _out is the Kennicott River discharge (Fig. 2e)),∇z _w is the hydraulic gradient across a cavity and A is Glen’s flow-law parameter. We take ∇z _w as the glacier surface slope divided by the cavity network sinuosity, assumed to be 3 (Reference Kessler and AndersonKessler and Anderson, 2004). The first term in Equation (4) represents the growth due to sliding, the second term represents the growth due to melting of cavity walls and the third term represents collapse due to creep.

5.3. Model results

To assess how well this model captures the expectations outlined above, we drive it with two different water budgets: a simple, synthetic case, constructed to mimic the character of the Kennicott Glacier water budget (Fig. 6a), and the positive-degree-day based 2006 Kennicott Glacier water budget (Reference Bartholomaus, Anderson and AndersonBartholomaus and others, 2008; Fig. 7a). In the synthetic case, two pairs of superimposed sine waves, representing diurnal and synoptic-scale fluctuations, make up both Q _in and Q _out. After several synoptic cycles, both the water inputs and outputs taper smoothly to zero to simulate the onset of winter (Fig. 6a).

Fig. 6.

Input for and results of our water-storage and basal-motion model: synthetic case. (a) Input and output components of the prescribed water balance. (b) ∂S/∂t, the difference between the inputs and outputs. (c) Volumes of water storage: subglacial storage, englacial storage and the total water storage (sum of subglacial and englacial). (d) Percent of local glacier flotation level reached by englacial water storage. (e) Glacier surface speed, predicted by Equation (3). The diurnal velocity ratios for days with peak u _b (i.e. 47 and 67) are 0.63.

Fig. 7.

Input for and results of our water-storage and basal-motion model: data-driven case. Panels are the same as Figure 6, but the model is driven by water inputs and outputs calculated at Kennicott Glacier (Reference Bartholomaus, Anderson and AndersonBartholomaus and others, 2008).

Both the synthetic and the data-driven model outputs capture some of the common characteristics of glacier motion (Figs 6e and 7e), particularly periods of rapid and diurnally variable motion. As determined by our choice of sliding law, the most rapid motion occurs when the englacial water pressure is closest to ice flotation. Englacial water responds more directly than subglacial water to ∂S/∂t (Figs 6c and 7c). Diurnal-scale fluctuations, amounting to 2–3% of the ice flotation level, are present throughout both records. When the englacial storage reaches a significant portion of the volume required for flotation, effective pressure is lowest, and additional, modest increases in the englacial water storage lead to diurnal variations in basal motion. In both model runs, cavity melt (second term of Equation (4)), contributed insignificantly to cavity enlargement.

During the simulated winter, after day 100, the rate of change of storage, ∂S/∂t, approaches zero and the volumes of englacial and subglacial storage stabilize; about 45% of the total storage is accounted for in the subglacial system, within cavities (Fig. 6c). Water pressure stabilizes at about 75% of ice flotation, intermediate between the highs and lows of the ‘summer’. The high, stable winter water pressure predicted by the model is close to winter values observed at other glaciers (e.g. 84% of flotation measured at Bench Glacier, Alaska (Reference Harper, Humphrey, Pfeffer, Fudge and O’NeelHarper and others, 2005)). Sliding speeds are maintained near their minimum rate, at 1.4 m d⁻¹, far below the 3.2 m d⁻¹ diurnal sliding maxima reached when ∂S/∂t was most positive.

The lag between maximum 24 hour averaged ∂S/∂t and maximum rates of basal motion varies between the synthetic and data-driven model runs in a subtle but important way. Despite having been run with identical parameter values, the amplitude of the synoptic englacial water-pressure fluctuations varies between the two runs. In the synthetic run, negative 24 hour averaged ∂S/∂t leads to a water-pressure drop from the flotation level to 45% of flotation, before rising again with positive averaged ∂S/∂t values (Fig. 6d). Due to the 10 days required for englacial water storage to recover from these synoptic falls in ∂S/∂t, neither maximum diurnal velocity ratios nor peak rates of basal motion occur simultaneously with maximum averaged ∂S/∂t, as has been described by Reference Bartholomaus, Anderson and AndersonBartholomaus and others (2008) and as is modeled using the Kennicott water balance (Fig. 7). Instead, the diurnally varying basal motion lags the switch to positive ∂S/∂t by 7–8 days and nearly coincides with the maximum total water storage. A somewhat similar lag is apparent in the Kennicott GPS data around day 160–165. At this time, the observed transition from Mode B to Mode A motion lags the transition to positive 24 hour average ∂S/∂t by 3 days (Fig. 4). In the data-driven model (Fig. 7), the water pressure varies less (55–90%), so that water pressures quickly reach the level at which they are sensitive to additional water inputs, and the lag between the onset of positive average ∂S/∂t and diurnal velocity fluctuations is shorter, typically 1–3 days. Therefore, the amplitude of synoptic and seasonal swings in englacial water pressure can strongly influence the strength of the link between the rate of change in water storage, ∂S/∂t, and basal motion. This influence is discussed further in Section 6.2.

5.4. Comparison between model results and observed motion at Kennicott Glacier

5.5. Model sensitivity

Our model requires the selection of three parameters that govern the relationship between subglacial hydrology and basal motion. Each of these three parameters, γ, and φ, is only weakly constrained by observations reported in the literature. Four reported values for γ, the exponent on the dependence of sliding speed on effective pressure, range from 0.18 to 0.4 (Reference JanssonJansson, 1995; Reference Sugiyama and GudmundssonSugiyama and Gudmundsson, 2004). has units that depend on the exponent γ, and varies from 1.5 × 10⁻⁵ m s⁻¹ Pa^0.18 to 1.1 × 10⁻³ m s⁻¹ Pa^0.4 (Reference JanssonJansson, 1995; Reference Sugiyama and GudmundssonSugiyama and Gudmundsson, 2004). Estimates of englacial macroporosity (crevasses/fractures and moulins/conduits) that allow the free exchange of water with the bed (our φ) are similarly limited. Estimates from Reference Bradford, Nichols, Mikesell and HarperBradford and others (2009) place an upper bound on the water content (micro-and macro-porosity) of Bench Glacier, 150 km southwest of Kennicott Glacier. Using ice-penetrating radar, they estimate a typical water content for the deeper portions of the glacier of 1.0–2.5%. On Storglaciären, Sweden, analyses of borehole video reveal that englacial voids make up 1.3% of the observed glacier volume (Reference PohjolaPohjola, 1994).

To test the sensitivity of our model, we ran the above synthetic experiment (Fig. 6) with values of γ between 0.05 and 0.7, between 3 × 10⁻⁶ and 3 × 10⁻³, and φ between 0.1% and 3%. As we are more interested in the timing and pattern of diurnal basal motion than the specific speeds, we present modeled diurnal velocity ratios as a metric of model success. In turn, we varied the values of these three parameters and ran our water-storage and basal-motion model with the synthetic water budget (Fig. 6a). For each set of the three parameters, the maximum diurnal velocity ratio of a post-spin-up synoptic cycle (e.g. days 40–60 of Fig. 6e) is presented in Figure 10. In each slice of the three-dimensional parameter space (i.e. each panel of Fig. 10), two parameters are varied, and the third is held at the value used to produce Figures 6 and 7 (γ = 0.22, and φ = 0.7%). As discussed above, we consider that diurnal velocity ratios of 0.1–5 indicate present and plausible diurnal velocity variations; this range is bracketed in bold black lines in Figure 10.

Fig. 10.

Sensitivity of diurnal velocity ratios to three parameters: γ, and φ. The diurnal velocity ratio ((max. diurnal velocity–min. diurnal velocity)/min. diurnal velocity) associated with the maximum synthetic velocity for the synthetic water balance shown in Figure 6a. Shading identifies the value of the diurnal velocity ratio; bold black lines bound the diurnal velocity ratios between 0.1 and 5, realistic values for diurnal velocity ratios when present. Three panels are slices through the three-dimensional parameter space; straight, white lines show the location of two other slices. Diurnal velocity ratio (a) as a function of γ and φ, ; (b) as a function of φ and , γ = 0.22; and (c) as a function of and γ, φ = 7 × 10⁻³. In (c), the values of and γ from other researchers are plotted with polygons. In each panel, the intersection of the thick, white lines represents the set of parameters used to create Figures 6 and 7. Small white dots show measurements that are the source of the contouring. Crenulated, intricate nature of the pattern of diurnal velocity ratios greater than 30 reflects limits of model stability, not anticipated dependence of the ratio on model parameters.

Parameter sets yielding diurnal velocity ratios of 0.1–5 span the range of observed parameter values. When provided with the synthetic water balance of Figure 6, three of the four pairs of and γ found in the literature will yield diurnal velocity ratios close to our bounds at an effective macroporosity of 0.7% (polygons in Fig. 10c). To bring the fourth pair of values, from Storglaciären, to a diurnal velocity ratio within our preferred range, φ would have to be >0.7%.

The depiction of parameter space in Figure 10 allows us to generalize the behavior of our model and discuss the general relationship between glacier hydrology and basal motion. The model is less sensitive to the selection of either or γ at low values of φ than at high values. For small or γ, the model is less sensitive to the selection of the other parameter while φ = 0.7%. Effective macroporosity greater than 2%, or γ not between 0.1 and 0.6, are combinations unlikely to allow a glacier to display diurnally fluctuating velocities, and would perhaps eliminate enhanced summer basal motion entirely. As enhanced summer basal motion is common to almost all glaciers, our sensitivity testing places some bounds on the expected range of these parameters, bounds that enclose the few values presented in previous studies (e.g. Reference JanssonJansson, 1995; Reference Sugiyama and GudmundssonSugiyama and Gudmundsson, 2004; Fig. 10).

6.1. Hydrology and basal motion of the glacier leading up to Hidden Creek Lake drainage

The drainage of HCL was preceded by the most rapid ice flow recorded outside the 1.5 days directly impacted by the passage of floodwaters (Fig. 2b). Beginning on day 179, one of the transitions to Mode A behavior, each day the glacier reach from GPS2 to GPS5 increased its speed en masse, with diurnal velocity ratios of approximately 2–3. Simultaneously, bed separation increased (Fig. 2c). Throughout most of our field season, the furthest up-glacier GPS sites exhibited the highest speeds, although during the interval from day 179 to 185 GPS4 both maintained a higher average speed and attained higher daily maxima (Fig. 3b) than any other GPS site, including GPS5, 3.5 km up-glacier of GPS4. Thus, longitudinal strain rates were compressional between GPS3 and GPS4 but, unusually, extensional between GPS4 and GPS5. Following HCL drainage, GPS5 resumed speeds in excess of those at GPS4 for the duration of our measurements; that reach of glacier therefore went back into compression.

The rate of lake stage rise did not begin to decline until day 183. Given that the hypsometry of HCL can be roughly approximated as a triangular trough, we see no evidence that the filling rate of HCL decreased, or that HCL began to drain, until several days after the start of GPS4’s particularly rapid motion on day 179. Thus, we do not associate the motion during days 179–183 (nor similarly during days 183–185) with early leakage of water from HCL. Instead, we interpret the motion during days 179–185 as likely being the result of an increase in air temperature beginning on day 179. The increased input of meltwater likely encountered a subglacial hydrologic system that was substantially collapsed after the preceding low temperatures (Figs 2a and 9).

The ingredients that favor initiation of a jökulhlaup include the presence of a nearby conduit, a path for water to leak between the lake and the conduit, and a high potential gradient between lake and conduit. Prior to warm weather beginning on day 179, cavities beneath the reach of glacier from GPS2 to GPS5, which include the reach adjacent to HCL, had substantially collapsed; bed separation was low (Fig. 9). Strong diurnal melt during days 179–185 then overwhelmed these cavities, leading to a period with strong diurnal velocity fluctuations and rapid basal motion (Fig. 7). We suggest that this sliding will increase the likelihood of a connection between lake and conduit (Reference Knight and TweedKnight and Tweed, 1991); this effect might have been augmented by crevassing observed between GPS4 and GPS5. Extensional strain rates visible between these stations in Figure 3b are consistent with crevassing. This increase in connectivity is seen in other analogous systems; for example, sliding promotes increases in lateral connectivity and effective permeability in fault systems (Reference Mallikamas and RajaramMallikamas and Rajaram, 2005). By this time in the melt season, it is likely that one or more conduits (Röthlisberger (R) channels) would have propagated to this elevation on the glacier, draining side-glacier lakes (e.g. Jumbo, then DFL) as it passes by. The few days of enhanced surface melt just prior to the jökulhlaup would have increased discharge of meltwater connecting to and then draining through this conduit system, which in turn would have lowered the pressure of the nearest channel as it grew to accept additional water (Reference RöthlisbergerRöthlisberger, 1972). Decreasing the R-channel water pressure increases the head gradient between it and the lake. Eventually, the hydrologic connection is established and water flow begins, inevitably leading to a jökulhlaup (Reference NyeNye, 1976).

6.2. Drawdown feedback leads to regulation of the glacier sliding system

Reference Iken and TrufferIken and Truffer (1997) described water-limited cavities as potential ‘sticky spots’ on the glacier bed. We may assess the ‘stickiness’ of isolated cavities by exploring the conceptual model outlined in Section 5.1. The growth of cavities through sliding increases the volume of subglacial storage and simultaneously draws down englacial water. The associated drop in water pressure will slow sliding, acting as a negative feedback or a regulator on the hydro-sliding system. Our intent in this section is to evaluate the relative importance of the various components involved in this feedback.

For a sliding law in the form of Equation (3), the rate of change of sliding speed with respect to water pressure is

where P _w is the water pressure (=ρ _w gz _w) and .

We must evaluate the rate at which new cavity space is generated by sliding. Recall the picture of the glacier bed with evenly spaced cavities and a defined cavity-generating capacity of the bed, f. Ignoring for the moment the viscous collapse of cavities, the rate of generation of basal cavity volume per unit area of bed is simply

where s, bed separation, has units of length. The basal sliding serves to lengthen all cavities. Again, we have ignored the loss of water to the channelized portion of the subglacial system.

If this newly generated subglacial cavity volume is immediately filled with englacial water from connected porosity, the volume of englacial water must decrease as cavities grow, but increase as new meltwater is added to the glacial system. The rate of change of englacial water volume per unit area of bed is therefore

where is the rate of water input to the glacier system (here assumed to be dominated by melt, with units of length over time). We may now assess how water pressure is reduced by sliding:

Inspection of this equation reveals that the greater the rate of melt input, the greater the rate of pressurization of the system (first term in parentheses). But the higher the sliding speed, the greater the rate at which water pressure is drawn down (second term in parentheses). Low porosity serves to amplify both effects: the smaller the porosity, the more sensitive the system. The cavity-generating capacity of the bed, f, also serves to amplify the drawdown effect.

The rate of change of sliding speed may now be cast as If we ignore melt inputs , and consider a closed system with only feedbacks between cavity opening and englacial water drawdown (as might be the case in winter), this becomes simply

The solution of this ordinary differential equation for the history of sliding is a simple exponential decay function with timescale

The shorter this timescale, the stronger the drawdown feedback and the more sensitive the sliding/hydrology system. Inspection of this equation illuminates the roles of the various components of the system. In the final expression, we isolate physical constants in the first parentheses, the variables that govern the strengths of the englacial and subglacial feedbacks in the second parentheses, and the state of the hydrologic system and sensitivity of sliding to this state in the third parentheses. In general, the more cavity volume that is generated per unit of sliding, the stronger the feedback; the greater the cavity-forming capacity of the bed, f (the taller and wider the cavities, and the smaller the spacing between cavities), the greater the cavity volume change per unit of sliding, and hence the greater the drawdown of the water table per unit of cavity volume. The smaller the porosity, the greater the drawdown, the amplifier of which we have already written.

The remaining dependence is on the effective pressure itself. The smaller the effective pressure, N (i.e. the greater the water pressure), the shorter the timescale and hence the greater the feedback. The more nonlinear the dependence of sliding speed on effective pressure, as represented by γ, the stronger the effect. This formulation yields additional insight into the switches between Modes A and B, in particular the presence or absence of diurnal velocity variations. In our model, N decreases towards zero (P _w approaches flotation; Fig. 6d) while 24 hour averaged ∂S/∂t is positive (Fig. 6b). Eventually, N may be sufficiently small such that τ falls below 1 day. At that point we can expect diurnal variations in sliding speed. In other words, the transition from Mode B to Mode A motion while water storage within the glacier is increasing (Reference Bartholomaus, Anderson and AndersonBartholomaus and others, 2008) is a natural result of an effective pressure-dependent ‘sliding law’ (Equation (3)).

What are plausible values for these numbers, and how quickly might sliding rates change? If we find that the timescale is of order 1 year, the feedback will not be strong, and the regulation will instead fall on other governors of the system that we have ignored in this analysis: variations in melt rate, cavity collapse, and rates of passage of water to and from the conduit system. If, on the other hand, the timescales are much shorter than this, the drawdown feedback becomes important in understanding seasonal to daily sliding histories.

Consider the following values, used to generate Figures 6 and 7: φ = 0.007 (0.7% porosity), g = 9.8 m s⁻², ρ _w = 1000 kg m⁻³, h = 1 m, w _c = 5 m, λ_x = λ_y = 10 m (hence f = 0.05), γ = 0.22, and . Say, for example, the water pressure is 95% flotation, or P _w = 0.95 P _i, hence N = 0.05 P _i. If the glacier thickness is 900 m, the timescale τ is then 17 days. If the glacier is thinner, 500 m, τ = 8 days. If the glacier is 900 m, but P _w = 0.99 P _i, τ = 2.4 days. If the glacier porosity is halved, τ = 1.2 days. As the timescale τ is inversely dependent on f, an order-of-magnitude increase in f from one glacier bed to another will result in order-ofmagnitude decrease in the expected response time. While we have not included in this discussion the addition of new melt from the glacier surface or the viscous collapse of cavities, we point out that their inclusion in a more comprehensive analysis would further shorten the response timescale, τ. It is clear from this scaling that the drawdown feedback may indeed be important on synoptic to daily timescales. The hydrologic state of the glacier and the nature of both bed and glacier should therefore matter.

With plausible values of the controlling parameters, this drawdown effect may yield very ‘sticky’ cavities, in which cavity growth can retard glacial motion, and therefore come into play as a regulator of the sliding system on temperate glaciers. Therefore, it ought not to be ignored in any hydro-sliding analysis. Depending on the geometry of the sub-glacial system, on the thickness of the glacier and on the state of the hydrologic system at the time (governing the connectedness of conduit-fast flow routes and cavity-slow flow routes), this drawdown feedback may become the dominant player.

The clearest example of the importance of this drawdown feedback occurs during the winter, when surface melt is shut off. In the absence of continued water addition to the glacier’s hydrologic system, the ice roof of the subglacial hydrologic system will collapse (cf. the last term on the right-hand side of Equation (4)). The winter subglacial hydrologic system becomes increasingly disarticulated, resulting in a ‘stickier’ glacier bed, despite the fact that basal water pressure in these cavities will rise to equal the ice overburden pressure. Any small amount of sliding associated with this high P _w will increase the cavity volume and result in drawdown of englacial water pressure. The negative feedback during winter is therefore identical to that described in the episodic shearing of granular materials, in which dilation upon shearing pulls down the water pressure that caused the shearing event in the first place (Reference Moore and IversonMoore and Iverson, 2002).

6.3. Weak connections between glacier hydrology and tidewater/outlet glacier dynamics

Although diurnal velocity fluctuations and seasonal speedups are nearly ubiquitous below the equilibrium line of terrestrially terminating glaciers and ice sheets far from marine margins (Reference WillisWillis, 1995; Reference Bartholomew, Nienow, Mair, Hubbard, King and SoleBartholomew and others, 2010), clear connections between glacier hydrology and glacier speed are rarely observed near the termini of lake-and marine-terminating glaciers (e.g. Reference Echelmeyer and HarrisonEchelmeyer and Harrison, 1990; Reference Howat, Tulaczyk, Waddington and BjörnssonHowat and others, 2008; Reference Joughin, Das, King, Smith, Howat and MoonJoughin and others, 2008). Our model suggests one potential explanation for this observation. As ice approaches the termini of lake-and marine-terminating glaciers, it most often experiences longitudinal strain rates that far exceed those observed on terrestrial glaciers. This is manifest in the last kilometers to tens of kilometers of glacier being chopped up by crevasses. We expect that this extensive crevassing ought to increase the macroporosity of the glacier above the threshold for which diurnal velocity fluctuations or seasonal melt-enhanced ice flow may occur. As stated above, that threshold may be as low as 2%, not far above the macroporosity measured at terrestrially terminating glaciers with less extreme strain rates (e.g. Reference PohjolaPohjola, 1994; Reference Bradford, Nichols, Mikesell and HarperBradford and others, 2009). Thus, while the mean water levels within gently sloping tidewater termini are at least as high as sea level (Reference PfefferPfeffer, 2007), the porosity of these glaciers may be so high that daily variations in water inputs are insufficient to change the subglacial water pressure substantially. If water pressures are insensitive to these water input fluctuations, then the rates of basal motion that they govern (Equation (3)) should also not vary on daily timescales.

6.4. Importance of incorporating hydrologic transience in ‘sliding laws’

The results of our study indicate that efforts to predict changes in basal motion from temperature parameterizations alone, or from predicted changes in meltwater modeled using such methods (e.g. Reference Parizek and AlleyParizek and Alley, 2004), will be unreliable unless they account for the evolution of the glacial hydrologic system (e.g. as in the present study, or Reference SchoofSchoof, 2010). The hydrologic system is transient both because the meltwater additions vary in time and space and because the hydrologic system varies in its efficiency in transmitting the meltwater. The strong annual sliding cycle that results from the combining seasonal variation in meltwater inputs and evolving glacial-hydrologic system, described here for an alpine glacier system, has now been documented on outlet glaciers from the Greenland ice sheet (Reference Bartholomew, Nienow, Mair, Hubbard, King and SoleBartholomew and others 2010; Reference SoleSole and others, 2011). The key point is that meltwater inputs alone do not control basal motion. Positive correlations between glacier acceleration and meltwater inputs (such as inferred from partial degree-day models) will occur in some cases, but in general such correlations oversimplify the system because they fail to incorporate the important role of dynamics of the subglacial hydrologic system.