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).     

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.     

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.     

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.     

Focal mechanism solution of the Himalayan earthquake of February 14, 1977

Through a closely spaced local network of seismic stations in Himachal Pradesh, India, supplemented by worldwide P-wave first-motion data, the source mechanism of the February 14, 1977 earthquake which occurred very close to the Rawalpindi area in Pakistan has been determined. The fault-plane solution as reported earlier for this event by Seeber and Armbruster (1979) showed thrust faulting. The reliability of their solution has been tested using more P-wave first-motion data from near Indian stations within the epicentral distance of 4°–7°, as well as from distant stations. The inclusion of data from these stations completely changed the type of faulting from thrust to normal type. The new solution parameters have been briefly discussed in relation to the local geological faults/thrusts.     

On the P-wave velocity and plate-tectonics implications for the Tyrrhenian deep-earthquake zone

P-wave velocities in the Tyrrhenian mantle have been determined for the 230–480 km depth range. Analysis of P-wave travel times for a set of Tyrrhenian deep earthquakes gives a velocity-distribution law which shows different behaviours in the 230–300 km and 300–480 km depth intervals. For the first interval the velocity gradient is 0.64 · 10⁻² sec⁻¹ and for the second one it is 0.59 · 10⁻² sec⁻¹. At a depth of 300 km the velocity decreases rapidly from 8.75 to 8.43 km/sec.The results have been analyzed in the framework of a Tyrrhenian structural model characterized by a lithospheric slab dipping 55–60° in the WNW direction.It is also pointed out that the analysis of some geodynamic features of the slabs of Pacific island arcs carried out by Oliver et al. (1973) and Sleep (1973) can be applied to the Tyrrhenian mantle geodynamic features.     

The earthquake of March 4, 1977, its recording peculiarities and source parameters

A detailed analysis of recording peculiarities at seismic stations of the Uniform System of Seismic Observations (USSO) is presented a complicated nature of the source being shown. Consideration is given to parameters of the earthquake source, including the seismic moment and the length of the rupture.Comparison of magnitudes M_LH and M_PV indicates an anomalous attenuation in surface waves, itis is 3–4 times weaker than it had been noticed in case of other intermediate-depth Carpathian earthquakes.On the basis of comparison of the logarithm of the ratio of P-wave spectra at different epicentral distances (30° –70° ), the fac tor characterizing the absorption of P wave is found to remain practically unchanged.Average value of the seismic moment is estimated to be 2.6 × 10²⁷ dyne × cm, the most reasonable length of the rupture 58 km, and its focus 100 –130 km. The source parameters of the earthquake in question are compared with those of the earthquake of November 10, 1940.     

Teleseismic p-wave traveltime residuals and deep structure of the Aegean region

Teleseismic P-wave traveltime residuals have been measured at the Greek seismic stations with respect to the Herrin 68 tables. In spite of the large scatter, some insight into crustal and upper-mantle structure of the Aegean region can be gained. The average absolute residuals (observed minus Herrin traveltimes) are of the order of + 2 s. The most plausible interpretation is an efficient low-velocity zone in the upper mantle. Simple estimates of densities and subsequently gravity with the aid of Birch's law suggest that the Aegean region is underlain by hot expanded upper mantle, perhaps involving partial melting. The relative P residuals (observed minus Herrin traveltime differences between a station and Athens) are generally positive and can be interpreted with lateral variations of the LVZ or of the crust. The latter interpretation is supported in some cases by seismic refraction data. The azimuthal variation of the relative residuals at stations on the non-volcanic arc bears a distinct relation with the arc orientation. At Archangelos (Rhodes) where we “see” through the Benioff zone, the residuals from N to W are between -1 and -2 s and indicate a high-velocity slab sinking below the Aegean sea. At Vamos (Crete), Valsamata (Kephallenia), and Joanina (Pindus Range) the largest (smallest) residuals are along directions parallel (perpendicular) to the arc. This can be interpreted by crustal thickening under the sedimentary arc and/or by velocity anisotropy with the maximum perpendicular to the arc. On the whole, our study supports the hypothesis that the Aegean region is a trench—island-arc—marginal-sea system.     

Seismicity and structural investigations of the Romanian Vrancea Region: Evidence for azimuthal variations of p-wave velocity and poisson's ratio

A study of the shallow and intermediate depth seismicity of the Romanian Vrancea region in the period 1964–1981 has been performed. The seismic events have been relocated by a standard location procedure using a regional velocity model. From the temporal and spatial distribution of the seismic activity, aspects of the seismicity before the large March 4, 1977 earthquake are treated, in particular the seismic gap in space and time prior to this event, found by Mârze (1979), which is critically discussed and revised. The concept of the precursor time/magnitude relationships of different authors is applied and its validity to the Vrancea region assessed. The hypocentral distribution shows that the intermediate depth seismic activity is confined to a small volume with dimensions of only some tens of kilometers. The results are interpreted in terms of the tectonics of the region. From an analysis of the travel-time residuals at different local stations, evidence for lateral velocity heterogeneities beneath the region is obtained e.g. a high velocity zone southeastwards of the Carpathian chain. Finally mean ratios, (i.e. Poisson's ratios), for various stations are calculated from P- and S-wave travel times. They show azimuthal variations of up to 6% for stations within the area where the intermediate seismic activity occurs in comparison with the station Focsani, situated eastwards in the Carpathian foredeeps. All these results are compatible with the plate tectonic concept for the Vrancea region, that is the subduction of an oceanic lithospheric slab under the Carpathian mountain arc, giving rise to such a highly active seismic zone.     

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.     

Effects of subducting lithosphere on p-wave amplitudes and residuals

In studying the structure of the lithosphere and asthenosphere using the kinematic and dynamic parameters of seismic waves, in any area, the tectonics of focal regions should also be considered. Subduction zones represent large-scale inhomogeneities which affect the propagation of seismic waves both at small and teleseismic epicentral distances.A study on the magnitude corrections of seismic stations in Central Europe and Scandinavia revealed that even in the case of close seismic stations the observed differences of the amplitudes of teleseismic P waves depend on the strike and dip of the sinking plates in the Northwestern Pacific. The smaller the angle between a seismic ray and the bottom of the lithospheric plate, the larger the numerical decrease of the magnitude corrections. A similar dependence was found for the P-wave residuals; the rays propagating more along the slab are accelerated more. Besides that, we have observed systematic differences in the P-wave residuals of about 2 s for foci with different positions within the subducting lithosphere of the Kurile arc in investigating the dependence of the residuals on the epicentral distance in the direction of subduction. The differences in the residuals disappear in this region at an epicentral distance of about 52°.     

Crustal anisotropy in southwest Ireland from analysis of controlled source shear-wave data

Results from a travel-time analysis of three-component shear-wave (S-wave) data recorded in southwest Ireland during a controlled source seismic experiment have been used to investigate the magnitude of crustal anisotropy. The data used were recorded from 20 in-line shots on three-component short-period stations deployed at approximately 1-km spacing along two parallel profiles. Analysis of the travel-time differences between vertically and horizontally polarised S-waves recorded on vertical, radial and transverse seismometer components was undertaken using seismic phases travelling near the Earth's surface (S_g) and reflected from the Moho (S_mS). Travel-time differences between the components for both phases scatter largely within the range ± 0.2 s, which is about the uncertainty in the measurements, with no observed coherent variation with shot-receiver offset. Synthetic S-wave seismograms were also computed from 1-D S-wave velocity models with varying degrees of anisotropy in the upper and in the lower crusts. Travel-time differences of S_g and S_mS phases picked from these synthetic seismograms confirm that for anisotropies with probable symmetries of magnitude 1–2% in either the upper or lower crust should result in an observable variation of the travel-time differences between the transverse and radial, and transverse and vertical components with source–receiver offset. The study shows that crustal anisotropy does not contribute significantly to the marked anisotropy recently deduced from SKS and SKKS measurements in Ireland, which is therefore confirmed to reside at sub-crustal and deeper mantle levels.     

Teleseismic P-wave residual investigations at Shillong,India

The northeast India region is seismically very active and it has experienced two large earthquakes of magnitude 8.7 during the last eight decades (1897 and 1950). We have analysed teleseismic P-wave residuals at Shillong, the only reliable seismic station operating in the region, to investigate a possible association of travel-time residual anomaly with earthquake occurrence. The period covered is from October 1964 through March 1976. The total number of events is 9479, including 1767 events with depth >/ 100 km. Six-monthly average residuals have been calculated. The standard deviations are less than 0.10 sec for these data sets. During the period of investigations, no major earthquake took place close to Shillong. The earthquake of June 1, 1969 with a magnitude (M_b) of 5.0, at an epicentral distance of 20 km from Shillong is the only significant event. This earthquake is found to be associated with a travel-time increase with a maximum amplitude of 0.4 sec. It appears that, in general, the P-wave velocity has decreased in the neighbourhood of Shillong since 1969. A quadrant-wise analysis of residuals indicates that the residual anomaly is most prominent in the SE quadrant from Shillong.     

Lithospheric structure of the Tornquist Zone resolved by nonlinear P and S teleseismic tomography along the TOR array

The main aim of the TOR project is to study the lithospheric–asthenospheric boundary structure under the Sorgenfrei–Tornquist Zone, across northern Germany, Denmark and southern Sweden. Relative arrival-time residuals of teleseismic P and S phases from 51 earthquakes, recorded by 150 seismic stations along the TOR array, were used to delineate the transition zone in the studied area. The effects of crustal structures were investigated by correcting the teleseismic residuals for travel-time variations in the crust based on a 3D crustal model derived from other data. The inversion was carried out for S phases. The results were then compared with the corresponding P-wave models. As expected, the derived models show that the relatively old and cold Baltic Shield has higher velocity at depth than the younger lithosphere farther South. The models show two sharp and distinct increases in depth to velocities which are low compared to our reference model, as we move from South to North. The location and sharpness of these boundaries suggests that the features resolved are, at least partially, compositional in origin, presumably related to mantle depletion. A sharp and steep subcrustal boundary is found roughly coincident with the southern edge of Sweden. This is below where the edge of the Baltic Shield is usually placed, based on surface geological evidence (the Sorgenfrei–Tornquist Zone). Another less significant transition is recognised more or less beneath the Elbe-lineament. Relatively high d(V_p / V_s) ratios under the central part of the profile (Denmark) indicate relatively low S-velocity in an area where a gravity high supports the hypothesis of extensive mafic intrusions.     

P-wave travel-time anomalies: Aleutian-Alaska region

The P-wave travel-time residuals in the Aleutian-Alaska region have been obtained by combining data from three underground nuclear tests so far carried out by the Atomic Energy Commission of the United States in Amchitka Island. The travel-time residuals show close correlation with the tectonism of the area. Attempts have been made to interpret qualitatively the pattern of the residuals in regard to the plate tectonics of the area. Areas of negative residuals appear interpretable as due to the underthrusted high-velocity lithospheric plate while areas with positive residuals seem to be associated with normal crust and upper mantle.     

Investigation of the lithosphere beneath the Vogelsberg volcanic complex with P-wave travel time residuals

With the aim of investigating the P-wave velocity structure below the Tertiary volcano Vogelsberg, a network of 10 mobile short period seismograph stations was installed in May 1987 for a period of 20 months. P-Wave travel time residuals relative to the station Kleiner Feldberg/Taunus (TNS) were determined for 168 seismic events using the Jeffreys - Bullen travel time tables. At all stations the relative residuals showed a positive sign, indicating a low velocity zone beneath the Vogelsberg. Maxima were found in the northern part of the Vogelsberg (station VAD +0.5 s) and in the region of the Amöneburger Basin (station RAU +0.28 s).The travel time residuals were inverted using the tomographic inversion method of Aki et al. (1977). The slowness perturbations of the single blocks were calculated relative to a crustal and upper mantle model of the Rhenish Massif. The results show an intracrustal low velocity body (about –9%) striking in a Variscan direction and underlying the north-eastern part of the Vogelsberg, and another velocity minimum (about – 6%) in the region of the Am6neburger Basin. In the lower crust and the upper mantle the velocities are reduced by about 4% relative to the starting model.The Variscan alignment of the low velocity zone under the Vogelsberg correlates with results of other geological studies. It can be assumed that during the rifting phase of the Upper Rhinegraben Variscan lineations have been reactivated, favouring uprising of magma along these old structures. The position and extension of the low velocity zone correlate with the assumed sediment distributions in the area of investigation. This may account for about one-half of the observed anomaly. The reason for the velocity reduction of about 4% in the entire underground region of the Vogelsberg down to a depth of about 70 km can be explained by the intensive fracturing of the lithosphere, caused by thermal and pressure gradients during the magma eruption process.     

Crustal structure of the Eastern Alps along the TRANSALP profile from wide-angle seismic tomography

The objective of the TRANSALP project is an investigation of the Eastern Alps with regard to their deep structure and dynamic evolution. The core of the project is a 340-km-long seismic profile at 12°E between Munich and Venice. This paper deals with the P-wave velocity distribution as derived from active source travel time tomography. Our database consists of Vibroseis and explosion seismic travel times recorded at up to 100 seismological stations distributed in a 30-km-wide corridor along the profile. In order to derive a velocity and reflector model, we simultaneously inverted refractions and reflections using a derivative of a damped least squares approach for local earthquake tomography. 8000 travel time picks from dense Vibroseis recordings provide the basis for high resolution in the upper crust. Explosion seismic wide-angle reflection travel times constrain both deeper crustal velocities and structure of the crust–mantle boundary with low resolution. In the resulting model, the Adriatic crust shows significantly higher P-wave velocities than the European crust. The European Moho is dipping south at an angle of 7°. The Adriatic Moho dips north with a gentle inclination at shallower depths. This geometry suggests S-directed subduction. Azimuthal variations of the first-break velocities as well as observations of shear wave splitting reveal strong anisotropy in the Tauern Window. We explain this finding by foliations and laminations generated by lateral extrusion. Based on the P-wave model we also localized almost 100 local earthquakes recorded during the 2-month acquisition campaign in 1999. Seismicity patterns in the North seem related to the Inn valley shear zone, and to thrusting of Austroalpine units over European basement. The alignment of deep seismicity in the Trento-Vicenza region with the top of the Adriatic lower crust corroborates the suggestion of a deep thrust fault in the Southern Alps.     

A dynamic model of the Hellenic Arc deduced from geophysical data

Combined gravity and seismic data from Greece and the adjacent areas have been used to explain the high seismicity and tectonic activity of this area. Computed 2-D gravity models revealed that below the Aegean region a large “plume” of hot upper-mantle material is rising, causing strong attenuation of the crust. The hot “plume” extends to the base of the lithosphere and has very probably been mobilized through compressional processes that forced the lithosphere to sink into the asthenosphere. The above model is supported by: high heat flow in the Aegean region; low velocity of the compressional waves of 7.7 km/sec for the upper mantle; lower density than normal extending to the base of the lithosphere; teleseismic P-wave travel-time residuals of the order of +2 sec for seismic events recorded at the Greek seismic stations; volcanics in the Aegean area with a chemical composition which can be explained by assuming an assimilation of oceanic crust by the upper mantle; deep seismicity (200 km) which has been interpreted by various authors as a Benioff zone.     

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.     

P-wave velocity structure under the Changbaishan volcanic region, NE China, from wide-angle reflection and refraction data

We present velocity models determined by inverting refracted and reflected arrivals along two active source lines in the Changbaishan volcanic region, NE China. We resolve a prominent low-velocity zone (LVZ) in the crust, with velocities as low as 5.4 km/s. Away from the LVZ, the velocity gradients in the crust are relatively smooth, with average P-wave velocities of about 6.0–6.5 km/s. The Moho is at about 35 km depth, thickening to about 40 km under the Tianchi volcano, and thinning to about 30 km under the LVZ. The LVZ is located about 30–60 km to the north of the summit of the Tianchi volcano (the most recently active volcano in the region), is about 30–75 km in north–south extent, is at most 35 km in east–west extent, and is in the depth range of about 10–25 km below the surface. We use these results to constrain receiver function inversions, and show that the receiver functions in the region are compatible with our findings. With these data alone, the significance of the LVZ in non-unique, although we do not see any evidence to support the presence of partial melt in the crust, and favor the interpretation that the LVZ indicates a residual crustal magma chamber.