Abstract Waveform tomography with very large datasets reveals the upper-mantle structure of the Arctic in unprecedented detail. Using tomography jointly with computational petrology, we estimate temperature in the lithosphere–asthenosphere depth range and infer lithospheric structure and evolution. Most of the boundaries of the mantle roots of cratons in the Arctic are coincident with their geological boundaries at the surface. The thick lithospheres of the Greenland and North American cratons are separated by a corridor of thin lithosphere beneath Baffin Bay and through the middle of the Canadian Arctic Archipelago; the southern archipelago is part of the North American Craton. The mantle root of the cratonic block beneath northern Greenland may extend westwards as far as central Ellesmere Island. The Barents and Kara seas show high velocities indicative of thick lithosphere, similar to cratons. The locations of intraplate basaltic volcanism attributed to the High Arctic Large Igneous Province are all on thin, non-cratonic lithosphere. The lithosphere beneath the central part of the Siberian Traps is warmer than elsewhere beneath the Siberian Craton. This observation is consistent with lithospheric erosion associated with the large igneous province volcanism. A corridor of relatively low seismic velocities cuts east–west across central Greenland. This indicates lithospheric thinning, which appears to delineate the track of the Iceland hotspot. Supplementary material: Figures with comparisons of different tomographic models at 50 and 200 km depths are available at https://doi.org/10.6084/m9.figshare.c.3817810

Abstract: In this paper we illustrate the geochemical diversity of the basalts of Panjal ( c. 289 Ma), Rajmahal ( c. 117 Ma) and Deccan ( c. 65 Ma) continental large igneous provinces (continental LIPs) from India as controlled by multiple mantle sources. The Panjal primitive basalts are formed by 5–20%, Rajmahal basalts by 10–20% and Deccan picritic basalts by 1–20% melting of the mantle. The depth of melting beneath Deccan is greater than beneath Rajmahal, which greater than beneath Panjal. We have categorized the Panjal, Rajmahal and Deccan basalts into Type I and Type II. The Type I basalts with lower La/Nb, Th/Nb and Th/Yb ratios and higher ε Nd i values ( ε Nd i = −3 to +8) are derived from the sublithospheric sources; some of the Type I basalts represent lithosphere-contaminated melts. The Type II basalts with higher La/Nb, Th/Nb and Th/Yb ratios and lower ε Nd i values ( ε Nd i = −20 to −3) are produced from the subcontinental lithosphere. We estimate that lithosphere and sublithosphere (plume and/or E-MORB patches within asthenosphere) have contributed, respectively, approximately 67 and 33% to the Panjal, 52 and 48% to Rajmahal, and 28 and 72% to Deccan LIPs. Fertility of the Indian subcontinental lithosphere may have increased with decreasing age related to the break-up of Gondwana. Supplementary material: References for the geochemical data, extracted from the GEOROC database, on the basaltic rocks from Panjal, Rajmahal and Deccan Large Igneous Provinces are available at: https://doi.org/10.6084/m9.figshare.c.3821881

ABSTRACT In this paper, we present a comprehensive review of literature data (~2600 analyses), including major and trace elements and Sr-Nd isotopes, on continental flood basalts from the Ethiopia-Yemen, Deccan (India), and Karoo (southern Africa) volcanic provinces in order to evaluate whether they can be attributable to similar tectonomagmatic processes that occurred during the past 200 m.y. in central Gondwana. Results indicate that the three investigated provinces share fundamental features, such as the following: (1) Major and trace element compositions are closely comparable, in terms of parental magmas and fractionation trends, for the various continental flood basalt suites recognized in the provinces, namely, low Ti (LT, TiO 2 0.5–3 wt%), high Ti (HT1, TiO 2 1–4 wt%), and very high Ti (HT2, TiO 2 2.5–7 wt%). (2) There is a clear zonal arrangement of continental flood basalts, with the hottest (potential temperature T p up to ~1600 °C) and deepest (up to 5 GPa) HT picrite-basalt magmas in the central area and cooler and shallower LT basalts (T p down to 1450 °C, pressure [ P ] = 2–3 GPa) at the periphery, corresponding to a maximum thermal difference of 60–110 °C from the inner to the outer zones in each province. This conforms to continental flood basalt generation from a lenticular melting region, plausibly reflecting thermo-compositionally zoned plume heads, with maximum excess temperature T ex = 250–300 °C with respect to the notional mid-ocean-ridge basalt (MORB) ambient mantle. (3) The central area of all provinces is characterized by the nearly exclusive occurrence of superheated HT picrite-basalt (and nearly coeval alkaline-carbonatite complexes) at the intersection of multiple extensional lineaments (faulting, rifting, and dike swarms), reflecting the focus of the tectonomagmatic activities. (4) The common occurrence of rhyolitic differentiates at the top of picrite-basalt lavas (e.g., Lalibela suite, northern Ethiopia; Pavagadh suite, Deccan; Lebombo suite, African Karoo) has to be considered an effect of the inversion of the stress regime, from generalized regional extension (continental flood basalt eruption) to localized continental rifting accompanying magma differentiation to rhyolites; activity at some of these rift and dike systems, e.g., the Western Afar Escarpment, the coastal dikes of Western Deccan, and the Rooi Rand dikes, was protracted until continental breakup and the opening of new oceanic branches of the Red Sea, the central Indian Ocean, and southwestern Indian Ocean, respectively. (5) The Sr-Nd isotope distributions of continental flood basalts show HT picrite-basalt magmas mostly recording mantle values unaffected by continental lithospheric signatures, and LT basalts mainly reflecting either mixed source components located at the lithosphere-asthenosphere transition or continental crust contamination, particularly in the Karoo and Deccan provinces. Overall, results from this review provide compelling evidence that hot mantle plumes impinged diachronously on the central Gondwana lithosphere, causing similar tectonomagmatic events and continental flood basalt zonal arrangements that reflect a common thermocompositional zonation of the plume head in the three investigated provinces.

The Siberian craton was affected by flood basalt volcanism at least twice during the Devonian (Yakutsk-Vilyui province) and Permian–Triassic (Siberian province) periods. In both cases volcanism appeared as brief pulses of flood basalt eruptions, followed by kimberlitic (and lamproitic) emplacement. Pressure estimations for the kimberlite-entrained mantle xenoliths reflect that the lithosphere was 190–230 km thick at the time of the Devonian flood basalt volcanism. Differently from Devonian kimberlites, the majority of Triassic kimberlites are diamond free, but at least one Triassic kimberlite pipe and some lamproites are diamondiferous, suggesting that the Siberian lithosphere remained thick during the Permian–Triassic flood basalt volcanic activity. If both the lithosphere and the asthenosphere were volatile poor, thick cratonic lithosphere prevented melting even at an elevated geotherm. During the Paleozoic, Siberia was surrounded by subduction systems. The water deep cycle in association with fast subduction and slab stagnation in the mantle transition zone is proposed to cause fluxing of the asthenosphere by water plus other fluids via wet diapir formation in the mantle transition zone. Such diapirs started to melt in the asthenosphere beneath thick cratonic lithosphere, producing voluminous melts. Mafic melts probably accumulated beneath cratonic lithosphere and rapidly erupted on the surface in response to stress-induced drainage events, as assumed for some other cratonic flood basalts.

The Tuareg Shield, to which Hoggar (southern Algeria) belongs, has a swell-shaped morphology of lithospheric scale of ~1000 km in diameter linked to Cenozoic volcanism occurring in several regions, including Atakor, the center of the swell, which reaches nearly 3000 m in altitude. The lack of high-resolution geophysical data for constraining its deep structure is at the origin of a controversy about its innermost nature and about the origin of the Cenozoic volcanism. During the course of this study, magnetotelluric (MT) broadband data were collected at 18 sites forming a northeast-southwest profile 170 km long within the Atakor region. The electrical resistivity model obtained by inverting the magnetotelluric data reveals lithospheric structure down to a depth of ~100 km. From this depth to the surface, the model does not show any regional anomaly that may result from a metasomatized lithosphere or from an asthenospheric upwelling, including a mantle plume. MT data reveal rather a lithosphere affected by a set of rather thin subvertical conductors that can be attributed to the electrical signature of some known shear zones resulting from the Pan-African evolution of the LATEA metacraton, which globally corresponds to the uplifted Central Hoggar swell. The main anomaly is deeply rooted in the lithosphere and underlies the Atakor-Manzaz volcanic districts. As a whole, MT data are therefore properly integrated within the hypothesis of the reactivation of shear zones due the intraplate deformation related to the collision between Africa and Europe since the Eocene, applied onto the metacratonic region.

Abstract Based on present knowledge of mantle peridotites from the Ligurian Tethys ophiolites, this paper presents new ideas and a new model for passive rifting to ocean spreading in the slow–ultraslow rifting Europe–Adria realm. Relevant points include: (i) the positive feedback between deformation and melt percolation during passive magmatic rifting; (ii) the positive feedback between natural evidence and experimental data on the behaviour of the mantle lithosphere during passive rifting; (iii) the significance of hidden magmatism and the associated melt thermal advection; (iv) the role of the wedge-shaped weakened and softened axial zone; and (v) the evidence of a transition from passive to active rifting in the Ligurian Tethys. Passive rifting induced passive asthenospheric upwelling and the onset of partial melting. Fractional melts migrated through the mantle lithosphere and stagnated at shallow levels (the hidden magmatism). Melt thermal advection heated the mantle lithosphere to temperatures ( T ) of ≥1200°C and formed a wedge-shaped axial zone of rheological softened/weakened mantle peridotites that served as the future locus of continental break-up. The hotter/deeper asthenosphere ascended within this axial zone, underwent partial melting and formed aggregated mid-ocean ridge basalts (MORBs) that migrated within dunite channels to form olivine gabbro intrusions and basaltic lava flows. Rifting evolved from passive to active, and the actively upwelling asthenosphere established a ridge-type system and thermal regime.

Investigations of lithospheric foundering and related magmatism have long focused on the central Andes, where there are postulated links between the eruption of mantle-derived lavas and periodic loss of the lower lithosphere. Whole-rock elemental and Nd-Sr-Pb isotopic results from a suite of late Miocene–Quaternary mafic lavas erupted onto the Puna Plateau clarify the relationship between this hypothesized process and lava composition. Zinc and Fe provide a critical perspective because they are partitioned differently during the melting of asthenospheric and lithospheric mantle. All Puna lavas have Zn/Fe T (×10 4 ) values >13, which requires clinopyroxene and perhaps garnet to be the dominant phase(s) in the melt source; this precludes a melt source of typical mantle asthenosphere. This result is contrary to classic models of delamination magmatism that suggest asthenospheric peridotite melts to generate these lavas. Pyroxenite (±garnet)–bearing lithospheric materials in the central Andes are likely common and heterogeneous in age, volatile content, and mineralogical composition, and if they are the melt source, this can explain the diversity in the elemental (La/Yb = 11–45; La/Ta = 22–40) and isotopic ( 87 Sr/ 86 Sr = 0.7055–0.7080; ε Nd = −1 to −7) compositions of these mafic magmas (MgO > 8%, Mg number > 60). We propose that compositionally diverse, gravitationally unstable pyroxenites both drive “dripping” of the lower lithosphere and are the source of the resulting melt. We also postulate that mantle-derived lavas erupted on the Puna Plateau were generated during localized foundering and melting of these materials. The cumulative effect of these drip events is a modern Puna Plateau with geodynamic anomalies including thin lithosphere and anomalously high surface elevation.

Abstract Seaward dipping reflector or SDR complexes comprise piles of individual basaltic flows and interbedded sediments that are thought to have formed subaerially at the flanks of tholeiitic shield volcanoes like those in the Afar or at larger magmatic complexes like Iceland. During the rift process, these flows subsequently acquire very steep true dips (up to 25°) almost always in the seaward direction. Past explanations for the acquisition of these dips involve progressive burial, loading and flexure by subsequent flows ( e.g ., Pálmason, 1980 ), or listric, landward-dipping faulting and magmatic dilation ( Geoffroy, 2005 ). These factors no doubt play a role, but we feel that such models fall short of a full explanation by exceeding reasonable amounts of flexure and the amount by which huge blocks of continental crust can rotate by faulting alone. Such models also do not provide an explanation for how the topmost SDR layer subsides rapidly to the depth of normal oceanic crust as the latter begins to form, a problem that has been apparent since Mutter et al . (1982) and Hinz (1981) . Here, we present and discuss some of the main observations visible in long-offset, depth imaged seismic reflection records of SDR complexes at magma rich passive margins and propose a new but simple magmatic evacuation model for their production. Having conducted a global review of seismic data imaging SDR complexes, we present a line from the southern Brazilian margin to serve as a template example that shows most of the important criteria worldwide. After identifying the main observations, we propose a simple model of magmatic evacuation, similar to salt or mud evacuation in sedimentary sequences, to explain the observations. We hope that this simple proposal will spawn new avenues of research to refine and support the general model.

The Grande Ronde Basalt lavas constitute ~63% of the Columbia River Basalt Group, a large igneous province in the NW United States. The lavas are aphyric or contain less than 5% phenocrysts of plagioclase, augite, pigeonite, and olivine (altered). Plagioclase hygrometry shows that the erupted lavas generally contained less than 0.3% dissolved H 2 O; however, the presence of rare disequilibrium An 96 plagioclase phenocrysts suggests that some magmas may have originally had 4.5 wt% dissolved H 2 O at depth, but they all degassed during ascent and eruption. The size of plagioclase phenocrysts suggests an average plagioclase phenocryst residence time in the magmas of 160 yr. Ignoring hiatuses between eruptions, we estimate that the ~110 flows of the Grande Ronde Basalt erupted over a cumulative time of 17,600 yr, with an average eruption rate of ~8.6 km 3 /yr. The average interval between eruptions is estimated to be 3658 yr. It is envisaged that a shallow intrusive dike-sill complex, rather than large kilometer-sized magma chamber(s), fed the Grande Ronde basalt lavas. We performed model simulations using the COMAGMAT software to retrace the pre-eruption histories of the Grande Ronde Basalt lavas. Based on such simulations and petrological reasoning, we find that the primary melts could have been generated from a spinel peridotite source at 1.5 GPa, perhaps under hydrous conditions. Extensive melting of lithospheric eclogite may have played an important role as well; however, this is not constrained by our simulations. All lavas were contaminated by the crust, and they were last processed (mixing, crystallization) during their short residence within shallow (6 km) intrusive rocks prior to eruption. Our petrologic conclusions lead us to present a petrotectonic model that supports the hypothesis that the Columbia River Basalt Group magma generation was greatly aided by a thinned lithosphere and H 2 O that may have come off the asthenospheric wedge.

A variety of xenoliths from the lower crust to mantle transition occur in Quaternary mafic intraplate lavas of the Bayuda volcanic field of northern Sudan. The lower-crust xenoliths include plagioclase- and garnet-bearing mafic granulite. Ultra-mafic garnet-bearing pyroxenite, websterite, hornblendite, and distinct peridotite xenoliths are from the upper lithospheric mantle. Sr, Nd, and Pb isotope signatures distinguish between ultramafic and granulite xenoliths. The latter show a strong compositional affinity to juvenile Neoproterozoic crust. The Pb isotope composition of the ultramafic xenoliths resembles the distinct high-μ signature ( 206 Pb/ 204 Pb >19.5) of their host basanite. These xenoliths may represent cumulates of late Mesozoic to Quaternary mafic intraplate magmatism. The felsic upper crust in a schematic lithospheric profile of the Bayuda area includes predominantly granitoids, migmatites, and metasedimentary rocks that represent reworked old cratonic or juvenile Neoproterozoic rocks. The deep lower crust is represented by mafic granulite, likely cumulate rocks from Neoproterozoic juvenile magmatism. The crust-mantle transition is characterized by ultramafic cumulate rocks possibly from the late Mesozoic to Quaternary magmatism. The peridotites of the same xenolith suites represent typical lithospheric mantle with variable degrees of depletion by melt extraction.

The Atakor massif is a part of the Hoggar volcanic province, which was emplaced on top of a basement swell initiated during the Cretaceous. There have been three main episodes of volcanic activity since the Miocene, separated by long periods of quiescence. The lava flows and domes were emitted along lithosphere-scale fault zones. With its famous scenery, the Atakor massif is one of the largest (2150 km 2 ) volcanic districts of the province. Mafic volcanic rocks are abundant in the center of the massif, but become scarce to the south, where only few scarps are observed. Phonolites occur only in the Assekrem area, whereas trachytes occur everywhere, with a marked enrichment in quartz to the south and the southeast (Tahifet area), where rhyolites are also exposed. Two magmatic groups have been identified based on field and petrological observations. The mafic group has a basanite-phonotephrite association, forming uplifted plateaus, scoria cones, and valley-filling lava flows. The presence of mantle-derived amphibole ± biotite megacrysts and peridotite mantle xenoliths together with the nonprimary chemical compositions of the magmatic rocks suggest that magmatic differentiation may have occurred within the upper mantle. The felsic group is composed of two diverging trends, a silica-saturated benmoreite-trachyte-rhyolite trend and a silica-undersaturated trachyte-phonolite trend. The primary magmas are considered to have been produced as a consequence of lithospheric mantle delamination along linear megashear zones inducing low degrees of decompression partial melting at variable depths (110–40 km) in the upwelling asthenosphere. The discrete volcanic episodes correspond to periods of reactivation of the major fault zones in response to discrete Neogene extensional tectonic events associated with Alpine orogenesis in the Western Mediterranean region induced by Africa-Eurasia collision.

The Taconic orogeny in Newfoundland consisted of three accretionary events (Taconic 1, 2, and 3). Taconic 1 is represented by ca. 495 Ma, west-directed obduction of the infant-arc Lushs Bight oceanic tract (510–501 Ma) onto the peri-Laurentian Dashwoods microcontinent. Subduction is inferred to have initiated at a spreading center abandoned during an inboard ridge jump responsible for separation of Dash-woods from Laurentia and opening of the Humber seaway. Clogging of the subduction zone by Dashwoods forced subduction to step back into the Humber seaway. Inception of the new subduction zone led to formation of the ca. 490 Ma Baie Verte oceanic tract. Closure of the Humber seaway formed the Notre Dame arc (489–477 Ma) built on Dashwoods and the coeval Snooks Arm arc built on the Baie Verte oceanic tract. Sea-way closure led to collision (Taconic 2) between the arcs and Laurentia, which caused significant shortening of the Notre Dame arc. After a magmatic gap of 7–10 m.y., the Notre Dame arc records a voluminous flare-up of predominantly tonalite magmatism (464–459 Ma) during the waning stages of Taconic 2. Magmatism overlaps with deformation and includes both arc and non-arc-like tonalite. This flare-up was related to break-off of the oceanic lithosphere of the downgoing slab. The rapidly upwelling asthenosphere that replaced the broken-off slab induced melting in the subarc mantle and arc infrastructure. Taconic 3 is represented by 455–450 Ma accretion of a peri-Laurentian arc that had formed after the ca. 480 Ma initiation of west-directed subduction in the Iapetus Ocean outboard of the Dashwoods microcontinent.

A conceptual shift is overdue in geodynamics. Popular models that present plate tectonics as being driven by bottom-heated whole-mantle convection, with or without plumes, are based on obsolete assumptions, are contradicted by much evidence, and fail to account for observed plate interactions. Subduction-hinge rollback is the key to viable mechanisms. The Pacific spreads rapidly yet shrinks by rollback, whereas the subduction-free Atlantic widens by slow mid-ocean spreading. These and other first-order features of global tectonics cannot be explained by conventional models. The behavior of arcs and the common presence of forearc basins on the uncrumpled thin leading edges of advancing arcs and continents are among features indicating that subduction provides the primary drive for both upper and lower plates. Subduction rights the density inversion that is produced when asthenosphere is cooled to oceanic lithosphere: plate tectonics is driven by top-down cooling but is enabled by heat. Slabs sink more steeply than they dip and, if old and dense, are plated down on the 660 km discontinuity. Broadside-sinking slabs push all sublithosphere oceanic upper mantle inward, forcing rapid spreading in shrinking oceans. Down-plated slabs are overpassed by advancing arcs and plates, and thus transferred to enlarging oceans and backarc basins. Plate motions make sense in terms of this subduction drive in a global framework in which the ridge-bounded Antarctic plate is fixed: most subduction hinges roll back in that frame, plates move toward subduction zones, and ridges migrate to tap fresh asthenosphere. This self-organizing kinematic system is driven from the top. Slabs probably do not subduct into, nor do plumes rise to the upper mantle from, the sluggish deep mantle.

Asthenosphere plume-to-ridge flow has often been proposed to explain both the existence of geochemical anomalies at the mid-ocean ridge segments nearest an off-axis hotspot and the existence of apparent geochemical provinces within the global mid-ocean spreading system. We have constructed a thin-spherical-shell finite element model to explore the possible structure of global asthenosphere flow and to determine whether plume-fed asthenosphere flow is compatible with present-day geochemical and geophysical observations. The assumptions behind the physical flow model are described in the companion paper to this study. Despite its oversimplifications (especially the steady-state assumption), Atlantic, Indian, and Pacific mid-ocean ridge isotope geochemistry can be fit well at medium and long wavelengths by the predicted global asthenosphere flow pattern from distinct plume sources. The model suggests that the rapidly northward-moving southern margin of Australia, not the Australia-Antarctic discordance, is the convergence zone for much plume material in the southern hemisphere. It also suggests a possible link between the strike of asthenosphere flow with respect to a ridge axis and along-axis isotopic peaks.

Anorogenic magmatism of the circum-Mediterranean area (the Tyrrhenian Sea, Sardinia, Sicily Channel, and the Middle East) and of continental Europe (the French Massif Central, Eifel, the Bohemian Massif, and the Pannonian basin) has been proposed to be related to the presence of one or more mantle plumes. Such conclusions based on geochemical data and seismic tomography are not fully justified because (1) a given chemical and isotopic composition of a magma can be explained by different petrogenetic models, (2) a given petrogenetic process can produce magmas with different chemical and isotopic composition, (3) tomographic studies do not furnish unique results (i.e., different models give different results), and (4) the commonly adopted interpretation of seismic wave velocity anomalies exclusively in terms of temperature is not unique; velocities are also dependent on other parameters, such as composition, melting, anisotropy, and anelasticity. Tomography and geochemistry are powerful tools but must be used in an interdisciplinary way, in combination with geodynamics and structural geology. Alone they cannot provide conclusive evidence for or against the existence of mantle plumes. The existence of large and/or extensive thermal anomalies under Europe is here considered unnecessary, because other models, based on the existence of upper-mantle heterogeneity, can explain the major-element, trace-element, and isotopic variability of the magmas. Volcanism in central Europe (the French Massif Central, Germany, and the Bohemian Massif) is concentrated in Cenozoic rifted areas and is here interpreted as the result of passive asthenosphere upwelling driven by decompression. Similarly, anorogenic magmatism in Sardinia, the Tyrrhenian Sea, and the Pannonian basin is explained as the result of lithospheric stretching in a back-arc geodynamic setting. The most important factors determining the locus and, in part, the geochemical characteristics of magmatic activity are the Moho and the lithosphere-asthenosphere boundary depths. Where both are shallowed by tectonic processes (e.g., in rift zones or back-arc basins), passive upwelling of asthenospheric mantle can explain the magmatic activity.

Neogene to Quaternary volcanism of the Carpathian-Pannonian region is part of the extensive volcanic activity in the Mediterranean and surrounding regions. Using the spatial and temporal distribution of the magmatic rocks, their major- and trace-element features, and Sr-Nd-Pb isotope characteristics, we suggest that lithospheric extension in the Pannonian Basin had a major role in the generation of the magmas. Dehydration of subducting slab should have resulted in thorough metasomatism in the mantle wedge during Cretaceous to early Miocene that would have lowered the melting temperature, therefore playing an indirect role in the generation of magmas later on. Mixing between mantle-derived magmas and lower-crustal melts was an important process at the first stage of the silicic and calc-alkaline magmatism in the Northern Pannonian Basin. However, the crustal component gradually decreased with time, which is consistent with magmatic activity in a continuously thinning continental plate. Calc-alkaline volcanism along the Eastern Carpathians was mostly postcollisional and could have been related to slab break-off processes. However, the fairly young (<1.5 Ma) potassic magmatism at the southeasternmost segment of the Carpathian volcanic arc could be explained by lithospheric delamination under the Vrancea zone. Alkaline basaltic volcanism began at the end of rifting of the Pannonian Basin (11 Ma) and continued until recently. We suggest that a mantle plume beneath the Pannonian Basin is highly unlikely and the mafic magmas were formed by small degree partial melting in a heterogeneous asthenospheric mantle, which has been close to the solidus temperature due to the lithospheric extension in the Miocene. Magmatism appears to have been in a waning phase for the last 2 m.y., but recent volcanic eruptions (<200 k.y.) indicate that future volcanic activity cannot be unambiguously ruled out.