Natural Trace Element Distribution Coefficients for Garnet, Clino- and Orthopyroxene: Variations with Temperature and Pressure

Van Achterbergh, E.1, Griffin, W.L.1,2, Shee, S.R.3, Wyatt, B.A.3 and Sharma, A.L.1

1. GEMOC National Key Centre, School of Earth Sciences, Macquarie University, NSW 2109, Australia

2. CSIRO Exploration and Mining, P.O. Box 136, North Ryde, NSW 2113, Australia

3. Stockdale Prospecting Ltd., P.O. Box 126, South Yarra, VIC 3141, Australia

Knowledge of the distribution of elements in minerals is an essential tool for modelling anatectic and metasomatic processes in the upper mantle. Using the high sensitivity Laser Ablation ICP-MS microprobe, new data on the partitioning of trace and rare earth elements (REE) between garnet and pyroxene have been obtained This study is based on a suite of peridotite xenoliths from the Wesselton kimberlite pipe and reports preliminary results for the partitioning of these elements between garnet, clinopyroxene (cpx) and orthopyroxene (opx). Petrographic investigation and detailed electron microprobe analysis of the major elements established that the minerals are homogeneous and likely to be in chemical equilibrium.

Compositional variation (for the major elements) within the selected suite is minimal (with one exception), whereas a range in temperature (T) and pressure (P) is represented. This allows an investigation of the effects that these physical parameters have on the distribution of the elements, independent of compositional effects. In accordance with Watson (1985), deviation from Henry's Law is assumed to be negligible. Details of the major element compositions are given in Table 1 and the trace and REE distribution coefficients are listed in Table 2.

Table 1:Major element data for the Wesselton xenoliths. All values in weight%, ND = not detected. T from Griffin et al. (1989) and Ryan et al., (1996), P from Finnerty and Boyd (1984).

Sample SiO2 TiO2 Al2O3 Cr2O3 FeO MnO MgO CaO Na2O K2O T (oC) P (kbar)
959 gt 40.3 0.03 20.9 3.73 8.34 0.47 19.4 5.25 0.03 ND 872 36
960 gt 39.8 0.05 20.6 3.79 8.34 0.55 19.9 5.01 0.04 ND 809 32
965 gt 39.9 0.07 21.1 3.24 8.26 0.44 19.9 4.90 0.03 ND 819 33
966 gt 40.0 0.03 18.6 6.69 7.18 0.36 19.7 5.89 0.03 ND 997 46
968 gt 40.6 0.04 20.6 3.89 7.30 0.42 20.2 5.24 0.03 ND 1013 46
959 cpx 53.8 0.09 2.23 1.83 2.15 0.07 15.6 21.0 1.93 0.01 872 36
960 cpx 53.8 0.13 3.03 2.17 2.50 0.02 14.6 19.8 2.65 0.01 809 32
965 cpx 54.4 0.14 2.48 2.07 2.29 0.06 15.3 20.6 2.11 ND 819 33
966 cpx 53.8 0.03 2.23 2.44 2.28 0.07 16.0 19.4 2.26 0.01 997 46
968 cpx 54.0 0.02 1.73 1.40 2.26 0.10 17.0 20.9 1.44 0.02 1013 46
994 cpx 53.1 0.04 2.51 2.18 2.25 0.05 15.0 20.6 2.25 0.01 821 34
995 cpx 52.6 0.17 2.36 1.66 2.43 0.03 15.4 20.8 2.02 ND 851 36
996 cpx 53.3 0.06 2.24 1.86 2.25 0.04 15.4 21.0 1.91 0.01 862 36
1001 cpx 52.7 0.45 2.28 2.41 3.71 0.10 16.8 17.7 1.96 0.01 1080 51
1007 cpx 54.0 0.25 2.31 2.40 2.38 0.07 16.2 18.5 2.29 0.05 1047 49
994 opx 57.0 0.02 0.72 0.26 4.94 0.12 36.6 0.25 0.06 0.002 821 34
995 opx 57.5 0.06 0.72 0.22 5.27 0.12 36.1 0.26 0.07 0.004 851 36
996 opx 57.5 0.03 0.74 0.26 5.05 0.11 36.3 0.27 0.04 0.002 862 36
1001 opx 56.9 0.06 0.69 0.42 6.45 0.12 34.6 0.62 0.16 0.006 1080 51
1007 opx 57.7 0.06 0.70 0.37 4.39 0.14 36.1 0.58 0.17 ND 1047 49

 

Garnet - Clinopyroxene Partitioning: A comparison of the new distribution coefficient data from Wesselton with data available in the literature (and some unpublished data), shows generally good agreement (Figure 1). This is despite the fact that the trace element composition of both the garnet and cpx reflects extensive metasomatism, which is commonly accompanied by the introduction of phlogopite. This suggests that regardless of the secondary processes involved, equilibration had been achieved before entrainment in the kimberlite.

The positive Nb and Zr anomalies in the otherwise smooth curves of the distribution coefficient plot are consistent with the literature data and with the experimental findings of Green et al. (1989) that the high field strength elements (HFSE) tend to concentrate in garnet. The distribution is controlled dominantly by the crystal chemistry, and the better fit of these ions into the garnet structure. As the temperature increases, thermal expansion of the lattice allows more Zr to enter garnet and DZr (garnet/cpx) therefore increases with increasing temperature. In accordance with the results of O'Reilly and Griffin (1995), no pressure effect was observed on the distribution of these elements between cpx and garnet. The distribution of Ti also appears to be temperature dependent. However, data from peridotite xenoliths in basalts from eastern China (Xu et al., 1998), suggest that pressure plays a more important role in the partitioning of Ti. The calculated pressures for these samples are lower than those measured for the Wesselton suite (19-20 kbar vs 32-46 kbar), while the temperature is higher (1117oC-1146oC vs 809oC-1013oC). The differences between the Chinese and South African samples indicate that at lower pressure, Ti prefers the cpx structure (DTi (garnet/cpx) <1), and at higher pressure the garnet structure (DTi (garnet/cpx) >1). Although the two suites are broadly similar in composition (calcic lherzolitic garnet and diopside), the Chinese cpx contain slightly less Ca and the garnet less Cr, which may have some effect on this trend.

Higher temperature and pressure also increase the D (garnet/cpx) for the light-(L) and medium(M) REE, but have no observable effect on the heavy REE. Therefore, the slope of the distribution coefficient curve changes, becoming shallower at higher temperature. This is consistent with the data of Mazzucchelli et al. (1992) on crustal garnet-cpx pairs, which are at the lowest temperature (around 800 oC) and show the steepest slope, and with the flatter slope of Hauri et al.'s (1994) experimental data at 1430 oC (Figure1).

Table 2:Trace and rare earth element distribution coefficients for garnet-cpx and cpx-opx pairs
 

Element gt/cpx

959

gt/cpx

960

gt/cpx

965

gt/cpx

966

gt/cpx

968

cpx/opx

994

cpx/opx

995

cpx/opx

996

cpx/opx

1001

cpx/opx

1007

Sr 0.0007 0.0007 0.0004 0.0023 0.0013 305 156 354 161 392
Nb 0.564 0.284 0.190 0.559 1.15          
La 0.008 0.007 0.003 0.022 0.013 151 149 186 59.8 13.3
Ce 0.009 0.008 0.006 0.034 0.030 129 218 456 50.3  
Pr 0.032 0.016 0.022 0.078 0.081          
Nd 0.064 0.052 0.047 0.191 0.172          
Zr 0.760 0.540 0.563 2.28 3.89 102 140 114 46.3 75.6
Hf 0.230 0.153 0.165 1.00 1.54          
Sm 0.344 0.319 0.273 0.832 0.722 86.8 12.3   7.87  
Eu 0.486 0.643 0.597 1.00 0.955          
Ti 0.947 0.640 0.524 1.56 7.93 2.25 2.39 2.15 12.6 2.74
Gd 1.01 0.952 1.39 1.42 1.63 11.5 6.00   9.37  
Dy 3.51 4.21 6.11 4.57 2.48          
Ho 6.50 6.72 10.2   1.16          
Y 9.25 9.29 9.78 8.84 8.34 75.3 55.1   52.0 30.5
Er 6.70 8.56 14.5   1.71 0.919     1.63  
Yb   9.8                
Lu 3.47 4.87     3.33          
Sc 4.10 3.23 3.67 5.17 6.89 9.47 10.3 10.0 6.47 7.67
Ga 1.27 0.982 1.09 1.89 2.27 1.31 1.51 1.52 2.14 1.58
V 0.411 0.318 0.350 0.712 0.875 8.38 8.35 7.88 5.72 6.26
Co 2.78 2.66 3.14 2.25 2.21 0.242 0.240 0.251 0.333 0.309
Ni 0.090 0.083 0.094 0.125 0.120 0.332 0.332 0.351 0.429 0.373
CaO 0.250 0.252 0.237 0.304 0.250 52.2 60.4 70.8 22.7 28.2

Figure 1: Garnet-cpx distribution coefficients for Wesselton xenoliths and comparative literature data

Clinopyroxene - Orthopyroxene Partitioning: Published data on cpx-opx partitioning are scarce, and a comprehensive comparison is not possible. However, where data are available (eg. McDonough et al., 1992) good agreement exists between datasets. The LREE abundances in opx are very low, and the data are characterised by large standard deviations. Of the elements analysed (Table 1), only Co and Ni prefer the opx structure. Ga and Er partition approximately equally and all other elements fit more easily into cpx. The well-documented T dependent partitioning of Ca into opx (Brey and Köhler, 1990) is evident in these data: as T increases, so the Ca in opx increases. The LREE, Zr, Sc and V follow the same trend. A single sample (1001) has lower XMg (0.81) than the rest of the suite (0.84-0.86) and as Mg decreases, more LREE, and fewer HREE enter the crystal structure. The Ti distribution coefficient also increases markedly as XMgopx decreases.

References Cited

Brey, G.P. and Köhler, T. (1990). J. Petrol., 31, 1353-1378

Finnerty, A.A. and Boyd, F.R. (1984). Geochim. Cosmochim. Acta, 48, 15-27

Green, T.H., Sie, S.H., Ryan, C.G. and Cousens, D.R. (1989). Chemical Geology, 74, 201-216

Griffin,W.L., Cousens,D.R., Ryan,C.G., Sie,S.H. and Suter,G.F.(1989). Contrib. Mineral. Petrol. 103, 199-202

Harte, B. and Kirkley, M.B. (1997). Chemical Geology, 136, 1-24

Hauri, E.H., Wagner, T.P. and Grove, T.L. (1994). Chemical Geology, 117, 149-166

Mazzucchelli, M., Rivalenti, G., Vannucci, R., Bottazzi, P., Ottolini, L., Hofmann, A.W., Sinigoi, S. and Demarchi, G. (1992). Geochim. Cosmochim. Acta, 56, 2371-2385

McDonough, W.F., Stosch, H.G. and Ware, N.G. (1992). Contrib. Mineral. Petrol., 110, 321-328

O'Reilly, S.Y. and Griffin, W.L. (1995). Chemical Geology, 121, 105-130

Ryan, C.G., Griffin, W.L. and Pearson, N.J. (1996). J. Geophys. Res. 101, 5611-5625.

Shimizu, N. and Allègre, C.J. (1978). Contrib. Mineral. Petrol. 67, 41-50

Watson, E.B. (1985). Geochim. Cosmochim. Acta, 49, 917-923

Xu, X., O'Reilly, S.Y., Griffin, W.L., Zhou, X and Huang, X. (1997).In: Mantle Dynamics and Plate Interactions in East Asia (Flower, M., Chung, S.L., Lo, C.H. and Lee, T.Y., eds). Amer. Geophys. Union Spec. Publ., in press