Diffusion of tetravalent cations in zircon

Diffusion rates for the three tetravalent cations U, Th and Hf have been measured in synthetic zircon. Diffusant sources included oxide powders and ground pre-synthesized silicates. Rutherford backscattering spectrometry (RBS) was used to measure depth profiles. Over the temperature range 1400–1650 °C, the following Arrhenius relations were obtained (diffusion coefficients in m²sec⁻¹): log D _Th = (1.936 ± 0.9820) + (− 792 ± 34 kJ mol⁻¹ /2.303 RT) log D _U = (0.212 ± 2.440) + (− 726 ± 83 kJ mol⁻¹ /2.303 RT) log D _Hf = (3.206 ± 1.592) + (− 812 ± 54 kJ mol⁻¹ /2.303 RT) The data show a systematic increase in diffusivity with decreasing ionic radius (i.e., faster diffusion rates for Hf than for U or Th), a trend also observed in our earlier study of rare earth diffusion in zircon. Diffusive fractionation may be a factor in the Lu-Hf system given the much slower diffusion rates of tetravalent cations when compared with the trivalent rare earths. The very slow diffusion rates measured for these tetravalent cations suggest that they are essentially immobile under most geologic conditions, permitting the preservation of fine-scale chemical zoning and isotopic signatures of inherited cores. Received: 12 July 1996 / Accepted: 2 December 1996     

Experimental determination of the diffusion coefficient for calcium in olivine between 900°C and 1500°C

Experiments have been carried out to determine the temperature, oxygen fugacity (fO₂) and compositional dependence of the tracer diffusion coefficient (D) of calcium in olivine. These data constrain the diffusion coefficient over the temperature range 900 to 1500°C for the three principal crystallographic axes. Well constrained linear relationships between the reciprocal of the absolute temperature and log(D) exist at any given oxygen fugacity. There is a strong dependence of the diffusion coefficient on oxygen fugacity with D ∝ fO₂^(1/3). This makes a knowledge of the T-fO₂ path followed by geological samples a prerequisite for modelling Ca diffusion in olivine. The best fitting preexponential factor (Do) and activation energy (E) to the Arrhenius equation log (D) = log [Do exp(−E/RT)] + 0.31Δ log fO₂ for Ca diffusion in olivine at a given oxygen fugacity (fO₂*) are given by:diffusion along [100]: log [Do (m²/s)] = −10.78 ± 0.43; E = 193 ± 11 kJ/moldiffusion along [010]: log [Do (m²/s)] = −10.46 ± 0.37; E = 201 ± 10 kJ/moldiffusion along [001]: log [Do (m²/s)] = −10.02 ± 0.29; E = 207 ± 8 kJ/molwhere Δ log fO₂ = log[fO₂*] − log[10⁻¹²] with fO₂* in units of bars. There is no measurable compositional dependence of the diffusion coefficient between Fo₈₃ and Fo₉₂. Diffusion in Fo₁₀₀ has a much higher activation energy than in Fe-bearing olivine and has a weaker fO₂ dependence.     

Water diffusivity in rhyolitic glasses as determined by in situ IR spectroscopy

The dehydration rate of hydrous rhyolitic glasses at 475–875 °C was measured by in situ infrared (IR) spectroscopy in order to determine the diffusion coefficient of water in rhyolitic glasses. The IR spectra of glass thin sections were obtained at 90-s intervals during 90 min at high temperatures, and the change in absorbance at 3550 cm^–1 corresponding to total water was monitored. The diffusion coefficients obtained from dehydration rates of the rhyolitic glasses are considered to be averaged value over the water-concentration profile in the sample. The averaged apparent diffusion coefficients increase with the initial total water content from 0.20 m² s^–1 for 0.7 wt% to 0.37 m² s^–1 for 2.8 wt% at 700 °C. The apparent activation energy for the diffusion of total water decreases with increasing initial water content from 112 ± 6 kJ mol^–1 for 0.7 wt% to 60 ± 17 kJ mol^–1 for 4.1 wt%. Assuming a linear relation between the diffusion coefficient of total water and the total water content, the diffusion coefficients at each initial total water content were also determined. The diffusion coefficients of total water at the water contents of 0.7 and 1.9 wt% and at 0.1 MPa were best fitted by ln D=[(12.9 ± 0.8) – (111 500 ± 6400)/RT] and ln D=[(10.6 ± 0.4) – (86 800 ± 2800)/RT], respectively, and are in agreement with previous data (D in m² s^–1, T in K). The present in situ IR dehydration experiment is a rapid and effective method for the determination of water diffusivity at high temperatures.     

Oxygen diffusion in basalt and andesite melts: experimental results and discussion of chemical versus tracer diffusion

Chemical diffusion coefficients for oxygen in melts of Columbia River basalt (Ice Harbor Dam flow) and Mt. Hood andesite have been determined at 1 atm. The diffusion model is that of sorption or desorption of oxygen into a sphere of uniform initial concentration from a constant and semi-infinite atmosphere. The experimental design utilizes a thermogravimetric balance to monitor the rate of weight change arising from the response of the sample redox state to an imposed fO₂. Oxygen diffusion coefficients are approximately an order-ofmagnitude greater for basaltic melt than for andesitic melt. At 1260° C, the oxygen diffusion coefficients are: D=1.65×10^–6cm²/s and D=1.43×10^–7cm²/s for the basalt and andesite melts, respectively. The high oxygen diffusivity in basaltic melt correlates with a high ratio of nonbridging oxygen/tetrahedrally coordinated cations, low melt viscosity, and high contents of network-modifying cations. The dependence of the oxygen diffusion coefficient on temperature is: D=36.4exp(–51,600±3200/RT)cm²/s for the basalt and D=52.5exp(–60,060±4900/RT)cm²/s for the andesite (R in cal/deg-mol; T in Kelvin). Diffusion coefficients are independent of the direction of oxygen diffusion (equilibrium can be approached from extremely oxidizing or reducing conditions) and thus, melt redox state. Characteristic diffusion distances for oxygen at 1260° C vary from 10^-2 to 10² m over the time interval of 1 to 10⁶ years. A compensation diagram shows two distinct trends for oxygen chemical diffusion and oxygen tracer diffusion. These different linear relationships are interpreted as supporting distinct oxygen transport mechanisms. Because oxygen chemical diffusivities are generally greater than tracer diffusivities and their Arrhenius activation energies are less, transport mechanisms involving either molecular oxygen or vacancy diffusion are favored.     

An ion microprobe study of anhydrous oxygen diffusion in anorthite: a comparison with hydrothermal data and some geological implications

The diffusion rate of ¹⁸O tracer atoms in anorthite (An₉₇Ab₀₃) under anhydrous conditions has been measured using SIMS techniques. The tracer source was ¹⁸O₂ 98.4% gas at 1 bar, in the temperature range 1300° C–850° C. The measured diffusion constants are D ₀=1 _–0.6 ⁺¹ ×10^–9 m²s^–1 Q=236±8 kJ mol^–1 Comparison of these values with published data for ¹⁸O diffusion in anorthite under hydrothermal conditions shows that dry oxygen diffusivities are orders of magnitude lower than equivalent wet values at similar temperatures. The effect of these differences on oxygen isotope equilibration during cooling is discussed.     

Si diffusion in zircon

Self-diffusion of Si under anhydrous conditions at 1 atm has been measured in natural zircon. The source of diffusant for experiments was a mixture of ZrO₂ and ³⁰Si-enriched SiO₂ in 1:1 molar proportions; experiments were run in crimped Pt capsules in 1-atm furnaces. ³⁰Si profiles were measured with both Rutherford backscattering spectrometry (RBS) and nuclear reaction analysis with the resonant nuclear reaction ³⁰Si(p,γ)³¹P. For Si diffusion normal to c over the temperature range 1,350–1,550°C, we obtain an Arrhenius relation D = 5.8 exp(−702 ± 54 kJ mol⁻¹/RT) m² s⁻¹ for the NRA measurements, which agrees within uncertainty with an Arrhenius relation determined from the RBS measurements [62 exp(−738 ± 61 kJ mol⁻¹/RT) m² s⁻¹]. Diffusion of Si parallel to c appears slightly faster, but agrees within experimental uncertainty at most temperatures with diffusivities for Si normal to c. Diffusion of Si in zircon is similar to that of Ti, but about an order of magnitude faster than diffusion of Hf and two orders of magnitude faster than diffusion of U and Th. Si diffusion is, however, many orders of magnitude slower than oxygen diffusion under both dry and hydrothermal conditions, with the difference increasing with decreasing temperature because of the larger activation energy for Si diffusion. If we consider Hf as a proxy for Zr, given its similar charge and size, we can rank the diffusivities of the major constituents in zircon as follows: D _Zr < D _Si << D _{O, dry} < D _{O, ‘wet’}.     

Rare earth element diffusion in diopside: influence of temperature, pressure, and ionic radius, and an elastic model for diffusion in silicates

Volume diffusion rates for five rare earth elements (La, Ce, Nd, Dy, and Yb) have been measured in single crystals of natural diopside at pressures of 0.1 MPa to 2.5 GPa and temperatures of 1,050 to1,450 °C. Polished, pre-annealed crystals were coated with a thin film of rare earth element oxides, then held at constant temperature and pressure for times ranging from 20 to 882 h. Diffusion profiles in quenched samples were measured by SIMS (secondary ion mass spectrometry) depth profiling. At 1 atm pressure, with the oxygen fugacity controlled near the quartz-fayalite-magnetite buffer, the following Arrhenius relations were obtained for diffusion normal to (001) (diffusion coefficient D in m²/s): log₁₀D_Yb=(-4.64ǂ.42)-(411ᆠ kJ/mol/2.303RT); log₁₀D_Dy=(-3.31ǃ.44)-(461ᆽ kJ/mol/2.303RT); log₁₀D_Nd=(-2.95DŽ.64)-(496ᇡ kJ/mol/2.303RT); log₁₀D_Ce=(-4.10ǃ.08)-(463ᆳ kJ/mol/2.303RT); log₁₀D_Lu=(-4.22DŽ.66)-(466ᇢ kJ/mol/2.303RT). Diffusion rates decrease significantly with increasing ionic radius, with La a factor of ~35 slower than Yb. The relationship between diffusivity and ionic radius is consistent with a model in which elastic strain plays a critical role in governing the motion of an ion through the crystal lattice. Activation volumes for Yb and Ce diffusion, at constant temperature and oxygen fugacity, are 9.0DŽ.0 cm³/mol and 8.9Dž.2 cm³/mol, respectively, corresponding to an order of magnitude decrease in diffusivity as pressure is increased from 0 to 3 GPa at 1,200 °C. Diffusion of Nd is such that grain-scale isotopic equilibrium in the mantle can be achieved in ~1 My under conditions near the peridotite solidus (~1,450 °C at 2.5 GPa). The equilibration time is much longer under P, T conditions of the lithospheric mantle or at the eclogite solidus (~1 Gy at 1.5 GPa and 1,150 °C). Because of the relatively strong decrease in diffusivity with pressure (two orders of magnitude between 2.5 and 15 GPa along an adiabatic temperature gradient), Nd transport in clinopyroxene will be effectively frozen at pressures approaching the transition zone, on time scales less than 100 My. Rare earth element diffusion rates are slow enough that significant disequilibrium uptake of REE by growing clinopyroxene phenocrysts may be preserved under natural conditions of basalt crystallization. The relative abundances and spatial distributions of REE in such crystals may provide a sensitive record of the cooling and crystallization history of the host lava.     

Rare earth element diffusion in a natural pyrope single crystal at 2.8 GPa

Volume diffusion rates of Ce, Sm, Dy, and Yb have been measured in a natural pyrope-rich garnet single crystal (Py₇₁Alm₁₆Gr₁₃) at a pressure of 2.8 GPa and temperatures of 1,200-1,450 °C. Pieces of a single gem-quality pyrope megacryst were polished, coated with a thin layer of polycrystalline REE oxide, then annealed in a piston cylinder device for times between 2.6 and 90 h. Diffusion profiles in the annealed samples were measured by SIMS depth profiling. The dependence of diffusion rates on temperature can be described by the following Arrhenius equations (diffusion coefficients in m²/s): % MathType!MTEF!2!1!+- % feaaeaart1ev0aaatCvAUfKttLearuavTnhis1MBaeXatLxBI9gBam % XvP5wqSXMqHnxAJn0BKvguHDwzZbqegm0B1jxALjhiov2DaeHbuLwB % Lnhiov2DGi1BTfMBaebbfv3ySLgzGueE0jxyaibaieYlf9irVeeu0d % Xdh9vqqj-hEeeu0xXdbba9frFj0-OqFfea0dXdd9vqaq-JfrVkFHe9 % pgea0dXdar-Jb9hs0dXdbPYxe9vr0-vr0-vqpWqaaeaabiGaciaaca % qabeaadaabauaaaOqaauaabeqaeeaaaaqaaiGbcYgaSjabc+gaVjab % cEgaNnaaBaaaleaacqaIXaqmcqaIWaamaeqaaOGaemiraq0aaSbaaS % qaaiabbMfazjabbkgaIbqabaGccqGH9aqpcqGGOaakcqGHsislcqaI % 3aWncqGGUaGlcqaI3aWncqaIZaWmcqGHXcqScqaIWaamcqGGUaGlcq % aI5aqocqaI3aWncqGGPaqkcqGHsisldaqadaqaaiabiodaZiabisda % 0iabiodaZiabgglaXkabiodaZiabicdaWiaaysW7cqqGRbWAcqqGkb % GscaaMe8UaeeyBa0Maee4Ba8MaeeiBaW2aaWbaaSqabeaacqqGTaql % cqqGXaqmaaGccqGGVaWlcqaIYaGmcqGGUaGlcqaIZaWmcqaIWaamcq % aIZaWmcqWGsbGucqWGubavaiaawIcacaGLPaaaaeaacyGGSbaBcqGG % VbWBcqGGNbWzdaWgaaWcbaGaeGymaeJaeGimaadabeaakiabdseaen % aaBaaaleaacqqGebarcqqG5bqEaeqaaOGaeyypa0JaeiikaGIaeyOe % I0IaeGyoaKJaeiOla4IaeGimaaJaeGinaqJaeyySaeRaeGimaaJaei % Ola4IaeGyoaKJaeG4naCJaeiykaKIaeyOeI0YaaeWaaeaacqaIZaWm % cqaIWaamcqaIYaGmcqGHXcqScqaIZaWmcqaIWaamcaaMe8Uaee4AaS % MaeeOsaOKaaGjbVlabb2gaTjabb+gaVjabbYgaSnaaCaaaleqabaGa % eeyla0IaeeymaedaaOGaei4la8IaeGOmaiJaeiOla4IaeG4mamJaeG % imaaJaeG4mamJaemOuaiLaemivaqfacaGLOaGaayzkaaaabaGagiiB % aWMaei4Ba8Maei4zaC2aaSbaaSqaaiabigdaXiabicdaWaqabaGccq % WGebardaWgaaWcbaGaee4uamLaeeyBa0gabeaakiabg2da9iabcIca % OiabgkHiTiabiMda5iabc6caUiabikdaYiabigdaXiabgglaXkabic % daWiabc6caUiabiMda5iabiEda3iabcMcaPiabgkHiTmaabmaabaGa % eG4mamJaeGimaaJaeGimaaJaeyySaeRaeG4mamJaeGimaaJaaGjbVl % abbUgaRjabbQeakjaaysW7cqqGTbqBcqqGVbWBcqqGSbaBdaahaaWc % beqaaiabb2caTiabbgdaXaaakiabc+caViabikdaYiabc6caUiabio % daZiabicdaWiabiodaZiabdkfasjabdsfaubGaayjkaiaawMcaaaqa % aiGbcYgaSjabc+gaVjabcEgaNnaaBaaaleaacqaIXaqmcqaIWaamae % qaaOGaemiraq0aaSbaaSqaaiabboeadjabbwgaLbqabaGccqGH9aqp % cqGGOaakcqGHsislcqaI5aqocqGGUaGlcqaI3aWncqaI0aancqGHXc % qScqaIYaGmcqGGUaGlcqaI4aaocqaI0aancqGGPaqkcqGHsisldaqa % daqaaiabikdaYiabiIda4iabisda0iabgglaXkabiMda5iabigdaXi % aaysW7cqqGRbWAcqqGkbGscaaMe8UaeeyBa0Maee4Ba8MaeeiBaW2a % aWbaaSqabeaacqqGTaqlcqqGXaqmaaGccqGGVaWlcqaIYaGmcqGGUa % GlcqaIZaWmcqaIWaamcqaIZaWmcqWGsbGucqWGubavaiaawIcacaGL % Paaaaaaaaa!0C76!

Glassy orthopyroxene granodiorites of the Pannonian Basin: tracers of ultra-high-temperature deep-crustal anatexis triggered by Tertiary basaltic volcanism

Glassy orthopyroxene granodiorite-tonalite (named pincinite after type locality) was described from basaltic lapilli tuffs of the Pliocene maar near Pinciná village in the Slovakian part of the Pannonian Basin. Two pincinite types exhibit a qualitatively similar mineral composition (quartz, An_20–55 plagioclase, intergranular silicic glass with orthopyroxene and ilmenite, ±K-feldspar), but strongly different redox potential and formation PT conditions. Peraluminous pincinite is reduced (6–7% of total iron as Fe³⁺ in corundum-normative intergranular dacitic glass) and contains ilmenite with 8–10 mol% Fe₂O₃ and orthopyroxene dominated by ferrosilite. High-density (up to 0.85 g/cm³) primary CO₂ inclusions with minor H₂, CH₄, H₂S, CO and N₂ (<2 mol% total) are present in Qtz and Plg. Equilibrium PT conditions inferred from the intergranular Opx–Ilm–Glass assemblage and fluid density correspond to 1,170±50°C, 5.6±0.4 kbar, respectively. Metaluminous pincinite is more oxidised (25–27% of total iron as Fe³⁺ in diopside-normative intergranular glass of rhyolite–trachyte–dacite composition) and contains Fe₂O₃-rich ilmenite (17–29 mol%) associated with enstatite. Fluid inclusions are composed of CO₂–H₂O mixtures with up to 38 mol% H₂O. Raman spectroscopy revealed H₂S along with dominant CO₂ in the carbonic phase. Equilibrium PT parameters for the intergranular Opx–Ilm–Glass assemblage correspond to 740±15°C, 2.8±0.1 kbar, respectively. Reducing gas species (<2 mol% total) in the CO₂-inclusions of the peraluminous pincinite resulted from hydrogen diffusion due to fH₂ gradient imposed during decrease of redox potential from the log fO₂ values near QFM during Qtz + Plg growth, to QFM-2 incidental to the superimposed Opx + Ilm assemblage in the intergranular melt. The decrease in oxygen fugacity was recorded also in the metaluminous pincinite, where log fO₂ values changed from ~QFM + 2.6 to QFM + 0.4, but hydrogen diffusion did not occur. Absence of OH-bearing minerals, major and trace element abundances (e.g. REE 300–320, Nb 55–57, Th 4–31, Zr 240–300 ppm, FeO_tot/MgO up to 11), and Sr–O isotope ratios in the pincinites are diagnostic of high-temperature anorogenic magmas originated by dehydration melting of biotite in quartz-feldspathoid crust (⁸⁷Sr/⁸⁶Sr>0.705–0.706, ¹⁸O>9 V-SMOW) around alkali basalt reservoir in depths between 17 and 20 km, and around late stage derivatives of the basalt fractionation, intruding the crust up to depths of 10–11 km. Low water activity in the pincinite parental melt was caused by CO₂-flux from the Tertiary basaltic reservoirs and intrusions. The anatexis leads to generation of a melt-depleted granulitic crust beneath the Pannonian Basin, and the pincinites are interpreted as equivalents of igneous charnockites and enderbites quenched at temperatures above solidus and unaffected by sub-solidus re-equilibration and metamorphic overprint.     

Mg–Fe interdiffusion rates in wadsleyite and the diffusivity jump at the 410-km discontinuity

Mg–Fe interdiffusion rates have been measured in wadsleyite aggregates at 16.0–17.0 GPa and 1230–1530 °C by the diffusion couple method. Oxygen fugacity was controlled using the NNO buffer, and water contents of wadsleyite were measured by infrared spectroscopy. Measured asymmetric diffusion profiles, analyzed using the Boltzmann–Matano equation, indicate that the diffusion rate increases with increasing iron concentration and decreasing grain size. In the case of wadsleyite containing 50–90 weight ppm H₂O, the Mg–Fe interdiffusion coefficients at compositions of Mg/(Mg + Fe)=0.95 in the coarse-grained region (about 60 m) and 0.90 in the fine-grained region (about 6 m) were determined to be a D_Xmg = 0.95 (m² s^–1)=1.24 × 10^–9 exp[–172 (kJ mol^–1)/RT] and D_Xmg = 0.90 (m² s^–1)=1.77 × 10^–9 exp[–143 (kJ mol^–1)/RT], respectively. Grain-boundary diffusion rates were estimated to be about 4 orders of magnitude faster than the volume diffusion rate. Grain-boundary diffusion dominates when the grain size is less than a few tens of microns. Results for the nominally dry diffusion couple in the present study are roughly consistent with previous studies, taking into account differences in pressure and grain size, although water contents of samples were not clear in previous studies. We observed that the diffusivity is enhanced by about 1 order of magnitude in wadsleyite containing 300–2100 wt. ppm H₂O at 1230 °C, which is almost identical to the enhancement associated with a 300 °C increase in temperature. It is still not conclusive that a jump in diffusivity exists between olivine and wadsleyite because water contents of olivine in previous diffusion studies and effects of water on the olivine diffusivity are uncertain.     

Self-diffusion of Si and O in diopside-anorthite melt at high pressures

Self-diffusion coefficients for Si and O in Di₅₈An₄₂ liquid were measured from 1 to 4 GPa and temperatures from 1510 to 1764°C. Glass starting powders enriched in ¹⁸O and ²⁸Si were mated to isotopically normal glass powders to form simple diffusion couples, and self-diffusion experiments were conducted in the piston cylinder device (1 and 2 GPa) and in the multianvil apparatus (3.5 and 4 GPa). Profiles of ¹⁸O/¹⁶O and ^29,30Si/²⁸Si were measured using secondary ion mass spectrometry. Self-diffusion coefficients for O (D(O)) are slightly greater than self-diffusion coefficients for Si (D(Si)) and are often the same within error. For example, D(O) = 4.20 ± 0.42 × 10⁻¹¹ m²/s and D(Si) = 3.65 ± 0.37 × 10⁻¹¹ m²/s at 1 GPa and 1662°C. Activation energies for self-diffusion are 215 ± 13 kJ/mol for O and 227 ± 13 kJ/mol for Si. Activation volumes for self-diffusion are −2.1 ± 0.4 cm³/mol and −2.3 ± 0.4 cm³/mol for O and Si, respectively. The similar self-diffusion coefficients for Si and O, similar activation energies, and small, negative activation volumes are consistent with Si and O transport by a cooperative diffusion mechanism, most likely involving the formation and disassociation of a high-coordinated intermediate species. The small absolute magnitudes of the activation volumes imply that Di₅₈An₄₂ liquid is close to a transition from negative to positive activation volume, and Adam-Gibbs theory suggests that this transition is linked to the existence of a critical fraction (∼0.6) of bridging oxygen.     

The partitioning of Fe and Mg between olivine and carbonate and the stability of carbonate under mantle conditions

We have investigated the effect of Fe on the stabilities of carbonate (carb) in lherzolite assemblages by determining the partitioning of Fe and Mg between silicate (olivine; ol) and carbonates (magnesite, dolomite, magnesian calcite) at high pressures and temperatures. Fe enters olivine preferentially relative to magnesite and ordered dolomite, but Fe and Mg partition almost equally between disordered calcic carbonate and olivine. Measurement of K _d (X _Fe ^carb X _Mg ^ol /X _Fe ^ol X _Mg ^carb ) as a function of Fe/ Mg ratio indicates that Fe–Mg carbonates deviate only slightly from ideality. Using the regular solution parameter for olivine W _FeMg ^ol of 3.7±0.8 kJ/mol (Wiser and Wood 1991) we obtain for (FeMg)CO₃ a W _FeMg ^carb of 3.05±1.50 kJ/mol. The effect of Ca–Mg–Fe disordering is to raise K _d substantially enabling us to calculate W _CaMg ^carb -W _CaFe ^carb of 5.3±2.2 kJ/mol. The activity-composition relationships and partitioning data have been used to calculate the effect of Fe/Mg ratio on mantle decarbonation and exchange reactions. We find that carbonate (dolomite and magnesian calcite) is stable to slightly lower pressures (by 1 kbar) in mantle lherzolitic assemblages than in the CaO–MgO–SiO₂(CMS)–CO₂ system. The high pressure breakdown of dolomite + orthopyroxene to magnesite + clinopyroxene is displaced to higher pressures (by 2 kbar) in natural compositions relative to CMS. CO₂. We also find a stability field of magnesian calcite in lherzolite at 15–25 kbar and 750–1000°C.     

Oxygen diffusion in three silicate melts along the join diopside-anorthite

The self-diffusion of oxygen has been measured for three silicate melts along the join diopsideanorthite. The experiments were done by isotope exchange between an “infinite” reservoir of oxygen gas and spheres of melt. The oxygen self-diffusion coefficients for the three melts are given as: C-1(diopside): D = 1.64 × 10¹ exp(?(63.2 ± 20)(kcal/mole)/RT) cm²/sec C-2(Di₅₈An₄₂): D = 1.35 × 10^?1 exp(?(46.8 ± 9)(kcal/mole)/RT) cm²/sec C-3(Di₄₀An₆₀): D = 1.29 × 10^?2 exp(?(44.2 ± 6)(kcal/mole)/RT) cm²/secThe self-diffusion coefficients do not agree with the Eyring equation unless mean ionic jump distances (λ) considerably larger than the diameter of oxygen anion are assumed. However, the sense of variation of the actual diffusivities is as the Eyring equation predicts.Consideration of the results of this study and the bulk of previous work shows that oxygen appears to conform to the compensation law for cationic diffusion in silicate melts and glasses. The range of oxygen diffusivities was also found to encompass the field of divalent cation diffusivities in silicate melts.Those results imply that the diffusion of oxygen in silicate melts may involve a contribution from a cation-like diffusion mechanism (discrete O^2? anions) as well as contributions from the diffusion of larger structural units.     

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.     

Thermodynamic data from redox reactions at high temperatures. III. Activity-composition relations in Ni-Pd alloys from EMF measurements at 850–1250 K,and calibration of the NiO+Ni-Pd assemblage as a redox sensor

The assemblage NiO+Ni-Pd alloy has been calibrated as a precise oxygen fugacity sensor in the temperature range 850–1250 K at 1 bar, using an electrochemical technique with oxygen-specific CSZ electrolytes, and Ni+NiO and Cu+Cu₂O as the reference electrodes. Nine compositions were studied, ranging from 0.12 to 0.83 X _Ni ^alloy . Steady EMFs, implying equilibrium, were rapidly achieved in all cells, and were found to be reversible on increasing and decreasing temperature with a precision approaching 0.1 mV. The estimated accuracy of the measurements on each cell is ±0.2 mV (1, corresponding to ±0.003 log-bar units in fo₂ at 1273 K). Compositions of the Ni-Pd alloys were measured after each run by electron microprobe, and these compositions were then checked for internal consistency by measuring the lattice parameter by X-ray diffraction. Nickel-rich alloys show positive deviations from ideality and endothermic enthalpies of mixing, but palladium-rich compositions have exothermic enthalpies of mixing and strong negative deviations from ideality. The excess entropies of mixing are positive for all compositions, and correlate approximately with the excess volumes of mixing. The highly asymmetrical deviations from ideality are well described by a polynomial expression of the Redlich-Kister form, with three terms for the enthalpies, and two for the excess entropies and volumes of mixing. The experimental data from this study have been used to re-formulate the Ni-Pd oxygen fugacity sensor to give an expression; _O2 ^ss = _O2 ^NNO – 2RT ln X _Ni ^alloy – [2 · (1 – X _Ni ^alloy )² · [(–2165–7.958 · T) + (9409 – 0.888 · T) · (4 X _Ni ^alloy – 1) + 2089 · (6 X _Ni ^alloy – 1) · (2 X _Ni ^alloy – 1)]](850<T<1300) where _O2 ^ss is in J mol^-1, T is in kelvins, and the expression for _O2 ^NNO is that given by O'Neill and Pownceby (1993). Values in terms of log fo₂ may be obtained from the above by dividing by RT ln 10. The estimated standard error in _O2 ^ss is on the order of ±200 J mol^-1, which is approximately ±0.01 log-bar units in fo₂ at 1273 K.     

Energy associated with dislocations: A calorimetric study using synthetic quartz

High temperature oxide melt solution calorimetry was used to study the energy associated with dislocations in quartz by comparing undeformed and deformed single crystals of synthetic quartz. Samples were deformed at 698 K, 1000–1500 MPa at a strain rate of 10^–5 sec^–1. Two sets of calorimetric measurements were made: (i) using a Pt capsule as a container for powdered sample, and (ii) using pellets made from sample powder without any container. For the first set of measurements, the undeformed sample with a dislocation density of enthalpy is sum of heat content H ₉₇₃-H ₂₉₅ and enthalpy of solution in molten lead borate at 973 K of 39.22 ± 1.00 kJ mol^–1, while the sample deformed in the dislocation creep regime with a dislocation density of 6 × 10¹⁰ to 1 × 10¹¹ cm^–2 gave an enthalpy of 38.59 ± 0.78 kJ mol^–1. For the second set of measurements the measured enthalpy of the undeformed sample was 38.87 ± 0.31 kJ mol^–1, and that of a deformed sample with a dislocation density of 3 × 10¹⁰ to 1 × 10¹¹ cm^–2 was 38.24 ± 0.58 kJ mol^–1.The present study and previous theoretical calculations and estimates are consistent and suggest that the energy associated with dislocations in quartz is 0.6 ± 0.6 kJ mol^–1 for a dislocation density of 10¹¹ cm^–2; a precise value is difficult to determine because of the overlapping errors. These results indicate that for geologically realistic dislocation densities, the maximum excess energy due to dislocations would be 0.5 kJ mol^–1 for most minerals; the exact value would depend on the Burgers vector as well as the shear modulus.     

Thermodynamic properties,low-temperature heat-capacity anomalies,and single-crystal X-ray refinement of hydronium jarosite, (H<Subscript>3</Subscript>O)Fe<Subscript>3</Subscript>(SO<Subscript>4</Subscript>)<Subscript>2</Subscript>(OH)<Subscript>6</Subscript>

Crystals of hydronium jarosite were synthesized by hydrothermal treatment of Fe(III)–SO₄ solutions. Single-crystal XRD refinement with R₁=0.0232 for the unique observed reflections (|F_o| > 4_F) and wR₂=0.0451 for all data gave a=7.3559(8) Å, c=17.019(3) Å, V^o=160.11(4) cm³, and fractional positions for all atoms except the H in the H₃O groups. The chemical composition of this sample is described by the formula (H₃O)_0.91Fe_2.91(SO₄)₂[(OH)_5.64(H₂O)_0.18]. The enthalpy of formation (H^o_f) is –3694.5 ± 4.6 kJ mol^–1, calculated from acid (5.0 N HCl) solution calorimetry data for hydronium jarosite, -FeOOH, MgO, H₂O, and -MgSO₄. The entropy at standard temperature and pressure (S^o) is 438.9±0.7 J mol^–1 K^–1, calculated from adiabatic and semi-adiabatic calorimetry data. The heat capacity (C_p) data between 273 and 400 K were fitted to a Maier-Kelley polynomial C_p(T in K)=280.6 + 0.6149T–3199700T^–2. The Gibbs free energy of formation is –3162.2 ± 4.6 kJ mol^–1. Speciation and activity calculations for Fe(III)–SO₄ solutions show that these new thermodynamic data reproduce the results of solubility experiments with hydronium jarosite. A spin-glass freezing transition was manifested as a broad anomaly in the C_p data, and as a broad maximum in the zero-field-cooled magnetic susceptibility data at 16.5 K. Another anomaly in C_p, below 0.7 K, has been tentatively attributed to spin cluster tunneling. A set of thermodynamic values for an ideal composition end member (H₃O)Fe₃(SO₄)₂(OH)₆ was estimated: G^o_f= –3226.4 ± 4.6 kJ mol^–1, H^o_f=–3770.2 ± 4.6 kJ mol^–1, S^o=448.2 ± 0.7 J mol^–1 K^–1, C_p (T in K)=287.2 + 0.6281T–3286000T^–2 (between 273 and 400 K).     

F-OH exchange equilibria between mica-amphibole mineral pairs

Fluoride-hydroxyl exchange equilibria between phlogopite-pargasite and phlogopite-tremolite mineral pairs were experimentally determined at 1,173K, 500 bars and 1,073–1,173 K, 500 bars respectively. The distribution of fluorine between phlogopite and pargasite was found to favor phlogopite slightly, G _ex ^. (1,173 K)=–1.71 kJ anion^–1, while in the case of phlogopite-tremolite, fluorine was preferentially incorporated in the mica, G _ex ^. (1,073)=– 5.67 kJ anion^–1 and G _ex ^. (1,173K)=–5.84 kJ anion^–1. These results have yielded new values of entropy and Gibbs energy of formation for fluortremolite, S _f =–2,293.4±16.0JK^–1 mol^–1 and G _f = –11,779.3±25.0 kJ mol^–1, respectively. In addition, F-OH mineral exchange equilibria support a recent molten oxide calorimetric value for the Gibbs energy of fluorphlogopite, G _f =–6,014.0±7.0 kJ mol^–1, which is approximately 40 kJ mol^–1 more exothermic than the tabulated value.This work performed in part at Sandia National Laboratories supported by the U.S. Department of Energy, DOE, under contract number DE-AC04-76DP00789     

Mid-Pleistocene lavas from the Seguam volcanic center,central Aleutian arc: closed-system fractional crystallization of a basalt to rhyodacite eruptive suite

In contrast to adjacent volcanic centers of the modern central Aleutian arc, Seguam Island developed on strongly extended arc crust. K-Ar dates indicate that mid-Pleistocene, late-Pleistocene, and Holocene eruptive phases constitute Seguam. This study focuses on the petrology of the mid-Pleistocene, 1.07–07 Ma, Turf Point Formation (TPF) which is dominated by an unusual suite of porphyritic basalt and basaltic andesite lavas with subordinate phenocryst-poor andesite to rhyodacite lavas. Increasing whole-rock FeO^*/MgO from basalt to dacite, the anhydrous Plag+Ol+Cpx±Opx±Mt phenocryst assemblage, groundmass pigeonite, and the reaction Ol+Liq=Opx preserved in the mafic lavas indicate a tholeiitic affinity. Thermometry and comparison to published phase equilibria suggests that most TPF basalts crystallized Plag+Ol+Cpx±Mt at 1160°C between about 3–5 kb (±1–2% H₂O), andesites crystallized Plag+Cpx+Opx±Mt at 1000°C between 3–4 kb with 3–5% H₂O, and dacites crystallized Plag +Cpx±Opx±Mt at 1000°C between 1–2 kb with 2–3% H₂O. All lavas crystallized at f _{o 2} close to the NNO buffer. Mineral compositions and textures indicate equilibrium crystallization of the evolved lavas; petrographic evidence of open-system mixing or assimilation is rare. MgO, CaO, Al₂O₃, Cr, Ni, and Sr abundances decrease and K₂O, Na₂O, Rb, Ba, Zr, and Pb increase with increasing SiO₂ (50–71%). LREE enrichment [(Ce/Yb)_n=1.7±0.2] characterizes most TPF lavas; total REE contents increase and Eu anomalies become more negative with increasing SiO₂. Relative to other Aleutian volcanic centers, TPF basalts and basaltic andesites have lower K₂O, Na₂O, TiO₂, Rb, Ba, Sr, Zr, Y, and LREE abundances. ⁸⁷Sr/⁸⁶Sr ratios (0.70361–0.70375) and ratios of ²⁰⁶Pb/²⁰⁴Pb (18.88–18.97), ²⁰⁷Pb/²⁰⁴Pb (15.58–15.62), ²⁰⁸Pb/²⁰⁴Pb (38.46–38.55) are the highest measured for any suite of lavas in the oceanic portion of the Aleutian arc. Conversely, Nd values (+5.8 to+6.7) are among the lowest from the Aleutians. Sr, Nd, and Pb ratios are virtually constant from basalt through rhyodacite, whereas detectable isotopic heterogenity is observed at most other Aleutian volcanic centers. Major and trace element, REE, and Sr, Nd, and Pb isotopic compositions are consistent with the basaltic andesitic, andesitic, dacitic, and rhyodacitic liquids evolving from TPF basaltic magma via closed-system fractional crystallization alone. Fractionation models suggest that removal of 80 wt% cumulate (61% Plag, 17% Cpx, 12% Opx, 7% Ol, and 3% Mt) can produce 20 wt% rhyodacitic residual liquid per unit mass of parental basaltic liquid. Petrologic and physical constraints favor segregation of small batches of basalt from a larger mid-crustal reservoir trapped below a low-density upper crustal lid. In these small magma batches, the degree of cooling, crystallization, and fractionation are functions of the initial mass of basaltic magma segregated, the thermal state of the upper crust, and the magnitude of extension. Tholeiitic magmas erupted at Seguam evolved by substantially different mechanisms than did calc-alkaline lavas erupted at the adjacent volcanic centers of Kanaga and Adak on unextended arc crust. These variable differentiation mechanisms and liquid lines of descent reflect contrasting thermal and mechanical conditions imposed by the different tectonic environments in which these centers grew. At Seguam, intra-arc extension promoted eruption of voluminous basalt and its differentiates, unmodified by interaction with lower crustal or upper mantle wallrocks.     

Sub-volcanic infiltration and syn-eruptive quenching of liquids in cumulate wall-rocks: the example of the gabbroic nodules of Stromboli (Aeolian Islands, Italy)

Summary ¶Fine- to coarse-grained plutonic nodules within the Petrazza pyroclastics (Paleo-Stromboli I period) consist of gabbroic rocks with variable amounts of interstitial material. They are characterised by cumulate textures and low pressure modal mineralogy formed by plagioclase (An_96–87)+clinopyroxene (Mg-v 82–94)+olivine (Fo_83–74)±amphibole±opaque minerals; the interstitial material consists of newly crystallised microlites (quenching) of plagioclase (An_73–55)+amphibole+clinopyroxene±olivine±biotite±opaques and highly variable amounts of residual glasses that range in composition from shoshonite and high-K basaltic andesite to high-K andesite and latite. The interstitial material has a relatively high but variable degree of vesicularity. The whole rock incompatible element abundances are lower than – but the patterns are typical of – in subduction related magmas and the incompatible trace-elements are well correlated with the amount of the interstitial material. The Sr, Pb and Nd isotopic ratios resemble those of the extrusive rocks of Stromboli older series and the mineral chemistry of the gabbros is similar to that of the HKCA Paleo-Stromboli lavas. Modal mineralogy, mineral chemistry and chemical-isotopic whole rock compositions suggest that the cumulus portions of the gabbroic nodules crystallised from basaltic magmas compositionally compatible with those erupted by Stromboli volcano. The interstitial material does not represent the residual liquid after in situ crystallisation of the gabbros; it is also distinct from the juvenile host andesite magma. Textural evidence, Fe–Mg mineral/liquid partioning and mass balance calculations indicate that the interstitial material (quench crystals and vesicular glass) derived from infiltrated hydrous basaltic liquid undercooling and vesiculation of which occurred during the eruption of the Petrazza pyroclastics.Received April 17, 2002; revised version accepted November 14, 2002 Published online June 2, 2003