Effect of oxygen fugacity on the H2O storage capacity of forsterite in the carbon-saturated systems

High pressure experiments have been performed in the systems Mg₂SiO₄-C-O-H and Mg₂SiO₄-K₂CO₃-C at 6.3 GPa and 1200 to 1600 °C using a split-sphere multi-anvil apparatus. In the Mg₂SiO₄-C-O-H system the composition of fluid was modeled by adding different amounts of water and stearic acid. The fO₂ was controlled by the Mo-MoO₂ or Fe-FeO oxygen buffers. Several experiments in the Mg₂SiO₄-C-O-H system and all experiments in the Mg₂SiO₄-K₂CO₃-C system have been conducted without buffering the fO₂. Forsterite in the system Mg₂SiO₄-K₂CO₃-C does not reveal OH absorption bands in the IR spectra, while forsterite coexisting with carbon-bearing fluid and silicate melt at logfO₂ from FMQ-2 to FMQ-5 (from 2 to 5 log units below fayalite-magnetite-quartz oxygen buffer) contains 800-1850 wt. ppm H₂O. The maximum concentrations were detected at 1400 °C and FMQ-3.5. We observed an increase in the solidus temperature in the system Mg₂SiO₄-C-O-H from 1200 to above 1600 °C with log fO₂ decreasing from FMQ-2 to FMQ-5. The increase of the solidus temperature and the broadening of the stability field of the H₂O-H₂-CH₄ subsolidus fluid phase at 1400-1600 °C explain the high H₂O storage capacity of forsterite relative to that crystallized from carbon-free, oxidized, hydrous, silicic melt. At temperatures above 1400 °C liquidus forsterite precipitated along with diamond from oxidized (FMQ-1) carbonate-silicate melt and from silicate melt dissolving the moderately reduced C-O-H fluid (from FMQ-2 to FMQ-3.5). Formation of diamond was not detected under ultra-reduced conditions (FMQ-5) at 1200-1600 °C. Olivine co-precipitating with diamond from dry carbonate-silicate or hydrous-silicic fluid/melt can provide information on the H₂O contents and speciation of the diamond-forming media in the mantle. The conditions for minimum post-crystallization alteration of olivine and its hydrogen content are discussed.     

Thermodynamics, structure, dynamics, and freezing of Mg2SiO4 liquid at high pressure

We perform first principles molecular dynamics simulations of Mg₂SiO₄ liquid and crystalline forsterite. On compression by a factor of two, we find that the Grüneisen parameter of the liquid increases linearly from 0.6 to 1.2. Comparison of liquid and forsterite equations of state reveals a temperature-dependent density crossover at pressures of ∼12-17 GPa. Along the melting curve, which we calculate by integration of the Clapeyron equation, the density crossover occurs within the forsterite stability field at P = 13 GPa and T = 2550 K. The melting curve obtained from the root mean-square atomic displacement in forsterite using the Lindemann law fails to match experimental or calculated melting curves. We attribute this failure to the liquid structure that differs significantly from that of forsterite, and which changes markedly upon compression, with increases in the degree of polymerization and coordination. The mean Si coordination increases from 4 in the uncompressed system to 6 upon twofold compression. The self-diffusion coefficients increase with temperature and decrease monotonically with pressure, and are well described by the Arrhenian relation. We compare our equation of state to the available highpressure shock wave data for forsterite and wadsleyite. Our theoretical liquid Hugoniot is consistent with partial melting along the forsterite Hugoniot at pressures 150-170 GPa, and complete melting at 170 GPa. The wadsleyite Hugoniot is likely sub-liquidus at the highest experimental pressure to date (200 GPa).     

Pressure dependence of the speciation of dissolved water in rhyolitic melts

Water speciation in rhyolitic melts with dissolved water ranging from 0.8 to 4 wt% under high pressure was investigated. Samples were heated in a piston-cylinder apparatus at 624-1027 K and 0.94-2.83 GPa for sufficient time to equilibrate hydrous species (molecular H₂O and hydroxyl group, H₂O_m + O ? 2OH) in the melts and then quenched roughly isobarically. The concentrations of both hydrous species in the quenched glasses were measured with Fourier transform infrared (FTIR) spectroscopy. For the samples with total water content less than 2.7 wt%, the equilibrium constant (K) is independent of total H₂O concentration. Incorporating samples with higher water contents, the equilibrium constant depends on total H₂O content, and a regular solution model is used to describe the dependence. K changes with pressure nonmonotonically for samples with a given water content at a given temperature. The equilibrium constant does not change much from ambient pressure to 1 GPa, but it increases significantly from 1 to 3 GPa. In other words, more molecular H₂O reacts to form hydroxyl groups as pressure increases from 1 GPa, which is consistent with breakage of tetrahedral aluminosilicate units due to compression of the melt induced by high pressure. The effect of 1.9 GPa (from 0.94 to 2.83 GPa) on the equilibrium constant at 873 K is equivalent to a temperature effect of 49 K (from 873 K to 922 K) at 0.94 GPa. The results can be used to evaluate the role of speciation in water diffusion, to estimate the apparent equilibrium temperature, and to infer viscosity of hydrous rhyolitic melts under high pressure.     

High pressure fluids in the system MgO–SiO2–H2O under upper mantle conditions

Fluids and melts have been trapped and analysed in high pressure experiments in the model mantle system MgO-SiO₂-H₂O at 6 to 10.5 GPa and 900 to 1,200 °C. The fluid/melt traps consisted of a diamond layer that was added to the experimental charge and was separate from the silicate phases. The recovered diamond traps were analysed by laser ablation - ICP - MS. Starting materials were synthetic mixtures of brucite, talc and silica with variable Mg/Si containing 11-31 wt% H₂O. Experiments on a serpentine starting composition [Mg₃Si₂O₅(OH)₄] result in MgO/SiO₂ weight ratios in the subsolidus fluids close to 1 at 6 GPa and close to 2 at 9 GPa. Melt compositions at 6 and 9 GPa have MgO/SiO₂ ratios close to that of forsterite. At a single pressure the amount of dissolved silicate in the fluid increases steadily with increasing temperature up to 1,150 °C, where a sudden increase of both SiO₂ and MgO is observed. This discrete step marks the solidus, which is more clearly developed at 6 than at 9 GPa. Thus, hydrous melts within the model mantle subsystem Mg₂SiO₄-Mg₂Si₂O₆-H₂O are chemically distinct from aqueous fluids up to at least 9 GPa, corresponding to 300 km depth. Extrapolation of the current data set implies that total convergence between fluid and melt along the solidus probably occurs at 12-13 GPa (~400 km), i.e. close to the Earth's mantle transition zone. Beneath cratons, interactions of hydrous fluids with upper mantle lithologies cause relative silica depletion (olivine enrichment) at depths greater than 200 km and silica (orthopyroxene) enrichment at shallower depths.     

Acoustic velocity measurements on Na2O-TiO2-SiO2 liquids: Evidence for a highly compressible TiO2 component related to five-coordinated Ti

Longitudinal acoustic velocities were measured at 1 bar in 10 Na₂O-TiO₂-SiO₂ (NTS) liquids for which previous density and thermal expansion data are reported in the literature. Data were collected with a frequency-sweep acoustic interferometer at centered frequencies of 4.5, 5, and 6 MHz between 1233 and 1896 K; in all cases, the sound speeds decrease with increasing temperature. Six of the liquids have a similar TiO₂ concentration (∼25 mol %), so that the effect of varying Na/Si ratio on the partial molar compressibility of the TiO₂ component can be evaluated. Theoretically based models for β_T and (∂V/∂P)_T as a function of composition and temperature are presented. As found previously for the partial molar volume of TiO₂ in sodium silicate melts, values of (13.7-18.8 × 10⁻²/GPa) vary systematically with the Na/Si and Na/(Si + Ti) ratio in the liquid. In contrast values of for the SiO₂ and Na₂O components (6.6 and 8.0 × 10⁻²/GPa, respectively, at 1573 K) are independent of composition. Na₂O is the only component that contributes to the temperature dependence of the compressibility of NTS liquids (1.13 ± 0.04 × 10⁻⁴/GPa K). The results further indicate that the TiO₂ component is twice as compressible as the Na₂O and SiO₂ components. The enhanced compressibility of TiO₂ appears to be related to the abundance of five-coordinated Ti (^[5]Ti) in these liquids, but not with a change in Ti coordination. Instead, it is proposed that the asymmetric geometry of ^[5]Ti in a square pyramidal site promotes different topological rearrangements in alkali titanosilicate liquids, which lead to the enhanced compressibility of TiO₂.     

Equation of state of the H2O, CO2, and H2O-CO2 systems up to 10 GPa and 2573.15 K: Molecular dynamics simulations with ab initio potential surface

Based on our previous development of the molecular interaction potential for pure H₂O and CO₂ [Zhang, Z.G., Duan, Z.H. 2005a. Isothermal-isobaric molecular dynamics simulations of the PVT properties of water over wide range of temperatures and pressures. Phys. Earth Planet Interiors149, 335-354; Zhang, Z.G., Duan, Z.H. 2005b. An optimized molecular potential for carbon dioxide. J. Chem. Phys.122, 214507] and the ab initio potential surface across CO₂-H₂O molecules constructed in this study, we carried out more than one thousand molecular dynamics simulations of the PVTx properties of the CO₂-H₂O mixtures in the temperature-pressure range from 673.15 to 2573.15 K up to 10.0 GPa. Comparison with extensive experimental PVTx data indicates that the simulated results generally agree with experimental data within 2% in density, equivalent to experimental uncertainty. Even the data under the highest experimental temperature-pressure conditions (up to 1673 K and 1.94 GPa) are well predicted with the agreement within 1.0% in density, indicating that the high accuracy of the simulation is well retained as the temperature and pressure increase. The consistent and stable predictability of the simulation from low to high temperature-pressure and the fact that the molecular dynamics simulation resort to no experimental data but to ab initio molecular potential makes us convinced that the simulation results should be reliable up to at least 2573 K and 10 GPa with errors less than 2% in density. In order to integrate all the simulation results of this study and previous studies [Zhang and Duan, 2005a, 2005b] and the experimental data for the calculation of volumetric properties (volume, density, and excess volume), heat properties, and chemical properties (fugacity, activity, and possibly supercritical phase separation), an equation of state (EOS) is laboriously developed for the CO₂, H₂O, and CO₂-H₂O systems. This EOS reproduces all the experimental and simulated data covering a wide temperature and pressure range from 673.15 to 2573.15 K and from 0 to 10.0 GPa within experimental or simulation uncertainty.     

Stability of clinohumite in the system MgO-SiO2-H2O

The low-pressure stability of clinohumite has been investigated in phase-equilibrium experiments on the reaction forsterite + brucite = clinohumite. The reaction was bracketed between 2.45 and 2.84 GPa at 650 °C, extending to between 1.37 and 1.57 GPa at 850 °C. At temperatures above the reaction brucite = periclase + vapour, the reaction clinohumite = forsterite + vapour was bracketed between 1.27 and 1.52 GPa at 900 °C, rising to between 1.90 and 2.00 GPa at 1000 °C. The position of the reaction forsterite + brucite = clinohumite is ∼0.5 GPa below the position determined in previous work, the difference arising either from pressure uncertainties in both studies, from enhanced reaction to clinohumite in this study due to the presence of excess brucite in the starting material, or from different concentrations of defects in the two samples. The brackets on the reaction were combined with other measured and estimated thermodynamic data for clinohumite to determine its enthalpy of formation and entropy, in a revised version of the thermodynamic dataset of Holland and Powell (1998). The values obtained were ΔH _f =−9607.29±3.05 kJ mol⁻¹, S=445 J mol⁻¹ K⁻¹. These data were used to calculate positions of other reactions involving clinohumite. The calculations suggest a larger stability field for clinohumite than implied by the results of previous experimental studies, indicating a need for more high-pressure phase-equilibrium studies to provide better thermodynamic data. Received: 30 April 1999 / Accepted: 8 November 1999     

Solubility of niobium in the system CaCO<Subscript>3</Subscript>–CaF<Subscript>2</Subscript>–NaNbO<Subscript>3</Subscript> at 0.1 GPa pressure: implications for the crystallization of pyrochlore from carbonatite magma

Experimental data are presented for the solubility of NaNbO₃ in the ternary system CaCO₃–CaF₂–NaNbO₃ (or calcite–fluorite–lueshite) over the temperature range 500–1,000°C at 0.1 GPa pressure. Liquidus to solidus phase relationships are given for the pseudo-binary join ([CaCO₃]₆₀[CaF₂]₄₀)_100-x–(NaNbO₃)_x (0<x<60 wt%). These data show that the maximum solubility of NaNbO₃ in these liquids is about 17 wt% (or 13.8 wt% Nb₂O₅) at approximately 930°C, and is represented by the appearance of pyrochlore as the primary liquidus phase. The sub-liquidus assemblages with decreasing temperature for NaNbO₃ contents of 20–50 wt% are: pyrochlore + liquid; pyrochlore + CaF₂ + liquid; pyrochlore + CaF₂ + CaCO₃ + liquid. The solidus assemblage is pyrochlore + CaF₂ + CaCO₃ at temperatures of approximately 700°C (20 wt% NaNbO₃) and 600°C (40 wt% NaNbO₃). NaNbO₃ is present only in sub-solidus assemblages. These data show that in this fluorine-bearing anhydrous system pyrochlore is the principal Nb-hosting supra-solidus phase, in contrast to fluorine-free hydrous melts from which perovskite-structured compounds crystallize. The crystallization of pyrochlore and/or perovskite-structured compounds from haplocarbonatite liquids is thus considered to be dependent upon the F/OH ratio of the melt.     

Effect of KCl on melting in the Mg<Subscript>2</Subscript>SiO<Subscript>4</Subscript>–MgSiO<Subscript>3</Subscript>–H<Subscript>2</Subscript>O system at 5 GPa

To examine the effect of KCl-bearing fluids on the melting behavior of the Earth’s mantle, we conducted experiments in the Mg₂SiO₄–MgSiO₃–H₂O and Mg₂SiO₄–MgSiO₃–KCl–H₂O systems at 5 GPa. In the Mg₂SiO₄–MgSiO₃–H₂O system, the temperature of the fluid-saturated solidus is bracketed between 1,200–1,250°C, and both forsterite and enstatite coexist with the liquid under supersolidus conditions. In the Mg₂SiO₄–MgSiO₃–KCl–H₂O systems with molar Cl/(Cl + H₂O) ratios of 0.2, 0.4, and 0.6, the temperatures of the fluid-saturated solidus are bracketed between 1,400–1,450°C, 1,550–1,600°C, and 1,600–1,650°C, respectively, and only forsterite coexists with liquid under supersolidus conditions. This increase in the temperature of the solidus demonstrates the significant effect of KCl on reducing the activity of H₂O in the fluid in the Mg₂SiO₄–MgSiO₃–H₂O system. The change in the melting residues indicates that the incongruent melting of enstatite (enstatite = forsterite + silica-rich melt) could extend to pressures above 5 GPa in KCl-bearing systems, in contrast to the behavior in the KCl-free system.     

Pressure dependence of viscosity of rhyolitic melts

Viscosity of silicate melts is a critical property for understanding volcanic and igneous processes in the Earth. We investigate the pressure effect on the viscosity of rhyolitic melts using two methods: indirect viscosity inference from hydrous species reaction in melts using a piston cylinder at pressures up to 2.8 GPa and direct viscosity measurement by parallel-plate creep viscometer in an internally-heated pressure vessel at pressures up to 0.4 GPa. Comparison of viscosities of a rhyolitic melt with 0.8 wt% water at 0.4 GPa shows that both methods give consistent results. In the indirect method, viscosities of hydrous rhyolitic melts were inferred based on the kinetics of hydrous species reaction in the melt upon cooling (i.e., the equivalence of rheologically defined glass transition temperature and chemically defined apparent equilibrium temperature). The cooling experiments were carried out in a piston-cylinder apparatus using hydrous rhyolitic samples with 0.8-4 wt% water. Cooling rates of the kinetic experiments varied from 0.1 K/s to 100 K/s; hence the range of viscosity inferred from this method covers 3 orders of magnitude. The data from this method show that viscosity increases with increasing pressure from 1 GPa to 3 GPa for hydrous rhyolitic melts with water content ?0.8 wt% in the high viscosity range. We also measured viscosity of rhyolitic melt with 0.13 wt% water using the parallel-plate viscometer at pressures 0.2 and 0.4 GPa in an internally-heated pressure vessel. The data show that viscosity of rhyolitic melt with 0.13 wt% water decreases with increasing pressure. Combining our new data with literature data, we develop a viscosity model of rhyolitic melts as a function of temperature, pressure and water content.     

Molecular dynamics studies of CaAl2Si2O8 liquid. Part II: Equation of state and a thermodynamic model

A thermodynamic model and equation of state (EOS) is developed from the molecular dynamics simulation experiments of Spera et al. (2009) for CaAl₂Si₂O₈ liquid over the temperature range 3500-6000 K and pressure interval 0-125 GPa. The model is constructed utilizing the isothermal Universal EOS of Vinet et al. (1986) combined with an expression for the temperature-dependence of the internal energy derived from density functional theory (Rosenfeld and Tarazona, 1998). It is demonstrated that this model is more successful at reproducing the data than the temperature-dependent Universal EOS (Vinet et al., 1987) or the volume-explicit EOS of Ghiorso (2004a). Distinct parameterizations are required to model low (<20 GPa) and high (>20 GPa) pressure regimes. This result is ascribed to the affect of liquid structure on macroscopic thermodynamic properties, specifically the interdependence of average cation-oxygen coordination number on the bulk modulus. The thermodynamic transition between the high- and low-pressure parameterizations is modeled as second order, although the nature of the transition is open to question and may well be first order or lambda-like in character.Analysis of the thermodynamic model reveals a predicted region of liquid-liquid un-mixing at low-temperatures (<1624 K) and pressures (<1.257 GPa). These pressure-temperature conditions are above the glass transition temperature but within the metastable liquid region. They represent the highest temperatures yet suggested for liquid-liquid un-mixing in a silicate bulk composition. A shock wave Hugoniot curve is calculated for comparison with the experimental data of Rigden et al. (1989) and of Asimow and Ahrens (2008). The comparison suggests that the model developed in this paper underestimates the density of the liquid by roughly 10% at pressures greater than ∼20 GPa.     

The system CaO-MgO-SiO2-CO2 at 1 GPa, metasomatic wehrlites, and primary carbonatite magmas

New experimental data in CaO-MgO-SiO₂-CO₂ at 1 GPa define the vapor-saturated silicate-carbonate liquidus field boundary involving primary minerals calcite, forsterite and diopside. The eutectic reaction for melting of model calcite (1% MC)-wehrlite at 1 GPa is at 1100 °C, with liquid composition (by weight) 72% CaCO₃ (CC), 9% MgCO₃ (MC), and 18% CaMgSi₂O₆ (Di). These data combined with previous results permit construction of the isotherm-contoured vapor-saturated liquidus surface for the calcite/dolomite field, and part of the adjacent forsterite and diopside fields. Nearly pure calcite crystals in mantle xenoliths cannot represent equilibrium liquids. We recently determined the complete vapor-saturated liquidus surface between carbonates and model peridotites at 2.7 GPa; the peritectic reaction for dolomite (25% MC)-wehrlite at 2.7 GPa occurs at 1300 °C, with liquid composition 60% CC, 29% MC, and 11% Di. The liquidus field boundaries on these two surfaces provide the road-map for interpretation of magmatic processes in various peridotite-CO₂ systems at depths between the Moho and about 100 km. Relationships among kimberlites, melilitites, carbonatites and the liquidus phase boundaries are discussed. Experimental data for carbonatite liquid protected by metasomatic wehrlite have been reported. The liquid trends directly from dolomitic towards CaCO₃ with decreasing pressure. The 1.5 GPa liquid contains 87% CC and 4% Di, much lower in silicate components than our phase boundary. However, the liquids contain approximately the same CaCO₃ (90 ± 1 wt%) in terms of only carbonate components. For CO₂-bearing mantle, all magmas at depth must pass through initial dolomitic compositions. Rising dolomitic carbonatite melt will vesiculate and may erupt as primary magmas through cracks from about ˜70 km. If it percolates through metasomatic wehrlite from 70 km toward the Moho at 35–40 km, primary calcic siliceous carbonatite magma can be generated with silicate content at least 11–18% (70–40 km) on the silicate-carbonate boundary. Received: 22 June 1998 / Accepted: 7 July 1999     

Melts in the mantle modeled in the system CaO-MgO-SiO2-CO2 at 2.7 GPa

The effect of CO₂ on mantle peridotites is modeled by experimental data for the system CaO-MgO-SiO₂-CO₂ at 2.7 GPa. The experiments provide isotherms for the vapor-saturated liquidus surface, bracket piercing points for field boundaries on the surface, and define the positions and compositions of isobaric invariant liquids on the boundaries (eutectics and peritectics). CO₂-saturated carbonatitic liquids (>80% carbonate) exist through approximately 200 °C above the solidus, with a transition to silicate liquids (>80% silicate) within ∼75 °C across a plateau on the liquidus. Carbonate-rich magmas cannot cross the silicate-carbonate liquidus field boundary, so the carbonate liquidus field is therefore a forbidden volume for liquid magmas. This confirms the fact that rounded, pure carbonates in mantle xenoliths cannot represent original liquids. A P-T diagram is constructed for the carbonation and melting reactions for mineral assemblages corresponding to lherzolite, harzburgite, websterite and wehrlite, with carbonate, CO₂ vapor (V), or both. The changing compositions of liquids in solidus reactions on the P-T diagram are illustrated by the changing compositions of eutectic and peritectic liquids on the liquidus surface. At an invariant point Q (∼2.8 GPa/1230 °C), all peridotite assemblages coexist with a calcite-dolomite solid solution (75 ± 5% CaCO₃) and a dolomitic carbonatite melt [57% CaCO₃ (CC), 33% MgCO₃ (MC), 10% CaMgSi₂O₆ (Di)], with 63% CC in the carbonate component. At higher pressures, dolomite-lherzolite, dolomite-harzburgite-V, and dolomite-websterite-V melt to yield similar liquids. Magnesian calcite-wehrlite is the only peridotite melting to carbonatitic liquids (more calcic) at pressures below Q (∼70 km). Dolomitic carbonatite magma rising through mantle to the near-isobaric solidus ledge near Q will begin to crystallize, releasing CO₂ (enhancing crack propagation), and metasomatizing lherzolite toward wehrlite. Received: 20 March 1998 / Accepted: 7 July 1999     

Interaction of model peridotite with H2O-KCl fluid: Experiment at 1.9 GPa and its implications for upper mantle metasomatism

Mineralogical and geochemical data suggest that chloride components play an important role in the transformation and partial melting of upper mantle peridotites. The effect of KCl on the transformation of hydrous peridotite rich in Al₂O₃, CaO, and Na₂O was examined in experiments aimed at studying interaction between model NCMAS peridotite with H₂O-KCl fluid under a pressure of 1.9 GPa, temperatures of 900–1200°C, and various initial H₂O/KCl ratios. The experimental results indicate that KCl depresses the solidus temperature of the hydrous peridotite: this temperature is <900°C at 1.9 GPa, which is more than 100°C lower than the solidus temperature (1000–1025°C) of hydrous peridotite in equilibrium with KCl-free fluid. The reason for the decrease in the melting temperature is that the interaction of KCl with silicates prevails over the effect of chloride on the water activity in the fluid. Experimental data highlight the key role of Al₂O₃ as a component controlling the whole interaction process between peridotite and H₂O-KCl fluid. Garnet, spinel, and pargasite-edenite amphibole in association with aluminous orthopyroxene are unstable in the presence of H₂O-KCl fluid at a chloride concentration in the fluid as low as approximately 2 wt % and are replaced by Cl-bearing phlogopite (0.4–1.1 wt % Cl). Interaction with H₂O-KCl fluid does not, however, affect clinopyroxene and forsterite, which are the Al poorest phases of the system. Chlorine stabilizes phlogopite at relatively high temperatures in equilibrium with melt at temperatures much higher than the solidus (>1200°C). The compositional evolution of melt generated during the melting of model peridotite in the presence of H₂O-KCl fluid is controlled, on the one hand, by the solubility of the H₂O-KCl fluid in the melt and, on the other hand, by phlogopite stability above the solidus. At temperatures below 1050°C, at which phlogopite does not actively participate in melting reactions, fluid dissolution results in SiO₂-undersaturated (35–40 wt %) and MgO-enriched (up to 45 wt %) melts containing up to 4–5 wt % K₂O and 2–3 wt % Cl. At higher temperatures, active phlogopite dissolution and, perhaps, also the separation of immiscible aqueous chloride liquid give rise to melts containing >10 wt % K₂O and 0.3–0.5 wt % Cl. Our experimental results corroborate literature data on the transformation of upper mantle peridotites into phlogopite-bearing associations and the formation of ultrapotassic and highly magnesian melts.     

Bromine cycle in subduction zones through in situ Br monitoring in diamond anvil cells

The geochemical partitioning of bromine between hydrous haplogranitic melts, initially enriched with respect to Br and aqueous fluids, has been continuously monitored in situ during decompression. Experiments were carried out in diamond anvil cells from 890 °C to room temperature and from 1.7 GPa to room pressure, typically from high P, T conditions corresponding to total miscibility (presence of a supercritical fluid). Br contents were measured in aqueous fluids, hydrous melts and supercritical fluids. Partition coefficients of bromine were characterized at pressure and temperature between fluids, hydrous melts and/or glasses, as appropriate: D^Br_fluid/melt = (Br)_fluid/(Br)_melt, ranges from 2.18 to 9.2 ± 0.5 for conditions within the ranges 0.66-1.7 GPa, 590-890 °C; and D^Br_fluid/glass = (Br)_fluid/(Br)_glass ranges from 60 to 375 at room conditions. The results suggest that because high pressure melts and fluids are capable of accepting high concentrations of bromine, this element may be efficiently removed from the slab to the mantle source of arc magmas. We show that Br may be highly concentrated in subduction zone magmas and strongly enriched in subduction-related volcanic gases, because its mobility is strongly correlated with that of water during magma degassing. Furthermore, our experimental results suggest that a non negligible part of Br present in the subducted slab may remain in the down-going slab, being transported toward the transition zone. This indicates that the Br cycle in subduction zones is in fact divided in two related but independent parts: (1) a shallower one where recycled Br may leave the slab with a water and silica-bearing “fluid” leading to enriched arc magmas that return Br to the atmosphere. (2) A deeper cycle where Br may be recycled back to the mantle maybe to the transition zone, where it may be present in high pressure water-rich metasomatic fluids.     

An experimental study of the solubility of molybdenum in H2O and KCl-H2O solutions from 500 °C to 800 °C, and 150 to 300 MPa

The solubility of molybdenum (Mo) was determined at temperatures from 500 °C to 800 °C and 150 to 300 MPa in KCl-H₂O and pure H₂O solutions in cold-seal experiments. The solutions were trapped as synthetic fluid inclusions in quartz at experimental conditions, and analyzed by laser ablation inductively coupled plasma mass spectrometry (LA ICPMS).Mo solubilities of 1.6 wt% in the case of KCl-bearing aqueous solutions and up to 0.8 wt% in pure H₂O were found. Mo solubility is temperature dependent, but not pressure dependent over the investigated range, and correlates positively with salinity (KCl concentration). Molar ratios of ∼1 for Mo/Cl and Mo/K are derived based on our data. In combination with results of synchrotron X-ray absorption spectroscopy of individual fluid inclusions, it is suggested that Mo-oxo-chloride complexes are present at high salinity (>20 wt% KCl) and ion pairs at moderate to low salinity (<11 wt% KCl) in KCl-H₂O aqueous solutions. Similarly, in the pure H₂O experiments molybdic acid is the dominant species in aqueous solution. The results of these hydrothermal Mo experiments fit with earlier studies conducted at lower temperatures and indicate that high Mo concentrations can be transported in aqueous solutions. Therefore, the Mo concentration in aqueous fluids seems not to be the limiting factor for ore formation, whereas precipitation processes and the availability of sulfur appear to be the main controlling factors in the formation of molybdenite (MoS₂).     

Molecular dynamics simulation of the CH4 and CH4-H2O systems up to 10 GPa and 2573 K

Based on our previous study of the intermolecular potential for pure H₂O and the strict evaluation of the competitive potential models for pure CH₄ and the ab initio fitting potential surface across CH₄-H₂O molecules in this study, we carried out more than two thousand molecular dynamics simulations for the PVTx properties of pure CH₄ and the CH₄-H₂O mixtures up to 2573 K and 10 GPa. Comparison of 1941 simulations with experimental PVT data for pure CH₄ shows an average deviation of 0.96% and a maximum deviation of 2.82%. The comparison of the results of 519 simulations of the mixtures with the experimental measurements reveals that the PVTx properties of the CH₄-H₂O mixtures generally agree with the extensive experimental data with an average deviation of 0.83% and 4% in maximum, which is equivalent to the experimental uncertainty. Moreover, the maximum deviation between the experimental data and the simulation results decreases to about 2% as temperature and pressure increase, indicating that the high accuracy of the simulation is well retained in the high temperature and pressure region.After the validation of the simulation method and the intermolecular potential models, we systematically simulated the PVTx properties of this binary system from 673 K and 0.05 GPa to 2573 K and 10 GPa. In order to integrate all the simulation results and the experimental data for the calculation of thermodynamic properties, an equation of state (EOS) is developed for the CH₄-H₂O system covering 673-2573 K and 0.01-10 GPa. Isochores for compositions <4 mol% CH₄ up to 773 K and 600 MPa are also determined in this paper. The program for the EOS can be downloaded from www.geochem-model.org/programs.htm.     

Structure, thermodynamic and transport properties of liquid MgSiO3: Comparison of molecular models and laboratory results

Liquid MgSiO₃ is a model for the Earth’s magma ocean and of remnant melt present near the core-mantle boundary. Here, models for molten MgSiO₃ are computed employing empirical potential molecular dynamics (EPMD) and results are compared to published results including two EPMD studies and three first-principles molecular dynamics (FPMD) models and to laboratory data. The EPMD results derived from the Oganov (OG) potential come closest to the density of MgSiO₃ liquid at the 1-bar melting point inferred from the melting curve. At higher P, EPMD densities calculated from the OG potential and FPMD broadly match shock wave studies, with the OG potential yielding the better comparison. Matsui (M) potential results deviate from other studies above ∼50 GPa. Overall, results based on the OG potential compare best to experimental densities over the P-T range of the mantle. Isothermally, upon increasing P the mean coordination numbers () of oxygen around Si and Mg monotonically increase with pressure. Tetrahedral Si and octahedral Si monotonically increase and decrease, respectively, whereas pentahedral Si maximizes at 10-20 GPa. Tetrahedral Mg decreases monotonically as P increases whereas pentahedral, octahedral and higher coordination polyhedra each show similar behavior first increasing and then decreasing after attaining a maximum; the P of the maximum for each polyhedra type migrates to higher P as the CN increases. Free oxygen and oxygen with one nearest neighbor of either Si or Mg decreases whereas Si or Mg with two or three nearest oxygens (i.e., tricluster oxygen) increases with increasing P isothermally. The increase of tricluster oxygen is consistent with spectroscopy on MgSiO₃ glass quenched from 2000 K and 0-40 GPa and high-energy X-ray studies constraining the coordination of O around Mg and around Si at 2300 K and 1 bar. Coordination statistics from FPMD studies for O around Si and Si around O are in agreement with the EPMD results based on the M and OG potentials. Mg self-diffusivity is greater than O and Si self-diffusivities for both the M and OG potentials. All D values monotonically decrease with increasing pressure isothermally and all atoms are more diffusive in the M liquid compared to the OG liquid except at T > ∼5000 K and P > 100 GPa. Previously published EPMD diffusivities fall between values given by the M and OG potentials, at least up to 45 GPa. The M liquid is generally less viscous than the OG liquid except at P > ∼80 GPa. Activation energy and volume are around 96 kJ/mol and 1.5 cm³/mol, respectively. The FPMD viscosity results at 120 GPa and 4000 and 4500 K are essentially identical to the values from the M and OG potentials. FPMD viscosity results are similar to the OG results for P < 60 GPa; at higher P, the FPMD viscosities are higher. At 4000 K and 100 GPa the shear viscosity of liquid MgSiO₃ is ∼0.1 Pa s. More extensive laboratory results are required to better define the thermodynamic, transport and structural properties of MgSiO₃ liquids and for comparison with computational studies.     

The stability of hydrous phases beyond antigorite breakdown for a magnetite-bearing natural serpentinite between 6.5 and 11 GPa

Phase relations for a natural serpentinite containing 5 wt% of magnetite have been investigated using a multi-anvil apparatus between 6.5 and 11 GPa and 400–850 °C. Post-antigorite hydrous phase assemblages comprise the dense hydrous magnesium silicates (DHMSs) phase A (11.3 wt% H₂O) and the aluminous phase E (Al-PhE, 11.9 wt% H₂O). In addition, a ferromagnesian hydrous silicate (11.1 wt% H₂O) identified as balangeroite (Mg,Fe)₄₂Si₁₆O₅₄(OH)₄₀, typically described in low pressure natural serpentinite, was found coexisting with Al-PhE between 650 and 700 °C at 8 GPa. In the natural antigorite system, phase E stability is extended to lower pressures (8 GPa) than previously reported in simple chemical systems. The reaction Al-phase E?=?garnet?+?olivine?+?H₂O is constrained between 750 and 800 °C between 8 and 11 GPa as the terminal boundary between hydrous mineral assemblages and nominally anhydrous assemblages, hence restricting water transfer into the deep mantle to the coldest slabs. The water storage capacity of the assemblage Al-PhE?+?enstatite (high-clinoenstatite)?+?olivine, relevant for realistic hydrated slab composition along a relatively cold temperature path is estimated to be ca. 2 wt% H₂O. Attempts to mass balance run products emphasizes the role of magnetite in phase equilibria, and suggests the importance of ferric iron in the stabilization of hydrous phases such as balangeroite and aluminous phase E.