Petrology and U–Pb geochronology of lower crustal xenoliths and the development of a craton, Slave Province, Canada

Lower crustal xenoliths recovered from Eocene to Cambrian kimberlites in the central and southern Slave craton are dominated by mafic granulites (garnet, clinopyroxene, plagioclase±orthopyroxene), with subordinate metatonalite and peraluminous felsic granulites. Geothermobarometry indicates metamorphic conditions of 650–800 °C at pressures of 0.9–1.1 GPa. The metamorphic conditions are consistent with temperatures expected for the lower crust of high-temperature low-pressure (HT-LP) metamorphic belts characteristic of Neoarchean metamorphism in the Slave craton. U–Pb geochronology of zircon, rutile and titanite demonstrate a complex history in the lower crust. Mesoarchean protoliths occur beneath the central Slave supporting models of an east-dipping boundary between Mesoarchean crust in the western and Neoarchean crust in the eastern Slave. At least, two episodes of igneous and metamorphic zircon growth occurred in the interval 2.64–2.58 Ga that correlate with the age of plutonism and metamorphism in the upper crust, indicating magmatic addition to the lower crust and metamorphic reworking during this period. In addition, discrete periods of younger zircon growth at ca. 2.56–2.55 and 2.51 Ga occurred 20–70 my after the cessation of ca. 2.60–2.58 Ga regional HT-LP metamorphism and granitic magmatism in the upper crust. This pattern of younger metamorphic events in the deep crust is characteristic of the Slave as well as other Archean cratons (e.g., Superior). The high temperature of the lower crust immediately following amalgamation of the craton, coupled with evidence for continued metamorphic zircon growth for >70 my after ‘stabilization’ of the upper crust, is difficult to reconcile with a thick (200 km), cool lithospheric mantle root beneath the craton prior to this event. We suggest that thick tectosphere developed synchronously or after these events, most likely by imbrication of mantle beneath the craton at or after ca. 2.6 Ga. The minimum age for establishing a cratonic like geotherm is given by lower crustal rutile ages of ca. 1.8 Ga in the southern Slave. Transient heating and possible magmatic additions to the lower crust continued through the Proterozoic, with possible additional growth of the tectosphere.     

The electrical structure of the Slave craton

The Slave craton in northwestern Canada, a relatively small Archean craton (600×400 km), is ideal as a natural laboratory for investigating the formation and evolution of Mesoarchean and Neoarchean sub-continental lithospheric mantle (SCLM). Excellent outcrop and the discovery of economic diamondiferous kimberlite pipes in the centre of the craton during the early 1990s have led to an unparalleled amount of geoscientific information becoming available.

Over the last 5 years deep-probing electromagnetic surveys were conducted on the Slave, using the natural-source magnetotelluric (MT) technique, as part of a variety of programs to study the craton and determine its regional-scale electrical structure. Two of the four types of surveys involved novel MT data acquisition; one through frozen lakes along ice roads during winter, and the second using ocean-bottom MT instrumentation deployed from float planes.

The primary initial objective of the MT surveys was to determine the geometry of the topography of the lithosphere–asthenosphere boundary (LAB) across the Slave craton. However, the MT responses revealed, completely serendipitously, a remarkable anomaly in electrical conductivity in the SCLM of the central Slave craton. This Central Slave Mantle Conductor (CSMC) anomaly is modelled as a localized region of low resistivity (10–15 Ω m) beginning at depths of 80–120 km and striking NE–SW. Where precisely located, it is spatially coincident with the Eocene-aged kimberlite field in the central part of the craton (the so-called “Corridor of Hope”), and also with a geochemically defined ultra-depleted harzburgitic layer interpreted as oceanic or arc-related lithosphere emplaced during early tectonism. The CSMC lies wholly within the NE–SW striking central zone defined by Grütter et al. [Grütter, H.S., Apter, D.B., Kong, J., 1999. Crust–mantle coupling; evidence from mantle-derived xenocrystic garnets. Contributed paper at: The 7th International Kimberlite Conference Proceeding, J.B. Dawson Volume, 1, 307–313] on the basis of garnet geochemistry (G10 vs. G9) populations.

Deep-probing MT data from the lake bottom instruments infer that the conductor has a total depth-integrated conductivity (conductance) of the order of 2000 Siemens, which, given an internal resistivity of 10–15 Ω m, implies a thickness of 20–30 km. Below the CSMC the electrical resistivity of the lithosphere increases by a factor of 3–5 to values of around 50 Ω m. This change occurs at depths consistent with the graphite–diamond transition, which is taken as consistent with a carbon interpretation for the CSMC.

Preliminary three-dimensional MT modelling supports the NE–SW striking geometry for the conductor, and also suggests a NW dip. This geometry is taken as implying that the tectonic processes that emplaced this geophysical–geochemical body are likely related to the subduction of a craton of unknown provenance from the SE (present-day coordinates) during 2630–2620 Ma. It suggests that the lithospheric stacking model of Helmstaedt and Schulze [Helmstaedt, H.H., Schulze, D.J., 1989. Southern African kimberlites and their mantle sample: implications for Archean tectonics and lithosphere evolution. In Ross, J. (Ed.), Kimberlites and Related Rocks, Vol. 1: Their Composition, Occurrence, Origin, and Emplacement. Geological Society of Australia Special Publication, vol. 14, 358–368] is likely correct for the formation of the Slave's current SCLM.     

Two anisotropic layers in the Slave craton

Four years of recording global earthquakes using a broadband seismometer located at the Ekati diamond mine revealed variations with earthquake azimuth in the arrival of SKS phases. These variations can be modeled assuming two distinct layers of anisotropy in the lithosphere. The lower layer probably lies in the mantle, and the anisotropy aligns with both North American plate motion and the strike of mantle structures identified by previous conductivity and geochemical analyses, at ˜N50°E. The upper layer is hypothesized to result from regional structures in the uppermost mantle and the crust; these trends are distinct from the mantle trends.     

Multiple-mineral inclusions in diamonds from the Snap Lake/King Lake kimberlite dike, Slave craton, Canada: a trace-element perspective

Multiple inclusions of minerals in diamonds from the Snap Lake/King Lake kimberlites of the southeastern Slave craton in Canada have been analyzed for trace elements to elucidate the petrogenetic history of these inclusions, and of their host diamonds. As observed worldwide, the harzburgitic-garnet diamond inclusions (DIs) possess sinusoidal REE patterns that indicate an early depletion event, followed by metasomatism by LREE-enriched, HREE-depleted fluids. Furthermore, these fluids appear to contain appreciable concentrations of LILE and HFSE, based on the increasing abundances of these elements in the olivine inclusion that occurs at the outer portion of a diamond compared to that near the core. The compositions of these fluids are probably a mixture of hydrous-silicic melt, carbonatitic melt, and brine, similar to the compositions of micro-inclusions in diamonds reported by Navon et al. (2003). Comparison between the compositions of majoritic and normal harzburgitic garnets shows that the former are more depleted in terms of major/minor elements (higher Cr#) but significantly more enriched in the REE (up to 10×). This characteristic may indicate the higher susceptibility for metasomatic enrichment of previously more depleted garnets. Garnets of eclogitic paragenesis show strong LREE-depleted patterns, whereas the coexisting omphacite inclusion has relatively flat light- and middle-REE but depleted HREE. Whole-rock reconstruction from coexisting garnet and omphacite inclusions indicates that the protolith of these inclusions was probably the extrusive section of an oceanic crust, subducted beneath the Slave craton.     

华北克拉通太古宙地壳演化和最古老的岩石

简要总结了华北克拉通太古宙地壳演化的规律、认识及存在的问题,特别是近年来在最古老的岩石和锆石年代学研究方面取得的进展。华北克拉通构造热事件主要发生在约2.5Ga,导致大规模陆壳的形成。仅在鲁西、胶东等少数地区有较大规模的约2.7Ga的地质体存在。古元古代晚期的陆陆碰撞使华北克拉通最终形成统一的整体。最近对鞍山地区的研究发现大量3.6～3.8Ga的岩石和锆石。在冀东、信阳、焦作及其他地区也有始太古代—古太古代的岩石和锆石存在。华北克拉通可能存在几个不同的古太古代—始太古代陆核。     

Mantle Xenoliths from the Southeastern Slave Craton: Evidence for Chemical Zonation in a Thick, Cold Lithosphere

We present the first data on the petrology of the mantle lithosphereof the Southeastern (SE) Slave craton, Canada. These are basedon petrographic, mineralogical and geochemical studies of mantlexenoliths in Pipe 5034 of the Cambrian Gahcho Kué kimberlitecluster. Major types of mantle xenoliths include altered eclogite,coarse garnet or spinel peridotite, and deformed garnet peridotite.The peridotites belong to the low-temperature suite and formedat T=600–1300°C and P= 25–80 kbar in a thick(at least 220–250 km), cool lithosphere. The SE Slavemantle is cooler than the mantle of other Archaean cratons andthat below other terranes of the Slave craton. The thick lithosphereand the relatively cool thermal regime provide favourable conditionsfor formation and preservation of diamonds beneath the SE Slaveterrane. Similar to average Archaean mantle worldwide, the SESlave peridotite is depleted in magmaphile major elements andcontains olivine with forsterite content of 91–93·5.With respect to olivine composition and mode, all terranes ofthe Slave mantle show broadly similar compositions and are relativelyorthopyroxene-poor compared with those of the Kaapvaal and Siberiancratons. The SE Slave spinel peridotite is poorer in Al, Caand Fe, and richer in Mg than deeper garnet peridotite. Thegreater chemical depletion of the shallow upper mantle is typicalof all terranes of the Slave craton and may be common for thesubcontinental lithospheric peridotitic mantle in general. Peridotiticxenoliths of the SE Slave craton were impregnated by kimberliticfluids that caused late-stage recrystallization of primary clinopyroxene,spinel, olivine and spinel-facies orthopyroxene, and formationof interstitial clinopyroxene. This kimberlite-related recrystallizationdepleted primary pyroxenes and spinel in Al. The kimberliticfluid was oxidizing, Ti-, Fe- and K-rich, and Na-poor, and introducedserpentine, chlorite, phlogopite and spinel into peridotitesat P < 35 kbar. KEY WORDS: kimberlite xenolith; lithosphere; mantle terrane; chemical zoning; thermobarometry; Slave craton     

Peridotitic mantle xenoliths from kimberlites on the Ekati Diamond Mine property, N.W.T., Canada: major element compositions and implications for the lithosphere beneath the central Slave craton

The composition, structure and thermal state of the lithosphere beneath the Slave craton have been studied by analysing over 300 peridotitic mantle xenoliths or multiphase xenocrysts entrained within kimberlites in the Lac de Gras area. These xenoliths are derived from seven kimberlites located on the Ekati Diamond Mine™ property and define a detailed stratigraphic profile through the central Slave lithosphere from less than 120 km down to 200 km. Two dominant peridotite types are present, namely garnet-bearing harzburgite and lherzolite with rare occurrences of chromite-facies peridotite, websterite and wehrlite. The pressures and temperatures (P–T's) defined by the entire data-set range from 28 to 62 kbar and 650 to 1250 °C, respectively, and approximately intersect the diamond stability field at 900 °C and 42 kbar. There is no apparent change in the geotherm with depth that is discernable beyond the resolution of the various thermobarometers. The peridotites can be divided into two compositional zones—a shallow layer dominated by garnet harzburgite that straddles the diamond–graphite boundary and a deeper layer that is strongly dominated by garnet lherzolite. Compositionally, the harzburgites (and to a lesser extent, the shallow lherzolites) are ultra-depleted relative to the more fertile deeper layer, irrespective of whether they reside within the graphite or diamond stability field. This ultra-depleted layer beneath Ekati continues to 150 km.     

Petrological constraints on seismic properties of the Slave upper mantle (Northern Canada)

Modes and compositions of minerals in Slave mantle xenoliths, together with their pressures and temperatures of equilibrium were used to derive model depth profiles of P- and S-wave velocities (Vp, Vs) for composites equivalent to peridotite, pyroxenite and eclogite. The rocks were modeled as isotropic aggregates with uniform distribution of crystal orientations, based on single-crystal elastic moduli and volume fractions of constituent minerals. Calculated seismic wave velocities are adjusted for in situ pressure and temperature conditions using (1) experimental P- and T- derivatives for bulk rocks' Vp and Vs, and (2) calculated P- and T- derivatives for bulk rocks' elastic moduli and densities. The peridotite seismic profiles match well with the globally averaged IASP91 model and with seismic tomography results for the Slave mantle. In peridotite, an observed increase of seismic wave velocities with depth is controlled by lower degrees of chemical depletion in the deeper upper mantle. In eclogite, seismic velocities increase more rapidly with depth than in peridotite. This follows from contrasting first-order pressure derivatives of bulk isotropic moduli for eclogite and peridotite, and from the lower compressibility of eclogite at high pressures. Our calculations suggest that depletion in cratonic mantle has a distinct seismic signature compared to non-cratonic mantle. Depleted mantle on cratons should have slower Vp, faster Vs and should show lower Poisson's ratios due to an orthopyroxene enrichment. For the modelled Slave craton xenoliths, the predicted effect on seismic wave velocities would be up to 0.05 km/s.     

Crustal structure of the Kaapvaal craton and its significance for early crustal evolution

High-quality seismic data obtained from a dense broadband array near Kimberley, South Africa, exhibit crustal reverberations of remarkable clarity that provide well-resolved constraints on the structure of the lowermost crust and Moho. Receiver function analysis of Moho conversions and crustal multiples beneath the Kimberley array shows that the crust is 35 km thick with an average Poisson's ratio of 0.25. The density contrast across the Moho is 15%, indicating a crustal density about 2.86 gm/cc just above the Moho, appropriate for felsic to intermediate rock compositions. Analysis of waveform broadening of the crustal reverberation phases suggests that the Moho transition can be no more than 0.5 km thick and the total variation in crustal thickness over the 2400 km² footprint of the array no more than 1 km. Waveform and travel time analysis of a large earthquake triggered by deep gold mining operations (the Welkom mine event) some 200 km away from the array yield an average crustal thickness of 35 km along the propagation path between the Kimberley array and the event. P- and S-wave velocities for the lowermost crust are modeled to be 6.75 and 3.90 km/s, respectively, with uppermost mantle velocities of 8.2 and 4.79 km/s, respectively. Seismograms from the Welkom event exhibit theoretically predicted but rarely observed crustal reverberation phases that involve reflection or conversion at the Moho. Correlation between observed and synthetic waveforms and phase amplitudes of the Moho reverberations suggests that the crust along the propagation path between source and receiver is highly uniform in both thickness and average seismic velocity and that the Moho transition zone is everywhere less than about 2 km thick. While the extremely flat Moho, sharp transition zone and low crustal densities beneath the region of study may date from the time of crustal formation, a more geologically plausible interpretation involves extensive crustal melting and ductile flow during the major craton-wide Ventersdorp tectonomagmatic event near the end of Archean time.     

Multi-Stage Modification of the Northern Slave Mantle Lithosphere: Evidence from Zircon- and Diamond-Bearing Eclogite Xenoliths Entrained in Jericho Kimberlite, Canada

The Jericho kimberlites are part of a small Jurassic kimberlitecluster in the northern Slave craton, Canada. A variety of datingtechniques were applied to constrain the nature and age of twoJericho kimberlites, JD-1 (170·2 ± 4·3Ma Rb–Sr phlogopite megacrysts, 172·8 ±0·7 Ma U–Pb eclogite rutile, 178 ± 5 MaU–Pb eclogite zircon lower intercept) and JD-3 (173 ±2 Ma Rb–Sr phlogopite megacryst; 176·6 ±3·2 Ma U–Pb perovskite), and all yielded identicalresults within analytical uncertainty. As there is no discernibledifference in the radiometric ages obtained for these two pipes,the composite Rb–Sr phlogopite megacryst date of 173·1± 1·3 Ma is interpreted as the best estimate forthe emplacement age of both Jericho pipes. The initial Sr isotopecomposition of 0·7053 ± 0·0003 derivedfrom phlogopite megacrysts overlaps the range (0·7043–0·7084)previously reported for Jericho whole-rocks. These strontiumisotope data, combined with the radiogenic initial ²⁰⁶Pb/²⁰⁴Pbratio of 18·99 ± 0·33 obtained in thisstudy, indicate that the Jericho kimberlites are isotopicallysimilar to Group 1 kimberlites as defined in southern Africa.The Jericho kimberlites are an important new source of mantlexenoliths that hold clues to the nature of the Slave cratonsubcontinental mantle. A high proportion (30%) of the Jerichomantle xenolith population consists of various eclogite typesincluding a small number (2–3%) of apatite-, diamond-,kyanite- and zircon-bearing eclogites. The most striking aspectof the Jericho zircon-bearing eclogite xenoliths is their peculiargeochemistry. Reconstructed whole-rock compositions indicatethat they were derived from protoliths with high FeO, Al₂O₃and Na₂O contents, reflected in the high-FeO (22·6–27·5wt %) nature of garnet and the high-Na₂O (8·47–9·44wt %) and high-Al₂O₃ (13·12–14·33 wt %)character of the clinopyroxene. These eclogite whole-rock compositionsare highly enriched in high field strength elements (HFSE) suchas Nb (133–1134 ppm), Ta (5–28 ppm), Zr (1779–4934ppm) and Hf (23–64 ppm). This HFSE enrichment is linkedto growth of large (up to 2 mm) zircon and niobian rutile crystals(up to 3 modal %) near the time of eclogite metamorphism. Thediamond-bearing eclogites on the other hand are characterizedby high-MgO (19·6–21·3 wt %) garnet andultralow-Na₂O (0·44–1·50 wt %) clinopyroxene.Paleotemperature estimates indicate that both the zircon- anddiamond-bearing eclogites have similar equilibration temperaturesof 950–1020°C and 990–1030°C, respectively,corresponding to mantle depths of 150–180 km. Integrationof petrographic, whole-rock and mineral geochemistry, geochronologyand isotope tracer techniques indicates that the Jericho zircon-bearingeclogite xenoliths have had a complex history involving Paleoproterozoicmetamorphism, thermal perturbations, and two or more episodesof Precambrian mantle metasomatism. The oldest metasomatic event(Type 1) occurred near the time of Paleoproterozoic metamorphism(1·8 Ga) and is responsible for the extreme HFSE enrichmentand growth of zircon and high-niobian rutile. A second thermalperturbation and concomitant carbonatite metasomatism (Type2) is responsible for significant apatite growth in some xenolithsand profound light rare earth element enrichment. Type 2 metasomatismoccurred in the period 1·0–1·3 Ga and isrecorded by relatively consistent whole-rock eclogite modelNd ages and secondary U–Pb zircon upper intercept dates.These eclogite xenoliths were derived from a variety of protoliths,some of which could represent metasomatized pieces of oceaniccrust, possibly linked to east-dipping subduction beneath theSlave craton during construction of the 1·88–1·84Ga Great Bear continental arc. Others, including the diamond-bearingeclogites, could be cumulates from mafic or ultramafic sillcomplexes that intruded the Slave lithospheric mantle at depthsof about 150–180 km. KEY WORDS: zircon- and diamond-bearing eclogites; Jericho kimberlite, geochronology; Precambrian metasomatism, northern Slave Craton     

Granulite sulphides as tracers of lower crustal origin and evolution: An example from the Slave craton, Canada

We carried out a detailed study of sulphide minerals, a ubiquitous mineral group in lower crustal mafic to peraluminous granulite xenoliths from the Diavik kimberlites, to assess their use in constraining the origin and tectonothermal evolution of the deep crust, and to obtain additional data on the composition of lower crust beneath ancient continents. Sulphides are overwhelmingly pyrrhotite with minor Ni (0.7-3.9 at.%), Co (0.1-0.7 at.%), and Cu contents (0.4-3.9 at.%). Sulphide modes in mafic granulites range from 0.14 to 0.55 vol%, translating into bulk rock S contents from ∼600 to 2000 ppm, similar to S contents in other mafic igneous rocks and indicating preservation of primary igneous S contents. In mafic granulites, Re and Os abundances in sulphides range from 42.5 to 726 ppb and 3.2 to 180 ppb, respectively, whereas those in peraluminous granulites are distinctly lower (36.1-282 ppb and 1.8-7.2 ppb, respectively), suggestive of Re and Os loss to fractionating sulphides in the more evolved precursors of these rocks.The significant within-sample variability of ¹⁸⁷Os/¹⁸⁸Os and correlation with ¹⁸⁷Re/¹⁸⁸Os indicates the preservation of primary Re-Os isotope systematics and time-integrated decay of the measured ¹⁸⁷Re. Within the large uncertainties inherent in the nature of the samples and technique, sulphides in some granulites may record major tectonothermal events in the central Slave craton spanning several billion years of evolution. Multiple generations of sulphide can occur in a single sample. These data attest to the heterogeneous composition and complex history of the Slave craton lower crust.     

Spatial distribution of kimberlite in the Slave craton, Canada: a geometrical approach

Exploration within the Slave craton has revealed clusters of kimberlite intrusions, commonly with internally consistent geochemical and temporal characteristics. Translation diagrams (“Fry analysis”) allow an unbiased geometrical examination of the distance and direction between each kimberlite occurrence and all others in the database. Recurrent patterns are visually accentuated due to the square function in data density. Circular histograms quantify the azimuthal density of kimberlite at various distances. For this study, the database comprises the geographic position of 212 kimberlite occurrences of which 70% are from the Lac de Gras field (LDG). Analyses are presented separately for the LDG data and for all non-LDG data in order to test for regional variations and to avoid overwhelming the craton-scale studies by the high density of LDG data.

Empirical grouping of kimberlite locations results in delineation of five elliptical clusters that encompass all but four kimberlite occurrences. Clusters within the western part of the craton are elongate to the north–northeast and align within a narrow zone (“Western Corridor”). Elsewhere, the clusters are elongate to the northwest or west–northwest and appear to be arranged en echelon within a poorly defined north–northwest trending zone (“Central Corridor”). Geometrical spatial analyses of kimberlite locations highlight the craton-scale pattern of emplacement within the two main corridors. At regional and local scales, individual intrusions are preferentially located towards the west–northwest (ca. 280°) and north–northeast (ca. 015°) of other intrusions, and these orientations are interpreted to reflect upper mantle trends in magma generation. At local scales (10–25 km), kimberlite of the central and southern craton tends to be located to the northeast (ca. 045°), and possibly weakly to the east–northeast (ca. 070°), of other intrusions, and these orientations correspond to major crustal fractures systems. It is proposed that kimberlite emplacement is controlled primarily by the interaction of elongate 280° and 015° source regions with near-surface deviations influenced by crustal fracture systems.

The 015° trend evident at craton, regional, and local scales is parallel to a swarm of alkaline diabase dykes that are concentrated in a ca. 30-km-wide corridor passing through Lac de Gras. A profound spatial association between significantly diamondiferous kimberlite and the margins of the dyke corridor suggests the corridor is the surface expression of a mantle-depth structure. It remains unclear whether the proposed mantle structure coincides with a diamond-rich zone near the base of the lithosphere, or delineates pathways favorable for diamond preservation during emplacement. The linear array of kimberlite within the western craton forms a parallel corridor that may be an analogous mantle structure, but which to date has failed to yield economic diamond concentrations.     

A review of Palaeoarchaean felsic volcanism in the eastern Kaapvaal craton:Linking plutonic and volcanic records

In the Kaapvaal craton of southern Africa, as well as other Archaean cratons worldwide, the progression from dominant tonalite-trondhjemite-granodiorite(TTG) to granite-monzogranite-syenogranite(GMS)rock types is interpreted to reflect progressive reworking and differentiation of the continental crust.Here we re-evaluate the early Archaean evolution of the Kaapvaal craton and propose a unified view of the plutonic and volcanic records based on elemental and isotopic(Nd, Hf) data and zircon U-Pb ages.We also report new whole-rock major and trace element analyses, zircon U-Pb ages and Hf-in-zircon analyses of igneous clasts from a conglomerate of the 3.2 Ga Moodies Group of the Barberton Greenstone Belt. Many of these clasts are derived from shallow intrusive rocks of granitic composition, which are scarcely represented in outcrop. Despite alteration, the volcanic rocks can be classified based on their trace element contents into two main groups by comparison with plutonic rocks. One group has characteristics resembling TTGs: relatively low and fractionated rare earth element concentrations with no Eu anomaly and relatively low concentrations of high field strength elements(Nb mostly ≤12 ppm). The second group has GMS-like characteristics: less fractionated REE, marked negative Eu anomalies and HFSE-increasing trends with progressing fractionation(Nb ≤ 50 ppm or more, Th up to 30-40 ppm). In addition, igneous clasts of Moodies Group conglomerate have chemical, mineralogical and isotopic characteristics that link them to GMS. New analyses of some of these clasts indicate elevated high field strength elements(Nb up to 20 ppm) and_(εHf)(t)of zircon down to -3.5. These rocks imply the presence of an already differentiated felsic crust at 3.5 Ga, which has Nd and Hf model ages indicating mantle extraction ages extending back to the Eoarchaean. The combined record of plutonic and volcanic rocks of the Kaapvaal craton provides a more complex scenario than previously suggested and indicates that TTG and GMS-like felsic magmas were emplaced broadly coevally in multiple pulses between ~3.5 Ga and 3.2 Ga.     

Dynamic Earth: crustal and mantle heterogeneity

The dynamic processes within the Earth leave their record in geophysical and geochemical variation about the general stratification with depth. A snapshot of current structure is provided by geophysical evidence, whereas geochemical information provides a perspective over the age of the Earth. The combination of information on the distribution of heterogeneity from geophysical and geochemical sources provides enhanced insight into likely geodynamic processes. A variety of techniques can be used to examine crustal structure, but the major source of information on seismic heterogeneity within the Earth comes from tomographic studies, exploiting surface waves for the lithosphere and body waves for the bulk Earth. A powerful tool for examining the character of mantle heterogeneity is the comparison of images of bulk-sound and shear-wave speed extracted in a single inversion, since this isolates the dependencies on the elastic moduli. Such studies are particularly effective when a common path coverage is achieved for P and S as, for example, when common source and receiver pairs are extracted for arrival times of the phases. The relative behaviour of the bulk-sound and shear-wave speeds allows the definition of heterogeneity regimes. For subduction zones, a large part of the imaged structure comes from S-wave speed variations. The narrow segments of fast wave speeds in the lower mantle, in the depth range 900 – 1500 km, are dominated by S variations, with very little bulk-sound contribution, so images of P-wave speed are controlled by shear-wave variability. Deep in the mantle, there are many features with high seismic-wave speed without an obvious association with subduction in the last 100 million years, which suggests long-lived preservation of components of the geodynamic cycle. The base of the Earth's mantle is a complex zone with widespread indications of heterogeneity on many scales, discontinuities of variable character, and shear-wave anisotropy. Discordance between P- and S-wave speed anomalies suggests the presence of chemical heterogeneity rather than just the effect of temperature.     

Peridotitic diamonds from the Slave and the Kaapvaal cratons—similarities and differences based on a preliminary data set

A comparison of the diamond productions from Panda (Ekati Mine) and Snap Lake with those from southern Africa shows significant differences: diamonds from the Slave typically are un-resorbed octahedrals or macles, often with opaque coats, and yellow colours are very rare. Diamonds from the Kaapvaal are dominated by resorbed, dodecahedral shapes, coats are absent and yellow colours are common. The first two features suggest exposure to oxidizing fluids/melts during mantle storage and/or transport to the Earth's surface, for the Kaapvaal diamond population.

Comparing peridotitic inclusions in diamonds from the central and southern Slave (Panda, DO27 and Snap Lake kimberlites) and the Kaapvaal indicates that the diamondiferous mantle lithosphere beneath the Slave is chemically less depleted. Most notable are the almost complete absence of garnet inclusions derived from low-Ca harzburgites and a generally lower Mg-number of Slave inclusions.

Geothermobarometric calculations suggest that Slave diamonds originally formed at very similar thermal conditions as observed beneath the Kaapvaal (geothermal gradients corresponding to 40–42 mW/m² surface heat flow), but the diamond source regions subsequently cooled by about 100–150 °C to fall on a 37–38 mW/m² (surface heat flow) conductive geotherm, as is evidenced from touching (re-equilibrated) inclusions in diamonds, and from xenocrysts and xenoliths. In the Kaapvaal, a similar thermal evolution has previously been recognized for diamonds from the De Beers Pool kimberlites. In part very low aggregation levels of nitrogen impurities in Slave diamonds imply that cooling occurred soon after diamond formation. This may relate elevated temperatures during diamond formation to short-lived magmatic perturbations.

Generally high Cr-contents of pyrope garnets (inside and outside of diamonds) indicate that the mantle lithosphere beneath the Slave originally formed as a residue of melt extraction at relatively low pressures (within the stability field of spinelperidotites), possibly during the extraction of oceanic crust. After emplacement of this depleted, oceanic mantle lithosphere into the Slave lithosphere during a subduction event, secondary metasomatic enrichment occurred leading to strong re-enrichment of the deeper (>140 km) lithosphere. Because of the extent of this event and the occurrence of lower mantle diamonds, this may be related to an upwelling plume, but it may equally just reflect a long term evolution with lower mantle diamonds being transported upwards in the course of “normal” mantle convection.     

A redox profile of the Slave mantle and oxygen fugacity control in the cratonic mantle

The authors report a redox profile based on Mössbauer data of spinel and garnet to a depth of 210 km from mantle xenoliths of the northern (N) and southeastern (SE) Slave craton (northern Canada). The profile transects three depth facies of peridotites that form segments of different bulk composition, represented by spinel peridotite, spinel–garnet peridotite, low-temperature garnet peridotite, high-temperature garnet peridotite, and pyroxenite. The shallow, more depleted N Slave spinel peridotite records lower oxygen fugacities compared to the deeper, less depleted N Slave spinel–garnet peridotite, consistent with their different spinel Fe³⁺ concentrations. Garnet peridotites show a general reduction in log fO₂ (FMQ)s with depth, where values for garnet peridotites are lower than those for spinel–garnet peridotites. There is a strong correlation between depletion and oxygen fugacity in the spinel peridotite facies, but little correlation in the garnet peridotite facies. The strong decrease in log fO₂ (FMQ) with depth that arises from the smaller partial molar volume of Fe³⁺ in garnet, and the observation of distinct slopes of log fO₂ (FMQ) with depth for spinel peridotite compared to spinel–garnet peridotite strongly suggest that oxygen fugacity in the cratonic peridotitic mantle is intrinsically controlled by iron equilibria involving garnet and spinel.

---

Sr evolution in the West Greenland-Labrador craton: a model for early Rb depletion in the mantle

A model for Sr evolution in the West Greenland-Labrador region is proposed based on the available data for granitic gneisses ranging in age from 3.6 AE to 1.9 AE and for recent continental tholeiites from West Greenland. The evolution of initial Sr isotopic composition is consistent with a two-stage model for the mantle in this area with a constant Rb/Sr ratio of 0.014 ± 0.002 from ～ 4.45 AE through the present. The model suggests the mantle may have been depleted in Rb ～ 4.45 AE ago if the Earth was originally chondritic with respect to Rb and Sr. This age of differentiation is consistent with recently proposed terrestrial differentiation based on Pb Pb analyses of ancient Amîtsoq gneiss feldspars (～4.47 AE, Gancarz and Wasserburg, 1977). A linear regression of the gneisses, if extrapolated to the present, predicts the present value of Sr⁸⁷/Sr⁸⁶ of the mantle in this area to be 0.7034 ± 0.0003 which is the measured Sr isotopic composition of the recent Svartenhuk tholeiites. Various implications of the model with respect to early Earth history, mantle evolution and pre-emplacement histories of gneissic suites are discussed.     

The Siberian lithosphere traverse: mantle terranes and the assembly of the Siberian Craton

The kimberlite fields scattered across the NE part of the Siberian Craton have been used to map the subcontinental lithospheric mantle (SCLM), as it existed during Devonian to Late Jurassic time, along a 1000-km traverse NE–SW across the Archean Magan and Anabar provinces and into the Proterozoic Olenek Province. 4100 garnets and 260 chromites from 65 kimberlites have been analysed by electron probe (major elements) and proton microprobe (trace elements). These data, and radiometric ages on the kimberlites, have been used to estimate the position of the local (paleo)geotherm and the thickness of the lithosphere, and to map the detailed distribution of specific rock types and mantle processes in space and time. A low geotherm, corresponding approximately to the 35 mW/m² conductive model of Pollack and Chapman [Tectonophysics 38, 279–296, 1977], characterised the Devonian lithosphere beneath the Magan and Anabar crustal provinces. The Devonian geotherm beneath the northern part of the area was higher, rising to near a 40 mW/m² conductive model. Areas intruded by Mesozoic kimberlites are generally characterised by this higher, but still ‘cratonic' geotherm. Lithosphere thickness at the time of kimberlite intrusion varied from ca. 190 to ca. 240 km beneath the Archean Magan and Anabar provinces, but was less (150–180 km) beneath the Proterozoic Olenek Province already in Devonian time. Thinner Devonian lithosphere (140 km) in parts of this area may be related to Riphean rifting. Near the northern end of the traverse, differences in geotherm, lithosphere thickness and composition between the Devonian Toluopka area and the nearby Mesozoic kimberlite fields suggest thinning of the lithosphere by ca. 50–60 km, related to Devonian rifting and Triassic magmatism. A major conclusion of this study is that the crustal terrane boundaries defined by geological mapping and geophysical data (extended from outcrops in the Anabar Shield) represent major lithospheric sutures, which continue through the upper mantle and juxtapose lithospheric domains that differ significantly in composition and rock-type distribution between 100 and 250 km depth. The presence of significant proportions of harzburgitic and depleted lherzolitic garnets beneath the Magan and Anabar provinces is concordant with their Archean surface geology. The lack of harzburgitic garnets, and the chemistry of the lherzolitic garnets, beneath most of the other fields are consistent with the Proterozoic surface rocks. Mantle sections for different terranes within the Archean portion of the craton show pronounced differences in bulk composition, rock-type distribution, metasomatic overprint and lithospheric thickness. These observations suggest that individual crustal terranes, of both Archean and Proterozoic age, had developed their own lithospheric roots, and that these differences were preserved during the Proterozoic assembly of the craton. Data from kimberlite fields near the main Archean–Proterozoic suture (the Billyakh Shear Zone) suggest that reworking and mixing of Archean and Proterozoic mantle was limited to a zone less than 100 km wide.     

Kimberlite AT-56: a mantle sample from the north central Superior craton, Canada

Kimberlite AT-56, discovered in February 2001, represents the most recent addition to the Attawapiskat kimberlite cluster, located in the James Bay Lowlands of Ontario, Canada. AT-56 is a small kimberlite body with a surface diameter of approximately 40 m and a steep southeastern plunge. It consists of a medium to coarse-grained matrix supported kimberlite with abundant olivine, clinopyroxene, garnet, ilmenite and mica macrocrysts in a green-black to orange-black matrix. The kimberlite is classified as a hypabyssal facies sparsely macrocrystic calcite kimberlite. Heavy mineral concentrates from two representative samples of AT-56 have been analyzed to characterize the mantle sampled by the kimberlite. Both samples yielded large heavy mineral concentrates comprised of roughly equal proportions of Mg-ilmenite, Cr-diopside, high-Cr garnet and low-Cr garnet. Mg-chromite is also present in quantities an order of magnitude less than the other constituents.

The high-Cr peridotitic garnet macrocrysts are only slightly more abundant than the low-Cr varieties, the population being dominated by G9 (lherzolitic) types with only a few (less than 10%) weakly sub-calcic G10 (probable harzburgitic) garnets present. Ni thermometry results for a representative selection of G9 and G10 garnets indicate that the majority equilibrated at temperatures ranging from 1000 to 1250 °C. A significant proportion of the low-Cr garnet population derived from AT-56 is characterized by relatively low-Ti (0.2 to 0.4 wt.% TiO₂) and elevated Na (0.07 to 0.13 wt.% Na₂O) contents characteristic of Group 1, diamond inclusion type eclogite garnets. These sodic garnets have elevated Cr₂O₃ contents (typically 1 to 2 wt.% Cr₂O₃), suggesting they may be websteritic in origin rather than eclogitic. Comparison of AT-56 garnet compositions with published data available for other Attawapiskat kimberlites suggests websteritic mantle has also been sampled by kimberlite bodies elsewhere in the Attawapiskat cluster and it may be an important diamond reservoir in this area.