Introduction

Enhanced surface meltwater runoff since the 1960s (Reference Hanna, Huybrechts, Janssens, Cappelen, Steffen and StephensHanna and others, 2005) has been implicated in ice flow acceleration across the margins of the Greenland ice sheet (GrIS) (Reference Rignot and KanagaratnamRignot and Kanagaratnam, 2006). This association has led to the concept of a direct hydromechanical link between the ice-sheet surface and its bed (Reference Zwally, Abdalati, Herring, Larson, Saba and SteffenZwally and others, 2002), a process well understood on valley glaciers but that until the last decade was undocumented on ice sheets. Covariations in ice velocity and surface uplift have been shown to lag peak daily melt in summer (Reference Van de WalVan de Wal and others, 2008; Reference Shepherd, Hubbard, Nienow, McMillan and JoughinShepherd and others, 2009) and occur during discrete pulses of rapid supraglacial lake drainage (Reference DasDas and others, 2008; Reference DoyleDoyle and others, 2012).

RADARSAT interferometric synthetic aperture radar (In-SAR) and GPS measurements covering the western margin of the GrIS confirm patterns of seasonal speed-up (Reference Joughin, Das, King, Smith, Howat and MoonJoughin and others, 2008; Reference BartholomewBartholomew and others, 2011; Reference Palmer, Shepherd, Nienow and JoughinPalmer and others, 2011). However, differences exist regarding the magnitude as well as the spatial and temporal patterns of these flow perturbations. Short-lived high-magnitude diurnal and seasonal increases in velocity (∼220% above the winter mean) have been recorded (Reference Shepherd, Hubbard, Nienow, McMillan and JoughinShepherd and others, 2009; Reference Bartholomew, Nienow, Mair, Hubbard, King and SoleBartholomew and others, 2010, Reference Bartholomew2011); however, a strong contrast between the duration of peak summer flow and the strength of the melt season (Reference Sundal, Shepherd, Nienow, Hanna, Palmer and HuybrechtsSundal and others, 2011) suggests the relationship between melt intensity and annual ice flux is more complex than originally thought. Long-term GPS measurements in this region have revealed a significant decrease in mean surface velocity despite increased ablation over a 17 year time-span (Reference Van de WalVan de Wal and others, 2008). A decade’s worth of satellite derived velocity data encompassing nearly all of Greenland’s outlet glaciers (200+) show complex spatial and temporal patterns of ice flow, highlighting the contrasting ice dynamics of land- and marine-terminating outlets (Reference Moon, Joughin, Smith and HowatMoon and others, 2012). Although investigations of individual glaciers using GPS surveys reveal a great deal about glacier behavior on short and interannual timescales, understanding the complex spatiotemporal behavior of the GrIS cannot be achieved without the help of satellite imagery.

Here we present high-resolution, 11 day repeat interval ice surface velocity maps covering Russell Glacier and three adjacent outlet glaciers (Fig. 1) during two melt seasons of contrasting intensity: the record high melt year of 2010, and the comparatively low melt year of 2009 (termed average by Reference Van As, Hubbard, Hasholt, Mikkelsen, Van den Broeke and FaustoVan As and others, 2012).

Fig. 1. Winter mean (2009 and 2010) velocity map for our study area derived from TSX data overlaid onto a Landsat 7 image from June 2008. The central flowline of Russell Glacier visible is shown in red, with flowlines of neighboring glaciers Isungata Sermia, outlet A and outlet B shown in black. The winter mean velocity map was derived from TSX data with mid-dates 11 November 2009, 7 February 2010, 9 November 2010 and 20 November 2010. Elevation contours were derived from a 2008 SPOT DEM.

Study Area

Russell Glacier is a land-terminating outlet glacier in West Greenland (67° N, 50°W; Fig. 1) that has received considerable scientific attention given its ease of access from nearby Kangerlussuaq. This land-terminating outlet is isolated from marine influences and exhibits changes in ice dynamics that are remote from tidewater influences and are directly attributable to surface melt forcing.

Methods

TerraSAR-X (TSX) data were acquired in StripMap mode producing images with 30 × 50 km swaths at 3 m resolution, with repeat orbits every 11 days (Reference Joughin, Smith, Howat, Scambos and MoonJoughin and others, 2010a). Two overlapping TSX image pairs were used covering an area of 3047 km² including an ice extent of 2400 km² (Fig. 1; Table 1a). Envisat advanced synthetic aperture radar (ASAR) data from 2009 were acquired in image mode, with a swath of 56−100 km, a repeat orbit time of 35 days and a 30 m spatial resolution (Table 1b). Additionally, panchromatic Landsat Thematic Mapper (TM)/Enhanced TM Plus (ETM+) and multispectral Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) scenes were acquired. Both Landsat and ASTER have a 15 m spatial resolution and 16 day repeat interval (Table 1c).

Table 1. Temporal coverage, image IDs and details of (a) the TSX dataset, (b) Envisat ASAR data and (c) ASTER and Landsat data. Date format is day/month/year

Velocity data were derived from TSX datasets using a combination of conventional interferometry and speckle-tracking techniques (Reference JoughinJoughin, 2002). Interferograms were produced using Gamma software (Reference Werner, Wegmüller, Strozzi Tand Wiesmann and Sawaya-LacosteWerner and others, 2001) where baselines were estimated by least-square fit from off-ice control points. Cross-correlation (speckle-tracking) procedures were used to estimate the range and azimuth offsets between interferometric pairs to determine displacements with sub-pixel accuracy. Error margins in the TSX velocity data were estimated to be ±5ma⁻¹ by measuring the mean difference from zero of stationary (off-ice) areas in each of the images. Additional velocity data from cross-correlation feature-tracking techniques using Envisat ASAR data (Table 1b) in Gamma, and panchromatic Landsat and ASTER scenes (Table 1c) in COSI-CORR software (Reference Leprince, Barbot, Ayoub and AvouacLeprince and others, 2007) were derived to fill gaps in the temporal coverage of TSX data during 2009. Feature tracking relies on finding the best-fit positions between search windows from two images allowing surface displacement to be derived. This method is described in detail by Reference Strozzi, Luckman, Murray, Wegmuller and WernerStrozzi and others (2002) and has been successfully applied to outlet glaciers in East Greenland (Reference Luckman, Murray, Jiskoot, Pritchard and StrozziLuckman and others, 2003). In this study, search window sizes of 64 pixels were used, which strike a balance between accommodating for the maximum expected displacement between image pairs and exploiting the spatial coverage. Although the spatial resolution is 15 m, we co-register the images with sub-pixel accuracy to produce velocity results with uncertainties in the order of 30 m a⁻¹, derived by measuring static, off-ice displacement.

Surface runoff is derived from a surface energy-balance model described in detail by Reference Van As, Hubbard, Hasholt, Mikkelsen, Van den Broeke and FaustoVan As and others (2012), utilizing data from a network of automated weather stations. The model uses temperature, air pressure, humidity, wind speed, and the downwelling components of shortwave radiation and longwave radiation to calculate the surface energy budget components (absorbed shortwave radiation, net longwave radiation, sensible heat flux, latent heat flux and subsurface heat flux). The net residual positive energy flux excess is applied to melting snow and ice, allowing surface melt to be calculated, with error margins of ±15% (Reference Van As, Hubbard, Hasholt, Mikkelsen, Van den Broeke and FaustoVan As and others, 2012). We compare the modeled surface runoff to the proglacial discharge measured at Watson River bridge (Reference Hasholt, Mikkelsen, Nielsen and LarsenHasholt and others, 2012).

Results

Surface velocity maps produced by the TSX data reveal complex spatial and temporal patterns of ice flow within Russell Glacier catchment (Fig. 2). To illustrate these patterns, surface velocity profiles for Russell Glacier, Isungata Sermia and two unnamed outlets, referred to as outlets A and B herein (Fig. 1), were extracted. The flowlines were determined from a combination of (1) structural mapping from high-resolution Satellite Pour l’Observation de la Terre (SPOT) images acquired in 2008, (2) velocity maps from previous studies (Reference Palmer, Shepherd, Nienow and JoughinPalmer and others, 2011; Reference Sundal, Shepherd, Nienow, Hanna, Palmer and HuybrechtsSundal and others, 2011) and (3) a preliminary basal digital elevation model (DEM) (personal communication from R. Pettersson, 2010). Velocities extracted along the flowline of Russell Glacier exhibit great temporal variation, with seasonal changes extending at least as far inland as 57 km, the limit of our velocity coverage (Fig. 2). Mean winter velocities for Russell Glacier are ∼100m a⁻¹ (Fig. 3), with summer speeds exceeding 200 m a⁻¹. All four of the outlet glaciers have similar hypsometry (Fig. 3); however, the extent and magnitude of accelerated flow exhibited by each outlet varies significantly. High-magnitude perturbations in ice flow experienced by outlets A and B (∼200% and ∼300% respectively) sharply contrast with the relatively modest (100−150%) increase in ice flow from Isungata Sermia and Russell glaciers in May 2010 (Fig. 4). By June 2010, mean surface velocities of outlet A had returned to or dropped below the winter mean, while the other outlets were still experiencing accelerated flow. In contrast to Russell Glacier, ice flow within the first 400 m elevation of Isungata Sermia continued at elevated rates above the winter mean for a longer period of the summer (Fig. 4).

Fig. 2. Surface velocities expressed as a percentage relative to the winter mean throughout summer 2010 over 11 day periods with midpoints (a) 6 May, (b) 19 June, (c) 30 June and (d) 22 July. Elevation contours derived from a SPOT DEM are also shown.

Fig. 3. Surface elevation (m) and mean winter velocities (m a⁻¹) of the four outlet glaciers Isungata Sermia, Russell Glacier, outlet A and outlet B, derived from the flowlines shown in Figure 1. Winter mean was derived from TSX data with mid-dates 11 November 2009, 7 February 2010, 9 November 2010 and 20 November 2010.

Fig. 4. Mean ice speed-up relative to winter mean of Russell Glacier’s neighboring outlets during 2010 in April, May, June, July and August. The flowlines of each of the outlets are shown in Figure 1, and the winter background is shown in Figure 3. Note the different scale in May to accommodate the considerable speed-up (∼300%) of outlet B.

To highlight spatial and temporal changes between 2009 and 2010, the flowline of Russell Glacier is split into three zones, the lower (400–600 m/0–10 km), middle (600–800 m/ 10–21 km) and upper (800–1000 m/21–34 km). In early May 2010, surface velocities in the lower portion of the ice-sheet margin doubled relative to the winter mean (Figs 2 and 5). Following the initial perturbation, velocities reduced and ice flow increased in the middle and upper marginal zones. Ice surface velocities derived from TSX data reveal a sustained (11 day average) peak of 249 m a⁻¹ around 19 June 2010, 26 km from the ice margin. Following this, ice flow decreased and returned to or dropped below the winter mean during July and August. A similar pattern occurred in 2009; however, the lack of temporal coverage of satellite observations during early summer prevents us from determining the exact time of peak flow. Owing to the paucity of data in the velocity record during May 2009, the initiation of accelerated flow (defined as ice flow 50% winter mean; Fig. 2) was derived from geodetic GPS stations located within the equivalent elevation zones (S. Doyle, unpublished data) and is in agreement with Reference BartholomewBartholomew and others (2011). Velocity results from mid- and late June 2009 capture the tail end of summer acceleration as ice flow steadily returned to winter values. Changes also occurred throughout the winter months along the Russell Glacier flowline, and mean velocities increased by 6% over the course of the 2009 winter (November 2009 to February 2010). In the upper zone, velocities in February 2010 were 22% greater at the end of winter compared to values from the start of winter the same year (November 2010), and in 2009, over the entire 57 km transect, mean velocities in March were 9% greater than in November the same year. In both years a linear trend exists of decreasing magnitude of winter speed-up with increasing elevation.

The energy-balance model reveals that, over the lower zone of the catchment, surface melting initiated on 16 May 2009 and 2 May 2010 and continued for 125 and 177 days respectively. The extended period of melt (a difference of 52 days) is due to the record-breaking temperatures during 2010 (2.7 standard deviations above the 1974–2010 mean), combined with low winter accumulation and low surface albedo (Reference Van As, Hubbard, Hasholt, Mikkelsen, Van den Broeke and FaustoVan As and others, 2012). The year 2010 is herein referred to as a ‘high melt year’, and is compared to the low melt year of 2009 (termed near-average by Reference Van As, Hubbard, Hasholt, Mikkelsen, Van den Broeke and FaustoVan As and others, 2012). To illustrate the contrasting character in ice flow between years, the results are expressed as percentage change in surface velocity relative to the winter mean, since the start of the respective melt seasons (Fig. 5). The winter mean is common to both years to allow absolute comparison, and was derived from TSX data with mid-dates 11 November 2009, 7 February 2010, 9 November 2010 and 20 November 2010.

Fig. 5. Ice speed-up relative to winter mean along the Russell Glacier transect since the start of the melt season in 2009 (blue) and 2010 (red) within three zones: lower (400–600 m/0–10 km), middle (600–800 m/10–21 km) and upper (800–1000 m/21–34 km). Owing to paucity of data in the TSX record, the initiation of the speed-up in 2009 was derived from geodetic GPS stations (S. Doyle, unpublished data) in agreement with Reference BartholomewBartholomew and others (2011). Modeled daily runoff derived from a network of automated weather stations is also shown within each zone for 2009 (light blue) and 2010 (light red), and overlapping areas are colored light purple. Calculated annual ice flux (Δv) is shown for each zone.

Accelerated ice flow occurs in 2010 just 4 days after the onset of surface melt at the ice-sheet margin (Figs 2 and 5). Ice in the lower ablation zone experiences the greatest flow increase, doubling in early May, after which the lower zone slows to winter mean levels by mid-July. Ice flow in the middle and upper zones also exhibits an initial period of accelerated flow shortly after the onset of melt, albeit of a reduced magnitude (to 50%) and up to 48 days after the onset of melt. The increase in annual ice flux (Δv) from summer ice velocity variations is expressed as the percentage by which the annual surface displacement exceeds that which would occur if winter values remained constant throughout the year (Fig. 5; Table 2). During 2010, Δv values vary from 12% in lower zones to 15% in the upper zones of the catchment. It is clear that the magnitude and extent of accelerated flow in 2010 exceeded that of 2009 in the lower and upper zones of the catchment (Fig. 5). In both 2009 and 2010, surface velocity in the lower zone of the catchment slowed and returned to overwinter levels ∼90 and ∼80 days after melt onset respectively (Table 2). In the upper zone there are further differences in the seasonal extent of accelerated flow between the years. In 2010, ice in the upper zone continued to flow 20% greater than the winter mean 114 days after the initiation of the melt season, dropping to winter values after ∼125 days. By contrast, in 2009, ice in the same elevation band had already returned to winter values within 102 days of melt onset (Fig. 5).

Table 2. Summary of seasonal characteristics between 2009 and 2010, illustrating the first day of speed-up (50% more than winter mean) after the start of the melt season, calculated annual ice flux and the number of days for velocity to return to background levels following accelerated flow. Due to paucity of velocity data at the start of the 2009 melt season, data denoted with an asterisk are based on S. Doyle (unpublished data) and Reference BartholomewBartholomew and others’ (2011) 2009 GPS data, which fall into the same elevation bands

Although flow in early May 2010 doubled relative to the winter mean, this was restricted to 12% of the 2400 km² catchment (Fig. 6), with the majority of the survey area (35%) flowing 25–50% faster than the winter mean. During June 2010 the majority of the catchment (42%) continued to flow at speeds 25–50% greater than the winter mean, with 7% and 12% of the area flowing at speeds less than the winter mean (19 and 30 June 2010 respectively). Throughout July 2010 the areal extent of ice flowing slower than the winter mean increased significantly (Table 3), and slowest velocities are obtained at the end of July when 48% of the catchment flowed at velocities less than the winter mean.

Fig. 6. Changes in the areal extent of the ice speed-up relative to the winter mean.

To highlight the differences in the magnitude and extent of accelerated flow between 2009 and 2010, mean ice speed-up relative to the winter mean along the entire Russell transect is compared with modeled daily runoff and proglacial discharge (Fig. 7). Despite the fact that in 2010 melt continued 53 days (42%) longer than during 2009, the time taken to return winter velocities is not proportional to this, with only a slight extension in the upper zone (Table 2). Furthermore the magnitude of the speed-up (mean ice flow relative to the winter mean across the entire transect) peaked at 60% in both years (Fig. 7); however, this persisted for a longer period in 2010 compared to the previous year. In both years accelerated flow does not coincide with peak proglacial discharge, which occurs in mid- to late July when ice flow is returning to winter background levels.

Fig. 7. Mean ice speed-up relative to winter mean along the Russell Glacier transect (400–1000 m/0–34 km) throughout the melt season during 2009 (bottom) and 2010 (top). Error bars represent the standard deviation of velocity covering all elevation bands. Modeled daily runoff (blue bars) (Reference Van As, Hubbard, Hasholt, Mikkelsen, Van den Broeke and FaustoVan As and others, 2012) and discharge (black line) (Reference Hasholt, Mikkelsen, Nielsen and LarsenHasholt and others, 2012) at the Watson River are also shown for both years.

Discussion

The time series of residual velocity patterns illustrate a strong seasonal trend extending up to at least 57 km (1200m elevation) from the ice margin (the limit of our data). The spatial extent of seasonal accelerated ice flow recorded here is in general agreement with previous studies (Reference Van de WalVan de Wal and others, 2008; Reference BartholomewBartholomew and others, 2011; Reference Palmer, Shepherd, Nienow and JoughinPalmer and others, 2011) and extends beyond the data limits of others (e.g. Reference Joughin, Das, King, Smith, Howat and MoonJoughin and others (2008) and Reference Sundal, Shepherd, Nienow, Hanna, Palmer and HuybrechtsSundal and others (2011), who identified seasonal patterns up to 1000 m and 20 km inland respectively). Consistent with Reference Palmer, Shepherd, Nienow and JoughinPalmer and others (2011), the velocity patterns (Fig. 2) of this study highlight spatial and temporal heterogeneity indicating that surface velocities are primarily controlled by the local-scale bed topography and the time-varying subglacial drainage network governed by surface meltwater inputs. The extent and magnitude of accelerated flow varies considerably between adjacent outlets (Figs 2 and 4) despite similar hypsometry (Fig. 3) and climatic control. This indicates that subglacial routing of water may, through the influence of basal topography, play a key role in governing individual outlet flow regime. The observed high temporal and spatial variability in ice flow in this study cautions against upscaling dynamic properties of glaciers to larger regions.

A notable feature of our results is the absolute magnitude of the velocity increase in 2010. Just 4 days after the onset of melt, 12% of the study area doubled in speed (based on 11 day averaged TSX data; Figs 2, 5 and 6). Previous studies identify summer velocity increases of 30% (Reference Van de WalVan de Wal and others, 2008), 50–100% (Reference Joughin, Das, King, Smith, Howat and MoonJoughin and others, 2008), 50–125% (Reference Sundal, Shepherd, Nienow, Hanna, Palmer and HuybrechtsSundal and others, 2011) and 48% (Reference Palmer, Shepherd, Nienow and JoughinPalmer and others, 2011). These transient accelerations combine to increase the annual ice flux. Accordingly values of Δv for 2010 in the lower zone (12%) are similar to previous records in this region (17%; Reference BartholomewBartholomew and others, 2011), and are considerably greater than in 2009 (7%). In the upper zone of the catchment, our results indicate that sustained and widespread melt in 2010 contributed more to the mean annual ice flux (15%) than in 2009: 10% in this study compared to 8% by Reference BartholomewBartholomew and others’ (2011) GPS measurements). Melt-induced acceleration involves surface meltwater accessing the subglacial system and flowing outwards under pressure, initially through an inefficient system of irregularly distributed channels or linked cavities (Reference SchoofSchoof, 2010). When subglacial water pressures approach ice-overburden pressure the ice sheet decouples from its bed, reducing basal friction, permitting enhanced flow via basal sliding. The flow perturbation propagates up-glacier, though the timing and pattern varies each year. In the upper zone the magnitude of melt-induced accelerated flow is reduced and a lag time exists as higher summer temperatures extend to increasingly higher elevations on the ice sheet. Soon after the initial melt-induced acceleration, the speed declines and ice flow in the lower zone returns to, or drops below, the winter mean ∼90 and ∼80 days after the onset of melt in 2009 and 2010 respectively. By mid-July 2010, ice in the lower marginal zone, which 2 months earlier experienced the greatest increase in velocity, flowed up to 75% slower than the winter mean (Figs 2 and 6).

The spatial extent of the deceleration observed after the initial increase in flow (Table 3) suggests a rapid evolution of the subglacial system from an inefficient distributed drainage system into an efficient one. The observed rapid return to winter velocities halfway through the melt season before peak discharge was reached (Fig. 7) supports this hypothesis and indicates that the direct correlation between melt and velocity is short-lived. Previous work suggests that in high-melt years the subglacial hydrology evolves to an efficient channelized state with such rapidity that marginal ice flux is actually reduced in comparison to low-melt years (Reference Sundal, Shepherd, Nienow, Hanna, Palmer and HuybrechtsSundal and others, 2011). Although rapid deceleration is observed for both years in the lower zone, the results here do not support this pattern. We identified that although the magnitude of the speed-up across the entire transect was similar in both years, accelerated flow persisted for a longer period in 2010 compared to 2009 (Fig. 7). This difference indicates that, in general, the extent and magnitude of melt-induced accelerated flow and therefore the calculated contribution to annual ice flux is greater during a high melt year (2010) than in a comparatively low melt year (2009) (Fig. 5; Table 2). Furthermore, the melt season of 2010 persisted 42% (53 days) longer than 2009; however, the time taken to return to winter velocities is not proportional to this, with only a slight extension in the upper zone (Table 2). This illustrates that despite the differences between the years, the magnitude and extent of accelerated flow over the entire transect during 2010 is not directly proportional to bulk meltwater runoff and discharge (Fig. 7), indicating that the relationship between meltwater availability and accelerated ice flow is short-lived and nonlinear. Numerical modeling studies have identified a dynamic switching threshold between subglacial drainage modes (Reference SchoofSchoof, 2010), with subsequent effect on ice surface velocities, recorded in previous work to be ∼1.4 cm d⁻¹ (Reference Sundal, Shepherd, Nienow, Hanna, Palmer and HuybrechtsSundal and others, 2011). Here we have presented evidence of continued accelerated ice flow during periods of melt supply well in excess of this value (Figs 5 and 7), illustrating the complexities of attempting to use simple relationships that fail to take into account the seasonal evolution of the drainage system.

Table 3. Summary of areal extent of ice flow greater or less than the winter mean. Date format is day/month/year

Previous studies illustrate that flow across the ice-sheet margin increases over the winter months (e.g. Reference Joughin, Smith, Howat, Scambos and MoonJoughin and others, 2010b). In the upper zone our data show mean velocities increased by 8 m a⁻¹ (6%) over the course of the 2009/10 winter. We also observed that surface ice flow was 22% greater at the end of winter (February 2010) compared to the start of winter the same year (November 2010), indicating increased ice flow during freezing surface conditions. Furthermore, in the current study the increase in the recorded velocity over winter was restricted to areas nearer the ice margin (850 m elevation), and little difference was observed at greater elevations, which matches the previously observed trend in this region (Reference Joughin, Smith, Howat, Scambos and MoonJoughin and others, 2010b). Work on temperate glaciers has shown a marked increase in seasonality with distance down-glacier (Reference Hooke, Calla, Holmlund, Nilsson and StroevenHooke and others, 1989), with velocities decreasing over the course of the winter. As winter approaches, the input of meltwater from the supraglacial system subsides as the surface freezes. Over winter, meltwater can be generated at the ice sheet bed through energy supplied by basal shear and geothermal heat (Reference Joughin, Das, King, Smith, Howat and MoonJoughin and others, 2008). Migrating deep pore water from the percolation zone, which persists despite cold conditions in the firn during winter, may also provide an additional source of meltwater (Reference Harper, Humphrey, Pfeffer, Brown and FettweisHarper and others, 2012; Reference Humphrey, Harper and PfefferHumphrey and others, 2012). This water becomes trapped within a distributed drainage system as basal cavities collapse, which over the winter will increase the water pressure at the ice–bed interface, permitting gradual enhancement of motion.

Conclusions

The results of this study provide snapshots of the highly variable ice flow distribution observed at the marginal zone of the GrIS extending at least 57 km up glacier. The observed seasonal flow variations are in broad agreement with previous studies (Reference Joughin, Das, King, Smith, Howat and MoonJoughin and others, 2008, Reference Joughin, Smith, Howat, Scambos and Moon2010b; Reference Bartholomew, Nienow, Mair, Hubbard, King and SoleBartholomew and others, 2010), though the large initial perturbation followed by rapid deceleration of the ice within the 2400 km² catchment indicates the basal drainage system responds rapidly to changes in supraglacial inputs. The large spatial variability of the observed velocity trends in this region, highlighted by the differences in ice flow between neighboring outlets, suggests a complex relationship between the routing of water from the surface to the bed and illustrates the importance of considering other forces such as bedrock geometry. Previous work proposes the concept that a positive feedback exists between meltwater production and ice velocity, which would lead to a more rapid response to projected future warming trends (Reference Zwally, Abdalati, Herring, Larson, Saba and SteffenZwally and others, 2002). However the hydraulic evolution of the basal drainage system acts to regulate ice flow, therefore self-limiting this positive feedback. Recent studies (Reference Sundal, Shepherd, Nienow, Hanna, Palmer and HuybrechtsSundal and others, 2011) support this hypothesis; in warmer years the period of melt-induced summer acceleration is shorter than in cooler years, as equilibrium subglacial hydrological conditions are reached earlier in the season. Modelling studies including a ‘dynamic switching’ threshold determine that above a critical rate of meltwater input, channelization and surface deceleration occurs, indicating that sustained high meltwater rates would suppress dynamic response (Reference SchoofSchoof, 2010). However, our analysis during two contrasting melt seasons does not completely concur with these findings, suggesting that a more complex relationship between meltwater availability and accelerated flow exists. The spatial and temporal heterogeneity of ice flow presented in this study suggests caution should be exercised when upscaling dynamic properties of glaciers to larger regions.