Partitioning behavior of chlorine and fluorine in the system apatite-silicate melt-fluid

The partitioning behavior of Cl among apatite, mafic silicate melt, and aqueous fluid and of F between apatite and melt have been determined in experiments conducted at 1066 to 1150 °C and 199-205 MPa. The value of D_Cl^apatite/melt (wt. fraction of Cl in apatite/Cl in melt) ≈0.8 for silicate melt containing less than ∼3.8 wt.% Cl. At higher melt Cl contents, small increases in melt Cl concentration are accompanied by large increases in apatite Cl concentration, forcing D_Cl^apatite/melt to increase as well. Melt containing less than 3.8% Cl coexists with water-rich vapor; that containing more Cl coexists with saline fluid, the salinity of which increases rapidly with small increases in melt Cl content, analogous to the dependency of apatite composition on melt Cl content. This behavior is due to the fact that the solubility of Cl in silicate melt depends strongly on the composition of the melt, particularly its Mg, Ca, Fe, and Si contents. Once the melt becomes “saturated” in Cl, additional Cl must be accommodated by coexisting fluid, apatite, or other phases rather than the melt itself. Because Cl solubility depends on composition, the Cl concentration at which D_Cl^apatite/melt and D_Cl^fluid/melt begin to increase also depends on composition. The experiments reveal that D_F^apatite/melt ≈3.4. In contrast to Cl, the concentration of F in silicate melt is only weakly dependent on composition (mainly on melt Ca contents), so D_F^apatite/melt is constant for a wide range of composition.The experimental data demonstrate that the fluids present in the waning stages of the solidification of the Stillwater and Bushveld complexes were highly saline. The Cl-rich apatite in these bodies crystallized from interstitial melt with high Cl/(F + OH) ratio. The latter was generated by the combined processes of fractional crystallization and dehydration by its reaction with the relatively large mass of initially anhydrous pyroxene through which it percolated.     

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.     

S-rich apatite-hosted glass inclusions in xenoliths from La Palma: constraints on the volatile partitioning in evolved alkaline magmas

The composition of S-rich apatite, of volatile-rich glass inclusions in apatite, and of interstitial glasses in alkaline xenoliths from the 1949 basanite eruption in La Palma has been investigated to constrain the partitioning of volatiles between apatite and alkali-rich melts. The xenoliths are interpreted as cumulates from alkaline La Palma magmas. Apatite contains up to 0.89 wt% SO₃ (3560 ppm S), 0.31 wt% Cl, and 0.66 wt% Ce₂O₃. Sulfur is incorporated in apatite via several independent exchange reactions involving (P⁵⁺, Ca²⁺) vs. (S⁶⁺, Si⁴⁺, Na⁺, and Ce³⁺). The concentration of halogens in phonolitic to trachytic glasses ranges from 0.15 to 0.44 wt% for Cl and from <0.07 to 0.65 wt% for F. The sulfur concentration in the glasses ranges from 0.06 to 0.23 wt% SO₃ (sulfate-saturated systems). The chlorine partition coefficients (D_Cl^{apatite/glass}) range from 0.4 to 1.3 (average D_Cl^{apatite/glass} = 0.8), in good agreement with the results of experimental data in mafic and rhyolitic system with low Cl concentrations. With increasing F in glass inclusions D_F^{apatite/glass} decreases from 35 to 3. However, most of our data display a high partition coefficient (~30) close to D_F^{apatite/glass} determined experimentally in felsic rock. D_S^{apatite/glass} decreases from 9.1 to 2.9 with increasing SO₃ in glass inclusions. The combination of natural and experimental data reveals that the S partition coefficient tends toward a value of 2 for high S content in the glass (>0.2 wt% SO₃). D_S^{apatite/glass} is only slightly dependent on the melt composition and can be expressed as: SO_{3 apatite} (wt%) = 0.157 * ln SO_{3 glass} (wt%) + 0.9834. The phonolitic compositions of glass inclusions in amphibole and haüyne are very similar to evolved melts erupted on La Palma. The lower sulfur content and the higher Cl content in the phonolitic melt compared to basaltic magmas erupted in La Palma suggest that during magma evolution the crystallization of haüyne and pyrrhotite probably buffered the sulfur content of the melt, whereas the evolution of Cl concentration reflects an incompatible behavior. Trachytic compositions similar to those of the (water-rich) glass inclusions analyzed in apatite and clinopyroxene are not found as erupted products. These compositions are interpreted to be formed by the reaction between water-rich phonolitic melt and peridotite wall-rock.     

Experimental measurements of zircon/melt trace-element partition coefficients

Zircon was grown from trace-element doped hydrous peralkaline rhyolite melts with buffered oxygen fugacities in cold-seal experiments at 0.1 and 0.2 GPa and 800 °C and piston-cylinder experiments at 1.5 GPa and 900-1300 °C. Zircon and glass were present in all run products, and small monazite crystals were present in eight of the 12 experiments. Average diameters of zircon crystals ranged from 5 to 20 μm at 800 °C to 30-50 μm at 1300 °C. Zircon crystals have thin rims, and adjacent glass has a narrow (∼1 μm thick) compositional boundary layer. Concentrations obtained through in-situ analysis of cores of run product zircon crystals and melt pools were used to calculate trace-element partition coefficients D^zircon/melt for P, Sc, Ti, V, Y, La, Ce, Pr, Nd, Eu, Gd, Ho, Yb, Lu, Hf, Th, and U. In most cases Lu was the most (D 12-105) and La the least (0.06-0.95) compatible elements. D values from this study fall within the range of previously measured values for Rare Earth Elements (REE). However, D values measured experimentally show less fractionation than those recently measured using natural phenocryst/matrix pairs. For example, D_Lu/D_La measured experimentally in this study range between 27 and 206 compared to a value of 706,522 for a natural zircon/dacite pair [Sano, Y., Terada, K., and Fukuoka, T. 2002 High mass resolution ion microprobe analysis of rare earth elements in silicate glass, apatite and zircon: lack of matrix dependency. Chem. Geol.184, 217-230]. Although D values from this study show good agreement with the lattice strain model, D values from natural phenocryst/matrix pairs combined with measured zircon compositions better reproduce host-rock (magma) compositions of igneous rocks. They also yield more reasonable estimates of magma compositions when combined with compositions of ‘‘out-of-context” zircons. For example, compositions of the Hadean detrital zircons from Jack Hills, Australia yield LREE-enriched magmas when combined with D values from phenocryst/matrix pairs yields, but yield LREE-depleted magmas when experimentally determined D values are used. We infer that experimentally measured D^zircon/melt values represent disequilibrium partitioning resulting from rapid zircon growth during short laboratory timescales. Rapid growth causes development of observed diffusive boundary layers in the melt adjacent to zircon crystals. D values from phenocryst/matrix pairs are therefore recommended for petrogenetic modeling.     

The distribution of sulfur between haplogranitic melts and aqueous fluids

The distribution of sulfur between haplogranitic melt and aqueous fluid has been measured as a function of oxygen fugacity (Co-CoO-buffer to hematite-magnetite buffer), pressure (0.5-3 kbar), and temperature (750-850 °C). Sulfur always strongly partitions into the fluid. At a given oxygen fugacity, pressure and temperature, the distribution of sulfur between melt and fluid can be described by one constant partition coefficient over a wide range of sulfur concentrations. Oxygen fugacity is the most important parameter controlling sulfur partitioning. While the fluid/melt partition coefficient of sulfur is 468 ± 32 under Co-CoO buffer conditions at 2 kbar and 850 °C, it decreases to 47 ± 4 at an oxygen fugacity 0.5-1 log unit above Ni-NiO at the same pressure and temperature. A further increase in oxygen fugacity to the hematite-magnetite buffer has virtually no effect on the partition coefficient (D^fluid/melt = 49 ± 2). The dependence of D^fluid/melt on temperature and pressure was systematically explored at an oxygen fugacity 0.5-1 log units above Ni-NiO. At 850 °C, the effect of pressure on the partition coefficient is small (D^fluid/melt = 58 ± 3 at 0.5 kbar; 94 ± 9 at 1 kbar; 47 ± 4 at 2 kbar and 68 ± 5 at 3 kbar) and temperature also has only a minor effect on partitioning.The data show the “sulfur excess” observed in many explosive volcanic eruptions can easily be explained by the presence of a small fraction of hydrous fluid in the magma chamber before the eruption. The sulfur excess can be calculated as the product of the fluid/melt partition coefficient of sulfur and the mass ratio of fluid over melt in the erupted material. For a plausible fluid/melt partition coefficient of 47 under oxidizing conditions, a 10-fold sulfur excess corresponds to a 17.6 wt.% of fluid in the erupted material. Large sulfur excesses (10-fold or higher) are only to be expected if only a small fraction of the magma residing in the magma chamber is erupted.The behavior of sulfur, which seems to be largely independent of pressure and temperature under oxidizing conditions is very different from chlorine, where the fluid/melt partition coefficient strongly increases with pressure. Variations in the SO₂/HCl ratio of volcanic gases, if they reflect primary processes in the magma chamber, therefore provide an indicator of pressure variations in a magma. In particular, major increases in the S/Cl ratio of an aqueous fluid coexisting with a felsic magma suggest a pressure reduction in the magma chamber and/or magma rising to the surface.     

Mg-Fe partitioning between olivine and ultramafic melts at high pressures

Precise determination of the partitioning of Mg and Fe²⁺ between olivine and ultramafic melt has been made at pressures from 5 to 13 GPa using a MA-8 type multi-anvil high-pressure apparatus (PREM) installed at Earthquake Research Institute, University of Tokyo. A very short rhenium capsule (<100 μm sample thickness) was adopted to minimize temperature variation within the sample container. Synthetic gels with the composition of the upper mantle peridotite were used as starting materials to promote the homogeneity. Analyses of quenched melts and coexisting olivines were made with an electron probe microanalyzer. The obtained partition coefficient, K_D [=(FeO/MgO)^ol/(FeO/MgO)^melt], decreases from 0.35 to 0.25 with increasing pressure from 5 to 13 GPa, suggesting a negative correlation between pressure and K_D above 5 GPa. Our result is consistent with a parabolic relationship between K_D and degree of polymerization (NBO/T) of melts reported by previous studies at lower pressures. The negative correlation between pressure and K_D suggests that olivine crystallizing in a magma ocean becomes more Mg-rich with depth and that primary magmas generated in the upper mantle become more Fe-rich with depth than previously estimated.     

Solubility of Au in Cl- and S-bearing hydrous silicate melts

The solubility of Au in Cl- and S-bearing hydrous rhyodacitic and andesitic melts has been experimentally investigated at 1050 °C, 200 MPa and log fO₂ close to the Ni/NiO solid oxygen buffer (NNO). The concentrations of Au in the experimental glasses have been determined using Laser Ablation ICP-MS (LA) with special efforts to avoid incorporation of Au micronuggets in the analysis. It is concluded that metal micronuggets are an experimental artefact and produced by Au partitioning into the fluids during heating with consequent precipitation on fluid dissolution in the melting glass powder. Hence, the micronuggets do not represent quench phases and must be excluded from the analysis. The micro-analytical data obtained by LA show that Au concentrations vary from ∼0.2 to ∼2.5 ppm by weight, generally consistent with the literature data for other melt compositions. The measured Au concentrations increase with increasing amounts of Cl and S dissolved in the silicate melt and show a correlation with the apparent activities of Cl and S in the system. The apparent activities of Cl and S are defined by the simplified linear relationship between volatile concentrations in the melt and activity of volatiles. The maximum activity (a^∗ = 1) is assumed to be reached at the saturation of the systems in respect of Cl-rich brine or FeS liquid for Cl and S, respectively. The dependence of Au solubility on the concentrations/activities of Cl and S at the fixed redox conditions shows that Au may form not only oxide- but also Cl- and S-bearing complexes in silicate melts. Furthermore, it indicates that exsolution of S and Cl from the melt by degassing/segregation/crystallization processes may lead to mobilization and extraction of Au into the fluid, liquid and/or mineral phase(s).     

Determination of fluid/melt partition coefficients by LA-ICPMS analysis of co-existing fluid and silicate melt inclusions: Controls on element partitioning

Analyses of co-existing silicate melt and fluid inclusions, entrapped in quartz crystals in volatile saturated magmatic systems, allowed direct quantitative determination of fluid/melt partition coefficients. Investigations of various granitic systems (peralkaline to peraluminous in composition, log fO₂ = NNO−1.7 to NNO+4.5) exsolving fluids with various chlorinities (1-14 mol/kg) allowed us to assess the effect of these variables on the fluid/melt partition coefficients (D). Partition coefficients for Pb, Zn, Ag and Fe show a nearly linear increase with the chlorinity of these fluid (D_Pb ∼ 6 ∗ m_Cl, D_Zn ∼ 8 ∗ m_Cl, D_Ag ∼ 4 ∗ m_Cl, D_Fe ∼ 1.4 ∗ m_Cl, where m_Cl is the molinity of Cl). This suggests that these metals are dissolved primarily as Cl-complexes and neither oxygen fugacity nor the composition of the melt affects significantly their fluid/melt partitioning. By contrast, partition coefficients for Mo, B, As, Sb and Bi are highest in low salinity (1-2 mol/kg Cl) fluids with maximum values of D_Mo ∼ 20, D_B ∼ 15, D_As ∼ 13, D_Sb ∼ 8, D_Bi ∼ 15 indicating dissolution as non-chloride (e.g., hydroxy) complexes. Fluid/melt partition coefficients of copper are highly variable, but highest between vapor like fluids and silicate melt (D_Cu ? 2700), indicating an important role for ligands other than Cl. Partition coefficients for W generally increase with increasing chlorinity, but are exceptionally low in some of the studied brines which may indicate an effect of other parameters. Fluid/melt partition coefficients of Sn show a high variability but likely increase with the chlorinity of the fluid (D_Sn = 0.3-42, D_W = 0.8-60), and decrease with decreasing oxygen fugacity or melt peraluminosity.     

Experimental determination of F and Cl partitioning between lherzolite and basaltic melt

We experimentally determined F and Cl partition coefficients together with that of 19 trace elements (including REE, U-Th, HFSE and LILE) between basaltic melt and lherzolite minerals: olivine, orthopyroxene, clinopyroxene, plagioclase and garnet. Under conditions from 8 to 25 kbars and from 1,265 to 1,430°C, compatibilities of F and Cl are globally ordered as D ^Cpx/melt > D ^Opx/melt > D ^Grt/melt > D ^Ol/melt > D ^Plag/melt, and D _F ^mineral/melt is larger than D _Cl^mineral/melt. Four other major results were brought to light. (1) Chlorine partition coefficients positively correlate with the jadeite component in orthopyroxene, and increase of the CaTs component promotes Cl incorporation in clinopyroxene. (2) Variations of fluorine partition coefficients correlate strongly with melt viscosity. (3) F and Cl partition coefficients correlate with the Young’s modulus (E ₀) of pyroxene octahedral sites (M sites) and with Raman vibrational modes of pyroxenes. This demonstrates a fundamental interaction between the M site of pyroxenes and the incorporation of F and Cl. (4) We also determined the parameters of the lattice-strain model applied to 3+ cation trace elements for the two M sites in orthopyroxene and clinopyroxene: D ₀^M1, D ₀^M2, r ₀^M1, r ₀^M2, E ₀^M1 and E ₀^M2.     

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.     

Trace element partitioning in plagioclase feldspar

Compilation and interpretation of experimental and natural Nernst partition coefficient (^{plagioclase/melt}D) data show that, with a few exceptions, increases in ^{plagioclase/melt}D correlate with decreasing anorthite-content of plagioclase. In contrast, increases of ^{plagioclase/melt}D for Ga, Sc, Cu, Zn, Zr, Hf and Ti, are better correlated against decreasing melt MgO or increasing melt SiO₂ contents. ^{plagioclase/melt}D for Ti and the rare earth elements (REE) show little dependence on temperature, but increase as the melt water content increases. ^{plagioclase/melt}D for K and Sr are sensitive to pressure. Variations of D₀ (the strain compensated partition coefficient), r₀ (the size of the site into which REE substitute), and E (Young’s Modulus of this site) were parameterized against variations of melt SiO₂, the An-content of plagioclase, and other combinations of variables, allowing ^{plagioclase/melt}D_REE-Y to be calculated from a variety of input parameters. The interrelations of temperature, melt MgO and SiO₂ content, and plagioclase anorthite-content for wet and dry systems were also parameterized to facilitate interpolation where such data are lacking. When combined, these semi-empirical parameterizations yield ^{plagioclase/melt}D results comparable to available experimental and natural data.     

Melting of hydrous pyroxenites with alkali amphiboles in the continental mantle: 2. Trace element compositions of melts and minerals

The trace element compositions of melts and minerals from high-pressure experiments on hydrous pyroxenites containing K-richterite are presented. The experiments used mixtures of a third each of the natural minerals clinopyroxene, phlogopite and K-richterite, some with the addition of 5% of an accessory phase ilmenite, rutile or apatite. Although the major element compositions of melts resemble natural lamproites, the trace element contents of most trace elements from the three-mineral mixture are much lower than in lamproites. Apatite is required in the source to provide high abundances of the rare earth elements, and either rutile and/or ilmenite is required to provide the high field strength elements Ti, Nb, Ta, Zr and Hf. Phlogopite controls the high levels of Rb, Cs and Ba.Since abundances of trace elements in the various starting mixtures vary strongly because of the use of natural minerals, we calculated mineral/melt partition coefficients (D_Min/melt) using mineral modes and melting reactions and present trace element patterns for different degrees of partial melting of hydrous pyroxenites. Rb, Cs and Ba are compatible in phlogopite and the partition coefficient ratio phlogopite/K-richterite is high for Ba (1 3 6) and Rb (12). All melts have low contents of most of the first row transition elements, particularly Ni and Cu ((0.1–0.01) × primitive mantle). Nickel has high D_Min/melt for all the major minerals (12 for K-richterite, 9.2 for phlogopite and 5.6 for Cpx) and so behaves at least as compatibly as in melting of peridotites. Fluorine/chlorine ratios in melts are high and D_Min/melt for fluorine decreases in the order apatite (2.2) > phlogopite (1.5) > K-richterite (0.87). The requirement for apatite and at least one Ti-oxide in the source of natural lamproites holds for mica pyroxenites that lack K-richterite. The results are used to model isotopic ageing in hydrous pyroxenite source rocks: phlogopite controls Sr isotopes, so that lamproites with relatively low ⁸⁷Sr/⁸⁶Sr must come from phlogopite-poor source rocks, probably dominated by Cpx and K-richterite. At high pressures (>4 GPa), peritectic Cpx holds back Na, explaining the high K₂O/Na₂O of lamproites.     

Trace element partitioning between apatite and silicate melts

We present new experimental apatite/melt trace element partition coefficients for a large number of trace elements (Cs, Rb, Ba, La, Ce, Pr, Sm, Gd, Lu, Y, Sr, Zr, Hf, Nb, Ta, U, Pb, and Th). The experiments were conducted at pressures of 1.0 GPa and temperatures of 1250 °C. The rare earth elements (La, Ce, Pr, Sm, Gd, and Lu), Y, and Sr are compatible in apatite, whereas the larger lithophile elements (Cs, Rb, and Ba) are strongly incompatible. Other trace elements such as U, Th, and Pb have partition coefficients close to unity. In all experiments we found D_Hf > D_Zr, D_Ta ≈ D_Nb, and D_Ba > D_Rb > D_Cs. The experiments reveal a strong influence of melt composition on REE partition coefficients. With increasing polymerisation of the melt, apatite/melt partition coefficients for the rare earth elements increase for about an order of magnitude. We also present some results in fluorine-rich and water-rich systems, respectively, but no significant influence of either H₂O or F on the partitioning was found. Furthermore, we also present experimentally determined partition coefficients in close-to natural compositions which should be directly applicable to magmatic processes.     

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.     

Olivine/melt transition metal partitioning, melt composition, and melt structure—Melt polymerization and Q-speciation in alkaline earth silicate systems

The two most abundant network-modifying cations in magmatic liquids are Ca²⁺ and Mg²⁺. To evaluate the influence of melt structure on exchange of Ca²⁺ and Mg²⁺ with other geochemically important divalent cations (m-cations) between coexisting minerals and melts, high-temperature (1470-1650 °C), ambient-pressure (0.1 MPa) forsterite/melt partitioning experiments were carried out in the system Mg₂SiO₄-CaMgSi₂O₆-SiO₂ with ?1 wt% m-cations (Mn²⁺, Co²⁺, and Ni²⁺) substituting for Ca²⁺ and Mg²⁺. The bulk melt NBO/Si-range (NBO/Si: nonbridging oxygen per silicon) of melt in equilibrium with forsterite was between 1.89 and 2.74. In this NBO/Si-range, the NBO/Si(Ca) (fraction of nonbridging oxygens, NBO, that form bonds with Ca²⁺, Ca²⁺-NBO) is linearly related to NBO/Si, whereas fraction of Mg²⁺-NBO bonds is essentially independent of NBO/Si. For individual m-cations, rate of change of K_D(m−Mg) with NBO/Si(Ca) for the exchange equilibrium, m^melt + Mg^olivine ? m^olivine + Mg^melt, is linear. K_D(m−Mg) decreases as an exponential function of increasing ionic potential, Z/r² (Z: formal electrical charge, r: ionic radius—here calculated with oxygen in sixfold coordination around the divalent cations) of the m-cation. The enthalpy change of the exchange equilibrium, ΔH, decreases linearly with increasing Z/r² [ΔH = 261(9)-81(3)·Z/r² (Å⁻²)]. From existing information on (Ca,Mg)O-SiO₂ melt structure at ambient pressure, these relationships are understood by considering the exchange of divalent cations that form bonds with nonbridging oxygen in individual Qⁿ-species in the melts. The negative ∂K_D(m−Mg)/∂(Z/r²) and ∂(ΔH)/∂(Z/r²) is because increasing Z/r² is because the cations forming bonds with nonbridging oxygen in increasingly depolymerized Qⁿ-species where steric hindrance is decreasingly important. In other words, principles of ionic size/site mismatch commonly observed for trace and minor elements in crystals, also govern their solubility behavior in silicate melts.     

Magmatic Evolution of Changbaishan Tianchi Volcano,China–North Korea: Evidence from Mineral-Hosted Melt and Fluid Inclusions

Data on mineral-hosted melt, fluid, and crystalline inclusions were used to study the composition and evolution of melts that produced rocks of Changbaishan Tianchi volcano, China–North Korea, and estimate their crystallization parameters. The melts crystallized within broad ranges of temperature (1220–700°C) and pressure (3100–1000 bar), at a drastic change in the redox potential: Δ log \(f_{O_2}\) from NNO + 0.92 to +1.42 for the basalt melts, NNO –1.61 to –2.09 for the trachybasaltic andesite melts, NNO –2.63 to –1.89 for the comendite melts, and NNO –1.55 to –3.15 for the pantellerite melts. The paper reports estimates of the compositions of melts that produced the continuous rock series from trachybasalt to comendite and pantellerite. In terms of trace-element concentrations, all of the mafic melts are comparable with OIB magmas. The silicic melts are strongly enriched in trace elements and REE. The most strongly enriched melts contain concentrations of certain elements almost as high as in ores of these elements. The paper reports data on H2O concentrations in melts of different composition. It is demonstrated that the variations in the H₂O concentrations were controlled by magma degassing. Data are reported on the Sr and Nd composition of the rocks. The deviations in the Sr isotopic composition are proportional to the ⁸⁷Sr/⁸⁶Sr ratio and could be produced in a melt with a high enough ⁸⁷Sr/⁸⁶Sr ratio during a geologically fairly brief time period. The evolution of melts that produced rocks of the volcano was controlled by crystallization differentiation of the parental basalt magmas at insignificant involvement of melt mixing and liquid immiscibility of silicate and sulfide melts. The alkaline salic rocks were generated in shallow-sitting (13–3.5 km) magmatic chambers in which the melts underwent profound differentiation that gave rise to pantellerites and comendites strongly enriched in trace elements (Th, Nb, Ta, Zr, and REE). Data on the composition of the magmas and parameters of their derivation are used to develop a generalized petrologic–geodynamic model for the origin of Changbaishan Tianchi volcano.     

Structure of H2O-saturated peralkaline aluminosilicate melt and coexisting aluminosilicate-saturated aqueous fluid determined in-situ to 800 °C and ∼800 MPa

The structure of H₂O-saturated silicate melts and of silicate-saturated aqueous solutions, as well as that of supercritical silicate-rich aqueous liquids, has been characterized in-situ while the sample was at high temperature (to 800 °C) and pressure (up to 796 MPa). Structural information was obtained with confocal microRaman and with FTIR spectroscopy. Two Al-bearing glasses compositionally along the join Na₂O•4SiO₂-Na₂O•4(NaAl)O₂-H₂O (5 and 10 mol% Al₂O₃, denoted NA5 and NA10) were used as starting materials. Fluids and melts were examined along pressure-temperature trajectories of isochores of H₂O at nominal densities (from PVT properties of pure H₂O) of 0.85 g/cm³ (NA10 experiments) and 0.86 g/cm³ (NA5 experiments) with the aluminosilicate + H₂O sample contained in an externally-heated, Ir-gasketed hydrothermal diamond anvil cell.Molecular H₂O (H₂O°) and OH groups that form bonds with cations exist in all three phases. The OH/H₂O° ratio is positively correlated with temperature and pressure (and, therefore, fugacity of H₂O, f_H2O) with (OH/H₂O°)^melt > (OH/H₂O°)^fluid at all pressures and temperatures. Structural units of Q³, Q², Q¹, and Q⁰ type occur together in fluids, in melts, and, when outside the two-phase melt + fluid boundary, in single-phase liquids. The abundance of Q⁰ and Q¹ increases and Q² and Q³ decrease with f_H2O. Therefore, the NBO/T (nonbridging oxygen per tetrahedrally coordination cations), of melt is a positive function of f_H2O. The NBO/T of silicate in coexisting aqueous fluid, although greater than in melt, is less sensitive to f_H2O.The melt structural data are used to describe relationships between activity of H₂O and melting phase relations of silicate systems at high pressure and temperature. The data were also combined with available partial molar configurational heat capacity of Qⁿ-species in melts to illustrate how these quantities can be employed to estimate relationships between heat capacity of melts and their H₂O content.     

Crystal-melt partitioning of noble gases (helium, neon, argon, krypton, and xenon) for olivine and clinopyroxene

Mineral-melt partition coefficients of all noble gases (^min/meltD_i) have been obtained for olivine (ol) and clinopyroxene (cpx) by UV laser ablation (213 nm) of individual crystals grown from melts at 0.1 GPa mixed noble gas pressure. Experimental techniques were developed to grow crystals virtually free of melt and fluid inclusions since both have been found to cause profound problems in previous work. This is a particularly important issue for the analysis of noble gases in crystals that have very low partition coefficients relative to coexisting melt and fluid phases. The preferred partitioning values obtained for the ol-melt system for He, Ne, Ar, Kr, and Xe are 0.00017(13), 0.00007(7), 0.0011(6), 0.00026(16), and , respectively. The respective cpx-melt partition coefficients are 0.0002(2), 0.00041(35), 0.0011(7), 0.0002(2), and . The data confirm the incompatible behaviour of noble gases for both olivine and clinopyroxene but unlike other trace elements these values show little variation for a wide range of atomic radius. The lack of dependence of partitioning on atomic radius is, however, consistent with the partitioning behaviour of other trace elements which have been found to exhibit progressively lower dependence of ^min/meltD_i on radius as the charge decreases. As all noble gases appear to exhibit similar ^min/meltD_i values we deduce that noble gases are not significantly fractionated from each other by olivine and clinopyroxene during melting and fractional crystallisation. Although incompatible, the partitioning values for noble gases also suggest that significant amounts of primordial noble gases may well have been retained in the mantle despite intensive melting processes. The implication of our data is that high primordial/radiogenic noble gas ratios (³He/⁴He, ²²Ne/²¹Ne, and ³⁶Ar/⁴⁰Ar) characteristic of plume basalt sources can be achieved by recycling a previously melted (depleted) mantle source rather than reflecting an isolated, non-degassed primordial mantle region.     

Sulfur partition coefficient between apatite and rhyolite: the role of bulk S content

Experimental studies have been performed to constrain sulfur behavior during apatite crystallization and to determine sulfur partition coefficient between apatite and melt (Kd_S^apatite/melt) at oxidizing conditions. Crystallization experiments have been conducted with a hydrous rhyolitic melt and different bulk sulfur contents (0.15 to 2 wt.% S) at 900 and 1,000°C, 200 MPa and Δlog =NNO+3.6. The sulfur content in the glass increases with increasing amount of added S. Anhydrite crystallizes for S added = 0.75 wt.% (0.10 and 0.13 wt.% SO₃ in glass at 900 and 1,000°C, respectively). The amount of anhydrite increases and the amount of apatite decreases with increasing amount of added sulfur. The sulfur exchange reaction in apatite is influenced by the bulk composition of the melt (e.g., P content). However, changing melt composition has only little effect on Kd_S^apatite/melt for the investigated rhyolitic composition. The Kd_S^apatite/melt does not depend directly on temperature, decreases from 14.2 to 2.7 with increasing S content in glass from SO₃=0.03 to 0.19 wt.%, respectively, and can be predicted from the following equation: ln Kd = −0.0025×S in melt (in ppm)+2.9178. The combination of experimental data obtained for rhyolitic and andesitic melts reveals that the sulfur partition coefficient tends toward a value of 2 for high-sulfur content in the glass (>0.2 wt.% SO₃). Using S in apatite as proxy for determining S content in melt is promising but additional experimental data are needed to clarify the individual effects of T, , and P and Ca content in the melt on S partitioning.     

Partitioning of Ru, Rh, Pd, Re, Ir, and Au between Cr-bearing spinel, olivine, pyroxene and silicate melts

A series of high temperature experiments was undertaken to study partitioning of several highly siderophile elements (HSE; Ru, Rh, Pd, Re, Os, Ir, Pt and Au) between Cr-rich spinel, olivine, pyroxene and silicate melt. Runs were carried out on a Hawaiian ankaramite, a synthetic eucrite basalt, and a DiAn eutectic melt, at one bar, 19 kbar, and 20 kbar, respectively, in the temperature range of 1200 to 1300°C, at oxygen fugacities between the nickel-nickel oxide (NNO) and hematite-magnetite (HM) oxygen buffers. High oxygen fugacities were used to suppress the formation of HSE-rich “nuggets” in the silicate melts. The resulting oxide and silicate crystals (<100 μm) were analyzed using both SIMS and LA-ICP-MS, with a spatial resolution of 15 to 50 μm. Rhenium, Au and Pd were all found to be incompatible in Cr-rich spinel (D_Re^sp/melt = 0.0012-0.21, D_Au^sp/melt = 0.076, D_Pd^sp/melt = 0.14), whereas Rh, Ru and Ir were all found to be highly compatible (D_Rh^sp/melt = 41-530, D_Ru^sp/melt = 76-1143, D_Ir^sp/melt = 5-22000). Rhenium, Pd, Au and Ru were all found to be incompatible in olivine (D_Re^oliv/melt = 0.017-0.073, D_Pd^oliv/melt = 0.12, D_Au^oliv/melt = 0.12, D_Ru^oliv/melt = 0.23), Re is incompatible in orthopyroxene and clinopyroxene (D_Re^opx/melt = 0.013, D_Re^cpx/melt = 0.18-0.21), and Pt is compatible in clinopyroxene (D_Pt^cpx/melt = 1.5). The results are compared to and combined with previous work on HSE partitioning among spinel-structured oxides, and applied to some natural magmatic suites to demonstrate consistency.