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.     

Seismotectonic re-evaluation of the 1976 Friuli, Italy, seismic sequence

The 1976 Friuli seismic sequence is known for itscomplexity, with several unusually large aftershocksoccurring in the epicentral area. The source region, where the northeastern part of the Alpine chain meetsthe northern Dinaric chain, is characterized by a complexcompressional tectonic regime. Previous studies have not clearly identified which of the two main tectonic systems was activated, in part due to the limited precision of derived earthquake parameters, such as hypocentral coordinates and focal mechanisms. We review the locations for the 10 largest events of the sequence, including theSeptember 16, 1977 and April 18, 1979 earthquakes andwe compute their centroid moment tensors. Source parameters are calculated using intermediate period surfacewaves and the modification of the Harvard centroid-moment tensoralgorithm proposed by Arvidsson and Ekström (1998). A summary of all available geological, geodetic and seismological data show that most of the earthquakes may be associated with the Periadriatic overthrust and other related thrust faults. Based on their locations and focal mechanisms only two ofthe early aftershocks (May 7, 1976 and May 9, 1976) appear instead to have occurred on Dinaric structures.     

Empirical seismic fragility for the precast RC industrial buildings damaged by the 2012 Emilia (Italy) earthquakes

The paper analyses the seismic fragility of precast reinforced concrete buildings using observational damage data gathered after the 2012 Emilia earthquakes that struck Northern Italy. The damage level in 1890 buildings was collected, classified and examined. Damage matrices were then evaluated, and finally, empirical fragility curves were fitted using Bayesian regression. Building damage was classified using a six‐level scale derived from EMS‐98. The completeness of the database and the spatial distribution of the buildings investigated were analysed using cadastral data as a reference. The intensity of the ground motion was quantified by the maximum horizontal peak ground acceleration, which was obtained from ShakeMaps. Copyright © 2017 John Wiley & Sons, Ltd.     

Simultaneous inversion of velocity structures and hypocentral locations: Application to the Friuli seismic area NE Italy

A layeredP- andS-wave velocity model is obtained for the Friuli seismic area using the arrival time data ofP- andS-waves from local earthquakes. A damped least-squares method is applied in the inversion.The data used are 994P-wave arrival times for 177 events which have epicenters in the region covered by the Friuli seismic network operated by Osservatorio Geofisico sperimentale (OGS) di Trieste, which are jointly inverted for the earthquake hypocenters andP-wave velocity model. TheS-wave velocity model is estimated on the basis of 978S-wave arrival times and the hypocenters obtained from theP-wave arrival time inversion. We also applied an approach thatP- andS-wave arrival time data are jointly used in the inversion (Roecker, 1982). The results show thatS-wave velocity structures obtained from the two methods are quite consistent, butP-wave velocity structures have obvious differences. This is apparent becauseP-waves are more sensitive to the hypocentral location thanS-waves, and the reading errors ofS-wave arrival times, which are much larger than those ofP-waves, bring large location errors in the joint inversion ofP- andS-wave arrival time. The synthetic data tests indicated that when the reading errors ofS-wave arrivals are larger than four times that ofP-wave arrivals, the method proposed in this paper seems more valid thanP- andS-wave data joint inversion. Most of the relocated events occurred in the depth range between 7 and 11 km, just above the biggest jump in velocity. This jump might be related to the detachment line hypothesized byCarulli et al. (1982). From the invertedP- andS-wave velocities, we obtain an average value 1.82 forV _p /V _s in the first 16 km depth.     

Response spectral attenuation relationships for Sumatran-subduction earthquakes and the seismic hazard implications to Singapore and Kuala Lumpur

Singapore and Kuala Lumpur, the capital of Malaysia, may well represent the classic examples of area with low seismic hazard but with high consequence. Both cities are located in a low-seismicity region of Southeast Asia, where active seismic sources are located more than 300 km away. Seismic designs have not been implemented in this seemingly low-hazard region though distant earthquakes in Sumatra had frequently shaken high-rise structures in the two cities. Several studies have been conducted to systematically assess the seismic hazards of Singapore and the Malay Peninsula. The present research particularly addresses issues in deriving a new set of attenuation relationships of peak ground acceleration (PGA), peak ground velocity (PGV) and response spectral acceleration (RSA) for the Sumatran-subduction earthquakes. To be relevant for the seismic hazard assessment of the remote metropolises, the derived attenuation relationships cover a long distance range from 150 to 1500 km. The attenuation relationships are derived using synthetic seismograms that account for source and path effects. The uncertainties in rupture parameters, such as stress drop, strike, dip and rake angles, have been defined according to the regional geological and tectonic settings as well as the ruptures of previous earthquakes. The seismic potential of the Sumatran subduction zone are high in the region from 2°N to 5°S as there has been no recurrence of great thrust events since 1861. A large event with M_w greater than 7.8 in this particular subduction zone may be capable of generating destructive ground motions in Singapore and Kuala Lumpur, even at a distance of 700 km.     

河南平顶山平煤矿区天然地震、爆破、塌陷时频特征分析

选取2019年3月—8月河南平顶山市宝丰县平煤矿区发生的ML 2.0—2.9天然地震、爆破、塌陷等9次震动事件,在区域地质构造背景和波形特征分析基础上,采用短时傅里叶变换(STFT)方法开展时频波谱分析,提取不同类型事件的时频特征。结果显示:(1)天然地震频率成分丰富,且高、低频分布均匀,P波在约3 Hz和8 Hz处存在2个峰值,S波存在多个峰值;(2)爆破事件的时频谱相对集中,以低频为主,P波频率峰值约5Hz,信号主频随时间变化,大致呈线性降低至1—2 Hz;(3)塌陷事件频率成分以4 Hz以下的低频为主,P波无明显峰值且频率成分单一,主频出现在2 Hz左右的面波。本文结果可为今后该矿区震动事件类型判断提供参考依据。     

Improved 2-D attenuation analysis for Northern Italy using a merged dataset from selected regional seismic networks

A merged, high-quality waveform dataset from different seismic networks has been used to improve our understanding of lateral seismic attenuation for Northern Italy. In a previous study on the same region, Morasca et al. (Bull Seismol Soc Am 98:1936–1946, 2008) were able to resolve only a small area due to limited data coverage. For this reason, the interpretation of the attenuation anomalies was difficult given the complexity of the region and the poor resolution of the available data. In order to better understand the lateral changes in the crustal structure and thickness of this region, we selected 770 earthquakes recorded by 54 stations for a total of almost 16,000 waveforms derived from seismic networks operating totally or partially in Northern Italy. Direct S-wave and coda attenuation images were obtained using an amplitude ratio technique that eliminates source terms from the formulation. Both direct and early-coda amplitudes are used as input for the inversions, and the results are compared. Results were obtained for various frequency bands ranging between 0.3 and 25.0 Hz and in all cases show significant improvement with respect to the previous study since the resolved area has been extended and more crossing paths have been used to image smaller scale anomalies. Quality-factor estimates are consistent with the regional tectonic structure exhibiting a general trend of low attenuation under the Po Plain basin and higher values for the Western Alps and Northern Apennines. The interpretation of the results for the Eastern Alps is not simple, possibly because our resolution for this area is still not adequate to resolve small-scale structures.     

Investigation on regional attenuation of Vrancea (Romania) intermediate-depth earthquakes

This paper aims at investigating possible regional attenuation patterns in the case of Vrancea(Romania) intermediate-depth earthquakes.Almost 500 pairs of horizontal components recorded during 13 intermediate-depth Vrancea earthquakes are employed in order to evaluate the regional attenuation patterns.The recordings are grouped according to the azimuth with regard to the Vrancea seismic source and subsequently,Q models are computed for each azimuthal zone assuming similar geometrical spreading.Moreover,the local soil amplification which was disregarded in a previous analysis performed for Vrancea intermediate-depth earthquakes is now clearly evaluated.The results show minor differences between the four regions situated in front of the Carpathian Mountains and considerable differences in attenuation of seismic waves between the forearc and backarc regions(with regard to the Carpathian Mountains).Consequently,an average Q model of the type Q(f) = 115×f~(1.25) is obtained for the four forearc regions,while a separate Q model of the type Q(f) = 70×f~(0.90) is computed for the backarc region.These results highlight the need to evaluate the seismic hazard of Romania by using ground motion models which take into account the different attenuation between the forearc/backarc regions.     

A geological model for the Colfiorito earthquakes (September-October 1997, central Italy)

In this study, surface and subsurface geologicaldata are integrated with seismological data in orderto reconstruct a structural model for theSeptember-October 1997 Colfiorito earthquakes. Theseismic sequence is mainly controlled by two majorSW-dipping normal faults outcropping in the area (M.Pennino-M. Prefoglio and M.Civitella-Preci faults).The activated faults detach, at depth, on a commoneast-dipping low-angle normal fault, the AltotiberinaFault (AF). The AF is interpreted as the base of anactive hangingwall block which is stretching towardNE. The decrease in maximum depth of the earthquakefoci from the Colfiorito area (about 8 km) to theSellano area (about 6 km), suggested by the available seismological data, could be related to the eastward-deepening geometry of the AFdetachment. The seismic fault planes, inferred fromfocal mechanisms and aftershock distributions, arecharacterised by a moderate dip (average 40^°)toward SW, which appears to be independent from thepresence of pre-existing thrust planes.     

美国Cove Fort-Sulphurdale地热区地震震源机制和各向异性研究

Cove Fort-Sulphurdale位于美国犹他州中西部, 是一个地质特征独特的地热资源富集区． 为了有效开发地热资源, 有必要了解该地区的应力分布和裂缝分布． 本文在前人工作的基础上使用地震资料, 利用波形匹配、 P波极性和横纵波振幅比联合约束反演的方法, 研究该地区地震的震源机制解; 同时结合该地区的横波分裂分析研究, 确定该地区的应力分布和断层/裂隙走向． 波形反演结果显示, 实际数据与模拟数据的波形匹配度非常高, 且二者的P波初动极性和横纵波振幅比也都非常相近． 由震源机制分析得出的断层面走向大部分趋于南北向, 与美国区域应力图显示的最大水平主应力指向(南北方向), 具有较好的一致性, 符合目前研究对该地区的构造认知． 横波分裂分析结果也表明如果各向异性主要是由该地区的断裂构造所引起, 那么该地区的主要断层或裂缝的方向可能为南北向．      

Coda Q c Attenuation and Source Parameter Analysis in Friuli (NE Italy) and its Vicinity

—?The digital data acquired by 16 short-period seismic stations of the Friuli-Venezia-Giulia seismic network for 56 earthquakes of magnitude 2.3–4.7 which occurred in and near NE Italy have been used to estimate the coda attenuation Q _c and seismic source parameters. The entire area under study has been divided into five smaller regions, following a criterion of homogeneity in the geological characteristics and the constrains imposed by the distribution of available events. Standard IASPEI routines for coda Q _c determination have been used for the analysis of attenuation in the different regions showing a marked anomaly in the values measured across the NE border between Friuli and Austria for Q ₀ value. A large variation exists in the coda attenuation Q _c for different regions, indicating the presence of great heterogeneities in the crust and upper mantle of the region. The mean value of Q _c (f) increases from 154–203 at 1.5?Hz to 1947–2907 at 48?Hz frequency band with large standard deviation estimates.¶Using the same earthquake data, the seismic-moment, M ₀, source radius, r and stress-drop, Δσ for 54 earthquakes have been estimated from P- and S-wave spectra using the Brune's seismic source model. The earthquakes with higher stress-drop (greater than 1?Kbar) occur at depths ranging from 8 to 14?km.     

Magnitudes and travel-time residuals of the Friuli earthquakes between May and August 1976

Summary The records of the first series of the North-Italian aftershock series from the period May–August 1976 at the stations Prhonice and Kaperské Hory were analyzed from the point of view of the verification of the travel-time tables used for the interpretation of near earthquake seismograms and of the determination of the short-period body wave magnitudes for near shocks. Also the time-energy distribution was investigated for a 90 days period.     

Strong motion recorded during the Emilia 2012 thrust earthquakes (Northern Italy): a comprehensive analysis

A complex seismic sequence characterised by two thrust earthquakes of magnitudes M \(_\mathrm{L}\) 5.9 and M \(_\mathrm{L}\) 5.8 occurred on May 20 and 29, 2012, respectively, and activated the central portion of the Ferrara Arc structure beneath the Po Plain in northern Italy. The sequence, referred to as Emilia 2012, was recorded by the Italian Strong Motion Network, the Irpinia Network, the Friuli Venezia Giulia Network and 15 temporary stations installed by the Civil Protection Department. In this study, we compile and analyse a large dataset that contains 3,273 waveforms from 37 \(M_\mathrm{L} \ge 4.0\) seismic events. The main aim of this paper is to characterise the ground motion induced by the Emilia 2012 seismic sequence and compare it with other data in the Italian strong motion database and to the recent Ground Motion Prediction Equations (GMPEs) developed for northern Italy, all of Italy and Europe. This is achieved by (1) the computation and analysis of the strong motion parameters of the entire Emilia Strong Motion Dataset (ESMD) and (2) a comprehensive investigation of the May 29 event recordings in terms of time–frequency analysis, the ground motion parameters and the response spectra. This detailed analysis was made possible by the temporary Civil Protection Department stations that were installed soon after the May 20 event at several municipalities in the epicentral area. Most of the recordings are characterised by low-frequency content and long durations, which is a result of the thick sedimentary cover that is typical of the Po Plain. The distributions of the observed horizontal peak ground accelerations and velocities (PGAs and PGVs) with distance are generally consistent with the GMPEs. This is particularly true for the data from M \(_\mathrm{L} \ge \) 5.0 (M \(_\mathrm{W}\ge \) 5.0) events, though the data are scattered at distances beyond approximately 60–70 km and show faster attenuation than the European GMPEs. The horizontal components for the May 29 event at two near-fault stations (Mirandola and San Felice sul Panaro) are overestimated by all of the analysed GMPEs. In contrast, the vertical components, which played an important role in the shaking near the source, are underestimated. The May 29 event produced intense velocity pulses on the horizontal components and the highest peak ground acceleration ever recorded in Italy on the vertical component of the Mirandola near-fault station. The ground motion recordings contained in the ESMD significantly enrich the Italian strong motion database. They contribute new information about (1) the possibility of exceeding the largest recorded PGA in Italy, (2) the development of a spectral design that takes into account the role of the vertical component and the extreme variability of the near-fault ground shaking, and (3) the characterisation of the ground motions in deep sedimentary basins.     

我国大陆地区常用浅壳地震的地震动参数衰减关系

对地震动参数衰减关系应用时涉及的地震动参数、地震参数、基本模型、回归方法等方面的要点进行了阐述。以华北地区为例,简要介绍了数十年我国工程抗震设防应用中常用的几组浅壳地震峰值加速度、加速度反应谱衰减关系的主要特点,并对其差异进行了比较分析。最后总结了我国现阶段衰减关系研究的不足,对其未来发展方向进行了展望。     

Shakemaps uncertainties and their effects in the post-seismic actions for the 2012 Emilia (Italy) earthquakes

The May 20, 2012, Emilia Ml 5.9 earthquake was followed by some major aftershocks, well recorded by a large number of temporary stations that were installed to monitor the sequence. These additional recordings allowed us a thorough testing of the performance of the ShakeMap—a procedure designed to provide rapid information on the experienced ground motion. We found that the shakemaps for the May 29, 2012, Ml 5.8 earthquake, obtained using the permanent stations only, underestimate significantly the ground motion computed with the highest station density, especially for PSA at long periods (T \(=\) 3.0 s). This low-frequency motion is controlled primarily by the surface waves recorded in the Po plain: the observed site effects are likely not accounted properly by the site correction coefficient based on Vs30 as implemented in the ShakeMap procedure. The shakemaps determined during the seismic sequence have been included in an Italian national law that was passed after the 2012 earthquake. According to this law, the factories safety verifications were bound to the comparison between the shakemaps and the design acceleration required by the current national seismic code. We then decide to appraise the impact of the shakemap accuracy on the law provisions. Following the law recommendations, we have estimated the possible errors resulting from the incomplete evaluation of the ground shaking: our results show that, if the complete dataset were available at the time of the law approval, the number of buildings for which the safety check was required would have been significantly smaller.     

Remarks about magnitudes and maximum intensities of the earthquakes in Friuli, 1976

Summary The magnitudesM and maximum intensitiesI ₀ of 66 earthquakes which occurred in the period from 1901 to 1973 in the eastern alpine area are compared with those of 62 earthquakes occurred in Friuli, northern Italy between 6 May and 15 September 1976. The average focal depth of the Friuli events is about 10 km. The existing empirical relation betweenM andI ₀ for this area is improved by many new data.     

Local earthquakes seismic tomography in the Betic Cordillera (southern Spain)

A crustal tomographic image, from the surface down to 35 km depth beneath the Betic Cordillera (southern Spain), is obtained using data on local earthquakes recorded at stations from the National and Andalusian Seismic Networks. The velocity structure and the hypocentre locations are derived from the inversion of P first arrival times, using an iterative simultaneous inversion method. The reliability of the results is assessed using different control parameters. The inverted velocity field in the uppermost layers shows a significant lateral variability which reflects most of the large-scale geological features of the Betic Cordillera. Well determined local surface anomalies allow to constrain the location and geometry of the most prominent Neogene sedimentary basins. The upper crust is well resolved throughout the whole region, and is characterized by relatively high velocities in the Internal Betics and in the South Iberian Massif and lower velocities within the External Betics. A relatively well constrained event cluster displays a NNE–SSW trend, and outlines the contact zone between the Internal and the External domains. The middle and lower crustal levels show reliable results beneath the central part of the Betic Cordillera. High averaged velocities are obtained within the South Iberian and the Alboran domains, in contrast to a relatively low velocity anomaly which characterizes the boundary between them. These findings support the hypothesis of the lack of well differentiated crustal levels below the contact zone, while crustal layering is better defined beneath the Alboran and the Iberian domains.     

Reliability of the automatic procedures for locating earthquakes in southwestern Alps and northern Apennines (Italy)

Reliable automatic procedure for locating earthquake in quasi-real time is strongly needed for seismic warning system, earthquake preparedness, and producing shaking maps. The reliability of an automatic location algorithm is influenced by several factors such as errors in picking seismic phases, network geometry, and velocity model uncertainties. The main purpose of this work is to investigate the performances of different automatic procedures to choose the most suitable one to be applied for the quasi-real-time earthquake locations in northwestern Italy. The reliability of two automatic-picking algorithms (one based on the Characteristic Function (CF) analysis, CF picker, and the other one based on the Akaike’s information criterion (AIC), AIC picker) and two location methods (“Hypoellipse” and “NonLinLoc” codes) is analysed by comparing the automatically determined hypocentral coordinates with reference ones. Reference locations are computed by the “Hypoellipse” code considering manually revised data and tested using quarry blasts. The comparison is made on a dataset composed by 575 seismic events for the period 2000–2007 as recorded by the Regional Seismic network of Northwestern Italy. For P phases, similar results, in terms of both amount of detected picks and magnitude of travel time differences with respect to manual picks, are obtained applying the AIC and the CF picker; on the contrary, for S phases, the AIC picker seems to provide a significant greater number of readings than the CF picker. Furthermore, the “NonLinLoc” software (applied to a 3D velocity model) is proved to be more reliable than the “Hypoellipse” code (applied to layered 1D velocity models), leading to more reliable automatic locations also when outliers (wrong picks) are present.