Ephesite,Na(LiAl2) [Al2Si2O10] (OH)2: I. thermal stability and standard state thermodynamic properties

Ephesite, Na(LiAl₂) [Al₂Si₂O₁₀] (OH)₂, has been synthesized for the first time by hydrothermal treatment of a gel of requisite composition at 300≦T(° C)≦700 and \(P_{H_2 O}\) upto 35 kbar. At \(P_{H_2 O}\) between 7 and 35 kbar and above 500° C, only the 2M₁ polytype is obtained. At lower temperatures and pressures, the 1M polytype crystallizes first, which then inverts to the 2M₁ polytype with increasing run duration. The X-ray diffraction patterns of the 1M and 2M₁ poly types can be indexed unambiguously on the basis of the space groups C2 and Cc, respectively. At its upper thermal stability limit, 2M₁ ephesite decomposes according to the reaction (1) $$\begin{gathered} {\text{Na(LiAl}}_{\text{2}} {\text{) [Al}}_{\text{2}} {\text{Si}}_{\text{2}} {\text{O}}_{{\text{10}}} {\text{] (OH)}}_{\text{2}} \hfill \\ {\text{ephesite}} \hfill \\ {\text{ = Na[AlSiO}}_{\text{4}} {\text{] + LiAl[SiO}}_{\text{4}} {\text{] + }}\alpha {\text{ - Al}}_{\text{2}} {\text{O}}_{\text{3}} {\text{ + H}}_{\text{2}} {\text{O}} \hfill \\ {\text{nepheline }}\alpha {\text{ - eucryptite corundum}} \hfill \\ \end{gathered}$$ Five reversal brackets for (1) have been established experimentally in the temperature range 590–750° C, at \(P_{H_2 O}\) between 400 and 2500 bars. The equilibrium constant, K, for this reaction may be expressed as (2) $$log K{\text{ = }}log f_{{\text{H}}_{\text{2}} O}^* = 7.5217 - 4388/T + 0.0234 (P - 1)T$$ where \(f_{H_2 O}^* = f_{H_2 O} (P,T)/f_{H_2 O}^0\) (1,T), with T given in degrees K, and P in bars. Combining these experimental data with known thermodynamic properties of the decomposition products in (1), the following standard state (1 bar, 298.15 K) thermodynamic data for ephesite were calculated: H _f,298.15 ⁰ =-6237372 J/mol, S _298.15 ⁰ =300.455 J/K·mol, G _298.15 ⁰ =-5851994 J/mol, and V _298.15 ⁰ =13.1468 J/bar·mol.     

An experimental study of Ti and Zr partitioning among zircon, rutile, and granitic melt

In order to evaluate the effect of trace and minor elements (e.g., P, Y, and the REEs) on the high-temperature solubility of Ti in zircon (zrc), we conducted 31 experiments on a series of synthetic and natural granitic compositions [enriched in TiO₂ and ZrO₂; Al/(Na + K) molar ~1.2] at a pressure of 10 kbar and temperatures of ~1,400 to 1,200 °C. Thirty of the experiments produced zircon-saturated glasses, of which 22 are also saturated in rutile (rt). In seven experiments, quenched glasses coexist with quartz (qtz). SiO₂ contents of the quenched liquids range from 68.5 to 82.3 wt% (volatile free), and water concentrations are 0.4–7.0 wt%. TiO₂ contents of the rutile-saturated quenched melts are positively correlated with run temperature. Glass ZrO₂ concentrations (0.2–1.2 wt%; volatile free) also show a broad positive correlation with run temperature and, at a given T, are strongly correlated with the parameter (Na + K + 2Ca)/(Si·Al) (all in cation fractions). Mole fraction of ZrO₂ in rutile $ \left( {\mathop X\nolimits_{{{\text{ZrO}}_{ 2} }}^{\text{rt}} } \right) $ in the quartz-saturated runs coupled with other 10-kbar qtz-saturated experimental data from the literature (total temperature range of ~1,400 to 675 °C) yields the following temperature-dependent expression: $ {\text{ln}}\left( {\mathop X\nolimits_{{{\text{ZrO}}_{ 2} }}^{\text{rt}} } \right) + {\text{ln}}\left( {a_{{{\text{SiO}}_{2} }} } \right) = 2.638(149) - 9969(190)/T({\text{K}}) $ , where silica activity $ a_{{{\text{SiO}}_{2} }} $ in either the coexisting silica polymorph or a silica-undersaturated melt is referenced to α-quartz at the P and T of each experiment and the best-fit coefficients and their uncertainties (values in parentheses) reflect uncertainties in T and $ \mathop X\nolimits_{{{\text{ZrO}}_{2} }}^{\text{rt}} $ . NanoSIMS measurements of Ti in zircon overgrowths in the experiments yield values of ~100 to 800 ppm; Ti concentrations in zircon are positively correlated with temperature. Coupled with values for $ a_{{{\text{SiO}}_{2} }} $ and $ a_{{{\text{TiO}}_{2} }} $ for each experiment, zircon Ti concentrations (ppm) can be related to temperature over the range of ~1,400 to 1,200 °C by the expression: $ \ln \left( {\text{Ti ppm}} \right)^{\text{zrc}} + \ln \left( {a_{{{\text{SiO}}_{2} }} } \right) - \ln \left( {a_{{{\text{TiO}}_{2} }} } \right) = 13.84\left( {71} \right) - 12590\left( {1124} \right)/T\left( {\text{K}} \right) $ . After accounting for differences in $ a_{{{\text{SiO}}_{2} }} $ and $ a_{{{\text{TiO}}_{2} }} $ , Ti contents of zircon from experiments run with bulk compositions based on the natural granite overlap with the concentrations measured on zircon from experiments using the synthetic bulk compositions. Coupled with data from the literature, this suggests that at T ≥ 1,100 °C, natural levels of minor and trace elements in “granitic” melts do not appear to influence the solubility of Ti in zircon. Whether this is true at magmatic temperatures of crustal hydrous silica-rich liquids (e.g., 800–700 °C) remains to be demonstrated. Finally, measured $ D_{\text{Ti}}^{{{\text{zrc}}/{\text{melt}}}} $ values (calculated on a weight basis) from the experiments presented here are 0.007–0.01, relatively independent of temperature, and broadly consistent with values determined from natural zircon and silica-rich glass pairs.     

A multi-isotope approach for understanding sources of water,carbon and sulfur in natural springs of the Central Appalachian region

Natural springs have been reliable sources of domestic water and have allowed for the development of recreational facilities and resorts in the Central Appalachians. The structural history of this area is complex and it is unknown whether these natural springs receive significant recharge from modern precipitation or whether they discharge old water recharged over geological times scales. The main objective of this study was to use stable isotopes of water ( $\delta^{18} {\text{O}}_{{{\text{H}}_{2} {\text{O}}}}$ and $\delta^{2} {\text{H}}_{{{\text{H}}_{2} {\text{O}}}}$ ), dissolved inorganic carbon ( $\delta^{13} {\text{C}}_{\text{DIC}}$ ) and dissolved sulfate ( $\delta^{34} {\text{S}}_{{{\text{SO}}_{4} }}$ and $\delta^{18} {\text{O}}_{{{\text{SO}}_{4} }}$ ) to delineate sources of water, carbon and sulfur in several natural springs of the region. Our preliminary isotope data indicate that all springs are being recharged by modern precipitation. The oxygen isotope composition indicates that waters in thermal springs did not encounter the high temperatures required for O isotope exchange between the water and silicate/carbonate minerals, and/or the residence time of water in the aquifers was short due to high flow rates. The carbon isotopic composition of dissolved inorganic carbon and sulfur/oxygen isotopic composition of dissolved sulfate provide evidence of low-temperature water–rock interactions and various biogeochemical transformations these waters have undergone along their flow path.     

Redox-controlled iron isotope fractionation during magmatic differentiation: an example from the Red Hill intrusion, S. Tasmania

This study presents accurate and precise iron isotopic data for 16 co-magmatic rocks and 6 pyroxene–magnetite pairs from the classic, tholeiitic Red Hill sill in southern Tasmania. The intrusion exhibits a vertical continuum of compositions created by in situ fractional crystallisation of a single injection of magma in a closed igneous system and, as such, constitutes a natural laboratory amenable to determining the causes of Fe isotope fractionation in magmatic rocks. Early fractionation of pyroxenes and plagioclase, under conditions closed to oxygen exchange, gives rise to an iron enrichment trend and an increase in $ f_{{{\text{O}}_{2} }} $ of the melt relative to the Fayalite–Magnetite–Quartz (FMQ) buffer. Enrichment in Fe³⁺/ΣFe_melt is mirrored by δ⁵⁷Fe, where ^VIFe²⁺-bearing pyroxenes partition ⁵⁷Fe-depleted iron, defining an equilibrium pyroxene-melt fractionation factor of $ \Updelta^{57} {\text{Fe}}_{{{\text{px}} - {\text{melt}}}} \le - 0.25\,\permille \times 10^{6} /T^{2} $ . Upon magnetite saturation, the $ f_{{{\text{O}}_{2} }} $ and δ⁵⁷Fe of the melt fall, commensurate with the sequestration of the oxidised, ⁵⁷Fe-enriched iron into magnetite, quantified as $ \Updelta^{57} {\text{Fe}}_{{{\text{mtn}} - {\text{melt}}}} = + 0.20\,\permille \times 10^{6} /T^{2} $ . Pyroxene–magnetite pairs reveal an equilibrium fractionation factor of $ \Updelta^{57} {\text{Fe}}_{{{\text{mtn}} - {\text{px}}}} \approx + 0.30\,\permille $ at 900–1,000?°C. Iron isotopes in differentiated magmas suggest that they may act as an indicator of their oxidation state and tectonic setting.     

P,T, f_{{\text{CO}}_{\text{2}} } and f_{{\text{H}}_{\text{2}} {\text{O}}} during metamorphism of calcareous sediments in the Waterville-Vassalboro area,south-central Maine

In a regional metamorphic terrain where six isograds have been mapped based on mineral reactions that are observed in metacarbonate rocks, the P-T conditions and fugacities of CO₂ and H₂O during metamorphism were quantified by calculations involving actual mineral compositions and experimental data. Pressure during metamorphism was near 3,500 bars. Metamorphic temperatures ranged from 380° C (biotite-chlorite isograd) to 520° C (diopside isograd). \(f_{{\text{CO}}_{\text{2}} }\) and \(f_{{\text{CO}}_{\text{2}} }\) / \(f_{{\text{H}}_{\text{2}} {\text{O}}}\) in general is higher in metacarbonate rocks below the zoisite isograd than in those above the zoisite isograd. Calculated \(f_{{\text{CO}}_{\text{2}} }\) and \(f_{{\text{H}}_{\text{2}} {\text{O}}}\) are consistent with carbonate rocks above the zoisite isograd having equilibrated during metamorphism with a bulk supercritical fluid in which \(P_{{\text{CO}}_{\text{2}} }\) + \(P_{{\text{H}}_{\text{2}} {\text{O}}}\) = P _total. Calculations indicate that below the zoisite isograd, however, \(P_{{\text{CO}}_{\text{2}} }\) + \(P_{{\text{H}}_{\text{2}} {\text{O}}}\) was less than P_total, and that this condition is not due to the presence of significant amounts of species other than CO₂ and H₂O in the system C-O-H-S. Calculated \(P_{{\text{CO}}_{\text{2}} }\) /( \(P_{{\text{CO}}_{\text{2}} }\) + \(P_{{\text{H}}_{\text{2}} {\text{O}}}\) ) is low (0.06–0.32) above the zoisite isograd. The differences in conditions above and below the zoisite isograd may indicate that the formation of zoisite records the introduction of a bulk supercritical H₂O-rich fluid into the metacarbonates. The results of the study indicate that \(f_{{\text{CO}}_{\text{2}} }\) and \(f_{{\text{H}}_{\text{2}} {\text{O}}}\) are constant on a thin section scale, but that gradients in \(f_{{\text{CO}}_{\text{2}} }\) and \(f_{{\text{H}}_{\text{2}} {\text{O}}}\) existed during metamorphism on both outcrop and regional scales.     

Oxygen isotope fractionation between rutile and water

Synthetic rutile-water fractionations (1000 ln α) at 775, 675, and 575° C were found to be ?2.8, ?3.5, and ?4.8, respectively. Partial exchange experiments with natural rutile at 575° C and with synthetic rutile at 475° C failed to yield reliable fractionations. Isotopic fractionation within the range 575–775° C may be expressed as follows: 1 $$1000\ln \alpha ({\rm T}i{\rm O}_{2 } - H_2 O) = - 4.1 \frac{{10^6 }}{{T_{k^2 } }} + 0.96$$ . Combined with previously determined quartz-water fractionations, the above data permit calibration of the quartz-rutile geothermometer: 1 $$1000\ln \alpha ({\text{S}}i{\rm O}_{2 } - Ti{\rm O}_{2 } ) = 6.6 \frac{{10^6 }}{{T_{k^2 } }} - 2.9$$ . When applied to B-type eclogites from Europe, as an example, the latter equation yields a mean equilibration temperature of 565° C.     

(Ca,Mg)2Si2O6 clinopyroxenes: A solution model based on nonconvergent site-disorder

Experiments at high pressure and temperature indicate that excess Ca may be dissolved in diopside. If the (Ca, Mg)₂Si₂O₆ clinopyroxene solution extends to more Ca-rich compositions than CaMgSi₂O₆, macroscopic regular solution models cannot strictly be applied to this system. A nonconvergent site-disorder model, such as that proposed by Thompson (1969, 1970), may be more appropriate. We have modified Thompson's model to include asymmetric excess parameters and have used a linear least-squares technique to fit the available experimental data for Ca-Mg orthopyroxene-clinopyroxene equilibria and Fe-free pigeonite stability to this model. The model expressions for equilibrium conditions \(\mu _{{\text{Mg}}_{\text{2}} {\text{Si}}_{\text{2}} {\text{O}}_{\text{6}} }^{{\text{opx}}} = \mu _{{\text{Mg}}_{\text{2}} {\text{Si}}_{\text{2}} {\text{O}}_{\text{6}} }^{{\text{cpx}}} \) (reaction A) and \(\mu _{{\text{Ca}}_{\text{2}} {\text{Si}}_{\text{2}} {\text{O}}_{\text{6}} }^{{\text{opx}}} = \mu _{{\text{Ca}}_{\text{2}} {\text{Si}}_{\text{2}} {\text{O}}_{\text{6}} }^{{\text{cpx}}} \) (reaction B) are given by: 1 $$\begin{gathered} \Delta \mu _{\text{A}}^{\text{O}} = {\text{RT 1n}}\left[ {\frac{{(X_{{\text{Mg}}}^{{\text{opx}}} )^2 }}{{X_{{\text{Mg}}}^{{\text{M1}}} \cdot X_{{\text{Mg}}}^{{\text{M2}}} }}} \right] - \frac{1}{2}\{ W_{21} [2(X_{{\text{Ca}}}^{{\text{M2}}} )^3 - (X_{{\text{Ca}}}^{{\text{M2}}} ] \hfill \\ {\text{ + 2W}}_{{\text{22}}} [X_{{\text{Ca}}}^{{\text{M2}}} )^2 - (X_{{\text{Ca}}}^{{\text{M2}}} )^3 + \Delta {\text{G}}_{\text{*}}^{\text{0}} (X_{{\text{Ca}}}^{{\text{M1}}} \cdot X_{{\text{Ca}}}^{{\text{M2}}} )\} \hfill \\ {\text{ + W}}^{{\text{opx}}} (X_{{\text{Wo}}}^{{\text{opx}}} )^2 \hfill \\ \Delta \mu _{\text{B}}^{\text{O}} = {\text{RT 1n}}\left[ {\frac{{(X_{{\text{Ca}}}^{{\text{opx}}} )^2 }}{{X_{{\text{Ca}}}^{{\text{M1}}} \cdot X_{{\text{Ca}}}^{{\text{M2}}} }}} \right] - \frac{1}{2}\{ 2W_{21} [2(X_{{\text{Mg}}}^{{\text{M2}}} )^2 - (X_{{\text{Mg}}}^{{\text{M2}}} )^3 ] \hfill \\ {\text{ + W}}_{{\text{22}}} [2(X_{{\text{Mg}}}^{{\text{M2}}} )^3 - (X_{{\text{Mg}}}^{{\text{M2}}} )^2 + \Delta {\text{G}}_{\text{*}}^{\text{0}} (X_{{\text{Mg}}}^{{\text{M1}}} \cdot X_{{\text{Mg}}}^{{\text{M2}}} )\} \hfill \\ {\text{ + W}}^{{\text{opx}}} (X_{{\text{En}}}^{{\text{opx}}} )^2 \hfill \\ \hfill \\ \end{gathered} $$ where 1 $$\begin{gathered} \Delta \mu _{\text{A}}^{\text{O}} = 2.953 + 0.0602{\text{P}} - 0.00179{\text{T}} \hfill \\ \Delta \mu _{\text{B}}^{\text{O}} = 24.64 + 0.958{\text{P}} - (0.0286){\text{T}} \hfill \\ {\text{W}}_{{\text{21}}} = 47.12 + 0.273{\text{P}} \hfill \\ {\text{W}}_{{\text{22}}} = 66.11 + ( - 0.249){\text{P}} \hfill \\ {\text{W}}^{{\text{opx}}} = 40 \hfill \\ \Delta {\text{G}}_*^0 = 155{\text{ (all values are in kJ/gfw)}}{\text{.}} \hfill \\ \end{gathered} $$ . Site occupancies in clinopyroxene were determined from the internal equilibrium condition 1 $$\begin{gathered} \Delta G_{\text{E}}^{\text{O}} = - {\text{RT 1n}}\left[ {\frac{{X_{{\text{Ca}}}^{{\text{M1}}} \cdot X_{{\text{Mg}}}^{{\text{M2}}} }}{{X_{{\text{Ca}}}^{{\text{M2}}} \cdot X_{{\text{Mg}}}^{{\text{M1}}} }}} \right] + \tfrac{1}{2}[(2{\text{W}}_{{\text{21}}} - {\text{W}}_{{\text{22}}} )(2{\text{X}}_{{\text{Ca}}}^{{\text{M2}}} - 1) \hfill \\ {\text{ + }}\Delta G_*^0 (X_{{\text{Ca}}}^{{\text{M1}}} - X_{{\text{Ca}}}^{{\text{M2}}} ) + \tfrac{3}{2}(2{\text{W}}_{{\text{21}}} - {\text{W}}_{{\text{22}}} ) \hfill \\ {\text{ (1}} - 2X_{{\text{Ca}}}^{{\text{M1}}} )(X_{{\text{Ca}}}^{{\text{M1}}} + \tfrac{1}{2})] \hfill \\ \end{gathered} $$ where δG _E ⁰ =153+0.023T+1.2P. The predicted concentrations of Ca on the clinopyroxene Ml site are low enough to be compatible with crystallographic studies. Temperatures calculated from the model for coexisting ortho- and clinopyroxene pairs fit the experimental data to within 10° in most cases; the worst discrepancy is 30°. Phase relations for clinopyroxene, orthopyroxene and pigeonite are successfully described by this model at temperatures up to 1,600° C and pressures from 0.001 to 40 kbar. Predicted enthalpies of solution agree well with the calorimetric measurements of Newton et al. (1979). The nonconvergent site disorder model affords good approximations to both the free energy and enthalpy of clinopyroxenes, and, therefore, the configurational entropy as well. This approach may provide an example for Febearing pyroxenes in which cation site exchange has an even more profound effect on the thermodynamic properties.     

The Hydrolysis of Al(III) in NaCl solutions: A Model for M(II), M(III), and M(IV) Ions

Understanding the identity and stability of the hydrolysis products of metals is required in order to predict their behavior in natural aquatic systems. Despite this need, the hydrolysis constants of many metals are only known over a limited range of temperature and ionic strengths. In this paper, we show that the hydrolysis constants of 31 metals [i.e. Mn(II), Cr(III), U(IV), Pu(IV)] are nearly linearly related to the values for Al(III) over a wide range of temperatures and ionic strengths. These linear correlations allow one to make reasonable estimates for the hydrolysis constants of +2, +3, and +4 metals from 0 to 300°C in dilute solutions and 0 to 100°C to 5 m in NaCl solutions. These correlations in pure water are related to the differences between the free energies of the free ion and complexes being almost equal $$ \Updelta {\text{G}}^\circ \left( {{\text{Al}}^{3 + } } \right) - \Updelta {\text{G}}^\circ \left( {{\text{Al}}\left( {\text{OH}} \right)_{j}^{{\left( {3 - j} \right)}} } \right) \cong \Updelta {\text{G}}^\circ \left( {{\text{M}}^{n + } } \right) - \Updelta {\text{G}}^\circ \left( {{\text{M}}\left( {\text{OH}} \right)_{j}^{{\left( {n - j} \right)}} } \right) $$ The correlation at higher temperatures is a result of a similar relationship between the enthalpies of the free ions and complexes $$ \Updelta {\text{H}}^\circ \left( {{\text{Al}}^{3 + } } \right) - \Updelta {\text{H}}^\circ \left( {{\text{Al}}\left( {\text{OH}} \right)_{j}^{3 - j} } \right) \cong \Updelta {\text{H}}^\circ \left( {{\text{M}}^{n + } } \right) - \Updelta {\text{H}}^\circ \left( {{\text{M}}\left( {\text{OH}} \right)_{j}^{n - j} } \right) $$ The correlations at higher ionic strengths are the result of the ratio of the activity coefficients for Al(III) being almost equal to that of the metal. $$ \gamma \left( {{\text{M}}^{n + } } \right)/\gamma \left( {{\text{M}}\left( {\text{OH}} \right)_{j}^{n - j} } \right) \cong \gamma \left( {{\text{Al}}^{3 + } } \right)/\gamma \left( {{\text{Al}}\left( {\text{OH}} \right)_{j}^{3 - j} } \right) $$ The results of this study should be useful in examining the speciation of metals as a function of pH in natural waters (e.g. hydrothermal fresh waters and NaCl brines).     

The system Fe-Si-O: Oxygen buffer calibrations to 1,500K

The five solid-phase oxygen buffers of the system Fe-Si-O, iron-wuestite (IW), wuestite-magnetite (WM), magnetite-hematite (MH), quartz-iron-fayalite (QIF) and fayalite-magnetite-quartz (FMQ) have been recalibrated at 1 atm pressure and temperatures from 800°–1,300° C, using a thermogravimetric gas mixing furnace. The oxygen fugacity, \(f_{{\text{O}}_{\text{2}} }\) was measured with a CaO-doped ZrO₂ electrode. Measurements were made also for wuestite solid solutions in order to determine the redox behavior of wuestites with O/Fe ratios varying from 1.05 to 1.17. For FMQ, additional determinations were carried out at 1 kb over a temperature range of 600° to 800° C, using a modified Shaw membrane. Results agree reasonably well with published data and extrapolations. The reaction parameters K, ΔG _r ^o , ΔH _r ^o , and ΔS _r ^o were calculated from the following log \(f_{{\text{O}}_{\text{2}} }\) /T relations (T in K): $$\begin{gathered} {\text{IW }}\log f_{{\text{O}}_{\text{2}} } = - 26,834.7/T + 6.471\left( { \pm 0.058} \right) \hfill \\ {\text{ }}\left( {{\text{800}} - 1,260{\text{ C}}} \right), \hfill \\ {\text{WM }}\log f_{{\text{O}}_{\text{2}} } = - 36,951.3/T + 16.092\left( { \pm 0.045} \right) \hfill \\ {\text{ }}\left( {{\text{1,000}} - 1,300{\text{ C}}} \right), \hfill \\ {\text{MH }}\log f_{{\text{O}}_{\text{2}} } = - 23,847.6/T + 13.480\left( { \pm 0.055} \right) \hfill \\ {\text{ }}\left( {{\text{1,040}} - 1,270{\text{ C}}} \right), \hfill \\ {\text{QIF }}\log f_{{\text{O}}_{\text{2}} } = - 27,517.5/T + 6.396\left( { \pm 0.049} \right) \hfill \\ {\text{ }}\left( {{\text{960}} - 1,140{\text{ C}}} \right), \hfill \\ {\text{FMQ }}\log f_{{\text{O}}_{\text{2}} } = - 24,441.9/T + 8.290\left( { \pm 0.167} \right) \hfill \\ {\text{ }}\left( {{\text{600}} - 1,140{\text{ C}}} \right). \hfill \\ \end{gathered}$$ These experimentally determined reaction parameters were combined with published 298 K data to determine the parameters G_f, H_f, and S_f for the phases wuestite, magnetite, hematite, and fayalite from 298 K to the temperatures of the experiments. The T? \(f_{{\text{O}}_{\text{2}} }\) data for wuestite solid solutions were used to obtain activities, excess free energies and Margules mixing parameters. The new data provide a more reliable, consistent and complete reference set for the interpretation of redox reactions at elevated temperatures in experiments and field settings encompassing the crust, mantle and core as well as extraterrestrial environments.     

Experimental investigation and application of the equilibrium rutile + orthopyroxene = quartz + ilmenite

Equilibria in the Sirf (Silica-Ilmenite-Rutile-Ferrosilite) system: $${\text{SiO}}_{\text{2}} + ({\text{Mg,Fe}}){\text{TiO}}_{\text{3}} {\text{ + (Mg,Fe)SiO}}_{\text{3}} $$ have been calibrated in the range 800–1100° C and 12–26 kbar using a piston-cylinder apparatus to assess the potential of the equilibria for geobarometry in granulite facies assemblages that lack garnet. Thermodynamic calculations indicate that the two end-member equilibria involving quartz + geikielite = rutile + enstatite, and quartz + ilmenite = rutile + ferrosilite, are metastable. We therefore reversed equilibria over the compositional range Fs_40–70, using Ag₈₀Pd₂₀ capsules with \(f_{{\text{O}}_{\text{2}} } \) buffered at or near iron-wüstite. Ilmenite compositions coexisting with orthopyroxene are \(X_{{\text{MgTiO}}_{\text{3}} }^{{\text{Ilm}}} \) of 0.06 to 0.15 and \(X_{{\text{Fe}}_{\text{2}} {\text{O}}_{\text{3}} }^{{\text{Ilm}}} \) of 0.00 to 0.01, corresponding toK _D values of 13.3, 10.2, 9.0 and 8.0 (±0.5) at 800, 900, 1000 and 1100° C, respectively, whereK _D =(XMg/XFe)^Opx/(XMg/XFe)^Ilm. Pressures have been calculated using equilibria in the Sirf system for granulites from the Grenville Province of Ontario and for granulite facies xenoliths from central Mexico. Pressures are consistent with other well-calibrated geobarometers for orthopyroxeneilmenite pairs from two Mexican samples in which oxide textures appear to represent equilibrium. Geologically unreasonable pressures are obtained, however, where oxide textures are complex. Application of data from this study on the equilibrium distribution of iron and magnesium between ilmenite and orthopyroxene suggests that some ilmenite in deep crustal xenoliths is not equilibrated with coexisting pyroxene, while assemblages from exposed granulite terranes have reequilibrated during retrogression. The Sirf equilibria are sensitive to small changes in composition and may be used for determination of activity/composition (a/X) relations of orthopyroxene if an ilmenite model is specified. A symmetric regular solution model has been used for orthopyroxene in conjunction with activity models for ilmenite available from the literature to calculatea/X relations in orthopyroxene of intermediate composition. Data from this study indicate that FeSiO₃?MgSiO₃ orthopyroxene exhibits small, positive deviations from ideality over the range 800–1100°C.     

A note on the significance of the assemblage calcite-quartz-plagioclase-paragonite-graphite

The biotite zone assemblage: calcite-quartz-plagioclase (An₂₅)-phengite-paragonite-chlorite-graphite, is developed at the contact between a carbonate and a pelite from British Columbia. Thermochemical data for the equilibrium paragonite+calcite+2 quartz=albite+ anorthite+CO₂+H₂O yields: $$\log f{\text{H}}_{\text{2}} {\text{O}} + \log f{\text{CO}}_{\text{2}} = 5.76 + 0.117 \times 10^{ - 3} (P - 1)$$ for a temperature of 700°K and a plagioclase composition of An₂₅. By combining this equation with equations describing equilibria between graphite and gas species in the system C-H-O, the following partial pressures: \(P{\text{H}}_2 {\text{O}} = 2572{\text{b, }}P{\text{CO}}_2 = 3162{\text{b, }}P{\text{H}}_2 = 2.5{\text{b, }}P{\text{CH}}_4 = 52.5{\text{b, }}P{\text{CO}} = 11.0{\text{b}}\) are obtained for \(f{\text{O}}_2 = 10^{ - 26}\) . If total pressure equals fluid pressure, then the total pressure during metamorphism was approximately 6 kb. The total fluid pressure calculated is extremely sensitive to the value of \(f{\text{O}}_2\) chosen.     

Experimental determination of the reaction paragonite=jadeite+kyanite+H2O,and internally consistent thermodynamic data for part of the system Na2O−Al2O3−SiO2−H2O,with applications to eclogites and blueschists

Hydrothermal reversal experiments have been performed on the upper pressure stability of paragonite in the temperature range 550–740 ° C. The reaction $$\begin{gathered} {\text{NaAl}}_{\text{3}} {\text{Si}}_{\text{3}} {\text{O}}_{{\text{1 0}}} ({\text{OH)}}_{\text{2}} \hfill \\ {\text{ paragonite}} \hfill \\ {\text{ = NaAlSi}}_{\text{2}} {\text{O}}_{\text{6}} + {\text{Al}}_{\text{2}} {\text{SiO}}_{\text{5}} + {\text{H}}_{\text{2}} {\text{O}} \hfill \\ {\text{ jadeite kyanite vapour}} \hfill \\ \end{gathered}$$ has been bracketed at 550 ° C, 600 ° C, 650 ° C, and 700 ° C, at pressures 24–26 kb, 24–25.5 kb, 24–25 kb, and 23–24.5 kb respectively. The reaction has a shallow negative slope (? 10 bar °C^?1) and is of geobarometric significance to the stability of the eclogite assemblage, omphacite+kyanite. The experimental brackets are thermodynamically consistent with the lower pressure reversals of Chatterjee (1970, 1972), and a set of thermodynamic data is presented which satisfies all the reversal brackets for six reactions in the system Na₂O-Al₂O₃-SiO₂-H₂O. The Modified Redlich Kwong equation for H₂O (Holloway, 1977) predicts fugacities which are too high to satisfy the reversals of this study. The P-T stabilities of important eclogite and blueschist assemblages involving omphacite, kyanite, lawsonite, Jadeite, albite, chloritoid, and almandine with paragonite have been calculated using thermodynamic data derived from this study.     

Talc paragenesis in some siliceous dolomitic rocks,and its sedimentologic significance

Thin (0.5–4 mm), contorted stringers of talc, associated with apatite and minor pyrite, are containdy Formation in eastern Alabama. The form, position and lithologic distribed within generally saccharoidal dolomite-quartz marbles of the Cambrian Shaution of the stringers strongly suggest an algalstromatolitic origin, with interlaminar trapped dolomitic muds. Metamorphic formation of talc plus apatite proceeded only within the stringers, whereas surrounding marble remained as unreacted dolomite plus quartz. Talc generation is best explained by the reaction $${\text{dolomite}} + {\text{silica}} + {\text{water}} + {\text{P}}_{\text{2}} {\text{O}}_{\text{5}} = {\text{talc}} + {\text{apatite}} + {\text{CO}}_{\text{2}}$$ in which the phosphate was supplied to the reaction from organic matter contained within the stromatolitic layers. The system was probably open to CO₂ during metamorphism, and \(P_{{\text{CO}}_{\text{2}} }\) remained relatively low.     

The reversed alumina contents of orthopyroxene in equilibrium with spinel and forsterite in the system MgO-Al2O3-SiO2

Equilibrium alumina contents of orthopyroxene coexisting with spinel and forsterite in the system MgO-Al₂O₃-SiO₂ have been reversed at 15 different P-T conditions, in the range 1,030–1,600° C and 10–28 kbar. The present data and three reversals of Danckwerth and Newton (1978) have been modeled assuming an ideal pyroxene solid solution with components Mg₂Si₂O₆ (En) and MgAl₂SiO₆ (MgTs), to yield the following equilibrium condition (J, bar, K): $$\begin{gathered} RT{\text{ln(}}X_{{\text{MgTs}}} {\text{/}}X_{{\text{En}}} {\text{) + 29,190}} - {\text{13}}{\text{.42 }}T + 0.18{\text{ }}T + 0.18{\text{ }}T^{1.5} \hfill \\ + \int\limits_1^P {\Delta V_{T,P}^{\text{0}} dP = 0,} \hfill \\ \end{gathered} $$ where $$\begin{gathered} + \int\limits_1^P {\Delta V_{T,P}^{\text{0}} dP} \hfill \\ = [0.013 + 3.34 \times 10^{ - 5} (T - 298) - 6.6 \times 10^{ - 7} P]P. \hfill \\ \end{gathered} $$ The data of Perkins et al. (1981) for the equilibrium of orthopyroxene with pyrope have been similarly fitted with the result: $$\begin{gathered} - RT{\text{ln(}}X_{{\text{MgTs}}} \cdot X_{{\text{En}}} {\text{) + 5,510}} - 88.91{\text{ }}T + 19{\text{ }}T^{1.2} \hfill \\ + \int\limits_1^P {\Delta V_{T,P}^{\text{0}} dP = 0,} \hfill \\ \end{gathered} $$ where $$\begin{gathered} + \int\limits_1^P {\Delta V_{T,P}^{\text{0}} dP} \hfill \\ = [ - 0.832 - 8.78{\text{ }} \times {\text{ 10}}^{ - {\text{5}}} (T - 298) + 16.6{\text{ }} \times {\text{ 10}}^{ - 7} P]{\text{ }}P. \hfill \\ \end{gathered} $$ The new parameters are in excellent agreement with measured thermochemical data and give the following properties of the Mg-Tschermak endmember: $$H_{f,970}^0 = - 4.77{\text{ kJ/mol, }}S_{298}^0 = 129.44{\text{ J/mol}} \cdot {\text{K,}}$$ and $$V_{298,1}^0 = 58.88{\text{ cm}}^{\text{3}} .$$ The assemblage orthopyroxene+spinel+olivine can be used as a geothermometer for spinel lherzolites, subject to a choice of thermodynamic mixing models for multicomponent orthopyroxene and spinel. An ideal two-site mixing model for pyroxene and Sack's (1982) expressions for spinel activities provide, with the present experimental calibration, a geothermometer which yields temperatures of 800° C to 1,350° C for various alpine peridotites and 850° C to 1,130° C for various volcanic inclusions of upper mantle origin.     

Experimental investigation of the reaction forsterite + H2O ⇌ serpentine + brucite

A new determination of the equilibrium reaction: $$\begin{gathered} 2{\text{ Mg}}_{\text{2}} [{\text{SiO}}_{\text{4}} ] + 3{\text{ H}}_{\text{2}} {\text{O}} \rightleftharpoons {\text{1 Mg}}_{\text{3}} [({\text{OH)}}_{\text{4}} |{\text{Si}}_{\text{2}} {\text{O}}_{\text{5}} ] + 1{\text{ Mg(OH)}}_{\text{2}} \hfill \\ \hfill \\ {\text{ forsterite serpentine brucite}} \hfill \\ \end{gathered} $$ yielded equilibrium temperatures which lie (at identical H₂O-pressures) about 60° C lower than all previously published data (Bowen and Tuttle, 1949; Yoder, 1952; Kitahara et al., 1966; Kitahara and Kennedy, 1967). It has been shown that the above authors have determined not the stable equilibrium curve but instead a metastable “synthesis boundary”. The actual (stable) equilibrium curve is located at 0,5 kb and 350° C 2,0 kb and 380° C 3,5 kb and 400° C 5,0 kb and 420° C 6,5 kb and 430° C.     

The role of the minor substitutions in the crystal structure of natural challacolloite, KPb2Cl5, and hephaistosite, TlPb2Cl5, from Vulcano (Aeolian Archipelago, Italy)

Crystals of challacolloite, KPb₂Cl₅, and hephaistosite, TlPb₂Cl₅, from volcanic sublimates formed on the crater rim of the “La Fossa Crater” at Vulcano, Aeolian Archipelago, Italy, were investigated. Chemical compositions were ${\left( {{\text{K}}_{{0.93}} {\text{Tl}}_{{0.02}} } \right)}_{{\Sigma = 0.95}} {\text{Pb}}_{{2.04}} {\left( {{\text{Cl}}_{{4.90}} {\text{Br}}_{{0.11}} } \right)}_{{\Sigma = 5.01}} $ and ${\text{Tl}}_{{0.94}} {\text{Pb}}_{{2.01}} {\left( {{\text{Cl}}_{{4.91}} {\text{Br}}_{{0.14}} } \right)}_{{\Sigma = 5.05}} $ , respectively. Single crystal X-ray measurements showed monoclinic symmetry for both phases, space group P2₁/c, with the following unit-cell parameters: a = 8.8989(4), b = 7.9717(5), c = 12.5624(8) Å, β = 90.022(4)°, V = 891.2(1) ų, Z = 4 (challacolloite) and a = 9.0026(6), b = 7.9723(6), c = 12.5693(9) Å, β = 90.046(4)°, V = 902.1(1) ų, Z = 4 (hephaistosite). The structure refinements converge to R = 3.99% and R = 3.86%, respectively. The effects of Br?Cl and K?Tl substitutions on the structure of these natural compounds have been discussed.     

Djerfisherite: nebular source of refractory potassium

Djerfisherite is an important carrier of potassium in highly reduced enstatite chondrites, where it occurs in sub-round metal-sulfide nodules. These nodules were once free-floating objects in the protoplanetary nebula. Here, we analyze existing and new data to derive an equation of state (EOS) for djerfisherites of $ {\text{K}}_{ 6} ({\text{Cu}},{\text{Fe}},{\text{Ni}})^{B} ({\text{Fe}},{\text{Ni}},{\text{Cu}})_{24}^{C} {\text{S}}_{ 2 6} {\text{Cl}} $ K 6 ( Cu , Fe , Ni ) B ( Fe , Ni , Cu ) 24 C S 2 6 Cl structural formula. We use this EOS to calculate the thermal stability of djerfisherite coexisting in equilibrium with a cooling vapor of solar composition enriched in a dust analogous to anhydrous, chondritic interplanetary dust (C-IDP). We find that condensed mineral assemblages closely match those found in enstatite chondrites, with djerfisherite condensing above 1,000 K in C-IDP dust-enriched systems. Results may have implications for the volatile budgets of terrestrial planets and the incorporation of K into early formed, highly reduced, planetary cores. Previous work links enstatite chondrites to the planet Mercury, where the surface has a terrestrial K/Th ratio, high S/Si ratio, and very low FeO content. Mercury’s accretion history may yield insights into Earth’s.     

Gehlenite stability in the system CaO-Al2O3-SiO2-H2O-CO2

P, T, \(X_{{\text{CO}}_{\text{2}} }\) relations of gehlenite, anorthite, grossularite, wollastonite, corundum and calcite have been determined experimentally at P _f =1 and 4 kb. Using synthetic starting minerals the following reactions have been demonstrated reversibly

1. 2 anorthite+3 calcite=gehlenite+grossularite+3 CO₂.
2. anorthite+corundum+3 calcite=2 gehlenite+3 CO₂.
3. 3anorthite+3 calcite=2 grossularite+corundum+3CO₂.
4. grossularite+2 corundum+3 calcite=3 gehlenite+3 CO₂.
5. anorthite+2 calcite=gehlenite+wollastonite+2CO₂.
6. anorthite+wollastonite+calcite=grossularite+CO₂.
7. grossularite+calcite=gehlenite+2 wollastonite+CO₂.

In the T, \(X_{{\text{CO}}_{\text{2}} }\) diagram at P _f =1 kb two isobaric invariant points have been located at 770±10°C, \(X_{{\text{CO}}_{\text{2}} }\) =0.27 and at 840±10°C, \(X_{{\text{CO}}_{\text{2}} }\) =0.55. Formation of gehlenite from low temperature assemblages according to (4) and (2) takes place at 1 kb and 715–855° C, \(X_{{\text{CO}}_{\text{2}} }\) =0.1–1.0. In agreement with experimental results the formation of gehlenite in natural metamorphic rocks is restricted to shallow, high temperature contact aureoles.     

Crystal chemistry of Fe3+-bearing (Mg,Fe)SiO3 perovskite: a single-crystal X-ray diffraction study

Magnesium silicate perovskite is the predominant phase in the Earth’s lower mantle, and it is well known that incorporation of iron has a strong effect on its crystal structure and physical properties. To constrain the crystal chemistry of (Mg, Fe)SiO₃ perovskite more accurately, we synthesized single crystals of Mg_0.946(17)Fe_0.056(12)Si_0.997(16)O₃ perovskite at 26 GPa and 2,073 K using a multianvil press and investigated its crystal structure, oxidation state and iron-site occupancy using single-crystal X-ray diffraction and energy-domain Synchrotron Mössbauer Source spectroscopy. Single-crystal refinements indicate that all iron (Fe²⁺ and Fe³⁺) substitutes on the A-site only, where \( {\text{Fe}}^{ 3+ } /\Upsigma {\text{Fe}}\sim 20\,\% \) based on Mössbauer spectroscopy. Charge balance likely occurs through a small number of cation vacancies on either the A- or the B-site. The octahedral tilt angle (Φ) calculated for our sample from the refined atomic coordinates is 20.3°, which is 2° higher than the value calculated from the unit-cell parameters (a = 4.7877 Å, b = 4.9480 Å, c = 6.915 Å) which assumes undistorted octahedra. A compilation of all available single-crystal data (atomic coordinates) for (Mg, Fe)(Si, Al)O₃ perovskite from the literature shows a smooth increase of Φ with composition that is independent of the nature of cation substitution (e.g., \( {\text{Mg}}^{ 2+ } - {\text{Fe}}^{ 2+ } \) or \( {\text{Mg}}^{ 2+ } {\text{Si}}^{ 4+ } - {\text{Fe}}^{ 3+ } {\text{Al}}^{ 3+ } \) substitution mechanism), contrary to previous observations based on unit-cell parameter calculations.     

Partition of Ni between olivine and sulfide: the effect of temperature,f_{{\text{O}}_{\text{2}} } and f_{{\text{S}}_{\text{2}} }

The experimental distribution coefficient for Ni/ Fe exchange between olivine and monosulfide (K_D3) is 35.6±1.1 at 1385° C, \(f_{{\text{O}}_{\text{2}} } = 10^{ - 8.87} ,f_{{\text{S}}_{\text{2}} } = 10^{ - 1.02} \) , and olivine of composition Fo96 to Fo92. These are the physicochemical conditions appropriate to hypothesized sulfur-saturated komatiite magma. The present experiments equilibrated natural olivine grains with sulfide-oxide liquid in the presence of a (Mg, Fe)-alumino-silicate melt. By a variety of different experimental procedures, K _D3 is shown to be essentially constant at about 30 to 35 in the temperature range 900 to 1400° C, for olivine of composition Fo97 to FoO, monosulfide composition with up to 70 mol. % NiS, and a wide range of \(f_{{\text{O}}_{\text{2}} } \) and \(f_{{\text{S}}_{\text{2}} } \) .