Seismic response of deep Quaternary sediments in historical center of L’Aquila City (central Italy)

On 6 April 2009 a M_w=6.1 earthquake produced severe destruction and damage over the historic center of L’Aquila City (central Italy), in which the accelerometer stations AQK and AQU recorded a large amount of near-fault ground motion data. This paper analyzes the recorded ground motions and compares the observed peak accelerations and the horizontal to vertical response spectral ratios with those revealed from numerical simulations. The finite element method is considered herein to perform dynamic modeling on the soil profile underlying the seismic station AQU. The subsurface model, which is based on the reviewed surveys that were carried out in previous studies, consists of 200–400 m of Quaternary sediments overlying a Meso-Cenozoic carbonate bedrock. The Martin-Finn-Seed's pore-water pressure model is used in the simulations. The horizontal to vertical response spectral ratio that is observed during the weak seismic events shows three predominant frequencies at about 14 Hz, 3 Hz and 0.6 Hz, which may be related to the computed seismic motion amplification occurring at the shallow colluvium, at the top and base of the fluvial-lacustrine sequence, respectively. During the 2009 L’Aquila main shock the predominant frequency of 14 Hz shifts to lower values probably due to a peculiar wave-field incidence angle. The predominant frequency of 3 Hz shifts to lower values when the earthquake magnitude increases, which may be associated to the progressive softening of soil due to the excess pore-water pressure generation that reaches a maximum value of about 350 kPa in the top of fluvial-lacustrine sequence. The computed vertical peak acceleration underestimates the experimental value and the horizontal to vertical peak acceleration ratio that is observed at station AQU decreases when the earthquake magnitude increases, which reveals amplification of the vertical component of ground motion probably due to near-source effects.     

A spatial correlation model of peak ground acceleration and response spectra based on data of the Istanbul Earthquake Rapid Response and Early Warning System

Ground motion intensity measures such as the peak ground acceleration (PGA) and the pseudo-spectral acceleration (PSA) at two sites due to the same seismic event are correlated. The spatial correlation needs to be considered when modeling ground-motion fields for seismic loss assessments, since it can have a significant influence on the statistical moments and probability distribution of aggregated seismic loss of a building portfolio.Empirical models of spatial correlation of ground motion intensity measures exist only for a few seismic regions in the world such as Japan, Taiwan and California, since for this purpose a dense observation network of earthquake ground motion is required. The Istanbul Earthquake Rapid Response and Early Warning System (IERREWS) provides one such dense array with station spacing of typically 2 km in the urban area of Istanbul. Based on the records of eight small to moderate (M_w3.5–M_w5.1) events, which occurred since 2003 in the Marmara region, we establish a model of intra-event spatial correlation for PGA and PSA up to the natural period of 1.0 s.The results indicate that the correlation coefficients of PGA and short-period PSA decay rapidly with increasing interstation distance, resulting in correlation lengths of approximately 3–4 km, while correlation lengths at longer natural periods (above 0.5 s) exceed 6 km. Finally, we implement the correlation model in a Monte Carlo simulation to evaluate economic loss in Istanbul's district Zeytinburnu due to a M_w7.2 scenario earthquake.     

Experimental investigation into ground vibrations induced by very high speed trains on a non-ballasted track

A field measurement of ground vibration was performed on the Beijing−Shanghai high-speed railway in China. In this paper, the experimental results of vertical ground vibration accelerations induced by very high speed trains running over a non-ballasted track on embankment with speeds from 300 to 410 km/h are reported and analyzed in detail for the first time. Characteristics of ground vibration accelerations in both time and frequency domains are analyzed based on the test data. It is shown that the periodic exciting action of high-speed train bogies can be identified in time histories of vertical accelerations of the ground within the range of 50 m from the track centerline. The first dominant sensitive frequency of the ground vibration acceleration results from the wheelbase of the bogie, and the center distance of two neighboring cars plays an important role in the significant frequencies of the ground vibration acceleration. Variations of time–response peak value and frequency-weighted vertical acceleration level of ground vibration in relation with train speed as well as the distance from the track centerline are also investigated. Results show that the time-domain peak value of ground vibration acceleration exhibits an approximately linear upward tendency with the increase of train speed. With the increasing distance from the track centerline, the frequency-weighted vertical acceleration level of the ground vibration attenuates more slowly than the time-domain peak value of the ground vibration acceleration does. Severe impact of high-speed railway ground vibration on human body comfort on the ground occurs at the speed of 380–400 km/h. The results given in the paper are also valuable for validating the numerical prediction of train induced ground vibrations.     

面向地震烈度仪的现地地震动预测模型的构建与验证

利用初期P波预警参数构建现地地震动预测模型,使其在达到设定阈值时快速发出报警信息,是现地地震预警系统面临的一个关键问题,直接关系到发布信息的准确性和及时性。针对地震烈度仪基于微机电系统传感器记录到的数据质量较差,通过两次积分获取的位移存在较大偏差,会引起更多的误报和漏报,本文采用不同阶数(1—4阶)的巴特沃斯滤波器,分别构建了基于P波3 s和全P波段数据的位移幅值PD、速度幅值PV和加速度幅值PA与地震动峰值速度PGV和峰值加速度PGA的现地地震动预测模型,然后利用收集到的川滇示范预警网地震事件记录进行验证。结果表明,对于地震烈度仪微机电系统传感器的记录,采用1阶巴特沃斯滤波器处理、基于全P波段波形拟合获取到的PV与PGV的相关性和PA与PGA的相关性为两种最优现地地震动预测模型。具体应用时,应同时利用两种或两种以上的统计关系进行现地地震动预测,并将实际地震动观测值作为额外的判定条件,以降低误报率和漏报率。     

Earthquake induced liquefaction features in the Karewas of Kashmir Valley North-West Himalayas,India: Implication to paleoseismicity

Earthquake induced liquefaction features found in the Karewas of Kashmir Valley are potential tools for estimating energy center, magnitude and peak ground acceleration of the paleoearthquakes. Size, pattern and spatial distribution of liquefaction features and in-situ geotechnical data were collected at the paleoliquefaction sites to assess the magnitude and peak ground acceleration of the paleoearthquake. The magnitude of paleoearthquakes was estimated to be of the order of (Mw = 6.2) using “Magnitude bound method” and “In-situ testing of geotechnical properties and assessment of liquefaction potential of liquefied beds (source stratum) using cyclic stress method”. The peak ground acceleration computed were in the range of 0.27 g to 0.83 g using cyclic stress method and 0.30 g using attenuation equation NDMA [21].     

Crustal velocity structure in the southern edge of the Central Alborz (Iran)

Inversion of local earthquake travel times and joint inversion of receiver functions and Rayleigh wave group velocity measurements were used to derive a simple model for the velocity crustal structure beneath the southern edge of the Central Alborz (Iran), including the seismically active area around the megacity of Tehran. The P and S travel times from 115 well-located earthquakes recorded by a dense local seismic network, operated from June to November 2006, were inverted to determine a 1D velocity model of the upper crust. The limited range of earthquake depths (between 2 km and 26 km) prevents us determining any velocity interfaces deeper than 25 km. The velocity of the lower crust and the depth of the Moho were found by joint inversion of receiver functions and Rayleigh wave group velocity data. The resulting P-wave velocity model comprises an upper crust with 3 km and 4 km thick sedimentary layers with P wave velocities (V_p) of ～5.4 and ～5.8 km s^?1, respectively, above 9 km and 8 km thick layers of upper crystalline crust (V_p ～6.1 and ～6.25 km s^?1 respectively). The lower crystalline crust is ～34 km thick (V_p ～ 6.40 km s^?1). The total crustal thickness beneath this part of the Central Alborz is 58 ± 2 km.     

Impact and cyclic shaking on loose sand properties in laminar box using gap sensors

This paper focuses on using high-frequency GAP-SENSORs (GSs), accelerometers, and load cells in a laminar shear box (LSB) filled with loose Toyoura sand to understand the effects of impact loads and cyclic shaking at 1-G on soil properties. The shear wave velocity at small strain (V_s) was calculated directly from first arrival reference using displacement time-history of two GSs under impact loading. Moreover, from first peak using the reduced deformation amplitude technique, damping ratio was calculated. In addition, shaking table tests were performed under harmonic loading with amplitude of acceleration inside the model ground varying from 0.02 g to 1 g. The frequencies of excitation varied from 1 Hz to 10 Hz. GSs and inside accelerometers were used to directly measure the outside lateral deformation and shear stress at different elevations of LSB, respectively. Results show that the shear modulus (G) and the damping ratio (D) behavior of model sand are generally consistent with the behavior presented by similar tests using only accelerometers. In addition, damping ratio increases as frequency loading increases. Characteristic changes in two shear stress components in shaking loading conditions were also investigated using high precision inside load cells.     

Sensitivity of ground-motion scenarios to earthquake source parameters in the Tehran metropolitan area,Iran

The northern Tehran fault (NTF) is a principal active fault of the Alborz mountain belt in the northern Iran. The fault is located north of the highly populated Metropolitan Area of Tehran. Historical records and paleoseismological studies have shown that the NTF poses a high seismic risk for the Tehran region and the surrounding cities (e.g. Karaj). A series of ground-motion simulations are carried out using a hybrid kinematic-stochastic model to calculate broadband (0.1–20 Hz) ground-motion time histories for deterministic earthquake scenarios (M7.2) on the NTF. We will describe the source characteristics of the target event to develop a list of scenario earthquakes that are probably similar to a large earthquake on the NTF. The effect of varying different rupture parameters such as rupture velocity and rise time on the resulting broadband strong motions has been investigated to evaluate the range of uncertainty in seismic scenarios. The most significant parameters in terms of ground-shaking level are the rise time and the value of the rupture velocity. For the worst-case scenario, the maximum expected horizontal acceleration, and velocity at rock sites in Tehran range between 128 and 1315 cm/s/s and 11–191 cm/s, respectively. For the lowest scenario, the corresponding values range between 102 and 776 cm/s/s and 12 to 81 cm/s. Nonlinear soil effects may change these results but are not accounted for in this study. The largest variability of ground motion is observed in neighborhood of asperity and also in the direction of rupture propagation. The calculated standard deviation of all ground-motion scenarios is less than 30% of the mean. The capability of the simulation method to synthesize expected ground motions and the appropriateness of the key parameters used in the simulations are confirmed by comparing the synthetic peak ground motions (PGA, PGV and response spectra) with empirical ground-motion prediction equations.     

Exploring the feasibility of earthquake early warning using records of the 2008 Wenchuan earthquake and its aftershocks

Earthquake early warning system (EEWS) is one of the effective ways to mitigate earthquake damage and can provide few seconds to tens of seconds of advanced warning time of impending ground motions, allowing for mitigation measures to be taken in the short term. After the devastating Ms8.0 Wenchuan earthquake, hundreds of M4-6 earthquakes occurred with depth range of 2–24 km. We explore a practical approach to earthquake early warning in Wenchuan area by determining a ground-motion period parameter τ_c and a high-pass filtered vertical displacement amplitude parameter Pd from the initial 3 s of the P waveforms of these aftershocks with M≥4.0. The empirical relationships both between τ_c and M, and between Pd and peak ground velocity PGV for the Wenchuan area are presented. The τ_c result shows that it is systematically greater for slow earthquakes, leading to a possible false alarm. The moment rate function is used to handle the fact that the Pd parameter alone miss the M=8.0 mainshock. These two relationships can be used to detect the occurrence of a major earthquake and provide onsite warning in the area around the station where onset of strong ground motions is expected within seconds after the arrival of the P wave. The robustness of onsite early warning can be increased by using multistation data when the station density is high or by combing τ_c and Pd as a single indicator.     

Seismic microzonation studies for the city of Ragusa (Italy)

The seismic history of the city of Ragusa (Italy), the geotechnical characterisation of the subsoil and the site response analysis should be correctly evaluated for the definition of the Seismic Geotechnical Hazard of the city of Ragusa, through geo-settled seismic microzoning maps. Basing on the seismic history of the city of Ragusa, the following earthquake scenarios have been considered: the “Val di Noto” earthquake of January 11, 1693 (with intensity X–XI on MCS scale, magnitude M_W=7.41 and epicentral distance of about 53 km); the “Etna” earthquake of February 20, 1818 (with intensity IX on MCS scale, magnitude M_W=6.23 and epicentral distance of about 64 km); the Vizzini earthquake of April 13, 1895 (with intensity I=VII–VIII on MCS scale, magnitude M_W=5.86 and epicentral distance of about 26 km); the “Modica” earthquake of January 23, 1980 (with intensity I=V–VI on MCS scale, magnitude M_W=4.58 and epicentral distance of about 10 km); the “Sicilian” earthquake of December 13, 1990 (with intensity I=VII on MCS scale, magnitude M_W=5.64 and epicentral distance of about 50 km). Geotechnical characterisation has been performed by in situ and laboratory tests, with the definition of shear wave velocity profiles in the upper 30 m of soil. Soil response analyses have been evaluated for about 120 borings location by some non-linear 1-D models. Finally the seismic microzonation of the city of Ragusa has been obtained in terms of maps with different peak ground acceleration at the surface; shaking maps for the central area of the city of Ragusa were generated via GIS for the earthquake scenarios.     

Characteristics of the strong ground motions from the 6 April 2009 L’Aquila earthquake,Italy

An M_w 6.25 earthquake occurred on April 6, 2009 at 03:33 a.m. local time, in the Abruzzo region (Central Italy), close to the city of L’Aquila. The earthquake ruptured a North-West (NW)–South-East (SE) oriented normal fault dipping toward the South-West (SW), with the city of L’Aquila lying a few kilometers away on the hanging wall.The main shock has been recorded by fifty-eight accelerometric stations: the highest number of digital recordings ever obtained in Italy for a single earthquake, one of the best-recorded earthquakes with a normal fault mechanism. Very high values of peak ground acceleration (0.3–0.65 g) were observed close to the center of L’Aquila (6 stations at zero JB distance from the fault). The earthquake caused severe loss of lives (299 victims and 1500 injured) and damage (about 18000 unusable buildings) in the epicentral area.In this study we analyze the ground motion characteristics of both the main shock in terms of peak ground acceleration (PGA), peak ground velocity (PGV), and pseudo-acceleration response spectra (5% of damping ratio). In particular, we compare the pseudo-acceleration response spectra for horizontal directions with the EC8 design spectrum and the new Italian building code (NTC08). In order to understand the characteristics of the ground motions induced by L’Aquila earthquake, we also study the source-related effects and site response of the strong motion stations that recorded the seismic sequence. A novel method is used for the analysis of inter-station and site-specific H/V spectral ratios for the main event and for 12 aftershocks.     

Seismic hazard evaluation in Anjar city area of western India: Microtremor array measurement

In western India during the Bhuj earthquake (Mw 7.6) on January 26, 2001, the Anjar City at ~30 km southwest of Bhuj experienced three types of damage scenario: severely damaged, less damaged and non-damaged. Similar damage patterns were also observed for the 1819 (Mw 7.8) and the 1956 (Mw 6.0) earthquakes. Microtremor array measurements were conducted in and around the Anjar city to examine the strength of soil structures and damage pattern. Significant differences are observed in frequencies and amplitudes in horizontal-to-vertical spectral ratio (HVSR) using microtremor measurements. The severely- damaged site shows two peak amplitudes: 2.8 at 1.2 Hz; and 4.0 at 8.0 Hz. The less-damaged site also shows two amplitudes: 2.5 and 2.1 at 1.4 Hz; and 2.0 Hz, respectively. The non-damaged site, on the other hand, shows that the HVSR curves become almost flatter. Similar results for three types of damage scenario based on analyses of earthquake records are also observed for the study area. The microtremor array measurements has revealed shear wave velocity V_s≥400 m/s at 18 m depth in the non-damaged, at 40 m in the less-damaged and at 60 m depth in the severely-damaged sites. The site amplitudes and the V_s values show a good correlation with the soil characteristics and damage pattern, suggesting that strength of soil layers at varying depths is a dictating factor for the estimate of the earthquake risk evaluation of the area under study.     

Shake table tests of a new steel–asphalt composite layer system for the seismic base isolation of housing units

Base isolation is an effective way to reduce earthquake energy transfer from ground to structure, but existing seismic isolation systems are not very suitable for rural buildings for some reasons. A new steel–asphalt composite layer for the seismic base isolation of housing units is present in this paper. Its dynamic characteristics and isolation effect are studied by shake table tests of two full-scale specimens. Different earthquake waves with different peak ground accelerations (PGA, from 0.1 g to 0.4 g) are input. Test results show that the isolation layer could efficiently reduce the input acceleration. Moreover, as the PGA increase, the isolation layer shows good function of the displacement limit.     

Artificial earthquake test of gas supply networks

In this paper, seismic behavior of gas supply networks is investigated by constructing a 24 m×24 m buried network and conducting an artificial earthquake test by detonating TNT explosives. The pipe network consists of typical welded steel pipes used in gas supply networks. Ground motion is produced by detonating 30 kg of TNT explosives buried 5 m below ground surface and 15 m away from the network. The test preparations are presented, including site selection, network layout, ground motion generation, and arrangements of different sensors. The measured ground acceleration, strain, and acceleration of welded steel pipes, and pipe–soil relative slippage are discussed. The deformation patterns of welded steel pipes are also analyzed and explained, including axial and bending deformations, systematic network characteristics, deformation relationship between ground and pipes, and dynamic response of buried pipes.     

Seismic microzoning from synthetic ground motion earthquake scenarios parameters: The case study of the city of Catania (Italy)

The city of Catania (Italy) in the South-Eastern Sicily has been affected in past times by several destroying earthquakes with high values of estimated magnitude. The seismogenic area to the south of Volcano Etna, known as Iblean Area, is placed between the African and the Euro-Asiatic plates on the west of the Ibleo-Maltese escarpment, to the south of the Graben of the Sicilian channel and on the east of the overlapping front of Gela. Basing on the seismic history of Catania, the following earthquake scenarios have been considered: the “Val di Noto” earthquake of January 11, 1693 (with intensity X-XI on MCS scale, magnitude M_W=7.41 and epicentral distance of about 13 km); the “Etna” earthquake of February 20, 1818 (with intensity IX on MCS scale, magnitude M_W=6.23 and epicentral distance of about 10 km). The soil response analysis at the surface, in terms of time history and response spectra, has been obtained by 1-D equivalent linear models for about 1200 borings location available in the data-bank of the central area of Catania of about 50 km², using deterministic design scenario earthquakes as input at the conventional bedrock.Seismic microzoning maps of the city of Catania have been obtained in terms of different peak ground acceleration at the surface and in terms of amplification ratios for given values of frequency.     

A new fault rupture scenario for the 2003 Mw 6.6 Bam earthquake,SE Iran: Insights from the high-resolution QuickBird imagery and field observations

The December 26, 2003 M_w 6.6 Bam earthquake is one of the most disastrous earthquakes in Iran. QuickBird panchromatic and multispectral satellite imagery with 61 cm and 2.4 m ground resolution, respectively provide new insights into the surface rupturing process associated with this earthquake. The results indicate that this earthquake produced a 2–5 km-wide surface rupture zone with a complex geometric pattern. A 10-km-long surface rupture zone developed along the pre-existing Bam fault trace. Two additional surface rupture zones, each 2–5 km long, are oblique to the pre-existing Bam fault in angles of 20–35°. An analysis of geometric and geomorphic features also shows that movement on the Bam fault is mainly right-lateral motion with some compressional component. This interpretation is consistent with field investigations, analysis of aftershocks as well as teleseismic inversion. Therefore, we suggest that the 2003 Bam earthquake occurred on the Bam fault, and that the surface ruptures oblique to the Bam fault are caused by secondary faulting such as synthetic shears (Reidel shears). Our fault model for the Bam earthquake provides a new tectonic scenario for explaining complex surface deformations associated with the Bam earthquake.     

Ground response analysis of Kanpur soil along Indo-Gangetic Plains

The present work deals with 1D and 2D ground response analysis and liquefaction analysis of alluvial soil deposits from Kanpur region along Indo-Gangetic plains. Standard penetration tests and seismic down hole tests have been conducted at four locations namely IITK, Nankari village, Mandhana and Bithoor at 1.5 m interval up to a depth of 30 m below the ground surface to find the variation of penetration blows and the shear wave velocity along the depth. From the selected sites undisturbed as well as representative soil samples have been collected for detailed soil classification. The soil profiles from four sites have been considered for 1D and 2D ground response analysis by applying the free field motions of three Himalayan earthquakes namely Chamba earthquake (M_w—5.1), Chamoli earthquake (M_w—6.4) and Uttarkashi earthquake (M_w—6.5). An average value of Peak Ground Acceleration (PGA) obtained from 1D and 2D analysis is considered for liquefaction analysis and post-liquefaction settlement. The excess pore water pressure ratio is greater than 0.8 at a depth of 24 m from ground surface for IITK, Nankari village, Bithoor sites. More than 50% of post liquefaction settlement is contributed by layers from 21–30 m for all sites. In general, the soil deposits in Kanpur region have silty sand and sand deposits and are prone to liquefaction hazards due to drastic decrease of cyclic resistance ratio (CRR) at four chosen sites in Kanpur.     

Earthquake ground motion predictive equations for Garhwal Himalaya,India

Predictive equations based on the stochastic approach are developed for earthquake ground motions from Garhwal Himalayan earthquakes of 3.5≤M_w≤6.8 at a distance of 10≤R≤250 km. The predicted ground motion parameters are response spectral values at frequencies from 0.25 to 20 Hz, and peak ground acceleration (PGA). The ground motion prediction equations (GMPEs) are derived from an empirically based stochastic ground motion model. The GMPEs show a fair agreement with the empirically developed ground motion equations from Himalaya as well as the NGA equation. The proposed relations also reasonably predict the observed ground motion of two major Himalayan earthquakes from Garhwal Himalayan region. For high magnitudes, there is insufficient data to satisfactorily judge the relationship; however it reasonably predicts the 1991 Uttarkashi earthquake (M_w=6.8) and 1999 Chamoli earthquake (M_w=6.4) from Garhwal Himalaya region.     

Analysis of observations backing up the existence of VLF and ionospheric TEC anomalies before the Mw6.1 earthquake in Greece,January 26, 2014

The present work integrates ground-based ionosphere measurements using very-low-frequency radio transmissions with satellite measurements of the total electron content to draw common conclusions about the possible impact that the Mw6.1 earthquake that took place in Greece on January 26, 2014, had on the ionosphere.Very-low-frequency radio signals reveal the existence of an ∼4-day anomaly in the wavelet spectra of the signals received inside the earthquake preparation zone and a significant increase in the normalized variance of the signals prior to the earthquake (approximately 1 day before).Through total electron content analysis, it was possible to identify a clear anomaly from 15:00 until 20:00 UT on the day before the earthquake that appears again on the day of the earthquake between 07:00 UT and 08:00 UT. The anomalous values reach TEC1Sigma ∼4.36 and 3.11, respectively. Their spatial and temporal distributions give grounds to assume a possible link with the earthquake preparation. The geomagnetic, solar and weather conditions during the considered period are presented and taken into account.This work is an initial and original step towards a multi-parameter approach to the problem of the possible earthquake-related effects on the ionosphere joining observations made from both ground stations and satellites. A well-founded knowledge of these phenomena is clearly necessary before dealing with their application to earthquake prediction purposes.