Experimental study of the incorporation of Li, Sc, Al and other trace elements into olivine

We performed a series of synthesis experiments at 1 atm pressure to investigate the substitution mechanisms of 1+ and 3+ ions into olivine. Forsterite crystals were grown from bulk compositions that contained the element of interest (e.g. Li) and different amounts of additional single trace elements. By working at constant (major element) liquid composition and temperature we eliminated all compositional effects other than those due to the trace elements. Mineral-melt pairs were then analysed to determine the compositional-dependence of the partition coefficient (D), which corresponds to , and where [element] refers to weight concentration of the element in the respective phase.We find that Li forms a stable coupled substitution with Sc and, at above ∼500 ppm Sc in the crystal, Li⁺ and Sc³⁺ ions form an ordered neutral complex ([LiSc]). This complex dissociates at lower trace element concentrations and a second, concentration-independent, mechanism begins to dominate. This second solution mechanism is most likely 2Li⁺ ⇔ Mg²⁺ where one of the Li atoms is in an interstitial position in the crystal lattice. Natural olivines show Li contents slightly greater than Sc (on an atomic basis), indicating that both substitution mechanisms are significant. Unlike Sc, Al does not appear to form a stable complex with Li in the olivine structure.Sodium is highly incompatible in olivine with of ∼0.00015-0.03. Olivine-liquid partitioning of Na⁺ is independent of Sc³⁺ or Al³⁺ concentration. This indicates that the coupled substitution of Na⁺ with any 3+ ions is unlikely. Instead, the relevant substitution mechanism appears to be 2Na⁺ ⇔ Mg²⁺. Although independent of 3+ ion concentration, is inversely correlated with the Li concentration of both melts and crystals, implying that Na competes (unsuccessfully) with Li to replace Mg in the olivine structure.Aluminium is highly incompatible in forsterite . Values of are similar for all phase pairs synthesised from starting materials containing between 10 and 100,000 ppm Al. This suggests that Al is principally incorporated in forsterite by replacing one Mg and one Si atom , where the Al atoms on octahedral (Mg) and tetrahedral (Si) sites are dissociated from one another.The incorporation of gallium into forsterite is influenced by the presence of Li. Where Li concentration in the crystal is much greater than that of Ga (on an atomic basis) we find an excellent correlation between and melt Li content. This relationship indicates that Ga³⁺ and Li⁺ replace 2Mg²⁺ on octahedral sites and that the Ga and Li atoms are, like Sc and Li, strongly associated in the crystal structure.The mechanism by which scandium is incorporated into forsterite is strongly governed by the presence Li. As discussed above, ordered complexes form readily in forsterite in Li-rich experiments. Under Li-absent but Sc-rich conditions (Sc in the crystal >∼500 ppm), is proportional to the concentration of Sc in the melt. This indicates that Sc incorporation is charge-balanced by the formation of magnesium vacancies , and that both species are associated . At lower Sc concentrations (<500 ppm in the crystal), the concentration-dependence of partitioning indicates that the complexes dissociate.Our results demonstrate that partitioning of 1+ and 3+ ions into olivine is complex and involves a range of point defects which yield strongly composition-dependent crystal-melt partition coefficients. Since physical and chemical properties of natural olivine, such as diffusion of ⁶Li and ⁷Li and H₂O solubility, depend on the concentrations of the defects identified in this study, our results provide an important insight into how determining substitution mechanisms can improve our understanding of large-scale mantle processes and properties.     

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.     

First-principles estimates of equilibrium magnesium isotope fractionation in silicate, oxide, carbonate and hexaaquamagnesium(2+) crystals

Equilibrium mass-dependent magnesium isotope fractionation factors are estimated for a range of crystalline compounds including oxides, silicates, carbonates, and salts containing the solvation complex. Fractionation factors for the gas-phase species Mg and MgO are also presented. Fractionation factors are calculated with density functional perturbation theory (DFPT), using norm-conserving pseudopotentials. The results suggest that there will be substantial inter-mineral fractionation, particularly between tetrahedrally coordinated Mg²⁺ in spinel (MgAl₂O₄) and the more common octahedrally coordinated Mg²⁺-sites in silicate and carbonate minerals. Isotope fractionations calculated for Mg²⁺ in hexaaquamagnesium(2+) salts are in good agreement with previous fractionation models of based on large molecular clusters (Black et al., 2007), but show possibly more significant disagreement with a more recent study (Rustad et al., 2010). These models further suggest that solvated , in the form of , will have higher ²⁶Mg/²⁴Mg than coexisting magnesite and dolomite. Calculated fractionations are consistent with Mg-isotope fractionations observed in peridotite mineral separates and inorganic carbonate precipitates. Predicted large, temperature-sensitive spinel-silicate fractionations, in particular, may find use in determining equilibration temperatures of peridotites and other high-temperature rock types.     

The ionic strength dependence of lead (II) carbonate complexation in perchlorate media

Lead speciation in many aqueous geochemical systems is dominated by carbonate complexation. However, direct observations of Pb²⁺ complexation by carbonate ions are few in number. This work represents the first investigation of the equilibrium over a range of ionic strength. Through spectrophotometric observations of formation at 25 °C in NaHCO₃-NaClO₄ solutions, formation constants of the form were determined between 0.001 and 5.0 molal ionic strength. Formation constant results were well represented by the equation:

     

A new value for the stable oxygen isotope fractionation between dissolved sulfate ion and water

Although the stable oxygen isotope fractionation between dissolved sulfate ion and H₂O (hereafter ) is of physico-chemical and biogeochemical significance, no experimental value has been established until present. The primary reason being that uncatalyzed oxygen exchange between and H₂O is extremely slow, taking ∼10⁵ years at room temperature. For lack of a better approach, values of 16‰ and 31‰ at 25 °C have been assumed in the past, based on theoretical ‘gas-phase’ calculations and extrapolation of laboratory results obtained at temperatures >75 °C that actually pertain to the bisulfate system. Here I use novel quantum-chemistry calculations, which take into account detailed solute-water interactions to establish a new value for of 23‰ at 25 °C. The results of the corresponding calculations for the bisulfate ion are in agreement with observations. The new theoretical values show that sediment -data, which reflect oxygen isotope equilibration between sulfate and ambient water during microbial sulfate reduction, are consistent with the abiotic equilibrium between and water.     

Ammonium loss and nitrogen isotopic fractionation in biotite as a function of metamorphic grade in metapelites from western Maine, USA

Ammonium fixed in micas of metamorphic rocks is a sensitive indicator both of organic-inorganic interactions during diagenesis as well as of the devolatilization history and fluid/rock interaction during metamorphism. In this study, a collection of geochemically well-characterized biotite separates from a series of graphite-bearing Paleozoic greenschist- to upper amphibolite-facies metapelites, western Maine, USA, were analyzed for ammonium nitrogen () contents and isotopic composition (δ¹⁵N_NH4) using the HF-digestion distillation technique followed by the EA-IRMS technique. Biotite separates, sampled from 9 individual metamorphic zones, contain 3000 to 100 ppm with a wide range in δ¹⁵N from +1.6‰ to +9.1‰. Average contents in biotite show a distinct decrease from about 2750 ppm for the lowest metamorphic grade (∼500 °C) down to 218 ppm for the highest metamorphic grade (∼685 °C). Decreasing abundances in are inversely correlated in a linear fashion with increasing K⁺ in biotite as a function of metamorphic grade and are interpreted as a devolatilization effect. Despite expected increasing δ¹⁵N_NH4 values in biotite with nitrogen loss, a significant decrease from the Garnet Zones to the Staurolite Zones was found, followed by an increase to the Sillimanite Zones. This pattern for δ¹⁵N_NH4 values in biotite inversely correlates with Mg/(Mg + Fe) ratios in biotite and is discussed in the framework of isotopic fractionation due to different exchange processes between or , reflecting devolatilization history and redox conditions during metamorphism.     

New indices and calibrations derived from the distribution of crenarchaeal isoprenoid tetraether lipids: Implications for past sea surface temperature reconstructions

Several studies have shown that there is a strong relationship between the distribution of crenarchaeotal isoprenoid glycerol dibiphytanyl glycerol tetraethers (GDGTs) and sea surface temperature (SST). Based on this, a ratio of certain GDGTs, called TEX₈₆ (TetraEther indeX of tetraethers consisting of 86 carbon atoms), was developed as a SST proxy. In this study, we determined the distribution of crenarchaeotal isoprenoid GDGTs in 116 core-top sediments mostly from (sub)polar oceans and combined these data with previously published core-top data. Using this extended global core-top dataset (n = 426), we re-assessed the relationship of crenarchaeal isoprenoid GDGTs with SST. We excluded data from the Red Sea from the global core-top dataset to define new indices and calibration models, as the Red Sea with its elevated salinity appeared to behave differently compared to other parts of the oceans. We tested our new indices and calibration models on three different paleo datasets, representing different temperature ranges. Our results indicate that the crenarchaeol regio-isomer plays a more important role for temperature adaptation in (sub)tropical oceans than in (sub)polar oceans, suggesting that there may be differences in membrane adaptation of the resident crenarchaeotal communities at different temperatures. We, therefore, suggest to apply two different calibration models. For the whole calibration temperature range (−3 to 30 °C), a modified version of TEX₈₆ with a logarithmic function which does not include the crenarchaeol regio-isomer, called , is shown to correlate best with SST: (r² = 0.86, n=396, p <0.0001). Application of on sediments from the subpolar Southern Ocean results in realistic absolute SST estimates and a similar SST trend compared to a diatom SST record from the same core. , which is defined as the logarithmic function of TEX₈₆, yields the best correlation with SST, when the data from the (sub)polar oceans are removed: (r² = 00.87, n = 255, p < 0.0001). Furthermore, gives the best correlation for mescosm data with temperatures ranging between 10 and 46 °C. For Quaternary sediments from the tropical Arabian Sea, both and yield similar trends and SST estimates. However, the extrapolation of calibration on a sediment record from a greenhouse world ocean predicts more reliable absolute SST estimates and relative SST changes in agreement with estimates based on the δ¹⁸O of planktonic foraminifera. Based on the comparison of and derived SSTs using the core top data, we recommend applying above 15 °C and below 15 °C. In cases where paleorecords encompass temperatures both below and above 15 °C, we suggest to use .     

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.     

Comparative carbonate and hydroxide complexation of cations in seawater

Comparative concentrations of carbonate and hydroxide complexes in natural solutions can be expressed in terms of reactions with bicarbonate that have no explicit pH dependence (). Stability constants for this reaction with n = 1 were determined using conventional formation constant data expressed in terms of hydroxide and carbonate. Available data indicate that stability constants appropriate to seawater at 25 °C expressed in the form are on the order of 10^4.2 for a wide range of cations (M^z+) with z = +1, +2 and +3. Φ₁ is sufficiently large that species appear to substantially dominate MOH^z−1 species in seawater. Evaluations of comparative stepwise carbonate and hydroxide stability constant behavior leading to the formation of n = 2 and n = 3 complexes suggest that carbonate complexes generally dominate hydroxide complexes in seawater, even for cations whose inorganic speciation schemes in seawater are currently presumed to be strongly dominated by hydrolyzed forms (). Calculated stability constants, and , indicate that the importance of carbonate complexation is sufficiently large that carbonate and hydroxide complexes would be generally comparable even if calculated Φ₂ and Φ₃ values are overestimated by two or more orders of magnitude. Inclusion of mixed ligand species in carbonate-hydroxide speciation models allows cation complexation intensities (M_T/[M^z+]) to be expressed in the following form:

     

Effects of magnesium ions on near-equilibrium calcite dissolution: Step kinetics and morphology

Dissolution kinetics at the aqueous solution-calcite interface at 50 °C were investigated using in situ atomic force microscopy (AFM) to reveal the influence of magnesium concentration and solution saturation state on calcite dissolution kinetics and surface morphology. Under near-equilibrium conditions, dissolved Mg²⁺ displayed negligible inhibitory effects on calcite dissolution even at concentrations of . Upon the introduction of , the solution saturation state with respect to calcite, , acted as a “switch” for magnesium inhibition whereby no significant changes in step kinetics were observed at Ω_calcite<0.2, whereas a sudden inhibition from Mg²⁺ was activated at Ω_calcite?0.2. The presence of the Ω-switch in dissolution kinetics indicates the presence of critical undersaturation in accordance with thermodynamic principles. The etch pits formed in solutions with exhibited a unique distorted rhombic profile, different from those formed in Mg-free solutions and in de-ionized water. Such unique etch pit morphology may be associated with the anisotropy in net detachment rates of counter-propagating kink sites upon the addition of Mg²⁺.     

Sorption of phosphate onto calcite; results from batch experiments and surface complexation modeling

The adsorption of phosphate onto calcite was studied in a series of batch experiments. To avoid the precipitation of phosphate-containing minerals the experiments were conducted using a short reaction time (3 h) and low concentrations of phosphate (?50 μM). Sorption of phosphate on calcite was studied in 11 different calcite-equilibrated solutions that varied in pH, P_CO2, ionic strength and activity of Ca²⁺, and . Our results show strong sorption of phosphate onto calcite. The kinetics of phosphate sorption onto calcite are fast; adsorption is complete within 2-3 h while desorption is complete in less than 0.5 h. The reversibility of the sorption process indicates that phosphate is not incorporated into the calcite crystal lattice under our experimental conditions. Precipitation of phosphate-containing phases does not seem to take place in systems with ?50 μM total phosphate, in spite of a high degree of super-saturation with respect to hydroxyapatite (SI_HAP ? 7.83). The amount of phosphate adsorbed varied with the solution composition, in particular, adsorption increases as the activity decreases (at constant pH) and as pH increases (at constant activity). The primary effect of ionic strength on phosphate sorption onto calcite is its influence on the activity of the different aqueous phosphate species. The experimental results were modeled satisfactorily using the constant capacitance model with >CaPO₄Ca⁰ and either >CaHPO₄Ca⁺ or > as the adsorbed surface species. Generally the model captures the variation in phosphate adsorption onto calcite as a function of solution composition, though it was necessary to include two types of sorption sites (strong and weak) in the model to reproduce the convex shape of the sorption isotherms.     

The influence of temperature and seawater composition on calcite crystal growth mechanisms and kinetics: Implications for Mg incorporation in calcite lattice

The composition of carbonate minerals formed in past and present oceans is assumed to be significantly controlled by temperature and seawater composition. To determine if and how temperature is kinetically responsible for the amount of Mg incorporated in calcite, we quantified the influence of temperature and specific dissolved components on the complex mechanism of calcite precipitation in seawater. A kinetic study was carried out in artificial seawater and NaCl-CaCl₂ solutions, each having a total ionic strength of 0.7 M. The constant addition technique was used to maintain [Ca²⁺] at 10.5 mmol kg⁻¹ while [] was varied to isolate the role of this variable on the precipitation rate of calcite.Our results show that the overall reaction of calcite precipitation in both seawater and NaCl-CaCl₂ solutions is dominated by the following reaction:

     

Sphalerite oxidation pathways detected by oxygen and sulfur isotope studies

Sphalerite oxidation is a common process under acid-mine drainage (AMD) conditions and results in the release of , Zn and potentially toxic trace metals, which can pollute rivers and oceans. However, there are only a few studies on the mechanisms of aerobic sphalerite oxidation. Oxygen and S isotope investigations of the produced may contribute to the understanding of sphalerite oxidation mechanisms so helping to interpret field data from AMD sites. Therefore, batch oxidation experiments with an Fe-rich sphalerite were performed under aerobic abiotic conditions at different initial pH values (2 and 6) for different lengths of time (2–100 days). The O and S isotope composition of the produced indicated changing oxidation pathways during the experiments. During the first 20 days of the experiments at both initial pH values, molecular O₂ was the exclusive O source of . Furthermore, the lack of S isotope enrichment processes between and sphalerite indicated that O₂ was the electron acceptor from sphalerite S. As the oxidation proceeded, a sufficient amount of released Fe(II) was oxidized to Fe(III) by O₂. Therefore, electrons could be transferred from sphalerite S sites to adsorbed hydrous Fe(III) and O from the hydration sphere of Fe was incorporated into the produced as indicated by decreasing δ¹⁸O_SO4 values which became more similar to the δ¹⁸O_H2O values. The enrichment of ³²S in relative to the sphalerite may also result from sphalerite oxidation by Fe(III).The incorporation of O₂ into during the oxidation of sphalerite was associated with an O isotope enrichment factor ε_SO4–O2 of ca. −22‰. The O isotope enrichment factor ε_SO4–H2O was determined to be ?4.1‰. A comparison with O and S studies of other sulfides suggests that there is no general oxidation mechanism for acid-soluble sulfides.     

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

Confirmation of mass-independent Ni isotopic variability in iron meteorites

We report high-precision analyses of internally-normalised Ni isotope ratios in 12 bulk iron meteorites. Our measurements of ⁶⁰Ni/⁶¹Ni, ⁶²Ni/⁶¹Ni and ⁶⁴Ni/⁶¹Ni normalised to ⁵⁸Ni/⁶¹Ni and expressed in parts per ten thousand (?) relative to NIST SRM 986 as and , vary by 0.146, 0.228 and 0.687, respectively. The precision on a typical analysis is 0.03?, 0.05? and 0.08? for , and , respectively, which is comparable to our sample reproducibility. We show that this ‘mass-independent’ Ni isotope variability cannot be ascribed to interferences, inaccurate correction of instrumental or natural mass-dependent fractionation, fractionation controlled by nuclear field shift effects, nor the influence of cosmic ray spallation. These results thus document the presence of mass-independent Ni isotopic heterogeneity in bulk meteoritic samples, as previously proposed by Regelous et al. (2008) (EPSL 272, 330-338), but our new analyses are more precise and include determination of ⁶⁴Ni. Intriguingly, we find that terrestrial materials do not yield homogenous internally-normalised Ni isotope compositions, which, as pointed out by Young et al. (2002) (GCA 66, 1095-1104), may be the expected result of using the exponential (kinetic) law and atomic masses to normalise all fractionation processes. The certified Ni isotope reference material NIST SRM 986 defines zero in this study, while appropriate ratios for the bulk silicate Earth are given by the peridotites JP-1 and DTS-2 and, relative to NIST SRM 986, yield deviations in , and of −0.006?, 0.036? and 0.119?, respectively. There is a strong positive correlation between and in iron meteorites analyses, with a slope of 3.03 ± 0.71. The variations of Ni isotope anomalies in iron meteorites are consistent with heterogeneous distribution of a nucleosynthetic component from a type Ia supernova into the proto-solar nebula.     

Molybdenum, tungsten and manganese partitioning in the system pyrrhotite-Fe-S-O melt-rhyolite melt: Impact of sulfide segregation on arc magma evolution

Experiments were performed to determine the partitioning of molybdenum, tungsten and manganese among a rhyolitic melt (melt), pyrrhotite (po), and an immiscible Fe-S-O melt (Fe-S-O). Sulfide phases such as these may be isolated from a silicate melt along with other crystallizing phases during the evolution of arc magma, and partition coefficients are required to model the effect of this process on molybdenum and tungsten budgets.We developed an experimental design to take advantage of properties of the phases under study. Careful control of temperature allowed pyrrhotite and magnetite to be stable along with an Fe-S-O melt, and this phase assemblage allowed the composition of run-product pyrrhotite to be used to calculate both fS₂ and fO₂ for the experiments. At run temperature, (1042 ± 2 °C), a rhyolitic melt can be formed at low pressure, under nominally dry conditions, which removed the need for confining pressure as well as externally imposed fugacities. The silica-saturated melt allowed the charges to be contained in sealed evacuated silica tubes without danger of reaction, and with closed system behavior for molybdenum and tungsten.Experiments were run for durations up to 2000 min. Molybdenite (mb) and wolframite (wo) were added to the experiments as sources for molybdenum and tungsten, respectively. Manganese was added to the system as both a component of the starting rhyolitic pumice, and of Mn-bearing wolframite. Oxygen fugacity in these experiments was fixed at the Ni-NiO oxygen fugacity buffer. Sulfur fugacity was 10⁻¹ bar. Run products were analyzed by EPMA and LA-ICP-MS. Analysis of the run products yielded ( standard deviation of the mean): , , , and . The partition coefficients for manganese in this system are and .Simple Rayleigh fractionation modeling suggests that oxidized felsic melts produced through fractional crystallization may have lost as much as 14% of their initial molybdenum, but only 2% of their initial tungsten, through the removal of an Fe-S-O melt along with crystalline phases. Modeling consistent with conditions of oxygen and sulfur fugacity influenced by assimilation of sulfide (with low concentrations of molybdenum and tungsten) from, for example, sedimentary rock, results in evolved magmas significantly depleted in molybdenum, but only moderately depleted in tungsten. The molybdenum:tungsten ratio can vary by two orders of magnitude. These systematics may help to explain some of the variability in metal ratios of intrusion-related hydrothermal ore deposits.     

Kinetics of gypsum crystal growth from high ionic strength solutions: A case study of Dead Sea - seawater mixtures

Gypsum precipitation kinetics were examined from a wide range of chemical compositions , ionic strengths (4.75-10 m) and saturation state with respect to gypsum (1.16-1.74) in seeded batch experiments of mixtures of Ca²⁺-rich Dead Sea brine and -rich seawater. Despite the variability in the experimental solutions, a single general rate law was formulated to describe the heterogeneous precipitation rate of gypsum from these mixtures:

     

Elemental and isotopic abundances of lithium in chondrule constituents in the Allende meteorite

Elemental and isotopic abundances of lithium in chondrule constituents in the Allende CV3 meteorite were determined using secondary ion mass spectrometry. Olivines and mesostasis dominated by a feldspathic phase are depleted in Li ( and , respectively). In contrast, low-Ca pyroxenes and mesostasis dominated by a Na-rich phase are enriched in Li ( and , respectively) and the interchondrule matrix is generally enriched in Li ( on average). The Li isotopic abundance of olivine ranges from to 21. The spatial distributions of elemental and isotopic abundances of Li in olivines within individual chondrules exhibit no systematic pattern. This suggests that the distribution of Li in olivine was not disturbed during aqueous alteration or thermal metamorphism on the Allende meteorite parent body. Although mesostasis is the last crystallizing phase from a chondrule melt and is expected to be enriched in Li, in the Allende meteorite it is generally depleted in Li. We suggest that during aqueous alteration on the CV asteroid, Li in mesostasis was leached out by aqueous fluids. The Li-enriched Na-rich mesostasis was probably produced later by infiltration of Na-rich fluids. It seems likely that aqueous fluids sequestered alkali elements from the Allende-chondrite region in the CV parent asteroid, although significant amounts of Li are preserved in ferrous olivine in the interchondrule matrix.     

An experimental study of Cobalt (II) complexation in Cl and H2S-bearing hydrothermal solutions

The speciation of cobalt (II) in Cl⁻ and H₂S-bearing solutions was investigated spectrophotometrically at temperatures of 200, 250, and 300 °C and a pressure of 100 bars, and by measuring the solubility of cobaltpentlandite at temperatures of 120-300 °C and variable pressures of H₂S. From the results of these experiments, it is evident that CoHS⁺ and predominate in the solutions except at 150 °C, for which the dominant chloride complex is CoCl₃⁻. The logarithms of the stability constant for CoHS⁺ show moderate variation with temperature, decreasing from 6.24 at 120 °C to 5.84 at 200 °C, and increasing to 6.52 at 300 °C. Formation constants for chloride species increase smoothly with temperature and at 300°C their logarithms reach 8.33 for , 6.44 for CoCl₃⁻, 4.94 to 5.36 for , and 2.42 for CoCl⁺. Calculations based on the composition of a model hydrothermal fluid (Ksp-Mu-Qz, KCl = 0.25 m, NaCl = 0.75 m, ΣS = 0.3 m) suggest that at temperatures ?200 °C, cobalt occurs dominantly as CoHS⁺, whereas at higher temperatures the dominant species is .     

Boron isotopes as pH proxy: A new look at boron speciation in deep-sea corals using B MAS NMR and EELS

Dissolved boron in modern seawater occurs in the form of two species, trigonal boric acid B(OH)₃ and tetrahedral borate ion . One of the key assumption in the use of boron isotopic compositions of carbonates as pH proxy is that only borate ions, , are incorporated into the carbonate. Here, we investigate the speciation of boron in deep-sea coral microstructures (Lophelia pertusa specimen) by using high field magic angle spinning nuclear magnetic resonance (¹¹B MAS NMR) and electron energy-loss spectroscopy (EELS). We observe both boron coordination species, but in different proportions depending on the coral microstructure, i.e. centres of calcification versus fibres. These results suggest that careful sampling is necessary before performing boron isotopic measurements in deep-sea corals. By combining the proportions of B(OH)₃ and determined by NMR and our previous ion microprobe boron isotope measurements, we propose a new equation for the relation between seawater pH and boron isotopic composition in deep-sea corals.