Specific conductivity and hydrogen ions measured on two ice core (NEEM and NGRIP)

A stratigraphy-based chronology for the North Greenland Eemian Ice Drilling (NEEM) ice core has been derived by transferring the annual layer counted Greenland Ice Core Chronology 2005 (GICC05) and its model extension (GICC05modelext) from the NGRIP core to the NEEM core using 787 match points of mainly volcanic origin identified in the electrical conductivity measurement (ECM) and dielectrical profiling (DEP) records. Tephra horizons found in both the NEEM and NGRIP ice cores are used to test the matching based on ECM and DEP and provide five additional horizons used for the timescale transfer. A thinning function reflecting the accumulated strain along the core has been determined using a Dansgaard-Johnsen flow model and an isotope-dependent accumulation rate parameterization. Flow parameters are determined from Monte Carlo analysis constrained by the observed depth-age horizons. In order to construct a chronology for the gas phase, the ice age-gas age difference (Delta age) has been reconstructed using a coupled firn densification-heat diffusion model. Temperature and accumulation inputs to the Delta age model, initially derived from the water isotope proxies, have been adjusted to optimize the fit to timing constraints from d15N of nitrogen and high-resolution methane data during the abrupt onset of Greenland interstadials. The ice and gas chronologies and the corresponding thinning function represent the first chronology for the NEEM core, named GICC05modelext-NEEM-1. Based on both the flow and firn modelling results, the accumulation history for the NEEM site has been reconstructed. Together, the timescale and accumulation reconstruction provide the necessary basis for further analysis of the records from NEEM.

GISP2 ice core N2O data

Paleoatmospheric records of trace-gas concentrations recovered from ice cores provide important sources of information on many biogeochemical cycles involving carbon, nitrogen, and oxygen. Here, we present a 106,000-year record of atmospheric nitrous oxide (N2O) along with corresponding isotopic records spanning the last 30,000 years, which together suggest minimal changes in the ratio of marine to terrestrial N2O production. During the last glacial termination, both marine and oceanic N2O emissions increased by 40 ± 8%. We speculate that our records do not support those hypotheses that invoke enhanced export production to explain low carbon dioxide values during glacial periods.

Synchronization of the NGRIP, GRIP, and GISP2 ice cores

We here present a synchronization of the NGRIP, GRIP, and GISP2 ice cores based mainly on volcanic events over the period 14.9-32.45 ka b2k (before AD 2000), corresponding to Marine Isotope Stage 2 (MIS 2) and the end of MIS 3. The matching provides a basis for applying the recent NGRIP-based Greenland Ice Core Chronology 2005 (GICC05) time scale to the GRIP and GISP2 ice cores, thereby making it possible to compare the synchronized palaeoclimate profiles of the cores in detail and to identify relative accumulation differences between the cores. Based on the matching, a period of anomalous high accumulation rates in the GISP2 ice core is detected within the period 16.5-18.3 ka b2k. The d18O and [Ca2+] profiles of the three cores are presented on the common GICC05 time scale and generally show excellent agreement across the stadial-interstadial transitions and across the two characteristic dust events in Greenland Stadial 3. However, large differences between the d18O and [Ca2+] profiles of the three cores are seen in the same period as the 7-9% increase in the GISP2 accumulation rate. We conclude that changes of the atmospheric circulation are likely to have occurred in this period, altering the spatial gradients in Greenland and resulting in larger variations between the records.

Mean air temperatures of Greenland

Melt rate and surface temperature on the Greenland ice sheet are parameterized in terms of snow accumulation, mean annual air temperatur and mean July air temperature. Melt rates are calculated using positive degree-days, and firn warming (i.e. the positive deviation of the temperature at 10-15 m depth from the mean annual air temperature) is estimated from the calculated amount of refrozen melt water in the firn. A comparison between observed and calculated melt rates shows that the parameterization provides a reasonable estimate of the present ablation rates in West Greenland between 61°N and 76°N. The average equilibrium line elevation is estimated to be about 1150 m and 1000 m for West and East Greenland respectively, which is several hundred meter lower than previous estimates. However, the total annual ablation from the ice sheet is found to be about 280 km**3 of water per year which agrees well with most other estimates. The melt-rate model predicts significant melting and consequently significant firn warming even at the highest elevations of the South Greenland ice sheet, whereas a large region of central Greenland north of 70° N experiences little or no summer melting. This agrees with the distribution of the dry snow facics as given by BENSON (1962).

(Table 1) Gas tie points used to correlate the ice cores NGRIP and EDML

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 d18O of O2 (d18Oatm) and methane (CH4) 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, CH4 records alone are not sufficient to construct a common gas timescale between the two cores. Millennial-scale variations of d18Oatm are evidenced over MIS 5 both on the Antarctic and Greenland ice cores and are coupled to CH4 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.

Ice rafted debris in the Nordic Seas

A continuous 3.5 Myr IRD record was produced from Ocean Drilling Program (ODP) Site 907. A timescale based on magnetic polarity chrons, oxygen isotope stratigraphy (for the last 1Myr) and orbital tuning was developed. The record documents a stepwise inception of large-scale glacial cycles in the Nordic Seas region, the first being a marked expansion of the Greenland ice sheet at 3.3 Ma. A second step occurred at 2.74 Ma by an expansion of large scale ice sheets in the Northern Hemisphere. Ice sheet variability around the Nordic Seas was tightly coupled to global ice volume over the past 3.3 Myr. Between 3 and 1 Ma, most of the variance of the IRD signal is in the 41 kyr band, whereas the last 1 Myr is characterized by stronger 100 kyr variance. The Gamma Ray Porosity Evaluator (GRAPE) density record is closely linked with IRD variations and documents sub orbital variability resembling the late Quaternary Heinrich/Bond cycles.

(Table 1) Two-way travel times and depths of picked internal reflection horizons (IRHs) of the NGRIP ice core

During the summer of 2003, a ground-penetrating radar survey around the North Greenland Icecore Project (NorthGRIP) deep ice-core drilling site (75°06' N, 42°20' W; 2957 m a.s.l.) was carried out using a shielded 250 MHz radar system. The drill site is located on an ice divide, roughly 300 km north-northwest of the summit of the Greenland ice sheet. More than 430 km of profiles were measured, covering a 10 km by 10 km area, with a grid centered on the drilling location, and eight profiles extending beyond this grid. Seven internal horizons within the upper 120 m of the ice sheet were continuously tracked, containing the last 400 years of accumulation history. Based on the age-depth and density-depth distribution of the deep core, the internal layers have been dated and the regional and temporal distribution of accumulation rate in the vicinity of NorthGRIP has been derived. The distribution of accumulation shows a relatively smoothly increasing trend from east to west from 145 kg/m**2/a to 200 kg/m**2/a over a distance of 50 km across the ice divide. The general trend is overlain by small-scale variations on the order of 2.5 kg/m**2/a/km, i.e. around 1.5% of the accumulation mean. The temporal variations of the seven periods defined by the seven tracked isochrones are on the order of +-4% of the mean of the last 400 years, i.e. at NorthGRIP ±7 kg/m**2/a. If the regional accumulation pattern has been stable for the last several thousand years during the Holocene, and ice flow has been comparable to today, advective effects along the particle trajectory upstream of NorthGRIP do not have a significant effect on the interpretation of climatically induced changes in accumulation rates derived from the deep ice core over the last 10 kyr.

Proglacial river dataset from the Akuliarusiarsuup Kuua River northern tributary, Southwest Greenland, 2008 - 2013, version 2.0

The dataset contains measurements of river stage and discharge for one sites along the Akuliarusiarsuup Kuua River's northern tributary, with 30 minute temporal resolution between June 2008 and August 2013 This river is a tributary to the Watson River discharging into Kangerlussuaq Fjord by the town of Kangerlussuaq, Southwest Greenland. Additional data of water temperature, air pressure are also provided. Compared to version 1.0 of the dataset, this dataset used a total of 36 in situ discharge observations collected between 2008 and 2012 to construct the rating curve. Furthermore, data of Station AK-004-001 between 2010-09-06T11:30 to 2010-09-07T13:30 have been removed from version 2.0 because these values were likely caused by backflow when a jokulhlaup from a large glacier dammed lake caused increased water levels in the downstreams lake. Thus, data measured at AK-004-001 between 2010-09-06T11:30 to 2010-09-07T13:30 are not representative for the AK-004 catchment.

Northeast Greenland - Digital Bathymetric Compilation, with links to DBM in ArcGIS format

A new digital bathymetric model (DBM) for the Northeast Greenland (NEG) continental shelf (74°N - 81°N) is presented. The DBM has a grid cell size of 250 m × 250 m and incorporates bathymetric data from 30 multibeam cruises, more than 20 single-beam cruises and first reflector depths from industrial seismic lines. The new DBM substantially improves the bathymetry compared to older models. The DBM not only allows a better delineation of previously known seafloor morphology but, in addition, reveals the presence of previously unmapped morphological features including glacially derived troughs, fjords, grounding-zone wedges, and lateral moraines. These submarine landforms are used to infer the past extent and ice-flow dynamics of the Greenland Ice Sheet during the last full-glacial period of the Quaternary and subsequent ice retreat across the continental shelf. The DBM reveals cross-shelf bathymetric troughs that may enable the inflow of warm Atlantic water masses across the shelf, driving enhanced basal melting of the marine-terminating outlet glaciers draining the ice sheet to the coast in Northeast Greenland. Knolls, sinks, and hummocky seafloor on the middle shelf are also suggested to be related to salt diapirism. North-south-orientated elongate depressions are identified that probably relate to ice-marginal processes in combination with erosion caused by the East Greenland Current. A single guyot-like peak has been discovered and is interpreted to have been produced during a volcanic event approximately 55 Ma ago.