Grain coarsening in contact metamorphic carbonates: effects of second-phase particles, fluid flow and thermal perturbations

Under contact metamorphic conditions, carbonate rocks in the direct vicinity of the Adamello pluton reflect a temperature-induced grain coarsening. Despite this large-scale trend, a considerable grain size scatter occurs on the outcrop-scale indicating local influence of second-order effects such as thermal perturbations, fluid flow and second-phase particles. Second-phase particles, whose sizes range from nano- to the micron-scale, induce the most pronounced data scatter resulting in grain sizes too small by up to a factor of 10, compared with theoretical grain growth in a pure system. Such values are restricted to relatively impure samples consisting of up to 10 vol.% micron-scale second-phase particles, or to samples containing a large number of nano-scale particles. The obtained data set suggests that the second phases induce a temperature-controlled reduction on calcite grain growth. The mean calcite grain size can therefore be expressed in the form D 1⁄4 C2 eQ*/RT(dp/fp)m*, where C2 is a constant, Q* is an activation energy, T the temperature and m* the exponent of the ratio dp/fp, i.e. of the average size of the second phases divided by their volume fraction. However, more data are needed to obtain reliable values for C2 and Q*. Besides variations in the average grain size, the presence of second-phase particles generates crystal size distribution (CSD) shapes characterized by lognormal distributions, which differ from the Gaussian-type distributions of the pure samples. In contrast, fluid-enhanced grain growth does not change the shape of the CSDs, but due to enhanced transport properties, the average grain sizes increase by a factor of 2 and the variance of the distribution increases. Stable d18O and d13C isotope ratios in fluid-affected zones only deviate slightly from the host rock values, suggesting low fluid/rock ratios. Grain growth modelling indicates that the fluid-induced grain size variations can develop within several ka. As inferred from a combination of thermal and grain growth modelling, dykes with widths of up to 1 m have only a restricted influence on grain size deviations smaller than a factor of 1.1.To summarize, considerable grain size variations of up to one order of magnitude can locally result from second-order effects. Such effects require special attention when comparing experimentally derived grain growth kinetics with field studies.

Differences in grain growth of calcite: a field-based modeling approach

Normal grain growth of calcite was investigated by combining grain size analysis of calcite across the contact aureole of the Adamello pluton, and grain growth modeling based on a thermal model of the surroundings of the pluton. In an unbiased model system, i.e., location dependent variations in temperature-time path, 2/3 and 1/3 of grain growth occurs during pro- and retrograde metamorphism at all locations, respectively. In contrast to this idealized situation, in the field example three groups can be distinguished, which are characterized by variations in their grain size versus temperature relationships: Group I occurs at low temperatures and the grain size remains constant because nano-scale second phase particles of organic origin inhibit grain growth in the calcite aggregates under these conditions. In the presence of an aqueous fluid, these second phases decay at a temperature of about 350 °C enabling the onset of grain growth in calcite. In the following growth period, fluid-enhanced group II and slower group III growth occurs. For group II a continuous and intense grain size increase with T is typical while the grain growth decreases with T for group III. None of the observed trends correlate with experimentally based grain growth kinetics, probably due to differences between nature and experiment which have not yet been investigated (e.g., porosity, second phases). Therefore, grain growth modeling was used to iteratively improve the correlation between measured and modeled grain sizes by optimizing activation energy (Q), pre-exponential factor (k0) and grain size exponent (n). For n=2, Q of 350 kJ/mol, k0 of 1.7×1021 μmns−1 and Q of 35 kJ/mol, k0 of 2.5×10-5 μmns−1 were obtained for group II and III, respectively. With respect to future work, field-data based grain growth modeling might be a promising tool for investigating the influences of secondary effects like porosity and second phases on grain growth in nature, and to unravel differences between nature and experiment.

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.

Orbital elements of the material surrounding comet 67P/Churyumov-Gerasimenko

Context. We investigate the dust coma within the Hill sphere of comet 67P/Churyumov-Gerasimenko. Aims. We aim to determine osculating orbital elements for individual distinguishable but unresolved slow-moving grains in the vicinity of the nucleus. In addition, we perform photometry and constrain grain sizes. Methods. We performed astrometry and photometry using images acquired by the OSIRIS Wide Angle Camera on the European Space Agency spacecraft Rosetta. Based on these measurements, we employed standard orbit determination and orbit improvement techniques. Results. Orbital elements and effective diameters of four grains were constrained, but we were unable to uniquely determine them. Two of the grains have light curves that indicate grain rotation. Conclusions. The four grains have diameters nominally in the range 0.14-0.50 m. For three of the grains, we found elliptic orbits, which is consistent with a cloud of bound particles around the nucleus. However, hyperbolic escape trajectories cannot be excluded for any of the grains, and for one grain this is the only known option. One grain may have originated from the surface shortly before observation. These results have possible implications for the understanding of the dispersal of the cloud of bound debris around comet nuclei, as well as for understanding the ejection of large grains far from the Sun.

On the initiation of boudinage and folding instabilities in layered elasto-visco-plastic solids - Insights from natural microstructures and numerical simulations

The present understanding of the initiation of boudinage and folding structures is based on viscosity contrasts and stress exponents, considering an intrinsically unstable state of the layer. The criterion of localization is believed to be prescribed by geometry-material interactions, which are often encountered in natural structures. An alternative localization phenomenon has been established for ductile materials, in which instability emerges for critical material parameters and loading rates from homogeneous conditions. In this thesis, conditions are sought under which this type of instability prevails and whether localization in geological materials necessarily requires a trigger by geometric imperfections. The relevance of critical deformation conditions, material parameters and the spatial configuration of instabilities are discussed in a geological context. In order to analyze boudinage geometries, a numerical eigenmode analysis is introduced. This method allows determining natural frequencies and wavelengths of a structure and inducing perturbations on these frequencies. In the subsequent coupled thermo-mechanical simulations, using a grain size evolution and end-member flow laws, localization emerges when material softening through grain size sensitive viscous creep sets in. Pinch-and-swell structures evolve along slip lines through a positive feedback between the matrix response and material bifurcations inside the layer, independent from the mesh-discretization length scale. Since boudinage and folding are considered to express the same general instability, both structures should arise independently of the sign of the loading conditions and for identical material parameters. To this end, the link between material to energy instabilities is approached by means of bifurcation analyses of the field equations and finite element simulations of the coupled system of equations. Boudinage and folding structures develop at the same critical energy threshold, where dissipative work by temperature-sensitive creep overcomes the diffusive capacity of the layer. This finding provides basis for a unified theory for strain localization in layered ductile materials. The numerical simulations are compared to natural pinch-and-swell microstructures, tracing the adaption of grain sizes, textures and creep mechanisms in calcite veins. The switch from dislocation to diffusion creep relates to strain-rate weakening, which is induced by dissipated heat from grain size reduction, and marks the onset of continuous necking. The time-dependent sequence uncovers multiple steady states at different time intervals. Microstructurally and mechanically stable conditions are finally expressed in the pinch-and-swell end members. The major outcome of this study is that boudinage and folding can be described as the same coupled energy-mechanical bifurcation, or as one critical energy attractor. This finding allows the derivation of critical deformation conditions and fundamental material parameters directly from localized structures in the field.

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.