The seismic microzonation of the city of Catania (Italy) for the maximum expected scenario earthquake of January 11, 1693

Safety against earthquake hazards presents two aspects: structural safety against potentially destructive dynamic forces and site safety related to geotechnical phenomena, such as amplification, landsliding and soil liquefaction. The correct evaluation of seismic hazard is, therefore, highly affected by risk factors due to geological nature and geotechnical properties of soils. In response to these new developments, several attempts have been made to identify and appraise geotechnical hazards and to represent them in the form of zoning maps, in which locations or zones with different levels of hazard potential are identified. The geotechnical zonation of the subsoil of the city of Catania (Italy) suggests a high vulnerability of the physical environment added to site amplification of the ground motion phenomena. The ground response analysis at the surface, in terms of time history and response spectra, has been obtained by some 1D equivalent linear models and by a 2D linear model, using a design scenario earthquake as input at the conventional bedrock. In particular, the study has regarded the evaluation of site effects in correspondence of the database of about 1200 boreholes and water-wells available in the data-bank of the Catania area. According to the response spectra obtained through the application of the 1D and 2D models, the city of Catania has been divided into some zones with different peak ground acceleration at the surface, to which corresponds a different value of the Seismic Geotechnical Hazard. A seismic microzoning map of the urban area of the city of Catania has been obtained. The map represents an important tool for the seismic improvement of the buildings, indispensable for the mitigation of the seismic risk.     

2-D spectral element simulations of destructive ground shaking in Catania (Italy)

This study wants to estimate the strong ground motion in the municipal area of Catania (Italy) for a catastrophic earthquake scenario. It is part of a larger research program funded by the National Research Council – National Group for the Defence Against Earthquakes (CNR-GNDT), The Catania Project, devoted to evaluating the seismic risk of a highly urbanised area, such as that of Catania, located in a seismically active region. The reference earthquake simulates the catastrophic event (M 7.2) of 1693. The ground shaking is computed solving the 2-D full-wave equation by the Chebyshev spectral element method (SPEM). Particular emphasis is given to the construction of realistic structural models, also including the finest local detail, obtained from the geophysical, geological and geotechnical data available. Simulations are performed for several sources, to account for both a change in source position and orientation, and the finite extension of the fault along its dip. Synthetic seismograms and peak ground acceleration (PGA) envelopes, calculated at the surface for four transects across the Catania area, constitute the main result of this study which can be used for practical purposes. Simulations show that ground motion is strongly influenced by both source characteristics and crustal structure. We have found that PGA values range between 0.1 g and 0.5 g, although particular site conditions strongly affect these values locally. For example, the frequencies of maximum interest in civil engineering (1.5–4 Hz) are enhanced selectively by a thick portion of surface sediments (i.e., 30–100 m for an average shear wave velocity of 500–600 m/s). An unexpected feature is the appreciable increase of PGA at large epicentral distances, which contradicts classical attenuation relations. All the results are examined through an analysis of the propagating wavefield.     

Empirical prediction of the ground shaking scenario for the Catania area

Simplified approaches for prediction of strong ground shaking, that can be efficiently implemented through a Geographical Information System (GIS), are applied to constructing earthquake scenarios for the Catania municipal area. Ground shaking maps in terms of peak ground acceleration and response spectral ordinates are developed both for a repetition of the catastrophic 1693 earthquake (M7.3), and for a lower energy event, like that which occurred in 1818 (M6.2). The maps account for amplification due to local site conditions. According to the commonly used attenuation relations, the local amplification is taken into account either in a simpler way or in a more detailed approach, derived from an extensive geotechnical zonation. The approximation introduced by the simpler methodology, to be used whenever accurate data are not available, is estimated by a comparison with damage distribution.     

Geostatistical assessment for the regional zonation of seismic site effects in a coastal urban area using a GIS framework

Earthquake-induced hazards are profoundly affected by site effects related to the amplification of ground motions, which are strongly influenced by local geologic conditions such as soil thickness, bedrock depth, and soil stiffness. Seismic disasters are often more severe over soft soils than over stiff soils or rocks due to differences in local site effects. In this study, on the basis of a geotechnical information system (GTIS) framework, we developed an advanced geostatistical assessment for the regional zonation of seismic site effects. In particular, to reliably predict spatial geotechnical information, we developed a procedural methodology for building an advanced GTIS within a geographic information system framework and applied it to the Busan area in Korea. The systemized GTIS comprised four functional components: database, geostatistical analysis, geotechnical analysis, and visualization. First, to build the GTIS, we collected pre-existing geotechnical data in and around the study area, and then conducted a walk-over site survey to acquire surface geo-knowledge data. Second, we determined the optimum geostatistical estimation method using a cross-validation-based verification test, considering site conditions. The advanced GTIS was used in a practical application to estimate the site effects in the study area. We created seismic zoning maps of geotechnical earthquake parameters, such as the depth to bedrock and the site period, and present them as part of a regional synthetic strategy for earthquake risk assessment.     

The role of alternating outcrops of sediments and basaltic lavas on seismic urban scenario: the study case of Catania, Italy

Experimental data and numerical modelling were used to study the effect of local geology on the seismic response of the Catania area. The town extends on a marly clays bedrock and terraced deposits made up by coastal sands and alluvial conglomerates. This sedimentary substratum is deeply entrenched by paleo-valleys filled by lava flows and pyroclastics. Available borehole data and elastic parameters were used to reconstruct a geotechnical model in order to perfome 1D numerical modeling. Seismic urban scenarios were simulated considering destructive (M _w = 7.0), strong (M _w = 6.2) and moderate (M _w = 5.7) earthquakes to assess the shaking level of the different outcropping formations. For each scenario seven real accelerograms were selected from the European Strong Motion Database to assess the expected seismic input at the bedrock. PGA and spectral acceleration at different periods were obtained in the urban area through the equivalent linear numerical code EERA, and contour maps of different levels of shaking were drawn. Standard and horizontal-to-vertical spectral ratios were achieved making use of a dataset of 172 seismic events recorded at ten sites located on the main outcropping lithotypes. Spectral ratios inferred from earthquake data were compared with theoretical transfer functions. Both experimental and numerical results confirm the role of the geological and morphologic setting of Catania. Amplification of seismic motion mainly occurs in three different stratigraphic conditions: (a) sedimentary deposits mainly diffused in the south of the study area; (b) spots of soft sediments surrounded by lava flows; (c) intensely fractured and scoriaceous basaltic lavas.     

Spatial zonations for regional assessment of seismic site effects in the Seoul metropolitan area

Earthquake-induced hazards are profoundly affected by site effects related to the amplification of ground motions, which are strongly influenced by local geologic conditions such as soil thickness or bedrock depth and soil stiffness. In this study, an integrated geographic information system (GIS)-based system for geotechnical data, called the geotechnical information system (GTIS), was developed to establish a regional counterplan against earthquake ground motions in the Seoul metropolitan area. In particular, to reliably predict spatial geotechnical information, a procedural methodology for building the GTIS within a GIS framework was developed and applied to the Seoul area in Korea. To build the GTIS, pre-existing geotechnical data were collected in and around the study area, and then a walk-over site survey was conducted to acquire surface geo-knowledge data. In addition, the representative shear wave velocities for geotechnical layers were derived by statistically analyzing many seismic test data in Korea. The GTIS was used in a practical application to estimate site effects in the study area; seismic zoning maps of geotechnical earthquake parameters, such as the depth to bedrock and the site period, were created and presented as a regional synthetic strategy for earthquake risk assessment. Furthermore, seismic zonation of site classification was also performed to determine the site amplification coefficients for seismic design and seismic performance evaluation at any site and administrative sub-unit in the study area. The methodology and results of the case study of seismic zonations in the Seoul area verified that the GIS-based GTIS can be very useful for the regional estimation of seismic risk and also to support decisions regarding seismic hazard mitigation, particularly in the metropolitan area.     

ESTIMATION OF SITE AMPLIFICATION FROM DAHAN DOWNHOLE RECORDINGS

Four-level Dahan downhole recordings are used to investigate site amplification due to the near-surface structure. The downhole array is installed at the depths of 0 m (surface), 50, 100 and 200 m. Eight sets of data presenting well-separated incident and reflected waves (from the free surface) are selected to analyse the site amplification. Five techniques of spectral ratios have been tested to determine the most reliable way to estimate site response. The superior one is the ratio of Root-Mean-Square (RMS) Fourier amplitude spectra, because it has the minimum number of pseudopeaks. However, the pseudopeaks can be further reduced by omitting those ratios with small denominators. Results show that the amplification factors are nearly uniform at the depths of 50 and 100 m, where their values fall between one and two. This factor is frequency-dependent at the surface and in the range from two to five under 6 Hz. These results agree well with those from time-domain analyses and the transmission-coefficient estimation except that the latter overestimates the amplification factor of the surface station. This suggests that the shallow velocity structure provided by the geotechnical survey might be underestimated. Numerical modelling also supports the underestimation of velocity at near surface since the synthetic SH waves based on a modified velocity model can match these observations well.     

Site amplifications in the city of Benevento (Italy): comparison of observed and estimated ground motion from explosive sources

A hybrid approach is applied to construct 2-D synthetic seismograms for explosive sources. The computation of the signals has been performed with a hybrid technique, that couples the modal summation method, describing the P-SV propagation from the source to a target area, to the finite difference method, allowing the computation of the 2-D response of the target area. We use this technique to study the site response in the city of Benevento (Southern Italy). The results are compared both with those obtained for the same model, but considering an extended double-couple source and SH wave propagation and with those derived from a direct analysis of experimental data obtained during the Benevento Seismic Risk Project, funded by the Commission of European Communities. The type and depth of the source and the source-receiver distance can provide different excitations, but the response of the target area, in terms of ratios of response spectra, remains in general stable in the high frequency range (>2 Hz). The explosive source does not excite lower frequencies and is therefore unable to provide estimates of site effects in this frequency range. We also show that it is important to consider both P-SV and SH waves when assessing site effects. Finally, a very good definition of a structural model from geological data is fundamental in order to explain the observed data in the time domain. The observed data validate our synthetic modelling for evaluating the site effects and for a tentative seismic microzonation of the city of Benevento.     

ON THE POSSIBILITY OF SEISMIC EXPLORATION USING SURFACE TORQUE SOURCE AND RELATED TOPICS*

In this study we derive expressions for particle displacement or particle velocity anywhere inside a stratified earth and at its surface due to horizontal torque source located in the top layer. Equivalently, invoking Green's function reciprocity theorem, the solution applies also to the case of a surface or subsurface source when the resulting displacement or velocity is measured within the top layer. In order to evaluate the closed-form analytical solution economically and accurately it is advisable to introduce inelastic attenuation. Causal inelastic attenuation also lends the necessary realism to the computed seismic trace. To provide proof that the analytical solution is indeed correct and applicable to the multilayer case, a thick uniform overburden was assumed to consist of many thin layers. The correctness of the computed particle velocity response can be very simply verified by inspection. The computed response can also serve as a check on other less accurate methods of producing synthetic seismograms, such as the techniques of finite differences, finite elements, and various sophisticated ray-tracing techniques. It is not difficult to construct horizontal surface torque source. It appears that such source is well suited for seismic exploration in areas with a high-velocity surface layer. A realistic source function is analyzed in detail and normalized displacement response evaluated at different incidence angles in the near and the far fields. In an effort to distinguish the features of an SH torque seismogram from a pressure seismogram two models with identical layerings and layer parameters have been set up. As expected the torque seismogram is very different from the compressional seismogram. One desirable feature of a torque seismogram is the fast decay of multiples. Exact synthetic seismograms have many uses; some of them, such as the study of complex interference phenomena, phase change at wide angle reflection, channeling effects, dispersion (geometrical and material), absolute gain, and inelastic attenuation, can be carried out accurately and effortlessly. They can also be used to improve basic processing techniques such as deconvolution and velocity analysis. The numerical evaluation of the analytical solution of the wave equation as described in this paper has a long history. Most of the work leading to this paper was carried out by one of us (M. J. K.) in the years 1957 to 1968 at the Geophysical Research Corporation. However, the full testing of the various computer codes was carried out only very recently at the Phillips Petroleum Company.     

The road map for seismic risk analysis in a Mediterranean city

The seismic risk analysis evaluation in the Mediterranean area is one of the main tasks for the preservation of Cultural Heritage and for the sustainable development of Mediterranean cities. The Mediterranean area is characterised by a medium–high level of seismic risk, so that earthquakes are the major cause for the destruction of monuments, residential and industrial buildings. A case history regarding the seismic risk analysis for the city of Catania (Italy) is presented, since the city has been heavy damaged in the past by strong earthquakes such as the 1169 earthquake (XI MCS), the 1542 earthquake (IX MCS), the 1693 earthquake (XI MCS) and the 1818 earthquake (VIII MCS) etc., which caused several thousands of deaths. Fault modelling, attenuation laws, synthetic accelerograms, recorded accelerograms and site effects are considered for the evaluation of the seismic action. Vulnerability of physical environment, related to the presence of cavities and to seismic-induced landslides and liquefaction has been analysed, with special reference to the new modelling of such phenomena and to the application of models to given areas. Soil–structure Interaction has been analysed for some geotechnical works, such as shallow foundation and retaining wall, by means of physical and numerical modelling. The paper deals with the vulnerability of physical environment (landslides, liquefaction, etc.), while the road map continues with the analysis of vulnerability of monuments and buildings, with the aim of the estimation of the seismic resistance required to defend against the seismic action given by the scenario earthquake. For the mitigation of seismic risk, structural improvements of R.C. buildings with different methodology and techniques have been analysed, as well as the guideline for the strengthening of buildings. The work shows that the seismic risk of the city is not a summation of the seismic risk of each building, because the vulnerability of the urban system plays an important role on the seismic risk evaluation of a given city. To this aim the vulnerability of the road infrastructures, lifelines, and urban framework have been also analysed in the project.     

Virtual California: Fault Model, Frictional Parameters, Applications

Virtual California is a topologically realistic simulation of the interacting earthquake faults in California. Inputs to the model arise from field data, and typically include realistic fault system topologies, realistic long-term slip rates, and realistic frictional parameters. Outputs from the simulations include synthetic earthquake sequences and space-time patterns together with associated surface deformation and strain patterns that are similar to those seen in nature. Here we describe details of the data assimilation procedure we use to construct the fault model and to assign frictional properties. In addition, by analyzing the statistical physics of the simulations, we can show that that the frictional failure physics, which includes a simple representation of a dynamic stress intensity factor, leads to self-organization of the statistical dynamics, and produces empirical statistical distributions (probability density functions: PDFs) that characterize the activity. One type of distribution that can be constructed from empirical measurements of simulation data are PDFs for recurrence intervals on selected faults. Inputs to simulation dynamics are based on the use of time-averaged event-frequency data, and outputs include PDFs representing measurements of dynamical variability arising from fault interactions and space-time correlations. As a first step for productively using model-based methods for earthquake forecasting, we propose that simulations be used to generate the PDFs for recurrence intervals instead of the usual practice of basing the PDFs on standard forms (Gaussian, Log-Normal, Pareto, Brownian Passage Time, and so forth). Subsequent development of simulation-based methods should include model enhancement, data assimilation and data mining methods, and analysis techniques based on statistical physics.     

Seismic characterization of the shore sand at Catania

In the framework of the Catania Project supported by GNDT (Gruppo Nazionale per la Difesa dai Terremoti) seismic measurements have been carried out in the shore sands of Catania to define their detailed shear wave velocity (V) profiles and attenuation properties (Q). The method used, FTAN, is based on the analysis of the dispersion curve extracted from signals generated by controlled sources. The computation of complete synthetic seismograms enables the estimation of the attenuation properties of the surface layers. Hereinafter we analyse the shear seismic characteristics of the shore sands as indicators of potential liquefaction phenomena.     

黄土地区地铁地下结构抗震研究综述

介绍了地铁地下结构抗震研究的发展历史、研究现状及黄土地区地铁地下结构抗震的相关研究,总结归纳出地震作用下地铁地下结构的地震破坏形式与地震反应特征;对地铁地下结构的不同抗震研究方法进行对比研究,指出各种研究方法的优缺点与适用情况。针对我国黄土地区地铁地下结构抗震分析理论研究不足,基础资料积累欠缺,抗震设计尚无可靠依据的实际情况,探讨了黄土动力本构模型、黄土与结构动力相互作用、离心机振动台试验研究及抗震分析方法等因素对黄土地区地铁地下结构抗震研究可靠性的影响,提出黄土地区地铁地下结构地震动力反应分析、地震破坏机理分析、抗震性能分析及抗震设计中需要深入研究的关键问题。     

Topographic effects on the Monte Po Hill in Catania (Italy)

The study regards an evaluation of site effects on and near the Monte Po hill, located in the north-eastern part of the city of Catania (Italy), an area at high seismic risk. At the beginning of 2007 a seismic station was located in a school building, situated at the slope toe, but no seismic events have been recorded as yet. Therefore, synthetic seismograms have been used to evaluate the ground response analysis at the surface. Because the average slope is moderate (less than 15°), 1-D computer codes have been used to model the equivalent-linear earthquake site response analyses of layered hill deposits, as generally performed by professionals. However, the slope response has also been analysed in greater detail, using a 2-D computer code and the soil characterisation has been evaluated accurately by means of borings, Down-Hole tests, SDMT tests and laboratory tests. Comparing 1-D with 2-D results the stratigraphic site amplification and the Topographic Aggravation Factor (TAF) have also been computed. The aim of the study is that it will form a basis for the design of works to remediate the damage caused by a landslide reactivated by the earthquake in Eastern Sicily on December 13, 1990 (M_L=5.4).     

Revisiting the 1669 Etnean eruptive crisis using a cellular automata model and implications for volcanic hazard in the Catania area

Cellular Automata provide an alternative approach to standard numerical methods for modelling some complex natural systems, the behaviour of which can be described in terms of local interactions of their constituent parts. SCIARA is a 2-D Cellular Automata model which simulates lava flows. It was tested on, validated by, and improved on several Etnean lava events such as the 1986–1987 eruption and the first and last phase of the 1991–1993 event. With respect to forecasting the surface covered by the lava flows, the best results were acceptable. The model has been used to determine hazard zones in the inhabited areas of Nicolosi, Pedara, S. Alfio and Zafferana (Sicily, Italy). The main goal of the current work in the Etnean area from Nicolosi to Catania has been the verification of the volcanic hazard effects of an eruptive crisis similar to the event that occurred in 1669. The simulation uses the volcanic data of the 1669 eruption with present-day morphology. Catania has been affected by some historical Etnean events, the most famous one being the 1669 eruption, involving 1 km³ of lava erupted over the course of 120 days. The simulation of ephemeral vents and the use of different histories within the experiments have been crucial in the determination of a new hazard area for Catania. In fact, during the simulation the city was never affected without the introduction of ephemeral vents, proving the fact that lava tubes played a fundamental role in the 1669 Catania lava crisis.     

Local seismic site amplification: effects of obliquely incident antiplane motions

Seismic site amplification studies are generally used to assess the effects of local geology and soil conditions on ground motion characteristics. Although extensive reviews on site amplification phenomena associated with stratigraphic effects can be found in the specialized literature, it should be pointed out that most of the practical applications have been limited to the study of vertically propagating shear horizontal (SH) waves, i.e., to the 1-D soil amplification problem. Furthermore, little attention, if any, has been devoted to the study of the effects of non-vertically incident SH waves on surface accelerograms and on the earthquake response of structures. In the present work, the study is extended to an investigation of 2-D site amplification of non-vertically propagating seismic shear waves in multilayered viscoelastic soil deposits. Sensitivity analyses of the effects of non-vertical incidence on site amplification functions are performed based on site geotechnical data collected from post-seismic investigations of the 1980 El-Asnam earthquake. Analytical results are discussed in terms of seismic site transfer functions, spectral ratios, surface acceleration time histories, and structural response spectra for different values of wave incidence angle. Both bedrock and rock outcropping cases are examined.     

Thessaloniki's Detailed Microzoning: Subsurface Structure as Basis for Site Response Analysis

—?The city of Thessaloniki is located in northern Greece, close to a large seismogenic area and it has experienced several destructive earthquakes during the present century. In this paper, we focus on the definition of the subsurface structure of the city from a site response analysis perspective. The paper presents, together with a summary of geology and seismicity, the results of a large-scale geophysical and geotechnical survey, in order to determine and validate geometry and dynamic properties of the main soil formations. The synthesis and combination of recent results regarding the dynamic properties with those obtained from the elaboration of a large database of classical geotechnical tests led to the design of a detail geotechnical map and of various 1-D profiles, 2-D cross sections and 3-D thematic maps for the main soil formations. These soil profiles and maps are oriented to site effect studies and provide a comprehensive picture easily adapted to geographic information systems (GIS) for planning and design purposes. Moreover, the results of this study were correlated with macroseismic observations reported in many earlier published microzonation studies of Thessaloniki. These comparative observations revealed the complexity of surface geology of the urban area, a fact which is expected to provoke additional amplification with respect to 1-D resonance.     

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.     

Mapping soil effects in Barcelona, Spain, through an integrated GIS environment

Maps of soil response for the city of Barcelona have been for the first time obtained through a GIS environment that integrates the different procedures for soil response estimation within a single tool. These maps constitute part of the results of the local scale application of the computer prototype for seismic risk assessment that was developed within the European project SERGISAI. The approach involves: collection of available data relevant to local geology, implementation of geotechnical models, estimation of the reference seismic action, generation of synthetic strong-motion time histories, and soil response calculation through 1D analytical method. The resulting predictive hazard maps of predominant period and amplification ratio delineate potential variations on ground shaking and constitute a first approximation towards an integrated approach to Barcelona urban area microzonation. Analysis of the observed differences, when comparing the analytical results in this study with previous empirical studies, provides a useful feedback to establish site dependence suitability and reliability of methods, to extract information on at-present inaccessible parameters needed for the characterisation of physical properties of soil, and also to delimit those areas where further in-depth survey research is needed for a proper seismic hazard assessment.     

Delineation of Near-Surface Structure in the Southern Part of 15th of May City, Cairo, Egypt Using Geological, Geophysical and Geotechnical Techniques

The integration of geological, geophysical, and geotechnical interpretation at the southern part of 15th of May City, have been used to evaluate the subsurface stratigraphy, especially the clay layer which may cause serious danger to construction. Those techniques have been used to delineate the subsurface structures as normal faults, which play a critical role on the stability of buildings. Geological setting of the area has been evaluated through the construction of a geological map from different geological sections and samples obtained from more than 30 observation points. Geophysical tools such as vertical electrical soundings (12 VESs), 2-D dipole–dipole array (7 sections), P-wave shallow seismic refraction (31 profiles) and multiple channel analysis of surface waves (31 MASW profiles) have been carried out to image the subsurface situation. Geotechnical evaluation using 26 boreholes, samples, laboratory tests and geotechnical parameters has been done at the area of interest. The geological setting demonstrates that the city had been constructed on the second and third members of Qurn Formation (Upper Eocene) composed of argillaceous limestone, marl and shale. Two normal faults are passing through the area were observed. The resistivity (VES and dipole–dipole) and seismic (P-waves and MASW) results reflect the presence of the two normal faults cross the study area, affecting the obtained section of marl, clayey marl and limestone layers. The geotechnical information indicate the presence of the normal faults and the existence of clay layer with swelling ability reaching 140%, which may cause cracks in the upper layers and/or subsidence.