Deep structure beneath the central part of the Balkan Peninsula was studied using P and S receiver function technique. Data from seismic stations from the Bulgarian National Seismological Network and several stations from neighbouring countries were used. Depth of Mohorovicic discontinuity has been estimated between 28–30 km in northern and central Bulgaria to 50 km in southwestern of Bulgaria. The 410 km mantle boundary is uplifted by 10 km relative to nominal depth in the area of Rhodopean Massif. In northern Bulgaria, the boundary is lowered by 10 km. Indications of a low-velocity layer are present at a depth exceeding 410 km. The thickness of the asthenosphere is estimated as 50 km and the depth of lithosphere-asthenosphere (LAB) boundary varies between 40 and 60 km. The results of this study have been published in Vinnik et. al., Izvestiya, Physics of the Solid Earth, 2021, Vol. 57, No. 6, pp. 849–863. This research has been carried out as part of a joint project supported by the National Science Foundation of Bulgaria (grant no. KP-06-RUSIA/27.09.2019) and the Russian Foundation for Basic Research (RFBR, grant no. 19-55-18008 Bolg_a).

By using P and S wave receiver functions and P and S wave travel time residuals, we have found velocity models for 16 seismograph stations in Eastern Anatolia. Our study is focused mainly on the mantle lithosphere, asthenosphere and transition zone. The volcanism and uplift of the Eastern Anatolia Plateau are thought to be related to the Bitlis slab break off and delamination of the continental lithosphere. Sinking cold slab and lithospheric drips can reduce temperature in the mantle transition zone (MTZ) by up to a few hundred degrees C. However, our analysis of seismic data provides no robust evidence of significant cooling of the transition zone. In the mantle immediately above the 410-km discontinuity there is a pronounced low S wave velocity layer that may be a source of the volcanism in the study region. Another low velocity layer is present at the base of the MTZ. The obtained S wave velocity models of the upper mantle can be divided into three groups. In the first group, the lithosphere—asthenosphere boundary (LAB) is at a depth of ~ 60 km. In the second group, the LAB is at a depth from 90 to 100 km. In the third group, the mantle lithosphere is practically absent. On a scale of our analysis there is no clear correspondence between the obtained mantle velocity models and the volcanism (< 23 Ma) exposed at the surface. Only the models of the first group are well represented in the neighboring Central Anatolian Plateau.

In spite of numerous active and passive seismological investigations, the existence of continuous or interrupted continental subduction below the Western Alps is still open to debate. Many of the observations focus on the Moho or the deeper part of the mantle, while reliable information on the Lithosphere‐Asthenosphere Boundary (LAB) below the Alpine region is scarce. Exploiting the data from the dense, broadband AlpArray Seismic Network we present a set of Receiver Function (RF) measurements on the Moho and LAB of a region encompassing the Western Alps, which includes the Ivrea Geophysical Body (IGB), a fragment of mantle placed at a few kilometers depth at the collision margin between Eurasia and Adria plates. We derive seismic velocity profiles of the crust‐uppermost mantle below each station down to about 250 km, through the joint inversion of P and S RF. We constrain the lateral variations of the Moho and LAB topographies across the colliding plates, and quantify the errors related to our measurements. Our results allow us to considerably expand the published data of the Moho depth and to add a unique set of new measurements of the LAB. Our observations show that Eurasia and Adria lithospheres have a comparable thickness (on average 90–100 km), and are colliding below the IGB, and that Eurasia is not presently continuously subducting below Adria. These observations suggest that there is a gap between the superficial (continental) European lithosphere and the deep (oceanic) lithosphere, confirming the discontinuous structure imaged by some seismic tomography models.

We investigate the upper mantle discontinuities in the central Mediterranean region by applying the P and S receiver function techniques on waveforms recorded at broadband stations located around the Tyrrhenian basin. P and S wave velocity profiles (down to 300-km depth) are calculated with joint inversion of P and S receiver functions. We could identify the Moho, lithosphere-asthenosphere boundary, and an underlying low-velocity layer between similar to 60- and similar to 200-km depth. The low-velocity layer is interpreted as asthenospheric material, and its lower boundary is identified below the western Ionian and Tyrrhenian basins as a sharp Lehmann discontinuity. Although the stations are located on different lithospheric domains we find a strong correlation between Moho and the lithosphere-asthenosphere boundary depths, which suggests ubiquitous coupling of the crust and lithospheric mantle, consistently with the southward opening of the Tyrrhenian basin. The Tyrrhenian and western Ionian basins present thinning of the transition zone of similar to 14 km, as inferred from a reduced P660s-P410s differential time. Below the southern Apennines we observe a standard differential time that implies an average mantle transition zone thickness. We explain these mantle transition zone thickness variations as due to temperature heterogeneity linked to the area's subduction history. Finally, under central Europe (the location of the deep S-to-P conversion points) two strong signals from nonstandard discontinuities within the mantle transition zone are observed. These signals can be explained as being generated at the boundaries of high seismic velocity layers that are spatially correlated with stagnant slabs in the transition zone detected by seismic tomography.

Summary Sharpness of the 410-km boundary is of interest because it is sensitive to water content in the transition zone. We evaluate the width of the 410-km discontinuity with a new seismic method. Our estimates are inferred from the amplitude ratio of the P2p410s and P410s seismic phases that are detected in P-wave receiver functions. We applied this method to seismic recordings from arrays of broad-band stations deployed in central Fennoscandia, southern Africa and southern China. The obtained estimates of width of the 410-km discontinuity range from 10 to 22 km and always exceed the width of 7 km which is expected for anhydrous conditions. The enlarged width may be interpreted in terms of hydrous conditions, but we have found only one region (the eastern Yangtze Craton in China) where the broad 410-km discontinuity, as expected, is accompanied by a broad transition zone. Water in the transition zone may be a kind of a global phenomenon, but evidence of the enlarged width of the transition zone may be missing in most of our data because the reference seismic model is affected by water, as well.

In this study, we use teleseismic P and S receiver functions (i.e., S-to-P and P-to-S converted signals) and their joint inversions to determine seismic discontinuities in the crust and upper mantle. Eight permanent broadband stations from the KOERI (Kandilli Observatory and Earthquake Research Institute) that are distributed along the eastern Pontides orogenic belt (EPOB), NE, Turkey, comprise our database. Inversion is performed by using a simulated annealing technique with and without travel time residuals. Our inversion results reveal the Moho depth, a high S velocity lid, a low-velocity zone, and the underlying upper mantle layer. The studied area is divided into two regions based on the station locations: (a) a northern region and (b) a southern region. The inversion results from the northern area produce crustal models that indicate that the uppermost crust is represented by a low P and S wave velocity (Vp = ~ 5.0 km/s and Vs = ~ 2.8 km/s). These velocities are clear evidence of (1) the sedimentary and volcanic rocks that widely crop out in the region and (2) a thinner uppermost crust, whereas the velocities of the southern region (Vp = ~ 6.0 km/s and Vs = ~ 3.1 km/s) indicate a thicker uppermost crust (~ 7 km). Our calculated Vp/Vs velocity ratio in the lower crust is approximately 1.90 and 1.80 in the northern and southern regions, respectively. These ratios are generally attributed to mafic rocks. Beneath the northern stations, the crustal thicknesses are 30, 33, 37, and 40 km from east to west, while the depths of the Moho are 46, 42, 39, and 44 km beneath the southern stations. Some velocity histograms show a transition from the high S velocity mantle lid to the low-velocity zone, which is known as the lithosphere–asthenosphere boundary. A representative value of the boundary’s depth is around 83 km in the north and ~ 88 km in the south. The hypothesis of some researchers regarding the existence of southward subduction beneath the eastern Pontides orogenic belt during the Late Mesozoic–Cenozoic is supported by our 2-D and 3-D velocity-depth models.

We investigate the mantle of central-eastern Greenland by using recordings with data from 24 local broad-band seismograph stations. We apply P wave receiver function technique and evaluate the difference in the arrival times of seismic phases that are formed by P to SV mode conversion at the 410-km and 660-km seismic discontinuities. These boundaries mark the top and bottom of the mantle transition zone (MTZ). The difference in the arrival time of the phases from the 410-km and 660-km discontinuities is sensitive to the thickness of the MTZ and relatively insensitive to volumetric velocity anomalies above the 410-km discontinuity. Near the east coast of Greenland in the region of the Skaergaard basalt intrusions we find two regions where the differential time is reduced by more than 2 s. The 410-km discontinuity in these regions is depressed by more than 20 km. The depression may be explained by a temperature elevation of similar to 150 degrees C. We hypothesize that the basaltic intrusions and the temperature anomalies at a depth of similar to 400 km are, at least partly, effects of the passage of Greenland over the Iceland hotspot at about 55 Ma. This explanation is consistent with the concept of tectosphere and implies that the upper mantle to a depth of similar to 400 km translates coherently with the Greenland plate.

We present models of the seismic structure of the crust and upper mantle in the interior of Greenland based on new seismological data from the TopoGreenland experiment. Until this experiment, all seismic data in Greenland was acquired close to the coast, where the crustal structure is affected by oceanic break-up. The TopoGreenland data acquisition programme in central-eastern Greenland included the first controlled source seismic experiment in interior Greenland and deployment of 24 broadband (BB) onshore stations for 3 years, partly on the ice cap. The 320 km long seismic refraction/wide-angle reflection profile was acquired on the ice cap by a team of six people during two-months in summer of 2011.

Receiver functions for the central Tien Shan and northern Tarim in central Asia reveal a pronounced depression on the 410 km discontinuity beneath the Permian basalts in Tarim. The depression may be caused by elevated temperature. The striking spatial correlation between the anomaly of the MTZ and the Permian basalts suggests that both may be effects of the same plume. This relation can be reconciled with the possible motion of Tarim on the order of 1000 km by assuming that the mantle layer, which has moved coherently with the plate since the Permian, extends to a depth of 410 km or more. Alternatively, the lithosphere and underlying mantle are decoupled at a depth of ∼ 200 km, but a cumulative effect of the Tarim plate motion since the Permian is less by an order of magnitude. A similar explanation is applicable to the Siberian traps.

Vp/Vs ratio where Vp and Vs are P- and S-wave velocities is an indicator of rock composition, but estimates of Vp/Vs for the lower continental crust remain sparse. We present estimates of Vs, Vp and Vp/Vs in the crust of the Archean-Paleoproterozoic Siberian craton that are obtained by simultaneous inversion of P and S receiver functions from GSN seismograph stations NRIL, YAK and TIXI. These stations are located in the region of the Siberian traps (NRIL), close to the Laptev Sea Rift (TIXI) and the Viluy rift system (YAK). The most conspicuous result of our analysis is a high Vp/Vs ratio (>= 2.0) at depths from 20-30 to 40 km. A very high Vp in this layer (from 7 to 8 km s(-1)) is indicative of magmatic underplating. We find broadly similar data in the western Mediterranean and in India. In a dry lower crust the Vp/Vs ratio is similar to 1.8, which is hard to reconcile with the estimates >2.0. A coincidence in depths of zones of high electric conductivity and of anomalously high Vp/Vs in Siberia suggests that both may have the same origin: fluid-filled porosity. The porosity which is required by our seismic observations is on the order of 1 per cent. Origins of the fluids may be linked with processes of magmatic underplating.

We investigate structure of the mantle transition zone (MTZ) under the central Tien Shan in central Asia by using recordings of seismograph stations in Kyrgyzstan, Kazakhstan and adjacent northern China. We apply P-wave receiver functions techniques and evaluate the differential time between the arrivals of seismic phases that are formed by P to SV mode conversion at the 410-km and 660-km seismic boundaries. The differential time is sensitive to the thickness of the MTZ and insensitive to volumetric velocity anomalies above the 410-km boundary. Under part of the southern central Tien Shan with the lowest S wave velocity in the uppermost mantle and the largest thickness of the crust, the thickness of the MTZ increases by 15–20km relative to the ambient mantle and the reference model IASP91. The increased thickness is a likely effect of low (about −150K) temperature. This anomaly is indicative of delamination and sinking of the mantle lithosphere. The low temperature in the MTZ might also be a relic of subduction of the oceanic lithosphere in the Paleozoic, but this scenario requires strong coupling and coherence between structures in the MTZ and in the lithosphere during plate motions in the last 300Myr. Our data reveal a reduction of thickness of the MTZ of 10–15km under the Fergana basin, in the neighborhood of the region of small-scale basaltic volcanism at the time near the Cretaceous-Paleogene boundary. The reduced thickness of the MTZ is the effect of a depressed 410-km discontinuity, similar to that found in many hotspots. This depression suggests a positive temperature anomaly of about 100–150K, consistent with the presence of a thermal mantle plume. A similar depression on the 410-km discontinuity is found underneath the Tarim basin.

P- and S-wave velocity profiles (Vp and Vs) from the Earth's surface to a depth of 300km are obtained by simultaneous inversion of P- and S-wave receiver functions for the temporary POLENET/LAPNET array in northern Finland and 5 permanent stations in southern Finland. The obtained Vp/Vs velocity ratio in the uppermost mantle is anomalously low (1.65–1.70 versus the standard 1.8). This ratio can be explained by a high (~30%) fraction of orthopyroxene in the depleted upper mantle. An increase of ~4% in the Vs values is detected at a depth of 110–130km. Under southern Finland the high-velocity mantle keel may extend beyond the depth of 300km, but under northern Finland we detect a low-Vs zone (LVZ) with the top at a depth of 160km. This depth corresponds to the intersection of the geotherm with the wet solidus of peridotite at a temperature near ~1100°C, and suggests partial melting as an explanation for the origin of the LVZ. The bottom of the LVZ (the Lehmann discontinuity) is found at a depth of 240km. The fast direction of seismic azimuthal anisotropy beneath northern Finland in the depth interval from 200–240km to 320km is parallel to the current APM direction (60°) but in the interval from 160km to 200–240km the fast direction (150°) appears to be normal to the APM direction (Vinnik et al., 2014). The present study indicates that the fast direction of 150° confines to the LVZ where the flow in the direction normal to the APM direction is unlikely. More likely, the flow direction is parallel to the APM direction and normal to the fast direction of anisotropy, as suggested by experiments with sheared peridotite-type rocks containing melt.