Compositional zoning of calcite in a high-pressure metamorphic calc-schist: clues to heterogeneous grain-scale fluid distribution during exhumation

Calcite in former aragonite–dolomite-bearing calc-schists from the ultrahigh-pressure metamorphic (UHPM) oceanic complex at Lago di Cignana, Valtournanche, Italy, preserved different kinds of zoning patterns at calcite grain and phase boundaries. These patterns are interpreted in terms of lattice diffusion and interfacial mass transport linked with a heterogeneous distribution of fluid and its response to a changing state of stress. The succession of events that occurred during exhumation is as follows: As the rocks entered the calcite stability field at T=530–550 °C, P ca. 1.2 GPa, aragonite occurring in the matrix and as inclusions in poikilitic garnet was completely transformed to calcite. Combined evidence from microstructures and digital element distribution maps (Mn-, Mg-, Fe- and Ca–Kα radiation intensity patterns) indicates that transformation rates have been much higher than rates of compositional equilibration of calcite (involving resorption of dolomite and grain boundary transport of Mg, Fe and Ca). This rendered the phase transformation an isochemical process. During subsequent cooling to T ca. 490 °C (where lattice diffusion effectively closed), grains of matrix calcite have developed diffusion-zoned rims, a few hundred micrometres thick, with Mg and Fe increasing and Ca decreasing towards the phase boundary. Composition profiles across concentrically zoned, large grains in geometrically simple surroundings can be successfully modelled with an error function describing diffusion into a semi-infinite medium from a source of constant composition. The diffusion rims in matrix calcite are continuous with quartz, phengite, paragonite and dolomite in the matrix. This points to an effective mass transport on phase boundaries over a distance of several hundred micrometres, if matrix dolomite has supplied the Mg and Fe needed for incorporation in calcite. In contrast, diffusion rims are lacking at calcite–calcite and most calcite–garnet boundaries, implying that only very minor mass transport has occurred on these interfaces over the same T–t interval. From available grain boundary diffusion data and experimentally determined fluid–solid grain boundary structures, inferred large differences in transport rates can be best explained by the discontinuous distribution of aqueous fluid along grain/phase boundaries. Observed patterns of diffusion zoning indicate that fluid was distributed not only along grain-edge channels, but spread out along most calcite–white mica and calcite–quartz two-grain junctions. On the other hand, the inferred non-wetting of calcite grain boundaries in carbonate-rich domains is compatible with fluid–calcite–calcite dihedral angles >60° determined by Holness and Graham (1995) for a wide range of fluid compositions under the P–T conditions of interest. Whereas differential stress has been very low at the stage of diffusion zoning (T > 490 °C), it increased as the rocks were cooling below 440 °C (at 0.3–0.5 GPa). Dislocation creep and the concomitant increase of strain energy in matrix calcite induced migration recrystallisation of high-angle grain boundaries. For that stage, the compositional microstructure of recrystallised calcite grain boundary domains indicates significant mass transport along calcite two-grain junctions, which at the established low temperatures is likely to have been accomplished by ionic diffusion within a hydrous grain boundary fluid film (“dynamic wetting” of migrating grain boundaries). Received: 10 January 2000 / Accepted: 10 April 2000     

Chemically induced grain boundary migration in calcite: temperature dependence,phenomenology, and possible applications to geologic systems

Chemically induced grain boundary migration (CIGM) is a solid-state reaction mechanism which involves grain boundary migration. Initially straight grain boundaries in calcite bicrystals exhibit CIGM when exposed to a SrCO₃ or BaCO₃-rich melt at temperatures of 680–840° C in a 1-bar CO₂ atmosphere. Under these conditions, CIGM is apparently a much more efficient way to form solid-solutions than lattice diffusion. The maximum observed rate of migration, 1×10^–9 m/s, occurred at about 760° C. Below that temperature, migration was thermally activated, but above it, the process was thermally inhibited. Some portions of the boundary swept through a given crystal volume more than once. In the (Ca, Sr) CO₃ system, the solute concentrations incorporated during a single pass of the boundary were far less than values expected at chemical equilibrium, and departed further from equilibrium at higher temperatures. CIGM may be geologically important when the distances characteristic of grain boundary migration and diffusion are much larger than those characteristic of lattice diffusion. The great similarity between results in ionic calcite and metals systems suggests CIGM is a very general process which may be expected to occur in a large number of geologic systems.     

Asymmetrically zoned reaction rims: assessment of grain boundary diffusivities and growth rates related to natural diffusion-controlled mineral reactions

This study explores garnet coronas around hedenbergite, which were formed by the reaction plagioclase + hedenbergite→garnet + quartz, to derive information about diffusion paths that allowed for material redistribution during reaction progress. Whereas quartz forms disconnected single grains along the garnet/hedenbergite boundaries, garnet forms ～20‐μm‐wide continuous polycrystalline rims along former plagioclase/hedenbergite phase boundaries. Individual garnet crystals are separated by low‐angle grain boundaries, which commonly form a direct link between the reaction interfaces of the plagioclase|garnet|hedenbergite succession. Compositional variations in garnet involve: (i) an overall asymmetric compositional zoning in Ca, Fe²⁺, Fe³⁺ and Al across the garnet layer; and (ii) micron‐scale compositional variations in the near‐grain boundary regions and along plagioclase/garnet phase boundaries. These compositional variations formed during garnet rim growth. Thereby, transfer of the chemical components occurred by a combination of fast‐path diffusion along grain boundaries within the garnet rim, slow diffusion through the interior of the garnet grains, and by fast diffusion along the garnet/plagioclase and the garnet/hedenbergite phase boundaries. Numerical simulation indicates that diffusion of Ca, Al and Fe²⁺ occurred about three to four, four and six to seven orders of magnitude faster along the grain boundaries than through the interior of the garnet grains. Fast‐path diffusion along grain boundaries contributed substantially to the bulk material transfer across the growing garnet rim. Despite the contribution of fast‐path diffusion, bulk diffusion through the garnet rim was too slow to allow for chemical equilibration of the phases involved in garnet rim formation even on a micrometre scale. Based on published garnet volume diffusion data the growth interval of a 20‐μm‐wide garnet rim is estimated at ～10³–10⁴ years at the inferred reaction conditions of 760 ± 50 °C at 7.6 kbar. Using the same parameterization of the growth law, 100‐μm‐ and 1‐mm‐thick garnet rims would grow within 10⁵–10⁶ and 10⁶–10⁷ years respectively.     

The role of the grain boundary at chemical and isotopic fronts in marble during contact metamorphism

Carbon and oxygen isotopic profiles around a low pressure metasomatic wollastonite reaction front in a marble of the Hida metamorphic terrain, central Japan, display typical metamorphic fluid-enhanced isotopic zonations. Isotopic profiles obtained from detailed microscale analyses perpendicular to the chemical reaction front in calcite marble show that diffusion-enhanced isotopic exchange may control these profiles. Carbon and oxygen isotopic behaviour in grain boundaries is remarkably different. Oxygen isotopic troughs (¹⁸O depleted rims) around the calcite-grain boundaries are widely observed in this contact aureole, demonstrating that diffusion of oxygen in calcite grain boundary dominates over lattice diffusion in calcite. In contrast, no difference is observed in carbon isotopic profiles obtained from grain cores and rims. There is thus no specific role of the grain boundary for diffusion of carbonic species in the metamorphic fluid during transportation. Carbon chemical species such as CO₂ and CO₃ ions in metamorphic fluid migrate mainly through lattice diffusion. The carbon and oxygen isotope profiles may be modelled by diffusion into a semi-infinite medium. Empirically lattice diffusion of oxygen isotopes is almost six times faster than that of carbon isotopes, and oxygen grain-boundary diffusion is ten times faster than oxygen lattice diffusion. Oxygen isotopic results around the wollastonite vein indicate that migration of the metamorphic fluid into calcite marble was small and was parallel to the aquifer. From the stability of wollastonite and the attainment of oxygen isotopic equilibrium, we suggest that diffusion of oxygen occurred through an aqueous fluid phase. The timescale of formation of the oxygen isotopic profile around the wollastonite vein is calculated to be about 0.76 × 10⁶ years using the experimentally determined diffusion constant. Received: 14 January 1997 / Accepted: 23 April 1998     

Volume and grain boundary diffusion of calcium in natural and hot-pressed calcite aggregates

 Calcium self-diffusion rates in natural calcite single crystals were experimentally determined at 700 to 900° C and 0.1 MPa in a stream of CO₂. Diffusion coefficients (D) were determined from ⁴²Ca concentration profiles measured with an ion microprobe. The Arrhenius parameters yield an activation energy (Q)=382±37 kJ/mol and pre-exponential factor (D₀)=0.13 m²/s, and there is no measurable anisotropy. Calcium grain boundary diffusion rates were experimentally determined in natural (Solnhofen) limestone and hot-pressed calcite aggregates at 650° to 850° C and 0.1 to 100 MPa pressure. The Solnhofen limestone was first pre-annealed for 24 h at 700° C and 100 MPa confining pressure under anhydrous conditions to produce an equilibrium microstructure for the diffusion experiments. Values for the product of the grain boundary diffusion coefficient (D′) and the effective grain boundary diffusion width (δ) were determined from ⁴²Ca concentration profiles measured with an ion microprobe. The results show that there is no measurable difference between D′δ values obtained for pre-annealed Solnhofen samples at 0.1 and 100 MPa or between hot-pressed calcite aggregates and pre-annealed Solnhofen samples. The temperature dependence for calcium grain boundary diffusion in Solnhofen samples annealed at 0.1 MPa is described by the Arrhenius parameters D^′ ₀δ=1.5×10⁻⁹ m³/s and Q=267±47 kJ/mol. Comparison of the results of this study with previously published data show that calcium is the slowest volume diffusing species in calcite. The calcium diffusivities measured in this study place constraints on several geological processes that involve diffusive mass transfer including diffusion-accommodated mechanisms in the deformation of calcite rocks. Received: 19 December 1994/Accepted: 30 June 1995     

Diffusion-controlled growth of wollastonite rims between quartz and calcite: comparison between nature and experiment

Growth rates of wollastonite reaction rims between quartz and calcite were experimentally determined at 0.1 and 1 GPa and temperatures from 850 to 1200 °C. Rim growth follows a parabolic rate law indicating that this reaction is diffusion‐controlled. From the rate constants, the D′δ‐values of the rate‐limiting species were derived, i.e. the product of grain boundary diffusion coefficient D′ and the effective grain boundary width, δ. In dry runs at 0.1 GPa, wollastonite grew exclusively on quartz surfaces. From volume considerations it is inferred that (D′_CaOδ)/(D′_SiO2δ)≥1.33, and that SiO₂ diffusion controls rim growth. D′_SiO2δ increases from about 10^?25 to 10^?23 m³ s^?1 as temperature increases from 850 to 1000 °C, yielding an apparent activation energy of 330±36 kJ mol^?1. In runs at 1 GPa, performed in a piston‐cylinder apparatus, there were always small amounts of water present. Here, wollastonite rims always overgrew calcite. Rims around calcite grains in quartz matrix are porous and their growth rates are controlled by a complex diffusion‐advection mechanism. Rim growth on matrix calcite around quartz grains is controlled by grain boundary diffusion, but it is not clear whether CaO or SiO₂ diffusion is rate‐limiting. D′δ increases from about 10^?21 to 10^?20 m³ s^?1 as temperature increases from 1100 to 1200 °C. D′_SiO2δ or D′_CaOδ in rims on calcite is c. 10 times larger than D′_SiO2δ in dry rims at the same temperature. Growth structures of the experimentally produced rims are very similar to contact‐metamorphic wollastonite rims between metachert bands and limestone in the Bufa del Diente aureole, Mexico, whereby noninfiltrated metacherts correspond to dry and brine‐infiltrated metacherts to water‐bearing experiments. However, the observed diffusivities were 4 to 5 orders of magnitude larger during contact‐metamorphism as compared to our experimental results.     

Sr/Ca and Ca/Ca fractionation during inorganic calcite formation: II. Ca isotopes

Ca isotope fractionation during inorganic calcite formation was experimentally studied by spontaneous precipitation at various precipitation rates (1.8 < log R < 4.4 μmol/m²/h) and temperatures (5, 25, and 40 °C) with traces of Sr using the CO₂ diffusion technique.Results show that in analogy to Sr/Ca [see Tang J., Köhler S. J. and Dietzel M. (2008) Sr²⁺/Ca²⁺ and ⁴⁴Ca/⁴⁰Ca fractionation during inorganic calcite formation: I. Sr incorporation. Geochim. Cosmochim. Acta] the ⁴⁴Ca/⁴⁰Ca fractionation during calcite formation can be followed by the Surface Entrapment Model (SEMO). According to the SEMO calculations at isotopic equilibrium no fractionation occurs (i.e., the fractionation coefficient α_calcite-aq = (⁴⁴Ca/⁴⁰Ca)_s/(⁴⁴Ca/⁴⁰Ca)_aq = 1 and Δ^44/40Ca_calcite-aq = 0‰), whereas at disequilibrium ⁴⁴Ca is fractionated in a primary surface layer (i.e., the surface entrapment factor of ⁴⁴Ca, F_44Ca < 1). As a crystal grows at disequilibrium, the surface-depleted ⁴⁴Ca is entrapped into the newly formed crystal lattice. ⁴⁴Ca depletion in calcite can be counteracted by ion diffusion within the surface region. Our experimental results show elevated ⁴⁴Ca fractionation in calcite grown at high precipitation rates due to limited time for Ca isotope re-equilibration by ion diffusion. Elevated temperature results in an increase of ⁴⁴Ca ion diffusion and less ⁴⁴Ca fractionation in the surface region. Thus, it is predicted from the SEMO that an increase in temperature results in less ⁴⁴Ca fractionation and the impact of precipitation rate on ⁴⁴Ca fractionation is reduced.A highly significant positive linear relationship between absolute ⁴⁴Ca/⁴⁰Ca fractionation and the apparent Sr distribution coefficient during calcite formation according to the equation

Δ^44/40Ca_calcite-aq=(−1.90±0.26)·logD_Sr−2.83±0.28     

Intragrain oxygen isotope zoning in titanite by SIMS: Cooling rates and fluid infiltration along the Carthage‐Colton Mylonite Zone,Adirondack Mountains,NY, USA

Oxygen isotopes are an attractive target for zoning studies because of the ubiquity of oxygen‐bearing minerals and the dependence of mineral ¹⁸O/¹⁶O ratios on temperature and fluid composition. In this study, subtle intragrain oxygen isotope zoning in titanite is resolved at the 10‐μm scale by secondary ion mass spectrometry. The patterns of δ¹⁸O zoning differ depending on microstructural context of individual grains and reflect multiple processes, including diffusive oxygen exchange, partial recrystallization, grain‐size reduction, and grain growth. Using the chronological framework provided by structural relations, these processes can be related to specific events during the Grenville orogeny. Titanite was sampled from two outcrops within the Carthage‐Colton Mylonite Zone (CCMZ), a long‐lived shear zone that ultimately accommodated exhumation of the Adirondack Highlands from beneath the Adirondack Lowlands during the Ottawan phase (1090–1020 Ma) of the Grenville orogeny. Titanite is hosted in the Diana metasyenite complex, which preserves three sequentially developed fabrics: an early NW‐dipping protomylonitic fabric (S₁) is crosscut by near‐vertical ultramylonitic shear zones (S₂), which are locally reoriented by a NNW‐dipping mylonitic fabric (S₃). Texturally early titanite (pre‐S₂) shows diffusion‐dominated δ¹⁸O zoning that records cooling from peak Ottawan, granulite‐facies conditions. Numerical diffusion models in the program Fast Grain Boundary yield good fits to observed δ¹⁸O profiles for cooling rates of 50 ± 20 °C Ma^?1, which are considerably faster than the 1–5 °C Ma^?1 cooling rates previously inferred for the Adirondack Highlands from regional thermochronology. High cooling rates are consistent with an episode of rapid shearing and exhumation along the CCMZ during gravitational collapse of the Ottawan orogen at c. 1050 Ma. Texturally later titanite (syn‐S₂) has higher overall δ¹⁸O and shows a transition from diffusion‐dominated to recrystallization‐dominated δ¹⁸O zoning, indicating infiltration of elevated‐δ¹⁸O fluids along S₂ shear zones and continued shearing below the blocking temperature for oxygen (~≤500 °C for grain sizes at the study site). The texturally latest titanite (post‐S₃) has growth‐dominated δ¹⁸O zoning, consistent with porphyroblastic grain growth following cessation of shearing along the Harrisville segment of the CCMZ.     

Grain growth in synthetic marbles with added mica and water

Evolution of grain size in synthetic marbles was traced from compaction of unconsolidated powder, through primary recrystallization and normal grain growth, to a size stabilized by second phases. To form the marbles, reagent grade CaCO₃ was mixed with 0, 1 and 5 volume% mica and heat-treated under pressure with added water. Densification with negligible recrystallization occurred within one hour at 500° C and 500 MPa confining pressure. Primary recrystallization occurred at 500–550° C, causing increases of grain size of factors of 2–5. Resulting samples had uniform grain size, gently curved grain boundaries, and near-equilibrium triple junctions; they were used subsequently for normal grain growth studies. Normal grain growth occurred above 550° C; at 800° C, grain size (D) increased from 7 m (D ₀) to 65 m in 24 hours. Growth rates fit the equation, D ⁿ -D ₀ ⁿ =Kt, where K is a constant and n2.6. Minor amounts of pores or mica particles inhibit normal grain growth and lead to a stabilized grain size, D _max, which depends on the size of the second phases and the inverse of their volume fraction raised to a power between 0.3 and 1. Once D _max is reached, normal growth continues only if second phases are mobile or coarsen, or if new driving forces are introduced that cause unpinning of boundaries. Normal grain growth in Solnhofen limestone was significantly slower than in pure synthetic marble, suggesting that migration is also inhibited by second phases in the limestone.     

Grain boundary and volume diffusion experiments in yttrium aluminium garnet bicrystals at 1,723 K: a miniaturized study

Yb-Y inter-diffusion along a single grain boundary of a synthetic yttrium aluminium garnet (YAG) bicrystal has been studied using analytical transmission electron microscopy (ATEM). To investigate the diffusion, a thin-film containing Yb as the diffusant was deposited perpendicular to the bicrystal grain boundary by pulsed laser deposition (PLD). Structural properties and their change with time in both the diffusant source and the grain boundary are reported. The diffusion profiles are incorporated in a numerical diffusion model, which is applied to determine the grain boundary diffusion coefficient, D _gb , at 1.723 K it is equal to 3 × 10⁻¹⁵ m²/s. We find that grain boundary diffusion is 4.85 orders of magnitude faster than volume diffusion, which was determined from the same diffusion experiment. This result is discussed in the context of special versus general grain boundaries. Finally, we successfully tested the capability of synchrotron-based nano-X-ray fluorescence analysis to map micro-chemical patterns in two dimensions with sub-micrometre resolution.     

Grain growth kinetics of dolomite,magnesite and calcite: a comparative study

The rates of grain growth of stoichiometric dolomite [CaMg(CO₃)₂] and magnesite (MgCO₃) have been measured at temperatures T of 700–800°C at a confining pressure P _c of 300 MPa, and compared with growth rates of calcite (CaCO₃). Dry, fine-grained aggregates of the three carbonates were synthesized from high purity powders by hot isostatic pressing (HIP); initial mean grain sizes of HIP-synthesized carbonates were 1.4, 1.1, and 17 μm, respectively, for CaMg(CO₃)₂, MgCO₃, and CaCO₃, with porosities of 2, 28, and 0.04% by volume. Grain sizes of all carbonates coarsened during subsequent isostatic annealing, with mean values reaching 3.9, 5.1, and 27 μm for CaMg(CO₃)₂, MgCO₃, and CaCO₃, respectively, in 1 week. Grain growth of dolomite is much slower than the growth rates of magnesite or calcite; assuming normal grain growth and n = 3 for all three carbonates, the rate constant K for dolomite (≃5 × 10⁻⁵ μm³/s) at T = 800°C is less than that for magnesite by a factor of ~30 and less than that for calcite by three orders of magnitude. Variations in carbonate grain growth may be affected by differences in cation composition and densities of pores at grain boundaries that decrease grain boundary mobility. However, rates of coarsening correlate best with the extent of solid solution; K is the largest for calcite with extensive Mg substitution for Ca, while K is the smallest for dolomite with negligible solid solution. Secondary phases may nucleate at advancing dolomite grain boundaries, with implications for deformation processes, rheology, and reaction kinetics of carbonates.     

Mechanisms of oxygen isotopic exchange and isotopic evolution of <Superscript>18</Superscript>O/<Superscript>16</Superscript>O-depleted periclase zone marbles in the Alta aureole,Utah: insights from ion microprobe analysis of calcite

Oxygen isotope ratios have been measured by ion microprobe and millimeter-scale dental drill along detailed sampling traverses across the boundary between periclase-bearing (δ¹⁸O = 11.8‰) and periclase-free (δ¹⁸O = 17.2‰) marble layers in the periclase (Per) zone of the Alta Stock aureole, Utah. These data define a steep, coherent gradient in δ¹⁸O that is displaced a short distance (~4 cm) into the periclase-free (Cal + Fo) layer. SEM and ion microprobe analyses show two isotopically and texturally distinct types of calcite at the grain scale. Clear (well polished) calcite grains are isotopically homogeneous (within analytical uncertainty; ±0.27‰, 2SD). More poorly polished (pitted), texturally retrograde ‘turbid’-looking calcite has lower and more variable δ¹⁸O values, and replaces clear calcite along fractures, cleavage traces or grain boundaries. Despite significant lowering of the δ¹⁸O values in calcite throughout both layers during prograde metamorphism, ion microprobe analyses indicate that individual clear calcite grains are now isotopically homogeneous across the entire gradient in δ¹⁸O. Diffusion calculations indicate that conservative time scales required for isotopic homogenization of calcite grains by volume diffusion, 30,000–62,000 years at 575–600°C, exceed significantly the timescale (~1,250 years) estimated for the prograde development of the δ¹⁸O gradient at the boundary between these two marble layers. The ion microprobe data and these diffusion calculations suggest instead that surface reaction mechanisms accompanying recrystallization are responsible for the observed oxygen isotope homogeneity of these calcite grains. Thus, the ion microprobe data are consistent with the formation of calcite in oxygen isotope exchange equilibrium with infiltrating fluid during prograde reaction and recrystallization. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Ca isotopes in carbonate sediment and pore fluid from ODP Site 807A: The Ca(aq)-calcite equilibrium fractionation factor and calcite recrystallization rates in Pleistocene sediments

The calcium isotopic compositions (δ⁴⁴Ca) of 30 high-purity nannofossil ooze and chalk and 7 pore fluid samples from ODP Site 807A (Ontong Java Plateau) are used in conjunction with numerical models to determine the equilibrium calcium isotope fractionation factor (α_s−f) between calcite and dissolved Ca²⁺ and the rates of post-depositional recrystallization in deep sea carbonate ooze. The value of α_s−f at equilibrium in the marine sedimentary section is 1.0000 ± 0.0001, which is significantly different from the value (0.9987 ± 0.0002) found in laboratory experiments of calcite precipitation and in the formation of biogenic calcite in the surface ocean. We hypothesize that this fractionation factor is relevant to calcite precipitation in any system at equilibrium and that this equilibrium fractionation factor has implications for the mechanisms responsible for Ca isotope fractionation during calcite precipitation. We describe a steady state model that offers a unified framework for explaining Ca isotope fractionation across the observed precipitation rate range of ∼14 orders of magnitude. The model attributes Ca isotope fractionation to the relative balance between the attachment and detachment fluxes at the calcite crystal surface. This model represents our hypothesis for the mechanism responsible for isotope fractionation during calcite precipitation. The Ca isotope data provide evidence that the bulk rate of calcite recrystallization in freshly-deposited carbonate ooze is 30-40%/Myr, and decreases with age to about 2%/Myr in 2-3 million year old sediment. The recrystallization rates determined from Ca isotopes for Pleistocene sediments are higher than those previously inferred from pore fluid Sr concentration and are consistent with rates derived for Late Pleistocene siliciclastic sediments using uranium isotopes. Combining our results for the equilibrium fractionation factor and recrystallization rates, we evaluate the effect of diagenesis on the Ca isotopic composition of marine carbonates at Site 807A. Since calcite precipitation rates in the sedimentary column are many orders of magnitude slower than laboratory experiments and the pore fluids are only slightly oversaturated with respect to calcite, the isotopic composition of diagenetic calcite is likely to reflect equilibrium precipitation. Accordingly, diagenesis produces a maximum shift in δ⁴⁴Ca of +0.15‰ for Site 807A sediments but will have a larger impact where sedimentation rates are low, seawater circulates through the sediment pile, or there are prolonged depositional hiatuses.     

An empirical estimate of the diffusion rate of oxygen in diopside

The intracrystalline diffusion rate of oxygen in diopside was constrained based on natural isotopic variations from a granulite facies marble from Cascade Slide, Adirondacks (New York, USA). The oxygen isotope compositions of the diopsides, measured as a function of grain size, are nearly constant (20.9 ± 0.3‰ vs. SMOW) over the entire measured size range (0.3–3.2 mm diameter). The δ¹⁸O values of the cores of calcite grains are 23.0‰. Temperature estimates based on the Δ¹⁸O(calcite-diopside) are 800d?C, in agreement with the highest previous thermometric estimates for these rocks. The lack of isotopic variation in the diopsides as a function of grain size requires that the oxygen intracrystalline diffusion rate in diopside from the Adirondack samples was very slow. The maximum diffusion rates (D_800d?C parallel to the c-axis) were calculated with an infinite reservoir model (IRM) and a finite reservoir model (FRM) that incorporates mineral modal abundances and initial isotopic variations. For an assumed activation energy (Q) = 100 kJ/mol, the IRM diffusion rate estimate of 1.6 times 10^-20cm²/s is two orders of magnitude faster than from the FRM; at Q=500kJ/mol, the D_800d?C estimate for both methods is c. 5.6 times 10^-20 cm²/s. The present results require that a hydrothermal fluid significantly enhances the diffusion rate of oxygen in diopside if previous data are correct. The δ¹⁸O(SMOW) and δ¹³C(PDB) values of the calcite, measured in situ with a CO₂ laser, are 22.9 ± 0.3, 0.1±0.3‰ in the grain cores, 22.1 ±0.3, 0.2 ±0.1‰ at the grain boundaries and 21.7 ±0.4, -0.6±0.1‰ abutting diopside grains. The δ¹⁸O and δ¹³δC values measured conventionally are: crystal cores, 22.96, -0.95‰; abutting diopside grains, 22.38, -0.93‰; bulk, 22.79, -0.95%. Use of the bulk δ¹⁸O(calcite) values for thermometry yields unreasonably high temperatures. The lower δ¹⁸O values at the calcite grain boundaries are not due to retrograde diffusional exchange with the diopside, they are thought to be a result of a late retrograde fluid infiltration.     

Diffusion along interphase boundaries and its effect on retrograde zoning patterns of metamorphic minerals

In this study we use two dimensional chemical patterns and numerical modeling to estimate the relative rates of chemical transport along interphase boundaries (ib) and through grain (s) interiors during retrograde Fe–Mg exchange between garnet and biotite at a garnet–biotite–quartz triple junction. We demonstrate that systematic variations in garnet–rim compositions, which are frequently observed along garnet–quartz interfaces, and deviations from concentric retrograde zoning patterns start to develop when chemical transport along the interphase boundaries becomes slow during cooling. The capacities for chemical transport along an interphase boundary depend on the product D ^ib K ^ib/s a, where D ^ib is the diffusion coefficient of the exchangeable components within the interphase boundary medium, K ^ib/s is the equilibrium partitioning coefficient between the cation exchange partners and the interphase boundary medium and a is the interphase boundary width. The model is applied to analyze the retrograde zoning patterns in garnets from the Mozambique belt (SE-Kenya), which cooled from 820°C at a rate of ca. 2°C/my. It is found that non-equilibrated compositions in garnet along garnet/quartz interphase boundaries started to develop below 700°C due to insufficient rates of chemical transport along these boundaries. The transport capacities of garnet/quartz interphase boundaries was estimated to have been between about 1E-23 cm³/s (575°C) and 1E-20 cm³/s (700°C) from modeling the observed X _Fe pattern in garnet close to a garnet–quartz–biotite triple junction and relying on published data on the diffusivity of Fe²⁺ in garnet. Similar transport capacities are obtained; when the interphase boundary is assumed to be filled with a material that has the transport properties and chemical composition of a free melt in equilibrium with garnet, biotite and quartz at the respective conditions. In contrast, if the transport properties of the interphase boundary medium are related to the diffusivities and solubility of Fe²⁺ and FeOH⁺ within a free aqueous solution, chemical transport along the interphase boundaries would be much more efficient, and exchange equilibrium would have been maintained during the entire cooling history of the rocks. The observation of systematic deviations from local equilibrium along the garnet–quartz interphase boundaries leads us to exclude the presence of an aqueous fluid along the interphase boundary at any time during cooling.     

Timescales and mechanisms of fluid infiltration in a marble: an ion microprobe study

Using a recently developed ion microprobe technique, a detailed oxygen isotope map of calcite grains in a coarse-grained marble has been constructed, supported by trace element (Mn, Sr, Fe) analysis and cathodoluminescence (CL) imaging, in order to constrain scales of oxygen isotope equilibrium, timescales and mechanisms of metamorphic fluid infiltration, and fluid sources and pathways. Results are compared with a previous study of this sample (Wada 1988) carried out using a cryo-microtome technique and conventional oxygen isotope analysis. The marble, from the high temperature/low pressure Hida metamorphic belt in north-central Japan, underwent granulite facies followed by amphibolite facies metamorphic events, the latter associated with regional granite intrusion. The CL imaging indicates two types of calcite, a yellow luminescing (YLC) and a purple luminescing (PLC) variety. The YLC, which occupies grain boundaries, fractures, replacement patches, and most of the abundant deformation twin lamellae, post-dates the dominant PLC calcite and maps out fluid pathways. Systematic relationships were established between oxygen isotope and trace element composition, calcite type and texture, based on 74 ¹⁸O/¹⁶O and 17 trace element analyses with 20–30 μ m spatial resolution. The YLC is enriched in Mn and Fe, and depleted in ¹⁸O and Sr compared to PLC, and is much more ¹⁸O depleted than is indicated from conventional analyses. Results are interpreted to indicate infiltration of ¹⁸O-depleted (metamorphic or magmatic) fluid (initial δ¹⁸O = 9‰–10.5‰) along grain boundaries, fractures and deformation twin lamellae, depleting calcite grains in Sr and enriching them in Mn and Fe. The sample is characterised by gross isotopic and elemental disequilibrium, with important implications for the application of chromatographic theory to constrain fluid fluxes in metacarbonate rocks. Areas of PLC unaffected by “short-circuiting” fluid pathways contain oxygen diffusion profiles of ∼10‰/∼200 μm in grain boundary regions or adjacent to fractures/patches. When correction is made for estimated grain boundary/fracture and profile orientation in 3D, profiles are indistinguishable within error. Modelling of these profiles gives consistent estimates of Dt (where D is the diffusion coefficient and t is time) of ∼0.8 × 10⁻⁸ m², from which, using experimental data for oxygen diffusion in calcite, timescales of fluid transport along grain boundaries at amphibolite facies temperatures of ∼10³ to ∼10⁴ years are obtained. These short timescales, which are much shorter than plausible durations of metamorphism, imply that rock permeabilities may be transiently much higher during fluid flow than those calculated from time integrated fluid fluxes or predicted from laboratory measurements. The preservation of ¹⁸O/¹⁶O profiles requires either rapid cooling rates (∼100–600 °C/million years), or, more plausibly, loss of grain boundary fluid such that a dry cooling history followed the transient passage of fluid. The δ¹⁸O/trace element correlations are also consistent with volume diffusion-controlled transport in the PLC. Fluid transport and element exchange occurred by two inter-related mechanisms on short timescales and on different lengthscales – long-distance flow along cracks, grain boundaries and twin lamellae coupled to ∼200 μm-scale volume diffusion of oxygen. Received: 8 December 1997 / Accepted: 18 May 1998     

Modelling grain‐recycling zoning during metamorphism

Chemical zoning, recorded by grain growth during metamorphism, is a key source of information about P–T–t paths. Interpretation of these data must be carried out using appropriate models and recognizing their inherent assumptions. To assist with defining how zoned minerals form, a set of geometric criteria for three types of chemical zoning developed in minerals (diffusion, growth and grain recycling) is outlined. Re‐equilibration of minerals by lattice diffusion causes zoning if the re‐equilibration is incomplete. Growth of porphyroblasts is commonly considered in pelites, but in metagranitoids, large monophase domains undergo coarsening by recycling of material from one grain to another as grain boundaries migrate driven by surface energy. This type of grain size increase is termed here ‘grain recycling’. Zoning developed during grain recycling due to equilibration of the recycled material with grain‐boundary chemistry is termed ‘grain‐recycling zoning’. Furthermore, short lattice diffusion lengths relative to grain sizes cause metamorphic fractionation because material in the grain cores is not in communication thermodynamically with the rest of the rock. A new model is derived for this sort of grain size increase coupled with metamorphic reactions using Theriak–Domino. An example is given of plagioclase undergoing an increase in anorthite content as epidote breaks down during amphibolite facies metamorphism of a metagranitoid. Agreement between naturally occurring zoning profiles and those derived from modelled P–T–t paths shows that this model can be used to extract metamorphic conditions from rocks which are not accessible using conventional thermobarometry.     

Compositional changes of minerals associated with dynamic recrystallizatin

The rate of compositional and isotopic exchange between minerals may be enhanced significantly if the rock is deformed simultaneously. The enhanced exchange rate may result from a reduction in grain size (shorter distance for volume diffusion), dissolution and growth of grains by diffusion creep (pressure solution), or the movement of high-angle grain boundaries through strained grains during recrystallization in the dislocation creep regime. The migration of high-angle grain boundaries provides high diffusivity paths for the rapid exchange of components during recrystallization. The operation of the latter process has been demonstrated by deforming aggregates consisting of two plagioclases (An₁ and An₇₉) at 900°C, 1 GPa confining pressure, and a strain rate of ∼2x10^-6s^-1. The polygonal, recrystallized grains were analyzed using an analytical transmission electron microscope and have a variable but often intermediate composition. At the conditions of these experiments, the volume interdiffusion rate of NaSi/CaAl is too slow to produce any observable chemical change, and microstructural-chemical relations indicate that the contribution from diffusion creep was insignificant except for initially fine-grained (2–10 μm) aggregates. These results indicate that strain-induced recrystallization can be an effective mechanism for enhancing the kinetics of metamorphic reactions and for resetting the isotope systematics of minerals such as feldspars, pyroxenes, and amphiboles.     

Sr isotopes and pore fluid chemistry in carbonate sediment of the Ontong Java Plateau: Calcite recrystallization rates and evidence for a rapid rise in seawater Mg over the last 10 million years

The ⁸⁷Sr/⁸⁶Sr ratios and Sr concentrations in sediment and pore fluids are used to evaluate the rates of calcite recrystallization at ODP Site 807A on the Ontong Java Plateau, an 800-meter thick section of carbonate ooze and chalk. A numerical model is used to evaluate the pore fluid chemistry and Sr isotopes in an accumulating section. The deduced calcite recrystallization rate is 2% per million years (%/Myr) near the top of the section and decreases systematically in older parts of the section such that the rate is close to 0.1/age (in years). The deduced recrystallization rates have important implications for the interpretation of Ca and Mg concentration profiles in the pore fluids. The effect of calcite recrystallization on pore fluid chemistry is described by the reaction length, L, which varies by element, and depends on the concentration in pore fluid and solid. When L is small compared to the thickness of the sedimentary section, the pore fluid concentration is controlled by equilibrium or steady-state exchange with the solid phase, except within a distance L of the sediment-water interface. When L is large relative to the thickness of sediment, the pore fluid concentration is mostly controlled by the boundary conditions and diffusion. The values of L for Ca, Sr, and Mg are of order 15, 150, and 1500 meters, respectively. L_Sr is derived from isotopic data and modeling, and allows us to infer the values of L_Ca and L_Mg. The small value for L_Ca indicates that pore fluid Ca concentrations, which gradually increase down section, must be equilibrium values that are maintained by solution-precipitation exchange with calcite and do not reflect Ca sources within or below the sediment column. The pore fluid Ca measurements and measured alkalinity allow us to calculate the in situ pH in the pore fluids, which decreases from 7.6 near the sediment-water interface to 7.1 ± 0.1 at 400-800 mbsf. While the calculated pH values are in agreement with some of the values measured during ODP Leg 130, most of the measurements are artifacts. The large value for L_Mg indicates that the pore fluid Mg concentrations at 807A are not controlled by calcite-fluid equilibrium but instead are determined by the changing Mg concentration of seawater during deposition, modified by aqueous diffusion in the pore fluids. We use the pore fluid Mg concentration profile at Site 807A to retrieve a global record for seawater Mg over the past 35 Myr, which shows that seawater Mg has increased rapidly over the past 10 Myr, rather than gradually over the past 60 Myr. This observation suggests that the Cenozoic rise in seawater Mg is controlled by continental weathering inputs rather than by exchange with oceanic crust. The relationship determined between reaction rate and age in silicates and carbonates is strikingly similar, which suggests that reaction affinity is not the primary determinant of silicate dissolution rates in nature.     

The legacy of crystal-plastic deformation in olivine: high-diffusivity pathways during serpentinization

Crystal-plastic olivine deformation to produce subgrain boundaries composed of edge dislocations is an inevitable consequence of asthenospheric mantle flow. Although crystal-plastic deformation and serpentinization are spatio-temporally decoupled, we identified compositional readjustments expressed on the micrometric level as a striped Fe-enriched ( [`(X)]<sub>\textFe</sub> \bar{X}_{\text{Fe}}  = 0.24 ± 0.02 (zones); 0.12 ± 0.02 (bulk)) or Fe-depleted ( [`(X)]_\textFe \bar{X}_{\text{Fe}}  = 0.10 ± 0.01 (zones); 0.13 ± 0.01 (bulk)) zoning in partly serpentinized olivine grains from two upper mantle sections in Norway. Focused ion beam sample preparation combined with transmission electron microscopy (TEM) and aberration-corrected scanning TEM, enabling atomic-level resolved electron energy-loss spectroscopic line profiling, reveals that every zone is immediately associated with a subgrain boundary. We infer that the zonings are a result of the environmental Fe²⁺Mg₋₁ exchange potential during antigorite serpentinization of olivine and the drive toward element exchange equilibrium. This is facilitated by enhanced solid-state diffusion along subgrain boundaries in a system, which otherwise re-equilibrates via dissolution-reprecipitation. Fe enrichment or depletion is controlled by the silica activity imposed on the system by the local olivine/orthopyroxene mass ratio, temperature and the effect of magnetite stability. The Fe-Mg exchange coefficients K_\textD^{\textAtg/\textOl} K_{\text{D}}^{{{\text{Atg}}/{\text{Ol}}}} between both types of zoning and antigorite display coalescence toward exchange equilibrium. With both types of zoning, Mn is enriched and Ni depleted compared with the unaffected bulk composition. Nanometer-sized, heterogeneously distributed antigorite precipitates along olivine subgrain boundaries suggest that water was able to ingress along them. Crystallographic orientation relationships gained via electron backscatter diffraction between olivine grain domains and different serpentine vein generations support the hypothesis that serpentinization was initiated along olivine subgrain boundaries.