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

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

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.     

Water diffusion in dacitic melt

H₂O diffusion in dacitic melt was investigated at 0.48-0.95 GPa and 786-893 K in a piston-cylinder apparatus. The diffusion couple design was used, in which a nominally dry dacitic glass makes one half and is juxtaposed with a hydrous dacitic glass containing up to ∼8 wt.% total water (H₂O_t). H₂O concentration profiles were measured on quenched glasses with infrared microspectroscopy. The H₂O diffusivity in dacite increases rapidly with water content under experimental conditions, similar to previous measurements at the same temperature but at pressure <0.15 GPa. However, compared with the low-pressure data, H₂O diffusion at high pressure is systematically slower. H₂O diffusion profiles in dacite can be modeled by assuming molecular H₂O (H₂O_m) is the diffusing species. Total H₂O diffusivity D_H2Ot within 786-1798 K, 0-1 GPa, and 0-8 wt.% H₂O_t can be expressed as: where D_H2Ot is in m²/s, T is temperature in K, P is pressure in GPa, K = exp(1.49 − 2634/T) is the equilibrium constant of speciation reaction (H₂O_m+O?2OH) in the melt, X = C/18.015/[C/18.015 + (100 − C)/33.82], C is wt.% of H₂O_t, and 18.015 and 33.82 g/mol correspond to the molar masses of H₂O and anhydrous dacite on a single oxygen basis. Compared to H₂O diffusion in rhyolite, diffusivity in dacite is lower at intermediate temperatures but higher at superliquidus temperatures. This general H₂O diffusivity expression can be applied to a broad range of geological conditions, including both magma chamber processes and volcanic eruption dynamics from conduit to the surface.     

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.     

Trace/minor element:calcium ratios in cultured benthic foraminifera. Part II: Ontogenetic variation

Ontogenetic (developmental stage) measurements of Mg/Ca and Sr/Ca were made on the benthic foraminifer Bulimina aculeata, which were cultured under controlled physicochemical conditions of temperature, pH, alkalinity, salinity, and trace- and minor-element concentrations. We utilized two methods of ontogenetic sampling—whole specimens progressively increasing in length and laser microdissection of a single specimen with subsequent analysis of dissected portions. A novel high-resolution laser-microdissection (HRLM) method allowed for precise (10 μm) cuts of the foraminiferal tests (shells) along the geometrically complex sutures distinguishing individual chambers. This new microdissection method limited sample loss and cross-contamination between foraminiferal chambers. Little or no variation in D_Sr was observed at different foraminiferal developmental stages. Conversely, D_Mg was enriched during a mid-developmental stage of whole-specimen samples (150-225 μm D_Mg = 1.6 × 10⁻³) compared to earlier and later stages (<150 μm, >225 μm D_Mg = 8.3 × 10⁻⁴). Further analysis of HRLM ontogenetic samples showed a larger, age-dependent D_Mg signature variation. This increase in shell Mg/Ca may contribute substantially to the measured inter-individual variability in Mg/Ca temperature prediction for cultured B. aculeata. Due to relatively large Mg/Ca inter- and intra-individual variability, measuring similar-size foraminiferal samples may improve the precision of paleotemperature prediction. Additionally, partial dissolution of the highest ontogenetically Mg-enriched calcite (D_Mg = 1.3 × 10⁻²-1.6 × 10⁻²) may occur in undersaturated bottom-water environments or during reductive cleaning procedures. Thus, the calcite phases remaining after partial dissolution by either natural or laboratory cleaning processes may not accurately represent the calcification environment.     

Confined pyrolysis of Tertiary lacustrine source rocks in the Western Qaidam Basin,Northwest China: Implications for generative potential and oil maturity evaluation

The parameter S₁ + S₂ (genetic potential) of Rock-Eval analysis is widely used as an evaluation of the genetic potential for the source rocks. Oligocene–Miocene saline lacustrine source rocks in the western Qaidam basin have low total organic C contents (TOC), most around 0.5% with a few exceptions >1.0%. Mineral matrix effects are substantial for source rocks with low TOC, resulting in relatively low S₁ and S₂ peaks. Based on the results of confined pyrolyses (sealed Au capsules) on 6 Oligocene–Miocene source rocks from the western Qaidam basin, with TOC ranging between 0.48% and 2.22%, the relationship between the S₁ + S₂ parameter and the maximum amount of extracted bitumen or saturated and aromatic hydrocarbons (SA) after the confined pyrolysis has been established as follows: bitumen (mg/g rock) = 1.4924 × (S₁ + S₂) + 0.3201 (r = 0.987), or SA (saturates + aromatics) (mg/g rock) = 0.7083 × (S₁ + S₂) + 0.4045 (r = 0.992). Based on these formulas, the amounts of hydrocarbons generated from source rocks can be reasonably estimated. The typical crude oils with low biomarker maturities in this region appear substantially different to the pyrolysates of these six rocks at 180–300 °C but comparable to the pyrolysates at 320 °C and higher temperatures based on molecular parameters. This result, in combination with the physical and gross compositions of the crude oils, suggests that the majority of these crude oils were generated from the source rocks during the main oil-generative stage, possibly at a maturity higher than R_o 0.74%.     

Seawater temperature proxies based on DSr, DMg, and DU from culture experiments using the branching coral Porites cylindrica

In order to investigate the incorporation of Sr, Mg, and U into coral skeletons and its temperature dependency, we performed a culture experiment in which specimens of the branching coral (Porites cylindrica) were grown for 1 month at three seawater temperatures (22, 26, and 30 °C). The results of this study showed that the linear extension rate of P. cylindrica has little effect on the skeletal Sr/Ca, Mg/Ca, and U/Ca ratios. The following temperature equations were derived: Sr/Ca (mmol/mol) = 10.214(±0.229) − 0.0642(±0.00897) × T (°C) (r² = 0.59, p < 0.05); Mg/Ca (mmol/mol) = 1.973(±0.302) + 0.1002(±0.0118) × T (°C) (r² = 0.67, p < 0.05); and U/Ca (μmol/mol) = 1.488(±0.0484) − 0.0212(±0.00189) × T (°C) (r² = 0.78, p < 0.05). We calculated the distribution coefficient (D) of Sr, Mg, and U relative to seawater temperature and compared the results with previous data from massive Porites corals. The seawater temperature proxies based on D calibrations of P. cylindrica established in this study are generally similar to those for massive Porites corals, despite a difference in the slope of D_U calibration. The calibration sensitivity of D_Sr, D_Mg, and D_U to seawater temperature change during the experiment was 0.64%/°C, 1.93%/°C, and 1.97%/°C, respectively. These results suggest that the skeletal Sr/Ca ratio (and possibly the Mg/Ca and/or U/Ca ratio) of the branching coral P. cylindrica can be used as a potential paleothermometer.     

Solubility of cyclooctasulfur in pure water and sea water at different temperatures

The solubility of cyclooctasulfur in water and sea water at various temperatures in the range between 4 and 80 °C was determined. Cyclooctasulfur in equilibrium with rhombic sulfur reacted with hot acidic aqueous potassium cyanide to form thiocyanate anion which was measured by anion chromatography. Sulfur solubility in pure water was found to increase with temperature by more than 78 times: from 6.1 nM S₈ at 4 °C to 478 nM S₈ at 80 °C. The following thermodynamic values for solubilisation of S₈ in water were calculated from the experimental data: K° = 3.01 ± 1.04 × 10⁻⁸, ΔG_r° = 42.93 ± 0.73 kJ mol⁻¹, ΔH_r° = 47.4 ± 3.6 kJmol⁻¹, ΔS_r° = 15.0 ± 11.7 J mol⁻¹ K⁻¹). Solubility of cyclooctasulfur in sea water was found to be 61 ± 13% of the solubility in pure water regardless of the temperature.     

Silicon control of strontium and cesium partitioning in hydroxide-weathered sediments

Cation partitioning and speciation in an aqueous soil suspension may depend on the coupling of reaction time, sorbate amount and mineral weathering reactions. These factors were varied in sediment suspension experiments to identify geochemical processes that affect migration of Sr²⁺ and Cs⁺ introduced to the subsurface by caustic high level radioactive waste (HLRW). Three glacio-fluvial and lacustrine sediments from the Hanford Site (WA, USA) were subjected to hyperalkaline (pH > 13), Na-Al-NO₃-OH solution conditions within a gradient field of (i) sorptive concentration (10⁻⁵-10⁻³ m) and (ii) reaction time (0-365 d). Strontium uptake (q_Sr) exceeded that of cesium at nearly all reaction times. Sorbent affinity for both Cs⁺ and Sr²⁺ increased with clay plus silt content at early times, but a prolonged slow uptake process was observed over the course of sediment weathering that erased the texture effect for Sr²⁺; all sediments showed similar mass normalized uptake after several months of reaction time. Strontium became progressively recalcitrant to desorption after 92 d, with accumulation and aging of neoformed aluminosilicates. Formation of Cs⁺ and Sr²⁺-containing cancrinite and sodalite was observed after 183 d by SEM and synchrotron μ-XRF and μ-XRD. EXAFS data for q_Sr ≈ 40 mmol kg⁻¹ showed incorporation of Sr²⁺ into both feldspathoid and SrCO₃(s) coordination environments after one year. Adsorption was predominant at early times and low sorbate amount, whereas precipitation, controlled largely by sediment Si release, became increasingly important at longer times and higher sorbate amount. Kinetics of contaminant desorption at pH 8 from one year-weathered sediments showed significant dependence on background cation (Ca²⁺ versus K⁺) composition. Results of this study indicate that co-precipitation and ion exchange in neoformed aluminosilicates may be an important mechanism controlling Sr²⁺ and Cs⁺ mobility in siliceous sediments impacted by hyperalkaline HLRW.     

Solubility of plutonium(VI) carbonate in saline solutions

Among the plutonium oxidation states found to form in the environment, mobile plutonium(VI) can exist under oxidizing conditions and in waters with high chloride content due to radiolysis effects. We are investigating the solubility and speciation of plutonium(VI) carbonate under conditions relevant to natural waters and brines such as those found near some geologic radioactive waste repositories. The solid Pu(VI) phase PuO₂CO₃(s) was prepared and its solubility was measured in NaCl and NaClO₄ solutions in a CO₂ atmosphere as a function of pH and ionic strength (0.1-5.6 m). The concentration of soluble plutonium in solution was calculated from spectroscopic data and liquid scintillation counting. Spectroscopic measurements also revealed the plutonium oxidation state. The apparent solubility product of PuO₂CO₃(s) was determined at selected electrolyte concentrations to be, log K_s,0 = −13.95 ± 0.07 (0.1 m NaCl), log K_s,0 = −14.07 ± 0.13 (5.6 m NaCl), and log K_s,0 = −15.26 ± 0.11 (5.6 m NaClO₄). Specific ion interaction theory was used to calculate the solubility product at zero ionic strength, .     

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.     

Trace element partitioning between apatite and silicate melts

We present new experimental apatite/melt trace element partition coefficients for a large number of trace elements (Cs, Rb, Ba, La, Ce, Pr, Sm, Gd, Lu, Y, Sr, Zr, Hf, Nb, Ta, U, Pb, and Th). The experiments were conducted at pressures of 1.0 GPa and temperatures of 1250 °C. The rare earth elements (La, Ce, Pr, Sm, Gd, and Lu), Y, and Sr are compatible in apatite, whereas the larger lithophile elements (Cs, Rb, and Ba) are strongly incompatible. Other trace elements such as U, Th, and Pb have partition coefficients close to unity. In all experiments we found D_Hf > D_Zr, D_Ta ≈ D_Nb, and D_Ba > D_Rb > D_Cs. The experiments reveal a strong influence of melt composition on REE partition coefficients. With increasing polymerisation of the melt, apatite/melt partition coefficients for the rare earth elements increase for about an order of magnitude. We also present some results in fluorine-rich and water-rich systems, respectively, but no significant influence of either H₂O or F on the partitioning was found. Furthermore, we also present experimentally determined partition coefficients in close-to natural compositions which should be directly applicable to magmatic processes.     

Determination of fluid/melt partition coefficients by LA-ICPMS analysis of co-existing fluid and silicate melt inclusions: Controls on element partitioning

Analyses of co-existing silicate melt and fluid inclusions, entrapped in quartz crystals in volatile saturated magmatic systems, allowed direct quantitative determination of fluid/melt partition coefficients. Investigations of various granitic systems (peralkaline to peraluminous in composition, log fO₂ = NNO−1.7 to NNO+4.5) exsolving fluids with various chlorinities (1-14 mol/kg) allowed us to assess the effect of these variables on the fluid/melt partition coefficients (D). Partition coefficients for Pb, Zn, Ag and Fe show a nearly linear increase with the chlorinity of these fluid (D_Pb ∼ 6 ∗ m_Cl, D_Zn ∼ 8 ∗ m_Cl, D_Ag ∼ 4 ∗ m_Cl, D_Fe ∼ 1.4 ∗ m_Cl, where m_Cl is the molinity of Cl). This suggests that these metals are dissolved primarily as Cl-complexes and neither oxygen fugacity nor the composition of the melt affects significantly their fluid/melt partitioning. By contrast, partition coefficients for Mo, B, As, Sb and Bi are highest in low salinity (1-2 mol/kg Cl) fluids with maximum values of D_Mo ∼ 20, D_B ∼ 15, D_As ∼ 13, D_Sb ∼ 8, D_Bi ∼ 15 indicating dissolution as non-chloride (e.g., hydroxy) complexes. Fluid/melt partition coefficients of copper are highly variable, but highest between vapor like fluids and silicate melt (D_Cu ? 2700), indicating an important role for ligands other than Cl. Partition coefficients for W generally increase with increasing chlorinity, but are exceptionally low in some of the studied brines which may indicate an effect of other parameters. Fluid/melt partition coefficients of Sn show a high variability but likely increase with the chlorinity of the fluid (D_Sn = 0.3-42, D_W = 0.8-60), and decrease with decreasing oxygen fugacity or melt peraluminosity.     

Zinc complexation in aqueous sulfide solutions: Determination of the stoichiometry and stability of complexes via ZnS(cr) solubility measurements at 100 °C and 150 bars

The solubility of ZnS_(cr) was measured at 100 °C, 150 bars in sulfide solutions as a function of sulfur concentration (m(S_total) = 0.02-0.15) and acidity (pH_t = 2-11). The experiments were conducted using a Ti flow-through hydrothermal reactor enabling the sampling of large volumes of solutions at experimental conditions, with the subsequent concentration and determination of trace quantities of Zn. Prior to the experiments, a long-term in situ conditioning of the solid phase was performed in order to attain the reproducible Zn concentrations (i.e. solubilities). The ZnS_(cr) solubility product was monitored in the course of the experiment. The following species were found to account for Zn speciation in solution: Zn²⁺ (pH_t < 3), (pH_t 3-4.5), (pH_t 5-8), and ZnS(HS)⁻ (pH_t > 8) (pH_t predominance regions are given for m(S_total) = 0.1). Solubility data collected in this study at pH_t > 3 were combined with the ZnS_(cr) solubility product determined at lower pH to yield the following equilibrium constants (t = 100 °C, P = 150 bars):

     

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.     

Raman spectroscopic measurements of CO2 density: Experimental calibration with high-pressure optical cell (HPOC) and fused silica capillary capsule (FSCC) with application to fluid inclusion observations

Raman spectroscopy is a powerful method for the determination of CO₂ densities in fluid inclusions, especially for those with small size and/or low fluid density. The relationship between CO₂ Fermi diad split (Δ, cm⁻¹) and CO₂ density (ρ, g/cm³) has been documented by several previous studies. However, significant discrepancies exist among these studies mainly because of inconsistent calibration procedures and lack of measurements for CO₂ fluids having densities between 0.21 and 0.75 g/cm³, where liquid and vapor phases coexist near room temperature.In this study, a high-pressure optical cell and fused silica capillary capsules were used to prepare pure CO₂ samples with densities between 0.0472 and 1.0060 g/cm³. The measured CO₂ Fermi diad splits were calibrated with two well established Raman bands of benzonitrile at 1192.6 and 1598.9 cm⁻¹. The relationship between the CO₂ Fermi diad split and density can be represented by: ρ = 47513.64243 − 1374.824414 × Δ + 13.25586152 × Δ² − 0.04258891551 × Δ³ (r² = 0.99835, σ = 0.0253 g/cm³), and this relationship was tested by synthetic fluid inclusions and natural CO₂-rich fluid inclusions. The effects of temperature and the presence of H₂O and CH₄ on this relationship were also examined.     

Water speciation and diffusion in haploandesitic melts at 743-873 K and 100 MPa

Water is an important volatile component in andesitic eruptions and deep-seated andesitic magma chambers. We report an investigation of H₂O speciation and diffusion by dehydrating haploandesitic melts containing ?2.5 wt.% water at 743-873 K and 100 MPa in cold-seal pressure vessels. FTIR microspectroscopy was utilized to measure species [molecular H₂O (H₂O_m) and hydroxyl group (OH)] and total H₂O (H₂O_t) concentration profiles on the quenched glasses from the dehydration experiments. The equilibrium constant of the H₂O speciation reaction H₂O_m+O?2OH, K = (X_OH)²/(X_H2OmX_O) where X means mole fraction on a single oxygen basis, in this Fe-free andesite varies with temperature as ln K = 1.547-2453/T where T is in K. Comparison with previous speciation data on rhyolitic and dacitic melts indicates that, for a given water concentration, Fe-free andesitic melt contains more hydroxyl groups. Water diffusivity at the experimental conditions increases rapidly with H₂O concentration, contrary to previous H₂O diffusion data in an andesitic melt at 1608-1848 K. The diffusion profiles are consistent with the model that molecular H₂O is the diffusion species. Based on the above speciation model, H₂O_m and H₂O_t diffusivity (in m²/s) in haploandesite at 743-873 K, 100 MPa, and H₂O_t ? 2.5 wt.% can be formulated as

     

Retention of biosignatures and mass-independent fractionations in pyrite: Self-diffusion of sulfur

Self-diffusion of sulfur in pyrite (FeS₂) was characterized over the temperature range ∼500-725 °C (∼1 bar pressure) by immersing natural specimens in a bath of molten elemental ³⁴S and characterizing the resulting diffusive-exchange profiles by Rutherford backscattering spectroscopy (RBS). The temperature dependence of the sulfur diffusivity (D_S) conforms to D_S = D_o exp(−E_a/RT), where the pre-exponential constant (D_o) and the activation energy (E_a) are constrained as follows:

     

Internally consistent thermodynamic data for magnesium sulfate hydrates

Experimental studies on the stability of several Mg-sulfate hydrates including epsomite (MgSO₄·7H₂O), hexahydrite (MgSO₄·6H₂O), starkeyite (MgSO₄·4H₂O), and kieserite (MgSO₄·H₂O) as a function of temperature and relative humidity are in poor agreement with calculations based on thermodynamic properties of these substances taken from the literature. Therefore, we synthesized four different MgSO₄ hydrates and measured their enthalpies of formation by solution calorimetry at T = 298.15 K. The resulting enthalpies of formation from the elements are:

•

Δ_fH⁰₂₉₈ (epsomite) = −3387.7 ± 1.3 kJmol⁻¹

•

Δ_fH⁰₂₉₈ (hexahydrite) = −3088.1 ± 1.1 kJmol⁻¹

•

Δ_fH⁰₂₉₈ (sanderite, MgSO₄·2H₂O) = −1894.9 ± 1.3 kJmol⁻¹

•

Δ_fH⁰₂₉₈ (kieserite) = −1612.4 ± 1.3 kJmol⁻¹

Using mathematical programming (MAP) techniques, standard thermodynamic values consistent both with our calorimetric data and previously published humidity brackets could be derived:

•

Epsomite: Δ_fH⁰₂₉₈ = −3388.7 kJmol⁻¹, S⁰₂₉₈ = 371.3 Jmol⁻¹ K⁻¹, Δ_fG⁰₂₉₈ = −2871.0 kJmol⁻¹

•

Hexahydrite: Δ_fH⁰₂₉₈ = −3087.3 kJmol⁻¹, S⁰₂₉₈ = 348.5 Jmol⁻¹ K⁻¹, Δ_fG⁰₂₉₈ = −2632.3 kJmol⁻¹

•

Starkeyite: Δ_fH⁰₂₉₈ = −2496.1 kJmol⁻¹, S⁰₂₉₈ = 259.9 Jmol⁻¹ K⁻¹, Δ_fG⁰₂₉₈ = −2153.8 kJmol⁻¹

•

Kieserite: Δ_fH⁰₂₉₈ = −1611.5 kJmol⁻¹, S⁰₂₉₈ = 126.0 Jmol⁻¹ K⁻¹, Δ_fG⁰₂₉₈ = −1437.9 kJmol⁻¹

Additionally, heat capacity measurements and standard entropy determinations of several magnesium sulfate hydrate minerals from the literature are analyzed and judged against estimates obtained from a linear combination of the heat capacities of MgSO₄ and hexagonal ice. The results of the MAP analysis are compared to these estimates to conclude that heat capacity and entropy correlate well with the number of waters of hydration. However, even the good correlation is not good enough to capture the fine variations in these properties. Consequently, their experimental measurement is inevitable if reliable thermodynamic data are sought. Our MAP thermodynamic data show that epsomite, hexahydrite, and kieserite have stability fields in the T-%RH space. Starkeyite is metastable. Although no MAP data could have been derived for pentahydrite (MgSO₄·5H₂O) and sanderite, their transient existence suggest that both of them are metastable as well.