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.     

Physicochemical parameter charts for gases in fluid inclusions

A great wealth of analytical data for fluid inclusions in minerals indicate that the major species of gases in fluid inclusions are H₂O, CO₂, CO, CH₄, H₂ and O₂. Three basic chemical reactions are supposed to prevail in rock-forming and ore-forming fluids: $$\begin{gathered} H_2 + 1/2{\text{ O}}_{\text{2}} = H_2 O, \hfill \\ CO + 1/2{\text{ O}}_{\text{2}} = CO_2 , \hfill \\ CH_4 + 2{\text{O}}_{\text{2}} = CO_2 + 2H_2 O, \hfill \\ \end{gathered} $$ and equilibria are reached among them. \(\lg f_{O_2 } - T,{\text{ }}\lg f_{CO_2 } - T\) and Eh-T charts for petrogenesis and minerogenesis in the supercritical state have been plotted under different pressures. On the basis of these charts \(f_{O^2 } ,{\text{ }}f_{CO_2 } \) , Eh, equilibrium temperature and equilibrium pressure can be readily calculated. In this paper some examples are presented to show their successful application in the study of the ore-forming environments of ore deposits.     

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.     

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.     

Oxygen isotope fractionation between rutile and water

Synthetic rutile-water fractionations (1000 ln α) at 775, 675, and 575° C were found to be ?2.8, ?3.5, and ?4.8, respectively. Partial exchange experiments with natural rutile at 575° C and with synthetic rutile at 475° C failed to yield reliable fractionations. Isotopic fractionation within the range 575–775° C may be expressed as follows: 1 $$1000\ln \alpha ({\rm T}i{\rm O}_{2 } - H_2 O) = - 4.1 \frac{{10^6 }}{{T_{k^2 } }} + 0.96$$ . Combined with previously determined quartz-water fractionations, the above data permit calibration of the quartz-rutile geothermometer: 1 $$1000\ln \alpha ({\text{S}}i{\rm O}_{2 } - Ti{\rm O}_{2 } ) = 6.6 \frac{{10^6 }}{{T_{k^2 } }} - 2.9$$ . When applied to B-type eclogites from Europe, as an example, the latter equation yields a mean equilibration temperature of 565° 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.     

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

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

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.     

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.     

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.     

Petrological cannibalism: the chemical and textural consequences of incremental magma body growth

The textures of minerals in volcanic and plutonic rocks testify to a complexity of processes in their formation that is at odds with simple geochemical models of igneous differentiation. Zoning in plagioclase feldspar is a case in point. Very slow diffusion of the major components in plagioclase means that textural evidence for complex magmatic evolution is preserved, almost without modification. Consequently, plagioclase affords considerable insight into the processes by which magmas accumulate in the crust prior to their eventual eruption or solidification. Here, we use the example of the 1980–1986 eruptions of Mount St. Helens to explore the causes of textural complexity in plagioclase and associated trapped melt inclusions. Textures of individual crystals are consistent with multiple heating and cooling events; changes in total pressure (P) or volatile pressure ( $P_{{{\text{H}}_{ 2} {\text{O}}}}$ P H 2 O ) are less easy to assess from textures alone. We show that by allying textural and chemical analyses of plagioclase and melt inclusions, including volatiles (H₂O, CO₂) and slow-diffusing trace elements (Sr, Ba), to published experimental studies of Mount St. Helens magmas, it is possible to disambiguate the roles of pressure and temperature to reconstruct magmatic evolutionary pathways through temperature–pressure–melt fraction (T– $P_{{{\text{H}}_{ 2} {\text{O}}}}$ P H 2 O –F) space. Our modeled crystals indicate that (1) crystallization starts at $P_{{{\text{H}}_{ 2} {\text{O}}}}$ P H 2 O  > 300 MPa, consistent with prior estimates from melt inclusion volatile contents, (2) crystal cores grow at $P_{{{\text{H}}_{ 2} {\text{O}}}}$ P H 2 O  = 200–280 MPa at F = 0.65–0.7, (3) crystals are transferred to $P_{{{\text{H}}_{ 2} {\text{O}}}}$ P H 2 O  = 100–130 MPa (often accompanied by 10–20 °C of heating), where they grow albitic rims of varying thicknesses, and (4) the last stage of crystallization occurs after minor heating at $P_{{{\text{H}}_{ 2} {\text{O}}}}$ P H 2 O  ~ 100 MPa to produce characteristic rim compositions of An₅₀. We hypothesize that modeled $P_{{{\text{H}}_{ 2} {\text{O}}}}$ P H 2 O decreases in excess of ~50 MPa most likely represent upward transport through the magmatic system. Small variations in modeled $P_{{{\text{H}}_{ 2} {\text{O}}}}$ P H 2 O , in contrast, can be effected by fluxing the reservoir with CO₂-rich vapors that are either released from deeper in the system or transported with the recharge magma. Temperature fluctuations of 20–40 °C, on the other hand, are an inevitable consequence of incremental, or pulsed, assembly of crustal magma bodies wherein each pulse interacts with ancestral, stored magmas. We venture that this “petrological cannibalism” accounts for much of the plagioclase zoning and textural complexity seen not only at Mount St. Helens but also at arc magmas generally. More broadly we suggest that the magma reservoir below Mount St. Helens is dominated by crystal mush and fed by frequent inputs of hotter, but compositionally similar, magma, coupled with episodes of magma ascent from one storage region to another. This view both accords with other independent constraints on the subvolcanic system at Mount St. Helens and supports an emerging view of many active magmatic systems as dominantly super-solidus, rather than subliquidus, bodies.     

Phase relationships in granulitic metapelites from the Ivrea-Verbano zone (Northern Italy)

Paragneisses of the Ivrea-Verbano zone exhibit over a horizontal distance of 5 km mineralogical changes indicative of the transition from amphibolite to granulite facies metamorphism. The most obvious change is the progressive replacement of biotite by garnet via the reaction: a $${\text{Biotite + sillimanite + quartz }} \to {\text{ Garnet + K - feldspar + H}}_{\text{2}} {\text{O}}$$ which results in a systematic increase in the modal ratio g = (garnet)/(garnet + biotite) with increasing grade. The systematic variations in garnet and biotite contents of metapelites are also reflected by the compositions of these phases, both of which become more magnesian with increasing metamorphic grade. The pressure of metamorphism has been estimated from the Ca₃Al₂Si₃O₁₂ contents of garnets coexisting with plagioclase, sillimanite and quartz. These phases are related by the equilibrium: b $$\begin{gathered} 3 CaAl_2 {\text{Si}}_{\text{2}} {\text{O}}_{\text{8}} \rightleftharpoons Ca_3 Al_2 {\text{Si}}_{\text{3}} {\text{O}}_{{\text{12}}} + 2 Al_2 {\text{SiO}}_{\text{5}} + {\text{SiO}}_{\text{2}} \hfill \\ plagioclase garnet sillimanite quartz \hfill \\ \end{gathered} $$ which has been applied to these rocks using the available data on the mixing properties of plagioclase and garnet solid solutions. Temperature and f _{H 2O} estimates have been made in a similar way using thermodynamic data on the biotite-garnet reaction (a) and the approximate solidus temperatures of paragneisses. Amphibolite to granulite facies metamorphism in the Ivrea-Verbano zone took place in the P-T ranges 9–11 kb and 700–820 °C. The differences in temperature and pressure of metamorphism between g= 0 and g = 1 (5 kms horizontal distance) were less than 50° C and approximately 1 kb. Retrogression and re-equilibration of garnets and biotites in the metapelites extended to temperatures more than 50° C below and pressures more than 1.5 kb below the peak of metamorphism, the degree of retrogression increasing with decreasing grade of the metamorphic “peak”. The pressure and temperature of the peak of metamorphism are not inconsistent with the hypothesis that the Ivrea-Verbano zone is a slice of upthrusted lower crust from the crust-mantle transition region, although it appears that the thermal gradient was too low for the zone to represent a near-vertical section through the crust. The most reasonable explanation of the granulite facies metamorphism is that it arose through intrusion of mafic rocks into a region already undergoing recrystallisation under amphibolite facies conditions.     

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.     

A multi-isotope approach for understanding sources of water,carbon and sulfur in natural springs of the Central Appalachian region

Natural springs have been reliable sources of domestic water and have allowed for the development of recreational facilities and resorts in the Central Appalachians. The structural history of this area is complex and it is unknown whether these natural springs receive significant recharge from modern precipitation or whether they discharge old water recharged over geological times scales. The main objective of this study was to use stable isotopes of water ( $\delta^{18} {\text{O}}_{{{\text{H}}_{2} {\text{O}}}}$ and $\delta^{2} {\text{H}}_{{{\text{H}}_{2} {\text{O}}}}$ ), dissolved inorganic carbon ( $\delta^{13} {\text{C}}_{\text{DIC}}$ ) and dissolved sulfate ( $\delta^{34} {\text{S}}_{{{\text{SO}}_{4} }}$ and $\delta^{18} {\text{O}}_{{{\text{SO}}_{4} }}$ ) to delineate sources of water, carbon and sulfur in several natural springs of the region. Our preliminary isotope data indicate that all springs are being recharged by modern precipitation. The oxygen isotope composition indicates that waters in thermal springs did not encounter the high temperatures required for O isotope exchange between the water and silicate/carbonate minerals, and/or the residence time of water in the aquifers was short due to high flow rates. The carbon isotopic composition of dissolved inorganic carbon and sulfur/oxygen isotopic composition of dissolved sulfate provide evidence of low-temperature water–rock interactions and various biogeochemical transformations these waters have undergone along their flow path.     

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.     

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.     

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.     

The Activity of H2O,Fugacity of O2 and Variation of Their Chemical Potential During Metamorphic Peak of the Granulite Complex Near Taipingzhai—Louzishan Region,East Hebei Province

The thermodynamic calculation of dehydration reacton suggests very low activity of H2O during metamorphic peak of the Archaean granulite complex in the region studied.The αH2O values for Al-rich gneiss and hypersthene biotite gneiss-granulite in the Taipingzhai region are usually between 0.10 and 0.20,and those in the Louzishan region are 0.15-0.25.The fugacity of O2 in terms of lgf O2 in whole region ranges form-8to-14.The average coefficients of (δμH2O/δHMg^Bt)and(δμO2/δXMg^Bt)in the Taipingzhai region are-0.293 and-1.60 respectively,and those in the Louzishan region are-0.364and-1.420.The activity of H2O is very low in the whole region,but its values and other data mentioned above are considerably constant from place to place within a given region,even in rocks of dirrerent lithological characters.However,they show a certain gradient between different regions.Such characteristics are compatible with the genetic mechanism known as“carbonic metamorphism” put forward by Newton et al.,i.e.,the α H2O during the peak stage is controlled by permeation of pervasive CO2 influx of the mantle source,and shows features of external buffering.     

Redox-controlled iron isotope fractionation during magmatic differentiation: an example from the Red Hill intrusion, S. Tasmania

This study presents accurate and precise iron isotopic data for 16 co-magmatic rocks and 6 pyroxene–magnetite pairs from the classic, tholeiitic Red Hill sill in southern Tasmania. The intrusion exhibits a vertical continuum of compositions created by in situ fractional crystallisation of a single injection of magma in a closed igneous system and, as such, constitutes a natural laboratory amenable to determining the causes of Fe isotope fractionation in magmatic rocks. Early fractionation of pyroxenes and plagioclase, under conditions closed to oxygen exchange, gives rise to an iron enrichment trend and an increase in $ f_{{{\text{O}}_{2} }} $ of the melt relative to the Fayalite–Magnetite–Quartz (FMQ) buffer. Enrichment in Fe³⁺/ΣFe_melt is mirrored by δ⁵⁷Fe, where ^VIFe²⁺-bearing pyroxenes partition ⁵⁷Fe-depleted iron, defining an equilibrium pyroxene-melt fractionation factor of $ \Updelta^{57} {\text{Fe}}_{{{\text{px}} - {\text{melt}}}} \le - 0.25\,\permille \times 10^{6} /T^{2} $ . Upon magnetite saturation, the $ f_{{{\text{O}}_{2} }} $ and δ⁵⁷Fe of the melt fall, commensurate with the sequestration of the oxidised, ⁵⁷Fe-enriched iron into magnetite, quantified as $ \Updelta^{57} {\text{Fe}}_{{{\text{mtn}} - {\text{melt}}}} = + 0.20\,\permille \times 10^{6} /T^{2} $ . Pyroxene–magnetite pairs reveal an equilibrium fractionation factor of $ \Updelta^{57} {\text{Fe}}_{{{\text{mtn}} - {\text{px}}}} \approx + 0.30\,\permille $ at 900–1,000?°C. Iron isotopes in differentiated magmas suggest that they may act as an indicator of their oxidation state and tectonic setting.