Introduction
The effects of dispersed particles on the creep behavior of polycrystalline ice have been investigated by a number of workers (Reference Holdsworth and BullHoldsworth and Bull, 1970; Reference Nayar, Lenel and AnsellNayar and others, 1971; Reference Baker and GerberichBaker and Gerberich, 1979; Reference Lange and AhrensLange and Ahrens, 1983; Reference Shoji and LangwayShoji and Langway, 1985; Reference Durham, Kirby and SternDurham and others, 1992; Reference Song, Cole and BakerSong and others, 2004). These studies include ice containing ultrafine (15 nm diameter) amorphous silica particles (Reference Nayar, Lenel and AnsellNayar and others, 1971; Reference Lange and AhrensLange and Ahrens, 1983), which show similar mechanical behavior to particle-strengthened metals, i.e. at low concentrations the strength increases as the particle content increases. In contrast, when the particle size is similar to that of silt (several hundred microns), the particles can decrease the strength (Reference Holdsworth and BullHoldsworth and Bull, 1970; Reference Baker and GerberichBaker and Gerberich, 1979; Reference Shoji and LangwayShoji and Langway, 1985). For example, Reference Baker and GerberichBaker and Gerberich (1979) found that at temperatures from −20°C to −5°C, solid inclusions in concentrations from 1.3 to 6.6 vol.%, increased the creep strain rate; Reference Shoji and LangwayShoji and Langway (1985) found an increased flow rate for dirty basal ice from Camp Century, Greenland, compared with clean basal ice; Reference Holdsworth and BullHoldsworth and Bull (1970) observed an enhanced flow rate in the basal amber-colored ice of Meserve Glacier, Antarctica; Reference Song, Cole and BakerSong and others (2004) observed an increased flow rate for laboratory-grown granular ice containing 0.1–4 wt.% particles distributed along the grain boundaries.
Microstructural evolution will strongly affect the deformation process: both the grain-size and the orientations of the grains will affect the deformation rates and mechanisms. Decreasing the grain-size can change the creep mechanism from dislocation glide to a diffusion-based grain boundary sliding mechanism (Reference Goldsby and KohlstedtGoldsby and Kohlstedt, 1997). Reference Jacka and MaccagnanJacka and Maccagnan (1984) indicated that the creep rate can increase by up to a factor of three times the minimum creep rate, and grains show a weak girdle-shaped c-axis orientation after 7.3% strain in uniaxial compression tests. Reference Cuffey, Conway, Gades, Hallet, Raymond and WhitlowCuffey and others (2000a, Reference Cuffey, Thorsteinsson and Waddingtonb) found an enhanced shear strain rate of ice-age ice in southern Greenland; they suggested that this arose from the development of c-axis fabric with a strong shear single maximum pattern and decrease in grain-size, with these two phenomena accounting for roughly 70% and 30% of the average enhancement. Thus, studying microstructural evolution is important to fully understand creep processes.
There are a number of studies (Reference DuvalDuval, 1981; Reference JackaJacka, 1984a; Reference Jacka and MaccagnanJacka and Maccagnan, 1984; Reference Azuma and HigashiAzuma and Higashi, 1985; Reference Wilson and SimWilson and Sim, 2002) concerning fabric development in laboratory-prepared and glacier ice. These studies examined fabric development under compression and shear. The studies indicate that a girdle-shaped polar figure will be formed during low-stress symmetric or nearly symmetric uniaxial compression. Generally, these small girdle-shaped circles are seen when the strain is above 10%. However, Reference Jacka and MaccagnanJacka and Maccagnan (1984) showed that a girdle-shaped circle can be generated even when the strain is as low as 7.3%. Although the stresses (0.1–0.6 MPa) used in the above fabric studies are higher than those observed in almost all natural ice, they might still not be high enough to affect the recrystallization mechanisms. However, compared with fabric development at those stress levels, fabric development under even higher stress may be different because the internal stress generated along the grain boundaries and particles during creep may not have enough time to be released (at high deformation rates). Reference Van der Veen and WhillansVan der Veen and Whillans (1994) developed a simulation for fabric development at a stress of 0.2 MPa and assumed that dynamic recrystallization occurred via strain-induced grain boundary migration (SIBM). They found that the simulation only worked well for short times when the lattice rotations were small. However, under high stress, recrystallization may not occur by SIBM but through the nucleation and growth of many small grains along the grain boundaries (Reference Humphreys and HatherlyHumphreys and Hatherly, 1996). This is due to the high strain energy that is built up along the grain boundaries under high stresses. This process will inevitably increase the number of grains. Furthermore, because strain energy can build up around the particles during deformation, the particles can stimulate the nucleation of new grains. Since the particles are distributed both along the grain boundaries and in the grain interiors, these new grains appear throughout the specimen, not just along the grain boundaries. Studies on metals and alloys (Reference Habiby and HumphreysHabiby and Humphreys, 1999; Reference Huang and HumphreysHuang and Humphreys, 1999; Reference Somerday and HumphreysSomerday and Humphreys, 2003) show that large second-phase particles enhance dynamic recrystallization during deformation by acting as preferred sites for nucleation; this is called particle-stimulated nucleation (PSN). In some previous PSN studies, Reference Bleck and BungeBleck and Bunge (1981) found a nearly random orientation distribution in fully recrystallized grains in an Al–Mn alloy, while Reference Habiby and HumphreysHabiby and Humphreys (1999) showed that in an Al–Si alloy most of the recrystallized nuclei forming around the particles are misoriented from the adjacent deformed matrix by 30–45°.
Our research is motivated by the need to understand how fine particles affect the mechanical behavior of ice. Since the flow processes are inevitably connected with microstructural evolution, the present paper focuses on the effect of dispersed particles, located both in grain boundaries and inside grains, on the creep, dynamic recrystallization and fabric development in laboratory-prepared granular ice specimens.
Experimental
Deionized, distilled and degassed water was used to grow thin plates of large-grained ice. These ice plates were then broken up, and the grains were used to seed granular freshwater ice specimens. The final cylindrical specimen dimensions were 127 mm in length and 50.4 mm in diameter.
For ice specimens with particles, water with silt-sized soil particles (obtained from Hanover, NH, and sieved to 50 ±10 µm) was frozen to a thickness of about 5 mm from bottom to top using a cooling plate. Another layer of ice with particles was then grown on the top of the first layer. This procedure was repeated several times until the thickness of the ice plate was about 30 mm. (It was necessary to grow a multilayer ice plate to distribute the particles throughout the ice grains.) The ice plate was subsequently broken up, and a 3–5 mm sieve fraction was obtained. These grains were then used to seed granular fresh-water ice specimens prepared by the method of Reference ColeCole (1979). Both the particle-free ice and the ice with particles had an average initial grain-size of approximately 5 mm. Figure 1 shows a thin section and an enlarged view of ice with 1 wt.% (~0.43 vol.%) particles. It can be seen that the particles are uniformly distributed (Fig. 1b). Unlike a previous study (Reference Song, Cole and BakerSong and others, 2004), in which the particles were located only along the grain boundaries, here the particles are distributed both along the grain boundaries and inside the grain interior.
Constant-load creep tests with an initial stress level of 1.45 MPa at −10°C were used to study the effect of particles on the creep behavior. After the creep strain reached ~10%, the stress had decreased to ~1.3 MPa, i.e. about 10% lower than the initial stress, due to the increase in the specimen’s cross-sectional area. A load cell mounted in line with the specimens monitored the axial load, and two displacement transducers mounted on circumferential rings (located at one-third and two-thirds along the length of the specimens) provided the deformation measurements.
The experiments employ cyclic loading to determine the dislocation density based on a dislocation-drag model (Reference Cole and DurellCole, 1995; Reference Cole and DurellCole and Durell, 2001). During cyclic tests, a zero-mean-stress sinusoidal waveform was applied to the specimens by an electro-hydraulic actuator. The experiments employed the reversed direct-stress testing method developed by Cole and Durell (1995), and the procedures described by Reference Cole, Dempsey, Bazant, Rajapakse and SunderCole (1993). These tests also show the grain boundary relaxation peak that occurs at high frequencies. Thus, to study the effects of particles on the grain boundary sliding, cyclic loading was applied to the specimens to determine the internal friction. The internal friction (the ratio of the hysteresis loop area to the peak strain energy during a load cycle) is a dimensionless measure of the energy dissipated per cycle.
After creep (strains of ~1–10%) the specimens were thin-sectioned in order to study the effects of particles on both dynamic recrystallization and fabric development. Microstructural analysis was performed using optical microscopy with polarized light. Fabric studies used the method developed by Reference LangwayLangway (1958) to determine the c-axis orientation to ±5°. The orientation of each specimen was measured by two workers (Song and E. Andreas) to ensure the accuracy of the fabric studies.
Results and Discussion
Creep
Figure 2 shows the relationship between strain rate and viscous strain for both the particle-free ice and the ice with 1 wt.% particles. It should be noted that the viscous strain does not include the elastic and anelastic components, which have been removed from the total strain because these two components will not affect the microstructure. Note that the creep rate of the ice with particles was always higher than that of the particle-free specimens, even at strains that are too low for appreciable recrystallization to occur. This agrees well with a previous study (Reference Song, Cole and BakerSong and others, 2004), in which particles are distributed only along the grain boundaries and also show increased creep rates, compared with those of particle-free ice. The higher creep rate of ice with particles is possibly due to an increased dislocation density caused by the particles. Since under these conditions (−10°C and 1.45 MPa) the main creep mechanism is dislocation glide on the basal plane (Reference Duval, Ashby and AndermanDuval and others, 1983), an increase in dislocation density will inevitably increase the creep rate. The curves in Figure 2 show typical creep behavior: the creep rate initially decreases (primary creep) to a minimum value and then increases (tertiary creep). During tertiary creep, dynamic recrystallization leads to increases in the creep rate. It can be seen that the strain at the minimum creep rate is about 0.01 for both particle-free ice and ice with particles, in agreement with Reference JackaJacka’s data (1984b).
Microstructural evolution
Figure 3 shows thin sections of the particle-free ice observed using polarized light before and after creep to the indicated strains. Note the refinement, i.e. reduction, in grain-size with increasing strain. Figure 4 shows higher-magnification images, where it is evident that dynamic recrystallization has occurred along the grain boundaries (see arrows). In the early stages of creep (strain ~1%), few new grains nucleate. With increasing strain, the effect of dynamic recrystallization increases until ~10% strain, at which point the specimen is substantially recrystallized. New grains replace the initial microstructure, resulting in a decrease in the average grain-size from ~5 mm to ~2 mm after 10% strain.
Figure 5 shows the thin sections of ice with 1 wt.% particles observed using polarized light both before and after creep to the indicated strains. Again, the grain-size decreases with increasing strain, but more rapidly than in the particle-free ice. Note also the heterogeneity of the microstructure, with some large, probably original, grains in a largely recrystallized matrix. Figure 6 shows enlarged images of the microstructure. Similar to the particle-free ice, new grains nucleate along the grain boundaries, but at a higher nucleation rate (cf. Figs 3–6). Even after only ~1% strain, the extent of dynamic recrystallization is very high in the ice with particles, with many new grains nucleated along the grain boundaries (see arrows in Fig. 6). At ~10% strain, the whole microstructure has largely recrystallized, and almost all of the old grains are replaced by the new, small, recrystallized grains. This particle-induced dynamic recrystallization was also observed in a previous study (Reference Song, Cole and BakerSong and others, 2004), which focused on the effect of particles on creep in microstructures with particles distributed only along the grain boundaries.
The particles (both located along the grain boundaries and inside the grain interiors) produce greater strain energy due to local lattice rotations (Reference Humphreys, Hansen, Jones and LeffersHumphreys, 1980) than occurs in the particle-free ice, and thus will increase the nucleation rate. This produces a higher recrystallization rate and leads to a greater reduction in grain-size with increasing strain than in the particle-free ice. Table 1 shows the relationship between average grain-size and strain level for particle-free ice, ice with 1 wt.% particles distributed both along the grain boundaries and in the grain interior, and ice with 1 wt.% particles distributed only along the grain boundaries for comparison. It can be seen that as the strain increases, the average grain-size decreases, with the ice with particles showing a much larger grain-size reduction than the particle-free ice. Ice with 1 wt.% particles distributed both along the grain boundaries and in the grain interiors, and ice with 1 wt.% particles distributed only along the grain boundaries show similar grain-size reduction, indicating that the particles in the grain interiors may have effects on dynamic recrystallization similar to those of the particles along the grain boundaries.
Typically, pure polycrystalline ice exhibits an internal friction peak (Reference Cole and DurellCole, 1995) due to grain boundary sliding located at ~5 Hz at −10°C. Figure 7 shows how particles affect the internal friction before creep. It is evident that the internal friction of the ice with particles is similar to that of particle-free ice at frequencies less than 0.1 Hz, but it is at least 30% lower at 1 Hz. This result is in accordance with a previous study (Reference Song, Cole and BakerSong and others, 2004), which focused on the effect of particles that are located only along the grain boundaries (also see Fig. 7). When particles are uniformly distributed throughout the microstructure, some fraction of them will be located along the grain boundaries. The decrease in internal friction at 1 Hz for ice with particles indicates that there were a sufficient number of particles along the grain boundaries to block grain boundary sliding and thus eliminate the grain boundary relaxation peak. Reference Baker, Liu, Jia, Hu, Cullen and DudleyBaker and others (2000) showed that for polycrystalline fresh-water ice, dislocation nucleation always occurs at the grain boundaries in response to the stresses built up by grain boundary sliding. When particles inhibit grain boundary sliding, the strain energy will increase along the grain boundaries, making them favorable sites for the nucleation of new grains.
Fabric development
Figure 8 shows the c-axis orientations for particle-free ice before and after creep to the indicated strains. It is evident that the random distribution of c-axis orientations before creep slowly changes to a weak, small girdle-shaped circle when the strain reaches 10%. Actually, during the compressive creep tests, the slip direction rotates away from the compression axis if one considers each crystal individually. Compared with some previous data, which show a strong girdle-shaped orientation at strains above 15% (Reference DuvalDuval, 1981; Reference Azuma and HigashiAzuma and Higashi, 1985), the ~10% strain is rather small to generate such a girdle-shaped fabric in our tests. The reason for this phenomenon is not clear at this stage. However, since our tests use a high stress level, the strain energy developed is also very high, which might help form the girdle-shaped orientation at a lower strain. Reference Jacka and MaccagnanJacka and Maccagnan (1984) also show that a girdle shape can form even with only 7.3% strain during uniaxial compression tests with a stress level of 0.2 MPa.
Figure 9 shows the c-axis orientations for ice with 1 wt.% particles before and after creep to the indicated strains. It can be seen that the grains show random orientation at all strain levels, and the fabric after ~10% strain is much weaker than that for particle-free ice with ~10% strain. Since the microstructure (Figs 5 and 6) shows that almost all the deformed grains were replaced by the new, recrystallized, small grains, the random orientation at ~10% strain indicates that the new grains do not have preferred orientations. The rather random c-axis orientations for ice with particles also suggests that dynamic recrystallization occurs through nucleation of new grains rather than through strain-induced boundary migration (Reference Gottstein and ShvindlermanGottstein and Shvindlerman, 1999) or through polygonization by edge dislocation movement, because the latter two processes will only change the orientation of grains by a few degrees. Thus, if dynamic recrystallization occurs through either of those two processes, a weak small-circle girdle, as found for deformed particle-free ice, should be observed.
Previous studies (Reference Gow and WilliamsonGow and Williamson, 1971, Reference Gow and Williamson1976) indicated the existence of dust bands at 1200–1800 m depth in the Byrd (Antarctica) ice cores. These debris bands are estimated to have been deposited in the ice sheet 43 500–7500 years ago and consist of volcanic ash. Reference Gow and WilliamsonGow and Williamson (1971) observed that the grain-size decreases dramatically and the c axes of the fine grains show a tight single-pole fabric. In our laboratory results on particle-containing ice, grain-size decreases dramatically compared with fresh-water ice during creep. However, since our ice with particles did not have enough time and deformation to develop the strong fabric that was evident in the dust-band ice cores (Reference Gow and WilliamsonGow and Williamson, 1971), the difference in fabric development is not surprising.
Conclusions
The effect of particles on the creep, dynamic recrystallization and fabric development in granular fresh-water ice has been investigated and the following conclusions are drawn:
Ice with 1 wt.% particles shows a higher creep rate at all stages of creep than particle-free ice.
New, small, grains nucleate along the grain boundaries for both the particle-free ice and particle-laden ice.
Particles inhibit grain boundary sliding and enhance dynamic recrystallization through the development of high local strain energies, leading to particle-stimulated nucleation and dramatically decreasing the grain-size.
The fabric for particle-free ice shows a weak girdle shape. Particles inhibit the development of this fabric by nucleating grains with nearly random orientations.
Acknowledgements
This research was supported by the US National Science Foundation Office of Polar Programs, Arctic Natural Sciences Program (OPP 011737). We thank G. Durell for his valuable assistance in developing the creep loading equipment and E. Andreas for her valuable help with the orientation measurements.