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     

High-temperature study and thermal expansion of phlogopite

 Unit-cell dimensions of a natural phlogopite from Pargas, Finland, have been determined in the temperature interval of 27–1050 °C by X-ray powder diffraction technique. Expansion rates vary discontinuously with temperature with a break at 412 °C. Below this temperature, the linear expansions (α) for a, b and c axis lengths are 3.74 × 10⁻⁵ K⁻¹, 1.09 × 10⁻⁵ K⁻¹, and 1.19 × 10⁻⁵ K⁻¹, respectively, and above that they are 0.86 × 10⁻⁵ K⁻¹, 0.80 × 10⁻⁵ K⁻¹, and 1.93 × 10⁻⁵ K⁻¹. The volume thermal expansion coefficients are 6.26 × 10⁻⁵ K⁻¹ and 3.71 × 10⁻⁵ K⁻¹ for low-temperature and high-temperature intervals, respectively. The observed kink in the rate of thermal expansions with temperature could be due to the different mode of structural changes. Thermogravimetric analysis of the sample indicates the oxidation of iron in the temperature range of 500–600 °C and dehydroxylation as well as decomposition of phlogopite in the temperature range of 900–1200 °C. Received: 8 September 1998 / Accepted: 28 February 2000     

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.     

Oxygen diffusion in anhydrous sanidine feldspar

Oxygen exchange experiments have been performed between single crystals of sanidine feldspar and oxygen gas enriched in ¹⁸O, at temperatures in the range 869–1053 °C, total pressure 1 atmosphere, for times up to 28 days. Oxygen isotope diffusion profiles in a direction perpendicular to (001) were determined with an ion microprobe. The experimental data obey a single Arrhenius relationship of the form D = 8.4 × 10⁻¹¹ exp. (−245 ± 15 kJ mol⁻¹/RT) m²s⁻¹. The results indicate that oxygen diffusion in anhydrous sanidine feldspar is marginally slower than oxygen diffusion in anhydrous anorthite. Comparison with published atomistic simulation studies suggests that oxygen transport in feldspar is by an interstitial mechanism. Received: 17 October 1997 / Accepted: 6 July 1998     

Partitioning of strontium between calcite,dolomite and liquids: An experimental study under higher temperature diagenetic conditions,and a model for the prediction of mineral pairs for geothermometry

The partitioning of Sr between calcite, dolomite and liquids is essentially independent of temperature between 150° and 350° C. The partition coefficients corrected for number of cation sites are b _calc=0.096 and b _dol= 0.048 for 1 mol cations/6 mol H₂O liquid. Upon dilution the partition coefficients increase, but their ratio stays constant at about 2∶1. This ratio is due to the fact that calcite has twice as many Ca-sites for Sr-substitution as dolomite. The 2∶1 relationship is also observed in natural calcite and dolomite which have undergone diagenesis. The temperature independence of partitioning is caused by the relatively small thermal expansion of calcite and dolomite. Thermal expansion between 25° and 400° C was found to follow the equations V _calc=7.0·10⁻⁴ T(°C)+36.95 and V _dol=6.9·10⁻⁴ T(°C)+32.24, V: cm³/mol. Therefore calcite and dolomite cannot serve as a temperature indicator. To have an ideal geothermometer a mineral pair with high and low thermal expansion is required. Literature date demonstrate that wurtzite, sphalerite, and galena are such minerals.     

Structural thermal behavior of paragonite and its dehydroxylate: a high-temperature single-crystal study

 The lattice constants of paragonite-2M₁, NaAl₂(AlSi₃)O₁₀(OH)₂, were determined to 800 °C by the single-crystal diffraction method. Mean thermal expansion coefficients, in the range 25–600 °C, were: α_a = 1.51(8) × 10⁻⁵, α_b = 1.94(6) × 10⁻⁵, α_c = 2.15(7) ×  10⁻⁵ °C⁻¹, and α_V = 5.9(2) × 10⁻⁵ °C⁻¹. At T higher than 600 °C, cell parameters showed a change in expansion rate due to a dehydroxylation process. The structural refinements of natural paragonite, carried out at 25, 210, 450 and 600 °C, before dehydroxylation, showed that the larger thermal expansion along the c parameter was mainly due to interlayer thickness dilatation. In the 25–600 °C range, Si,Al tetrahedra remained quite unchanged, whereas the other polyhedra expanded linearly with expansion rate proportional to their volume. The polyhedron around the interlayer cation Na became more regular with temperature. Tetrahedral rotation angle α changed from 16.2 to 12.9°. The structure of the new phase, nominally NaAl₂ (AlSi₃)O₁₁, obtained as a consequence of dehydroxylation, had a cell volume 4.2% larger than that of paragonite. It was refined at room temperature and its expansion coefficients determined in the range 25–800 °C. The most significant structural difference from paragonite was the presence of Al in fivefold coordination, according to a distorted trigonal bipyramid. Results confirm the structural effects of the dehydration mechanism of micas and dioctahedral 2:1 layer silicates. By combining thermal expansion and compressibility data, the following approximate equation of state in the PTV space was obtained for paragonite: V/V ₀ = 1 + 5.9(2) × 10⁻⁵ T(°C) − 0.00153(4) P(kbar). Received: 12 July 1999 / Revised, accepted: 7 December 1999     

Fluid source and thermal history of an epithermal vein deposit,Owen Lake,central British Columbia: evidence from bitumen and fluid inclusions

Fluid source and thermal history are determined for the barite and bitumen-bearing, early Eocene (ca. 50 Ma) polymetallic, epithermal veins of the Owen Lake deposit, central British Columbia, Canada. Carbon isotopic values for the bitumen are highly negative (δ ¹³ C ca −29%) indicating a probable terrigenous source, which may be 1 no older than Late Cretaceous or, 2 Eocene plant-fossil-bearing units stratigraphically above the Owen Lake deposit. Heat generated by suspected magmatic activity resulted in downflow of meteoric water and upflow of hydrothermal water, mixing at the site of deposition. Aqueous and hydrocarbon fluid inclusions occur within barite; Th of both types of inclusions indicate a temperature range of approximately 100° to 180 °C. Tm(ice) of aqueous inclusions range from −4.5 to −0.2 °C indicating a range of 7.2 to 0.4 equivalent weight percent NaCl. Parageneticaly younger bitumen has a vitrinite reflectance of 0.6% indicating maturation level in the temperature range of 80° to 120 °C, strongly suggesting a cooling thermal regime during barite and bitumen deposition, consistent with a late stage paragenesis.     

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.     

Mg tracer diffusion in aluminosilicate garnets at 750-850° C, 1 atm. and 1300° C, 8.5 GPa

 Tracer diffusion coefficients of Mg in natural aluminosilicate garnets of composition Alm₃₈Pyr₅₀Gr₁₀Sp₂ and Alm₇₃Pyr₂₁Gr₅Sp₁ have been measured at 1 bar, 750-850° C and at 8.5 GPa, 1300° C by chemically depositing a salt layer enriched in ²⁶Mg on the specially prepared surface of a garnet single crystal. Diffusion anneals at 1 atmosphere (101325 Pa) were carried out at a controlled f _{O 2} of ∼10^−17.5 bars maintained by a flowing gas mix of CO-CO₂. Annealing conditions were carefully chosen to avoid decomposition of garnet by redox reactions. High pressure anneals were carried out in a multianvil apparatus. Induced diffusion profiles (0.1–0.6 μm) were measured by an ion-microprobe with SIMS attachment. Diffusion coefficients at 1 atmosphere are in excellent agreement with extrapolation of data from high P-T experiments (Loomis et al. 1985; Chakraborty and Ganguly 1992) and also with the low temperature (750–900° C) dataset of Cygan and Lasaga (1985) if the diffusion coefficients are assumed to be proportional to f _{O 2} ^1/6. Such an f _{O 2} dependence, however, makes this dataset inconsistent with the recent dataset of Schwandt et al. (1995) on garnets of composition (Alm₁₅Pyr₇₂Gr₁₃Sp₀) unless a strong compositional dependence of Mg tracer diffusivity for Mg-rich garnets is invoked. The present experimental results show that such a compositional dependence is weak to non-existent for garnets with >38 mole percent almandine component. It is emphasized that the temperature dependence of diffusion coefficients at constant oxygen fugacities (activation energy ≈54 kcal/mol) are different from that along an oxygen fugacity buffer (activation energy ≈64.5 kcal/mol), as already pointed out by Chakraborty and Ganguly (1991). This distinction is of importance for modelling natural processes. The measurements at low temperatures either eliminate the need for, or greatly reduce the uncertainty of, extrapolation of laboratory data for modelling metamorphic processes. The high pressure results combined with those from Chakraborty and Ganguly (1992) and Loomis et al. (1985) indicate that pressure dependence of Mg tracer diffusivity in garnets is much stronger than that in forsterite (Chakraborty et al. 1994). This difference in pressure dependence of diffusivity may be caused by the difference in compressibility of the coordination polyhedra of Mg between olivines and garnets. Activation volumes of Mg tracer diffusion as high as 8 cm³/mol may be estimated using the present data in combination with earlier results. These data suggest that at a temperature of 1300° C, Mg tracer diffusion rates in garnets will decrease by an order of magnitude for every 100 km depth. The pressure effect will be stronger at lower temperatures. For calculations involving diffusion coefficients of garnets at high pressures (e.g. mantle xenoliths, eclogites) the pressure dependence of diffusivity must be taken into account. Received: 21 December 1994 / Accepted: 22 September 1995     

Pb diffusion in rutile

Diffusion of Pb was measured in natural and synthetic rutile under dry, 1 atmosphere conditions, using mixtures of Pb titanate or Pb sulfide and TiO₂ as the sources of diffusant. Pb depth profiles were then measured with Rutherford Backscattering Spectrometry (RBS). Over the temperature range 700–1100 °C, the following Arrhenius relation was obtained for the synthetic rutile: D=3.9 × 10⁻¹⁰exp(−250 ± 12 kJ mol⁻¹/RT) m²s⁻¹. Results for diffusion in natural and synthetic rutile were quite similar, despite significant differences in trace element compositions. Mean closure temperatures calculated from the diffusion parameters are around 600 °C for rutile grains of ∼100 μm size. This is about 100 °C higher than rutile closure temperature determinations from past field-based studies, suggesting that rutile is more resistant to Pb loss through volume diffusion than previously thought. Received: 28 June 1999 / Accepted: 29 December 1999     

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.     

The breakdown of potassium feldspar at high water pressures

The equilibrium position of the reaction between sanidine and water to form “sanidine hydrate” has been determined by reversal experiments on well characterised synthetic starting materials in a piston cylinder apparatus. The reaction was found to lie between four reversed brackets of 2.35 and 2.50 GPa at 450 °C, 2.40 and 2.59 GPa at 550 °C, 2.67 and 2.74 GPa at 650 °C, and 2.70 and 2.72 GPa at 680 °C. Infrared spectroscopy showed that the dominant water species in sanidine hydrate was structural H₂O. The minimum quantity of this structural H₂O, measured by thermogravimetric analysis, varied between 4.42 and 5.85 wt% over the pressure range of 2.7 to 3.2 GPa and the temperature range of 450 to 680 °C. Systematic variation in water content with pressure and temperature was not clearly established. The maximum value was below 6.07 wt%, the equivalent of 1 molecule of H₂O per formula unit. The water could be removed entirely by heating at atmospheric pressure to produce a metastable, anhydrous, hexagonal KAlSi₃O₈ phase (“hexasanidine”) implying that the structural H₂O content of sanidine hydrate can vary. The unit cell parameters for sanidine hydrate, measured by powder X-ray diffraction, were a = 0.53366 (±0.00022) nm and c = 0.77141 (±0.00052) nm, and those for hexasanidine were a = 0.52893 (±0.00016) nm and c = 0.78185 (±0.00036) nm. The behaviour and properties of sanidine hydrate appear to be analogous to those of the hydrate phase cymrite in the equivalent barium system. The occurrence of sanidine hydrate in the Earth would be limited to high pressure but very low temperature conditions and hence it could be a potential reservoir for water in cold subduction zones. However, sanidine hydrate would probably be constrained to granitic rock compositions at these pressures and temperatures. Received: 6 May 1997 / Accepted: 2 October 1997     

Hydrogeochemistry,environmental isotopes and the origin of the Hamamayagi-Ladik thermal spring (Samsun,Turkey)

Hamamayagi thermal spring (HTS) is located along the North Anatolian Fault Zone. The thermal spring has a temperature of 36°C, with total dissolved solids ranging from 485.6 to 508.5 mg/L. Hard, brittle, and gray limestones Permian aged are the reservoir rocks of the HTS. δ¹⁸O–δ²H isotope ratios clearly indicate a meteoric origin for the waters. The δ³⁴S value of sulfate in the thermal water is nearly 4.1‰ and implies a diagenetic environment characterized by reduced sulfur compounds. The δ¹³C ratio for dissolved inorganic carbonate in the HTS lies between −1.78 and −1.62‰, showing that it originates from the dissolution of fresh-water carbonates. Quartz geothermometry suggests a reservoir temperature of 52–85°C for the Hamamayagi geothermal field, but chalcedony geothermometers suggest reservoir temperatures between 30 and 53°C.     

Transport properties of low-sanidine single-crystals, glasses and melts at high temperature

Thermal diffusivity (D) was measured using laser-flash analysis from oriented single-crystal low-sanidine (K_0.92Na_0.08Al_0.99Fe³⁺ _0.005Si_2.95O₈), and three glasses near KAlSi₃O₈. Viscosity measurements of the three supercooled liquids, in the range 10^6.8 to 10^12.3 Pa s, confirm near-Arrhenian behavior, varying subtly with composition. For crystal and glass, D decreases with T, approaching a constant near 1,000 K: D _sat ∼ 0.65 ± 0.3 mm² s⁻¹ for bulk crystal and ∼0.53 ± 0.03 mm² s⁻¹ for the glass. A rapid decrease near 1,400 K is consistent with crossing the glass transition. Melt behavior is approximated by D = 0.475 ± 0.01 mm² s⁻¹. Thermal conductivity (k _lat) of glass, calculated using previous heat capacity (C _P) and new density data, increases with T because C _P strongly increases with T. For melt, k _lat reaches a plateau near 1.45 W m⁻¹ K⁻¹, and is always below k _lat of the crystal. Melting of potassium feldspars impedes heat transport, providing positive thermal feedback that may promote further melting in continental crust.     

Point defects and cation tracer diffusion in (Ti<Subscript><Emphasis Type="Italic">x</Emphasis></Subscript> Fe<Subscript><Emphasis Type="Italic">1−x</Emphasis></Subscript>)<Subscript>3−δ</Subscript>O<Subscript>4</Subscript>. II. Cation tracer diffusion

 Cation tracer diffusion coefficients, D_Me ^*, for Me=Fe, Mn, Co and Ti, were measured using radioactive isotopes in the spinel solid solution (Ti _x Fe _1−x )_3−δO₄ as a function of the oxygen activity. Experiments were performed at different cationic compositions (x=0, 0.1, 0.2 and 0.3) at 1100, 1200, 1300 and 1400 °C. The oxygen activity dependence of all data for D_Me ^* at constant temperature and cationic composition can be described by equations of the type D_Me ^*=D^∘ _Me[V]. C_V·a _O2 ^2/3+D_Me[I] ^∘·a _O2 ^−2/3·D_Me[V] ^∘ and D_Me[I] ^∘ are constants and C_V is a factor of the order of unity which decreases with increasing δ. All log D_Me ^* vs. loga _O2 curves obtained for different values of x and for different temperatures go through a minimum due to a change in the type of point defects dominating the cation diffusion with oxygen activity. Cation vacancies prevail for the cation diffusion at high oxygen activities while cation interstitials become dominant at low oxygen activities. At constant values of x, D_Me[V] ^∘ decreases with increasing temperature while D_Me[I] ^∘ increases.     

Formation and kinetics of porous cordierite from fly ash

Porous cordierite ceramics were prepared from a mixture of coal fly ash and basic magnesium carbonate at 1100-1350℃. Porosity, flexural strength and thermal expansion coefficient of the samples sintered at 1300℃ were estimated to be 26%, 65 MPa and 4.21×10^-6/℃, respectively. The kinetics of the formation progress was investigated by stepwise isothermal dilatometry （SID） accompanied with XRD, SEM and porosity measurement. It was found that the isothermal shrinkage data from SID could be well analyzed to get kinetic parameters according to the erapirical rate equation developed by Makipirtti-Meng, dY/dt=nk（T）Y（1-Y）（Y/1-Y）^（1/n）,where Y is the fractional shrinkage during the sintering process and n is a dimensionless component. The apparent activation energy △E values for 900-1000℃ and 1050-1 150℃ were 1294 and 1778 kJ/mol, respectively.     

Retrograde fluids in the Archean Shawmere anorthosite, Kapuskasing Structural Zone, Ontario, Canada

The Archean Shawmere anorthosite lies within the granulite facies portion of the Kapuskasing Structural Zone (KSZ), Ontario, and is crosscut by numerous linear alteration veins containing calcite + quartz ± dolomite ± zoisite ± clinozoisite ± margarite ±paragonite ± chlorite. These veins roughly parallel the trend of the Ivanhoe Lake Cataclastic Zone. Equilibria involving clinozoisite + margarite + quartz ± calcite ± plagioclase show that the vein minerals were stable at T < 600 °C, XCO₂ < 0.4 at P ≈ 6 kbar. The stabilities of margarite and paragonite in equilibrium with quartz are also consistent with T < 600 °C and XCO₂ < 0.4 at 6 kbar. Additional assemblages consisting of calcite + clinochlore + quartz + talc + margarite indicate T < 500 °C with XCO₂ > 0.9. Thus, vein formation, while clearly retrograde, spanned a range of temperatures, and fluid compositions evolved from H₂O-rich to CO₂-rich. The calcite in the retrograde veins has δ¹⁸O values that range from 8.4 to 11.2‰ (average = +9.7 ± 0.9‰) and δ¹³C values that range from −3.9 to −1.6‰ (average = −3.1 ± 0.6‰). These values indicate that the fluids from which calcite precipitated underwent extensive exchange with the anorthosite and other crustal lithologies. The fluids may have been initially derived either from devolatilization of metamorphic rocks or crystallization of igneous rocks in the adjacent Abitibi subprovince. Vein quartz contains CO₂-rich fluid inclusions (final melting T = −57.0 to −58.7 °C) that range in size from 5 to 17 μm. Measured homogenization temperatures (T h) range from −44.0 to 14.5 °C, however for most inclusions (46 of S1), T h = −44.0 to −21.1 °C (ρCO₂ ≈ 1.13 to 1.05 g/cm³). At 400 to 600 °C, these densities correspond to pressures of 3.5 to 7 kbar, which is the best estimate of pressures of vein formation. It has been argued that some high density CO₂-rich fluid inclusions found in the KSZ were formed during peak metamorphism and thus document the presence of a CO₂-rich fluid during peak granulite facies metamorphism (Rudnick et al. 1984). The association of high density CO₂-rich fluid inclusions with clearly retrograde veins documents the formation of similar composition and density inclusions after the peak of metamorphism. Thus, the coincidence of entrapment pressures calculated from fluid inclusion density measurements with peak metamorphic pressures alone should not be considered strong evidence for peak metamorphic inclusion entrapment. All fluid inclusion results are consistent with an initially semi-isobaric retrograde P–T path. Received: 2 April 1996 / Accepted: 15 November 1996     

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.     

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’}.     

The thermal expansion of ScAlO3 — A silicate perovskite analogue

The crystal structure of ScAlO₃ has been refined at temperatures up to 1100° C on the basis of x-ray powder diffraction data. The thermal expansion is adequately described by a Grüneisen-Debye model with the elastic Debye temperature and an effective Grüneisen parameter of 1.6. The volumetric thermal expansion of 3.0% between 10 and 1100° C, corresponding to a mean thermal expansion coefficient of 2.7 × 10⁻⁵ K⁻¹, is entirely attributable to the expansion of the AlO₆ octahedra. The interoctahedral angles, though not fixed by symmetry, do not vary significantly with temperature —indicating that the expansivities of the constituent AlO₆ and distorted ScO₈ polyhedra are well matched. Similar considerations of polyhedral expansivity suggest thermal expansion coefficients of ∼2 × 10⁻⁵K⁻¹ for cubic CaSiO₃ perovskite and a value between 2 × 10⁻⁵ K⁻¹ and 4 × 10⁻⁵ K⁻¹ for MgSiO₃ perovskite. The lower value is consistent with the reconnaissance measurements for Mg_0.9Fe_0.1SiO₃ (Knittle et al. 1986) below 350° C, with low-temperature measurements of single-crystal MgSiO₃ (Ross and Hazen 1989), and with the results of some recent calculations. The markedly greater expansivity ∼4 × 10⁻⁵ K⁻¹ measured at higher temperatures (350–570° C) by Knittle et al. is inconsistent with the simple Grüneisen-Debye quasiharmonic model and may reflect the marginal metastability of the orthorhombic perovskite phase. Under these circumstances, extrapolation of the measured expansivity is hazardous and may result in the under-estimation of lower mantle densities and the drawing of inappropriate inferences concerning the need for chemical stratification of the Earth's mantle.