Source parameters of the 2004 Kaliningrad earthquakes

An analysis of source parameters of the two unexpected earthquakes from the Kaliningrad (Russia) area is presented. The earthquakes occurred on 21 September 2004 at 11:05:01 and 13:32:31 UT, respectively. The first event was located at the latitude φ = 54.924°N and the longitude λ = 20.120°E, with a depth h = 16 km, and the second event at φ = 54.876°N, λ = 20.120°E and h = 20 km. Magnitudes Mw of the two events were very similar: 5.1 and 5.2. The magnitude values reported by various international data centers have been meaningfully different. The reason is the presence of high-frequency components in Z velocity component of the S wavefield. They were observed along the direction defined by two stations, BLEU in Sweden and SUW in Poland, located in opposite sides of the source. Along the direction perpendicular to it, the effects are relatively very small. The high-frequency waves are understood to mean components in the 6–8 Hz band for event 1 and 2-4 Hz for event 2. The effects in question are also clearly visible on displacement spectrograms. The magnitude values calculated at such stations from S-wave amplitudes or from seismic spectra are clearly overestimated and are close to 6. Therefore, we made a careful selection of channels in order to determine the spectral parameters and, on this basis, the source parameters. The size of the source is relatively small, of about 2 km. The closest seismic station is at 100 source radii from the source. One can clearly see the effect of the TT zone which markedly reduces the seismic moment value for seismic stations laying on the opposite sides of the source. Both events have very similar spatial distributions of the source parameters: magnitude, seismic moment and radius.     

A Source Study of the 9 September 1998 (Mw 5.6) Castelluccio Earthquake in Southern Italy using Teleseismic and Strong Motion Data

Teleseismic and strong motion data are used to derive the source parameters of the September 9, 1998 Castelluccio earthquake (M5.6). Teleseismic body-wave modeling indicates normal faulting (Plane 1: strike = 328°, dip = 50°, rake = −75°; Plane 2: strike = 126°, dip = 42°, rake = −107°) along NW–SE striking planes. Both nodal planes of the computed focal mechanism are tested for their capability of reproducing the mainshock acceleration time histories at three strong motion stations. Synthetic accelerograms are estimated using the stochastic method for finite sources in combination with the H/V ratios technique for the incorporation of the site effect. Our preferred model, which provided the best fit between synthetic and observed waveforms and corresponding Fourier amplitude spectra, consists of a NE-dipping normal fault with dimensions 8 × 7 km. The rupture nucleation point, which is assumed to coincide with the hypocenter location, was confined to the southeasternmost, deepest part of the fault. Our results are in good agreement with the so far released information regarding the aftershock sequence of the examined event, as well as the general seismotectonic knowledge on the broader epicentral area.     

数据域特征波反射走时反演

地震波的运动学信息(走时、斜率等)通常用于宏观速度建模.针对走时反演方法,一个基本问题是走时拾取或反射时差的估计.对于成像域反演方法,可以通过成像道集的剩余深度差近似计算反射波时差.在数据域中,反射地震观测数据是有限频带信号,如果不能准确地确定子波的起跳时间,难以精确地确定反射波的到达时间.另一方面,如果缺乏关于模型的先验信息,则很难精确测量自地下同一个反射界面的观测数据同相轴和模拟数据同相轴之间的时差.针对走时定义及时差测量问题,首先从叠前地震数据的稀疏表达出发,利用特征波场分解方法,提取反射子波并估计局部平面波的入射和出射射线参数.进一步,为了实现自动和稳定的走时拾取,用震相的包络极值对应的时间定义反射波的到达时,实现了立体数据中间的自动生成.理论上讲,利用包络极值定义的走时大于真实的反射波走时,除非观测信号具有无限带宽(即delta脉冲).然而,走时反演的目的是估计中-大尺度的背景速度结构,因此走时误差导致的速度误差仍然在可以接受的误差范围内.利用局部化传播算子及特征波聚焦成像条件将特征波数据直接投影到地下虚拟反射点,提出了一种新的反射时差估计方法.既避免了周期跳跃现象以及串层等可能性,又消除了振幅因素对时差测量的影响.最后,在上述工作基础之上,提出了一种基于特征波场分解的新型全自动反射走时反演方法(CWRTI).通过对泛函梯度的线性化近似,并用全变差正则化方法提取梯度的低波数部分,实现了背景速度迭代反演.在理论上,无需长偏移距观测数据或低频信息、对初始模型依赖性低且计算效率高,可以为后续的全波形反演提供可靠的初始速度模型.理论和实际资料的测试结果证明了本文方法的有效性.     

重叠区域伪谱法计算非均匀地球介质地震波传播

伪谱法是一种高效、高精度计算非均匀介质地震波传播的数值算法,由于它的微分算子的全局性,一般认为该方法不适于并行计算. 本文介绍了并行计算非均匀介质中地震波传播的重叠区域分解算法,给出了一种基于傅里叶伪谱法的并行算法. 文中给出的算法将介质划分为相互重叠的若干区域,在各个子域上单独求解,利用重叠部分的解的传递,将各个子域连接起来,实现了伪谱法在分布式并行处理机上的计算. 文中给出了一个将二维区域分解的算例,比较了并行算法和整体算法的结果,分析了并行算法的计算精度. 结果表明,并行算法会有效降低计算时间,并且保证计算精度. 该方法在大规模三维非均匀介质的地震波场模拟方面有应用价值.      

A review of the influencing factors on teleseismic traveltime tomography

Teleseismic traveltime tomography is an important tool for investigating the crust and mantle structure of the Earth. The imaging quality of teleseismic traveltime tomography is affected by many factors, such as mantle heterogeneities, source uncertainties and random noise. Many previous studies have investigated these factors separately. An integral study of these factors is absent. To provide some guidelines for teleseismic traveltime tomography, we discussed four main influencing factors:the ...     

Application of teleseismic tomography to the study of shallow structure beneath Shizigou in the western Qaidam basin

Teleseismic body wave traveltime tomography is used to inverse the three-dimensional seismic velocity structure beneath Shizigou in the western Qaidam basin.The travel time are picked from the continuous observation data on a small seismic array of stations deployed during 2004-2007.The tomographic results obtained indicate that a NW-trending low velocity anomaly just beneath the target region insert northeastwards with a high dip angle.In the north, northeast and east of the low velocity anomaly, some high...