Reversed experimental calibration of the garnet-clinopyroxene Fe — Mg exchange thermometer

The temperature-sensitive Fe,Mg exchange equilibrium,

Transition enthalpies and entropies of high pressure zinc metasilicates and zinc metagermanates

The enthalpies of transition at T= 298 K between zinc metasilicate assemblages, measured by molten oxide solution calorimetry, are:

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.     

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.     

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 carbon dioxide solubility in alkali basalts: an experimental study

Experiments were conducted to determine CO₂ solubilities in alkali basalts from Vesuvius, Etna and Stromboli volcanoes. The basaltic melts were equilibrated with nearly pure CO₂ at 1,200°C under oxidizing conditions and at pressures ranging from 269 to 2,060 bars. CO₂ solubility was determined by FTIR measurements. The results show that alkalis have a strong effect on the CO₂ solubility and confirm and refine the relationship between the compositional parameter Π devised by Dixon (Am Mineral 82:368–378, 1997) and the CO₂ solubility. A general thermodynamic model for CO₂ solubility in basaltic melts is defined for pressures up to 2 kbars. Based on the assumption that O²⁻ and CO₃²⁻ mix ideally, we have:

The crystal structure of a triclinic roscherite

Summary The crystallography of roscherite is more complicated than previously thought. Single crystal X-ray work on material from Foote Mine (California) gave triclinic symmetry. The unit cell corresponding to the one adopted for monoclinic roscherite hasa=15.921,b=11.965,c=6.741 Å, =91°04, =94°21, =89°59 1/2, space group . The least-squares refinement of the structure using 2380 non zero reflections with anisotropic temperature factors resulted in a conventional reliability factorR=0.060.The X-ray study indicates the formula while that proposed for monoclinic roscherite is The atomic arrangements of both varieties of roscherite are very similar. The lowering of symmetry is caused by the segregation of the trivalent cations into only half of the sites of a monoclinic point position. Crystallochemical considerations suggest that the symmetry of roscherite does not depend on the kind of trivalent cations occupying the 6-coordinated position, but rather by the ratio between trivalent and divalent metal ions.

---

Die Kristallstruktur eines triklinen Roscherites

Zusammenfassung Die Kristallographie des Roscherites ist komplizierter als man bisher annahm. Einkristall-Röntgenuntersuchungen an Material von Foote Mine (Kalifornien) ergaben trikline Symmetrie. Die Elementarzelle, welche der für monoklinen Roscherit angenommenen entspricht, hata=15,921,b=11,965,c=6,741 Å, =91°04, =94°21, =89°59 1/2, Raumgruppe . Die Verfeinerung der Struktur mit der Methode der kleinsten Quadrate ergab unter Verwendung anisotroper Temperaturfaktoren für 2380 beobachtete Reflexe einen konventionellen ZuverlässigkeitsindexR=0,060.Die Röntgenuntersuchung weist auf die Formel , während für monoklinen Roscherit vorgeschlagen wurde. Die Atomanordnungen beider Abarten des Roscherites sind sehr ähnlich. die Symmetrieerniedrigung wird dadurch hervorgerufen, daß die dreiwertigen Kationen nur die Hälfte der Positionen einer monoklinen Punktlage besetzen. Kristallchemische Überlegungen weisen darauf hin, daß die Symmetrie nicht von der Art der dreiwertigen Kationen, welche eine 6-koordinierte Punktlage besetzen, abhängt, sondern vielmehr von dem Mengenverhältnis zwischen 3-wertigen und 2-wertigen Metallionen.

---

---

With 1 Figure     

Dislocation dissociation in CaGeO3 perovskite

Dissociated dislocations have been observed for the first time by transmission electron microscopy in the perovskite-structure compound CaGeO₃. Dislocations with Burgers vectors \(\left[ {1\bar 10} \right]\) and [001] (in pseudo-cubic index) are dissociated into collinear partials on the (110) plane: $$\left[ {1\bar 10} \right] = {1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-\nulldelimiterspace} 2}\left[ {1\bar 10} \right] + {1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-\nulldelimiterspace} 2}\left[ {1\bar 10} \right]$$ and [001] = 1/2[001] + 1/2[001]. The partials react to form octagonal extended nodes. The stacking fault ribbons with displacement vector \(\left[ {1\bar 10} \right]\) have a width of 350 A, which corresponds to a stacking fault energy of 35 erg/cm² (or mJ/m²).     

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

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

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

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.     

Thermochemical and pressure-volume-temperature systematics of data on solids,examples: Tungsten and MgO

Data systematization using the constraints from the equation $$Cp = Cv + \alpha _P {}^2V_T K_T T$$ where C _p, C _v, α _p, K _T and V are respectively heat capacity at constant pressure, heat capacity at constant volume, isobaric thermal expansion, isothermal bulk modulus and molar volume, has been performed for tungsten and MgO. The data are $$K_T (W) = 1E - 5/(3.1575E - 12 + 1.6E - 16T + 3.1E - 20T^2 )$$ $$\alpha _P (W) = 9.386E - 6 + 5.51E - 9T$$ $$C_P (W) = 24.1 + 3.872E - 3T - 12.42E - 7T^2 + 63.96E - 11T^3 - 89000T^{ - 2} $$ $$K_T (MgO) = 1/(0.59506E - 6 + 0.82334E - 10T + 0.32639E - 13T^2 + 0.10179E - 17T^3 $$ $$\alpha _P (MgO) = 0.3754E - 4 + 0.7907E - 8T - 0.7836/T^2 + 0.9148/T^3 $$ $$C_P (MgO) = 43.65 + 0.54303E - 2T - 0.16692E7T^{ - 2} + 0.32903E4T^{ - 1} - 5.34791E - 8T^2 $$ For the calculation of pressure-volume-temperature relation, a high temperature form of the Birch-Murnaghan equation is proposed $$P = 3K_T (1 + 2f)^{5/2} (1 + 2\xi f)$$ Where $$K_T = 1/(b_0 + b_1 T + b_2 T^2 + b_3 T^3 )$$ $$f = (1/2)\{ [V(1,T)/V(P,T)]^{2/3} - 1\} $$ $$\xi = ({3 \mathord{\left/ {\vphantom {3 4}} \right. \kern-\nulldelimiterspace} 4})[K'_0 + K'_1 \ln ({T \mathord{\left/ {\vphantom {T {300}}} \right. \kern-\nulldelimiterspace} {300}}) - 4]$$ where in turn $$V(1,T) = V_0 [\exp (\int\limits_{300}^T {\alpha dT)]} $$ . The temperature dependence of the pressure derivative of the bulk modulus (K′₁) is estimated by using the shock-wave data. For tungsten the data are K′₀ = 3.5434, K′₁ = 0.032; for MgO K′₀ = 4.17 and K′₁ = 0.1667. For calculating the Gibbs free energy of a solid at high pressure and at temperatures beyond that of melting at 1 atmosphere, it is necessary to define a high-temperature reference state for the fictive solid.     

Stochastic analysis of bed-thickness distribution of sediments

On formation of a bed and distribution of bed thickness, A. N. Kolmogorov presented a mathematical explanation that if repetitive alternations of material accumulation and erosion form a sequence of beds, the resultant bed-thickness distribution curve takes a shape truncated by the ordinate at zero thickness. In this truncated distribution curve, its continuation and extension from positive to negative thickness represents the distribution of beds with negative thickness, that is, the depth of erosion. When a distribution curve, including both positive and negative parts, is expressed by a function f(x),the ratio \(\int_0^\infty {f(x)dx to} \int_{ - \infty }^\infty {f(x)dx} \) ,called Kolmogorov's coefficient and designated as p,is a parameter representing the degree of accumulation in the depositional environment. On the assumption that f(x)is described by the Gaussian distribution function, the coefficient pfor Permian and Pliocene sequences in central Japan was calculated. The coefficients also were obtained from published data for different types of sediments from other areas. It was determined that they are more or less different depending on their depositional environments. The calculated results are summarized as follows: $$\begin{gathered} p = 0.80 - 1.0for{\text{ }}alluvial{\text{ }}or{\text{ }}fluvial{\text{ }}deposits \hfill \\ p = 0.65 - 0.95for{\text{ }}nearshore{\text{ }}sediments \hfill \\ p = 0.55 - 0.95for{\text{ }}geosynclinal{\text{ }}sediments \hfill \\ p = 0.90 - 1.0for{\text{ }}varves \hfill \\ \end{gathered} $$ In addition, a ratio \(q = \int_0^\infty {xf(x)dx/} \int_{ - \infty }^\infty {|x|f(x)dx} \) ,called Kolmogorov's ratio in this paper, is introduced for estimating a degree of total thickness actually observed in the field relative to total thickness once present in a basin. The calculated results of Kolmogorov's ratio are as follows: $$\begin{gathered} q = 0.88 - 1.0for{\text{ }}alluvial{\text{ }}or{\text{ }}fluvial{\text{ }}deposits \hfill \\ q = 0.68 - 0.98for{\text{ }}nearshore{\text{ }}sediments \hfill \\ q = 0.55 - 0.96for{\text{ }}geosynclinal{\text{ }}sediments \hfill \\ q = 0.92 - 1.0for{\text{ }}varves \hfill \\ \end{gathered} $$ The sedimentological significance of these values is discussed.     

Spinel-pyroxene-garnet relationships and their dependence on Cr/Al ratio

The partitioning of Cr and Al between coexisting spinel and clinopyroxene and the dependence of spinel-cpxgarnet equilibria on Cr/Al ratio have been investigated by a combination of phase equilibrium experiments, high temperature solution calorimetry and thermodynamic calculations.The exchange equilibrium: has a measured enthalpy change for pure phases of –2,100±500 cal at 970 K and 1 atm. Experimental reversals of Cr-Al partitioning between the spinel and clinopyroxene phases yield the following partitioning relationship: where X _i ^j refers to atomic fraction of i in the octahedral sites of phase j. The compositional dependence of partitioning implies that Al-Cr mixing in spinel is nonideal with, on the symmetrical model, a W _Cr-Al ^Sp of 2,700±500 cal/gm. atom. In contrast, aluminum-chromium mixing in clinopyroxene is close to ideal.The measured stability field of knorringite (Mg₃Cr₂Si₂O₁₂) and mixing properties of garnet have been used in conjunction with our experimental data to calculate the influence of Cr/Al ratio on the important reaction: orthopyroxene+clinopyroxene+spinel=olivine+garnetThe stability field of spinel lherzolite increases by about 2.8 Kb for every increase of 0.1 in Cr/(Cr+Al) ratio up to Cr/(Cr+Al) of 0.7. The calculated stabilization is in very good agreement with the experimental results of O'Neill (1981). The partitioning relationships are such that, at the low ratios of Cr/Al (0.07) of primitive lherzolite, clinopyroxene buffers spinel composition and sharpens the spinelgarnet reaction interval from 10 Kb (little or no clinopyroxene) down to about 2 Kb in pyroxene-rich pyrolite.     

Intrinsic oxygen Fugacity measurements on chromites from the Bushveld Complex and their petrogenetic significance

Oxygen Fugacity measurements were carried out on chromites from the Eastern Bushveld Complex (Maandagshoek) and are compared with former measurements on chromites from the western Bushveld Complex (Zwartkop Chrome Mine). These results together with those of Hill and Roeder (1974) yield the following conditions of formation for the massive chromitite layers: Western Bushveld Complex (Zwartkop Chrome Mine) $$\begin{gathered} Layer{\text{ }}T(^\circ C) p_{O_2 } (atm) \hfill \\ LG3{\text{ 1160}} - {\text{1234 10}}^{ - {\text{5}}} - 10^{ - 7.6} \hfill \\ LG4{\text{ 1175}} - {\text{1200 10}}^{ - 6.35} - 10^{ - 7.20} \hfill \\ LG6{\text{ 1162}} - {\text{1207 10}}^{ - 6.20} - 10^{ - 7.50} \hfill \\ \hfill \\ \end{gathered} $$ Eastern Bushveld Complex (Farm Maandagshoek) $$\begin{gathered} {\text{LXI 1115}} - {\text{1150 10}}^{ - 7.80} - 10^{ - 8.80} \hfill \\ ( = {\text{Steelpoort Seam)}} \hfill \\ {\text{LX 1125 10}}^{ - 8.25} \hfill \\ {\text{V 1120 10}}^{ - 8.55} \hfill \\ {\text{LII 1120 10}}^{ - 8.0} - 10^{ - 8.60} \hfill \\ \end{gathered} $$ The comparison of the data shows, that the chronitite layers within each particular sequence were formed under approximately identicalp _{o 2}- andT-conditions. The chromites from the western Bushveld Complex, however, were formed at higher temperatures and higher oxygen fugacities than the chromites from the eastern Bushveld Complex. Fromp _{o 2}-T-curves of disseminated chromites and the temperatures derived above, the following conditions of formation for the host rocks were obtained: Western Bushveld Complex $$T = 1200^\circ {\text{C; }}p_{{\text{o}}_{\text{2}} } = 10^{ - 7.25} - 10^{ - 7.50} $$ Eastern Bushveld Complex $$T = 1125^\circ {\text{C; }}p_{{\text{o}}_{\text{2}} } = 10^{ - 8.50} - 10^{ - 9.25} $$ Consequently, the host rocks in the Zwartkop-Chrome-Mine, were formed under higher temperatures and higher oxygen fugacities than the host rocks at Maandagshoek. The rock sequence in the Zwartkop-Chrome-Mine therefore originated in an earlier stage of the differentiation of the Bushveld magma. Comparison of the chromites from the host rocks with the chromites from massive layers supports Ulmer's (1969) thesis that an increase of the oxygen fugacity is responsible for the formation of massive chromitite layers. The values in this investigation show that increases of only about 0.5–1.0 log units are necessary to enhance chromitite layer formation.     

Experimental study of Fe-Mg distribution between biotite and orthopyroxene at P=490 MPa

Partitioning of Mg and Fe between coexisting biotite and orthopyroxene has been experimentally determined at temperatures 700, 750 and 800° C and 490 MPa total pressure in the system KAlO₂-MgO-FeO-SiO₂-H₂O. Oxygen fugacity was controlled by the QFM buffer. Starting materials were synthetic minerals of differing Fe/(Fe+Mg) values. Run products were analyzed for partitioning of components by a microprobe. Orthopyroxene was established to be notably inhomogeneous, whereas biotite was essentially homogeneous. To establish equilibrium relations, statistical treatment of the results of each experiment in addition to the whole complex of experimental data was applied. The regression equations for isotherms of the Fe-Mg partitioning between the minerals studied have been obtained. As a result, the equation for a two-dimensional regression may be written as: $$\begin{gathered} Y = (A + A_1 t + A_2 t^2 )(X - X^4 ) + (B + B_1 t + B_1 t^2 )(X^2 - X^4 ) + \hfill \\ (C + C_1 t + C_1 t^2 )(X^3 - X^4 ) + X^4 {\text{ where }}Y = X_{{\text{Opx}}}^{{\text{Fe}}} ;{\text{ X}} = {\text{X}}_{{\text{Bi}}}^{{\text{Fe}}} ; \hfill \\ t = 1000/T,K, \hfill \\ \begin{array}{*{20}c} {A = {\text{ }}4.59398,} & {A_1 = - {\text{ }}8.29838,} & {A_2 = {\text{ }}4.97316,} \\ {B = - 11.13731,} & {B_1 = {\text{ }}28.19304,} & {B_2 = - 20.98240,} \\ {A = {\text{ }}8.25072,} & {C_1 = - 20.80485,} & {C_2 = {\text{ }}15.35967} \\ \end{array} \hfill \\ {\text{ }}\sigma = 0.0143{\text{ }} \hfill \\ \end{gathered}$$ . This equation enables extrapolation of partitioning isotherms over a wide range of temperatures.     

From the simplest chemical system to the natural one: garnet peridotite barometry

The available experimental data on garnet-bearing-assemblages for synthetic chemical systems (MAS, FMAS, CMAS) have been used to calibrate consistent models for the Al-solubility in orthopyroxene coexisting with garnet, on the basis of equilibrium reaction Py(opx) ? Py(gt). The alternative reaction En(opx)+MgTs(opx) ? Py(gt) is discarded as it yields larger a-posteriori uncertainties. To provide a reliable equation, directly applicable to natural garnet lherzolites, each successive synthetic-system calibration is tested against Mori and Green's (1978) natural-system reequilibration data. For the MAS system, an ideal solution model with constant ΔH°, ΔV° and ΔS° based on 12-oxygen structural formulae for aluminous pyroxenes yields the best fit (GPa, K), $${\text{25,134 + 9,941 }}P - 23.177{\text{ }}T{\text{ + }}RT{\text{ ln (}}X_{{\text{Al}}}^{TB'} {\text{) = 0}}$$ . The MAS synthetic-system calibration can be directly applied to the FMAS system by adding an empirical correction term (20,835 [X _Fe ^gt ]²) independent of either pressure and temperature. However, this correction term is not important because of the limited Fe content of mantle peridotites. When calcium is added to the MAS system, the equilibrium constant is calculated as: $$K_{{\text{CMAS}}} = {{[(1 - X_{{\text{Ca}}}^{M2} )^2 (X_{{\text{Al}}}^{TB'} )]} \mathord{\left/ {\vphantom {{[(1 - X_{{\text{Ca}}}^{M2} )^2 (X_{{\text{Al}}}^{TB'} )]} {[(1 - X_{{\text{Ca}}}^X )^3 (X_{{\text{Al}}}^Y )^2 ]}}} \right. \kern-\nulldelimiterspace} {[(1 - X_{{\text{Ca}}}^X )^3 (X_{{\text{Al}}}^Y )^2 ]}}$$ where M2 and TB′ are pyroxene sites and X and Y are garnet sites. Up to 5 GPa, X _Ca ^X ～ and the CMAS experimental data agree well with the MAS model, but for Yamada and Takahashi's (1983) higher pressure experiments (up to 10 GPa), this no longer holds. Indeed, the garnet solid solution does not behave ideally and an asymmetric regular solution model is needed for application to the deepest natural samples available (>7GPa). Calibration based on new high pressure data yields, $$\begin{gathered} \Delta G_{{\text{CMAS}}}^{XS} = (X_{{\text{Ca}}}^X )(1 - X_{{\text{Ca}}}^X )(0.147 - X_{{\text{Ca}}}^X ) \hfill \\ {\text{ }} \cdot {\text{(6,440,535 - 1,490,654 }}P{\text{)}} \hfill \\ \end{gathered}$$ . According to tests of the inferred solution model, the CFMAS system is a good analogue of natural systems in the pressure, temperature and composition ranges covered by the natural-system reequilibration data (up to 1,500° C and 4 GPa). Simultaneous application of this thermobarometer and of the two-pyroxene mutual solubility thermometer (Bertrand and Mercier 1985) to the phases of the garnet-peridotite xenoliths from Thaba Putsoa, Lesotho, yields a refined paleogeotherm for southern Africa strongly contrasting with previous results. The “granular” nodules yield a thermal gradient of about 8 K/km characteristic of a lithospheric-type environment, whereas the “sheared” ones show a lower gradient of about 1 K/km. This is a typical geotherm expected for a steady thermal state with an inflexion point at the depth of about 160 km corresponding to the lithosphere/asthenosphere boundary.     

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

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