Aristotle University of Thessaloniki,School of Geology,Geophysical Laboratory,Thessaloniki,Greece

Parameters of seismic hazard are estimated by the application of the maximum likelihood method. The technique is based on a procedure which utilizes data of different quality, e.g., the ones where the uncertainty in the assessment of the magnitudes is great and those where the magnitudes are computed with great precision. In other words, the data were extracted from both historical (incomplete) and recorded (complete) files. The historical part of the catalogue contains only the strongest events, whereas the complete part can be divided into several subcatalogues each one assumed to be complete above a specified threshold magnitude. Uncertainty in the determination of magnitudes has also been taken into account. The method allow us to estimate the seismic hazard parameters which are the maximum regional magnitude, Mmax , the activity rate, lˆ, of the seismic events and the well known b-value, the slope of the magnitude-frequency relationship. The parameter b, which is interrelated to b (b = bloge), is also obtained. All these parameters are of physical significance. The mean Return Periods, RP, of earthquakes with a certain lower magnitude M ³ m are also determined. The method is applied in some regions of the circum-Pacific belt, which includes various tectonic features, and where catastrophic earthquakes are known from the historical era. The seismic hazard level is also calculated as a function of the form q(Mmax , RP7.5 ) and a relative hazard scale (defined as an index K) is defined for each seismic region. According to this, the investigated regions are classified into five groups of very low, low, intermediate, high and very high seismic hazard levels. This classification is useful for both theoretical and practical reasons and provides a picture of quantitative seismicity.

Global data have been widely used for seismicity and seismic hazard assessment by seismologists. In the present study we evaluate worldwide seismicity in terms of maps of maximum observed magnitude (Mmax), seismic moment (M 0 ) and seismic moment rate (M 0S). The data set used consists of a complete and homogeneous global catalogue of shallow (h £ 60 km) earthquakes of magnitude MS ³ 5.5 for the time period 1894-1992. In order to construct maps of seismicity and seismic hazard the parameters a and b derived from the magnitude-frequency relationship were estimated by both: a) the least squares, and b) the maximum likelihood, methods. The values of a and b were determined considering circles centered at each grid point 1° (of a mesh 1° ´1°) and of varying radius, which starts from 30 km and moves with a step of 10 km. Only a and b values which fulfill some predefined conditions were considered in the further procedure for evaluating the seismic hazard maps. The obtained worldwide M max distribution in general delineates the contours of the plate boundaries. The highest values of M max observed are along the circum-Pacific belt and in the Himalayan area. The subduction plate boundaries are characterized by the largest amount of M 0 , while areas of continental collision are next. The highest values of seismic moment rate (per 1 year and per equal area of 10 000 km 2) are found in the Southern Himalayas. The western coasts of U.S.A., Northwestern Canada and Alaska, the Indian Ocean and the eastern rift of Africa are characterized by high values of M 0 , while most of the Pacific subduction zones have lower values of seismic moment rate. Finally we analyzed the seismic hazard in South America comparing the predicted by the NUVEL1 model convergence slip rate between Nazca and South America plates with the average slip rate due to earthquakes. This consideration allows for distinguishing between zones of high and low coupling along the studied convergence plate boundary.

A probabilistic approach was applied to map the seismic hazard in Greece and the surrounding region. The procedure does not require any specification of seismic sources or/and seismic zones and allows for the use of the whole seismological record, comprising both historical and instrumental data, available for the region of interest. The new seismic hazard map prepared for Greece and its vicinity specifies a 10% probability of exceedance of the given Peak Ground Acceleration (PGA) values for shallow seismicity and intermediate soil conditions for an exposure time of 50 years. When preparing the map, the new PGA attenuation relation given by Margaris et al. (2001) was employed. The new map shows a spatial distribution of the seismic hazard that corresponds well with the features of shallow seismicity within the examined region. It depicts the level of seismic hazard in which the exceedance of the PGA value of 0.25 g may be expected to occur within limited areas. The highest estimated levels of seismic hazard inside the territory of Greece are found in the Northern Sporades Islands, where PGA values in excess of 0.50 g are reached at individual sites, and in the Zante Island in Western Greece, where PGA values in the range of 0.35 g to 0.40 g are obtained at more numerous localities. High values are also observed in the sea between the Karpathos and Rhodes islands, near the Island of Amorgos (Cyclades Archipelago) and in the Southwestern Peloponnesus. The levels of seismic hazard at the sites of seven Greek cities (Athens, Jannena, Kalamata, Kozani, Larisa, Rhodes and Thessaloniki) were also estimated in terms of probabilities that a given PGA value will be exceeded at least once during a time interval of 1, 50 and 100 years at those sites. These probabilities were based on the maximum horizontal PGA values obtained by applying the design earthquake procedure, and the respective median values obtained were 0.24 g for Athens, 0.28 g for Jannena, 0.30 g for Kalamata, 0.21 g for Kozani, 0.24 g for Larisa, 0.43 g for Rhodes and 0.35 g for Thessaloniki. The probabilities of exceedance of the estimated maximum possible PGA value were also calculated for the cities to illustrate the uncertainty of maximum PGA assessment.

Seismic hazard in terms of probability of exceedance of a given intensity in a given time span,was assessed for 12 sites in Japan.The method does not use any attenuation law.Instead,the dependence of local intensity on epicentral intensity I 0 is calculated directly from the data,using a Bayesian model.According to this model (Meroni et al., 1994),local intensity follows the binomial distribution with parameters (I 0 ,p ).The parameter p is considered as a random variable following the Beta distribution.This manner of Bayesian estimates of p are assessed for various values of epicentral intensity and epicentral distance.In order to apply this model for the assessment of seismic hazard,the area under consideration is divided into seismic sources (zones)of known seismicity.The contribution of each source on the seismic hazard at every site is calculated according to the Bayesian model and the result is the combined effect of all the sources.High probabilities of exceedance were calculated for the sites that are in the central part of the country,with hazard decreasing slightly towards the north and the south parts.