ALBERT

Description: The majority of Antarctic ice is discharged via long and narrow fast-flowing ice streams. At ice stream margins, the rapid transition from the vertical shearing flow in the ice ridges surrounding the stream to a rapidly sliding plug flow in the stream itself leads to high stress concentrations and a velocity field whose form is non-trivial to determine. In this paper, we develop a boundary layer theory for this narrow region separating a lubrication-type ice ridge flow and a membrane-type ice stream flow. This allows us to derive jump conditions for the outer models describing ridge and stream self-consistently. Much of our focus is, however, on determining the velocity and shear heating fields in the margin itself. Ice stream margins have been observed to change position over time, with potentially significant implications for ice stream discharge. Our boundary layer model allows us to extend previous work that has determined rates of margin migration from a balance between shear heating in the margin and the cooling effect of margin migration into the colder ice of the surrounding ice ridge. Solving for the transverse velocity field in the margin allows us to include the effect of advection due to lateral inflow of ice from the ridge on margin migration, and we demonstrate that this reduces the rate of margin migration, as previously speculated. © 2015 Cambridge University Press.

Description: The dynamical contribution of marine ice sheets to sea level rise is largely controlled by grounding line (GL) dynamics. Seroussi et al. (2014) emphasised the sensitivity of numerical ice flow model results to the practical implementation of the friction of the ice on its bed in the very close vicinity of the GL. Elmer/Ice is a reference finite element (FE) ice flow model used in recent marine ice sheet model intercomparison (MISMIP) exercises. In the model, the GL is defined as the nodes where the ice is in contact with the bedrock but belong to both grounded and floating elements. Inherently to the FE method, computing the contribution of the friction by element requires evaluating the friction at the integration points. In Elmer/Ice, this is done by interpolating the values of the friction parameter C prescribed at the nodes. In this brief communication, we discuss and compare three alternative ways to prescribe the friction at the GL: (i) C is prescribed and non null at the GL nodes, (ii) C is set to zero at the GL nodes, and (iii) C is discontinuous at the GL nodes (i.e. is prescribed and non null for grounded elements and otherwise null). So far, all published results using Elmer/Ice were obtained with the first method. Using the MISMIP3d diagnostic experiment, we first show that, although the change in the total force at the base is insignificant, the three methods lead to significantly different velocity fields. We then show that these methods also lead to different steady state GL positions and different transient behaviours. Such model sensitivity to the methods discussed here is certainly specific to the high friction prescribed in the MISMIP experiments and should be smaller in real setups where friction in the vicinity of the GL would be expected to be lower. Results obtained with the three methods are available as Supplement for future comparisons.

Description: Two mechanisms are generally proposed to explain seasonal variations in the calving front of tidewater glaciers: submarine melting of the calving face and the mechanical back-force applied by the ice mélange. However, the way these processes affect the calving rate and the glacier dynamics remains uncertain. In this study, we used the finite element model Elmer/Ice to simulate the impact of these forcings on more than 200 two dimensional theoretical flowline glacier configurations. The model, which includes calving processes, suggests that frontal melting affects the position of the terminus only slightly (〈 a few hundred meters) and does not affect the pluriannual glacier mass balance at all. However, the ice mélange has a greater impact on the advance and retreat cycles of the glacier front (more than several 1000 m) and its consequences for the mass balance are not completely negligible, stressing the need for better characterization of forcing properties. We also show that ice mélange forcing against the calving face can mechanically prevent crevasse propagation at sea level and hence prevent calving. Results also revealed different behaviors in grounded and floating glaciers: in the case of a floating extension, the heaviest forcings can disrupt the glacier equilibrium by modifying its buttressing and ice flux at the grounding line.

Description: Radar inference of the bulk properties of glacier beds, most notably identifying basal melting, is, in general, derived from the basal reflection coefficient. On the scale of an ice sheet, unambiguous determination of basal reflection is primarily limited by uncertainty in the englacial attenuation of the radio wave, which is an Arrhenius function of temperature. Existing bed-returned power algorithms for deriving attenuation assume that the attenuation rate is regionally constant, which is not feasible at an ice-sheet-wide scale. Here we introduce a new semi-empirical framework for deriving englacial attenuation, and, to demonstrate its efficacy, we apply it to the Greenland Ice Sheet. A central feature is the use of a prior Arrhenius temperature model to estimate the spatial variation in englacial attenuation as a first guess input for the radar algorithm. We demonstrate regions of solution convergence for two input temperature fields and for independently analysed field campaigns. The coverage achieved is a trade-off with uncertainty and we propose that the algorithm can be "tuned" for discrimination of basal melt (attenuation loss uncertainty  ∼ 5 dB). This is supported by our physically realistic ( ∼ 20 dB) range for the basal reflection coefficient. Finally, we show that the attenuation solution can be used to predict the temperature bias of thermomechanical ice sheet models and is in agreement with known model temperature biases at the Dye 3 ice core.

Description: In ice flow modelling, the use of control methods to assimilate the dynamic and geometric state of an ice body has become common practice. These methods have primarily focussed on inverting for one of the two least known properties in glaciology, namely the basal friction coefficient or the ice viscosity parameter. Here, we present an approach to infer both properties simultaneously for the whole of the Antarctic ice sheet. During the assimilation, the root-mean-square deviation between modelled and observed surface velocities is reduced to 12.3 m a−1, with a value of 16.4 m a−1 for the ice shelves. An exception in terms of the velocity mismatch is the Thwaites Glacier ice shelf, where the RMS value attains almost 80 m a−1. The reason is that the underlying BEDMAP2 geometry ignores the presence of an ice rise, that exerts major control on the dynamics of the eastern part of the ice shelf. On these grounds, we suggest an approach to account for pinning points not included in BEDMAP2 by locally allowing an optimisation of basal friction during the inversion. In this way, the velocity mismatch on the Thwaites ice shelf is more than halved. A characteristic velocity mismatch pattern emerges for unaccounted pinning points close to the marine shelf front. This pattern is exploited to manually identify 7 uncharted features around Antarctica that exert significant resistance to the shelf flow. Potential pinning points are detected on Fimbul, West, Shakelton, Nickerson and Venable ice shelves. As pinning points can provide substantial resistance to shelf flow, with considerable consequences if they became ungrounded in the future, the model community is in need of detailed bathymetry there. Our data assimilation points to some of these dynamically important features, not present in BEDMAP2, and implicitly quantifies their relevance.

Description: Submarine melting of the calving face of tidewater glaciers and the mechanical back force applied by the ice mélange layer are two mechanisms generally proposed to explain seasonal variations at the calving front of tidewater glaciers. However, the way these processes affect the calving rate and glacier dynamics remains uncertain. In this study, we used a finite element-based model that solves the full Stokes equations to simulate the impact of these forcings on two-dimensional theoretical flow line glacier configurations. The model, which includes calving processes, suggests that frontal melting affects the position of the terminus only slightly (less than a few hundred metres) and does not affect the multiannual glacier mass balance at all. However, the ice mélange has a greater impact on the advance and retreat cycles of the glacier front (more than several kilometres) and its consequences for the mass balance are not completely negligible, stressing the need for better characterization of forcing properties. We also show that ice mélange forcing against the calving face can mechanically prevent crevasse propagation at sea level and hence prevent calving. Results also reveal different behaviours in grounded and floating glaciers: in the case of a floating extension, the strongest forcings can disrupt the glacier equilibrium by modifying its buttressing and ice flux at the grounding line.

Description: In ice flow modelling, the use of control methods to assimilate the dynamic and geometric state of an ice body has become common practice. These methods have primarily focussed on inverting for one of the two least known properties in glaciology, namely the basal friction coefficient or the ice viscosity parameter. Here, we present an approach to infer both properties simultaneously for the whole of the Antarctic ice sheet. After the assimilation, the root-mean-square deviation between modelled and observed surface velocities attains 8.7 m a−1 for the entire domain, with a slightly higher value of 14.0 m a−1 for the ice shelves. An exception in terms of the velocity mismatch is the Thwaites Glacier Ice Shelf, where the RMS value is almost 70 m a−1. The reason is that the underlying Bedmap2 geometry ignores the presence of an ice rise, which exerts major control on the dynamics of the eastern part of the ice shelf. On these grounds, we suggest an approach to account for pinning points not included in Bedmap2 by locally allowing an optimisation of basal friction during the inversion. In this way, the velocity mismatch on the ice shelf of Thwaites Glacier is more than halved. A characteristic velocity mismatch pattern emerges for unaccounted pinning points close to the marine shelf front. This pattern is exploited to manually identify seven uncharted features around Antarctica that exert significant resistance to the shelf flow. Potential pinning points are detected on Fimbul, West, Shackleton, Nickerson and Venable ice shelves. As pinning points can provide substantial resistance to shelf flow, with considerable consequences if they became ungrounded in the future, the model community is in need of detailed bathymetry there. Our data assimilation points to some of these dynamically important features not present in Bedmap2 and implicitly quantifies their relevance.

Description: The dynamical contribution of marine ice sheets to sea level rise is largely controlled by grounding line (GL) dynamics. Two marine ice sheet model intercomparison exercises, namely MISMIP and MISMIP3d, have been proposed to the community to test and compare the ability of models to capture the GL dynamics. Both exercises are known to present a discontinuity of the friction at the GL, which is believed to increase the model sensitivity to mesh resolution. Here, using Elmer/Ice, the only Stokes model which completed both intercomparisons, the sensitivity to the mesh resolution is studied from an extended MISMIP experiment in which the friction continuously decreases over a transition distance and equals zero at the GL. Using this MISMIP-like setup, it is shown that the sensitivity to the mesh resolution is not improved for a vanishing friction at the GL. For the original MISMIP experiment, i.e. for a discontinuous friction at the GL, we further show that the results are moreover very sensitive to the way the friction is interpolated in the close vicinity of the GL. In the light of these new insights, and thanks to increased computing resources, new results for the MISMIP3d experiments obtained for higher resolutions than previously published are made available for future comparisons as the Supplement.

Description: The deadliest outburst flood from an englacial cavity occurred on Glacier de Tête Rousse in the Mont Blanc area, French Alps, in 1892. A subglacial reservoir was discovered in the same glacier in 2010 and drained artificially in 2010, 2011 and 2012 to protect the 3000 inhabitants downstream. The mechanism leading to the spontaneous refilling of the cavity following these pumping operations has been analyzed. For this purpose, the subglacial water volume changes between 2010 and 2013 were reconstructed. The size of the cavity following the pumping was found to have decreased from 53500 m3 in 2010 to 12 750 m3 in 2013. Creep and the partial collapse of the cavity roof explain a large part of the volume loss. Analysis of cavity filling showed a strong relationship between measured surface melting and the filling rate, with a time delay of 4–6 hours. A permanent input of 15 m3 d−1, not depending on surface melt, was also found. The meltwater and rain from the surface is conveyed to bedrock through crevasses and probably through a permeable layer of rock debris at the glacier bed. The drainage pathway permeability was estimated at 0.054 ms−1 from water discharge measurements and dye-tracing experiments.