Annual Surface Velocities Research Articles

Terminus thinning drives recent acceleration of a Greenlandic lake-terminating outlet glacier

Abstract Ice-contact proglacial lakes affect ice dynamics and the transition of glacier termini from land- to lake-terminating has been shown to cause ice flow acceleration. In recent decades, the number and size of Greenlandic ice-marginal lakes has increased, highlighting the need to further understand these lake-terminating ice-margins as their influence on ice sheet mass balance increases. Here, time series of satellite-derived observations of ice velocity, surface elevation, and terminus position were generated at a lake-terminating outlet glacier, Isortuarsuup Sermia, and the nearby land-terminating Kangaasarsuup Sermia in south-west Greenland. At Isortuarsuup Sermia, annual surface velocity at the terminus increased by a factor of 2.5 to 214 ± 4 m yr−1 (2013–2021), with the magnitude of this acceleration declining with distance up-glacier. Meanwhile, near-terminus surface elevation changed at a rate of −2.3 ± 1.1 m yr−1 (2012–2021). Conversely, velocity change at Kangaasarsuup Sermia was minimal, while surface elevation change was approximately half at comparable elevations (−1.2 ± 0.3 m yr−1). We attribute these dynamic differences to thinning at Isortuarsuup Sermia and subsequent retreat from a stabilising sublacustrine moraine, and emphasise the potential of proglacial lakes to enhance future rates of mass loss from the Greenland Ice Sheet.

The Expanding of Proglacial Lake Amplified the Frontal Ablation of Jiongpu Co Glacier since 1985

In High Mountain Asia, most glaciers and glacial lakes have undergone rapid variations throughout changes in the climate. Unlike land-terminating glaciers, lake-terminating glaciers show rapid shrinkage due to dynamic interactions between proglacial lakes and glacier dynamics. In this study, we conducted a detailed analysis of the changes in the surface elevation, velocity, and especially frontal ablation on Jiongpu Co lake-terminating glacier. The results show that the Jiongpu Co glacier has twice as much negative mass balance compared to other glaciers, and the annual surface velocity has anomalously increased (3.6 m a−1 per decade) while other glaciers show a decreased trend. The frontal ablation fraction in the net mass loss of the Jiongpu Co glacier increased from 26% to 52% with the accelerated expansion of the proglacial lake. All available evidence indicates the presence of positive feedback between the proglacial lake and its host glacier. Our findings highlight the existence of proglacial lake affects the spatial change patterns of the lake-terminating glacier. Furthermore, the ongoing enlargement of the lake area amplifies the changes associated with the evolution of the lake-terminating glacier.

Dynamic Monitoring of Laohugou Glacier No. 12 with a Drone, West Qilian Mountains, West China

Laohugou glacier No. 12 (LHG12), located in the northeast of the Qinghai–Tibet Plateau, is the largest valley glacier in the Qilian mountains. Since 1957, LHG12 has shrunk significantly. Due to the limitations of in situ observations, simulations and investigations of LHG12 have higher levels of uncertainty. In this study, consumer-level, low-altitude microdrones were used to conduct repeated photogrammetry at the lower part of LHG12, and a digital orthophoto map (DOM) and a digital surface model (DSM) with a resolution at the centimeter scale were generated, from 2017 to 2021. The dynamic parameters of the glacier were detected by artificial and automatic extraction methods. Using a combination of GNSS and drone-based data, the dynamic process of LHG12 was analyzed. The results show that the terminus of LHG12 has retreated by 194.35 m in total and by 19.44 m a−1 on average during 2008–2021. The differential ablation leading to terminus retreat distance markedly increased during the study period. In 2019–2021, the maximum annual surface velocity was 6.50 cm day−1, and during ablation season, the maximum surface velocity was 13.59 cm day−1, 52.17% higher than it is annually. The surface parameters, motion, and mass balance characteristics of the glacier had significant differences between the west and east branches. The movement in the west branch is faster than it is in the east branch. Because of the extrusion of the two ice flows, there is a region with a faster surface velocity at the ablation area. The ice thickness of LHG12 is decreasing due to intensified ablation, leading to a deceleration in the surface velocity. In large glaciers, this phenomenon is more obvious than it is in small glaciers in the Qilian mountains.

SURFACE VELOCITY DYNAMICS OF SAMUDRA TAPU GLACIER, INDIA FROM 2013 TO 2017 USING LANDSAT-8 DATA

Abstract. In glacier dynamics, surface velocity of glacier is an important parameter to understand the behaviour of glacier in absence of mass balance and long-term glacier area change information. In present study, surface velocity of Samudra Tapu Glacier, India is estimated using freely available Landsat-8 OLI (PAN) images during 2013–2017. To estimate surface velocity, open source COSI-Corr tool is used which is based on cross-correlation algorithm. Maximum annual surface velocity estimated is 55.68 ± 4.01 m/year during 2015–2016 while the minimum surface velocity being 44.99 ± 4.67 m/year in 2016–2017. The average annual velocity during 2013–2017 was 50.51 ± 4.49 m/year which is higher than other glaciers in Chandra basin. The variation in annual surface velocity is analysed which not only depends on mass loss but also on temperature, pressure and internal drainage. Further, as one moves opposite to glacier terminus, the surface velocity increases with the increase in glacier elevation and slope. The higher surface velocity can be attributed to the fact that Samudra Tapu is a top-heavy glacier based on HI index analysis having larger accumulation area along with high glacier ice-thickness.

Spatiotemporal variation in surface velocity in Chandra basin glacier between 1999 and 2017 using Landsat-7 and Landsat-8 imagery

Glacier surface velocity is one of the main parameters of the glaciers used to interpret the glacier dynamics. In the present research paper, surface velocity of four glaciers, namely Bara Shigri, Chhota Shigri, Gepang Gath and Hamtah in Chandra basin between 1999 and 2017 during two time spans of 1999–2002 and 2013–2017, is estimated to analyze spatiotemporal characteristics of surface velocity of these glaciers using a cross-correlation algorithm applied on panchromatic images of Landsat-7 Enhanced Thematic Mapper Plus (ETM+) and Landsat-8 Operational Land Imager (OLI) sensors. A decrease of ∼34% and ∼33% in average annual surface velocity is observed during 2013–2017 in comparison with the average annual surface velocity during 1999–2002 for Bara Shigri and Hamtah glaciers mainly due to lowering of glacier surface. However, Gepang Gath has shown a very low reduction of ∼3% during the same period due to the presence of proglacial lake in its terminus. The analysis of transverse velocity profile of all the four glaciers is also carried out to determine the dominant factor responsible for glacier movement.

Modelling seasonal meltwater forcing of the velocity of land-terminating margins of the Greenland Ice Sheet

Abstract. Surface runoff at the margin of the Greenland Ice Sheet (GrIS) drains to the ice-sheet bed, leading to enhanced summer ice flow. Ice velocities show a pattern of early summer acceleration followed by mid-summer deceleration due to evolution of the subglacial hydrology system in response to meltwater forcing. Modelling the integrated hydrological–ice dynamics system to reproduce measured velocities at the ice margin remains a key challenge for validating the present understanding of the system and constraining the impact of increasing surface runoff rates on dynamic ice mass loss from the GrIS. Here we show that a multi-component model incorporating supraglacial, subglacial, and ice dynamic components applied to a land-terminating catchment in western Greenland produces modelled velocities which are in reasonable agreement with those observed in GPS records for three melt seasons of varying melt intensities. This provides numerical support for the hypothesis that the subglacial system develops analogously to alpine glaciers and supports recent model formulations capturing the transition between distributed and channelized states. The model shows the growth of efficient conduit-based drainage up-glacier from the ice sheet margin, which develops more extensively, and further inland, as melt intensity increases. This suggests current trends of decadal-timescale slowdown of ice velocities in the ablation zone may continue in the near future. The model results also show a strong scaling between average summer velocities and melt season intensity, particularly in the upper ablation area. Assuming winter velocities are not impacted by channelization, our model suggests an upper bound of a 25 % increase in annual surface velocities as surface melt increases to 4× present levels.

Changing contribution of peak velocity events to annual velocities following a multi-decadal slowdown at White Glacier

ABSTRACTAs the focus of intensive glaciological studies in the 1960–70s, White Glacier on Axel Heiberg Island, Canada, has played an important role in understanding the dynamics of mostly-cold polythermal glaciers in the high Arctic. In this study, we examine the magnitude, duration and timing of peak velocity events in the summers of 2013–15 using continuous dual-frequency GPS observations, and compare them with similar measurements made in 1968. Summer speed-up events in 1968 and 2014, in which ice velocities reached 200% above winter values, were found to occur in conjunction with formation and drainage of an ice-marginal lake. Despite thinning of the glacier by >20 m and a decrease in annual surface velocities of 15–35% since the 1960s, the relative magnitude and duration of these peak events has increased, particularly at lower elevations, in comparison with the observations at the same locations many decades ago. Given the long-term slowdown of the glacier, the relative contribution of summer displacement to the net annual motion has therefore increased significantly, with summer motion over the span of <2 months now accounting for nearly half of the total annual displacement.

Sliding velocity fluctuations and subglacial hydrology over the last two decades on Argentière glacier, Mont Blanc area

ABSTRACTThe subglacial observatory beneath the Argentière glacier provides a rare opportunity to study the interactions between glacier sliding velocity and subglacial runoff. The sliding velocity has been monitored in this cavity almost continuously since 1997 and the resulting data indicate a decrease in annual sliding velocities over the last two decades. We found close relationships between annual surface velocity, sliding velocity and ice thickness. These relationships indicate that the ice-flow velocity changes do not depend on subglacial water runoff changes at the annual timescale. The seasonal magnitudes of sliding also show a decrease over the last two decades. At the seasonal timescale, sliding velocity increases before or simultaneously with the large runoff increase in May, indicating a distributed drainage system. Conversely, at the end of the melt season, sliding velocity continues to decrease after the runoff returns to low winter values. The simultaneous increases of runoff and sliding velocity occur mainly before the spring transition. Later, sliding velocity generally appears not to be related to water inputs coming from the surface, except for some large accelerations after midAugust that are always associated with periods of rapidly increasing water inputs to the subglacial drainage system.

Dynamics of an alpine cirque glacier

Alpine cirques are excavated by glacial erosion, a process that depends in turn on the movement of ice by basal sliding. Cirque glacier flow is usually depicted as rotational sliding of a rigid block, but this model is based on little evidence and implies unorthodox glacier behavior given typical cirque dimensions. The small (∼1 km2), temperate West Washmawapta Glacier occupies an archetypal overdeepened and “armchair-shaped” cirque in the Canadian Rockies. We measured (1) the annual surface velocity field, (2) ice thickness, (3) sliding and internal deformation at one borehole, and (4) sliding in a marginal cavity. The glacier moves slowly, with surface velocities of 3 to 10 m/yr. The maximum ice thickness (∼185 m) occurs in the center of the cirque basin and roughly coincides with the position of greatest ice flux. Using our field measurements, a standard constitutive relation for ice, and simplifying assumptions related to the depth distribution of strain rates, we approximated the driving and resisting forces acting on sections of the glacier, and inferred the general pattern of basal sliding. Sliding is minimum in the center of the cirque and increases toward the margins, especially up the stoss side of the riegel. Internal deformation accounts for all motion in the cirque center, even if an unusually low viscosity for temperate ice is assumed. Basal shear stresses tend toward 105 Pa everywhere, a typical value for mountain glaciers. Transverse and longitudinal straining are significant in some parts of the glacier. Although a component of rotational flow must occur internally, the glacier does not conform to the rotational sliding model in any essential respect.

Iceberg calving rates from northern Ellesmere Island ice caps, Canadian Arctic, 1999–2003

AbstractOptical satellite imagery was used to estimate glacier surface velocities and iceberg calving rates from Agassiz and western Grant Ice Caps, Nunavut, Canada, between 1999 and 2003. The largest mean annual surface velocities ranged from ∼400 to 700 m a−1, but velocities in the ∼100–200 m a−1 range were common. Summer velocities were up to an order of magnitude larger than annually averaged velocities. The highest velocity (∼1530 m a−1) was measured on the floating tongue of Lake Tuborg Glacier between 19 July and 19 August 2001. Calving rates from individual glaciers varied by up to a factor of two between successive years. Summer calving rates were ∼2–8 times larger than annual average rates. The average ratio of the calving flux due to terminus-volume change to that due to ice flow through the glacier terminus was ∼0.81 for the annual rates and ∼1.71 for summer rates. The estimated mean annual calving rate from Agassiz Ice Cap in the period 1999–2002 was 0.67 ± 0.15 km3 a−1, of which ∼54% emanated from Eugenie Glacier alone. This total rate is similar to a previously estimated calving rate from Devon Ice Cap.

Monitoring of the Taconnaz ice fall (French Alps) using measurements of mass balance, surface velocities and ice cliff position

Glacier de Taconnaz is a hanging glacier in the Mont-Blanc area with an upper accumulation area stretching from Dôme du Goûter (4300 m a.s.l.) down to a wide ice cliff at an altitude of about 3300 m a.s.l., over which flows most of the ice accumulated upstream. During winter, ice blocks breaking off the cliff often trigger off large avalanches constituted of a mixture of snow and ice which can be devastating for the downward inhabited areas in the Chamonix valley as was the case in 1988 and 1999. The approach proposed here is twofold. It first consists of quantifying the flow regime over the accumulation area and consequently the yearly amount of ice falling off the cliff. For this purpose, a network of large stakes was first installed in June 2001 and regularly updated since. From the measured height above surface corrected with time series mass balance of two nearby glaciers, a yearly net mass balance map over the accumulation area has been established and converted into a mean discharge of 1.2 ± 0.3 × 106 m3 year− 1 at the cliff. Moreover, monthly remote topographic surveys of the same stakes provided annual surface velocities, which after some assumptions led to a similar independent amount of discharge. The second part consists of trying to assess the frequency at which the collapses occur from similar monthly topographic measurements of the cliff edge position. Although the time series are still short, a single period of about 180 days seems to control the major ice falls. These data also allow for estimation of ice volumes breaking off which, associated with flow values, could help constraining the frequency of falling events.

Intra-annual and intra-seasonal flow dynamics of a High Arctic polythermal valley glacier

AbstractMeasurements of surface dynamics on polythermal John Evans Glacier, Nunavut, Canada, over two winter periods and every 7–10 days throughout two melt seasons (June–July 2000, 2001) provide new insight into spatio-temporal patterns of High Arctic glacier dynamics. In the lower ablation zone, mean annual surface velocities are 10–21 m a–1, but peak velocities up to 50% higher are attained during late June/early July. In the upper ablation zone and lower accumulation zone, mean annual surface velocities are typically 10–18 m a–1, and peak velocities up to 40% higher occur during late July. In the upper accumulation zone, mean annual surface velocities are 2–9 m a–1, and motion in mid- to late July exceeds this by up to 10%. Rapid drainage of ponded supraglacial water in the upper ablation zone to an initially distributed subglacial drainage system in mid-June may force excess surface motion in the warm-based lower glacier. The data indicate that the duration of the velocity response may be related to the rate of channelization of the basal drainage, and the velocity response may be transmitted up-glacier by longitudinal coupling. An increase in surface velocities in the middle glacier in late July occurs in conjunction with the opening of two further moulins in the accumulation zone.

Velocities, strain rates, stresses, crevassing and faulting on Glacier de Saint-Sorlin, French Alps, 1957–76

AbstractStake surveys on Glacier de Saint-Sorlin, French Alps, during the period 1957–76 show that annual surface velocities fit a linear vectorial model, with a term depending on the site, another on the year, and an important random component. Strain rates, viscosities and stresses at shallow depth are computed using strain triangles of hectometric (102 m) size. Between 1961/62 and 1972/73 the isotropic point, where streamlines cease to converge, moved downstream about 200 m. This trend may explain increasing velocities. Crevasses appear where annual strain is > 1.2%. Faults limit the effective shear stress at the surface to about 0.38 bar. Ten modes of flow are distinguished, instead of only two for the two-dimensional model (compressive and extensive). The gradients of shallow stresses, which ensure extra driving forces, are computed with another mosaic of triangles of similar size (stress triangles). There are also important extra driving forces at the bottom, which force the flow to deviate from the direction of the steepest surface slope. Two criteria allow elimination of stress triangles where these unknown basal extra driving forces are important. Even so, no sliding law in terms of mean annual values can be obtained. This study shows that the classical perturbation theory, which explains advances and retreats by the arrival of kinematic waves, is unsuitable for glaciers of kilometric size.