Growth of multilayered polycrystalline reaction rims in the MgO–SiO<Subscript>2</Subscript> system,part II: modelling

Part I of this contribution (Gardés et al. in Contrib Mineral Petrol, 2010) reported time- and temperature-dependent experimental growth of polycrystalline forsterite-enstatite double layers between single crystals of periclase and quartz, and enstatite single layers between forsterite and quartz. Both double and single layers displayed growth rates decreasing with time and pronounced grain coarsening. Here, a model is presented for the growth of the layers that couples grain boundary diffusion and grain coarsening to interpret the drop of the growth rates. It results that the growth of the layers is such that (Δx)² ∝ t ^1−1/n , where Δx is the layer thickness and n the grain coarsening exponent, as experimentally observed. It is shown that component transport occurs mainly by grain boundary diffusion and that the contribution of volume diffusion is negligible. Assuming a value of 1 nm for the effective grain boundary width, the following Arrhenius laws for MgO grain boundary diffusion are derived: log D _gb,0^Fo (m²/s) = −2.71 ± 1.03 and E _gb^Fo = 329 ± 30 kJ/mol in forsterite and log D _gb,0^En (m²/s) = 0.13 ± 1.31 and E _gb^En = 417 ± 38 kJ/mol in enstatite. The different activation energies are responsible for the changes in the enstatite/forsterite thickness ratio with varying temperature. We show that significant biases are introduced if grain boundary diffusion-controlled rim growth is modelled assuming constant bulk diffusivities so that differences in activation energies of more than 100 kJ/mol may arise. It is thus important to consider grain coarsening when modelling layered reaction zones because they are usually polycrystalline and controlled by grain boundary transport.     

Grain boundary diffusion of Si, Mg, and O in enstatite reaction rims: a SIMS study using isotopically doped reactants

Diffusion-controlled growth rates of polycrystalline enstatite reaction rims between forsterite and quartz were determined at 1,000 °C and 1 GPa in presence of traces of water. Iron-free, pure synthetic forsterite with normal oxygen and silicon isotopic compositions and quartz extremely enriched in ¹⁸O and ²⁹Si were used as reactants. The relative mobility of ¹⁸O and ²⁹Si in reactants and rims were determined by SIMS step scanning. The morphology of the rim shows that enstatite grows by a direct replacement of forsterite. Rim growth is modelled within a mass-conserving reference frame that implies advancement of reaction fronts from the initial forsterite-quartz interface in both directions. The isotopic compositions at the two reaction interfaces are controlled by the partial reactions Mg₂SiO₄=0.5 Mg₂Si₂O₆+MgO at the forsterite-enstatite, and MgO+SiO₂=0.5 Mg₂Si₂O₆ at the enstatite-quartz interface, implying that grain boundary diffusion of MgO is rate-controlling. Isotopic profiles show no silicon exchange across the propagating reaction interfaces. This propagation, controlled by MgO diffusion, is faster than the homogenisation of Si by self-diffusion behind the advancing fronts. From this, and using % MathType!MTEF!2!1!+- % feaaeaart1ev0aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn % hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr % 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9 % vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x % fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiramaaDa % aaleaacaWGtbGaamyAaiaacYcacaWGfbGaamOBaaqaaiaadAfacaWG % VbGaamiBaaaaaaa!3DD2! D_Si,En^VolD_{Si,En}^{Vol} at dry conditions from the literature, results a % MathType!MTEF!2!1!+- % feaaeaart1ev0aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn % hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr % 4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq-Jc9 % vqaqpepm0xbba9pwe9Q8fs0-yqaqpepae9pg0FirpepeKkFr0xfr-x % fr-xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGabmirayaafa % Waa0baaSqaaiaadofacaWGPbGaaiilaiaadweacaWGUbaabaaaaOGa % eqiTdqgaaa!3CCD! D￠_Si,En dD'_{Si,En}^{} \delta value of 3᎒^-24 m³ s^-1 at 1,000 °C. The isotopic profiles for oxygen are more complex. They are interpreted as an interplay between the propagation of the interfaces, the homogenisation of the isotope concentrations by grain boundary self-diffusion of O within the rim, and the isotope exchange across the enstatite-quartz interface, which was open to ¹⁸O influx from quartz. Because of overlapping diffusion processes, boundary conditions are unstable and D´_Ox,En' cannot be quantified. Using measured rim growth rates, the grain boundary diffusivity D´_MgO' of MgO in iron-free enstatite is 8᎒^-22 m³ s^-1 at 1,000 °C and 1 GPa. Experiments with San Carlos olivine (fo₉₂) as reactant reveal lower rates by a factor of about 4. Our results show that isotope tracers in rim growth experiments allow identification of the actual interface reactions, recognition of the rate-controlling component and further calculation of D´' values for specific components.     

Diffusion-controlled growth of albite and pyroxene reaction rims

The growth rates of albite and pyroxene (enstatite + diopside + spinel) reaction rims were measured at 1000°C and ˜700 MPa and found to be parabolic indicating diffusion-controlled growth. The parabolic rate constants for the pyroxene (+ spinel) rims in samples with 0.5 wt% H₂O added or initially vacuum dried at 25°C and 250°C are 1.68 ± 0.09, 0.54 ± 0.05 and 0.25 ± 0.06 μm²/h, respectively. The values for albite rim growth in samples initially dried at 60°C and with 0.1 wt% H₂O added are 0.25 ± 0.04 and 0.33 ± 0.03 μm²/h, respectively. The latter values were used to derive the product of the grain boundary diffusion coefficient D′_A, where A = SiO₂, NaAlO₂, or NaAlSi₋₁, and the grain boundary thickness δ in albite. The calculated D′_SIO2δ in the albite aggregate for the situations of two different water contents are about 9.9 × 10⁻²³ and 1.4 × 10⁻²² m³ s⁻¹, respectively. Both the rate constants and the calculated D′_Aδ demonstrate that the effect of water content on the grain boundary diffusion rate in monomineralic albite and polymineralic pyroxene (+ spinel) aggregates is small, consistent with recent studies of monomineralic enstatite and forsterite rims. Received: 1 July 1995 / Accepted: 1 August 1996     

Rates of grain boundary diffusion through enstatite and forsterite reaction rims

 The growth rates of enstatite rims produced by reaction of Fo₉₂ and SiO₂ were determined at 250–1500 MPa and 900–1100°C for a wide range of water contents. Growth rates were also determined for forsterite rims between MgO and Mg₂Si₂O₆ and between MgO and SiO₂. Rim growth rates are parabolic indicating diffusion-controlled growth of the polycrystalline rims which are composed of ˜ 2 μm diameter grains. Rim growth rates were used to calculate the product of the grain boundary diffusion coefficient (D'_A) times the effective grain boundary thickness (δ) assuming in turn that MgO, SiO₂, and Mg₂Si₋₁ are the diffusing components (coupled diffusion of a cation and oxygen or interdiffusion of Mg and Si). The values for D'_MgOδ, D', and D' for enstatite at 1000°C and 700 MPa confining pressure with about 0.1 wt %  water are about five times larger than the corresponding D'_Aδ values for samples initially vacuum dried at 250°C. Most of the increase in D'_Aδ occurs with the first 0.1 wt %  water. The activation energy for diffusion through the enstatite rims (1100–950°C) is 162 ± 30 kJ/mole. The diffusion rate through enstatite rims is essentially unchanged for confining pressures from 210–1400 MPa, but the nucleation rate is greatly reduced at low confining pressure (for  ≤ 1.0 wt % water present) and limits the conditions at which rim growth can be measured. The corresponding values for D'_Aδ through forsterite rims are essentially identical for the two forsterite-producing reactions when 0.1 wt % water is added and similar to the D'_Aδ values for enstatite at the same conditions. The D'_Aδ values for forsterite are ˜ 28 times larger for samples starting with 0.1 wt %  water compared to samples that were first vacuum dried. Thus water enhances these grain boundary diffusion rates by a factor of 5–30 depending on the mineralogy, but the total range in D'_Aδ is only slightly more than an order of magnitude for as wide a range of water contents as expected for most crustal conditions. Received: 1 July 1995 / Accepted: 1 August 1996     

Zr diffusion in titanite

Chemical diffusion of Zr under anhydrous, pO₂-buffered conditions has been measured in natural titanite. The source of diffusant was either zircon powder or a ZrO₂–Al₂O₃–titanite mixture. Experiments were run in sealed silica glass capsules with solid buffers (to buffer at NNO or QFM). Rutherford Backscattering Spectrometry (RBS) was used to measure diffusion profiles. The following Arrhenius parameters were obtained for Zr diffusion parallel to c over the temperature range 753–1,100°C under NNO-buffered conditions: D _Zr = 5.33 × 10⁻⁷ exp(−325 ± 30 kJ mol⁻¹/RT) m² s⁻¹ Diffusivities are similar for experiments buffered at QFM. These data suggest that titanite should be moderately retentive of Zr chemical signatures, with diffusivities slower than those for O and Pb in titanite, but faster than those for Sr and the REE. When applied in evaluation of the relative robustness of the recently developed Zr-in-titanite geothermometer (Hayden and Watson, Abstract, 16th V.M. Goldschmidt Conference 2006), these findings suggest that Zr concentrations in titanite will be less likely to be affected by later thermal disturbance than the geothermometer based on Zr concentrations in rutile (Zack et al. in Contrib Mineral Petrol 148:471–488, 2004; Watson et al. in Contrib Mineral. Petrol, 2006), but much less resistant to diffusional alteration subsequent to crystallization than the Ti-in-Zircon geothermometer (Watson and Harrison in Science 308:841–844, 2005).     

Paragonite stability at 700°C in the presence of H<Subscript>2</Subscript>O–NaCl fluids: constraints on H<Subscript>2</Subscript>O activity and implications for high pressure metamorphism

We carried out reversed piston-cylinder experiments on the equilibrium paragonite = jadeite + kyanite + H₂O at 700°C, 1.5–2.5 GPa, in the presence of H₂O-NaCl fluids. Synthetic paragonite and jadeite and natural kyanite were used as starting materials. The experiments were performed on four different nominal starting compositions: X(H₂O)=1.0, 0.90, 0.75 and 0.62. Reaction direction and extent were determined from the weight change in H₂O in the capsule, as well as by optical and scanning electron microscopy (SEM). At X(H₂O)=1.0, the equilibrium lies between 2.25 and 2.30 GPa, in good agreement with the 2.30–2.45 GPa reversal of Holland (Contrib Miner Petrol 68:293–301, 1979). Lowering X(H₂O) decreases the pressure of paragonite breakdown to 2.10–2.20 GPa at X(H₂O)=0.90 and 1.85–1.90 GPa at X(H₂O)=0.75. The experiments at X(H₂O) = 0.62 yielded the assemblage albite + corundum at 1.60 GPa, and jadeite + kyanite at 1.70 GPa. This constrains the position of the isothermal paragonite–jadeite–kyanite–albite–corundum–H₂O invariant point in the system Na₂O–Al₂O₃–SiO₂–H₂O to be at 1.6–1.7 GPa and X(H₂O)~0.65±0.05. The data indicate that H₂O activity, a(H₂O), is 0.75–0.86, 0.55–0.58, and <0.42 at X(H₂O)=0.90, 0.75, and 0.62, respectively. These values approach X(H₂O)², and agree well with the a(H₂O) model of Aranovich and Newton (Contrib Miner Petrol 125:200–212, 1996). Our results demonstrate that the presence or absence of paragonite can be used to place limits on a(H₂O) in high-pressure metamorphic environments. For example, nearly pure jadeite and kyanite from a metapelite from the Sesia Lanzo Zone formed during the Eo-Alpine metamorphic event at 1.7–2.0 GPa, 550–650°C. The absence of paragonite requires a fluid with low a(H₂O) of 0.3–0.6, which could be due to the presence of saline brines.     

Constraints on phase diagram topology for the system CaO−MgO−SiO2−CO2−H2O

Published phase diagrams for the siliceous carbonate system CaO–MgO–SiO₂–CO₂–H₂O are contradictory because of different estimates of the relative stability of magnesite. Experimental data on magnesite are too ambiguous to determine the validity of these estimates. Therefore, field evidence is used to select the correct phase diagram topology for siliceous carbonate and carbonate ultramafic rocks at pressures of about 2–5 kbar. The primary selection criterion is provided by the existence of the stable assemblage talc+dolomite+forsterite+tremolite+antigorite, which occurs in the Bergell contact aureole and Swiss Central Alps. Field evidence also is used to argue that the reaction magnesite+quartz=enstatite must occur at lower temperature than the reaction dolomite+quartz=diopside. T-X _{CO 2} and P _{CO 2}-T phase diagrams consistent with these observations are calculated from experimental and thermo-dynamic data. For antigorite ophicarbonate rocks, remarkable agreement is obtained between the spatial distribution of low variance mineral assemblages and the calculated diagrams.     

Reaction rim growth in the system MgO-Al2O3-SiO2 under uniaxial stress

We synthesize reaction rims between thermodynamically incompatible phases in the system MgO-Al₂O₃-SiO₂ applying uniaxial load using a creep apparatus. Synthesis experiments are done in the MgO-SiO₂ and in the MgO-Al₂O₃ subsystems at temperatures ranging from 1150 to 1350 °C imposing vertical stresses of 1.2 to 29 MPa at ambient pressure and under a constant flow of dry argon. Single crystals of synthetic and natural quartz and forsterite, synthetic periclase and synthetic corundum polycrystals are used as starting materials. We produce enstatite rims at forsterite-quartz contacts, enstatite-forsterite double rims at periclase-quartz contacts and spinel rims at periclase-corundum contacts. We find that rim growth under the “dry” conditions of our experiments is sluggish compared to what has been found previously in nominally “dry” piston cylinder experiments. We further observe that the nature of starting material, synthetic or natural, has a major influence on rim growth rates, where natural samples are more reactive than synthetic ones. At a given temperature the effect of stress variation is larger than what is anticipated from the modification of the thermodynamic driving force for reaction due to the storage of elastic strain energy in the reactant phases. We speculate that this may be due to modification of the physical properties of the polycrystals that constitute the reaction rims or by deformation under the imposed load. In our experiments rim growth is very sluggish at forsterite-quartz interfaces. Rim growth is more rapid at periclase-quartz contacts. The spinel rims that are produced at periclase-corundum interfaces show parabolic growth indicating that reaction rim growth is essentially diffusion controlled. From the analysis of time series done in the MgO-Al₂O₃ subsystem we derive effective diffusivities for the Al₂O₃ and the MgO components in a spinel polycrystal as ${\rm D}_{MgO} = 1.4 \pm 0.2 \cdot 10^{-15}$  m²/s and ${\rm D}_{Al_2O_3} = 3.7 \pm 0.6 \cdot 10^{-16}$  m²/s for T?=?1350 °C and a vertical stress of 2.9 MPa.     

Diffusive transport of water in porous feldspars from granitic saprolites: In situ experiments using FTIR spectroscopy

A novel experimental cell was developed for in situ measurements of transport phenomena in porous media using Fourier-Transform Infrared (FTIR) Spectroscopy. The technique was employed at ambient pressure in the temperatures range of 11–44 °C to study the H₂O → D₂O exchange between water-saturated weathered feldspars (bulk porosity of 5–19 vol% for feldspar) from granitic saprolites and a surrounding aqueous liquid. Such measurements are an important step for understanding internal weathering reactions of feldspars in soils and aquifers. Effective diffusion coefficients D_eff for water in water-saturated porous feldspars were determined assuming one-dimensional diffusion in a quasi-homogeneous medium. The values of D_eff vary from 7.2 × 10⁻¹⁰ to 1.9 × 10⁻¹¹ m²/s and are 1–2 orders of magnitude lower than the diffusion coefficients (D) of protons and molecular H₂O in liquid water. The activation energy for the H₂O → D₂O exchange process in porous feldspars ranges from 7.8 to 18.8 kJ/mol.The results imply that the effective diffusivity of water is mainly controlled by physical properties of the feldspars like porosity, pore connectivity, pore geometry and distribution. Perthitic feldspars with homogeneous pore distribution in the albitic lamellas have diffusional tortuosity factors X = D/D_eff between 3 and 10 while alkali feldspars with inhomogeneously distributed and disconnected pores have much higher X values up to 129. Diffusion anisotropy has been verified for a vein perthite with diffusion perpendicular to the lamellas being faster by 0.3–0.5 log units than within the lamellas. It has to be emphasized that the study is based only on few selected feldspars, including perthitic feldspar, and additional work on samples with different weathering stages is needed to test the importance of the different parameters controlling diffusive transport in the pore system.     

A preliminary stable isotope study on Mogok Ruby, Myanmar

The primary occurrence of ruby in the Mogok area, northern Myanmar is exclusively found in marble along with spinel–forsterite-bearing marble and phlogopite–graphite marble. These marble units are enclosed within banded biotite–garnet–sillimanite–oligoclase gneisses. Samples of these marbles collected for C–O stable isotope analysis show two trends of δ¹³C–δ¹⁸O variation resulting most likely from fluid–rock interactions. Ruby-bearing marble and phlogopite–graphite marble follow a trend with coupled C–O depletion, whereas spinel–forsterite-bearing marble follows a δ¹⁸O depletion trend with relatively constant δ¹³C values. Ruby formation might have resulted from CO₂-rich fluid–rock interaction, while spinel–forsterite-bearing marble was genetically related to CO₂-poor fluid–rock interaction. Both fluids may have arisen from external sources. Based on graphite Raman spectral thermometry, the estimated temperature for phlogopite–graphite marble, and probably ruby-bearing marble, was lower than 607 °C, and for spinel–forsterite-bearing marble, lower than 710 °C. Contrasting C/O diffusion between graphite/ruby/spinel/forsterite and calcite, local variations of isotopic compositions of newly formed minerals as a result of non-pervasive fluid infiltration, and open-system isotopic disturbance during cooling may have affected C-/O-isotopic fractionations between minerals. The estimated high formation temperatures for ruby and spinel/forsterite imply that the parental fluids may have been related to nearby igneous intrusions and/or metamorphic processes. Whether these two types of fluid were genetically related is unclear based on the present data.     

Experimental study of incongruent evaporation kinetics of enstatite in vacuum and in hydrogen gas

Variations in bulk Mg/Si ratios in the various groups of chondritic meteorites indicate that Mg/Si fractionation occurred in the primitive solar nebula. Enstatite (MgSiO₃) evaporates incongruently forming forsterite (Mg₂SiO₄) as an evaporation residue; therefore, evaporation of enstatite produces Mg/Si variations in solid (Mg-rich) and gas (Si-rich) and must be considered as a probable process responsible for Mg/Si fractionation recorded in chondrites. To understand the evaporation kinetics of enstatite, incongruent evaporation experiments on enstatite single crystals have been carried out in vacuum and in hydrogen gas at temperatures of 1300 to 1500°C. A polycrystalline forsterite layer is formed on the surface of enstatite by preferential evaporation of the SiO₂ component, both in vacuum and in hydrogen gas. The thickness of the forsterite layer in vacuum increases with time in the early stage of evaporation and later the thickness of the forsterite layer remains constant (several microns). This is due to the change in the rate limiting process from surface reaction plus nucleation and growth to diffusion in the surface forsterite layer. The activation energy of the diffusion-controlled evaporation rate constant of enstatite is 457 (±58) kJ/mol. A thinner forsterite layer is formed on the surface of enstatite in hydrogen gas than in vacuum. Evaporation of enstatite in hydrogen gas is also considered to be controlled by diffusion of ions through the forsterite layer. The thin forsterite layer formed in hydrogen gas is ascribed to the enhanced evaporation rate of forsterite in the presence of hydrogen gas.The results are applied to incongruent evaporation under the solar nebular conditions. The steady thickness of the forsterite of nebular pressure-temperature conditions is estimated to be submicron because of the enhanced evaporation rate of forsterite under hydrogen-rich nebular conditions if evaporated gases are taken away immediately and no back reaction occurs (an open system). Because enstatite grains in the solar nebula would be comparable to the estimated steady thickness of forsterite, evaporation of such enstatite grains under kinetic conditions could play an important role in producing variations in Mg/Si ratios between solid and gas in the solar nebula.     

Zr-in-rutile thermometry in HP/UHP eclogites from Western China

Four Zr-in-rutile thermometry calibrations are applied to eclogites from Western China. Here, we show that if rutile grows in equilibrium with Qtz and Zrn, and is isolated inside garnet, it preserves its Zr composition and does not undergo compositional change due to cation exchange with the host garnet. It thus preserves the composition for the P–T conditions of its formation and the growth zoning of the host garnet. For the HP/UHP metamorphic temperature, the Tomkins et al. (J Metamorph Geol 25:703–713, 2007) calibration yields temperatures that agree well with previous studies, whereas the other three calibrations (Zack et al. in Contrib Mineral Petrol 148:471–488, 2004; Watson et al. in Contrib Mineral Petrol 151:413–433, 2006; Ferry and Watson in Contrib Mineral Petrol in 154:429–437, 2007), which do not include a pressure correction, give systematically lower temperatures. Zr contents of rutile inclusions within garnet show systematic decrease from garnet core to rim. The rutile inclusions in garnet rims contain the lowest Zr content, similar to that in the matrix. Analyses confirm that the pressure plays a significant role in modifying the primary temperature dependence of the Zr content of rutile. Rutiles trapped in garnets are unable to re-equilibrate easily during retrogression, but those in the matrix can do so, providing retrograde P–T path information.     

Grain boundary diffusion in enstatite

We have investigated grain boundary diffusion rates in enstatite by heating single crystals of quartz packed in powdered San Carlos olivine (Mg_0.90Fe_0.10)₂SiO₄ at controlled oxygen fugacities in the range 10^?5.7 to 10^?8.7?atm and temperatures from 1350° to 1450?°C for times from 5 to 100?h at 1?atm total pressure. Following the experiments, the thickness of the coherent polycrystalline reaction rim of pyroxene that had formed between the quartz and olivine was measured using backscatter scanning imaging in the electron microprobe. Quantitative microprobe analysis indicated that the composition of this reaction phase is (Mg_0.92Fe_0.08)₂Si₂O₆. The rate of growth of the pyroxene increases with increasing temperature, is independent of the oxygen fugacity, and is consistent with a parabolic rate law, indicating that the growth rate is controlled by ionic diffusion through the pyroxene rim. Microstructural observations and platinum marker experiments suggest that the reaction phase is formed at the olivine-pyroxene interface, and is therefore controlled by the diffusion of silicon and oxygen. The parabolic rate constants determined from the experiments were analyzed in terms of the oxide activity gradient across the rim to yield mean effective diffusivities for the rate-limiting ionic species, assuming bulk transport through the pyroxene layer. These effective diffusivities are faster than the lattice diffusivities for the slowest species (silicon) calculated from creep experiments, but slower than measured lattice diffusivities for oxygen in enstatite. Thus, silicon grain boundary diffusion is most likely to be the rate-limiting process in the growth of the pyroxene rims. Also, as oxygen transport through the pyroxene rims must be faster than silicon transport, diffusion of oxygen along the grain boundaries must be faster than through the lattice. The grain boundary diffusivity for silicon in orthopyroxenite is then given by D¯gbSiδ=(3.3±3.0)×10^?9f0.0O₂e^?400±65/RT?m³s^?1, where the activation energy for diffusion is in kJ/mol, and δ is the grain boundary width in m. Calculated growth rates for enstatite under these conditions are significantly slower than predicted by an extrapolation from similar experiments performed at 1000?°C under high pressure (hydrous) conditions by Yund and Tullis (1992), perhaps due to water-enhancement of diffusion in their experiments.     

The upper thermal stability of clinochlore,Mg5Al[AlSi3O10](OH)8, at 10–35 kb P_{{\text{H}}_{\text{2}} {\text{O}}}

Clinochlore, which is, within the limits of error, the thermally most stable member of the Mg-chlorites, breaks down at = P _tot to the assemblage enstatite+forsterite+spinel+H₂O along a univariant curve located at 11 kb, 838 ° C; 15kb, 862 ° C; and 18 kb, 880 ° C (±1 kb ±10 ° C). At water pressures above that of an invariant point at 20.3 kb and 894 ° C involving the phases clinochlore, enstatite, forsterite, spinel, pyrope, and hydrous vapor, clinochlore disintegrates to pyrope+forsterite+spinel+H₂O. The resulting univariant curve has a steep, negative dP/dT slope of –930 bar/ °C at least up to 35 kb.Thus, given the proper chemical environment, Mg-chlorites have the potential of appearing as stable phases within the earth's upper mantle to maximum depths between about 60 and 100 km depending on the prevailing undisturbed geotherm, and to still greater depths in subduction zones. However, unequivocal criteria for mantle derived Mg-chlorites are difficult to find in ultrabasic rocks.     

Transformation of enstatite — diopside — jadeite pyroxenes to garnet

The high-pressure stability of enstatite(En)-diopside(Di)-jadeite(Jd) pyroxenes has been investigated experimentally with a split-sphere anvil apparatus (USSA-2000). On the enstatite-pyrope join, the compositions of garnet coexisting with enstatite were determined at 100–165 kbar and 1450–1850° C. The results indicate complete solubility between enstatite and pyrope. In the system CaO-MgO-Al₂O₃-SiO₂ (CMAS), the compositions of coexisting pyroxenes and garnet were determined at 100–165 kbar and 1250–1750° C. At 157 kbar, 1650° C, garnet with the composition En₇₉Di₂₁ (mol%) forms on the En-Di join. In the system Na₂O-MgO-Al₂O₃-SiO₂ (NMAS), the compositions of coexisting pyroxenes and garnet were determined at 60–160 kbar and 1200–1850° C. On the En-Jd join, the first garnet has the composition En₄₈Jd₅₂ at 135 kbar, 1650° C, and En₅₃Jd₄₇ at 140 kbar, 1500° C. On the Di-Jd join, the first garnet with the composition Di₆₃Jd₃₇ forms around 170 kbar, 1650° C. In the En-Di-Jd system, the first appearance of garnet with the composition En₄₂Di₉Jd₄₉ is estimated at 133 kbar, 1650° C. The new pyroxene with the composition NaMg_0.5Si_2.5O₆ (NaPx) transforms to garnet at 154 kbar, 1650° C. The experimental results indicate that the transformation of a twopyroxene assemblage to garnet and residual pyroxene in the En-Di-Jd system could occur at pressures consistent with the 400 km seismic discontinuity and in a pressure interval of 0–3 kbar.     

Reply to “Oxygen isotope heterogeneity of the mantle beneath the Canary Islands: a discussion of the paper of Gurenko et al.”

The comment by Day et al. (Contrib Mineral Petrol, 2012) (1) discusses the validity of the previously obtained oxygen isotope data for El Hierro and La Palma (Canary Island) olivines, (2) questions the approach by Gurenko et al. (Contrib Mineral Petrol 162:349–363, 2011) of using weakly correlated variations of δ¹⁸O_olivine values with X _px (proportion of pyroxenite-derived melt in the parental magma), and (3) provides reasons why oxygen isotope data by secondary ion mass spectrometry (SIMS) “offer sensitive means for detecting melt-crust interactions.” We respond these comments and report a new set of oxygen isotope measurements performed by SIMS and single-grain laser fluorination methods. These measurements confirm our previous data and conclusions and demonstrate the ability of the SIMS technique to analyze O isotopes in terrestrial samples with 2-sigma uncertainty better than ±0.25 ‰.     

Development of orthopyroxene-Fe/Mg ferrite symplectites by continuous olivine oxidation

The development of orthopyroxene-Fe/Mg ferrite symplectites associated with olivine is discussed with respect to the chemical reactions by which they form. Previously proposed reactions are presented graphically and the differences between them are reviewed. With the exception of exsolution, these are all discontinuous reactions in the sense that olivine is replaced by the two-phase symplectite assemblage.Olivine-hosted symplectites developed in the margins of lherzolite xenoliths from Kauai, Hawaii, demonstrate a reaction mechanism which has not been previously documented from natural samples. Original Fo₉₀ olivine in these samples oxidized to a new assemblage consisting of orthopyroxene (En_92–95)-Fe/Mg ferrite (Mf_35–50) symplectites developed within more magnesian olivine (Fo_92–96) hosts. Thus, by this mechanism, olivine of a different composition persists as part of a final three-phase assemblage. As oxidation advanced, the compositions of all three product phases became continuously more magnesian and the stoichiometric coefficients of the orthopyroxene and Fe/Mg ferrite continuously increased, whereas those of the product olivine decreased in the mass-balance equations. These characteristics demonstrate that the reaction was controlled by oxygen diffusion into the xenoliths from the highly oxidized alkali picrite melt in which they were entrained. Thermodynamic calculations suggest that a gradient in oxygen fugacity of 10^0.9 bars existed across the xenolith rims and resulted in compositional gradients of 4 mol% fayalite and ferrosilite and 15 mol% magnetite.     

The low-pressure stability of phase A,Mg7Si2O8(OH)6

 Phase A, Mg₇Si₂O₈(OH)₆, is a dense hydrous magnesium silicate whose importance as a host of H₂O in the Earth’s mantle is a subject of debate. We have investigated the low-pressure stability of phase A in experiments on the reaction phase A=brucite+forsterite. Experiments were conducted in piston-cylinder and multi-anvil apparatus, using mixtures of synthetic phase A, brucite and forsterite. The reaction was bracketed between 2.60 and 2.75 GPa at 500° C, between 3.25 and 3.48 GPa at 600° C and between 3.75 and 3.95 GPa at 650^° C. These pressures are much lower than observed in the synthesis experiments of Yamamoto and Akimoto (1977). At 750^° C the stability field of brucite + chondrodite was entered. The enthalpy of formation and entropy of phase A at 1 bar (10⁵ Pa), 298 K, were derived from the experimental brackets on the reaction phase A=brucite+forsterite using a modified version of the thermodynamic dataset THERMOCALC of Holland and Powell (1990), which includes a new equation of state of H₂O derived from the molecular dynamics simulations of Brodholt and Wood (1993). The data for phase A are: ΔH ^o _f =−7126±8 kJ mol^-1, S ^o=351 J K^-1 mol^-1. Incorporating these data into THERMOCALC allows the positions of other reactions involving phase A to be calculated, for example the reaction phase A + enstatite=forsterite+vapour, which limits the stability of phase A in equilibrium with enstatite. The calculated position of this reaction (753° C at 7 GPa to 937° C at 10 GPa) is in excellent agreement with the experimental brackets of Luth (1995) between 7 and 10 GPa, supporting the choice of equation of state of H₂O used in THERMOCALC. Comparison of our results with calculated P-T paths of subducting slabs (Peacock et al. 1994) suggests that, in the system MgO–SiO₂–H₂O, phase A could crystallise in compositions with Mg/Si>2 at pressures as low as 3 GPa. In less Mg rich compositions phase A could crystallise at pressures above approximately 6 GPa. Received: 3 July 1995/Accepted: 14 December 1995     

The temperature dependence of water solubility in enstatite

The solubility of water in pure enstatite was measured on samples synthesized under water-saturated conditions at 15 kbar and temperatures ranging from 700 to 1,100°C. Polarized FTIR measurements on millimetre-sized, clear crystals showed that water solubility increases strongly with temperature, from 101 ppm by weight at 700°C to 269 ppm by weight at 1,100°C. The position and shape of the infrared bands hardly changes with temperature, with one notable exception: a band close to 3,380 cm^–1 is present in samples synthesized between 700 and 1,000°C, while this band is absent from samples synthesized at 1,100°C. This effect appears to be very reproducible and points towards a slight change in the crystal structure of enstatite between 1,000 and 1,100°C at 15 kbar. The water solubility data of this study as well as those of Rauch and Keppler (Contrib Mineral Petrol 143:525–536, 2002) can be reproduced by the equation where K is water solubility, is water fugacity, A is 0.01354 ppm/bar, V^solid=12.1 cm³/mol is the volume change of enstatite during incorporation of water, and H^{1 bar}=-4,563 J/mol is the reaction enthalpy at 1 bar. This equation predicts the following behaviour of water solubility in enstatite as a function of pressure and temperature: (1) water solubility increases with pressure up to a maximum around 80 kbar; (2) water solubility decreases with temperature at 1 bar; and (3) water solubility increases with temperature between 10 and 100 kbar. If the observed temperature dependence for enstatite were representative for other upper mantle minerals as well, it would have the following implications: (1) Lateral temperature gradients in the upper mantle could cause major variations in water contents at the same depth; in particular, hot mantle plumes may scavenge water from the surrounding shallow upper mantle. (2) The scavenging of water by hot plumes could be a major factor in increasing the mobility of plumes. (3) The predicted temperature dependence of water solubility at the base of the upper mantle may allow plumes to bypass the transition zone water filter postulated by Bercovici and Karato (Nature 425:39–44, 2003).