Description: An overview of the deformation and recrystallization mechanisms that are active in the North Greenland Eemian Ice Drilling (NEEM) ice core is given, based on microscale models, light microscopy and cryogenic electron backscatter diffraction (cryo-EBSD). The Holocene ice (0-1419 m depth) deforms by dislocation creep with basal slip accommodated by non-basal slip. The amount of non-basal slip is controlled by the extent of strain induced boundary migration (SIBM). The most important recrystallization mechanisms and processes in the Holocene ice are grain dissection, strain induced boundary migration (SIBM), and bulging nucleation. In the glacial ice (1419-2207 m of depth) basal slip is accommodated by both non-basal slip and grain boundary sliding (GBS). Rotation recrystallization is more important, while SIBM is less important in the glacial ice compared to the Holocene ice. In the Eemian ice (2207-2540 m depth), which is at high temperature, different microstructures occur depending on the impurity content of the ice. The difference in microstructure and deformation mechanisms, between interglacial and glacial ice can have important consequences for ice rheology and ice sheet dynamics.

Description: The ice microstructure in the lower part of the North Greenland Eemian Ice Drilling (NEEM) ice core consists of relatively fine-grained ice with a single maximum crystallographic preferred orientation (CPO) alternated by much coarser-grained ice with a partial (great circle) girdle or multi-maxima CPO. In this study, the grain-size-sensitive (GSS) composite flow law of Goldsby and Kohlstedt (2001) was used to study the effects of grain size and premelting (liquid-like layer along the grain boundaries) on strain rate in the lower part of the NEEM ice core. The results show that the strain rates predicted in the fine-grained layers are about an order of magnitude higher than in the much coarser-grained layers. The dominant deformation mechanisms, based on the flow relation of Goldsby and Kohlstedt (2001), between the layers is also different, with basal slip rate limited by grain boundary sliding (GBS-limited creep) being the dominant deformation mechanism in the finer-grained layers, while GBS-limited creep and dislocation creep (basal slip rate limited by non-basal slip) contribute both roughly equally to bulk strain in the coarse-grained layers. Due to the large difference in microstructure between finer-grained ice and the coarse-grained ice at premelting temperatures (T〉262 K), it is expected that the fine-grained layers deform at high strain rates, while the coarse-grained layers are relatively stagnant. The difference in microstructure, and consequently in viscosity, between impurity-rich and low-impurity ice can have important consequences for ice dynamics close to the bedrock.

Description: The effect of grain size on strain rate of ice in the upper 2207 m in the North Greenland Eemian Ice Drilling (NEEM) deep ice core was investigated using a rheological model based on the composite flow law of Goldsby and Kohlstedt (1997, 2001). The grain size was described by both a mean grain size and a grain size distribution, which allowed the strain rate to be calculated using two different model end-members: (i) the microscale constant stress model where each grain deforms by the same stress and (ii) the microscale constant strain rate model where each grain deforms by the same strain rate. The model results predict that grain-size-sensitive flow produces almost all of the deformation in the upper 2207 m of the NEEM ice core, while dislocation creep hardly contributes to deformation. The difference in calculated strain rate between the two model end-members is relatively small. The predicted strain rate in the fine-grained Glacial ice (that is, ice deposited during the last Glacial maximum at depths of 1419 to 2207 m) varies strongly within this depth range and, furthermore, is about 4–5 times higher than in the coarser-grained Holocene ice (0–1419 m). Two peaks in strain rate are predicted at about 1980 and 2100 m depth. The prediction that grain-size-sensitive creep is the fastest process is inconsistent with the microstructures in the Holocene age ice, indicating that the rate of dislocation creep is underestimated in the model. The occurrence of recrystallization processes in the polar ice that did not occur in the experiments may account for this discrepancy. The prediction of the composite flow law model is consistent with microstructures in the Glacial ice, suggesting that fine-grained layers in the Glacial ice may act as internal preferential sliding zones in the Greenland ice sheet.

Description: Apart from evaluating the crystallographic orientation, focus of microstructural analysis of natural ice during the last decades has been to create depth-profiles of mean grain size. Several ice flow models incorporated mean grain size as a variable. Although such a mean value may coincide well with the size of a large proportion of the grains, smaller/larger grains are effectively ignored. These smaller/larger grains, however, may affect the ice flow modeling. Variability in grain size is observed on centimeter, meter and kilometer scale along deep polar ice cores. Composite flow laws allow considering the effect of this variability on rheology, by weighing the contribution of grain-size-sensitive (GSS, diffusion/grain boundary sliding) and grain-size-insensitive (GSI, dislocation) creep mechanisms taking the full grain size distribution into account [Herwegh et al., 2005, J. Struct. Geol., 27, 503-521]. Extraction of hundreds of grain size distributions for different depths along an ice core has become relatively easy by automatic image processing techniques [T. Binder et al., 2013, J. Microsc., 250, 130-141]. The shallow ice approximation is widely adopted in ice sheet modeling and approaches the full-Stokes solution for small ratios of vertical to horizontal characteristic dimensions. In this approximation shear stress in the vertical plain dominates the strain. This assumption is not applicable at ice divides or dome structures, where most deep ice core drilling sites are located. Within the upper two thirds of the ice column longitudinal stresses are not negligible and ice deformation is dominated by vertical strain. The Dansgaard-Johnsen model [W. Dansgaard & S.J. Johnsen, 1969, J. Glaciol., 8, 215-223] predicts a dominating, constant vertical strain rate for the upper two thirds of the ice sheet, whereas in the lower ice column vertical shear becomes the main driver for ice deformation. We derived vertical strain rates from the upper NEEM ice core (North-West Greenland) and compared them to classical estimates of strain rates at the NEEM site. Assuming intervals of constant accumulation rates, we found a variation of vertical strain rates by a factor 2-3 in the upper ice column. We discuss the current applicability of composite flow laws to grain size distributions extracted from ice cores drilled at sites where the flow direction rotates by 90 degrees with depth (i.e. ice divide). An interesting finding is that a transition to a glacial period in future would be associated with a decrease in vertical strain rate (due to a reduced accumulation rate) and an increase of the frequency of small grains (due to an enhanced impurity content). Composite flow laws assign an enhanced contribution of GSS creep to this transition. It is currently unclear which factor would have a greater influence.