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 englacial stratigraphic architecture of internal‐reflection horizons (IRHs) as imaged by ice‐penetrating radar (IPR) across ice sheets reflects the cumulative effects of surface mass balance, basal melt and ice flow. IRHs, considered isochrones, have typically been traced in interior, slow‐flowing regions. Here, we identify three distinctive IRHs spanning the Institute and Möller catchments that cover 50% of West Antarctica's Weddell Sea Sector and are characterised by a complex system of ice‐stream tributaries. We place age constraints on IRHs through their intersections with previous geophysical surveys tied to Byrd Ice Core, and by age‐depth modeling. We further show where the oldest ice likely exists within the region; and that Holocene ice‐dynamic changes were limited to the catchment's lower reaches. The traced IRHs from this study have clear potential to nucleate a wider continental‐scale IRH database for validating ice‐sheet models.

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.

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] We present ice velocities from a land-terminating transect extending >115 km into the western Greenland Ice Sheet during three contrasting melt years (2009–2011) to determine whether enhanced melting accelerates dynamic mass loss. We find no significant correlation between surface melt and annual ice flow. There is however a positive correlation between melt and summer ice displacement, but a negative correlation with winter displacement. This response is consistent with hydro-dynamic coupling; enhanced summer ice flow results from longer periods of increasing surface melting and greater duration ice surface to bed connections, while reduced winter motion is explicable by drainage of high basal water pressure regions by larger more extensive subglacial channels. Despite mean interannual surface melt variability of up to 70%, mean annual ice velocities changed by

Predictions of the Greenland Ice Sheet’s response to climate change are limited in part by uncertainty in the coupling between meltwater lubrication of the ice-sheet bed and ice flow1, 2, 3. This uncertainty arises largely from a lack of direct measurements of water flow characteristics at the bed of the ice sheet. Previous work has been restricted to indirect observations based on seasonal and spatial variations in surface ice velocities4, 5, 6, 7 and on meltwater flux8. Here, we employ rhodamine and sulphur hexafluoride tracers, injected into the drainage system over three melt seasons, to observe subglacial drainage properties and evolution beneath the Greenland Ice Sheet, up to 57 km from the margin. Tracer results indicate evolution from a slow, inefficient drainage system to a fast, efficient channelized drainage system over the course of the melt season. Further inland, evolution to efficient drainage occurs later and more slowly. An efficient routing of water was established up to 41 km or more from the margin, where the ice is approximately 1 km thick. Overall, our findings support previous interpretations of drainage system characteristics, thereby validating the use of surface observations as a means of investigating basal processes.

[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 Pleistocene ice sheets left a clear signature of erosion, but the rate at which ice sheets erode is difficult to determine from either paleolandscapes or observations of contemporary processes. Here we use two years of sediment flux data, derived from meltwaters emerging from an outlet glacier in west Greenland, to calculate an average rate of subglacial erosion across a catchment extending >50 km inland from the ice margin. Erosion in this zone occurs at 4.8 ± 2.6 mm a −1 , a rate 1–2 orders of magnitude greater than previous estimates of erosion rate beneath the Greenland Ice Sheet. Our results suggest that where surface meltwaters are able to access the bed, the rate of erosion by ice sheets is in keeping with the rapid erosion observed at temperate alpine glaciers. During deglacial phases, when meltwater was abundant, ice sheet margins should therefore have acted as highly efficient agents of erosion.

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

We present global positioning system observations that capture the full inland extent of ice motion variations in 2009 along a transect in the west Greenland Ice sheet margin. In situ measurements of air temperature and surface ablation, and satellite monitoring of ice surface albedo and supraglacial lake drainage are used to investigate hydrological controls on ice velocity changes. We find a strong positive correlation between rates of annual ablation and changes in annual ice motion along the transect, with sites nearest the ice sheet margin experiencing greater annual variations in ice motion (15–18%) than those above 1000 m elevation (3–8%). Patterns in the timing and rate of meltwater delivery to the ice–bed interface provide key controls on the magnitude of hydrologically-forced velocity variations at each site. In the lower ablation zone, the overall contribution of variations in ice motion to annual flow rates is limited by evolution in the structure of the subglacial drainage system. At sites in the upper ablation zone, a shorter period of summer melting and delayed establishment of a hydraulic connection between the ice sheet surface and its bed limit the timeframe for velocity variations to occur. Our data suggest that land-terminating sections of the Greenland Ice Sheet will experience increased dynamic mass loss in a warmer climate, as the behaviour that we observe in the lower ablation zone propagates further inland. Findings from this study provide a conceptual framework to understand the impact of hydrologically-forced velocity variations on the future mass balance of land-terminating sections of the Greenland Ice Sheet.

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

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 Canadian Arctic contains the largest area of land ice (~150,000 km2) on Earth outside the ice sheets of Greenland and Antarctica and is a potentially significant contributor to global sea level change. The current ice cover includes large ice caps that are remnants of the Wisconsinan Laurentide and Innuitian ice sheets, and many smaller ice caps and valley glaciers that formed during the late Holocene. Most of these ice masses have decreased in area over the past century as a result of climate warming in the first half of the 20th century and since the mid-1980s. In general, smaller ice masses have lost a higher proportion of their area, but the largest total area losses have come from the larger ice caps. Both iceberg calving and negative surface mass balances have contributed to this episode of glacier shrinkage. Long-term calving rates are not well known, however, and many tidewater glaciers exhibit velocity variability on a range of timescales that may affect calving rates. Floating ice shelves in northern Ellesmere Island have lost over 90 % of their area in the 20th century, with the most recent phase of disintegration occurring since 2000. Some fjords in the region are now ice free for the first time in over 3000 years. Regional rates of mass loss have accelerated strongly since 2005, and Canadian Arctic glaciers and ice caps have emerged as the most significant non–ice sheet contributor to the nonsteric component of global sea level rise.

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.

Changes to the dynamics of the Greenland ice sheet can be forced by various mechanisms including surface-melt-induced ice acceleration and oceanic forcing of marine-terminating glaciers. We use observations of ice motion to examine the surface melt-induced dynamic response of a land-terminating outlet glacier in southwest Greenland to the exceptional melting observed in 2012. During summer, meltwater generated on the Greenland ice sheet surface accesses the ice sheet bed, lubricating basal motion and resulting in periods of faster ice flow. However, the net impact of varying meltwater volumes upon seasonal and annual ice flow, and thus sea level rise, remains unclear. We show that two extreme melt events (98.6% of the Greenland ice sheet surface experienced melting on July 12, the most significant melt event since 1889, and 79.2% on July 29) and summer ice sheet runoff ~3.9 σ above the 1958-2011 mean resulted in enhanced summer ice motion relative to the average melt year of 2009. However, despite record summer melting, subsequent reduced winter ice motion resulted in 6% less net annual ice motion in 2012 than in 2009. Our findings suggest that surface melt-induced acceleration of land-terminating regions of the ice sheet will remain insignificant even under extreme melting scenarios.

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.

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

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

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.

Radar altimeters are one of the main tools for measuring elevation changes across the Antarctic and Greenland ice sheets and larger ice caps. A ground-based radar was deployed in autumn 2004 and spring 2006 in the percolation zone of the Greenland Ice Sheet. This radar is a high bandwidth system operating in the Ku band, the same frequency as several satellite altimeters. Measurements were made over an elevation range of 1795 to 2350 m, along with snow pit and shallow core studies. These measurements demonstrate the spatial and temporal variations in the backscatter. Relative strengths of surface and volume reflections change dramatically between spring and autumn and there is also high spatial variability across the percolation zone. The extent of percolation will affect elevation estimates made by radar altimeters.

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

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.