Hostname: page-component-78c5997874-94fs2 Total loading time: 0 Render date: 2024-11-19T05:18:20.430Z Has data issue: false hasContentIssue false

Formation Processes of Ice Fabric Pattern in Ice Sheets

Published online by Cambridge University Press:  20 January 2017

N. Azuma
Affiliation:
Department of Applied Physics, Faculty of Engineering, Hokkaido University, Sapporo, Japan 060
A. Higashi
Affiliation:
Division of Natural Sciences, International Christian University, Mitaka-shi, Tokyo, Japan 181
Rights & Permissions [Opens in a new window]

Abstract

Experimental studies on the Formation process of the single maximum fabric of ice from its random fabric were carried out using both Greenland Dye-3 ice cores and artificially grown polycrystalline ice. From the observations of the rotation of c-axes of individual grains in thin specimens, an empirical equation was obtained for describing the rotation process of c-axes. Combined with another empirical equation on the relationship between the uniaxial strain of whole specimen and the uniaxial strain in individual grains of different c-axis direction, an equation was derived which describes the formation of single maximum fabric in bulk ice samples.

Computer simulation of the formation process was performed based on an equation starting from an artificially made random fabric pattern. From the results obtained by a stepwise simulation of 2% strain intervals, the degree of concentration expressed by θm, the angle which halves the number of grains in the fabric diagram is correlated with the accumulated strain in bulk ice (Fig. 8). The strain deduced by applying this relationship to the obtained θm for samples of Dye-3 ice core coincided well with actual strain which the samples had experienced. This implies that the equation is valid and the mechanism of c-axis rotation prevails in polar ice sheets.

Type
Research Article
Copyright
Copyright © International Glaciological Society 1985

INTRODUCTION

The formation mechanism of the characteristic fabric pattern of ice in various ice masses on the earth is of vital importance for interpretations of the flow history of ice masses, and for prediction of changes in properties and behavior of ice. It is generally recognized that the fabric pattern of ice in polar ice sheets gradually changes from a random pattern near the surface to the single maxima at depth.

Laboratory experiments have revealed the kinds of characteristic fabrics that can be obtained under certain conditions of the deformation of ice; the small girdle fabric pattern was obtained by the uniaxial compression test (Steinemann 195S, Watanabe and Ôura 1968, Tanaka 1972 and Huang, Ohtomo and Wakahama 1983), and a single maximum or two-maximum fabric was developed by simple shear experiments (Shumskii 1958, Rigsby I960, Kamb 1972 and Duval 1981). However, no detailed experiments or analyses have been carried out clarifying mechanisms of the development of preferred orientation of c-axes in various experiments.

In general, the following can be considered as the three main mechanisms for the development of the ice fabric:

  • 1) rotation of basal plane acompanied by deformation

  • 2) grain growth

  • 3) recrystallization

Although processes in natural ice masses should be governed by combinations of these mechanisms, it is desirable to study experimentally how these mechanisms contribute the formation of characteristic fabric pattern. In polar ice sheets where ice temperatures are very low, mechanisms (1) and (2) are considered to prevail over (3). With respect to (1), it has been reported that the c-axis rotated towards the compression axis when a single crystal of ice was subjected to uniaxial compression (Jones and Brunet 1978), and also when a thin specimen of polycrystalline ice was compressed (Wakahama 1964). However, the rotation angle of the c-axis was not correlated to the amount of uniaxial strain and this phenomenon was not considered to be related to the process of the fabric formation.

In the present work, the quantitative relationship between the rotation angle of c-axis and the amount of uniaxial strain was obtained from uniaxial compressive deformation experiments using thin specimens of polycrystalline ice. Deformation experiments of the same type were carried out for bulk specimens cut from an ice core taken from shallow depth in Greenland ice sheet. Random fabric of this core changed to the broad single maximum when the uniaxial strain reached 10%, and it is inferred that the same mechanism of c-axis rotation takes place in bulk polycrystalline ice when neither grain growth nor recrystallization take place in the process. Computer simulations were performed for the production of single maximum fabric using a constitutive equation describing the relationship between the rotation angle of c-axis and the amount of uniaxial strain.

EXPERIMENTAL PROCEDURES

Specimen preparation

In order to observe the rotation process of the C-axis of individual grains during the compressive deformation, thin specimens 2 mm thick and 30 x 20 mm2 in area were prepared from both Dye-3 ice core and artificially grown polycrystalline ice. The latter was grown from the surface of pure water at a moderate rate, so that the c-axes became horizontal under a certain depth. When thin specimens were cut horizontally from this ice, it had two-dimensional random orientation fabric of c-axes. Its average grain size was 2 mm.

In the case of Dye-3 ice core (sample 1851, 1900 m depth and sample 1950, 2000 m depth), thin specimens were cut from the core so that the plane of the wide area was parallel to the core axis and its long axis made 45° with the core axis. Since both cores have the single maximum fabric, the c-axis of each grain lies on the plane of the thin specimen and angles between the compressive axis and the mean direction of c-axes of both cores are 40° and 60° respectively, with a standard error of 10°. Average grain sizes for the cores of sample 1851 and sample 1950 were 1.5 mm and 2.7 mm respectively.

For a study of the effect of compressive deformation on the random fabric, normal compressive tests were carried out with bulk specimens 20 x 20 x 70 mm2. They were cut from Dye-3 ice core no. 120 (201 m depth) so that the long axis of the specimen was parallel to the core axis, and therefore, the specimen’s compressive axis in the test coincided with that of the natural compressive stress in its original depth in the ice sheet. The c-axes, fabric of this core is almost random and its average grain size is approximately 2 mm.

Fig. 1. Photographs of a thin specimen through crossed polaroids during compressive deformation test no. 2. Axial strain for each is: 1) 0%, 2) 3.4%, 3) 6.5%, 4) 10.0%, 5) 11.5%, 6) 12.4%. Compression axis is in the top and bottom direction of the figure.

Uniaxial compression tests of thin specimens

Processes of the c-axis rotation were recorded on photographs of thin specimens intermittently taken on the process of uniaxial compression tests by transmitted light through the crossed polaroids which sandwiched the specimen. Thin sections of approximately 2 mm revealed clear images of individual grains and even of slip bands when they appeared by stress, as can be seen in Figure 1. Every compression test was carried out using a testing machine (for details of which see Higashi and Shoji 1979) under an atmospheric pressure or 20 MPa hydrostatic pressure with a constant strain rate in the order of 10-7s-1. Experimental conditions including temperatures and the characteristics of the used specimens are given in Table 1. The thin specimen was also sandwiched by transparent acryllic-resin boards to avoid the buckling which may take place when such a specimen is compressed perpendicular to the thin section. This arrangement gave plane strain conditions.

Since it is expected that the rotation of the c-axis in an individual grain is related to internal strain, the internal local uniaxial strain of each grain was determined by the relative movement of several air bubbles in it (black dots in Figure 1). The least-square method of evaluating the plane strain using several data is developed from the standard method (Nye 1957a).

Uniaxial compression lest of bulk specimens

Uniaxial compression tests with bulk specimens for the study of fabric change were carried out with constant strain rate of 4.38 x 10-7s-1, under an atmospheric pressure. Experimental temperatures were -13 C (test 5) and -20 °C (test 6). When the total strain reached 10%, the test was terminated; the specimen was immediately removed from the testing machine and horizontal thin section specimens were prepared . Such procedures were carried out in a -20 °C cold room.

EXPERIMENTAL RESULTS

Rotation of c-axes

Figure 1 shows six successive photographs of a thin section specimen taken before the deformation and when the total strain of the specimen reached 3.4% 6.5%, 10.0%, 11.5% and 12.4% respectively in no. 2 test (Table 1). The slip bands appear in individual grain even when the total strain reached only 2% and they became clear as the strain increased. Since the c-axis must be perpendicular to the slip plane in each grain and it lay parallel to the plane of the wide area or to the plane of the figure, we can directly determine the c-axis of individual grains from these photographs. Angle x between the specimen’s compressive axis and the direction of the c-axis in each grain was measured using these determined c-axis direction at various increasing value of strain. We will call x the misorientation angle of a grain. Decreasing tendency of x with increase of the strain can be clearly seen in Figure 2 in which several examples of the change of x from different values x of various grains are plotted against the total strain E of specimens. They indicate that the c-axis rotates in individual grains as the deformation proceeds and the larger the value χ0 the more rapid the rotation rate.

Table 1.

Fig. 2. Changes of the misorientation angle x with increasing strain ε. Several grains of different initial angle x0 are designated by different symbols. Solid lines show the relationship x vs. ε by the equation 3.

Fig. 3. Change of the nearly random fabric of Dye-3 core no.120 a) before the compression test, b) after test no. 5, c) after test no.6, and d) the simulated fabric diagram by equation 4, starting from a) and strained 10%. The center of each diagram coincides with the compression axis.

Change of the random fabric

Fabric diagrams for the originally undeformed, and deformed specimens in the tests at -13 C and -20 C are illustrated in circles of three diagrams of Figure 3. This figure shows clearly that the c-axes become concentrated after the deformation, a tendency quantitatively expressed in the following way. Accumulated fraction curves as shown in Figure 4 are drawn by plotting the accumulated fraction F(e) of number of c-axes within certain misorientation angle θ against the values of θ divided into S intervals in this case. Three curves in Figure 4 exhibit F(e) for three data sets corresponding to Figure 3(a), (b) and (c). The almost straight solid line for the initial random fabric change to those broken one for test no 5 and dotted one for test no 6. They exhibit a steep slope below 40° and gentle slope at larger angles than 60°, indicating the tendency that the fabric becomes concentrated toward x approaches zero as the total strain increases. This process will be discussed and simulated in the next section with the aid of a formula obtained from the c-axes rotation process in the preceding subsection.

DISCUSSION

italic>Rotation of c-axes

In the case of the compression test of a single crystal under restraint, it is well known that the normal to a slip plane rotates as it tends to become parallel to the compression axis (Schmid and Boas 1935). Since the basal plane is a dominant slip plane in ice crystal, the angle x between the c-axis and the compression axis

Fig. 4. Cumulative fraction F(θ) curve for the fabrics of Dye-3 #120 before the test (----- ), after test S (- - - ) and test 6 ( . . . ).

decreases from its initial value compressive strain e equation:

(1)

Individual grains in our thin specimens are considered to be restrained by their neighbours when they are deformed. Therefore, it is expected that the rotation of the c-axis of each grain could be correlated to the local compressive strain e_ in the grain in such a way that ε is inserted for εS in eq. (1). εg in individual grains may be proportional to the gross axial strain ε of the polycrystalline plane specimen with a proportional constant as a function of Xo Figure 5 shows relationships between the measured εg of individual grains and their initial orientation angles Xo at two stages when the total axial strain reached 3% and 7%. Although data scatter considerably, they fit to the curve of an equation of the form:

Fig. 5. Relationship between the unixaial compressive strain εg of individual grains which had initial misorientotior. angle X0 after a certain total strain, 3% and 7%. Solid lines show equation 2.

(2)

This equation implies that the arithmetic mean of εg over; Xo from O° to 90° is equal to the total axial strain £. Details of this implication together with its expansion to the three-dimensional case will be given elsewhere.

Inserting εg of this equation for εs in eq. (1), we obtain the changed orientation χ related to the total axial strain e and the initial angle ε of the orientation in the following form:

(3)

The solid lines in Figure 2 are the curves calculated by this equation (3) started from several certain initial χ0 s. Good coincidence of the curves on the plotted empirical data of x vs. strain suggests that the rotation mechanism of c-axes envisaged to obey equation (I) is valid.

Formation of the broad single maximum fabric

As was described in the preceding subsection, the nearly random fabric of ice at comparatively shallow core of Dye-3 tends to concentrate into the broad single maximum when its bulk specimen was uniaxially compressed. In view of the fact that the c-axes of individual grains rotate so that the misorientation angle tends to decrease to zero degree in a thin polycrystalline specimen during the compressive deformation as described above, it is inferable that the same mechanism takes place for the concentration process of the c-axes in the bulk polycrystalline ice. Therefore, it can be expected that the broad single maximum fabric will be produced in nature when polycrystalline ice body is subjected to the compressive deformation as the case in the upstream accumulation zone of an ice sheet.

For predicting the degree of c-axes concentration owing to the deformation, two computer simulations on the formation process of the broad single maximum fabric were carried out. The constitutive equation used for these simulations is as follows:

(4)

This equation is of the same type as equation (3) because it should express the decreasing tendency of the misorientation angle θ with increasing strain. The angle θ is equivalent to x in equation (3) but express the angle between the c-axis of individual grain and the compressive axis in a bulk polycrystalline specimen. Since the fabric pattern varies with increasing strain, the simulation should be performed stepwise with a certain strain interval Δε. The constant Π in equation (3) is substituted by a parameter A which also varies with change of the fabric pattern.

The value of A can be determined by the following equation which is derived from an empirical fact that the uniaxial strain of a bulk specimen e is equal to the arithmetic mean of the uniaxial strain e of whole individual grains:

(5)

where N is the total number of grains of which misorientation angle e0, are measured. Angle Θ in equation (4) is the initial value of θ at each step of the simulation: the real initial value for each grain at the start and the resultant of the preceding step for the next step. Different values of A near 3 according to equation (5) were used in successive steps. In practice for the present, Δε was taken as 2%,

The one simulation started from the fabric of Dye-3 no. 120 core which exhibited nearly random fabric as illustrated in Figure 3a. Fabric diagram shown in Figure 3d is the resultant fabric when ε reached 10%. It is comparable to those obtained by compression tests as illustrated in Figures 3b and 3c. Another simulation was carried out starting from an artificially-made random fabric as shown in Figure 6a in which 200 orientation of c-axes were selected by the use of random number generated in the computer. Results of the simulation at 10%, 30% and 50% of the total compressive strain are given in Figures 6b, 6c and 6d.

Fig. 6. Simulated fabric diagrams started from artificially made random fabric with various degrees of uniaxial compressive strain. The compressive axis is in the center of each diagram.

Fig. 7. F(θ) - θ curve for the simulated fabrics shown in Fig. 6.

F(θ) curves, the accumulated fractions of the c-axes within θ are drawn in Figure 7 from such calculations as performed for the above simulation. The figure clearly shows that the broadness of the c-axes concentration decreases as the total strain increases. If we draw a horizontal line at a level of 50% of F(a) in Figure 7, its intersections with the predicted curve of F(θ) at various total strain ε give us the angle θ or the inclination from the compression axis at which the number of grains are halved in the fabric diagram. This is called the median inclination θm and a good indicator of the degree of concentration of the c-axes in the fabric. The strain dependence of θm obtained from results of the above simulation is given as a solid line in Figure 8.

Fig. 8. Relationship between the median inclination θm and the vertical strain for the Dye-3 core ice and specimens of compression tests. The solid line shows the strain dependence of θm obtained from results of simulated fabrics.

Individual relationships between the median inclination θm for the observed fabrics of specimens of ice core from various depths at Dye 3, Greenland and their given vertical strain ε are plotted on this figure (closed circles). The data for specimens of compression test no 5 and no 6 are also plotted (open circles). The vertical strain given to the ‘ice at a certain depth was estimated by a simple calculation based on Nye’s flow model (Nye 1957b), assuming that the ice thickness and the accumulation rate has not changed during the past few thousand years. The accumulation rate used for the calculation is the measured value of 0.54 ma-1 at Dye 3 (Reeh and others 1978). For specimens of the uniaxial compression tests, the vertical strain was given by adding the compressive strain of 10% to the calculated vertical strain of 10% corresponding to the depth 201 m of specimen no. ¡20.

Field data and empirical results plotted on the figure fit fairly well with the solid line. This proves that the formation mechanism of the broad single maximum of the fabrics of ice cores from the depth of plateau area of Greenland ice sheet is principally the rotation of c-axes in individual grains composing the ice in it. This is believed to be generally true in polar ice sheets in which recrystallization does not occur intensively because of low temperatures.

ACKNOWLEDGEMENTS

We thank Professor C C Langway, Jr for providing ice core samples of Dye-3, and Dr T Hondoh and Dr M Nakawo for useful discussion while preparing this paper,

References

Duval, P 1981 Creep and Fabrics of polycrystalline ice under shear and compression. Journal of Glaciology 27(95): 129-140 Google Scholar
Higashi, A, Shoji, H 1979 [Mechanical tests of Antarctic deep core ice under hydrostatic pressure; instrumentation and preliminary results]. Oyo-butsuri 48(1): 41-17 Google Scholar
Huang, M, Ohtomo, M, Wakahama, G 1982 [Repeated recrystallization of polycrystalline ice]. Low Temperature Science A41: 13-23 Google Scholar
Jones, S J, Brunet, J G 1978 Deformation of ice single crystals close to the melting point. Journal of Glaciology 21(85): 445-155 Google Scholar
Kamb, B 1972 Experimental recrystallization of ice under stress. Geophysical Monograph 16: 211-352 Google Scholar
Nye, J F 1957a Physical properties of crystals. Oxford, Oxford University Press Google Scholar
Nye, J F 1957b The distribution of stress and velocity in glaciers and ice sheets. Proceedings of the Royal Society of London A, 239(1216): 113-133 Google Scholar
Reeh, N, Clausen, H B, Dasgaard, W, Gundestrup, N, Hammer, C U, Johnsen, S J 1978 Secular trends of accumulation rates at three Greenland stations. Journal of Glaciology 20(82): 27-30 CrossRefGoogle Scholar
Rigsby, G P 1960 Crystal orientation in glacier and experimentally deformed ice. Journal of Glaciology 3(27): 589-606 Google Scholar
Schmid, E, Boas, W 1935 Kristallplastizilat. Berlin, Springer-Verlag Google Scholar
Shumskii, P A 1958 The mechanism of ice straining and its recrystallization. Symposium of Chamonix. IASH Pubi 47: 244-248 Google Scholar
Steinemann, S 1958 Experimentalle Untersuchungen zur Piastizitat von Eis. Beitrage zur Geologie der Schweiz. Geotechnische Serie. Hydrologie, 10Google Scholar
Tanaka, H 1972 On preferred orientation of glacier and experimentally deformed ice. Journal of Geological Society of Japan 78(12): 659-675 Google Scholar
Wakahama, G 1964 [On the plastic deformation of ice. V. Plastic deformation of polycrystalline ice). Low Temperature Science. A22: 1-24 Google Scholar
Watanabe, O, Oura, H 1968 [Experimental studies on orientation of polycrystalline ice by unconfined compression 1]. Low Temperature Science, A26: 1-28 Google Scholar
Figure 0

Fig. 1. Photographs of a thin specimen through crossed polaroids during compressive deformation test no. 2. Axial strain for each is: 1) 0%, 2) 3.4%, 3) 6.5%, 4) 10.0%, 5) 11.5%, 6) 12.4%. Compression axis is in the top and bottom direction of the figure.

Figure 1

Table 1.

Figure 2

Fig. 2. Changes of the misorientation angle x with increasing strain ε. Several grains of different initial angle x0 are designated by different symbols. Solid lines show the relationship x vs. ε by the equation 3.

Figure 3

Fig. 3. Change of the nearly random fabric of Dye-3 core no.120 a) before the compression test, b) after test no. 5, c) after test no.6, and d) the simulated fabric diagram by equation 4, starting from a) and strained 10%. The center of each diagram coincides with the compression axis.

Figure 4

Fig. 4. Cumulative fraction F(θ) curve for the fabrics of Dye-3 #120 before the test (----- ), after test S (- - - ) and test 6 ( . . . ).

Figure 5

Fig. 5. Relationship between the unixaial compressive strain εg of individual grains which had initial misorientotior. angle X0 after a certain total strain, 3% and 7%. Solid lines show equation 2.

Figure 6

Fig. 6. Simulated fabric diagrams started from artificially made random fabric with various degrees of uniaxial compressive strain. The compressive axis is in the center of each diagram.

Figure 7

Fig. 7. F(θ) - θ curve for the simulated fabrics shown in Fig. 6.

Figure 8

Fig. 8. Relationship between the median inclination θm and the vertical strain for the Dye-3 core ice and specimens of compression tests. The solid line shows the strain dependence of θm obtained from results of simulated fabrics.