Annual layering in the NGRIP ice core during the Eemian

The Greenland NGRIP ice core continuously covers the period from present day back to 123 ka before present, which includes several thousand years of ice from the previous interglacial period, MIS 5e or the Eemian. In the glacial part of the core, annual layers can be identified from impurity records and visual stratigraphy, and stratigraphic layer counting has been performed back to 60 ka. In the deepest part of the core, however, the ice is close to the pressure melting point, the visual stratigraphy is dominated by crystal boundaries, and annual layering is not visible to the naked eye. In this study, we apply a newly developed setup for high-resolution ice core impurity analysis to produce continuous records of dust, sodium and ammonium concentrations as well as conductivity of melt water. We analyzed three 2.2 m sections of ice from the Eemian and the glacial inception. In all of the analyzed ice, annual layers can clearly be recognized, most prominently in the dust and conductivity profiles. Part of the samples is, however, contaminated in dust, most likely from drill liquid. It is interesting that the annual layering is preserved despite a very active crystal growth and grain boundary migration in the deep and warm NGRIP ice. Based on annual layer counting of the new records, we determine a mean annual layer thickness close to 11 mm for all three sections, which, to first order, confirms the modeled NGRIP time scale (ss09sea). The counting does, however, suggest a longer duration of the climatically warmest part of the NGRIP record (MIS5e) of up to 1 ka as compared to the model estimate. Our results suggest that stratigraphic layer counting is possible basically throughout the entire NGRIP ice core, provided sufficiently highly-resolved profiles become available.

Design-based stereology to quantify structural properties of artificial and natural snow using thin sections

The quantification of the structural properties of snow is traditionally based on model-based stereology. Model-based stereology requires assumptions about the shape of the investigated structure. Here, we show how the density, specific surface area, and grain boundary area can be measured using a design-based method, where no assumptions about structural properties are necessary. The stereological results were also compared to X-ray tomography to control the accuracy of the method. The specific surface area calculated with the stereological method was 19.8 ± 12.3% smaller than with X-ray tomography. For the density, the stereological method gave results that were 11.7 ± 12.1% larger than X-ray tomography. The statistical analysis of the estimates confirmed that the stereological method and the sampling used are accurate. This stereological method was successfully tested on artificially produced ice beads but also on several snow types. Combining stereology and polarisation microscopy provides a good estimate of grain boundary areas in ice beads and in natural snow, with some limitatio

Characterization of granite matrix porosity and pore-space geometry by in situ and laboratory methods

Most available studies of interconnected matrix porosity of crystalline rocks are based on laboratory investigations; that is, work on samples that have undergone stress relaxation and were affected by drilling and sample preparation. The extrapolation of the results to in situ conditions is therefore associated with considerable uncertainty, and this was the motivation to conduct the ‘in situ Connected Porosity’ experiment at the Grimsel Test Site (Central Swiss Alps). An acrylic resin doped with fluorescent agents was used to impregnate the microporous granitic matrix in situ around an injection borehole, and samples were obtained by overcoring. The 3-D structure of the porespace, represented by microcracks, was studied by U-stage fluorescence microscopy. Petrophysical methods, including the determination of porosity, permeability and P -wave velocity, were also applied. Investigations were conducted both on samples that were impregnated in situ and on non-impregnated samples, so that natural features could be distinguished from artefacts. The investigated deformed granites display complex microcrack populations representing a polyphase deformation at varying conditions. The crack population is dominated by open cleavage cracks in mica and grain boundary cracks. The porosity of non-impregnated samples lies slightly above 1 per cent, which is 2–2.5 times higher than the in situ porosity obtained for impregnated samples. Measurements of seismic velocities (Vp ) on spherical rock samples as a function of confining pressure, spatial direction and water saturation for both non-impregnated and impregnated samples provide further constraints on the distinction between natural and induced crack types. The main conclusions are that (1) an interconnected network of microcracks exists in the whole granitic matrix, irrespective of the distance to ductile and brittle shear zones, and (2) conventional laboratory methods overestimate the matrix porosity. Calculations of contaminant transport through fractured media often rely on matrix diffusion as a retardation mechanism.

Microfabric memory of vein quartz for strain localization in detachment faults: A case study on the Simplon fault zone

This manuscript deals with the adaptation of quartz-microfabrics to changing physical deformation conditions, and discusses their preservation potential during subsequent retrograde deformation. Using microstructural analysis, a sequence of recrystallization processes in quartz, ranging from Grain-Boundary Migration Recrystallization (GBM) over Subgrain-Rotation Recrystallization (SGR) to Bulging Nucleation (BLG) is detected for the Simplon fault zone (SFZ) from the low strain rim towards the internal high strain part of the large-scale shear zone. Based on: (i) the retrograde cooling path; (ii) estimates of deformation temperatures; and (iii) spatial variation of dynamic recrystallization processes and different microstructural characteristics, continuous strain localization with decreasing temperature is inferred. In contrast to the recrystallization microstructures, crystallographic preferred orientations (CPO) have a longer memory. CPO patterns indicative of prism and rhomb glide systems in mylonitic quartz veins, overprinted at low temperatures (�400 �C), suggest inheritance of a high-temperature deformation. In this way, microstructural, textural and geochemical analyses provide information for several million years of the deformation history. The reasons for such incomplete resetting of the rock texture is that strain localization is caused by change in effective viscosity contrasts related to temporal large- and small-scale temperature changes during the evolution of such a long-lived shear zone. The spatially resolved, quantitative investigation of quartz microfabrics and associated recrystallization processes therefore provide great potential for an improved understanding of the geodynamics of large-scale shear zones.

Small quantity but large effect — How minor phases control strain localization in upper mantle shear zones

Low viscosity domains such as localized shear zones exert an important control on the geodynamics of the uppermost mantle. Grain size reduction and subsequent strain localization related to a switch from dislocation to diffusion creep is one mechanism to form low viscosity domains. To sustain strain localization, the grain size of mantle minerals needs to be kept small over geological timescales. One way to keep olivine grain sizes small is by pinning of mobile grain boundaries during grain growth by other minerals (second phases). Detailed microstructural studies based on natural samples from three shear zones formed at different geodynamic settings, allowed the derivation of the olivine grain-size dependence on the second-phase content. The polymineralic olivine grain-size evolution with increasing strain is similar in the three shear zones. If the second phases are to pin the mobile olivine grain boundary the phases need to be well mixed before grain growth. We suggest that melt-rock and metamorphic reactions are crucial for the initial phase mixing in mantle rocks. With ongoing deformation and increasing strain, grain boundary sliding combined with mass transfer processes and nucleation of grains promotes phase mixing resulting in fine-grained polymineralic mixtures that deform by diffusion creep. Strain localization due to the presence of volumetrically minor minerals in polymineralic mantle rocks is only important at high strain deformation (ultramylonites) at low temperatures (<~800°C). At smaller strain and stress conditions and/or higher temperatures other parameters like overall energy available to deform a given rock volume, the inheritance of mechanical anisotropies or the presence of water or melts needs to be considered to explain strain localization in the upper mantle.

Granular flow in polymineralic rocks bearing sheet silicates: New evidence from natural examples

Natural deformation in carbonate mylonites bearing sheet silicates occurs via a complex interaction of granular flow and solution transfer processes and involves continuous cycles of dissolution, grain boundary diffusion, nucleation and growth. In this way, new sheet silicates (a) nucleate within voids formed by grain boundary sliding of calcite grains. (b) grow, and (c) rotate towards the shear plane. As a consequence, small mica grains show a wide range of orientations with respect to the shear plane, but moderate to large grains are subparallel both to each other and to the shear plane. Increases of average grain sizes with increasing temperature of sheet silicates in mica-rich layers is more pronounced than in mica-poor layers. In the calcitic matrix however, sheet silicates can only grow via solution-precipitation and mass transfer processes. Therefore, the observed grain size variability indicates drastic differences in mass transfer behavior between the individual layers, which might be related to differences in the fluid flux. Based on these observations, a conceptual model for the microfabric evolution in sheet silicate bearing mylonites is presented. © 2001 Elsevier Science B.V. All rights reserved.

Temperature- and strain-rate-dependent microfabric evolution in monomineralic mylonite: evidence from in situ deformation of norcamphor

Norcamphor (C7H10O) was subjected to plane strain simple shear in a see-through deformation rig at four different strain rate and temperature conditions. Two transient stages in the microfabric evolution to steady state are distinguished. The grain scale mechanisms associated with the microstructural and textural evolution vary with the applied temperature, strain rate and strain. In high-temperature-low-strain-rate experiments, computer integrated polarization microscopy reveals that the texture evolution is closely related to the crystallographic rotation paths and rotation rates of individual grains. High c-axis rotation rates at low to intermediate shear strains are related to the development of a symmetrical c-axis cross girdle by the end of the first transient stage (γ = 1.5 to 2). During the second transient stage (γ = 1.5 to 6), the cross girdle yields to an oblique c-axis single girdle as c-axis rotation rates decrease and the relative activity of grain boundary migration recrystallization increases. Steady state (γ > 8) is characterized by a stable end orientation of the sample texture and the cyclic growth, rotation and consumption of individual grains within the aggregate.

The evolution of high-temperature mylonitic microfabrics: evidence from simple shearing of a quartz analogue (norcamphor)

Plane strain simple shearing of norcamphor (C7H10O) in a see-through deformation rig to a shear strain of γ = 10.5 at a homologous temperature of Th = 0.81 yields a microfabric similar to that of quartz in amphibolite facies mylonite. Synkinematic analysis of the norcamphor microfabric reveals that the development of a steady-state texture is linked to changes in the relative activities of several grain-scale mechanisms. Three stages of textural and microstructural evolution are distinguished: (1) rotation and shearing of the intracrystalline glide planes are accommodated by localized deformation along three sets of anastomozing microshears. A symmetrical c-axis girdle reflects localized pure shear extension along the main microshear set (Sa) oblique to the bulk shear zone boundary (abbreviated as SZB); (2) progressive rotation of the microshears into parallelism with the SZB increases the component of simple shear on the Sa microshears. Grain-boundary migration recrystallization favours the survival of grains with slip systems oriented for easy glide. This is associated with a textural transition towards two stable c-axis point maxima whose skeletal outline is oblique with respect to the Sa microshears and the SZB; and (3) at high shear strains (γ > 8), the microstructure, texture and mechanism assemblage are strain invariant, but strain continues to partition into rotating sets of microshears. Steady state is therefore a dynamic, heterogeneous condition involving the cyclic nucleation, growth and consumption of grains.

Deformation mechanisms in second-phase affected microstructures and their energy balance

Based on the relationship Zener parameter (Z=second-phase size/second-phase volume fraction) vs. calcite grain size (dg), second-phase controlled aggregates and microstructures that are weakly affected by second-phases are discriminated. The latter are characterized by large but constant grain sizes, high calcite grain boundary fractions and crystallographic preferred orientations (CPO), while calcite grain size and calcite grain boundary fraction decrease continuously and CPO weakens with decreasing Z in second-phase controlled microstructures. These observations suggest that second-phase controlled microstructures predominantly deform via granular flow because pinning of calcite grain boundaries reduces the efficiency of dynamic recrystallization favoring mass transfer processes and grain boundary sliding. In contrast, the balance of grain size reduction and growth by dynamic recrystallization maintains a steady state grain size in microstructures that are only weakly affected by second-phases promoting a predominance of dislocation creep. With increasing temperature, the relationship between Z and dg persists but the calcite grain size increases continuously. Based on microstructures, the energy of each modifying process is calculated and its relative contribution is compared with energies of the competing processes (surface energy, dragging energy, dynamic recrystallization energy). The steady state microstructures result from a temperature-dependent energy minimization procedure of the system.

The effect of dissolved magnesium on diffusion creep in calcite

We experimentally tested a series of synthetic calcite marbles with varying amounts of dissolved magnesium in a standard triaxial deformation machine at 300 MPa confining pressure, temperatures between 700 and 850°C, stresses between 2 and 100 MPa, and strain rates between 10−7 and 10−3 s−1. The samples were fabricated by hot isostatic pressing of a mixture of calcite and dolomite at 850°C and 300 MPa. The fabrication protocol resulted in a homogeneous, fine-grained high-magnesian calcite aggregate with minimal porosity and with magnesium contents between 0.07 and 0.17 mol% MgCO3. At stresses below 40 MPa the samples deformed with linear viscosity that depended inversely on grain size to the 3.26±0.51 power, suggesting that the mechanisms of deformation were some combination of grain boundary diffusion and grain boundary sliding. Because small grain sizes tended to occur in the high-magnesium calcite, the strength also appeared to vary inversely with magnesium content. However, the strength at constant grain size does not depend on the amount of dissolved magnesium, and thus, the impurity effect seems to be indirect. At stresses higher than 40 MPa, the aggregates become non-linearly viscous, a regime we interpret to be dislocation creep. The transition between the two regimes depends on grain size, as expected. The activation energy for diffusion creep is 200±30 kJ/mol and is quite similar to previous measurements in natural and synthetic marbles deformed at similar conditions with no added magnesium.

The influence of nano-scale second-phase particles on deformation of fine grained calcite mylonites

Grey and white carbonate mylonites were collected along thrust planes of the Helvetic Alps. They are characterised by very small grain sizes and non-random grain shape (SPO) and crystallographic preferred orientation (CPO). Presumably they deformed in the field of grain size sensitive flow by recrystallisation accommodated intracrystalline deformation in combination with granular flow. Both mylonites show a similar mean grain size, but in the grey mylonites the grain size range is larger, the grain shapes are more elongate and the dynamically recrystallised calcite grains are more often twinned. Grey mylonites have an oblique CPO, while the CPO in white mylonites is symmetric with respect to the shear plane. Combustion analysis and TEM investigations revealed that grey mylonites contain a higher amount of highly structured kerogens with particle sizes of a few tens of nanometers, which are finely dispersed at the grain boundaries. During deformation of the rock, nano-scale particles reduced the migration velocity of grain boundaries by Zener drag resulting in slower recrystallisation rates of the calcite aggregate. In the grey mylonites, more strain increments were accommodated by individual grains before they became refreshed by dynamic recrystallisation than in white mylonites, where grain boundary migration was less hindered and recrystallisation cycles were faster. Consequently, grey mylonites represent ‘deformation’ microfabrics while white mylonites are characterised by ‘recrystallisation’ microfabrics. Field geologists must utilise this different deformation behavior when applying the obliquity in CPO and SPO of the respective mylonites as reliable shear sense indicators.

Evolution of mylonitic microfabrics (EMM), a computer application for educational purposes

`Evolution of mylonitic microfabrics' (EMM) is an interactive Filemaker Pro 3.0 application that documents a series of see-through deformation experiments on polycrystalline norcamphor. The application comprises computer animations, graphics and text explanations designed to give students and researchers insight into the interaction and dynamic nature of small-scale, mylonitic processes like intracrystalline glide, dynamic recrystallization and strain localization (microshearing). EMM shows how mylonitic steady state is achieved at different strain rates and temperatures. First, rotational mechanisms like glide-induced vorticity, subgrain rotation recrystallization and rigid-body rotation bring grains' crystal lattices into orientations that are favorable for intracrystalline glide. In a second stage, selective elimination of grains whose lattices are poorly oriented for glide involves grain boundary migration. This strengthens the texture. Temperature and strain rate affect both the relative activity of different strain accommodation mechanisms and the rate of microfabric change. Steady-state microfabrics are characterized by stable texture, grain size and shape-preferred orientations of grains and domains. This involves the cyclical generation and elimination of dynamically recrystallized grains and microshear zones.

A new technique to automatically quantify microstructures of fine grained carbonate mylonites: two-step etching combined with SEM imaging and image analysis

A two-step etching technique for fine-grained calcite mylonites using 0.37% hydrochloric and 0.1% acetic acid produces a topographic relief which reflects the grain boundary geometry. With this technique, calcite grain boundaries become more intensely dissolved than their grain interiors but second phase minerals like dolomite, quartz, feldspars, apatite, hematite and pyrite are not affected by the acid and therefore form topographic peaks. Based on digital backscatter electron images and element distribution maps acquired on a scanning electron microscope, the geometry of calcite and the second phase minerals can be automatically quantified using image analysis software. For research on fine-grained carbonate rocks (e.g. dolomite calcite mixtures), this low-cost approach is an attractive alternative to the generation of manual grain boundary maps based on photographs from ultra-thin sections or orientation contrast images.