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.     

Lithospheric mantle structure of the Siberian craton inferred from the superlong Meteorite and Rift seismic profiles

Modeling of the seismic, thermal, and density structure of the Siberian craton lithospheric mantle at depths of 100-300 km has been performed along the superlong Meteorite and Rift seismic profiles. The 2D velocity sections reflect the specific features of the internal structure of the craton: lateral inhomogeneities, seismic-boundary relief at depths of ~ 100, 150, 240, and 300 km, velocities of 8.3-8.7 km/s, and the lack of low-velocity zone in the lower lithosphere. Mapping of the thermal state along the Meteorite and Rift profiles shows a significant temperature decrease in the cratonic mantle as compared with the average temperatures of the surrounding Phanerozoic mantle (> 300 °C) estimated from the global reference model AK135. Lateral temperature variations, reflecting the thermal anomalies in the cratonic keel, are observed at depths of < 200 km (with some decrease in temperature in the central part of the craton), whereas at depths of > 200 km, temperature variations are negligible. This suggests the preservation of residual thermal perturbations at the base of the lithosphere, which must lead to the temperature equalization in the transition zone between the lithosphere and the asthenosphere. Variations in chemical composition have a negligible effect on the thermal state but affect strongly the density structure of the mantle. The results of modeling admit a significant fertilization of matter at depths more than 180-200 km and stratification of the cratonic mantle by chemical composition. The thicknesses of chemical (petrologic) and thermal boundary layers beneath the Siberian craton are estimated. The petrologic lithosphere is localized at depths of ~ 200 km. The bottom of the thermal boundary layer is close to the 1450 °C isotherm and is localized at a depth of 300 km, which agrees with heat flow and seismic-tomography data.     

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.     

S-wave velocity structure of the North China from inversion of Rayleigh wave phase velocity

We constructed the S-wave velocity structure of the crust and uppermost mantle (10–100 km) beneath the North China based on the teleseismic data recorded by 187 portable broadband stations deployed in this region. The traditional two-step inversion scheme was adopted. Firstly, we measured the interstation fundamental Rayleigh wave phase velocity of 10–60 s and imaged the phase velocity distributions using the Tarantola inversion method. Secondly, we inverted the 1-D S-wave velocity structure with a grid spacing of 0.25° × 0.25° and constructed the 3-D S-wave velocity structure of the North China. The 3-D S-wave velocity model provides valuable information about the destruction mechanism and geodynamics of the North China Craton (NCC). The S-wave velocity structures in the northwestern and southwestern sides of the North–South Gravity Lineament (NSGL) are obviously different. The southeastern side is high velocity (high-V) while the northeastern side is low velocity (low-V) at the depth of 60–80 km. The upwelling asthenosphere above the stagnated Pacific plate may cause the destruction of the Eastern Block and form the NSGL. A prominent low-V anomaly exists around Datong from 50 to 100 km, which may due to the upwelling asthenosphere originating from the mantle transition zone beneath the Western Block. The upwelling asthenosphere beneath the Datong may also contribute to the destruction of the Eastern Block. The Zhangjiakou-Penglai fault zone (ZPFZ) may cut through the lithosphere and act as a channel of the upwelling asthenosphere. A noticeable low-V zone also exists in the lower crust and upper mantle lid (30–50 km) beneath the Beijing–Tianjin–Tangshan (BTT) region, which may be caused by the upwelling asthenosphere through the ZPFZ.     

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 Saharan Metacraton

This article introduces the name “Saharan Metacraton” to refer to the pre-Neoproterozoic––but sometimes highly remobilized during Neoproterozoic time––continental crust which occupies the north-central part of Africa and extends in the Saharan Desert in Egypt, Libya, Sudan, Chad and Niger and the Savannah belt in Sudan, Kenya, Uganda, Congo, Central African Republic and Cameroon. This poorly known tract of continental crust occupies 5,000,000 km² and extends from the Arabian-Nubian Shield in the east to the Tuareg Shield to the west and from the Congo craton in the south to the Phanerozoic cover of the northern margin of the African continent in southern Egypt and Libya. The term “metacraton” refers to a craton that has been remobilized during an orogenic event but is still recognizable dominantly through its rheological, geochronological and isotopic characteristics. Neoproterozoic remobilization of the Saharan Metacraton was in the forms of deformation, metamorphism, emplacement of igneous bodies, and probably local episodes of crust formation related to rifting and oceanic basin development. Relics of unaffected or only weakly remobilized old lithosphere are present as exemplified by the Archean to Paleoproterozoic charnockites and anorthosites of the Uweinat massif at the Sudanese/Egyptian/Libyan boarder. The article explains why the name “Saharan Metacraton” should be used, defines the boundaries of the metacraton, reviews geochronological and isotopic data as evidence for the presence of pre-Neoproterozoic continental crust, and discusses what happened to the Saharan Metacraton during the Neoproterozoic. A model combining collisional processes, lithospheric delamination, regional extension, and post-collisional dismembering by horizontal shearing is proposed.     

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.     

Extent,kinematics and tectonic origin of the Precambrian Aswa Shear Zone in eastern Africa

The Aswa Shear Zone (ASZ) is a fundamental Precambrian lithospheric structure that has been shaped by many tectonic events in eastern Africa. It separates the Saharan Metacraton in the northeast from the Northern Uganda terrane (which represents part of the Northeastern Congo block of the Congo craton) to the southwest. Nonetheless, its tectonic evolution is not fully understood. We used high-resolution airborne magnetic and radiometric data over Uganda integrated with Shuttle Radar Topography Mission (SRTM) Digital Elevation Model (DEM) in South Sudan to assess the extent, kinematics and contribute to the understanding of the tectonic origin of the ASZ. (1) Our results showed that the ASZ extends in a NW–SE direction for ~ 550 km in Uganda and South Sudan. (2) The airborne magnetic and radiometric data revealed a much wider (~ 50 km) deformation belt than the 5-10 km of the exposed surface expression of the ASZ. The deformation belt is defined by three NW-trending sinistral strike-slip shear zones bounding structural domains with magnetic fabric showing splays of secondary shear zones and shear-related folds. These folds are tighter close to the discrete shear zones with their axial traces becoming sub-parallel to the shear zones. A similar fold pattern is observed in South Sudan from the SRTM DEM. We interpreted these folds as due to ENE–WSW contraction associated with the sinistral strike-slip movement. (3) To the northeast, the magnetic patterns and radiometric signatures suggest the presence of a series of W-verging nappes indicative of strong E–W to NE–SW contraction deformation. (4) We relate the evolution of the ASZ to E–W to NE–SW Neoproterozoic oblique collision between East and West Gondwana. The deformation related to this collision was partitioned into E–W to NE–SW contraction resulting in W-verging thrusts in the east and a sinistral strike-slip movement along the NW-trending ASZ with the strain localized at the boundary between the Saharan Metacraton and the Northern Uganda terrane.     

Long wavelength gravity anomalies over India: Crustal and lithospheric structures and its flexure

Long wavelength gravity anomalies over India were obtained from terrestrial gravity data through two independent methods: (i) wavelength filtering and (ii) removing crustal effects. The gravity fields due to the lithospheric mantle obtained from two methods were quite comparable. The long wavelength gravity anomalies were interpreted in terms of variations in the depth of the lithosphere–asthenosphere boundary (LAB) and the Moho with appropriate densities, that are constrained from seismic results at certain points. Modeling of the long wavelength gravity anomaly along a N–S profile (77°E) suggest that the thickness of the lithosphere for a density contrast of 0.05 g/cm³ with the asthenosphere is maximum of ∼190 km along the Himalayan front that reduces to ∼155 km under the southern part of the Ganga and the Vindhyan basins increasing to ∼175 km south of the Satpura Mobile belt, reducing to ∼155–140 km under the Eastern Dharwar craton (EDC) and from there consistently decreasing south wards to ∼120 km under the southernmost part of India, known as Southern Granulite Terrain (SGT).The crustal model clearly shows three distinct terrains of different bulk densities, and thicknesses, north of the SMB under the Ganga and the Vindhyan basins, and south of it the Eastern Dharwar Craton (EDC) and the Southern Granulite Terrain (SGT) of bulk densities 2.87, 2.90 and 2.96 g/cm³, respectively. It is confirmed from the exposed rock types as the SGT is composed of high bulk density lower crustal rocks and mafic/ultramafic intrusives while the EDC represent typical granite/gneisses rocks and the basement under the Vindhyan and Ganga basins towards the north are composed of Bundelkhand granite massif of the lower density. The crustal thickness along this profile varies from ∼37–38 km under the EDC, increasing to ∼40–45 km under the SGT and ∼40–42 km under the northern part of the Ganga basin with a bulge up to ∼36 km under its southern part. Reduced lithospheric and crustal thicknesses under the Vindhyan and the Ganga basins are attributed to the lithospheric flexure of the Indian plate due to Himalaya. Crustal bulge due to lithospheric flexure is well reflected in isostatic Moho based on flexural model of average effective elastic thickness of ∼40 km. Lithospheric flexure causes high heat flow that is aided by large crustal scale fault system of mobile belts and their extensions northwards in this section, which may be responsible for lower crustal bulk density in the northern part. A low density and high thermal regime in north India north of the SMB compared to south India, however does not conform to the high S-wave velocity in the northern part and thus it is attributed to changes in composition between the northern and the southern parts indicating a reworked lithosphere. Some of the long wavelength gravity anomalies along the east and the west coasts of India are attributed to the intrusives that caused the breakup of India from Antarctica, and Africa, Madagascar and Seychelles along the east and the west coasts of India, respectively.     

Chronology and tectonic implications of Neoproterozoic blocks in the South Qinling Orogenic Belt,Central China

The Qinling Orogenic Belt (QOB) located between the North China Craton (NCC) and the Yangtze Craton (YZC) is composed of the North Qinling Belt (NQB), the South Qinling Belt (SQB) and the northern margin of the YZC. Detailed geological and geochronological investigations have revealed distinct Neoproterozoic blocks of various scales in the middle and western segments of the SQB, including the Madao block (MDB), Mihunzhen intrusion (MHI), Zhenggou block (ZGB), and Lengshuigou block (LSB) which constitute an east-west trending Neoproterozoic uplift zone of the basement continental blocks. These blocks are mainly composed of four lithological groups. Group #1 consists mainly of diorites in the LSB, the zircons from which yield a weighted mean ²⁰⁶Pb/ ²³⁸U age of ca. 941 Ma. Group #2 is chiefly composed of hornblende gabbros and diorites in the MHI and LSB, which were formed at ca. 885 Ma. Group #3 comprises massive diorites, quartz diorite, tonalites, granodiorites, and monzogranites in the MDB, MHI, ZGB and LSB, which were emplaced during ca. 785–740 Ma. Group #4 is composed of hornblende gabbros with an emplacement age of ca. 667 Ma in the ZGB.Detailed whole-rock geochemical and zircon Hf isotopic studies reveal the following: (1) The diorites of Group #1 were produced by partial melting of depleted mantle which was enriched by slab-derived melts, with the parental magmas contaminated by crustal materials. (2) The gabbros of Group #2 were derived from the partial melting of depleted mantle enriched by slab-derived melts and the diorites are the fractional crystallization products of the gabbroic magmas. (3) Group #3 which can be further sub-divided based on lithological assemblages and zircon Hf isotopic features into two subgroups, one representing massive diorites, quartz diorite, tonalites, granodiorites, and monzogranites (DTGMs) and the other composed of gneissic quartz diorites and granodiorites. Among these, the DTGMs were derived through magma mixing between melts derived from the depleted mantle wedge altered by slab-derived fluids and melts from juvenile sources, which subsequently underwent amphibole-dominated fractionation, whereas the gneissic granitoids formed through partial melting of thickened lower crust contaminated by depleted mantle melts. (4) The gabbros of Group #4 originated from a depleted lithospheric mantle that was enriched by slab-derived melts and fluids with contribution of asthenospheric mantle-derived materials. In conjunction with data from previous studies on the Neoproterozoic blocks in the SQB and basement blocks in the northern margin of the YZC, our new geological, geochronological and geochemical data suggest a large Neoproterozoic uplift zone in the SQB, which was destructed by Paleozoic to Mesozoic magmatism and deformation. The Neoproterozoic uplift zone of the SQB might have been separated from the northern margin of the YZC during the formation of the Mianlue Ocean, and might have evolved under an active continental margin setting and subsequent continental rift setting accompanied by significant crustal growth. The magmatism also resulted in the formation of important Neoproterozoic ore deposits and supplied the material sources for some of the major Mesozoic ore deposits.     

Metacraton: Nature,genesis and behavior

In this paper, we show with examples that cratons involved in intercontinental collisions in a lower plate position are often affected by orogenic events, leading to the transformation of their margins. In some cases, craton interiors can also be shaped by intense collisional processes, leading to the generation of intracratonic orogenic belts. We propose to call these events “metacratonization” and the resulting lithospheric tract “metacraton”. Metacratons can appear similar to typical orogenic belts (i.e. active margin transformed by collisional processes) but are actually sharply different. Their main distinctive characteristics (not all are present in each metacraton) are: (1) absence of pre-collisional events; (2) absence of lithospheric thickening, high-pressure metamorphism being generated by subduction, leading to high gradient in strain and metamorphic intensity; (3) preservation of allochthonous pre-collisional oceanic terranes; (4) abundant post-collisional magmatism associated with shear zones but not with lithospheric thickening; (5) presence of high-temperature–low-pressure metamorphism associated with post-collisional magmatism; (6) intracontinental orogenic belts unrelated to subduction and oceanic basin closures. Reactivation of the rigid but fractured metacratonic lithosphere will cause doming, asthenospheric volcanism emplacement, and mineralizations due to repetitive mineral enrichments. This paper provides several geological cases exemplifying these different metacratonic features in Scandinavia, Sahara, Central Africa and elsewhere. A special focus is given to the Saharan Metacraton because it is where the term “metacraton” originated and it is a vastly expanded tract of continental crust (5,000,000 km²). Metacratonization is a common process in the Earth's history. Considering the metacraton concept in geological studies is crucial for understanding the behavior of cratons and their partial destruction.     

He–Ar isotope geochemistry of iron and gold deposits reveals heterogeneous lithospheric destruction in the North China Craton

The North China Craton (NCC) provides a classic example for extensive destruction of the cratonic lithosphere. The Mesozoic magmatism which contributed to the decratonization of the NCC was also accompanied by the formation of a variety of mineral deposits. In order to gain further insights into the cratonic destruction process, typical iron and gold deposits are investigated here. Helium–argon isotopic data on pyrite, from typical skarn iron deposits of the Beiminghe and Fushan in the Han-Xing district of the central NCC, and the Linglong and Canzhuang gold deposits in the Jiaodong district in the eastern NCC, are presented in this paper. The ³He/⁴He, ⁴⁰Ar/³⁶Ar and ⁴⁰Ar/⁴He ratios show generally uniform patterns within the individual deposits and reveal a complex evolutionary history of the ore-forming fluids with varying degree of crust–mantle interaction. The ore-forming fluids associated with the gold mineralization at the Jiaodong mine have higher content of fluids of mantle origin with mantle helium ranging from 1.24% to 18.02% (average 6.73%; N = 18). In contrast, the ore-forming fluids related to the iron ore deposits contain less mantle contribution with mantle helium ranging from 0.12% to 4.96% (average 1.29%; N = 10). Our results suggest complex and heterogeneous crust–mantle processes associated with the magmatism and metallogeny, where the lithosphere of the eastern NCC was subjected to more extensive thinning and destruction as compared with that in the western part, consistent with the observations from geophysical studies in the region. Our study demonstrates that fluids associated with the Mesozoic metallogenic processes in the NCC provide useful insights into the geodynamics of destruction and refertilization of the cratonic lithosphere.     

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.     

Lithospheric structure and crust–mantle decoupling in the southeast edge of the Tibetan Plateau

Convergence between the Indian plate and the Eurasian plate has resulted in the uplift of the Tibetan Plateau, and understanding the associated dynamical processes requires investigation of the structures of the crust and the lithosphere of the Tibetan Plateau. Yunnan is located in the southwest edge of the plateau and adjacent to Myanmar to the west. Previous observations have confirmed that there is a sharp transition in mantle anisotropy in this area, as well as clockwise rotations of the surface velocity, surface strain, and fault orientation. We use S receiver functions from 54 permanent broad-band stations to investigate the structures of the crust and the lithosphere beneath Yunnan. The depth of the Moho is found to range from 36 to 40 km beneath southern Yunnan and from 55 to 60 km beneath northwestern Yunnan, with a dramatic variation across latitude 25–26°N. The depth of the lithosphere–asthenosphere boundary (LAB) ranges from 180 km to less than 70 km, also varying abruptly across latitude 25–26°N, which is consistent with the sudden change of the fast S-wave direction (from NW–SE to E–W across 26–28°N). In the north of the transition belt, the lithosphere is driven by asthenospheric flow from Tibet, and the crust and the upper mantle are mechanically coupled and moving southward. Because the northeastward movement of the crust in the Burma micro-plate is absorbed by the right-lateral Sagaing Fault, the crust in Yunnan keeps the original southward movement. However, in the south of the transition belt, the northeastward mantle flow from Myanmar and the southward mantle flow from Tibet interact and evolve into an eastward flow (by momentum conservation) as shown by the structure of the LAB. This resulting mantle flow has a direction different from that of the crustal movement. It is concluded that the Sagaing Fault causes the west boundary condition of the crust to be different from that of the lithospheric mantle, thus leading to crust–mantle decoupling in Yunnan.     

Upper mantle P-wave tomography across the Longmenshan fault belt from passive-source seismic observations along Aba-Longquanshan profile

A passive seismic experiment across the Longmenshan (LMS) fault belt had been conducted between August 2006 and July 2007 for the understanding of geodynamic process between the Eastern Tibet and Sichuan basin. We herein collected 3677 first P arrival times with high precision from seismograms of 288 teleseismic events so as to reconstruct the upper mantle velocity structure. Our results show that the depth of the Lithosphere–asthenosphere boundary (LAB) changes from 70 km beneath Eastern Tibet to about 110 km beneath Longquanshan, Sichuan Basin, which is consistent with the receiver function imaging results. The very thin mantle part of the lithosphere beneath Eastern Tibet may suggest the lithosphere delamination due to strong interaction between the Tibetan eastward escaping flow and the rigid resisting Sichuan basin, which can be further supported by the existences of two high-velocity anomalies beneath LAB in our imaging result. We also find there are two related low-velocity anomalies beneath the LMS fault belt, which may indicate magmatic upwelling from lithosphere delamination and account for the origin of tremendous energy needed by the devastating Wenchuan earthquake.     

Basaltic dykes of the Eastern Belt of Peninsular Malaysia: The effects of the difference in crustal thickness of Sibumasu and Indochina

Basaltic dykes of Peninsular Malaysia are confined to the Eastern Belt (Indochina/East Malaya block) as compared with the Western Belt (Sibumasu Block). The dyke intruded through a crustal fracture formed by stress developed from the evolution of two offshore basins (Malay and Penyu basins) east of Peninsular Malaysia. The Ar–Ar dating from the present study combined with the previous geochronological data indicate that the ages of dykes range from 79 ± 2 Ma to 179 ± 2 Ma. Thus it is difficult to correlate the dykes with the closure of Tethys during Permo-Triassic time because of the younger age of the dykes. The majority of the dykes exposed in the Eastern Belt may have been attributed to the difference of crustal thickness between the Eastern and Western belt of Peninsular Malaysia. A thicker Western Belt crust (13 km more than both Eastern and Central belts) is difficult to rupture with normal plate tectonic stress and therefore serves to contain the rise of a mantle derived melt. The chemistry indicates the basalts are olivine to quartz normative and are of the continental within-plate category.     

Density heterogeneity of the cratonic lithosphere: A case study of the Siberian Craton

Using free-board modeling, we examine a vertically-averaged mantle density beneath the Archean–Proterozoic Siberian Craton in the layer from the Moho down to base of the chemical boundary layer (CBL). Two models are tested: in Model 1 the base of the CBL coincides with the LAB, whereas in Model 2 the base of the CBL is at a 180 km depth. The uncertainty of density model is < 0.02 t/m³ or < 0.6% with respect to primitive mantle. The results, calculated at in situ and at room temperature (SPT) conditions, indicate a heterogeneous density structure of the Siberian lithospheric mantle with a strong correlation between mantle density variations and the tectonic setting. Three types of cratonic mantle are recognized from mantle density anomalies. ‘Pristine’ cratonic regions not sampled by kimberlites have the strongest depletion with density deficit of 1.8–3.0% (and SPT density of 3.29–3.33 t/m³ as compared to 3.39 t/m³ of primitive mantle). Cratonic mantle affected by magmatism (including the kimberlite provinces) has a typical density deficit of 1.0–1.5%, indicative of a metasomatic melt-enrichment. Intracratonic sedimentary basins have a high density mantle (3.38–3.40 t/m³ at SPT) which suggests, at least partial, eclogitization. Moderate density anomalies beneath the Tunguska Basin imply that the source of the Siberian LIP lies outside of the Craton. In situ mantle density is used to test the isopycnic condition of the Siberian Craton. Both CBL thickness models indicate significant lateral variations in the isopycnic state, correlated with mantle depletion and best achieved for the Anabar Shield region and other intracratonic domains with a strongly depleted mantle. A comparison of synthetic Mg# for the bulk lithospheric mantle calculated from density with Mg# from petrological studies of peridotite xenoliths from the Siberian kimberlites suggests that melt migration may produce local patches of metasomatic material in the overall depleted mantle.     

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.     

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.     

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.