NEEM Ice Core Research Articles

Millennial-scale variability of Greenland dust provenance during the last glacial maximum as determined by single particle analysis

Greenland ice core records exhibited 100-fold higher dust concentrations during the Last Glacial Maximum (LGM) than during the Holocene, and dust input temporal variability corresponded to different climate states in the LGM. While East Asian deserts, the Sahara, and European loess have been suggested as the potential source areas (PSAs) for Greenland LGM dust, millennial-scale variability in their relative contributions within the LGM remains poorly constrained. Here, we present the morphological, mineralogical, and geochemical characteristics of insoluble microparticles to constrain the provenance of dust in Greenland NEEM ice core samples covering cold Greenland Stadials (GS)-2.1a to GS-3 (~ 14.7 to 27.1 kyr ago) in the LGM. The analysis was conducted on individual particles in microdroplet samples by scanning electron microscopy with energy dispersive X-ray spectroscopy and Raman microspectroscopy. We found that the kaolinite-to-chlorite (K/C) ratios and chemical index of alteration (CIA) values were substantially higher (K/C: 1.4 ± 0.7, CIA: 74.7 ± 2.9) during GS-2.1a to 2.1c than during GS-3 (K/C: 0.5 ± 0.1, CIA: 65.8 ± 2.8). Our records revealed a significant increase in Saharan dust contributions from GS-2.1a to GS-2.1c and that the Gobi Desert and/or European loess were potential source(s) during GS-3. This conclusion is further supported by distinctly different carbon contents in particles corresponding to GS-2.1 and GS-3. These results are consistent with previous estimates of proportional dust source contributions obtained using a mixing model based on Pb and Sr isotopic compositions in NEEM LGM ice and indicate millennial-scale changes in Greenland dust provenance that are probably linked to large-scale atmospheric circulation variabilities during the LGM.

Basal debris of the NEEM ice core, Greenland: a window into sub-ice-sheet geology, basal ice processes and ice-sheet oscillations

AbstractWe present new data from the debris-rich basal ice layers of the NEEM ice core (NW Greenland). Using mineralogical observations, SEM imagery, geochemical data from silicates (meteoric10Be, εNd,87Sr/86Sr) and organic material (C/N, δ13C), we characterize the source material, succession of previous glaciations and deglaciations and the paleoecological conditions during ice-free episodes. Meteoric10Be data and grain features indicate that the ice sheet interacted with paleosols and eroded fresh bedrock, leading to mixing in these debris-rich ice layers. Our analysis also identifies four successive stages in NW Greenland: (1) initial preglacial conditions, (2) glacial advance 1, (3) glacial retreat and interglacial conditions and (4) glacial advance 2 (current ice-sheet development). C/N and δ13C data suggest that deglacial environments favored the development of tundra and taiga ecosystems. These two successive glacial fluctuations observed at NEEM are consistent with those identified from the Camp Century core basal sediments over the last 3 Ma. Further inland, GRIP and GISP2 summit sites have remained glaciated more continuously than the western margin, with less intense ice-substratum interactions than those observed at NEEM.

Sea ice fluctuations in the Baffin Bay and the Labrador Sea during glacial abrupt climate changes

Sea ice decline in the North Atlantic and Nordic Seas has been proposed to contribute to the repeated abrupt atmospheric warmings recorded in Greenland ice cores during the last glacial period, known as Dansgaard-Oeschger (D-O) events. However, the understanding of how sea ice changes were coupled with abrupt climate changes during D-O events has remained incomplete due to a lack of suitable high-resolution sea ice proxy records from northwestern North Atlantic regions. Here, we present a subdecadal-scale bromine enrichment (Brenr) record from the NEEM ice core (Northwest Greenland) and sediment core biomarker records to reconstruct the variability of seasonal sea ice in the Baffin Bay and Labrador Sea over a suite of D-O events between 34 and 42 ka. Our results reveal repeated shifts between stable, multiyear sea ice (MYSI) conditions during cold stadials and unstable, seasonal sea ice conditions during warmer interstadials. The shift from stadial to interstadial sea ice conditions occurred rapidly and synchronously with the atmospheric warming over Greenland, while the amplitude of high-frequency sea ice fluctuations increased through interstadials. Our findings suggest that the rapid replacement of widespread MYSI with seasonal sea ice amplified the abrupt climate warming over the course of D-O events and highlight the role of feedbacks associated with late-interstadial seasonal sea ice expansion in driving the North Atlantic ocean-climate system back to stadial conditions.

Volcanism and the Greenland ice cores: A new tephrochronological framework for the last glacial-interglacial transition (LGIT) based on cryptotephra deposits in three ice cores

Chemical profiles from Greenland ice cores show that the frequency of volcanism was higher during the last glacial-interglacial transition (LGIT) and early Holocene, (17–9 ka b2k) than in any other period during the last 110 kyr. This increased frequency has partly been linked to climate-driven melting of the Icelandic ice sheet during the last deglaciation, with regional isostatic changes thought to alter mantle viscosity and lead to more eruptions. Our study is the first to construct a comprehensive tephrochronological framework from Greenland ice cores over the LGIT to aid in the reconstruction of volcanic activity over this period. The framework is based on extensive high-resolution sampling of three Greenland ice cores between 17.4 and 11.6 ka b2k and comprises a total of 64 cryptotephra deposits from the NGRIP, GRIP and NEEM ice cores. We show that many of these tephras are preserved within the core without an associated chemical signature in the ice, which implies that reconstructions of volcanism based solely on glacio-chemical indicators might underestimate the number of events. Single glass shards from each deposit were geochemically characterised to trace the volcanic source and many of these deposits could be correlated between cores. We show that the 64 deposits represent tephra deposits from 42 separate volcanic events, and of these, 39 are from Iceland, two from the north Pacific region (Japan and USA) and one has an unknown source. Six deposits can be correlated to terrestrial and/or marine tephra deposits in the Northern Hemisphere and the remaining 36 are unreported in other archives. We did not locate tephra from the compositionally distinctive Laacher See eruption (∼13 ka b2k) in our records. Combining our new discoveries with the previously published tephra framework, raises the number of individual tephra horizons found in Greenland ice over this interval to 50. This significantly improves the regional tephrochronological framework, our knowledge of the eruptive history of Iceland during the LGIT and provides new tephra constraints over key LGIT climate events. Consequentially, this framework can guide sampling strategies of future tephra studies in the terrestrial and marine realms aiming to link these records to the Greenland ice cores to assess regional climate synchroneity.

High-resolution aerosol concentration data from the Greenland NorthGRIP and NEEM deep ice cores

Abstract. Records of chemical impurities from ice cores enable us to reconstruct the past deposition of aerosols onto polar ice sheets and alpine glaciers. Through this they allow us to gain insight into changes of the source, transport and deposition processes that ultimately determine the deposition flux at the coring location. However, the low concentrations of the aerosol species in the ice and the resulting high risk of contamination pose a formidable analytical challenge, especially if long, continuous and highly resolved records are needed. Continuous flow analysis, CFA, the continuous melting, decontamination and analysis of ice-core samples has mostly overcome this issue and has quickly become the de facto standard to obtain high-resolution aerosol records from ice cores after its inception at the University of Bern in the mid-1990s. Here, we present continuous records of calcium (Ca2+), sodium (Na+), ammonium (NH4+), nitrate (NO3-) and electrolytic conductivity at 1 mm depth resolution from the NGRIP (North Greenland Ice Core Project) and NEEM (North Greenland Eemian Ice Drilling) ice cores produced by the Bern Continuous Flow Analysis group in the years 2000 to 2011 (Erhardt et al., 2021). Both of the records were previously used in a number of studies but were never published in full 1 mm resolution. Alongside the 1 mm datasets we provide decadal averages, a detailed description of the methods, relevant references, an assessment of the quality of the data and its usable resolution. Along the way we will also give some historical context on the development of the Bern CFA system. The data is available in full 1 mm and 10-year-averaged resolution on PANGAEA (https://doi.org/10.1594/PANGAEA.935838, Erhardt et al., 2021)

Discovery of argon in air-hydrate crystals in a deep ice core using scanning electron microscopy and energy-dispersive X-ray spectroscopy

AbstractTiny samples of ancient atmosphere in air bubbles within ice cores contain argon (Ar), which can be used to reconstruct past temperature changes. At a sufficient depth, the air bubbles are compressed by the overburden pressure under low temperature and transform into air-hydrate crystals. While the oxygen (O2) and nitrogen (N2) molecules have indeed been identified in the air-hydrate crystals with Raman spectroscopy, direct observational knowledge of the distribution of Ar at depth within ice sheet and its enclathration has been lacking. In this study, we applied scanning electron microscopy (SEM) and energy-dispersive X-ray spectroscopy (EDS) to five air-hydrate crystals in the Greenland NEEM ice core, finding them to contain Ar and N. Given that Ar cannot be detected by Raman spectroscopy, the method commonly used for O2 and N2, the SEM-EDS measurement method may become increasingly useful for measuring inert gases in deep ice cores.

The Influence of Orbital Forcing on 10Be Deposition in Greenland Over the Glacial Period

Understanding the transport and deposition of the cosmogenic isotope10Be is vital for the application of the isotope data to infer past changes of solar activity, to reconstruct past Earth’s magnetic field intensity and climate change. Here, we use data of the cosmogenic isotope10Be from the Greenland ice cores, namely the NEEM and GRIP ice cores, to identify factors controlling its distribution. After removing the effects of the geomagnetic field on the cosmogenic radionuclide production rate, the results expose imprints of the 20–22 ka precession cycle on the Greenland10Be records of the last glacial period. This finding can further improve the understanding of10Be variability in ice sheets and has the prospect of providing better reconstructions of geomagnetic and solar activity based on cosmogenic radionuclide records.

Atmospheric iron supply and marine productivity in the glacial North Pacific Ocean

Abstract. Iron (Fe) is a key element in the Earth climate system, as it can enhance marine primary productivity in the high-nutrient low-chlorophyll (HNLC) regions where, despite a high concentration of major nutrients, chlorophyll production is low due to iron limitation. Eolian mineral dust represents one of the main Fe sources to the oceans; thus, quantifying its variability over the last glacial cycle is crucial to evaluate its role in strengthening the biological carbon pump. Polar ice cores, which preserve detailed climate records in their stratigraphy, provide a sensitive and continuous archive for reconstructing past eolian Fe fluxes. Here, we show the Northern Hemisphere Fe record retrieved from the NEEM ice core (Greenland), which offers a unique opportunity to reconstruct the past Fe fluxes in a portion of the Arctic over the last 108 kyr. Holocene Fe fluxes (0.042–11.7 ka, 0.5 mg m−2 yr−1) at the NEEM site were 4 times lower than the average recorded over the last glacial period (11.7–108 ka, 2.0 mg m−2 yr−1), whereas they were greater during the Last Glacial Maximum (LGM; 14.5–26.5 ka, 3.6 mg m−2 yr−1) and Marine Isotope Stage 4 (MIS 4; 60–71 ka, 5.8 mg m−2 yr−1). Comparing the NEEM Fe record with paleoceanographic records retrieved from the HNLC North Pacific, we found that the coldest periods, characterized by the highest Fe fluxes, were distinguished by low marine primary productivity in the subarctic Pacific Ocean, likely due to the greater sea ice extent and the absence of major nutrients upwelling. This supports the hypothesis that Fe fertilization during colder and dustier periods (i.e., LGM and MIS 4) was more effective in other regions, such as the midlatitude North Pacific, where a closer relationship between marine productivity and the NEEM Fe fluxes was observed.

New technique for high-precision, simultaneous measurements of CH<sub>4</sub>, N<sub>2</sub>O and CO<sub>2</sub> concentrations; isotopic and elemental ratios of N<sub>2</sub>, O<sub>2</sub> and Ar; and total air content in ice cores by wet extraction

Abstract. Air in polar ice cores provides unique information on past climatic and atmospheric changes. We developed a new method combining wet extraction, gas chromatography and mass spectrometry for high-precision, simultaneous measurements of eight air components (CH4, N2O and CO2 concentrations; δ15N, δ18O, δO2∕N2 and δAr∕N2; and total air content) from an ice-core sample of ∼ 60 g. The ice sample is evacuated for ∼ 2 h and melted under vacuum, and the released air is continuously transferred into a sample tube at 10 K within 10 min. The air is homogenized in the sample tube overnight at room temperature and split into two aliquots for mass spectrometric and gas chromatographic measurements. Care is taken to minimize (1) contamination of greenhouse gases by using a long evacuation time, (2) consumption of oxygen during sample storage by a passivation treatment on sample tubes, and (3) fractionation of isotopic ratios with a long homogenization time for splitting. Precision is assessed by analyzing standard gases with artificial ice and duplicate measurements of the Dome Fuji and NEEM ice cores. The overall reproducibility (1 SD) of duplicate ice-core analyses are 3.2 ppb, 2.2 ppb and 2.9 ppm for CH4, N2O and CO2 concentrations; 0.006 ‰, 0.011 ‰, 0.09 ‰ and 0.12 ‰ for δ15N, δ18O, δO2∕N2 and δAr∕N2; and 0.63 mLSTP kg−1 for total air content, respectively. Our new method successfully combines the high-precision, small-sample and multiple-species measurements, with a wide range of applications for ice-core paleoenvironmental studies.

A first chronology for the East Greenland Ice-core Project (EGRIP) over the Holocene and last glacial termination

Abstract. This paper provides the first chronology for the deep ice core from the East Greenland Ice-core Project (EGRIP) over the Holocene and the late last glacial period. We rely mainly on volcanic events and common peak patterns recorded by dielectric profiling (DEP) and electrical conductivity measurement (ECM) for the synchronization between the EGRIP, North Greenland Eemian Ice Drilling (NEEM) and North Greenland Ice Core Project (NGRIP) ice cores in Greenland. We transfer the annual-layer-counted Greenland Ice Core Chronology 2005 (GICC05) from the NGRIP core to the EGRIP ice core by means of 381 match points, typically spaced less than 50 years apart. The NEEM ice core has previously been dated in a similar way and is only included to support the match-point identification. We name our EGRIP timescale GICC05-EGRIP-1. Over the uppermost 1383.84 m, we establish a depth–age relationship dating back to 14 967 years b2k (years before the year 2000 CE). Tephra horizons provide an independent validation of our match points. In addition, we compare the ratio of the annual layer thickness between ice cores in between the match points to assess our results in view of the different ice-flow patterns and accumulation regimes of the different periods and geographical regions. For the next years, this initial timescale will be the basis for climatic reconstructions from EGRIP high-resolution proxy data sets, e.g. stable water isotopes, chemical impurity or dust records.

Using a composite flow law to model deformation in the NEEM deep ice core, Greenland – Part 1: The role of grain size and grain size distribution on deformation of the upper 2207 m

Abstract. The effect of grain size on strain rate of ice in the upper 2207 m in the North Greenland Eemian Ice Drilling (NEEM) deep ice core was investigated using a rheological model based on the composite flow law of Goldsby and Kohlstedt (1997, 2001). The grain size was described by both a mean grain size and a grain size distribution, which allowed the strain rate to be calculated using two different model end-members: (i) the microscale constant stress model where each grain deforms by the same stress and (ii) the microscale constant strain rate model where each grain deforms by the same strain rate. The model results predict that grain-size-sensitive flow produces almost all of the deformation in the upper 2207 m of the NEEM ice core, while dislocation creep hardly contributes to deformation. The difference in calculated strain rate between the two model end-members is relatively small. The predicted strain rate in the fine-grained Glacial ice (that is, ice deposited during the last Glacial maximum at depths of 1419 to 2207 m) varies strongly within this depth range and, furthermore, is about 4–5 times higher than in the coarser-grained Holocene ice (0–1419 m). Two peaks in strain rate are predicted at about 1980 and 2100 m depth. The prediction that grain-size-sensitive creep is the fastest process is inconsistent with the microstructures in the Holocene age ice, indicating that the rate of dislocation creep is underestimated in the model. The occurrence of recrystallization processes in the polar ice that did not occur in the experiments may account for this discrepancy. The prediction of the composite flow law model is consistent with microstructures in the Glacial ice, suggesting that fine-grained layers in the Glacial ice may act as internal preferential sliding zones in the Greenland ice sheet.

Using a composite flow law to model deformation in the NEEM deep ice core, Greenland – Part 2: The role of grain size and premelting on ice deformation at high homologous temperature

Abstract. The ice microstructure in the lower part of the North Greenland Eemian Ice Drilling (NEEM) ice core consists of relatively fine-grained ice with a single maximum crystallographic preferred orientation (CPO) alternated by much coarser-grained ice with a partial (great circle) girdle or multi-maxima CPO. In this study, the grain-size-sensitive (GSS) composite flow law of Goldsby and Kohlstedt (2001) was used to study the effects of grain size and premelting (liquid-like layer along the grain boundaries) on strain rate in the lower part of the NEEM ice core. The results show that the strain rates predicted in the fine-grained layers are about an order of magnitude higher than in the much coarser-grained layers. The dominant deformation mechanisms, based on the flow relation of Goldsby and Kohlstedt (2001), between the layers is also different, with basal slip rate limited by grain boundary sliding (GBS-limited creep) being the dominant deformation mechanism in the finer-grained layers, while GBS-limited creep and dislocation creep (basal slip rate limited by non-basal slip) contribute both roughly equally to bulk strain in the coarse-grained layers. Due to the large difference in microstructure between finer-grained ice and the coarse-grained ice at premelting temperatures (T>262 K), it is expected that the fine-grained layers deform at high strain rates, while the coarse-grained layers are relatively stagnant. The difference in microstructure, and consequently in viscosity, between impurity-rich and low-impurity ice can have important consequences for ice dynamics close to the bedrock.

Iron in the NEEM ice core relative to Asian loess records over the last glacial-interglacial cycle.

Mineral dust can indirectly affect the climate by supplying bioavailable iron (Fe) to the ocean. Here, we present the records of dissolved Fe (DFe) and total Fe (TDFe) in North Greenland Eemian Ice Drilling (NEEM) ice core over the past 110 kyr BP. The Fe records are significantly negatively correlated with the carbon-dioxide (CO2) concentrations during cold periods. The results suggest that the changes in Fe fluxes over the past 110 kyr BP in the NEEM ice core are consistent with those in Chinese loess records because the mineral-dust distribution is controlled by the East Asian deserts. Furthermore, the variations in the dust input on a global scale are most likely driven by changes in solar radiation during the last glacial–interglacial cycle in response to Earth's orbital cycles. In the last glacial–interglacial cycle, the DFe/TDFe ratios were higher during the warm periods (following the post-Industrial Revolution and during the Holocene and last interglacial period) than during the main cold period (i.e. the last glacial maximum (LGM)), indicating that the aeolian input of iron and the iron fertilization effect on the oceans have a non-linear relationship during different periods. Although the burning of biomass aerosols has released large amounts of DFe since the Industrial Revolution, no significant responses are observed in the DFe and TDFe variations during this period, indicating that severe anthropogenic contamination has no significant effect on the DFe (TDFe) release in the NEEM ice core.

Solar and climate signals revealed by seasonal 10Be data from the NEEM ice core project for the neutron monitor period

10Be in ice cores has been instrumental for reconstructing past changes in solar activity prior to direct observations. For a robust use of these records, it is pivotal to understand the 10Be transport and deposition. However, there are only few high-resolution seasonal 10Be data longer than one full solar cycle (11 years) that could enable a quantification of the influences of atmospheric circulation and deposition processes on the 10Be signal in ice. Here we present a seasonally resolved 10Be data set covering the neutron monitor period (1951–2002) from a firn core connected to the NEEM (North Greenland Eemian Ice Drilling) project. The results suggest that both summer and winter 10Be reflect the production signal induced by solar modulation of galactic cosmic rays. However, superimposed on this solar signal we find additional meteorologically driven influences on 10Be transport and deposition. We found that the tropopause pressure over 30°N represents an important factor influencing NEEM 10Be concentrations on seasonal and annual scales. 10Be deposited in summer also correlates significantly with the tropopause pressure over Greenland suggesting a direct contribution of stratospheric intrusions during summer to the 10Be deposition in Greenland. To correct for these transport/deposition influences, we apply a first-order correction to the 10Be data using a multi-linear regression model. The “climate-corrected” 10Be data shows a comparable skill for reconstructing production rate changes as the 10Be composite record from five different ice cores in Greenland. The results suggest that the correction approach can be a complementary method to the stacking to better isolate the production rate signal from the 10Be data when only limited data are available.

Excess methane in Greenland ice cores associated with high dust concentrations

Ice core records of atmospheric methane (CH4) and its isotopic composition provide important information about biogeochemical cycles in the past. Interpreting these data requires that they faithfully record the composition of the atmosphere. In this study, we describe anomalies of up to 30–40 ppb CH4 that are only observed in dust-rich (>∼60 ng Ca/g ice), glacial-period ice measured with standard melt-refreeze methods. The stable isotopic composition of CH4 is also significantly affected. Results from the GISP2 and NEEM ice cores from Greenland show that excess CH4 is either released or produced in the presence of liquid water in amounts which are highly correlated with the abundance of Ca2+ and mineral dust in the sample. Additional experiments show that excess CH4 is unaffected by the addition of HgCl2 (a microbial inhibitor) and is not related to ice core storage time. Dust concentrations in Antarctic ice cores are an order of magnitude lower than in Greenlandic ice cores and no excess CH4 was observed in samples from the Antarctic WAIS Divide (WD) and South Pole (SPICE) ice cores. While the overall structure of the ice core atmospheric methane history is minimally impacted by excess CH4, the impacts on the isotopic record and on inverse models used to reconstruct CH4 sources are greater. We propose three potential mechanisms to explain the presence of excess CH4: (1) that CH4 is adsorbed on dust particles prior to deposition on the ice sheet and is slowly desorbed during the melt-extraction step of sample analysis; (2) that dust acts as a micro-environment within the ice sheet for methanogenic extremophiles; or (3) that excess CH4 is a product of abiotic degradation of organic compounds during the melt-extraction step of sample analysis.

A Polarimetric Coherence Method to Determine Ice Crystal Orientation Fabric From Radar Sounding: Application to the NEEM Ice Core Region

Ice crystal orientation fabric (COF) records information about past ice-sheet deformation and influences the present-day flow of ice. Polarimetric radar sounding provides a means to infer anisotropic COF patterns due to the associated birefringence of polar ice. Here, we develop a polarimetric coherence (phase-based) method to determine horizontal properties of the COF. The method utilizes the azimuth and depth dependence of the vertical gradient of the hhvv coherence phase to infer the dielectric principal axes and birefringence, which are then related to the second-order fabric orientation tensor. Specifically, under the assumption that one of the orientational eigenvectors is vertical, we can determine the horizontal eigenvectors and the difference between the horizontal eigenvalues (a measure of horizontal fabric asymmetry). The method exploits single-polarized data acquired with varying antenna orientation. It applies to ground-based “multi-polarization” surveys and is demonstrated using data acquired by Center for Remote Sensing of Ice Sheets (CReSIS) using Multi-Channel Coherent Radar Depth Sounder (MCRDS) from the North Greenland Eemian Ice Drilling (NEEM) ice core region in Greenland. The analysis is validated using a combination of polarimetric matrix backscatter simulations and comparison with COF data from the NEEM ice core. The results are consistent with a conventional model of ice deformation at an ice divide where a lateral tension component is present, with minor horizontal COF asymmetry and the greatest horizontal concentration of crystallographic axes orientated near parallel to the ice divide.

High-resolution mass spectrometric characterization of dissolved organic matter from warm and cold periods in the NEEM ice core

Dissolved organic matter (DOM) is an important component of ice cores but is currently poorly characterized. DOM from one Holocene sample (HS, aged at 1600–4500 B.P.) and one Last Glacial Maximum sample (LS, aged at 21000–25000 B.P.) from the North Greenland Eemian Ice Drilling (NEEM) ice core were analyzed by ultra-high resolution Fourier transform ion cyclotron resonance mass spectrometry (FT-ICR-MS). CHO compounds contributed 50% of the compounds identified in negative-ionization mode in these two samples, with significant contributions from organic N, S, and P compounds, likely suggesting that marine DOM was an important source in these samples. Overall, the chemical compositions are similar between these two samples, suggesting their consistent DOM sources. However, subtle differences in the DOM between these two samples are apparent and could indicate differences in source strength or chemistry occurring through both pre- and post-depositional processes. For example, higher relative amounts of condensed carbon compounds in the HS DOM (5%), compared to the LS DOM (2%), suggest potentially important contributions from terrestrial sources. Greater incorporation of P in the observed DOM in the LS DOM (22%), compared to the HS DOM (13%), indicate more active microbiological processes that likely contribute to phosphorus incorporation into the DOM pool. Although these two samples present only a preliminary analysis of DOM in glacial/interglacial periods, the data indicate a need to expand the analysis into a broader range of ice-core samples, geographical locations, and glacial/interglacial periods.

High-resolution isotopic evidence for a potential Saharan provenance of Greenland glacial dust

Dust concentrations in Greenland ice show pronounced glacial/interglacial variations with almost two orders of magnitude increase during the Last Glacial Maximum. Greenland glacial dust was previously sourced to two East Asian deserts: the Taklimakan and Gobi deserts. Here we report the first high-resolution Pb and Sr isotopic evidence for a significant Saharan dust influence in Greenland during the last glacial period, back to ~31 kyr ago, from the Greenland NEEM ice core. We find that during Greenland Stadials 3–5.1 (~31 to 23 kyr ago), the primary dust provenance was East Asia, as previously proposed. Subsequently, the Saharan isotopic signals emerge during Greenland Stadials 2.1a–2.1c (~22.6 to 14.7 kyr ago) and from the late Bølling-Allerød to the Younger Dryas periods (~13.6 to 12 kyr ago), coincident with increased aridity in the Sahara and efficient northward transport of dust during these cold periods. A mixing isotopic model proposes the Sahara as an important source, accounting for contribution to Greenland glacial dust of up to 50%, particularly during Greenland Stadial 2.1b and the late Bølling-Allerød to the Younger Dryas periods. Our findings provide new insights into climate-related dust provenance changes and essential paleoclimatic constraints on dust-climate feedbacks in northern high latitudes.

Simulating the Last Interglacial Greenland stable water isotope peak: The role of Arctic sea ice changes

Last Interglacial (LIG), stable water isotope values (δ18O) measured in Greenland deep ice cores are at least 2.5‰ higher compared to the present day. Previous isotopic climate simulations of the LIG do not capture the observed Greenland δ18O increases. Here, we use the isotope-enabled HadCM3 (UK Met Office coupled atmosphere-ocean general circulation model) to investigate whether a retreat of Northern Hemisphere sea ice was responsible for this model-data disagreement. Our results highlight the potential significance of sea ice changes on the LIG Greenland isotopic maximum. Sea ice loss in combination with increased sea surface temperatures, over the Arctic, affect δ18O: water vapour enriched in heavy isotopes and a shorter distillation path may both increase δ18O values over Greenland. We show, for the first time, that simulations of the response to Arctic sea ice reduction are capable of producing the likely magnitude of LIG δ18O increases at NEEM, NGRIP, GIPS2 and Camp Century ice core sites. However, we may underestimate δ18O changes at the Renland, DYE3 and GRIP ice core locations. Accounting for possible ice sheet changes is likely to be required to produce a better fit to the LIG ice core δ18O values.

Climate information preserved in seasonal water isotope at NEEM: relationships with temperature, circulation and sea ice

Abstract. Analyzing seasonally resolved δ18O ice core data can aid the interpretation of the climate information in ice cores, also providing insights into factors governing the δ18O signal that cannot be deciphered by investigating the annual δ18O data only. However, the seasonal isotope signal has not yet been investigated in northern Greenland, e.g., at the NEEM (North Greenland Eemian Ice Drilling) ice core drill site. Here, we analyze seasonally resolved δ18O data from four shallow NEEM ice cores covering the last 150 years. Based on correlation analysis with observed temperature, we attribute about 70 and 30 % of annual accumulation to summer and winter, respectively. The NEEM summer δ18O signal correlates strongly with summer western Greenland coastal temperature and with the first principal component (PC1) of summer δ18O from multiple seasonally resolved ice cores from central/southern Greenland. However, there are no significant correlations between NEEM winter δ18O data and western Greenland coastal winter temperature or southern/central Greenland winter δ18O PC1. The stronger correlation with temperature during summer and the dominance of summer precipitation skew the annual δ18O signal in NEEM. The strong footprint of temperature in NEEM summer δ18O record also suggests that the summer δ18O record rather than the winter δ18O record is a better temperature proxy at the NEEM site. Despite the dominant signal of the North Atlantic Oscillation (NAO) and the Atlantic Multidecadal Oscillation (AMO) in the central–southern ice core data, both NAO and AMO exert weak influences on NEEM seasonal δ18O variations. The NEEM seasonal δ18O is found to be highly correlated with Baffin Bay sea ice concentration (SIC) in the satellite observation period (1979–2004), suggesting a connection of the sea ice extent with δ18O at NEEM. NEEM winter δ18O significantly correlates with SIC even for the period prior to satellite observation (1901–1978). The NEEM winter δ18O may reflect sea ice variations of Baffin Bay rather than temperature itself. This study shows that seasonally resolved δ18O records, especially for sites with a seasonal precipitation bias such as NEEM, provide a better understanding of how changing air temperature and circulation patterns are associated with the variability in the δ18O records.