Search

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.

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

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.

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

The puzzling fact that Jakobshavn Isbræ, West Greenland, is flowing very fast but without any significant seasonal velocity changes, despite big amounts of surface-derived meltwater entering the ice stream (Echelmeyer and Harrison, 1990), has been explained by a combination of different types of measurements. It is now well established from seismic measurements and radio-echo sounding that Jakobshavn Isbræ flows through a deeply eroded subglacial trench that, even 50 km inland of the grounding line, extends as far as 1500 m below sea level (Clarke and Echelmeyer, 1996; Legarsky and Huang, 2006). Temperature measurements in boreholes down to 65% of the 2500 m thick ice stream at site B, some 50 km upstream of the calving front (Fig. 1b; Iken and others, 1993), were used to infer the presence of a substantial layer of temperate ice, the thickness of which was estimated to be at least 300 m by modeling and matching internal layering structures (Funk and others, 1994; Lüthi and others, 2002). The presence of a thick layer of temperate ice under very high driving stress allows for high ice-deformation rates, which contribute substantially to the observed fast flow velocities. Basal motion, while certainly important, seems to be barely influenced by the seasonal meltwater input (Echelmeyer and Harrison, 1990).

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.

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.

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.

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.

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

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.