Site response study in Abruzzo (Central Italy): underground array versus surface stations

In this site response study we examined local earthquakes recorded at surface stations of a local seismic network and at a temporary underground seismic array installed in a tunnel underneath the Gran Sasso Massif in Abruzzo (central Italy). This allowed us to compare the seismic site response beneath the mountain and on the surface in similar geological environment (soft rock sites). We applied spectral ratios method on different segments of the seismograms and used different reference spectra in the 1–20 Hz frequency band. We found little or no amplification effects at most of the surface stations whereas site transfer functions evaluated with respect to underground sites show an amplification factor up to 6 in the 1–8 Hz frequency range. Coda spectral ratios estimated at soft rock sites are confirmed as good estimates of shear wave transfer function.     

Broadband frequency-dependent amplification of seismic waves across Bucharest, Romania

The determination of seismic amplitude amplification is a fundamental contribution to seismic hazard assessment. While often only high-frequency amplitude variations (>1 Hz) are taken into account, we analyse broadband waveforms from 0.14 to 8.6 Hz using a temporary network of 32 stations in and around the earthquake-prone city of Bucharest. Spectral amplitudes are calculated with an adaptive multiple-taper approach. Across our network (aperture 25 km × 25 km), we find a systematic northwest/southeast-oriented structural influence on teleseismic P-wave amplitudes from 0.14 to 0.86 Hz that can be explained by constructive interference in the dipping Cenozoic sedimentary layers. For higher frequencies (1.4–8.75 Hz), more local site effects prevail and can be correlated partly among neighbouring stations. The transition between systematic and localised amplitude variations occurs at about 1 Hz.     

Evidence for fault-related directionality and localized site effects from strong motion recordings of the 2003 Boumerdes (Algeria) earthquake: Consequences on damage distribution and the Algerian seismic code

The Algiers–Boumerdes region has been struck by a destructive magnitude 6.8 (M_w) earthquake on May 21, 2003. The study presented in this paper is based on main shock strong motions from 13 stations of the Algerian accelerograph network. A maximum 0.58g peak ground acceleration (PGA) has been recorded at 20 km from the epicenter, only about 150 m away from a PGA of 0.34g, with both a central frequency around 5 Hz, explained by a strong very localized site effect, confirmed by receiver function technique results showing peaks at 5 Hz with amplitudes changing by a factor of 2. Soil amplifications are also evidenced at stations located in the quaternary Mitidja basin, explaining the higher PGA values recorded at these stations than at stations located on firm soil at similar distances from the epicenter. A fault-related directionality effect observed on the strong motion records and confirmed by the study of the seismic movement anisotropy, in agreement with the N65 fault plan direction, explains the SW–NE orientation of the main damage zone. In the near field, strong motions present a high-frequency content starting at 3 Hz with a central frequency around 8 Hz, while in the far field their central frequency is around 3 Hz, explaining the high level of damage in the 3- to 4-story buildings in the epicentral zone. The design spectra overestimate the recorded mean response spectra, and its high corner frequency is less than the recorded one, leading to a re-examination of the seismic design code that should definitively integrate site-related coefficient, to account for the up to now neglected site amplification, as well as a re-modeling of the actual design spectra. Finally, both the proposed Algerian attenuation law and the worldwide laws usually used in Algeria underestimate the recorded accelerations of the 6.8 (M_w) Boumerdes earthquake, clearly showing that it is not possible to extrapolate the proposed Algerian law to major earthquakes.     

祁连山主动源观测台站场地响应分析

场地响应是精确测定震源参数和稳健估计出高精度地壳介质衰减变化的关键影响因素之一,采用H/V尾波谱比法,对祁连山主动源观测台网中40个短周期观测台站的场地响应进行了分析.结果表明,祁连山主动源观测台网中40个短周期观测台站均存在不同程度的放大(衰减)作用,根据40个台站观测波形的尾波场地响应曲线特征,可将其分为平坦类、放...     

Characterization of regional and local scattering effects from small-aperture seismic array recordings

In 1998, four small-aperture arrays separated by 20 km have been deployed in the southern French Alps in order to record the natural seismicity during two consecutive months. One of the main objectives of this experiment was to characterize the heterogeneities that influence the wavefield propagation, by analysis of the coda characteristics recorded by each of the seismic arrays.The time-azimuth-velocity characteristics of the coda phases have been obtained using a high-resolution wavenumber decomposition method.A statistical analysis, using the coda characteristics of the whole data set (20 regional earthquakes) recorded by the four arrays has been performed and lead to the calculation of the density of scattered energy within the medium. Three regional heterogeneous areas (> 10 km) have been characterized, and are located in the N–NE, W and S directions from the four arrays. Scattered energy is also located at local distances (< 10 km) from the four arrays.The comparison of the waveforms recorded (i) with one of the array, for different groups of earthquakes, and (ii) with the four arrays, for one group of earthquakes, show that the geometrical configuration of the source–scattering area–recording site system can strongly influence the energetic distributions related to the secondary phases of the seismograms.In particular cases, we experimentally show that the interactions of the direct wavefield with the heterogeneous structures of the medium can be sufficiently energetic to induce secondary scattered phases that dominate the seismic motions recorded at a given site. In such case, these effects should be taken into account for the quantification of the expected ground motion recorded during an earthquake.     

Calculating the Station Corrections to Earthquake Energy Class and the Acoustic Impedance for Kamchatka Stations

This paper describes a method for calculating the station corrections to earthquake energy class and the acoustic impedance based on records of earthquake codas in the 0.8–1.8 Hz frequency band. The calculation is based on measurements of the relative coda level at two stations, with one of the two being used as the reference. The data were recorded by digital seismometers: CMG-3TB (GURALP), CMG-5TD (GURALP), CMG-6TD (GURALP), and GMS-AC-73i HHV (GeoSIG). We present estimates of class corrections for the Kamchatka network of seismic stations.     

Site amplification effect of coda waves in West Tunnan Experimental Site

The site effect of coda waves for four types of basement under 6 stations was studied by using the local earthquakes from digital seismograph network of West Yunnan Experimental Site. Results from the band 0.75 to 24 Hz showed that the site effect was abviously interrelated with not only the frequency but also the site condition of the station. For those stations which sited on the hard-rock, such as Jianchuan (limestone), Shuanglang (limestone), Qiaohou (sandstone) and Eryuan (metamorghic rock), the site amplification as a function of frequency was similar in shape and the changes in amplification were flat. A typical behavior of sediment was measured at Qina station. The site amplification remained high at low frequency, but turned to a low level as frequency increased. However, in contrast to Qina, the variations of site amplification at Yangbi station were quite different. Such different behavior between the two sites could be caused by the difference in sediment thickness under valley basins. It could be seen that the site amplification pattern of horizontal motion approached approximately to the vertical despite of some differences. The Fourier spectra of small earthquakes also showed the character of site. It can be expected that the site response has an important effect on the shape of the spectra and that the size of the magnitude is not likely responsible for the spectra recorded at the same site. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,15, Supp., 1993.     

Regional features of coda attenuation of aftershocks of Wuding, Yunnan earthquake in near field

IntroductionSeismic wave is an important information source to analyze and study seismic source and internal structure of the earth. Studies of wave on different frequency ranges can provide many features of seismic source and internal structure of the earth. In recent years in quantitative study of seismic wave seismologists at home and abroad have developed a wide study field on the coda. Especially seismologists have done a lot of work in coda geneses, mechanism and its attenuation to explo…     

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.     

云南武定地震余震近场尾波衰减的区域特征

根据地方震尾波散射理论,利用云南武定MS6.5主震后近场田心和云龙两个流动数字化台站记录的余震序列地震资料,以不同中心频率（1.5～20Ｈｚ）取样分频研究了近场尾波振幅衰减率（f ）的空间分布特征．结果表明,两个台站给出的尾波振幅衰减率（f ）存在明显的小区域范围非均匀性变化．位于主发震断层两侧的云龙台（东侧）和田心台（西侧）测得的尾波振幅衰减率（f ）值,当 f ＜6.0Ｈｚ时基本一致．相应的Qc（f ）值为66～120的低值;但在高频范围（f 6.0Ｈｚ）,由于穿过主断层的近场尾波遭到断层破碎带强烈的吸收,因而云龙台测得的尾波振幅衰减率（f ）值,比不穿过主断层的田心台给出的（f ）降低30％,并讨论了产生近场尾波振幅衰减空间分布差异的原因．本文的研究结果还表明,波源因子A0（f）不仅与震源强度相关,也与频率犳相关,用1.5～20Ｈｚ求出的Ｓ波平均自由程犔在主断层东侧为30ｋｍ、西侧为40ｋｍ ,由此证实了地震波穿过主断层破碎带时有更强的散射．     

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.     

Estimation of topographical effects at Narni ridge (Central Italy): comparisons between experimental results and numerical modelling

In the present work the seismic site response of Narni ridge (Central Italy) is evaluated by comparing experimental results and numerical simulations. The inhabited village of Narni is located in central Apennines at the top of a steep massive limestone ridge. From March to September 2009 the site was instrumented with 10 weak-motion stations, 3 of which located at the base of the ridge and 7 at the top. The velocimetric network recorded 642 events of M_L up to 5.3 and hypocentral distance up to about 100 km. The great amount of data are related to the April 2009 L’Aquila sequence. The site response was analyzed using both reference (standard spectral ratio, SSR) and non reference spectral techniques (horizontal to vertical spectral ratio, HVSR). Moreover directional analyses were performed in order to evaluate the influence of the ridge orientation with respect to the selected source-site paths. In general the experimental results show amplification factors for frequencies between 4 and 5 Hz for almost all stations installed along the crest. The SSR technique provides amplification factors up to 4.5 in a direction perpendicular to the main elongation of the ridge. The results obtained from the data analyses were used as a target for bidimensional and tridimensional numerical simulations, performed using a hybrid finite-boundary element method and a boundary element method for 2D and 3D modelling, respectively. In general, the results obtained through numerical simulation fit well the experimental data in terms of range of amplified frequencies, but they underestimate by a factor of about 2 the observed amplifications.     

基于台阵记录的土层山体场地效应分析

选取2008年5月25日至8月7日期间由甘肃省文县上城山地形效应台阵获取的12次汶川地震余震事件(M_S≥4.0),在分析其地震动基本参数的基础上,采用参考场地谱比(RSSR)法和水平-竖向谱比(HVSR)法,研究了不同地震作用下上城山地形台阵的场地效应.分析结果显示:随着高程和覆盖层的增加,记录台站地震动的PGA呈增大趋势,地震频谱形状由宽变窄;上城山台阵记录到的地震波在地形基阶频段(2—4 Hz)和浅部土层频段(7—9 Hz)的幅值明显放大,RSSR曲线显示山顶NS向的土层频段谱比大于山体地形频段谱比;由于土层山体竖向地震动在中高频段放大,使得HVSR方法谱比结果在中高频段较RSSR方法所得结果明显偏低,而在山体基阶频段附近两种方法的谱比值接近.松散土层山体的台阵记录特征体现了地形和土层对地震动的联合作用,揭示了强震区起伏地形场地震害加重及地震滑坡集中发生的原因所在.      

1985年云南禄劝地震余震尾波 Q 值的某些研究

建立在用来解释地方震和区域地震的尾波散射模型的基础上,用云南短周期区域台网12个台站的 VGK 地震仪记录的68个禄劝地震的余震序列尾波观测资料,分六个时间段估算了尾波 Q 值.发现在所观测的频率范围内(0.40——1.65Hz),Q 值与频率有关,估算的 Q 值在80——240之间,频率相关值=0.45.对表示散射强度的波源因子 B(fp)的估算结果多数在10-23——10-24量级,低 Q 值地区表现出激发强度高.值得注意的是禄劝地震前后,离震中区较近的台站所测得的 Q 值发生了明显变化,有的变化达2.0倍以上.Q 值随时间的区域变化,这在地震预报中也许有重要的实际意义.此外,测定了余震震源参数.对震源参数间的定量关系的讨论结果给出:IgE=1.59ML+11.335;E=(2.1010-5M0;在3.0ML5Pa,一个显著的特征是应力降与地震震级不存在明显的线性关系.     

Attenuation of Coda Waves at the Changbaishan Tianchi Volcanic Area in Northeast China

21 earthquakes recorded by a temporary seismic network in the Changbaishan Tianchi volcanic area in Northeast China operated during the summer of 2002 and 2003 were analyzed to estimate the S coda attenuation. The attenuation quality factor Q_c was estimated using the single scattering attenuation model of Sato (1977) in the frequency band from 4 to 24 Hz. All the events studied in this paper occurred at depths from 2 to 6 km with ML of 1.4–2.8. The epicentral distances are less than 25 km. For all events which occurred near the Tianchi Lake (caldera), the Q_c patterns obtained at the stations near the lake are similar, and the Q_c values are relatively small. At the stations located about 15 km east of the Tianchi Lake, however, the average Q_c is significantly higher. For an event which occurred 25km from the lake to the west, Q_c patterns derived at the stations near the lake are quite similar to the above mentioned Q_c for stations located in the east. Further study shows that Q_c value in the north and central areas of the volcano is relatively lower than that in the surrounding area. Compared to other volcanic areas in the world, the average Q_c of the Changbaishan Tianchi volcanic area is obviously lower. The deep seismic sounding and teleseismic receiver function studies indicated more than one lower velocity layer in the crust. The MT studies suggested the presence of high conductive bodies beneath the area. We interpret the strong attenuation of coda waves near the Changbaishan Tianchi volcano as being possibly related to high temperature medium caused by shallow magma chambers.     

Variations of local seismic response in Benevento (Southern Italy) using earthquakes and ambient noise recordings

The city of Benevento (Southern Italy) has been repeatedly struck by large historical earthquakes. A heterogeneous geologic structure and widespread soft soil conditions make the estimation of site effects crucial for the seismic hazard assessment of the city. From 2000 until 2004, we installed seismic stations to collect earthquake data over zones with different geological conditions. Despite the high level of urban noise, we recorded more than 150 earthquakes at twelve sites. This data set yields the first, well documented experimental evidence for weak to moderate local amplifications. We investigated site effects primarily by the classical spectral ratio technique (CSR) using a rock station placed on the Benevento hill as reference. All sites in the Calore river valley and in the eastern part of the Benevento hill show a moderate high-frequency (f > 4 Hz) amplification peak. Conversely, sites in the Sabato river valley share weak-to-moderate amplification in a wide frequency band (from 1–2 to 7–10 Hz), without evident frequency peaks. Application of no-reference-site techniques to earthquake and noise data confirms the results of the CSRs in the sites of the Calore river valley and of the eastern part of the Benevento hill, but fails in providing indications for site effects in the Sabato river valley, being the H/V ratios nearly flat. One-dimensional modeling indicates that the ground motion amplification can be essentially explained in terms of a vertically varying geologic structure. High-frequency narrow peaks are caused by the strong impedance contrast existing between near-surface soft deposits and stiff cemented conglomerates. Conversely, broad-band amplifications in the Sabato river valley are likely due to a more complex layering with weak impedance contrasts both in the shallow and deep structure of the valley.     

尾波Qc值随时间变化在地震预测中应用的研究

利用云南地区22个区域数字地震台站记录到的,自1999年下半年至2004年下半年的1403条数字地震波形资料,采用Sato单次散射模型,计算了各台站周围50km内的尾波Qc值;并给出了其中14个资料相对较为丰富的台站的尾波Qc值随时间变化结果。依据这14个台的3Hz、10Hz、18Hz频率点的2002年1月1日前后平均Qc值的变化,研究认为在未来几年内,云南地区比较危险的区域可能位于滇中块体的东南边界附近、滇中块体的西南部及其西边界附近。     

Comparison of amplitude decay rates in reflection,refraction, and local earthquake records

Three types of seismic data recorded near Coalinga, California were analyzed to study the behavior of scattered waves: 1) aftershocks of the May 2, 1983 earthquake, recorded on verticalcomponent seismometers deployed by the USGS; 2) regional refraction profiles using large explosive sources recorded on essentially the same arrays above; 3) three common-midpoint (CMP) reflection surveys recorded with vibrator sources over the same area. Records from each data set were bandpassed filtered into 5 Hz wide passbands (over the range of 1–25 Hz), corrected for geometric spreading, and fit with an exponential model of amplitude decay. Decay rates were expressed in terms of inverse codaQ (Q _c ^–1 ).Q _c ^–1 values for earthquake and refraction data are generally comparable and show a slight decrease with increasing frequency. Decay rates for different source types recorded on proximate receivers show similar results, with one notable exception. One set of aftershocks shows an increase ofQ _c ^–1 with frequency.Where the amplitude decay rates of surface and buried sources are similar, the coda decay results are consistent with other studies suggesting the importance of upper crustal scattering in the formation of coda. Differences in the variation ofQ _c ^–1 with frequency can be correlated with differences in geologic structure near the source region, as revealed by CMP-stacked reflection data. A more detailed assessment of effects such as the depth dependence of scattered contributions to the coda and the role of intrinsic attenuation requires precise control of source-receiver field geometry and the study of synthetic seismic data calculated for velocity models developed from CMP reflection data.     

Temporal Variations of Seismic Wave Velocity Associated with 1998 M6.1 Shizukuishi Earthquake

—?We attempt to detect temporal variations of seismic wave velocity before and after 1998 M6.1 Shizukuishi, northeastern Japan, earthquake by using waveform data from explosions and earthquake doublets spanning the period immediately before and after the earthquake. Direct P waves of the second explosion are delayed by ～20 ms at observation stations with epicentral distances less than 15 km. This tendency does not change if the analysis frequency band is changed. Our result suggests that the P-wave velocity decreased by at least 1% in the extremely shallow region of the hanging wall of the M6.1 thrust event after its occurrence. On the other hand, there was the frequency dependence of the coda wave delays for both artificial sources and for natural events. At 5–10 Hz, immediate sharp increases by more than 20 ms in time delays and lower coherency were observed at several stations. We estimated the region in which P-wave velocity might have decreased after the M6.1 earthquake. Maximum depth of the region is 13 km. The region includes the aftershock area of the M6.1 earthquake, but is eccentric to the earthquake in the west. Considering the frequency band analyzed (5–10 Hz), the scale of the spatial inhomogeneity which led to the coda wave delay is several hundreds meters. We investigated candidates for the cause of the direct P-wave and coda wave delay. Observed direct P-wave delay can be partly explained by the stress changes caused by coseismic fault slip. However, the coda wave delay cannot be explained by the stress changes that are limited to the superficial area. Crustal heterogeneity should have changed at coseismic time in the deeper area where aftershocks of the M6.1 earthquake occurred.     

Earthquake activity during the 1983 Etna eruption

Seismic stations, with automatic P-picking and satellite retransmission were set up on Mount Etna following the eruption started on March 1983. Positions of the stations were chosen in order to complement the permanent telemetered network of Catania University.Comparison between locations obtained by both networks were made for earthquakes recorded by at least 5 ARGOS DCP (Data Collection Platform) stations. We observed a satisfactory agreement for events inside both networks.By merging data of both networks, it has been possible to locate more than 50 earthquakes for which separate computation was not possible due to the low number of arrivals.On 3rd-4th June a swarm of deep seismic events was observed. Hypocenters of these earthquakes are clearly located in a NNW-SSE-trending vertical zone of 5 km width at a depth of 7–36 km.Changes in the distribution of shallow seismic activity, before and after this swarm, have been observed.