Ground motion estimation in the Eastern-Southern Alps: from Ground Motion Predictive Equations to real-time ShakeMaps

The aim of this PhD thesis is to estimate ground motions in the South-Eastern Alps area. For this purpose we purposed empirical relationships that estimate the ground motion parameters as function of the magnitude, the epicentral distance and the soil geological characterization. Later on these relationships are used to calibrate the ShakeMaps software to generate ground motion shake maps in real time for the Friuli-Venezia Giulia region. The GMPEs (Ground Motion Predictive Equations) for PGA, PGV and SA are computed in the South-Eastern Alps area using strong motion observations. 900 accelerometric waveforms are selected and filtered between 0.1 Hz and 30 Hz; the epicentral distance varies from 1 km to 100 km, while the local magnitude, calibrated by comparison with various catalogues, varies in a relatively wide range (3.0 £ ML £ 6.3). Various attenuation models are tested and the best result is selected by the adoption of specific evaluation criteria derived from statistical considerations (R2 value, ANOVA test, residuals analysis). The coefficients of the final model are determined from ML, the epicentral distance, the site effects, the magnitude saturation, the correlation between the distance and the magnitude and the near-source effects. The coefficients of the GMPEs are computed from vertical and horizontal components (the latter represented both as the largest horizontal component and the vectorial addiction); the cluster analysis reduces the final uncertainties on the empirical relations. The comparison with the previous results evidences that the obtained relationships are characterized by a strong attenuation at low magnitudes and large distances. Similar results are obtained for the relationships derived from data recorded all over Northern Italy. The recent evolution of the seismic networks provides a large number of data, available in real time, so it is possible to quickly estimate shake maps. The “ShakeMap” software has been adapted to the South-Eastern Alps region and implemented to obtain a stable interface with the “Antelope” acquisition system in order to extract the ground motion parameters from the waveforms and the generation of the shake maps within 5 minutes from the earthquake occurrence. This procedure requires a dense and uniform spatial distribution of the recording instruments in the field and a geological classification of the soil derived from the average velocities of the S waves in the first 30m below the recording instruments (Vs30). In the geological classification the soil is divided into three classes (bedrock, stiff soil and soft soil), and the amplification coefficients are computed using the relationships proposed by Borcherdt (1994). The related maps are generated using the GMPEs and the empirical relations that predict the macroseismic intensity from the ground motion, both derived from data observed in the Alpine region. The GMPEs that are obtained in this thesis are inserted in the ShakeMap software to generate shake maps in real time or quasi real time in North-Eastern Italy. To evaluate the effects of the station coverage on the shake maps, synthetic seismograms are computed for the Bovec 2004 earthquake by varying the grid size and the network geometry. The results indicate that a dense and uniform spatial distribution in the field and a careful choice of the grid size are crucial to calibrate the shake maps in a given geographical area. The shake maps are generated for eight earthquakes occurred in the studied area in the last 30 years. Furthermore, the finite-fault model is utilized for the seismic events of the Friuli 1976 and Bovec 1998 selecting the source parameters proposed in previous studies. The model validation is done computing the misfit value between the observed macroseismic data (DBMI04 catalogue) and the “instrumental” intensities that are obtained from the synthetic seismograms using empirical relationships between the ground motion and intensity. This analysis has been done for the earthquakes of Cansiglio (1936), Friuli (1976) and Bovec (1998). The synthetic seismograms are calculated for an upper cutoff frequency of 10 Hz applying the reflectivity model. The ground motion parameters are extracted from synthetic signals computed at the presently operating seismic stations and the shake maps are generated. The macroseismic intensity is derived from various relationships; the lowest misfit is obtained using the relation proposed by Kästli and Fäh (2006) for all considered seismic events and this seem to validate our Shake Maps model.

Lo scopo di questa tesi di dottorato è la stima del moto forte del suolo nell’area delle Alpi Sud-Orientali. A tal fine sono state proposte delle relazioni empiriche che stimano i parametri del moto in funzione della magnitudo, della distanza dall’epicentro e della classificazione geologica del suolo; successivamente tali relazioni sono state usate per calibrare il software ShakeMaps con il fine di generare in tempo reale le mappe di scuotimento del terreno per la regione Friuli-Venezia Giulia. Le GMPEs (Ground Motion Predictive Equations) per PGA, PGV e SA sono state calcolate nell’area delle Alpi Sud-Orientali utilizzando registrazioni del moto forte del terreno. Sono state selezionate 900 forme d’onde accelerometriche filtrate tra 0.1 Hz e 30 Hz; la distanza epicentrale varia tra 1 km a 100 km, mentre la magnitudo locale, opportunamente calibrata confrontando diversi cataloghi, varia in un intervallo relativamente ampio (3.0 £ ML £ 6.3). Sono stati testati diversi modelli di attenuazione e il miglior risultato è stato individuato utilizzando specifici criteri di valutazione derivanti da considerazioni di carattere statistico (valore di R2, uso dell’ANOVA test, analisi dei residui). I coefficienti del modello finale sono stati determinati oltre che da ML, dalla distanza epicentrale e dagli effetti dovuti al sito, anche dalla saturazione della magnitudo, dalla correlazione tra magnitudo e distanza e dagli effetti di “near-source”. I coefficienti delle GMPEs sono stati calcolati per le componenti verticali ed orizzontali (rappresentata sia con la componente maggiore sia con la somma vettoriale delle due componenti); la tecnica dell’analisi dei gruppi ha permesso di ridurre l’incertezza finale sulle relazioni empiriche. Il confronto con i risultati ottenuti precedentemente evidenzia come le relazioni ottenute in questa tesi abbiano una maggiore attenuazione a basse magnitudo e a grandi distanze; risultati analoghi sono stati ottenuti per le relazioni ricavate dai dati registrati in tutta l’Italia Settentrionale. L’evoluzione recente delle reti sismiche rende oggi disponibile una grossa mole di dati acquisiti in tempo reale, per cui risulta fattibile stimare velocemente lo scuotimento del terreno tramite mappe; il software “ShakeMap” è stato adattato alle Alpi Sud-Orientali implementato allo scopo di ottenere una stabile interfaccia con il sistema di acquisizione dati “Antelope” che garantisca l’estrazione dei parametri del moto dalle forme d’onda e la creazione delle mappe di scuotimento entro 5 minuti dall’evento sismico. Questa procedura richiede una fitta e uniforme distribuzione spaziale degli strumenti di registrazione sul territorio e una classificazione geologica del suolo fatta usando le velocita’ medie, Vs30, dei primi 30m del mezzo immediatamente sotto gli strumenti. La classificazione geologica del suolo prevede la suddivisione in tre categorie (suolo rigido, suolo addensato e suolo soffice) mentre i coefficienti di amplificazione sono stati calcolati usando le relazioni proposte da Borcherdt (1994). Le relative mappe vanno calcolate usando le GMPEs e le relazioni empiriche che legano il moto del terreno all’intensità macrosismica, basate ambedue su dati registrati nella regione alpina. Le GMPEs discusse in precedenza sono state inserite nel software “ShakeMap” per la produzione delle mappe di scuotimento in tempo reale e quasi-reale nell’Italia Nord-Orientale. Per valutare l’effetto della densità di stazioni sulle mappe di scuotimento sono stati calcolati dei sismogrammi sintetici relativi al terremoto di Bovec 2004 variando il passo di griglia e la geometria dei ricevitori. I risultati ottenuti indicano come una distribuzione fitta e uniforme di strumenti sul territorio e una scelta accurata delle dimensioni della griglia dei ricevitori siano cruciali per calibrare le mappe di scuotimento in una ben determinata area geografica. Le mappe di scuotimento del suolo sono state generate per otto terremoti avvenuti nell’area considerata negli ultimi 30 anni; inoltre per gli eventi del Friuli 1976 e Bovec 1998 è stato utilizzato il modello di faglia finita con i parametri di sorgente stimati in precedenti studi. La validazione del modello è stata fatta calcolando il misfit tra le intensità macrosismiche osservate (catalogo DBMI04) e quelle “strumentali” che sono state ottenute dai sismogrammi sintetici tramite relazioni empiriche tra moto del suolo ed intensità. L’analisi è stata fatta per i terremoti del Cansiglio (1936), del Friuli (1976) e di Bovec (1998). I sismogrammi sintetici sono stati calcolati ad una frequenza massima di 10 Hz applicando il modello della riflettività; i parametri del moto sono stati estratti dai segnali sintetici calcolati nelle attuali stazioni di registrazione e successivamente sono state generate le mappe di scuotimento. L’intensità macrosismica “strumentale” è stata ricavata applicando diverse relazioni; il minor misfit è stato ottenuto usando le relazioni proposte da Kästli and Fäh (2006) per tutti e tre i terremoti considerati, il che sembra validare il nostro modello di Shake Maps.