Equations of state of Plagioclase Feldspars

The volume variation with pressure of seven intermediate plagioclase feldspars has been determined by high-pressure single-crystal X-ray diffraction. The bulk moduli of plagioclases for a 3rd-order Birch-Murnaghan EoS can be described by the following pair of equations:

The magnesiowüstite: iron equilibrium and its implications for the activity-composition relations of (Mg,Fe)<Subscript>2</Subscript>SiO<Subscript>4</Subscript> olivine solid solutions

The chemical potential of oxygen (µO₂) in equilibrium with magnesiowüstite solid solution (Mg, Fe)O and metallic Fe has been determined by gas-mixing experiments at 1,473 K supplemented by solid-cell EMF experiments at lower temperatures. The results give:

Experimental study of a low-alkali tholeiite at 1–5 kbar: optimal condition for the crystallization of high-An plagioclase in hydrous arc tholeiite

We conducted melting experiments on a low-alkali tholeiite (SiO₂ ~52 wt%, MgO ~6.5 wt%, CaO/Na₂O~4.4, Al₂O₃/SiO₂ ~0.33) under both H₂O-undersaturated and H₂O-saturated conditions to investigate the effect of H₂O on the Ca–Na partitioning between plagioclase and melt. Experiments were performed in the temperature and pressure ranges of 1,000–1,300°C and 1–5 kbar, respectively, with varying H₂O contents of 0–12wt%. Redox condition was 0–2 log unit above NNO (nickel–nickel oxide) buffer. Temperature-bulk H₂O diagrams for the low-alkali tholeiite are constructed at 1, 2, and 5 kbar, and compositions of near-liquidus plagioclase and coexisting melt are determined. To exclude the effect of melt composition (CaO/Na₂O and Al₂O₃/SiO₂ ratios) on plagioclase composition and to reveal the effect of H₂O on An (=100×Ca/(Ca+Na)) content and (=(Ca/Na)_pl/(Ca/Na)_melt), we focused on the composition of near-liquidus plagioclases which crystallized from melts with nearly constant CaO/Na₂O and Al₂O₃/SiO₂ ratios. Our experimental results show that, at each experimental pressure, An content of the near-liquidus plagioclase and the K_D^Ca-Na almost linearly increases as H₂O content in melt increases. Each of the An content and the variations in a low-alkali tholeiitic system (CaO/Na₂O~4.0–4.5, Al₂O₃/SiO₂ ~0.27–0.33) can be described by one equation using temperature, pressure, and melt H₂O content as parameters. An content and of liquidus plagioclase increases with increasing melt H₂O and with decreasing pressure, elucidating that nearly H₂O-saturated conditions of 2–3 kbar is optimal for the crystallization of the most An-rich plagioclase (>An₈₈). We suggest this pressure condition of 2–3 kbar, corresponding to depth of 7–11 km, plays an important role for the origin of An-rich plagioclase in H₂O-rich low-alkali tholeiite. At pressures more than ca. 4 kbar, crystallization of liquidus Ca-rich clinopyroxene decreases the CaO/Na₂O ratio of liquid, thus prohibiting the crystallization of high-An plagioclase from hydrous tholeiite.     

Experimental and computational study of trace element distribution between orthopyroxene and anhydrous silicate melt: substitution mechanisms and the effect of iron

Although orthopyroxene (Opx) is present during a wide range of magmatic differentiation processes in the terrestrial and lunar mantle, its effect on melt trace element contents is not well quantified. We present results of a combined experimental and computational study of trace element partitioning between Opx and anhydrous silicate melts. Experiments were performed in air at atmospheric pressure and temperatures ranging from 1,326 to 1,420°C in the system CaO–MgO–Al₂O₃–SiO₂ and subsystem CaO–MgO–SiO₂. We provide experimental partition coefficients for a wide range of trace elements (large ion lithophile: Li, Be, B, K, Rb, Sr, Cs, Ba, Th, U; rare earth elements, REE: La, Ce, Nd, Sm, Y, Yb, Lu; high field strength: Zr, Nb, Hf, Ta, Ti; transition metals: Sc, V, Cr, Co) for use in petrogenetic modelling. REE partition coefficients increase from $ D_{\text{La}}^{{\text{Opx}} {\hbox{-}} {\text{melt}}} \sim 0.0005 Although orthopyroxene (Opx) is present during a wide range of magmatic differentiation processes in the terrestrial and lunar mantle, its effect on melt trace element contents is not well quantified. We present results of a combined experimental and computational study of trace element partitioning between Opx and anhydrous silicate melts. Experiments were performed in air at atmospheric pressure and temperatures ranging from 1,326 to 1,420°C in the system CaO–MgO–Al₂O₃–SiO₂ and subsystem CaO–MgO–SiO₂. We provide experimental partition coefficients for a wide range of trace elements (large ion lithophile: Li, Be, B, K, Rb, Sr, Cs, Ba, Th, U; rare earth elements, REE: La, Ce, Nd, Sm, Y, Yb, Lu; high field strength: Zr, Nb, Hf, Ta, Ti; transition metals: Sc, V, Cr, Co) for use in petrogenetic modelling. REE partition coefficients increase from $ D_{\text{La}}^{{\text{Opx}} {\hbox{-}} {\text{melt}}} \sim 0.0005 $ D_{\text{La}}^{{\text{Opx}} {\hbox{-}} {\text{melt}}} \sim 0.0005 to $ D_{\text{Lu}}^{{{\text{Opx}} {\hbox{-}} {\text{melt}}}} \sim 0.109 $ D_{\text{Lu}}^{{{\text{Opx}} {\hbox{-}} {\text{melt}}}} \sim 0.109 , D values for highly charged elements vary from $ D_{\text{Th}}^{{{\text{Opx}} {\hbox{-}} {\text{melt}}}} \sim 0.0026 $ D_{\text{Th}}^{{{\text{Opx}} {\hbox{-}} {\text{melt}}}} \sim 0.0026 through $ D_{\text{Nb}}^{{{\text{Opx}} {\hbox{-}} {\text{melt}}}} \sim 0.0033 $ D_{\text{Nb}}^{{{\text{Opx}} {\hbox{-}} {\text{melt}}}} \sim 0.0033 and $ D_{\text{U}}^{{{\text{Opx}} {\hbox{-}} {\text{melt}}}} \sim 0.0066 $ D_{\text{U}}^{{{\text{Opx}} {\hbox{-}} {\text{melt}}}} \sim 0.0066 to $ D_{\text{Ti}}^{{\text{Opx}} {\hbox{-}} {\text{melt}}} \sim 0.058 $ D_{\text{Ti}}^{{\text{Opx}} {\hbox{-}} {\text{melt}}} \sim 0.058 , and are all virtually independent of temperature. Cr and Co are the only compatible trace elements at the studied conditions. To elucidate charge-balancing mechanisms for incorporation of REE into Opx and to assess the possible influence of Fe on Opx-melt partitioning, we compare our experimental results with computer simulations. In these simulations, we examine major and minor trace element incorporation into the end-members enstatite (Mg₂Si₂O₆) and ferrosilite (Fe₂Si₂O₆). Calculated solution energies show that R²⁺ cations are more soluble in Opx than R³⁺ cations of similar size, consistent with experimental partitioning data. In addition, simulations show charge balancing of R³⁺ cations by coupled substitution with Li⁺ on the M1 site that is energetically favoured over coupled substitution involving Al–Si exchange on the tetrahedrally coordinated site. We derived best-fit values for ideal ionic radii r ₀, maximum partition coefficients D ₀, and apparent Young’s moduli E for substitutions onto the Opx M1 and M2 sites. Experimental r ₀ values for R³⁺ substitutions are 0.66–0.67 ? for M1 and 0.82–0.87 ? for M2. Simulations for enstatite result in r ₀ = 0.71–0.73 ? for M1 and ~0.79–0.87 ? for M2. Ferrosilite r ₀ values are systematically larger by ~0.05 ? for both M1 and M2. The latter is opposite to experimental literature data, which appear to show a slight decrease in $ r_{0}^{{{\text{M}}2}} $ r_{0}^{{{\text{M}}2}} in the presence of Fe. Additional systematic studies in Fe-bearing systems are required to resolve this inconsistency and to develop predictive Opx-melt partitioning models for use in terrestrial and lunar magmatic differentiation models.     

Fluorite stability in silicic magmas

Recent experimental evidence is used to assess the conditions under which fluorite forms an early crystallising phase in silicic magmas. Fluorite solubility primarily depends on the (Na + K)/Al balance in the coexisting silicic melt, reaching a minimum in metaluminous melts. It can display reaction relationships with topaz and titanite, depending on changes in melt composition during crystallisation. An empirical model of fluorite stability in Ca-poor peralkaline rhyolite melts is derived and applied to selected rocks:

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).     

H<Subscript>2</Subscript>O diffusion in peralkaline to peraluminous rhyolitic melts

The diffusion of water in a peralkaline and a peraluminous rhyolitic melt was investigated at temperatures of 714–1,493 K and pressures of 100 and 500 MPa. At temperatures below 923 K dehydration experiments were performed on glasses containing about 2 wt% H₂O _t in cold seal pressure vessels. At high temperatures diffusion couples of water-poor (<0.5 wt% H₂O _t ) and water-rich (~2 wt% H₂O _t ) melts were run in an internally heated gas pressure vessel. Argon was the pressure medium in both cases. Concentration profiles of hydrous species (OH groups and H₂O molecules) were measured along the diffusion direction using near-infrared (NIR) microspectroscopy. The bulk water diffusivity () was derived from profiles of total water () using a modified Boltzmann-Matano method as well as using fittings assuming a functional relationship between and Both methods consistently indicate that is proportional to in this range of water contents for both bulk compositions, in agreement with previous work on metaluminous rhyolite. The water diffusivity in the peraluminous melts agrees very well with data for metaluminous rhyolites implying that an excess of Al₂O₃ with respect to alkalis does not affect water diffusion. On the other hand, water diffusion is faster by roughly a factor of two in the peralkaline melt compared to the metaluminous melt. The following expression for the water diffusivity in the peralkaline rhyolite as a function of temperature and pressure was obtained by least-squares fitting:

Ca-Fe-Mg olivines: phase relations and a solution model

Reversed phase equilibrium experiments in the system (Ca, Mg, Fe)₂SiO₄ provide four tielines at P?1 bar and 1 kbar and 800° C–1,100° C. These tielines have been used to model the solution properties of the olivine quadrilateral following the methods described by Davidson et al. (1981) for quadrilateral clinopyroxenes. The discrepancy between the calculated phase relations and the experimentally determined tielines is within the uncertainty of the experiments. The solution properties of quadrilateral olivines can be described by a non-convergent site-disorder model that allows for complete partitioning of Ca on the M2 site, highly disordered Fe-Mg cation distributions and limited miscibility between high-Ca and low-Ca olivines. The ternary data presented in this paper together with binary solution models for the joins Fo-Mo and Fa-Kst have been used to evaluate two solution parameters: $$\begin{gathered} F^0 \equiv 2(\mu _{{\rm M}o}^0 - \mu _{{\rm K}st}^0 ) + \mu _{Fa}^0 - \mu _{Fo}^0 = 12.660 (1.6) kJ, \hfill \\ \Delta G_*^0 \equiv \mu _{{\rm M}gFe}^0 + \mu _{FeMg}^0 - \mu _{Fo}^0 - \mu _{Fa}^0 = 7.030 (3.9) kJ. \hfill \\ \end{gathered} $$ . Ternary phase quilibrium data for olivines tightly constrain the value of F⁰, but not that for ΔG _* ⁰ which describes nonideality in Fe-Mg mixing. From this analysis, we infer a function for the apparent standard state energy of Kst: $$\begin{gathered} \mu _{{\rm K}st}^0 = - 102.79 \pm 0.8 - (T - 298)(0.137026) \hfill \\ + (T - 298 - T1n(T/298))(0.155519) \hfill \\ + (T - 298)^2 (2.8242E - 05)/2 \hfill \\ + (T - 298)^2 (2.9665E + 03)/(2T(298)^2 ) kJ \hfill \\ \end{gathered} $$ where T is in Kelvins and the 298 K value is relative to oxides.     

A Coupled Riverine-Marine Fractionation Model for Dissolved Rare Earths and Yttrium

Fractionation of yttrium (Y) and the rare earth elements (REEs) begins in riverine systems and continues in estuaries and the ocean. Models of yttrium and rare earth (YREE) distributions in seawater must therefore consider the fractionation of these elements in both marine and riverine systems. In this work we develop a coupled riverine/marine fractionation model for dissolved rare earths and yttrium, and apply this model to calculations of marine YREE fractionation for a simple two-box (riverine/marine) geochemical system. Shale-normalized YREE concentrations in seawater can be expressed in terms of fractionation factors ( _ij ) appropriate to riverine environments ( ) and seawater ( ):

The effect of liquid composition on the partitioning of Ni between olivine and silicate melt

We report the results of experiments designed to separate the effects of temperature and pressure from liquid composition on the partitioning of Ni between olivine and liquid, \(D_{\text{Ni}}^{\text{ol/liq}}\). Experiments were performed from 1300 to 1600 °C and 1 atm to 3.0 GPa, using mid-ocean ridge basalt (MORB) glass surrounded by powdered olivine in graphite–Pt double capsules at high pressure and powdered MORB in crucibles fabricated from single crystals of San Carlos olivine at one atmosphere. In these experiments, pressure and temperature were varied in such a way that we produced a series of liquids, each with an approximately constant composition (~12, ~15, and ~21 wt% MgO). Previously, we used a similar approach to show that \(D_{\text{Ni}}^{\text{ol/liq}}\) for a liquid with ~18 wt% MgO is a strong function of temperature. Combining the new data presented here with our previous results allows us to separate the effects of temperature from composition. We fit our data based on a Ni–Mg exchange reaction, which yields \(\ln \left( {D_{\text{Ni}}^{\text{molar}} } \right) = \frac{{ -\Delta _{r(1)} H_{{T_{\text{ref}} ,P_{\text{ref}} }}^{ \circ } }}{RT} + \frac{{\Delta _{r(1)} S_{{T_{\text{ref}} ,P_{\text{ref}} }}^{ \circ } }}{R} - \ln \left( {\frac{{X_{\text{MgO}}^{\text{liq}} }}{{X_{{{\text{MgSi}}_{ 0. 5} {\text{O}}_{ 2} }}^{\text{ol}} }}} \right).\) Each subset of constant composition experiments displays roughly the same temperature dependence of \(D_{\text{Ni}}^{\text{ol/liq}}\) (i.e.,\(-\Delta _{r(1)} H_{{T_{\text{ref}} ,P_{\text{ref}} }}^{ \circ } /R\)) as previously reported for liquids with ~18 wt% MgO. Fitting new data presented here (15 experiments) in conjunction with our 13 previously published experiments (those with ~18 wt% MgO in the silicate liquid) to the above expression gives \(-\Delta _{r(1)} H_{{T_{\text{ref}} ,P_{\text{ref}} }}^{ \circ } /R\) = 3641 ± 396 (K) and \(\Delta _{r(1)} S_{{T_{\text{ref}} ,P_{\text{ref}} }}^{ \circ } /R\) = ? 1.597 ± 0.229. Adding data from the literature yields \(-\Delta _{r(1)} H_{{T_{\text{ref}} ,P_{\text{ref}} }}^{ \circ } /R\) = 4505 ± 196 (K) and \(\Delta _{r(1)} S_{{T_{\text{ref}} ,P_{\text{ref}} }}^{ \circ } /R\) = ? 2.075 ± 0.120, a set of coefficients that leads to a predictive equation for \(D_{\text{Ni}}^{\text{ol/liq}}\) applicable to a wide range of melt compositions. We use the results of our work to model the melting of peridotite beneath lithosphere of varying thickness and show that: (1) a positive correlation between NiO in magnesian olivine phenocrysts and lithospheric thickness is expected given a temperature-dependent \(D_{\text{Ni}}^{\text{ol/liq}} ,\) and (2) the magnitude of the slope for natural samples is consistent with our experimentally determined temperature dependence. Alternative processes to generate the positive correlation between NiO in magnesian olivines and lithospheric thickness, such as the melting of olivine-free pyroxenite, are possible, but they are not required to explain the observed correlation of NiO concentration in initially crystallizing olivine with lithospheric thickness.     

Thermodynamic properties of andradite and application to skarn with coexisting andradite and hedenbergite

The enthalpy of formation of andradite (Ca₃Fe₂Si₃O₁₂) has been estimated as-5,769.700 (±5) kJ/mol from a consideration of the calorimetric data on entropy (316.4 J/mol K) and of the experimental phaseequilibrium data on the reactions: 1 $$\begin{gathered} 9/2 CaFeSi_2 O_6 + O_2 = 3/2 Ca_3 Fe_2 Si_3 O_{12} + 1/2 Fe_3 O_4 + 9/2 SiO_2 (a) \hfill \\ Hedenbergite andradite magnetite quartz \hfill \\ \end{gathered} $$ 1 $$\begin{gathered} 4 CaFeSi_2 O_6 + 2 CaSiO_3 + O_2 = 2 Ca_3 Fe_2 Si_3 O_{12} + 4 SiO_2 (b) \hfill \\ Hedenbergite wollastonite andradite quartz \hfill \\ \end{gathered} $$ 1 $$\begin{gathered} 18 CaSiO_3 + 4 Fe_3 O_4 + O_2 = 6Ca_3 Fe_2 Si_3 O_{12} (c) \hfill \\ Wollastonite magnetite andradite \hfill \\ \end{gathered} $$ 1 $$\begin{gathered} Ca_3 Fe_2 Si_3 O_{12} = 3 CaSiO_3 + Fe_2 O_3 . (d) \hfill \\ Andradite pseudowollastonite hematite \hfill \\ \end{gathered} $$ and $$log f_{O_2 } = E + A + B/T + D(P - 1)/T + C log f_{O_2 } .$$ Oxygen-barometric scales are presented as follows: $$\begin{gathered} E = 12.51; D = 0.078; \hfill \\ A = 3 log X_{Ad} - 4.5 log X_{Hd} ; C = 0; \hfill \\ B = - 27,576 - 1,007(1 - X_{Ad} )^2 - 1,476(1 - X_{Hd} )^2 . \hfill \\ \end{gathered} $$ For the assemblage andradite (Ad)-hedenbergite (Hd)-magnetite-quartz: $$\begin{gathered} E = 13.98; D = 0.0081; \hfill \\ A = 4 log(X_{Ad} / X_{Hd} ); C = 0; \hfill \\ B = - 29,161 - 1,342.8(1 - X_{Ad} )^2 - 1,312(1 - X_{Hd} )^2 . \hfill \\ \end{gathered} $$ For the assemblage andradite-hedenbergite-wollastonite-quartz: 1 $$\begin{gathered} E = 13.98;{\text{ }}D = 0.0081; \hfill \\ A = 4\log (X_{Ad} /X_{Hd} );{\text{ C = 0;}} \hfill \\ B = - 29,161 - 1,342.8(1 - X_{Ad} )^2 - 1,312(1 - X_{Hd} )^2 . \hfill \\ \end{gathered} $$ For the assemblage andradite-hedenbergite-calcitequartz: 1 $$\begin{gathered} E = - 1.69;{\text{ }}D = - 0.199; \hfill \\ A = 4\log (X_{Ad} /X_{Hd} );{\text{ C = 2;}} \hfill \\ B = - 20,441 - 1,342.8(1 - X_{Ad} )^2 - 1,312(1 - X_{Hd} )^2 . \hfill \\ \end{gathered} $$ For the assemblage andradite-hedenbergite-wollastonite-calcite: 1 $$\begin{gathered} E = - 17.36;{\text{ }}D = - 0.403; \hfill \\ A = 4\log (X_{Ad} /X_{Hd} );{\text{ C = 4;}} \hfill \\ B = - 11,720 - 1,342.8(1 - X_{Ad} )^2 - 1,312(1 - X_{Hd} )^2 \hfill \\ \end{gathered} $$ The oxygen fugacity of formation of those skarns where andradite and hedenbergite assemblage is typical can be calculated by using the above equations. The oxygen fugacity of formation of this kind of skarn ranges between carbon dioxide/graphite and hematite/magnetite buffers. It increases from the inside zones to the outside zones, and appears to decrease with the ore-types in the order Cu, Pb?Zn, Fe, Mo, W(Sn) ore deposits.     

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.     

The Influence of Redox Reactions on the Uptake of Dissolved Ce by Suspended Fe and Mn Oxide Particles

Laboratory experiments were conducted to evaluate the partitioning ofrare earth elements (REE) between solution and suspended particles. Becauseof their strong tendency to complex, the REE can be used to study a varietyof marine processes and in particular particle scavenging. In this study, anemphasis was placed on examining abiotic redox processes that influence theuptake of dissolved Ce by particles. Batch sorption experiments wereconducted with REE and synthetic mineral phases over the range of pH4–9. The solutions varied in ionic strength between 0 and 0.7 M andconsisted of individual solutes (NaNO₃, NaCl, andNa₂SO₄), ionic mixtures that duplicate theseawater composition, and natural seawater. The uptake of REE from solutionwas also studied at a Pt electrode coated with using cyclic voltametry. Experimental results are consistent with uptake of dissolved Ce onto occurring by a combination of oxidativescavenging and surface complexation. The contribution of oxidativescavenging to the removal of Ce from solution is most pronounced at acidicpH, where the strictly trivalent REE exhibit little propensity for sorptiononto . Sorption of dissolved Ce onto FeOOH occursin a manner analogous to that of the other strictly trivalent REE and nocontribution from oxidative scavenging is observed on this mineral phase atlow pH. Our work also substantiates the hypothesis that anions in solution,particularly and Cl^-, aswell as those adsorbed on the surface of the particles, influence the extentof Ce uptake by . This revised version was published online in July 2006 with corrections to the Cover Date.     

High temperature calorimetry of sulfide systems

The standard enthalpies of formation of FeS (troilite), FeS₂ (pyrite), Co_0.9342S, Co₃S₄ (linnaeite), Co₉S₈ (cobalt pentlandite), CoS₂ (cattierite), CuS (covellite), and Cu₂S (chalcocite) have been determined by high temperature direct reaction calorimetry at temperatures between 700 K and 1021 K. The following results are reported: $$\Delta {\rm H}_{f,FeS}^{tr} = - 102.59 \pm 0.20kJ mol^{ - 1} ,$$ $$\Delta {\rm H}_{f,FeS}^{py} = - 171.64 \pm 0.93kJ mol^{ - 1} ,$$ $$\Delta {\rm H}_{f,Co_{0.934} S} = - 99.42 \pm 1.52kJ mol^{ - 1} ,$$ $$\Delta {\rm H}_{f,Co_9 S_8 }^{ptl} = - 885.66 \pm 16.83kJ mol^{ - 1} ,$$ $$\Delta {\rm H}_{f,Co_3 S_4 }^{In} = - 347.47 \pm 7.27kJ mol^{ - 1} ,$$ $$\Delta {\rm H}_{f,CoS_2 }^{ct} = - 150.94 \pm 4.85kJ mol^{ - 1} ,$$ $$\Delta {\rm H}_{f,Cu_2 S}^{cc} = - 80.21 \pm 1.51kJ mol^{ - 1} ,$$ and $$\Delta {\rm H}_{f,CuS}^{cv} = - 53.14 \pm 2.28kJ mol^{ - 1} ,$$ The enthalpy of formation of CuFeS₂ (chalcopyrite) from (CuS+FeS) and from (Cu+FeS₂) was determined by solution calorimetry in a liquid Ni_0.60S_0.40 melt at 1100 K. The results of these measurements were combined with the standard enthalpies of formation of CuS, FeS, and FeS₂, to calculate the standard enthalpy of formation of CuFeS₂. We found \(\Delta {\rm H}_{f,CuFeS_2 }^{ccp} = - 194.93 \pm 4.84kJ mol^{ - 1}\) . Our results are compared with earlier data given in the literature; generally the agreement is good and our values agree with previous estimates within the uncertainties present in both.     

Non-ideal Mg-Fe binary mixing in cordierite: Constraints from experimental data on Mg-Fe partitioning in garnet and cordierite and a reformulation of garnet-cordierite geothermometer

The non-ideal regular Mg-Fe binary in cordierite has been derived through multivariate linear regression of the expressionRT InK_D +(P- 1)ΔVK ₁ ⁰ , 298 along with updated subfegular mixing parameter of almandine-pyrope solution (Hackler and Wood 1989; Berman 1990). The data base used for multivariate analyses consists of published experimental data (n = 177) on Mg-Fe partitioning between garnet and cordierite in theP-T range 650–1050°C and 4–12 K bar. The non-ideality can be approximated by temperature-dependent Margules parameters. The retrieved values of ΔH_<T> ^o and ΔH_<T> ^o of exchange reaction between garnet and cordierite and enthalpy and entropy of mixing of Mg-Fe cordierite were combined with recent quaternary (Fe-Mg-Ca-Mn) mixing data in garnet to obtain the geothermometric expressions to determine temperature (T Kelvin): $$\begin{gathered} T(WH) = 6832 + 0.031(P - 1) - \{ 166(X_{Mg}^{Gt} )^2 - 506(X_{Fe}^{Gt} )^2 + 680X_{Fe}^{Gt} X_{Mg}^{Gt} + 336(X_{Ca} + X_{Mn} ) \hfill \\ (X_{Mg} - X_{Fe} )^{Gt} - 3300X_{Ca}^{Gt} - 358X_{Mn}^{Gt} \} + 954(X_{Fe} - X_{Mg} )^{Crd} /1.987\ln K_D + 3.41 + 1.5X_{Ca}^{Gt} \hfill \\ + 1.23(X_{Fe} - X_{Mg} )^{Crd} \hfill \\ \end{gathered} $$ $$\begin{gathered} T(Br) = 6920 + 0.031(p - 1) - \{ 18(X_{Mg}^{Gt} )^2 - 296(X_{Fe}^{Gt} )^2 + 556X_{Fe}^{Gt} X_{Mg}^{Gt} - 6339X_{Ca}^{Gt} X_{Mg}^{Gt} \hfill \\ - 99(X_{Ca}^{Gt} )^2 + 4687X_{Ca}^{Gt} (X_{Mg} - X_{Fe}^{Gt} ) - 4269X_{Ca}^{Gt} X_{Fe}^{Gt} - 358X_{Mn}^{Gt} \} + 640(X_{Fe} - X_{Mg} )^{Crd} \hfill \\ + 1.90X_{Ca}^{Gt} (X_{Mg} - X_{Ca} )^{Gt} . \hfill \\ \end{gathered} $$      

Opening and resetting temperatures in heating geochronological systems

We present a theoretical model for diffusive daughter isotope loss in radiochronological systems with increasing temperature. It complements previous thermochronological models, which focused on cooling, and allows for testing opening and resetting of radiochronometers during heating. The opening and resetting temperatures are, respectively,

Thermal expansion and spontaneous strain associated with the normal-incommensurate phase transition in melilites

The linear thermal expansions of åkermanite (Ca₂MgSi₂O₇) and hardystonite (Ca₂ZnSi₂O₇) have been measured across the normal-incommensurate phase transition for both materials. Least-squares fitting of the high temperature (normal phase) data yields expressions linear in T for the coefficients of instantaneous linear thermal expansion, $$\alpha _1 = \frac{1}{l}\frac{{dl}}{{dT}}$$ for åkermanite: $$\begin{gathered} \alpha _{[100]} = 6.901(2) \times 10^{ - 6} + 1.834(2) \times 10^{ - 8} T \hfill \\ \alpha _{[100]} = - 2.856(1) \times 10^{ - 6} + 11.280(1) \times 10^{ - 8} T \hfill \\ \end{gathered} $$ for hardystonite: $$\begin{gathered} \alpha _{[100]} = 15.562(5) \times 10^{ - 6} - 1.478(3) \times 10^{ - 8} T \hfill \\ \alpha _{[100]} = - 11.115(5) \times 10^{ - 6} + 11.326(3) \times 10^{ - 8} T \hfill \\ \end{gathered} $$ Although there is considerable strain for temperatures within 10° C of the phase transition, suggestive of a high-order phase transition, there appears to be a finite ΔV of transition, and the phase transition is classed as “weakly first order”.     

The dependence of the Fe2+-Mg cation-partitioning between olivine and basaltic liquid on pressure,temperature and composition

The equilibrium partitioning of Fe²⁺ and Mg between olivine and liquid along a liquid line of descent has been determined for a calc-alkaline system, ranging in composition from picritic to andesitic. Experiments were conducted between 1000–1450° C and 1 bar to 30 kbar. Within the compositional range investigated and , the compositional dependence of the Fe²⁺ and Mg partitioning is a function of the Mg-content of the liquid. The Mg-content of the liquid correlates strongly with temperature. The variation of the Fe²⁺ and Mg partitioning were therefore evaluated individualy as functions of composition and temperature alone. The composition dependence of the cation-partitioning coefficients (Kd) is given by the following two equations: