Adsorption of Sr and Ba at the cleaved mica-water interface: Free energy profiles and interfacial structure

Monte Carlo simulations show that the adsorption position of the Sr²⁺ or Ba²⁺ ion on the cleaved muscovite surface is determined by the radius of the ion’s first hydration shell, hydrogen bonding of the first shell water molecules with the basal oxygens of muscovite as well as a requirement of minimization of the number of muscovite’s lattice cations in the ion’s first coordination shell. Accordingly, Sr²⁺ or Ba²⁺ adsorbs in ditrigonal cavities at a distance of 1.12 Å or 1.35 Å, respectively, from the basal surface on dehydrated mica and above tetrahedral substitutions at a height of 1.93 ± 0.02 Å or 2.15 ± 0.03 Å, respectively, at the highest simulated water coverage of 28 H₂O per ion. The ion’s displacement from a ditrigonal cavity occurs upon adsorption of 2 H₂O per ion for Sr²⁺ and 3 H₂O per ion for Ba²⁺. At a coverage of 28 H₂O per ion, outer-sphere adsorption of Sr²⁺ or Ba²⁺ at a height of 3.9 ± 0.2 Å or 4.17 ± 0.07 Å, respectively, is possible albeit unfavorable on the free energy scale by 107 ± 7 kJ/mol or 89 ± 13 kJ/mol, respectively, as compared to inner-sphere adsorption. Activation energies for the transformation between inner-sphere and outer-sphere adsorptions are calculated to be 121 ± 3 kJ/mol for Sr²⁺ and 99 ± 10 kJ/mol for Ba²⁺. A comparison of these values with those reported recently for Mg²⁺ and Ca²⁺ results in an adsorption affinity sequence Mg²⁺ > Ca²⁺ > Sr²⁺ > Ba²⁺ in agreement with the sequence predicted recently for low dielectric constant solids (which include mica) (Sverjensky, 2006). A recent resonant anomalous X-ray reflectivity study of Sr²⁺ adsorption on muscovite (Park et al., 2006) has questioned the common assumption (Stumm, 1992), which is supported by the present simulation results, that inner-sphere adsorption is stronger than outer-sphere adsorption. A modification of the cleaved muscovite surface as a result of Sr²⁺ adsorption in muscovite’s ditrigonal cavities and related destabilization of muscovite’s hydroxyl groups is proposed as a possible reason for this controversy.     

The estimation of the pK1HA of acids in seawater using the Pitzer equations

The equations of Pitzer have been used to calculate the stoichiometric ionization constants, pK¹_HA, for acids in NaCl media at 25°C. The calculated results for the ionization of HAc, H₂O, B(OH)₃, H₂CO^?₃, H₃PO₄, H₂PO^?₄, HPO^2?₄, H₃AsO₄, H₂AsO^?₄ and HAsO^2?₄ are in good agreement with the measured values, providing higher order interaction terms (θ and ψ) are used. The pK¹_HA measurements of these acids in NaCl media containing Mg²⁺ and Ca²⁺ were used to determine Pitzer specific interaction parameters at I = 0.7. With these Pitzer coefficients, it was possible to make reliable estimates for the activity coefficients of anions in seawater (S = 35) that form strong interactions with Mg²⁺ and Ca²⁺. The calculated activity coefficients yield reliable estimates for the pK¹_HA of acids in seawater.     

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.     

The activity coefficients of Fe(III) hydroxide complexes in NaCl and NaClO4 solutions

The osmotic coefficients of FeCl₃ at 25 °C from 0.15 to 1.7 m [Rumyantsev et al., Z. Phys. Chem., 218, 1089-1127, 2004] have been used to determine the Pitzer parameters (β⁽⁰⁾, β⁽¹⁾ and C^?) for FeCl₃. Since the differences in the Pitzer coefficients of rare earths in NaCl and NaClO₄ are small, the values of Fe(ClO₄)₃ have been estimated using the differences between La(ClO₄)₃ and LaCl₃. The Pitzer coefficients for FeCl₃ combined with enthalpy and heat capacity data for the rare earths can be used to estimate the activity coefficients of Fe³⁺ in NaCl over a wide range of temperatures (0 to 50 °C) and ionic strength (0 to 6 m).The activity coefficients of Fe³⁺ in NaCl and NaClO₄ solutions have been used to determine the activity coefficients of Fe(OH)²⁺ in these solutions from the measured first hydrolysis constants of Fe³⁺ [Byrne et al., Mar. Chem., 97, 34-48, 2005]. The activity coefficients of , Fe(OH)₃ and from 0 to 50 °C have also been determined from the solubility measurements of Fe(III) in NaCl solutions [Liu and Millero, Geochim. Cosmochim Acta, 63, 3487-3497, 1999]. These activity coefficients have been fitted to the Pitzer equations. These results can be used to estimate the speciation of Fe(III) with OH⁻ in natural waters with high concentrations of NaCl from 0 to 50 °C.     

Trace/minor element:calcium ratios in cultured benthic foraminifera. Part I: Inter-species and inter-individual variability

Trace/minor element signatures (D_Cd, D_Ba, D_Mg, and D_Sr) were measured in the tests (shells) of benthic foraminifera cultured in a trace-metal-concentration-controlled system. The culture system was constructed of inert materials and designed to limit microhabitat effects. This system ensured that variation observed in cultured foraminiferal element:calcium (TE/Ca) signatures was due to biologically mediated (vital) effects only. Two species, Bulimina aculeata and Rosalina vilardeboana, reproduced prolifically during two 4-to-8-month culture periods. In every case (i.e., for both species and each element), the inter-individual variability was larger than the analytical precision. Mean (±1 standard deviation) D_E signatures for B. aculeata were: D_Cd: 1.5 ± 0.4, D_Ba × 10: 2.1 ± 0.7, D_Mg × 1000: 0.62 ± 0.15, and D_Sr × 10: 1.5 ± 0.1. Cultured B. aculeata D_Mg, calibrated from culture and core-top (live) field specimens, predicted temperatures within ±2.0 °C. The observed inter-individual variability from culture specimens was as large or larger than comparable results from core-top investigations. R. vilardeboana D_Cd signatures were significantly lower, while D_Ba, D_Mg, and D_Sr signatures were significantly higher than B. aculeata values. Since our culture system minimizes microhabitat variability, the variation in measured TE/Ca ratios suggests that biological processes are a significant factor in inter-individual and inter-species variability. Comparison of cultured and field-collected foraminiferal D_Ba signatures supports previous findings that pore-water chemistry is a major environmental influence on foraminiferal test chemistry.     

Mg2+ removal in the system Mg2+—amorphous SiO2—H2O by adsorption and Mg-hydroxysilicate precipitation

Two chemical processes can remove Mg²⁺ from suspensions containing amorphous silica (am-SiO₂) at low temperatures: adsorption and precipitation of a Mg-hydroxysilicate resembling sepiolite. Mg²⁺ removal from am-SiO₂ suspensions was investigated, and the relative role of the two removal processes evaluated, as a function of: pH, ionic strength, Mg²⁺ concentration, and temperature.The extent of Mg²⁺ adsorption onto am-SiO₂ decreases with increasing NaCl concentration due to displacement of Mg²⁺ by Na⁺. At NaCl concentrations of 0.05 M and above, adsorption occurs only at pH values above 8.5, where rapid dissolution of am-SiO₂ gives rise to high concentrations of dissolved silica, resulting in supersaturation with respect to sepiolite. Removal of Mg²⁺, at concentrations of 40 to 650 μM, from am-SiO₂ suspensions in 0.70 M NaCl at 25 °C occurs at pH 9.0 and above. Experiments show that under these conditions adsorption and Mg-hydroxysilicate precipitation remove Mg²⁺ at similar rates. For 0.05 M Mg²⁺, at 0.70 M ionic strength and 25 °C, measurable Mg²⁺ removal occurs down to ca. pH 7.5 but is primarily due to Mg-hydroxysilicate precipitation. For the same solution conditions at 5°C, Mg²⁺ removal occurs above pH 8.0 and is primarily due to adsorption.Assuming that increasing pressure does not greatly enhance adsorption, Mg²⁺ adsorption onto am-SiO₂ is an insignificant process in sea water. The surface charge of pristine am-SiO₂ in sea water is primarily controlled by interactions with Na⁺. The principal reaction between Mg²⁺ and am-SiO₂ in marine sediments is sepiolite precipitation.The age distribution of sepiolite in siliceous pelagic sediments is influenced by temperatures of bottom waters and by geothermal gradients.     

Isopiestic measurement of the osmotic and activity coefficients for the NaOH-NaAl(OH)4-H2O system at 313.2 K

By using a specially designed and constructed isopiestic apparatus, we measured the osmotic coefficients at 313.2 K for the NaOH-NaAl(OH)₄-H₂O system with the total alkali molality, m_NaOHT (m_NaOH + m_NaAl[OH]4), from 0.05 mol/kg H₂O to 12 mol/kg H₂O and α_K (m_NaOHT/m_NaAl(OH)4) from 1.64 to 5.53. The mean standard deviation of the measurements is 0.0038. Several sets of the Pitzer model parameters for NaOH-NaAl(OH)₄-H₂O system were then obtained by regressing the measured osmotic coefficients with the Pitzer model and the Pitzer model parameters for NaOH(aq). One set of the results is as follows: β⁽⁰⁾_NaOH: 0.08669, β⁽¹⁾_NaOH: 0.31446, β⁽²⁾_NaOH: −0.00007367, C^Φ_NaOH: 0.003180, β⁽⁰⁾_NaAl(OH)4: 0.03507, β⁽¹⁾_NaAl(OH)4: 0.02401, C^Φ_NaAl(OH)4: −0.001066, θ_{OH⁻Al(OH)4⁻}: 0.08177, Ψ_Na⁺OH⁻_Al(OH)4⁻: −0.01162. The mean standard difference between the calculated and the measured osmotic coefficients is 0.0088. With the obtained Pitzer model parameters, we calculated the values of K = (γ_{NaAl(OH)4,cal²} · m_{Al(OH)4⁻,exp})/(γ_NaOH,cal² · m_OH⁻,exp) for the gibbsite solubility. The results show that the obtained Pitzer model parameters are reliable, and the relative error of the calculated activity coefficients should be < 2.1%. We also compared the calculated gibbsite solubility data among several activity coefficients models over a range of m_NaOHT at various temperatures. The comparison indicates that our activity coefficients model may be approximately applied in the ranges of temperature from 298.2 to 323.2 K and m_NaOHT from 0 to 8 mol/kg H₂O. We also calculated the stoichiometric activity coefficients of NaOH and NaAl(OH)₄ and the activity of H₂O for the NaOH-NaAl(OH)₄-H₂O system, and these calculations establish their variations with m_NaOHT and α_K. These variations imply that the strengths of the repulsive interactions among various anions are in the following sequence: Al(OH)₄⁻-Al(OH)₄⁻ < Al(OH)₄⁻-OH⁻ < OH⁻-OH⁻, and the attractive interaction between Al(OH)₄⁻ and H₂O is weaker than that between OH⁻ and H₂O.     

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.     

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

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

Ra/Ra and Ra/Ba ratios to track barite formation and transport in the water column

We measured ²²⁸Ra_ex/²²⁶Ra_ex and ²²⁶Ra_ex/Ba_ex ratios in suspended and sinking particles collected at the Oceanic Flux Program (OFP) time-series site in the western Sargasso Sea and compared them to seawater ratios to provide information on the origin and transport of barite (BaSO₄) in the water column. The ²²⁸Ra_ex/²²⁶Ra_ex ratios of the suspended particles down to 2000 m are nearly identical to those of seawater at the same water depth. These ratios are much lower than expected if suspended barite was produced in surface waters and indicate that barite is produced throughout the mesopelagic layer. The ²²⁸Ra_ex/²²⁶Ra_ex activity ratios of sinking particles collected at 1500 and 3200 m varied mostly between 0.1 and 0.2, which is intermediate between the seawater ratio at these depths (<0.03) and the seawater ratios found in the upper 250 m (0.31-0.42). This suggests that excess Ba (i.e., Ba_ex = Ba_total − Ba_lithogenic), considered to be mainly barite, present in the sinking flux is a mixture of crystals formed recently in the upper water column, formed several years earlier in the upper water column, or formed recently in deeper waters. We observe a sizeable temporal variability in the ²²⁸Ra_ex/²²⁶Ra_ex ratios of sinking particles, which indicates temporal variability in the relative proportion of barite crystals originating from surface (with a high ²²⁸Ra_ex/²²⁶Ra_ex ratio) and mesopelagic (with a low ²²⁸Ra_ex/²²⁶Ra_ex ratio) sources. However, we could not discern a clear pattern that would elucidate the factors that control this variability. The ²²⁶Ra/Ba ratios measured in seawater are consistent with the value reported from the GEOSECS expeditions (2.3 dpm μmol⁻¹) below 500 m depth, but are significantly lower in the upper 500 m. High ²²⁶Ra_ex/Ba_ex ratios and elevated Sr concentrations in suspended particles from the upper water column suggest preferential uptake of ²²⁶Ra over Ba during formation of SrSO₄ skeletons by acantharians, which must contribute to barite formation in shallow waters. Deeper in the water column the ²²⁶Ra_ex/Ba_ex ratios of suspended particles are lower than those of seawater. Since ²²⁸Ra_ex/²²⁶Ra_ex ratios demonstrate that suspended barite at these depths has been produced recently and in situ, their low ²²⁶Ra_ex/Ba_ex ratios indicate preferential uptake of Ba over Ra in barite formed in mesopelagic water.     

Empirical prediction of the Pitzer’s interaction parameters for cationic Al species with both SiO2(aq) and CO2(aq): Implications for the geochemical modelling of very saline solutions

The Pitzer’s interaction parameters, λ_N–M, involving the Mth cationic Al species Al³⁺ or AlOH²⁺ or AlO⁺ and the Nth neutral species SiO_2(aq) (at temperatures of 25–300 °C) or CO_2(aq) (at temperatures of 25–150 °C), have been evaluated through empirical linear relationships between λ_N–M and the surface electrostatic field of the ionic species of interest. These relationships have been obtained starting from the known λ_N–M for both SiO_2(aq) and CO_2(aq) with the main dissolved cations. The Pitzer’s interaction parameter thus estimated for the pair CO_2(aq)–Al³⁺ at 25 °C, 0.327, is 20–40% higher than the corresponding values obtained from CO₂ solubilities in concentrated solutions of AlCl₃, 0.272 ± 0.010 (2σ), and Al₂(SO₄)₃, 0.232 ± 0.002 (2σ), partly corroborating the empirical approach adopted in this study. To test the Pitzer’s interaction coefficients for cationic Al species with aqueous SiO₂, the log K values of the kaolinite dissolution reaction have been computed starting from available experimental data at 23–25 °C and ionic strengths of 0.0001–0.12 mol/kg adopting, alternatively, the Pitzer’s equations and the Debye–Hückel equation. A satisfactory agreement has been found between the log K values obtained through these two approaches, with maximum deviations of 0.11–0.12 log units. This good convergence of results is encouraging as it represents a necessary condition to prove the reliability of the Pitzer’s interaction coefficients estimated in this work. These results are a first step to take into account specific interactions among solutes in concentrated electrolyte solutions, such as those hosted in sedimentary basins or geothermal waters, for instance through the Pitzer’s equations. However, experimental or field data at higher ionic strengths are absolutely necessary to validate the reliability of the Pitzer’s interaction coefficients determined in this study.     

Calcite dissolution kinetics in the system CaCO3-H2O-CO2 at high undersaturation

Dissolution rates of limestone covered by a water film open to a CO₂-containing atmosphere are controlled by the chemical composition of the CaCO₃-H₂O-CO₂ solution at the water-mineral interface. This composition is determined by the Ca²⁺-concentration at this boundary, conversion of CO₂ into H⁺ and in the solution, and by diffusional mass transport of the dissolved species from and towards the water-limestone interface. A system of coupled diffusion-reaction equations for Ca²⁺, , and CO₂ is derived. The Ca²⁺ flux rates at the surface of the mineral are defined by the PWP-empirical rate law. These flux rates by the rules of stoichiometry must be equal to the flux rates of CO₂ across the air-water interface. In the solution, CO₂ is converted into H⁺ and . At low water-film thickness this reaction becomes rate limiting. The time dependent diffusion-reaction equations are solved for free drift dissolution by a finite-difference scheme, to obtain the dissolution rate of calcite as a function of the average calcium concentration in the water film. Dissolution rates are obtained for high undersaturation. The results reveal two regimes of linear dissolution kinetics, which can be described by a rate law F = α_i(m_ic_eq − c), where c is the calcium concentration in the water film, c_eq the equilibrium concentration with respect to calcite. For index i = 0, a fast rate law, which here is reported for the first time, is found with α₀ = 3 × 10⁻⁶ m s⁻¹ and m₀ = 0.3. For c > m₀c_eq, a slow rate law is valid with α₁ = 3 × 10⁻⁷ m  s⁻¹ and m₁ = 1, which confirms earlier work. The numbers given above are valid for film thickness of several tenths of a millimetre and at 20 °C. These rates are proven experimentally, using a flat inclined limestone plate covered by a laminar flowing water film injected at an input point with known flow rate Q and calcium concentration. From the concentration measured after flow distance x the dissolution rates are determined. These experiments have been performed at a carbon-dioxide pressure of 0.00035 atm and also of 0.01 atm. The results are in good agreement to the theoretical predictions.     

Complex formation in the ternary system Mg(II)-CO2-H2O

The carbonato and hydrogencarbonato complexes of Mg²⁺ were investigated at 25 and 50° in solutions of the constant ClO₄^? molality (3 M) consisting preponderantly of NaClO₄. The experimental data could be explained assuming the following equilibria: Mg²⁺ + CO_2B + H₂O ag MgHCO⁺₃ + H⁺, $\text{log}{}^1\text{β}{}_1\text{= ?7.64}{}_4\text{± 0.01}{}_7\text{(25°), ?7.46}{}_2\text{± 0.01}{}_1\text{(50°)}$, Mg²⁺ + 2 CO_2g + 2 H₂Oag Mg(HCO₃)⁰₂ ± 2 H⁺, $\text{log}{}^1\text{β}{}_2\text{= ?15.00 ± 0.14 (25°)}$, ?15.37 ± 0.39 (50°), Mg²⁺ + CO_2g + H₂Oag MgCO⁰₃ + 2 H⁺, $\text{log}{}^1\text{k}{}_1\text{= ?15.64 ± 0.06 (25°)}$,?15.23 ± 0.02 (50°), with the assumption γ_MgCO3⁰ = γ_Mg(HCO3)⁰2, ΔG⁰(I = 0) for the reaction MgCO⁰₃ + CO_2g + H₂O = Mg(HCO₃)⁰₂ was estimated to be ?3.9₁ ± 0.8₆ and 0.6 ± 2.4 kJ/mol at 25 and 50°C, respectively. The abundance of carbonate linked Mg(II) species in fresh water systems is discussed.     

Partitioning of Ni2+, Co2+, Fe2+, Mn2+ and Mg2+ between olivine and silicate melts: compositional dependence of partition coefficient

Partitioning of Ni²⁺, Co²⁺, Fe²⁺, Mn²⁺ and Mg²⁺ between olivine and silicate melts has been determined near the join (Mg_0.5·-Fe_0.5)₂SiO₄-K₂O 4SiO₂ and for seven different basaltic compositions. The experiments were made at 1 atm total pressure, 1500-1100°C, and under moderate to reducing oxygen fugacities. The concentration factor, defined as K^MO = (MO)^ol/(MO)^liq (molar ratio), increases markedly for all the cations studied as the olivine component of the liquid decreases. Much of the increase in K^MO is considered as due to the compositional effect of the coexisting liquid: the temperature effect on K^MO is probably opposite to the compositional effect (K^MO decreases as temperature decreases).The partition coefficient K^MO-MgO = (MO/MgO)^ol/(MO/MgO)^liq for the reaction, M_ol²⁺ + Mg_liq²⁺ = M_liq²⁺ + Mg_ol²⁺. is relatively constant over a wide range of SiO₂ content of the liquid, except in the case of Ni²⁺. The partition coefficients have similar ranges both in synthetic and natural rock systems: K^NiO-MgO = 1.8–3.0, K^CoO-MgO = 0.6–0.8, K^FeO-MgO = 0.27–0.38, and K^MnO-MgO = 0.23–0.32. There is a systematic variation in the partition coefficient K^MO-MgO with the composition of liquid; K^MO-MgO increases with increasing SiO₂ content of melt. The partition coefficients for Co²⁺, Fe²⁺ and Mn²⁺ are useful to test the equilibration of olivine with magma of a wide compositional range.     

Distribution of barium and fulvic acid at the mica-solution interface using in-situ X-ray reflectivity

The interfacial structures of the basal surface of muscovite mica in solutions containing (1) 5 × 10⁻³ m BaCl₂, (2) 500 ppm Elliott Soil Fulvic Acid I (ESFA I), (3) 100 ppm Elliott Soil Fulvic Acid II (ESFA II), (4) 100 ppm Pahokee Peat Fulvic Acid I (PPFA), and (5) 5 × 10⁻³ m BaCl₂ and 100 ppm ESFA II were obtained with high resolution in-situ X-ray reflectivity. The derived electron-density profile in BaCl₂ shows two sharp peaks near the mica surface at 1.98(2) and 3.02(4) Å corresponding to the heights of a mixture of Ba²⁺ ions and water molecules adsorbed in ditrigonal cavities and water molecules coordinated to the Ba²⁺ ions, respectively. This pattern indicates that most Ba²⁺ ions are adsorbed on the mica surface as inner-sphere complexes in a partially hydrated form. The amount of Ba²⁺ ions in the ditrigonal cavities compensates more than 90% of the layer charge of the mica surface. The electron-density profiles of the fulvic acids (FAs) adsorbed on the mica surface, in the absence of Ba²⁺, had overall thicknesses of 4.9-10.8 Å and consisted of one broad taller peak near the surface (likely hydrophobic and positively-charged groups) followed by a broad humped pattern (possibly containing negatively-charged functional groups). The total interfacial electron density and thickness of the FA layer increased as the solution FA concentration increased. The sorbed peat FA which has higher ash content showed a higher average electron density than the sorbed soil FA. When the muscovite reacted with a pre-mixed BaCl₂-ESFA II solution, the positions of the two peaks nearest the surface matched those in the BaCl₂ solution. However, the occupancy of the second peak decreased by about 30% implying that the hydration shell of surface-adsorbed Ba²⁺ was partially substituted by FA. The two surface peaks were followed by a broad less electron-dense layer suggesting a sorption mechanism in which Ba²⁺ acts dominantly as a bridging cation between the mica surface and FA. When the muscovite reacted first with FA and subsequently with BaCl₂, more Ba²⁺ could be adsorbed on the FA-coated mica surface. The peak closest to the mica included Ba²⁺ ions adsorbed directly on the mica in an amount similar to that in the BaCl₂ solution but more broadly distributed. A second peak observed within the FA layer suggests that the FA coating provides additional sites for Ba²⁺ sorption. The results indicate that enhanced uptake of heavy metals can occur when an organic coating already exists on a mineral surface.     

An isopiestic study of aqueous NaBr and KBr at 50 °C: Chemical equilibrium model of solution behavior and solubility in the NaBr-H2O, KBr-H2O and Na-K-Br-H2O systems to high concentration and temperature

The isopiestic method has been used to determine the osmotic coefficients of the binary solutions NaBr-H₂O (from 0.745 to 5.953 mol kg⁻¹) and KBr-H₂O (from 0.741 to 5.683 mol kg⁻¹) at the temperature t = 50 °C. Sodium chloride solutions have been used as isopiestic reference standards. The isopiestic results obtained have been combined with all other experimental thermodynamic quantities available in literature (osmotic coefficients, water activities, bromide mineral’s solubilities) to construct a chemical model that calculates solute and solvent activities and solid-liquid equilibria in the NaBr-H₂O, KBr-H₂O and Na-K-Br-H₂O systems from dilute to high solution concentration within the 0-300 °C temperature range. The Harvie and Weare [Harvie C., and Weare J. (1980) The prediction of mineral solubilities in naturalwaters: the Na-K-Mg-Ca-Cl-SO₄-H₂O system from zero to high concentration at 25 °C. Geochim. Cosmochim. Acta44, 981-997] solubility modeling approach, incorporating their implementation of the concentration-dependent specific interaction equations of Pitzer [Pitzer K. (1973) Thermodynamics of electrolytes. I. Theoretical basis and general equations. J. Phys. Chem.77, 268-277] is employed. The model for binary systems is validated by comparing activity coefficient predictions with those given in literature, and not used in the parameterization process. Limitations of the mixed solutions model due to data insufficiencies are discussed. This model expands the variable temperature sodium-potassium model of Greenberg and Moller [Greenberg J., and Moller N. (1989) The prediction of mineral solubilities in natural waters: a chemical equilibrium model for the Na-K-Ca-Cl-SO₄-H₂O system to high concentration from 0 to 250 °C. Geochim. Cosmochim. Acta53, 2503-2518] by evaluating Br⁻ pure electrolyte and mixing solution parameters and the chemical potentials of three bromide solid phases: NaBr-2H₂O (cr), NaBr (cr) and KBr (cr).     

An atomic force microscopy study of calcite dissolution in saline solutions: The role of magnesium ions

In situ Atomic Force Microscopy, AFM, experiments have been carried out using calcite cleavage surfaces in contact with solutions of MgSO₄, MgCl₂, Na₂SO₄ and NaCl in order to attempt to understand the role of Mg²⁺ during calcite dissolution. Although previous work has indicated that magnesium inhibits calcite dissolution, quantitative AFM analyses show that despite the fact that Mg²⁺ inhibits etch pit spreading, it increases the density and depth of etch pits nucleated on calcite surfaces and, subsequently, the overall dissolution rates: i.e., from 10^−11.75 mol cm⁻² s⁻¹ (in deionized water) up to 10^−10.54 mol cm⁻² s⁻¹ (in 2.8 M MgSO₄). Such an effect is concentration-dependent and it is most evident in concentrated solutions ([Mg²⁺] >> 50 mM). These results show that common soluble salts (especially Mg sulfates) may play a critical role in the chemical weathering of carbonate rocks in nature as well as in the decay of carbonate stone in buildings and statuary.     

Ra/Ra and Ra/Ba ratios in the Western Mediterranean Sea: Barite formation and transport in the water column

²²⁶Ra, ²²⁸Ra and Ba distributions as well as ²²⁸Ra/²²⁶Ra and ²²⁶Ra/Ba ratios were measured in seawater, suspended and sinking particles at the DYFAMED station in the Western Mediterranean Sea at different seasons of year 2003 in order to track the build-up and fate of barite through time. The study of the ²²⁸Ra_ex/²²⁶Ra_ex ratios (Ra_ex = Ra activities corrected for the lithogenic Ra) of suspended particles suggests that Ba_ex (Ba_ex = Ba concentrations corrected for the lithogenic Ba, mostly barite) formation takes place not only in the upper 500 m of the water column but also deeper (i.e. throughout the mesopelagic layer). Temporal changes in the ²²⁸Ra_ex/²²⁶Ra_ex ratios of sinking particles collected at 1000 m depth likely reflect changes in the relative proportion of barite originating from the upper water column (with a high ²²⁸Ra/²²⁶Ra ratio) and formed in the mesopelagic layer (with a low ²²⁸Ra/²²⁶Ra ratio). ²²⁸Ra_ex/²²⁶Ra_ex ratios measured in sinking particles collected in the 1000 m-trap in April and May suggest that barite predominantly formed in the upper water column during that period, while barite found outside the phytoplankton bloom period (February and June) appears to form deeper in the water column. Combining ratios of both the suspended and sinking particles provides information on aggregation/disaggregation processes. High ²²⁶Ra_ex/Ba_ex ratios were also found in suspended particles collected in the upper 500 m of the water column. Because celestite is expected to be enriched in Ra [Bernstein R. E., Byrne R. H. and Schijf J. (1998) Acantharians: a missing link in the oceanic biogeochemistry of barium. Deep-Sea Res. II45, 491-505], acantharian skeletons may contribute to these high ratios in shallow waters. The formation of both acantharian skeletons and barite enriched in ²²⁶Ra may thus contribute to the decrease in the dissolved ²²⁶Ra activity and ²²⁶Ra/Ba ratios of surface waters observed between February and June 2003 at the DYFAMED station.     

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.     

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