ALBERT

Description: Cold glaciers at the highest locations of the European Alps have been investigated with great success by drilling ice cores to retrieve their stratigraphic climate records. Findings like the Oetztal ice man have demonstrated that small ice bodies at summit locations of comparatively lower altitudes may also contain old ice if comprising ice frozen to the underlying bedrock. In this case, constraining the maximum age of their basal ice part may help to identify past periods with minimum ice extent in the Alps. Facing ongoing warming and recent years with extremely negative glacier mass balance, these sites may not preserve their unique climate information for much longer, however. Since sampling and dating the lowermost ice is essential, and usually requires substantial logistical (drilling) effort, we utilize here the direct access to basal ice offered by an existing ice cave at Chli Titlis (3030 m), Central Switzerland. Our dedicated approach comprises a combination of standard glaciological tools with the analysis of the isotopic and physical properties and sophisticated radiocarbon dating techniques. By this means we demonstrate that, in comparison to an earlier study at Chli Titlis, stagnant cold basal ice conditions still exist fairly unchanged more than 25 years after the pioneering exploration. Our radiocarbon analysis constrains the maximum age of the ice at Chli Titlis to about 5000 years before present. By this means, the approach presented here will contribute to a future systematic investigation of cold-based summit glaciers also targeting the Eastern Alps.

Description: One of the great challenges in glaciology is the ability to estimate the bulk ice anisotropy in ice sheets and glaciers, which is needed to improve our understanding of ice-sheet dynamics. We investigate the effect of crystal anisotropy on seismic velocities in glacier ice and revisit the framework which is based on fabric eigenvalues to derive approximate seismic velocities by exploiting the assumed symmetry. In contrast to previous studies, we calculate the seismic velocities using the exact c axis angles describing the orientations of the crystal ensemble in an ice-core sample. We apply this approach to fabric data sets from an alpine and a polar ice core. Our results provide a quantitative evaluation of the earlier approximative eigenvalue framework. For near-vertical incidence our results differ by up to 135 m s−1 for P-wave and 200 m s−1 for S-wave velocity compared to the earlier framework (estimated 1 % difference in average P-wave velocity at the bedrock for the short alpine ice core). We quantify the influence of shear-wave splitting at the bedrock as 45 m s−1 for the alpine ice core and 59 m s−1 for the polar ice core. At non-vertical incidence we obtain differences of up to 185 m s−1 for P-wave and 280 m s−1 for S-wave velocities. Additionally, our findings highlight the variation in seismic velocity at non-vertical incidence as a function of the horizontal azimuth of the seismic plane, which can be significant for non-symmetric orientation distributions and results in a strong azimuth-dependent shear-wave splitting of max. 281 m s−1 at some depths. For a given incidence angle and depth we estimated changes in phase velocity of almost 200 m s−1 for P wave and more than 200 m s−1 for S wave and shear-wave splitting under a rotating seismic plane. We assess for the first time the change in seismic anisotropy that can be expected on a short spatial (vertical) scale in a glacier due to strong variability in crystal-orientation fabric (±50 m s−1 per 10 cm). Our investigation of seismic anisotropy based on ice-core data contributes to advancing the interpretation of seismic data, with respect to extracting bulk information about crystal anisotropy, without having to drill an ice core and with special regard to future applications employing ultrasonic sounding.

Description: Cold glaciers at the highest locations of the European Alps have been investigated by drilling ice cores to retrieve their stratigraphic climate records. Findings like the Oetztal ice man have demonstrated that small ice bodies at summit locations of comparatively lower altitudes may also contain old ice if locally frozen to the underlying bedrock. In this case, constraining the maximum age of their lowermost ice part may help to identify past periods with minimum ice extent in the Alps. However, with recent warming and consequent glacier mass loss, these sites may not preserve their unique climate information for much longer. Here we utilized an existing ice cave at Chli Titlis (3030 m), central Switzerland, to perform a case study for investigating the maximum age of cold-based summit glaciers in the Alps. The cave offers direct access to the glacier stratigraphy without the logistical effort required in ice core drilling. In addition, a pioneering exploration had already demonstrated stagnant cold ice conditions at Chli Titlis, albeit more than 25 years ago. Our englacial temperature measurements and the analysis of the isotopic and physical properties of ice blocks sampled at three locations within the ice cave show that cold ice still exists fairly unchanged today. State-of-the-art micro-radiocarbon analysis constrains the maximum age of the ice at Chli Titlis to about 5000 years before present. By this means, the approach presented here will contribute to a future systematic investigation of cold-based summit glaciers, also in the Eastern Alps.

Description: One of the greatest challenges in glaciology, with respect to sea level predictions, is the ability to gain information on bulk ice anisotropy in ice sheets and glaciers, which is urgently needed to improve our understanding of ice-sheet dynamics. Therefore, we investigate the effect of crystal anisotropy on seismic velocities in a glacier. We revisit the framework which is based on fabric eigenvalues to derive approximate seismic velocities by exploiting the assumed symmetry. In contrast to previous studies, we calculate the seismic velocities using the exact c-axis angles describing the orientations of the crystal ensemble in an ice-core sample. We apply this approach to fabric data sets from an Alpine (KCC) and a polar (EDML) ice core. The results allow a quantitative evaluation of the earlier approximative eigenvalue framework. Additionally, our findings highlight the variation in seismic velocity as a function of the horizontal azimuth of the seismic plane, which can be significant in case of non-symmetric orientation distributions and results in a strong azimuth-dependent shear-wave splitting. For the first time, we assess the change in seismic anisotropy that can be expected on a short spatial scale in a glacier due to a strong variability in crystal-orientation fabric. Our investigation of seismic anisotropy based on ice-core data contributes to advancing the interpretation of seismic data, with respect to extracting bulk information about crystal anisotropy without having to drill an ice core and with special regard to future applications employing ultrasonic sounding.

Description: Crystal orientation fabric (COF) analysis provides information about the c-axis orientation of ice grains and the associated anisotropy and microstructural information about deformation and recrystallisation processes within the glacier. This information can be used to introduce modules that fully describe the microstructural anisotropy or at least direction-dependent enhancement factors for glacier modelling. The COF was studied at an ice core that was obtained from the temperate Rhonegletscher, located in the central Swiss Alps. Seven samples, extracted at depths between 2 and 79 m, were analysed with an automatic fabric analyser. The COF analysis revealed conspicuous four-maxima patterns of the c-axis orientations at all depths. Additional data, such as microstructural images, produced during the ice sample preparation process, were considered to interpret these patterns. Furthermore, repeated high-precision global navigation satellite system (GNSS) surveying allowed the local glacier flow direction to be determined. The relative movements of the individual surveying points indicated longitudinal compressive stresses parallel to the glacier flow. Finally, numerical modelling of the ice flow permitted estimation of the local stress distribution. An integrated analysis of all the data sets provided indications and suggestions for the development of the four-maxima patterns. The centroid of the four-maxima patterns of the individual core samples and the coinciding maximum eigenvector approximately align with the compressive stress directions obtained from numerical modelling with an exception for the deepest sample. The clustering of the c axes in four maxima surrounding the predominant compressive stress direction is most likely the result of a fast migration recrystallisation. This interpretation is supported by air bubble analysis of large-area scanning macroscope (LASM) images. Our results indicate that COF studies, which have so far predominantly been performed on cold ice samples from the polar regions, can also provide valuable insights into the stress and strain rate distribution within temperate glaciers.

Description: 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.