Experimental study of liquid immiscibility in the fluid-magmatic silicate systems containing Ti,Nb, Sr,REE, and Zr

Aluminosilicate alkaline systems containing Ti, REE (La, Ce), Y, Sr, and Nb were studied experimentally at T = 1200 and 1250°C and P = 2 kbar in the presence of aqueous and alkaline fluid. In the fluid-absent experiments with systems containing the same elements, loparite crystals in a silicate matrix were obtained. These systems in the presence of aqueous and alkaline fluid pressure conditions demonstrate an immiscible splitting into two liquids: (1) aluminosilicate matrix and (2) droplets enriched in Ti, REE (La, Ce), Y, Sr, and Nb, which contain silicate admixture and are compositionally close to loparites. According to approximate estimates, the partition coefficients between the melts of the droplets and aluminosilicate matrix (K = C _dr/C _sil) are more than 5 for TiO₂, < 0.35 for SiO₂, 10–20 for Nb₂O₅, > 15 for REE, and from 2.3 to 7.6 for SrO. Obtained unmixing may be of great significance for explaining the genesis of REE-Nb (loparite) deposits. In addition, the experiments demonstrate that Zr may be accumulated together with Ti and REE due to liquid immiscibility of this type. The partition coefficient of ZrO₂ between the melts of droplets and aluminosilicate matrix (K = C _dr/C _sil) according to approximate estimates varies from ∼ 3.5 to 9.     

Liquid Immiscibility and the Evolution of Basaltic Magma

This experimental study examines relationships between alternativeevolution paths of basaltic liquids (the so-called Bowen andFenner trends), and silicate liquid immiscibility. Syntheticanalogues of natural immiscible systems exhibited in volcanicglasses and melt inclusions were used as starting mixtures.Conventional quench experiments in 1 atm gas mixing furnacesproved unable to reproduce unmixing of ferrobasaltic melts,yielding instead either turbid, opalescent glasses, or crystallizationof tridymite and pyroxenes. In contrast, experiments involvingin situ high-temperature centrifugation at 1000g (g = 9·8m/s²) did yield macroscopic unmixing and phase separation. Centrifugationfor 3–4 h was insufficient to complete phase segregation,and resulted in sub-micron immiscible emulsions in quenchedglasses. For a model liquid composition of the Middle Zone ofthe Skaergaard intrusion at super-liquidus temperatures of 1110–1120°C,centrifugation produced a thin, silicic layer (64·5 wt%SiO₂ and 7·4 wt% FeO) at the top of the main Fe-richglass (46 wt% SiO₂ and 21 wt% FeO). The divergent compositionsat the top and bottom were shown in a series of static runsto crystallize very similar crystal assemblages of plagioclase,pyroxene, olivine, and Fe–Ti oxides. We infer from theseresults that unmixing of complex aluminosilicate liquids maybe seriously kinetically hampered (presumably by a nucleationbarrier), and thus conventional static experiments may not correctlyreproduce it. In the light of our centrifuge experiments, immiscibilityin the Skaergaard intrusion could have started already at thetransition from the Lower to the Middle Zone. Thus, magma unmixingmight be an important factor in the development of the Fe-enrichmenttrend documented in the cumulates of the Skaergaard LayeredSeries. KEY WORDS: liquid immiscibility; Skaergaard; layered intrusions; experimental petrology     

Thermochemistry of sulfide liquids III: Ni-bearing liquids at 1 bar

Experiments were performed in the system O–S–Fe–Ni designed to extend our understanding of the chemistry of sulfide liquids. Results indicate that adding nickel to Fe-rich sulfide liquids in equilibrium with silicate liquids extends their stability field to much higher oxygen fugacities and lower sulfur fugacities. Increasing Ni/Fe at a given temperature and sulfur and oxygen fugacity is accompanied by a significant decrease in the oxygen content of the sulfide liquid. Results of these experiments are combined with data from the literature to calibrate an associated regular solution model for O–S–Fe–Ni liquids. This model represents a complete refit of the associated regular solution model of Kress (Contrib Mineral Petrol 139:316–325, 2000). The resulting model is combined with the olivine solution model of Hirschmann (Am Mineral 76:1232–1248, 1991) to explore the effect of variations in oxygen and sulfur fugacities on the distribution of Fe and Ni between olivine and sulfide liquid. Predicted olivine–sulfide distribution trends parallel those observed by Gaetani and Grove (Geochim Cosmochim Acta 61:1829–1846, 1997), Gaetani and Grove (Earth Planet Sci Lett 169:147–163, 1999), Brenan and Caciagli (Geochim Cosmochim Acta 64:307–320, 2000) and Brenan (Geochim Cosmochim Acta 67:2663–2681, 2003), but are systematically offset toward lower predicted Ni in the sulfide. Nevertheless our results are consistent with the assertion that low K _D ^os values in magmatic ore deposits such as the J-M Reef reflect high iron contents in the sulfides combined with relatively high oxygen fugacities.

---

Thermochemistry of sulfide liquids. I. the system O-S-Fe at 1 bar

Recrystallized globules representing former immiscible sulfide liquids are found in a variety of igneous environments. Relatively little is known about the physical properties and thermochemistry of sulfide liquids, despite their importance in igneous systems. This study presents results of a series of experiments designed to calibrate a thermodynamic model for sulfide liquids in the system O-S-Fe at one atmosphere pressure. Sulfide liquids were equilibrated under controlled oxygen and sulfur fugacities at temperatures between 1100 and 1350 ° C in equilibrium with a silica mineral and a silicate melt. Experiments were quenched in a high-speed double-roller “splat” quencher in order to assure that measured compositions were as close to equilibrium liquid values as possible. Sulfide liquids are not stable in equilibrium with a silica-saturated silicate melt at log₁₀(f _O2) > FMQ-1 at 1250 °C and log₁₀(f _S2)=−3. Iron content of the sulfide changes little with variations in oxygen and sulfur fugacity at a given temperature. Consequently, oxygen and sulfur contents are inversely correlated in these liquids. Sulfur is present entirely as sulfide. Iron appears to be present in both its ferric and ferrous states. Data from this study were combined with data compiled from the literature to calibrate an asymmetric regular solution thermodynamic mixing model for O-S-Fe liquids. This model reproduces miscibility gaps and data from this study quite well, but exhibits minor but systematic errors at the O-Fe binary. The observed inverse correlation between sulfur and oxygen is reflected in the predicted free-energy surface by a sharp energy valley running along a line of constant Fe content. Received: 10 April 1996 / Accepted: 15 November 1996     

The influence of crystal field stabilization energy on Fe2+ partitioning in paragenetic minerals

In order to explore possible quantitative relations between crystal field stabilization energy, CFSE, and partitioning behaviour of the 3d ⁶-configured Fe²⁺ ion, a suite of 29 paragenetic rock-forming minerals from 12 high-grade metamorphic rock samples of the Ukrainian shield, including the parageneses garnet/orthopyroxene/clinopyroxene (2x), orthopyroxene/clinopyroxene, garnet/clinopyroxene, garnet/orthopyroxene/biotite, garnet/biotite, garnet/cordierite, garnet/cordierite/biotite, garnet/orthopyroxene/clinopyroxene/Ca-amphibole, Ca-amphibole/biotite (retrograde), was studied by electron microprobe analysis to obtain the respective K _D Fe^{2+ (Ph1/Ph2)} values and by polarized single crystal electronic absorption spectroscopy to evaluate the respective CFSE_Fe2+ values. Other than in the case of Cr³⁺, a clear quantitative relation between K _D ^(Ph1/Ph2) and the ΔCFSE^(Ph1/Ph2) was only observed when geometrical factors, mainly the volume of crystallographic sites and ionic radii of ions competing in the partitioning process, are similar in the respective two paragenetic phases to within 15–20%. In such cases, the ΔCFSE_Fe2+ contribution to K _D ^(Ph1/Ph2) amounts to 0.1 to 0.2 log K _D per 100 cm⁻¹ΔCFSE. The conclusion is that ΔCFSE_Fe2+ plays only a secondary role after geometrical factors, in the partitioning behaviour of Fe²⁺. The reason for this is seen in the facts that, compared to the 3d ³-configured Cr³⁺ ion, CFSE of the 3d ⁶-configured Fe²⁺ amounts only to 20–25%, and that the former ion enters only octahedral sites with similar geometrical properties in the paragenetic mineral phases. Received: 17 November 1998 / Accepted: 28 June 1999     

Mutual replacement reactions in alkali feldspars II: trace element partitioning and geothermometry

Perthitic alkali feldspar primocrysts in layered syenites in the Klokken intrusion in South Greenland, underwent dissolution–reprecipitation reactions in a circulating post-magmatic aqueous fluid at ~450°C, and are to a large degree pseudomorphs. These ‘mutual replacement’ reactions provide a perfect natural experiment with which to study trace element partitioning between sodium and potassium feldspars growing simultaneously. The reactant ‘phase’ was a cryptoperthitic feldspar consisting of low albite and low microcline in a coherent sub-μm ‘braid’ intergrowth and the product phases were ‘strain-free’ incoherent subgrains of low albite and low microcline forming microporous patch perthites on scales up to 200 μm. The driving force for the reaction was reduction of coherency strain energy. The mechanisms of this process are described in Part I. Five mixed braid perthite–patch perthite crystals were analysed for major and trace elements using laser ablation-inductively coupled plasma mass spectrometry with a 19 μm beam diameter. This gave bulk analyses of the braid texture, which were in the range Ab_73–54Or_45–27An_4.3–0.8, but could resolve Ab- and Or-rich patches in patch perthite. The major element bulk compositions of the crystals were retained during the replacement reactions. Major components in patches plot on tielines in the Ab–Or–An ternary system that pass through or very close to the parent braid perthite composition and indicate local equilibrium on the scale of a few tens of mm. Many trace elements, including REE, were lost to the fluid during the deuteric reactions, but the effect is large only for Fe and Ti. Cs, Pb and Sr were added to some crystals. Plots of log distribution coefficient D for Rb, Ba, Pb, Eu²⁺, La and Ce between Or- and Ab-rich patches against ionic radius are straight lines, assuming eightfold coordination, and to a first approximation are independent of ionic charge. K also lies on these lines, and the smaller ions Na and Ca lie close to them. The best linear fits were obtained using ionic radii for ^[8]K and ^[8]Ca, but there is ambiguity as to whether ^[7]Na or ^[5]Na is most appropriate. The linear relationship shows that the listed trace elements are in the feldspar M-site rather than in inclusions. Tl is in M although an exact D could not be obtained. The very large Cs ion partitions strongly into the Or-rich phase but its D value appears to be less than predicted by extrapolation. The near-linearity arises because partitioning is occurring between two solids into sites which have similar Young’s moduli, so that the parabolas that normally represent trace element partitioning between crystals and liquids (which have negligible shear strength) approximately cancel out. Ga and Be are in T-sites, as well as some of the Fe and Ti present, although part is in oxide inclusions. The site of Sc is unclear, but if structural it is likely to be T. Partitioning on M-sites is a potential geothermometer but because the effective size of the irregular M-site is defined by its K and (Na + Ca) contents, which are controlled by ternary solvus relationships, its calibration is not independent of conventional two-feldspar geothermometers. Trace elements may however provide a useful means of confirming that feldspar pairs are in equilibrium, and of recognising feldspar intergrowths produced by non-isochemical replacement rather than exsolution. Two-feldspar geothermometry for the ternary phases in the low-albite microcline patch perthites gives temperatures above the stability range of microcline, markedly so if a correction is made for Si–Al ordering. This is probably because current geothermometers are too sensitive to low concentrations of An in ordered Or-rich feldspars. This interpretation is supported by two-feldspar assemblages growing at known temperatures in geothermal systems and sedimentary basins. This paper and the earlier Part I are dedicated in the memory of J. V. Smith and W. L. Brown, both of whom died in 2007, in acknowledgement of their unrivalled contributions to the study of the feldspar minerals over more than half a century. An erratum to this article can be found at     

Carbonate-rich melt infiltration in peridotite xenoliths from the Eurasian–North American modern plate boundary (Chersky Range,Yakutia)

A suite of mainly spinel peridotite and subordinate pyroxenite xenoliths and megacrysts were studied in detail, enabling us to characterize upper mantle conditions and processes beneath the modern North American–Eurasian continental plate boundary. The samples were collected from 37-Ma-old basanites cropping out in the Main Collision Belt of the Chersky Range, Yakutia Republic (Russian Far East). The spinel lherzolites reflect a mantle sequence, equilibrated at temperatures of 890–1,025 °C at pressures of 1.1–2 GPa, with melt extraction estimated to be around 2–6 %. The spinel harzburgites are characterized by lower P–T equilibration conditions and estimated melt extraction up to 12 %. Minor cryptic metasomatic processes are recorded in the clinopyroxene trace elements, revealing that percolating hydrous fluid-rich melts and basaltic melts affected the peridotites. One of the lherzolites preserves a unique melt droplet with primary dolomite in perfect phase contact with Na-rich aluminosilicate glass and sodalite. On the basis of the well-constrained P–T frame of the xenolith suite, as well as the rigorously documented melt extraction and metasomatic history of this upper mantle section, we discuss how a carbonated silicate melt infiltrated the lherzolite at depth and differentiated into an immiscible carbonate and silicate liquid shortly before the xenolith was transported to the surface by the host basalt. Decreasing temperatures triggered crystallization of primary dolomite from the carbonate melt fraction and sodalite as well as quenched glass from the Na-rich aluminosilicate melt fraction. Rapid entrainment and transport to the Earth’s surface prevented decarbonatization processes as well as reaction phenomena with the host lherzolite, preserving this exceptional snapshot of upper mantle carbonatization and liquid immiscibility.     

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     

Trace element partitioning between orthopyroxene and anhydrous silicate melt on the lherzolite solidus from 1.1 to 3.2 GPa and 1,230 to 1,535°C in the model system Na<Subscript>2</Subscript>O–CaO–MgO–Al<Subscript>2</Subscript>O<Subscript>3</Subscript>–SiO<Subscript>2</Subscript>

We determined experimentally the Nernst distribution coefficient between orthopyroxene and anhydrous silicate melt for trace elements i in the system Na₂O–CaO–MgO–Al₂O₃–SiO₂ (NCMAS) along the dry model lherzolite solidus from 1.1 GPa/1,230°C up to 3.2 GPa/1,535°C in a piston cylinder apparatus. Major and trace element composition of melt and orthopyroxene were determined with a combination of electron microprobe and ion probe analyses. We provide partitioning data for trace elements Li, Be, B, K, Sc, Ti, V, Cr, Co, Ni, Rb, Sr, Y, Zr, Nb, Cs, Ba, La, Ce, Sm, Nd, Yb, Lu, Hf, Ta, Pb, U, and Th. The melts were chosen to be boninitic at 1.1 and 2.0 GPa, picritic at 2.3 GPa and komatiitic at 2.7 and 3.2 GPa. Orthopyroxene is Tschermakitic with 8 mol% Mg-Tschermaks MgAl[AlSiO₆] at 1.1 GPa while at higher pressure it has 18–20 mol%. The rare earth elements show a continuous, significant increase in compatibility with decreasing ionic radius from D _La^opx−melt ∼ 0.0008 to D _Lu^opx−melt ∼ 0.15. For the high-field-strength elements compatibility increases from D _Th^opx−melt ∼ 0.001 through D _Nb^opx−melt ∼ 0.0015, D _U^opx−melt ∼ 0.002, D _Ta^opx−melt ∼ 0.005, D _Zr^opx−melt ∼ 0.02 and D _Hf^opx−melt ∼ 0.04 to D _Ti^opx−melt ∼ 0.14. From mathematical and graphical fits we determined best-fit values for D ₀^M1, D ₀^M2, r ₀^M1, r ₀^M2, E ₀^M1, and E ₀^M2 for the two different M sites in orthopyroxene according to the lattice strain model and calculated the intracrystalline distribution between M1 and M2. Our data indicate extreme intracrystalline fractionation for most elements in orthopyroxene; for the divalent cations D _i ^M2−M1 varies by three orders of magnitude between D _Co^M2−M1 = 0.00098–0.00919 and D _Ba^M2−M1 = 2.3–28. Trivalent cations Al and Cr almost exclusively substitute on M1 while the other trivalent cations substitute on M2; D _La^M2−M1 reaches extreme values between 6.5 × 10⁷ and 1.4 × 10¹⁶. Tetravalent cations Ti, Hf, and Zr almost exclusively substitute on M1 while U and Th exclusively substitute on M2. Our new comprehensive data set can be used for polybaric-polythermal melting models along the Earth’s mantle solidus. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Experimental and computational study of trace element distribution between orthopyroxene and anhydrous silicate melt: substitution mechanisms and the effect of iron

Although orthopyroxene (Opx) is present during a wide range of magmatic differentiation processes in the terrestrial and lunar mantle, its effect on melt trace element contents is not well quantified. We present results of a combined experimental and computational study of trace element partitioning between Opx and anhydrous silicate melts. Experiments were performed in air at atmospheric pressure and temperatures ranging from 1,326 to 1,420°C in the system CaO–MgO–Al₂O₃–SiO₂ and subsystem CaO–MgO–SiO₂. We provide experimental partition coefficients for a wide range of trace elements (large ion lithophile: Li, Be, B, K, Rb, Sr, Cs, Ba, Th, U; rare earth elements, REE: La, Ce, Nd, Sm, Y, Yb, Lu; high field strength: Zr, Nb, Hf, Ta, Ti; transition metals: Sc, V, Cr, Co) for use in petrogenetic modelling. REE partition coefficients increase from $ D_{\text{La}}^{{\text{Opx}} {\hbox{-}} {\text{melt}}} \sim 0.0005 Although orthopyroxene (Opx) is present during a wide range of magmatic differentiation processes in the terrestrial and lunar mantle, its effect on melt trace element contents is not well quantified. We present results of a combined experimental and computational study of trace element partitioning between Opx and anhydrous silicate melts. Experiments were performed in air at atmospheric pressure and temperatures ranging from 1,326 to 1,420°C in the system CaO–MgO–Al₂O₃–SiO₂ and subsystem CaO–MgO–SiO₂. We provide experimental partition coefficients for a wide range of trace elements (large ion lithophile: Li, Be, B, K, Rb, Sr, Cs, Ba, Th, U; rare earth elements, REE: La, Ce, Nd, Sm, Y, Yb, Lu; high field strength: Zr, Nb, Hf, Ta, Ti; transition metals: Sc, V, Cr, Co) for use in petrogenetic modelling. REE partition coefficients increase from $ D_{\text{La}}^{{\text{Opx}} {\hbox{-}} {\text{melt}}} \sim 0.0005 $ D_{\text{La}}^{{\text{Opx}} {\hbox{-}} {\text{melt}}} \sim 0.0005 to $ D_{\text{Lu}}^{{{\text{Opx}} {\hbox{-}} {\text{melt}}}} \sim 0.109 $ D_{\text{Lu}}^{{{\text{Opx}} {\hbox{-}} {\text{melt}}}} \sim 0.109 , D values for highly charged elements vary from $ D_{\text{Th}}^{{{\text{Opx}} {\hbox{-}} {\text{melt}}}} \sim 0.0026 $ D_{\text{Th}}^{{{\text{Opx}} {\hbox{-}} {\text{melt}}}} \sim 0.0026 through $ D_{\text{Nb}}^{{{\text{Opx}} {\hbox{-}} {\text{melt}}}} \sim 0.0033 $ D_{\text{Nb}}^{{{\text{Opx}} {\hbox{-}} {\text{melt}}}} \sim 0.0033 and $ D_{\text{U}}^{{{\text{Opx}} {\hbox{-}} {\text{melt}}}} \sim 0.0066 $ D_{\text{U}}^{{{\text{Opx}} {\hbox{-}} {\text{melt}}}} \sim 0.0066 to $ D_{\text{Ti}}^{{\text{Opx}} {\hbox{-}} {\text{melt}}} \sim 0.058 $ D_{\text{Ti}}^{{\text{Opx}} {\hbox{-}} {\text{melt}}} \sim 0.058 , and are all virtually independent of temperature. Cr and Co are the only compatible trace elements at the studied conditions. To elucidate charge-balancing mechanisms for incorporation of REE into Opx and to assess the possible influence of Fe on Opx-melt partitioning, we compare our experimental results with computer simulations. In these simulations, we examine major and minor trace element incorporation into the end-members enstatite (Mg₂Si₂O₆) and ferrosilite (Fe₂Si₂O₆). Calculated solution energies show that R²⁺ cations are more soluble in Opx than R³⁺ cations of similar size, consistent with experimental partitioning data. In addition, simulations show charge balancing of R³⁺ cations by coupled substitution with Li⁺ on the M1 site that is energetically favoured over coupled substitution involving Al–Si exchange on the tetrahedrally coordinated site. We derived best-fit values for ideal ionic radii r ₀, maximum partition coefficients D ₀, and apparent Young’s moduli E for substitutions onto the Opx M1 and M2 sites. Experimental r ₀ values for R³⁺ substitutions are 0.66–0.67 ? for M1 and 0.82–0.87 ? for M2. Simulations for enstatite result in r ₀ = 0.71–0.73 ? for M1 and ~0.79–0.87 ? for M2. Ferrosilite r ₀ values are systematically larger by ~0.05 ? for both M1 and M2. The latter is opposite to experimental literature data, which appear to show a slight decrease in $ r_{0}^{{{\text{M}}2}} $ r_{0}^{{{\text{M}}2}} in the presence of Fe. Additional systematic studies in Fe-bearing systems are required to resolve this inconsistency and to develop predictive Opx-melt partitioning models for use in terrestrial and lunar magmatic differentiation models.     

Trace-element partitioning between garnet,clinopyroxene and Fe-rich picritic melts at 3 to 7 GPa

We have determined mineral-melt partition coefficients (D values) for 20 trace elements in garnet-pyroxenite run products, generated in 3 to 7 GPa, 1,425–1,750°C experiments on a high-Fe mantle melt (97SB68) from the Paraná-Etendeka continental-flood-basalt (CFB) province. D values for both garnet (∼Py₆₃Al₂₅Gr₁₂) and clinopyroxene (∼Ca_0.2Mg_0.6Fe_0.2Si₂O₆) show a large variation with temperature but are less dependent on pressure. At 3 GPa, D ^cpx/liq values for pyroxenes in garnet-pyroxenite run products are generally lower than those reported from Ca-rich pyroxenes generated in melting experiments on eclogites and basalts (∼Ca_0.3–0.5Mg_0.3–0.6Fe_0.07–0.2Si₂O₆) but higher than those for Ca-poor pyroxenes from peridotites (∼Ca_0.2Mg_0.7Fe_0.1Si₂O₆). D ^grt/liq values for light and heavy rare-earth elements are ≤0.07 and >0.8, respectively, and are similar to those for peridotitic garnets that have comparable grossular but higher pyrope contents (Py_70–88Al_l7–20Gr_8–14). 97SB68 D _LREE^grt/liq values are higher and D _HREE^grt/liq values lower than those for eclogitic garnets which generally have higher grossular contents but lower pyrope contents (Py_20–70Al_10–50Gr_10–55). D values agree with those predicted by lattice strain modelling and suggest that equilibrium was closely approached for all of our experimental runs. Correlations of D values with lattice-strain parameters and major-element contents suggest that the wollastonite component and pyrope:grossular ratio exert major controls on 97SB68 clinopyroxene and garnet partitioning, respectively. These are controlled by the prevailing pressure and temperature conditions for a given bulk-composition. The composition of co-existing melt was found to have a relatively minor effect on 97SB68 D values. The variations in D values displayed by different mantle lithologies are subtle and our study confirms previous investigations which have suggested that the modal proportions of garnet and clinopyroxene are by far the most influential factor in determining incompatible trace-element concentrations in mantle melts. The trace-element partition coefficients we have determined may be used to place high-pressure constraints on garnet-pyroxenite melting models.     

Carbonatite melt inclusions in coexisting magnetite,apatite and monticellite in Kerimasi calciocarbonatite,Tanzania: melt evolution and petrogenesis

Kerimasi calciocarbonatite consists principally of calcite together with lesser apatite, magnetite, and monticellite. Calcite hosts fluid and S-bearing Na–K–Ca-carbonate inclusions. Carbonatite melt and fluid inclusions occur in apatite and magnetite, and silicate melt inclusions in magnetite. This study presents statistically significant compositional data for quenched S- and P-bearing, Ca-alkali-rich carbonatite melt inclusions in magnetite and apatite. Magnetite-hosted silicate melts are peralkaline with normative sodium-metasilicate. On the basis of our microthermometric results on apatite-hosted melt inclusions and forsterite–monticellite phase relationships, temperatures of the early stage of magma evolution are estimated to be 900–1,000°C. At this time three immiscible liquid phases coexisted: (1) a Ca-rich, P-, S- and alkali-bearing carbonatite melt, (2) a Mg- and Fe-rich, peralkaline silicate melt, and (3) a C–O–H–S-alkali fluid. During the development of coexisting carbonatite and silicate melts, the Si/Al and Mg/Fe ratio of the silicate melt decreased with contemporaneous increase in alkalis due to olivine fractionation, whereas the alkali content of the carbonatite melt increased with concomitant decrease in CaO resulting from calcite fractionation. Overall the peralkalinity of the bulk composition of the immiscible melts increased, resulting in a decrease in the size of the miscibility gap in the pseudoquaternary system studied. Inclusion data indicate the formation of a carbonatite magma that is extremely enriched in alkalis with a composition similar to that of Oldoinyo Lengai natrocarbonatite. In contrast to the bulk compositions of calciocarbonatite rocks, the melt inclusions investigated contain significant amount of alkalis (Na₂O + K₂O) that is at least 5–10 wt%. The compositions of carbonatite melt inclusions are considered as being better representatives of parental magma composition than those of any bulk rock.     

Laser Ablation ICPMS study of trace element partitioning between plagioclase and basaltic melts: an experimental approach

Plagioclase-melt partition coefficients (D) for 34 trace elements at natural concentration levels were determined experimentally in a natural MORB composition at atmospheric pressure using thin Pt-wire loops. Experiments were carried out at three temperatures (1,220, 1,200, and 1,180°C), and at three different oxygen fugacities (fO₂ = IW, QFM, air) in order to assess the effect of fO₂ on the partitioning of elements with multiple valence (Fe, Eu, Cr). Run products were analyzed by laser-ablation ICP-MS. Most trace element Ds increase slightly as temperature decreases, except for D _Zr, D _Fe, D _Eu and D _Cr that vary systematically with fO₂. Applying the Lattice Strain Model to our data suggests the presence of Fe²⁺ entirely in the octahedral site at highly to moderate reducing conditions, while Fe³⁺ was assigned wholly to the tetrahedral site of the plagioclase structure. Furthermore, we provide a new quantitative framework for understanding the partitioning behaviour of Eu, which occurs as both 2+ and 3+ cations, depending on fO₂and confirm the greater compatibility of Eu²⁺, which has an ionic radius similar to Sr, relative to Eu³⁺ in plagioclase and the higher Eu²⁺/ Eu³⁺ under reducing conditions. For petrogenetic basaltic processes, a combined fractionation of Eu²⁺–Sr and Fe–Mg by plagioclase has considerable potential as an oxybarometer for natural magmatic rocks.     

H<Subscript>2</Subscript>O storage capacity of MgSiO<Subscript>3</Subscript> clinoenstatite at 8–13 GPa, 1,100–1,400°C

We present H₂O analyses of MgSiO₃ pyroxene crystals quenched from hydrous conditions in the presence of olivine or wadsleyite at 8–13.4 GPa and 1,100–1,400°C. Raman spectroscopy shows that all pyroxenes have low clinoenstatite structure, which we infer to indicate that the crystals were high clinoenstatite (C2/c) during conditions of synthesis. H₂O analyses were performed by secondary ion mass spectrometry and confirmed by unpolarized Fourier transform infrared spectroscopy on randomly oriented crystals. Measured H₂O concentrations increase with pressure and range from 0.08 wt.% H₂O at 8 GPa and 1,300°C up to 0.67 wt.% at 13.4 GPa and 1,300°C. At fixed pressure, H₂O storage capacity diminishes with increasing temperature and the magnitude of this effect increases with pressure. This trend, which we attribute to diminishing activity of H₂O in coexisting fluids as the proportion of dissolved silicate increases, is opposite to that observed previously at low pressure. We observe clinoenstatite 1.4 GPa below the pressure stability of clinoenstatite under nominally dry conditions. This stabilization of clinoenstatite relative to orthoenstatite under hydrous conditions is likely owing to preferential substitution of H₂O into the high clinoenstatite polymorph. At 8–11 GPa and 1,200–1,400°C, observed H₂O partitioning between olivine and clinoenstatite gives values of D ^ol/CEn between 0.65 and 0.87. At 13 GPa and 1,300°C, partitioning between wadsleyite and clinoenstatite, D ^wd/CEn, gives a value of 2.8 ± 0.4.     

Quantifying garnet-melt trace element partitioning using lattice-strain theory: assessment of statistically significant controls and a new predictive model

As a complement to our efforts to update and revise the thermodynamic basis for predicting garnet-melt trace element partitioning using lattice-strain theory (van Westrenen and Draper in Contrib Mineral Petrol, this issue), we have performed detailed statistical evaluations of possible correlations between intensive and extensive variables and experimentally determined garnet-melt partitioning values for trivalent cations (rare earth elements, Y, and Sc) entering the dodecahedral garnet X-site. We applied these evaluations to a database containing over 300 partition coefficient determinations, compiled both from literature values and from our own work designed in part to expand that database. Available data include partitioning measurements in ultramafic to basaltic to intermediate bulk compositions, and recent studies in Fe-rich systems relevant to extraterrestrial petrogenesis, at pressures sufficiently high such that a significant component of majorite, the high-pressure form of garnet, is present. Through the application of lattice-strain theory, we obtained best-fit values for the ideal ionic radius of the dodecahedral garnet X-site, r ₀(3+), its apparent Young’s modulus E(3+), and the strain-free partition coefficient D ₀(3+) for a fictive REE element J of ionic radius r ₀(3+). Resulting values of E, D ₀, and r ₀ were used in multiple linear regressions involving sixteen variables that reflect the possible influence of garnet composition and stoichiometry, melt composition and structure, major-element partitioning, pressure, and temperature. We find no statistically significant correlations between fitted r ₀ and E values and any combination of variables. However, a highly robust correlation between fitted D ₀ and garnet-melt Fe–Mg exchange and D _Mg is identified. The identification of more explicit melt-compositional influence is a first for this type of predictive modeling. We combine this statistically-derived expression for predicting D ₀ with the new expressions for predicting E and r ₀ outlined in the first of our pair of companion papers into an updated set of formulae that use easy-to-measure quantities (e.g. garnet composition, pressure, temperature) to predict variations in E, r ₀, and D ₀. These values are used in turn to calculate D values for those garnets. The updated model substantially improves upon a previous model (van Westrenen et al. in Contrib Mineral Petrol 142:219–234, 2001), and accounts well for trivalent cation partitioning in nominally anhydrous systems up to at least 15 GPa, including for eclogitic bulk compositions and for Fe-rich systems appropriate to magmagenesis on the Moon and Mars. The new model is slightly less successful in predicting partitioning with strongly majoritic garnets, although the mismatch is much less than with the original 2001 model. Although it also improves upon the 2001 model in predicting partitioning in hydrous systems, the mismatch between model and observation is still unacceptably large. The same statistical tools were applied in an attempt to predict tetravalent partitioning as well, because lattice-strain based techniques are not applicable to such partitioning. However, no statistically significant predictive relationships emerged from that effort. Our analyses show that future efforts should focus on filling the gap in partitioning data between ∼10 and 25 GPa to evaluate more closely the gradual transition of garnet to majorite, and on systematically expanding the hydrous partitioning database to allow extension of our model to water-bearing systems.     

The Partitioning of Tungsten bwtween Aqueous Fluids and Silicate Melts

An experimental study has been carried out to determine the partition coefficients of tungsten between aqueous fluids and granitic melts at 800 °C and 1.5 kb with natural granite as the starting material. The effects of the solutions on the partition coefficients of tungsten show a sequence of P > CO ₃ ²⁻ > B > H₂O. The effects are limited (generallyK _D < 0.3) and the tungsten shows a preferential trend toward the melt over the aqueous fluid. The value ofK _D increases with increasing concentration of phosphorus; theK _D increases first and then reduces with the concentration of CO ₃ ²⁻ when temperature decreases, theK _D between the solution of CO ₃ ²⁻ and the silicate melt increases, and that between the solution of B₄O ₇ ²⁻ and the silicate melt decreases. The partition coefficients of phosphorus and sodium between fluids and silicate melts have been calculated from the concentrations of the elements in the melts. TheK _D value for phosphorus is 0.38 and that for sodium is 0.56. Evidence shows that the elements tend to become richer and richer in the melts.     

Fluid transfer of gold,palladium, and rare earth elements and genesis of ore occurrences in the Subpolar Urals

To estimate the behavior of Au, Pd, REE, and Y in magmatic and postmagmatic processes, a series of experimental studies on the solubility of noble metals and REE in magma, magmatic fluid, and hydrothermal solutions has been performed in wide temperature and pressure ranges (300–400°C, 860–1350°C; 1–14 kbar). The coefficients of Au and Pd partitioning (D ^F/L) between fluid and tholeiitic melt have been determined. Depending on P, T, and the composition of the system, they vary from 1 to 11 for Au and 0.02 to 1 for Pd. The phase solubility technique was used to determine Au and Pd solubility in hydrothermal fluid. The effects of temperature, composition, and fluid acidity on Au and Pd solubility have been estimated. The high solubility of these metals in aqueous chloride solutions has been established for both Au (28–803 mg/kg at T = 300°C, 305–1123 mg/kg at T = 350°C, and 330–1400 mg/kg at T = 400°C) and Pd (40–126 mg/kg at T = 300°C, 62–152 mg/kg at T = 350°C, and 20–210 mg/kg at T = 400°C). The coefficients of REE and Y partitioning (D ^F/L) between fluid and tholeiitic or alkaline melts have been determined. They vary from 0.00n to 2 depending on P, T, and fluid composition. The experimental data on Au and Pd solubility in solutions and magmatic fluids and the wide variation of REE D ^F/L between fluid and melt show that magmatic and hydrothermal fluids are efficient agents of Au, Pd, and REE transfer and fractionation. The obtained experimental data were used for elucidating sources of fluids and their role in the genesis of Au-Pd-REE occurrences in the Subpolar Urals.     

Garnet + liquid equilibrium

New experiments were performed to determine saturation conditions for garnet and silicate liquid. Starting compositions were natural basalt powders ranging from komatiite to nephelinite, which were partially melted at pressures between 25 and 100 kbar. Rounded grains of natural pyrope or grossular were added to some experiments to induce garnet saturation, and to aid the segregation of liquid pools for microprobe analysis. Simple expressions describing K _eq as a function of P, T and liquid composition were calibrated by linear least squares analysis of the data from this, and other, studies. Since garnets do not often occur as phenocrysts, equations were designed to predict garnet compositions when P, T and a silicate liquid composition are given. The regression data have a pressure range of 20–270 kbar, and compositions as diverse as nephelinite and komatiite. These models should thus apply to a broad range of geological problems. The majorite component in garnet was found to increase with increasing P, but compositional effects are also important. A garnet saturation surface applied to liquids with chondritic compositions shows that such liquids crystallize garnet with Mj contents of 0.27–0.42 at 200 kbar. Models of Earth differentiation thus need to account not only for fractionation of majorite, but also for Fe-, Ca-, Na- and Ti-bearing garnet components, which occur in non-trivial quantities at high pressure. Since many models of igneous petrogenesis rely on mineral-melt partition coefficients for the minor elements Na, Ti, and Cr, partition coefficients for these elements were also examined. The K _d ^gar/liq for Na was found to be P-sensitive; Na contents of basalts may thus potentially yield information regarding depths of partial melting. Received : 28 May 1997 / 25 November 1997     

Experimental determination of REE partition coefficients between zircon,garnet and melt: a key to understanding high‐T crustal processes

The partitioning of rare earth elements (REE) between zircon, garnet and silicate melt was determined using synthetic compositions designed to represent partial melts formed in the lower crust during anatexis. The experiments, performed using internally heated gas pressure vessels at 7 kbar and 900–1000 °C, represent equilibrium partitioning of the middle to heavy REE between zircon and garnet during high‐grade metamorphism in the mid to lower crust. The D_REE (zircon/garnet) values show a clear partitioning signature close to unity from Gd to Lu. Because the light REE have low concentrations in both minerals, values are calculated from strain modelling of the middle to heavy REE experimental data; these results show that zircon is favoured over garnet by up to two orders of magnitude. The resulting general concave‐up shape to the partitioning pattern across the REE reflects the preferential incorporation of middle REE into garnet, with D_Gd (zircon/garnet) ranging from 0.7 to 1.1, D_Ho (zircon/garnet) from 0.4 to 0.7 and D_Lu (zircon/garnet) from 0.6 to 1.3. There is no significant temperature dependence in the zircon–garnet REE partitioning at 7 kbar and 900–1000 °C, suggesting that these values can be applied to the interpretation of zircon–garnet equilibrium and timing relationships in the ultrahigh‐T metamorphism of low‐Ca pelitic and aluminous granulites.     

Experimental investigation of large-ion-lithophile-element-, high-field-strength-element- and rare-earth-element-partitioning between calcic amphibole and basaltic melt: the effects of pressure and oxygen fugacity

The effects of pressure and oxygen fugacity (fO₂) on trace element partitioning between pargasitic amphibole and alkali-basaltic melts have been determined at pressures from 1.5 to 2.5 GPa and oxygen fugacities at 2 log units above and below the nickel–nickel oxide buffer. Amphibole crystallization experiments were performed in a piston cylinder apparatus and partition coefficients between amphibole and quenched melt of large-ion-lithophile elements (LILE: Rb, Sr, Ba), high-field-strength elements (HFSE: Zr, Nb, Ta, Hf, U, Th) and rare-earth elements (REE: La to Lu; +Y) were measured with a LASER ablation inductively coupled plasma – mass spectrometer. Increasing pressure from 1.5 to 2.5 GPa at similar temperatures and approximately constant fO₂ increases D _Rb but decreases D _Zr and D _Hf and D _REE (D _La, D _Ce, D _Pr). An empirical relationship was observed between D _Zr and (Ti/Al)_M2 in the amphibole, which can be described by: