# A 120,000-year long climate record from a NW-Greenland deep ice core at ultra-high resolution

Vasileios Gkinis, Bo Vinther, Trevor Popp, Thea Quistgaard, Anne-Katrine Faber, Christian Holme, Camilla-Marie Jensen, Mika Lanzky, Anine-Maria Lütt, Vasileios Mandrakis, et al.

► To cite this version:

Vasileios Gkinis, Bo Vinther, Trevor Popp, Thea Quistgaard, Anne-Katrine Faber, et al. A 120,000-year long climate record from a NW-Greenland deep ice core at ultra-high resolution. Scientific Data, 2021, 8, pp.141. 10.1038/s41597-021-00916-9. hal-03238378

HAL Id: hal-03238378

https://hal.science/hal-03238378

Submitted on 27 May 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
scientific **data**

<figure><img src="image_2.png" /></figure>

Check for updates

OPEN

DATA DESCRIPTOR

A 120,000-year long climate record
from a NW-Greenland deep ice core
at ultra-high resolution

Vasileios Gkinis<fnref n="1" />, Bo M. Vinther<fnref n="1" />, Trevor J. Popp<fnref n="1" />, Thea Quistgaard<fnref n="1" />,
Anne-Katrine Faber<fnref n="1" />,², Christian T. Holme<fnref n="1" />, Camilla-Marie Jensen<fnref n="1" />,³, Mika Lanzky<fnref n="1" />,⁴,
Anine-Maria Lütt<fnref n="1" />, Vasileios Mandrakis<fnref n="1" />,³, Niels-Ole Ørum<fnref n="1" />, Anna-Sofie Pedersen<fnref n="1" />,
Nikol Vaxevani<fnref n="1" />, Yongbiao Weng<fnref n="1" />,², Emilie Capron<fnref n="1" />,⁵, Dorthe Dahl-Jensen<fnref n="1" />,⁶, Maria Hörhold<fnref n="1" />,⁷,
Tyler R. Jones<fnref n="1" />,⁸, Jean Jouzel⁹, Amaëlle Landais⁹, Valérie Masson-Delmotte⁹, Hans Oerter⁷,
Sune O. Rasmussen<fnref n="1" />, Hans Christian Steen-Larsen<fnref n="1" />,⁵, Jørgen-Peder Steffensen<fnref n="1" />,⁵,
Árný-Erla Sveinbjörnsdóttir¹⁰, Anders Svensson<fnref n="1" />, Bruce Vaughn⁸ & James W. C. White⁸

We report high resolution measurements of the stable isotope ratios of ancient ice (δ¹⁸O, δD) from the North Greenland Eemian deep ice core (NEEM, 77.45° N, 51.06° E). The record covers the period 8–130 ky b2k (y before 2000) with a temporal resolution of ≈0.5 and 7 y at the top and the bottom of the core respectively and contains important climate events such as the 8.2 ky event, the last glacial termination and a series of glacial stadials and interstadials. At its bottom part the record contains ice from the Eemian interglacial. Isotope ratios are calibrated on the SMOW/SLAP scale and reported on the GICC05 (Greenland Ice Core Chronology 2005) and AICC2012 (Antarctic Ice Core Chronology 2012) time scales interpolated accordingly. We also provide estimates for measurement precision and accuracy for both δ¹⁸O and δD.

**Background & Summary**

The isotopic composition of ice from deep cores drilled on Greenland and Antarctica has traditionally been used as a proxy for past temperatures<fnref n="1" />⁻⁴ offering a picture of past climate that extends as far as 120,000 y in the past⁵ for the case of Greenland and 800,000 y for the case of Antarctica⁶. Typically expressed with the δ notation (isotopic abundances are reported as deviations of a sample's isotopic ratio relative to that of a reference water (e.g. VSMOW) expressed in per mille (‰) as: $δ^i = \left( \frac{^{iR}_{sample} - 1}{^{iR}_{VSMOW}} \right) \times 1000$ where $^{2R} = \frac{2H}{1H}$ and $^{18R} = \frac{^{18O}}{^{16O}}$), the isotopic composition records hydrological cycle changes spanning centennial to millennial scales. At sufficiently high resolution, water isotope records can be used in order to resolve climatic signals at higher frequencies, investigate abrupt climate events, count annual layers for the purpose of ice core chronologies and extract paleoclimate information stored in their spectral signature⁷⁻¹⁴.

Stable water isotope analysis has traditionally been performed with Isotope Ratio Mass Spectrometry (IRMS) utilising an array of different implementations, each one presenting advantages and disadvantages with respect to precision, accuracy, sample consumption and throughput, as well as ease of use and labour intensity¹⁵⁻¹⁹. With the advent of fast electronics and room-temperature, turn-key, solid state laser sources, there emerged various optical

<footnote n="1">¹Physics of Ice Climate and Earth, Niels Bohr Institute, University of Copenhagen, Tagensvej 16, 2200, Copenhagen, Denmark. ²Geophysical Institute and Bjerknes Centre for Climate Research, University of Bergen, Allégaten 70, 5007, Bergen, Norway. ³Climate and Environmental Physics, Physics Institute, Sidlerdtrasse 5, 3012, Bern, Switzerland. ⁴Department of Geosciences, University of Oslo, Sem Sæ lands Vei 1, 0371, Oslo, Norway. ⁵Institut des Géosciences de l'Environnement, Université Grenoble Alpes, CNRS, IRD, G-INP, 38000, Grenoble, France. ⁶Centre for Earth Observation Science, 535 Wallace Building, University of Manitoba, Winnipeg, MB, R3T 2N2, Canada. ⁷Alfred-Wegener-Institut Helmholtz-Zentrum für Polar- und Meeresforschung, Am Handelshafen 12, 27570, Bremerhaven, Germany. ⁸Institute of Arctic and Alpine Research, University of Colorado, Boulder, Colorado, 80309-0450, USA. ⁹Laboratoire des Sciences du Climat et de l'Environnement (LSCE), Institut Pierre Simon Laplace (CEA-CNRS-UVSQ UMR 8212), Université Paris Saclay, Gif-sur-Yvette, France. ¹⁰Institute of Earth Sciences, University of Iceland, Sturlugata 7, Reykjavik, Iceland. e-mail: v.gkinis@nbi.ku.dk</footnote>
<figure><img src="image_3.png" /><figcaption>**Fig. 1** Sampling scheme for the NEEM ice core. CFA refers to the section sampled for Continuous Flow Analysis and Iso refers to the section sampled for water isotope analysis.</figcaption></figure>

spectroscopic techniques, typically in the near Infra-Red spectral region that allowed for quasi-simultaneous analysis of the three isotope ratios $^{2}$H/$^{1}$H, $^{17}$O/$^{16}$O and $^{18}$O/$^{16}$O on the same water sample$^{20,21}$. Commercial instruments using Cavity Enhanced Spectroscopic techniques allowed for routine simultaneous measurements of δ$^{18}$O and δD using a greatly reduced amount of sample and with precision and accuracy comparable, if not better than that of IRMS$^{22,23}$. The new techniques offered the possibility for in-situ measurements of the isotopic composition of vapour$^{24–26}$ while in combination with continuous melter systems, and micro-volume flash vaporisers, they yielded ice core data sets of unprecedented resolution$^{23,27–31}$.

Here, we present a high-resolution (0.05 m) water isotope (δ$^{18}$O and δD) record$^{32}$ from the Greenland NEEM ice core, drilled from 2007 to 2012 at 77.45°N–51.06°W. The record covers the segment of the ice core below 1210.5 m - a depth that corresponds to an age of 8000 y b2k (y before 2000 AD)- and extends to a final depth of 2536.5 m. The oldest age of the record is 129,258 y b2k at the depth of 2432.15 m. We report ages on the GICC05 (ages 8–121.446 ky b2k) and the AICC2012 (ages 60–129.258 ky b2k) time scales interpolated accordingly. Beyond that point and for the deepest 104 m, an age scale cannot be reconstructed as the record is too disturbed due to ice folds$^{33}$.

Isotopic analysis has been performed using Laser Spectroscopy and in particular three different versions of the Cavity Ring Down Spectrometers (CRDS) from Picarro Inc. (L2120-i, L2130-i, L2140-i). All measurements are reported on the international VSMOW-SLAP (VSMOW and SLAP refer to the International Atomic Energy reference water materials and stand for Vienna Mean Ocean Water and Standard Light Antarctic Precipitation) isotope scale after careful calibration using a triplet of local standards calibrated against the primary International Atomic Energy Agency reference materials. We describe the methodology of the measurement including information on instrument calibration procedures and an accuracy and precision assessment.

The dataset will be useful for future studies of past climate variability on long (millennial) and short (decadal to annual) time scales. Studies focusing on the spectral properties of the isotopic signal will also benefit from the high resolution apparent in the record. The record includes significant past-climate events such as the 8.2 ky climatic transition, the last Glacial Termination, as well as the series of Glacial Stadial-Interstadials of the last 120 ky (GI-1, GS-1 to GI-25c, GS-26)0$^{34}$.

## Methods

**Discrete Ice Core Sampling.** The sampling of the ice core for water isotope analysis was performed in the field during the field seasons 2009–2012. A wedge with a cross section of approximately 10 cm$^2$ (roughly 12% of the cross section of the ice core) was cut using a band saw with a resolution of 0.025 m for the top 601.7 m and 0.05 m for the interval 601.75–2536.5 (Fig. 1). The dataset presented here, has a resolution of 0.05 m. The ice samples were stored and transported to Copenhagen, Denmark frozen in plastic bags for further preparation and isotopic analysis.

Subsequently, and after arrival to the Stable Isotope Lab at the Niels Bohr Institute, the samples were melted at room temperature in air-tight containers and transferred to 10 ml PE-LD narrow mouth bottles (VITLAB-138093). From that point, the samples were stored frozen until the time of isotopic analysis using CRDS. Prior to analysis, 200 µl of sample were transferred from the PE-LD bottles to 2 ml short thread glass vials with integrated 300 µl micro-inserts (La-Pha-Pack 11092357) sealed with triple layer Teflon/Silicon/Teflon caps (La-Pha-Pack 09150480).

**Stable isotope analysis.** *Cavity ring down spectroscopy.* Isotope analysis is performed with Cavity Ring Down Spectroscopy (CRDS) in the near Infra-Red region. For this study, we have used three different models of
<table><caption>**Table 1.** Types of instruments used and their relative contribution.</caption><thead><tr><td>Instrument</td><td>Points measured</td><td>Coverage [%]</td></tr></thead><tbody><tr><td>L2120-i</td><td>6857</td><td>25.9</td></tr><tr><td>L2130-i</td><td>14384</td><td>54.2</td></tr><tr><td>L2140-i</td><td>5027</td><td>19</td></tr><tr><td>Interpolated</td><td>251</td><td>0.9</td></tr><tr><td>Total</td><td>26519</td><td>100.0</td></tr></tbody></table>

<table><caption>**Table 2.** PAL Autosampler Injection method.</caption><thead><tr><td>Property</td><td>Value</td></tr></thead><tbody><tr><td>Syringe Volume</td><td>10 µl</td></tr><tr><td>Sample Volume</td><td>3.8 µl</td></tr><tr><td>Fill Speed</td><td>400 nl/s</td></tr><tr><td>Fill Strokes</td><td>1</td></tr><tr><td>Pullup Delay</td><td>2.0 s</td></tr><tr><td>Inject Speed</td><td>50 nl/s</td></tr><tr><td>Pre Injection Delay</td><td>0 s</td></tr><tr><td>Post Injection Delay</td><td>10 s</td></tr></tbody></table>

the Picarro L21XX-i analyser namely the L2120-i, L2130-i and L2140-i. The record contains a total of 26,519 data points acquired with the three instruments. The exact number of data points per instrument is given in Table 1.

All three models follow the same spectroscopic principle for the measurement of isotope ratios and are very similar with respect to their performance. As we show in section 4, the L2130-i and L2140-i models reach higher levels of precision and accuracy when compared to the older generation L2120-i instrument. CRDS measurements are performed in the vapour phase with the sample continuously flowing through a high finesse optical cavity. The cavity consists of three high reflectivity mirrors in a V-shape configuration²². A continuous emission diode laser is coupled to the optical cavity in which the temperature and the pressure are accurately controlled (75 Torr and 85 °C). The light leaking out of one of the cavity mirrors is measured by a solid state photodiode. When the light intensity reaches a certain threshold, the laser is abruptly switched off and the light intensity signal on the photodiode decays exponentially. The time constant of the exponential decay, commonly referred to as 'ring-down time' depends on the concentration of the absorbant molecules present in the optical cavity. For this study we have performed measurements of δ¹⁸O and δD, thus the CRDS protocol used, takes into account absorption lines of the ¹H₂¹⁶O, ¹H₂¹⁸O and ¹H²H¹⁶O isotopologues.

*Sample vaporisation and injection method.* Vaporisation of the sample is performed using a high throughput, low volume vaporiser (Picarro-A0212 – discontinued model as of 2016). The A0212 vaporiser utilises a continuous vaporisation approach similar to that developed for the Continuous Flow δ¹⁸O / δD analysis²³. A nano-flow of the liquid sample generated by a micro-volume syringe, is vaporised continuously on a hot stainless steel plate at a temperature of 170 °C within a stream of dry air ([H₂O] < 30 ppm). The CRDS analyser pumps the water vapour sample into the optical cavity with a flow rate of 30 cc/min (STP) via an open split connection with the inlet of the vaporiser.

A total of 3.8 µl of liquid water sample is injected with a liquid flow of 50 nl/s using a 10 µl syringe (SGE 10R-C/F5). The stability of the sample injection is a critical prerequisite for achieving high quality measurements with the continuous vaporisation method. An injection cleaning procedure is performed on a daily basis using an ultrasonic acetic acid bath followed by several manual rinses with deionised water. The thorough daily cleaning is vital for achieving measurements of high quality. The injection is performed using a CTC GC-PAL autosampler (PAL hereafter) with an injection method that results in a pulse duration of ≈80 sec and a pulse amplitude of 20,000 ppm. The amount of time required for one injection, including the filling and rinsing steps of the injection method is ≈135 sec. In Table 2 we present the properties of the injection method and in Fig. 2 we plot the water concentration and the δD signal over a sequence of four injection pulses.

*Tray configuration - PAL Autosampler jobs sequence.* In Table 3 we describe the sequence of the PAL autosampler jobs and the tray configuration. The measurement run begins with the VSMOW-SLAP calibration block (Vials #1 and #3) followed by the unknown (with the term “unknown” here we refer to the ice core samples whose isotopic composition is to be measured) samples run sequentially with increasing depth (Vials #4–56). A measurement of every local water standard consists of 20 injections of which the first 12 are discarded, in order to account for memory effects caused by the large isotopic step between the standards. The measurement of the unknowns is based on a block of four injections of which the first one is discarded. Three blocks containing a “check” standard that is treated as an unknown (Vial #2) are introduced in the middle and at the end of the VSMOW-SLAP calibration as well as at the end of the run after all unknowns have been measured. The “check” standard provides a measure of the obtained accuracy as well as the long term stability of the run. In total, 300
<figure><img src="image_4.png" /><figcaption>Fig. 2 Vaporiser pulse example and run tray configuration. (a) A sequence of four vaporiser pulses focused at the transition from the “−40” to the “NEEM” standard water. (b) Vial tray configuration and raw δD data of a full measurement run as a function of time.</figcaption></figure>

<table><caption>Table 3. Description of the sequence of PAL Autosampler jobs throughout a full measurement run.</caption><thead><tr><th>Job</th><th>From Vial</th><th>To Vial</th><th>Inj per Vial</th><th>Discard</th><th>Description</th></tr></thead><tbody><tr><td>1</td><td>1</td><td>3</td><td>20</td><td>12</td><td>VSMOW/SLAP calibration block</td></tr><tr><td>2</td><td>2</td><td>2</td><td>8</td><td>5</td><td>Standard 2 “check”</td></tr><tr><td>3</td><td>4</td><td>4</td><td>8</td><td>5</td><td>1st unknown</td></tr><tr><td>4</td><td>5</td><td>56</td><td>4</td><td>1</td><td>Unknowns 2–53</td></tr><tr><td>5</td><td>2</td><td>2</td><td>16</td><td>10</td><td>Standard 2 “check”</td></tr></tbody></table>

<table><caption>Table 4. Copenhagen local standards used for the NEEM record. All values are given in ‰ with respect to VSMOW.</caption><thead><tr><th>Name</th><th>δ<sup>18</sup>O</th><th>δD</th></tr></thead><tbody><tr><td>−22</td><td>−21.88</td><td>−168.3</td></tr><tr><td>NEEM</td><td>−33.5</td><td>−257.1</td></tr><tr><td>−40</td><td>−39.93</td><td>−310.7</td></tr><tr><td>DC-02</td><td>−54.07</td><td>−428.2</td></tr></tbody></table>

injections are performed, resulting to a total run duration of approximately 11 hours and 40 minutes. A complete run is illustrated in Fig. 2.

VSMOW-SLAP calibration protocols. The raw isotope measurements are calibrated on the VSMOW-SLAP scale using a two fixed-point calibration similar to²⁷ and following the IAEA recommended procedures. For the record presented here, we have used 4 water isotopic standards calibrated against the primary IAEA reference waters VSMOW and SLAP (Table 4). The standards were combined in two triplets [−22, NEEM, −40] and [NEEM, −40, DC-02]. The first triplet [−22, NEEM, −40] is used for sections of the Holocene, whereas the more isotopically depleted samples during the Glacial are calibrated using the second triplet [NEEM, −40, DC-02].

The two extreme standards of every triplet are used in order to define the slope and the intercept of the calibration line

$$ \delta_{\text{vsmow}} = \alpha \delta_{\text{raw}} + \beta \qquad (1) $$

where the slope α and the intercept β are for the example of the low water standards triplet equal to:

$$ \alpha = \frac{\delta_{\text{NEEMvsmow}} - \delta_{\text{DC02Mvsmow}}}{\delta_{\text{NEEMraw}} - \delta_{\text{DC02raw}}} \qquad (2) $$

$$ \beta = \delta_{\text{NEEMvsmow}} - \frac{\delta_{\text{NEEMvsmow}} - \delta_{\text{DC02vsmow}}}{\delta_{\text{NEEMraw}} - \delta_{\text{DC02raw}}} \delta_{\text{NEEMraw}} \qquad (3) $$

From Eqs. 1–3 we calculate the VSMOW calibrated value of the middle “check” standard, which is then compared to the assigned VSMOW value of the “check” standard water in Table 4.
<table><caption>**Table 5.** Name and description of the variables included in the data file.</caption><thead><tr><td>Column</td><td>Name</td><td>Unit</td><td>Description</td></tr></thead><tbody><tr><td>1</td><td>Depth</td><td>m</td><td>Bottom Depth of each discrete sample</td></tr><tr><td>2</td><td>Age</td><td>y b2k</td><td>Age on the GICC05 counted time scale</td></tr><tr><td>3</td><td>Age</td><td>y b2k</td><td>Age on the GICC05 model extended time scale</td></tr><tr><td>4</td><td>Age</td><td>y b2k</td><td>Age on the AICC2012 time scale</td></tr><tr><td>5</td><td>MCE</td><td>y</td><td>Maximum Counting Error for the GICC05 time scale</td></tr><tr><td>6</td><td>δ<sup>18</sup>O</td><td>%o</td><td>δ<sup>18</sup>O of the sample</td></tr><tr><td>7</td><td>δD</td><td>%o</td><td>δD of the sample</td></tr><tr><td>8</td><td>δ<sup>18</sup>O std dev</td><td>%o</td><td>Run precision estimate for δ<sup>18</sup>O</td></tr><tr><td>9</td><td>δD std dev</td><td>%o</td><td>Run precision estimate for δD</td></tr><tr><td>10</td><td>Offset</td><td>%o</td><td>Run accuracy estimate for δ<sup>18</sup>O</td></tr><tr><td>11</td><td>Offset</td><td>%o</td><td>Run accuracy estimate for δD</td></tr></tbody></table>

**Data analysis.** *Data analysis for individual runs.* Upon completion, individual runs are processed using a spreadsheet-based data analysis scheme. The scheme performs the calculations of the VSMOW-SLAP calibration slope and intercept and subsequently applies them to the raw data, yielding a VSMOW-SLAP calibrated dataset. The data analysis sheet also calculates and reports the uncertainties of each measurement block (water standards and unknown samples). The main quality-control metric of the protocol is the calibrated value of the “check” standard, as well as the standard deviation calculated based on the three “check” standard blocks. The first value gives an estimate of the measurement’s accuracy, while the second is an estimate of the individual run’s noise level. The data analysis scheme selects the appropriate number of injections to be considered as described in Sec. 2 and generates a list of the analysed samples with their unique ID’s and isotopic composition expressed in per mille. Possible outliers in the run can occur in the rare case of faulty sample injections. These outliers are identified by visual inspection of the run results, as well as an inspection of the injection level and stability of the suspected unknowns. These data points are marked and then filled by means of linear interpolation. The record contains a total of 251 interpolated points accounting for 0.9% of the record (Table 1).

*Post processing of the record.* Compilation of the full record based on the individual run's spreadsheets is performed using a Python routine that reads the data of interest sequentially. The Python routine reads the samples ID's and their respective isotopic composition and assigns a section bottom depth to every sample. Additionally, metadata of interest are collected from each individual run and compiled on a separate file. These data concern sample injection characteristics, metrics related to the check standard as well as the values related to the VSMOW-SLAP calibration line.

*The age scale.* The GICC05 (Greenland Ice Core Chronology 2005)³⁵, GICC05modelext³⁶ and AICC2012 (Antarctic Ice Core Chronology 2012)³⁷,³⁸ timescales have all been transferred to NEEM high resolution depths using the approaches of³³,³⁹,⁴⁰.

The GICC05 NEEM depths down to 1335.85 m have been transformed to NGRIP1 depths by linear interpolation between ECM match points and tephra horizons, and from these, GICC05 ages were obtained using the NGRIP1 GICC05 depth-age relationship. For the section from 1335.90 m to 1955.90 m, the NEEM depths have been transformed to NGRIP2 depths by linear interpolation between ECM match points and tephras, and from these, GICC05 ages were obtained using the NGRIP2 GICC05 depth-age relationship (see³⁹ for matching details between NGRIP1, NGRIP2 and NEEM). Note that GICC05 comes with an estimate of the Maximum Counting Error (MCE) of the time scale⁹. In a standard deviation context, the maximum counting error should be regarded as 2-σ⁰⁴¹.

The GICC05modelext time scale has been transferred to NEEM depths in the following ways: Between 1955.95 m to 2203.55 m, the NEEM depths have been transformed to NGRIP2 depths by linear interpolation between ECM match points and tephras, and from these, GICC05modelext ages were obtained using the NGRIP2 GICC05modelext depth-age relationship (see³⁹ for matching details). Below 2203.55 m NEEM depths, the ice in the NEEM core contains time reversals due to foldings. A detailed discussion of the matching and GICC05modelext timescale transfer between NEEM and NGRIP in this folded part of the core can be found in³³. The depth-age relationship for this section is interpolated to 5 cm resolution from the lower resolution relationship given in the supplementary data file in³³.

The AICC2012 time scale has been transferred to NEEM depths in the following ways: For ages younger than 60.2ka (1955.95 m NEEM depth) AICC12 is constructed to be practically identical to GICC05³⁸, so no separate age transfer is needed. Between 1955.95 m to 2196.95 m, the NEEM depths have been transformed to NGRIP2 by linear interpolation between ECM match points and tephras, and from these, AICC2012 ages were obtained using the NGRIP2 AICC2012 depth-age relationship (see³⁹ for matching details, and³⁸ for the NGRIP2 AICC2012 timescale). For the section below 2196.95 m the ice in the NEEM core contains time reversals due to foldings. The depth-age relationship for this section is derived from a matching of the NEEM and EDML cores³³. The lower resolution NEEM-EDML1 timescale given in the³³ supplementary data has been interpolated to 5 cm resolution, and EDML1 ages have subsequently been converted to AICC2012 ages using the relationship between EDML1
<figure><img src="image_5.png" /><figcaption>Fig. 3 The compiled NEEM record from the depth of 1210.5 (corresponding to 8000 y b2k) and below.</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>Fig. 4 δ¹⁸O record and the coverage obtained for every instrument type.</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>Fig. 5 Quality-control metrics based on the measurement of the “check” standard. The top row presents the precision level of every measurement run while at the bottom row we present the offset of the “check” standard measurement from its assigned value- a measure of accuracy for every measurement run.</figcaption></figure>
<figure><img src="image_8.png" /><figcaption>Fig. 6 Histograms of the quality-control metrics of Fig. 5.</figcaption></figure>

<table><caption>**Table 6.** Statistics of the “check” standard quality-control metrics as presented in Fig. 6 for the three instrument types.</caption><thead><tr><td>Instrument</td><td>δ<sup>18</sup>O precision [%]</td><td>δ<sup>18</sup>O accuracy [%]</td><td>δD precision [%]</td><td>δD accuracy [%]</td></tr></thead><tbody><tr><td>L2120-i</td><td>0.065 ± 0.028</td><td>0.023 ± 0.059</td><td>0.524 ± 0.163</td><td>0.151 ± 0.232</td></tr><tr><td>L2130-i</td><td>0.031 ± 0.017</td><td>-0.001 ± 0.035</td><td>0.301 ± 0.098</td><td>0.034 ± 0.146</td></tr><tr><td>L2140-i</td><td>0.034 ± 0.015</td><td>-0.011 ± 0.033</td><td>0.390 ± 0.118</td><td>0.012 ± 0.172</td></tr></tbody></table>

and AICC2012 (see$^{38,42}$ and$^{40}$). It should be noted, that transferring AICC2012 to NEEM from NGRIP2 and from EDML does not yield the exact same ages for all NEEM depths due to matching uncertainties. Between 2196.95 m and 2197.00 m the two transfers of AICC2012 do, however, yield intersecting ages. Hence, to ensure continuity of the NEEM AICC2012 time scale, the AICC2012 ages are transferred from NGRIP2 down to 2196.95 m NEEM depth, while AICC2012 ages are transferred from EDML starting from 2197.00 m NEEM depth.

## Data Records

The dataset$^{32}$ can be found in the PANGAEA data repository (https://doi.org/10.1594/PANGAEA.925552). In Table 5 we list the variables included in the data record together with a short description of each one of them. The full record is plotted in Fig. 3 (δD and δ$^{18}$O) whereas in Fig. 4 we present the contribution of every CRDS model to the measurement of the total record.

## Technical Validation

**Measurement precision and accuracy.** The measurement precision of every run, is estimated from the calculation of the standard deviation of the valid injections of the “check” standard blocks (Vial #2). Measurement precision is given for δ$^{18}$O and δD in columns 8 and 9 of the data set. Similarly, the difference between the assigned and the post-calibration value of the “check” standard is used as a measure of the accuracy of the run. This metric is given in columns 10 and 11 for δ$^{18}$O and δD respectively.

In Fig. 5 we plot the two metrics as a function of depth categorised by the instrument type and in Fig. 6 we bin the data points in histograms. From these two figures it is apparent that the performance of the L2130-i and L2140-i versions of the CRDS spectrometer is superior to that of the L2120-i version. The difference in performance concerns both the precision and the accuracy of the measurement and it appears to be more profound for the δ$^{18}$O measurement. This behaviour is explained by the more accurate spectroscopic corrections, as well as the more precise control of the optical cavity’s temperature and pressure of the newer L2130-i and L2140-i models. The statistics of the distributions for the two metrics are outlined in Table 6. Based on the data of Table 6 we can commend that the obtained accuracy as expressed by the “check” standard offset, is well within the 2σ (95%) interval as defined by the precision of the “check” standard measurement for all the runs of the record.
## Code availability

The Python code used for the transfer, organising of the data, estimation of the precision and accuracy metrics as well as the plots included in this manuscript can be found in https://github.com/vgkinis/neem_isotope_dataDescriptor_code. In the repository, we also provide auxiliary code with basic routines for post-processing of the PANGAEA data file.

Received: 26 January 2021; Accepted: 7 April 2021;

Published online: 26 May 2021

## References

1. Dansgaard, W. Stable isotopes in precipitation. *Tellus* **16**, 436–468, https://doi.org/10.3402/tellusa.v16i4.8993 (1964).

2. Ciais, P. & Jouzel, J. Deuterium and oxygen-18 in precipitation - isotopic model, including mixed cloud processes. *J. Geophys. Res.-Atmospheres* **99**, 16793–16803, https://doi.org/10.1029/94JD00412 (1994).

3. Jouzel, J. et al. Validity of the temperature reconstruction from water isotopes in ice cores. *J. Geophys. Res.-Oceans* **102**, 26471–26487, https://doi.org/10.1029/97JC01283 (1997).

4. Jouzel, J. A brief history of ice core science over the last 50 yr. *Clim. Past* **9**, 2525–2547, https://doi.org/10.5194/cp-9-2525-2013 (2013).

5. NGRIP members. High-resolution record of Northern Hemisphere climate extending into the last interglacial period. *Nature* **431**, 147–151, https://doi.org/10.1038/nature02805 (2004).

6. EPICA Community Members. Eight glacial cycles from an Antarctic ice core. *Nature* **429**, 623–628, https://doi.org/10.1038/nature02599 (2004).

7. Steffensen, J. P. et al. High-resolution Greenland ice core data show abrupt climate change happens in few years. *Science* **321**, 680–684, https://doi.org/10.1126/science.1157707 (2008).

8. Vinther, B. M. et al. A synchronized dating of three Greenland ice cores throughout the Holocene. *J. Geophys. Res.-Atmospheres* **111**, D13102, https://doi.org/10.1029/2005JD006921 (2006).

9. Rasmussen, S. O. et al. A new Greenland ice core chronology for the last glacial termination. *J. Geophys. Res.* **111**, D06102, https://doi.org/10.1029/2005JD006079 (2006).

10. Masson-Delmotte, V. et al. GRIP Deuterium excess reveals rapid and orbital-scale changes in Greenland moisture origin. *Science* **309**, 118–121, https://doi.org/10.1126/science.1108575 (2005).

11. Gkinis, V., Simonsen, S. B., Buchardt, S. L., White, J. W. C. & Vinther, B. M. Water isotope diffusion rates from the NorthGRIP ice core for the last 16,000 years - Glaciological and paleoclimatic implications. *Earth and Planetary Science Letters* **405**, 132–141, https://doi.org/10.1016/j.epsl.2014.08.022 (2014).

12. Buizert, C. et al. Greenland temperature response to climate forcing during the last deglaciation. *Science* (New York, N.Y.) **345**, 1177–1180, https://doi.org/10.1126/science.1254961 (2014).

13. Jones, T. R. et al. Southern Hemisphere climate variability forced by Northern Hemisphere ice-sheet topography. *Nature* **554**, 351–355, https://doi.org/10.1038/nature24669 (2018).

14. Holme, C., Gkinis, V. & Vinther, B. M. Molecular diffusion of stable water isotopes in polar firn as a proxy for past temperatures. *Geochim. Cosmochim. Acta* **225**, 128–145, https://doi.org/10.1016/j.gca.2018.01.015 (2018).

15. Bigeleisen, J., Perlman, M. L. & Prosser, H. C. Conversion of Hydrogenic Materials to Hydrogen for Isotopic Analysis. *Analytical Chemistry* **24**, 1356–1357, https://doi.org/10.1021/ac60068a025 (1952).

16. Epstein, T. & Mayeda, S. Variations of ¹⁸O content of waters from natural sources. *Geochim. Cosmochim. Acta* **4**, 213–224, https://doi.org/10.1016/0016-7037(53)90051-9 (1953).

17. Vaughn, B. H. et al. An automated system for hydrogen isotope analysis of water. *Chemical Geology* **152**, 309–319, https://doi.org/10.1016/S0009-2541(98)00117-X (1998).

18. Begley, I. S. & Scrimgeour, C. M. High-precision δ²H and δ¹⁸O measurement for water and volatile organic compounds by Continuous-Flow Pyrolysis Isotope Ratio Mass Spectrometry. *Analytical Chemistry* **69**, 1530–1535, https://doi.org/10.1021/ac960935r (1997).

19. Gehre, M., Geilmann, H., Richter, J., Werner, R. A. & Brand, W. A. Continuous flow ²H/²¹O and δ¹⁸O analysis of water samples with dual inlet precision. *Rapid Communications In Mass Spectrometry* **18**, 2650–2660, https://doi.org/10.1002/rcm.1672 (2004).

20. Kerstel, E. R. T., van Trigt, R., Dam, N., Reuss, J. & Meijer, H. A. J. Simultaneous determination of the ²H/¹H, ¹⁷O/¹⁶O and ¹⁸O/¹⁶O isotope abundance ratios in water by means of laser spectrometry. *Analytical Chemistry* **71**, 5297–5303, https://doi.org/10.1021/ac990621e (1999).

21. Van Trigt, R., Meijer, H. A. J., Sveinbjörnsdóttir, A. E., Johnsen, S. J. & Kerstel, E. R. T. Measuring stable isotopes of Hydrogen and Oxygen in ice by means of laser spectrometry: the Bølling transition in the Dye-3 (south Greenland) ice core. *Ann. Glaciol.* **35**, 125–130, https://doi.org/10.3189/172756402781816906 (2002).

22. Crosson, E. R. A cavity ring-down analyzer for measuring atmospheric levels of methane, carbon dioxide, and water vapor. *Applied Physics B-Lasers And Optics* **92**, 403–408, https://doi.org/10.1063/1.1139895 (2008).

23. Gkinis, V., Popp, T. J., Johnsen, S. J. & Blunier, T. A continuous stream flash evaporator for the calibration of an IR cavity ring-down spectrometer for the isotopic analysis of water. *Isotopes In Environmental and Health Studies* **46**, 463–475, https://doi.org/10.1080/10256016.2010.538052 (2010).

24. Aemisegger, F. et al. Measuring variations of δ¹⁸O and δ²H in atmospheric water vapour using two commercial laser-based spectrometers: an instrument characterisation study. *Atmos. Meas. Tech.* **5**, 1491–1511, https://doi.org/10.5194/amt-5-1491-2012 (2012).

25. Steen-Larsen, H. C. et al. Continuous monitoring of summer surface water vapor isotopic composition above the Greenland Ice Sheet. *Atmos. Chem. Phys.* **13**, 4815–4828, https://doi.org/10.5194/acp-13-4815-2013 (2013).

26. Benetti, M. et al. Stable isotopes in the atmospheric marine boundary layer water vapour over the Atlantic Ocean, 2012–2015. *Sci. Data* **160128**, https://doi.org/10.1038/sdata.2016.128 (2017).

27. Gkinis, V. et al. Water isotopic ratios from a continuously melted ice core sample. *Atmos. Meas. Tech.* **4**, 2531–2542, https://doi.org/10.5194/amt-4-2531-2011 (2011).

28. Maselli, O. J., Fritzsche, D., Layman, L., McConnell, J. R. & Meyer, H. Comparison of water isotope-ratio determinations using two cavity ring-down instruments and classical mass spectrometry in continuous ice-core analysis. *Isotopes Environ. Health Stud.* **387**–398, https://doi.org/10.1080/10256016.2013.781598 (2013).

29. Emanuelsson, B. D., Baisden, W. T., Bertler, N. A. N., Keller, E. D. & Gkinis, V. High-resolution continuous-flow analysis setup for water isotopic measurement from ice cores using laser spectroscopy. *Atmos. Meas. Tech.* **8**, 2869–2883, https://doi.org/10.5194/amt-8-2869-2015 (2015).

30. Jones, T. R. et al. Improved methodologies for continuous-flow analysis of stable water isotopes in ice cores. *Atmos. Meas. Tech.* **10**, 617–632, https://doi.org/10.5194/amt-10-617-2017 (2017).

31. Jones, T. R. et al. Water isotope diffusion in the WAIS Divide ice core during the Holocene and last glacial. *J. Geophys. Res. Earth Surf.* **122**, 290–309, https://doi.org/10.1002/2016JF003938 (2017).
32. Gkinis, V. et al. NEEM ice core High Resolution (0.05 m) Water Isotope Ratios (18 O/16 O, 2H/1H) covering 8-129 ky b2k, https://doi.org/10.1594/PANGAEA.925552 (2020).

33. NEEM Community Members. Eemian interglacial reconstructed from a Greenland folded ice core. *Nature* **493**, 489–494, https://doi.org/10.1038/nature11789 (2013).

34. Rasmussen, S. O. et al. A stratigraphic framework for abrupt climatic changes during the Last Glacial period based on three synchronized Greenland ice-core records: refining and extending the INTIMATE event stratigraphy. *Quat. Sci. Rev.* **106**, 14–28, https://doi.org/10.1016/j.quascirev.2014.09.007 (2014).

35. Svensson, A. et al. A 60 000 year Greenland stratigraphic ice core chronology. *Clim. Past* **4**, 47–57, https://doi.org/10.5194/cp-4-47-2008 (2008).

36. Wolff, E. W., Chappellaz, J., Blunier, T., Rasmussen, S. O. & Svensson, A. Millennial-scale variability during the last glacial: The ice core record. *Quat. Sci. Rev.* **29**, 2828–2838, https://doi.org/10.1016/j.quascirev.2009.10.013 (2010).

37. Bazin, L. et al. An optimized multi-proxy, multi-site Antarctic ice and gas orbital chronology (AICC2012): 120-800 ka. *Clim. Past* **9**, 1715–1731, https://doi.org/10.5194/cp-9-1715-2013 (2013).

38. Veres, D. et al. The Antarctic ice core chronology (AICC2012): an optimized multi-parameter and multi-site dating approach for the last 120 thousand years. *Clim. Past* **9**, 1733–1748, https://doi.org/10.5194/cp-9-1733-2013 (2013).

39. Rasmussen, S. O. et al. A first chronology for the North Greenland Eemian Ice Drilling (NEEM) ice core. *Clim. Past* **9**, 2713–2730, https://doi.org/10.5194/cp-9-2713-2013 (2013).

40. Govin, A. et al. Sequence of events from the onset to the demise of the Last Interglacial: Evaluating strengths and limitations of chronologies used in climatic archives. *Quat. Sci. Rev.* **129**, 1–36, https://doi.org/10.1016/j.quascirev.2015.09.018 (2015).

41. Andersen, K. K. et al. The Greenland Ice Core Chronology 2005, 15-42 ka. Part 1: constructing the time scale. *Quat. Sci. Rev.* **25**, 3246–3257, https://doi.org/10.1016/j.quascirev.2006.08.002 (2006).

42. Ruth, U. et al. "EDML1": a chronology for the EPICA deep ice core from Dronning Maud Land, Antarctica, over the last 150 000 years. *Clim. Past* **3**, 475–484, https://doi.org/10.5194/cp-3-475-2007 (2007).

## Acknowledgements

We are grateful to the drilling engineers, operators and field assistants that made the retrieval, processing and distribution of the NEEM ice core possible. We would like to acknowledge financial support from the European Research Council under European Union's Seventh Framework Programme for research and innovation (FP7-ENVIRONMENT grant agreement ID: 243908; FP7-IDEAS-ERC grant agreement ID: 246815, 610055; Marie Sklodowska-Curie grant agreement ID: 600207), the Villum Foundation (Grant 00022995, 000280161, 16572), the US National Science Foundation (Award 0806387), the French National Research Agency (grants ANR-07-VULN-0009, ANR-19-MPGA-0001). NEEM is directed and organised by the Centre of Ice and Climate at the Niels Bohr Institute, funded by the Danish National Research Foundation and US NSF, Office of Polar Programs. It is supported by funding agencies and institutions in Belgium (FNRS-CFB and FWO), Canada (NRCan/GSC), China (CAS), Denmark (FIST), France (IPEV, CNRS/INSU, CEA and ANR), Germany (AWI), Iceland (RannIs), Japan (NIPR), South Korea (KOPRI), the Netherlands (NWO/ALW), Sweden (VR), Switzerland (SNF), the United Kingdom (NERC) and the US (US NSF, Office of Polar Programs) and the EU Seventh Framework Programmes Past4Future and WaterundertheIce.

## Author contributions

D.D.J., J.P.S., B.M.V., J.J., V.M.D., A.L., H.O., A.E.S. and J.W.W. directed and planned the water isotope strategy for the NEEM ice core project. V.G., T.J.P. and B.M.V. designed, tested and implemented the sample preparation, laser spectroscopy analysis and data processing protocols utilised in this study with valuable input from A.L., H.C.S.L., B.V., T.J. and M.H. V.G., T.J.P., T.Q., A.K.F., C.T.H., C.M.J., M.L., A.M.L., V.M., N.O.Ø., A.S.P., N.V. and Y.W. performed sample preparation, CRDS and data analysis. B.M.V., E.C., S.O.R. and A.S. contributed to the adaptation of the time scales to the NEEM high resolution isotope record.

## Competing interests

The authors declare no competing interests.

## Additional information

Correspondence and requests for materials should be addressed to V.G.

Reprints and permissions information is available at www.nature.com/reprints.

Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_9.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

The Creative Commons Public Domain Dedication waiver http://creativecommons.org/publicdomain/zero/1.0/ applies to the metadata files associated with this article.

© The Author(s) 2021