Critical evaluation and geological correlation of teleseismic phase data from Indian seismological stations

We evaluated the quality of seismic phase data from Indian seismological stations through the analysis of teleseismic travel times reported during 1976–83 and infer that only WWSSN stations (NDI, SHL, POO, KOD) apart from GBA and HYB can be rated satisfactory while the majority of stations (more than 40) produce very poor quality data sets. Detailed analysis of teleseismic P-wave travel time residuals shows that while the average structure of the upper mantle beneath India has high velocity (negative residuals) there are marked lateral variations. In particular, three zones of anomalous positive residuals (low velocity) are observed: one beneath the north western part of the Deccan trap, the second covering the southernmost peninsula (granulite terrain) and a third rather localized one, to the north of Delhi coinciding with Delhi-Haridwar ridge. New Delhi exhibits strong negative residuals in the E-SE quadrant along with negative station anomaly, implying that it is underlain by an anomalous high velocity crust/upper mantle. The negative residuals observed over India, continue beneath the Himalaya till the south of Lhasa but change sign further northward, suggesting the northern limit of the Indian upper mantle structure.     

Three-dimensional intraplate stress distributions associated with topography and crustal density inhomogeneities beneath the Deccan Volcanic Province

The generation of intraplate earthquakes has been attributed to perturbations in the stress regime, either due to surface and sub-surface loading or strength weakening of the rock mass. The present work aims at estimating the intraplate stresses associated with topography and crustal density inhomogeneities beneath the Deccan Volcanic Province (DVP). A layered crustal model with irregular interfaces of small amplitude has been used for elastostatic stress calculations. The computed principal stress differences show a significant concentration at 5–20 km depths beneath the western side of the region. The maximum magnitude of principal stress difference occurs beneath the Karad at a depth of 10 km with a value of 60 MPa. The deviatoric stress estimates are further superposed on inferred stresses due to the regional plate tectonic forces. These results show principal stress difference concentrations beneath the Koyna, Poona and Karad regions which may thus be more vulnerable to brittle failure. It is also seen that the principal total stress directions point to the strike slip motion at Koyna, similar to that which is associated with the 1967 Koyna earthquake.     

The azimuthal variation of teleseismic P-residuals in southern California: Implications for upper-mantle structure

The azimuthal variation of teleseismic P-delays has been investigated for stations of the USGS-Caltech Southern California Seismographic Network. Normalized residuals show azimuthal variations as large as 1.2 s, and must be explained in terms of upper mantle structure. The observed azimuthal dependence implies the presence of a region of depressed velocity beneath the Imperial Valley, and regions of increased velocity below the Sierra Nevada, southwest Arizona, and much of the Transverse Ranges. The last is a major high velocity ridge-like structure, extending from a depth of ~40 km to over 100 km, which crosses, but is not offset by, the San Andreas Fault. This suggests that the plate boundary at depth may diverge from its surface expression. The horizontal shear resulting from the divergence of crust and mantle plate boundaries may be accommodated by a zone of decoupling associated with the regionally observed 7.8 km/s (P_n) layer.     

Apparent velocity and azimuths of short period P-Waves recorded at indian WWSSN stations network

Over 200 earthquakes in the distance range 30°–90° and azimuthal range 0°–360°, recorded at Indian WWSSN stations, have been used in the present study. We have treated the four WWSSN recording stations i.e. New Delhi, Poona, Shillong and Kodai-Kanal, as Super Large Aperture Seismic Array (SLASA) Network with New Delhi being its cross-over point. Short period P-wave data as obtained from these stations have been analysed using a least square technique. Slowness and azimuthal anomalies have been computed for all these events. Relative time residuals have also been calculated. A velocity model has been derived on the basis of the slowness and travel-time data. The results do not indicate presence of any triplication in the travel-time curve. Variations in the relative residuals refer to the tectonic features beneath the recording stations. The P-wave velocity increases continuously in the lower mantle region and there is no indication for the presence of any appreciable velocity gradient.     

Basement configuration of the Jhagadia–Rajpipla profile in the western part of Deccan syneclise,India from travel-time inversion of seismic refraction and wide-angle reflection data

2-D velocity structure up to the basement is derived by travel-time inversion of the first arrival seismic refraction and wide-angle reflection data along the SW–NE trending Jhagadia–Rajpipla profile, located on the western part of Deccan syneclise in the Narmada–Tapti region. The study region is mostly covered by alluvium. Inversion of refraction and wide-angle reflection data reveals four layered velocity structure above the basement. The first two layers with P-wave velocities of 1.95–2.3 km s^?1 and 2.7–3.05 km s^?1 represent the Recent and Quaternary sediments respectively. The thickness of these sediments varies from 0.15 km to 3.4 km. The third layer with a P-wave velocity of 4.8–5.1 km s^?1 corresponds to the Deccan volcanics, whose thickness varies from 0.5 km to 1.0 km. Presence of a low velocity zone (LVZ) below the high velocity volcanic rocks in the study area is inferred from the travel-time ‘skip’ and amplitude decay of the first arrival refraction data and the wide-angle reflection from top of the LVZ present immediately after the first arrival refraction from Deccan Trap layer. The thickness of the low velocity Mesozoic sediments varies from 0.3 km to 1.7 km. The basement with a P-wave velocity of 5.9–6.15 km s^?1 lies at a depth of 4.9 km near Jhagadia and shallows to 1.2 km towards northeast near Rajpipla. The results indicate presence of low velocity Mesozoic sediments hidden below the Deccan Trap layer in the western part of the Deccan syneclise.     

Anomalous granulite crust of South India — signatures from converted teleseismic waves

Analysis of teleseismic waves usingS-P converted phases, travel time-terms and residual travel times point to the presence of an anomalous thick (4—5 km thicker) low velocity (−3%) crust beneath Kodaikanal (KOD) on granulites characterized by an oriented inhomogeneity inferred as possibly due to Mylonites in contrast to the nature of crust beneath the adjoining precambrian granite-gneiss terrain. The observed seismic signatures in the South India granulites (represented by KOD) offer an opportunity to discriminate between the competing hypotheses of tectonic thickening and magmatic underplating to explain the origin of the granulites of South India. This analysis lends support to the hypothesis of a continent-continent collision origin for the granulites in the study region.     

Structure of the Earth’s crust in the northern part of the Barents–Kara region along the 4-AR DSS profile

The 1370 km long 4-AR reference profile crosses the North Barents Basin, the northern end of the Novaya Zemlya Rise, and the North Kara Basin. Integrated geophysical studies including common deep point (CDP) survey and deep seismic sounding (DSS) were carried out along the profiles. The DSS was performed using autonomous bottom seismic stations (ABSS) spaced 10–20 km apart and a powerful air gun producing seismic signals with a step size of 250 m. As a result, detailed P- and S-wave velocity structures of the crust and upper mantle were studied. The basic method was ray-tracing modeling. The Earth’s crust along the entire profile is typically continental with compressional wave velocities of 5.8–7.2 km/s in the consolidated part. Crustal thickness increases from 30 km near the islands of Franz Josef Land to 35 km beneath the North Barents Basin, 50 km beneath the Novaya Zemlya Rise, and 40 km beneath the North Kara Basin. The North Barents Basin 15 km deep is characterized by unusually low velocities in the consolidated crust: The upper crust layer with velocities of 5.8–6.4 km/s has a thickness of about 15 km beneath the basin (usually, this layer wedges beneath deep sedimentary basins). Another special property of the crust in the North Barents Basin is the destroyed structure of the Moho.     

Relative teleseismic travel-time residuals, North Island, New Zealand, and their relation to upper-mantle structure

Relative travel-time residuals computed from clear P-wave arrivals at fourteen seismograph stations in the North Island, New Zealand, from five deep-focus events in the Banda Sea region, show large spatial variations of up to 3 sec. The variations can be explained by higher than normal velocities in the oceanic lithosphere which is underthrust to depths of 350 km beneath the North Island. After correction for crustal structure, the residuals imply an average P-wave velocity about 11% higher than in the surrounding mantle. The lack of suitable source events at azimuths other than northwest prevents a more detailed investigation by this means.     

Structure of the crust and uppermost mantle of the Yanshan Belt and adjacent regions at the northeastern boundary of the North China Craton from Rayleigh Wave Dispersion Analysis

We have used Rayleigh wave dispersion analysis and inversion to produce a high resolution S-wave velocity imaging profile of the crust and uppermost mantle structure beneath the northeastern boundary regions of the North China Craton (NCC). Using waveform data from 45 broadband NCISP stations, Rayleigh wave phase velocities were measured at periods from 10 to 48 s and utilized in subsequent inversions to solve for the S-wave velocity structure from 15 km down to 120 km depth. The inverted lower crust and uppermost mantle velocities, about 3.75 km/s and 4.3 km/s on average, are low compared with the global average. The Moho was constrained in the depth range of 30–40 km, indicating a typical crustal thickness along the profile. However, a thin lithosphere of no more than 100 km was imaged under a large part of the profile, decreasing to only ~ 60 km under the Inner Mongolian Axis (IMA) where an abnormally slow anomaly was observed below 60 km depth. The overall structural features of the study region resemble those of typical continental rift zones and are probably associated with the lithospheric reactivation and tectonic extension widespread in the eastern NCC during Mesozoic–Cenozoic time. Distinctly high velocities, up to ~ 4.6 km/s, were found immediately to the south of the IMA beneath the northern Yanshan Belt (YSB), extending down to > 100-km depth. The anomalous velocities are interpreted as the cratonic lithospheric lid of the region, which may have not been affected by the Mesozoic–Cenozoic deformation process as strongly as other regions in the eastern NCC. Based on our S-wave velocity structural image and other geophysical observations, we propose a possible lithosphere–asthenosphere interaction scenario at the northeastern boundary of the NCC. We speculate that significant undulations of the base of the lithosphere, which might have resulted from the uneven Mesozoic–Cenozoic lithospheric thinning, may induce mantle flows concentrating beneath the weak IMA zone. The relatively thick lithospheric lid in the northern YSB may serve as a tectonic barrier separating the on-craton and off-craton regions into different upper mantle convection systems at the present time.     

The seismic anisotropy beneath Southern Sakhalin: Evidence from the S-wave parameters of deep-focus earthquakes

The measurements of the parameters of split shear (S) waves from local deep-focus earthquakes recorded in 2005–2007 by a network of 12 seismic stations in Southern Sakhalin are presented. The results revealed the heterogeneous distribution of the anisotropic properties beneath Southern Sakhalin. The azimuths of the fast S-wave polarization beneath the stations in the central part of the peninsula are oriented along the NNW and NNE-NE directions normal to and along the Kuril Trench. Beneath the stations located along the western and eastern coasts, the azimuths of the fast S-wave polarization change their direction from NNW in the northern area to E-SE in the southern area. The highest anisotropy degree (up to 0.9–1.5%) is recorded beneath the central part of Southern Sakhalin. The maximum values of the discrepancy in the arrival time of the split S-waves are observed when the azimuth of the fast S-wave is oriented along the NNE beneath the active fault zones. The analysis of the variations of the S-wave lag time shows their weak depth dependence. The highest anisotropy is assumed in the upper layers of the medium (down to a depth of about 250 km). The results obtained for the dominating wave frequency of 1–5 Hz represent mainly the medium-scale anisotropy of the top of the studied region.     

Tomographic studies of the upper mantle beneath the Italian region*

Teleseismic P arrivals at seismological stations are inverted into a model of velocity perturbations down to a depth of about 470 km. Directionally independent average residuals, computed from steeply inciding waves, are transformed into a model of lithospheric thickness. Both models show a good correspondence with the main tectonic features of the Italian Peninsula. Positive velocity perturbations are observed beneath the Alps and in depths over 200 km also beneath the Po Basin. A high-velocity anomaly of the Tyrrhenian subduction is less pronounced, probably due to a directional dependence of P velocities in the mantle. Negative velocity perturbations indicate several low-velocity regions, e.g. beneath the Northern Apennines, the Sicily region and in the upper 100 km beneath the Po Basin. The amplitudes of velocity perturbations beneath the depth of 200 km are smaller on the average than those in the upper two layers. The whole region is characterized by large undulations of the lithosphere base which reaches depths from less than 60 km to more than 150 km. The most prominent lithospheric root beneath the Alps is a product of the collision between the European and the Adriatic plates while the lithospheric thickening beneath the Calabrian coast is likely to be connected with the eastern wing of the Tyrrhenian subduction. The dramatic changes of lithosphere thickness between the northern and the southern Apenninic arcs and northern Calabria as well as the thinnings at the western closure of the Po Basin, indicate important deep-seated boundaries of lithospheric blocks of autonomous geodynamic development.     

Geochemical investigation of a semi-continuous extrusive basaltic section from the Deccan Volcanic Province, India: implications for the mantle and magma chamber processes

Spatial and temporal variations in the geochemistry of an extrusive basaltic section of Deccan traps record progressive changes in mantle melting and crustal filtration and are relevant to understand continental flood basalt (CFB) magmatism. In the present work we have carried out detailed field, petrographic, density and magnetic susceptibility, and geochemical investigations on a small, semi-continuous extrusive section in the eastern Deccan Volcanic Province (DVP) to understand the role of shallow magma chambers in CFB magmatism. Four formations, Ajanta, Chikhli, Buldhana and Karanja crop out in the Gangakhed–Ambajogai area with increasing elevation. Our studies indicate that: (1) the Karanja Formation represents a major magma addition, as indicated by abrupt change in texture, increases in MgO, CaO, Ni, Cr, and Sr, and drastic decreases in Al₂O₃, Na₂O, K₂O, Rb, Ba, REE, bulk-rock density and magnetic susceptibility; (2) assimilation fractional crystallization, crystal-laden magmas, and accessory cumulus phases influence the trace element chemistry of Deccan basalts; (3) the predicted cumulate sequence of olivine gabbro–leucogabbro–oxide-apatite gabbro is supported by the observed layered series in a shallow magma chamber within the DVP; (4) the initial magma was saturated with olivine, plagioclase, and augite, and final the pressure of equilibration for the Gangakhed–Ambajogai section basalts is ~2 kbar (~6 km depth); (5) petrophysical parameters act as proxies for magmatic processes; (6) a small layer of oxide-rich basalts may represent the latest erupted pulse in a given magmatic cycle in the DVP; (7) parental basalts to some of the red boles, considered as formation boundaries, might represent small degree partial melts of the mantle; (8) SW Deccan basaltic-types continue into the eastern DVP; and (9) in addition to the magma chamber processes, dynamic melting of the mantle may have controlled DVP geochemistry. The present study underscores the importance of mapping specific stratigraphic intervals in limited areas to understand mantle and magma chamber processes relevant to CFB magmatism.     

Shear wave velocity structure beneath the Archaean granites around Hyderabad, inferred from receiver function analysis

Broadband receiver functions abstracted from teleseismicP waveforms recorded by a 3-component Streckeisen seismograph at Hyderabad, have been inverted to constrain the shear velocity structure of the underlying crust. Receiver functions obtained from the Hyderabad records of both shallow and intermediate focus earthquakes lying in different station-event azimuths, show a remarkable coherence in arrival times and shapes of the significant shear wave phases:Ps, PpPs, PsPs/PpSs, indicating horizontal stratification within the limits of resolution. This is also supported by the relatively small observed amplitudes of the tangential component receiver functions which are less than 10% of the corresponding radial component. Results of several hundred inversions of stacked receiver functions from closely clustered events (within 2°), show that the crust beneath the Hyderabad granites has a thickness of 36 ± 1 km, consisting of a 10 km thick top layer in which shear wave velocity is 3.54 ± 0.07 km/sec, underlain by a 26 ± 1 km thick lower crust in which the shear wave velocity varies uniformly with a small gradient of 0.02 km/sec/km. The shear wave velocity at its base is 4.1 ± 0.05 km/sec, just above the moho transition zone which is constrained to be less than 4 km thick, overlying a 4.74 ±0.1 km/sec half space.     

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.     

Deep structure of the baikal and other continental rift zones from seismic data

Abyssal variations beneath the Baikal rift zone are revealed in an irregular seismic stratification of the crust, the presence of an intracrust waveguide and by the vast (> 200,000 km²) underlying area of anomalously low velocity (P_n = 7.6−7.8 km/sec) uppermost mantle. In its abyssal structure the Baikal rift is heterogeneous along the strike, with sharp changes in crustal thickness (35–50 km).Comparison of first-arrival seismic-velocity curves and also the respective velocity columns reveals the essential similarity of upper-mantle seismic cross-sections for all continental rift zones. The anomalous upper layer of the mantle (ca. 7.7 km/sec) is relatively thin (15-13 km) and can be linked with the mantle waveguide only locally.     

2D model of the crust and uppermost mantle along rift profile, Siberian craton

A two-dimensional model of the crust and uppermost mantle for the western Siberian craton and the adjoining areas of the Pur-Gedan basin to the north and Baikal Rift zone to the south is determined from travel time data from recordings of 30 chemical explosions and three nuclear explosions along the RIFT deep seismic sounding profile. This velocity model shows strong lateral variations in the crust and sub-Moho structure both within the craton and between the craton and the surrounding region. The Pur-Gedan basin has a 15-km thick, low-velocity sediment layer overlying a 25-km thick, high-velocity crystalline crustal layer. A paleo-rift zone with a graben-like structure in the basement and a high-velocity crustal intrusion or mantle upward exists beneath the southern part of the Pur-Gedan basin. The sedimentary layer is thin or non-existent and there is a velocity reversal in the upper crust beneath the Yenisey Zone. The Siberian craton has nearly uniform crustal thickness of 40–43 km but the average velocity in the lower crust in the north is higher (6.8–6.9 km/s) than in the south (6.6 km/s). The crust beneath the Baikal Rift zone is 35 km thick and has an average crustal velocity similar to that observed beneath the southern part of craton. The uppermost mantle velocity varies from 8.0 to 8.1 km/s beneath the young West Siberian platform and Baikal Rift zone to 8.1–8.5 km/s beneath the Siberian craton. Anomalous high Pn velocities (8.4–8.5 km/s) are observed beneath the western Tunguss basin in the northern part of the craton and beneath the southern part of the Siberian craton, but lower Pn velocities (8.1 km/s) are observed beneath the Low Angara basin in the central part of the craton. At about 100 km depth beneath the craton, there is a velocity inversion with a strong reflecting interface at its base. Some reflectors are also distinguished within the upper mantle at depth between 230 and 350 km.     

Crust–upper mantle seismic velocity structure across Southeastern China

One in-line wide-angle seismic profile was conducted in 1990 in the course of the Southeastern China Continental Dynamics project aimed at the study of the contact between the Cathaysia block and the Yangtze block. This 380-km-long profile extended in NW–SE direction from Tunxi, Anhui Province, to Wenzhou, Zhejiang Province. Five in-line shots were fired and recorded at seismic stations with spacing of about 3 km along the recording line. We have used two-dimensional ray tracing to model P- and S-wave arrivals and provide constraints on the velocity structure of the upper crust, middle crust, lower crust, Moho discontinuity, and the top part of the lithospheric mantle. P-wave velocity, S-wave velocity and V_P/V_S ratio are mapped. The crust is 36-km thick on average, albeit it gradually thins from the northwest end to the southeast end (offshore) of the profile. The average crustal velocity is 6.26 km/s for P-waves but 3.6 km/s for S-waves. A relatively narrow low-velocity layer of about 4 km of thickness, with P- and S-wave velocities of 6.2 km/s and 3.5 km/s, respectively, marks the bottom of the middle crust at a depth of 23-km northwest and 17-km southeast. At the crust–mantle transition, the P- and S-wave velocity change quickly from 7.4 to 7.8 km/s (northwest) and 8.0 to 8.2 km/s (southeast) and from 3.9 to 4.2 km/s (northwest) and 3.9 to 4.5 km/s (southeast), respectively. This result implies a lateral contrast in the upper mantle velocity along the 140 km sampled by the profile approximately. The average V_P/V_S ratio ranges from 1.68–1.8 for the upper crust to 1.75 for the middle and 1.75–1.85 for lower crust. With the interpretation of the wide-angle seismic data, Jiangshan–Shaoxin fault is considered as the boundary between the Yangtze and the Cathaysia block.     

Crustal velocity and Vp/Vs structures beneath central Anatolia from local seismic tomography

We applied a seismic tomography technique to arrival time data generated by local crustal earthquakes in central Anatolia in order to study the three-dimensional velocity and Vp/Vs structures and their relation with the complex tectonic processes and seismic activity occurring in the study region. The relatively equal and large number of both P- and S-wave arrival times comprising a total of 51,650 arrivals and the relatively uniform distribution of the recording stations imply that the obtained velocity anomalies are reliable features down to a depth of 40 km. This is also evident from the results of the checkerboard resolution test, hit count, and the ray-path coverage. The inversion results indicate the existence of strong lateral heterogeneities in the crust and uppermost mantle beneath central Anatolia. Prominent low-velocity anomalies are clearly imaged at all layers especially beneath existing volcanoes and the active fault segments. Higher-than-average Vp/Vs ratios are widely distributed, indicating the possible existence of over-pressurized fluids that may be responsible for the triggering of the large crustal earthquakes along the north and east Anatolian fault zones. We noticed that the seismic activity occurs mainly at the low-velocity areas and to a lesser extent in some high-velocity zones, perhaps because of the complex tectonics and geological structures. These observations imply that all the zones with velocity anomalies—either low or high—are potential sites for strain energy accumulation and subsequent release. The obtained velocity and Vp/Vs models are consistent with previous geophysical measurements conducted beneath central Anatolia and give much deeper understanding of the current seismotectonic processes occurring in the region.     

Crustal accretion beneath the Koyna coastal region (India) and Late Cretaceous geodynamics

In 1967 a major earthquake in the Koyna region attracted attention to the hitherto considered stable Indian shield. The region is covered by a thick pile of Deccan lava flows and characterized by several hidden tectonic features and complex geophysical signatures. Although deep seismic sounding studies have provided vital information regarding the crustal structure of the Koyna region, much remains unknown. The two available DSS profiles in the region have been combined along the trend of Bouguer gravity anomalies. Unified 2-D density modelling of the Koyna crust/mantle suggests a ca. 3 km thick and 40 km wide high velocity/high density anomalous layer at the base of the crust along the coastline. The thickness of this anomalous layer decreases gradually towards the east and ahead of the Koyna gravity low the layer ceases to be visible. Based on the seismic and gravity data interpretation in the geodynamical/rheological boundary conditions the anomalous layer is attributed to igneous crustal accretion at the base of the crust. It is suggested that the underplated layer is the imprint of the magmatism caused by the deep mantle plume when the northward migrating Indian plate passed over the Reunion hotspot.     

Tomographic Imaging of the India-Asia Plate Collisional Tectonics and Mantle Upwelling Beneath Western Tibet

To better understand the lithosphere mantle collision tectonics between the India plate and Asia plate, we determine three dimensional P wave velocity structure beneath western Tibet using 27,439 arrival times from 2,174 teleseismic events recorded by 182 stations of Hi-CLIMB Project and 16 stations in the north of Hi-CLMB. Our tomographic images show the velocity structure significantly difference beneath northern and southern Qiangtang, which can further prove that the Longmu Co-Shuanghu ophiolitic belt is a significant tectonic boundary fault zone. There are two prominent high velocity anomalies and two prominent low velocity anomalies in our images. One obvious high velocity anomalies subduct beneath the Tibet at the long distance near 34°N, whereas it is broke off by an obvious low velocity anomaly under the IYS. We interpret them as northward subducting Indian lithosphere mantle and the low velocity anomanly under IYS likely reflects mantle material upwelling triggered by tearing of the northward subduction Indian lithosphere. The other prominent high velocity anomaly was imaged at a depth from 50 km to 200 km horizontal and up to the northern Qiangtang with its southern edge extending to about 34°N through Hoh Xil block. We infer it as the southward subducting Asia lithosphere mantle. The other widely low velocity anomaly beneath the Qiangtang block lies in the gap between the frontier of India plate and Asia plate, where is the channel of mantle material upwelling.