scholarly journals Deriving seismic velocities on the micro-scale from c-axis orientations in ice cores

2018 ◽  
Author(s):  
Johanna Kerch ◽  
Anja Diez ◽  
Ilka Weikusat ◽  
Olaf Eisen

Abstract. 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.

2015 ◽  
Vol 9 (1) ◽  
pp. 367-384 ◽  
Author(s):  
A. Diez ◽  
O. Eisen

Abstract. A preferred orientation of the anisotropic ice crystals influences the viscosity of the ice bulk and the dynamic behaviour of glaciers and ice sheets. Knowledge about the distribution of crystal anisotropy is mainly provided by crystal orientation fabric (COF) data from ice cores. However, the developed anisotropic fabric influences not only the flow behaviour of ice but also the propagation of seismic waves. Two effects are important: (i) sudden changes in COF lead to englacial reflections, and (ii) the anisotropic fabric induces an angle dependency on the seismic velocities and, thus, recorded travel times. A framework is presented here to connect COF data from ice cores with the elasticity tensor to determine seismic velocities and reflection coefficients for cone and girdle fabrics. We connect the microscopic anisotropy of the crystals with the macroscopic anisotropy of the ice mass, observable with seismic methods. Elasticity tensors for different fabrics are calculated and used to investigate the influence of the anisotropic ice fabric on seismic velocities and reflection coefficients, englacially as well as for the ice–bed contact. Hence, it is possible to remotely determine the bulk ice anisotropy.


2018 ◽  
Vol 12 (5) ◽  
pp. 1715-1734 ◽  
Author(s):  
Johanna Kerch ◽  
Anja Diez ◽  
Ilka Weikusat ◽  
Olaf Eisen

Abstract. 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.


2014 ◽  
Vol 8 (4) ◽  
pp. 4349-4395 ◽  
Author(s):  
A. Diez ◽  
O. Eisen

Abstract. A preferred orientation of the anisotropic ice crystals influences the viscosity of the ice bulk and the dynamic behaviour of glaciers and ice sheets. Knowledge about the distribution of crystal anisotropy, to understand its contribution to ice dynamics, is mainly provided by crystal orientation fabric (COF) data from ice cores. However, the developed anisotropic fabric does not only influence the flow behaviour of ice, but also the propagation of seismic waves. Two effects are important: (i) sudden changes in COF lead to englacial reflections and (ii) the anisotropic fabric induces an angle dependency on the seismic velocities and, thus, also recorded traveltimes. A framework is presented here to connect COF data with the elasticity tensor to determine seismic velocities and reflection coefficients for cone and girdle fabrics from ice-core data. We connect the microscopic anisotropy of the crystals with the macroscopic anisotropy of the ice mass, observable with seismic methods. Elasticity tensors for different fabrics are calculated and used to investigate the influence of the anisotropic ice fabric on seismic velocities and reflection coefficients, englacially as well as for the ice-bed contact. Our work, therefore, provides a contribution to remotely determine the state of bulk ice anisotropy.


2020 ◽  
Author(s):  
Felix Kästner ◽  
Simona Pierdominici ◽  
Judith Elger ◽  
Christian Berndt ◽  
Alba Zappone ◽  
...  

<p>Deeply rooted thrust zones are key features of tectonic processes and the evolution of mountain belts. Exhumed and deeply-eroded orogens like the Scandinavian Caledonides allow to study such systems from the surface. Previous seismic investigations of the Seve Nappe Complex have shown indications for a strong but discontinuous reflectivity of this thrust zone, which is only poorly understood. The correlation of seismic properties measured on borehole cores with surface seismic data can help to constrain the origin of this reflectivity. In this study, we compare seismic velocities measured on cores to in situ velocities measured in the borehole. The core and downhole velocities deviate by up to 2 km/s. However, velocities of mafic rocks are generally in close agreement. Seismic anisotropy increases from about 5 to 26 % at depth, indicating a transition from gneissic to schistose foliation. Differences in the core and downhole velocities are most likely the result of microcracks due to depressurization of the cores. Thus, seismic velocity can help to identify mafic rocks on different scales whereas the velocity signature of other lithologies is obscured in core-derived velocities. Metamorphic foliation on the other hand has a clear expression in seismic anisotropy. To further constrain the effects of mineral composition, microstructure and deformation on the measured seismic anisotropy, we conducted additional microscopic investigations on selected core samples. These analyses using electron-based microscopy and X-ray powder diffractometry indicate that the anisotropy is strongest for mica schists followed by amphibole-rich units. This also emphasizes that seismic velocity and anisotropy are of complementary importance to better distinguish the present lithological units. Our results will aid in the evaluation of core-derived seismic properties of high-grade metamorphic rocks at the COSC-1 borehole and elsewhere.</p>


2020 ◽  
Author(s):  
Daniele FIlippi ◽  
Chiara Giorio

<p>The Beyond EPICA Oldest Ice (BEOI) project will drill an ice core dating back to 1.5 million-years (1.5 Myr) ago. This ice core is of particular interest to the scientific community as it will be the only one covering the climate history of the Mid Pleistocene Transition, when glacial-interglacial cycles changed from a 40 Kyr to 100 Kyr cyclicity, and for which causes are not well understood currently. Obtaining useful climatic information beyond 800 Kyr represents an analytical challenge due to the fact that the deepest section of the ice core is very compact and the amount of sample available is very low.</p><p>Current analytical methods for the determination of organics in ice are characterized by a large number of steps that requires large amounts of sample for a single analysis. This results in the loss of the high time resolution desired from ice cores which is particularly problematic for deeper (i.e. older) records where the ice is more compact.</p><p>This work aims at combining the growing field of microfluidics with improvements to conventional mass spectrometry to allow for continuous analysis of organics in ice cores, melted in continuous on a melting-head. In fact, microfluidic is a powerful technology in which, only a small amount of liquid (10<sup>-9</sup>-10<sup>-18</sup> liters) is manipulated and controlled with an extremely high precision. The method invokes a three-step process: (1) the melted ice core sample is sent to a nebulizer to produce aerosol, then (2) the aerosol is dried to remove water content and concentrate the sample, and (3) the aerosol is sent to a mass spectrometer for continuous analysis through a modified electrospray ionization (ESI) probe.</p><p>This novel system, once operational, can be applied to a range of ice cores but is especially useful for older ice cores given the stratification of deeper segments. It will allow the research community to measure organic compounds with a high time resolution, even in the oldest of ice, to retrieve paleoclimatic information that would otherwise be lost using traditional methods.</p>


2014 ◽  
Vol 2 (1) ◽  
pp. SB69-SB77 ◽  
Author(s):  
Niven Shumaker ◽  
Daniel Haymond ◽  
Joe Martin

A geopressure interpretation technique known as the seismic velocity method is a common workflow in which shale compaction functions are characterized at offset control wells, matched to interval seismic velocities, and then used to predictively calculate geopressure away from well control. The seismic velocity method is used to interpret the expected geopressure profile at the Deep Blue subsalt exploration well in Green Canyon 723 in the deep water Gulf of Mexico. The Deep Blue prospect is distinct from other prospects in the play fairway in that the prospective section is overlain by a salt withdrawal minibasin, whereas the offsetting fields are positioned either along the flanks of minibasins or under a thick allochthonous salt canopy. Predrill geopressure interpretations using numerous tomographic imaging velocity data sets shows a large degree of consistency with the magnitude of geopressure encountered in offsetting supra salt and subsalt fields. Results from the Deep Blue 1 exploration well indicate the predrill geopressure interpretation from interval seismic velocities failed to anticipate the extreme degree overpressure encountered in the subsalt section of the well due to poor deep velocity resolution and an “unloaded” compaction signature. The magnitude of overpressure in the primary section is attributed to the emplacement of an unconformable halokinetic sequence over the primary subsalt basin. An interpretive paradigm is described in which the Deep Blue pressure cell is created through two halokinetic episodes: (1) rapid progradation of a salt canopy followed by (2) subsequent salt withdrawal and emplacement of an overlying minibasin. The linkage between halokinetic sequences, burial history, and the development of overpressure can be used to predictively characterize subsalt geopressure environments.


2019 ◽  
Author(s):  
Seyedhamidreza Mojtabavi ◽  
Frank Wilhelms ◽  
Eliza Cook ◽  
Siwan Davies ◽  
Giulia Sinnl ◽  
...  

Abstract. This paper provides the first chronology for the deep ice core from the East GReenland Ice-core Project (EGRIP) over the Holocene and late last glacial period. We rely mainly on volcanic events and common patterns of peaks in dielectric profiling (DEP), electrical conductivity measurements (ECM) and tephra records for the synchronization between the EGRIP, NEEM and NGRIP ice cores in Greenland. We transfer the annual-layer-counted Greenland Ice Core Chronology 2005 (GICC05) timescale from the NGRIP core to the EGRIP ice core by means of 373 match points. The NEEM ice core is only used for supporting match-point identification. We name our EGRIP time scale GICC05-EGRIP-1. Over the uppermost 1383.84 m, we establish a depth–age relationship dating back to 14,965 a b2k (years before the year 2000 CE). Tephra horizons provide an independent validation of our match points. In addition, we compare the ratio of annual layer thicknesses 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. This initial timescale is the basis of interpretation and refinement of the presently derived EGRIP high-resolution data sets of chemical impurities.


1988 ◽  
Vol 10 ◽  
pp. 209-209
Author(s):  
C.C. Langway ◽  
K. Goto-Azuma

Measurements of the chemical constituents in polar-snow deposits translate into chronological records representing a history of atmospheric composition for the periods involved. The 2037 m deep continuous and undisturbed ice core recovered at Dye 3, Greenland between 1979 and 1981 contains a temporal record of sequential snow deposits for the past 9 × 104 a B.P. (Dansgaard and others 1985). The upper 90 m of the deep core were unsuitable for chemistry studies, but stratigraphic continuity with present-day accumulation was obtained by hand-excavating a 5.4 m deep pit and augering two shallow cores to 138 and 113 m depths. The pit and shallow cores represent the last two centuries of snow precipitation.To date, over 6000 individual samples of the pit, shallow and deep ice cores have been measured by ion chromatography for Cl−, NO3−, and SO42− in the field and laboratory (Herron and Langway 1985, Finkel and Langway 1985, Finkel and others 1986), under clean-room conditions. All pit and shallow-core samples were prepared in a continuous sequence of eight samples per year, as identified by other stable and radioactive isotope-dating methods. The deep ice-core samples were selected and prepared from core intervals spaced over the 2037 m profile from time units which showed evidence of abrupt or transitory periods in climate change or volcanic disturbances, as defined by stable isotopes (Dansgaard and others 1985), atmospheric gases (Oeschger and others 1985) and dust (Hammer and others 1985).Approximately 1700 new measurements from the Dye 3 samples are included in this study. Variability in the chemical constituents and their concentration levels is present and meaningful on a short-term and long-term basis. The time units measured represent seasons, years, decades, centuries and longer geological periods. Particular attention is given to two new high- and low-frequency detailed chronological data sets from (1) a continuous 26 m core profile, representing 3000 years, extending from the Holocene/Wisconsin boundary back into the late Wisconsin and (2) measurements made on 106 samples spaced every 2 m over the Wisconsin-age ice from 1786 to 2008 m.


2003 ◽  
Vol 37 ◽  
pp. 317-324 ◽  
Author(s):  
Beatriz Benjumea ◽  
Yury Ya. Macheret ◽  
Francisco J. Navarro ◽  
Teresa Teixidό

AbstractRadio-wave velocity measurements in temperate and polythermal glaciers, combined with dielectric mixture formulae by Looyenga or Paren, have been used during the last decade to estimate the water content in temperate ice. We have used a similar mixture formula by Riznichenko, but based on elastic properties of the material, to estimate the water content from seismic velocity data. To compare the suitability of the two methods, we have used seismic and radar data from a temperate glacier on an Antarctic island. The estimated water contents are within 0.4–2.3% (average 1.2 ±0.6%) when radio-wave velocities are used, and within 0.9–3.2% (average 2.2±0.9%) when seismic velocities are used. These results are similar to those directly measured from ice cores and to those estimated from radar data on other temperate glaciers. The water-content estimates from seismic data are higher than those from radar data, which we attribute to the different behaviour of seismic and radar velocities as functions of density. Near-surface conditions (ice–firn conditions, presence of crevasses, etc.) have a strong influence on the propagation of elastic and electromagnetic waves, and thus on the accuracy of the velocity determinations and water-content estimates, and so should not be disregarded.


2020 ◽  
Author(s):  
Ling Fang ◽  
Theo Jenk ◽  
Thomas Singer ◽  
Shugui Hou ◽  
Margit Schwikowski

Abstract. High-alpine glaciers are valuable archives of past climatic and environmental conditions. The interpretation of the preserved signal requires a precise chronology. Radiocarbon (14C) dating of the water-insoluble organic carbon (WIOC) fraction has become an important dating tool to constrain the age of ice cores from mid-latitude and low-latitude glaciers. However, in some cases this method is restricted by the low WIOC concentration in the ice. In this work, we report first 14C dating results using the dissolved organic carbon (DOC) fraction, which is present at concentrations of at least a factor of two higher than the WIOC fraction. We evaluated this new approach by comparison to the established WIO14C dating based on parallel ice core sample sections from four different Eurasian glaciers covering an age range of several hundred to around 20’000 years. 14C dating of the two fractions yielded comparable ages with WIO14C revealing a slight, barely significant, systematic offset towards older ages. Our data suggests this to be caused by incompletely removed carbonate from mineral dust (14C depleted) contributing to the WIOC fraction. While in the DOC extraction procedure inorganic carbon is monitored to ensure complete removal, the average removal efficiency for WIOC samples was here estimated to be ~96%. We did not find any indication of in-situ production systematically contributing to DO14C as suggested in a previous study. By using the DOC instead of the WIOC fraction for 14C dating, the required ice mass can be reduced to typically ~250 g, yielding a precision of ±200 years or even better if sample sizes typically required for WIO14C dating are used. This study shows the potential of pushing radiocarbon dating of ice forward even to remote and Polar Regions, where the carbon content in the ice is particularly low, when applying the DOC fraction for 14C dating.


Sign in / Sign up

Export Citation Format

Share Document