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

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

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

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³⁺.     

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

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

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.     

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

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

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.     

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.     

Experimental study of subsolidus phase relations and mixing properties of pyroxene and plagioclase in the system Na2O-CaO-Al2O3-SiO2

In the system Na₂O-CaO-Al₂O₃-SiO₂ (NCAS), the equilibrium compositions of pyroxene coexisting with grossular and corundum were experimentally determined at 40 different P-T conditions (1,100–1,400° C and 20.5–38 kbar). Mixing properties of the Ca-Tschermak — Jadeite pyroxene inferred from the data are (J, K): $$\begin{gathered} G_{Px}^{xs} = X_{{\text{CaTs}}} X_{{\text{Jd}}} [14,810 - 7.15T - 5,070(X_{{\text{CaTs}}} - X_{{\text{Jd}}} ) \hfill \\ {\text{ }} - 3,350(X_{{\text{CaTs}}} - X_{{\text{Jd}}} )^2 ] \hfill \\ \end{gathered} $$ The excess entropy is consistent with a complete disorder of cations in the M2 and the T site. Compositions of coexisting pyroxene and plagioclase were obtained in 11 experiments at 1,190–1,300° C/25 kbar. The data were used to infer an entropy difference between low and high anorthite at 1,200° C, corresponding to the enthalpy difference of 9.6 kJ/mol associated with the C \(\bar 1\) =I \(\bar 1\) transition in anorthite as given by Carpenter and McConnell (1984). The resulting entropy difference of 5.0 J/ mol · K places the transition at 1,647° C. Plagioclase is modeled as ideal solutions, C \(\bar 1\) and I \(\bar 1\) , with a non-first order transition between them approximated by an empirical expression (J, bar, K): $$\Delta G_T = \Delta G_{1,473} \left[ {1 - 3X_{Ab} \tfrac{{T^4 - 1,473^4 }}{{\left( {1,920 - 0.004P} \right)^4 - 1,473^4 }}} \right],$$ where $$\Delta G_{1,473} = 9,600 - 5.0T - 0.02P$$ The derived mixing properties of the pyroxene and plagioclase solutions, combined with the thermodynamic properties of other phases, were used to calculate phase relations in the NCAS system. Equilibria involving pyroxene+plagioclase +grossular+corundum and pyroxene+plagioclase +grossular+kyani te are suitable for thermobarometry. Albite is the most stable plagioclase.     

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.     

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.     

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.     

The combined effect of F and H2O on interdiffusion between peralkaline dacitic and rhyolitic melts

The effective binary diffusion coefficient (EBDC) of silicon has been measured during the interdiffusion of peralkaline, fluorine-bearing (1.3 wt% F), hydrous (3.3 and 6 wt% H₂O), dacitic and rhyolitic melts at 1.0 GPa and temperatures between 1100°C and 1400°C. From Boltzmann-Matano analysis of diffusion profiles the diffusivity of silicon at 68 wt% SiO₂ can be described by the following Arrhenius equations (with standard errors): $$\begin{gathered} {\text{with 1}}{\text{.3 wt\% F and 3}}{\text{.3\% H}}_{\text{2}} {\text{O:}} \hfill \\ {\text{D}}_{{\text{Si}}} = \begin{array}{*{20}c} { + {\text{3}}{\text{.59}}} \\ {{\text{3}}{\text{.66}} \times {\text{10}}^{ - {\text{9}}} } \\ { - {\text{1}}{\text{.86}}} \\ \end{array} {\text{exp}}\left( {{{ - {\text{86}}{\text{.1}} \pm {\text{8}}{\text{.9}}} \mathord{\left/ {\vphantom {{ - {\text{86}}{\text{.1}} \pm {\text{8}}{\text{.9}}} {{\text{RT}}}}} \right. \kern-\nulldelimiterspace} {{\text{RT}}}}} \right) \hfill \\ {\text{with 1}}{\text{.3 wt\% F and 6}}{\text{.0\% H}}_{\text{2}} {\text{O:}} \hfill \\ {\text{D}}_{{\text{Si}}} = \begin{array}{*{20}c} { + {\text{3}}{\text{.59}}} \\ {{\text{3}}{\text{.51}} \times {\text{10}}^{ - {\text{8}}} } \\ { - {\text{1}}{\text{.77}}} \\ \end{array} {\text{exp}}\left( {{{ - {\text{109}}{\text{.5}} \pm {\text{8}}{\text{.9}}} \mathord{\left/ {\vphantom {{ - {\text{109}}{\text{.5}} \pm {\text{8}}{\text{.9}}} {{\text{RT}}}}} \right. \kern-\nulldelimiterspace} {{\text{RT}}}}} \right) \hfill \\ \end{gathered} $$ where D is in m²s^?1 and activation energies are in kJ/mol. Diffusivities measured at 64 and 72 wt% SiO₂ are only slightly different from those at 68 wt% SiO₂ and frequently all measurements are within error of each other. Silicon, aluminum, iron, magnesium, and calcium EBDCs were also calculated from diffusion profiles by error function inversion techniques assuming constant diffusivity. With one exception, silicon EBDCs calculated by error function techniques are within error of Boltzmann-Matano EBDCs. Average diffusivities of Fe, Mg, and Ca were within a factor of 2.5 of silicon diffusivities whereas Al diffusivities were approximately half those of silicon. Alkalies diffused much more rapidly than silicon and non-alkalies, however their diffusivities were not quantitatively determined. Low activation energies for silicon EBDCs result in rapid diffusion at magmatic temperatures. Assuming that water and fluorine exert similar effects on melt viscosity at high temperatures, the viscosity can be calculated and used in the Eyring equation used to determine diffusivities, typically to within a factor of three of those measured in this study. This correlation between viscosity and diffusivity can be inverted to calculate viscosities of fluorine- and water-bearing granitic melts at magmatic temperatures; these viscosities are orders of magnitude below those of hydrous granitic melts and result in more rapid and effective separation of granitic magmas from partially molten source rocks. Comparison of Arrhenius parameters for diffusion measured in this study with Arrhenius parameters determined for diffusion in similar compositions at the same pressure demonstrates simple relationships between Arrhenius parameters, activation energy-E_a, kJ/mol, pre-exponential factor-D_o, m²s^?1, and the volatile, X=F or OH^?, to oxygen, O, ratio of the melt {(X/X+O)}: $$\begin{gathered} {\text{E}}a = - {\text{1533\{ }}{{\text{X}} \mathord{\left/ {\vphantom {{\text{X}} {\left( {{\text{X}} + {\text{O}}} \right)}}} \right. \kern-\nulldelimiterspace} {\left( {{\text{X}} + {\text{O}}} \right)}}{\text{\} }} + {\text{213}}{\text{.3}} \hfill \\ {\text{D}}_{\text{O}} = {\text{2}}{\text{.13}} \times {\text{10}}^{ - {\text{6}}} {\text{exp}}\left[ { - {\text{6}}{\text{.5\{ }}{{\text{X}} \mathord{\left/ {\vphantom {{\text{X}} {\left( {{\text{X}} + {\text{O}}} \right)}}} \right. \kern-\nulldelimiterspace} {\left( {{\text{X}} + {\text{O}}} \right)}}{\text{\} }}} \right] \hfill \\ \end{gathered} $$ These relationships can be used to estimate diffusion in various melts of dacitic to rhyolitic composition containing both fluorine and water. Calculations for the contamination of rhyolitic melts by dacitic enclaves at 800°C and 700°C provide evidence for the virtual inevitability of diffusive contamination in hydrous and fluorine-bearing magmas if they undergo magma mixing of any form.     

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.     

Kinetic rate laws as derived from order parameter theory II: Interpretation of experimental data by laplace-transformation,the relaxation spectrum,and kinetic gradient coupling between two order parameters

A unifying theory of kinetic rate laws, based on order parameter theory, is presented. The time evolution of the average order parameter is described by $$\langle Q\rangle \propto \smallint P(x)e^{^{^{^{^{^{^{ - xt} } } } } } } dx = L(P)$$ where t is the time, x is the effective inverse susceptibility, and L indicates the Laplace transformation. The probability function P(x) can be determined from experimental data by inverse Laplace transformation. Five models are presented:

1. Polynomial distributions of P(x) lead to Taylor expansions of 〈Q〉 as $$\langle Q\rangle = \frac{{\rho _1 }}{t} + \frac{{\rho _2 }}{{t^2 }} + ...$$
2. Gaussian distributions (e.g. due to defects) lead to a rate law $$\langle Q\rangle = e^{ - x_0 t} e^{^{^{^{^{\frac{1}{2}\Gamma t^2 } } } } } erfc\left( {\sqrt {\frac{\Gamma }{2}} t} \right)$$ where x ₀ is the most probable inverse time constant, Γ is the Gaussian line width and erfc is the complement error integral.
3. Maxwell distributions of P are equivalent to the rate law 〈Q〉∝e^?k√t .
4. Pseudo spin glasses possess a logarithmic rate law 〈Q〉∝lnt.
5. Power laws with P(x)=x ^a lead to a rate law: ln〈Q〉=-(α + 1) ln t.

The power spectra of Q are shown for Gaussian distributions and pseudo spin glasses. The mechanism of kinetic gradient coupling between two order parameters is evaluated.     

Thermodynamic properties of glaucophane new data from calorimetric and spectroscopic measurements

New data concerning glaucophane are presented. New high temperature drop calorimetry data from 400 to 800 K are used to constrain the heat capacity at high temperature. Unpublished low temperature calorimetric data are used to estimate entropy up to 900 K. These data, corrected for composition, are fitted for C _p and S to the polynomial expressions (J · mol^?1 · K^?2) for T> 298.15 K: $$\begin{gathered} C_p = 11.4209 * 10^2 - 40.3212 * 10^2 /T^{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-\nulldelimiterspace} 2}} - 41.00068 * 10^6 /T^2 \hfill \\ + 52.1113 * 10^8 /T^3 \hfill \\ \end{gathered} $$ $$\begin{gathered} S = 539 + 11.4209 * 10^2 * \left( {\ln T - \ln 298.15} \right) - 80.6424 * 10^2 \hfill \\ * \left( {T^{ - {1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-\nulldelimiterspace} 2}} - 1/\left( {298.15} \right)^{{1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-\nulldelimiterspace} 2}} } \right) + 20.50034 * 10^6 \hfill \\ * \left( {T^{ - 2} - 1/\left( {298.15} \right)^2 } \right) - 17.3704 * 10^8 * \left( {T^{ - 3} - \left( {1/298.15} \right)^3 } \right) \hfill \\ \end{gathered} $$ IR and Raman spectra from 50 to 3600 cm^?1 obtained on glaucophane crystals close to the end member composition are also presented. These spectroscopic data are used with other data (thermal expansion, acoustic velocities etc.) in vibrational modelling. This last method provides an independent way for the determination of the thermodynamic properties (C_p and entropy). The agreement between measured and calculated properties is excellent (less than 2% difference between 100 and 1000 K). It is therefore expected that vibrational modelling could be applied to other amphiboles for which spectroscopic data are available. Finally, the enthalpy of formation of glaucophane is calculated.     

Dislocation-assisted diffusion of oxygen in albite

Oxygen diffusion in albite has been determined by the integrating (bulk ¹⁸O) method between 750° and 450° C, for a P _H2O of 2 kb. The original material has a low dislocation density (<10⁶ cm^?2), and its lattice diffusion coefficient (D ₁), given below, agrees well with previous determinations. A sample was deformed at high temperature and pressure to produce a uniform dislocation density of 5 × 10⁹ cm^?2. The diffusion coefficient (D _a) for this deformed material, given below, is about 0.5 and 0.7 orders of magnitude larger than D ₁ at 700° and 450° C, respectively. This enhancement is believed due to faster diffusion along the cores of dislocations. Assuming a dislocation core radius of 4 Å, the calculated pipe diffusion coefficient (D _p), given below, is about 5 orders of magnitude larger than D ₁. These results suggest that volume diffusion at metamorphic conditions may be only slightly enhanced by the presence of dislocations. $$\begin{gathered} D_1 = 9.8 \pm 6.9 \times 10^{ - 6} (cm^2 /\sec ) \hfill \\ {\text{ }} \cdot \exp [ - 33.4 \pm 0.6(kcal/mole)/RT] \hfill \\ \end{gathered} $$ $$\begin{gathered} D_a = 7.6 \pm 4.0 \times 10^{ - 6} (cm^2 /\sec ) \hfill \\ {\text{ }} \cdot \exp [ - 30.9 \pm 1.1(kcal/mole)/RT] \hfill \\ \end{gathered} $$ $$\begin{gathered} D_p \approx 1.2 \times 10^{ - 1} (cm^2 /\sec ) \hfill \\ {\text{ }} \cdot \exp [ - 29.8(kcal/mole)/RT]. \hfill \\ \end{gathered} $$