Search

Subduction processes play a pivotal role in facilitating material exchange between the crust and mantle, contributing to the growth of continents. However, the onset and evolution of subduction remain hotly debated. Here, we developed a high‐dimensional machine learning (ML) model to use multiple compositional data (e.g., Nb/La, Nb, Ti, Nb/U, Pb/Nd, and Nb/Th) to distinguish arc‐type from non‐arc basalts worldwide, then applied this well‐trained ML model to identify and delineate the secular occurrence of Archean arc‐type basalts. Our findings indicate that basaltic arc magmatism can be traced back to only a few early Archean cratons. A noticeable increase in the prevalence of arc‐type basalts during the Mesoarchean and Neoarchean suggests the synchronous formation of Archean subduction across different cratons. This observation hints at transitioning from the predominantly stagnant geodynamic regime to the early, more mobile tectonic regime. Plain Language Summary: Deep beneath the Earth's surface, subduction moves material from the crust to the mantle, generating magma and igneous activity. This process is crucial for building the continental crust. Despite its importance, how the subduction process evolved over billions of years remains debated. To better understand when and where subduction began, we trained a machine learning model to distinguish between arc (subduction‐related) and non‐arc (not subduction‐related) basalts to recognize the chemical "fingerprints" of subduction in these rocks. By applying the same model to ancient basalts, we could trace basaltic arc‐like magmatism back to several Archean eon cratons, known as cratons, early in Earth's history (during the Archean Eon). The number of arc‐type basalts increased during the Mesoarchean and late Archean, suggesting that subduction expanded to different regions for about 700 million years. The transition suggests a significant shift in how the Earth's crust moved and changed, from a primarily steady state to a more mobile tectonics as observed today. Our study sheds new light on the ancient history of our planet and helps us understand how the Earth's crust has evolved over billions of years. Key Points: Machine learning can accurately distinguish between arc and non‐arc basalts worldwideSubduction‐related basaltic magmatism may have produced on several cratons during the early ArcheanSubduction‐related processes may have started in cratons worldwide during the late Archean [ABSTRACT FROM AUTHOR]

The mechanisms governing a commonly observed seismic velocity drop in the cratonic lithosphere, referred to as the mid‐lithospheric discontinuity (MLD), have been widely debated. To identify the composition and seismic structure of MLDs, we have analyzed Sp receiver functions (SRF) and mantle xenocrysts for six regions across Australia. We utilize locations where seismic stations and kimberlite‐hosted mantle xenocrysts are both available, allowing for comparison between seismological and petrological constraints. Our results show negative SRF phases indicative of the MLD coincide with clinopyroxene‐depleted zones at 60–140 km depth. Clinopyroxenes with different chemical compositions across the MLD define a litho‐chemical discontinuity. Modeling and experimental data show that MLDs may be explained by modified lherzolite with 10%–20% modal pargasite. Pargasite MLDs may form when rising H2O‐bearing melts cross the amphibole dehydration curve and react with clinopyroxene in lherzolite. Because the amphibole dehydration curve is isobaric at 80–120 km, pargasite will be precipitated as horizontal channels. Plain Language Summary: Seismic imaging of Earth's lithosphere has revealed a seismic velocity drop at 60–120 km depth, termed the mid‐lithosphere discontinuity (MLD). The mechanisms governing the MLD have been extensively debated. To advance the understanding of the MLD we carried out an interdisciplinary study to investigate the seismic structure and composition of MLDs identified in the cratonic lithosphere of Australia. We integrated analysis of seismic converted signals from discontinuities recorded by permanent seismometers with geochemical data for kimberlite‐hosted mantle xenocrysts collected in the vicinity. The method allows direct comparison between seismological and petrological constraints. Our results suggest that the MLD comprises anomalously low abundances of clinopyroxene and separates geochemically distinct layers within the lithospheric mantle. Experimental results and modeling suggests that the observed decrease in seismic velocity and absence of clinopyroxene may relate to the formation of pargasite channels in modified lherzolite. Key Points: Sp receiver functions and mantle xenocrysts are used to study mid‐lithosphere discontinuities beneath Australian cratonsMid‐lithosphere discontinuity at 60–120 km depth corresponds with mantle xenocryst populations that are depleted in clinopyroxeneModeling and experimental data suggests the mid‐lithosphere discontinuity is caused by pargasite‐bearing lherzolite [ABSTRACT FROM AUTHOR]

The evolution of the Gondwana along the flank of the West African Craton was complex and is far from understood. Subduction-related activity along this margin spanned between c. 750 and 500 Ma. Sections close to African cratons record the earliest stages, while Autochthonous and Allochthonous domains of the Variscan Belt preserve the latest stages of the arc system, essentially between c. 540 and 500 Ma. The geochemistry of the Ediacaran-early Cambrian siliciclastic series deposited along this Cadomian active margin preserves the evolutionary history of their sources, which are related to activity in the arc and nearby continental areas. In this sense, the SW Iberian Massif (Ossa-Morena Complex) preserves a section of this Ediacaran-early Cambrian peri-Gondwanan arc. Its evolution can be tracked through the characterization of the subduction-related magmatism (including the Mérida Massif) and coeval metasedimentary record (Serie Negra Group and Malcocinado Formation) during a time interval spanning almost 100 m.y., from pre-602 Ma to at least c. 534 Ma. This study reveals that arc magmatism is closely linked with synorogenic deposition in a complex way so far unexplored. Arc recycling is revealed by the isotopic equivalence of synorogenic strata to the first magmatic event (pre-602 Ma), and by geochronological data of the arc-building pulses. The earliest magmatic pulses (c. 602–550 Ma) are characterized by significant crustal input, likely favoured by subduction erosion. Subsequently, magmatism evolved towards larger mantle involvement (c. 540–534 Ma), likely associated with progressive variation in the slab angle. These slab-mantle-upper plate interactions generated changes in the arc dynamics leading to an extensional setting with alkaline magmatism during the Cambrian. This review proposes a model of petrogenetic and geodynamic arc evolution between the Ediacaran and the Early Cambrian. The gathered data could improve the accuracy of future palaeogeographic reconstructions for the northern margin of Gondwana. [ABSTRACT FROM AUTHOR]

Syncollisional magmatism plays an important but underappreciated role in continental crust growth and maturation. However, the origin of syncollisional intermediate magmas in continental subduction zones is controversial, with some models suggesting they form by arc-related processes, and others indicating they form by later slab breakoff-induced melting. Diorite porphyry dikes intruding granitic gneiss in the Paleo-Tethyan Sulu ultrahigh-pressure (UHP) continental collisional orogen have inherited zircon grains with 206Pb/238U ages of ca. 749-238 Ma, and magmatic zircons with weighted mean ages of 216-215 Ma, falling within the well-constrained time range (ca. 235-208 Ma) tracking exhumation of the Sulu UHP rocks from UHP peak conditions to amphibolite facies; they are thus syncollisional. The dikes have high Cr (330-402 ppm), Ni (84.5-103 ppm), and Mg# (64-66) values, showing a mantle origin. The porphyries have relatively high Sm/Yb, Nb/Y, La/Yb, and Gd/Yb ratios, representing a classic signature of slab breakoff magmatism. Together with the arc-like traceelement patterns and enriched Sr-Nd isotope compositions, ages, and εHf(t) values (-19.5 to -17.0) of magmatic zircons and their tectonic setting, we propose a syncollisional slab breakoff model in which the melts were initially generated from asthenospheric upwelling in the gap created when the oceanic slab attached to the Yangtze craton detached underneath the North China craton during Late Triassic collision following Paleo-Tethys Ocean closure. The diorite porphyry dikes have consistent Sr-Nd isotope compositions and spatiotemporal relationships with the nearby Shidao gabbro-syenite-granite complex, for which the tectonic affinity is controversial. Thus, we argue that the diorite porphyries and Shidao complex were sourced from two cratons, including the enriched subcontinental lithospheric mantle of the North China craton, which interacted with abundant felsic melts derived from the sinking slab breaking away from the subducted crust of the Yangtze continental-ocean transitional margin. This study sheds new light on crustal recycling versus continental growth in collisional orogens and implies that considerable syncollisional intermediate magmas can be generated by slab breakoff in continental subduction zones, representing hybrid additions to continental growth that are different and more evolved than arc magmas and have a composition similar to that of the bulk continental crust. [ABSTRACT FROM AUTHOR]

There are numerous Neoarchean metasedimentary rock packages in the Slave craton, and they serve as important archives of tectonic processes. Little is documented on the Neoarchean sedimentary packages of the Winter Lake greenstone belt of the central Slave craton, however, and their interpretation can aid in the understanding of the final stages of Slave craton amalgamation. This project investigates the depositional environments and tectonic settings of the Itchen Formation and Sherpa Formation of the Winter Lake greenstone belt. Our study provides constraints for reconstructing the Neoarchean evolution of the central Slave craton through bedrock mapping and facies analysis. The Itchen Formation consists of submature mudstone, siltstone, and sandstone, with preserved graded bedding, planar bedding, and flame structures. Unconformably overlying the Itchen Formation is the Sherpa Formation, which is dominated by polymictic conglomerates and coarse-grained sandstones with preserved cross-bedding, imbricated clasts, and scour surfaces. The Itchen Formation is interpreted to have been deposited in a convergent basin (i.e., retro-arc foreland basin), where two facies associations outline turbidite and suspension sedimentation consistent with submarine fan deposition on a continental slope and a basin floor environment. By contrast, the Sherpa Formation has three facies associations representing dominantly alluvial–fluvial environments in terrestrial–marine–lacustrine settings deposited in pull-apart basins resulting from transtensional forces associated with the Beniah fault zone. [ABSTRACT FROM AUTHOR]

Comprehending the temperature distribution within Earth's lithospheric mantle is of paramount importance for understanding the dynamics of Earth's interior. Traditional mineral‐based thermobarometers effectively constrain temperature and pressure for particular compositions, but their application is limited at the global scale. Here, we trained machine‐learning (ML) algorithms on 985 published high‐temperature and high‐pressure experiments for use as thermometers and barometers to overcome the limitations of classic methods. We compared our ML models to classic thermobarometers to assess the accuracies of predicted pressures and temperatures. The comparison shows that the ML models outperform classic methods and better fit various mineral pairs. Global application of the ML models unveils mantle conditions beneath cratons. Furthermore, depths to the lithosphere‐asthenosphere boundary (LAB) calibrated based on the ML thermobarometry results are generally deeper by ∼40 km than those derived geophysically, implying the existence of melt‐bearing or hydrated mineral zones at the LAB. Plain Language Summary: The temperature of the lithospheric mantle is pivotal to processes such as the production of magmas and the structural stability of cratons. Researchers have traditionally employed mineral‐based thermobarometers to lithospheric mantle temperatures, but those conventional techniques have limited applications and are susceptible to inaccuracies. Here, we explored a new approach using machine learning, a powerful tool for identifying complex patterns in high‐dimensional data space. We trained machine‐learning models using data from high‐pressure and high‐temperature experiments, then compared our machine‐learning models to traditional methods. Our results show that the machine‐learning models better predict pressure and temperature for specific mineral combinations. We also used our optimized model to predict mantle temperatures and pressures based on a global data set of xenolith analyses, revealing insights about the temperature of the lithospheric mantle and the depth to the lithosphere‐asthenosphere boundary beneath continents. This research shows that machine learning can greatly improve our understanding of Earth's deep processes, providing more accurate insights into its dynamics and evolution. Key Points: Thermometers and barometers calibrated and validated for xenolithic mantle minerals by machine learningGlobal lithospheric mantle temperatures predicted based on xenolithic mineral compositions calculated by machine‐learning thermobarometryThermal continental lithosphere‐asthenosphere boundarys refined by machine‐learning thermobarometry are consistently ∼40 km deeper than those determined geophysically [ABSTRACT FROM AUTHOR]

The study of the volcano–sedimentary sequence in the lower part of the Zaonega formation in the Paleoproterozoic Onega structure (Karelian craton, Fennoscandian Shield) has shown that tuffaceous and high-silica rocks predominate in its composition. High-silica rocks (SiO2 up to 94 wt %) are depleted in all elements and probably represent chemogenic siliceous silts. Tuffaceous rocks vary broadly in composition and were formed during mixing of tuffogenous basalt material and high-silica chemogenic sediments. In terms of levels of contents and the distribution character of trace elements, tuffogenous rocks are close to N-MORB volcanic rocks. This rock association is typical of the early stages of continental rifting in the Phanerozoic and may indicate the formation of volcano–sedimentary complexes of the Zaonega formation in the environment of continental rifting. The tuffaceous rocks in the lower part of the Zaonega formation are geochemically identical to dolerite dikes and MORB-type volcanic rocks 2.10–2.14 Ga in age. Their formation was probably related to this episode of large-scale extension and thinning of the continental lithosphere of the Karelian craton in the Middle Paleoproterozoic. In this case, the age limit of the Zaonega and underlying Tulomozero formations should be somewhat older than the 2.06–2.10 Ga interval accepted in the modern regional stratigraphic schemes of the Paleoproterozoic Fennoscandian shield. [ABSTRACT FROM AUTHOR]

Surface topography is an important yet largely neglected aspect of the early evolution of cratons. The lateral accretion of cratonic nuclei inevitably forms orogenic belts that subsequently provide a sediment source for large, resource‐rich intracratonic basins, but to date, geodynamic models have focused exclusively on lithospheric root processes. Here we use two‐dimensional thermal‐mechanical models to study the topography and lithospheric deformation during 50 Myr of compression of a cratonic nucleus, to simulate the lateral accretion phase of craton growth in the Neoarchean. Although the cratonic nucleus thickens slightly during the compression phase, most of the deformation occurs in the regions adjacent to the nucleus that have weaker lithosphere. Here, crustal thickness triples developing high topography in excess of 10 km without active erosion. Models with different initial rheological parameters will have different final topography and lithosphere geometry, but in general it is difficult to shorten and deform the depleted cratonic nucleus, unless there are significantly weak heterogeneities in the mantle lithosphere. We apply two quantitative analysis techniques to objectively evaluate a multitude of model outputs. Cross‐correlation clustering (CCC) measures the degree of similarity between topography profiles and categorizes models based on the general topographic character. Six different topography families are possible in the context of our models and crustal strength is the most important parameter affecting the shape. From principal component analysis (PCA) we identify four dominant lithosphere geometries. When used together, these two methods provide distinct yet complementary information about the surface and subsurface deformation features in our models. Plain Language Summary: Cratons, the ancient building blocks of continents, have existed on our planet, virtually unchanged, for the last 2 to 3 billion years. These regions contain the world's largest deposits of economic and precious metals, and they preserve the earliest evidence of life on Earth. The longevity of cratons is attributed to their unique chemical composition and structure that modern day plate tectonics is unable to reproduce. Therefore, understanding how cratons form may help unlock some of the mystery of Earth's early history. While many previous geodynamic studies on this topic have been published, the evolution of surface topography during craton formation via lateral compression has been ignored until now. We present 54 numerical models to test a wide range of initial conditions, showing that topographic relief is largely controlled by crustal strength. We also demonstrate how quantitative analysis techniques can effectively cluster models based on their large‐scale behavior and therefore aid in the interpretation of model results. Key Points: Quantitative clustering techniques effectively summarize results from numerous numerical models and facilitate objective comparisonSignificant surface topography formed during cratonic nuclei compression in supercontinent cycles is strongly controlled by crustal strengthHomogenous cratonic lithosphere is resistant to pure shear thickening when subjected to horizontal compression [ABSTRACT FROM AUTHOR]

Metamorphosed dacitic porphyry dikes were first found in the western part of the Vorontsovka terrane, which is located in the Paleoproterozoic Volga–Don orogen at the margin of Archean Sarmatia and Volga–Ural cratons. The magmatic protolith age for the metadacites is ca. 2.07 Ga. These are ferrous, metaluminous calc-alkali I-type granitoids. The sodium specialization of the rocks and their low concentrations of Mg, Cr, Ni, and incompatible elements, with significant REE fractionation, the absence of Eu* anomalies, high Sr/Y ratio, remarkably high (Gd/Yb)n values (>10), and the radiogenic Nd isotopic composition indicate that the dacitic melts were derived from a juvenile mafic source. According to petrogenetic estimations, such conditions could be caused by the partial melting of depleted N-MORB basites in equilibrium with an eclogitic residue. The dacitic magmas were likely generated by the partial melting of mafic rocks at lower levels of the significantly thickened crust (>60 km) in relation to collision processes. [ABSTRACT FROM AUTHOR]

The thermal structure of the lithosphere is key to understanding its thickness, properties, evolution, and geothermal resources. Cratons are known for their low heat flow and deep lithospheric roots. However, present-day cratons in East China have geothermal characteristics that are highly complex, with variable heat flow values, diverging from the typical thermal state of cratons. In this study, we conducted a detailed analysis of the geothermal geological background of the cratons in East China, summarizing the thermal state and tectono-thermal processes of different tectonic units, calculating the temperature at various depths, and discussing differences in temperature and thermal reservoirs at different depths. The observed lithospheric thermal thickness within the North Jiangsu Basin and the Bohai Bay Basin is notably reduced in comparison to that of the Jianghan Basin and the Southern North China Basin. The phenomenon of craton destruction during the Late Mesozoic emerges as a pivotal determinant, enhancing the geothermal resource prospects of both the Bohai Bay Basin and the North Jiangsu Basin. Our findings contribute significantly to the augmentation of theoretical frameworks concerning the origins of heat sources in global cratons. Furthermore, they offer invaluable insights for the methodical exploration, evaluation, advancement, and exploitation of geothermal resources. [ABSTRACT FROM AUTHOR]

The gravity anomalies reflect density perturbations at different depths, which control the physical state and dynamics of the lithosphere and sub-lithospheric mantle. However, the gravity effect of the crust masks the mantle signals. In this study, we develop two frameworks (correction with density contrasts and actual densities) to calculate the gravity anomalies generated by the layered crust. We apply the proposed approaches to evaluate the global mantle gravity disturbances based on the new crustal models. Consistent patterns and an increasing linear trend of the mantle gravity disturbances with lithospheric thickness and Vs velocities at 150 km depth are obtained. Our results indicate denser lithospheric roots in most cratons and lighter materials in the oceanic mantle. Furthermore, our gravity map corresponds well to regional geological features, providing new insights into mantle structure and dynamics. Specifically, (1) reduced anomalies associated with the Superior and Rae cratons indicate more depleted roots compared with other cratons of North America. (2) Negative anomalies along the Cordillera (western North America) suggest mass deficits owing to the buoyant hot mantle. (3) Positive anomalies in the Baltic, East European, and Siberian cratons support thick, dense lithosphere with significant density heterogeneities, which could result from thermo-chemical modifications of the cratonic roots. (4) Pronounced positive anomalies correspond to stable blocks, e.g., Arabian Platform, Indian Craton, and Tarim basin, indicating a thick, dense lithosphere. (5) Low anomalies in the active tectonic units and back-arc basins suggest local mantle upwellings. (6) The cold subducting/detached plates may result in the high anomalies observed in the Zagros and Tibet. [ABSTRACT FROM AUTHOR]

The Bloubergstrand Member, consisting of mafic to intermediate tuffs and flows within the Tygerberg Formation of the Malmesbury Group, is a crucial but understudied component of the western Saldania Belt in southwestern Africa. With a U-Pb zircon age of 555 ± 5 Ma, these volcanic rocks provide valuable insights into the origins and, more broadly, the tectonic setting of the Saldania Belt during the final stages of Malmesbury Group sedimentation in the context of the construction of southwestern Gondwana at approximately 560 to 540 Ma. The Bloubergstrand Member amygdaloidal volcanic rocks vary from basaltic andesitic to andesitic and have experienced varying degrees of alteration. The volcanic rocks are enriched in the large ion lithophile elements (LILE) relative to the high field strength elements (HFSE). When normalised to primitive mantle values, they show enrichment in Cs, Rb, K and Pb and mild enrichment in Th, U, Zr and Hf, with the light rare earth elements (LREE) enriched relative to the heavy rare earth elements (HREE). They are depleted in Ba, Nb, Ta, Sr, Eu, P and Ti. eNd(t) values are mildly negative ranging from -3.60 to -2.39, and Nd TDM model ages range from 1.4 to 1.7 Ga. Initial 87Sr/86Sr ratios vary from 0.70478 to 0.70620, with one higher value of 0.72270, the latter likely due to extensive alteration. The Bloubergstrand Member volcanic rocks exhibit characteristics suggesting their origin from partial melting of lithospheric to transitional asthenospheric upper mantle, influenced by sediment-derived melts and variable degrees of crustal contamination. With continental arc basalt compositions and similarities to shoshonites, the samples reflect variable degrees of partial melting and source heterogeneities. Fractional crystallisation of pyroxenes, hornblende, and plagioclase contributes to compositional variability. Erupting in a continental arc margin, likely part of the Arachania block of the Kalahari craton, separated during the Tonian break-up of Rodinia, the volcanics were associated with the Marmora basin formed by eastward-directed subduction below the western Kaapvaal Craton margin. Extruded in a distal position relative to the Cuchilla Dionisio Arc in present-day southern Brazil and Uruguay, these volcanic rocks preceded the closure of the southern Marmora basin. The mantle melting, possibly a result of slab roll-back, break-off, ridge subduction, or a combination, served as a precursor to the lower crustal heating that generated the granitic magmas of the Cape Granite Suite. [ABSTRACT FROM AUTHOR]

The start of the Paleoproterozoic supercontinent cycle is typically taken as the initiation of orogenesis at ca. 2.1 Ga leading to the assembly of Earth's first supercontinent, Columbia. However, the dearth of ca. 2.5–2.2 Ga geological records makes it difficult to deduce tectonic factors during the onset of the Paleoproterozoic supercontinent cycle. The petrogenesis of tonalite–trondhjemite–granodiorite (TTG) provides useful proxies for tracing prevailing geodynamic regimes of early continental evolution. However, marked decreases of TTG and other magmatism occurred across the Archean–Paleoproterozoic transition and have previously precluded forming testable hypotheses. Early Paleoproterozoic TTGs have been identified in the North China Craton (NCC) and other cratons, which may represent the last major pulse of TTGs globally. Here we present low δ18O and δ30Si ca. 2.3 Ga TTGs from the NCC, together with thermodynamic modeling and compilation of stable O and Si isotopes for TTGs globally through time. The ca. 2.3 Ga TTGs were derived from the partial melting of Archean basaltic crust and give lighter average zircon δ18O (3.15 ± 0.35‰) and whole‐rock δ30Si values (−0.17 ± 0.08‰) than most Archean TTGs. Considering coeval mafic‐felsic igneous rocks, and lithospheric thinning since ca. 2.5 Ga based on estimated crustal thickness through the Neoarchean–Paleoproterozoic, we posit the onset of an intraplate rifting consistent with the anomalous low‐δ18O magmatism. Continental rifting of Archean cratons/supercratons plausibly hints at the formation of rifts driving subduction initiation as the veritable onset of the Paleoproterozoic supercontinent cycle. Plain Language Summary: The Earth's first supercontinent, Paleoproterozoic–Mesoproterozoic Columbia, likely began to form during ca. 2.1 Ga (billion years ago) orogenesis. However, the limited record of ca. 2.5–2.2 Ga magmatism makes it difficult to establish how exactly the Paleoproterozoic supercontinent cycle commenced. The tonalite–trondhjemite–granodiorite (TTG) series of felsic rocks forms the backbone of Archean cratons and can provide key constraints on the Earth's early tectonic evolution. Some early Paleoproterozoic TTGs occur in North China and other cratons, which may represent the last major TTG event globally, presenting the possibility to fulfill this role. Here we carry out zircon O and bulk‐rock Si isotopic analyses and thermodynamic modeling for these ca. 2.3 Ga TTGs, together with Si‐O‐Hf data compilation of TTGs globally through time. These ca. 2.3 Ga TTGs show lighter zircon O and bulk‐rock Si isotopes than those of Archean TTGs. They might be formed within‐plate, associated with lithospheric thinning and mantle‐crust remelting to form mafic‐felsic magmatism. We propose an intraplate rift setting for this early ca. 2.3 Ga magmatism, which drove the mobilism that led to the subsequent 2.1–1.8 Ga orogenesis and should be posited as the real start of the Paleoproterozoic supercontinent cycle. Key Points: Ca. 2.3 Ga tonalite–trondhjemite–granodiorite (TTG) in the North China Craton (NCC) shows the affinity of low‐pressure (LP) TTGs and contains low zircon δ18O and bulk‐rock δ30Si valuesCa. 2.3 Ga low δ18O and δ30Si TTG and coeval low‐δ18O magmatism herald the initiation of an intraplate rift settingSuch intraplate rifting hints at the formation of rifts driving subduction as the real onset of the Paleoproterozoic supercontinent cycle [ABSTRACT FROM AUTHOR]

The article explores the concept of eclogite viscous drainage and its impact on the evolution of the continental lithosphere. Eclogites, which are high-pressure metamorphic rocks, can flow within the lithospheric mantle through deformation. This process has been observed in the Pamir Plateau and may have occurred in other accretion zones throughout Earth's history. The article discusses various mechanisms associated with eclogite, such as delamination, dripping, and subduction acceleration, and their effects on the continental lithosphere's evolution. It also highlights the role of eclogite viscous drainage in stabilizing and cratonizing the lithosphere. [Extracted from the article]

We present high-resolution 3-D images of isotropic P -wave velocity (Vp), azimuthal anisotropy (AAN) and radial anisotropy (RAN) down to 700 km depth beneath the North China Craton (NCC) and adjacent areas, which are obtained by inverting a great number of high-quality arrival time data recorded at 1374 portable seismic stations and 635 permanent stations in the study region. Our results reveal new and detailed features of the upper mantle structure beneath the NCC. Varying structural heterogeneities are revealed beneath different tectonic blocks, and differences also exist between northern and southern parts of each block. The fast velocity directions (FVDs) of azimuthal anisotropy are mainly NW–SE under the Alaxa block, and NE–SW beneath the Tibetan Plateau. The FVDs present an arc transition along the boundary faults separating the Tibetan Plateau, the Alaxa block, the western NCC, and the Sichuan basin. Low- Vp anomalies with positive RANs (i.e. horizontal Vp  > vertical Vp) are revealed at 100–200 km depths under the Tibetan Plateau, reflecting frozen-in anisotropy in the thick lithosphere. Significant low- Vp anomalies with a circular AAN pattern exist at 0–700 km depths beneath the Datong volcano. In addition, negative RAN occurs right below the volcano, whereas positive RANs appear around it, suggesting that the Datong volcano is fed by hot upwelling flow from the lower mantle associated with collapsing of subducted slab materials down to the lower mantle. The eastern NCC shows complex Vp AANs and RANs. Seismic anisotropy exhibits east–west variations in the upper mantle across the Tanlu fault zone. The west of the Tanlu fault shows negative RANs (vertical Vp  > horizontal Vp), whereas its east shows positive RANs at 300–500 km depths. The low- Vp anomaly under the Datong volcano is connected with a large low- Vp anomaly beneath the eastern NCC above ∼250 km depth, suggesting that the hot upwelling flow under Datong may migrate laterally to the asthenosphere under the eastern NCC and contribute to the lithospheric delamination and destruction there. [ABSTRACT FROM AUTHOR]

Inconsistencies between observations from long and short period seismic waves and geochemical data mean craton formation and evolution remains enigmatic. Specifically, internal layering and radial anisotropy are poorly constrained. Here, we show that these inconsistencies can be reconciled by inverting cratonic Rayleigh and Love surface wave dispersion curves for shear‐wave velocity and radial anisotropy using a flexible Bayesian scheme. This approach requires no explicit vertical smoothing and only adds anisotropy to layers where required by the data. We show that all cratonic lithospheres are comprised of a positively radially anisotropic upper layer, best explained by Archean underplating, and an isotropic layer beneath, indicative of two‐stage formation. Within the positively radially anisotropic upper layer, we find a variable amplitude low velocity zone within 9 of 12 cratons studied, that is well correlated with observed Mid‐Lithospheric Discontinuities (MLDs). The MLD is best explained by metasomatism after craton formation. Plain Language Summary: The ancient cores of the continents, or cratons, are a treasure‐trove of >2.5 billion years of Earth's history. However, scientists disagree on the processes that led to their formation, or whether they have evolved significantly through time. This is because the geological and geophysical methods used to investigate cratons often yield conflicting results. By capitalizing on an up‐to‐date global long‐wavelength seismic data set, we image the cores of 12 cratons using an advanced statistical method, Bayesian inference. The flexible method requires few choices to be made a priori, is driven by the quality of the data itself and measures uncertainties on results. By detecting velocity differences between horizontally and vertically vibrating seismic waves, we show that all cratons likely comprise an upper layer formed during the hot early Earth, by a process that strongly aligns the constituent minerals of the rocks in the horizontal plane. Below this ∼125 km thick upper layer, a lower layer (∼150 km thick) shows no clear alignment of minerals and so was likely formed by a different process, at a later time. Variable slow wavespeed zones exist within the upper layer that match previous results from short‐wavelength seismic data. Key Points: Bayesian surface wave inversion reconciles existing differences between seismic images of craton structure from long and short period dataCratons show a shallow low velocity zone (LVZ) within a layer of positive radial anisotropy and a high velocity isotropic layer beneathCratons are formed in two stages shown by anisotropic structure and are later modified producing a LVZ and mid lithosphere discontinuities [ABSTRACT FROM AUTHOR]

We present the first magnetotelluric sounding (MTS) data acquired in the southern part of the Lena River delta which is a junction zone of the Siberian Craton and the Verkhoyansk fold-and-thrust belt, in the transition zone from the Eurasian Continent to the shelf of the Laptev Sea. The MTS data were used to construct a vertical section of the bulk electrical resistivity (ER) structure of the Earth's crust down to a depth of 8 km. The upper high-resistivity layer (320–1000 Ohm m) corresponds to the permafrost rocks ranging in thickness from 1 km to 400 m under the channels of the Lena River. The underlying deformed Mesoproterozoic–Lower Triassic rocks can vary in composition, but they are poorly differentiated in ER values (110–240 Ohm m). We detected three low-resistivity anomalies (10–60 Ohm m) related to the fluid-saturated core and damage zones of active faults. A deep-lying high-resistivity anomaly (320–350 Ohm m) identified at the northeastern part of the section can correspond to the Lower Proterozoic or Archean metamorphic rocks. [ABSTRACT FROM AUTHOR]

This study presents results from apatite fission track (AFT) thermochronology to investigate the thermal history and exhumation dynamics of the Rio Negro–Juruena basement, situated within the western Guiana Shield of the Amazonian Craton. AFT dating and associated thermal history modeling in South America has largely been restricted to the plate's margins (e.g., Andean active margin, Brazilian passive margin, and others). Our paper reports on low-temperature thermochronological data from the internal part of the western Guiana Shield for the first time. This area is part of a vast cratonic lithosphere that is generally thought to be stable and little influenced by Mesozoic and Cenozoic tectonics. Our data, however, show AFT central ages ranging from 79.1 ± 3.2 to 177.1 ± 14.8 Ma, with mean confined track lengths of ca. 12 µm. Contrary to what might be expected of stable cratonic shields, inverse thermal history modeling indicates a rapid basement cooling event in the early Cretaceous. This cooling is interpreted as a significant exhumation event of the basement that was likely driven by the coeval extensional tectonics associated with back-arc rifts in the Llanos and Putumayo–Oriente–Maranon basins. The extensional tectonics facilitated both basement uplift and subsidence of the adjoining basins, increasing erosional dynamics and consequent exhumation of the basement rocks. The tectonic setting shifted in the late Cretaceous from extensional to contractional, resulting in reduced subsidence of the basins and consequential diminishing cooling rates of the Guiana Shield basement. Throughout the Cenozoic, only gradual, slow subsidence occurred in the study area due to regional flexure linked to the Andean orogeny. Comparative analysis with low-temperature thermochronology data from other west Gondwana cratonic segments highlights that exhumation episodes are highly controlled by tectonic inheritance, lithospheric strength, and proximity to rift zones. This study underscores the complex interplay between tectonic events and the response of cratonic lithosphere over geological timescales and highlights extensional settings as an important geological context for craton exhumation. [ABSTRACT FROM AUTHOR]

Seismic interrogation of the upper mantle from the base of the crust to the top of the mantle transition zone has revealed discontinuities that are variable in space, depth, lateral extent, amplitude and lack a unified explanation for their origin. Improved constraints on the detectability and properties of mantle discontinuities can be obtained with P-to-S receiver function (Ps-RF) where energy scatters from P to S as seismic waves propagate across discontinuities of interest. However, due to the interference of crustal multiples, uppermost mantle discontinuities are more commonly imaged with lower resolution S-to-P receiver function (Sp-RF). In this study, a new method called CRISP-RF (Clean Receiver-function Imaging using SParse Radon Filters) is proposed, which incorporates ideas from compressive sensing and model-based image reconstruction. The central idea involves applying a sparse Radon transform to effectively decompose the Ps-RF into its underlying wavefield contributions, that is direct conversions, multiples, and noise, based on the phase moveout and coherence. A masking filter is then designed and applied to create a multiple-free and denoised Ps-RF. We demonstrate, using synthetic experiment, that our implementation of the Radon transform using a sparsity-promoting regularization outperforms the conventional least-squares methods and can effectively isolate direct Ps conversions. We further apply the CRISP-RF workflow on real data, including single station data on cratons, common-conversion-point stack at continental margins and seismic data from ocean islands. The application of CRISP-RF to global data sets will advance our understanding of the enigmatic origins of the upper mantle discontinuities like the ubiquitous mid-lithospheric discontinuity and the elusive X-discontinuity. [ABSTRACT FROM AUTHOR]

Distributions of long-lived radioelements (Th, U and K) in the upper crust play a vital role in lithospheric thermal modelling and understanding of underlying geodynamic processes. In the present study, we report abundances of Th, U, K and radiogenic heat production (A) in Paleoarchaean and early Neoarchaean granitoids from the Singhbhum Craton, located in the eastern part of the Indian shield, for the first time in a systematic way, using laboratory gamma-ray spectrometric set-up. We have studied 204 samples consisting of Paleoarchaean gneiss, three phases of Paleoarchaean Singhbhum Granite and early Neoarchaean granitoid. Th, U, K and A in the Paleoarchaean Singhbhum Granite are found to be marginally higher (10.8 ppm, 1.4 ppm, 2.4 per cent and 1.3 μWm−3) than the Paleoarchaean gneiss (9.6 ppm, 1.7 ppm, 1.8 per cent and 1.3 μWm−3). In comparison, such values are much higher in the early Neoarchaean granitoid (31.1 ppm, 4.2 ppm, 3.9 per cent, and 3.6 μWm−3). The surface heat production (1.36 μWm−3) of this craton is comparatively low compared to most cratons worldwide. In addition, it has the lowest average crustal heat production (0.42 µWm−3), due to which it will contribute little to the observed surface heat flow. Moreover, Th is a dominant component in heat production, followed by U and K in these granitoids. Spatially, radioelement ratios Th/U and K/U show similar trends as the radioelements, indicating systematic depletion of Th and U, compared to K. The study also indicates that the spatially separated coeval granitoids are formed from different magma sources in the heterogeneous crust that prevailed during the Archaeans. The lower radioelement abundances in the Paleoarchaean granitoids mostly resulted from the partial melting of mafic sources, whereas higher radioelement abundances in the early Neoarchaean granitoids resulted from the partial melting of the felsic source. [ABSTRACT FROM AUTHOR]

Ancient orogens within the supercontinent like Columbia can remain stable evolution as long as the cratons. What kind of lithospheric mantle was beneath those orogens and how it evolved into a stable state are still enigmatic. The Trans‐North China orogen (TNCO) is one of the typical collisional orogens within the Columbia supercontinent and was formed at ca. 1.85 Ga. Our work reveals that a cluster of kimberlites intruded the orogenic belt at ca. 1.54 Ga. These rocks were originally generated under a thick lithosphere (>200 km). Their entrained olivine cores show a composition of overlapping olivines from refractory mantle peridotites. The results suggest a thick and refractory lithospheric mantle beneath the TNCO at ca. 1.54 Ga. Such craton‐like property may result from large volume melt extraction from the lithospheric mantle, possibly caused by the ca. 1.78 Ga large igneous event, which eventually induces the long‐term stability of the TNCO during the subsequent supercontinent cycle. [ABSTRACT FROM AUTHOR]

The Cenozoic sediments are very thick in the southwest Tarim Basin and very thin in the northwest, but what controls these variations is unclear. Here, we use two‐dimensional thermo‐mechanical models to investigate how the lateral variations in rheological strength and depletion density of cratonic lithosphere mantle affect the cratonic basin deformation. Model results show that the basin basement uplift occurs above either the region with crustal thickening or high depletion in the mantle. A model with a stronger and density‐depleted northern half of cratonic lithosphere mantle in the context of compression matches the differential Cenozoic subsidence and deformation observed in the Tarim Basin well. We propose that a Permian plume led to the lateral heterogeneity of the lithosphere mantle under the Tarim craton, and the modified lithosphere mantle characteristics caused the differential Cenozoic sediment thickness in the Tarim Basin. Plain Language Summary: Tarim Basin is a typical intraplate cratonic basin in northwest China. It is bounded by the West Kunlun Shan to the south and the Tian Shan to the north, respectively. The Cenozoic sediment thickness in the Tarim Basin displays significant lateral variation, which cannot be explained by flexural bending of the Tarim lithosphere under orogenic loads of the West Kunlun Shan and Tian Shan. Recent geophysical data outlined the extent of the Permian plume head, which coincides with the region covered by thinnest Cenozoic sediments. Here, we conduct a series of 2‐D thermo‐mechanical numerical simulations to examine the correlation of the cratonic lithosphere mantle characteristics with cratonic basin's basement geometry. Model results show that a craton with a stronger and density‐depleted northern half of lithosphere mantle in the context of compression matches well with the south‐north differential Cenozoic subsidence and deformation in the Tarim Basin. Accordingly, we suggest that the differential Cenozoic sediment thickness in the Tarim Basin was likely caused by the Permian plume‐modified lithosphere mantle. Our study indicates that the plume‐modified mantle characteristics play a critical role in the basement deformation of the above sediment basin. Key Points: The strength and depletion density of cratonic lithosphere mantle control the deformation of cratonic lithosphere underlying sediment basinThe differential Cenozoic sediment thickness in the Tarim Basin is related to the heterogeneity of cratonic lithosphere mantleThe lithosphere mantle heterogeneity under the Tarim craton was likely caused by post‐Permian plume‐driven recratonization [ABSTRACT FROM AUTHOR]

The continental mid‐lithosphere discontinuity (MLD) is widely detected within cratons, with the dominant depth range of 70–100 km and a significant reduction of shear‐wave velocity of 2%–12%. However, the formation mechanism and corresponding strength of the MLD are widely debated, which may strongly affect the roles of MLD in craton evolution. The comparisons among variable mechanisms indicate that the strength of the MLD varies from the relatively high viscosity of wet olivine to the rather low viscosity of antigorite. Thus, systematic numerical modeling has been conducted with the MLD of contrasting strengths, that is, the wet olivine‐induced MLD or antigorite‐induced MLD, to investigate the roles of MLD in the craton instability under variable tectonic regimes (stable, extension, compression, mantle flow traction, or mantle plume). The models show that the cratonic lithosphere with wet olivine‐induced MLD maintains its stability under all the tectonic regimes. In contrast, the antigorite‐induced MLD with lowest viscosity could significantly promote the decoupling of lithosphere, and facilitate the lithospheric deformation. However, lithospheric delamination only occurs with the rather weak MLD interacting with the sub‐plate asthenosphere upwelling during craton extension or mantle plume activity. The sufficient amount of melts is essential for this process, which requires a large amount of extension or a mantle plume with rather high temperature anomaly and large size. Therefore, craton destruction is still difficult, and requires additional strict conditions. This may explain the general stability of most cratons with widespread MLDs. Plain Language Summary: Geophysical observations have revealed a seismic discontinuity within cratons, which is defined as continental mid‐lithosphere discontinuity (MLD). The MLD is often considered to be rather weak, which could facilitate the lithospheric delamination, and even the craton destruction. However, most cratons preserve the stability during long‐term evolution. Thus, the viscosity of the MLD and its roles in craton evolution require further investigations. In this study, we have conducted systematic numerical modeling to investigate the effects of the MLD on craton instability under variable tectonic regimes. The wet olivine‐induced MLD and antigorite‐induced MLD are two end‐members of MLD with contrasting strengths, which are compared in the present models. The model results indicate that only the connection and interaction between the weak antigorite‐induced MLD and the sub‐plate asthenosphere could cause craton destruction during tectonic extension or upwelling of mantle plume. The sufficient amount of melts is essential for this process to destroy the lower part of lithosphere. Thus, craton destruction requires additional strict conditions, which may explain the stability of most cratons on the Earth, although with the presence of MLD. Key Points: The viscosity of mid‐lithosphere discontinuity varies from relatively high viscosity of wet olivine to the rather low value of antigoriteThe wet olivine‐induced mid‐lithosphere discontinuity could not lead to craton destruction in any tectonic regimesCraton destruction occurs with the presence of antigorite‐induced mid‐lithosphere discontinuity combing additional strict conditions [ABSTRACT FROM AUTHOR]

Correlations between the Neoproterozoic belts surrounding the West African Craton and northern Brazilian cratons have long been a subject of interest and controversies. Due to the splitting of African and South American continents by the Atlantic oceanic domains, no direct links are preserved, requiring relying on various geological or geophysical characteristics to propose such correlations. In addition to the opening of the Atlantic oceanic domains, another difficulty arises from the covering of northern Brazilian belts by upper Paleozoic and Mesozoic basins, making these correlations speculative. Here, we propose new correlations based on the comparison between African and Brazilian belts. Recently, new geochronological and geodynamical data obtained in West Africa have evidenced four orogens in the Neoproterozoic belts of the western margin of the West Africa Craton, while the belts on the eastern side underwent only one orogeny. Similarities with the Pan-African I (900–650 Ma) and with the Pan-African II (650–480 Ma) orogenic events have been evidenced in the western Brazilian belts (Araguay and Paraguay). The first two orogens on the western margin (Pan-African I and Pan-African II) can thus be extended to the western Brazilian belts and can be considered as parts of a single geodynamic system running from the Mauritania to the Paraguay including the "Gurupi rift" as an aulacogen connected to the NNW-SSE Panafrican I and II oceanic domains. Consequently, the eastern Brazilian belt should rather be linked the Eastern Trans-Saharan belts. [ABSTRACT FROM AUTHOR]

Modern ideas about the tectonics of the Eurasian continental margin are considered. Cratons, fold belts, platforms, and sedimentary basins are briefly characterized. A significant part of the Arctic continental margins is represented by folded structures of different ages. In the course of geological evolution, tectonic structures successively formed, modified and died off, on which the processes of rifting and ocean formation were superimposed. Multidirectional tectonic movements, both vertical and horizontal, which led to the formation of the contours of oceanic basins and the ridges separating them, also influenced the continental margin. The destruction of the continental crust of the outer part of the continental shelf resulted in the formation of graben-like and rift troughs. The history of the Arctic oceanic basin began with the opening of the straits that connected the once isolated basin, the waters of which were largely desalinated. The latest stage of development was characterized by processes in which the transgressions and regressions of the Arctic Basin, the development and degradation of terrestrial and underground glaciation, and other processes played a leading role. [ABSTRACT FROM AUTHOR]

Subduction‐induced metasomatism may increase heterogeneity in the subcontinental lithospheric mantle (SCLM) over time. Syenite and its relevant varieties are typically derived from the most metasomatized portions of the SCLM and thus provide a key to decipher its heterogeneity. The only circum‐cratonic syenite province in the world, the Late Triassic belt around North China craton, is unique and offers an opportunity to investigate heterogeneity of the SCLM right before its destruction in the Mesozoic. New and previous data of whole‐rock Sr‐Nd and zircon in situ Hf‐O isotopes of syenite intrusions encircling the craton indicate two distinct types of metasomatism characterized by an alteration of juvenile or ancient materials. The spatiotemporally distinct sources can be discriminated by melts from oceanic slabs subducting beneath much of North China followed by deep subduction of South China craton only in the south. This unique setting of a circum‐cratonic syenite belt proves in a single case study that subduction‐induced metasomatism may cause craton‐scale SCLM heterogeneity. Plain Language Summary: The subcontinental lithospheric mantle (SCLM), comprising the roots of cratons, is depleted, rigid, and buoyant, all of which contribute to maintaining the stability of a craton. However, the SCLM can be replenished and destroyed by melts from recycled subducting slabs, known as mantle metasomatism. SCLM heterogeneity can be caused by metasomatism by carbonate and silicate melts derived from subducted oceanic and continental plates. It can then be recorded by mantle xenoliths or mantle‐derived rocks, but such records are typically either highly localized or asynchronous with other such records. In this study, Triassic syenites of North China craton that are distributed craton‐wide and coeval provide a chance to investigate SCLM heterogeneity through whole‐rock Sr‐Nd isotopes and zircon in situ Hf‐O isotopes. We propose that the contrasting subduction systems related to the Paleo‐Asian and Paleo‐Qinling oceanic plates subducting on the northern and southern margins of North China craton, respectively, as well as the deep continental subduction of South China craton to the southeast of North China, gave rise to the observed SCLM heterogeneity during the Permo–Triassic, which provided a preconditioning for the later Mesozoic lithospheric destruction of the cratonic root. Key Points: A unique circum‐cratonic Triassic syenite belt from North China craton (NCC) reveals heterogeneity of the subcontinental lithospheric mantleDifferential subduction systems surrounding NCC during the Permo–Triassic caused the heterogeneity through mantle metasomatismEarly Mesozoic mantle metasomatism served as a preconditioning for lithospheric destruction of NCC before the effect of the Pacific Ocean [ABSTRACT FROM AUTHOR]

We presented two high‐resolution, three‐dimensional (3D) global P‐ and joint S‐wave velocity models of the whole mantle with a particular focus on the North Pacific realm using our new global tomographic inversion method. The 3D ray tracing based on data‐adaptive block parameterization was implemented to calculate ∼17 million body‐wave absolute travel times and ray paths of multiple classes (e.g., P, S, PP, SSS, PcP, ScSScS, Pdiff, and PKP) recorded by local networks in China, Alaska, and global stations. Rayleigh wave data at 40–250 s was also incorporated for resolution improvement. Compared to previous models, especially at long spatial wavelengths, our new models better exhibited rich variations in the properties of slabs and cratons subducting/subducted into the mantle. The Kurile‐Kamchatka‐Aleutian arc exhibits a distinct and well‐defined pattern of long linear bands characterized by high velocity (high‐Vp, high‐Vs) anomalies. These anomalies descend continuously from the subduction zone and then extend horizontally at the mantle transition zone (MTZ) beneath the Bering Sea, Japan Sea, or alternatively penetrate into the lower mantle below the Okhotsk Sea. We observed there was no slab‐related gap trapped on top of the MTZ beneath the Alaska Peninsula, which was possibly relevant for the subduction of the extinct Kula and Pacific Plates described in previous studies. A further unexpected discovery, deduced from the slab's location and tectonic history, was the identification of a high‐Vp anomaly running parallel to the plate boundary in the Aleutian forearc region, which was considered as an integral part of the Kula slab. Plain Language Summary: The Circum‐North Pacific contains an active North‐Northwestward subduction zone of the ancient Izanagi and Kula Plates and the present Pacific Plate during the Late Mesozoic to Cenozoic. This successive subduction system leads to complex plate tectonics characterized by substantial topographic variations, widespread subduction zones, abundant intraplate volcanic activities, and ancient craton lithospheric deformation. Here we image high‐resolution P‐ and S‐wave velocity models, which were inverted by our updated global tomographic method using data containing more than 17 million body and surface wave measurements. Our models indicate no gaps at the top of the mantle transition zone, and it differs from the previous studies beneath the Aleutian Arc. In addition, two parallel high‐Vp belts are imaged on both sides of the Aleutian Arc between ∼300 and 1,000 km depth, which were considered as part of the remnant Kula slab and present Pacific slab. These findings provide both a good understanding of the slab geometry and evolution associated with multiple plate subduction in the North Pacific subduction zone and an example to understand the way how the Earth's interior works. Key Points: Global P‐ and S‐wave tomography models were determined for the North Pacific by the joint inversion of body wave and Rayleigh wave dataMassive‐slab‐related high‐velocity anomalies revealed the geometries of the subduction slabs under the North Pacific realmTwo subparallel high‐velocity belts have presumably related to the subduction of the Pacific and Kula slabs [ABSTRACT FROM AUTHOR]

The Satpura Mountain Belt, or the Central Indian Tectonic Zone (CITZ), is a Paleoproterozoic mountain belt of Central India, with Mesoproterozoic and Neoproterozoic rejuvenation history. Structural studies carried out in the southern margin of this ancient orogenic belt have identified a two-phase evolution. The documentation of mapped fold interference of Type 2 'boomerang-shaped' pattern and the spread of early lineations over the hinge area of second-generation folds provide an estimate of ~40-43° angle between the regional axial directions of the folds of the two orogenic episodes. The map patterns indicate the presence of basement-involved recumbent fold nappes placed above the sedimentary cover sequence. The identified décollement structure on the surface provides information regarding the deformation of the upper crust in ductile manner. The patterns of the deep crustal reflectors mapped through previous deep seismic reflection study interpreted during the present investigation brought out the thick-skinned deformation in the region as well as identification of crustal low-velocity zones of lower strength and a domal geometry of the lower crust. The regional thrusts are found to deform deeper crustal layers, one thrust affecting the crust-mantle boundary. The behaviours of the upper and lower crust in terms of rheological layering have been used to interpret the regional structural style of the lithosphere. The structural features of the region do not support the argument regarding the status of the Central Indian Shear as a suture, but indicate it to be the outermost emergent thrust towards the foreland side. The patterns of dip of the thrusts and thrust transport directions indicate a northward dip of the South Indian Block during the subduction stage and the position of the collisional suture of the Satpura orogen to the north of the Sausar belt. A tectonic model has been proposed here to explain the formation of the Satpura Mountain Belt or CITZ involving subduction-collision at ~2250 Ma, followed by erosion, stretching and basin development between 2100 and 1500 Ma. A collisional event in the eastern and western craton margins of India, at ~1500 Ma, caused the second uplift of the Satpura Mountain Belt and its northern margin, developing a temporal unconformity of ~300 Ma gap between the Semri and Kaimur Groups (Vindhyan Supergroup). A third event in the region at ~950 Ma was similar to the second event. No continental amalgamation took place along the CITZ during the tectonothermal events at ~1500 Ma and ~950 Ma; these involved reactivation of old thrusts and basal accretion which were the effect of far-field stresses emanating from the orogenic events at the eastern margin of the Bastar Craton and western margin of the Aravalli Craton. [ABSTRACT FROM AUTHOR]

The Svecofennian gabbro intrusions coincide temporally with the global 2100–1800 Ma orogens related to the amalgamation of the Mesoproterozoic supercontinent Nuna. We provide a new reliable 1891–1875 Ma palaeomagnetic pole for Fennoscandia based on rock magnetic and palaeomagnetic studies on the Svecofennian intrusions in central Finland to fill gaps in the Palaeoproterozoic palaeomagnetic record. By using the new pole together with other global high-quality data, we propose a new palaeogeographic reconstruction at 1885 Ma. This, together with previous data, supports a long-lived relatively stable position of Fennoscandia at low to moderate latitudes at 1890–1790 Ma. Similar stable pole positions have also been obtained for Kalahari at 1880–1830 Ma, Siberia at 1880–1850 Ma, and possibly India at 1980–1775 Ma. A new reconstruction at the beginning of this period indicates the convergence of several cratons at 1885 Ma in the initial stages of the amalgamation of the Nuna supercontinent at low to moderate latitudes. The close proximity of cratons at low to moderate latitudes is further supported by global and regional palaeoclimatic indicators. Stable position of several cratons could indicate a global period of minimal apparent drift at ca. 1880–1830 Ma. Before this period, the global palaeomagnetic record indicates large back-and-forth swings, most prominently seen in the high-resolution 2020–1870 Ma Coronation loops of the Slave craton. These large back-and-forth movements have been explained as resulting from an unstable geomagnetic field or basin- or local-scale vertical-axis rotations. However, the most likely explanation is inertial interchange true polar wander (IITPW) events, which is in line with the suggestion of large amplitude true polar wander events during the formation of the supercontinent. [ABSTRACT FROM AUTHOR]

Five main stages of the formation and reworking of continental crust in the region are recognized on the basis of correlation between the U–Pb isotopic ages of detrital zircon grains from present-day alluvial sediments of the largest rivers of West Transbaikalia and North Mongolia and primary rocks of the Mongolian–Transbaikalian sector of the Central Asian Foldbelt: (i) Neoarchean–Paleoproterozoic (2.7–1.7 Ga), corresponding to global crust-forming processes at the basement of the Precambrian cratons; (ii) Neoproterozoic (Early and Late Baikalian) and (iii) Early Paleozoic (Caledonian), which reflect the processes of the formation of juvenile crust by means of island arc magmatism and accretion of island arcs leading to the origination of fold structures of the southern margin of the Siberian Craton; and (iv) Late Paleozoic and (v) Early to Middle Mesozoic, which reflect the processes of reworking of the Early Precambrian and Caledonian crust by intraplate magmatism and riftogenesis. It is shown that the maximum ages on the probability density curve of detrital zircon grains depend on various (often unrelated) factors, e.g., the absolute age of the rock, the area of the eroded surface, the amount of zircon in rocks of various compositions, the transportation distance, etc., and the statistics of age distribution therefore could not be used as an indicator of intensity of endogenic events. The absence of Cenozoic detrital zircon grains in alluvial sediments of the region with abundant basaltic volcanism of this period (Mongolia, Baikal region) indicates that the endogenic events, which mostly produced basalts, are omitted in the "zircon record." [ABSTRACT FROM AUTHOR]

This study presents results from apatite fission track (AFT) thermochronology to investigate the thermal history and exhumation dynamics of the Rio Negro-Juruena basement, situated within the western Guiana Shield of the Amazonian Craton. AFT dating and associated thermal history modelling in South America has largely been restricted to the plate's margins (e.g. Andean active margin, Brazilian passive margin and others). Our paper reports on low-temperature thermochronological data from the internal part of the western Guiana Shield for the first time. This area is part of a vast cratonic lithosphere that is generally thought to be stable and little influenced by Mesozoic and Cenozoic tectonics. Our data however show AFT central ages ranging from 79.1 ± 3.2 Ma to 177.1 ± 14.8 Ma with mean confined track lengths of ca. 12 µm. Contrary to what might be expected of stable cratonic shields, inverse thermal history modeling indicates a rapid basement cooling event in the early Cretaceous. This cooling is interpreted as a significant exhumation event of the basement that was likely driven by the coeval extensional tectonics associated with back-arc rifts in the Llanos and Putumayo-Oriente-Maranon basins. The extensional tectonics facilitated both basement uplift and subsidence of the adjoining basins, increasing erosional dynamics and consequent exhumation of the basement rocks. The tectonic setting shifted in the late Cretaceous from extensional to contractional, resulting in reduced subsidence of the basins and consequential diminishing cooling rates of the Guiana Shield basement. Throughout the Cenozoic, only gradual, slow subsidence occurred in the study area due to regional flexure linked to the Andean orogeny. Comparative analysis with low-temperature thermochronology data from other West Gondwana cratonic segments highlights that exhumation episodes are highly controlled by tectonic inheritance, lithospheric strength, and proximity to rift zones. This study underscores the complex interplay between tectonic events and the response of cratonic lithosphere over geological time scales and highlights extensional settings as an important geological context for craton exhumation. [ABSTRACT FROM AUTHOR]

The Dom Feliciano Belt of southern Brazil and Uruguay represents part of a larger Neoproterozoic orogenic system formed during the amalgamation of Western Gondwana. The hinterland and foreland domains in parts of the belt preserve deformation structures and metamorphic assemblages that developed during early crustal thickening from c. 650 Ma. However, the metamorphic history of the southern foreland, in Uruguay, and its relationship with the hinterland, is not so well understood. We show that metamorphism in the southern hinterland is characterized by near‐isothermal decompression from ~10 kbar (~770°C) down to ~6 kbar, reflecting exhumation from depths of ~40 km during convergent thrusting and crustal thickening. This metamorphic event and associated magmatism is constrained by garnet Lu–Hf and zircon U–Pb dating to c. 655–640 Ma, supporting age and P–T constraints from previous studies. In contrast, prograde metamorphism in the foreland supracrustal rocks reached maximum lower‐amphibolite facies conditions (~6–7 kbar and ~550–570°C) and is constrained by garnet Lu–Hf dating to 582 ± 23 Ma. An exposed sheet of imbricated foreland basement rocks reached partial melting at upper‐amphibolite facies conditions, and metamorphism is similarly constrained to c. 585–570 Ma by monazite U–Pb dating. The data indicate that metamorphism in the foreland occurred during a sinistral transpressional event c. 55–85 Ma after the start of crustal thickening recorded in the hinterland, whereby strain partitioning during sinistral transpression led to imbrication in the foreland and oblique thrusting of the basement over more distal supracrustal rocks. This event is coeval with transpressional deformation in the Kaoko and Gariep belts, indicating a distinct two‐stage tectonic history driven by the three‐way convergence between the Congo, Kalahari, and South American cratons. [ABSTRACT FROM AUTHOR]

One of the most striking geological features of the Pamir is the south‐dipping lithospheric slab beneath the orogen characterized by an intracontinental Wadati‐Benioff zone. A widely accepted hypothesis over the past 40 years interprets the slab to represent southward subducted cratonic Asian continental lithosphere, which predicts significant cratonic Asia‐sourced crustal materials (e.g., Tarim Basin) beneath the Pamir. Alternatively, recent studies have interpreted the slab to be lithosphere delaminated from the base of the Pamir. To test these hypotheses, depth–tectonic affinity relations of crustal xenoliths carried by Miocene volcanic rocks in the eastern Pamir, interpreted to be sourced from the Pamir deep lithosphere, are used to determine whether they represented Asian affinity cratonic crust. Thermodynamic calculations, zircon U–Pb geochronology combined with rare earth element analysis, and whole‐rock major‐trace element and Sr–Nd isotopic analyses document that (1) eclogite and pyroxenite xenoliths (~31–43 kbar/~960–1170°C) are the deepest sourced portions of the lithosphere from ~100 to 140 km depth, the protoliths of which represent the mid‐lower crustal rocks of the Cretaceous Pamir magmatic arc, rather than material from cratonic Asia, and (2) granulite xenoliths (~20 kbar/~900°C) represent the Cenozoic lower crustal rocks of Pamir terranes from ~70 km depth. These results indicate the south‐dipping slab represents delaminated Pamir lower crust and mantle lithosphere, rather than intracontinental subduction of Asian lithosphere, and further support the hypothesis of minimal Cenozoic northward translation of the Pamir. [ABSTRACT FROM AUTHOR]

Archean gneiss domes are a very common structure in old cratons, these are different from the structural patterns in the post‐Archean linear orogenic belts related to subduction and collision. There are two tectonic models in Archean tectonics: horizontal tectonism and vertical tectonism. The formation mechanisms of Archean gneiss domes are still debated. Here we present data from the Anziling gneiss dome, a typical gneiss dome structure from the eastern North China Craton, which was formed dominantly by vertical tectonism during Neoarchean times (2554–2452 Ma). The Anziling gneiss dome is located in the east of Shuangshanzi, East Heibei Province and consists of the Anziling gneiss (a diopside‐monzonitic gneiss) in the center outwards surrounded by the Jiguanshan gneiss (hornblende‐plagioclase gneiss) and the Niuxinshan gneiss (trondhjemitic gneiss), together showing an irregular elliptical shape with a NE‐ward extending tail on the surface. But in the cross‐section, it represents a deflective dome, both north‐western and south‐eastern limbs dip approximately to the west. The dome suffered inhomogeneous ductile shear deformation, much stronger with more displacement in the north‐western half (especially close to the margin) than that in the eastern half. The core area shows only a relatively weak deformation. Therefore, we suggest that the Anziling gneiss dome was formed by the TTG (trondhjemite, tonalite, granodiorite) pluton rising obliquely up towards the east by diapirism in the first stage, followed by sub‐horizontal emplacement towards the north and then to the north‐east associated with clockwise rotation in the second stage. The Anziling gneiss dome and its neighbouring Shuangshanzi shear zone were formed simultaneously, indicating that the Anziling gneiss dome was driven by both gravitational forces and magma diapirism. [ABSTRACT FROM AUTHOR]

Outcrops of coarse‐grained, high‐angle clinoform sets are mainly thought to represent Gilbert‐type deltas. Superficially similar clinoform sets may, however, form in marine, sub‐wave base settings. They are interpreted to have formed as a result of storms where downwelling, seaward‐directed currents transported sand from the coastal area and shoreface across the shelf in suspension or as bedload to be deposited as clinothems. An additional transport of sand took place by strong coast‐parallel currents. The clinoform sets appear to be associated with rift events, the creation of accommodation space and an increasing supply of coarse‐grained sediment. A major protracted rift phase was initiated in East Greenland in Middle Jurassic times and intensified through the Late Jurassic to reach a climax close to the Jurassic–Cretaceous boundary. Rifting caused uplift of borderlands, creation of accommodation space and development of shallow marine shelves passing offshore into submarine slopes and deeper basinal areas. Each rift event was accompanied by the formation of clinoform sets prograding seawards towards the east and southeast away from the cratonic coastline in water depths below the wave base. The clinoform sets are interpreted as typical motifs for rift events in the relatively shallow epeiric Jurassic seaway between East Greenland and Norway. Outcrops of such sets represent a little‐known, commonly misinterpreted sedimentary system and may serve as motifs for rifting in shallow marine areas elsewhere in the geological record. Similar sets have been recorded in outcrop from the Mediterranean, and elsewhere and are probably more common than hitherto realized. [ABSTRACT FROM AUTHOR]

Seismic anisotropy tomography is the updated geophysical imaging technology that can reveal 3-D variations of both structural heterogeneity and seismic anisotropy, providing unique constraints on geodynamic processes in the Earth's crust and mantle. Here we introduce recent advances in the theory and application of seismic anisotropy tomography, thanks to abundant and high-quality data sets recorded by dense seismic networks deployed in many regions in the past decades. Applications of the novel techniques led to new discoveries in the 3-D structure and dynamics of subduction zones and continental regions. The most significant findings are constraints on seismic anisotropy in the subducting slabs. Fast-velocity directions (FVDs) of azimuthal anisotropy in the slabs are generally trench-parallel, reflecting fossil lattice-preferred orientation of aligned anisotropic minerals and/or shape-preferred orientation due to transform faults produced at the mid-ocean ridge and intraslab hydrated faults formed at the outer-rise area near the oceanic trench. The slab deformation may play an important role in both mantle flow and intraslab fabric. Trench-parallel anisotropy in the forearc has been widely observed by shear-wave splitting measurements, which may result, at least partly, from the intraslab deformation due to outer-rise yielding of the incoming oceanic plate. In the mantle wedge beneath the volcanic front and back-arc areas, FVDs are trench-normal, reflecting subduction-driven corner flows. Trench-normal FVDs are also revealed in the subslab mantle, which may reflect asthenospheric shear deformation caused by the overlying slab subduction. Toroidal mantle flow is observed in and around a slab edge or slab window. Significant azimuthal and radial anisotropies occur in the big mantle wedge beneath East Asia, reflecting hot and wet upwelling flows as well as horizontal flows associated with deep subduction of the western Pacific plate and its stagnation in the mantle transition zone. The geodynamic processes in the big mantle wedge have caused craton destruction, back-arc spreading, and intraplate seismic and volcanic activities. Ductile flow in the middle-lower crust is clearly revealed as prominent seismic anisotropy beneath the Tibetan Plateau, which affects the generation of large crustal earthquakes and mountain buildings. [ABSTRACT FROM AUTHOR]

Extension tectonics responsible for intracratonic rift basin formation are often the consequences of active or passive tectonic regimes. The present work puts forth a plume-related rifting mechanism for the creation and evolution of two Proterozoic sedimentary basins outlining the Bundelkhand craton, namely the Bijawar and Vindhyan basins. Using global gravity data, a regional scale study is performed over the region encompassing the southern boundary of the Bundelkhand craton consisting of Bijawar basin, Vindhyan basin and Deccan basalt outcrops. The gravity highs in the central part of the observed Bouguer gravity anomaly as well as the upward continued regional anomaly, derived from global gravity grid data, suggests that the Vindhyan sedimentary basin overlies a deeper high-density crustal source. The deepest interface as obtained from the radially averaged power spectrum analysis is observed to occur at a depth of ~30.3 km, indicating that the sources responsible for the observed gravity signatures occur at larger depths. 3D inversion of Bouguer gravity anomaly data based on Parker-Oldenburg's algorithm revealed the Moho depth of ~32 km below the Vindhyan basin, i.e., south of the craton. 2D crustal models along two selected profiles showcase a thick underplated layer with maximum thickness of ~12 km beneath the southern part of the Bundelkhand craton. The inferred large E-W trending underplating and deciphered shallower Moho beneath the regions south of the exposed Bundelkhand craton points to crustal thinning compensated by magmatic emplacement due to a Paleoproterozoic plume activity below the craton margin. [ABSTRACT FROM AUTHOR]

Shear wave splitting measurement of teleseismic data has been used to determine the fast polarization directions and delay times for 38 temporary stations and 15 permanent stations from a NW linear seismic array across the eastern Tarim basin (ETB) and the northern Tibetan Plateau (NTP), and 10 permanent stations on both sides of the array. We present an image of upper mantle anisotropy in the ETB and NTP using the 63 new measurements. The results show that the fast directions and delay times have complex spatial distribution characteristics. The delay times within the interior of the Tarim basin are very small, with an average value of 0.6 s, which is not only smaller than that in the Altyn Tagh fault and Tianshan on the southern and northern margins of the basin, but also smaller than that in the NTP, reflecting that the delay time of stable blocks is smaller than that of active blocks. Along the array, from east to west, the fast directions contrarotate from NNW in the southern Songpan-Garze terrane to NW in the northern Songpan-Garze terrane, to near E-W or ENE in the north of the East Kunlun fault and southern margin of the Qaidam basin, then first abruptly rotate to NW in the Qiman Tagh fault on the northwestern margin of the Qaidam basin, second abruptly rotate to ENE in the Altyn Tagh fault and south of the ETB, and third abruptly rotate to NW in the north of the ETB, then finally rotate to WNW in the Tianshan. The comparative analysis between the fast wave directions measured by shear wave splitting and predicted from the surface deformation field shows that, with the exclusion of the five observations with larger misfits within the interior of the ETB (with an average misfit of 27°), the misfits in the NTP and northern and southern margins of the Tarim basin are relatively small (with an average misfit of 9°). In addition, the fast wave directions of the tectonic units such as the Altyn Tagh fault, East Kunlun fault, and Tianshan are parallel to the strikes of faults and mountains in the region, which indicates that the deep and shallow deformations of the NTP and northern and southern margins of the ETB are consistent, where the crust-mantle coupling extent of lithospheric deformation is higher, according with the vertical coherent deformation of the lithosphere. Conversely, the crust-mantle coupling extent within the interior of the Tarim Basin is weak, and it is characterized by weak anisotropy, stable rigidity, and thick lithosphere, which may remain the "fossil" anisotropy of ancient craton. [ABSTRACT FROM AUTHOR]

KLR‐1 scientific deep borehole, drilled in the 1993 Killari earthquake region of Maharashtra (India), penetrated 8 m thick infratrappean sediments below a 338 m‐thick column of Deccan volcanics that rested over Neoarchean amphibolite to granulite facies mid‐crustal basement. Detailed palynological studies revealed the presence of characteristic Gondwana palynomorphs in sedimentary succession that include Plicatipollenites indicus, P. trigonalis, Parasaccites korbaensis, P. diffusus, Microbaculispora indica, Cyclobaculisporites minimus, Cyclogranisporites barakarensis, C. gondwanensis, Jayantisporites pseudozonatus, Leiosphaeridia bokaroensis, Latosporites sp., Verrucosisporites sp., cf. V. ambiplicatus and Striatopodocarpites gondwanensis. Qualitative analyses of palynomorphs would suggest an Asselian age (298–295 Ma) for their formation and burial in dysoxic to anoxic conditions. Recorded palynological assemblages correlate well with continental Lower Gondwana (Early Permian) sediments of India, as well as other Gondwana continents, indicating a larger extent of Gondwana sedimentation than hitherto known today. Petrological studies indicate the presence of two types of carbonate units that overly the crystalline basement. A relatively coarse‐grained sparry carbonate rock immediately overlies the basement (Sample KIL‐7). The size of the carbonate crystals in the sparite vary from 100 to 500 μm. This unit is overlain by an extremely fine‐grained micrite resembling lime mud (Sample KIL‐6). Micrite is made of carbonate grains with an average size of 10 μm. X‐ray diffraction studies indicate that calcite is the dominant carbonate mineral, while dolomite and quartz occur in minor amounts. Geochemical analyses indicate that in these two samples of carbonate rock that overlie the basement, CaO (33.38 to 38.16 wt %) is higher than MgO (16.0–23.85 wt%). These infratrappeans sediments are also characterized by stable carbon (δ13C) and oxygen (δ18O) isotopic ratios ranging from −1.59 to −8.74‰ and −9.12 to −13.50‰ versus PDB, respectively, thereby conforming to non‐marine fresh water deposits. These findings provide first evidences of reactivation and rifting of Dharwar intracratonic regions during earliest Permian Gondwana Period, which led to severe lithospheric mantle destruction, and massive regional uplift and erosion of the granitic‐gneissic upper crust. Large scale inland‐rifting during the Gondwana sedimentation period was possibly a prelude to India‐Antarctic breakup during the early Cretaceous. [ABSTRACT FROM AUTHOR]

'Transitional kimberlite' is a collective term previously used to classify rocks occurring in southern Africa that show bulk rock geochemical and Sr–Nd isotope features intermediate between (cratonic) lamproites and kimberlites. However, it is now well established that detailed petrographic and mineral chemical criteria represent a more robust approach towards the classification of kimberlites, lamproites and related rocks. Here, we re-assess the classification of southern African 'transitional kimberlites' by combining new petrographic observations and mineral compositional results for samples from six localities (Leicester, Frank Smith, Wimbledon, Melton Wold, Droogfontein, and Silvery Home) straddling the southwestern margin of the Kaapvaal Craton. These new data indicate that Leicester and Frank Smith are archetypal kimberlites, whereas Wimbledon, Melton Wold, Droogfontein, and Silvery Home represent bona fide olivine lamproites. We combine the mineral chemical results with new (Wimbledon) and existing bulk rock trace element and Nd–Hf isotope compositions, and emplacement ages, to assess whether the previously documented trends in Nd–Hf isotope vs time for these 'transitional kimberlites' constrain their petrological evolution. Modal groundmass mineralogy, bulk rock K/La and chromite compositions, the latter being a proxy for primitive melt composition, are linearly correlated with emplacement age and initial Nd–Hf isotope compositions. These observations suggest derivation of both older lamproites (181–115 Ma) and younger kimberlites (114–93 Ma), from a common evolving source. The temporal evolution of Nd–Hf isotope compositions in these rocks converge to values typical of archetypal Cretaceous kimberlites elsewhere in the Kaapvaal Craton, but are clearly different from the isotopic compositions of on-craton Kaapvaal lamproites (previously known as orangeites). This observation distinguishes the petrogenesis of the Wimbledon, Melton Wold, Droogfontein, and Silvery Home lamproites from those of 'typical' Kaapvaal lamproites. We hypothesize that progressive consumption of enriched and hence fertile K-bearing components in a common sub-lithospheric (i.e. convective mantle) source beneath the southwestern margin of the Kaapvaal Craton might represent a plausible scenario to explain the temporal evolution of petrographic and geochemical traits of the examined lamproites and kimberlites. A source in the lithospheric mantle is considered at odds with the contrasting location of the current localities as they occur both off- and on-craton. Migration of the African plate between 180 and 90 Ma over a relatively stationary convective mantle (plume?) source is not compatible with the spatial–temporal distribution of 'transitional kimberlites'. Instead, we invoke viscous coupling between an upper asthenospheric source and the lithosphere to reconcile a single evolving source with the geographic and age distribution of these rocks. This work supports the hypothesis that olivine lamproites occurring in intra-continental settings share similar genetic features with kimberlites. [ABSTRACT FROM AUTHOR]

Paleoproterozoic banded iron formation (BIF) iron ore of the Zhengjiapo region of the Changyi metallogenic belt, Eastern Block of North China Craton contains abundant coexisting antiperthite and mesoperthite textures. The antiperthite and mesoperthite occur in felsic domains of the Zhengjiapo BIF ore and enable derivation of peak temperature metamorphic conditions. Thermodynamic phase modeling shows that equilibrium conditions of corresponding textures, considering the related mineral assemblage of Pl + Qz + Kfs + Mag + Opx + L, are in the range of 870–940 °C and 5.0–6.8 kbar. Ternary feldspar thermometry using reintegrated compositions of antiperthite and mesoperthite in the felsic domain of the studied BIF iron ore reveals even higher peak metamorphic temperatures of 1045–1080 °C. The ultra-high temperature–low pressure conditions of Precambrian BIF have not yet been reported from the North China Craton. The documented ultra-high temperature metamorphism of the Zhengjiapo BIF iron ore in the Changyi metallogenic belt indicates that the BIF was involved in the collision-related tectonic process during Paleoproterozoic to have occurred in the Jiao-Liao-Ji orogenic belt. [ABSTRACT FROM AUTHOR]

The closure timing of the Paleo-Asian Ocean and the terminal stage of the Central Asian Orogenic Belt have been widely debated in the geological community. It's known that the gradual scissor-like closure of the Paleo-Asian Ocean occurred from west to east during the Paleozoic period. However, it was not until the Triassic period that the complete closure of the ocean occurred at the northeastern margin of the North China Craton. Nevertheless, there is still much uncertainty regarding the Late Triassic tectonic setting in Northeast China. In this study, we focused on the Upper Triassic Dajianggang Formation, located in the Shuangyang area of central Jilin Province, which is situated on the northeastern margin of the North China Plate. Our aim was to determine the formation age of the Dajianggang Formation by analyzing the detrital particle composition, petrogeochemistry, detrital zircon U-Pb isotope dating, and apatite fission track thermochronology. Our results indicated that the primary sandstone provenance area of the Dajianggang Formation in the Shuangyang area is the island arc orogenic belt. The tectonic background of the sandstone provenance area is mainly a continental island arc environment. The provenance area is mostly composed of felsic rocks with sedimentary tendencies, and some of its material may have originated from the northern margin of the North China Craton or the eroded recycle orogenic belt. LA-ICP-MS U-Pb dating of detrital zircons shows that the Dajianggang Formation formed after 226.8 ± 5 Ma. Moreover, analysis of the thermal evolution history modelling shows that the Dajianggang Formation in the Shuangyang area continued to be deposited and heated in the early stage, and then experienced rapid exhumation around 30 Ma. This suggests that the study area underwent an orogenic process during the early stage of formation, but then transitioned into a post-orogenic extension period, which constrained the final closure of the Paleo-Asian Ocean prior to the Late Triassic period. In addition, our study indicates that the remote effect of the Pacific subduction did not reach the study area until 30 Ma. The central age of the detrital apatite fission track of sample 19DJ-1 is 94.2 ± 8.3 Ma, which is younger than its corresponding stratigraphic age. The two peak ages of the fission track analysis are 62.9 ± 5.4 Ma and 126 ± 11 Ma. These findings provide new evidence for the tectonic evolution of Northeast China and shed light on the Late Triassic tectonic setting, as well as the influence time of subsequent tectonic domains in the southern part of Northeast China. [ABSTRACT FROM AUTHOR]

The volcanic activity of the Xiong'er Group and its concomitant sedimentation are related to the stretching–breakup of the Columbia supercontinent. The Dagushi Formation overlies the Paleoproterozoic Shuangfang Formation with an angular unconformity. The Dagushi Formation, as the earliest clastic strata of the Xiong'er Group and the first stable sedimentary cover overlying the Archean crystalline basement in the southern margin of the North China Craton, provides tectonic evolution information that predates Xiong'er volcanic activity. By distinguishing lithologic characteristics and sedimentary structures, we identified that the sedimentary facies of the Dagushi Formation were braided river delta lake facies from bottom to top. The U–Pb ages of the detrital zircons of the Dagushi Formation can be divided into four groups: ~1905–1925, ~2154–2295, ~2529–2536, and ~2713–2720 Ma, indicating the provenance from the North China Craton basement. Based on the geochemical characteristics of the Dagushi Formation, we suggest that the sediments accumulated rapidly near the source, which were principally felsic in nature, and were supplemented by recycled materials. The provenance area pointed to the underlying metamorphic crystalline basement of the North China Craton as the main source area with an active tectonic background. The Chemical Index of Alteration (CIA) values of the Dagushi Formation sandstone samples ranged from 60.8 to 76.7, indicating that the source rocks suffered from slight to moderate chemic chemical weathering. The Index of Composition Variability (ICV) values ranged from 0.8 to 1.3, which indicates the first cyclic sediments. The vertical facies and provenance changes of the Dagushi Formation reflect a continuous crust fracturing process that occurred in the North China Craton. [ABSTRACT FROM AUTHOR]

The Jiao-Liao-Ji Belt (JLJB) is the most representative Paleoproterozoic orogenic belt in the North China Craton (NCC). The sedimentation, metamorphism and magmatism of the Ji'an Group and associated granites provide significant insights into the tectonic evolution of the JLJB. In this study, we have synthesized published geochemistry and geochronology data on metasedimentary, metavolcanic and igneous rocks. According to the available data, the protoliths of the metasedimentary rocks are sets of shale, wacke, arkose, quartz sandstone and carbonate, while the protoliths of the metavolcanic rocks are calc-alkaline basalt, basaltic andesite, andesite, dacite and rhyolite. The rock assemblages indicate a transformation of the tectonic environment from a passive margin to an active continental margin following the onset of plate convergence and subduction. The A2-type gneissic granite (Qianzhuogou pluton) is formed in a subsequent back-arc basin extension setting at 2.20–2.14 Ga. The Ji'an Group was finally deposited in an active continental margin during the closure of a back-arc basin at 2.14–2.0 Ga. Then, the sediments were involved in a continent–arc–continent collision between the Longgang and Nangrim blocks at ~1.95 Ga. This process was accompanied by HP granulite-facies metamorphism at ~1.90 Ga. The subsequent exhumation and regional extension resulted in decompression melting during 1.90–1.86 Ga, producing metamorphism with an isothermal decompression clockwise P–T path. The resulting metapelites are characterized by perthite + sillimanite, and mafic granulites are characterized by orthopyroxene + clinopyroxene. The S-type porphyritic granite (Shuangcha pluton) is formed during the crustal anatexis. Meanwhile, extensive anatexis produced significant heating and triggered prograde to peak metamorphism with an anticlockwise P–T path. Cordierite-bearing symplectites around the garnet in the metapelites indicate a superposed isobaric cooling metamorphism. The ages of monazites and anatectic zircons suggest that the post-exhumation cooling occurred at 1.86–1.80 Ga. The Paleoproterozoic magmatism, sedimentation and metamorphism suggest a process of subduction back-arc basin extension and closure, collision and exhumation for the tectonic evolution of the JLJB. [ABSTRACT FROM AUTHOR]

We present a new S-wave velocity tomographic model of the upper mantle beneath the Australian Plate and its boundaries that we call Aus22. It includes azimuthal anisotropy and was constrained by waveforms from 0.9 million vertical-component seismograms, with the densest data sampling in the hemisphere centred on the Australian continent, using all available data covering this hemisphere. Waveform inversion extracted structural information from surface waves, S- and multiple S-waves and constrained S- and P-wave speeds and S-wave azimuthal anisotropy of the crust and upper mantle, down to the 660-km discontinuity. The model was validated by resolution tests and, for particular locations in Australia with notable differences from previous models, independent inter-station measurements of surface-wave phase velocities. Aus22 can be used to constrain the structure and evolution of the Australian Plate and its boundaries in fine detail at the regional scale. Thick, high-velocity (and, by inference, cold) cratonic lithosphere occupies nearly all of western and central Australia but shows substantial lateral heterogeneity. It extends up to the northern edge of the plate, where it collides with island arcs, without subducting. Diamondiferous kimberlites and lamproite deposits are underlain by cratonic lithosphere, except for the most recent diamondiferous lamproites in the King Leopold Orogen. The rugged eastern boundary of the cratonic lithosphere resolved by the model provides a lithospheric definition of the Tasman Line. Just east of the Tasman Line, an area of intermediate-thick lithosphere is observed in the southern part of the continent. The eastern part of Australia is underlain by thin, warm lithosphere, evidenced by low seismic velocities. All the sites of Cenozoic intraplate volcanism in eastern Australia are located on thin lithosphere. A low-velocity anomaly is present in the mantle transition zone (410–660 km depths) beneath the Lord Howe and Tasmanid hotspots, indicative of anomalously high temperature and consistent with a deep mantle upwelling feeding these hotspots and, possibly, also the East Australia hotspot. High seismic velocities at 200–410 km depth below New Guinea indicate the presence of slab fragments, probably linked to the subduction of the Australian Plate. High seismic velocities are observed in the transition zone below northeast Australia and indicate the presence of subducted lithospheric fragments trapped in the transition zone, possibly parts of the former northern continental margin of Australia. [ABSTRACT FROM AUTHOR]

We identified a ca. 180 Ma diamondiferous lamproite event in Zambia, establishing a link between ultrapotassic volcanism and the early Jurassic Karoo flood basalt province of sub-Saharan Africa. The cratonic lamproites erupted through the Permo–Triassic Luangwa Rift structure, but MgO-rich ultrapotassic magma formation was unrelated to rifting and triggered by plume–lithosphere interactions during the Karoo LIP event. Elevated Li–Zn–Ti concentrations in magmatic olivine (up to 18.5 ppm Li at 86–90 mol.% forsterite) and strong Sr–Nd–Hf–Pb isotopic enrichment of the host lamproites (87Sr/86Sr = 0.70701–0.70855, εNd = −10.8 to −10, εHf = −20.3 to −19.1, 206Pb/204Pb = 16.8–17.5) suggest partial melting of phlogopite-metasomatized lithospheric mantle domains, at approximately 180–200 km depth. The mantle-like δ7Li values (+2.8 to +5.7‰) of the most pristine lamproite samples are compatible with source enrichment by asthenosphere-derived melts, without significant involvement of recycled sedimentary components. This geochemical fingerprint stands in sharp contrast to the negative δ7Li compositions of primitive K-rich volcanic rocks from collision zone settings, where the shallow mantle sources contain recycled sediment. Isotope modelling demonstrates that the sub-Saharan lamproites originate from a MARID-style metasomatized peridotitic mantle source that underwent incompatible element enrichment at ca. 1 Ga, during tectonic activity associated with Rodinia supercontinent formation. Plume-sourced basaltic and picritic magmas of the 180 Ma Karoo LIP interacted with such K-rich hydrous lithospheric mantle domains, thereby attaining enriched incompatible element and radiogenic isotope compositions. Nd–Hf isotope mass balance suggests that up to 25% of MARID-sourced lamproite melt component contributed to some of the high-Ti flood volcanic units. Although large quantities of volatiles can be transferred from Earth's mantle to the atmosphere via plume–lithosphere interactions, it is unlikely that outgassing of mantle-sourced sulphur can exceed the climatic impact caused by the release of much more abundant carbon from thick continental roots. Thus, the excess SO 2 required to account for transient atmospheric cooling during the early Jurassic, coincident with the Karoo LIP event, must have had a thermogenic origin near the surface of Earth. [ABSTRACT FROM AUTHOR]

Metasomatized peridotites entrained by kimberlites are of great importance in understanding the geochemical changes occurring in the sub-continental lithospheric mantle (SCLM) during percolation and interaction of various fluids with the harzburgitic mantle. To simulate a scenario where slab-derived silicic fluid infiltrates the SCLM, high-pressure, high-temperature experiments were performed on a harzburgite + silicic fluid system at an 80:20 rock:fluid ratio (bulk ∼4.5 wt% H 2 O and ∼3 wt% CO 2). Experiments were performed at 3–6 GPa and 900–1200 °C on a rocking multi-anvil apparatus, conditions corresponding to the deep SCLM. The compositions of the mineral assemblages were analyzed in order to evaluate the formation of the metasomatic rocks at these conditions. At 3–6 GPa between 900–1000 °C, phlogopite-K-richterite- peridotites (PKP) were formed, at 3–5 GPa and 1100 °C, phlogopite-peridotites (PP), and at 6 GPa between 1100–1200 °C, garnet-phlogopite peridotites (GPP) were formed. Similar mineral assemblages and mineral chemistry are found in metasomatic xenoliths from Kimberly, South Africa, and surrounding localities, registering similar P-T conditions. The results presented in this study, together with previous studies on various peridotite + H 2 O ± CO 2 systems, clarify the relationship between bulk composition characteristics (e.g., K 2 O/H 2 O, peralkalinity) and the stability of different aluminum-bearing hydrous phases (e.g., phlogopite and amphibole). Fluxing harzburgite with silicic high-density fluids (HDFs) resulted in a peralkaline system, stabilizing both phlogopite and K-richterite. The stabilization of these Al-bearing hydrous phases led to the formation of a Cr-rich spinel, shifting the spinel to garnet transition to higher pressures. In addition, fluxing the harzburgite with silicic HDF led to the oxidation of the SCLM by one log unit, similar to the f O 2 registered by natural metasomatic xenoliths. [ABSTRACT FROM AUTHOR]

The North China Craton (NCC) is one of the most complex cratons in the world. It underwent a series of tectonothermal events during the Neoarchean–Paleoproterozoic. The petrogenesis of potassic granitoids, the timing, and the style of Archean crustal growth are still debated. Systematic field and petrological stdudies on the potassic meta-granites from the Guyang-Chayouzhongqi region were carried out. New U-Pb ages, zircon Lu-Hf isotopic analyses, and whole-rock geochemical data were obtained. Two groups (∼2.7 Ga and ∼2.5 Ga) of potassic meta-granites were recognized. The ∼2.7 Ga meta-granites are mainly A2-type, with variable εHf(t) values (−8.4 to +3.3) and Archean one stage model ages (TDM = ∼3.0 Ga), indicating that their source was derived from ancient anatectic TTG-like granite and depleted mantle, which suggests that thin crust had formed in the Guyang-Chayouzhongqi region by ∼3.0 Ga. Similar to the K-rich granites in the NCC, most of the ∼2.5 Ga potassic meta-granites are typical of A1-type granite, and are enriched in Sm and Gd and depleted in Nb, Ta, P, and Ti. The ages and isotopic data indicate that the ∼2.5 Ga meta-granites were generated from juvenile crustal sources with Neoarchean TTGs. The overall zircon U-Pb and Hf isotopic data furthermore suggested that the ∼2.7 Ga event is the most important stage of magmatic accretion in the NCC, similar to other cratons. In contrast, reworking or metamorphic alteration was the main crustal process in the NCC at ∼2.5 Ga. [ABSTRACT FROM AUTHOR]

The chemical composition of the deep continental crust is key to understanding the formation and evolution of the continental crust. Constraining the chemical composition of present‐day deep continental crust is, however, limited by indirect accessibility. This paper presents a modeling method for constraining deep crustal chemical structures from observed crustal seismic structures. We compiled a set of published composition models for the continental crust to construct functional relationships between seismic wave speed and major oxide content in the crust. Phase equilibria and compressional wave speeds (VP) for each composition model were calculated over a range of depths and temperatures of the deep crust. For conditions within the alpha(α)‐quartz stability field, robust functional relationships were obtained between VP and major oxide contents of the crust. Based on these relationships, observed VP of the deep crust can be inverted to chemical compositions for regions with given geotherms. We provide a MATLAB code for this process (CalcCrustComp). We apply this method to constrain compositions from deep crustal VP of global typical tectonic settings and the North China Craton (NCC). Our modeling results suggest that the lower crust in subduction‐related and rifting‐related tectonic settings may be more mafic than platforms/shields and orogens. The low VP signature in the deep crust of the NCC can be explained by intermediate crustal compositions, higher water contents, and/or higher temperatures. The chemical structure obtained by this method can serve as a reference model to further identify deep crustal features. Plain Language Summary: Understanding the formation and evolution of the continental crust requires knowledge of the chemical composition of the deep continental crust. Due to limited accessibility, it is difficult to determine the composition of deep crust. Here we present a method for estimating the composition of deep continental crust by combining published crust composition models, thermodynamic phase equilibrium calculations, and seismic wave speed calculations. Published model crustal compositions were used to establish functional relationships between crustal composition and seismic wave speed. Through the application of these functional relationships, it is possible to construct a chemical structure of the deep continental crust in regions subjected to seismic wave constraints. We apply this method to the global lower continental crust in typical tectonic settings and the North China Craton. The influence of crustal composition, geotherm, and water content on crustal VP is evaluated. Key Points: We present a geophysical‐geochemical modeling method to derive a chemical structure of deep continental crust from its VP structureThe global lower continental crust and the North China Craton are examined using this methodHeterogeneity in the crust is highlighted by the inferred chemical structure [ABSTRACT FROM AUTHOR]