Search

The increasing volume and spatio-temporal resolution of satellite-derived ice velocity data has created new exploratory opportunities for the quantitative analysis of glacier dynamics. One potential technique, Proper Orthogonal Decomposition (POD), also known as Empirical Orthogonal Functions, has proven to be a powerful and flexible technique for revealing coherent structures in a wide variety of environmental flows. In this study we investigate the applicability of POD to an openly available TanDEM-X/TerraSAR-X derived ice velocity dataset from Sermeq Kujalleq (Jakobshavn Isbræ), Greenland. We find three dominant modes with annual periodicity that we argue are explained by glaciological processes. Mode 1 is interpreted as relating to the stress-reconfiguration at the glacier terminus, known to be an important control on the glacier’s dynamics. Modes 2 and 3 together relate to the development of the spatially heterogenous glacier hydrological system and are primarily driven by the pressurisation and efficiency of the subglacial hydrological system. During the melt season, variations in the velocity shown in Modes 2 and 3 are explained by the drainage of nearby supraglacial melt ponds, as identified with a Google Earth Engine MODIS dynamic thresholding technique. By isolating statistical structures within velocity datasets, and through their comparison to glaciological theory and complementary datasets POD indicates which glaciological processes are responsible for the changing bulk velocity signal, as observed from space. With the proliferation of optical and radar derived velocity products (e.g. MEaSUREs/ESA CCI/PROMICE) we suggest POD, and potentially other modal decomposition techniques, will become increasingly useful in future studies of ice dynamics.

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 increasing volume and spatio-temporal resolution of satellite-derived ice velocity data has created new exploratory opportunities for the quantitative analysis of glacier dynamics. One potential technique, Proper Orthogonal Decomposition (POD), also known as Empirical Orthogonal Functions, has proven to be a powerful and flexible technique for revealing coherent structures in a wide variety of environmental flows: mapping hydraulic vortex shedding patterns, the dynamics of fluidised granular beds, and the magnetohydrodynamics of sunspots.POD exactly describes a series of snapshots from a flow field with the product of ranked spatially orthogonal Eigenfunctions, or “modes” of spatial weighting, and one-dimensional “temporal” coefficients (Eigenvectors). In many cases the variance of the flow field is well described by just a few dominant modes. The orthogonal nature of each mode, by definition, means that the relative contribution of independent forcing mechanisms on the flow can, in theory, be separated.In this study we investigate the applicability of POD to freely available TanDEM-X/TerraSAR-X derived ice velocity datasets of Sermeq Kujalleq (Jakobshavn Glacier), Greenland. We outline the POD procedure using the singular value decomposition of a rearranged and resampled velocity matrix and investigate the factors responsible for the dominant modes. We find dominant modes interpreted as relating to the stress-reconfiguration at the glacier terminus and the development of the glacier hydrological system, but also find that the POD is sensitive to data resampling and quality. With the proliferation of publicly available optical and radar derived velocity products (e.g. MEaSUREs/ESA CCI) we suggest POD, and potentially other modal decomposition techniques, will become increasingly useful in future studies of ice dynamics.

The Greenland Ice Sheet (GrIS) is known to have experienced widespread retreat over the last century. Information on outlet glacier dynamics, prior to this, are limited due to both a lack of observations and a paucity of mapped or mappable deglacial evidence which restricts our understanding of centennial to millennial timescale dynamics of the GrIS. Here we present glacial geomorphological mapping, for upper Godthåbsfjord, covering 5800 km2 at a scale of 1:92,000, using a combination of ASTER GDEM V2, a medium-resolution DEM (error

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.

This progress report is the first of three overviews of important fields of study within contemporary glaciology, illustrated with reference to exemplar research contributions and highlighting key developments in each field since 2005, the date of the last review of ‘Glaciology’ in this journal. The topics covered in this paper are: (1) the mass balance of glaciers and ice sheets and their contributions to global sea level change; and (2) atmospheric melt-induced influences on the dynamics of the Greenland Ice Sheet. The paper aims to provide an overview of key developments over the last six years from each topic, with discussion of the continuing debates and unresolved issues which constitute the foci for ongoing and future research.

An acceleration of ice-mass loss has been observed near the margin of the Greenland ice sheet, partly as a result of faster ice motion. Observations by GPS receivers reveal high seasonal variability in ice motion, with summer motion up to 220% higher than winter background levels. The Greenland ice sheet contains enough water to raise sea levels by 7 m. However, its present mass balance and future contribution to sea level rise is poorly understood1. Accelerated mass loss has been observed near the ice sheet margin, partly as a result of faster ice motion2,3,4. Surface melt waters can reach the base of the ice sheet and enhance basal ice motion5,6. However, the response of ice motion to seasonal variations in meltwater supply is poorly constrained both in space and time. Here we present ice motion data obtained with global positioning system receivers located along a ∼35 km transect at the western margin of the Greenland ice sheet throughout a summer melt season. Our measurements reveal substantial increases in ice velocity during summer, up to 220% above winter background values. These speed-up events migrate up the glacier over the course of the summer. The relationship between melt and ice motion varies both at each site throughout the melt season and between sites. We suggest that these patterns can be explained by the seasonal evolution of the subglacial drainage system similar to hydraulic forcing mechanisms for ice dynamics that have been observed at smaller glaciers.

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

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

Statistical classification of hydrograph form is used to elucidate the controls on diurnal runoff cycle evolution at Haut Glacier d'Arolla, Switzerland during the 1998 and 1999 melt seasons. Hydrographs are objectively grouped using statistical techniques into four principal types that are qualitatively interpreted as rising, falling, peaked-falling and peaked hydrographs. Peaked hydrographs are further grouped on the basis of the magnitude of their bulk flow, baseflow and diurnal flow components. Comparison with the evolution of meltwater sources and pathways demonstrates that runoff cycles evolve systematically during the melt season in response to removal of the seasonal snowpack from the ablation area. Peaked hydrographs predominate following the onset of snowpack removal and demonstrate an increasing and progressively earlier diurnal peak, but also an unusually low baseflow component that is probably due to surface melt mainly contributing direct to subglacial channels. Runoff cycle evolution has potentially significant geomorphic implications because peaked surface runoff cycles result in the formation of hydraulically efficient, channelised subglacial drainage and a significant increase in the gradient of the relationship between suspended sediment transport and discharge. Increasingly peaked diurnal cycles also result in increased basal sediment availability, most likely related to high diurnal water pressure variation within subglacial channels that may also have enhanced rates of basal sliding and hence subglacial erosion. Differences in runoff cycle form and evolution therefore have the potential to significantly influence glacial erosion rates and sediment yields.

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, ...

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.

The flow of glaciers and ice-streams is strongly influenced by the presence of water at the interface between ice and bedrock. 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.

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.

[1] We measure hydrological parameters in meltwater draining from an outlet glacier in west Greenland to investigate seasonal changes in the structure and behaviour of the hydrological system of a large catchment in the Greenland ice sheet (GrIS). Our data reveal seasonal upglacier expansion and increase in hydraulic efficiency of the subglacial drainage system, across a catchment >600 km2, to distances >50 km from the ice-sheet margin. This expansion occurs episodically in response to the drainage of surface meltwaters into a hitherto inefficient subglacial drainage system as new input locations become active progressively further upglacier; this system is similar to Alpine glaciers. These observations provide the first synopsis of seasonal hydrological behaviour in the ablation zone of the GrIS.

[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] 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.

Substrate properties and hydrological conditions at the base of Haut Glacier d'Arolla, Switzerland, are nonuniform. The thickness and grain size of subglacial sediment vary on length scales less than one ice thickness, whereas hydrological conditions vary seasonally and on length scales of about one ice thickness transverse to ice motion. There is a close relationship between the annually averaged velocity field and this nonuniformity of bed conditions. Basal motion dominates in an area with large water-level variations in boreholes, whereas ice deformation contributes significantly to total movement elsewhere. Localized enhanced basal motion occurs primarily in summer, especially during a July “spring event,” but this motion is barely discernible in annually averaged surface velocity measurements, because transverse coupling suppresses surface expression of the basal motion discontinuity. These results highlight the need to include representations of bed nonuniformity in models of glacier flow and to consider ice deformation and basal motion as interdependent processes.

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.

The influence of supraglacial debris on the rate and spatial distribution of glacier surface melt is well established, but its potential impact on the structure and evolution of the drainage system of extensively debris-covered glaciers has not been previously investigated. Forty-eight dye injections were conducted on Miage Glacier, Italian Alps, throughout the 2010 and 2011 ablation seasons. An efficient conduit system emanates from moulins in the mid-part of the glacier, which are downstream of a high melt area of dirty ice and patchy debris. High melt rates and runoff concentration by intermoraine troughs encourages the early-season development of a channelized system downstream of this area. Conversely, the drainage system beneath the continuously debris-covered lower ablation area is generally inefficient, with multi-peaked traces suggesting a distributed network, which likely feeds into the conduit system fed by the upglacier moulins. Drainage efficiency from the debris-covered area increased over the season but trace flow velocity remained lower than from the upper glacier moulins. Low and less-peaked melt inputs combined with the hummocky topography of the debris-covered area inhibits the formation of an efficient drainage network. These findings are relevant to regions with extensive glacial debris cover and where debris cover is expanding.

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.