Low solubility of He and Ar in carbonatitic liquids: Implications for decoupling noble gas and lithophile isotope systems

In order to asses the importance of carbonatitic liquids in transporting noble gases in the mantle, the solubilities of He and Ar in carbonatitic liquids were estimated from analyses of calcium-potassium carbonate glasses that had been synthesized at 1 bar and temperatures between 850 and 950 °C under He or Ar enriched atmospheres. Despite poor reproducibility related to difficulties synthesizing carbonatite glass, we have been able to estimate He and Ar solubilities in carbonatite liquids to be 1 × 10⁻⁸ and 2 × 10⁻⁹ mol g⁻¹ at 1 bar respectively (with ?50% uncertainty). Despite the significant uncertainties on these estimates, it is clear that the noble gases are not massively soluble in carbonatite liquids (within error, these solubilities are identical to their equivalent solubilities in tholeiitic melts). Assuming the results of these low pressure experiments can be applied to mantle conditions, it seems unlikely that carbonatite metasomatism per se transports noble gases within the mantle. It is nevertheless possible that partitioning of lithophile trace elements (including the important radioelements, U, K and Th) and noble gases between a carbonatitic melt and a silicate melt could effectively decouple lithophile and noble gas isotope systematics because the carbonatitic melt expressedly does not transport noble gases, yet is known to efficiently transport incompatible trace elements.     

Gold solubility in oxidized and reduced, water-saturated mafic melt

We have performed experiments to evaluate Au solubility in natural, water-saturated basaltic melts as a function of oxygen fugacity. Experiments were carried out at 1000 °C and 200 MPa, and oxygen fugacity was controlled at the fayalite-magnetite-quartz (FMQ) oxygen fugacity buffer and FMQ + 4. All experiments were saturated with a metal-chloride aqueous solution loaded initially as a 10 wt% NaCl eq. fluid. The stable phase assemblage at FMQ consists of basalt melt, olivine, clinopyroxene, a single-phase aqueous fluid, and metallic Au. The stable phase assemblage at FMQ + 4 consists of basalt melt, clinopyroxene, magnetite-spinel solid solution, a single-phase aqueous fluid, and metallic Au. Silicate glasses (i.e., quenched melt) and their contained crystalline material were analyzed by using both electron probe microanalysis (EPMA) and laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS). Measured Au concentrations in the quenched melt range from 4.8 μg g⁻¹ to 0.64 μg g⁻¹ at FMQ + 4, and 0.54 μg g⁻¹ to 0.1 μg g⁻¹ at FMQ. The measured solubility of Au in olivine and clinopyroxene was consistently below the LA-ICP-MS limit of detection (i.e., 0.1 μg g⁻¹). These melt solubility data place important limitations on the dissolved Au content of water-saturated, Cl- and S-bearing basaltic liquids at geologically relevant fO₂ values. The new data are compared to published, experimentally-determined values for Au solubility in dry and hydrous silicate liquids spanning the compositional range from basalt to rhyolite, and the effects of melt composition, oxygen fugacity, pressure and temperature are discussed.     

Na2O solubility in CaO-MgO-SiO2 melts

The sodium solubility in silicate melts in the CaO-MgO-SiO₂ (CMS) system at 1400 °C has been measured by using a closed thermochemical reactor designed to control alkali metal activity. In this reactor, Na_(g) evaporation from a Na₂O-xSiO₂ melt imposes an alkali metal vapor pressure in equilibrium with the molten silicate samples. Because of equilibrium conditions in the reactor, the activity of sodium-metal oxide in the molten samples is the same as that of the source, i.e., aNa₂O_(sample) = aNa₂O_(source). This design also allows to determine the sodium oxide activity coefficient in the samples. Thirty-three different CMS compositions were studied. The results show that the amount of sodium entering from the gas phase (i.e., Na₂O solubility) is strongly sensitive to silica content of the melt and, to a lesser extent, the relative amounts of CaO and MgO. Despite the large range of tested melt compositions (0 < CaO and MgO < 40; 40 < SiO₂ < 100; in wt%), we found that Na₂O solubility is conveniently modeled as a linear function of the optical basicity (Λ) calculated on a Na-free basis melt composition. In our experiments, γNa₂O_(sample) ranges from 7 × 10⁻⁷ to 5 × 10⁻⁶, indicating a strongly non-ideal behavior of Na₂O solubility in the studied CMS melts (γNa₂O_(sample) ? 1). In addition to showing the effect of sodium on phase relationships in the CMS system, this Na₂O solubility study brings valuable new constraints on how melt structure controls the solubility of Na in the CMS silicate melts. Our results suggest that Na₂O addition causes depolymerization of the melt by preferential breaking of Si-O-Si bonds of the most polymerized tetrahedral sites, mainly Q⁴.     

Temperature independent thermal expansivities of calcium aluminosilicate melts between 1150 and 1973 K in the system anorthite-wollastonite-gehlenite (An-Wo-Geh): A density model

The thermal expansivities of 10 compositions from within the anorthite-wollastonite-gehlenite (An-Wo-Geh) compatibility triangle have been investigated using a combination of calorimetry and dilatometry on the glassy and liquid samples. The volumes at room temperature were derived from densities measured using the Archimedean buoyancy method. For each sample, density was measured at 298 K using glass that had a cooling-heating history of 10-10 K min⁻¹. The thermal expansion coefficient of the glass from 298 K to the glass transition interval was measured by a dilatometer and the heat capacity was measured using a differential scanning calorimeter from 298 to 1135 K. The thermal expansion coefficient and the heat flow were determined at a heating rate of 10 K min⁻¹ on glasses which were previously cooled at 10 K min⁻¹. Supercooled liquid density, molar volume and molar thermal expansivities were indirectly determined by combining differential scanning calorimetric and dilatometric measurements assuming that the kinetics of enthalpy and shear relaxation are equivalent. The data obtained on supercooled liquids were compared to high-temperature predictions from the models of (Lange, R.A., Carmichael, I.S.E., 1987. Densities of Na₂O-K₂O-CaO-MgO-FeO-Fe₂O₃-Al₂O₃-TiO₂-SiO₂ liquids: New measurements and derived partial molar properties. Geochim. Cosmochim. Acta51, 2931-2946; Courtial, P., Dingwell, D.B., 1995. Nonlinear composition dependence of molar volume of melts in the CaO-Al₂O₃-SiO₂ system. Geochim. Cosmochim. Acta59 (18), 3685-3695; Lange, R.A., 1997. A revised model for the density and thermal expansivity of K₂O-Na₂O-CaO-MgO-Al₂O₃-SiO₂ liquids from 700 to 1900 K: extension to crustal magmatic temperatures. Contrib. Mineral. Petrol.130, 1-11). The best linear fit combines the supercooled liquid data presented in this study and the high temperature data calculated using the Courtial and Dingwell (1995) model. This dilatometric/calorimetric method of determining supercooled liquid molar thermal expansivity greatly increases the temperature range accessible for thermal expansion. It represents a substantial increase in precision and understanding of the thermodynamics of calcium aluminosilicate melts. This enhanced precision demonstrates clearly the temperature independence of the melt expansions in the An-Wo-Geh system. This contrasts strongly with observations for neighboring system such as anorthite-diopside and raises the question of the compositional/structural origins of temperature dependence of thermal expansivity in multicomponent silicate melts.     

Olivine/melt transition metal partitioning, melt composition, and melt structure—Melt polymerization and Q-speciation in alkaline earth silicate systems

The two most abundant network-modifying cations in magmatic liquids are Ca²⁺ and Mg²⁺. To evaluate the influence of melt structure on exchange of Ca²⁺ and Mg²⁺ with other geochemically important divalent cations (m-cations) between coexisting minerals and melts, high-temperature (1470-1650 °C), ambient-pressure (0.1 MPa) forsterite/melt partitioning experiments were carried out in the system Mg₂SiO₄-CaMgSi₂O₆-SiO₂ with ?1 wt% m-cations (Mn²⁺, Co²⁺, and Ni²⁺) substituting for Ca²⁺ and Mg²⁺. The bulk melt NBO/Si-range (NBO/Si: nonbridging oxygen per silicon) of melt in equilibrium with forsterite was between 1.89 and 2.74. In this NBO/Si-range, the NBO/Si(Ca) (fraction of nonbridging oxygens, NBO, that form bonds with Ca²⁺, Ca²⁺-NBO) is linearly related to NBO/Si, whereas fraction of Mg²⁺-NBO bonds is essentially independent of NBO/Si. For individual m-cations, rate of change of K_D(m−Mg) with NBO/Si(Ca) for the exchange equilibrium, m^melt + Mg^olivine ? m^olivine + Mg^melt, is linear. K_D(m−Mg) decreases as an exponential function of increasing ionic potential, Z/r² (Z: formal electrical charge, r: ionic radius—here calculated with oxygen in sixfold coordination around the divalent cations) of the m-cation. The enthalpy change of the exchange equilibrium, ΔH, decreases linearly with increasing Z/r² [ΔH = 261(9)-81(3)·Z/r² (Å⁻²)]. From existing information on (Ca,Mg)O-SiO₂ melt structure at ambient pressure, these relationships are understood by considering the exchange of divalent cations that form bonds with nonbridging oxygen in individual Qⁿ-species in the melts. The negative ∂K_D(m−Mg)/∂(Z/r²) and ∂(ΔH)/∂(Z/r²) is because increasing Z/r² is because the cations forming bonds with nonbridging oxygen in increasingly depolymerized Qⁿ-species where steric hindrance is decreasingly important. In other words, principles of ionic size/site mismatch commonly observed for trace and minor elements in crystals, also govern their solubility behavior in silicate melts.     

Canadian Shield brine from the Con Mine, Yellowknife, NT, Canada: Noble gas evidence for an evaporated Palaeozoic seawater origin mixed with glacial meltwater and Holocene recharge

Dissolved noble gas concentrations were measured in high salinity (270 g/L) Ca(Na)-Cl groundwaters from the Con Mine, Yellowknife, Canada in an effort to discriminate between two possible origins, as either a brine generated by evaporative enrichment in a Paleozoic inland sea, or marine water concentrated by freezing during glacial times. Major ion and isotope geochemistry indicate that brines from the deepest level remain relatively undisturbed by mixing with modern water introduced by mining. Mixing calculations are used to quantify fractions of brine, glacial meltwater and modern water. From this, noble gas concentrations were corrected for excess air with Ne and normalized to 100% brine solution. Over-pressuring of helium and argon in the brine provide age constraints based on the accumulation of geogenic ⁴He and ⁴⁰Ar. Radiogenic age calculations together with the local geological history suggest brine emplacement during early Palaeozoic time, likely during the Devonian when evaporitic inland seas existed in this region. The concentrations of the atmospherically derived noble gases in the brine fraction (Kr = 1.4E-8, Xe = 8.5E-10 cc_STP/cc_H2O) are close to atmospheric equilibrium for brine at 25 °C (Kr = 7.3E-9, Xe = 8.0E-10 cc_STP/cc_H2O), but are far lower than would be expected for closed-system concentration of seawater by freezing (Kr = 2.8E-6, Xe = 4.2E-7 cc_STP/cc_H2O). Thus, despite the complicated mixing history of the brine, the atmospheric and geogenic noble gases provide strong evidence for an origin as air-equilibrated brine from evaporated Paleozoic seawater, which infiltrated via density displacement through existing fractures and faults into the Canadian Shield.     

H2S fluxes from Mt. Etna, Stromboli, and Vulcano (Italy) and implications for the sulfur budget at volcanoes

We present here new measurements of sulfur dioxide and hydrogen sulfide emissions from Vulcano, Etna, and Stromboli (Italy), made by direct sampling at vents and by filter pack and ultraviolet spectroscopy in downwind plumes. Measurements at the F0 and FA fumaroles on Vulcano yielded SO₂/H₂S molar ratios of ≈0.38 and ≈1.4, respectively, from which we estimate an H₂S flux of 6 to 9 t · d⁻¹ for the summit crater. For Mt. Etna and Stromboli, we found SO₂/H₂S molar ratios of ≈20 and ≈15, respectively, which combined with SO₂ flux measurements, suggest H₂S emission rates of 50 to 113 t · d⁻¹ and 4 to 8 t · d⁻¹, respectively. We observe that “source” and plume SO₂/H₂S ratios at Vulcano are similar, suggesting that hydrogen sulfide is essentially inert on timescales of seconds to minutes. This finding has important implications for estimates of volcanic total sulfur budget at volcanoes since most existing measurements do not account for H₂S emission.     

Solubility products of amorphous ferric arsenate and crystalline scorodite (FeAsO4 · 2H2O) and their application to arsenic behavior in buried mine tailings

Published solubility data for amorphous ferric arsenate and scorodite have been reevaluated using the geochemical code PHREEQC with a modified thermodynamic database for the arsenic species. Solubility product calculations have emphasized measurements obtained under conditions of congruent dissolution of ferric arsenate (pH < 3), and have taken into account ion activity coefficients, and ferric hydroxide, ferric sulfate, and ferric arsenate complexes which have association constants of 10^4.04 (FeH₂AsO₄²⁺), 10^9.86 (FeHAsO₄⁺), and 10^18.9 (FeAsO₄). Derived solubility products of amorphous ferric arsenate and crystalline scorodite (as log K_sp) are −23.0 ± 0.3 and −25.83 ± 0.07, respectively, at 25 °C and 1 bar pressure. In an application of the solubility results, acid raffinate solutions (molar Fe/As = 3.6) from the JEB uranium mill at McClean Lake in northern Saskatchewan were neutralized with lime to pH 2-8. Poorly crystalline scorodite precipitated below pH 3, removing perhaps 98% of the As(V) from solution, with ferric oxyhydroxide (FO) phases precipitated starting between pH 2 and 3. Between pH 2.18 and 7.37, the apparent log K_sp of ferric arsenate decreased from −22.80 to −24.67, while that of FO (as Fe(OH)₃) increased from −39.49 to −33.5. Adsorption of As(V) by FO can also explain the decrease in the small amounts of As(V)(aq) that remain in solution above pH 2-3. The same general As(V) behavior is observed in the pore waters of neutralized tailings buried for 5 yr at depths of up to 32 m in the JEB tailings management facility (TMF), where arsenic in the pore water decreases to 1-2 mg/L with increasing age and depth. In the TMF, average apparent log K_sp values for ferric arsenate and ferric hydroxide are −25.74 ± 0.88 and −37.03 ± 0.58, respectively. In the laboratory tests and in the TMF, the increasing crystallinity of scorodite and the amorphous character of the coexisting FO phase increases the stability field of scorodite relative to that of the FO to near-neutral pH values. The kinetic inability of amorphous FO to crystallize probably results from the presence of high concentrations of sulfate and arsenate.     

Geochemical controls on sediment reactivity and buffering processes in a heterogeneous aquifer

The injection and recovery of oxic water into deep anoxic aquifers may help to alleviate short- and long-term imbalance between freshwater supply and demand. The extent and structure of physical and geochemical heterogeneity of the aquifer will impact the water quality evolution during injection, storage and recovery. Water–sediment interactions within the most permeable parts of the aquifer, where the bulk of the injectant will penetrate, may dominate, however, water quality may also be impacted by interactions within the finer-grained, less permeable but potentially highly reactive media. In this study, the heterogeneity of the reductive capacity of an aquifer selected for water reuse projects was characterised, the amount, type and reactivity of the sedimentary reductants present determined, and the relationship between reductive capacity and sedimentary lithologies quantified. The average potential reductive capacities (PRC_TOT), based on total organic C and pyrite concentrations of the sediment, were quantified for sands (382 μmol O₂ g⁻¹), clays (1522 μmol O₂ g⁻¹), and silts (1957 μmol O₂ g⁻¹). Twenty-seven samples, spanning the three different lithologies, were then incubated for 50 days and the measured reductive capacities (MRC) determined for the sands (29.2 μmol O₂ g⁻¹), silts (136 μmol O₂ g⁻¹), and clays (143 μmol O₂ g⁻¹). On average, the MRC were 10% of the PRC_TOT. The main consumers of O₂ were pyrite (20–100%), sedimentary organic matter (SOM; 3–56%), siderite (3–28%) and Fe(II)-aluminosilicates (8–55%). The incubation data plus hydrogeochemical modelling, indicated that pH-buffering was controlled firstly by dissolution of trace level carbonates, followed by dissolution of feldspars. Zinc, Co, Ni, Cd and Pb were readily mobilized during incubation.     

Surface chemistry of disordered mackinawite (FeS)

Disordered mackinawite, FeS, is the first formed iron sulfide in ambient sulfidic environments and has a highly reactive surface. In this study, the solubility and surface chemistry of FeS is described. Its solubility in the neutral pH range can be described by K_s^app = {Fe²⁺} · {H₂S(aq)} · {H⁺}⁻² = 10^+4.87±0.27. Acid-base titrations show that the point of zero charge (PZC) of disordered mackinawite lies at pH ∼7.5. The hydrated disordered mackinawite surface can be best described by strongly acidic mono-coordinated and weakly acidic tricoordinated sulfurs. The mono-coordinated sulfur site determines the acid-base properties at pH < PZC and has a concentration of 1.2 × 10⁻³ mol/g FeS. At higher pH, the tricoordinated sulfur, which has a concentration of 1.2 × 10⁻³ mol/g FeS, determines surface charge changes. Total site density is 4 sites nm⁻². The acid-base titration data are used to develop a surface complexation model for the surface chemistry of FeS.     

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.     

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.     

Hydration of rhyolitic glass during weathering as characterized by IR microspectroscopy

The mechanism and rate of hydration of rhyolitic glass during weathering were studied. Doubly polished thin sections of two rhyolites with different duration of weathering (Ohsawa lava: 26,000 yr, Awanomikoto lava: 52,000 yr) were prepared. Optical microscope observation showed that altered layers had developed along the glass surfaces. IR spectral line profile analysis was conducted on the glass sections from the surface to the interior for a length of 250 μm and the contents of molecular H₂O (H₂O_m), OH species (OH) and total water (H₂O_t) were determined. The diffusion profile of H₂O_m in Ohsawa lava extends beyond the layer observed by optical microscope. The content of H₂O_m in the hydrated region is much higher than that of OH species. Thus, the reaction from H₂O_m to OH appears to be little and H₂O_m is the dominant water species moving into the glass during weathering. Based on the concentration profiles, the diffusion coefficients of H₂O_m(D_H2Om) and H₂O_t(D_H2Ot) were determined to be 2.8 × 10⁻¹⁰ and 3.4 × 10⁻¹⁰ μm² s⁻¹ for Ohsawa lava, and 5.2 × 10⁻¹¹ and 4.1 × 10⁻¹¹ μm² s⁻¹ for Awanomikoto lava, respectively. The obtained D_H2Om during weathering are more than 2-3 orders of magnitude larger than the diffusion coefficient at ∼20 °C that is extrapolated from the diffusivity data for >400 °C. This might suggest that the mechanism of water transport is different at weathering conditions and >400 °C.     

Diffusion of HTO,Cl and I in Upper Toarcian argillite samples from Tournemire: Effects of initial iodide concentration and ionic strength

Diffusion parameters for HTO, ³⁶Cl⁻_, and ¹²⁵I⁻ were determined on Upper Toarcian argillite samples from the Tournemire Underground Research Laboratory (Aveyron, France) using the through diffusion technique. The direction of diffusion was parallel to the bedding plane. The purpose of the present study was 3-fold; it was intended (i) to confirm the I⁻ interaction with Upper Toarcian argillite and to verify the effects of initial I⁻ concentration on this affinity, as previously observed by means of radial diffusion experiments, (ii) to highlight any discrepancy between Cl⁻ and I⁻ diffusivity, and (iii) to investigate the effect of an increase of the ionic strength of the solution on the anionic tracers’ diffusive behaviour. The results show that the effective diffusion coefficient (D_e) and diffusion accessible porosity (ε_a) values obtained with an ionic strength (I.S.) synthetic pore water of 0.01 eq L⁻¹ are: D_e = 2.35–2.50 × 10⁻¹¹ m² s⁻¹ and ε_a = 12.0–15.0% for HTO, and D_e = 14.5–15.5 × 10⁻¹³ m² s⁻¹ and ε_a = 2.5–2.9% for ³⁶Cl⁻. Because of anionic exclusion effects, anions diffuse slower and exhibit smaller diffusion accessible porosities than HTO, taken as a water tracer. The associated effective diffusion coefficient (D_e) and rock capacity factor (α) obtained for ¹²⁵I are: D_e = 7.00–8.60 × 10⁻¹³ m² s⁻¹ and α = 4.3–7.2%. Such values make it possible to calculate low ¹²⁵I distribution ratios (0.0057 < R_D < 0.0192 mL g⁻¹) which confirm the trend indicating that the ¹²⁵I rock capacity factor increases with the decrease of the initial I⁻ concentration. Additional through-diffusion experiments were carried out with a higher ionic strength synthetic pore water (I.S. = 0.11 eq L⁻¹). No evolution of HTO diffusion parameters was observed. The anionic tracers’ effective diffusion coefficient increased by a factor of two but no clear evolution of their accessible porosity was observed. Such a paradox could be related to the particularly small mean pore size of the Upper Toarcian argillite of Tournemire. The most significant finding of this study is the large discrepancy (factor of two) between the values of the effective diffusion coefficient for ¹²⁵I and ³⁶Cl. Whatever the ionic strength of the synthetic solution used, ¹²⁵I exhibited D_e values two times lower than those of ³⁶Cl. A detailed explanation for this difference cannot be given at present even if a hypothesis based on ion-pairing or on steric-exclusion cannot be excluded. This makes questionable the assumption usually made for quantifying ¹²⁵I sorption and postulating that ³⁶Cl and ¹²⁵I would diffuse in the same porosity. In other terms, at Tournemire, ¹²⁵I sorption could be more pronounced than previously indicated.     

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.     

The oxidation state of sulfur in synthetic and natural glasses determined by X-ray absorption spectroscopy

Sulfur K-edge X-ray absorption near edge structure (XANES) spectra were recorded for experimental glasses of various compositions prepared at different oxygen fugacities (fO₂) in one-atmosphere gas-mixing experiments at 1400 °C. This sample preparation method only results in measurable S concentrations under either relatively reduced (log fO₂ < −9) or oxidised (log fO₂ > −2) conditions. The XANES spectra of the reduced samples are characterised by an absorption edge crest at 2476.4 eV, typical of S²⁻. In addition, spectra of Fe-bearing compositions exhibit a pronounced absorption edge shoulder. Spectra for all the Fe-free samples are essentially identical, as are the spectra for the Fe-bearing compositions, despite significant compositional variability within each group. The presence of a sulfide phase, such as might exsolve on cooling, can be inferred from a pre-edge feature at 2470.5 eV.The XANES spectra of the oxidised samples are characterised by an intense transition at 2482.1 eV, typical of the sulfate anion SO₄²⁻. Sulfite (SO₃²⁻) has negligible solubility in silicate melts at low pressures. The previous identification of sulfite species in natural glass samples is attributed to an artefact of the analysis (photoreduction of S⁶⁺). S⁴⁺ does, however, occur unambiguously with S⁶⁺ in Fe-free and Fe-poor compositions prepared in equilibrium with CaSO₄ at 4-16 kbar, and when buffered with Re/ReO₂ at 10 kbar. Solubility of S⁴⁺ thus requires partial pressures of SO₂ considerably in excess of 1 bar. A number of experiments were undertaken in an attempt to access intermediate fO₂s more applicable to terrestrial volcanism. Although these were largely unsuccessful, S²⁻ and S⁶⁺ were found to coexist in some samples that were not in equilibrium with the imposed fO₂.The XANES spectra of natural olivine-hosted melt inclusions and submarine glasses representative of basalts at, or close to, sulfide saturation show mainly dissolved S²⁻, but with minor sulfate, and additionally a peak at 2469.5 eV, which, although presumably due to immiscible sulfide, is 1 eV lower than that typical of FeS. These sulfate and sulfide-related peaks disappear with homogenisation of the inclusions by heating to 1200 °C followed by rapid quenching, suggesting that both these features are a result of cooling under natural conditions. The presence of small amounts of sulfate in otherwise reduced basaltic magmas may be explained by the electron exchange reaction: S²⁻ + 8Fe³⁺ = S⁶⁺ + 8Fe²⁺, which is expected to proceed strongly to the right with decreasing temperature. This reaction would explain why S²⁻ and S⁶⁺ are frequently found together despite the very limited fO₂ range over which they are thermodynamically predicted to coexist. The S XANES spectra of water-rich, highly oxidised, basaltic inclusions hosted in olivine from Etna and Stromboli confirm that nearly all S is dissolved as sulfate, explaining their relatively high S contents.     

Nitrogen solubility in basaltic melt. Part I. Effect of oxygen fugacity

The role of the oxygen fugacity on the incorporation of nitrogen in basaltic magmas has been investigated using one atmosphere high temperature equilibration of tholeiitic-like compositions under controlled nitrogen and oxygen partial pressures in the [C-N-O] system. Nitrogen was extracted with a CO₂ laser under high vacuum and analyzed by static mass spectrometry. Over a redox range of 18 oxygen fugacity log units, this study shows that the incorporation of nitrogen in silicate melts follows two different behaviors. For log fO₂ values between −0.7 and −10.7 (the latter corresponding to IW − 1.3), nitrogen dissolves as a N₂ molecule into cavities of the silicate network (physical solubility). Nitrogen presents a constant solubility (Henry’s) coefficient of 2.21 ± 0.53 × 10⁻⁹ mol g⁻¹ atm⁻¹ at 1425°C, identical within uncertainties to the solubility of argon. Further decrease in the oxygen fugacity (log fO₂ between −10.7 and −18 corresponding to the range from IW − 1.3 to IW − 8.3) results in a drastic increase of the solubility of nitrogen by up to 5 orders of magnitude as nitrogen becomes chemically bounded with atoms of the silicate melt network (chemical solubility). The present results strongly suggest that under reducing conditions nitrogen dissolves in silicate melts as N³⁻ species rather than as CN⁻ cyanide radicals. The nitrogen content of a tholeiitic magma equilibrated with N₂ is computed from thermochemical processing of our data set as

     

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.     

Solubility and solution mechanism of H2O in alkali silicate melts and glasses at high pressure and temperature

The solubility behavior of H₂O in melts in the system Na₂O-SiO₂-H₂O was determined by locating the univariant phase boundary, melt = melt + vapor in the 0.8-2 GPa and 1000°-1300°C pressure and temperature range, respectively. The NBO/Si-range of the melts (0.25-1) was chosen to cover that of most natural magmatic liquids. The H₂O solubility in melts in the system Na₂O-SiO₂-H₂O (X_H2O) ranges between 18 and 45 mol% (O = 1) with (∂X_H2O/∂P)_T∼14-18 mol% H₂O/GPa. The (∂X_H2O/∂P)_T is negatively correlated with NBO/Si (= Na/Si) of the melt. The (∂X_H2O/∂T)_P is in the −0.03 to +0.05 mol% H₂O/°C range, and is negatively correlated with NBO/Si. The [∂X_H2O/∂(NBO/Si)]_P,T is in the −3 to −8 mol% H₂O/(NBO/Si) range. Melts with NBO/Si similar to basaltic liquids (∼0.6-∼1.0) show (∂X_H2O/∂T)_P<0, whereas more polymerized melts exhibit (∂X_H2O/∂T)_P>0. Complete miscibility between hydrous melt and aqueous fluid occurs in the 0.8-2 GPa pressure range for melts with NBO/Si ≤0.5 at T >1100°C. Miscibility occurs at lower pressure the more polymerized the melt.     

Thermodynamic properties of chlorite CCa-2. Heat capacities, heat contents and entropies

The heat capacities of the international reference clay mineral chlorite CCa-2 from Flagstaff Hill, California, were measured by low temperature adiabatic calorimetry and differential scanning calorimetry, from 5 to 520 K (at 1 bar). The studied chlorite is a Fe-bearing trioctahedral chlorite with an intermediary composition between ideal clinochlore (Si₃Al)(Mg₅Al)O₁₀(OH)₈ and chamosite (Si₃Al)(Fe₅Al)O₁₀(OH)₈. Only few TiO₂ impurities were detected in the natural chlorite sample CCa-2. Its structural formula, obtained after subtracting the remaining TiO₂ impurities, is (Si_2.633Al_1.367)(Al_1.116Mg_2.952Mn_0.012Ca_0.011)O₁₀(OH)₈. From the heat capacity results, the entropy, standard entropy of formation and heat content of the chlorite were deduced. At 298.15 K, the heat capacity of the chlorite is 547.02 (±0.27) J mol⁻¹ K⁻¹ and the molar entropy is 469.4 (±2.9) J mol⁻¹ K⁻¹. The standard molar entropy of formation of the clay mineral from the elements is −2169.4 (±4.0) J mol⁻¹ K⁻¹.