Mesozoic lithospheric mantle of the northeastern Siberian craton (evidence from inclusions in kimberlite)

Several thousand clinopyroxene, garnet, and phlogopite inclusions of mantle rocks from Jurassic and Triassic kimberlites in the northeastern Siberian craton have been analyzed and compared with their counterparts from Paleozoic kimberlites, including those rich in diamond. The new and published mineral chemistry data make a basis for an updated classification of kimberlite-hosted clinopyroxenes according to peridotitic and mafic (eclogite and pyroxenite) parageneses. The obtained results place constraints on the stability field of high-Na lherzolitic clinopyroxenes, which affect the coexisting garnet and decrease its Ca contents. As follows from analyses of the mantle minerals from Mesozoic kimberlites, the cratonic lithosphere contained more pyroxenite and eclogite in the Mesozoic than in the Paleozoic. It virtually lacked ultradepleted harzburgite-dunite lithologies and contained scarce eclogitic diamonds. On the other hand, both inclusions in diamond and individual eclogitic minerals from Mesozoic kimberlites differ from eclogitic inclusions in diamond from Triassic sediments in the northeastern Siberian craton. Xenocrystic phlogopites from the D’yanga pipe have ⁴⁰Ar/³⁹Ar ages of 384.6, 432.4, and 563.4 Ma, which record several stages of metasomatic impact on the lithosphere. These phlogopites are younger than most of Paleozoic phlogopites from the central part of the craton (Udachnaya kimberlite). Therefore, hydrous mantle metasomatism acted much later on the craton periphery than in the center. Monomineral clinopyroxene thermobarometry shows that Jurassic kimberlites from the northeastern craton part trapped lithospheric material from different maximum depths (170 km in the D’yanga pipe and mostly < 130 km in other pipes). The inferred thermal thickness of cratonic lithosphere decreased progressively from ~ 260 km in the Devonian-Carboniferous to ~ 225 km in the Triassic and to ~ 200 km in the Jurassic, while the heat flux (Hasterok-Chapman model) was 34.9, 36.7, and 39.0 mW/m², respectively. Dissimilar PT patterns of samples from closely spaced coeval kimberlites suggest different emplacement scenarios, which influenced both the PT variations across the lithosphere and the diamond potential of kimberlites.     

Petrological-geophysical models of the internal structure of the lithospheric mantle of the Siberian Craton

Based on the simultaneous inversion of unique ultralong-range seismic profiles Craton, Kimberlite, Meteorite, and Rift, sourced by peaceful nuclear and chemical explosions, and petrological and geochemical data on the composition of xenoliths of garnet peridotite and fertile primitive mantle material, the first reconstruction was obtained for the thermal state and density of the lithospheric mantle of the Siberian craton at depths of 100–300 km accounting for the effects of phase transformation, anharmonicity, and anelasticity. The upper mantle beneath Siberia is characterized by significant variations in seismic velocities, relief of seismic boundaries, degree of layering, and distribution of temperature and density. The mapping of the present-day lateral and vertical variations in the thermal state of the mantle showed that temperatures in the central part of the craton at depths of 100–200 km are somewhat lower than those at the periphery and 300–400°C lower than the mean temperature of tectonically younger mantle surrounding the craton. The temperature profiles derived from the seismic models lie between the 32.5 and 35 mW/m² conductive geotherms, and the mantle heat flow was estimated as 11–17 mW/m². The depth of the base of the cratonic thermal lithosphere (thermal boundary layer) is close to the 1450 ± 100°C isotherm at 300 ± 30 km, which is consistent with published heat flow, thermobarometry, and seismic tomography data. It was shown that the density distribution in the Siberian cratonic mantle cannot be described by a single homogeneous composition, either depleted or enriched. In addition to thermal anomalies, the mantle density heterogeneities must be related to variations in chemical composition with depth. This implies significant fertilization at depths greater than 180–200 km and is compatible with the existence of chemical stratification in the lithospheric mantle of the craton. In the asthenosphere-lithosphere transition zone, the craton root material is not very different in chemical composition, thermal regime, and density from the underlying asthenosphere. It was shown that minor variations in the chemical composition of the cratonic mantle and position of chemical (petrological) boundaries and the lithosphere-asthenosphere boundary cannot be reliably determined from the interpretation of seismic velocity models only.     

Inhibited eclogitization and consequences for geophysical rock properties and delamination models: Constraints from cratonic lower crustal xenoliths

Studies on lower crustal and mantle xenoliths as well as geophysical data provide important information on the cratonic lithosphere. While geothermobarometric calculations of a majority of mantle xenoliths are in agreement with the typically low surface heat flow values of a craton (~ 40 mW/m²), P–T estimates for lower crustal xenoliths deviate significantly from the cratonic geotherms. Independent from the individual cratonic history, the temperatures are ~ 200–300 °C higher than what is expected at the base of the lower crust (~ 500–600 °C at ~ 1.3–1.6 GPa). Possible explanations may be a lack of equilibration to the cratonic geotherm or a relatively recent localized heat input. The presence of granulitic rocks under eclogite-facies conditions which are expected to prevail in the lower cratonic crust has consequences for the interpretation of geophysical rock properties. A mafic granulite which has been preserved under eclogite-facies conditions has densities and P-wave velocities similar to a felsic composition equilibrated to eclogite-facies conditions. Furthermore, phase diagrams calculated from xenolith bulk compositions demonstrate that eclogitization at relatively high temperatures as required for delamination of continental crust can only be triggered at significantly higher pressures than lithostatic at the base of the lower crust. As long as P–T conditions and the rock composition entail the assemblage to be granulitic, the addition of fluid at temperatures above 800 °C will not result in eclogitization, but rather in melt generation. This can also lead to an increase in density of up to 3%, however, this is strongly dependent on the amount of water saturation.     

Origin of high-velocity anomalies beneath the Siberian craton: A fingerprint of multistage magma underplating since the Neoarchean

Despite the violent eruption of the Siberian Traps at ~ 250 Ma, the Siberian craton has an extremely low heat flow (18–25 mW/m²) and a very thick lithosphere (300–350 km), which makes it an ideal place to study the influence of mantle plumes on the long-term stability of cratons. Compared with seismic velocities of rocks, the lower crust of the Siberian craton is composed mainly of mafic granulites and could be rather heterogeneous in composition. The very high V_p (> 7.2 km/s) in the lowermost crust can be fit by a mixture of garnet granulites, two-pyroxene granulites, and garnet gabbro due to magma underplating. The high-velocity anomaly in the upper mantle (V_p = 8.3-8.6 km/s) can be interpreted by a mixture of eclogites and garnet peridotites. Combined with the study of lower crustal and mantle xenoliths, we recognized multistage magma underplating at the crust-mantle boundary beneath the Siberian craton, including the Neoarchean growth and Paleoproterozoic assembly of the Siberian craton beneath the Markha terrane, the Proterozoic collision along the Sayan-Taimyr suture zone, and the Triassic Siberian Trap event beneath the central Tunguska basin. The Moho becomes a metamorphism boundary of mafic rocks between granulite facies and eclogite facies rather than a chemical boundary that separates the mafic lower crust from the ultramafic upper mantle. Therefore, multistage magma underplating since the Neoarchean will result in a seismic Moho shallower than the petrologic Moho. Such magmatism-induced compositional change and dehydration will increase viscosity of the lithospheric mantle, and finally trigger lithospheric thickening after mantle plume activity. Hence, mantle plumes are not the key factor for craton destruction.     

Upper mantle structure of the Saharan Metacraton

The ∼500,000 km² Saharan Metacraton in northern Africa (metacraton refers to a craton that has been mobilized during an orogenic event but that is still recognisable through its rheological, geochronological and isotopic characteristics) is an Archean–Paleoproterozoic cratonic lithosphere that has been destabilized during the Neoproterozoic. It extends from the Arabian–Nubian Shield in the east to the Trans-Saharan Belt in the west, and from the Oubanguides Orogenic Belt in the south to the Phanerozoic cover of North Africa. Here, we show that there are high S-wave velocity anomalies in the upper 100 km of the mantle beneath the metacraton typical of cratonic lithosphere, but that the S-wave velocity anomalies in the 175–250 km depth are much lower than those typical of other cratons. Cratons have possitive S-wave velocity anomalies throughout the uppermost 250 km reflecting the presence of well-developed cratonic root. The anomalous upper mantle structure of the Saharan Metacraton might be due to partial loss of its cratonic root. Possible causes of such modification include mantle delamination or convective removal of the cratonic root during the Neoproterozoic due to collision-related deformation. Partial loss of the cratonic root resulted in regional destabilization, most notably in the form of emplacement of high-K calc-alkaline granitoids. We hope that this work will stimulate future multi-national research to better understand this part of the African Precambrian. Specifically, we call for efforts to conduct systematic geochronological, geochemical, and isotopic sampling, deploy a reasonably-dense seismic broadband seismic network, and conduct systematic mantle xenoliths studies.     

Lithospheric structure of the North China Craton: Integrated gravity,geoid and topography data

The lithospheric structure of ancient cratons provides important constraints on models relating to tectonic evolution and mantle dynamics. Here we present the 3D lithospheric structure of the North China Craton (NCC) from a joint inversion of gravity, geoid and topography data. The NCC records a prolonged history of Archean and Paleoproterozoic accretion of crustal blocks through subduction and collision building the cratonic architecture, which was subsequently differentially destroyed during Mesozoic through extensive magmatism. The thermal structure obtained in our study is considered to define the lithosphere-asthenosphere boundary (LAB) of the NCC, and reflects the density variations within the mantle lithosphere. Employing the Moho depths from deep seismic sounding profiles for the inversion, and based on repeated computations using different parameters, we estimate the Moho depth, LAB depth and average crustal density of the craton. The Moho depth varies from 28 to 50 km and the LAB depth varies from 105 to 205 km. The LAB and Moho show concordant thinning from West to East of the NCC. The average crustal density is 2870 kg m^− 3 in the western part of the NCC, higher than that in the eastern part (2750 kg m^− 3). The results of joint inversion in our study yielded LAB depth and lithospheric thinning features similar to those estimated from thermal and seismic studies, although our results show different depth and variations in the thickness. The lithosphere gently thins from 145 to 105 km in the eastern NCC, where as the thinning is much less pronounced in the western NCC with average depth of about 175 km. The joint inversion results in this study provide another perspective on the lithospheric structure from the density properties and corresponding geophysical responses in an ancient craton.     

Local modification of the lithosphere beneath the central and western North China Craton: 3-D constraints from Rayleigh wave tomography

We have imaged the lithospheric structure beneath the central and western North China Craton (NCC) with Rayleigh wave tomography. The Rayleigh waveforms of 100 teleseismic events recorded by 208 broadband stations are used to yield high-resolution phase velocity maps at 13 periods from 20 s to 143 s. A 3-D S-wave velocity model is constructed based on the phase velocity maps. Our S-wave velocity model is broadly consistent with the results of previous tomography studies, but shows more detailed variations within the lithosphere. The Trans-North China Orogen (TNCO) is generally characterized by low-velocity anomalies but exhibits great heterogeneities. Two major low-velocity zones (LVZs) are observed in the north and south, respectively. The northern LVZ laterally coincides with sites of Cenozoic magmatism and extends to depths greater than 200 km. We propose that a small-scale mantle upwelling is present, confined to the north of the TNCO. A high-velocity patch in the uppermost mantle is also observed between the two LVZs adjacent to the narrow transtensional zone of the Cenozoic Shanxi–Shaanxi Rift (SSR). We interpret this as the remnant of a cratonic mantle root. The Ordos Block in the western NCC is associated with high-velocity anomalies, similarly reflecting the existence of cratonic mantle root, but a discernible low-velocity layer is observed at depths of 100–150 km in this location. We interpret that this mid-lithospheric structure was probably formed by metasomatic processes during the early formation of the NCC. Based on the observations from our S-wave velocity model, we conclude that the current highly heterogeneous lithospheric structure beneath the TNCO is the result of multiphase reworking of pre-existing mechanically weak zones since the amalgamation of the craton. The latest Cenozoic lithospheric reworking is dominated by the far-field effects of both Pacific plate subduction and the India–Eurasia collision.     

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

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

Thermochemical structure of the North China Craton from multi-observable probabilistic inversion: Extent and causes of cratonic lithosphere modification

We present a 3D thermochemical model of the North China Craton (NCC) from the surface down to 350 km by jointly inverting surface wave phase velocity data, geoid height, surface heat flow and absolute elevation with a multi-observable probabilistic inversion method. Our model reveals a thin (~ 65–100 km) and chemically fertile lithosphere (87 < Mg# < 90) beneath the Eastern NCC, consistent with independent results from mantle xenoliths, and supports the idea that the Eastern NCC experienced significant lithospheric destruction and refertilization during the Phanerozoic. In contrast, beneath the Trans-North China Orogen, Inner Mongolia Suture Zone and Yinshan belt, we observe a more heterogeneous (chemically and thermally) lithosphere, indicating that these areas have been partly involved in lithospheric modification and mechanical erosion at multiple scales. A cold and chemically refractory (Mg# > 90) lithospheric mantle is imaged beneath the central TNCO and Ordos Block, reaching depths > 260 km. This lithospheric “keel” is surrounded to the east by a high-temperature sublithospheric anomaly that originates at depths > 280 km. The spatial distribution of this anomaly and its correlation with the location of recent volcanism in the region suggest that the anomaly represents a deep mantle upwelling being diverted by the cratonic keel and spreading onto regions of shallow lithosphere. Our results indicate that the present-day thermochemical structure beneath the NCC is the result of a complex interaction between a large-scale return flow associated with the subduction of the Pacific slab and the shallow lithospheric structure.     

Seismic structure and rheology of the crust under mainland China

The crust and upper mantle in mainland China were relatively densely probed with wide-angle seismic profiling since 1958, and the data have provided constraints on the amalgamation and lithosphere deformation of the continent. Based on the collection and digitization of crustal P-wave velocity models along related wide-angle seismic profiles, we construct several crustal transects across major tectonic units in mainland China. In our study, we analyzed the seismic activity, and seismic energy releases during 1970 and 2010 along them. We present seismogenic layer distribution and calculate the yield stress envelopes of the lithosphere along the transects, yielding a better understanding of the lithosphere rheology strength beneath mainland China. Our results demonstrate that the crustal thicknesses of different tectonic provinces are distinctively different in mainland China. The average crustal thickness is greater than 65 km beneath the Tibetan Plateau, about 35 km beneath South China, and about 36–38 km beneath North China and Northeastern China. For the basins, the thickness is ~ 55 km beneath Qaidam, ~ 50 km beneath Tarim, ~ 40 km beneath Sichuan and ~ 35 km beneath Songliao. Our study also shows that the average seismic P-wave velocity is usually slower than the global average, equivalent with a more felsic composition of crust beneath the four tectonic blocks of mainland China resulting from the complex process of lithospheric evolution during Triassic and Cenozoic continent–continent and Mesozoic ocean–continent collisions. We identify characteristically different patterns of seismic activity distribution in different tectonic blocks, with bi-, or even tri-peak distribution of seismic concentration in South Tibet, which may suggest that crustal architecture and composition exert important control role in lithosphere deformation. The calculated yield stress envelopes of lithosphere in mainland China can be divided into three groups. The results indicate that the lithosphere rheology structure can be described by jelly sandwich model in eastern China, and crème brulee models with weak and strong lower crust corresponding to lithosphere beneath the western China and Kunlun orogenic belts, respectively. The spatial distribution of lithospheric rheology structure may provide important constraints on understanding of intra- or inter-plate deformation mechanism, and more studies are needed to further understand the tectonic process(es) accompanying different lithosphere rheology structures.     

The Mesoproterozoic mantle plume beneath the northern part of the Siberian craton

The study of the Mesoproterozoic (1473 ± 24 Ma) dolerites of the Olenek uplift of the Siberian craton basement has shown their petrologic and geochemical similarity to typical OIB produced with participation of a mantle plume. The dolerites are characterized by variations in the geochemical composition explained by different degrees of melting of the same source. A conclusion is drawn that the parental melts of the rocks were slightly modified by crustal contamination, as evidenced from their Nd isotope composition (£_Nd(T) = + 0.6 to − 0.8) and the presence of inherited zircons of four ages (2564, 2111, 2053, and 1865 Ma). Since the Siberian craton in the structure of the Nuna supercontinent (Columbia) was located relatively close to the Baltic continent and the Congo and Sao Francisco cratons, we assume that the Early Mesoproterozoic mafic intrusions (1500–1470 Ma) of all these cratons belong to the same large igneous province (LIP). The province formation was related to the activity of superplume (or mantle hot field), which supplied mantle matter to the lithosphere basement. The superplume core was probably located beneath the northern part of the Siberian craton, where basites are compositionally most similar to the primary mantle source.     

Amount of Asian lithospheric mantle subducted during the India/Asia collision

Body wave seismic tomography is a successful technique for mapping lithospheric material sinking into the mantle. Focusing on the India/Asia collision zone, we postulate the existence of several Asian continental slabs, based on seismic global tomography. We observe a lower mantle positive anomaly between 1100 and 900 km depths, that we interpret as the signature of a past subduction process of Asian lithosphere, based on the anomaly position relative to positive anomalies related to Indian continental slab. We propose that this anomaly provides evidence for south dipping subduction of North Tibet lithospheric mantle, occurring along 3000 km parallel to the Southern Asian margin, and beginning soon after the 45 Ma break-off that detached the Tethys oceanic slab from the Indian continent. We estimate the maximum length of the slab related to the anomaly to be 400 km. Adding 200 km of presently Asian subducting slab beneath Central Tibet, the amount of Asian lithospheric mantle absorbed by continental subduction during the collision is at most 600 km. Using global seismic tomography to resolve the geometry of Asian continent at the onset of collision, we estimate that the convergence absorbed by Asia during the indentation process is ~ 1300 km. We conclude that Asian continental subduction could accommodate at most 45% of the Asian convergence. The rest of the convergence could have been accommodated by a combination of extrusion and shallow subduction/underthrusting processes. Continental subduction is therefore a major lithospheric process involved in intraplate tectonics of a supercontinent like Eurasia.     

Isotope (U–Pb,Sm–Nd,Rb–Sr) geochronology of alkaline basic plutons of the Kuznetsk Alatau

On the northeastern slope of the Kuznetsk Alatau, small differentiated alkaline basic intrusive massifs form an isometric area ~ 100 km across. They are composed of subalkalic and alkali gabbroids, basic and ultrabasic foidolites, nepheline and alkali syenites, and carbonatites. Results of complex (U–Pb, Sm–Nd, and Rb–Sr) isotope dating suggest that alkaline basic magmatism developed at two stages, in the Middle Cambrian–Early Ordovician (~ 510–480 Ma) and in the Early–Middle Devonian (~ 410–385 Ma). Finding of accessory zircons (age 1.3–2.0 Ga) in alkaline rocks suggests that the ascent of mantle plume was accompanied by the melting of fragments of Proterozoic mature continental crust composing the basement of the Caledonian orogen of the Kuznetsk Alatau. Probably, parental Cambrian–Ordovician alkaline mafic melts initiated metasomatism and lithosphere erosion. During the next melting of lithosphere substrate in ~ 100 Myr, this caused the generation of magmas of similar composition with inherited isotope parameters (ε_Nd(T) ≈ + 4.8 to + 5.7, T_Nd(DM) ≈ 0.8–0.9 Ga) pointing to the similar nature of their matter sources in the moderately depleted mantle.     

Change in the viscosity of kimberlite and basaltic magmas during their origin and evolution (prediction)

Based on the analysis of experimental data on the viscosity of mafic to ultramafic magmatic melts with the use of our structure-chemical model for the calculation and prediction of the viscosity of magmas, we have first predicted that diamond-carryihg kimberlite magma must ascend from mantle to crust with considerable acceleration. The viscosity of kimberlite magma decreases by more than three times during its genesis, evolution, and ascent from mantle to crust despite the significant decrease in the temperature of the ascending kimberlite magma (~ 150 °C) and its partial crystallization and degassing. In the case of partial melting (< 1 wt.%) of carbonated peridotite in the mantle at depths of 250-350 km, high-viscosity (~ 35 Pas) kimberlite melts can be generated at ~ 8.5 GPa and ~ 1350 °C, the water content in the melt being up to ~ 8 wt.%, C(OH^-) = 0-2 wt.%, and C(H₂O) = 0-6 wt.%. On the other hand, during the formation of kimberlite pipes, dikes, and sills, the viscosity of near-surface kimberlite melts is much lower (~ 10 Pa s) at ~ 50 MPa and 1200 °C, the volume contents of crystals (V_cr) and the fluid phase (bubbles) (Vfl) are 35 and 5 vol.%, respectively, and the water content in magma, C(OH^-), is 0.5 wt.%. On the contrary, the viscosity of basaltic magmas increases by more than two orders of magnitude during their ascent from mantle to crust. The basaltic magmas which can be generated in the asthenosphere at depths of ~ 100 km have the minimum viscosity (up to ~ 2.3 Pas) at ~ 4.0 GPa, 1350 °C, C(OH^-) - 3 wt.%, and C(H₂O) - 5 wt.%. However, at the final stage of evolution (e.g., during volcanic eruptions), the viscosity of basaltic magma is considerably higher (600 Pa s) at ~ 10 MPa, 1180 °C, V_cr - 30 vol.%, Vf - 15 vol.%, and C(OH^-) - 0.5 wt.%.     

Tectonics of the North Qilian orogen,NW China

The Qilian Orogen at the northern margin of the Tibetan Plateau is a type suture zone that recorded a complete history from continental breakup to ocean basin evolution, and to the ultimate continental collision in the time period from the Neoproterozoic to the Paleozoic. The Qilian Ocean, often interpreted as representing the “Proto-Tethyan Ocean”, may actually be an eastern branch of the worldwide “Iapetus Ocean” between the two continents of Baltica and Laurentia, opened at ≥ 710 Ma as a consequence of breakup of supercontinent Rodinia.Initiation of the subduction in the Qilian Ocean probably occurred at ~ 520 Ma with the development of an Andean-type active continental margin represented by infant arc magmatism of ~ 517–490 Ma. In the beginning of Ordovician (~ 490 Ma), part of the active margin was split from the continental Alashan block and the Andean-type active margin had thus evolved to western Pacific-type trench–arc–back-arc system represented by the MORB-like crust (i.e., SSZ-type ophiolite belt) formed in a back-arc basin setting in the time period of ~ 490–445 Ma. During this time, the subducting oceanic lithosphere underwent LT-HP metamorphism along a cold geotherm of ~ 6–7 °C/km.The Qilian Ocean was closed at the end of the Ordovician (~ 445 Ma). Continental blocks started to collide and the northern edge of the Qilian–Qaidam block was underthrust/dragged beneath the Alashan block by the downgoing oceanic lithosphere to depths of ~ 100–200 km at about 435–420 Ma. Intensive orogenic activities occurred in the late Silurian and early Devonian in response to the exhumation of the subducted crustal materials.Briefly, the Qilian Orogen is conceptually a type example of the workings of plate tectonics from continental breakup to the development and evolution of an ocean basin, to the initiation of oceanic subduction and formation of arc and back-arc system, and to the final continental collision/subduction and exhumation.     

3D thermal structure of the continental lithosphere beneath China and adjacent regions

Based on the Crust2.0 model and the topography data of Chinese continent and its adjacent regions, a three-dimensional finite element model is constructed in terms of the spherical coordinate system. In our numerical model, the average annual ground temperature from 195 meteorological stations and temperature of upper mantle derived from the seismic velocities are adopted as the top and bottom boundary conditions, respectively. The observed thermal conductivity and heat production from P wave velocity based on empirical formula are employed in our numerical model as well. The comparison between the calculated and observed surface heat flow proved that our results are reliable. The temperature beneath the Precambrian cratons is lower than that of other areas for 100–300 °C also. The typical temperature rang at the Moho is estimated to be 800–1000 °C beneath the Tibetan plateau and 500–700 °C beneath the Precambrian cratons (such as Indian plate, Sichuan basin, South China, North China and Tarim), respectively. The thermal state in the eastern part of Sino-Korean craton at the depth deeper than 60 km indicates that it was destructed. The thermal structure in center of Tibetan plateau (especially beneath Qiangtang area) supports the proposed flow of lower crustal or upper mantle material to the east. Generally, the distribution of volcanoes in Chinese continent is consistent with the high temperature areas in the crust or upper mantle. There are many obvious thermal transition zones across the orogenic belts. The thermal transition zone between eastern and western parts in the crust of Chinese continent is consistent with the north–south seismic zone.     

Thermo-mechanical modeling of the obduction process based on the Oman Ophiolite case

Obduction emplaces regional-scale fragments of oceanic lithosphere (ophiolites) over continental lithosphere margins of much lower density. For this reason, the mechanisms responsible for obduction remain enigmatic in the framework of plate tectonics. We present two-dimensional (2D) thermo-mechanical models of obduction and investigate the possible dynamics and physical controls of this process. Model geometry and boundary conditions are based on available geological and geochronological data and numerical modeling results are validated against petrological and structural observations of the Oman (Semail) Ophiolite. Our model reproduces the stages of oceanic subduction initiation away from the Arabian margin, the emplacement of the Oman Ophiolite on top of it, and the domal exhumation of the metamorphosed margin through the ophiolitic nappe. A systematic study indicates that 350–400 km of bulk shortening provides the best fit for both maximum pressure–temperature conditions of the metamorphosed margin (1.5–2.5 GPa/450–600 °C) and the dimension of the ophiolitic nappe (~ 170 km width). Our results confirm that a thermal anomaly located close to the Arabian margin (~ 100 km) is needed to initiate obduction. We further suggest that a strong continental basement rheology is a prerequisite for ophiolite emplacement.     

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

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

The origin and geochemical evolution of the Woodlark Rift of Papua New Guinea

Geochemical, isotopic, and geochronologic data for exhumed rocks in the Woodlark Rift of Papua New Guinea (PNG) allow a tectonic link to be established with the Late Cretaceous Whitsunday Volcanic Province (WVP) of northeastern Australia. Most of the metamorphic rocks in the Woodlark Rift have Nd isotopic compositions (ε_Nd = + 1.7 to + 6.2) similar to the Nd isotopic compositions of rocks in the WVP (ε_Nd = + 1.3 to + 6.6; Ewart et al., 1992), and contain inherited zircons with 90 to 100 Ma U–Pb ages that overlap the timing of magmatism in the WVP. None of the metamorphic rocks in the Woodlark Rift have the highly evolved Hf and Nd isotopic compositions expected of ancient continental crust. Magmas were erupted in the WVP during the middle Cretaceous as eastern Gondwana was rifted apart. The protoliths of felsic and intermediate metamorphic rocks in the Woodlark Rift are interpreted to be related to the magmatic products produced during this Cretaceous rifting event. Some mafic metamorphic rocks exposed in the western Woodlark Rift (eclogites and amphibolites) are not related to the WVP and instead could have originated as basaltic lavas crystallized from mantle melts at (U)HP depths in the Late Cenozoic, or as fragments of Mesozoic aged oceanic lithosphere.Isotopic and elemental comparisons between basement gneisses and Quaternary felsic volcanic rocks demonstrate that felsic lavas in the D'Entrecasteaux Islands did not form solely from partial melting of metamorphic rocks during exhumation. Instead, the isotopic compositions and geochemistry of Quaternary felsic volcanic rocks indicate a significant contribution from the partial melting of the mantle in this region. When combined with geophysical data for the western Woodlark Rift, this suggests that future seafloor spreading will commence south of Fergusson Island, and west of the present-day active seafloor spreading rift tip.     

Dynamic subsidence of Eastern Australia during the Cretaceous

During the Early Cretaceous Australia's eastward passage over sinking subducted slabs induced widespread dynamic subsidence and formation of a large epeiric sea in the eastern interior. Despite evidence for convergence between Australia and the paleo-Pacific, the subduction zone location has been poorly constrained. Using coupled plate tectonic–mantle convection models, we test two end-member scenarios, one with subduction directly east of Australia's reconstructed continental margin, and a second with subduction translated ~ 1000 km east, implying the existence of a back-arc basin. Our models incorporate a rheological model for the mantle and lithosphere, plate motions since 140 Ma and evolving plate boundaries. While mantle rheology affects the magnitude of surface vertical motions, timing of uplift and subsidence depends on plate boundary geometries and kinematics. Computations with a proximal subduction zone result in accelerated basin subsidence occurring 20 Myr too early compared with tectonic subsidence calculated from well data. This timing offset is reconciled when subduction is shifted eastward. Comparisons between seismic tomography and model temperature cross-sections, and an absence of subduction zone volcanism in eastern Australia in the Early Cretaceous provide support for the back-arc basin scenario.