Trace-element partitioning between garnet and clinopyroxene in mantle-derived pyroxenites and eclogites: P-T-X controls

CHEMICAL GEOLOGY INCLUDING

ISOTOPE GEOSCIENCE

ELSEVIER

Chemical Geology 121 (1995) 105-130

Trace-element partitioning between garnet and clinopyroxene in mantle-derived pyroxenites and eclogites: P - T - X controls Suzanne Y. O'Reilly

W.L. Griffin b

aSchool of Earth Sciences, Macquarie University, Sydney, N.S. W. 2109, Australia bCSIRO Division of Exploration and Mining, Box 136, North Ryde, Sydney, N.S. W. 2113, Australia

Received 26 May 1994; revision accepted 25 October 1994

Abstract Trace-element abundances in coexisting clinopyroxenes and garnets from five suites of mantle-derived garnet pyroxenite and eclogite xenoliths and two suites of eclogite inclusion pairs in diamonds have been determined by electron microprobe (Ti) and proton microprobe (Ni, Zn, Ga, Sr, Y, Zr). The sample sets provide garnet--clinopyroxene pairs from a range of depths on contrasting geothermal gradients (40-90 mW m 2) from South Africa, and western and eastern Australia. These data are used to assess the effects of phase composition, pressure (P) and temperature (T) on the partitioning of each element between garnet and clinopyroxene. Ni partitioning is moderately dependent on XMg and strongly dependent on T. Dzn (Zncpx/Znc,,t) is ~ 0.9 + 0.1 and may show a slight increase with T, but is independent of P. Ga partitioning is strongly dependent on X cpx as well as P and T. Sr and Y distributions are mainly controlled"oy aca " c,nt and are weakly dependent on Tbut not on P. Zr partitioning is mainly dependent on T; Dzr decreases with increasing X c°x and posslbly ' Gnt, but shows no P effect. DT~ decreases markedly with increasing P (and Xca possibly T) and has a weak tendency to increase with increasing X cox. These results on trace-element partitioning in natural, equilibrated systems provide data for the modelling of melting and metasomatic processes in the upper mantle. They also provide a test dataset to evaluate experimental determinations of partition coefficients, and some relationships (e.g., Dzr, DNI) may be used to cross-check thermobarometry calculations where spurious calculated Fe 3+/Fe 2÷ ratios may give invalid temperature estimates.

1. I n t r o d u c t i o n Element partitioning between coexisting mineral pairs provides data fundamental to crystal chemistry, geothermobarometry and the mathematical modelling of geological processes such as partial melting, fractional crystallization and metasomatism. Considerable effort therefore has been expended on the experimental measurement of major- and minor-element partitioning as a function o f pressure ( P ) and temperature (T) [EW I 0009-2541/95/$09.50 © 1995 Elsevier Science B.V. All rights reserved SSD10009-254 1 (94)00 147-2

under mantle conditions. However, these studies commonly have problems with the attainment of equilibrium, and the necessity of doping the charge to reach analytical detection limits leads to the use of unnatural compositions. Furthermore, few experiments on minorelement partitioning have assessed the effect of compositional variability in the major-element composition of the phases involved. Studies of element distribution in mantle-derived xenoliths can provide critical data to support experimental studies, although they also must face the problem of assessing the degree of equilibrium between

106

s. E O'Red/y,W.L. Griffin/ Chemical Geology 121 (1995) 105-130

phases. The natural samples have other limitations; many xenolith suites represent samples taken along a single geotherm, where P and T vary sympathetically in some regular fashion. The natural processes of sampling in any one xenolith suite therefore provide a relatively narrow range of well-correlated P - T whereas in the laboratory, P and Tcan be varied independently, and P - T combinations unlikely to occur in nature are routinely used to study P-Teffects. However, selection of xenolith samples from areas with different geothermal regimes can provide both a useful range of P - T combinations, and many xenoliths appear to be derived from depths beyond the reach of routine experimental studies. Finally, a wider range of complex compositions is available in xenoliths, while only a small range can be studied in a time-consuming laboratory research program. Analysis of element distribution in such xenolith assemblages therefore can provide valuable data which can guide experimental work and which in turn can be used to evaluate the validity of experimental results. In this study, we have determined trace-element abundances by electron microprobe (Ti) and proton microprobe (Ni, Zn, Ga, Sr, Y, Zr) in coexisting clinopyroxenes and garnets from mantle-derived garnet pyroxenite (i.e. low-Jd Cpx) and eclogite (i.e. Jd-rich Cpx) xenoliths and from inclusion pairs in diamonds. All of these minerals are low in Cr; partitioning between Cr-rich garnet, orthopyroxene and chrome diopside in garnet peridotite xenoliths will be discussed elsewhere (W.L. Griffin et al., in prep.). The live suites of xenoliths and two diamond-inclusion suites used here represent a variety of mantle pressure (P) and temperature (T) regimes from a variety of geothermal profiles and therefore can be used to assess the effects of P, T and mineral composition on trace-element partitioning. The samples have been selected to be as free of alteration and disequilibrium microstructures (e.g., zoning) as possible. The aims of this work are: (1) to establish the concentration ranges of selected trace elements in mantlederived clinopyroxene+low-Cr garnet pairs; (2) to assess as far as possible from these natural samples the effects of pressure, temperature and composition on clinopyroxene/garnet partition coefficients; (3) to provide realistic constraints on partition coefficients for use in modelling of melting and metasomatic processes in the upper mantle; and (4) to establish a framework

for evaluation of experimental determinations of partition coefficients.

2. Samples Xenoliths in kimberlites and basalts are (in principle) useful for partition-coefficient studies because they commonly record the ambient T at the depth of entrainment, quenched in by rapid ascent (Boyd, 1976; O'Reilly, 1989) and eruption. Garnet-clinopyroxene pairs encapsulated in diamonds (and not in contact) record T at the time of their entrapment in the growing diamond; this may be significantly different from the ambient mantle temperature at the time of entrainment in the kimberlitic magma (Griffin et al., 1992, 1993). The sample suites used here were chosen to provide a wide range of P and Tconditions (from single pipes as well as collectively), and include samples that come from different geotherms, i.e. different combinations of P - T conditions (Fig. 1). 2.1. Xenoliths from South African kimberlites

These samples are xenoliths from the kimberlites at Roberts Victor [ including some described by Viljoen et al. (1994)], Newlands and Bobbejahn mines T:=C 200

10

600

O~,~,~.

1000

~.

1400

1800

J grit pyroxenites 5O

2O 3O

O0

40 Q,,

i50 ~. 50

8

60 70 _ _

' ~ ~ diamond Argyleeclogite- sure _ ~k,~lncluslons ? Inclusions In d i a m o ~ 250

80 90

Fig. 1. P and T fields represented by the the xenolithsuites studied. References are given in the text.

S.Y. O'Reilly, W.L. Griffin/ Chemical Geology 121 (1995) 105-130 (grouped below as " R o v i c " ) , Kaalvallei (K.S. Viljoen, 1994 and in prep.) and Monastery Mine (Moore et al., 1991). Several extreme compositions, including grospydites and kyanite- or corundum-bearing eclogites, are included in these suites. Diamond-bearing eclogites come from the Roberts Victor, Orapa and Excelsior mines (Robinson et al., 1984; Hatton and Gurney, 1987; McCandless and Gurney, 1989). The area in which many of these kimberlites erupted had a well-constrained xenolith-based geotherm, corresponding approximately to Pollack and Chapman's (1964) model of conductive geotherm for 40-mWm -2 heat flow (Finnerty and Boyd, 1987) to which clinopyroxene/garnet temperatures may be referred to obtain approximate pressure estimates (Table 1). Diamond eclogites are stable only at P > 40 kbar on this geotherm. 2.2. South African and Australian diamond inclusions These samples represent clinopyroxene-garnet pairs, not in contact, from single diamonds. The diamond inclusions (hereafter referred to as DI) from Argyle diamonds were described by Griffin et al. (1988) and the data on the South African DI are given by Gurney et al. ( 1979, 1986a, b), Moore and Gurney (1989), Rickard et al. (1989) and J.J. Gurney (unpublished data). Both areas had a 40-mW-m -2 geotherm at the time of eruption (Finnerty and Boyd, 1987; Jaques et al., 1990), but the high temperatures recorded by some DI pairs may reflect growth of the diamonds during essentially isobaric thermal pulses, rather than derivation from extreme depths (Griffin et al., 1992; Fig. 1). 2.3. Garnet clinopyroxenites from Bullenmerri and Gnotuk Maars (B / G), Victoria, Australia The petrology of these xenoliths (hereafter called the B / G suite) has been described by Griffin et al. (1984), and the elevated geotherm defined by these xenoliths is discussed by O'Reilly and Griffin (1985). The importance of this suite is that the xenoliths reflect significantly lower P at any given T than either of the other suites. Thus a T of 1000°C will be reached at ~ 40-kbar pressure in the South African xenolith suite, but at ~ 12 kbar in the B / G suite (Fig. l ). Comparison of these suites therefore provides a means of separating

107

the effects of pressure and temperature on element partitioning. Trace-element abundances and residence sites for mantle-derived spinel lherzolite xenoliths from this locality were reported previously (O'Reilly et al., 1991). Two xenoliths from another locality in eastern Australia are also included within this group as they lie on the same geotherm and extend the P range (samples 69-27 and 38926 from the Delegate locality; Griffin and O'Reilly, 1986).

3. Problems Problems encountered in interpreting the trace-element data include the usual decisions on appropriate thermobarometry calculations (Finnerty and Boyd, 1987) and the possible effects of disequilibrium. We have used the Ellis and Green (1979) garnetclinopyroxene thermometer to estimate temperatures for each mineral pair. This thermometer, and its various modifications (Krogh, 1988) are extremely sensitive to the ferric/ferrous iron ratio in the clinopyroxene, which is in turn sensitive to analytical error, especially in Si (Neumann, 1976). The Fe 3÷/Fe 2 ÷ ratio has been estimated for each analysis by normalization to four cations. In most of the xenoliths used here, the calculated Fe 3 ÷ content of both clinopyroxene and garnet is within the analytical errors; Fe 3 + therefore is taken as zero and the T estimates are accepted as reliable. In some cases the calculated Fe 3÷ contents of the pyroxenes are significant, and lead to improbably low temperatures, as noted by numerous earlier workers. Re-analysis does not remove this effect, and we interpret at least part of the high Fe 3 + contents as the result of late-stage (deuteric?) oxidation of the pyroxenes. In these cases (underlined in Table 1) we have arbitrarily taken an average value between the temperatures calculated assuming zero and maximum Fe 3 + contents; this procedure typically changes the estimated T by 100-150°C and may produce more scatter in several of the plots shown below. However, removal of these samples does not produce much improvement in the regressions, and the Tuncertainty does not significantly affect the deliberately broad conclusions drawn from the data. An approximate pressure has been estimated for each sample (Table 1) by referring T to the appropriate

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

108

Table 1 Representative data for clinopyroxenes and garnets Eastern Australia pyroxenites Sample Locality Comment

Ti/Gnt (ppm) Ti/Cpx Cr/Gnt Cr/Cpx Ni/Gnt Ni/Cpx Zn/Gnt Zn/Cpx Ga/Gnt Ga/Cpx Sr/Gnt Sr/Cpx Y/Gnt Y/Cpx Zr/Gnt Zr/Cpx D, Ti D, Cr D, Ni D, Zn D, Ga D, Sr D, Y D, Zr X~'~ X~rg"

X~2' %Jd, Cpx %K20, Cpx %Na20, Gnt D(con'), Ni D(COIT), Ga D(corr), Sr D (corr), Zr T (°C) P (kbar)

Diamond eclogites

GNSI GNWI D R I O I 6 2 BMI71 Lakes Bullenmerri/Gnotuk, western Victoria

GN47

580 465 3515 2640 1200 998 1240 684 46_+10 35_+2 332+6 415_+6 32_+1 37_+1 31-+1 39_+1 8-+0.5 6_+0.5 18_+1 15_+0.6 <1 ~ <1 155_+3 83_+2 24_+1 20_+1 3_+0.5 2_+0.5 12_+0.5 4.5-+0.5 20_+1 8-+1

b.d. b 1620 1850 1230 18±2 298_+3 4_+1 3_+1 4_+0.5 3_+0.5 <1 135-+2 12_+1 4_+0.3 5-+1 8+1

6.1 1.0 7.2 0.97 2.3 #

5.7 0.7 11.9 1.05 2.5 #

0.13 1.7 0.623 0.802 0.20 2.5 b.d. b.d. 14.0 2.2 # 1.7 1090 16

0.10 1.6 0.675 0.823 0.13 11 b.d. b.d. 17.1 2.3 # 1.8 1045 15

330 2260 958 1440 49_+2 571 _+5 24±0.5 26-+1 7±0.4 14-+0.5 1_+0.3 153-+2 8+_0.4 1_+0.3 5_+0.3 7_+0.8

805 4715 2940 2530 46_+2 565_+6 36_+1 38_+1 7_+0.5 15_+1 1_+0.5 51-+1 23+1 4_+0.5 6+_1 13_+1

6.9 1.5 11.7 1.08 2.0 153 0.13 1.4

5.9 0.9 12.3 1.06 2.1 51 0.17 2.2

0.702 0.843 0,15 4.6 b.d. b.d. 17.8 1.9 165 1.5 1065 15

0.678 0.81 0.20 6.4 b.d. b.d. 17.5 2.0 # 2.3 1085 16

#~ 0.7 16.6 0.75 0.8 # 0.33 1.6 0.707 0.868 0.14 0.5 b.d. b.d. 23.1 0.7 # 1.6 955 12

G N 6 4 9 ( 1 ) 69-27

XMI I

JJG892

Delegate, N.S.W. Sapph ~

Orapa graphite

Orapa

940 3215 6770 4720 44_+3 367_+5 21_+1 21+_1 6_+1 10_+1 2+0.5 20_+0.5 68_+2 7+_0.5 19_+1 13_+1

390 1460 440 684 65_+3 308_+4 24_+1 16+_1 4_+0.5 10_+1 2_+0.5 108_+2 33_+1 5-+1 32_+2 17+1

900 2400 1710 2460 62+_2 315_+6 9_+1 27_+2 9.5±0.5

3780 3900 616 342 31 +3 179_+3 75+2 65+ 1 12.5+0.6 19_+1 1.4_+0.4 166_+2 22-+1 2.5_+0.4 21-+1 15_+ I

5.9 1.4 13.0 1.00 0.2 24 0.17 1.4

3.4 0.7 8.3 1.00 1.7 10 0.10 0.7

3.7 1.6 4.7 0.67 2.5 54 0.15 0.5

2.7 1.4 5.1 3.00

0.673 0.857 0.14 7.2 b.d. b.d. 18.7 0.0 38 1.5 955 12

0.758 0.862 0.15 1.9 b.d. b.d. 14.6 1.6

0.821 0.88 0.21 17 b.d. b.d. 10.3 2.2 75 0.8 1340 50

420 2480 958 1370 26_+2 338_+4 7_+1 7_+1 4.5-+0.4 7_+1 3_+0.5 71-+1 18+1 3±0.3 9.5_+0.5 14_+1

# 0.7 I 100 17

38926

16_+1 <2 340_+3 17-+1 3_+1 21_+3 65-+3

1.0 0.6 5.8 0.87

1.7

1.5

0.18 3.1

119 0.11 1).7

#

0.88 0.921 0.10 17 b.d 0.16 I 1.5 1.3 # 3.3 1000 45

0.672 0.823 0.17 25 0.06 0.16 8.0 1.0 136 1.1 1245 67

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

Diamond eclogites (cont.)

109

African diamond inclusions

JJG889

AKI-IOc

H R V 2 3 4 7 b HRV247F

CBS17327

JJG144

Orapa

Orapa

Roberts Victor graphite

Roberts Victor graphite

Newlands

Newlands

PREM45 Premier

1620 1440 411 479 48±3 486+_6 39+_1 38+_1 15+_0.5 21+_1 1.4+_0.5 353+_3 11+_0.6 <1.8 32+_ I 33+_2

2280 1860 342 479 59±2 499+_5 48+_2 38±1 13.5+_0.5 22+_0.5 1.2±0.4 288+_3 10+_0.7
2400 2400 547 411 82±2 611 +_6 36+_I 37+_I 13.5+_0.5 16.5+_0.5 <1.3 255+_3 14+_0.7 2.5+_0.5 21 ± 1 31 ± 2

1740 2340 274 479 26±3 217±2 104+_3 90+_2 13+_0.5 25+_0.5 <1.3 198_+4 31+_I 2.2+_0.4 16.5± 1 22.5 ± 1

4740 3540 2580 3540 3360 4860 3240 2880 3720 3600 2950 2730 3220 2260 1290 700+270 a 810+ 150 1060± 170 955 ± 175 994_+200 71±5 50+_4 18+_3 16+_3 16+_4 242+_8 168+_5 84+_5 84+_5 100+_4 73+_2 99±3 103+_3 148±4 112±3 84+_3 85±3 89±2 131_+3 110_+3 16±1.5 1 1 . 5 + _ 1 1 4 + _ 1 . 5 19+_1.5 16+_I 21+_2 21+_I 15±1 27_+2 21_+1 4±1 12±1 10.5+_I 7.5+_1.5 7±1 86±3 139±4 940+_24 392_+6 278+_4 47±3 51±3 25+_1.5 57+_4 47+_2 4.5+_I 4+_1 4±t 3.5+1 4±1 33 ± 2 39+_2 28+_3 33 ± 2 27+2 9±2 15+_3 18+_5 32_+3 32+_3

2280 1320 2820 600 70 753 470 684 20__.3 104±3 79+_8 377+_3 78+_3 35+_2 67±5 27+_1 12+_0.5 12.5+_0.5 22+_2 14.5+_0.5 <1.5 3.3+_0.5 305+_6 202+_4 30+_1 8.5+_0.6 2_+0.5 <1.3 20+_2 8± 1 27+_2 3.5+_ 1 1.2 6.8 4.0 0.86 1.8 # 0.10 1.4 0.53 0.8 0.19 38 0.30 0.10 8.0 I. 1 # 1.9 1055 48

0.5 0.9 3.6 0.77 1.2 61 # 0.4 0.79 0.912 0.33 29 0.11 0.12 6.7 0.6 94 0.8 1305 71

0.9 1.2 10.1 0.97 1.4 252 # 1.0 0.741 0.913 0.19 17 0.05 0.09 16.3 h I

271 1.3 1030 47

0.8 1.4 8.5 0.79 1.6 2~ # 0.9 0.677 0.896 0.26 23 0.06 0.11 14.4 1.2 266 1.2 1050 48

h0 0.8 7.5 1.03 1.2 #

1.3 1.8 8.4 0.87 1.9 #

0.18 1.5 0.699 0.832 0.09 25 0.02 0.11 12.1 0.7 # 1.8 1175 57

0.07 1.4 0.517 0.788 0.19 39 0.17 b.d. 12.2 1. I # 1.9 1055 48

PREMI05 ORAPA4

JWA20

Premier

Jwaneng

Orapa

JWA24 Jwaneng

1.0 0.2 3.4 1.15 1.3 22 0.10 0.3

0.9 0.3 3.4 0.86 1.8 12 0.08 0.4

I.! 0.3 4.7 0.86 I.I 90 0.16 0.6

1.1 0.4 5.3 0.89 1.4 52 0.06 1.0

1.1 0.8 6.3 0.98 1.3 40 0.09 1.2

0.544 0.707 0.23 27 0.18 0.20 4.5 08 45 0.6 1360 75

0.505 0.71 0.23 43 0.19 0.25 4.5 h0 35 1.0 1265 68

0.603 0.808 0.19 24 1.12 0.19 8.9 0.6 109 1.0 l 150 56

0.388 0.645 0.26 46 0.06 0.21 6.0 0.5 # 1.6 1060 49

0.447 0.695 0.23 39 0.07 0.19 8.2 0.6 # 1.7 1090 51

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

110

Table 1 (continued) Rovic eclogites Sample Locality Comment T i / G n t (ppm) Ti/Cpx Cr/Gnt Cr/Cpx Ni/Gnt Ni/Cpx Zn/Gnt Zn/Cpx Ga/Gnt Ga/Cpx Sr/Gnt Sr/Cpx Y/Gnt Y/Cpx Zr/Gnt Zr/Cpx D, Ti D, Cr D, Ni D, Zn D, Ga D, Sr D,Y D, Zr Xr~ X~ X~)2' %Jd, Cpx %K20, Cpx %Na20, Gnt D(con-), Ni D(con'), Ga D(corr), Sr D(corr), Zr T (°C) P (kbar)

RV-A Rovic Type 1/Ky

RVIG Rovic Type l

RV2G Rovic Type l

RV3G Rovic Type I

BD3699

KA64-6

74506

BDI I75

BDI I88

Rovic Type 1

Rovic Type 1

Bobbejahn Type 1

Rovic Type II

Rovic Type II

2220 1920 820 0 41 +_2 239+_12 71 _+3 55 _+2 12!1 46+_2 1.8 +_0.8 99 + 4 22_+2 <2 28+_2 14_+3

1620 2340 410 750 33_+2 179+_5 70_+3 60 _+2 10_+1 17_+1 4.2 _+0.6 207 + 3 22_+2 4_+ 1 21 + 2 27_+2

1560 1680 610 820 45+_2 295_+6 78+3 66 _+ 1 11_+1 18+_1 4.1 _+0.7 236 + 3 15_+1.5 <2 21 _+2 33+2

1920 2040 1090 1020 42 298_+5 60_+3 48 _+ 1 10_+1 18_+0.6 3.6 _+0.6 267 _+ 3 17_+2 < 1.5 29_+2 30_+ I

2100 2220 750 815 85_+4 468+7 117_+6 73 + 2 17_+2 31_+1 5 _+ I 223 + 3 21+2 <3 84_+6 50_+3

1680 2160 410 680 47_+2 321 + 2 2 83+2 36 _+2 11.3_+0.6 14.5_+1.5 2.5 _+0.5 200 _+5 33_+2 <2.5 83_+4 18+2

2160 3660 610 950 59_+6 351 _+22 29_+ 1 23 _+ 1 10+1 19.5+_0.7 3.8 +_0.5 250 _+4 27_+2 <2 31 _+2 52_+2

540 1380 1710 1916 75_+2 639+_5 77_+3 53 _+ 1 12_+1 23_+1 7.5 _+0.5 464 _+5 10.5+0.5 <2 20_+ I 28+2

240 300 2740 1370 66_+3 613+6 47_+2 30 + I 11 14.5_+1 2 _+0.5 17 + 1 12_+1 1 _+0.5 5 +0.5 10+0.5

0.9

0.5

1.4 1.8 5.4 0.86 1.7 49 0.18 1.3

0.604 0.73 0.14 69 0.10 0.15 8.4 2.4 69 1.5 1175 57

0.568 0.746 0.35 35 0.12 0.12 8.3 1.0 84 1.8 1210 60

# 5.8 0.77 3.8 55 #

1.1 1.3 6.6 0.85 1.6 58 #

1.1 0.9 7.1 0.80 1.8 74 #

1.6 0.640 0.819 0.12 34 0.13 0.1 I I 1.0 1.0 70 2.0 1175 57

1.0 0.641 0.829 0.16 34 0.17 0.12 11. l 1.1 90 ].5 1150 56

1.1 1.1 5.5 0.62 1.8 45 # 0.6 0.563 0.76 0.30 55 0.30 0.19 8.7 0.7 75 1.4 1100 52

1.3 1.7 6.8 0.43 1.3 80 #

1.7 1.6 6.0 0.79 2.0 66 #

0.2 0.638 0.851 0.15 35 0.15 0.11 11.8 0.6 95 0.7 1175 57

2.6 1.1 8.5 0.69 1.9 62 #

1.7 0.688 0.829 0.12 23 0.02 0.12 10.5 1.5 78 2.0 1060 ~ 49

1.4 0.706 0.895 1/.20 33 0.07 0.07 -3.1 1.3 82 1.8 1040 47

Abbreviations: Sapph = sapphirine-bearing; Cor = corundum-bearing; Ky = kyanite-bearing; San = sanidi ne-bearing; Coes = coesite-bearing. "Sapphirine-bearing pyroxenite (see Griffin and O'Reilly, 1986). bb.d. = below detection limit. " < defines minimum detection limit at 99% confidence limit. ' % gives one standard deviation. *#denotes that one element of the pair is below the minimum detection limit. ~Underlined Tcalculated using Fe 3+ (see text).

1.3 0.5 9.3 I).64 1.3 9 0.08 2.0 0.706 0.928 0.26 29 0.02 0.05 2.7 I).6 35 2.4 920 40

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

111

Rocvia eclogites (cont.)

BD1186

BDII91

74507

7,.1510

Rovic Type 11

Rovic Type 1I

Bobbejahn Type I1

Bobbejahn Type 11

NL-3 Newlands Type II

NL-4 Newlands Type II

Bultfontein Type II

960 I 140 205 205 23_+3 332_+5 74_+3 62_+1 12_+1 24_+1 3.7_+1 19_+ 1 41_+2 <2 9_+1 12_+1

960 I 140 545 340 22_+2 292+10 85_+3 69_+3 12_+1 18_+2 2_+1 18_+2 66_+4 <2 5+I 8_+2

1860 2040 20940 12510 43_+3 284_+5 19_+1 16_+1 15_+1 10.5_+1 9.5_+1 369_+4 31.5_+t.5 6.7_+1.5 44_+2 57_+4

1440 1440 1360 2175 56_+3 390_+5 17_+1 14_+1 14_+1 10_+1 <2.5 174_+3 13_+1 <2.5 31_+2 23_+2

960 1320 2720 2450 69_+3 622_+14 42_+2 38_+2 11_+1 12_+1 <2.5 340_+5 14_+1 <3 19_+1.5 40_+3

1020 900 70 b.d. 38_+4 278_+6 57_+3 40_+2 14.5_+1 16_+1 4.7_+0.7 170_+3 9_+1 <3 12_+1 7.5_+2

420 1500 270 130 20_+5 128_+7 52_+5 65_+3 13_+2 22_+1 <1 61 _+2 31_+3 <2 <4 9_+1

1.2 1,0 14.4 0.84 2.0 5

1.2 0.6 13.3 0.81 1.5 9

#

#

1.3 0.48 0.887 0.28 37 0.01 0.08 20.1 1.3 33 1.8 980 42

1.6 0.469 0.836 0.24 36 0.02 0.09 18.0 0.8 33 2.1 1015 46

1.1 0.6 6.6 0.84 0.7 39 0.21 1.3 0.838 0.93 0.1 I 14 0.05 0.07 13.2 0.4 50 1.5 1025 45

1.0 1.6 7.0 0.82 0.7 # #

1.4 0.9 9.0 0.90 I.I # #

0.9 7.3 0.70 1.1 36

0.7

2.1

0.6

0.79 0.91 0.15 11 b.d. 0.05 13.2 0.5 # 0.9 925 40

0.718 0.866 0.09 16 0.03 0.05 14.3 0.8

0.668 0.845 0.22 27 0.12 0.08 12.2 0.6 58 1.0 1170 57

2.3 1000 45

NL-/ Newlands Cor 840 840 205 475 96_+3 365_+7 52+3 19_+1 11+1 26_+1 <2 53_+2 4_+0.5 <2 20+1.5 3.9_+1

3.6 0.5 6.4 1.25 1.7

#

#

#

KA64-5

# # #

1.0 2.3 3.8 0.37 2.4

74501

74505

Bobbejahn Cor/Ky

Bobbejahn Cor

720 540 b.d. b.d. 64_+2 509_+4 48_+2 28_+1 I6_+1 25_+1 4_+1 15-+0.5 19_+2 <1.5 7.5_+1 10_+1 0.8 #

#

0.2 0.438 0.79 0.30 14 0.02 0.09 10.2 1.4

# # 1060 49

0.638 0.88 0.39 63 0.05 0.11 9.4 1.1 # 1.1 1180 57

0.599 0.913 0.54 58 b.d. 0.05 14.3 0.4 58 2.1 1090 51

0.6 #

7.1 0.48 1.8 9 #

1.3

840 540 135 b.d. 43_+3 163_+3 83+3 48_+1 15_+1 25_+1 11_+2 188_+3 7_+2 <2 42_+4 12+1.5

0.9 #

8.0 0.58 1.6 4

# #

660 610 70 b.d. 54_+2 385_+6 50-+2 24_+1 15_+1 27_+2 3.5_+0.5 30_+2 44_+4 <2 7_+1 5_+1

SRV-/ Rovic San/Coes/Ky

3.8 0.58 1.7 17 #

0.7

0.3

0.563 0.873 0.51 59 0.01 0.05 12.6 0.6 60 1.5 I 185 58

0.549 0.866 0.50 47 0.20 0.08 9.1 0.7 67 0.9 1 180 58

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

112

All of these features also may contribute to scatter in the plots shown below.

geotherm (Fig. 1) ; these have at least as much uncertainty as the T estimates; they are included primarily to emphasise the significant differences between the B / G suite and the kimberlite-derived suites. Potential causes of disequilibrium include metasomatism by the host magma (e.g., basalt or kimberlite) or associated mantle events, which could preferentially change the trace-element concentrations in one phase, or change Fe 3 +/Fe 2 + (e.g., B allhaus et al., 1991 ) and hence lead to spurious temperature estimates. In addition, short-term heating in the mantle or in the magma may result in incomplete trace-element diffusion, although the zoning expected from this process has not been observed here. DI pairs may not have been precipitated at the same growth stage (Griffin et al., 1992, 1993) and thus may not represent equilibrium pairs.

4. Methods The proton microprobe analyses for Ni, Zn, Ga, Sr, Y and Zr have been carried out using the CSIRO's HIAF facility, as described by Griffin et al. (1992, 1993) and in more detail by Ryan et al. (1990). These procedures provide rapid, standardless quantitative analysis of volumes ca. 30 × 30 × 30/xm. Typically 310-point analyses were collected on each phase, any outliers (judged to result from analysis of inclusions or alteration material below the surface of the grain) were discarded, and the selected analyses were digitally summed to form a single spectrum with improved sta-

Table 2 Electron microprobe data for clinopyroxenes and garnets in Table 1 Diamond eclogites

Eastern Australia pyroxenites GNS1 GNWI DRIOI62 BMI71 Lakes Bullenmerri/Gnotuk, western Victoria

GN47

GN649(1)

69-27 38926 Delegate, N.S.W.

XMll Orapa

JJG892 Orapa

Cpx

Cpx

Cpx

Cpx

Cpx

Cpx

Cpx

Cpx

Cpx

Cpx

SiO2 TiO2 AI203 Cr203 FeO MnO MgO CaO Na20 K20

49.00 0.61 9.81 b.d. 5.08 b.d. 12.30 20.90 1.96 b.d.

51.99 0.42 7.27 0.10 4.84 0.08 13.74 19.52 1.93 b.d.

51.29 0.31 7.51 0.21 4.40 b.d. 14.09 21.36 1.39 b.d.

50.64 0.70 6.57 0.37 5.49 0.08 13.45 20.08 1.65 b.d.

51.72 0.23 5.29 0.18 3.79 0.07 15.25 22.70 0.76 b.d.

52.24 0.45 5.97 0.20 4.09 b.d. 14.24 2 1.80 1.34 b.d.

51.10 0.56 5.30 0.69 3.99 0.09 15.2I 21.74 0.84 b.d.

47.72 0.29 16.97 0.10 2.41 b.d. 10.01 20.20 2.41 b.d.

54.95 0.40 4.81 0.36 2.34 0.06 15.22 18.71 2.89 b.d.

53.99 0.65 7.83 0.05 4.49 0.08 I 1.73 16.26 4.47 0.06

Total

99.66

99.89

100.56

99.03

99.99

100.33

99.52

100.11

99.74

99.61

Sample Locality

Gnt

Gnt

Gnt

Gnt

Gnt

Gnt

Gnt

Gnt

Gnt

Gnt

SiO2 TiO2 A120~ CDO3 FeO MnO MgO CaO Na20

40.30 0.10 23.10 b.d. 14.50 0.34 14.40 7.81 b.d.

41.21 0.08 23.13 b.d. 13.59 0.29 16.89 5.20 b.d.

42.45 0.06 23.95 0.14 12.44 0.26 17.36 6.07 b.d.

40.85 0.13 22.90 0.43 14.85 0.40 15.30 5.32 b.d.

41.29 b.d. 23.59 0.27 12.35 0.51 17.20 5.65 b.d.

41.27 0.07 23.33 0.14 13.39 0.42 16.36 5.41 b.d.

41.52 0.16 23.14 0.99 10.14 0.31 18.06 5.91 b.d.

42.17 0.07 24.28 b.d. 7.20 0.20 18.38 8.38 b.d.

41.64 0.15 23.91 0.25 9.78 0.35 19.95 3.82 0.16

41.15 0.63 22.99 0.09 13.08 0.36 15.05 6.52 0.16

Total

100.55

100.39

102.73

100.18

100.86

100.39

100.23

100.68

100.01

100.03

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

113

5. Results

tistics. The analytical uncertainties s h o w n in T a b l e 1 are g i v e n as one standard deviation m e a s u r e d f r o m the counting statistics on the s u m m e d spectrum, whereas the m i n i m u m detection limits are g i v e n at the 99% c o n f i d e n c e limit. E l e c t r o n m i c r o p r o b e ( E M P ) analyses for m a j o r elements, K, N a and Ti h a v e b e e n carried out at M a c q u a r i e University, the Australian National U n i v e r s i t y and the U n i v e r s i t y o f C a p e T o w n , using standard techniques. Basic E M P data for the samples s h o w n in T a b l e 1 are g i v e n in T a b l e 2, to facilitate calculation o f other parameters o f potential interest. E M P data for other suites are contained in the cited references.

5.1. Major-element parameters In analyzing the t r a c e - e l e m e n t data, we have found that the m a j o r - e l e m e n t c o m p o s i t i o n o f the garnet ( G n t ) and c l i n o p y r o x e n e ( C p x ) c o m m o n l y has a strong effect on t r a c e - e l e m e n t partitioning e v e n at very low concentrations. T h e s e effects are c o m p l e x , but apparently can be represented adequately by reference to three parameters (Figs. 2 - 4 ) . xCpx ranges f r o m 0.55 to 0.98, Mg [ M g / ( M g + F e ) ] with m o s t values falling b e t w e e n 0.7 and 0.95. This range reflects variation both in bulk c o m p o s i t i o n , and in temperature. Within a limited range o f bulk composition, the temperature effect on F e / M g partitioning

Diamond eclogites (cont.)

African diamond inclusions

JJG889 Orapa

AKI-IOc Orapa

Cpx

Cpx

55.82 0.47 9.47 b.d. 4.68 0.07 10.48 13.07 5.93 0.03

54.80 0.09 9.90 0.08 2.15 b.d. 12.50 16.60 4.18 0.11

55.00 0.24 5.13 0.07 2.51 0.03 14.80 19.20 3.01 0.05

55.20 0.31 7.80 0.07 2.61 0.05 12.60 17.50 4.65 0.06

55.78 0.40 5.19 0.06 5.22 0.11 14.56 15.10 3.63 0.06

56.19 0.39 10.24 0.10 4.77 0.04 9.78 13.56 5.55 0.18

100.02

100.41

100.04

100.85

100,11

Gnt

Gnt

Gnt

Gnt

Gnt

HRV247b Roberts Victor Cpx

HRV247F Roberts Victor Cpx

CBS17327 Newlands

JJG144 Newlands

PREM45 Premier

PREMI05 ORAPA4 Premier Orapa

JWA20 Jwaneng

JWA24 Jwaneng

Cpx

Cpx

Cpx

Cpx

Cpx

Cpx

Cpx

54.26 0.47 9.19 0.06 7.63 0.10 10.34 14.40 3.93 0.18

54.99 0.54 12.39 b.d. 5.84 0.08 8.04 11.83 6.18 0.19

54.48 0.48 7.92 0.12 4.92 0.12 11.64 15.05 3.47 1.12

54.79 0.62 11.32 0.05 6.27 0.03 7.77 12.21 6.22 0.06

54.28 0.60 9.68 0.03 6.22 0.05 9.28 13.63 5.18 0.07

100.80

100.56

100.08

99.32

99.34

99.02

Gnt

Gnt

Gnt

Gnt

Gnt

Gnt

40.16 0.38 22.78 0,01 18,26 0,43 11.69 7.18 0.13

41.60 0.19 23.30 0.12 6.97 0.15 14.70 12.70 0.12

41.60 0.27 23.60 0.06 10.40 0.26 16.70 7.45 0.09

40.60 0.38 23.00 0.05 11.80 0.22 13.90 9.80 0.11

41.52 0.40 23.43 0.08 13.59 0.30 17.71 3.33 0.11

39.75 0.28 22.81 0.05 t8.98 0.32 11.19 7.06 0.11

39.68 0.79 21.70 0.09 16.96 0.28 11.34 8.50 0.20

39.64 0.59 22.34 0.04 17.69 0.31 10.13 8.16 0.25

40.08 0.43 22.49 0.13 15.29 0.47 13.04 7.18 0.19

38.98 0.59 21.80 0.03 21.14 0.41 7.51 9.29 0.21

39.50 0.56 22.13 0.07 19.85 0.43 9.02 8.46 0.19

101.02

99.85

100.43

99.86

100.47

100.55

99.54

99.15

99.30

99.96

100.21

S.Y. O'Reilly, W.L. Griffin/Chemical Geology 12l (1995) 105-130

114

between Grit and Cpx (Ellis and Green, 1979) will produce a general decrease in XC~" with increasing T; this is in fact observed in the South African and Australian xenolith suites (Fig. 2), despite the scatter induced by the variation in bulk M g / ( M g + Fe) within each suite. However, the diamond-inclusion suites show no obvious correlation between vL'Mg cpx and T, which suggests they are derived from a wide range of bulk compositions. X cp× ranges from zero to 0.7 (Fig. 3); the highest values are found in the pyroxenes of some grospydites and kyanite- or corundum-bearing eclogites. In the absence of a buffering assemblage (e.g., plagioclase + quartz) the bulk composition of the rock determines the maximum jadeite content of the pyroxene,

but the jadeite content of clinopyroxenes still is broadly dependent on pressure, because of the changes in Grit/ Cpx ratio with increasing P. The African xenoliths show a broad positive correlation between X cpX and T (Fig. 3) ; this reflects the rapid increase in P relative to T along the cratonic geotherm. The correlation is best displayed by the Kaalvallei eclogites and the "Rovic" suite, where two distinct trends are visible, corresponding to the Group-I (high-Jd) and Group-II (Iow-Jd) eclogites of MacGregor and Carter (1970). The B/G pyroxenes show lower X cp× at any T than those from the South African xenoliths; this reflects the higher geotherm (lower P at any T) of the B / G suite. Although the bulk Na content of the B / G xenoliths is comparable to those of many African eclogites, the

Table 2 (continued) Rovic eclogites Sample Locality

RV-A Rovic

RVIG Rovic

RV2G Rovic

RV3G Rovic

BD3699 Rovic

KA64-6 Rovic

74506 Bobbejahn

BDII75 Rovic

BDI188 Rovic

BDl186 Rovic

Cpx

Cpx

Cpx

Cpx

Cpx

Cpx

Cpx

Cpx

Cpx

Cpx

SiO2 TiO2 A1203 Cr203 FeO MnO MgO CaO Na20 K~O

54.83 0.29 18.76 0.00 1.86 0.00 4.41 10.45 9.98 0.11

55.94 0.39 8.62 0.11 6.50 0.13 10.73 12.88 5.44 0.12

55.72 0.35 8.96 0.11 4.76 0.08 11.28 13.64 5.35 0.13

56.06 0.34 9.15 0.15 4.03 0.05 10.99 14.74 5.07 0.17

55.27 0.37 14.37 0.12 2.47 0.03 7.15 11.70 7.37 0.30

55.32 0.36 8.78 0.10 4,19 0,07 10.60 14,05 5,59 0,15

55.10 0.61 6,79 0.14 4.58 0.04 12.47 15.89 4.66 0.02

56.51 0.23 9.62 0.28 2.37 b.d. 11.38 16.10 5.21 0.07

56.15 0.05 9.41 0.20 1,65 b.d. 12,02 18,09 4.22 0.02

54.46 0.19 10.94 0.03 4.73 0.04 8.21 14.09 6.62 0.01

Total

100.69

100.86

100.38

100.75

99.15

99.21

100.30

101.77

101,81

99.32

Grit

Gnt

Gnt

Grit

Gnt

Gnt

Gnt

Grit

Gnt

Gnt

SiO2 TiO_~ AI203 Cr203 FeO MnO MgO CaO Na20

40.03 0.42 22.43 0.00 14.84 0.32 11.03 10.02 0.13

39.77 0.27 22.94 0.06 19.00 0.58 14.02 4.38 0.12

41.51 0.26 22.45 0.21 15.89 0.42 15.85 4.49 0.11

41.11 0.32 22.62 0.16 15.01 0.47 15.00 6.27 0.12

40.42 0.35 22.51 0.11 14.56 0.28 10.51 11.09 0.19

42.00 0.28 22.83 0.06 15.00 0.48 14.83 5.52 0.11

41.53 0.36 23.04 0.09 14.02 0.41 17.33 4.79 0.12

42.08 0.09 23.69 0.25 12.09 0.23 16.29 8.16 0.07

42.14 0.04 23.92 0.23 I 1.08 0.26 14.95 10.67 0.05

39.12 0.16 22,02 0.04 17.76 0,32 9.20 10.23 0,08

Total

99.22

101.14

101.07

101.08

100.02

101.11

101.69

102.95

103.34

98.93

b.d. below detection limit.

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

higher T at any P results in a more aluminous Cpx, a higher Cpx/Gnt ratio, and a lower X cpx on average. The African DI pyroxenes show a clear decrease in X cox with increasing T; this suggests that T increased more rapidly than P for this suite, and is consistent with other indications that the diamonds grew during thermal pulses which affected rocks distributed over a relatively limited range of P (Griffin et al., 1992, 1993). The Argyle DI pyroxenes show no correlation between X cpx and T, but have generally higher X cpx at any T than the African DI pyroxenes; this presumably is a bulk-composition effect. xGnt Ca appears to show no correlation with T (Fig. 4) ; the fact that the B / G samples are indistinguishable from the African eclogites in this plot suggests that xG,t Ca also is independent of P, and is controlled primarily by bulk composition.

115

5.2. Nickel

Ni contents range from < 10 to 200 ppm in garnet, and from < 10 to 1500 ppm in clinopyroxene; the broad fan on the correlation plot (Fig. 5a) indicates a significant range in DNi [D;= (ppm of element i in C p x ) / (ppm of element i in Gnt)]. In both phases, there is a broad positive correlation between Ni content and XMg (Fig. 5b and c). The plot o f DNi VS. 1000/T (Fig. 5d) shows a broad negative correlation between DNi and T, showing that Ni preferentially enters the Cpx, and that this tendency increases as T decreases. A plot of DNi against ,l. vcpx Mg (Fig. 5e) shows an overall positive correlation; however, closer inspection of this plot reveals that some individual suites actually show a negative correlation between DNi andX cpx (most clearly shown by the diamond--eclogite suite), so that the plot consists

Rovic eclogites (cont.)

BDII91

74507

74510 Bobbejahn

NL-4 Newlands

Bultfontein

NL-1 Newlands

74505

Bobbejahn

NL-3 Newlands

74501

Rovic

Bobbejahn

Bobbejahn

SRV-I Rovic

Cpx

Cpx

Cpx

Cpx

Cpx

Cpx

Cpx

Cpx

Cpx

Cpx

54.99 0.19 9.97 0.05 5.43 0.03 8.74 14.39 6.21 0.02

54.94 0.34 3.04 1.84 2.14 0.06 15.93 18.89 2.59 b.d.

54.85 0.24 3.71 0.32 2.76 0.07 15.34 20.32 2.47 b.d.

55.44 0.22 4.31 0.36 4.12 0.02 14.90 18.11 2.99 0.03

55.28 0.15 7.50 b.d. 3.89 b.d. 11.91 16.96 4.38 0.12

52.00 0.25 6.91 0.02 6.32 0.02 10.97 20.32 2.49 0.02

56.26 0.14 19.41 0.07 1.24 0.00 5.11 8.85 9.45 0.05

55.43 0.09 17.95 b.d. 1.06 b.d. 6.28 10.94 8.55 b.d.

55.59 0.09 18.26 0.00 1.50 0.00 5.76 10.27 8.78 0.01

56.61 0.09 17.45 0.00 1.70 0.00 6.19 11.80 6.98 0.20

100.02

99.77

100.08

100.50

100.19

99.32

100.58

100.30

100.26

101.02

Gnt

Gnt

Gnt

Gnt

Gnt

KA64-5

Gnt

Gnt

Gnt

39.39 0.16 21.84 0.08 18.92 0.39 9.38 8.94 0.09

42.49 0.31 21.63 3.08 7.38 0.35 21.47 4.48 0.07

42.46 0.24 23.69 0.20 9.20 0.35 19.46 5.78 0.05

42.16 0.16 23.52 0.40 12.97 0.28 18.51 3.67 0.05

41.57 0.17 23.32 0.01 12.99 0.29 14.69 8.40 0.08

40.29 0.07 21.87 0.04 18.72 0.41 8.19 11.39 0.09

40.94 0.14 22.75 0.03 11.06 0.19 10.92 14.96 0.11

40.40 0.12 23.25

99.39

101.26

101.43

101.72

101.52

101.07

101.10

Gnt

Gnt

9.02 0.19 7.56 20.83 0.05

40.17 0.11 23.01 0.01 10.77 0.27 7.77 19.34 0.05

40.46 0.14 22.87 0.02 11.11 0.26 7.60 18.98 0.08

101.42

101.50

101.52

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

116

1500 • 7

1.0

1300 r

90o

1100 r

•

a

0.9

oa

I

.o

o• ~•~"~Oo ++o • ~x~

~.. 0.8 o

X

X

(°c)

T

0.7

X

~oo~

l

X +

~

XX

A.

X

06 X

0.5

~

T

0.5

0.6

0.7

X

Afr diam incl

•

Argyle diam incl

o A • O +

B/G Diam ecIog Kaalvallei Monastery Rovic r

i

X

0.8

0.9

5.3. Zinc 1.0

IO00/T,K Fig. 2. X~," plotted against 1000/T for all the suites of clinopyroxene-garnet pairs. Descriptions of these suites and a discussion of the plot are given in the text. B/G refers to garnet pyroxenites from the Bullenmerri/Gnotuk locality in eastern Australia and includes one sample from a New South Wales locality (Delegate). Rotfic refers to xenolith suites from the kimberlites at Roberts Victor, Newlands and Bobbejahn mines. Diamond-bearing eclogites come from the Roberts Victor, Orapa and Excelsior mines. Subsequent diagrams use the same symbols.

1500

80

1300

1100

900

(°C)

+ +

+

60

• x

x ++

•

o

i•

•~ ~0

40

~o

~ ~

o•

X X A I~X~

~

X

1500

60

oo~

1300 ,

1100 ,

•

+

#. + ~ . CII

o• °

i

i

0.6

0.7

i 0 0.8

900

(°C)

50

+

o

0.5

Zn contents range from < 3 to 150 ppm in both Gnt and Cpx; the correlation plot (Fig. 6a) shows a broad scatter around a 1:1 line. The highest Zn contents are found in the African diamond inclusions. The Argyle DI data are not used here, because the surfaces of many grains had been contaminated by Zn from the surgical tape used for scanning electron microscopy (SEM) mounting (Griffin et al., 1988). Zn contents show a broad correlation with X M g in both phases (Fig. 6b and c). There is no obvious correlation of Dz,, with either ~.Gnt (Fig. 6e and f). However, the high-Ca X M g o r z,~Ca X Ca G"t >0.40) garnets of grospydites tend to take up slightly more Zn than other garnets, so that these rocks typically have anomalously low Dzn values. The plot of Dz, vs. 1000/T (Fig. 6d) shows a broad band between Dz° = 0.7 and Dz, = 1.5, regardless of T, with the grospydites (and two low-Ca garnets) extending to lower Dz,. In detail, the eclogites from Monas-

4- 4-

O

0

the corrected regression shows a significantly reduced scatter. We conclude that the partitioning of Ni between Gnt and Cpx in these samples shows a moderate dependence on XMg,and a strong dependence on T. The B/G samples remain on the high side of the scatter even after correction for XMg Cp×, which might indicate a small inverse correlation between DNi and P at constant T. However, the possible effect is too small to be reliably demonstrated by this dataset.

II

, 0.9

4O

A

o

30

,~

¢0 O

20

1,0

1000/T,K

Fig. 3. Jadeite component in clinopyroxene plotted against 1000/T for all samples. Symbols used are the same as defined for Fig. 2.

Ib

X

~lI+ + I + [] Ii,~ x x xWr~X . x ~ + O

X

I

of several broad en 6chelon negative trends. In detail, samples with a similar T range fall along trends that can be described by the equation DNi(corr)= DN~+y(XC~"-0.6), where y = 0 for T = 1 4 0 0 1500°C, y = 10 for T = 1200-1300°C and y = 2 0 for T= < 1200°C. The corrected DNi-Values (omitting three values that became < 0) are plotted in Fig. 5f;

+

10

O

0 0.5

I~

X

+

'lI~.lO

+

~0~__

Z~

u

I

~

i

i

~

I

0.6

0.7

0.8

0.9

1.0

1000/T,K Fig. 4. Xc, in garnet plotted against 1000/T for all samples. Svmbt:.,ls used are the same as defined for Fig. 2.

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

117

200

1800 a

• 150

>¢

o. o t=

X

1200

_•100

•

"

-7

0

@ co ~ + ~ • ~o/~ . ~A+o o+=

600

xx

&I

50"

X

0

&

•

0

~o

x ~'-~o ,-1. mA= 'U ~O

100 NI In gnt

~x~°~

0

I

50

•

°

x x

0

.~• 0 & i=+o.~_+=

b

I

0.3

200

150

1800

x

100

•

+

= o

I

I

I

0.4

0.5

0.6 0.7 XMg, grit

1500

1300

1100

900

I

i

I

I

0.8

0.9

1.0

(°C) d

C

X

1200

10

~o

tJ a

c z

o ,w

600

: o

+

•'A°

t o •=

I I.nDNI - - 2.3 + 5.7(1000/T)

X

0

.1

I

0.5

0.6

0.7

0.8

0.9

1.0

I

0.5

I

0.6

XMg, cpx 30

00 ~

1500 ,

1300 i

I

0.7 0.8 1000~,K 1100 ,

R"2 = 0.433 I

0.9 900 (°C) ,

1.0

f

0

I,

2O +

t~

o

+6>° o ,, 10 x .

0 0.5

, 0.6

a,

,Q 0.8 XMg, cpx

0.7

, 0.9

1.0

0.5

0.6

0.7

0.8

0.9

1.0

1 O00/'F,K

Fig. 5. Variation o f Ni p p m a n d / o r DNi with key chemical variables a n d 1000/T. DNi(COIT) =DNi +y(XMg Cpx - - 0 . 6 ) , where y = 0 for T = 1 4 0 0 1500°C v = 10 for T = 1 2 0 0 - 1 3 0 0 ° C and y = 20 for T = < 1200°C, and is discussed in the text. Symbols used are the s a m e as defined for Fig. 2.

118

S. E 0 'Reilly. W.L. Grig'fin/ Chemical Geology 12I (1995) 105 130

150

150

X X x

O

~.100

X

X

X

x~

~ 100'

tX

O C

# X

~

D

A

~×

C

C

N

so

5O

, 0

50

0

150

100

200

0.5

",

,

, 0.6

0.7

.

o

0.8

~,r.

0.9

XMg, cpx

Zn In gnt 1500

10

1300

1100

900 (°C i

d

150

X X X

X

X x

+

X

~ 100

,5

¢.,, 1

C

X

N

+O

113 O ~ 1 ~

5O

0

I

0.3

+

+O+

0.4

I

0.5

"~

A

I

i

0.6

0.7

0.8

0.9

.1 0.5

I

I

0.6

i

0.7

XMg, gnt

0.8

0.9

1000/T,K

10

10:

~

t-

#1

mn

++ +

¢#

A ,,

•

+

t

.1

i

0.3

0.4

I

0.5

I

0.6

XMg, gnt

I

0.7

.1

i

0.8

0.9

i

0

.10

+

I

.20

i

.30

i

.40

i

.50

.60

XCa, gnt

Fig. 6. Plots of Zn ppm in clinopyroxene and garnet and/or D z" against key chemical parameters and 1000/T. Symbols used are the same as defined for Fig. 2.

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

tery Mine, the Group-I eclogites from Kaalvallei and the African DI pairs show weak positive correlations between Dzn and T. The B / G samples fall in with the other suites, and show no dependence of Dzn on T; this suggests that the effect of P on Dzn is small. We conclude that Dzn is typically ~ 0.9 + 0.1 ; it may increase slightly with T, but is independent of P. The only identifiable bulk-composition effect is the greater uptake of Zn by some high-Ca garnets. These regularities make the partitioning of Zn a sensitive test for contamination and for disequilibrium induced by metasomatism. 5.4. Gallium

Ga contents range from 3 to 46 ppm in Cpx and from 4 to 25 ppm in garnet. Cpx typically contains more Ga than coexisting Gnt (Fig. 7a), but in detail some individual suites scatter across the general trend, reflecting a large range in DGa. A sapphirine granulite xenolith from Delegate (Griffin and O'Reilly, 1986) is anomalous in that the garnet contains 45 ppm Ga, whereas the coexisting Cpx contains only 8 ppm. It is interesting to note that Ga does not "follow" AI, in the sense of being more abundant in the more aluminous phase, even though the chemical "coherence" of Ga and AI is widely accepted in the literature (e.g., McKay and Mitchell, 1988). Ga contents of Cpx are well correlated with X cpx both overall and within individual suites (Fig. 7b), suggesting that Ga is involved in the chargebalance mechanism for Na substitution, as noted by Griffin et al. (1988). The B / G samples fall above the general trend (i.e. high Ga relative to Jd), which suggests that Ga also may be involved in the Tschermak's substitution together with A1. The plot of DG~ vs. 1000/T (Fig. 7c) shows a broad scatter with an overall negative slope (i.e. DGa increases with T). However, this trend is largely defined by some Argyle DI data at high T, and some Kaalvallei data at low T; most of the data lie in a band that shows a weak positive slope. In particular, the data from Monastery and Roberts Victor and the diamond eclogites show positive slopes (although there are many outliers). As might be expected, there is a strong overall correlation ofDG= withX cpx (Fig. 7d); the data from B / G lie well off this trend, and the data from the Argyle DI lie at a steep angle to it, since they show a wide variation in Dca at uniformly high X cpx. Separa-

119

tion of the data on this plot by T shows that DGabroadly increases with decreasing T at constant X cpX. We have used the data for the South African eclogites, in the range T = 1100-1200°C, to derive a correction for the effect of X cpx on D~a at constant T: DGa(corr) = D ~ a - 0 . 0 2 ( % J d ) . This simple procedure corrects all D~a to X cpX = 0; only 4 values give Dca < 0. A plot of D ~ ( c o r r ) vs. 1000/T (Fig. 7e and f) actually shows more overall scatter than a plot of the uncorrected data (Fig. 7d). However, the scatter within each suite is reduced, and several interesting points stand out. The B / G data show no clear T dependence, but 6 of the 7 points lie well above the African eclogite trend. The diamond-inclusion data from Argyle lie along a trend normal to that defined by the African eclogites, while the African DI data show no clear correlation of DGa with T (Fig. 7f). These contrasting trends can be reconciled by the following interpretation (Fig. 7f): (1) The T dependence over a narrow P range is given by the Argyle DI data, i.e. D~a increases with T; (2) the P dependence at constant T is shown by the difference between the African eclogites and the B / G suite, i.e. DGa increases as P decreases; (3) the broad decrease in Doa with increasing T observed in the African eclogite suites reflects the covariation of T with P along a cratonic geotherm, i.e. P increases relatively more rapidly than T, as noted above. The large degree of scatter within the African xenolith and DI data may reflect differences between the geotherms in each locality at the time of eruption, a larger degree of disequilibrium than is suggested by the other data, or compositional controls that have not been recognized. We conclude that D ~ is strongly dependent on xCpx Jd , because Ga is involved in the charge-balance mechanism for Na (and for tetrahedral AI-Si substitution in high-Tpyroxenes). Doa also shows strong P T dependence; Ga preferentially enters Cpx with increasing T and Gnt with increasing P. 5.5. Strontium

Sr is strongly concentrated in Cpx, where contents range from 3 to 1000 ppm (Fig. 8a); few garnets contain > 10 ppm Sr except for the Argyle DI garnets, in which Sr concentrations can reach 60 ppm. There obviously is a large range in Dsr, but the uncertainty in the ratio is large in many cases due to the analytical uncer-

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

120

50

50

a

+ b 40 -~

40

M

~

+

~- 30

~a,. M X

t•+

30

X •

+&

m

oe~ •I

¢~ 20

10

oo~d

~-~-a

m

l

~

+,

20

"

10

-O

+ ~ I

maD

• l

I

I

10

20

30

•

1500

1300

1100

I

I

I

+

+

•

•

•

0

I

20

I

40 %Jd, cpx

Ga In gnt 10

X

I

50

+

X

+

•

0

I

40

O /m ~.==.

Oooo 0

x

60

80

900 (°C) I

+

d +

,

3-

,+

o

•

al "~ 1

0

1 : ~ 2 - °°°

+

O

ml~ ¢ ~ X = , •

I

I

0.6

I

I

0.7

0.8

0.9

+!~

1500

1300

x+

:

x

I

0

I

20

I

40

60

80

%Jd, cpx

I~0~,K 10

+

n

0

1.0

a+J=+~ x~

+ & J=l&~lB --~P,m O +

1-

.1 0.5

•

0

900 (°C)

1100

O

1500

1300

1100

I

I

I

900 (°C) I

Argyle •

b 0 0

o

+ o~O,o

1 _ x _= •

o ~

~ ~ ( ~ G s A F

o

•

+

n .I

I

0.5

0.6

I

0.7

1000/F,K

I

0.8

I

0.9

1.0

.1 0.5

I

0.6

I

0.7

I

0.8

I

0.9

1.0

1000/F,K

Fig. 7. Plots o f G a ( p p m ) a n d / o r Dc~, in c l i n o p y r o x e n e a g a i n s t Ga ( p p m ) in garnet, %Jd in c l i n o p y r o x e n e and 1000/T. Dc,,(corr) = D m . - 0 . 0 2 ( % J d ) is plotted against 1000/T in (e), and (f) clarifies that plot by showing the trends of selected suites. Symbols used are the same as defined for Fig. 2.

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

1500

1000

121

1300

900 (°C)

1100

1000

i

b

a •&

800 o

&

100 x

t,,.x

600 +

t,,.,,. (,q

•

200

+

X

10

•

X

•A~+

x

A A

x ÷

400

+

!

0

|

10

I

20

'

!

I

30 40 Sr in gnt

-

' .....

50

1

!

!

60

0.5

!

0.7

0.6

0.8

0.9

1000rr,K 300

300

d O

2O0

200 C3

o

0

T < 1 0 2 5 °C

•

T 1025-1175 °C

•

T > 1 1 7 5 oc

&

100

o

100

++ ×

o~+

~

x~¢~o" +

0 • v;~

0

AA•

a

÷ I

.20

.10

0

A~ 8=•

•

I

.30 .40 XCa, gnt

+

.50

+.

•

n

i

m

•

|

•

Oqg

0

.60

A

I

0

•

.10

I

•

.20

i

k•

i

I

.30

.40

.50

.60

XCa, gnt 1500

300

1300

900 (oc

1100

i &

f

1000 1

/,,

&A

P 200

/,

o

X

o

I

x &

¢h

x

lOO

=.,"i,,i.+.. o -

0

+ -'

I

20

• " ,

+

I

+

+0

+

,... i

40 %Jd, cpx

10 1

i

60

80

0.5

i

0.6

"

•

i

0.7 1000/F,K

•

1

0.8

0.9

Fig. 8. Plots of Sr ppm in clinopyroxene and garnet and Dsr against relevant parameters ~LSdiscussed in the text. :~vmbols used are the same as defined for Fig. 2. (d) shows the "Rovic" dataset subdivided on the basis ofT.

122

s.Y. 0 'Reilly. W.L. Gri/fin / Chemical Geology 121 (1995) 105-130

tainty in the analysis of Sr in Gnt (Table 1). Excluding some outliers (especially the low-temperature members of the " R o v i c " suite), Ds,. shows a broad overall negative correlation with T (Fig. 8b). There also is a strong negative correlation between Ds~• and ~vcmc~,,at least up to ~ 25% Gross (Fig. 8c), reflecting a greater uptake of Sr by Ca-rich garnets over this range. In general, the high-T samples within the " R o v i c " eclogite suite have higher Ds,. at a given X~.',; ~ than the lowerT garnets (Fig. 8d). Using the medium-T group from this suite, we derive a correction for the effect oi'"-'~ .,~m c,, nt on Ds~: D s r ( c o r r ) = D s , . + 1." 3. ~. y. G c~, • After removal of this effect, there is still a weak negative correlation between Ds,. and X cv~ at low Jd contents (Fig. 8e) but the data are too scattered to provide a reliable correction for this effect. The plot ofDs~(corr ) vs. 1000/T (Fig. 8f) shows a reduced scatter (compared to Fig. 8b), and hence a somewhat more pronounced negative correlation of Dsr with T; however, this trend is not definitive within any one suite, with the possible exception of the diamond eclogites (n = 4) and the African diamond inclusions. We conclude that Ds~ is strongly controlled by X~'~tc+,,especially at low grossular contents. There also appears to be a weak T effect on Ds~, with more Sr entering Gnt as T increases. No P effect is discernible in the limited data available from the B / G suite. 5.6. Y t t r i u m

The Y contents of Gnt range from < 2 - 1 2 0 ppm; most Cpx contain < l0 ppm and many contain < 1 ppm (Fig. 9a). There clearly is a large range of Dr, but as for St, the low Y contents of the Cpx produce relatively large uncertainties on some Dv-values. There is a good correlation between the Y content of Gnt and XC,,t in some suites, notably the African DI and the Ca Monastery eclogites, but no correlation is obvious in the other suites (Fig. 9b). There is no apparent correlation between Y and Na shown in the plot of Y vs. %Jd (Fig. 9d). The behaviour of Dy with T varies widely among the sample suites (Fig. 9c) ; the African DI show a broad positive correlation between Dy and 7", the Monastery suite shows a negative correlation, and the other suites show none. Overall, the data define a large scatter with a slight negative correlation between Dy and 7+. However, there is an overall negative correlation between Dy and xC'mc~,(Fig. 9d), which

is essentially independent of T. We have used the African D1 data to derive a correction for the effect of x Gnt Ca on Dy: D v ( c o r r ) = D y + 140X(C~p~. Removal of this effect greatly reduces the scatter in the data, and leaves a very small positive correlation between Dy and T (Fig. 9e). The B / G data/'all within the range of the other suites. We conclude that the partitioning of Y is controlled almost entirely by ~c~, vG,,t , and that there is only a weak tendency for Y to enter Cpx with increasing T; no P effect on partitioning is discernible. 5. 7. Z i r c o n i u m

Zr contents generally range from < 1 to 50 ppm in Cpx and from < 1 to 80 ppm in Gnt, with generally higher contents found in Gnt (Fig. 10a). Several Monastery eclogites have unusually high-Zr Cpx, and the Argyle DI garnets tend to contain anomalously high Zr. Dzr for most samples ranges from 0.1 to 3, and shows an overall negative correlation with T; a number of samples, especially from Monastery, scatter widely around this main trend (Fig. 10b). There is a broad overall negative correlation between Dz,. and L,vCpXjd (Fig. 10c), and this correlation is pronounced within several individual suites, notably the Argyle DI, B/G, and the eclogites from Kaalvallei and Roberts Victor. Examination of the African eclogite suite shows that in general Dz,. increases with Tat any X cox . We have used the intermediate-T group from this suite to derive a correction for the effect of vq,x z . Jd on Dz,-: Dzr(corr) = Dzr + 14X cpX . After removal of this effect, there is no clear overall correlation between Dz,.(corr) and l,vc;"c:~+(Fig. 10d). No correction has been applied for this effect; similarly, no correlation could be seen between Dz, (corr) and vcp× (not shown). ~Mg The plot ofDzr(corr ) vs. 1000/T (Fig. 10e) shows a greatly reduced scatter compared with the raw data (Fig. 10b). Some of the high-T, high-Jd pairs from Argyle may have been overcorrected by the procedure applied. Two Monastery eclogites still fall well above the trend; these may reflect disequilibrium caused by late-stage metasomatic uptake of Zr by Cpx. The B / G samples are not distinguishable from the other suites. We conclude that Dz~ is mainly controlled by T, with more Zr entering Gnt as T increases. Dz,. is lowered by Gt+t increased ,.vcP×ju, and possibly also by increased V~c,, •

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

10

123

a

120 O

8

~

El

x X

n

U r-

ra

o

>"

m"'m

4

X X)NK XO

~

~

+

o

n

D

m

0

Q

o

•~

X

4

I~&

+

2 O+

~'+

o ~

0

I

I

I

1500

1300

i

t

"÷

A I

.20

+ I

.30

•

I

.40

.50

.60

50

60

XCa, gnt

Y In gnt

10

m

•

!

.10

120

80

4O

I

900 (°C)

1100

10

I

c n

8 -

>
+

1

~.6U e-

>-

x

o

x

oC~O

¥&xx

0O

.1

+.~

~

0

4-~

Ox~

o O

El

X X

0o

•

I

O

o

X

X

[3

X

^

~

~

z$ + I

0

0.9

0.8

X IP,-

ZX

0

I

0.7

X

X[:3~

&

2-

i

0.6

0.5

El

r'l O >0<

O

.01

X

O

10

20

30

40

% J d , cpx

1000/T,K

0.5

1500

1300

1100

I

I

]

900 (°C) I

f 0.4-

o

x

>. 0.3

+x M oox XOo

0.2

0

0.1

Oa

0 ~

0

[]

1

^ft.. ~_ a ¢~_a~~ "-""~"~d'~o~%

o =

n

x +

.1

X EIA ~

D&

0.0

'°1

[]

)0(

130

~0

LnDy = - 0.59 - 0.43(1000/1") R^2 = 0.005

,

,

,

,

,

.10

.20

.30

.40

.50

XCa,

gnt

.01 .60

I

0.5

0.6

i

0.7 1 000/T,K

i

0.8

0.9

Fig. 9. Plots of Y (ppm) in clinopyroxene and garnet against key chemical parameters, and Dy against 1000/T. Dv(corr) = Dv + 140X~'," (as discussed in the text). Symbols used are the same as defined for Fig. 2.

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

124

100 :

200

1500

1300

1100

I

I

I

900(°C) I

0 150

[]

[]

t21

rl

[]

g

-=loo z~

50

[]

10:

X

~mo

t

1 +

+ x .o

~

[]

[] •

!

I

0

.1

I

I

50

150

100 Zr In gnt

A-O-

I

I

0.6

0.5

200

X

1

0.7

0.8

0.9

1000/T,K

20

5 4-

15

& •

3-

•

& + +

0

4+

o +

•

m~ •

o,oo .

++ + I

I

20

40 %Jd, cpx

1500

100

1300

+ I

I

'

60

8O

0

.10

I

.20

I

.30 XCa, gnt

I

.40

I

.50

.60

900 (°C)

11 O0

I

O

!= 10

[]

8

.

Im

L,,

N e~

i

0.5

xVO

0.6

- +

o_

0.7 1000/T,K

_

•

0.8

0.9

Fig. 10. Variation o f Z r ( p p m ) and Dzr with key chemical parameters and 1000/T. Dzr(Corr) = D z , + 14X~I as discussed in the text. Symbols used are the same as defined for Fig. 2. TM

S. E 0 'Reilly, W.L. Griffin / Chemical Geology 121 (1995) 105-130

No pressure effect on Dz,- can be recognized from these data.

125

6. Discussion and conclusions

6.1. Compositional dependence 5.8. Titanium

Ti contents range up to 5500 ppm in Cpx and 8000 ppm in Gnt; in most samples, the two phases have similar Ti contents (Fig. 1la). However, the Ti contents of the B / G garnets are uniformly low, despite the relatively high Ti in Cpx. The Kaalvallei eclogites show two distinct groups on the basis of Ti content. There appears in general to be little dependence of D-ri on T (Fig. 1lb), but the diamond-eclogite suite displays a pronounced negative correlation between DTi and T. Dvi appears in general to be independent of-Xc~ cnt (Fig. 1 ld), but a group of relatively high-Ca samples in the Argyle and diamond-eclogite suites have anomalously low Dye. There is no overall correlation between Dw and X cpx, but a strong positive correlation is seen within the diamond--eclogite, Monastery and Kaalvallei suites, suggesting that X cpx may be an important factor. Inspection of the data shows that for a given Jd content, D-r~ decreases with increasing T. Using the medium-T subset of the African eclogite suite, we have derived a correction for the effect of X cpx on D-ri: DTi (corr) = D w - 9 2 Y co x + 0.55. This algorithm cor~" ~ J d ycp~ r e c t s O y i t o a constant ,~ j~ = 0.20. However, application of this correction produces no improvement in the overall correlation between DTi and 1000/T; several low-Ca samples from Argyle still fall offthe main trend to lower Dw (Fig. 1 le). Both raw and corrected Dw show a broad negative correlation with T, suggesting an increased preference of Ti for Grit at higher T. The extremely low Ti contents of the B/G garnets (leading to high /)TO suggest that D-r~ decreases significantly with increasing P at constant T, although this effect cannot be quantified with the present dataset. It therefore is probable that the apparent dependence of Dv~ on T observed in Fig. 11 e simply reflects the relationship of T to P along a typical cratonic geotherm. We conclude that Dv~ decreases markedly with increasing P, and also may decrease with increasing T. These effects are partially offset by a weak tendency for Dr~ to increase with increasing X cp~ .

Among the elements studied here, Zn shows the least compositional effect on partitioning. This probably reflects the fact that Zn substitutes simply for Fe in both phases, and the Zn contents of coexisting Cpx and Gnt tend to be very similar. There is some indication that Dz, decreases as Xca,~t increases. The partitioning of Ni shows a weak dependence on XMg, reflecting the similarity of ionic radii between Ni and Mg. The partitioning of Y and Sr appears to be controlled essentially by the Ca content of Gnt. The substitution of Sr into garnet probably is controlled primarily by ionic-radius effects, where Sr enters the larger site provided by Ca. Caporuscio and Smyth (1990) showed that the larger X site of grossular-rich garnets favours the uptake of intermediate rare-earth elements (MREE) relative to heavy rare-earth elements (HREE), and the ionic radius of Y~+ is most similar to that of Ho -~+. The substitution of Y into clinopyroxene must involve a charge balance, but our data do not identify the mechanism; defect substitution may be involved (Wood, 1974; B. Hensen, pers. commun., 1993). The data of Caporuscio and Smyth (1990) suggest that partitioning of Y also should be affected by X cpX, but this effect is not obvious in our data, probably because of the overriding importance of XcC~, "t. The partitioning of Ga between Gnt and Cpx is very strongly affected by --J0 yCpx , as well as being controlled by T and P. The P - T dependence ,~f within each ~ l i xycp× jd bulk composition makes these effects difficult to separate. The strong correlation between the Na and Ga contents of Cpx suggests that Ga enters the Cpx lattice as part of the charge-balance mechanism. Both the Ca ~ Na and Na ~ [2]substitutions on the M2 site are compensated by trivalent cations substituting on the MI site, a contraction of the M1 polyhedron, and a more strongly negative M1 site potential (Caporuscio and Smyth, 1990; Oberti and Caporuscio, 1991 ) ; these features may favour the uptake of the Ga 3 ~ ion on M 1. The partitioning of Zr and Ti between Cpx and Grit v x . , this suggests that shows a weak dependence on vl,cjo a similar mechanism affects the uptake of these elements, even at low concentrations.

S. E O 'Reilly, W.L. Griffin / Chemical Geology 12 l (I 995) 105-130

126

6000

10 a X

5000

X

~00 0

x x

4000 O

+ A,

•

0

n

o ~ ~x., ~ • x

+

+

o

3000

0

[3 i-i

o

0

XX

P

,,il'

gl

x" ×

~+.~

A

2000 1000

++ i

0

2000

i

4000

I

6000 in gnt

10000

8000

.1 0.5

I

10 I

o o o

0,9

d

O

O

+

+ A~

+

I= a 4

08

o

O &

I

07 1000/T,K

oo°

4

0

I

06

[

+

•,r

+ 0

0

+

"+

, x t ++ + I

I

20

+

I

40 °/~1d, ¢px

,1

60

80

I

0

10

I

20

I

I

30 XCa, gnt

40

50

60

10 0

000

0

0

+

o +

O

&

~. ±x~O~+ ++

I.- 1 o o

Ilaoc

+ A

e~ .1

.01

all. A

i

0.5

0.6

i

0.7 1000/T,K

i

0.8

0.9

Fig. I 1. Ti (ppm) in clinopyroxene plotted against key chemical parameters and 1000/T. D n ( c o r r ) the effect of X.~'~Px on D.r,. Symbols used are the same as defined for Fig. 2.

= Dii-2.2XJ~

TM

+ 0.55 and accounts for

s. E 0 'Reilly, W.L. Griffin/ Chemical Geology 121 (1995) 105-130 6.2. Temperature dependence

There appears to be little temperature dependence on the partitioning of Zn, Sr and Y between Gnt and Cpx; Dzn is essentially constant, and Dsr and Dy are controlled primarily by crystal-chemical factors. The partitioning of Ni and Zr shows a significant T effect, which overrides the modest effects of phase composition; both elements partition preferentially into Grit at higher T. The partitioning of Ti may reflect a weak T effect, with D decreasing as T increases, but it is difficult to separate this from the pressure effect with the available data. The behaviour of the Argyle DI suite suggests that partitioning of Ga between Gnt and Cpx is strongly dependent on T, as well as on P and X cp~. At constant P and X cvX, D appears to increase with increasing T, as Ga preferentially enters Cpx. 6.3. Pressure dependence

The analysis in this study of the pressure dependence of partitioning is based on the difference in ambient geotherms for the African and B / G sample suites. For most elements, including Ni, Zn, Y, Sr and Zr, the B / G samples do not stand out from the African ones once compositional dependence has been removed. We therefore conclude that the partitioning of these elements between Grit and Cpx in eclogites and pyroxenites shows no measurable P dependence. The B / G garnets contain much less Ti than the coexisting Cpx, in contrast to the African samples with similar T, and this indicates the existence of a dependence of Dw on P. This P effect may also account for the apparent decrease in Dw with increasing T within the African suite, since P increases rapidly relative to Talong the ~ 40-mW-m 2 paleogeotherm from which these samples were probably derived. Smith et al. ( 1991 ) showed that DT~ is higher in peridotite xenoliths sampled from elevated geotherms than in those from cratonic kimberlites; this suggests that a similar P effect is present in ultramafic rocks. A similar situation applies in the case of Ga, after the dependence Of DG, on X cpX has been removed; the B / G samples show higher DGa than the African eclogites at similar T, implying that DGa decreases with increasing P. As for Ti, this P dependence may account for the general decrease of D with increasing T in the

127

African eclogite suites. The true T dependence of DG~ probably is shown by the Argyle DI pairs, where DG,~ increases sharply with increasing T. Scatter in these data may be due to lattice effects such as defect substitution, as strain may also be a significant factor but cannot be assessed by this dataset. 6.4. Applications of partitioning data

None of these elements promise simply applicable alternative geothermometers or geobarometers for eclogitic systems. The corrections derived here for compositional effects are necessarily crude, the scatters about the D-Tregressions are still large, and in the case of Ga and Ti the strong P effects are not quantified. However, these data do provide a series of tests for equilibrium between phases in xenolith samples, and hence can be used to improve the choice of appropriate samples for other studies, especially including isotopic analysis. Samples that fall well outside the trends shown here, after compositional effects are considered, probably are in disequilibrium, and the nature of the anomalies may provide guides to the processes that have produced the disequilibrium. The broad patterns of trace-element partitioning described here also provide information on the environment of formation of Gnt-Cpx rocks. The partitioning behaviour of Ga and Ti, in particular, is useful for recognizing eclogitic material from different P - T regimes. For example, the data on Ga partitioning in the Argyle DI pairs are consistent with the formation of these diamonds in thermal pulses occurring over relatively restricted ranges of P. Studies of other DI suites may reveal similar patterns. Eclogites or eclogitic DI pairs from asthenospheric P - T regimes, in which T increases adiabatically with P, should show anomalously low DG~and DTi at high T, relative to the' 'African eclogite" trends. Conversely, samples from strongly advective thermal regimes should stand out because of high DGa and Dw at relatively low T, as do the B / G samples. These data also can be used to evaluate the validity of experimental determinations of partition coefficients, and to guide the choice of experimental conditions and the compositions used in the experiments. For example, Green et al. (1989) report garnet/liquid and clinopyroxene/liquid partition coefficients for Zr and Y, derived from experiments on a basaltic composition

128

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

at 25 kbar and 1100°C. The garnets in the experimental charge are zoned, whereas the pyroxenes are homogeneous; Dzr (Cpx/Gnt) -values calculated from their data range from 0.14 to 0.17; by comparison with the data shown in Fig. 10d and e, these values are low by perhaps a factor of 10. The values of Dv calculated for the same samples range from 0.08 to 0.12; after correction for XcG2t these values appear to be low by perhaps a factor of 5 (cf. Fig. 9). It seems likely that traceelement equilibrium was not attained in these experiments, even in the rims of grains, or that the analyses are flawed. A reanalysis of one of these samples by laser ablation microprobe-inductively coupled plasma mass spectrometry ( L A M - I C P - M S ) (Jenner et al., 1993) gave higher values of Dzr = 0.51 and Dy = 0.17, which still are lower than the data from natural samples at this T. Fujimaki et al. (1984) give Dzr=0.67 for coexisting garnet + clinopyroxene produced in an olivine tholeiite composition at 20 kbar and 1150°C. No X cpx is given, but it is expected to be low, and the raw value of Dzr falls within the range of our data (Fig. 10). Shimizu (1980) measured D s r in clinopyroxenegarnet pairs produced by experiments on high-A1 basalt and alkali basalt at 30 kbar and 1400-1500°C. His three Ds~-values range from 3.6 to 6.7 and are inversely correlated with XcG2t, which ranges from 0.16 to 0.27; these values are consistent with the data shown in Fig. 8d and f. Similarly, Johnson (1993) reports data from a Kilauea basalt (Hawaii) composition at 1300-1470°C, 20-30 kbar (conditions well above even the southeastern Australian geotherm of Fig. 1): DTi = 1.28; Dsr= 29: Dy = 0.13; Dzr=0.18. While these data cannot be corrected for phase composition with the data given in the abstract, the individual values fall within the ranges of those reported here. More data from such experiments will help to quantify the P dependence of

indication of whether Fe-Mg temperatures derived by calculating Fe 3+ in Gnt and Cpx are realistic. The general ranges of D~ established here for given values of P and T, and the qualitative dependence of D on bulk composition, provide important constraints on the values used in the mathematical modelling of partial melting, fractional crystallization and metasomatic processes in the mantle. Experimental data on clinopyroxene/liquid partition coefficients are much more abundant than data on garnet/liquid coefficients; our data may be used to calculate appropriate garnet/liquid trace-element coefficients from the clinopyroxene/liquid data, and thus greatly extend the potential use of trace elements in modelling involving basaltic compositions at high pressures.

Acknowledgements John Gurney generously provided samples of diamondiferous eclogites and diamond-inclusion pairs which contributed significantly to the database, and Fanus Viljoen kindly provided samples from Roberts Victor and Kaalvallei. We have benefited from comments on the early drafts by Norm Pearson, Bas Hensen, Jing Guo, Ming Zhang, Heinz Stosch, Chris Ryan and Dmitri Ionov, and we thank them for their time and helpful discussions. Gordon Medaris and Doug Smith provided constructive reviews of the penultimate version. Jing Guo assisted with the data handling (a large task) and the diagrams. Chris Ryan provided invaluable help with the proton microprobe techniques and Tin Tin Win assisted with the proton microprobe analyses. This research was supported financially by the Australian Research Council and the Macquarie University Research Grants Scheme.

DTi.

As noted earlier, the estimation of temperature in Gnt-Cpx rocks, using Fe-Mg geothermometry, requires the calculation of Fe3+/Fe 2+ from stoichiometry. Since this calculation accumulates all of the analytical errors, the F e 3 + values typically are overestimated; Fe 3 + also may be affected by late-stage alteration processes. The trace-element data may be used to cross-check the Fe 3÷ calculations. Specifically, the DNi-T and DZr-T relations described here can help to constrain the equilibration temperatures, and give an

References Ballhaus, C., Berry, R.F. and Green, D.H., 1991. High pressure experimental calibration of the olivine-orthopyroxene-spinel oxygen barometer - - implications for redox conditions in the upper mantle. Contrib. Mineral. Petrol., 85: 85-94. Boyd, F.R., 1976. Inflected and non-inflected geotherms. Carnegie Inst. Washington Yearbk., 75: 521-523. Caporuscio, F.A. and Smyth, J.R., 1990. Trace element crystal chemistry of mantle eclogites. Contrib. Mineral. Petrol., 105: 550561.

S.Y. O'Reilly, W.L. Griffin / Chemical Geology 121 (1995) 105-130 Ellis, D.J. and Green, D.H., 1979. An experimental study of the effect of Ca upon garnet-clinopyroxene Fe-Mg exchange equilibria. Contrib. Mineral. Petrol., 71: 13-22. Finnerty, A.A. and Boyd, F.R., 1987. Thermobarometry for garnet peridotite xenoliths: a basis for mantle stratigraphy. In: P.H. Nixon (Editor), Mantle Xenoliths. Wiley, New York, N.Y., pp. 381-402. Fujimaki, H., Tatsumoto, M. and Aoki, K., 1984. Partition coefficients of Hf, Zr and REE between phenocrysts and groundmasses. Proc. 14th Lunar Conf., Part 2. J. Geophys. Res., 89(Suppl.): B662-B672. Hatton, C.J. and Gurney, J.J., 1987. Roberts Victor eclogites and their relation to the mantle. In: P.H. Nixon (Editor), Mantle Xenoliths. Wiley, New York, N.Y., pp. 453-463. Green, T.H., Sie, S.H., Ryan, C.G. and Cousens, D.R., 1989. Proton microprobe determined partitioning of Nb, Ta, Zr, Sr and Y between garnet, clinopyroxene and basaltic magma at high pressure and temperature. Chem. Geol., 74: 201-216. Griffin, W.L. and O'Reilly, S.Y., 1986. Sapphirine in a mantlederived xenolith from Delegate, Australia. Mineral. Mag., 50: 635-640. Griffin, W.L., Wass, S.Y. and Hollis, J.D., 1984. Ultramafic xenoliths from Bullenmerri and Gnotuk maars, Victoria, Australia: petrology of a sub-continental crust-mantle transition. J. Petrol., 25: 53-87. Griffin, W.L., Jaques, A.L., Sie, S.H., Ryan, C.G., Cousens, D.R. and Suter, G.F., 1988. Conditions of diamond growth: a proton microprobe study of inclusions in West Australian diamonds. Contrib. Mineral. Petrol., 99: 143-152. Griffin, W.L., Gurney, J.J. and Ryan, C.G., 1992. Variations in trapping temperatures and trace elements in peridotite-suite inclusions from African diamonds: evidence for two inclusion suites, and implications for lithosphere stratigraphy. Contrib. Mineral. Petrol., 110: 1-15. Griffin, W.L., Sobolev, N.V., Ryan, C.G., Pokhilenko, N.P., Win, T.T. and Yefimov, Y., 1993. Trace elements in garnets and chromites: diamond formation in the Siberian lithosphere. Lithos, 29: 235-256. Gurney, J.J., Hams, J.W. and Rickard, R.S., 1979. Silicate and oxide inclusions in diamonds from the Finsch kimberlite pipe. In: F.R. Boyd and H.O.A. Meyer (Editors), Kimberlites, Diatremes and Diamonds. Am. Geophys. Union, Washington, D.C., pp. 1-15. Gurney, J.J., Harris, J.W. and Rickard, R.S., 1986a. Premier Mine diamond inclusions. Univ. of Cape Town, Stellenbosch, Kimberlite Res. Group, Intern. Rep. No. 2. Gurney, J.J., Hams, J.W., Rickard, R.S. and Moore, R.O., 1986b. Silicate and oxide inclusions in diamonds from the Orapa mine, Botswana. Univ. of Cape Town, Stellenbosch, Kimberlite Res. Group, Intern. Rep. No. 4. Jaques, A.L., Hall, A.E., Sheraton, J.D., Smith, C.B., Sun, S.-s., Drew, R.M., Foudalis, C. and Ellingsen, K., 1990. Composition of crystalline inclusions and C-isotopic composition of Argyle and Ellendale diamonds. In: Kimberlites and Related Rocks, Vol. 2: Their Mantle/Crust Setting. Geol. Soc. Aust., Spec. Publ. No. 14. Blackwells, Melbourne, Vic., pp. 966-989. Jenner, G.A., Foley, S.F., Jackson, S.E., Green, T.H., Fryer, B.J. and Longerich, H.P., 1993. Determination of partition coefficients

129

for trace elements in high pressure-temperature experimental run products by laser ablation microprobe-inductively coupled plasma-mass spectrometry (LAM-ICP-MS). Geochim. Cosmochim. Acta, 57:5099-5104. Johnson, K.T.M., 1993. Experimental cpx/ and garnet/melt partitionong of REE and other trace elements. Eos (Trans. Am. Geophys. Union), 1993 Fall Meet., p.658. Krogh, E.J., 1988. The garnet-clinopyroxene Fe-Mg geothermom e t e r - - a reinterpretation of existing experimental data. Contrib. Mineral. Petrol., 99: 44-48. MacGregor, I.H. and Carter, J.L., 1970. The chemistry of clinopyroxenes and garnets of eclogite and peridotite xenoliths from the Roberts Victor mine, South Africa. Phys. Earth Planet. Inter., 3: 391-397. McCandless, T.E. and Gurney, J.J., 1989. Sodium in garnet and potassium in clinopyroxene: criteria for classifying mantle eclogites. In: Kimberlites and Related Rocks, Vol. 2: Their Mantle/ Crust Setting. Geol. Soc. Aust., Spec. Publ., No. 14. Blackwells, Melbourne, Vic., pp. 827-832. McKay, D.B. and Mitchell, R.H., 1988. Abundance and distribution of gallium in some spinel and garnet lherzolites. Geochim. Cosmochim. Acta, 52: 2867-2870. Moore, R.O. and Gurney, J.J., 1989. Mineral inclusions in diamond from the Monastery kimberlite, South Africa. In: Kimberlites and Related Rocks, Vol. 2: Their Mantle/Crust Setting. Geol. Soc. Aust., Spec. Publ. No. 14. Blackwells, Melbourne, Vic., pp. 1029-1041. Moore, R.O., Gumey, J.J., Griffin, W.L. and Shimizu, N., 1991. Ultra-high pressure garnet inclusions in Monastery diamonds: trace element abundance patterns and conditions of origin. European. J. Mineral., 3: 213-230. Neumann, E.-R., 1976. Two refinements for the calculation of structural formulae for pyroxenes and amphiboles. Nor. Geol. Tidsskr., 56: 1-6. Oberti, R. and Caporuscio, F.A., 1991. Crystal chemistry of clinopyroxenes from mantle eclogites: a study of the key role of the M2 site population by means of crystal-structure refinement. Am. Mineral., 76:1141-1152. O'Reilly, S.Y., 1989. Australian xenolith types, distribution and transport. In: R.W. Johnson (Editor), lntraplate Volcanism in Eastern Australia and New Zealand. Cambridge University Press, Cambridge, pp. 249-253. O'Reilly, S.Y. and Griffin, W.L., 1985. A xenolith-derived geotherm for southeastern Australia and its geophysical implications. Tectonophysics, 111: 41-63. O'Reilly, S.Y., Griffin, W.L. and Ryan, C.G., 1991. Residence of trace elements in metasomatized spinel lherzolite xenoliths: a proton microprobe study. Contrib. Mineral. Petrol., 109: 98-~ 13. Pollack, H.N. and Chapman, D.S., 1964. On the regional variation of heat flow, geotherms and lithospheric thickness. Tectonophysics, 38: 279-296. Rickard, R.S., Hams, J.W., Gumey, J.J. and Cardoso, P., 1989. Mineral inclusions in diamond from the Koffiefontein Mine. In: Kimberlites and Related Rocks, Vol. 2: Their Mantle/Crust Setting. Geol. Soc. Aust., Spec. Publ. No. 14. Blackwells, Melbourne, Vic., pp. 1054-1062.

130

S. E O 'Reill_v, W.L. Gri~lin / Chemical Geology 121 (1995) 105-130

Robinson, D.N., Gurney, J.J. and Shee, S.R., 1984. Diamond eclogite and graphite eclogite xenoliths from Orapa, Botswana. In: J. Kornprobst (Editor), The Mantle and Crust-Mantle Relationships. Elsevier, Amsterdam, pp. 11-24. Ryan, C.G., Cousens, D.R., Sie, S.H., Griffin, W.L. and Suler, G.F., 1990. Quantitative PIXE microanalysis of geological material using the CSIRO proton microprobe. Nucl. lnstrum. Methods, B47:55-7 I. Shimizu, H., 1980. Experimental study on rare-earth element partitioning in minerals formed at 20 and 30 kb for basaltic systems. Geochem. J., 14: 185-202. Smith, D., Griffin, W.L., Ryan, C.G. and Sie, S.H., 1991. Traceelement zonation in garnets from The Thumb: Heating and melt

infiltration below the Colorado Plateau. Contrib. Mineral. Petrol., 107:60 79. Viljoen, K.S., 1994~ The petrology and geochemistry of a suite of mantle-derived eclogite xenoliths l~om the Kaalvallei kimberlite, South Africa. Ph.D. Thesis, University of the Witwatersrand, Johannesburg, 160 pp. Viljoen, K.S., Robinson, D.N., Swash, P.M., Griffin, W.L., Otter, M.L., Ryan, C.G. and Win, T.T., 1994. Diamond- and graphitebearing peridotite xenoliths from the Roberts Victor Kimberlite, South Africa. In: H.O.A. Meyer and O.H. Leonardos (Editors), Kimberlites, Related Rocks and Mantle Xenoliths. Comp. Pesq. Recurs. Miner., Spec. Publ., IA/93: 285-303. Wood, B.J., 1974. The solubility of alumina in orthopyroxene coexisting with garnet. Contrib. Mineral. Petrol., 46: 1-15.