Scattering of surface waves modelled by the integral equation method

The integral equation method is used to model the propagation of surface waves in 3-D structures. The wavefield is represented by the Fredholm integral equation, and the scattered surface waves are calculated by solving the integral equation numerically. The integration of the Green's function elements is given analytically by treating the singularity of the Hankel function at   R = 0  , based on the proper expression of the Green's function and the addition theorem of the Hankel function. No far-field and Born approximation is made. We investigate the scattering of surface waves propagating in layered reference models imbedding a heterogeneity with different density, as well as Lamé constant contrasts, both in frequency and time domains, for incident plane waves and point sources.     

Anomalous high-frequency wave propagation from the Tonga–Kermadec seismic zone to New Zealand

Summary. Bulletins of the International Seismological Centre (ISC) show very large residuals, up to 15 s early, for arrivals from events in the Tonga–Kermadec subduction zone to the New Zealand network of seismometers. The very early arrivals are confined to events south of about 22°S, and shallower than about 350 km. The waveforms show two distinct phases: an early, emergent, first phase with energy in the high-frequency band 2–10 Hz, and a distinct second phase, containing lower frequency energy, arriving at about the time predicted by JB tables.
The residuals are attributed to propagation through the cold, subducted lithosphere, which has a seismic velocity 5 per cent faster, on average, than normal. Ray tracing shows that the ray paths lie very close to the slab for events south of 22°S, but pass well beneath the slab for events further north, corresponding to the change in residual pattern. This characteristic of the ray paths is due to the curved shape of the seismic zone, and in particular to the bend in the zone where the Louisville ridge intersects the trench at 25°S.
The residuals can only be explained if the high velocity anomaly extends to a depth of 450 km in the region of the gap in deep seismicity from 32 to 36°S. The very high-frequency character of the first phase requires the path from the bottom of the slab to the stations to be of high Q , and to transmit 2–10 Hz energy with little attenuation.
The absence of low-frequency energy in the first phase is due to the narrowness of the high-velocity slab, which transmits only short-wavelength waves. The second phase, which contains low frequencies, is identified as a P -wave travelling beneath the subducted slab in normal mantle. There is no need to invoke any special structures, such as low-velocity waveguides or reflectors, to explain any of the observations. The S -wave arrivals show similar effects.     

Intrusive Mechanisms at Mt. Etna Forerunning the July-August 2001 Eruption from Seismic and Ground Deformation Data

In this work we present seismological and ground deformation evidence for the phase preparing the July 18 to August 9, 2001 flank eruption at Etna. The analysis performed, through data from the permanent seismic and ground deformation networks, highlighted a strong relationship between seismic strain release at depth and surface deformation. This joint analysis provided strong constraints on the magma rising mechanisms. We show that in the last ten years, after the 1991–1993 eruption, an overall accumulation of tension has affected the volcano. Then we investigate the months preceding the 2001 eruption. In particular, we analyse the strong seismic swarm on April 20–24, 2001, comprising more than 200 events (M_max = 3.6) with prevalent dextral shear fault mechanisms in the western flank. The swarm showed a ca. NE-SW earthquake alignment which, in agreement with previous cases, can be interpreted as the response of the medium to an intrusive process along the approximately NNW-SSE volcano-genetic trend. These mechanisms, leading to the July 18 to August 9, 2001 flank eruption, are analogous to ones observed some months before the 1991–1993 flank eruption and, more recently, in January 1998 before the February-November 1999 summit eruption.     

山西洪洞ML5.0地震前测震学异常指标分析

对山西洪洞2003年ML5．0地震进行了分析，发现震前1a左右，在临汾盆地周边地区相继出现断层总面积、η值、A(b)值、D值、5值、YH值等地震活动性指标异常，这些指标异常反映了5．0级地震前地震活动呈现增强的特征。根据非线性理论中的尺度依赖性原理提出一种对特定范围的中强地震进行检验性预测的方法。     

近期国际地震预测预报研究进展的几个侧面

本文从国际组织、国际计划、科学出版物、主要代表性人物及其学术观点、实际地震预测与地震前兆检验等方面，介绍了20世纪90年代后期以来国际上地震预测预报研究进展的几个侧面，这些进展在一定程度上反映了此间地震预测预报研究从“低谷”走向新的发展的思想脉络。     

再论各向异性介质转换波地震学,第二部分:参数计算研究

在具有垂直对称轴横向各向同性介质中，利用四种参数来确定中间至远偏移距转换波（C-波）动校正。它们是C-波叠加速度Vc2，垂直速度比和有效速度比γ0和γeff以及各向异性参数χeff。我们将这四种参数作为C波叠加速度模型。C-波速度分析的目的就是确定这种叠加速度模型。C-波叠加速度模型Vc2，γ0，γeff，和χeff可以由P-波和C-波反射动校正资料获得。然而错误的传播是C-波反射动校正反演中的严重问题。当前短排列叠加速度由于是从双曲线动校正推算而得，因而其精度不足以为各向异性参数提供有意义的反演值。中间偏移非双曲线动校正不再被人们所勿略，而是可以用一个背景γ加以量化。非双曲线分析通过中间偏移距的γ校正量可以产生Vc2，若数据不含燥音，其误差小于1％。方法稳健，允许γ启始假定值的误差达20％。该方法也适用垂直非均匀各向异性介质。精度的提高使能够用4分量地震资料计算各向异性参数。为此提出了两种工作流程：双扫描和单扫描流程。理论数据和实际数据的应用表明这两种流程得出的结果其精度相似，但是单扫描流程比双扫描更有效。     

可控源音频大地电磁测量法在砂岩型铀成矿地质环境研究中的应用

本文简要介绍了可控源音频大地电磁法测量的基本原理、数据处理方法,并结合新疆伊犁、内蒙古等地区的实测成果,论述了该方法在砂岩型铀成矿地质环境研究中的应用的效果。