Search

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

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.

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.

Current understanding of ice dynamics predicts that increasing availability and variability of meltwater will have an impact on basal motion and therefore on the evolution and future behavior of the Greenland ice sheet. We present measurements of ice deformation, subglacial water pressure, and surface velocity that show periodic and episodic variations on several time scales (seasonal, multiday, and diurnal). These variations, observed with GPS and sensors at different depths throughout the ice column, are not synchronous but show delayed responses of ice deformation with increasing depth and basal water pressure in antiphase with surface velocity. With the help of a Full-Stokes ice flow model, these observations are explained as ice motion in a caterpillar-like fashion. Caused by patches of different basal slipperiness, horizontal stress transfer through the stiff central part of the ice body leads to spatially varying surface velocities and ice deformation patterns. Variation of this basal slipperiness induces characteristic patterns of ice deformation variability that explain the observed behavior. Ice flow in the ablation zone of the Greenland ice sheet is therefore controlled by activation of basal patches by varying slipperiness in the course of a melt season, leading to caterpillar-like ice motion superposed on the classical shear deformation.

Simultaneous observations of moulins and boreholes in western Greenland show that water delivery to the base of the ice sheet via moulins affects short-term ice velocity fluctuations, but not late-season ice velocity decelerations. Increased meltwater delivery to the base of the Greenland Ice Sheet will increase the ice sheet velocity, accelerating its inevitable rush towards the ocean and subsequent sea level rise. Or will it? Debate on this topic is at the forefront of cryospheric research, but has been hampered by the lack of simultaneous observations of hydraulic head in moulins, vertical shafts that deliver water to the base of the ice sheet, and boreholes, which monitor basal water pressure. Lauren Andrews and colleagues now provide these observations from a small region in western Greenland and reveal that water delivery by moulins into channelized basal flow is indeed linked to short-term fluctuations in ice velocity. Late season decelerations in ice velocity, however, seem to be controlled by changes in unchannelized flow, rather than any shifts in the moulin system. Seasonal acceleration of the Greenland Ice Sheet is influenced by the dynamic response of the subglacial hydrologic system to variability in meltwater delivery to the bed1,2 via crevasses and moulins (vertical conduits connecting supraglacial water to the bed of the ice sheet). As the melt season progresses, the subglacial hydrologic system drains supraglacial meltwater more efficiently1,2,3,4, decreasing basal water pressure4 and moderating the ice velocity response to surface melting1,2. However, limited direct observations of subglacial water pressure4,5,6,7 mean that the spatiotemporal evolution of the subglacial hydrologic system remains poorly understood. Here we show that ice velocity is well correlated with moulin hydraulic head but is out of phase with that of nearby (0.3–2 kilometres away) boreholes, indicating that moulins connect to an efficient, channelized component of the subglacial hydrologic system, which exerts the primary control on diurnal and multi-day changes in ice velocity. Our simultaneous measurements of moulin and borehole hydraulic head and ice velocity in the Paakitsoq region of western Greenland show that decreasing trends in ice velocity during the latter part of the melt season cannot be explained by changes in the ability of moulin-connected channels to convey supraglacial melt. Instead, these observations suggest that decreasing late-season ice velocity may be caused by changes in connectivity in unchannelized regions of the subglacial hydrologic system. Understanding this spatiotemporal variability in subglacial pressures is increasingly important because melt-season dynamics affect ice velocity beyond the conclusion of the melt season8,9,10.

The geological nature of glacier beds plays a key role in ice sheet dynamics. Whereas little is known about Greenland’s subglacial geology, the presence of basal sediments is a necessary condition for fast-flowing Antarctic ice streams. Such sediments sustain subglacial till layers, which if water saturated and under high pore-water pressure provide little resistance to ice flow. Using receiver function modeling of teleseismic P-waves, we report a thick (at least tens of meters) sediment layer beneath a site in Greenland’s ablation zone, ∼15 km away from the western ice sheet margin. Although we do not discuss the origin or detailed nature of these subglacial sediments, we suggest that they are capable of sustaining a till layer. Due to the prevalence of an inefficient, pressurized subglacial drainage system, this till layer would typically be under high pore pressures and fail at the shear stresses transmitted from the overlaying ice. Ice flow resistance is thus focused on regions where till is consolidated or absent, or where form drag over obstacles takes place. In contrast to Antarctica, Greenland’s surface melt now affects practically all latitudes, and efficient hydraulic connections between ice sheet surface and bed are pervasive throughout the ablation zone. This implies that rapid, melt-induced changes in subglacial water pressure are possible. The time required for basal till strength to react to these changes depends on the till’s properties. We estimate that formation and destruction of flow resistance can occur on time scales of less than a few years. This could lead to changes in ice flow that are currently difficult to predict.

Mountain glaciers sample a combination of climate fields – temperature, precipitation and radiation – by accumulation and melting of ice. Flow dynamics acts as a transfer function that maps volume changes to a length response of the glacier terminus. Long histories of terminus positions have been assembled for several glaciers in the Alps. Here I analyze terminus position histories from an ensemble of seven glaciers in the Alps with a macroscopic model of glacier dynamics to derive a history of glacier equilibrium line altitude (ELA) for the time span 400–2010 C.E. The resulting climatic reconstruction depends only on records of glacier variations. The reconstructed ELA history is similar to recent reconstructions of Alpine summer temperature and Atlantic Multidecadal Oscillation (AMO) index, but bears little resemblance to reconstructed precipitation variations. Most reconstructed low-ELA periods coincide with large explosive volcano eruptions, hinting at a direct effect of volcanic radiative cooling on mass balance. The glacier advances during the LIA, and the retreat after 1860, can thus be mainly attributed to temperature and volcanic radiative cooling., The Cryosphere, 8 (2), ISSN:1994-0416, ISSN:1994-0424

Dynamic changes on Greenland outlet glaciers are a primary driver for increases in ice-sheet mass loss and its contribution to sea-level rise. One dramatic example of such change has been observed at Jakobshavn Isbræ, which has thinned, retreated and doubled in speed since the early 2000s. Complementary to large changes on decadal scales, we observe the glacier response on shorter timescales, driven by tidal forcing and calving events. During a 14 day period in August 2009, we documented changes in geometry and speed near the terminus. On this timescale, ice flow responds to forcing at the front from iceberg calving and ocean tides. We observe a step-increase in velocity near the terminus during a large calving event, with transient deceleration in the days following the event. A simple calving-response model explains 94–99% of variations in detrended positions at the four sites considered. During each day, variability due to tidal forcing covers 10–90% of the variability that remains after removing effects accounted for by the calving-response model. The influence of the tidal forcing on flow decays upstream with a characteristic length scale of 2 km, comparable with about two ice thicknesses.

During summer 2011, seismic activity in the ablation zone of the western Greenland ice sheet (GrIS) was monitored using a network of three-component seismometers. The seismic record includes a large variety of icequakes and seismic tremors that demonstrate a clear correlation with subglacial water flow. We verified the existence of well-known shallow icequakes (related to surface crevasse formation), deep icequakes (located at 100–160 m depth) and narrow-banded short-term seismic tremors (tens of seconds in duration). In addition, we present previously unreported long-term tremors lasting several hours. Using attenuation of the measured tremor amplitude, we locate the epicentre of this long-term tremor to a large moulin within our study area. Between 3 and 11 Hz, our continuous seismic record is dominated by this ‘moulin tremor’ and shows strong correlation with the water level of the generating moulin. We argue that monitoring of icequake and glacial tremor sources bears high potential for investigating glacier hydraulics and dynamics, and is thus an ideal supplement to traditional glaciological measurements.

[1] Cold ice within a polythermal ice body controls its flow dynamics through the temperature dependence of viscosity, and affects glacier hydrology by blocking water flow paths. Lakes on the surface, linked by persistent, deeply incised meltwater streams, are hallmark features of cold ice in the ablation zone of a glacier or ice sheet. Ice radar is a convenient method to map scattering from internal water bodies present in ice at the pressure melting temperature (PMT). Consequently, lack of internal scatters is indicative of cold ice. We use a helicopter-borne 30 MHz ice radar to delineate the extent of cold ice within Grenzgletscher (Zermatt, Swiss Alps). The inferred thermal structure is validated with temperature measurements in 15 deep boreholes, showing excellent agreement. The cold ice occupies 80–90 % of the total ice thickness in a 400 m wide flow band along the central flow line. Quantitative interpretation of ice radar scattering power indicates a decrease of ice water content between PMT and 0.5 K below PMT, as predicted by theory, and observed in the laboratory. The cold ice which emerges at the surface in the lower ablation zone is impermeable to water, and is thus devoid of moulins if not crevassed. The surface water from melt and rain is routed through deeply incised, persistent streams and lakes, and cryoconite holes are frequent, in stark contrast to the adjacent temperate ice from other tributaries. The cold ice thus has a strong control on glacier hydrology, but is likely to change due to continued warming.

A range of englacial temperature measurements was acquired in the Monte Rosa area at the border of Switzerland and Italy in the years 1982, 1991, 1994, 1995, 1999, 2000, 2003, 2007 and 2008. Englacial temperatures revealed no evidence of warming at the firn saddle of Colle Gnifetti at 4452 m a.s.l. between 1982 and 1991, the 1991 to 2000 period showed an increase of 0.05 °C per year at a depth of 20 m. From 2000 to 2008 a further increase of 1.3 °C or 0.16 °C per year was observed, indicating that the amount of infiltrating and refreezing meltwater at Colle Gnifetti has probably increased since 2000. The measured temperatures give clear evidence of firn warming since 1991. This is confirmed by five existing boreholes with measured temperature down to bedrock, which were drilled in 1982, 1995, 2003 and 2005. All the observed temperature profiles show a slight bending to warmer temperatures in their uppermost part indicating a warming of the firn, which can be related to the observed atmospheric warming in the 20th century. However, the drilling sites on Colle Gnifetti are still located in the recrystallisation-infiltration zone. A much stronger warming of 6.8 °C or 0.4 °C per year was found at locations beneath Colle Gnifetti on Grenzgletscher from 1991 to 2008. This warming is one order of magnitude greater than the atmospheric warming and can be explained only by a strong increase in the latent heat input by infiltrating and refreezing meltwater. The observations indicate that since 1991, an important firn area beneath Colle Gnifetti has already undergone a firn facies change from the recrystallisation-infiltration to the cold infiltration zone due to an increasing supply of surface melt energy., The Cryosphere, 5 (1), ISSN:1994-0416, ISSN:1994-0424

The length change record of 91 glaciers in the Swiss Alps was analyzed with a novel macroscopic glacier model (LV-model). Based on a history of equilibrium line variations, synthetic length change data were calculated. From the LV-models matching best the measured length changes, characteristic parameters were obtained. The volume time scale thus determined ranges from 5 to 170 years for glaciers of different slope and length. The analysis shows that the observed glacier length changes cannot be reproduced with an equilibrium line variation based on temperature and precipitation alone. The equilibrium line has to be lowered by 100 to 200 meters during several phases of the Little Ice Age (in the time span 1650 to 1850) to obtain observed glacier responses. Such an effect might be attributable to either higher winter precipitation in the Alps, or to radiation forcing.

Recently, increasing concentrations of persistent organic pollutants (POPs) have been observed in the sediment of glacier-fed Lake Oberaar, Switzerland. Melting glaciers have been suggested as a secondary source of POPs released to Alpine lakes. Here we further investigate whether climate warming may accelerate the release of POPs previously deposited to Alpine glaciers ("glacier hypothesis"). To test this hypothesis, a dynamic multimedia mass balance model is developed for the catchment area of Lake Oberaar and is applied to polychlorinated biphenyls (PCBs), polychlorinated dibenzo-p-dioxins and dibenzofurans (PCDD/Fs), and dichlorodiphenyl trichloroethane (DDT). This lake model is combined with two other models. The first is a dynamic multimedia mass balance model parametrized for the Swiss lowlands that is used to calculate (on the basis of historical emission data) the atmospheric concentrations that are an advective input into the model of the lake catchment. The second is a flow model of Oberaar Glacier that determines the residence time of persistent chemicals in the glacier after their deposition to the glacier surface. According to results from these three models in combination, the release of POPs by the glacier is currently increasing and accounts for the observed increase in concentrations in the lake sediment. The models indicate that approximately half of the amount of PCBs, PCDD/Fs, and DDT initially incorporated into the glacier ice is still stored in the glacier. Under the assumption that the climate is warming, accelerated release of POPs is to be expected for the future; in a model run where no climate warming is assumed, the period of time required for release of the same amount of chemicals is longer by several decades than in the scenario with a changing climate.

The transient response of glaciers to climate variations is investigated with a novel low-order model of an idealized glacier resting on uniformly inclined bedrock. The model consists of two volumes representing the accumulation and ablation areas, which are joined by a flux gate controlling mass flow according to the shallow-ice approximation. Under the assumption of a constant vertical mass-balance gradient, a volume–length scaling relation is derived which depends explicitly on mass-balance gradient and bedrock slope. Analytic expressions for the volume and area timescales are given, which are inversely proportional to the mass-balance gradient and a geometric factor which is the ratio between the vertical extent of the ablation area and the ice thickness at the equilibrium line. From the low-order model, a dynamical system in length and volume is obtained. Results from this system are in good agreement with solutions obtained from a transient finite-element model solving the full force-balance and mass-conservation equations. Under periodic forcing there are significant deviations of the length response, which show that the usual relaxation-type parameterization of length change is not well suited for short-term reactions. Switching on and off periodic climate forcing, the model glaciers show surprisingly large initial and final transient responses that have not been investigated before. These results are of significance for the interpretation of length and thickness changes observed on glaciers.

We present the thermal distribution in the confluence area of Gorner- and Grenzgletscher, Valais, Switzerland. The area was mapped by ice-penetrating radar at 1–5 and 40 MHz. The higher-frequency data reveal a thick surface layer of low backscatter in the center of the Grenzgletscher branch. Based on datasets of borehole-temperature measurements and flow velocity, we interpret this as a thick layer of cold ice, advected from the accumulation region of Grenzgletscher. Along seven profiles the base of the low-backscatter zone can be found at a maximum depths between approximately 100 and 200 m. Laterally, the layer extends some 400 m, ∼1/3 of the width of the Grenzgletscher branch. The lower boundary of the low-backscatter zone is systematically higher than the cold–temperate transition surface found in the boreholes. This discrepancy is attributed to the direct sensitivity of radar backscatter to liquid-water inclusions, rather than to the temperature distributions as observed in boreholes. We present the current state of the cold layer and discuss its influence on other glacier characteristics.

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.

Glaciological and related observations from 1961 to 2005 at the summit of Mt Wrangell (62.00° N, 144.02°W; 4317m a.s.l.), a massive glacier-covered shield volcano in south-central Alaska, show marked changes that appear to have been initiated by the Great Alaska Earthquake (Mw = 9.2) of 27 March 1964. The 4×6 km diameter, ice-filled Summit Caldera with several post-caldera craters on its rim, comprises the summit region where annual snow accumulation is 1–2m of water equivalent and the mean annual temperature, measured 10 m below the snow surface, is –20°C. Precision surveying, aerial photogrammetry and measurements of temperature and snow accumulation were used to measure the loss of glacier ice equivalent to about 0.03 km3 of water from the North Crater in a decade. Glacier calorimetry was used to calculate the associated heat flux, which varied within the range 20–140Wm–2; total heat flow was in the range 20–100MW. Seismicity data from the crater’s rim show two distinct responses to large earthquakes at time scales from minutes to months. Chemistry of water and gas from fumaroles indicates a shallow magma heat source and seismicity data are consistent with this interpretation.

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.

Observations of surface motion and ice deformation from 2002–03 were used to infer mean stress fields in a cross-section of Black Rapids Glacier, Alaska, USA, over seasonal timescales. Basal shear stresses in a well-defined zone north of the center line (orographic left) were approximately 7% and 16% lower in spring and summer, respectively, than in winter. Correspondingly higher stresses were found near the margins. These changes in the basal shear stress distribution were sufficiently large to cause mean surface velocities to be 1.2 and 1.5 times larger in spring and summer than in winter. These results were inferred with a simple inverse finite-element flow model that can successfully reproduce bulk surface velocities and tiltmeter data. Stress redistribution between the well-defined zone and the margins may also occur over much shorter time periods as a result of rapidly changing basal conditions (ice–bed decoupling or enhanced till deformation), thereby causing large variations in surface velocity and strongly influencing the glacier’s net motion during summer.

In the frame of the Swiss Alpine Airborne SAR Experiment 2003 (SASARE), multitemporal/multibaseline polarimetric synthetic aperture radar (SAR) interferometry (Pol-InSAR) datasets were acquired with the German Experimental SAR (E-SAR) platform at L- and P-band in the region of the Aletsch glacier in Switzerland. For SASARE, several sub-test-sites had been selected below and significantly above the equilibrium line. We provide results from the analysis of the multibaseline Pol-InSAR L- and P-band measurements of the polarization-dependent complex interferometric coherence of an Alpine glacier in the cold infiltration zone. Interferometric volume decorrelation is demonstrated as a function of polarization. Conventional polarimetric decomposition techniques confirm the polarization-dependent backscattering properties of the analyzed snow and ice surfaces.

Continuum damage mechanics describes the progressive deterioration of material subjected to loading. Jointly used with a level-set method, it proves to be a promising approach to computing the interface motion of a damaged material. For polycrystalline ice, a local isotropic damage evolution law (generalized Kachanow’s law) applied to Glen’s flow law allows the description of tertiary creep and facilitates the modeling of crevasse opening using a failure criterion based on damage accumulation. The use of a level-set method permits the description, in a continuum approach, of the motion of a fractured glacier surface. Using these methods, a model is developed. The ability of this model to describe phenomena connected to crevasse opening is presented. The rupture of a large ice block from a hanging glacier is computed and analyzed. The regular acceleration of such an unstable ice block prior to its collapse is calculated and compared to the acceleration function obtained from observations. A good agreement between the two acceleration functions was found.

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.

Modelling the heat flow in small, cold high-altitude glaciers is important for the interpretation of paleoclimatic data from ice cores. Coupled glacier-flow and heat-flow models are presented that incorporate the densification, heat advection and possible phase transitions at the permafrost boundaries within the bedrock. Marked bends observed in the temperature profiles from two recent boreholes on Colle Gnifetti, Swiss Alps, are interpreted with the help of a transient heat-flow model, driven with a temperature history. The conclusion is that substantial warming of the mean firn temperature at shallow depths has taken place over the last few decades. This has not been observed before in cold-firn regions of the Alps. Modelled heat fluxes in the Monte Rosa massif are strongly influenced by the mountain topography. This leads to a spatial variability of the temperature gradient near the glacier base which has been observed in boreholes to the bedrock. In order to match the measured temperature profiles in the glacier, the vertical heat flux at great depth must be set to an extremely low value. It is shown with the help of the transient heat-flow model that this is a paleoclimatic effect, possibly enhanced by a degrading permafrost base.

A flow model for cold firn proves to be very successful in calculating the age-depth relation of several deep ice cores drilled on Colle Gnifetti, Monte Rosa, Swiss Alps. The compressibility of firn is taken into account by an appropriate constitutive equation, first employed in glaciology by Gagliardini and Meyssonnier (1997), which is implemented in a finite-element code. Flow models of the Colle Gnifetti saddle glaciation in two and three dimensions are based on digital elevation maps of the surface σnd the bedrock, based on radio-echo soundings of the ice thickness. Firn density and the englacial temperature fields are either prescribed or calculated in coupled models. Measured surface velocities, density profiles, the ages of chemically dated layers in ice cores and the closure of a 100 m deep borehole provide benchmarks for the models. The good agreement of modeled and measured quantities confirms that the model includes the relevant physical processes and particularly that the firn flow law is well suited for this type of glacier. The study provides new constraints on the age of the ice near the base as well as the source regions of the ice in the cores.An exceptional flow behavior of the basal ice layer was detected in measurements of borehole closure and inclination. Measurtxl6ed deformation rates exceed upper bounds derived from the flow models, and are thus attributed to altered rheological properties.

Ice deformation and basal motion characterize the dynamical behavior of the Greenland ice sheet (GrIS). We evaluate the contribution of basal motion from ice deformation measurements in boreholes drilled to the bed at two sites in the western marginal zone of the GrIS. We find a sustained high amount of basal motion contribution to surface velocity of 44–73% in winter, and up to 90% in summer. Measured ice deformation rates show an unexpected variation with depth that can be explained with the help of an ice-flow model as a consequence of stress transfer from slippery to sticky areas. This effect necessitates the use of high-order ice-flow models, not only in regions of fast-flowing ice streams but in all temperate-based areas of the GrIS. The agreement between modeled and measured deformation rates confirms that the recommended values of the temperature-dependent flow rate factor A are a good choice for ice-sheet models., Journal of Glaciology, 60 (222), ISSN:0022-1430, ISSN:1727-5652

[1] A novel method to reconstruct glacier volume changes from a measured length record is presented and tested for 13 glaciers in the Swiss Alps. The response of a glacier to changes in climate is modeled with a two-parameter dynamical system in the variables “length” and “volume”. Driven by a history of equilibrium line altitude (ELA), the model yields variations of glacier length and volume. A dynamically equivalent simple model (DESM) is determined for each glacier by matching modeled and measured length changes. The volume changes predicted with the DESM agree well with measurements for 12 glaciers, whereas agreement is poor for one glacier with topographic breaks in the terminus area. For all glaciers, which are located in different climate regions, the length and volume changes are reproduced with the same ELA history. This agreement shows that the macroscopic glacier response to the climate history is well correlated over a whole mountain range. Modeling the future evolution of the glaciers under a constant present-day climate reveals that fast-reacting glaciers are close to equilibrium, whereas length and volume of the large valley glaciers would be reduced during the next century by an amount similar to the volume lost during the last 150 years.

[1] We used time-lapse imagery, seismic and audio recordings, iceberg and glacier velocities, ocean wave measurements, and simple theoretical considerations to investigate the interactions between Jakobshavn Isbrae and its proglacial ice melange. The melange behaves as a weak, granular ice shelf whose rheology varies seasonally. Sea ice growth in winter stiffens the melange matrix by binding iceberg clasts together, ultimately preventing the calving of full-glacier-thickness icebergs (the dominant style of calving) and enabling a several kilometer terminus advance. Each summer the melange weakens and the terminus retreats. The melange remains strong enough, however, to be largely unaffected by ocean currents (except during calving events) and to influence the timing and sequence of calving events. Furthermore, motion of the melange is highly episodic: between calving events, including the entire winter, it is pushed down fjord by the advancing terminus (at ∼40 m d−1), whereas during calving events it can move in excess of 50 × 103 m d−1 for more than 10 min. By influencing the timing of calving events, the melange contributes to the glacier's several kilometer seasonal advance and retreat; the associated geometric changes of the terminus area affect glacier flow. Furthermore, a force balance analysis shows that large-scale calving is only possible from a terminus that is near floatation, especially in the presence of a resistive ice melange. The net annual retreat of the glacier is therefore limited by its proximity to floatation, potentially providing a physical mechanism for a previously described near-floatation criterion for calving.

[1] The recent loss of Jakobshavn Isbrae's extensive floating ice tongue has been accompanied by a change in near terminus behavior. Calving currently occurs primarily in summer from a grounded terminus, involves the detachment and overturning of several icebergs within 30–60 min, and produces long-lasting and far-reaching ocean waves and seismic signals, including “glacial earthquakes”. Calving also increases near-terminus glacier velocities by ∼3% but does not cause episodic rapid glacier slip, thereby contradicting the originally proposed glacial earthquake mechanism. We propose that the earthquakes are instead caused by icebergs scraping the fjord bottom during calving.

The supply, storage, and flow of water under an ice sheet are crucial for its overall dynamics. The density difference between ice and water limits subglacial water pressure to the pressure from overlying ice. Near this limit, the ice gradually becomes decoupled from the base (see the figure), water-filled cavities form, and sediments become weakened. The velocity of ice moving over the bedrock increases due to reduced friction and faster sediment deformation. Subglacial water pressure is therefore the key parameter that controls basal motion, ice flux, and the future evolution of an ice sheet. On page 777 of this issue ( 1 ), Meierbachtol et al. present a detailed record of water pressure variability under one region of the Greenland ice sheet.

[1] During summer 2001, a 140 m long ice core was recovered from the Belukha glacier (49°48′26″N, 86°34′43″E, 4062 m a.s.l.) in the Siberian Altai. The englacial temperature of −17°C indicates that this unique glaciochemical record is well preserved and suitable for the reconstruction of air pollution levels in this previously unexplored region. The upper 86 m were dated by 210Pb and cover the period 1862–2001. A lack of strong winter minima was observed in the δ18O record and attributed to the small amount of precipitation during that season and to wind erosion. The ion concentrations are comparable to those observed in Swiss glaciers, except for ammonium and formate, where enhanced concentrations indicate biogenic emissions from Siberian forests. Sulfate, ammonium and nitrate records all show anthropogenic impacts despite the remoteness of this site.

[1] A zone of highly irregular ice motion has been detected within the ice sheet adjacent to Jakobshavn Isbrae, a fast moving ice stream in West Greenland. Two tilt sensors in different boreholes and frozen into the ice show angle variations of up to 0.1° and 0.3° within a 1200 s interval. The amplitude of these flutter oscillations increases toward the transition from Holocene to Wisconsin ice at 682 m depth, 150 m above bedrock. All other tilt sensors above and below show low background noise. A numerical model confirms that the stress configuration is exceptional in the neighborhood of the ice stream. Stress transfer from the ice stream to the margins causes an extra forcing of the basal ice at the drill site. Several possible explanations for the observed flutter oscillations are given, the most likely being dynamic rupture events or stick-slip motion along the Holocene-Wisconsin transition. If true, this would correspond to an active overthrust fault with very high loading rates.

Sliding of the Greenland ice sheet is affected by the production of surface meltwater. A new theory shows that whether the result is a long-term speed-up or slow-down of ice motion depends on the variability in melt input. See Letter p.803 Recent observations in Greenland have led to the notion that surface meltwater supply to the ice-sheet bed lubricates ice flow, suggesting that climate warming could lead to runaway glacial acceleration. Christian Schoof uses a physically based model that captures drainage channelization under the ice to challenge this view. The model shows that increased melt supply leads to channelization and a drop in water pressure. This causes ice flow to slow down rather than speed up. However, this effect can be overcome if melt supply is variable over short timescales, when temporary pressure spikes can lead to accelerated flow. The positive melt/dynamic thinning feedback is therefore still viable, but is heavily dependent on daily temperature variations and rain events that are largely ignored in current ice sheet models.

At a site on the ice sheet adjacent to the Jakobshavn ice stream in West Greenland, ice deformation rates and temperatures have been measured in boreholes to the bedrock at 830 m depth. Enhanced deformation rates were recorded just below the Holocene–Wisconsin transition at 680 m depth. A 31 m layer of temperate ice and the temperature minimum of −22°C at 520 m depth were detected. The good agreement of these data with results of a two-dimensional thermomechanically coupled flow model implies that the model input is adequate. Discrepancies between modelled and measured temperature profiles on a flowline at the ice-stream centre have been attributed to effects not accounted for by the model. We have suggested that the convergent three-dimensional flow leads to a vertical extension of the basal ice entering the stream. A thick basal layer of temperate and Wisconsin ice would explain the fast flow of this ice stream. As a test of this hypothesis, the new core-borehole conductivity (CBC) method has been used to compare conductivity sequences from the ice stream to those of the adjacent ice sheet. The correlation thus inferred suggests that the lowest 270 m of the ice sheet correspond to the lowermost 1700 m of the stream, and, consequently, that the lower part of the ice stream has experienced a very large vertical extension.