Multi-stage garnet formation and destruction in Kimberley harzburgitic xenoliths, South Africa

Abstract

Thirty-nine garnet harzburgites from Kimberley in the Kaapvaal Craton (South Africa) were studied to constrain the origin, age and evolution of sub-cratonic lithospheric mantle (SCLM). In order to avoid chemical overprinting by recent metasomatism, only garnet harzburgites that appeared clinopyroxene-free to the naked eye were sampled. The majority of garnets were, however, in equilibrium with clinopyroxene (24 of 39). Whole rock and mineral major–trace element geochemistry and garnet Sr–Nd–Hf isotope data are presented. Equilibration pressures range from 3.8–6.1 GPa, indicating the harzburgites were derived from a large portion of the SCLM (~115–185 km). High olivine Mg# (~93.4, n = 39) and low whole rock heavy rare earth elements (HREE) contents are consistent with large degrees of partial melting (>45%) and garnet exhaustion leaving a dunitic residue with olivine ≥90%, orthopyroxene ≤10% and HREE <0.01 times chondrite. Mineral modes, whole rock Al2O3 (0.5–3.2 wt%) and SiO2 (43.1–49.1 wt%), however, indicate heterogeneous re-introduction of garnet (≤13%) and orthopyroxene (≤50%). Harzburgites with high garnet and relatively low orthopyroxene modes (mostly ~7–13% and ~ 9–30%; n = 6) are characterised by mildly sinusoidal garnet REE patterns (Tb-Dy minimum and high HREE) and Archaean depleted Hf TDM ages (2.7–3.3 Ga; εHfe: +190 to +709). In contrast, harzburgites with high orthopyroxene and relatively low garnet and modes (~1.5–7.5% and ~ 25–50%; n = 19) are characterised by highly sinuous REE patterns (Ho-Yb minimum and low HREE) and Proterozoic enriched Hf TDM ages (0.7–1.6 Ga; εHfe: −16 to +6). It is inferred that Archaean G10 garnet re-introduction caused a significant increase in HREE, making melt depletion models based on HREE inaccurate. Orthopyroxene addition, a few hundred million years later, most likely at ~2.7 Ga and associated with Ventersdorp magmatic activity, caused partial consumption of garnet and olivine, and changed garnet compositions leading to: 1) Cr/Al ratio increase; 2) HREE decrease; 3) more sinusoidal REE patterns; and 4) un-radiogenic 176Hf/177Hf. Garnets define a Lu-Hf isochron age of 2702 ± 64 Ma (εHfi = +44, n = 31), which is interpreted as a consequence of partial isotopic equilibrium within the SCLM and mixing of the garnet- and orthopyroxene-rich metasomatic components. The low LILE contents and absence of Nb-Ta anomalies are consistent with modal metasomatism caused by intra-plate magmatism. In addition, the REE signatures of metasomatic agents in equilibrium with the garnets suggest that carbonatitic melts and SiO2-rich hydrous melts were responsible for re-introduction of garnet and orthopyroxene, respectively. Sr-Nd isotope systematics were disrupted associated with kimberlite magmatism (Nd isochron: 217 ± 58 Ma, εNdi = +4; n = 34), consistent with recent G10 garnet transformation into G9 garnets (Ca + Fe-enriched). This event may have caused garnet addition (up to 1%), suggesting that garnet was formed or destroyed in at least 4 different events: i) initial extensive polybaric melting, ii) asthenospheric melts re-introducing the bulk of the garnet, iii) orthopyroxene addition and garnet loss, all in the Archaean, and iv) minor garnet addition possibly related to recent kimberlite magmatism prior to eruption.