On the relation of variations in the travel times of seismic waves with the changes in the speed of the Earth’s rotation

The results of long-term sounding of the lithosphere by seismic waves from the deep-focus Hindu Kush earthquakes are presented. The travel time series of the first longitudinal wave on a fixed base are constructed for six seismic observation stations (SS) located on the Russian Platform (the Obninsk SS), on the Siberian Platform (the Eltsovka SS), on the Cis-Ural Trough (the Arti SS), in the Central-Ural Megazone (the Sverdlovsk SS), in the Transbaikalia (the Bodaibo SS), and in the Northern Tien Shan (the Przhevalsk SS). The time series duration in years for these stations ranged from 1964–1970 to 2007. The travel time series of seismic waves for the stations indicated are characterized by multi-slope negative linear trends caused by changes in the stress-strain state of rocks. From the comparison of the trend slopes at different stations it follows that the changes in the stress conditions within the lithosphere are relatively weak in the aseismic regions of the Russian and Siberian Platforms and in the Ural Megazone, whereas in the seismically active regions of Tien Shan, Transbaikalia and the Cis-Ural Trough they are more pronounced. The correlation has been observed between the time series trends of the average annual travel times of seismic waves and the time series of the Earth’s rotation speed. The strongest correlation between the series can be seen for the stations, located on the platforms with weak manifestations of both seismicity and active geodynamic processes. Within the long-term periods of deceleration and acceleration of the Earth’s rotation, travel times of seismic waves are decreased and increased, respectively.     

利用瑞利面波相速度和方位各向异性研究鄂尔多斯块体的岩石圈变形特征

利用在鄂尔多斯块体内部布设的45个宽频带流动台站和固定台站的资料,用双平面波方法反演了20~143 s共12个周期的基阶瑞利面波的平均相速度和方位各向异性,并反演了一维S波速度结构.反演结果显示50~100 s中长周期的瑞利面波相速度高于AK135速度模型的相速度,为高速异常,S波速度显示高速异常主要位于180 km深度范围内,表明鄂尔多斯块体保留有厚的高速岩石圈.20~111 s周期的方位各向异性强度小于1%,较小的各向异性表明鄂尔多斯块体岩石圈变形较弱.20~50 s周期的平均快波方向为近EW向,67~143 s周期的平均快波方向为NW-SE向,相对发生了整体改变,快波方向的转变约开始于80~100 km深度范围,这表明岩石圈上下部存在着由不同变形机制导致的各向异性.上部岩石圈中各向异性可能主要为残留的“化石”各向异性,而下部岩石圈各向异性可能是现今板块构造运动导致的变形而形成.鄂尔多斯块体岩石圈垂向上的变形差异可能主要与岩石圈温度随深度的变化以及青藏高原NE-NNE向挤压引起的上部岩石圈逆时针旋转有关.     

P wave amplitude measurements at the Kaapvaal network: do they help in resolving deep earth structure?

Measurements of the amplitudes of seismic body waves at teleseismic distances have rarely resolved significant features mainly because of the large scatter of the data. However, amplitudes are easy to measure and may provide additional constraints on structure to supplement times and waveforms. A new approach to analysing body wave amplitudes at a regional network of similar instruments seeks to minimize scatter by first deriving amplitude station corrections analogous to station corrections for times. After correction for station effects, amplitudes from several events can be combined to give regional amplitude–distance curves without using information on event magnitudes. However, the earthquakes providing the observations must lie in a restricted range of azimuths from the stations of the network and provide considerable overlap in the range of distances between adjacent events, with no gaps in distance coverage. The advantages of the method are explored using P wave amplitudes from two sets of earthquakes in the Indonesian and South American regions recorded by the Kaapvaal network deployed across southern Africa. In the first example, high amplitudes near 88° distance suggest the presence of a small discontinuity at the top of D″ that causes constructive interference between the closely separated arrivals of a small triplication in the travel times. The second example, supplemented by calculations using synthetic data, shows how long-wavelength regional variations in amplitudes can be resolved to assist the interpretation of times and waveforms. However, the limited range of distances in the observations and lateral heterogeneities at any depths can result in bias or tilt of the amplitude–distance relationships. Constraining the depths of the structure causing the long-wavelength variations is a subject for future research.     

Travel-time residuals and three-dimensional velocity structure of Italy

Jeffreys-Bullen P and PKP travel-time residuals observed at more than 50 seismic stations distributed along Italy and surrounding areas in the time interval 1962–1979, indicate the complex velocity pattern of this region. Strong lateral velocity inhomogeneities and low velocity zones are required to explain the observed pattern of residuals. In particular, late arrivals of about 1 sec are observed in the Apenninic mountain range, requiring both greater crustal thickness and low velocity layers, coherent with seismic refraction data and surface wave dispersion measurements. The seismic stations located in the Western and Eastern Alps indicate the presence of high velocities. In the Western Alps the strong azimuthal variation of residuals and the high values of early arrivals have a close relationship to the Ivrea body, an intrusive crustal complex characterized by a velocity as high as 7–7.2 km/sec.A travel-time inversion performed with theAki et al. (1977) block model, confirms the peculiar characteristics and the sharp variations in the lithosphere of the whole Italian region, with values of velocity perturbations between many adjacent blocks, ranging in size from 50 to 100 km, and independent from the earth parametrization chosen, reaching values up to 10% in the lithospheric part and 5% in the asthenosphere. 3-D inversion requires also high velocity along the Tyrrhenian coastal margin, equivalent to an uprise of major crustal and lithospheric discontinuities along this part of the Italian peninsula. Moreover low velocity material must be present in the northern part of the Adriatic foreland, in the lithosphere-asthenosphere system, closely related to the stress and seismicity pattern, and the lateral bending of the lithosphere in the same region.     

Heterogeneities in the field of short period seismic wave attenuation in the lithosphere of central Tien Shan

Records of deep-focus Hindu Kush earthquakes in the depth ranges 70–110 and 190–230 km made by 45 digital and analogue seismic stations were analyzed to study the attenuation field of short period seismic waves in the lithosphere of central Tien Shan. The dynamic characteristics studied include the ratio of peak amplitudes in S and P waves (S/P) and the ratio of the S-wave maximum to the coda level in the range t = 400 ± 5 s, where t is the lapse time (S/c400) for 1.25 Hz. Comparatively high values of S/P are shown to prevail in most of the area, corresponding to lower S-wave attenuation. Upon this background is a band of high and intermediate attenuation in the west of the area extending along the Talas-Fergana fault in the south and afterwards turning north-northeast. The rupture areas of the two largest (M ≥ 7.0) earthquakes which have occurred in Tien Shan during the last 25 years are confined to this band. Abnormally high values of S/c400 were obtained for stations situated in the rupture zone of the August 19, 1992, magnitude 7.3 Suusamyr earthquake and around it. For two of the stations we found considerable time variations in the coda envelope before the earthquake. The effective Q was derived from compressional and shear wave data for the entire area, as well as for the band of high attenuation. Comparison with previous data shows that the attenuation field in the area has changed appreciably during 20–25 years, which can only be due to a rearrangement of the fluid field in the crust and uppermost mantle. It is hypothesized that a large earthquake is very likely to occur in the northern part of the attenuating band.     

CORRELATION PATTERNS IN CHARACTERISTICS OF SPATIALLY VARIABLE SEISMIC GROUND MOTIONS

A methodology for the investigation of the spatial variation of seismic ground motions is presented; data recorded at the SMART-1 dense instrument array in Lotung, Taiwan, during Events 5 and 39 are used in the analysis. The seismic motions are modeled as superpositions of sinusoidal functions, described by their amplitude, frequency, wavenumber and phase. For each event and direction (horizontal or vertical) analysed, the approach identifies a coherent, common component in the seismic motions at all recording stations, and variabilities in amplitudes and phases around the common component sinusoidal characteristics, that are particular for each recording station. It is shown that the variations in both the amplitudes and the phases of the motions at the station locations around the common component characteristics contribute significantly to the spatially variable nature of the motions, and, furthermore, they are correlated: increase in the variability of the amplitudes of the motions recorded at individual stations around the common amplitude implies increase in the variability of the phases around the common phase. The dispersion range of the amplitude and phase variability around their corresponding common components appear also to be associated with physical parameters. The spatially variable arrival time delays of the waveforms at the stations due to their upward travelling through the site topography, in addition to the wave passage delays identified from signal processing techniques, constitute another important cause for the spatial variation of the motions; their consideration in the approach facilitates also the identification of the correlation patterns in the amplitudes and phases. © 1997 by John Wiley & Sons, Ltd.     

Variations in short-period shear-wave attenuation in the Baikal Rift Zone and their relationship to seismicity

This paper presents results from a study of variations in short-period shear-wave attenuation in the lithosphere of the Baikal Rift Zone (BRZ). We used earthquake records made at the Ulan-Bator station (ULN) at epicentral distances of ～400–1300 km. The ratios of maximum amplitudes in the Sn and Pn waves were considered. We show that these ratios are on the whole considerably higher than those in other areas of Central Asia. It was found that zones of low shear wave attenuation coincide with the rupture zones of large earthquakes that occurred during the 19th and 20th centuries. We identified zones of high attenuation where no large (M～ > 7.0) seismic events have occurred during at least 180 years. The hypothesis we propose is that precursory processes before future large earthquakes may be occurring in these zones. We discuss the question of whether wave attenuation characteristics may be related to seismicity.     

青藏高原中部INDEPTH-Ⅲ剖面上地幔远震S波层析成像

印度板块向欧亚俯冲前缘位于班公—怒江缝合带附近,但是印度岩石圈地幔的俯冲形态和形变过程仍然缺乏共识,在不同地区使用不同方法获得的结果之间存在明显差异.本文使用青藏高原中部INDEPTH-Ⅲ剖面远震S波波形数据,提取走时信息,通过层析成像方法获得剖面下方S波速度扰动图像.结果显示:在班公—怒江缝合带下方100至300km深度范围内存在一个高角度(约65°)北倾的S波高速体,推测可能是回退的印度岩石圈板片或/和小规模对流引起的岩石圈拆沉后残留的印度大陆岩石圈板片.     

T waves from deep earthquakes generated exactly at the bottom of deep-sea trenches

T waves (seismic water waves), which were generated by deep-focused earthquakes, have been found by an array of sensitive ocean-bottom seismographic observations depolyed on the western Pacific basin. The points of generation of T waves have been exactly located by use of the accurate velocity of water waves which were known from explosions. The positions obtained are at the bottom of deep-sea trenches; however, the positions are slightly (10–35 km) ocean-side of the trench. T waves have been known to be generated by seismic waves which were transmitted from the focus to the trench bottom along the descending lithosphere. The intensity of the observed T waves implies that the Q value along the descending lithosphere is more than 4000. The positions of T-wave generation are consistent with the 8.2- to 8.6-km/s stratified structure of the oceanic lithosphere. T waves from shallow earthquakes beneath the lower continental slope are also clearly observed by bottom seismography.     

Seismological evidence for a lithospheric normal faulting — the Sanriku earthquake of 1933

The focal process of the Sanriku earthquake of March 2, 1933, is discussed in relation to the bending mechanism of the lithosphere. On the basis of the P times obtained at more than 200 stations, it is confirmed that the hypocenter of this earthquake is within the lithosphere beneath the Japan trench. The P wave fault plane solution, the amplitude of long-period (100 s) Love and Rayleigh waves and two near-field observations suggest, almost definitely, that the Sanriku earthquake represents a predominantly normal faulting on a plane dipping 45° towards N 90° W. A fault size of 185 × 100 km², in agreement with the size of the aftershock area, is required to yield a slip dislocation of 3.3 m, a value consistent with the tsunami data. This result suggests that the fracture took place over the entire thickness of the lithosphere, thereby precluding the possibility that the Sanriku earthquake merely represents a surface tensile crack due to the bending of the lithosphere. This large scale lithospheric faulting is presumably due to a gravitational pull exerted by the cold sinking lithosphere. The fracture probably took place on an old fault plane which had once fractured and healed up. The existence of this fracture zone which decouples, to some extent, the oceanic lithosphere from the sinking lithosphere accounts for the sharp bend of the lithosphere beneath oceanic trenches and also the abrupt disappearance of seismic activity across oceanic trenches. The sharp bend of the lithosphere is therefore a result, not the cause, of great earthquakes beneath oceanic trenches.     

Broadband frequency-dependent amplification of seismic waves across Bucharest, Romania

The determination of seismic amplitude amplification is a fundamental contribution to seismic hazard assessment. While often only high-frequency amplitude variations (>1 Hz) are taken into account, we analyse broadband waveforms from 0.14 to 8.6 Hz using a temporary network of 32 stations in and around the earthquake-prone city of Bucharest. Spectral amplitudes are calculated with an adaptive multiple-taper approach. Across our network (aperture 25 km × 25 km), we find a systematic northwest/southeast-oriented structural influence on teleseismic P-wave amplitudes from 0.14 to 0.86 Hz that can be explained by constructive interference in the dipping Cenozoic sedimentary layers. For higher frequencies (1.4–8.75 Hz), more local site effects prevail and can be correlated partly among neighbouring stations. The transition between systematic and localised amplitude variations occurs at about 1 Hz.     

Stress-strain state of the lithosphere in the southern Baikal region and northern Mongolia from data on seismic moments of earthquakes

Investigation and understanding of the present-day geodynamic situation are of key importance for the elucidation of the laws and evolution of the seismic process in a seismically active region. In this work, seismic moments of nearly 26000 earthquakes with K _p ≥ 7 (M _LH ≥ 2) that occurred in the southern Baikal region and northern Mongolia (SBNM) (48°–54°N, 96°–108°E) from 1968 through 1994 are determined from amplitudes and periods of maximum displacements in transverse body waves. The resulting set of seismic moments is used for spatial-temporal analysis of the stress-strain state of the SBNM lithosphere. The stress fields of the Baikal rift and the India-Asia collision zone are supposed to interact in the region studied. Since the seismic moment of a tectonic earthquake depends on the type of motion in the source, seismic moments and focal mechanisms of earthquakes belonging to four long-term aftershock and swarm clusters of shocks in the Baikal region were used to “calibrate” average seismic moments in accordance with the source faulting type. The study showed that the stress-strain state of the SBNM lithosphere is spatially inhomogeneous and nonstationary. A space-time discrepancy is observed in the formation of faulting types in sources of weak (K _p = 7 and 8) and stronger (K _p ≥ 9) earthquakes. This discrepancy is interpreted in terms of rock fracture at various hierarchical levels of ruptures on differently oriented general, regional, and local faults. A gradual increase and an abrupt, nearly pulsed, decrease in the vertical component of the stress field S _v is a characteristic feature of time variations. The zones where the stress S _v prevails are localized at “singular points” of the lithosphere. Shocks of various energy classes in these zones are dominated by the normal-fault slip mechanism. For earthquakes with K _p = 9, the source faulting changes with depth from the strike-slip type to the normal-strike-slip and normal types, suggesting an increase in S _v . On the whole, the results of this study are well consistent with the synergism of open unstable dissipative systems and are usable for interpreting the main observable variations in the stress-strain state of the lithosphere in terms of spatiotemporal variations in the vertical component of the stress field S _v . This suggests the influence of rifting on the present-day geodynamic processes in the SBNM lithosphere.     

Creepex of Kamchatka shallow earthquakes

Results of investigation of the lithosphere in the Kamchatka seismic focal zone from dynamic characteristics of earthquake records obtained at regional stations are presented. It is assumed that the specificity of the source zone can be estimated by the relation Cr = K _P ? bK _S ? c characterizing relative energies (energy classes, according to [Fedotov, 1972]) of short period transverse and longitudinal waves in the source. Azimuthal, spatial, and temporal variations in Cr and their relation to focal mechanisms are examined. Spatiotemporal variations in this parameter are shown to be caused by the influence of variations in the conditions in the source zone (its substance or process) on the radiation of P and S waves.     

Temporary Array Data for Studying Seismic Anisotropy of Variscan Massifs – The Armorican Massif,French Massif Central And Bohemian Massif

We present the first results of a comparison of deep lithosphere structure of three Variscan massifs - the Armorican Massif, French Massif Central and Bohemian Massif, as revealed by recent tomographic studies of seismic anisotropy. The data originate from several field measurements made in temporary arrays of stations equipped with both short-period and broadband seismometers with digital recording. The study is based on teleseismic body waves and a joint inversion of anisotropic data (P-residual spheres, the fast shear-wave polarizations and split times) and demonstrates that the three Variscan massifs appear to consist of at least two parts with different orientation of large-scale fabric derived from seismic anisotropy. The boundaries of anisotropic lithospheric domains are related to prominent tectonic features recognised on the surface as sutures, shear zones or transfer fault zones, as well as grabens, thus indicating that some of them extend deep through the entire lithosphere.     

Shallow geology characterization using Rayleigh and Love wave dispersion curves derived from seismic noise array measurements

The local geology and shallow S-wave velocity structure of a site are recognized to be key factors for the increase in the damaging potential of seismic waves. Indeed, seismic amplitudes may be amplified in frequency ranges unfavorable for building stock by the presence of soft sedimentary covers over lying hard bedrock. Hence, microzonation activities, which aim at assessing the site response as accurately as possible, have become a fundamental task for the seismic risk reduction of urbanized areas. Methods based on the measurement of seismic noise, which typically are fast, non-invasive, and low cost, have become a very attractive option in microzonation studies.Using observations derived from seismic noise recordings collected by two-dimensional arrays of seismic stations, we present a novel joint inversion scheme for surface wave curves. In particular, the Love wave, the Rayleigh wave dispersion and the HVSR curves are innovatively combined in a joint inversion procedure carried out following a global search approach (i.e., the Genetic Algorithm).The procedure is tested using a data set of seismic noise recordings collected at the Bevagna (Italy) test-site. The results of the novel inversion scheme are compared with the inversion scheme proposed by Parolai et al. (2005), where only Rayleigh wave dispersion and HVSR curves are used, and with a cross-hole survey.     

A discussion of world-wide measurements of tidal gravity with respect to oceanic interactions,lithosphere heterogeneities,Earth's flattening and inertial forces

By the end of 1981 the International Center of Earth Tides (ICET) had collected and evaluated a considerable amount of data from 180 stations, including those of the Trans-World Profile which ensure for the first time a world-wide distribution including the tropical areas and the southern hemisphere. In 1979–80, new oceanic cotidal maps of high quality, established by Schwiderski, became available. These maps fit “on-land” tidal-gravity measurements quite successfully. A new theoretical approach developed by Wahr in 1981 has resulted in a set of theoretical formulae establishing the latitude-dependence of the classical elastic amplitude factors for tidal deformations.We calculate here, for six tidal waves, the correlations between the observed gravity variations and those resulting from a calculation based upon the Schwiderski maps. These correlations are highly significant. After subtraction of these oceanic effects we calculate the latitude-dependences of the experimentally determined amplitude factors, which are found to fit Wahr's theoretical formulae. There remains, however, a serious discrepancy in the constant terms of the various formulae. The effects of heterogeneities in the lithosphere on tidal deformations are also clearly identified.     

Geodynamic processes in seismically active areas of the Tien Shan from monitoring data with the use of nuclear explosions

Recent geodynamic processes in the Tien Shan region are studied by the analysis of time series of effective velocities and traveltime delays relative to the IASPEI-91 traveltime curve of the weakly refracted wave P _n from nuclear explosions at the Semipalatinsk test site over the period of 1968–1989. The time series were constructed for 10 seismic stations located at distances of 800–1200 km from the test site in the regions of the Northern, Central, and Southern Tien Shan. The twenty-year period of observations at stations in the North Tien Shan showed a significant decrease in traveltime delays by 0.20–0.76 s, which corresponds to a 0.2–0.7% increase in seismic velocities. An opposite pattern is observed at stations of the Central and Southern Tien Shan: traveltime delays increased by 0.2–0.5 s and, accordingly, seismic velocities dropped by 0.2–0.5%. These results suggest the predominance of compression processes in the crust and upper mantle during the period of observations in Northern Tien Shan and extension processes in the Central and Southern Tien Shan. The series of velocities and traveltime residuals are characterized by the presence of rhythmic oscillations of various amplitudes and periods against a linear trend. A correlation between variations in kinematic parameters and yearly numbers of earthquakes is observed at all stations. Diagrams of the spectral time analysis reveal rhythms with periods of 2–3 and 5–7 yr. The data obtained in this study are consistent with results of studying the stress-strain state of the Tien Shan crust from focal mechanisms of earthquakes and the velocities of recent crustal movements from GPS data. It is found that the amplitude of variations in kinematic parameters of the P wave at stations located in seismically active regions (the Tien Shan, Kopet Dagh, the Caucasus, Altai, and Sayany) is two to five times higher compared to aseismic regions (the Russian and Kazakh plains).     

Seismic Anisotropy and Velocity Variations in the Mantle beneath the Saxothuringicum-Moldanubicum Contact in Central Europe

—We report on results of a passive seismic experiment undertaken to study the 3-D velocity structure and anisotropy of the upper mantle around the contact zone of the Saxothuringicum and Moldanubicum in the western margin of the Bohemian Massif in central Europe. Spatial variations of P-wave velocities and lateral variations of the particle motion of split shear waves over the region monitor changes of structure and anisotropy within the deep lithosphere and the asthenosphere. A joint interpretation of P-residual spheres and shear-wave splitting results in an anisotropic model of the lithosphere with high velocities plunging divergently from the contact of both tectonic units. Lateral variations of the mean residuals are related to a southward thickening of the lithosphere beneath the Moldanubicum.     

Seismomagnetic response of a fault zone

Based on the results of instrumental observations of geomagnetic variations caused by the propagation of seismic waves through a fault zone, the dependences between the amplitudes of the induced seismomagnetic effect and seismic signal as a function of distance r to the midline of the fault are obtained. For the first time, it is shown that the amplitude of the seismomagnetic effect is maximal in the fault damage zone. The phenomenological model describing the generation of magnetic signals by seismic waves propagating through the crushed rock in the tectonic fault zone is suggested. It is assumed that geomagnetic variations are generated by the changes in the electrical conductivity of the fragmented rocks as a result of the deformation of the rock pieces contacts. The amplitudes of the geomagnetic variations calculated from the model agree with the instrumental observations.     

Stacking of records of nuclear explosions registered by the world network of seismic stations

The methods and results of stacking of seismograms of nuclear explosions registered by the world network of seismic stations from 1980 until 1988 are presented. Stacking of seismograms obtained at one seismic station from the sources located in the same place allows us substantially to increase the useful signal/noise ratio. This, in turn, allows more confident correlation of various seismic waves coming to the seismic stations after the first arrivals, particularly, converted and reflected waves from different boundaries inside the Earth.