Influence of anisotropy on velocity and age distribution at Scharffenbergbotnen blue ice area

We use a full-Stokes thermo-mechanically coupled ice-flow model to study the dynamics of the glacier inside Scharffenbergbotnen valley, Dronning Maud Land, Antarctica. The domain encompasses a high accumulation rate region and, downstream a sublimation-dominated bare ice ablation area. The ablation ice area is notable for having old ice at its surface since the vertical velocity is upwards, and horizontal velocities are almost stagnant there. We compare the model simulation with field observations of velocities and the age distribution of the surface ice. A satisfactory match with simulations using an isotropic flow law was not found because of too high horizontal velocities and too slow vertical ones. However, the existence of a pronounced ice fabric may explain the present day surface velocity distribution in the inner Scharffenbergbotnen blue ice area. Near absence of data on the temporal evolution of Scharffenbergbotnen since the Late Glacial Maximum necessitates exploration of the impact of anisotropy using prescribed ice fabrics: isotropic, single maximum, and linear variation with depth, in both two-dimensional and three dimensional flow models. The realistic velocity field simulated with a non-collinear orthotropic flow law, however produced surface ages in significant disagreement with the few reliable age measurements and suggests that the age field is not in a steady state and that the present distribution is a result of a flow reorganization at about 15 000 yr BP. In order to fully understand the surface age distribution a transient simulation starting from the Late Glacial Maximum including the correct initial conditions for geometry, age, fabric and temperature distribution would be needed. It is the first time that the importance of anisotropy has been demonstrated in the ice dynamics of a blue ice area. This is useful to understand ice flow in order to better interpret archives of ancient ice for paleoclimate research

A double continuum hydrological model for glacier applications

Peer reviewedPublisher PD

A three-dimensional full Stokes model of the grounding line dynamics: effect of a pinning point beneath the ice shelf

The West Antarctic ice sheet is confined by a large area of ice shelves, fed by inland ice through fast flowing ice streams. The dynamics of the grounding line, which is the line-boundary between grounded ice and the downstream ice shelf, has a major influence on the dynamics of the whole ice sheet. However, most ice sheet models use simplifications of the flow equations, as they do not include all the stress components, and are known to fail in their representation of the grounding line dynamics. Here, we present a 3-D full Stokes model of a marine ice sheet, in which the flow problem is coupled with the evolution of the upper and lower free surfaces, and the position of the grounding line is determined by solving a contact problem between the shelf/sheet lower surface and the bedrock. Simulations are performed using the open-source finite-element code Elmer/Ice within a parallel environment. The model's ability to cope with a curved grounding line and the effect of a pinning point beneath the ice shelf are investigated through prognostic simulations. Starting from a steady state, the sea level is slightly decreased to create a contact point between a seamount and the ice shelf. The model predicts a dramatic decrease of the shelf velocities, leading to an advance of the grounding line until both grounded zones merge together, during which an ice rumple forms above the contact area at the pinning point. Finally, we show that once the contact is created, increasing the sea level to its initial value does not release the pinning point and has no effect on the ice dynamics, indicating a stabilising effect of pinning points

Modelling environmental influences on calving at Helheim Glacier in eastern Greenland

Calving is an important mass-loss process for many glaciers worldwide, and has been assumed to respond to a variety of environmental influences. We present a grounded, flowline tidewater glacier model using a physically-based calving mechanism, applied to Helheim Glacier, eastern Greenland. By qualitatively examining both modelled size and frequency of calving events, and the subsequent dynamic response, the model is found to realistically reproduce key aspects of observed calving behaviour. Experiments explore four environmental variables which have been suggested to affect calving rates: water depth in crevasses, basal water pressure, undercutting of the calving face by submarine melt and backstress from ice mélange. Of the four variables, only crevasse water depth and basal water pressure were found to have a significant effect on terminus behaviour when applied at a realistic magnitude. These results are in contrast to previous modelling studies, which have suggested that ocean temperatures could strongly influence the calving front. The results raise the possibility that Greenland outlet glaciers could respond to the recent trend of increased surface melt observed in Greenland more strongly than previously thought, as surface ablation can strongly affect water depth in crevasses and water pressure at the glacier bed

Immigració i integració

Abstract not availabl

Effects of undercutting and sliding on calving: a global approach applied to Kronebreen, Svalbard

In this paper, we study the effects of basal friction, sub-aqueous undercutting and glacier geometry on the calving process by combining six different models in an offline-coupled workflow: a continuum–mechanical ice flow model (Elmer/Ice), a climatic mass balance model, a simple subglacial hydrology model, a plume model, an undercutting model and a discrete particle model to investigate fracture dynamics (Helsinki Discrete Element Model, HiDEM). We demonstrate the feasibility of reproducing the observed calving retreat at the front of Kronebreen, a tidewater glacier in Svalbard, during a melt season by using the output from the first five models as input to HiDEM. Basal sliding and glacier motion are addressed using Elmer/Ice, while calving is modelled by HiDEM. A hydrology model calculates subglacial drainage paths and indicates two main outlets with different discharges. Depending on the discharge, the plume model computes frontal melt rates, which are iteratively projected to the actual front of the glacier at subglacial discharge locations. This produces undercutting of different sizes, as melt is concentrated close to the surface for high discharge and is more diffuse for low discharge. By testing different configurations, we show that undercutting plays a key role in glacier retreat and is necessary to reproduce observed retreat in the vicinity of the discharge locations during the melting season. Calving rates are also influenced by basal friction, through its effects on near-terminus strain rates and ice velocity

Subglacial hydrology from high-resolution ice-flow simulations of the Rhine Glacier during the Last Glacial Maximum: a proxy for glacial erosion

At the Last Glacial Maximum (LGM), the Rhine Glacier complex (Rhine and Linth glaciers) formed large piedmont lobes extending north into the Swiss and German Alpine forelands. Numerous overdeepened valleys there were formed by repeated glaciations. A characteristic of these overdeepened valleys is their location close to the LGM ice margin, away from the Alps. Numerical models of ice flow of the Rhine Glacier indicate a poor fit between the sliding distance, a proxy for glacial erosion, and the location of these overdeepenings. Calculations of the hydraulic potential based on the computed time-dependent ice surface elevations of the Rhine Glacier lobe obtained from a high-resolution thermo-mechanically coupled Stokes flow model are used to estimate the location of subglacial water drainage routes. Results indicate that the subglacial water discharge is high and focused along glacial valleys and overdeepenings when water pressure is equal to the ice overburden pressure. These conditions are necessary for subglacial water to remove basal sediments, expose fresh bedrock, and favor further erosion by quarrying and abrasion. Knowledge of the location of paleo-subglacial water drainage routes may be useful to understand patterns of subglacial erosion beneath paleo-ice masses that do not otherwise relate to the sliding of ice. Comparison of the erosion pattern from subglacial meltwater with those from quarrying and abrasion shows the importance of subglacial water flow in the formation of distal overdeepenings in the Swiss lowlands.</p

Coupling of ice-shelf melting and buttressing is a key process in ice-sheets dynamics

Increase in ice-shelf melting is generally presumed to have triggered recent coastal ice-sheet thinning. Using a full-Stokes finite element model which includes a proper description of the grounding line dynamics, we investigate the impact of melting below ice shelves. We argue that the influence of ice-shelf melting on the ice-sheet dynamics induces a complex response, and the first naive view that melting inevitably leads to loss of grounded ice is erroneous. We demonstrate that melting acts directly on the magnitude of the buttressing force by modifying both the area experiencing lateral resistance and the ice-shelf velocity, indicating that the decrease of back stress imposed by the ice-shelf is the prevailing cause of inland dynamical thinning. We further show that feedback from melting and buttressing forces can lead to nontrivial results, as an increase in the average melt rate may lead to inland ice thickening and grounding line advance. Citation: Gagliardini, O., G. Durand, T. Zwinger, R. C. A. Hindmarsh, and E. Le Meur (2010), Coupling of ice-shelf melting and buttressing is a key process in ice-sheets dynamics, Geophys. Res. Lett., 37, L14501, doi:10.1029/2010GL043334

Sensitivity of basal conditions in an inverse model : Vestfonna ice cap, Nordaustlandet/Svalbard

Peer reviewe