Deterministic ground-motion scenarios for engineering applications: the case of thessaloniki, Greece.

In this paper we present a deterministic study to estimate seismic ground motions expected in urban areas located near active faults. The purpose was to generate bedrock synthetic time series to be used as seismic input into site effects evaluation analysis and loss estimates for the urban area and infrastructures of Thessaloniki (Northern Greece). Two simulation techniques (a full wave method to generate low frequency,~< 1Hz, and a hybrid deterministic-stochastic technique to simulate high-frequency seismograms, ~> 1 Hz) were used to compute time series associated with four different reference earthquakes having magnitude from 5.9 to 6.5 and located within 30 km of Thessaloniki. The propagation medium and different source parameters were tested through the modeling of the 1978 Thessaloniki earthquake (M 6.5). Moreover two different nucleation points were considered for each fault in order to introduce additional variability in the ground motion estimates. Between the two cases, the quasi-unilateral rupture propagation toward the city produces both higher median PGA and PGV values and higher variability than bilateral ones. Conversely, the low-frequency ground motion (PGD) is slightly influenced by the position of the nucleation point and its variability is related to the final slip distribution on the faults of the reference earthquakes and to the location of the sites with respect to the nodal planes of the radiation pattern. To validate our deterministic shaking scenarios we verified that the synthetic peak ground motions (PGA, PGV) and spectral ordinates are within one standard deviation of several ground-motion prediction equations valid for the region. At specific sites we combined the low- and high-frequency synthetics to obtain broadband time series that cover all the frequency band of engineering interest (0-25 Hz). The use of synthetic seismograms instead of empirical equations in the hazard estimates provides a complete evaluation of the expected ground motions both in frequency and time domains, including predictions at short distances from the fault (0 – 10 km) and at periods larger than 2 – 3 seconds

Modelling directivity effects of the October 21, 2002 (Mw = 5.7), Molise, Southern Italy, earthquake

Acceleration time series recorded by the Italian Strong Motion Network (RAN) during the October 31, 2002 (Mw=5.8), Molise earthquake, are employed in order to investigate source effects on the ground motion in the epicentral area. We consider two different seismogenic sources: a fault model inferred from inversion of teleseismic, regional and local seismic signals [Vallée and Di Luccio, 2005], and a fault model based on seismotectonic data [Basili and Vannoli, 2005]. Both source studies suggest a deep location of the earthquake fault plane (ranging from 6.0 to 20.1 km and from 12.0 to 19.9 km, respectively), however, with considerably different fault lengths (5.2 and 10.5 km, respectively), and widths (14.2 and 8 km, respectively). Due to these differences, only the second model allows for effective horizontal unilateral rupture propagation. Finite fault effects are modelled by the Deterministic-Stochastic-Method (DSM) [Pacor et al., 2005], and the Hybrid Integral-Composite source model (HIC) [Gallovic and Brokesova, 2006]. In both methods k-square slip distributions on the faults are considered. We simulate the October 31, 2002 earthquake considering: 1) Vallée and Di Luccio [2005] faultwith a bilateral rupture propagation, and 2) Basili and Vannoli [2005] fault with unilateral directions of the rupture propagation. The spectral attenuation is modelled using a regional estimate of the quality factor [Castro et al., 2004] and k values estimated from acceleration records. Comparison between synthetic and recorded data at nearby stations (hypocentral distances < 60 km) performed in terms of frequency content and peak ground motion, favours the model with unilateral propagation of the rupture. Assuming the source model with unilateral rupture propagation, we utilize both asymptotic and full wave field methods in order to simulate ground shaking scenarios for an area extending up to 150 km epicentral distance. These results are then subjected to comparison with peak ground accelerations recorded in the far field

Approcci deterministici per la stima del moto del suolo: vantaggi e limiti

Nel Progetto DPC-INGV S3 “Scenari di scuotimento in aree di interesse prioritario e/o strategico”, le stime del moto del suolo sono state ottenute attraverso l’applicazione di diverse tecniche di simulazione di sismogrammi sintetici. Le esperienze maturate nel corso del progetto hanno condotto alla stesura di linee guida per il calcolo degli scenari di scuotimento al bedrock attraverso approcci deterministici. In questo ambito è stata introdotta una classificazione degli scenari deterministici secondo tre diversi livelli di complessità

Variability of kinematic source parameters and its implication on the choice of the design scenario

Near-fault seismic recordings for recent earthquakes (Chi Chi earthquake, 1999, and Parkfield earthquake, 2004) show the high spatial heterogeneity of ground motion. This variability is controlled by fault geometry, rupture complexity, and also by wave propagation and site effects. Nowadays, the number of available records in the near-source region is still not enough to infer a robust parameterization of the ground motion and to retrieve multiparametric predictive equations valid at close distances from the fault. The use of a synthetic approach may help to overcome this limitation and to study the strong ground motion variability. In this article we focus on ground-motion dependence on different earthquakes breaking the same fault, as it has been rarely recorded by instruments. We model seismic scenarios from different rupture models of a fault similar to the 1980 Irpinia, Italy, earthquake source (Mw 6.9). A discrete wavenumber/finite element technique is used to compute fullwave displacement and velocity time series in the low-frequency band (up to 2 Hz). We investigate the variability of the ground motion as a function of different source parameters (rupture velocity, slip distribution, nucleation point, and source time function), whose values depend on the state of knowledge of the physical model driving the process. The probability density functions of the simulated ground-motion parameters, such as displacement response spectrum and peak ground velocity, are used to identify particular scenarios that match specific engineering requests

Variability of kinematic source parameters and its implication on the choice of the design scenario

Near-fault seismic recordings for recent earthquakes (Chi Chi earthquake, 1999; Parkfield earthquake, 2004) show the high spatial heterogeneity of ground motion. This variability is controlled by fault geometry, rupture complexity, and also by wave propagation and site effects. Nowadays, the number of available records in near-source region is still not enough to infer a robust parameterization of the ground motion and to retrieve multi-parametric predictive equations valid at close distances from the fault. The use of a synthetic approach may help to overcome this limitation and to study the strong ground motion variability. In this paper we focus on ground-motion dependence on different earthquakes breaking the same fault, as it has been rarely recorded by instruments. We model seismic scenarios from different rupture models of a fault similar to the 1980 Irpinia, Italy, earthquake source (Mw 6.9). A discrete wavenumber-finite element technique is used to compute full-wave displacement and velocity time series in the low-frequency band (up to 2 Hz). We investigate the variability of the ground motion as a function of different source parameters (rupture velocity, slip distribution, nucleation point, source time function), whose values depend on the state of knowledge of the physical model driving the process. The probability density functions of the simulated ground motion parameters, such as displacement response spectrum (SD) and peak ground velocity (PGV), have been used to identify particular scenarios that match specific engineering requests

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

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

Ground motion scenarios for the 1997 Colfiorito, central Italy, earthquake

In this paper we report the results of several investigations aimed at evaluating ground motion scenarios for the September 26th, 1997 Colfiorito earthquake (Mw 6.0, 09:40 UTC). We model the observed variability of ground motions through synthetic scenarios which simulate an earthquake rupture propagating at constant rupture velocity (2.7 km/s) and the inferred directivity. We discuss the variability of kinematic source parameters, such as the nucleation position and the rupture velocity, and how it influences the predicted ground motions and it does not account for the total standard deviation of the empirical predictive model valid for the region. Finally, we used the results from the scenario studies for the Colfiorito earthquake to integrate the probabilistic and deterministic approaches for seismic hazard assessment

Ground motion shaking scenarios for the 1997 Colfiorito earthquake

In the recent years, two Italian research projects have been devoted to the simulation of ground shaking scenarios in different areas. A large part of the activities has been performed in the Umbria region and was in particular related to the 1997 Colfiorito earthquake. In general the statistical-deterministic approach was adopted for evaluating the scenarios for strong motion parameters (peak values, spectral ordinates, signal integral quantities, and so on) associated with the occurrence of a characteristic earthquake on a given fault. This approach is based on the realistic occurrence of a single earthquake related to the fracture of an a priori well identified active fault. According to the characteristic earthquake model, an earthquake rupture can repeatedly occurs along the same fault (or fault system) with an almost constant geometry, mechanism and seismic moment, these parameters being mainly related to the direction and intensity of the large scale tectonic stress regime. These ideas are supported by numerous paleoseismic studies of active faults in different tectonic environments [e.g., Pantosti and Valensise, 1990]. On the other hand, each faulting process may not repeat the same style of nucleation, propagation and arrest during successive rupture episodes occurring along a given fault zone, depending these characteristics on the pre-fracturing conditions of rock strength and/or yielding stress along the fault zone. It is therefore assumed that the large scale source characteristics (i.e., fault size and position, focal mechanism and seismic moment) are a priori known as the result of previous geological, geophysical and historical seismicity investigations. The variability of the rupture process is expected to produce variable strong ground motions at the earth surface, depending on the distribution of the kinematic parameters (final slip distribution, rupture velocity, slip duration …) along the faulting surface. In order to account for the possible variation of the source process from one rupture event to another, a large number of synthetic seismograms should be computed for different (and possible) rupture histories occurring along the characteristic fault selected, so to provide a representative set of strong motion records to be used for hazard estimation. By this strategy, the massive computation of synthetics for different possible rupture models does not provide a single earthquake scenario (as for the standard deterministic approach) but a set of possible scenarios whose variability substantially reflects the heterogeneity of the source process. The advantage of this approach is that the variability of the selected strong ground motion parameter at a given site can be described by the statistical quantities inferred from the large number of simulations available. The earthquake scenario can then be represented, for example, by a couple of maps, one describing the spatial distribution of the mean value of the considered ground motion parameter and the other representing the associated variability for example in terms of standard deviation