Dual speciation of nitrogen in silicate melts at high pressure and temperature: An experimental study

Solubility and speciation of nitrogen in silicate melts have been investigated between 1400 and 1700 °C and at pressures ranging from 10 to 30 kbar for six different binary alkali and alkaline-earth silicate liquids and a Ca-Mg-alumino silicate. Experiments were performed in a piston-cylinder apparatus. The nitrogen source is silver azide, which breaks down to Ag and molecular N₂ below 300 °C. At high pressure and temperature, the nitrogen content may be as high as 0.7 wt% depending on the melt composition, pressure, and temperature. It increases with T, P and the polymerization state of the liquid. Characterization by Raman spectroscopy and ¹⁵N solid state MAS NMR indicates that nitrogen is not only physically dissolved as N₂ within the melt structure like noble gases, but a fraction of nitrogen interacts strongly with the silicate network. The most likely nitrogen-bearing species that can account for Raman and NMR results is nitrosyl group. Solubility data follow an apparent Henry’s law behavior and are in good agreement with previous studies when the nitrosyl content is low. On the other hand, a significant departure from a Henry’s law behavior is observed for highly depolymerized melts, which contain more nitrosyl than polymerized melts. Possible solubility mechanisms are also discussed. Finally, a multi-variant empirical relation is given to predict the relative content of nitrosyl and molecular nitrogen as a function of P, T, and melt composition and structure. This complex speciation of nitrogen in melts under high pressure may have significant implication concerning crystal-melt partitioning of nitrogen as well as for potential elemental and isotopic fractionation of nitrogen in the deep Earth.     

The solubility of gold in H2O-H2S vapour at elevated temperature and pressure

This experimental study sheds light on the complexation of gold in reduced sulphur-bearing vapour, specifically, in H₂O-H₂S gas mixtures. The solubility of gold was determined in experiments at temperatures of 300, 350 and 365 °C and reached 2.2, 6.6 and 6.3 μg/kg, respectively. The density of the vapour varied from 0.02 to 0.22 g/cm³, the mole fraction of H₂S varied from 0.03 to 0.96, and the pressure in the cell reached 263 bar. Statistically significant correlations of the amount of gold dissolved in the fluid with the fugacity of H₂O and H₂S permit the experimental data to be fitted to a solvation/hydration model. According to this model, the solubility of gold in H₂O-H₂S gas mixtures is controlled by the formation of sulphide or bisulphide species solvated by H₂S or H₂O molecules. Formation of gold sulphide species is favoured statistically over gold bisulphide species and thus the gold is interpreted to dissolve according to reactions of the form:

(A1)     

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

高温高压下NaCl-H2O体系的Raman光谱分析: 一种新的流体包裹体盐度测定方法

利用水热金刚石压腔研究了高温高压下0.0%、5.0%、10.0%和20.0%的NaCl-H2O体系的Raman光谱,研究结果发现:对于盐度一定的NaCl-H2O体系,当压力不变时,水的Raman伸缩谱峰的波数变化量△v(O-H)随温度的升高而不断增加;当压力不同时,△v(O-H)随温度的变化关系(偏移斜率S)不受压力的影响,而只与盐度相关.利用该性质可以用来确定流体包裹体的盐度(W),其计算公式如下:W=123.25S2-128.11S+32.40,误差为0.46%.     

Nature of polymerization and properties of silicate melts and glasses at high pressure

The structures of sodium silicate and aluminosilicate glasses quenched from melts at high pressure (6-10 GPa) with varying degrees of polymerization (fractions of nonbridging oxygen) were explored using solid-state NMR [¹⁷O and ²⁷Al triple-quantum magic-angle spinning (3QMAS) NMR]. The bond connectivity in melts among four and highly coordinated network polyhedra, such as ^[4]Al, ^[5,6]Al, ^[4]Si, and ^[5,6]Si, at high pressure is shown to be significantly different from that at ambient pressure. In particular, in the silicate and aluminosilicate melts, the proportion of nonbridging oxygen (NBO) generally decreases with increasing pressure, leading to the formation of new oxygen clusters that include 5- and 6-coordinated Si and Al in addition to 4-coordinated Al and Si, such as ^[4]Si-O-^[5,6]Si, ^[4]Si-O-^[5,6]Al and Na-O-^[5,6]Si. While the fractions of ^[5,6]Al increase with pressure, the magnitude of this increase diminishes with increasing degrees of ambient-pressure polymerization under isobaric conditions. Incorporating the above structural information into models of melt properties reproduces the anomalous pressure-dependence of O²⁻ diffusivity and viscosity often observed in silicate melts.     

Review of the systematics of CO2–H2O fluid inclusions

Aqueous solutions that contain volatile (gas) components are one of the most important types of fluid in the Earth's crust. The record that such fluids have left in the form of fluid inclusions in minerals provides a wealth of insight into the geochemical and petrologic processes in which the fluids participated. This article reviews the systematics of CO₂–H₂O fluid inclusions as a starting point for interpreting the chemically more complex systems. The phase relations of the binary are described with respect to a qualitative P–T–X model, and isoplethic–isochoric paths through this model are used to explain the equilibrium and non-equilibrium behaviour of fluid inclusions during microthermometric heating and cooling. The P–T–X framework is then used to discuss the various modes of fluid inclusion entrapment, and how the resulting assemblage textures can be used to interpret the P–T conditions, phase states, and evolution paths of the parent solutions. Finally, quantitative methods are reviewed by which bulk molar volume and composition of CO₂–H₂O fluid inclusions can be determined from microthermometric observations of phase transitions.     

Kinetics of H2-O2-H2O redox equilibria and formation of metastable H2O2 under low temperature hydrothermal conditions

Hydrothermal experiments were conducted to evaluate the kinetics of H_2(aq) oxidation in the homogeneous H₂-O₂-H₂O system at conditions reflecting subsurface/near-seafloor hydrothermal environments (55-250 °C and 242-497 bar). The kinetics of the water-forming reaction that controls the fundamental equilibrium between dissolved H_2(aq) and O_2(aq), are expected to impose significant constraints on the redox gradients that develop when mixing occurs between oxygenated seawater and high-temperature anoxic vent fluid at near-seafloor conditions. Experimental data indicate that, indeed, the kinetics of H_2(aq)-O_2(aq) equilibrium become slower with decreasing temperature, allowing excess H_2(aq) to remain in solution. Sluggish reaction rates of H_2(aq) oxidation suggest that active microbial populations in near-seafloor and subsurface environments could potentially utilize both H_2(aq) and O_2(aq), even at temperatures lower than 40 °C due to H_2(aq) persistence in the seawater/vent fluid mixtures. For these H₂-O₂ disequilibrium conditions, redox gradients along the seawater/hydrothermal fluid mixing interface are not sharp and microbially-mediated H_2(aq) oxidation coupled with a lack of other electron acceptors (e.g. nitrate) could provide an important energy source available at low-temperature diffuse flow vent sites.More importantly, when H_2(aq)-O_2(aq) disequilibrium conditions apply, formation of metastable hydrogen peroxide is observed. The yield of H₂O_2(aq) synthesis appears to be enhanced under conditions of elevated H_2(aq)/O_2(aq) molar ratios that correspond to abundant H_2(aq) concentrations. Formation of metastable H₂O₂ is expected to affect the distribution of dissolved organic carbon (DOC) owing to the existence of an additional strong oxidizing agent. Oxidation of magnetite and/or Fe⁺⁺ by hydrogen peroxide could also induce formation of metastable hydroxyl radicals (•OH) through Fenton-type reactions, further broadening the implications of hydrogen peroxide in hydrothermal environments.     

Molecular H2O in armenite, BaCa2Al6Si9O30·2H2O, and epididymite, Na2Be2Si6O15·H2O: Heat capacity, entropy and local-bonding behavior of confined H2O in microporous silicates

Armenite, ideal formula BaCa₂Al₆Si₉O₃₀·2H₂O, and its dehydrated analog BaCa₂Al₆Si₉O₃₀ and epididymite, ideal formula Na₂Be₂Si₆O₁₅·H₂O, and its dehydrated analog Na₂Be₂Si₆O₁₅ were studied by low-temperature relaxation calorimetry between 5 and 300 K to determine the heat capacity, C_p, behavior of their confined H₂O. Differential thermal analysis and thermogravimetry measurements, FTIR spectroscopy, electron microprobe analysis and powder Rietveld refinements were undertaken to characterize the phases and the local environment around the H₂O molecule.The determined structural formula for armenite is Ba_0.88(0.01)Ca_1.99(0.02)Na_0.04(0.01)Al_5.89(0.03)Si_9.12(0.02)O₃₀·2H₂O and for epididymite Na_1.88(0.03)K_0.05(0.004)Na_0.01(0.004)Be_2.02(0.008)Si_6.00(0.01)O₁₅·H₂O. The infrared (IR) spectra give information on the nature of the H₂O molecules in the natural phases via their H₂O stretching and bending vibrations, which in the case of epididymite only could be assigned. The powder X-ray diffraction data show that armenite and its dehydrated analog have similar structures, whereas in the case of epididymite there are structural differences between the natural and dehydrated phases. This is also reflected in the lattice IR mode behavior, as observed for the natural phases and the H₂O-free phases. The standard entropy at 298 K for armenite is S° = 795.7 ± 6.2 J/mol K and its dehydrated analog is S° = 737.0 ± 6.2 J/mol K. For epididymite S° = 425.7 ± 4.1 J/mol K was obtained and its dehydrated analog has S° = 372.5 ± 5.0 J/mol K. The heat capacity and entropy of dehydration at 298 K are Δ = 3.4 J/mol K and ΔS^rxn = 319.1 J/mol K and Δ = −14.3 J/mol K and ΔS^rxn = 135.7 J/mol K for armenite and epididymite, respectively. The H₂O molecules in both phases appear to be ordered. They are held in place via an ion-dipole interaction between the H₂O molecule and a Ca cation in the case of armenite and a Na cation in epididymite and through hydrogen-bonding between the H₂O molecule and oxygen atoms of the respective silicate frameworks. Of the three different H₂O phases ice, liquid water and steam, the C_p behavior of confined H₂O in both armenite and epididymite is most similar to that of ice, but there are differences between the two silicates and from the C_p behavior of ice. Hydrogen-bonding behavior and its relation to the entropy of confined H₂O at 298 K is analyzed for various microporous silicates.The entropy of confined H₂O at 298 K in various silicates increases approximately linearly with increasing average wavenumber of the OH-stretching vibrations. The interpretation is that decreased hydrogen-bonding strength between a H₂O molecule and the silicate framework, as well as weak ion-dipole interactions, results in increased entropy of H₂O. This results in increased amplitudes of external H₂O vibrations, especially translations of the molecule, and they contribute strongly to the entropy of confined H₂O at T < 298 K.     

CO2 solubility and speciation in intermediate (andesitic) melts: the role of H2O and composition

We determined total CO₂ solubilities in andesite melts with a range of compositions. Melts were equilibrated with excess C-O(-H) fluid at 1 GPa and 1300°C then quenched to glasses. Samples were analyzed using an electron microprobe for major elements, ion microprobe for C-O-H volatiles, and Fourier transform infrared spectroscopy for molecular H₂O, OH⁻, molecular CO₂, and CO₃²⁻. CO₂ solubility was determined in hydrous andesite glasses and we found that H₂O content has a strong influence on C-O speciation and total CO₂ solubility. In anhydrous andesite melts with ∼60 wt.% SiO₂, total CO₂ solubility is ∼0.3 wt.% at 1300°C and 1 GPa and total CO₂ solubility increases by about 0.06 wt.% per wt.% of total H₂O. As total H₂O increases from ∼0 to ∼3.4 wt.%, molecular CO₂ decreases (from 0.07 ± 0.01 wt.% to ∼0.01 wt.%) and CO₃²⁻ increases (from 0.24 ± 0.04 wt.% to 0.57 ± 0.09 wt.%). Molecular CO₂ increases as the calculated mole fraction of CO₂ in the fluid increases, showing Henrian behavior. In contrast, CO₃²⁻ decreases as the calculated mole fraction of CO₂ in the fluid increases, indicating that CO₃²⁻ solubility is strongly dependent on the availability of reactive oxygens in the melt. These findings have implications for CO₂ degassing. If substantial H₂O is present, total CO₂ solubility is higher and CO₂ will degas at relatively shallow levels compared to a drier melt. Total CO₂ solubility was also examined in andesitic glasses with additional Ca, K, or Mg and low H₂O contents (<1 wt.%). We found that total CO₂ solubility is negatively correlated with (Si + Al) cation mole fraction and positively correlated with cations with large Gibbs free energy of decarbonation or high charge-to-radius ratios (e.g., Ca). Combining our andesite data with data from the literature, we find that molecular CO₂ is more abundant in highly polymerized melts with high ionic porosities (>∼48.3%), and low nonbridging oxygen/tetrahedral oxygen (<∼0.3). Carbonate dominates most silicate melts and is most abundant in depolymerized melts with low ionic porosities, high nonbridging oxygen/tetrahedral oxygen (>∼0.3), and abundant cations with large Gibbs free energy of decarbonation or high charge-to-radius ratio. In natural silicate melt, the oxygens in the carbonate are likely associated with tetrahedral and network-modifying cations (including Ca, H, or H-bonds) or a combinations of those cations.     

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

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

Solution mechanisms of H2O in depolymerized peralkaline melts

Solubility mechanisms of water in depolymerized silicate melts quenched from high temperature (1000°-1300°C) at high pressure (0.8-2.0 GPa) have been examined in peralkaline melts in the system Na₂O-SiO₂-H₂O with Raman and NMR spectroscopy. The Na/Si ratio of the melts ranged from 0.25 to 1. Water contents were varied from ∼3 mol% and ∼40 mol% (based on O = 1). Solution of water results in melt depolymerization where the rate of depolymerization with water content, ∂(NBO/Si)/∂X_H2O, decreases with increasing total water content. At low water contents, the influence of H₂O on the melt structure resembles that of adding alkali oxide. In water-rich melts, alkali oxides are more efficient melt depolymerizers than water. In highly polymerized melts, Si-OH bonds are formed by water reacting with bridging oxygen in Q⁴-species to form Q³ and Q² species. In less polymerized melts, Si-OH bonds are formed when bridging oxygen in Q³-species react with water to form Q²-species. In addition, the presence of Na-OH complexes is inferred. Their importance appears to increase with Na/Si. This apparent increase in importance of Na-OH complexes with increasing Na/Si (which causes increasing degree of depolymerization of the anhydrous silicate melt) suggests that water is a less efficient depolymerizer of silicate melts, the more depolymerized the melt. This conclusion is consistent with recently published ¹H and ²⁹Si MAS NMR and ¹H-²⁹Si cross polarization NMR data.     

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.     

Solution behavior of reduced COH volatiles in silicate melts at high pressure and temperature

The solubility and solution mechanisms of reduced COH volatiles in Na₂OSiO₂ melts in equilibrium with a (H₂ + CH₄) fluid at the hydrogen fugacity defined by the iron-wüstite + H₂O buffer [f_H2(IW)] have been determined as a function of pressure (1-2.5 GPa) and silicate melt polymerization (NBO/Si: nonbridging oxygen per silicon) at 1400 °C. The solubility, calculated as CH₄, increases from ∼0.2 wt% to ∼0.5 wt% in the melt NBO/Si-range ∼0.4 to ∼1.0. The solubility is not significantly pressure-dependent, probably because f_H2(IW) in the 1-2.5 GPa range does not vary greatly with pressure. Carbon isotope fractionation between methane-saturated melts and (H₂ + CH₄) fluid varied by ∼14‰ in the NBO/Si-range of these melts.The (C..H) and (O..H) speciation in the quenched melts was determined with Raman and ¹H MAS NMR spectroscopy. The dominant (C..H)-bearing complexes are molecular methane, CH₄, and a complex or functional group that includes entities with CCH bonding. Minor abundance of complexes that include SiOCH₃ bonding is tentatively identified in some melts. There is no spectroscopic evidence for SiC or SiCH₃. Raman spectra indicate silicate melt depolymerization (increasing NBO/Si). The [CH₄/CCH]^melt abundance ratio is positively correlated with NBO/Si, which is interpreted to suggest that the (CCH)-containing structural entity is bonded to the silicate melt network structure via its nonbridging oxygen. The ∼14‰ carbon isotope fractionation change between fluid and melt is because of the speciation changes of carbon in the melt.     

The solution behavior of H2O in peralkaline aluminosilicate melts at high pressure with implications for properties of hydrous melts

Solubility and solution mechanisms of H₂O in depolymerized melts in the system Na₂O-Al₂O₃-SiO₂ were deduced from spectroscopic data of glasses quenched from melts at 1100 °C at 0.8-2.0 GPa. Data were obtained along a join with fixed nominal NBO/T = 0.5 of the anhydrous materials [Na₂Si₄O₉-Na₂(NaAl)₄O₉] with Al/(Al+Si) = 0.00-0.25. The H₂O solubility was fitted to the expression, X_H2O=0.20+0.0020f_H2O-0.7X_Al+0.9(X_Al)², where X_H2O is the mole fraction of H₂O (calculated with O = 1), f_H2O the fugacity of H₂O, and X_Al = Al/(Al+Si). Partial molar volume of H₂O in the melts, , calculated from the H₂O-solulbility data assuming ideal mixing of melt-H₂O solutions, is 12.5 cm³/mol for Al-free melts and decreases linearly to 8.9 cm³/mol for melts with Al/(Al+Si) ∼ 0.25. However, if recent suggestion that is composition-independent is applied to constrain activity-composition relations of the hydrous melts, the activity coefficient of H₂O, , increases with Al/(Al+Si).Solution mechanisms of H₂O were obtained by combining Raman and ²⁹Si NMR spectroscopic data. Degree of melt depolymerization, NBO/T, increases with H₂O content. The rate of NBO/T-change with H₂O is negatively correlated with H₂O and positively correlated with Al/(Al+Si). The main depolymerization reaction involves breakage of oxygen bridges in Q⁴-species to form Q² species. Steric hindrance appears to restrict bonding of H⁺ with nonbridging oxygen in Q³ species. The presence of Al³⁺ does not affect the water solution mechanisms significantly.     

H2O diffusion in dacitic and andesitic melts

The diffusion of water in dacitic and andesitic melts was investigated at temperatures of 1458 to 1858 K and pressures between 0.5 and 1.5 GPa using the diffusion couple technique. Pairs of nominally dry glasses and hydrous glasses containing between 1.5 and 6.3 wt.% dissolved H₂O were heated for 60 to 480 s in a piston cylinder apparatus. Concentration profiles of hydrous species (OH groups and H₂O molecules) and total water (C_H2Ot = sum of OH and H₂O) were measured along the cylindrical axis of the diffusion sample using IR microspectroscopy. Electron microprobe traverses show no significant change in relative proportions of anhydrous components along H₂O profiles, indicating that our data can be treated as effective binary interdiffusion between H₂O and the rest of the silicate melt. Bulk water diffusivity (D_H2Ot) was derived from profiles of total water using a modified Boltzmann-Matano method as well as using fittings assuming a functional relationship between D_H2Ot and C_H2Ot. In dacitic melts D_H2Ot is proportional to C_H2Ot up to 6 wt.%. In andesitic melts the dependence of D_H2Ot on C_H2Ot is less pronounced. A pressure effect on water diffusivity could not be resolved for either dacitic or andesitic melt in the range 0.5 to 1.5 GPa. Combining our results with previous studies on water diffusion in rhyolite and basalt show that for a given water content D_H2Ot increases monotonically with increasing melt depolymerization at temperatures >1500 K. Assuming an Arrhenian behavior in the whole compositional range, the following formulation was derived to estimate D_H2Ot (m²/s) at 1 wt.% H₂O_t in melts with rhyolitic to andesitic composition as a function of T (K), P (MPa) and S (wt.% SiO₂):

     

Solubility behavior of alkali aluminosilicate components in aqueous fluids and silicate melts at high pressure and temperature

The solubility behavior of K₂O, Na₂O, Al₂O₃, and SiO₂ in silicate-saturated aqueous fluid and coexisting H₂O-saturated silicate melts in the systems K₂O-Al₂O₃-SiO₂-H₂O and Na₂O-Al₂O₃-SiO₂-H₂O has been examined in the 1- to 2-GPa pressure range at 1100°C. Glasses of Na- and K-tetrasilicate compositions with 0, 3, and 6 mol% Al₂O₃ were used as starting materials. In both systems, the oxides dissolve incongruently in aqueous fluid and silicate melt. When recalculated to an anhydrous basis, the aqueous fluids are enriched in alkalis and depleted in silica and alumina relative to their proportions in the starting materials. The extent of incongruency is more pronounced in the Na₂O-Al₂O₃-SiO₂-H₂O system than in the K₂O-Al₂O₃-SiO₂-H₂O system.The partition coefficients of the oxides, D_oxide^fluid/melt, are linear and positive functions of the oxide concentration in the fluid for each composition. There is a slight dependence of the partition coefficients on bulk composition. No effect of pressure could be discerned. For alkali metals, the fluid/melt partition coefficients range from 0.06 to 0.8. For Al₂O₃ this range is 0.01 to 0.2, and for SiO₂, it is 0.01 to 0.32. For all compositions, D_K2O^fluid/melt∼D_Na2O^fluid/melt>D_SiO2^fluid/melt>D_Al2O3^fluid/melt for the same oxide concentration in the fluid. D_K2O^fluid/melt, D_Na2O^fluid/melt, and D_SiO2^fluid/melt correlate negatively with the Al₂O₃ content of the systems. This correlation is consistent with a solubility model of alkalis that involve associated KOH°, NaOH°, silicate, and aluminate complexes.     

CO2 solubility and solubility mechanisms in silicate melts at high pressures

CO₂ solubility has a slight negative temperature dependence in olivine melilitite at 30 kb with 9% CO₂ dissolved at 1,450 °C, 8.5% at 1,550 °C and 8.3% at 1,650° C. CO₂ is dissolved as the carbonate molecule (CO ₃ ^2– ) only. Feldspar melts (albite-anorthite) dissolve much less CO₂ at 30 kb (around 2%) with a slight increase with increasing anorthite content. A CO₂ absorption peak in infrared spectra of albite-rich glasses diappears in favour of the CO ₃ ^2– peak with increasing anorthite content. It is inferred that CO₂ was present as CO ₃ ^2– in albite-rich melts also, but reverts to CO₂ during quenching because of bonding differences related to Ca²⁺ and Na⁺ in the melts.     

Fractionation of major elements between coexisting H2O-saturated silicate melt and silicate-saturated aqueous fluids in aluminosilicate systems at 1-2 GPa

From experimental data in the systems Na₂O-Al₂O₃-SiO₂-H₂O, K₂O-Al₂O₃-SiO₂-H₂O at 1100°C, and CaO-Al₂O₃-SiO₂-H₂O at 1200°C in the 1-2 GPa pressure range, the solution behavior of the individual oxides in coexisting H₂O-saturated silicate melts and silicate-saturated aqueous fluids appears to be incongruent. Recalculated on an anhydrous basis, in the CaO-Al₂O₃-SiO₂-H₂O system, CaO^fluid/CaO^melt < 1, whereas in the Na₂O-Al₂O₃-SiO₂-H₂O and K₂O-Al₂O₃-SiO₂-H₂O systems, K₂O^fluid/K₂O^melt and Na₂O^fluid/Na₂O^melt both are greater than 1. The aqueous fluids are depleted in alumina relative to silicate melt.In the Na₂O-Al₂O₃-SiO₂-H₂O, K₂O-Al₂O₃-SiO₂-H₂O, and CaO-Al₂O₃-SiO₂-H₂O systems, fluid/melt partition coefficients for the individual oxides range between ∼0.005 and 0.35 depending on oxide, bulk composition and pressure. The alkali partition coefficients are about an order of magnitude higher than that of CaO. Alumina and silica partition coefficient values in the CaO-Al₂O₃-SiO₂-H₂O system are 10-20% of the values for the same oxides in the Na₂O-Al₂O₃-SiO₂-H₂O and K₂O-Al₂O₃-SiO₂-H₂O systems.Positive correlations among individual partition coefficients and oxide concentrations in the aqueous fluids are consistent with complexing in the fluid that involves silicate polymers associated with alkalis and alkaline earths and aluminosilicate complexes where alkalis and alkaline earths may serve to charge-balance Al³⁺, which is, perhaps, in tetrahedral coordination. Alkali aluminosilicate complexes in aqueous fluid appear more stable than Ca-aluminosilicate complexes.     

H2O diffusion models in rhyolitic melt with new high pressure data

Water diffusion in silicate melts is important for understanding bubble growth in magma, magma degassing and eruption dynamics of volcanos. Previous studies have made significant progress on water diffusion in silicate melts, especially rhyolitic melt. However, the pressure dependence of H₂O diffusion is not constrained satisfactorily. We investigated H₂O diffusion in rhyolitic melt at 0.95–1.9 GPa and 407–1629 °C, and 0.2–5.2 wt.% total water (H₂O_t) content with the diffusion-couple method in a piston-cylinder apparatus. Compared to previous data at 0.1–500 MPa, H₂O diffusivity is smaller at higher pressures, indicating a negative pressure effect. This pressure effect is more pronounced at low temperatures. Assuming H₂O diffusion in rhyolitic melt is controlled by the mobility of molecular H₂O (H₂O_m), the diffusivity of H₂O_m (D_H2Om) at H₂O_t ≤ 7.7 wt.%, 403–1629 °C, and ≤ 1.9 GPa is given by

D_H2Om=D₀exp(aX),