Journal cover Journal topic
The Cryosphere An interactive open-access journal of the European Geosciences Union
Journal topic
TC | Articles | Volume 13, issue 5
The Cryosphere, 13, 1495–1511, 2019
© Author(s) 2019. This work is distributed under
the Creative Commons Attribution 4.0 License.
The Cryosphere, 13, 1495–1511, 2019
© Author(s) 2019. This work is distributed under
the Creative Commons Attribution 4.0 License.

Research article 27 May 2019

Research article | 27 May 2019

Recrystallization processes, microstructure and crystallographic preferred orientation evolution in polycrystalline ice during high-temperature simple shear

Recrystallization processes, microstructure and crystallographic preferred orientation evolution in polycrystalline ice during high-temperature simple shear
Baptiste Journaux1,2, Thomas Chauve2,3, Maurine Montagnat2, Andrea Tommasi4, Fabrice Barou4, David Mainprice4, and Léa Gest2 Baptiste Journaux et al.
  • 1Department of Earth and Space Sciences, University of Washington, NASA Astrobiology Institute, Seattle, USA
  • 2Université Grenoble Alpes, CNRS, IRD, G-INP, IGE, 38000 Grenoble, France
  • 3PGP, Department of Geoscience, University of Oslo, Oslo, Norway
  • 4Géosciences Montpellier, Université de Montpellier/CNRS, 34095 Montpellier, France

Correspondence: Baptiste Journaux ( and Maurine Montagnat (

Back to toptop

Torsion experiments were performed in polycrystalline ice at high temperature (0.97 Tm) to reproduce the simple shear kinematics that are believed to dominate in ice streams and at the base of fast-flowing glaciers. As clearly documented more than 30 years ago, under simple shear ice develops a two-maxima c axis crystallographic preferred orientation (CPO), which evolves rapidly into a single cluster CPO with a c axis perpendicular to the shear plane. Dynamic recrystallization mechanisms that occur in both laboratory conditions and naturally deformed ice are likely candidates to explain the observed CPO evolution. In this study, we use electron backscatter diffraction (EBSD) and automatic ice texture analyzer (AITA) to characterize the mechanisms accommodating deformation, the stress and strain heterogeneities that form under torsion of an initially isotropic polycrystalline ice sample at high temperature, and the role of dynamic recrystallization in accommodating these heterogeneities. These analyses highlight an interlocking microstructure, which results from heterogeneity-driven serrated grain boundary migration, and sub-grain boundaries composed of dislocations with a [c]-component Burgers vector, indicating that strong local stress heterogeneity develops, in particular, close to grain boundaries, even at high temperature and high finite shear strain. Based on these observations, we propose that nucleation by bulging, assisted by sub-grain boundary formation and followed by grain growth, is a very likely candidate to explain the progressive disappearance of the c axis CPO cluster at low angle to the shear plane and the stability of the one normal to it. We therefore strongly support the development of new polycrystal plasticity models limiting dislocation slip on non-basal slip systems and allowing for efficient accommodation of strain incompatibilities by an association of bulging and formation of sub-grain boundaries with a significant [c] component.

1 Introduction
Back to toptop

Ice deforms by shear in many natural conditions such as glaciers and ice sheets, and in particular along ice streams. In most of the deep ice cores studied, compression and extension are the dominant deformation geometries owing to the core location close to geographical domes or along ridges. Simple shear occurs in the deeper portions of these cores. However, the strongest simple shear occurs in fast-flowing glaciers and ice streams (see review by Hudleston2015). During large-scale ice flow, deformation induces development of strong crystallographic preferred orientation (CPO) or texture. Since the late 1970s a growing number of studies have provided measurements of increasing accuracy of the CPO evolution along ice cores (refer to Gow and Williamson1976; Alley1988; Lipenkov et al.1989, for pioneer work, and Faria et al.2014b, for a review). Pioneer work of Hudleston (1977) analyzed the CPO evolution in a narrow shear zone of the Barnes Ice Cap (Canada), providing well-documented observations of CPO development owing to shear in natural conditions.

CPO development in ice is related to the strong viscoplastic anisotropy of this material (Alley1992). Ice Ih is hexagonal with Space group P63mmc, no. 194 and Laue class 6/mmm. Plastic strain in ice Ih is mainly accommodated at high temperature by dislocation glide on the basal plane through three 1/3<1120>(0001) equivalent slip systems (Duval et al.1983; Hondoh2000). The viscoplastic anisotropy of the ice crystal results in strong strain heterogeneities at the inter- and intragranular scale when polycrystalline ice is submitted to a macroscopic stress (Duval et al.1983; Grennerat et al.2012, and references therein) and hence in activation of dynamic recrystallization processes, even at very low strains (Duval1979; Duval et al.1983; De La Chapelle et al.1998; Chauve et al.2017a).

In return, the CPO-induced anisotropy at the macroscale leads to a bulk anisotropic viscosity of ice that impacts ice rheology. This macroscale mechanical anisotropy is responsible for modifications in ice stratigraphy and age–depth relationship (see, e.g., Martín et al.2009; Buiron et al.2011). It may also produce large-scale flow heterogeneities, such as basal folding (see, e.g., Dahl-Jensen et al.1997; Thorsteinsson and Waddington2002; Gillet-Chaulet et al.2006; Ma et al.2010; Bons et al.2016). In simple shear the evolution of CPO (as well as other factors like the grain size) is associated with an increase in shear strain rate by a factor of 10 at constant stress (Cuffey and Paterson2010; Hudleston2015).

This phenomenon is well-studied in the Earth mantle, where olivine CPOs produced by similar deformation and recrystallization mechanisms (see reviews by Ismaïl and Mainprice1998; Tommasi and Vauchez2015) result in elastic and viscoplastic anisotropy. The elastic anisotropy in the mantle may be directly observed in seismic data (see Tommasi and Vauchez2015). The viscoplastic anisotropy is only observed indirectly, via the changes it produces in the deformation patterns. It nevertheless plays a major role in the development of plate tectonics (Tommasi et al.2009).

In addition to the obvious implications of ice rheology studies for ice caps and glacier dynamics, ice is also often used as a reference material to study anisotropic mechanical behaviors in rocks and alloys. Since ductile deformation experiments in rocks are challenging, ice is a good analogue to polycrystalline geo-materials since “high-temperature” deformation in ice at TTm close to 1 can be achieved in a cold room and at ambient pressure. Mantle rocks deformed close to the melting point in the asthenosphere beneath mid-ocean ridges, sampled in ophiolites, display strong olivine CPO due to deformation and dynamic recrystallization with fast grain boundary migration (e.g., Cassard et al.1981; Boudier and Coleman1981; Higgie and Tommasi2012). The mechanical anisotropy due to these CPOs may significantly change the flow patterns beneath the ridge (Blackman et al.2017).

Pioneering experimental studies have shown that at low shear strains (γ<0.6), an initially random CPO polycrystalline ice sample develops a bimodal distribution of the grain's c axes (Kamb1972; Duval1981; Bouchez and Duval1982). These are noted as the M1 sub-maximum, which remains normal to the macroscopic shear plane, and the M2 sub-maximum, which rotates from near parallel to the shear plane towards the M1 maximum with increasing finite shear strain (Bouchez and Duval1982). For important shear strains (γ≥2), only the M1 maximum persists. This orientation corresponds to crystals with their basal planes well-oriented for dislocation glide parallel to the shear plane. A similar evolution was observed in more recent shear experiments on artificial ice polycrystals by Li et al. (2000) and Budd et al. (2013), as well as by Wilson and Peternell (2012), who analyzed the influence of the initial CPO and the importance of recrystallization processes on the CPO evolution. The vanishing of the M2 maximum has been hypothesized by Bouchez and Duval (1982) to result from rigid body rotation and recrystallization of grains unfavorably oriented for slip.

A large part of the knowledge on the microscopic processes occurring in polycrystalline ice under simple shear deformation is still derived from data published over 30 years ago (Kamb1959, 1972; Duval1981; Bouchez and Duval1982; Burg et al.1986). The tools and methods used in these studies to analyze the CPO were often manual and highly dependent on the operator experience. Electron backscatter diffraction (EBSD) and automatic ice CPO analyzer (AITA) can now provide high-spatial- and high-angular-resolution quantitative data, enabling a global and statistical study of the processes accommodating strain at the micro-scale. Recent experiments (Qi et al.2017, 2019) using these new characterization techniques have shed new light on some aspects of the question. They have, for instance, disproved the hypothesis by Kamb (1972) that CPO evolution in ice mainly depends on the finite shear strain and is not sensitive to temperature, strain rate or stress. Indeed, Qi et al. (2017) showed that during axial compression the final CPO is sensitive to stress or strain rate and Qi et al. (2019) that the rate of evolution of the CPO in simple shear is sensitive to temperature.

A modern and comprehensive update of the observations of ice deformation in simple shear seems more and more necessary as the increase in calculation power allows us to run mean-field or full-field polycrystal plasticity models routinely (see Montagnat et al.2014, for a review and Llorens et al.2017, for a recent application). These new numerical approaches require accurate experimental constraints on the microstructure evolution. It is worth noting that these models are still unable to accurately reproduce the double sub-maximum of the c axis preferred orientation evolution of polycrystalline ice observed in simple shear experiments (Bouchez and Duval1982; Burg et al.1986; Wilson and Peternell2012) and natural shear zone (Hudleston1977).

The complex interplay of dislocation slip system activities, their interactions, and the role of the various mechanisms accommodating the local strain and stress heterogeneities, such as cross-slip, kink-banding and dynamic recrystallization, is still either not or poorly represented in most of the existing modeling approaches, very likely because of the complexity of the interactions between these mechanisms. Concerning the olivine CPO development in the mantle, Signorelli and Tommasi (2015) succeeded in reproducing the evolution of olivine CPO in simple shear with a mean-field modeling approach, by making use of a two-level mechanical interaction scheme that integrates the formation of low-angle grain boundaries as expected during dynamic recrystallization.

The present work aims at providing new constraints on the processes controlling CPO evolution in simple shear in polycrystalline ice by presenting a comprehensive description of the evolution of CPO and grain shape fabric in polycrystalline ice during hot torsion experiments. We use state of the art analytical techniques: EBSD equipped with a cryostage and an automatic ice texture analyzer (AITA), to retrieve quantitative data of unprecedented spatial and angular resolution on crystallographic orientations, grain size and grain shape fabrics as a function of shear strain. These data then are used to estimate the active dislocation systems and to investigate the role of sub-grain-scale processes on the recrystallization regime that allows the formation of a CPO and microstructure more favorable to shear strain.

2 Experimental methods
Back to toptop

2.1 Polycrystalline ice preparation

Unstrained equiaxial polycrystalline ice was prepared by evenly packing 200 µm sieved ice particles in a mold. The container was sealed and the air contained in the porosity was removed using a primary vacuum rotary vane oil pump. Outgassed water at 0 C was quickly introduced afterwards to fill the porosity. The resulting mush was frozen bottom up at −5C for 48 h. To enable a more homogenous grain size and recovery of any small strain accumulated during grain growth, the granular ice cylinders were annealed for at least 120 h at −7C. This method produces untextured homogenous samples with negligible porosity and ∼1.5 mm mean grain size. Ultra pure Milli-Q water (18 MΩcm) was used for all sample preparations. Cylinders around 35 mm in diameter and 60 mm in height were carved from the granular ice using a lathe located in a −10C cold room. Specific sample sizes are indicated in Table 1.

Table 1Torsion experiment details for each sample. The sample that was annealed is indicated with an asterisk. Mmax is the calculated torque, τmax is the maximum calculated shear stress, εmax is the maximum measured strain, and γ˙f and ε˙f are the final measures shear strain rate and strain rate.

Download Print Version | Download XLSX

2.2 Torsion experiments

Deformation experiments were performed using a torsion apparatus placed inside a cold room at a temperature of -7±0.5C (0.97Tm, with Tm the melting temperature at 0 C). Samples were carefully fixed in the apparatus ensuring no axial load. The design of the torsion apparatus does not allow for displacements parallel to the rotation axis; the imposed deformation is therefore perfect simple shear. During the experiments, the evolution of the CPO under these fixed-end boundary conditions might produce axial stresses (Swift1947; Li et al.2000). The latter cannot be measured in the present apparatus, but polycrystal plasticity models indicate that these axial stresses may attain values similar to those of the shear stresses when the CPO is oblique to the imposed shear (Castelnau et al.1996). A more precise description of the apparatus is provided in the Supplement. Straight lines parallel to the rotation axis were traced on the surface of the sample to allow for observation of any strain heterogeneity along the cylinder (Fig. 1). The surface of the sample was coated with silicon grease and wrapped gently with a cellophane foil to ensure negligible loss by sublimation.

We imposed a constant torque M on one site of the sample corresponding to a maximum stress τmax between 0.4 and 0.6 MPa (Table 1). The τmax was calculated from the torque M using a stress–strain power law as described in Paterson and Olgaard (2000). Since most of the present experiments recorded secondary creep conditions (Fig. 1), a stress exponent of n=3 was chosen based on results from Duval et al. (1983) and more recent work from Treverrow et al. (2012) that observed stress exponents between 2.9 and 3.1 during the secondary creep regime, and of 3.5 for tertiary creep in compression and shear tests. The actual torque applied to each sample is indicated in Table 1. These values of maximum tangential stress have been chosen because they allow high enough strain rates while preventing the formation of cracks.

Macroscopic strain (Fig. 1) was followed using a linear variable differential transformer (LVDT) during the primary creep regime (i.e., small strains) and a homemade optical rotary encoder device (5000 counts per revolution) for larger strains. Maximum macroscopic shear strain γ, the maximum natural extensional strain ε, and the strain rate ε˙, calculated through


(Paterson and Olgaard2000), are reported in Table 1.

2.3 Analytical methods

CPO, microstructure and sub-grain strain heterogeneity analyses were performed using both AITA and EBSD, as the combination of the two techniques allows the acquisition of a comprehensive dataset. AITA enables the analyses of large samples (up to 10×10 cm2), but only recovers the c axis (optical axis) crystallographic orientations. Conversely, the cryo-EBSD setup used in the present study (CamScan X500FE CrystalProbe at Géosciences Montpellier) is only capable of mapping smaller ice samples of 1×2 cm2, but is able to recover the full orientation of the crystal (all crystallographic axes) with an expected angular resolution better than 0.7 (Randle1992) and with a practical spatial resolution for ice down to 5 µm. Sections tangential to the experimental cylinders were cut out for both AITA and EBSD analyses. The general geometry of these sections is illustrated in Fig. 2. To provide flat surfaces for AITA thin sections and EBSD samples, the sample sections were microtomed by less than 1 mm, which gives an error on γmax of less than 5 %.

AITA has been performed on thin sections typically 300 µm thick, as is routinely carried out on ice (Russell-Head and Wilson1999; Montagnat et al.2011; Peternell et al.2011; Chauve et al.2015; Wilson et al.2015). We worked at a spatial resolution of 25 µm and were able to obtain an angular resolution of about ∼3. All optical CPO measurements were conducted at −7C. Grain boundaries (GBs) were extracted using the segmentation technique described in Montagnat et al. (2015), which automatically defines grain boundaries using color change detection in AITA maps, followed by manual corrections based on coherent microstructure shapes and discernible misorientations. All manual corrections are operator dependent, so all analyses have been performed by the first author for consistency. This allowed us to differentiate, with a good level of confidence, sub-grain boundaries (SGBs) from grain boundaries, even if optical measurements do not give access to the full crystallographic orientation, as EBSD measurements do. From the extracted skeleton images we derived grain size and shape statistics using the PolyLX MATLAB® toolbox developed by Lexa (2003). Grain size frequency plots are reported for every finite shear strain in the Supplement. From the grain shapes, we calculated the PARIS grain shape factor to quantify grain boundary interlocking as PARIS =2(P-PE)/PE100 %, with P the perimeter and PE the convex hull perimeter (Heilbronner and Barrett2014). This parameter increases with increasing grain boundary sinuosity.

Electron backscatter diffraction (EBSD) patterns were obtained with a CamScan X500FE CrystalProbe equipped with a liquid nitrogen cryogenic stage (GATAN®) at Géosciences Montpellier (France). Sample surfaces were prepared by carefully flattening the surface with a microtome blade at −60C inside a freezer placed near the EBSD just before loading it quickly on the cold stage. This method provides a very good surface finish with no significant frost during measurements and limits to a minimum the heat transfer to the sample, ensuring limited microstructure evolution due to annealing. The CrystalProbe–EBSD working conditions were 15 kV, 3.5 nA, and a low vacuum with 1 Pa of gaseous nitrogen and 25 mm working distance. More details on the EBSD geometry can be found in Demouchy et al. (2011). To avoid ice sample sublimation in the vacuum chamber, which happens at approximately −60C at 1 Pa for ice Ih (Fernicola et al.2011), the temperature of the cold stage was kept at -100±10C during the analyses. EBSD Kikuchi patterns were automatically indexed using the HKL Channel5 software suite (Oxford Instruments®). Maps were acquired with a spatial resolution of 25 µm identified as a good compromise between resolution and time constrains, with an indexation rate higher than 85 %. Analyses of the EBSD maps were made using the MTEX MATLAB® toolbox (Hielscher and Schaeben2008; Bachmann et al.2010, 2011; Mainprice et al.2014).

To observe the geometry of sub-grain boundaries and extract statistical quantities, local misorientation analyses, such as kernel average misorientation (KAM) or the local misorientation relative to the mean orientation (Mis2Mean), were performed on the EBSD data (see review by Wright et al.2011). Both analyses enable visualization of the spatial distribution of intragranular stored strain structures, generally in the form of sub-grain boundaries (Wright et al.2011), but do not provide any information about the type of dislocations accommodating the observed misorientation. Conversely, the weighted Burgers vector (WBV) analysis allows us to extract a minimum value for the dislocation density in sub-grain boundaries and provides an estimate of the slip systems necessary to explain the observed misorientations (Wheeler et al.2009). As conventional EBSD maps are 2-D, only five components of the Nye tensor α can be extracted (α12, α21, α13, α23, α33). This allows us to obtain a projection of the Nye tensor on the EBSD map plane which partially constrains the Burgers vector. It provides a lower bound for the dislocation density and constrains the Burgers vector of geometrically necessary dislocations (GNDs) accommodating the observed misorientations. Even if WBV analysis does not recover the full Burger vectors coordinates, it does not provide “false” components. This means that if a [c] component is deduced from WBV analysis of a 2-D EBSD map, the true Burgers vector does actually contain at least this amount of [c] component. For the WBV analysis we used the MATLAB® toolbox developed by Wheeler et al. (2009) and previously used on ice in Chauve et al. (2017b). To provide the relative number of [c]-component dislocations, the calculated WBV for each pixel is projected on the non-independent four hexagonal lattice directions ([1120], [2110], [1210], [0001]), respectively called WBVa1, WBVa2, WBVa3 and WBVc. The ratio of [c] component (WBVc) is then given by the relation rWBVc=(WBVc/WBV), with WBV the Euclidian norm of the WBV. A high-pass value of 0.4 (WBV>2.8×10-4µm−1) was imposed for the misorientation analysis to filter the noise resulting from the angular resolution of the EBSD data. We suggest referring to Appendix A of Chauve et al. (2017b) for more details about this analysis. Nevertheless, it should be recalled that the WBV analysis does not directly provide information about mobile dislocations, responsible for the majority of the plastic deformation. It only measures the GND, as it is based upon the Nye tensor.

3 Results
Back to toptop

Five torsion experiments were performed up to different finite shear strains γmax, from 0.012 to 1.94. Data for all experiments are presented in Table 1. One sample (TGI0.71), with a γmax of 0.71, was annealed for 72 h at −7 after the deformation experiment to study the effect of annealing on both CPO and microstructure. The microstructure of one of these samples has been partially described in Chauve et al. (2017b) to compare the statistical representation of WBV (or Burgers vector) with [c] component in GNDs in samples submitted to uniaxial unconfined compression and torsion.

3.1 Macroscopic strain

Retrieved samples presented homogeneous strain with no apparent localization or cracks, as illustrated in Fig. 1. TGI0.012 and TGI0.2 shear strain data were very noisy due to signal transmission issues on the acquisition setup, which were fixed for the TGI0.42 and TGI1.96 runs. Typical primary creep regime (i.e., hardening) is observed for TGI0.42 and TGI1.96 in Fig. 1 until approximately 20 000 s (≈6 h). The noise in TGI0.012 and TGI0.2 data was too strong to distinguish any primary creep hardening. After this phase of hardening, a minimum in strain rate is reached (secondary creep regime). This minimum is clearly visible in experiments TGI0.2, TGI0.42, TGI0.71 and TGI1.96 at 200 000 s (≈56 h). For the longer TGI1.96 experiment, after the secondary creep phase, strain rate continuously increased until achieving an almost steady-state strain rate of 8.5×10-7 s−1, suggesting that tertiary creep regime was reached for this run (Fig. 1). TGI0.012 stayed mostly in the primary creep regime with a γmax=0.012. Samples TGI0.2 and TGI0.42 experienced a maximum shear strain of 0.2 and 0.42, respectively. The final strain rate for these two experiments is illustrated by the tangent dashed line in Fig. 1, which shows that softening was still occurring when the samples were unloaded (only the last data points lie on the tangent). We therefore consider that those experiments were still in a transient regime when stopped. A strain step is visible for TGI0.42 between 550 000 s (≈153 h) and 650 000 s (≈180 h); it corresponds to a period where the weight used to apply constant torque was in contact with the metallic frame of the torsion apparatus, consequently decreasing the applied stress. After 650 000 s (≈180 h), the weight was put back in a good configuration and the strain curve recovered its pre-step slope. The significant variations observed in strain rate evolution with time between the different runs cannot be attributed to a variation in initial grain size, CPOs or in the applied torque alone (Table 1), but rather to coarse-grained microstructure of the samples, which resulted in fewer than 25 grains per radii and, hence, in a not perfectly isotropic macroscopic behavior. The strain–time curves presented in Fig. 1 are therefore useful to characterize each run creep regime independently, but should be used with care in comparison between different samples or with other experiments.

Figure 1Creep curves for TGI0.2, TGI0.42, TGI0.71 and TGI1.96 experiments with the tangent to end of curve as black dashed lines, and corresponding final strain rate ε˙f. The blank parts of some of the curves correspond to a loss of signal between the acquisition setup and the sensor located in the cold room. In the top left frame is a zoom on the primary creep regime (<10 h) for better readability. The photograph in the bottom right shows the TGI1.96 sample after unloading, with the inclined ink color stripes, initially straight, used as strain markers.


3.2 Microstructure evolution

AITA maps of tangential sections for each sample and associated analyses are presented in Fig. 2 and Table 2. The evolution of CPO and microstructure (grain size and shape) as a function of finite shear strain is clearly visible. The color wheel gives the orientation of the <0001> c axis. The annealed sample TGI0.71 is discussed separately at the end of this section. As the shear strain increases, the distribution of colors changes from random color toward a predominance of red and blue grains, until the map becomes dominated by red grains for TGI1.96 (γmax=1.96).

Figure 2(a) AITA analyses of tangential sections for every sample showing the evolution of CPO and microstructure, with x parallel to the shear plane and y parallel to the torsion axis. The scheme on the upper right shows the geometrical relation between the tangential section and the retrieved cylindrical samples. Arrows indicate the shear directions. Color-coded c axis orientation is indicated by the color wheel. (b) Pole figures (equal area) showing the c axis orientations evolution with shear strain. Shear direction, direction of finite extension (ED) and compression direction (CD) are indicated with a black line, dotted black line and grey dashed line, respectively. Direction angles for the M1 and M2 sub-maxima are reported on the pole figures with red and green triangles, respectively, along with their angles with the ED direction. Theoretical perfectly symmetric M2 poles to M1 with ED direction (M2S) are represented in orange outline open triangles. J index (Bunge1982) and φ angle between M1 and M2 sub-maxima are indicated below for γmax≥0.2. (c) Rose diagram of grain shape fabrics. Mean orientation is indicated as a red line and the 95 % confidence interval as a red arc. Finite shear strain ellipses are plotted onto each rose diagram. The total number of grains segmented is reported as ng.


Table 2Grain size and shape analysis results based on AITA measurements. The sample that was annealed is indicated by an asterisk. ng is the number of grains used for the analyses and dg the median of the grain size (with associated interquartile range). Standard deviations are indicated in parentheses.

Download Print Version | Download XLSX

This evolution records the progressive change in orientation of the c axis with increasing shear strain, displayed as pole figures in Fig. 2b. For the undeformed ice sample (γmax=0) we observe approximately random orientations of c axes with a texture (CPO) index Ji=1.07. The c axis texture (CPO) J index is derived by calculating the second moment of the Euler angles' orientation distribution function (ODF) of discrete crystal orientation data (Bunge1982). It is commonly used to quantify CPO strength in geology and in material science. The strain reached by TGI0.012 sample is not high enough to produce any significant change in the CPO. As the strain increases we observe the formation of two sub-maxima, as reported in previous studies (Kamb1959, 1972; Duval1981; Bouchez and Duval1982; Burg et al.1986), one almost perpendicular to the shear plane (M1) and the other one at a low angle (initially 11–17) from the shear plane (M2) as indicated in Fig. 2. With increasing shear, M2 rotates into higher angles to the shear plane (up to 47) while M1 remains stable at 90 from the macroscopic shear plane. More importantly, the maxima M1 and M2 increase and decrease in intensity, respectively. This evolution towards a single maximum CPO results in an increase in the J index (Fig. 3). The plane containing the finite extension direction (ED) and normal to the compression direction (CD) has a very low [c] axis CPO density for γ≥0.2. One should note that for 0.2γ0.42, the bimodal pattern formed by the two M1 and M2 sub-maxima is almost symmetric relative to the maximum finite extension axis (ED).

Nevertheless the symmetry between M1 and M2 around the finite extension direction is not perfect. The angle between M2 and ED is generally larger than the angle between M1 and ED by 1 to 3. The exception is the annealed sample TGI0.71, where the difference between the two is 17.4, with M2 closer to ED than it should have been for a perfect symmetry. This change may be due to a post-deformation CPO evolution, due to grain growth during annealing. A small lag in the reorientation of the M2 sub-maximum relative to the M1 sub-maximum is also observed in other simple shear experiments (Bouchez and Duval1982; Qi et al.2019). The limited number of experiments performed in the present study precludes a statistical analysis of this behavior. The evolution of the angle between M1 and M2 (φ) with increasing shear strain is discussed in more detail in Qi et al. (2019), who compared observations from ice shear experiments carried out at different temperatures and numerical modeling.

AITA-extracted grain statistics are presented in Table 2. With increasing γ the median grain size (with associated interquartile range) diminishes from 0.7(4) to 0.12(16) mm. The median is preferred relative to the mean because of the non-normal distribution of grain sizes. The coarser grain size in the TGI0.42 sample compared to TGI0.2 and TGI0.012 is probably due to a difference in the initial grain size, but the grain sizes remain comparable if interquartile ranges are considered. Grain shape fabrics calculated using moment-based ellipse fitting technics (Lexa2003) are illustrated as rose diagrams in Fig. 2c. As the shear strain increases, we observe the development of a grain shape fabric following the maximum finite elongation direction indicated with the corresponding overlying strain ellipse in Fig. 2c. The grain shape fabric clearly evolves until γ=0.42 and then does not vary significantly up to γ=1.96

Grain boundary interlocking estimated through the PARIS factor (reported in Table 2) increases during the transient strain regime (from 2.2 % at γ=0 to 11.5 % at γ=0.42) and then decreases. At the highest shear strain of γ=1.96, the PARIS (5.9 %) is half of the one at a shear strain of 0.42.

The annealed sample TGI0.71 shows a median grain size of 1.7 mm, above the initial unstrained value of 0.7 (Table 2). The CPO does not seem to be modified by the annealing with the M1 and M2 maxima maintained and a J index =2.96, which lies in between the values for γmax 0.42 and 1.96, as would be expected from a non-annealed sample with the same finite shear strain (Fig. 2b). The grain shape fabric in the annealed sample is nevertheless more marked than in the strained but non-annealed samples. The PARIS factor decreases significantly, as much of the interlocking of grains is smoothed by annealing.

In summary, the crystallographic orientation evolution during simple shear strain of polycrystalline ice can be described as follows. M1 orientations rapidly develop (by γ=0.2). They remain stable in orientation and increase in intensity until they become predominant at high strains (γ>1.96). The M2 maximum also rapidly develops; it is roughly symmetric to M1 with respect to the maximum finite shortening direction. However the M2 maximum continuously decreases in intensity during the transient regime. The nearly complete disappearance of M2 correlates with achievement of steady state (tertiary creep). All other crystallographic orientations start to disappear as early as γ=0.2 and this tendency is accentuated with increasing shear strain. Finally, if annealing occurs after the strain stops, the CPO is maintained even if the grain shape fabric and grain size are modified.

Figure 3c axis texture (CPO) index Ji determined from AITA and EBSD data for all the shear strains investigated. Annealed sample TGI0.71 is highlighted for clarity.


3.3 EBSD

3.3.1 Maps and textures

Multiple EBSD maps were acquired for every final shear strain. Representative examples of EBSD maps for each sample are shown in Fig. 4a, with the number of maps nm acquired for each γmax indicated at the bottom. The access to the full crystallographic orientations enables more complete grain segmentation based on misorientation data. Grain boundaries are characterized by misorientations ≥5 and sub-grain boundaries between 1 and 5 (Chauve et al.2017b).

The orientations of <0001><1010> and <1120> axes relative to the shear plane and direction are shown in Fig. 4b. All maps acquired for each γmax were combined for better statistical analysis. An excellent level of agreement with CPOs from AITA is observed for <0001> CPOs, which rapidly develop M1 and M2 sub-maxima, with disappearance of the M2 maximum and reinforcement of the M1 one at high shear strain. The increase in c axis CPO J index confirms this evolution, with very consistent values with those based on the AITA measurements for a same γmax (Fig. 3). For the high shear strain sample at γmax=1.96, the <1010> and <1120> axes (<a> and <m> axis) form a girdle, which tends to align in the shear plane. Within this girdle, there is a preferred orientation of both <a> and <m> directions toward the shear direction. The present CPO is similar to those formed in direct shear experiments (Qi et al.2019). It is consistent with equivalent contribution of the three <a> axes in accommodating shear on the (0001) plane.

Some elongation of the distribution of the M1 and M2 sub-maxima towards the z direction (Fig. 4b), which is the normal to the shear direction in the shear plane, is visible in our results. This elongation is best expressed for the M1 maximum in the highest-strain sample TG1.96, for which pole figures for <0001><1010> and <1120> lattice vectors are now represented in two perpendicular reference frames in Fig. 4b for better readability. Similar elongated distributions of <c> axes have been reported in direct shear experiments by Qi et al. (2019). Some elongation of the M1 maximum is also observed in the highest shear strain sample (γ=2) of Bouchez and Duval (1982) as well as in other shear experiments in Li et al. (2000), Wilson and Peternell (2012), and Budd et al. (2013). However, most naturally sheared ice samples do not have elongated <c> maxima (Hudleston1977).

Figure 4(a) Example of EBSD patterns obtained from tangential sections for each maximal shear strain, at −7C. Map colors correspond to the orientation of the crystallographic lattice toward to the y axis (torsion axis), as indicated in the inverse pole figure on the TGI0.012 EBSD map (top left). Red grains have <0001> (c axis) parallel to the torsion axis, green grains have <1010> (a axis) parallel to the torsion axis and blue grains have <1120> (m axis). Boundaries with misorientation greater than 5 are represented by black lines and boundaries between 1 and 5 are depicted as white lines. (b) Lower hemisphere pole figures (equal area) representing the orientations of <0001>, <1010> and <1120> directions of the EBSD data were combined for each γmax. Additional pole figures projected normal to the x direction are represented for γmax=1.96. Orientation density intensities are scaled to their maximum value for better readability. Associated c axis texture (CPO) J index is indicated above the pole figure. The shear plane (SP) is indicated as a black line, and the finite extension (ED) and compression (CD) directions are represented by black and grey dashed lines, respectively. The number of maps used to obtain the pole figures is indicated as nm and the total number of grains segmented as ng.


3.3.2 Weighted Burgers vector analysis and intragranular strain

During plastic strain, some of the stress on the ice lattice is relieved by the formation and glide of dislocations. Some of those dislocations are trapped in the crystal lattice, forming stable arrays of geometrically necessary dislocations (GNDs), which may define sub-grain boundaries. Such intragranular strain heterogeneities are observed in numerous experimental studies on polycrystalline ice and in natural environments, (Wilson et al.1986; Mansuy et al.2000; Piazolo et al.2008; Montagnat et al.2011; Weikusat et al.2011; Faria et al.2014a; Chauve et al.2017b). As described in the methods section, WBV analysis allows the calculation of a minimum density of GNDs and definition of the relative proportion of dislocations with <a> and [c]-component Burgers vector within this GND population. Figure 5 provides representative WBV maps at every finite shear strain investigated.

Figure 5Maps of the norm of the weighted Burgers vector, WBV. (a) Map examples for every shear strain. Zoomed in region in (b) delimited by black dashed rectangles. (b) Details on the TGI0.2, TGI0.42 and TGI1.96. Grains boundaries (>5) are reported as black lines. Non-indexed regions reported in grey. Pole figures representing c axis orientations are reported at the top right.


In these WBV maps (Fig. 5) we observe a large number of linear sub-grain boundaries already developed at γmax=0.2. Evolution of the GNDs density with increasing strain is hard to assess by visually comparing the maps for different finite shear strains, but it can be assessed by the analysis of the quantitative statistical distribution of <a> and [c] GNDs (Fig. 6). The sub-grain boundaries are mostly organized in walls and the direction of the trace of these planes on the thin section typically corresponds to the parent grain [c]-axis direction. This kind of sub-grain boundary, very common in both natural and experimentally deformed ice, has been identified as a basal tilt boundary composed of basal edge dislocations (b=a, with b the Burgers vector) through EBSD and X-ray Laue diffraction (Piazolo et al.2008; Weikusat et al.2011; Chauve et al.2017b). These sub-grain boundaries are visible in Fig. 5b as purple and blue lines. Other types of low-angle boundaries were also identified in both natural and experimentally deformed ice; these are basal twist boundaries with basal screw dislocation and non-basal tilt boundaries with variable numbers of dislocations with a [c] axis Burgers vector component (i.e., <a+c>