Publications

Abstract: Coastal areas are highly diverse, ecologically rich, regions of key socio-economic activity, and are particularly sensitive to sea-level change. Over most of the 20th century, global mean sea level has risen mainly due to warming and subsequent expansion of the upper ocean layers as well as the melting of glaciers and ice caps. Over the last three decades, increased mass loss of the Greenland and Antarctic ice sheets has also started to contribute significantly to contemporary sea-level rise. The future mass loss of the two ice sheets, which combined represent a sea-level rise potential of similar to 65 m, constitutes the main source of uncertainty in long-term (centennial to millennial) sea-level rise projections. Improved knowledge of the magnitude and rate of future sea-level change is therefore of utmost importance. Moreover, sea level does not change uniformly across the globe and can differ greatly at both regional and local scales. The most appropriate and feasible sea level mitigation and adaptation measures in coastal regions strongly depend on local land use and associated risk aversion. Here, we advocate that addressing the problem of future sea-level rise and its impacts requires (i) bringing together a transdisciplinary scientific community, from climate and cryospheric scientists to coastal impact specialists, and (ii) interacting closely and iteratively with users and local stakeholders to co-design and co-build coastal climate services, including addressing the high-end risks.

Abstract: We present the result of the third Marine Ice Sheet Model Intercomparison Project, MISMIP+. MISMIP+ is intended to be a benchmark for ice-flow models which include fast sliding marine ice streams and floating ice shelves and in particular a treatment of viscous stress that is sufficient to model buttressing, where upstream ice flow is restrained by a downstream ice shelf. A set of idealized experiments first tests that models are able to maintain a steady state with the grounding line located on a retrograde slope due to buttressing and then explore scenarios where a reduction in that buttressing causes ice stream acceleration, thinning, and grounding line retreat. The majority of participating models passed the first test and then produced similar responses to the loss of buttressing. We find that the most important distinction between models in this particular type of simulation is in the treatment of sliding at the bed, with other distinctions – notably the difference between the simpler and more complete treatments of englacial stress but also the differences between numerical methods – taking a secondary role.

Abstract: Predictions for sea-level rise this century due to melt from Antarctica range from zero to more than one metre. The highest predictions are driven by the controversial marine ice-cliff instability (MICI) hypothesis, which assumes that coastal ice cliffs can rapidly collapse after ice shelves disintegrate, as a result of surface and sub-shelf melting caused by global warming. But MICI has not been observed in the modern era and it remains unclear whether it is required to reproduce sea-level variations in the geological past. Here we quantify ice-sheet modelling uncertainties for the original MICI study and show that the probability distributions are skewed towards lower values (under very high greenhouse gas concentrations, the most likely value is 45 centimetres). However, MICI is not required to reproduce sea-level changes due to Antarctic ice loss in the mid-Pliocene epoch, the last interglacial period or 1992-2017; without it we find that the projections agree with previous studies (all 95th percentiles are less than 43 centimetres). We conclude that previous interpretations of these MICI projections over-estimate sea-level rise this century; because the MICI hypothesis is not well constrained, confidence in projections with MICI would require a greater range of observationally constrained models of ice-shelf vulnerability and ice-cliff collapse.

Abstract: The sensitivity of Antarctic sea-ice to increasing glacial freshwater release into the Southern Ocean is studied in a series of 31-year ocean/sea-ice/iceberg model simulations. Glaciological estimates of ice-shelf melting and iceberg calving are used to better constrain the spatial distribution and magnitude of freshwater forcing around Antarctica. Two scenarios of glacial freshwater forcing have been designed to account for a decadal perturbation in glacial freshwater release to the Southern Ocean. For the first time, this perturbation explicitly takes into consideration the spatial distribution of changes in the volume of Antarctic ice shelves, which is found to be a key component of changes in freshwater release. In addition, glacial freshwater-induced changes in sea ice are compared to typical changes induced by the decadal evolution of atmospheric states. Our results show that, in general, the increase in glacial freshwater release increases Antarctic sea ice extent. But the response is opposite in some regions like the coastal Amundsen Sea, implying that distinct physical mechanisms are involved in the response. We also show that changes in freshwater forcing may induce large changes in sea-ice thickness, explaining about one half of the total change due to the combination of atmospheric and freshwater changes. The regional contrasts in our results suggest a need for improving the representation of freshwater sources and their evolution in climate models.

Abstract: Basal slip accounts for a large part of the flow of ice streams draining ice from Antarctica and Greenland into the ocean. Therefore, an appropriate representation of basal slip in ice flow models is a prerequisite for accurate sea level rise projections. Various friction laws have been proposed to describe basal slip in models. Here, we compare the influence on grounding line (GL) dynamics of four friction laws: the traditional Weertman law and three effective pressure-dependent laws, namely the Schoof, Tsai and Budd laws. It turns out that, even when they are tuned to a common initial reference state, the Weertman, Budd and Schoof laws lead to thoroughly different steady-state positions, although the Schoof and Tsai laws lead to much the same result. In particular, under certain circumstances, it is possible to obtain a steady GL located on a reverse slope area using the Weertman law. Furthermore, the predicted transient evolution of the GL as well as the projected contributions to sea level rise over a 100-year time horizon vary significantly depending on the friction law. We conclude on the importance of choosing an appropriate law for reliable sea level rise projections and emphasise the need for a coupling between ice flow models and physically based subglacial hydrological models.

Abstract: Numerical modeling of the Antarctic ice sheet has gone through a paradigm shift over the last decade. While initially models focussed on long-time diffusive response to surface mass balance changes, processes occurring at the marine boundary of the ice sheet are progressively incorporated in newly developed state-of-the-art ice-sheet models. These models now exhibit fast, short-term volume changes, in line with current observations of mass loss. Coupling with ocean models is currently on its way and applied to key areas of the Antarctic ice sheet. New model intercomparisons have been launched, focusing on ice/ocean interaction (MISMIP+, MISOMIP) or ice-sheet model initialization and multi-ensemble projections (ISMIP6). Nevertheless, the inclusion of new processes pertaining to ice-shelf calving, evolution of basal friction, and other processes, also increase uncertainties in the contribution of the Antarctic ice sheet to future sea-level rise.

Abstract: Coupled ice sheet-ocean models capable of simulating moving grounding lines are just becoming available. Such models have a broad range of potential applications in studying the dynamics of marine ice sheets and tidewater glaciers, from process studies to future projections of ice mass loss and sea level rise. The Marine Ice Sheet-Ocean Model Intercomparison Project ( MISOMIP) is a community effort aimed at designing and coordinating a series of model intercomparison projects ( MIPs) for model evaluation in idealized setups, model verification based on observations, and future projections for key regions of the West Antarctic Ice Sheet ( WAIS). Here we describe computational experiments constituting three interrelated MIPs for marine ice sheet models and regional ocean circulation models incorporating ice shelf cavities. These consist of ice sheet experiments under the Marine Ice Sheet MIP third phase ( MISMIP+), ocean experiments under the Ice Shelf-Ocean MIP second phase ( ISOMIP+) and coupled ice sheet-ocean experiments under the MISOMIP first phase ( MISOMIP1). All three MIPs use a shared domain with idealized bedrock topography and forcing, allowing the coupled simulations ( MISOMIP1) to be compared directly to the individual component simulations ( MISMIP+ and ISOMIP+). The experiments, which have qualitative similarities to Pine Island Glacier Ice Shelf and the adjacent region of the Amundsen Sea, are designed to explore the effects of changes in ocean conditions, specifically the temperature at depth, on basal melting and ice dynamics. In future work, differences between model results will form the basis for the evaluation of the participating models.

Abstract: 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.

Abstract: Recent increase in Antarctic freshwater release to the Southern Ocean is suggested to contribute to change in water masses and sea ice. However, climate models differ in their representation of the freshwater sources. Recent improvements in altimetry-based detection of small icebergs and in estimates of the mass loss of Antarctica may help better constrain the values of Antarctic freshwater releases. We propose a model-based seasonal climatology of iceberg melt over the Southern Ocean using state-of-the-art observed glaciological estimates of the Antarctic mass loss. An improved version of a Lagrangian iceberg model is coupled with a global, eddy-permitting ocean/sea ice model and compared to small icebergs observations. Iceberg melt increases sea ice cover, about 10% in annual mean sea ice volume, and decreases sea surface temperature over most of the Southern Ocean, but with distinctive regional patterns. Our results underline the importance of improving the representation of Antarctic freshwater sources. This can be achieved by forcing ocean/sea ice models with a climatological iceberg fresh-water flux. (C) 2016 Elsevier Ltd. All rights reserved.

Abstract: 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 ma(-1) for the entire domain, with a slightly higher value of 14.0 ma(-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 ma(-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.

Abstract: Locally grounded features in ice shelves, called ice rises and rumples, play a key role buttressing discharge from the Antarctic Ice Sheet and regulating its contribution to sea level. Ice rises typically rise several hundreds of meters above the surrounding ice shelf; shelf flow is diverted around them. On the other hand, shelf ice flows across ice rumples, which typically rise only a few tens of meters above the ice shelf. Ice rises contain rich histories of deglaciation and climate that extend back over timescales ranging from a few millennia to beyond the last glacial maximum. Numerical model results have shown that the buttressing effects of ice rises and rumples are significant, but details of processes and how they evolve remain poorly understood. Fundamental information about the conditions and processes that cause transitions between floating ice shelves, ice rises and ice rumples is needed in order to assess their impact on ice-sheet behavior. Targeted high-resolution observational data are needed to evaluate and improve prognostic numerical models and parameterizations of the effects of small-scale pinning points on grounding-zone dynamics. (C) 2015 The Authors. Published by Elsevier B.V. This is an open access article under the CC BY license

Abstract: An important share of paleoclimatic information is buried within the lowermost layers of deep ice cores. Because improving our records further back in time is one of the main challenges in the near future, it is essential to judge how deep these records remain unaltered, since the proximity of the bedrock is likely to interfere both with the recorded temporal sequence and the ice properties. In this paper, we present a multiparametric study (delta D-delta O-18(ice), delta O-18(atm), total air content, CO2, CH4, N2O, dust, high-resolution chemistry, ice texture) of the bottom 60 m of the EPICA (European Project for Ice Coring in Antarctica) Dome C ice core from central Antarctica. These bottom layers were subdivided into two distinct facies: the lower 12 m showing visible solid inclusions (basal dispersed ice facies) and the upper 48 m, which we will refer to as the “basal clean ice facies”. Some of the data are consistent with a pristine paleoclimatic signal, others show clear anomalies It is demonstrated that neither large-scale bottom refreezing of subglacial water, nor mixing (be it internal or with a local basal end term from a previous/initial ice sheet configuration) can explain the observed bottom-ice properties. We focus on the high-resolution chemical profiles and on the available remote sensing data on the subglacial topography of the site to propose a mechanism by which relative stretching of the bottom-ice sheet layers is made possible, due to the progressively confining effect of subglacial valley sides. This stress field change, combined with bottom-ice temperature close to the pressure melting point, induces accelerated migration recrystallization, which results in spatial chemical sorting of the impurities, depending on their state (dissolved vs. solid) and if they are involved or not in salt formation. This chemical sorting effect is responsible for the progressive build-up of the visible solid aggregates that therefore mainly originate “from within”, and not from incorporation processes of debris from the ice sheet's substrate. We further discuss how the proposed mechanism is compatible with the other ice properties described. We conclude that the paleoclimatic signal is only marginally affected in terms of global ice properties at the bottom of EPICA Dome C, but that the timescale was considerably distorted by mechanical stretching of MIS20 due to the increasing influence of the subglacial topography, a process that might have started well above the bottom ice. A clear paleoclimatic signal can therefore not be inferred from the deeper part of the EPICA Dome C ice core. Our work suggests that the existence of a flat monotonic ice bedrock interface, extending for several times the ice thickness, would be a crucial factor in choosing a future “oldest ice” drilling location in Antarctica.

Abstract: Over the past 40 years Pine Island Glacier in West Antarctica has thinned at an accelerating rate(1-3), so that at present it is the largest single contributor to sea-level rise in Antarctica(4). In recent years, the grounding line, which separates the grounded ice sheet from the floating ice shelf, has retreated by tens of kilometres(5). At present, the grounding line is crossing a retrograde bedrock slope that lies well below sea level, raising the possibility that the glacier is susceptible to the marine ice-sheet instability mechanism(6-8). Here, using three state-of-the-art ice-flow models(9-11), we show that Pine Island Glacier's grounding line is probably engaged in an unstable 40 km retreat. The associated mass loss increases substantially over the course of our simulations from the average value of 20 Gt yr(-1) observed for the 1992-2011 period(4), up to and above 100 Gt yr(-1), equivalent to 3.5-10mm eustatic sea-level rise over the following 20 years. Mass loss remains elevated from then on, ranging from 60 to 120 Gt yr(-1).

Abstract: Calving of icebergs is a major negative component of polar ice-sheet mass balance. Here we present a new calving model relying on both continuum damage mechanics and linear elastic fracture mechanics. This combination accounts for both the slow sub-critical surface crevassing and the rapid propagation of crevasses when calving occurs. First, damage to the ice occurs over long timescales and enhances the viscous flow of ice. Then brittle fractures propagate downward, at very short timescales, when the ice body is considered as an elastic medium. The model was calibrated on Helheim Glacier, Southeast Greenland, a well-monitored glacier with fast-flowing outlet. This made it possible to identify sets of model parameters to enable a consistent response of the model and to produce a dynamic equilibrium in agreement with the observed stable position of the Helheim ice front between 1930 and today.

Abstract: Marine ice-sheet stability is mostly controlled by the dynamics of the grounding line, i.e. the junction between the grounded ice sheet and the floating ice shelf. Grounding line migration has been investigated within the framework of MISMIP (Marine Ice Sheet Model Intercomparison Project), which mainly aimed at investigating steady state solutions. Here we focus on transient behaviour, executing short-term simulations (200 yr) of a steady ice sheet perturbed by the release of the buttressing restraint exerted by the ice shelf on the grounded ice upstream. The transient grounding line behaviour of four different flowline ice-sheet models has been compared. The models differ in the physics implemented (full Stokes and shallow shelf approximation), the numerical approach, as well as the grounding line treatment. Their overall response to the loss of buttressing is found to be broadly consistent in terms of grounding line position, rate of surface elevation change and surface velocity. However, still small differences appear for these latter variables, and they can lead to large discrepancies (> 100 %) observed in terms of ice sheet contribution to sea level when cumulated over time. Despite the recent important improvements of marine ice-sheet models in their ability to compute steady state configurations, our results question the capacity of these models to compute short-term reliable sea-level rise projections.

Abstract: The Fourth IPCC Assessment Report concluded that ice sheet flow models, in their current state, were unable to provide accurate forecast for the increase of polar ice sheet discharge and the associated contribution to sea level rise. Since then, the glaciological community has undertaken a huge effort to develop and improve a new generation of ice flow models, and as a result a significant number of new ice sheet models have emerged. Among them is the parallel finite-element model Elmer/Ice, based on the open-source multi-physics code Elmer. It was one of the first full-Stokes models used to make projections for the evolution of the whole Greenland ice sheet for the coming two centuries. Originally developed to solve local ice flow problems of high mechanical and physical complexity, Elmer/Ice has today reached the maturity to solve larger-scale problems, earning the status of an ice sheet model. Here, we summarise almost 10 yr of development performed by different groups. Elmer/Ice solves the full-Stokes equations, for isotropic but also anisotropic ice rheology, resolves the grounding line dynamics as a contact problem, and contains various basal friction laws. Derived fields, like the age of the ice, the strain rate or stress, can also be computed. Elmer/Ice includes two recently proposed inverse methods to infer badly known parameters. Elmer is a highly parallelised code thanks to recent developments and the implementation of a block preconditioned solver for the Stokes system. In this paper, all these components are presented in detail, as well as the numerical performance of the Stokes solver and developments planned for the future.

Abstract: Predictions of marine ice-sheet behaviour require models able to simulate grounding-line migration. We present results of an intercomparison experiment for plan-view marine ice-sheet models. Verification is effected by comparison with approximate analytical solutions for flux across the grounding line using simplified geometrical configurations (no lateral variations, no buttressing effects from lateral drag). Perturbation experiments specifying spatial variation in basal sliding parameters permitted the evolution of curved grounding lines, generating buttressing effects. The experiments showed regions of compression and extensional flow across the grounding line, thereby invalidating the boundary layer theory. Steady-state grounding-line positions were found to be dependent on the level of physical model approximation. Resolving grounding lines requires inclusion of membrane stresses, a sufficiently small grid size (<500 m), or subgrid interpolation of the grounding line. The latter still requires nominal grid sizes of <5 km. For larger grid spacings, appropriate parameterizations for ice flux may be imposed at the grounding line, but the short-time transient behaviour is then incorrect and different from models that do not incorporate grounding-line parameterizations. The numerical error associated with predicting grounding-line motion can be reduced significantly below the errors associated with parameter ignorance and uncertainties in future scenarios.

Abstract: Over the last two decades, the Greenland ice sheet (GrIS) has been losing mass at an increasing rate, enhancing its contribution to sea-level rise (SLR). The recent increases in ice loss appear to be due to changes in both the surface mass balance of the ice sheet and ice discharge (ice flux to the ocean). Rapid ice flow directly affects the discharge, but also alters ice-sheet geometry and so affects climate and surface mass balance. Present-day ice-sheet models only represent rapid ice flow in an approximate fashion and, as a consequence, have never explicitly addressed the role of ice discharge on the total GrIS mass balance, especially at the scale of individual outlet glaciers. Here, we present a new-generation prognostic ice-sheet model which reproduces the current patterns of rapid ice flow. This requires three essential developments: the complete solution of the full system of equations governing ice deformation; a variable resolution unstructured mesh to resolve outlet glaciers and the use of inverse methods to better constrain poorly known parameters using observations. The modelled ice discharge is in good agreement with observations on the continental scale and for individual outlets. From this initial state, we investigate possible bounds for the next century ice-sheet mass loss. We run sensitivity experiments of the GrIS dynamical response to perturbations in climate and basal lubrication, assuming a fixed position of the marine termini. We find that increasing ablation tends to reduce outflow and thus decreases the ice-sheet imbalance. In our experiments, the GrIS initial mass (im)balance is preserved throughout the whole century in the absence of reinforced forcing, allowing us to estimate a lower bound of 75 mm for the GrIS contribution to SLR by 2100. In one experiment, we show that the current increase in the rate of ice loss can be reproduced and maintained throughout the whole century. However, this requires a very unlikely perturbation of basal lubrication. From this result we are able to estimate an upper bound of 140 mm from dynamics only for the GrIS contribution to SLR by 2100.

Abstract: Recent glaciological surveys have revealed a significant increase of ice discharge from polar ice caps into the ocean. In parallel, ice flow models have been greatly improved to better reproduce current changes and forecast the future behavior of ice sheets. For these models, surface topography and bedrock elevation are crucial input parameters that largely control the dynamics and the ensuing overall mass balance of the ice sheet. For obvious reasons of inaccessibility, only sparse and uneven bedrock elevation data is available. This raw data is processed to produce Digital Elevation Models (DEMs) on a regular 5 km grid. These DEMs are used to constrain the basal boundary conditions of all ice sheet models. Here, by using a full-Stokes finite element code, we examine the sensitivity of an ice flow model to the accuracy of the bedrock description. In the context of short-term ice sheet forecast, we show that in coastal regions, the bedrock elevation should be known at a resolution of the order of one kilometer. Conversely, a crude description of the bedrock in the interior of the continent does not affect modeling of the ice outflow into the ocean. These findings clearly indicate that coastal regions should be prioritized during future geophysical surveys. They also indicate that a paradigm shift is required to change the current design of DEMs describing the bedrock below the ice sheets: they must give users the opportunity to incorporate high-resolution bedrock elevation data in regions of interest. Citation: Durand, G., O. Gagliardini, L. Favier, T. Zwinger, and E. le Meur (2011), Impact of bedrock description on modeling ice sheet dynamics, Geophys. Res. Lett., 38, L20501, doi: 10.1029/2011GL048892.

Abstract: 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.

Abstract: Polar ice is known to be one of the most anisotropic natural materials. For a given fabric the polycrystal viscous response is strongly dependent on the actual state of stress and strain rate. Within an ice sheet, grounded-ice parts and ice shelves have completely different stress regimes, so one should expect completely different impacts of ice anisotropy on the flow. The aim of this work is to quantify, through the concept of enhancement factors, the influence of ice anisotropy on the flow of grounded ice and ice shelves. For this purpose, a full-Stokes anisotropic marine ice-sheet flowline model is used to compare isotropic and anisotropic diagnostic velocity fields on a fixed geometry. From these full-Stokes results, we propose a definition of enhancement factors for grounded ice and ice shelves, coherent with the asymptotic models used for these regions. We then estimate realistic values for the enhancement factors induced by ice anisotropy for grounded ice and ice shelves.

Abstract: [1] The stability of marine ice sheets and outlet glaciers is mostly controlled by the dynamics of their grounding line, i.e., where the bottom contact of the ice changes from bedrock or till to ocean water. The last report of the Intergovernmental Panel on Climate Change has clearly underlined the poor ability of models to capture the dynamics of outlet glaciers. Here we present computations of grounding line dynamics on the basis of numerical solutions of the full Stokes equations for ice velocity, coupled with the evolution of the air ice- and sea ice-free interfaces. The grounding line position is determined by solving the contact problem between the ice and a rigid bedrock using the finite element code Elmer. Results of the simulations show that marine ice sheets are unstable on upsloping beds and undergo hysteresis under perturbation of ice viscosity, confirming conclusions from boundary layer theory. The present approach also indicates that a 2-D unconfined marine ice sheet sliding over a downsloping bedrock does not exhibit neutral equilibrium. It is shown that mesh resolution around the grounding line is a crucial issue. A very fine grid size (< 100 m spacing) is needed in order to achieve consistent results.

Abstract: While there are no indications of mixing back to 800 000 years in the EPICA Dome C ice core record, comparison with marine sediment records shows significant differences in the timing and duration of events prior to stage 11 (similar to 430 ka, thousands of years before 1950). A relationship between the isotopic composition of atmospheric oxygen (delta O-18 of O-2, noted delta O-18(atm))and daily northern hemisphere summer insolation has been observed for the youngest four climate cycles. Here we use this relationship with new delta O-18 of O-2 measurements to show that anomalous flow in the bottom 500 m of the core distorts the duration of events by up to a factor of 2. By tuning delta O-18(atm) to orbital precession we derive a corrected thinning function and present a revised age scale for the interval corresponding to Marine Isotope Stages 11-20 in the EPICA Dome C ice core. Uncertainty in the phasing of delta O-18(atm) with respect to insolation variations in the precession band limits the accuracy of this new agescale to +/- 6 kyr (thousand of years). The previously reported similar to 30 kyr duration of interglacial stage 11 is unchanged. In contrast, the duration of stage 15.1 is reduced by a factor of 2, from 31 to 16 kyr.

Abstract: One-dimensional ( 1-D) ice flow models are used to construct the age scales at the Dome C and Dome Fuji drilling sites ( East Antarctica). The poorly constrained glaciological parameters at each site are recovered by fitting independent age markers identified within each core. We reconstruct past accumulation rates, that are larger than those modelled using the classical vapour saturation pressure relationship during glacial periods by up to a factor 1.5. During the Early Holocene, changes in reconstructed accumulation are not linearly related to changes in ice isotopic composition. A simple model of past elevation changes is developed and shows an amplitude variation of 110-120m at both sites. We suggest that there is basal melting at Dome C ( 0.56 +/- 0.19 mm/yr). The reconstructed velocity profile is highly non-linear at both sites, which suggests complex ice flow effects. This induces a non-linear thinning function in both drilling sites, which is also characterized by bumps corresponding to variations in ice thickness with time.

Abstract: [1] On the basis of a detailed study of the ice microstructure of the European Project for Ice Coring in Antarctica (EPICA) ice core at Dome Concordia, Antarctica, we analyze the effect of impurities (solubles, and insolubles, that is, dust particles) on the grain growth process in cold ice sheets. As a general trend, the average grain size increases with depth. This global increase, induced by the normal grain growth process, is punctuated by several sharp decreases that can be associated with glacial-interglacial climatic transitions. To explain the modifications of the microstructure with climatic changes, we discuss the role of soluble and insoluble impurities on the grain growth process, coupled with an analysis of the pinning of grain boundaries by microparticles. Our data indicate that high soluble impurity content does not necessarily imply a slowdown of grain growth kinetics, whereas the pinning of grain boundaries by dust explains all the observed modifications of the microstructure. We propose a numerical model of the evolution of the average grain size in deep ice cores that takes into account recrystallization processes such as normal grain growth and rotation recrystallization as well as the pinning effect induced by dust particles, bubbles, and clathrates on the grain boundaries. Applied to the first 2135 m of the Dome Concordia core, the model reproduces accurately the measured mean grain radius. This indicates a major role of dust in the modification of polar ice microstructure and shows that the average grain size is not a true paleothermometer, as it is correlated with climatic transitions through the dust content of the ice.

Abstract: By using recently developed automatic instruments for fabric and texture measurements on ice, both the c-axis orientation and area of the individual crystals can be determined. Each grain can then be associated with its volume fraction, defined as a function of its measured cross-sectional area, to describe the microstructure of a polycrystal. The relevance of this approach is studied using a three-dimensional microstructure obtained from the Potts model. In particular, the area weighting is compared to the classical implicit equal weighting used by glaciologists, which assumes that all the grains have the same volume fraction (discrete uniform distribution). Then, using the measurements of c-axis orientation and crystal size performed on the North Greenland Icecore Project (NorthGRIP) ice core, we compare area-weighted and equal-weighted fabrics. All these comparisons are made with respect to the orientation tensor. According to the ability of the Potts model to reproduce the ice microstructure, it is shown that using the grain cross-sectional area to infer its volume fraction improves the description of the actual polycrystal fabric.