Sliding dominates slow-flowing margin regions, Greenland Ice Sheet

Complex motion of Greenland Ice Sheet outlet glaciers with basal temperate ice

Uncertainty associated with ice sheet motion plagues sea level rise predictions. Much of this uncertainty arises from imperfect representations of physical processes including basal slip and internal ice deformation, with ice sheet models largely incapable of reproducing borehole-based observations. Here, we model isolated three-dimensional domains from fast-moving (Sermeq Kujalleq/Store Glacier) and slow-moving (Isunnguata Sermia) ice sheet settings in Greenland. By incorporating realistic geostatistically simulated topography, we show that a spatially highly variable layer of temperate ice (much softer ice at the pressure-melting point) forms naturally in both settings, alongside ice motion patterns which diverge substantially from those obtained using smoothly varying BedMachine topography. Temperate ice is vertically extensive (>100 meters) in deep troughs but thins notably (<5 meters) over bedrock highs, with basal slip rates reaching >90 or <5% of surface velocity dependent on topography and temperate layer thickness. Developing parameterizations of the net effect of this complex motion can improve the realism of predictive ice sheet models.

Threshold response to melt drives large-scale bed weakening in Greenland

Ice speeds in Greenland are largely set by basal motion¹, which is modulated by meltwater delivery to the ice base^2-4. Evidence suggests that increasing melt rates enhance the subglacial drainage network's capacity to evacuate basal water, increasing bed friction and causing the ice to slow^5-10. This limits the potential of melt forcing to increase mass loss as temperatures increase¹¹. Here we show that melt forcing has a pronounced influence on dynamics, but factors besides melt rates primarily control its impact. Using a method to examine friction variability across the entirety of western Greenland, we show that the main impact of melt forcing is an abrupt north-to-south change in bed strength that cannot be explained by changes in melt production. The southern ablation zone is weakened by 20-40 per cent compared with regions with no melt, whereas in northern Greenland the ablation zone is strengthened. We show that the weakening is consistent with persistent basal water storage and that the threshold is linked to differences in sliding and hydropotential gradients, which exert primary control on the pressures within drainage pathways that dewater the bed. These characteristics are mainly set by whether a margin is land or marine terminating, suggesting that dynamic changes that increase mass loss are likely to occur in northern Greenland as temperatures increase. Our results point to physical representations of these findings that will improve simulated ice-sheet evolution at centennial scales.

A first constraint on basal melt-water production of the Greenland ice sheet

The Greenland ice sheet has been one of the largest sources of sea-level rise since the early 2000s. However, basal melt has not been included explicitly in assessments of ice-sheet mass loss so far. Here, we present the first estimate of the total and regional basal melt produced by the ice sheet and the recent change in basal melt through time. We find that the ice sheet's present basal melt production is 21.4 +4.4/-4.0 Gt per year, and that melt generated by basal friction is responsible for about half of this volume. We estimate that basal melting has increased by 2.9 ± 5.2 Gt during the first decade of the 2000s. As the Arctic warms, we anticipate that basal melt will continue to increase due to faster ice flow and more surface melting thus compounding current mass loss trends, enhancing solid ice discharge, and modifying fjord circulation.

Changing friction at the base of an Alpine glacier

Repeating earthquakes are a global phenomenon of tectonic faults. Multiple ruptures on the same fault asperities lead to nearly identical waveforms characteristic for these seismic events. We identify their microseismic counterparts beneath an Alpine glacier, where basal sliding accounts for a significant amount of ice flow. In contrast to tectonic faults, Alpine glacier beds are subject to large variations in sliding velocity and effective normal stresses. This leads to inter- and sub-seasonal variations in released seismic moment from stick–slip asperities, which we explain with the rate-and-state friction formalism. During summer, numerically modelled effective normal stresses at asperities are three times higher than in winter, which increases the local shear resistance by the same factor. Stronger summer asperities therefore tend to form in bed regions well connected to the efficient subglacial drainage system. Moreover, asperities organise themselves into a state of subcriticality, transferring stresses between each other. We argue that this seismic stick–slip behavior has potentially far-reaching consequences for glacier sliding and in particular for catastrophic failure of unstable ice masses.

Thermodynamics of a fast-moving Greenlandic outlet glacier revealed by fiber-optic distributed temperature sensing

Measurements of ice temperature provide crucial constraints on ice viscosity and the thermodynamic processes occurring within a glacier. However, such measurements are presently limited by a small number of relatively coarse-spatial-resolution borehole records, especially for ice sheets. Here, we advance our understanding of glacier thermodynamics with an exceptionally high-vertical-resolution (~0.65 m), distributed-fiber-optic temperature-sensing profile from a 1043-m borehole drilled to the base of Sermeq Kujalleq (Store Glacier), Greenland. We report substantial but isolated strain heating within interglacial-phase ice at 208 to 242 m depth together with strongly heterogeneous ice deformation in glacial-phase ice below 889 m. We also observe a high-strain interface between glacial- and interglacial-phase ice and a 73-m-thick temperate basal layer, interpreted as locally formed and important for the glacier's fast motion. These findings demonstrate notable spatial heterogeneity, both vertically and at the catchment scale, in the conditions facilitating the fast motion of marine-terminating glaciers in Greenland.

A slip law for hard-bedded glaciers derived from observed bed topography

Ice-sheet responses to climate warming and associated sea-level rise depend sensitively on the form of the slip law that relates drag at the beds of glaciers to their slip velocity and basal water pressure. Process-based models of glacier slip over idealized, hard (rigid) beds with water-filled cavities yield slip laws in which drag decreases with increasing slip velocity or water pressure (rate-weakening drag). We present results of a process-based, three-dimensional model of glacier slip applied to measured bed topographies. We find that consideration of actual glacier beds eliminates or makes insignificant rate-weakening drag, thereby uniting process-based models of slip with some ice-sheet model parameterizations. Computed slip laws have the same form as those indicated by experiments with ice dragged over deformable till, the other common bed condition. Thus, these results may point to a universal slip law that would simplify and improve estimations of glacier discharges to the oceans.

Data Reduction Using Statistical and Regression Approaches for Ice Velocity Derived by Landsat-8, Sentinel-1 and Sentinel-2

During the last decade, the number of available satellite observations has increased significantly, allowing for far more frequent measurements of the glacier speed. Appropriate methods of post-processing need to be developed to efficiently deal with the large volumes of data generated and relatively large intrinsic errors associated with the measurements. Here, we process and combine together measurements of ice velocity of Russell Gletscher in Greenland from three satellites—Sentinel-1, Sentinel-2, and Landsat-8, creating a multi-year velocity database with high temporal and spatial resolution. We then investigate post-processing methodologies with the aim of generating corrected, ordered, and simplified time series. We tested rolling mean and median, cubic spline regression, and linear non-parametric local regression (LOWESS) smoothing algorithms to reduce data noise, evaluated the results against ground-based GPS in one location, and compared the results between two locations with different characteristics. We found that LOWESS provides the best solution for noisy measurements that are unevenly distributed in time. Using this methodology with these sensors, we can robustly derive time series with temporal resolution of 2–3 weeks and improve the accuracy on the ice velocity to about 10 m/yr, or a factor of three compared to the initial measurements. The presented methodology could be applied to the entire Greenland ice sheet with an aim of reconstructing comprehensive sub-seasonal ice flow dynamics and mass balance.