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     

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     

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     

Iron Diffusion in Single-Crystal Diopside

 Iron tracer diffusion experiments in diopside have been performed using natural and synthetic single crystals of diopside, and stable iron tracers enriched in ⁵⁴Fe, at temperatures in the range 950–1100 °C, total pressure 1 atm, for times up to 29 days. Iron isotope diffusion profiles were determined with an ion microprobe. For experiments performed at log pO₂ = −13, in directions parallel to the c axis and the b axis of two natural, low iron (Fe ∼ 1.8 at %) diopsides, the data obey a single Arrhenius relationship of the form D = 6.22_−5.9 ^+49.6×10⁻¹⁵ exp(−161.5 ± 35.0 kJ mol⁻¹/RT) m² s⁻¹. A single datum for iron diffusion in iron-free, single-crystal diopside at 1050 °C, is approximately 1 order of magnitude slower than in the natural crystals. The pO₂ dependence of iron diffusion in natural crystals at 1050 °C (power exponent = 0.229 ± 0.036) indicates a vacancy mechanism; this is consistent with the results of unpublished atomistic simulation studies. There is no evidence of anisotropy for iron diffusion in diopside. Received: 16 March 1999 / Accepted: 10 April 2000     

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     

Diffusion of dissolved SiO2 in H2O at 1 GPa, with implications for mass transport in the crust and upper mantle

Cation diffusion rates at 690 ± 30 °C have been calculated by inverse modelling of observed manganese (Mn) zonation profiles in 40 garnets from two kyanite-bearing metapelite samples from the High Himalayan Crystalline Series, Zanskar, northwest India. Knowledge of the initial growth profile of Mn in garnet is a pre-requisite for this technique. Following previous workers we model Mn partitioning into growing garnet in terms of a Rayleigh fractionation process, and demonstrate that the Mn_{garnet:whole rock} partition coefficient is 60–100. Three-dimensional zonation profiles were obtained by successively grinding and polishing ∼1 cm slabs of each sample at 0.1–0.2 mm intervals and analysing the garnets at each stage, thus ensuring that core sections were measured. The diffusion model assumes that garnet has a spherical geometry and behaves as a closed system, and simulates diffusive modification of the hypothetical Mn Rayleigh growth profile for each garnet. The derived measure of the time-integrated diffusion history for each garnet is then combined with radiometric and field-relation constraints for the duration of the Himalayan metamorphic event to calculate cation diffusion rates. The average cation interdiffusion rate calculated for garnets in the two samples examined is (6 ± 3.2) × 10⁻²³ m²s⁻¹. This interdiffusion rate pertains to a temperature of 690 ± 30 °C, which is 0.97 × T _PEAK, the peak temperature conditions experienced by the samples estimated using standard thermobarometric techniques. Garnet compositions are P_y2–17Alm_65–77Gro_6–16Sp_1–17. These new diffusion data are consistent with, and more precise than, existing high-temperature (>1000 °C) experimentally determined diffusion data, although some uncertainties remain difficult to constrain. Qualitative comparison between diffusively modified Mn growth profiles in garnets from the Scottish Dalradian and the Himalayan garnets suggests that the duration of metamorphism affecting the Dalradian garnets was 10–20 times longer than that endured by the Himalayan garnets. Received: 5 June 1996 / Accepted: 29 January 1997     

Columbite solubility in granitic melts: consequences for the enrichment and fractionation of Nb and Ta in the Earth's crust

The behaviour of niobium and tantalum in magmatic processes has been investigated by conducting MnNb₂O₆ and MnTa₂O₆ solubility experiments in nominally dry to water-saturated peralkaline (aluminium saturation index, A.S.I. 0.64) to peraluminous (A.S.I. 1.22) granitic melts at 800 to 1035 °C and 800 to 5000 bars. The attainment of equilibrium is demonstrated by the concurrence of the solubility products from dissolution, crystallization, Mn-doped and Nb- or Ta-doped experiments at the same pressure and temperature. The solubility products of MnNb₂O₆ (K_sp ^Nb) and MnTa₂O₆ (K_sp ^Ta) at 800 °C and 2 kbar both increase dramatically with alkali contents in water-saturated peralkaline melts. They range from 1.2 × 10⁻⁴ and 2.6 × 10⁻⁴ mol²/kg², respectively, in subaluminous melt (A.S.I. 1.02) to 202 × 10⁻⁴ and 255 × 10⁻⁴ mol²/kg², respectively, in peralkaline melt (A.S.I. 0.64). This increase from the subaluminous composition can be explained by five non-bridging oxygens being required for each excess atom of Nb⁵⁺ or Ta⁵⁺ that is dissolved into the melt. The K_sp ^Nb and K_sp ^Ta also increase weakly with Al content in peraluminous melts, ranging up to 1.7 × 10⁻⁴ and 4.6 × 10⁻⁴ mol²/kg², respectively, in the A.S.I. 1.22 composition. Columbite-tantalite solubilities in subaluminous and peraluminous melts (A.S.I. 1.02 and 1.22) are strongly temperature dependent, increasing by a factor of 10 to 20 from 800 to 1035 °C. By contrast columbite-tantalite solubility in the peralkaline composition (A.S.I. 0.64) is only weakly temperature dependent, increasing by a factor of less than 3 over the same temperature range. Similarly, K_sp ^Nb and K_sp ^Ta increase by more than two orders of magnitude with the first 3 wt% H₂O added to the A.S.I. 1.02 and 1.22 compositions, whereas there is no detectable change in solubility for the A.S.I. 0.64 composition over the same range of water contents. Solubilities are only slightly dependent on pressure over the range 800 to 5000 bars. The data for water-saturated sub- and peraluminous granites have been extrapolated to 600 °C, conditions at which pegmatites and highly evolved granites may crystallize. Using a melt concentration of 0.05 wt% MnO, 70 to 100 ppm Nb or 500 to 1400 ppm Ta are required for manganocolumbite and manganotantalite saturation, respectively. The solubility data are also used to model the fractionation of Nb and Ta between rutile and silicate melts. Predicted rutile/melt partition coefficients increase by about two orders of magnitude from peralkaline to peraluminous granitic compositions. It is demonstrated that the γNb₂O₅/γTa₂O₅ activity coefficient ratio in the melt phase depends on melt composition. This ratio is estimated to decrease by a factor of 4 to 5 from andesitic to peraluminous granitic melt compositions. Accordingly, all the relevant accessory phases in subaluminous to peraluminous granites are predicted to incorporate Nb preferentially over Ta. This explains the enrichment of Ta over Nb observed in highly fractionated granitic rocks, and in the continental crust in general. Received: 9 August 1996 / Accepted: 26 February 1997     

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.     

Fractal analysis of magnetite grains — implications for interpreting deformation mechanism

In the present study, the grain size (d) and shape of 225 magnetite grains, that crystallized at T>600°C in a syntectonic granite (Godhra Granite, India) are evaluated and implications of data to decipher deformation mechanism of magnetite are discussed. Fractal (ruler) dimension (D) analysis of magnetite grains is performed and it is demonstrated that they show fractal behaviour. Smaller magnetite grains tend to be more serrated than the larger ones, which is manifested in the higher fractal (ruler) dimension (D) of the former. Assuming a natural strain rate ranging between 10⁻¹⁰ s⁻¹ and 10⁻¹⁴ s⁻¹, the grain size data fall dominantly in the dislocation creep field of the existing deformation mechanism map of magnetite for 630°C. However, SEM-EBSD studies reveal that subgrains are absent in the magnetite grains and they did not undergo dislocation creep. Thus it is inferred that the shape of magnetite grains was not controlled by dislocation creep. It is concluded that the higher serration and increased fractal dimension of finer magnetite grains implies the importance of diffusion creep as an important deformation mechanism at high-T for magnetite in polymineralic rocks.     

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     

Cl/Br ratios and stable chlorine isotope analysis of magmatic–hydrothermal fluid inclusions from Butte,Montana and Bingham Canyon,Utah

A bulk geochemical study has been carried out on fluid inclusion leachates extracted from quartz veins from porphyry Cu deposits in Butte, Montana, USA and Bingham Canyon, Utah, USA. The leachates mostly represent low-salinity magmatic–hydrothermal fluid inclusions. Their halogen ratios (Br/Cl) of fluid inclusion leachates were determined by ion chromatography, and δ³⁷Cl values of the leachates were measured by continuous-flow isotope ratio mass spectrometry. Br/Cl ratios from early pre-Main stage and later Main stage veins at Butte range from 0.60 to 1.88 × 10⁻³ M. Ratios are similar in pre-Main stage veins with sericite bearing selvages and Main stage samples ranging from 0.81 to 1.08 × 10⁻³ and from 0.92 to 1.88 × 10⁻³ M, respectively, clustering below seawater (1.54 × 10⁻³ M) and overlapping mantle values (~1–2 × 10⁻³ M). Two samples associated with early pre-Main stage potassic alteration yield distinctly lower Br/Cl ratios of 0.60 and 0.64 × 10⁻³ M. Butte δ³⁷Cl values range from −0.8‰ to −2.3‰ with no significant difference between pre-Main stage and Main stage samples. Br/Cl ratios for quartz veins from Bingham Canyon range from 0.18 to 3.68 × 10⁻³ M. Br/Cl ratios from Bingham range above and below previously reported for porphyry copper deposits. In contrast to Butte, δ³⁷Cl values for Bingham are lower, ranging from −0.9‰ to −4.1‰. In the absence of any processes which can significantly fractionate chlorine isotopes at high temperatures, we suggest that the porphyry system at Bingham, and to a lesser extent at Butte, have inherited negative chlorine isotopic signatures from the subducting slab generated at low temperatures.     

Temperature evolution between 50 K and 320 K of the thermal expansion tensor of gypsum derived from neutron powder diffraction data

The unit cell parameters, extracted from Rietveld analysis of neutron powder diffraction data collected between 4.2 K and 320 K, have been used to calculate the temperature evolution of the thermal expansion tensor for gypsum for 50 ≤ T ≤ 320 K. At 300 K the magnitudes of the principal axes are α ₁₁  = 1.2(6) × 10⁻⁶ K⁻¹, α ₂₂  = 36.82(1) × 10⁻⁶ K⁻¹ and α ₃₃  = 25.1(5) × 10⁻⁶ K⁻¹. The maximum axis, α ₂₂ , is parallel to b, and using Institution of Radio Engineers (IRE) convention for the tensor orthonormal basis, the axes α ₁₁ and α ₃₃ have directions equal to (−0.979, 0, 0.201) and (0.201, 0, 0.979) respectively. The orientation and temperature dependent behaviour of the thermal expansion tensor is related to the crystal structure in the I2/a setting. Received 12 February 1998 / Revised, accepted 19 October 1998     

The enthalpy of transformation of Ca(OH)2-I (portlandite) to Ca(OH)2-II (EuI2 structure) by low-temperature DSC

 Thermodynamic properties of high-pressure minerals that are not recoverable from synthesis experiments by conventional quenching methods (“unquenchable” phases) usually are calculated from equation of state data and phase diagram topologies. The present study shows that, with cryogenic methods of recovery and sample treatment, phases with a suitable decomposition rate can be made accessible to direct thermodynamic measurements. A set of samples of Ca(OH)₂-II has been synthesized in a multianvil device and subsequently recovered by cooling the high-pressure assembly with liquid nitrogen. Upon heating from liquid nitrogen to room temperature, the material transformed back to Ca(OH)₂-I. The heat effect of this backtransformation was measured by differential scanning calorimetry. A commercial differential scanning calorimeter (Netzsch DSC 404), modified to allow sample loading at liquid nitrogen temperature was used to heat the material from −150 to +200 °C at rates varying between 5 and 15 °C min⁻¹. The transformation started around −50 °C very gradually, and peaked at about 0 °C. To obtain a baseline correction, each sample was scanned under exactly the same conditions after the backtransformation was complete. Because of the relative sluggishness, onset and offset temperatures were not well defined as compared to fast (e.g., melting) reactions. To aid in integration, the resulting signals were successfully fitted using a generic asymmetric peak model. The enthalpy of backtransformation was determined to be ΔH =−10.37 ± 0.50 kJ mol⁻¹. From previous in situ X-ray diffraction experiments, the location of the direct transformation in P-T space has been constrained to 5.7 ± 0.4 GPa at 500 °C (Kunz et al. 1996). With the reaction volume known from the same study, and assuming that ΔC _p of the transformation remains negligible between the conditions of our measurements and 500 °C, our result gives an estimate of the entropy of transition and the P-T slope of the reaction curve. To a first approximation, the values ΔS = −16.00 ± 0.65 J(mol · K)⁻¹ and dP/dT = 0.0040 ± 0.0002 GPa/K have been determined. These results need to be refined by equation of state data for Ca(OH)₂-II. Received: 30 December 1999 / Accepted: 10 April 2000     

Origin of the Zálesí U–Ni–Co–As–Ag/Bi deposit,Bohemian Massif,Czech Republic: fluid inclusion and stable isotope constraints

The Zálesí vein-type deposit is hosted by Early Paleozoic high-grade metamorphic rocks on the northern margin of the Bohemian Massif. The mineralization is composed of three main stages: uraninite, arsenide, and sulfide. The mineral assemblages formed at low temperatures (~80 to 130°C, locally even lower) and low pressures (<100 bars). The salinity of the aqueous hydrothermal fluids (0 to 27 wt.% salts) and their chemical composition vary significantly. Early fluids of the oldest uraninite stage contain a small admixture of a clathrate-forming gas, possibly CO₂. Salinity correlates with oxygen isotope signature of the fluid and suggests mixing of brines [δ ¹⁸O around +2‰ relative to standard mean ocean water (SMOW)] with meteoric waters (δ ¹⁸O around −4‰ SMOW). The fluid is characterized by highly variable halogen ratios (molar Br/Cl = 0.8 × 10⁻³ to 5.3 × 10⁻³; molar I/Cl = 5.7 × 10⁻⁶ to 891 × 10⁻⁶) indicating a dominantly external origin for the brines, i.e., from evaporated seawater, which mixed with iodine-enriched halite dissolution brine. The cationic composition of these fluids indicates extensive interaction of the initial brines with their country rocks, likely associated with leaching of sulfur, carbon, and metals. The brines possibly originated from Permian–Triassic evaporites in the neighboring Polish Basin, infiltrated into the basement during post-Variscan extension and were finally expelled along faults giving rise to the vein-type mineralization. Cenozoic reactivation by low-salinity, low-δ ¹⁸O (around −10‰ SMOW) fluids of mainly meteoric origin resulted in partial replacement of primary uraninite by coffinite-like mineral aggregates.     

Strain-rate estimation using fractal analysis of quartz grains in naturally deformed rocks

The area-perimeter fractal dimension (D) of quartz grains has earlier been proposed as a strain-rate gauge based on experimental deformation of quartz aggregates. To test the application in naturally deformed rocks, D is calculated in (a) three quartzites belonging to the Lunavada Group of rocks (Aravalli Mountain Belt, NW India) that developed textures between 420–600°C and (b) one quartz reef sample from the Malanjkhand Granite (Central India), which underwent dynamic recrystallization between 250–400°C. Using the above T ranges and calculated D values, strain-rates are calculated for the two sets of samples. A 10^−12.7 s⁻¹ strain rate at 250°C is calculated for the quartz reef sample. However, at higher temperatures the calculated strain-rate is >10⁻¹⁰ s⁻¹ for the quartz reef and the quartzite samples. The quartzites show evidence of dynamic recrystallization by grain boundary migration (GBM) and subgrain rotation (SGR), while the quartz reef is replete with evidence of bulging (BLG) recrystallization. T and calculated strainrates are plotted on available recrystallization map of quartz. It is demonstrated that whilst the T/strain-rate of the quartzites does not fall in the region of GBM and SGR, the T/strain-rate of the quartz reef falls in the BLG region. The problems with strain-rate calculations using area perimeter fractal dimension are discussed. It is concluded that the method of strain-rate calculation can be used only for lower T.     

Supercooled diopside melt: confirmation of temperature-dependent expansivity using container-based dilatometry

The expansivity of supercooled diopside liquid has been determined using techniques of container-based dilatometry. Two thermal strategies have been employed, one in which the sample is brought to volumetric equilibrium by long-duration dwells at low temperatures (817 °C) and one in which scanning dilatometry of the sample has been performed at somewhat higher temperatures (890–913 °C). The results of both experiments yield a supercooled liquid expansivity for diopside liquid in the temperature range of 817–913 °C of 84.4 ± 2.8 × 10⁻⁴ cm³/mol K. The expansivity is 65% higher than that obtained for diopside melt obtained at superliquidus temperatures using the double bob Archimedean method. Combined fitting of the new low temperature, volume–temperature data from the present study and the superliquidus data from the literature has been performed. The combined fit yields the following equations for the volume–temperature relationship of diopside liquid (T=temperature in °C):

Magnesium self-diffusion in orthoenstatite

Magnesium self-diffusion coefficients were determined experimentally for diffusion parallel to each of the three crystallographic directions in natural orthoenstatite (En₈₈Fs₁₂). Experiments were conducted at 1 atm in CO-CO₂ gas mixing furnaces, which provided oxygen fugacities equivalent to the iron-wüstite buffer. Diffusion of ²⁵Mg was induced in polished samples of oriented orthoenstatite using a film of isotopically enriched ²⁵MgO as the source material. Very short (<0.15 μm) diffusional penetration profiles were measured by ion microprobe depth profiling. The diffusion coefficients determined for four temperatures (900, 850, 800, 750 °C) provide the activation energies, E _a , and frequency factors, D _o, where D = D _o exp (−E _a /RT) for Mg self-diffusion parallel to each crystallographic direction: a-axis, E _a  = 360 ± 52 kJ/mole and D _o = 1.10 × 10⁻⁴ m²/s; b-axis, E _a  = 339 ± 77 kJ/mole and D _o = 6.93 × 10⁻⁶ m²/s and c-axis, E _a  = 265 ± 66 kJ/mole and D _o = 4.34 × 10⁻⁹ m²/s. In this temperature range, any possible anisotropy of cation diffusion is very small, however the activation energy for diffusion parallel to the c-axis (001) is the lowest and the activation energies for diffusion parallel to the a-axis (100) and b-axis (010) are higher. Application of these diffusion results to the silicate phases of the Lowicz mesosiderite meteorite provides cooling rates for the silicate portion of the meteorite (4–11 °C/100 years) that are similar, although slower, to previous estimates. These silicate cooling rates are still several orders of magnitude faster than the cooling rates (0.1 °C/10⁶ years) for the metal portions. Received: 22 January 1997 / Accepted: 2 October 1997     

An experimental study of effect of pre‐existing fabric on deformation of foliated mylonite at high temperature and pressure

We performed deformation experiments on a foliated mylonite under high temperature and pressure conditions in this study. To investigate the effect of pre‐existing fabric on the rheology of rocks, our samples were drilled from natural mylonite with the cylinder axis parallel to the foliation (PAR) and perpendicular to the foliation (PER). We performed 25 tests on seven PAR samples and 21 tests on seven PER samples at temperatures ranging from 600 to 890 °C, confining pressures ranging from 800 to 1400 MPa, and steady‐state strain rates of 1 × 10⁻⁴, 1 × 10⁻⁵ and 2.5 × 10⁻⁶ s⁻¹. In the temperatures of 600–700 °C, the deformation is accommodated by semi‐brittle flow, with the average stress exponent being 6–7 assuming power law flow; in the temperature range of 800–890 °C, deformation is mainly by plastic flow, with an average stress exponent of n = 3 and activation energies of Q = 354 ± 52 kJ/mol (PER and PAR samples). The experimental results show that the strengths of PER samples are higher than those of PAR samples. Deformation microstructures have been studied by optical and electron microscopy. The original foliation of PER samples is destroyed by deformation and replaced by a new foliation, but the deformation of PAR samples followed the original foliation. Electron backscatter diffraction (EBSD) measurements show a strong lattice preferred orientation (LPO) of the quartz c axis fabrics of the starting samples and deformed PER and PAR samples. However, the c axis fabric of quartz in experimentally deformed PER and PAR samples varied with temperature and strain rate is different from that seen in the starting mylonite sample. The initial quartz c axis fabric of the starting mylonite sample has been transformed into a new fabric during experimental deformation. Dehydration melting of biotite and hornblende occurred in both PER and PAR samples at temperatures of 800–890 °C. Copyright © 2014 John Wiley & Sons, Ltd.     

Thermal decomposition of Ferrian chamosite: an infrared emission spectroscopic study

The thermal behaviour of ripidolite, an iron-rich chlorite, has been studied in situ by infrared emission spectroscopy up to 800 °C. The more di,trioctahedral nature due to significant amounts of Fe³⁺ is reflected, in addition to the two bands around 3420 and 3560 cm⁻¹, by an extra band around 3345 cm⁻¹. This extra band is absent in pure dioctahedral chlorites without Fe³⁺. These bands have been assigned to (AlAl)O-OH, (SiAl)O-OH and (SiSi)O-OH stretching modes with increasing frequencies. The bands disappear upon dehydroxylation around 650 °C. A similar behaviour is observed for the corresponding libration modes around 716, 759 and 802 cm⁻¹. The stretching and bending modes of the inner-OH of the octahedral sheet in the 2:1 clay-like layer are observed around 3645, 943 and 904 cm⁻¹. Although the bands decrease in intensity, they remain present up to 800 °C as dehydroxylation of the octahedral sheet is not yet complete at this temperature. The presence of two bending modes is explained as being due to a differentiation between Mg-OH and Fe-OH modes. At 650 °C a new sharp band is observed around 502 cm⁻¹ assigned to a (Fe,Mg)-O-Al bending mode caused by the formation of a spinel-like interlayer phase after dehydroxylation. Received: 4 June 1999 / Accepted: 6 August 1999