Partitioning coefficients between olivine and silicate melts

Variation of Nernst partition coefficients (D) between olivine and silicate melts cannot be neglected when modeling partial melting and fractional crystallization. Published natural and experimental ^{olivine/liquid}D data were examined for covariation with pressure, temperature, olivine forsterite content, and melt SiO₂, H₂O, MgO and MgO/MgO + FeO^total. Values of ^{olivine/liquid}D generally increase with decreasing temperature and melt MgO content, and with increasing melt SiO₂ content, but generally show poor correlations with other variables. Multi-element ^{olivine/liquid}D profiles calculated from regressions of D REE–Sc–Y vs. melt MgO content are compared to results of the Lattice Strain Model to link melt MgO and: D₀ (the strain compensated partition coefficient), E_M³⁺ (Young's Modulus), and r₀ (the size of the M site). Ln D₀ varies linearly with Ln MgO in the melt; E_M³⁺ varies linearly with melt MgO, with a dog-leg at ca. 1.5% MgO; and r₀ remains constant at 0.807 Å. These equations are then used to calculate ^{olivine/liquid}D for these elements using the Lattice Strain Model. These empirical parameterizations of ^{olivine/liquid}D variations yield results comparable to experimental or natural partitioning data, and can easily be integrated into existing trace element modeling algorithms. The ^{olivine/liquid}D data suggest that basaltic melts in equilibrium with pure olivine may acquire small negative Ta–Hf–Zr–Ti anomalies, but that negative Nb anomalies are unlikely to develop. Misfits between results of the Lattice Strain Model and most light rare earth and large ion lithophile partitioning data suggest that kinetic effects may limit the lower value of D for extremely incompatible elements in natural situations characterized by high cooling/crystallization rates.     

Distribution of Mg2+ between olivine and silicate melt,and its implications regarding melt structure

Consideration of experimental data on the distribution of Mg²⁺ between olivine and silicate liquid clearly demonstrates that the distribution coefficient (K_Mg) is dependent upon variations in temperature, pressure and melt composition, largely because these variables control the solubility of Mg²⁺ in the melt phase. Attempts to minimize composition dependence of K_Mg, utilizing various activity-composition models for silicate melts, have been partially successful. Composition-related effects do not appear to be large, however, for melts of restricted range in composition (e.g., tholeiitic or lunar basalts) as long as the contents of alkalis and the alkali/alumina ratio are relatively small (on a molar basis). For such melts, K_Mg may be used as a reliable geothermometer. By analogy, these conclusions can be extended to the distribution of other divalent cation such as Fe²⁺, Mn²⁺, Ni²⁺ and Co²⁺.     

An electrochemical study of Ni2+, Co2+, and Zn2+ ions in melts of composition CaMgSi2O6

Cyclic voltammetry has been done for Ni²⁺, Co²⁺, and Zn²⁺ in melts of diopside composition in the temperature range 1425 to 1575°C. Voltammetric curves for all three ions excellently match theoretical curves for uncomplicated, reversible charge transfer at the Pt electrode. This implies that the neutral metal atoms remain dissolved in the melt. The reference electrode is a form of oxygen electrode. Relative to that reference assigned a reduction potential of 0.00 volt, the values of standard reduction potential for the ions are $\text{E}{}^1\text{(}\text{Ni}{}^{\mathrm{2+}}\text{Ni}{}^0\text{,}\text{diopside}\text{, 1500°C) = ?0.32 ± .01}\text{V}$, $\text{E}{}^1\text{(}\text{Co}{}^{\mathrm{2+}}\text{Co}{}^0\text{,}\text{diopside}\text{, 1500°C) = ?0.45 ± .02}\text{V}$, and $\text{E}{}^1\text{(}\text{Zn}{}^{\mathrm{2+}}\text{Zn}{}^0\text{,}\text{diopside}\text{, 1500°C) = ?0.53 ± .01}\text{V}$. The electrode reactions are rapid, with first order rate constants of the order of 10^?2 cm/sec. Diffusion coefficients were found to be 2.6 × 10^?6 cm²/sec for Ni²⁺, 3.4 × 10^?6 cm²/sec for Co²⁺, and 3.8 × 10^?6 cm²/sec for Zn²⁺ at 1500°C. The value of $\text{E1}$ ($\text{Ni}{}^{\mathrm{2+}}\text{Ni}{}^0$, diopside) is a linear function of temperature over the range studied, with values of ?0.35 V at 1425°C and ?0.29 V at 1575°C. At constant temperature the value of $\text{E}{}^1\text{(}\text{Ni}{}^{\mathrm{2+}}\text{Ni}{}^0\text{, 1525°C}$) was not observed to vary with composition over the range CaO · MgO · 2SiO₂ to CaO·MgO·3SiO₂ or from 1.67 CaO·0.33MgO·2SiO₂ to 0.5 CaO·1.5MgO·2SiO₂. The value for the diffusion coefficient for Ni²⁺ decreased by an order of magnitude at 1525°C over the compositional range CaO · MgO · 1.25SiO₂ to CaO · MgO · 3SiO₂. This is consistent with a mechanism by which Ni²⁺ ions diffuse by moving from one octahedral coordination site to another in the melt, with the same Ni²⁺ species discharging at the cathode regardless of the SiO₂ concentration in the melt.     

Microscope-photometric methods for non-destructive Fe2+ – Fe3+ determination in chloritoids (Fe2+, Mn2+, Mg)2 (Al,Fe3+)4 Si2O10(OH)4

Six natural chloritoid crystals with Fe²⁺ and Fe³⁺ contents ranging from 4.15 to 12.81 and from 0.411 to 0.849g-atoms/l, respectively, as determined by means of microprobe and Mössbauer techniques, served as reference material to develop non-destructive microscope-spectrophotometric methods for quantitative Fe²⁺ – Fe³⁺ determinations in chloritoids from unpolarized spectra of (001) platelets. Fe²⁺ concentrations in g-atom/l can be obtained from [ [Fe³⁺]=C₁xD₁/t where D₁ = log₁₀(I_0/I at 28,000 cm^-1 and t=crystal thickness in cm; C₁ is a conttant that may be influenced somewhat by experimental conditions and is found to be 0.002289 with the experimental set-up used in this study. Fe²⁺ concentrations in g-atom/l can be obtained from [Fe²⁺]=C₁xD₁/D₁-C₃ with D₂=log₁₀(I₀/I) at 16,300 cm^?1 and constants C₄ = 45.36 and C₅ = 3.540. Due to the uncertainties in absorbance measurements, D₁ and D₂ and the thickness measurements, the accuracies are ±0.05 and ±0.15 g-atom/l for [Fe³⁺] and [Fe²⁺], respectively. The determinations may be carried out on chloritoid grains in normal thin sections with an areal resolution of ～10 μm.     

Polymerization of silicate and aluminate tetrahedra in glasses,melts and aqueous solutions—II. The network modifying effects of Mg2+ , K+, Na+, Li+, H+, OH−, F−, Cl−, H2O,CO2 and H3O+ on silicate polymers

The effect of the group IA and VIIA ions, as well as Mg²⁺, and the molecules H₂O, CO₂, H₃O⁺ and OH^? on the energy of the Si-O bond in a H₆Si₂O₇ cluster has been calculated using semiempirical molecular orbital calculations (CNDO/2). Three types of elementary processes, i.e. substitution, addition, and polymerization reactions have been used to interpret data on the dynamic viscosity, surface tension and surface charge, hydrolytic weakening, diffusivity, conductivity, freezing point depression, and degree of polymerization of silicates in melts, glasses, and aqueous solutions. As a test of our calculational procedure, observed X-ray emission spectra of binary alkali silicate glasses were compared with calculated electronic spectra. The well known bondlength variations between the bridging bond [Si-O(br)] and the non-bridging bond [Si-O(nbr)] in alkali silicates are shown to be due to the propagation of oscillating bond-energy patterns through the silica framework. A kinetic interpretation of some results of our calculations is given in terms of the Bell-Evans-Polanyi reaction principle.     

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.     

Pressure dependence on partition coefficients for trace elements between olivine and the coexisting melts

Partition coefficients between olivine and melt at upper mantle conditions, 3 to 14 GPa, have been determined for 27 trace elements (Li, Be, B, Na, Mg, Al, Si, P, K, Ca, Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, Rb, Sr, Y, Zr, Cs, Ba, La and Ce) using secondary-ion mass-spectrometry (SIMS) and electron-probe microanalysis (EPMA). The general pattern of olivine/melt partitioning on Onuma diagrams resembles those reported previously for natural systems. This agreement strongly supports the argument that partitioning is under structural control of olivine even at high pressure. The partition coefficients for mono- and tri-valent cations show significant pressure dependence, both becoming larger with pressure, and are strongly correlated with coupled substitution into cation sites in the olivine structure. The dominant type of trace element substitution for mono- and tri-valent cations into olivine changes gradually from (Si, Mg)↔(Al, Cr) at low pressure to (Si, Mg)↔(Al, Al) and (Mg, Mg)↔(Na, Al) at high pressure. The change in substitution type results in an increase in partition coefficients of Al and Na with pressure. An inverse correlation between the partition coefficients for divalent cations and pressure has been observed, especially for Ni, Co and Fe. The order of decreasing rate of partition coefficient with pressure correlates to strength of crystal field effect of the cation. The pressure dependence of olivine/melt partitioning can be attributed to the compression of cation polyhedra induced by pressure and the compensation of electrostatic valence by cation substitution. Received: March 6, 1997 / Revised, accepted: March 12, 1998     

水锰矿与Fe2+的相互作用与转化过程

作为表生土壤环境中易生成且分布广泛的氧化锰矿物,水锰矿(γ-MnOOH)能参与铁氧化物的生成过程,影响Fe_(2+)的迁移、转化和归趋。本文考察了pH值为3.0~7.0的模拟水溶液体系中水锰矿与Fe_(2+)的相互作用及其生成铁氧化物的过程,分析了Fe_(2+)浓度、pH值和空气(氧气)对Fe(Ⅲ)氧化物晶体结构类型、化学组成和反应速率的影响。研究结果表明,水锰矿氧化Fe_(2+)产物主要为针铁矿和纤铁矿;pH值为3.0~5.0时产物为针铁矿,而pH值为7.0时产物为针铁矿与纤铁矿的混合物,且高浓度Fe_(2+)会促使纤铁矿生成;引入空气利于针铁矿生成;反应速率随着pH值升高、氧气分压的增大而加快。本工作丰富了对铁氧化物在常见锰氧化物表面形成和转化过程的认识。     

The incorporation of Mg2+ and Sr2+ into calcite overgrowths: influences of growth rate and solution composition

The concentrations of Mg²⁺ and Sr²⁺ incorporated within calcite overgrowths precipitated from seawater and related solutions, determined at 25°C, were independent of the precipitation rate over approximately an order of magnitude. The saturation states used to produce this range of precipitation rates varied from 3 to 17 depending on the composition of the solution.The amount of Mg²⁺ incorporated in the overgrowths was not directly proportional to $\text{Mg}{}^{\mathrm{2+}}\text{Ca}{}^{\mathrm{2+}}$ in solution over the entire range (1–20) of ratios studied. Below a ratio of 7.5, the overgrowth was enriched in MgCO₃ relative to what is predicted by the constant distribution coefficient measured above a ratio of 7.5. This increased MgCO₃ correlates with the relative enrichment of adsorbed Mg²⁺. Above a ratio of 7.5 the concentration of MgCO₃ in the calcite overgrowths followed a classical thermodynamic behavior characterized by a constant distribution coefficient of 0.0123 (±0.008 std dev).The concentration of SrCO₃ incorporated in the overgrowths was linearly related to the MgCO₃ content of the overgrowths, and is attributed to increased solubility of SrCO₃ in calcite due to the incorporation of the smaller Mg²⁺ ions.The kinetic data indicate that the growth mechanism involves the adsorption of the cations on the surface of the calcite prior to dehydration and final incorporation. It is suggested that dehydration of cations at the surface is the rate controlling step.     

Thermodynamic model of natural brines accounting for the presence of trace components: II. System Na+, K+, Mg2+ ‖ Cl?, Br?-H2O

This paper reports the calculations of parameters for the Pitzer equation and thermodynamic potentials of solid phases crystallizing in water-salt systems modeling chloride brines taking into account the presence of bromide ions in them. Solubility diagrams were calculated for corresponding ternary and quaternary systems containing chlorides and bromides of sodium, potassium, and magnesium at 25°C. The results of calculations are in adequate agreement with the available published experimental data on solubility and can be used to model salt crystallization during the concentration of seawater and brines. Original Russian Text ? M.V. Charykova, N.A. Charykov, 2007, published in Geokhimiya, 2007, No. 10, pp. 1129–1138.     

Interaction between sulphide and H2O in silicate melts

Reaction between dissolved water and sulphide was experimentally investigated in soda-lime-silicate (NCS) and sodium trisilicate (NS3) melts at temperatures from 1000 to 1200 °C and pressures of 100 or 200 MPa in internally heated gas pressure vessels. Diffusion couple experiments were conducted at water-undersaturated conditions with one half of the couple being doped with sulphide (added as FeS or Na₂S; 1500-2000 ppm S by weight) and the other with H₂O (∼3.0 wt.%). Additionally, two experiments were performed using a dry NCS glass cylinder and a free H₂O fluid. Here, the melt was water-saturated at least at the melt/fluid interface. Profiling by electron microprobe (sulphur) and infrared microscopy (H₂O) demonstrate that H₂O diffusion in the melts is faster by 1.5-2.3 orders of magnitude than sulphur diffusion and, hence, H₂O can be considered as a rapidly diffusing oxidant while sulphur is quasi immobile in these experiments.In Raman spectra a band at 2576 cm⁻¹ appears in the sulphide - H₂O transition zone which is attributed to fundamental S-H stretching vibrations. Formation of new IR absorption bands at 5025 cm⁻¹ (on expense of the combination band of molecular H₂O at 5225 cm⁻¹) and at 3400 cm⁻¹ was observed at the front of the in-diffusing water in the sulphide bearing melt. The appearance and intensity of these two IR bands is correlated with systematic changes in S K-edge XANES spectra. A pre-edge excitation at 2466.5 eV grows with increasing H₂O concentration while the sulphide peak at 2474.0 eV decreases in intensity relative to the peak at 2477.0 eV and the feature at 2472.3 eV becomes more pronounced (all energies are relative to the sulphate excitation, calibrated to 2482.5 eV). The observations by Raman, IR and XANES spectroscopy indicate a well coordinated S²⁻ - H₂O complex which was probably formed in the glasses during cooling at the glass transition. No oxidation of sulphide was observed in any of the diffusion couple experiments. On the contrary, XANES spectra from experiments conducted with a free H₂O fluid show complete transformation of sulphide to sulphate near the melt surface and coexistence of sulphate and sulphide in the center of the melt. This can be explained by a lower H₂O activity in the diffusion couple experiments or by the need of a sink for hydrogen (e.g., a fluid which can dissolve high concentration of hydrogen) to promote oxidation of sulphide by H₂O via the reaction S²⁻ + 4H₂O = SO₄²⁻ + 4H₂. Sulphite could not be detected in any of the XANES spectra implying that this species, if it exists in the melt, it is a subordinate or transient species only.     

溶出珊瑚砂中Ca~(2+)、Mg~(2+)的动力学实验研究及其影响因素

通过改变固液比、摇床转速、珊瑚砂粒径、温度、溶液pH值及溶液含盐量等参数,对珊瑚砂在水溶液中溶出Ca~(2+)、Mg~(2+)进行了实验,以探讨溶出过程中的动力学规律和影响因素。实验结果表明,珊瑚砂中Ca~(2+)、Mg~(2+)溶出量随反应时间逐渐增大;摇床转速越快、固液比越大、温度越高、溶液pH值越低,Ca~(2+)、Mg~(2+)溶出量越大;当珊瑚砂粒径为2.36~4.75 mm、溶液含盐量为100 mg/L时,Ca~(2+)、Mg~(2+)溶出量最大。统计分析表明,摇床转速、温度及溶液p H值均对珊瑚砂溶出有显著影响,但溶液p H值影响最大。珊瑚砂在水中的溶出过程符合收缩核内扩散模型,表明控制整个溶出过程反应速率的决定因素是内扩散速率;在15~40℃时,珊瑚砂在纯水中溶出Ca~(2+)、Mg~(2+)的活化能分别为78.07和74.91 k J/mol。     

多硅白云母中Fe2+的高压有序效应

21个多硅白云母来自4个低温高压变质带。用电子探针、X射线粉末衍射及穆斯堡尔谱测定其化学成分、b0值及Fe²⁺占位。建立Fe²⁺(M1)/Fe²⁺(M2)对b0值相关图,发现Fe²⁺在八面体晶位有序化并解释其有序机理。     

Mg2+调控低温环境下单水方解石的合成

在4℃的低温条件下,通过模拟纳木错湖水中主要离子种类,控制反应体系中的Mg/Ca摩尔比,探究纳木错单水方解石(MHC)的无机成因,并利用X-射线衍射(XRD)、傅里叶变换红外光谱(FTIR)和电感耦合等离子体发射光谱(ICP-OES)对实验沉淀产物进行了分析。结果表明,单水方解石是通过无定形碳酸钙(ACC)转化而来的,且单水方解石的形成和稳定受溶液中的Mg2+调控。当反应溶液中Mg/Ca 2时,生成方解石和镁方解石;当反应溶液中Mg/Ca 2时,单水方解石生成,并且在溶液中的存在时间随着Mg/Ca摩尔比的增加而增加。Mg/Ca=2时,单水方解石在溶液中仅存在48 h,随后转化为文石和少量镁方解石。单水方解石形成过程中Mg2+的作用主要是抑制方解石成核和生长,使单水方解石能够在溶液中形成并稳定存在。此结果为纳木错单水方解石无机成因提供一定的研究基础。     

Oxide solution mechanisms in silicate melts: Systematic variations in the activity coefficient of SiO2

Phase equilibria and spectroscopic data are used to develop a simple model for the interaction of various oxide components and molten SiO₂. Network modifying oxides, M_xO_y produce nonbridging oxygens thereby depolymerizing the SiO₂ network. The energetics of nonbridging oxygen formation are least favorable when the field strength of the metal cation is high. This produces relatively strong M-O and Si-O-Si bridging bonds at the expense of weaker Si-O-M bonds (De Jonget al., 1980). This relationship is manifested by an increase in positive deviations from ideality with increasing cation field strength in M_xO_y-SiO₂ systems; the activity coefficient of SiO₂ is inversely correlated with Si-O-M bond strength. Network forming oxides (aluminates, phosphates, titanates. zirconates, etc.) may copolymerize with the SiO₂ network. Mixing on the same quasi-lattice produces solutions which approach ideality. Deviations from ideality in such solutions can be linked to distortions in the SiO₂ network. Discrete anion formers (phosphates, titanates, chromates, zirconates) complex with metal oxides other than SiO₂ to form discrete structural units which do not copolymerize with SiO₂. The SiO₂ network is essentially shielded from the high charge density cations in such systems and unmixing is common. As a result, the relative deviations from ideality in such melts are high. It is important to recognize that oxides such as P₂O₅, TiO₂ and ZrO₂ may act as either network-formers or discrete anion formers depending upon melt composition, and are probably distributed between these two “sites” in most geologically important liquids. The latter structural role is favored in more basic compositions.     

Partitioning of trace elements between rutile and silicate melts: Implications for subduction zones

Mineral/melt trace element partition coefficients were determined for rutile (TiO₂) for a large number of trace elements (Zr, Hf, Nb, Ta, V, Co, Cu, Zn, Sr, REE, Cr, Sb, W, U, Th). Whilst the high field strength elements (Zr, Hf, Nb, Ta) are compatible in rutile, other studied trace elements are incompatible (Sr, Th, REE). In all experiments we found D_Ta > D_Nb, D_Hf > D_Zr and D_U > D_Th. Partition coefficients for some polyvalent elements (Sb, W, and Co) were sensitive to oxygen fugacity. Melt composition exerts a strong influence on HFSE partition coefficients. With increasing polymerization of the melt, rutile/melt partition coefficients for the high field strength elements Zr, Hf, Nb and Ta increase about an order of magnitude. However, D_Nb/D_Ta and D_Hf/D_Zr are not significantly affected by melt composition. Because D_U ? D_Th, partial melting of rutile-bearing eclogite in subducted lithosphere may cause excesses of ²³⁰Th over ²³⁸U in some island arc lavas, whereas dehydration of subducted lithosphere may cause excesses of ²³⁸U over ²³⁰Th. From our partitioning results we infer partition coefficients for protactinium (Pa) which we predict to be much lower than previously anticipated. Contrary to previous studies, our data imply that rutile should not significantly influence observed ²³¹Pa-²³⁵U disequilibria in certain volcanic rocks.     

热处理褐铁矿去除水中的Mn2+

本文利用褐铁矿中针铁矿经热脱水相变获得以纳米晶赤铁矿为主要物相的纳米-微米多级孔结构材料,并用于模拟净化富Mn~(2+)地下水。同时考察了热处理温度、初始pH值、初始Mn~(2+)浓度、吸附反应时间等对材料去除溶液中Mn~(2+)的影响。XRD、TEM、BET表征结果表明,300℃热处理产物中赤铁矿孔径最小为2.7 nm,比表面积最大达到107.4 m~2/g。吸附实验结果表明,在pH值5~10的范围内,p H值对煅烧褐铁矿颗粒对Mn~(2+)去除效果影响较小;材料在贫氧条件下对水中低浓度Mn~(2+)的最大吸附量为6.45 mg/g;吸附动力学符合准二级动力学模型;褐铁矿热处理形成的纳米晶赤铁矿对Mn~(2+)具有吸附和催化氧化作用,其中的杂质锰氧化物对Mn~(2+)的吸附和催化氧化具有增强作用。     

The effect of melt composition on the wetting angle between silicate melts and olivine

The wetting angle between silicate melts containing Ca, Li, Na, or K and olivine single crystals have been measured as part of an investigation of the dependence of the solid-liquid interfacial energy on melt composition and olivine orientation. The wetting angle increases with increasing silica content of the melt on (100) surfaces, but decreases with increasing silica content on (010) and (001) surfaces. For a given silica content, the wetting angle on (100) decreases in going from Ca to Li to Na to K, while the wetting angle on (010) and (001) increases in going from Ca to K-bearing melts. Based on published values for liquid-vapor interfacial energies, the observed changes in wetting angle with changes in melt composition indicate that the solid-liquid interfacial energy increases with increasing silica content of the melt for the (100) surface. However, for (010) and (001) surfaces, the variation of the solid-liquid interfacial energy with silica content depends upon whether Ca or K is present in the melt. In addition, the solid-liquid interfacial energy depends upon the orientation of the olivine in the following manner: _sl ⁽⁰¹⁰⁾ _sl ⁽⁰⁰¹⁾ _sl ⁽¹⁰⁰⁾ .