Spatial gradients of temperature, accumulation and delta18O-ice in Greenland over a series of Dansgaard-Oeschger events

Air and water stable isotope measurements from four Greenland deep ice cores (GRIP, GISP2, NGRIP and NEEM) are investigated over a series of Dansgaard–Oeschger events (DO 8, 9 and 10), which are representative of glacial millennial scale variability. Combined with firn modeling, air isotope data allow us to quantify abrupt temperature increases for each drill site (1σ = 0.6 °C for NEEM, GRIP and GISP2, 1.5 °C for NGRIP). Our data show that the magnitude of stadial–interstadial temperature increase is up to 2 °C larger in central and North Greenland than in northwest Greenland: i.e., for DO 8, a magnitude of +8.8 °C is inferred, which is significantly smaller than the +11.1 °C inferred at GISP2. The same spatial pattern is seen for accumulation increases. This pattern is coherent with climate simulations in response to reduced sea-ice extent in the Nordic seas. The temporal water isotope (δ18O)–temperature relationship varies between 0.3 and 0.6 (±0.08) ‰ °C−1 and is systematically larger at NEEM, possibly due to limited changes in precipitation seasonality compared to GISP2, GRIP or NGRIP. The gas age−ice age difference of warming events represented in water and air isotopes can only be modeled when assuming a 26% (NGRIP) to 40% (GRIP) lower accumulation than that derived from a Dansgaard–Johnsen ice flow model.

Eemian interglacial reconstructed from a Greenland folded ice core

Efforts to extract a Greenland ice core with a complete record of the Eemian interglacial (130,000 to 115,000 years ago) have until now been unsuccessful. The response of the Greenland ice sheet to the warmer-than-present climate of the Eemian has thus remained unclear. Here we present the new North Greenland Eemian Ice Drilling ('NEEM') ice core and show only a modest ice-sheet response to the strong warming in the early Eemian. We reconstructed the Eemian record from folded ice using globally homogeneous parameters known from dated Greenland and Antarctic ice-core records. On the basis of water stable isotopes, NEEM surface temperatures after the onset of the Eemian (126,000 years ago) peaked at 8±4 degrees Celsius above the mean of the past millennium, followed by a gradual cooling that was probably driven by the decreasing summer insolation. Between 128,000 and 122,000 years ago, the thickness of the northwest Greenland ice sheet decreased by 400±250 metres, reaching surface elevations 122,000 years ago of 130±300 metres lower than the present. Extensive surface melt occurred at the NEEM site during the Eemian, a phenomenon witnessed when melt layers formed again at NEEM during the exceptional heat of July 2012. With additional warming, surface melt might become more common in the future.

Fire in ice: two millennia of boreal forest fire history from the Greenland NEEM ice core

Biomass burning is a major source of greenhouse gases and influences regional to global climate. Pre-industrial fire-history records from black carbon, charcoal and other proxies provide baseline estimates of biomass burning at local to global scales spanning millennia, and are thus useful to examine the role of fire in the carbon cycle and climate system. Here we use the specific biomarker levoglucosan together with black carbon and ammonium concentrations from the North Greenland Eemian (NEEM) ice cores (77.49° N, 51.2° W; 2480 m a.s.l) over the past 2000 years to infer changes in boreal fire activity. Increases in boreal fire activity over the periods 1000–1300 CE and decreases during 700–900 CE coincide with high-latitude NH temperature changes. Levoglucosan concentrations in the NEEM ice cores peak between 1500 and 1700 CE, and most levoglucosan spikes coincide with the most extensive central and northern Asian droughts of the past millennium. Many of these multi-annual droughts are caused by Asian monsoon failures, thus suggesting a connection between low- and high-latitude climate processes. North America is a primary source of biomass burning aerosols due to its relative proximity to the Greenland Ice Cap. During major fire events, however, isotopic analyses of dust, back trajectories and links with levoglucosan peaks and regional drought reconstructions suggest that Siberia is also an important source of pyrogenic aerosols to Greenland.

Modern solar maximum forced late twentieth century Greenland cooling

The abrupt Northern Hemispheric warming at the end of the twentieth century has been attributed to an enhanced greenhouse effect. Yet Greenland and surrounding subpolar North Atlantic remained anomalously cold in 1970s to early 1990s. Here we reconstructed robust Greenland temperature records (North Greenland Ice Core Project and Greenland Ice Sheet Project 2) over the past 2100 years using argon and nitrogen isotopes in air trapped within ice cores and show that this cold anomaly was part of a recursive pattern of antiphase Greenland temperature responses to solar variability with a possible multidecadal lag. We hypothesize that high solar activity during the modern solar maximum (approximately 1950s–1980s) resulted in a cooling over Greenland and surrounding subpolar North Atlantic through the slowdown of Atlantic Meridional Overturning Circulation with atmospheric feedback processes.

Chemical compositions of solid particles present in the Greenland NEEM ice core over the last 110,000 years

This study reports the chemical composition of particles present along Greenland’s North Greenland Eemian Ice Drilling (NEEM) ice core, back to 110,000 years before present. Insoluble and soluble particles larger than 0.45 μm were extracted from the ice core by ice sublimation, and their chemical composition was analyzed using scanning electron microscope and energy dispersive X-ray spectroscopy and micro-Raman spectroscopy. We show that the dominant insoluble components are silicates, whereas NaCl, Na₂SO₄, CaSO ₄, and CaCO₃ represent major soluble salts. For the first time, particles of CaMg(CO₃)₂ and Ca(NO₃)₂ 4H₂O are identified in a Greenland ice core. The chemical speciation of salts varies with past climatic conditions. Whereas the fraction of Na salts (NaCl + Na₂SO₄) exceeds that of Ca salts (CaSO₄+ CaCO₃) during the Holocene (0.6–11.7 kyr B.P.), the two fractions are similar during the Bølling-Allerød period (12.9–14.6 kyr B.P.). During cold climate such as over the Younger Dryas (12.0–12.6 kyr B.P.) and the Last Glacial Maximum (15.0–26.9 kyr B.P.), the fraction of Ca salts exceeds that of Na salts, showing that the most abundant ion generally controls the salt budget in each period. High-resolution analyses reveal changing particle compositions: those in Holocene ice show seasonal changes, and those in LGM ice show a difference between cloudy bands and clear layers, which again can be largely explained by the availability of ionic components in the atmospheric aerosol body of air masses reaching Greenland.

Water isotopic ratios from a continuously melted ice core sample

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

Continuous measurements of methane mixing ratios from ice cores

This work presents a new, field-deployable technique for continuous, high-resolution measurements of methane mixing ratios from ice cores. The technique is based on a continuous flow analysis system, where ice core samples cut along the long axis of an ice core are melted continuously. The past atmospheric air contained in the ice is separated from the melt water stream via a system for continuous gas extraction. The extracted gas is dehumidified and then analyzed by a Wavelength Scanned-Cavity Ring Down Spectrometer for methane mixing ratios. We assess the performance of the new measurement technique in terms of precision (±0.8 ppbv, 1σ), accuracy (±8 ppbv), temporal (ca. 100 s), and spatial resolution (ca. 5 cm). Using a firn air transport model, we compare the resolution of the measurement technique to the resolution of the atmospheric methane signal as preserved in ice cores in Greenland. We conclude that our measurement technique can resolve all climatically relevant variations as preserved in the ice down to an ice depth of at least 1980 m (66 000 yr before present) in the North Greenland Eemian Ice Drilling ice core. Furthermore, we describe the modifications, which are necessary to make a commercially available spectrometer suitable for continuous methane mixing ratio measurements from ice cores.

The response of atmospheric nitrous oxide to climate variations during the last glacial period

Detailed insight into natural variations of the greenhouse gas nitrous oxide (N2O) in response to changes in the Earth's climate system is provided by new measurements along the ice core of the North Greenland Ice Core Project (NGRIP). The presented record reaches from the early Holocene back into the previous interglacial with a mean time resolution of about 75 years. Between 11 and 120 kyr BP, atmospheric N2O concentrations react substantially to the last glacial-interglacial transition (Termination 1) and millennial time scale climate variations of the last glacial period. For long-lasting Dansgaard/Oeschger (DO) events, the N2O increase precedes Greenland temperature change by several hundred years with an increase rate of about 0.8-1.3 ppbv/century, which accelerates to about 3.8-10.7 ppbv/century at the time of the rapid warming in Greenland. Within each bundle of DO events, the new record further reveals particularly low N2O concentrations at the approximate time of Heinrich events. This suggests that the response of marine and/or terrestrial N2O emissions on a global scale are different for stadials with and without Heinrich events.

High-resolution glacial and deglacial record of atmospheric methane by continuous-flow and laser spectrometer analysis along the NEEM ice core

The Greenland NEEM (North Greenland Eemian Ice Drilling) operation in 2010 provided the first opportunity to combine trace-gas measurements by laser spectroscopic instruments and continuous-flow analysis along a freshly drilled ice core in a field-based setting. We present the resulting atmospheric methane (CH4) record covering the time period from 107.7 to 9.5 ka b2k (thousand years before 2000 AD). Companion discrete CH4 measurements are required to transfer the laser spectroscopic data from a relative to an absolute scale. However, even on a relative scale, the high-resolution CH4 data set significantly improves our knowledge of past atmospheric methane concentration changes. New significant sub-millennial-scale features appear during interstadials and stadials, generally associated with similar changes in water isotopic ratios of the ice, a proxy for local temperature. In addition to the midpoint of Dansgaard–Oeschger (D/O) CH4 transitions usually used for cross-dating, sharp definition of the start and end of these events brings precise depth markers (with ±20 cm uncertainty) for further cross-dating with other palaeo- or ice core records, e.g. speleothems. The method also provides an estimate of CH4 rates of change. The onsets of D/O events in the methane signal show a more rapid rate of change than their endings. The rate of CH4 increase associated with the onsets of D/O events progressively declines from 1.7 to 0.6 ppbv yr−1 in the course of marine isotope stage 3. The largest observed rate of increase takes place at the onset of D/O event #21 and reaches 2.5 ppbv yr−1.

NGRIP CH4 concentration from 120 to 10 kyr before present and its relation to a δ15N temperature reconstruction from the same ice core

During the last glacial cycle, Greenland temperature showed many rapid temperature variations, the so-called Dansgaard–Oeschger (DO) events. The past atmospheric methane concentration closely followed these temperature variations, which implies that the warmings recorded in Greenland were probably hemispheric in extent. Here we substantially extend and complete the North Greenland Ice Core Project (NGRIP) methane record from the Preboreal Holocene (PB) back to the end of the last interglacial period with a mean time resolution of 54 yr. We relate the amplitudes of the methane increases associated with DO events to the amplitudes of the local Greenland NGRIP temperature increases derived from stable nitrogen isotope (δ15N) measurements, which have been performed along the same ice core (Kindler et al., 2014). We find the ratio to oscillate between 5 parts per billion (ppb) per °C and 18 ppb °C−1 with the approximate frequency of the precessional cycle. A remarkably high ratio of 25.5 ppb °C−1 is reached during the transition from the Younger Dryas (YD) to the PB. Analysis of the timing of the fast methane and temperature increases reveals significant lags of the methane increases relative to NGRIP temperature for DO events 5, 9, 10, 11, 13, 15, 19, and 20. These events generally have small methane increase rates and we hypothesize that the lag is caused by pronounced northward displacement of the source regions from stadial to interstadial. We further show that the relative interpolar concentration difference (rIPD) of methane is about 4.5% for the stadials between DO events 18 and 20, which is in the same order as in the stadials before and after DO event 2 around the Last Glacial Maximum. The rIPD of methane remains relatively stable throughout the full last glacial, with a tendency for elevated values during interstadial compared to stadial periods.

Temperature reconstruction from 10 to 120 kyr b2k from the NGRIP ice core

In order to reconstruct the temperature of the North Greenland Ice Core Project (NGRIP) site, new measurements of δ15N have been performed covering the time period from the beginning of the Holocene to Dansgaard–Oeschger (DO) event 8. Together with previously measured and mostly published δ15N data, we present for the first time a NGRIP temperature reconstruction for the whole last glacial period from 10 to 120 kyr b2k (thousand years before 2000 AD) including every DO event based on δ15N isotope measurements combined with a firn densification and heat diffusion model. The detected temperature rises at the onset of DO events range from 5 °C (DO 25) up to 16.5 °C (DO 11) with an uncertainty of ±3 °C. To bring measured and modelled data into agreement, we had to reduce the accumulation rate given by the NGRIP ss09sea06bm timescale in some periods by 30 to 35%, especially during the last glacial maximum. A comparison between reconstructed temperature and δ18Oice data confirms that the isotopic composition of the stadial was strongly influenced by seasonality. We evidence an anticorrelation between the variations of the δ18Oice sensitivity to temperature (referred to as α) and obliquity in agreement with a simple Rayleigh distillation model. Finally, we suggest that α might be influenced by the Northern Hemisphere ice sheet volume.

Chlorine isotope composition in chlorofluorocarbons CFC-11, CFC-12 and CFC-113 in firn, stratospheric and tropospheric air

The stratospheric degradation of chlorofluorocarbons (CFCs) releases chlorine, which is a major contributor to the destruction of stratospheric ozone (O3). A recent study reported strong chlorine isotope fractionation during the breakdown of the most abundant CFC (CFC-12, CCl2F2, Laube et al., 2010a), similar to effects seen in nitrous oxide (N2O). Using air archives to obtain a long-term record of chlorine isotope ratios in CFCs could help to identify and quantify their sources and sinks. We analyse the three most abundant CFCs and show that CFC-11 (CCl3F) and CFC-113 (CClF2CCl2F) exhibit significant stratospheric chlorine isotope fractionation, in common with CFC-12. The apparent isotope fractionation (ϵapp) for mid- and high-latitude stratospheric samples are respectively −2.4 (0.5) and −2.3 (0.4) ‰ for CFC-11, −12.2 (1.6) and −6.8 (0.8) ‰ for CFC-12 and −3.5 (1.5) and −3.3 (1.2) ‰ for CFC-113, where the number in parentheses is the numerical value of the standard uncertainty expressed in per mil. Assuming a constant isotope composition of emissions, we calculate the expected trends in the tropospheric isotope signature of these gases based on their stratospheric 37Cl enrichment and stratosphere–troposphere exchange. We compare these projections to the long-term δ (37Cl) trends of all three CFCs, measured on background tropospheric samples from the Cape Grim air archive (Tasmania, 1978–2010) and tropospheric firn air samples from Greenland (North Greenland Eemian Ice Drilling (NEEM) site) and Antarctica (Fletcher Promontory site). From 1970 to the present day, projected trends agree with tropospheric measurements, suggesting that within analytical uncertainties, a constant average emission isotope delta (δ) is a compatible scenario. The measurement uncertainty is too high to determine whether the average emission isotope δ has been affected by changes in CFC manufacturing processes or not. Our study increases the suite of trace gases amenable to direct isotope ratio measurements in small air volumes (approximately 200 mL), using a single-detector gas chromatography–mass spectrometry (GC–MS) system.