Upper mantle heterogeneities beneath Eastern Europe

P-wave travel-time residuals for seismograph stations in eastern Europe as reported by ISC for the years 1964–1977 were used for constructing a seismic image of upper mantle heterogeneities in the network region. For the depth range 0–100 km, dominant tectonic features like the Pannonian Basin and the Aegean Sea and western Turkey correlate well with pronounced velocity lows which a ppear to extenddown to a 300 km depth. The velocity anomaly patterns in the depth intervals 300–500 km and 500–600 km are broadly similar but quite different from those of shallower depths. The observed seismic heterogeneities are briefly discussed in terms of large-scale tectonic and geophysical (heat-flow) characteristics of eastern Europe.     

Seismotectonics of the Calabrian arc

Relocation of intermediate and deep earthquakes of Tyrrhenian Sea area through joint hypocenter determination for the period 1962–1979, has allowed a more detailed definition of the geometry of this peculiar Benioff zone. Earthquakes dip along a quasi-vertical plane to 250 km depth; there is a 50° dip in the 250–340 km depth range, and a low dip angle to 480 km depth. The structure sketched from the hypocenters is almost continuous, but most energy has been released in the 230–340 km depth interval. An evaluation of fault plane solutions of intermediate earthquakes in this area indicates predominance of down-dip compressions in the central part of the slab. At the border, strike-slip motion occurs independent of depth. Some earthquakes that occurred at intermediate depth (less than 100 km) along the Ionian margin of Calabria show predominance of reverse faulting, with the P-axis oriented SE-NW. However, shallow earthquakes in the Calabria-Sicily region indicate a more complex motion, with predominance of normal faulting. A possible interpretation of these features according to the available geological history, which involves subduction of continental lithosphere, is discussed.     

Structural characteristics off Miyagi forearc region, the Japan Trench seismogenic zone, deduced from a wide-angle reflection and refraction study

The Japan Trench subduction zone, located east of NE Japan, has regional variation in seismicity. Many large earthquakes occurred in the northern part of Japan Trench, but few in the southern part. Off Miyagi region is in the middle of the Japan Trench, where the large earthquakes (M > 7) with thrust mechanisms have occurred at an interval of about 40 years in two parts: inner trench slope and near land. A seismic experiment using 36 ocean bottom seismographs (OBS) and a 12,000 cu. in. airgun array was conducted to determine a detailed, 2D velocity structure in the forearc region off Miyagi. The depth to the Moho is 21 km, at 115 km from the trench axis, and becomes progressively deeper landward. The P-wave velocity of the mantle wedge is 7.9–8.1 km/s, which is typical velocity for uppermost mantle without large serpentinization. The dip angle of oceanic crust is increased from 5–6° near the trench axis to 23° 150 km landward from the trench axis. The P-wave velocity of the oceanic uppermost mantle is as small as 7.7 km/s. This low-velocity oceanic mantle seems to be caused by not a lateral anisotropy but some subduction process. By comparison with the seismicity off Miyagi, the subduction zone can be divided into four parts: 1) Seaward of the trench axis, the seismicity is low and normal fault-type earthquakes occur associated with the destruction of oceanic lithosphere. 2) Beneath the deformed zone landward of the trench axis, the plate boundary is characterized as a stable sliding fault plain. In case of earthquakes, this zone may be tsunamigenic. 3) Below forearc crust where P-wave velocity is almost 6 km/s and larger: this zone is the seismogenic zone below inner trench slope, which is a plate boundary between the forearc and oceanic crusts. 4) Below mantle wedge: the rupture zones of thrust large earthquakes near land (e.g. 1978 off Miyagi earthquake) are located beneath the mantle wedge. The depth of the rupture zones is 30–50 km below sea level. From the comparison, the rupture zones of large earthquakes off Miyagi are limited in two parts: plate boundary between the forearc and oceanic crusts and below mantle wedge. This limitation is a rare case for subduction zone. Although the seismogenic process beneath the mantle wedge is not fully clarified, our observation suggests the two possibilities: earthquake generation at the plate boundary overridden by the mantle wedge without serpentinization or that in the subducting slab.     

Stress and viscosity in the asthenosphere

Stresses and effective viscosities in the asthenosphere to a depth of 400 km are calculated on the basis of Weertmans “temperature method” i.e., on relating viscosity to the ratio of the temperature to the melting point (=homologous temperature). Some oceanic and continental geotherms and two melting point—depth curves, the dry pyrolite solidus and the forsterite₉₀ melting curve are used for the conversion of the homologous temperature to the effective viscosity. Two creep laws are considered, the linear, grain-size-dependent Nabarro—Herring (NH) creep law, and a power creep law, in which the creep rate is proportional to the third power of the stress. A plate tectonic model yields creep rates of 2 · 10⁻¹⁴ s⁻¹ for the oceanic and 3 · 10⁻¹⁵ s⁻¹ for the continental asthenosphere. These values are held constant for the calculations and may be valid for regions inside plates.The dry pyrolite mantle model results in high homologous temperatures in the asthenosphere below oceans (0.9), very low stresses (a few bars and lower) and shows a low viscosity “layer” of about 200-km thickness. Below continental shields the homologous temperature has a maximum value of 0.73, stresses are around 5–20 bar and the low-viscosity region is thicker and less pronounced than in the oceanic case. The Fo₉₀ mantle model generally gives lower homologous temperatures (maximum value below oceans beside active ridges 0.75). The stresses in the asthenosphere beneath oceans vary from a few bars to about 50 bar and below continents to about 100 bar. The low-viscosity region seems to reach great depths without forming a “channel”. The Figs. 1 and 2 show the approximate viscosity—depth distribution for the two mantle models under study.Assuming a completely dry mantle and a mean grain size of 5 mm, power law creep will be the dominating creep process in the asthenosphere. However, grains may grow in a high-temperature—low-stress regime (i.e., below younger oceans), an effect which will further diminish the influence of NH creep. In the upper 100–150 km of the earth some fluid phases may affect considerably creep processes.     

Seismic velocities of granulite-facies xenoliths from Central Ireland: Implications for lower crustal composition and anisotropy

V_p and V_s values have been measured experimentally and calculated for granulite-facies lower crustal xenoliths from central Ireland close to the Caledonian Iapetus suture zone. The xenoliths are predominantly foliated and lineated metapelitic (garnet–sillimanite–K-feldspar) granulites. Their metapelitic composition is unusual compared with the mostly mafic composition of lower crustal xenoliths world-wide. Based on thermobarometry, the metapelitic xenoliths were entrained from depths of c. 20–25 ± 3.5 km and rare mafic granulites from depths of 31–33 ± 3.4 km. The xenoliths were emplaced during Lower Carboniferous volcanism and are considered to represent samples of the present day lower crust.V_p values for the metapelitic granulites range between 6.26 and 7.99 km s^{− 1} with a mean value of 7.09 ± 0.4 km s^{− 1}. Psammite and granitic orthogneiss samples have calculated V_p values of 6.51 and 6.23 km s^{− 1}, respectively. V_s values for the metapelites are between 3.86 and 4.34 km s^{− 1}, with a mean value of 4.1 ± 0.15 km s^{− 1}. The psammite and orthogneiss have calculated V_s values of 3.95 and 3.97 km s^{− 1}, respectively.The measured seismic velocities correlate with density and with modal mineralogy, especially the high content of sillimanite and garnet. V_p anisotropy is between 0.15% and 13.97%, and a clear compositional control is evident, mainly in relation to sillimanite abundance. Overall V_s anisotropy ranges from 1% to 11%. Poisson's ratio (σ) lies between 0.25 and 0.35 for the metapelitic granulites, mainly reflecting a high V_p value due to abundant sillimanite in the sample with the highest σ. Anisotropy is probably a function of deformation associated with the closure of the Iapetus ocean in the Silurian as well as later extension in the Devonian. The orientation of the bulk strain ellipsoid in the lower crust is difficult to constrain, but lineation is likely to be NE–SW, given the strike-slip nature of the late Caledonian and subsequent Acadian deformation.When corrected for present-day lower crustal temperature, the experimentally determined V_p values correspond well with velocities from the ICSSP, COOLE I and VARNET seismic refraction lines. Near the xenolith localities, the COOLE I line displays two lower crustal layers with in situ V_p values of 6.85–6.9 and 6.9–8.0 km s^{− 1}, respectively. The upper (lower velocity) layer corresponds well with the metapelitic granulite xenoliths while the lower (higher velocity) layer matches that of the basic granulite xenoliths, though their metamorphic pressures suggest derivation from depths corresponding to the present-day upper mantle.     

Three-dimensional P-wave velocity image under the Carpathian arc

An inversion of P-wave travel time residuals from selected earthquakes in the distance range 30°–98° to two seismic station networks was used to model P-wave velocity anomalies down to 250 km depth. In the first inversion experiment a region between 43.5°–47.5°N and 21°–29°E was modelled, using 35 seismic stations, while in the second one a region between 44°–47°N and 25°–29°E was modelled, using 19 seismic stations. The 4-layer block model of the first inversion offers 19% reduction in residual variance, while the 5-layer block model of the second one offers 26% reduction, the rest being explained by noise and smaller scale heterogeneities. The obtained velocity anomalies correlate remarkably well with the gravity anomalies and with the tectonic model for the Vrancea region of Fuchs et al. (1979).     

Temperature-depth profiles in Czechoslovakia and some adjacent areas derived from heat-flow measurements, deep seismic sounding and other geophysical data

The results of seismic measurements along three deep seismic sounding (DSS) profiles on the territory of Czechoslovakia and in adjacent countries have provided sufficient material about the crustal structure and the depth of the Moho discontinuity. These data, together with gravity and aeromagnetic data and the determinations of heat-flow values, were used to select several locations where the temperature—depth profiles were calculated. The Moho temperature of about 500 C beneath the Bohemian Massif increases to 800–1000 C and even more beneath the inner Neogene depressions of the Carpathian system. The regional differences in mantle heat-flow contribution between both these provinces may reach 1 μcal. cm⁻² sec⁻¹; such a variation in energy inflow may then be the driving force for the geological evolution. The geophysical implications of different thermal structure of the crust are discussed. Because of high subsurface temperatures in the Hungarian basin, partial melting at a depth of about 30 km may not be excluded.     

Results of Scripps long range profiles in the Pacific

During 1976 the first installment of a long range seismic profile was conducted in the North Pacific to a range of 600 km using shots to two tons in size. The line was shot to a closely-spaced array of Scripps ocean bottom seismographs and was parallel to magnetic anomaly 32 at an age of approximately 70 · 10⁶ yr. The line extended between the Clarion and Molokai Fracture Zones and did not cross any major topographic features. Linearized and extremal travel-time inversions were conducted to provide bounds on the compressional velocity as a function of depth. The velocity does not exceed 8.4 km s⁻¹ to a depth of 60 km at which point the data no longer provide any resolution. The constraints on the acceptable models were improved by using array processing methods to measure phase velocity and synthetic seismogram techniques to model phase and amplitude information. The oceanic crust is composed of a series of gradients with no first order discontinuities. The “Moho” is smeared out over a depth of 1.5–2.0 km even though “wide-angle reflections” from the Moho, the phase P_MP, are clearly seen in the data. The upper lithosphere is characterized by a general tendency for the velocity to decrease with depth and the tendency is occasionally overwhelmed (at about 27 and 52 km depth) by rapid velocity changes perhaps associated with phase or compositional changes.     

A tentative 2D thermal model of central India across the Narmada-Son Lineament (NSL)

This work deals with 2D thermal modeling in order to delineate the crustal thermal structure of central India along two Deep Seismic Sounding (DSS) profiles, namely Khajuriakalan–Pulgaon and Ujjan–Mahan, traversing the Narmada-Son-Lineament (NSL) in an almost north–south direction. Knowledge of the crustal structure and P-wave velocity distribution up to the Moho, obtained from DSS studies, has been used for the development of the thermal model. Numerical results reveal that the Moho temperature in this region of central India varies between 500 and 580 °C. The estimated heat flow density value is found to vary between 46 and 49 mW/m². The Curie depth varies between 40 and 42 km and is in close agreement with the Curie depth (40±4 km) estimated from the analysis of MAGSAT data. Based on the present work and previous work, it is suggested that the major part of peninsular India consisting of the Wardha–Pranhita Godavari graben/basin, Bastar craton and the adjoining region of the Narmada Son Lineament between profiles I and III towards the north and northwest of the Bastar craton are characterized with a similar mantle heat flow density value equal to 23 mW/m². Variation in surface heat flow density values in these regions are caused by variation in the radioactive heat production and fluid circulation in the upper crustal layer.     

On the formation of continental silicic melts in thermochemical mantle convection models: implications for early Earth

Seafloor magnetotelluric (MT) data were collected at seven sites across the Hawaiian hot spot swell, spread approximately evenly between 120 and 800 km southwest of the Hawaiian-Emperor island chain. All data are consistent with an electrical strike direction of 300°, aligned along the seamount chain, and are well fit using two-dimensional (2D) inversion. The major features of the 2D electrical model are a resistive lithosphere underlain by a conductive lower mantle, and a narrow, conductive, ‘plume’ connecting the surface of the islands to the lower mantle. This plume is required; without it the swell bathymetry produces a large divergence of the along-strike and across-strike components of the MT fields, which is not seen in the data. The plume radius appears to be less than 100 km, and its resistivity of around 10 Ωm, extending to a depth of 150 km, is consistent with a bulk melt fraction of 5–10%.A seismic low velocity region (LVR) observed by Laske et al. [Laske, G., Phipp Morgan, J., Orcutt, J.A., 1999. First results from the Hawaiian SWELL experiment, Geophys. Res. Lett. 26, 3397–3400] at depths centered around 60 km and extending 300 km from the islands is not reflected in our inverse model, which extends high lithospheric resistivities to the edge of the conductive plume. Forward modeling shows that resistivities in the seismic LVR can be lowered at most to 30 Ωm, suggesting a maximum of 1% connected melt and probably less. However, a model of hot subsolidus lithosphere of 10² Ωm (1450–1500 °C) within the seismic LVR increasing to an off-swell resistivity of >10³ Ωm (<1300 °C) fits the MT data adequately and is also consistent with the 5% drop in seismic velocities within the LVR. This suggests a ‘hot, dry lithosphere’ model of thermal rejuvination, or possibly underplated lithosphere depleted in volatiles due to melt extraction, either of which is derived from a relatively narrow mantle plume source of about 100 km radius. A simple thermal buoyancy calculation shows that the temperature structure implied by the electrical and seismic measurements is in quantitative agreement with the swell bathymetry.     

Configuration of subducting Philippine Sea plate and crustal structure in the central Japan region

A seismic experiment with six explosive sources and 391 seismic stations was conducted in August 2001 in the central Japan region. The crustal velocity structure for the central part of Japan and configuration of the subducting Philippine Sea plate were revealed. A large lateral variation of the thickness of the sedimentary layer was observed, and the P-wave velocity values below the sedimentary layer obtained were 5.3–5.8 km/s. P-wave velocity values for the lower part of upper crust and lower crust were estimated to be 6.0–6.4 and 6.6–6.8 km/s, respectively. The reflected wave from the upper boundary of the subducting Philippine Sea plate was observed on the record sections of several shots. The configuration of the subducting Philippine Sea slab was revealed for depths of 20–35 km. The dip angle of the Philippine Sea plate was estimated to be 26° for a depth range of about 20–26 km. Below this depth, the upper boundary of the subducting Philippine Sea plate is distorted over a depth range of 26–33 km. A large variation of the reflected-wave amplitude with depth along the subducting plate was observed. At a depth of about 20–26 km, the amplitude of the reflected wave is not large, and is explained by the reflected wave at the upper boundary of the subducting oceanic crust. However, the reflected wave from reflection points deeper than 26 km showed a large amplitude that cannot be explained by several reliable velocity models. Some unique seismic structures have to be considered to explain the observed data. Such unique structures will provide important information to know the mechanism of inter-plate earthquakes.     

Upper mantle structure of the Northern Eurasia from peaceful nuclear explosion data

Several long-range seismic profiles were carried out in Russia with Peaceful Nuclear Explosions (PNE). The data from 25 PNEs recorded along these profiles were used to compile a 3-D upper mantle velocity model for the central part of the Northern Eurasia. 2-D crust and upper mantle models were also constructed for all profiles using a common methodology for wavefield interpretation. Five basic boundaries were traced over the study area: N1 boundary (velocity level, V = 8.35 km/s; depth interval, D = 60–130 km), N2 (V = 8.4 km/s; D = 100–140 km), L (V = 8.5 km/s; D = 180–240 km) and H (V = 8.6 km/s; D = 300–330 km) and structural maps were compiled for each boundary. Together these boundaries describe a 3-D upper mantle model for northern Eurasia. A map characterised the velocity distribution in the uppermost mantle down to a depth of 60 km is also presented. Mostly horizontal inhomogeneity is observed in the uppermost mantle, and the velocities range from the average 8.0–8.1 km/s to 8.3–8.4 km/s in some blocks of the Siberian Craton. At a depth of 100–200 km, the local high velocity blocks disappear and only three large anomalies are observed: lower velocities in West Siberia and higher velocities in the East-European platform and in the central part of the Siberian Craton. In contrast, the depths to the H boundary are greater beneath the craton and lower beneath in the West Siberian Platform. A correlation between tectonics, geophysical fields and crustal structure is observed. In general, the old and cold cratons have higher velocities in the mantle than the young platforms with higher heat flows.Structural peculiarities of the upper mantle are difficult to describe in form of classical lithosphere–asthenosphere system. The asthenosphere cannot be traced from the seismic data; in contrary the lithosphere is suggested to be rheologically stratified. All the lithospheric boundaries are not simple discontinuities, they are heterogeneous (thin layering) zones which generate multiphase reflections. Many of them may be a result of fluids concentrated at some critical P–T conditions which produce rheologically weak zones. The most visible rheological variations are observed at depths of around 100 and 250 km.     

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.     

Anatolian surface wave evaluated at GEOFON Station ISP Isparta, Turkey

We have studied seismic surface waves of 255 shallow regional earthquakes recently recorded at GEOFON station ISP (Isparta, Turkey) and have selected these 52 recordings with high signal-to-noise ratio for further analysis. An attempt was made by the simultaneous use of the Rayleigh and Love surface wave data to interpret the planar crust and uppermost mantle velocity structure beneath the Anatolian plate using a differential least-square inversion technique. The shear-wave velocities near the surface show a gradational change from approximately 2.2 to 3.6 km s^− 1 in the depth range 0–10 km. The mid-crustal depth range indicating a weakly developed low velocity zone has shear-wave velocities around 3.55 km s^− 1. The Moho discontinuity characterizing the crust–mantle velocity transition appears somewhat gradual between the depth range  25–45 km. The surface waves approaching from the northern Anatolia are estimated to travel a crustal thickness of  33 km whilst those from the southwestern Anatolia and part of east Mediterranean Sea indicate a thicker crust at  37 km. The eastern Anatolia events traveled even thicker crust at  41 km. A low sub-Moho velocity is estimated at  4.27 km s^− 1, although consistent with other similar studies in the region. The current velocities are considerably slower than indicated by the Preliminary Reference Earth Model (PREM) in almost all depth ranges.     

Lithosphere structure of European sedimentary basins from regional three-dimensional gravity modelling

A three-dimensional (3D) density model, approximated by two regional layers—the sedimentary cover and the crystalline crust (offshore, a sea-water layer was added), has been constructed in 1° averaging for the whole European continent. The crustal model is based on simplified velocity model represented by structure maps for main seismic horizons—the “seismic” basement and the Moho boundary. Laterally varying average density is assumed inside the model layers. Residual gravity anomalies, obtained by subtraction of the crustal gravity effect from the observed field, characterize the density heterogeneities in the upper mantle. Mantle anomalies are shown to correlate with the upper mantle velocity inhomogeneities revealed from seismic tomography data and geothermal data. Considering the type of mantle anomaly, specific features of the evolution and type of isostatic compensation, the sedimentary basins in Europe may be related into some groups: deep sedimentary basins located in the East European Platform and its northern and eastern margins (Peri-Caspian, Dnieper–Donets, Barents Sea Basins, Fore–Ural Trough) with no significant mantle anomalies; basins located on the activated thin crust of Variscan Western Europe and Mediterranean area with negative mantle anomalies of −150 to −200×10⁻⁵ ms⁻² amplitude and the basins associated with suture zones at the western and southern margins of the East European Platform (Polish Trough, South Caspian Basin) characterized by positive mantle anomalies of 50–150×10⁻⁵ ms⁻² magnitude. An analysis of the main features of the lithosphere structure of the basins in Europe and type of the compensation has been carried out.     

Heat flow, heat production and fission track data from the Hercynian basement around the Provençal Basin (Western Mediterranean)

In the complex structural framework of the Western Mediterranean. Hercynian areas are expected to be thermally preserved from the recent tectonic evolution. The thermal regime of these areas is studied using heat flow, heat production and fission track data. The surface heat flow is significantly higher in Corsica (76 ± 10 mW m⁻²) than in the Maures and Estérel (58 ± 2 mW m⁻²). Neither heat production nor erosion subsequent to the Alpine orogeny in Corsica can explain such a difference. It is suggested that a deep thermal source related to the asymmetric evolution of the Provençal basin could explain the higher heat flow in Corsica. A model of thermal structure based on the present day thermal regime of the Maures and Estérei is proposed for the stable Hercynian crust in this area. The mantle heat flow is 20–25 mW m⁻² and the temperature at Moho level is 375–500°C, depending on the thermal parameter distribution with depth.     

Crustal structure and transform margin development south of Svalbard based on ocean bottom seismometer data

The Barents Sea is located in the northwestern corner of the Eurasian continent, where the crustal terrain was assembled in the Caledonian orogeny during Late Ordovician and Silurian times. The western Barents Sea margin developed primarily as a transform margin during the early Tertiary. In the northwestern part south of Svalbard, multichannel reflection seismic lines have poor resolution below the Permian sequence, and the early post-orogenic development is not well known here. In 1998, an ocean bottom seismometer (OBS) survey was collected southwest to southeast of the Svalbard archipelago. One profile was shot across the continental transform margin south of Svalbard, which is presented here. P-wave modeling of the OBS profile indicates a Caledonian suture in the continental basement south of Svalbard, also proposed previously based on a deep seismic reflection line coincident with the OBS profile. The suture zone is associated with a small crustal root and westward dipping mantle reflectivity, and it marks a boundary between two different crystalline basement terrains. The western terrain has low (6.2–6.45 km s⁻¹) P-wave velocities, while the eastern has higher (6.3–6.9 km s⁻¹) velocities. Gravity modeling agrees with this, as an increased density is needed in the eastern block. The S-wave data predict a quartz-rich lithology compatible with felsic gneiss to granite within and west of the suture zone, and an intermediate lithological composition to the east. A geological model assuming westward dipping Caledonian subduction and collision can explain the missing lower crust in the western block by subduction erosion of the lower crust, as well as the observed structuring. Due to the transform margin setting, the tectonic thinning of the continental block during opening of the Norwegian-Greenland Sea is restricted to the outer 35 km of the continental block, and the continent–ocean boundary (COB) can be located to within 5 km in our data. Distinct from the outer high commonly observed on transform margins, the upper part of the continental crust at the margin is dominated by two large, rotated down-faulted blocks with throws of 2–3 km on each fault, apparently formed during the transform margin development. Analysis of the gravity field shows that these faults probably merge to one single fault to the south of our profile, and that the downfaulting dominates the whole margin segment from Spitsbergen to Bjørnøya. South of Bjørnøya, the faulting leaves the continental margin to terminate as a graben 75 km south of the island. Adjacent to the continental margin, there is no clear oceanic layer 2 seismic signature. However, the top basement velocity of 6.55 km s⁻¹ is significantly lower than the high (7 km s⁻¹) velocity reported earlier from expanding spread profiles (ESPs), and we interpret the velocity structure of the oceanic crust to be a result of a development induced by the 7–8-km-thick sedimentary overburden.     

Lower crustal intrusions beneath the southern Baikal Rift Zone: Evidence from full-waveform modelling of wide-angle seismic data

The Cenozoic Baikal Rift Zone (BRZ) is situated in south-central Siberia in the suture between the Precambrian Siberian Platform and the Amurian plate. This more than 2000-km long rift zone is composed of several individual basement depressions and half-grabens with the deep Lake Baikal at its centre. The BEST (Baikal Explosion Seismic Transect) project acquired a 360-km long, deep seismic, refraction/wide-angle reflection profile in 2002 across southern Lake Baikal. The data from this project is used for identification of large-scale crustal structures and modelling of the seismic velocities of the crust and uppermost mantle. Previous interpretation and velocity modelling of P-wave arrivals in the BEST data has revealed a multi layered crust with smooth variation in Moho depth between the Siberian Platform (41 km) and the Sayan-Baikal fold belt (46 km). The lower crust exhibits normal seismic velocities around the rift structure, except for beneath the rift axis where a distinct 50–80-km wide high-velocity anomaly (7.4–7.6 ± 0.2 km/s) is observed. Reverberant or “ringing” reflections with strong amplitude and low frequency originate from this zone, whereas the lower crust is non-reflective outside the rift zone. Synthetic full-waveform reflectivity modelling of the high-velocity anomaly suggests the presence of a layered sequence with a typical layer thickness of 300–500 m coinciding with the velocity anomaly. The P-wave velocity of the individual layers is modelled to range between 7.4 km/s and 7.9 km/s. We interpret this feature as resulting from mafic to ultra-mafic intrusions in the form of sills. Petrological interpretation of the velocity values suggests that the intrusions are sorted by fractional crystallization into plagioclase-rich low-velocity layers and pyroxene- and olivine-rich high-velocity layers. The mafic intrusions were probably intruded into the ductile lower crust during the main rift phase in the Late Pliocene. As such, the intrusive material has thickened the lower crust during rifting, which may explain the lack of Moho uplift across southern BRZ.     

Laboratory measurement of P-wave velocity in crustal and upper mantle xenoliths from Ichino-megata, NE Japan: ultrabasic hydrous lower crust beneath the NE Honshu arc

Ultrasonic laboratory measurements of P-wave velocity (Vp) were carried out up to 1.0 GPa in a temperature range of 25–400 °C for crustal and mantle xenoliths of Ichino-megata, northeast Japan. The rocks used in the present study cover a nearly entire range of lithological variation of the Ichino-megata xenoliths and are considered as representative rock samples of the lower crust and upper mantle of the back arc side of the northeast (NE) Honshu arc. The Vp values measured at 25 °C and 1.0 GPa are 6.7–7.2 km/s for the hornblende gabbros (38.6–46.9 wt.% SiO₂), 7.2 km/s for the hornblende-pyroxene gabbro (43.8 wt.% SiO₂), 6.9–7.3 km/s for the amphibolites (36.1–44.3 wt.% SiO₂), 8.0–8.1 km/s for the spinel lherzolites (46.2–47.2 wt.% SiO₂) and 6.30 km/s for the biotite granite (72.1 wt.% SiO₂). Combining the present data with the Vp profile of the NE Honshu arc [Iwasaki, T., Kato, W., Moriya, T., Hasemi, A., Umino, N., Okada, T., Miyashita, K., Mizogami, T., Takeda, T., Sekine, S., Matsushima, T., Tashiro, K., Miyamachi, H. 2001. Extensional structure in northern Honshu Arc as inferred from seismic refraction/wide-angle reflection profiling. Geophys. Res. Lett. 28 (12), 2329–2332], we infer that the 15 km thick lower crust of the NE Honshu arc is composed of amphibolite and/or hornblende (±pyroxene) gabbro with ultrabasic composition. The present study suggests that the Vp range of the lower crustal layer (6.6–7.0 km/s) in the NE Honshu arc, which is significantly lower than that obtained from various seismic measurements (e.g. the northern Izu-Bonin-Mariana arc: 7.1–7.3 km/s), is due to the thick hydrous lower crustal layer where hornblende, plagioclase and magnetite are dominant.     

Ray path interpretation of the crustal structure beneath Saudi Arabia

Interpretation of a long-range seismic refraction line in Saudi Arabia has shown that beneath the Arabian Shield velocity generally increases with depth, from about 6 km s⁻¹ at the surface to about 7 km s⁻¹ at the top of the crust-mantle transition zone. The base of this transition zone (Moho) occurs at 37–44 km in depth. Intracrustal discontinuities can also be recognized, the most important being in the 10–20 km-depth range and separating the upper from the lower crust. Laterally, the variations in the intracrustal discontinuities and the total crustal thickness can be correlated with previously defined tectonic regions. Beneath the Red Sea shelf and coastal plain the crust, including 4 km of sediments, is only 15–17.5 km thick. With the aid of both seismic and gravity data an abrupt, steeply dipping transition from the crust of the Red Sea shelf and coastal plain to that of the Arabian Shield has been derived. With a jump of more than 20 km in Moho depth, this appears to be the major discontinuity between the Red Sea depression and the Arabian continental shield.