Search

Understanding glaciers evolution is of major concern to evaluate their contribution to sea level rise in the context of global warming. Among the various processes involved in glacier dynamics, fractures like the calving of ice at the front of marine terminating glacier and crevasse formation affect the stress state, potentially modifying the glacier’s velocity. Crevasses also impact the melting rate. The fractures alter the roughness of the ground, increasing the amount of absorbed radiation, and open new networks in which the meltwater is likely to penetrate deeper toward the glacier bed.In this work we propose to model fractures in glacier based on finite strain elastoplasticity, using the Material Point Method (Wolper et al. 2021): we solve the classical governing equations for ice deformation in an Eulerian-Lagrangian framework and we use a strain softening Drucker-Prager constitutive model to simulate plasticity. Thanks to the Lagrangian part of the model, the fractures appear explicitly where high levels of total plastic deformation are reached. The behaviour of a glacier flowing over a step in the bedrock is investigated. The simulations show that crevasse patterns appear with regular spacing between the fractures. We perform a parametric study to determine which parameters affect the length of these patterns and potential dimensionless numbers are outlined.

Outlet glaciers and ice streams transport ice from the ice sheets to the ocean, where the glaciers lose mass by iceberg calving. Sermeq Kujalleq (SKK, Jakobshavn Isbræ) is one of the largest and most dynamic ice streams of the Greenland Ice Sheet with ice flow velocities up to 40 m/day. With extensive fieldwork and detailed repeated UAV surveys we aim at understanding the complex processes occuring at the ice stream margins and at the calving front of SKK. Such processes are often neclected in numerical models inducing uncertainties in projections of the ice sheet evolution.Within the framework of the COEBELI project we conducted drone photogrammetry surveys in July 2022 at SKK along other field measurements including in-situ GPS, GPRI, seismometers, and time-lapse imagery. Despite challenging weather conditions and constraints due to flying restrictions, we acquired 17 repeated flight surveys over the calving front of SKK during two weeks. As a result, we produced a large imagery data set, which was processed to infer high-resolution ortho-images and digital elevation models (DEM). Comparing the different products enables us to estimate changes in surface topography and ice dynamics. During the observation period several large calving events occurred allowing us to investigate the interaction between frontal processes and ice flow dynamics. With the very detailed data we can study crevasse opening, acceleration at the front, weaknesses in the ice and their origin, and the reaction of the glacier to large calving events.

Outlet glaciers and ice streams are the main channels through which ice sheets transport their mass towards the ocean. One of Greenland’s largest outlet glaciers Sermeq Kujalleq in Kangia (Jakobshavn Isbrae) has been broadly researched after experiencing a rapid retreat of the terminus and accompanying speedup to up to 40 m/day in the early 2000’s. However, such short-term ice dynamic variations remain poorly understood making numerical models difficult to constrain and predictions on future sea-level rise uncertain.The short-term ice dynamics of Sermeq Kujalleq consists in transient states and can only be captured by in-situ measurements of high spatial and temporal resolution. Glacier seismology has proven to be a valuable tool to study these dynamics, it provides data with a high temporal resolution and can provide information on processes happening below the ice surface. Within the COEBELI project we combine passive glacier seismology with global navigation satellite system (GNSS) receivers, long-range drones, time-lapse cameras and terrestrial radar interferometry to capture processes such as calving and basal sliding at their respective timescales.Here, we present results from a multi-array seismic deployment at Sermeq Kujalleq in Summer 2022. From May until September two arrays were deployed in the upstream part of the fast-flowing ice stream (>22 km from calving front) and one array on slower moving ice North of the main trunk. For a 3-week period in July, four more arrays were deployed on the fast-flowing ice stream closer to the calving front (50 Hz) icequakes. By studying these different signals, we are able to better constrain the processes and forces that control fluctuating ice-flow velocity and calving events.

Sermeq Kujalleq in Kangia (Jakobshavn Isbræ), Greenland is one of the most studied glaciers in the world mainly due to its recent retreat associated with extremely fast ice stream flow and high solid ice discharge. However, large limitations remain in the understanding of its short-term ice dynamics as the study of sub-daily variations, generally undetectable in spaceborne observations, requires high-rate field measurements that are challenging to acquire. Here, we present glacier surface velocities determined in Post-Processed Kinematic (PPK) mode from eight autonomous Global Navigation Satellite System (GNSS) stations deployed in July 2022 along the ice stream at a distance of 4 to 30 kilometers from the calving front. During this field campaign, we identified an 8-hour-long glacier speedup which was recorded at all GNSS stations and reached up to 11% of the pre-event velocity, followed by a 12-hour-long slowdown of similar magnitude. We further found the peak velocity was first measured at a GNSS station 16 kilometers away from the calving front, then recorded consecutively at each of the three other downstream GNSS stations with a positive time lag corresponding to a ~3 km/h wave propagation speed. At the station closest to the calving front, the timing of peak velocity corresponded to the occurrence of large-scale calving events. We further present line-of-sight glacier surface velocities measured along three shear margin transects with a terrestrial radar interferometer deployed simultaneously with the GNSS array. Across all profiles, we observed a widespread and simultaneous response of fast- and slow-moving ice suggesting a strong coupling between the main trunk and the shear margins of the ice stream.

Here, we present a compilation of 85 ice temperature profiles from 79 boreholes from the Greenland Ice Sheet and peripheral ice caps, as well as local ice caps in the Canadian Arctic. Only 25 profiles (32 %) were previously available in open-access data repositories. The remaining 54 profiles (68 %) are being made digitally available here for the first time. These newly available profiles, which are associated with pre-2010 boreholes, have been submitted by community members or digitized from published graphics and/or data tables. All 85 profiles are now made available in both absolute (meters) and normalized (0 to 1 ice thickness) depth scales, and are accompanied by extensive metadata. This metadata includes a transparent description of data provenance. The ice temperature profiles span 70 years, with the earliest profile being from 1950 at Camp VI, West Greenland. To highlight the value of this database in evaluating ice flow simulations, we compare the ice temperature profiles from the Greenland Ice Sheet with an ice flow simulation by the Parallel Ice Sheet Model (PISM). We find a cold bias in modeled near-surface ice temperatures within the ablation area, a warm bias in modeled basal ice temperatures at inland cold-bedded sites, and an apparent underestimation of deformational heating in high-strain settings. These biases provide process-level insight on simulated ice temperatures., The Cryosphere Discussions, ISSN:1994-0432, ISSN:1994-0440

Sermeq Kujalleq in Kangia (Jakobshavn Isbræ), Greenland has been extensively investigated over the past decades due to its recent retreat associated with extremely fast ice stream flow and high solid ice discharge. However, its short-term dynamics still remain poorly understood as they consist in transient states that can only be captured by high spatial and temporal in situ measurements. In the new COEBELI project, we aim at combining high resolution field data sets from seismic arrays, global navigation satellite system (GNSS) receivers, long-range uncrewed aerial vehicles and terrestrial radar interferometers (TRI) to achieve a comprehensive and detailed study of the short-term ice stream dynamics. Here, we present TRI and GNSS retrievals of surface velocity and elevation acquired during the first, exploratory field campaign of the COEBELI project in summer 2021. Seven kilometers away from the calving front, we specifically identified a slowdown of 1.12 m d-1 within a single day in the main trunk of the ice stream. While the absolute slowdown is larger in the main trunk than in the outer area of the shear margin (1.12 m d-1 versus 0.75 m d-1), it corresponds to a larger fraction of the pre-slowdown velocity in the latter zone (-4.48% versus -7.94%). We further discuss the challenges associated with the acquisition, processing and analysis of high-resolution data sets for the study of such complex and dynamic environments.

International audience; The complexity of karst groundwater flow modelling is reflected by the amount of simulation approaches. The goal of the Karst Modelling Challenge (KMC) is comparing different approaches on one single system using the same data set. Thirteen teams with different computational models for simulating discharge variations at karst springs have applied their respective models on one single data set coming from the Milandre Karst Hydrogeological System (MKHS). The approaches include neural networks, reservoir models, semi-distributed models and fully distributed groundwater models. Four and a half years of hourly or daily meteorological input and hourly discharge data were provided for model calibration. The validation comprised forecasting one year of discharge, without the observed discharge data. The model performance was evaluated using the volume conservation, Nash-Sutcliffe efficiency (NSE) and the Kling-Gupta efficiency (KGE) applied on the total discharge and individual flow components. As a result, the comparison of model performances is a challenging task due to the differences in the model architecture but also required time steps: some of the models require aggregated daily steps while others could be run using hourly data, which provided some interesting differences depending on how the data was transformed. The use of instantaneous data (e.g. value at noon) produces less bias that averaging hourly data over one day. The transformation of hourly into daily data produces a decrease of Nash and KGE of 0.05 to 0.08 (i.e. from 1 to ~0.93). The resulting simulations (forecasted values for year 2016) produced KGEs ranging between 0.83 and 0.37 (0.83 to −0.24 for NSE). Although the simulations matched the monitored flows reasonably well, most models struggled to simulate baseflow conditions accurately. In general, the models that performed the best for this exercise were the global ones (Gardenia and Varkarst), with a limited number of parameters, which can be calibrated using automatic calibration procedures. The neural network models also showed a fair potential, with one providing reasonable results despite the relatively short dataset available for warming-up (4.5 years). Semi-and fully distributed models also suggested that with some more effort they could perform well. The accuracy of model predictions does not seem to increase by using models with more than 9–12 calibration parameters. An evaluation of the relative errors between the forecasted and the observed values revealed that for most models, 50% of the forecasted values contained more than 50% of difference against the observed discharge rate, with 25% having a difference larger than 100%. A significant part of the poorly forecasted values corresponded to base-flow which was surprising given that as base-flow is generally much easier to predict than peak flow. Hence, this shows that modelling approaches and criteria for the calibration are too oriented towards peak-flow sections of the hydrographs, and that improvements could be gained by more focus on the base-flow.

Glaciers calving icebergs into the ocean significantly contribute to sea-level rise and can trigger tsunamis, posing severe hazards for coastal regions. Computational modeling of such multiphase processes is a great challenge involving complex solidfluid interactions. Here, a new continuum damage Material Point Method has been developed to model dynamic glacier fracture under the combined effects of gravity and buoyancy, as well as the subsequent propagation of tsunami-like waves induced by released icebergs. We reproduce the main features of tsunamis obtained in laboratory experiments as well as calving characteristics, the iceberg size, tsunami amplitude and wave speed measured at Eqip Sermia, an ocean-terminating outlet glacier of the Greenland ice sheet. Our hybrid approach constitutes important progress towards the modeling of solidfluid interactions, and has the potential to contribute to refining empirical calving laws used in large-scale earth-system models as well as to improve hazard assessments and mitigation measures in coastal regions, which is essential in the context of climate change., Communications Earth & Environment, 2, ISSN:2662-4435

Glacier calving plays a key role in the recently observed dynamic massloss of the Greenland ice sheet. Calving waves, generated by thesudden detachment of ice from the glacier terminus, can reach tens ofmeters of height and have devastating effects upon impact onsurrounding shores. In this study, we describe a new method for thedetection of source location and timing of calving waves, and theanalysis of their magnitude and spreading properties using aterrestrial radar interferometer (TRI). This method was applied to11,500 minute-interval TRI acquisitions from Eqip Sermia, Greenland.More than 2,000 calving waves were detected within sevendays. Quantitative assessment with a Wave Power Index (WPI) showedspatially distinctive patterns: the sector of the calving front endingin deep water shows a higher wave activity (+49%) with highercumulative WPI (+34%) than the shallow sector. In combination witha detection of meltwater plume locations, we highlighted a 2.3 timeshigher occurrence of visible meltwater plumes in the deep sector than theshallow one. We found both the cumulated WPI and the number of wavesto increase by more than 80% in the presence of a meltwater plumein the deep sector while only by 30% in the shallow sector. Wetherefore explain the higher calving activity in the deep sector to be strongly related to a combination of higher occurrence of meltwater plumes and more efficient calving enhancement linked to better connections to deep warm waters.

The recent rapid retreat of many Arctic outlet glaciers has been attributed to increased oceanic melt, but the relationship between oceanic melt and iceberg calving remains poorly understood. Here, we employ a transient finite‐element model that simulates oceanic melt and ice break‐off at the terminus. The response of an idealized tidewater glacier to various submarine melt rates and seasonal variations is investigated. Our modeling shows that for zero to low oceanic melt, the rate of volume loss at the front is similar or higher than for intermediate oceanic melt rates. Only very high melt rates lead to increasing volume losses. These results highlight the complex interplay between oceanic melt and calving and question the general assumption that increased submarine melt leads to higher calving fluxes and enhanced retreat. Models for tidewater glacier evolution should therefore consider calving and oceanic melt as tightly coupled processes rather than as simple, additive parametrizations.

We observed two outlet glaciers in West- and Nordwest-Greenland with a terrestrial radar interferometer (TRI), pressure sensors and time-lapse cameras over six and two years, respectively. The resulting detailed dataset provides us with insights on the calving process and the changes in front geometry over the last years. Since the two glaciers are characterised by different geometries and velocity fields, the influence of those parameters on the calving process can be investigated. The combination of the three different observation methods enable us to overcome their individual disadvantages. With the time-lapse camera taking pictures of the glacier front every 10 seconds, we detect all calving events of different sizes and styles but cannot quantify the volume. We used the TRI to quantify the volumes of aerial calving events by DEM differentiation. Further, calving waves measured with pressure sensors are used to distinguish between different calving types. We develop a relationship between calving volumes and wave heights and use this as an additional indirect method to estimate calving volumes. We find that the calving style and size as well as the front geometry is mainly controlled by the bed topography and the presence of a subglacial discharge plume. The location of the plume is observed to migrate from year to year, which leads also to changes in the calving pattern. Calving style and pattern as well as glacier velocity fields and geometry changes are additionally compared with environmental conditions such as the air temperature and the presence of ice-mélange in the proglacial fjord. In years with an early spring we find different front characteristics and calving patterns than for years with colder conditions.

Calving is a crucial process for the recently observed dynamic mass loss changes of the Greenland ice sheet. Despite its importance for global sea level change, major limitations in understanding the process of calving remain. This study presents high-resolution calving event data and statistics recorded with a terrestrial radar interferometer at the front of Eqip Sermia, a marine-terminating outlet glacier in Greenland. The derived digital elevation models with a spatial resolution of several metres recorded at 1 min intervals were processed to provide source areas and volumes of 906 individual calving events during a 6 d period. The calving front can be divided into sectors ending in shallow and deep water with different calving statistics and styles. For the shallow sector, characterized by an inclined and very high front, calving events are more frequent and larger than for the vertical ice cliff of the deep sector. We suggest that the calving volume deficiency of 90 % relative to the estimated ice flux in our observations of the deep sector is removed by oceanic melt, subaquatic calving, and small aerial calving events. Assuming a similar ice thickness for both sectors implies that subaqueous mass loss must be substantial for this sector with a contribution of up to 65 % to the frontal mass loss. The size distribution of the shallow sector is represented by a log-normal model, while for the deep sector the log-normal and power-law model fit well, but none of them are significantly better. Variations in calving activity and style between the sectors seem to be controlled by the bed topography and the front geometry. Within the short observation period no simple relationship between environmental forcings and calving frequency or event volume could be detected., The Cryosphere, 14 (3), ISSN:1994-0416, ISSN:1994-0424

Large, plume-like internal ice-layer-structures have been observed in radar images from both Antarctica and Greenland, rising from the ice-sheet base to up to half of the ice thickness. Their origins are not yet understood. Here, we simulate their genesis by basal freeze-on using numerical ice-flow modelling and analyse the transient evolution of the emerging ice-plume and the surrounding ice-layer structure as a function of both freeze-on rate and ice flux. We find good agreement between radar observations, modelled ice-plume geometry and internal layer structure, and further show that plume height relates primarily to ice-flux and only secondarily to freeze-on. An in-depth analysis, performed for Northern Greenland of observed spatial plume distribution related to ice flow, basal topography and water availability supports our findings regarding ice flux and suggests freeze-on is controlled by ascending subglacial water flow. Our results imply that widespread basal freeze-on strongly affects ice stratigraphy and consequently ice-core interpretations., Subsurface ice-sheet radar images reveal large plume-shaped bodies rising from the base, with their origin not yet understood. Here, the authors show that freeze-on of water at the ice-sheet base combined with ice-flux explains the vertical extent, shape and structure of the observed plumes.

The impact of summer surface melt on the dynamics of the Greenland Ice Sheet is modulated by the state of the subglacial hydrologic system. Studies of ice motion indicate that efficiency of the subglacial system increases over the melt season, decreasing the sensitivity of ice motion to surface melt. However, these inferences are based on limited indirect observations of the subglacial hydrologic system that leave many factors poorly constrained, particularly the presence and stability of subglacial channels. Here we use observations from 11 GPS stations, from which we derive ice velocity, longitudinal strain rates, and basal uplift, alongside observations of surface ablation and supraglacial lake drainage events, to explore the coevolution of ice motion and the subglacial hydrologic system in the Pakitsoq region of western Greenland during the 2011 melt season. We observe ice acceleration after the onset of local surface melting, followed by gradual ice deceleration, consistent with the pattern expected from increased subglacial drainage efficiency. Supraglacial lake drainages appear to precipitate ice deceleration and increased basal traction, suggesting that lake drainages effectively reorganize the local subglacial hydrologic system into a more efficient state that persists through the remainder of the melt season. At high elevations, ice velocity and inferred basal uplift suggest that continued cavity growth or sediment behavior, not subglacial channelization, drive the apparent increase in subglacial efficiency. Our results further indicate that these transient perturbations are critical in the seasonal evolution of ice motion.

Ocean terminating glaciers in Arctic regions have undergone rapid dynamic changes in recent years, which have been related to a dramatic increase in calving rates. Iceberg calving is a dynamical process strongly influenced by the geometry at the terminus of tidewater glaciers. We investigate the effect of varying water level, calving front slope and basal sliding on the stress state and flow regime for an idealized grounded ocean-terminating glacier and scale these results with ice thickness and velocity. Results show that water depth and calving front slope strongly affect the stress state while the effect from variations in basal sliding is much smaller. An increased relative water level or a reclining calving front slope strongly decrease the stresses and velocities in the vicinity of the terminus and hence have a stabilizing effect on the calving front. We find that surface stress magnitude and distribution are determined by solely the water depth relative to ice thickness for simple geometries. Based on this scaled relationship for the stress peak at the surface, and assuming a critical stress for damage initiation, we propose a simple and new parametrization for calving rates for grounded tidewater glaciers that is in good agreement with observations.

Moulins permit access of surface meltwater to the glacier bed, causing basal lubrication and ice speedup in the ablation zone of western Greenland during summer. Despite the substantial impact of moulins on ice dynamics, the conditions under which they form are poorly understood. We assimilate a time-series of ice surface velocity from a network of eleven Global Positioning System receivers into an ice sheet model to estimate ice sheet stresses during winter, spring, and summer in an approx. 30 x 10 km region. Surface-parallel von Mises stress increases slightly during spring speedup and early summer, sufficient to allow formation of 16% of moulins mapped in the study area. In contrast, 63% of moulins experience stresses over the tensile strength of ice during a short (hours) supraglacial lake drainage event. Lake drainages appear to control moulin density, which is itself a control on subglacial drainage efficiency and summer ice velocities.

Measuring the ice flow motion accurately is essential to better understand the time evolution of glaciers and ice sheets, and therefore to better anticipate the future consequence of climate change in terms of sea-level rise. Although there exist a variety of remote sensing methods to fill this task, in-situ measurements are always needed for validation or to capture high temporal resolution movements. Yet glaciers are in general hostile environments where the installation of instruments might be tedious and risky when not impossible. Here we report the first-ever in-situ measurements of ice flow motion using a remotely controlled Unmanned Aerial Vehicle (UAV). We used a multicopter UAV to land on a highly crevassed area of Eqip Sermia Glacier, West Greenland, to measure the displacement of the glacial surface with the aid of an on-board differential GNSS receiver. Despite the unfortunate loss of the UAV, we measured approximately 70 cm of displacement over 4.36 hours without setting foot onto the glacier – a result validated by applying UAV photogrammetry and template matching techniques. Our study demonstrates that UAVs are promising instruments for in-situ monitoring, and have a great potential for capturing short-term ice flow variations in inaccessible glaciers – a task that remote sensing techniques can hardly achieve.

In previous studies, the incorporation of polychlorinated biphenyls (PCBs) has been quantified in the accumulation areas of Alpine glaciers. Here, we introduce a model framework that quantifies mass fluxes of PCBs in glaciers and apply it to the Silvretta glacier (Switzerland). The models include PCB incorporation into the entire surface of the glacier, downhill transport with the flow of the glacier ice, and chemical fate in the glacial lake. The models are run for the years 1900-2100 and validated by comparing modeled and measured PCB concentrations in an ice core, a lake sediment core, and the glacial streamwater. The incorporation and release fluxes, as well as the storage of PCBs in the glacier increase until the 1980s and decrease thereafter. After a temporary increase in the 2000s, the future PCB release and the PCB concentrations in the glacial stream are estimated to be small but persistent throughout the 21st century. This study quantifies all relevant PCB fluxes in and from a temperate Alpine glacier over two centuries, and concludes that Alpine glaciers are a small secondary source of PCBs, but that the aftermath of environmental pollution by persistent and toxic chemicals can endure for decades.

Our ability to measure, quantify and assimilate hydrological properties and processes of snow in operational models is disproportionally poor compared to the significance of seasonal snowmelt as a global water resource and major risk factor in flood and avalanche forecasting. We show here that strong electrical self-potential fields are generated in melting in situ snowpacks at Rhone Glacier and Jungfraujoch Glacier, Switzerland. In agreement with theory, the diurnal evolution of self-potential magnitudes ( ∼ 60–250mV) relates to those of bulk meltwater fluxes (0–1.2 × 10−6m3s−1) principally through the permeability and the content, electrical conductivity and pH of liquid water. Previous work revealed that when fresh snow melts, ions are eluted in sequence and electrical conductivity, pH and self-potential data change diagnostically. Our snowpacks had experienced earlier stages of melt, and complementary snow pit measurements revealed that electrical conductivity ( ∼ 1–5 × 10−6Sm−1) and pH ( ∼ 6.5–6.7) as well as permeabilities (respectively ∼ 9.7 × 10−5 and ∼ 4.3 × 10−5m2 at Rhone Glacier and Jungfraujoch Glacier) were invariant. This implies, first, that preferential elution of ions was complete and, second, that our self-potential measurements reflect daily changes in liquid water contents. These were calculated to increase within the pendular regime from ∼ 1 to 5 and ∼ 3 to 5.5% respectively at Rhone Glacier and Jungfraujoch Glacier, as confirmed by ground truth measurements. We conclude that the electrical self-potential method is a promising snow and firn hydrology sensor owing to its suitability for (1) sensing lateral and vertical liquid water flows directly and minimally invasively, (2) complementing established observational programs through multidimensional spatial mapping of meltwater fluxes or liquid water content and (3) monitoring autonomously at a low cost. Future work should focus on the development of self-potential sensor arrays compatible with existing weather and snow monitoring technology and observational programs, and the integration of self-potential data into analytical frameworks., The Cryosphere, 10 (1), ISSN:1994-0416, ISSN:1994-0424

Following the general warming trend in Greenland, an increase in calving rates, retreat and ice flow has been observed at ocean-terminating outlet glaciers. These changes contribute substantially to the current mass loss of the Greenland Ice Sheet. In order to constrain models of ice dynamics as well as estimates of mass change, detailed knowledge of geometry and ice-flow are needed, in particular on the rapidly changing tongues of ocean-terminating outlet glaciers. In this study, we validate velocity estimates and spatial patterns close to the calving terminus of such an outlet derived from an iterative offset tracking method based on SAR intensity data with a collection of three independent reference measurements of glacier flow. These reference data sets are comprised of measurements from differential GPS, a Terrestrial Radar Interferometer (TRI) and repeated UAV surveys. Our approach for the SAR-velocity processing aims achieving at high spatial and temporal resolution in order to best resolve the steep velocity gradients in the terminus area and to exploit the 12 day repeat interval of the single-satellite Sentinel-1A sensor. Results from images of the medium-sized ocean terminating outlet glacier Eqip Sermia acquired by Sentinel-1A and RADARSAT-2 exhibit a mean difference of 8.7 % when compared to the corresponding GPS measurements. An areal comparison of our SAR velocity-fields with independently generated velocity maps from TRI and UAV showed a good agreement in magnitude and spatial patterns, with mean differences smaller than 0.7 md−1. In comparison with existing operational velocity products, our SAR-derived velocities showed a strongly improved spatial velocity pattern near the margins and calving front. There 10 % to 20 % higher surface ice velocities are produced, which has substantial implications on ice fluxes and on mass budget estimates of ice sheets. Further, we showed that offset tracking from SAR intensity data at a high spatio-temporal resolution is a valid method to derive glacier flow fields for fast-flowing glacier termini of outlet glaciers and, given the repeat interval of 12 days of the Sentinel-1A sensor (6 days with Sentinel-1B), has the potential to be applied operationally in a quasi-continuous mode.

The rapid drainage of supraglacial lakes introduces large pulses of meltwater to the subglacial environment and creates moulins, surface-to-bed conduits for future melt. Introduction of water to the subglacial system has been shown to affect ice flow, and modeling suggests that variability in water supply and delivery to the subsurface play an important role in the development of the subglacial hydrologic system and its ability to enhance or mitigate ice flow. We developed a fully automated method for tracking meltwater and rapid drainages in large (> 0.125 km2) perennial lakes and applied it to a 10 yr time series of ETM+ and MODIS imagery of an outlet glacier flow band in West Greenland. Results indicate interannual variability in maximum coverage and spatial evolution of total lake area. We identify 238 rapid drainage events, occurring most often at low (< 900 m) and middle (900–1200 m) elevations during periods of net filling or peak lake coverage. We observe a general progression of both lake filling and draining from lower to higher elevations but note that the timing of filling onset, peak coverage, and dissipation are also variable. Lake coverage is sensitive to air temperature, and warm years exhibit greater variability in both coverage evolution and rapid drainage. Mid-elevation drainages in 2011 coincide with large surface velocity increases at nearby GPS sites, though the relationships between ice-shed-scale dynamics and meltwater input are still unclear.

Ice temperature profiles from the Greenland Ice Sheet contain information on the deformation history, past climates and recent warming. We present full-depth temperature profiles from two drill sites on a flow line passing through Swiss Camp, West Greenland. Numerical modeling reveals that ice temperatures are considerably higher than would be expected from heat diffusion and dissipation alone. The possible causes for this extra heat are evaluated using a Lagrangian heat flow model. The model results reveal that the observations can be explained with a combination of different processes: enhanced dissipation (strain heating) in ice-age ice, temperate paleo-firn, and cryo-hydrologic warming in deep crevasses.

Penetration of surface meltwater to the bed of the Greenland Ice Sheet each summer causes an initial increase in ice speed due to elevated basal water pressure, followed by slowdown in late summer that continues into fall and winter. While this seasonal pattern is commonly explained by an evolution of the subglacial drainage system from an inefficient distributed to efficient channelized configuration, mounting evidence indicates that subglacial channels are unable to explain important aspects of hydrodynamic coupling in late summer and fall. Here we use numerical models of subglacial drainage and ice flow to show that limited, gradual leakage of water and lowering of water pressure in weakly connected regions of the bed can explain the dominant features in late and post melt season ice dynamics. These results suggest that a third weakly connected drainage component should be included in the conceptual model of subglacial hydrology., Surface meltwater draining to the bed of the Greenland Ice Sheet each summer causes ice flow changes inconsistent with the prevailing theory of channelizing subglacial drainage. Here, the authors show this is caused by limited, gradual leakage of water from previously ignored weakly connected regions of the bed.

Penetration of surface meltwater to the bed of the Greenland Ice Sheet each summer causes an initial increase in ice speed due to elevated basal water pressure, followed by slowdown in late summer that continues into fall and winter. While this seasonal pattern is commonly explained by an evolution of the subglacial drainage system from an inefficient distributed to efficient channelized configuration, mounting evidence indicates that subglacial channels are unable to explain important aspects of hydrodynamic coupling in late summer and fall. Here we use numerical models of subglacial drainage and ice flow to show that limited, gradual leakage of water and lowering of water pressure in weakly connected regions of the bed can explain the dominant features in late and post melt season ice dynamics. These results suggest that a third weakly connected drainage component should be included in the conceptual model of subglacial hydrology.

Glacier calving can cause violent tsunami waves which, upon landfall, can cause severe destruction. Here we present data acquired during a calving event from Eqip Sermia, an ocean-terminating glacier in west Greenland. During an exceptionally well-documented event, the collapse of 9 × 105 m3 ice from a 200 m high ice cliff caused a tsunami wave of 50 m height, traveling at a speed of 25–33 m s−1. This wave was filmed from a tour boat at 800 m distance from the calving face, and simultaneously measured with a terrestrial radar interferometer and a tide gauge. Tsunami wave run-up height on the steep opposite shore at a distance of 4 km was 10–15 m, destroying infrastructure and eroding old vegetation. These observations indicate that such high tsunami waves are a recent phenomenon in the history of this glacier. Analysis of the data shows that only moderately bigger tsunami waves are to be expected in the future, even under rather extreme scenarios.

We present results from a chemical fate model quantifying incorporation of polychlorinated biphenyls (PCBs) into the Silvretta glacier, a temperate Alpine glacier located in Switzerland. Temperate glaciers, in contrast to cold glaciers, are glaciers where melt processes are prevalent. Incorporation of PCBs into cold glaciers has been quantified in previous studies. However, the fate of PCBs in temperate glaciers has never been investigated. In the model, we include melt processes, inducing elution of water-soluble substances and, conversely, enrichment of particles and particle-bound chemicals. The model is validated by comparing modeled and measured PCB concentrations in an ice core collected in the Silvretta accumulation area. We quantify PCB incorporation between 1900 and 2010, and discuss the fate of six PCB congeners. PCB concentrations in the ice core peak in the period of high PCB emissions, as well as in years with strong melt. While for lower-chlorinated PCB congeners revolatilization is important, for higher-chlorinated congeners, the main processes are storage in glacier ice and removal by particle runoff. This study gives insight into PCB fate and dynamics and reveals the effect of snow accumulation and melt processes on the fate of semivolatile organic chemicals in a temperate Alpine glacier.

Our ability to measure, quantify and assimilate hydrological properties and processes of snow in operational models is disproportionally poor compared to the significance of seasonal snowmelt as a global water resource and major risk factor in flood and avalanche forecasting. Encouraged by recent theoretical, modelling and laboratory work, we show here that the diurnal evolution of aerially-distributed self-potential magnitudes closely track those of bulk meltwater fluxes in melting in-situ snowpacks at Rhone and Jungfraujoch glaciers, Switzerland. Numerical modelling infers temporally-evolving liquid water contents in the snowpacks on successive days in close agreement with snow-pit measurements. Muting previous concerns, the governing physical and chemical properties of snow and meltwater became temporally invariant for modelling purposes. Because measurement procedure is straightforward and readily automated for continuous monitoring over significant spatial scales, we conclude that the self-potential geophysical method is a highly-promising non-intrusive snow-hydrological sensor for measurement practice, modelling and operational snow forecasting.

Glacier calving can cause violent impulse waves which, upon landfall, can lead to destructive tsunami-like waves. Here we present data acquired during a calving event from Eqip Sermia, an ocean-terminating glacier in West Greenland. During an exceptionally well documented event, the collapse of 9 × 105 m3 ice from a 200 m high ice cliff caused an impulse wave of 50 m height, traveling at a speed of 25–30 m s-1. This wave was filmed from a tour boat in 800 m distance from the calving face, and simultaneously measured with a terrestrial radar interferometer and a tide gauge. Tsunami wave run-up height on the steep opposite shore in 4 km distance was 10–15 m, destroying infrastructure and eroding old vegetation. These observations indicate that such high tsunami waves are a recent phenomenon in the history of this glacier. Analysis of the data shows that only moderately bigger tsunami waves are to be expected in the future, even under rather extreme scenarios.