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.     

Re-evaluation of lead isotopic data,southern Massif Central,France

Three independent Pb isotope homogenizing processes operating on large volumes of rock material during limited intervals in the Phanerozoic have been used to define a unique evolutionary curve for rock and ore lead isotopic compositions of the southern Massif Central, France. The model is

The Reaction Titanite+Kyanite=Anorthite+Rutile and Titanite-Rutile Barometry in Eclogites

Titanite and rutile are a common mineral pair in eclogites, and many equilibria involving these phases are potentially useful in estimating pressures of metamorphism. We have reversed one such reaction,

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

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

(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 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 role of the minor substitutions in the crystal structure of natural challacolloite, KPb2Cl5, and hephaistosite, TlPb2Cl5, from Vulcano (Aeolian Archipelago, Italy)

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

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

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

Solubility of alumina in orthopyroxene plus spinel as a geobarometer in complex systems. Applications to spinel-bearing alpine-type peridotites

The addition of Fe and Cr to the simple system MgO-SiO₂-Al₂O₃ markedly affects the activities of phases involved in the equilibrium

The origin of garnet in the anorthosite-charnockite suite of the Adirondacks

Detailed analysis of textural and chemical criteria in rocks of the anorthosite-charnockite suite of the Adirondack Highlands suggests that development of garnet in silica-saturated rocks of the suite occurs according to the reaction: $$\begin{gathered} {\text{Anorthite}} {\text{Orthopyroxene}} {\text{Quartz}} \hfill \\ {\text{2CaAl}}_{\text{2}} {\text{Si}}_{\text{2}} {\text{O}}_{\text{8}} + (6 - \alpha )({\text{Fe,Mg}}){\text{SiO}}_{\text{3}} + \alpha {\text{Fe - Oxide + (}}\alpha {\text{ - 2)SiO}}_{\text{2}} \hfill \\ {\text{Garnet}} {\text{Clinopyroxene}} \hfill \\ = {\text{Ca(Fe,Mg)}}_{\text{5}} {\text{Al}}_{\text{4}} {\text{Si}}_{\text{6}} {\text{O}}_{{\text{24}}} + {\text{Ca(Fe,Mg)Si}}_{\text{2}} {\text{O}}_{\text{6}} \hfill \\ \end{gathered} $$ , where α is a function of the distribution of Fe and Mg between the several coexisting ferromagnesian phases. Depending upon the relative amounts of Fe and Mg present, quartz may be either a reactant or a product. Using an aluminum-fixed reference frame, this reaction can be restated in terms of a set of balanced partial reactions describing the processes occurring in spatially separated domains within the rock. The fact that garnet invariably replaces plagioclase as opposed to the other reactant phases indicates that the aluminum-fixed model is valid as a first approximation. This reaction is univariant and produces unzoned garnet. It differs from a similar equation proposed by de Waard (1965) for the origin of garnet in Adirondack metabasic rocks, i.e. 6 Orthopyroxene+2 Anorthite = Clinopyroxene+Garnet+2 Quartz, the principle difference being that iron oxides (ilmenite and/or magnetite) are essential reactant phases in the present reactions. The product assemblage (garnet+clinopyroxene+plagioclase ± orthopyroxene ± quartz) is characteristic of the clinopyroxene-almandine subfacies of the granulite facies.     

Thermodynamics of coexisting cummingtonite-hornblende Pairs

The activity-composition relations for calcium-rich and calcium-poor amphiboles are calculated from the composition of coexisting cummingtonite-hornblende pairs from a suite of New Zealand rhyolites. The activities are formulated in terms of site occupancies and the regular solution model is used to represent non-ideal mixing of the cations on each site. The regular solution parameters for each site are calculated from the compositions of the coexisting amphiboles. The resulting activity-composition relations are used to calibrate the standard Gibbs energy change for the reaction $${\text{7MgSiO}}_{\text{3}} {\text{ + SiO}}_{\text{2}} {\text{ + H}}_{\text{2}} {\text{O = Mg}}_{\text{7}} {\text{Si}}_{\text{8}} {\text{O}}_{{\text{22}}} {\text{(OH)}}_{\text{2}} $$ assuming that the lowest temperature rhyolites in this suite crystallised at \(P_{{\text{H}}_2 {\text{O}}} = P_{{\text{total}}} \)      

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.     

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

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

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

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

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.     

Chemistry, Raman and infrared spectroscopic characterization of the phosphate mineral reddingite: (MnFe)3(PO4)2(H2O,OH)3, a mineral found in lithium-bearing pegmatite

Detailed investigation of an intermediate member of the reddingite–phosphoferrite series, using infrared and Raman spectroscopy, scanning electron microcopy and electron microprobe analysis, has been carried out on a homogeneous sample from a lithium-bearing pegmatite named Cigana mine, near Conselheiro Pena, Minas Gerais, Brazil. The determined formula is $ ({\text{Mn}}_{1.60} {\text{Fe}}_{1.21} {\text{Ca}}_{0.01} {\text{Mg}}_{0.01} )_{\sum 2.83} ({\text{PO}}_{4} )_{2.12} \cdot ({\text{H}}_{2} {\text{O}}_{2.85} {\text{F}}_{0.01} )_{\sum 2.86} $ , indicating predominance in the reddingite member. Raman spectroscopy coupled with infrared spectroscopy supports the concept of phosphate, hydrogen phosphate and dihydrogen phosphate units in the structure of reddingite-phosphoferrite. Infrared and Raman bands attributed to water and hydroxyl stretching modes are identified. Vibrational spectroscopy adds useful information to the molecular structure of reddingite–phosphoferrite.     

Weathering Processes in the Min Jiang: Major Elements, {}^{87}{\text{Sr/}}{}^{86}{\text{Sr}},\;\delta {}^{34}{\text{S}}_{{\text{SO}}_{\text{4}} } ,\;{\text{and}}\;\delta {}^{18}{\text{O}}_{{\text{SO}}_{\text{4}} }

We investigated the dissolved major elements, $ {}^{87}{\text{Sr/}}{}^{86}{\text{Sr}},\;\delta {}^{34}{\text{S}}_{{\text{SO}}_{\text{4}} } ,\;{\text{and}}\;\delta {}^{18}{\text{O}}_{{\text{SO}}_{\text{4}} } $ composition of the Min Jiang, a headwater tributary of the Chang Jiang (Yangtze River). A forward calculation method was applied to quantify the relative contribution to the dissolved load from rain, evaporite, carbonate, and silicate reservoirs. Input from carbonate weathering dominated the major element composition (58–93%) and that from silicate weathering ranged from 2 to 18% in unperturbed Min Jiang watersheds. Most samples were supersaturated with respect to calcite, and the CO₂ partial pressures were similar to or up to ～5 times higher than atmospheric levels. The Sr concentrations in our samples were low (1.3–2.5 μM) with isotopic composition ranging from 0.7108 to 0.7127, suggesting some contribution from felsic silicates. The Si/(Na^* + K) ratios ranged from 0.5 to 2.5, which indicate low to moderate silicate weathering intensity. The $ \delta {}^{34}{\text{S}}_{{\text{SO}}_{\text{4}} } \;{\text{and}}\;\delta {}^{18}{\text{O}}_{{\text{SO}}_{\text{4}} } $ for five select samples showed that the source of dissolved sulfate was combustion of locally consumed coal. The silicate weathering rates were 23–181 × 10³ mol/km²/year, and the CO₂ consumption rates were 31–246 × 10³ mol/km²/year, which are moderate on a global basis. Upon testing various climatic and geomorphic factors for correlation with the CO₂ consumption rate, the best correlation coefficients found were with water temperature (r ² = 0.284, p = 0.009), water discharge (r ² = 0.253, p = 0.014), and relief (r ² = 0.230, p = 0.019).     

Thermodynamic properties of almandine-grossular garnet solid solutions

The mixing properties of Fe₃Al₂Si₃O₁₂-Ca₃Al₂Si₃O₁₂ garnet solid solutions have been studied in the temperature range 850–1100° C. The experimental method involves measuring the composition of garnet in equilibrium with an assemblage in which the activity of the Ca₃Al₂Si₃O₁₂ component is fixed. Experiments on the assemblage garnet solid solution, anorthite, Al₂SiO₅ polymorph and quartz at known pressure and temperature fix the activity of the Ca₃Al₂Si₃O₁₂ component through the equilibrium: 1 $$\begin{gathered} {\text{3CaAl}}_{\text{2}} {\text{Si}}_{\text{2}} {\text{O}}_{\text{8}} \rightleftarrows {\text{Ca}}_{\text{3}} {\text{Al}}_{\text{2}} {\text{Si}}_{\text{3}} {\text{O}}_{{\text{12}}} \hfill \\ {\text{Anorthite garnet}} \hfill \\ {\text{ + 2Al}}_{\text{2}} {\text{SiO}}_{\text{5}} {\text{ + SiO}}_{\text{2}} \hfill \\ {\text{ sillimanite/kyanite quartz}}{\text{.}} \hfill \\ \end{gathered}$$ This equilibrium, with either sillimanite or kyanite as the aluminosilicate mineral, was used to control \({\text{a}}_{{\text{Ca}}_{\text{3}} {\text{Al}}_{\text{2}} {\text{Si}}_{\text{3}} {\text{O}}_{{\text{12}}} }^{{\text{gt}}} \) . The compositions of the garnet solutions produced were determined by measurement of their unit cell edges. At 1 bar Fe₃Al₂Si₃O₁₂-Ca₃Al₂Si₃O₁₂ garnets exhibit negative deviations from ideality at the Fe-rich end of the series and positive deviations at the calcium end. With increasing pressure the activity coefficients for the Ca₃Al₂Si₃O₁₂ component increase because the partial molar volume of this component is greater than the molar volume of pure grossular. Previous studies indicate that the activity coefficients for the Ca₃Al₂Si₃O₁₂ component also increase with increasing (Mg/Mg+Fe) ratio of the garnet. The region of negative deviation from ideality implies a tendency towards formation of a stable Fe-Ca garnet component. Evidence in support of this conclusion has been found in a natural Fe-rich garnet which was found to contain two different garnet phases of distinctly different compositions.     

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.     

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

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