Description: It has been debated for decades whether rigid inclusions, such as porphyroclasts and porphyroblasts, do or do not rotate in a softer matrix during deformation. Experiments and numerical simulations with viscous matrix rheologies show ongoing rotation of circular inclusions, whereas using Mohr-Coulomb plasticity results in nonrotation. Because the rocks in which inclusions are found normally undergo deformation by dislocation creep, we applied a full-field crystal plasticity approach to investigate the rotation behavior of rigid circular inclusions. We show that the inclusion's rotation strongly depends on the anisotropy of the matrix minerals. Strongly anisotropic minerals will develop shear bands that reduce the rotation of inclusions. Inhibition of rotation can only occur after a significant amount of strain. Our results may help to explain why geologic rigid objects often show evidence for rotation, but not necessarily in accordance with the viscous theory that is usually applied to these systems.

Description: Results of numerical simulations of co-axial deformation of pure ice up to high-strain, combining full-field modelling with recrystallisation are presented. Grain size and lattice preferred orientation analysis and comparisons between simulations at different strain-rates show how recrystallisation has a major effect on the microstructure, developing larger and equi-dimensional grains, but a relatively minor effect on the development of a preferred orientation of c-axes. Although c-axis distributions do not vary much, recrystallisation appears to have a distinct effect on the relative activities of slip systems, activating the pyramidal slip system and affecting the distribution of a-axes. The simulations reveal that the survival probability of individual grains is strongly related to the initial grain size, but only weakly dependent on hard or soft orientations with respect to the flow field. Dynamic recrystallisation reduces initial hardening, which is followed by a steady state characteristic of pure-shear deformation.

Description: We performed numerical simulations on the micro-dynamics of ice with air inclusions as a second phase. This provides first results of a numerical approach to model dynamic recrystallisation in polyphase crystalline aggregates. Our aim was to investigate the rheological effects of air inclusions and explain the onset of dynamic recrystallisation in the permeable firn. The simulations employ a full field theory crystal plasticity code coupled to codes simulating dynamic recrystallisation processes and predict time-resolved microstructure evolution in terms of lattice orientations, strain distribution, grain sizes and grain boundary network. Results show heterogeneous deformation throughout the simulations and indicate the importance of strain localisation controlled by air inclusions. This strain localisation gives rise to locally increased energies that drive dynamic recrystallisation and induce heterogeneous microstructures that are coherent with natural firn microstructures from EPICA Dronning Maud Land ice coring site in Antarctica. We conclude that although overall strains and stresses in firn are low, strain localisation associated with locally increased strain energies can explain the occurrence of dynamic recrystallisation.

Description: Understanding the flow of ice on the microstructural scale is essential for improving our knowledge of large-scale ice dynamics, and thus our ability to predict future changes of ice sheets. Polar ice behaves anisotropically during flow, which can lead to strain localisation. In order to study how dynamic recrystallisation affects to strain localisation in deep levels of polar ice sheets, we present a series of numerical simulations of ice polycrystals deformed under simple-shear conditions. The models explicitly simulate the evolution of microstructures using a full-field approach, based on the coupling of a viscoplastic deformation code (VPFFT) with dynamic recrystallisation codes. The simulations provide new insights into the distribution of stress, strain rate and lattice orientation fields with progressive strain, up to a shear strain of three. Our simulations show how the recrystallisation processes have a strong influence on the resulting microstructure (grain size and shape), while the development of lattice preferred orientations (LPO) appears to be less affected. Activation of non-basal slip systems is enhanced by recrystallisation and induces a strain hardening behaviour up to the onset of strain localisation and strain weakening behaviour. Simulations demonstrate that the strong intrinsic anisotropy of ice crystals is transferred to the polycrystalline scale and results in the development of strain localisation bands than can be masked by grain boundary migration. Therefore, the finite-strain history is non-directly reflected by the final microstructure. Masked strain localisation can be recognised in ice cores, such as the EDML, from the presence of stepped boundaries, microshear and grains with zig-zag geometries.

Description: Deformation localisation can lead to a variety of structures, such as shear zones and bands that range from grain to crustal scale, from discrete zones to anastomosing networks, and shear zone related folds. We present numerical simulations of the deformation of an intrinsically anisotropic material with a single maximum crystal preferred orientation (CPO) in simple shear. We use the Viscoplastic Full-Field Transform (VPFFT) crystal plasticity code coupled with the modelling platform ELLE to achieve very high strains. The VPFFT-approach simulates deformation by dislocation glide, taking into account the different available slip systems and their critical resolved shear stresses. We vary the anisotropy of the material from isotropic to highly anisotropic, as well as the orientation of the initial CPO. To visualize deformation structures, we use passive markers, for which we also systematically vary the initial orientation. At low strains the amount of strain rate localisation and resulting deformation structures depend on the initial CPO in all anisotropic models. Three regimes can be recognised: distributed shear localisation, synthetic shear bands and antithetic shear bands. However, at very high strains localisation behaviour always tends to converge to a similar state, independent of the initial CPO. Shear localisation is often detected by folded layers, which may be parallel to the anisotropy (e.g. cleavage formed by aligned mica), or the deformation of passive layering, such as original sedimentary layers. The resulting fold patterns vary strongly, depending on the original layer orientation. This can result in misleading structures that seem to indicate the opposite sense of shear.

Description: Deformation localisation can lead to a variety of structures, such as shear zones and shear bands (from grain to crustal scale), from isolated zones to anastomosing networks. The heterogeneous strain field can furthermore result in a wide range of highly diverse fold geometries. In anisotropic materials the deformation behaviour is controlled by the viscoplastic anisotropy of the material and their ability to form a crystallographic preferred orientation (CPO). Here we present a selection of a series of numerical simulations which aim to investigate (1) the influence of an initial CPO in an anisotropic material on localization behaviour and (2) the role of layering/passive markers on the development of deformation structures.

Description: We present a series of simple shear numerical simulations of dynamic recrystallization of two‐phase non‐linear viscous materials that represent temperate ice. Firstly, we investigate the effect of the presence of water on the resulting microstructures and, secondly, how water influences on P‐wave (Vp) and fast S‐wave (Vs) velocities. Regardless the water percentage, all simulations evolve from a random fabric to a vertical single maximum. For a purely solid aggregate, the highest Vp quickly aligns with the maximum c‐axis orientation. At the same time, the maximum c‐axis development reduces Vs in this orientation. When water is present, the developed maximum c‐axis orientation is less intense, which results in lower Vp and Vs. At high percentage of water, Vp does not align with the maximum c‐axis orientation. If the bulk modulus of ice is assumed for the water phase (i.e., implying that water is at high pressure), we find a remarkable decrease of Vs while Vp remains close to the value for purely solid ice. These results suggest that the decrease in Vs observed at the base of the ice sheets could be explained by the presence of water at elevated pressure, which would reside in isolated pockets at grain triple junctions. Under these conditions water would not favor sliding between ice grains. However, if we consider that deformation dominates over recrystallization water pockets get continuously stretched, allowing water films to be located at grain boundaries. This configuration would modify and even overprint the maximum c‐axis‐dependent orientation and the magnitude of seismic anisotropy.

Description: The ice mass balances of Antarctic and Greenland ice sheets represent the largest uncertainty for predicting future sea-level rise. Understanding how ice flows from the accumulation to the ablation zone is therefore crucial for correctly estimating the changing mass in polar ice-sheets. On Earth, ice crystals have a hexagonal symmetry (ice lh) with a strong anisotropy favouring basal slip. This results in a progressive development of a vertical c-axis preferred orientation (LPO) of ice polycrystalline aggregates during deformation. In depth, the elastic anisotropy of polycrystalline ice gradually increases due the development of a vertical LPO. Observations of P-wave (Vp) and S-wave (Vs) velocities in ice sheets reveal a strong decrease of ~25% of Vs in depth, while Vp remains approximately constant. According to Wittlinger and Farra (2015) the low Vs may be due to the presence of unfrozen liquids resulting from pre-melting at grain joints and/or melting of chemical solutions buried in ice. Although previous studies of two-phase rocks (including melt and water) show that seismic velocities depend on both LPO and water content, studies on the effect of melt on polar ice seismic velocity are scarce. In this contribution we investigate the changes in P- and faster S-wave velocities during deformation of polycrystalline ice with different melt fractions.

Description: We performed numerical simulations on the microdynamics of ice with air inclusions as a second phase. Our aim was to investigate the rheological effects of air inclusions and explain the onset of dynamic recrystallization in the permeable firn. The simulations employ a full-field theory crystal plasticity code coupled to codes simulating dynamic recrystallization processes and predict time-resolved microstructure evolution in terms of lattice orientations, strain distribution, grain sizes and grain-boundary network. Results show heterogeneous deformation throughout the simulations and indicate the importance of strain localization controlled by air inclusions. This strain localization gives rise to locally increased energies that drive dynamic recrystallization and induce heterogeneous microstructures that are coherent with natural firn microstructures from EPICA Dronning Maud Land ice coring site in Antarctica. We conclude that although overall strains and stresses in firn are low, strain localization associated with locally increased strain energies can explain the occurrence of dynamic recrystallization.