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.     

Estimation of Strong Ground Motions for 2001 Bhuj (M w 7.6), India Earthquake

The strong ground motions for the 2001 Bhuj (M _w 7.6) India earthquake have been estimated on hard rock and B/C boundary (NEHRP) levels using a recently modified version of stochastic finite fault modeling based on dynamic corner frequency (Motazedian and Atkinson in Bull Seismol Soc Am 95, 995–1010 2005). Incorporation of dynamic corner frequency removes the limitations of earlier stochastic methods. Simulations were carried out at 13 sites in Gujarat where structural response recorder (SRR) recordings are available. In addition, accelerograms were simulated at the B/C boundary at a large number of points distributed on a grid. The corresponding response spectra have also been estimated. The values of peak ground accelerations and spectral accelerations at three periods (0.4, 0.75 and 1.25 s) are presented in the form of contour maps. The maximum value of peak ground acceleration (PGA) in the center of meizoseismal zone is 550 cm/s². The response spectral acceleration in same zone is 900 cm/s² (T = 0.4 s), 600 cm/s² (T = 0.75 s) and 300 cm/s² (T = 1.25 s). The innermost PGA contour is on the fault plane. A comparison of the PGA values obtained at 13 sites in this study with those obtained in earlier studies on the same sites, but employing different methods, show that the present PGA values are comparable at most of the sites. The rate of decay of PGA values is fast at short distances as compared to that at longer distances. The PGA values obtained here put some constraints on the expected values from a similar earthquake in the region. A synthetic intensity map has been prepared from the estimated values of PGA using an empirical relation. A comparison with the reported intensity map of the earthquake shows the synthetic MMI values, as expected, are lower by 1 unit compared to reported intensity map. The contour map of PGA along with the contour maps of spectral acceleration at various periods permit the assessment of damage potential to various categories of houses and other structures. Such information will be quite important in planning of mitigation and disaster management programs in the region.     

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.     

Probability of actuating a seismic shutoff device

Automatic seismic shutoff devices are used to reduce the risk to gas and liquid distribution systems from earthquakes. In the USA, the gas shutoff devices are tested and certified according to the American Society of Civil Engineers' Standard ASCE 25. During tests, devices are shaken by simple harmonic (sinusoidal) motions of different frequencies and checked for actuation. Because earthquake motions are not sinusoidal, the amplitude of earthquake motions that will actuate these devices is not clearly understood. This paper determines the probability of actuating devices by earthquake motions of different amplitudes. The probability of actuation increases with increase in the resultant peak horizontal ground acceleration (PGA). The probability of actuation is 50% for PGA = 0.23g and 90% for PGA = 0.31g, where g = 9.81 m/s² = acceleration due to gravity. On a ‘stiff soil’ site in San Francisco, CA, the mean recurrence interval of actuation is 51 years. On a similar site in Boston, MA, the mean recurrence interval of actuation is 3000 years. ASCE 25 compliant devices are actuated by high frequencies in ground motions. There is greater uncertainty in the actuation of these devices by ground motions that are damaging to very flexible systems. Copyright © 2016 John Wiley & Sons, Ltd.     

Simulation of strong earthquake characteristics of a scenario earthquake (MS7.5) based on the enlightenment of 2022 MS6.9 earthquake in Menyuan

The Menyuan area is an important transportation hub in the Hexi Corridor. The Menyuan M_S6.9 earthquake that occurred on January 8, 2022 had a major impact on the local infrastructure and transportation of this region. Due to the high possibility of similar strong earthquakes occurring in this area in the future, preliminary assessment of the seismic intensity characteristics of destructive earthquakes in this region is essential for effective disaster control. This paper uses the empirical Green′s function (EGF) method as a numerical simulation tool to predict the ground motion intensity of Datong Autonomous County under the action of the scenario earthquake (M_S7.5). Seismic records of aftershocks of the 2016 Menyuan M_S6.4 earthquake were used as Green’s functions for this simulation. The uncertainties associated with various source parameters were considered, and 36 possible earthquake scenarios were simulated to obtain 72 sets of horizontal ground motions in Datong County. The obtained peak ground acceleration (PGA) vs. time histories of the horizontal ground motion were screened using the attenuation relationships provided by the fifth-edition of China's Seismic Ground Motion Parameter Zoning Map and the NGA-West2 dataset. Ultimately, 32 possible acceleration-time histories were selected for further analysis. The screened PGA values ranged from 78.8 to 153 cm/s². The uncertainty associated with the initial rupture point was found to greatly affect the results of the earthquake simulation. The average acceleration spectrum of the selected acceleration-time history exceeded the expected spectrum of a intermediate earthquake, which means that buildings in Datong County might sustain some damage should the scenario earthquake occur. This research can provide reliable ground motion input for urban earthquake damage simulation and seismic design in Datong County. Growing the dataset of small earthquakes recorded in this region will facilitate the large-scale simulation of ground motions under different earthquake scenarios.     

Seismic fragility of traction elevators

Second‐generation performance‐based earthquake engineering (PBEE‐2) requires a library of component fragility functions to estimate probabilistic damage to a wide variety of building components. The present work draws on a large body of (mostly) post‐earthquake reconnaissance and (some) post‐earthquake survey observations of traction elevators to create fragility functions useful in PBEE‐2. Two surveys provide detailed observations of 115 representative elevators at 12 hospitals shaken in the 1989 Loma Prieta and 1994 Northridge earthquakes and selected without regard to or foreknowledge of damage. Of these, 55 failed and 60 did not. Approximately half were installed after an important code change of 1972, so one can distinguish the performance of pre‐1973 and post‐1973 elevator construction. They experienced a range of strong motion: 22 with peak ground acceleration (PGA) < 0.25 g, 93 with 0.25 g < PGA < 0.85 g. The hospitals had elevator failure rates as low as 0% and as high as 100%. A third survey describes damage qualitatively for six sites with PGA ≤ 0.25 and per‐site failure rates of 0% to perhaps 30%. Fragility functions are offered where the damage state is the loss of functionality of the elevator. The elevators in these surveys exhibit a median capacity of PGA ≈ 0.35 g with a logarithmic standard deviation of 0.40. Capacity is modestly sensitive to whether the elevator was installed before or after 1973. Using building‐specific intensity measures such as S_a(T₁) does not improve the fragility functions. Copyright © 2015 John Wiley & Sons, Ltd.     

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.     

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.     

Peak and integral seismic parameters of L��Aquila 2009 ground motions: observed versus code provision values

The M _w 6.3 L’Aquila earthquake of April 6, 2009 hit a wide area of the Abruzzo region (Central Italy). The epicentre of the main shock was very close to the urban centre of L’Aquila, the regional capital, with an epicentral distance less than 10 km. It was the strongest earthquake ever recorded in Italy which has provided ground motion recordings from accelerometric stations located in close proximity to the epicentre. Because of this, several remarkable results can be achieved by analysing the strong motion recorded signals in terms of peak (PGA, PGV and PGD) and integral (Housner Intensity, I _H) seismic parameters. Additionally, an alternative time-domain representation of recorded signals has been used to furnish a rapid comparison of traces recorded at different stations and along different directions. Some comparisons between the response spectra derived from the recordings and the elastic demand spectra provided in the new seismic Italian code have also been performed. PGA recorded values are very high and generally higher than code values for seismic actions with return period T _R = 475 years. In some cases, this also happens for seismic actions with T _R up to 2,475 years. With regard to I _H, recorded values are generally higher for T _R = 475 years, whilst they are remarkably lower for T _R = 2,475 years. Accurate analyses have been carried out in the article to better understand the above differences and their significance and implications.     

A note on the M w 6.3 earthquake in Iceland on 29 May 2008 at 15:45 UTC

This preliminary study aims to investigate a M _w 6.3 earthquake that occurred in South Iceland on Thursday 29 May 2008 at 15:45 UTC. The epicentre was in the Olfus District between the towns of Selfoss and Hveragerdi. This study examines the data recorded and the damage observed immediately after the event. Horizontal accelerations of up to 80%g were recorded in the epicentral region and there is visual evidence that the vertical acceleration exceeded 1 g. The PGA data is compared to a ground motion estimation model developed for the South Iceland earthquakes in June 2000. In general the basic properties of this event are found to be similar to the characteristics of the South Iceland earthquakes in June 2000. The duration of strong-motion is short and the intensity attenuates rapidly with increasing distance. The earthquake action resisted by buildings in the near fault area is inspected through evaluation of elastic as well as inelastic response spectra. The vast majority of structures seemed to withstand the strong-motion fairly competently and without significant visual damage due firstly to the low-rise, predominantly reinforced concrete or timber, style of buildings. Secondly, the short duration of strong-motion contributed to the endurance of structures.     

Topographic site effects of Xishan Park ridge in Zigong city,Sichuan considering epicentral distance

Ground motion amplification induced by topography plays a vital role in engineering seismology. A topographic array of 8 accelerographs has been operating along the ridge in Xishan Park since 2007. The topographic site effects in Zigong city are studied based on the strong motion data of 2008 M_s 8.0 Wenchuan earthquake (the epicentral distance?=?225 km) and 2019 M_s 5.2 Zizhong earthquake (the epicentral distance?=?29 km). We compare the peak ground acceleration (PGA) of the two earthquakes and find that the PGA of Station 7#, which locates on a relatively steep slope, is amplified by 4.41 times comparing with the reference station in Zizhong earthquake, while this value is only 1.62 in Wenchuan earthquake. Fourier amplitude spectrum shows that the high frequency content of Zizhong earthquake is more abundant because of its smaller epicentral distance. By using the standard spectral ratio (SSR) method, we conclude that the amplification occurs because high-frequency ground motion is likely to resonate at small-scale features. Finally, the 3D numerical simulations are used to verify these conclusions. Our work indicates that more sophisticated numerical models need to be established for more accurate topographic site effects quantification. In addition, the influence of nearby topographic features should be considered when selecting reference stations.

     

Reconnaissance report and preliminary ground motion simulation of the 12 May 2008 Wenchuan earthquake

The M _w = 8.0 Wenchuan earthquake of May 12, 2008, caused destruction over a wide area. The earthquake cost more than 69,000 lives and the damage is reported to have left more than 5 million people homeless. It is estimated that 5.36 million buildings were destroyed and 21 million buildings were damaged in Sichuan and the nearby provinces. Economic losses due to the event are estimated to be 124 billion USD. From a field reconnaissance trip conducted in October 2008, it is evident that the combination of several factors, including mountainous landscape, strong ground shaking, extensive landslides and rock-falls, has exacerbated the human and economic consequences of this earthquake. Extensive damage occurred over a wide area due to the shear size of the earthquake rupture combined with poor quality building construction. In order to investigate the ground shaking during the earthquake, we have conducted a strong ground motion simulation study, applying a hybrid broadband frequency technique. The preliminary results show large spatial variation in the ground shaking, with the strongest ground motions along the fault plane. The simulation results have been calibrated against the recorded ground motion from several near-field stations in the area, and acceleration values of the order of 1 g are obtained, similar to what was recorded during the event. Comparison with the damage distribution observed in the field confirms that the effect of fault rupture complexity on the resulting ground motion distribution also controls to a large extent the damage distribution. The applied simulation technique provides a promising platform for predictive studies.     

Development of damage probability matrices based on Greek earthquake damage data

A comprehensive study is presented for empirical seismic vulnerability assessment of typical structural types, representative of the building stock of Southern Europe, based on a large set of damage statistics. The observational database was obtained from post-earthquake surveys carried out in the area struck by the September 7, 1999 Athens earthquake. After analysis of the collected observational data, a unified damage database has been created which comprises 180,945 damaged buildings from/after the near-field area of the earthquake. The damaged buildings are classified in specific structural types, according to the materials, seismic codes and construction techniques in Southern Europe. The seismic demand is described in terms of both the regional macroseismic intensity and the ratio αg/ao, where αg is the maximum peak ground acceleration (PGA) of the earthquake event and ao is the unique value PGA that characterizes each municipality shown on the Greek hazard map. The relative and cumulative frequencies of the different damage states for each structural type and each intensity level are computed in terms of damage ratio. Damage probability matrices (DPMs) and vulnerability curves are obtained for specific structural types. A comparison analysis is fulfilled between the produced and the existing vulnerability models.     

Quantitative archaeoseismological investigation of the Great Theatre of Larissa, Greece

Larissa, the capital of Thessaly, is located in the eastern part of Central Greece, at the southern border of a Late Quarternary graben, the Tyrnavos Basin. Palaeoseismological, morphotectonic, and geophysical investigations as well as historical and instrumental records show evidences for seismic activity in this area. Previous investigations documented the occurrence of several moderate to strong earthquakes during Holocene time on active faults with recurrence intervals of a few thousand years. The historical and instrumental records suggest a period of seismic quiescence during the last 400–500 years. The present archaeoseismological research, based on a multidisciplinary approach is devoted to improve the knowledge on past earthquakes, which occurred in the area. This study focuses on damages observed on the walls of the scene building of the Great Theatre of Larissa. The Theatre was built at the beginning of the third century BC and consists of a semicircular auditorium, an almost circular arena and a main scene building. Archaeological and historical investigations document a partial destruction of the theatre during the second to first century BC. Recent excavations show that the building complex after it was repaired suffered additional structural damages, probably from seismic loading. The damages investigated in detail are displacements, rotations and ruptures of numerous blocks at the walls of the scene building. In order to test the earthquake hypothesis as cause of the damages a simplified seismotectonic model of the Tyrnavos Basin and its surroundings was used with a composite earthquake source model to calculate synthetic seismograms at the Larissa site for various earthquake scenarios. Horizontal to vertical seismic ratio (HVSR) measurements in the theatre and its vicinity were used to estimate local site effects. The synthetic seismograms are then used as input accelerations for a finite element model of the walls, which simulates seismically induced in-plane sliding within the walls. Results show that some of the surrounding faults have the potential to produce seismic ground motion that can induce in-plane sliding of blocks. Model calculations were used to constrain the characteristics of the ground acceleration and infer the causative fault and earthquake by comparing the calculated and observed distribution of the displacements of the blocks. Ground motions with a PGA at the site of 0.62–0.82 g, which could be induced by an M 5.8–6.0 earthquake on the nearby Larissa Fault, would be sufficient to explain the damage.     

Estimation of Strong Ground Motion from a Great Earthquake Mw 8.5 in Central Seismic Gap Region, Himalaya (India) Using Empirical Green’s Function Technique

In the present study ground motions for a Mw 8.5 scenario earthquake are estimated at 13 sites in Kumaun-Garhwal region using the empirical Green’s function technique. The recordings of 1991 Uttarkashi earthquake of Mw 6.8 at these sites are used as an element earthquake. A heterogeneous source model consisting of two asperities is considered for simulating the ground motions. The entire central seismic gap (CSG) can expect acceleration in excess of 100 cm/s² with NW portion in excess of 400 cm/s² and SE between 100 and 200 cm/s². The central portion can expect peak ground acceleration (PGA) between 200 and 400 cm/s². It has been observed from simulation of strong ground motion that sites located near the rupture initiation point can expect accelerations in excess of 1g. In the present analysis, Bhatwari and Uttarkashi can expect ground accelerations in excess of 1g. The estimates of the PGA are compared with earlier studies in the same region using different methodologies and it was found that the results are comparable. This has put constrains on the expected PGAs in this region. The obtained PGA values can be used in identifying the vulnerable areas in the central Himalaya, thereby facilitating the planning, design and construction of new structures and strengthening of the existing structures in the region.     

Ground-motion Attenuation Relation from Strong-motion Records of the 2001 Mw 7.7 Bhuj Earthquake Sequence (2001–2006), Gujarat, India

Predictive relations are developed for peak ground acceleration (PGA) from the engineering seismoscope (SRR) records of the 2001 M_w 7.7 Bhuj earthquake and 239 strong-motion records of 32 significant aftershocks of 3.1 ≤ M_w ≤ 5.6 at epicentral distances of 1 ≤ R ≤ 288 km. We have taken advantage of the recent increase in strong-motion data at close distances to derive new attenuation relation for peak horizontal acceleration in the Kachchh seismic zone, Gujarat. This new analysis uses the Joyner-Boore’s method for a magnitude-independent shape, based on geometrical spreading and anelastic attenuation, for the attenuation curve. The resulting attenuation equation is,

Strong motion simulation by the composite source modeling： A case study of 1679 M8.0 Sanhe-Pinggu earthquake

In this study, a composite source model has been used to calculate the realistic strong ground motions in Beijing area, caused by 1679 MS8.0 earthquake in Sanhe-Pinggu. The results could provide us the useful physical parame-ters for the future seismic hazard analysis in this area. Considering the regional geological/geophysical background, we simulated the scenario earthquake with an associated ground motions in the area ranging from 39.3°N to 41.1°N in latitude and from 115.35°E to 117.55°E in longitude. Some of the key factors which could influence the characteristics of strong ground motion have been discussed, and the resultant peak ground acceleration (PGA) distribution and the peak ground velocity (PGV) distribution around Beijing area also have been made as well. A comparison of the simulated result with the results derived from the attenuation relation has been made, and a suf-ficient discussion about the advantages and disadvantages of composite source model also has been given in this study. The numerical results, such as the PGA, PGV, peak ground displacement (PGD), and the three-component time-histories developed for Beijing area, have a potential application in earthquake engineering field and building code design, especially for the evaluation of critical constructions, government decision making and the seismic hazard assessment by financial/insurance companies.     

System identification estimation of soil properties at the Lotung site

Dynamic properties of the soils at the Lotung test site, Lotung, Taiwan, are estimated from seismic vertical array measurements (input–output data sets) using both time-invariant and time-variant parametric modeling methods (system identification). Soil properties are directly mapped from model parameters to an equivalent lumped mass model of the soil interval. Shear stiffness and damping ratios were calculated for 8 events with M_L ranging from 4.5 to 7.0. Shear stiffness ranged between 0.5 and 6 MN/m, inversely proportional to PGA. The equivalent viscous damping ratio varied from 2 to 30% of critical damping, proportional to PGA. Degradation of soil behavior, while less pronounced with increasing depth, consistently occurs above a peak input acceleration of 0.07 g. Although “non-linear” behavior is evident above 0.17 g, Event 7 (0.21 g) is accurately predicted using a linear constant parameter model estimated from the smaller Event 8 aftershock ground motions.     

6th April 2009 L��Aquila earthquake, Italy: reinforced concrete building performance

On 6th April 2009 an earthquake of magnitude M _w =  6.3 occurred in the Abruzzo region; the epicentre was very close to the city of L’Aquila (about 6 km away). The event produced casualties and damage to buildings, lifelines and other infrastructures. An analysis of the main damage that reinforced concrete (RC) structures showed after the event is presented in this study. In order to isolate the main causes of structural and non-structural damage, the seismological characteristics of the event are examined, followed by an analysis of the existing RC building stock in the area. The latter issue came under scrutiny after the release of official data about structural types and times of construction, combined with a detailed review of the most important seismic codes in force in the last 100 years in Italy. Comparison of the current design provisions of the Italian and European codes with previous standards allows the main weaknesses of the existing building stock to be determined. Damage to structural and non-structural elements is finally analyzed thanks to photographic material collected in the first week after the event; the main causes of damage are then inferred.     

The role of alternating outcrops of sediments and basaltic lavas on seismic urban scenario: the study case of Catania, Italy

Experimental data and numerical modelling were used to study the effect of local geology on the seismic response of the Catania area. The town extends on a marly clays bedrock and terraced deposits made up by coastal sands and alluvial conglomerates. This sedimentary substratum is deeply entrenched by paleo-valleys filled by lava flows and pyroclastics. Available borehole data and elastic parameters were used to reconstruct a geotechnical model in order to perfome 1D numerical modeling. Seismic urban scenarios were simulated considering destructive (M _w = 7.0), strong (M _w = 6.2) and moderate (M _w = 5.7) earthquakes to assess the shaking level of the different outcropping formations. For each scenario seven real accelerograms were selected from the European Strong Motion Database to assess the expected seismic input at the bedrock. PGA and spectral acceleration at different periods were obtained in the urban area through the equivalent linear numerical code EERA, and contour maps of different levels of shaking were drawn. Standard and horizontal-to-vertical spectral ratios were achieved making use of a dataset of 172 seismic events recorded at ten sites located on the main outcropping lithotypes. Spectral ratios inferred from earthquake data were compared with theoretical transfer functions. Both experimental and numerical results confirm the role of the geological and morphologic setting of Catania. Amplification of seismic motion mainly occurs in three different stratigraphic conditions: (a) sedimentary deposits mainly diffused in the south of the study area; (b) spots of soft sediments surrounded by lava flows; (c) intensely fractured and scoriaceous basaltic lavas.