A. Ice Deforms Primarily by Dislocation Glide on the Basal Plane {4}

Many mechanisms can contribute to deformation, and there almost certainly are situations in which ice enters other deformation fields. If grain-size were small enough, diffusional creep (Nabarro-Herring or Coble creep) would become important (Reference Duval, Ashby and AndermanDuval and others, 1983). In low-density firn, grain-boundary sliding may be dominant (Reference AlleyAlley, 1987). However, for conditions thought to occur in densified, ice-sheet ice, direct laboratory observations and theory strongly support dislocation glide on the basal plane (Reference Weertman, Whalley, Jones and GoldWeertman, 1973, Reference Weertman1983; Reference Duval, Ashby and AndermanDuval and others, 1983; Reference Fukuda, Hondoh and HigashiFukuda and others, 1987). Other glide systems are 100 to 1000 times “harder” than basal glide.

The basal plane has only two independent slip systems, but continuum deformation of a polycrystal requires four or five independent slip systems. The basal glide must be accommodated in some way: by hard, prismatic glide; by dislocation climb; by grain-boundary sliding; by Nabarro-Herring type diffusion; or in some other manner. Diffusional processes (grain-boundary migration, dislocation climb, Nabarro-Herring creep) seem likely considering the low strain rates and relatively high temperatures (almost always > 70% of the absolute melting temperature, and often >90%) of most natural ice deformation.

B. If Strain Rate, Temperature, Stress and Cumulative Strain are Large Enough, Strain Rate Increases Approximately with the Third Power of the Stress {4}

The (or a?) flow law for ice usually is written as

(1)

Here is the strain rate, A_o is a prefactor that may depend on such factors as fabric, grain-size or impurity concentrations (and can be written as aA, where a is an enhancement factor on the “normal” behavior A, and may be a tensor quantity), Q is the activation energy, k is Boltzmann’s constant, T is the absolute temperature, σ is the stress deviator and n is the power-law constant. For uniaxial deformation, and σ refer to that deformation; more generally, they are the second invariants of the strain-rate and deviatoric-stress tensors. A weak dependence on pressure often is ignored, or is allowed for by referring T to the in situ melting point (Reference Weertman and WeertmanWeertman, 1968; Reference Doake and WolffDoake and Wolff, 1985).

For dislocation glide, n should be equal to 3 based on the following analysis, summarized from the review by Reference Weertman, Whalley, Jones and GoldWeertman (1973). The strain rate is simply the product of the rate at which dislocations move through a material, their abundance and the deformation achieved per dislocation, or

(2)

where α is a geometric factor of the order of 1, ρ is the dislocation density, b is the Burgers vector of the dislocations (approximately the molecular spacing) and ν is the dislocation velocity. The dislocation velocity is, both theoretically and observationally,

(3)

where γ is a constant and σ is the applied stress causing dislocation glide. The dislocation density, ρ, for steady-state deformation is approximately that which gives an internal stress from the dislocations equal to the applied stress. The internal stress is μb/(βd), where β is a constant of the order of 1, μ is the shear modulus, and d is the dislocation spacing and is related to the dislocation density by ρ ≈ 1/d ². Then

(4)

Combining Equations (2)–(4) yields

(5)

which is essentially Equation (1) with a more physical statement of A₀ and with n = 3.

This result has broad support from laboratory and field observations, especially at high stresses (σ > ≈ 1 bar) (Reference Weertman, Whalley, Jones and GoldWeertman, 1973; Reference HobbsHobbs, 1974; Reference Budd and JackaBudd and Jacka, 1989); results at lower stresses (σ < ≈ 1 bar) are equivocal as to whether n = 3 or n < 3 (including n = 1), probably because of the great difficulty in conducting such tests (see Reference Weertman, Whalley, Jones and GoldWeertman, 1973; Reference Budd and JackaBudd and Jacka, 1989).

C. If Temperature, Stress and Cumulative Strain are Small Enough but Dislocation Glide Still Dominates, the Strain Rate Increases Approximately with the First Power Of the Stress {2}

This is the regime of Harper-Dorn creep. From the equations above, it is evident that, if dislocation density is constant, then linear-viscous behavior is predicted. Constant (stress-independent) dislocation density may apply at very small strains before dislocation multiplication has begun, or at very small stresses and strain rates, which may not be large enough to cause dislocation multiplication above the initial annealed concentration (see Energetics, below) or which may allow recovery (principally diffusional) processes to “heal” dislocations as fast as they are produced.

There is little doubt that creep exists with n = 1. There is greater uncertainty whether n = 1 creep is steady, and whether it occurs in ice under natural conditions (see Reference WeertmanWeertman, 1985). A number of workers have reported n = 1 for natural ice (or at least n < 3) (see references in Reference Weertman, Whalley, Jones and GoldWeertman, 1973; also Reference Doake and WolffDoake and Wolff, 1985; Reference Pimienta and DuvalPimienta and Duval, 1987), but this point is by no means settled (Reference PatersonPaterson, 1985; Reference Budd and JackaBudd and Jacka, 1989).

A. The Average Grain-Size in an Existing Population of Grains Always Increases (Grain Growth Always Occurs), Although Creation of New, Smaller Grains may Cause the Average Grain-Size to Decrease or Remain Constant Over Time {4}

Grain growth in ice is driven primarily by two factors: the curvature of energetic grain boundaries, and differences in stored strain energy between grains. The curvature of grain boundaries that possess a surface tension causes what amounts to a pressure difference between adjacent grains, with the higher pressure on the concave side of the boundary (usually the smaller-grain side). Molecules tend to jump from the high-pressure side to the low-pressure side, causing the boundary to move toward the high-pressure side. The process usually accelerates with time; as the boundary moves toward the center of a small grain, that grain becomes smaller, the curvature of its boundary increases, and the migration speeds up until the smaller grain is consumed entirely. The average grain-size in the sample is increased every time one grain consumes another in this fashion.

Boundary migration can be stopped or slowed by second-phase particles, pores or dissolved impurities at high concentrations; however, concentrations found in natural glacier ice rarely or never are high enough to stop grain growth (see review by Reference Alley, Perepezko and BentleyAlley and others, 1986a,Reference Alley, Perepezko and Bentleyb). Grain growth also can be slowed significantly if boundaries assume special orientations that cause them to have low energy (see Grain Growth B, below), but this also does not stop grain growth. The tendency for boundary migration from curvature and tension exists and will always tend to cause grain growth until arrested by second-phase particles and other factors, which in typical ice-sheet ice would require average grains about 1 m across or larger (Reference Alley, Perepezko and BentleyAlley and others, 1986b). If grain growth in a polycrystal with ordinary, high-angle grain boundaries is driven by the curvature and surface tension of boundaries and if the effects of contaminant materials are not large, then the cross-sectional area of grains increases linearly with time; this is the regime of normal grain growth (e.g. Reference GowGow, 1969; Reference Alley, Perepezko and BentleyAlley and others, 1986a).

Boundary migration also can be caused by differences in stored strain energy across grain boundaries. Dislocations and other defects introduced into a crystal lattice by strain increase its energy. Because boundary migration reconstitutes the ice through which the boundary moves, the ice just behind a moving boundary has not been deformed and thus is unstrained and has low energy. (For very rapid migration, excess vacancies may enter the lattice behind a moving grain boundary, raising its energy and slowing the boundary (Reference GleiterGleiter, 1979), but this is unlikely to be a significant factor under typical conditions in ice sheets (Reference Alley, Perepezko and BentleyAlley and others, 1986a).) When a boundary separates strained from unstrained ice, the free energy of the strained side is higher, and the total energy can be lowered if the boundary migrates into the strained grain. Once started, the migration of a boundary in deforming ice will tend to continue because the ice immediately behind it is unstrained and that in front is older and thus more strained. Eventually, a strained grain being consumed in this fashion will become small enough that curvature effects also are important, and it will be consumed rapidly. (Initially, sufficiently large differences in stored strain energy can cause a boundary to move away from its center of curvature, away from a small unstrained grain into a large strained grain, in the opposite direction of that expected from curvature alone. As the strained grain becomes small, both strain energy and curvature work in the same direction, causing the strained grain to be consumed. Consumption of any grain always increases the average grain-size.)

It thus appears that if grain boundaries are not pinned by second-phase particles, they always migrate. And, if new grains are not created, migration of boundaries inevitably leads to grain growth. The only ways to avoid this are if each grain consumes a neighbor at the exact same rate that another neighbor consumes it, or if all grains have the same strain energy and are perfectly the same size and shape and perfectly packed, so that all boundaries are exactly straight. Real materials always are imperfect and have some randomness, so grain growth always occurs.

B. If the Average Grain-Size is Not Increasing (or is Not Increasing as Rapidly as Expected for Normal Grain Growth), Then Existing Grains are Being Subdivided (Polygonization) or Entirely New Grains are Being Created (Recrystallization) {4}

Consumption of some grains by others without increase in the average grain-size in a constant-density material requires that new grains are formed to replace those lost. Two mechanisms are believed to dominate such new-grain formation: polygonization and recrystallization (“rotation recrystallization” and “migration recrystallization”, respectively, in the terminology of Reference PoirierPoirier (1985, chapter 6)).

Polygonization is the formation of new grains by subdividing the old grains. If different parts of a grain are subjected to different stress states (virtually guaranteed by grain–grain interactions in a polycrystal), then the grain can become bent or twisted. Dislocations tend to organize between relatively undeformed regions called sub-grains to form sub-grain boundaries that relieve this bending or twisting and lower the energy of the system. If a sub-grain boundary becomes sufficiently strong (or a sub-grain becomes sufficiently rotated), then the boundary becomes a full grain boundary; the division between sub-grains and grains is somewhat arbitrary. Low-angle grain boundaries produced by polygonization initially cause a small mismatch between grain lattices, which means that the boundaries have low energy and migrate slowly under curvature-induced stress. Continued deformation can increase the angle of mismatch. Polygonization often can be identified by occurrence of more nearest-neighbor grains with small misorientations than would occur randomly in a sample. (Poirier (1985) distinguished between the formation of new grains by subdividing, which he called “creep polygonization”, and the refinement of this structure by further deformation, which he called “rotation recrystallization”. Here I use “polygonization” to refer to both.)

Recrystallization (“migration recrystallization” of Reference PoirierPoirier (1985)) is the production of new grains at initially high angles to their neighbors. The new grains typically nucleate from existing grain or sub-grain boundaries separating regions with different stored strain energy, with the boundary bowing out toward the more strained side. The stored strain energy in the region ahead of the migrating boundary must be large enough to overcome the increase in grain-boundary energy caused by this bowing out of the boundary (see Energetics, below). The new grains in recrystallization usually form in energetically favored orientations with their c-axes at high angles (≈45° ) to compressional axes (e.g. Reference Kamb, Heard, Borg and RaleighKamb, 1972).

C. Growth Selection Does Not Form Fabrics Over Broad Ranges of Normal Grain Growth {3}

The energy of grain boundaries is nearly independent of relative grain orientations, except for certain favored orientations (e.g. twin boundaries, small misorientations) which may reduce the boundary energy by more than an order of magnitude (see review in Reference VerhoevenVerhoeven, 1975, chapter 7). This opens the possibility that, if two grains with similar orientations grow together, one will not consume the other rapidly because of the low boundary energy. Repetition of such growing together of similarly oriented grains could yield regions in which grains exhibit a preferred orientation. However, simple calculations show that the favored orientations are sufficiently rare that this is unlikely to matter in natural ice, and indeed, no strong fabrics are observed during normal grain growth.

Growth selection may contribute to fabrics when stored strain energy is higher than in normal grain growth. For example, polygonization produces low-angle boundaries, and the grains across such boundaries are consumed only slowly. However, the new grains will occur within an area that started as a single grain, and boundaries between such areas are high-angle grain boundaries that can grow. In this way, relatively large areas of grains may develop with similar orientations. Also, it is possible that grains in some orientations will store more strain energy than those in other orientations. The high-energy grains would be consumed in grain growth, leaving grains with a preferred orientation. In the first case, any fabric measurement that averaged over a large enough area would be unaffected by growth selection. The second case might matter, but is essentially unexplored in ice.

A. During Ice Deformation, c-axes Always Rotate Toward Compressional Axes and Away From Tensional Axes {4}

Deformation parallel to the basal plane in ice is 100–1000 times easier than in any other direction (Reference Weertman, Whalley, Jones and GoldWeertman, 1973). An applied stress thus causes basal glide for almost any orientation, which tends to cause grain deformation. However, in a polycrystal, the probability that adjacent grains will deform in exactly the same direction is vanishingly small (except in the limit of large strains in simple shear) (Fig. 1). Even if the deformation is in the same direction in adjacent grains, the active slip planes are unlikely to match exactly across the grain boundary. Adjacent grains then interfere with the deformation of any given grain. This interference can be avoided if the grain rotates at the same time that it deforms in basal glide. The adjustment between two grains is taken up by slip along the grain boundary, which allows sliding easily (Reference Ignat and FrostIgnat and Frost, 1987). The result is grain deformation on the basal plane, grain elongation, and rotation of the normal to the glide plane (the c-axis in ice) toward compressional and away from tensional axes (Figs 1 and 2). This result is general for essentially any crystal deforming by dislocation glide and constrained by a testing apparatus or adjacent grains in a polycrystal (e.g. Reference VerhoevenVerhoeven, 1975, chapter 3).

Fig. 1.

Cartoon showing grain rotation from deformation of a polycrystal. In (a), uniaxial tension σ_x causes resolved shear stresses on the basal planes, ba, of two grains, which cause mutual interference for basal glide. Deformation is achieved by basal glide plus rotation of c-axes, c, away from the tensional axis, leading to the situation in (b).

Fig. 2.

Schematic evolution of c-axis fabrics as a function of stress state and active processes. An initial random distribution (top) yields solid-cone distributions (left) if rotation alone is active, and girdle/surface-cone distributions (right) if recrystallization also is active.

B. Rotation of Existing Grains in Polycrystals is Described to Good Approximation by the Results of Reference Azuma and HigashiAzuma and Higasbi (1985) {2}

Reference Azuma and HigashiAzuma and Higashi (1985) deformed thin ice specimens uniaxially and found that the strain in each grain was approximately proportional to the resolved shear stress along its basal plane, with the bulk strain equal to the strain of individual grains averaged over the entire sample. They then found that this model predicts the observed fabric development at Dye 3, Greenland, in the region of normal grain growth, where such a model for an existing population of grains should apply. More recently, Reference Lipenkov, Barkov, Duval and PimientaLipenkov and others (1989) have found that the Reference Azuma and HigashiAzuma and Higashi (1985) model, modified for uniaxial tension following Reference Fujita, Nakawo, Mae and MatsudaFujita and others (1987), accurately describes fabric observations in the Vostok, East Antarctica, core.

Briefly, under uniaxial tension or compression, the resolved shear stress on the basal plane of an ice crystal is proportional to the geometric factor S = cos ϕ₀ sin ϕ₀, where ϕ₀ is the angle between the c-axis of a grain and the unique stress axis. A grain that elongates from length l_o to length l has incremental strain Δϵ _g = (l – l _o)/l _o and its c-axis rotates from angle ϕ_o to angle ϕ according to cosϕ = cosϕ_o l _o/l for uniaxial tension and cosϕ = cosϕ _o l/l_o for uniaxial compression. An incremental strain Δε averaged over the entire sample causes strain in an individual grain Δϵ _g = ΔϵS/ S where S is the average of the geometric factor S over all grains in the sample. Combining these relations yields

(6a)

for rotation of a grain in uniaxial tension, and

(6b)

for rotation in uniaxial compression; pure shear is the superposition of these, and simple shear is pure shear plus a rigid-body rotation. The parameter S is the geometric softness of an ice sample to deformation, arising from the average resolved shear stress on basal planes. Numerical integration of Equations (6a) or (6b) for uniaxial deformation shows that S generally decreases with increasing strain, but that strains in excess of 100% are needed to halve S (Reference AlleyAlley, 1988); we will find later that uniaxial deformation without recrystallization usually does not persist to strains significantly greater than 100%, so that the fabric effects from uniaxial deformation are not large, but that simple shear does reach larger strains and affect ice hardness significantly through fabric development.

The main principle underlying Equations (6) is that the strain in any grain is proportional to the resolved shear stress on its basal plane; the rest of the derivation is simply geometry. This is not a necessary result; in some materials at some temperatures, the strain in individual grains may be nearly independent of orientation, or the strain may be such as to minimize differences in deformation between adjacent grains, or some intermediate behavior (e.g. Reference EtchecoparEtchecopar, 1977; Reference Lister, Paterson and HobbsLister and others, 1978). Strain proportional to resolved shear stress would be most likely at relatively high temperatures and low stresses, as in ice. Under these conditions, diffusional processes can accommodate differences in strain between adjacent grains (Reference Duval and LliboutryDuval and Lliboutry, 1985), allowing grain elongations to be nearly independent of neighbors, although still with constraint from neighbors sufficient to force rotation. This adds some theoretical support to the laboratory and field results that Equations (6) do describe grain rotation in ice. However, further data on this point would be desirable.

C. In Polygonization, New Grains Form within A Few Degrees of Old Grains, So That Fabrics are Not Affected Greatly {4}

This is almost true by definition. If a new grain forms by progressive sub-grain rotation, and is defined to have formed when it has reached an angle of a few degrees, then the new grain will plot very close to the grain from which it formed on a c-axis fabric diagram, causing little change in the overall fabric. Polygonization does not remove the old grains, and the boundary between a new grain and its parent may have low enough energy to slow consumption of one by the other.

D. New Grains in Recrystallization Form with c-axes Midway Between Compressive and Tensile Axes, Causing Formation of New Grains to Change Fabrics {3}

During recrystallization, new grains are believed to form in energetically favorable postions with c-axes approximately bisecting the angle between compressional and tcnsional axes, so that resolved shear stress on basal planes is maximized (see Reference Kamb, Heard, Borg and RaleighKamb, 1972). This produces a girdle about the unique axis in compression and tension, and two concentrations in the plane containing the compressional and tensional axes in pure shear and simple shear (and usually with greater scatter of positions perpendicular to this plane than parallel to it. For simple shear, this means that one of the concentrations is normal to the shear plane, and the other is almost, although generally not exactly, in the shear plane) (Fig. 2).

The new grains are relatively strain-free and consume older grains by boundary migration. The new grains also accumulate strain energy while rotating and deforming, and eventually are consumed by newer grains. Recrystallization thus tends to replace an existing population of grains with a new population of grains which have c-axes that on average are rotated a little toward compressional axes from nucleation positions midway between tensional and compressional axes.

A. If Normal Grain Growth is Occurring (That is, If the Cross-Sectional Area of Grains is Increasing Linearly with Time, the Grain-Size Distribution Is Unimodal and Grains are Equant and Convex, with No Strain Shadows at Least Initially) Then:

B. If c-axis Fabrics are Tightening, Strain Shadows and Sub-Grains Appear, and Grain Growth is Slower Than Normal (Including Zero or Negative Growth), Then:

C. If a Strong Single-Maximum Fabric Develops in Simple Shear, Together with Slow, Zero or Negative Grain Growth and With Strain Shadows or Sub-Grains, Then:

D. If a Girdle Fabric Develops During Active Deformation, Together With Zero Grain Growth and With Strain Shadows, Then:

Normal Grain-Growth Region — Hypothesis: n = 1, but 1 ≤ n ≤ 3 Possible. Fabric Enhancement from Reference AlleyAlley (1988)

In the upper 100–1000 m of the central regions of cold ice sheets, stored strain energy is small compared to grain-boundary energy. Deformation occurs by dislocation glide on basal planes and probably is accommodated by diffusional processes, c-axes of grains rotate toward compressional axes during deformation. No new grains are produced. The c-axis fabric then records the cumulative strain (by its tightness about the compressional axis) and the stress state (by its shape and orientation), and affects the ice hardness (through the fabric effect on the prefactor A₀ in Equation (1)), probably as sketched by Reference AlleyAlley (1988).

The best estimate of the dislocation density needed for steady-state creep would not produce enough strain energy to affect normal grain growth; hence, we lack a direct probe to distinguish between the steady-state dislocation density (which gives n = 3 in Equation (1)) and the case where diffusional recovery processes or other processes balance or prevent dislocation multiplication and yield n = 1. It is noteworthy, however, that recovery processes must be active or dislocation multiplication must be suppressed, because the normal grain-growth regime can persist to ≈ 10–100% strain or possibly more (Reference Gow and WilliamsonGow and Williamson, 1976; Reference Herron, Langway, Brugger, Langway, Oeschger and DansgaardHerron and others, 1985; Reference Lipenkov, Barkov, Duval and PimientaLipenkov and others, 1989), whereas at higher stresses or temperatures in the laboratory ≈1% strain initiates recrystallization, reaching a steady state of deformation by ≈ 10% strain (Reference JackaJacka, 1984b). Recovery processes remove dislocations and reduce dislocation density. This, combined with reasonably good evidence for n = 1 at low stresses in laboratory experiments on other materials (see Reference WeertmanWeertman, 1985) and a suggestion of n = 1 from laboratory and field observations on ice (Reference Weertman, Whalley, Jones and GoldWeertman, 1973; Reference Doake and WolffDoake and Wolff, 1985; Reference Pimienta and DuvalPimienta and Duval, 1987), make n = 1 the leading candidate for the flow law during normal grain growth. However, 1 ≤ n ≤ 3 is possible, including the possibility that n lacks a physical meaning, but only is a curve-fitting parameter for creep controlled by a changing balance of processes.

Polygonization Region – Hypothesis: n = 1, but 1 ≤ n ≤ 3 Possible. Fabric Enhancement from Reference AlleyAlley (1988)

Deeper than the normal grain-growth regime, but still in approximately the upper one-half to two-thirds of the central regions of large ice sheets, grain-size from at least some sites (Reference Gow and WilliamsonGow and Williamson, 1976; Reference Herron, Langway, Brugger, Langway, Oeschger and DansgaardHerron and others, 1985) is nearly independent of depth and age, abundant sub-grains and strain shadows occur, but fabrics are similar to and slightly tighter than those in the normal grain-growth region. This behavior probably arises from polygonization. That is, dislocations organize into sub-grain boundaries to reduce their energy, and with increasing deformation these sub-grain boundaries strengthen until they cross an arbitrary threshold to become grain boundaries. The new grains have similar orientations to old grains, so fabrics can continue to tighten with deformation. Lack of girdle fabrics shows that old grains are not being consumed by new, high-angle nuclei. The easiest explanation is that recrystallization is slow or absent. Cumulative strains can approach or slightly exceed 100% in this region without causing rapid recrystallization, whereas 1–10% strain causes recrystallization during more rapid deformation in the laboratory (Reference JackaJacka, 1984b; Reference Jacka, Gao Xiangqun. and Gao Kun,Jacka and Gao, 1989), showing that recovery processes continue to be important; however, the sub-grains and polygonization attest to increasing stored strain energy compared to normal grain growth. Because recovery processes are active, n = 1 remains possible, although the increasing stored strain energy makes n = 3 more likely than for normal grain growth.

Laboratory experiments (e.g. Reference JackaJacka, 1984b) show that uniaxial deformation causes ice to harden by a factor of roughly 2–10 with increasing strain up to about 1%, prior to initiation of recrystallization, with this hardening probably caused by development of dislocation tangles (Fig. 3). If n = 1, then such tangles should not be significant in the interiors of sub-grains, and the ice hardness should be unchanged from the normal grain-growth regime, except for fabric changes that can be treated in the same way as in that regime. If n > 1, then it is likely that dislocations are building up within the sub-grains, and they may become tangled and harden the ice. Then it is possible that A₀ should be decreased somewhat independent of fabric.

Simple-Shear Region – Hypothesis: n = 3, but 1 ≤ n ≤ 3 Possible. Fabric Enhancement from Reference Budd and JackaBudd and Jacka (1989, fig. 9)

Still deeper, the basal shear stress comes to dominate over the normal stresses, so simple shear occurs. The fabric is dominated by a strong vertical single maximum, with grain-size small and relatively constant. A critical question is whether a second, weaker maximum exists.

As discussed above, laboratory experiments (e.g. Reference Kamb, Heard, Borg and RaleighKamb, 1972) show that recrystallization in simple shear produces nuclei in two positions: one strong maximum with c-axes normal to the shear plane, and one weaker maximum with c-axes nearly parallel to the shear direction and nearly parallel to the shear plane, c-axes in the stronger maximum will tend to stay there and c-axes in the weaker maximum will tend to rotate toward the stronger one.

If many nuclei are being formed near the horizontal position, they should be observable as a second maximum, even if weak and spread out through rotation. However, if polygonization is occurring without recrystallization, then only one maximum should exist. As discussed below, occurrence of two maxima argues for η = 3, whereas one maximum leaves open the question of η between 1 and 3, inclusive.

Available data are equivocal. Some observations show no clear evidence of a second maximum (e.g. 1302 m in the Byrd core; Reference Gow and WilliamsonGow and Williamson, 1976), but others do seem to show a weak second maximum (e.g. 1210, 1487 and 1634 m at Byrd). There seems to be no way to achieve such a secondary maximum without recrystallization. Then either recrystallization occurs in some ice undergoing simple shear deep in ice sheets but not in adjacent ice (possibly from differences in ice or from waves of recrystallization), or else recrystallization occurs in all such simple-shear ice but varies in strength relative to rotation. (Or, the observing or contouring conventions chosen by field observers failed to highlight some weak secondary maxima that did exist.)

It thus appears that at least a little recrystallization occurs in deep simple-shear zones of ice sheets, probably along with much polygonization, favoring n = 3 (see below). This would agree with analyses of ice-sheet shape and velocity that yield n = 3, because the shear zone near the bed largely controls the ice-sheet shape and velocity (Reference Hamley, Smith and YoungHamley and others, 1985). However, we cannot absolutely rule out n < 3. The best data on fabric-enhancement factors in simple shear are probably those summarized by Reference Budd and JackaBudd and Jacka (1989, fig. 9), so they should be used initially.

Warm, Basal Multiple Fabrics – Hypothesis: No Good Hypothesis

It is not clear what flow law to use for the very coarse, interlocking grains with a multiple-pole, girdle-type fabric found in warm ice near the bed at Byrd Station (Reference Gow and WilliamsonGow and Williamson, 1976) and elsewhere. The very large grains argue against rapid deformation. The annealing hypothesis of Reference Budd and JackaBudd and Jacka (1989; and earlier work) seems to provide the best available explanation. If correct, then this texture and fabric indicate that the ice was deformed rapidly but now is nearly stagnant and at high temperature. When deformation stops, recrystallization in twin orientations is initiated and then “runs away”, with grains growing very large. The suggestion of stagnant, or nearly so, ice indicates that there is no need for a flow law for this ice. If a flow law is needed, the girdle fabric shows that recrystallization has occurred, but the very large grains show that recovery processes have been active since the initiation of recrystallization. A recovery-dominated n = 1 might be slightly more likely, but the evidence is weak.

Ice Shelves and Ice Streams – Hypothesis: n = 3. Fabric Effects From Jacka and Maccagnan (1984)

With ice shelves, and probably with their inland extensions as ice streams, we finally reach a firmer basis for the flow law. Laboratory experiments show clearly that relatively high normal strain rates (high stresses or temperatures) cause girdle-type fabrics (i.e. no c-axes near compressional axes), a true steady state of deformation (a constant applied stress yields a constant strain rate, grain-size and c-axis fabric, and thus almost certainly a steady-state average dislocation density) and n = 3 (e.g. Reference JackaJacka, 1984a, Reference Jackab; Reference Jacka and MaccagnanJacka and Maccagnan, 1984). Ice shelves under temperatures and normal stresses approaching laboratory values develop girdle-type fabrics and are observed to spread with n = 3 (e.g. Reference ThomasThomas, 1973; Reference Jezek, Alley and ThomasJezek and others, 1985). Reference Doake and WolffDoake and Wolff (1985) argued that ice-shelf spreading is consistent with n = 1 or n = 3, but required that other workers significantly underestimated their errors and that the good, high-stress data from Amery Ice Shelf, Antarctica, be ignored to avoid the natural conclusion from the data that n = 3. A girdle-type fabric has been observed seismically on Ice Stream Β in West Antarctica (Reference Blankenship, Alley and BentleyBlankenship and others, 1989).

The girdle-type fabrics of the ice shelves and ice streams show that recrystallization occurs. Because this requires high stored strain energy in the grain or sub-grain into which a new boundary migrates, rather than just in sub-grain boundaries, this shows that grains contain a large dislocation density at least equal to the steady-state value, and thus argues for n = 3. Then n = 3 is a good hypothesis wherever recrystallization occurs. (At higher stresses than occur in natural ice sheets, n > 3 is possible.)

Reference Jacka and MaccagnanJacka and Maccagnan (1984) showed how the angle of the girdle fabric varies between onset of recrystallization and approach to steady state, and how the ice hardness depends on this angle. The values recommended by Jacka and Maccagnan provide a good starting point for modeling.

Summary of Hypotheses

We thus sketch a crude picture of macroscopic and microscopic properties in cold natural ice. Rapid deformation under normal stresses in ice shelves and ice streams causes recrystallization and n = 3. Simple-shear deformation in deep ice probably causes recrystallization and n = 3, but the evidence is not as strong as for ice shelves and does not rule out n = 1 or some intermediate value. Polygonization probably dominates in somewhat shallower ice with normal stresses, and recovery processes must be important. The recovery processes suggest n < 3, including n = 1, but the evidence is weak. In very shallow ice with normal grain growth, recovery processes are most rapid relative to storage of strain energy, dislocation multiplication may be suppressed, and n = 1 (or at least n < 3) is most likely, but proof again is lacking. A reasonable approach to testing the flow law would be to divide the ice sheet into layers based on macroscopic observations of fabric and texture, and to assume n has a constant value (or a constant slope of variation) within a layer.

The onset or re-orientation of recrystallization during progressive burial of ice or initiation of streaming flow require that the ice first harden from storage of strain energy and then soften from recrystallization removal of this strain energy, as shown in Figure 3. This should cause recrystallizating zones to be surrounded by bands of stiff ice, which may have important implications for ice flow (which I (Reference AlleyAlley, 1990) have discussed briefly and hope to explore in a future manuscript).

Progressive strain causes fabric formation, c-axes rotate toward compressive axes in normal grain growth and polygonization, in a process that cannot obtain a true steady state under normal stresses, c-axes rotate to the normal to the shear plane in simple shear, and form a girdle around compressive axes in recrystallization, and can reach a steady state in these cases. The fabric effect on creep rates is a tensor property but reduces to a simple enhancement factor for deformation in one stress state. Laboratory experiments constrain the enhancement factor well for simple shear and recrystallization, and theory constrains the smaller variation in enhancement factor under normal grain growth and polygonization.