Flow-induced anisotropy in polar ice and related ice-sheet flow modelling

• Published in: Journal of Non-Newtonian Fluid Mechanics
• Main Authors: Gillet-Chaulet, Fabien, GAGLIARDINI, Olivier, Meyssonnier, Jacques, Zwinger, Thomas, Ruokolainen, Juha
• Other Authors: Laboratoire de glaciologie et géophysique de l'environnement (LGGE), Centre National de la Recherche Scientifique (CNRS)-Institut national des sciences de l'Univers (INSU - CNRS)-Observatoire des Sciences de l'Univers de Grenoble (OSUG), Université Joseph Fourier - Grenoble 1 (UJF)-Institut polytechnique de Grenoble - Grenoble Institute of Technology (Grenoble INP )-Institut national des sciences de l'Univers (INSU - CNRS)-Institut national de recherche en sciences et technologies pour l'environnement et l'agriculture (IRSTEA)-Université Savoie Mont Blanc (USMB Université de Savoie Université de Chambéry )-Centre National de la Recherche Scientifique (CNRS)-Université Joseph Fourier - Grenoble 1 (UJF)-Institut polytechnique de Grenoble - Grenoble Institute of Technology (Grenoble INP )-Institut national de recherche en sciences et technologies pour l'environnement et l'agriculture (IRSTEA)-Université Savoie Mont Blanc (USMB Université de Savoie Université de Chambéry )-Centre National de la Recherche Scientifique (CNRS), Scientific Computing Ltd (CSC), European Project for Ice Coring in Antarctica (EPICA), by European Commission under the Environment and Climate Program (contract ENV4-CT95-0074)
• Format: Article in Journal/Newspaper
• Language: English
• Published: HAL CCSD 2006
• Subjects: Polar ice anisotropy; Ice-sheet flow modelling; Invariant-based optimal fitting; [SDU.STU.GL]Sciences of the Universe [physics]/Earth Sciences/Glaciology; Ice Sheet
• Online Access: https://hal-insu.archives-ouvertes.fr/insu-00375483; https://doi.org/10.1016/j.jnnfm.2005.11.005

International audience As fibers or other crystalline materials exhibiting hexagonal symmetry, the crystal of ice can be orientated by using only one single vector, i.e. its c-axis. Such a characteristic allows to apply specific methods to deal with the properties of the polycrystalline aggregate. Among others, the fabric (texture) of the ice polycrystal can be described by an ODF, i.e. a scalar function of two angles that gives the distribution of the orientation of all the constituents (grains). This paper presents a strain-induced anisotropic flow law for polycrystalline ice and the associated equations describing the evolution of its fabric. This constitutive law is formulated at the polycrystal scale and tabulated using a micro–macro model. The fabric is defined by the second- and fourth-order orientation tensors for the c-axes, assuming the so-called “invariant-based optimal fitting closure approximation”. Both the anisotropic constitutive law and the fabric evolution equations have been implemented in a finite element code in order to solve large scale ice flow problem. As an application, the flow of an idealized ice sheet over a bumpy bed is studied.