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.     

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.     

Speciation of reduced C-O-H volatiles in coexisting fluids and silicate melts determined in-situ to ∼1.4 GPa and 800 °C

The structure of silicate melts in the system Na₂O·4SiO₂ saturated with reduced C-O-H volatile components and of coexisting silicate-saturated C-O-H solutions has been determined in a hydrothermal diamond anvil cell (HDAC) by using confocal microRaman and FTIR spectroscopy as structural probes. The experiments were conducted in-situ with the melt and fluid at high temperature (up to 800 °C) and pressure (up to 1435 MPa). Redox conditions in the HDAC were controlled with the reaction, Mo + H₂O = MoO₂ + H₂, which is slightly more reducing than the Fe + H₂O = FeO + H₂ buffer at 800 °C and less.The dominant species in the fluid are CH₄ + H₂O together with minor amounts of molecular H₂ and an undersaturated hydrocarbon species. In coexisting melt, CH₃ - groups linked to the silicate melt structure via Si-O-CH₃ bonding may dominate and possibly coexists with molecular CH₄. The abundance ratio of CH₃ - groups in melts relative to CH₄ in fluids increases from 0.01 to 0.07 between 500 and 800 °C. Carbon-bearing species in melts were not detected at temperatures and pressures below 400 °C and 730 MPa, respectively. A schematic solution mechanism is, Si-O-Si + CH₄?Si-O-CH₃+H-O-Si. This mechanism causes depolymerization of silicate melts. Solution of reduced (C-O-H) components will, therefore, affect melt properties in a manner resembling dissolved H₂O.     

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.     

Partitioning of F between H2O and CO2 fluids and topaz rhyolite melt

Fluid/melt distribution coefficients for F have been determined in experiments conducted with peraluminous topaz rhyolite melts and fluids consisting of H₂O and H₂O+CO₂ at pressures of 0.5 to 5 kbar, temperatures of 775°–1000°C, and concentrations of F in the melt ranging from 0.5 to 6.9 wt%. The major element, F, and Cl concentrations of the starting material and run product glasses were determined by electron microprobe, and the concentration of F in the fluid was calculated by mass balance. The H₂O concentrations of some run product glasses were determined by ion microprobe (SIMS). The solubility of melt in the fluid phase increases with increasing F in the system; the solubility of H₂O in the melt is independent of the F concentration of the system with up to 6.3 wt% F in the melt. No evidence of immiscible silica- and fluoriderich liquids was detected in the hydrous but water-undersaturated starting material glasses (8.5 wt% F in melt) or in the water-saturated run product glasses. F concentrates in topaz rhyolite melts relative to coexisting fluids at most conditions studied; however, D_F (wt% F in fluid/wt% F in melt) increases strongly with increasing F in the system. Maximum values of D_F in this study are significantly larger than those previously reported in the literature. Linear extrapolation of the data suggests that D_F is greater than one for water-saturated, peraluminous granitic melts containing 8 wt% F at 800° C and 2 kbar. D_F increases as temperature and as (H₂O/H₂O+CO₂) of the fluid increase. For topaz rhyolite melts containing 1 wt% F and with H₂O-rich fluids, D_F is independent of changes in pressure from 2 to 5 kbar at 800° C; for melts containing 1 wt% F and in equilibrium with CO₂-bearing fluids the concentrations of F in fluid increases with increasing pressure. F-and lithophile element-enriched granites may evolve to compositions containing extreme concentrations of F during the final stages of crystallization. If F in the melt exceeds 8 wt%, D_F is greater than one and the associated magmatic-hydrothermal fluid contains >4 molal F. Such F-enriched fluids may be important in the mass transport of ore constituents, i.e., F, Mo, W, Sn, Li, Be, Rb, Cs, U, Th, Nb, Ta, and B, from the magma.     

Partitioning behavior of chlorine and fluorine in the system apatite-melt-fluid. II: Felsic silicate systems at 200 MPa

Hydrothermal experiments were conducted to determine the partitioning of Cl between rhyolitic to rhyodacitic melts, apatite, and aqueous fluid(s) and the partitioning of F between apatite and these melts at ca. 200 MPa and 900-924 °C. The number of fluid phases in our experiments is unknown; they may have involved a single fluid or vapor plus saline liquid. The partitioning behavior of Cl between apatite and melt is non-Nernstian and is a complex function of melt composition and the Cl concentration of the system. Values of D_Cl^apat/melt (wt. fraction of: Cl in apatite/Cl in melt) vary from 1 to 4.5 and are largest when the Cl concentrations of the melt are at or near the Cl-saturation value of the melt. The Cl-saturation concentrations of silicate melts are lowest in evolved, silica-rich melts, so with elevated Cl concentrations in a system and with all else equal, the maximum values of D_Cl^apat/melt occur with the most felsic melt. In contrast, values of D_F^apat/melt range from 11 to 40 for these felsic melts, and many of these are an order of magnitude greater than those applying to basaltic melts at 200 MPa and 1066-1150 °C. The Cl concentration of apatite is a simple and linear function of the concentration of Cl in fluid. Values of D_Cl^fluid/apat for these experiments range from 9 to 43, and some values are an order of magnitude greater than those determined in 200-MPa experiments involving basaltic melts at 1066-1150 °C.In order to determine the concentrations and interpret the behavior of volatile components in magmas, the experimental data have been applied to the halogen concentrations of apatite grains from chemically evolved rocks of Augustine volcano, Alaska; Krakatau volcano, Indonesia; Mt. Pinatubo, Philippines; Mt. St. Helens, Washington; Mt. Mazama, Oregon; Lascar volcano, Chile; Santorini volcano, Greece, and the Bishop Tuff, California. The F concentrations of these magmas estimated from apatite-melt equilibria range from 0.06 to 0.12 wt% and are generally equivalent to the concentrations of F determined in the melt inclusions. In contrast, the Cl concentrations of the magmas estimated from apatite-melt equilibria (e.g., ca. 0.3-0.9 wt%) greatly exceed those determined in the melt inclusions from all of these volcanic systems except for the Bishop Tuff where the agreement is good. This discrepancy in estimated Cl concentrations of melt could result from several processes, including the hypothesis that the composition of apatite represents a comparatively Cl-enriched stage of magma evolution that precedes melt inclusion entrapment prior to the sequestration of Cl by coexisting magmatic aqueous and/or saline fluid(s).     

The behavior of trace elements during the chemical evolution of the H2O-, B-, and F-rich granite–pegmatite–hydrothermal system at Ehrenfriedersdorf, Germany: a SXRF study of melt and fluid inclusions

Detailed melt and fluid inclusion studies in quartz hosts from the Variscan Ehrenfriedersdorf complex revealed that ongoing fractional crystallization of the highly evolved H₂O-, B-, and F-rich granite magma produced a pegmatite melt, which started to separate into two immiscible phases at about 720°C, 100 MPa. With cooling and further chemical evolution, the immiscibilty field expanded. Two conjugate melts, a peraluminous one and a peralkaline one, coexisted down to temperatures of about 490°C. Additionally, high-salinity brine exsolved throughout the pegmatitic stage, along with low-density vapor. Towards lower temperatures, a hydrothermal system gradually developed. Boiling processes occurred between 450 and 400°C, increasing the salinities of hydrothermal fluids at this stage. Below, the late hydrothermal stage is dominated by low-salinity fluids. Using a combination of synchrotron radiation-induced X-ray fluorescence analysis and Raman spectroscopy, the concentration of trace elements (Mn, Fe, Zn, As, Sb, Rb, Cs, Sr, Zr, Nb, Ta, Ag, Sn, Ta, W, rare earth elements (REE), and Cu) was determined in 52 melt and 8 fluid inclusions that are representative of distinct stages from 720°C down to 380°C. Homogenization temperatures and water contents of both melt and fluid inclusions are used to estimate trapping temperatures, thus revealing the evolutionary stage during the process. Trace elements are partitioned in different proportions between the two pegmatite melts, high-salinity brines and exsolving vapors. Concentrations are strongly shifted by co ncomitant crystallization and precipitation of ore-forming minerals. For example, pegmatite melts at the initial stage (700°C) have about 1,600 ppm of Sn. Concentrations in both melts decrease towards lower temperatures due to the crystallization of cassiterite between 650 and 550°C. Tin is preferentially fractionated into the peralkaline melt by a factor of 2–3. While the last pegmatite melts are low in Sn (64 ppm at 500°C), early hydrothermal fluids become again enriched with about 800 ppm of Sn at the boiling stage. A sudden drop in late hydrothermal fluids (23 ppm of Sn at 370°C) results from precipitation of another cassiterite generation between 400 and 370°C. Zinc concentrations in peraluminous melts are low (some tens of parts per million) and are not correlated with temperature. In coexisting peralkaline melts and high-T brines, they are higher by a factor of 2–3. Zinc continuously increases in hydrothermal fluids (3,000 ppm at 400°C), where the precipitation of sphalerite starts. The main removal of Zn from the fluid system occurs at lower temperatures. Similarly, melt and fluid inclusion concentrations of many other trace elements directly reflect the crystallization and precipitation history of minerals at distinctive temperatures or temperature windows.     

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.     

SIMS determination of trace element partition coefficients between garnet, clinopyroxene and hydrous basaltic liquids at 2–7.5 GPa and 1080–1200°C

Trace element partition coefficients (D's) for up to 13 REE, Nb, Ta, Zr, Hf, Sr and Y have been determined by SIMS analysis of seven garnets, four clinopyroxenes, one orthopyroxene and one phlogopite crystallized from an undoped basanite and a lightly doped (200 ppm Nb, Ta and Hf) quartz tholeiite. Experiments were conducted at 2–7.5 GPa, achieving near-liquidus crystallization at relatively low temperatures of 1080–1200°C under strongly hydrous conditions (5–27 wt.% added water). Garnet and pyroxene D_REE show a parabolic pattern when plotted against ionic radius, and conform closely to the lattice strain model of Blundy and Wood (Blundy, J.D., Wood, B.J., 1994. Prediction of crystal–melt partition coefficients from elastic moduli. Nature 372, 452–454). Comparison, at constant pressure, between hydrous and anhydrous values of the strain-free partition coefficient (D₀) for the large cation sites of garnet and clinopyroxene reveals the relative importance of temperature and melt water content on partitioning. In the case of garnet, the effect of lower temperature, which serves to increase D₀, and higher water content, which serves to decrease D₀, counteract each other to the extent that water has little effect on garnet–melt D₀ values. In contrast, the effect of water on clinopyroxene–melt D₀ overwhelms the effect of temperature, such that D₀ is significantly lower under hydrous conditions. For both minerals, however, the lower temperature of the hydrous experiments tends to tighten the partitioning parabolas, increasing fractionation of light from heavy REE compared to anhydrous experiments.

Three sets of near-liquidus clinopyroxene–garnet two-mineral D values increase the range of published experimental determinations, but show significant differences from natural two-mineral D's determined for subsolidus mineral pairs. Similar behaviour is observed for the first experimental data for orthopyroxene–clinopyroxene two-mineral D's when compared with natural data. These differences are in large part of a consequence of the subsolidus equilibration temperatures and compositions of natural mineral pairs. Great care should therefore be taken when using natural mineral–mineral partition coefficients to interpret magmatic processes.

The new data for strongly hydrous compositions suggest that fractionation of Zr–Hf–Sm by garnet decreases with increasing depth. Thus, melts leaving a garnet-dominated residuum at depths of about 200 km or greater may preserve source Zr/Hf and Hf/Sm. This contrasts with melting at shallower depths where both garnet and clinopyroxene will cause Zr–Hf–Sm fractionation. Also, at shallower depths, clinopyroxene-dominated fractionation may produce a positive Sr spike in melts from spinel lherzolite, but for garnet lherzolite melting, no Sr spike will result. Conversely, clinopyroxene megacrysts with negative Sr spikes may crystallize from magmas without anomalous Sr contents when plotted on mantle compatibility diagrams. Because the characteristics of strongly hydrous silicate melt and solute-rich aqueous fluid converge at high pressure, the hydrous data presented here are particularly pertinent to modelling processes in subduction zones, where aqueous fluids may have an important metasomatic role.     

C�CO�CH�CS fluids and granitic magma: how S partitions and modifies CO2 concentrations of fluid-saturated felsic melt at 200?MPa

Hydrothermal volatile-solubility and partitioning experiments were conducted with fluid-saturated haplogranitic melt, H₂O, CO₂, and S in an internally heated pressure vessel at 900°C and 200?MPa; three additional experiments were conducted with iron-bearing melt. The run-product glasses were analyzed by electron microprobe, FTIR, and SIMS; and they contain ??0.12 wt% S, ??0.097 wt% CO₂, and ??6.4 wt% H₂O. Apparent values of log f _O2 for the experiments at run conditions were computed from the [(S⁶⁺)/(S⁶⁺+S^2?)] ratio of the glasses, and they range from NNO ?0.4 to NNO?+?1.4. The C?CO?CH?CS fluid compositions at run conditions were computed by mass balance, and they contained 22?C99?mol% H₂O, 0?C78?mol% CO₂, 0?C12?mol% S, and <3 wt% alkalis. Eight S-free experiments were conducted to determine the H₂O and CO₂ concentrations of melt and fluid compositions and to compare them with prior experimental results for C?CO?CH fluid-saturated rhyolite melt, and the agreement is excellent. Sulfur partitions very strongly in favor of fluid in all experiments, and the presence of S modifies the fluid compositions, and hence, the CO₂ solubilities in coexisting felsic melt. The square of the mole fraction of H₂O in melt increases in a linear fashion, from 0.05 to 0.25, with the H₂O concentration of the fluid. The mole fraction of CO₂ in melt increases linearly, from 0.0003 to 0.0045, with the CO₂ concentration of C?CO?CH?CS fluids. Interestingly, the CO₂ concentration in melts, involving relatively reduced runs (log f _O2????NNO?+?0.3) that contain 2.5?C7?mol% S in the fluid, decreases significantly with increasing S in the system. This response to the changing fluid composition causes the H₂O and CO₂ solubility curve for C?CO?CH?CS fluid-saturated haplogranitic melts at 200?MPa to shift to values near that modeled for C?CO?CH fluid-saturated, S-free rhyolite melt at 150?MPa. The concentration of S in haplogranitic melt increases in a linear fashion with increasing S in C?CO?CH?CS fluids, but these data show significant dispersion that likely reflects the strong influence of f _O2 on S speciation in melt and fluid. Importantly, the partitioning of S between fluid and melt does not vary with the (H₂O/H₂O?+?CO₂) ratio of the fluid. The fluid-melt partition coefficients for H₂O, CO₂, and S and the atomic (C/S) ratios of the run-product fluids are virtually identical to thermodynamic constraints on volatile partitioning and the H, S, and C contents of pre-eruptive magmatic fluids and volcanic gases for subduction-related magmatic systems thus confirming our experiments are relevant to natural eruptive systems.     

Partitioning of lanthanides and Y between immiscible silicate and fluoride melts, fluorite and cryolite and the origin of the lanthanide tetrad effect in igneous rocks

Some F-rich granitic rocks show anomalous, nonchondritic ratios of Y/Ho, extreme negative Eu anomalies, and unusual, discontinuous, segmented chondrite-normalised plots of rare earth elements (REE). The effects of F-rich fluids have been proposed as one of the explanations for the geochemical anomalies in the evolved granitic systems, as the stability of nonsilicate complexes of individual rare earths may affect the fluid-melt element partitioning. The lanthanide tetrad effect, related to different configurations of 4f-electron subshells of the lanthanide elements, is one of the factors affecting such complexing behaviour. We present the first experimental demonstration of the decoupling of Y and Ho, and the tetrad effect in the partitioning of rare earths between immiscible silicate and fluoride melts. Two types of experiments were performed: dry runs at atmospheric pressure in a high-temperature centrifuge at 1100 to 1200°C, and experiments with the addition of H₂O at 700 to 800°C and 100 MPa in rapid-quench cold-seal pressure vessels. Run products were analysed by electron microprobe (major components), solution-based inductively coupled plasma mass spectrometry (ICP-MS) (REE in the centrifuged runs), and laser ablation ICP-MS (REE and Li in the products of rapid-quench runs). All the dry centrifuge runs were performed at super-liquidus, two-phase conditions. In the experiments with water-bearing mixtures, minor amounts of aqueous vapour were present in addition to the melts. We found that lanthanides and Y concentrated strongly in the fluoride liquids, with two-melt partition coefficients reaching values as high as 100-220 in water-bearing compositions. In all the experimental samples, two-melt partition coefficients of lanthanides show subtle periodicity consistent with the tetrad effect, and the partition coefficient of Y is greater than that of Ho. One of the mixtures also produced abundant fluorite (CaF₂) and cryolite (Na₃AlF₆) crystals, which enabled us to study fluorite-melt and cryolite-melt REE partitioning. REE concentrations in fluorite are high and comparable to those in the fluoride melt. However, fluorite-melt partition coefficients appear to depend mostly on ionic radii and show neither significant tetrad anomalies, nor differences in Y and Ho partitioning. In contrast, REE concentrations in cryolite are low (∼5-10 times lower than in the silicate melt), and cryolite-melt REE partitioning shows very strong tetrad and Y-Ho anomalies. Our results imply that Y-Ho and lanthanide tetrad anomalies are likely to be caused mainly by aluminofluoride complexes, and the tetrad REE patterns in natural igneous rocks can result from fractionation of F-rich magmatic fluids.     

Determination of melt influence on divalent element partitioning between anorthite and CMAS melts

We propose a theory for crystal-melt trace element partitioning that considers the energetic consequences of crystal-lattice strain, of multi-component major-element silicate liquid mixing, and of trace-element activity coefficients in melts. We demonstrate application of the theory using newly determined partition coefficients for Ca, Mg, Sr, and Ba between pure anorthite and seven CMAS liquid compositions at 1330 °C and 1 atm. By selecting a range of melt compositions in equilibrium with a common crystal composition at equal liquidus temperature and pressure, we have isolated the contribution of melt composition to divalent trace element partitioning in this simple system. The partitioning data are fit to Onuma curves with parameterizations that can be thermodynamically rationalized in terms of the melt major element activity product (a_Al2O3)(a_SiO2)² and lattice strain theory modeling. Residuals between observed partition coefficients and the lattice strain plus major oxide melt activity model are then attributed to non-ideality of trace constituents in the liquids. The activity coefficients of the trace species in the melt are found to vary systematically with composition. Accounting for the major and trace element thermodynamics in the melt allows a good fit in which the parameters of the crystal-lattice strain model are independent of melt composition.     

Solubility of water-hydrogen fluid in haplogranite,albite, and sodium disilicate melts

Solubility curves of water-hydrogen fluid were studied using a high-pressure gas apparatus at a pressure of 200 MPa under variable fluid composition in haplogranite (Ab ₃₉ Or ₃₂ Qtz ₂₉, 950°C), Na-disilicate (Na₂Si₂O₅, 950°C), and albite melts (1200°C). The mole fraction of hydrogen in experiments was controlled directly by Ar-H₂ mixtures using a specially designed cell with a Shaw membrane. $ X_{H_2 }^{Ar - H_2 } $ X_{H_2 }^{Ar - H_2 } ranged from 0 to 1. In some experiments with haplogranite and Na-disilicate melts under oxidizing conditions, in order to increase the accuracy of experimental parameters, the fugacities of oxygen and hydrogen were controlled using the double-capsule technique and the solid-phase buffer mixtures Ni-NiO (NNO) and Co-CoO (CCO). The addition of H₂ to the H₂O-saturated systems ($ X_{H_2 }^{H_2 O - H_2 } $ X_{H_2 }^{H_2 O - H_2 } ≥ 0.012) results in the appearance of a distinct maximum on the solubility curves at $ X_{H_2 }^{H_2 O - H_2 } $ X_{H_2 }^{H_2 O - H_2 } = 0.05–0.07 (H₂ mole fractions were calculated for real H₂O-H₂ mixtures of real gases), and the maximum content of H₂O-H₂ fluid increases relative to the H₂O-saturated melts by 1.51 wt % for haplogranite melt at $ X_{H_2 } $ X_{H_2 } = 0.063, 2.68 wt % for albite melt at $ X_{H_2 } $ X_{H_2 } = 0.066, and 3.54 wt % for Na-disilicate melt at $ X_{H_2 } $ X_{H_2 } = 0.067. A further increase in H₂ content in the gas mixture decreases the solubility of H₂O-H₂ fluid in the melts, and under pure H₂ pressure, the contents of fluid components are 0.08 wt % in haplogranite melt and 0.06 wt % in albite melt. The ¹H NMR study of aluminosilicate and Na-silicate glasses obtained under the pressure of H₂O and H₂O-H₂ fluids suggests different mechanisms of the dissolution of H₂O and H₂O-H₂ fluids in magmatic melts. In addition to the spectra of dissolved water fluid, the spectra of quenched glasses synthesized under H₂O-H₂ fluid pressure exhibited a narrow line of molecular hydrogen with a width at half height of 1.8–2.0 kHz at $ X_{H_2 } $ X_{H_2 } ≥ 0.653 for albite and $ X_{H_2 } $ X_{H_2 } ≥ 0.063 for Na-disilicate and two lines at $ X_{H_2 } $ X_{H_2 } ≥ 0.063 for the haplogranite composition.     

Aqueous fluids and hydrous melts in high-pressure and ultra-high pressure rocks: Implications for element transfer in subduction zones

High-pressure (HP) and ultra-high pressure (UHP) terranes are excellent natural laboratories to study subduction-zone processes. In this paper we give a brief theoretical background and we review experimental data and observations in natural rocks that constrain the nature and composition of the fluid phase present in HP and UHP rocks. We argue that a fluid buffered by a solid residue is compositionally well defined and is either an aqueous fluid (total amount of dissolved solids < 30 wt.%) or a hydrous melt (H₂O < 35 wt.%). There is only a small temperature range of approximately 50–100 °C, where transitional solute-rich fluids exist. A review of available experimental data suggest that in felsic rocks the second critical endpoint is situated at 25–35 kbar and  700 °C and hence must be considered in the study of UHP rocks. Despite this, the nature of the fluid phase can be constrained by relating the peak metamorphic conditions of rocks to the position of the wet solidus even if the peak pressure exceeds the pressure where the wet solidus terminates at the second critical endpoint. Transitional solute-rich fluids are expected in UHP terrains (P > 30 kbar) with peak temperatures of about 700 ± 50 °C. At higher temperatures, hydrous granitic melts occur whereas at lower temperatures aqueous fluids coexists with eclogite-facies minerals. This argument is complemented by evidence on the nature of the fluid phase from high-pressure terrains. We show that in the diamond-bearing, high-temperature UHP rocks from the Kokchetav Massif there are not only hydrous felsic melts, but probably also carbonate and sulfide melts present.

Hydrous quartzo-feldspathic melts are mainly produced in high temperature UHP rocks and their composition is relatively well constrained from experiments and natural rocks. In contrast, constraining the composition of aqueous fluids is more problematic. The combined evidence from experiments and natural rocks indicates that aqueous fluids liberated at the blueschist to eclogite facies transition are dilute. They contain only moderate amounts of LILE, Sr and Pb and do not transport significant amounts of key trace elements such as LREE, U and Th. This indicates that there is a decoupling of water and trace element release in subducted oceanic crust and that aqueous fluids are unable to enrich the mantle wedge significantly. Instead we propose that fluid-present melting in the sediments on top of the slab is required to transfer significant amounts of trace elements from the slab to the mantle wedge. For such a process to be efficient, top slab temperature must be at least 700–750 °C at sub-arc depth. Slab melting is likely to be triggered by fluids that derive from dehydration of mafic and ultramafic rocks in colder (deeper) portions of the slab.     

1kbar、800℃下REE在富磷过铝质熔体/流体相间分配的实验研究

利用"RQV-快速内冷淬火"(或称之为"外加热冷封式")高温高压实验装置,实验研究了1kbar、800℃条件下12个REE+Y在富磷过铝质熔体/含水流体相间的分配,并利用EMP、LA-ICPMS和ICP-MS分析技术分别测定了实验初始物、实验产物玻璃中主要化学组成以及熔体相和流体相中REE含量。实验结果表明,REE元素(La,Nd,Sm,Eu,Gd,Tb,Dy,Ho,Er,Tm,Yb和Lu)在流体/熔体相间的分配系数(Dfluid/melt)在(0.1～19.9)×10-4范围,DfYluid/melt在(0.2～7.8)×10-4范围,指示REE和Y强烈趋向于在熔体中富集。REE在流体/熔体相间的分配系数(Dfluid/melt)与体系中P2O5含量变化呈近抛物线状分布,其最大值对应于残余熔体中w(P2O5)为1.44%处。REE在流体/熔体相间的分配系数(Dfluid/melt)随REE的原子序数增大而逐渐降低,构成右倾的平滑曲线,总体上显示出DLREE>DMREE>DHREE的趋势。Y与Ho在流体/熔体相间分配系数的比值(DY/DHo)约为1(0.91～1.28),不随体系中P2O5变化而变化的特征。上述特征表明熔体-流体作用不会导致Y-Ho及REE间的分异,因此,可推断熔体-流体作用过程不可能是过铝质岩浆体系中产生稀土"四重效应"机制。     

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.     

Role of fluid in the genesis of carbonatites and alkaline rocks: Geochemical evidence

In most alkaline-ultrabasic-carbonatite ring complexes, the distribution of trace elements in the successive derivatives of mantle magmas is usually controlled by the Rayleigh equation of fractional crystallization in accordance with their partition coefficients, whereas, that of late derivatives, nepheline syenites and carbonatites, is usually consistent with trends characteristic of silicate-carbonate liquid immiscibility. In contrast to the carbonatites of ring complexes, carbonatites from deep-seated linear zones have no genetic relation with alkaline-ultrabasic magmatism, and the associated alkaline rocks are represented only by the nepheline syenite eutectic association. The geochemical study of magmatic rocks from the Vishnevye Gory nepheline syenite-carbonatite complex (Urals), which is assigned to the association of deep-seated linear zones, showed that neither differentiation of a parental melt nor liquid immiscibility could produce the observed trace element distribution (Sr, Rb, REE, and Nb) in miaskites and carbonatites. Judging from the available fragmentary experimental data, the distribution patterns can be regarded as possible indicators of element fractionation between alkaline carbonate fluid and alkaline melt. Such trace element distribution is presumably controlled by a fluid-melt interaction; it was also observed in carbonatites and alkaline rocks of some ring complexes, and its scarcity can be explained by the lower density of aqueous fluid released from magma at shallower depths.     

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.     

Gold solubility in oxidized and reduced, water-saturated mafic melt

We have performed experiments to evaluate Au solubility in natural, water-saturated basaltic melts as a function of oxygen fugacity. Experiments were carried out at 1000 °C and 200 MPa, and oxygen fugacity was controlled at the fayalite-magnetite-quartz (FMQ) oxygen fugacity buffer and FMQ + 4. All experiments were saturated with a metal-chloride aqueous solution loaded initially as a 10 wt% NaCl eq. fluid. The stable phase assemblage at FMQ consists of basalt melt, olivine, clinopyroxene, a single-phase aqueous fluid, and metallic Au. The stable phase assemblage at FMQ + 4 consists of basalt melt, clinopyroxene, magnetite-spinel solid solution, a single-phase aqueous fluid, and metallic Au. Silicate glasses (i.e., quenched melt) and their contained crystalline material were analyzed by using both electron probe microanalysis (EPMA) and laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS). Measured Au concentrations in the quenched melt range from 4.8 μg g⁻¹ to 0.64 μg g⁻¹ at FMQ + 4, and 0.54 μg g⁻¹ to 0.1 μg g⁻¹ at FMQ. The measured solubility of Au in olivine and clinopyroxene was consistently below the LA-ICP-MS limit of detection (i.e., 0.1 μg g⁻¹). These melt solubility data place important limitations on the dissolved Au content of water-saturated, Cl- and S-bearing basaltic liquids at geologically relevant fO₂ values. The new data are compared to published, experimentally-determined values for Au solubility in dry and hydrous silicate liquids spanning the compositional range from basalt to rhyolite, and the effects of melt composition, oxygen fugacity, pressure and temperature are discussed.     

Apatite solubility in carbonatitic liquids and trace element partitioning between apatite and carbonatite at high pressure

We have measured apatite solubility in calcic carbonatitic liquids and determined apatite/melt partition coefficients for a series of trace elements, including the rare earth elements (REE), high field strength elements (HFSE), Rb, Sr, U-Th-Pb. Experiments were performed between 4 and 6 GPa, from 1200 to 1380 °C, using the multianvil apparatus. Trace element concentrations were determined by laser ablation ICP-MS and electron microprobe. In addition, a specific protocol was designed to measure carbon concentration in the apatites, using the electron microprobe. Two starting apatite samples were used in order to test for the effect of apatite chemistry on partitioning behavior.Apatite solubility is lower in calcitic melts by a factor 3-5 compared to dolomitic melts (3-5.5 vs. 10-18 wt.% P₂O₅ in melt). We interpret this difference in terms of solubility product in the liquid and propose an empirical model for apatite saturation that takes into account melt calcium content. We conclude that calcitic melts that may form by melting of carbonated eclogites could be saturated with residual apatite, contrary to dolomitic melts formed in carbonated peridotites.Compatibility behavior of the REE depends on apatite silica content: REE are compatible in apatites containing 3.5-5 wt.% SiO₂, with values between 1.5 and 4, whereas REE are incompatible in apatites containing 0.2 wt.% SiO₂. HFSE, U, Th, and Y are compatible in silica-rich apatite, with while . Strontium is always retained in the melt, with of the order of 0.5. Lead appears to be incompatible in apatite, although this finding is weakened by almost complete Pb loss to sample container. High silica concentration favors REE incorporation in apatite by allowing for charged balanced coupled substitution. Sulfur and carbonate may also favor REE incorporation in apatite. Our results allow to reconcile previously published experimental determinations of REE partitioning. We use our experimentally determined partition coefficients to investigate the impact of residual apatite during partial melting of recycled carbonated material (eclogite + sediments) and discuss how the chemical characteristics of the produced liquids can be affected by residual apatite.