3D S-wave velocity structure of the Ningdu basin in Jiangxi province inferred from ambient noise tomography with dense array

The Ningdu basin, located in southern Jiangxi province of southwest China, is one of the Mesozoic basin groups which has exploration prospects for geothermal energy. A study on the detailed velocity structure of the Ningdu basin can provide important information for geothermal resource exploration. In this study, we deployed a dense seismic array in the Ningdu basin to investigate the 3D velocity structure and discuss implications for geothermal exploration and geological evolution. Based on the dense seismic array including 35 short-period (5 s-100 ?Hz) seismometers with an average interstation distance of ～5 ?km, Rayleigh surface wave dispersion curves were extracted from the continuous ambient noise data for surface wave tomographic inversion. Group velocity tomography was conducted and the 3D S-wave velocity structure was inverted by the neighborhood algorithm. The results revealed obvious low-velocity anomalies in the center of the basin, consistent with the low-velocity Cretaceous sedimentary rocks. The basement and basin-controlling fault can also be depicted by the S-wave velocity anomalies. The obvious seismic interface is about 2 ?km depth in the basin center and decreases to 700 ?m depth near the basin boundary, suggesting spatial thickness variations of the Cretaceous sediment. The fault features of the S-wave velocity profile coincide with the geological cognition of the western boundary basin-controlling fault, which may provide possible upwelling channels for geothermal fluid. This study suggests that seismic tomography with a dense array is an effective method and can play an important role in the detailed investigations of sedimentary basins.     

重庆浅层S波速度结构研究

基于重庆台网及邻区固定台站自2011年1月至2018年10月观测到的远震波形记录,利用接收函数直达P波振幅约束浅层结构的方法计算了重庆及周边地区台站下方浅层S波速度结构,结果表明：重庆浅层S波速度结构与盆山构造明显相关,盆地内表现为低速异常,与沉积层特征相对应,大巴山、大娄山区域则表现为相对高速的特征;华蓥山断裂S波速度高于断裂两侧的沉积层速度,表现出四川盆地沉积层中间薄两边厚的特点;华蓥山以东的川东地区多为套滑脱层构造,其浅层S波速度表现为低速异常。最后讨论了2010年以来重庆5次显著地震的孕震环境：川东滑脱构造地区的垫江M_S4.4地震和石柱M_S4.5地震与该区明显的低速异常有关;荣昌M_S4.7和M_S4.8地震的震源区无明显的高低速特征,可能与注水相关;武隆M_S5.0地震发生在高低速交界部位的有利于积累应变的高速体一侧。     

Estimation of site effects using strong motion data of BYTNet array in Turkey

Simultaneous estimation of effects of source, propagation path, and local site amplification was carried out using observed strong motion records in a frequency range from 0.8 to 20 Hz for the purpose of empirical evaluation of the local site effects in different geological conditions in the northwestern part of Turkey. The analyzed data are S-wave portions of 162 accelerograms from 39 shallow events observed at 14 sites of BYTNet array. A spectral separation method was applied to the observed S-wave spectra. The solutions for source spectra, inelasticity factor of propagation path for S-waves (Q _s-value), and factor of site amplification at each site were obtained in a least squares sense. In the analysis, we assumed that the factor of the site amplification at a reference site is the same as that of theoretical amplification of S-waves to the soil model whose bottom layer has an S-wave velocity around 2.15 km/s. The estimated Q _s-value of the propagation path is modeled as Q _s(f)?=?87.4f^0.78. The estimated site amplifications are characterized into three groups. The sites in the first group belong to rock site with no dominant peaks at a frequency range of 2 to 10 Hz. The second group of hard soil sites is characterized with moderately dominant peaks at a frequency of 5 Hz. The last group for soft soil sites has common peaks at a frequency of 4 Hz with larger amplitudes than those in the hard soil group. We, then, compare the amplifications with average S-wave velocity in top 30 m of the shallow S-wave profiles and proposed linear empirical formula between them at each frequency. We, furthermore, inverted the observed amplification factors into S-wave velocity and Q _s-value profiles of the deep soil over the basement.     

Broadband ground-motion simulation of the 24 May 2014 Gokceada (North Aegean Sea) earthquake (Mw 6.9) in NW Turkey considering local soil effects

On 24 May 2014, a Mw 6.9 earthquake occurred in the west of Gokceada Island, northern Aegean Sea. The earthquake was close to Canakkale, Enez, Tekirdag cities, and damaged 300 buildings in the Marmara Region, NW Turkey. We simulated its broadband (0.1–10 Hz) ground motions including 1D deep and shallow structures soil amplification effects at the 12 strong ground motion stations in the western Marmara Region. The 1D deep velocity structures from the focal layer to the engineering bedrock with an S-wave velocity of 0.78 km/s in different azimuthal directions were tuned by comparing the observed group-velocity dispersion curves of Rayleigh and Love waves from the mainshock with theoretical ones. We also added the shallow parts from previous surveys into the 1D models. Synthetic seismograms on the engineering bedrock were generated using the discrete wave number method with a source model and the 1D deep velocity structures. Then the surface motion was generated considering shallow soil amplification. The synthetic seismograms are generally in good agreement with the observed low and high-frequency parts at most of the stations indicating an appropriateness of the source model and the 1D structural model.     

松辽盆地沉积层结构的短周期地震背景噪声成像研究

使用位于松辽盆地内部的NECESSArray台阵连续两年背景噪声数据,通过波形互相关和多重滤波方法提取到2~14 s较短周期的Rayleigh波群速度和相速度频散曲线,基于快速行进（FMM）面波成像方法得到群速度和相速度成像结果,并采用最小二乘迭代线性方法反演获得了松辽盆地深至12 km的三维S波速度结构.本文成像结果显示：松辽盆地内部S波速度分布的横向不均匀性与该区域的构造单元呈现出良好的空间对应关系.从地表至下方的6 km深度,盆地北部比南部表现出更加强烈的低速异常,这一特征可能与盆地南北的沉积构造差异有关.中央坳陷区低速异常的边界与嫩江断裂走向相互平行,表明盆地基底断裂对盆地形成演化具有一定的控制作用.在垂直速度结构剖面中,2.9 km·s^-1的S波速度等值线与地震反射剖面显示的盆地基底深度大致对应.基于S波速度模型和盆地基底速度（2.9 km·s^-1）,我们获得精细的松辽盆地沉积层厚度模型,结果表明松辽盆地的沉积层厚度分布呈现出中间厚、四周薄的特征,中央坳陷区的沉积层厚度范围大约在3~6 km.     

S-wave envelope broadening characteristics of microearthquakes in the Canary Islands

This study analyzes the S-wave envelope broadening characteristics of 290 earthquakes recorded by 14 stations of the Spanish National Seismograph Network in the Canary Islands region. The S-wave peak delay time (t _p ) and envelope duration (t _q ) parameters are evaluated phenomenologically to infer the strength of velocity inhomogeneities of the medium along each seismic ray path. Crustal (0?≤?h?≤?18 km) and upper mantle (18?<?h?≤?80 km) events are analyzed separately. Results in the frequency range 1 to 12 Hz for hypocentral distances from 30 to 600 km show that both t _p and t _q increase according to a power of hypocentral distance and they are independent of frequency. The spatial distribution of the peak delay time reveals weak strength of heterogeneity in most of the region at shallow depths. Relatively strong inhomogeneous zones are generated under the island of Tenerife and Gran Canaria at depths of 11–22 km.     

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.     

S-wave velocity structures of sediments estimated from array microtremor records and site responses in the near-fault region of the 1999 Chi-Chi,Taiwan earthquake

The 1999 Chi-Chi, Taiwan earthquake, M_W = 7.6, caused severe damage in the near-fault region of the earthquake. In order to evaluate site effects in the near-field strong motions we estimate S-wave velocity structures of sediments at four sites using array records of microtremors. We also recalculated S-wave velocity structures at other four sites previously reported. To show the validity of the estimated S-wave velocity structures we separate empirical site responses from aftershock records using the generalized inversion method and show the agreement between empirical and theoretical site responses. We also show an observed fact that suggests soil nonlinearity during the Chi-Chi earthquake by comparing horizontal-to-vertical spectral ratios (HVRs) for main shock records with HVRs for aftershock records. Then we calculate one-dimensional equivalent-linear site responses using the estimated S-wave velocity structures and the main shock records observed on the surface. It is found that site amplification due to thick (about 6 km) sediments is one of the important factors for explaining the long-period velocity pulses of about 5 to 10 sec observed at sites in the footwall during the Chi-Chi earthquake. It is also found that the theoretical site responses of shallow soft sediments at sites that sustained severe damage in the hanging wall shows significant amplification around 1 sec. As the amplitude of velocity pulses with period around 1 sec is most critical in causing damage to ordinary buildings of moderate heights, our results suggest that the 1-sec period velocity pulses, amplified by the site response of shallow sediments should contribute to the severe damage during the Chi-Chi earthquake.     

羌塘盆地中部地区地壳S波速度结构及构造意义

羌塘盆地是我国最大的海相盆地.本文根据在羌塘盆地内布设的27个宽频带地震观测台站记录的远震波形数据,利用非线性复谱比反演算法得到各台站下方100 km深度范围内S波速度结构.结果表明.羌塘地区Moho深度较为平缓,平均深度为61 km;北羌塘地壳内低速层广泛分布;北羌塘具有两个较大的沉积盆地,龙尾错和白滩湖坳陷,沉积厚度分别有10 km和15 km.尽管北羌塘下地壳受到强烈的新生代火山岩作用改造,但是这种深部岩浆热作用会加速烃源岩中有机质的热演化历程.北羌塘两个盆地具有很好的油气前景.与北羌塘低速层分布相比,南羌塘下低速层更深,可能与班公怒江洋于中生代的俯冲消减及拉萨地体北向俯冲有关.     

基于地震背景噪声的四川威远地区浅层速度结构成像研究

地震S波速度是防震减灾中场地分类、强地面震动模拟建模等的重要参数.地震背景噪声成像方法可以重建地下浅层三维结构以及探测浅层速度结构变化,为页岩气开采提供参考,有效规避地震灾害风险.本文收集了2015年11月至2016年2月间四川威远地区50个流动台站记录到的垂直分量连续波形记录,利用波形互相关提取2~6 s的基阶Rayleigh波经验格林函数,采用一维地壳结构模型拟合相速度频散曲线,获得了该区域近地表5 km以内的三维S波速度模型.结果表明,威远地区深度2 km以内的S波速度为2.0~2.7 km·s^-1;2.0~5.0 km深度的S波速度横向分布不均匀,西北地区大于2.7 km·s^-1,东南地区在2.3~2.8 km·s^-1之间,分别与背斜构造和较厚的盆地沉积层相对应.3个剖面图均表明S波速度随着深度的增加而逐渐增大,并且在台站S29附近的页岩气田钻井底部与2016年1月7日四川威远M_L3.9地震震源位置较为接近,均处于松散的沉积盖层与较为坚硬的花岗岩基底的分界处,推测它们之间可能存在一定的联系.     

利用接收函数方法确定青藏高原东北缘近地表S波速度

公式推导和数值实验都表明远震接收函数直达波幅度对近地表S波速度比较敏感,可以用于确定近地表（如沉积层）的S波速度值.本文利用该方法分析了青藏高原东北缘区域数字地震台网36个台站的远震资料,得到了研究区近地表S波速度分布特征.研究结果表明：在深度5±2.5 km范围内,祁连山—河西走廊区域近地表S波速度均值较小;位于地形梯度带的甘东南区域近地表S波速度呈现较为复杂的分布.通过与该区域沉积层厚度结果对比,发现近地表S波速度与沉积层厚度呈现较好的负相关.本研究是对目前接收函数方法的有效补充,在探测研究近地表速度结构方面,相比人工地震方法其成本几乎可忽略且环保.在目前地震波观测数据海量增长的背景下,该方法对于确定不同区域近地表结构具有较大潜力.     

S-wave velocity structure in Tangshan earthquake region and its adjacent areas from joint inversion of receiver functions and surface wave dispersion*

Using the seismic records of 83 temporary and 17 permanent broadband seismic stations deployed in Tangshan earthquake region and its adjacent areas (39°N–41.5°N, 115.5°E–119.5°E), we conducted a nonlinear joint inversion of receiver functions and surface wave dispersion. We obtained some detailed information about the Tangshan earthquake region and its adjacent areas, including sedimentary thickness, Moho depth, and crustal and upper mantle S-wave velocity. Meanwhile, we also obtained the v_P/v_S structure along two sections across the Tangshan region. The results show that: (1) the Moho depth ranges from 30 km to 38 km, and it becomes shallower from Yanshan uplift area to North China basin; (2) the thickness of sedimentary layer ranges from 0 km to 3 km, and it thickens from Yanshan uplift region to North China basin; (3) the S-wave velocity structure shows that the velocity distribution of the upper crust has obvious correlation with the surface geological structure, while the velocity characteristics of the middle and lower crust are opposite to that of the upper crust. Compared with the upper crust, the heterogeneity of the middle and lower crust is more obvious; (4) the discontinuity of Moho on the two sides of Tangshan fault suggests that Tangshan fault cut the whole crust, and the low v_S and high v_P/v_S beneath the Tangshan earthquake region may reflect the invasion of mantle thermal material through Tangshan fault.     

Site response in Tecoman,Colima, Mexico—II. Determination of subsoil structure and comparison with observations

In a companion paper local transfer functions were estimated at Tecoman using earthquake and microtremor data. In this paper, the subsoil structure at this city is investigated using seismic refraction and cross-correlation of noise records as a case study. P- and S-wave refraction profiles were measured at five sites within the city. Standard analysis constrained only very shallow layers. The P-wave refraction deployment was also used to record ambient vibration. These data were processed using an extension of the SPAC (SPatial AutoCorrelation (Aki, 1957) [1]) method; cross-correlation is computed between station pairs and the results are inverted to obtain a phase velocity dispersion curve. Penetration depth was larger than that from the refraction experiments but the shear-wave velocity of the basement could not be determined. For this reason, additional microtremor measurements were made using broad band seismometers with a larger spacing between stations. The results allowed to constrain the shear-wave velocity of the basement. Site amplification computed for the final profiles compare well with observed ground motion amplification at Tecoman. The case of Tecoman illustrates that even a simple subsoil structure may require crossing data from different experiments to correctly constrain site effects.     

Geophysical Images of the North Anatolian Fault Zone in the Erzincan Basin, Eastern Turkey, and their Tectonic Implications

The collision between the Arabian and Eurasian plates in eastern Turkey causes the Anatolian block to move westward. The North Anatolian Fault (NAF) is a major strike-slip fault that forms the northern boundary of the Anatolian block, and the Erzincan Basin is the largest sedimentary basin on the NAF. In the last century, two large earthquakes have ruptured the NAF within the Erzincan Basin and caused major damage (M _s = 8.0 in 1939 and M _s = 6.8 in 1992). The seismic hazard in Erzincan from future earthquakes on the NAF is significant because the unconsolidated sedimentary basin can amplify the ground motion during an earthquake. The amount of amplification depends on the thickness and geometry of the basin. Geophysical constraints can be used to image basin depth and predict the amount of seismic amplification. In this study, the basin geometry and fault zone structure were investigated using broadband magnetotelluric (MT) data collected on two profiles crossing the Erzincan Basin. A total of 24 broadband MT stations were acquired with 1–2 km spacing in 2005. Inversion of the MT data with 1D, 2D and 3D algorithms showed that the maximum thickness of the unconsolidated sediments is ~3 km in the Erzincan Basin. The MT resistivity models show that the northern flanks of the basin have a steeper dip than the southern flanks, and the basin deepens towards the east where it has a depth of 3.5 km. The MT models also show that the structure of the NAF may vary from east to west along the Erzincan Basin.     

基于地震探测的盐源盆地构造特征及动力学意义

基于覆盖盐源盆地的短周期天然地震台阵和布设的一条人工地震测线所获得的地震数据,从中提取地震能量属性,并通过地震层析成像获得该地区的浅部地震速度结构,继而对短周期地震台阵一个月的噪声数据进行互相关得到经验格林函数,再通过时频分析获得相速度频散曲线,反演获得不同深度的S波速度分布.研究结果显示,盐源盆地地震特征主要分为三层...     

利用SPAC法估算地壳S波速度结构

S波速度结构能够反映地球介质的物性差异,是地壳内低速区结构特征判别的重要依据.本文尝试利用空间自相关法(SPAC法)从地震台站微动信号的垂直分量中提取瑞利波相速度频散曲线,通过对频散曲线的反演获得地下介质的S波速度结构.以国家数字测震台网8个宽频带地震台站的实测微动数据为例,采用SPAC方法获得了首都圈地区北京附近约30 km 深度范围内的一维S波速度结构.结果表明,该区结晶基底埋深较浅约2 km;分别在5~8 km 和12~16 km 深处发育S波低速层;8 km 和 20 km 处是S波速度差异较大的速度分界面.这一结果与以往地震学及人工地震探测结果较为吻合,表明SPAC法估算地壳S波速度结构是可行、有效的.     

华北地区地壳上地幔S波三维速度结构

利用华北地区大型流动地震台阵的记录资料,采用近震和远震联合成像方法,得到了水平分辨率0.5°×0.5°、深至600km的S波速度结构.研究结果表明,上地壳S波速度结构与地表地质构造基本一致,燕山—太行山山脉均呈现高速异常,延庆—怀来盆地、大同盆地表现为低速异常,华北盆地内部的拗陷和隆起分别呈现低速和高速.唐山地区中地壳、山西裂陷盆地中下地壳存在明显的低速异常,可能分别与流体和热物质作用有关,有利于形成孕育强震的地质构造环境.90km的速度结构图像依然与地表的构造特征有较大的相关性,可能说明深部结构对地表构造有一定的控制作用.燕山隆起区岩石圈的厚度可达120～150km左右,华北盆地的岩石圈厚度可能在80km左右,太行山地区的岩石圈厚度介于两者之间.山西裂陷盆地上地幔低速层较厚,反映了该区不稳定的构造环境造成了地幔热物质的上涌.华北盆地下方220～320km出现的高速异常体,可能揭示了华北盆地上地幔仍然存在拆沉后残留的难熔、高密度的古老岩石圈地幔.研究区东部地幔转换带呈低速异常,推测可能与太平洋板块俯冲至该区下方地幔转换带前缘120°E左右的俯冲板块相变脱水有关.     

基于广义反演法的不同位置地震引起盆地放大效应的研究——以仙台盆地为例

基于强震观测记录,采用广义反演法和HVSR法分析日本仙台盆地48个强震台站以及场地v_S30和盆地深度对放大效应的影响。研究表明,相比于广义反演法,HVSR法能够较好地给出场地主频,但会明显低估放大效应的幅值;处于盆地外不同位置的浅源地震引起的盆地放大效应差异明显,仙台盆地南部海域地震引起的放大效应最大,盆地北部陆地地震引起的放大效应最小;盆地南部和东部地震引起的S波的放大效应与v_S30的相关性较强,北部地震的放大效应与v_S30基本不相关;盆地S波的放大效应与盆地深度在0.5～5 Hz频段内的相关性较强,在0.25～0.5 Hz和5～10 Hz频段内基本不相关。     

Velocity Structure of the West-Bengal Sedimentary Basin, India along the Palashi-Kandi Profile Using a Travel-time Inversion of Wide-angle Seismic data and Gravity Modeling – An Update

The 2-D shallow velocity structure along the north-south Palashi-Kandi profile in the West Bengal sedimentary basin has been updated by travel-time inversion of seismic refraction, wide-angle reflection and gravity data. A six-layer shallow model up to a depth of about 7 km has been derived. The first layer, which has an average velocity of 2.0 kms^?1, represents the alluvium deposit, which rests over the shale formation with average velocity of 3.0 kms^?1. The thin (200 m) Sylhet limestone, observed at a nearby Palashi well, remains hidden in the present data set. Hence a 200-m thin layer with a velocity of 3.7 kms^?1, corresponding to the Sylhet limestone, has been assumed to be present throughout the profile. The fourth layer with a velocity of 4.5–4.7 kms^?1 at a depth of 1.7–2.4 km represents the Rajmahal traps. The ‘skip’ phenomenon and rapid amplitude decay of first arrivals indicate a low-velocity zone (LVZ) in the study area. Using the ‘skip’ phenomena and wide-angle reflection data, identified on seismograms, the LVZ with a velocity of 4.0 kms^?1, indicating the Gondwana sediments, has been delineated below the Rajmahal traps. The next layer with a velocity 5.4–5.6 kms^?1 overlying the crystalline basement (5.8–6.25 kms^?1) may be associated with the Singhbhum group of meta volcanic rock that has been exposed in the western part of the basin. The basement lies at a variable depth of 4.9 to 6.8 km. The overall uncertainties of various velocity and boundary nodes are ± 0.15 kms^?1 and ± 0.5 km, respectively. The elevated basement feature in the north might have acted as a structural barrier for the deposition of Sylhet limestone during the Eocene epoch. The seismically derived shallow structure correctly explains the observed Bouguer gravity anomaly along the profile. The addition of reflections in the present analysis provides a stronger control on the depths and velocities of basement and overlying sedimentary formations, compared to the earlier model derived mainly by the first arrival seismic data.     

汶川地震余震震源参数及地震动衰减与场地影响参数反演分析

基于强震动观测台钻孔数据, 统计得到了四川、 甘肃地区20 m深的浅硬土层场地平均剪切波速模型． 通过与美国加州地区的波速模型对比, 结合四川、 甘肃地区地震预报的地壳模型, 延拓剪切波速模型至40 km． 应用四分之一波长法计算了浅硬土层场地的平均场地放大系数． 利用场地放大系数, 消除13次余震中浅硬场地的台站场地响应, 反演了龙门山断层上、 下盘的介质品质因子以及13次余震的有效应力降． 通过与普通土层场地的强震动记录对比, 对浅硬土层场地的放大系数进行调整, 得到了深厚土层场地的平均放大系数．