incorporation in magmas at depth: changes in melt local environment and the influence partitioning

The structure of two Lu doped (4000 ppm) model end member silicate liquids, a highly polymerised haplogranite (Si-Al-Na-K-O) and a less polymerised anorthite-diopside (Si-Al-Mg-Ca-O), have been studied up to 8 GPa using in situ x-ray diffraction techniques. The results are the first to identify trace rare Earth element incorporation in silicate melts at high pressure. At pressures below 5 GPa, the bonding environment of Lu-O was found to be dependent on composition with coordination number CN Lu − O = 8 and bond distance r Lu − O = 2 . 36 ˚A in the haplogranite melt, decreasing to CN Lu − O = 6 and r Lu − O = 2 . 29 ˚A in the anorthite-diopside melt. This compositional variance in coordination number at low pressure is consistent with observations made for Y-O in glasses at ambient conditions and is coincident with a dramatic increase in the partition coefficients previously observed for rare Earth elements with increasing melt polymerisation. With increasing pressure we find that CN Lu − O and r Lu − O remain constant in the haplogranite melt. However, an abrupt change in both Lu-O coordination and bond distance is observed at 5 GPa in the anorthite-diopside melt, with CN Lu − O increasing from 6 to 8-fold and r Lu − O from 2.29 to 2.39 ˚A. This occurs over a similar pressure range where a change in the P-dependence in the reported rare Earth element partition coefficients is observed for garnet-, clinopyroxene-, and olivine-melt systems. This work shows that standard models for predicting trace elements at depth must incorporate the effect of pressure-induced structural transformations in the melt in order to realistically predict partitioning behaviour. element studies due to the low concentrations of trace elements ( < 0.1 wt%) required to represent natural systems. X-ray diffraction has the benefit of providing information on both the short and medium range structure of the melt, includ- ing absolute distance and average coordination numbers, without relying on structural models. In this study we employ two end member melt compositions, silicic haplogran- ite and a ‘model basalt’ iron-free anorthite-diopside to monitor structural changes at high pressure and temperature that may affect the incorporation of Lu at 4000 ppm concentration into the

[1]  C. Sanloup Density of magmas at depth , 2016 .

[2]  G. Shen,et al.  Toward comprehensive studies of liquids at high pressures and high temperatures: Combined structure, elastic wave velocity, and viscosity measurements in the Paris–Edinburgh cell , 2014 .

[3]  Mark L. Rivers,et al.  Atomistic insight into viscosity and density of silicate melts under pressure , 2014, Nature Communications.

[4]  Donna Morton,et al.  Structural change in molten basalt at deep mantle conditions , 2013, Nature.

[5]  M. Salanne,et al.  Molecular dynamics simulations of Y in silicate melts and implications for trace element partitioning , 2013 .

[6]  G. Morard,et al.  Magmas trapped at depth and the continental lithosphere-asthenosphere boundary , 2012 .

[7]  T. Hirata,et al.  Element partitioning between olivine and melt up to 10 GPa: Implications for the effect of pressure , 2012 .

[8]  T. Hirata,et al.  Pressure effect on element partitioning between minerals and silicate melt: Melting experiments on basalt up to 20 GPa , 2012 .

[9]  H. Morita,et al.  Relationship between topological order and glass forming ability in densely packed enstatite and forsterite composition glasses , 2011, Proceedings of the National Academy of Sciences.

[10]  L. Stixrude,et al.  First principles molecular dynamics simulations of diopside (CaMgSi2O6) liquid to high pressure , 2011 .

[11]  H. E. Maynard-Casely,et al.  Compressibility change in iron-rich melt and implications for core formation models , 2011 .

[12]  Nobumasa Funamori,et al.  High-pressure structural transformation of SiO 2 glass up to 100 GPa , 2010 .

[13]  A. Barnhoorn,et al.  PVT relation of MgO derived by simultaneous elastic wave velocity and in situ X-ray measurements: A new pressure scale for the mantle transition region , 2010 .

[14]  B. Guillot,et al.  Computer modeling of natural silicate melts: What can we learn from ab initio simulations , 2009, 0907.2530.

[15]  H. O’Neill,et al.  The influence of melt composition on the partitioning of REEs, Y, Sc, Zr and Al between forsterite and melt in the system CMAS , 2008 .

[16]  W. Westrenen,et al.  Quantifying garnet-melt trace element partitioning using lattice-strain theory: new crystal-chemical and thermodynamic constraints , 2007 .

[17]  B. Guillot,et al.  A computer simulation study of natural silicate melts. Part I: Low pressure properties , 2007 .

[18]  S. Klemme,et al.  Rare earth element partitioning between titanite and silicate melts: Henry's law revisited , 2006 .

[19]  G. Gaetani The influence of melt structure on trace element partitioning near the peridotite solidus , 2004 .

[20]  S. Klemme,et al.  Effect of melt composition on the partitioning of trace elements between titanite and silicate melt , 2003 .

[21]  J. Longhi,et al.  Near mantle solidus trace element partitioning at pressures up to 3.4 GPa , 2002 .

[22]  W. Westrenen,et al.  Atomistic simulation of trace element incorporation into garnets—comparison with experimental garnet-melt partitioning data , 2000 .

[23]  D. Dingwell THE GLASS TRANSITION IN HYDROUS GRANITIC MELTS , 1998 .

[24]  B. Wood,et al.  A predictive model for rare earth element partitioning between clinopyroxene and anhydrous silicate melt , 1997 .

[25]  F. Farges Does ZrF “complexation” occur in magmas? , 1996 .

[26]  J. Yarger,et al.  Al Coordination Changes in High-Pressure Aluminosilicate Liquids , 1995, Science.

[27]  H. Mao,et al.  In Situ Determination of the NiAs Phase of FeO at High Pressure and Temperature , 1994, Science.

[28]  B. Wood,et al.  Prediction of crystal–melt partition coefficients from elastic moduli , 1994, Nature.

[29]  H. Keppler,et al.  Pressure-induced coordination changes of transition-metal ions in silicate melts , 1993, Nature.

[30]  P. Barbey,et al.  Effects of H 2 O on liquidus phase relations in the haplogranite system at 2 and 5 kbar , 1992 .

[31]  S. Hart,et al.  Nickel partitioning between olivine and silicate melt , 1978 .

[32]  P. C. Hess,et al.  Implications of liquid-liquid distribution coefficients to mineral-liquid partitioning , 1978 .

[33]  R. D. Shannon Revised effective ionic radii and systematic studies of interatomic distances in halides and chalcogenides , 1976 .

[34]  J. H. Hubbell,et al.  Atomic form factors, incoherent scattering functions, and photon scattering cross sections , 1975 .

[35]  F. Hajdu Revised parameters of the analytic fits for coherent and incoherent scattered X‐ray intensities of the first 36 atoms , 1972 .

[36]  Vnco,et al.  Solubility mechanisms of H 2 O in silicate melts at high pressures and temperatures : a Raman spectroscopic study , 2007 .

[37]  P. Beattie On the occurrence of apparent non-Henry's Law behaviour in experimental partitioning studies , 1993 .

[38]  John Ziman,et al.  A theory of the electrical properties of liquid metals , 1965 .

[39]  R. Hosemann,et al.  Direct analysis of diffraction by matter , 1962 .