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.     

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.     

Partition of Ni between olivine and sulfide: the effect of temperature,f_{{\text{O}}_{\text{2}} } and f_{{\text{S}}_{\text{2}} }

The experimental distribution coefficient for Ni/ Fe exchange between olivine and monosulfide (K_D3) is 35.6±1.1 at 1385° C, \(f_{{\text{O}}_{\text{2}} } = 10^{ - 8.87} ,f_{{\text{S}}_{\text{2}} } = 10^{ - 1.02} \) , and olivine of composition Fo96 to Fo92. These are the physicochemical conditions appropriate to hypothesized sulfur-saturated komatiite magma. The present experiments equilibrated natural olivine grains with sulfide-oxide liquid in the presence of a (Mg, Fe)-alumino-silicate melt. By a variety of different experimental procedures, K _D3 is shown to be essentially constant at about 30 to 35 in the temperature range 900 to 1400° C, for olivine of composition Fo97 to FoO, monosulfide composition with up to 70 mol. % NiS, and a wide range of \(f_{{\text{O}}_{\text{2}} } \) and \(f_{{\text{S}}_{\text{2}} } \) .     

The ferric-ferrous ratio of natural silicate liquids equilibrated in air

Results of chemical analyses of glasses produced in 46 melting experiments in air at 1,350° C and 1,450° C on rocks ranging in composition from nephelinite to rhyolite have been combined with other published data to obtain an empirical equation relating in \((X_{{\text{Fe}}_{\text{2}} {\text{O}}_{\text{3}} }^{{\text{liq}}} /X_{{\text{FeO}}}^{{\text{liq}}} )\) to T, \(\ln f_{{\text{O}}_{\text{2}} } \) and bulk composition. The whole set of experimental data range over 1,200–1,450° C and oxygen fugacities of 10^?9.00 to 10^?0.69 bars, respectively. The standard errors of temperature and \(\log _{10} f_{{\text{O}}_{\text{2}} } \) predictions from this equation are 52° C and 0.5 units, respectively, for 186 experiments.     

Use of FePt alloys to eliminate the iron loss problem in 1 atmosphere gas mixing experiments: Theoretical and practical considerations

Theoretical and practical considerations are combined to place limits on the iron content of an FePt alloy that is in equilibrium with silicate melt, olivine and a gas phase of known \(f_{{\text{O}}_{\text{2}} }\) . Equilibrium constants are calculated for the reactions: (1) $$2{\text{Fe}}^{\text{o}} + {\text{SiO}}_{\text{2}} + {\text{O}}_{\text{2}} \rightleftharpoons {\text{Fe}}_{\text{2}} {\text{SiO}}_{\text{4}}$$ (2) $${\text{Fe}}^{\text{o}} + \frac{1}{2}{\text{O}}_{\text{2}} \rightleftharpoons {\text{FeO}}$$ . These equilibria may be used to choose an appropriate iron activity for the FePt alloy of an experiment. The temperature dependence of the equilibrium constants is calculated from experimental data. The Gibbs free energy of reaction (1) obtained using thermochemical data is in close agreement with ΔG_rxn calculated from the experimental data. Reaction (1) has the advantage that it is independent of the Fe²⁺/Fe³⁺ ratio of the melt, but is limited to applications where olivine is a crystallizing phase and requires a formulation for \(a_{{\text{SiO}}_{\text{2}} }^{{\text{liq}}}\) . Reaction (2) uses an empirical approximation for the FeO/Fe₂O₃ ratio of the liquid, and is independent of olivine saturation. However, it requires a formulation for a _FeO ^liq . Either equilibrium constant may be used to calculate the appropriate FePt alloy in equilibrium with a silicate melt. If experiments are conducted at an \(f_{{\text{O}}_{\text{2}} }\) parallel that of a buffer assemblage, a small range of FePt alloys may be used over a large temperature interval. For example, an alloy containing from 6 % to 9 % Fe by weight is in equilibrium with olivine-saturated tholeiites and komatiites at the quartzfayalite-magnetite buffer over the temperature interval 1,400° C to 1,100° C. Lunar basalt liquids in equilibrium with olivine at 1/2 log unit below the iron-wüstite buffer require an FePt alloy that contains 30–50 wt. % iron over a similar temperature interval.     

Experimentelle Untersuchung der Reaktion Dolomit + Kalifeldspat + H2O ⇌ Phlogopit + Calcit + CO2

The equilibrium curve for the reaction 3 dolomite + 1 K-feldspar + 1 H₂O=1 phlogopite + 3 calcite + 3 CO₂ was determined experimentally at a total gas pressure of 2000 bars using two different methods.

1. In the first case water alone was added to the reactants. The CO₂ component of the gas phase was producted solely by the reaction under favourable P-T conditions. This manner of carrying out the reaction is called the “water method”. With this method sufficient time must be allowed for the gas phase to attain a constant composition (see Fig. 1). Reverse reactions were carried out using reaction products of the forward reaction.
2. In the second case silver oxalate + water were added to the reactants. Breakdown of the silver oxalate leads to formation of a CO₂-H₂O gasphase of definite composition. At constant temperature and gas pressure the \(X_{{\text{CO}}_{\text{2}} } \) determines whether the reaction products will be phlogopite + calcite or dolomite + K-feldspar. In this case it is not necessary to wait for equilibrium to be attained. This method is abbreviated the “oxalate method”. Reactants for reverse reactions are not identical with the products of the forward reaction.

At high temperatures the results of the two different methods agree well (see Tables 1 and 2). Equilibrium was attained in one case at 490° C and \(X_{{\text{CO}}_{\text{2}} } \) of approximately 0.77, and in the other case at 520° C and \(X_{{\text{CO}}_{\text{2}} } \) of 0.90. At lower temperatures there are considerable differences in the results. With the water method an \(X_{{\text{CO}}_{\text{2}} } \) of about 0.25 was reached at 450° C. With the oxalate method dolomite K-feldspar and water still react with each other at even higher \(X_{{\text{CO}}_{\text{2}} } \) values. Phlogopite, calcite and CO₂ are formed together with metastable talc. There are no criteria to indicate which of the methods is the correct one at lower temperatures and in Fig. 2, therefore, both equilibrium curves are plotted.     

Gehlenite stability in the system CaO-Al2O3-SiO2-H2O-CO2

P, T, \(X_{{\text{CO}}_{\text{2}} }\) relations of gehlenite, anorthite, grossularite, wollastonite, corundum and calcite have been determined experimentally at P _f =1 and 4 kb. Using synthetic starting minerals the following reactions have been demonstrated reversibly

1. 2 anorthite+3 calcite=gehlenite+grossularite+3 CO₂.
2. anorthite+corundum+3 calcite=2 gehlenite+3 CO₂.
3. 3anorthite+3 calcite=2 grossularite+corundum+3CO₂.
4. grossularite+2 corundum+3 calcite=3 gehlenite+3 CO₂.
5. anorthite+2 calcite=gehlenite+wollastonite+2CO₂.
6. anorthite+wollastonite+calcite=grossularite+CO₂.
7. grossularite+calcite=gehlenite+2 wollastonite+CO₂.

In the T, \(X_{{\text{CO}}_{\text{2}} }\) diagram at P _f =1 kb two isobaric invariant points have been located at 770±10°C, \(X_{{\text{CO}}_{\text{2}} }\) =0.27 and at 840±10°C, \(X_{{\text{CO}}_{\text{2}} }\) =0.55. Formation of gehlenite from low temperature assemblages according to (4) and (2) takes place at 1 kb and 715–855° C, \(X_{{\text{CO}}_{\text{2}} }\) =0.1–1.0. In agreement with experimental results the formation of gehlenite in natural metamorphic rocks is restricted to shallow, high temperature contact aureoles.     

Experimental investigation of the mechanism of the reaction: 1 tremolite+11 dolomite ⇌ 8 forsterite+13 calcite+9 CO2+1H2O

Stoichiometric mixtures of tremolite and dolomite were heated to 50° C above equilibrium temperatures to form forsterite and calcite. The pressure of the CO₂-H₂O fluid was 5 Kb and \(X_{{\text{CO}}_{\text{2}} }\) varied from 0.1 to 0.6. The extent of the conversion was determined by the amount of CO₂ produced. The resulting mixtures of unreacted tremolite and dolomite and of newly-formed forsterite and calcite were examined with a scanning electron microscope. All tremolite and dolomite grains showed obvious signs of dissolution. At fluid compositions with \(X_{{\text{CO}}_{\text{2}} }\) less than about 0.4, the forsterite and calcite crystals are randomly distributed throughout the charges, indicating that surfaces of the reactants are not a controlling factor with respect to the sites of nucleation of the products. A change is observed when \(X_{{\text{CO}}_{\text{2}} }\) is greater than about 0.4; the forsterite and calcite crystals now nucleate and grow at the surface of the dolomite grains, thus indicating a change in mechanism at medium CO₂ concentrations. As the reaction progresses, the dolomite grains become more and more surrounded by forsterite and calcite, finally forming armoured relics of dolomite. Under experimental conditions this characteristic texture can only be formed if the CO₂-concentration is greater than about 40 mole %. These findings make it possible to estimate the CO₂-concentration from the texture of the dolomite+tremolite+forsterite+calcite assemblage. The results suggest a dissolution-precipitation mechanism for the reaction investigated. In a simplified form it consists of the following 4 steps:

1. Dissolution of the reactants tremolite and dolomite.
2. Diffusion of the dissolved constituents in the fluid.
3. Heterogeneous nucleation of the product minerals.
4. Growth of forsterite and calcite from the fluid.

Two possible explanations are discussed for the development of the amoured texture at \(X_{{\text{CO}}_{\text{2}} }\) above 0.4. The first is based upon the assumption that dolomite has a lower rate of dissolution than tremolite at high \(X_{{\text{CO}}_{\text{2}} }\) values resulting in preferential calcite and forsterite nucleation and growth on the dolomite surface. An alternative explanation is the formation of a raised CO₂ concentration around the dolomite grains at high \(X_{{\text{CO}}_{\text{2}} }\) values, leading to product precipitation on the dolomite crystals.     

An empirical Redlich-Kwong-type equation of state for water to 1,000 °C and 200 Kbar

An empirically derived Redlich-Kwong type of equation of state (ERK) is proposed for H₂O, expressing a, the term related to the attraction between the molecules, as a pressure-independent function of temperature, and b, the covolume, as a temperature-independent function of pressure. The coefficients of a(T) and b(P) were derived by least squares non-linear regression, using P-V-T data given by Burnham et al. (1969b) and Rice and Walsh (1957) in conjunction with more precise recent data obtained by Tanishita et al. (1976), Hilbert (1979) and Schmidt (1979): $$a(T) = 1.616 x 10^8 - 4.989 x 10^4 T - 7.358 x 10^9 T^{ - 1} $$ and $$ = \frac{{1 + 3.4505x 10^{--- 4} P + 3.8980x 10^{--- 9} P^2 - 2.7756x 10^{--- 15} P^3 }}{{6.3944x 10^{--- 2} + 2.3776x 10^{--- 5} + 4.5717x 10^{--- 10} P^2 }}$$ , where T is expressed in Kelvin and P in bars. The ERK works very well at upper mantle conditions, at least up to 200 kbar and 1,000 °C. At subcritical conditions and those somewhat above the critical point, it still reproduces the molar Gibbs energy, \(\tilde G_{{\text{H}}_{\text{2}} {\text{O}}} \) , with a maximum deviation of 400 joules. Thus, for the purpose of calculation of geologically interesting heterogeneous equilibria, it predicts the thermodynamic properties of H₂O well enough. The values of molar volume, \(\tilde V_{{\text{H}}_{\text{2}} {\text{O}}} \) , and \(\tilde G_{{\text{H}}_{\text{2}} {\text{O}}} \) are tabulated in the appendix over a considerable P-T range. A FORTRAN program generating these functions as well as a FORTRAN subroutine for calculating the fugacity values, \(f_{{\text{H}}_{\text{2}} {\text{O}}} \) for incorporation into existing programs, are available upon request.     

Quaternary acid magma in New Zealand

The voluminous Pleistocene—Recent Taupo rhyolites typically contain phenocrysts of plagioclase (oligoclase-labradorite), quartz, titanomagnetite, ilmenite, and ferromagnesian silicates. Ferromagnesian assemblages correlate with well defined Fe-Ti oxide equilibration temperature ranges and allow the rhyolites to be subdivided as follows: (1) Cummingtonite (c)—calcic hornblende (hb)—orthopyroxene (opx); 725–755°C, (2) Hb-opx, 750–825°C, (3) Biotite-hb-(c-opx), 720–765°C, (4) Opx-clinopyroxene (cpx), 860–915°C, (5) Fe olivine-opx-cpx, one sample with temperature of 900°C. Plagioclase and orthopyroxene phenocryst compositions typically exhibit a range of composition up to ～20 mol.%. Calculated average phenocryst equilibration pressures (P _total) range between 0.5–4.9 kb, and average 2.2 kb (～7–8 km depth), indicating upper crustal crystallization. These calculations are very sensitive to variations in phenocryst composition. Calculated \(/_{{\text{H}}_2 {\text{O}}} \) for the amphibole and biotite-bearing rhyolites indicate phenocryst equilibration under \(P_{{\text{H}}_2 {\text{O}}} \simeq P_{{\text{total}}} \) , with \(X_{{\text{H}}_2 {\text{O}}} \) ～0.17–0.24 (5–8 wt. %). The precipitation of cummingtonite is thus temperature dependent, the upper limit being close to 760°C. Eruptive mechanisms of the lavas, pumices, and ash-flow deposits are evidently not primarily controlled by temperature, P _total, \(P_{{\text{H}}_2 {\text{O}}} \) , or crystal content of the magmas, and explanations must lie in kinetic and fluid dynamic behavior of the magmas. For the Taupo rhyolites, it is suggested that the critical size of a magma body (i.e. Rayleigh number) is a controlling factor in that it will influence the convective regime; fully turbulent convection is deduced to have occurred within the larger magma bodies. One consequence is intense vesiculation, prior to eruption, within the uppermost zones of these magma chambers, and the voluminous pumice deposits are believed to emanate from such chambers. Oscillatory compositional zoning within pyroxene phenocrysts is consistent with magma convection.     

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.     

Activity-composition relations of the magnesiowüstite solid solution series in equilibrium with metallic iron in the temperature range 1050–1400 K

The equilibrium $${\text{(1}} - y{\text{)Fe}}_{(s)} + \tfrac{{\text{1}}}{{\text{2}}}{\text{O}}_{{\text{2(g)}}} \rightleftarrows {\text{Fe}}_{{\text{1}} - y} {\text{O}}_{{\text{(}}s,{\text{ in MW)}}} $$ was studied by measuring oxygen potentials for a range of different magnesiowüstite compositions relative to those of the iron-wüstite system in an oxygen concentration cell involving yttria stabilized zirconia as the solid electrolyte. The temperature range covered was 1050 to 1400 K. Separate measurements of the mole fraction of trivalent iron in magnesiowüstite (x(Fe³⁺)) were made and the composition dependence ofx(Fe³⁺) was taken into account in calculations of the activity-composition relations of FeO, Fe_2/3O and MgO.     

Garnet-clinopyroxene geothermometer

A garnet-clinopyroxene geothermometer based on the available experimental data on compositions of coexisting phases in the system MgO-FeO-MnO-Al₂O₃-Na₂O-SiO₂ is as follows: $$T({\text{}}K) = \frac{{8288 + 0.0276 P {\text{(bar)}} + Q1 - Q2}}{{1.987 \ln K_{\text{D}} + 2.4083}}$$ where P is pressure, and Q1, Q2, and K _D are given by the following equations $$Q1 = 2,710{\text{(}}X_{{\text{Fe}}} - X_{{\text{Mg}}} {\text{)}} + 3,150{\text{ }}X_{{\text{Ca}}} + 2,600{\text{ }}X_{{\text{Mn}}} $$ (mole fractions in garnet) $$\begin{gathered}Q2 = - 6,594[X_{{\text{Fe}}} {\text{(}}X_{{\text{Fe}}} - 2X_{{\text{Mg}}} {\text{)]}} \hfill \\{\text{ }} - 12762{\text{ [}}X_{{\text{Fe}}} - X_{{\text{Mg}}} (1 - X_{{\text{Fe}}} {\text{)]}} \hfill \\{\text{ }} - 11,281[X_{{\text{Ca}}} (1 - X_{{\text{Al}}} ) - 2X_{{\text{Mg}}} 2X_{{\text{Ca}}} ] \hfill \\{\text{ + 6137[}}X_{{\text{Ca}}} (2X_{{\text{Mg}}} + X_{{\text{Al}}} )] \hfill \\{\text{ + 35,791[}}X_{{\text{Al}}} (1 - 2X_{{\text{Mg}}} )] \hfill \\{\text{ + 25,409[(}}X_{{\text{Ca}}} )^2 ] - 55,137[X_{{\text{Ca}}} (X_{{\text{Mg}}} - X_{{\text{Fe}}} )] \hfill \\{\text{ }} - 11,338[X_{{\text{Al}}} (X_{{\text{Fe}}} - X_{{\text{Mg}}} )] \hfill \\\end{gathered} $$ [mole fractions in clinopyroxene Mg = MgSiO₃, Fe = FeSiO₃, Ca = CaSiO₃, Al = (Al₂O₃-Na₂O)] K _D = (Fe/Mg) in garnet/(Fe/Mg) in clinopyroxene. Mn and Cr in clinopyroxene, when present in small concentrations are added to Fe and Al respectively. Fe is total Fe²⁺+Fe³⁺.     

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.     

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}}} \)      

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

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.     

An experimental study of Ti and Zr partitioning among zircon, rutile, and granitic melt

In order to evaluate the effect of trace and minor elements (e.g., P, Y, and the REEs) on the high-temperature solubility of Ti in zircon (zrc), we conducted 31 experiments on a series of synthetic and natural granitic compositions [enriched in TiO₂ and ZrO₂; Al/(Na + K) molar ~1.2] at a pressure of 10 kbar and temperatures of ~1,400 to 1,200 °C. Thirty of the experiments produced zircon-saturated glasses, of which 22 are also saturated in rutile (rt). In seven experiments, quenched glasses coexist with quartz (qtz). SiO₂ contents of the quenched liquids range from 68.5 to 82.3 wt% (volatile free), and water concentrations are 0.4–7.0 wt%. TiO₂ contents of the rutile-saturated quenched melts are positively correlated with run temperature. Glass ZrO₂ concentrations (0.2–1.2 wt%; volatile free) also show a broad positive correlation with run temperature and, at a given T, are strongly correlated with the parameter (Na + K + 2Ca)/(Si·Al) (all in cation fractions). Mole fraction of ZrO₂ in rutile $ \left( {\mathop X\nolimits_{{{\text{ZrO}}_{ 2} }}^{\text{rt}} } \right) $ in the quartz-saturated runs coupled with other 10-kbar qtz-saturated experimental data from the literature (total temperature range of ~1,400 to 675 °C) yields the following temperature-dependent expression: $ {\text{ln}}\left( {\mathop X\nolimits_{{{\text{ZrO}}_{ 2} }}^{\text{rt}} } \right) + {\text{ln}}\left( {a_{{{\text{SiO}}_{2} }} } \right) = 2.638(149) - 9969(190)/T({\text{K}}) $ , where silica activity $ a_{{{\text{SiO}}_{2} }} $ in either the coexisting silica polymorph or a silica-undersaturated melt is referenced to α-quartz at the P and T of each experiment and the best-fit coefficients and their uncertainties (values in parentheses) reflect uncertainties in T and $ \mathop X\nolimits_{{{\text{ZrO}}_{2} }}^{\text{rt}} $ . NanoSIMS measurements of Ti in zircon overgrowths in the experiments yield values of ~100 to 800 ppm; Ti concentrations in zircon are positively correlated with temperature. Coupled with values for $ a_{{{\text{SiO}}_{2} }} $ and $ a_{{{\text{TiO}}_{2} }} $ for each experiment, zircon Ti concentrations (ppm) can be related to temperature over the range of ~1,400 to 1,200 °C by the expression: $ \ln \left( {\text{Ti ppm}} \right)^{\text{zrc}} + \ln \left( {a_{{{\text{SiO}}_{2} }} } \right) - \ln \left( {a_{{{\text{TiO}}_{2} }} } \right) = 13.84\left( {71} \right) - 12590\left( {1124} \right)/T\left( {\text{K}} \right) $ . After accounting for differences in $ a_{{{\text{SiO}}_{2} }} $ and $ a_{{{\text{TiO}}_{2} }} $ , Ti contents of zircon from experiments run with bulk compositions based on the natural granite overlap with the concentrations measured on zircon from experiments using the synthetic bulk compositions. Coupled with data from the literature, this suggests that at T ≥ 1,100 °C, natural levels of minor and trace elements in “granitic” melts do not appear to influence the solubility of Ti in zircon. Whether this is true at magmatic temperatures of crustal hydrous silica-rich liquids (e.g., 800–700 °C) remains to be demonstrated. Finally, measured $ D_{\text{Ti}}^{{{\text{zrc}}/{\text{melt}}}} $ values (calculated on a weight basis) from the experiments presented here are 0.007–0.01, relatively independent of temperature, and broadly consistent with values determined from natural zircon and silica-rich glass pairs.     

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.     

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.