Observations of teleseismicP wave coda for underground explosions

The earlyP wave coda (5–15 sec after the first arrival) of underground explosions at the Nevada Test Site is studied in the time domain using 2082 teleseismic short-period recordings, with the intent of identifying near-source contributions to the signals in the frequency range 0.2–2.0 Hz. Smaller magnitude events tend to have relatively high coda levels in the 0.4–0.8 Hz frequency band for both Yucca Flat and Pahute Mesa explosions. Coda complexity in this low-frequency passband is negatively correlated with burial depth for Pahute Mesa events but is only weakly correlated with depth for Yucca Flat events. Enhanced excitation of relatively long-period scattered waves for smaller, less deeply buried events is required to explain this behavior. Coda complexity in the 0.8–1.1 Hz band is positively correlated with magnitude and depth for Pahute Mesa events, but has no such dependence for Yucca Flat events. This may result from systematic variations between the spectra of direct signals and coda arrivals caused bypP interference for the largest events, all of which were detonated at Pahute Mesa. Another possible explanation is a frequency-dependent propagation effect on the direct signals of the larger events, most of which were located in the center of the mesa overlying strong lateral velocity gradients in the crust and upper mantle. Event average complexity varies spatially for both test sites, particularly in the 0.8–1.1 Hz band, providing evidence for frequency-dependent focussing or scattering by near-source structure. Both the direct arrivals and the early coda have strong azimuthal amplitude patterns that are produced by defocussing by mantle heterogeneity. The direct arrivals have stronger coherent azimuthal patterns than the early coda for Pahute Mesa events, indicating more pronounced deep crustal and shallow mantle defocussing for the direct signals. However, for Yucca Flat events the direct arrivals have less coherent azimuthal patterns than the coda, suggesting that a highly variable component of near-source scattering preferentially affecting the downgoing energy is superimposed on a pattern produced by mantle heterogeneity that affects the entire signal. This complicated behavior of the direct arrivals may be the result of triplications and caustics produced by the complex basement structure known to underlie the Yucca Flat test site. The presence of strong azimuthal patterns in the early coda indicates that source strength estimates based on early coda are subject to biases similar to those affecting estimates based on direct arrivals.     

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.     

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.     

Seismic tomography constrained by bouguer gravity anomalies: Applications in western Washington

Tomographic inversions for velocity variations in western Washington indicate a high correlation with surface geology and geophysical measurements, including gravity observations. By assuming a simple linear relationship between density and velocity (Birch's law) it is possible to calculate the gravity field predicted from the velocity perturbations obtained by local tomographic inversion. While the predicted gravity matches observations in parts of the model, the overall correlation is not satisfactory. In this paper we suggest a method of constraining the tomographic inversion to fit the gravity observations simultaneously with the seismic travel time data. The method is shown to work well with synthetic data in 3 dimensions where the assumption of Birch's law holds strictly. If the sources of the gravity anomalies are assumed to be spatially localized, integration can be carried out over a relatively small volume below the observation points and sparse matrix techniques can be applied. We have applied the constrained inversion method to western Washington using 4,387 shallow earthquakes, to depths of 40.0 km, (36,865 raypaths) convering a 150×250 km region and found that the gravitational constraints may be satisfied with minor effect on the degree of misfit to the seismic data.     

Spatio-temporal variations in the structure of the attenuation field of the S-wave in the region of Nevada nuclear test site

The characteristics of the attenuation field of short-period shear waves in the region of Nevada nuclear test site (NNTS) are studied. The seismograms of underground nuclear explosions (UNEs) and earthquakes recorded by three seismic stations in 1975–2012 at the epicentral distances of up to 1000 km are processed by the methods based on the analysis of the amplitude ratios of Sn to Pn and Lg to Pg waves, as well as the S-coda envelopes for close events. It is shown that the structure of the attenuation field in the Earth’s crust and upper mantle in the NNTS region experienced significant temporal variations during the interval of nuclear operations. The strongest variations were associated with UNEs conducted in the Pahute Mesa area, which held about two-thirds of the most intense explosions. Our data indicate that temporal variations in the structure of the attenuation field are related to the migration of deep fluids. A comparison of the general characteristics of the attenuation field in the regions of the three large nuclear test sites is presented.     

Three-dimensional Q structure in Jiashi earthquake region of Xinjiang

Introduction It is well known that the occurrence of moderate and strong earthquakes and the distribution of seismically active belts are most directly related to deep crustal faults. Therefore, studying geometrical spreading and physical state of medium of deep crustal faults has always been a major subject in earthquake science. It is not only of great significance to the understanding of deep structural background and seismogenic mechanism of strong earthquakes, but also of great help to se…     

Seismic structure of the upper mantle beneath the western Philippine Sea

We investigate the seismic structure of the western Philippine Sea using two sets of seismological observations: ScS reverberations, which provide the layering framework for a regional upper mantle model, and observations of frequency-dependent phase delays for direct S waves, surface-reflected phases (sS, SS, sSS), and surface waves (R₁, G₁), which constrain the velocity and anisotropy structure within the layers. The combined data set, comprising 17 discontinuity amplitudes and layer travel times from the ScS-reverberation stack and more than 1000 frequency-dependent phase delays, was inverted for a path-averaged, radially anisotropic model. Mineralogical estimates of the bulk sound velocity and density are incorporated as complementary constraints. The final model, PHB3, is characterized by a 11.5-km thick crust, an anisotropic lid bounded by a sharp negative G discontinuity at 89 km, an anisotropic low-velocity layer extending to 166 km, a subjacent high-gradient region, and transition-zone discontinuities at depths of 408 km, 520 km and 664 km. The lid is slower than in a comparable model for the Tonga–Hawaii corridor (PA5), but also significantly thicker, requiring a compositional variation between the two regions. We explore the hypothesis that the thickness of the oceanic lid is controlled by the melting depth at the spreading centers during crust formation, and that the thicker crust and lid in the Philippine Sea results from deeper melting owing to a higher potential temperature and perhaps a higher water content in the upper mantle.     

The effect of spall on Lg waves

Introduction The Lg phase plays a central role in nuclear test detection, discrimination, and yield estima-tion because of its stability in propagation, but the excitation mechanisms of Lg phase are not fully understood. Regional Lg is a short-period guided wave composed mainly of a sequence of multi-ply reflected post-critical S waves trapped in the crustal wave guide. Several studies carried out recently (Gupta et al, 1992; Patton and Taylor, 1995; Gupta et al, 1997; Wallace, 1991) have …     

Three-dimensional Q structure in Jiashi earthquake region of Xinjiang

3-D S-waveQ structure in Jiashi earthquake region is inverted based on the attenuation of seismic waves recorded from earthquakes in this region in 1998 by the Research Center of Exploration Geophysics (RCEG), CSB, and a rough configuration of deep crustal faults in the earthquake region is presented. First, amplitude spectra of S-waves are extracted from 450 carefully-chosen earthquake records, called observed amplitude spectra. Then, after instrumental and site effect correction, theoretical amplitude spectra are made to fit observed amplitude spectra with nonlinear damped least-squares method to get the observed travel time overQ, provided that earthquake sources conform to Brune’s disk dislocation model. Finally, by 3-D ray tracing method, theoretical travel time overQ is made to fit observed travel time overQ with nonlinear damped least-squares method. In the course of fitting, the velocity model, which is obtained by 3-D travel time tomography, remains unchanged, while onlyQ model is modified. When fitting came to the given accuracy, the ultimateQ model is obtained. The result shows that an NE-trending lowQ zone exists at the depths of 10–18 km, and an NW-trending lowQ zone exists at the depths of 12–18 km. These roughly coincide with the NE-trending and the NW-trending low velocity zones revealed by other scientists. The difference is that the lowQ zones have a wider range than the low velocity zones. Foundation item: Joint Seismological Science Foundation of China (957-07-414) and State Key Basic Research Development and Programming Project (95-13-02-02). Contribution No. RCEG200105, Research Center of Exploration Geophysics, China Seismological Bureau.     

层裂对区域震相Lg波的影响

地下核爆炸Lg波的激发机制, 目前尚未十分清楚. 普遍接受的观点是: 伴随层裂源的补偿线性矢量偶极源激发的Rg波的散射形成的S波, 是Lg波的最主要贡献因素. 本文利用频谱比对方法, 结合合成理论地震图方法, 基于东哈萨克斯坦地区速度模型, 分析了东哈萨克斯坦地壳速度模型下的层裂对区域震相Lg波的影响, 并进一步分析了该地区的实际观测资料. 结果表明, 层裂时间函数对Lg波具有扇贝形调制作用; Lg波振幅谱的低谷点源自于补偿线性矢量偶极（CLVD）源激发的Rg波. 这一结果也表明, 东哈萨克斯坦速度模型下的Lg波受层裂的调制作用, 与内华达试验区的Lg波是一样的, 进一步支持了Patton和Taylor的观点, 对于进一步理解及更好地利用Lg波具有重要理论指导意义.      

Theoretical study on the influence of CMB topography on the core reflection ScS

Computing synthetic seismograms for media with localized heterogeneous regions can be performed using hybrid methods. Here, a combination of a finite-difference (FD) technique and a frequency-wavenumber (ω − k) filtering is applied to model wave reflection at different kinds of core-mantle boundary (CMB) topography. The FD method is only applied in the neighbourhood of the CMB, while the ω − k filter is used to continue the reflected wavefield to the Earth's surface. Synthetic SH-seismograms for ScS with a dominant frequency of 0.5 Hz are computed at epicentral distances from 44° to 69°. The topography varies in amplitude (maximum amplitude of 1.0–2.7 km) and in its wavenumber spectrum; it is either monochromatic (wavelengths from 55 to 270 km) or statistical (coloured noise). The seismograms for a CMB with topography are compared with those for a plane CMB. We observe that monochromatic topography with short wavelengths (less than 100 km) results in amplitude reduction and shorter travel times than in the case of a plane CMB, but no variations with epicentral distance appear, whereas greater wavelengths exhibit amplitude variations with distance as well as travel time residuals, which both correlate with the CMB topography. Statistical models show amplitude variations with epicentral distance, while the travel time residuals are very small (less than 0.1 s). All synthetics illustrate that wavefront healing occurs along the ray path from the CMB to the Earth's surface. While the seismograms at the CMB exhibit strong fluctuations, the fluctuations at the surface are smoothed and reduced. This demonstrates that it is necessary to use wave theoretical methods for computing synthetic seismograms for complicated structures at greater depth. It also follows that travel times are less sensitive to the structure than the amplitudes.     

重庆武隆M5.0地震前地磁异常特征分析

选取恩施台、涪陵江东台、石柱黄水台以及武隆仙女山台2012年1月至2017年11月地磁数据进行异常分析,得到各个台站的地磁谐波振幅比和加卸载响应比计算结果,并结合2017年11月23日重庆武隆M5.0地震震中附近地质构造背景分析地震前后地磁异常现象。分析结果认为大部分地磁谐波振幅比变化特征类似于地电阻率的变化特征,表现为"下降-转折-上升"过程。本次地震大致发生在转折上升的初期阶段,其中震中距较小的台站地磁谐波振幅比变化趋势不同步,但变化幅度大致相同,而震中距较大的台站地磁谐波振幅比变化趋势较为同步,变化幅度较大。此外,研究还表明震前地磁加卸载响应比异常高值在时空分布上与震中位置具有很好的对应关系。     

2020年新疆于田MS6.4地震前地磁谐波振幅比变化特征分析

以2020年6月26日于田M_S6.4地震为例, 采用地磁谐波振幅比法, 选用距震中370 km和360 km范围内的且末台、 狮泉河台所记录到的秒值连续波形数据, 对震前地磁异常进行了提取和回溯性分析, 发现地磁谐波振幅比反映了观测点深部电阻率变化, 2020年于田M_S6.4地震前且末台和狮泉河台地磁谐波振幅比出现了较明显的异常现象。 异常特征有: ① 地磁谐波振幅比Y_ZHx和Y_ZHy主要表现为“下降—转折—恢复”的形态, 地震发生在下降型变化恢复之后上升阶段; ② 震前约2年起, 且末台和狮泉河台谐波振幅比表现长短周期变化不同步或NS向与EW向变化不同步的现象; ③ 且末台Y_ZHx(NS)和狮泉河台Y_ZHy(EW)谐波振幅比异常变化中的转折上升形态具有自较长周期向较短周期的迁移现象; ④ 狮泉河台谐波振幅比趋势性异常起始时间早于且末台, 持续时间平均2年左右, 最大异常幅度达到0.04; ⑤ 且末台Y_ZHx(NS)和狮泉河台Y_ZHy(EW)异常相比对应其他分量更为显著, 这也许震中和台站所在位置有关系; ⑥ 地震前和田大地电场(距震中约265 km)岩体裂隙优势方位角(α)也显示异常变化, 可见该地震前多地球物理探测参量的异常响应, 呈现多地球物理量耦合, 增加了异常可信度。     

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

利用在鄂尔多斯块体内部布设的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波速度结构

采用作者最近新提出的频谱参数化方法,对北京电信传输台网2447个远震P波到时残差数据进行了反演,说明新的参数化方法能有效地应用于实际数据.反演结果表明,京津唐张地区深部100多km内存在明显的速度横向不均匀性;50 km以上部分的结构与地表的地质构造有较好的对应关系:在燕山隆起区下部表现为高速异常区,而凹陷区下部则呈现低速异常区;深度100 km以下,地表构造的痕迹逐步消失,异常蝠度也相应降低.研究区域范围内存在三个大的低速异常块体,其中以天津、唐山附近的水平尺度为100 km左右的低速异常区幅度最大,向下一直延伸到约100 km的深度,而且在100 km以下存在一低速异常构造与之对应.研究区域内发生的强震大部分都位于低速异常区边缘向高速区过渡的梯度带上.     

Linearized rays,amplitude and inversion

In this paper, ray theoretical amplitudes and travel times are calculated in slightly perturbed velocity models using perturbation analysis. Also, test inversions using travel time and amplitude are computed. The pertubation method is tested using a 3-D velocity model for NORSAR having velocity variations up to 8.0 percent. The perturbed amplitudes are found to be in excellent agreement with the calculated ray amplitudes. Velocity inversions based on travel time and amplitude are next investigated. Perturbation analysis using linearized ray equations is efficiently used to compute amplitude derivatives with respect to model parameters. The trial linearized inversions use smaller velocity variations of 1.7 percent to avoid possible effects due to ray shift, even though the perturbation analysis is valid for larger variations. The trial 2-D inversion results show that linearized amplitude inversions are complementary and not redundant to travel time inversions, even in smoothly varying models.     

吉林松原M≥5.0震群前地磁谐波振幅比异常特征

以东北地区(40°—50°N,120°—130°E)为研究区域,利用黑龙江、吉林和辽宁地区地磁观测数据,进行谐波振幅比计算,对2013年10月31日至11月23日吉林松原接连发生的5个M≥5.0震群前该研究区地下电性结构进行分析,发现:(1)在该区域内显著地震发生前,望奎、通河、三岗、通化和铁岭地震台均出现下降—转折—恢复的异常变化特征,持续时间1—4.5年,变化幅度为0.04—0.07,而通河台Y_(ZH_y)项的不同周期间存在不同步现象,且异常持续2年以上;(2)通河台分钟值谐波振幅比具有由深至浅的周期迁移特征,三岗台秒值数据结果 Y_(ZH_y)向具有由浅入深的迁移特征;(3)2014年3月至2016年4月通河台水平方向的EW向数据不同深度介质电性出现差异性异常,可能预示着该台周边具有发生中强地震的危险,应重点关注该异常变化。     

2017年8月8日九寨沟MS7.0地震前成都台地磁谐波振幅比异常分析

2017年8月8日四川省九寨沟县发生M_S7.0地震,成都地震基准台距此次地震震中约255km,震前该台地磁谐波振幅比出现趋势性异常。对成都台GM4磁通门磁力仪秒数据进行谐波振幅比计算,结果显示,谐波振幅比在九寨沟M_S7.0地震前表现出下降—转折—恢复上升的异常变化形态,地震发生在异常恢复期,同时,异常表现出由长周期向短周期迁移、SN向与EW向变化不同步的特征,芦山M_S7.0地震前成都台地磁谐波振幅比亦呈现了相似的异常变化特征。     

2017年11月23日重庆武隆5.0级地震前地磁谐波振幅比异常研究

近年来,地磁谐波振幅比方法在国内得到较广泛的应用,并成为一种较为有效的地震预测方法。本文利用谐波振幅比方法,对2017年11月23日重庆武隆5.0级地震震中附近台站的地磁三分量数据进行计算,并分析了各台站10~60min周期谐波振幅比值异常特征。研究发现,异常台站几乎均位于震源机制解的拉张区(P波向下的区域),具有象限分布特征;异常台站NS向谐波振幅比值变化存在分层现象,短周期(浅部)的变化与长周期(深部)变化成相反相位,从电磁学机理看,浅部与深部之间存在EW向的面电流;近震中台站的异常低点时间存在由长周期(深部)向短周期(浅部)迁移的现象。