Implication of seismic noise for determining the structure of the upper Earth Rock Mass

An assertion that the cross-correlation function of seismic noise, considered as a result of the superposition of the surface waves, excited by the sources, randomly distributed over the Earth’s surface, determines the Green function of the surface wave is verified by numerical modeling. The maximum wave periods, for which this assertion is correct, are estimated and the errors in determination of the phase and group velocity of the surface waves are evaluated. The procedure for the determination of the correlation function and estimation from it of the group velocity are tested thoroughly based on the example of the pair of the BJT and TLY stations in Asia. This procedure is used for obtaining the group and phase velocities of the Rayleigh waves on the traces between the OBN-ARU and PUL-ARU stations. The velocity sections of the transverse waves are built based on the dispersion curves of the phase and group velocities of averages along these traces. The region of the lowered velocity in the upper mantle at depths of 150–300 km is revealed on both traces. From the analysis of correlation functions, which are subjected to narrow-band filtering, it is shown that the frequency composition of noise varies from the East and from the West from the profiles between the stations: in the East (Siberia) the noise has an appreciably lower-frequency than in the West (Western Europe).     

A new stacking technique in ambient noise tomography

We have developed a new stacking technique in ambient noise tomography to obtain high-quality dispersion curves of Rayleigh waves.This technique is used to stack the vertical components of the Estimated Green Functions(EGFs) obtained respectively from cross correlation of the ambient noise data recorded by a remote seismic station and one of the short distance seismic stations of a seismic array.It is based on a phase-matched filter and is implemented by a four-step iterative process:signal compression,stacking,signal extraction and signal decompression.The iterative process ends and gives the dispersion curve of Rayleigh wave when the predicted one and the processing result converge.We have tested the method using the vertical components of synthetic Rayleigh wave records.Results show that this new stacking method is stable and it can improve the quality of dispersion curves.In addition,we have applied this method to real data.We see that the results given by our new technique are obviously better than the ones employing the traditional method which is a three-step process:signal compression,signal extraction and signal decompression.In conclusion,the new method proposed in this paper can improve the signal to noise ratio of EGFs,and can therefore potentially improve the resolution of ambient noise tomography.     

Phase and Correlation in `Random' Seismic Fields and the Reconstruction of the Green Function

We first present a summary of recent results on coda interpretation. We emphasize the observation of the stabilization of P to S energy ratio indicating the modal equipartition of the wavefield. This property clearly shows that the coda waves are in the regime of multiple scattering. Numerical solutions of the elastic radiative transfer equation are used to illustrate the evolution of the wave-field towards P-to-S energy stabilization, and asymptotically to complete isotropy. The energy properties of the coda have been widely studied but the phase properties have often been neglected. The recently observed coherent backscattering enhancement, an expression of the so-called `weak localization', demonstrates that interference effects still persist for multiple diffracted waves. Another manifestation of the persistence of the phase is the possibility to reconstruct the Green function between two stations by averaging the cross correlation of coda waves produced by distant earthquakes and recorded at those two stations. This reconstruction is directly related to the properties of reciprocity and time reversal of any wavefield. Using broadband seismic coda waves, we show that the dominant phases of the Green function in the band 2 s–10 s, namely fundamental mode Rayleigh and Love waves, are reconstructed. We analyze the time symmetry of the cross correlation and show how the level of symmetry evolves with the isotropization of the diffuse field with lapse time. Similarly we investigate the correlation in continuous ambient noise records. Whereas the randomness of the coda results from multiple scattering by randomly distributed scatterers, we assume that the seismic noise is random mostly because of the distribution of sources at the surface of the Earth. Surface waves can be extracted from long time series. The dispersion curves of Rayleigh waves are deduced from the correlations. On paths where measurements from earthquake data are also available, we show that they are in good agreement with those deduced from noise correlation. The measurement of velocities from correlation of noise along paths crossing different crustal structures opens the way for a `passive imaging' of the Earth's structure.     

利用背景噪声数据提取地震台站间 面波的可靠性分析 以中国大陆中东部地区的宽频带台站为例

利用中国大陆中东部地区以国家台网为主的100个分布均匀的宽频带地震台记录到的21个月的连续波形数据, 经过单台数据处理和互相关叠加计算后, 由时频分析法提取了研究区各台站对间瑞雷波的格林函数. 为了检验经验格林函数的可靠性和稳定性, 对沿部分路径的经验格林函数和频散曲线进行了质量评估. 检测结果表明, 自21个月叠加的台站对间背景噪声中提取的经验格林函数与实际的地震面波一致, 提取的格林函数可靠. 此外, 统计了使用从3—21个月不同长度数据叠加后, 经验格林函数信噪比大于10的频散曲线数目. 结果表明, 至少要使用12个月的数据才能提取到信噪比足够大, 数目足够多, 可用于反演面波速度结构的经验格林函数; 12个月的叠加时长, 可以保证30 s以下周期的频散曲线在时间上稳定.      

基于背景噪声的晋冀鲁豫交界地区瑞利波群速度层析成像

晋冀鲁豫交界地区位于华北克拉通的东南部, 研究该区地壳结构有助于进一步了解河北南部及其邻区地质构造演化特征、 发震机制、 发震规律。 本文收集了晋冀鲁豫交界地区两年84个地震台连续波形观测资料, 通过台站对波形互相关, 得到各台站对间的经验格林函数, 再利用图像分析技术提取瑞利波频散曲线。 本文基于面波直接成像方法得到不同深度剪切波速度结构。 结果显示, 该区部分速度异常与该区地表地质构造单元存在相关性, 速度异常呈区域分布特征, 山区和盆地呈现非常显著的速度差异, 浅层上地壳结构速度变化特征与沉积层厚度、 结晶基底埋深等区域构造密切相关, 中下地壳及上地幔存在不规则的速度异常, 与该区的地震发生存在相关性。     

Crustal S-wave velocity structure of the Yellowstone region using a seismic ambient noise method

The Yellowstone volcano is one of the largest active volcanoes in the world, and its potential hazards demand detailed seismological and geodetic studies. Previous studies with travel time tomography and receiver functions have revealed a low-velocity layer in the crust beneath the Yellowstone volcano, suggesting the presence of a magma chamber at depth. We use ambient seismic noise from regional seismic stations to retrieve short-period surface waves and then study the shallow shear velocity structure of the Yellowstone region by surface wave dispersion analysis. We first obtained a crustal model of the area outside of the Yellowstone volcano and then constructed an absolute shear wave velocity structure in combination with receiver function results for the crust beneath the Yellowstone volcano. The velocity model shows a low-velocity layer with shear velocity at around 1.3 km/s, suggesting that a large-scale magma chamber exists at shallow levels within the crust of the Yellowstone volcanic region.     

Estimated Green's function extracted from superconducting gravimeters and seismometers

Estimated Green's function (EGF) between stations has been extracted from ambient seismic noise, direct surface wave and coda waves. It is also confirmed by laboratory experiments on ultrasonics and theoretical derivations assuming diffusive wave field, equi-partition of modes or random sources on an enclosed surface. This method provides a new approach to study the crust and mantle structure at regional scale, continental scale and global scale. Following the achievements with seismometer records, the records of infrasonic station, hydrophone and microphone were also used to obtain the EGFs of different wave fields. Since superconducting gravimeter is a better long period instrument than regular seismometer, EGF at longer period is expected to be obtained with the cross correlation of gravity data. In this paper, we will show the EGFs extracted by cross-correlations between the superconducting gravimeters and the seismometers. Both the travel times and dispersion curves obtained from different data types are consistent. The result shows that it is possible to retrieve the deep structure by the cross correlation of gravity data.     

分布式光纤声波传感系统在近地表成像中的应用Ⅱ:背景噪声成像

随着城市化的发展,城市地球物理日益成为地球物理研究的重要方向,地震成像是构建城市地下空间三维/四维图像的重要手段,但面临观测成本高的困难.近年来国际上新发展的分布式光纤声波传感器作为高密度地震观测系统已经在地震层析成像方面得到了应用,在提高成像分辨率的同时,又降低了观测成本.本研究使用国产分布式光纤声波传感器开展了观测实验,利用480m埋地光缆记录了13h背景噪声,计算得到噪声互相关函数,获得了高频Rayleigh面波信号.采用多道面波分析方法提取相速度频散曲线,其结果与传统检波器记录和主动源结果较为一致.采用遗传算法反演得到了研究区内二维S波速度剖面,获得了下方沉积物横向变化特征.通过本次实验,初步验证了国产设备开展地震背景噪声成像研究、构建地下浅层结构模型的可行性.     

利用噪声研究地下介质速度结构及其变化在中国大陆的应用

背景噪声层析成像技术已被广泛应用于地壳和上地幔速度结构的研究,这种方法不依靠地震的发生和人工源爆破,只需记录连续的噪声信号而无需产生信号,因为噪声穿过地下介质时会携带信息,然后通过利用台站记录到的连续背景噪声数据进行互相关计算和叠加,即可得到台站间的经验格林函数,从而获取对地下结构的认识。这种方法已经很好地应用于中国的东北地区、华北克拉通、青藏高原以及华南地区,并成功地揭示了这些地区地壳与上地幔顶部的速度结构。此外近年来,一些学者还利用噪声互相关技术研究地下介质地震波速度随时间的变化,通过对汶川大地震前后连续噪声记录的研究发现,大震发生后呈现同震波速降低和震后波速逐渐恢复的特点,这表明可以通过观测地震波特性的变化来监测地下应力的变化,从而为大震的预测预防工作提供科学依据。本文主要综述了近些年来背景噪声技术及其在中国大陆地区的应用。     

Probabilistic inversion of Rayleigh-wave dispersion data: an application to Mt. Etna, Italy

We present a methodology for determining the elastic properties of the shallow crust from inversion of surface wave dispersion characteristics through a fully nonlinear procedure. Using volcanic tremor data recorded by a small-aperture seismic array on Mount Etna, we measured the surface waves dispersion curves with the multiple signal classification technique. The large number of measurements allows the determination of an a priori probability density function without the need of making any assumption about the uncertainties on the observations. Using this information, we successively conducted the inversion of phase velocities using a probabilistic approach. Using a wave-number integration method, we calculated the predicted dispersion function for thousands of 1-D models through a systematic grid search investigation of shear-wave velocities in individual layers. We joined this set of theoretical dispersion curves to the experimental probability density function (PDF), thus obtaining the desired structural model in terms of an a posteriori PDF of model parameters. This process allowed the representation of the objective function, showing the non-uniqueness of the solutions and providing a quantitative view of the uncertainties associated with the estimation of each parameter. We then compared the solution with the surface wave group velocities derived from diffuse noise Green’s functions calculated at pairs of widely spaced (~5–10 km) stations. In their gross features, results from the two different approaches are comparable, and are in turn consistent with the models presented in several earlier studies.     

利用背景噪声研究福建金钟库区地壳介质波速变化

基于2011年10月—2015年12月福建测震台网金钟库区周边台站的连续波形数据,利用地震背景噪声互相关技术提取台站对间的经验格林函数,分周期(1~2 s,2~4 s和4~8 s)测量每天经验格林函数与参考格林函数的走时偏移,进而计算了库区各台站对之间的相对地震波速变化。结果表明:2~4 s周期内地下介质相对波速变化表现出夏高冬低的年变规律,与库区水位变化及当地降雨量存在明显的正相关性,说明在此周期范围内,水位变化及降雨量对波速变化的影响最大。2013年8—10月,仙游地区接连发生了4次M_L≥4.0地震,其中在2013年9月4日仙游震群最大M_L5.0地震发生后,波速变化出现趋势性下降,幅度约为0.1%,随后波速逐渐恢复,在整个震群活动后期,对库区介质波速的影响不明显。     

甘肃华亭地区矿震记录特征分析

利用"十五"数字网络项目兰州宽频带地震台网资料,对甘肃华亭地区矿震记录波形采取多种数据处理技术,剔除噪声干扰、滤波后,分析解释了甘肃华亭矿震的典型记录特征。结果表明:华亭矿震的宽频带记录表现为近台记录大周期面波发育,持续时间较长;平凉台(PLT)场地放大作用明显,所计算的震级偏大。     

Rayleigh wave phase velocity tomography and strong earthquake activity on the southeastern front of the Tibetan Plateau

To investigate the relationship between velocity structure and earthquake activity on the southeastern front of the Tibetan Plateau, we make use of continuous observations of seismic ambient noise data obtained at 55 broadband stations from the regional Yunnan Seismic Network. These data are used to compute Rayleigh wave Green's Functions by cross-correlating between two stations, extracting phase velocity dispersion curves, and finally inverting to image Rayleigh wave phase velocity with periods between 5 and 34 s by ambient noise tomography. The results show significant lateral variations in crustal and uppermost mantle structures in the studied region. Phase velocity anomalies at short periods(5–12 s) are closely related to regional tectonic features such as sediment thickness and the depth of the crystalline basement. The Sichuan-Yunnan rhombic block, enclosed by the Honghe, Xiaojiang and Jianchuan faults, emerges as a large range of low-velocity anomalies at periods of 16–26 s, that inverts to high-velocity anomalies at periods of 30–34 s. The phase velocity variation in the vicinity of the Sichuan-Yunnan rhombic block suggests that the low-velocity anomaly area in the middle-lower crust may correspond to lower crustal channelized flow of the Tibetan Plateau. The spatial distribution of strong earthquakes since 1970 reveals that the Yunnan region is inhomogeneous and shows prominent characteristics of block motion. However, earthquakes mostly occur in the upper crust, with the exception of the middle-Yunnan block where earthquakes occur at the interface zone between high and low velocity as well as in the low-velocity zones, with magnitudes being generally less than 7. There are few earthquakes of magnitude 5 at the depths of 15–30 km, where gather earthquakes of magnitude 7 or higher ones which mainly occur in the interface zone between high and low velocities with others extending to the high-velocity abnormal zone.     

秦岭及周边地区瑞雷波方位各向异性

除了使用前人提取的面波频散曲线外, 还从秦岭及周边地区布设的59个宽频流动地震台的数据和中国地震局及各省局台网的数据中筛选出的地震事件波形和台站间噪声互相关格林函数中提取瑞雷波群速度频散曲线. 利用二维面波层析成像反演获得了瑞雷波周期为10—50 s的各向同性群速度及方位各向异性分布. 结果显示: 周期为10 s的各向同性群速度和方位各向异性分布与各构造单元存在明显的对应关系; 周期为10—50 s的面波在四川盆地和鄂尔多斯盆地内的快波方向多为近NS向. 与前人研究结果不同的是, 本文得到的秦岭、 大巴山构造带周期为10—50 s的面波快波方向均与山脉走向近似平行, 且与SKS波分析得到的快波方向一致. 这表明秦岭和大巴山之下整个岩石圈的快波方向都与山脉走向平行, 预示着秦岭和大巴山整个地壳, 甚至岩石圈发生了类似的变形. 由于四川盆地和鄂尔多斯盆地面波快波方向与SKS波结果差别较大, 推测青藏高原隆升扩展对这两个盆地的地壳基本无影响, 但对其岩石圈上地幔却产生了重大影响.      

Velocity model of the Hronov-Poříčí Fault Zone from Rayleigh wave dispersion

The Hronov-Po?í?í Fault Zone (HPFZ) is an active tectonic area with regularly occurring shallow earthquakes up to magnitude 5. For their exact locations, at least an average velocity model of the area is needed. A method of measuring local phase velocities of surface waves using the array of stations deployed permanently in the HPFZ is introduced. Seven regional and teleseismic events are selected to represent different backazimuths of propagation. Applicable range of periods is estimated for each event. The coherency of the waves reaching the array is constraining the short period range. The dimension of the array is a limiting factor for the long-periods. A dispersion curve of Rayleigh wave phase velocity measured at the vertical component and characterizing 1D properties of the target area is determined using the seven measurements for the interval from 1 to 40 s. An isometric method is used to invert the determined dispersion curve for shear and longitudinal velocity distribution from the surface to the depth of 65 km.     

鄂西地区宽频地震台阵背景噪声特征

利用鄂西地区长时间段宽频地震台站的三分量背景噪声记录,采用波形互相关方法得到台站对间的互相关函数,并通过聚束分析获得瑞雷波和勒夫波的慢度谱,研究鄂西地区背景噪声源的时空分布特征。结果表明,5~10 s周期范围,背景噪声来源于南太平洋且没有季节变化;10~20 s周期范围,慢度谱上显示明显的能量环,表明噪声源来源于多个方向,且表现出强烈和急剧的季节变化;20~40 s周期范围,慢度谱上也存在明显的能量环,其产生机制可能与此周期下提出的次重力波机制相似。在不同的周期范围内,噪声源分布方位有所不同,但在周期10~40 s范围噪声源在各方向均有分布。因此,利用长时间段连续噪声数据计算的互相关函数在周期10~40 s范围内满足背景噪声面波层析成像的理论前提。     

利用面波和接收函数联合反演滇西地区壳幔速度结构

考虑到面波频散对介质S波速度、接收函数对界面深度的各自敏感性优势,综合利用面波和接收函数资料实现联合反演,求取滇西地区壳幔速度结构. 本文利用适配滤波频时分析技术处理覆盖滇西地区的长周期面波资料,获得105～1050s周期范围内的面波群速度频散,进而利用分格反演方法提取研究区内1°×1°网格纯路径频散;基于滇西地区宽频带三分量远震记录,经反褶积后得到台站下方的远震P波接收函数. 联立面波纯路径频散信息和接收函数资料建立系统方程,利用阻尼最小二乘法实现联合反演,从而获得滇西地区壳幔S波速度结构. 结果表明,滇西地区以红河断裂为界,东西两侧壳幔结构存在明显差异,断裂西侧约20km深度处存在一厚度为10km左右的低速层,而东侧并不明显;滇缅泰块体上的畹町、沧源一带属于上地幔低速区,而另一个地幔低速区则位于滇中块体上的康滇古隆起上,两处地幔低速区与大地高热流分布、强震活动具有较好的对应关系.     

Structure of the upper crust in the Middle Asia by dispersion of surface waves from near earthquakes

Summary The records of Rayleigh surface waves from near earthquakes recorded at the seismic stations of Middle Asia have been analysed. The structure of the upper crust from dispersion of phase and group velocities of Rayleigh waves has been obtained.Communication presented to the XIII General Assembly of the European Seismological Commission, Brashov, 1972.     

<Emphasis Type="Italic">S</Emphasis>-wave velocity model for several regions of the Kamchatka Peninsula from the cross correlations of ambient seismic noise

The data from the seismic networks of the Kamchatka Branch of the Geophysical Survey of the Russian Academy of Sciences are used for calculating the cross correlations of seismic noise for the stationary digital stations over 2013 and for radio telemetric stations (RTS) in the region of the Klyuchevskoy volcano over the period from January 1, 2009 to May 31, 2013. Four hundred and two correlations overall are calculated. The fundamental-mode group velocities of the Rayleigh waves are calculated in the periods ranging from 5 to 50 s. The calculations for the region of the Klyuchevskaya group of volcanoes are based on the RTS data and cover the periods from 2 to 8 s. The two-dimensional (2D) maps of group velocity distributions in different periods are constructed with the use of the algorithm of surface wave tomography (Barmin, 2001). The velocity sections for the selected Kamchatka regions are reconstructed by the dispersion curve inversion technique (Mordret, 2014). For each region, the structure of the Earth’s crust and upper mantle down to a depth of 50 km was obtained.