Influence of temperature,composition, silica activity and oxygen fugacity on the H<Subscript>2</Subscript>O storage capacity of olivine at 8 GPa

Olivine crystals were grown in the presence of a hydrous silicate fluid during multi-anvil experiments at 8 GPa and 1,000–1,600°C. Experiments were conducted both in a simple system (FeO–MgO–SiO₂–H₂O) and in a more complex system containing additional elements (CaO–Na₂O–Al₂O₃–Cr₂O₃–TiO₂–FeO–MgO–SiO₂–H₂O). Silica activity was buffered by the presence of either pyroxene (high a _SiO2) or ferropericlase (low a _SiO2), and was buffered by the presence of Ni + NiO or Fe + FeO, or constrained by the presence of Fe₂O₃. Raman spectroscopy was used to identify pyroxene polymorphs in the run products. Clinoenstatite was present in the 1,000°C experiment, and enstatite in experiments at 1,400–1,520°C. The H₂O content of olivine was measured using secondary ion mass spectroscopy, and infrared spectroscopy was used to investigate the nature of hydrous defects. The H₂O storage capacity of olivine decreases with increasing temperature at 8 GPa. In contrast to previous experimental results at ≤2 GPa, no significant effect of varying oxygen fugacity is evident, but H₂O storage capacity is enhanced under conditions of low silica activity. No significant growth of low wavenumber (<3,400 cm⁻¹) peaks, generally associated with high at low pressure, was observed in the FTIR spectra of olivine from the high experiments. Our experiments show that previous high pressure H₂O storage capacity measurements for olivine synthesized under more oxidizing conditions than the Earth’s mantle are not likely to be compromised by the of the experiments. However, the considerable effect of temperature on H₂O storage capacity in olivine must be taken into account to avoid overestimation of the bulk upper mantle H₂O storage capacity.     

THERIA_G: a software program to numerically model prograde garnet growth

We present the software program THERIA_G, which allows for numerical simulation of garnet growth in a given volume of rock along any pressure–temperature–time (P–T–t) path. THERIA_G assumes thermodynamic equilibrium between the garnet rim and the rock matrix during growth and accounts for component fractionation associated with garnet formation as well as for intracrystalline diffusion within garnet. In addition, THERIA_G keeps track of changes in the equilibrium phase relations, which occur during garnet growth along the specified P–T–t trajectory. This is accomplished by the combination of two major modules: a Gibbs free energy minimization routine is used to calculate equilibrium phase relations including the volume and composition of successive garnet growth increments as P and T and the effective bulk rock composition change. With the second module intragranular multi-component diffusion is modelled for spherical garnet geometry. THERIA_G allows to simulate the formation of an entire garnet population, the nucleation and growth history of which is specified via the garnet crystal size frequency distribution. Garnet growth simulations with THERIA_G produce compositional profiles for the garnet porphyroblasts of each size class of a population and full information on equilibrium phase assemblages for any point along the specified P–T–t trajectory. The results of garnet growth simulation can be used to infer the P–T–t path of metamorphism from the chemical zoning of garnet porphyroblasts. With a hypothetical example of garnet growth in a pelitic rock we demonstrate that it is essential for the interpretation of the chemical zoning of garnet to account for the combined effects of the thermodynamic conditions of garnet growth, the nucleation history and intracrystalline diffusion. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.

An experimental study of the grain-scale processes of peridotite melting: implications for major and trace element distribution during equilibrium and disequilibrium melting

The grain-scale processes of peridotite melting were examined at 1,340°C and 1.5 GPa using reaction couples formed by juxtaposing pre-synthesized clinopyroxenite against pre-synthesized orthopyroxenite or harzburgite in graphite and platinum-lined molybdenum capsules. Reaction between the clinopyroxene and orthopyroxene-rich aggregates produces a melt-enriched, orthopyroxene-free, olivine + clinopyroxene reactive boundary layer. Major and trace element abundance in clinopyroxene vary systematically across the reactive boundary layer with compositional trends similar to the published clinopyroxene core-to-rim compositional variations in the bulk lherzolite partial melting studies conducted at similar P–T conditions. The growth of the reactive boundary layer takes place at the expense of the orthopyroxenite or harzburgite and is consistent with grain-scale processes that involve dissolution, precipitation, reprecipitation, and diffusive exchange between the interstitial melt and surrounding crystals. An important consequence of dissolution–reprecipitation during crystal-melt interaction is the dramatic decrease in diffusive reequilibration time between coexisting minerals and melt. This effect is especially important for high charged, slow diffusing cations during peridotite melting and melt-rock reaction. Apparent clinopyroxene-melt partition coefficients for REE, Sr, Y, Ti, and Zr, measured from reprecipitated clinopyroxene and coexisting melt in the reactive boundary layer, approach their equilibrium values reported in the literature. Disequilibrium melting models based on volume diffusion in solid limited mechanism are likely to significantly underestimate the rates at which major and trace elements in residual minerals reequilibrate with their surrounding melt. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Oxidation of pyrite during extraction of carbonate associated sulfate

The sulfur isotopic composition of carbonate associated sulfate (CAS) has been used to investigate the geochemistry of ancient seawater sulfate. However, few studies have quantified the reliability of δ³⁴S of CAS as a seawater sulfate proxy, especially with respect to later diagenetic overprinting. Pyrite, which typically has depleted δ³⁴S values due to authigenic fractionation associated with bacterial sulfate reduction, is a common constituent of marine sedimentary rocks. The oxidation of pyrite, whether during diagenesis or sample preparation, could thus adversely influence the sulfur isotopic composition of CAS. Here, we report the results of CAS extractions using HCl and acetic acid with samples spiked with varying amounts of pyrite. The results show a very strong linear relationship between the abundance of fine-grained pyrite added to the sample and the resultant abundance and δ³⁴S value of CAS. This data represents the first unequivocal evidence that pyrite is oxidized during the CAS extraction process. Our mixing models indicate that in samples with much less than 1 wt.% pyrite and relatively high δ³⁴S_pyrite values, the isotopic offset imparted by oxidation of pyrite should be much less than ? 4‰. A wealth of literature exists on the oxidation of pyrite by Fe³⁺ and we believe this mechanism drives the oxidation of pyrite during CAS extraction, during which the oxygen used to form sulfate is taken from H₂O, not O₂. Consequently, extracting CAS under anaerobic conditions would only slow, but not halt, the oxidation of pyrite. Future studies of CAS should attempt to quantify pyrite abundance and isotopic composition.     

Controls on porphyroblast size along a regional metamorphic field gradient

Garnet-bearing schists from the Waterville Formation of south-central Maine provide an opportunity to examine the factors governing porphyroblast size over a range of metamorphic grade. Three-dimensional sizes and locations for all garnet porphyroblasts were determined for three samples along the metamorphic field gradient spanning lowest garnet through sillimanite grade, using high-resolution X-ray computed tomography. Comparison of crystal size distributions to previous data sets obtained by stereological methods for the same samples reveals significant differences in mode, mean, and shape of the distributions. Quantitative textural analysis shows that the garnets in each rock crystallized in a diffusion-controlled nucleation and growth regime. In contrast to the typical observation of a correlation between porphyroblast size and position along a metamorphic field gradient, porphyroblast size of the lowest-grade specimen is intermediate between the high- and middle-grade specimens’ sizes. Mean porphyroblast size does not correlate with peak temperatures from garnet-biotite Fe-Mg exchange thermometry, nor is post-crystallization annealing (Ostwald Ripening) required to produce the observed textures, as was previously proposed for these rocks. Robust pseudosection calculations fail to reproduce the observed garnet core compositions for two specimens, suggesting that these calc-pelites experienced metasomatism. For each of these two specimens, Monte Carlo calculations suggest potential pre-metasomatism bulk compositions that replicate garnet core compositions. Pseudosection analyses allow the estimation of the critical temperatures for garnet growth: ∼481, ∼477, and ∼485°C for the lowest-garnet-zone, middle-garnet-zone, and sillimanite-zone specimens, respectively. Porphyroblast size appears to be determined in this case by a combination of the heating rate during garnet crystallization, the critical temperature for the garnet-forming reaction and the kinetics of nucleation. Numerical simulations of thermally accelerated, diffusion-controlled nucleation, and growth for the three samples closely match measured crystal size distributions. These observations and simulations suggest that previous hypotheses linking the garnet size primarily to the temperature at the onset of porphyroblast nucleation can only partially explain the observed textures. Also important in determining porphyroblast size are the heating rate and the distribution of favorable nucleation sites.     

Surface complexation of U(VI) on goethite (α-FeOOH)

Sorption of U(VI) to goethite is a fundamental control on the mobility of uranium in soil and groundwater. Here, we investigated the sorption of U on goethite using EXAFS spectroscopy, batch sorption experiments and DFT calculations of the energetics and structures of possible surface complexes. Based on EXAFS spectra, it has previously been proposed that U(VI), as the uranyl cation , sorbs to Fe oxide hydroxide phases by forming a bidentate edge-sharing (E2) surface complex, >Fe(OH)₂UO₂(H₂O)_n. Here, we argue that this complex alone cannot account for the sorption capacity of goethite (α-FeOOH). Moreover, we show that all of the EXAFS signal attributed to the E2 complex can be accounted for by multiple scattering. We propose that the dominant surface complex in CO₂-free systems is a bidentate corner-sharing (C2) complex, (>FeOH)₂UO₂(H₂O)₃ which can form on the dominant {101} surface. However, in the presence of CO₂, we find an enhancement of UO₂ sorption at low pH and attribute this to a (>FeO)CO₂UO₂ ternary complex. With increasing pH, U(VI) desorbs by the formation of aqueous carbonate and hydroxyl complexes. However, this desorption is preceded by the formation of a second ternary surface complex (>FeOH)₂UO₂CO₃. The three proposed surface complexes, (>FeOH)₂UO₂(H₂O)₃, >FeOCO₂UO₂, and (>FeOH)₂UO₂CO₃ are consistent with EXAFS spectra. Using these complexes, we developed a surface complexation model for U on goethite with a 1-pK model for surface protonation, an extended Stern model for surface electrostatics and inclusion of all known UO₂-OH-CO₃ aqueous complexes in the current thermodynamic database. The model gives an excellent fit to our sorption experiments done in both ambient and reduced CO₂ environments at surface loadings of 0.02-2.0 wt% U.     

Low-oxygen and chemical kinetic constraints on the geochemical niche of neutrophilic iron(II) oxidizing microorganisms

Neutrophilic iron oxidizing bacteria (FeOB) must actively compete with rapid abiotic processes governing Fe(II) oxidation and as a result have adapted to primarily inhabit low-O₂ environments where they can more successfully compete with abiotic Fe(II) oxidation. The spatial distribution of these microorganisms can be observed through the chemical gradients they affect, as measured using in situ voltammetric analysis for dissolved Fe(II), Fe(III), O₂, and FeS_(aq). Field and laboratory determination of the chemical environments inhabited by the FeOB were coupled with detailed kinetic competition studies for abiotic and biotic oxidation processes using a pure culture of FeOB to quantify the geochemical niche these organisms inhabit. In gradient culture tubes, the maximum oxygen levels, which were associated with growth bands of Sideroxydans lithotrophicus (ES-1, a novel FeOB), were 15-50 μM. Kinetic measurements made on S. lithotrophicus compared biotic/abiotic (killed control) Fe oxidation rates. The biotic rate can be a significant and measurable fraction of the total Fe oxidation rate below O₂ concentrations of approximately 50 μM, but biotic Fe(II) oxidation (via the biotic/abiotic rate comparison) becomes difficult to detect at higher O₂ levels. These results are further supported by observations of conditions supporting FeOB communities in field settings. Variablity in cell densities and cellular activity as well as variations in hydrous ferrous oxide mineral quantities significantly affect the laboratory kinetic rates. The microbial habitat (or geochemical niche) where FeOB are active is thus largely controlled by the competition between abiotic and biotic kinetics, which are dependent on Fe(II) concentration, P_O2, temperature and pH in addition to the surface area of hydrous ferric oxide minerals and the cell density/activity of FeOB. Additional field and lab culture observations suggest a potentially important role for the iron-sulfide aqueous molecular cluster, FeS_(aq), in the overall cycling of iron associated with the environments these microorganisms inhabit.     

The Fe-C system at 5 GPa and implications for Earth’s core

Earth’s core may contain C, and it has been suggested that C in the core could stabilize the formation of a solid inner core composed of Fe₃C. We experimentally examined the Fe-C system at a pressure of 5 GPa and determined the Fe-C phase diagram at this pressure. In addition, we measured solid metal/liquid metal partition coefficients for 17 trace elements and examined the partitioning behavior between Fe₃C and liquid metal for 14 trace elements. Solid metal/liquid metal partition coefficients are similar to those found in one atmosphere studies, indicating that the effect of pressure to 5 GPa is negligible. All measured Fe₃C/liquid metal partition coefficients investigated are less than one, such that all trace elements prefer the C-rich liquid to Fe₃C. Fe₃C/liquid metal partition coefficients tend to decrease with decreasing atomic radii within a given period. Of particular interest, our 5 GPa Fe-C phase diagram does not show any evidence that the Fe-Fe₃C eutectic composition shifts to lower C contents with increasing pressure, which is central to the previous reasoning that the inner core may be composed of Fe₃C.     

Coupled and complex: Human-environment interaction in the Greater Yellowstone Ecosystem, USA

Complexity theory has received considerable attention over the past decade from a wide variety of disciplines. Some who write on this topic suggest that complexity theory will lead to a unifying understanding of complex phenomena; others dismiss it as a passing and disruptive fad. We suggest that for the analysis of coupled natural/human systems, the truth emerges from the middle ground. As an approach focused as much on the connections among system elements as the elements themselves, we argue that complexity theory provides a useful conceptual framework for the study of coupled natural/human systems. It is, if nothing else, a framework that leads us to ask interesting questions about, for example, sustainability, resilience, threshold events, and predictability.In this paper we attempt to demystify the ongoing discussions on complexity theory by linking its evocative and overloaded terminology to real-world processes. We illustrate how a shift in focus from system elements to connections among elements can lead to meaningful insight into human-environment interactions that might otherwise be overlooked. We ground our discussion in ongoing interdisciplinary research surrounding Yellowstone National Park’s northern elk winter range; a tightly coupled natural/human system that has been the center of debate, conflict, and compromise for more than 135 years.