Shear-wave velocity structure of Africa from Rayleigh-wave analysis

The 3D S-velocity structure beneath Africa is shown by means of a 2D S-velocity mapping for depths raging from zero to 500 km, determined by the regionalization and inversion of Rayleigh-wave dispersion. The traces of 94 earthquakes, occurred from 1990 to 2009 in the study area, have been used to obtain the Rayleigh-wave dispersion. These earthquakes were registered by 61 seismic stations located on Africa and the surrounding area. The dispersion curves were obtained for periods between 5 and 300 s, by digital filtering with a combination of MFT and TVF filtering techniques. After that, all seismic events (and some stations) were grouped to obtain a dispersion curve for each source-station path. These dispersion curves were regionalized and after inverted according to generalized inversion theory, to obtain shear-wave velocity models for rectangular blocks with a size of 5° × 5°. The 3D S-velocity structure obtained through this procedure is shown in the 2D S-velocity maps plotted for several depths. These results agree well with the geology and other geophysical results previously obtained. The obtained S-velocity models suggest the existence of lateral and vertical heterogeneity. The zones with consolidated and old structures (as cratons) present greater S-velocity values than the other younger zones. Nevertheless, in the depth range from 20 to 40 km, the different Moho depths present in the study area generate the principal variation of S-velocity. A similar behaviour is found for the depth range from 60 to 230 km, in which the lithosphere–asthenosphere boundary generates the principal variations of S-velocity. Finally, it should be highlighted a new and interesting feature obtained in this study: the definition of the base of the asthenosphere, for depths ranging from 160 to 280 km, in the whole African continent.     

Shear velocity and Vp/Vs ratio structure of the crust beneath the southern margin of South China continent

We here present the results of the inverse modeling of crustal S-phases recorded from a 400-km-long seismic profile, with azimuth nearly N30W, from Lianxian, near Hunan Province, to Gangkou Island, near Guangzhou City, Guangdong Province, in the southern margin of South China continent. The finding in this case is that many shot gathers provided by this wide-angle seismic experiment show relatively strong reflected and refracted S-phases, in particular some crustal refractions (Sg waves) and Moho reflections (SmS waves or simply Sm waves). The P-wave velocity structure of the crust and uppermost mantle was already obtained through the interpretation of vertical-component shot gathers. Now, with constraints introduced by the P-wave velocity architecture and after picking up S-wave traveltime data on the seismograms, we have obtained the S-velocity model of the crust by adjusting these traveltimes but keeping the geometry of the crustal reflectors. Our results demonstrate: (1) the average crustal S-velocity is about 3.64 km/s to the northwest of the Wuchuan-Sihui fault, and 3.62 km/s to the southeast of this fault; (2) relatively constant S-velocity of about 3.42 km/s for the upper crust, 3.55 km/s for the middle crust and laterally varying shear velocity around 3.82 km/s for the lower crust; (3) correspondingly, Vp/Vs ratio is 1.73 for the upper crust, 1.71 for the middle crust and 1.74 for the lower crust. Both shear velocities and Vp/Vs ratio correlate well with the major active faults that break the study area, and show significant changes especially in the upper crust. High Poisson’s ratio (1.8) is observed at shallow depth beneath the Minzhong depression to the southeast of the Wuchuan-Sihui fault and the Huiyuan depression in the southern margin of South China continent. In contrast, a very low Vp/Vs ratio (1.68) is observed between 8 and 14 km depth beneath Huiyuan. At deeper depth, a high Vp/Vs ratio (1.76) is observed in the lower crust beneath the Minzhong depression.     

Shear-wave velocity structure of the south-eastern part of the Iberian Peninsula from Rayleigh wave analysis

In this study, we present the lithospheric structure of the south-eastern part of the Iberian Peninsula by means of a set of 2D images of shear velocity, for depths ranging from 0 to 50 km. This goal will be attained by means of the inversion of the Rayleigh wave dispersion. For it, the traces of 25 earthquakes occurred on the neighbouring of the study area, from 2001 to 2003, will be considered. These earthquakes have been registered by 11 broadband stations located on Iberia. All seismic events have been grouped in source zones to get an average dispersion curve for each source-station path. The dispersion curves have been measured for periods between 2 and 45 s, by combination of two digital filtering techniques: Multiple Filter Technique and Time Variable Filtering. The resulting set of source-station averaged dispersion curves has been inverted according to the generalized inversion theory, to get S-wave velocity models for each source-station path. Later, these models have been interpolated using the method of kriging, to obtain a 2D mapping of the S-wave velocity structure for the south-eastern part of Iberia. The results presented in this paper show that the techniques used here are a powerful tool to investigate the crust and upper mantle structure, through the dispersion analysis and its inversion to obtain shear velocity distributions with depth. By means of this analysis, principal structural features of the south-eastern part of Iberia, such as the existence of lateral and vertical heterogeneity in the whole study area, or the location of the Moho discontinuity at 30 km of depth (with an average S-velocity of uppermost mantle of 4.7 km/s), have been revealed. Other important structural features revealed by this analysis have been that the uppermost of Iberian massif shows higher velocity values than the uppermost of the Alpine domain, indicating that the massif is old and tectonically stable. The average velocity of the crust in Betic cordillera is of 3.5 km/s, while in the Iberian massif is 3.7 km/s. All these features are in agreement with the geology and other previous geophysical studies.     

Crust and upper mantle structure and its tectonic implications in the South China Sea and adjacent regions

We present a 3D S-velocity model for the crust and upper mantle of the South China Sea and the surrounding regions, constrained from the analysis of over 12,000 of fundamental Rayleigh wave dispersion curves between 10 s and 150 s periods. The lateral resolution was found to vary from 2° to 4° with the increasing period over the study region. A robust scheme of Debayle and Sambridge allowed us to conduct the tomographic inversion efficiently for massive datasets. Group velocity maps varying with period show lateral heterogeneities, well related to the geological and tectonic features in the study region. The 3D S-velocity model was constructed from the 1D structure inversion of the tomographic group velocity dispersion curves at each node. The obtained average crustal structure is similar to the PREM model, while the average mantle velocity is typically lower than the global average. The complicated 3D structures reveal three prominent features correlated with geological divisions: sea basin regions, island and arc regions, and continental regions. The derived crustal and lithospheric thicknesses range from ∼15 to >50 km and from ∼60 to >140 km, respectively, with the thinnest in the South China Sea, the thickest in eastern Tibet and the Yangtze Block, and the medium in the South China Fold Belt, Indochina, and island arc regions. Our results further confirm that (1) a Mesozoic subduction zone, which is interpreted as the tectonic weak zone during the Paleogene, exists along the South China margin; (2) the influence of the Indochina extrusion along the Red River Fault is limited for the South China Sea region; (3) there is a slab remnant of the proto-South China Sea beneath Borneo. New findings suggest that the Mesozoic subduction zone should be built into any evolution model for the region, as well as the other two major tectonic boundaries of the Red River Fault and proto-South China Sea subduction zone.     

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.     

S-wave velocity structure of the crust and upper mantle under southeastern China by surface wave dispersion analysis

Group velocity dispersion data of fundamental-mode Rayleigh and Love waves for 12 wave paths within southeastern China have been measured by applying the multiple-filter technique to the properly rotated three-component digital seismograms from two Seismic Research Observatory stations, TATO and CHTO. The generalized surface wave inversion technique was applied to these group velocity dispersion data to determine the S-wave velocity structures of the crust and upper mantle for various regions of southeastern China. The results clearly demonstrate that the crust and upper mantle under southeastern China are laterally heterogeneous. The southern China region south of 25°N and the eastern China region both have a crustal thickness of 30 km. The eastern Tibet plateau along the 100°E meridian has a crustal thickness of 60 km. Central China, consisting mainly of the Yangtze and Sino-Korean platforms, has a crustal thickness of 40 km. A distinct S-wave low-velocity layer at 10–20 km depth in the middle crust was found under wave paths in southeastern China. On the other hand, no such crustal low-velocity layer is evident under the eastern Tibet plateau. This low-velocity layer in the middle crust appears to reflect the presence of a sialic low-velocity layer perhaps consisting of intruded granitic laccoliths, or possibly the remnant of the source zone of widespread magmatic activities known to have taken place in these regions since the late Carboniferous.     

Crustal and upper mantle S-wave velocity structure beneath the Bransfield Strait (West Antarctica) from regional surface wave tomography

Regional surface wave tomography in the sub-Antarctic Scotia Sea is helpful in revealing the nature of the crust and the S-wave seismic velocity profile beneath the Bransfield Strait. The joint use of our regional network, global seismographic network stations and local temporary arrays provide better lateral resolution than that obtained in our previous studies concerning the Scotia Sea region.Tomographic analysis of data obtained using 10 broad band seismic stations and more than 300 regional events, shows that the Bransfield Basin is characterised by a strong group velocity reduction of 8% with respect to the surrounding areas, in the period range from 15 s to 50 s.The crustal and upper mantle models of the eastern, central and western Bransfield Basin are obtained by joint inversion of Rayleigh and Love local dispersion curves from 15 s to 50 s. In addition our data set is expanded to a broader period interval (1–80 s), in central Bransfield Strait in order to better constrain the upper mantle and shallow crust.The main results can be summarized as follows: (a) the crust thins distinctly from W toward E; the variation is consistent with the type of volcanism, earthquake distribution and bathymetric observations, (b) low upper mantle velocities (soft lid) extend down to depths exceeding 70 km as a consequence of elevated temperatures, (c) the crust beneath the central Bransfield Basin displays continental characteristics with a gradually increasing S-wave velocity distribution versus depth analogous to the East African Rift structure of Kenya, (d) negative velocity gradients are present in the lower crust beneath the eastern Bransfield Basin; these could be interpreted as magmatic bodies originating from decompression melting of the mantle.     

Fission track and fault kinematics analyses for new insight into the Late Cenozoic tectonic regime changes in West-Central Sulawesi (Indonesia)

A 3-D density model for the Cretan and Libyan Seas and Crete was developed by gravity modelling constrained by five 2-D seismic lines. Velocity values of these cross-sections were used to obtain the initial densities using the Nafe–Drake and Birch empirical functions for the sediments, the crust and the upper mantle. The crust outside the Cretan Arc is 18 to 24 km thick, including 10 to 14 km thick sediments. The crust below central Crete at its thickest section, has values between 32 and 34 km, consisting of continental crust of the Aegean microplate, which is thickened by the subducted oceanic plate below the Cretan Arc. The oceanic lithosphere is decoupled from the continental along a NW–SE striking front between eastern Crete and the Island of Kythera south of Peloponnese. It plunges steeply below the southern Aegean Sea and is probably associated with the present volcanic activity of the southern Aegean Sea in agreement with published seismological observations of intermediate seismicity. Low density and velocity upper mantle below the Cretan Sea with ρ  3.25 × 10³ kg/m³ and V_p velocity of compressional waves around 7.7 km/s, which are also in agreement with observed high heat flow density values, point out at the mobilization of the upper mantle material here. Outside the Hellenic Arc the upper mantle density and velocity are ρ ≥ 3.32 × 10³ kg/m³ and V_p = 8.0 km/s, respectively. The crust below the Cretan Sea is thin continental of 15 to 20 km thickness, including 3 to 4 km of sediments. Thick accumulations of sediments, located to the SSW and SSE of Crete, are separated by a block of continental crust extended for more than 100 km south of Central Crete. These deep sedimentary basins are located on the oceanic crust backstopped by the continental crust of the Aegean microplate. The stretched continental margin of Africa, north of Cyrenaica, and the abruptly terminated continental Aegean microplate south of Crete are separated by oceanic lithosphere of only 60 to 80 km width at their closest proximity. To the east and west, the areas are floored by oceanic lithosphere, which rapidly widens towards the Herodotus Abyssal plain and the deep Ionian Basin of the central Mediterranean Sea. Crustal shortening between the continental margins of the Aegean microplate and Cyrenaica of North Africa influence the deformation of the sediments of the Mediterranean Ridge that has been divided in an internal and external zone. The continental margin of Cyrenaica extends for more than 80 km to the north of the African coast in form of a huge ramp, while that of the Aegean microplate is abruptly truncated by very steep fractures towards the Mediterranean Ridge. Changes in the deformation style of the sediments express differences of the tectonic processes that control them. That is, subduction to the northeast and crustal subsidence to the south of Crete. Strike-slip movement between Crete and Libya is required by seismological observations.     

Lithospheric structure beneath the East China Sea revealed by Rayleigh-wave phase velocities

We explore the variations of Rayleigh-wave phase-velocity beneath the East China Sea in a broad period range (5–200 s). Rayleigh-wave dispersion curves are measured by the two-station technique for a total of 373 interstation paths using vertical-component broad-band waveforms at 32 seismic stations around the East China Sea from 6891 global earthquakes.The resulting maps of Rayleigh-wave phase velocity and azimuthal anisotropy provide a high resolution model of the lithospheric mantle beneath the East China Sea. The model exhibits four regions with different isotropic and anisotropic patterns: the Bohai Sea, belonging to the North China Craton, displays a continental signature with fast velocities at short periods; the Yellow Sea, very stable unit associated with low deformation, exhibits fast velocities and limited anisotropy; the southern part of the East China Sea, with high deformation and many fractures and faults, is related to slow velocities and high anisotropic signature; and the Ryukyu Trench shows high-velocity perturbations and slab parallel anisotropy.     

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.     

The Levantine Basin—crustal structure and origin

The origin of the Levantine Basin in the Southeastern Mediterranean Sea is related to the opening of the Neo-Tethys. The nature of its crust has been debated for decades. Therefore, we conducted a geophysical experiment in the Levantine Basin. We recorded two refraction seismic lines with 19 and 20 ocean bottom hydrophones, respectively, and developed velocity models. Additional seismic reflection data yield structural information about the upper layers in the first few kilometers. The crystalline basement in the Levantine Basin consists of two layers with a P-wave velocity of 6.0–6.4 km/s in the upper and 6.5–6.9 km/s in the lower crust. Towards the center of the basin, the Moho depth decreases from 27 to 22 km. Local variations of the velocity gradient can be attributed to previously postulated shear zones like the Pelusium Line, the Damietta–Latakia Line and the Baltim–Hecateus Line. Both layers of the crystalline crust are continuous and no indication for a transition from continental to oceanic crust is observed. These results are confirmed by gravity data. Comparison with other seismic refraction studies in prolongation of our profiles under Israel and Jordan and in the Mediterranean Sea near Greece and Sardinia reveal similarities between the crust in the Levantine Basin and thinned continental crust, which is found in that region. The presence of thinned continental crust under the Levantine Basin is therefore suggested. A β-factor of 2.3–3 is estimated. Based on these findings, we conclude that sea-floor spreading in the Eastern Mediterranean Sea only occurred north of the Eratosthenes Seamount, and the oceanic crust was later subducted at the Cyprus Arc.     

The lithosphere-asthenosphere system in Fennoscandia

We present a summary of the available information on Rayleigh-wave dispersion data for the Fennoscandian region. The observations have been combined to produce regional dispersion relations which have then been subjected to the “hedgehog” inversion procedure. The results are presented on a map outlining the thickness of the lid and the shear velocities in both the lid and the asthenosphere channel. Lid thickness up to around 135 km is found in the Bothnia-north-central Finland area with, if any, weak shear velocity contrast to the underlying layer. The surrounding areas are characterized by lid thickness up to around 75 km; a stronger low-velocity zone to lid contrast may be found in the Caledonian and Baltic Sea area (0.25÷0.45 km/s). Taking into account Moho depth data and the aforementioned results, a map of the lithosphere-asthenosphere system was derived.     

Crustal thinning beneath the Rwenzori region, Albertine rift, Uganda, from receiver-function analysis

The Rwenzori mountains in western Uganda, with a maximum elevation of more than 5,000 m, are located within the Albertine rift valley. We have deployed a temporary seismic network on the Ugandan side of the mountain range to study the seismic velocity structure of the crust and upper mantle beneath this section of the rift. We present results from a receiver-function study revealing a simple crustal structure along the eastern rift flank with a more or less uniform crustal thickness of about 30 km. The complexity of inner-crustal structures increases drastically within the Rwenzori block. We apply different inversion techniques to obtain reliable results for the thickness of the crust. The observations expose a significantly thinner crust beneath the Rwenzori range with thickness values ranging from about 20–28 km beneath northern and central parts of the mountains. Our study therefore indicates the absence of a crustal root beneath the Rwenzori block. Beneath the Lake Edward and Lake George basins we detect the top of a layer of significantly reduced S-wave velocity at 15 km depth. This low-velocity layer may be attributed to the presence of partial melt beneath a region of recent volcanic activity.     

Lower-crust S-wave velocity beneath western Yunnan Province from waveform inversion of dense seismic observations

We use teleseismic body waveforms to explore S-wave layered velocity structures beneath 30 portable digital seismic stations deployed around western Yunnan Province. Results show that the Moho depth in this region is ∼40 km and decreases in general from north to south, consistent with previous geophysical studies. Associated with this lateral variation of the Moho depth, the lower crust above the Moho discontinuity has a 15–25 km thick zone with an S-wave velocity lower than that of the upper crust. This lower velocity zone might be interpreted as a lower crust weak channel, which may mechanically partially decouple the upper-crust deformation from the underlying mantle. Thus, the inverted S-wave velocity structure could provide new evidence for the lateral flow of lower crust in the build-up of the south-eastern Tibetan plateau.     

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.     

Rayleigh wave dispersion in Papua New Guinea

Phase and group velocity dispersions of fundamental mode Rayleigh waves for the path Port Moresby-Rabaul have been computed. The structural interpretation of dispersion curves for the Solomon Sea reveals the channel shear velocity as about 4.35 km/s and the channel starts around 95 km below the earth's surface. In other words, the elastic velocities-depth structure beneath the Solomon Sea is very similar to that beneath the northern part of the Western Mediterranean Basin. Further, group velocity dispersion of fundamental mode Rayleigh waves at the PMG and the RAB seismic stations have been computed. The structural interpretation of a dispersion curve for the orogenic belt of Papua New Guinea reveals a crust of about 62.5 km thick and uppermost mantle shear velocity of about 4.3 km/s.     

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.     

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.     

Moho depth variation beneath southwestern Japan revealed from the velocity structure based on receiver function inversion

The Philippine Sea plate is subducting under the Eurasian plate beneath the Chugoku-Shikoku region, southwestern Japan. We have constructed depth contours for the continental and oceanic Mohos derived from the velocity structure based on receiver function inversion. Receiver functions were calculated using teleseismic waveforms recorded by the high-density seismograph network in southwestern Japan. In order to determine crustal velocity structure, we first improved the linearized time-domain receiver function inversion method. The continental Moho is relatively shallow ( 30 km) at the coastline of the Sea of Japan and at the Seto Inland Sea, and becomes deeper–greater than 40 km–around 35°N and 133.8°E. Near the Seto Inland Sea, a low-velocity layer of thickness 10 km lies under the continental Moho. This low-velocity layer corresponds to the subducting oceanic crust of the Philippine Sea plate. The oceanic Moho continues to descend from south to northwest and exhibits complicated ridge and valley features. The oceanic Moho runs around 25 km beneath the Pacific coast and 45 km beneath the Seto Inland Sea, and it extends to at least to 34.5°N. The depth variation of the Moho discontinuities is in good qualitative agreement with the concept of isostasy. From the configurations of both the continental and oceanic Mohos, we demonstrate that the continental lower crust and the subducting oceanic crust overlap beneath the southern and central part of Shikoku and that a mantle wedge may exist beneath the western and eastern part of Shikoku. The southern edge of the overlapping region coincides with the downdip limit of the slip area of a megathrust earthquake.     

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.