Investigating changes in basal conditions of Variegated Glacier prior to and during its 1982–1983 surge

Variegated Glacier (Alaska) is known to surge periodically after a sufficient amount of cumulative mass balance is reached, but this observation is difficult to link with changes in the basal conditions. Here, using a 10-yr dataset, consisting of surface topography and surface velocity observations along a flow line for 25 dates, we have reconstructed the evolution of the basal conditions prior to and during the 1982–1983 surge. The model solves the full-Stokes problem along the central flow line using the finite element method. For the 25 dates of the dataset, the basal friction parameter distribution is inferred using the inverse method proposed by Arthern and Gudmundsson (2010). This method is here slightly modified by incorporating a regularisation term in the cost function to avoid short wavelength changes in the friction parameter. Our results indicate that dramatic changes in the basal conditions occurred between 1973 to 1983. Prior to the surge, periodic changes can be observed between winter and summer, with a regular increase of the sliding from 1973 to 1982. During the surge, the basal friction decreased dramatically and an area of very low friction moved from the upper part of the glacier to its terminus. Using a more complex friction law, these changes in basal sliding are then interpreted in terms of basal water pressure. Our results support that dramatic changes took place in the subglacial drainage system of Variegated Glacier, moving from a relatively efficient drainage system prior to the surge to an inefficient one during the surge. By reconstructing the water pressure evolution at the base of the glacier it is possible to propose a scenario for the hydrological history leading to the occurrence of a surge

Sensitivity of centennial mass loss projections of the Amundsen basin to the friction law

Reliable projections of ice sheets' future contributions to sea-level rise require models that are able to accurately simulate grounding-line dynamics, starting from initial states consistent with observations. Here, we simulate the centennial evolution of the Amundsen Sea Embayment in response to a prescribed perturbation in order to assess the sensitivity of mass loss projections to the chosen friction law, depending on the initialisation strategy. To this end, three different model states are constructed by inferring both the initial basal shear stress and viscosity fields with various relative weights. Then, starting from each of these model states, prognostic simulations are carried out using a Weertman, a Schoof and a Budd friction law, with different parameter values. Although the sensitivity of projections to the chosen friction law tends to decrease when more weight is put on viscosity during initialisation, it remains significant for the most physically acceptable of the constructed model states. Independently of the considered model state, the Weertman law systematically predicts the lowest mass losses. In addition, because of its particular dependence on effective pressure, the Budd friction law induces significantly different grounding-line retreat patterns than the other laws and predicts significantly higher mass losses.</p

Impact of mesh resolution for MISMIP and MISMIP3d experiments using Elmer/Ice

International audienceThe 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

Change of the ice rheology with climatic transitions – implication on ice flow modelling and dating of the EPICA Dome C core

International audienceThe study of the distribution of the crystallographic orientations (the fabric) along ice cores supplies information on the past and current ice flows of ice-sheets. Beside the usually observed formation of a vertical single maximum fabric, the EPICA Dome Concordia ice core (EDC) shows an abrupt and unexpected strenghtening of its fabric during termination II around 1750 m depth. Such strengthenings were already observed for sites located on an ice-sheet. This suggests that horizontal shear could occur along the EDC core. Moreover, the change in the fabric leads to a modification of the viscosity between neighbouring ice layers. Through the use of an anisotropic ice flow model, we quantify the change in viscosity and investigate its implication on ice flow and dating

Change in ice rheology during climate variations &ndash; implications for ice flow modelling and dating of the EPICA Dome C core

The study of the distribution of crystallographic orientations (i.e., the fabric) along ice cores provides information on past and current ice flow in ice-sheets. Besides the usually observed formation of a vertical single maximum fabric, the EPICA Dome C ice core (EDC) shows an abrupt and unexpected strengthening of its fabric during termination II around 1750 m depth. Such strengthening has already been observed for sites located on an ice-sheet flank. This suggests that horizontal shear could occur along the EDC core. Moreover, the change in the fabric leads to a modification of the effective viscosity between neighbouring ice layers. Through the use of an anisotropic ice flow model, we quantify the change in effective viscosity and investigate its implication for ice flow and dating

Brief communication: Candidate sites of 1.5 Myr old ice 37 km southwest of the Dome C summit, East Antarctica

The search for ice as old as 1.5 Myr requires the identification of places that maximize our chances to retrieve old, well-resolved, undisturbed and datable ice. One of these locations is very likely southwest of the Dome C summit, where elevated bedrock makes the ice thin enough to limit basal melting. A 3-D ice flow simulation is used to calculate five selection criteria, which together delineate the areas with the most appropriate glaciological properties. These selected areas (a few square kilometers) lie on the flanks of a bedrock high, where a balance is found between risks of basal melting, stratigraphic disturbances and sufficient age resolution. Within these areas, several sites of potential 1.5 Myr old ice are proposed, situated on local bedrock summits or ridges. The trajectories of the ice particles towards these locations are short, and the ice flows over a smoothly undulating bedrock. These sites will help to choose where new high-resolution ground radar surveys should be conducted in upcoming field seasons

Assimilation of Antarctic velocity observations provides evidence for uncharted pinning points

International audienceIn 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

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

Peer reviewe

Application of a two-step approach for mapping ice thickness to various glacier types on Svalbard

The basal topography is largely unknown beneath most glaciers and ice caps, and many attempts have been made to estimate a thickness field from other more accessible information at the surface. Here, we present a two-step reconstruction approach for ice thickness that solves mass conservation over single or several connected drainage basins. The approach is applied to a variety of test geometries with abundant thickness measurements including marine- and land-terminating glaciers as well as a 2400-km2 ice cap on Svalbard. The input requirements are kept to a minimum for the first step. In this step, a geometrically controlled, non-local flux solution is converted into thickness values relying on the shallow ice approximation (SIA). In a second step, the thickness field is updated along fast-flowing glacier trunks on the basis of velocity observations. Both steps account for available thickness measurements. Each thickness field is presented together with an error-estimate map based on a formal propagation of input uncertainties. These error estimates point out that the thickness field is least constrained near ice divides or in other stagnant areas. Withholding a share of the thickness measurements, error estimates tend to overestimate mismatch values in a median sense. We also have to accept an aggregate uncertainty of at least 25-% in the reconstructed thickness field for glaciers with very sparse or no observations. For Vestfonna ice cap (VIC), a previous ice volume estimate based on the same measurement record as used here has to be corrected upward by 22-%. We also find that a 13-% area fraction of the ice cap is in fact grounded below sea level. The former 5-% estimate from a direct measurement interpolation exceeds an aggregate maximum range of 6-23-% as inferred from the error estimates here.This study received primary funding from the German Research Foundation (DFG) under grant number FU1032/1-1. Results presented in this publication are based on numerical simulations conducted at the high-performance computing centre of the Regionales Rechenzentrum Erlangen (RRZE) of the University of Erlangen-Nuremberg. The reconstruction approach also benefits from co-development work of the Elmer/Ice team at the CSC-IT Center for Science Ltd. (Finland). The velocity analysis on Svalbard was funded by DFG within the priority programme 1158 Antarctic Research with Comparable Investigations in Arctic Sea Ice Areas under contract number BR2105/9-1 and received financial support from the Helmholtz Association of the German Research Centres (HGF) Alliance on Remote Sensing and Earth System Dynamics. Thickness data collection in Wedel Jarlsberg Land was funded by the Spanish R&D projects C11093001 and C150954001, NCBiR/PolarCLIMATE-2009/2-2/2010 from the Polish National Centre for R&D, by IPY/269/2006 from the Polish Ministry of Science and Higher Education, by Polish-Norwegian funding through the AWAKE (PNRF-22-AI-1/07) project, by the EU FP7 ice2sea programme (grant number 226375) and by funds of the Leading National Research Centre (KNOW) received by the Centre for Polar Studies of the University of Silesia, Poland. The DEM generation inWedel Jarlsberg Land received financial support from the European Research Council (grant 320816) and from ESA (project Glaciers CCI, 4000109873/14/I-NB). TanDEM-X data were provided under AO XTIGLAC6770. The WRF-SMB field was produced within the PERMANOR project funded by the Norwegian Research Council (255331)

Ice fabric in an Antarctic ice stream interpreted from seismic anisotropy

Here we present new measurements of an anisotropic ice fabric in a fast moving (377 ma−1) ice stream in West Antarctica. We use ∼6000 measurements of shear wave splitting observed in microseismic signals from the bed of Rutford Ice Stream, to show that in contrast to large-scale ice flow models, which assume that ice is isotropic, the ice in Rutford Ice Stream is dominated by a previously unobserved type of partial girdle fabric. This fabric has a strong directional contrast in mechanical properties, shearing 9.1 times more easily along the ice flow direction than across flow. This observed fabric is likely to be widespread and representative of fabrics in other ice streams and large glaciers, suggesting it is essential to consider anisotropy in data-driven models to correctly predict ice loss and future flow in these regions. We show how passive microseismic monitoring can be effectively used to provide these data