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

We report high resolution measurements of the stable isotope ratios of ancient ice (δ18O, δ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 δ18O 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 1-4 offering a picture of past climate that extends as far as 120,000 y in the past 5 for the case of Greenland and 800,000 y for the case of Antarctica 6 . 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: , 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 [7][8][9][10][11][12][13][14] .
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 [15][16][17][18][19] . With the advent of fast electronics and room-temperature, turn-key, solid state laser sources, there emerged various optical (2021) 8:141 | https://doi.org/10.1038/s41597-021-00916-9 www.nature.com/scientificdata www.nature.com/scientificdata/ 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 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][25][26] while in combination with continuous melter systems, and micro-volume flash vaporisers, they yielded ice core data sets of unprecedented resolution 23,[27][28][29][30][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 www.nature.com/scientificdata www.nature.com/scientificdata/ 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 22 . 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 δ 18  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  www.nature.com/scientificdata www.nature.com/scientificdata/ 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 27 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) where the slope α and the intercept β are for the example of the low water standards triplet equal to: NEEMraw D C02raw NEEMraw 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.    www.nature.com/scientificdata www.nature.com/scientificdata/ 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 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 39 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 9 . In a standard deviation context, the maximum counting error should be regarded as 2-σ0 41 .
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 39 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 33 . 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 33 .
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 38 , 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 39 for matching details, and 38 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 33 . The lower resolution NEEM-EDML1 timescale given in the 33 supplementary data has been interpolated to 5 cm resolution, and EDML1 ages have subsequently been converted to AICC2012 ages using the relationship between EDML1  Table 5. Name and description of the variables included in the data file.

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.  Table 6. Statistics of the "check" standard quality-control metrics as presented in Fig. 6 for the three instrument types.