NEEM Deep Ice Core Research Articles

Using a composite flow law to model deformation in the NEEM deep ice core, Greenland – Part 1: The role of grain size and grain size distribution on deformation of the upper 2207 m

Abstract. 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.

Using a composite flow law to model deformation in the NEEM deep ice core, Greenland – Part 2: The role of grain size and premelting on ice deformation at high homologous temperature

Abstract. 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.

EBSD analysis of subgrain boundaries and dislocation slip systems in Antarctic and Greenland ice

Abstract. Ice has a very high plastic anisotropy with easy dislocation glide on basal planes, while glide on non-basal planes is much harder. Basal glide involves dislocations with the Burgers vector b = 〈a〉, while glide on non-basal planes can involve dislocations with b = 〈a〉, b = [c], and b = 〈c + a〉. During the natural ductile flow of polar ice sheets, most of the deformation is expected to occur by basal slip accommodated by other processes, including non-basal slip and grain boundary processes. However, the importance of different accommodating processes is controversial. The recent application of micro-diffraction analysis methods to ice, such as X-ray Laue diffraction and electron backscattered diffraction (EBSD), has demonstrated that subgrain boundaries indicative of non-basal slip are present in naturally deformed ice, although so far the available data sets are limited. In this study we present an analysis of a large number of subgrain boundaries in ice core samples from one depth level from two deep ice cores from Antarctica (EPICA-DML deep ice core at 656 m of depth) and Greenland (NEEM deep ice core at 719 m of depth). EBSD provides information for the characterization of subgrain boundary types and on the dislocations that are likely to be present along the boundary. EBSD analyses, in combination with light microscopy measurements, are presented and interpreted in terms of the dislocation slip systems. The most common subgrain boundaries are indicative of basal 〈a〉 slip with an almost equal occurrence of subgrain boundaries indicative of prism [c] or 〈c + a〉 slip on prism and/or pyramidal planes. A few subgrain boundaries are indicative of prism 〈a〉 slip or slip of 〈a〉 screw dislocations on the basal plane. In addition to these classical polygonization processes that involve the recovery of dislocations into boundaries, alternative mechanisms are discussed for the formation of subgrain boundaries that are not related to the crystallography of the host grain.The finding that subgrain boundaries indicative of non-basal slip are as frequent as those indicating basal slip is surprising. Our evidence of frequent non-basal slip in naturally deformed polar ice core samples has important implications for discussions on ice about plasticity descriptions, rate-controlling processes which accommodate basal glide, and anisotropic ice flow descriptions of large ice masses with the wider perspective of sea level evolution.

Recent changes in north-west Greenland climate documented by NEEM shallow ice core data and simulations, and implications for past-temperature reconstructions

Abstract. Combined records of snow accumulation rate, δ18O and deuterium excess were produced from several shallow ice cores and snow pits at NEEM (North Greenland Eemian Ice Drilling), covering the period from 1724 to 2007. They are used to investigate recent climate variability and characterise the isotope–temperature relationship. We find that NEEM records are only weakly affected by inter-annual changes in the North Atlantic Oscillation. Decadal δ18O and accumulation variability is related to North Atlantic sea surface temperature and is enhanced at the beginning of the 19th century. No long-term trend is observed in the accumulation record. By contrast, NEEM δ18O shows multidecadal increasing trends in the late 19th century and since the 1980s. The strongest annual positive δ18O values are recorded at NEEM in 1928 and 2010, while maximum accumulation occurs in 1933. The last decade is the most enriched in δ18O (warmest), while the 11-year periods with the strongest depletion (coldest) are depicted at NEEM in 1815–1825 and 1836–1846, which are also the driest 11-year periods. The NEEM accumulation and δ18O records are strongly correlated with outputs from atmospheric models, nudged to atmospheric reanalyses. Best performance is observed for ERA reanalyses. Gridded temperature reconstructions, instrumental data and model outputs at NEEM are used to estimate the multidecadal accumulation–temperature and δ18O–temperature relationships for the strong warming period in 1979–2007. The accumulation sensitivity to temperature is estimated at 11 ± 2 % °C−1 and the δ18O–temperature slope at 1.1 ± 0.2 ‰ °C−1, about twice as large as previously used to estimate last interglacial temperature change from the bottom part of the NEEM deep ice core.

Extraction and parametrization of grain boundary networks in glacier ice, using a dedicated method of automatic image analysis

Microstructure analysis of polar ice cores is vital to understand the processes controlling the flow of polar ice on the microscale. This paper presents an automatic image processing framework for extraction and parametrization of grain boundary networks from images of the NEEM deep ice core. As cross-section images are acquired using controlled surface sublimation, grain boundaries and air inclusions appear dark, whereas the inside of grains appears grey. The initial segmentation step of the software is to separate possible boundaries of grains and air inclusions from background. A Machine learning approach is utilized to gain automatic, reliable classification, which is required for processing large data sets along deep ice cores. The second step is to compose the perimeter of section profiles of grains by planar sections of the grain surface between triple points. Ultimately, grain areas, grain boundaries and triple junctions of the later are diversely parametrized. High resolution is achieved, so that small grain sizes and local curvatures of grain boundaries can systematically be investigated.

Water isotopic ratios from a continuously melted ice core sample

Abstract. A new technique for on-line high resolution isotopic analysis of liquid water, tailored for ice core studies is presented. We built an interface between a Wavelength Scanned Cavity Ring Down Spectrometer (WS-CRDS) purchased from Picarro Inc. and a Continuous Flow Analysis (CFA) system. The system offers the possibility to perform simultaneuous water isotopic analysis of δ18O and δD on a continuous stream of liquid water as generated from a continuously melted ice rod. Injection of sub μl amounts of liquid water is achieved by pumping sample through a fused silica capillary and instantaneously vaporizing it with 100% efficiency in a~home made oven at a temperature of 170 °C. A calibration procedure allows for proper reporting of the data on the VSMOW–SLAP scale. We apply the necessary corrections based on the assessed performance of the system regarding instrumental drifts and dependance on the water concentration in the optical cavity. The melt rates are monitored in order to assign a depth scale to the measured isotopic profiles. Application of spectral methods yields the combined uncertainty of the system at below 0.1‰ and 0.5‰ for δ18O and δD, respectively. This performance is comparable to that achieved with mass spectrometry. Dispersion of the sample in the transfer lines limits the temporal resolution of the technique. In this work we investigate and assess these dispersion effects. By using an optimal filtering method we show how the measured profiles can be corrected for the smoothing effects resulting from the sample dispersion. Considering the significant advantages the technique offers, i.e. simultaneuous measurement of δ18O and δD, potentially in combination with chemical components that are traditionally measured on CFA systems, notable reduction on analysis time and power consumption, we consider it as an alternative to traditional isotope ratio mass spectrometry with the possibility to be deployed for field ice core studies. We present data acquired in the field during the 2010 season as part of the NEEM deep ice core drilling project in North Greenland.

Understanding the climatic signal in the water stable isotope records from the NEEM shallow firn/ice cores in northwest Greenland

[1] Samples of precipitation and atmospheric water vapor were collected together with shallow firn/ice cores as part of the new deep drilling project in northwest Greenland: the NEEM project. These samples were analyzed for their isotope composition to understand the processes affecting the climatic signal archived in the water stable isotope records from the NEEM deep ice core. The dominant moisture source for the snow deposited at the NEEM-site may be originating as far south as 35°N from the western part of the Atlantic Ocean. The surface atmospheric water vapor appears in isotopic equilibrium with the snow surface indicating a large water exchange between the atmosphere and snowpack. The interannual variability of NEEM shallow firn/ice cores stable isotope data covering the last ∼40 years shows an unexpectedly weak NAO signal. Regional to global atmospheric models simulate a dominant summer precipitation in the NEEM area, suggesting that the intermittency of modern winter precipitation is responsible for the lack of a strong NAO imprint. The interannual variability of NEEM isotope data however shows a strong correlation with interannual variations of Baffin Bay sea ice cover, a relationship consistent with air mass trajectories. NEEM deep ice core isotopic records may therefore provide detailed information on past Baffin Bay sea ice extent. NEEM stable water isotope content increasing trend points to a local warming trend of ∼3.0°C over the last 40 years.

Seasonal variations of snow chemistry at NEEM, Greenland

AbstarctWe conducted a pit study in July 2009 at the NEEM (North Greenland Eemian Ice Drilling) deep ice-coring site in northwest Greenland. To examine the seasonal variations of snow chemistry and characteristics of the drill site, we collected snow/firn samples from the wall of a 2 m deep pit at intervals of 0.03 m and analyzed them for electric conductivity, pH, Cl–, NO3–, SO42–, CH3SO3– (MSA), Na+, K+, Mg2+, Ca2+ and stable isotopes of water (δ18O and δD). Pronounced seasonal variations in the stable isotopes of water were observed, which indicated that the snow had accumulated regularly during the past 4 years. Concentrations of Na+, Cl– and Mg2+, which largely originate from sea salt, peaked in winter to early spring, while Ca2+, which mainly originates from mineral dust, peaked in late winter to spring, slightly later than Na+, Cl– and Mg2+. Concentrations of NO3– showed double peaks, one in summer and the other in winter to spring, whereas those of SO42– peaked in winter to spring. The winter-to-spring concentrations of NO3– and SO42– seem to have been strongly influenced by anthropogenic inputs. Concentrations of MSA showed double peaks, one in spring and the other in late summer to autumn. Our study confirms that the NEEM deep ice core can be absolutely dated to a certain depth by counting annual layers, using the seasonal variations of stable isotopes of water and those of ions. We calculated the annual surface mass balance for the years 2006–08. The mean annual balance was 176 mm w.e., and the balances for winter-to-summer and summer-to-winter halves of the year were 98 and 78 mm, respectively. Snow deposition during the winter-to-summer half of the year was greater than that during the summer-to-winter half by 10–20mm for all three years covered by this study.