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.     

Melting of the Indarch meteorite (EH4 chondrite) at 1 GPa and variable oxygen fugacity: Implications for early planetary differentiation processes

In order to derive constraints on planetary differentiation processes, and ultimately the formation of the Earth, it is required to study a variety of meteoritic materials and to investigate their melting relations and elemental partitioning at variable pressures, temperatures, and oxygen fugacities (f_O2). This study reports the first high pressure (HP) and high temperature (HT) investigation of an enstatite chondrite (Indarch). Four series of experiments exploring various f_O2 conditions have been carried out at 1 GPa in a piston-cylinder apparatus using the EH4 chondrite Indarch. We show that temperature and redox conditions have important effects on the phase equilibria of the meteorite: the solidus and liquidus temperatures of the silicate portion increase with decreasing f_O2, and the stability fields of various phases are modified. Olivine and pyroxene are stable around 1.5 log f_O2 unit below the iron-wüstite buffer (IW−1.5), whereas quartz and pyroxene is the stable assemblage under the most reducing conditions, between IW−5.0 and IW−4.0, due to reduction of the silicate. While these changes are occurring in the silicate, the metal gains Si from the silicate, (Fe, Mg, Mn, Ca, Cr)-bearing sulfides are observed at f_O2 less than IW−4, and the partitioning of Ni and Mo are both affected by the presence of Si in Fe-S-C liquids. The f_O2 has also a significant effect on the liquid metal-liquid silicate partitioning behavior of Si and S, two possible light elements in planetary cores, and of the slightly siderophile elements Cr and Mn. With decreasing f_O2, S becomes increasingly lithophile, Si becomes increasingly siderophile, and Cr and Mn both become strongly siderophile and chalcophile. The partitioning behavior of these elements places new constraints on models of core segregation for the Earth and other differentiated bodies.     

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.     

Effect of pressure, temperature, and oxygen fugacity on the metal-silicate partitioning of Te, Se, and S: Implications for earth differentiation

We have measured liquid Fe metal-liquid silicate partitioning (D_i) of tellurium, selenium, and sulfur over a range of pressure, temperature, and oxygen fugacity (1-19 GPa, 2023-2693 K, fO₂ −0.4 to −5.5 log units relative to the iron-wüstite buffer) to better assess the role of metallic melts in fractionating these elements during mantle melting and early Earth evolution. We find that metal-silicate partitioning of all three elements decreases with falling FeO activity in the silicate melt, and that the addition of 5-10 wt% S in the metal phase results in a 3-fold enhancement of both D_Te and D_Se. In general, Te, Se, and S all become more siderophile with increasing pressure, and less siderophile with increasing temperature, in agreement with previous work. In all sulfur-bearing experiments, D_Te is greater than D_Se or D_S, with the latter two being similar over a range of P and T. Parameterized results are used to estimate metal-silicate partitioning at the base of a magma ocean which deepens as accretion progresses, with the equilibration temperature fixed at the peridotite liquidus. We show that during accretion, Te behaves like a highly siderophile element, with expected core/mantle partitioning of >10⁵, in contrast to the observed core/mantle ratio of ∼100. Less extreme differences are observed for Se and S, which yielded core/mantle partitioning 100- to 10 times higher, respectively, than the observed value. Addition of ∼0.5 wt% of a meteorite component (H, EH or EL ordinary chondrite) is sufficient to raise mantle abundances to their current level and erase the original interelement fractionation of metal-silicate equilibrium.     

Concordant Rb-Sr, Sm-Nd, and Ar-Ar ages for Northwest Africa 1460: A 346 Ma old basaltic shergottite related to “lherzolitic” shergottites

Multiple lines of evidence show that the Rb-Sr, Sm-Nd, and Ar-Ar isotopic systems individually give robust crystallization ages for basaltic (or diabasic) shergottite Northwest Africa (NWA) 1460. In contrast to other shergottites, NWA 1460 exhibits minimal evidence of excess ⁴⁰Ar, thus allowing an unambiguous determination of its Ar-Ar age. The concordant Rb-Sr, Sm-Nd, and Ar-Ar results for NWA 1460 define its crystallization age to be 346 ± 17 Ma (2σ). In combination with petrographic and trace element data for this specimen and paired meteorite NWA 480, these results strongly refute the suggestion by others that the shergottites are ∼4.1 Ga old. Current crystallization and cosmic-ray exposure (CRE) age data permit identification of a maximum of nine ejection events for Martian meteorites (numbering more than 50 unpaired specimens as of 2008) and plausibly as few as five such events. Although recent high resolution imaging of the Martian surface has identified limited areas of sparsely cratered terrains, the meteorite data suggest that either these areas are representative of larger areas from which the meteorites might come, or that the cratering chronology needs recalibration. Time-averaged ⁸⁷Rb/⁸⁶Sr = 0.16 for the mantle source of the parent magma of NWA 1460/480 over the ∼4.56 Ga age of the planet is consistent with previously estimated values for bulk silicate Mars in the range 0.13-0.16, and similar to values of ∼0.18 for the “lherzolitic” shergottites. Initial ε_Nd for NWA 1460/480 at 350 ± 16 Ma ago was +10.6 ± 0.5, which implies a time-averaged ¹⁴⁷Sm/¹⁴⁴Nd of 0.217 in the Martian mantle prior to mafic melt extraction, similar to values of 0.211-0.216 for the “lherzolitic” shergottites. These time-averaged values do not imply a simple two-stage mantle/melt evolution, but must result from multiple episodes of melt extractions from the source regions. Much higher “late-stage” ε_Nd values for the depleted shergottites imply similar processes carried to a greater degree. Thus, NWA 1460/480, the “lherzolitic” shergottites and perhaps EET 79001 give the best (albeit imperfect) estimate of the Sr- and Nd-isotopic characteristics of bulk silicate Mars.     

Ammonium in aqueous fluids to 600 °C, 1.3 GPa: A spectroscopic study on the effects on fluid properties, silica solubility, and K-feldspar to muscovite reactions

The behavior of ammonium, NH₄⁺, in aqueous systems was studied based on Raman spectroscopic experiments to 600 °C and about 1.3 GPa. Spectra obtained at ambient conditions revealed a strong reduction of the dynamic three-dimensional network of water with addition of ammonium chloride, particularly at small solute concentrations. The differential scattering cross section of the ν₁-NH₄⁺ Raman band in these solutions was found to be similar to that of salammoniac.The Raman band of silica monomers at ∼780 cm⁻¹ was present in all spectra of the fluid at high temperatures in hydrothermal diamond-anvil cell experiments with H₂O ± NH₄Cl and quartz or the assemblage quartz + kyanite + K-feldspar ± muscovite/tobelite. However, these spectra indicated that dissolved silica is less polymerized in ammonium chloride solutions than in comparable experiments with water. Quantification based on the normalized integrated intensity of the H₄SiO₄⁰ band showed that the silica solubility in experiments with H₂O + NH₄Cl was significantly lower than that in equimolal NaCl solutions. This suggests that ammonium causes a stronger decrease in the activity of water in chloridic solutions than sodium.The Raman spectra of the fluid also showed that a significant fraction of ammonium was converted to ammonia, NH₃, in all experiments at temperatures above 300 °C. This indicates a shift towards acidic conditions for experiments without a buffering mineral assemblage. The estimated pH of the fluid was ∼2 at 600 °C, 0.26 GPa, 6.6 m initial NH₄Cl, based on the ratio of the integrated ν₁-NH₃ and ν₁-NH₄⁺ intensities and the HCl⁰ dissociation constant. The NH₃/NH₄⁺ ratio increased with temperature and decreased with pressure. This implies that more ammonium should be retained in K-bearing minerals coexisting with chloridic fluids upon high-P low-T metamorphism. At 500 °C, 0.73 GPa, ammonium partitions preferentially into the fluid, as constrained from infrared spectroscopy on the muscovite and from mass balance.The conversion of K-feldspar to muscovite proceeded much faster in experiments with NH₄Cl solutions than in comparable experiments with water. This is interpreted as being caused by enhancement of the rate-limiting alumina solubility, suggesting complexation of Al with NH₄. Nucleation and growth of mica at the expense of K-feldspar and NH₄⁺/K⁺ exchange between fluid and K-feldspar occurred simultaneously, but incorporation of NH₄⁺ into K-feldspar was distinctly faster than K-feldspar consumption.     

Controls on gold solubility in arc magmas: An experimental study at 1000 °C and 4 kbar

In order to (1) explain the worldwide association between epithermal gold-copper-molybdenum deposits and arc magmas and (2) test the hypothesis that adakitic magmas would be Au-specialized, we have determined the solubility of Au at 4 kbar and 1000 °C for three intermediate magmas (two adakites and one calc-alkaline composition) from the Philippines. The experiments were performed over a fO₂ range corresponding to reducing (∼NNO−1), moderately oxidizing (∼NNO+1.5) and strongly oxidizing (∼NNO+3) conditions as measured by solid Ni-Pd-O sensors. They were carried out in gold containers, the latter serving also as the source of gold, in presence of variable amounts of H₂O and, in a few additional experiments, of S. Concentrations of Au in glasses were determined by LA-ICPMS. Gold solubility in melt is very low (30-240 ppb) but increases with fO₂ in a way consistent with the dissolution of gold as both Au¹⁺ and Au³⁺ species. In the S-bearing experiments performed at ∼NNO−1, gold solubility reaches much higher values, from ∼1200 to 4300 ppb, and seems to correlate with melt S content. No systematic difference in gold solubility is observed between the adakitic and the non-adakitic compositions investigated. Oxygen fugacity and the sulfur concentration in melt are the main parameters controlling the incorporation and concentration of gold in magmas. Certain adakitic and non-adakitic magmas have high fO₂ and magmatic S concentrations favorable to the incorporation and transport of gold. Therefore, the cause of a particular association between some arc magmas and Au-Cu-Mo deposits needs to be searched in the origin of those specialized magmas by involvement of Au- and S-rich protoliths. The subducted slab, which contains metal-rich massive sulfides, may constitute a potentially favorable protolith for the genesis of magmas specialized with respect to gold.     

CO2 solubility in Martian basalts and Martian atmospheric evolution

To understand possible volcanogenic fluxes of CO₂ to the Martian atmosphere, we investigated experimentally carbonate solubility in a synthetic melt based on the Adirondack-class Humphrey basalt at 1-2.5 GPa and 1400-1625 °C. Starting materials included both oxidized and reduced compositions, allowing a test of the effect of iron oxidation state on CO₂ solubility. CO₂ contents in experimental glasses were determined using Fourier transform infrared spectroscopy (FTIR) and Fe³⁺/Fe^T was measured by Mössbauer spectroscopy. The CO₂ contents of glasses show no dependence on Fe³⁺/Fe^T and range from 0.34 to 2.12 wt.%. For Humphrey basalt, analysis of glasses with gravimetrically-determined CO₂ contents allowed calibration of an integrated molar absorptivity of 81,500 ± 1500 L mol⁻¹ cm⁻² for the integrated area under the carbonate doublet at 1430 and 1520 cm⁻¹. The experimentally determined CO₂ solubilities allow calibration of the thermodynamic parameters governing dissolution of CO₂ vapor as carbonate in silicate melt, K_II, (Stolper and Holloway, 1988) as follows: , ΔV⁰ = 20.85 ± 0.91 cm³ mol⁻¹, and ΔH⁰ = −17.96 ± 10.2 kJ mol⁻¹. This relation, combined with the known thermodynamics of graphite oxidation, facilitates calculation of the CO₂ dissolved in magmas derived from graphite-saturated Martian basalt source regions as a function of P, T, and f_O2. For the source region for Humphrey, constrained by phase equilibria to be near 1350 °C and 1.2 GPa, the resulting CO₂ contents are 51 ppm at the iron-wüstite buffer (IW), and 510 ppm at one order of magnitude above IW (IW + 1). However, solubilities are expected to be greater for depolymerized partial melts similar to primitive shergottite Yamato 980459 (Y 980459). This, combined with hotter source temperatures (1540 °C and 1.2 GPa) could allow hot plume-like magmas similar to Y 980459 to dissolve 240 ppm CO₂ at IW and 0.24 wt.% of CO₂ at IW + 1. For expected magmatic fluxes over the last 4.5 Ga of Martian history, magmas similar to Humphrey would only produce 0.03 and 0.26 bars from sources at IW and IW + 1, respectively. On the other hand, more primitive magmas like Y 980459 could plausibly produce 0.12 and 1.2 bars at IW and IW + 1, respectively. Thus, if typical Martian volcanic activity was reduced and the melting conditions cool, then degassing of CO₂ to the atmosphere may not be sufficient to create greenhouse conditions required by observations of liquid surface water. However, if a significant fraction of Martian magmas derive from hot and primitive sources, as may have been true during the formation of Tharsis in the late Noachian, that are also slightly oxidized (IW + 1.2), then significant contribution of volcanogenic CO₂ to an early Martian greenhouse is plausible.     

Constraints on core formation from Pt partitioning in mafic silicate liquids at high temperatures

A new high temperature piston cylinder design has enabled the measurement of platinum solubility in mafic melts at temperatures up to 2500 °C, 2.2 GPa pressure, and under reducing conditions for 1-10 h. These high temperature and low f_O2 conditions may mimic a magma ocean during planetary core formation. Under these conditions, we measured tens to hundreds of ppm Pt in the quenched silicate glass corresponding to , 4-12 orders of magnitude lower than extrapolations from high f_O2 experiments at 1 bar and at temperatures no higher than 1550 °C. Moreover, the new experiments provide coupled textural and compositional evidence that noble metal micro-nuggets, ubiquitous in experimental studies of the highly siderophile elements, can be produced on quench: we measure equally high Pt concentrations in the rapidly quenched nugget-free peripheral margin of the silicate as we do in the more slowly quenched nugget-bearing interior region. We find that both temperature and melt composition exercise strong control on and that Pt⁰ and Pt¹⁺ may contribute significantly to the total dissolved Pt such that low f_O2 does not imply low Pt solubility. Equilibration of metal alloy with liquid silicate in a hot primitive magma might not have depleted platinum to the extent previously believed.     

Oxygen fugacity dependence of Os solubility in haplobasaltic melt

Os equilibrium solubilities were determined at 1350 °C over a wide range of oxygen fugacities (−12 < log fO₂ < −7) applying the mechanically assisted equilibration technique (MAE) at 10⁵ Pa (= 1 bar). Os concentrations in the glass samples were analysed using ID-NTIMS. Additional LA-ICP-MS and SEM analyses were performed to detect, visualize and analyse the nature and chemistry of “nanonuggets.” Os solubilities determined range at a constant temperature of 1350 °C from 0.63 ± 0.04 to 37.4 ± 1.16 ppb depending on oxygen fugacity. At the highest oxygen fugacities, Os³⁺ can be confirmed as the main oxidation state of Os. At low oxygen fugacities (below log fO₂ = −8), samples are contaminated by nanonuggets which, despite the MAE technique, were still not removed entirely from the melt. However, the present results indicate that applying MAE technology does reduce the amount of nanonuggets present significantly, resulting in the lowest Os solubility results reported to date under these experimental conditions, and extending the experimentally accessible range of fO₂ for these studies to lower values. Calculated metal/silicate melt partition coefficients are therefore higher compared to previous studies, making Os more siderophile. Neglecting the as yet unknown temperature dependence of the Os metal/silicate melt partition coefficient, extrapolation of the obtained Os solubilities to conditions for core-mantle equilibrium, results in a , while metallic alloy/silicate melt partition coefficients range from 1.4 × 10⁶ to 8.6 × 10⁷, in agreement with earlier findings. Therefore remains too high by 2-4 orders of magnitude to explain the Os abundance in the Earth’s mantle as result of core-mantle equilibrium during core formation.     

Characterization of elemental release during microbe-basalt interactions at T = 28 °C

This study used batch reactors to characterize the rates and mechanisms of elemental release during the interaction of a single bacterial species (Burkholderia fungorum) with Columbia River Flood Basalt at T = 28 °C for 36 days. We primarily examined the release of Ca, Mg, P, Si, and Sr under a variety of biotic and abiotic conditions with the aim of evaluating how actively metabolizing bacteria might influence basalt weathering on the continents. Four days after inoculating P-limited reactors (those lacking P in the growth medium), the concentration of viable planktonic cells increased from ∼10⁴ to 10⁸ CFU (Colony Forming Units)/mL, pH decreased from ∼7 to 4, and glucose decreased from ∼1200 to 0 μmol/L. Mass-balance and acid-base equilibria calculations suggest that the lowered pH resulted from either respired CO₂, organic acids released during biomass synthesis, or H⁺ extrusion during uptake. Between days 4 and 36, cell numbers remained constant at ∼10⁸ CFU/mL and pH increased to ∼5. Purely abiotic control reactors as well as control reactors containing inert cells (∼10⁸ CFU/mL) showed constant glucose concentrations, thus confirming the absence of biological activity in these experiments. The pH of all control reactors remained near-neutral, except for one experiment where the pH was initially adjusted to 4 but rapidly rose to 7 within 2 days. Over the entire 36 day period, P-limited reactors containing viable bacteria yielded the highest Ca, Mg, Si, and Sr release rates. Release rates inversely correlate with pH, indicating that proton-promoted dissolution was the dominant reaction mechanism. Both biotic and abiotic P-limited reactors displayed low P concentrations. Chemical analyses of bacteria collected at the end of the experiments, combined with mass-balances between the biological and fluid phases, demonstrate that the absence of dissolved P in the biotic reactors resulted from microbial P uptake. The only P source in the basalt is a small amount of apatite (∼1.2%), which occurs as needles within feldspar grains and glass. We therefore conclude that B. fungorum utilized apatite as a P source for biomass synthesis, which stimulated elemental release from coexisting mineral phases via pH lowering. The results of this study suggest that actively metabolizing bacteria have the potential to influence elemental release from basalt in continental settings.     

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.     

Pt-Re-Os and Sm-Nd isotope and HSE and REE systematics of the 2.7 Ga Belingwe and Abitibi komatiites

High-precision Pt-Re-Os and Sm-Nd isotope and highly siderophile element (HSE) and rare earth element (REE) abundance data are reported for two 2.7 b.y. old komatiite lava flows, Tony’s flow (TN) from the Belingwe greenstone belt, Zimbabwe, and the PH-II flow (PH) from Munro Township in the Abitibi greenstone belt, Canada. The emplaced lavas are calculated to have contained ∼25% (TN) and ∼28% (PH) MgO. These lavas were derived from mantle sources characterized by strong depletions in highly incompatible lithophile trace elements, such as light REE (Ce/Sm_N = 0.64 ± 0.02 (TN) and 0.52 ± 0.01 (PH), ε¹⁴³Nd(T) = +2.9 ± 0.2 in both sources). ¹⁹⁰Pt-¹⁸⁶Os and ¹⁸⁷Re-¹⁸⁷Os isochrons generated for each flow yield ages consistent with respective emplacement ages obtained using other chronometers. The calculated precise initial ¹⁸⁶Os/¹⁸⁸Os = 0.1198318 ± 3 (TN) and 0.1198316 ± 5 (PH) and ¹⁸⁷Os/¹⁸⁸Os = 0.10875 ± 17 (TN) and 0.10873 ± 15 (PH) require time-integrated ¹⁹⁰Pt/¹⁸⁸Os and ¹⁸⁷Re/¹⁸⁸Os of 0.00178 ± 11 and 0.407 ± 8 (TN) and 0.00174 ± 18 and 0.415 ± 5 (PH). These parameters, which by far represent the most precise and accurate estimates of time-integrated Pt/Os and Re/Os of the Archean mantle, are best matched by those of enstatite chondrites. The data also provide evidence for a remarkable similarity in the composition of the sources of these komatiites with respect to both REE and HSE. The calculated absolute HSE abundances in the TN and PH komatiite sources are within or slightly below the range of estimates for the terrestrial Primitive Upper Mantle (PUM). Assuming a chondritic composition of the bulk silicate Earth, the strong depletions in LREE, yet chondritic Re/Os in the komatiite sources are apparently problematic because early Earth processes capable of fractionating the LREE might also be expected to fractionate Re/Os. This apparent discrepancy could be reconciled via a two-stage model, whereby the moderate LREE depletion in the sources of the komatiites initially occurred within the first 100 Ma of Earth’s history as a result of either global magma ocean differentiation or extraction and subsequent long-term isolation of early crust, whereas HSE were largely added subsequently via late accretion. The komatiite formation, preceded by derivation of basaltic magmas, was a result of second-stage, large-degree dynamic melting in mantle plumes.     

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.     

Structure of H2O-saturated peralkaline aluminosilicate melt and coexisting aluminosilicate-saturated aqueous fluid determined in-situ to 800 °C and ∼800 MPa

The structure of H₂O-saturated silicate melts and of silicate-saturated aqueous solutions, as well as that of supercritical silicate-rich aqueous liquids, has been characterized in-situ while the sample was at high temperature (to 800 °C) and pressure (up to 796 MPa). Structural information was obtained with confocal microRaman and with FTIR spectroscopy. Two Al-bearing glasses compositionally along the join Na₂O•4SiO₂-Na₂O•4(NaAl)O₂-H₂O (5 and 10 mol% Al₂O₃, denoted NA5 and NA10) were used as starting materials. Fluids and melts were examined along pressure-temperature trajectories of isochores of H₂O at nominal densities (from PVT properties of pure H₂O) of 0.85 g/cm³ (NA10 experiments) and 0.86 g/cm³ (NA5 experiments) with the aluminosilicate + H₂O sample contained in an externally-heated, Ir-gasketed hydrothermal diamond anvil cell.Molecular H₂O (H₂O°) and OH groups that form bonds with cations exist in all three phases. The OH/H₂O° ratio is positively correlated with temperature and pressure (and, therefore, fugacity of H₂O, f_H2O) with (OH/H₂O°)^melt > (OH/H₂O°)^fluid at all pressures and temperatures. Structural units of Q³, Q², Q¹, and Q⁰ type occur together in fluids, in melts, and, when outside the two-phase melt + fluid boundary, in single-phase liquids. The abundance of Q⁰ and Q¹ increases and Q² and Q³ decrease with f_H2O. Therefore, the NBO/T (nonbridging oxygen per tetrahedrally coordination cations), of melt is a positive function of f_H2O. The NBO/T of silicate in coexisting aqueous fluid, although greater than in melt, is less sensitive to f_H2O.The melt structural data are used to describe relationships between activity of H₂O and melting phase relations of silicate systems at high pressure and temperature. The data were also combined with available partial molar configurational heat capacity of Qⁿ-species in melts to illustrate how these quantities can be employed to estimate relationships between heat capacity of melts and their H₂O content.     

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

Microbial dissolution of calcite at T = 28 °C and ambient pCO2

This study used batch reactors to quantify the mechanisms and rates of calcite dissolution in the presence and absence of a single heterotrophic bacterial species (Burkholderia fungorum). Experiments were conducted at T = 28°C and ambient pCO₂ over time periods spanning either 21 or 35 days. Bacteria were supplied with minimal growth media containing either glucose or lactate as a C source, NH₄⁺ as an N source, and H₂PO₄⁻ as a P source. Combining stoichiometric equations for microbial growth with an equilibrium mass-balance model of the H₂O-CO₂-CaCO₃ system demonstrates that B. fungorum affected calcite dissolution by modifying pH and alkalinity during utilization of ionic N and C species. Uptake of NH₄⁺ decreased pH and alkalinity, whereas utilization of lactate, a negatively charged organic anion, increased pH and alkalinity. Calcite in biotic glucose-bearing reactors dissolved by simultaneous reaction with H₂CO₃ generated by dissolution of atmospheric CO₂ (H₂CO₃ + CaCO₃ → Ca²⁺ + 2HCO₃⁻) and H⁺ released during NH₄⁺ uptake (H⁺ + CaCO₃ → Ca²⁺ + HCO₃⁻). Reaction with H₂CO₃ and H⁺ supplied ∼45% and 55% of the total Ca²⁺ and ∼60% and 40% of the total HCO₃⁻, respectively. The net rate of microbial calcite dissolution in the presence of glucose and NH₄⁺ was ∼2-fold higher than that observed for abiotic control experiments where calcite dissolved only by reaction with H₂CO₃. In lactate bearing reactors, most H⁺ generated by NH₄⁺ uptake reacted with HCO₃⁻ produced by lactate oxidation to yield CO₂ and H₂O. Hence, calcite in biotic lactate-bearing reactors dissolved by reaction with H₂CO₃ at a net rate equivalent to that calculated for abiotic control experiments. This study suggests that conventional carbonate equilibria models can satisfactorily predict the bulk fluid chemistry resulting from microbe-calcite interactions, provided that the ionic forms and extent of utilization of N and C sources can be constrained. Because the solubility and dissolution rate of calcite inversely correlate with pH, heterotrophic microbial growth in the presence of nonionic organic matter and NH₄⁺ appears to have the greatest potential for enhancing calcite weathering relative to abiotic conditions.     

The influence of sulfur on partitioning of siderophile elements

Sulfur is a potential light element in the liquid outer core of the Earth. Its presence in segregating metal may have had an influence in distribution of metal-loving (siderophile) elements during early accretion and core formation events in the Earth. The observed “excess” abundance of siderophile elements in the terrestrial mantle, relative to an abundance expected from simple core-mantle equilibrium at low temperature and pressure, may indicate a reduction in the iron-loving tendency of siderophile elements in the presence of sulfur in the metallic phase. The present experimental partitioning study between iron-carbon-sulfur-siderophile element bearing liquid metal and liquid silicate shows that for some siderophile elements this sulfur effect may be significant enough to even change their character to lithophile. Large and intricate variations in metal-silicate partition coefficients (D^met/sil) have been observed for many elements, e.g., Ni, Co, Ge, W, P, Au, and Re as a function of sulfur content. Moderately siderophile elements Ge, P, and W show the most significant response (sulfur-avoidance) by an enhanced segregation into the associated sulfur-deficient phases. Highly siderophile elements Ir, Pt, and Re show a different style of sulfur-avoidance (alloy-preference) by segregating as sulfur-poor, siderophile element-rich alloys. Both groups are chalcophobic. D^met/sil for Ni, Co, and Au moderately decreases with increasing sulfur-content in the liquid metal. D^met/sil for chalcophile element, Cr, in contrast, increases with sulfur. Irrespective of the sulfur-content, in the presence of a carbon-saturated liquid metal, P is always lithophile. The general nonmetal-avoidance tendency of siderophile elements (and acceptance of chalcophile elements) in the liquid metal, postulated by Jones and Malvin (1990) in the FeNiS(sulfur)M (siderophile) system is found to be present in the metal-silicate system as well. A sulfur-bearning liquid metal segregation can potentially reduce the metal-loving nature of many elements to explain the excess paradox. Sulfur-bearing core segregation, however, might require an efficient draining of exsolved immiscible sulfide liquids from the molten silicate, or an increasing siderophility of sulfur at high pressure to reduce the mantle sulfur content to the observed (<300 ppm) value. Moreover, the chondritic relative abundance pattern of many moderately or highly siderophile elements in the upper mantle is not explained by the presence of sulfur in the segregating metals. Core formation is more complex and intricate than equilibrium segregation.     

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.     

Radium uptake during barite recrystallization at 23 ± 2 °C as a function of solution composition: An experimental Ba and Ra tracer study

High-purity synthetic barite powder was added to pure water or aqueous solutions of soluble salts (BaCl₂, Na₂SO₄, NaCl and NaHCO₃) at 23 ± 2 °C and atmospheric pressure. After a short pre-equilibration time (4 h) the suspensions were spiked either with ¹³³Ba or ²²⁶Ra and reacted under constant agitation during 120-406 days. The pH values ranged from 4 to 8 and solid to liquid (S/L) ratios varied from 0.01 to 5 g/l. The uptake of the radiotracers by barite was monitored through repeated sampling of the aqueous solutions and radiometric analysis. For both ¹³³Ba and ²²⁶Ra, our data consistently showed a continuous, slow decrease of radioactivity in the aqueous phase.Mass balance calculations indicated that the removal of ¹³³Ba activity from aqueous solution cannot be explained by surface adsorption only, as it largely exceeded the 100% monolayer coverage limit. This result was a strong argument in favor of recrystallization (driven by a dissolution-precipitation mechanism) as the main uptake mechanism. Because complete isotopic equilibration between aqueous solution and barite was approached or even reached in some experiments, we concluded that during the reaction all or substantial fractions of the initial solid had been replaced by newly formed barite.The ¹³³Ba data could be successfully fitted assuming constant recrystallization rates and homogeneous distribution of the tracer into the newly formed barite. An alternative model based on partial equilibrium of ¹³³Ba with the mineral surface (without internal isotopic equilibration of the solid) could not reproduce the measured activity data, unless multistage recrystallization kinetics was assumed. Calculated recrystallization rates in the salt solutions ranged from 2.8 × 10⁻¹¹ to 1.9 × 10⁻¹⁰ mol m⁻² s⁻¹ (2.4-16 μmol m⁻² d⁻¹), with no specific trend related to solution composition. For the suspensions prepared in pure water, significantly higher rates (∼5.7 × 10⁻¹⁰ mol m⁻² s⁻¹ or ∼49 μmol m⁻² d⁻¹) were determined.Radium uptake by barite was determined by monitoring the decrease of ²²⁶Ra activity in the aqueous solution with alpha spectrometry, after filtration of the suspensions and sintering. The evaluation of the Ra uptake experiments, in conjunction with the recrystallization data, consistently indicated formation of non-ideal solid solutions, with moderately high Margules parameters (W_AB = 3720-6200 J/mol, a₀ = 1.5-2.5). These parameters are significantly larger than an estimated value from the literature (W_AB = 1240 J/mol, a₀ = 0.5).In conclusion, our results confirm that radium forms solid solutions with barite at fast kinetic rates and in complete thermodynamic equilibrium with the aqueous solutions. Moreover, this study provides quantitative thermodynamic data that can be used for the calculation of radium concentration limits in environmentally relevant systems, such as radioactive waste repositories and uranium mill tailings.