The influence of melt structure on trace element partitioning near the peridotite solidus

This experimental study examines the mineral/melt partitioning of Na, Ti, La, Sm, Ho, and Lu among high-Ca clinopyroxene, plagioclase, and silicate melts analogous to varying degrees of peridotite partial melting. Experiments performed at a pressure of 1.5 GPa and temperatures of 1,285 to 1,345 °C produced silicate melts saturated with high-Ca clinopyroxene, plagioclase and/or spinel, and, in one case, orthopyroxene and garnet. Partition coefficients measured in experiments in which clinopyroxene coexists with basaltic melt containing ~18 to 19 wt% Al₂O₃ and up to ~3 wt% Na₂O are consistent with those determined experimentally in a majority of the previous studies, with values of ~0.05 for the light rare earths and of ~0.70 for the heavy rare earths. The magnitudes of clinopyroxene/melt partition coefficients for the rare earth elements correlate with pyroxene composition in these experiments, and relative compatibilities are consistent with the effects of lattice strain energy. Clinopyroxene/melt partition coefficients measured in experiments in which the melt contains ~20 wt% Al₂O₃ and ~4 to 8 wt% Na₂O are unusually large (e.g., values for Lu of up to 1.33±0.05) and are not consistent with the dependence on pyroxene composition found in previous studies. The magnitudes of the partition coefficients measured in these experiments correlate with the degree of polymerization of the melt, rather than with crystal composition, indicating a significant melt structural influence on trace element partitioning. The ratio of non-bridging oxygens to tetrahedrally coordinated cations (NBO/T) in the melt provides a measure of this effect; melt structure has a significant influence on trace element compatibility only for values of NBO/T less than ~0.49. This result suggests that when ascending peridotite intersects the solidus at relatively low pressures (~1.5 GPa or less), the compatibility of trace elements in the residual solid varies significantly during the initial stages of partial melting in response to the changing liquid composition. It is unlikely that this effect is important at higher pressures due to the increased compatibility of SiO₂, Na₂O, and Al₂O₃ in the residual peridotite, and correspondingly larger values of NBO/T for small degree partial melts.Editorial responsibility: T.L. Grove     

Partial melting of secondary pyroxenite at 1 and 1.5 GPa,and its role in upwelling heterogeneous mantle

We performed partial melting experiments at 1 and 1.5 GPa, and 1180–1400 °C, to investigate the melting under mantle conditions of an olivine-websterite (GV10), which represents a natural proxy of secondary (or stage 2) pyroxenite. Its subsolidus mineralogy consists of clinopyroxene, orthopyroxene, olivine and spinel (+garnet at 1.5 GPa). Solidus temperature is located between 1180 and 1200 °C at 1 GPa, and between 1230 and 1250 °C at 1.5 GPa. Orthopyroxene (±garnet), spinel and clinopyroxene are progressively consumed by melting reactions to produce olivine and melt. High coefficient of orthopyroxene in the melting reaction results in relatively high SiO₂ content of low melt fractions. After orthopyroxene exhaustion, melt composition is controlled by the composition of coexisting clinopyroxene. At increasing melt fraction, CaO content of melt increases, whereas Na₂O, Al₂O₃ and TiO₂ behave as incompatible elements. Low Na₂O contents reflect high partition coefficient of Na between clinopyroxene and melt (\(D_{{{\text{Na}}_{ 2} {\text{O}}}}^{{{\text{cpx}}/{\text{liquid}}}}\)). Melting of GV10 produces Quartz- to Hyperstene-normative basaltic melts that differ from peridotitic melts only in terms of lower Na₂O and higher CaO contents. We model the partial melting of mantle sources made of different mixing of secondary pyroxenite and fertile lherzolite in the context of adiabatic oceanic mantle upwelling. At low potential temperatures (T _P < 1310 °C), low-degree melt fractions from secondary pyroxenite react with surrounding peridotite producing orthopyroxene-rich reaction zones (or refertilized peridotite) and refractory clinopyroxene-rich residues. At higher T _P (1310–1430 °C), simultaneous melting of pyroxenite and peridotite produces mixed melts with major element compositions matching those of primitive MORBs. This reinforces the notion that secondary pyroxenite may be potential hidden components in MORB mantle source.     

Phase relations and melting of carbonated peridotite between 10 and 20 GPa: a proxy for alkali- and CO<Subscript>2</Subscript>-rich silicate melts in the deep mantle

We determined the melting phase relations, melt compositions, and melting reactions of carbonated peridotite on two carbonate-bearing peridotite compositions (ACP: alkali-rich peridotite + 5.0 wt % CO₂ and PERC: fertile peridotite + 2.5 wt % CO₂) at 10–20 GPa and 1,500–2,100 °C and constrain isopleths of the CO₂ contents in the silicate melts in the deep mantle. At 10–20 GPa, near-solidus (ACP: 1,400–1,630 °C) carbonatitic melts with < 10 wt % SiO₂ and > 40 wt % CO₂ gradually change to carbonated silicate melts with > 25 wt % SiO₂ and < 25 wt % CO₂ between 1,480 and 1,670 °C in the presence of residual majorite garnet, olivine/wadsleyite, and clinoenstatite/clinopyroxene. With increasing degrees of melting, the melt composition changes to an alkali- and CO₂-rich silicate melt (Mg# = 83.7–91.6; ~ 26–36 wt % MgO; ~ 24–43 wt % SiO₂; ~ 4–13 wt % CaO; ~ 0.6–3.1 wt % Na₂O; and ~ 0.5–3.2 wt % K₂O; ~ 6.4–38.4 wt % CO₂). The temperature of the first appearance of CO₂-rich silicate melt at 10–20 GPa is ~ 440–470 °C lower than the solidus of volatile-free peridotite. Garnet + wadsleyite + clinoenstatite + carbonatitic melt controls initial carbonated silicate melting at a pressure < 15 GPa, whereas garnet + wadsleyite/ringwoodite + carbonatitic melt dominates at pressure > 15 GPa. Similar to hydrous peridotite, majorite garnet is a liquidus phase in carbonated peridotites (ACP and PERC) at 10–20 GPa. The liquidus is likely to be at ~ 2,050 °C or higher at pressures of the present study, which gives a melting interval of more than 670 °C in carbonated peridotite systems. Alkali-rich carbonated silicate melts may thus be produced through partial melting of carbonated peridotite to 20 GPa at near mantle adiabat or even at plume temperature. These alkali- and CO₂-rich silicate melts can percolate upward and may react with volatile-rich materials accumulate at the top of transition zone near 410-km depth. If these refertilized domains migrate upward and convect out of the zone of metal saturation, CO₂ and H₂O flux melting can take place and kimberlite parental magmas can be generated. These mechanisms might be important for mantle dynamics and are potentially effective metasomatic processes in the deep mantle.     

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

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

The influence of water on melting of mantle peridotite

This experimental study examines the effects of variable concentrations of dissolved H₂O on the compositions of silicate melts and their coexisting mineral assemblage of olivine + orthopyroxene ± clinopyroxene ± spinel ± garnet. Experiments were performed at pressures of 1.2 to 2.0 GPa and temperatures of 1100 to 1345 °C, with up to ∼12 wt% H₂O dissolved in the liquid. The effects of increasing the concentration of dissolved H₂O on the major element compositions of melts in equilibrium with a spinel lherzolite mineral assemblage are to decrease the concentrations of SiO₂, FeO, MgO, and CaO. The concentration of Al₂O₃ is unaffected. The lower SiO₂ contents of the hydrous melts result from an increase in the activity coefficient for SiO₂ with increasing dissolved H₂O. The lower concentrations of FeO and MgO result from the lower temperatures at which H₂O-bearing melts coexist with mantle minerals as compared to anhydrous melts. These compositional changes produce an elevated SiO₂/(MgO + FeO) ratio in hydrous peridotite partial melts, making them relatively SiO₂ rich when compared to anhydrous melts on a volatile-free basis. Hydrous peridotite melting reactions are affected primarily by the lowered mantle solidus. Temperature-induced compositional variations in coexisting pyroxenes lower the proportion of clinopyroxene entering the melt relative to orthopyroxene. Isobaric batch melting calculations indicate that fluid-undersaturated peridotite melting is characterized by significantly lower melt productivity than anhydrous peridotite melting, and that the peridotite melting process in subduction zones is strongly influenced by the composition of the H₂O-rich component introduced into the mantle wedge from the subducted slab. Received: 7 April 1997 / Accepted: 9 January 1998     

Experimental comparison of trace element partitioning between clinopyroxene and melt in carbonate and silicate systems, and implications for mantle metasomatism

Experiments in the systems diopside-albite (Di-Ab) and diopside-albite-dolomite (Di-Ab-Dmt), doped with a wide range of trace elements, have been used to characterise the difference between clinopyroxene-silicate melt and clinopyroxene-carbonate melt partitioning. Experiments in Di-Ab-Dmt yielded clinopyroxene and olivine in equilibrium with CO₂-saturated dolomitic carbonate melt at 3 GPa, 1375 °C. The experiments in Di-Ab were designed to bracket those conditions (3 GPa, 1640 °C and 0.8 GPa, 1375 °C), and so minimise the contribution of differential temperature and pressure to partitioning. Partition coefficients, determined by SIMS analysis of run products, differ markedly for some elements between Di-Ab and Di-Ab-Dmt systems. Notably, in the carbonate system clinopyroxene-melt partition coefficients for Si, Al, Ga, heavy REE, Ti and Zr are higher by factors of 5 to 200 than in the silicate system. Conversely, partition coefficients for Nb, light REE, alkali metals and alkaline earths show much less fractionation (<3). The observed differences compare quantitatively with experimental data on partitioning between immiscible carbonate and silicate melts, indicating that changes in melt chemistry provide the dominant control on variation in partition coefficients in this case. The importance of melt chemistry in controlling several aspects of element partitioning is discussed in light of the energetics of the partitioning process. The compositions of clinopyroxene and carbonate melt in our experiments closely match those of near-solidus melts and crystals in CMAS-CO₂ at 3 GPa, suggesting that our partition coefficients have direct relevance to melting of carbonated mantle lherzolite. Melts so produced will be characterised by elevated incompatible trace element concentrations, due to the low degrees of melting involved, but marked depletions of Ti and Zr, and fractionated REE patterns. These are common features of natural carbonatites. The different behaviour of trace elements in carbonate and silicate systems will lead to contrasted styles of trace element metasomatism in the mantle. Received: 15 July 1999 / Accepted: 18 February 2000     

Peridotite–eclogite–carbonatite systems at 7.0–8.5 GPa: concentration barrier of diamond nucleation and syngenesis of its silicate and carbonate inclusions

Diamond crystallization in multicomponent melts of variable composition is studied. The melt carbonates are K₂CO₃, CaCO₃?MgCO₃, and K-Na-Ca-Mg-Fe-carbonatites, and the melt silicates are model peridotite (60 wt.% olivine, 16 wt.% orthopyroxene, 12 wt.% clinopyroxene, and 12 wt.% garnet) and eclogite (50 wt.% garnet and 50 wt.% clinopyroxene). In the experiments carried out under the PT-conditions of diamond stability, the carbonate-silicate melts behave like completely miscible liquid phases. The concentration barriers of diamond nucleation (CBDN) in the melts with variable proportions of silicates and carbonates have been determined at 8.5 GPa. In the system peridotite–K₂CO₃–CaCO₃?MgCO₃–carbonatite they correspond to 30, 25, and 30 wt.% silicates, respectively, and in the analogous eclogite–carbonate system, 45, 30, and 35 wt.%. In the silicate-carbonate melts with higher silicate contents seed diamond growth occurs, which is accompanied by the crystallization of thermodynamically unstable graphite phase. In the experiments with melts compositionally corresponding to the CBDN at 7.0 GPa and 1200–1700 °C, a full set of silicate minerals of peridotite (olivine, orthopyroxene, clinopyroxene, garnet) and eclogite (garnet, clinopyroxene) parageneses was obtained. The minerals occur as syngenetic inclusions in natural diamonds; moreover, the garnets contain an impurity of Na, and the pyroxenes, K. The experimental data indicate that peridotite-carbonate and eclogite-carbonate melts are highly effective for the formation of diamond (or unstable graphite) together with syngenetic minerals and melts, which agrees with the carbonate-silicate (carbonatite) model for the mantle diamond formation.     

Melting of carbonated pelites at 8–13 GPa: generating K-rich carbonatites for mantle metasomatism

The melting behaviour of three carbonated pelites containing 0–1 wt% water was studied at 8 and 13 GPa, 900–1,850°C to define conditions of melting, melt compositions and melting reactions. At 8 GPa, the fluid-absent and dry carbonated pelite solidi locate at 950 and 1,075°C, respectively; >100°C lower than in carbonated basalts and 150–300°C lower than the mantle adiabat. From 8 to 13 GPa, the fluid-present and dry solidi temperatures then increase to 1,150 and 1,325°C for the 1.1 wt% H₂O and the dry composition, respectively. The melting behaviour in the 1.1 wt% H₂O composition changes from fluid-absent at 8 GPa to fluid-present at 13 GPa with the pressure breakdown of phengite and the absence of other hydrous minerals. Melting reactions are controlled by carbonates, and the potassium and hydrous phases present in the subsolidus. The first melts, which composition has been determined by reverse sandwich experiments, are potassium-rich Ca–Fe–Mg-carbonatites, with extreme K₂O/Na₂O wt ratios of up to 42 at 8 GPa. Na is compatible in clinopyroxene with D_\textNa^{\textcpx/\textcarbonatite} = 10-18 D_{\text{Na}}^{{{\text{cpx}}/{\text{carbonatite}}}} = 10{-}18 at the solidus at 8 GPa. The melt K₂O/Na₂O slightly decreases with increasing temperature and degree of melting but strongly decreases from 8 to 13 GPa when K-hollandite extends its stability field to 200°C above the solidus. The compositional array of the sediment-derived carbonatites is congruent with alkali- and CO₂-rich melt or fluid inclusions found in diamonds. The fluid-absent melting of carbonated pelites at 8 GPa contrasts that at ≤5 GPa where silicate melts form at lower temperatures than carbonatites. Comparison of our melting temperatures with typical subduction and mantle geotherms shows that melting of carbonated pelites to 400-km depth is only feasible for extremely hot subduction. Nevertheless, melting may occur when subduction slows down or stops and thermal relaxation sets in. Our experiments show that CO₂-metasomatism originating from subducted crust is intimately linked with K-metasomatism at depth of >200 km. As long as the mantle remains adiabatic, low-viscosity carbonatites will rise into the mantle and percolate upwards. In cold subcontinental lithospheric mantle keels, the potassic Ca–Fe–Mg-carbonatites may freeze when reacting with the surrounding mantle leading to potassium-, carbonate/diamond- and incompatible element enriched metasomatized zones, which are most likely at the origin of ultrapotassic magmas such as group II kimberlites.     

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.     

Stability of phlogopite in ultrapotassic kimberlite-like systems at 5.5–7.5 GPa

Hydrous K-rich kimberlite-like systems are studied experimentally at 5.5–7.5 GPa and 1200–1450?°C in terms of phase relations and conditions for formation and stability of phlogopite. The starting samples are phlogopite–carbonatite–phlogopite sandwiches and harzburgite–carbonatite mixtures consisting of Ol?+?Grt?+?Cpx?+?L (±Opx), according to the previous experimental results obtained at the same P–T parameters but in water-free systems. Carbonatite is represented by a K- and Ca-rich composition that may form at the top of a slab. In the presence of carbonatitic melt, phlogopite can partly melt in a peritectic reaction at 5.5 GPa and 1200–1350?°C, as well as at 6.3–7.0 GPa and 1200?°C: 2Phl?+?CaCO₃ (L)?Cpx?+?Ol?+?Grt?+?K₂CO₃ (L)?+?2H₂O (L). Synthesis of phlogopite at 5.5 GPa and 1200–1350?°C, with an initial mixture of H₂O-bearing harzburgite and carbonatite, demonstrates experimentally that equilibrium in this reaction can be shifted from right to left. Therefore, phlogopite can equilibrate with ultrapotassic carbonate–silicate melts in a?≥?150?°C region between 1200 and 1350?°C at 5.5 GPa. On the other hand, it can exist but cannot nucleate spontaneously and crystallize in the presence of such melts in quite a large pressure range in experiments at 6.3–7.0 GPa and 1200?°C. Thus, phlogopite can result from metasomatism of peridotite at the base of continental lithospheric mantle (CLM) by ultrapotassic carbonatite agents at depths shallower than 180–195 km, which creates a mechanism of water retaining in CLM. Kimberlite formation can begin at 5.5 GPa and 1350?°C in a phlogopite-bearing peridotite source generating a hydrous carbonate–silicate melt with 10–15 wt% SiO₂, Ca# from 45 to 60, and high K enrichment. Upon further heating to 1450?°C due to the effect of a mantle plume at the CLM base, phlogopite disappears and a kimberlite-like melt forms with SiO₂ to 20 wt% and Ca#?=?35–40.     

Phase relations and chemistry of Ti-rich K-richterite-bearing mantle assemblages: an experimental study to 8.0 GPa in a Ti-KNCMASH system

Experiments have been conducted in a peralkaline Ti-KNCMASH system representative of MARID-type bulk compositions to delimit the stability field of K-richterite in a Ti-rich hydrous mantle assemblage, to assess the compositional variation of amphibole and coexisting phases as a function of P and T, and to characterise the composition of partial melts derived from the hydrous assemblage. K-richterite is stable in experiments from 0.5 to 8.0 GPa coexisting with phlogopite, clinopyroxene and a Ti-phase (titanite, rutile or rutile + perovskite). At 8.0 GPa, garnet appears as an additional phase. The upper T stability limit of K-richterite is 1200–1250 °C at 4.0 GPa and 1300–1400 °C at 8.0 GPa. In the presence of phlogopite, K-richterite shows a systematic increase in K with increasing P to 1.03 pfu (per formula unit) at 8.0 GPa/1100 °C. In the absence of phlogopite, K-richterite attains a maximum of 1.14 K pfu at 8.0 GPa/1200 °C. Titanium in both amphibole and mica decreases continuously towards high P with a nearly constant partitioning while Ti in clinopyroxene remains more or less constant. In all experiments below 6.0 GPa ΣSi + Al in K-richterite is less than 8.0 when normalised to 23 oxygens+stoichiometric OH. Rutiles in the Ti-KNCMASH system are characterised by minor Al and Mg contents that show a systematic variation in concentration with P(T) and the coexisting assemblage. Partial melts produced in the Ti-KNCMASH system are extremely peralkaline [(K₂O+Na₂O)/Al₂O₃ = 1.7–3.7], Si-poor (40–45 wt% SiO₂), and Ti-rich (5.6–9.2 wt% TiO₂) and are very similar to certain Ti-rich lamproite glasses. At 4.0 GPa, the solidus is thought to coincide with the K-richterite-out reaction, the first melt is saturated in a phlogopite-rutile-lherzolite assemblage. Both phlogopite and rutile disappear ca. 150 °C above the solidus. At 8.0 GPa, the solidus must be located at T≤1400 ^°C. At this temperature, a melt is in equilibrium with a garnet- rutile-lherzolite assemblage. As opposed to 4.0 GPa, phlogopite does not buffer the melt composition at 8.0 GPa. The experimental results suggest that partial melting of MARID-type assemblages at pressures ≥4.0 GPa can generate Si-poor and partly ultrapotassic melts similar in composition to that of olivine lamproites. Received: 23 December 1996 / Accepted: 20 March 1997     

Experiments on melt–rock reaction in the shallow mantle wedge

This experimental study simulates the interaction of hotter, deeper hydrous mantle melts with shallower, cooler depleted mantle, a process that is expected to occur in the upper part of the mantle wedge. Hydrous reaction experiments (~6 wt% H₂O in the melt) were conducted on three different ratios of a 1.6 GPa mantle melt and an overlying 1.2 GPa harzburgite from 1060 to 1260 °C. Reaction coefficients were calculated for each experiment to determine the effect of temperature and starting bulk composition on final melt compositions and crystallizing assemblages. The experiments used to construct the melt–wall rock model closely approached equilibrium and experienced <5% Fe loss or gain. Experiments that experienced higher extents of Fe loss were used to critically evaluate the practice of “correcting” for Fe loss by adding iron. At low ratios of melt/mantle (20:80 and 5:95), the crystallizing assemblages are dunites, harzburgites, and lherzolites (as a function of temperature). When the ratio of deeper melt to overlying mantle is 70:30, the crystallizing assemblage is a wehrlite. This shows that wehrlites, which are observed in ophiolites and mantle xenoliths, can be formed by large amounts of deeper melt fluxing though the mantle wedge during ascent. In all cases, orthopyroxene dissolves in the melt, and olivine crystallizes along with pyroxenes and spinel. The amount of reaction between deeper melts and overlying mantle, simulated here by the three starting compositions, imposes a strong influence on final melt compositions, particularly in terms of depletion. At the lowest melt/mantle ratios, the resulting melt is an extremely depleted Al-poor, high-Si andesite. As the fraction of melt to mantle increases, final melts resemble primitive basaltic andesites found in arcs globally. An important element ratio in mantle lherzolite composition, the Ca/Al ratio, can be significantly elevated through shallow mantle melt–wall rock reaction. Wall rock temperature is a key variable; over a span of <80 °C, reaction with deeper melt creates the entire range of mantle lithologies from a depleted dunite to a harzburgite to a refertilized lherzolite. Together, the experimental phase equilibria, melt compositions, and reaction coefficients provide a framework for understanding how melt–wall rock reaction occurs in the natural system during melt ascent in the mantle wedge.     

Fluid and melt inclusions in the Mesozoic Fangcheng basalt from North China Craton: implications for magma evolution and fluid/melt-peridotite reaction

Melt inclusions and fluid inclusions in the Fangcheng basalt were investigated to understand the magma evolution and fluid/melt-peridotite interaction. Primary silicate melt inclusions were trapped in clinopyroxene and orthopyroxene phenocrysts in the Fangcheng basalt. Three types of melt inclusions (silicate, carbonate, and sulfide) coexisting with fluid inclusions occur in clinopyroxene xenocrysts and clinopyroxene in clinopyroxenite xenoliths. In situ laser-ablation ICP-MS analyses of major and trace element compositions on individual melt inclusions suggest that the silicate melt inclusions in clinopyroxene and orthopyroxene phenocrysts were trapped from the same basaltic magma. The decoupling of major and trace elements in the melt inclusions indicates that the magma evolution was controlled by melt crystallization and contamination from entrapped ultramafic xenoliths. Trace element patterns of melt inclusions are similar to those of the average crust of North China Craton and Yangtze Craton, suggesting a considerable crustal contribution to the magma source. Calculated parental melt of the Fangcheng basalt has features of low MgO (5.96 wt%), high Al₂O₃ (16.81 wt%), Sr (1,670 ppm), Y (>35 ppm), and high Sr/Y (>40), implying that subducted crustal material was involved in the genesis of the Fangcheng basalt. The coexisting fluid and melt inclusions in clinopyroxene xenocrysts and in clinopyroxene of xenoliths record a rare melt-peridotite reaction, that is olivine + carbonatitic melt₁ (rich in Ca) = clinopyroxene + melt₂ ± CO₂. The produced melt₂ is enriched in LREE and CO₂ and may fertilize the mantle significantly, which we consider to be the cause for the rapid replacement of lithospheric mantle during the Mesozoic in the region.     

Source and H2O content of high-MgO magmas in island arc settings: an experimental study of a primitive calc-alkaline basalt from St. Vincent, lesser antilles arc

Liquidus phase relationships have been determined for a high-MgO basalt (STV301: MgO=12.5 wt%, Ni=250 ppm, Cr=728 ppm) from Black Point, St Vincent (Lesser Antilles arc). Piston-cylinder experiments were conducted between 7.5 and 20 kbar under both hydrous and oxidizing conditions. AuPd capsules were used as containers. Compositions of supraliquidus glasses and mass-balance calculations show that Fe loss is < 10% in the majority of experiments. Two series of water concentrations in melt were investigated: (i) 1.5 wt% and (ii) 4.5 wt% H₂O, as determined by SIMS analyses on quenched glasses and with the by difference technique. The Fe³⁺/Fe²⁺ partitioning between Cr-Al spinel and melt and olivine-spinel equilibria show that oxidizing fO₂ were imposed (NNO + 1.5 for the 1.5 wt% H₂O series, NNO + 2.3 for the 4.5 wt% H₂O series). For both series of water concentrations, the liquid is multiply-saturated with a spinel lherzolite phase assemblage on its liquidus, at 1235°C, 11.5 kbar (1.5 wt% H₂O) and 1185°C, 16 kbar (4.5 wt% H₂O). Liquidus phases are homogeneous and comparable to typical mantle compositions. Mineral-melt partition coefficients are generally identical to values under anhydrous conditions. The modal proportion cpx/opx on the liquidus decreases from the 1.5 wt% to the 4.5 wt% H₂O series. The experimental data are consistent with STV301 being a product of partial melting of lherzolitic mantle. Conditions of multiple saturation progressively evolve toward lower temperatures and higher pressures with increasing melt H₂O concentration. Phase equilibria constraints, i.e., the necessity of preserving the mantle signature seen in high-MgO and picritic arc basalts, and glass inclusion data suggest that STV301 was extracted relatively dry (∼ 2 wt% H₂O) from its mantle source. However, not all primary arc basalts are extracted under similarly dry conditions because more hydrous melts will crystallize during ascent and will not be present unmodified at the surface. From degrees of melting calculated from experiments on KLB-1, extraction of a 12.5 wt% MgO melt with ∼ 2 wt% H₂O would require a H₂O concentration of 0.3 wt% in the sub-arc mantle. For mantle sources fluxed with a slab-derived hydrous component, extracted melts may contain up to ∼ 5.5 wt% H₂O.     

Distribution of REE between clinopyroxene and basaltic melt along a mantle adiabat: effects of major element composition, water, and temperature

The distribution of rare earth elements (REE) between clinopyroxene (cpx) and basaltic melt is important in deciphering the processes of mantle melting. REE and Y partition coefficients from a given cpx-melt partitioning experiment can be quantitatively described by the lattice strain model. We analyzed published REE and Y partitioning data between cpx and basaltic melts using the nonlinear regression method and parameterized key partitioning parameters in the lattice strain model (D ₀, r ₀ and E) as functions of pressure, temperature, and compositions of cpx and melt. D ₀ is found to positively correlate with Al in tetrahedral site (Al ^T ) and Mg in the M2 site (Mg^M2) of cpx and negatively correlate with temperature and water content in the melt. r ₀ is negatively correlated with Al in M1 site (Al^M1) and Mg^M2 in cpx. And E is positively correlated with r ₀. During adiabatic melting of spinel lherzolite, temperature, Al ^T , and Mg^M2 in cpx all decrease systematically as a function of pressure or degree of melting. The competing effects between temperature and cpx composition result in very small variations in REE partition coefficients along a mantle adiabat. A higher potential temperature (1,400°C) gives rise to REE partition coefficients slightly lower than those at a lower potential temperature (1,300°C) because the temperature effect overwhelms the compositional effect. A set of constant REE partition coefficients therefore may be used to accurately model REE fractionation during partial melting of spinel lherzolite along a mantle adiabat. As cpx has low Al and Mg abundances at high temperature during melting in the garnet stability field, REE are more incompatible in cpx. Heavy REE depletion in the melt may imply deep melting of a hydrous garnet lherzolite. Water-dependent cpx partition coefficients need to be considered for modeling low-degree hydrous melting.     

Mineral equilibria of diamond-forming carbonate-silicate systems

Based on experimental and mineralogical data, the model of mantle carbonate-silicate (carbonatite) melts as dominating parental media for natural diamonds was substantiated. It was demonstrated that the compositions of silicate constituents of parental melts were variable and saturated with respect to mantle rocks, namely pyrope peridotite, garnet pyroxenite, and eclogite. Based on concentration contributions and role in diamond genesis, major (carbonate and silicate) and minor (admixture) components were distinguished. The latter components may be both soluble (oxides, phosphates, chlorides, carbon dioxide, and water) and insoluble (sulfides, metals, and carbides) in silicate-carbonate melts. This paper presents the results of a study of diamond crystallization in multicomponent melts of variable composition with carbonate components (K₂CO₃, CaCO₃ · MgCO₃, and K-Na-Ca-Mg-Fe carbonatite) and silicate components represented by model peridotite (60 wt % olivine, 16 wt % orthopyroxene, 12 wt % clinopyroxene, and 12 wt % garnet) and eclogite (50 wt % garnet and 50 wt % clinopyroxene). Carbonate-silicate melts behave like completely miscible liquid phases in experiments performed under the P-T conditions of diamond stability. The concentration barriers of diamond nucleation (CBDN) in melts with variable proportions of silicates and carbonates were determined at 8.5 GPa. In the peridotite system with K₂CO₃, CaCO₃ · MgCO₃, and carbonatite, they correspond to 30, 25, and 30 wt % silicates, respectively, and in the eclogite system, the CBDN is shifted to 45, 30, and 35 wt % silicates. In the silicate-carbonate melts with higher silicate contents, diamond grows on seeds, which is accompanied by the crystallization of thermodynamically unstable graphite. At P = 7.0 GPa and T = 1200−1800°C, we studied and constructed phase diagrams for the multicomponent peridotite-carbonate and eclogite-carbonate systems as a physicochemical basis for revealing the syngenetic relationships between diamond and its silicate (olivine, ortho- and clinopyroxene, and garnet) and carbonate (aragonite and magnesite) inclusions depending on the physicochemical conditions of growth media. The results obtained allowed us to reconstruct the evolution of diamond-forming systems. The experiments revealed similarity between the compositions of synthetic silicate minerals and inclusions in natural diamonds (high concentrations of Na in garnets and K in clinopyroxenes). It was experimentally demonstrated that the formation of Na-bearing majoritic garnets is controlled by the P-T parameters and melt alkalinity. Diamonds with inclusions of such garnets can be formed in alkalic carbonate-silicate (aluminosilicate) melts. A mechanism was suggested for sodic end-member dissolution in majoritic garnets, and garnet with the composition Na₂MgSi₅O₁₂ and tetragonal symmetry was synthesized for the first time.     

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.     

Dissolution of orthopyroxene in basanitic magma between 0.4 and 2 GPa: further implications for the origin of Si-rich alkaline glass inclusions in mantle xenoliths

A large body of recent work has linked the origin of Si-Al-rich alkaline glass inclusions to metasomatic processes in the upper mantle. This study examines one possible origin for these glass inclusions, i.e., the dissolution of orthopyroxene in Si-poor alkaline (basanitic) melt. Equilibrium dissolution experiments between 0.4 and 2 GPa show that secondary glass compositions are only slightly Si enriched and are alkali poor relative to natural glass inclusions. However, disequilibrium experiments designed to examine dissolution of orthopyroxene by a basanitic melt under anhydrous, hydrous and CO₂-bearing conditions show complex reaction zones consisting of olivine, ± clinopyroxene and Si-rich alkaline glass similar in composition to that seen in mantle xenoliths. Dissolution rates are rapid and dependent on volatile content. Experiments using an anhydrous solvent show time dependent dissolution rates that are related to variable diffusion rates caused by the saturation of clinopyroxene in experiments longer than 10 minutes. The reaction zone glass shows a close compositional correspondence with natural Si-rich alkaline glass in mantle-derived xenoliths. The most Si-and alkali-rich melts are restricted to pressures of 1 GPa and below under anhydrous and CO₂-bearing conditions. At 2 GPa glass in hydrous experiments is still Si-␣and alkali-rich whereas glass in the anhydrous and CO₂-bearing experiments is only slightly enriched in SiO₂ and alkalis compared with the original solvent. In the low pressure region, anhydrous and hydrous solvent melts yield glass of similar composition whereas the glass from CO₂-bearing experiments is less SiO₂ rich. The mechanism of dissolution of orthopyroxene is complex involving rapid incongruent breakdown of the orthopyroxene, combined with olivine saturation in the reaction zone forming up to 60% olivine. Inward diffusion of CaO causes clinopyroxene saturation and uphill diffusion of Na and K give the glasses their strongly alkaline characteristics. Addition of Na and K also causes minor SiO₂ enrichment of the reaction glass by increasing the phase volume of olivine. Olivine and clinopyroxene are transiently stable phases within the reaction zone. Clinopyroxene is precipitated from the reaction zone melt near the orthopyroxene crystal and redissolved in the outer part of the reaction zone. Olivine defines the thickness of the reaction zone and is progressively dissolved in the solvent as the orthopyroxene continues to dissolve. Although there are compelling reasons for supporting the hypothesis that Si-rich alkaline melts are produced in the mantle by orthopyroxene – melt reaction in the mantle, there are several complications particularly regarding quenching in of disequilibrium reaction zone compositions and the mobility of highly polymerized melts in the upper mantle. It is considered likely that formation of veins and pools of Si-rich alkaline glass by orthopyroxene – melt reaction is a common process during the ascent of xenoliths. However, reaction in situ within the mantle will lead to equilibration and therefore secondary melts will be only moderately siliceous and alkali poor. Received: 24 August 1998 / Accepted: 2 December 1998     

Cu diffusivity in granitic melts with application to the formation of porphyry Cu deposits

We report new experimental data of Cu diffusivity in granite porphyry melts with 0.01 and 3.9 wt% H₂O at 0.15–1.0 GPa and 973–1523 K. A diffusion couple method was used for the nominally anhydrous granitic melt, whereas a Cu diffusion-in method using Pt₉₅Cu₅ as the source of Cu was applied to the hydrous granitic melt. The diffusion couple experiments also generate Cu diffusion-out profiles due to Cu loss to Pt capsule walls. Cu diffusivities were extracted from error function fits of the Cu concentration profiles measured by LA-ICP-MS. At 1 GPa, we obtain \({D_{{\text{Cu, dry, 1 GPa}}}}=\exp \left[ {( - {\text{13.89}} \pm {\text{0.42}}) - \frac{{{\text{12878}} \pm {\text{540}}}}{T}} \right],\) and \({D_{{\text{Cu, 3}}{\text{.9 wt\% }}{{\text{H}}_{\text{2}}}{\text{O}},{\text{ 1 GPa}}}}=\exp \left[ {( - 16.31 \pm 1.30) - \frac{{{\text{8148}} \pm {\text{1670}}}}{T}} \right],\) where D is Cu diffusivity in m²/s and T is temperature in K. The above expressions are in good agreement with a recent study on Cu diffusion in rhyolitic melt using the approach of Cu₂S dissolution. The observed pressure effect over 0.15–1.0 GPa can be described by an activation volume of 5.9 cm³/mol for Cu diffusion. Comparison of Cu diffusivity to alkali diffusivity and its variation with melt composition implies fourfold-coordinated Cu⁺ in silicate melts. Our experimental results indicate that in the formation of porphyry Cu deposits, the diffusive transport of magmatic Cu to sulfide liquids or fluid bubbles is highly efficient. The obtained Cu diffusivity data can also be used to assess whether equilibrium Cu partitioning can be reached within certain experimental durations.