Cr and Al diffusion in chromite spinel: experimental determination and its implication for diffusion creep

Diffusion coefficients of Cr and Al in chromite spinel have been determined at pressures ranging from 3 to 7 GPa and temperatures ranging from 1,400 to 1,700°C by using the diffusion couple of natural single crystals of MgAl₂O₄ spinel and chromite. The interdiffusion coefficient of Cr–Al as a function of Cr# (=Cr/(Cr + Al)) was determined as D _Cr–Al = D ₀ exp {−(Q′ + PV*)/RT}, where D ₀ = exp{(10.3 ± 0.08) × Cr#^0.54±0.02} + (1170 ± 31.2) cm²/s, Q′ = 520 ± 81 kJ/mol at 3 GPa, and V* = 1.36 ± 0.25 cm³/mol at 1,600°C, which is applicable up to Cr# = 0.8. The estimation of the self-diffusion coefficients of Cr and Al from Cr–Al interdiffusion shows that the diffusivity of Cr is more than one order of magnitude smaller than that of Al. These results are in agreement with patterns of multipolar Cr–Al zoning observed in natural chromite spinel samples deformed by diffusion creep.     

Compositional limits and analogs of monoclinic triple-chain silicates

Growing recognition of triple-chain silicates in nature has prompted experimental research into the conditions under which they can form and the extent of solid solution that is feasible for some key chemical substitutions. Experiments were done primarily in the range of 0.1–0.5 GPa and 200–850 °C for durations of 18–1,034 h. A wide range of bulk compositions were explored in this study that can be classified broadly into two groups: those that are Na free and involve various possible chemical substitutions into jimthompsonite (Mg₁₀Si₁₂O₃₂(OH)₄), and those that are Na bearing and involve chemical substitutions into the ideal end-member Na₄Mg₈Si₁₂O₃₂(OH)₄. Numerous attempts to synthesize jimthompsonite or clinojimthompsonite were unsuccessful despite the type of starting material used (reagent oxides, magnesite + SiO₂, talc + enstatite, or anthophyllite). Similarly, the chemical substitutions of F⁻ for OH⁻, Mn²⁺, Ca²⁺, or Fe²⁺ for Mg²⁺, and 2Li⁺ for Mg²⁺ and a vacancy were unsuccessful at nucleating triple-chain silicates. Conversely, nearly pure yields of monoclinic triple-chain silicate could be made at temperatures of 440–630 °C and 0.2 GPa from the composition Na₄Mg₈Si₁₂O₃₂(OH)₄, as found in previous studies, though its composition is most likely depleted in Na as evidenced by electron microprobe and FTIR analysis. Pure yields of triple-chain silicate were also obtained for the F-analog composition Na₄Mg₈Si₁₂O₃₂F₄ at 550–750 °C and 0.2–0.5 GPa if a flux consisting of Na-halide salt and water in a 2:1 ratio by weight was used. In addition, limited chemical substitution could be documented for the substitutions of 2 Na⁺ for Na⁺ + H⁺ and of Mg²⁺ + vacancy for 2Na⁺. For the former, the Na content appears to be limited to 2.5 cations giving the ideal composition of Na_2.5Mg₈Si₁₂O_30.5(OH)_5.5, while for the latter substitution the Na content may go as low as 1.1 cations giving the composition Na_1.1Mg_9.4Si₁₂O_31.9(OH)_4.1 based on a fixed number of Si cations. Further investigation involving Mg for Na cation exchange may provide a pathway for the synthesis of Na-free clinojimthompsonite. Fairly extensive solid solution was also observed for triple-chain silicates made along the compositional join Na₄Mg₈Si₁₂O₃₂(OH)₄–Ca₂Mg₈Si₁₂O₃₂(OH)₄ where the limit of Ca substitution at 450 °C and 0.2 GPa corresponds to Na_0.7Ca_1.8Mg_7.8Si₁₂O_31.9(OH)_4.1 (with the OH content adjusted to achieve charge balance). Aside from the Na content, this composition is similar to that observed as wide-chain lamellae in host actinolite. The relative ease with which Na-rich triple chains can be made experimentally suggests that these phases might exist in nature; this study provides additional insights into the range of compositions and formation conditions at which they might occur.     

The high PT stability of apatite and Cl partitioning between apatite and hydrous potassic phases in peridotite: an experimental study to 19 GPa with implications for the transport of P,Cl and K in the upper mantle

High PT experiments were performed in the range 2.5–19 GPa and 800–1,500°C using a synthetic peridotite doped with trace elements and OH-apatite or with Cl-apatite + phlogopite. The aim of the study was (1) to investigate the stability and phase relations of apatite and its high PT breakdown products, (2) to study the compositional evolution with P and T of phosphate and coexisting silicate phases and (3) to measure the Cl-OH partitioning between apatite and coexisting calcic amphibole, phlogopite and K-richterite. Apatite is stable in a garnet-lherzolite assemblage in the range 2.5–8.7 GPa and 800–1,100°C. The high-P breakdown product of apatite is tuite γ-Ca₃ (PO₄)₂, which is stable in the range 8–15 GPa and 1,100–1,300°C. Coexisting apatite and tuite were observed at 8 GPa/1,050°C and 8.7 GPa/1,000°C. MgO in apatite increases with P from 0.8 wt% at 2.5 GPa to 3.2 wt% at 8.7 GPa. Both apatite and tuite may contain significant Na, Sr and REE with a correlation indicating 2 Ca²⁺=Na⁺ + REE³⁺. Tuite has always higher Sr and REE and lower Fe and Mg than apatite. Phosphorus in the peridotite phases decreases in the order P_melt ≫ P_grt ≫ P_Mg2SiO4 > P_cpx > P_opx. The phosphate-saturated P₂O₅ content of garnet increases from 0.07 wt% at 2.5 GPa to 1.5 wt% at 12.8 GPa. Due to the low bulk Na content of the peridotite, ^[8]Na^[4]P^[8]M²⁺ ₋₁ ^[4]Si₋₁ only plays a minor role in controlling the phosphorus content of garnet. Instead, element correlations indicate a major contribution of ^[6]M^2+[4]P^[6]M³⁺ ₋₁ ^[4]Si₋₁. Pyroxenes contain ~200–500 ppm P and olivine has 0.14–0.23 wt% P₂O₅ in the P range 4–8.7 GPa without correlation with P, T or X_Mg. At ≥12.7 GPa, all Mg₂SiO₄ polymorphs have <200 ppm P. Coexisting olivine and wadsleyite show an equal preference for phosphorus. In case of coexisting wadsleyite and ringwoodite, the latter fractionates phosphorus. Although garnet shows by far the highest phosphorus concentrations of any peridotite silicate phase, olivine is no less important as phosphorus carrier and could store the entire bulk phosphorus budget of primitive mantle. In the Cl-apatite + phlogopite-doped peridotite, apatite contains 0.65–1.35 wt% Cl in the PT range 2.5–8.7 GPa/800–1,000°C. Apatite coexists with calcic amphibole at 2.5 GPa, phlogopite at 2.5–5 GPa and K-richterite at 7 GPa, and all silicates contain between 0.2 and 0.6 wt% Cl. No solid potassic phase is stable between 5 and 8.7 GPa. Cl strongly increases the solubility of K in hydrous fluids. This may lead to the breakdown of phlogopite and give rise to the local presence in the mantle of fluids strongly enriched in K, Cl, P and incompatible trace elements. Such fluids may get trapped as micro-inclusions in diamonds and provide bulk compositions suitable for the formation of unusual phases such as KCl or hypersilicic Cl-rich mica.     

High-temperature elasticity of magnesioferrite spinel

The elastic moduli of magnesioferrite spinel, MgFe₂O₄, and their temperature dependence have been determined for the first time by ultrasonic measurements on a polycrystalline specimen. The measurements were carried out at 300 MPa and to 700°C in a gas-medium high-pressure apparatus. On heating, both the elastic bulk (K _S) and shear (G) moduli decrease linearly to 350°C. By combining with extant thermal-expansion data, the values for the room-temperature K _S and G, and their temperature derivatives are as follows: K ₀ = 176.3(7) GPa, G ₀ = 80.1(2) GPa, (∂K _S/∂T) _P = −0.032(3) GPa K⁻¹ and (∂G/∂T) _P = −0.012(1) GPa K⁻¹. Between 350 and 400°C, there are abrupt increases of 1.4% in both of the elastic moduli; these closely coincide with the magnetic Curie transition that was observed by thermal analyses at about 360°C.     

The thermoelastic behavior of clintonite up to 10?GPa and 1,000°C

The thermoelastic behavior of a natural clintonite-1M [with composition: Ca_1.01(Mg_2.29Al_0.59Fe_0.12)_Σ3.00(Si_1.20Al_2.80)_Σ4.00O₁₀(OH)₂] has been investigated up to 10 GPa (at room temperature) and up to 960°C (at room pressure) by means of in situ synchrotron single-crystal and powder diffraction, respectively. No evidence of phase transition has been observed within the pressure and temperature range investigated. P–V data fitted with an isothermal third-order Birch–Murnaghan equation of state (BM-EoS) give V ₀ = 457.1(2) ?³, K _T0 = 76(3)GPa, and K′ = 10.6(15). The evolution of the “Eulerian finite strain” versus “normalized stress” shows a linear positive trend. The linear regression yields Fe(0) = 76(3) GPa as intercept value, and the slope of the regression line leads to a K′ value of 10.6(8). The evolution of the lattice parameters with pressure is significantly anisotropic [β(a) = 1/3K _T0(a) = 0.0023(1) GPa⁻¹; β(b) = 1/3K _T0(b) = 0.0018(1) GPa⁻¹; β(c) = 1/K _T0(c) = 0.0072(3) GPa⁻¹]. The β-angle increases in response to the applied P, with: β_P = β₀ + 0.033(4)P (P in GPa). The structure refinements of clintonite up to 10.1 GPa show that, under hydrostatic pressure, the structure rearranges by compressing mainly isotropically the inter-layer Ca-polyhedron. The bulk modulus of the Ca-polyhedron, described using a second-order BM-EoS, is K _T0(Ca-polyhedron) = 41(2) GPa. The compression of the bond distances between calcium and the basal oxygens of the tetrahedral sheet leads, in turn, to an increase in the ditrigonal distortion of the tetrahedral ring, with ∂α/∂P ≈ 0.1°/GPa within the P-range investigated. The Mg-rich octahedra appear to compress in response to the applied pressure, whereas the tetrahedron appears to behave as a rigid unit. The evolution of axial and volume thermal expansion coefficient α with temperature was described by the polynomial α(T) = α₀ + α₁ T ^−1/2. The refined parameters for clintonite are as follows: α₀ = 2.78(4) 10⁻⁵°C⁻¹ and α₁ = −4.4(6) 10⁻⁵°C^1/2 for the unit-cell volume; α₀(a) = 1.01(2) 10⁻⁵°C⁻¹ and α₁(a) = −1.8(3) 10⁻⁵°C^1/2 for the a-axis; α₀(b) = 1.07(1) 10⁻⁵°C⁻¹ and α₁(b) = −2.3(2) 10⁻⁵°C^1/2 for the b-axis; and α₀(c) = 0.64(2) 10⁻⁵°C⁻¹ and α₁(c) = −7.3(30) 10⁻⁶°C^1/2for the c-axis. The β-angle appears to be almost constant within the given T-range. No structure collapsing in response to the T-induced dehydroxylation was found up to 960°C. The HP- and HT-data of this study show that in clintonite, the most and the less expandable directions do not correspond to the most and the less compressible directions, respectively. A comparison between the thermoelastic parameters of clintonite and those of true micas was carried out.     

Synchrotron X-ray analysis of norbergite, Mg2.98Fe0.01Ti0.02Si0.99O4(OH0.31F1.69) structure at high pressure up to 8.2?GPa

The natural norbergite, Mg_2.98Fe_0.01Ti_0.02Si_0.99O₄(OH_0.31F_1.69) is examined by synchrotron X-ray diffraction analysis at pressures up to 8.2 GPa. The measured linear compressibilities of the crystallographic axes are β _a  = 2.18(4) × 10⁻³, β _b  = 2.93(7) × 10⁻³, and β _c  = 2.77(7) × 10⁻³ (GPa⁻¹), respectively and the calculated isothermal bulk modulus of the norbergite is K _T = 113(2) GPa based on the Birch–Murnaghan equation of state assuming a pressure derivative of K′ = 4. The crystal structures of norbergite are refined at room temperature and pressures of 4.7, 6.3, and 8.2 GPa, yielding R values for the structure refinements of 4.6, 5.3, and 5.3%, respectively. The bulk moduli of the polyhedral sites are 293(15) GPa for the tetrahedron, 106(5) GPa for the M2 octahedron, 113(2) GPa for the M3 octahedron, and 113(3) GPa for the total void space. The bulk modulus exhibits a good linear correlation with the filling factor for polyhedral sites in structures of the humite minerals and forsterite, reflecting the Si⁴⁺ + 4O²⁻ ⇔ □ + 4(OH, F)⁻ substitution in the humite minerals. Moreover, two simply linear trends were observed in the relationship between bulk modulus and packing index for natural minerals and dense hydrous magnesium silicate minerals. This relationship would reflect that the differences in compression mechanism were involved with hydrogen bonding in these minerals. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Behavior of epidote at high pressure and high temperature: a powder diffraction study up to 10 GPa and 1,200 K

The thermo-elastic behavior of a natural epidote [Ca_1.925 Fe_0.745Al_2.265Ti_0.004Si_3.037O₁₂(OH)] has been investigated up to 1,200 K (at 0.0001 GPa) and 10 GPa (at 298 K) by means of in situ synchrotron powder diffraction. No phase transition has been observed within the temperature and pressure range investigated. P–V data fitted with a third-order Birch–Murnaghan equation of state (BM-EoS) give V ₀ = 458.8(1)ų, K _T0 = 111(3) GPa, and K′ = 7.6(7). The confidence ellipse from the variance–covariance matrix of K _T0 and K′ from the least-square procedure is strongly elongated with negative slope. The evolution of the “Eulerian finite strain” vs “normalized stress” yields Fe(0) = 114(1) GPa as intercept values, and the slope of the regression line gives K′ = 7.0(4). The evolution of the lattice parameters with pressure is slightly anisotropic. The elastic parameters calculated with a linearized BM-EoS are: a ₀ = 8.8877(7) Å, K _T0(a) = 117(2) GPa, and K′(a) = 3.7(4) for the a-axis; b ₀ = 5.6271(7) Å, K _T0(b) = 126(3) GPa, and K′(b) = 12(1) for the b-axis; and c ₀ = 10.1527(7) Å, K _T0(c) = 90(1) GPa, and K’(c) = 8.1(4) for the c-axis [K _T0(a):K _T0(b):K _T0(c) = 1.30:1.40:1]. The β angle decreases with pressure, β_P(°) = β_P0 −0.0286(9)P +0.00134(9)P ² (P in GPa). The evolution of axial and volume thermal expansion coefficient, α, with T was described by the polynomial function: α(T) = α₀ + α₁ T ^−1/2. The refined parameters for epidote are: α₀ = 5.1(2) × 10⁻⁵ K⁻¹ and α₁ = −5.1(6) × 10⁻⁴ K^1/2 for the unit-cell volume, α₀(a) = 1.21(7) × 10⁻⁵ K⁻¹ and α₁(a) = −1.2(2) × 10⁻⁴ K^1/2 for the a-axis, α₀(b) = 1.88(7) × 10⁻⁵ K⁻¹ and α₁(b) = −1.7(2) × 10⁻⁴ K^1/2 for the b-axis, and α₀(c) = 2.14(9) × 10⁻⁵ K⁻¹ and α₁(c) = −2.0(2) × 10⁻⁴ K^1/2 for the c-axis. The thermo-elastic anisotropy can be described, at a first approximation, by α₀(a): α₀(b): α₀(c) = 1 : 1.55 : 1.77. The β angle increases continuously with T, with β_T(°) = β_T0 + 2.5(1) × 10⁻⁴ T + 1.3(7) × 10⁻⁸ T ². A comparison between the thermo-elastic parameters of epidote and clinozoisite is carried out.     

Experimental study of Pd hydrolysis in aqueous solutions at 25–70°C

The hydrolysis of the Pd²⁺ ion in HClO₄ solutions was examined at 25–70°C, and the thermodynamic constants of equilibrium K ₍₁₎⁰ and K ₍₂₎⁰were determined for the reactions Pd²⁺ + H₂O = PdOH⁺ + H⁺ and Pd²⁺ + 2H₂O = Pd(OH)₂⁰ + 2H⁺, respectively. The values of log K ₍₁₎⁰ = −1.66 ± 0.5 (25°C) and −0.65 ± 0.25 (50°C) and log K ₍₂₎⁰ = −4.34 ± 0.3 (25°C) and −3.80 ± 0.3 (50°C) were derived using the solubility technique at 0.95 confidence level. The values of log K ₍₁₎⁰ = −1.9 ± 0.6 (25°C), −1.0 ± 0.4 (50°C), and −0.5 ± 0.3 (70°C) were obtained by spectrophotometric techniques. The palladium ion is significantly hydrolyzed at elevated temperatures (50–70°C) even in strongly acidic solutions (pH 1–1.5), and its hydrolysis is enhanced with increasing temperature.     

High-pressure transitions and thermochemistry of MGeO3 (M=Mg, Zn and Sr) and Sr-silicates: systematics in enthalpies of formation of A2+B4+O3 perovskites

Phase transitions in MgGeO₃ and ZnGeO₃ were examined up to 26 GPa and 2,073 K to determine ilmenite–perovskite transition boundaries. In both systems, the perovskite phases were converted to lithium niobate structure on release of pressure. The ilmenite–perovskite boundaries have negative slopes and are expressed as P(GPa)=38.4–0.0082T(K) and P(GPa)=27.4−0.0032T(K), respectively, for MgGeO₃ and ZnGeO₃. Enthalpies of SrGeO₃ polymorphs were measured by high-temperature calorimetry. The enthalpies of SrGeO₃ pseudowollasonite–walstromite and walstromite–perovskite transitions at 298 K were determined to be 6.0±8.6 and 48.9±5.8 kJ/mol, respectively. The calculated transition boundaries of SrGeO₃, using the measured enthalpy data, were consistent with the boundaries determined by previous high-pressure experiments. Enthalpy of formation (ΔH _f°) of SrGeO₃ perovskite from the constituent oxides at 298 K was determined to be −73.6±5.6 kJ/mol by calorimetric measurements. Thermodynamic analysis of the ilmenite–perovskite transition boundaries in MgGeO₃ and ZnGeO₃ and the boundary of formation of SrSiO₃ perovskite provided transition enthalpies that were used to estimate enthalpies of formation of the perovskites. The ΔH _f° of MgGeO₃, ZnGeO₃ and SrSiO₃ perovskites from constituent oxides were 10.2±4.5, 33.8±7.2 and −3.0±2.2 kJ/mol, respectively. The present data on enthalpies of formation of the above high-pressure perovskites were combined with published data for A²⁺B⁴⁺O₃ perovskites stable at both atmospheric and high pressures to explore the relationship between ΔH _f° and ionic radii of eightfold coordinated A²⁺ (R _A) and sixfold coordinated B⁴⁺ (R _B) cations. The results show that enthalpy of formation of A²⁺B⁴⁺O₃ perovskite increases with decreasing R _A and R _B. The relationship between the enthalpy of formation and tolerance factor ( R _o: O²⁻ radius) is not straightforward; however, a linear relationship was found between the enthalpy of formation and the sum of squares of deviations of A²⁺ and B⁴⁺ radii from ideal sizes in the perovskite structure. A diagram showing enthalpy of formation of perovskite as a function of A²⁺ and B⁴⁺ radii indicates a systematic change with equienthalpy curves. These relationships of ΔH _f° with R _A and R _B can be used to estimate enthalpies of formation of perovskites, which have not yet been synthesized.     

Thermodynamic data of the high-pressure phase Mg5Al5Si6O21(OH)7 (Mg-sursassite)

 Calorimetric and P–V–T data for the high-pressure phase Mg₅Al₅Si₆O₂₁(OH)₇ (Mg-sursassite) have been obtained. The enthalpy of drop solution of three different samples was measured by high-temperature oxide melt calorimetry in two laboratories (UC Davis, California, and Ruhr University Bochum, Germany) using lead borate (2PbO·B₂O₃) at T=700 ^∘C as solvent. The resulting values were used to calculate the enthalpy of formation from different thermodynamic datasets; they range from −221.1 to −259.4 kJ mol⁻¹ (formation from the oxides) respectively −13892.2 to −13927.9 kJ mol⁻¹ (formation from the elements). The heat capacity of Mg₅Al₅Si₆O₂₁(OH)₇ has been measured from T=50 ^∘C to T=500 ^∘C by differential scanning calorimetry in step-scanning mode. A Berman and Brown (1985)-type four-term equation represents the heat capacity over the entire temperature range to within the experimental uncertainty: C _P (Mg-sursassite) =(1571.104 −10560.89×T ^−0.5−26217890.0 ×T ⁻²+1798861000.0×T ⁻³) J K⁻¹ mol⁻¹ (T in K). The P V T behaviour of Mg-sursassite has been determined under high pressures and high temperatures up to 8 GPa and 800 ^∘C using a MAX 80 cubic anvil high-pressure apparatus. The samples were mixed with Vaseline to ensure hydrostatic pressure-transmitting conditions, NaCl served as an internal standard for pressure calibration. By fitting a Birch-Murnaghan EOS to the data, the bulk modulus was determined as 116.0±1.3 GPa, (K ^′=4), V _T,0 =446.49 ³ exp[∫(0.33±0.05) × 10⁻⁴ + (0.65±0.85)×10⁻⁸ T dT], (K _T/T) _P  = −0.011± 0.004 GPa K⁻¹. The thermodynamic data obtained for Mg-sursassite are consistent with phase equilibrium data reported recently (Fockenberg 1998); the best agreement was obtained with Δ_f H ⁰ ₂₉₈ (Mg-sursassite) = −13901.33 kJ mol⁻¹, and S ⁰ ₂₉₈ (Mg-sursassite) = 614.61 J K⁻¹ mol⁻¹. Received: 21 September 2000 / Accepted: 26 February 2001     

The transition from blueschist to greenschist facies modeled by the reaction glaucophane + 2 diopside + 2 quartz = tremolite + 2 albite

The reaction glaucophane + 2 diopside + 2 quartz = tremolite + 2 albite is proposed to model the transition from the blueschist to greenschist facies. This reaction was investigated experimentally over the range of 1.0–2.1 GPa and 500–800°C using synthetic phases in the chemical system Na₂O–CaO–MgO–Al₂O₃–SiO₂–H₂O. Reversals of this reaction were possible at 500 and 550°C and growth of the low-pressure assemblage at 600°C; however, at temperatures of 600°C and higher and at pressures above 1.6 GPa omphacite nucleation (at the expense of diopside and albite) became quite strong and prevented attaining clear reversals of this reaction. Compositional changes in the amphiboles were determined by both electron microprobe analyses and correlations between unit-cell dimensions and composition. Glaucophane and particularly tremolite showed clear signs of compositional re-equilibration and merged to a single amphibole of winchite composition by about 754°C. These data were used to model the miscibility gap between glaucophane and tremolite using either the asymmetric multicomponent formulism parameters of W _TR,GL of 68 kJ with α_TR of 1.0 and α_GL of 0.75 or a simple two-site asymmetric thermodynamic mixing expression with Margules parameters W _NaCa of 13.4 kJ and W _CaNa of 19.3 kJ. Combination of these thermodynamic models of the miscibility gap with extant thermodynamic data for the other phases yields a calculated location of the above reaction, involving pure diopside and albite, that is in good agreement with the observed experimental reversals and amphibole compositions over the range of 0.94–1.93 GPa and 400–754°C. The calculated effect of jadeite solid solution into diopside is to reduce the dP/dT slope from 0.0028 to 0.0021 GPa/°C and decrease the pressure by 0.28 GPa at 754°C. The dP/dT slope of this reaction boundary lies close to a linear geotherm of 13°C/km and is consistent with the slopes of other solid–solid reactions that have been used to model the blueschist-to-greenschist facies transition.     

Location and quantitative analysis of OH in coesite

 The incorporation of hydrogen (deuterium) into the coesite structure was investigated at pressures from 3.1 to 7.5 GPa and temperatures of 700, 800, and 1100 °C. Hydrogen could only be incorporated into the coesite structure at pressures greater 5.0 GPa and 1100 °C . No correlation between the concentration of trace elements such as Al and B and the hydrogen content was observed based on ion probe analysis (1335 ± 16 H ppm and 17 ± 1 Al ppm at 7.5 GPa, 1100 °C). The FTIR spectra show three relatively intense bands at 3575, 3516, and 3459 cm⁻¹ (ν₁ to ν₃, respectively) and two very weak bands at 3296 and 3210 cm⁻¹ (ν₄ and ν₅, respectively). The band at 3516 cm⁻¹ is strongly asymmetric and can be resolved into two bands, 3528 (ν_2a) and 3508 (ν_2b) cm⁻¹, with nearly identical areas. Polarized infrared absorption spectra of coesite single-crystal slabs, cut parallel to (0 1 0) and (1 0 0), were collected to locate the OH dipoles in the structure and to calibrate the IR spectroscopy for quantitative analysis of OH in coesite (ɛ _{i ,tot}=190 000 ± 30 000 l mol⁻¹ _H2O cm⁻²). The polarized spectra revealed a strong pleochroism of the OH bands. High-pressure FTIR spectra at pressures up to 8 GPa were performed in a diamond-anvil cell to gain further insight into incorporation mechanism of OH in coesite. The peak positions of the ν₁, ν₂, and ν₃ bands decrease linearly with pressure. The mode Grüneisen parameters for ν₁, ν₂, and ν₃ are −0.074, −0.144 and −0.398, respectively. There is a linear increase of the pressure derivatives with band position which follows the trend proposed by Hofmeister et al. (1999). The full widths at half maximum (FWHM) of the ν₁, ν₂, and ν₃ bands increase from 35, 21, and 28 cm⁻¹ in the spectra at ambient conditions to 71, 68, and 105 in the 8 GPa spectra, respectively. On the basis of these results, a model for the incorporation of hydrogen in coesite was developed: the OH defects are introduced into the structure by the substitution Si^4+(Si2)+4O²⁻= ^[4]□^(Si2) + 4OH⁻, which gives rise to four vibrations, ν₁, ν_2a, ν_2b, and ν₃. Because the OH(D)-bearing samples do contain traces of Al and B, the bands ν₄ and ν₅ may be coupled to Al and/or B substitution. Received: 19 December 2000 / Accepted: 23 April 2001     

Thermal expansion behavior of Ce2Zr2O7 up to 898 K in conjunction with structural analyses by neutron diffraction

The thermal expansion of cubic pyrochlore Ce₂Zr₂O₇ has been measured from room temperature to 898 K on polycrystalline material in conjunction with structural analyses using neutron diffraction. This compound has a thermal expansion coefficient in line with the other comparable lanthanoide pyrochlore oxides. The coefficient can be expressed as α(T) = 8.418 × 10⁻⁶ + 0.9861 × 10⁻⁹ × T. The structural refinements performed for each measured temperature showed a comparable linear evolution of the Ce–O/Zr–O distances (within 0.57%).     

Hydrogen partitioning between melt,clinopyroxene, and garnet at 3 GPa in a hydrous MORB with 6 wt.% H<Subscript>2</Subscript>O

To understand partitioning of hydrogen between hydrous basaltic and andesitic liquids and coexisting clinopyroxene and garnet, experiments using a mid-ocean ridge basalt (MORB) + 6 wt.% H₂O were conducted at 3 GPa and 1,150–1,325°C. These included both isothermal and controlled cooling rate crystallization experiments, as crystals from the former were too small for ion microprobe (SIMS) analyses. Three runs at lower bulk water content are also reported. H₂O was measured in minerals by SIMS and in glasses by SIMS, Fourier Transform infrared spectroscopy (FTIR), and from oxide totals of electron microprobe (EMP) analyses. At 3 GPa, the liquidus for MORB with 6 wt.% H₂O is between 1,300 and 1,325°C. In the temperature interval investigated, the melt proportion varies from 100 to 45% and the modes of garnet and clinopyroxene are nearly equal. Liquid composition varies from basaltic to andesitic. The crystallization experiments starting from above the liquidus failed to nucleate garnets, but those starting from below the liquidus crystallized both garnet and clinopyroxene. SIMS analyses of glasses with >7 wt.% H₂O yield spuriously low concentrations, perhaps owing to hydrogen degassing in the ultra-high vacuum of the ion microprobe sample chamber. FTIR and EMP analyses show that the glasses have 3.4 to 11.9 wt.% water, whilst SIMS analyses indicate that clinopyroxenes have 1,340–2,330 ppm and garnets have 98–209 ppm H₂O. D _H ^cpx−gt is 11 ± 3, D _H ^cpx−melt is 0.023 ± 0.005 and D _H ^gt−melt is 0.0018 ± 0.0006. Most garnet/melt pairs have low values of D _H ^gt−melt, but D _H ^gt−melt increases with TiO₂ in the garnet. As also found by previous studies, values of D _H ^cpx−melt increase with Al₂O₃ of the crystal. For garnet pyroxenite, estimated values of D _H ^{pyroxenite−melt} decrease from 0.015 at 2.5 GPa to 0.0089 at 5 GPa. Hydration will increase the depth interval between pyroxenite and peridotite solidi for mantle upwelling beneath ridges or oceanic islands. This is partly because the greater pyroxene/olivine ratio in pyroxenite will tend to enhance the H₂O concentration of pyroxenite, assuming that neighboring pyroxenite and peridotite bodies have similar H₂O in their pyroxenes. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Heat capacity,entropy, and phase equilibria of dmitryivanovite

The heat capacity (C _p ) of dmitryivanovite synthesized with a cubic press was measured in the temperature range of 5–664 K using the heat capacity option of a physical properties measurement system and a differential scanning calorimeter. The entropy of dmitryivanovite at standard temperature and pressure (STP) was calculated to be 110.1 ± 1.6 J mol⁻¹ K⁻¹ from the measured C _p data. With the help of new phase equilibrium experiments done at 1.5 GPa, the phase transition boundary between krotite and dmitryivanovite was best represented by the equation: P (GPa) = −2.1825 + 0.0025 T (K). From the temperature intercept of this phase boundary and other available thermodynamic data for krotite and dmitryivanovite, the enthalpy of formation and Gibbs free energy of formation of dmitryivanovite at STP were calculated to be −2326.7 ± 2.1 and −2,208.1 ± 2.1 kJ mol⁻¹, respectively. It is also inferred that dmitryivanovite is the stable CaAl₂O₄ phase at STP and has a wide stability field at high pressures whereas the stability field of krotite is located at high temperatures and relatively low pressures. This conclusion is consistent with natural occurrences (in Ca–Al-rich inclusions) of dmitryivanovite and krotite, where the former is interpreted as the shock metamorphic product of originally present krotite.     

Thermal and redox state of the subcontinental lithospheric mantle of NE Spain from thermobarometric data on mantle xenoliths

Mantle xenoliths in within-plate Cenozoic alkaline mafic lavas from NE Spain are used to assess the local subcontinental lithospheric mantle geotherm and the influence of melting and metasomatism on its oxidation state. The xenoliths are mainly anhydrous spinel lherzolites and harzburgites and gradations between, with minor pyroxenites. Most types show protogranular textures, but transitional protogranular–porphyroclastic and equigranular lherzolites also exist. Different thermometers used in the estimates provide higher subsolidus equilibrium temperatures for harzburgites (1,062 ± 29°C) than for lherzolites (972 ± 89°C), although there is overlap; the lowest temperatures correspond to porphyroclastic lherzolites, whereas pyroxenites give the highest temperatures (up to 1,257°C). Maximum pressures for subsolidus equilibrium of peridotites are at 2.0–1.8 GPa. Later they followed adiabatic decompression and harzburgites registered lower pressures (1.02 ± 0.19 GPa) than lherzolites (1.41 ± 0.27 GPa). One pyroxenite gives values consistent with the spinel lherzolite field (1.08 GPa). The shallowest barometric data are in agreement with the highest local conductive geotherms, which implies that the lithosphere–asthenosphere boundary is at 70–60 km minimum depth. Higher equilibrium temperatures for the harzburgites could be explained by the existence of mafic magma bodies or dykes at the lower crust–mantle boundary. Paleo-fO₂ conditions during partial melting as inferred from the covariation between V and MgO concentrations are mainly between QFM−1 and QFM−2 in log units. However, most thermobarometric fO₂ estimates are between QFM−1 and QFM+1, suggesting oxidation caused by later metasomatism during uplift and cooling.     

Correlation between electron paramagnetic resonance and thermoluminescence in natural sodalite

Samples of natural sodalite, Na₈Al₆Si₆O₂₄Cl₂, submitted to gamma irradiation and to thermal treatments, have been investigated using the thermoluminescence (TL) and electron paramagnetic resonance (EPR) techniques. Both, natural and heat-treated samples at 500°C in air for 30 min, present an EPR signal around g = 2.01132 attributed to oxygen hole centers. The EPR spectra of irradiated samples show an intense line at g = 2.0008 superimposed by a hyperfine multiplet of 11 lines due to an O⁻ ion in an intermediate position with respect to two adjacent Al nuclei. In the TL measurements, the samples were annealed at 500°C for 30 min and then irradiated with γ doses varying from 0.001 to 20 kGy. All the samples have shown TL peaks at 110, 230, 270, 365, and 445°C. A correlation between the EPR g = 2.01132 line and the 365°C TL peak was observed. A TL model is proposed in which a Na⁺ ion acts as a charge compensator when an Al³⁺ ion replaces a Si⁴⁺ lattice ion. The γ ray destruction of the Al–Na complex provides an electron trapped at the Na and a hole trapped at a non-bridging oxygen ion adjacent to the Al³⁺ ion.     

Stability of clinohumite in the system MgO-SiO2-H2O

The low-pressure stability of clinohumite has been investigated in phase-equilibrium experiments on the reaction forsterite + brucite = clinohumite. The reaction was bracketed between 2.45 and 2.84 GPa at 650 °C, extending to between 1.37 and 1.57 GPa at 850 °C. At temperatures above the reaction brucite = periclase + vapour, the reaction clinohumite = forsterite + vapour was bracketed between 1.27 and 1.52 GPa at 900 °C, rising to between 1.90 and 2.00 GPa at 1000 °C. The position of the reaction forsterite + brucite = clinohumite is ∼0.5 GPa below the position determined in previous work, the difference arising either from pressure uncertainties in both studies, from enhanced reaction to clinohumite in this study due to the presence of excess brucite in the starting material, or from different concentrations of defects in the two samples. The brackets on the reaction were combined with other measured and estimated thermodynamic data for clinohumite to determine its enthalpy of formation and entropy, in a revised version of the thermodynamic dataset of Holland and Powell (1998). The values obtained were ΔH _f =−9607.29±3.05 kJ mol⁻¹, S=445 J mol⁻¹ K⁻¹. These data were used to calculate positions of other reactions involving clinohumite. The calculations suggest a larger stability field for clinohumite than implied by the results of previous experimental studies, indicating a need for more high-pressure phase-equilibrium studies to provide better thermodynamic data. Received: 30 April 1999 / Accepted: 8 November 1999     

Trace element partitioning between orthopyroxene and anhydrous silicate melt on the lherzolite solidus from 1.1 to 3.2 GPa and 1,230 to 1,535°C in the model system Na<Subscript>2</Subscript>O–CaO–MgO–Al<Subscript>2</Subscript>O<Subscript>3</Subscript>–SiO<Subscript>2</Subscript>

We determined experimentally the Nernst distribution coefficient between orthopyroxene and anhydrous silicate melt for trace elements i in the system Na₂O–CaO–MgO–Al₂O₃–SiO₂ (NCMAS) along the dry model lherzolite solidus from 1.1 GPa/1,230°C up to 3.2 GPa/1,535°C in a piston cylinder apparatus. Major and trace element composition of melt and orthopyroxene were determined with a combination of electron microprobe and ion probe analyses. We provide partitioning data for trace elements Li, Be, B, K, Sc, Ti, V, Cr, Co, Ni, Rb, Sr, Y, Zr, Nb, Cs, Ba, La, Ce, Sm, Nd, Yb, Lu, Hf, Ta, Pb, U, and Th. The melts were chosen to be boninitic at 1.1 and 2.0 GPa, picritic at 2.3 GPa and komatiitic at 2.7 and 3.2 GPa. Orthopyroxene is Tschermakitic with 8 mol% Mg-Tschermaks MgAl[AlSiO₆] at 1.1 GPa while at higher pressure it has 18–20 mol%. The rare earth elements show a continuous, significant increase in compatibility with decreasing ionic radius from D _La^opx−melt ∼ 0.0008 to D _Lu^opx−melt ∼ 0.15. For the high-field-strength elements compatibility increases from D _Th^opx−melt ∼ 0.001 through D _Nb^opx−melt ∼ 0.0015, D _U^opx−melt ∼ 0.002, D _Ta^opx−melt ∼ 0.005, D _Zr^opx−melt ∼ 0.02 and D _Hf^opx−melt ∼ 0.04 to D _Ti^opx−melt ∼ 0.14. From mathematical and graphical fits we determined best-fit values for D ₀^M1, D ₀^M2, r ₀^M1, r ₀^M2, E ₀^M1, and E ₀^M2 for the two different M sites in orthopyroxene according to the lattice strain model and calculated the intracrystalline distribution between M1 and M2. Our data indicate extreme intracrystalline fractionation for most elements in orthopyroxene; for the divalent cations D _i ^M2−M1 varies by three orders of magnitude between D _Co^M2−M1 = 0.00098–0.00919 and D _Ba^M2−M1 = 2.3–28. Trivalent cations Al and Cr almost exclusively substitute on M1 while the other trivalent cations substitute on M2; D _La^M2−M1 reaches extreme values between 6.5 × 10⁷ and 1.4 × 10¹⁶. Tetravalent cations Ti, Hf, and Zr almost exclusively substitute on M1 while U and Th exclusively substitute on M2. Our new comprehensive data set can be used for polybaric-polythermal melting models along the Earth’s mantle solidus. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

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