A deterministic seismic hazard analysis for shallow earthquakes in Greece

The maximum expected ground motion in Greece is estimated for shallow earthquakes using a deterministic seismic hazard analysis (DSHA). In order to accomplish this analysis the input data include an homogeneous catalogue of earthquakes for the period 426 BC–2003, a seismogenic source model with representative focal mechanisms and a set of velocity models. Because of the discrete character of the earthquake catalogue and of errors in location of single seismic events, a smoothing algorithm is applied to the catalogue of the main shocks to get a spatially smoothed distribution of magnitude. Based on the selected input parameters synthetic seismograms for an upper frequency content of 1 Hz are computed on a grid of 0.2° × 0.2°. The resultant horizontal components for displacement, velocity, acceleration and DGA (Design Ground Acceleration) are mapped. The maps which depict these results cannot be compared with previously published maps based on probabilistic methodologies as the latter were compiled for a mean return period of 476 years. Therefore, in order to validate our deterministic analysis, the final results are compared with PGA estimated from the maximum observed macroseismic intensity in Greece during the period 426 BC–2003.Since the results are obtained for point sources, with the frequency content scaled with moment magnitude, some sensitivity tests are performed to assess the influence of the finite extent of fault related to large events. Sensitivity tests are also performed to investigate the changes in the peak ground motion quantities when varying the crustal velocity models in some seismogenic areas. The ratios and the relative differences between the results obtained using different models are mapped and their mean value computed. The results highlight the importance in the deterministic approach of using good and reliable velocity models.     

Attenuation laws and seismic hazard assessment

In the present study, the seismic hazard in the city of Patras, central Greece, is estimated. The computations are based on a slightly modified version of the method proposed by Cornell, allowing the introduction of individual attenuation laws for each seismic source.The obtained results emphasize the dependence of hazard on attenuation and the importance of the use of local attenuation laws in seismic hazard assessment.Paper presented at the 21st General Assembly of the European Seismological Commission, held in Sofia, 1988.     

Ground zoning against seismic hazard in Athens, Greece

The 1999, Ms=5.9, Athens earthquake caused serious structural damage to buildings in the western part of Athens, Greece. This paper presents the ground zoning against seismic hazard proposed shortly after the earthquake in order to aid reconstruction of the area. Existing engineering geological and geotechnical data were combined with local observations to provide a unified set of classification criteria, consistent with provisions of the Greek Seismic Code EAK. The accuracy and the possible limitations of this zoning procedure are addressed through comparison with observed damage distribution as well as results from seismic ground response analyses performed at sites with well established soil profiles. There is clear evidence that the proposed zones correspond to geological formations exhibiting grossly different seismic response with regard to the design of common engineering structures. However, the mostly qualitative nature of the guidelines for ground categorisation provided by EAK and the general lack of systematic, site-specific geotechnical data for the whole area induce uncertainties in the definition of the seismic design actions for the different zones. These objective uncertainties certainly demand increased conservatism but do not limit application of the proposed methodology for first aid, preliminary planning in the event of destructive earthquakes.     

A Moore's cellular automaton model to get probabilistic seismic hazard maps for different magnitude releases: A case study for Greece

Cellular automata are simple mathematical idealizations of natural systems and they supply useful models for many investigations in natural science. Examples include sandpile models, forest fire models, and slider block models used in seismology. In the present paper, they have been used for establishing temporal relations between the energy releases of the seismic events that occurred in neighboring parts of the crust. The catalogue is divided into time intervals, and the region is divided into cells which are declared active or inactive by means of a threshold energy release criterion. Thus, a pattern of active and inactive cells which evolves over time is determined. A stochastic cellular automaton is constructed starting with these patterns, in order to simulate their spatio-temporal evolution, by supposing a Moore's neighborhood interaction between the cells. The best model is chosen by maximizing the mutual information between the past and the future states. Finally, a Probabilistic Seismic Hazard Map is given for the different energy releases considered. The method has been applied to the Greece catalogue from 1900 to 1999. The Probabilistic Seismic Hazard Maps for energies corresponding to m = 4 and m = 5 are close to the real seismicity after the data in that area, and they correspond to a background seismicity in the whole area. This background seismicity seems to cover the whole area in periods of around 25–50 years. The optimum cell size is in agreement with other studies; for m > 6 the optimum area increases according to the threshold of clear spatial resolution, and the active cells are not so clustered. The results are coherent with other hazard studies in the zone and with the seismicity recorded after the data set, as well as provide an interaction model which points out the large scale nature of the earthquake occurrence.     

Seismicity, seismotectonics and seismic hazard of Italy

A first generation of probabilistic seismic hazard maps of the Italian country are presented. They are based on seismogenic zoning deriving from a kinematic model of the structural tectonic units and on an earthquake catalogue with the foreshock and aftershock events filtered out. The following ground motion parameters have been investigated and mapped using attenuation equations based on strong-motion recordings of Italian earthquakes: peak ground acceleration and velocity; Arias intensity; strong motion duration; and the pseudovelocity and pseudoacceleration spectral values at 14 fixed frequencies both for the vertical and the largest horizontal component. A Poissonian model of earthquake occurrence is assumed as a default and the hazard maps are presented in terms of ground motion values expected to be exceeded at a 10% probability level in 50 years (return period 475 years) according to the requirement of Eurocode 8 for the seismic classification of national territories, as well as in terms of exceedance probabilities of selected ground motion values. Finally, as a tentative study, the use of hybrid methods (implementing both seismogenic zones and structures), renewal processes (including earthquake forecasting) and the influence of site effects (as the basis for the planning of earthquake scenarios) were explored.     

An estimate of probabilistic seismic hazard for five cities in Greece by using the parametric-historic procedure

A probabilistic procedure was applied to assess seismic hazard for the sites of five Greek cities (Athens, Heraklion, Patras, Thessaloniki and Volos) using peak ground acceleration as the hazard parameter. The methodology allows the use of either historical or instrumental data, or a combination of both. It has been developed specifically for the estimation of seismic hazard at a given site and does not require any specification of seismic sources or/and seismic zones. A new relation for the attenuation of peak ground acceleration was employed for the shallow seismicity in Greece. The computations involved the area- and site-specific parts. When assessing magnitude recurrence for the areas surrounding the five cities, the maximum magnitude, m_max, was estimated using a recently derived equation. The site-specific results were expressed as probabilities that a given peak ground acceleration value will be exceeded at least once during a time interval of 1, 50 and 100 years at the sites of the cities. They were based on the maximum peak ground acceleration values computed by assuming the occurrence of the strongest possible earthquake (of magnitude m_max) at a very short distance from the site and using the mean value obtained with the help of the attenuation law. This gave 0.24 g for Athens, 0.53 g for Heraklion (shallow) and 0.39 g Heraklion (intermediate-depth seismicity), 0.30 g for Patras, 0.35 g for Thessaloniki and 0.30 g for Volos. In addition, the probabilities of exceedance of the estimated maximum peak ground acceleration values were calculated for the sites. The standard deviation of the new Greek attenuation law demonstrates the uncertainty and large variation of predicted peak ground acceleration values.     

Seismic hazard in the Sannio-Matese area of Italy

The assessment of seismic hazard at five selected sites in the Sannio-Matese region is based on the computer program SRAMSC. Owing to the extensive historical data base for the output parameter, the MSK intensiy is chosen. The seismicity model is made up of five narrow area seismic sources. Circular or elliptical macroseismic fields are assigned to individual sources. A generalized Kövesligethy equation is used for this purpose as the attenuation relationship. The study reveals similar and a rather high hazard at the sites at Benevento, Boiano, and Melfi, which are located in the zone of highest seismic activity. At the Pomigliano and Lucera sites, the assessed hazard is much lower.     

A seismic source zone model for the seismic hazard assessment of the Italian territory

We designed a new seismic source model for Italy to be used as an input for country-wide probabilistic seismic hazard assessment (PSHA) in the frame of the compilation of a new national reference map.

We started off by reviewing existing models available for Italy and for other European countries, then discussed the main open issues in the current practice of seismogenic zoning.

The new model, termed ZS9, is largely based on data collected in the past 10 years, including historical earthquakes and instrumental seismicity, active faults and their seismogenic potential, and seismotectonic evidence from recent earthquakes. This information allowed us to propose new interpretations for poorly understood areas where the new data are in conflict with assumptions made in designing the previous and widely used model ZS4.

ZS9 is made out of 36 zones where earthquakes with M_w > = 5 are expected. It also assumes that earthquakes with M_w up to 5 may occur anywhere outside the seismogenic zones, although the associated probability is rather low. Special care was taken to ensure that each zone sampled a large enough number of earthquakes so that we could compute reliable earthquake production rates.

Although it was drawn following criteria that are standard practice in PSHA, ZS9 is also innovative in that every zone is characterised also by its mean seismogenic depth (the depth of the crustal volume that will presumably release future earthquakes) and predominant focal mechanism (their most likely rupture mechanism). These properties were determined using instrumental data, and only in a limited number of cases we resorted to geologic constraints and expert judgment to cope with lack of data or conflicting indications. These attributes allow ZS9 to be used with more accurate regionalized depth-dependent attenuation relations, and are ultimately expected to increase significantly the reliability of seismic hazard estimates.     

苏南地区地质灾害区划评价

苏南是江苏省经济发展水平较高的地区,也是全国人口密度较大、经济较发达和城镇化程度较高的地区,地质灾害是影响苏南快速发展的重要因素之一。研究采用层次分析法,尝试以重要性、危害性和易发性3个要素作为指标建立评价体系,对苏南地区进行地质灾害区域划分,并在区划基础上考虑社会经济因素对地质灾害危害程度的影响,将苏南地区划分出地面沉降重点防治区和滑坡塌陷重点防治区2个重点防治区,其他区域为一般防治区。评价结果对于苏南地区地质灾害防治工作具有较大的参考价值。     

Probabilistic seismic hazard zonation of Syria

Earthquake hazard maps for Syria are presented in this paper. The Peak Ground Acceleration (PGA) and the Modified Mercalli Intensity (MMI) on bedrock, both with 90% probability of not being exceeded during a life time of 50, 100 and 200 years, respectively are developed. The probabilistic PGA and MMI values are evaluated assuming linear sources (faults) as potential sources of future earthquakes. A new attenuation relationship for this region is developed. Ten distinctive faults of potential earthquakes are identified in and around Syria. The pertinent parameters of each fault, such as theb-parameter in the Gutenberg-Richter formula, the annual rate ₄ and the upper bound magnitudem ₁ are determined from two sets of seismic data: the historical earthquakes and the instrumentally recorded earthquake data (AD 1900–1992). The seismic hazard maps developed are intended for preliminary analysis of new designs and seismic check of existing civil engineering structures.     

Probabilistic seismic hazard maps for Panama

Probabilistic seismic hazard maps in term of Modified Mercalli (MM) intensity are derived by applying the Cornell-McGuire method to four earthquake source zones in Panama and adjacent areas. The maps contain estimates of the maximum MM intensity for return periods of 5, 25 and 100 yr. The earthquake phenomenon is based on the point source model. The probabilistic iso-intensity map for a return period of 50 yr indicates that the Panama Suture Zone (PSZ) could experience a maximum (MM) intensity IX, and the Panama Fracture Zone (PFZ) an MM intensity VIII, for the rest of the area this varies from IV up to VIII. The present study intends to serve as a reference for more advanced approaches, to stimulate discussions and suggestions on the data base, assumptions and inputs, and path for the risk based assessment of the seismic hazard in the site selection and in the design of common buildings and engineering.     

苏南地区地质灾害区划评价

苏南是江苏省经济发展水平较高的地区,也是全国人口密度较大、经济较发达和城镇化程度较高的地区,地质灾害是影响苏南快速发展的重要因素之一。研究采用层次分析法,尝试以重要性、危害性和易发性3个要素作为指标建立评价体系,对苏南地区进行地质灾害区域划分,并在区划基础上考虑社会经济因素对地质灾害危害程度的影响,将苏南地区划分出地面沉降重点防治区和滑坡塌陷重点防治区2个重点防治区,其他区域为一般防治区。评价结果对于苏南地区地质灾害防治工作具有较大的参考价值。     

Some aspects of the seismic hazard in Panama City

The definition of earthquake sources in the Panama region on the basis of both tectonics and average seismicity rates, have recently led to the concept of a microplate surrounded by seismically active areas. The effects of these earthquakes on the place where the most important concentration of investments and population is located, the capital city of Panama, are analyzed in this paper using statistical approaches.The parameters of Gumbel's Type-I distribution of extreme values for a continuous interval of 60 yr annual maximum magnitudes were used to make probabilistic estimations of the seismic hazard in Panama City. An earthquake with magnitude 7.5 is capable of producing a modified Mercalli intensity VII in Panama City, provided the source distance is of the order of 100 km. This earthquake has a probability of occurrence of 69% in 50 yr.     

On the reliability of different methods of seismic hazard assessment in Greece

Historical and present century instrumental data have been used to determine seismic hazard in 35 sites of Greece by the application of Cornell's method (Cornell, 1968) and the mean value method. The macroseismic intensity has been considered as a measure of seismic hazard. Comparison of the results of the two methods showed that, in general, the mean value method gives higher values, particularly for low probabilities of exceedance. In addition, for some sites, the differences of the expected intensities resulting from the two methods, indicate that finer tuning of the seismogenic souce model is required, or suggest time dependence. Although each one of these methods has its own merits, the method based on seismic zonation (Cornell's method) has several advantages and must be preferred when an accurate zonation is possible by the use of macroseismic and instrumental seismic data, together with geological and geomorphological information. However, reliable estimates of seismic hazard at a particular site require work on a microzoning scale, incorporating historical, archaeological, and recent geological data.Paper presented at the 21st General Assembly of the European Seismological Commission held in Sofia, 1988.     

Influence of attenuation models in seismic hazard evaluation for the North-Eastern Sicily area

From the influence already revealed by attenuation models on the values of expected intensity for North-Eastern Sicily, the necessity arises to quantify the weight of these models and of their respective coefficients on the projection of intensity.A first evaluation is presented in this paper using the Sponheuer, Blake, and Grandori models.A comparison of the expected intensity maps allows a first critical estimate, showing the greater adaptability of the Grandori model to describe the attenuation of intensity for the investigated area.     

Different approaches to the seismic hazard of Sannio-Matese (Southern Italy)

The seismic hazard of Sannio-Matese is calculated using a new seismogenetic zoning of southern Italy of seven areas and different, azimuth-dependent, attenuation laws. Various approaches (Gumbel's first and third asymptotic distribution, Cornell) lead to similar results for the different exposure times considered (100 and 200 years) and probability levels (37 and 68%). The present seismic regulation proposed in 1980 by the Italian National Council for Researches, and based on a different approach, is, in general, confirmed by the results.     

The Olmsted fault zone, southernmost Illinois: A key to understanding seismic hazard in the northern new Madrid seismic zone

Geological deformation in the northern New Madrid seismic zone, near Olmsted, Illinois (USA), is analyzed using integrated compressional-wave (P) and horizontally polarized-wave (SH) seismic reflection and regional and dedicated borehole information. Seismic hazards are of special concern because of strategic facilities (e.g., lock and dam sites and chemical plants on the Ohio River near its confluence with the Mississippi River) and because of alluvial soils subject to high amplification of earthquake shock. We use an integrated approach starting with lower resolution, but deeper penetration, P-wave reflection profiles to identify displacement of Paleozoic bedrock. Higher resolution, but shallower penetration, SH-wave images show deformation that has propagated upward from bedrock faults into Pleistocene loess. We have mapped an intricate zone more than 8 km wide of high-angle faults in Mississippi embayment sediments localized over Paleozoic bedrock faults that trend north to northeast, parallel to the Ohio River. These faults align with the pattern of epicenters in the New Madrid seismic zone. Normal and reverse offsets along with positive flower structures imply a component of strike-slip; the current stress regime favors right-lateral slip on northeast-trending faults. The largest fault, the Olmsted fault, underwent principal displacement near the end of the Cretaceous Period 65 to 70 million years ago. Strata of this age (dated via fossil pollen) thicken greatly on the downthrown side of the Olmsted fault into a locally subsiding basin. Small offsets of Tertiary and Quaternary strata are evident on high-resolution SH-wave seismic profiles. Our results imply recent reactivation and possible future seismic activity in a critical area of the New Madrid seismic zone. This integrated approach provides a strategy for evaluating shallow seismic hazard-related targets for engineering concerns.     

The importance of paleoseismology in seismic hazard studies for critical facilities

This paper is intended to provide a perspective on the use of paleoseismological studies in the seismic hazard assessment of critical facilities, such as dams, chemical/petrochemical facilities and nuclear power plants. In particular, the use of data obtained from paleoseismological studies for probabilistic seismic hazard analyses, when the required probabilities of exceedance are very low (e.g. 10^{− 6}–10^{− 7}) is considered. Recent revisions to the IAEA Safety Standards that provide guidance to Member States in their work related to the seismic safety of nuclear power plants are presented to illustrate the importance of this emerging discipline.     

地震危险性、地震危害性和地震易损性

做好城市防震减灾工作 ,关键是对未来可能遭遇地震灾害定量化预测。本文论述地震灾害定量化的三要素 :地震危险性 ,地震危害性和地震易损性 ,它们概念完全不同 ,且很容易混淆 ,但又存在因与果的关系。本文还简要介绍厦门市地震科技工作者开展的“闽南地区综合防震减灾示范工程”。