Love wave contribution to the ambient vibration H/V amplitude peak observed with array measurements

This study applies array methods to measure the relative proportions of Love and Rayleigh waves in the ambient vibration wavefield. Information on these properties is of special relevance for frequencies around the horizontal-to-vertical (H/V) spectral amplitude ratio peak. The analysis of H/V curves, a popular technique in site characterisation, commonly assumes that the curves represent the frequency-dependent Rayleigh wave ellipticity. For the detailed interpretation of amplitudes or the inversion of the curves, it is therefore necessary to estimate and correct for the contribution of other wave types to the ambient vibration wavefield. I use available ambient vibration array measurements to determine the relative amount of Love and Rayleigh waves on the horizontal components by frequency-dependent analysis of the main propagation and polarisation directions, with a special emphasis on the H/V peak frequency as determined from the same recordings. Tests with synthetic data demonstrate the feasibility of this approach, at least in the presence of dominant source regions. Analysis of the data from 12 measurements at nine European sites, which include shallow as well as deep locations that span a wide range of impedance contrasts at the sediment-bedrock interface, indicates that the relative contribution of Rayleigh waves varies widely with frequency, from close to 0% to more than 70%. While most data sets show relative Rayleigh wave contributions between 40% and 50% around the H/V peak, there are also examples where Love waves clearly dominate the wavefield at the H/V peak, even for a site with a low impedance contrast. Longer-term measurements at one site indicate temporal variations in the relative Rayleigh wave content between day- and nighttime. Results calculated with the method introduced herein generally compare well with results of modified spatial autocorrelation analysis. These two methods might be used in a complimentary fashion, as both rely on different properties of the ambient vibration wavefield. This study illustrates that it is possible to measure the relative Rayleigh wave content of the noise wavefield from array data. Furthermore, the examples presented herein indicate it is important to estimate this property, as the assumption that there are an equal proportion of Love and Rayleigh waves is not always correct.     

Experiments on Exploration of Shallow Fine Structures and the Construction of the 1-D Velocity Model in the Pingtan Island, Fujian

112 short-period seismographs were set up in the 400km² area of Pingtan Island and its surrounding areas in Fujian. The combined observations of the airgun source and ambient noise source were carried out using a dense array to receive the 387 airgun signals excited around the island and one month of continuous ambient noise recording. The 1-D P-wave and S-wave shallow velocity model of Pingtan Island is obtained by the inversion of the airgun body wave''s first arrival time data, and the reliability of the velocity model is verified by using the surface wave phase velocity dispersion curve, which can provide initial model for subsequent 3-D imaging. The experimental results show that this experiment is a successful demonstration of local scale green non-destructive detection, which can provide basic data for shallow surface structure research and strong vibration simulation of the Pingtan Island.     

Factors Governing the Impact of Emerged Salt Diapirs on Water Resources

Salt diapirs in southern Iran are typically in contact with karstic and alluvial aquifers and consequently they are the most likely sources of groundwater salinization in this arid region. However, there are some salt diapirs that have no significant degradation effect on adjacent aquifers. Assessments of 62 of 122 Iranian‐emerged salt diapirs based on geological, geomorphological, hydrogeological, and hydrochemical investigations indicated that 45% of the studied salt diapirs did not have a negative impact on surrounding water resources, whereas 55% of the salt diapirs have degraded water quality of adjacent aquifers. The impacts ranged from low‐ to high‐grade salinization. We characterize here four major factors that control the impact of salt diapirs on surrounding water resources: (1) the evolutionary stage of the diapir, (2) the geological and (3) hydrogeological setting of the diapir, and (4) human activities. Identification of the major factors governing the influence of salt diapirs on the adjacent aquifers is necessary to understand the mechanism of salt diapir impact on adjacent aquifers, and subsequently to decide how to mitigate the deteriorating effect of the diapirs on the surrounding water resources.     

1D velocity structure of the Po River plain (Northern Italy) assessed by combining strong motion and ambient noise data

Strong ground motions recorded on the sedimentary deposits of the Po River alluvial plain during the Emilia (Northern Italy) Mw 5.7 earthquake of May 29, 2012 are used to assess the vertical profile of shear-wave velocity above the limestone basement. Data were collected by a linear array installed for site effect studies after the Mw 5.9 mainshock of May 20, 2012. The array stations, equipped with both strong and weak motion sensors, are aligned in the South–North direction, at distances ranging from 1 to 26 km from the epicenter. The vertical components of ground motion show very distinctive, large-amplitude, low-frequency dispersive wave trains. Wavelet analysis yields group-velocity dispersion curve in the 0.2–0.7 Hz frequency band. The availability of a long ambient noise record allows estimates of the site resonance frequency along with its stability among stations. The joint inversion of dispersion of surface waves and ellipticity curves derived from ambient noise H/V allows extending investigations down to the sediment-limestone interface, at a depth of about 5,000 m. Our results add new information about the velocity structure at a scale that is intermediate between the local scale already investigated by other authors with small-aperture arrays using ambient noise and the regional scale inferred from modeling of seismogram waveforms recorded at hundreds of kilometers from the source.     

CALCULATIONS ON THE THERMAL CONDITIONS AROUND A SALT DIAPIR*

Temperature modeling around the Vejrum salt structure in Denmark shows that the heat flow near the top is twice the regional heat flow. The temperature in the vicinity of the top of the diapir is thus anomalously high. Overlying permeable formations could therefore be attractive for geothermal energy exploitation. A negative temperature anomaly of ?20°C is calculated around the root of the salt body. Comparisons between calculated and measured temperatures are used to test the model. It is shown that measurements of heat flow (or temperature gradients) in wells penetrating the region of groundwater circulation could be used in identification of salt diapirs.     

Control of preexisting faults and near-surface diapirs on geometry and kinematics of fold-and-thrust belts (Internal Prebetic,Eastern Betic Cordillera)

We have determined, for the first time, the 3D geometry of a sector of the eastern Internal Prebetic comprised between Parcent and Altea diapirs, combining structural, borehole and multichannel seismic reflection data. The tectonic structure of the Jurassic-Cretaceous carbonate series is characterized by regional ENE-WSW fold-and-thrusts that interact with oblique N-S and WNW-ESE folds, detached over Triassic evaporites and clays. The structural style comprises box-shape anticlines, and N-vergent anticlines with vertical to overturned limbs frequently bordered by reverse and strike-slip faults. The anticlines surround a triangular broad synclinal structure, the Tárbena basin, filled by a late Oligocene to Tortonian sedimentary sequence that recorded folding and thrusting history. The location and geometrical characteristics of fold-and-thrusts may be controlled by the positive inversion of pre-existing Mesozoic normal faults, and by the position and shape of near-surface diapirs composed of Triassic rocks. Therefore, we propose an initial near-surface diapir emplacement of Triassic evaporitic rocks driven by late Jurassic to early Cretaceous rifting of the southern Iberian paleomargin. Thrusting and folding started during the latest Oligocene (∼28–23 Ma) roughly orthogonal to the NW-directed shortening. Deformation migrated to the south during Aquitanian (∼23–20 Ma), when tectonic inversion implied the left-lateral transpressive reactivation of N-S striking former normal faults and right-lateral/reverse reactivation of inherited WNW-ESE faults. We show two mechanisms driving the extrusion of the diapirs during contraction: lateral migration of a pre-existing near-surface diapir associated with dextral transpression; and squeezing of a previous near-surface diapir at the front of an anticline. Our study underlines the value of 3D geological modeling to characterize geometry and kinematics of complex fold-and-thrust belts influenced by preexisting faults and near-surface diapirs.     

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.     

Inversion of local S-wave velocity structures from average H/V ratios, and their use for the estimation of site-effects

H/V spectral ratios from microtremors areused to retrieve the S-velocity structurefrom a single ambient vibration record, byusing its relation to the ellipticity ofthe fundamental mode Rayleigh wave and theamplitude of observed H/V ratio.Constraints are needed in order to restrictthe possible range of solutions, and theinversion is applied to sites where thethickness of the unconsolidated sedimentsis approximately known from boreholeinformation. Within the uncertainty, theinverted structures agree well with theresults from other S-wave measuringtechniques such as downhole and cross-holemeasurements, and the analysis of ambientvibrations measured on an array.The influence of the inversion uncertaintyon site-amplification estimates forearthquakes is then investigated. For allinverted models, site response is computedfor a large number of events, which allowsto define the uncertainty by the aprioriunknown source position and mechanism of afuture earthquake. In most cases thevariability between the results obtainedfor the different models is much smallerthan the variability introduced by theunknown source position. The accuracy withwhich S-wave velocity structures can beretrieved from observed H/V ratios istherefore sufficient for an application ofthe method in seismic hazard analysis for aspecific site.     

Empirical site response evaluations: case studies in Israel

A closed-loop, time-efficient technique for site effect evaluation for seismic microzonation in urban areas is developed and implemented. It includes the following successive steps: microtremor measurements using a triangular array of three-component seismometers; estimation of the Nakamura horizontal-vertical spectral ratio, calculation of Rayleigh wave dispersion using the Aki spatial correlation method; inversion of the dispersion curve to determine a shear wave velocity model of the shallow subsurface; analytical computation of Rayleigh wave spectral ellipticity and 1-D SH transfer function for this model; and a comparison of empirical and analytical ratios and correction of the model if required. The technique has been tested at several locations in a number of towns near the active Jordan Rift valley. A reasonable agreement is obtained when comparing empirical spectral ratios and analytical transfer functions.     

基于微动方法研究五大连池火山区尾山火山锥浅层剪切波速度结构

为了研究五大连池火山区尾山火山锥浅层三维波速结构特征,在尾山火山锥附近区域布设了无线地震检波器密集台阵,记录连续地震背景噪声波形数据.基于微动方法(拓展空间自相关方法)提取了台站间2~5 Hz频率范围的Rayleigh面波相速度频散曲线.利用面波层析成像方法反演获得2~5 Hz Rayleigh面波二维相速度图像,基于每一个网格节点的频散曲线,进一步反演获得了尾山火山锥附近区域地表至700m深度的三维剪切波速度结构.成像结果显示:在0~150m较浅深度,靠近尾山火山锥区域显示为相对高速异常,远离火山锥区域则显示为相对低速异常.而至150~700m较深深度,波速异常特征与浅部相反,靠近尾山火山锥的区域显示为相对低速异常,远离火山锥的区域显示为相对高速异常.在远离尾山火山锥区域,浅层的相对低速异常可能与松散沉积层有关,深部的高速异常则反映了结晶变质岩的影响.在靠近尾山火山锥区域,浅部的相对高速异常应该反映了出露地表的玄武岩,而深部的相对低速异常则可能反映了火山通道周围广泛发育的破碎裂隙结构及其火山喷发后孔隙流体填充的影响.     

Diapirism of depleted peridotite — a model for the origin of hot spots

A model is proposed for the origin of hot spots that depends on the existence of major-element heterogeneities in the mantle. Generation of basaltic crust at spreading centers produces a layer of residual peridotite ～20–25 km thick directly beneath the crust which is depleted in Fe/Mg, TiO₂, CaO, Al₂O₃, Na₂O and K₂O, and which has a slightly lower density than undepleted peridotite beneath it. Upon recycling of this depleted peridotite back into the deep mantle at subduction zones, it becomes gravitationally unstable, and tends to rise as diapirs through undepleted peridotite. For a density contrast of 0.05 g cm^?3, a diapir 60 km in diameter would rise at roughly 8 cm y^?1, and could transport enough heat to the base of the lithosphere to cause melting and volcanism at the surface. Hot spots are thus viewed as a passive consequence of mantle convection and fractionation at spreading centers rather than a plate-driving force.It is suggested that depleted diapirs exist with varying amounts of depletion, diameters, upward velocities and source volumes. Such variations could explain the occurrence of hot spots with widely varying lifetimes and rates of lava production. For highly depleted diapirs with very low Fe/Mg, the diapir would act as a heat source and the asthenosphere and lower lithosphere drifting across the diapir would serve as the source region of magmas erupted at the surface. For mildly depleted diapirs with Fe/Mg only slightly less than in normal undepleted mantle, the diapir could provide not only the source of heat but also most or all of the source material for the erupted magmas. The model is consistent with isotopic data that require two separate and ancient source regions for mid-ocean ridge and oceanic island basalts. The source for mid-ocean ridge basalts is considered to be material upwelling at spreading centers from the deep mantle. This material forms the oceanic lithosphere. Oceanic island basalts are considered to be derived from varying mixtures of sublithospheric and lower lithospheric material and the rising diapir itself.     

基于贝叶斯理论的接收函数与环境噪声联合反演

基于Bayes反演理论(Tarantola,1987,2005),在接收函数非线性复谱比反演方法基础上(刘启元等,1996),本文讨论了接收函数与地震环境噪声Rayleigh波相速度频散的联合反演.本文采用修正后的快速广义反射/透射系数方法(Pei et al., 2008,2009) 计算Rayleigh波相速度频散, 并引入地壳泊松比的全局性搜索.数值检验表明:(1)接收函数与环境噪声的联合反演能够有效地解决反演结果对初始模型依赖的问题,即使对地壳速度结构仅有非常粗略的初始估计(例如,垂向均匀模型),本文方法仍能给出模型参数的可靠估计;(2)由于环境噪声与接收函数在频带上的适配性明显优于地震面波,接收函数与环境噪声的非线性联合反演能更好地约束台站下方近地表的速度结构;对于周期范围为2~40s的环境噪声相速度频散,利用本文方法能够可靠推测台站下方0~80 km深度范围的S波速度结构, 其浅表速度结构的分辨率可达到1 km; (3)本文方法能够可靠地估计地壳泊松比,泊松比的全局性搜索有助于合理解释接收函数和环境噪声的面波频散数据.利用本文方法对川西台阵KWC05台站观测的接收函数与环境噪声的联合反演表明,该台站下方地壳厚度为44 km,上地壳具有明显的高速结构,24~42 km范围的中下地壳具有低速结构.该台站下方地壳的平均泊松比为0.262,壳内低速带的泊松比为0.27.     

冷底辟与热底辟地震信息

冷底辟是指盐、膏、泥岩等塑性地层向上侵入或刺穿所产生的构造样式;热底辟主要指岩浆侵入或喷发所产生的底辟构造样式。不论哪种底辟构造,在人工地震剖面上反映均非常明显,其底辟内部呈杂乱反射,而围岩反射界面在此突然中断。在岩浆底辟处常出现磁力高异常,并有较高的层速度,冷底辟具有较低的层速度,盐膏底辟也具有较高的层速度。构造与火山地震关系密切,而冷底辟也可以形成破坏性地震,是今后值得注意的研究课题     

Improved salt imaging in a basin context by high resolution potential field data: Nordkapp Basin,Barents Sea

The seismic imaging of salt diapirs in the Nordkapp Basin gave rise to considerable problems in defining their shape and volume. Independent information was added by integrating the interpretation with high resolution gravity and magnetic data. We developed a novel, iterative workflow, separated into sub‐categories: sediments, salt structures, basement and Moho. Distinctions between the sources of the anomalies from different depths was achieved by utilizing the different decay characteristics of gravity, gravity gradiometry and high resolution magnetic anomalies. The workflow was applied to the southern part of the Nordkapp Basin. It started with the sedimentary model derived from seismics, populated with measured densities and magnetic susceptibilities and a starting model for the base salt. The residual after the removal of this model was interpreted in terms of a crustal model, including flexural isostatic calculations for the Moho with the sedimentary load. The residual after the removal of crustal and early sedimentary model was used to tune the salt model. As these major and minor modelling steps depend on each other, an iterative process was applied to stepwise improve the density and magnetic susceptibility model. The first vertical gradient of gravity and the magnetic field were found to give most information about the cap rock of the diapirs. The improvement in salt imaging, integrated with results from controlled‐source electromagnetic and magneto‐telluric modelling is shown for the salt diapir Uranus, where a well, terminated in the salt, constrains the minimum of the depth to base salt.     

天然源面波勘探台阵对比试验

为了对比天然源面波勘探不同台阵布局的探测效果, 筛选出探测成果可靠、 效率高和便于野外施工的天然源面波勘探台阵阵形, 在天水市黄土覆盖区的同一场地分别用4种常见的阵形进行数据采集试验, 并对各种阵形数据使用空间自相关法或扩展空间自相关法提取相应的频散曲线, 通过反演得到了试验点地下的浅层速度结构模型. 分析对比试验结果表明: 4种台阵提取的频散曲线数值很相近; 频散谱能量集中度较高的是嵌套式等边三角形和圆形台阵, L形和直线形台阵相对分散; L形台阵低频段(4—8 Hz)比直线形台阵差, 其高频段(8—40 Hz)比直线形台阵好. 针对直线形台阵在高频段信噪比较低的情况, 在确保探测成果可靠性的前提下, 为了提高探测效率, 提出了在同一直线形台阵开展天然源与人工源面波联合勘探的数据采集方法. 实验结果证实, 这种联合勘探方法不仅可弥补直线形台阵高频段的不足, 确保探测精度和结果的可靠性, 而且还能实现“高低”频兼顾, 即“深浅”兼顾.      

基于环境噪声测试的隐伏断层场地H/V谱比特征分析

基于四川石棉县安顺场隐伏断层穿越区的60余次环境噪声测试,探讨了隐伏断层对环境噪声H/V谱比特征的影响,并且通过研究区所获卓越频率对其覆盖层厚度进行估算,以覆盖层厚度的变化揭示了隐伏断层的行迹。结果显示:① 在排除可能受河流影响引发的高频区后,隐伏断层带处的H/V谱比曲线多呈现高频、低放大系数和多峰值频率等特点;② 覆盖层厚度在部分强风化区和河流高频区存在5—10 m的差异;③ 环境噪声阵列反演所得部分覆盖层厚度变化较大区域与隐伏断层的形迹吻合。      

Mechanism of development of inversion folding in the subsalt

The numerical modeling of the Archimedean upwelling of a salt bed predicts that, with mature diapirs (from the finger stage and later), an about 2-km-thick zone of inversion (mirror) folding, where the suprasalt diapirs correspond to the subsalt synclines while the suprasalt inter-diapir sags correspond to the anticlines in the subsalt, should be expected immediately below the bottom of the salt deposits. These deformations decay with depth. The development of the inversion folding is exceptionally due to the flow caused by the rising of the unstable layer.     

基于邻域算法的瑞利面波垂直-水平振幅比及频散曲线联合反演及应用

瑞利面波垂直-水平振幅比(或ZH振幅比)是一个随频率变化的函数,对于台站下方浅层地壳结构非常敏感,且具有和频散资料不同的深度敏感核,是传统频散反演方法的一个很好的补充,从而可以将基阶瑞利面波的ZH振幅比和面波频散数据联合起来更好地反演获得观测台站下方的速度结构.本文提出了基于邻域算法的面波频散曲线与ZH振幅比联合反演方法,我们进行了基于理论模型的模拟测试,证明了联合反演是一种更为可靠的反演方法,且能更好地约束浅层地壳结构.相比于频散曲线单独反演,联合反演不仅可以精确反演获得地壳的Vs结构,对分层地壳的Vp/Vs也能很好地约束.然后我们将联合反演算法应用于实际测量数据,获得了中国西南昆明台(KMI)下方更为准确的地壳横波速度结构及Vp/Vs模型.     

地脉动空间自相关方法反演 浅层S波速度结构

在近地表地球物理领域, 基于地脉动(或称背景噪声)提取的面波频散曲线反演地下S波速度结构是一种简单经济的工程勘察方法. 本文基于地脉动的空间自相关方法对一个微型台阵观测的背景噪声记录进行处理, 介绍了一种简单易行的提取频散曲线的数据处理方法, 获得了6.7—23 Hz频段的可靠频散曲线; 通过对该观测频散曲线与预测模型的频散曲线进行拟合, 反演得到S波速度结构. 结果表明, 该速度结构与钻孔直接测试的结果相吻合.      

On the repeatability and consistency of three-component ambient vibration array measurements

Ambient vibration measurements with small, temporary arrays that produce estimates of surface wave dispersion have become increasingly popular as a low-cost, non-invasive tool for site characterisation. An important requirement for these measurements to be meaningful, however, is the temporal consistency and repeatability of the resulting dispersion and spatial autocorrelation curve estimates. Data acquired within several European research projects (NERIES task JRA4, SESAME, and other multinational experiments) offer the chance to investigate the variability of the derived data products. The dataset analysed here consists of repeated array measurements, with several years of time elapsed between them. The measurements were conducted by different groups in different seasons, using different instrumentations and array layouts, at six sites in Greece and Italy. Ambient vibration amplitude spectra and locations of dominant sources vary between the two measurements at each location. Still, analysis indicates that this does not influence the derived dispersion information, which is stable in time and neither influenced by the instrumentation nor the analyst. The frequency range over which the dispersion curves and spatial autocorrelation curves can be reliably estimated depends on the array dimensions (minimum and maximum aperture) used in the specific deployment, though, and may accordingly vary between the repeated experiments. The relative contribution of Rayleigh and Love waves to the wavefield can likewise change between repeated measurements. The observed relative contribution of Rayleigh waves is generally at or below 50%, with especially low values for the rural sites. Besides, the visibility of higher modes depends on the noise wavefield conditions. The similarity of the dispersion and autocorrelation curves measured at each site indicates that the curves are stable, mainly determined by the sub-surface structure, and can thus be used to derive velocity information with depth. Differences between velocity models for the same site derived from independently determined dispersion and autocorrelation curves—as observed in other studies—are consequently not adequately explained by uncertainties in the measurement part.