Archean geotherms and supracrustal assemblages

Metamorphic mineral assemblages suggest the existence of variable geotherms and lithospheric thicknesses beneath late Archean continental crust. Archean granite-greenstone terranes reflect steep geotherms (50–70°C/km) while high-grade terranes reflect moderate geotherms similar to present continental crust with high heat flow (25–40°C/km). Corresponding lithosphere thicknesses for each terrane during the late Archean are 35–50 km and 50–75 km, respectively.Early Archean ( 3.0 b.y.) greenstones differ from late Archean ( 2.7 b.y.) greenstones by the rarity or absence of andesite and graywacke and the relative abundance of pelite, quartzite, and komatiite. Mature clastic sediments in early greenstones reflect shallow-water, stable-basin deposition. Such rocks, together with granite-bearing conglomerate and felsic volcanics imply the existence of still older granitic source terranes. The absence or rarity of andesite in early greenstones reflects the absence of tectonic conditions in which basaltic and tonalitic magmas are modified to produce andesite.A model is presented in which early Archean greenstones form at the interface between tonalite islands and oceanic lithosphere, over convective downcurrents; high-grade supracrustals form on stable continental edges or interiors; and late Archean greenstones form in intracontinental rifts over mantle plumes.     

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.     

Earth''s earliest continental lithosphere, hydrothermal flux and crustal recycling

The Kaapvaal craton in southern Africa and the Pilbara craton of northwestern Australia are the largest regions on Earth to have retained relatively pristine mid-Archaean rocks (3.0–4.0 Ga).The Kaapvaal craton covers about 1.2×10⁶ km², and varies in lithospheric thickness between 170 and 350 km. At surface, the craton can be subdivided into a number of Archaean sub-domains; some of the subdomains are also well defined at depth, and local variations in tomography of the lithosphere correspond closely with subdomain boundaries at surface.The Archaean history of the Kaapvaal craton spans about 1 Gyr and can be conveniently subdivided into two periods, each of about the same length as the Phanerozoic. The first period, from circa 3.7-3.1 Ga, records the initial separation of the cratonic lithosphere from the asthenosphere, terminating with a major pulse of accretion tectonics between 3.2 and 3.1 Ga, which includes the formation of “paired metamorphic belts”. This period of continental growth can be compared to plate tectonic processes occurring in modern-day oceanic basins. However, the difference is that in the mid-Archaean, these oceanic processes appear to have occurred in shallower water depths than the modern ocean basins. The second period, from circa 3.1-2.6 Ga, records intra-continental and continental-edge processes: continental growth during this period occurred predominantly through a combination of tectonic accretion of crustal fragments and subduction-related igneous processes, in much the same way as has been documented along the margins of the Pacific and Tethys oceans since the Mesozoic.The intra-oceanic processes resulted in small, but deep-rooted continental nucleii; the first separation of this early continental lithosphere could only have occurred when the mean elevation of mid-oceanicridges sank below sea-level. Substantial recycling of continental lithosphere into the mantle must have occurred during this period of Earth history. During the second period, at least two large continental nucleii amalgamated during collisional processes which, together with internal chemical differentiation processes, created the first stable continental landmass. This landmass, which is known to have been substantially bigger than its present outline, may have been part of the Earth's first supercontinent.The oldest known subdomains of the craton include the oceanic-like rocks of the Barberton greenstone belt. The comagmatic mafic-ultramafic rocks (3.48–3.49 Ga) of this belt represent a remnant of very early oceanic-like lithosphere (known as the Jamestown Ophiolite Complex), which was obducted, approximately 45 Ma after its formation, onto a volcanic arc-like terrain by processes similar to those which have emplaced modern ophiolites at convergent margins of Phanerozoic continents. The early metamorphic history, metamorphic mineralogy, oxygen isotope profiles and degree of hydration of the 3.49 Ga Jamestown Ophiolite Complex are similar to present day subseafloor hydrothermal systems. The ratio of ΔMg to ΔSi for hydrothermally altered igneous rocks, both present day and Archaean, are remarkably uniform at −5(±0.9) and the same as that of hydrothermal fluids venting on the present-day East Pacific Rise. This observation suggests that the process of Mg exchange for Si in hydrothermal systems was commonplace throughout Earth's history.The chemistry of vent fluids and hydrothermally altered igneous rocks was combined with an inventory of ³He in the mantle to model Earth's total hydrothermal flux. An Archaean flux (at 3.5 Ga) of about 10 times present day was accompanied by a correspondingly greater abundance of Mg(OH), SiO₂, carbonate and Fe---Mn metasomatic rock types as well as massive sulphides. Assuming a constant column of seawater since the Archaean, the average residence time of seawater in the oceanic crust was 1.65−8.90×10⁵ years in the Archaean. Assuming that ³He and heat are transported from the mantle in silicate melts in uniform proportions, the model stipulates that accretion of oceanic crust decreased from about 3.43−6.5×10¹⁷ g/yr to a present-day rate of 0.52−0.8×10¹⁷ g/yr, with a drop in heat flow from 1.4−2.6×10²⁰ cal/yr to 2.1−3.2×10¹⁹ cal/year.The total amounts of SiO₂ and Fe mobilised in marine hydrothermal systems since 3.5 Ga is less than their masses in the present exosphere reservoirs (crust, hydrosphere, atmosphere). The total amounts of Mg, K, CO₂, Ca and Mn are greater than their respective masses in exosphere reservoirs; therefore, they must have been recycled into mantle. The total mass of recycled hydrothermal components is small compared to the mass of the mantle. The flux of volatiles in hydrothermal systems is large compared to their volume in the atmosphere suggesting that the CO₂ and O₂ budgets of the atmosphere have been influenced by hydrothermal processes, especially in the Archaean.     

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.     

Plume head–lithosphere interactions near intra-continental plate boundaries

Plume–lithosphere interactions (PLI) have important consequences both for tectonic and mineralogical evolution of the lithosphere: for example, Archean metallogenic crises at the boundaries of the West African and Australian cratons coincide with postulated plume events. In continents, PLI are often located near boundaries between younger plates (e.g., orogenic) and older stable plates (e.g., cratons), which represent important geometrical, thermal and rheological barriers that interact with the emplacement of the plume head (e.g., Archean West Africa, East Africa, Pannonian–Carpathian system). The observable PLI signatures are conditioned by plume dynamics but also by lithosphere rheology and structure. We address the latter problem by considering a free-surface numerical model of PLI with two stratified elasto-viscous–plastic (EVP) lithospheric plates, one of which is older and thicker than another. The results show that: (1) plume head flattening is asymmetric, it is blocked from one side by the cold vertical boundary of the older plate, which leads to the mechanical decoupling of the crust from the mantle lithosphere, and to localized faulting at the cratonic margin; (2) the return flow from the plume head results in sub-vertical down-thrusting (delamination) of the lithosphere at the margin, producing sharp vertical cold boundary down to the 400 km depth; (3) plume head flattening and migration towards the younger plate results in concurrent surface extension above the centre of the plume and in compression (pushing), down-thrusting and magmatic events at the cratonic margin (down-thrusting is also produced at the opposite border of the younger plate); these processes may result in continental growth at the “craton side”; (4) topographic signatures of PLI show basin-scale uplifts and subsidences preferentially located at cratonic margins. Negative Rayleigh–Taylor instabilities in the lithosphere above the plume head provide a mechanism for crustal delamination. Inferred consequences of PLI near intra-continental plate boundaries, such as faulting at cratonic edges and enhanced magmatic activity, could explain plume-related metallogenic crises, as suggested for West Africa and Australia.     

Formation and evolution of Precambrian continental lithosphere in South China

An overview is presented for the formation and evolution of Precambrian continental lithosphere in South China. This is primarily based on an integrated study of zircon U–Pb ages and Lu–Hf isotopes in crustal rocks, with additional constraints from Re–Os isotopes in mantle-derived rocks. Available Re–Os isotope data on xenolith peridotites suggest that the oldest subcontinental lithospheric mantle beneath South China is primarily of Paleoproterozoic age. The zircon U–Pb ages and Lu–Hf isotope studies reveal growth and reworking of the juvenile crust at different ages. Both the Yangtze and Cathaysia terranes contain crustal materials of Archean U–Pb ages. Nevertheless, zircon U–Pb ages exhibit two peaks at 2.9–3.0 Ga and ~ 2.5 Ga in Yangtze but only one peak at ~ 2.5 Ga in Cathaysia. Both massive rocks and crustal remnants (i.e., zircon) of Archean U–Pb ages occur in Yangtze, but only crustal remnants of Archean U–Pb ages occur in Cathaysia. Zircon U–Pb and Lu–Hf isotopes in the Kongling complex of Yangtze suggest the earliest episode of crustal growth in the Paleoarchean and two episodes of crustal reworking at 3.1–3.3 Ga and 2.8–3.0 Ga. Both negative and positive ε_Hf(t) values are associated with Archean U–Pb ages of zircon in South China, indicating both the growth of juvenile crust and the reworking of ancient crust in the Archean. Paleoproterozoic rocks in Yangtze exhibit four groups of U–Pb ages at 2.1 Ga, 1.9–2.0 Ga, ~ 1.85 Ga and ~ 1.7 Ga, respectively. They are associated not only with reworking of the ancient Archean crust in the interior of Yangtze, but also with the growth of the contemporaneous juvenile crust in the periphery of Yangtze. In contrast, Paleoproterozoic rocks in Cathaysia were primarily derived from reworking of Archean crust at 1.8–1.9 Ga. The exposure of Mesoproterozoic rocks are very limited in South China, but zircon Hf model ages suggest the growth of juvenile crust in this period due to island arc magmatism of the Grenvillian oceanic subduction. Magmatic rocks of middle Neoproterozoic U–Pb ages are widespread in South China, exhibiting two peaks at about 830–800 Ma and 780–740 Ma, respectively. Both negative and positive ε_Hf(t) values are associated with the middle Neoproterozoic U–Pb ages of zircon, suggesting not only growth and reworking of the juvenile Mesoproterozoic crust but also reworking of the ancient Archean and Paleoproterozoic crust in the middle Neoproterozoic. The tectonic setting for this period of magmatism would be transformed from arc–continent collision to continental rifting with reference to the plate tectonic regime in South China.     

Proterozoic–Paleozoic development of the basement of the Central Andes (18–26°S) — a mobile belt of the South American craton

The Late Precambrian–Early Paleozoic metamorphic basement forms a volumetrically important part of the Andean crust. We investigated its evolution in order to subdivide the area between 18 and 26°S into crustal domains by means of petrological and age data (Sm–Nd isochrons, K–Ar). The metamorphic crystallization ages and t_DM ages are not consistent with growth of the Pacific margin north of the Argentine Precordillera by accretion of exotic terranes, but favor a model of a mobile belt of the Pampean Cycle. Peak metamorphic conditions in all scattered outcrop areas between 18 and 26°S are similar and reached the upper amphibolite facies conditions indicated by mineral paragensis and the occurrence of migmatite. Sm–Nd mineral isochrons yielded 525±10, 505±6 and 509±1 Ma for the Chilean Coast Range, the Chilean Precordillera and the Argentine Puna, and 442±9 and 412±18 Ma for the Sierras Pampeanas. Conventional K–Ar cooling age data of amphibole and mica cluster around 400 Ma, but are frequently reset by Late Paleozoic and Jurassic magmatism. Final exhumation of the Early Paleozoic orogen is confirmed by Devonian erosional unconformities. Sm–Nd depleted mantle model ages of felsic rocks from the metamorphic basement range from 1.4 to 2.2 Ga, in northern Chile the average is 1.65±0.16 Ga (1σ; n=12), average t_DM of both gneiss and metabasite in NW Argentina is 1.76±0.4 Ga (1σ; n=22), and the isotopic composition excludes major addition of juvenile mantle derived material during the Early Paleozoic metamorphic and magmatic cycle. These new data indicate a largely similar development of the metamorphic basement south of the Arequipa Massif at 18°S and north of the Argentine Precordillera at 28°S. Variations of metamorphic grade and of ages of peak metamorphism are of local importance. The protolith was derived from Early to Middle Proterozoic cratonic areas, similar to the Proterozoic rocks from the Arequipa Massif, which had undergone Grenvillian metamorphism at ca. 1.0 Ga.     

Effects of buoyancy and mechanical layering on collisional deformation of continental lithosphere: Results from physical modeling

We use two suites of lithospheric-scale physical experiments to investigate the manner in which deformation of the continental lithosphere is affected by both (1) variations of lithospheric density (quantified by the net buoyant mass per area in the lithospheric mantle layer, M_B), and (2) the degree of coupling between the crust and lithospheric mantle (characterized by a modified Ampferer ratio, Am). The dynamics of the experiments can be characterized with a Rayleigh–Taylor type ratio, _CLM. Models with a positively buoyant lithospheric mantle layer (M_B > 0 and _CLM > 0) result in distributed root formation and a wide deformation belt. In contrast, models with a negatively buoyant lithospheric mantle layer strongly coupled to the crust (M_B < 0, 0 > _CLM > ≈ − 0.2, and Am > ≈ 10^{− 3}) exhibit localized roots and narrow deformation belts. Syncollisional delamination of the model lithospheric mantle layer and a wide deformation belt is exhibited in models with negatively buoyant lithospheric mantle layers weakly coupled to the crust (M_B < 0, _CLM < 0, and Am < ≈ 10^{− 3}). Syncollisional delamination of the continental lithosphere may initiate due to buoyancy contrasts within the continental plate, instead of resulting from wedging by the opposing plate. Rayleigh–Taylor instabilities dominate the style of deformation in models with a negatively buoyant lithospheric mantle layer strongly coupled to the crust and a slow convergence rate (M_B < 0 and _CLM > ≈ − 0.2). The degree of coupling (Am) between the model crust and lithospheric mantle plays a lesser role in both the style of lower-lithospheric deformation and the width of the crustal deformed zone with increasing density of the lithospheric mantle layer.     

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.     

Continental rifting as a function of lithosphere mantle strength

The role of the uppermost mantle strength in the pattern of lithosphere rifting is investigated using a thermo-mechanical finite-element code. In the lithosphere, the mantle/crust strength ratio (S_M/S_C) that decreases with increasing Moho temperature T_M allows two strength regimes to be defined: mantle dominated (S_M > S_C) and crust dominated (S_M < S_C). The transition between the two regimes corresponds to the disappearance of a high strength uppermost mantle for T_M > 700 °C. 2D numerical simulations for different values of S_M/S_C show how the uppermost mantle strength controls the style of continental rifting. A high strength mantle leads to strain localisation at lithosphere scale, with two main patterns of narrow rifting: “coupled crust–mantle” at the lowest T_M values and “deep crustal décollement” for increasing T_M values, typical of some continental rifts and non-volcanic passive margins. The absence of a high strength mantle leads to distributed deformations and wide rifting in the upper crust. These numerical results are compared and discussed in relation with series of classical rift examples.     

Regional patterns in the paragenesis and age of inclusions in diamond, diamond composition, and the lithospheric seismic structure of Southern Africa

The Archean lithospheric mantle beneath the Kaapvaal–Zimbabwe craton of Southern Africa shows ±1% variations in seismic P-wave velocity at depths within the diamond stability field (150–250 km) that correlate regionally with differences in the composition of diamonds and their syngenetic inclusions. Seismically slower mantle trends from the mantle below Swaziland to that below southeastern Botswana, roughly following the surface outcrop pattern of the Bushveld-Molopo Farms Complex. Seismically slower mantle also is evident under the southwestern side of the Zimbabwe craton below crust metamorphosed around 2 Ga. 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 Proterozoic and show little correspondence with these lithospheric variations. However, silicate inclusions in diamonds and their host diamond compositions for the above kimberlites, Finsch, Jagersfontein, Roberts Victor, Premier, Venetia, and Letlhakane do show some regional relationship to the seismic velocity of the lithosphere. Mantle lithosphere with slower P-wave velocity correlates with a greater proportion of eclogitic versus 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 whereas the converse is true for diamonds from higher velocity mantle. The oldest formation ages of diamonds indicate that the mantle keels which became continental nuclei were created by middle Archean (3.2–3.3 Ga) mantle depletion events with high degrees of melting and early harzburgite formation. The predominance of sulfide inclusions that are eclogitic in the 2.9 Ga age population links late Archean (2.9 Ga) subduction-accretion events involving an oceanic lithosphere component to craton stabilization. These events resulted in a widely distributed younger Archean generation of eclogitic diamonds in the lithospheric mantle. 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.     

Archean greenstone-tonalite duality: Thermochemical mantle convection models or plate tectonics in the early Earth global dynamics?

Mantle convection and plate tectonics are one system, because oceanic plates are cold upper thermal boundary layers of the convection cells. As a corollary, Phanerozoic-style of plate tectonics or more likely a different version of it (i.e. a larger number of slowly moving plates, or similar number of faster plates) is expected to have operated in the hotter, vigorously convecting early Earth. Despite the recent advances in understanding the origin of Archean greenstone–granitoid terranes, the question regarding the operation of plate tectonics in the early Earth remains still controversial. Numerical model outputs for the Archean Earth range from predominantly shallow to flat subduction between 4.0 and 2.5 Ga and well-established steep subduction since 2.5 Ga [Abbott, D., Drury, R., Smith, W.H.F., 1994. Flat to steep transition in subduction style. Geology 22, 937–940], to no plate tectonics but rather foundering of 1000 km sectors of basaltic crust, then “resurfaced” by upper asthenospheric mantle basaltic melts that generate the observed duality of basalts and tonalities [van Thienen, P., van den Berg, A.P., Vlaar, N.J., 2004a. Production and recycling of oceanic crust in the early earth. Tectonophysics 386, 41–65; van Thienen, P., Van den Berg, A.P., Vlaar, N.J., 2004b. On the formation of continental silicic melts in thermochemical mantle convection models: implications for early Earth. Tectonophysics 394, 111–124]. These model outputs can be tested against the geological record. Greenstone belt volcanics are composites of komatiite–basalt plateau sequences erupted from deep mantle plumes and bimodal basalt–dacite sequences having the geochemical signatures of convergent margins; i.e. horizontally imbricated plateau and island arc crust. Greenstone belts from 3.8 to 2.5 Ga include volcanic types reported from Cenozoic convergent margins including: boninites; arc picrites; and the association of adakites–Mg andesites- and Nb-enriched basalts.Archean cratons were intruded by voluminous norites from the Neoarchean through Proterozoic; norites are accounted for by melting of subduction metasomatized Archean continental lithospheric mantle (CLM). Deep CLM defines Archean cratons; it extends to  350 km, includes the diamond facies, and xenoliths signify a composition of the buoyant, refractory, residue of plume melting, a natural consequence of imbricated plateau-arc crust. Voluminous tonalites of Archean greenstone–granitoid terranes show a secular trend of increasing Mg#, Cr, Ni consistent with slab melts hybridizing with thicker mantle wedge as subduction angle steepens. Strike-slip faults of 1000 km scale; diachronous accretion of distinct tectonostratigraphic terranes; and broad Cordilleran-type orogens featuring multiple sutures, and oceanward migration of arcs, in the Archean Superior and Yilgarn cratons, are in common with the Altaid and Phanerozoic Cordilleran orogens. There is increasing geological evidence of the supercontinent cycle operating back to  2.7 Ga: Kenorland or Ur  2.7–2.4 Ga; Columbia  1.6–1.4 Ga; Rodinia  1100–750 Ma; and Pangea  230 Ma. High-resolution seismic reflection profiling of Archean terranes reveals a prevalence of low angle structures, and evidence for paleo-subduction zones. Collectively, the geological–geochemical–seismic records endorse the operation of plate tectonics since the early Archean.     

Secular evolution of the lithosphere beneath the eastern North China Craton: evidence from Mesozoic basalts and high-Mg andesites

Geochemical and isotopic data from Mesozoic lavas from the Jianguo, Niutoushan, Wulahada, and Guancaishan volcanic fields on the northern margin of the North China Craton provide evidence for secular lithospheric evolution of the region. Jianguo lavas are alkaline basalts with LILE- and LREE-enrichment ((La/Yb)_N=12.2-13.2) and MORB-like Sr-Nd-Pb isotopic ratios ((⁸⁷Sr/⁸⁶Sr)_i<0.704; ε_Nd=3.9-4.8; (²⁰⁶Pb/²⁰⁴Pb)_i≈18). Niutoushan basalts are similar but show evidence of olivine fractionation. Wulahada lavas are high-Mg andesites (Mg#∼67) with EM1 Sr-Nd-Pb isotopic signatures. Geochemical data suggest that the basalts originated from MORB-type asthenosphere whereas the high-Mg andesites were derived an EM1 mantle source, i.e., a refractory lithospheric mantle modified by a previously subducted slab. The result, combined with the available data of the Mesozoic basalts from the southern portion of the NCC (Zhang et al., 2002), manifests a vast secular evolution of the lithospheric mantle beneath the eastern NCC from the Paleozoic refractory continental lithosphere to this Mesozoic modified lithosphere. Compared with the cratonic margin, the lithospheric mantle beneath the center of the craton was less extensively modified, implying the secular evolution was related to the subduction processes surrounding the NCC. Therefore, we suggest that the interaction of the slab-derived silicic melt with the old refractory lithospheric mantle converted the Paleozoic cratonic lithospheric mantle into the late Mesozoic fertile mantle, which was also different from the Cenozoic counterpart. A geodynamic model is proposed to illustrate such a secular lithosphere evolution.     

http://www.sciencedirect.com/science/article/pii/S1674987111001125

<正>Greenstone belts of the eastern Dharwar Craton,India are reinterpreted as composite tectonostratigraphic terranes of accreted plume-derived and convergent margin-derived magmatic sequences based on new high-precision elemental data.The former are dominated by a komatiile plus Mg-tholeiitic basalt volcanic association,with deep water siliciclastic and banded iron formation(BIF) sedimentary rocks.Plumes melted at90 km under thin rifted continental lithosphere to preserve inlraoceanic and continental margin aspects.Associated alkaline basalts record subduction-recycling of Mesoarchean oceanic crust,incubated in the asthenosphere.and erupted coevally with Mg basalts from a heterogeneous mantle plume.Together.komaliites-Mg basalts-alkaline basalts plot along the Phanerozoic mantle array in Th/Yb versus Nb/Yb coordinate space,representing zoned plumes,establishing that these reservoirs were present in the Neoarchean mantle. Convergent margin magmatic associations are dominated by tholeiitic to calc-alkaline basalts eompositionally similar to recent intraoceanic arcs.As well,boninitic flows sourced in extremely depleted mantle are present,and the association of arc basalts with Mg-andesites-Nb enriched basalts-adakites documented from Cenozoic arcs characterized by subduction of young(20 Ma),hot,oceanic lithosphere. Consequently.Cenozoic style "hot" subduction was operating in the Neoarchean.These diverse volcanic associations were assembled to give composite terranes in a subduction-accretion orogen at～2.1 Ga,coevally with a global accretionary orogen at ~2.7 Ga,and associated orogenic gold mineralization. Archean lithospheric mantle,distinctive in being thick,refractory,and buoyant,formed complementary to the accreted plume and convergent margin terranes.as migrating arcs captured thick plumeplateaus. and the refractory,low density.residue of plume melting coupled with accreted imbricated plume-arc crust.     

The building blocks of continental crust: Evidence for a major change in the tectonic setting of continental growth at the end of the Archean

Oceanic arcs are commonly cited as primary building blocks of continents, yet modern oceanic arcs are mostly subducted. Also, lithosphere buoyancy considerations show that oceanic arcs (even those with a felsic component) should readily subduct. With the exception of the Arabian–Nubian orogen, terranes in post-Archean accretionary orogens comprise < 10% of accreted oceanic arcs, whereas continental arcs compose 40–80% of these orogens. Nd and Hf isotopic data suggest that accretionary orogens include 40–65% juvenile crustal components, with most of these (> 50%) produced in continental arcs.Felsic igneous rocks in oceanic arcs are depleted in incompatible elements compared to average continental crust and to felsic igneous rocks from continental arcs. They have lower Th/Yb, Nb/Yb, Sr/Y and La/Yb ratios, reflecting shallow mantle sources in which garnet did not exist in the restite during melting. The bottom line of these geochemical differences is that post-Archean continental crust does not begin life in oceanic arcs. On the other hand, the remarkable similarity of incompatible element distributions in granitoids and felsic volcanics from continental arcs is consistent with continental crust being produced in continental arcs.During the Archean, however, oceanic arcs may have been thicker due to higher degrees of melting in the mantle, and oceanic lithosphere would be more buoyant. These arcs may have accreted to each other and to oceanic plateaus, a process that eventually led to the production of Archean continental crust. After the Archean, oceanic crust was thinner due to cooling of the mantle and less melt production at ocean ridges, hence, oceanic lithosphere is more subductable. Widespread propagation of plate tectonics in the late Archean may have led not only to rapid production of continental crust, but to a change in the primary site of production of continental crust, from accreted oceanic arcs and oceanic plateaus in the Archean to primarily continental arcs thereafter.     

Stabilization and reactivation of cratonic lithosphere from the lower crustal record in the western Canadian shield

New U–Pb geochronology for an extensive exposure of high-pressure granulites in the East Lake Athabasca region of the western Canadian shield is consistent with a history characterized by 2.55 Ga stabilization of cratonic lithosphere, 650 million years of lower crustal residence and cratonic stability, and 1.9 Ga reactivation of the craton during lithospheric attenuation and asthenospheric upwelling. High precision single-grain and fragment zircon data define distinctive discordia arrays between 2.55 and 1.9 Ga. U–Pb ion microprobe spot analyses yield a similar range of U–Pb dates with no obvious correlation between date and cathodoluminescence zonation. We attribute the complex U–Pb zircon systematics to growth of the primary populations during a 2.55 Ga high-pressure granulite facies event (~1.3 GPa, 850°C) recorded by the dominant mineral assemblage of the mafic granulite gneisses, with subsequent zircon recrystallization and minor secondary zircon growth during a second high-pressure granulite facies event (1.0 GPa, ~800°C) at 1.9 Ga. The occurrence of two discrete granulite facies metamorphic events in the lower crust, separated by an interval of 650 million years that included isobaric cooling for at least some of this time, suggests that the rocks resided at lower crustal depths until 1.9 Ga. We infer that this phase of lower crustal residence and little tectonic activity is coincident with an extended period of cratonic stability. Detailed structural and thermochronological datasets indicate that multistage unroofing of the lower crustal rocks occurred in the following 200 million years. Extended lower crustal residence would logically be the history inferred for lower crust in most cratonic regions, but the unusual aspect of the history in the East Lake Athabasca region is the subsequent lithospheric reactivation that initiated transport of the lower crust to the surface. We suggest that a weakened strength profile related to the 1.9 Ga heating left the lithosphere susceptible to far-field tectonic stresses from bounding orogens that drove the lower crustal exhumation. An ultimate return to cratonic stability is responsible for the preservation of this extensive lower crustal exposure since 1.7 Ga.     

Density structure of the cratonic mantle in Southern Africa: 2. Correlations with kimberlite distribution,seismic velocities,and Moho sharpness

We present a new regional model for the depth-averaged density structure of the cratonic lithospheric mantle in southern Africa constrained on a 30′ × 30′ grid and discuss it in relation to regional seismic models for the crust and upper mantle, geochemical data on kimberlite-hosted mantle xenoliths, and data on kimberlite ages and distribution. Our calculations of mantle density are based on free-board constraints, account for mantle contribution to surface topography of ca. 0.5–1.0 km, and have uncertainty ranging from ca. 0.01 g/cm³ for the Archean terrains to ca. 0.03 g/cm³ for the adjacent fold belts. We demonstrate that in southern Africa, the lithospheric mantle has a general trend in mantle density increase from Archean to younger lithospheric terranes. Density of the Kaapvaal mantle is typically cratonic, with a subtle difference between the eastern, more depleted, (3.31–3.33 g/cm³) and the western (3.32–3.34 g/cm³) blocks. The Witwatersrand basin and the Bushveld Intrusion Complex appear as distinct blocks with an increased mantle density (3.34–3.35 g/cm³) with values typical of Proterozoic rather than Archean mantle. We attribute a significantly increased mantle density in these tectonic units and beneath the Archean Limpopo belt (3.34–3.37 g/cm³) to melt-metasomatism with an addition of a basaltic component. The Proterozoic Kheis, Okwa, and Namaqua–Natal belts and the Western Cape Fold Belt with the late Proterozoic basement have an overall fertile mantle (ca. 3.37 g/cm³) with local (100–300 km across) low-density (down to 3.34 g/cm³) and high-density (up to 3.41 g/cm³) anomalies. High (3.40–3.42 g/cm³) mantle densities beneath the Eastern Cape Fold belt require the presence of a significant amount of eclogite in the mantle, such as associated with subducted oceanic slabs.We find a strong correlation between the calculated density of the lithospheric mantle, the crustal structure, the spatial pattern of kimberlites, and their emplacement ages. (1) Blocks with the lowest values of mantle density (ca. 3.30 g/cm³) are not sampled by kimberlites and may represent the “pristine” Archean mantle. (2) Young (< 90 Ma) Group I kimberlites sample mantle with higher density (3.35 ± 0.03 g/cm³) than the older Group II kimberlites (3.33 ± 0.01 g/cm³), but the results may be biased by incomplete information on kimberlite ages. (3) Diamondiferous kimberlites are characteristic of regions with a low-density cratonic mantle (3.32–3.35 g/cm³), while non-diamondiferous kimberlites sample mantle with a broad range of density values. (4) Kimberlite-rich regions have a strong seismic velocity contrast at the Moho, thin crust (35–40 km) and low-density (3.32–3.33 g/cm³) mantle, while kimberlite-poor regions have a transitional Moho, thick crust (40–50 km), and denser mantle (3.34–3.36 g/cm³). We explain this pattern by a lithosphere-scale (presumably, pre-kimberlite) magmatic event in kimberlite-poor regions, which affected the Moho sharpness and the crustal thickness through magmatic underplating and modified the composition and rheology of the lithospheric mantle to make it unfavorable for consequent kimberlite eruptions. (5) Density anomalies in the lithospheric mantle show inverse correlation with seismic Vp, Vs velocities at 100–150 km depth. However, this correlation is weaker than reported in experimental studies and indicates that density-velocity relationship in the cratonic mantle is strongly non-unique.     

Density structure of the cratonic mantle in southern Africa: 1. Implications for dynamic topography

The origin of high topography in southern Africa is enigmatic. By comparing topography in different cratons, we demonstrate that in southern Africa both the Archean and Proterozoic blocks have surface elevation 500–700 m higher than in any other craton worldwide, except for the Tanzanian Craton. An unusually high topography may be caused by a low density (high depletion) of the cratonic lithospheric mantle and/or by the dynamic support of the mantle with origin below the depth of isostatic compensation (assumed here to be at the lithosphere base). We use free-board constraints to examine the relative contributions of the both factors to surface topography in the cratons of southern Africa. Our analysis takes advantage of the SASE seismic experiment which provided high resolution regional models of the crustal thickness.We calculate the model of density structure of the lithospheric mantle in southern Africa and show that it has an overall agreement with xenolith-based data for lithospheric terranes of different ages. Density of lithospheric mantle has significant short-wavelength variations in all tectonic blocks of southern Africa and has typical SPT values of ca. 3.37–3.41 g/cm³ in the Cape Fold and Namaqua–Natal fold belts, ca. 3.34–3.35 g/cm³ in the Proterozoic Okwa block and the Bushveld Intrusion Complex, ca. 3.34–3.37 g/cm³ in the Limpopo Belt, and ca. 3.32–3.33 g/cm³ in the Kaapvaal and southern Zimbabwe cratons.The results indicate that 0.5–1.0 km of surface topography, with the most likely value of ca. 0.5 km, cannot be explained by the lithosphere structure within the petrologically permitted range of mantle densities and requires the dynamic (or static) contribution from the sublithospheric mantle. Given a low amplitude of regional free air gravity anomalies (ca. + 20 mGal on average), we propose that mantle residual (dynamic) topography may be associated with the low-density region below the depth of isostatic compensation. A possible candidate is the low velocity layer between the lithospheric base and the mantle transition zone, where a temperature anomaly of 100–200 °C in a ca. 100–150 km thick layer may explain the observed reduction in Vs velocity and may produce ca. 0.5–1.0 km to the regional topographic uplift.     

Fission track and fault kinematics analyses for new insight into the Late Cenozoic tectonic regime changes in West-Central Sulawesi (Indonesia)

A 3-D density model for the Cretan and Libyan Seas and Crete was developed by gravity modelling constrained by five 2-D seismic lines. Velocity values of these cross-sections were used to obtain the initial densities using the Nafe–Drake and Birch empirical functions for the sediments, the crust and the upper mantle. The crust outside the Cretan Arc is 18 to 24 km thick, including 10 to 14 km thick sediments. The crust below central Crete at its thickest section, has values between 32 and 34 km, consisting of continental crust of the Aegean microplate, which is thickened by the subducted oceanic plate below the Cretan Arc. The oceanic lithosphere is decoupled from the continental along a NW–SE striking front between eastern Crete and the Island of Kythera south of Peloponnese. It plunges steeply below the southern Aegean Sea and is probably associated with the present volcanic activity of the southern Aegean Sea in agreement with published seismological observations of intermediate seismicity. Low density and velocity upper mantle below the Cretan Sea with ρ  3.25 × 10³ kg/m³ and V_p velocity of compressional waves around 7.7 km/s, which are also in agreement with observed high heat flow density values, point out at the mobilization of the upper mantle material here. Outside the Hellenic Arc the upper mantle density and velocity are ρ ≥ 3.32 × 10³ kg/m³ and V_p = 8.0 km/s, respectively. The crust below the Cretan Sea is thin continental of 15 to 20 km thickness, including 3 to 4 km of sediments. Thick accumulations of sediments, located to the SSW and SSE of Crete, are separated by a block of continental crust extended for more than 100 km south of Central Crete. These deep sedimentary basins are located on the oceanic crust backstopped by the continental crust of the Aegean microplate. The stretched continental margin of Africa, north of Cyrenaica, and the abruptly terminated continental Aegean microplate south of Crete are separated by oceanic lithosphere of only 60 to 80 km width at their closest proximity. To the east and west, the areas are floored by oceanic lithosphere, which rapidly widens towards the Herodotus Abyssal plain and the deep Ionian Basin of the central Mediterranean Sea. Crustal shortening between the continental margins of the Aegean microplate and Cyrenaica of North Africa influence the deformation of the sediments of the Mediterranean Ridge that has been divided in an internal and external zone. The continental margin of Cyrenaica extends for more than 80 km to the north of the African coast in form of a huge ramp, while that of the Aegean microplate is abruptly truncated by very steep fractures towards the Mediterranean Ridge. Changes in the deformation style of the sediments express differences of the tectonic processes that control them. That is, subduction to the northeast and crustal subsidence to the south of Crete. Strike-slip movement between Crete and Libya is required by seismological observations.     

Permian basaltic rocks in the Tarim basin, NW China: Implications for plume–lithosphere interaction

There are large areas of Permian basaltic rocks in the Tarim basin (PBRT) in northwestern China. Precise Ar–Ar dating of these rocks revealed an eruption age span of 262 to 285 Ma. Most of the PBRT is composed of alkaline basaltic rocks with high TiO₂ (2.43%–4.59%, weight percent), high Fe₂O₃ + FeO (12.63%–17.83%) and P₂O₅ (0.32%–1.38%) contents. Trace elements of these rocks have affinities with oceanic island basalts (OIB), as shown in chondrite normalized rare earth elements (REE) diagrams and primitive mantle normalized incompatible elements diagrams. The rocks show complex Sr–Nd isotopic character based on which they can be subdivided into two distinct groups: group 1 has relatively small initial (t = 280 Ma)⁸⁷Sr/⁸⁶Sr ratio ( 0.7048) and positive εNd(t) (3.42–4.66) values. Group 2 has relatively large initial ⁸⁷Sr/⁸⁶Sr ratio (0.7060–0.7083) and negative εNd(t) (from − 2.79 to − 2.16) values. Lead isotopes are even more complex with variations of (²⁰⁶Pb/²⁰⁴Pb)t, (²⁰⁷Pb/²⁰⁴Pb)t and (²⁰⁸Pb/²⁰⁴Pb)t ranging from 17.9265 to 18.5778, 15.4789 to 15.6067 and 37.2922 to 38.1437, respectively. Moreover, these two groups have different trace elements ratios such as Nb/La, Ba/Nb, Zr/Nb, Nb/Ta and Zr/Hf, implying different magmatic processes. Based on the geochemistry of basaltic rocks and an evaluation of the tectonics, deformation, and the compositions of crust and lithospheric mantle in Tarim, we conclude that these basaltic rocks resulted from plume–lithosphere interaction. Permian mantle plume caused an upwelling of the Tarim lithosphere leading to melting of the asthenospheric mantle by decompression. The magma ascended rapidly to the base of lower crust, where different degrees of assimilation of OIB-like materials and fractionation occurred. Group 1 rocks formed where the upwelling is most pronounced and the assimilation was negligible. In other places, different degrees of assimilation and fractionation account for the geochemical traits of group 2.