Layout 6 ANNALS OF GEOPHYSICS, 57, 1, 2014, S0189; doi:10.4401/ag-6363 S0189 The seismic microzonation of level 3 of Sant’Agata Fossili (northern Italy) based on a multidisciplinary approach Giuseppe Di Capua1,*, Massimo Compagnoni2, Giuseppe Di Giulio3, Marco Marchetti1, Giuliano Milana1, Silvia Peppoloni1, Floriana Pergalani2, Vincenzo Sapia1 1 Istituto Nazionale di Geofisica e Vulcanologia, Rome, Italy 2 Politecnico di Milano, Dipartimento di Ingegneria Civile e Ambientale, Milan, Italy 3 Istituto Nazionale di Geofisica e Vulcanologia, L’Aquila, Italy ABSTRACT In this paper the results of a detailed seismic microzonation, per- formed at Sant’Agata Fossili (Piemonte region, northern Italy) are presented. We study the local seismic response of this small village using a level 3, that is the most accurate level following the Italian code of seismic microzonation. The activity steps consist in a gradual widening of knowledge of the different aspects of the amplification phenomena. A multidisciplinary approach has been performed to ob- tain the local seismic response: including a study of local geology, geo- physical and geotechnical characterization of the lithologies, and numerical and experimental analyses. We finally compare the ob- tained elastic response spectra to the prescribed spectra of the Italian Building Code (in Italian: Norme Tecniche per le Costruzioni). Our re- sults show the geologic and geophysical differences of the subsoil, that produce different local seismic response in terms of amplification fac- tors and acceleration response spectra. 1. Introduction The seismic hazard is usually connected to a high building vulnerability and to significant historic values of the urban centres, especially in Italy. In the activities finalized to the assessment of the local seismic hazard and the seismic risk, the geological and geophysical studies can be considered as fundamental to perform ac- tivities of prevention and reconstruction. The aim of a study of seismic microzonation is the evaluation of the expected local seismic hazard, related to the expected seismic input and to geologic and geo- morphologic properties of the subsoil. On the basis of these considerations the Piemonte region (northern Italy) has performed, in a Strategic European Project RISKNAT (http://webgis.arpa.piemonte.it/risknat/), a study of seismic microzonation for the small munici- pality of Sant’Agata Fossili (approximately 500 citizens). This municipality was affected by the low magnitude seismic event of April 11, 2003 (Mw=4.9), and has suf- fered the most damage in the epicentral area, reaching a macroseismic intensity of IMCS=VI-VII. Starting from the recent “Indirizzi e Criteri per la Microzonazione Sismica” (Guidelines and Criteria for the Seismic Microzonation) (ICMS) [GdL MS 2008], a study of level 3 has been carried out. The level 3 con- sists of an accurate analysis of seismic microzonation based on geophysical surveys and numerical simula- tion. A multidisciplinary group has been created and the fundamental steps of the project can be summa- rized as follows: (1) study of local geology; (2) geophysical and geotechnical characterization of the lithologies; (3) numerical and experimental analyses to evalu- ate the local seismic response. In our study, the geological and geophysical het- erogeneities of the subsoil properties have been de- fined. These heterogeneities can be considered as responsible for differences in the local seismic re- sponse, in terms of amplification factors and accelera- tion response spectra. Furthermore, the resulting acceleration response spectra have been compared to those proposed by the Italian seismic code “Norme Tecniche per le Costruzioni” (Technical Standards for Construction) [NTC 2008] in order to have a correct Article history Received May 30, 2013; accepted November 21, 2013. Subject classification: Seismic microzoning, Amplification factors, Elastic acceleration response spectra, Seismic refraction, Electrical resistivity tomography, HVSR, MASW, Down-Hole, Italian Building Code. level of seismic protection in the project design phases. Finally, our study gives some results in term of comparison among different geophysical techniques (noise measurements, Down-Hole tests - DH, Multi- channel Analysis of Surface Waves Method-MASW, 2D array, electrical tomography). 2. Historical and recent seismicity The information on the historical seismicity for Sant’Agata Fossili is scarce. This could be due to the small dimension of the municipality and therefore its little political and economical importance, that can lead to poor consideration in the historical bibliography. An alternative explanation could be the low local seismic- ity and consequently the absence of damages; proba- bly all the two causes are present. Considering the seismic events in the “Catalogo Parametrico dei Terremoti Italiani” (Parametric Cata- logue of Italian Earthquakes: GdL CPTI [2004] and Rovida et al. [2011]), the area has been hit by low en- ergy earthquakes. Particularly the strongest and nearest events are: 1541 (IMCSmax=8), 1680 (IMCSmax=7), 1828 (IMCSmax=8) and 1913 (IMCSmax=5) earthquakes, with magnitude between 4.9 and 5.8, producing minor dam- age in the municipality (Figure 1). The recent seismicity shows events characterized by magnitude lesser than 4 and only the event of April 11, 2003 has reached magnitude producing a damage equal to IMCS=VI-VII. 3. Experimental setting To investigate the subsoil, 19 pre-existent bore- holes (with thicknesses between 15 and 20 m) have been collected, 2 new boreholes (S1 and S2, with thick- nesses between 40 and 30 m) and 2 Standard Penetra- tion Test (SPT) have been performed, with the collec- tion of 6 undisturbed samples to carry out the dynamic and static laboratory tests. Different geophysical inves- tigations, as 2 Down-Hole surveys (DH), 27 noise meas- urements, 1 passive 2D array, 1 MASW, 1 seismic refraction line, 2 electrical tomography, have been con- ducted (Figure 2). The investigations have pointed out a homogeneous lithological sequence, characterized by different thickness of the deposits. The sequence is generally formed by an- thropic filling material or eluvial-colluvial deposits ap- proximately 2 m thick over-imposed to silts, sands and clays. Below these deposits, marls can be considered as a geologic bedrock. In the northern-western area, gravels and cemented gravels have been found below the an- thropic filling deposits or silts. In the central area an arti- ficially stabilized landslide has been recognized. On the 6 samples, collected in the silts sands and clays, the following geotechnical analyses have been performed: classification tests, oedometric tests, shear and resonant column tests. A series of seismic noise measurements have been carried out (Figure 2), finalized to obtain information on the resonant frequencies (f0). DI CAPUA ET AL. 2 Figure 1. Location of the strongest and closest earthquakes to San- t’Agata Fossili. The orange points show the most important towns in the area closest to Sant’Agata Fossili, the squares indicate the epi- centres of historical and recent earthquakes, the stars show the epi- centre of the April 11, 2003 earthquake. Figure 2. Location of the existing and new geotechnical and geo- physical investigations. 3 The average spectral ratios H/V have been calculated [HVNSR; Nakamura 1989], and the curves are repre- sented in Figure 3. The results show clear resonance peaks, characterized by a f0 values varying between 2.1 Hz and 7.4 Hz in the northern and western sector of Sant’Agata Fossili, whereas, in the eastern sector, the HVNSR curves show the absence of clear resonance peaks, indicating likely the absence of a strong Vs con- trast in the subsoil [Ibs-von Seht and Wohlenberg 1999, Delgado et al. 2000a,b, Parolai et al. 2002, Bonnefoy- Claudet et al. 2006, Gosar and Lenart 2010, Benjumea et al. 2011]. In these two identified sectors, two DH (S1 and S2) have been performed and the results in term of Vp and Vs velocity profiles are reported in Figure 4 (location in Figure 2) together with the borehole stratigraphy. The Vs velocity profiles show some differences: S1 is char- acterized by a gradually increment of the Vs with depth, and the seismic bedrock (Vs > 800 m/s) is reached at a depth of about 40 m (eastern sector), whereas S2 shows an high increment of Vs at the contact between the silts and marls (depth 11 m), a second velocity contrast is present in the lower part of the profile at a depth of ap- proximately 30 m (western sector). In the area (blue polygon in Figure 2) between the sectors characterized by the different behavior of Vs profile, some integrative investigations have been per- formed, finalized to the study of the limit between the two sectors; particularly a 2D small-aperture array recording seismic noise and a 1D array of 72 geophones aimed to perform a MASW analysis [Park et al. 2005], using both active (minigun-like) and passive source to derive surface-wave dispersion curves. The HVNSR curves of the 11 stations of the 2D array indicate a f0 varying spatially, in a short distance, from 4 Hz to 6 Hz. The inversion of surface-wave dispersion curves, measured with the 2D and with the 1D arrays, allows to derive a Vs profile. However, the surface wave analysis THE SEISMIC MICROZONATION OF SANT'AGATA FOSSILI Figure 3. Location of the seismic noise measurements. The H/V curves and the resonant frequencies are reported. In the gray polygon, in the northern area, the passive measures in array configuration have been performed. indicate a heterogeneity of the Vs between the western and eastern sectors of the array area (Figure 5), con- firming also the variability of the HVNSR curves. Sum- marizing, it is possible to notice from the two Vs profiles of Figure 5 a first thin layer with a Vs of about 100 m/s; then a second layer characterized by different Vs and thickness values in the two sectors: western sector of the array area with a thickness of about 10 m and Vs of about 340 m/s, eastern sector with a thickness of about 15 m and Vs of about 480 m/s and a third layer charac- terized by a Vs greater than 800 m/s. In the studied area 2 electrical tomography lines (Figure 6, location in Figure 2) have been conducted: the line X-X’ shows the presence of 3 intervals of resis- tivity values: 5-15 Ohm*m, 15-25 Ohm*m and 25-45 Ohm*m and the line Y-Y’, performed on the stabilized landslide, shows a large resistivity contrast between sur- face layers (50-80 Ohm*m) and the bottom layers (40 Ohm*m) and indicates that the landslide thickness is about 8-10 m. Summarizing, on the basis of borehole stratigra- phies and results of geophysical measurements, we have identified 4 sectors for the studied area (Figure 7). Sector 1 in the eastern side of the village shows nearly flat HVNSR curves. This agrees with the results of the DH-S1, which show increasing shear velocity with depth. The Vs profile does not indicate a strong velocity contrast between the anthropic filling material and silt deposits and the underlying clayey marl. The lack of a high contrast can be noticed, considering the results of the electrical tomography (thickness of 10-12 m). In sector 2, the noise measurements indicated a clear resonance peak with f0 varying from 3.6 to 4.9 Hz. The DH-S2 in this sector shows a thickness of 12 m for a low-velocity deposit (Vs < 200 m/s) overlay- ing a stiffer basement and the electrical tomography shows that this deposit can be correlated to a low re- sistivity layer (5-45 Ohm*m). Sector 3 is characterized by the presence of con- glomerates outcropping in the northwest portion of DI CAPUA ET AL. 4 Figure 4. Vp and Vs velocity profiles obtained by Down-Holes in the S1 and S2 boreholes. The different colours identify the man-made de- posits and the eluvial-colluvial deposits, the silts and clays as well as the marls. Figure 5. Inversion of surface-wave dispersion curves, measured with the 2D and with the 1D arrays. 5 the village with a resonance frequency between 2.1 and 3.6 Hz. Sector 4 is in the central part of the village where a stabilized landslide is present. In sector 4 the two noise measurements show f0 of 7.4 and 5.2 Hz, considering the Vs values of 200 m/s as characteristic of the landslide ma- terial, these values allow to define the thickness of the landslide to about 8 m. The landslide thickness is consis- tent with the electrical tomography data (Figure 6). In conclusion, the comparison of the results shows a good agreement among the different investigations, al- lowing a valid reconstruction of the subsoil. 4. Geolithological models and sections The elaboration of all data (geologic, geophysi- cal and geotechnical) give, as mentioned, 4 main geo- lithological sectors (Figure 7) and different lithological models with the individuation of the geo- physical characteristics of each layer. We derive sev- eral geo-lithological sections (Figure 8), used in the 1D numerical analyses of the local seismic response. The layers consist in man-made and eluvial-colluvial de- posits (1-2 m) over-imposed to silts and clays and locally gravels, while the bedrock is formed by marls and, in one area (sector 3), conglomerates. The silt and clay deposits show different charac- teristics in the 4 sectors: (1) In sector 1 the silt and clay deposits have Vs of 250 m/s (Figure 4) and a unit weight of 20 kN/m3. In this sector the H/V curves are characterized by a flat behaviour (Figure 3); (2) In sectors 2, 3 and 4 the Vs value of the deposits is 200 m/s (Figure 4) and the unit weight is 19 kN/m3. The contrast of the Vs between the deposits and the bedrock is significant, as shown by the clear resonance peaks in the H/V curves in the range 4-5 Hz (Figure 3). The marls present a different behaviour of the Vs, as deduced by the results of the DH, considering the THE SEISMIC MICROZONATION OF SANT'AGATA FOSSILI Figure 6. Results of the electric tomography for the 2 analyzed sections. Location in Figure 2. In the section Y-Y’, it is traced the possible sliding surface of the landslide and is reported the stratigraphy of the closest borehole. Figure 7. Areas of the municipality where groups of geophysical investigation, using different techniques, have been performed. two sectors 1 and 2 (Figure 7): (1) In sector 1 the S1 survey shows that marls are characterized by a gradual increase of the Vs values; the unit weight is about 20-21 kN/m3. The bedrock depth is reached at about 35 m; (2) In sector 2 and 4 the marls are compact with a thickness of 17 m and Vs equal to 700 m/s and a unit weight of 21 kN/m3 in the first layer, then the marls reach the Vs values of 800 m/s. In sector 3 the sequence is characterized by silts and clays with a thickness of 3 m, over-imposed to ce- mented gravels and conglomerates characterized by Vs values of 600 m/s and unit weights about 21 kN/m3, over-imposed to bedrock (conglomerate). The dynamic geotechnical tests have given the curves of the behaviour of the normalized shear mod- ulus (G/G0) and the damping ratio (D) relating to the shear stain (γ%) for the silt and clay deposits (Figure 9). The curves relating to the sample of depth 4.0-4.5 m are applied to the silts and clays of the sector 1, the curves relating to the sample of depth 9.0-9.5 m are ap- plied to the silts and clays of the sectors 2, 3, 4 and to marls. The curves for gravels are taken from the scien- tific literature [Rollins 1998]. In the sector 4, characterized by a stabilized land- slide, the subsoil model takes into account the litho- logical and geophysical characteristics discussed in the previous paragraph. DI CAPUA ET AL. 6 Figure 8. Geolithological sections (location in Figure 2): note the man-made and the eluvial-colluvial deposits, the silts and clays, and the geologic bedrock formed by marls. Figure 9. Behaviour of the G/G0 and D correlating to γ. 7 5. Numerical analysis The numerical and experimental analyses are fi- nalized to the evaluation of the expected seismic re- sponse in term of elastic acceleration response spectra and amplification factors, in order to derive a seismic microzonation map of level 3 [GdL MS 2008]. The numerical analyses, as prescribed by level 3 of microzoning, require realistic accelerograms and ac- celeration response spectra as seismic inputs. Based on the seismic hazard of Italy [GdL MPS 2004] the ex- pected maximum acceleration values (amax) and the val- ues of the spectral ordinates of the acceleration response spectra are calculated for different return pe- riods within the national territory. The return period chosen is 475 years, considered as conventional for the local seismic response analysis. The values of amax for the Piemonte region are between 0.036 g and 0.150 g; for Sant’Agata Fossili the closer 4 points of the grid to the coordinates (in latitude and longitude) of the mu- nicipality have been selected: the obtained amax is the average of the 4 points with the inverse of the distance between the municipality point and the grid points and it is 0.096 g. [GdL MPS 2004]. Recorded accelerograms have been extracted from the Database ITACA [Luzi and Sabetta 2006]: following the rules prescribed by the NTC [2008]. In particular for the selection of the seismic input we consider the: (1) historical seismicity and seismogenetic charac- teristic (strike fault) [DISS Working Group 2010]; (2) couple magnitude-distance (from the data in GdL MPS 2004, the characteristic event has a magnitude between 4.5 and 5.0 and a distance between 0-10 km); (3) maximum expected acceleration (0.096 g). Only the accelerograms recorded on bedrock have been selected (subsoil category “A”; NTC [2008]). The selected accelerograms have been scaled in order to have their acceleration peak similar to the ex- pected amax. The characteristics of the 5 accelerograms are reported in Table 1, and are namely: the code of the station, the latitude and longitude of the station, the dis- tance of the station from the epicentre, the name of the event, the name of the station, the component of the motion, the lithology of the station, and the peak ground acceleration (PGA). The accelerograms are plot- ted in Figure 10. The geologic and geomorphologic analyses of the area show the presence of horizontal layers and the presence of a ridge, so we use a mono-dimensional (1D) code [Idriss et al. 1992] for the study of the litho- logical amplifications and a two-dimensional (2D) code [Callerio et al. 2000] for the study of the topographic amplifications. In the 1D code the soil profile is idealized as a sys- tem of homogeneous, visco-elastic sub-layers of infi- nite horizontal extent. The response of this system is calculated considering vertically propagating shear waves. The bedrock is considered deformable, to avoid the reflection of waves into the model, in fact a rigid layer reflects all the reflected waves from the surface, instead in the case of the deformable layer, the waves are spread into the bedrock. The code adopts the equiv- alent linear analysis using an iterative procedure to ob- tain the characteristics of the soil compatible with the effective strain in each layer in each iteration. There- fore, the process is iterative and the code works in the frequency domain, using the Fourier analysis. On the basis of the total amount of the geologi- cal, geotechnical and geophysical data and of the 5 geo-lithological sections (Figure 8), 18 1D geolitho- logical columns have been pointed out (Figure 11). The local seismic response has been computed for each column, using the 1D code and applying the 5 input accelerograms, G/G0 and D curves. The results are expressed as the average of the 5 analyses in terms of elastic acceleration response spectra (using a 5% of the critical damping) and in terms of amplification fac- tors Fa0.1-0.5 and Fa0.5-1.5. These parameters were cal- culated considering the ratio of the elastic velocity THE SEISMIC MICROZONATION OF SANT'AGATA FOSSILI Code Lat. (°) Long. (°) Epicentral distance (km) Earthquake Station Comp. Lithology PGA (g) PNR 44.877 7.344 16.0 VAL CHISONE Pinerolo N-S Rock 0.044 NZZ 44.782 8.357 6.0 ZONA ALESSANDRIA Nizza Monferrato W-E Rock 0.090 TRT 44.892 8.882 17.0 OLTREPO PAVESE Tortona N-S W-E Rock Rock 0.042 0.086 Figure 8. Geolithological sections (location in Figure 2): note the man-made and the eluvial-colluvial deposits, the silts and clays, and the geologic bedrock formed by marls. response spectra between the integral of output and input in the period ranges between 0.1-0.5 s and 0.5- 1.5 s [Pergalani et al. 1999]. The obtained average acceleration response spec- tra have been compared with the response spectra pre- scribed by the national code [NTC 2008], assuming the correspondent category of subsoil for each analysed column (using the Vs30 value). To perform this comparison, the calculated spectra have been adapted to the behaviour of the code spectra (using DI CAPUA ET AL. 8 Figure 10. Accelerograms used in the numerical analyses of the local seismic response. Figure 11. Location of the 18 lithostratigraphic columns used for the analyses of local seismic response 1D. 9 the procedure in GdL MS [2008]) and characterized by a segment of constant spectral acceleration, a segment of constant spectral velocity and a segment of con- stant spectral displacement. In Figure 12 the graphs are reported. Some analysed columns show the same re- sponse spectra, so they have been joined into one graph. As shown, some national code spectra are lower than the 1D calculated spectra, particularly for low periods. The results in terms of average amplification factors are reported in Table 2. The Fa0.1-0.5 values are included be- tween the range of 1.5-2.2 and the Fa0.5-1.5 values are included between the range of 1.1-1.6, typical values of a medium-high amplification. The 2D analyses were performed using a bound- ary elements method. The method considers, with lin- ear segments, only the boundary of the real structure, reducing the computational time, using the Green function and a linear elastic analysis. The structure is characterized by an elastic and homogeneous material. The code considers a continuous displacement field and works in the frequency domain applying Fourier’s trans- formation to the motion equation. The topographic am- plifications due to the presence of a ridge have been cal- culated by the application of a 2D code. The profile of the ridge is characterized by a width (L) of 1700 m and height (H) of 170 m. This morphology can be classified as a pointed ridge (ratio H/L equal to 0.1). In the nu- merical analysis a homogeneous, elastic bedrock char- acterized by a Vs of 800 m/s and a unit weight of 22 kN/m3 has been considered. The average response spec- tra and the average amplification factors do not show morphological effects. 6. Conclusions The level of seismic hazard of Sant’Agata Fossili was computed for homogeneous sectors of the munici- pality, based on geological and geophysical investiga- tions. The values of the amplification factors can be used in planning, for identifying the areas in which the miti- gation actions are necessary and the areas suitable for new constructions. The acceleration response spectra can be used during the project design for new buildings THE SEISMIC MICROZONATION OF SANT'AGATA FOSSILI Figure 12. Acceleration response spectra of the analyzed columns (location in Figure 11) and relative national code spectra. and for the evaluation of the safety of existent buildings. On the basis of the numerical and experimental analyses of the local seismic response, the seismic mi- crozonation map of level 3 has been performed, in which the municipality has been subdivided in 4 main sectors (Figure 13). In particular, for each sector the different response spectra are reported and the different sectors show dif- ferent behaviour. In general, the results of the numerical analyses are higher than the relative national code spec- tra. Only in the sector 3, in presence of the cemented gravels, the amplification phenomena are negligible. The values of the average amplification factors (Fa0.1-0.5 and Fa0.5-1.5), associated to each sector, are re- ported in Table 3; also in this case, the variability of the amplification factors is evident. Only in sector 3, the values of the amplification factors are equal to 1. The HVNSR results show clear resonance peaks, characterized by f0 values between 2.1 Hz and 7.4 Hz in the northern and western sector of Sant’Agata Fossili, whereas, in the eastern sector, the HVNSR curves show the absence of clear resonance peaks. In any case, the interaction between soil resonance and building reso- nance effects needs further investigation, considering the height of reinforced-concrete or masonry buildings and of 2 prevailing storeys, that have likely frequencies > 7 Hz [NTC 2008]. In conclusion, an integrated approach based on geological and geophysical data allows investigating the heterogeneities in the subsoil properties, also consider- DI CAPUA ET AL. 10 PA1 PA2 PA3 PA4 PB2 PD1 PB3 PC3 PC2 PD2 PE1 PE2 DH1 PC1 PE3 PE4 DH2 PB1 Fa0.1-0.5 2.12 1.04 1.98 2.19 2.06 1.55 1.46 1.94 2.15 1.75 2.10 Fa0.5-1.5 1.26 1.01 1.55 1.19 1.40 1.12 1.10 1.13 1.27 1.16 1.36 Table 2. Average amplification factors of the analyzed columns, considering the ratio of the elastic velocity response spectra between the integral of output and input in the period ranges between 0.1-0.5 s and 0.5-1.5 s. Figure 13. Seismic microzonation map of level 3 of Sant’Agata Fossili. For each sector, similar response spectra as presented in Figure 12 have been jointed in one graph. 11 ing a small spatial scale of investigation. This study shows that, in urban contexts, the seismic microzona- tion analyses are valid and that they offer an instrument of knowledge of the territory, used for the seismic pre- vention of buildings, if the seismic vulnerability is known. In fact a realistic evaluation of the expected am- plifications is necessary to perform the reduction of the seismic risk and the safety of the buildings. Acknowledgements. This project has been developed by the financial contribution of Regione Piemonte and by the technical contribution of Mr. Riccardo Conte. References Benjumea B., A. Macau, A. Gabàs, F. Bellmunt, S. Figueras and J. Cirés (2011). Integrated geophysical profiles and H/V microtremor measurements for subsoil characterization. Near Surface Geophysics, 9. Bonnefoy-Claudet S., F. Cotton and P.Y. Bard (2006). The nature of noise wavefield and its applications for site effects studies. A literature review. Earth-Sci- ence Reviews, 79, 205-227. Callerio A., V. Petrini and F. Pergalani (2000). ELCO, A program for two dimensional analyses using bound- ary element method, Technical report, IRRS, Milano. Delgado J., C. Lòpez Casado, A. Estévez, J. Giner, A. Cuenca and S. Molina (2000a). Mapping soft soils in the Segura river valley (SE Spain): a case studi of mi- crotremors as an exploration tool. Journal of Ap- plied Geophysics (45), 19-32. Delgado J., C. Lòpez Casado, J. Giner, A. Estévez, A. Cuenca and S. Molina (2000b). Microtremors as a geophysical exploration tool: application and limita- tion. Pure and applied geophysics (157), 1445-1462. DISS Working Group (2010). Database of Individual Seismogenic Sources (DISS), Version 3.1.1: A com- pilation of potential sources for earthquakes larger than M 5.5 in Italy and surrounding areas, http://diss.rm.ingv.it/diss/. GdL CPTI (2004). Catalogo Parametrico dei Terremoti Italiani, versione 2004 (CPTI04). INGV, Bologna, http://emidius.mi.ingv.it/CPTI04/index.html. GdL MPS (2004). Redazione della mappa di pericolosità sismica prevista dall’Ordinanza PCM 3274 del 20 marzo 2003. Rapporto Conclusivo per il Diparti- mento della Protezione Civile, INGV, Milano-Roma, aprile 2004, 65 pp.+5 appendici, http://zonesis- miche.mi.ingv.it/ (in Italian). GdL MS (2008). Indirizzi e Criteri per la Microzon- azione Sismica. Conferenza delle Regioni e delle Province autonome-Dipartimento della Pro- tezione Civile, Roma, 3 vol. e DVD (in Italian). Gosar A. and A.Lenart (2010). Mapping the thickness of sediments in the Ljubljana Moor basin (Slovenia) using microtremors. Bulletin of earthquake enge- neering (8), 501-518. Ibs-von Seht M., and J. Wohlenberg (1999). Microtremor Measurements Used to Map Thickness of Soft Sedi- ments. Bulletin of the Seismological Society of America , 89 (1). Idriss I.M. and J.I. Sun (1992). User’s manual for SHAKE91, A computer program for conducting equivalent linear seismic response analyses of hori- zontally layered soil deposits, Report of Dip. of Civil and Environmental Engineering, University of Cal- ifornia, Davis. Luzi L.and F. Sabetta (2006). Data base dei dati ac- celerometrici italiani relative al periodo 1972-2004, Istituto Nazionale di Geofisica e Vulcanologia 2006. Nakamura Y. (1989). A method for dynamic character- istics estimation of subsurface using microtremor on the ground surface: Quarterly Report of the Railway Technical Research Institute of Japan, 30, no. 1, 25-33. NTC (2008). Norme Tecniche per le Costruzioni, DM 14 gennaio 2008, Gazzetta Ufficiale n. 29 del 4 feb- braio 2008, Supplemento Ordinario n. 30. Istituto Poligrafico e Zecca dello Stato, Roma (www.cslp.it) (in Italian). Park C.B., R.D. Miller, N. Ryden, J. Xia and J. Ivanov (2005). Combined use of active and passive surface waves. Journal of Engineering and Environmental Geophysics (JEEG), 10, (3), 323-334. Parolai S., P. Bormann and C. Milkereit (2002). New Re- lationships between Vs , Thickness of sediments, and resonance frequency calculated by the H/V ratio of seismic noise for the Cologne area (Ger- many). Bulletin of the seismological society of THE SEISMIC MICROZONATION OF SANT'AGATA FOSSILI Sector 1 Sector 2 Sector 3 Sector 4 Fa0.1-0.5 1.7 2.1 2.1–1.0 2.2 Fa0.5-1.5 1.2 1.4 1.3–1.0 1.2 Table 3. Average amplification factors for each geolithological sector. America, 92 (6), 2521-2527. Pergalani F., R. Romeo, L. Luzi, V. Petrini, A. Pugliese and T. Sanò (1999). Seismic microzoning of the area struck by Umbria-Marche (central Italy) Ms 5.9 earthquake of the 26 september 1997. Soil Dynam- ics and Earthquake Engineering, Elsevier Science, vol. 18, 4, 279-296. Regione Piemonte, Progetto RISKNAT, http://www.risknat-alcotra.org/rna/index.cfm/b- 6.html (in Italian). Rollins K.M. (1998). Shear modulus and damping rela- tionship for gravels. Journal of Geotechnical and Geoenvironmental Engineering. Rovida A., R. Camassi, P. Gasperini and M. Stucchi (eds.) (2011). CPTI11, the 2011 version of the Para- metric Catalogue of Italian Earthquakes. Milano, Bologna, http://emidius.mi.ingv.it/CPTI, doi: 10.6092/INGV.IT-CPTI11. *Corresponding author: Giuseppe Di Capua, Istituto Nazionale di Geofisica e Vulcanologia, Rome, Italy; email: giuseppe.dicapua@ingv.it. © 2014 by the Istituto Nazionale di Geofisica e Vulcanologia. All rights reserved. DI CAPUA ET AL. 12