Heterogeneous modification and reactivation of a craton margin beneath the Korean Peninsula from teleseismic travel time tomography

Margins of old continental lithosphere are likely prone to ongoing modification processes. Therefore, constraining detailed structures beneath the margin can be essential in understanding the evolution of the continental lithosphere. The eastern margin of the Eurasian plate is a natural laboratory that allows us to study the strong effects from multiple episodes of continental collision and subduction of different oceanic plates since their formation. To reveal the detailed evolution of cratons at their margins, we describe, for the first time, the upper mantle structures beneath the southern Korean Peninsula (SKP) based strictly on teleseismic relative arrival time data from densely deployed local seismic arrays, which allows us to constrain the details of the lithospheric structures beneath the Archean-Proterozoic basement. We imaged a thick (~150 km) high-velocity anomaly mainly beneath the Proterozoic Yeongnam Massif with large velocity contrasts (dlnVp ≈ 4.0% and dlnVs ≈ 6.0%) at its boundaries, suggesting the presence of a long-lasting cratonic root in the southwestern SKP. On the other hand, low-velocity anomalies were found beneath the Proterozoic Gyeonggi Massif, Gyeongsang arc-back-arc basin, and along the eastern margin of the SKP, indicating significantly modified regions. The possible existence of a remnant cratonic root beneath the SKP and contrasting lithospheric structures across the different Precambrian massifs suggests the highly heterogeneous modification of cratonic lithosphere at the eastern Eurasian plate margin. Strong velocity reductions, which indicate a thermally elevated upper mantle potentially with partial melts, correspond to areas of Cenozoic basalts, high surface heat flow, and high topography along the eastern KP margin. We interpret this coincidence as a result of recent reactivation of a craton margin, which is controlled by intense interaction between the convective upper mantle and heterogeneous continental lithosphere.     

Global 1° × 1° thermal model TC1 for the continental lithosphere: Implications for lithosphere secular evolution

This paper reports a new 1° × 1° global thermal model for the continental lithosphere (TC1). Geotherms for continental terranes of different ages (> 3.6 Ga to present) constrained by reliable data on borehole heat flow measurements (Artemieva, I.M., Mooney, W.D. 2001. Thermal structure and evolution of Precambrian lithosphere: a global study. J. Geophys. Res 106, 16387–16414.), are statistically analyzed as a function of age and are used to estimate lithospheric temperatures in continental regions with no or low-quality heat flow data (ca. 60% of the continents). These data are supplemented by cratonic geotherms based on electromagnetic and xenolith data; the latter indicate the existence of Archean cratons with two characteristic thicknesses, ca. 200 and > 250 km. A map of tectono-thermal ages of lithospheric terranes complied for the continents on a 1° × 1° grid and combined with the statistical age relationship of continental geotherms (z = 0.04  t + 93.6, where z is lithospheric thermal thickness in km and t is age in Ma) formed the basis for a new global thermal model of the continental lithosphere (TC1). The TC1 model is presented by a set of maps, which show significant thermal heterogeneity within continental upper mantle, with the strongest lateral temperature variations (as large as 800 °C) in the shallow mantle. A map of the depth to a 550 °C isotherm (Curie isotherm for magnetite) in continental upper mantle is presented as a proxy to the thickness of the magnetic crust; the same map provides a rough estimate of elastic thickness of old (> 200 Ma) continental lithosphere, in which flexural rigidity is dominated by olivine rheology of the mantle.Statistical analysis of continental geotherms reveals that thick (> 250 km) lithosphere is restricted solely to young Archean terranes (3.0–2.6 Ga), while in old Archean cratons (3.6–3.0 Ga) lithospheric roots do not extend deeper than 200–220 km. It is proposed that the former were formed by tectonic stacking and underplating during paleocollision of continental nuclei; it is likely that such exceptionally thick lithospheric roots have a limited lateral extent and are restricted to paleoterrane boundaries. This conclusion is supported by an analysis of the growth rate of the lithosphere since the Archean, which does not reveal a peak in lithospheric volume at 2.7–2.6 Ga as expected from growth curves for juvenile crust.A pronounced peak in the rate of lithospheric growth (10–18 km³/year) at 2.1–1.7 Ga (as compared to 5–8 km³/year in the Archean) well correlates with a peak in the growth of juvenile crust and with a consequent global extraction of massif-type anorthosites. It is proposed that large-scale variations in lithospheric thickness at cratonic margins and at paleoterrane boundaries controlled anorogenic magmatism. In particular, mid-Proterozoic anorogenic magmatism at the cratonic margins was caused by edge-driven convection triggered by a fast growth of the lithospheric mantle at 2.1–1.7 Ga. Belts of anorogenic magmatism within cratonic interiors can be caused by a deflection of mantle heat by a locally thickened lithosphere at paleosutures and, thus, can be surface manifestations of exceptionally thick lithospheric roots. The present volume of continental lithosphere as estimated from the new global map of lithospheric thermal thickness is 27.8 (± 7.0) × 10⁹ km³ (excluding submerged terranes with continental crust); preserved continental crust comprises ca. 7.7 × 10⁹ km³. About 50% of the present continental lithosphere existed by 1.8 Ga.     

Integrated models of diamond formation and craton evolution

Two decades of diamond research in southern Africa allow the age, average N content and carbon composition of diamonds, and the dominant paragenesis of their syngenetic silicate and sulfide inclusions to be integrated on a cratonwide scale with a model of craton formation. Individual eclogitic sulfide inclusions in diamonds from the Kimberley area kimberlites, Koffiefontein, Orapa and Jwaneng have Re–Os isotopic ages that range from circa 2.9 Ga to the mid-Proterozoic and display little correspondence with the prominent variations in the P-wave velocity (±1%) that the mantle lithosphere shows at depths within the diamond stability field (150–225 km). Silicate inclusions in diamonds and their host diamond compositions for the above kimberlites, Finsch, Jagersfontein, Roberts Victor, Premier, Venetia, and Letlhakane show a regional relationship to the seismic velocity of the lithosphere. Mantle lithosphere with slower P-wave velocity relative to the craton average correlates with a greater proportion of eclogitic vs. peridotitic silicate inclusions in diamond, a greater incidence of younger Sm–Nd ages of silicate inclusions, a greater proportion of diamonds with lighter C isotopic composition, and a lower percentage of low-N diamonds. The oldest formation ages of diamonds support a model whereby mantle that became part of the continental keel of cratonic nuclei first was created by middle Archean (3.2–3.3 Ga or older) mantle depletion events with high degrees of melting and early harzburgite formation. The predominance of eclogitic sulfide inclusions in the 2.9 Ga age population links late Archean (2.9 Ga) subduction–accretion events to craton stabilization. These events resulted in a widely distributed, late Archean generation of eclogitic diamonds in an amalgamated craton. Subsequent Proterozoic tectonic and magmatic events altered the composition of the continental lithosphere and added new lherzolitic and eclogitic diamonds to the already extensive Archean diamond suite. Similar age/paragenesis systematics are seen for the more limited data sets from the Slave and Siberian cratons.     

Lithospheric growth at margins of cratons

Deep seismic reflection profiles collected across Proterozoic–Archean margins are now sufficiently numerous to formulate a consistent hypothesis of how continental nuclei grow laterally to form cratonic shields. This picture is made possible both because the length of these regional profiles spans all the tectonic elements of an orogen on a particular cratonic margin and because of their great depth range. Key transects studied include the LITHOPROBE SNORCLE 1 transect and the BABEL survey, crossing the Slave and Baltic craton margins, respectively. In most cases, the older (Archean) block appears to form a wedge of uppermost mantle rock embedded into the more juvenile (Proterozoic) block by as much as 100–200 km at uppermost mantle depths and Archean lithosphere is therefore more laterally extensive at depth than at the surface. Particularly bright reflections along the Moho are cited as evidence of shear strain within a weak, low-viscosity lower crustal channel that lies along the irregular top of the indenting wedge. The bottom of the wedge is an underthrust/subduction zone, and associated late reversal in subduction polarity beneath the craton margin emerges as a common characteristic of these margins although related arc magmatism may be minor.     

Crustal structure and properties of Archean cratons of Gondwanaland: similarity and difference

This is a synopsis of available data the on crustal structure and properties of thirteen Archean cratons of Gondwanaland (the cratons of Africa, Australia, Antarctica, South America, and the Indian subcontinent). The data include estimates of surface area, rock age and lithology, Moho depth, thickness of lithosphere and sediments, as well as elevations, all summarized in a table. The cratons differ in size from 0.05 x 10⁶ km² (Napier craton) to 4 x 10⁶ km² (Congo craton) and span almost the entire Archean period from 3.8 to 2.5 Ga. Sediments are mostly thin, though reach 7 km in the Congo and West African cratons. Elevations above sea level are from 0 to 2 km; some relatively highland cratons (Kaapvaal, Zimbabwe, and Tanzanian) rise to more than 1 km. On the basis of regional seismic data, the Moho map for cratons has been improved. The Moho diagrams for each craton are constructed. The analysis of the available new data shows that the average Moho depth varies from 33 to 44 km: Pilbara (33 km), Grunehogna (35 km), Sao Francisco (36 km), Yilgarn (37 km), Dharwar (38 km), Tanzanian (39 km), Zimbabwe (39 km), Kaapvaal (40 km), Gawler (40 km), Napier (40 km), West Africa (40 km), Congo (42 km), and Amazon (44 km) cratons. The Moho depth within the cratons is less uniform than it was assumed before: from 28 to 52 km. The new results differ significantly from the earlier inference of a relatively flat Moho geometry beneath Archean cratons. According to the new data, early and middle Archean undeformed crust is characterized by a shallow Moho depth (28-38 km), while late Archean or deformed crust may be as thick as 52 km.     

Constraints on the thermal evolution of continental lithosphere from U-Pb accessory mineral thermochronometry of lower crustal xenoliths,southern Africa

U-Pb isotopic thermochronometry of rutile, apatite and titanite from kimberlite-borne lower crustal granulite xenoliths has been used to constrain the thermal evolution of Archean cratonic and Proterozoic off-craton continental lithosphere beneath southern Africa. The relatively low closure temperature of the U-Pb rutile thermochronometer (~400-450 °C) allows its use as a particularly sensitive recorder of the establishment of "cratonic" lithospheric geotherms, as well as subsequent thermal perturbations to the lithosphere. Contrasting lower crustal thermal histories are revealed between intracratonic and craton margin regions. Discordant Proterozoic (1.8 to 1.0 Ga) rutile ages in Archean (2.9 to 2.7 Ga) granulites from within the craton are indicative of isotopic resetting by marginal orogenic thermal perturbations influencing the deep crust of the cratonic nucleus. In Proterozoic (1.1 to 1.0 Ga) granulite xenoliths from the craton-bounding orogenic belts, rutiles define discordia arrays with Neoproterozoic (0.8 to 0.6 Ga) upper intercepts and lower intercepts equivalent to Mesozoic exhumation upon kimberlite entrainment. In combination with coexisting titanite and apatite dates, these results are interpreted as a record of postorogenic cooling at an integrated rate of approximately 1 °C/Ma, and subsequent variable Pb loss in the apatite and rutile systems during a Mesozoic thermal perturbation to the deep lithosphere. Closure of the rutile thermochronometer signals temperatures of 𙠂 °C in the lower crust during attainment of cratonic lithospheric conductive geotherms, and such closure in the examined portions of the "off-craton" Proterozoic domains of southern Africa indicates that their lithospheric thermal profiles were essentially cratonic from the Neoproterozoic through to the Late Jurassic. These results suggest similar lithospheric thickness and potential for diamond stability beneath both Proterozoic and Archean domains of southern Africa. Subsequent partial resetting of U-Pb rutile and apatite systematics in the cratonic margin lower crust records a transient Mesozoic thermal modification of the lithosphere, and modeling of the diffusive Pb loss from lower crustal rutile constrains the temperature and duration of Mesozoic heating to 𙡦 °C for ₞ ka. This result indicates that the thermal perturbation is not simply a kimberlite-related magmatic phenomenon, but is rather a more protracted manifestation of lithospheric heating, likely related to mantle upwelling and rifting of Gondwana during the Late Jurassic to Cretaceous. The manifestation of this thermal pulse in the lower crust is spatially and temporally correlated with anomalously elevated and/or kinked Cretaceous mantle paleogeotherms, and evidence for metasomatic modification in cratonic mantle peridotite suites. It is argued that most of the geographic differences in lithospheric thermal structure inferred from mantle xenolith thermobarometry are likewise due to the heterogeneous propagation of this broad upper mantle thermal anomaly. The differential manifestation of heating between cratonic margin and cratonic interior indicates the importance of advective heat transport along pre-existing lithosphere-scale discontinuities. Within this model, kimberlite magmatism was a similarly complex, space- and time-dependent response to Late Mesozoic lithospheric thermal perturbation.     

Mesozoic thermal evolution of the southern African mantle lithosphere

The thermal structure of Archean and Proterozoic lithospheric terranes in southern Africa during the Mesozoic was evaluated by thermobarometry of mantle peridotite xenoliths erupted in alkaline magmas between 180 and 60 Ma. For cratonic xenoliths, the presence of a 150–200 °C isobaric temperature range at 5–6 GPa confirms original interpretations of a conductive geotherm, which is perturbed at depth, and therefore does not record steady state lithospheric mantle structure.

Xenoliths from both Archean and Proterozoic terranes record conductive limb temperatures characteristic of a “cratonic” geotherm (40 mW m⁻²), indicating cooling of Proterozoic mantle following the last major tectonothermal event in the region at 1 Ga and the probability of thick off-craton lithosphere capable of hosting diamond. This inference is supported by U–Pb thermochronology of lower crustal xenoliths [Schmitz and Bowring, 2003. Contrib. Mineral. Petrol. 144, 592–618].

The entire region then suffered a protracted regional heating event in the Mesozoic, affecting both mantle and lower crust. In the mantle, the event is recorded at 150 Ma to the southeast of the craton, propagating to the west by 108–74 Ma, the craton interior by 85–90 Ma and the far southwest and northwest by 65–70 Ma. The heating penetrated to shallower levels in the off-craton areas than on the craton, and is more apparent on the southern margin of the craton than in its western interior. The focus and spatial progression mimic inferred patterns of plume activity and supercontinent breakup 30–100 Ma earlier and are probably connected.

Contrasting thermal profiles from Archean and Proterozoic mantle result from penetration to shallower levels of the Proterozoic lithosphere by heat transporting magmas. Extent of penetration is related not to original lithospheric thickness, but to its more fertile character and the presence of structurally weak zones of old tectonism. The present day distribution of surface heat flow in southern Africa is related to this dynamic event and is not a direct reflection of the pre-existing lithospheric architecture.     

Granitoid events in space and time: Constraints from igneous and detrital zircon age spectra

The goal of this study is to evaluate the global age distribution of granitoid magmatism and juvenile continental crust production with U/Pb isotopic ages from igneous and detrital zircons, and with Nd isotopic data. Granitoid age peaks, which are largely defined by TIMS data, are narrow and precise in contrast to detrital peaks that are often broad and hump-shaped due to the larger uncertainties of SHRIMP and LAM-ICPMS data. Granitic age peaks do not always have detrital counterparts and vice versa. Possible contributing factors to this mismatch are removal of crustal sources by erosion, inadequate sampling of granitoids because of cover by younger rocks, or small age peaks hidden by large age peaks in detrital spectra.Seven igneous peaks are found on five or more cratons or continents (3300, 2700, 2680, 2500, 2100, 1900 and 1100 Ma) and seven detrital peaks occur on three or more continents (2785, 2700, 2600, 2500, 1900, 1650 and 1200 Ma). Nd isotope distributions suggest important additions of juvenile continental crust at 2700, 2500, 2120, 1900, 1700, 1650, 800, 570 and 450 Ma. Tight clusters of craton ages occur for Superior–Karelia, Sao Francisco–Nain, and Kaapvaal–Siberia in the early Archean and for Wyoming–Kaapvaal–Slave, Superior–Nain, and West Africa–Amazonia in the late Archean. The global 2700-Ma peak is not a simple spike, but involves several peaks between 2760 and 2650 Ma. Events older than 3700 Ma are limited to the Yilgarn, Slave, Nain and North China cratons, and events between 2600 and 2500 Ma are widespread only in East Asia, Central and East Africa, and India.Single, short-lived mantle plume events at 2700 and 1900 Ga (or any other time) cannot easily account for prolonged episodes of granitoid magmatism during the Precambrian. The causes of geographically widespread and geographically restricted events are probably not the same.     

A top to bottom lithospheric study of Africa and Arabia

We study the lithospheric structure of Africa, Arabia and adjacent oceanic regions with fundamental-mode surface waves over a broad period range. Including group velocities with periods shorter than 35 s allows us to examine shallower features than previous studies of the whole continent. In the process, we have developed a crustal thickness map of Africa. Main features include crustal thickness increases under the West African, Congo, and Kalahari cratons. We find crustal thinning under Mesozoic and Cenozoic rifts, including the Benue Trough, Red Sea, and East, Central, and West African rift systems, along with less abrupt crustal thickness changes at passive continental margins. We also find crustal thickness differences in North Africa between the West African Craton and East Saharan Shield. Crustal shear wave velocities are generally faster in oceanic regions and cratons, and slower in more recent crust and in active and remnant orogenic regions. Deeper structure, related to the thickness of cratons and modern rifting, is generally consistent with previous work. Under cratons we find thick lithosphere and fast upper mantle velocities, while under rifts we find thinned lithosphere and slower upper mantle velocities. However, we also find the lack of a thick cratonic keel beneath the central portion of the Congo Craton. There are no consistent effects in areas classified as hotspots, indicating that there seem to be numerous origins for these features. Finally, it appears that the African Superswell, which is responsible for high elevation and uplift over large portions of Africa, has had a significantly different impact (as indicated by features such as temperature, time of influence, etc.) in the north and the south. This is consistent with episodic activity at shallow depths, which is well-expressed in northeastern Africa and Arabia today.     

Heat flow and destabilized cratons: a comparative study of the North China,Siberian, and Wyoming cratons

Heat has a close relation to cratonic stability. We reconstructed the thermal history of the eastern North China, Wyoming, and Siberian cratons by using xenolith and thermal indicator data, and then the thermal lithospheric thickness in geological time was calculated by combining the crustal structural and thermal parameters. Researches from three cratons demonstrate that heat-flow values of the eastern North China Craton (eNCC) and Wyoming Craton (WC) have increased exceeding 20 mW/m² in the Meso-Cenozoic and the thermal lithospheric thickness has decreased by >100 km because of the cratonic destabilization. However, the thermal history variation of the Siberian Craton (SC) is small, ~5 mW/m², the lithosphere has remained stable and it has experienced relatively little lithospheric thinning <100 km. The differences of the cratonic thermal regime are related to the breakup of the Pangea supercontinent. Destabilized cratons on active continental margin are directly affected by plate subduction. Stable cratons in the interior of the supercontinents are directly affected by the upwelling of superplumes and that may result in lithospheric thinning. Differences in the thermal evolutions of cratons provide new geothermal evidence for the different thermal regimes associated with the breakup of the Pangea supercontinent.     

Evolution of the Australian lithosphere

The evolution of the Australian plate can be interpreted in a plate‐tectonic paradigm in which lithospheric growth occurred via vertical and horizontal accretion. The lithospheric roots of Archaean lithosphere developed contemporaneously with the overlying crust. Vertical accretion of the Archaean lithosphere is probably related to the arrival of large plumes, although horizontal lithospheric accretion was also important to crustal growth. The Proterozoic was an era of major crustal growth in which the components of the North Australian, West Australian and South Australian cratons were formed and amalgamated during a series of accretionary events and continent‐continent collisions, interspersed with periods of lithospheric extension. During Phanerozoic accretionary tectonism, approximately 30% of the Australian crust was added to the eastern margin of the continent in a predominantly supra‐subduction environment. Widespread plume‐driven rifting during the breakup of Gondwana may have contributed to the destruction of Archaean lithospheric roots (as a result of lithospheric stretching). However, lithospheric growth occurred at the same time due to mafic underplating along the eastern margin of the plate. Northward drift of Australia during the Tertiary led to the development of a complex accretionary margin at the leading edge of the plate (Papua New Guinea).     

Structure and evolution of an obliquely sheared continental margin: Rio Muni, West Africa

Oblique-shear margins are divergent continental terrains whose breakup and early drift evolution are characterized by significant obliquity in the plate divergence vector relative to the strike of the margin. We focus on the Rio Muni margin, equatorial West Africa, where the ca. 70-km-wide Ascension Fracture Zone (AFZ) exhibits oblique–slip faulting and synrift half-graben formation that accommodated oblique extension during the period leading up to and immediately following whole lithosphere failure and continental breakup (ca. 117 Ma). Oblique extension is recorded also by strike–slip and oblique–slip fault geometry within the AFZ, and buckling of Aptian synrift rocks in response to block rotation and local transpression. Rio Muni shares basic characteristics of both rifted and transform margins, the end members of a spectrum of continental margin kinematics. At transform margins, continental breakup and the onset of oceanic spreading (drifting) are separate episodes recorded by discrete breakup and drift unconformities. Oceanic opening will proceed immediately following breakup on a rifted margin, whereas transform and oblique-shear margins may experience several tens of millennia between breakup and drift. Noncoeval breakup and drift have important consequences for the fit of the equatorial South American and African margins because, in reconstructing the configuration of conjugate continental margins at the time of their breakup, it cannot be assumed that highly segmented margins like the South Atlantic will match each other at their ocean–continent boundaries (OCBs). Well known ‘misfits’ in reconstructions of South Atlantic continental margins may be accounted for by differential timing of breakup and drifting between oblique-shear margins and their adjacent rifted segments.     

1891–1883 Ma Southern Bastar–Cuddapah mafic igneous events, India: A newly recognized large igneous province

A newly recognized remnant of a Paleoproterozoic Large Igneous Province has been identified in the southern Bastar craton and nearby Cuddapah basin from the adjacent Dharwar craton, India. High precision U–Pb dates of 1891.1 ± 0.9 Ma (baddeleyite) and 1883.0 ± 1.4 Ma (baddeleyite and zircon) for two SE-trending mafic dykes from the BD2 dyke swarm, southern Bastar craton, and 1885.4 ± 3.1 Ma (baddeleyite) for a mafic sill from the Cuddapah basin, indicate the existence of 1891–1883 Ma mafic magmatism that spans an area of at least 90,000 km² in the south Indian shield.This record of 1.9 Ga mafic/ultramafic magmatism associated with concomitant intracontinental rifting and basin development preserved along much of the south-eastern margin of the south Indian shield is a widespread geologic phenomenon on Earth. Similar periods of intraplate mafic/ultramafic magmatism occur along the margin of the Superior craton in North America (1.88 Ga Molson large igneous province) and in southern Africa along the northern margin of the Kaapvaal craton (1.88–1.87 Ga dolerite sills intruding the Waterberg Group). Existing paleomagnetic data for the Molson and Waterberg 1.88 Ga large igneous provinces indicate that the Superior and Kalahari cratons were at similar paleolatitudes at 1.88 Ga but a paleocontinental reconstruction at this time involving these cratons is impeded by the lack of a robust geological pin such as a Limpopo-like 2.0 Ga deformation zone in the Superior Province. The widespread occurrence of 1.88 Ga intraplate and plate margin mafic magmatism and basin development in numerous Archean cratons worldwide likely reflects a period of global-scale mantle upwelling or enhanced mantle plume activity at this time.     

Heat flow and the structure of Precambrian lithosphere

Studies of heat flow from Precambrian terrains have demonstrated three empirical relationships; a temporal relationship between heat flow and tectonic age, a spatial pattern between heat flow and the proximity of Archean cratons, and a temporal relationship between heat flow and the age of lithosphere stabilization. In the first relationship, heat flow is inversely related to tectonic age. The second pattern is characterized by low heat flow from Archean cratons and Proterozoic terrains adjacent to cratonic margins (pericratonic terrains), and higher heat flow from Proterozoic terrains that are more than a few hundred kilometers from a craton. In the third pattern, heat flow decreases as the age of stabilization of the lithosphere increases. A number of interpretations of Precambrian heat flow have been offered to explain one or more of these relationships. The simple cooling of a thermal boundary layer predicts essentially no change in heat flow in terrains older than 1.5 Ga, and therefore does not likely provide a comprehensive framework for the interpretation of Precambrian heat flow. By contrast, two other interpretations, (1) thicker lithosphere beneath Archean terrains than beneath Proterozoic terrains, and (2) greater heat production in Proterozoic crust than in Archean crust, when combined with the special structural configuration of sutures, can both contribute to the spatial and temporal heat flow distributions. Xenolith thermobarometry constraints on lithospheric temperatures, however, limit the contribution of age-dependent crustal heat production, and therefore at least part of the heat flow distributions derive from variations in lithosphere thickness.     

Speculations on evolution of the terrestrial lithosphere–asthenosphere system—Plumes and plates

In the early Earth, accretionary impact heating, including collision with a large, Mars-sized object, decay of short-lived radioisotopes, and (after an initial thermal run-up) continuous segregation of the liquid Fe–Ni core resulted in extensive the melting of the silicate mantle and in the formation of a near-surface magma mush ocean. Progressive, continuous degassing and chemical–gravitational differentiation of the crust–mantle system accompanied this Hadean stage, and has gradually lessened during the subsequent cooling of the planet. Mantle and core overturn was vigorous in the Hadean Earth, reflecting deep-seated chemical heterogeneities and concentrations of primordial heat. Hot, bottom-up mantle convection, including voluminous plume ascent, efficiently rid the planet of much thermal energy, but gradually decreased in importance with the passage of time. Formation of lithospheric scum began when planetary surface temperatures fell below those of basalt and peridotite solidi. Thickening and broadening of lithospheric plates are inferred from the post-Hadean rock record. Developmental stages of mantle circulation included: (a) 4.5–4.4 Ga, early, chaotic magma ocean circulation involving an incipient or pre-plate regime; (b) 4.4–2.7 Ga, growth of small micro-oceanic and microcontinental platelets, all returned to the mantle prior to 4.0 Ga, but increasing in size and progressively suturing sialic crust-capped lithospheric amalgams at and near the surface over time; (c) 2.7–1.0 Ga, assembly of cratons surmounting larger, supercontinental plates; and (d) 1.0 Ga–present, modern, laminar-flowing asthenospheric cells capped by gigantic, Wilson-cycle lithospheric plates. Restriction of komatiitic lavas to the Archean, and of ophiolite complexes ± alkaline igneous rocks, high-pressure and ultrahigh-pressure metamorphic terranes to progressively younger Proterozoic–Phanerozoic orogenic belts supports the idea that planetary thermal relaxation promoted the increasingly negative buoyancy of cooler oceanic lithosphere. The Thickening of oceanic plates enhanced the gravitational instability and the consequent overturn of the outer Earth as cold, top-down oceanic mantle convection. The scales and dynamics of deep-seated asthenospheric circulation, and of lithospheric foundering + shallow asthenospheric return flow evidently have evolved gradually over geologic time in response to the progressive cooling of the Earth.     

Non‐uniform splitting of a single mantle plume by double cratonic roots: Insight into the origin of the central and southern East African Rift System

Using numerical thermo‐mechanical experiments we analyse the role of an active mantle plume and pre‐existing lithospheric thickness differences in the structural development of the central and southern East African Rift system. The plume‐lithosphere interaction model setup captures the essential features of the studied area: two cratonic bodies embedded into surrounding lithosphere of normal thickness. The results of the numerical experiments suggest that localization of rift branches in the crust is mainly defined by the initial position of the mantle plume relative to the cratons. We demonstrate that development of the Eastern branch, the Western branch and the Malawi rift can be the result of non‐uniform splitting of the Kenyan plume, which has been rising underneath the southern part of the Tanzanian craton. Major features associated with Cenozoic rifting can thus be reproduced in a relatively simple model of the interaction between a single mantle plume and pre‐stressed continental lithosphere with double cratonic roots.     

Permo–Triassic intraplate magmatism and rifting in Eurasia: implications for mantle plumes and mantle dynamics

At the transition from the Permian to the Triassic, Eurasia was the site of voluminous flood-basalt extrusion and rifting. Major flood-basalt provinces occur in the Tunguska, Taymyr, Kuznetsk, Verkhoyansk–Vilyuy and Pechora areas, as well as in the South Chinese Emeishen area. Contemporaneous rift systems developed in the West Siberian, South Kara Sea and Pyasina–Khatanga areas, on the Scythian platform and in the West European and Arctic–North Atlantic domain. At the Permo–Triassic transition, major extensional stresses affected apparently Eurasia, and possibly also Pangea, as evidenced by the development of new rift systems. Contemporaneous flood-basalt activity, inducing a global environmental crisis, is interpreted as related to the impingement of major mantle plumes on the base of the Eurasian lithosphere. Moreover, the Permo–Triassic transition coincided with a period of regional uplift and erosion and a low-stand in sea level. Permo–Triassic rifting and mantle plume activity occurred together with a major reorganization of plate boundaries and plate kinematics that marked the transition from the assembly of Pangea to its break-up. This plate reorganization was possibly associated with a reorganization of the global mantle convection system. On the base of the geological record, we recognize short-lived and long-lived plumes with a duration of magmatic activity of some 10–20 million years and 100–150 million years, respectively. The Permo–Triassic Siberian and Emeishan flood-basalt provinces are good examples of “short-lived” plumes, which contrast with such “long lived” plumes as those of Iceland and Hawaii. The global record indicates that mantle plume activity occurred episodically. Purely empirical considerations indicate that times of major mantle plume activity are associated with periods of global mantle convection reorganization during which thermally driven mantle convection is not fully able to facilitate the necessary heat transfer from the core of the Earth to its surface. In this respect, we distinguish between two geodynamically different scenarios for major plume activity. The major Permo–Triassic plume event followed the assembly Pangea and the detachment of deep-seated subduction slabs from the lithosphere. The Early–Middle Cretaceous major plume event, as well as the terminal–Cretaceous–Paleocene plume event, followed a sharp acceleration of global sea-floor spreading rates and the insertion of new subduction zone slabs deep into the mantle. We conclude that global plate kinematics, driven by mantle convection, have a bearing on the development of major mantle plumes and, to a degree, also on the pattern of related flood-basalt magmatism.     

Lithosphere mapping beneath the North American plate

Major- and trace-element analyses of garnets from heavy-mineral concentrates have been used to derive the compositional and thermal structure of the subcontinental lithospheric mantle (SCLM) beneath 16 areas within the core of the ancient Laurentian continent and 11 areas in the craton margin and fringing mobile belts. Results are presented as stratigraphic sections showing variations in the relative proportions of different rock types and metasomatic styles, and the mean Fo content of olivine, with depth. Detailed comparisons with data from mantle xenoliths demonstrate the reliability of the sections.

In the Slave Province, the SCLM in most areas shows a two-layer structure with a boundary at 140–160 km depth. The upper layer shows pronounced lateral variations, whereas the lower layer, after accounting for different degrees of melt-related metasomatism, shows marked uniformity. The lower layer is interpreted as a subcreted plume head, added at ca. 3.2 Ga; this boundary between the layers rises to <100 km depth toward the northern and southern edges of the craton. Strongly layered SCLM suggests that plume subcretion may also have played a role in the construction of the lithosphere beneath Michigan and Saskatchewan.

Outside the Slave Province, most North American Archon SCLM sections are less depleted than similar sections in southern Africa and Siberia; this may reflect extensive metasomatic modification. In E. Canada, the degree of modification increases toward the craton margin, and the SCLM beneath the Kapuskasing Structural Zone is typical of that beneath Proterozoic to Phanerozoic mobile belts.

SCLM sections from several Proterozoic areas around the margin of the Laurentian continental core (W. Greenland, Colorado–Wyoming district, Arkansas) show discontinuities and gaps that are interpreted as the effects of lithosphere stacking during collisional orogeny. Some areas affected by Proterozoic orogenesis (Wyoming Craton, Alberta, W. Greenland) appear to retain buoyant, modified Archean SCLM. Possible juvenile Proterozoic SCLM beneath the Colorado Plateau is significantly less refractory. The SCLM beneath the Kansas kimberlite field is highly melt-metasomatised, reflecting its proximity to the Mid-Continent Rift System.

A traverse across the continent shows that the upper part of the cratonic SCLM is highly magnesian; the decrease in mg# with depth is interpreted as the cumulative effect of metasomatic modification through time. The relatively small variations in seismic velocity within the continental core largely reflect the thickness of this depleted layer. The larger drop in seismic velocity in the surrounding Proton and Tecton belts reflects the closely coupled changes in SCLM composition and geotherm.     

Sm–Nd isotopic mapping of lithospheric growth and stabilization in the eastern Kaapvaal craton

Whole-rock Sm–Nd isotope systematics of 79 Archean granitoids from the eastern Kaapvaal craton, southern Africa, are used to delineate lithospheric boundaries and to constrain the timescale of crustal growth, assembly and geochemical differentiation c. 3.66–2.70 Ga. Offsets in ε_Nd values for 3.2–3.3 Ga granitoids across the Barberton greenstone belt (BGB) are consistent with existing models for c. 3.23 Ga accretion of newly formed lithosphere north of the BGB onto pre-existing c. 3.66 Ga lithosphere south of the BGB along a doubly verging subduction margin. The Nd isotopic signature of c. 3.3–3.2 Ga magmatic rocks show that significant crustal growth occurred during subduction–accretion. After c. 3.2 Ga, however, the Nd signature of intrusive rocks c. 3.1 and 2.7 Ga is dominated by intracrustal recycling rather than by new additions from the mantle, signalling cratonic stability.     

A Kaapvaal craton debate: Nucleus of an early small supercontinent or affected by an enhanced accretion event?

Incorporation of the Kaapvaal craton within a speculative Neoarchaean–Palaeoproterozoic supercontinent has long been debated, and this idea provides a potential solution to solving the apparently enigmatic provenance of the huge quantities of gold within the famous Witwatersrand auriferous deposits of Kaapvaal. Within a framework of a postulated Neoarchaean “Kenorland” (“northern”; present-day reference) supercontinent, we examine possible “southern” cratons that may have been contiguous with Kaapvaal: Pilbara, Zimbabwe, Dharwar, São Francisco, Amazon, Congo. Brief reviews of their basic geology and inferred evolution in syn-Witwatersrand basin times (c. 3.1–2.8 Ga) show no obvious support for any such supercontinental amalgamations. An alternative idea to explain a measure of gross similarity amongst several Neoarchaean cratons is through global events, such as a c. 3125–3000 Ma cratonic-scale erosive event interpreted for both Pilbara and Kaapvaal, and a much more widespread magmatic event at c. 2760–2680 Ma. We postulate that a global superplume event at c. 3.0 Ga included a plume beneath the Kaapvaal cratonic nucleus, thus halting any subduction around that terrane due to the thermal anomaly. Such a speculative global magmatic event is assumed to have enhanced production of juvenile oceanic crust at mid-ocean ridges, including those “offshore” of the thermally elevated Kaapvaal nucleus. Intra-oceanic obduction complexes may have built up fairly rapidly under such conditions, globally, and once the plume event had abated, “normal” plate tectonics would have resulted in composite (greenstone-tonalite, possibly also including granite) terranes accreting with nuclei such as Kaapvaal. This enhanced plume-related cratonic growth can be seen as a rapid accretion event. Formation of the envisaged ophiolite complexes possibly encompassed deformation-related first-order concentration of gold, and once accretion occurred around Kaapvaal's nucleus, from north and west (present-day frame of reference), a second-order (deformation-related) gold concentration may have resulted. The third order of gold concentration would logically have occurred once placer systems reworked detritus derived from the orogens along the N and W margins of Kaapvaal. Such conditions and placer gold deposits are known from many Neoarchaean cratons. The initial source of gold was presumably from the much hotter Mesoarchaean mantle and may have been related to major changes in Earth's tectonic regime at c. 3.0 Ga. The unique nature of Kaapvaal is probably its early stabilization, enabling formation of a complex flexural foreland basin system, in which vast quantities of placer sediments and heavy minerals could be deposited, and preserved from younger denudation through a unique post-Witwatersrand history.