Basal Crevasses Research Articles

Sensitivity of a calving glacier to ice–ocean interactions under climate change: new insights from a 3-D full-Stokes model

Abstract. Iceberg calving accounts for between 30 % and 60 % of net mass loss from the Greenland Ice Sheet, which has intensified and is now the single largest contributor to global sea level rise in the cryosphere. Changes to calving rates and the dynamics of calving glaciers represent a significant uncertainty in projections of future sea level rise. A growing body of observational evidence suggests that calving glaciers respond rapidly to regional environmental change, but predictive capacity is limited by the lack of suitable models capable of simulating calving mechanisms realistically. Here, we use a 3-D full-Stokes calving model to investigate the environmental sensitivity of Store Glacier, a large outlet glacier in West Greenland. We focus on two environmental processes: undercutting by submarine melting and buttressing by ice mélange, and our results indicate that Store Glacier is likely to be able to withstand moderate warming perturbations in which the former is increased by 50 % and the latter reduced by 50 %. However, severe perturbation with a doubling of submarine melt rates or a complete loss of ice mélange destabilises the calving front in our model runs. Furthermore, our analysis reveals that stress and fracture patterns at Store's terminus are complex and varied, primarily due to the influence of basal topography. Calving style and environmental sensitivity vary greatly, with propagation of surface crevasses significantly influencing iceberg production in the northern side, whereas basal crevasses dominate in the south. Any future retreat is likely to be initiated in the southern side by a combination of increased submarine melt rates in summer and reduced mélange strength in winter. The lateral variability, as well as the importance of rotational and bending forces at the terminus, underlines the importance of using the 3-D full-Stokes stress solution when modelling Greenland's calving glaciers.

Seismic Noise Interferometry Reveals Transverse Drainage Configuration Beneath the Surging Bering Glacier

AbstractSubglacial drainage systems are known to critically control ice flows, but their spatial configuration and temporal evolution are poorly constrained due to inaccessibility. Here we report a 12‐year‐long monitoring of the drainage underneath Bering Glacier, Alaska, by correlating ambient noise recorded at two seismic stations on the sides of the glacier. We find that the seismic surface waves traveling across Bering Glacier slowed down by 1–2% during its latest 2008–2011 surge, likely due to the switch of the subglacial drainage from a channelized system to a distributed system. In contrast to current models, the relative amplitude of velocity reductions for Rayleigh and Love waves requires the distributed drainage to be highly anisotropic and aligned perpendicular to the ice flow direction. We infer that the subglacial water flow is mainly through a network of transverse basal crevasses during surges and thus can sustain the high water pressure and ice flow speed.

On the evaluation of the stress intensity factor in calving models using linear elastic fracture mechanics

ABSTRACTWe investigate the appropriateness of calving or crevasse models from the literature using linear elastic fracture mechanics (LEFM). To this end, we compare LEFM model-predicted stress intensity factors (SIFs) against numerically computed SIFs using the displacement correlation method in conjunction with the finite element method. We present several benchmark simulations wherein we calculate the SIF at the tips of water-filled surface and basal crevasses penetrating through rectangular ice slabs under different boundary conditions, including grounded and floating conditions. Our simulation results indicate that the basal boundary condition significantly influences the SIF at the crevasse tips. We find that the existing calving models using LEFM are not generally accurate for evaluating SIFs in grounded glaciers or floating ice shelves. We also illustrate that using the ‘single edge crack’ weight function in the LEFM formulations may be appropriate for predicting calving from floating ice shelves, owing to the low fracture toughness of ice; whereas, using the ‘double edge crack’ or ‘central through crack’ weight functions is more appropriate for predicting calving from grounded glaciers. To conclude, we recommend using the displacement correlation method for SIF evaluation in real glaciers and ice shelves with complex geometries and boundary conditions.

Modelled fracture and calving on the Totten Ice Shelf

Abstract. The Totten Ice Shelf (IS) has a large drainage basin, much of which is grounded below sea level, leaving the glacier vulnerable to retreat through the marine ice sheet instability mechanism. The ice shelf has also been shown to be sensitive to changes in calving rate, as a very small retreat of the calving front from its current position is predicted to cause a change in flow at the grounding line. Therefore understanding the processes behind calving on the Totten IS is key to predicting its future sea level rise contribution. Here we use the Helsinki Discrete Element Model (HiDEM) to show that not all of the fractures visible at the front of the Totten IS are produced locally, but that the across-flow basal crevasses, which are part of the distinctive cross-cutting fracture pattern, are advected into the calving front area from upstream. A separate simulation of the grounding line shows that re-grounding points may be key areas of basal crevasse production, and can produce basal crevasses in both an along- and across-flow orientation. The along-flow basal crevasses at the grounding line may be a possible precursor to basal channels, while we suggest the across-flow grounding-line fractures are the source of the across-flow features observed at the calving front. We use two additional models to simulate the evolution of basal fractures as they advect downstream, demonstrating that both strain and ocean melt have the potential to deform narrow fractures into the broad basal features observed near the calving front. The wide range of factors which influence fracture patterns and calving on this glacier will be a challenge for predicting its future mass loss.

A Full‐Stokes 3‐D Calving Model Applied to a Large Greenlandic Glacier

AbstractIceberg calving accounts for around half of all mass loss from both the Greenland and Antarctic ice sheets. The diverse nature of calving and its complex links to both internal dynamics and climate make it challenging to incorporate into models of glaciers and ice sheets. Here we present results from a new open‐source 3‐D full‐Stokes calving model developed in Elmer/Ice. The calving model implements the crevasse depth criterion, which states that calving occurs when surface and basal crevasses penetrate the full thickness of the glacier. The model also implements a new 3‐D rediscretization approach and a time‐evolution scheme which allow the calving front to evolve realistically through time. We test the model in an application to Store Glacier, one of the largest outlet glaciers in West Greenland, and find that it realistically simulates the seasonal advance and retreat when two principal environmental forcings are applied. These forcings are (1) submarine melting in distributed and concentrated forms and (2) ice mélange buttressing. We find that ice mélange buttressing is primarily responsible for Store Glacier's seasonal advance and retreat. Distributed submarine melting prevents the glacier from forming a permanent floating tongue, while concentrated plume melting has a disproportionately large and potentially destabilizing effect on the calving front position. Our results also highlight the importance of basal topography, which exerts a strong control on calving, explaining why Store Glacier has remained stable during a period when neighboring glaciers have undergone prolonged interannual retreat.

On the interpretation of ice-shelf flexure measurements

ABSTRACTTidal flexure in ice shelf grounding zones has been used extensively in the past to determine grounding line position and ice properties. Although the rheology of ice is viscoelastic at tidal loading frequencies, most modelling studies have assumed some form of linear elastic beam approximation to match observed flexure profiles. Here we use density, radar and DInSAR measurements in combination with full-Stokes viscoelastic modelling to investigate a range of additional controls on the flexure of the Southern McMurdo Ice Shelf. We find that inclusion of observed basal crevasses and density dependent ice stiffness can greatly alter the flexure profile and yet fitting a simple elastic beam model to that profile will still produce an excellent fit. Estimates of the effective Young's modulus derived by fitting flexure profiles are shown to vary by over 200% depending on whether these factors are included, even when the local thickness is well constrained. Conversely, estimates of the grounding line position are relatively insensitive to these considerations for the case of a steep bed slope in our study region. By fitting tidal amplitudes only, and ignoring phase information, elastic beam theory can provide a good fit to observations in a wide variety of situations. This should, however, not be taken as an indication that the underlying rheological assumptions are correct.

An englacial hydrologic system of brine within a cold glacier: Blood Falls, McMurdo Dry Valleys, Antarctica

ABSTRACTTaylor Glacier hosts an active englacial hydrologic system that feeds Blood Falls, a supraglacial outflow of iron-rich subglacial brine at the terminus, despite mean annual air temperatures of −17°C and limited surface melt. Taylor Glacier is an outlet glacier of the East Antarctic ice sheet that terminates in Lake Bonney, McMurdo Dry Valleys. To image and map the brine feeding Blood Falls, we used radio echo sounding to delineate a subhorizontal zone of englacial brine upstream from Blood Falls and elongated in the ice flow direction. We estimate volumetric brine content in excess of 13% within 2 m of the central axis of this zone, and likely much higher at its center. Brine content decreases, but remains detectable, up to 45 m away along some transects. Hence, we infer a network of subparallel basal crevasses allowing injection of pressurized subglacial brine into the ice. Subglacial brine is routed towards Blood Falls by hydraulic potential gradients associated with deeply incised supraglacial valleys. The brine remains liquid within the subglacial and englacial environments through latent heat of freezing coupled with elevated salt content. Our findings suggest that cold glaciers could support freshwater hydrologic systems through localized warming by latent heat alone.

Bounds on the calving cliff height of marine terminating glaciers

AbstractIncreased calving and rapid retreat of glaciers can contribute significantly to sea level rise, but the processes controlling glacier retreat remain poorly understood. We seek to improve our understanding of calving by investigating the stress field controlling tensile and shear failure using a 2‐D full‐Stokes finite element model. Using idealized rectangular geometries, we find that when rapidly sliding glaciers thin to near buoyancy, full thickness tensile failure occurs, similar to observations motivating height‐above‐buoyancy calving laws. In contrast, when glaciers are frozen to their beds, basal crevasse penetration is suppressed and calving is minimal. We also find that shear stresses are largest when glaciers are thickest. Together, the tensile and shear failure criteria map out a stable envelope in an ice‐thickness‐water‐depth diagram. The upper and lower bounds on cliff height can be incorporated into numerical ice sheet models as boundary conditions, thus bracketing the magnitude of calving rates in marine‐terminating glaciers.

The Bothnian Sea ice stream: early Holocene retreat dynamics of the south‐central Fennoscandian Ice Sheet

The Gulf of Bothnia hosted a variety of palaeo‐glaciodynamic environments throughout the growth and decay of the last Fennoscandian Ice Sheet, from the main ice‐sheet divide to a major corridor of marine‐ and lacustrine‐based deglaciation. Ice streaming through the Bothnian and Baltic basins has been widely assumed, and the damming and drainage of the huge proglacial Baltic Ice Lake has been implicated in major regional and hemispheric climate changes. However, the dynamics of palaeo‐ice flow and retreat in this large marine sector have until now been inferred only indirectly, from terrestrial, peripheral evidence. Recent acquisition of high‐resolution multibeam bathymetry opens these basins up, for the first time, to direct investigation of their glacial footprint and palaeo‐ice sheet behaviour. Here we report on a rich glacial landform record: in particular, a palaeo‐ice stream pathway, abundant traces of high subglacial meltwater volumes, and widespread basal crevasse squeeze ridges. The Bothnian Sea ice stream is a narrow flow corridor that was directed southward through the basin to a terminal zone in the south‐central Bothnian Sea. It was activated after initial margin retreat across the Åland sill and into the Bothnian basin, and the exclusive association of the ice‐stream pathway with crevasse squeeze ridges leads us to interpret a short‐lived stream event, under high extension, followed by rapid crevasse‐triggered break‐up. We link this event with a c. 150‐year ice‐rafted debris signal in peripheral varved records, at c. 10.67 cal. ka BP. Furthermore, the extensive glacifluvial system throughout the Bothnian Sea calls for considerable input of surface meltwater. We interpret strongly atmospherically driven retreat of this marine‐based ice‐sheet sector.

Accelerated ice shelf rifting and retreat at Pine Island Glacier, West Antarctica

AbstractPine Island Glacier has undergone several major iceberg calving events over the past decades. These typically occurred when a rift at the heavily fractured shear margin propagated across the width of the ice shelf. This type of calving is common on polar ice shelves, with no clear connection to ocean‐ice dynamic forcing. In contrast, we report on the recent development of multiple rifts initiating from basal crevasses in the center of the ice shelf, resulted in calving further upglacier than previously observed. Coincident with rift formation was the sudden disintegration of the ice mélange that filled the northern shear margin, resulting in ice sheet detachment from this margin. Examination of ice velocity suggests that this internal rifting resulted from the combination of a change in ice shelf stress regime caused by disintegration of the mélange and intensified melting within basal crevasses, both of which may be linked to ocean forcing.

Modeling hydraulic fracture of glaciers using continuum damage mechanics

ABSTRACTThe presence of water-filled crevasses is known to increase the penetration depth of crevasses and this has been hypothesized to play an important role controlling iceberg calving rate. Here, we develop a continuum-damage-based poro-mechanics formulation that enables the simulation of water-filled basal and surface crevasse propagation. The formulation incorporates a scalar isotropic damage variable into a Maxwell-type viscoelastic constitutive model for glacial ice, and the effect of the water pressure on fracture propagation using the concept of effective solid stress. We illustrate the model by simulating quasi-static hydrofracture in idealized rectangular slabs of ice in contact with the ocean. Our results indicate that water-filled basal crevasses only propagate when the water pressure is sufficiently large, and that the interaction between simultaneously propagating water-filled surface and basal crevasses can have a mutually positive influence leading to deeper crevasse propagation, which can critically affect glacial stability. Therefore, this study supports the hypothesis that hydraulic fracture is a plausible mechanism for the accelerated breakdown of glaciers.

Rhombohedral crevasse-fill ridges at the marine margin of a surging Svalbard ice cap

The ice cap of Austfonna in eastern Svalbard is the largest in the Eurasian Arctic at 8000 km2 and has about 200 km of marine-terminating ice cliffs (Dowdeswell et al. 2008). Several of its drainage basins are of surge type (Meier & Post 1969) and, between 1936 and 1938, one of these basins (Brasvellbreen; 1100 km2) increased its velocity rapidly and underwent an advance of about 20 km along a 30 km wide front (Schytt 1969). Since that time the ice-cap terminus has stagnated and retreated, losing mass by a combination of surface melting, thinning and iceberg production. Retreat has revealed several distinctive and well-preserved submarine landforms (Fig. 1) linked to this recent surge activity (Solheim & Pfirman 1985; Solheim 1991). Fig. 1. ( a ) Rhombohedral ridges and a terminal moraine ridge on the seafloor offshore of Brasvellbreen, an outlet of the Austfonna ice cap in eastern Svalbard. Image courtesy of the Norwegian Hydrographic Service (Permission no. 14/G754). Acquisition system Kongsberg EM1002. Frequency 97 kHz. Grid-cell size 3 m. The ridges are interpreted to have been formed by the filling of basal crevasses with deforming diamictic sediment. ( b ) Oblique aerial photograph of the surge and rapid advance of Brasvellbreen in 1938 (Photo S38-1958; courtesy of Norsk Polarinstitutt). ( c ) Location of study area (red box; map from IBCAO v. 3.0). ( d …

Debris entrainment and landform genesis during tidewater glacier surges

AbstractThe englacial entrainment of basal debris during surges presents an opportunity to investigate processes acting at the glacier bed. The subsequent melt‐out of debris‐rich englacial structures during the quiescent phase produces geometrical ridge networks on glacier forelands that are diagnostic of surge activity. We investigate the link between debris entrainment and proglacial geomorphology by analyzing basal ice, englacial structures, and ridge networks exposed at the margins of Tunabreen, a tidewater surge‐type glacier in Svalbard. The basal ice facies display clear evidence for brittle and ductile tectonic deformation, resulting in overall thickening of the basal ice sequence. The formation of debris‐poor dispersed facies ice is the result of strain‐induced metamorphism of meteoric ice near the bed. Debris‐rich englacial structures display a variety of characteristics and morphologies and are interpreted to represent the incorporation and elevation of subglacial till via the squeezing of till into basal crevasses and hydrofracture exploitation of thrust faults, reoriented crevasse squeezes, and preexisting fractures. These structures are observed to melt‐out and form embryonic geometrical ridge networks at the base of a terrestrially grounded ice cliff. Ridge networks are also located at the terrestrial margins of Tunabreen, neighboring Von Postbreen, and in a submarine position within Tempelfjorden. Analysis of network characteristics allows these ridges to be linked to different formational mechanisms of their parent debris‐rich englacial structures. This in turn provides an insight into variations in the dominant tectonic stress regimes acting across the glacier during surges.

Dynamics of glacier calving at the ungrounded margin of Helheim Glacier, southeast Greenland.

During summer 2013 we installed a network of 19 GPS nodes at the ungrounded margin of Helheim Glacier in southeast Greenland together with three cameras to study iceberg calving mechanisms. The network collected data at rates up to every 7 s and was designed to be robust to loss of nodes as the glacier calved. Data collection covered 55 days, and many nodes survived in locations right at the glacier front to the time of iceberg calving. The observations included a number of significant calving events, and as a consequence the glacier retreated ~1.5 km. The data provide real‐time, high‐frequency observations in unprecedented proximity to the calving front. The glacier calved by a process of buoyancy‐force‐induced crevassing in which the ice downglacier of flexion zones rotates upward because it is out of buoyant equilibrium. Calving then occurs back to the flexion zone. This calving process provides a compelling and complete explanation for the data. Tracking of oblique camera images allows identification and characterisation of the flexion zones and their propagation downglacier. Interpretation of the GPS data and camera data in combination allows us to place constraints on the height of the basal cavity that forms beneath the rotating ice downglacier of the flexion zone before calving. The flexion zones are probably formed by the exploitation of basal crevasses, and theoretical considerations suggest that their propagation is strongly enhanced when the glacier base is deeper than buoyant equilibrium. Thus, this calving mechanism is likely to dominate whenever such geometry occurs and is of increasing importance in Greenland.

Temperature distribution and thermal anomalies along a flowline of the Greenland ice sheet

AbstractEnglacial and basal temperature data for the Greenland ice sheet (GrIS) are sparse and mostly limited to deep interior sites and ice streams, providing an incomplete representation of the thermal state of ice within the ablation zone. Here we present 11 temperature profiles at five sites along a 34 km east–west transect of West Greenland. These profiles depict ice temperatures along a flowline and local temperature variations between closely spaced boreholes. A temperate basal layer is present in all profiles, increasing in thickness in the flow direction, where it expands from ∼3% of ice height furthest inland to 100% at the margin. Temperate thickness growth is inconsistent with modeled heat contributions from strain heating, heat conduction, and vertical extension of the temperate layer. We suggest that basal crevassing, facilitated by water pressures at or near ice overburden pressure, is responsible for the large temperate ice thicknesses observed. High-temperature kinks at 51–85 m depth are likely remnants from the thermal influence of partially water-filled crevasses up ice sheet. Steep horizontal temperature gradients between closely grouped boreholes suggest the recent thermal influence of a moulin. These profiles demonstrate the ability of meltwater to rapidly alter ice temperatures at all depths within the ablation zone.

A model of viscoelastic ice-shelf flexure

AbstractWe develop a formal thin-plate treatment of the viscoelastic flexure of floating ice shelves as an initial step in treating various problems relevant to ice-shelf response to sudden changes of surface loads and applied bending moments (e.g. draining supraglacial lakes, iceberg calving, surface and basal crevassing). Our analysis is based on the assumption that total deformation is the sum of elastic and viscous (or power-law creep) deformations (i.e. akin to a Maxwell model of viscoelasticity, having a spring and dashpot in series). The treatment follows the assumptions of well-known thin-plate approximation, but is presented in a manner familiar to glaciologists and with Glen’s flow law. We present an analysis of the viscoelastic evolution of an ice shelf subject to a filling and draining supraglacial lake. This demonstration is motivated by the proposition that flexure in response to the filling/drainage of meltwater features on the Larsen B ice shelf, Antarctica, contributed to the fragmentation process that accompanied its collapse in 2002.

Evolution of basal crevasses links ice shelf stability to ocean forcing

Iceberg calving is one of the primary mechanisms responsible for transferring ice from the Antarctic ice shelves to the ocean, but remains poorly understood. Previous theories of calving have sought to explain the calving process as a brittle phenomenon that occurs rapidly when surface or basal crevasses penetrate the entire ice thickness. Here we show that the strain-rate-weakening nature of ice permits initially narrow basal crevasses to seed an instability that gives rise to locally enhanced ductile deformation and thinning over length scales that are large compared to the ice thickness. This ductile failure process, called necking, amplifies long wavelength features of bottom topography and allows basal crevasses to penetrate an increasing fraction of the ice thickness as they advect downstream. Application of the model to the four largest Antarctic ice shelves shows that necking occurs downstream of pinning points and sharp protrusions in the ice shelf embayment where stress is highly concentrated. However, model predictions are sensitive to assumptions about basal melting and refreezing within crevasses, suggesting that the combination of mechanical instability and ice–ocean interaction on the scale of an individual crevasse may play a leading role in controlling ice shelf stability.

Buoyant flexure and basal crevassing in dynamic mass loss at Helheim Glacier

Despite the role that calving plays in Greenland mass loss, the mechanisms of calving are poorly constrained. Observations of Greenland’s Helheim Glacier suggest that buoyant flexure at the glacier terminus leads to the propagation of basal crevasses and iceberg calving. Iceberg calving accounts for a significant proportion of annual mass loss from marine-terminating glaciers1,2 and may have been a factor in the rapid demise of ancient ice sheets3. The largest contributions from the main outlet glaciers of the Greenland ice sheet to sea-level rise over the next two centuries have been projected to be dynamic in origin, that is, driven by glacier flow and calving4. However, present physical models remain a coarse approximation of real calving mechanisms because models are poorly constrained by sparse glacier geometry observations5. Here we present a record of daily digital elevation models from the calving margin of Greenland’s Helheim Glacier at a high spatial resolution. Our digital elevation models are derived from stereo terrestrial photography taken over the summers of 2010 and 2011. We find that during these two summers dynamic mass loss at Helheim Glacier was dominated by calving events exceeding 1 km3 that were the result of buoyant flexure and the propagation of basal crevasses. We suggest that this buoyancy-driven mechanism for calving may be common elsewhere in Greenland and could be a first-order control on the ice sheet’s future contribution to sea-level rise.

Oscillatory subglacial drainage in the absence of surface melt

Abstract. The presence of strong diurnal cycling in basal water pressure records obtained during the melt season is well established for many glaciers. The behaviour of the drainage system outside the melt season is less well understood. Here we present borehole observations from a surge-type valley glacier in the St Elias Mountains, Yukon Territory, Canada. Our data indicate the onset of strongly correlated multi-day oscillations in water pressure in multiple boreholes straddling a main drainage axis, starting several weeks after the disappearance of a dominant diurnal mode in August 2011 and persisting until at least January 2012, when multiple data loggers suffered power failure. Jökulhlaups provide a template for understanding spontaneous water pressure oscillations not driven by external supply variability. Using a subglacial drainage model, we show that water pressure oscillations can also be driven on a much smaller scale by the interaction between conduit growth and distributed water storage in smaller water pockets, basal crevasses and moulins, and that oscillations can be triggered when water supply drops below a critical value. We suggest this in combination with a steady background supply of water from ground water or englacial drainage as a possible explanation for the observed wintertime pressure oscillations.

Modeling ice‐ocean interaction in ice‐shelf crevasses

AbstractOcean freezing within ice‐shelf basal crevasses could potentially act as a stabilizing influence on ice shelves; however, ice‐ocean interaction and ocean dynamics within these crevasses are as yet poorly understood. To this end, an idealized 2‐D model of an ice‐shelf basal crevasse has been developed using Fluidity, a finite‐element ocean model using an unstructured mesh. A simple model of frazil ice formation and deposition has been incorporated into Fluidity to better represent the freezing process. Model results show two different flow regimes, dependent on the amount of freezing in the crevasse: one driven by freezing at the top of the crevasse and the other by the ingress of meltwater from outside the crevasse. In the first, freezing at the top of the crevasse leads to the formation of an unstable overturning circulation due to the rejection of dense, salty water. In the second, a buoyant layer is formed along the sides and roof of the crevasse, stratifying the water column. Frazil ice precipitation is found to be the dominant freezing process at the top of the basal crevasse in the freeze‐driven case, with direct freezing being dominant in the melt‐driven case. In both cases, melting occurs lower down on the walls of the crevasse due to the strong overturning circulation. The freezing in ice‐shelf crevasses and rifts is found to be highly dependent upon ocean temperature, providing a stabilizing influence on ice shelves underlain by cold waters that is not present elsewhere.