Seismic Hazard Assessment of Iran

The development seismic hazard map of Iran is based on probabilistic seismic hazard computation using the historical earthquakes data, geology ,tectonics, fault activity and seismic source models in Iran. These maps have been prepared to indicate the earthquake hazard of Iran in the forms of iso-acceleration contour lines, and seismic hazard zonations by utilizing to current probabilistic procedures. They contain the probabilistic estimates of Peak Ground Acceleration for the return periods of 75 and 475 years. The map has been divided into intervals of 0.25 degrees in both latitudinal and longitudinal directions for calculating the peak ground acceleration values at each grid point and drawing seismic hazard curves. The results presented in this study are to provide a base for preparing of seismic risk map, the estimation of earthquake insurance premiums, and the preliminary site evaluation of critical facilities.


Introduction
The probabilistic assessment of the seismic hazard in Iran involves specifying the likelihood ,magnitude, location, and nature of earthquakes that might have damaging effect in the region or at the site and estimating the peak acceleration of the ground shaking.The basis for all seismic hazard assessment is the analysis of seismicity or the occurrence of earthquake in space and time.The term seismicity in this study has been used to describe different parts of earthquake record.These parts include : 1 ( the pre-historic and historical seismicity data (pre-1900 ( ,

(
The instrumental seismicity data ) 1900 -1997 ( , a n d 3 ( Seismic sources of the historical and instrumental earthquakes, rate of activity, and fault interaction.Accurate knowledge of seismicity is an important tool for understanding active tectonics .Seismic hazard analysis require an assessment of the future earthquake potential in Iran.It is, therefore, necessary to estimate the maximum earthquake magnitude and recurrence character that might be generated by a particular active fault.The most common uses of more detailed geologic data have been to constrain maximum earthquake magnitudes using empirical relationships between earthquake rupture dimensions and magnitude.
One of the basic elements in assessing seismic hazards is to recognize seismic sources that could affect the particular location at which the hazard is being evaluated.These sources are often called seismotectonic sources.Defining and understanding seismotectonic sources is often the major part of a seismic hazard analysis and requires knowledge of the regional and local geology, seismicity ,and tectonics.Two types of seismotectonic sources were identified in this study, fault or line sources, which may be sources of future earthquakes, and localizing structures and seismotectonic provinces which appear as area sources .
In this study, the geological and seismological data were integrated into a probabilistic seismic hazard model for Iran which located between 25 to 40 degrees N and 44 to 63 degrees E .The main steps and parameters of the probabilistic seismic hazard assessment conducted in this study have been summarized as follows : 1. Macroseismic ,historical and instrumental seismic data that provided information on the regional stress and the location and frequency of earthquakes occurrence during the past several hundred years.
2. Seismic source model that based on geological and seismological evidence are identified.These sources described the potential locations of future earthquakes within an area.

3.
Source seismicity information which were estimated recurrence relationships and the assessment of the maximum magnitudes for each individual seismic source.
4. Attenuation of the strong ground motion that provided information on the variation of ground motion levels with earthquake magnitude, and source-to-site distance.

5.
Recurrence forecasting which includes statistical methods used in the estimation of source seismicity information and stochastic model (Homogeneous Poisson Process) used in the estimation of seismic occurrences.

Earthquake Data Base
A review of Iran's historical earthquakes (Pre-1900) is provided by Ambraseys and Melville ( 1982 .(To improve the quality and accuracy of these earthquake data the new Earthquake Catalogue of Iran has been compiled by Berberian in 1994.For present century ,IIEES catalogue based on reports from international seismological institutes have been used.With the inception of seismological stations, the location accuracy of the earthquakes has considerably improved after mid 1960's.Thus, especially for Iran, it can be stated that the reliable earthquake data base exists only for the last several decades.Both of these catalogues have provided a base for the spatial correlation of the seismicity with the seismotectonic sources in Iran .From 1930 to 1985, the seismicity analysis of Iran was conducted by Wilson(1930), Niazi and Basford (1968), Nowroozi (1972Nowroozi ( , 1976))

Seismotectonic Provinces of Iran
The Iranian plateau can be characterized by active faults, recent volcanic and high surface elevation along the Alpied earthquake belt.Tectonic studies indicate that the Iranian plateau has a very high density of active and recent faults.Earthquake data of Iran show that most activity is concentrated along the Zagros fold thrust belt while less activity is observed in central and eastern Iran.Thus ,several regions are vulnerable to destructive earthquakes.The preparation of an earthquake hazard map is the delineation of the seismotectonic province and the assessment of the associated maximum earthquake potentials.They can be defined as a geographic region of some geological, geophysical and seismological similarity with the assumption of uniform earthquake potential.Earthquakes are assumed to occur randomly throughout the seismotectonic provinces event though the earthquake record may indicate some clustering at preferred locations.The seismotectonic province of Iran is defined as an area bounded by geological features which mark a difference in seismic characteristics of one province from its neighboring provinces.Each province has equal seismic potential and uniform geological structure and trends .
The seismotectonic structure of Iran was studied by several investigators in the past .Stoklin(1968), Takin(1972), and Berberian (1976) have suggested simplified divisions consisting of only nine, four, and four regions, respectively.A more elaborated division, consisting of twenty-three seismotectonic provinces was suggested by Nowroozi (1976).In the present work, a modified and updated catalogue of large and damaging Iranian earthquakes has been prepared to provide a based for the correlation between seismic activity and seismotectonic provinces.On this basis, the new seismotectonic provinces have been proposed by Tavakoli ( 1996 .( The boundaries of the provinces are established through analysis of seismic history, relocated epicenter for the past several decades, tectonic environments, active faults ,regional geomorphology, and plate boundaries.Tavakoli (1996) has divided Iran into 20 seismotectonic provinces, as it has been shown in Figure (2).The activity faults are located in parts of twenty known seismotectonic provinces of Iran.

Earthquake Hazard Parameters
Earthquake hazard analysis requires and assessment of earthquake hazard parameters and the future of earthquake potential in a region.Earthquake hazard parameters such as ,maximum expected magnitude, Mmax, activity rate , λ ,and b value of the Gutenberg-Richter relation have been evaluated for each seismotectonic province of Iran.The applied maximum likelihood method used in this study ) Kijko and Sellevoll, 1992) permits the combination of both historical and instrumental data.The maximum likelihood estimation for fitting the seismicity parameter ) β ( has a good applicability correlation with the Iranian earthquake data.This is due to the fact that earthquake magnitudes have always been reported by some uncertainty.In this method also artificially homogeneous data set is simulated through the determination of the period over which the data in a given time span are completely reported.In the present study, unlike previous works, seismic gaps, i.e., when records are missing or the seismic network was not in operation, and uncertainties of earthquake magnitudes are considered in analysis.This consideration is necessary for regions like Iran where few Earthquake data base are available .Therefore, based on the collected data, the earthquake hazard parameters, as shown in Table 1, were estimated for each seismotectonic province of Iran.It should be noted that in this work, the Gutenberg-Richter parameters was assumed to be constant for each province.

Attenuation Relationships
Magnitude ,distance and site conditions are the principal variables used in predicting future ground motions.A number of predictive relationships derived from regression analysis of strong motion data are available for peak horizontal ground acceleration.The analysis requires a seismic source model of the region and appropriate attenuation relationships.It was decided to use the Campbell ) 1990 ( and Campbell-Bozorgnia (1993) attenuation relationships in this study .These relationships were derived empirically from recorded accelerograms due to earthquakes in different parts of the world.Most of these accelerograms used to drive these relationships are recorded in Western North America, and some other part of the world (including the Tabas and Manjil earthquake in Iran).However , Chandra et. al (1979) show that attenuation of acceleration is more rapid in Iran than Western North America.Therefore, the use of these relationships may be conservative.It should be noted that Zare, Bard, and Ghafory-Ashtiany are developing an attenuation law for Iran based on a new corrected strong motion accelerometer data base.However, it was not ready to be used for this hazard assessment.Therefore, these relationships have been accepted as appropriate models for evaluation of the ground motion parameters in the seismic zonations study.

3.
The least acceleration contours are expected in two regions.One includes a narrow band which has a NW-SE trend and extends from Urumiyeh to Esfahan.The other one is in the central Lut zone in eastern Iran.Corresponding maximum acceleration values for both of these regions are less than 0.20 g, and 0.35g for return periods of 75 and 475 years, respectively.

4.
The peak ground acceleration attenuation relationships considered in this study are given for the mean maximum horizontal acceleration in bedrock.In soil deposits of appreciable depth, the ground accelerations will be somewhat larger than those indicated on these hazard maps.

5.
The maps are intended as a pattern for the assessment of the seismic hazard concerning engineering structures.Moreover, they may also be used for other purposes such as the preparation of seismic risk maps, the estimation of earthquake insurance premiums, and the preliminary site evaluation of critical facilities.

6.
Since only major known sources have been considered in preparation of these maps ,it is recommended that for important type of structures, site specific studies which require deterministic hazard investigations, and identified local seismic sources on large scale maps to be performed.

Figure ( 1 )
indicates the recent seismicity of Iran.It shows the high inhomogeneouty and seismic activity dispersion of the Iranian plateau.No earthquakes larger than Ms= 7.0 have been experienced in the Zagros region during the 20th century, but shocks of magnitude over Ms= 7.0 have occurred in central and eastern Iran.
4, is now agreed by several investigators that the seismicity in Iran is related to the local surface geology and tectonics.Many destructive earthquakes such as Silakhor(Ms=7.4,1909),Salmas (Ms=7.4,1930), Ms=7.3, 1997) confirm this phenomenon in Iran.Therefore, the Iranian plateau is one of the seismically active areas of the world and frequently suffers destructive and catastrophic earthquakes that cause heavy loss of human life and widespread damage.