The effect of H<Subscript>2</Subscript>O on the olivine liquidus of basaltic melts: experiments and thermodynamic models

We designed and carried out experiments to investigate the effect of H₂O on the liquidus temperature of olivine-saturated primitive melts. The effect of H₂O was isolated from other influences by experimentally determining the liquidus temperatures of the same melt composition with various amounts of H₂O added. Experimental data indicate that the effect of H₂O does not depend on pressure or melt composition in the basaltic compositional range. The influence of H₂O on melting point lowering can be described as a polynomial function This expression can be used to account for the effect of H₂O on olivine-melt thermometers, and can be incorporated into fractionation models for primitive basalts. The non-linear effect of H₂O indicates that incorporation of H₂O in silicate melts is non-ideal, and involves interaction between H₂O and other melt components. The simple speciation approach that seems to account for the influence of H₂O in simple systems (albite-H₂O, diopside-H₂O) fails to describe the mixing behavior of H₂O in multi-component silicate melts. However, a non-ideal solution model that treats the effect of H₂O addition as a positive excess free energy can be fitted to describe the effect of melting point lowering.     

Thermodynamic mixing properties of Mg(Al,Cr)2O4 spinel crystalline solution at high temperatures and pressures

Three Al-Cr exchange isotherms at 1,250°, 1,050°, and 796° between Mg(Al, Cr)₂O₄ spinel and (Al, Cr)₂O₃ corundum crystalline solutions have been studied experimentally at 25 kbar pressure. Starting from gels of suitable bulk compositions, close approach to equilibrium has been demonstrated in each case by time studies. Using the equation of state for (Al, Cr)₂O₃ crystalline solution (Chatterjee et al. 1982a) and assuming that the Mg(Al, Cr)₂O₄ can be treated in terms of the asymmetric Margules relation, the exchange isotherms were solved for Δ G ^*, and . The best constrained data set from the 1,250° C isotherm clearly shows that the latter two quantities do not overlap within three standard deviations, justifying the choice of asymmetric Margules relation for describing the excess mixing properties of Mg(Al, Cr)₂O₄ spinels. Based on these experiments, the following polybaric-polythermal equation of state can be formulated: , P expressed in bars, T in K, G _m ^ex and W _G,i ^Sp in joules/mol. Temperature-dependence of G _m ^ex is best constrained in the range 796–1,250° C; extrapolation beyond that range would have to be done with caution. Such extrapolation to lower temperature shows tentatively that at 1 bar pressure the critical temperature, T _c, of the spinel solvus is 427° C, with dT_c/dP≈1.3 K/kbar. The critical composition, X _c, is 0.42 , and changes barely with pressure. Substantial error in calculated phase diagrams will result if the significant positive deviation from ideality is ignored for Al-Cr mixing in such spinels.     

The H<Subscript>2</Subscript>O solubility of alkali basaltic melts: an experimental study

Experiments were conducted to determine the water solubility of alkali basalts from Etna, Stromboli and Vesuvius volcanoes, Italy. The basaltic melts were equilibrated at 1,200°C with pure water, under oxidized conditions, and at pressures ranging from 163 to 3,842 bars. Our results show that at pressures above 1 kbar, alkali basalts dissolve more water than typical mid-ocean ridge basalts (MORB). Combination of our data with those from previous studies allows the following simple empirical model for the water solubility of basalts of varying alkalinity and fO₂ to be derived: \textH ₂ \textO( \textwt% ) = \text H ₂ \textO_\textMORB ( \textwt% ) + ( 5.84 ×10 ^{- 5} *\textP - 2.29 ×10 ^{- 2} ) ×( \textNa₂ \textO + \textK₂ \textO )( \textwt% ) + 4.67 ×10 ^{- 2} ×\Updelta \textNNO - 2.29 ×10 ^{- 1} {\text{H}}_{ 2} {\text{O}}\left( {{\text{wt}}\% } \right) = {\text{ H}}_{ 2} {\text{O}}_{\text{MORB}} \left( {{\text{wt}}\% } \right) + \left( {5.84 \times 10^{ - 5} *{\text{P}} - 2.29 \times 10^{ - 2} } \right) \times \left( {{\text{Na}}_{2} {\text{O}} + {\text{K}}_{2} {\text{O}}} \right)\left( {{\text{wt}}\% } \right) + 4.67 \times 10^{ - 2} \times \Updelta {\text{NNO}} - 2.29 \times 10^{ - 1} where H₂O_MORB is the water solubility at the calculated P, using the model of Dixon et al. (1995). This equation reproduces the existing database on water solubilities in basaltic melts to within 5%. Interpretation of the speciation data in the context of the glass transition theory shows that water speciation in basalt melts is severely modified during quench. At magmatic temperatures, more than 90% of dissolved water forms hydroxyl groups at all water contents, whilst in natural or synthetic glasses, the amount of molecular water is much larger. A regular solution model with an explicit temperature dependence reproduces well-observed water species. Derivation of the partial molar volume of molecular water using standard thermodynamic considerations yields values close to previous findings if room temperature water species are used. When high temperature species proportions are used, a negative partial molar volume is obtained for molecular water. Calculation of the partial molar volume of total water using H₂O solubility data on basaltic melts at pressures above 1 kbar yields a value of 19 cm³/mol in reasonable agreement with estimates obtained from density measurements.     

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

Thermodynamic Properties of Brucite Determined by Solubility Studies and Their Significance to Nuclear Waste Isolation

Solubility experiments were conducted for the dissolution reaction of brucite, Mg(OH)₂ (cr): Experiments were conducted from undersaturation in deionized (DI) water and 0.010–4.4 m NaCl solutions at 22.5°C. In addition, brucite solubility was measured from supersaturation in an experiment in which brucite was precipitated via dropwise addition of 0.10 m NaOH into a 0.10 m MgCl₂ solution also at 22.5°C. The attainment of the reversal in equilibrium was demonstrated in this study. The solubility constant at 22.5°C at infinite dilution calculated from the experimental results from the direction of supersaturation by using the specific interaction theory (SIT) is: with a corresponding value of 17.0 ± 0.2 (2σ) when extrapolated to 25°C. The dimensionless standard chemical potential (μ°/RT) of brucite derived from the solubility data in 0.010 m to 4.4 m NaCl solutions from undersaturation extrapolated to 25°C is −335.76 ± 0.45 (2σ), with the corresponding Gibbs free energy of formation of brucite, , being −832.3 ± 1.1 (2σ) kJ mol⁻¹. In combination with the auxiliary thermodynamic data, the is calculated to be 17.1 ± 0.2 (2σ), based on the above Gibbs free energy of formation for brucite. This study recommends an average value of 17.05 ± 0.2 in logarithmic unit as solubility constant of brucite at 25°C, according to the values from both supersaturation and undersaturation. Sandia National Laboratories is a multiprogram laboratory operated by Sandia Corporation, a Lockheed Martin Company, for the United States Department of Energy’s National Nuclear Security Administration under Contract DE-AC04-94AL85000.     

Aluminum partitioning between olivine and ultrabasic silicate liquid to 6 GPa

Trace element analyses of 1-atm and high-pressure experiments show that in komatiite and peridotite, the olivine (OL)/liquid (L) distribution coefficient for Al₂O₃ ( ) increases with pressure and temperature. Olivine in equilibrium with liquid accepts as much as 0.2 wt% Al₂O₃ in solution at 6 GPa. Convergence to equilibrium compositions at this high level is shown by cation diffusion of Al into synthetic forsterite crystals of low-Al contents in the presence of melt. Convergence to low-Al equilibrium compositions at lower P and T is shown by diffusion of Al out of synthetic forsterite with high initial Al content. Isobaric and isothermal experimental data subsets reveal that temperature and pressure variations both have real effects on . Variation in silicate melt composition has no detectable effect on within the limited range of experimentally investigated mixtures. Least-squares regression for 24 experiments, using komatiite and peridotite, performed at 1 atm to 6 GPa and 1300 to 1960°C, gives the best fit equation: Increase in with increasingly higher-pressure melting is consistent with incorporation of a spinel-like component of low molar volume into olivine, although other substitutions possibly involving more complex coupling cannot be ruled out. High P-T ultrabasic melting residues, if pristine, may be recognized by the high calculated from microprobe analyses of Al₂O₃ concentrations in residual olivines and estimated Al₂O₃ concentration in the last liquid removed. In general the low levels of Al in natural olivine from mantle xenoliths suggest that pristine residues are rarely recovered.     

Energetics of hydration of cordierite and water barometry in cordierite-granulites

The Gibbs free energy and volume changes attendant upon hydration of cordierites in the system magnesian cordierite-water have been extracted from the published high pressure experimental data at \(P_{{\text{H}}_{\text{2}} {\text{O}}} \) =P _total, assuming an ideal one site model for H₂O in cordierite. Incorporating the dependence of ΔG and ΔV on temperature, which was found to be linear within the experimental conditions of 500°–1,000°C and 1–10,000 bars, the relation between the water content of cordierite and P, T and \(f_{{\text{H}}_{\text{2}} {\text{O}}} \) has been formulated as $$\begin{gathered} X_{{\text{H}}_{\text{2}} {\text{O}}}^{{\text{crd}}} = \hfill \\ \frac{{f_{{\text{H}}_{\text{2}} {\text{O}}}^{{\text{P, T}}} }}{{\left[ {{\text{exp}}\frac{1}{{RT}}\left\{ {64,775 - 32.26T + G_{{\text{H}}_{\text{2}} {\text{O}}}^{{\text{1, }}T} - P\left( {9 \times 10^{ - 4} T - 0.5142} \right)} \right\}} \right] + f_{{\text{H}}_{\text{2}} {\text{O}}}^{{\text{P, T}}} }} \hfill \\ \end{gathered} $$ The equation can be used to compute H₂O in cordierites at \(P_{{\text{H}}_{\text{2}} {\text{O}}} \) <1. Our results at different P, T and partial pressure of water, assuming ideal mixing of H₂O and CO₂ in the vapour phase, are in very good agreement with the experimental data of Johannes and Schreyer (1977, 1981). Applying the formulation to determine \(X_{{\text{H}}_{\text{2}} {\text{O}}}^{{\text{crd}}} \) in the garnet-cordierite-sillimanite-plagioclase-quartz granulites of Finnish Lapland as a test case, good agreement with the gravimetrically determined water contents of cordierite was obtained. Pressure estimates, from a thermodynamic modelling of the Fe-cordierite — almandine — sillimanite — quartz equilibrium at \(P_{{\text{H}}_{\text{2}} {\text{O}}} = 0\) and \(P_{{\text{H}}_{\text{2}} {\text{O}}} \) =P_total, for assemblages from South India, Scottish Caledonides, Daly Bay and Hara Lake areas are compatible with those derived from the garnetplagioclase-sillimanite-quartz geobarometer.     

Non-ideal mixing in the phlogopite-annite binary: constraints from experimental data on Mg−Fe partitioning and a reformulation of the biotite-garnet geothermometer

The existing experimental data [Ferry and Spear 1978; Perchuk and Lavrent'eva 1983] on Mg?Fe partitioning between garnet and biotite are disparate. The underlying assumption of ideal Mg?Fe exchange between the minerals has been examined on the basis of recently available thermochemical data. Using the updated mixing parameters for the pyrope-almandine asymmetric regular solution as inputs [Ganguly and Saxena 1984; Hackler and Wood 1984], thermodynamic analysis points to non-ideal mixing in the phlogopite-annite binary in the temperature range of 550°C–950°C. The non-ideality can be approximated by a temperature-independent, one constant Margules parameter. The retrieved values for enthalpy of mixing for Mg?Fe biotites and the standard state enthalpy and entropy changes of the exchange reaction were combined with existing thermochemical data on grossular-pyrope and grossular-almandine binaries to obtain geothermometric expressions for Mg?Fe fractionation between biotite and garnet. [T in K] $$\begin{gathered} {\text{T(HW) = [20286 + 0}}{\text{.0193P - \{ 2080(X}}_{{\text{Mg}}}^{{\text{Gt}}} {\text{)}}^{\text{2}} {\text{ - 6350(X}}_{{\text{Fe}}}^{{\text{Gt}}} {\text{)}}^{\text{2}} \hfill \\ {\text{ - 13807(X}}_{{\text{Ca}}}^{{\text{Gt}}} {\text{)(1 - X}}_{{\text{Mn}}}^{{\text{Gt}}} {\text{) + 8540(X}}_{{\text{Fe}}}^{{\text{Gt}}} {\text{)(X}}_{{\text{Mg}}}^{{\text{Gt}}} {\text{)(1 - X}}_{{\text{Mn}}}^{{\text{Gt}}} {\text{)}} \hfill \\ {\text{ + 4215(X}}_{{\text{Ca}}}^{{\text{Gt}}} {\text{)(X}}_{{\text{Mg}}}^{{\text{Gt}}} {\text{ - X}}_{{\text{Fe}}}^{{\text{Gt}}} {\text{)\} + 4441}}{{{\text{(2X}}_{{\text{Mg}}}^{{\text{Bt}}} {\text{ - 1)]}}} \mathord{\left/ {\vphantom {{{\text{(2X}}_{{\text{Mg}}}^{{\text{Bt}}} {\text{ - 1)]}}} {{\text{[13}}{\text{.138}}}}} \right. \kern-\nulldelimiterspace} {{\text{[13}}{\text{.138}}}} \hfill \\ {\text{ + 8}}{\text{.3143 InK}}_{\text{D}} {\text{ + 6}}{\text{.276(X}}_{{\text{Ca}}}^{{\text{Gt}}} ){\text{(1 - X}}_{{\text{Mn}}}^{{\text{Gt}}} )] \hfill \\ {\text{T(GS) = [13538 + 0}}{\text{.0193P - \{ 837(X}}_{{\text{Mg}}}^{{\text{Gt}}} )^{\text{2}} {\text{ - 10460(X}}_{{\text{Fe}}}^{{\text{Gt}}} )^2 \hfill \\ {\text{ - 13807(X}}_{{\text{Ca}}}^{{\text{Gt}}} )(1{\text{ - X}}_{{\text{Mn}}}^{{\text{Gt}}} {\text{) + 19246(X}}_{{\text{Fe}}}^{{\text{Gt}}} ){\text{(X}}_{{\text{Mg}}}^{{\text{Gt}}} ){\text{(1 - X}}_{{\text{Mn}}}^{{\text{Gt}}} ) \hfill \\ {\text{ }}{{{\text{ + 5649(X}}_{{\text{Ca}}}^{{\text{Gt}}} ){\text{(X}}_{{\text{Mg}}}^{{\text{Gt}}} {\text{ - X}}_{{\text{Fe}}}^{{\text{Gt}}} ){\text{\} + 7972(2X}}_{{\text{Mg}}}^{{\text{Bt}}} {\text{ - 1)]}}} \mathord{\left/ {\vphantom {{{\text{ + 5649(X}}_{{\text{Ca}}}^{{\text{Gt}}} ){\text{(X}}_{{\text{Mg}}}^{{\text{Gt}}} {\text{ - X}}_{{\text{Fe}}}^{{\text{Gt}}} ){\text{\} + 7972(2X}}_{{\text{Mg}}}^{{\text{Bt}}} {\text{ - 1)]}}} {{\text{[6}}{\text{.778}}}}} \right. \kern-\nulldelimiterspace} {{\text{[6}}{\text{.778}}}} \hfill \\ {\text{ + 8}}{\text{.3143InK}}_{\text{D}} {\text{ + 6}}{\text{.276(X}}_{{\text{Ca}}}^{{\text{Gt}}} )(1{\text{ - X}}_{{\text{Mn}}}^{{\text{Gt}}} )] \hfill \\ \end{gathered} $$ The reformulated geothermometer is an improvement over existing biotite-garnet geothermometers because it reconciles the experimental data sets on Fe?Mg partitioning between the two phases and is based on updated activity-composition relationship in Fe?Mg?Ca garnet solid solutions.     

Amphibole composition in tonalite as a function of pressure: an experimental calibration of the Al-in-hornblende barometer

The Al-in-hornblende barometer, which correlates Al^tot content of magmatic hornblende linearly with crystallization pressure of intrusion (Hammarstrom and Zen 1986), has been calibrated experimentally under water-saturated conditions at pressures of 2.5–13 kbar and temperatures of 700–655°C. Equilibration of the assemblage hornlende-biotite-plagioclase-orthoclasequartz-sphene-Fe-Ti-oxide-melt-vapor from a natural tonalite 15–20° above its wet solidus results in hornblende compositions which can be fit by the equation: P(±0.6 kbar) = –3.01 + 4.76 Al _hbl ^tot r ²=0.99, where Al^tot is the total Al content of hornblende in atoms per formula unit (apfu). Al^tot increase with pressure can be ascribed mainly to a tschermak-exchange ( ) accompanied by minor plagioclase-substitution ( ). This experimental calibration agrees well with empirical field calibrations, wherein pressures are estimated by contact-aureole barometry, confirming that contact-aureole pressures and pressures calculated by the Al-in-hornblende barometer are essentially identical. This calibration is also consistent with the previous experimental calibration by Johnson and Rutherford (1989b) which was accomplished at higher temperatures, stabilizing the required buffer assemblage by use of mixed H₂O-CO₂ fluids. The latter calibration yields higher Al^tot content in hornblendes at corresponding pressures, this can be ascribed to increased edenite-exchange ( ) at elevated temperatures. The comparison of both experimental calibrations shows the important influence of the fluid composition, which affects the solidus temperature, on equilibration of hornblende in the buffering phase assemblage.     

Experimental determination of activities in Fe−Mg olivine at 1400 K

The exchage equilibrium has been used to measure activity-composition relations along the olivine join FeSi_0.5O₂−MgSi_0.5O₂ at 1400 K and 1 atm pressure. Equilibrium Fe−Mg partitioning between the two phases was determined by reversing the compositions of olivine coexisting with oxide and matallic iron over the composition range Fo₂₃ to Fo₉₂. A detailed study of the thermodynamic properties of the oxide phase has recently been made by Srečec et al. and we have confirmed their results in the composition range of interest. Application of the oxide data to the exchange equilibrium enables the properties of olivine to be determined. Within experimental uncertainly (Fe, Mg)Si_0.5O₂ olivine can, at 1400 K, be treated as a symmetric solution with W _Fe-Mg ^o1 of 3.7±0.8 kJ/mol. The data permit the presence of only very slight asymmetry in the series. The data do not support recent assertions that olivine is highly non-ideal (W≈10 kJ/mol) under these conditions.     

IR calibrations for water determination in olivine,r-GeO<Subscript>2</Subscript>, and SiO<Subscript>2</Subscript> polymorphs

Mineral-specific IR absorption coefficients were calculated for natural and synthetic olivine, SiO₂ polymorphs, and GeO₂ with specific isolated OH point defects using quantitative data from independent techniques such as proton–proton scattering, confocal Raman spectroscopy, and secondary ion mass spectrometry. Moreover, we present a routine to detect OH traces in anisotropic minerals using Raman spectroscopy combined with the “Comparator Technique”. In case of olivine and the SiO₂ system, it turns out that the magnitude of ε for one structure is independent of the type of OH point defect and therewith the peak position (quartz ε = 89,000 ± 15,000  \textl \textmol_{\textH2\textO}^-1 \textcm^-2\text{l}\,\text{mol}_{{\text{H}_2}\text{O}}^{-1}\,\text{cm}^{-2}), but it varies as a function of structure (coesite ε = 214,000 ± 14,000  \textl \textmol_{\textH2\textO}^-1 \textcm^-2\text{l}\,\text{mol}_{{\text{H}_2}\text{O}}^{-1}\,\text{cm}^{-2}; stishovite ε = 485,000 ± 109,000  \textl \textmol_{\textH2\textO}^-1 \textcm^-2\text{l}\,\text{mol}_{{\text{H}_2}\text{O}}^{-1}\,\text{cm}^{-2}). Evaluation of data from this study confirms that not using mineral-specific IR calibrations for the OH quantification in nominally anhydrous minerals leads to inaccurate estimations of OH concentrations, which constitute the basis for modeling the Earth’s deep water cycle.     

Some geobarometers involving cordierite in the FeO-Al2O3-SiO2(±H2O) system: refinements,thermodynamic calibration,and applicability in granulite facies rocks

Five geobarometers involving cordierite have been formulated for quantitative pressure sensing in high grade metapelites. The relevant reactions in the FeO-Al₂O₃-SiO₂ (±H₂O) system are based on the assemblages (A) cordierite-garnet-sillimanite-quartz, (B) cordierite-spinel-quartz, (C) cordierite-garnet-spinel-sillimanite, (D) cordierite-garnet-orthopyroxene-quartz and (E) cordierite-orthopyroxene-sillimanite-quartz. Application of the barometric formulations to a large number of granulite grade rocks indicates that the cordierite-garnet-sillimanite-quartz equilibrium is widely applicable and registers pressures which are in good agreement with the “consensus” pressure estimates. The dispersion in the computed P values, expressed as one standard deviation, is within ±1.2 kbar. The geobarometers (B) and (C) also yield pressures which are reasonable and compare well with those computed from equilibrium (A). The estimated pressures from (D) and (E), both involving orthopyroxene, are at variance with these estimates. It has been argued that the discrepancy in pressures obtained from these geobarometers stems from an inadequate knowledge of activity-composition relations and/or errors in input thermodynamic data of aluminous orthopyroxene. The convergence of pressure values estimated from the barometric formulations, especially (A), (B) and (C), implies that the present formulations are more dependable than the existing formulations and are also capable of setting limits on P values in response to varying $$\begin{gathered} {\text{1/2Fe}}_{\text{2}} {\text{Al}}_{\text{4}} {\text{Si}}_{\text{5}} {\text{O}}_{{\text{18}}} \hfill \\ {\text{ = 1/3Fe}}_{\text{3}} {\text{Al}}_{\text{2}} {\text{Si}}_{\text{3}} {\text{O}}_{{\text{12}}} {\text{ + 2/3Al}}_{\text{2}} {\text{SiO}}_{\text{5}} {\text{ + 5/6SiO}}_{\text{2}} {\text{. (A)}} \hfill \\ {\text{1/2Fe}}_{\text{2}} {\text{Al}}_{\text{4}} {\text{Si}}_{\text{5}} {\text{O}}_{{\text{18}}} {\text{ = FeAl}}_{\text{2}} {\text{O}}_{\text{4}} {\text{ + 5/2SiO}}_{\text{2}} {\text{. (B)}} \hfill \\ {\text{Fe}}_{\text{2}} {\text{Al}}_{\text{4}} {\text{Si}}_{\text{5}} {\text{O}}_{{\text{18}}} {\text{ + FeAl}}_{\text{2}} {\text{O}}_{\text{4}} \hfill \\ = {\text{Fe}}_{\text{3}} {\text{Al}}_{\text{2}} {\text{Si}}_{\text{3}} {\text{O}}_{{\text{12}}} {\text{ + 2Al}}_{\text{2}} {\text{SiO}}_{\text{5}} {\text{. (C)}} \hfill \\ {\text{1/2Fe}}_{\text{2}} {\text{Al}}_{\text{4}} {\text{Si}}_{\text{5}} {\text{O}}_{{\text{18}}} {\text{ + Fe}}_{\text{2}} {\text{Si}}_{\text{2}} {\text{O}}_{\text{6}} \hfill \\ = {\text{Fe}}_{\text{3}} {\text{Al}}_{\text{2}} {\text{Si}}_{\text{3}} {\text{O}}_{{\text{12}}} {\text{ + 3/2SiO}}_{\text{2}} .{\text{ (D)}} \hfill \\ {\text{1/2Fe}}_{\text{2}} {\text{Al}}{}_{\text{4}}{\text{Si}}_{\text{5}} {\text{O}}_{{\text{18}}} \hfill \\ = 1/2{\text{Fe}}_{\text{2}} {\text{Si}}_{\text{2}} {\text{O}}_{\text{6}} {\text{ + Al}}_{\text{2}} {\text{SiO}}_{\text{5}} {\text{ + 1/2SiO}}_{\text{2}} .{\text{ (E)}} \hfill \\ \end{gathered}$$ . The present communication addresses the calibration, applicability and reliability of these barometers with reference to granulite facies metapelites.     

The combined effect of F and H2O on interdiffusion between peralkaline dacitic and rhyolitic melts

The effective binary diffusion coefficient (EBDC) of silicon has been measured during the interdiffusion of peralkaline, fluorine-bearing (1.3 wt% F), hydrous (3.3 and 6 wt% H₂O), dacitic and rhyolitic melts at 1.0 GPa and temperatures between 1100°C and 1400°C. From Boltzmann-Matano analysis of diffusion profiles the diffusivity of silicon at 68 wt% SiO₂ can be described by the following Arrhenius equations (with standard errors): $$\begin{gathered} {\text{with 1}}{\text{.3 wt\% F and 3}}{\text{.3\% H}}_{\text{2}} {\text{O:}} \hfill \\ {\text{D}}_{{\text{Si}}} = \begin{array}{*{20}c} { + {\text{3}}{\text{.59}}} \\ {{\text{3}}{\text{.66}} \times {\text{10}}^{ - {\text{9}}} } \\ { - {\text{1}}{\text{.86}}} \\ \end{array} {\text{exp}}\left( {{{ - {\text{86}}{\text{.1}} \pm {\text{8}}{\text{.9}}} \mathord{\left/ {\vphantom {{ - {\text{86}}{\text{.1}} \pm {\text{8}}{\text{.9}}} {{\text{RT}}}}} \right. \kern-\nulldelimiterspace} {{\text{RT}}}}} \right) \hfill \\ {\text{with 1}}{\text{.3 wt\% F and 6}}{\text{.0\% H}}_{\text{2}} {\text{O:}} \hfill \\ {\text{D}}_{{\text{Si}}} = \begin{array}{*{20}c} { + {\text{3}}{\text{.59}}} \\ {{\text{3}}{\text{.51}} \times {\text{10}}^{ - {\text{8}}} } \\ { - {\text{1}}{\text{.77}}} \\ \end{array} {\text{exp}}\left( {{{ - {\text{109}}{\text{.5}} \pm {\text{8}}{\text{.9}}} \mathord{\left/ {\vphantom {{ - {\text{109}}{\text{.5}} \pm {\text{8}}{\text{.9}}} {{\text{RT}}}}} \right. \kern-\nulldelimiterspace} {{\text{RT}}}}} \right) \hfill \\ \end{gathered} $$ where D is in m²s^?1 and activation energies are in kJ/mol. Diffusivities measured at 64 and 72 wt% SiO₂ are only slightly different from those at 68 wt% SiO₂ and frequently all measurements are within error of each other. Silicon, aluminum, iron, magnesium, and calcium EBDCs were also calculated from diffusion profiles by error function inversion techniques assuming constant diffusivity. With one exception, silicon EBDCs calculated by error function techniques are within error of Boltzmann-Matano EBDCs. Average diffusivities of Fe, Mg, and Ca were within a factor of 2.5 of silicon diffusivities whereas Al diffusivities were approximately half those of silicon. Alkalies diffused much more rapidly than silicon and non-alkalies, however their diffusivities were not quantitatively determined. Low activation energies for silicon EBDCs result in rapid diffusion at magmatic temperatures. Assuming that water and fluorine exert similar effects on melt viscosity at high temperatures, the viscosity can be calculated and used in the Eyring equation used to determine diffusivities, typically to within a factor of three of those measured in this study. This correlation between viscosity and diffusivity can be inverted to calculate viscosities of fluorine- and water-bearing granitic melts at magmatic temperatures; these viscosities are orders of magnitude below those of hydrous granitic melts and result in more rapid and effective separation of granitic magmas from partially molten source rocks. Comparison of Arrhenius parameters for diffusion measured in this study with Arrhenius parameters determined for diffusion in similar compositions at the same pressure demonstrates simple relationships between Arrhenius parameters, activation energy-E_a, kJ/mol, pre-exponential factor-D_o, m²s^?1, and the volatile, X=F or OH^?, to oxygen, O, ratio of the melt {(X/X+O)}: $$\begin{gathered} {\text{E}}a = - {\text{1533\{ }}{{\text{X}} \mathord{\left/ {\vphantom {{\text{X}} {\left( {{\text{X}} + {\text{O}}} \right)}}} \right. \kern-\nulldelimiterspace} {\left( {{\text{X}} + {\text{O}}} \right)}}{\text{\} }} + {\text{213}}{\text{.3}} \hfill \\ {\text{D}}_{\text{O}} = {\text{2}}{\text{.13}} \times {\text{10}}^{ - {\text{6}}} {\text{exp}}\left[ { - {\text{6}}{\text{.5\{ }}{{\text{X}} \mathord{\left/ {\vphantom {{\text{X}} {\left( {{\text{X}} + {\text{O}}} \right)}}} \right. \kern-\nulldelimiterspace} {\left( {{\text{X}} + {\text{O}}} \right)}}{\text{\} }}} \right] \hfill \\ \end{gathered} $$ These relationships can be used to estimate diffusion in various melts of dacitic to rhyolitic composition containing both fluorine and water. Calculations for the contamination of rhyolitic melts by dacitic enclaves at 800°C and 700°C provide evidence for the virtual inevitability of diffusive contamination in hydrous and fluorine-bearing magmas if they undergo magma mixing of any form.     

The Equation of State for Caspian Sea Waters

The density ρ of Caspian Sea waters was measured as a function of temperature (273.15–343.15) K at conductivity salinities of 7.8 and 11.3 using the Anton-Paar Densitometer. Measurements were also made on one of the samples (S = 11.38) diluted with water as a function of temperature (T = 273.15–338.15 K) and salinity (2.5–11.3). These latter results have been used to develop an equation of state for the Caspian Sea (σ = ±0.007 kg m⁻³)

Clinoeulite (magnesian clinoferrosilite) in a eulysite of a metamorphosed iron formation in the vredefort structure,South Africa

A unique clinopyroxene (En₁₉Fs₇₈Wo₃), clinoeulite, space group P2₁/c, $${\text{(Fe}}_{{\text{1}}{\text{.48}}} {\text{Mg}}_{{\text{0}}{\text{.37}}} {\text{Mn}}_{{\text{0}}{\text{.08}}}^{{\text{2 + }}} {\text{Ca}}_{{\text{0}}{\text{.05}}} {\text{Al}}_{{\text{0}}{\text{.01}}} {\text{)}}_{{\text{1}}{\text{.99}}} {\text{ [Si}}_{{\text{2}}{\text{.01}}} {\text{O6],}}$$ contains sharp exsolution lamellae of ferroaugite (En₁₇Fs₄₃Wo₄₀) from which the former presence of a ferropigeonite near En₁₇Fs₇₀Wo₁₃ can be calculated. This two-pyroxene intergrowth is the main component of a eulysite containing also magnetite, olivine (Fo₉Fa₈₆Te₅), quartz, oligoclase-K feldspar inter-growth, and retrograde cummingtonite with about 76 % grunerite end member. The occurrence of this most unusual rock type in the center of the Vredefort structure is attributed to a period of high-temperature metamorphism (at least 800 °–850 °C) which was followed by hot deformation of the rock during the Vredefort event thus probably preventing the common formation of orthopyroxene through pigeonite exsolution and inversion upon cooling. After this tectonic deformation, the rock recrystallized within the low-temperature stability range of clinoeulite to yield fine annealing textures. Late-stage equilibria at temperatures well below 500 °C include the complete unmixing of a former high-temperature anorthoclase, a Mg/Fe redistribution in the clinoeulite and olivine and, with the introduction of water, the partial formation of cummingtonite through reaction of clinoeulite, olivine, and quartz. During weathering the olivine was transformed to a nearly opaque, anhydrous ferrisilicate which, except for the change of Fe²⁺ to Fe³⁺ and the oxygen introduction, largely retained its original chemistry.     

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:

Optical absorption and mössbauer spectra of purple and green yoderite,a kyanite-related mineral

Mössbauer and polarized optical absorption spectra of the kyanite-related mineral yoderite were recorded. Mössbauer spectra of the purple (PY) and green yoderite (GY) from Mautia Hill, Tanzania, show that the bulk of the iron is Fe³⁺ in both varieties, with Fe²⁺/(Fe²⁺+Fe³⁺) ratios near 0.05. Combining this result with new microprobe data for PY and with literature data for GY gives the crystallochemical formulae: $$\begin{gathered} ({\text{Mg}}_{{\text{1}}{\text{.95}}} {\text{Fe}}_{{\text{0}}{\text{.02}}}^{{\text{2 + }}} {\text{Mn}}_{{\text{0}}{\text{.01}}}^{{\text{2 + }}} {\text{Fe}}_{{\text{0}}{\text{.34}}}^{{\text{3 + }}} {\text{Mn}}_{{\text{0}}{\text{.07}}}^{{\text{3 + }}} {\text{Ti}}_{{\text{0}}{\text{.01}}} {\text{Al}}_{{\text{3}}{\text{.57}}} )_{5.97}^{[5,6]} \hfill \\ {\text{Al}}_{{\text{2}}{\text{.00}}}^{{\text{[5]}}} [({\text{Si}}_{{\text{3}}{\text{.98}}} {\text{P}}_{{\text{0}}{\text{.03}}} ){\text{O}}_{{\text{18}}{\text{.02}}} ({\text{OH)}}_{{\text{1}}{\text{.98}}} ] \hfill \\ \end{gathered}$$ and PY and $$\begin{gathered} ({\text{Mg}}_{{\text{1}}{\text{.98}}} {\text{Fe}}_{{\text{0}}{\text{.02}}}^{{\text{2 + }}} {\text{Mn}}_{{\text{< 0}}{\text{.001}}}^{{\text{2 + }}} {\text{Fe}}_{{\text{0}}{\text{.45}}}^{{\text{3 + }}} {\text{Ti}}_{{\text{0}}{\text{.01}}} {\text{Al}}_{{\text{3}}{\text{.56}}} )_{6.02}^{[5,6]} \hfill \\ {\text{Al}}_{{\text{2}}{\text{.00}}}^{{\text{[5]}}} [({\text{Si}}_{{\text{3}}{\text{.91}}} {\text{O}}_{{\text{17}}{\text{.73}}} {\text{(OH)}}_{{\text{2}}{\text{.27}}} ] \hfill \\ \end{gathered}$$ for GY. The Mössbauer spectra at room temperature contain one main doublet with isomer shifts and quadrupole splittings of 0.36 (PY), 0.38 (GY) and 1.00 (PY), 0.92 (GY) mm s^?1, respectively. These values correspond to Fe³⁺ in six or five-fold coordination. The doublet components have anomalously large half widths indicating either accomodation of Fe³⁺ in more than one position (e.g., octahedraA1 and five coordinatedA2) or the yet unresolved superstructure. Besides strong absorption in the ultraviolet (UV) starting from about 25,000 cm^?1, the polarized optical absorption spectra are dominated by strong bands around 16,500 and 21,000 cm^?1 (PY) and a medium strong band at around 13,800 cm^?1 (GY). Position and polarization of these bands, in combination with the UV absorption, explain the colour and pleochroism of the two varieties. The bands in question are assigned to homonuclear metal-to-metal charge transfer transitions: Mn²⁺(A1) Mn³⁺(A1′) ? Mn³⁺(A1) Mn²⁺(A1′) and Mn²⁺(A1) Mn³⁺(A2 ? Mn³⁺(A1) Mn²⁺(A2) in PY and Fe²⁺(A1) Fe³⁺(A1′) ? Fe³⁺(A1) Fe²⁺(A1′) in GY. The evidence for homonuclear Mn²⁺ Mn³⁺ charge transfer (CTF) is not quite clear and needs further study. Heteronuclear FeTi CTF does not contribute to the spectra. In PY, additional weak bands were resolved at energies around 17,700, 18,700, 21,000, and 21,900 cm^?1 and assigned to Mn³⁺ in two positions. Weak bands around 10,000 cm^?1 in both varieties are assigned to Fe²⁺ spin-alloweddd-transitions. Very weak and sharp bands, around 15,400, 16,400, 21,300, 22,100, 23,800, and 25,000 cm^?1 are identified in GY and assigned to Fe³⁺ spin-forbiddendd-transitions.     

Solubility of Ca and Al in orthopyroxene from spinel peridotite: an improved version of an empirical geothermometer

Geothermometric equations for spinel peridotites by Fujii (1976), Gasparik and Newton (1984), and Chatterjee, and Terhart (1985) based on the reaction enstatite (en)+spinel (sp)Mg–Tschermaks (mats)+forsterite (fo) were tested using a nearly isothermal suite of mantle xenoliths from the Eifel, West Germany. In spite of using activities of MgAl₂O₄, en, and mats to allow for the non-ideal solution behaviour of the constituent phases, temperatures calculated from these equations systematically change as a function of Cr/(Cr+AL+Fe³⁺) in spinel. We propose an improved version of the empirical geothermometer for spinel peridotites of Sachtleben and Seck (1981) derived from the evaluation of the solubilities of Ca and Al in orthopyroxene from more than 100 spinel peridotites from the Rhenish Volcanic Province. A least squares regression yielded a smooth correlation between

Genesis of secondary uranium minerals associated with jasperoid veins,El Erediya area,Eastern Desert,Egypt

Uranium mineralization in the El Erediya area, Egyptian Eastern Desert, has been affected by both high temperature and low temperature fluids. Mineralization is structurally controlled and is associated with jasperoid veins that are hosted by a granitic pluton. This granite exhibits extensive alteration, including silicification, argillization, sericitization, chloritization, carbonatization, and hematization. The primary uranium mineral is pitchblende, whereas uranpyrochlore, uranophane, kasolite, and an unidentified hydrated uranium niobate mineral are the most abundant secondary uranium minerals. Uranpyrochlore and the unidentified hydrated uranium niobate mineral are interpreted as alteration products of petscheckite. The chemical formula of the uranpyrochlore based upon the Electron Probe Micro Analyzer (EPMA) is . It is characterized by a relatively high Zr content (average ZrO₂ = 6.6 wt%). The average composition of the unidentified hydrated uranium niobate mineral is , where U and Nb represent the dominant cations in the U and Nb site, respectively. Uranophane is the dominant U⁶⁺ silicate phase in oxidized zones of the jasperoid veins. Kasolite is less abundant than uranophane and contains major U, Pb, and Si but only minor Ca, Fe, P, and Zr. A two-stage metallogenetic model is proposed for the alteration processes and uranium mineralization at El Erediya. The primary uranium minerals were formed during the first stage of the hydrothermal activity that formed jasperoid veins in El Eradiya granite (130–160 Ma). This stage is related to the Late Jurassic–Early Cretaceous phase of the final Pan-African tectono-thermal event in Egypt. After initial formation of El Erediya jasperoid veins, a late stage of hydrothermal alteration includes argillization, dissolution of iron-bearing sulfide minerals, formation of iron-oxy hydroxides, and corrosion of primary uranium minerals, resulting in enrichment of U, Ca, Pb, Zr, and Si. During this stage, petscheckite was altered to uranpyrochlore and oxy-petscheckite. Uranium was likely transported as uranyl carbonate and uranyl fluoride complexes. With change of temperature and pH, these complexes became unstable and combined with silica, calcium, and lead to form uranophane and kasolite. Finally, at a later stage of low-temperature supergene alteration, oxy-petscheckite was altered to an unidentified hydrated uranium niobate mineral by removal of Fe.