Explore the vast range of titles available.

Radio‐echo sounding (RES) shows large‐scale englacial stratigraphic folds are ubiquitous in Greenland's ice sheet. However, there is no consensus yet on how these folds form. Here, we use the full‐Stokes code Underworld2 to simulate ice movements in three‐dimensional convergent flow, mainly considering ice anisotropy due to a crystallographic preferred orientation, vertical viscosity and density gradients in ice layers, and bedrock topography. Our simulated folds show complex patterns and are classified into: large‐scale folds (>100 m amplitude), small‐scale folds (<<100 m) and basal‐shear folds. The amplitudes of large‐scale folds tend to be at their maximum in the middle of the ice column or just below, in accordance with observations in RES data. We conclude that ice anisotropy amplifies the perturbations in ice layers (mainly due to bedrock topography) into large‐scale folds during flow. Density differences between the warm deep ice and cold ice above may enhance fold amplification. Plain Language Summary Polar ice sheets are composed of compacted former snow layers deposited at the ice surface. If not distorted or deformed, these layers are flat or adapt to the underlying bedrock topography. However, vertical radar scans of Greenland's ice sheet show large‐scale folds of up to hundreds of meters in height. To investigate how these large‐scale folds form, we set up a three‐dimensional numerical ice‐sheet model and simulate fold growth. Our modeling emphasizes the distinctive physical properties of ice required for fold formation, notably its anisotropy (the direction dependency of the flow strength) and power‐law rheology (when ice becomes softer with increasing strain rate). These findings help to better explain ice flow dynamics. Key Points Large‐scale fold formation in polar ice sheets is mainly controlled by ice anisotropy and bedrock topography Buoyancy of deep warm ice can further enhance fold amplification The implementation of ice anisotropy should be included in large‐scale ice flow modeling

Understanding post‐depositional processes altering the layer sequence in ice cores is especially needed to avoid misinterpretation of the oldest and most highly thinned layers. The record of soluble and insoluble impurities represents an important part of the paleoclimate proxies in ice cores but is known to be affected through interaction with the ice matrix, diffusion, and chemical reactions. Laser ablation inductively coupled plasma mass spectrometry (LA‐ICP‐MS) has been recognized for its micron‐scale resolution and micro‐destructiveness in ice core impurity analysis. Employing LA‐ICP‐MS for 2D chemical imaging has already revealed a close relationship between the ice grain boundary network and impurity signals with a significant soluble component, such as Na and Mg. Here we show the latest improvements in chemical imaging with LA‐ICP‐MS, by increasing the spatial resolution to 20 μm and extending the simultaneous analysis to also mostly insoluble impurities, such as Al and Fe. All analytes reveal signals of dispersed spots in a sample of an East Greenland ice core. Based on their average size around 50–60 times larger than an average particle and their heterogeneous elemental ratios these spots are interpreted as particle clusters. To distinguish their origin, a simple colocalization classification reveals elemental ratios consistent with marine and mineral dust aerosol. Based on already existing data from cryo‐Raman spectroscopy, we discuss potential ways to integrate the two methods in a future comparison. Such a combined approach may help constraining post‐depositional changes to the dust‐related insoluble impurity components, such as cluster formation and chemical reactions at grain boundaries. Plain Language Summary Aerosols of marine and terrestrial origin delivered to the polar ice sheets are archived in the ice and can be studied via the analysis of ice cores. The chemical composition and size of mineral dust can deliver important information about past climatic changes and atmospheric transport. However, it has already been shown that this insoluble material is not always passively archived in the ice but can undergo changes in its chemical composition and size, for example, by forming particle aggregates. To investigate these processes, it is preferable to study the chemical composition of insoluble particles and their localization within the ice matrix. Here we show how this can be done by a new chemical imaging method for ice using laser ablation inductively coupled plasma mass spectrometry. In a sample of a Greenland ice core we find clear signals of particle clusters, 50–60 times larger than a single particle. Based on their chemical composition, it is possible to differentiate between marine and terrestrial material. We discuss the results against findings previously obtained for the same sample from a different method, cryo‐Raman spectroscopy. Bringing together several methods may provide important added value for a more comprehensive understanding of these important indicators of past climate. Key Points Laser ablation inductively coupled plasma mass spectrometry imaging investigates the geochemical composition and localization of particles in samples of the East Greenland Ice Core Project ice core, East Greenland The maps reveal clusters of insoluble particles 50–60 times larger than an average particle, at intra‐grain locations and grain boundaries Geochemical signals of the particles are consistent with cryo‐Raman spectroscopy and with known dust sources and sea salt aerosol

The flow of glaciers and polar ice sheets is controlled by the highly anisotropic rheology of ice crystals that have hexagonal symmetry (ice lh). To improve our knowledge of ice sheet dynamics, it is necessary to understand how dynamic recrystallization (DRX) controls ice microstructures and rheology at different boundary conditions that range from pure shear flattening at the top to simple shear near the base of the sheets. We present a series of two-dimensional numerical simulations that couple ice deformation with DRX of various intensities, paying special attention to the effect of boundary conditions. The simulations show how similar orientations of c-axis maxima with respect to the finite deformation direction develop regardless of the amount of DRX and applied boundary conditions. In pure shear this direction is parallel to the maximum compressional stress, while it rotates towards the shear direction in simple shear. This leads to strain hardening and increased activity of non-basal slip systems in pure shear and to strain softening in simple shear. Therefore, it is expected that ice is effectively weaker in the lower parts of the ice sheets than in the upper parts. Strain-rate localization occurs in all simulations, especially in simple shear cases. Recrystallization suppresses localization, which necessitates the activation of hard, non-basal slip systems. This article is part of the themed issue ‘Microdynamics of ice’.

Microstructures from deep ice cores reflect the dynamic conditions of the drill location as well as the thermodynamic history of the drill site and catchment area in great detail. Ice core parameters (crystal lattice-preferred orientation (LPO), grain size, grain shape), mesostructures (visual stratigraphy) as well as borehole deformation were measured in a deep ice core drilled at Kohnen Station, Dronning Maud Land (DML), Antarctica. These observations are used to characterize the local dynamic setting and its rheological as well as microstructural effects at the EDML ice core drilling site (European Project for Ice Coring in Antarctica in DML). The results suggest a division of the core into five distinct sections, interpreted as the effects of changing deformation boundary conditions from triaxial deformation with horizontal extension to bedrock-parallel shear. Region 1 (uppermost approx. 450 m depth) with still small macroscopic strain is dominated by compression of bubbles and strong strain and recrystallization localization. Region 2 (approx. 450-1700 m depth) shows a girdle-type LPO with the girdle plane being perpendicular to grain elongations, which indicates triaxial deformation with dominating horizontal extension. In this region (approx. 1000 m depth), the first subtle traces of shear deformation are observed in the shape-preferred orientation (SPO) by inclination of the grain elongation. Region 3 (approx. 1700-2030 m depth) represents a transitional regime between triaxial deformation and dominance of shear, which becomes apparent in the progression of the girdle to a single maximum LPO and increasing obliqueness of grain elongations. The fully developed single maximum LPO in region 4 (approx. 2030-2385 m depth) is an indicator of shear dominance. Region 5 (below approx. 2385 m depth) is marked by signs of strong shear, such as strong SPO values of grain elongation and strong kink folding of visual layers. The details of structural observations are compared with results from a numerical ice sheet model (PISM, isotropic) for comparison of strain rate trends predicted from the large-scale geometry of the ice sheet and borehole logging data. This comparison confirms the segmentation into these depth regions and in turn provides a wider view of the ice sheet. This article is part of the themed issue ‘Microdynamics of ice’.

Deciphering the localisation of solid and dissolved impurities on the micron-scale in glacial ice remains a challenge, but is critical to understand the integrity of ice core records and internal deformation. Here we report on the state-of-the-art in microstructural impurity research by highlighting recent progress in bringing together cryo-Raman spectroscopy and laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS). We show the potential of both methods and discuss possibilities to improve inter-method approaches aiming for a more holistic understanding of the evolution of impurity localisation throughout the ice column, including post-depositional processes. In this framework, we elaborate on future research priorities such as LA-ICP-MS imaging on firn samples and integrating a large cryo-cell with imaging capabilities.

The dynamic mass loss of ice sheets constitutes one of the biggest uncertainties in projections of ice-sheet evolution. One central, understudied aspect of ice flow is how the bulk orientation of the crystal orientation fabric translates to the mechanical anisotropy of ice. Here we show the spatial distribution of the depth-averaged horizontal anisotropy and corresponding directional flow-enhancement factors covering a large area of the Northeast Greenland Ice Stream onset. Our results are based on airborne and ground-based radar surveys, ice-core observations, and numerical ice-flow modelling. They show a strong spatial variability of the horizontal anisotropy and a rapid crystal reorganisation on the order of hundreds of years coinciding with the ice-stream geometry. Compared to isotropic ice, parts of the ice stream are found to be more than one order of magnitude harder for along-flow extension/compression while the shear margins are potentially softened by a factor of two for horizontal-shear deformation. This study uses radio-echo sounding measurements, ice-core data and models to map the spatial variation in ice-crystal orientation in the northeast Greenland Ice Stream and shows how it potentially affects the ice-flow dynamics in this region.

Only a few localised ice streams drain most of the ice from the Greenland Ice Sheet. Thus, understanding ice stream behaviour and its temporal variability is crucially important to predict future sea-level change. The interior trunk of the 700 km-long North-East Greenland Ice Stream (NEGIS) is remarkable due to the lack of any clear bedrock channel to explain its presence. Here, we present a 3-dimensional analysis of the folding and advection of its stratigraphic horizons, which shows that the localised flow and shear margins in the upper NEGIS were fully developed only ca 2000 years ago. Our results contradict the assumption that the ice stream has been stable throughout the Holocene in its current form and show that upper NEGIS-type development of ice streaming, with distinct shear margins and no bed topography relationship, can be established on time scales of hundreds of years, which is a major challenge for realistic mass-balance and sea-level rise projections. The flow of ice streams leaves traces in the stratigraphy of the ice sheets. Made visible by radar, they reveal the history of the upper North East Greenland Ice Stream. The ice stream is found to have existed in its current form for only about the last 2000 years.

Ever since the first deep ice cores were drilled, it has been a challenge to determine their original, in-situ orientation. In general, the orientation of an ice core is lost as the drill is free to rotate during transport to the surface. For shallow ice cores, it is usually possible to match the adjacent core breaks, which preserves the orientation of the ice column. However, this method fails for deep ice cores, such as the EastGRIP ice core in Northeast Greenland. We provide a method to reconstruct ice core orientation using visual stratigraphy and borehole geometry. As the EastGRIP ice core is drilled through the Northeast Greenland Ice Stream, we use information about the directional structures to perform a full geographical re-orientation. We compared the core orientation with logging data from core break matching and the pattern of the stereographic projections of the crystals’ c-axis orientations. Both comparisons agree very well with the proposed orientation method. The method works well for 441 out of 451 samples from a depth of 1375–2120 m in the EastGRIP ice core. It can also be applied to other ice cores, providing a better foundation for interpreting physical properties and understanding the flow of ice.

The increasing catalogue of high-quality ice-penetrating radar data provides a unique insight in the internal layering architecture of the Greenland ice sheet. The stratigraphy, an indicator of past deformation, highlights irregularities in ice flow and reveals large perturbations without obvious links to bedrock shape. In this work, to establish a new conceptual model for the formation process, we analysed the radar data at the onset of the Petermann Glacier, North Greenland, and created a three-dimensional model of several distinct stratigraphic layers. We demonstrate that the dominant structures are cylindrical folds sub-parallel to the ice flow. By numerical modelling, we show that these folds can be formed by lateral compression of mechanically anisotropic ice, while a general viscosity contrast between layers would not lead to folding for the same boundary conditions. We conclude that the folds primarily form by converging flow as the mechanically anisotropic ice is channelled towards the glacier. A range of mechanisms has been proposed for large-scale folding in polar ice sheets. Here, using new three-dimensional reconstructions of such folds in the onset region of the Greenland Petermann Glacier, the authors show that these formed due to flow convergence and the high mechanical anisotropy of ice.

During the past 20 years, multi-channel radar emerged as a key tool for deciphering an ice sheet's internal architecture. To assign ages to radar reflections and connect them over large areas in the ice sheet, the layer genesis has to be understood on a microphysical scale. Synthetic radar trace modelling based on the dielectric profile of ice cores allows for the assignation of observed physical properties’ variations on the decimetre scale to radar reflectors extending from the coring site to a regional or even whole-ice-sheet scale. In this paper we rely on the available dielectric profiling data of the northern Greenland deep ice cores: NGRIP, NEEM and EGRIP. The three records are well suited for assigning an age model to the stratigraphic radar-mapped layers, and linking up the reflector properties to observations in the cores. Our modelling results show that the internal reflections are mainly due to conductivity changes. Furthermore, we deduce fabric characteristics at the EGRIP drill site from two-way-travel-time differences of along and across-flow polarized radarwave reflections of selected horizons (below 980 m). These indicate in deeper parts of the ice column an across-flow concentrated c-axis fabric.