Platinum solubility and partitioning in a felsic melt-vapor-brine assemblage

The solubility and partitioning of Pt in a S-free vapor - brine - rhyolite melt - Pt metal assemblage has been quantified at 800 °C, f_O2=NNO and pressures of 100 and 140 MPa. Vapor and brine were sampled at run conditions by trapping these phases as glass-hosted fluid inclusions as the melt cooled through the glass transition temperature. The vapor and brine were in equilibrium with the melt at the time of trapping and, thus, represent fluids which were sampled at the termination of each experimental run. The microthermometrically determined salinities of vapor and brine are ∼2 and ∼63 wt.% NaCl eq. and ∼9 and ∼43 wt.% NaCl eq. at 100 and 140 MPa, respectively. Platinum solubilities in vapor, brine and glass (i.e., quenched melt) were quantified by using laser ablation - inductively coupled plasma - mass spectrometry (LA-ICP-MS). Equilibrium is discussed with reference to the major and trace element concentrations of glass-hosted fluid inclusions as well as the silicate melt over run times that varied from 110 to 377 h at 140 MPa and 159 to 564 h at 100 MPa. Platinum solubility values (±1σ) in H₂O-saturated felsic melt are 0.28 ± 0.13 μg/g and 0.38 ± 0.06 μg/g at 140 and 100 MPa, respectively. Platinum solubility values () at 140 and 100 MPa, respectively, in aqueous vapor are 0.91 ± 0.29 μg/g and 0.37 ± 0.17 μg/g and in are brine 16 ± 10 μg/g and 3.3 ± 1.0 μg/g. The measured solubility data were used to calculate Nernst-type partition coefficients for Pt between vapor/melt, brine/melt and vapor/brine. The partition coefficient values () for vapor/melt, brine/melt and vapor/brine at 140 MPa are 2.9 ± 1.0, 67 ± 27, and 0.13 ± 0.05 and at 100 MPa are 1.0 ± 0.2, 6.8 ± 2.4, and 0.15 ± 0.05. The partitioning data were used to model the Pt-scavenging capacity of vapor and brine during the crystallization-driven degassing (i.e., second boiling) of a felsic silicate melt over a depth range (i.e., 3-6 km) consistent with the evolution of magmatic-hydrothermal ore deposits. Model calculations suggest that aqueous vapor and brine can scavenge sufficient quantities of Pt, and by analogy other platinum group elements (PGE), to produce economically important PGE-rich magmatic-hydrothermal ore deposits in Earth’s upper continental crust.     

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.     

Copper partitioning in a melt-vapor-brine-magnetite-pyrrhotite assemblage

The effect of sulfur on the partitioning of Cu in a melt-vapor-brine ± magnetite ± pyrrhotite assemblage has been quantified at 800 °C, 140 MPa, f_O2 = nickel-nickel oxide (NNO), logf_S2=-3.0 (i.e., on the magnetite-pyrrhotite curve at NNO), logf_H2S=-1.3 and logf_SO2=-1. All experiments were vapor + brine saturated. Vapor and brine fluid inclusions were trapped in silicate glass and self-healed quartz fractures. Vapor and brine are dominated by NaCl, KCl and HCl in the S-free runs and NaCl, KCl and FeCl₂ in S-bearing runs. Pyrrhotite served as the source of sulfur in S-bearing experiments. The composition of fluid inclusions, glass and crystals were quantified by laser-ablation inductively coupled plasma mass spectrometry. Major element, chlorine and sulfur concentrations in glass were quantified by using electron probe microanalysis. Calculated Nernst-type partition coefficients (±2σ) for Cu between melt-vapor, melt-brine and vapor-brine are , , and , respectively, in the S-free system. The partition coefficients (±2σ) for Cu between melt-vapor, melt-brine and vapor-brine are , , and , respectively, in the S-bearing system. Apparent equilibrium constants (±1σ) describing Cu and Na exchange between vapor and melt and brine and melt were also calculated. The values of are 34 ± 21 and 128 ± 29 in the S-free and S-bearing runs, respectively. The values of are 33 ± 22 and60 ± 5 in the S-free and S-bearing runs, respectively. The data presented here indicate that the presence of sulfur increases the mass transfer of Cu into vapor from silicate melt. Further, the nearly threefold increase in suggests that Cu may be transported as both a chloride and sulfide complex in magmatic vapor, in agreement with hypotheses based on data from natural systems. Most significantly, the data demonstrate that the presence of sulfur enhances the partitioning of Cu from melt into magmatic volatile phases.     

The partitioning behavior of As and Au in S-free and S-bearing magmatic assemblages

The partitioning of As and Au between rhyolite melt and low-salinity vapor (2 wt% NaCl eq.) in a melt-vapor-Au metal ± magnetite ± pyrrhotite assemblage has been quantified at 800 °C, 120 MPa and f_O2=NNO. The S-bearing runs have calculated values for the fugacities of H₂S, SO₂ and S₂ of logf_H2S=1.1, logf_SO2=-1.5, and logf_S2=-3.0. The ratio of H₂S to SO₂ is on the order of 400. The experiments constrain the effect of S on the partitioning behavior of As and Au at magmatic conditions. Calculated average Nernst-type partition coefficients (±1σ) for As between vapor and melt, , are 1.0 ± 0.1 and 2.5 ± 0.3 in the S-free and S-bearing assemblages, respectively. These results suggest that sulfur has a small, but statistically meaningful, effect on the mass transfer of As between silicate melt and low-salinity vapor at the experimental conditions. Efficiencies of removal, calculated following Candela and Holland (1986), suggest that the S-free and S-bearing low-salinity vapor can scavenge approximately 41% and 63% As from water-saturated rhyolite melt, respectively, during devolatilization assuming that As is partitioned into magnetite and pyrrhotite during second boiling. The S-free data are consistent with the presence of arsenous acid, As(OH)₃ in the vapor phase. However, the S-bearing data suggest the presence of both arsenous acid and a As-S complex in S-bearing magmatic vapor. Apparent equilibrium constants, , describing the partitioning of As between melt and vapor are −1.3 (0.1) and −1.1 (0.1) for the S-free and S-bearing runs, respectively. The increase in the value of with the addition of S suggests a role for S in complexing and scavenging As from the melt during degassing.The calculated vapor/melt partition coefficients (±1σ) for Au between vapor and melt, , in S-free and S-bearing assemblages are 15 ± 2.5 and 12 ± 0.3, respectively. Efficiencies of removal (Candela and Holland, 1986) for the S-free melt, calculated assuming that magnetite is the dominant Au-sequestering solid phase during crystallization (Simon et al., 2003), suggest that magmatic vapor may scavenge on the order of 72% Au from a water-saturated melt. Efficiencies of removal calculated for the S-bearing assemblage, assuming pyrrhotite and magnetite are the dominant Au-sequestering solid phases, indicate that vapor may scavenge on the order of 60% Au from the melt. These model calculations suggest that the loss of pyrrhotite and magnetite from a melt, owing to punctuated differentiation during ascent and emplacement, does not prohibit the ability of a rhyolite melt to generate a large-tonnage Au deposit. Apparent equilibrium constants describing the partitioning of Au between melt and vapor were calculated using the mean values for the S-free and S-bearing assemblages; only S-bearing data from runs longer than 400 h were used as shorter runs may not have reached equilibrium with respect only to vapor/melt partitioning of Au. The values for are −4.4 (0.1) and −4.2 (0.2) for the S-free and S-bearing runs, respectively. These data suggest that the presence of S does not affect the mass transfer of Au from degassing silicate melt to an exsolved, low-salinity vapor in a low-f_S2 assemblage (i.e., pyrrhotite-magnetite at NNO) at the experimental conditions reported here. Efficiencies of removal are calculated and used to model the mass transfer of Au from a crystallizing silicate melt to an exsolved, low-salinity vapor phase. The calculations suggest that the model, absolute tonnage of Au scavenged and transported by S-free and S-bearing vapors, from a crystallizing melt, would be comparable and that the time-integrated flux of low-salinity vapor could be responsible for a significant quantity of the Au in magmatic-hydrothermal ore deposits.     

Experimental constraints on the mobility of Rhenium in silicate liquids

The volatization of Rhenium (Re) from melts of natural basalt, dacite and a synthetic composition in the CaO-MgO-Al₂O₃-SiO₂ system has been investigated at 0.1 MPa and 1250-1350 °C over a range of fO₂ conditions from log fO₂ = −10 to −0.68. Experiments were conducted using open top Pt crucibles doped with Re and Yb. Analysis of quenched glasses by laser ablation-inductively coupled plasma mass spectrometry (LA-ICP-MS) normal to the melt/gas interface showed concentration profiles for Re, to which a semi-infinite one-dimensional diffusion model could be applied to extract diffusion coefficients (D). The results show Re diffusivity in basalt at 1300 °C in air is log D_Re = −7.2 ± 0.3 cm²/s and increases to log D_Re = −6.6 ± 0.3 cm²/s when trace amounts of Cl were added to the starting material. At fO₂ conditions below the nickel-nickel oxide (NNO) buffer Re diffusivity decreases to and to in dacitic melt. In the CMAS composition, . The diffusivity of Re is comparable to Ar and CO₂ in basalt at 500 MPa favoring its release as a volatile. Our results support the contention that subaerial degassing is the cause of lower Re concentrations in arc-type and ocean island basalts compared to mid-ocean ridge basalts.     

Trace element partitioning between ilmenite, armalcolite and anhydrous silicate melt: Implications for the formation of lunar high-Ti mare basalts

We performed a series of experiments at high pressures and temperatures to determine the partitioning of a wide range of trace elements between ilmenite (Ilm), armalcolite (Arm) and anhydrous lunar silicate melt, to constrain geochemical models of the formation of titanium-rich melts in the Moon. Experiments were performed in graphite-lined platinum capsules at pressures and temperatures ranging from 1.1 to 2.3 GPa and 1300-1400 °C using a synthetic Ti-enriched Apollo ‘black glass’ composition in the CaO-FeO-MgO-Al₂O₃-TiO₂-SiO₂ system. Ilmenite-melt and armalcolite-melt partition coefficients (D) show highly incompatible values for the rare earth elements (REE) with the light REE more incompatible compared to the heavy REE ( 0.0020 ± 0.0010 to 0.069 ± 0.010 for ilmenite; 0.0048 ± 0.0023 to 0.041 ± 0.008 for armalcolite). D values for the high field strength elements vary from highly incompatible for Th, U and to a lesser extent W (for ilmenite: 0.0013 ± 0.0008, 0.0035 ± 0.0015 and 0.039 ± 0.005, and for armalcolite 0.008 ± 0.003, 0.0048 ± 0.0022 and 0.062 ± 0.03), to mildly incompatible for Nb, Ta, Zr, and Hf (e.g. 0.28 ± 0.05 and : 0.76 ± 0.07). Both minerals fractionate the high field strength elements with D_Ta/D_Nb and D_Hf/D_Zr between 1.3 and 1.6 for ilmenite and 1.3 and 1.4 for armalcolite. Armalcolite is slightly more efficient at fractionating Hf from W during lunar magma ocean crystallisation, with D_Hf/D_W = 12-13 compared to 6.7-7.5 for ilmenite. The transition metals vary from mildly incompatible to compatible, with the highest compatibilities for Cr in ilmenite (D ∼ 7.5) and V in armalcolite (D ∼ 8.1). D values show no clear variation with pressure in the small range covered.Crystal lattice strain modelling of D values for di-, tri- and tetravalent trace elements shows that in ilmenite, divalent elements prefer to substitute for Fe while armalcolite data suggest REE replacing Mg. Tetravalent cations appear to preferentially substitute for Ti in both minerals, with the exception of Th and U that likely substitute for the larger Fe or Mg cations. Crystal lattice strain modelling is also used to identify and correct for very small (∼0.3 wt.%) melt contamination of trace element concentration determinations in crystals.Our results are used to model the Lu-Hf-Ti concentrations of lunar high-Ti mare basalts. The combination of their subchondritic Lu/Hf ratios and high TiO₂ contents requires preferential dissolution of ilmenite or armalcolite from late-stage, lunar magma ocean cumulates into low-Ti partial melts of deeper pyroxene-rich cumulates.     

Thermodynamic properties of H4SiO4 in the ideal gas state as evaluated from experimental data

Solid phases of silicon dioxide react with water vapor with the formation of hydroxides and oxyhydroxides of silica. Recent transpiration and mass-spectrometric studies convincingly demonstrate that H₄SiO₄ is the predominant form of silica in vapor phase at water pressure in excess of 10⁻² MPa. Available literature transpiration and solubility data for the reactions of solid SiO₂ phases and low-density water, extending from 424 to 1661 K, are employed for the determination of Δ_fG⁰, Δ_fH⁰ and S⁰ of H₄SiO₄ in the ideal gas state at 298.15 K, 0.1 MPa. In total, there are 102 data points from seven literature sources. The resulting values of the thermodynamic functions of H₄SiO₄(g) are: Δ_fG⁰ = −1238.51 ± 3.0 kJ mol⁻¹, Δ_fH⁰ = −1340.68 ± 3.5 kJ mol⁻¹ and S⁰ = 347.78 ± 6.2 J K⁻¹ mol⁻¹. These values agree quantitatively with one set of ab initio calculations. The relatively large uncertainties are mainly due to conflicting data for H₄SiO₄(g) from various sources, and new determinations of would be helpful. The thermodynamic properties of this species, H₄SiO₄(g), are necessary for realistic modeling of silica transport in a low-density water phase. Applications of this analysis may include the processes of silicates condensation in the primordial solar nebula, the precipitation of silica in steam-rich geothermal systems and the corrosion of SiO₂-containing alloys and ceramics in moist environments.     

Experimental studies of metal-silicate partitioning of Sb: Implications for the terrestrial and lunar mantles

The terrestrial mantle has a well defined Sb depletion of ∼7 ± 1 (Jochum and Hofmann, 1997), and the lunar mantle is depleted relative to the Earth by a factor of ∼50 ± 5 (Wolf and Anders, 1980). Despite these well defined depletions, there are few data upon which to evaluate their origin—whether due to volatility or core formation. We have carried out a series of experiments to isolate several variables such as oxygen fugacity, temperature, pressure, and silicate and metallic melt compositions, on the magnitude of . The activity of Sb in FeNi metal is strongly composition dependent such that solubility of Sb as a function of fO₂ must be corrected for the metal composition. When the correction is applied, Sb solubility is consistent with 3+ valence. Temperature series (at 1.5 GPa) shows that decreases by a factor of 100 over 400 °C, and a pressure series exhibits an additional decrease between ambient pressure (100 MPa) and 13 GPa. A strong dependence upon silicate melt composition is evident from a factor of 100 decrease in between nbo/t values of 0.3 and 1.7. Consideration of all these variables indicates that the small Sb depletion for the Earth’s mantle can be explained by high PT equilibrium partitioning between metal and silicate melt . The relatively large lunar Sb depletion can also be explained by segregation of a small metallic core, at lower pressure conditions where is much higher (2500).     

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.     

Solubility of NaCl in CO2 at high pressure and temperature: First experimental measurements

NaCl solubility in gaseous carbon dioxide has been measured in the pressure range from 30 to 70 MPa at 623 and 673 K. Our originally-designed high pressure apparatus allows in situ sampling of a portion of the fluid phase for chemical analysis. The results indicate that the solubility of NaCl increases with both temperature and pressure, and is about 4-5 orders of magnitude higher than saturated NaCl pressure values at the same temperature conditions (6.02 × 10⁻¹² at 623 K and 1.51 × 10⁻¹⁰ at 673 K). It is also 1-2 orders of magnitude greater than predictions according to the Equation of State of the ternary H₂O-CO₂-NaCl system by Duan, Moeller and Weare [Duan, Z., Moller, N., and Weare, J. H. (1995) Equation of state for the NaCl-H₂O-CO₂ system: prediction of phase equilibria and volumetric properties. Geochim. Cosmochim. Acta59, 2869] and has the opposite pressure dependence. The activity values of NaCl in the vapor phase, calculated from the experiments (with pure molten NaCl as a standard state in the vapor), have been fitted to the Darken Quadratic Formalism: , where, x_NaCl,v is mole the fraction of NaCl in the vapor phase, , , where P is the pressure in MPa and T the absolute temperature. Caution should be exerted while extrapolating this empirical equation far beyond the experimental P-T-compositional range.     

Olivine/melt transition metal partitioning, melt composition, and melt structure—Influence of Al for Si substitution in the tetrahedral network of silicate melts

The influence on olivine/melt transition metal (Mn, Co, Ni) partitioning of substitution in the tetrahedral network of silicate melt structure has been examined at ambient pressure in the 1450-1550 °C temperature range. Experiments were conducted in the systems NaAlSiO₄-Mg₂SiO₄- SiO₂ and CaAl₂Si₂O₈-Mg₂SiO₄-SiO₂ with about 1 wt% each of MnO, CoO, and NiO added. These compositions were used to evaluate how, in silicate melts, substitution and ionization potential of charge-balancing cations affect activity-composition relations in silicate melts and mineral/melt partitioning.The exchange equilibrium coefficient, , is a positive and linear function of melt Al/(Al + Si) at constant degree of melt polymerization, NBO/T. The is negatively correlated with the ionic radius, r, of the M-cation and also with the ionization potential (Z/r², Z = electrical charge) of the cation that serves to charge-balance Al³⁺ in tetrahedral coordination in the melts. The activity coefficient ratio, (γ_M/γ_Mg)^melt, is therefore similarly correlated.These melt composition relationships are governed by the distribution of Al³⁺ among coexisting Q-species in the peralkaline (depolymerized) melts coexisting with olivine. This distribution controls Q-speciation abundance, which, in turn, controls (γ_M/γ_Mg)^melt and . The relations between melt structure and olivine/melt partitioning behavior lead to the suggestion that in natural magmatic systems mineral/melt partition coefficients are more dependent on melt composition and, therefore, melt structure the more alkali-rich and the more felsic the melt. Moreover, mineral/melt partition coefficients are more sensitive to melt composition the more highly charged or the smaller the ionic radius of the cation of interest.     

Silicon isotopic fractionation during adsorption of aqueous monosilicic acid onto iron oxide

The quantification of silicon isotopic fractionation by biotic and abiotic processes contributes to the understanding of the Si continental cycle. In soils, light Si isotopes are selectively taken up by plants, and concentrate in secondary clay-sized minerals. Si can readily be retrieved from soil solution through the specific adsorption of monosilicic acid () by iron oxides. Here, we report on the Si-isotopic fractionation during adsorption on synthesized ferrihydrite and goethite in batch experiment series designed as function of time (0-504 h) and initial concentration (ic) of Si in solution (0.21-1.80 mM), at 20 °C, constant pH (5.5) and ionic strength (1 mM). At various contact times, the δ²⁹Si vs. NBS28 compositions were determined in selected solutions (ic = 0.64 and 1.06 mM Si) by MC-ICP-MS in dry plasma mode with external Mg doping with an average precision of ±0.08‰ (±2σ_SEM). Per oxide mass, ferrihydrite (74-86% of initial Si loading) adsorbed more Si than goethite (37-69%) after 504 h of contact over the range of initial Si concentration 0.42-1.80 mM. Measured against its initial composition (δ²⁹Si = +0.01 ± 0.04‰ (±2σ_SD)), the remaining solution was systematically enriched in ²⁹Si, reaching maximum δ²⁹Si values of +0.70 ± 0.07‰ for ferrihydrite and +0.50 ± 0.08‰ for goethite for ic 1.06 mM. The progressive ²⁹Si enrichment of the solution fitted better a Rayleigh distillation path than a steady state model. The fractionation factor ²⁹ε (±1σ_SD) was estimated at −0.54 ± 0.03‰ for ferrihydrite and −0.81 ± 0.12‰ for goethite. Our data imply that the sorption of onto synthetic iron oxides produced a distinct Si-isotopic fractionation for the two types of oxide but in the same order than that generated by Si uptake by plants and diatoms. They further suggest that the concentration of light Si isotopes in the clay fraction of soils is partly due to sorption onto secondary clay-sized iron oxides.     

The effect of sulfur on the partitioning of Ni and other first-row transition elements between olivine and silicate melt

The effect of sulfur dissolved as sulfide (S²⁻) in silicate melts on the activity coefficients of NiO and some other oxides of divalent cations (Ca, Cr, Mn, Fe and Co) has been determined from olivine/melt partitioning experiments at 1400 °C in six melt compositions in the system CaO-MgO-Al₂O₃-SiO₂ (CMAS), and in derivatives of these compositions at 1370 °C, obtained from the six CMAS compositions by substituting Fe for Mg (FeCMAS). Amounts of S²⁻ were varied from zero to sulfide saturation, reaching 4100 μg g⁻¹ S in the most sulfur-rich silicate melt. The sulfide solubilities compare reasonably well with those predicted from the parameterization of the sulfide capacity of silicate melts at 1400 °C of O’Neill and Mavrogenes (2002), although in detail systematic deviations indicate that a more sophisticated model may improve the prediction of sulfide capacities.The results show a barely discernible effect of S²⁻ in the silicate melt on Fe, Co and Ni partition coefficients, and also surprisingly, a tiny but resolvable effect on Ca partitioning, but no detectable effect on Cr, Mn or some other lithophile incompatible elements (Sc, Ti, V, Y, Zr and Hf). Decreasing Mg# of olivine (reflecting increasing FeO in the system) has a significant influence on the partitioning of several of the divalent cations, particularly Ca and Ni. We find a remarkably systematic correlation between and the ionic radius of M²⁺, where M = Ca, Cr, Mn, Fe, Co or Ni, which is attributable to a simple relationship between size mismatch and excess free energies of mixing in Mg-rich olivine solid solutions.Neither the effect of S²⁻ nor of Mg#^ol is large enough by an order of magnitude to account for the reported variations of obtained from electron microprobe analyses of olivine/glass pairs from mid-ocean ridge basalts (MORBs). Comparing these MORB glass analyses with the Ni-MgO systematics of MORB from other studies in the literature, which were obtained using a variety of analytical techniques, shows that these electron microprobe analyses are anomalous. We suggest that the reported variation of with S content in MORB is an analytical artifact.Mass balance of melt and olivine compositions with the starting compositions shows that dissolved S²⁻ depresses the olivine liquidus of haplobasaltic silicate melts by 5.8 × 10⁻³ (±1.3 × 10⁻³) K per μg g⁻¹ of S²⁻, which is negligible in most contexts. We also present data for the partitioning of some incompatible trace elements (Sc, Ti, Y, Zr and Hf) between olivine and melt. The data for Sc and Y confirm previous results showing that and decrease with increasing SiO₂ content of the melt. Values of average 0.01 with most falling in the range 0.005-0.015. Zr and Hf are considerably more incompatible than Ti in olivine, with and about 10⁻³. The ratio / is well constrained at 0.611 ± 0.016.     

The problem of inherited Ar* in dating impact glass by the Ar/Ar method: Evidence from the Tswaing impact crater (South Africa)

The Tswaing meteorite impact crater is a 1.13 km diameter structure located in the 2.05 Ga Nebo granite of the Bushveld Complex. The impact age had previously been determined by fission track dating to 220 ± 104 ka. ⁴⁰Ar/³⁹Ar step-heating and total fusion experiments performed on single- and multi-grain impact glass aliquots gave apparent ages ranging from 1.0 ± 0.3 Ma to 204 ± 6 Ma. These “ages” indicate that the radiogenic Ar derived from the target rocks has not been completely degassed as a result of the impact process, despite fusion of the target material. Results of step-heating experiments imply that the trapped within the glass is located in two distinct reservoirs thought to be the glass matrix and fluid/vapor inclusions (or un-melted residual clasts). Calculations assuming an age of 0.2 ± 0.1 Ma for Tswaing (fission track data) reveal that the amount of inherited ⁴⁰Ar*() relative to the pre-impact concentration varies from 0.015% to 4.15%. The spread defined by likely reflects the various quench rates experienced by the glass, most certainly due to the pre-impact position of the sample relative to the center of the crater. We compare the influence of on the apparent ⁴⁰Ar/³⁹Ar age determination of five impact structures. Our calculations show that the main characteristic controlling the age offset (for a given proportion of ) is the age difference between the impact and the target rocks (i.e., the ⁴⁰Ar* concentration in the target rock). The buffer effect for a given crater structure can be predicted knowing the age of the basement and having a rough estimation of the age of the crater structure itself. The occurrence of is likely influenced by (1) the degree of polymerization (i.e., silicate structure complexity) of the target rock and presumably related to the diffusivity of Ar in the melt and glass, (2) the Ar partial pressure at the grain boundary, (3) the quantity of energy involved in the impact, and (4) the porosity of the target rocks. For glass that inevitably suffers inherited and/or excess ⁴⁰Ar*, the use of the inverse isochron technique can be appropriate but should be applied with careful statistical treatment.     

Carbon solubility in core melts in a shallow magma ocean environment and distribution of carbon between the Earth’s core and the mantle

The solubility of carbon in Fe and Fe-5.2 wt.% Ni melts, saturated with graphite, determined by electron microprobe analysis of quenched metal melts was 5.8 ± 0.1 wt.% at 2000 °C, 6.7 ± 0.2 wt.% at 2200 °C, and 7.4 ± 0.2 wt.% at 2410 °C at 2 GPa, conditions relevant for core/mantle differentiation in a shallow magma ocean. These solubilities are slightly lower than low-pressure literature values and significantly beneath calculated values for even higher pressures [e.g., Wood B. J. (1993) Carbon in the core. Earth Planet. Sci. Lett.117, 593-607]. The trend of C solubility versus temperature for Fe-5.2 wt.% Ni melt, within analytical uncertainties, is similar to or slightly lower (∼0.2-0.4 wt.%) than that of pure Fe. Carbon content of core melts and residual mantle silicates derived from equilibrium batch or fractional segregation of core liquids and their comparison with our solubility data and carbon content estimate of the present day mantle, respectively, constrain the partition coefficient of carbon between silicate and metallic melts, in a magma ocean. For the entire range of possible bulk Earth carbon content from chondritic to subchondritic values, of 10⁻⁴ to 1 is derived. But for ∼1000 ppm bulk Earth carbon, is between 10⁻² and 1. Using the complete range of possible for a magma ocean at ∼2200 °C, we predict maximum carbon content of the Earth’s core to be ∼6-7 wt.% and a preferred value of 0.25 ± 0.15 wt.% for a bulk Earth carbon concentration of ∼1000 ppm.     

Thermodynamic properties of anhydrous smectite MX-80, illite IMt-2 and mixed-layer illite-smectite ISCz-1 as determined by calorimetric methods. Part I: Heat capacities, heat contents and entropies

The heat capacities of the anhydrous international reference clay minerals, smectite MX-80, illite IMt-2 and mixed-layer illite-smectite ISCz-1, were measured by low temperature adiabatic calorimetry and differential scanning calorimetry, from 6 to 520 K (at 1 bar). The samples were chemically purified and Na-saturated. Dehydrated clay fractions <2 μm were studied. The structural formulae of the corresponding clay minerals, obtained after subtracting the remaining impurities, are K_0.026Na_0.435Ca_0.010(Si_3.612Al_0.388) (Al_1.593Mg_0.228Ti_0.011)O₁₀(OH)₂ for smectite MX-80, K_0.762Na_0.044(Si_3.387Al_0.613) (Al_1.427Mg_0.241O₁₀(OH)₂ for illite IMt-2 and K_0.530Na_0.135(Si_3.565Al_0.435)(Al_1.709Mg_0.218Ti_0.005)O₁₀(OH)₂for mixed-layer ISCz-1. From the heat capacity values, we determined the molar entropies, standard entropies of formation and heat contents of these minerals. The following values were obtained at 298.15 K and 1 bar:

	(J mol−1 K−1)	S0 (J mol−1 K−1)
Smectite MX-80	326.13 ± 0.10	280.56 ± 0.16
Illite IMt-2	328.21 ± 0.10	295.05 ± 0.17
Mixed-layer ISCz-1	320.79 ± 0.10	281.62 ± 0.15

Full-size table

     

A combined Sm-Nd, Rb-Sr, and U-Pb isotopic study of Mg-suite norite 78238: Further evidence for early differentiation of the Moon

Lunar Mg-suite norite 78238 was dated using the Sm-Nd, Rb-Sr, and U-Pb isotopic systems in order to constrain the age of lunar magma ocean solidification and the beginning of Mg-suite magmatism, as well as to provide a direct comparison between the three isotopic systems. The Sm-Nd isotopic system yields a crystallization age for 78238 of 4334 ± 37 Ma and an initial value of −0.27 ± 0.74. The age-initial (T-I) systematics of a variety of KREEP-rich samples, including 78238 and other Mg-suite rocks, KREEP basalts, and olivine cumulate NWA 773, suggest that lunar differentiation was completed by 4492 ± 61 Ma assuming a Chondritic Uniform Reservoir bulk composition for the Moon. The Rb-Sr isotopic systematics of 78238 were disturbed by post-crystallization processes. Nevertheless, selected data points yield two Rb-Sr isochrons. One is concordant with the Sm-Nd crystallization age, 4366 ± 53 Ma. The other is 4003 ± 95 Ma and is concordant with an Ar-Ar age for 78236. The ²⁰⁷Pb-²⁰⁶Pb age of 4333 ± 59 Ma is concordant with the Sm-Nd age. The U-Pb isotopic systematics of 78238 yield linear arrays equivalent to younger ages than the Pb-Pb system, and may reflect fractionation of U and Pb during sample handling. Despite the disturbed nature of the U-Pb systems, a time-averaged μ (²³⁸U/²⁰⁴Pb) value of the source can be estimated at 27 ± 30 from the Pb-Pb isotopic systematics. Because KREEP-rich samples are likely to be derived from source regions with the highest U/Pb ratios, the relatively low μ value calculated for the 78238 source suggests the bulk Moon does not have an exceedingly high μ value.     

Determination of Sr and Ba partition coefficients between apatite from fish (Sparus aurata) and seawater: The influence of temperature

The Sr/Ca and Ba/Ca ratios in inorganic apatite are strongly dependent on the temperature of the aqueous medium during precipitation. If valid in biogenic apatite, these thermometers would offer the advantage of being more resistant to diagenesis than those calibrated on biogenic calcite and aragonite. We have reared seabreams (Sparus aurata) in tanks with controlled conditions during experiments lasting for more than 2 years at 13, 17, 23 and 27 °C, in order to determine the variations in Sr and Ba partitioning relative to Ca (D^Sr and D^Ba, respectively) between seawater and fish apatitic hard tissues (i.e. teeth and bones), as a function of temperature. The sensitivity of the Sr and Ba thermometers (i.e. ∂D^Sr/∂T and ∂D^Ba/∂T, respectively), are similar in bone (/∂T = 0.0036 ± 0.0003 and /∂T = 0.0134 ± 0.0026, respectively) and enamel (/∂T = 0.0037 ± 0.0005 and /∂T = 0.0107 ± 0.0026, respectively). The positive values of ∂D^Sr/∂T and ∂D^Ba/∂T in bone and enamel indicate that D^Sr and D^Ba increase with increasing temperature, a pattern opposite to that observed for inorganic apatite. This distinct thermodependent trace element partitioning between inorganic and organic apatite and water highlights the contradictory effects of the crystal-chemical and biological controls on the partitioning of Ca, Sr and Ba in vertebrate organisms. Taking into account the diet Sr/Ca and Ba/Ca values, it is shown that the bone Ba/Ca signature of fish can be explained by Ca-biopurification and inorganic apatite precipitation, whereas both of these processes fail to predict the bone Sr/Ca values. Therefore, the metabolism of Ca as a function of temperature still needs to be fully understood. However, the biogenic Sr thermometer is used to calculate an average seawater temperature of 30.6 °C using the Sr/Ca compositions of fossil shark teeth at the Cretaceous/Tertiary boundary, and a typical seawater Sr/Ca ratio of 0.02. Finally, while the present work should be completed with data obtained in natural contexts, it is clear that Sr/Ca and Ba/Ca ratios in fossil biogenic apatite already constitute attractive thermometers for marine paleoenvironments.     

Co-precipitation of radium in high ionic strength systems: 2. Kinetic and ionic strength effects

High concentrations of naturally occurring radium pose environmental and health concerns in natural and industrial systems. The adsorption of Ra²⁺ in saline water is limited compared to its adsorption in fresh water, but the process of co-precipitation may be effective in decreasing its concentration. However, despite its importance, Ra co-precipitation has rarely been studied in high ionic strength environments such as those in evaporitic systems.The fate of Ra in the reject brine of a desalination plant was studied via evaporation batch experiments at ionic strengths (I) ranging from 0.7 to 7.0 mol kg⁻¹. Precipitation sequences revealed that Ra co-precipitated with barite, even though the latter was a trace mineral compared to the precipitated gypsum. The concentration-based effective partition coefficient, , for the co-precipitation reaction was 1.04 ± 0.01. This value of is significantly lower than the value for relatively diluted solutions (1.8 ± 0.1). This low value of is mainly the result of a kinetic effect but is also slightly affected by the ionic strength.Both effects are quantitatively examined in the present paper. It is suggested that a kinetic effect influences the nucleation of (Ra,Ba)SO₄, reducing the value of the partition coefficient. This kinetic effect is caused by the favorable nucleation of a more soluble phase (i.e., a phase with a higher BaSO₄ fraction). An additional decrease in the partition coefficient results from the ionic strength effect. Considering the activity of Ra²⁺ and Ba²⁺ in the solution (rather than their concentration) makes it possible to determine the activity-based partition coefficient (), which accounts for the ionic strength effect. was calculated empirically from the experiments and theoretically via a kinetic model. The two derived values are consistent with one another and indicate the combined effect of ionic strength and precipitation kinetics.Finally, the common assumption that γ_Ra²⁺/γ_Ba²⁺=1 was re-examined using a numerical model to predict the experimental results. As the ionic strength increases, this assumption becomes less appropriate for predicting the change in as calculated in the experiments. Understanding the co-precipitation of Ra in such systems is crucial for risk assessments in which both Ra concentration and ionic strength are relatively high.