The Fe-rich liquidus in the Fe-FeS system from 1 bar to 10 GPa

The composition and evolution of a metallic planetary core is determined by the behavior with pressure of the eutectic and the liquidus on the Fe-rich side of the Fe-FeS eutectic. New experiments at 6 GPa presented here, along with existing experimental data, inform a thermodynamic model for this liquidus from 1 bar to at least 10 GPa. Fe-FeS has a eutectic that becomes more Fe-rich but remains constant in T up to 6 GPa. The 1 bar, 3 GPa, and 6 GPa liquidi all cross at a pivot point at 1640 ± 5 K and FeS_37 ± 0.5. This liquid/crystalline metal equilibrium is T-x-fixed and pressure independent through 6 GPa. Models of the 1 bar through 10 GPa experimental liquidi show that with increasing P there is an increase in the T separation between the liquidus and the crest of the metastable two-liquid solvus. The solvus crest decreases in T with increasing P. The model accurately reproduces all the experimental liquidi from 1 bar to 10 GPa, as well as reproducing the 0-6 GPa pivot point. The 14 GPa experimental liquidus ( [Chen et al., 2008a] and Chen et al., 2008b) deviates sharply from the lower pressure trends indicating that the 0-10 GPa model no longer applies to this 14 GPa data.     

No mushy zones in the Earth’s core

Mushy zones, assemblages of crystals and their pore-space liquids, have been invoked for both the upper and lower boundaries of the liquid outer core. The timescale of very slow accumulation compared with solidification at either of these interfaces militates against such zones, where instead hard ground should be expected to form by solidification at the interface. Such adcumulus growth involves isothermal, isocompositional solidification by successful exchange of evolving solute with fresh melt from an infinite reservoir. At both boundaries of the outer core, the removal of rejected material is significantly aided by compositional convection. The accumulation rates at the outer core boundaries are orders of magnitude slower than required for adcumulus growth, as calibrated both by field and experimental evidence in silicate melts. A conceptual phase diagram for the core-mantle boundary helps to visualize the relevant equilibria. Capture of core metal into the mantle has been suggested to occur via a mushy zone, to explain a high electrical conductivity there, as plausibly required by the secular behavior of the Earth’s nutation. One conjecture is that the rejected light elements from the freezing of the inner core might be able to congregate as a porous flotation sediment at the top of the core. The idea of porosity in such a mushy zone must be rejected from experience with solidification of cumulates from magmas.A high electrical conductivity might instead be caused by solution of core metal by mantle, followed by exsolution. The hottest part of the mantle lies in contact with the molten outer core, where the maximum solubility of Fe must occur in the major mantle phases. On leaving the core-mantle boundary, the mantle must cool and may exsolve metal on the metal-silicate solvus. If the iron-rich metal resides chiefly in the rheologically weaker metal oxide phase, which coats the deforming perovskite grains, it may furnish a short circuit for mantle conductivity in the basal mantle. At still cooler and higher levels, the mantle encounters more normal mantle redox conditions, and any exsolved Fe metal should oxidize to FeO in the metal oxide and perovskite phases, ceasing to be a conductor.     

Experimental study of platinum solubility in silicate melt to 14 GPa and 2273 K: Implications for accretion and core formation in Earth

We determined the solubility limit of Pt in molten haplo-basalt (1 atm anorthite-diopside eutectic composition) in piston-cylinder and multi-anvil experiments at pressures between 0.5 and 14 GPa and temperatures from 1698 to 2223 K. Experiments were internally buffered at ∼IW + 1. Pt concentrations in quenched-glass samples were measured by laser-ablation inductively coupled-plasma mass spectrometry (LA-ICPMS). This technique allows detection of small-scale heterogeneities in the run products while supplying three-dimensional information about the distribution of Pt in the glass samples. Analytical variations in ¹⁹⁵Pt indicate that all experiments contain Pt nanonuggets after quenching. Averages of multiple, time-integrated spot analyses (corresponding to bulk analyses) typically have large standard deviations, and calculated Pt solubilities in silicate melt exhibit no statistically significant covariance with temperature or pressure. In contrast, averages of minimum ¹⁹⁵Pt signal levels show less inter-spot variation, and solubility shows significant covariance with pressure and temperature. We interpret these results to mean that nanonuggets are not quench particles, that is, they were not dissolved in the silicate melt, but were part of the equilibrium metal assemblage at run conditions. We assume that the average of minimum measured Pt abundances in multiple probe spots is representative of the actual solubility. The metal/silicate partition coefficients (D^met/sil) is the inverse of solubility, and we parameterize D^met/sil in the data set by multivariate regression. The statistically robust regression shows that increasing both pressure and temperature causes D^met/silto decrease, that is, Pt becomes more soluble in silicate melt. D^met/sil decreases by less than an order of magnitude at constant temperature from 1 to 14 GPa, whereas isobaric increase in temperature produces a more dramatic effect, with D^met/sil decreasing by more than one order of magnitude between 1623 and 2223 K. The Pt abundance in the Earth’s mantle requires that D^met/sil is ∼1000 assuming core-mantle equilibration. Geochemical models for core formation in Earth based on moderately and slightly siderophile elements are generally consistent with equilibrium metal segregation at conditions generally in the range of 20-60 GPa and 2000-4000 K. Model extrapolations to these conditions show that the Pt abundance of the mantle can only be matched if oxygen fugacity is high (∼IW) and if Pt mixes ideally in molten iron, both very unlikely conditions. For more realistic values of oxygen fugacity (∼IW − 2) and experimentally-based constraints on non-ideal mixing, models show that D^met/sil would be several orders of magnitude too high even at the most favorable conditions of pressure and temperature. These results suggest that the mantle Pt budget, and by implication other highly siderophile elements, was added by late addition of a ‘late veneer’ phase to the accreting proto-Earth.     

Experimental study and modeling of fluid reaction paths in the quartz-kyanite ± muscovite-water system at 0.7 GPa in the 350-550 °C range: Implications for Al selective transfer during metamorphism

In order to quantify Al transfer in response to fluid-mineral equilibration under evolving metamorphic conditions, isobaric (0.7 GPa) experiments were conducted in the 350-550 °C range. Disequilibrium was induced (1) by holding initially pure water and natural minerals (kyanite + quartz ± muscovite enclosed in a perforated inner capsule) under isothermal conditions and (2) by stepwise temperature variations. In all experiments, secondary Al-bearing phases crystallized in the external tube of a “tube-in-tube” setup (SEM characterization); they are interpreted as witnesses of the evolution of the fluid composition (fluid reaction path). These reaction paths and the subsequent amount of secondary crystallizations were modeled using thermodynamic data from SUPCRT92 and estimates of both starting-mineral dissolution rates and elemental diffusion coefficients from the literature. A major result is that the amount of aluminum transferred to secondary phases is a thousand times larger than the calculated Al concentration in the fluid. Although the crystallization of Al-bearing phases was expected as a response to a temperature decrease, the stepwise temperature increase (20 °C/day) also led to aluminum transfer towards secondary phases. In the course of re-equilibration, the fluid first becomes saturated with respect to aluminosilicates and then reaches silica saturation, due to the low solubility of Al-minerals. Consequently, aluminosilicates partly recrystallize in response to a temperature increase. Crystallization of secondary Al-phases in the external tube implies that aqueous aluminum was efficiently transported from the inner capsule, even in the pure Al₂O₃-SiO₂-H₂O system. Therefore, mass balance calculations considering a constant Al reference frame, i.e., postulating Al immobility, should be regarded with caution.     

Experimental determination of the solubility of natural wollastonite in pure water up to pressures of 5 GPa and at temperatures of 400-800 °C

The solubility of natural, near-end-member wollastonite-I (>99.5% CaSiO₃) has been determined at temperatures from 400 to 800 °C and pressures between 0.8 and 5 GPa in piston-cylinder apparatus with the weight-loss method. Chemical analysis of quench products and optical monitoring in a hydrothermal diamond anvil cell demonstrates that no additional phases form during dissolution. Wollastonite-I, therefore, dissolves congruently in the pressure-temperature range investigated. The solubility of CaSiO₃ varies between 0.175 and 13.485 wt% and increases systematically with both temperature and pressure up to 3.0 GPa. Above 3.0 GPa wollastonite-I reacts rapidly to the high-pressure modification wollastonite-II. No obvious trends are evident in the solubility of wollastonite-II, with values between 1.93 and 10.61 wt%. The systematics of wollastonite-I solubility can be described well by a composite polynomial expression that leads to isothermal linear correlation with the density of water. The molality of dissolved wollastonite-I in pure water is then

log(m_woll)=2.2288-3418.23×T^-1+671386.84×T^-2+logρ_H2O×(5.4578+2359.11×T^-1).     

Systematics of metal-silicate partitioning for many siderophile elements applied to Earth’s core formation

Superliquidus metal-silicate partitioning was investigated for a number of moderately siderophile (Mo, As, Ge, W, P, Ni, Co), slightly siderophile (Zn, Ga, Mn, V, Cr) and refractory lithophile (Nb, Ta) elements. To provide independent constrains on the effects of temperature, oxygen fugacity and silicate melt composition, isobaric (3 GPa) experiments were conducted in piston cylinder apparatus at temperature between 1600 and 2600 °C, relative oxygen fugacities of IW−1.5 to IW−3.5, and for silicate melt compositions ranging from basalt to peridotite. The effect of pressure was investigated through a combination of piston cylinder and multi-anvil isothermal experiments between 0.5 and 18 GPa at 1900 °C. Oxidation states of siderophile elements in the silicate melt as well as effect of carbon saturation on partitioning are also derived from these results. For some elements (e.g. Ga, Ge, W, V, Zn) the observed temperature dependence does not define trends parallel to those modeled using metal-metal oxide free energy data. We correct partitioning data for solute interactions in the metallic liquid and provide a parameterization utilized in extrapolating these results to the P-T-X conditions proposed by various core formation models. A single-stage core formation model reproduces the mantle abundances of several siderophile elements (Ni, Co, Cr, Mn, Mo, W, Zn) for core-mantle equilibration at pressures from 32 to 42 GPa along the solidus of a deep peridotitic magma ocean (∼3000 K for this pressure range) and oxygen fugacities relevant to the FeO content of the present-day mantle. However, these P-T-fO₂ conditions cannot produce the observed concentrations of Ga, Ge, V, Nb, As and P. For more reducing conditions, the P-T solution domain for single stage core formation occurs at subsolidus conditions and still cannot account for the abundances of Ge, Nb and P. Continuous core formation at the base of a magma ocean at P-T conditions constrained by the peridotite liquidus and fixed fO₂ yields concentrations matching observed values for Ni, Co, Cr, Zn, Mn and W but underestimates the core/mantle partitioning observed for other elements, notably V, which can be reconciled if accretion began under reducing conditions with progressive oxidation to fO₂ conditions consistent with the current concentration of FeO in the mantle as proposed by Wade and Wood (2005). However, neither oxygen fugacity path is capable of accounting for the depletions of Ga and Ge in the Earth’s mantle. To better understand core formation, we need further tests integrating the currently poorly-known effects of light elements and more complex conditions of accretion and differentiation such as giant impacts and incomplete equilibration.     

Experimental investigation of the solubility of albite and jadeite in H2O, with paragonite + quartz at 500 and 600 °C, and 1-2.25 GPa

The solubilities of the assemblages albite + paragonite + quartz and jadeite + paragonite + quartz in H₂O were determined at 500 and 600 °C, 1.0-2.25 GPa, using hydrothermal piston-cylinder methods. The three minerals are isobarically and isothermally invariant in the presence of H₂O, so fluid composition is uniquely determined at each pressure and temperature. A phase-bracketing approach was used to achieve accurate solubility determinations. Albite + quartz and jadeite + quartz dissolve incongruently in H₂O, yielding residual paragonite which could not be retrieved and weighed. Solution composition fixed by the three-mineral assemblage at a given pressure and temperature was therefore bracketed by adding NaSi₃O_6.5 glass in successive experiments, until no paragonite was observed in run products. Solubilities derived from experiments bounding the appearance of paragonite thus constrain the equilibrium fluid composition. Results indicate that, at a given pressure, Na, Al, and Si concentrations are higher at 600 °C than at 500 °C. At both 500 and 600 °C, solubilities of all three elements increase with pressure in the albite stability field, to a maximum at the jadeite-albite-quartz equilibrium. In the jadeite stability field, element concentrations decline with continued pressure increase. At the solubility maximum, Na, Al, and Si concentrations are, respectively, 0.16, 0.05, and 0.48 molal at 500 °C, and 0.45, 0.27, and 1.56 molal at 600 °C. Bulk solubilities are 3.3 and 10.3 wt% oxides, respectively. Observed element concentrations are everywhere greater than those predicted from extrapolated thermodynamic data for simple ions, monomers, ion pairs, and the silica dimer. The measurements therefore require the presence of additional, polymerized Na-Al-Si-bearing species in the solutions. The excess solubility is >50% at all conditions, indicating that polymeric structures are the predominant solutes in the P-T region studied. The solubility patterns likely arise from combination of the large solid volume change associated with the albite-jadeite-quartz equilibrium and the rise in Na-Al-Si polymerization with approach to the hydrothermal melting curves of albite + quartz and jadeite + quartz. Our results indicate that polymerization of Na-Al-Si solutes is a fundamental aspect of fluid-rock interaction at high pressure. In addition, the data suggest that high-pressure metamorphic isograds can impose unexpected controls on metasomatic mass transfer, that significant metasomatic mass transfer prior to melting should be considered in migmatitic terranes, and that polymeric complexes may be an important transport agent in subduction zones.     

The effect of CO2 on the speciation of RbBr in solution at temperatures to 579 °C and pressures to 0.26 GPa

Carbon dioxide- and salt-bearing solutions are common in granulite, ore-forming and magmatic environments. The presence of CO₂ affects mineral solubilities, fluid miscibility, and viscosity and wetting properties, and is expected to affect salt speciation. EXAFS measurements of RbBr-H₂O-CO₂ fluids contained in corundum-osed synthetic fluid inclusions (SFLINCs) have been used to investigate the effect of CO₂ on salt speciation at temperatures to 579 °C and pressures to around 0.26 GPa.Forward modelling indicates that solute dehydration is difficult to distinguish from up to around 40% of Rb-Br ion-pairing, so results refer to the total number of nearest neighbours, which are likely to be mostly O present in waters of hydration, but may also include Br, if ion pairing is present. Additionally, results relate to the number of well-ordered neighbours in the first shell, because nearest neighbours with a high degree of disorder may be present but contribute minimally to the EXAFS signal. Analysis of the EXAFS results at the Rb edge for the CO₂-free solution is consistent with previous work and shows that the number of nearest neighbours for Rb in CO₂-free solutions decreases from 6 ± 0.6 to 1.4 ± 0.1 as temperature increases from 20 to 534 °C. The decrease is accompanied by a decrease in Rb-x bondlengths of 0.05 Å, where x is the first shell scatterer. Results for the CO₂-bearing solution are different to those for the CO₂-free solution. The number of nearest neighbours is 16 and 22% less than for the CO₂-bearing solution at 312 and 445 °C respectively. Changes in the numbers of nearest neighbours correlate well with calculated changes in the bulk solution dielectric constant; CO₂-bearing and CO₂-free solutions lie on the same trend, which suggests that it may be possible to calculate the number of nearest neighbours from dielectric constant. Rb-x bondlengths for the CO₂-bearing solution are statistically indistinguishable to those for the CO₂-free inclusions. Results for Br are worse quality than for Rb so EXAFS analysis could not be completed, however XANES spectra for CO₂-free and CO₂-bearing solutions are consistent with solute dehydration similar to that recorded by the Rb spectra. The conclusions of this study provide support for the notion that CO₂ has a fundamental effect on the mechanics of solubility, and that these effects should be incorporated into conceptual and quantitative thermodynamic models.     

A threshold in soil formation at Earth’s arid-hyperarid transition

The soils of the Atacama Desert in northern Chile have long been known to contain large quantities of unusual salts, yet the processes that form these soils are not yet fully understood. We examined the morphology and geochemistry of soils on post-Miocene fans and stream terraces along a south-to-north (27° to 24° S) rainfall transect that spans the arid to hyperarid transition (21 to ∼2 mm rain y⁻¹). Landform ages are ? 2 My based on cosmogenic radionuclide concentrations in surface boulders, and Ar isotopes in interbedded volcanic ash deposits near the driest site indicate a maximum age of 2.1 My. A chemical mass balance analysis that explicitly accounts for atmospheric additions was used to quantify net changes in mass and volume as a function of rainfall. In the arid (21 mm rain y⁻¹) soil, total mass loss to weathering of silicate alluvium and dust (−1030 kg m⁻²) is offset by net addition of salts (+170 kg m⁻²). The most hyperarid soil has accumulated 830 kg m⁻² of atmospheric salts (including 260 kg sulfate m⁻² and 90 kg chloride m⁻²), resulting in unusually high volumetric expansion (120%) for a soil of this age. The composition of both airborne particles and atmospheric deposition in passive traps indicates that the geochemistry of the driest soil reflects accumulated atmospheric influxes coupled with limited in-soil chemical transformation and loss. Long-term rates of atmospheric solute addition were derived from the ion inventories in the driest soil, divided by the landform age, and compared to measured contemporary rates. With decreasing rainfall, the soil salt inventories increase, and the retained salts are both more soluble and present at shallower depths. All soils generally exhibit vertical variation in their chemistry, suggesting slow and stochastic downward water movement, and greater climate variability over the past 2 My than is reflected in recent (∼100 y) rainfall averages. The geochemistry of these soils shows that the transition from arid to hyperarid rainfall levels marks a fundamental geochemical threshold: in wetter soils, the rate and character of chemical weathering results in net mass loss and associated volumetric collapse after 10⁵ to 10⁶ years, while continuous accumulation of atmospheric solutes in hyperarid soils over similar timescales results in dramatic volumetric expansion. The specific geochemistry of hyperarid soils is a function of atmospheric sources, and is expected to vary accordingly at other hyperarid sites. This work identifies key processes in hyperarid soil formation that are likely to be independent of location, and suggests that analogous processes may occur on Mars.     

The Earth’s tungsten budget during mantle melting and crust formation

During silicate melting on Earth, W is one of the most incompatible trace elements, similar to Th, Ba or U. As W is also moderately siderophile during metal segregation, ratios of W and the lithophile Th and U in silicate rocks have therefore been used to constrain the W abundance of the Earth’s mantle and the Hf-W age of core formation. This study presents high-precision W concentration data obtained by isotope dilution for samples covering important silicate reservoirs on Earth. The data reveal significant fractionations of W from other highly incompatible lithophile elements such as Th, U, and Ta. Many arc lavas exhibit a selective enrichment of W relative to Th, U, and Nb-Ta, reflecting W enrichment in the sub-arc mantle via fluid-like components derived from subducting plates. In contrast, during enrichment by melt-like subduction components, W is generally slightly depleted relative to Th and U, but is still enriched relative to Ta. Hence, all arc rocks and the continental crust exhibit uniformly low Ta/W (ca. 1), whereas W/Th and W/U may show opposite fractionation trends, depending on the role of fluid- and melt-like subduction components. Further high-precision W data for OIBs and MORBs reveal a systematic depletion of W in both rock types relative to other HFSE, resulting in high Ta/W that are complementary to the low Ta/W observed in arc rocks and the continental crust. Similar to previous interpretations based on Nb/U and Ce/Pb systematics, our Ta/W data confirm a depletion of the depleted upper mantle (DM) in fluid mobile elements relative to the primitive mantle (PRIMA). The abundance of W in the depleted upper mantle relative to other immobile and highly incompatible elements such as Nb and Ta is therefore not representative of the bulk silicate Earth. Based on mass balance calculations using Ta-W systematics in the major silicate reservoirs, the W abundance of the Earth’s primitive mantle can be constrained to 12 ppb, resulting in revised ratios of W-U and W-Th of 0.53 and 0.14, respectively. The newly constrained Hf-W ratio of the silicate Earth is 25.8, significantly higher than previously estimated (18.7) and overlaps within error the Hf-W ratio proposed for the Moon (ca. 24.9). The ¹⁸²Hf-¹⁸²W model age for the formation of the Earth’s core that is inferred from the ¹⁸²W abundance and the Hf/W of the silicate Earth is therefore younger than previously calculated, by up to 5 Myrs after solar system formation depending on the accretion models used. The similar Hf/W ratios and ¹⁸²W compositions of the Earth and the silicate Moon suggest a strong link between the Moon forming giant impact and final metal-silicate equilibration on the Earth.     

Adsorption of Rb and Sr at the orthoclase (0 0 1)-solution interface

Adsorption of Rb⁺ and Sr²⁺ at the orthoclase (0 0 1)-solution interface is probed with high-resolution X-ray reflectivity and resonant anomalous X-ray reflectivity. Specular X-ray reflectivity data for orthoclase in contact with 0.01 m RbCl solution at pH 5.5 reveal a systematic increase in electron density adjacent to the mineral surface with respect to that observed in contact with de-ionized water (DIW). Quantitative analysis indicates that Rb⁺ adsorbs at a height of 0.83 ± 0.03 Å with respect to the bulk K⁺ site with a nominal coverage of 0.72 ± 0.10 ions per surface unit mesh (55.7 Å²). These results are consistent with an ion-exchange reaction in which Rb⁺ occupies an inner-sphere adsorption (IS) site. In contrast, X-ray reflectivity data for orthoclase in contact with 0.01 m Sr(NO₃)₂ solution at pH 5.3 reveal few significant changes with respect to DIW. Resonant anomalous X-ray reflectivity was used to probe Sr²⁺ adsorption and to image its vertical distribution. This element-specific measurement reveals that Sr²⁺ adsorbs with a total coverage of 0.37 ± 0.02 ions per surface unit mesh, at a substantially larger height (3.28 ± 0.05 Å) than found for Rb⁺, and with a relatively broad density distribution (having a root-mean-square width of 1.88 ± 0.08 Å for a single-peak model), implying that Sr²⁺ adsorbs primarily as a fully-hydrated outer-sphere (OS), species. Comparison to a two-height model suggests that 13 ± 5% of the adsorbed Sr²⁺ may be present as an IS species. This partitioning implies a ∼5 kJ/mol difference in free energy between the IS and OS Sr²⁺ on orthoclase. Differences in the partitioning of Sr²⁺ between IS and OS species for orthoclase (0 0 1) and muscovite (0 0 1) suggest control by the geometry of the IS adsorption site. Results for the OS distribution are compared to predictions of the Poisson-Boltzmann equation in the strong coupling regime, which predicts an intrinsically narrow vertical diffuse ion distribution; the OS distribution might thus be thought of as the diffuse ion profile in the limit of high surface charge.     

A model for C-O-H fluid in the Earth’s mantle

A model is presented for predicting the composition (H₂O, CO₂, CH₄, H₂, CO, O₂ and C₂H₆) in the C-O-H fluid system under high temperatures and pressures found in the Earth’s mantle. The model is based on a molecular dynamic equation of state, statistical mechanics calculations and non-stoichiometric global free-energy minimization. Although the model is not fitted to experimental data on C-O-H speciation, it does accurately reproduce these datasets and should extrapolate at least to the depths of ∼80-220 km. The model results suggest that (1) in the upper cratonic mantle, H₂O is the dominant fluid species in the C-O-H fluid system; (2) the abundance of CO₂ increases with decreasing depth, the trend of CH₄ is just the opposite; (3) the boundary between lithosphere and asthenosphere generally divides fluid systems into H₂O-CH₄+ minor species and H₂O-CO₂+ minor species, respectively; (4) it is entirely possible to generate methane and ethane and possibly other hydrocarbons under mantle conditions, confirming previously experimental results.     

Molecular dynamics simulations of the orthoclase (0 0 1)- and (0 1 0)-water interfaces

Molecular dynamics simulations of water in contact with the (0 0 1) and (0 1 0) surfaces of orthoclase (KAlSi₃O₈) were carried out to investigate the structure and dynamics of the feldspar-water interface, contrast the intrinsic structural properties of the two surfaces, and provide a basis for future work on the diffusion of ions and molecules in microscopic mineral fractures. Electron density profiles were computed from the molecular dynamics trajectories and compared with those derived experimentally from high-resolution X-ray reflectivity measurements by Fenter and co-workers [Fenter P., Cheng L., Park C., Zhang H. and Sturchio N. C. (2003a) Structure of the orthoclase (0 0 1)- and (0 1 0)-water interfaces by high-resolution X-ray reflectivity. Geochim. Cosmochim. Acta67, 4267-4275]. For each surface, three scenarios were considered whereby the interfacial species is potassium, water, or a hydronium ion. Excellent agreement was obtained for the (0 0 1) surface when potassium is the predominant interfacial species; however, some discrepancies in the position of the interfacial peaks were obtained for the (0 1 0) surface. The two surfaces showed similarities in the extent of water ordering at the interface, the activation energies for water and potassium desorption, and the adsorption localization of interfacial species. However, there are also important differences between the two surfaces in the coordination of a given adsorbed species, adsorption site densities, and the propensity for water molecules in surface cavities and those in the first hydration layer to coordinate to surface bridging oxygen atoms. These differences may have implications for the extent of dissolution in the low-pH regime since hydrolysis of Si(Al)OSi(Al) bonds is a major dissolution mechanism.     

A synthetic fluid inclusion study of copper solubility in hydrothermal brines from 525 to 725 °C and 0.3 to 1.7 GPa

Fluid inclusions were synthesized in a piston-cylinder apparatus under mineral-buffered conditions over a range of Cl concentration (0.29 to 11.3 mol kg⁻¹), temperature (525 to 725 °C), and pressure (0.3 to 1.7 GPa). All fluids were buffered by the mineral assemblage native copper + cuprite + talc + quartz. In situ fluid composition was determined by analysing individual fluid inclusions by LA-ICPMS and independently analysing the quench solution. The solubility data provide basic information necessary to model the high temperature behaviour of Cu in magmatic-hydrothermal systems. Copper concentrations up to ∼15 wt% were measured at 630 °C and 0.34 GPa. These results give an upper limit for Cu in natural fluids and support field-based observations of similar high Cu concentrations in fluids at near-magmatic conditions. Experimental evidence indicates that Cu⁺ may form neutral chloride complexes with the general stoichiometry with n up to 4, though n ? 2 is typical for the majority of the experimental conditions. At high pressure (>∼0.5 GPa) there is evidence that hydroxide species, e.g., CuOH⁰, become increasingly important and may predominate over copper(I)-chloride complexes. The roles of fluid mixing, cooling and decompression in ore-forming environments are also discussed.     

Speciation of reduced C-O-H volatiles in coexisting fluids and silicate melts determined in-situ to ∼1.4 GPa and 800 °C

The structure of silicate melts in the system Na₂O·4SiO₂ saturated with reduced C-O-H volatile components and of coexisting silicate-saturated C-O-H solutions has been determined in a hydrothermal diamond anvil cell (HDAC) by using confocal microRaman and FTIR spectroscopy as structural probes. The experiments were conducted in-situ with the melt and fluid at high temperature (up to 800 °C) and pressure (up to 1435 MPa). Redox conditions in the HDAC were controlled with the reaction, Mo + H₂O = MoO₂ + H₂, which is slightly more reducing than the Fe + H₂O = FeO + H₂ buffer at 800 °C and less.The dominant species in the fluid are CH₄ + H₂O together with minor amounts of molecular H₂ and an undersaturated hydrocarbon species. In coexisting melt, CH₃ - groups linked to the silicate melt structure via Si-O-CH₃ bonding may dominate and possibly coexists with molecular CH₄. The abundance ratio of CH₃ - groups in melts relative to CH₄ in fluids increases from 0.01 to 0.07 between 500 and 800 °C. Carbon-bearing species in melts were not detected at temperatures and pressures below 400 °C and 730 MPa, respectively. A schematic solution mechanism is, Si-O-Si + CH₄?Si-O-CH₃+H-O-Si. This mechanism causes depolymerization of silicate melts. Solution of reduced (C-O-H) components will, therefore, affect melt properties in a manner resembling dissolved H₂O.     

Redox evolution and mass transfer during serpentinization: An experimental and theoretical study at 200 °C, 500 bar with implications for ultramafic-hosted hydrothermal systems at Mid-Ocean Ridges

Highly reducing and high-pH vent fluids characterize moderately low temperature ultramafic-hosted hydrothermal systems, such as the recently discovered Lost City hydrothermal field at 30°N Mid-Atlantic Ridge Ridge (MAR). To better understand the role of mineral reaction rates on changes in fluid chemistry and mineralization processes in these and similar systems, we conducted an experimental study involving seawater and peridotite at 200 °C, 500 bar. Time series changes in fluid chemistry were monitored and compared with analogous data predicted using experimental and theoretical data for mineral dissolution rates. Although there was qualitative agreement between predicted and measured changes in the chemical evolution of the fluid for some species, the rate and magnitude of increase in pH, dissolved chloride and H₂ did not agree well with predictions based on theoretical modeling results. Experimental data indicate that dissolved H₂ abruptly and intermittently increased, reaching a value only approximately 20% of that predicted assuming magnetite as the primary Fe-bearing alteration phase. The distribution and valence of Fe in primary and secondary minerals reveal that the most abundant secondary mineral, serpentine, contained significant amounts of both ferric and ferrous Fe, with the less abundant brucite, also being Fe-rich (X_Fe = 0.3). Surprisingly, magnetite was present in only trace amounts, indicating that H₂ generation was largely accommodated by the formation of Fe-chrysotile. Accordingly, the diversity of Fe-bearing secondary minerals together with rates of serpentinization less than theoretically predicted, account best for the relatively low dissolved H₂ concentrations produced. Thus, the experimental data can be used to obtain provisional estimates of thermodynamic data for Fe-bearing minerals, enhancing the application of reaction path models depicting mass transfer processes during serpentinization at mid-ocean ridges. Similarly, the observed differences between theoretically predicted and experimentally measured pH values result from constraints imposed by complex patterns of mass transfer inherent to the experimental system. In particular, the experimental observation of a late stage increase in Na/Cl ratio likely results from the dissolution of a Na₂O component of clinopyroxene, which causes pH to increase sufficiently to induce precipitation of a Ca-bearing phase, perhaps portlandite. As with the redox variability observed during the experiment, this event could not be predicted, underscoring the need to use caution when modeling alteration processes in the chemically complex ultramafic-hosted hydrothermal systems at elevated temperatures and pressures.     

Water ordering and surface relaxations at the hematite (1 1 0)-water interface

Structural characterization of iron oxide-water interfaces provides insight into the mechanisms through which these minerals control contaminant fate and element cycling in soil, sedimentary, and groundwater systems. Ordering of interfacial water and structural relaxations at the hematite (1 1 0) surface have been investigated in situ using high-resolution specular X-ray reflectivity. These measurements demonstrate that relaxations are constrained to primarily the top ∼5 Å of the surface. Near-surface iron atoms do not relax substantially, although the uppermost layer displays an increased distribution width, while the undercoordinated oxygens on the surface uniformly relaxed outward. Two sites of adsorbed water and additional layering of water farther from the surface were observed. Water fully covers the (1 1 0) surface and appears to form a continuous network extending into bulk solution, with positional order decreasing to that of a disordered bulk fluid within 1 nm. The arrangement of water is similar to that on the hematite (0 1 2) surface, which has a similar surface topography, although these surfaces display different vibrational amplitudes or positional disorder of adsorbed water molecules and average spacings of near-surface layered water. Comparison between these surfaces suggests that interfacial water ordering on hematite is controlled primarily by surface structure and steric constraints and that highly ordered water is likely common to most hematite-water interfaces.     

Experimental partitioning of uranium between liquid iron sulfide and liquid silicate: Implications for radioactivity in the Earth’s core

Measurable uranium (U) is found in metal sulfide liquids in equilibrium with molten silicate at conditions appropriate for a planetary magma ocean: 1-10 GPa, 1750-2100 °C, 0-28 wt% S, and fO₂ 2 log units below IW. However, the transfer of U from metal sulfide to silicate under our experimental conditions is so complete that insufficient U would remain so as to be of any importance to the core’s heat budget. U content in the sulfide phase increases strongly with S content but shows no significant variability with either pressure or temperature. Maximum is 0.001 while most values are considerably lower.     

Gold and copper partitioning in magmatic-hydrothermal systems at 800 °C and 100 MPa

Porphyry-type ore deposits sometimes contain fluid inclusion compositions consistent with the partitioning of copper and gold into vapor relative to coexisting brine at the depositional stage. However, this has not been reproduced experimentally at magmatic conditions. In an attempt to determine the conditions under which copper and gold may partition preferentially into vapor relative to brine at temperatures above the solidus of granitic magmas, we performed experiments at 800 °C, 100 MPa, oxygen fugacity () buffered by Ni-NiO, and fixed at either 3.5 × 10⁻² by using intermediate solid solution-pyrrhotite, or 1.2 × 10⁻⁴ by using intermediate solid solution-pyrrhotite-bornite. The coexisting vapor (∼3 wt.% NaCl eq.) and brine (∼68 wt.% NaCl eq.) were composed initially of NaCl + KCl + HCl + H₂O, with starting HCl set to <1000 μg/g in the aqueous mixture. Synthetic vapor and brine fluid inclusions were trapped at run conditions and subsequently analyzed by laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS). Our experiments demonstrate that copper and gold partitioned strongly into the magmatic volatile phase(s) (MVP) (i.e., vapor or brine) relative to a silicate melt over the entire imposed range of . Nernst style partition coefficients between coexisting brine (b) and melt (m), D^b/m (±1σ), range from 3.6(±2.2) × 10¹ to 4(±2) × 10² for copper and from 1.2(±0.6) × 10² to 2.4(±2.4) × 10³ for gold. Partition coefficients between coexisting vapor (v) and melt, D^v/m range from 2.1 ± 0.7 to 18 ± 5 and 7(±3) × 10¹ to 1.6(±1.6) × 10² for copper and gold, respectively. Partition coefficients for all experiments between coexisting brine and vapor, D^b/v (±1σ), range from 7(±2) to 1.0(±0.4) × 10² and 1.7(±0.2) to 15(±2) for copper and gold, respectively. Observed average D^b/v at an of 1.2 × 10⁻⁴ were elevated, 95(±5) and 15 ± 1 for copper and gold, respectively, relative to those at the higher of 3.5 × 10⁻² where D^b/v were 10(±5) for copper and 7(±6) for gold. Thus, there is an inverse relationship between the and the D^b/v for both copper and gold with increasing resulting in a decrease in the D^b/v signifying increased importance of the vapor phase for copper and gold transport. This suggests that copper and gold may complex with volatile S-species as well as Cl-species at magmatic conditions, however, none of the experiments of our study at 800 °C and 100 MPa had a D^b/v ? 1. We did not directly determine speciation, but infer the existence of some metal-sulfur complexes based on the reported data. We suggest that copper and gold partition preferentially into the brine in most instances at or above the wet solidus. However, in most systems, the mass of vapor is greater than the mass of brine, and vapor transport of copper and gold may become more important in the magmatic environment at higher , lower , or near the critical point in a salt-water system. A D^b/v ? 1 at subsolidus hydrothermal conditions may also occur in response to changes in temperature, , , and/or acidity.Additionally, both copper and gold were observed to partition into intermediate solid solution and bornite much more strongly than into vapor, brine or silicate melt. This suggests that, although vapor and brine are both efficient at removing copper and gold from a silicate melt, the presence of Cu-Fe sulfides can sequester a substantial portion of the copper and gold contained within a silicate melt if the Cu-Fe sulfides are abundant.     

In situ Atomic Force Microscopy study of dissolution of the barite (0 0 1) surface in water at 30 °C

The dissolution behavior of the barite (0 0 1) surface in pure water at 30 °C was investigated using in situ Atomic Force Microscopy (AFM), to better understand the dissolution mechanism and the microtopographical changes that occur during the dissolution, such as steps and etch pits. The dissolution of the barite (0 0 1) surface started with the slow retreat of steps, after which, about 60 min later, the <hk0> steps of one unit cell layer or multi-layers became two-step fronts (fast “f” and slow “s” steps) with a half-unit cell layer showing different retreat rates. The “f” step had a fast retreat rate (≈(14 ± 1) × 10⁻² nm/s) and tended to have a jagged step edge, whereas the “s” step (≈(1.8 ± 0.1) × 10⁻² nm/s) had a relatively straight front. The formation of the “f” steps led to the formation of a new one-layer step, where the front of the “s” step was overtaken by that of the immediate underlying “f” step. The “f” steps also led to the decrease of the <hk0> steps and the increase in the percentage of stable steps parallel to the [0 1 0] direction during the dissolution.Etch pits, which could be observed after about 90 min, were of three types: triangular etch pits with a depth of a half-unit cell, shallow etch pits, and deep etch pits. The triangular etch pits were bounded by the step edges parallel to [0 1 0], [1 2 0], and [] and had opposite orientations in the upper half and lower half layers. Shallow etch pits that had a depth of two or more half-unit cell layers had any two consecutive pits pointing in the opposite direction of each other. The triangular etch pit appeared to grow by simultaneously removal of a row of ions parallel to each direction from the three step edges. At first, deep etch pits were elongated in the [0 1 0] direction with a curved outline and then gradually developed to an angular form bounded by the {1 0 0}, {3 1 0}, and (0 0 1) faces. The retreat rate of the (0 0 1) face was much slower than those of the {1 0 0} and {3 1 0} and tended to separate into two rates ((0.13 ± 0.01) × 10⁻² nm/s for the deep etch pits derived from a screw dislocation and (0.07 ± 0.01) × 10⁻² nm/s for those from other line defects).The changes in the dissolution rate of a barite (0 0 1) surface during the dissolution were estimated using the retreat rates and densities of the various steps as well as the growth rates, density, and areas of the lateral faces of the deep etch pits that were obtained from this AFM analysis. Our results showed that the dissolution rate of the barite (0 0 1) surface gradually increased and approached the bulk dissolution rate because of the change in the main factor determining the dissolution rate from the density of the steps to the growth and the density of the deep etch pits on the surface.