SEISMIC HAZARD ASSESSMENT OF METROPOLITAN TEHRAN,IRAN

This paper presents a probabilistic seismic hazard assessment of Tehran, the capital of Iran. Two maps have been prepared to indicate the earthquake hazard of Tehran and its vicinity in the form of iso-acceleration contour lines. They display the probabilistic estimate of Peak Ground Acceleration (PGA) over bedrock for the return periods of 475 and 950 years. Tehran is a densely populated metropolitan in which more than 10 million people live. Many destructive earthquakes happened in Iran in the last centuries. It comes from historical references that at least 6 times, Tehran has been destroyed by catastrophic earthquakes. The oldest one happened in the 4th century BC. A collected catalogue, containing both historical and instrumental events and covering the period from the 4th century BC to 1999 is then used. Seismic sources are modelled and recurrence relationship is established. For this purpose the method proposed by Kijko [2000] was employed considering uncertainty in magnitude and incomplete earthquake catalogue. The calculations were performed using the logic tree method and three weighted attenuation relationships; Ramazi [1999], 0.4, Ambraseys and Bommer [1991], 0.35, and Sarma and Srbulov [1996], 0.25. Seismic hazard assessment is then carried out for 12×11 grid points using SEISRISK III. Finally, two seismic hazard maps of the studied area based on Peak Ground Acceleration (PGA) over bedrock for 10% probability of exceedance in two life cycles of 50 and 100 years are presented. The results showed that the PGA ranges from 0.27(g) to 0.46(g) for a return period of 475 years and from 0.33(g) to 0.55(g) for a return period of 950 years. Since population is very dense in Tehran and vulnerability of buildings is high, the risk of future earthquakes will be very significant.     

Seismic Hazard Analysis Using Discrete Faults in Northwestern Pakistan: Part II – Results of Seismic Hazard Analysis

This article is the second of two companion articles that evaluate seismic hazard in northwestern (NW) Pakistan. Using the properties and characteristics of discrete faults in NW Pakistan described in the first article, probabilistic and deterministic seismic hazard analyses for 11 major cities in NW Pakistan were conducted. The results from both probabilistic and deterministic seismic hazard analyses exhibit good agreement. Median deterministic spectra compare favorably with uniform hazard spectra (UHS) for 475- or 975-year return periods, while the 84^th-percentile deterministic spectra compare favorably with the UHS for a 2475-year return period. Peak ground accelerations (PGAs) for 2475-year return periods exceed 1.0 g for the cities of Kaghan and Muzaffarabad, which are surrounded by major faults. The PGAs for a 475-year return period for these cities are approximately 0.6g — 3 to 4 times greater than estimates by previous studies using diffuse areal source zones. The PGAs for some cities located farther from faults (including Astor, Malakand, Mangla, Peshawar, and Talagang) are similar to those predicted using diffuse areal source zones. Seismic hazard maps for PGA and spectral accelerations at periods of 0.2 s and 1.0 s corresponding to three return period (2475, 975, and 475 years) were produced. Based on deaggregation results, a discussion of the conditional mean spectra for engineering applications is presented.     

Stochastic Modeling of Earthquake Scenarios in Greater Tehran

Tehran, the capital of Iran, with millions of inhabitants, has been affected several times by historical and recent earthquakes that confirm the importance of seismic hazard assessment for the area. The main objective of this article is to present a probabilistic procedure to construct time series compatible with the source-path and site reflecting the influence of different magnitude events at different distances that may occur during a specified time period. A Monte Carlo approach is used to generate numerous synthetic catalogs for the evaluation of the probabilistic seismic hazard in greater Tehran over hard rock site for a return period of 475 years. The disaggregation of the seismic hazard is carried out to identify hazard-dominating events and to associate them with one or more specific faults, rather than a given distance. The stochastic finite-fault technique based on region specific seismic parameters is used to generate time series of earthquake scenario.     

Earthquake Risk Scenarios in Urban Areas: A Review with Applications to the Ciutat Vella District in Barcelona,Spain

ABSTRACT

The vulnerability index and capacity spectrum-based methods used to assess the expected physical damage in dwellings are revisited, and their main advantages and shortcomings are analyzed and discussed. Then, a vulnerability index-based method is used to assess the expected damage in the oldest district of Barcelona, Spain. The risk analysis is performed for earthquake scenarios defined by macroseismic intensities of V, V-VI, VI, VI-VII, and VII to design seismic emergency plans. According to the Spanish Seismic Code, the basic peak ground acceleration for a 475-year return period is 0.04 g, which corresponds to an intensity of VI. Thus, the expected physical damage is quantified and its impact on population and other quantities, such as debris and economic cost, are studied. In spite of the low-to-moderate seismic hazard in Barcelona, the results show that the risk is high, due to the high exposure and vulnerability of the built environment.     

Design Earthquake Ground Motion Prediction for Perth Metropolitan Area with Microtremor Measurements for Site Characterization

Perth is the largest city in Western Australia and home to three-quarters of the state's residents. In recent decades, there have been a lot of earthquake activities just east of Perth in an area known as the South-West Seismic Zone. Previous numerical results of site response analyses based on limited available geology information for PMA indicated that Perth Basin might amplify the bedrock motion by more than 10 times at some frequencies and at some sites. Hence, more detailed studies on site characterization and amplification are necessary. The microtremor method using spatial autocorrelation (SPAC) processing is a useful tool for gaining thickness and shear wave velocity (SWV) of sediments and has been adopted in many previous studies. In this study, the response spectrum of rock site corresponding to the 475-year return period for PMA is defined according to the probabilistic seismic hazard analysis (PSHA) based on the latest ground motion attenuation model of Southwest Western Australia. Site characterization in PMA is performed using two microtremor measurements, namely SPAC technique and H/V method. The clonal selection algorithm (CSA) is introduced to perform direct inversion of SPAC curves to determine the soil profiles of representative PMA sites investigated in this study. Using the simulated bedrock motion as input, the responses of the soil sites are estimated using numerical method based on the shear-wave velocity vs. depth profiles determined from the SPAC technique. The response spectrum of the earthquake ground motion on surface of each site is derived from the numerical results of the site response analysis, and compared with the respective design spectrum defined in the Australian Earthquake Loading Code. The comparison shows that the code spectra are conservative in the short period range, but may slightly underestimate the response spectrum at some long period range.     

Seismic Reliability of Code-Conforming Italian Buildings

ABSTRACT

This paper presents and discusses some research results related to the seismic failure risk of standard, residential and industrial, buildings designed for damage, and life-safety according to the Italian seismic code, which is somewhat similar to Eurocode 8. The five considered structural typologies are as follows: masonry, cast-in-place reinforced concrete, precast reinforced concrete, steel, and base-isolated buildings. The archetype structures have been designed according to standard practice at three sites, representative of the seismic hazard across the country. Seismic risk is defined here as the annual rate of earthquakes able to cause structural failure in terms of usability-preventing damage and global collapse. For each structure, the failure rates have been evaluated in the framework of performance-based earthquake engineering, that is, via integration of site’s probabilistic hazard and structural fragility. The former has been computed consistently with the official hazard model for Italy that is also used to define design actions in the code. The latter has been addressed via nonlinear dynamic analysis of three-dimensional numerical structural models. Results indicate that, generally, design procedures are such that seismic structural reliability tends to decrease with increasing seismic hazard of the building site, despite the homogeneous return period of exceedance of the design seismic ground-motion.     

Probabilistic Seismic Hazard Model of West Bengal,India

We deliver a next-generation Probabilistic Seismic Hazard model of West Bengal based on improved seismogenic source characterization considering both the Layered Polygonal sources & Tectonic sources in the hypocentral depth range of 0–25 km & 25–70 km, seismic local-specific site condition, and adaption of appropriate region specific ground motion prediction equations in a logic tree framework. The surface consistent Probabilistic seismic hazard distribution in terms of Peak Ground Acceleration (PGA) & 5% damped Pseudo Spectral Acceleration (PSA) at different time periods for 10% probability of exceedance in 50 years have been generated and the design response spectra computed.     

EXPERIMENTAL STUDY OF TOPOGRAPHIC AMPLIFICATION USING THE ISRAEL SEISMIC NETWORK

Earthquakes and microtremor records are used for estimating the site response of hard rock sites comprising four three-component seismic stations which operate as part of the Israel Seismic Network. The response functions are determined by implementing the horizontal-to-vertical component spectral ratio of earthquake shear-waves (receiver function estimates) and microtremors (Nakamura's estimate) observed simultaneously at the site. The sites of seismic stations ATZ (Mt. Atzmon), MBH (Mt. Berech) and MRNI (Mt. Meron) exhibit amplification attributed to topography effects. At ATZ, within the 1.3–2.0 Hz range, the amplification is in the order of factor 4. At MBH amplification levels of 3.0–3.5 are observed in the frequency range 1.5–4.0 Hz. Station MRNI exhibits a relatively strong amplification effect (up to 4) in the frequency range of about 2.5 to 3.5 Hz. Slight amplification around 5 Hz is observed at ATR (the proposed site for a nuclear power plant). These effects were correlated with the thickness of the weathered layer above unweathered chalk. A comparison between the amplification factor observed during earthquakes and those inferred from microtremors shows that these are, in general, in agreement. However, details of the spectral ratios from different microtremor recordings are not exactly the same. Differences appear mainly in the frequency at which the maximum amplification occurs. These observations demonstrate the usefulness of non-reference technique in estimating the topographical effects of ground shaking. These methods may be used in the process of seismic hazard assessment for ridges and mountain tops, common sites for settlements, communication relay stations, bridges, rope-drive and power transmission towers.     

MICROZONATION OF THE EXPECTED SEISMIC SITE EFFECTS ACROSS PORT VILA,VANUATU

The city of Port Vila, Vanuatu, is located in one of the most active seismic regions on earth. Earthquakes are felt frequently and, due to very rapid plate convergence rates, return period of large earthquakes (M>6) in the New Hebrides Benioff zone can be less than 10 years. Even though Port Vila does not lie on an identified seismic fault zone, strong motions by nearby earthquakes have to be expected due to the city's geographical location close to the plate boundary of New Hebrides convergence zone. An accurate estimation of the seismic ground motion across the city is of prime importance for urban developments and mitigation of earthquake risk. Following many examples of monitored strong earthquakes in the current century, it is evident that the local site effects may have a dominant contribution to the intensity of damage and destruction. In this study we focussed on the first stage of associating site effects and seismic hazard by preparing a microzonation map for Port Vila. The seismic microzonation of the city has been carried out to provide a detailed map of the zones that exhibit site effects in terms of resonance frequencies and approximated amplification of the ground shaking. Having in mind that these data will be used in improving building design to sustain strong ground motions, our analysis is limited to the frequency band of 1–10 Hz, corresponding to the expected resonance of different types of buildings in Port Vila. The Nakamura technique has been used to estimate site amplification effects from single station noise recordings. Interestingly, excluding one site located on an old dump zone, the amplification factors at about the 100 sites surveyed in Port Vila remain below 3 with an average well below 2 in the 1 to 10 Hz frequency band. These results suggest that there is no significant Vs velocity change in consequently layered material and that the uppermost sedimentary layers in the surveyed down town area are relatively thin. These observations are in agreement with the mapping of limestone terraces throughout Port Vila area. However, both the surface geology and results from seismic zonation indicate a thicker (up to several tens of meters) sedimentary cover around the Bauerfield airport and in the Mele terrace zone. Low resonance frequencies (around and below 1 Hz) and amplification factor of the order of 5 were observed over this large area, immediately outside Port Vila. Any building development in this area should take these results into account.     

Earthquake-Induced Damage Detection in a Monumental Masonry Bell-Tower Using Long-Term Dynamic Monitoring Data

ABSTRACT

This work investigates the use of an advanced long-term vibration-based structural health monitoring tool to automatically detect earthquake-induced damages in heritage structures. Damage produced in a monumental bell-tower at increasing values of the Peak Ground Acceleration (PGA) of the seismic input is predicted by incremental nonlinear dynamic analysis, using a Finite Element model calibrated on the basis of experimentally identified natural modes. Then, predicted damage effects are artificially introduced in the monitoring data to check for their detectability. The results demonstrate that a very small damage, associated to a low intensity and low return period earthquake, is clearly detected by the monitoring system.     

GROUND MOTION ANALYSIS IN THE ANCHORAGE BASIN: 1-D APPROACH

Understanding the dynamic behaviour of soil in the Anchorage basin in southcentral Alaska is essential for seismic hazard assessment of this highly seismically active region. The analysis of site responses for 40 sites from weak-motion and strong-motion data with amplitudes less than 0.1 g showed a strong influence of subsurface geological conditions on the characteristics of ground motion. Particularly, the sites in the central part of the city, including the downtown area, showed prominent resonance peaks around 1 Hz with amplification values up to about 4. The numerical analysis, based on one-dimensional multi-layer soil models shows that site response characteristics, and especially aforementioned peaks, are largely related to the thick, soft layer of Quaternary deposits, particularly cohesive fades of Bootlegger Cove Formation. The computed transfer functions for soil profiles of six representative sites are in accord with the site responses in the frequency range from 1 to 5 Hz. There is no significant change in amplification values below 2 Hz corresponding to large-amplitude (up to 0.38 g) ground motions; however, above 2-3 Hz the amplification values are greatly reduced in this case.     

SEISMIC WAVE ATTENUATION FOR VRANCEA EVENTS REVISITED

New aspects of the frequency-dependent attenuation of the seismic waves travelling from Vrancea subcrustal sources toward NW (Transylvanian Basin) and SE (Romanian Plain) are evidenced by the recent experimental data made available by the CALIXTO'99 tomography experiment. The observations validate the previous theoretical computations performed for the assessment, by means of a deterministic approach, of the seismic hazard in Romania. They reveal an essential aspect of the seismic ground motion attenuation that has important implications on the probabilistic assessment of seismic hazard from Vrancea intermediate-depth earthquakes. The attenuation toward NW is shown to be a much stronger frequency-dependent effect than the attenuation toward SE and the seismic hazard computed by the deterministic approach fits satisfactorily well the observed ground motion distribution in the low-frequency band (<1Hz). The apparent contradiction with the historically-based intensity maps arises mainly from a systematic difference in the eigenperiods (type and size) of the buildings in the intra- and extra-Carpathians regions, thus the existing macroseismic data, based on buildings of small dimensions, i.e. with high eigenfrequency (5–10 Hz), can hardly be representative of the real hazard for new and large dimension, tall buildings, with eigenfrequency above 1 Hz.     

EFFECTS OF SOIL LAYERING ON THE CHARACTERISTICS OF BASIN-EDGE INDUCED SURFACE WAVES AND DIFFERENTIAL GROUND MOTION

This paper studies the effects of soil layering in the basin based on the characteristics of basin-edge induced surface waves and associated differential ground motion. Seismic responses of various basin-edge models were simulated using software based on parsimonious finite difference staggered grid approximation of 2.5D eiastodynamic wave equation. Seismic responses of various models with different number of soil/soft rock layers but for a fixed thickness of deposit, fundamental frequency and impedance contrast revealed a decrease of surface wave amplitude with an increase in the number of layers in the basin. Shifting of dominant frequency towards the higher values was obtained with an increase of number of layers. An increase of dispersion of surface waves with an increase of number of soil layers in the basin was observed. A minor increase of Rayleigh wave velocity with an increase of number of soil layers was also obtained, but in the case of Love wave it was almost negligible.

Spectral analysis of the edge-induced surface waves revealed that the anomalous earthquake intensity may arise in a zone of width of 2.5–3.0 km, parallel to basin-edge and at an offset of 0.5–0.7 km from the edge. Maximum horizontal differential ground motion (HDGM) developed by Love wave (≈4.9×10^?2) was more than that of Rayleigh wave (≈9.4×10^?3). Large variation in HDGM caused by the surface waves was obtained with a change in the number of layers in the basin and maximum HDGM was observed when there were only two layers in the basin. It was inferred that the effect of soil layering in the basin was more on the Rayleigh wave as compared to the Love wave. Development of large HDGM near the basin-edge and its dependency on the number of soil layers reveals that basin-edge induced surface waves need special attention during seismic microzonation or seismic hazard prediction.     

Modeling and Seismic Response Analysis of RC Precast Italian Code-Conforming Buildings

ABSTRACT

In this study, industrial single-story RC precast buildings are investigated. Twenty-four case studies have been considered, in which the column height, the beam spans and the seismic hazard level are varied. The seismic design of the selected case studies is performed according to the Italian building code and additional technical documentation. Three-dimensional nonlinear models are defined to perform static and dynamic analyses for the seismic assessment of the selected case studies. Demand/capacity ratios in terms of the selected engineering demand parameters are computed for ten increasing values of the seismic input return period.     

SEISMIC HAZARD IN GREECE BASED ON DIFFERENT STRONG GROUND MOTION PARAMETERS

The available Greek strong ground motion records to date are used in order to study the duration of strong-motion in Greece, covering magnitudes between 4.5 and 6.9 and distances from 1 km to 128 km. An attenuation relation of strong-motion duration is calculated and compared to earlier existing similar relations proposed for Greece and Japan. Furthermore, the seismic hazard for the area of Greece is assessed, using the strong-motion parameters of duration and peak ground acceleration. The results are presented in the form of a map according to which Greece is classified in four different categories of equal seismic hazard.     

SEISMIC HAZARD ASSESSMENT IN FLORENCE CITY ITALY

A seismic hazard analysis of Florence city was performed in the frame of a project concerning the dynamic behaviour of cable-stayed bridges. Both a probabilistic approach and a methodology based on the use of a local macroseismic catalogue were applied. A local catalogue was expressly compiled for this purpose, to collect the macroseismic intensities actually observed at the site as a result of past earthquakes. This sort of catalogue is an independent tool to verify the assumptions of the probabilistic approach (seismic zoning, earthquake recurrence relation, attenuation model), though it can supply results in terms of macroseismic intensity only and reflects the effective seismic history at the site, without taking into account any variability. The Cornell' methodology was used to assess probabilistic hazard in terms of macroseismic intensity, peak ground acceleration, peak ground velocity, and pseudovelocity uniform response spectra. The local catalogue points out level VII of the Mercalli-Cancani-Sieberg scale (MCS) as the maximum intensity historically observed in Florence. The probabilistic approach leads to the consideration of intensity VIII MCS as the maximum credible for the city. The probabilistic analysis in terms of ground motion was performed using attenuation relations estimated for alluvium sites, since the geology of Florence area is represented by fluvial and lacustrine deposits of various thickness. Peak ground acceleration values with 90% non exceedence probability in 50 and 500 years are respectively 145 and 219 cm/s's for a shallow alluvium site, and 95 and 157 cm/s's for a deep alluvium site; the corresponding peak ground velocity values for sites located on alluvium are 6.41 and 11.76 cm/s. Uniform response spectra are provided for shallow and deep alluvium sites, according to frequency-dependent attenuation relations estimated from strong Italian earthquakes.     

ATTENUATION RELATIONSHIP FOR LOW MAGNITUDE EARTHQUAKES USING STANDARD SEISMOMETRIC RECORDS

Northwestern Italian weak-motion data were used to study attenuation characteristics of horizontal peak ground acceleration (PGA) and horizontal peak ground velocity (PGV) from earthquakes of local magnitudes (M _l ) up to 5.1. Data have been provided by the RSNI (Regional seismic network of Northwestern Italy) and RSLG (Regional seismic network of Lunigiana-Garfagnana) waveform database. The database consists of more than 14000 horizontal components recorded in the period 1999-2002 by both broadband and enlarged band seismometers. The accuracy of the procedure used to extract PGA values from the velocity recordings was verified comparing observed and derived PGA values at station STV2, which was equipped with both a temporary K2 Kinemctrtcs accelerometer and Guralp CMG40 broadband sensor. The attenuation of both peak ground acceleration and peak ground velocity was found to be logarithmically distributed with a strong attenuation for low distances (less than 50 km) and low M _l values (<3.0). The resulting equations are:

Log(PGA)=?3.19+0.87M?0.042M ²?1.92 Log(R)+0.249S,

Log(PGA)=?4.23+0.76M?0.018M²?1.56 Log(R)+0.230S,

where PGA is expressed in g, PGV is expressed in m/s, M is local magnitude, R is the hypocentral distance in kilometers and S is a dummy variable assuming values of 0 and 1 for rock and soil respectively. For increasing distance and magnitude, both PGA and PGV values show a linear distribution. The validity range of the obtained attenuation relationships is 0–200 km for distances and M _l up to 4.5. Sensitivity studies performed by analysis of residuals, showed that predicted PGA and PGV values are stable with respect to reasonable variations of the model and distances providing the data. Comparisons with attenuation relationships proposed for Italian region, derived from strong motion records, are also presented.     

OBTAINING PROBABILISTIC ESTIMATES OF DISPLACEMENT ON A LANDSLIDE DURING FUTURE EARTHQUAKES

The assessment of earthquake triggered landslide hazard may be undertaken using both deterministic and probabilistic techniques. Probabilistic methods have been developed because much of the data can be considered as random variables where parameters such as the angle of internal friction and moisture content do not have a single fixed value but may assume any number of values across a range. This random variability can be modelled by a probability density function (PDF) which describes the relative likeli-hood that a random variable will assume a particular value. Instead of using just the average or expected value of an input parameter, the complete range of possible values can be used to estimate a range of possible outcomes. Thus the probability of a slope being unstable can be obtained rather than a single indicator of stability. Such proba-bilistic analyses allow for the incorporation of the likely variability of each parameter and therefore allow a more intimate assessment of slope stability to be derived. Utilising empirical relationships for calculating earthquake ground motions and associated slope displacement, an investigation was undertaken to identify the contribution that modern simulation techniques could make to the assessment of earthquake-triggered landslides. To achieve this, geotechnical and earthquake data obtained from a deep-seated landslide triggered during the M _w 7.0 Loma Prieta earthquake was used. By incorporating the variability of the geotechnical parameters and the uncertainty in earthquake location the model derived the probabilities associated with increasing amounts of slope displacement during future probable earthquakes. Analysis was undertaken for four of the principal fault segments in the San Francisco Bay area. These estimates were then combined with the occurrence probabilities of the earthquakes to provide temporal estimates of dis-placement for a 30 year period. Results indicated that a M _w 7.0 earthquake located on the Peninsula Segment of the San Andreas fault was most hazardous with a 11% chance of minor slope displacement (≥0.10 m) and a 6% chance of moderate slope displacement (≥0.30 m) within the next 30 years.     

Full-Scale Shaking Table Tests on a Substandard RC Building Repaired and Strengthened with Post-Tensioned Metal Straps

The effectiveness of a novel Post-Tensioned Metal Strapping (PTMS) technique at enhancing the seismic behavior of a substandard RC building was investigated through full-scale, shake-table tests during the EU-funded project BANDIT. The building had inadequate reinforcement detailing in columns and joints to replicate old construction practices. After the bare building was initially damaged significantly, it was repaired and strengthened with PTMS to perform additional seismic tests. The PTMS technique improved considerably the seismic performance of the tested building. While the bare building experienced critical damage at an earthquake of PGA = 0.15 g, the PTMS-strengthened building sustained a PGA = 0.35 g earthquake without compromising stability.     

Empirical Correlations of Frequency-Content Parameters of Ground Motions with Other Intensity Measures

Empirical correlations between the frequency-content parameters of earthquake ground motions and amplitude-, cumulative-, and duration-based intensity measures (IMs) are examined in this study. Three commonly used scalar frequency-content parameters are considered, namely the mean period (T_m), the average spectral period (T_avg), and the smoothed spectral predominant period (T₀). It is found that the frequency-content parameters have weakly negative correlations with high-frequency IMs such as peak ground acceleration (PGA) and spectral accelerations (SAs) at periods smaller than 0.3 s, low-to-moderate positive correlations with peak ground velocity (PGV) and SA within a period range of 0.5 s–10 s, negligible correlations with cumulative-based IMs, and weakly positive correlations (in the vicinity of 0.1–0.3) with significant durations. Simple piecewise parametric equations are proposed to fit the empirical correlations of T_m, T_avg, and T₀ with SA over the entire period range. The presented correlation results and parametric models enable the frequency-content parameters to be easily used in various applications such as ground-motion selection and vector-based probabilistic seismic hazard analysis.