Search

Gaining knowledge of tidewater glacier (TWG) margin evolution, solid ice flux and their responses to climate over large spatio-temporal scales provides valuable context for the projection of future Greenland ice sheet (GrIS) change. Although studies of sector-wide responses of TWGs exist, studies at an ice-sheet-wide scale have only just become feasible. Here, we present a dataset of 224 annual TWG margins for 1984–2017 (n= 3801), showing that averaged over regional scales, normalised terminus change is linear. Regionally linear retreat trends were identified across most sectors of the GrIS starting in the mid-1990s, although in contrast to previous studies, the northeastern sector is shown to have experienced sustained retreat since the mid-1980s. Through cointegration analyses, individual glaciers are shown to have differing sensitivities to potential climate drivers, though on a sector-wide scale the northwest and southeast are shown to be especially sensitive to annual sea surface temperature and June–July–August air temperature, respectively. Although 92% of the analysed glaciers experience retreat across the GrIS, observed increases in absolute flux for the entire ice sheet can be explained by changes in just 11 of these TWGs.

The retention of meltwater in the accumulation area of the Greenland ice sheet and other Arctic ice masses buffers their contribution to sea level change. However, sustained warming also results in impermeable ice layers or ‘ice slabs’ that seal the underlying pore space. Here, we use a 1-D, physically based, high-resolution model to simulate the surface mass balance (SMB), percolation, refreezing, ice layer formation and runoff from across the high-elevation area of Devon Ice Cap, Canada, from 2001 to 2016. We vary the thickness of the ‘impermeable’ ice layer at which underlying firn becomes inaccessible to meltwater. Thick near-surface ice layers are established by an initial deep percolation, the formation of decimetre ice layers and the infilling of interleaving pore space. The cumulative SMB increases by 48% by varying impermeable layer thickness between 0.01 and 5 m. Within this range we identify narrower range (0.25–1 m) that can simulate both the temporal variability in SMB and the observed near-surface density structure. Across this range, cumulative SMB variation is limited to 6% and 45–49% of mass retention takes place within the annually replenished snowpack. Our results indicate cooler summers after intense mid-2000s warming have led to a partial replenishment of pore space.

Understanding the controls on the amount of surface meltwater that refreezes, rather than becoming runoff, over polar ice masses is necessary for modeling their surface mass balance and ultimately for predicting their future contributions to global sea level change. We present a modified version of a physically based model that includes an energy balance routine and explicit calculation of near-surface meltwater refreezing capacity, to simulate the evolution of near-surface density and temperature profiles across Devon Ice Cap in Arctic Canada. Uniquely, our model is initiated and calibrated using high spatial resolution measurements of snow and firn densities across almost the entire elevation range of the ice cap for the summer of 2004 and subsequently validated with the same type of measurements obtained during the very different meteorological conditions of summer 2006. The model captures the spatial variability across the transect in bulk snowpack properties although it slightly underestimates the flow of meltwater into the firn of previous years. The percentage of meltwater that becomes runoff is similar in both years; however, the spatial pattern of this melt-runoff relationship is different in the 2 years. The model is found to be insensitive to variation in the depth of impermeable layers within the firn but is very sensitive to variation in air temperature, since the refreezing capacity of firn decreases with increasing temperature. We highlight that the sensitivity of the ice cap's surface mass balance to air temperature is itself dependent on air temperature.

The flow of glaciers and ice streams is strongly influenced by the presence of water at the interface between ice and bed. In this paper, a hydrological model evaluating the subglacial water pressure is developed with the final aim of estimating the sliding velocities of glaciers. The global model fully couples the subglacial hydrology and the ice dynamics through a water-dependent friction law. The hydrological part of the model follows a double continuum approach which relies on the use of porous layers to compute water heads in inefficient and efficient drainage systems. This method has the advantage of a relatively low computational cost that would allow its application to large ice bodies such as Greenland or Antarctica ice streams. The hydrological model has been implemented in the finite element code Elmer/Ice, which simultaneously computes the ice flow. Herein, we present an application to the Haut Glacier d'Arolla for which we have a large number of observations, making it well suited to the purpose of validating both the hydrology and ice flow model components. The selection of hydrological, under-determined parameters from a wide range of values is guided by comparison of the model results with available glacier observations. Once this selection has been performed, the coupling between subglacial hydrology and ice dynamics is undertaken throughout a melt season. Results indicate that this new modelling approach for subglacial hydrology is able to reproduce the broad temporal and spatial patterns of the observed subglacial hydrological system. Furthermore, the coupling with the ice dynamics shows good agreement with the observed spring speed-up.

Several different methodologies have previously been employed in the tracking of glacier terminus change, though a systematic comparison of these has not been undertaken. The frequent application of single methods to multiple glaciers over large geographical areas such as Greenland, raises the question of whether individual methodologies are robust. In this study we evaluate three existing methodologies that have been widely used to track terminus change (the centre-line, bow and box methods) against a full range of idealized glaciological scenarios and six examples of real glaciers. We also evaluate two new methodologies that aim to reduce measurement error compared with the existing methodologies. The first is a modification to the box method that can account for termini retreating through fjords that change orientation (termed the curvilinear box method), while the second determines the average terminus position relative to the glacier centre line using an inverse distance weighting extrapolation (termed the extrapolated centre-line method). No single method tested achieved complete accuracy for all scenarios, though the extrapolated centre-line method was able to successfully account for variable fjord orientation, width and terminus geometry with the least error.

[1] The influence of meltwater on the dynamics and geomorphic impact of the Greenland Ice Sheet is strongly controlled by the morphology of the ice sheet's drainage system. However, this system and its evolution through the melt season remain poorly understood. Here we present the results of an intensive programme of dye tracing experiments undertaken along the lower 14 km of a land-terminating Greenlandic outlet glacier over a period of four months during the 2010 melt season. These data are interpreted in conjunction with observations of proglacial discharge, englacial water storage, surface melt rates and ice velocity to produce a detailed picture of the changing hydrology of the glacier. Following the onset of melt in the spring, inputs to the drainage system regularly exceed outputs, causing the englacial water level to rise to the ice sheet surface. During this time there is a rapid transition from distributed to channelized drainage in those parts of the drainage system closed by ice deformation over winter. As the melt season progresses, channel efficiency increases and englacial storage and ice velocity decrease. High-velocity events continue to be observed following the channelization of the drainage system however, indicating that hydrological forcing of ice velocity occurs despite the existence of channels during periods when meltwater inputs exceed the capacity of the subglacial drainage system.

The transfer of surface-generated meltwater to the subglacial drainage system through full ice thickness crevassing may lead to accelerated glacier velocities, with implications for ice motion under future climatic scenarios. Accurate predictions of where surface meltwater accesses the ice/bed interface are therefore needed in fully coupled hydrodynamic ice-sheet models. We present a spatially distributed modelling routine for predicting the location and timing of delivery of surface-derived meltwater to the ice/bed interface through moulins and supraglacial lake drainage. The model is explained as it is applied to the Croker Bay glacial catchment of Devon Ice Cap, Canada. The formation of moulins, drainage of lakes, and the transfer of meltwater through the full ice thickness are modelled for the 2004 and 2006 ablation seasons. Through this case study we assess the model’s sensitivity to degree-day factors, fracture toughness, tensile strength and crevasse width, and confirm that parameters influencing the rate at which water fills a crevasse are the most significant controls on the ability of a crevasse to reach the bed. Increased surface melt production, therefore, has the potential to significantly influence the spatial and temporal transfer of meltwater through surface-to-bed connections in a warmer climate

Mair, D., Hubbard, B., Nienow, P., Willis, I., Fischer, U. H. (2008). Diurnal fluctuations in glacier ice deformation: Haut Glacier d'Arolla, Switzerland. Earth Surface Processes and Landforms, 33(8), 1272-1284. Sponsorship: Carnegie Trust for the Universities of Scotland

Interpretation of ice mass elevation changes observed by satellite altimetry demands quantification of the proportion of elevation change which is attributable to variations in firn densification. Detailed stratigraphic logging of snowpack structure and density was carried out at ~1km intervals along a 47 km transect on Devon Ice Cap, Canada, in spring (pre-melt) and autumn (during/ after melt) 2004 and 2006 to characterize seasonal snowpack variability across the full range of snow facies. Simultaneous meteorological measurements were gathered. Spring (pre-melt) snowpacks show low variability over large spatial scales, with low-magnitude changes in density. The end-of-summer/ autumn density profiles show high variability in both 2004 and 2006, with vastly different melt regimes generating dissimilar patterns of ice-layer formation over the two melt seasons. Dye-tracing experiments from spring to autumn 2006 reveal that vertical and horizontal distribution of meltwater flow within and below the annual snowpack is strongly affected by the pre-existing, often subtle stratigraphic interfaces in the snowpack, rather than its bulk properties. Strong interannual variability suggests that using a simple relationship between air temperature, elevation and snowpack densification to derive mass change from measurements of elevation change across High Arctic ice caps may be misguided. Melt timing and duration are important extrinsic factors governing snowpack densification and ice-layer formation in summer, rather than averaged air temperatures.

Meltwater delivered to the bed of the Greenland Ice Sheet is a driver of variable ice-motion through changes in effective pressure and enhanced basal lubrication. Ice surface velocities have been shown to respond rapidly both to meltwater production at the surface and to drainage of supraglacial lakes, suggesting efficient transfer of meltwater from the supraglacial to subglacial hydrological systems. Although considerable effort is currently being directed towards improved modelling of the controlling surface and basal processes, modelling the temporal and spatial evolution of the transfer of melt to the bed has received less attention. Here we present the results of spatially distributed modelling for prediction of moulins and lake drainages on the Leverett Glacier in Southwest Greenland. The model is run for the 2009 and 2010 ablation seasons, and for future increased melt scenarios. The temporal pattern of modelled lake drainages are qualitatively comparable with those documented from analyses of repeat satellite imagery. The modelled timings and locations of delivery of meltwater to the bed also match well with observed temporal and spatial patterns of ice surface speed-ups. This is particularly true for the lower catchment (< 1000 m a.s.l.) where both the model and observations indicate that the development of moulins is the main mechanism for the transfer of surface meltwater to the bed. At higher elevations (e.g. 1250-1500 m a.s.l.) the development and drainage of supraglacial lakes becomes increasingly important. At these higher elevations, the delay between modelled melt generation and subsequent delivery of melt to the bed matches the observed delay between the peak air temperatures and subsequent velocity speed-ups, while the instantaneous transfer of melt to the bed in a control simulation does not. Although both moulins and lake drainages are predicted to increase in number for future warmer climate scenarios, the lake drainages play an increasingly important role in both expanding the area over which melt accesses the bed and in enabling a greater proportion of surface melt to reach the bed., Postprints der Universität Potsdam : Mathematisch-Naturwissenschaftliche Reihe, 513

This paper describes the structural glaciology of Haut Glacier d'Arolla, a small valley glacier fed by two distinct accumulation basins in the Swiss Alps. A considerable body of field data is presented alongside observations from ground and aerial photographs. Suites of structures identified in the field and from aerial photographs are first described in nongenetic terms before being assigned regular structural terms. Haut Glacier d'Arolla is dominated by primary stratification, which is progressively folded and eventually transposed into longitudinal foliation as it moves into the glacier tongue. Crevasses and crevasse traces cross-cut and in places displace primary stratification and longitudinal foliation. Crevasse traces are formed by the closure of crevasses or may represent tensional veins. On their journey downglacier, crevasse traces become increasingly rotated. Close to the snout, some crevasse traces become reactivated as thrust faults. Strain ellipses, derived from the velocity field, ...

The surface velocity field of Devon Ice Cap, Nunavut, Canada, was mapped using interferometric synthetic aperture radar (InSAR). Ascending European Remote-sensing Satellite 1 and 2 (ERS-1/-2) tandem mode data were used for the western and southeast sectors, and 3 day repeat pass ERS-1 imagery for the northeast sector. Speckle-tracking procedures were used with RADARSAT 1 imagery to obtain surface velocities over the terminus of Belcher Glacier (a major calving front) where decorrelation between ERS data occurred. The InSAR data highlight a significant contrast in ice-flow dynamics between the east and west sides of the ice cap. Ice movement west of the main north–south divide is dominated by relatively uniform ‘sheet’ flow, but three fast-flowing outlet glaciers that extend 14–23km beyond the ice-cap margin also drain this region. Several outlet glaciers that extend up to 60 km inland from the eastern margin drain the eastern side of the ice cap. The dominant ice-flow regimes were classified based on the relationship between the driving stress (averaged over a length scale of ten ice thicknesses) and the ratio of surface velocity to ice thickness. The mapped distribution of flow regimes appears to depict the spatial extent of basal sliding across the ice cap. This is supported by a close relationship between the occurrence of flow stripes on the ice surface and flow regimes where basal sliding was found to be an important component of the glacier motion. Iceberg calving rates were computed using measured surface velocities and ice thicknesses derived from airborne radio-echo sounding. The volume of ice calved between 1960 and 1999 was estimated to be 20.5 ± 4.7 km3 (or 0.57 km3 a–1). Approximately 89% of this loss occurred along the eastern margin. The largest single source is Belcher Glacier, which accounts for ~50% of the total amount of ice calved.

High rates of surface uplift and horizontal velocities were measured during a hydrologically induced spring speed-up event. Spatial patterns of surface uplift are analyzed to estimate components of vertical motion due to flow along an inclined bed and vertical strain. Areas are identified where surface uplift was most likely due in part to the opening or enlargement of subglacial cavities by bed separation. Results suggest a widespread enlargement of subglacial cavities during the event, and survival of residual cavities after the event. The spatial pattern of cavity enlargement closely matches previously identified axes of preferential subglacial drainage. It is suggested that localized cavity opening along axes of preferential drainage may constitute the initial stage in the seasonal development of channelized subglacial drainage. It is concluded that spatial and temporal variations in glacier motion may play an active role in determining the structure and rate of development of subglacial drainage during the summer melt season.

Many tidewater glaciers in Greenland are known to have undergone significant retreat during the last century following their Little Ice Age maxima. Where it is possible to reconstruct glacier change over this period, they provide excellent records for comparison to climate records, as well as calibration/validation for numerical models. These glacier change records therefore allow for tests of numerical models that seek to simulate tidewater glacier behaviour over multi-decadal to centennial timescales. Here we present a detailed record of behaviour from Kangiata Nunaata Sermia (KNS), SW Greenland, between 1859 and 2012, and compare it against available oceanographic and atmospheric temperature data between 1871 and 2012. We also use these records to evaluate the ability of a well-established one-dimensional flow-band model to replicate behaviour for the observation period. The record of terminus change demonstrates that KNS has advanced/retreated in phase with atmosphere and ocean climate anomalies averaged over multi-annual to decadal timescales. Results from an ensemble of model runs demonstrate that observed dynamics can be replicated. Model runs that provide a reasonable match to observations always require a significant atmospheric forcing component, but do not necessarily require an oceanic forcing component. Although the importance of oceanic forcing cannot be discounted, these results demonstrate that changes in atmospheric forcing are likely to be a primary driver of the terminus fluctuations of KNS from 1859 to 2012. We propose that the detail and length of the record presented makes KNS an ideal site for model validation exercises investigating links between climate, calving rates, and tidewater glacier dynamics.

Tidewater glaciers in Greenland experienced widespread retreat during the last century. Information on their behaviour prior to this is often poorly constrained due to lack of observations, while determining the drivers prior to instrumental records is also problematic. Here we present a record of the dynamics of Kangiata Nunaata Sermia (KNS), southwest Greenland, from its Little Ice Age maximum (LIAmax) to 1859 – the period before continuous air temperature observations began at Nuuk in 1866. Using glacial geomorphology, historical accounts, photographs and GIS analyses, we provide evidence KNS was at its LIAmax by 1761, had retreated by ~5 km by 1808 and a further 7 km by 1859. This predates retreat at Jakobshavn Isbræ by 43–113 years, demonstrating the asynchroneity of tidewater glacier terminus response following the LIA. We use a one-dimensional flowband model to determine the relative sensitivity of KNS to atmospheric and oceanic climate forcing. Results demonstrate that terminus forcing rather than surface mass balance drove the retreat. Modelled glacier sensitivity to submarine melt rates is also insufficient to explain the retreat observed. However, moderate increases in crevasse water depth, driving an increase in calving, are capable of causing terminus retreat of the observed magnitude and timing.

The surface motion of Haut Glacier d’Arolla, Switzerland, was monitored at a high spatial and temporal resolution. Data are analyzed to calculate surface velocities, surface strain rates and the components of the glacier force budget before, during and after an early melt season speed-up or “spring event”. We investigate the extent to which variations in glacier motion can be attributed to hydrologically induced local forcing or to non-local forcing transmitted via horizontal stress gradients. Enhanced glacier motion is dependent on a change in the spatial distribution of areas of high drag across the glacier.

An indirect methodology for determining the distribution of mass balance at high spatial resolution using remote sensing and ice-flow modelling is presented. The method, based on the mass-continuity equation, requires two datasets collected over the desired monitoring interval: (i) the spatial pattern of glacier surface-elevation change, and (ii) the mass-flux divergence field. At Haut Glacier d’Arolla, Valais, Switzerland, the mass-balance distribution between September 1992 and September 1993 is calculated at 20 m resolution from the difference between the pattern of surface-elevation change derived from analytical photogrammetry and the mass-flux divergence field determined from three-dimensional, numerical flow modelling constrained by surface-velocity measurements. The resultant pattern of mass balance is almost totally negative, showing a strong dependence on elevation, but with large localized departures. The computed distribution of mass balance compares well (R2 = 0.91) with mass-balance measurements made at stakes installed along the glacier centre line over the same period. Despite the highly optimized nature of the flow-modelling effort employed in this study, the good agreement indicates the potential this method has as a strategy for deriving high spatial and temporal-resolution estimates of mass balance.

A three-dimensional, finite-difference model based on a first-order solution of the ice-flow equations is applied to Haut Glacier d’Arolla, Switzerland. The numerical model successfully converges at horizontal resolutions down to 70 m, so a number of detailed comparisons with field data can be made. Modelled surface velocities with no basal sliding component are compared with surface velocities observed on the glacier over four different time periods. The best fit is achieved with over-winter surface velocities (R2 = 0.75) using a rate factor, A, in Glen’s flow law of 0.063 a−1 bar−3. Surface zones of maximum computed effective stress display a high level of coincidence with observed crevassing, the orientation of which is successfully predicted by the direction of the tensile component of the computed principal surface stress. Comparison of the relative magnitude and direction of computed principal stresses with principal strains measured at the ice surface also correspond closely. In an attempt to simulate the observed annual velocity distribution within a cross-section of the glacier tongue, we incorporate two basal-motion patterns into the model. By treating net annual ice motion as a time-weighted composite of three separate flow situations: normal sliding, enhanced sliding and no sliding, we are able to reproduce the key features of the observed cross-sectional ice and basal slip velocity distributions. These experiments indicate there may be substantial decoupling taking place along an elongated narrow zone at the bed of Haut Glacier d’Arolla and that this decoupling interacts in a complex manner with the englacial stress and strain field.

[1] High resolution measurements of ice motion along a ∼120 km transect in a land-terminating section of the GrIS reveal short-term velocity variations (

We present subdaily ice flow measurements at four GPS sites between 36 and 72 km from the margin of a marine‐terminating Greenland outlet glacier spanning the 2009 melt season. Our data show that >35 km from the margin, seasonal and shorter–time scale ice flow variations are controlled by surface melt–induced changes in subglacial hydrology. Following the onset of melting at each site, ice motion increased above background for up to 2 months with resultant up‐glacier migration of both the onset and peak of acceleration. Later in our survey, ice flow at all sites decreased to below background. Multiple 1 to 15 day speedups increased ice motion by up to 40% above background. These events were typically accompanied by uplift and coincided with enhanced surface melt or lake drainage. Our results indicate that the subglacial drainage system evolved through the season with efficient drainage extending to at least 48 km inland during the melt season. While we can explain our observations with reference to evolution of the glacier drainage system, the net effect of the summer speed variations on annual motion is small (∼1%). This, in part, is because the speedups are compensated for by slowdowns beneath background associated with the establishment of an efficient subglacial drainage system. In addition, the speedups are less pronounced in comparison to land‐terminating systems. Our results reveal similarities between the inland ice flow response of Greenland marine‐ and land‐terminating outlet glaciers.

We present estimates of accumulation rate along a 200 km transect ranging in elevation from 2750 to 3150 m in the dry snow zone on the western slope of the Greenland Ice Sheet. An airborne radar altimeter is used to estimate the thickness of annual internal layers and, in conjunction with ground based snow/firn density profiles, annual accumulation rates between 1998 and 2003 are derived. A clear gradient in the thickness of each layer observed by the radar altimeter and in the associated estimates of annual accumulation is seen along the transect, with a 33.6% ± 16% mean decrease in accumulation from west to east. The observed inter-annual variability is high, with the annual mean accumulation rate estimated at 0.359 m.w.e. yr−1 (s.d. ± 0.049 m.w.e. yr−1). Mean accumulation rates modelled using meteorological models overestimate our results by 16% on average, but by 32% and 42% in the years 2001 and 2002. The methodology presented here demonstrates the potential to obtain accurate and spatially extensive accumulation rates from radar altimeters in regions of ice sheets where field observations are sparse, and accumulation rates greater than several tens of cm.

[1] This paper estimates the mass balance of the Prince of Wales Icefield, Ellesmere Island, Canada, averaged over four decades, from measurements of surface mass balance (SMB) and iceberg calving. Shallow ice core net accumulation measurements and annual mass balance stake measurements are used in conjunction with a digital elevation model and knowledge of the location of the dominant moisture source for precipitation over the ice cap to interpolate and extrapolate spatial patterns of SMB across the Prince of Wales Icefield. The contribution of iceberg calving to the mass balance is calculated from estimates of (1) the annual volume of ice discharged at the major tidewater glacier termini and (2) the annual volume loss or gain due to terminus fluctuations. Two different approaches to determining the SMB conclude that the SMB of the ice field is approximately in balance (average equals −0.1 ± 0.4 km3 w.e. a−1, where w.e. means water equivalent) largely because of its proximity to the main year-round moisture source that is the Smith Sound portion of the North Open Water polynya. Iceberg calving is a highly significant component of mass loss (−1.9 ± 0.2 km3 w.e. a−1) and is sufficient to make the overall mass balance of the ice field averaged over the period 1963–2003 clearly negative (−2 ± 0.45 km3 w.e. a−1, equivalent to a mean-specific mass balance across the ice field of −0.1 m w.e. a−1). The Prince of Wales Icefield contributes ∼0.005 mm a−1 to global eustatic sea level rise.

With global land and sea temperatures rising, the importance of accurate monitoring of the world's ice sheets is increasing. Satellite radar altimetry can be used to measure elevation changes of ice sheets from which mass balance can be derived. In the percolation zone of ice sheets, summer melt which percolates into the snowpack and refreezes causes a re-distribution of mass through densification, which can result in elevation changes which are not representative of changes in mass balance. We present data collected in the percolation zone of the Greenland Ice Sheet prior to and post the processes of summer melt, percolation and refreezing. Data from nine sites along two 1-km transects show that in 2004 there was a 31.6% increase in accumulation over the summer, but due to surface melting percolation and refreezing the average snowpack density increased by 26.2%, resulting in only a 5.3% increase in elevation. Our results indicate that in areas of substantial seasonal melt and refreezing, densification rates must be accurately quantified if mass balance estimates are to be usefully derived from surface elevation change.

The mass balance of polythermal ice masses is critically dependent on the proportion of surface-generated meltwater that subsequently refreezes in the snowpack and firn. In order to quantify this effect and to characterize its spatial variability, we measured near-surface (26%, resulting in a 32% increase in net accumulation. This ‘seasonal densification’ increased at lower elevations, rising to 47% 10 km closer to the ice-sheet margin at 1860ma.s.l. Density/depth profiles from nine sites within 1 km2 at ~1945ma.s.l. reveal complex stratigraphies that change over short spatial scales and seasonally. We conclude that estimates of mass-balance change cannot be calculated solely from observed changes in surface elevation, but that near-surface densification must also be considered. However, predicting spatial and temporal variations in densification may not be straightforward. Further, the development of complex firn-density profiles both masks discernible annual layers in the near-surface firn and ice stratigraphy and is likely to introduce error into radar-derived estimates of surface elevation.

[1] This paper uses a combination of field data and three-dimensional modeling to investigate the spatial variability in basal conditions required to induce observed fluctuations in diurnal ice velocity at Haut Glacier d'Arolla, Switzerland. A network of surface velocity markers was observed at intervals of as little as four hours over diurnal cycles in both winter and late summer. Winter motion showed limited diurnal variability, presumably due to the absence of supraglacial meltwater inputs. By contrast, diurnal fluctuations in ice motion were recorded in summer across the lower and upper glacier. In the lower glacier, surface velocities were intimately linked to hydrological forcing in the vicinity of a subglacial channel. Previously observed diurnal excursions of meltwater away from the channel should reduce areas of basal drag adjacent to the channel thereby impacting on ice dynamics. Using a first-order ice flow approximation, we investigated the distribution of basal shear traction adjacent to the channel necessary to replicate the observed surface velocity field during periods of rapid ice motion. The modeling suggests that the observed variations in diurnal velocity will only occur with extensive reductions in basal drag across a transverse zone of up to 560 m across, well beyond the immediate vicinity and previously observed extent of diurnal excursions of meltwater away from the subglacial channel.

[1] In April–May 2000, eight boreholes were drilled to ∼15–20 m depth on the Devon Ice Cap. 137Cs γ activity profiles of each borehole showed a peak count rate at depth that is associated with fallout from atmospheric thermonuclear weapons testing in 1963. Snow, firn, and ice densities were measured at each core site and were used to estimate the average pattern of mass balance across the accumulation zone of the ice cap over the period 1963–2000. The average mass balance across the entire ice cap for the period 1960–2000 was also estimated using a degree-day model driven by data derived from on-ice temperature sensors and long-term measurements at Resolute Bay. Best fitting degree-day factors were determined for different sectors of the ice cap by comparing model output with repeated annual mass balance measurements made along two transects (Koerner, 1970). The results suggest that the ice cap has lost ∼1.6 km3 water per year, equivalent to a mean net mass balance of approximately −0.13 m We a−1. Estimates of the mean mass balance for individual drainage basins reveal regions of positive and negative mass balance that are consistent with remotely sensed observations of advancing and retreating ice cap margins, respectively.

[1] The relationship between the evolution of subglacial drainage system morphology and spatial patterns of glacier surface velocity was investigated using dye tracing experiments and ground surveying throughout the 1995 melt season at Haut Glacier d'Arolla, Switzerland. With the onset of high and variable melt season discharges, subglacial drainage changed from a predominantly distributed system to a predominantly channelized system. The change occurred later farther up glacier. During the period of drainage evolution the glacier was subjected to three periods of rapidly rising meltwater discharge. The magnitude and spatial pattern of the glacier's velocity response differed between these periods and can be explained in terms of the impact of the evolving drainage system morphology on the amplitude and spatial distribution of basal hydrological forcing. Increasing discharge through a distributed drainage system caused widespread basal forcing and high glacier velocity. Increasing discharge through incipient channels below moulins, not yet connected to the main channel, caused more localized basal forcing and slightly increased glacier velocity. Increasing discharge through a fully channelized drainage system caused no significant basal forcing and glacier velocity was not significantly different from the annual deformation flow pattern. Empirical orthogonal function analysis of flow patterns defined two distinct spatial modes of surface velocity which corresponded closely with the drainage system morphologies inferred to be present during each event. The relative importance of these modes changed through the melt season, suggesting a temporal change in the spatial pattern of hydrologically induced basal forcing.

Three early-melt-season high-velocity events (or “spring events”) occurred on Haut Glacier d’Arolla, Switzerland, during the melt seasons of 1998 and 1999. The events involve enhanced glacier velocity during periods of rapidly increasing bulk discharge in the proglacial stream and high subglacial water pressures. However, differences in spatial patterns of surface velocity, internal ice deformation rates, the spatial extent of high subglacial water pressures and in rates of subglacial sediment deformation suggest different hydrological and mechanical controls. The data from two of the events suggest widespread ice–bed decoupling, particularly along a subglacial drainage axis creating the highest rates of basal motion and “plug flow” in the overlying ice. The other event showed evidence of less extensive ice–bed decoupling and sliding along the drainage axis with more mechanical support for ice overburden transferred to areas adjacent to decoupled areas. We suggest that: (1) plug flow may be a common feature on glaciers experiencing locally induced reductions in basal drag; (2) under certain circumstances, enhanced surface motion may be due in part to non-locally forced enhanced bed deformation; and (3) subglacial sediment deformation is confined to a depth of the order of centimetres to decimetres.

Records of surface motion, englacial tilt and repeat inclinometry are used to determine patterns of surface, internal and basal motion across the tongue of Haut Glacier d’Arolla, Switzerland, over temporal scales ranging from days to months. Findings are interpreted with reference to contemporaneous measurements of subglacial water pressures, and prior knowledge of the glacier’s subglacial drainage-system structure. Long-term inclinometry results show pronounced extrusion flow over a subglacial drainage axis, with basal velocities up to twice those measured at the glacier surface. Deformation profiles are more conventional away from the drainage axis, with basal velocities ∼60–70% of surface velocities. Comparison of long-term tilt rates from repeat inclinometry and englacial tiltmeters shows close correspondence. Englacial tiltmeter data are used to reconstruct internal velocity profiles and to split surface velocities into internal deformation and basal motion contributions over spring, summer and autumn/winter periods. Although, spatial patterns of surface movement are similar between periods, patterns of internal and basal motion are not. Results are interpreted in terms of the location of sticky and slippery spots, with temporally changing patterns of basal drag reflecting changing distributions of water pressure.