Mantle Plume Research Articles

Upper mantle shear velocity structure of the Cathaysia Block and surrounding areas: New insight into deep geodynamics

The Cathaysia Block (CAB) and its surrounding areas experienced intensive magmatism and mineralization in the Yanshanian period (ca. 200–90 Ma), but their mechanism and deep geodynamics are still debated. In addition, the origin and structures of the Hainan mantle plume are still unclear. To resolve these issues and investigate the large-scale lithospheric thinning and extension in the eastern South China Block, we determine a detailed 3-D S-wave velocity (Vs) model down to 700 km depth by collecting 24,190 S-wave teleseismic data recorded at 164 permanent stations and 125 portable stations deployed in the CAB and surrounding areas. Our results show that high-Vs anomalies exist separately in the study volume. Two high-Vs anomalies appear in the shallow upper mantle and the mantle transition zone, which may reflect the present thin lithosphere and the stagnant Paleo-Pacific slab, respectively. Two other high-Vs anomalies exist in the upper mantle, which may reflect the detached lithosphere and subducting slabs. In contrast, low-Vs anomalies appear widely beneath the CAB, which reflect a tilting magmatic conduit beneath the Wuyishan metallogenic belt (WYMB) and magmatic chambers beneath the Nanling metallogenic belt (NLMB). In addition, our results show that the Hainan plume has a double-layered appearance. Combining our tomographic results with previous multidisciplinary findings, we consider that (1) the subduction and rollback of the Paleo-Pacific Plate may have played different roles in the Yanshanian mineralization of the WYMB and the NLMB; (2) the double-layered appearance of the Hainan plume may be formed due to the influence of plume self-evolution dynamics and pre-existing deep structures; and (3) lithospheric delamination triggered by gravity instability may have occurred beneath the Xuefengshan Mountain in the late Mesozoic.

Zircon Xenocrysts From Easter Island (Rapa Nui) Reveal Hotspot Activity Since the Middle Jurassic

AbstractWe report the finding of mantle‐derived zircon grains retrieved from red soils, regoliths, and beach sands from Easter Island, that are much older than the island volcanism (0–2.5 Ma) and its underlying lithosphere (Pliocene, 3–4.8 Ma). A large population of 0–165 Ma old zircons have coherent oxygen (δ18O 3.8–5.9‰) and hafnium (εHf(t)+3.5–+12.5) mantle isotopic signatures. These results are consistent with the crystallization of zircon from plume‐related melts. In addition, a chemically diverse population with ages from the Paleozoic to the Archean was found. These older populations are enigmatic but they could represent remnants of ancient subducted sediments. Meanwhile, the ∼0–165 Ma population is interpreted as plume‐derived, suggesting that the hotspot started at least ∼165 Ma ago. A spike of ∼164–160 Ma zircons could represent a Large Igneous Province (LIP) stage upon the first arrival of the plume. We use plate reconstructions to show that such a LIP would have formed on the Phoenix Plate and would have subducted below the Antarctic Peninsula around 100–105 Ma. There, LIP subduction would offer a solution for the enigmatic Palmer Land deformation event, previously proposed to result from a collision with an unknown indenter. The here‐reported “ghost” of a prolonged hotspot activity suggests that fragments of the Easter plume and of the ambient sub‐lithospheric mantle stored and re‐sampled zircon xenocrysts due to convective (re)circulation at the scale of the plume head. Our study demonstrates how zircon geochronology and geochemistry provide novel insights into global‐scale geodynamics, offering new perspectives on the dynamics of mantle plumes and hotspot activity.

Deflected Mantle Flow and Shearing‐Aligned Lithospheric Melt Under the Strike‐Slip Dead Sea Rift

AbstractContinental rifting is one of the fundamental tectonics of the Earth evolution while our current understandings on the dynamic mechanism of the strike‐slip ones are relatively limited. Here, we have utilized three kinds of core‐refracted shear waves (including PKS, SKKS, and SKS) and employed the shear‐wave splitting technique to systematically investigate the azimuthal anisotropy of the upper mantle under the typical strike‐slip Dead Sea rift. There are a total of 1,855 well‐determined anisotropic measurements from 187 stations with dominantly N‐S fast orientations. We have proposed a new model from a joint analysis of multiple newly available geophysical observations to interpret the resulting anisotropy as mainly due to the absolute‐plate‐motion‐driven mantle flow deflected by the thick lithosphere of the eastern Arabian plate. Relatively larger splitting times are locally revealed at the rift zone and attributed to additional lithospheric anisotropy from the shearing‐oriented melt pockets whose existence further induces complex anisotropy with slight difference of fast orientations between the upper and lower layer anisotropy. The overall rift‐parallel fast orientations, when combined with the absence of low‐velocity and hot thermal anomalies in the mantle transition zone, rule out the role of mantle plume or edge‐driven convection in the rift development and further infer the Dead Sea rift to evolve in a passive mode.

The geodynamics of plume-influenced mid-ocean ridges: insights from the Foundation Segment of the Pacific-Antarctic Ridge

The intersection of the Foundation Plume with the Pacific-Antarctic Ridge is a key location in global geodynamics where a mantle plume is approached by and interacting with a fast-spreading mid-ocean ridge. Here, we discuss a comprehensive major and trace element and Sr-Nd-Pb isotope dataset of new and existing samples from the young Foundation Seamount Chain (<5 Ma) and adjacent section of the Pacific-Antarctic Ridge. We use the geochemistry of axial, off-axis and intraplate lavas to map the spatial extent of plume dispersal underneath the ridge as well as the internal zonation of the upwelling plume. We show that the unusual length, increased crustal thickness and occurrence of silicic rocks on the axis of the Foundation Segment are the direct result of plume being tapped by the axial melting zone. We demonstrate that the plume is not homogeneous but shows a HIMU-like (high time-integrated 238U/204Pb) OIB (Ocean Island Basalt) component characterized by 206Pb/204Pb of up to 20.5 in its center and a more EM1-like (Enriched Mantle one) OIB component characterized by low U/Pb and 206Pb/204Pb but high Rb/Nb and 87Sr/86Sr towards its edges. Plume entrainment leads to a high magma supply rate that fosters the formation of silicic rocks and triggers the lengthening of the segment over time. However, plume dispersal is not symmetric as the geochemical tracers for the OIB component are extending <100 km northwards but >300 km southwards. We relate this to the current plate tectonic framework in which the obliquity between the migrating ridge and the absolute plate motions induces a sub-axial asthenospheric flow that preferentially channels plume material southwards.

Petrogenesis and an Evaluation of the Melting Conditions of the Late Permian ELIP Picrites, SW China: Constraints Due to Primary Magma and Olivine Composition

The late Permian Emeishan large igneous province (ELIP) in SW China is a melting product of the Emeishan mantle plume. Recently, it has been debated whether peridotite or pyroxenite is the dominant lithology of the mantle source in the ELIP. To address this, systematic analyses of bulk-rock and coexisting spinel and olivine compositions were conducted on picrites from Lijiang–Yongsheng, Dali–Binchuan, Yumen, Muli, and Ertan. The ELIP picrites exhibit positive TiO2–CaO and negative MgO–CaO correlations, as well as low FC3MS values (−0.24–0.1), supporting a peridotite-dominated mantle source. This lithology of the mantle source is also supported by the high 100 × Mn–Fe (1.43–1.73) and Mn–Zn (13.6–18.4) values but low 10,000 × Zn–Fe (8.0–12.7) ratios of the olivine phenocrysts. The estimated mantle potential temperature for Lijiang, Yongsheng, Yumen–Ertan, Muli, and Dali–Binchuan picrites decreased away from Lijiang and Yongsheng, suggesting that the Lijiang and Yongsheng areas were the center of the ELIP. The Lijiang–Yongsheng primary magma shows similar SiO2 content but lower Al2O3 contents (average of 8.24 wt.%) and higher MgO contents (average of 21.42 wt.%) than those of Dali–Binchuan primary magma (Al2O3: 9.86 wt.%; MgO: 19.02 wt.%). Also considering the high Gd–Yb (average of 3.05) and La–Yb (average of 14.61) ratios and mantle potential temperature (average of 1599 °C), we proposed that Lijiang–Yongsheng lavas are produced via the melting of a garnet–peridotitic mantle. In contrast, the Dali–Binchuan lavas with low Gd–Yb (average of 1.91) and La–Yb (average of 5.88) ratios can be explained by their formation in the garnet–spinel transition zone of a peridotitic mantle. The Yumen–Ertan primary magma displays similar mantle potential temperature (average of 1600 °C), Al2O3 and FeO content, and Gd–Yb ratios to those of Lijiang–Yongsheng lavas, indicating that YumenvErtan primary magma may be attributed to the partial melting of garnet with minor peridotite. Therefore, heterogeneous plume-head mantle sources lead to the evaluation of melting conditions of the late Permian ELIP picrites.

МОРФОЛОГИЯ ОБЛАСТИ СОЧЛЕНЕНИЯ ХРЕБТА РЕЙКЬЯНЕС С РИФТОВЫМИ ЗОНАМИ ИСЛАНДИИ

The joint zone between subaquatic Reykjanes Ridge and Icelandic rift zones known as Reykjanes Rift Zone is a unique transform zone formed in conditions of interaction between the mantle plume and the mid-oceanic ridge. Its morphostructures include a row of en-echelon volcanic systems with oblique extension. The modern tectonic and volcanic relief was formed almost completely during Holocene: the previous landforms were reworked or destroyed by erosion of the Late Pleistocene glaciation. Morphology of the Reykjanes Rift Zone is transitional between the Reykjanes Ridge and the adjacent Western Rift Zone from the northeast. Unlike the Reykjanes Ridge, volcanic activity is not concentrated with in axial volcanic ridges, but is confined by eruptive centers represented by fissure volcanoes. Central volcanoes, in contrast to the Western Rift Zone, are not expressed in topography. Oppositely, tectonic activity increases with distance from the Reykjanes Ridge. Fault scarps reach their largest size within the Western Rift Zone, where Thingvellir graben is located. At the same time, the tectonic and volcanic landforms of the Reykjanes rift zone are in paragenetic association, appearing sequentially during the episodes of riftogenesis.

The origin and implications of primordial helium depletion in the Afar mantle plume

Mantle plumes are responsible for the Earth’s largest volcanic provinces. In the prevailing paradigm, the deep mantle is less degassed than convecting shallow mantle, implying that plume-derived lavas have higher concentrations of primordial volatiles such as helium (He). Demonstrating this has led to explanations that question the established Earth model. Here, we show that the 3He/4He of basalts from the Red Sea display coherent relationships with trace elements, allowing the helium concentration of the Afar plume to be calculated. Contrary to the prevailing model it appears the helium concentration of the Afar plume is 10-25% of the upper mantle. This contradiction is resolved if the plume material itself is a mixture of helium-rich high-3He/4He deep mantle with helium-depleted low-3He/4He recently subducted oceanic crust. This implies that helium-depleted domains may exist in convecting mantle and that moderately high 3He/4He plumes likely do not contain a notable contribution of the deep mantle.

Geophysical characterization and 2D gravimetric modeling: Application to tectonic control of the Vitória-Trindade Ridge (VTR–ES, Brazil)

There are numerous inquiries regarding the origin and development of the Vitória-Trindade Ridge (VTR) in eastern Brazil, a major tectonic feature composed of seamounts, guyots, banks, and islands, prompting extensive studies in the region. These investigations aim to create a more precise crustal model to aid in understanding various geological phenomena, such as the separation of South America from Africa, the interaction between continental and oceanic crusts, and the behavior of oceanic fracture zones. These zones may have served as conduits for magma feeding into the ridge and adjacent areas, or the seamounts could have formed due to a fixed mantle plume related to the West-East trend. This study focuses on exploring the tectonic structure of the region, particularly the subsurface of the VTR, utilizing geophysical data modeling techniques. This includes data from topography, gravimetric and magnetic anomalies, and analysis of reflection seismic lines. The research aims to accurately characterize the geological and geophysical context of the continental margin and some igneous features in the offshore region of the Espírito Santo and Bahia states. The anomalous hotspots correlate with pre-existing fracture zones resulting from the Mesozoic rifting of South American and African continents, which due to later tectonic episodes, have masked traces of deep-seated tectonic structures. The integration of geophysical methods, particularly gravimetry and magnetometry, when combined with reflection seismic data and 2D modeling techniques, offer insights into the geodynamic processes involved in the formation of the volcanic ridge. The results obtained showed the interference generated by magmatism in the oceanic crust and the transition zone. The seismic data allowed the visualization of the approximate geometries and interrelationships of some volcanic edifices that are compared with the Davis Bank, located in the central region of the VTR.

Petrogenesis of the Deccan high-Mg basalts and picrites

Tholeiitic basalts and picrites from the Deccan Traps were used to constrain the pressure and temperature conditions of mantle melting for their origin. Clinopyroxene thermobarometry indicates that all Deccan tholeiites crystallized at low pressures in the upper crust (< 6 kbar/1047–1221 °C). In comparison, the Deccan alkalic rocks crystallized at pressures up to ~ 12.7 kbar. Rare samples of the tholeiites plot on their low-pressure olivine-plagioclase-clinopyroxene (Ol-Pl-Cpx) cotectic boundaries or olivine control lines in phase diagrams. These samples represent unmodified magmatic liquids. Primary magmas of the basalts that plot on their cotectic boundaries were modeled through reverse fractionation by incrementally adding equilibrium Ol + Pl + Cpx, Ol + Pl and Ol ± spinel, until the liquid was multiply saturated with lherzolite at a high pressure. The high-Mg basalts are contaminated with continental crust. Hence, a crustal partial melt was simultaneously subtracted according to energy constraints at each reverse fractionation step for these samples. The results show that the high-Mg basalts are 41–53% fractionated and 1–6% contaminated, and the low-Mg basalts are 63–67% fractionated. Their primary magmas were last equilibrated with spinel lherzolite at 10–13 kbar/1289–1333 °C. A picrite and two very high-Mg basalts plot on their olivine control lines. So, their primary magmas were calculated by adding only equilibrium olivine. These samples are 9–25% fractionated, and their primary magmas were last equilibrated with garnet lherzolite at 25–36 kbar/1452–1531 °C. The estimated mantle potential temperatures are 1400–1500 °C for the Deccan tholeiites, consistent with their origin from a mantle plume.

Geochemistry and conditions of formation of Mesoarchean banded iron formations (BIF-1) from the Kostomuksha Greenstone belt, Karelian Craton

Three variably old groups of banded iron formation (BIF) are known in the Kostomuksha Greenstone belt (KGB) of the Karelian Craton. This paper deals with the earliest of them, Mesoarchean (2.87–2.81 Ga) – BIF-1. BIF-1 occurs among the komatiite-basalt unit of the KGB. BIF-1 consists mainly of quartz and magnetite, with varying amounts of amphibole, biotite, and garnet; they contain 48.3-58.6 SiO2 and 21.34–33.82 wt. %, Fe2O3T, suggesting that the rocks are BIF. BIF-1 of the KGB, as well as most Archean BIFs, contain high Fe2O3T, concentration, display a contrasting positive Eu anomaly, lost of Ce anomaly, the depletion of LREE relative to HREE. However, they stand out among other BIFs with high Al2O3, TiO2, MgO, K2O, Cr, Ni, Zr, Ba, Cu and Zn concentrations. BIF-1 was formed in a marine basin in an anoxic atmosphere due to hydrothermal fluids, the proportion of which varies from 20 to 80 %, and a terrigenous component derived mainly from basalts, komatiites, and dacites in host rocks. Mesoarchean BIF-1 of the KGB s was formed in a small rift within an oceanic volcanic plateau, the formation of which is associated with the influence of a mantle plume on the oceanic lithosphere.

Marginal stability analyses for thermochemical convection and its implications for the dynamics of continental lithosphere and core–mantle boundary regions

SUMMARY Significant compositional differences may exist in the lithospheric mantle and above the core–mantle boundary (CMB) relative to the ambient mantle. The intrinsic density differences may affect the development of thermal boundary layer (TBL) instabilities associated with lithospheric delamination and formation of thermochemical plumes. In this study, we explored the instability of two-layer thermochemical fluid using two different techniques: marginal stability analysis with a propagator-matrix method and finite element modelling. We investigated both the instabilities in lithospheric mantle (i.e. lithospheric instability) and the mantle above the CMB (i.e. plume-forming instability) using a background temperature Tbg(z) with the TBL. For lithospheric instability, we found that two-layer fluid with free-slip boundary conditions mainly undergoes the same three different convective modes (i.e. two oscillatory convection modes and one layered convection regime) as that with no-slip boundary conditions reported in Jaupart et al. However, with free-slip boundary conditions, the transitions between these convection modes occur at larger values of buoyancy number B. Free-slip boundary conditions lead to smaller critical Rayleigh number Rac, but larger convective wavelength and oscillation frequency ωc, compared with those with no-slip boundary conditions. Our numerical modelling results demonstrate that Rac and ωc predicted from the classical marginal stability analyses using Tbg(z) with TBL temperature may have significant errors when the oscillatory period is comparable with or larger than the timescale of lithospheric thermal diffusion that causes Tbg(z) to vary with time significantly. In this case, using a more gently sloped background temperature profile ignoring the TBL temperature, the stability analysis predicts more accurate stability conditions, thus presenting an effective remedy to the stability analysis. For plume-forming instability, because of the reduced viscosity in the hot and compositionally dense bottom layer, the transition to the layered convection occurs at significantly smaller B values, and in the oscillatory convection regime, Rac is larger but ωc is smaller, compared with those for lithospheric instability. Finally, our study provides a successful benchmark of numerical models of thermochemical convection by comparing Rac and ωc from numerical models with those from the marginal stability analysis.

Review of morphotectonic studies in the East African Rift System: Insights from legacy research to modern advances

The East African Rift System (EARS) is a prime example of continental rifting, characterized by its iconic topographic features, which are significantly shaped by the interplay of tectonic and climatic forces. This review synthesizes past decades of morphotectonic studies (recent studies) on the plateau landscapes of the EARS, focusing particularly on their interiors and margins, and highlights their crucial role in advancing our understanding of landscape evolution. The paper revisits the legacy of geologic research in the EARS and delves into recent morphotectonic studies that quantify the rates and patterns of plateau incision, exhumation, and uplift through advanced thermochronological and topographic analyses. Additionally, the review identifies key unresolved questions, such as the origins of the plateau topography–whether from volcanic activity or mantle plume influence and explores how these plateaus maintained high elevations during the humid Oligocene Epoch without significant incision. The recent work described in this paper alongside legacy perspectives, continues to refine our geological understanding of the EARS, with broader implications for studies worldwide.

The Gulf of Suez rifting: implications from low-temperature thermochronology

ABSTRACT The Gulf of Suez is a continental rift basin that has been developed during the Oligocene-Miocene as a northern extension of the Red Sea rift system. The Oligo-Miocene rift-activated flanks are represented by exposures of the Neoproterozoic basement rocks of the Arabian-Nubian Shield. An additional thermal overprint (i.e. mantle plume) and the preceding tectono-thermal events that also affected the rifting dynamics and the uplift of its flanks remain uncertain. We here present thermochronological data for 20 basement samples collected from the Gabal Somr ElQaa area at the northwestern flank of the Gulf of Suez. The resulted zircon fission-track ages could be divided temporospatially into two groups: the Cambrian and the Silurian-Devonian with average ages of 503 ± 14 Ma and 407 ± 11 Ma, respectively. While apatite fission-track data reveal three different age groups of the Devonian-Carboniferous, the Permian-Triassic, and the Late Cretaceous, with average ages of 349 ± 10 Ma, 246 ± 6 Ma, and 81 ± 5 Ma, respectively. The modelled t-T histories indicate four cooling pulses during the Neoproterozoic, the Devonian-Carboniferous, the Cretaceous, and the Oligocene-Miocene. Integrating our findings with the regional tectonism and depositional records indicates that these cooling pulses could be attributed to synchronous activities of the post-accretion erosional event, the Variscan tectonic event, the Mid-Atlantic Ridge spreading, and the Gulf of Suez rifting. The Gulf of Suez is a passive rift, formed by a single and rapid rifting pulse within the framework of the Red Sea rift system. An additional east-southward far-field thermal overprint from the Arabian marginal plume partially affected the Suez rifting process. This thermal print affected only the rift’s southern portion of its eastern flank, causing syn-rift cooling ages, highly elevated rift flanks, and an increased heat flow. The western rift flank and the northern segment of the eastern flank are characterized by pre-rift cooling ages, modest rift flanks, and reduction in heat flow.

Mechanisms of Co enrichment in the Zhubu Ni–Co–PGE deposit in the Emeishan LIP, SW China: Constraints from mineral geochemistry and Fe–S isotopes of sulfides

The Zhubu Ni–Co–PGE sulfide deposit in the Permian Emeishan Large Igneous Province (LIP), SW China is hosted by a mafic–ultramafic intrusion that comprises two morphologically distinct ore-forming spaces: a central, layered sequence that represents a magma chamber surrounded by a marginal sub-vertical conduit zone. The marginal conduit zone contains significant Ni–Co–PGE sulfide mineralization, suggesting that high volumes of sulfide-bearing silicate melt flowed through the magma conduit system. Determining the hosts of Co is key to understanding the Co enrichment process in the Zhubu deposit, as well as in the Emeishan LIP as a whole. We assess this here using the geochemistry and Fe–S isotope compositions of sulfides in the Zhubu deposit. Cobalt is mainly concentrated in pentlandite sulfide in the websterite, with Co contents increasing in the order pyrrhotite < chalcopyrite < pyrite < pentlandite < violarite; the abundance of violarite, however, is notably lower than that of the other sulfides. The Co contents of pentlandite, pyrite, and chalcopyrite are negatively correlated with their Ni and Fe contents, suggesting that Co substitutes isomorphically for Ni and Fe in the sulfide structure. The δ56Fe values of sulfides increase in the order pyrrhotite < chalcopyrite < pentlandite < pyrite resulted from kinetic fractionation of Fe isotopes. The δ56Fe values of pentlandite and chalcopyrite in lherzolite are higher than those in websterite, implying that lherzolite and websterite formed in compositionally distinct magmas, Fe isotope compositions become heavier with crystal fractionation. The δ34S values of sulfides increase in the order pyrrhotite < pentlandite < chalcopyrite < pyrite, similar to the fractionation sequence of sulfide melts, suggesting that externally derived S was added to the magma and contributed to its sulfide saturation. The Co contents of chalcopyrite and pyrite vary positively with δ56Fe values, and negatively with δ34S values, indicating that the addition of externally derived S and resulting sulfide segregation were critical to Co enrichment.

Electrical conductivity of mantle minerals beneath East Asia revealed by geomagnetic observatory data

The electrical conductivity of the earth’s mantle can provide important information about geodynamic phenomena. East Asia is associated with complex tectonics and geodynamic processes. Hence, it is necessary to better understand the deep structure beneath East Asia. In this study, geomagnetic data obtained from East Asian observatories are employed to image the conductivity structure of the mantle at depths ranging from 410 to 900 km. First, the data are processed using the modified bounded influence remote reference processing (BIRRP) method and the ratio method is used to correct for the ocean effect. Thereafter, the stable C-response curves at the 27 observatories are estimated, and 1D electrical conductivity models for these observatories are established using the L-BFGS method. The conductivity-depth profiles reveal a heterogeneous distribution of the electrical conductivity beneath East Asia. The mantle transition zone (MTZ) beneath East China and Japan is found to be more conductive, whereas the MTZ beneath central and southern regions of China is more resistive. In East China, the dehydration of the stagnant Pacific slab may lead to an increase in the conductivity of the mantle minerals. There is also the possibility of upwelling of the thermal material from the lower mantle beneath the Japanese Island arc. In Northwest China, there exists a large high-conductive body beneath the Tarim area, which could indicate an upwelling of the Tarim mantle plume. Our results provide insights into the deep structure of the earth at the mantle scale.

Anatomy of the Emeishan Mantle Plume Head: Insights From New Geochronologic, Geochemical, and Geologic Data

AbstractThe link between mantle plumes and the formation of large igneous provinces (LIPs) is well established although the anatomy of these remains equivocal. Recent experimental studies and geophysical data suggest that the mantle plume head is more likely to be irregular and asymmetric, rather than an axisymmetric flattened disk. The Emeishan large igneous province (ELIP) provides a unique opportunity to test this hypothesis. According to robust petrographic, geochronologic, and geochemical evidence from the late Permian basalts in the Sichuan Basin, and in conjunction with a comprehensive compilation of geologic maps and published geochemical data from the ELIP, we identified several giant radial “fingering” structures. Based on the shallow mantle source from the center to margin in the ELIP and relief of the lithosphere‐asthenosphere boundary, we propose a new mantle plume model to explain the evolution of the Emeishan plume periphery, where narrow finger‐like protrusions and plumelets developed outwards from the main body of the plume to the edges of the flattened plume head. Dragged fingers might have been torn apart into some plumelets, which dispersed and were trapped beneath the thinnest lithosphere relief, and eventually erupted to form small‐scale flood basalt in the Outer Zone of the ELIP.

Tectonic evolution of the Early Paleozoic proto-East Asian continental margin: Inference from Paleozoic metamorphic rocks in the South Kitakami belt, northeast Japan

To constrain the incipient Pacific-type orogeny and tectonic processes in the Early Paleozoic proto-East Asian continental margin, the Motai–Matsugadaira–Yamagami (MMY) metamorphic rocks in the South Kitakami belt, northeast Japan were investigated. They are divided into two different types: amphibolite-facies rocks associated with serpentinite and blueschist-facies rocks associated with pelitic and psammitic schists. Three geochemical groups are identified from the MMY metamorphic rocks. Groups 1 and 2 resemble geochemical characteristics of mid-ocean ridge basalt and continental arc rocks, respectively. Group 3 exhibits considerable depletion of highly incompatible elements, which is caused by the high degree of partial melting of a hot mantle plume. The zircon U–Pb ages of Group 1 indicate that the protoliths experienced amphibolite-facies metamorphism soon after their formation in the Early Ordovician. Group 2 exhibits a coeval zircon U–Pb age with Group 1. The age distribution of detrital zircons in the MMY psammitic schists shows a peak of 500–400 Ma, the presence of Archean to Neoproterozoic zircons, and the youngest Late Devonian zircon. The following model is proposed for the tectonic evolution of the proto-East Asian continental margin: (1) the formation of an arc in the eastern margin of the South China craton in the Cambrian to Ordovician; (2) the subduction of a spreading axis and an oceanic plateau at the same time as the continental arc formation; (3) the consumption and subduction of arc materials by tectonic erosion; and (4) the formation of the Carboniferous accretionary complex and high-pressure metamorphic rocks under steady oceanic plate subduction. The proposed tectonic evolution model may also be applicable to equivalent Early Paleozoic rocks in southwest Japan.

Craton deformation from flat-slab subduction and rollback

The mechanisms underlying the deformation and eventual destruction of Earth’s cratons remain enigmatic, despite proposed links to subduction and deep mantle plume processes. Here we study the deformation of the North China Craton using four-dimensional mantle flow models of the plate–mantle system since the late Mesozoic, integrating constraints from lithospheric deformation, mantle seismic tomography and the evolution of surface topography. We find that flat-slab subduction induced landward shortening and lithospheric thickening, while subsequent flat-slab rollback caused seaward extension and lithospheric thinning. Both subduction phases resulted in substantial topographic changes in basin sediments. Rapid flat-slab rollback, coupled with a viscosity jump and phase change across the 660 km mantle discontinuity, was a key ingredient in shaping a large mantle wedge. We argue that craton deformation through lithospheric extension and thinning was triggered by the subduction of a flat slab and its subsequent rollback. The integration of data into mechanical models provides insights into the four-dimensional dynamic interplay involving subduction, mantle processes, craton deformation and topography.

CSRM‐1.0: A China Seismological Reference Model

AbstractHigh‐resolution seismic model is crucial for advancing our understandings on geological processes and enhancing seismic hazard mitigation programs. We construct a high‐resolution China Seismological Reference Model (CSRM‐1.0) in the top 100 km of the crust and uppermost mantle in continental China following a top‐down construction process. The employed seismic constraints include P‐wave polarization angle from tele‐seismic event, short‐period Rayleigh wave ellipticity from ambient noise, long‐period Rayleigh wave ellipticity from earthquake data, receiver function, empirical Green's function from ambient noise, Rayleigh wave phase/group velocity dispersion curves from regional earthquakes, and Pn‐wave travel time extracted from seismic data of 4,435 stations. CSRM‐1.0 has a spatial crustal resolution of ∼60 km beneath the north‐south seismic belt and trans‐North China orogen regions and ∼120 km beneath the rest of continental China, and a spatial mantle resolution of ∼300 km. CSRM‐1.0 exhibits prominent velocity heterogeneities in the crust and uppermost mantle and an eastward thinning of the crust, geographically correlating with geological settings. CSRM‐1.0 improvements include accurate estimation of shallow seismic structure, increased spatial resolution and improved model accuracy. Crustal composition inferred from CSRM‐1.0 exhibits a general transition from a felsic upper crust to a mafic lower crust. Mafic rocks in the lower crust are found predominantly along inter‐block boundaries and sporadically within the interiors of blocks, likely resulted from preferential inter‐block intrusions of magmas related to various oceanic plate subductions and the Emeishan mantle plume. This study contributes seismic constraints and CSRM‐1.0 to the CSRM product center (http://chinageorefmodel.org) as a backbone open‐access geophysical cyberinfrastructure.

Is the mantle transition zone uniform beneath Precambrian shields?

In the present study, we examine the mantle transition zone (MTZ) structure below the prominent Precambrian shields on the Earth to understand its thermal and compositional properties and depth distribution of MTZ discontinuities. For this study, we used data from Canadian, Brazilian, Baltic, African and Australian Shields. We migrate the receiver functions to the depth domain using 3D tomographic models to examine the topography of MTZ discontinuities that define the upper and lower boundaries of the MTZ. The depth migration results were obtained using two different 3D tomographic velocity models, LLNL_G3D_JPS and GyPSuM. The analysis shows that both models yield similar results. The findings indicate a thinner than usual MTZ beneath all the Precambrian shields with an average thickness of ∼238 ± 8 km. Notably, the present study reveals the upper boundary of the MTZ (410 km discontinuity) displays distinctive topography; in contrast, the lower boundary of the MTZ (660 km discontinuity) is situated at shallower depths. Therefore, the main factor causing the variation in MTZ thickness is the shallowing of the 660 km discontinuity, indicating that the post-spinel transition occurred at higher temperatures with a negative Clapeyron slope, supporting the theory of whole-mantle convection. This suggests that mantle plumes have some appreciable impact on the MTZ beneath Precambrian shields, which are limited to the base of the MTZ. Alternatively, the global mantle warming explanation could be invoked, as it best explains the similar characteristics of the 660 km discontinuity beneath the Precambrian shields. However, this phenomenon needs to be tested through numerical modelling.