High-pressure phase equilibria and element partitioning experiments on Apollo 15 green C picritic glass: Implications for the role of garnet in the deep lunar interior

We report results of nominally anhydrous near-liquidus experiments on a synthetic analog to very low-titanium Apollo 15 green C lunar picritic glass from ∼2 to 5 GPa. Apollo 15 green C glass (A15C) is saturated with garnet and pyroxene on the liquidus at ∼3 GPa. However, such an assemblage is unlikely to represent the lunar-mantle source region for this glass, and instead an olivine + orthopyroxene-dominated source is favored, in accord with earlier lower-pressure experiments on A15C. Near-liquidus garnet has a slight but significant majorite component at ∼5 GPa in this iron-rich bulk composition, as expected from our previous work in ordinary-chondritic bulk compositions. Ion microprobe measurements of partitioning of Sr, Ba, Sc, Nd, Sm, Dy, Yb, Y, Zr, Hf, and Th between garnet and coexisting melt in these experiments are the first garnet partition coefficients (D values) available that are directly relevant to lunar compositions. D values for these garnets differ significantly compared to D values for garnets grown in more magnesian, terrestrial bulk compositions, which until now are all that have been available in modeling the possible role of garnet in the lunar interior. For example, D values for heavy rare earth elements are lower than are those from terrestrial basaltic systems. These partitioning values are well-described by the lattice-strain partitioning model, but predictive relationships for garnet partitioning using that model fail to match the measured values, as was the case in our earlier work on chondritic compositions. Using our new D values in place of the “terrestrial” values in a variety of models of lunar petrogenesis, we suggest that garnet is unlikely to be present in the source regions for very titanium-poor lunar liquids despite its appearance on the liquidus of A15C.     

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.     

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

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

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.     

Trace element incorporation into pyrope–grossular solid solutions: an atomistic simulation study

 We have performed atomistic computer simulations on trace element incorporation into the divalent dodecahedral X-sites of pyrope (Py — Mg₃Al₂Si₃O₁₂) – grossular (Gr — Ca₃Al₂Si₃O₁₂) solid solutions. An ionic model and the Mott–Littleton two-region approach to defect energies were used to calculate the energetics of substitution by a range of divalent trace-elements and of charge-balanced substitution by trivalent ions in the static limit. Results are compared with experimental high-temperature, high-pressure garnet-melt trace element partitioning data obtained for the same garnet solid solution to refine our understanding of the factors controlling element partitioning into solid solutions. Defect energies (U _def,f), relaxation (lattice strain) energies (U _rel), and solution energies (U _sol) were derived using two different approaches. One approach assumes the presence of one type of hybrid X-site with properties intermediate between pure Mg and Ca sites, and the other assumes discrete Mg and Ca X-sites, and thus two distinct cation sublattices. The hybrid model is shown to be inadequate, since it averages out local distortions in the garnet structure. The discrete model results suggest trace elements are more soluble in Py₅₀Gy₅₀ than in either end-member compound. Physically this is due to small changes in size of the X-sites and the removal of unfavourable interactions between third nearest neighbours of the same size. Surprisingly, depending on the local order, large trace element cations may substitute for Mg²⁺ and small trace elements for Ca²⁺ in Py₅₀Gr₅₀. These solubilities provide an explanation for the anomalous trace-element partitioning behaviour along the pyrope–grossular join observed experimentally. Received: 27 January 2000 / Accepted: 14 February 2003     

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.     

High-precision high field strength element partitioning between garnet, amphibole and alkaline melt from Kakanui, New Zealand

The high field strength elements (HFSE: Zr, Hf, Nb, Ta, and W) are an important group of chemical tracers that are increasingly used to investigate magmatic differentiation processes. Successful modeling of these processes requires the availability of accurate mineral-melt partition coefficients (D). To date, these have largely been determined by ion microprobe or laser ablation-ICP-MS analyses of the run products of high-pressure, high-temperature experiments. Since HFSE are (highly) incompatible, relatively immobile, high-charge, and difficult to ionize, these experiments and their analysis are challenging. Here we explore whether high-precision analyses of natural mineral-melt systems can provide additional constraints on HFSE partitioning.The HFSE concentrations in natural garnet and amphibole and their alkaline host melt from Kakanui, New Zealand are determined with high precision isotope dilution on a multi-collector-ICP-MS. Major and trace element compositions combined with Lu-Hf isotopic systematics and detailed petrographic sample analysis are used to assess mineral-melt equilibrium and to provide context for the HFSE D measurements. The whole-rock nephelinite, ∼1 mm sized amphiboles in the nephelinite, and garnet megacrysts have similar initial Hf isotope ratios with a mean initial ¹⁷⁶Hf/¹⁷⁷Hf_{(34 Ma)} = 0.282900 ± 0.000026 (2σ). In contrast, the amphibole megacrysts are isotopically distinct (¹⁷⁶Hf/¹⁷⁷Hf_{(34 Ma)} = 0.282830 ± 0.000011). Rare earth element D values for garnet megacryst-nephelinite melt and ∼1 mm amphibole-nephelinite melt plotted as a function of ionic radii show classic near-parabolic trends that are in excellent agreement with crystal lattice-strain models. These observations are consistent with equilibrium between the whole-rock nephelinite, the ∼1 mm amphibole grains within the nephelinite and the garnet megacrysts.High-precision isotope dilution results for Zr and Hf in garnet (D_Zr = 0.220 ± 0.007 and D_Hf = 0.216 ± 0.005 [2σ]), and for all HFSE in amphibole are consistent with previous experimental findings. However, our measurements for Nb and Ta in garnet (D_Nb = 0.0007 ± 0.0001 and D_Ta = 0.0011 ± 0.0006 [2σ]) show that conventional methods may overestimate Nb and Ta concentrations, thereby overestimating both Nb and Ta absolute D values for garnet by up to 3 orders of magnitude and underestimating D_Nb/D_Ta by greater than a factor of 100. As a consequence, the role of residual garnet in imposing Nb/Ta fractionation may be less important than previously thought. Moreover, garnet D_Hf/D_W = 17 and D_Nb/D_Zr = 0.003 imply fractionation of Hf from W and Nb from Zr upon garnet crystallization, which may have influenced short-lived ¹⁸²Hf-¹⁸²W and ⁹²Nb-⁹²Zr isotopic systems in Hadean time.     

Effect of Fe on garnet–melt trace element partitioning: experiments in FCMAS and quantification of crystal-chemical controls in natural systems

Garnet–melt trace element partitioning experiments were performed in the system FeO–CaO–MgO–Al₂O₃–SiO₂ (FCMAS) at 3 GPa and 1540°C, aimed specifically at studying the effect of garnet Fe²⁺ content on partition coefficients (D^Grt/Melt). D^Grt/Melt, measured by SIMS, for trivalent elements entering the garnet X-site show a small but significant dependence on garnet almandine content. This dependence is rationalised using 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], which describes partitioning of an element i with radius r_i and valency Z in terms of three parameters: the effective radius of the site r₀(Z), the strain-free partition coefficient D₀(Z) for a cation with radius r₀(Z), and the apparent compressibility of the garnet X-site given by its Young's modulus E_X(Z). Combination of these results with data in Fe-free systems [Van Westrenen, W., Blundy, J.D., Wood, B.J., 1999. Crystal-chemical controls on trace element partitioning between garnet and anhydrous silicate melt. Am. Mineral. 84, 838–847] and crystal structure data for spessartine, andradite, and uvarovite, leads to the following equations for r₀(3+) and E_X(3+) as a function of garnet composition (X) and pressure (P):

r₀(3+) [Å]=0.930X_Py+0.993X_Gr+0.916X_Alm+0.946X_Spes+1.05(X_And+X_Uv)−0.005(P [GPa]−3.0)(±0.005 Å)

E_X(3+) [GPa]=3.5×10¹²(1.38+r₀(3+) [Å])^−26.7(±30 GPa)

Accuracy of these equations is shown by application to the existing garnet–melt partitioning database, covering a wide range of P and T conditions (1.8 GPa<P<5.0 GPa; 975°C<T<1640°C). D^Grt/Melt for all 3+ elements entering the X-site (REE, Sc and Y) are predicted to within 10–40% at given P, T, and X, when D^Grt/Melt for just one of these elements is known. In the absence of such knowledge, relative element fractionation (e.g. D_Sm^Grt/Melt/D_Nd^Grt/Melt) can be predicted. As an example, we predict that during partial melting of garnet peridotite, group A eclogite, and garnet pyroxenite, r₀(3+) for garnets ranges from 0.939±0.005 to 0.953±0.009 Å. These values are consistently smaller than the ionic radius of the heaviest REE, Lu. The above equations quantify the crystal-chemical controls on garnet–melt partitioning for the REE, Y and Sc. As such, they represent a major advance en route to predicting D^Grt/Melt for these elements as a function of P, T and X.     

Interpreting granulite facies events through rare earth element partitioning arrays

The use of rare earth element (REE) partition coefficients is an increasingly common tool in metamorphic studies, linking the growth or modification of accessory mineral geochronometers to the bulk silicate mineral assemblage. The most commonly used mineral pair for the study of high‐grade metamorphic rocks is zircon and garnet. The link from U–Pb ages provided by zircon to the P–T information recorded by garnet can be interpreted in relation to experimental data. The simplistic approach of taking the average REE abundances for zircon and garnet and comparing them directly to experimentally derived partition coefficients is imperfect, in that it cannot represent the complexity of a natural rock system. This study describes a method that uses all the zircon analyses from a sample, and compares them to different garnet compositions in the same rock. Using the most important REE values, it is possible to define zircon–garnet equilibrium using an array rather than an average. The array plot describes partitioning between zircon and garnet using D_Yb and D_Yb/D_Gd as the defining features of the relationship. This approach provides far more sensitivity to mineral reactions and diffusional processes, enabling a more detailed interpretation of metamorphic history of the sample.     

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     

An experimental investigation of the partitioning of REE between garnet and liquid with reference to the role of defect equilibria

The partitioning of samarium and thulium between garnets and melts in the systems Mg₃Al₂-Si₃O₁₂-H₂O and Ca₃Al₂Si₃O₁₂-H₂O has been studied as a function of REE concentration in the garnets at 30 kbar pressure. Synthesis experiments of variable time under constant P, T conditions indicate that garnet initially crystallizes rapidly to produce apparent values of D _Sm (D _Sm=concentration of Sm in garnet/concentration of Sm in liquid) which are too large in the case of pyrope and too small in the case of grossular. As the experiment proceeds, Sm diffuses out of or into the garnet and the equilibrium value of D _Sm is approached. Approximate values of diffusion coefficients for Sm in pyrope garnet obtained by this method are 6 × 10^–13 cm² s^–1 at 1,300 ° C and 2 × 10^–12 cm² s^–1 at 1,500 ° C, and for grossular, 8.3 × 10^–12 cm² s^–1 at 1,200 ° C and 4.6 × 10^–11 cm² s^–1 at 1,300 ° C. The equilibrium values of D _Sm have been reversed by experiments with Sm-free pyrope and Sm-bearing glass, and with Sm-bearing grossular and Sm-free glass.Between 12 ppm and 1,000 ppm Sm in pyrope at 1,300 ° C and between 80 ppm and >2 wt.% Tm in pyrope at 1,500 ° C, partition coefficients are constant and independent of REE concentration. Above 100 ppm of Sm in garnet at 1,500 ° C, partition coefficients are independent of Sm concentration. At lower concentrations, however, D _Sm is dependent upon the Sm content of the garnet. The two regions may be interpreted in terms of charge-balanced substitution of Sm₃Al₅O₁₂ in the garnet at high Sm concentrations and defect equilibria involving cation vacancies at low concentrations. At very low REE concentrations (< 1 ppm Tm in grossular at 1,300 ° C) D_REE garnet/liquid again becomes constant with an apparent Henry's Law value greater than that at high concentrations. This may be interpreted in terms of a large abundance of cation vacancies relative to the number of REE ions.The importance of defects in the low concentration region has been confirmed by adding other REE (at 80 ppm level) to the system Mg₃Al₂Si₃O₁₂-H₂O at low Sm concentrations. These change D _Sm in the defect region, demonstrating their role in the production of vacancies.Experiments on a natural pyropic garnet indicate that defect equilibria are of importance to REE partitioning within the concentration ranges found in nature.     

An experimental calibration of the “Nickel in Garnet” geothermometer with applications

The partitioning of Ni between olivine and garnet in two multi-component compositions was studied as a function of temperature and pressure in a 6–8 type multi-anvil apparatus. Weighted least squares regression of the experimental results for both compositions give a temperature dependence for Ni partitioning between olivine (Ol) and garnet (Gt) of the form: T=–10210(±114)/[lnD _Gt/Ol ^Ni –3.59(±67)] where DGt/Ol/Ni=Ni content of Gt/Ni content of olivine (by weight) and Tis in Kelvin. The pressure effect on Ni exchange between olivine and garnet is estimated to be 0.13 J/bar, and should therefore have a negligible effect on the geothermometer, as is demonstrated in the experiments. The experimentally derived Ni in Gt geothermometer is in agreement with an extant empirical version between 1100 to 900°C, but the two thermometers produce significantly different results outside of this temperature range.     

Hydrogen partitioning between melt,clinopyroxene, and garnet at 3 GPa in a hydrous MORB with 6 wt.% H<Subscript>2</Subscript>O

To understand partitioning of hydrogen between hydrous basaltic and andesitic liquids and coexisting clinopyroxene and garnet, experiments using a mid-ocean ridge basalt (MORB) + 6 wt.% H₂O were conducted at 3 GPa and 1,150–1,325°C. These included both isothermal and controlled cooling rate crystallization experiments, as crystals from the former were too small for ion microprobe (SIMS) analyses. Three runs at lower bulk water content are also reported. H₂O was measured in minerals by SIMS and in glasses by SIMS, Fourier Transform infrared spectroscopy (FTIR), and from oxide totals of electron microprobe (EMP) analyses. At 3 GPa, the liquidus for MORB with 6 wt.% H₂O is between 1,300 and 1,325°C. In the temperature interval investigated, the melt proportion varies from 100 to 45% and the modes of garnet and clinopyroxene are nearly equal. Liquid composition varies from basaltic to andesitic. The crystallization experiments starting from above the liquidus failed to nucleate garnets, but those starting from below the liquidus crystallized both garnet and clinopyroxene. SIMS analyses of glasses with >7 wt.% H₂O yield spuriously low concentrations, perhaps owing to hydrogen degassing in the ultra-high vacuum of the ion microprobe sample chamber. FTIR and EMP analyses show that the glasses have 3.4 to 11.9 wt.% water, whilst SIMS analyses indicate that clinopyroxenes have 1,340–2,330 ppm and garnets have 98–209 ppm H₂O. D _H ^cpx−gt is 11 ± 3, D _H ^cpx−melt is 0.023 ± 0.005 and D _H ^gt−melt is 0.0018 ± 0.0006. Most garnet/melt pairs have low values of D _H ^gt−melt, but D _H ^gt−melt increases with TiO₂ in the garnet. As also found by previous studies, values of D _H ^cpx−melt increase with Al₂O₃ of the crystal. For garnet pyroxenite, estimated values of D _H ^{pyroxenite−melt} decrease from 0.015 at 2.5 GPa to 0.0089 at 5 GPa. Hydration will increase the depth interval between pyroxenite and peridotite solidi for mantle upwelling beneath ridges or oceanic islands. This is partly because the greater pyroxene/olivine ratio in pyroxenite will tend to enhance the H₂O concentration of pyroxenite, assuming that neighboring pyroxenite and peridotite bodies have similar H₂O in their pyroxenes. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

An experimental study of Fe-Mg partitioning between garnet and olivine and its calibration as a geothermometer

The partitioning of Fe and Mg between coexisting garnet and olivine has been studied at 30 kb pressure and temperatures of 900 ° to 1,400 °C. The results of both synthesis and reversal experiments demonstrate that K _D (= (Fe/Mg)_gt/(Fe/Mg)_OI) is strongly dependent on Fe/Mg ratio and on the calcium content of the garnet. For example, at 1,000 °C/30 kb, K _D varies from about 1.2 in very iron-rich compositions to 1.9 at the magnesium end of the series. Increasing the mole fraction of calcium in the garnet from 0 to 0.3 at 1,000 ° C increases K _D in magnesian compositions from 1.9 to about 2.5.The observed temperature and composition dependence of K _D has been formulated into an equation suitable for geothermometry by considering the solid solution properties of the olivine and garnet phases. It was found that, within experimental error, the simplest kind of nonideal solution model (Regular Solution) fits the experimental data adequately. The use of more complex models did not markedly improve the fit to the data, so the model with the least number of variables was adopted.Multiple linear regression of the experimental data (72 points) yielded, for the exchange reaction: 3Fe₂SiO₄+2Mg₃Al₂Si₃O₁₂ olivine garnet 2Fe₂Al₂Si₃O₁₂+3Mg₂SiO₄ garnet olivine H ° (30kb) of –10,750 cal and S ° of –4.26 cal deg^–1 mol^–1. Absolute magnitudes of interaction parameters (W _ij ) derived from the regression are subject to considerable uncertainty. The partition coefficient is, however, strongly dependent on the following differences between solution parameters and these differences are fairly well constrained: W _FeMg ^ol -W _FeMg ^gt 800 cal W _CaMg ^gt -W _CaFe ^gt 2,670 cal.The geothermometer is most sensitive in the temperature and composition regions where K _D is substantially greater than 1. Thus, for example, peridotitic compositions at temperatures less than about 1,300 ° C should yield calculated temperatures within 60 °C of the true value. Iron rich compositions (at any temperature) and magnesian compositions at temperatures well above 1,300 °C could not be expected to yield accurate calculated temperatures.For a fixed K _D the influence of pressure is to raise the calculated temperature by between 3 and 6 °C per kbar.     

Experimental determination of Mn-Mg mixing properties in garnet,olivine and oxide

We have measured the mixing properties of Mn-Mg olivine and Mn-Mg garnet at 1300° C from a combination of interphase partitioning experiments involving these phases, Pt-Mn alloys and Mn-Mg oxide solid solutions. Activity coefficients of Mn dilute in Pt-Mn alloys at 1300° C/1 atm were measured by equilibrating the alloy with MnO at known f _{O 2}. Assuming that the log f _{O 2} of the Mn-MnO equilibrium under these conditions is-17.80 (Robie et al. 1978), we obtain for _Mn: log_Mn = –5.25 + 3.67 X_Mn + 11.41X² _Mn Mixing properties of (Mn,Mg)O were determined by reversing the Mn contents of the alloys in equilibrium with oxide at known f _{O 2}. Additional constraints were obtained by measuring the maximum extent of immiscibility in (Mn,Mg)O at 800 and 750° C. The data are adequately described by an asymmetric (Mn,Mg)O solution with the following upper and lower limits on nonideality: (a) W_Mn = 19.9kj/Mol; W_Mg = 13.7kj/Mol; (b) W_Mn = 19.9kj/Mol; W_Mg = 8.2kj/Mol; Olivine-oxide partitioning was tightly bracketed at 1300° C and oxide properties used to obtain activity-composition relations for Mn-Mg olivine. Despite strong M2 ordering of Mn in olivine, the macroscopic properties are adequately described by a symmetric model with: W^ol = 5.5 ± 2.5 kj/mol (1-site basis) Using these values for olivine, garnet-olivine partitioning at 27 kbar/1300° C leads to an Mn-Mg interaction parameter in garnet given by: W^gt = 1.5 ± 2.5kJ/mol (1-site basis) Garnet-olivine partitioning at 9 kbar/1000° C is consistent with the same extent of garnet nonideality and the apparent absence of excess volume on the pyrope-spessartine join indicates that any pressure-dependence of W^Gt must be small. If olivine and garnet properties are both treated as unknown and the garnet-olivine partitioning data alone used to derive W^Ol and W^Gt, by multiple linear regression, best-fit values of 6.16 and 1.44 kJ/mol. are obtained. These are in excellent agreement with the values derived from metal-oxide, oxide-olivine and olivine-garnet equilibria.     

The influence of low aluminum concentrations on the composition and conditions of crystallization of majorite–knorringite garnets: Experiment at 7.0 GPa and 1500–1700°C

Crystallization of garnet in high-chromium restite formed under the conditions of partial melting in the spinel facies and subsequently subducted into the garnet depth facies was studied experimentally in the MgO–Al₂O₃–Cr₂O₃–SiO₂ system. The crystallization of garnet and the dependence of its composition on the temperature and bulk composition of the system with low Al concentration were studied as well. Experiments in the knorringite–majorite–pyrope system with 5, 10, and 20 mol % Prp were carried out at 7 GPa. The phase associations for the starting composition of pure knorringite Mg₃Cr₂Si₃O₁₂ included chromiumbearing enstatite MgSiO₃ (up to 3.2 wt % Cr₂O₃) and eskolaite Cr₂O₃. Addition of Al resulted in crystallization of high-chromium majoritic garnet. The portion of garnet in the samples always exceeded the concentration of pyrope in the starting composition owing to the formation of the complex majorite–knorringite–pyrope series of solid solutions. With increasing content of pyrope (from 5 to 20 mol %) and increasing temperature, the modal concentration of garnet increased significantly (from 6–12 to 22–37%). The garnet was characterized by high concentrations of the pyrope (23–80 mol %) and knorringite (22–70 mol %) components. The excess of Si (>3 f.u.) with decreasing Cr concentration provided evidence for the contribution of the majorite–knorringite trend to the variation in garnet composition. On the basis of the natural data, most of the garnets composing xenoliths of ultrabasic rocks in kimberlites and occurring as inclusions in diamonds are low-chromium; i.e., their protolith was not subjected to partial melting, at least in the spinel depth facies.     

Partition coefficients for REE between garnets and liquids: implications of non-Henry's Law behaviour for models of basalt origin and evolution

Partition coefficients for the rare earth elements (REE) Ce, Sm and Tm between coexisting garnets and hydrous liquids have been determined at high pressure and temperatures (30 kbar and 1300 and 1500°C). Two synthetic systems were studied, Mg₃Al₂Si₃O₁₂-H₂O and Ca₃Al₂Si₃O₁₂-H₂O, in addition to a natural pyrope-bearing system.Deviations from Henry's Law behaviour occur at geologically relevant REE concentrations. At concentrations < 3 ppm Ce, < 12 ppm Sm, < 80 ppm Tm in pyrope and < 100 ppm Ce, < 250 ppm Sm, < 1000 ppm Tm in grossular (at 30 kbar and 1300°C), D^{garnet liquid}_REE increases as the REE concentration in the garnet decreases. At higher concentrations, D_REE is constant. D^{grossular liquid}_REE also constant when the garnet contains less than about 2 ppm Sm or Tm. The REE concentration at which D_REE becomes constant increases with increasing temperature, decreasing REE ionic radius and increasing Ca content of the garnet.Partitioning behaviour of Ce, Sm and Tm between a natural pyrope-rich garnet and hydrous liquid is analogous to that in the synthetic systems and substantiates the substitution model proposed by Harrison and Wood (1980).Values of D_REE^{garnet/liquid} for which Henry's Law is obeyed are systematically higher for grossular than for pyrope (D^{pyrope/liquid} = 0.067(Ce), 0.108(Sm), 0.155(Tm) and D^{grossular/Liquid} = 0.65(Ce), 0.75(Sm), 4.55(Tm).The implications of non-Henry's Law partitioning of REE for models of basalt petrogenesis involving garnet are far-ranging. Deviations from Henry's Law permit refinements to be made to calculated REE abundances once basic model parameters have been defined.