Lithospheric petrology of the eastern Arabian Plate: Constraints from Al-Ashkhara (Oman) xenoliths

Mafic granulite and spinel lherzolite xenoliths from Cenozoic alkaline basalts near Al-Ashkhara, eastern Oman, have been selected for a systematic mineralogical, geochemical and Sr–Nd–Pb isotopic study. This is the only place in E Arabia where samples of both lower crust and upper mantle can be examined. Lower crustal xenoliths consist of two mineralogically and chemically distinct groups: gabbronorite (subequal abundances of ortho- and clino-pyroxene and plagioclase) and plagioclase pyroxenite (dominant pyroxene and subordinate plagioclase). Temperature estimates for lower crustal xenoliths using the two pyroxene geothermometer (T-Wells) yield 810–865 °C. The mineral assemblage (spinel–pyroxene–plagioclase) and Al content in pyroxene indicate that plagioclase-bearing xenoliths equilibrated at 5–8 kbar (13 and 30 km depth) in the lower crust. εNd and ⁸⁷Sr/⁸⁶Sr calculated at 700 Ma for Al-Ashkhara lower crustal xenoliths (+ 6.4 to + 6.6; ⁸⁷Sr/⁸⁶Sr = 0.7028 to 0.7039) are consistent with the interpretation that juvenile, mafic melts were added to the lower crust during Neoproterozoic time and that there was no discernible contribution from pre-Neoproterozoic crust. Upper mantle xenoliths consist of both dry and hydrous (phlogopite-bearing) lherzolites. These peridotites are more Fe-rich than expected for primitive mantle or melt residues and probably formed by pervasive circulation of melts that have refertilized pre-existing mantle peridotites. Mineral equilibration temperatures range from 990 to 1070 °C. Isotopic compositions calculated at 700 Ma are εNd = + 6.8 to + 7.8 and ⁸⁷Sr/⁸⁶Sr = 0.7016 to 0.7025, indicating depleted upper mantle. Pb isotopic compositions indicate that the metasomatism was relatively recent, perhaps related to Paleogene tectonics and basanite igneous activity. Nd model ages for the spinel peridotite xenoliths range between 0.59 and 0.65 Ga. The xenolith data suggest that eastern Arabian lower crust is of hotspot origin, in contrast to western Arabian lower crust, which mostly formed at a convergent plate margin. Geochemical and isotopic differences between lower crust and upper mantle indicate that these are unrelated, possibly because delamination replaced the E Arabian mantle root in Neoproterozoic time.     

Lower crustal xenoliths from Junan,Shandong province and their bearing on the nature of the lower crust beneath the North China Craton

Geochronological, petrological and geochemical studies were performed on the granulite xenoliths from a Late Cretaceous basaltic breccia dike in Junan, Shandong province, eastern China. These xenoliths show close similarities to the Nushan granulite xenoliths from the southern margin of the North China Craton (NCC) and the Archean granulite terrains in terms of mineralogy and bulk rock compositions, but are quite different from the Hanuoba mafic granulite xenoliths from the northern NCC. In-situ zircon U–Pb age and Hf isotopic analyses, together with geochemical data reveal that the protolith of these xenoliths was formed around 2.3 Ga ago, through assimilation–fractional crystallization of a mafic magma. P–T conditions of these xenoliths suggest that the lower crust beneath the Junan region reaches to a depth of 35 km, which agree well with the result deduced from various geophysical methods. The consistent petrological and seismic Moho depths, the observed velocity structure and calculated velocity of these xenoliths imply the absence of underplating induced crust–mantle transition zone, which was well formed in the northern NCC. Compared to 40–50 km depth of the lower crust in Early Jurassic, the lower crust beneath Junan extended to a depth of 30 km in Late Cretaceous, suggesting that the lower crust of NCC was significantly thinned during Late Mesozoic.     

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.     

Crustal and lithospheric thinning beneath the seismogenic Kachchh rift zone,Gujarat (India): Its implications toward the generation of the 2001 Bhuj earthquake sequence

A high-resolution passive seismic experiment in the Kachchh rift zone of the western India has produced an excellent dataset of several thousands teleseismic events. From this network, 500 good teleseismic events recorded at 14 mobile broadband sites are used to estimate receiver functions (for the 30–310° back-azimuth ranges), which show a positive phase at 4.5–6.1 s delay time and a strong negative phase at 8.0–11.0 s. These phases have been modeled by a velocity increase at Moho (i.e. 34–43 km) and a velocity decrease at 62–92 km depth. The estimation of crustal and lithospheric thicknesses using the inversion of stacked radial receiver functions led to the delineation of a marked thinning of 3–7 km in crustal thickness and 6–14 km in lithospheric thickness beneath the central rift zone relative to the surrounding un-rifted parts of the Kachchh rift zone. On an average, the Kachchh region is characterized by a thin lithosphere of 75.9 ± 5.9 km. The marked velocity decrease associated with the lithosphere–asthenoshere boundary (LAB), observed over an area of 120 km × 80 km, and the isotropic study of xenoliths from Kachchh provides evidence for local asthenospheric updoming with pockets of partial melts of CO₂ rich lherzolite beneath the Kachchh seismic zone that might have caused by rifting episode (at 88 Ma) and the associated Deccan thermal-plume interaction (at 65 Ma) episodes. Thus, the coincidence of the area of the major aftershock activity and the Moho as well as asthenospheric upwarping beneath the central Kachchh rift zone suggests that these pockets of CO₂-rich lherzolite partial melts could perhaps provide a high input of volatiles containing CO₂ into the lower crust, which might contribute significantly in the seismo-genesis of continued aftershock activity in the region. It is also inferred that large stresses in the denser and stronger lower crust (at 14–34 km depths) induced by ongoing Banni upliftment, crustal intrusive, marked lateral variation in crustal thickness and related sub-crustal thermal anomaly play a key role in nucleating the lower crustal earthquakes beneath the Kachchh seismic zone.     

Mineral chemistry and zircon geochronology of xenocrysts and altered mantle and crustal xenoliths from the Aries micaceous kimberlite: Constraints on the composition and age of the central Kimberley Craton,Western Australia

The Neoproterozoic (∼ 820 Ma) Aries micaceous kimberlite intrudes the central Kimberley Basin, northern Western Australia, and has yielded a suite of 27 serpentinised ultramafic xenoliths, including spinel-bearing and rare, metasomatised, phlogopite–biotite and rutile-bearing types, along with minor granite xenoliths. Proton-microprobe trace-element analysis of pyrope and chromian spinel grains derived from heavy mineral concentrates from the kimberlite has been used to define a ∼ 35–40 mW/m² Proterozoic geotherm for the central Kimberley Craton. Lherzolitic chromian pyrope highly depleted in Zr and Y, and Cr-rich magnesiochromite xenocrysts (class 1), probably were derived from depleted garnet peridotite mantle at ∼ 150 km depth. Sampling of shallower levels of the lithospheric mantle by kimberlite magmas in the north and north-extension lobes entrained high-Fe chromite xenocrysts (class 2), and aluminous spinel-bearing xenoliths, where both spinel compositions are anomalously Fe-rich for spinels from mantle xenoliths. This Fe-enrichment may have resulted from Fe–Mg exchange with olivine during slow cooling of the peridotite host rocks. Fine exsolution rods of aluminous spinel in diopside and zircon in rutile grains in spinel- and rutile-bearing serpentinised ultramafic xenoliths, respectively, suggest nearly isobaric cooling of host rocks in the lithospheric mantle, and indicate that at least some aluminous spinel in spinel-facies peridotites formed through exsolution from chromian diopside. Fe–Ti-rich metasomatism in the spinel-facies Kimberley mantle probably produced high-Ti phlogopite–biotite + rutile and Ti, V, Zn, Ni-enriched aluminous spinel ± ilmenite associations in several ultramafic xenoliths. U–Pb SHRIMP ²⁰⁷Pb/²⁰⁶Pb zircon ages for one granite (1851 ± 10 Ma) and two serpentinised ultramafic xenoliths (1845 ± 30 Ma; 1861 ± 31 Ma) indicate that the granitic basement and lower crust beneath the central Kimberley Basin are at least Palaeoproterozoic in age. However, Hf-isotope analyses of the zircons in the ultramafic xenoliths suggest that the underlying lithospheric mantle is at least late Archean in age.     

Crustal shear-wave velocity structure beneath northeast India from teleseismic receiver function analysis

We investigated the seismic shear-wave velocity structure of the crust beneath nine broadband seismological stations of the Shillong–Mikir plateau and its adjoining region using teleseismic P-wave receiver function analysis. The inverted shear wave velocity models show ∼34–38 km thick crust beneath the Shillong Plateau which increases to ∼37–38 km beneath the Brahmaputra valley and ∼46–48 km beneath the Himalayan foredeep region. The gradual increase of crustal thickness from the Shillong Plateau to Himalayan foredeep region is consistent with the underthrusting of Indian Plate beyond the surface collision boundary. A strong azimuthal variation is observed beneath SHL station. The modeling of receiver functions of teleseismic earthquakes arriving the SHL station from NE backazimuth (BAZ) shows a high velocity zone within depth range 2–8 km along with a low velocity zone within ∼8–13 km. In contrast, inversion of receiver functions from SE BAZ shows high velocity zone in the upper crust within depth range ∼10–18 km and low velocity zone within ∼18–36 km. The critical examination of ray piercing points at the depth of Moho shows that the rays from SE BAZ pierce mostly the southeast part of the plateau near Dauki fault zone. This observation suggests the effect of underthrusting Bengal sediments and the underlying oceanic crust in the south of the plateau facilitated by the EW-NE striking Dauki fault dipping 30⁰ toward northwest.     

Destruction of the North China Craton: Delamination or thermal/chemical erosion? Mineral chemistry and oxygen isotope insights from websterite xenoliths

The destruction (or reactivation) of the North China Craton (NCC) is one of the important issues related to the Phanerozoic evolution of eastern China, although the processes of destruction remain debated. Two main mechanisms – delamination and thermal/chemical erosion – have been proposed based on the geochemistry of Mesozoic–Cenozoic basalts and entrained deep–seated xenoliths. A key criterion in distinguishing between these mechanisms is the nature of the melt, derived from delaminated crust or the asthenosphere, that modified the lithospheric mantle. Here we investigate the mechanism of destruction of the NCC based on mineral compositions and oxygen isotopic data from olivines, and strontium isotopic data for clinopyroxenes within websterite xenoliths from the Early Cretaceous Feixian basalts in the eastern NCC. Olivines in websterite xenoliths have higher Mg# (86–86.4), Ni content (2187–2468 ppm), and lower Ca (983–1134 ppm), Ti (58.1–76.1 ppm), and Mn (1478–1639 ppm) contents than olivine phenocrysts (Mg# = 71.0–77.3, Ni = 233–1038 ppm, Ca = 1286–2857 ppm, Ti = 120–300 ppm, and Mn = 2092–4106 ppm) from Late Cretaceous basalts. Additionally, olivines in websterite xenolith have δ¹⁸O values of 7.10 ± 0.21‰ to 8.40 ± 0.21‰, evidently higher than those of typical mantle-derived olivines. Similarly, orthopyroxenes (Opx) and clinopyroxenes (Cpx) in the websterite xenoliths have much higher Mg# (86.3–89.2 and 87.5–90.3, respectively), and Ni contents (1097–1491 ppm and 581–809 ppm, respectively) than orthopyroxene- and clinopyroxene-phenocrysts (Opx: Mg# = 82.2–83.9, Ni = 730–798 ppm; Cpx: Mg# = 74.2–84.6, Ni = 117–277 ppm) from Late Cretaceous basalts. The ⁸⁷Sr/⁸⁶Sr ratios of clinopyroxenes in the websterite xenoliths range from 0.70862 to 0.70979, and are much higher than those of clinopyroxene grains from peridotite xenoliths and basalts. These data indicate that olivines are the residue of ancient lithospheric mantle that was modified intensively by a melt derived from recycled continental crust, and that the silicic and calcic metasomatic melt might have been derived from the partial melting of the subducted Yangtze slab and delaminated lower crust of the NCC. The existence of recycled continental crust in the Mesozoic lithospheric mantle implies that delamination was an important mechanism of destruction of the NCC.     

Nature and processes of the lithospheric mantle beneath the western Qinling: Evidence from deformed peridotitic xenoliths in Cenozoic kamafugite from Haoti,Gansu Province,China

The Cenozoic Haoti kamafugite field (23 Ma) is situated at the western Qinling Orogen, Gansu Province in China, which is a conjunction region of the North China Craton, the Yangtze Craton and the Tibetan Plateau. Fresh peridotitic xenoliths entrained in these volcanic rocks provide an opportunity to study the nature and processes of the lithospheric mantle beneath the western Qinling. These xenoliths can be divided into two groups based on the petrological features and mineral compositions, type 1 and type 2. Type 1 xenoliths with strongly deformed texture have higher Fo (90–92.5) contents in olivines, Mg# (91–94) and Cr# (15–35) of clinopyroxenes, and Cr# (36–67) of spinels than the weakly deformed type 2 xenoliths, which have the corresponding values of 89–90, 89–91.5, 10–15 and 5–15 in minerals, respectively. CaO contents in fine-grained olivines are slightly higher than 0.10 wt% compared with coarse-grained ones (less than 0.10 wt%). Fine-grained clinopyroxenes have low Al₂O₃ + CaO contents (generally <23 wt%) relative to coarse-grained ones (>23 wt%). Fo contents in fine-grained olivines mainly in the melt pocket of the type 1 xenoliths are higher than those in coarse-grained ones, which is somewhat contrary to the type 2 xenoliths without melt pocket. Clinopyroxenes of the type 2 display higher Na₂O contents (1.7–1.9 wt%) than those of the type 1 (<1.4 wt%). P–T estimations reveal that the type 1 xenoliths give temperature in range of 1106–1187 °C and pressure of 21–26 kbar and that relatively low temperature (907 and 1022 °C) and pressure (19.0 and 18.5 kbar) for the type 2 xenoliths. The type 1 xenoliths are characterized by depletion due to high degree of partial melting (>10%), modal metasomatic and deformed characteristics, and may represent the old refractory lithospheric mantle. In contrast, the type 2 peridotites show fertile features with low degree of partial melting (<5%) and may represent the newly-accreted lithospheric mantle. The lithospheric mantle beneath the western Qinling underwent partial melting, recrystallization, deformation and metasomatism due to asthenospheric upwelling and the latest decompression responding to the Cenozoic extensive tectonic environment. These processes perhaps are closely related to the evolution of Tibetan Plateau caused by the India-Asian collision.     

Highly heterogeneous lithospheric mantle beneath the Central Zone of the North China Craton evolved from Archean mantle through diverse melt refertilization

High-Mg# peridotite xenoliths in the Cenozoic Hebi basalts from the North China Craton have refractory mineral compositions (Fo > 91.5) and highly heterogeneous Sr–Nd isotopic compositions (⁸⁷Sr/⁸⁶Sr = 0.7031–0.7048, ¹⁴³Nd/¹⁴⁴Nd = 0.5130–0.5118) ranging from MORB-like to EM1-type mantle, which are similar to those of peridotites from Archean cratons. Thus, the high-Mg# peridotites may represent relics of the ancient lithospheric mantle. Published Re–Os isotopic data for Cenozoic basalt-borne xenoliths show T_RD ages of 3.0–1.5 Ga for the peridotites from Hebi (the center of the craton), 2.2–0 Ga for those from Hannuoba and Jining (north margin of the craton), and 2.6–0 Ga for those from Fanshi and Yangyuan (midway between the center and north margin of the craton). In situ Re–Os data of sulfides in Hannuoba peridotites suggest that whole-rock Re–Os model ages represent mixtures of multiple generations of sulfides with varying Os isotopic compositions. These observations indicate that initial lithospheric mantle beneath the Central Zone of the North China Craton formed during the Archean and was refertilized by multiple melt additions after its formation. The refertilization became more intensive from the interior to the margin of the craton, leading to the high heterogeneity of the lithospheric mantle: more ancient and refractory peridotites with highly variable Sr–Nd isotopic compositions in the interior, and more young and fertile peridotites with depleted Sr–Nd isotopic composition in the margin. Our data, coupled with published petrological and geochemical data of peridotites from the Central Zone of the North China Craton, suggest that the lithospheric mantle beneath this region is highly heterogeneous, likely produced by refertilization of Archean mantle via multiple additions of melts/fluids, which were closely related to the Paleoproterozoic collision between the Eastern and the Western Blocks and subsequent circum-craton subduction events.     

An overview of the crustal structure of the Tibetan plateau after 35 years of deep seismic soundings

Since the pioneer wide-angle seismic profile along the Yadong–Gulu rift acquired in 1974 by the ex-Institute of Geophysics, Chinese Academy of Sciences (CAS), several research programs aimed to deep geophysics, performed thanks to the participation of Chinese national and international institutions, have been developed during last 35 years, including 23 wide-angle seismic profiles with total length of about 6000 km. These profiles are unevenly distributed, most of them in eastern Tibet and few profiles in western Tibet. In this paper, we make a summarized presentation of all these wide-angle seismic profiles and provide an overall view of the seismic velocity structure of the crust beneath the broad Tibetan plateau, which is the product of the continuous convergence and collision of the Indian and Eurasian plates since about 50 Ma ago. Different patterns of crustal thickness variation related to the tectonic blocks and along suture zones of the region are displayed. The crust thickness is confirmed to be about 70–75 km under southern Tibet, and 60–65 km under northern, northeastern and southeastern Tibet. The leading edge of the subducted lithosphere reaches the northern margin of the plateau and directly contacts with Tarim Basin. Westward of the 90°E boundary, the Indian crust is moving towards the northern edge of the plateau and collides with Tarim Basin at 80°E while reach the Bangong–Nujiang suture belt at 88°E; eastward of the 90°E boundary, the northern edge of the crust should be at 50–100 km south of Bangong–Nujiang suture. The results supply helpful constrains to understand the mechanism of the continent–continent collision and its consequences in the plateau and neighbouring areas.     

Seismic structure of the crust and uppermost mantle of South America and surrounding oceanic basins

We present a new set of contour maps of the seismic structure of South America and the surrounding ocean basins. These maps include new data, helping to constrain crustal thickness, whole-crustal average P-wave and S-wave velocity, and the seismic velocity of the uppermost mantle (P_n and S_n). We find that: (1) The weighted average thickness of the crust under South America is 38.17 km (standard deviation, s.d. ±8.7 km), which is ∼1 km thinner than the global average of 39.2 km (s.d. ±8.5 km) for continental crust. (2) Histograms of whole-crustal P-wave velocities for the South American crust are bi-modal, with the lower peak occurring for crust that appears to be missing a high-velocity (6.9–7.3 km/s) lower crustal layer. (3) The average P-wave velocity of the crystalline crust (P_cc) is 6.47 km/s (s.d. ±0.25 km/s). This is essentially identical to the global average of 6.45 km/s. (4) The average P_n velocity beneath South America is 8.00 km/s (s.d. ±0.23 km/s), slightly lower than the global average of 8.07 km/s. (5) A region across northern Chile and northeast Argentina has anomalously low P- and S-wave velocities in the crust. Geographically, this corresponds to the shallowly-subducted portion of the Nazca plate (the Pampean flat slab first described by Isacks et al., 1968), which is also a region of crustal extension. (6) The thick crust of the Brazilian craton appears to extend into Venezuela and Colombia. (7) The crust in the Amazon basin and along the western edge of the Brazilian craton may be thinned by extension. (8) The average crustal P-wave velocity under the eastern Pacific seafloor is higher than under the western Atlantic seafloor, most likely due to the thicker sediment layer on the older Atlantic seafloor.     

Crustal structure in an onshore–offshore transitional zone near Hong Kong,northern South China Sea

To study the crustal structure beneath the onshore–offshore transitional zone, a wide-angle onshore–offshore seismic experiment was carried out in northern South China Sea near Hong Kong, using large volume airgun sources at sea and seismic stations on land. The crustal velocity model constructed from traveltime fitting shows that the sedimentary thickness abruptly increases seaward of the Dangan Islands based on the characteristics of Pg and Multiple Pg, and the crustal structure beneath the sedimentary layer is relatively simple. The Moho depth is about 25–28 km along the profile and the P-wave velocity increases gradually with depth. The velocities in the upper crust range from 5.5 to 6.4 km/s, while that in the lower crust is 6.4–6.9 km/s. It also reveals a low velocity zone with a width of more than 10 km crossing the crust at about 75–90 km distance, which suggests that the Littoral Fault Zone (LFZ) exists beneath the onshore–offshore transitional zone. The magnetism anomalies, bouguer gravity anomalies and active seismic zone along the coastline imply the LFZ is a main tectonic fault in the onshore–offshore area. Combined with two previously published profiles in the continental South China (L–G profile) and in the northern margin of South China Sea (OBS1993) respectively, we constructed a land-sea super cross-section about 1000 km long. The results show the onshore–offshore transitional zone is a border separating the unstretched and the stretched continental crust. The low velocity layer (LVL) in the middle crust was imaged along L–G profile. However, the high velocity layer (HVL) in the lower crust was detected along OBS1993. By analyzing the mechanisms of the LVL in the middle crust and HVL in the base of crust, we believe the crustal structures had distinctly different attributes in the continental South China and in the northern SCS, which indicates that the LFZ could be the boundary fault between them.     

P,S wave velocity model of the crust and upper most mantle of Albania region

This paper describes the one-dimensional (1D) velocity model computed by VELEST in the SEISAN seismic analysis system, inverting re-picked P-wave and S-wave arrival times recorded during 2002–2006 by the Albanian, Montenegro, Thessalonica and Macedonia seismic networks. The re-picked data yield P-wave and S-wave velocities proved to be more suitable compared to bulletin data for this detailed inversion study. Seismic phases recorded by the Albania seismic network and integrated with data from the Montenegro, Thessalonica and Macedonia networks are used to prepare the Albanian seismic bulletin. Earthquake hypocenters from the Albanian bulletins have also location errors that are negligible for civil protection purposes, large scale seismotectonic analyses and more accurate hypocentral determinations which are necessary for detailed seismotectonic and geodynamic studies.It was noted that the smoothness of the velocity variation increased with depth. A velocity of 5.5 km/s was calculated for the upper crust, 6.1 km/s was calculated for the middle crust and 6.9 km/s was computed for the lower crust. P wave velocity was 7.85 km/s at depth of 50 km and for the upper mantle it is 8.28 km/s. Using the improved velocity model, the earthquakes which occurred in Albania in the past 5 years were able to be relocated, achieving constrained hypocentral determinations for events in Albania. The interpretation of the 1 D velocity models infers interesting features of the deep structure of Albania. These results represent an important step towards more detailed seismotectonic analyses.     

Mid-Neoproterozoic crustal evolution of the northeastern Yangtze Block: Evidence from the felsic-gneiss xenoliths hosted in the Donghai Cenozoic basalts

Crustal xenoliths can provide new insights into the unexposed crust, and those from the northeastern Yangtze Block have rarely been studied. This paper reports U–Pb–Hf isotopes and trace-element compositions of zircons from six felsic xenoliths hosted by the Neogene alkali basalts in the Donghai region (i.e. Anfengshan and Pingmingshan) of the Sulu orogen in central eastern China. The xenoliths are mainly composed of orthoclase and quartz, or orthoclase and natrolite, with accessory minerals of Fe–Ti oxides and zircon. Most zircon grains show core-rim structures, with the cores and rims being magmatic and metamorphic in origin, respectively. The zircon cores mainly yield ages of ca. 827–794 Ma, while the zircon rims give ages of ca. 232–212 Ma. We interpret the zircon core ages as the time of an early Mid-Neoproterozoic magmatic event in the northeastern Yangtze Block and the zircon rim ages as the time of collision between the Yangtze and North China Blocks. Our data suggest that much more ca. 830–800 Ma magmatic records are possibly preserved in the unexposed deep crust, and the early Mid-Neoproterozoic is an important era for the crust evolution of the northeastern Yangtze Block. The new zircon Hf isotopic analyses show that the Anfengshan sample (south of Donghai) has zircon ε_Hf (820 Ma) values ranging from −15.3 to −9.4, and two-stage Hf model ages of 2.66–2.30 Ga; the Pingmingshan sample (southeast of Donghai) has zircon ε_Hf (820 Ma) values ranging from −1.4 to +3.8, and two-stage Hf model ages of 1.80–1.47 Ga. These data suggest that ancient crust as old as Neoarchean to Mesoproterozoic was involved in the early Mid-Neoproterozoic magmatism. Combined with the previously reported zircon U–Pb–Hf results of the exposed rocks, it is highlighted that crustal recycling was dominant in the early Mid-Neoproterozoic (ca. 830–800 Ma) magmatism, whereas both crustal recycling and addition of mantle-derived melts were significant in the late Mid-Neoproterozoic (ca. 800–720 Ma) magmatism in the northeastern Yangtze Block.     

Mapping crustal structure beneath southern Tibet: Seismic evidence for continental crustal underthrusting

Receiver function imaging along a temporary seismic array (ANTILOPE-2) reveals detailed information of the underthrusting of the Indian crust in southern Tibet. The Moho dips northward from ~ 50 km to 80 km beneath the Himalaya terrane, and locally reaches ~ 85 km beneath the Indus–Yalung suture. It remains at ~ 80 km depth across the Lhasa terrane, and shallows to ~ 70 km depth under the Qiangtang terrane. An intra-crustal interface at ~ 60 km beneath the Lhasa terrane can be clearly followed southward through the Main Himalaya Thrust and connects the Main Boundary Thrust at the surface, which represents the border of the Indian crust that is underthrusting until south of the Bangong–Nujiang Suture. A mid-crustal low velocity zone is observed at depths of 14–30 km beneath the Lhasa and Himalaya terranes probably formed by partial melt and/or aqueous fluids.     

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.     

Distinct variations of crustal shear wave velocity structure and radial anisotropy beneath the North China Craton and tectonic implications

This study presents the crustal shear wave velocity structure and radial anisotropy along two linear seismic arrays across the North China Craton (NCC) from ambient noise tomography. About a half to one year long ambient noise data from 87 stations were used for obtaining the inter-station surface wave empirical Green's functions (EGFs) from cross-correlation. Rayleigh and Love dispersion curves within the period band 5–30 s were measured from the EGFs of the vertical and transverse components, respectively. These dispersion data were then used to determine the crustal shear wave velocity structure (V_SV and V_SH) and radial anisotropy (2(V_SH ? V_SV) / (V_SH + V_SV)) from point-wise linear inversion with constraints from receiver function analysis. Our results reveal substantial structural variations among different parts of the NCC. The Bohai Bay Basin in the eastern NCC is underlain by a thin crust (~ 30 km) with relatively low velocities (particularly V_SV) and large positive radial anisotropy in the middle to lower crust. Such a crustal structure is no longer of a cratonic type and may have resulted from the widespread tectonic extension and intensive magmatism in this region since late Mesozoic. Beneath the Ordos Basin in the western NCC, the crust is relatively thicker (≥ 40 km) and well stratified, and presents a large-scale low velocity zone in the middle to lower crust and overall weak radial anisotropy except for a localized lower crust anomaly. The overall structural features of this region resemble those of typical Precambrian shields, in agreement with the long-term stability of the region. The crustal structure under the Trans North China Orogen (TNCO, central NCC) is more complicated and characterized by smaller scale velocity variations, strong positive radial anisotropy in the middle crust and rapid change to weak-to-negative anisotropy in the lower crust. These features may reflect complex deformations and crust–mantle interactions, probably associated with tectonic extension and magmatic underplating during the Mesozoic to Cenozoic evolution of the region. Our structural images in combination with previous seismic, geological and geochemical observations suggest that the Phanerozoic lithospheric reactivation and destruction processes may have affected the crust (especially the middle and lower crust) of the eastern NCC, and the effect probably extended to the TNCO, but may have minor influence on the crust of the western part of the craton.     

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.     

Crustal and upper mantle S-wave velocity structures across the Taiwan Strait from ambient seismic noise and teleseismic Rayleigh wave analyses

Although orogeny tapers off in western Taiwan large and small earthquakes do occur in the Taiwan Strait, a region largely untouched in previous studies owing mostly to logistical reasons. But the overall crustal structure of this region is of particular interest as it may provide a hint of the proto-Taiwan before the orogeny.By combining time domain empirical Green’s function (TDEGF) from ambient seismic noise using station-pairs and traditional surface wave two-station method (TS) we are able to construct Rayleigh wave phase velocity dispersion curves between 5 and 120 s. Using Broadband Array in Taiwan for Seismology (BATS) stations in Taiwan and in and across the Strait we are able to derive average 1-D Vs structures in different parts of this region. The results show significant shear velocity differences in the upper 15 km crust as expected. In general, the highest Vs in the upper crust observed in the coastal area of Mainland China and the lowest Vs appears along the southwest offshore of the Taiwan Island; they differ by about 0.6–1.1 km/s. For different parts of the Strait, the upper crust Vs structures are lower in the middle by about 0.1–0.2 km/s relative to those in the northern and southern parts. The upper mantle Vs structure (Moho – 150 km) beneath the Taiwan Strait is about 0.1–0.3 km/s lower than the AK135 model. The overall crustal thickness is approximately 30 km, much thinner and less variable than under the Taiwan Island. The inversion of seismic velocity structures using shorter period band dispersion data in the sea areas with water depth deeper than 1000 m should take water layer into consideration except for the continental shelves.     

An overview of the earth crust under China

We assess the results of a number of deep seismic soundings performed in China over the last few decades, and study the variations in crustal structure in 18 tectonic units comprised of three platforms and 15 fold systems. Thickness data on 344 Mesozoic–Cenozoic sedimentary basins, as well as data on Moho depth are collected in order to discuss the relationship between the thickness of the basins and the average thickness of the consolidated crust in each tectonic unit. The degree of mirror-image symmetry between Mesozoic–Cenozoic sedimentary basins and the uplifting topmost parts of the mantle is herein analyzed using deep geophysical data on sedimentary basins. By applying standard methods of least-squares analysis to both datasets, we have obtained both the average thickness of the consolidated crust and the mirror-image symmetry factor for every platform and tectonic fold system, thereby allowing us to explore the correlation between the depths of the bottom of the sedimentary basins and the top of the uplifting mantle. The thickness of the consolidated crust in China is found to be between 20 and 63 km, following a pattern of gradual thickening from east to west. Expressed in terms of spatial seismicity and the concentration of seismic energy, and according to the sharing-out of earthquake hypocenters in the top 80 km of the earth, the rheology of the area does not appear to suit the widely accepted “jelly–sandwich” model for the continents, which does not seem to be entirely valid in China. The findings on the mirror-image symmetry factor show that this parameter varies mainly in the range ? 0.5 to ? 1.8 compared with each other tectonic unit. Most of the tectonic systems in China appear to be isostatically compensated. Tibet is an exception, in that the symmetry factor is positive for the Gangdise–Nyainqentanglha (1.0) and Himalayan (0.1) fold systems, implying that these tectonic zones are far from being in the isostatic equilibrium of the other regions. We have also analyzed the zoning characteristics through the logarithmic relationship ln R = ln (h/|a|), using the average thickness of the consolidated crust and the absolute value of the symmetry factor. The key finding is that regardless of geographical location, all the values of R = h/|a| for those sedimentary basins in which oil/gas reservoirs have to date been found, fall into the narrow range of values between 19.38 and 37.40. There is some appeal in a possible relationship involving the ratio of crustal thickness to symmetry factor, more so when the results obtained appear to suggest a prognostic tool for exploratory practice in relation to oil/gas reservoirs.