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.     

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.     

Application of horizontal-to-vertical (H/V) Fourier spectral ratio for analysis of site effect on rock (NEHRP-class B) sites in Taiwan

The frequency-dependent amplification for rock (NEHRP-class B) sites was studied using earthquake ground-motion database collected in Taiwan during implementation of the Taiwan Strong Motion Instrumentation Program. The database used includes several hundred records from earthquakes of M_L 4.0–7.3 occurred between 1993 and 2004. The characteristics of amplification were evaluated using the well-known technique of horizontal-to-vertical Fourier spectral ratio (H/V) of the S-wave phase [Lermo J, Chavez-Garcia FJ. Site effect evaluation using spectral ratios with only one station. Bull Seism Soc Am 1993;83:1574–94]. The study allows us to analyze peculiarities of rock sites amplification in Northern and Eastern Taiwan. It was suggested to divide the NEHRP-class B site amplification into four types based on frequency of maximum amplification and the shape of amplification function. The applicability of the technique was also checked for a few stiff and soft soil sites (NEHRP-classes D and E).     

Estimation of Shallow S-Wave Velocity Structure in Damascus City,Syria, Using Microtremor Exploration

Array measurements of microtremors were carried out at thirty sites in Damascus city, Syria to estimate S-wave velocity structures of shallow soil formations for site effect analysis. The microtremor data were recorded by 6 vertical-component seismometers distributed along the circumferences of two circles as well as a 3-component seismometer deployed in the center. The phase velocities were estimated at each site from the vertical components of recorded microtremor data by using the Spatial Autocorrelation method. Then, Genetic Simulated Annealing Algorithm technique was applied for inversion of the phase velocities to estimate 1-D S-wave velocity structures beneath the sites. The inverted Vs profiles are not uniform in Damascus city and the results show that a shallow soft layer (∼200 m/s) appears in the eastern part of the city as well as the central part along Barada River. This layer controls the amplification distribution in the city with a high amplification mainly observed at the locations having this layer. The inversion results also show that the depth to the engineering bedrock (∼750 m/s) is very shallow along the foothills of Mt. Qasyoun in the north-west. Then the depth increases towards the east and the south. The maximum depth to the engineering bedrock (∼80 m) was observed in the southern part of Damascus. To validate the results of the inversions, the spectral ratios between the horizontal and vertical components (H/V) of the recorded microtremor data at the central seismometer were compared with the computed ellipticities of the fundamental-mode Rayleigh-waves based on the respective Vs structure. The results show a good agreement in a period range of 0.05 s to 0.5 s. In this period range, the dominant peaks of the H/V ratios are due to the overall effect of the velocity contrasts between the shallow layers representing the subsurface S-wave velocity structure. Moreover, the average S-wave velocity for the top 10 m of soils (V_S10) shows a better correlation with the averaged site amplification in a period range of 0.05 s to 0.5 s than V_S30 which indicates that V_S10 can be a better proxy for high-frequency site amplification in the case of Damascus city.     

Soil characteristics in Doon Valley (north west Himalaya,India) by inversion of H/V spectral ratios from ambient noise measurements

Past and recent observations have shown that the local site conditions significantly affect the behavior of seismic waves and its potential to cause destructive earthquakes. Thus, seismic microzonation studies have become crucial for seismic hazard assessment, providing local soil characteristics that can help to evaluate the possible seismic effects. Among the different methods used for estimating the soil characteristics, the ones based on ambient noise measurements, such as the H/V technique, become a cheap, non-invasive and successful way for evaluating the soil properties along a studied area.In this work, ambient noise measurements were taken at 240 sites around the Doon Valley, India, in order to characterize the sediment deposits. First, the H/V analysis has been carried out to estimate the resonant frequencies along the valley. Subsequently, some of this H/V results have been inverted, using the neighborhood algorithm and the available geotechnical information, in order to provide an estimation of the S-wave velocity profiles at the studied sites.Using all these information, we have characterized the sedimentary deposits in different areas of the Doon Valley, providing the resonant frequency, the soil thickness, the mean S-wave velocity of the sediments, and the mean S-wave velocity in the uppermost 30 m.     

安徽强震动台站场地响应研究

采用H/V谱比法,利用安徽强震台网记录的波形数据,计算了3个强震台的场地响应。结果表明,当台站场地为基岩或覆盖层较薄时,由地脉动和S波数据分别计算出的谱比曲线一致性很好。当覆盖层大于30.5 m的时候,两者的放大因子有差异,即用S波计算的场地放大响应要明显大于用脉动计算的场地响应。     

Ambient noise H/V spectral ratio in site effects estimation in Fateh jang area,Pakistan

Local geology or local site effect is a crucial component while conducting seismic risk assessment studies. Investigations made by utilization of ambient noise are an effective tool for local site estimation. The present study is conducted to perform site response analysis at 13 different sites within urban settlements of Fateh jang area(Pakistan). The aim of this study was achieved by utilizing Nakamura method or H/V spectral ratio method. Some important local site parameters, e.g., the fundamental frequencies f0 of soft sediments, amplitudes A0 of corresponding H/V spectral ratios, and alluvium thicknesses over 13 sites within the study area, were measured and analyzed. The results show that the study area reflects low fundamental frequency f0. The fundamental frequencies of the sediments are highly variable and lie in a range of 0.6–13.0 Hz. Similarly, amplification factors at these sites are in the range of 2.0–4.0.     

Site characterizations for the Tehran network (TDMMO) in Tehran Region Using Micro-Earthquake,Microtremor and Quarry Blast Data

We did a Study of Horizontal-to-Vertical Component Spectral Ratio in the Tehran seismic zone. Micro-earthquakes, microtremors and quarry blasts data were used as an estimation of the site response in the Tehran zone. Site effects were studied based on horizontal to vertical ratios by the Nakamura׳s technique. Also, we used the spectra of signals for three components with the lowest noise levels for spectral slope studies. The analysis used seismic events from a network of 13 seismic stations by the permanent local seismological network of the Tehran Disaster Mitigation and Management Organization (TDMMO) from 2004 to 2007. The number of events used were different for each station. Quarry blast events were with 1.2≤M_L≤2.2 and micro-earthquakes were with 1.1≤M_L≤4.1.By comparing results for earthquake, microtremor and quarry blast, we could see that there is a significant difference between them. The data showed clear observations, especially in high-frequencies. The H/V spectral ratios indicate dominant frequency for rock/soft site with a higher ratio level for quarry blast ratios, which are comparable to the earthquake results due to their difference sources. The results derived by spectral H/V ratios and spectral analysis may be used to distinguish between local earthquakes and quarry blasts.     

Empirical evaluation of microtremor H/V spectral ratio

The objective of this work is to perform a purely empirical assessment of the actual capabilities of the horizontal-to-vertical (H/V) spectral ratio technique to provide reliable and relevant information concerning site conditions and/or site amplification. This objective has been tackled through the homogeneous (re)processing of a large volume of earthquakes and ambient noise data recorded by different research teams in more than 200 sites located mainly in Europe, but also in the Caribbean and in Tehran. The original recordings were first gathered in a specific database with information on both the sites and recorded events. Then, for all sites close to an instrumented reference, average site-to-reference spectral ratios (“spectral ratio method” (SSR)) were derived in a homogeneous way (window selection, smoothing, signal-to-noise ratio threshold, averaging), as well as H/V ratios (“HVSRE–RF”) on earthquake recordings. H/V ratios were also obtained from noise recordings at each site (either specific measurements, or extracted from pre- or post-event noise windows). The spectral curves resulting from these three techniques were estimated reliable for a subset of 104 sites, and were thus compared in terms of fundamental frequency, amplitude and amplification bandwidth, exhibiting agreements and disagreements, for which interpretations are looked for in relation with characteristics of site conditions. The first important result consists in the very good agreement between fundamental frequencies obtained with either technique, observed for 81% of the analyzed sites. A significant part of the disagreements correspond to thick, low frequency, continental sites where natural noise level is often very low and H/V noise ratios do not exhibit any clear peak. The second important result is the absence of correlation between H/V peak amplitude and the actual site amplification measured on site-to-reference spectral ratios. There are, however, two statistically significant results about the amplitude of the H/V curve: the peak amplitude may be considered as a lower bound estimate of the actual amplification indicated by SSR (it is smaller for 79% of the 104 investigated sites), and, from another point of view, the difference in amplitude exhibits a questioning correlation with the geometrical characteristics of the sediment/basement interface: large SSR/HV differences might thus help to detect the existence of significant 2D or 3D effects.     

The Effect of Velocity Inversions on H/V

We analyzed the phenomenology of microtremor H/V curves under inversions in the shear-wave velocity (V_s) profile in the subsoil. Under no V_s inversion the spectral signature of the H/V peaks is found to be ‘eye-shaped’ with the horizontal components higher than the vertical. Conversely, under negative velocity gradients, numerous of differences emerge. I) A H/V ratio below 1 is observed for a wide range of frequencies, due to the decrease of the horizontal components below the vertical one. II) In the presence of persistent H/V < 1, small bumps in the H/V ratio given by local minima in the vertical spectral component may represent the relics of the peaks indicating resonances and stratigraphic discontinuities. As a consequence, in the presence of velocity inversions the H/V > 2 SESAME (2004) criterion fails but a stratigraphic interpretation may still be possible. III) The H/V curves should always be interpreted together with the single component spectra. IV) Microtremor H/V measurements for stratigraphic/microzonation purposes on stiff artificial soils, (asphalt, concrete, cement, pavements) should always be avoided since the latter often produce velocity inversions. This may have consequences in the intermediate to high frequency domain ( > 1 Hz) also in the application of reference site methods, like H_site/H_bedrock, to microtremor. Theoretical modeling confirms these experimental findings.     

H/V spectral ratio analysis and Rayleigh modelization in Eastern Thuringia,Germany

We use records from the East Thuringian Seismic Network (OTSN, Ostthüringer Seismisches Netzwerk) to characterize the site response for each station, and to analyze the scope and limits of the Rayleigh modeling for H/V spectral ratio. The stations considered in this work can be classified by their seismic response as hard rock sites or as sites with some site effect. From these results we propose velocity models based on Rayleigh modeling (theoretical Rayleigh wave ellipticity).Our results show that for locations affected by site effects the H/V spectral ratio can be modeled by the theoretical ellipticity of layered velocity models. For hard rock sites the spectral ratio is rather flat and the modeling with the theoretical ellipticity was not very clear. This may be explained by the fact that for hard rock sites the conditions for a clear fundamental frequency associated with S-wave resonance, and therefore with Rayleigh wave ellipticity, are not fulfilled.     

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.     

Applicability of the H/V spectral ratio of microtremors in assessing site effects on seismic motion

The authors examine the reliability of site response estimations obtained by the horizontal to vertical (H/V) spectral ratios of microtremors by means of cross‐validation with the ratio of the horizontal spectra of earthquake motion with respect to reference sites. The data comprise microtremor and ground motion records recorded at 150 sites of Yokohama strong motion array. The use of non‐supervised pattern recognition techniques aims to group the sites with more objectivity. Attributes defining the overall shape of the amplification spectra serve as input in the computation of Euclidean distance similarity coefficients amongst sites. The implementation of the Ward clustering scheme leads to the attainment of a meaningful tree diagram. Its analysis shows the possibility of summarizing the results into six general patterns. A good coincidence of site effects estimates at 80 per cent of the sites becomes apparent. However, this coincidence appears poor for sites characterized by H/V amplification ratios around 2 or smaller and predominant periods longer than 0.5 s. In such cases, the presence of stiff, sandy sediments in the soil profile proves common. To proscribe H/V estimations, relying solely on the small spectral ratios criterion seems inadequate. Copyright © 2001 John Wiley & Sons, Ltd.     

Response of the radiation belt electron flux to the solar wind velocity: Parameterization by radial distance and energy

The solar wind velocity is the primary driver of the electron flux variability in Earth's radiation belts. The response of the logarithmic flux (“log-flux”) to this driver has been determined at the geosynchronous orbit and at a fixed energy [Baker, D.N., McPherron, R.L., Cayton, T.E., Klebesadel, R.W., 1990. Linear prediction filter analysis of relativistic electron properties at 6.6 RE. Journal of Geophysical Research 95(A9), 15,133–15,140) and as a function of L shell and fixed energy [Vassiliadis, D., Klimas, A.J., Kanekal, S.G., Baker, D.N., Weigel, R.S., 2002. Long-term average, solar-cycle, and seasonal response of magnetospheric energetic electrons to the solar wind speed. Journal of Geophysical Research 107, doi:10.1029/2001JA000506). In this paper we generalize the response model as a function of particle energy (0.8–6.4 MeV) using POLAR HIST measurements. All three response peaks identified earlier figure prominently in the high-altitude POLAR measurements. The positive response around the geosynchronous orbit is peak P₁ (τ=2±1 d; L=5.8±0.5; E=0.8–6.4 MeV), associated with high-speed, low-density streams and the ULF wave activity they produce. Deeper in the magnetosphere, the response is dominated by a positive peak P₀ (0±1 d; 2.9±0.5R_E; 0.8–1.1 MeV), of a shorter duration and producing lower-energy electrons. The P₀ response occurs during the passage of geoeffective structures containing high IMF and high-density parts, such as ICMEs and other mass ejecta. Finally, the negative peak V₁ (0±0.5 d; 5.7±0.5R_E; 0.8–6.4 MeV) is associated with the “D_st effect” or the quasiadiabatic transport produced by ring-current intensifications. As energies increase, the P₁ and V₁ peaks appear at lower L, while the D_st effect becomes more pronounced in the region L<3. The P₀ effectively disappears for E>1.6 MeV because of low statistics, although it is evident in individual events. The continuity of the response across radial and energy scales supports the earlier hypothesis that each of the three modes corresponds to a qualitatively different type of large-scale electron acceleration and transport.     

Forecasting the ionospheric f0F2 parameter one hour ahead using a support vector machine technique

This paper proposes a method for forecasting the ionospheric critical frequency, f₀F₂, 1 h in advance, using the support vector machine (SVM) approach. The inputs to the SVM network are the time of day, seasonal information, 2 month running mean sunspot number (R2), 3 day running mean of the 3 h planetary magnetic ap index, the solar zenith angle, the present value f₀F₂(t) and its first and second increments, the observation of f₀F₂ at t?23 h, the 30-day mean value at time, t, f_mF₂ (t) and the previous 30 day running mean of f₀F₂ at t?23 h f_mF₂(t?23). The output is the predicted f₀F₂ 1 h ahead. The network is trained to use the ionospheric sounding data at Haikou, Guangzhou, Chongqing, Lanzhou, Beijing, Changchun and Manzhouli stations at high and low solar activities. The performance of the SVM model was verified with observed data. It is shown that the predicted f₀F₂ has good agreement with the observed f₀F₂. The performance of the SVM model is superior to that of the autocorrelation and persistence models, and that it is comparable to that of the neural network model.     

The Himalayan foreland basin crust and upper mantle

Modeling of multimode surface wave group velocity dispersion data sampling the eastern and the western Ganga basins, reveals a three layer crust with an average V_s of 3.7 km s^?1, draped by ～2.5 km foreland sediments. The Moho is at a depth of 43 ± 2 km and 41 ± 2 km beneath the eastern and the western Ganga basins respectively. Crustal V_p/V_s shows a felsic upper and middle crust beneath the eastern Ganga basin (1.70) compared to a more mafic western Ganga basin crust (1.77). Due to higher radiogenic heat production in felsic than mafic rocks, a lateral thermal heterogeneity will be present in the foreland basin crust. This heterogeneity had been previously observed in the north Indian Shield immediately south of the foreland basin and must also continue northward below the Himalaya. The high heat producing felsic crust, underthrust below the Himalayas could be an important cause for melting of midcrustal rocks and emplacement of leucogranites. This is a plausible explanation for abundance of leucogranites in the east-central Himalaya compared to the west. The uppermost mantle V_s is also significantly lower beneath the eastern Ganga basin (4.30 km s^?1) compared to the west (4.44 km s^?1).     

Elasticity of antigorite,seismic detection of serpentinites,and anisotropy in subduction zones

Serpentinization of the mantle wedge is an important process that influences the seismic and mechanical properties in subduction zones. Seismic detection of serpentines relies on the knowledge of elastic properties of serpentinites, which thus far has not been possible in the absence of single-crystal elastic properties of antigorite. The elastic constants of antigorite, the dominant serpentine at high-pressure in subduction zones, were measured using Brillouin spectroscopy under ambient conditions. In addition, antigorite lattice preferred orientations (LPO) were determined using an electron back-scattering diffraction (EBSD) technique. Isotropic aggregate velocities are significantly lower than those of peridotites to allow seismic detection of serpentinites from tomography. The isotropic V_P/V_S ratio is 1.76 in the Voigt–Reuss–Hill average, not very different from that of 1.73 in peridotite, but may vary between 1.70 and 1.86 between the Voigt and Reuss bonds. Antigorite and deformed serpentinites have a very high seismic anisotropy and remarkably low velocities along particular directions. V_P varies between 8.9 km s^? 1 and 5.6 km s^? 1 (46% anisotropy), and 8.3 km s^? 1 and 5.8 km s^? 1 (37%), and V_S between 5.1 km s^? 1 and 2.5 km s^? 1 (66%), and 4.7 km s^? 1 and 2.9 km s^? 1 (50%) for the single-crystal and aggregate, respectively. The V_P/V_S ratio and shear wave splitting also vary with orientation between 1.2 and 3.4, and 1.3 and 2.8 for the single-crystal and aggregate, respectively. Thus deformed serpentinites can present seismic velocities similar to peridotites for wave propagation parallel to the foliation or lower than crustal rocks for wave propagation perpendicular to the foliation. These properties can be used to detect serpentinite, quantify the amount of serpentinization, and to discuss relationships between seismic anisotropy and deformation in the mantle wedge. Regions of high V_P/V_S ratios and extremely low velocities in the mantle wedge of subduction zones (down to about 6 and 3 km.s^?1 for V_P and V_S, respectively) are difficult to explain without strong preferred orientation of serpentine. Local variations of anisotropy may result from kilometer-scale folding of serpentinites. Shear wave splittings up to 1–1.5 s can be explained with moderately thick (10–20 km) serpentinite bodies.     

Optical signatures of dissolved organic matter in the watershed of a globally large river (Yangtze River,China)

Parallel factor analysis of fluorescence excitation emission matrices of surface water samples of a globally large river (Yangtze River, China) watershed identified three classes of fluorescent dissolved organic matter (FDOM) that had ex/em = 280/330 nm, 305/385 nm and 350/450 nm respectively, resembling “peak T”, “peak M” and “peak C” commonly identified in natural water, respectively. Peak T (a tyrosine/tryptophan-like FDOM) did not show correlations to peak M or C which were humic-like substances, while a positive correlation (r = 0.935, p < 0.001) was present between the natural log-transformed maximum fluorescence intensity (F_max) of peaks T and M indicating a tight link during their production and processing. F_max values (in Raman unit nm^?1) normalized to dissolved organic carbon (DOC) concentration were low, varying in ranges 15.93–85.95, 29.83–83.54 and 19.73–51.05 × 10^?5 nm^?1 (μmol/L)^?1 for peaks T, M and C, respectively, in line with the history of strong photobleaching of the water samples as indicated by fairly high absorption spectral slope ratios (0.75–1.53 with a mean 1.03). Intermediate fluorescence index (FI) (1.46–1.83 with a mean 1.61) and small specific absorption at 254 nm (0.64–1.93 with a mean 1.15 m^?1 mg^?1 L) of the water samples, indicated the presence of both aquatic microbial DOM (e.g. peak T) and soil DOM (e.g. peak C). Peak C could be substantially removed by UV-A (320–400 nm) irradiation, while peak M was slightly increased when a microbe-containing water was exposed to the same UV-A irradiation. Taken together, peak C was attributed to diffuse soil source while peak M was likely attributed to joint effects of microbial activities and solar irradiation on the chromophores in the sample.     

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.