Introduction
Experimental deformation studies have documented that ice can produce distinct crystallographic preferred orientations (fabrics) very rapidly in regimes of compression at higher temperature (-10 to -1°C) (Reference WilsonWilson, 1982; Reference WilsonWilson and Russell-Head, 1982) or during simple shear (Reference Kamb, Heard, Borg, Carter and RaleighKamb, 1972; Reference BakerBaker, 1982; Reference Bouchez and DuvalBouchez and Duval, 1982). In polar ice sheets and in glaciers, where the temperature increases towards the bedrock contact (Reference Budd, Jenssen and RadokBudd and others, 1971; Reference RitzRitz, 1989), the deformation tends to localize into diffuse zones where the ice crystallographic properties change through a combination of compression and simple shear (Reference AlleyAlley, 1992). It is also well known that anisotropy in crystal properties and initial orientations govern the plasticity and microstructural evolution (Reference Kocks, Kocks, Tomé; and WenkKocks, 1998) of a material such as polycrystalline ice. Few experiments (e.g. Reference DuvalDuval, 1981; Reference Budd and JackaBudd and Jacka, 1989; references therein) document the evolution of ice fabrics in such regimes where the driving force for the localization of the deformation is shear stress. This produces a strain-induced anisotropy or foliation (Reference Hudleston, Saxena and BhattacharjiHudleston, 1977) where the ice grains develop a single-maximum c-axis fabric offset from the normal to the shear plane.
Such fabric development is produced by intracrystalline plastic deformation (Reference AlleyAlley, 1992) and achieved by slip on the basal (0001) plane (Reference SteinemannSteinemann, 1958) and along the (11 20) direction (Reference Glen, Jones and OuraGlen and Jones, 1967) and this depends on the orientation of the crystal. In randomly oriented polycrystalline ice (isotropic) deformed in compression, two-maxima c-axis fabrics are developed primarily in the plane of deformation at about 25-30° to the compression axis (Reference WilsonWilson, 1982). For simple shear after the initial formation of two maxima (Reference Burg, Wilson and MitchellBurg and others, 1986), the maximum perpendicular to the applied shear plane becomes dominant at higher shear strains (Reference Bouchez and DuvalBouchez and Duval, 1982). In strongly layered ice deformed at -2°C, Reference Wilson and SimWilson and Sim (2002) have shown that the layer orientation and pre-existing c-axis fabrics in a shear zone have little effect on the way dynamic recrystallization produces the final fabric, with little preservation of the earlier fabrics. In compression and simple shear experiments at -2°C, Reference Li and JackaLi and Jacka (1998) observed that the shear strain rate is in part dependent on the initial crystal-orientation fabric, with minimum strain rates in isotropic ice being 2.5 times lower than in ice with initial c-axis concentrations that lie parallel to the compression axis (i.e. easy-glide orientation).
Critical in the deformation of ice is the control of grain anisotropy, with an orientation dependence for the crystal to respond to the applied changes in the stress and strain rate. In this study, where ice crystals exhibit ‘easy’ slip only on the basal system (forming a set of two independent systems), it is expected that in common with other minerals (Reference Wenk, Canova, Molinari and KocksWenk and others, 1989) the flow stress will be lowest and strain rate highest when the dominant grain population lies close to an orientation where the resolved shear stress (Schmid factor) is a maximum with respect to the compressive stress (Reference Kocks, Kocks, Tomé; and WenkKocks, 1998). On the other hand, if the ice crystals lie in a ‘hard’- glide orientation it is expected that the strain rate will be lower. Similarly, for randomly oriented polycrystals with ‘equiaxed’ grains it would be expected that strain rates will lie between the grain aggregates dominated by large grain anisotropy with either hard or easy glide. Such intuitive models have been used to calculate the fabric evolution of a crystal aggregate (Reference AzumaAzuma, 1995; Reference Lebensohn and ToméLebensohn and Tomé, 1994; Reference ThorsteinssonThorsteinsson, 2002) and can lead to predictions for the macroscopic behaviour of an ice mass (Reference Budd and JackaBudd and Jacka, 1989; Reference PatersonPaterson, 1991).
Reference Li, Jacka and BuddLi and others (1996) discuss comparisons of the flow-law dependence of polycrystalline ice in terms of stress and strain rates. Their analysis is based on the assumption of a straightforward arithmetic average of the crystal properties to represent aggregate properties (e.g. Reference Li, Jacka and BuddLi and others, 1996, Equation (2)). In general, a straightforward average of the crystal properties is not correct, as averaging procedures ignore the local perturbations induced by the grain interaction with neighbours (Reference Tomé, Kocks, Tomé and WenkTomé, 1998). In addition, Reference Li and JackaLi and Jacka (1998) suggest the average bulk deformation of an ice aggregate is assumed to be a function of the ice’s initial preferred orientation. However, the compressibility of ice is a function of the stress components, the imposed strain field, temperature and the grain shape anisotropies.
Summary of experimental results for blocks of layered columnar ice and three samples of polycrystalline ice deformed under a combined compressive stress (σ = 0.22 MPa) and shear stress (τ = 0.4MPa). Table illustrates the variation in initial layer orientation, temperature, bulk compressive strain, bulk shear strain, minimum octahedral shear strain rate and increased octahedral shear strain rate observed at the termination of an experiment
We have conducted this study to further examine the combination of compression and shear in a range of temperatures from -20 to -2°C and compare our observations with the -2°C experiments of Reference Li and JackaLi and Jacka (1998) and the findings of Reference Budd and JackaBudd and Jacka (1989). As well as determining flow rates, we examine the accompanying changes in c-axis patterns produced in ice where there is a strong pre-existing grain structure and a population of elongate grains that parallel the glide plane (0001), which provides an anisotropy in the ice mass. The microfabric observations are based on results obtained using a fabric analyser (e.g. Reference Wilson, Russell-Head, Kunzi and ViolaWilson and others, 2007; Reference Wilson and PeternellWilson and Peternell, 2011). In conjunction with this, we investigate the differences in the evolution of microstructure and fabrics that can be attributed to different temperature regimes under the same combined stress configuration used by Reference Li, Jacka and BuddLi and others (1996). It is also shown that the shear strain rate is mainly temperature- dependent and not purely a consequence of the crystal- orientation fabric.
Experimental Procedure
The test samples (Table 1) were deformed and measured in a deformation apparatus similar to that used by Reference Li, Jacka and BuddLi and others (1996) and Reference Wilson, Russell-Head and SimWilson and others (2003). Experiments were performed under both constant load or compressive stress and constant shear stress conditions at -20 to -2°C (Table 1). The ice-test samples were laboratory-prepared and comprised (1) polycrystalline ice, with initially random (isotropic) crystal orientations and mean crystal area 1.6 mm2, or (2) layers (~10mm wide) of fine fibrous ice (2 mm × 5 mm) intergrown with elongate crystals (10mm × ~20mm) similar to that described by Reference Wilson and SimWilson and Sim (2002). The fine ice grains had c-axes that were randomly oriented and therefore offer strength isotropy, with multiple planes capable of undergoing easy glide, similar to a randomly oriented polycrystalline ice aggregate. This is in contrast to the anisotropy present in the elongate grains where the elongation parallels the basal plane, which is perpendicular to the c-axis. Layering in these samples was variably oriented with respect to the compression axis (x 1 in Table 1). This meant that the basal plane orientation in the elongate grains was in a hard-glide orientation when the layering was parallel to the compression axis (x 1), or close to the easy- glide orientation when oriented at ~40°. Polycrystalline ice samples were prepared in the manner described by Reference WilsonWilson and Russell-Head (1982).
The sample is frozen into two support boxes as illustrated by Reference Li, Jacka and BuddLi and others (1996), and the strain is not uniformly distributed, but localized to the centre of the sample (Fig. 1 zone of shear). There are transitional zones where the original layering is progressively rotated, and these are juxtaposed to relatively undeformed boundary zones. Deformation in such high-strain zones is usually noncoaxial, i.e. the principal axes of finite strain (x 1 x 2 x 3; Fig. 1) do not lie parallel to those of successive strain increments. In order to build on the results of Reference Li and JackaLi and Jacka (1998) and the experiments summarized by Reference Budd and JackaBudd and Jacka (1989), four series of experimental configurations were used: (1) A first set of tests was performed with initial layers coincident with the compression axis either parallel (Fig. 1a) or perpendicular (Fig. 1b) and perpendicular to the shear plane (x 2 x 3) over the temperature range -20 to -2°C (Figs 2 and 4, further below). (2) A set where the samples were then re-deformed (Fig. 1c) after a 90° rotation in the manner described by Reference Wilson, Russell-Head and SimWilson and others (2003) in order to evaluate the effect of a pre-existing anisotropy on the mechanical behaviour and fabric development. (3) A third set of tests at -2°C with randomly oriented polycrystalline ice and (4) layers inclined (Fig. 1d and e) to the compression axis (x 1) to evaluate the effect of the anisotropy on the mechanical behaviour (Fig. 3, further below). Many experiments were duplicated in order to verify the mechanical properties;however, some of the -10°C and -5°C samples melted or fractured before a proper microstructural analysis could be undertaken.
Initial sample configurations under a vertical compressive stress and anticlockwise shear stress. (a) Parallel columnar ice with layering inclined parallel to the x 1 x 2 plane. (b) Columnar ice with layering perpendicular to the x 1 x 2 plane as in experiment 2-16. (c) Composite sample with pre-deformed ice rotated 90º and inserted between ice wedges to have the same dimensions as in (a) in the manner described by Reference Wilson, Russell-Head and SimWilson and others (2003). (d) Parallel columnar ice with layering inclined 40º anticlockwise to the x 1 x 2 plane. (e) Parallel columnar ice with layering inclined 20º clockwise to the x 1 x 2 plane.
Creep curves comparing samples experimentally deformed at -2 to -20°C in combined compression and shear with a compressive stress of 0.22 MPa and shear stress of 0.4 MPa. (a) Creep curves for -2°C experiments 2-18, 2-36, 2-41 and the simple shear experiments of Reference Li and JackaLi and Jacka (1998) where A is a test on initially isotropic ice with random c-axis distributions, B is a test on a sample with a two-maxima c-axis distribution, and C is a test on a sample with single maximum of c-axes or grains that lie in an easy-glide orientation. (b) Creep curves for -5°C experiments 2-20 and 2-39, where layering is parallel to the x 1 x 2 plane, and 2-16 where layering is perpendicular to the x 1 x 2 plane. (c) Creep curves for -10°C parallel experiment 2-59 and composite samples 2-62 and 2-64. Also shown is composite -8.2°C experiment 2-62. (d) Creep curves for -15°C experiment 2-66 and -20°C experiment 2-71.
Creep curves for samples experimentally deformed at –5ºC (a) and –2ºC (a, b) in combined compression and shear with a compressive stress of 0.22MPa and shear stress of 0.4MPa. (a) Polycrystalline ice in experiments 2-19, 2-51 and 2-52. (b) Columnar ice with layering parallel (2-40), inclined 20º and 40º clockwise (2-46, 2-55) and 20º and 40º anticlockwise (2-43, 2-48) to the x1x2 plane.
Ice crystal statistics observed with the evolution of the recrystallized grain structure that accompanies a rise in temperature from –20ºC to –2ºC in four samples illustrated in the left column: (a) 2-71, (b) 2-66, (c), 2-64 and (d) 2-41. Ice grains from similar-sized areas of each sample were digitized manually and are labelled according to their grain sizes. Middle column shows grain-size histograms from each sample. Grain sizes are described by the crystal’s equal-area distance. n is number of measured crystals. A normal distribution function is plotted (red), and mean value and standard deviation are indicated. Right column shows a shape preferred-orientation analysis of the digitized grains as a rose diagram. n is number of measured crystals. Dotted circles: percentage of grains for each direction; red line: mean direction; red dotted line: circular standard deviation.
Orientation of c-axes was obtained using the fabric analysis apparatus described by Reference Peternell, Kohlmann, Wilson, Seiler and GleadowPeternell and others (2009), an instrument that determines the orientation of c-axes of uniaxial crystals at each pixel in the field of view (Reference Wilson, Russell-Head and SimWilson and others, 2003, Reference Wilson, Russell-Head, Kunzi and Viola2007). Retardation (highest birefringence colour), c-axis plunge, c-axis orientation trend, ‘Achsenverteilungsanalyse’ (AVA; axial distribution analysis) maps, grain boundary and quality images are produced in addition to the conventional crossed polar images (Fig. 5, further below). Data can be analysed individually using the G50 INVESTIGATOR software (Reference Wilson, Russell-Head, Kunzi and ViolaWilson and others, 2007). From this the operator can obtain either a bulk crystallographic preferred orientation, or manually select, using a cursor, the trend and plunge at any pixel site (2.8 pm maximum pixel resolution) in an individual grain. Based on the quality images (Reference Peternell, Kohlmann, Wilson, Seiler and GleadowPeternell and others, 2009), doubtful results can be rejected. This accounts for different numbers of c-axes seen in the automatically generated crystallographic preferred-orientation stereonets (Fig. 6, further below). A standard colour code (Fig. 6a) was used to identify the c-axis directions as azimuths and dips in the micrographs. In this scheme, north-trending axes are in red, east-trending axes are in green, west-trending axes are in blue, axes normal to the plane of the thin section are in black, and gradational shades represent angular differences between neighbouring pixels. This is a more precise visualization of orientation patterns than the 180° trends that were generated with earlier versions of the fabric analysers (Reference Wilson, Russell-Head and SimWilson and others, 2003).
Evolution of microstructure, c-axis fabric and grain statistics from samples 2-45 (a) and 2-46 (d) at small general strains and performed at –2ºC. (a) AVA image of sample 2-45 where the layering defined by the elongate grains was initially at 40º to x 1 and has undergone a rotation of 58 during which the elongate grains develop undulose extinction with minor migration recrystallization of the finer grains. Stereonets of measured c-axes are equal-area lower-hemisphere projections from the top and bottom outside areas of vertical sections and from a horizontal section through the centre of the deformed sample; the number, n, of measured c-axes is shown at the bottom of the stereonet. (b, e) Grain-size histogram with a normal distribution function (red) and mean value and standard deviation. (d) Experiment 2-46: the elongate grains were initially at 20º to x 1 and now lie at ~45º to x 1, and display undulose extinction with recrystallization initiated parallel to the grain anisotropy. Stereonets of measured c-axes are from top and bottom outside areas of vertical section and from a horizontal section through the centre of the deformed sample. (c, f) Manually digitized ice crystals (outlined by black rectangle) from sections from (a) and (d) respectively, and shape preferred-orientation analysis of the digitized grains shown in rose diagrams. Digitized crystals are labelled according to their sizes, in three categories for sample 2-45 (c) and two categories for 2-46 (f). For each category in (c) and (f) a shape preferred-orientation analysis has been performed. n is number of measured crystals. Dotted circles: percentage of grains for each direction; red line: mean direction; red dotted line: circular standard deviation.
AVA orientations and c-axis fabric in x 1 x 3 vertical sections. (a) The AVA colour code relates to the three-dimensional c-axis orientation. (b) Microstructure in –20ºC experiment 2-71 and stereonets show the volume distribution of c-axes in the outside of sample (c) and in central region (d). (e) Microstructure in –15ºC experiment 2-66 with stereonets (f, g) from the outside and central regions. (h) Microstructure in –10ºC experiment 2-64 with stereonets (i, j) from the outside and central regions. (k) Microstructure in –2ºC experiment 2-41 with stereonets (l, m) from the outside (zones 1 and 2) and central region (zone 3). (n, o) c-axes measured in a vertical section (x 1 x 3) (n) and a horizontal section (x 2 x 3) (o). (p) Microstructure in polycrystalline ice in –2ºC experiment 2-52 with (q, r, s, t) stereonets from the transitional regions (zones 1 and 2). (s) illustrates c-axes measured in the centre of the shear zone (zone 3) in a vertical section (x 1 x 3) and (t) in a horizontal section (x 2 x 3). All stereonets are equal-area lower-hemisphere projections; the number, n, of measured c-axes is shown at the bottom of the stereonet.
Mechanical Behaviour
Tests were carried out in confined vertical compression under a compressive stress of 0.22 MPa up to 17% shortening strain (Table 1). Bulk compressive strains (Table 1) were obtained by comparing the original versus final lengths of the sample in the x 1 direction. Bulk shear strains (Table 1) were calculated from the angular shear strain measured from the displacement of the upper piston relative to the base of the sample. Creep curves for representative tests are shown in Figures 2 and 3 as plots of octahedral shear strain rate as a function of octahedral shear strain, with the results of the shear tests of Reference Li and JackaLi and Jacka (1998) incorporated into Figure 2a.
Primary layering and hard-glide fabric orientations symmetrically related to the x 1 x 2 plane
At -2°C the secondary (minimum) octahedral shear strain rates observed in tests with initial layering parallel to the x 1 x 2 plane (Fig. 1a) range from 5.9 × 10-8s-1 to 2.4 × 10-7s-1 (Table 1). Compared with the polycrystalline ice deformed at -2°C by Reference Li and JackaLi and Jacka (1998), the creep curve (Fig. 2a) is similar, with the same accelerating stage in the creep curve (i.e. between the minimum and the tertiary strain rates). The only minor difference is the presence of an intermediate minimum at an octahedral shear strain of >10% that Reference Wilson, Russell-Head and SimWilson and others (2003) attributed to the initiation of shearing and the onset of dynamic recrystallization in the shear zone. Such an intermediate minimum is not observed in the uniaxial compression tests (Reference JackaJacka and Maccagnan, 1984; Reference Gao and JackaGao and Jacka, 1987). The position of the minima is also comparable with results presented by Reference Wilson, Russell-Head and SimWilson and others (2003) for more complex composite layered ice samples (Fig. 1c).
At -5°C where the layering was perpendicular (Fig. 1a) to the direction of shear, there was a rapid decrease in minimum octahedral shear strain rates (Fig. 2b;experiment 2-20) to 2.0 × 10-8s-1; whereas, where the ice layering was contained in the direction of shear (Fig. 1b), there was a gradual decrease in minimum octahedral shear strain rates to a minimum (Fig. 2b;experiment 2-16) of 1.4 × 10-8s-1.
At -10°C where the layering was perpendicular (Fig. 1a) to the direction of shear, there was a rapid decrease in minimum octahedral shear strain rates (Fig. 2c; experiment 2-59) to 4.0 × 10-8s-1. This is comparable to the creep curves for the composite layered samples (2-64 and 2-62; Fig. 2c) deformed at -10 and -8.2°C, which reached a minimum (at 3.2-4.0 × 10-8s-1) that is comparable to layered ice dominated by grains in a hard-glide orientation (experiment 2-59), with comparable tertiary strain-rate trends in samples deformed at higher temperatures.
With decreasing temperature, the minimum octahedral shear strain rates in the layered ice (Fig. 1a) vary from 6 × 10-8 s-1 (at -15°C) to ~4.5 × 10-9 s-1 (at -20°C; Fig. 2d). However, in the two -20°C experiments listed in Table 1, a true minimum was not reached. Similarly the bulk compressive (ε~ 0.5) and shear strains (γ~ 0.1) were also very small compared with the -15°C and higher-temperature experiments. In all experiments above -15°C, after 1% octahedral shear strain, the curves exhibit strain-rate increases to a greater tertiary strain rate and follow similar paths, with small decreases at octahedral shear strains of ~8% before undergoing a secondary increase.
Polycrystalline ice with random fabrics
Using the same stress configuration as the layered ice experiments, a set of tests on polycrystalline ice (mean grain area 1.6 mm2) was also undertaken (Fig. 3a). At -2°C the minimum octahedral shear strain rate varied from 2.1 to 2.4 × 10-8 s-1 and at -5°C it was 1.6 × 10-8 s-1. These results are at least an order of magnitude different from the results of Reference Li and JackaLi and Jacka (1998) and Reference Jacka, Li and HondohJacka and Li (2000) whose samples had a mean grain area of 1.2 mm2. Unfortunately these authors did not support their results with microstructural observations and failed to describe the localization of deformation in their samples. The results here also support the conclusions of other workers (e.g. Reference Duval and Le GacDuval and Le Gac, 1980; Reference Jones and ChewJones and Chew, 1983; Reference JackaJacka, 1984), who found that the creep rate was independent of crystal size.
In comparison, for the layered ice samples, all of which show the onset of zones containing recrystallized grains (Fig. 4) that transect the anisotropic ice, the minimum octahedral shear strain rates are very comparable. This suggests that after the onset of extensive recrystallization the minimum strain rate is not sensitive to variations in initial fabric intensity or grain morphology. However, the greatest difference in this comparison (Fig. 3) was that the minimum strain rate was reached at much smaller values of strain (~0.5%) in the polycrystalline ice than in the layered ice samples (~1%). This earlier onset of dynamic recrystallization in the polycrystalline ice suggests there may be an earlier attainment of tertiary creep if the ice mass is dominated by equiaxed polycrystalline grains.
Layering and dominant c-axis fabrics inclined to the x 1 x 2 plane
In order to evaluate the consequences of ‘easy’ versus ‘hard’ glide, another series of -2°C tests was undertaken using the anisotropic ice samples (Fig. 3b), where the orientation of layering in the latter tests was varied with respect to the compression axis and shear direction (Fig. 1c-e). Four tests (2-45, 2-47, 2-48 and 2-55; Table 1) were undertaken with an anisotropy initially oriented at 40° to evaluate the effect of easy glide on the basal plane as it was being rotated into a 45° orientation with respect to the compression axis. Minimum strain rates ranged from 1.6 to 5.9 × 10-8s-1 (Fig. 3b) and overlapped with the values obtained in the randomly oriented polycrystalline ice (Fig. 3a). As the dominant layering and hence the basal planes are initially inclined 20° and parallel to the x 1 x 2 plane, this minimum shear strain rate progressively increased to 3.2 and 5.9 × 10-8s-1, and hence overlapped with the values observed in the -5 and -10°C experiments (Fig. 2b and c).
Intuitively it would have been thought that grains close to the easy-glide orientation (~45°) with respect to the principal directions of compression and extension would be subjected to the largest shear stresses and should glide more easily. However, in experiments 2-48 and 2-55 where the layering was initially 40° the observed minimum octahedral shear strain rates (~1.6 × 10-8s-1) are lower and at greater octahedral shear strains (~1%) than in the experiments with the dominant grains lying close to the hard-glide orientation (e.g. 2-40;Fig. 3b). Similarly, in experiments 2-43 and 2-46 (Fig. 3b) where the layering was initially oriented 20° with respect to the direction of shear, the observed minimum octahedral shear strain rates (3.2-3.9 × 10-8 s-1) are higher and occur at lesser octahedral shear strains (~0.5%). The reason for this apparent anomaly can be seen in the evolution of the microstructure. With the progressive deformation, there was a passive rotation of the elongate grains into a ~45° orientation, but no further grain rotation of the c-axis (Fig. 5b). Accompanying this was the development of a new grain population initiated from the smaller randomly oriented grains as a zone that parallels the initial anisotropy (Fig. 5c). There is no evolution of a new grain population that transects the elongate grains until the octahedral shear strain exceeds 1% (Fig. 3b).
Microstructural and Textural Evolution
Role of grain anisotropy in recrystallization
The nature of recrystallization and the relationship to the shape anisotropy and crystallographic orientation is best illustrated in the microstructures observed in the tests on the anisotropic ice (Figs 4-6). Recrystallization has involved both modification of pre-existing grains and the preferential growth of new grains, rather than the rotation of grains with pre-existing shape anisotropy. The transition between these processes is clearly evident by observing the grain evolution with increasing temperature (Fig. 4) and in the lower strained samples (Fig. 5).
At -20°C, grain modification involves the preferential grain boundary migration of the randomly oriented small pre-existing grains (Fig. 4a; green and blue grains) and nucleation of very small (yellow) grains. Minimal change in the pre-existing grain geometry is reflected by their grain- size distribution and the shape-preferred orientation.
Between -15 and -10°C the anisotropic grains are progressively rotated (Fig. 4b and c) in the direction indicated by their shape-preferred orientation. There is a significant reduction in the area occupied by relics of the elongate grains with grain-size distribution that still reflects that observed at -20°C. Nevertheless, a slightly higher number of smaller grains indicates growth of the new nucleated grains and grain-size reduction of the larger pre- existing (blue) grains. At -15°C, new nucleated grains are arranged in more discrete planes parallel to the initial and rotated layers. With increasing temperature the grain texture becomes more granular, and grain boundary migration and grain growth become the dominant processes, causing a wider spread of the shape-preferred orientation.
By the time the ice has reached a tertiary creep regime at -2°C, there has been extensive grain boundary migration, an increase in grain size and a complete destruction of any pre-existing grain shape fabric (Fig. 4d). Initial stages of this process can be observed in Figure 5, where the elongate grains lie in an easy-glide orientation and have been subjected to comparatively small components of shear (γ = 0.05-0.07) and compressional strain (1.5-1.6%). All grains have experienced grain growth and there is undulose extinction in the relics of the pre-existing larger elongate grains (Fig. 5a). There is a preponderance of finer new grains (yellow grains in Fig. 5c) occurring on the boundaries of the pre-existing layering, and their mean direction parallels the compression direction. On the other hand, the area occupied by the older grains (0.2 to >0.14%; Fig. 5c) reflects the two initial grain shape distributions. Similarly the c-axis preferred orientation in both the vertical and horizontal sections (Fig. 5a) reflects the fabric distribution observed in the undeformed sample.
In experiment 2-46, there has been nucleation and growth of many more new grains, and the concentration can be explained by the localization of strain in the centre of the sample (Fig. 5d). There is elimination of the elongate grains that lie perpendicular to the layering (Fig. 5f), and the fraction of subgrain boundaries associated with undulose extinction also decreases. There is also the beginning of a new c-axis fabric pattern in both the vertical and horizontal sections (Fig. 5d).
We want to emphasize that (1) the initial layering of sample 2-45 (Fig. 4a; 40° clockwise) is in an easier glide orientation than sample 2-46 (Fig. 5d; 20° clockwise) and shows a larger minimal shear strain rate (Fig. 3b);and (2) shape-preferred orientation of the new nucleated grains (Fig. 5c and f; yellow grains) is always parallel to the compression direction and therefore at a higher angle to the pre-existing grain geometry and layering in sample 2-45, causing a microstructural interlocking of the initial layers. There is also minimal evidence of rotational recrystallization, as the aspect ratios of the larger preexisting grains are very large relative to the imposed strain. As a result, the shape modification cannot be caused by intracrystalline slip alone.
Textural evolution
At the termination of the -20°C tests, after applying the load for >50days, there are minimal changes in grain morphology from the start material, and bulk compressive and shear strains are very small (Table 1). In the elongate strongly anisotropic grains that lie in hard-glide orientations (i.e. c-axes lie perpendicular to the compression direction), there is only gentle warping (Fig. 6b). The small fibrous grains, however, that lie between these elongate grains in the outside region are identical to those of the start material. However, there is some minor grain boundary migration and growth in the central region of the sample after deformation. There is no evolution of any obvious zone of increasing shear strain, and the c-axis fabrics across the sample are comparable from the outside (Fig. 6c) to the centre (Fig. 6d), with a horizontal concentration that can be related to the initial elongate grains and a girdle pattern in the smaller grains. Therefore, in the -20°C tests, grain growth is insufficient to significantly modify the initial c-axis pattern, and grains are reasonably stable in this stress and temperature regime.
At the termination of the -15°C experiment (Fig. 6e) undertaken over the same time interval as the -20°C experiments, there was greater bulk shortening (4.9%) and minimal shear strain (γ = 0.09). Least deformed outside or transitional zones of the -15°C sample retain relics of the large elongate grains that display undulose extinction. The c-axes of these elongate grains reflect their initial orientation (Fig. 6f), which is perpendicular to the compression direction and hence parallel to the hard-glide orientation. Small fibrous grains are more equidimensional than the starting material and appear to have undergone grain growth, but still produce the scattered c-axis patterns seen in the stereonet (Fig. 6f). A central zone of recrystallized grains, which corresponds to a high-strain zone, now transects the sample, with a double-maximum concentration of c-axes (Fig. 6g) and the removal of the initial concentration of grains in the hard-glide orientation. It is inferred that this new grain population in the central zone has developed by misorientation and migration of boundaries as the deformation is progressively localized in the shear zone.
In the layered ice at -10°C, the bulk shortening of the sample is ~10%, with considerably higher shortenings (Table 1) in the composite samples (18%) for comparable bulk shear strains (γ~ 1). In the transitional zone, there was still a reflection of the original c-axis patterns (Fig. 6i). However, in the central zone of shearing in all -10°C and -8°C experiments there was a marked concentration of c-axes oriented oblique to the shortening direction (Fig. 6j) and sub-perpendicular to the shear zone, with sub-maxima lying at ~50° to the dominant maximum.
In the -2°C experiment illustrated here (Fig. 6k) and by Reference Wilson and SimWilson and Sim (2002), there is retention of the random c-axis pattern of the starting material in the outside transitional zones (Fig. 6l) and retention of the original grain size in relics of the elongate grains and extensive grain growth in the smaller grains (zone 1). The transition into a central region (Fig. 6m; zone 2), where extensive recrystallization equates with the zone of higher shear strain, is very rapid. In the higher-strained central area (zone 3), there are two asymmetric concentrations of c-axes (Fig. 6n) relative to the compression axis, with a strong single maximum rotated with the sense of shear. This pattern of preferred orientation has been reproduced in comparable experiments (Reference Wilson and SimWilson and Sim, 2002; Reference Wilson, Russell-Head and SimWilson and others, 2003). The c-axis pattern measured in horizontal sections, i.e. parallel to the shear zones, displays strong concentrations of c-axes that lie at high angles to the central shear zone (Fig. 6o). There are two asymmetric concentrations of c-axes relative to the compression axis with a strong single maximum, rotated in the direction of the sense of shear.
Polycrystalline ice (Fig. 6p) preserved in the outside zones (Fig. 6q; zone 1) shows the characteristic random c-axis orientation of the start material. This rapidly changes to a double-maximum fabric (Fig. 6r; zone 2) and an asymmetric double maximum in the central high-strain zone (Fig. 6s and t).
Discussion
As pointed out by earlier workers (e.g. Reference Barnes, Tabor and WalkerBarnes and others, 1971; Reference Budd and JackaBudd and Jacka, 1989), the temperature difference is the most important parameter in determining the position of the minimum creep rate, with very much higher values of shear stress required before the ice begins to deform at-20°C. However, it is clear from the current set of tests that at -20°C, where the minimum octahedral strain rate is ~10-9s-1, there is minimal evidence of recrystallization and change in crystallographic properties. At these low temperatures, there appears to be an absence of basal slip that would contribute to the generation of dislocations and the activation of recrystallization processes. At -15°C the minimum octahedral shear strain rate is 6 × 10-8s-1 and there is a marked onset of recrystallization and the development of anisotropic flow that corresponds to the development of a distinct two- maxima c-axis fabric pattern. For temperatures above -10°C, where the initial ice may have variable layer and crystallographic orientations, the minimum octahedral strain rate varies from 4 to 3.2 × 10-8 s-1. It is in this temperature regime (-10 to -2°C) that tertiary strain rates are comparable to those in natural ice where temperature-activated diffusion enhances the recrystallization processes (Reference Montagnat and DuvalMontagnat and Duval, 2000), with the ultimate development of strong c-axis maxima and sub-maxima that are asymmetric with respect to the normal to the shear zone boundary. This was particularly noticeable in the -2°C experiments where the central shear zone preserves no evidence of the pre-existing grain structures and is dominated by extensive recrystallization and grain growth.
Initiation of the grain growth observed involved two processes. (1) Growth of the small pre-existing randomly oriented grains during ongoing strain eventually leads to the formation of distinct new grains as described by Reference Duval, Ashby and AndermanDuval and others (1983). The pre-existing grain boundary morphology and grain shape appears to control the grain growth (Fig. 4). (2) Nucleation of small new grains at pre-existing layer boundaries is preferentially parallel to the compression direction (Fig. 5c and f). This is not considered to be static (normal) grain growth as described by Reference Alley, Perepezko and BentleyAlley and others (1986). Instead it is probably initiated by preferential growth by deformation-driven processes such as those described by Reference KambKamb (1959, Reference Kamb, Heard, Borg, Carter and Raleigh1972), which involve nucleating grains that minimize the difference between the elastic strain energy, the product of the mean stress and volumetric strain (see also Reference Paterson, Mayo, Kobayashi and WadsworthPaterson, 1990).
The c-axis patterns are similar to the observations of Reference Wilson, Russell-Head and SimWilson and others (2003) who suggested that the initial pattern of crystal-orientation fabrics and rate of rotation of the inherited layering govern the initial tertiary shear strain rate. However, from observing the differences in the creep curves (Figs 2 and 3) it is temperature that is critical in determining the initial minimum creep rate, and this masks the influence of the initial grain orientation or c-axis pattern. On the other hand, Reference Li and JackaLi and Jacka (1998) argue that at -2°C an initial small-circle fabric may enhance horizontal shear rates over an isotropic fabric. In considering these and the current set of experiments, it appears that the shear strain enhancement is intimately related to the degree of recrystallization and grain growth. Attainment of a tertiary creep rate of ~10-7s-1 to 10-6s-1 is independent of the initial crystal orientation, but is more dependent on the equilibrium crystal growth in the zone of higher shear strain and temperature (above -15°C).
Thus it appears that under a combined compression and shear, the onset of dynamic recrystallization is a gradual process that becomes more important above -15°C. This is supported by other laboratory investigations (Reference Kamb, Heard, Borg, Carter and RaleighKamb, 1972; Reference Bouchez and DuvalBouchez and Duval, 1982) and the observations in natural ice by Reference Montagnat and DuvalMontagnat and Duval (2000). At -15°C where there appears to be less compressive deformation (ε ~ 5%), a double-maximum fabric is produced. Above -10°C, where there is still a strong compressive component, doublemaximum fabrics develop (Reference WilsonWilson, 1982; Reference WilsonWilson and Russell-Head, 1982), with three- or four-maxima fabrics produced in the simple shear experiments of Reference DuvalDuval (1981). The coupling of these cumulative strains with temperature determines the final fabric evolution. This is clearly seen in the -2°C experiments, where a strong single-maximum fabric with a weaker subsidiary fabric is developed in both the anisotropic and polycrystalline ice samples.
The laboratory tests of Reference JackaJacka (1984) showed that secondary (minimum) creep rate and steady-state creep rate are independent of crystal size. In the tests by Reference Li, Jacka and BuddLi and others (1996) the polycrystalline ice had an initial mean crystal area of 1.2 mm2, whereas the polycrystalline ice in the current tests had a mean initial crystal area of 1.6 mm2. If the creep curves for polycrystalline ice at-2°C are compared, it is seen that the minimum strain rate for Reference Li, Jacka and BuddLi and others’ (1996) finer ice is higher (Fig. 2a, curves A-C) than for the coarser ice used in these experiments (Fig. 3a). However, as recrystallization and fabric development occurs, the steady-state creep rates converge between ~10-7s-1 and 10-6s-1. The current study also demonstrates that both temperature and initial variations in the orientation of both crystal fabric and anisotropy (Fig. 3b) appear to influence minimum strain rate.
The dominance of an initial fabric and the role of intracrystalline slip on the basal plane in controlling the flow can be observed in the anisotropic ice deformed at -2°C. Certainly grains close to the easy-glide orientation (~45°) with respect to the principal directions of compression and extension would be subjected to the largest shear stresses and should glide more easily. However, in the case of experiments 2-48 and 2-55 (Fig. 3b) the observed minimum octahedral shear strain rates (~1.6 × 10-8s-1) are lower than in the experiments with the dominant grains lying close to the hard-glide orientation (experiments 2-40, 2-43 and 2-46). As the deformation proceeds and a new grain population evolves, achieved by the migration of high-angle boundaries, there is a ‘loss of memory’ of initial fabric constraints, with the ice converging into a pattern of steady- state flow (Fig. 3b). In particular, the migration of the new grains in a high angle orientation into the pre-existing layering (Fig. 5c) causes a microstructural interlock, i.e. strain hardening, and therefore can lead to lower secondary (minimum) strain rates than expected when the initial layering is lying in an easy-glide orientation.
Migration of boundaries in the initial randomly oriented grains is the dominant feature at -20°C (experiments 2-68 and 2-71) and appears to be the main strain accommodation process during the stage of primary creep (Fig. 2d). Once a sample reaches the minimum (secondary creep), there appears to be nucleation of a discrete population of small new recrystallized grains (Fig. 4a; yellow grains). In all samples examined close to the minimum octahedral strain rate, there appears to be the nucleation of discrete new grains with a random c-axis preferred orientation (e.g. Fig. 5c; yellow grains). Above -15°C the distinction between a grain boundary migration and a discrete grain nucleation stage becomes more difficult (Fig. 4b and c). During the tertiary creep stage, there appears to be a further episode of grain boundary migration (Fig. 4d) that contributes to the final recrystallization history.
These results can be applied to naturally deformed ice where there may be a mesoscopic layering defined by seasonal stratification, or deformed ice with different grain aspect ratios and grain distributions. It has been shown that despite having initial strong fabrics these do not influence the evolution of the new fabric pattern where there is a general strain. This is particularly noticeable in grains with unfavourable orientations for intracrystalline slip; these grains appear to be removed by grain boundary migration of more favourably oriented grains and partially reoriented by selective recrystallization as described by Reference Wilson and ZhangWilson and Zhang (1996). A transition from a symmetrical double-maximum to an asymmetric maxima c-axis fabric pattern with increasing strain is not predicted by the fabric simulation models of Reference Van der Veen and WhillansVan der Veen and Whillans (1994), probably because the fabric simulation model does not involve combined crystallographic slip and dynamic recrystallization.
Conclusions
Results of this study indicate that the initiation of recrystallization, the strength and deformational behaviour of ice depend on the effects of temperature, which are linked to the localization of the deformation. At -20°C, there is practically no new grain growth when the deformation is close to a minimum creep rate of 10-9 s-1. For ice capable of sustaining temperatures higher than -15°C and differential stress over long periods of time, tertiary creep rates of 10-7 to 10-6s-1 are attained. This involves the ice deforming by intracrystalline gliding or through dynamic recrystallization. A pre-existing c-axis orientation has little effect once a steady-state creep rate above -15°C is attained. Hardening of the ice when the majority of the grains were in an easy- glide orientation only has implications for the short-term flow in an ice mass, as easy glide is a feature always assumed to facilitate increased strain. Similarly, the anisotropy in grain shape and orientation was found to only cause strain-rate fluctuations in the secondary (minimum) creep rate but have little effect on the steady-state creep rate, as recrystallization accompanied by the development of a new fabric controls the tertiary creep.
Moreover, the results of these tests suggest that finegrained polycrystalline ice may, as a result of more rapid recrystallization, approach steady-state tertiary creep more rapidly than coarser-grained ice and support earlier observations (e.g. Reference Duval and Le GacDuval and Le Gac, 1980; Reference JackaJacka, 1984). Both the polycrystalline and layered ice exhibited approximately the same ductility in the tertiary creep regime. At lower temperatures (e.g. -20°C and -15°C) the duration of the primary stage of creep was much longer, with minima at strains of ~0.5%. The position of the minima at -10 to -2°C at ~1% strain is also comparable with results presented by Reference Wilson, Russell-Head and SimWilson and others (2003) for more complex composite layered ice samples, which suggests any internal layering and a previous c-axis orientation pattern in an ice mass has minimal effect on the overall flow pattern (tertiary creep) as all longer-term tests reached tertiary strain rates of 10-7 to 10-6s-1. This is comparable to many bodies of polar ice which deform in a near-steady-state tertiary creep regime (Reference PatersonPaterson, 1991).
Acknowledgements
This work was supported by Australian Research Council grants (DP0344708 and DP0773097). D. Russell-Head is thanked for technical advice and assistance regarding the fabric analyser. Li Jun and other members of the Australian Antarctic Division Glaciology Program are thanked for discussions about the deformation rig and for supplying the software used to record the creep curves. This paper presents work initially undertaken while the authors were at the University of Melbourne. We thank that institution for its support, and H. Sim for help with the laboratory work and some of the initial data analysis. Two anonymous reviewers and T.H. Jacka are thanked for comments on an earlier version of the paper.
