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.     

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.     

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.     

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.     

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.     

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.     

The dilatant behaviour of sand–pile interface subjected to loading and stress relief

Property and behaviour of sand–pile interface are crucial to shaft resistance of piles. Dilation or contraction of the interface soil induces change in normal stress, which in turn influences the shear stress mobilised at the interface. Although previous studies have demonstrated this mechanism by laboratory tests and numerical simulations, the interface responses are not analysed systematically in terms of soil state (i.e. density and stress level). The objective of this study is to understand and quantify any increase in normal stress of different pile–soil interfaces when they are subjected to loading and stress relief. Distinct element modelling was carried out. Input parameters and modelling procedure were verified by experimental data from laboratory element tests. Parametric simulations of shearbox tests were conducted under the constant normal stiffness, constant normal load and constant volume boundary conditions. Key parameters including initial normal stress ( $ \sigma_{{{\text{n}}0}}^{\prime } $ ), initial void ratio (e ₀), normal stiffness constraining the interface and loading–unloading stress history were investigated. It is shown that mobilised stress ratio ( $ \tau /\sigma_{\text{n}}^{\prime } $ ) and normal stress increment ( $ \Updelta \sigma_{\text{n}}^{\prime } $ ) on a given interface are governed by $ \sigma_{{{\text{n}}0}}^{\prime } $ and e ₀. An increase in $ \sigma_{{{\text{n}}0}}^{\prime } $ from 100 to 400 kPa leads to a 30 % reduction in $ \Updelta \sigma_{\text{n}}^{\prime } $ . An increase in e ₀ from 0.18 to 0.30 reduces $ \Updelta \sigma_{\text{n}}^{\prime } $ by more than 90 %, and therefore, shaft resistance is much lower for piles in loose sands. A unique relationship between $ \Updelta \sigma_{\text{n}}^{\prime } $ and normal stiffness is established for different soil states. It can be applied to assess the shaft resistance of piles in soils with different densities and subjected to loading and stress relief. Fairly good agreement is obtained between the calculated shaft resistance based on the proposed relationship and the measured results in centrifuge model tests.     

Raman spectra of gypsum under pressure

Raman sprectra of a gypsum crystal were made at pressures between 0.001 and 7 kbar using He gas as the pressure medium. \(\frac{{{\text{d}}v}}{{dP}}\) values for bands in the range 3,600–100 cm^?1 were obtained. Comparison of results with \(\frac{{{\text{d}}v}}{{{\text{d}}T}}\) from the literature for temperatures of 77 and 300° K. shows that the internal modes of the SO₄ units are more sensitive to pressure than to temperature. The effect is small. Coupled H₂O-SO₄ translational modes are greatly affected by both pressure and temperature while coupled Ca-SO₄ mode are less so. It was found that stretching vibrations of water molecules were affected differently under pressure. The band at 3,500 cm^?1 is more greatly displaced by pressure \(\left( {\frac{{{\text{d}}v}}{{{\text{d}}P}} = {\text{2}}{\text{.11cm}}^{{\text{ - 1}}} /{\text{kbar}}} \right)\) than the band at 3,400 cm^?1 \(\left( {\frac{{{\text{d}}v}}{{{\text{d}}P}} \simeq {\text{2}}{\text{.11cm}}^{{\text{ - 1}}} /{\text{kbar}}} \right)\) . Assuming two different hydrogen bond intensities for the water molecules, one can attribute this difference in behavior of stretching modes to and increase in hydrogen bonding of one of the hydrogens which is exterior to the double H₂O planes in the gypsum structure. The great variety of pressure derivatives for the different types of vibrational modes observed indicates that each molecular unit readjusts internally to pressure induced volume changes and the some of the chemical bonds between the units are significantly affected.     

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.     

Interaction of B(OH)_3^0 and HCO_3^ - in Seawater: Formation of B(OH)_2 CO_3^ -

Boron is known to interact with a wide variety of protonated ligands(HL) creating complexes of the form B(OH)₂L^-.Investigation of the interaction of boric acid and bicarbonate in aqueoussolution can be interpreted in terms of the equilibrium $B(OH)_3^0 + HCO_3^ - \rightleftharpoons B(OH)_2 CO_3^ - + H_2 O$ The formation constant for this reaction at 25 °C and 0.7 molkg^-1 ionic strength is $K_{BC} = \left[ {B(OH)_2 CO_3^ - } \right]\left[ {B(OH)_3^0 } \right]^{ - 1} \left[ {HCO_3^ - } \right]^{ - 1} = 2.6 \pm 1.7$ where brackets represent the total concentration of each indicatedspecies. This formation constant indicates that theB(OH)₂ $CO_3^ - $ concentration inseawater at 25 °C is on the order of 2 μmol kg^-1. Dueto the presence of B(OH)₂ $CO_3^ - $ , theboric acid dissociation constant ( $K\prime _B $ ) in natural seawaterdiffers from $K\prime _B $ determined in the absence of bicarbonate byapproximately 0.5%. Similarly, the dissociation constants of carbonicacid and bicarbonate in natural seawater differ from dissociation constantsdetermined in the absence of boric acid by about 0.1%. Thesedifferences, although small, are systematic and exert observable influenceson equilibrium predictions relating CO₂ fugacity, pH, totalcarbon and alkalinity in seawater.     

The Hydrolysis of Al(III) in NaCl solutions: A Model for M(II), M(III), and M(IV) Ions

Understanding the identity and stability of the hydrolysis products of metals is required in order to predict their behavior in natural aquatic systems. Despite this need, the hydrolysis constants of many metals are only known over a limited range of temperature and ionic strengths. In this paper, we show that the hydrolysis constants of 31 metals [i.e. Mn(II), Cr(III), U(IV), Pu(IV)] are nearly linearly related to the values for Al(III) over a wide range of temperatures and ionic strengths. These linear correlations allow one to make reasonable estimates for the hydrolysis constants of +2, +3, and +4 metals from 0 to 300°C in dilute solutions and 0 to 100°C to 5 m in NaCl solutions. These correlations in pure water are related to the differences between the free energies of the free ion and complexes being almost equal $$ \Updelta {\text{G}}^\circ \left( {{\text{Al}}^{3 + } } \right) - \Updelta {\text{G}}^\circ \left( {{\text{Al}}\left( {\text{OH}} \right)_{j}^{{\left( {3 - j} \right)}} } \right) \cong \Updelta {\text{G}}^\circ \left( {{\text{M}}^{n + } } \right) - \Updelta {\text{G}}^\circ \left( {{\text{M}}\left( {\text{OH}} \right)_{j}^{{\left( {n - j} \right)}} } \right) $$ The correlation at higher temperatures is a result of a similar relationship between the enthalpies of the free ions and complexes $$ \Updelta {\text{H}}^\circ \left( {{\text{Al}}^{3 + } } \right) - \Updelta {\text{H}}^\circ \left( {{\text{Al}}\left( {\text{OH}} \right)_{j}^{3 - j} } \right) \cong \Updelta {\text{H}}^\circ \left( {{\text{M}}^{n + } } \right) - \Updelta {\text{H}}^\circ \left( {{\text{M}}\left( {\text{OH}} \right)_{j}^{n - j} } \right) $$ The correlations at higher ionic strengths are the result of the ratio of the activity coefficients for Al(III) being almost equal to that of the metal. $$ \gamma \left( {{\text{M}}^{n + } } \right)/\gamma \left( {{\text{M}}\left( {\text{OH}} \right)_{j}^{n - j} } \right) \cong \gamma \left( {{\text{Al}}^{3 + } } \right)/\gamma \left( {{\text{Al}}\left( {\text{OH}} \right)_{j}^{3 - j} } \right) $$ The results of this study should be useful in examining the speciation of metals as a function of pH in natural waters (e.g. hydrothermal fresh waters and NaCl brines).     

Ephesite,Na(LiAl2) [Al2Si2O10] (OH)2: I. thermal stability and standard state thermodynamic properties

Ephesite, Na(LiAl₂) [Al₂Si₂O₁₀] (OH)₂, has been synthesized for the first time by hydrothermal treatment of a gel of requisite composition at 300≦T(° C)≦700 and \(P_{H_2 O}\) upto 35 kbar. At \(P_{H_2 O}\) between 7 and 35 kbar and above 500° C, only the 2M₁ polytype is obtained. At lower temperatures and pressures, the 1M polytype crystallizes first, which then inverts to the 2M₁ polytype with increasing run duration. The X-ray diffraction patterns of the 1M and 2M₁ poly types can be indexed unambiguously on the basis of the space groups C2 and Cc, respectively. At its upper thermal stability limit, 2M₁ ephesite decomposes according to the reaction (1) $$\begin{gathered} {\text{Na(LiAl}}_{\text{2}} {\text{) [Al}}_{\text{2}} {\text{Si}}_{\text{2}} {\text{O}}_{{\text{10}}} {\text{] (OH)}}_{\text{2}} \hfill \\ {\text{ephesite}} \hfill \\ {\text{ = Na[AlSiO}}_{\text{4}} {\text{] + LiAl[SiO}}_{\text{4}} {\text{] + }}\alpha {\text{ - Al}}_{\text{2}} {\text{O}}_{\text{3}} {\text{ + H}}_{\text{2}} {\text{O}} \hfill \\ {\text{nepheline }}\alpha {\text{ - eucryptite corundum}} \hfill \\ \end{gathered}$$ Five reversal brackets for (1) have been established experimentally in the temperature range 590–750° C, at \(P_{H_2 O}\) between 400 and 2500 bars. The equilibrium constant, K, for this reaction may be expressed as (2) $$log K{\text{ = }}log f_{{\text{H}}_{\text{2}} O}^* = 7.5217 - 4388/T + 0.0234 (P - 1)T$$ where \(f_{H_2 O}^* = f_{H_2 O} (P,T)/f_{H_2 O}^0\) (1,T), with T given in degrees K, and P in bars. Combining these experimental data with known thermodynamic properties of the decomposition products in (1), the following standard state (1 bar, 298.15 K) thermodynamic data for ephesite were calculated: H _f,298.15 ⁰ =-6237372 J/mol, S _298.15 ⁰ =300.455 J/K·mol, G _298.15 ⁰ =-5851994 J/mol, and V _298.15 ⁰ =13.1468 J/bar·mol.     

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.     

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

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.     

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.     

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.     

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.