RESIDUAL AND MOBILE MANTLE SULFIDES: LAM-MC-ICPMS ANALYSIS OF RE-OS ISOTOPES

O.Alard, N.J. Pearson, W.L. Griffin, S. Graham and S.E. Jackson

GEMOC, Macquarie

Introduction:

The Re-Os isotopic system is now widely used to study the timing and nature of melting and metasomatism in the lithospheric mantle. Petrographic and microanalytical data show that essentially all of the Os, and most of the Re, in mantle-derived peridotites is concentrated in sulfide phases. These sulfides occur as two basic types: (1) residual, Ni-rich monosulfide solid solutions (now largely exsolved), that are generally enclosed inside primary silicates, oxides and diamonds; (2) interstitial sulfides, interpreted as deposited from metasomatic fluids or melts. The enclosed sulfides typically are enriched in Os+Ir and relatively depleted in Pt+Pd, while the interstitial sulfides typically are depleted in Os+Ir and relatively enriched in Pd±Pt (Alard et al., 1999a).

Many, and perhaps most, mantle-derived peridotites contain a mixture of these two sulfide types, and this mixing has major implications for the interpretation of whole-rock PGE patterns. It also raises serious questions about the interpretation of whole-rock Re-Os data. Burton et al. (1999) recently demonstrated a significant difference in the γOs of interstitial and enclosed sulfides in a peridotite xenolith from Kilbourne Hole, suggesting that a meaningful depletion age of the rock may be given only by the sulfide inclusions. The sulfide phases are small, rare and difficult to separate, and separation obliterates the microstructural context needed for interpretation. We therefore have investigated the potential of laser-ablation analysis of Os isotopes in such sulfides, using multi-collector ICPMS techniques.

Analytical Methods

We have used a Merchantek LUV266 nm laser microprobe, attached to a Nu Plasma MC-ICPMS. The Nu instrument's geometry features a variable dispersion ion optics and a fixed collector array with 12 Faraday cups and 3 ETP ion counters. Two configurations have been tested thus far:

· Masses 194, 192, 191, 190, 189, 188, 187, 186 and 185 in Faraday cups

· Masses 194, 192, 191, 190, 189, 188, 186 in Faraday cups; 187, 185 in ion counters.

The overlap of 187Re on 187Os has been corrected using a 187Re/185Re = 1.674; this value was measured using an Os standard spiked with Re, and iteratively determining the value of 187Re/185Re required to produce the correct 187Os/188Os. Os-isotope data are normalised to 189Os/192Os=0.39593.

Table 1 shows the accuracy and precision of the instrument, as determined by replicate analysis of several Os standards. The accuracy of the Re correction is demonstrated by repeated analysis of a mixed Re-Os standard supplied by V. Bennett (ANU-RSES).

Table 1. Solution Analyses of Os Standards


Sample
conc.
192/188
190/188
189/188
187/188
186/188
ANU
95, 10 ppb
mean
3.08010
1.98340
1.21961
0.12344
0.11999
(n=8)
2sd
0.00029
0.00100
0.00040
0.00019
0.00014
ANU
10ppb
mean
3.08383
1.98422
1.21987
0.12344
0.12008
(n=7)
w/2.5 ppb Re
2sd
0.00041
0.00056
0.00048
0.00123
0.00025
nominal
0.1233
±0.0003
DTM
10 ppb
mean
3.08204
1.98547
1.21992
0.17363
0.12110
(n=12)
2sd
0.00586
0.00879
0.00171
0.00067
0.00123
nominal
0.17367
±0.00058
JMC Os
30 ppb
mean
3.07997
1.98333
1.21960
0.10690
0.12044
Bern
2sd
0.00045
0.00115
0.00044
0.00017
0.00014
(n=8)
nominal
0.10694
±0.00004

Note: All analyses carried out using Cetac MCN6000.

The LAM analysis has been tested using a synthetic NiS bead (PGE-A) with 200 ppm Os of known isotopic composition and PGE pattern. Repeated analyses of PGE-A (Table 2) showed the precision obtainable as a function of signal intensity, and demonstrated that under the conditions (small spot size, short runs) expected for real samples with similar levels of Os, we could achieve precision equivalent to N-TIMS results using the Faraday cup array; better precision was achieved using the combined Faraday plus ion-counting setup. In both cases, a major limitation on the precision obtainable is the need for a signal high enough and long enough to get a precise measurement of the mass bias in the Faraday cups.

Table 2. LAM-MC-ICPMS analyses of PGE-A

Run
integ. Time
Os, volts
187Os/188Os
±2se
setup
TIMS
0.10645
0.00002
10
3100 sec
0.90
0.10668
0.00008
Faraday
16
3100 sec
0.37
0.10742
0.00016
Faraday
11
3100 sec
0.35
0.10634
0.00019
Faraday
15
3100 sec
0.26
0.10654
0.00024
Faraday
26
3100 sec
0.21
0.10700
0.00036
Faraday
13
3100 sec
0.06
0.10682
0.00104
Faraday
49
60 sec
0.61
0.10661
0.00024
Faraday
47
60 sec
0.26
0.10668
0.00052
Faraday
48
60 sec
0.21
0.10604
0.00058
Faraday
49
30 sec
0.65
0.10671
0.00034
Faraday
47
30 sec
0.24
0.10649
0.00070
Faraday
48
30 sec
0.24
0.10630
0.00080
Faraday
Mean
0.10664
0.00071
Faraday
39
60 sec
0.25
0.10645
0.00036
Farad.+IC
40
60 sec
0.06
0.10608
0.00020
Farad.+IC
41
60 sec
0.06
0.10619
0.00019
Farad.+IC
42
60 sec
0.13
0.10660
0.00015
Farad.+IC
43
60 sec
0.12
0.10634
0.00028
Farad.+IC
44
60 sec
0.37
0.10712
0.00013
Farad.+IC
45
60 sec
0.28
0.10690
0.00012
Farad.+IC
46
60 sec
0.24
0.10679
0.00018
Farad.+IC
Mean
60 sec
0.10656
0.00072
Farad.+IC
39
30 sec
0.29
0.10693
0.00048
Farad.+IC
40
30 sec
0.05
0.10587
0.00022
Farad.+IC
41
30 sec
0.05
0.10627
0.00040
Farad.+IC
42
30 sec
0.12
0.10646
0.00024
Farad.+IC
43
30 sec
0.10
0.10622
0.00038
Farad.+IC
44
30 sec
0.37
0.10742
0.00036
Farad.+IC
45
30 sec
0.27
0.10675
0.00034
Farad.+IC
46
30 sec
0.20
0.10628
0.00038
Farad.+IC
Mean
30 sec
0.10653
0.00098
Farad.+IC

Results

The results of sample analyses are shown in Table 3 and Figure 1.

Table 3. Results on Natural Samples


Sample
Os (Volts)
Os 

(ppm)

F/IC
187Os/188Os
±2se
187Re/188Os
TRD

(Ga)

±2SE
TMA

(Ga)

±2SE
Siberia: Udachnaya
UD-Ol61
1.5
849*
F
0.1072
0.0002
0.006
3.02
0.03
3.07
0.05
UD-Ol60
0.06
~54
IC
0.1033
0.0012
0.005
3.58
0.18
3.62
0.20
UD-Ol67
0.01
~19
IC
0.1118
0.0013
0.081
2.35
0.19
2.94
0.51
* data obtained by LAM-ICP-MS, confirmed by Proton probe analyses (Os=892±126 ppm, Pers. Comm. C.G. Ryan) 
Slave Craton
A-481
0.22
~250
IC
0.1080
0.0001
0.014
2.91
0.02
3.01
0.07
Kaapvaal Craton
FS-8 1
0.04
~23*
F
0.12146
0.0064
0.180
0.92
0.96
1.66
1.70
FS-8 2
0.02
~15*
IC
0.12494
0.0010
0.120
0.39
0.15
0.56
0.27
*LAM-ICP-MS analyses of FS-8 sulfides range from 7 to 20 ppm 
Mt Gambier, S. Australia, Newer Volcanics Province.
GAM-VL9-S3
0.02
42*
IC
0.1160
0.0007
0.016
1.74
0.10
1.81
0.14
GAM-VL9-S1(b)
0.002
1*
IC
0.1283
0.0056
0.035
-0.11
0.85
-0.13
0.91
Mt Quincan, (Queensland), Atherton Volcanic Province.
94-A4-37 S1
0.02
~45
IC
0.1242
0.0006
0.045
0.51
0.09
0.57
0.13
94-A4-39-s1(a)
0.003
~4
IC
0.1203
0.0036
0.071
1.09
0.54
1.32
0.72
Montferrier, Southern France (Montpellier Volcanic Province)
Pg2-s3a
0.004
~2
IC
0.1444
0.0028
0.041
-2.6
0.4
-2.9
0.3
PG 5 s1
0.039
~20
IC
0.1291
0.0007
0.222
-0.2
0.1
-0.5
0.0
MTF 37a
0.036
~20
F
0.1337
0.0034
0.060
-1.0
0.5
-1.1
0.5
MTF 37a
0.047
~20
F
0.1316
0.0042
0.072
-0.6
0.6
-0.8
0.6
MTF 37a
0.034
~20
F
0.1209
0.0052
0.211
1.0
0.8
2.1
1.7
MtfF37C S3
0.009
~5
IC
0.1361
0.0014
0.137
-1.3
0.2
-2.0
-0.1

Siberia: Three sulfide inclusions in olivine macrocrysts, derived from disaggregated peridotites, were extracted from the Udachnaya kimberlite. The largest has been analysed for PGEs by LAM-ICPMS and has a high Os content. A similar precision and a slightly older age was obtained from a smaller inclusion. A still smaller inclusion gave a lower precision, and a slightly younger model age. These ages are equivalent to the oldest TRD model ages (3.0-3.2 Ga) obtained by conventional analysis of Udachnaya peridotite xenoliths (Pearson et al., 1995), and to the TRD ages of three sulfide inclusions in a single diamond from the same pipe (3.1±0.3 to 3.5±0.3 Ga; Pearson et al., 1999).

Slave Craton: A single large inclusion was extracted from an olivine macrocryst (Fo92) in the A841 kimberlite, Lac de Gras area. The age is equivalent to the oldest model ages (3.1 Ga) reported from xenoliths in the Jericho pipe to the north (Irvine et al., 1999), and is significantly older than the oldest event (2.8 Ga) known in the overlying crust (Yamashita and Creaser, 1999).

Kaapvaal Craton: Two interstitial sulfide inclusions in a high-temperature sheared peridotite xenolith from the Frank Smith kimberlite (100 Ma eruption age) have modest Os contents and the analyses have relatively low precision. Other geochemical evidence (low mg#, high Ti, Zr and Y in rims of strongly zoned garnets; Griffin et al., 1989) suggests that these xenoliths were infiltrated by asthenosphere-derived melts close to the time of kimberlite eruption, and the sulfides may contain a mixture of old and young components.

Mt. Gambier (S. Australia): One enclosed sulfide gives a model age of 1.8±0.17 Ga, within uncertainty of the 1.9 Ga model ages obtained by Handler (1998) on three xenoliths from the same volcano. An interstitial sulfide from another Mt. Gambier xenolith gives a 187Os/188Os equivalent to modern asthenosphere, emphasising that these rocks contain two generations of sulfides with different Os isotopic compositions and Re/Os ratios.

Massif Central: MTF37 is a fertile peridotite (Al2O3 = 4.2%) with a PUM-like 187Os/188Os (0.1277; Meisel et al., 1996,1999). Individual grains of interstitial sulfides have 187Os/188Os ranging from 0.1209-0.1361, indicating that the whole-rock 187Os/188Os value represents mixing of older (Proterozoic) components and younger radiogenic Os (Fig. 1). Other xenoliths from the same locality display even more radiogenic Os, with 187Os/188Os = 0.1272-0.1756, and whole-rock 187Os/188Os increases with S contents and Pd/Ir ratios (Alard et al., 1999). The data suggest that the radiogenic Os and high S contents of this suite reflect high proportions of metasomatic sulfide. This sample suite highlights the resolving power of the in-situ analytical method.
 
 

Figure 1. Os data for Montferrier sulfides. PUM and CI values from Meisel et al. (1996).

Conclusions

LAM-MC-ICPMS analysis of single sulfide inclusions in mantle-derived peridotites can give useful data on the Re-Os system. With current technology, the best results are obtained from sulfide grains 350 µm in diameter, with Os contents 340 ppm. The data obtained from such grains have a precision equivalent to N-TIMS data on whole rock samples. The lower-precision data obtained on low-Os interstitial sulfides are still very useful in understanding the movement of Os within the lithosphere, as demonstrated by the Massif Central samples. The major advantage of the LAM-MC-ICPMS technique lies in providing in-situ analyses with the complete spatial and microstructural context required for meaningful interpretation.

Further improvements in sensitivity are expected to increase precision; these developments include multicycle ion-counter analysis, a shield torch for the MC-ICPMS and the use of He as the carrier gas in the ablation cell.

The application of this technique will make Re-Os analysis of mantle-derived rocks more rapid and less expensive, and enable larger-scale surveys of the age structure of the continental lithosphere.

References

Alard, O., Griffin, W.L., Lorand, J.P., Jackson, S.E. and O'Reilly, S.Y. 1999a. Highly siderophile elements in mantle sulfide by LAM-ICPMS. Ext. Abst. 9th Goldschmidt Conf.

Alard O., Lorand, J.P., Reisberg, L., Bodinier. J.L., O'Reilly, S.Y., Griffin, W.L. and Dautria, J.M., 1999b. Sulfur enrichment via volatile-richmetasomatism in montferrier xenoliths (Southern France): Consequences for highly siderophile elements in the sub-continental mantle. Orogenic lherzolites and mantle processes conference, Pavia, Italia, Ofioliti, 24:48-49.

Burton, K.W., Chiano, P, Birck, J.-L. and Allegre, C.J. 1999. Osmium isotope disequilibrium between minerals in a spinel-lherzolite. Earth Planet. Sci. Lett. 172, 311-322.

Griffin, W.L., Smith, D., Boyd, F.R., Cousens, D.R., Ryan, C.G., Sie, S.H. and Suter, G.F. 1989: Trace element zoning in garnets from sheared mantle xenoliths. Geochim. Cosmochim. Acta. 53, 561-567

Handler, M.R., Bennett, V.C. and Esat, T.M. 1997. The persistence of off-cratonic lithospheric mantle: Os isotopic systematics of variably metasomatised southeast Australian xenoliths. Earth Planet. Sci. Lett. 151, 61-75.

Irvine, G.J., Kopylova, M.G., Carlson, R.W., Pearson, D.G.,Shirey, S.H. and Kjarsgaard, B.A. 1999. Age of the lithospheric mantle beneath and around the Slave Craton: A Re-Os isotope study of peridotite xenoliths from the Jericho and Somerset Island kimberlites. Ext. Abst. 9th Goldschmidt Conf.

Meisel, T., Walker, R.J., Morgan, J.W., 1996, The osmium isotopic composition of the Earth's primitive mantle. Nature, 383, 517-520.

Meisel T., Walker R.J., Irving, A.J. and Lorand J.-P. 1999: Osmium isotopic composition of mantle xenoliths: a global perspective (submitted to GCA).

Pearson, D.G., Shirey, S.B., Carlson, R.W., Boyd, F.R., Pokhilenko, N.P. and Shimizu, N. 1995. Re-Os, Sm-Nd and Rb-Sr isotope evidence for thick Archaean lithospheric mantle beneath the Siberian craton modified by multistage metasomatism. Geochim. Cosmochim. Acta 59, 959-977.

Pearson, D.G., Shirey, S.B., Bulanova, G.P., Carlson, R.W. and Milledge, H.J. 1999. Re-Os isotope measurements of single sulfide inclusions in a Siberian diamond and its nitrogen aggregation systematics. Geochim. Cosmochim. Acta 63, 703-711.

Yamashita, K. and Creaser, R.A. 1999. Geochemical and isotopic constraints for the evolution of the Slave Province. Ext. Abst. 9th Goldschmidt Conf.