Skip Navigation


Journal of Petrology Advance Access originally published online on November 29, 2005
Journal of Petrology 2006 47(4):673-703; doi:10.1093/petrology/egi089
This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrow Supplementary data
Right arrow All Versions of this Article:
47/4/673    most recent
egi089v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrow Search for citing articles in:
ISI Web of Science (41)
Right arrowRequest Permissions
Google Scholar
Right arrow Articles by BECKER, M.
Right arrow Articles by ROEX, A. P. L.
Right arrow Search for Related Content
GeoRef
Right arrow GeoRef Citation
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?

© The Author 2005. Published by Oxford University Press. All rights reserved. For Permissions, please e-mail: journals.permissions@oxfordjournals.org

Geochemistry of South African On- and Off-craton, Group I and Group II Kimberlites: Petrogenesis and Source Region Evolution

MEGAN BECKER and ANTON P. LE ROEX*

DEPARTMENT OF GEOLOGICAL SCIENCES, UNIVERSITY OF CAPE TOWN, RONDEBOSCH, 7701, SOUTH AFRICA

RECEIVED MAY 17, 2005; ACCEPTED NOVEMBER 2, 2005


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 GEOLOGICAL SETTING AND SAMPLING
 ANALYTICAL TECHNIQUES
 PETROGRAPHY
 GEOCHEMISTRY
 PETROGENESIS
 DISCUSSION
 CONCLUSIONS
 SUPPLEMENTARY DATA
 REFERENCES
 
Bulk-rock geochemical compositions of hypabyssal kimberlites, emplaced through the Archaean Kaapvaal craton and Proterozoic Namaqua–Natal belt, are used to estimate close-to-primary magma compositions of Group I kimberlites (Mg-number = 0·82–0·87; ~22–28 wt % MgO; ~21–30 wt % SiO2; ~10–17 wt % CaO; ~0·2–1·7 wt % K2O) and Group II kimberlites (Mg-number = 0·86–0·89; ~23–29 wt % MgO; ~28–36 wt % SiO2; 8–13 wt % CaO; ~1·6–4·6 wt % K2O). Group I kimberlites are distinguished from Group II by their lower Ba/Nb (<12), Th/Nb (<1·1) and La/Nb (<1·1) but higher Ce/Pb (>22) ratios. The distinct rare earth element patterns of the two types of kimberlites indicate a more highly metasomatized source for Group II kimberlites, with more residual clinopyroxene and less residual garnet. The similarity of Sr and Nd isotope ratios and diagnostic trace element ratios (Ce/Pb, Nb/U, La/Nb, Ba/Nb, Th/Nb) of Group I kimberlites to ocean island basalts (OIB), but more refractory Mg-numbers and Ni contents, are consistent with derivation of Group I kimberlites from subcontinental lithospheric mantle (SCLM) that has been enriched by OIB-like melts or fluids. Source enrichment ages and plate reconstructions support a direct association of these melts or fluids with Mesozoic upwelling beneath southern Africa of a mantle plume(s), at present located beneath the southern South Atlantic Ocean. In contrast, the geochemical characteristics of both on- and off-craton Group II kimberlites show strong similarity to calc-alkaline magmas, particularly in their Nb and Ta depletion and Pb enrichment. It is suggested that Group II kimberlites are derived from both Archaean and Proterozoic lithospheric mantle source regions metasomatized by melts or fluids associated with ancient subduction events, unrelated to mantle plume upwelling. The upwelling of mantle plumes beneath southern Africa during the Mesozoic, at the time of Gondwana break-up, may have acted as a heat source for partial melting of the SCLM and the generation of both Group I and Group II kimberlite magmas.

KEY WORDS: kimberlite; geochemistry; petrogenesis; mantle plumes; South Africa


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 GEOLOGICAL SETTING AND SAMPLING
 ANALYTICAL TECHNIQUES
 PETROGRAPHY
 GEOCHEMISTRY
 PETROGENESIS
 DISCUSSION
 CONCLUSIONS
 SUPPLEMENTARY DATA
 REFERENCES
 
Studies of kimberlite magmatism provide a unique opportunity to investigate partial melting processes in the deeper parts of the Earth's upper mantle. Their unique ultramafic, potassic and volatile-rich character, with enrichment in both compatible and incompatible trace elements, suggests derivation from mantle source regions with a complex evolutionary history. Although there have been numerous studies focused on characterizing the major element geochemistry (e.g. Clement, 1982Go; Shee, 1985Go), trace element geochemistry (e.g. Wedepohl & Muramatsu, 1975Go; Smith et al., 1985bGo) or isotope geochemistry (e.g. Smith, 1983Go; Fraser et al., 1985–1986Go; Nowell et al., 2004Go), only relatively recently has trace element geochemistry been used to constrain petrogenetic processes and the evolution of kimberlite source regions (e.g. Fraser & Hawkesworth, 1992Go; Tainton & McKenzie, 1994Go; Beard et al., 1998Go; le Roex et al., 2003Go; Chalapathi Rao et al., 2004Go; Coe, 2004Go; Harris et al., 2004Go).

Semi-quantitative genetic models argue for derivation of kimberlite by very low degrees of partial melting of metasomatically enriched, previously melt-depleted, garnet peridotite facies lithospheric mantle (Tainton & McKenzie, 1994Go; le Roex et al., 2003Go; Harris et al., 2004Go). Such models are consistent with experimental studies suggesting carbonated garnet peridotite source regions (e.g. Dalton & Presnall, 1998aGo; Ulmer & Sweeney, 2002Go). More exotic source regions such as garnetite (Edgar & Charbonneau, 1993Go) and clinopyroxene and mica-rich veins (Foley, 1992bGo) have, however, also been suggested.

Two broad varieties of southern African kimberlites are recognized, and have been classified variously as ‘basaltic’ and ‘micaceous’ (on petrographic grounds; Wagner, 1914Go), ‘Group I’ and ‘Group II’ (on isotopic grounds; Smith, 1983Go) or ‘kimberlite’ and ‘orangeite’ (Mitchell, 1995Go). Group I kimberlites are defined as being slightly less radiogenic in terms of their Sr isotope and more radiogenic in terms of their Nd isotope compositions than present-day Bulk Earth, and in this regard show similarity to ocean island basalts (OIB). Smith (1983)Go associated such kimberlites with a sub-lithospheric mantle source. Group II kimberlites, which tend to be more micaceous, have highly radiogenic Sr isotope compositions relative to Bulk Earth and have been associated with direct derivation from the subcontinental lithospheric mantle (SCLM; Smith, 1983Go; Fraser et al., 1985–1986Go; Tainton & McKenzie, 1994Go; Coe, 2004Go). The interpretation of the results of Hf isotope studies are variable; Nowell et al. (2004)Go have proposed that both Group I and Group II kimberlites have sub-lithospheric sources, whereas Griffin et al. (2000)Go and Choukroun et al. (2005)Go preferred a model with contributions from both asthenospheric and sub-lithospheric mantle sources. Similarly, le Roex et al. (2003)Go have argued for a two-stage model for the generation of Group I kimberlites, involving both sub-lithospheric and lithospheric mantle source components.

The origin of both Group I and Group II Mesozoic kimberlites within southern Africa, from both on-craton and off-craton settings, is investigated in this study. This is done by combining major and trace element and isotope and subsequently developing semi-quantitative genetic models. These models are aimed at identifying petrogenetic processes and source region compositions that can account for the observed similarities and differences. Possible causes of kimberlite source region enrichment (metasomatism), and the origin of metasomatic melts or fluids, are also explored.

To achieve these objectives, close-to-primary magma compositions for the individual kimberlites analysed in this study are determined, after taking into account the effects of various secondary processes and fractional crystallization. Semi-quantitative modelling of kimberlite genesis and source geochemistry is used to evaluate mantle source region characteristics and evolutionary processes, as a basis for the development of an integrated model describing the petrogenesis of South African Group I and Group II Mesozoic kimberlites.


    GEOLOGICAL SETTING AND SAMPLING
 TOP
 ABSTRACT
 INTRODUCTION
 GEOLOGICAL SETTING AND SAMPLING
 ANALYTICAL TECHNIQUES
 PETROGRAPHY
 GEOCHEMISTRY
 PETROGENESIS
 DISCUSSION
 CONCLUSIONS
 SUPPLEMENTARY DATA
 REFERENCES
 
Kimberlites in southern Africa have been emplaced either through the Archaean Kaapvaal craton (de Wit et al., 1992Go; ‘on-craton kimberlites’) or through the surrounding Namaqua–Natal Mesoproterozoic metamorphic belt (off-craton kimberlites). Most of the 23 Group I and Group II kimberlites that form the basis of this study crop out within the shales, mudstones and sandstones of the Phanerozoic Karoo Supergroup, and range in age from ~74 Ma to ~150 Ma (Davis, 1977Go; Smith et al., 1985aGo, 1994Go); the Group II kimberlites are predominantly older than the Group I kimberlites.

On-craton kimberlites analysed in this study (Fig. 1; Table 1) are the Group I kimberlite occurrences at Koffiefontein, Jagersfontein and Monastery, and the Bellsbank, Newlands, New Elands and Finsch Group II kimberlites. Off-craton kimberlites comprise the Abbotsford East and Zeekoegat Group I kimberlites from East Griqualand (Nixon et al., 1983Go), as well as Group I and Group II kimberlites located within the vicinity of Victoria West, Prieska and Hanover (Robey, 1981Go; Skinner et al., 1992Go). Some of the kimberlites located near Prieska and Hanover are inferred to have been emplaced within Archaean crystalline basement (e.g. Goedehoop, Andriesfontein kimberlites) according to their relationship to the Brakbos fault, the generally accepted border between Proterozoic and Archaean Provinces (Pretorious, 1974Go; Cornell et al., 1986Go; Schmitz et al., 2004Go). In this study these kimberlites are considered to be on-craton.


Figure 1
View larger version (26K):
[in this window]
[in a new window]
 
Fig. 1. Location of on- and off-craton Group I and Group II kimberlites analysed in this study relative to the Kaapvaal craton (shaded in grey). Light dot–dash boundary represents the Namaqua–Natal mobile belt. Figure adapted from Carlson et al. (2000)Go. Filled symbols indicate locations of the Kimberley Group I kimberlites, the off-craton Uintjiesberg Group I kimberlite and the Star and Swartruggens Group II kimberlites. Numbers represent clusters of kimberlites; for individual kimberlites localities and coordinates are given in Table 1.

 

View this table:
[in this window]
[in a new window]
 
Table 1: Summary table of kimberlite coordinates and locations as illustrated in Fig. 1

 
Sampling was focused on hypabyssal kimberlite to avoid the excessive alteration and contamination present in most diatreme-facies kimberlite; material was selected to avoid obvious alteration, weathering and crustal contaminants. Where possible, multiple samples were collected from individual kimberlite bodies; the final sample suite includes material recovered from field outcrop, underground mine workings, collections housed at the University of Cape Town and collections provided by the De Beers Geoscience Centre.


    ANALYTICAL TECHNIQUES
 TOP
 ABSTRACT
 INTRODUCTION
 GEOLOGICAL SETTING AND SAMPLING
 ANALYTICAL TECHNIQUES
 PETROGRAPHY
 GEOCHEMISTRY
 PETROGENESIS
 DISCUSSION
 CONCLUSIONS
 SUPPLEMENTARY DATA
 REFERENCES
 
All samples were prepared for geochemical analysis at the University of Cape Town (UCT). Weathered surfaces were removed using a hydraulic splitter, after which samples were split and passed through a jaw crusher. The jaw crush chips were individually hand picked, with careful attention being paid to selecting kimberlite fragments free from weathering, crustal and mantle xenoliths, megacrysts and calcite veining. The clean samples were then powdered using a Sieb swing mill and carbon steel vessel. A subset of kimberlites was received from the De Beers Geoscience Centre as powders (K16/6; K51/2; K19/15; K5/P2; K6/11; K2/2 and K2/13) and, therefore, could not be prepared as desired.

Major element and Ni and Cr concentrations were determined by X-ray fluorescence (XRF) on a Philips X'Unique wavelength-dispersive spectrometer using a low dilution fusion technique (Willis, 1999Go). Selected trace elements (Co, V, Cu, Nb and Zr) were determined by XRF on 6 g powder briquettes with accuracy and precision values for XRF analyses similar to those reported by le Roex et al. (1981)Go. A more comprehensive set of trace elements including the rare earth elements (REE) was analysed by solution inductively coupled plasma mass spectrometry (ICP-MS) on a Perkin–Elmer ELAN 6000 system. A detailed description of the sample preparation techniques and analyses of procedural blanks and standards has been given by le Roex et al. (2001Go, 2003Go). The relative standard deviation for most elements was generally better than 3 %. Zr and Nb data by XRF and ICP-MS generally show good agreement, suggesting that incomplete dissolution of Fe–Ti-oxides or zircon was not a significant problem. CO2 was determined using a Karbonat Bombe according to the method described by Birch (1981)Go, with relative precision between 2 and 5 %.

Representative samples were prepared for Sr and Nd isotope analysis by thermal ionization mass spectrometry (TIMS) on a VG Sector seven-collector system run in multi-dynamic mode, using standard ion exchange techniques as described by Hart (1977)Go and Zindler (1980)Go. 87Sr/86Sr isotope ratios were normalized against a 86Sr/88Sr ratio of 0·1194 and 143Nd/144Nd against a 146Nd/144Nd ratio of 0·7219. International standards NBS 987 and La Jolla yielded 87Sr/86Sr and 143Nd/144Nd ratios of 0·710260 ± 17 (2{sigma}; n = 3) and 0·511829 ± 10 (2{sigma}; n = 1), respectively.


    PETROGRAPHY
 TOP
 ABSTRACT
 INTRODUCTION
 GEOLOGICAL SETTING AND SAMPLING
 ANALYTICAL TECHNIQUES
 PETROGRAPHY
 GEOCHEMISTRY
 PETROGENESIS
 DISCUSSION
 CONCLUSIONS
 SUPPLEMENTARY DATA
 REFERENCES
 
The analysed kimberlites can be divided according to texture into macrocrystic and aphanitic varieties (Clement & Skinner, 1979Go; Clement et al., 1984Go). The non-genetic textural term ‘macrocrystic’ is used to indicate the presence of predominantly large, anhedral, olivine crystals (>2 mm) that may have been derived from disaggregated mantle xenoliths (Clement et al., 1984Go; le Roex et al., 2003Go). All samples with >10 vol. % macrocrysts are classified as macrocrystic. More euhedral olivine and phlogopite crystals, >1 mm in size, are referred to as phenocrysts and if <1 mm, as microphenocrysts. Mineral abundances are summarized for representative kimberlite samples in Table 1.

Both on- and off-craton Group I kimberlites are mineralogically classified as serpentine calcite or calcite kimberlites (Table 1; using the classification scheme proposed by Clement & Skinner, 1979Go); most are sparsely macrocrystic to macrocrystic in texture, except for the aphanitic Lekkerfontein kimberlite. The subhedral to anhedral olivine macrocrysts, which may be up to 10 mm in size, are variably serpentinized and sometimes chloritized, and are the most abundant macrocryst phase (10–33 vol. %), whereas phlogopite (<10 vol. %) and ilmenite (<5 vol. %) are less abundant. Phlogopite macrocrysts tend to be corroded and show undulose extinction or kink banding, as well as, in some samples (Goedehoop, Klipgatsfontein, Hebron kimberlites), the development of calcite stringers along cleavage planes. Rare ilmenite macrocrysts, characterized by groundmass perovskite aggregations along their rims, are possibly correlated with the kimberlite ‘megacryst suite’, including zircon, known to occur in these kimberlites (e.g. Monastery, Lekkerfontein and Gansfontein kimberlites; Boyd & Nixon, 1975Go; Monastery, Lekkerfontein and Gansfontein kimberlites; Robey, 1981Go).

Euhedral to subhedral olivine phenocrysts and microphenocrysts are present in all the analysed Group I kimberlites, varying in abundance between 10 and 30 vol. %. In contrast, phlogopite phenocrysts occur in only two of the samples (JJG 4295, LKF 1). Fine-grained calcite and serpentine are the predominant groundmass phases, forming either a uniform or segregationary texture (Clement & Skinner, 1979Go), with some of the segregations showing very well-developed rhombohedral calcite crystals up to 500 µm in size. Calcite segregations can reach up to a few millimetres in size. Scattered, stubby-shaped groundmass phlogopite crystals occur in most kimberlites except for Goedehoop, Hebron and Kligatsfontein, whereas equant groundmass opaque oxide and lesser amounts of perovskite crystals (~100 µm) are ubiquitous. In comparison with other on-craton kimberlites (e.g. Kimberley kimberlites, Clement, 1982Go), the off-craton kimberlites in this study appear not to be characterized by the presence of groundmass monticellite. However, given that primary monticellite is prone to alteration, forming secondary serpentine and calcite (Clement, 1982Go), its paucity may reflect the more altered nature of the off-craton kimberlites.

Both on- and off-craton Group II kimberlites are classified as phlogopite kimberlites or calcite phlogopite/phlogopite calcite kimberlites and vary in texture between sparsely macrocrystic and macrocrystic (up to 45 vol. % macrocrysts). Anhedral olivine macrocrysts, up to 15 mm in size, are variably altered, by serpentinization along fractures, calcitization in patches or chloritization along their rims. Strained phlogopite macrocrysts (up to 20 vol. %) show strong pleochroism and are rimmed by tetraferriphlogopite or chlorite (e.g. Newlands kimberlite). In the off-craton Eendekuil kimberlite, the phlogopite macrocrysts show slight flow alignment. Garnet macrocrysts, rimmed by kelyphite or chlorite, were also identified in a few of the kimberlite samples (Newlands: JJG 24, KN 2; Markt: MRK 3).

Olivine phenocrysts are ubiquitous in all of the analysed Group II kimberlites and may constitute up to 35 vol. % of the rock. Tabular phlogopite phenocrysts, rimmed by tetraferriphlogopite, are also common, whereas subdhedral diopside phenocrysts were identified only in Eendekuil kimberlite sample EKL 1. Short, stubby, interlocking phlogopite crystals, rimmed by tetraferriphlogopite, form the dominant (15–40 vol. %) groundmass phase in most of the analysed Group II kimberlites; these vary in size from being very fine-grained (50–100 µm) to fine-grained (100–500 µm). Groundmass calcite is the next most common matrix phase (up to 30 vol. %) and is generally interstitial, although it may sometimes form calcite segregations up to 4 mm in size. Although matrix serpentine was recognized in a few kimberlite samples, it is likely that some is secondary in origin (e.g. Brandewynskuil kimberlite). Similarly, the chloritization of groundmass phlogopite in some kimberlites (Bellsbank and Newlands) and pervasive carbonatization in others (Roberts Victor, RVK 1) is inferred to be an alteration feature. Accessory opaque oxides and perovskite crystals are very fine-grained (20–50 µm) and relatively sparse; perovskite crystals may be poikilitically enclosed in phlogopite microphenocrysts (e.g. New Elands kimberlite). Trace amounts of apatite needles, 50–500 µm in size, have also been recognized in the Bellsbank, New Elands, Newlands, Roberts Victor and Markt kimberlites. No significant petrographic differences are recognized between on- and off-craton Group II kimberlites.


    GEOCHEMISTRY
 TOP
 ABSTRACT
 INTRODUCTION
 GEOLOGICAL SETTING AND SAMPLING
 ANALYTICAL TECHNIQUES
 PETROGRAPHY
 GEOCHEMISTRY
 PETROGENESIS
 DISCUSSION
 CONCLUSIONS
 SUPPLEMENTARY DATA
 REFERENCES
 
The major and trace element compositions of selected on- and off-craton kimberlites are reported in Table 2, with Nd–Sr isotope compositions in Table 3. The complete kimberlite dataset used in this study is available online (Electronic Appendix 1; http://www.petrology.oupjournals.org). The data are supplemented with analyses from the on-craton Group I Kimberley kimberlites (le Roex et al., 2003Go), Swartruggens and Star Group II kimberlites (Coe, 2004Go) shown as fields, and with analyses from the off-craton Group I Uintjiesberg kimberlite (Harris et al., 2004Go), shown as data points in Figs 26.


Figure 2
View larger version (32K):
[in this window]
[in a new window]
 
Fig. 2. Variation of (a) SiO2 vs MgO, (b) FeO* vs MgO, (c) TiO2 vs K2O and (d) CaO vs SiO2 in on- and off-craton Group I and Group II kimberlites. Symbols for off-craton Group I kimberlites include analyses from the Uintjiesberg kimberlite (Harris et al., 2004Go). Dashed field represents the evolved Group II Eendekuil kimberlite. Light shaded field represents the on-craton Group I Kimberley kimberlites (le Roex et al., 2003Go) and dark shaded field represents the on-craton Swartruggens and Star Group II kimberlites (Coe, 2004Go). Group I–Group II dividing lines in (c) are from Smith et al. (1985b)Go.

 

Figure 3
View larger version (35K):
[in this window]
[in a new window]
 
Fig. 3. Variation of (a) Ni vs MgO, (b) Rb vs La, (c) Nb vs La, (d) Pb vs Ce, (e) Ba vs Nb and (f) Th vs La in on- and off-craton Group I and Group II kimberlites. Dividing lines for Ce/Pb in (d) and Ba/Nb in (e) provide distinction between kimberlite groups. Symbols for off-craton Group I kimberlites include analyses from the Uintjiesberg kimberlite (Harris et al., 2004Go). Light shaded field represents the on-craton Group I Kimberley kimberlites (le Roex et al., 2003Go) and dark shaded field represents the on-craton Swartruggens and Star Group II kimberlites (Coe, 2004Go).

 

Figure 4
View larger version (28K):
[in this window]
[in a new window]
 
Fig. 4. Variation of (a) Ce/Pb vs Ba/Nb, (b) Ba/Nb vs La/Nb, (c) La/Sm vs La/Yb and (d) Gd/Yb vs La/Yb in on- and off-craton Group I and Group II kimberlites. Dividing lines for Ba/Nb, Ce/Pb and La/Nb in (a) and (b) provide distinction between kimberlite groups. Dashed field represents the Group IIb kimberlites. Light shaded field represents the on-craton Group I Kimberley kimberlites (le Roex et al., 2003Go) and dark shaded field represents the on-craton Swartruggens and Star Group II kimberlites (Coe, 2004Go).

 

Figure 5
View larger version (23K):
[in this window]
[in a new window]
 
Fig. 5. Chondrite-normalized REE diagrams for on- and off-craton (a) Group I and (b) Group II kimberlites. Data for labelled fields are from le Roex et al. (2003)Go, Coe (2004)Go and Harris et al. (2004)Go. Normalizing values are from Sun & McDonough (1989)Go.

 

Figure 6
View larger version (35K):
[in this window]
[in a new window]
 
Fig. 6. Primitive mantle-normalized trace element patterns for on- and off-craton (a) Group I and (b) Group II kimberlites, with highlighted trace element anomalies. Data for labelled fields are from le Roex et al. (2003)Go, Coe (2004)Go and Harris et al. (2004)Go. Normalizing values are from Sun & McDonough (1989)Go.

 

View this table:
[in this window]
[in a new window]
 
Table 2: Bulk-rock geochemical analyses of on- and off-craton Group I and Group II kimberlites, analysed by XRF and ICP-MS

 

View this table:
[in this window]
[in a new window]
 
Table 3: Measured Sr and Nd isotope ratios and calculated parent/daughter ratios for on- and off-craton, Group I and Group II kimberlites

 
On- and off-craton Group I kimberlites
SiO2 and MgO concentrations of Group I kimberlites show broad positive correlations (Fig. 2a) with SiO2 contents varying from 22·5 to 33·3 wt %, whereas MgO varies between 18·6 and 33·3 wt %. On the whole, the more macrocrystic samples tend to have higher SiO2 and MgO contents, and Mg-numbers (atomic Mg/(Mg + Fe2+) >0·85), than the less macrocrystic samples (Mg-numbers = 0·76–0·85). FeO* (all iron as Fe2+) concentrations vary from 7·5 to 15·3 wt % and tend to show a broad negative correlation with MgO (Fig. 2b). Samples from the Monastery, Lekkerfontein and Gansfontein kimberlites that have >11 wt % FeO* are interpreted to have been contaminated by ilmenite megacrysts, known to be abundant at these localities (Boyd & Nixon, 1975Go; Robey, 1981Go). It is also apparent that off-craton kimberlites tend to have lower SiO2 and MgO, but higher FeO* than on-craton kimberlites (disregarding those kimberlites affected by ilmenite megacryst contamination), although still broadly falling within the field of on-craton Group I kimberlites.

The range in TiO2 contents of the samples is greater than that for Al2O3 (1·0–5·1 wt % TiO2 and 1·8–4·4 wt % Al2O3; Fig. 2c, Table 2), and samples fall generally within the field of the on-craton Kimberley kimberlites. Ignoring the possible effect of ilmenite megacryst contamination (e.g. TiO2 >4 wt % in Gansfontein, Monastery and Lekkerfontein kimberlites), it is apparent that off-craton kimberlites have compositions that tend to plot towards the high TiO2 end of the field of on-craton Group I kimberlites (Fig. 2c). The majority of the Group I kimberlite samples have low K2O abundances (<1·5 wt %; Fig. 2c), except for a few outliers (1·6–2·4 wt % K2O; Lekkerfontein, Monastery, Gansfontein kimberlites). A well-defined negative correlation exists between CaO and SiO2 (Fig. 2d), whereas CO2 (1·5–11·7 wt %) and CaO (7·0–16·7 wt %) are positively correlated (Table 2). Most off-craton kimberlites tend to have higher CaO and CO2 than the on-craton kimberlites, although still broadly falling within the field of on-craton kimberlites.

Transition metals show good positive correlations with MgO (e.g. Ni in Fig. 3a), with Ni varying between 340 and 1500 ppm, and Cr varying between 450 and 2028 ppm (Table 2). No significant differences between on- and off-craton kimberlites are apparent within this group of elements, although concentrations do show some correlation with olivine macrocryst content. Large ion lithophile element (LILE) concentrations are highly variable (e.g. 9–112 ppm Rb; 280–2276 ppm Ba) and show no correlation with the more immobile incompatible elements (e.g. La in Fig. 3b). Inter-element correlations of high field strength elements (HFSE; Zr, Nb and Hf ) with other incompatible elements (e.g. La) are generally good over a wide range of concentrations (Fig. 3c). Inter-element ratios involving fluid-mobile elements are highly variable (e.g. K/Rb = 33–240; Ce/Pb = 14–49), whereas those involving the more immobile elements are more restricted (Figs 3 and 4; e.g. La/Nb = 0·7 ± 0·1; Th/Nb = 0·09 ± 0·02; Zr/Hf = 52 ± 6·6; Nb/Ta = 20 ± 5·6).

All of the analysed Group I kimberlites, both on- and off-craton, have broadly similar sub-parallel chondrite-normalized REE patterns (Fig. 5a), with considerable enrichment of the light REE (LREE; La = 300–1300 x chondrite) relative to the heavy REE (HREE; Lu = 3–8 x chondrite). In detail, La/YbN (44·7–267) and La/SmN (4·3–7·2) ratios are significantly variable. The analysed on- and off-craton kimberlites generally fall within the fields of the on-craton Kimberley and off-craton Uintjiesberg kimberlites, respectively (Fig. 5a), although it is evident from Fig. 4c and d that off-craton Group I kimberlites tend to have lower La/Sm, La/Yb and Gd/Yb ratios than some of their on-craton counterparts.

Primitive mantle-normalized trace element patterns are generally smooth, with strong enrichment in the more incompatible elements (e.g. Nb = 155–603 and Ba = 24–326 x primitive mantle), and relative depletions of Rb, K, Sr, Hf and Ti (Fig. 6a). These negative anomalies appear to be a characteristic feature of Group I kimberlites (le Roex et al., 2003Go; Harris et al., 2004Go). The magnitudes of the anomalies can be described by the X/X* notation (where X is the normalized value and X* the interpolated value calculated from the adjacent elements). Negative Rb and K anomalies tend to be larger (K/K* = 0·0–0·3; calculated using Th and La and not U because of its mobility) than the Ti anomalies (Ti/Ti* = 0·3–1·0). In some instances (e.g. Koffiefontein kimberlite), the magnitudes of the anomalies increase with the level of incompatible element enrichment in individual samples. Minor negative Sr anomalies, as well as occasional negative Hf anomalies, are also features of the analysed Group I kimberlites. Rare, strongly positive Zr–Hf anomalies (Fig. 6a) are attributed to zircon megacryst contamination and partial assimilation (Lekkerfontein, Gansfontein and Monastery kimberlites).

Calculated initial 87Sr/86Sr and 143Nd/144Nd ratios of the Group I kimberlites from this study, as well as those analysed to complement the datasets presented by le Roex et al. (2003)Go and Harris et al. (2004)Go, are reported in Table 3 and illustrated in Fig. 7. Initial 87Sr/86Sr ratios of the analysed Group I kimberlites vary from 0·70328 ± 2 to 0·70537 ± 2, whereas initial 143Nd/144Nd ratios vary between 0·51249 ± 1 and 0·51271 ± 1, and generally show good overlap with the field of southern African Group I kimberlites (Fig. 7). The strong displacement of the single sample COL 9 to high 87Sr/86Sr (Fig. 7) is attributed to alteration by crustal fluids. Although subtle differences exist in initial 87Sr/86Sr ratios between samples that may be related to the olivine macrocryst content (the more olivine macrocrystic samples tend to have higher 87Sr/86Sr(i)), no significant differences are evident between initial 87Sr/86Sr and 143Nd/144Nd ratios of on- and off-craton Group I kimberlites.


Figure 7
View larger version (23K):
[in this window]
[in a new window]
 
Fig. 7. Variation of initial 143Nd/144Nd and 87Sr/86Sr for on- and off-craton Group I and Group II analysed in this study (Table 3). Isotope data plotted for the Finsch and New Elands kimberlites analysed include data from Fraser & Hawkesworth (1992)Go and Smith (1983)Go. Error bars for 143Nd/144Nd and 87Sr/86Sr represent the magnitude of the change in isotope ratios over 100 Myr. Initial isotope ratios are calculated according to the age or inferred age of intrusion and 2{sigma} errors are less than the symbol size. Labelled shaded fields represent published data for South African Group I, Group II and Transitional kimberlites (Smith, 1983Go; Fraser & Hawkesworth, 1992Go; Tainton, 1992Go; Clark, 1994Go; Nowell et al., 1999Go, 2004Go; Coe, 2004Go; Becker & le Roex, in preparationGo) and South Atlantic OIB (O'Nions & Pankhurst, 1974Go; O'Nions et al., 1977Go; le Roex, 1985Go; le Roex et al., 1990Go).

 
On- and off-craton Group II kimberlites
Although there is considerable scatter in the data, for a given MgO content the analysed Group II kimberlites are characterized by higher SiO2 contents (29·0–43·1 wt %) than the Group I kimberlites (Fig. 2a); this feature is also shown by the on-craton Swartruggens and Star Group II kimberlites (Coe, 2004Go). Unlike the Group I kimberlites, there is no obvious correlation between MgO and SiO2, although some within-pipe correlations do exist (e.g. Markt and Eendekuil kimberlites). Variation in MgO concentrations and Mg-number is large (MgO = 17·1–36·3 wt %; Mg-number = 0·80–0·91), with the more MgO-rich samples being more macrocrystic in character (e.g. Bellsbank, Newlands and Markt kimberlites). The off-craton Eendekuil kimberlite (argued later to be highly fractionated) falls outside these ranges, having substantially lower MgO and Mg-number (Fig. 2a; Table 2). FeO* contents mostly range from 6·4 to 8·0 wt % (Fig. 2b), with the evolved off-craton Eendekuil kimberlite being richest in FeO* (up to 11·5 wt %).

Most samples show low abundances of TiO2 (0·56–1·4 wt %) and Al2O3 (0·79–4·5 wt %), but some (e.g. Eendekuil and Brandewynskuil kimberlites), argued later to show evidence for crustal contamination and crystal fractionation, tend to high TiO2 (>2 wt %) and Al2O3 (>5 wt %). On the whole, the Group II kimberlites show distinctively lower TiO2, but higher K2O (up to 5 wt %) than the Group I kimberlites (Fig. 2c), a characteristic feature noted by Smith et al. (1985b)Go. The on-craton Newlands and Bellsbank Group II kimberlites have unusually low K2O (0·53–1·9 wt %), and for reasons discussed later are referred to as Group IIb kimberlites. Good correlations exist between CaO (4·8–21·3 wt %), CO2 (up to 15·1 wt %) and SiO2 (Fig. 2d); for a given SiO2 content Group II kimberlites tend to have higher CaO than Group I kimberlites. CaO and CO2 concentrations in the analysed Group II kimberlites are more variable than in the Swartruggens and Star kimberlites (Fig. 2d). Other than off-craton Group II kimberlites having generally higher CaO for a given SiO2 content, no other obvious differences are evident in the major element geochemistry of on- and off-craton Group II kimberlites.

The transition metals Ni (570–2062 ppm) and Cr (485–2326 ppm) show a wide range in abundance and correlate well with MgO (Fig. 3a) and Mg-number. Although the off-craton Group II kimberlites tend to have lower Ni than their on-craton counterparts, this might simply reflect the small size of the off-craton dataset and the unusually evolved Eendekuil kimberlite. Of perhaps more significance is that Group II kimberlites tend to have higher Ni concentrations, for a given MgO content than Group I kimberlites (Fig. 3a). LILE and Pb concentrations of Group II kimberlites show poor correlation with the more immobile incompatible elements (e.g. La; Fig. 3b), but are distinctively higher (31–210 ppm Rb; 780–7345 ppm Ba; 17–64 ppm Pb (Sanddrift sample K5/P2 excepted) than in Group I kimberlites (Table 2; Fig. 3).

Incompatible element ratios are particularly useful in distinguishing between Group I and Group II kimberlites, with the latter having significantly lower Ce/Pb (16 ± 12) and higher La/Nb (1·6 ± 0·5), Ba/Nb (24 ± 10) and Th/Nb (0·21 ± 0·06) (e.g. Table 2; Fig. 4). The high values for ratios involving Nb reflect the characteristic depletion in Nb (38–192 ppm) compared with Group I kimberlites (~100–400 ppm; Fig. 3). In contrast, inter-element ratios of other HFSE pairs (e.g. Zr/Hf = 47 ± 5; Nb/Ta = 25 ± 10) are not dissimilar to those of Group I kimberlites (Zr/Hf = 52 ± 6·6; Nb/Ta = 20 ± 5·6).

The chondrite-normalized REE patterns of the analysed Group II kimberlites are sub-parallel and characterized by extreme LREE enrichment (La up to 1600 x chondrite) relative to the HREE (Lu = 2–10 x chondrite; Fig. 5b), with most of the analysed samples falling in the field defined by the on-craton Swartruggens and Star Group II kimberlites; exceptions are Finsch, Middlewater and Sanddrift. La/Yb and La/Sm are generally high and illustrative of the steep LREE slope in chondrite-normalized diagrams (La/YbN = 44–180, La/SmN = 5·3–9·3), with the two on-craton Group IIb kimberlites (Bellsbank and Newlands) having appreciably higher La/YbN (183–300) and La/SmN (11·1–12·3) ratios (Fig. 4d). Of significance is that the Group II kimberlites have distinctly higher La/Sm and La/Yb ratios, but lower Gd/Yb ratios than Group I kimberlites (Fig. 4c and d).

Steep patterns, similar to the chondrite-normalized REE patterns, occur on primitive mantle-normalized trace element diagrams (Fig. 6b) with the analysed samples showing strong enrichment in incompatible elements (e.g. Ba = 112–779 and Th = 87–677 x primitive mantle). Superimposed on the generally smooth patterns, there are distinct relative depletions in Rb, Ti and Sr, enrichment in Pb and a broad depletion in the group of elements K–Ta–Nb (Fig. 6b), features that are common to the Swartruggens and Star Group II kimberlites (Coe, 2004Go). Samples from the Bellsbank and Newlands kimberlites are distinct in having a strong relative depletion in K (K/K* <0·1) and Rb (see insert to Fig. 6b), more akin to the Group I kimberlites (Fig. 6a). Consequently, as mentioned above, Group II kimberlites are subdivided into two varieties, Group II with subdued relative depletion in K, and Group IIb with strong relative depletion in K (and Rb). The Group IIb kimberlites also tend to have larger negative Ti anomalies (Ti/Ti* = 0·2–0·3 vs 0·2–0·8). Positive Pb anomalies on primitive mantle-normalized diagrams are a feature of all but one (Sanddrift kimberlite, which also has a large negative P anomaly) of the analysed Group II kimberlites (Ce/PbN <1; Fig. 6b), with some samples showing fairly large anomalies. The Swartruggens and Star Group II kimberlites are likewise relatively enriched in Pb (Fig. 6b). Most samples also show a minor relative depletion in Hf (Sm/HfN <2·0), but in the Eendekuil kimberlite there is a much better developed distinct negative Zr and Hf anomaly (Sm/HfN ~2·9). Depletion in Ta and Nb is a ubiquitous feature of the analysed Group II kimberlites, as well as of the Swartruggens and Star Group II kimberlites (Coe, 2004Go) with La/NbN ratios generally being >1·3 (Table 2).

Initial 87Sr/86Sr ratios of the Group II kimberlites analysed in this study vary from 0·70728 ± 2 to 0·70837 ± 2, whereas initial 143Nd/144Nd ratios are characteristically unradiogenic and vary from 0·51186 ±1 to 0·51205 ± 1 (Table 3). No significant isotopic differences are evident between on- and off-craton Group II kimberlites. The co-variation of initial 143Nd/144Nd and 87Sr/86Sr ratios shows a scattered distribution that lies within the field previously defined for South African Group II kimberlites (Fig. 7), lying well within the ‘enriched quadrant’ (Smith, 1983Go; Fraser & Hawkesworth, 1992Go; Tainton, 1992Go; Nowell et al., 1999Go, 2004Go; Coe, 2004Go).


    PETROGENESIS
 TOP
 ABSTRACT
 INTRODUCTION
 GEOLOGICAL SETTING AND SAMPLING
 ANALYTICAL TECHNIQUES
 PETROGRAPHY
 GEOCHEMISTRY
 PETROGENESIS
 DISCUSSION
 CONCLUSIONS
 SUPPLEMENTARY DATA
 REFERENCES
 
The volatile-rich and hybrid nature of kimberlite magmas requires that considerable care be given to addressing contamination and alteration processes that might have compromised the kimberlite magma composition subsequent to leaving the source region (e.g. le Roex et al., 2003Go). Before considering possible genetic models, we first consider the possible impact of such secondary processes on the geochemistry of the analysed samples.

Alteration and crustal contamination
The brecciated and permeable nature of kimberlites makes them prone to alteration (Berg & Allsopp, 1972Go), particularly as primary kimberlite magmas are volatile-rich. Alteration of the kimberlite can be both a primary (movement of crystallizing deuteric fluids) and a secondary (movement of meteoric waters) process. Not only do the kimberlites in this study show various petrographic features indicative of alteration (e.g. serpentinization or carbonatization), but these effects can also be identified in the major element geochemistry (e.g. loss of MgO and gain of volatiles during serpentinization), trace element geochemistry (e.g. lack of correlation of the mobile elements Rb and K with more immobile elements such as La; Fig. 3b) and isotope geochemistry (e.g. increase of 87Sr/86Sr ratios, Barrett & Berg, 1975Go; seen most strongly in sample COL 9). The HFSE are generally unaffected by such processes.

Similarly, because kimberlites commonly entrain numerous crustal xenoliths, the degree of assimilation of country rock and its effect on bulk-rock geochemistry needs to be considered. The effects of crustal contamination (e.g. increase in SiO2, Al2O3 and Na2O, and decrease in MgO) can be partially quantified with the contamination index (C.I.) of Clement (1982)Go, or qualitatively recognized using parameters identified by le Roex et al. (2003Go; e.g. raised HREE patterns, raised Pb in Group I kimberlites). Most Group I and Group II kimberlites in this study have a C.I. near unity and less than 1·5, respectively, i.e. classify as being uncontaminated.

Through careful initial sampling and removal of all visibly altered or veined material and xenolith fragments, the effects of alteration and contamination have been minimized. Furthermore, by applying the geochemical criteria listed by le Roex et al. (2003)Go, in combination with Clement's (1982)Go C.I., those analysed compositions that appear to have been severely compromised either are excluded from the following petrogenetic discussions or are specifically mentioned. The interpretation of the geochemistry of the remaining samples is done with caution, with reference being made to occasions when alteration or contamination still needs to be borne in mind.

Macrocryst entrainment
Macrocrysts entrained by kimberlites are typically anhedral with strained physical characteristics that suggest they are xenocrysts, derived from disaggregated mantle peridotite xenoliths (Clement et al., 1984Go; Shee, 1985Go). Consequently, macrocrystic kimberlites are generally considered to be aggregates of mantle minerals as opposed to true liquid compositions (le Roex et al., 2003Go). A number of previous studies (e.g. Fraser & Hawkesworth, 1992Go; Beard et al., 2000Go; le Roex et al., 2003Go; Coe, 2004Go; Harris et al., 2004Go) have demonstrated that some of the variation in kimberlite bulk-rock geochemistry can be attributed to the entrainment of xenoliths of mantle peridotite; this is also true for the present study.

To consider only liquid compositions, the effect of macrocryst entrainment has to be addressed for each kimberlite sample. The regional nature of this study has resulted in few of the kimberlites being sampled in sufficient detail to evaluate evolutionary trends and the effect of peridotite entrainment, within a single kimberlite locality, as has been done in other kimberlite studies (e.g. Harris et al., 2004Go). Consequently, what we have done is to assume that all the olivine macrocrysts in a particular sample are derived from disaggregated garnet lherzolite xenoliths [using an average composition calculated from Grégoire et al. (2003)Go]. Using proportions based on the estimated modal olivine abundance in the kimberlite sample, an equivalent proportion of garnet lherzolite was accordingly subtracted. In view of the inherent errors in estimating modal abundances, this correction procedure was applied only to samples with >10% modal olivine macrocrysts. Compatible trace elements were corrected using an average composition of garnet lherzolite (Grégoire et al., 2003Go), whereas incompatible trace elements were effectively enriched in abundance by the proportion of lherzolite removed (assuming D = 0 for all elements). Although only an approximation, the very low partition coefficients and low incompatible element abundances in garnet lherzolite do not warrant a more complex approach. Given the broad range in 87Sr/86Sr and 143Nd/144Nd of mantle xenoliths entrained by South African kimberlites (e.g. Erlank et al., 1987Go; Walker et al., 1989Go; Pearson & Nowell, 2002Go), as well as the particularly low Nd concentrations in mantle peridotites (<5 ppm Nd), the isotopic compositions of the kimberlites have been left unchanged.

Fractional crystallization
Crystal fractionation en route to the surface or during shallow emplacement and cooling of the kimberlite magma will affect the composition of the primary kimberlite magma. Fractionation of olivine and phlogopite, and possibly calcite, have been proposed, based on detailed studies of individual kimberlites (e.g. Beard et al., 2000Go; le Roex et al., 2003Go; Coe, 2004Go; Harris et al., 2004Go). Evolved, aphanitic kimberlite, in particular, is a good example of the end product of such fractionation processes (le Roex et al., 2003Go). Crystal fractionation/accumulation may have occurred to variable degrees, for example via flow differentiation, in macrocryst-bearing kimberlites (e.g. the evolved Eendekuil kimberlite).

On the basis of a number of lines of evidence, le Roex et al. (2003)Go and Harris et al. (2004)Go have argued that the ubiquitous negative anomalies in K, Rb, Sr and Ti in primitive mantle-normalized trace element diagrams are not related to crystal fractionation processes, but are features of the primary magmas of the Group I Kimberley and Uintjiesberg kimberlites. Following the arguments presented by le Roex et al. (2003)Go and Harris et al. (2004)Go, we suggest that the ubiquitous negative Rb, K and Ti anomalies on primitive mantle-normalized trace element diagrams for the Group I kimberlites in this study are not the result of extensive phlogopite (Rb, K, Ti), ilmenite (Ti) or perovskite (Ti) fractionation (the last is supported also by the absence of an associated negative Nb and Ta anomaly), but are features of the primary kimberlite magmas. However, it is recognized that minor amounts of fractionation, as well as K and Rb mobility, may have an influence on the magnitudes of the anomalies.

Even though Nb and Ta are depleted in Group II kimberlites, the characteristic paucity of ilmenite and perovskite as groundmass phases (e.g. Skinner, 1989Go) argues against fractionation of a Ti-phase as the cause of the ubiquitous negative Ti anomaly and relative depletion in Nb and Ta in primitive mantle-normalized trace element patterns. Based on similar calculations to those used by le Roex et al. (2003)Go and Harris et al., (2004)Go, it can be argued that relative depletion in Rb, K, Sr, Ti, Nb and Ta is a primary feature of normal Group II and Group IIb kimberlites in this study, not caused by the fractional crystallization of phlogopite (K, Ti; Coe, 2004Go). Furthermore, the strong relative depletion in K and Rb present in the Group IIb kimberlites is argued here not to be due to phlogopite fractionation, but to be a feature inherited from the mantle source region.

Composition of close-to-primary kimberlite magmas
To allow meaningful and regional inter-group comparison of the geochemistry of on- and off-craton Group I and Group II kimberlites, it is important to identify a representative composition for each kimberlite that is close-to-primary in character. Because of the hybrid nature of kimberlites, and as no examples of quenched kimberlite exist (Mitchell, 2004Go), aphanitic kimberlite has been considered a proxy for a primary magma (Edgar & Charbonneau, 1993Go; Price et al., 2000Go). More recently, le Roex et al. (2003)Go and Harris et al. (2004)Go have constrained the compositions of primary kimberlite magmas by using the trajectories between macrocrystic and aphanitic kimberlite on major element variation diagrams. However, owing to the regional nature of this study, primary kimberlite magma compositions cannot be constrained here using similar methods, as sufficient samples are not available from each locality. Here we define close-to-primary kimberlite magma compositions as those being least affected by alteration, crustal contamination and fractional crystallization, and, specifically, having been corrected for macrocryst entrainment where necessary (as described earlier). Not every kimberlite represents a close-to-primary magma (e.g. aphanitic Lekkerfontein kimberlite) and in a few cases (e.g. Pampoenpoort, Jonkerwater, Middlewater and Sandrift kimberlites, where samples were made available as powders from De Beers) the absence of thin sections prohibited any correction for macrocryst content; these samples are omitted from further discussion.

Inferred close-to-primary magma compositions of Group I kimberlites analysed in this study have ~22–28 wt % MgO, ~21–30 wt % SiO2, 2·2–3·1 wt % Al2O3, ~10–17 wt % CaO, ~5·0–14 wt % CO2, ~0·2–1·7 wt % K2O and 660–1190 ppm Ni, with Mg-number varying between 0·82 and 0·87; these are broadly similar to primary Group I kimberlite compositions proposed in other studies (Price et al., 2000Go; le Roex et al., 2003Go; Harris et al., 2004Go). Group II kimberlites in this study have close-to-primary magma compositions with ~28–38 wt % SiO2, ~23–29 wt % MgO, ~1·0–4·3 wt % Al2O3, ~1·6–4·6 wt % K2O, 8·1–13 wt % CaO, 1·3–7·4 wt % CO2 and 590–1410 ppm Ni, with Mg-number between 0·86 and 0·89. The close-to-primary compositions of the analysed Group II kimberlites are similar to those of the Swartruggens and Star Group II kimberlites (Coe, 2004Go), but tend be slightly more Mg-rich and less K2O-rich.

Inferred close-to-primary Group I kimberlite magma compositions from this study are also broadly consistent with experimental ‘kimberlite’ compositions (20–36 wt % SiO2, 1·8–3·2 wt % Al2O3, 26–30 wt % MgO), obtained from low-degree melting of a carbonated lherzolite by Dalton & Presnall (1998a)Go, except that the CaO content of the experimental kimberlites is significantly higher (15–31 wt %). The MgO and SiO2 contents of inferred close-to-primary Group II kimberlites in this study also show some overlap with experimental compositions produced by <1% partial melting of a carbonated lherzolite at 6 GPa (Dalton & Presnall, 1998aGo).

Table 4 gives average Group I and Group II close-to-primary kimberlite magma compositions for major and trace elements and selected incompatible element ratios, calculated using data from this study and the literature (le Roex et al., 2003Go; Coe, 2004Go; Harris et al., 2004Go). There is considerable scatter around the averages (as reflected by the reported standard deviations), but nevertheless certain elements stand out as being significantly different between the two kimberlite varieties. These include higher TiO2, Fe2O3, MgO and CaO, and lower SiO2 and K2O in Group I kimberlites, and substantially higher Ba, Rb and Pb, and lower Nb and Ta in Group II kimberlites compared with Group I kimberlites. Although less marked, NiO appears to be slightly higher in Group II kimberlites on average than in Group I kimberlites (Table 4). It is argued later that these differences are inherited from their respective source regions.


View this table:
[in this window]
[in a new window]
 
Table 4: Average compositions and associated standard deviations of close-to-primary Group I and Group II kimberlites from South Africa

 
Partial melting models
Experimental studies have shown that Group I kimberlite-like magmas can be produced by low degrees of partial melting of a carbonated garnet peridotite mantle source at ~3–8 GPa (Eggler & Wendlandt, 1979Go; Canil & Scarfe, 1990Go; Dalton & Presnall, 1998aGo) and that a continuum exists between carbonatitic and kimberlitic melts, with the transition between these magma types taking place between 0 and 1% melting (Dalton & Presnall, 1998a). Ulmer & Sweeney (2002)Go have also shown that Group II kimberlites could be formed by partial melting of a carbonated garnet harzburgite or lherzolite at pressures up to 10 GPa.

The characteristic incompatible element enrichment and fractionated REE patterns of kimberlites are consistent with derivation by very low degrees of partial melting of metasomatically enriched sources (e.g. Wass & Rogers, 1980Go). Partial melting is generally inferred to take place in the presence of residual garnet because of the characteristic low HREE abundances and fractionated HREE patterns in both on- and off-craton kimberlites relative to typical mid-ocean ridge basalt (MORB; Rogers et al., 1992Go; Tainton & McKenzie, 1994Go; le Roex et al., 2003Go). Metasomatized lherzolite xenoliths from the Kaapvaal craton typically contain clinopyroxenes that have been enriched in the LREE (van Achterbergh et al., 2001Go; Grégoire et al., 2003Go) and therefore clinopyroxene is likely to be the major mineral host of these elements during partial melting. The semi-quantitative modelling of kimberlite source compositions outlined below is, therefore, based on a source comprising metasomatized garnet lherzolite and very low degrees of partial melting (F <2%). Some workers have argued that in view of the absence of olivine and orthopyroxene as liquidus phases during the experimental melting of kimberlites and lamproites, the source of such magmas is not a lherzolite (Foley, 1992aGo; Edgar & Charbonneau, 1993Go; Girnis et al., 1995Go; Mitchell, 2004Go), but rather some more exotic mineral assemblage (e.g. clinopyroxene–phlogopite veins; Foley, 1992bGo); this possibility is also explored.

In view of the assumptions required to estimate primary kimberlite magma compositions, complex inversion modelling (e.g. Tainton & McKenzie, 1994Go) is not justified. Furthermore, given the number of assumptions that need to be made when modelling mantle melting processes, and to allow possible differences between the sources of Group I and Group II kimberlites, as well as on- and off-craton kimberlites to be recognized, it is necessary to restrict at least some of the variables. First, close-to-primary magma compositions were used together with a fixed degree (F = 1%) of partial melting. Modal analyses of peridotite xenoliths entrained by on- and off-craton kimberlites (Mathias et al., 1970Go; Maaloe & Aoki, 1977Go; Boyd & Mertzman, 1987Go; Franz et al., 1996Go; Grégoire et al., 2003Go; Boyd et al., 2004Go) tend to be broadly similar, particularly with regard to the proportion of modal garnet (3–9 vol. %); consequently, we have assumed a fixed modal mineralogy for both on- and off-craton source regions. Semi-quantitative forward modelling of batch melting is based on the residual source mineralogy to avoid having to assume the reaction coefficients of the melting equation, using the partition coefficients summarized by le Roex et al. (2003)Go.

Considering the unusual composition of kimberlites, the choice of suitable partition coefficients is limited to a silica-undersaturated system such as those for basaltic, carbonatitic or ultrapotassic melts (e.g. Schmidt et al., 1999Go); none of which are ideal for kimberlites. Partition coefficients for basalts were chosen for modelling for both Group I and Group II kimberlites, because of their availability as a coherent dataset of elements. Differences in the CO2 content of close-to-primary Group I (0·2–1·7%) and Group II kimberlites (1·3–7·4%) were not considered large enough to justify the use of carbonatitic partition coefficients for Group II kimberlites.

The observation that Group II kimberlites have on average steeper LREE (higher La/Sm), but flatter HREE slopes (lower Gd/Yb) than Group I kimberlites (Fig. 4), could imply differences in the source mineralogy. However, calculated melting trajectories (Fig. 8) illustrate that increasing the degree of partial melting decreases La/Sm and Gd/Yb, but cannot account for the differences between the kimberlite groups. Rather, it appears that Group II kimberlites are derived from sources with slightly higher La/Sm but lower Gd/Yb than Group I kimberlites. Given the difference in partition coefficients between La, Sm, Gd and Yb, this is consistent with derivation of Group II kimberlites from a source with less residual garnet but more residual clinopyroxene than that giving rise to Group I kimberlites. The residual modal mineralogy used for subsequent modelling is 67% olivine, 26% orthopyroxene, 3% garnet, 4% clinopyroxene for the Group I kimberlites. The residual modal mineralogy used for Group II kimberlites has the same olivine and orthopyroxene proportions (67% ol; 26% opx) as for Group I kimberlites but with less garnet (1%) and more clinopyroxene (6%). Although the measured La/Sm and Gd/Yb ratios of Kaapvaal garnet lherzolite xenoliths are variable, and show little correlation with the proportion of modal garnet or clinopyroxene (van Achterbergh et al., 2001Go; Grégoire et al., 2003Go; Simon et al., 2003Go), the values used in the modelling and shown in Fig. 8 fall within the observed ranges.


Figure 8
View larger version (22K):
[in this window]
[in a new window]
 
Fig. 8. Gd/Yb vs La/Sm for on- and off-craton Group I and Group II close-to-primary kimberlite magmas, including the Uintjiesberg kimberlite (Harris et al., 2004Go). Dashed field represents the Group IIb kimberlites. Illustrated curves represent melting trajectories of inferred Group I and Group II kimberlite source regions (La/Sm = 6·0; Gd/Yb = 3·0 and La/Sm = 6·4; Gd/Yb = 2·5, respectively). Residual modal mineralogy as follows: Group I, ol:opx:cpx:gt = 0·67:0·26:0·04:0·03; Group II, ol:opx:cpx:gt = 0·67:0·26:0·06:0·01. Numbers shown represent the degree of melting. Shaded fields represent the on-craton Group I Kimberley (le Roex et al., 2003Go) and on-craton Group II Swartruggens and Star kimberlites (Coe, 2004Go).

 
Source region characteristics
Calculated compositions of mantle source regions in equilibrium with close-to-primary Group I kimberlites are all enriched in the LREE (4–25 x chondrite) and only slightly enriched in the HREE relative to chondrite (1·8–3·6 x chondrite; Fig. 9a). They have La/SmN = 2·3–3·9 and La/YbN = 1·9–11·5, but are depleted in the less incompatible elements relative to primitive mantle (Yb = 0·5–1 x primitive mantle; Fig. 9b). Similarly, calculated sources are also enriched in the more incompatible elements relative to primitive mantle (e.g. Th = 1–9 x primitive mantle). In the absence of residual phases capable of fractionating individual elements, the negative anomalies present in the close-to-primary magma compositions are transferred to the predicted source regions (Fig. 9b). Consequently, the calculated sources have negative Rb, Sr, K and Ti anomalies (K/K* <0·2; Ti/Ti = 0·2–0·7), as well as more subdued negative Hf anomalies. Predicted source compositions generally show good overlap with those calculated for the Group I on-craton Kimberley kimberlites and off-craton Uintjiesberg kimberlite (Harris et al., 2004Go). Calculated source compositions overlap the field of garnet lherzolites from the Kaapvaal craton for all elements other than K, which is significantly more depleted in the predicted sources than in peridotites (Fig. 9b). However, many of the peridotites from the study of Grégoire et al. (2003)Go contain modal phlogopite and are thus likely to have elevated K concentrations.


Figure 9
View larger version (37K):
[in this window]
[in a new window]
 
Fig. 9. Calculated chondrite- and primitive mantle-normalized source compositions of on- and off-craton Group I and Group II kimberlites, calculated by semi-quantitative forward modelling from close-to-primary kimberlite magmas assuming F = 1% and residual source region mineralogy as in Fig. 8. Calculated source compositions of the on-craton Group I Kimberley kimberlites (le Roex et al., 2003Go), off-craton Group I Uintjiesberg kimberlite (Harris et al., 2004Go) and on-craton Swartruggens and Star Group II kimberlites (Coe, 2004Go) are shown for comparison. Dark shaded field represents garnet lherzolite xenoliths from the Kaapvaal craton (Grégoire et al., 2003Go). Normalizing values are from Sun & McDonough (1989).

 
Because both on- and off-craton source compositions are calculated using the same parameters, any differences in the respective kimberlite magma compositions will be transferred to the sources. Consequently, predicted source compositions of off-craton Group I kimberlites appear to have lower La/Sm, La/Yb and Gd/Yb than on-craton Group I kimberlites. Alternatively, if the on- and off-craton sources were compositionally similar, then the off-craton kimberlites formed by slightly higher degrees of partial melting, as all Group I kimberlites fall roughly on the same melting trajectory (Fig. 8).

Predicted source regions in equilibrium with close-to-primary Group II kimberlite magmas are all enriched relative to chondrite, with the LREE (5–37 x chondrite) being considerably more enriched than the HREE (0·9–2·6 x chondrite; Fig. 9c). The sources of the Group IIb kimberlites are more enriched in the incompatible elements with higher La/Sm and La/Yb (La/SmN ~ 6·9) than those of normal Group II kimberlites (La/SmN = 3·2–5·4), suggesting that they are possibly lower degree partial melts than some of the other Group II kimberlites analysed in this study (Fig. 8). Calculated sources are also enriched in the highly incompatible LILE and HFSE and depleted in the less incompatible elements, relative to primitive mantle (e.g. Th = 1–11 x primitive mantle and Yb = 0·3–0·8 x primitive mantle, respectively; Fig. 9d). Negative Rb, K and Ti anomalies are larger for the sources of Group IIb kimberlites (K/K* ~0·1, Ti/Ti* ~0·2) than normal Group II kimberlites (K/K* = 0·3–1·0, Ti/Ti* = 0·2–0·6). All Group II kimberlite source regions show subdued depletion in Sr, Ta and Nb (La/NbN = 1·5–3·9), and slight enrichment in Pb (Ce/PbN = 0·2–1·0). Calculated source compositions (Fig. 9d) fall within the field of garnet lherzolite xenoliths from Grégoire et al. (2003)Go. There appear to be no significant differences in calculated source compositions between on- and off-craton Group II kimberlites.

A number of recent studies of kimberlite petrogenesis based on bulk-rock geochemistry (Tainton & McKenzie, 1994Go; Beard et al., 1998Go; Beard et al., 2000Go; le Roex et al., 2003Go; Chalapathi Rao et al., 2004Go; Coe, 2004Go; Harris et al., 2004Go) have all argued that the source regions of kimberlites have (at least) a two-stage evolutionary history, comprising initial melt-depletion with subsequent metasomatic enrichment in incompatible elements. Close-to-primary magmas of on- and off-craton Group I and Group II kimberlites all have high Mg-numbers (Mg-number = 0·82–0·89) and compatible trace element concentrations (e.g. 650–1400 ppm Ni) that suggest equilibration with residual olivine of Fo92–96, having more than 2000 ppm Ni [calculated using KDFe–Mg = 0·34–0·36 (Herzberg & O'Hara, 2002Go) and DNi from Beattie et al. (1991)Go]. It is notable that the calculated compositions of residual olivine in equilibrium with close-to-primary kimberlites tend to be more similar to refractory olivine included in diamonds than to olivine from Kaapvaal peridotite xenoliths (e.g. Phillips et al., 2004Go). All kimberlites have low Al2O3 and Sc concentrations, as well as relative depletion in the HREE and fractionated HREE patterns indicative of the role of residual garnet, and suggest derivation from previously depleted mantle source regions (Tainton & McKenzie, 1994Go). Studies of peridotite xenoliths from the Kaapvaal craton and surrounding Proterozoic mantle have recognized that the SCLM beneath southern Africa has a depleted or infertile character, deficient in Fe, Ca and Al in comparison with oceanic mantle; additionally the Archaean mantle is more refractory than the younger Proterozoic mantle (Boyd & McCallister, 1976Go; Boyd & Mertzman, 1987Go; Boyd et al., 2004Go). Given the apparent initial depleted character of kimberlite source regions, we concur with Tainton & McKenzie (1994)Go, le Roex et al. (2003)Go and Harris et al. (2004)Go, who have argued that they are more likely to be located within the SCLM than within the convecting asthenosphere.

Residual accessory phases
As forward modelling is based on a residual garnet lherzolite source with no accessory phases present, bulk partition coefficients for Rb, K, Ti, Sr and Nb during partial melting are all substantially less than unity; thus any anomalies in the close-to-primary kimberlite magmas are transferred to their source regions. Although some workers have argued for the presence of a residual K-bearing phase during partial melting to be the cause of negative K anomalies in some other alkaline magmas (e.g. Rogers et al., 1992Go; Späth et al., 2001Go; Williams et al., 2004Go), others such as le Roex et al. (2003)Go and Harris et al. (2004)Go have argued that these are primary features of kimberlite source regions.

Notwithstanding the possible superimposed effect of low-temperature element mobility, evaluation of the ubiquitous negative Rb and K anomalies in the Group I kimberlites analysed in this study leads us to the same conclusion regarding the origin of the apparent depletion in these elements, i.e. that it is a primary source feature, and does not result from buffering by residual phlogopite. The role of residual phlogopite in Group II kimberlites is more ambiguous. Unlike Group I kimberlites, the absolute K content of Group II kimberlites is consistent with the presence of residual phlogopite (based, for example, on stoichiometry considerations and melting modes; Späth et al., 2001Go), although most show only a minor relative depletion in K, suggesting that if present initially, phlogopite was fully or close to fully consumed during melting. This is consistent with the experimental evidence that phlogopite is not stable on the Group II kimberlite liquidus during melting in the presence of carbonate at pressures greater than 4–6·5 GPa, but rather breaks down to form garnet, orthopyroxene and a K-rich liquid (Yamashita et al., 1995Go; Ulmer & Sweeney, 2002Go). An exception to this generalization might be the Group IIb kimberlites that show strong relative depletion in K (Fig. 6b). The radiogenic 87Sr/86Sr character of Group II kimberlites is consistent with the presence of phlogopite in the source.

Group II kimberlites show a distinct, and characteristic, depletion in Nb and Ta relative to elements of similar incompatibility (e.g. Fig. 6b), which suggests possible buffering against a residual titanate phase, such as rutile. Because rutile is typically a metasomatic phase (Haggerty, 1987Go), it is likely to be preferentially melted. However, simple calculations involving the Ti budget and melt mode suggest that the presence of residual rutile during partial melting is unlikely, as it would be entirely consumed after <0·5% partial melting (Coe, 2004Go). It is, therefore, interpreted that the relative depletion in Nb and Ta (La/NbN >1·3; Fig. 4b) is a characteristic source feature of all Group II kimberlites analysed in this study. Furthermore, the presence of negative Nb and Ta anomalies in Karoo basalts (Marsh et al., 1997Go), argued by some workers to represent subcontinental lithospheric melts (e.g. Hawkesworth et al., 1984Go), as well as in some Kaapvaal craton xenoliths and their constituent minerals (van Achterbergh et al., 2001Go; Grégoire et al., 2003Go), suggests that this may be an important characteristic of the sub-Gondwana lithosphere.

Close-to-primary kimberlite magmas are, therefore, argued to have been unaffected by residual phlogopite or any titanate phases. Following the preferred models of le Roex et al. (2003)Go and Harris et al. (2004)Go for Group I and Coe (2004)Go for Group II close-to-primary kimberlite magmas, depletion in Rb, K, Ti and, to a lesser extent, Sr, Hf and Nb–Ta (Group II kimberlites) are characteristics inherited directly from their source regions.

Timing of source enrichment
Nd isotope signatures of on- and off-craton Group I kimberlites are depleted relative to present-day Bulk Earth; minimum Nd model ages calculated with respect to CHUR and Depleted Mantle suggest that metasomatic enrichment of their source regions occurred relatively recently (0–700 Ma), possibly just prior to the Jurassic–Cretaceous emplacement of Group I kimberlites. The Group II kimberlite sources, however, are characterized by time-integrated 143Nd/144Nd ratios indicative of ancient LREE source enrichment (this study; Smith, 1983Go; Fraser et al., 1985–1986Go; Coe, 2004Go). Calculated model ages suggest source enrichment at 600–900 Ma relative to CHUR, or 1000–1300 Ma relative to Depleted Mantle, but, most importantly, they indicate that the source regions of Group II kimberlites were metasomatically enriched considerably prior to those of Group I kimberlites. Similar geochemical compositions and Nd model ages for both on- and off-craton Group II kimberlite source regions suggest enrichment after the collision of the Proterozoic Namaqua–Natal belt and the Kaapvaal craton at ~1·2–1·0 Ga (Thomas et al., 1994Go).


    DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 GEOLOGICAL SETTING AND SAMPLING
 ANALYTICAL TECHNIQUES
 PETROGRAPHY
 GEOCHEMISTRY
 PETROGENESIS
 DISCUSSION
 CONCLUSIONS
 SUPPLEMENTARY DATA
 REFERENCES
 
The expanded geochemical dataset for both Group I and Group I kimberlites provided by this study affords the opportunity to investigate the cause of similarities and differences more rigorously than has been possible in the past, to investigate geodynamic relationships, to compare the influence of on- and off-craton settings, and to evaluate the possible origin of the metasomatic fluids that have enriched their mantle source regions.

On-craton vs off-craton differences
The role of tectonic setting, in terms of on- or off-craton locations, in controlling kimberlite geochemistry is ambiguous. Group II kimberlites show little systematic difference in major and trace element or Nd–Sr isotope geochemistry between on- and off-craton varieties, suggesting that both Proterozoic and Archaean mantle sources have experienced similar evolutionary histories leading up to Group II kimberlite magmatism. Both on- and off-craton Group II kimberlites appear to derive from within the garnet stability field, although the presence of diamonds in the on-craton kimberlites suggests derivation from greater depths. In contrast, off-craton Group I kimberlites tend to fall towards the low SiO2 and MgO, but high FeO*, TiO2, CaO and CO2 ends of the compositional fields of their on-craton counterparts (Fig. 2), suggesting possible derivation from more fertile mantle sources (Boyd & McCallister, 1976Go; Janney et al., 2001Go; Boyd et al., 2004Go). Alternatively, experimental studies suggest that partial melting at lower pressures results in a decrease in SiO2 and MgO, and increase in FeO, Al2O3, CaO and CO2 (e.g. Herzberg, 1992Go; Dalton & Presnall, 1998bGo; Gudfinnsson & Presnall, 2003Go), which would be consistent with the absence of diamonds in off-craton Group I kimberlites (Clifford, 1966Go), although both varieties are still derived from within the garnet stability field, based on similarly fractionated HREE patterns (Fig. 5a). LREE/HREE ratios also suggest that off-craton Group I kimberlites are higher degree melts than on-craton varieties (Fig. 8), in disagreement with the experimental studies of Dalton & Presnall (1998a)Go and Herzberg (1992)Go, which indicate that partial melts become progressively richer in SiO2 and MgO as the degree of melting increases.

Alternatives to a garnet lherzolite source
Given that the source regions of Group II kimberlites appear to be more highly enriched in incompatible elements than those of Group I kimberlites, a more exotic source mineralogy, corresponding to the ‘clinopyroxene and mica rich-vein’ model of Foley (1992b)Go needs to be considered. Foley (1992b)Go has argued that the source of ultrapotassic rocks is free of olivine, because of the absence of olivine as a liquidus phase in various experimental studies (e.g. Foley, 1992aGo; Girnis et al., 1995Go), and that the resulting magmas are hybrids of vein (e.g. cpx, phlog, apatite, spinel) and wall-rock components (e.g. ol, opx, cpx, garnet) with initial melting taking place within the vein. The presence of cross-cutting veins in mantle samples is well established (e.g. Erlank et al., 1987Go); these may represent the channelways for metasomatic fluids that then cryptically or modally metasomatize the wall rock, or are the final product of extensive modal metasomatism.

Grégoire et al. (2002)Go have recently published REE abundances in phlogopite and clinopyroxene from phlogopite-rich mantle xenoliths from the Kimberley region, which allows evaluation of the Foley (1992b)Go model. Clinopyroxene associated with phlogopite in such xenoliths appears to have chondrite-normalized REE patterns that are slightly convex upward with respect to the middle REE, i.e. La and Ce are slightly depleted relative to Nd and Sm, whereas phlogopite is virtually devoid of REE. The bulk composition thus essentially reflects uniform dilution of the clinopyroxene concentrations by the amount of phlogopite present. If one assumes a vein source comprising 50% each of clinopyroxene and phlogopite, then very low degrees of melting (1%) will give REE patterns reasonably similar with respect to the HREE, but very distinct with respect to the LREE from those seen in Group II kimberlites (Fig. 10), with the latter having very much steeper LREE patterns (La/SmN = 8 ± 2, Table 4, vs La/SmN ~2 in the vein melts). One per cent melting is probably low for melting within a metasomatic vein, and more realistic higher degrees of melting will only exacerbate the discrepancy.


Figure 10
View larger version (20K):
[in this window]
[in a new window]
 
Fig. 10. Calculated chondrite-normalized REE patterns for a 1% partial melt of a clinopyroxene–phlogopite vein compared with average close-to-primary Group II kimberlite compositions. Vein composition calculated assuming 50% clinopyroxene and 50% phlogopite, using published REE abundances from clinopyroxene and phlogopite in clinopyroxene–phlogopite xenoliths from the Kimberley region (Grégoire et al., 2002Go). Chondrite normalizing values are from Sun & McDonough (1989).

 
High Mg-number as well as Ni and Cr contents of close-to-primary Group II kimberlite magmas provide a further argument against a clinopyroxene–phlogopite vein source, in that they require equilibrium with refractory olivine at some stage during their petrogenesis. Similarly, Re–Os isotopes and Os population studies imply that Group II kimberlites are derived from melting of olivine-rich sources (Pearson et al., 1995Go, 2003Go). These observations together suggest that if vein melting was involved, then considerable wall-rock interaction took place to give the kimberlites their refractory character.

Metasomatic veining is a variant of disseminated modal metasomatism (e.g. Menzies et al., 1987Go), and the modelling technique used here is unable to discriminate the physical nature of metasomatic enrichment to this level of detail. Moreover, it is based on the modal residual mineralogy at the time of melt extraction; we therefore favour a model of melt generation from a modally metasomatized garnet lherzolite.

Role and origin of metasomatizing fluids
The geochemically enriched nature of the calculated Group I and Group II kimberlite source regions, coupled with the evidence for refractory compositions (high Mg-number and Ni; low Al2O3 and HREE), argues for metasomatic enrichment of a previously melt-depleted source region. Many peridotite xenoliths entrained by kimberlites show evidence for progressive metasomatism reflected in increasingly greater modal proportions of phlogopite, ± K-richterite, but with lower proportions of garnet (Erlank et al., 1987Go; Waters & Erlank, 1988Go). Textural and geochemical evidence indicates that phlogopite, clinopyroxene and Cr-spinel all grow at the expense of garnet (Erlank et al., 1987Go; van Achterbergh et al., 2001Go; Grégoire et al., 2003Go; Simon et al., 2003Go). These reactions are consistent with our earlier arguments that the distinct differences in REE ratios between Group I and Group II kimberlites, and the more Si, K and LILE-enriched character of the latter, are due to their sources having more clinopyroxene and phlogopite formed by metasomatism at the expense of garnet than those of Group I kimberlites. Group II kimberlites, are therefore interpreted to be derived from more highly metasomatized source regions than Group I kimberlites, and although there may be more phlogopite present in the source of the former, this phase is likely to be entirely consumed during partial melting (Ulmer & Sweeney, 2002Go).

It has been argued above that the source of both Group I and Group II kimberlites is the melt-depleted SCLM, enriched through percolation of metasomatic melts or fluids derived from sub-lithospheric sources. Hf isotope evidence from at least the MARID xenoliths in Group I kimberlites supports the involvement of the SCLM and asthenospheric mantle (Choukroun et al., 2005Go). Differences in diagnostic trace element ratios and Sr and Nd isotope compositions of Group I and Group II kimberlites allow a comment on the likely origin of the metasomatic fluids responsible for the enrichment of their respective source regions. Group I kimberlites have unradiogenic Sr and radiogenic Nd isotope ratios, similar to ocean island basalts (OIB; Fig. 7), which led Smith (1983)Go to propose that their sources contain an asthenospheric component. More recently, le Roex et al. (2003)Go and Harris et al. (2004)Go have shown that certain diagnostic incompatible trace element ratios in Group I kimberlites correspond to those characteristic of OIB (e.g. Nb/U = 47 ± 10; Ce/Pb = 25 ± 5; Hofmann et al., 1986Go). Figure 11 shows selected incompatible element ratios of Group I kimberlites from this study and the correlation with South Atlantic OIB compositions is clear, with Group I kimberlites having Ba/Nb <12, La/Nb <1·1, Ce/Pb >22 and Th/Nb <0·14. Given this strong correlation, and following le Roex et al. (2003)Go and Harris et al. (2004)Go, we argue that the enrichment of Group I kimberlite SCLM source regions is a consequence of upward percolation of alkaline melts or fluids sourced from upwelling Mesozoic mantle plumes that impinged on the base of the lithosphere just prior to continental break-up. Nd model ages are consistent with this relatively recent enrichment. In a geodynamic context, the plumes are inferred to have provided both material, in the form of metasomatizing alkaline fluids or melts, and heat to initiate melting to give rise to Group I kimberlite magmatism. In this regard, le Roex (1986)Go, le Roex et al. (2003)Go and Harris et al. (2004)Go have noted that the age of Cretaceous Group I kimberlite magmatism corresponds to the passage of southern Africa over the Shona and Bouvet mantle plumes, currently located in the South Atlantic. By analogy, and although not forming part of this study, source regions of older Group I southern African kimberlites (e.g. Premier, Kuruman cluster) would have been metasomatized by melts or fluids derived from sub-lithospheric upwellings beneath the craton during the Proterozoic.


Figure 11
View larger version (25K):
[in this window]
[in a new window]
 
Fig. 11. Variation of (a) Th/Nb vs La/Nb and (b) Ce/Pb vs Ba/Nb for on- and off-craton Group I and Group II kimberlites. Dividing lines for Th/Nb, La/Nb, Ce/Pb and Ba/Nb provide distinction between kimberlite groups. Field of South Atlantic ocean island basalts (compilation from PetDB; www.petdb.ldeo.columbia.edu) and schematic field for calc-alkaline magmas (compilation from PetDB) are shown for comparison. Data from Kimberley (le Roex et al., 2003Go), Uintjiesberg (Harris et al., 2004Go), Swartruggens and Star (Coe, 2004Go) are included.

 
Group II kimberlites have an isotope signature that requires ancient enrichment of their source regions, and subsequent isolation from the convecting mantle to allow development of their more extreme isotopic compositions that are unlike OIB (Fig. 7). In contrast to Group I kimberlites, diagnostic elemental ratios of Group II kimberlites are clearly very distinct from those characteristic of OIB (Fig. 11), largely as a result of their relative depletion in Nb and Ta, and enrichment in Pb. Group II kimberlites also have lower absolute TiO2 contents with larger negative Ti anomalies than Group I kimberlites. Depletion in Nb and other HFSE is a classic signature of calc-alkaline magmas in subduction zones (e.g. Elliott et al., 1997Go; Woodhead et al., 1998Go), and it is evident from Fig. 11 that the fields for Group II kimberlites strongly overlap with those for calc-alkaline magmas. Fluids derived from subduction zones are characterized by raised Ba/Th ratios, whereas melts derived from the down-going sediment pile tend to have raised Th/Nb ratios (e.g. Elliott et al., 1997Go; Woodhead et al., 1998Go). Group II kimberlites trend towards high Ba/Th and Th/Nb ratios, suggesting that both subduction zone fluids and sediment melt were involved in the metasomatic overprint of their lithospheric mantle source regions. This is consistent with the model of Helmstaedt & Gurney (1984)Go, who suggested that the sources of kimberlites were metasomatized by volatiles derived from subducted slabs.

Minimum Nd model ages and the similar geochemical characteristics of on- and off-craton Group II kimberlites suggest that metasomatic enrichment occurred after welding of the Proterozoic mantle underlying the Namaqua–Natal belt to the Archaean Kaapvaal craton. As no differences have been recognized between the isotopic and trace element signatures of on- and off-craton kimberlites, both Archaean and Proterozoic kimberlite source regions appear to have been enriched by the same metasomatic event. If multiple subduction and accretion events occurred during the 1·2–1·0 Ga formation and metamorphism of the Namaqua–Natal belt (Thomas et al., 1994Go), similar to that imaged by the Lithoprobe study in Canada (e.g. Cook et al., 1999Go), then upward percolating melts or fluids derived from the down-going slabs would have affected both overlying Proterozoic and Archaean mantle. Metasomatism by subduction zone fluids associated with the Pan-African event in southern Africa may also have contributed to the enrichment of on- and off-craton Group II kimberlite mantle source regions.

Independent support for this source region characteristic is found in the compositions of clinopyroxene in mantle xenoliths from both the Kaapvaal craton (e.g. Grégoire et al., 2003Go), and the Namaqua–Natal belt (A. P. le Roex & C. Class, unpublished data, 2005). In both tectonic settings, clinopyroxenes with two distinctly different geochemical signatures can be identified, both showing evidence for metasomatic enrichment, but one having strong depletion in Nb and Ta and greater LREE enrichment.

McKenzie & O'Nions (1983)Go have argued that the sources of mantle plumes may be delaminated SCLM and thus, in theory, such plumes could be the source of metasomatic fluids with the Group II kimberlite geochemical signature. However, this is deemed unlikely as no OIB within the South Atlantic or Southern Ocean shows depletion in Nb or isotopic signatures similar to Group II kimberlites (le Roex & Erlank, 1982Go; le Roex, 1985Go; Hart, 1988Go; le Roex et al., 1990Go). Although Tristan, Gough and Discovery OIB have high Ba/Nb ratios similar to Group II kimberlites, this is a feature of high Ba in EM-1 OIB, thought to reflect a subducted sediment component in the mantle plume source, rather than depletion in Nb (le Roex, 1985Go; Sun & McDonough, 1989Go; Weaver, 1991Go). The palaeo-tracks of these particular OIB plumes also suggest that they are unlikely to have influenced South African kimberlite SCLM source regions (Hartnady & le Roex, 1985Go; le Roex, 1986Go; O'Connor & le Roex, 1992Go).

It is noteworthy that the 180 Ma Karoo continental flood basalts show similar depletion in Nb, Ta and Ti, with strongly enriched isotope signatures (Duncan et al., 1984Go; Hawkesworth et al., 1984Go), suggesting that they may have been derived from a common source region located within the sub-Gondwana continental lithosphere. After the massive outpouring of the Karoo continental flood basalts at 180 Ma (Marsh et al., 1997Go), and the eruption of Group II kimberlites from 200 to 110 Ma (e.g. Skinner, 1989Go; Gurney et al., 2005Go), the more highly metasomatized source appears to have been largely exhausted. Consequently, the younger (~110–70 Ma; Smith et al., 1985aGo, 1994Go) on- and off-craton Group I kimberlites analysed in this study were derived from the more recently (Mesozoic) enriched lithospheric mantle carrying a plume signature. The off-craton Group I East Griqualand kimberlites (Abbotsford East and Zeekoegat) erupted at ~150 Ma (Davis, 1977Go; Smith et al., 1985a)Go, prior to many of the Group II kimberlites, were derived from a mantle source region apparently unaffected by subduction zone metasomatism, but enriched by melts or fluids that have been associated with the Bouvet mantle plume, inferred to have underlain East Griqualand at that time (le Roex, 1986Go). Although the geochemistry of Group II kimberlites shows no evidence for a plume component, it is possible that the Mesozoic mantle plumes lying beneath South Africa at the time of Gondwana break-up provided the heat source for partial melting within the overlying metasomatized lithospheric mantle that led to Group II kimberlite magmatism.


    CONCLUSIONS
 TOP
 ABSTRACT
 INTRODUCTION
 GEOLOGICAL SETTING AND SAMPLING
 ANALYTICAL TECHNIQUES
 PETROGRAPHY
 GEOCHEMISTRY
 PETROGENESIS
 DISCUSSION
 CONCLUSIONS
 SUPPLEMENTARY DATA
 REFERENCES
 
Detailed bulk-rock major and trace element and Nd–Sr isotope data of over 20 Group I and Group II kimberlites emplaced on and off the Kaapvaal craton allow recognition of clear similarities and differences in their geochemistry, and allow constraints to be placed on their genesis and the evolution of their mantle source regions.

(1) Group I kimberlites tend to have higher TiO2, CaO and CO2, but lower SiO2 and K2O, compared with Group II kimberlites, whereas Group II kimberlites have slightly higher Ni, and substantially higher Ba, Rb and Pb and lower Nb and Ta concentrations. Diagnostic trace element ratios of Group I kimberlites that distinguish them from Group II kimberlites include Ba/Nb <12, La/Nb <1·1, Ce/Pb >22 and Th/Nb <0·14. These differences are interpreted to reflect differences inherited from their respective source regions.

(2) Marginally lower SiO2 and MgO, but higher TiO2, CaO and CO2, in off-craton Group I kimberlites compared with their on-craton counterparts suggest derivation of the former from slightly less depleted source regions, or at slightly lower pressures, than on-craton kimberlites. The lower La/Sm, La/Yb and Gd/Yb of off-craton Group I kimberlites also are consistent with derivation by slightly higher degrees of melting than on-craton kimberlites. No significant differences are observed between on- and off-craton Group II kimberlites.

(3) Semi-quantitative modelling indicates that the high Mg-number (0·82–0·89) and Ni (600–1400 ppm) contents coupled with strong enrichment in incompatible elements of both Group I and Group II close-to-primary kimberlites requires derivation from a refractory source that has subsequently been metasomatically enriched. Calculated source compositions are not unlike those observed in garnet lherzolites from the Kaapvaal craton, with La = 4–37 x chondrite and Yb = 0·9–3·6 x chondrite.

(4) Contrasting REE ratios between Group I and Group II kimberlites (higher La/Sm and La/Yb, and lower Gd/Yb in Group II kimberlites) suggest that the source regions of Group II kimberlites have a greater modal proportion of clinopyroxene but less garnet, consistent with being more highly metasomatized than the Group I kimberlite source regions.

(5) The requirement for a refractory source protolith that was subsequently metasomatically enriched prior to melting is taken to imply that the source region of both Group I and Group II kimberlites is located within the SCLM. Characteristic trace element and Nd–Sr isotope ratios indicate that source regions of on- and off-craton Group I kimberlites have been enriched by OIB-like fluids or melts associated with the passage of Mesozoic mantle plume(s) beneath the SCLM. Diagnostic Nd–Sr isotope and trace element ratios of Group II kimberlites are unlike those of Group I kimberlites or OIB, but show more affinity to calc-alkaline magmas. We interpret this to suggest ancient metasomatic enrichment associated with the infiltration of subduction zone derived fluids or melts, which affected both the Archaean and Proterozoic lithospheric source regions of Group II kimberlites. Mesozoic mantle plumes may have subsequently provided the heat source to initiate partial melting and Group II kimberlite generation at the time of continental break-up.


    SUPPLEMENTARY DATA
 TOP
 ABSTRACT
 INTRODUCTION
 GEOLOGICAL SETTING AND SAMPLING
 ANALYTICAL TECHNIQUES
 PETROGRAPHY
 GEOCHEMISTRY
 PETROGENESIS
 DISCUSSION
 CONCLUSIONS
 SUPPLEMENTARY DATA
 REFERENCES
 
Supplementary data for this paper are available at Journal of Petrology online.


    ACKNOWLEDGEMENTS
 
We are grateful to Dr Jock Robey for his help in the field, and Dr Craig Smith for making certain samples in the De Beers sample store available to us. Finsch kimberlite analyses are from Noxolo Zwane's unpublished B.Sc.(Hons) project. Constructive reviews were received from Rick Carlson, Bill Griffin and Goonie Marsh. The National Research Foundation, De Beers Consolidated Mines (Ltd) and the University of Cape Town have all provided financial support for this research.


* Corresponding author. E-mail: aleroex{at}geology.uct.ac.za


    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 GEOLOGICAL SETTING AND SAMPLING
 ANALYTICAL TECHNIQUES
 PETROGRAPHY
 GEOCHEMISTRY
 PETROGENESIS
 DISCUSSION
 CONCLUSIONS
 SUPPLEMENTARY DATA
 REFERENCES
 
Allsopp, H. L. & Barrett, D. R. (1975). Rb–Sr age determinations on South African kimberlite pipes. Physics and Chemistry of the Earth 9, 605–617.[CrossRef]

Barrett, D. R. & Berg, G. W. (1975). Complementary petrographic and strontium-isotope ratio studies of South African kimberlite. Physics and Chemistry of the Earth 9, 619–635.[CrossRef]

Beard, A. D., Downes, H., Hegner, E., Sablukov, S. M., Vetrin, V. R. & Balogh, K. (1998). Mineralogy and geochemistry of Devonian ultramafic minor intrusions of the southern Kola Peninsula, Russia: implications for the petrogenesis of kimberlites and melilitites. Contributions to Mineralogy and Petrology 130, 288–303.[CrossRef][Web of Science]

Beard, A. D., Downes, H., Hegner, E. & Sablukov, S. M. (2000). Geochemistry and mineralogy of kimberlites from the Arkhangelsk Region, NW Russia: evidence for transitional kimberlite magma types. Lithos 51, 47–73.[CrossRef][Web of Science]

Beattie, P., Ford, C. & Russell, D. (1991). Partition coefficients for olivine–melt and orthopyroxene–melt systems. Contributions to Mineralogy and Petrology 109, 212–224.[CrossRef][Web of Science]

Becker, M. & le Roex, A. P. (in prep). Geochemistry and petrogenesis of South African isotopically transitional kimberlites located on and off the Kaapvaal Craton. Lithos.

Berg, G. W. & Allsopp, H. L. (1972). Low 87Sr/86Sr ratios in fresh South African kimberlites. Earth and Planetary Science Letters 16, 27–30.[CrossRef][Web of Science]

Birch, G. F. (1981). The karbonat-bombe; a precise, rapid and cheap instrument to determine calcium carbonate in sediments and rocks. Transactions of the Geological Society of South Africa 84, 199–203.

Boyd, F. R. & McCallister, R. H. (1976). Densities of fertile and sterile garnet peridotites. Geophysical Research Letters 3, 509–512.[Web of Science]

Boyd, F. R. & Mertzman, S. A. (1987). Composition and structure of the Kaapvaal lithosphere, southern Africa. Geochemical Society, Special Publications 1, 13–24.

Boyd, F. R. & Nixon, P. H. (1975). Origins of the ultramafic nodules from the kimberlites of northern Lesotho and the Monastery Mine, South Africa. Physics and Chemistry of the Earth's Interior 9, 431–453.

Boyd, F. R., Pearson, D. G., Hoal, K. O., Hoal, B. G., Nixon, P. H., Kingston, M. J. & Mertzman, S. A. (2004). Garnet lherzolites from Louwrensia, Namibia: bulk composition and P/T relations. Lithos 77, 573–592.[CrossRef][Web of Science]

Canil, D. & Scarfe, C. M. (1990). Phase relations in peridotite + CO2 systems to 12 GPa: implications for the origin of kimberlite and carbonate stability in the Earth's upper mantle. Journal of Geophysical Research 95, 15805–15816.[CrossRef][Web of Science]

Carlson, R. W., Boyd, F. R., Shirey, S. B., Janney, P. E., Grove, T. L., Bowring, S. A., Schmitz, M. D., Dann, J. C., Bell, D. R., Gurney, J. J., Richardson, S. H., Tredoux, M., Menzies, A. H., Pearson, D. G., Hart, R. J., Wilson, A. H. & Moser, D. (2000). Continental growth, preservation, and modification in southern Africa. GSA Today 10, 1–7.

Chalapathi Rao, N. V., Gibson, S. A., Pyle, D. M. & Dickin, A. P. (2004). Petrogenesis of Proterozoic lamproites and kimberlites from the Cuddapah Basin and Dharwar Craton, southern India. Journal of Petrology 45, 907–948.[Abstract/Free Full Text]

Choukroun, M., O'Reilly, S. Y., Griffin, W. L., Pearson, N. J. & Dawson, J. B. (2005). Hf isotopes of MARID (mica–amphibole–rutile–ilmenite–diopside) rutile trace metasomatic processes in the lithospheric mantle. Geology 33, 45–48.[Abstract/Free Full Text]

Clark, T. C. (1994). An integrated geochemical and isotopic study of the Prieska Province kimberlites from the Republic of South Africa. M.Sc. thesis, University of the Witwatersrand.

Clement, C. R. (1982). A comparative geological study of some major kimberlite pipes in the Northern Cape and Orange Free State. Ph.D. thesis, University of Cape Town.

Clement, C. R. & Skinner, M. W. (1979). Textural genetic classification of kimberlites. Transactions of the Geological Society of South Africa 88, 403–410.

Clement, C. R., Skinner, E. M. W. & Scott Smith, B. H. (1984). Kimberlite redefined. Journal of Geology 92, 223–228.[Web of Science]

Clifford, T. N. (1966). Tectono-metallogenic units and metallogenic provinces of Africa. Earth and Planetary Science Letters 1, 241–434.[CrossRef]

Coe, N. (2004). Petrogenesis of the Swartruggens and Star Group II kimberlite dyke swarms, South Africa. M.Sc. thesis, University of Cape Town, 146 pp.

Cook, F. A., van der Velden, A. R., Hall, K. W. & Roberts, B. J. (1999). Frozen subduction in Canada's Northwest Territories: Lithoprobe deep lithospheric reflection profiling of the western Canadian Shield. Tectonics 18, 1–24.[CrossRef][Web of Science]

Cornell, D. H., Hawkesworth, C. J., van Calsteren, P. & Scott, W. D. (1986). Sm–Nd study of Precambrian crustal development in the Prieska–Copperton region, Cape Province. Transactions of the Geological Society of South Africa 89, 17–28.

Dalton, J. A. & Presnall, C. (1998a). The continuum of primary carbonatitic–kimberlitic melt compositions in equilibrium with lherzolite: data from the system CaO–MgO–Al2O3–SiO2–CO2 at 6 GPa. Journal of Petrology 39, 1953–1964.[CrossRef][Web of Science]

Dalton, J. A. & Presnall, D. C. (1998b). Carbonatitic melts along the solidus of model lherzolite in the system CaO–MgO–Al2O3–SiO2–CO2 from 3 to 7 GPa. Contributions to Mineralogy and Petrology 131, 123–136.[CrossRef][Web of Science]

Davis, G. L. (1977). The ages and uranium contents of zircons from kimberlites and associated rocks. Carnegie Institution of Washington, Yearbook 76, 631–635.

de Wit, M. J., Roering, C., Hart, R. J., Armstrong, R. A., de Ronde, C. E. J., Green, R. W. E., Tredoux, M., Peberdy, E. & Hart, R. A. (1992). Formation of an Archean continent. Nature 357, 553–562.[CrossRef]

Duncan, A. R., Erlank, A. J. & Marsh, J. S. (1984). Regional geochemistry of the Karoo Igneous Province. Special Publication of the Geological Society of South Africa 13, 355–388.

Edgar, A. D. & Charbonneau, H. E. (1993). Melting experiments on a SiO2-poor, CaO-rich aphanitic kimberlite from 5–10 GPa and their bearing on sources of kimberlite magmas. American Mineralogist 78, 132–142.[Abstract]

Eggler, D. H. & Wendlandt, R. F. (1979). Experimental studies on the relationship between kimberlite magmas and partial melting of peridotite. In: Boyd, F. R. & Meyer, H. O. A. (eds) Kimberlites, Diatremes and Diamonds; their Geology, Petrology and Geochemistry. Washington, DC: American Geophysical Union, pp. 330–338.

Elliott, T., Plank, T., Zindler, A., White, W. & Bourdon, B. (1997). Element transport from slab to volcanic front at the Mariana arc. Journal of Geophysical Research 102, 14991–15019.[CrossRef]

Erlank, A. J., Waters, F. G., Hawkesworth, C. J., Haggerty, S. E., Allsopp, H. L., Rickard, R. S. & Menzies, M. (1987). Evidence for mantle metasomatism in peridotite nodules from the Kimberley pipes, South Africa. In: Menzies, M. & Hawkesworth, C. J. (eds) Mantle Metasomatism. London: Academic Press: pp. 221–312.

Foley, S. (1992a). Petrological characterisation of the source components of potassic magmas: geochemical and experimental constraints. Lithos 28, 187–204.[CrossRef][Web of Science]

Foley, S. (1992b). Vein-plus-wall-rock melting mechanisms in the lithosphere and the origin of potassic magmas. Lithos 28, 435–453.[CrossRef][Web of Science]

Franz, L., Brey, G. P. & Okrusch, M. (1996). Steady state geotherm, thermal disturbances, and tectonic development of the lower lithosphere underneath the Gibeon Kimberlite Province, Namibia. Contributions to Mineralogy and Petrology 126, 181–198.[CrossRef][Web of Science]

Fraser, K. J. & Hawkesworth, C. J. (1992). The petrogenesis of group 2 ultrapotassic kimberlite from Finsch mine, South Africa. Lithos 28, 327–345.[CrossRef][Web of Science]

Fraser, K. J., Hawkesworth, C. J., Erlank, A. J., Mitchell, R. H. & Scott-Smith, B. H. (1985–1986). Sr, Nd and Pb isotope and minor element geochemistry of lamproites and kimberlites. Earth and Planetary Science Letters 76, 57–70.

Girnis, A. V., Brey, G. P. & Ryabchikov, I. D. (1995). Origin of Group 1A kimberlites: fluid-saturated melting experiments at 45–55 kbar. Earth and Planetary Science Letters 134, 283–296.[CrossRef][Web of Science]

Grégoire, M., Bell, D. R. & le Roex, A. P. (2002). Trace element geochemistry of glimmerite and MARID mantle xenoliths: their relationship to kimberlite and to phlogopite-bearing peridotite revisited. Contributions to Mineralogy and Petrology 142, 603–625.[Web of Science]

Grégoire, M., Bell, D. R. & le Roex, A. P. (2003). Garnet lherzolites from the Kaapvaal Craton (South Africa): trace element evidence for a metasomatic history. Journal of Petrology 44, 629–657.[Abstract/Free Full Text]

Griffin, W. L., Pearson, N. J., Belousova, E., Jackson, S. E., van Achterberg, E., O'Reilly, S. Y. & Shee, S. R. (2000). The Hf isotope composition of cratonic mantle: LAM-MC-ICP-MS analysis of zircon megacrysts in kimberlites. Geochimica et Cosmochimica Acta 64, 133–147.[CrossRef][Web of Science]

Gudfinnsson, G. H. & Presnall, D. C. (2003). Continuous gradations among primary kimberlitic, carbonatitic, melilititic, and komatiitic melts in equilibrium with garnet lherzolite at 3–8 GPa. In: Mitchell, R. & Scott-Smith, B. (eds) VIIIth International Kimberlite Conference, Victoria, BC, Canada, Extended Abstracts.

Gurney, J. J., Helmstaedt, H. H., le Roex, A. P., Richardson, S. H. & Westerlund, K. J. (2005). Diamonds: crustal distribution and formation processes in time and space integrated into an advanced deposit model. Economic Geology (in press).

Haggerty, S. E. (1987). Metasomatic mineral titanates in upper mantle xenoliths. In: Nixon, P. H. (ed.) Mantle Xenoliths. New York: John Wiley, pp. 671–690.

Harris, M., le Roex, A. & Class, C. (2004). Geochemistry of the Uintjiesberg kimberlite, South Africa: petrogenesis of an off-craton, group I, kimberlite. Lithos 74, 149–165.[CrossRef][Web of Science]

Hart, S. R. (1977). The geochemistry and evolution of early Precambrian mantle. Contributions to Mineralogy and Petrology 61, 109–128.[CrossRef][Web of Science]

Hart, S. R. (1988). Heterogeneous mantle domains: signatures, genesis and mixing chronologies. Earth and Planetary Science Letters 90, 273–296.[CrossRef][Web of Science]

Hartnady, C. J. H. & le Roex, A. P. (1985). Southern Ocean hotspot tracks and the Cenozoic absolute motion of the African, Antarctic and South American plates. Earth and Planetary Science Letters 75, 245–257.[CrossRef][Web of Science]

Hawkesworth, C. J., Marsh, J. S., Duncan, A. R., Erlank, A. J. & Norry, M. J. (1984). The role of continental lithosphere in the generation of the Karoo volcanic rocks: evidence from combined Nd- and Sr-isotope studies. Special Publication of the Geological Society of South Africa 13, 341–354.

Helmstaedt, H. & Gurney, J. J. (1984). Kimberlites of Southern Africa—are they related to subduction processes? In: Kornprobst, J. (ed.) Kimberlites and Related Rocks: Kimberlites 1. Amsterdam: Elsevier, pp. 424–435.

Herzberg, C. (1992). Depth and degree of melting of komatiites. Journal of Geophysical Research 97, 4521–4540.

Herzberg, C. & O'Hara, M. J. (2002). Plume-associated ultramafic magmas of Phanerozoic age. Journal of Petrology 43, 1857–1883.[Abstract/Free Full Text]

Hofmann, A. W., Jochum, K. P., Seufert, M. & White, W. M. (1986). Nb and Pb in oceanic basalts: new constraints on mantle evolution. Earth and Planetary Science Letters 79, 33–45.[CrossRef][Web of Science]

Janney, P. E., Carlson, R. W., Shirey, S. B., Bell, D. R. & le Roex, A. P. (2001). Re–Os age and thermal structure of off-craton lithospheric mantle in western South Africa. EOS Transactions, American Geophysical Union 82, 20.

le Roex, A. P. (1985). Geochemistry, mineralogy and magmatic evolution of the basaltic and trachytic lavas from Gough Island, South Atlantic. Journal of Petrology 26, 149–186.[Abstract/Free Full Text]

le Roex, A. P. (1986). Geochemical correlation between southern African kimberlites and South Atlantic hotspots. Nature 324, 243–245.[CrossRef]

le Roex, A. P. & Erlank, A. J. (1982). Quantitative evaluation of fractional crystallization in Bouvet Island lavas. Journal of Volcanology and Geothermal Research 13, 309–338.[CrossRef][Web of Science]

le Roex, A. P., Erlank, A. J. & Needham, H. D. (1981). Geochemical and mineralogical evidence for the occurrence of at least three distinct magma types in the ‘Famous’ region. Contributions to Mineralogy and Petrology 77, 24–37.[CrossRef][Web of Science]

le Roex, A. P., Cliff, R. A. & Adair, B. J. I. (1990). Tristan da Cunha, South Atlantic: geochemistry and petrogenesis of a basanite–phonolite lava series. Journal of Petrology 31, 779–812.[Abstract/Free Full Text]

le Roex, A. P., Späth, A. & Zartman, R. E. (2001). Lithospheric thickness beneath the southern Kenya Rift: implications from basalt geochemistry. Contributions to Mineralogy and Petrology 142, 86–106.

le Roex, A. P., Bell, D. R. & Davis, P. (2003). Petrogenesis of group I kimberlites from Kimberley, South Africa: evidence from bulk-rock geochemistry. Journal of Petrology 44, 2261–2286.[Abstract/Free Full Text]

le Roux, P. J., le Roex, A. P., Schilling, J.-G., Shimizu, N., Perkins, W. W. & Pearce, N. J. G. (2002). Mantle heterogeneity beneath the southern Mid-Atlantic Ridge: trace element evidence for contamination of ambient asthenospheric mantle. Earth and Planetary Science Letters 203, 479–498.[CrossRef][Web of Science]

Maaloe, S. & Aoki, K.-i. (1977). The major element composition of the upper mantle estimated from the composition of lherzolites. Contributions to Mineralogy and Petrology 63, 161–173.[CrossRef][Web of Science]

Marsh, J. S., Hooper, P. R., Rehacek, J., Duncan, R. A. & Duncan, A. R. (1997). Stratigraphy and Age of Karoo Basalts of Lesotho and Implications for Correlations within the Karoo Igneous Province. Washington, DC: American Geophysical Union.

Mathias, M., Siebert, J. C. & Rickwood, P. C. (1970). Some aspects of the mineralogy and petrology of ultramafic xenoliths in kimberlite. Contributions to Mineralogy and Petrology 26, 75–123.[CrossRef][Web of Science]

McKenzie, D. & O'Nions, R. K. (1983). Mantle reservoirs and ocean island basalts. Nature 301, 229–231.[CrossRef]

Menzies, M., Rogers, N., Tindle, A. & Hawkesworth, C. J. (1987). Metasomatic and enrichment processes in lithospheric peridotites, and effect of asthenosphere–lithosphere interaction. In: Menzies, M. & Hawkesworth, C. J. (eds) Mantle Metasomatism. London: Academic Press, pp. 313–361.

Mitchell, R. (1995). Kimberlites, Orangites and Related Rocks. New York: Plenum.

Mitchell, R. H. (2004). Experimental studies at 5–12 GPa of the Ondermatjie hypabyssal kimberlite. Lithos 76, 551–564.[CrossRef][Web of Science]

Nixon, P. H., Boyd, F. R. & Boctor, N. Z. (1983). East Griqualand kimberlites. Transactions of the Geological Society of South Africa 86, 221–236.

Nowell, G. M., Pearson, D. G., Kempton, P. D., Noble, S. R. & Smith, C. B. (1999). Origins of kimberlites: a hafnium perspective. In: Gurney, J. J., Gurney, J. L., Pascoe, M. D. & Richardson, S. H. (eds) Proceedings of the VIIth International Kimberlite Conference. Cape Town: Red Roof Design, pp. 616–624.

Nowell, G. M., Pearson, D. G., Bell, D. R., Carlson, R. W., Smith, C. B., Kempton, P. D. & Noble, S. R. (2004). Hf isotope systematics of kimberlites and their megacrysts: new constraints on their source regions. Journal of Petrology 45, 1583–1612.[Abstract/Free Full Text]

O'Connor, J. M. & le Roex, A. P. (1992). South Atlantic hot spot-plume systems: 1. Distribution of volcanism in time and space. Earth and Planetary Science Letters 113, 343–364.[CrossRef][Web of Science]

O'Nions, R. K. & Pankhurst, R. J. (1974). Petrogenetic significance of isotope and trace element variations in volcanic rocks from the Mid-Atlantic. Journal of Petrology 15, 603–634.[Abstract/Free Full Text]

O'Nions, R. K., Hamilton, P. J. & Evensen, N. M. (1977). Variations in 143Nd/144Nd and 87Sr/86Sr ratios in oceanic basalts. Earth and Planetary Science Letters 34, 13–22.[CrossRef][Web of Science]

Pearson, D. G. & Nowell, G. M. (2002). The continental lithospheric mantle: characteristics and significance as a mantle reservoir. Philosophical Transactions of the Royal Society of London, Series A 360, 2383–2410.

Pearson, D. G., Rogers, N. W., Irving, A. J., Smith, C. B. & Hawkesworth, C. J. (1995). Source regions of kimberlites and lamproites: constraints from Re–Os isotopes. In: VIth International Kimberlite Conference, Extended Abstracts. Novosibirsk: Russian Academy of Sciences, United Institute of Geology, Geophysics and Mineralogy, pp. 430–432.

Pearson, D. G., Nowell, G. M., Dowall, D. P., Kjarsgaard, B. A., Kopylova, M. G. & Armstrong, J. A. (2003). The relative roles of lithosphere and convecting mantle in kimberlites from the Slave Province, NWT: constraints from Re–Os isotopes and olivine population studies. VIIIth International Kimberlite Conference, Extended Abstracts.

Phillips, D., Harris, J. W. & Viljoen, K. S. (2004). Mineral chemistry and thermobarometry of inclusions for De Beers pool diamonds, Kimberley, South Africa. Lithos 77, 155–179.[CrossRef][Web of Science]

Pretorious, D. A. (1974). The structural boundary between the Kaapvaal craton and the Sanama crustal provinces. In: Information Circular of Economic Geology Research Unit. Johannesburg: University of the Witwatersrand, 26 pp.

Price, S. E., Russell, J. K. & Kopylova, M. G. (2000). Primitive magma from the Jericho Pipe, N.W.T., Canada: constraints on primary kimberlite melt chemistry. Journal of Petrology 41, 789–808.[Abstract/Free Full Text]

Robey, J. (1981). Kimberlites in the Central Cape Province, R.S.A. Ph.D. thesis, University of Cape Town.

Rogers, N. W., Hawkesworth, C. J. & Palacz, Z. A. (1992). Phlogopite in the generation of olivine-melilitites from Namaqualand, South Africa, and its implications for element fractionation processes in the upper mantle. Lithos 28, 347–365.[CrossRef][Web of Science]

Schmidt, K. H., Bottazzi, P., Vannucci, R. & Mengel, K. (1999). Trace element partitioning between phlogopite, clinopyroxene and leucite lamproite melt. Earth and Planetary Science Letters 168, 287–300.[CrossRef][Web of Science]

Schmitz, M. D., Bowring, S. A., de Wit, M. J. & Gartz, V. (2004). Subduction and terrane collision stabilize the western Kaapvaal craton tectosphere 2·9 billion years ago. Earth and Planetary Science Letters 222, 363–376.

Shee, S. R. (1985). The petrogenesis of the Wesselton Mine Kimberlite, Kimberley, Cape Province, R.S.A. Ph.D. thesis, University of Cape Town.

Simon, N. S. C., Irvine, G. J., Davies, G. R., Pearson, D. G. & Carlson, R. W. (2003). The origin of garnet and clinopyroxene in ‘depleted’ Kaapvaal peridotites. Lithos 71, 289–322.[CrossRef][Web of Science]

Skinner, E. M. W. (1989). Contrasting Group I and Group II kimberlite petrology: towards a genetic model for kimberlites. In: Ross, J., Jacques, A. L., Ferguson, J., Green, D. H., O'Reilly, S. Y., Danchin, R. V. & Janse, A. J. A. (eds) IVth International Kimberlite Conference. Perth: Geological Society of Australia, pp. 528–544.

Skinner, E. M. W., Smith, C. B., Viljoen, K. S. & Clark, T. C. (1992). The petrography, tectonic setting and emplacement ages of kimberlites in the south western border region of the Kaapvaal Craton, Prieska area, South Africa. In: Meyer, H. O. A. & Leonardos, O. H. (eds) Kimberlites, Related Rocks and Mantle Xenoliths. Brasilia: CPRM, pp. 80–97.

Smith, C. B. (1983). Pb, Sr and Nd isotopic evidence for sources of southern African Cretaceous kimberlites. Nature 304, 51–54.[CrossRef]

Smith, C. B., Allsopp, H. L., Kramers, J. D., Hutchinson, G. & Roddick, J. C. (1985a). Emplacement ages of Jurassic–Cretaceous South African kimberlites by the Rb–Sr method on phlogopite and whole-rock samples. Transactions of the Geological Society of South Africa 88, 249–266.

Smith, C. B., Gurney, J. J., Skinner, E. M. W., Clement, C. R. & Ebrahim, N. (1985b). Geochemical character of the southern African kimberlites: a new approach based on isotopic constraints. Transactions of the Geological Society of South Africa 88, 267–280.

Smith, C. B., Clark, T. C., Barton, E. S. & Bristow, J. W. (1994). Emplacement ages of kimberlite occurrences in the Prieska region, southwest border of the Kaapvaal Craton, South Africa. Chemical Geology 113, 149–169.[CrossRef][Web of Science]

Späth, A., le Roex, A. P. & Opiyo-Akech, N. (2001). Plume–lithosphere interaction and the origin of continental rift-related alkaline volcanism—the Chyulu Hills Volcanic Province, southern Kenya. Journal of Petrology 42, 765–787.[Abstract/Free Full Text]

Sun, S.-s. & McDonough, W. F. (1989). Chemical and isotopic systematics of oceanic basalts: implications for mantle composition and processes. In: Saunders, A. D. & Norry, M. J. (eds) Magmatism in the Ocean Basins. Geological Society, London, Special Publications 42, 313–345.

Tainton, K. M. (1992). The petrogenesis of group-2 kimberlites and lamproites from the Northern Cape Province, South Africa. Ph.D. thesis, Cambridge.

Tainton, K. M. & McKenzie, D. (1994). The generation of kimberlites, lamproites, and their source rocks. Journal of Petrology 35, 787–817.[Abstract/Free Full Text]

Thomas, R. J., Agenbacht, A. L. D., Cornell, D. H. & Moore, J. M. (1994). The Kibaren of southern Africa: tectonic evolution and metallogeny. Ore Geology Reviews 9, 131–160.

Ulmer, P. & Sweeney, R. J. (2002). Generation and differentiation of group II kimberlites: constraints from a high-pressure experimental study to 10 GPa. Geochimica et Cosmochimica Acta 66, 2139–2153.[CrossRef][Web of Science]

van Achterbergh, E., Griffin, W. L. & Stiefenhofer, J. (2001). Metasomatism in mantle xenoliths from the Letlhakane kimberlites: estimation of element fluxes. Contributions to Mineralogy and Petrology 141, 397–414.[Web of Science]

Wagner, P. A. (1914). The Diamond Fields of Southern Africa. Johannesburg: Transvaal Leader.

Walker, R. J., Carlson, R. W., Shirey, S. B. & Boyd, F. R. (1989). Os, Sr, Nd, and Pb isotope systematics of southern African peridotite xenoliths; implications for the chemical evolution of subcontinental mantle. Geochimica et Cosmochimica Acta 53, 1583–1595.[CrossRef][Web of Science]

Wass, S. Y. & Rogers, N. W. (1980). Mantle metasomatism—precursor to continental alkaline volcanism. Geochimica et Cosmochimica Acta 44, 1811–1823.[CrossRef][Web of Science]

Waters, F. G. & Erlank, A. J. (1988). Assessment of the vertical extent and distribution of mantle metasomatism below Kimberley, South Africa. Journal of Petrology, 185–204.

Weaver, B. L. (1991). Trace element evidence for the origin of ocean-island basalts. Geology 19, 123–126.[Abstract/Free Full Text]

Wedepohl, K. H. & Muramatsu, Y. (1975). The Chemical Composition of Kimberlites Compared with the Average Composition of Three Basaltic Magma Types. Washington, DC: American Geophysical Union.

Williams, H. M., Turner, S. P., Pearce, J. A., Kelley, S. P. & Harris, N. B. W. (2004). Nature and source regions for post-collisional, potassic magmas in southern and northern Tibet from geochemical variation and inverse trace element modelling. Journal of Petrology 45, 555–607.[Abstract/Free Full Text]

Willis, J. P. (1999). Instrumental Parameters and Data Quality for Routine Major and Trace Element Determinations by WDXRFS. Cape Town: University of Cape Town.

Woodhead, J. D., Eggins, S. M. & Johnson, R. W. (1998). Magma genesis in the New Britain Island Arc: further insights into melting and mass transfer process. Journal of Petrology 39, 1641–1668.[CrossRef][Web of Science]

Yamashita, H., Arima, M. & Ohtani, E. (1995). High pressure melting experiments on group II kimberlite up to 8 GPa: implications for mantle metasomatism. In: VIth International Kimberlite Conference, Extended Abstracts. Novosibirsk: Russian Academy of Sciences, United Institute of Geology, Geophysics and Mineralogy, pp. 669–671.

Zindler, A. (1980). Geochemical processes in the Earth's mantle and the nature of crust–mantle interactions: evidence from studies of Nd and Sr isotope ratios in mantle derived igneous rocks and lherzolite nodules. Ph.D. thesis, MIT, Cambridge, MA.


Add to CiteULike CiteULike   Add to Connotea Connotea   Add to Del.icio.us Del.icio.us    What's this?


This article has been cited by other articles:


Home page
J PetrologyHome page
N. T. Arndt, M. Guitreau, A.-M. Boullier, A. Le Roex, A. Tommasi, P. Cordier, and A. Sobolev
Olivine, and the Origin of Kimberlite
J. Petrology, March 1, 2010; 51(3): 573 - 602.
[Abstract] [Full Text] [PDF]


Home page
J PetrologyHome page
R. G. Gonzaga, M. A. Menzies, M. F. Thirlwall, D. E. Jacob, and A. Leroex
Eclogites and Garnet Pyroxenites: Problems Resolving Provenance Using Lu-Hf, Sm-Nd and Rb-Sr Isotope Systems
J. Petrology, January 6, 2010; (2010) egp091v1.
[Abstract] [Full Text] [PDF]


Home page
Can MineralHome page
N. V. C. Rao and A. N. Dongre
MINERALOGY AND GEOCHEMISTRY OF KIMBERLITES NK-2 AND KK-6, NARAYANPET KIMBERLITE FIELD, EASTERN DHARWAR CRATON, SOUTHERN INDIA: EVIDENCE FOR A TRANSITIONAL KIMBERLITE SIGNATURE
Can Mineral, October 1, 2009; 47(5): 1117 - 1135.
[Abstract] [Full Text] [PDF]


Home page
J PetrologyHome page
S. A. Gibson, J. Malarkey, and J. A. Day
Melt Depletion and Enrichment beneath the Western Kaapvaal Craton: Evidence from Finsch Peridotite Xenoliths
J. Petrology, October 22, 2008; (2008) egn048v1.
[Abstract] [Full Text] [PDF]


Home page
Can MineralHome page
T. C. Birkett
FIRST-ROW TRANSITION ELEMENTS, Y AND Ga IN KIMBERLITE AND LAMPROITE: APPLICATIONS TO DIAMOND PROSPECTIVITY AND PETROGENESIS
Can Mineral, October 1, 2008; 46(5): 1269 - 1282.
[Abstract] [Full Text] [PDF]


Home page
Can MineralHome page
T. F.D. Nielsen and K. K. Sand
THE MAJUAGAA KIMBERLITE DIKE, MANIITSOQ REGION, WEST GREENLAND: CONSTRAINTS ON AN Mg-RICH SILICOCARBONATITIC MELT COMPOSITION FROM GROUNDMASS MINERALOGY AND BULK COMPOSITIONS
Can Mineral, August 1, 2008; 46(4): 1043 - 1061.
[Abstract] [Full Text] [PDF]


Home page
J PetrologyHome page
G. P. Brey, V. K. Bulatov, A. V. Girnis, and Y. Lahaye
Experimental Melting of Carbonated Peridotite at 6-10 GPa
J. Petrology, April 1, 2008; 49(4): 797 - 821.
[Abstract] [Full Text] [PDF]


Home page
Can MineralHome page
L. Melluso, M. Lustrino, E. Ruberti, P. Brotzu, C. de Barros Gomes, L. Morbidelli, V. Morra, D. P. Svisero, and F. d'Amelio
MAJOR- AND TRACE-ELEMENT COMPOSITION OF OLIVINE, PEROVSKITE, CLINOPYROXENE, Cr-Fe-Ti OXIDES, PHLOGOPITE AND HOST KAMAFUGITES AND KIMBERLITES, ALTO PARANAIBA, BRAZIL
Can Mineral, February 1, 2008; 46(1): 19 - 40.
[Abstract] [Full Text] [PDF]


Home page
South African Journal of GeologyHome page
M. Merry and A. le Roex
Megacryst suites from the Lekkerfontein and Uintjiesberg kimberlites, southern Africa: evidence for a non-cognate origin
South African Journal of Geology, December 1, 2007; 110(4): 597 - 610.
[Abstract] [Full Text] [PDF]


Home page
South African Journal of GeologyHome page
M. Becker, A. P. le Roex, and C. Class
Geochemistry and petrogenesis of South African transitional kimberlites located on and off the Kaapvaal Craton
South African Journal of Geology, December 1, 2007; 110(4): 631 - 646.
[Abstract] [Full Text] [PDF]


Home page
GeologyHome page
C. Paton, J. M. Hergt, D. Phillips, J. D. Woodhead, and S. R. Shee
New insights into the genesis of Indian kimberlites from the Dharwar Craton via in situ Sr isotope analysis of groundmass perovskite
Geology, November 1, 2007; 35(11): 1011 - 1014.
[Abstract] [Full Text] [PDF]


Home page
South African Journal of GeologyHome page
M. de Wit
The Kalahari Epeirogeny and climate change: differentiating cause and effect from core to space
South African Journal of Geology, September 1, 2007; 110(2-3): 367 - 392.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrow Supplementary data
Right arrow All Versions of this Article:
47/4/673    most recent
egi089v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrow Search for citing articles in:
ISI Web of Science (41)
Right arrowRequest Permissions
Google Scholar
Right arrow Articles by BECKER, M.
Right arrow Articles by ROEX, A. P. L.
Right arrow Search for Related Content
GeoRef
Right arrow GeoRef Citation
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?