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

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.     

The reaction garnet+clinopyroxene+quartz =2 orthopyroxene+anorthite: A potential geobarometer for granulites

A mineralogic geobarometer based on the reaction garnet+clinopyroxene+quartz=2 orthopyroxene+anorthite is proposed. The geobarometric formulations for the Fe- and Mg- end member equilibria are $$\begin{gathered} P_{({\text{Fe}})} {\text{ }}({\text{bars}}){\text{ = 32}}{\text{.097 }}T{\text{ }} - {\text{ 26385 }} - {\text{ 22}}{\text{.79 (}}T - 848 - T1{\text{n(}}T/848{\text{))}} \hfill \\ {\text{ }} - (3.655 + 0.0138T){\text{ }}\left( {\frac{{{\text{(}}T - 848{\text{)}}^{\text{2}} }}{T}} \right) \hfill \\ {\text{ }} - {\text{(3}}{\text{.123) }}T1{\text{n }}\frac{{(a_{a{\text{n}}}^{{\text{Plag}}} )(a_{{\text{fs}}}^{{\text{P}}\ddot u{\text{x}}} )^2 }}{{(a_{{\text{alm}}}^{{\text{Gt}}} )(a_{{\text{hed}}}^{{\text{Opx}}} )}} \hfill \\ P_{({\text{Mg}})} {\text{ (bars) = 9}}{\text{.270 }}T + 4006 - 0.9305{\text{ }}(T - 848 - T1{\text{n (}}T/848{\text{)}}) \hfill \\ {\text{ }} - (1.1963{\text{ }} - {\text{ }}6.0128{\text{ x 10}}^{ - {\text{3}}} T)\left( {\frac{{(T - 848)^2 }}{T}} \right) \hfill \\ {\text{ }} - 3.489{\text{ }}T1{\text{n }}\frac{{(a_{an}^{{\text{Plag}}} ){\text{ }}(a_{{\text{ens}}}^{{\text{Opx}}} )}}{{{\text{(}}a_{{\text{pyr}}}^{{\text{Gt}}} {\text{) (}}a_{{\text{diop}}}^{{\text{Cpx}}} {\text{)}}}}. \hfill \\ \end{gathered}$$ The end member thermodynamic data have been taken from the data base of Helgeson et al. (1978) and Saxena and Erikson (1983). The activities of pyroxene components and anorthite in plagioclase have been modelled after Wood and Banno (1973) and Newton (1983) respectively. The activities of pyrope and almandine are calculated from the binary interaction parameters for garnet solid solutions proposed by Saxena and Erikson (1983). Pressures computed from these equations for fifty sets of published mineral data from several granulite areas are comparable with those obtained from dependable geobarometers. The pressure values determined from the Fe-end member equilibrium appear to be more reasonable than those from the Mg-end member reaction. It is likely that the difference in pressures computed from the Fe- and Mg-end members, ΔP ^*, have been caused by non-ideal mixing in the phases, especially in garnets.     

A thermodynamic formulation of hydrous cordierite

A thermodynamic formulation of hydrous Mg-cordierite (Mg₂Al₄Si₅O₁₈·nH₂O) has been obtained by application of calorimetric and X-ray diffraction data for hydrous cordierite to the results of hydrothermal syntheses. The data include measurements of the molar heat capacity and enthalpy of hydration and the molar volume. The synthesis data are consistent with a thermodynamic formulation in which H₂O mixes ideally on a single crystallographic site in hydrous cordierite. The standard molar Gibbs free energy of hydration is-9.5±1.0 kJ/mol (an average of 61 syntheses). The standard molar entropy of hydration derived from this value is-108±3 J/mol-K. An equation providing the H₂O content of cordierite as a function of temperature and fugacity of H₂O is as follows (n moles of H₂O per formula unit, n<1): $$\begin{gathered}n = {{f_{{\text{ H}}_{\text{2}} O}^{\text{V}} } \mathord{\left/{\vphantom {{f_{{\text{ H}}_{\text{2}} O}^{\text{V}} } {\left( {f_{{\text{ H}}_{\text{2}} O}^{\text{V}} + {\text{exp}}\left[ { - {\text{3}}{\text{.8389}} - 5025.2\left( {\frac{1}{T} - \frac{1}{{298.15}}} \right)} \right.} \right.}}} \right.\kern-\nulldelimiterspace} {\left( {f_{{\text{ H}}_{\text{2}} O}^{\text{V}} + {\text{exp}}\left[ { - {\text{3}}{\text{.8389}} - 5025.2\left( {\frac{1}{T} - \frac{1}{{298.15}}} \right)} \right.} \right.}} \hfill \\{\text{ }}\left. {\left. { - {\text{ln}}\left( {\frac{T}{{{\text{298}}{\text{.15}}}}} \right) - \left( {\frac{{298.15}}{T} - 1} \right)} \right]} \right) \hfill \\\end{gathered}$$ Application of this formulation to the breakdown reaction of Mg-cordierite to an assemblage of pyrope-sillimanite-quartz±H₂O shows that cordierite is stabilized by 3 to 3.5 kbar under H₂O-saturated conditions. The thermodynamic properties of H₂O in cordierite are similar to those of liquid water, with a standard molar enthalpy and Gibbs free energy of hydration that are the same (within experimental uncertainty) as the enthalpy and Gibbs free energy of vaporization. By contrast, most zeolites have Gibbs free energies of hydration two to four times more negative than the corresponding value for the vaporization of water.     

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.     

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.     

Multiple-pumped-well aquifer test to determine the anisotropic properties of a karst limestone aquifer in Pasco County, Florida, USA

Groundwater-level data from an aquifer test utilizing four pumped wells conducted in the South Pasco wellfield in Pasco County, Florida, USA, were analyzed to determine the anisotropic transmissivity tensor, storativity, and leakance in the vicinity of the wellfield. A weighted least-squares procedure was used to analyze drawdowns measured at eight observation wells, and it was determined that the major axis of transmissivity extends approximately from north to south and the minor axis extends approximately from west to east with an angle of anisotropy equal to N4.54°W. The transmissivity along the major axis ${\left( {T_{{\xi \xi }} } \right)}$ is 14,019 m² day^–1, and the transmissivity along the minor axis ${\left( {T_{{\eta \eta }} } \right)}$ is 4,303 m² day^–1. The equivalent transmissivity $T_{e} = {\left( {T_{{\xi \xi }} T_{{\eta \eta }} } \right)}^{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-0em} 2}} = 7,767{{\text{m}}^{2} } \mathord{\left/ {\vphantom {{{\text{m}}^{2} } {{\text{day}}^{{ - {\text{1}}}} }}} \right. \kern-0em} {{\text{day}}^{{ - {\text{1}}}} }$ , and the ratio of anisotropy is 3.26. The storativity of the aquifer is 7.52?×?10^?4, and the leakance of the overlying confining unit is 1.37?×?10^?4 day^?1. The anisotropic properties determined for the South Pasco wellfield in this investigation confirm the results of previous aquifer tests conducted in the wellfield and help to quantify the NW–SE to NE–SW trends for regional fracture patterns and inferred solution-enhanced flow zones in west-central Florida.     

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

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:

Statistical mechanics of coupled solid solutions in the dilute limit

A statistical mechanical analysis of the limiting laws for coupled solid solutions shows that the random model, in which the configurational entropy is calculated as if atoms mix randomly on each crystallographic site, is correct as a first approximation. In coupled solid solutions, since atoms of different valence substitute on the same sites, significant short-range order which reduces the entropy can be expected. A first-order correction is rigorously obtained for the entropy in dilute binary short-range ordered coupled solid solutions: $$\bar S^{{\text{XS}}} {\text{/R = }}Q\left( {{\text{e}}^{--H_{\text{A}} /{\text{R}}T} \left( {\frac{{H_{\text{A}} }}{{{\text{R}}T}} + 1} \right) - 1} \right)N_2^a N_4^b ,$$ where Q is the number of positions an associated cation pair can assume per formula unit, H _A is the association energy per formula unit, and N ₂ ^a and N ₄ ^b are the site occupancy fractions for atoms 2 and 4 that are dilute on sites a and b. S ^XS is the configurational entropy minus the random model entropy. Aluminous pyroxenes on the joints diopside-jadeite and diopside-CaTs are examined as examples. A generalization for dilute multiple component solutions, including possible long-range ordering variations is given by: $$\frac{{\bar S^{{\text{XS}}} }}{{\text{R}}}{\text{ = }}\sum\limits_i {\sum\limits_j {\sum\limits_k {Q_i } } \left( {{\text{e}}^{--H_{\text{A}}^{j{\text{ }}k{\text{, }}i} /{\text{R}}T} \left( {\frac{{H_{\text{A}}^{j{\text{ }}k{\text{, }}i} }}{{{\text{R}}T}} + 1} \right) - 1} \right)N_j^l N_k^m ,} $$ where i labels each crystallographically distinct pair, j and k label atomic species, l and m label crystallographic sites, and the N's are site occupancy fractions for the solute atoms. A total association model is examined as well as the partial association and random models. Real solution behavior must lie between the total association model and the random model. Molecular models in which the ideal activity is proportional to a mole fraction, which in itself is not always unambiguously defined, do not lie in this range and furthermore have no physical justification.     

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.     

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.     

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.     

Assemblages among tephroite,pyroxmangite, rhodochrosite,quartz: Experimental data and occurrences in the Rhetic Alps

Reactions involving the phases quartz-rhodochrosite-tephroite-pyroxmangite-fluid have been studied experimentally in the system MnO-SiO₂-CO₂-H₂O at a pressure of 2 000 bars and resulted in the following expressions 1 $$\begin{gathered} {\text{Rhodochrosite + Quartz = Pyroxmangite + CO}}_2 \hfill \\ {\text{ log}}_{{\text{10}}} K^{{\text{2000 bars}}} = - \frac{{11.765}}{T} + 18.618. \hfill \\ {\text{Rhodochrosite + Pyroxmangite = Tephroite + CO}}_2 \hfill \\ {\text{ log}}_{{\text{10}}} K^{{\text{2000 bars}}} = - \frac{{7.083}}{T} + 11.870. \hfill \\ \end{gathered}$$ which can be used to derive data for the remaining two reactions among the phases under consideration. Field data from the Alps are in agreement with the metamorphic sequence resulting from the experiments.     

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 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:

Zinc isotope fractionation under vaporization processes and in aqueous solutions

Equilibrium Zn isotope fractionation was investigated using first-principles quantum chemistry methods at the B3LYP/6-311G* level. The volume variable cluster model method was used to calculate isotope fractionation factors of sphalerite, smithsonite, calcite, anorthite, forsterite, and enstatite. The water-droplet method was used to calculate Zn isotope fractionation factors of Zn²⁺-bearing aqueous species; their reduced partition function ratio factors decreased in the order \(\left[ {{\text{Zn}}\left( {{\text{H}}_{2} {\text{O}}} \right)_{6} } \right]^{2 + } > \left[ {{\text{ZnCl}}\left( {{\text{H}}_{2} {\text{O}}} \right)_{5} } \right]^{ + } > \left[ {{\text{ZnCl}}_{2} \left( {{\text{H}}_{2} {\text{O}}} \right)_{4} } \right] > \left[ {{\text{ZnCl}}_{3} \left( {{\text{H}}_{2} {\text{O}}} \right)_{2} } \right]^{ - } > {\text{ZnCl}}_{4} ]^{2 - }\). Gaseous ZnCl₂ was also calculated for vaporization processes. Kinetic isotope fractionation of diffusional processes in a vacuum was directly calculated using formulas provided by Richter and co-workers. Our calculations show that in addition to the kinetic isotope effect of diffusional processes, equilibrium isotope fractionation also contributed nontrivially to observed Zn isotope fractionation of vaporization processes. The calculated net Zn isotope fractionation of vaporization processes was 7–7.5‰, with ZnCl₂ as the gaseous species. This matches experimental observations of the range of Zn isotope distribution of lunar samples. Therefore, vaporization processes may be the cause of the large distribution of Zn isotope signals found on the Moon. However, we cannot further distinguish the origin of such vaporization processes; it might be due either to igneous rock melting in meteorite bombardments or to a giant impact event. Furthermore, isotope fractionation between Zn-bearing aqueous species and minerals that we have provided helps explain Zn isotope data in the fields of ore deposits and petrology.     

Enstatite-diopside solvus and geothermometry

The enstatite-diopside solvus presents certain interesting thermodynamic and crystal-structural problems. The solvus may be considered as parts of two solvi one with the ortho-structure and the other with clino-structure. By assuming the standard free energy change for the two reactions (MgMgSi₂O₆)opx ? (MgMgSi₂O₆)cpx and (CaMgSi₂O₆) opx ? (CaMgSi₂O₆) cpx as 500 and 1 000 to 3 000 cal/mol respectively, it is possible to calculate the regular solution parameter W for orthopyroxene and clinopyroxene. These W's essentially refer to mixing on M2 sites. The expression for the equilibrium constant by assuming ideal mixing for Fe-Mg, Fe-Ca and non-ideal mixing for Ca-Mg on binary M1 and ternary M2 sites is given by 1 $$K_a = \frac{{X_{{\text{Mg - cpx}}}^{{\text{M1}}} X_{{\text{Mg - cpx}}}^{{\text{M2}}} \exp \left[ {\frac{{W_{{\text{cpx}}} }}{{RT}}\left\{ {X_{{\text{Ca - cpx}}}^{{\text{M2}}} \left( {X_{{\text{Ca - cpx}}}^{{\text{M2}}} + X_{{\text{Fe - cpx}}}^{{\text{M2}}} } \right)} \right\}} \right]}}{{X_{{\text{Mg - cpx}}}^{{\text{M1}}} X_{{\text{Mg - opx}}}^{{\text{M2}}} \exp \left[ {\frac{{W_{{\text{cpx}}} }}{{RT}}\left\{ {X_{{\text{Ca - opx}}}^{{\text{M2}}} \left( {X_{{\text{Ca - opx}}}^{{\text{M2}}} + X_{{\text{Fe - opx}}}^{{\text{M2}}} } \right)} \right\}} \right]}}$$ where X's are site occupancies, R is 1.987 and T is temperature in ^oK. Temperature of pyroxene crystallization may be estimated by substituting for T in the above equation until the equation ?RT In K _a=500 is satisfied. The shortcomings of this method are the incomplete standard free energy data on the end member components and the absence of site occupancy data in pyroxenes at high temperatures. The assumed free energy data do, however, show the possible extent of inaccuracy in temperature estimates resulting from the neglect of Mg-Ca non ideality.