Site response and source spectra of S waves in the Zagros region, Iran

S wave amplitude spectra from shallow earthquakes with magnitudes ranging between 4.2 and 6.2 in the Zagros region of Iran that occurred between 1998 and 2008 are used to examine source parameters and site response of S waves. A generalized inversion scheme has been used to separate the source, propagation path, and local site effects from S wave spectra. For removing the trade-off between source and site terms and propagation effects (including geometric and anelastic attenuation), the spectral amplitudes of the records used were corrected for attenuation and geometrical spreading function using a path model proposed by Zafarani and Soghrat (Bull Seism Soc Am 102:2031–2045, 2012) for the region. We assume a Brune’s point source model to retrieve source parameters like corner frequency, moment magnitude, and high-frequency fall off coefficient, for each event. When the source spectra are interpreted in terms of Brune’s model, the average stress drops obtained are about 7.1 and 5.9 MPa (71 and 59 bars), respectively for the eastern and western Zagros regions. Stress drops range from 1.4 to 35.0 MPa (14 to 350 bars), with no clear dependence on magnitude. The results in terms of stress drop and S wave seismic energy indicate that the Zagros events are more similar to interplate earthquakes of western North America than to intraplate events of eastern North America. The method also provides us with site responses for all 40 stations individually and is an interesting alternative to other methods, such as the H/V method. A new empirical relationship between body-wave magnitudes and moment magnitude has been proposed for the Iranian plateau using derived seismic moment from the inversion.     

Application of Regional Phase Amplitude Tomography to Seismic Verification

v--vWe have applied tomographic techniques to amplitude data to quantify regional phase path effects for use in source discrimination studies. Tomography complements interpolation (kriging) methods by extending our ability to apply path corrections into regions devoid of calibration events, as well as raising levels of confidence in the corrections because of their more physical basis. Our tomography technique solves for resolvable combinations of attenuation, source-generation, site and spreading terms. First difference regularization is used to remove singularities and reduce noise effects.¶In initial tests the technique was applied to a data set of 1488, 1.0 Hz, P_g/L_g amplitude ratios from 13 stations for paths inside a 30° by 40° box covering western China and surrounding regions. Tomography reduced variance 60%, relative to the power-law distance correction traditionally applied to amplitude ratios. Relative P_g/L_g attenuation varied with geologic region, with low values in Tibet, intermediate values in basins and high values for platforms and older crust. Spatial patterns were consistent with previous path effect studies in Asia, especially local earthquake coda-Q. Relative spreading was consistent with expected values for P_g and L_g. Relative site terms were similar to one another, yet some tradeoff with attenuation was evident.¶Tomography residuals followed systematic trends with distance, which may result from the evolution from direct to coda phases, focusing, model tradeoff or data windowing effects. Examination of the residuals using a kriging interpolator showed coherent geographical variations, indicating unmodeled path effects. The residual patterns often follow geological boundaries, which could result from attenuating zones or minor blockages that are too thin to be resolved, or that have anisotropic effect on regional phases. These results will guide efforts to reparameterize tomography models to more effectively represent regional wave attenuation and blockage. The interpolated residuals also can be combined with predictions of the tomographic model to account for path effects in discrimination studies on a station by station basis.     

S wave inelastic attenuation and site effect in the northern Tianshan area

In this paper,according to the relationship of observation spectrum with source spectrum,travel-path attenuationfactor and site effect,we use the digital seismic data of moderate-small earthquakes to invert the soft rock site ef-fect and S wave inelastic attenuation under the constraint of site effect,and determine the parameters of path andsite.The parameters obtained by this method exclude the topographic effect of seismic stations,and truly reflectthe seismic effect of free soft rock sites and the S wave inelastic attenuation in the studied region.The average softrock site effect is about 1.5 times in the frequency domain of 2～4 Hz and the S wave quality factor is obtained asQS=278 f0.346,which can be directly applied to the stochastic modeling of ground motions in the studied region.     

北天山地区S波非弹性衰减和场地效应研究

采用中小地震的数字观测资料,依据地震观测谱和震源谱、传播路径衰减因子及场地效应的关系式,在对场地效应约束的条件下,反演S波非弹性衰减和软基岩场地效应,确定研究区路径、 场地参数. 该方法确定的参数, 避免了台站地形效应的影响,真实地反映了研究区自由软基岩场地的地震效应和传播路径非弹性衰减的影响. 研究区的平均软基岩场地效应, 在2~4 Hz频率范围为1.5倍左右; S波品质因子为QS=278 f 0.346.这些参数可以直接用于研究区的地面运动随机模拟.      

QLg tomography in Gujarat,Western India

We propose a novel Lg attenuation tomography model (Q_Lg tomography) for the state of Gujarat, Western India, using earthquake data recorded by the Gujarat Seismic Network, operated by the Institute of Seismological Research in Gandhinagar. The waveform dataset consist of 400 3-component recordings, produced by 60 earthquakes with magnitude (M_L) spanning from 3.6 to 5.1, recorded at 60 seismic stations having epicentral distances spanning between 200 and 500 km. Spectral amplitude decays for Lg wave displacement were obtained by generalized inversion at 17 frequencies spanning between 0.9 and 9 Hz. Lg wave propagation efficiency was measured by Lg/Pn spectral ratio categorizing as efficient ratio ≥6 for 86%, intermediate ratio of 3–6 for 10% and inefficient ratio <3 for 4% paths of total 400 ray paths. The earthquake size and quality of waveform recorded at dense network found sufficient to resolve lateral variation of Q_Lg in Gujarat.Average power-law attenuation relationship obtained for Gujarat as Q_Lg(f) = 234f^0.64, which corresponds to high attenuation in comparison to peninsular India shield region and other several regions around the world. Q_Lg tomography resolves the highly attenuating crust of extremely fractured Saurashtra region and tectonically active Kachchh region. The Gujarat average attenuation is also lying in between them. The low attenuation in Cambay and Narmada rift basins and extremely low attenuation in patch of Surendranagar area is identified. This study is the first attempt and can be utilized as pivotal criteria for scenario hazard assessment, as maximum hazard has been reported in highly attenuating tectonically active Kachchh region and in low attenuating Cambay, Narmada and Surendranagar regions. The site and source terms are also obtained along with the Q_Lg inversion. The estimated site responses are comparable with observed local geological condition and agree with the previously reported site amplifications at the same sites. The source terms are comparable with local magnitude estimated from Network. The M_w (Lg) is nearly equivalent to M_L (GSN) and the slight differences are noted for larger magnitude events.     

Crustal attenuation characteristics of S-waves beneath the Eastern Tohoku region,Japan

An inversion method was applied to crustal earthquakes dataset to find S-wave attenuation characteristics beneath the Eastern Tohoku region of Japan.Accelerograms from 85 shallow crustal earthquakes up to25 km depth and magnitude range between 3.5 and 5.5were analyzed to estimate the seismic quality factor Q_s. A homogeneous attenuation model Q_s for the wave propagation path was evaluated from spectral amplitudes, at 24 different frequencies between 0.5 and 20 Hz by using generalized inversion technique. To do this, non-parametric attenuation functions were calculated to observe spectral amplitude decay with hypocentral distance. Then, these functions were parameterized to estimate Q_s. It was found that in Eastern Tohoku region, the Q_s frequency dependence can be approximated with the function 33 f~(1.22) within a frequency range between 0.5 and 20 Hz. However,the frequency dependence of Q_s in the frequency range between 0.5 and 6 Hz is best approximated by Q_s(f) = 36f~(0.94) showing relatively weaker frequency dependence as compared to the relation Q_s(f) = 6 f~(2.09) for the frequency range between 6 and 15 Hz. These results could be used to estimate source and site parameters for seismic hazard assessment in the region.     

Determination of source,path and site parameters based on non-linear inversion of accelerometric data in Greece

Source, propagation path and site conditions are the factors affecting seismic ground motion. Consequently, recordings acquired at a seismic station are formed by the convolution of these three factors. In this work S-wave acceleration Fourier spectra of earthquakes recorded at local and regional scale, by the ITSAK accelerometric network for the period 2010–2016, are modeled as a product of source, propagation path (including geometric and anelastic attenuation) and site effects. The data set consists of 136 crustal earthquakes occurred in the broader Aegean area, with magnitudes 4.2?≤?M_w?≤?6.5 and epicentral distances 20 km?≤?R?≤?350 km, recorded at 112 broadband accelerometric stations installed at sites with various geologic conditions. Based on this data set, an iterative Gauss–Newton inversion method to solve the non-linear problem and retrieve the different terms of source, propagation path and site, is applied. This method uses an initial input model trying to find the best and at the same time a stable solution for the inverted parameters, which are, moment magnitude (M_w), corner frequency (f_c), attenuation quality factor (Q_s?=?Q_of^α), slope of the geometric attenuation (1/R^γ) and site transfer function (S(f)). The initial values of the starting model can be either known from other studies or inferred within a reasonable range. Depending on the level of knowledge on these input parameters, the associated standard deviation can be adjusted (large values for unknown parameters or small values for parameters which are well constrained). Results of the analyses exhibit satisfactory agreement of estimated source parameters with those proposed by seismological centers in Greece and propagation path properties similar to the ones determined in relevant previous studies for the same region. In addition, the site transfer functions obtained by the non-linear inversion are comparable with those calculated for the same sites using either standard spectral ratio or horizontal-to-vertical spectral ration (HVSR—receiver function) techniques. The aforementioned results are encouraging for reliable earthquake source parameters, propagation path properties and site effect assessment, in areas of intermediate to high seismicity.     

Empirical analysis of path effects on prediction equations of pseudo‐velocity response spectra in northern Japan

We study path effects on prediction equations of pseudo‐velocity response spectra (natural period of 0.1–5.0 s) in northern Japan, where heterogeneous attenuation structure exists. The path effects have been examined by comparing the regression analysis results for two different prediction equations. The first equation consists of a single term of anelastic attenuation conventionally. The second equation consists of two terms of anelastic attenuation in consideration of the heterogeneous attenuation structure. In the second equation, we divide a source‐to‐site distance into two distances at the attenuation boundary beneath the volcanic front. The boundary is considered to separate the relatively high Q fore‐arc side mantle wedge (FAMW) from the low Q back‐arc side mantle wedge (BAMW). Strong motion records (hypocentral distances less than 300 km) from interplate and intraslab events with Mw 5.1–7.3 are used. Regression analysis results show that the standard errors are significantly reduced by the second prediction equation at short periods (0.1–0.5 s), whereas the difference in standard errors from both prediction equations is negligible at intermediate and long periods. The Qs values (quality factor for S‐wave) converted from two anelastic attenuation coefficients for the second prediction equation are remarkably similar to the path‐averaged Qs values for the FAMW and BAMW by other studies using spectral inversion method. From these findings, we conclude that the path effects on the prediction equation of pseudo‐velocity response spectra are satisfactorily accomplished by the second prediction equation. Copyright © 2009 John Wiley & Sons, Ltd.     

Estimation of local site effects by a generalized inversion scheme using observed records of ‘Small-Titan’

A simple-extended array of observation for strong motions, named ‘Small-Titan’, was installed in Sendai City, Japan in 1997. This array system consists of 20 observation sites, which have an average of about 4 km for the station-to-station distance. Each observation station was located at various kinds of soils to effectively obtain local site effects. Since its completion, the system has obtained about 160 motion records. As is well known, ground motions are influenced by the three main factors such as source, path and local site effects. An accurate prediction of ground motions should thus be performed based on the three factors. In this context, theoretical methods used commonly lack enough information related with source and path effects. A generalized inversion scheme to separate the effects of source, propagation path and local site from S-wave spectra was applied to the records of Small-Titan obtained in the Sendai region, Japan. The objective of this paper is to estimate local site effects in Sendai City and to obtain meaningful information of input motions available to earthquake-resistant design, aimed at setting in terms of local soil conditions. This paper also describes the applicability of the generalized inversion scheme to seismic zoning.

Site amplification factors estimated from the generalized inversion scheme vary significantly, reflecting the effects of soil characteristics at each site. Strong site-dependency was found in the site amplification factors. In particular, the local site amplification factors in the western side are smaller than those in the eastern side, across the Nagamachi–Rifu Fault that exists in the central part of the Sendai area. Results from the generalized inversion scheme are compared with those obtained by the direct spectral ratio scheme, the H/V method and the bootstrap inversion. There is a good agreement in the amplification factors between the generalized inversion scheme, the direct spectral ratio scheme and the bootstrap inversion scheme. The H/V method for estimating site amplification is found to extract the same dominant periods as the other methods. This paper especially emphasizes a good agreement in not only dominant period, but also amplification factor.

The generalized inversion scheme is a useful method in separating the effects of source, propagation path and local site and in revealing the site amplification. In addition, the generalized inversion scheme provides meaningful information of input motions available to earthquake-resistant design, aimed at setting in terms of local soil conditions. Moreover, this paper indicates that the generalized inversion scheme is applicable to the seismic zoning toward disaster mitigation.     

Separation of Source, Propagation and Site Effects from S Waves¶of Local Earthquakes in Bursa Region, Northwestern Turkey

—?We have used micro-earthquake recordings (M= 1.8–4.1) of local events in the distance range of 5–60?km in order to quantify the attenuation and site effects in the vicinity of the Bursa city, Marmara region, Turkey. The data set consists of 120 three-component recorded accelograms from 69 earthquakes, recorded at six stations. Each station is deployed on different geologic units, such as massive limestone, slope deposit and Quaternary young sediments, in the framework of the Marmara Poly-Project.¶In this study a nonparametric inversion method was applied to acceleration records from the Bursa region to estimate source, site and path effects using a two-step inversion. At the first step, we determined attenuation functions by analyzing the distance dependence of the spectral amplitudes and retrieved values of Q _s (f) = 46.59f ^0.67. At the second step, the corrected S-waves spectral records for the attenuation function, including the geometrical spreading effect, were inverted to separate source and site response for 21 different frequencies selected between 0.5 and ～25?Hz. The near-surface attenuation, κ value, was also estimated by using the model proposed by Anderson and Hough (1984) at each site. We observed that κ₀ is smaller for stations located on rock site (I?dιr, SIGD, κ₀～0.004) compared to the one that is located on Neogene sediment (Çukurca, SCKR, κ₀～0.018).¶Site amplifications from inversion showed that the station located within the Bursa basin, Çukurca (SCKR), is the most important site with about 4.0 amplification value at 1.8?Hz. Demirta? (SDEM) amplifies the spectral amplitudes about 3.0 times at 2.0?Hz, SHMK about 3.5 times between 2.5 and 3.5?Hz and SHMT nearly reaching 3.5 times between 1.5 and 4.0?Hz. However, stations located on the Uluda? Mountain Massif (SKAY and SIGD), which correspond to a deep limestone geological unit, have the smallest amplification, that values between 0.6 and 1.4.     

时间域相邻道地震波衰减成像研究

在时间域中,利用地震记录中的振幅与上升时间信息可以计算出岩石介质的品质因子Q值.由于震源能量难以测量,本文通过共炮点道集相邻道循环对比的方法,消除了震源对计算过程的影响,推导出时间域相邻道振幅与上升时间衰减成像的计算公式,并在振幅衰减成像的计算过程中,去除了非粘滞性衰减成分.为了更全面地综合利用波形信息,并充分利用各种方法的长处,将波速成像、振幅衰减成像、上升时间衰减成像集成到一个处理流程中,通过在岩体结构研究中的实际应用,证明这种综合方法具有灵活、实用、可靠的特点.     

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…     

考虑介质对地震波吸收衰减作用的场地放大效应估计

当参考场地震源距小于研究场地震源距时,传统谱比法低估了场地放大效应,且参考场地与研究场地震源距相差越大,低估程度越大;当参考场地震源距大于研究场地震源距时,传统谱比法高估了场地放大效应,且参考场地与研究场地震源距相差越大,高估程度越大,计算结果依赖参考场地的选择。上述问题存在的原因是传统谱比法估计场地放大效应时未考虑地震波传播过程中介质的吸收衰减作用,故传统谱比法适用于参考场地与研究场地地质构造类似且震源距相差较小的情况,为此本文通过考虑介质对地震波的吸收衰减作用,提出改进谱比法,可适用于地质构造类似、任意震源距的情况。本文分别采用未改进方法和改进方法估计汶川地震中渭河盆地多个台站场地反应,证实了传统谱比法存在的上述问题。采用改进谱比法估计的场地放大倍数均在2 Hz以下出现1个峰值,且这些峰值多为最大值;采用未改进谱比法估计的场地放大倍数峰值出现的频率多大于2 Hz;改进谱比法可信度较高。     

2000年姚安MS65地震震源及近邻区的低QS

利用震源距23 km范围内观测的2000年姚安MS65地震余震记录，计算了震源及近邻区域的尾波规一直达S波在频率15～20 Hz之间的傅里叶谱振幅.结果显示谱振幅随震源距增大而增大， 在对谱振幅进行了震源辐射方向性校正之后, 才出现谱振幅随震源距衰减的现象.由此获得了震源及近邻区域S波的Q(f)值,可表示为QS(f)=89f098其值比由尾波得出的姚安地区的平均QC(f)值低得多，表明了震源破裂带的强烈非均匀性对QS(f)的重大影响.     

浅壳与上地幔地震的水平加速度谱阻尼修正系数模型

对日本K-net和KiK-net台网中6 466条浅壳与上地幔地震动记录进行了统计分析,考虑四种场地类别分别建立了浅壳与上地幔地震的加速度谱阻尼修正系数模型,并利用随机效应模型将模型总残差分离为事件间残差、事件内残差（分为场地间残差和场地内残差）,计算其相应的标准差,探究地震震源、路径、场地等因素对模型误差的影响。研究结果显示:所提出的模型能较好地拟合实际计算值,其随机误差整体随谱周期和阻尼比的增大而增大;震源效应导致的随机误差在谱周期小于2.0 s时小于其它效应;在大部分谱周期上,场地效应导致的随机误差小于路径效应及其它效应;事件间残差与震源和震源深度显著相关,场地内残差与震源距显著相关。      

朝鲜2009年和2013年两次核爆的地震学特征对比研究

利用布设在长白山地区临时地震台站接收到的朝鲜核爆的波形资料,对2009年5月25日和2013年2月12日两次朝鲜核试验的地震学特征进行比较.震中距范围从145km到420km.采用P/S型谱比值方法识别朝鲜核爆,通过与2009年3月20日长春地震和2013年1月23日沈阳地震事件的比较,表明在频率大于3 Hz时P/S型谱比值能够有效识别发生在中朝边境地区的地下核试验.选定参考台站,利用区域震相Pg波的振幅谱比值计算朝鲜核爆至各台站路径上的相对衰减.结合介质速度模型,在一定程度上反映了长白山地区衰减情况,为进一步研究长白山地区衰减层析成像提供初始模型.     

浙江地区Lg波路径衰减关系及 台站场地响应参数

基于浙江省数字地震台网2007—2012年记录的台湾地区355次M_S≥4.0地震的波谱比数据, 使用场地响应与路径衰减联合反演的方法, 对穿越台湾海峡特定路径的Lg波衰减关系及浙江地区36个测震台站的场地响应参数进行了计算. 计算中选取频率为1—7 Hz、 间隔为0.2 Hz, 计算得到了三分向(U--D, N--S, E--W)的衰减关系及场地响应参数. 三分向衰减参数关系为γ_UD(f)=0.0055f0.94, γ_NS(f)=0.0064f1.11和γ_EW(f)=0.0048f0.90. 浙江地区36个台站均获得了稳定的场地响应, 其中庆元台在高频端表现出较明显的放大效应, 其它台站在所研究的频率范围内各分向场地响应大致相当, 并未显示出明显的方向性特征. 最后对场地响应的异常部分进行了具体分析和验证, 表明本文计算结果正确可靠.      

Use of spectral acceleration data for determination of three-dimensional attenuation structure in the Pithoragarh region of Kumaon Himalaya

Three-dimensional attenuation structures are related to the subsurface heterogeneities present in the earth crust. An algorithm for estimation of three-dimensional attenuation structure in the part of Garhwal Himalaya, India has been presented by Joshi (Curr Sci 90:581–585, 2006b; Nat Hazards 43:129–146, 2007). In continuation of our earlier approach, we have presented a method in which strong motion data have been used to estimate frequency-dependent three-dimensional attenuation structure of the region. The border district of Pithoragarh in the Higher Himalaya, India, lies in the central seismic gap region of Himalaya. This region falls in the seismic zones IV and V of the seismic zoning map of India. A dense network consisting of eight accelerographs has been installed in this region. This network has recorded several local events. An algorithm based on inversion of strong motion digital data is developed in this paper to estimate attenuation structure at different frequencies using the data recorded by this network. Twenty strong motion records observed at five stations have been used to estimate the site amplification factors using inversion algorithm defined in this paper. Site effects obtained from inversion has been compared with that obtained using Nakamura (1988) and Lermo et al. (Bull Seis Soc Am 83:1574–1594, 1993) approach. The obtained site amplification term has been used for correcting spectral acceleration data at different stations. The corrected spectral acceleration data have been used as an input to the developed algorithm to avoid effect of near-site soil amplification term. The attenuation structure is estimated by dividing the entire area in several three-dimensional block of different frequency-dependent shear wave quality factor Q _β (f). The input to this algorithm is the spectral acceleration of S phase of the corrected accelerogram. The outcome of the algorithm is given in terms of attenuation coefficient and source acceleration spectra. In the present study, this region has been divided into 25 rectangular blocks with thickness of 10 km and surface dimension of 12.5 × 12.1 km, respectively. Present study gives three-dimensional attenuation model of the region which can be used for both hazard estimation and simulation of strong ground motion.     

高斯线调频连续小波变换的时频能量衰减因子方法及其应用

本文基于高斯线调频连续小波变换，提出了能够反映震源区或近场区小地震波形在震源深度、震源尺度、震源破裂机制、地震波传播途径、地震波衰减等方面的差异特征信息的特征指标，即小波变换的时频能量衰减因子方法。为了验证该方法的有效性，选取了震源破裂机制有明显差异特性的天然地震、爆破或塌方资料以及非强地震孕震区同一地点发生的多次小震资料。研究结果表明，天然地震与爆破或塌方记录的连续小波变换时频能量衰减因子有明显区别；正常地震活动背景下地震记录的连续小波变换时频能量衰减因子变化比较稳定；利用该方法，通过对某一地区的连续观测，可望为强震预测提供有效判据。