Search

The stability of the Greenland ice sheet through past glacial-interglacial cycles provides knowledge that can contribute to understanding the future mass loss from the Greenland ice sheet in a warmer climate. Paleo-climatic records from ice cores provide constraints on the past climate and ice sheet thickness in Greenland through the current interglacial, the Holocene, 11.7 kyr to present, but is limited to a few ice cores from the central areas. In the previous interglacial period, the Eemian, 130 kyr to 110 kyr before present, the ice core constraints are sparse, and beyond the Eemian, the climate evolution is known from Antarctic ice cores and marine sediments. Here we present initial results from an ice flow modelling study using the PISM ice flow model to simulate the evolution of the Greenland ice sheet in the Eemian and the Holocene periods. We discuss how paleo-climatic data from ice cores and marine sediments can be combined with ice flow modelling. We find that the Greenland ice sheet retreated to a minimum volume of up to ∼1.2 m sea-level equivalent smaller than present in the early or mid-Holocene, and that the ice sheet has continued to recover from this minimum up to present day. In all our runs, the ice sheet is approaching a steady state at the end of the 20th century. Our studies show that the Greenland ice sheet evolves on shorter and longer timescales, and that assessment of future mass loss must take into account the history and current state., The 28th IUGG General Assembly (IUGG2023) (Berlin 2023)

Thinning of the deep ice core layers must be considered when the water isotopic composition of the Oldest Ice Core is to be analyzed. From an experimental point of view, a novel instrument combining a micro-destructive cold femtosecond - Laser Ablation (LA) sampling system, that provides high spatial resolution together with minimal usage of ice sample, and a Cavity Ring Down Spectrometer is being built for high-quality water isotope measurements. Laser ablation results in crater formation and its morphology depends on the laser parameters used. Optical images that show crater morphology under different experimental conditions allow crater characterization towards an efficient cold LA sampling. An ablation chamber and a transfer line are both the connecting parts between the LA system and the CRDS instrument. They are to be designed and constructed in the optimal size and shape to collect the ablated mass and guarantee its smooth delivery to the CRDS analyzer with minimum disturbance. Coupling a Laser Ablation system with a CRDS analyzer has already been achieved using a laser operating at the nanosecond regime and a cryo-cell as the ablation chamber. Comparison of the two Laser Ablation systems, by the means of ice sampling and collection of the ablated material, will be of great importance to understand the ablation mechanism and post-ablation processes on ice and further develop a system dedicated to water isotope measurements.

Impurities in polar ice play a critical role in ice flow, deformation, and the integrity of the ice core record. Especially cloudy bands, visible layers with high impurity concentrations, are prominent features in ice from glacial periods. Their physical and chemical properties are poorly understood, highlighting the need to analyse them in more detail. We bridge the gap between decimetre and micrometre scales by combining the visual stratigraphy line scanner, fabric analyser, microstructure mapping, Raman spectroscopy, and laser ablation inductively coupled plasma mass spectrometry 2D impurity imaging. We classified approximately 1300 cloudy bands from glacial ice from the East Greenland Ice-core Project (EGRIP) ice core into seven different types. We determine the localisation and mineralogy of more than 1000 micro-inclusions at 13 depths. The majority of the minerals found are related to terrestrial dust, such as quartz, feldspar, mica, and hematite. We further found carbonaceous particles, dolomite, and gypsum in high abundance. Rutile, anatase, epidote, titanite, and grossular are infrequently observed. The 2D impurity imaging at 20 µm resolution revealed that cloudy bands are clearly distinguishable in the chemical data. Na, Mg, and Sr are mainly present at grain boundaries, whereas dust-related analytes, such as Al, Fe, and Ti, are located in the grain interior, forming clusters of insoluble impurities. We present novel vast micrometre-resolution insights into cloudy bands and describe the differences within and outside these bands. Combining the visual and chemical data results in new insights into the formation of different cloudy band types and could be the starting point for future in-depth studies on impurity signal integrity and internal deformation in deep polar ice cores.

We here describe, document, and make available a wide range of data sets used for annual layer identification in ice cores from DYE-3, GRIP, NGRIP, NEEM, and EGRIP. The data stem from detailed measurements performed both on the main deep cores and shallow cores over more than forty years using many different setups developed by research groups in several countries, and comprise both discrete measurements from cut ice samples and continuous-flow analysis data. The data series were used for the construction of the Greenland Ice-Core Chronology 2005 (GICC05) and/or the revised GICC21. Now that the underlying data are made available, we also release the individual annual layer positions of the GICC05 time scale which are based on these data sets. We hope that the release of the data sets will stimulate further studies of the past climate taking advantage of these highly resolved data series covering a large part of the interior of the Greenland ice sheet.

We present a 1D numerical model which calculates the age of ice around Dome C. It accounts either for melting or for a layer of stagnant ice above the bedrock, depending on the value of an inverted mechanical ice thickness. It is constrained by horizons picked from radar observations and dated using the EPICA Dome C (EDC) ice core age profile. We used 3 different radar datasets with the widest reaching airbourne radar system covering an area of 10,000 km2 and zooming in to 5 km transects over Little Dome C (LDC) with a ground based system. We find that stagnant ice exists in many places including above the LDC relief where the new Beyond EPICA drill site (BELDC) is located. The modelled thickness of this layer of stagnant ice roughly corresponds to the thickness of the basal unit observed in one of the radar surveys and observations made with Autonomous phase-sensitive Radio-Echo Sounder (ApRES). At BELDC, the modelled stagnant ice thickness is 182 ± 63 m and the modelled maximum age (that we define as the age at a maximum age density of 20 kyr m−1) is 1.49 ± 0.18 Ma at a depth of 2505 ± 34 m. This is very similar to all sites situated on the LDC relief such as that of the Million Year Ice Core project being conducted by the Australian Antarctic Division (AAD). The model was also applied to radar data in the area 10–20 km north of EDC (North Patch, NP), where we find either a thin layer of stagnant ice (generally < 60 m) or a very low melt rate (< 0.1 mm yr−1). The modelled maximum age at NP is over 2 Ma in most places, with ice at 1.5 Ma having a resolution of 9–12 kyr m−1 , making it an exciting prospect for a future oldest ice drill site.

An optimization problem is proposed for inferring physical properties of polycrystals given ultrasonic (elastic) wave velocity measurements, made across multiple sample orientations. The feasibility of the method is demonstrated by inferring both the effective grain elastic parameters and the grain c -axis orientation distribution function (ODF) of ice-core samples from Priestley glacier, Antarctica. The method relies on expanding the ODF in terms of a spherical harmonic series, which allows for a non-parametric estimation of the sample ODF. Moreover, any linear combination of the Voigt (strain) and Reuss (stress) homogenization scheme is allowed, although for glacier ice, the exact choice is found to matter little for bulk elastic behaviour, and thus for inferred physical properties, too. Finally, the accuracy of the inferred grain elastic parameters is discussed, including the well-posedness and shortcomings of the inverse problem, relevant for future adoptions in glaciology, geology and elsewhere.

The 13th Symposium on Polar Science/Ordinary sessions [OM] Polar meteorology and glaciology, Wed. 16 Nov.

9 pags., 5 figs., 2 tabs., 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., The research leading to these results has received funding from the European Research Council under the European Union’s Seventh Framework Programme (FP7/2007–2013) grant agreements #243908 “Past4Future. Climate change—Learning from the past climate” and #610055 “Ice2ice”. We thank our colleagues at Centre for Ice and Climate (Copenhagen, Denmark) for their generous contribution and providing the NEEM ice core samples. NEEM is directed and organized by the Center of Ice and Climate at the Niels Bohr Institute 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), Korea (KOPRI), The Netherlands (NWO/ALW), Sweden (VR), Switzerland (SNF), United Kingdom (NERC), and the USA (US NSF, Office of Polar Programs).

Here, we present a compilation of 85 ice temperature profiles from 79 boreholes from the Greenland Ice Sheet and peripheral ice caps, as well as local ice caps in the Canadian Arctic. Only 25 profiles (32 %) were previously available in open-access data repositories. The remaining 54 profiles (68 %) are being made digitally available here for the first time. These newly available profiles, which are associated with pre-2010 boreholes, have been submitted by community members or digitized from published graphics and/or data tables. All 85 profiles are now made available in both absolute (meters) and normalized (0 to 1 ice thickness) depth scales, and are accompanied by extensive metadata. This metadata includes a transparent description of data provenance. The ice temperature profiles span 70 years, with the earliest profile being from 1950 at Camp VI, West Greenland. To highlight the value of this database in evaluating ice flow simulations, we compare the ice temperature profiles from the Greenland Ice Sheet with an ice flow simulation by the Parallel Ice Sheet Model (PISM). We find a cold bias in modeled near-surface ice temperatures within the ablation area, a warm bias in modeled basal ice temperatures at inland cold-bedded sites, and an apparent underestimation of deformational heating in high-strain settings. These biases provide process-level insight on simulated ice temperatures., The Cryosphere Discussions, ISSN:1994-0432, ISSN:1994-0440

The dynamic mass loss of ice sheets due to ice flow constitutes one of the biggest uncertainties in projections of ice-sheet evolution and sea-level rise. One central, understudied aspect of ice flow is how the bulk orientation of the ice-crystal lattice (fabric) translates to the mechanical anisotropy of ice. Here, we present a comprehensive analysis of the depth-averaged fabric distribution and corresponding directional flow enhancement factors covering a large area of the Northeast Greenland Ice Stream onset region. Our results are based on a combination of methods applied to extensive airborne and ground-based radar surveys, ice-core observations, and numerical ice-flow modelling. They show a strong spatial variability of the horizontal fabric anisotropy and a rapid crystal reorganization on the order of hundreds of years coinciding with the ice-stream geometry. The ice in parts of the ice stream is found to be up to one order of magnitude harder for pure-shear deformation in flow direction compared to isotropic ice, while the shear margins are potentially softened for horizontal simple shear by a factor of two. The high spatial resolution and large-scale coverage of our results is unprecedented and significantly contributes towards more accurate ice-flow models and a better understanding of ice-stream dynamics. Moreover, our estimated time scale of crystal reorganization and methodological approaches can potentially also be transferred to the investigation of ice on extraterrestrial bodies and other geologic materials.

Bulk directional enhancement factors are determined for axisymmetric (girdle and single-maximum) orientation fabrics using a transversely isotropic grain rheology with an orientation-dependent non-linear grain fluidity. Compared to grain fluidities that are simplified as orientation independent, we find that bulk strain-rate enhancements for intermediate-to-strong axisymmetric fabrics can be up to a factor of ten larger, assuming stress homogenization over the polycrystal scale. Our work thus extends previous results based on simple basal slip (Schmid) grain rheologies to the transversely isotropic rheology, which has implications for large-scale anisotropic ice-flow modelling that relies on a transversely isotropic grain rheology. In order to derive bulk enhancement factors for arbitrary evolving fabrics, we expand the c-axis distribution in terms of a spherical harmonic series, which allows the rheology-required structure tensors through order eight to easily be calculated and provides an alternative to current structure-tensor-based modelling.

Future sea-level predictions require that the history of the Antarctic Ice Sheet is well understood and constrained by observations. Much of the ice sheets’ ice-dynamic properties are governed by processes at the ice-bed interface which can be imaged with radar sounding surveys. Here we use a combination of ultra-wideband radio-echo sounding data, satellite radar and laser altimetry data, as well as electromagnetic waveform modeling to characterize the properties of the ice base and the evolution of the subglacial morphology of the Jutulstraumen drainage basin (western Dronning Maud Land, Antarctica). Based on the classification of the bed topography, we reconstruct the step-by-step modifications the subglacial landscape has experienced since the beginning of the glaciation of Antarctica, 34 million years ago. Between 2017 and 2020, we find evidence of active episodic cascade-like subglacial water transport along the subglacial valley network. In addition, our high-resolution radio-echo sounding data reveal a cluster of anomalous basal ice units whose material properties we constrain by electromagnetic waveform modeling. Through this, we aim to derive the physical conditions at the ice base, and establish a link to the subglacial hydrology system. The combination of these observations will represent an important step towards a better understanding of large-scale ice-sheet dynamics in western Dronning Maud Land.

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.

An electrically large ultrawideband ultrahigh frequency (UHF) monopole antenna array has been designed to sound up to 3 km of ice and meet the logistical requirements of transportation to arctic regions. The monopole array is comprised of 16 planar subarray modules, which in combination form a 16 m by 17 m Mills cross array configuration to maximize sensitivity and spatial selectivity in both cross-track and along-track directions. Each planar subarray module is 1 m by 2 m in size with a 6.35 cm thick rigid insulation foam panel separating the individual monopole elements from metal foil ground plane on the top such that the maximum radiation is directed to nadir. Each subarray panel consists of 4 by 8 circular monopole antenna elements with a spacing of 0.25 m. Each monopole element is printed on a 130 mm by 80 mm 62 mil FR4 board. The total weight of each subarray panel is 9 kg, making for very lightweight and low-profile antenna construction. The antenna array together with the radar system was deployed to the East Greenland Ice-coring Project site in August 2018 for demonstrating surface-based ice sounding at the UHF band.

The Beyond EPICA project for retrieving the Oldest Ice Core (1.5 Myr) in Antarctica aims at obtaining high resolution climate records and water isotopes will be one of the most important parameters investigated. Given the extremely thin nature of the annual ice core layers, as we get deep down to the core, analysis of such an ice core requires new adopted techniques on water isotopes with high accuracy and precision. Laser ablation (LA) is an established powerful technique used in various fields and it can also be applied in ice sampling serving a dual purpose: a.direct solid-gas transition and b. the smallest amount of sample possible is used for analysis and that makes LA a micro-distructive process. A new instrument which couples LA sampling with the established Cavity Ring Down Spectroscopy (CRDS) for water isotopic analysis is developed. This novel design will allow both fast gas phase sample collection directly from the ice sample and high quality water isotopic measurements. Particular focus was given in the LA system which consists of a High Energy femtosecond IR LASER and the optical elements that focus the LASER beam into the ice surface. The focusing lens system is placed inside a freezer, up above a motorized stage that accomodates the ice sample. An enclosure supplied with dry air flow was build around the optics and tested by the means of humidity experiments. Subsequent series of experiments with varying laser ablation parameters: pulse energy, repetition rate, ablation time, together with the ablated crater characterization allow the evaluation of LA efficiency in ice and thus the optimization of the parameters controlling the ablation mechanism. Understanding the LA mechanism will provide the knowledge to further develop the sampling procedure and efficiently control and guide the vaporized ice into a CRDS instrument for detection.

The diffusion of stable water isotopes in firn and ice is a complex process which can smooth the measured isotope profile and thus remove the high frequency variations. While for the firn this process is relatively well understood and can be accounted for, this is not the case in the deepest parts of the ice core, where additional processes are introduced due to the increased temperatures near bedrock. Combined with the extreme thinning the deep ice has experienced over hundreds of thousands of years, variations up to centennial or even millennial timescales can be heavily attenuated. Understanding how to best recover this signal in the deep ice is crucial to get a reliable record in deep ice cores such as the Beyond EPICA Oldest Ice Core that is currently being drilled.In order to reconstruct the climate signal of this old ice, an accurate estimation of the diffusion length is necessary. Current estimation methods are mostly suitable for firn and shallow ice as they are assuming a rather stationary underlying climate signal. In this contribution, we present a method which approaches the issue without assuming the low frequency climate variability is negligible. Using this method on the high-resolution Dome-C isotope data (doi.pangaea.de/10.1594/PANGAEA.939445), we provide an improved estimate of the diffusion length of the Dome-C ice-core. Both the diffusion length estimate of the deep ice in Dome-C as well as the new method are useful for the interpretation of future deep ice coring projects such as Beyond EPICA Oldest Ice Core.

Large volcanic eruptions occurring in the last glacial period can be detected by their accompanying sulfuric acid deposition in continuous ice cores. Here we employ continuous sulfate and sulfur records from three Greenland and three Antarctic ice cores to estimate the emission strength, the frequency and the climatic forcing of large volcanic eruptions that occurred during the second half of the last glacial period and the early Holocene, 60–9 kyr before 2000 CE (b2k). Over most of the investigated interval the ice cores are synchronized, making it possible to distinguish large eruptions with a global sulfate distribution from eruptions detectable in one hemisphere only. Due to limited data resolution and large variability in the sulfate background signal, particularly in the Greenland glacial climate, we only list Greenland sulfate depositions larger than 20 kg km−2 and Antarctic sulfate depositions larger than 10 kg km−2. With those restrictions, we identify 1113 volcanic eruptions in Greenland and 737 eruptions in Antarctica within the 51 kyr period – for which the sulfate deposition of 85 eruptions is found at both poles (bipolar eruptions). Based on the ratio of Greenland and Antarctic sulfate deposition, we estimate the latitudinal band of the bipolar eruptions and assess their approximate climatic forcing based on established methods. A total of 25 of the identified bipolar eruptions are larger than any volcanic eruption occurring in the last 2500 years, and 69 eruptions are estimated to have larger sulfur emission strengths than the Tambora, Indonesia, eruption (1815 CE). Throughout the investigated period, the frequency of volcanic eruptions is rather constant and comparable to that of recent times. During the deglacial period (16–9 ka b2k), however, there is a notable increase in the frequency of volcanic events recorded in Greenland and an obvious increase in the fraction of very large eruptions. For Antarctica, the deglacial period cannot be distinguished from other periods. This confirms the suggestion that the isostatic unloading of the Northern Hemisphere (NH) ice sheets may be related to the enhanced NH volcanic activity. Our ice-core-based volcanic sulfate records provide the atmospheric sulfate burden and estimates of climate forcing for further research on climate impact and understanding the mechanism of the Earth system.