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

Thermochemistry of strontium incorporation in aragonite from atomistic simulations

We have investigated the thermodynamics of mixing between aragonite (orthorhombic CaCO₃) and strontianite (SrCO₃). In agreement with experiment, our simulations predict that there is a miscibility gap between the two solids at ambient conditions. All Sr_xCa_1−xCO₃ solids with compositions 0.12 < x < 0.87 are metastable with respect to separation into a Ca-rich and a Sr-rich phase. The concentration of Sr in coral aragonites (x ∼ 0.01) lies in the miscibility region of the phase diagram, and therefore formation of separated Sr-rich phases in coral aragonites is not thermodynamically favorable. The miscibility gap disappears at around 380 K. The enthalpy of mixing, which is positive and nearly symmetric with respect to x = 0.5, is the dominant contribution to the excess free energy, while the vibrational and configurational entropic contributions are small and of opposite sign. We provide a detailed comparison of our simulation results with available experimental data.     

New collectors for sphalerite flotation

A series of N-arylhydroxamic acids (31) were synthesized and tested as collectors to float sphalerite from a Canadian copper–zinc ore. The compounds were classified into four types namely, N-aryl-C-alkyl, N-aryl-C-aryl, N-aryl-C-aralkyl, and dihydroxamic acids based on the type of substitution. Dihydroxamic acids were found to be poor mineral collectors while the efficiency of sphalerite flotation increased in the order N-aryl-C-aryl < N-aryl-C-alkyl < N-aryl-C-aralkyl. Sphalerite was floated without activation by copper sulfate, and the best sphalerite recovery of about 80% (grade 32%) was achieved with N-hydrocinnamoyl-N-phenylhydroxylamine (HCNPHA) 67 g/t collector dosage. However, pyrite also floated along with sphalerite and this appeared as a major disadvantage to be addressed.     

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.     

A Raman spectroscopic investigation of Fe-rich sphalerite: effect of fe-substitution

A Raman spectroscopic study of Fe-rich sphalerite (Zn_{1 − x} Fe _x S) has been carried out for six samples with 0.10 ≤ x ≤ 0.24. Both the intensities and frequencies of the TO and LO modes of sphalerite are approximately independent of Fe concentration. However, the substitution of Zn by Fe results in five additional bands with frequencies between the TO (271 cm⁻¹) and LO (350 cm⁻¹) modes. Three of these bands are attributed to resonance modes (i.e. Y ₁, Y ₂ and Y ₃ modes). The fourth band (B mode) is assigned to a breathing mode of the nearest-neighbor sulfur atoms around the Fe atoms. The band at 337 cm⁻¹ is attributed to the presence of Fe³⁺. The excellent correlations between the normalized intensities of these five different modes and x _Fe show that these modes depend on Fe-content. Another extra mode at 287 cm⁻¹ is assigned to the presence of Cd in sphalerite.     

Experimental determination of the Raman CH4 symmetric stretching (ν1) band position from 1-650 bar and 0.3-22 °C: Application to fluid inclusion studies

The position of the Raman methane (CH₄) symmetric stretching band (ν₁) over the range 1-650 bar and 0.3-22 °C has been determined using a high-pressure optical cell mounted on a Raman microprobe. Two neon emission lines that closely bracket the CH₄ band were collected simultaneously with each CH₄ spectrum. The peak position was determined after least squares fitting using a summed Gaussian-Lorentzian method, resulting in a precision of ≈±0.02 cm⁻¹ in peak position determination. The CH₄ν₁ band position shifts to lower wave number with increasing pressure. At a given pressure, the band shifts to lower wave number with decreasing temperature, and the magnitude of the temperature shift increases with increasing pressure. The relationship between the Raman CH₄ν₁ band position and temperature and pressure determined here may be used to estimate the internal pressure in natural or synthetic CH₄-bearing fluid inclusions. This information, in turn, may be used to determine the density of pure CH₄ fluid inclusions and the salinity of CH₄-bearing aqueous inclusions.     

An experimental study on the preparation of tochilinite-originated intercalation compounds comprised of Fe1−xS host layers and various kinds of guest layers

Tochilinite represents a mineral group of ordered mixed-layer structures containing alternating Fe_1−xS layers with mackinawite-like structure and metal hydroxide layers with Mg(OH)₂-like structure. In this article, we report the preparation of a series of tochilinite-originated (or Fe_1−xS-based) intercalation compounds (ICs). According to their preparation procedures, these ICs can be divided into four kinds. The first kind of IC was sodium tochilinite (Na-tochilinite), which was prepared by the hydrothermal reaction of metallic Fe particles with concentrated Na₂S·9H₂O aqueous solutions. The hydroxide layer of the Na-tochilinite was a mixed hydroxide of Na⁺ ions along with a certain amount of Fe²⁺ ions. When the hydroxide layer of the Na-tochilinite completely dissolved in aqueous solutions, a Fe-deficient mackinawite-like phase Fe_1−xS was obtained, which was probably an electron-deficient p-type conductor. The second kind of ICs was prepared by ‘low-temperature direct intercalation in aqueous solutions, using Na-tochilinite as a parental precursor. When the Na-tochilinite was ultrasonicated in aqueous solutions containing Lewis basic complexing agents (like NH₃, N₂H₄, 2,2′-bipyridine (bipy), and 1,10-phenanthroline (phen)), the Na⁺ ions of the Na-tochilinite were removed and the Lewis basic complexing agents entered the hydroxide layer of the Na-tochilinite and became coordinated with the Fe²⁺ ions, and the second kind of ICs was thus produced. The second kind of ICs includes NH₃ IC, N₂H₄ IC, N₂H₄-NH₃ IC, [Fe(bipy)₃]²⁺-containing IC and [Fe(phen)₃]²⁺-containing IC. The third kind of ICs, which includes NH₃ IC, N₂H₄-NH₃ IC and N₂H₄-LiOH (NaOH) IC, was prepared by the hydrothermal reaction of metallic Fe particles with (NH₄)₂S aqueous solution, S (elemental) + N₂H₄·H₂O aqueous solution, and S + N₂H₄·H₂O + LiOH (NaOH) aqueous solution, respectively. The third kind of ICs has a close relationship with the second kind of ICs both in composition and structure. The fourth kind of ICs was prepared by the oxidation and reduction of some of the N₂H₄-containing ICs mentioned above, which include N₂H₂ (diazene or diimide) IC, N₂ (dinitrogen) IC and NH₃ IC. The N₂H₂ IC was prepared by mild air oxidation of the N₂H₄-LiOH IC. The N₂ IC was prepared by strong air oxidation of the N₂H₄-LiOH IC, however, we have not been able to separate the pure phase N₂ IC. Hydrothermal reduction of the N₂H₄ IC made by the direct intercalation method in strong reducing environment by H₂S + Fe (metal) led to the production of the NH₃ IC of the fourth kind of ICs. The NH₃ ICs prepared by the three methods had similar compositions and structures. As almost all the ICs reported in this paper were extremely sensitive both to air and to the electron beam, they were mainly characterized by XRD.The properties and interrelationships (or mutual transformations) of the Fe_1−xS-based ICs revealed novel chemistry occurring in the sub-nanoscopic space between the micrometer- to nanometer-sized electron-deficient Fe_1−xS layers. An important finding of this novel chemistry was that the Fe_1−xS-based ICs tended to oxidize or reduce the intercalated species when the redox state of their environments varied. The results of our experiments potentially have many cosmochemical implications. The most important implication is that our experimental results, along with previous studies, strongly suggested that some of the ammonium salts, ammonia and carbonates existing in the matrix of the CM carbonaceous chondrites may have been formed by abiotic reactions employing molecular nitrogen as the nitrogen source and carbon monoxide as the carbon source and iron sulfide and/or iron hydroxide as catalysts.     

Dissolution kinetics of diopside as a function of solution saturation state: Macroscopic measurements and implications for modeling of geological storage of CO2

Measurements of the dissolution rate of diopside (r) were carried out as a function of the Gibbs free energy of the dissolution reaction (ΔG_r) in a continuously stirred flow-through reactor at 90 °C and pH_90 °C = 5.05. The overall relation between r and ΔG_r was determined over a free energy range of −130.9 < ΔG_r < −47.0 kJ mo1⁻¹. The data define a highly non-linear, sigmoidal relation between r and ΔG_r. At far-from-equilibrium conditions (ΔG_r ? −76.2 kJ mo1⁻¹), a rate plateau is observed. In this free energy range, the rates of dissolution are constant, independent of [Ca], [Mg] and [Si] concentrations, and independent of ΔG_r. A sharp decrease of the dissolution rate (∼1 order of magnitude) occurs in the transition ΔG_r region defined by −76.2 < ΔG_r ? −61.5 kJ mo1⁻¹. Dissolution closer to equilibrium (ΔG_r > −61.5 kJ mo1⁻¹) is characterised by a much weaker inverse dependence of the rates on ΔG_r. Modeling the experimental r-ΔG_r data with a simple classical transition state theory (TST) law as implemented in most available geochemical codes is found inappropriate. An evaluation of the consequences of the use of geochemical codes where the r-ΔG_r relation is based on basic TST was carried out and applied to carbonation reactions of diopside, which, among other reactions with Ca- and Mg-bearing minerals, are considered as a promising process for the solid state sequestration of CO₂ over long time spans. In order to take into account the actual experimental r-ΔG_r relation in the geochemical code that we used, a new module has been developed. It reveals a dramatic overestimation of the carbonation rate when using a TST-based geochemical code. This points out that simulations of water-rock-CO₂ interactions performed with classical geochemical codes should be evaluated with great caution.     

Activity-composition relations in the system CaCO3-MgCO3 predicted from static structure energy calculations and Monte Carlo simulations

Thermodynamic mixing properties and subsolidus phase relations of the rhombohedral carbonate system, (1 − x) · CaCO₃ − x · MgCO₃, were modelled in the temperature range of 623-2023 K with static structure energy calculations based on well-parameterised empirical interatomic potentials. Relaxed static structure energies of a large set of randomly varied structures in a 4 × 4 × 1 supercell of calcite (a = 19.952 Å, c = 17.061 Å) were calculated with the General Utility Lattice Program (GULP). These energies were cluster expanded in a basis set of 12 pair-wise effective interactions. Temperature-dependent enthalpies of mixing were calculated by the Monte Carlo method. Free energies of mixing were obtained by thermodynamic integration of the Monte Carlo results. The calculated phase diagram is in good agreement with experimental phase boundaries.     

Spectroscopic analysis (FTIR, Raman) of water in mafic and intermediate glasses and glass inclusions

Micro-Raman spectroscopy, even though a very promising technique, is not still routinely applied to analyse H₂O in silicate glasses. The accuracy of Raman water determinations critically depends on the capability to predict and take into account both the matrix effects (bulk glass composition) and the analytical conditions on band intensities. On the other hand, micro-Fourier transform infrared spectroscopy is commonly used to measure the hydrous absorbing species (e.g., hydroxyl OH⁻ and molecular H₂O) in natural glasses, but requires critical assumptions for the study of crystal-hosted glasses. Here, we quantify for the first time the matrix effect of Raman external calibration procedures for the quantification of the total H₂O content (H₂O_T = OH⁻ + H₂O_m) in natural silicate glasses. The procedures are based on the calibration of either the absolute (external calibration) or scaled (parameterisation) intensity of the 3550 cm⁻¹ band. A total of 67 mafic (basanite, basalt) and intermediate (andesite) glasses hosted in olivines, having between 0.2 and 4.8 wt% of H₂O, was analysed. Our new dataset demonstrates, for given water content, the height (intensity) of Raman H₂O_T band depends on glass density, reflectance and water environment. Hence this matrix effect must be considered in the quantification of H₂O by Raman spectroscopy irrespective of the procedure, whereas the parameterisation mainly helps to predict and verify the self-consistency of the Raman results. In addition, to validate the capability of the micro-Raman to accurately determine the H₂O content of multicomponent aluminosilicate glasses, a subset of 23 glasses was analysed by both micro-Raman and micro-FTIR spectroscopy using the band at 3550 cm⁻¹. We provide new FTIR absorptivity coefficients (ε₃₅₅₀) for basalt (62.80 ± 0.8 L mol⁻¹ cm⁻¹) and basanite (43.96 ± 0.6 L mol⁻¹ cm⁻¹). These values, together with an exhaustive review of literature data, confirm the non-linear decline of the FTIR absorptivity coefficient (ε₃₅₅₀) as the glass depolymerisation increases. We demonstrate the good agreement between micro-FTIR and micro-Raman determination of H₂O in silicate glasses when the matrix effects are properly considered.     

Tissue N content and N natural abundance in epilithic mosses for indicating atmospheric N deposition in the Guiyang area,SW China

Tissue N contents and δ¹⁵N signatures in 175 epilithic mosses were investigated from urban to rural sites in Guiyang (SW China) to determine atmospheric N deposition. Moss N contents (0.85–2.97%) showed a significant decrease from the urban area (mean = 2.24 ± 0.32%, 0–5 km) to the rural area (mean = 1.27 ± 0.13%, 20–25 km), indicating that the level of N deposition decreased away from the urban environment, while slightly higher N contents re-occurred at sites beyond 30 km, suggesting higher N deposition in more remote rural areas. Moss δ¹⁵N ranged from −12.50‰ to −1.39‰ and showed a clear bimodal distribution (−12‰ to −6‰ and −5‰ to −2‰), suggesting that there are two main sources for N deposition in the Guiyang area. More negative δ¹⁵N (mean = −8.87 ± 1.65‰) of urban mosses mainly indicated NH₃ released from excretory wastes and sewage, while the less negative δ¹⁵N (from −3.83 ± 0.82‰ to −2.48 ± 0.95‰) of rural mosses were mainly influenced by agricultural NH₃. With more negative values in the urban area than in the rural area, the pattern of moss δ¹⁵N variation in Guiyang was found to be opposite to cities where N deposition is dominated by NO_x–N. Therefore, NH_x–N is the dominant N form deposited in the Guiyang area, which is supported by higher NH_x–N than NO_x–N in local atmospheric deposition. From the data showing that moss is responding to NH_x–N/NO_x–N in deposition it can be further demonstrated that the variation of moss δ¹⁵N from the Guiyang urban to rural area was more likely controlled by the ratio of urban-NH_x/agriculture-NH_x than the ratio of NH_x–N/NO_x–N. The results of this study have extended knowledge of atmospheric N sources in city areas, showing that urban sewage discharge could be important in cities co-generic to Guiyang.     

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

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

A unified equation for calculating methane vapor pressures in the CH4-H2O system with measured Raman shifts

A unified equation has been derived by using all available data for calculating methane vapor pressures with measured Raman shifts of C-H symmetric stretching band (υ₁) in the vapor phase of sample fluids near room temperature. This equation eliminates discrepancies among the existing data sets and can be applied at any Raman laboratory. Raman shifts of C-H symmetric stretching band of methane in the vapor phase of CH₄-H₂O mixtures prepared in a high-pressure optical cell were also measured at temperatures between room temperature and 200 °C, and pressures up to 37 MPa. The results show that the CH₄υ₁ band position shifts to higher wavenumber as temperature increases. We also demonstrated that this Raman band shift is a simple function of methane vapor density, and, therefore, when combined with equation of state of methane, methane vapor pressures in the sample fluids at elevated temperatures can be calculated from measured Raman peak positions. This method can be applied to determine the pressure of CH₄-bearing systems, such as methane-rich fluid inclusions from sedimentary basins or experimental fluids in hydrothermal diamond-anvil cell or other types of optical cell.     

Iron-hydroxide, iron-sulfate and hydrous-silica coatings in acid-mine tailings facilities: A comparative study of their trace-element composition

Surface alteration-layers often coat minerals in acid-mine drainage systems and the characterization of their chemical composition is required to understand the uptake or release of potentially toxic elements. Samples with micrometer-thick rock coatings were collected from bedrock in contact with three acidic tailings ponds and a small lake, all located within the Copper Cliff mine tailings disposal area in Sudbury, Ontario, Canada. Distribution and concentration of trace-metals in the rock coatings were characterized with Laser-Ablation Inductively-Coupled Plasma Mass Spectroscopy and Micro X-ray Fluorescence Spectroscopy. The rock coatings are composed of goethite, ferrihydrite, schwertmannite, jarosite and amorphous silica. The latter phase is a product of the non-stoichiometric weathering of the underlying siliceous rock. Layers within the coatings are distinguished on the basis of their atomic Fe:Si ratios: FeO_x coatings have Fe:Si > 4:1, Si-FeO_x coatings have Fe:Si = 4:1 to 1:1 and SiO_x coatings have Si > Fe. Iron-rich coatings (FeO_x) in contact with acidic tailings ponds (pH < 3.5) have lower trace-metal concentrations than their Si-rich counterparts, whereas FeO_x in contact with lake water at near neutral pH have similar trace-metal concentrations than Si-FeO_x and SiO_x, most likely the result of higher adsorption rates of metals at near neutral pH conditions. High trace-metal concentrations in Si-FeO_x and SiO_x are explained by the presence of jarosite-group minerals, which formed within Si-rich alteration layers through mixing of leached alkaline cations and trace elements from the underlying rock and Fe³⁺-sulfate solutions from the pond. Calculated enrichment factors for trace metals and metalloids in the coatings (relative to the pond) indicate that the mobility for Pb, As, Cr and Cu in the upper part of tailings ponds is commonly lower than the mobility for Zn, Mn, Co and Ni. The environmental significance of these findings is discussed in terms of the attenuation of trace metals in the coatings and the widespread occurrences of silica gels and jarosite-group minerals.     

The crystal and melt structure of spinel and alumina at high temperature: An in-situ XANES study at the Al and Mg K-edge

We present structural information obtained on spinel and alumina at high temperature (298-2400 K) using in-situ XANES at the Mg and Al K-edges. For spinel, ^[4](Al_x,Mg_1−x)^[6](Al_2−x,Mg_x)O_4, with increasing temperature, a substitution of Mg by Al and Al by Mg in their respective sites is observed. This substitution corresponds to an inversion of the Mg and Al sites. There is a significant change in the Al K-edge spectra between crystal and liquid, which can be attributed to a change of the ^[6]Al normally observed in corundum at room temperature, to a mixture of ^[6]Al-^[4]Al in the liquid state. This conclusion is in good agreement with previous ²⁷Al NMR experiments. Furthermore, both experiments at the Al and Mg K-edges are in good agreement with XANES calculation made using FDMNES code.     

Application of the Pitzer ion interaction model to isopiestic data for the Fe2(SO4)3-H2SO4-H2O system at 298.15 and 323.15 K

Recent isopiestic studies of the Fe₂(SO₄)₃-H₂SO₄-H₂O system at 298.15 K are represented with an extended version of Pitzer’s ion interaction model. The model represents osmotic coefficients for aqueous {(1 − y)Fe₂(SO₄)₃ + yH₂SO₄} mixtures from 0.45 to 3.0 m at 298.15 K and 0.0435 ? y ? 0.9370. In addition, a slightly less accurate representation of a more extended molality range to 5.47 m extends over the same y values, translating to a maximum ionic strength of 45 m. Recent isopiestic data for the system at 323.15 K are represented with the extended Pitzer model over a limited range in molality and solute fraction. These datasets are also represented with the usual “3-parameter” version of Pitzer’s model so that it may be incorporated in geochemical modeling software, but is a slightly less accurate representation of thermodynamic properties for this system. Comparisons made between our ion interaction model and available solubility data display partial agreement for rhomboclase and significant discrepancy for ferricopiapite. The comparisons highlight uncertainty remaining for solubility predictions in this system as well as the need for additional solubility measurements for Fe³⁺-bearing sulfate minerals. The resulting Pitzer ion interaction models provide an important step toward an accurate and comprehensive representation of thermodynamic properties in this geochemically important system.     

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.     

X-ray photoelectron spectroscopic studies of dolomite surfaces exposed to undersaturated and supersaturated aqueous solutions

Cleaved surfaces of dolomite were studied using ex-situ X-ray photoelectron spectroscopy (XPS) following exposure of the surfaces to various experimental conditions. Dolomite samples exposed to air, to a highly undersaturated solution (0.1 M NaCl, pH = 9), and to solution with a supersaturation (−Δμ/kT) of 5.5 (pH = 9) were investigated with semiquantitative methods of analysis to ascertain the degree of non-stoichiometry resulting at the dolomite surface from reactive conditions. It was found that the dolomite cleavage surface in undersaturated solution was not altered significantly from the stoichiometric surface termination. The composition of the cleaved surface after exposure to supersaturated solution, a surface known to have self-limiting growth characteristics under similar conditions, was found to be Ca²⁺ rich (Ca_xMg_2 − x(CO₃)₂, 1.7 > x > 1.3). The observations, while underscoring differences in hydration/dehydration kinetics of the two alkaline earth cations, suggest that achievement of equilibrium at dolomite-water interfaces may be subject to significant barriers from both undersaturated and supersaturated solutions.