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.

Linking temperature estimates and microstructures in deformed polymineralic mantle rocks

To constrain deformation temperatures of mantle shear zones, we studied a strike-slip shear zone (Hilti massif, Semail ophiolite, Oman) and focused on the interaction between microstructural mechanisms and chemical equilibration processes. Quantitative microfabric analysis on harzburgites with different deformation intensity (porphyroclastic tectonite, mylonite, and ultramylonite) was combined with orthopyroxene geothermometry. The average grain size of all phases decreases with decreasing shear zone thickness. Dynamic recrystallization of porphyroclasts in combination with dissolution-precipitation and nucleation result in small-sized, chemically equilibrated pyroxenes. The composition of orthopyroxene was used to calculate deformation temperatures. In the case of the porphyroclastic tectonites, the chemical composition of orthopyroxene has been reset by diffusion yielding temperature estimates of 880-900 degrees C. The mylonites were deformed by dislocation creep of olivine and show a broad range of calculated temperatures, which result from a combination of grain size reduction and inheritance of equilibrium compositions from earlier high-temperature events and diffusion. In mylonites, diffusion profiles combined with geothermometry and grain size analysis indicate a mylonitic deformation temperature of 800-900 degrees C possibly followed by diffusion. In ultramylonites, the smallest grains (<30 mu m) reveal equilibration at temperatures of similar to 700 degrees C during the last stages of ductile deformation, which was dominated by diffusion creep of olivine. Our results provide a crucial link between temperature and evolution of microstructures from dislocation creep to diffusion creep in mantle shear zones.

Fluid-rock-strain feedbacks in extensional shear zones

La présence de fluide météorique synchrone à l'activité du détachement (Farmin, 2003 ; Mulch et al., 2007 ; Gébelin et al., 2011), implique que les zones de cisaillement sont des systèmes ouverts avec des cellules de convections à l'échelle crustale et un intense gradient géothermique au sein du détachement (Morrison et Anderson, 1998, Gottardi et al., 2011). De plus, les réactions métamorphiques liées à des infiltrations fluides dans les zones de cisaillement extensionnel peuvent influencer les paramètres rhéologiques du système (White and Knipe, 1978), et impliquer la localisation de la déformation dans la croûte. Dans ce manuscrit, deux zones de cisaillement infiltrées par des fluides météoriques sont étudiées, l'une étant largement quartzitique, et l'autre de nature granitique ; les relations entre déformation, fluides, et roches s'appuient sur des approches structurales, microstructurales, chimiques et isotopiques. L'étude du détachement du Columbia river (WA, USA) met en évidence que la déformation mylonitique se développe en un million d'années. La phase de cisaillement principal s'effectue à 365± 30°C d'après les compositions isotopiques en oxygène du quartz et de la muscovite. Ces minéraux atteignent l'équilibre isotopique lors de leur recristallisation dynamique contemporaine à la déformation. La zone de cisaillement enregistre une baisse de température, remplaçant le mécanisme de glissement par dislocation par celui de dissolution- précipitation dans les derniers stades de l'activité du détachement. La dynamique de circulation fluide bascule d'une circulation pervasive à chenalisée, ce qui engendre localement la rupture des équilibres d'échange isotopiques. La zone de cisaillement de Bitterroot (MT, USA) présente une zone mylonitique de 600m d'épaisseur, progressant des protomylonites aux ultramylonites. L'intensité de la localisation de la déformation se reflète directement sur l'hydratation des feldspaths, réaction métamorphique majeure dite de « rock softening ». Une étude sur roche totale indique des transferts de masse latéraux au sein des mylonites, et d'importantes pertes de volume dans les ultramylonites. La composition isotopique en hydrogène des phyllosilicates met en évidence la présence (1) d'une source magmatique/métamorphique originelle, caractérisée par les granodiorites ayant conservé leur foliation magmatique, jusqu'aux protomylonites, et (2) une source météorique qui tamponne les valeurs des phyllosilicates des fabriques mylonitiques jusqu'aux veines de quartz non-déformées. Les compositions isotopiques en oxygène des minéraux illustrent le tamponnement de la composition du fluide météorique par l'encaissant. Ce phénomène cesse lors du processus de chloritisation de la biotite, puisque les valeurs des chlorites sont extrêmement négatives (-10 per mil). La thermométrie isotopique indique une température d'équilibre isotopique de la granodiorite entre 600-500°C, entre 500-300°C dans les mylonites, et entre 300 et 200°C dans les fabriques cassantes (cataclasites et veines de quartz). Basé sur les résultats issus de ce travail, nous proposons un modèle général d'interactions fluide-roches-déformation dans les zones de détachements infiltrées par des fluides météoriques. Les zones de détachements évoluent rapidement (en quelques millions d'années) au travers de la transition fragile-ductile ; celle-ci étant partiellement contrôlée par l'effet thermique des circulations de fluide météoriques. Les systèmes de détachements sont des lieux où la déformation et les circulations fluides sont couplées ; évoluant rapidement vers une localisation de la déformation, et de ce fait, une exhumation efficace. - The presence of meteoric fluids synchronous with the activity of extensional detachment zones (Famin, 2004; Mulch et al., 2007; Gébelin et al., 2011) implies that extensional systems involve fluid convection at a crustal scale, which results in high geothermal gradients within active detachment zones (Morrison and Anderson, 1998, Gottardi et al., 2011). In addition, the metamorphic reactions related to fluid infiltration in extensional shear zones can influence the rheology of the system (White and Knipe, 1978) and ultimately how strain localizes in the crust. In this thesis, two shear zones that were permeated by meteoric fluids are studied, one quartzite-dominated, and the other of granitic composition; the relations between strain, fluid, and evolving rock composition are addressed using structural, microstructural, and chemical/isotopic measurements. The study of the Columbia River detachment that bounds the Kettle core complex (Washington, USA) demonstrates that the mylonitic fabrics in the 100 m thick quartzite- dominated detachment footwall developed within one million years. The main shearing stage occurred at 365 ± 30°C when oxygen isotopes of quartz and muscovite equilibrated owing to coeval deformation and dynamic recrystallization of these minerals. The detachment shear zone records a decrease in temperature, and dislocation creep during detachment shearing gave way to dissolution-precipitation and fracturing in the later stages of detachment activity. Fluid flow switched from pervasive to channelized, leading to isotopic disequilibrium between different minerals. The Bitterroot shear zone detachment (Montana, USA) developed a 600 m thick mylonite zone, with well-developed transitions from protomylonite to ultramylonite. The localization of deformation relates directly to the intensity of feldspar hydration, a major rock- softening metamorphic reaction. Bulk-rock analyses of the mylonitic series indicate lateral mass transfer in the mylonite (no volume change), and significant volume loss in ultramylonite. The hydrogen isotope composition of phyllosilicates shows (1) the presence of an initial magmatic/metamorphic source characterized by the granodiorite in which a magmatic, and gneissic (protomylonite) foliation developed, and (2) a meteoric source that buffers the values of phyllosilicates in mylonite, ultramylonite, cataclasite, and deformed and undeformed quartz veins. The mineral oxygen isotope compositions were buffered by the host-rock compositions until chloritization of biotite started; the chlorite oxygen isotope values are negative (-10 per mil). Isotope thermometry indicates a temperature of isotopic equilibrium of the granodiorite between 600-500°C, between 500-300°C in the mylonite, and between 300 and 200°C for brittle fabrics (cataclasite and quartz veins). Results from this work suggest a general model for fluid-rock-strain feedbacks in detachment systems that are permeated by meteoric fluids. Phyllosilicates have preserved in their hydrogen isotope values evidence for the interaction between rock and meteoric fluids during mylonite development. Fluid flow generates mass transfer along the tectonic anisotropy, and mylonites do not undergo significant volume change, except locally in ultramylonite zones. Hydration of detachment shear zones attends mechanical grain size reduction and enhances strain softening and localization. Self-exhuming detachment shear zones evolve rapidly (a few million years) through the transition from ductile to brittle, which is partly controlled by the thermal effect of circulating surface fluids. Detachment systems are zones in the crust where strain and fluid flow are coupled; these systems. evolve rapidly toward strain localization and therefore efficient exhumation.

Evolution of a polymineralic mantle shear zone and the role of second phases in the localization of deformation

The influence of second phases (e.g., pyroxenes) on olivine grain size was studied by quantitative microfabric analyses of samples of the Hilti massif mantle shear zone (Semail ophiolite, Oman). The microstructures range from porphyroclastic tectonites to ultramylonites, from outside to the center of the shear zone. Starting at conditions of ridge-related flow, they formed under continuous cooling leading to progressive strain localization. The dependence of the average olivine grain size on the second-phase content can be split into a second-phase controlled and a dynamic recrystallization-controlled field. In the former, the olivine grain size is related to the ratio between the second-phase grain size and volume fraction (Zener parameter). In the latter, dynamic recrystallization manifested by a balance between grain growth and grain size reduction processes yields a stable olivine grain size. In both fields the average olivine and second-phase grain size decreases with decreasing temperature. Combining the microstructural information with deformation mechanism maps suggests that the porphyroclastic tectonites (similar to 1100 degrees C) and mylonites (similar to 800 degrees C) formed under the predominance of dislocation creep. Since olivine-rich layers are intercalated with layer parallel, polymineralic bands in the mylonites, nearly equiviscous conditions can be assumed. In the ultramylonites, diffusion creep represents the major deformation mechanism in the polymineralic layers. It is this switch in deformation mechanism from dislocation creep to diffusion creep that forces strain to localize in the fine-grained polymineralic domains at low temperatures (<similar to 700 degrees C), underlining the role of the second phases on strain localization in cooling mantle rocks.

Structural analysis of the Paint Lake Deformation Zone, Northern Ontario /

The Paint Lake Deformation Zone (PLDZ), located within the Superior Province of Canada, demarcates a major structural and lithological break between the Onaman-Tashota Terrane to the north and the Beardmore-Geraldton Belt to the south. The PLDZ is an east-west trending lineament, approximately 50 km in length and up to 1 km in width, comprised of an early ductile component termed the Paint Lake Shear Zone and a late brittle component known as the Paint Lake Fault. Structures associated with PLDZ development including S-, C- and C'-fabrics, stretching lineations, slickensides, C-C' intersection lineations, Z-folds and kinkbands indicate that simple shear deformation dominated during a NW-SE compressional event. Movement along the PLDZ was in a dextral sense consisting of an early differential motion with southside- down and a later strike-slip motion. Although the locus of the PLDZ may in part be lithologically controlled, mylonitization which accompanied shear zone development is not dependent on the lithological type. Conglomerate, intermediate and mafic volcanic units exhibit similar mesoscopic and microscopic structures where transected by the PLDZ. Field mapping, supported by thin section analysis, defines five strain domains increasing in intensity of deformation from shear zone boundary to centre. A change in the dominant microstructural deformation mechanism from dislocation creep to diffusion creep is observed with increasing strain during mylonitization. C'-fabric development is temporally associated with this change. A decrease in the angular relationship between C- and C'-fabrics is observed upon attaining maximum strain intensity. Strain profiling of the PLDZ demonstrates the presence of an outer primary strain gradient which exhibits a simple profile and an inner secondary strain gradient which exhibits a more complex profile. Regionally metamorphosed lithologies of lower greenschist facies outside the PLDZ were subjected to retrograde metamorphism during deformation within the PLDZ.

Quartz recrystallization regimes, c-axis texture transitions and fluid inclusion reequilibration in a prograde greenschist to amphibolite facies mylonite zone (Ribeira Shear Zone, SE Brazil)

The crystal-plastic behavior of quartz mylonites from the Ribeira Shear Zone (SE Brazil), a major strike-slip structure that was active during a prograde metamorphic phase related to the Neoproterozoic Brasiliano-Pan African Orogeny, was investigated using a multi-method approach. Geothermobarometry results indicate deformational conditions ranging from similar to 300 to similar to 630 degrees C and 500-700 MPa. A strong correlation between mapped metamorphic zones and a dominance of different dynamic recrystallization mechanisms of quartz occurs within the mylonite zone. Bulging recrystallization (BLG) dominates within the chlorite zone between 300 and 410 degrees C, subgrain rotation recrystallization (SGR) operates within the biotite zone from 410 to 520 degrees C, and grain boundary migration recrystallization (GBM) dominates in the garnet zone above 520 degrees C. The development of quartz c-axis textures is mainly governed by temperature and dynamic recrystallization mechanisms. Textures from BLG zone mylonites are characterized by maxima around Z; SGR zone mylonites display single girdles or asymmetric type I crossed girdles; and GBM zone mylonites comprise maxima around Y and intermediate between X and Z. The scarcity or absence of water-bearing fluid inclusions in quartz mylonites from the SGR and GBM zones, which are dominated by carbonic inclusions, suggests water-deficient conditions, whereas BLG zone mylonites are dominated by water-bearing inclusions. This evidence indicates that water was available in the protoliths but has been eliminated with increasing deformation and deformation temperature. No effect of the water content variation on the quartz microstructural and recrystallized grain size evolution was detected, and little influence on c-axis texture development was observed. Most of the fluid inclusion densities were reequilibrated during the shear zone exhumation history, recording a decompression in the range of 300-500 MPa, while microstructural reequilibration effects related to the prograde metamorphism are largely preserved. Fluid inclusion microstructures and densities from two SGR zone samples preserved evidence for a near isothermal compression within the interior of the Ribeira Shear Zone during the prograde metamorphism. (C) 2009 Elsevier B.V. All rights reserved.

A microstructure based model of the deformation mechanisms and flow stress during elevated temperature straining of a magnesium alloy

We show that the variation of flow stress with strain rate and grain size in a magnesium alloy deformed at a constant strain rate and 450 °C can be predicted by a crystal plasticity model that includes grain boundary sliding and diffusion. The model predicts the grain size dependence of the critical strain rate that will cause a transition in deformation mechanism from dislocation creep to grain boundary sliding, and yields estimates for grain boundary fluidity and diffusivity.

From transient to steady state deformation and grain size: A thermodynamic approach using elasto-visco-plastic numerical modeling

Numerical simulation experiments give insight into the evolving energy partitioning during high-strain torsion experiments of calcite. Our numerical experiments are designed to derive a generic macroscopic grain size sensitive flow law capable of describing the full evolution from the transient regime to steady state. The transient regime is crucial for understanding the importance of micro structural processes that may lead to strain localization phenomena in deforming materials. This is particularly important in geological and geodynamic applications where the phenomenon of strain localization happens outside the time frame that can be observed under controlled laboratory conditions. Ourmethod is based on an extension of the paleowattmeter approach to the transient regime. We add an empirical hardening law using the Ramberg-Osgood approximation and assess the experiments by an evolution test function of stored over dissipated energy (lambda factor). Parameter studies of, strain hardening, dislocation creep parameter, strain rates, temperature, and lambda factor as well asmesh sensitivity are presented to explore the sensitivity of the newly derived transient/steady state flow law. Our analysis can be seen as one of the first steps in a hybrid computational-laboratory-field modeling workflow. The analysis could be improved through independent verifications by thermographic analysis in physical laboratory experiments to independently assess lambda factor evolution under laboratory conditions.

Results from Purgatory Conglomerate investigations

Quartz Crystallographic Preferred Orientation (CPO) patterns are most commonly a result of deformation by dislocation creep. We investigated whether Dissolution-Precipitation Creep (DPC) a process that occur at lower differential stresses and temperatures, may result in CPO in quartz. Within the Purgatory Conglomerate, DPC led to quartz dissolution along cobble surfaces perpendicular to the shortening direction, and quartz precipitation in overgrowths at the ends of the cobbles (strain shadows), parallel to the maximum extension direction. The Purgatory Conglomerate is part of the SE Narragansett basin where strain intensity increases from west to east and is associated with top-to-the-west transport and folding during the Alleghanian orogeny. Quartz c-axis orientations as revealed by Electron Backscatter Diffraction (EBSD) methods, were random in all analyzed domains within the cobbles and strain shadows irrespective of the intensity of strain or metamorphic grade of the sample. Quartz dissolution probably occurred exclusively along the cobbles' margins, leaving the remaining grains unaffected by DPC. The fact that quartz precipitated in random orientations may indicate that the strain shadows were regions of little or no differential stress.

From point defects to plate tectonic faults

Understanding and explaining emergent constitutive laws in the multi-scale evolution from point defects, dislocations and two-dimensional defects to plate tectonic scales is an arduous challenge in condensed matter physics. The Earth appears to be the only planet known to have developed stable plate tectonics as a means to get rid of its heat. The emergence of plate tectonics out of mantle convection appears to rely intrinsically on the capacity to form extremely weak faults in the top 100 km of the planet. These faults have a memory of at least several hundred millions of years, yet they appear to rely on the effects of water on line defects. This important phenomenon was first discovered in laboratory and dubbed ``hydrolytic weakening''. At the large scale it explains cycles of co-located resurgence of plate generation and consumption (the Wilson cycle), but the exact physics underlying the process itself and the enormous spanning of scales still remains unclear. We present an attempt to use the multi-scale non-equilibrium thermodynamic energy evolution inside the deforming lithosphere to move phenomenological laws to laws derived from basic scaling quantities, develop self-consistent weakening laws at lithospheric scale and give a fully coupled deformation-weakening constitutive framework. At meso- to plate scale we encounter in a stepwise manner three basic domains governed by the diffusion/reaction time scales of grain growth, thermal diffusion and finally water mobility through point defects in the crystalline lattice. The latter process governs the planetary scale and controls the stability of its heat transfer mode.

Effect of water on olivine single crystals plasticity, deformed under high pressure and high temperature

The Earth's upper mantle, mainly composed of olivine, is seismically anisotropic. Seismic anisotropy attenuation has been observed at 220km depth. Karato et al. (1992) attributed this attenuation to a transition between two deformation mechanisms, from dislocation creep above 220km to diffusion creep below 220km, induced by a change in water content. Couvy (2005) and Mainprice et al. (2005) predicted a change in Lattice Preferred Orientation induced by pressure, which comes from a change of slip system, from [100] slip to [001] slip, and is responsible for the seismic anisotropy attenuation. Raterron et al. (2007) ran single crystal deformation experiments under anhydrous conditions and observed that the slip system transition occurs around 8GPa, which corresponds to a depth of 260Km. Experiments were done to quantify the effects of water on olivine single crystals deformed using D-DIA press and synchrotron beam. Deformations were carried out in uniaxial compression along [110]c, [011]c, and [101]c, crystallographic directions, at pressure ranging from 4 to 8GPa and temperature between 1373 and 1473K. Talc sleeves about the annulus of the single crystals were used as source of water in the assembly. Stress and specimen strain rates were calculated by in-situ X-ray diffraction and time resolved imaging, respectively. By direct comparison of single crystals strain rates, we observed that [110]c deforms faster than [011]c below 5GPa. However above 6GPa [011]c deforms faster than [110]c. This revealed that [100](010) is the dominant slip system below 5GPa, and above 6GPa [001](010) becomes dominant. According to our results, the slip system transition, which is induced by pressure, occurs at 6GPa. Water influences the pressure where the switch over occurs, by lowering the transition pressure. The pressure effect on the slip systems activity has been quantified and the hydrolytic weakening has also been estimated for both orientations. Data also shows that temperature affects the slip system activity. The regional variation of the depth for the seismic anisotropy attenuation, which would depend on local hydroxyl content and temperature variations and explains the seismic anisotropy attenuation occurring at about 220Km depth in the mantle, where the pressure is about 6GPa. Deformation of MgO single crystal oriented [100], [110] and [111] were also performed. The results predict a change in the slip system activity at 23GPa, again induced by pressure. This explains the seismic anisotropy observed in the lower mantle.

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.

Dislocation jamming and Andrade creep

We simulate the glide motion of an assembly of interacting dislocations under the action of an external shear stress and show that the associated plastic creep relaxation follows Andrades law. Our results indicate that Andrade creep in plastically deforming crystals involves the correlated motion of dislocation structures near a dynamic transition separating a flowing from a jammed phase. Simulations in the presence of dislocation multiplication and noise confirm the robustness of this finding and highlight the importance of metastable structure formation for the relaxation process.

The creep behaviour of pressure diecast zinc-aluminium based alloys

The creep behaviour of three pressure diecast commercial zinc-aluminium based alloys: Mazak 3, corresponding to BS 1004A, and the new alloys ZA.8 and ZA.27 with a series of alloys with compositions ranging from 0% to 30% aluminium was investigated. The total creep elongation of commercial alloys was shown to be well correlated using an empirical equation. Based on this a parametrical relationship was derived which allowed the total creep extension to be related to the applied stress, the temperature and the time of test, so that a quantitative assessment of creep of the alloys could be made under different conditions. Deviation from the normal creep kinetics occurred in alloys ZA.8 and ZA.27 at very low stresses, 150°C, due to structural coarsening combined with partial transformation of ε -phase into T' phase. The extent of primary creep was found to increase with aluminium content, but secondary creep rates decreased in the order Mazak 3, ZA.8 and ZA.27. Thus, based on the above equation, ZA.8 was found to have a substantially better total creep resistance than ZA.27, which in turn was marginally better than Mazak 3 for strains higher than 0.5%, but inferior for smaller strains, due to its higher primary creep extension. The superior creep resistance of ZA.8 was found to be due to the presence of strictly-orientated, thin plate-like precipitates of ε(CuZn4) phase in the zinc matrix of the eutectic and the lamellarly decomposed β phase, in which the precipitation morphology and orientation of ε in the zinc matrix was determined. Over broad ranges of temperature and stresses, the stress exponents and activation energies for creep were found to be consistent with some proposed creep rate mechanisms; i.e. viscous glide for Mazak 3, dislocation climb over second phase particles for ZA.8 and dislocation climb for ZA.27, controlled by diffusion in the zinc-rich phase. The morphology of aluminium and copper-rich precipitates formed from the solid solution of zinc was clearly revealed. The former were found to further increase the creep rate of inherently low creep resistant zinc, but the latter contributed significantly to the creep resistance. Excess copper in the composition, however, was not beneficial in improving the creep resistance. Decomposition of β in copper-containing alloys was found to be through a metastable Zn-Al phase which is strongly stabilised by copper, and the final products of the decomposition had a profound effect on the creep strength of the alloys. The poor creep resistance of alloy ZA.27 was due to the presence of particulate products derived from decomposed β-phase and a large volume of fine, equiaxed products of continuously decomposed α-dendrites.