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.     

Sulfur concentration at sulfide saturation (SCSS) in magmatic silicate melts

The sulfur concentration in silicate melts at sulfide saturation (SCSS) was experimentally investigated in a temperature range from 1150 to 1450 °C and a pressure range from 500 MPa to 1 GPa in a piston-cylinder apparatus. The investigated melt compositions varied from rhyolitic to basaltic and water concentrations varied from 0 to ∼9 wt%. All experiments were saturated with FeS melt or pyrrhotite crystals. Temperature was confirmed to have a positive effect on the SCSS. Experimental oxygen fugacities were either near the carbon-carbon monoxide buffer or one log unit above the nickel-nickel oxide buffer, and found to positively affect the SCSS. Combining our results with data from the literature we constructed a model to predict the SCSS in melts ranging in composition from komatiitic to rhyolitic, with water concentrations from 0 to 9 wt%, at pressures from 1 bar to 9 GPa and oxygen fugacities between ∼2 log units below the fayalite-magnetite-quartz buffer to ∼2 log units above it. The coefficients were obtained by multiple linear regression of experimental data and the best model found for the prediction of the SCSS is:

     

Viscosity of flux-rich pegmatitic melts

Viscosity experiments were conducted with two flux-rich pegmatitic melts PEG0 and PEG2. The Li₂O, F, B₂O₃ and P₂O₅ contents of these melts were 1.04, 4.06, 2.30 and 1.68 and 1.68, 5.46, 2.75 and 2.46 wt%, respectively. The water contents varied from dry to 9.04 wt% H₂O. The viscosity was determined in internally heated gas pressure vessels using the falling sphere method in the temperature range 873–1,373 K at 200 and 320 MPa pressure. At 1,073 K, the viscosity of water-rich (~9 wt% H₂O) melts is in the range of 3–60 Pa s, depending on the melt composition. Extrapolations to lower temperature assuming an Arrhenian behavior indicate that highly fluxed pegmatite melts may reach viscosities of ~30 Pa s at 773 K. However, this value is a minimum estimation considering the strongly non-Arrhenian behavior of hydrous silicate melts. The experimentally determined melt viscosities are lower than the prediction of current models taking compositional parameters into account. Thus, these models need to be improved to predict accurately the viscosity of flux-rich water bearing melts. The data also indicate that Li influences significantly the melt viscosity. Decreasing the molar Al/(Na + K + Li) ratio results in a strong viscosity decrease, and highly fluxed melts with low Al/(Na + K + Li) ratios (~0.8) have a rheological behavior which is very close to that of supercritical fluids.     

Molecular dynamics investigation of viscosity, chemical diffusivities and partial molar volumes of liquids along the MgO-SiO2 join as functions of pressure

Molecular dynamics simulations are carried out to systematically address the effects of composition and pressure on melts along the MgO-SiO₂ join and elucidate the role of structural modifier content on silicate melt properties. The MgO-SiO₂ system shows non-ideal mixing with a negative excess volume of mixing at low pressures, but the mixing becomes closer to ideal at higher pressures. At atmospheric pressure, the viscosities and diffusivities vary by more than 3 orders of magnitude as the composition is varied along this join, with the low SiO₂ melts characterized by lower viscosities and higher diffusivities; these results are in quantitative agreement with experimental results for the dependence of viscosity and diffusivity on structural modifier content in a wide range of silicate systems. The transport properties of melts in this system converge at higher pressures; at pressures greater than ∼15 GPa the viscosity and diffusivities vary by less than an order of magnitude across the entire range of compositions. The relevance of equations that relate the viscosity and diffusivity is also addressed.     

A long-duration experiment on hydrous species geospeedometer and hydrous melt viscosity

A 2.4-year controlled-cooling-rate experiment was carried out to investigate the dependence of hydrous species concentrations in rhyolitic melt on cooling rate. The experiment allows us to obtain speciation for a cooling rate of 1.68 × 10⁻⁶ K/s, extending previous experimental data by two orders of magnitude. Furthermore, a viscosity as high as 10^17.2 Pa s is inferred for this hydrous rhyolitic melt with 0.85 wt% total H₂O at 671 K. The results are applied to examine whether a geospeedometry model and four viscosity models may be extrapolated to slower cooling rates or lower temperatures. Two of the viscosity models and the geospeedometry model can be extrapolated by two orders of magnitude upwards in terms of viscosity or downwards in terms of cooling rate.     

Viscosity of silicate melts in CaMgSi2O6–NaAlSi2O6 system at high pressure

In situ X-ray viscometry of the silicate melts was carried out at high pressure and at high temperature. The viscosity of the silicate melts in the diopside(Di)–jadeite(Jd) system was determined in the pressure range from 1.88 GPa to 7.9 GPa and in the temperature range from 2,003 K to 2,173 K. The viscosity of the Di 25%–Jd 75% melt decreases continuously to 5.0 GPa, whereas the viscosity of the Di 50%–Jd 50% melt increases over 3.5 GPa. The viscosity of the Di50%–Jd 50% melt reaches a minimum around 3.5 GPa. Since the amounts of silicon in the two melts are the same, the difference in the pressure dependence of the viscosity may be controlled by another network-forming element, i.e., aluminum. The difference in the pressure dependence of the viscosities in the melts with two intermediate compositions in the Di–Jd system is estimated to be due to the difference in the melt structures at high pressures and high temperatures.     

A model for the viscosity of rhyolite as a function of H2O-content and pressure: A calibration based on centrifuge piston cylinder experiments

The Newtonian viscosity of synthetic rhyolitic liquids with 0.15-5.24 wt% dissolved water was determined in the interval between 580 and 1640 °C and pressures of 1 atm and 5-25 kbar. Measurements were performed by combining static and accelerated (up to 1000g) falling sphere experiments on water-bearing samples, with high temperature concentric cylinder experiments on 0.15 wt% H₂O melts. These methods allowed viscosity determinations between 10² and 10⁷ Pa s, and cover the complete range of naturally occurring magmatic temperatures, pressures, and H₂O-contents for rhyolites.Our viscosity data, combined with those from previous studies, were modeled by an expression based on the empirical Vogel-Fulcher-Tammann equation, which describes viscosities and derivative properties (glass transition temperature T_g, fragility m, and activation volume of viscous flow Va) of silicic liquids as a function of P-T-X_(H2O). The fitted expressions do not account for composition-dependent parameters other than X_(H2O) and reproduce the entire viscosity database for silicic liquids to within 3.0% average relative error on log η (i.e. std. error of estimate of 0.26 log units).The results yield the expected strong decrease of viscosity with temperature and water content, but show variable pressure dependencies. Viscosity results to be strongly affected by pressure at low pressures; an effect amplified at low temperatures and water contents. Fragility, as a measure for the deviation from Arrhenian behavior, decreases with H₂O-content but is insensitive to pressure. Activation volumes are always largely negative (e.g., less than −10 cm³/mol) and increase strongly with H₂O-content. Variations in melt structure that may account for the observed property variations are discussed.     

High-pressure melting relations in Fe-C-S systems: Implications for formation, evolution, and structure of metallic cores in planetary bodies

We present new high-pressure temperature experiments on melting phase relations of Fe-C-S systems with applications to metallic core formation in planetary interiors. Experiments were performed on Fe-5 wt% C-5 wt% S and Fe-5 wt% C-15 wt% S at 2-6 GPa and 1050-2000 °C in MgO capsules and on Fe-13 wt% S, Fe-5 wt% S, and Fe-1.4 wt% S at 2 GPa and 1600 °C in graphite capsules. Our experiments show that: (a) At a given P-T, the solubility of carbon in iron-rich metallic melt decreases modestly with increasing sulfur content and at sufficiently high concentration, the interaction between carbon and sulfur can cause formation of two immiscible melts, one rich in Fe-carbide and the other rich in Fe-sulfide. (b) The mutual solubility of carbon and sulfur increases with increasing pressure and no super-liquidus immiscibility in Fe-rich compositions is likely expected at pressures greater than 5-6 GPa even for bulk compositions that are volatile-rich. (c) The liquidus temperature in the Fe-C-S ternary is significantly different compared to the binary liquidus in the Fe-C and Fe-S systems. At 6 GPa, the liquidus of Fe-5 wt% C-5 wt% S is 150-200 °C lower than the Fe-5 wt% S. (d) For Fe-C-S bulk compositions with modest concentration of carbon, the sole liquidus phase is iron carbide, Fe₃C at 2 GPa and Fe₇C₃ at 6 GPa and metallic iron crystallizes only with further cooling as sulfur is concentrated in the late crystallizing liquid. Our results suggest that for carbon and sulfur-rich core compositions, immiscibility induced core stratification can be expected for planets with core pressure less than ∼6 GPa. Thus planetary bodies in the outer solar system such as Ganymede, Europa, and Io with present day core-mantle boundary (CMB) pressures of ∼8, ∼5, and 7 GPa, respectively, if sufficiently volatile-rich, may either have a stratified core or may have experienced core stratification owing to liquid immiscibility at some stage of their accretion. A similar argument can be made for terrestrial planetary bodies such as Mercury and Earth’s Moon, but no such stratification is predicted for cores of terrestrial planets such as Earth, Venus, and Mars with the present day core pressure in the order ?136 GPa, ?100 GPa, and ?23 GPa. (e) Owing to different expected densities of Fe-rich (and carbon-bearing) and sulfur-rich metallic melts, their settling velocities are likely different; thus core formation in terrestrial planets may involve rain of more than one metallic melt through silicate magma ocean. (f) For small planetary bodies that have core pressures <6 GPa and have a molten core or outer core, settling of denser carbide-rich liquid or flotation of lighter, sulfide-rich melt may contribute to an early, short-lived geodynamo.     

The composition of liquids coexisting with dense hydrous magnesium silicates at 11-13.5 GPa and the endpoints of the solidi in the MgO-SiO2-H2O system

High-pressure liquids in the MgO-SiO₂-H₂O (MSH) system have been investigated at 11 and 13.5 GPa and between 1000 and 1350 °C. A bulk composition more magnesian than the tie-line forsterite-H₂O was employed for the study. Rocking multi-anvil experiments were combined with a diamond trap set-up. After termination of the experiments, the liquid trapped in the diamond layer was analysed by laser ablation ICP-MS using the ‘freezing’ technique. At 11 GPa, liquids coexist with one or two of phase A, clinohumite, chondrodite, and forsterite. A marked discontinuity in the evolution of liquid compositions near 1100 °C is observed at 11 GPa. A step of ∼13 wt% H₂O and 13 wt% MgO is interpreted to result from overstepping the fluid-saturated solidus reaction mass balanced to 1.00(18) phase A + 1.07(4) fluid = 0.63(15) chondrodite + 1.44(2) melt. At 13.5 GPa liquids coexist with one or two of hydrous wadsleyite, clinohumite, superhydrous B, phase B, and forsterite. The discontinuity in liquid composition is no longer present, indicating that the second critical endpoint of the solidus has been overstepped. Thus, hydrous melts in the Mg-rich part of the MSH system (molar bulk Mg/Si > 2) are chemically distinct from aqueous fluids at pressure up to 11 GPa. Convergence of fluid and melt compositions along the solidus resulting in a supercritical liquid occurs between 11 and 13.5 GPa, at which pressure the entire MSH system becomes supercritical.     

The effect of water on the viscosity of a haplogranitic melt under P-T-X conditions relevant to silicic volcanism

 The viscosities of hydrous haplogranitic melts synthesized by hydrothermal fusion at 2 kbar pressure and 800 to 1040° C have been measured at temperatures just above the glass transition and at a pressure of 1 bar using micropenetration techniques. The micropenetration viscometry has been performed in the viscosity range of 10⁹ Pa s to 10¹² Pa s. The samples ranged in water content from 0.4 to 3.5 wt%. For samples with up to 2.5 wt% H₂O, the water contents have been determined using infrared spectroscopy obtained before and after each viscometry experiment to be constant over the duration of the measurements. Above this water content a measurable loss of water occurs during the viscometry. The viscosity data illustrate an extremely nonlinear decrease in viscosity with added water. The viscosity drops drastically with the addition of 0.5 wt% of water and then shallows out at water contents of 2 wt%. An additional viscosity datum point obtained from the analysis of fluid inclusions in a water-saturated HPG8 confirms a near invariance of the viscosity with the addition of water between 2 and 6 wt%. These measurements may be compared directly with the data of Hess et al. (1995, in press) for the effects of excess alkali and alkaline earth oxides on the viscosity of HPG8 (also obtained at 1 bar). The viscosity of the melts, compared on an equivalent molar basis, increases in the order H₂O<(Li₂O<Na₂O< K₂O<Rb₂O,Cs₂O<BaO<SrO<CaO<MgO< BeO). The extraordinary decrease in melt viscosity with added water is poorly reproduced by the calculation scheme of Shaw (1972) for the range of water contents investigated here. The speciation of water in the quenched glasses can be used to quantify the dependence of the viscosity on hydroxyl content. Considering only the hydroxyl groups as active fluidizers in the hydrous melts the nonlinearity of the viscosity decrease and the difference with the effects of the alkali oxides becomes larger. Consequences for degassing calcalkaline rhyolite are discussed. Received: 17 August 1995/Accepted: 8 January 1996     

The rheological evolution of alkaline Vesuvius magmas and comparison with alkaline series from the Phlegrean Fields, Etna, Stromboli and Teide

Somma-Vesuvius is considered one of the highest-risk volcanic systems in the world due to its high population density and record of highly destructive explosive activity. Eruptive style at Vesuvius varies greatly, alternating between effusive and explosive activities, and is likely strongly controlled by the evolution of the physical and chemical properties of the magma. Nevertheless, with the exception of the 1631 eruption, the rheological properties of Vesuvius magmas remain largely unconstrained. Here, we investigate the Newtonian shear viscosity (η) of dry and hydrous melts from the Mercato (plinian) and 1906 (violent strombolian) eruptions. These eruptions differ in size, eruptive style and magma chemistry (from phonolite to phono-tephrite). To evaluate the dry liquid viscosity variation covered by the eruptive products of the recent activity at Vesuvius, we measured the melt viscosities of bulk rock compositions and, for highly crystalline samples, of the separated groundmasses of tephras from the Pollena and 1906 eruptions. Hydrated samples with up to 4.24 wt% dissolved water were synthesised in a piston cylinder apparatus at confining pressure up to 10 kbar. The dry high temperature and the dry and hydrous low-temperature viscosities were obtained by combining the concentric cylinder and micropenetration techniques. The measured viscosities were parameterized by a modified Vogel-Fulcher-Tammann equation, accounting for the effect of water content, and were compared with previous measurements and models. At magmatic temperatures, the viscosities of Mercato samples are about four orders of magnitude higher than that of the least viscous investigated products from the 1906 eruption. Complex numerical models to forecast eruptive scenarios and their environmental impact are extremely sensitive to the accuracy of the input parameters and constitutive equations of magma properties. As a consequence, the numerical expressions obtained here are of particular relevance in the context of hazard assessment related to the different possible eruptive scenarios at Vesuvius through numerical simulation tools. The effect of composition on the liquid viscosities is compared to other high-Na (e.g., samples from Teide and Etna) and high-K (e.g., samples from Stromboli and Phlegrean Fields) alkaline magmas.     

The influence of pressure and composition on the viscosity of andesitic melts

The effect of pressure and composition on the viscosity of both anhydrous and hydrous andesitic melts was studied in the viscosity range of 10⁸ to 10^11.5 Pa · s using parallel plate viscometry. The pressure dependence of the viscosity of three synthetic, iron-free liquids (andesite analogs) containing 0.0, 1.06, and 1.96 wt.% H₂O, respectively, was measured from 100 to 300 MPa using a high-P-T viscometer. These results, combined with those from Richet et al. (1996), indicate that viscosities of anhydrous andesitic melts are independent of pressure, whereas viscosities of hydrous melts slightly increase with increasing pressure. This trend is consistent with an increased degree of depolymerization in the hydrous melts. Compositional effects on the viscosity were studied by comparing iron-free and iron-bearing compositions with similar degrees of depolymerization. During experiments at atmospheric and at elevated pressures (100 to 300 MPa), the viscosity of iron-bearing anhydrous melts preequilibrated in air continuously increased, and the samples became paramagnetic. Analysis of these samples by transmission electron microscopy showed a homogeneous distribution of crystals (probably magnetite) with sizes in the range of 10 to 50 nm. No significant difference in the volume fractions of crystals was found in samples after annealing for 170 to 830 min at temperatures ranging from 970 to 1122 K. An iron-bearing andesite containing 1.88 wt.% H₂O, which was synthesized at intrinsic fO₂ conditions in an internally heated pressure vessel, showed a similar viscosity behavior as the anhydrous melts. The continuous increase in viscosity at a constant temperature is attributed to changes of the melt structure due to exsolution of iron-rich phases. By extrapolating the time evolution of viscosity down to the time at which the run temperature was reached, for both the anhydrous (at 1055 K) and the hydrous (at 860 K) iron-bearing andesite, the viscosity is 0.7 log units lower than predicted by the model of Richet et al. (1996). This may be explained by differences in structural properties of Fe²⁺ and Fe³⁺ and their substitutes Mg²⁺, Ca²⁺, and Al³⁺, which were used in the analogue composition.The effect of iron redox state on the viscosity of anhydrous, synthetic andesite melts was studied at ambient pressure using a dilatometer. Reduced iron-bearing samples were produced by annealing melts in graphite crucibles in an Ar/CO atmosphere for different run times. In contrast to the oxidized sample, no variation of viscosity with time and no exsolution of iron oxide phases was observed for the most reduced glasses. This indicates that trivalent iron promotes the exsolution of iron oxide in supercooled melts. With decreasing Fe³⁺/ΣFe ratio from 0.58 to 0.34, the viscosity decreases by ∼1.6 log units in the investigated temperature range between 964 and 1098 K. A more reduced glass with Fe³⁺/ΣFe = 0.21 showed no additional decrease in viscosity. Our conclusion from these results is that the viscosity of natural melts may be largely overestimated when using data obtained from samples synthesized in air.     

Mineral/melt partitioning of trace elements during hydrous peridotite partial melting

This experimental study examines the mineral/melt partitioning of incompatible trace elements among high-Ca clinopyroxene, garnet, and hydrous silicate melt at upper mantle pressure and temperature conditions. Experiments were performed at pressures of 1.2 and 1.6 GPa and temperatures of 1,185 to 1,370 °C. Experimentally produced silicate melts contain up to 6.3 wt% dissolved H ₂O, and are saturated with an upper mantle peridotite mineral assemblage of olivine+orthopyroxene+clinopyroxene+spinel or garnet. Clinopyroxene/melt and garnet/melt partition coefficients were measured for Li, B, K, Sr, Y, Zr, Nb, and select rare earth elements by secondary ion mass spectrometry. A comparison of our experimental results for trivalent cations (REEs and Y) with the results from calculations carried out using the Wood-Blundy partitioning model indicates that H ₂O dissolved in the silicate melt has a discernible effect on trace element partitioning. Experiments carried out at 1.2 GPa, 1,315 °C and 1.6 GPa, 1,370 °C produced clinopyroxene containing 15.0 and 13.9 wt% CaO, respectively, coexisting with silicate melts containing ~1–2 wt% H ₂O. Partition coefficients measured in these experiments are consistent with the Wood-Blundy model. However, partition coefficients determined in an experiment carried out at 1.2 GPa and 1,185 °C, which produced clinopyroxene containing 19.3 wt% CaO coexisting with a high-H ₂O (6.26±0.10 wt%) silicate melt, are significantly smaller than predicted by the Wood-Blundy model. Accounting for the depolymerized structure of the H ₂O-rich melt eliminates the mismatch between experimental result and model prediction. Therefore, the increased Ca ²⁺ content of clinopyroxene at low-temperature, hydrous conditions does not enhance compatibility to the extent indicated by results from anhydrous experiments, and models used to predict mineral/melt partition coefficients during hydrous peridotite partial melting in the sub-arc mantle must take into account the effects of H ₂O on the structure of silicate melts.     

Electrical conductivity of dry and hydrous NaAlSi<Subscript>3</Subscript>O<Subscript>8</Subscript> glasses and liquids at high pressures

We report the first study of electrical conductivities of silicate melts at very high pressures (up to 10 GPa) and temperatures (up to 2,173 K). Impedance spectroscopy was applied to dry and hydrous albite (NaAlSi₃O₈) glasses and liquids (with 0.02–5.7 wt% H₂O) at 473–1,773 K and 0.9–1.8 GPa in a piston-cylinder apparatus, using a coaxial cylindrical setup. Measurements were also taken at 473–2,173 K and 6–10 GPa in two multianvil presses, using simple plate geometry. The electrical conductivity of albite melts is found to increase with temperature and water content but to decrease with pressure. However, at 6 GPa, conductivity increases rapidly with temperature above 1,773 K, so that at temperatures beyond 2,200 K, conductivity may actually increase with pressure. Moreover, the effect of water in enhancing conductivity appears to be more pronounced at 6 GPa than at 1.8 GPa. These observations suggest that smaller fractions of partial melt than previously assumed may be sufficient to explain anomalously high conductivities, such as in the asthenosphere. For dry melt at 1.8 GPa, the activation energy at T > 1,073 K is higher than that at T < 1,073 K, and the inflection point coincides with the rheological glass transition. Upon heating at 6–10 GPa, dry albite glass often shows a conductivity depression starting from ~1,173 K (due to crystallization), followed by rapid conductivity enhancement when temperature approaches the albite liquidus. For hydrous melts at 0.9–1.8 GPa, the activation energies for conductivity at ≥1,373 K are lower than those at <973 K, with a complex transition pattern in between. Electrical conductivity and previously reported Na diffusivity in albite melt are consistent with the Nernst–Einstein relation, suggesting the dominance of Na transport for electrical conduction in albite melts.     

Water diffusion in Mount Changbai peralkaline rhyolitic melt

Diffusion couple experiments with wet half (up to 4.6 wt%) and dry half were carried out at 789–1,516 K and 0.47–1.42 GPa to investigate water diffusion in a peralkaline rhyolitic melt with major oxide concentrations matching Mount Changbai rhyolite. Combining data from this work and a related study, total water diffusivity in peralkaline rhyolitic melt can be expressed as:

$ D_{{{\text{H}}_{ 2} {\text{O}}_{\text{t}} }} = D_{{{\text{H}}_{ 2} {\text{O}}_{\text{m}} }} \left( {1 - \frac{0.5 - X}{{\sqrt {[4\exp (3110/T - 1.876) - 1](X - X^{2} ) + 0.25} }}} \right), $

$ {\text{with}}\;D_{{{\text{H}}_{ 2} {\text{O}}_{\text{m}} }} = \exp \left[ { - 1 2. 7 8 9- \frac{13939}{T} - 1229.6\frac{P}{T} + ( - 27.867 + \frac{60559}{T})X} \right], $

where D is in m² s^?1, T is the temperature in K, P is the pressure in GPa, and X is the mole fraction of water and calculated as X = (C/18.015)/(C/18.015 + (100 ? C)/33.14), where C is water content in wt%. We recommend this equation in modeling bubble growth and volcanic eruption dynamics in peralkaline rhyolitic eruptions, such as the ~1,000-ad eruption of Mount Changbai in North East China. Water diffusivities in peralkaline and metaluminous rhyolitic melts are comparable within a factor of 2, in contrast with the 1.0–2.6 orders of magnitude difference in viscosities. The decoupling of diffusivity of neutral molecular species from melt viscosity, i.e., the deviation from the inversely proportional relationship predicted by the Stokes–Einstein equation, might be attributed to the small size of H₂O molecules. With distinct viscosities but similar diffusivity, bubble growth controlled by diffusion in peralkaline and metaluminous rhyolitic melts follows similar parabolic curves. However, at low confining pressure or low water content, viscosity plays a larger role and bubble growth rate in peralkaline rhyolitic melt is much faster than that in metaluminous rhyolite.

     

Rb and Sr partitioning between haplogranitic melts and aqueous solutions

Rubidium and strontium partitioning experiments between haplogranitic melts and aqueous fluids (water or 1.16-3.56 m (NaCl + KCl) ± HCl) were conducted at 750-950 °C and 0.2-1.4 GPa to investigate the effects of melt and fluid composition, pressure, and temperature. In addition, we studied if the applied technique (rapid and slow quench, and in-situ determination of trace element concentration in the fluid) has a bearing on the obtained data. There is good agreement of the data from different techniques for chloridic solutions, whereas back reactions between fluid and melt upon cooling have a significant effect on results from the experiments with water.The Rb fluid-melt partition coefficient shows no recognizable dependence on melt composition and temperature.For chloridic solutions, it is ∼0.4, independent of pressure. In experiments with water, it is one to two orders of magnitude lower and increases with pressure. The strontium fluid-melt partition coefficient does not depend on temperature. It increases slightly with pressure in Cl free experiments. In chloridic fluids, there is a sharp increase in the Sr partition coefficient with the alumina saturation index (ASI) from 0.003 at an ASI of 0.8 to a maximum of 0.3 at an ASI of 1.05. At higher ASI, it decreases slightly to 0.2 at an ASI of 1.6. It is one to two orders of magnitude higher in chloridic fluids compared to those found in H₂O experiments. The Rb/Sr ratio in non-chloridic solutions in equilibrium with metaluminous melts increases with pressure, whereas the Rb/Sr ratio in chloridic fluids is independent of pressure and decreases with fluid salinity.The obtained fluid-melt partition coefficients are in good agreement with data from natural cogenetic fluid and melt inclusions. Numerical modeling shows that although the Rb/Sr ratio in the residual melt is particularly sensitive to the degree of fractional crystallization, exsolution of a fluid phase, and associated fluid-melt partitioning is not a significant factor controlling Rb and Sr concentrations in the residual melt during crystallization of most granitoids.     

The experimental calibration of the iron isotope fractionation factor between pyrrhotite and peralkaline rhyolitic melt

A first experimental study was conducted to determine the equilibrium iron isotope fractionation between pyrrhotite and silicate melt at magmatic conditions. Experiments were performed in an internally heated gas pressure vessel at 500 MPa and temperatures between 840 and 1000 °C for 120-168 h. Three different types of experiments were conducted and after phase separation the iron isotope composition of the run products was measured by MC-ICP-MS. (i) Kinetic experiments using ⁵⁷Fe-enriched glass and natural pyrrhotite revealed that a close approach to equilibrium is attained already after 48 h. (ii) Isotope exchange experiments—using mixtures of hydrous peralkaline rhyolitic glass powder (∼4 wt% H₂O) and natural pyrrhotites (Fe_1 − xS) as starting materials— and (iii) crystallisation experiments, in which pyrrhotite was formed by reaction between elemental sulphur and rhyolitic melt, consistently showed that pyrrhotite preferentially incorporates light iron. No temperature dependence of the fractionation factor was found between 840 and 1000 °C, within experimental and analytical precision. An average fractionation factor of Δ ⁵⁶Fe/⁵⁴Fe_{pyrrhotite-melt} = −0. 35 ± 0.04‰ (2SE, n = 13) was determined for this temperature range. Predictions of Fe isotope fractionation between FeS and ferric iron-dominated silicate minerals are consistent with our experimental results, indicating that the marked contrast in both ligand and redox state of iron control the isotope fractionation between pyrrhotite and silicate melt. Consequently, the fractionation factor determined in this study is representative for the specific Fe²⁺/ΣFe ratio of our peralkaline rhyolitic melt of 0.38 ± 0.02. At higher Fe²⁺/ΣFe ratios a smaller fractionation factor is expected. Further investigation on Fe isotope fractionation between other mineral phases and silicate melts is needed, but the presented experimental results already suggest that even at high temperatures resolvable variations in the Fe isotope composition can be generated by equilibrium isotope fractionation in natural magmatic systems.     

Electrical conductivity of hydrous basaltic melts: implications for partial melting in the upper mantle

The Earth’s uppermost asthenosphere is generally associated with low seismic wave velocity and high electrical conductivity. The electrical conductivity anomalies observed from magnetotelluric studies have been attributed to the hydration of mantle minerals, traces of carbonatite melt, or silicate melts. We report the electrical conductivity of both H₂O-bearing (0–6 wt% H₂O) and CO₂-bearing (0.5 wt% CO₂) basaltic melts at 2 GPa and 1,473–1,923 K measured using impedance spectroscopy in a piston-cylinder apparatus. CO₂ hardly affects conductivity at such a concentration level. The effect of water on the conductivity of basaltic melt is markedly larger than inferred from previous measurements on silicate melts of different composition. The conductivity of basaltic melts with more than 6 wt% of water approaches the values for carbonatites. Our data are reproduced within a factor of 1.1 by the equation log σ = 2.172 − (860.82 − 204.46 w ^0.5)/(T − 1146.8), where σ is the electrical conductivity in S/m, T is the temperature in K, and w is the H₂O content in wt%. We show that in a mantle with 125 ppm water and for a bulk water partition coefficient of 0.006 between minerals and melt, 2 vol% of melt will account for the observed electrical conductivity in the seismic low-velocity zone. However, for plausible higher water contents, stronger water partitioning into the melt or melt segregation in tube-like structures, even less than 1 vol% of hydrous melt, may be sufficient to produce the observed conductivity. We also show that ~1 vol% of hydrous melts are likely to be stable in the low-velocity zone, if the uncertainties in mantle water contents, in water partition coefficients, and in the effect of water on the melting point of peridotite are properly considered.     

Trace element geochemistry of the 1991 Mt. Pinatubo silicic melts, Philippines: Implications for ore-forming potential of adakitic magmatism

The dacite pumice erupted from Mt. Pinatubo on June 15, 1991 (whole-rock, rhyolitic groundmass glasses and homogenized melt inclusions) has been analyzed using inductively coupled plasma-mass spectrometry (ICP-MS), nanosecond and femtosecond laser ablation ICP-MS and secondary ion mass spectrometry (SIMS) to evaluate its ore-forming potential. Data suggest that adakite magmas are metal-rich and concentrate ore metals during magmatic differentiation. Sulfides segregate in limited amounts under the hydrous, oxidizing conditions typical of adakitic magmas resulting in incompatible behavior for Au (6-22 ppb), Cu (26-77 ppm), and Pb, Mo, As, and Sb in melts of dacitic to rhyolitic compositions. Metal transfer from this adakite magma to the coexisting aqueous phase was favored by the peraluminous composition of the rhyolitic melt and high aqueous chloride concentrations. Mass balance calculations suggest that the pre-eruptive aqueous phase could have extracted a minimum of 100 t Au and 5 × 10⁵ t Cu from the Mt. Pinatubo magma. Our data suggest that intrusives having adakitic signatures are genetically associated with Au-Cu and Cu-Mo mineralization, auriferous porphyry copper deposits, and epithermal gold veins. High H₂O, Cl, Sr/Y, Pb/Ce, Mo/Ce, As/Ce and Sb/Ce in Mt. Pinatubo melts reflect the contribution of deep fluids derived from subducted sediments and altered MORBs in the dacite genesis. The slab-derived fluids carrying mobile elements are likely responsible for the enrichment of adakite magmas in gold, associated metals and H₂O, and may explain the exceptional ore-forming potential of adakite magmatism.     

The influence of bulk composition on the speciation of water in silicate glasses

Infrared spectroscopy was used to determine the concentrations of molecular water and hydroxyl groups in hydrous rhyolitic, orthoclasic, jadeitic, and Ca–Al-silicate glasses synthesized by quenching of melts from elevated presure and temperature. The rhyolitic glasses and some of the Ca–Al-silicate glasses were quenched from water-vapor-saturated melts and used to determine the solubility of water in melts of these compositions. For all compositions studied, hydroxyl groups are the dominant hydrous species at low total water contents, whereas molecular water dominates at elevated water contents. Although the trends in species concentrations in all these compositions are similar, the proportions of the two hydrous species are influenced by silicate chemistry: increasing silica content and K relative to Na both favor molecular water over hydroxyl. Results on rhyolitic glass demonstrate that molecular water is also favored by decreasing temperature at T<850°C. For rhyolitic glasses quenched from vapor-saturated melts, the mole fraction of molecular water is proportional to water fugacity for P(H₂O)1500 bars, demonstrating that the behavior of molecular water is approximately Henrian at total water contents up to at least several weight percent. Data on water solubility for albitic, orthoclasic, and Ca–Al-silicate melts to higher pressures can also be fit by assuming Henrian behavior for molecular water and can be used to set constraints on the partial molar volume of water in these melts. The demonstration of Henry's law for molecular water in these liquids provides a link between spectroscopic measurements of microscopic species concentrations and macroscopic thermodynamic properties.