The Mw 6.3, 2009 L’Aquila earthquake: source, path and site effects from spectral analysis of strong motion data

The strong motion data of 2009 April 6 L’Aquila (Central Italy) earthquake (Mw = 6.3) and of 12 aftershocks (4.1 ≤ Mw ≤ 5.6) recorded by 56 stations of the Italian strong motion network are spectrally analysed to estimate the source parameters, the seismic attenuation, and the site amplification effects. The obtained source spectra for S wave have stress drop values ranging from 2.4 to 16.8 MPa, being the stress drop of the main shock equal to 9.2 MPa. The spectral curves describing the attenuation with distance show the presence of shoulders and bumps, mainly around 50 and 150 km, as consequence of significant reflected and refracted arrivals from crustal interfaces. The attenuation in the first 50 km is well described by a quality factor equal to Q( f ) = 59 f 0.56 obtained by fixing the geometrical spreading exponent to 1. Finally, the horizontal-to-vertical spectral ratio provides unreliable estimates of local site effects for those stations showing large amplifications over the vertical component of motion

Frequency variation in site response as observed from strong motion data of the L’Aquila (2009) seismic sequence

Previous works based mainly on strong-motion recordings of large Japanese earthquakes showed that site amplification and soil fundamental frequency could vary over long and short time scales. These phenomena were attributed to non-linear soil behaviour: the starting fundamental frequency and amplification were both instantaneously decreasing and then recovering for a time varying from few seconds to several months. The recent April 6, 2009 earthquake (M W 6.3), occurred in the L’Aquila district (central Italy), gave us the possibility to test hypotheses on time variation of amplification function and soil fundamental frequency, thanks to the recordings provided by a pre-existing strong-motion array and by a large number of temporary stations. We investigated the intra- and inter-event soil frequency variations through different spectral analyses, including time-frequency spectral ratios and S-Transform (Stockwell et al. in IEEE Trans Signal Process 44:998–1001, 1996). Finally, analyses on noise recordings were performed, in order to study the soil behaviour in linear conditions. The results provided puzzling evidences. Concerning the long time scale, little variation was observed at the permanent stations of the Aterno Valley array. As for the short time-scale variation, the evidence was often contrasting, with some station showing a time-varying behavior, while others did not change their frequency with respect to the one evaluated from noise measurements. Even when a time-varying fundamental frequency was observed, it was difficult to attribute it to a classical, softening non-linear behaviour. Even for the strongest recorded shocks, with peak ground acceleration reaching 0.7 g, variations in frequency and amplitude seems not relevant from building design standpoint. The only exception seems to be the site named AQV, where the analyses evidence a fundamental frequency of the soil shifting from 3 Hz to about 1.5 Hz during the mainshock

Comparison between empirical predictive equations calibrated at global and national scale and the Italian strong-motion data

In Italy in the last years many ground motion prediction equations (hereinafter GMPEs) were calibrated both at national and regional scale using weak and strong motion data recorded in the last 30 years by several networks. Moreover many of the Italian strongest earthquakes were included in global datasets in order to calibrate GMPEs suitable to predict ground-motion at very large scale. In the last decade the Sabetta and Pugliese (1996) relationships represented a reference for the ground motion predictions in Italy. At present all Italian strong-motion data, recorded from 1972 by RAN (Italian Accelerometric Network), and more recently by other regional networks (e.g. RAIS, Strong motion network of Northern Italy), are collected in ITACA (ITalian ACcelerometric Archive). Considering Italian strong-motion data with Mw  4.0 and distance (Joyner-Boore or epicentral) up to 100 km, new GMPEs were developed by Bindi et al. (2009), aimed at replacing the older Italian relationships. The occurrence of the recent 23rd December 2008, Mw 5.4, Parma (Northern Italy) earthquake and the 6th April 2009, Mw 6.3, L’Aquila earthquake, allowed to upgrade the ITACA data set and gave us the possibility to validate the predictive capability of many GMPEs, developed using Italian, European and global data sets. The results are presented in terms of quality of performance (fit between recorded and predicted values) using the maximum likelihood approach as explained in Spudich et al. (1999). Considering the strong-motion data recorded during the L’Aquila sequence the considered GMPEs, in average, overestimate the observed data, showing a dependence of the residuals with distance in particular at higher frequencies. An improvement of fit is obtained comparing all Italian strong-motion data included in ITACA with the European GMPEs calibrated by Akkar and Bommer (2007 a,b) and the global models calibrated by Cauzzi and Faccioli (2008). In contrast, Italian data seem to attenuate faster than the NGA models calibrated by Boore and Atkinson (2008), in particular at higher frequencies

Frequency variation in site response over long and short time scales, as observed from strong motion data of the L’Aquila (2009) seismic sequence

Previous works based mainly on strong-motion recordings of large Japanese earthquakes showed that site amplification and soil fundamental frequency could vary over long and short time scales. These phenomena were attributed to non-linear soil behaviour due to inelastic, softening non-linearity: the starting fundamental frequency and amplification were both decreasing and not recovering for a time varying from few hours to several months. The recent April 6th 2009 earthquake (MW 6.3), occurred in the L'Aquila district (central Italy), gave us the possibility to test hypotheses on time variation of amplification function and soil fundamental frequency, thanks to the recordings provided by a preexisting strong-motion array and by a large number of temporary stations. We performed spectral ratio studies for the permanent stations of the Aterno Valley array where a reference station was available. The temporary stations and permanent ones were studied using time-frequency analyses through the S-Transform approach (Stockwell et al., 1996). Finally, analyses on noise recordings were performed, in order to study the soil behaviour in linear conditions. The results provided puzzling evidences. Concerning the long time scale, little variation was observed at the permanent stations of the Aterno Valley array. As for the short time-scale variation, the evidence was often contrasting, with some station showing a time-varying behavior, while others did not change their frequency with respect to the one evaluated from noise HVSR. Even when a time-varying fundamental frequency was observed, it was difficult to attribute it to a classical, softening non-linear behaviour. Even for the strongest recorded shocks, with PGA reaching 0.7 g, variations in frequency and amplitude seems not relevant from building design standpoint. The only exception seems to be the site named AQV, where the analyses evidence a fundamental frequency of the soil, shifting from 3 Hz to about 1.5 Hz during the mainshock

Ground Motion Prediction Equations Derived from the Italian Strong Motion Database

We present a set of ground motion prediction equations (GMPEs) derived for the geometrical mean of the horizontal components and the vertical, considering the latest release of the strong motion database for Italy. The regressions are performed over the magnitude range 4–6.9 and considering distances up to 200 km. The equations are derived for peak ground acceleration (PGA), peak ground velocity (PGV) and 5%-damped spectral acceleration at periods between 0.04 and 2 s. The total standard deviation (sigma) varies between 0.34 and 0.38 log10 unit, confirming the large variability of ground shaking parameters when regional data sets containing small to moderate magnitude events (M < 6) are used. The between-stations variability provides the largest values for periods shorter than 0.2 s while, for longer periods, the between-events and between-stations distributions of error provide similar contribution to the total variabilit

Overview of the Italian strong Motion database ITACA 1.0

The Italian Strong Motion Database, ITACA, was developed within projects 2 S6 and S4, funded in the framework of the agreements between the Italian Department of 3 Civil Protection (Dipartimento della Protezione Civile, DPC) and the Istituto Nazionale di 4 Geofisica e Vulcanologia (INGV), starting from 2005. The alpha version of the database 5 was released in 2007 and subsequently upgraded to version 1.0 after: (i) including the most 6 recent strongmotion data (from2005 to 2007) recorded in Italy, in addition to the 2008 Parma 7 earthquake, M 5.4, and the M 4.0, 2009 Abruzzo seismic events; (ii) processing the raw 8 strong motion data using an updated procedure; (iii) increasing the number of stations with a 9 measured shear wave velocity profile; (iv) improving the utilities to retrieve time series and 10 ground motion parameters; (v) implementing a tool for selecting time series in agreement 11 with design-response spectra; (vi) compiling detailed station reports containing miscella12 neous information such as photo, maps and site parameters; (vii) developing procedures for 13 the automatic generation of station reports and for the updating of the header files. After such 14 improvements, ITACA 1.0 was published at the web site http://itaca.mi.ingv.it, in 2010. It 15 presently contains 3,955 three-component waveforms, comprising the most complete cata16 logue of the Italian accelerometric records in the period 1972–2007 (3,562 records) and the 17 strongest events in the period 2008–2009. Records were mainly acquired by DPC through its 18 Accelerometric National Network (RAN) and, in few cases, by local networks and temporary 19 stations or networks. This paper introduces the published version of the Italian StrongMotion 20 database (ITACA version 1.0) together with main improvements and new functionalities

Pan-European ground-motion prediction equations for the average horizontal component of PGA, PGV, and 5 %-damped PSA at spectral periods up to 3.0 s using the RESORCE dataset

This article presents a set of Ground-Motion Prediction Equations (GMPEs) for Europe and the Middle East, derived from the RESORCE strong motion data bank, following a standard regression approach. The parametric GMPEs are derived for the peak ground acceleration, peak ground velocity, and 5 %-damped pseudo-absolute acceleration response spectra computed over 23 periods between 0.02 and 3 s, considering the average horizontal-component ground-motions. The GMPEs are valid for distances less than 300 km, hypocentral depth up to 35 km and over the magnitude range 4–7.6. Two metrics for the source-to-station distance (i.e. Joyner-Boore and hypocentral) are considered. The selected dataset is composed by 2,126 recordings (at a period of 0.1 s) related to 365 earthquakes, that includes strong-motion data from 697 stations.The EC8 soil classification (four classes from A to D) discriminates recording sites and four classes (normal, reverse, strike-slip, and unspecified) describe the style of faulting. A subset which contains only stations with measured Vs30 and earthquakes with specified focal mechanism (1,224 records from 345 stations and 255 earthquakes) is used to test of the accuracy of the median prediction and the variability associated to the broader data set. A random effect regression scheme is applied and bootstrap analyses are performed to estimate the 95 % confidence levels for the parameters. The total standard deviation sigma is decomposed into between-events and within-event components, and the site-to-site component is evaluated as well. The results show that the largest contribution to the total sigma is coming from the within-event component. When analyzing the residual distributions, no significant trends are observed that can be ascribed to the earthquake type (mainshock-aftershock classification) or to the non-linear site effects. The proposed GMPEs have lower median values than global models at short periods and large distances, while are consistent with global models at long periods (T>1) s. Consistency is found with two regional models developed for Turkey and Italy, as the considered dataset is dominated by waveforms recorded in these regions

Building damage scenarios based on exploitation of Housner intensity derived from finite faults ground motion simulations

In this paper earthquake damage scenarios for residential buildings (about 4200 units) in Potenza (Southern Italy) have been estimated adopting a novel probabilistic approach that involves complex source models, site effects, building vulnerability assessment and damage estimation through Damage Probability Matrices. Several causative faults of single seismic events, with magnitude up to 7, are known to be close to the town. A seismic hazard approach based on finite faults ground motion simulation techniques has been used to identify the sources producing the maximum expected ground motion at Potenza and to generate a set of ground motion time histories to be adopted for building damage scenarios. Additionally, site effects, evaluated in a previouswork through amplification factors of Housner intensity, have been combined with the bedrock values provided by hazard assessment. Furthermore, a new relationship between Housner and EMS-98 macroseismic intensity has been developed. This relationship has been used to convert the probability mass functions of Housner intensity obtained from synthetic seismograms amplified by the site effects coefficients into probability mass function of EMS-98 intensity. Finally, the Damage Probability Matrices have been applied to estimate the damage levels of the residential buildings located in the urban area of Potenza. The proposed methodology returns the full probabilistic distribution of expected damage, thus avoiding average damage index or uncertainties expressed in term of dispersion indexes

La difesa dai terremoti in Lombardia: stato dell’arte e prospettive

L’INGV (http://www.mi.ingv.it/) svolge in Lombardia ricerche nel campo della mitigazione del rischio sismico, mediante studi mirati al miglioramento delle conoscenze sulla storia sismica, sul modello strutturale legato al regime tettonico in atto e alla definizione del moto del suolo atteso. Accanto alle indagini necessarie alla caratterizzazione sismogenetica e dei possibili effetti dello scuotimento sismico, si pone l’attività di monitoraggio sismico, che la Sezione di Milano-Pavia espleta mediante la rete accelerometrica (RAIS, http://rais.mi.ingv.it/) che consta di 20 postazioni distribuite in prevalenza sul territorio regionale. In un recente lavoro di sintesi - di studi pluriennali su fonti storiche e di revisioni critiche di materiali pubblicati - si è tentato di colmare l’evidente carenza conoscitiva sulle caratteristiche della sismicità storica dell’area lombarda. Tale carenza è particolarmente evidente se si rapportano le informazioni oggi disponibili sui terremoti del passato in quest’area con quelle relative al settore veneto-friulano. La regione analizzata, compresa tra il bacino del fiume Adda e il Lago di Garda, è stata caratterizzata da alcuni terremoti con Mw>5.5 (es., 1117, Veronese; 1222, Brescia; 1901, Salò) e vari eventi con Mw compresa fra 4.8 e 5.5 (es., 1065, Brescia; 1396, Monza; 1642, Bergamo). Per molti degli eventi sismici fino al 1700 le informazioni sono desumibili soltanto da scarse fonti storiche. La determinazione epicentrale e l’attribuzione della magnitudo per tali eventi sono da considerarsi, pertanto, con una certa cautela al fine di definire le caratteristiche sismiche del territorio. I terremoti del 1117 e del 1222, in tale contesto, rappresentano un’eccezione rispetto ai dati generalmente disponibili. E tuttavia vari problemi tuttora aperti rendono assai difficile l’utilizzo della distribuzione del danno attribuibile a questi eventi nella prospettiva di un’affidabile parametrizzazione. Le indagini geologico-strutturali e di geologia del Quaternario, finalizzate a definire un quadro strutturale compatibile con il regime tettonico in atto, sono necessarie per giustificare la storia sismica e, in sostanza, per definire il comportamento sismogenetico della regione. Le geometrie dei sistemi di faglia attivi alpini sono ora sufficientemente noti. Le conoscenze permettono di formulare ipotesi sismotettoniche relative all’origine dei terremoti dell’area gardesana e del Bresciano. In via di definizione sono invece le geometrie dei fronti appenninici, cui sono attribuibili i terremoti al di sopra della soglia nel settore padano. Le ricerche attuali sono altresì indirizzate ad una migliore caratterizzazione del complesso settore compreso tra la parte meridionale del Lago di Garda (area di Sirmione), Verona e Mantova, all’interno del quale potrebbe collocarsi l’area epicentrale del terremoto del 1117 (o una delle aree epicentrali, qualora si considerasse questo evento come rappresentato da una sequenza sismica). Uno dei settori regionali con maggiore frequenza di eventi sismici è l’area gardesana occidentale, per questo motivo la rete di monitoraggio presenta una notevole densità di postazioni nel Bresciano e se ne è pianificato l’addensamento nel Veronese. Nel corso del 2007, la rete ha consentito la registrazione di 516 forme d’onda relative a 28 eventi locali e regionali (di cui una decina localizzati nell’area citata) con magnitudo da 1.3 a 4.2, di cui sono stati calcolati i parametri di interesse ingegneristico. L’analisi delle registrazioni ha permesso di ricavare informazioni utili per il calcolo di scenari di scuotimento. Un esempio di taali applicazioni è rappresentato dallo scenario realizzato utilizzando come terremoto di riferimento l'evento del 24 Novembre 2004 (M 5.2)

Building damage scenarios based on exploitation of Housner Intensity derived from finite faults ground motion simulations

In this paper earthquake damage scenarios for residential buildings (about 4200 units) in Potenza (Southern Italy) have been estimate adopting a probabilistic approach that involves complex source models, site effects, building vulnerability assessment and damage estimation through Damage Probability Matrices (DPMs). The studied area experienced several destructive earthquakes in historical and recent times. Several causative faults of single seismic events, with magnitude up to 7, are known to be close to the town. A seismic hazard approach based on finite faults ground motion simulation techniques has been used to identify the sources producing the maximum expected ground motion at Potenza and to generate a set of ground motion time histories to be used for building damage scenarios. Additionally, site effects, evaluated in the framework of the DPC-INGV S3 project through amplification factors of Housner intensity (IH), have been combined with the bedrock values provided by hazard assessment. Furthermore, a new relationship between IH and macroseismic intensity in terms of EMS98 has been developed. This relationship has been used to convert the Probability Density Functions (PDFs) for IH obtained from synthetic seismograms and convolved by the site effects coefficients into PDFs for EMS98 intensity. Finally, the DPMs approach has been applied to estimate the damage levels of the residential buildings in the urban area of Potenza