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.     

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.     

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.     

The crystal structure of waylandite from Wheal Remfry, Cornwall, United Kingdom

Sr- and Ca-rich waylandite, $ {\left( {{\hbox{B}}{{\hbox{i}}_{0.{54}}}{\hbox{S}}{{\hbox{r}}_{0.{31}}}{\hbox{C}}{{\hbox{a}}_{0.{25}}}{{\hbox{K}}_{0.0{1}}}{\hbox{B}}{{\hbox{a}}_{0.0{1}}}} \right)_{\Sigma 1.12}}{{\hbox{H}}_{0.{18}}}{\left( {{\hbox{A}}{{\hbox{l}}_{{2}.{96}}}{\hbox{C}}{{\hbox{u}}_{0.0{2}}}} \right)_{\Sigma 2.98}}{\left[ {{{\left( {{{\hbox{P}}_{0.{97}}}{{\hbox{S}}_{0.0{3}}}{\hbox{S}}{{\hbox{i}}_{0.0{1}}}} \right)}_{\Sigma 1.00}}{{\hbox{O}}_4}} \right]_2}{\left( {\hbox{OH}} \right)_6} $ , from Wheal Remfry, Cornwall, United Kingdom has been investigated by single-crystal X-ray diffraction and electron microprobe analyses. Waylandite crystallises in space group R $ \overline 3 $ ? m, with the cell parameters: a?=?7.0059(7) Å, c?=?16.3431(12) Å and V?=?694.69(11) ų. The crystal structure has been refined to R ₁?=?3.76%. Waylandite has an alunite-type structure comprised of a rhombohedral stacking of (001) composite layers of corner-shared AlO₆ octahedra and PO₄ tetrahedra, with (Bi,Sr,Ca) atoms occupying icosahedrally coordinated sites between the layers.     

Semi-empirical model for retrieval of soil moisture using RISAT-1 C-Band SAR data over a sub-tropical semi-arid area of Rewari district,Haryana (India)

We have estimated soil moisture (SM) by using circular horizontal polarization backscattering coefficient (\(\sigma ^{\mathrm{o}}_{\mathrm{RH}}\)), differences of circular vertical and horizontal \(\sigma ^{\mathrm{o}} \, (\sigma ^{\mathrm{o}}_{\mathrm{RV}} {-} \sigma ^{\mathrm{o}}_{\mathrm{RH}})\) from FRS-1 data of Radar Imaging Satellite (RISAT-1) and surface roughness in terms of RMS height (\({\hbox {RMS}}_{\mathrm{height}}\)). We examined the performance of FRS-1 in retrieving SM under wheat crop at tillering stage. Results revealed that it is possible to develop a good semi-empirical model (SEM) to estimate SM of the upper soil layer using RISAT-1 SAR data rather than using existing empirical model based on only single parameter, i.e., \(\sigma ^{\mathrm{o}}\). Near surface SM measurements were related to \(\sigma ^{\mathrm{o}}_{\mathrm{RH}}\), \(\sigma ^{\mathrm{o}}_{\mathrm{RV}} {-} \sigma ^{\mathrm{o}}_{\mathrm{RH}}\) derived using 5.35 GHz (C-band) image of RISAT-1 and \({\hbox {RMS}}_{\mathrm{height}}\). The roughness component derived in terms of \({\hbox {RMS}}_{\mathrm{height}}\) showed a good positive correlation with \(\sigma ^{\mathrm{o}}_{\mathrm{RV}} {-} \sigma ^{\mathrm{o}}_{\mathrm{RH}} \, (R^{2} = 0.65)\). By considering all the major influencing factors (\(\sigma ^{\mathrm{o}}_{\mathrm{RH}}\), \(\sigma ^{\mathrm{o}}_{\mathrm{RV}} {-} \sigma ^{\mathrm{o}}_{\mathrm{RH}}\), and \({\hbox {RMS}}_{\mathrm{height}}\)), an SEM was developed where SM (volumetric) predicted values depend on \(\sigma ^{\mathrm{o}}_{\mathrm{RH}}\), \(\sigma ^{\mathrm{o}}_{\mathrm{RV}} {-} \sigma ^{\mathrm{o}}_{\mathrm{RH}}\), and \({\hbox {RMS}}_{\mathrm{height}}\). This SEM showed \(R^{2}\) of 0.87 and adjusted \(R^{2}\) of 0.85, multiple R=0.94 and with standard error of 0.05 at 95% confidence level. Validation of the SM derived from semi-empirical model with observed measurement (\({\hbox {SM}}_{\mathrm{Observed}}\)) showed root mean square error (RMSE) = 0.06, relative-RMSE (R-RMSE) = 0.18, mean absolute error (MAE) = 0.04, normalized RMSE (NRMSE) = 0.17, Nash–Sutcliffe efficiency (NSE) = 0.91 (\({\approx } 1\)), index of agreement (d) = 1, coefficient of determination \((R^{2}) = 0.87\), mean bias error (MBE) = 0.04, standard error of estimate (SEE) = 0.10, volume error (VE) = 0.15, variance of the distribution of differences \(({\hbox {S}}_{\mathrm{d}}^{2}) = 0.004\). The developed SEM showed better performance in estimating SM than Topp empirical model which is based only on \(\sigma ^{\mathrm{o}}\). By using the developed SEM, top soil SM can be estimated with low mean absolute percent error (MAPE) = 1.39 and can be used for operational applications.     

Low dispersive modeling of Rayleigh waves on partly staggered grids

In elastic media, finite-difference (FD) implementations of free-surface (FS) boundary conditions on partly staggered grid (PSG) use the highly dispersive vacuum formulation (VPSG). The FS boundary is embedded into a “vacuum” grid layer (null Lame’s constants and negligible density values) where the discretized equations of motion allow computing surface displacements. We place a new set of compound (stress-displacement) nodes along a planar FS and use unilateral mimetic FD discretization of the zero-traction conditions for displacement computation (MPSG). At interior nodes, MPSG reduces to standard VPSG methods and applies fourth-order centered FD along cell diagonals for staggered differentiation combined with nodal second-order FD in time. We perform a dispersion analysis of these methods on a Lamb’s problem and estimate dispersion curves from the phase difference of windowed numerical Rayleigh pulses at two FS receivers. For a given grid sampling criterion (e.g., six or ten nodes per reference S wavelength λ ^S), MPSG dispersion errors are only a quarter of the VPSG method. We also quantify root-mean-square (RMS) misfits of numerical time series relative to analytical waveforms. MPSG RMS misfits barely exceed 10 % when nine nodes sample the minimum S wavelength $\lambda _{\text {MIN}}^{\mathrm {S}}$ in transit (along distances $\sim $ 145 $\lambda _{\text {MIN}}^{\mathrm {S}}$ ). In same tests, VPSG RMS misfits exceed 70 %. We additionally compare MPSG to a consistently fourth-order mimetic method designed on a standard staggered grid. The latter equates the former’s dispersion errors on grids twice denser and shows higher RMS precision only on grids with six or less nodes per $\lambda _{\text {MIN}}^{\mathrm {S}}$ .     

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.     

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.     

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.     

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.     

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

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.     

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.     

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.     

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.     

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.     

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

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

Experimental determination of the temperature and pressure dependence of the Fe-Mg partition coefficient for coexisting garnet and clinopyroxene

An experimental study initiated to calibrate the distribution coefficient \(K_D = \frac{{({\text{FeO}}/{\text{MgO}})_{{\text{ga}}} }}{{{\text{(FeO}}/{\text{MgO)}}_{{\text{cpx}}} }}\) in eclogites as a geothermometer has been done on (a) a mineral mis, (b) a glass of the typical tholeiite composition and (c) a series of glasses of tholeiite compositions with \(6.2 < \frac{{100{\text{Mg}}}}{{{\text{Mg}} + {\text{Fe}}^{ + + } }} < 93.\) The mineral mix was found to be unsuitable as reactant due to incomplete equilibration but the minimum K _D of the mineral mix and the K _D from glass of tholeiite composition are identical within experimental uncertainty. These data constitute a reversal of the garnet/clinopyroxene partition relationship and provide justification of the use of glass as a reactant for the project. To eliminate any uncertainty in interpretation of mineral analyses due to possible variation in Fe⁺⁺⁺/Fe⁺⁺ between runs, experiments were carried out in iron capsules on the nine tholeiite glasses, thus maintaining iron as Fe⁺⁺. Microprobe analytical techniques yielded mineral analyses of comparable accuracy to analyses of natural phases for experiments within the temperature range from 600° C to 1500° C and a pressure range from 20 kb to 40 kb. It has been shown that for \(6.2 < \frac{{100{\text{Mg}}}}{{{\text{Mg}} + {\text{Fe}}^{ + + } }} < 85\) , the bulk chemical composition does not perceptibly affect the K _D value. At 30 kb the K _D value ranges from 18.0 at 600° C to 1.45 at 1400° C, defining the linear relationship in a ln K _D vs 1/T(°K) plot. The pressure dependence of the K _D -value has been shown to be greater than previously predicted. There is a straight line relationship in ln K _D vs Pressure (Kb) between 20 and 40 kb at constant temperature (1100°C). This enables us to determine K _D =fn (T, P) and \(T(^\circ {\text{K}}) = \frac{{3686 + 28.35 \times P({\text{Kb}})}}{{\ln K_D + 2.33}}\) . This expression uniquely determines the temperature of equilibration of natural eclogites of basaltic bulk composition when the K _D ^ga,cpx is known and a pressure estimate can be given.