NorthGRIP Ice Core Research Articles

Late Quaternary (MIS 5a-5b) climate and environments of western Iberia inferred from the small-mammal assemblage of Gruta da Oliveira, Torres Novas, Portugal

Marine Isotope Stage 5 (MIS 5) is well represented in palynological studies of North Atlantic marine cores but in only a handful of archeological sequences from the adjacent Iberian landmass. In this paper, we undertake a multi-proxy analysis of small mammal assemblages (insectivores, bats, and rodents) from a 6 m-thick succession, dated to MIS 5a-5b (Layers 7–22), situated in the southeastern edge of the Central Limestone massif of Estremadura at Gruta da Oliveira, near Torres Novas, Portugal. Application of quantitative (Habitat Weighting, Simpson's Diversity Index) and qualitative (Mutual Ecogeographic Range, Bioclimatic Model, Quantified Ecology) methods suggest that open woodland habitats were dominant through the time of accumulation with reconstructed mean annual temperatures (MAT) tracking the 18O curve from the NorthGRIP ice core, from the end of MIS 5c to the beginning of MIS 4. Our findings are consistent with the general trends derived from other Gruta da Oliveira proxies (wood charcoal, coprolite pollen, large mammal associations), limited evidence from other MIS 5 terrestrial sites in Iberia, and offshore marine palynological records.

Initial Mediterranean response to major climate reorganization during the last interglacial-glacial transition

Millennial-scale Dansgaard Oeschger (DO) variability at northern high latitudes has influenced climatic and environmental conditions in the Mediterranean during the last glacial period. There is evidence that the hemispheric transmission of the DO variability occurred at the end of DO event 25; however, the exact timing and the trigger that activated the environmental response in the Mediterranean remains incompletely understood. Here, we provide evidence that the clear millennial-scale teleconnection between Greenland and the Mediterranean started at ∼111.4 ka BP and was initiated by a sub-millennial scale cooling in Greenland (GI-25b). High-resolution sediment proxies and the pollen record of Lago Grande di Monticchio (MON), Italy, reflect climatic instability during the last millennium of the last interglacial, which was characterised by a first and short cooling episode (MON 1) at 111.44 ± 0.69 ka BP, coinciding with the Greenland cold sub-event GI-25b in duration and timing (within dating uncertainties). MON and Greenland (NorthGRIP ice core) also agree in recording a subsequent warm rebound phase that abruptly culminated in the stadial MON 2/GS-25, marking the transition into the last glacial period. Our results show that the GI-25b triggered an early environmental response at MON to centennial-scale climate change in Greenland as a prelude to the millennial-scale teleconnection that was maintained during the glacial period.

Abrupt or gradual? Change point analysis of the late Pleistocene–Holocene climate record from Chew Bahir, southern Ethiopia

AbstractWe used a change point analysis on a late Pleistocene–Holocene lake-sediment record from the Chew Bahir basin in the southern Ethiopian Rift to determine the amplitude and duration of past climate transitions. The most dramatic changes occurred over 240 yr (from ~15,700 to 15,460 yr) during the onset of the African Humid Period (AHP), and over 990 yr (from ~4875 to 3885 yr) during its protracted termination. The AHP was interrupted by a distinct dry period coinciding with the high-latitude Younger Dryas stadial, which had an abrupt onset (less than ~100 yr) at ~13,260 yr and lasted until ~11,730 yr. Wet-dry-wet transitions prior to the AHP may reflect the high-latitude Dansgaard-Oeschger cycles, as indicated by cross-correlation of the potassium record with the NorthGRIP ice core record between ~45–20 ka. These findings may contribute to the debates regarding the amplitude, and duration and mechanisms of past climate transitions, and their possible influence on the development of early modern human cultures.

Predicting the Geothermal Heat Flux in Greenland: A Machine Learning Approach

AbstractGeothermal heat flux (GHF) is a crucial boundary condition for making accurate predictions of ice sheet mass loss, yet it is poorly known in Greenland due to inaccessibility of the bedrock. Here we use a machine learning algorithm on a large collection of relevant geologic features and global GHF measurements and produce a GHF map of Greenland that we argue is within ∼15% accuracy. The main features of our predicted GHF map include a large region with high GHF in central‐north Greenland surrounding the NorthGRIP ice core site, and hot spots in the Jakobshavn Isbræ catchment, upstream of Petermann Gletscher, and near the terminus of Nioghalvfjerdsfjorden glacier. Our model also captures the trajectory of Greenland movement over the Icelandic plume by predicting a stripe of elevated GHF in central‐east Greenland. Finally, we show that our model can produce substantially more accurate predictions if additional measurements of GHF in Greenland are provided.

Implementation of counted layers for coherent ice core chronology

Abstract. A recent coherent chronology has been built for four Antarctic ice cores and the NorthGRIP (NGRIP) Greenland ice core (Antarctic Ice Core Chronology 2012, AICC2012) using a Bayesian approach for ice core dating (Datice). When building the AICC2012 chronology, and in order to prevent any confusion with official ice core chronology, the AICC2012 chronology for NGRIP was forced to fit exactly the GICC05 chronology based on layer counting. However, such a strong tuning did not satisfy the hypothesis of independence of background parameters and observations for the NGRIP core, as required by Datice. We present here the implementation in Datice of a new type of markers that is better suited for constraints deduced from layer counting: the duration constraints. Estimating the global error on chronology due to such markers is not straightforward and implies some assumption on the correlation between individual counting errors for each interval of duration. We validate this new methodological implementation by conducting twin experiments and a posteriori diagnostics on the NGRIP ice core. Several sensitivity tests on marker sampling and correlation between counting errors were performed to provide some guidelines when using such a method for future dating experiments. Finally, using these markers for NGRIP in a five-core dating exercise with Datice leads to new chronologies that do not differ by more than 410 years from AICC2012 for Antarctic ice cores and 150 years from GICC05 for NGRIP over the last 60 000 years.

Water isotope diffusion rates from the NorthGRIP ice core for the last 16,000 years – Glaciological and paleoclimatic implications

A high resolution (0.05 m) water isotopic record (δO18) is available from the NorthGRIP ice core. In this study we look into the water isotope diffusion history as estimated by the spectral characteristics of the δO18 time series covering the last 16,000 years. The diffusion of water vapor in the porous medium of the firn pack attenuates the initial isotopic signal, predominantly having an impact on the high frequency components of the power spectrum. Higher temperatures induce higher rates of smoothing and thus the signal can be used as a firn paleothermometer. We use a water isotope diffusion model coupled to a steady-state densification model in order to infer the temperature signal from the site, assuming the accumulation and strain rate history as estimated using the GICC05 layer counted chronology and a Dansgaard–Johnsen ice flow model. The temperature reconstruction accurately captures the timing and magnitude of the Bølling–Allerød and Younger Dryas transitions. A Holocene climatic optimum is seen between 7 and 9 ky b2k with an apparent cooling trend thereafter. Our temperature estimate for the Holocene climatic optimum, points to a necessary adjustment of the ice thinning function indicating that the ice flow model overestimates past accumulation rates by about 10% at 8 ky b2k. This result, is supported by recent gas isotopic fractionation studies proposing a similar reduction for glacial conditions. Finally, the record presents a climatic variability over the Holocene spanning millennial and centennial scales with a profound cooling occurring at approximately 4000 years b2k. The new reconstruction technique is able to provide past temperature estimates by overcoming the issues apparent in the use of the classical δO18 slope method. It can at the same time resolve temperature signals at low and high frequencies.

Duration of Greenland Stadial 22 and ice-gas Δage from counting of annual layers in Greenland NGRIP ice core

Abstract. High-resolution measurements of chemical impurities and methane concentrations in Greenland ice core samples from the early glacial period allow the extension of annual-layer counted chronologies and the improvement of gas age-ice age difference (Δage) essential to the synchronization of ice core records. We report high-resolution measurements of a 50 m section of the NorthGRIP ice core and corresponding annual layer thicknesses in order to constrain the duration of the Greenland Stadial 22 (GS-22) between Greenland Interstadials (GIs) 21 and 22, for which inconsistent durations and ages have been reported from Greenland and Antarctic ice core records as well as European speleothems. Depending on the chronology used, GS-22 occurred between approximately 89 (end of GI-22) and 83 kyr b2k (onset of GI-21). From annual layer counting, we find that GS-22 lasted between 2696 and 3092 years and was followed by a GI-21 pre-cursor event lasting between 331 and 369 yr. Our layer-based counting agrees with the duration of stadial 22 as determined from the NALPS speleothem record (3250 ± 526 yr) but not with that of the GICC05modelext chronology (2620 yr) or an alternative chronology based on gas-marker synchronization to EPICA Dronning Maud Land ice core. These results show that GICC05modelext overestimates accumulation and/or underestimates thinning in this early part of the last glacial period. We also revise the possible ranges of NorthGRIP Δdepth (5.49 to 5.85 m) and Δage (498 to 601 yr) at the warming onset of GI-21 as well as the Δage range at the onset of the GI-21 precursor warming (523 to 654 yr), observing that temperature (represented by the δ15N proxy) increases before CH4 concentration by no more than a few decades.

Magnetization of polar ice: a measurement of terrestrial dust and extraterrestrial fallout

Laboratory-induced remanent magnetization of polar ice constitutes a measurement of the magnetization carried by the ferromagnetic dust particles in the ice. This non-destructive technique provides a novel kind of information on the dust deposited on the surface of polar ice sheets. Measurements made on ice samples from Greenland (North GRIP ice core) and Antarctica (Vostok and EPICA-Dome C ice cores) allowed the recognition of a fraction of magnetic minerals in ice whose concentration and magnetic properties are directly related to that of insoluble dust. The source of this fraction of magnetic minerals thus appears closely related to terrestrial dust transport and deposition and its magnetic properties are informative of the dust provenance areas. The rock-magnetic properties of the dust may reflect distinct changes of dust source areas from glacial to interglacial periods in agreement with and adding further information to the isotopic (87Sr/86Sr and 143Nd/144Nd) analyses. A second magnetic fraction consists of particles of nanometric size, which are superparamagnetic at freezer temperature and whose concentration is independent of the mass of aerosol dust found in the ice. The source of these nanometric-sized magnetic particles is ascribed to fallout of “meteoric smoke” and their concentration in ice was found to be compatible with the extraterrestrial fallout inferred from Ir concentrations. The diameter of the smoke particles as inferred from magnetic measurements is in the range of about 7–20 nm.

Past surface temperatures at the NorthGRIP drill site from the difference in firn diffusion of water isotopes

Abstract. A new ice core paleothermometer is introduced based on the temperature dependent diffusion of the stable water isotopes in the firn. A new parameter called differential diffusion length is defined as the difference between the diffusion length of the two stable water isotopologues 2H1H16O and 1H218O. A model treatment of the diffusion process of the firn and the ice is presented along with a method of retrieving the diffusion signal from the ice core record of water isotopes using spectral methods. The model shows how the diffusion process is highly dependent on the inter-annual variations in the surface temperatures. It results in a diffusion length longer than if the firn was isothermal. The longer diffusion length can be explained by the strong nonlinearly behaviour of the saturation pressure over ice in the range of the surface temperature fluctuations. The method has been tested on δ18O and δD measurements, spanning the transition from the last glacial to the holocene, from the NorthGRIP ice core. The surface temperature reconstruction based on the differential diffusion resembles other temperature reconstructions for the NorthGRIP ice core. However, the Allerød warming is seen to be significantly warmer than observed in other ice core based temperature reconstructions. The mechanisms behind this behaviour are not fully understood. The method shows the need of an expansion of high resolution stable water isotope datasets from ice cores. However, the new ice core paleothermometer presented here will give valuable insight into past climate, through the physical process of isotope diffusion in the firn column of ice sheets.

Sensitivity of interglacial Greenland temperature and δ<sup>18</sup>O: ice core data, orbital and increased CO<sub>2</sub> climate simulations

Abstract. The sensitivity of interglacial Greenland temperature to orbital and CO2 forcing is investigated using the NorthGRIP ice core data and coupled ocean-atmosphere IPSL-CM4 model simulations. These simulations were conducted in response to different interglacial orbital configurations, and to increased CO2 concentrations. These different forcings cause very distinct simulated seasonal and latitudinal temperature and water cycle changes, limiting the analogies between the last interglacial and future climate. However, the IPSL-CM4 model shows similar magnitudes of Arctic summer warming and climate feedbacks in response to 2 × CO2 and orbital forcing of the last interglacial period (126 000 years ago). The IPSL-CM4 model produces a remarkably linear relationship between TOA incoming summer solar radiation and simulated changes in summer and annual mean central Greenland temperature. This contrasts with the stable isotope record from the Greenland ice cores, showing a multi-millennial lagged response to summer insolation. During the early part of interglacials, the observed lags may be explained by ice sheet-ocean feedbacks linked with changes in ice sheet elevation and the impact of meltwater on ocean circulation, as investigated with sensitivity studies. A quantitative comparison between ice core data and climate simulations requires stability of the stable isotope – temperature relationship to be explored. Atmospheric simulations including water stable isotopes have been conducted with the LMDZiso model under different boundary conditions. This set of simulations allows calculation of a temporal Greenland isotope-temperature slope (0.3–0.4‰ per °C) during warmer-than-present Arctic climates, in response to increased CO2, increased ocean temperature and orbital forcing. This temporal slope appears half as large as the modern spatial gradient and is consistent with other ice core estimates. It may, however, be model-dependent, as indicated by preliminary comparison with other models. This suggests that further simulations and detailed inter-model comparisons are also likely to be of benefit. Comparisons with Greenland ice core stable isotope data reveals that IPSL-CM4/LMDZiso simulations strongly underestimate the amplitude of the ice core signal during the last interglacial, which could reach +8–10 °C at fixed-elevation. While the model-data mismatch may result from missing positive feedbacks (e.g. vegetation), it could also be explained by a reduced elevation of the central Greenland ice sheet surface by 300–400 m.

Radiocarbon wiggle matching of Swedish lake varves reveals asynchronous climate changes around the 8.2 kyr cold event

Snowball, I., Muscheler, R., Zillén, L., Sandgren, P., Stanton, T. & Ljung, K. 2010: Radiocarbon wiggle matching of Swedish lake varves reveals asynchronous climate changes around the 8.2 kyr cold event. Boreas, Vol. 39, pp. 720–733. 10.1111/j.1502‐3885.2010.00167.x. ISSN 0300‐9483.A series of nine radiocarbon measurements obtained at increments of 50 years was used to wiggle match varved sediments in Lake Kälksjön (west central Sweden) to the calendar year time scale provided by the tree‐ring‐derived radiocarbon calibration curve. The anchor point for a series of 400 varves known to span the ‘8.2 kyr cold event’ was derived from a combination of correlation analysis, Monte Carlo statistics and Bayesian age modelling. The GRIP and NorthGRIP ice‐cores were matched to the same absolute time scale by comparing 10Be data and tree‐ring Δ14C. The radiocarbon‐based wiggle match, organic carbon measurements, mineral magnetic parameters and XRF data reconstruct a distinct period of enhanced erosion in Kälksjön's catchment between 8066 ± 25 and 7920 ± 25 cal. yr BP. Prior to human impact, the increased erosion in the boreal environment is assigned to increased winter precipitation in the form of snow, which caused spring meltwater discharge to intensify. Our results suggest that the abrupt onset of increased winter precipitation in west central Sweden started at least 50 years after the onset of the ‘8.2 kyr cold event’ as defined by oxygen isotope data from Greenland. The study highlights the value of synchronized time scales in the reconstruction of abrupt climate changes based on proxy data.

Millennial and sub-millennial scale climatic variations recorded in polar ice cores over the last glacial period

Abstract. Since its discovery in Greenland ice cores, the millennial scale climatic variability of the last glacial period has been increasingly documented at all latitudes with studies focusing mainly on Marine Isotopic Stage 3 (MIS 3; 28–60 thousand of years before present, hereafter ka) and characterized by short Dansgaard-Oeschger (DO) events. Recent and new results obtained on the EPICA and NorthGRIP ice cores now precisely describe the rapid variations of Antarctic and Greenland temperature during MIS 5 (73.5–123 ka), a time period corresponding to relatively high sea level. The results display a succession of abrupt events associated with long Greenland InterStadial phases (GIS) enabling us to highlight a sub-millennial scale climatic variability depicted by (i) short-lived and abrupt warming events preceding some GIS (precursor-type events) and (ii) abrupt warming events at the end of some GIS (rebound-type events). The occurrence of these sub-millennial scale events is suggested to be driven by the insolation at high northern latitudes together with the internal forcing of ice sheets. Thanks to a recent NorthGRIP-EPICA Dronning Maud Land (EDML) common timescale over MIS 5, the bipolar sequence of climatic events can be established at millennial to sub-millennial timescale. This shows that for extraordinary long stadial durations the accompanying Antarctic warming amplitude cannot be described by a simple linear relationship between the two as expected from the bipolar seesaw concept. We also show that when ice sheets are extensive, Antarctica does not necessarily warm during the whole GS as the thermal bipolar seesaw model would predict, questioning the Greenland ice core temperature records as a proxy for AMOC changes throughout the glacial period.

Synchronising EDML and NorthGRIP ice cores using δ 18O of atmospheric oxygen (δ 18O atm) and CH 4 measurements over MIS5 (80–123 kyr)

Water isotope records from the EPICA Dronning Maud Land (EDML) and the NorthGRIP ice cores have revealed a one to one coupling between Antarctic Isotope Maxima (AIM) and Greenland Dansgaard-Oeschger (DO) events back to 50 kyr. In order to explore if this north–south coupling is persistent over Marine Isotopic Stage 5 (MIS 5), a common timescale must first be constructed. Here, we present new records of δ 18O of O 2 (δ 18O atm) and methane (CH 4) measured in the air trapped in ice from the EDML (68–147 kyr) and NorthGRIP (70–123 kyr) ice cores. We demonstrate that, through the period of interest, CH 4 records alone are not sufficient to construct a common gas timescale between the two cores. Millennial-scale variations of δ 18O atm are evidenced over MIS 5 both on the Antarctic and Greenland ice cores and are coupled to CH 4 profiles to synchronise the NorthGRIP and EDML records. They are shown to be a precious tool for ice core synchronisation. With this new dating strategy, we produce the first continuous and accurate sequence of the north–south climatic dynamics on a common ice timescale for the last glacial inception and the first DO events of MIS 5, reducing relative dating uncertainties to an accuracy of a few centuries at the onset of DO events 24 to 20. This EDML-NorthGRIP synchronisation provides new firm evidence that the bipolar seesaw is a pervasive pattern from the beginning of the glacial period. The relationship between Antarctic warming amplitudes and their concurrent Greenland stadial duration highlights the particularity of DO event 21 and its Antarctic counterpart. Our results suggest a smaller Southern Ocean warming rate for this long DO event compared to DO events of MIS 3.

Marine glacial and interglacial stratigraphy in Vendsyssel, northern Denmark: foraminifera and stable isotopes

The marine Quaternary of Vendsyssel has been studied in a series of new boreholes in the area, and the climatic development is discussed on the basis of foraminiferal assemblages and stable isotopes. The foraminiferal zones are correlated with previously published records from northern Denmark, and the spatial local and regional distribution is discussed in details based on the new evidence. The new data show that the marine sedimentation in Vendsyssel was not continuous from the Late Saalian to the Middle Weichselian, as previously thought. For example, there is indication of a hiatus at our key site, Åsted Vest in the central part of Vendsyssel, at the transition between regional foraminiferal zones N4 and N3, i.e. at the Late Saalian (MIS 6) – Eemian (MIS 5e) transition. The hitherto most complete Early Weichselian succession (zone N2) in Vendsyssel is presented from Åsted Vest. Deposits from the Early Weichselian sea‐level lowstands (MIS 5d and 5b) may, however, be missing in parts of the area. Two major breaks in the marine deposition during the Middle Weichselian represent glacial advances into northern Denmark. The first event occurred just after deposition of the regional foraminiferal zone N2 (late MIS 4), and the second event in the middle part of zone N1 (early MIS 3). Zone N1 is succeeded by a series of non‐marine units deposited during the sea‐level lowstand of the Weichselian maximum glaciation (late MIS 3 and MIS 2), including deeply incised tunnel valleys, which have been refilled with non‐marine sediments during the Late Weichselian. Vendsyssel was inundated by the sea again during the Late Weichselian, at c. 18 kyr BP. Subsequently, the marine conditions were gradually changed by forced regression caused by local isostatic uplift, and around the Weichselian–Holocene transition most of Vendsyssel was above sea level. A continuous deposition across the Late Weichselian–Holocene boundary only occurred at relatively deep sites such as Skagen. The environmental and climatic indications for Vendsyssel are in accordance with the global sea‐level curve, and the Quaternary record is correlated with the oxygen isotope record from the NorthGRIP ice core, as well as the marine isotope stages.

Formal definition and dating of the GSSP (Global Stratotype Section and Point) for the base of the Holocene using the Greenland NGRIP ice core, and selected auxiliary records

AbstractThe Greenland ice core from NorthGRIP (NGRIP) contains a proxy climate record across the Pleistocene–Holocene boundary of unprecedented clarity and resolution. Analysis of an array of physical and chemical parameters within the ice enables the base of the Holocene, as reflected in the first signs of climatic warming at the end of the Younger Dryas/Greenland Stadial 1 cold phase, to be located with a high degree of precision. This climatic event is most clearly reflected in an abrupt shift in deuterium excess values, accompanied by more gradual changes in δ18O, dust concentration, a range of chemical species, and annual layer thickness. A timescale based on multi‐parameter annual layer counting provides an age of 11 700 calendar yr b2 k (before AD 2000) for the base of the Holocene, with a maximum counting error of 99 yr. A proposal that an archived core from this unique sequence should constitute the Global Stratotype Section and Point (GSSP) for the base of the Holocene Series/Epoch (Quaternary System/Period) has been ratified by the International Union of Geological Sciences. Five auxiliary stratotypes for the Pleistocene–Holocene boundary have also been recognised. Copyright © 2008 John Wiley & Sons, Ltd.

The Global Stratotype Section and Point (GSSP) for the base of the Holocene Series/Epoch (Quaternary System/Period) in the NGRIP ice core

The Greenland ice core from NorthGRIP (NGRIP) contains a proxy climate record across the Pleistocene-Holocene boundary of unprecedented clarity and resolution. Analysis of an array of physical and chemical parameters within the ice enables the base of the Holocene, as reflected in the first signs of climatic warming at the end of the Younger Dryas/Greenland Stadial 1 cold phase, to be located with a high degree of precision. This climatic event is most clearly reflected in an. abrupt shaft in deuterium excess values, accompanied by more gradual changes in delta O-18, dust concentration, a range of chemical species, and annual layer thickness. A timescale based on multi-parameter annual layer counting provides an age of 11,700 yr b2k (before AD2000) for the base of the Holocene, with, an estimated 2 sigma uncertainty of 99 yr: It is proposed that an archived core from this unique sequence should constitute the Global Stratotype Section and Point (GSSP) for the base of the Holocene Series/Epoch (Quaternary System/Period).

A 60 000 year Greenland stratigraphic ice core chronology

Abstract. The Greenland Ice Core Chronology 2005 (GICC05) is a time scale based on annual layer counting of high-resolution records from Greenland ice cores. Whereas the Holocene part of the time scale is based on various records from the DYE-3, the GRIP, and the NorthGRIP ice cores, the glacial part is solely based on NorthGRIP records. Here we present an 18 ka extension of the time scale such that GICC05 continuously covers the past 60 ka. The new section of the time scale places the onset of Greenland Interstadial 12 (GI-12) at 46.9±1.0 ka b2k (before year AD 2000), the North Atlantic Ash Zone II layer in GI-15 at 55.4±1.2 ka b2k, and the onset of GI-17 at 59.4±1.3 ka b2k. The error estimates are derived from the accumulated number of uncertain annual layers. In the 40–60 ka interval, the new time scale has a discrepancy with the Meese-Sowers GISP2 time scale of up to 2.4 ka. Assuming that the Greenland climatic events are synchronous with those seen in the Chinese Hulu Cave speleothem record, GICC05 compares well to the time scale of that record with absolute age differences of less than 800 years throughout the 60 ka period. The new time scale is generally in close agreement with other independently dated records and reference horizons, such as the Laschamp geomagnetic excursion, the French Villars Cave and the Austrian Kleegruben Cave speleothem records, suggesting high accuracy of both event durations and absolute age estimates.

Relation between neighbouring grains in the upper part of the NorthGRIP ice core — Implications for rotation recrystallization

We propose a new method to investigate the relationships between neighbouring crystals and apply it to the textures measured along the upper 900 m of the NorthGRIP ice core. This method shows unexpected correlations between neighbours in the so-called normal grain growth regime, questioning the classical view on the onset of rotation recrystallization in ice-sheets. Moreover, the fractionation rate associated to the rotation recrystallization appears constant through time. Finally, grains with low-angle boundaries do not present a special distribution pattern of their c-axes. This suggests that rotation recrystallization is an isotropic process not affected by the direction of the macroscopic strain.

Estimating the basal melt rate at NorthGRIP using a Monte Carlo technique

AbstractFrom radio-echo sounding (RES) surveys and ice core data it can be seen that the ice sheet is melting at the base in a large area in Northern Greenland. The RES images reveal internal layers in the ice. The layers are former deposition surfaces and are thus isochrones. Undulations of the isochrones in regions where the base is smooth suggest that the basal melt rate changes over short distances. This indicates that the geothermal heat flux is very high and has large spatial variability in Northern Greenland. In this study, the basal melt rate at the NorthGRIP drill site in North-Central Greenland is calculated by inverse modelling. We use simple one- and two-dimensional flow models to simulate the ice flow along the NNW-trending ice ridge leading to NorthGRIP. The accumulation is calculated from a dynamical model. Several ice flow parameters are unknown and must be estimated along with the basal melt rate using a Monte Carlo method. The Monte Carlo inversion is constrained by the observed isochrones, dated from the timescale established for the NorthGRIP ice core. The estimates of the basal melt rates around NorthGRIP are obtained from both the one- and two-dimensional models. Combining the estimated basal melt rates with the observed borehole temperatures allows us to convert the basal melt rates to geothermal heat flow values. From the two-dimensional model we find the basal melt rate and geothermal heat flux at NorthGRIP to be 6.1 mma–1 and 129 mWm–2, respectively.

The Greenland Ice Core Chronology 2005, 15–42 ka. Part 2: comparison to other records

A new Greenland Ice Core Chronology (GICC05) based on multi-parameter counting of annual layers has been obtained for the last 42 ka. Here we compare the glacial part of the new time scale, which is based entirely on records from the NorthGRIP ice core, to existing time scales and reference horizons covering the same period. These include the GRIP and NorthGRIP modelled time scales, the Meese-Sowers GISP2 counted time scale, the Shackleton–Fairbanks GRIP time scale (SFCP04) based on 14C calibration of a marine core, the Hulu Cave record, three volcanic reference horizons, and the Laschamp geomagnetic excursion event occurring around Greenland Interstadial 10. GICC05 is generally in good long-term agreement with the existing Greenland ice core chronologies and with the Hulu Cave record, but on shorter time scales there are significant discrepancies. Around the Last Glacial Maximum there is a more than 1 ka age difference between GICC05 and SFCP04 and a more than 0.5 ka discrepancy in the same direction between GICC05 and the age of a recently identified tephra layer in the NorthGRIP ice core. Both SFCP04 and the tephra age are based on 14C-dated marine cores and fixed marine reservoir ages. For the Laschamp event, GICC05 agrees with a recent independent dating within the uncertainties.